void EvolDF1nlep::Df1threshold_nlep(double M, double nf){
matrix<double> drsT(dim,0.), dreT(dim,0.);
double alsM = model.Als(M) / 4. / M_PI;
double ale = model.getAle() / 4. / M_PI ;
drsT = alsM * Df1threshold_deltarsT(nf);
dreT = ale * Df1threshold_deltareT(nf);
switch(order_ew){
case NLO_ew:
*elem[NLO_ew] += (*elem[LO])*dreT + (*elem[LO_ew]) * drsT ;
break;
default:
throw std::runtime_error("Error in EvolDF1nlep::Df1threshold_nlep()");
}
switch(order){
case NNLO:
*elem[NNLO] = 0.;
case NLO:
*elem[NLO] += (*elem[LO])* drsT ;
break;
default:
throw std::runtime_error("Error in EvolDF1nlep::Df1threshold_nlep()");
}
}
void EvolDC1Buras::DC1PenguinThresholds(double M, orders order)
{
double alsM = model.Als(M) / 4. / M_PI;
gslpp::matrix<double> drsT(dim,0.);
drsT = alsM * StrongThresholds();
*elem[NLO] = (*elem[NLO]) + (*elem[LO]) * drsT;
}
