AmpDK2 Class Reference

A class for calculating the amplitudes contributing to \(\epsilon_K\) and \(\Delta M_K\). More...

#include <AmpDK2.h>

Inheritance diagram for AmpDK2:

Detailed Description

A class for calculating the amplitudes contributing to \(\epsilon_K\) and \(\Delta M_K\).

Author

HEPfit Collaboration

Copyright

GNU General Public License

This class is used to compute the theoretical value of amplitudes contributing to \(\epsilon_K\) and \(\Delta M_K\). The hadronic matrix elements are defined for the operators \(O_1, \ldots, O_5\) in the chiral limit as can be found in ...

Definition at line 30 of file AmpDK2.h.

Public Member Functions
	AmpDK2 (const StandardModel &SM_i)
	Constructor.

Protected Member Functions
gslpp::complex	AmpDK (orders order)
	compute the amplitude for kaon oscillations
gslpp::complex	AmpDMKNP (orders order)
	compute the NP part of the amplitude for kaon oscillations

Private Attributes
const StandardModel &	mySM

Constructor & Destructor Documentation

AmpDK2()

AmpDK2::AmpDK2

(

const StandardModel &

SM_i

)

Constructor.

Parameters

[in]

SM_i

a reference to an object of type StandardModel

Definition at line 12 of file AmpDK2.cpp.

: mySM(SM_i) 
{
    mySM.initializeBParameter("BK");
}

Member Function Documentation

AmpDK()

gslpp::complex AmpDK2::AmpDK ( orders  order )

protected

compute the amplitude for kaon oscillations

Parameters

[in]

order

the QCD order of the computation

Definition at line 18 of file AmpDK2.cpp.

{
    if (mySM.getFlavour().getHDF2().getCoeffK().getOrder() < order % 3)
        throw std::runtime_error("AmpDK::computeThValue(): requires cofficient of order not computed"); 
 
#if SUSYFIT_DEBUG & 4
    std::cout << "alsmu" << mySM.Als(2.) << std::endl;
#endif
    
    gslpp::vector<gslpp::complex> ** allcoeff = mySM.getFlavour().ComputeCoeffK(
            mySM.getBK().getMu(),
            mySM.getBK().getScheme());
            
    gslpp::vector<double> me(mySM.getBK().getBpars());
#if SUSYFIT_DEBUG & 2
    std::cout << "B-parameter: " << me(0) << std::endl;
#endif
#if SUSYFIT_DEBUG & 4
    std::cout << "alsmu" << mySM.Als(2.) << std::endl;
#endif
 
    double MK = mySM.getMesons(QCD::K_0).getMass();
    double Ms = mySM.Mrun(mySM.getBK().getMu(),
                mySM.getQuarks(QCD::STRANGE).getMass_scale(),
                mySM.getQuarks(QCD::STRANGE).getMass(), QCD::STRANGE, FULLNNLO);
    double Md = mySM.Mrun(mySM.getBK().getMu(),
                mySM.getQuarks(QCD::DOWN).getMass_scale(),
                mySM.getQuarks(QCD::DOWN).getMass(), QCD::DOWN, FULLNNLO);
    double KK = MK/(Ms+Md)*MK/(Ms+Md);
    double FK = mySM.getMesons(QCD::K_0).getDecayconst();
    double mm = MK*FK*FK;
    KK *= mm;
    me(0) *= 1./3.*mm;
    me(1) *= -5./24.*KK;
    me(2) *= 1./24.*KK;
    me(3) *= 1./4.*KK;
    me(4) *= 1./12.*KK;
 
#if SUSYFIT_DEBUG & 2
    std::cout << "matrix element: " << me(0) << std::endl;
#endif
#if SUSYFIT_DEBUG & 4
    std::cout << "alsmu" << mySM.Als(2.) << std::endl;
#endif
    switch(order) {
        case FULLNLO:
           return(((*(allcoeff[LO]) + *(allcoeff[NLO])) * me) / HCUT * 1.e3); // ns^-1
        case LO:
            return(((*(allcoeff[LO])) * me) / HCUT * 1.e3); // ns^-1
        default:
            throw std::runtime_error("AmpDK2::AmpDK(): order not implemented"); 
    }
}

AmpDMKNP()

gslpp::complex AmpDK2::AmpDMKNP ( orders  order )

protected

compute the NP part of the amplitude for kaon oscillations

Parameters

[in]

order

the QCD order of the computation

Definition at line 72 of file AmpDK2.cpp.

{
    if (mySM.getFlavour().getHDF2().getCoeffmK().getOrder() < order % 3)
        throw std::runtime_error("AmpDMKNP::computeThValue(): requires cofficient of order not computed");
 
    gslpp::vector<gslpp::complex> ** allcoeff = mySM.getFlavour().ComputeCoeffmK(
            mySM.getBK().getMu(),
            mySM.getBK().getScheme());
 
    gslpp::vector<double> me(mySM.getBK().getBpars());
    double MK = mySM.getMesons(QCD::K_0).getMass();
    double Ms = mySM.Mrun(mySM.getBK().getMu(),
                mySM.getQuarks(QCD::STRANGE).getMass_scale(),
                mySM.getQuarks(QCD::STRANGE).getMass(), QCD::STRANGE, FULLNNLO);
    double Md = mySM.Mrun(mySM.getBK().getMu(),
                mySM.getQuarks(QCD::DOWN).getMass_scale(),
                mySM.getQuarks(QCD::DOWN).getMass(), QCD::DOWN, FULLNNLO);
    double KK = MK/(Ms+Md)*MK/(Ms+Md);
    double FK = mySM.getMesons(QCD::K_0).getDecayconst();
    double mm = MK*FK*FK;
    KK *= mm;
    me(0) *= 1./3.*mm;
    me(1) *= -5./24.*KK;
    me(2) *= 1./24.*KK;
    me(3) *= 1./4.*KK;
    me(4) *= 1./12.*KK;
 
    switch(order) {
        case FULLNLO:
            return(((*(allcoeff[LO]) + *(allcoeff[NLO])) * me) / HCUT * 1.e3); // ns^-1
        case LO:
            return(((*(allcoeff[LO])) * me) / HCUT * 1.e3); // ns^-1
        default:
            throw "AmpDM2::AmpDMKNP(): order not implemented";
    }
}

Member Data Documentation

mySM

const StandardModel& AmpDK2::mySM

private

Definition at line 55 of file AmpDK2.h.

---

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

• AmpDK2.h
• AmpDK2.cpp