24{
25 double mu =
coeff.
getMu(), eta, M, alsM, kM, b0, b1, b0e, b1e;
26
27
28
29
30
31 if (mu == -1) throw std::runtime_error("Error in HeffDF1::LowScaleCoeff(): coeff not initialized.");
32 if (
model.
Nf(mu) != 5)
throw std::runtime_error(
"Error in HeffDF1::LowScaleCoeff(): defined for 5 flavours only.");
34 throw std::runtime_error("Error in HeffDF1::LowScaleCoeff(): order not computed at the high scale.");
35
39 alsM /= 4. * M_PI;
40 if (order_qed !=
QED0)
41 {
47 }
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64 switch (order_qed)
65 {
67 switch (order_qcd)
68 {
75 default:
76 throw std::runtime_error("Error in HeffDF1::LowScaleCoeff(): undefined order.");
77
78 }
80 switch (order_qcd)
81 {
88 default:
89 throw std::runtime_error("Error in HeffDF1::LowScaleCoeff(): undefined order.");
90 }
92 switch (order_qcd)
93 {
102 default:
103 throw std::runtime_error("Error in HeffDF1::LowScaleCoeff(): undefined order.");
104 }
105 default:
106 throw std::runtime_error("Error in HeffDF1::LowScaleCoeff(): undefined order.");
107 }
108}
const double Nf(const double mu) const
The number of active flavour at scale .
const double getMuw() const
A get method to retrieve the matching scale around the weak scale.
const double Ale(double mu, orders order, bool Nf_thr=true) const
The running electromagnetic coupling in the scheme.
const double Beta_s(int nm, unsigned int nf) const
QCD beta function coefficients including QED corrections - eq. (36) hep-ph/0512066.
const double Als(const double mu, const orders order, const bool Nf_thr, const bool qed_flag) const
The running QCD coupling in the scheme including QED corrections.
const double Beta_e(int nm, unsigned int nf) const
QED beta function coefficients - eq. (36) hep-ph/0512066.
qcd_orders getOrder_QCD() const
qed_orders getOrder_QED() const
master
Inheritance diagram for HeffDF1:
1.9.8 