StandardModelMatching Class Reference

A class for the matching in the Standard Model. More...

#include <StandardModelMatching.h>

Inheritance diagram for StandardModelMatching:

Detailed Description

A class for the matching in the Standard Model.

Author

HEPfit Collaboration

Copyright

GNU General Public License

Definition at line 26 of file StandardModelMatching.h.

Public Member Functions
virtual std::vector< WilsonCoefficient > &	CMdbd2 ()
	\( \Delta B = 2 \), \( B_{d} \) More...
virtual std::vector< WilsonCoefficient > &	CMdbs2 ()
	\( \Delta B = 2 \), \( B_{s} \) 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 > &	CMBMll (QCD::lepton lepton)=0
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

StandardModelMatching()

StandardModelMatching::StandardModelMatching

(

const StandardModel &

SM_i

)

Definition at line 16 of file StandardModelMatching.cpp.

: SM(SM_i),
        mcdbd2(5, NDR, NLO),
        mcdbs2(5, NDR, NLO),
        mcdd2(5, NDR, NLO),
        mcdk2(5, NDR, NLO),
        mck(10, NDR, NLO),
        mckcc(10, NDR, NLO),
        mcbsg(8, NDR, NNLO),
        mcprimebsg(8, NDR, NNLO),
        mcBMll(13, NDR, NLO),
        mcprimeBMll(13, NDR, NLO),
        mcbnlep(10, NDR, NLO, NLO_QED11),
        mcbnlepCC(10, NDR, NLO),
        mcd1(10, NDR, NLO),
        mcd1Buras(10, NDR, NLO),
        mckpnn(1, NDR, NLO, NLO_QED11),
        mckmm(1, NDR, NLO),
        mcbsnn(1, NDR, NLO),
        mcbdnn(1, NDR, NLO),
        mcbsmm(8, NDR, NNLO, NLO_QED22),
        mcbdmm(8, NDR, NNLO, NLO_QED22),
        mcbtaunu(3, NDR, LO),
        mcDLij(2, NDR, LO),
        mcDLi3j(20, NDR, LO),
        mcmueconv(8, NDR, LO),
        mcgminus2mu(2, NDR, NLO),
        Vckm(3,3,0.)
{    
    swa = 0.;
    swb = 0.;
    swc = 0.;
    xcachea = 0.;
    xcacheb = 0.;
    xcachec = 0.;
    
 
    for (int j=0; j<10; j++) {
        CWD1ArrayLO[j] = 0.; 
        CWD1ArrayNLO[j] = 0.;
        CWbnlepArrayLOqcd[j] = 0.;
        CWbnlepArrayNLOqcd[j] = 0.;
        CWbnlepArrayLOew[j] = 0.;
        CWbnlepArrayNLOew[j] = 0.;
    };
    
    
    for(int j=0; j<19; j++){
        CWBMllArrayLO[j] = 0.;
        CWBMllArrayNLO[j] = 0.;
    }
    
    for(int j=0; j<8; j++){
        CWbsgArrayLO[j] = 0.;
        CWbsgArrayNLO[j] = 0.;
        CWbsgArrayNNLO[j] = 0.;
        CWprimebsgArrayLO[j] = 0.;
        CWprimebsgArrayNLO[j] = 0.;
    }
    
    for(int j=0; j<8; j++){
        CWBsmmArrayNNLOqcd[j] = 0.; 
        CWBsmmArrayNLOqcd[j] = 0.;
        CWBsmmArrayLOqcd[j] = 0.;
        CWBsmmArrayNLOewt4[j] = 0.;
        CWBsmmArrayNLOewt2[j] = 0.;
        CWBsmmArrayNLOew[j] = 0.;
    }
    
    for(int j=0; j<8; j++){
        CWBdmmArrayNNLOqcd[j] = 0.; 
        CWBdmmArrayNLOqcd[j] = 0.;
        CWBdmmArrayLOqcd[j] = 0.;
        CWBdmmArrayNLOewt4[j] = 0.;
        CWBdmmArrayNLOewt2[j] = 0.;
        CWBdmmArrayNLOew[j] = 0.;
    }
    
    
    Nc = SM.getNc();
    CF = SM.getCF();
    gamma0 = 6. * (Nc - 1.) / Nc;
    J5 = SM.Beta1(5) * gamma0 / 2. / SM.Beta0(5) / SM.Beta0(5) - ((Nc - 1.)/(2. * Nc) * (-21. + 57./Nc - 19./3. * Nc + 20./3.)) / 2. / SM.Beta0(5);
    BtNDR = 5. * (Nc - 1.) / 2. / Nc + 3. * CF;
}

~StandardModelMatching()

StandardModelMatching::~StandardModelMatching ( )

virtual

Definition at line 102 of file StandardModelMatching.cpp.

{}

Member Function Documentation

CMdbd2()

std::vector< WilsonCoefficient > & StandardModelMatching::CMdbd2 ( )

virtual

\( \Delta B = 2 \), \( B_{d} \)

Returns

return the vector of SM Wilson coefficients

Implements ModelMatching.

Reimplemented in FlavourWilsonCoefficient_DF2Matching.

Definition at line 993 of file StandardModelMatching.cpp.

{
    double gammam = 6. * CF;
    double Bt;  
    
    
    double xt = x_t(Mut);
    gslpp::complex co = GF / 4. / M_PI * Mw * SM.getCKM().computelamt_d();
    
    vmcdb.clear();
 
    switch (mcdbd2.getScheme()) {
        case NDR:
            Bt = BtNDR;
            break;
        case HV:
        case LRI:
        default:
            std::stringstream out;
            out << mcdbd2.getScheme();
            throw std::runtime_error("StandardModel::CMdb2(): scheme " + out.str() + "not implemented"); 
    }
 
    mcdbd2.setMu(Mut);
    
    switch (mcdbd2.getOrder()) {
        case NNLO:
        case NLO:
            mcdbd2.setCoeff(0, co * co * 4. * (SM.Als(Mut, FULLNLO) / 4. / M_PI * (S1(xt) + //* CHECK ORDER *//
                    (Bt + gamma0 * log(Mut / Mw)) * S0(xt, xt) + 2. * gammam * S0p(xt) * log(Mut / Mw))), NLO);
#if SUSYFIT_DEBUG & 1
    std::cout << "Mw = " << Mw << " mt(mut=" << Mut << ")= " << Mt_mut << " xt(mut=" << Mut << ")= " << xt << "matching of DB=2: S0(xt) = " << S0(xt) << 
                ", S1(xt) = " << S1(xt) +
                    (Bt + gamma0 * log(Muw / Mw)) * S0(xt, xt) + 2. * gammam * S0p(xt) * log(Muw / Mw) 
            << ", lambdat_d^2 = " << SM.getCKM().computelamt_d()*SM.getCKM().computelamt_d() << std::endl;
#endif
        case LO:
            mcdbd2.setCoeff(0, co * co * 4. * (S0(xt, xt)), LO);
            break;
        default:
            std::stringstream out;
            out << mcdbd2.getOrder();
            throw std::runtime_error("StandardModelMatching::CMdbd2(): order " + out.str() + "not implemented"); 
    }
    
 

CMdbs2()

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

virtual

\( \Delta B = 2 \), \( B_{s} \)

Returns

return the vector of SM Wilson coefficients

Implements ModelMatching.

Reimplemented in GeneralTHDMMatching, THDMWMatching, FlavourWilsonCoefficient_DF2Matching, and THDMMatching.

Definition at line 1043 of file StandardModelMatching.cpp.

{
    double gammam = 6. * CF;
    double Bt;
    double xt = x_t(Mut);
    gslpp::complex co = GF / 4. / M_PI * Mw * SM.getCKM().computelamt_s();
 
    vmcds.clear();
 
    switch (mcdbs2.getScheme()) {
        case NDR:
            Bt = BtNDR;
            break;
        case HV:
        case LRI:
        default:
            std::stringstream out;
            out << mcdbs2.getScheme();
            throw std::runtime_error("StandardModel::CMdbs2(): scheme " + out.str() + "not implemented"); 
    }
 
    mcdbs2.setMu(Mut);
 
    switch (mcdbs2.getOrder()) {
        case NNLO:
        case NLO:          
            mcdbs2.setCoeff(0, co * co * 4. * (SM.Als(Mut, FULLNLO) / 4. / M_PI * (S1(xt) + //* CHECK ORDER *//
                    (Bt + gamma0 * log(Mut / Mw)) * S0(xt, xt) + 2. * gammam * S0p(xt) * log(Mut / Mw))), NLO);
         case LO:
            mcdbs2.setCoeff(0, co * co * 4. * (S0(xt, xt)), LO);
            break;
        default:
            std::stringstream out;
            out << mcdbs2.getOrder();
            throw std::runtime_error("StandardModelMatching::CMdbs2(): order " + out.str() + "not implemented"); 
    }    
 

CMdd2()

std::vector< WilsonCoefficient > & StandardModelMatching::CMdd2 ( )

virtual

\( \Delta C = 2 \),

Returns

return the vector of SM Wilson coefficients

Implements ModelMatching.

Reimplemented in FlavourWilsonCoefficient_DF2Matching.

Definition at line 1199 of file StandardModelMatching.cpp.

{
    vmcd2.clear();
 
    switch (mcdd2.getScheme()) {
        case NDR:
        case HV:
        case LRI:
          break; 
        default:
            std::stringstream out;
            out << mcdd2.getScheme();
            throw std::runtime_error("StandardModel::CMdd2(): scheme " + out.str() + "not implemented"); 
    }
 
    mcdd2.setMu(Muw);
 
    switch (mcdd2.getOrder()) {
        case NNLO:
        case NLO:
            for(int i=0; i<5; i++)
                mcdd2.setCoeff(i, 0., NLO);
        case LO:
            for(int j=0; j<5; j++)
                mcdd2.setCoeff(j, 0., LO);
            break;
        default:
            std::stringstream out;
            out << mcdd2.getOrder();
            throw std::runtime_error("StandardModelMatching::CMdd2(): order " + out.str() + "not implemented"); 
    }

updateSMParameters()

void StandardModelMatching::updateSMParameters

(

)

Updates to new Standard Model parameter sets.

Definition at line 105 of file StandardModelMatching.cpp.

{
    Mut = SM.getMut();
    Muw = SM.getMuw(); 
    Ale = SM.getAle();
    Mt_muw = SM.Mrun(Muw, SM.getQuarks(QCD::TOP).getMass_scale(), 
                        SM.getQuarks(QCD::TOP).getMass(), FULLNNLO);
    Mt_mut = SM.Mrun(Mut, SM.getQuarks(QCD::TOP).getMass_scale(), 
                        SM.getQuarks(QCD::TOP).getMass(), FULLNNLO);
    alstilde = SM.Als(Muw, FULLNNNLO, true) / 4. / M_PI; //  SM.Alstilde5(Muw) WHICH ONE TO USE?
    aletilde = SM.Ale(Muw, FULLNLO) / 4. / M_PI; // Ale / 4. / M_PI; // WHERE IS ale(mu)?
    GF = SM.getGF();
    Mw_tree = SM.Mw_tree();
    Mw = SM.Mw(); /* on-shell Mw */
    sW2 = SM.sW2(); /* on-shell sW2 */
    /* NP models should be added here after writing codes for Mw. */
    /* The following is commented out till further review.*/
//    if (SM.getModelName()=="StandardModel") {
//        Mw = SM.Mw(); /* on-shell Mw */
//        sW2 = SM.sW2(); /* on-shell sW2 */
//    } else {
//        Mw = Mw_tree;
//        sW2 = 1.0 - Mw*Mw/SM.getMz()/SM.getMz();
//    }
//    sW2 = (M_PI * Ale ) / ( sqrt(2.) * GF * Mw * Mw ); // WARNING: only for 
    Vckm = SM.getVCKM();
    lam_t = SM.getCKM().computelamt();
    L = 2*log(Muw/Mw);
    Lz = 2*log(Muw/SM.getMz());
    mu_b = SM.getMub();
}

---

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

• StandardModelMatching.h
• StandardModelMatching.cpp