THDMWDeltaS Class Reference

#include <THDMWSTU.h>

Inheritance diagram for THDMWDeltaS:

Detailed Description

Definition at line 60 of file THDMWSTU.h.

Public Member Functions
double	computeThValue ()
	THDMW contribution to \(S\). More...
	THDMWDeltaS (const StandardModel &SM_i)
	Constructor for DeltaS. More...
Public Member Functions inherited from THDMWSTU
double	F (const double m02, const double m12) const
	Function \(F(m02,m12)\) used for THDMW. Remember that this function is defined for THDMW and SUSY while for THDM we have a multiplicative factor 0.5. More...
	THDMWSTU (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...

Additional Inherited Members
Public Attributes inherited from THDMWSTU
const THDMW &	myTHDMW
Protected Attributes inherited from THDMWSTU
THDMWcache *	mycache
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

THDMWDeltaS()

THDMWDeltaS::THDMWDeltaS

(

const StandardModel &

SM_i

)

Constructor for DeltaS.

Definition at line 49 of file THDMWSTU.cpp.

: THDMWSTU(SM_i)
{}

Member Function Documentation

computeThValue()

double THDMWDeltaS::computeThValue ( )

virtual

THDMW contribution to \(S\).

Returns

\(\Delta S\)

Reimplemented from THDMWSTU.

Definition at line 53 of file THDMWSTU.cpp.

{
    double mSp2=myTHDMW.getMyTHDMWCache()->mSpsq;
    double mSR2=myTHDMW.getMyTHDMWCache()->mSRsq;
    double mSI2=myTHDMW.getMyTHDMWCache()->mSIsq;
    double s_W2 = myTHDMW.sW2();
    //std::cout<< "SinW_2=" << s_W2 << std::endl;//Why this value?
    double c_W2 = 1-myTHDMW.sW2();
    double MZ2=pow(myTHDMW.getMz(),2);
    //double MW2=pow(myTHDMW->Mw(),2);
    //We define the self energies removing the coupling because it will cancel
    
    
    
    //gslpp::complex delPiWW_MZ;
    gslpp::complex delPiZZ_MZ;
    //gslpp::complex delPiWW_0;
    gslpp::complex delPiZZ_0;
    gslpp::complex delPiAAprime_0;
    gslpp::complex delPiZAprime_0;
    //gslpp::complex delPiAA_MZ;
    //gslpp::complex delPiZA_MZ;
    
    //g1 on arxiv:0907.2696 is e/s_W
    //delPiWW_MZ=(1./(2.*pow(M_PI,2)*s_W2))*(mycache->B00_MZ2_MZ2_mSi2_mSp2(MZ2,mSI2,mSp2)+mycache->B00_MZ2_MZ2_mSr2_mSp2(MZ2,mSR2,mSp2)-0.5*mycache->A0_MZ2_mSp2(MZ2,mSp2)-0.25*mycache->A0_MZ2_mSp2(MZ2,mSR2)-0.25*mycache->A0_MZ2_mSp2(MZ2,mSI2));
    
        
    //delPiWW_0=(1./(8.*pow(M_PI,2)*s_W2))*(0.5*F(mSp2,mSR2)+0.5*F(mSp2,mSI2));
    
    
    //delPiAA_MZ=(2./pow(M_PI,2))*(mycache->B00_MZ2_MZ2_mSp2_mSp2(MZ2,mSp2)-0.5*mycache->A0_MZ2_mSp2(MZ2,mSp2));
    //delPiZA_MZ=((1.-2.*s_W2)/(pow(M_PI,2)*sqrt(c_W2)*sqrt(s_W2)))*(mycache->B00_MZ2_MZ2_mSp2_mSp2(MZ2,mSp2)-0.5*mycache->A0_MZ2_mSp2(MZ2,mSp2));
 
    
    //double nu2=myTHDMW->getTHDMW_nu2();
    //double nu3=myTHDMW->getTHDMW_nu3();
    //double v=246.2;
    
    //std::cout<< "prove10=" << delPiAAprime_0 << std::endl;
    //std::cout<< "prove11=" << delPiAA_MZ/MZ2 << std::endl;
    //std::cout<< "prove20=" << delPiZAprime_0 << std::endl;
    //std::cout<< "prove21=" << delPiZA_MZ/MZ2 << std::endl;
    //std::cout<< "DeltaS approx.=" << pow(v,2)*nu2/(6.*M_PI*mSp2) << std::endl;
    //std::cout<< "DeltaS old=" << (16.*M_PI*s_W2*c_W2)*(((delPiZZ_MZ.real()-delPiZZ_0.real())/MZ2)-(c_W2-s_W2)*delPiZAprime_0.real()/(sqrt(s_W2)*sqrt(c_W2))-delPiAAprime_0.real()) << std::endl;
 
    //std::cout<< "DeltaS fraction old=" << (16.*M_PI*s_W2*c_W2)*(((delPiZZ_MZ.real()-delPiZZ_0.real())/MZ2)-(c_W2-s_W2)*delPiZAprime_0.real()/(sqrt(s_W2)*sqrt(c_W2))-delPiAAprime_0.real())/(pow(v,2)*nu2/(6.*M_PI*mSp2)) << std::endl;
    //std::cout<< "DeltaS fraction new=" << ((16.*M_PI*s_W2*c_W2)*(((delPiZZ_MZ.real()-delPiZZ_0.real())/MZ2)-(c_W2-s_W2)*delPiZA_MZ.real()/(sqrt(s_W2)*sqrt(c_W2)*MZ2)-delPiAA_MZ.real()/MZ2))/(pow(v,2)*nu2/(6.*M_PI*mSp2)) << std::endl;
 
    
    //std::cout<< "B00_MZ2_MZ2_mSr2_mSi2(MZ2,mSR2,mSI2)=" << mycache->B00_MZ2_MZ2_mSr2_mSi2(MZ2,mSR2,mSI2) << std::endl;
    //std::cout<< "A0_MZ2_mSp2=" << mycache->A0_MZ2_mSp2(MZ2,mSp2) << std::endl;
    //std::cout<< "mSp2=" << mSp2 << std::endl;
    //std::cout<< "mSr2=" << mSR2 << std::endl;
    //std::cout<< "MZ2=" << MZ2 << std::endl;
    //std::cout<< "THDMWpositiveMassSquares=" << myTHDMW->getMyTHDMWCache()->setOtherParameters() << std::endl;
    //std::cout<< "THDMWpositiveMassSquares=" << THDMWpositiveMassSquares << std::endl;
    if(fabs(myTHDMW.getMyTHDMWCache()->setOtherParameters()-1.)<1.e-10)
    {
        delPiZZ_MZ=(1./(2.*pow(M_PI,2)*c_W2*s_W2))*(pow(1.-2.*s_W2,2)*(mycache->B00_MZ2_MZ2_mSp2_mSp2(MZ2,mSp2)-0.5*mycache->A0_MZ2_mSp2(MZ2,mSp2))+mycache->B00_MZ2_MZ2_mSr2_mSi2(MZ2,mSR2,mSI2)-0.25*mycache->A0_MZ2_mSp2(MZ2,mSR2)-0.25*mycache->A0_MZ2_mSp2(MZ2,mSI2));    
        delPiZZ_0=(1./(8.*pow(M_PI,2)*s_W2*c_W2))*(0.5*F(mSI2,mSR2));
        delPiAAprime_0=-(1/(6.*pow(M_PI,2)))*mycache->B0_MZ2_0_mSp2_mSp2(MZ2,mSp2);
        delPiZAprime_0=-((1-2*s_W2)/(12.*pow(M_PI,2)*sqrt(c_W2)*sqrt(s_W2)))*mycache->B0_MZ2_0_mSp2_mSp2(MZ2,mSp2);
 
        return (16.*M_PI*s_W2*c_W2)*(((delPiZZ_MZ.real()-delPiZZ_0.real())/MZ2)-(c_W2-s_W2)*delPiZAprime_0.real()/(sqrt(s_W2)*sqrt(c_W2))-delPiAAprime_0.real());
 
    }
    else
    {
        return std::numeric_limits<double>::quiet_NaN();
    }
    }

---

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

• THDMWSTU.h
• THDMWSTU.cpp