unitarityRp1 Class Reference

NLO over LO ratio for the first Z2-even 00 eigenvalue (corresponding to unitarityNLOev1). More...

#include <unitarity.h>

Inheritance diagram for unitarityRp1:

Detailed Description

NLO over LO ratio for the first Z2-even 00 eigenvalue (corresponding to unitarityNLOev1).

Definition at line 1128 of file unitarity.h.

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

Private Attributes
const THDM &	myTHDM

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

Constructor & Destructor Documentation

unitarityRp1()

unitarityRp1::unitarityRp1

(

const StandardModel &

SM_i

)

unitarityRp1 constructor.

Definition at line 1967 of file unitarity.cpp.

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

Member Function Documentation

computeThValue()

double unitarityRp1::computeThValue ( )

virtual

Returns

\(|a^{\text{even,NLO}}_{00+}/a^{\text{even,LO}}_{00+}|\)

Implements ThObservable.

Definition at line 1971 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 WFRc1 = myTHDM.getMyTHDMCache()->WFRcomb1;
    double WFRc2 = myTHDM.getMyTHDMCache()->WFRcomb2;
    double Rpeps = myTHDM.getMyTHDMCache()->Rpeps;
 
    double betalambda1 = 6.0*la1Q*la1Q + 2.0*la3Q*la3Q + 2.0*la3Q*la4Q + la4Q*la4Q + la5Q*la5Q
                          + 6.0*la1Q*Yb1Q*Yb1Q + 2.0*la1Q*Ytau1Q*Ytau1Q
                          - 6.0*Yb1Q*Yb1Q*Yb1Q*Yb1Q - 2.0*Ytau1Q*Ytau1Q*Ytau1Q*Ytau1Q;
    double betalambda2 = 6.0*la2Q*la2Q + 2.0*la3Q*la3Q + 2.0*la3Q*la4Q + la4Q*la4Q + la5Q*la5Q 
                          + 6.0*la2Q*Yb2Q*Yb2Q + 2.0*la2Q*Ytau2Q*Ytau2Q + 6.0*la2Q*YtQ*YtQ
                          - 6.0*Yb2Q*Yb2Q*Yb2Q*Yb2Q - 2.0*Ytau2Q*Ytau2Q*Ytau2Q*Ytau2Q - 6.0*YtQ*YtQ*YtQ*YtQ;
    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 betalambda4 = la1Q*la4Q + la2Q*la4Q + 4.0*la3Q*la4Q + 2.0*la4Q*la4Q + 4.0*la5Q*la5Q
                          + 3.0*la4Q*Yb1Q*Yb1Q + 3.0*la4Q*Yb2Q*Yb2Q + la4Q*Ytau1Q*Ytau1Q + la4Q*Ytau2Q*Ytau2Q + 3.0*la4Q*YtQ*YtQ
                          - 2.0*Ytau1Q*Ytau1Q*Ytau2Q*Ytau2Q - 6.0*Yb1Q*Yb1Q*(Yb2Q*Yb2Q - YtQ*YtQ);
 
    double uniNLO1a = -3.0*la1Q/(16.0*M_PI);
    double uniNLO1b = 9.0*betalambda1/(256.0*M_PI*M_PI*M_PI);
    gslpp::complex uniNLO1c = (gslpp::complex::i()*M_PI-1.0)*(9.0*la1Q*la1Q+(2.0*la3Q+la4Q)*(2.0*la3Q+la4Q))/(256.0*M_PI*M_PI*M_PI);
    gslpp::complex uniNLO1d = -3.0*la1Q/(512.0*M_PI*M_PI*M_PI) * WFRc1;
 
    double uniNLO2a = -3.0*la2Q/(16.0*M_PI);
    double uniNLO2b = 9.0*betalambda2/(256.0*M_PI*M_PI*M_PI);
    gslpp::complex uniNLO2c = (gslpp::complex::i()*M_PI-1.0)*(9.0*la2Q*la2Q+(2.0*la3Q+la4Q)*(2.0*la3Q+la4Q))/(256.0*M_PI*M_PI*M_PI);
    gslpp::complex uniNLO2d = -3.0*la2Q/(512.0*M_PI*M_PI*M_PI) * (-WFRc1+2.0*WFRc2);
 
    double uniNLO3a = -(2.0*la3Q+la4Q)/(16.0*M_PI);
    double uniNLO3b = 3.0*(2.0*betalambda3+betalambda4)/(256.0*M_PI*M_PI*M_PI);
    gslpp::complex uniNLO3c = 3.0*(gslpp::complex::i()*M_PI-1.0)*(la1Q+la2Q)*(2.0*la3Q+la4Q)/(256.0*M_PI*M_PI*M_PI);
    gslpp::complex uniNLO3d = -(2.0*la3Q+la4Q)/(512.0*M_PI*M_PI*M_PI) * WFRc2;
 
    gslpp::complex uniA=uniNLO1a+uniNLO1b+uniNLO1c+uniNLO1d;
    gslpp::complex uniB=uniNLO2a+uniNLO2b+uniNLO2c+uniNLO2d;
    gslpp::complex uniC=uniNLO3a+uniNLO3b+uniNLO3c+uniNLO3d;
 
    double Rp=0.01;
    double a0LO = uniNLO1a+uniNLO2a+sqrt((uniNLO1a-uniNLO2a)*(uniNLO1a-uniNLO2a)+4.0*uniNLO3a*uniNLO3a);
    if(fabs(a0LO)>Rpeps)
    {
        Rp = ((uniA+uniB+sqrt((uniA-uniB)*(uniA-uniB)+4.0*uniC*uniC))-a0LO).abs()/fabs(a0LO);
    }
 
    return Rp;
}

Member Data Documentation

myTHDM

const THDM& unitarityRp1::myTHDM

private

Definition at line 1141 of file unitarity.h.

---

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

• unitarity.h
• unitarity.cpp