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} \)

virtual std::vector< WilsonCoefficient > &	CMdbs2 ()
	\( \Delta B = 2 \), \( B_{s} \)

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 > &	CMBMll (QCD::lepton lepton)=0

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

◆ 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, NLO_QED11),
mckcc(10, NDR, NLO, NLO_QED11),
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(2, NDR, NLO, NLO_QED11),
mckmm(1, NDR, NLO),
mcbsnn(2, NDR, NLO, NLO_QED11),
mcbdnn(2, NDR, NLO, NLO_QED11),
mcbsmm(8, NDR, NNLO, NLO_QED22),
mcbdmm(8, NDR, NNLO, NLO_QED22),
mculeptonnu(5, NDR, LO),
mcDLij(2, NDR, LO),
mcDLi3j(20, NDR, LO),
mcmueconv(8, NDR, LO),
mcgminus2mu(2, NDR, NLO),
mcC(2, NDR, QCD2, QED2), // blocks have to have the same orders
mcP(4, NDR, QCD2, QED2),
mcM(2, NDR, QCD2, QED2),
mcL(2, NDR, QCD2, QED2),
mcQ(4, NDR, QCD2, QED2),
mcB(1, NDR, QCD2, QED2),
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 107 of file StandardModelMatching.cpp.

107 {

108}

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, NPSMEFTd6GeneralMatching, SUSYMatching, and SUSYMassInsertionMatching.

Definition at line 966 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("StandardModelMatching::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 & 2
            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");
    }
 
 
    vmcdb.push_back(mcdbd2);
    return (vmcdb);
}

◆ CMdbs2()

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

virtual

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

Returns: return the vector of SM Wilson coefficients

Implements ModelMatching.

Reimplemented in FlavourWilsonCoefficient_DF2Matching, LoopMediatorsMatching, GeneralTHDMMatching, NPSMEFTd6GeneralMatching, SUSYMatching, SUSYMassInsertionMatching, THDMMatching, and THDMWMatching.

Definition at line 1016 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("StandardModelMatching::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");
    }
 
    vmcds.push_back(mcdbs2);
    return (vmcds);
}

◆ CMdd2()

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

virtual

\( \Delta C = 2 \),

Returns: return the vector of SM Wilson coefficients

Implements ModelMatching.

Reimplemented in FlavourWilsonCoefficient_DF2Matching, NPSMEFTd6GeneralMatching, SUSYMatching, and SUSYMassInsertionMatching.

Definition at line 1169 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("StandardModelMatching::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");
    }
 
    vmcd2.push_back(mcdd2);
    return (vmcd2);
}

◆ updateSMParameters()

void StandardModelMatching::updateSMParameters ( )

Updates to new Standard Model parameter sets.

Definition at line 110 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(), QCD::TOP, FULLNNLO);
    Mt_mut = SM.Mrun(Mut, SM.getQuarks(QCD::TOP).getMass_scale(),
            SM.getQuarks(QCD::TOP).getMass(), QCD::TOP, FULLNNLO);
    alstilde = SM.Als(Muw, FULLNNNLO, true, 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:

Detailed Description

Public Member Functions

Constructor & Destructor Documentation

◆ StandardModelMatching()

◆ ~StandardModelMatching()

Member Function Documentation

◆ CMdbd2()

◆ CMdbs2()

◆ CMdd2()

◆ updateSMParameters()