unitarityR12 Class Reference

#include <unitarity.h>

Inheritance diagram for unitarityR12:

Detailed Description

Definition at line 1734 of file unitarity.h.

Public Member Functions
double	computeThValue ()
	unitarityR12 (const StandardModel &SM_i)
	unitarityR12 constructor. More...
Public Member Functions inherited from ThObservable
double	getBinMax ()
	A get method to get the maximum value of the bin. More...
double	getBinMin ()
	A get method to get the minimum value of the bin. More...
const StandardModel &	getModel ()
	A get method to get the model. More...
const std::vector< std::string >	getParametersForObservable ()
	A get method to get the parameters for the specific observable. More...
void	setBinMax (double max)
	A set method to set the maximum value of the bin. More...
void	setBinMin (double min)
	A set method to set the minimum value of the bin. More...
void	setParametersForObservable (std::vector< std::string > parametersForObservable_i)
	A set method to get the parameters for the specific observable. More...
	ThObservable (const StandardModel &SM_i)
	Constructor. More...
	ThObservable (const ThObservable &orig)
	The copy constructor. More...
virtual	~ThObservable ()
	The default destructor. More...

Private Attributes
const THDM &	myTHDM

Additional Inherited Members
Protected Attributes inherited from ThObservable
double	max
	the bin maximum. More...
double	min
	The bin minimum. More...
std::vector< std::string >	parametersForObservable
	a vector of parameter namesfor the specific observable More...
const StandardModel &	SM
	A reference to an object of StandardMode class. More...

Constructor & Destructor Documentation

unitarityR12()

unitarityR12::unitarityR12

(

const StandardModel &

SM_i

)

unitarityR12 constructor.

Definition at line 3696 of file unitarity.cpp.

: ThObservable(SM_i),myTHDM(static_cast<const THDM&> (SM_i))
{}

Member Function Documentation

computeThValue()

double unitarityR12::computeThValue ( )

virtual

Implements ThObservable.

Definition at line 3700 of file unitarity.cpp.

{
    double YtQ = myTHDM.getMyTHDMCache()->Ytop_at_Q;
    double Yb1Q = myTHDM.getMyTHDMCache()->Ybottom1_at_Q;
    double Yb2Q = myTHDM.getMyTHDMCache()->Ybottom2_at_Q;
    double Ytau1Q = myTHDM.getMyTHDMCache()->Ytau1_at_Q;
    double Ytau2Q = myTHDM.getMyTHDMCache()->Ytau2_at_Q;
    double la1Q = myTHDM.getMyTHDMCache()->lambda1_at_Q;
    double la2Q = myTHDM.getMyTHDMCache()->lambda2_at_Q;
    double la3Q = myTHDM.getMyTHDMCache()->lambda3_at_Q;
    double la4Q = myTHDM.getMyTHDMCache()->lambda4_at_Q;
    double la5Q = myTHDM.getMyTHDMCache()->lambda5_at_Q;
    double WFRc2 = myTHDM.getMyTHDMCache()->WFRcomb2;
    double WFRc3 = myTHDM.getMyTHDMCache()->WFRcomb3;
    double Rpeps = myTHDM.getMyTHDMCache()->Rpeps;
 
    double betalambda3 = 3.0*la1Q*la3Q + 3.0*la2Q*la3Q + 2.0*la3Q*la3Q + la1Q*la4Q + la2Q*la4Q + la4Q*la4Q + la5Q*la5Q
                          + 3.0*la3Q*Yb1Q*Yb1Q + 3.0*la3Q*Yb2Q*Yb2Q + la3Q*Ytau1Q*Ytau1Q + la3Q*Ytau2Q*Ytau2Q + 3.0*la3Q*YtQ*YtQ
                          - 6.0*Yb1Q*Yb1Q*(Yb2Q*Yb2Q + YtQ*YtQ) - 2.0*Ytau1Q*Ytau1Q*Ytau2Q*Ytau2Q;
    double betalambda5 = la5Q*(la1Q + la2Q + 4.0*la3Q + 6.0*la4Q)
                          + 3.0*la5Q*Yb1Q*Yb1Q + la5Q*Ytau1Q*Ytau1Q + la5Q*(3.0*Yb2Q*Yb2Q + Ytau2Q*Ytau2Q + 3.0*YtQ*YtQ)
                          - 6.0*Yb1Q*Yb1Q*Yb2Q*Yb2Q - 2.0*Ytau1Q*Ytau1Q*Ytau2Q*Ytau2Q;
 
    double uniNLO12a = -la3Q/(16.0*M_PI);
    double uniNLO12b = 3.0*betalambda3/(256.0*M_PI*M_PI*M_PI);
    gslpp::complex uniNLO12c = (gslpp::complex::i()*M_PI-1.0)*(la3Q*la3Q+la5Q*la5Q)/(256.0*M_PI*M_PI*M_PI);
    gslpp::complex uniNLO12d = -la3Q/(512.0*M_PI*M_PI*M_PI) * WFRc3;
 
    double uniNLO13a = -la5Q/(16.0*M_PI);
    double uniNLO13b = 3.0*betalambda5/(256.0*M_PI*M_PI*M_PI);
    gslpp::complex uniNLO13c = (gslpp::complex::i()*M_PI-1.0)*la3Q*la5Q/(128.0*M_PI*M_PI*M_PI);
    gslpp::complex uniNLO13d = -la4Q/(512.0*M_PI*M_PI*M_PI) * WFRc2;
 
    gslpp::complex uniA=uniNLO12a+uniNLO12b+uniNLO12c+uniNLO12d;
    gslpp::complex uniC=uniNLO13a+uniNLO13b+uniNLO13c+uniNLO13d;
 
    double R=0.01;
    double a0LO = uniNLO12a-uniNLO13a;
    gslpp::complex a0NLO = uniA-uniC;
    if(a0NLO.abs()>Rpeps)
    {
        R = (a0NLO-a0LO).abs()/(a0NLO.abs());
    }
 
    return R;
}

Member Data Documentation

myTHDM

const THDM& unitarityR12::myTHDM

private

Definition at line 1744 of file unitarity.h.

---

The documentation for this class was generated from the following files:

• unitarity.h
• unitarity.cpp