273 for (
unsigned int i = 0; i < mcbdm.size(); i++) {
275 for (
int j =
LO; j <= ordDF1; j++) {
276 for (
int k =
LO; k <= j; k++) {
277 for (
int l = 0; l < 8; l++) {
293 double eta = alsM / alsmu;
297 double B00E = 80./9., B01E = 176./9.;
299 double logeta =
log(eta);
300 double fatt = (B00E * B10S /B00S /B00S - B01E /B00S);
301 double app = B00E * (1. - eta)/ B00S;
310 for (
unsigned int i = 0; i < mcbdm.size(); i++){
312 for (j =
LO; j <= ordDF1; j++){
313 for (
int k =
LO; k <= j; k++){
322 for (j =
LO_ew; j <= ordDF1ew; j++){
334 (*(mcbdm[i].getCoeff(
orders(
LO)))) +
340 (*(mcbdm[i].getCoeff(
orders(
NLO))))) +
pow((logeta) * fatt,2) *
342 (*(mcbdm[i].getCoeff(
orders(
LO)))) ) +
pow( app, 2 ) *
346 (*(mcbdm[i].getCoeff(
orders(
LO)))) +
348 (*(mcbdm[i].getCoeff(
orders(
NLO)))) )+ logeta * fatt
365 (*(mcbdm[i].getCoeff(
orders(
LO))))) + ((logeta/eta) * fatt)
368 (*(mcbdm[i].getCoeff(
orders(
LO))))) + ( app * app/( eta) ) *
386 (*(mcbdm[i].getCoeff(
orders(
LO)))) +
389 eta * logeta * fatt *
398 (*(mcbdm[i].getCoeff(
orders(
LO)))) +
408 (*(mcbdm[i].getCoeff(
orders(
LO)))) + ( app
409 * app/(eta * eta ) ) *
419 (*(mcbdm[i].getCoeff(
orders(
LO)))) +
423 (*(mcbdm[i].getCoeff(
orders(
NLO)))) - logeta * fatt *
425 (*(mcbdm[i].getCoeff(
orders(
LO))))) - app *
438 (*(mcbdm[i].getCoeff(
orders(
LO)))) - ( app /( eta ) ) *
virtual void setMu(double mu)
const StandardModel & model
double Beta1(const double nf) const
The coefficient for a certain number of flavours .
virtual StandardModelMatching * getMyMatching() const
A get method to access the member pointer of type StandardModelMatching.
std::vector< WilsonCoefficient > Bdmumu_WC_cache
orders
An enum type for orders in QCD.
complex pow(const complex &z1, const complex &z2)
gslpp::vector< gslpp::complex > ** getCoeff() const
void setCoeff(const gslpp::vector< gslpp::complex > &z, orders order_i)
orders_ew getOrder_ew() const
orders_ew
An enum type for orders in electroweak.
void setScheme(schemes scheme)
std::vector< double > Bdmumu_Mu_cache
schemes Bdmumu_scheme_cache
complex log(const complex &z)
gslpp::matrix< double > & Df1Evol(double mu, double M, orders order, orders_ew order_ew, schemes scheme=NDR)
double Beta0(const double nf) const
The coefficient for a certain number of flavours .
double alphatilde_s(double mu)
WilsonCoefficient coeffdmumu