AmpDS1 Class Reference

#include <AmpDS1.h>

Inheritance diagram for AmpDS1:

Detailed Description

Definition at line 16 of file AmpDS1.h.

Public Member Functions
	AmpDS1 (const StandardModel &SM_i)
	compute the amplitude for \( K_L \) decay in 2 pion More...

Protected Member Functions
gslpp::complex	AmpDS1pp0 (orders order)
gslpp::complex	AmpDS1pp2 (orders order)

Private Attributes
const StandardModel &	mySM

Constructor & Destructor Documentation

AmpDS1()

AmpDS1::AmpDS1

(

const StandardModel &

SM_i

)

compute the amplitude for \( K_L \) decay in 2 pion

Parameters

Flavour

Definition at line 13 of file AmpDS1.cpp.

: mySM(SM_i) 
{
    mySM.initializeBParameter("BKd1");
    mySM.initializeBParameter("BKd3");
}

Member Function Documentation

AmpDS1pp0()

gslpp::complex AmpDS1::AmpDS1pp0 ( orders  order )

protected

Parameters

order

Returns

the amplitude for \( K_L \) decay in 2 pion with 0 isospin change

Definition at line 20 of file AmpDS1.cpp.

{
    if (mySM.getFlavour().getHDS1().getCoeffDS1PP().getOrder() < order){
        std::stringstream out;
        out << order;
        throw std::runtime_error("AmpDK1::computeThValue(): requires cofficient of order" 
                                 + out.str() + "not computed");
    }
 
    gslpp::vector<gslpp::complex> ** allcoeff = mySM.getFlavour().ComputeCoeffDS1PP(
            mySM.getBKd1().getMu(),
            mySM.getBKd1().getScheme());
 
    gslpp::vector<double> me1(mySM.getBKd1().getBpars());
    
    double MK = mySM.getMesons(QCD::K_0).getMass();
    double MP = mySM.getMesons(QCD::P_0).getMass();
    double FK = mySM.getMesons(QCD::K_0).getDecayconst();
    double FP = mySM.getMesons(QCD::P_0).getDecayconst();
    double Ms = mySM.Mrun(mySM.getBKd1().getMu(),
                mySM.getQuarks(QCD::STRANGE).getMass_scale(),
                mySM.getQuarks(QCD::STRANGE).getMass(), FULLNNLO);
    double Md = mySM.Mrun(mySM.getBKd1().getMu(),
                mySM.getQuarks(QCD::DOWN).getMass_scale(),
                mySM.getQuarks(QCD::DOWN).getMass(), FULLNNLO);
    
    double X = sqrt(1.5)*FP*(MK*MK-MP*MP);
    double KK = FP/(FK-FP);
    double Q6b = -4.*sqrt(1.5)*(MK*MK/(Ms+Md))*(MK*MK/(Ms+Md))*FP/KK;
    
    me1(0) *= -1./9.*X;
    me1(1) *= 5./9.*X;
    me1(2) *= 1./3.*X;
    me1(3) = me1(2)+me1(1)-me1(0);
    me1(4) *= Q6b/3.;
    me1(5) *= Q6b;
    me1(6) *= -(1./6.*Q6b*(KK+1)-0.5*X);
    me1(7) *= -(0.5*Q6b*(KK+1)-1./6.*X);
    me1(8) = 1.5*me1(0)-0.5*me1(2);
    me1(9) = me1(1)+0.5*me1(0)-0.5*me1(2);
    me1(10) *= 0.;
    me1(11) *= 0.;
 
    switch(order) {
        case NLO:
           return((*(allcoeff[LO]) + *(allcoeff[NLO])) * me1);
        case LO:
            return((*(allcoeff[LO])) * me1);
        default:
            std::stringstream out;
            out << order;
            throw std::runtime_error("AmpDK1::AmpDK(): order " + out.str() + "not implemented");
    }
}

AmpDS1pp2()

gslpp::complex AmpDS1::AmpDS1pp2 ( orders  order )

protected

Parameters

order

Returns

the amplitude for \( K_L \) decay in 2 pion with double isospin change

Definition at line 75 of file AmpDS1.cpp.

{
    if (mySM.getFlavour().getHDS1().getCoeffDS1PP().getOrder() < order){
        std::stringstream out;
        out << order;
        throw std::runtime_error("AmpDK1::computeThValue(): requires cofficient of "
                                 "order" + out.str() + "not computed");
    }
    
    gslpp::vector<gslpp::complex> ** allcoeff = mySM.getFlavour().ComputeCoeffDS1PP(
            mySM.getBKd3().getMu(),
            mySM.getBKd3().getScheme());
    
    gslpp::vector<double> me2(mySM.getBKd3().getBpars());
    
    double MK = mySM.getMesons(QCD::K_0).getMass();
    double MP = mySM.getMesons(QCD::P_0).getMass();
    double FK = mySM.getMesons(QCD::K_0).getDecayconst();
    double FP = mySM.getMesons(QCD::P_0).getDecayconst();
    double Ms = mySM.Mrun(mySM.getBKd3().getMu(),
                mySM.getQuarks(QCD::STRANGE).getMass_scale(),
                mySM.getQuarks(QCD::STRANGE).getMass(), FULLNNLO);
    double Md = mySM.Mrun(mySM.getBKd3().getMu(),
                mySM.getQuarks(QCD::DOWN).getMass_scale(),
                mySM.getQuarks(QCD::DOWN).getMass(), FULLNNLO);
    
    double X = sqrt(1.5)*FP*(MK*MK-MP*MP);
    double KK = FP/(FK-FP);
    double Q6b = -4.*sqrt(1.5)*(MK*MK/(Ms+Md))*(MK*MK/(Ms+Md))*FP/KK;
    
    me2(0) *= 4./9.*sqrt(2.)*X;
    me2(1) *= me2(0);
    me2(2) *= 0.;
    me2(3) *= 0.;
    me2(4) *= 0.;
    me2(5) *= 0.;
    me2(6) *= -(KK/6./sqrt(2.)*Q6b*(KK+1)-X/sqrt(2.));
    me2(7) *= -(KK/2./sqrt(2.)*Q6b*(KK+1)-sqrt(2.)/6.*X);
    me2(8) = 1.5*me2(0);
    me2(9) = 1.5*me2(0);
    me2(10) = 0.;
    me2(11) = 0.;
    
    switch(order) {
        case NLO:
           return((*(allcoeff[LO]) + *(allcoeff[NLO])) * me2);
        case LO:
            return((*(allcoeff[LO])) * me2);
        default:
            std::stringstream out;
            out << order;
            throw std::runtime_error("AmpDK1::AmpDK(): order " + out.str() + "not implemented");;
    }
}

Member Data Documentation

mySM

const StandardModel& AmpDS1::mySM

private

Definition at line 42 of file AmpDS1.h.

---

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

• AmpDS1.h
• AmpDS1.cpp