Rb0GTHDM Class Reference

#include <GeneralTHDMEWPO.h>

Inheritance diagram for Rb0GTHDM:

Detailed Description

Definition at line 30 of file GeneralTHDMEWPO.h.

Public Member Functions
double	computeThValue ()
	Rb0GTHDM (const StandardModel &SM_i)
	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 GeneralTHDM &	myGTHDM

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

Rb0GTHDM()

Rb0GTHDM::Rb0GTHDM

(

const StandardModel &

SM_i

)

Constructor.

Definition at line 26 of file GeneralTHDMEWPO.cpp.

: ThObservable(SM_i),myGTHDM(static_cast<const GeneralTHDM&> (SM_i)) 
{}

Member Function Documentation

computeThValue()

double Rb0GTHDM::computeThValue ( )

virtual

Returns

Rb0GTHDM

From: arXiv: 1002.1071v2 and "QCD corrections to e+e- cross section and Z decay rates: Concepts and results"

Masses of the top and bottom quarks at the scale of the W mass (mu)

Loop functions

Implements ThObservable.

Definition at line 33 of file GeneralTHDMEWPO.cpp.

{
    
    if (!myGTHDM.getATHDMflag())
    {
        throw std::runtime_error("Rb only aviable in the A2HDM.");
    }
    
    double GF = myGTHDM.getGF();
    double as = myGTHDM.getAlsMz();
    double a = myGTHDM.getAle();
 
    double MW=myGTHDM.Mw_tree(); //The Mw mass is used do define the electroweak scale, we can  
                                 //use the tree level value and forget about the NP contributions
 
    double MZ=myGTHDM.getMz();
    
    /*square mass of the charged Higgs*/
    
    double mHp2=myGTHDM.getmHp2();
    
    /*up and down couplings */
    gslpp::complex u = myGTHDM.getNu_11();
    gslpp::complex d = myGTHDM.getNd_11();
    
    /*The following inputs are needed but are not defined before*/
    
    double Cb = 1.0086;
 
    double mu = MW;
    double mtMZ = myGTHDM.Mrun(mu, myGTHDM.getQuarks(QCD::TOP).getMass_scale(),
                              myGTHDM.getQuarks(QCD::TOP).getMass(), QCD::TOP, FULLNNLO);   
    double mbMZ = myGTHDM.Mrun(mu, myGTHDM.getQuarks(QCD::BOTTOM).getMass_scale(),
                              myGTHDM.getQuarks(QCD::BOTTOM).getMass(), QCD::BOTTOM, FULLNNLO);
   /* std::cout << "mtMZ = " << mtMZ << std::endl;
    std::cout << "mbMZ = " << mbMZ << std::endl;
    std::cout << "GF = " << GF << std::endl;
    std::cout << "MZ = " << MZ << std::endl;
    std::cout << "MW = " << MW << std::endl;
    std::cout << "as = " << as << std::endl;
    std::cout << "a = " << a << std::endl;*/
 
 
 
    double xH = (mtMZ*mtMZ)/(mHp2);
    double pi = M_PI;
    
    
    double Qb = myGTHDM.getQuarks(QCD::BOTTOM).getCharge();
//    double Qt = myGTHDM.getQuarks(QCD::TOP).getCharge();
    
 
       
    /*double gAb = SM.getQuarks(QCD::BOTTOM).getIsospin();
    double gVb = SM.getQuarks(QCD::BOTTOM).getIsospin()-2.0*SM.getQuarks(QCD::BOTTOM).getCharge()*SM.sW2();*/
    double gbLSM = -0.42112;
    double gbRSM = 0.07744;
 
    double Sb = 1.3214;
    
  
        
        double f1 =(xH*xH - xH - xH*log(xH))/((1.0 - xH)*(1.0 - xH));
        
        double f2 =  - (6.0*xH*(xH-2.0)*gsl_sf_dilog(1.0-1./xH))/((xH-1.0)*(xH-1.0))+ xH*(-27.0+11.0*xH)/((xH-1.0)*(xH-1.0))+ xH*(25.0-9.0*xH)*log(xH)/((xH-1.0)*(xH-1.0)*(xH-1.0))+ (6.0*xH*(3.0-xH)/((xH-1.0)*(xH-1.0)) - (12.0*xH*log(xH))/((xH-1.0)*(xH-1.0)*(xH-1.0)))*log((mtMZ*mtMZ)/(MZ*MZ)) - 3.0*f1;
        
        /*Couplings*/
        
        double gbL= gbLSM + (sqrt(2.0)*GF*mtMZ*mtMZ*(u).abs2())/(16.0*pi*pi)*(f1 + (as/(3.0*pi))*f2);
        
        double gbR = gbRSM - (sqrt(2.0)*GF*mbMZ*mbMZ*(d).abs2())/(16.*pi*pi)*(f1+ (as/(3.0*pi))*f2);
        double sb = ((gbL- gbR)*(gbL- gbR)+ (gbL+ gbR)*(gbL+ gbR))*(1.0 + 3.0*a*Qb*Qb/(4.0*pi));
        
        
   /* std::cout << "f1 = " << f1 << std::endl;
    std::cout << "f2 = " << f2 << std::endl;
    std::cout << "gbL = " << gbL << std::endl;
    std::cout << "gbR = " << gbR << std::endl;
    std::cout << "sb = " << sb << std::endl;*/
 
        
        /*Rb*/
        
        double Rb0GGTHDM = 1.0/(1.0 + (Sb*Cb)/sb);
        // std::cout << "Rb0GGTHDM = " << Rb0GGTHDM << std::endl;
 
        return Rb0GGTHDM;
 
       
       /* double DeltaRb0=1.0;
        double Rb0SM=myGTHDM.R0_f(SM.getQuarks(SM.BOTTOM));
        return Rb0SM+DeltaRb0;*/
    
}

Member Data Documentation

myGTHDM

const GeneralTHDM& Rb0GTHDM::myGTHDM

private

Definition at line 43 of file GeneralTHDMEWPO.h.

---

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

• GeneralTHDMEWPO.h
• GeneralTHDMEWPO.cpp