78{
79
82 double mu = 0.;
83
84
85
86
87
88
89
90
92
95
96 switch (scheme) {
98 for (unsigned int i = 0; i < mcbCC.size(); i++) {
99 if (i == 0) {
100 for (
int j =
LO; j <= ordDF1; j++) {
101 for (
int k =
LO; k <= j; k++) {
102
104 u->Df1Evolnlep(mu, mcbCC[i].getMu(),
107 }
108 }
109 switch (ordqedDF1) {
112 LO,
LO_QED, mcbCC[i].getScheme()) * (*(mcbCC[i].getCoeff(
NLO))) +
113 u->Df1Evolnlep(mu, mcbCC[i].getMu(),
115 u->Df1Evolnlep(mu, mcbCC[i].getMu(),
116 LO,
NLO_QED11, mcbCC[i].getScheme()) * (*(mcbCC[i].getCoeff(
LO))) +
117 u->Df1Evolnlep(mu, mcbCC[i].getMu(),
121 LO,
LO_QED, mcbCC[i].getScheme()) * (*(mcbCC[i].getCoeff(
LO))) +
122 u->Df1Evolnlep(mu, mcbCC[i].getMu(),
124 break;
125 default:
126 throw std::runtime_error("Error in HeffDS1::ComputeCoeffDS1PPz()");
127 }
132 for (int l = 2; l < 10; l++){
137 }
139
140 if(muc != mu){
141 switch (ordDF1) {
147 break;
148 default:
149 throw std::runtime_error("Error in HeffDS1::ComputeCoeffDS1PPz()");
150 }
151 switch (ordqedDF1) {
160 break;
161 default:
162 throw std::runtime_error("Error in HeffDS1::ComputeCoeffDS1PPz()");
163 }
164 }
165 else{
166 switch (ordDF1) {
171 break;
172 default:
173 throw std::runtime_error("Error in HeffDS1::ComputeCoeffDS1PPz()");
174 }
175 switch (ordqedDF1) {
180 break;
181 default:
182 throw std::runtime_error("Error in HeffDS1::ComputeCoeffDS1PPz()");
183 }
184 }
185 }
186 }
188 default:
189 throw "HeffDS1::ComputeCoeffDS1PPz(double mu, schemes scheme): scheme not implemented";
190 }
191}
void CharmMatch()
compute the matching at the charm threshold within the SM, NDR scheme implemented
master
1.9.8