A class for the \(\Delta F = 2\) effective Hamiltonian. More...

#include <HeffDF2.h>

Detailed Description

A class for the \(\Delta F = 2\) effective Hamiltonian.

Author: HEPfit Collaboration

Copyright: GNU General Public License

This class is used to compute the contributions to the \(\Delta F = 2\) by taking the values of the Wilson coefficients from the model matching classes and runs them down to the relevant scale for the observables.

Definition at line 28 of file HeffDF2.h.

Public Member Functions
gslpp::matrix< double >	AnomalousDimension (orders order, unsigned int nf=0) const

void	ChangeScheme (schemes schout, WilsonCoefficient &c_in, orders order)
	change scheme for a Wilson Coefficient More...

gslpp::vector< gslpp::complex > **	ComputeCoeffBd (double mu, schemes scheme=NDR, bool SM=false)

gslpp::vector< gslpp::complex > **	ComputeCoeffBs (double mu, schemes scheme=NDR, bool SM=false)

gslpp::vector< gslpp::complex > **	ComputeCoeffdd (double mu, schemes scheme=NDR)

gslpp::vector< gslpp::complex > **	ComputeCoeffK (double mu, schemes scheme=NDR)

gslpp::vector< gslpp::complex > **	ComputeCoeffmK (double mu, schemes scheme=NDR)
	for Delta M_K the SM contribution is set to zero More...

WilsonCoefficient	getCoeffBd () const

WilsonCoefficient	getCoeffBs () const

WilsonCoefficient	getCoeffDD () const

WilsonCoefficient	getCoeffK () const

WilsonCoefficient	getCoeffmK () const

EvolDF2 &	getUDF2 () const

	HeffDF2 (const StandardModel &SM)
	constructor More...

virtual	~HeffDF2 ()
	destructor More...

Private Member Functions
gslpp::complex	S0tt (double mu) const

Private Attributes
WilsonCoefficient	coeffbd

WilsonCoefficient	coeffbs

WilsonCoefficient	coeffDd

WilsonCoefficient	coeffk

WilsonCoefficient	coeffmk

gslpp::matrix< double >	drNDRLRI

std::unique_ptr< EvolDF2 >	evolDF2

const StandardModel &	model

Constructor & Destructor Documentation

◆ HeffDF2()

HeffDF2::HeffDF2 ( const StandardModel & SM )

constructor

Parameters

SM
SM_Matching

Definition at line 12 of file HeffDF2.cpp.

:       model(SM),
        drNDRLRI(5, 5, 0),
        coeffbd(5, NDR, NLO),
        coeffbs(5, NDR, NLO),
        coeffDd(5, NDR, NLO),
        coeffk(5, NDR, NLO),
        coeffmk(5, NDR, NLO),
        evolDF2(new EvolDF2(5, NDR, NLO, SM))
{
    
    double Nc = SM.getNc();
    drNDRLRI(0,0) = -(((-1. + Nc) * (-7. + log(4096.))) / Nc);
    drNDRLRI(1,1) = (-2. * (-1. + 6. * Nc * Nc - 8. * log(2.) + Nc * (-13. + log(1024.)))) / (3. * Nc);
    drNDRLRI(1,2) = (-2. * (13. - 10. * log(2.) + Nc * (-5. + log(256.)))) / (3. * Nc);
    drNDRLRI(2,1) = (-8. + 6. * Nc * Nc + 20. * log(2.) - 8. * Nc * (1. + log(4.))) / (3. * Nc);
    drNDRLRI(2,2) = (2. * (4. + Nc - 10. * Nc * log(2.) + log(256.))) / (3. * Nc);
    drNDRLRI(3,3) = (2. - 4. * Nc * Nc + log(4.)) / Nc;
    drNDRLRI(3,4) = 2. - log(4.);
    drNDRLRI(4,3) = -2. * (1. + log(2.));
    drNDRLRI(4,4) = (2. + log(4.)) / Nc;
}

◆ ~HeffDF2()

HeffDF2::~HeffDF2 ( )

virtual

destructor

Definition at line 35 of file HeffDF2.cpp.

36{}

Member Function Documentation

◆ AnomalousDimension()

gslpp::matrix< double > HeffDF2::AnomalousDimension	(	orders	order,
		unsigned int	nf = `0`
	)		const

Parameters

order
nf	is the number of active flavours

Returns: the anomalous dimension for DeltaF=2 processes

◆ ChangeScheme()

void HeffDF2::ChangeScheme	(	schemes	schout,
		WilsonCoefficient &	c_in,
		orders	order
	)

change scheme for a Wilson Coefficient

Parameters

schout	is the renormalization scheme in output
c_in	is the Wilson Coefficient to be converted to scheme schout
order

Definition at line 205 of file HeffDF2.cpp.

{
    schemes schin = c_in.getScheme();
    if (schout == schin || order == LO) return;
    WilsonCoefficient c_out(5, schout, order);
    switch(schin) {
        case NDR:
            switch(schout) {
                case LRI:
                        c_out.setCoeff(*c_in.getCoeff(NLO) -
                            model.Als(c_in.getMu()) / 4. / M_PI * drNDRLRI.transpose()*
                            (*c_in.getCoeff(LO)), NLO);
                    c_in.setCoeff(*c_out.getCoeff(NLO),NLO);
                    c_in.setScheme(schout);
                    break;
                 default:
                    throw std::runtime_error("HeffDF2::ChangeScheme(): out scheme not implemented"); 
            }
            break;
        default:
            throw std::runtime_error("HeffDF2::ChangeScheme(): in scheme not implemented"); 
    }
}

◆ ComputeCoeffBd()

gslpp::vector< gslpp::complex > ** HeffDF2::ComputeCoeffBd	(	double	mu,
		schemes	scheme = `NDR`,
		bool	SM = `false`
	)

Parameters

mu	is the low energy scale
scheme	indicates the renormalization scheme (Default: NDR)

Returns: the effective hamiltonian at the scale mu for B_d oscillations

Definition at line 38 of file HeffDF2.cpp.

{
    std::vector<WilsonCoefficient> mc;
 
    if (SM)
        mc = model.getMatching().StandardModelMatching::CMdbd2();
    else
    mc = model.getMatching().CMdbd2();
    
    coeffbd.setMu(mu);
    
    coeffbd.setScheme(mc[0].getScheme());
    
    orders ordDF2 = coeffbd.getOrder();
    for (unsigned int i = 0; i < mc.size(); i++){
        ChangeScheme(mc[0].getScheme(),mc[i],ordDF2);
        for (int j = LO; j <= ordDF2; j++){
            for (int k = LO; k <= j; k++){                
                coeffbd.setCoeff(*coeffbd.getCoeff(orders(j)) +
                    evolDF2->Df2Evol(mu, mc[i].getMu(), orders(k), mc[i].getScheme()) *
                    (*(mc[i].getCoeff(orders(j - k)))), orders(j));
            }
        }
    }
    
    ChangeScheme(scheme, coeffbd, ordDF2);
 
    return coeffbd.getCoeff();
}

◆ ComputeCoeffBs()

gslpp::vector< gslpp::complex > ** HeffDF2::ComputeCoeffBs	(	double	mu,
		schemes	scheme = `NDR`,
		bool	SM = `false`
	)

Parameters

mu	is the low energy scale
scheme	indicates the renormalization scheme (Default: NDR)
SM	boolean to select only the SM contributions to the coefficients (Default: false)

Returns: the effective hamiltonian at the scale mu for B_s oscillations

Definition at line 68 of file HeffDF2.cpp.

{
 
    std::vector<WilsonCoefficient> mc;
    
    if (SM)
        mc = model.getMatching().StandardModelMatching::CMdbs2();
    else
        mc = model.getMatching().CMdbs2();
    
    coeffbs.setMu(mu);
    
    coeffbs.setScheme(mc[0].getScheme());
    
    orders ordDF2 = coeffbs.getOrder();
    for (unsigned int i = 0; i < mc.size(); i++){
        ChangeScheme(mc[0].getScheme(),mc[i],ordDF2);
        for (int j = LO; j <= ordDF2; j++){
            for (int k = LO; k <= j; k++){
                coeffbs.setCoeff(*coeffbs.getCoeff(orders(j)) +
                    evolDF2->Df2Evol(mu, mc[i].getMu(), orders(k), mc[i].getScheme()) *
                    (*(mc[i].getCoeff(orders(j - k)))), orders(j));
            }
        }
    }
 
    ChangeScheme(scheme, coeffbs, ordDF2);
 
    return coeffbs.getCoeff();
}

◆ ComputeCoeffdd()

gslpp::vector< gslpp::complex > ** HeffDF2::ComputeCoeffdd	(	double	mu,
		schemes	scheme = `NDR`
	)

Parameters

mu	is the low energy scale
scheme	indicates the renormalization scheme (Default: NDR)

Returns: the effective hamiltonian at the scale mu for D oscillations

Definition at line 99 of file HeffDF2.cpp.

{
 
    std::vector<WilsonCoefficient> mc = model.getMatching().CMdd2();
    
    coeffDd.setMu(mu);
 
    coeffDd.setScheme(mc[0].getScheme());
 
    orders ordDF2 = coeffDd.getOrder();
    for (unsigned int i = 0; i < mc.size(); i++){
       ChangeScheme(mc[0].getScheme(),mc[i],ordDF2);
       for (int j = LO; j <= ordDF2; j++){
            for (int k = LO; k <= j; k++){
                coeffDd.setCoeff(*coeffDd.getCoeff(orders(j)) +
                    evolDF2->Df2Evol(mu, mc[i].getMu(), orders(k), mc[i].getScheme()) *
                    (*(mc[i].getCoeff(orders(j - k)))), orders(j));
            }
        }
    }
 
    ChangeScheme(scheme, coeffDd, ordDF2);
    
    return coeffDd.getCoeff();
}

◆ ComputeCoeffK()

gslpp::vector< gslpp::complex > ** HeffDF2::ComputeCoeffK	(	double	mu,
		schemes	scheme = `NDR`
	)

Parameters

mu	is the low energy scale
scheme	indicates the renormalization scheme (Default: NDR)

Returns: the effective hamiltonian at the scale mu for K oscillations

Definition at line 125 of file HeffDF2.cpp.

{
 
    std::vector<WilsonCoefficient> mc = model.getMatching().CMdk2();
    gslpp::vector<gslpp::complex> zero(5,0.);
    
    coeffk.setScheme(mc[0].getScheme());
 
    coeffk.setMu(mu);
    coeffk.setCoeff(zero,LO);
    coeffk.setCoeff(zero,NLO);
 
    orders ordDF2 = coeffk.getOrder();
    for (unsigned int i = 0; i < mc.size(); i++){
        if (i == 0){
            coeffk.setCoeff(0, evolDF2->etatt(mu) * model.getMatching().S0tt() * (1. - model.getOptionalParameter("DeltattEpsK")) //EW correction from 2108.00017
                             + evolDF2->etacc(mu) * model.getMatching().S0c()
                             + evolDF2->etact(mu) * model.getMatching().S0ct(),
                            NLO);
#if SUSYFIT_DEBUG & 2
    std::cout << "mu = " << mu<< ", S0tt = " << model.getMatching().S0tt() << 
            ", S0cc = " << model.getMatching().S0c() << 
            ", S0ct = " << model.getMatching().S0ct() << std::endl;
    std::cout << ", etatt = " << evolDF2->etatt(mu) << 
            ", etacc = " << evolDF2->etacc(mu) << 
            ", etact = " << evolDF2->etact(mu) << std::endl;
    std::cout <<  "tt = " << evolDF2->etatt(mu)*model.getMatching().S0tt() << 
            ", cc = " << evolDF2->etacc(mu)*model.getMatching().S0c() << 
            ", ct = " << evolDF2->etact(mu)*model.getMatching().S0ct() << std::endl;
#endif
            
        }
        else {
            ChangeScheme(mc[0].getScheme(),mc[i],ordDF2);
            for (int j = LO; j <= ordDF2; j++){
                for (int k = LO; k <= j; k++){
                    coeffk.setCoeff(*coeffk.getCoeff(orders(j)) +
                        evolDF2->Df2Evol(mu, mc[i].getMu(), orders(k), mc[i].getScheme()) *
                        (*(mc[i].getCoeff(orders(j - k)))), orders(j));
                }
            }
        }
    }
    
    ChangeScheme(scheme, coeffk, ordDF2);
    
    return coeffk.getCoeff();
}

◆ ComputeCoeffmK()

gslpp::vector< gslpp::complex > ** HeffDF2::ComputeCoeffmK	(	double	mu,
		schemes	scheme = `NDR`
	)

for Delta M_K the SM contribution is set to zero

Parameters

mu	is the low energy scale
scheme	indicates the renormalization scheme (Default: NDR)

Returns: the effective hamiltonian at the scale mu for Delta M_K

Definition at line 175 of file HeffDF2.cpp.

{
 
    const std::vector<WilsonCoefficient>& mc = model.getMatching().CMdk2();
    gslpp::vector<gslpp::complex> zero(5,0.);
    
    coeffmk.setMu(mu);
 
    orders ordDF2 = coeffmk.getOrder();
    for (unsigned int i = 0; i < mc.size(); i++){
        if (i == 0){
            coeffmk.setCoeff(zero, NLO);
            coeffmk.setCoeff(zero, LO);
        }
        else {
            for (int j = LO; j <= ordDF2; j++){
                for (int k = LO; k <= j; k++){
                    coeffmk.setCoeff(*coeffmk.getCoeff(orders(j)) +
                        evolDF2->Df2Evol(mu, mc[i].getMu(), orders(k), mc[i].getScheme()) *
                        (*(mc[i].getCoeff(orders(j - k)))), orders(j));
                }
            }
        }
    }
    
    ChangeScheme(scheme, coeffmk, ordDF2);
    
    return coeffmk.getCoeff();
}

◆ getCoeffBd()

WilsonCoefficient HeffDF2::getCoeffBd ( ) const

inline

Definition at line 100 of file HeffDF2.h.

                                         {
        return coeffbd;
    }

◆ getCoeffBs()

WilsonCoefficient HeffDF2::getCoeffBs ( ) const

inline

Definition at line 104 of file HeffDF2.h.

                                         {
        return coeffbs;
    }

◆ getCoeffDD()

WilsonCoefficient HeffDF2::getCoeffDD ( ) const

inline

Definition at line 108 of file HeffDF2.h.

                                         {
        return coeffDd;
    }

◆ getCoeffK()

WilsonCoefficient HeffDF2::getCoeffK ( ) const

inline

Definition at line 112 of file HeffDF2.h.

                                        {
        return coeffk;
    }

◆ getCoeffmK()

WilsonCoefficient HeffDF2::getCoeffmK ( ) const

inline

Definition at line 116 of file HeffDF2.h.

                                          {
        return coeffmk;
    }

◆ getUDF2()

EvolDF2 & HeffDF2::getUDF2 ( ) const

inline

Definition at line 120 of file HeffDF2.h.

                             {
        return *evolDF2;
    }

◆ S0tt()

gslpp::complex HeffDF2::S0tt ( double mu ) const

private

Member Data Documentation

◆ coeffbd

WilsonCoefficient HeffDF2::coeffbd

private

Definition at line 129 of file HeffDF2.h.

◆ coeffbs

WilsonCoefficient HeffDF2::coeffbs

private

Definition at line 130 of file HeffDF2.h.

◆ coeffDd

WilsonCoefficient HeffDF2::coeffDd

private

Definition at line 131 of file HeffDF2.h.

◆ coeffk

WilsonCoefficient HeffDF2::coeffk

private

Definition at line 132 of file HeffDF2.h.

◆ coeffmk

WilsonCoefficient HeffDF2::coeffmk

private

Definition at line 133 of file HeffDF2.h.

◆ drNDRLRI

gslpp::matrix<double> HeffDF2::drNDRLRI

private

Definition at line 128 of file HeffDF2.h.

◆ evolDF2

std::unique_ptr<EvolDF2> HeffDF2::evolDF2

private

Definition at line 135 of file HeffDF2.h.

◆ model

const StandardModel& HeffDF2::model

private

Definition at line 127 of file HeffDF2.h.

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

Detailed Description

Public Member Functions

Private Member Functions

Private Attributes

Constructor & Destructor Documentation

◆ HeffDF2()

◆ ~HeffDF2()

Member Function Documentation

◆ AnomalousDimension()

◆ ChangeScheme()

◆ ComputeCoeffBd()

◆ ComputeCoeffBs()

◆ ComputeCoeffdd()

◆ ComputeCoeffK()

◆ ComputeCoeffmK()

◆ getCoeffBd()

◆ getCoeffBs()

◆ getCoeffDD()

◆ getCoeffK()

◆ getCoeffmK()

◆ getUDF2()

◆ S0tt()

Member Data Documentation

◆ coeffbd

◆ coeffbs

◆ coeffDd

◆ coeffk

◆ coeffmk

◆ drNDRLRI

◆ evolDF2

◆ model