#include <THDMWMatching.h>

Inheritance diagram for THDMWMatching:

Detailed Description

Author: HEPfit Collaboration

Copyright: GNU General Public License

Definition at line 26 of file THDMWMatching.h.

Public Member Functions
virtual gslpp::complex	C10NP (double xt, double xS, gslpp::complex etaU)

virtual std::vector< WilsonCoefficient > &	CMBMll (QCD::lepton lepton)

virtual std::vector< WilsonCoefficient > &	CMdbs2 ()

virtual gslpp::complex	CNPVLLcc (double xc, double xb, double xS, gslpp::complex etaU, gslpp::complex etaD)

virtual double	CNPVLLccEtadEtasu (double xc, double xb, double xS)

virtual double	CNPVLLccEtau2 (double xc, double xb, double xS)

virtual gslpp::complex	CNPVLLct (double xc, double xb, double xt, double xS, gslpp::complex etaU, gslpp::complex etaD)

virtual double	CNPVLLctEtadEtasu (double xc, double xb, double xt, double xS)

virtual double	CNPVLLctEtau2 (double xc, double xb, double xt, double xS)

virtual double	CNPVLLctEtau4 (double xc, double xb, double xt, double xS)

virtual gslpp::complex	CNPVLLtt (double xc, double xb, double xt, double xS, gslpp::complex etaU, gslpp::complex etaD)

virtual double	CNPVLLttEtadEtasu (double xc, double xb, double xt, double xS)

virtual double	CNPVLLttEtau2 (double xc, double xb, double xt, double xS)

virtual double	CNPVLLttEtau4 (double xc, double xb, double xt, double xS)

virtual gslpp::complex	CPNP (double xt, double xS, gslpp::complex etaU, gslpp::complex etaD)

virtual gslpp::complex	CSNP (double nu1, double xh, double xt, double xS, gslpp::complex etaU, gslpp::complex etaD)

	THDMWMatching (const THDMW &THDMW_i)

Public Member Functions inherited from StandardModelMatching
virtual std::vector< WilsonCoefficient > &	CMdbd2 ()
	\( \Delta B = 2 \), \( B_{d} \) More...

virtual std::vector< WilsonCoefficient > &	CMdd2 ()
	\( \Delta C = 2 \), More...

	StandardModelMatching (const StandardModel &SM_i)

void	updateSMParameters ()
	Updates to new Standard Model parameter sets. More...

virtual	~StandardModelMatching ()

Public Member Functions inherited from ModelMatching
virtual std::vector< WilsonCoefficient > &	CMbnlep (const int a)=0

virtual std::vector< WilsonCoefficient > &	CMbnlepCC (const int a)=0

virtual std::vector< WilsonCoefficient > &	CMbsg ()=0

virtual std::vector< WilsonCoefficient > &	CMd1 ()=0

virtual std::vector< WilsonCoefficient > &	CMd1Buras ()=0

virtual std::vector< WilsonCoefficient > &	CMprimeBMll (QCD::lepton lepton)=0

virtual std::vector< WilsonCoefficient > &	CMprimebsg ()=0

virtual	~ModelMatching ()

Constructor & Destructor Documentation

◆ THDMWMatching()

THDMWMatching::THDMWMatching ( const THDMW & THDMW_i )

Definition at line 16 of file THDMWMatching.cpp.

                                                   :
  
     StandardModelMatching(THDMW_i),
     myTHDMW(THDMW_i),
     myCKM(3, 3, 0.),
     mcBMll(13, NDR, NLO),
     mcbsg(8, NDR, NNLO),
     mcdbs2(5, NDR, NLO)
 {
 }

Member Function Documentation

◆ C10NP()

gslpp::complex THDMWMatching::C10NP	(	double	xt,
		double	xS,
		gslpp::complex	etaU
	)

virtual

Returns: C10 Wilson coefficient of \( B_q \to l \bar{l}\) according to [53].

Definition at line 83 of file THDMWMatching.cpp.

                                                                           {
     
     gslpp::complex C10 = etaU.abs2()*(pow(xt,2)/6)*(1/(xS-xt)+(xS/pow(xS-xt,2))*(log(xt)-log(xS)));
     return C10;
     }

◆ CMBMll()

std::vector< WilsonCoefficient > & THDMWMatching::CMBMll ( QCD::lepton lepton )

virtual

Implements ModelMatching.

Definition at line 31 of file THDMWMatching.cpp.

 {
     //From 1504.00839 Note that in this paper the Effective Hamiltonian has a typo and a factor 1/sw^2 it's missing
     double Mut = myTHDMW.getMut();
     double Muw = myTHDMW.getMuw();
     double MW = myTHDMW.Mw();
     double xt = x_t(Mut);
     double mhsq = myTHDMW.getMyTHDMWCache()->mhsq;
     double xh = mhsq/pow(MW,2);
     double sW2 = myTHDMW.sW2();
     double mSp2 =myTHDMW.getMyTHDMWCache()->mSpsq;
     double xS = mSp2/pow(MW,2);
     double nu1 = myTHDMW.getTHDMW_nu1();
     gslpp::complex etaU = myTHDMW.getTHDMW_etaU();
     gslpp::complex etaD = myTHDMW.getTHDMW_etaD();
     
     double mb=myTHDMW.getQuarks(QCD::BOTTOM).getMass();
    
     double ml=myTHDMW.getLeptons(lepton).getMass();
     
     
     
     vmcBMll = StandardModelMatching::CMBMll(lepton);
       switch (mcbsg.getScheme()) {
         case NDR:
  
             break;
         default:
         std::stringstream out;
        out << mcBMll.getScheme();
      throw std::runtime_error("THDMWMatching::CMBMll(): scheme " + out.str() + "not implemented");
     }
     mcBMll.setMu(Muw);
     
      switch (mcBMll.getOrder()) {
         case NNLO:
         case NLO:           
         case LO:           
             mcBMll.setCoeff(9 , C10NP(xt, xS, etaU)/(sW2), LO);
             mcBMll.setCoeff(10 , CSNP(nu1, xh, xt, xS, etaU, etaD)*mb*ml/(MW*MW*sW2), LO);
             mcBMll.setCoeff(11 , CPNP(xt, xS, etaU, etaD)*(mb*ml)/(MW*MW*sW2), LO);
             break;
         default:
             std::stringstream out;
             out << mcBMll.getOrder();
             throw std::runtime_error("THDMWMatching::CMBMll(): order " + out.str() + "not implemented"); 
             }
     vmcBMll.push_back(mcBMll);
     return (vmcBMll);
 }

◆ CMdbs2()

std::vector< WilsonCoefficient > & THDMWMatching::CMdbs2 ( )

virtual

Returns: THDMW Wilson coefficients for \( B_s \to \bar{B_s}\) according to [53]

Reimplemented from StandardModelMatching.

Definition at line 118 of file THDMWMatching.cpp.

                                                     {
  
     double Mut = myTHDMW.getMut();
     double Muc = myTHDMW.getMuc();
     double Mub = myTHDMW.getMub();
     double xt = x_t(Mut);//ratio mass of the top square at the scale of its mass over mass of the W square. Why not Mw scale???
     double xc= x_c(Muc);
     double mb =myTHDMW.Mrun(Mub, myTHDMW.getQuarks(QCD::BOTTOM).getMass_scale(),
                               myTHDMW.getQuarks(QCD::BOTTOM).getMass(), FULLNNLO);
     double mSp2=myTHDMW.getMyTHDMWCache()->mSpsq;
     double MW=myTHDMW.Mw();
     double xb=pow(mb,2)/(MW*MW);
     
     double xS=mSp2/(MW*MW);
     gslpp::complex etaU = myTHDMW.getTHDMW_etaU();
     gslpp::complex etaD = myTHDMW.getTHDMW_etaD();
     double GF=myTHDMW.getGF();
     double co = GF / 4. / M_PI * MW;
  
  
     vmcds = StandardModelMatching::CMdbs2();
     mcdbs2.setMu(Mut);
  
     switch (mcdbs2.getOrder()) {
         case NNLO:
         case NLO:
         case LO:
             mcdbs2.setCoeff(0, co * co * CNPVLL(xc, xb, xt, xS, etaU, etaD).conjugate(), LO);
             break;
         default:
             std::stringstream out;
             out << mcdbs2.getOrder();
             throw std::runtime_error("THDMMatching::CMdbs2(): order " + out.str() + "not implemented");
     }
  
     vmcds.push_back(mcdbs2);
     //The following are the primed coefficients.
     mcdbs2.setMu(Mut);
  
     switch (mcdbs2.getOrder()) {
         case NNLO:
         case NLO:
         case LO:
             mcdbs2.setCoeff(1, co * co * CNPSRR1(xc, xb, xt, xS, etaU, etaD).conjugate()+4*co * co * CNPSRR2(xc, xb, xt, xS, etaU, etaD).conjugate(), LO);
             mcdbs2.setCoeff(2, 8*co * co * CNPSRR2(xc, xb, xt, xS, etaU, etaD).conjugate(), LO);
             //OSRR2=4*O2+8*O3 -> From Buras basis to SUSY basis
             break;
         default:
             std::stringstream out;
             out << mcdbs2.getOrder();
             throw std::runtime_error("THDMMatching::CMdbs2(): order " + out.str() + "not implemented");
     }
  
     vmcds.push_back(mcdbs2);
  
     return(vmcds);
 }

◆ CNPVLLcc()

gslpp::complex THDMWMatching::CNPVLLcc	(	double	xc,
		double	xb,
		double	xS,
		gslpp::complex	etaU,
		gslpp::complex	etaD
	)

virtual

Returns: \( C^{NP,cc}_{VLL} \) according to [53]

Definition at line 278 of file THDMWMatching.cpp.

                                                                                                           {
       return etaD*etaU.conjugate()*CNPVLLccEtadEtasu(xc, xb, xS) +
               etaU.abs2()*CNPVLLccEtau2(xc, xb ,xS);
   }

◆ CNPVLLccEtadEtasu()

double THDMWMatching::CNPVLLccEtadEtasu	(	double	xc,
		double	xb,
		double	xS
	)

virtual

Returns: \( C^{NP,cc}_{VLL, \eta_d \eta^*_u} \) according to [53]

Definition at line 270 of file THDMWMatching.cpp.

                                                                       {
       return (-2*xb*xc)/(3.*(-xS + 2*xc*(1 + xS)));
   }

◆ CNPVLLccEtau2()

double THDMWMatching::CNPVLLccEtau2	(	double	xc,
		double	xb,
		double	xS
	)

virtual

Returns: \( C^{NP,cc}_{VLL, \eta_u^2} \) according to [53]

Definition at line 274 of file THDMWMatching.cpp.

                                                                   {
       return (-4*xb*xc)/(9.*(-xS + 3*xc*(1 + xS)));
   }

◆ CNPVLLct()

gslpp::complex THDMWMatching::CNPVLLct	(	double	xc,
		double	xb,
		double	xt,
		double	xS,
		gslpp::complex	etaU,
		gslpp::complex	etaD
	)

virtual

Returns: \( C^{NP,ct}_{VLL} \) according to [53]

Definition at line 241 of file THDMWMatching.cpp.

                                                                                                                      {
       return etaD*etaU.conjugate()*CNPVLLctEtadEtasu(xc ,xb ,xt , xS)+pow(etaU.abs2(),2)*CNPVLLctEtau4(xc ,xb ,xt , xS) +
               etaU.abs2()*CNPVLLctEtau2(xc ,xb ,xt , xS);
   }

◆ CNPVLLctEtadEtasu()

double THDMWMatching::CNPVLLctEtadEtasu	(	double	xc,
		double	xb,
		double	xt,
		double	xS
	)

virtual

Returns: \( C^{NP,ct}_{VLL, \eta_d \eta^*_u} \) according to [53]

Definition at line 222 of file THDMWMatching.cpp.

                                                                                  {
       return (-2*xb*xc*xt*(-1 - pow(xS,2) + xt + xS*xt))/(3.*(-1 + xS)*(-xS + xc*(1 + xS))*(xS - xt)*(-1 + xt)) + 
    (2*xb*xc*xt*log(xc))/(3.*(-(xS*xt) + xc*(2*xt + xS*(1 + 2*xt)))) + 
    (2*xb*xc*xt*(-2*pow(xS,2) + xt + xS*xt)*log(xS))/(3.*pow(-1 + xS,2)*(-2*xc + xS)*pow(xS - xt,2)) + 
    (2*xb*xc*xt*(xS + (-2 + xt)*xt)*log(xt))/(3.*pow(xS - xt,2)*pow(-1 + xt,2)*(-xc + xt));
   }

◆ CNPVLLctEtau2()

double THDMWMatching::CNPVLLctEtau2	(	double	xc,
		double	xb,
		double	xt,
		double	xS
	)

virtual

Returns: \( C^{NP,ct}_{VLL, \eta_u^2} \) according to [53]

Definition at line 235 of file THDMWMatching.cpp.

                                                                              {
       return -((f10(xt,xS) + f11(xt,xS) + f12(xt,xS))*f9(xc,xb,xt,xS)) + 
    f1(xc,xt,xS)*(f2(xb,xt,xS) + f3(xb,xt,xS) + f4(xb,xt,xS))*log(xS) + 
    f5(xc,xt,xS)*(f6(xb,xt,xS) + f7(xb,xt,xS) + f8(xb,xt))*log(xt);
   }

◆ CNPVLLctEtau4()

double THDMWMatching::CNPVLLctEtau4	(	double	xc,
		double	xb,
		double	xt,
		double	xS
	)

virtual

Returns: \( C^{NP,ct}_{VLL, \eta_u^4} \) according to [53]

Definition at line 229 of file THDMWMatching.cpp.

                                                                              {
       return (11*xc*pow(xt,2)*(-18*xS*pow(xS - xt,2) + xb*(5*pow(xS,2) - 22*xS*xt + 5*pow(xt,2))))/(324.*pow(xS - xt,3)*(-(xS*xt) + xc*(2*xS + 3*xt))) - 
    (11*xc*xS*xt*(xb*xS*(xS - 3*xt) - 3*pow(xS - xt,2)*xt)*log(xS))/(54.*(4*xc - xS)*pow(xS - xt,4)) + 
    (11*xc*pow(xt,2)*(xb*xS*(xS - 3*xt) - 3*pow(xS - xt,2)*xt)*log(xt))/(54.*(3*xc - xt)*pow(xS - xt,4));
   }

◆ CNPVLLtt()

gslpp::complex THDMWMatching::CNPVLLtt	(	double	xc,
		double	xb,
		double	xt,
		double	xS,
		gslpp::complex	etaU,
		gslpp::complex	etaD
	)

virtual

Returns: \( C^{NP,tt}_{VLL} \) according to [53]

Definition at line 265 of file THDMWMatching.cpp.

                                                                                                                      {
       return etaD*etaU.conjugate()*CNPVLLttEtadEtasu(xc ,xb ,xt , xS)+pow(etaU.abs2(),2)*CNPVLLttEtau4(xc ,xb ,xt , xS) +
               etaU.abs2()*CNPVLLttEtau2(xc ,xb ,xt , xS);
   }

◆ CNPVLLttEtadEtasu()

double THDMWMatching::CNPVLLttEtadEtasu	(	double	xc,
		double	xb,
		double	xt,
		double	xS
	)

virtual

Returns: \( C^{NP,tt}_{VLL, \eta_d \eta^*_u} \) according to [53]

Definition at line 246 of file THDMWMatching.cpp.

                                                                                  {
       return (2*xb*xt*(pow(xS,2)*(1 + xt) - xS*pow(1 + xt,2) + xt*(3 - 3*xt + 2*pow(xt,2))))/(3.*(-1 + xS)*pow(xS - xt,2)*pow(-1 + xt,2)) + 
    (4*xb*pow(xt,2)*(-pow(xS,2) + xt)*log(xS))/(3.*pow(-1 + xS,2)*pow(xS - xt,3)) + 
    (4*xb*pow(xt,2)*(-pow(xS,2) + xt + 2*xS*xt - 3*pow(xt,2) + pow(xt,3))*log(xt))/(3.*pow(xS - xt,3)*pow(-1 + xt,3));
   }

◆ CNPVLLttEtau2()

double THDMWMatching::CNPVLLttEtau2	(	double	xc,
		double	xb,
		double	xt,
		double	xS
	)

virtual

Returns: \( C^{NP,tt}_{VLL, \eta_u^2} \) according to [53]

Definition at line 259 of file THDMWMatching.cpp.

                                                                              {
       return  f24(xt,xS)*(f30(xb,xt,xS) + f31(xb,xt,xS)) + 
    f13(xt,xS)*(f14(xb,xt,xS) + f15(xb,xt,xS) + f16(xb,xt,xS))*log(xS) - 
    f17(xt,xS)*(f22(xb,xt,xS) + f23(xb,xt,xS))*log(xt);
   }

◆ CNPVLLttEtau4()

double THDMWMatching::CNPVLLttEtau4	(	double	xc,
		double	xb,
		double	xt,
		double	xS
	)

virtual

Returns: \( C^{NP,tt}_{VLL, \eta_u^4} \) according to [53]

Definition at line 252 of file THDMWMatching.cpp.

                                                                              {
       return (-11*pow(xt,2)*(-9*pow(xS - xt,2)*(xS + xt) + xb*(5*pow(xS,2) - 22*xS*xt + 5*pow(xt,2))))/(162.*pow(xS - xt,4)) + (11*pow(xt,2)*(-6*xS*pow(xS - xt,2)*xt + 
         xb*(pow(xS,3) - 3*pow(xS,2)*xt - 3*xS*pow(xt,2) + pow(xt,3)))*log(xS))/(54.*pow(xS - xt,5)) - 
    (11*pow(xt,2)*(-6*xS*pow(xS - xt,2)*xt + 
         xb*(pow(xS,3) - 3*pow(xS,2)*xt - 3*xS*pow(xt,2) + pow(xt,3)))*log(xt))/(54.*pow(xS - xt,5));
   }

◆ CPNP()

gslpp::complex THDMWMatching::CPNP	(	double	xt,
		double	xS,
		gslpp::complex	etaU,
		gslpp::complex	etaD
	)

virtual

Returns: CP Wilson coefficient of \( B_q \to l \bar{l}\) according to [53].

Definition at line 100 of file THDMWMatching.cpp.

                                                                                            {
     double sW2 = myTHDMW.sW2();
     gslpp::complex CP =(xt/(3*pow(xS-xt,2)))*(etaD*etaU.conjugate()*(-(xt+xS)/2+xt*xS*(log(xS)-log(xt))/(xS-xt))+
     (etaU.abs2()/(6*(xS-xt)))*((pow(xS,2)-8*xS*xt-17*pow(xt,2))/6+pow(xt,2)*(3*xS+xt)*(log(xS)-log(xt))/(xS-xt)))+
     (2*sW2*xt/(9*pow(xS-xt,2)))*(etaD*etaU.conjugate()*((5*xt-3*xS)/2+xS*(2*xS-3*xt)*(log(xS)-log(xt))/(xS-xt))
     +(etaU.abs2()/(6*(xS-xt)))*((17*pow(xS,2)-64*xS*xt+71*pow(xt,2))/6)
     -(4*pow(xS,3)-12*pow(xS,2)*xt+9*xS*pow(xt,2)+3*pow(xt,3))*(log(xS)-log(xt))/(xS-xt))
     +etaU.abs2()*(1-sW2)*pow(xt,2)*(xS*(log(xS)-log(xt))+xt-xS)/(3*pow(xS-xt,2));
     return CP;
     }

◆ CSNP()

gslpp::complex THDMWMatching::CSNP	(	double	nu1,
		double	xh,
		double	xt,
		double	xS,
		gslpp::complex	etaU,
		gslpp::complex	etaD
	)

virtual

Returns: CS Wilson coefficient of \( B_q \to l \bar{l}\) according to [53].

Definition at line 89 of file THDMWMatching.cpp.

                                                                                                                  {
     double mhsq= myTHDMW.getMyTHDMWCache()->mhsq;
     double vev= myTHDMW.v();
     
     gslpp::complex CS =(pow(vev,2)*xt*nu1/(mhsq*6*(xS-xt)))*(etaD*etaU.conjugate()*(xt*(log(xS)-log(xt))/(xS-xt)-1)+
             etaU.abs2()*(pow(xt,2)*(log(xS)-log(xt))/(2*pow(xS-xt,2))+(xS-3*xt)/(4*(xS-xt))))+
             (2*xt/(3*xh))*(etaD*etaU.conjugate()*xt*(1-xS*(log(xS)-log(xt))/(xS-xt))/(xS-xt)+
             etaU.abs2()*(xt/(2*pow(xS-xt,2)))*((xS+xt)/2-xS*xt*(log(xS)-log(xt))/(xS-xt)));
     return CS;
     }

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

Detailed Description

Public Member Functions

Constructor & Destructor Documentation

◆ THDMWMatching()

Member Function Documentation

◆ C10NP()

◆ CMBMll()

◆ CMdbs2()

◆ CNPVLLcc()

◆ CNPVLLccEtadEtasu()

◆ CNPVLLccEtau2()

◆ CNPVLLct()

◆ CNPVLLctEtadEtasu()

◆ CNPVLLctEtau2()

◆ CNPVLLctEtau4()

◆ CNPVLLtt()

◆ CNPVLLttEtadEtasu()

◆ CNPVLLttEtau2()

◆ CNPVLLttEtau4()

◆ CPNP()

◆ CSNP()