THDMWMatching Class Reference

#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} \)
virtual std::vector< WilsonCoefficient > &	CMdd2 ()
	\( \Delta C = 2 \),
	StandardModelMatching (const StandardModel &SM_i)
void	updateSMParameters ()
	Updates to new Standard Model parameter sets.
virtual	~StandardModelMatching ()
Public Member Functions inherited from ModelMatching
virtual std::vector< WilsonCoefficient > &	CMbnlep (int a)=0
virtual std::vector< WilsonCoefficient > &	CMbnlepCC (const int a)=0
virtual std::vector< WilsonCoefficient > &	CMbsg ()=0
virtual std::vector< WilsonCoefficient > &	CMBXsnn (QCD::lepton lepton)=0
virtual std::vector< WilsonCoefficient > &	CMd1 ()=0
virtual std::vector< WilsonCoefficient > &	CMd1Buras ()=0
virtual std::vector< WilsonCoefficientNew > &	CMDF1 (std::string blocks, unsigned int nops)=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 [Cheng:2015lsa].

Definition at line 106 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 [Cheng:2015lsa]

Reimplemented from StandardModelMatching.

Definition at line 151 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(), QCD::BOTTOM, 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 [Cheng:2015lsa]

Definition at line 311 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 [Cheng:2015lsa]

Definition at line 303 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 [Cheng:2015lsa]

Definition at line 307 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 [Cheng:2015lsa]

Definition at line 274 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 [Cheng:2015lsa]

Definition at line 255 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 [Cheng:2015lsa]

Definition at line 268 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 [Cheng:2015lsa]

Definition at line 262 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 [Cheng:2015lsa]

Definition at line 298 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 [Cheng:2015lsa]

Definition at line 279 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 [Cheng:2015lsa]

Definition at line 292 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 [Cheng:2015lsa]

Definition at line 285 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 [Cheng:2015lsa].

Definition at line 123 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 [Cheng:2015lsa].

Definition at line 112 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:

• THDMWMatching.h
• THDMWMatching.cpp