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

#include <unitarity.h>

Inheritance diagram for unitarityRp13:

Detailed Description

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

Definition at line 1364 of file unitarity.h.

Public Member Functions
double	computeThValue ()

	unitarityRp13 (const StandardModel &SM_i)
	unitarityRp13 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

◆ unitarityRp13()

unitarityRp13::unitarityRp13 ( const StandardModel & SM_i )

unitarityRp13 constructor.

Definition at line 2629 of file unitarity.cpp.

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

2631 {}

Member Function Documentation

◆ computeThValue()

double unitarityRp13::computeThValue ( )

virtual

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

Implements ThObservable.

Definition at line 2633 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 WFRc3 = myTHDM.getMyTHDMCache()->WFRcomb3;
     double WFRc4 = myTHDM.getMyTHDMCache()->WFRcomb4;
     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 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 uniNLO15a = -la1Q/(16.0*M_PI);
     double uniNLO15b = 3.0*betalambda1/(256.0*M_PI*M_PI*M_PI);
     gslpp::complex uniNLO15c = (gslpp::complex::i()*M_PI-1.0)*(la1Q*la1Q+la5Q*la5Q)/(256.0*M_PI*M_PI*M_PI);
     gslpp::complex uniNLO15d = -la1Q/(512.0*M_PI*M_PI*M_PI) * (WFRc1-2.0*WFRc2+WFRc3+2.0*WFRc4);
  
     double uniNLO16a = -la2Q/(16.0*M_PI);
     double uniNLO16b = 3.0*betalambda2/(256.0*M_PI*M_PI*M_PI);
     gslpp::complex uniNLO16c = (gslpp::complex::i()*M_PI-1.0)*(la2Q*la2Q+la5Q*la5Q)/(256.0*M_PI*M_PI*M_PI);
     gslpp::complex uniNLO16d = -la2Q/(512.0*M_PI*M_PI*M_PI) * (-WFRc1+2.0*WFRc2-WFRc3+2.0*WFRc4);
  
     double uniNLO17a = -la5Q/(16.0*M_PI);
     double uniNLO17b = 3.0*betalambda5/(256.0*M_PI*M_PI*M_PI);
     gslpp::complex uniNLO17c = (gslpp::complex::i()*M_PI-1.0)*(la1Q+la2Q)*la5Q/(256.0*M_PI*M_PI*M_PI);
     gslpp::complex uniNLO17d = -2.0*la5Q/(512.0*M_PI*M_PI*M_PI) * WFRc4;
  
     gslpp::complex uniA=uniNLO15a+uniNLO15b+uniNLO15c+uniNLO15d;
     gslpp::complex uniB=uniNLO16a+uniNLO16b+uniNLO16c+uniNLO16d;
     gslpp::complex uniC=uniNLO17a+uniNLO17b+uniNLO17c+uniNLO17d;
  
     double Rp=0.01;
     double a0LO = uniNLO15a+uniNLO16a+sqrt((uniNLO15a-uniNLO16a)*(uniNLO15a-uniNLO16a)+4.0*uniNLO17a*uniNLO17a);
     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& unitarityRp13::myTHDM

private

Definition at line 1377 of file unitarity.h.

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