#include <EvolDF1_diujlknu.h>

Inheritance diagram for EvolDF1_diujlknu:

Detailed Description

Definition at line 15 of file EvolDF1_diujlknu.h.

Public Member Functions
gslpp::matrix< double >	AnomalousDimension (orders order, unsigned int nf) const
	ADM in the JMS basis, in the order CnueduVLLkkij, CnueduVLRkkij, CnueduSRRkkij, CnueduSRLkkij, CnueduTRRkkij, from 1706.00410. More...

gslpp::matrix< double >	AnomalousDimension (orders_qed order_qed, unsigned int nf) const
	QED ADM in the JMS basis, in the order CnueduVLLkkij, CnueduVLRkkij, CnueduSRRkkij, CnueduSRLkkij, CnueduTRRkkij, from 1706.00410. More...

gslpp::matrix< double > &	Df1diujlknuEvol (double mu, double M, orders order, schemes scheme=NDR)

	EvolDF1_diujlknu (unsigned int dim_i, schemes scheme, orders order, const StandardModel &model)
	constructor More...

virtual	~EvolDF1_diujlknu ()

Public Member Functions inherited from RGEvolutor
gslpp::matrix< double > *	Evol (orders order)
	Evolution matrix set at a fixed order of QCD coupling. More...

gslpp::matrix< double > *	Evol (orders_qed order_qed)
	Evolution matrix set at a fixed order of Electroweak coupling. More...

gslpp::matrix< double > **	getEvol () const

double	getM () const
	Retrieve the upper scale of the Wilson Coefficients. More...

	RGEvolutor (unsigned int dim, schemes scheme, orders order)
	constructor More...

	RGEvolutor (unsigned int dim, schemes scheme, orders order, orders_qed order_qed)
	constructor More...

void	setEvol (const gslpp::matrix< double > &m, orders order_i)

void	setEvol (const gslpp::matrix< double > &m, orders_qed order_qed_i)

void	setEvol (unsigned int i, unsigned int j, double x, orders order_i)

void	setEvol (unsigned int i, unsigned int j, double x, orders order_i, orders_qed order_qed)

void	setM (double M)
	Sets the upper scale for the running of the Wilson Coefficients. More...

void	setMu (double mu)
	Sets the lower scale for the running of the Wilson Coefficients. More...

void	setScales (double mu, double M)
	Sets the upper and lower scale for the running of the Wilson Coefficients. More...

virtual	~RGEvolutor ()
	destructor More...

Public Member Functions inherited from WilsonTemplate< gslpp::matrix< double > >
double	getMu () const

orders	getOrder () const

orders_qed	getOrder_qed () const

schemes	getScheme () const

unsigned int	getSize () const

virtual void	resetCoefficient ()

void	setScheme (schemes scheme)

	WilsonTemplate (const WilsonTemplate< gslpp::matrix< double > > &orig)

	WilsonTemplate (unsigned int dim, schemes scheme_i, orders order_i, orders_qed order_qed_i=NO_QED)

virtual	~WilsonTemplate ()

Private Member Functions
void	Df1diujlknuEvol (double mu, double M, double nf, schemes scheme)

Private Attributes
double	a [5]

double	alsMZ_cache

double	b [5][5][5]

double	c [3][5][5][5]

double	d [3][5][5][5]

unsigned int	dim

const StandardModel &	model

double	Mz_cache

Additional Inherited Members
Protected Member Functions inherited from WilsonTemplate< gslpp::matrix< double > >
gslpp::matrix< double > *	Elem (orders order) const

gslpp::matrix< double > *	Elem (orders_qed order_qed) const

void	setElem (const gslpp::matrix< double > &v, orders order_i)

void	setElem (const gslpp::matrix< double > &v, orders_qed order_qed_i)

Protected Attributes inherited from RGEvolutor
double	M

Protected Attributes inherited from WilsonTemplate< gslpp::matrix< double > >
gslpp::matrix< double > *	elem [MAXORDER_QED+1]

double	mu

orders	order

orders_qed	order_qed

schemes	scheme

unsigned int	size

Constructor & Destructor Documentation

◆ EvolDF1_diujlknu()

EvolDF1_diujlknu::EvolDF1_diujlknu	(	unsigned int	dim_i,
		schemes	scheme,
		orders	order,
		const StandardModel &	model
	)

constructor

Parameters

dim
scheme
order
model

Definition at line 13 of file EvolDF1_diujlknu.cpp.

                                                                                                               :   
RGEvolutor(dim_i, scheme, order),
    model(model),
    dim(dim_i){
 
    gslpp::matrix<double> g0t(AnomalousDimension(LO, 3).transpose());
 
    gslpp::matrix<gslpp::complex>vv(dim, dim, 0.);
    gslpp::vector<gslpp::complex>ee(dim, 0.);
 
    g0t.eigensystem(vv, ee);
 
    gslpp::matrix<double> v(vv.real());
    gslpp::vector<double> e(ee.real());
 
    gslpp::matrix<double> vi = v.inverse();
    for (unsigned int k = 0; k < dim; k++) {
        a[k] = e(k);
//        std::cout << "a[" << k << "] = " << a[k] << std::endl;
        for (unsigned int j = 0; j < dim; j++) {
            for (unsigned int i = 0; i < dim; i++) {
                b[i][j][k] = v(i, k) * vi(k, j);
//                if (b[i][j][k] != 0.)
//                    std::cout << "b[" << i << "][" << j << "][" << k << "] = " << b[i][j][k] << std::endl;
            }
        }
    }
 
    gslpp::matrix<double> h(dim, dim, 0.);
    for (int l = 0; l < 2; l++) {
        gslpp::matrix<double> gg = vi * (AnomalousDimension(NLO, 5 - l).transpose()) * v;
        double b0 = model.Beta0(5 - l);
        for (unsigned int i = 0; i < dim; i++)
            for (unsigned int j = 0; j < dim; j++)
                h(i, j) = (i == j) * e(i) * model.Beta1(5 - l) / (2. * b0 * b0) -
                gg(i, j) / (2. * b0 + e(i) - e(j));
        gslpp::matrix<double> j = v * h * vi;
        gslpp::matrix<double> jv = j*v;
        gslpp::matrix<double> vij = vi*j;
        for (unsigned int i = 0; i < dim; i++)
            for (unsigned int j = 0; j < dim; j++)
                for (unsigned int k = 0; k < dim; k++) {
                    c[l][i][j][k] = jv(i, k) * vi(k, j);
//                    if (c[l][i][j][k] != 0.)
//                        std::cout << "c[" << l << "][" << i << "][" << j << "][" << k << "] = " << c[l][i][j][k] << std::endl;
                    d[l][i][j][k] = -v(i, k) * vij(k, j);
//                    if (d[l][i][j][k] != 0.)
//                        std::cout << "d[" << l << "][" << i << "][" << j << "][" << k << "] = " << d[l][i][j][k] << std::endl;
                }
    }
 
}

◆ ~EvolDF1_diujlknu()

EvolDF1_diujlknu::~EvolDF1_diujlknu ( )

virtual

Definition at line 66 of file EvolDF1_diujlknu.cpp.

66 {

67}

Member Function Documentation

◆ AnomalousDimension() [1/2]

gslpp::matrix< double > EvolDF1_diujlknu::AnomalousDimension	(	orders	order,
		unsigned int	nf
	)		const

ADM in the JMS basis, in the order CnueduVLLkkij, CnueduVLRkkij, CnueduSRRkkij, CnueduSRLkkij, CnueduTRRkkij, from 1706.00410.

Parameters

order
nf	number of active flavours

Returns: Anomalous dimension for \( \bar{d}_i u_j \bar{\nu} \ell_k \) operators in the JMS basis ordered as CnueduVLLkkij, CnueduVLRkkij, CnueduSRRkkij, CnueduSRLkkij, CnueduTRRkkij

Definition at line 69 of file EvolDF1_diujlknu.cpp.

{
    gslpp::matrix<double> ad(dim, dim, 0.);
            switch (order) {
                case LO:
                    ad(2, 2) = -4.;
                    ad(3, 3) = -4.;
                    ad(4, 4) = 4./3.;
                    break;
                case NLO:
                    ad(2, 2) = 2./9.*(-303.+10.*nf);
                    ad(3, 3) = ad(2,2);
                    ad(4, 4) = 2./27.*(543.-26.*nf);
                    break;
                default:
                    throw std::runtime_error("EvolDF1_diujlknu::AnomalousDimension(): order not implemented");
            }
    return (ad);
}

◆ AnomalousDimension() [2/2]

gslpp::matrix< double > EvolDF1_diujlknu::AnomalousDimension	(	orders_qed	order_qed,
		unsigned int	nf
	)		const

QED ADM in the JMS basis, in the order CnueduVLLkkij, CnueduVLRkkij, CnueduSRRkkij, CnueduSRLkkij, CnueduTRRkkij, from 1706.00410.

Parameters

order
nf	number of active flavours

Returns: QED Anomalous dimension for \( \bar{d}_i u_j \bar{\nu} \ell_k \) operators in the JMS basis ordered as CnueduVLLkkij, CnueduVLRkkij, CnueduSRRkkij, CnueduSRLkkij, CnueduTRRkkij

Definition at line 89 of file EvolDF1_diujlknu.cpp.

{
    gslpp::matrix<double> ad(dim, dim, 0.);
            switch (order_qed) {
                case LO_QED:
                    ad(0, 0) = -2.;
                    ad(1, 1) = -1.;
                    ad(2, 2) = 2./3.;
                    ad(3, 3) = 2./3.;
                    ad(3, 4) = 1./12.;
                    ad(4, 3) = 4.;
                    ad(4, 4) = -20./9.;
                    break;
                default:
                    throw std::runtime_error("EvolDF1_diujlknu::AnomalousDimension(): order_qed not implemented");
            }
    return (ad);
}

◆ Df1diujlknuEvol() [1/2]

void EvolDF1_diujlknu::Df1diujlknuEvol	(	double	mu,
		double	M,
		double	nf,
		schemes	scheme
	)

private

Definition at line 157 of file EvolDF1_diujlknu.cpp.

{
 
    gslpp::matrix<double> resLO(dim, 0.), resNLO(dim, 0.), resNNLO(dim, 0.);
 
    int l = 5 - (int) nf;
    double alsM = model.Als(M, FULLNLO) / 4. / M_PI;
    double alsmu = model.Als(mu, FULLNLO) / 4. / M_PI;
 
    double eta = alsM / alsmu;
 
    for (unsigned int k = 0; k < dim; k++) {
        double etap = pow(eta, a[k] / 2. / model.Beta0(nf));
        for (unsigned int i = 0; i < dim; i++)
            for (unsigned int j = 0; j < dim; j++) {
                resNNLO(i, j) += 0.;
                resNLO(i, j) += c[l][i][j][k] * etap * alsmu;
                resNLO(i, j) += d[l][i][j][k] * etap * alsM;
                resLO(i, j) += b[i][j][k] * etap;
            }
    }
    switch (order) {
        case NNLO:
            *elem[NNLO] = 0.; // Marco can implement it if he wishes to!
        case NLO:
            *elem[NLO] = (*elem[LO]) * resNLO + (*elem[NLO]) * resLO;
        case LO:
            *elem[LO] = (*elem[LO]) * resLO;
            break;
        case FULLNNLO:
        case FULLNLO:
        default:
            throw std::runtime_error("Error in EvolDF1_diujlknu::Df1diujlknuEvol()");
    }
}

◆ Df1diujlknuEvol() [2/2]

gslpp::matrix< double > & EvolDF1_diujlknu::Df1diujlknuEvol	(	double	mu,
		double	M,
		orders	order,
		schemes	scheme = `NDR`
	)

Parameters

mu	low energy scale
M	matching scale
order
scheme

Returns: the Wilson coefficients evolved from the scale M to the scale mu

Definition at line 108 of file EvolDF1_diujlknu.cpp.

{
    switch (scheme) {
        case NDR:
            break;
        case LRI:
        case HV:
        default:
            std::stringstream out;
            out << scheme;
            throw std::runtime_error("EvolDF1_diujlknu::Df1diujlknuEvol(): scheme " + out.str() + " not implemented ");
    }
 
    double alsMZ = model.getAlsMz();
    double Mz = model.getMz();
    if(alsMZ == alsMZ_cache && Mz == Mz_cache) {
        if (mu == this->mu && M == this->M && scheme == this->scheme)
            return (*Evol(order));        
    }
    alsMZ_cache = alsMZ;
    Mz_cache = Mz;
 
    if (M < mu) {
        std::stringstream out;
        out << "M = " << M << " < mu = " << mu;
        throw out.str();
    }
 
    setScales(mu, M); // also assign evol to identity
    if (M != mu) {
        double m_down = mu;
        double m_up = model.AboveTh(m_down);
        double nf = model.Nf(m_down);
 
        while (m_up < M) {
            Df1diujlknuEvol(m_down, m_up, nf, scheme);
            m_down = m_up;
            m_up = model.AboveTh(m_down);
            nf += 1.;
        }
        Df1diujlknuEvol(m_down, M, nf, scheme);
    }
    
    //add the leading QED log without resumming it for the moment, to be improved later
    *elem[LO] = (*elem[LO]) + AnomalousDimension(LO_QED,3).transpose()*model.getAle()/2./M_PI*log(mu/M);
    
    return (*Evol(order));
}

Member Data Documentation

◆ a

double EvolDF1_diujlknu::a[5]

private

Definition at line 61 of file EvolDF1_diujlknu.h.

◆ alsMZ_cache

double EvolDF1_diujlknu::alsMZ_cache

private

Definition at line 67 of file EvolDF1_diujlknu.h.

◆ b

double EvolDF1_diujlknu::b[5][5][5]

private

Definition at line 62 of file EvolDF1_diujlknu.h.

◆ c

double EvolDF1_diujlknu::c[3][5][5][5]

private

Definition at line 63 of file EvolDF1_diujlknu.h.

◆ d

double EvolDF1_diujlknu::d[3][5][5][5]

private

Definition at line 64 of file EvolDF1_diujlknu.h.

◆ dim

unsigned int EvolDF1_diujlknu::dim

private

Definition at line 66 of file EvolDF1_diujlknu.h.

◆ model

const StandardModel& EvolDF1_diujlknu::model

private

Definition at line 65 of file EvolDF1_diujlknu.h.

◆ Mz_cache

double EvolDF1_diujlknu::Mz_cache

private

Definition at line 68 of file EvolDF1_diujlknu.h.

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

Detailed Description

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Constructor & Destructor Documentation

◆ EvolDF1_diujlknu()

◆ ~EvolDF1_diujlknu()

Member Function Documentation

◆ AnomalousDimension() [1/2]

◆ AnomalousDimension() [2/2]

◆ Df1diujlknuEvol() [1/2]

◆ Df1diujlknuEvol() [2/2]

Member Data Documentation

◆ a

◆ alsMZ_cache

◆ b

◆ c

◆ d

◆ dim

◆ model

◆ Mz_cache