Hyperiso 1.0.3
Modular flavour-physics calculations, Wilson coefficients and statistical inference
Loading...
Searching...
No Matches
BWilson.cpp
Go to the documentation of this file.
1#include "BWilson.h"
2
3// ---------- C1 ----------
4
5C1::C1() : WilsonCoefficient("C1", GroupMapper::str(WGroup::B, ScaleType::MATCHING)) {
6
7 matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));
8 matching_info[QCDOrder::NLO] = {
9 {{"WPARAM_MATCH_SM", 3}}, // L
10 [](const auto& src) {
11 auto L = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 3);
12 return 15. + 6. * L;
13 },
14 get_lhaid_from_name(QCDOrder::NLO)
15 };
16
17 matching_info[QCDOrder::NNLO] = {
18 {
19 {"WPARAM_MATCH_SM", 3}, // L
20 {"WPARAM_MATCH_SM", LhaID(2, 1)} // xt
21 },
22 [](const auto& src) {
23 auto L = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 3);
24 auto xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
25 return -T(xt) + 7987./72. + 17. * PI2 / 3. + 475./6. * L + 17. * L * L;
26 },
27 get_lhaid_from_name(QCDOrder::NNLO)
28 };
29}
30
31// ---------- C2 ----------
32
33C2::C2() : WilsonCoefficient("C2", GroupMapper::str(WGroup::B, ScaleType::MATCHING)) {
34 matching_info[QCDOrder::LO] = {
35 {},
36 compute_LO,
37 get_lhaid_from_name(QCDOrder::LO)
38 };
39
40 matching_info[QCDOrder::NLO] = MatchingInfo(LhaID(3040405, 4141, 1, 0));
41
42 matching_info[QCDOrder::NNLO] = {
43 {
44 {"WPARAM_MATCH_SM", 3} // L
45 },
46 compute_NNLO,
47 get_lhaid_from_name(QCDOrder::NNLO)
48 };
49}
50
51double C2::compute_LO(const ParamSrc&) {
52 return 1.;
53}
54
55double C2::compute_NNLO(const ParamSrc& src) {
56 auto L = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 3);
57 return 127. / 18. + 4. / 3. * PI2 + 46. / 3. * L + 4. * L * L;
58}
59
60
61// ---------- C3 ----------
62
63C3::C3() : WilsonCoefficient("C3", GroupMapper::str(WGroup::B, ScaleType::MATCHING)) {
64
65 matching_info[QCDOrder::NLO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));
66
67 matching_info[QCDOrder::NLO] = MatchingInfo(get_lhaid_from_name(QCDOrder::NLO));
68
69 matching_info[QCDOrder::NNLO] = {
70 {
71 {"WPARAM_MATCH_SM", 3}, // L
72 {"WPARAM_MATCH_SM", 6}, // mass_top_muW
73 {"WPARAM_MATCH_SM", LhaID(2, 1)}, // xt
74 {"EW_SCALE", 1} // Q_match
75 },
76 compute_NNLO,
77 get_lhaid_from_name(QCDOrder::NNLO)
78 };
79}
80
81double C3::compute_NNLO(const ParamSrc& src) {
82 double L = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 3);
83 double mass_top_muW = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 6);
84 double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
85 double Q_match = src.get_val(ParameterType::WILSON, "EW_SCALE", 1);
86
87 double coeff_temp = G1t(xt, log(Q_match * Q_match / (mass_top_muW * mass_top_muW)))
88 - 680. / 243.
89 - 20. / 81. * PI2
90 - 68. / 81. * L
91 - 20. / 27. * L * L;
92
93 return coeff_temp;
94}
95
96// ---------- C4 ----------
97
98C4::C4() : WilsonCoefficient("C4", GroupMapper::str(WGroup::B, ScaleType::MATCHING)) {
99
100 matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));
101
102 // NLO
103 matching_info[QCDOrder::NLO] = {
104 {
105 {"WPARAM_MATCH_SM", 3}, // L
106 {"WPARAM_MATCH_SM", LhaID(2, 1)} // x_t
107 },
108 compute_NLO,
109 get_lhaid_from_name(QCDOrder::NLO)
110 };
111
112 // NNLO
113 matching_info[QCDOrder::NNLO] = {
114 {
115 {"WPARAM_MATCH_SM", 3}, // L
116 {"WPARAM_MATCH_SM", 6}, // mass_top_muW
117 {"WPARAM_MATCH_SM", LhaID(2, 1)}, // x_t
118 {"EW_SCALE", 1} // Q_match
119 },
120 compute_NNLO,
121 get_lhaid_from_name(QCDOrder::NNLO)
122 };
123}
124
125double C4::compute_NLO(const ParamSrc& src) {
126 double L = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 3);
127 double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
128 return E0t(xt) - 7. / 9. + 2. / 3. * L;
129}
130
131double C4::compute_NNLO(const ParamSrc& src) {
132 double L = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 3);
133 double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
134 double mass_top_muW = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 6);
135 double Q_match = src.get_val(ParameterType::WILSON, "EW_SCALE", 1);
136 return E1t(xt, log(Q_match * Q_match / (mass_top_muW * mass_top_muW)))
137 + 950. / 243.
138 + 10. / 81. * PI2
139 + 124. / 27. * L
140 + 10. / 27. * L * L;
141}
142
143// ---------- C5 ----------
144
145C5::C5() : WilsonCoefficient("C5", GroupMapper::str(WGroup::B, ScaleType::MATCHING)) {
146
147 matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));
148
149 matching_info[QCDOrder::NLO] = MatchingInfo(get_lhaid_from_name(QCDOrder::NLO));
150
151 matching_info[QCDOrder::NNLO] = {
152 {
153 {"WPARAM_MATCH_SM", 3}, // L
154 {"WPARAM_MATCH_SM", 6}, // mass_top_muW
155 {"WPARAM_MATCH_SM", LhaID(2, 1)}, // xt
156 {"EW_SCALE", 1} // Q_match
157 },
158 compute_NNLO,
159 get_lhaid_from_name(QCDOrder::NNLO)
160 };
161}
162
163double C5::compute_NNLO(const ParamSrc& src) {
164 double L = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 3);
165 double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
166 double mass_top_muW = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 6);
167 double Q_match = src.get_val(ParameterType::WILSON, "EW_SCALE", 1);
168
169 double coeff_temp =
170 -G1t(xt, log(Q_match * Q_match / (mass_top_muW * mass_top_muW))) / 10.
171 + 2. / 15. * E0t(xt)
172 + 68. / 243.
173 + 2. / 81. * PI2
174 + 14. / 81. * L
175 + 2. / 27. * L * L;
176
177 return coeff_temp;
178}
179
180// ---------- C6 ----------
181
182C6::C6() : WilsonCoefficient("C6", GroupMapper::str(WGroup::B, ScaleType::MATCHING)) {
183
184 matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));
185
186 matching_info[QCDOrder::NLO] = MatchingInfo(get_lhaid_from_name(QCDOrder::NLO));
187
188 matching_info[QCDOrder::NNLO] = {
189 {
190 {"WPARAM_MATCH_SM", 3}, // L
191 {"WPARAM_MATCH_SM", 6}, // mass_top_muW
192 {"WPARAM_MATCH_SM", LhaID(2, 1)}, // xt
193 {"EW_SCALE", 1} // Q_match
194 },
195 compute_NNLO,
196 get_lhaid_from_name(QCDOrder::NNLO)
197 };
198}
199
200double C6::compute_NNLO(const ParamSrc& src) {
201 double L = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 3);
202 double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
203 double mass_top_muW = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 6);
204 double Q_match = src.get_val(ParameterType::WILSON, "EW_SCALE", 1);
205
206 double coeff_temp =
207 -3. / 16. * G1t(xt, log(Q_match * Q_match / (mass_top_muW * mass_top_muW)))
208 + E0t(xt) / 4.
209 + 85. / 162.
210 + 5. / 108. * PI2
211 + 35. / 108. * L
212 + 5. / 36. * L * L;
213
214 return coeff_temp;
215}
216
217
218// ---------- C7 ----------
219
220C7::C7() : WilsonCoefficient("C7", GroupMapper::str(WGroup::B, ScaleType::MATCHING)) {
221 // LO
222 matching_info[QCDOrder::LO] = {
223 {
224 {ParameterType::WILSON, "WPARAM_MATCH_SM", LhaID(2, 1)} // x_t
225 },
226 compute_LO,
227 get_lhaid_from_name(QCDOrder::LO)
228 };
229
230 // std::cout << "matching_info LO : " << get_lhaid_from_name(QCDOrder::LO) << std::endl;
231 // NLO
232 matching_info[QCDOrder::NLO] = {
233 {
234 {"WPARAM_MATCH_SM", 3}, // L
235 {"WPARAM_MATCH_SM", 6}, // mass_top_muW
236 {"WPARAM_MATCH_SM", LhaID(2, 1)}, // x_t
237 {"EW_SCALE", 1} // Q_match
238 },
239 compute_NLO,
240 get_lhaid_from_name(QCDOrder::NLO)
241 };
242
243 // NNLO
244 matching_info[QCDOrder::NNLO] = {
245 {
246 // {"WPARAM_MATCH_SM", 3}, // L
247 {"WPARAM_MATCH_SM", 6}, // mass_top_muW
248 {"WPARAM_MATCH_SM", LhaID(2, 1)}, // x_t
249 {"WPARAM_MATCH_SM", 7}, // xtW (top pole)
250 {"WPARAM_MATCH_SM", 8}, // xtt
251 {"EW_SCALE", 1}, // Q_match
252 {ParameterType::SM, "MASS", 24} // m_W
253 },
254 compute_NNLO,
255 get_lhaid_from_name(QCDOrder::NNLO)
256 };
257}
258
259double C7::compute_LO(const ParamSrc& src) {
260 double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
261 return -0.5 * A0t(xt) - 23. / 36.;
262}
263
264double C7::compute_NLO(const ParamSrc& src) {
265 double L = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 3);
266 double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
267 double mass_top_muW = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 6);
268 double Q_match = src.get_val(ParameterType::WILSON, "EW_SCALE", 1);
269
270 return -0.5 * A1t(xt, log(Q_match * Q_match / (mass_top_muW * mass_top_muW)))
271 + 713. / 243. + 4. / 81. * L
272 - 4. / 9. * (E0t(xt) - 7. / 9. + 2. / 3. * L);
273}
274
275double C7::compute_NNLO(const ParamSrc& src) {
276 // double L = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 3);
277 double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
278 double xtW = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 7);
279 double xtt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 8);
280 double mtop = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 6);
281 double Q = src.get_val(ParameterType::WILSON, "EW_SCALE", 1);
282 double mW = src.get_val(ParameterType::SM, "MASS", 24);
283
284 double logqt = log(Q * Q / (mtop * mtop));
285 double logqW = log(Q * Q / (mW * mW));
286
287 double coeff_temp = C7t2mt(xtt)
288 + logqt * (
289 ((-592. * pow(xt, 5.) - 22. * pow(xt, 4.) + 12814. * pow(xt, 3.) - 6376. * xt * xt + 512. * xt)
290 / 27. / pow(xt - 1., 5.)) * Li2(1. - 1. / xt)
291 + ((-26838. * pow(xt, 5.) + 25938. * pow(xt, 4.) + 627367. * pow(xt, 3.) - 331956. * xt * xt + 16989. * xt - 460.)
292 / 729. / pow(xt - 1., 6.)) * log(xt)
293 + ((34400. * pow(xt, 5.) + 276644. * pow(xt, 4.) - 2668324. * pow(xt, 3.) + 1694437. * xt * xt - 323354. * xt + 53077.)
294 / 2187. / pow(xt - 1., 5.))
295 + logqt * (
296 ((-63. * pow(xt, 5.) + 532. * pow(xt, 4.) + 2089. * pow(xt, 3.) - 1118. * xt * xt)
297 / 9. / pow(xt - 1., 6.)) * log(xt)
298 + ((1186. * pow(xt, 5.) - 2705. * pow(xt, 4.) - 24791. * pow(xt, 3.) - 16099. * xt * xt + 19229. * xt - 2740.)
299 / 162. / pow(xt - 1., 5.))
300 )
301 )
302 - (C7c2MW(xtW) + 13763. / 2187. * logqW + 814. / 729. * pow(logqW, 2.));
303
304 return coeff_temp;
305}
306
307// ---------- C8 ----------
308
309C8::C8() : WilsonCoefficient("C8", GroupMapper::str(WGroup::B, ScaleType::MATCHING)) {
310 // LO
311 matching_info[QCDOrder::LO] = {
312 {
313 {"WPARAM_MATCH_SM", LhaID(2, 1)} // x_t
314 },
315 compute_LO,
316 get_lhaid_from_name(QCDOrder::LO)
317 };
318
319 // NLO
320 matching_info[QCDOrder::NLO] = {
321 {
322 {"WPARAM_MATCH_SM", 3}, // L
323 {"WPARAM_MATCH_SM", 6}, // mass_top_muW
324 {"WPARAM_MATCH_SM", LhaID(2, 1)}, // x_t
325 {"EW_SCALE", 1} // Q_match
326 },
327 compute_NLO,
328 get_lhaid_from_name(QCDOrder::NLO)
329 };
330
331 // NNLO
332 matching_info[QCDOrder::NNLO] = {
333 {
334 // {"WPARAM_MATCH_SM", 3}, // L
335 {"WPARAM_MATCH_SM", 6}, // mass_top_muW
336 {"WPARAM_MATCH_SM", LhaID(2, 1)}, // x_t
337 {"WPARAM_MATCH_SM", 7}, // xtW
338 {"WPARAM_MATCH_SM", 8}, // xtt
339 {"EW_SCALE", 1}, // Q_match
340 {ParameterType::SM, "MASS", 24} // m_W
341 },
342 compute_NNLO,
343 get_lhaid_from_name(QCDOrder::NNLO)
344 };
345}
346
347double C8::compute_LO(const ParamSrc& src) {
348 double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
349 return -0.5 * F0t(xt) - 1. / 3.;
350}
351
352double C8::compute_NLO(const ParamSrc& src) {
353 double L = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 3);
354 double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
355 double mass_top_muW = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 6);
356 double Q_match = src.get_val(ParameterType::WILSON, "EW_SCALE", 1);
357
358 return -0.5 * F1t(xt, log(Q_match * Q_match / (mass_top_muW * mass_top_muW)))
359 + 91. / 324.
360 - 4. / 27. * L
361 - (E0t(xt) - 7. / 9. + 2. / 3. * L) / 6.;
362}
363
364double C8::compute_NNLO(const ParamSrc& src) {
365 // double L = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 3);
366 double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
367 double xtW = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 7);
368 double xtt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 8);
369 double mtop = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 6);
370 double Q = src.get_val(ParameterType::WILSON, "EW_SCALE", 1);
371 double mW = src.get_val(ParameterType::SM, "MASS", 24);
372
373 double logqt = log(Q * Q / (mtop * mtop));
374 double logqW = log(Q * Q / (mW * mW));
375
376 double coeff_temp = C8t2mt(xtt)
377 + logqt * (
378 ((-148. * pow(xt, 5.) + 1052. * pow(xt, 4.) - 4811. * pow(xt, 3.) - 3520. * xt * xt - 61. * xt)
379 / 18. / pow(xt - 1., 5.)) * Li2(1. - 1. / xt)
380 + ((-15984. * pow(xt, 5.) + 152379. * pow(xt, 4.) - 1358060. * pow(xt, 3.) - 1201653. * xt * xt - 74190. * xt + 9188.)
381 / 1944. / pow(xt - 1., 6.)) * log(xt)
382 + ((109669. * pow(xt, 5.) - 1112675. * pow(xt, 4.) + 6239377. * pow(xt, 3.) + 8967623. * xt * xt + 768722. * xt - 42796.)
383 / 11664. / pow(xt - 1., 5.))
384 + logqt * (
385 ((-139. * pow(xt, 4.) - 2938. * pow(xt, 3.) - 2683. * xt * xt)
386 / 12. / pow(xt - 1., 6.)) * log(xt)
387 + ((1295. * pow(xt, 5.) - 7009. * pow(xt, 4.) + 29495. * pow(xt, 3.) + 64513. * xt * xt + 17458. * xt - 2072.)
388 / 216. / pow(xt - 1., 5.))
389 )
390 )
391 - (C8c2MW(xtW) + 16607. / 5832. * logqW + 397. / 486. * pow(logqW, 2.));
392
393 return coeff_temp;
394}
395
396// ---------- C9 ----------
397
398
399C9::C9() : WilsonCoefficient("C9", GroupMapper::str(WGroup::B, ScaleType::MATCHING)) {
400 // LO
401 matching_info[QCDOrder::LO] = {
402 {
403 {"WPARAM_MATCH_SM", LhaID(2, 1)}, // x_t
404 {"WPARAM_MATCH_SM", 3}, // L
405 {"WPARAM_SI_SM", 4} // sw2
406 },
407 compute_LO,
408 get_lhaid_from_name(QCDOrder::LO)
409 };
410
411 // NLO
412 matching_info[QCDOrder::NLO] = {
413 {
414 {"WPARAM_MATCH_SM", 3}, // L
415 {"WPARAM_MATCH_SM", 6}, // mass_top_muW
416 {"WPARAM_MATCH_SM", LhaID(2, 1)}, // x_t
417 {"WPARAM_SI_SM", 4}, // sw2
418 {"EW_SCALE", 1} // Q_match
419 },
420 compute_NLO,
421 get_lhaid_from_name(QCDOrder::NLO)
422 };
423
424 matching_info[QCDOrder::NNLO] = MatchingInfo(get_lhaid_from_name(QCDOrder::NNLO));
425
426}
427
428double C9::compute_LO(const ParamSrc& src) {
429 double L = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 3);
430 double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
431 double sw2 = src.get_val(ParameterType::WILSON, "WPARAM_SI_SM", 4);
432
433 // printf("L = %.5f\n", L);
434 // printf("x_t = %.5f\n", xt);
435 // printf("sw2 = %.5f\n", sw2);
436
437 // printf("C9_LO = %.5f\n", (1. - 4. * sw2) / sw2 * C0t(xt)
438 // - B0t(xt) / sw2
439 // - D0t(xt)
440 // + 1. / (4. * sw2)
441 // + 38. / 27.
442 // - 4. / 9. * L);
443
444 return (1. - 4. * sw2) / sw2 * C0t(xt)
445 - B0t(xt) / sw2
446 - D0t(xt)
447 + 1. / (4. * sw2)
448 + 38. / 27.
449 - 4. / 9. * L;
450}
451
452double C9::compute_NLO(const ParamSrc& src) {
453 double L = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 3);
454 double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
455 double mtop = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 6);
456 double Q = src.get_val(ParameterType::WILSON, "EW_SCALE", 1);
457 double sw2 = src.get_val(ParameterType::WILSON, "WPARAM_SI_SM", 4);
458
459 double logqt = log(Q * Q / (mtop * mtop));
460
461 // printf("L = %.5f\n", L);
462 // printf("x_t = %.5f\n", xt);
463 // printf("mtop = %.5f\n", mtop);
464 // printf("Q = %.5f\n", Q);
465 // printf("sw2 = %.5f\n", sw2);
466 // printf("logqt = %.5f\n", logqt);
467
468 // printf("C9_NLO = %.5f\n", (1. - 4. * sw2) / sw2 * C1t(xt, logqt)
469 // - B1t(xt, logqt) / sw2
470 // - D1t(xt, logqt)
471 // + 1. / sw2
472 // + 524. / 729.
473 // - 128. / 243. * PI2
474 // - 16. / 3. * L
475 // - 128. / 81. * L * L);
476
477 return (1. - 4. * sw2) / sw2 * C1t(xt, logqt)
478 - B1t(xt, logqt) / sw2
479 - D1t(xt, logqt)
480 + 1. / sw2
481 + 524. / 729.
482 - 128. / 243. * PI2
483 - 16. / 3. * L
484 - 128. / 81. * L * L;
485}
486
487// ---------- C10 ----------
488
489C10::C10() : WilsonCoefficient("C10", GroupMapper::str(WGroup::B, ScaleType::MATCHING)) {
490 // LO
491 matching_info[QCDOrder::LO] = {
492 {
493 {"WPARAM_MATCH_SM", LhaID(2, 1)}, // x_t
494 {"WPARAM_SI_SM", 4} // sw2
495 },
496 compute_LO,
497 get_lhaid_from_name(QCDOrder::LO)
498 };
499
500 // NLO
501 matching_info[QCDOrder::NLO] = {
502 {
503 {"WPARAM_MATCH_SM", 6}, // mtop at muW
504 {"WPARAM_MATCH_SM", LhaID(2, 1)}, // x_t
505 {"WPARAM_SI_SM", 4}, // sw2
506 {"EW_SCALE", 1} // Q_match
507 },
508 compute_NLO,
509 get_lhaid_from_name(QCDOrder::NLO)
510 };
511
512 // NNLO
513 matching_info[QCDOrder::NNLO] = {
514 {
515 {"WPARAM_MATCH_SM", 6}, // mtop_muW
516 {"WPARAM_MATCH_SM", LhaID(2, 1)}, // x_t
517 {"WPARAM_MATCH_SM", 7}, // xtW
518 {"WPARAM_MATCH_SM", 8}, // xtt
519 {"WPARAM_SI_SM", 4}, // sw2
520 {"EW_SCALE", 1}, // Q_match
521 {ParameterType::SM, "MASS", 24} // m_W
522 },
523 compute_NNLO,
524 get_lhaid_from_name(QCDOrder::NNLO)
525 };
526}
527
528double C10::compute_LO(const ParamSrc& src) {
529 double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
530 double sw2 = src.get_val(ParameterType::WILSON, "WPARAM_SI_SM", 4);
531
532 // printf("x_t = %.5f\n", xt);
533 // printf("sw2 = %.5f\n", sw2);
534
535 return (B0t(xt) - C0t(xt) - 0.25) / sw2;
536}
537
538double C10::compute_NLO(const ParamSrc& src) {
539 double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
540 double mtop = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 6);
541 double Q = src.get_val(ParameterType::WILSON, "EW_SCALE", 1);
542 double sw2 = src.get_val(ParameterType::WILSON, "WPARAM_SI_SM", 4);
543
544 // printf("x_t = %.5f\n", xt);
545 // printf("mtop = %.5f\n", mtop);
546 // printf("mu_W = %.5f\n", Q);
547 // printf("sw2 = %.5f\n", sw2);
548
549 double logqt = log(Q * Q / (mtop * mtop));
550 return (B1t(xt, logqt) - C1t(xt, logqt)) / sw2 - 1. / sw2;
551}
552
553double C10::compute_NNLO(const ParamSrc& src) {
554 double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});
555 double xtW = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 7);
556 double xtt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 8);
557 double mtop = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", 6);
558 double Q = src.get_val(ParameterType::WILSON, "EW_SCALE", 1);
559 double mW = src.get_val(ParameterType::SM, "MASS", 24);
560 double sw2 = src.get_val(ParameterType::WILSON, "WPARAM_SI_SM", 4);
561
562 // printf("x_t = %.5f\n", xt);
563 // printf("mtop = %.5f\n", mtop);
564 // printf("mu_W = %.5f\n", Q);
565 // printf("sw2 = %.5f\n", sw2);
566 // printf("xtW = %.5f\n", xtW);
567 // printf("xtt = %.5f\n", xtt);
568 // printf("mW = %.5f\n", mW);
569
570 double logqt = log(Q * Q / (mtop * mtop));
571 double logqW = log(Q * Q / (mW * mW));
572
573 double coeff_temp =
574 (C10Wt2mt(xtt)
575 + logqt * (
576 (69. + 1292. * xt - 209. * xt * xt) / 18. / pow(xt - 1., 3.)
577 - (521. * xt + 105. * xt * xt - 50. * pow(xt, 3.)) / 9. / pow(xt - 1., 4.) * log(xt)
578 - (47. * xt + xt * xt) / 3. / pow(xt - 1., 3.) * Li2(1. - 1. / xt)
579 + logqt * (
580 (61. * xt + 11. * xt * xt) / 3. / pow(xt - 1., 3.)
581 - (49. * xt + 96. * xt * xt - pow(xt, 3.)) / 6. / pow(xt - 1., 4.) * log(xt)
582 )
583 )
584 - (C10Wc2MW(xtW) - 23. / 6. * logqW)
585 + C10Zt2mt(xtt)
586 + logqt * (
587 (188. * xt + 4. * xt * xt + 95. * pow(xt, 3.) - 47. * pow(xt, 4.)) / 6. / pow(xt - 1., 3.) * Li2(1. - 1. / xt)
588 + (1468. * xt + 1578. * xt * xt - 25. * pow(xt, 3.) - 141. * pow(xt, 4.)) / 18. / pow(xt - 1., 4.) * log(xt)
589 - (4622. * xt + 1031. * xt * xt + 582. * pow(xt, 3.) - 475. * pow(xt, 4.)) / 36. / pow(xt - 1., 3.)
590 + logqt * (
591 (49. * xt + 315. * xt * xt - 4. * pow(xt, 3.)) / 6. / pow(xt - 1., 4.) * log(xt)
592 - (440. * xt + 257. * xt * xt + 72. * pow(xt, 3.) - 49. * pow(xt, 4.)) / 12. / pow(xt - 1., 3.)
593 )
594 )
595 + C10Z2tri(xtt)
596 ) * (-2. / sw2);
597
598 return coeff_temp;
599}
600
601
ScaleType
ScaleType
Definition GeneralEnum.h:666
WGroup
WGroup
Definition GeneralEnum.h:618
C10::compute_NNLO
static double compute_NNLO(const ParamSrc &src)
Definition BWilson.cpp:553
C10::compute_NLO
static double compute_NLO(const ParamSrc &src)
Definition BWilson.cpp:538
C10::C10
C10()
Definition BWilson.cpp:489
C10::compute_LO
static double compute_LO(const ParamSrc &src)
Definition BWilson.cpp:528
C1::C1
C1()
Definition BWilson.cpp:5
C2::C2
C2()
Definition BWilson.cpp:33
C2::compute_LO
static double compute_LO(const ParamSrc &src)
Definition BWilson.cpp:51
C2::compute_NNLO
static double compute_NNLO(const ParamSrc &src)
Definition BWilson.cpp:55
C3::C3
C3()
Definition BWilson.cpp:63
C3::compute_NNLO
static double compute_NNLO(const ParamSrc &src)
Definition BWilson.cpp:81
C4::compute_NLO
static double compute_NLO(const ParamSrc &src)
Definition BWilson.cpp:125
C4::C4
C4()
Definition BWilson.cpp:98
C4::compute_NNLO
static double compute_NNLO(const ParamSrc &src)
Definition BWilson.cpp:131
C5::compute_NNLO
static double compute_NNLO(const ParamSrc &src)
Definition BWilson.cpp:163
C5::C5
C5()
Definition BWilson.cpp:145
C6::compute_NNLO
static double compute_NNLO(const ParamSrc &src)
Definition BWilson.cpp:200
C6::C6
C6()
Definition BWilson.cpp:182
C7::C7
C7()
Definition BWilson.cpp:220
C7::compute_LO
static double compute_LO(const ParamSrc &src)
Definition BWilson.cpp:259
C7::compute_NLO
static double compute_NLO(const ParamSrc &src)
Definition BWilson.cpp:264
C7::compute_NNLO
static double compute_NNLO(const ParamSrc &src)
Definition BWilson.cpp:275
C8::compute_LO
static double compute_LO(const ParamSrc &src)
Definition BWilson.cpp:347
C8::compute_NLO
static double compute_NLO(const ParamSrc &src)
Definition BWilson.cpp:352
C8::C8
C8()
Definition BWilson.cpp:309
C8::compute_NNLO
static double compute_NNLO(const ParamSrc &src)
Definition BWilson.cpp:364
C9::C9
C9()
Definition BWilson.cpp:399
C9::compute_NLO
static double compute_NLO(const ParamSrc &src)
Definition BWilson.cpp:452
C9::compute_LO
static double compute_LO(const ParamSrc &src)
Definition BWilson.cpp:428
GroupMapper
Mapper for WGroup <-> WGroupId <-> optional external string.
Definition wgroup_ids.hpp:65
ParamSrc
Lightweight view over a set of source parameters keyed by ParamId.
Definition SourcesView.h:249
ParamSrc::get_val
scalar_t get_val(const ParamId &id) const
Retrieves the current value of a parameter.
Definition SourceView.cpp:111
WilsonCoefficient
Abstract base class representing a Wilson coefficient and its matching information.
Definition Wilson.h:153
WilsonCoefficient::get_lhaid_from_name
LhaID get_lhaid_from_name(QCDOrder order)
Computes the LhaID directly from the coefficient base name and mapping conventions.
Definition Wilson.cpp:61
WilsonCoefficient::matching_info
std::map< QCDOrder, MatchingInfo > matching_info
Matching metadata indexed by QCD order.
Definition Wilson.h:330
PI2
constexpr double PI2
Definition constants.h:8
Li2
double Li2(double x)
Computes the dilogarithm function Li2(x).
Definition polylog.cpp:1257
pow
scalar_t pow(const scalar_t &base, const scalar_t &exp)
Definition scalar.cpp:75
B1t
double B1t(double x, double l)
Wilson coefficient B1t depending on x and scale l.
Definition special_SM.cpp:796
C8c2MW
double C8c2MW(double x)
Wilson coefficient C8c at M_W scale.
Definition special_SM.cpp:1090
T
double T(double x)
Wilson coefficient T(x).
Definition special_SM.cpp:866
C10Wt2mt
double C10Wt2mt(double x)
Wilson coefficient C10Wt at 2m_t scale.
Definition special_SM.cpp:1104
G1t
double G1t(double x, double l)
Wilson coefficient G1t depending on x and scale l.
Definition special_SM.cpp:845
A1t
double A1t(double x, double l)
Wilson coefficient A1t depending on x and scale l.
Definition special_SM.cpp:765
C0t
double C0t(double x)
Scalar one-loop Wilson coefficient C0t.
Definition special_SM.cpp:782
D1t
double D1t(double x, double l)
Wilson coefficient D1t depending on x and scale l.
Definition special_SM.cpp:812
C1t
double C1t(double x, double l)
Wilson coefficient C1t depending on x and scale l.
Definition special_SM.cpp:805
C7c2MW
double C7c2MW(double x)
Wilson coefficient C7c at M_W scale.
Definition special_SM.cpp:1057
E0t
double E0t(double x)
Scalar one-loop Wilson coefficient E0t.
Definition special_SM.cpp:838
C8t2mt
double C8t2mt(double x)
Wilson coefficient C8t at 2m_t scale.
Definition special_SM.cpp:1073
A0t
double A0t(double x)
Scalar one-loop Wilson coefficient A0t.
Definition special_SM.cpp:758
C10Z2tri
double C10Z2tri(double x)
Wilson coefficient C10Z from triangle diagrams.
Definition special_SM.cpp:1138
B0t
double B0t(double x)
Scalar one-loop Wilson coefficient B0t.
Definition special_SM.cpp:775
F0t
double F0t(double x)
Scalar one-loop Wilson coefficient F0t.
Definition special_SM.cpp:820
D0t
double D0t(double x)
Scalar one-loop Wilson coefficient D0t.
Definition special_SM.cpp:789
F1t
double F1t(double x, double l)
Wilson coefficient F1t depending on x and scale l.
Definition special_SM.cpp:828
C10Zt2mt
double C10Zt2mt(double x)
Wilson coefficient C10Zt at 2m_t scale.
Definition special_SM.cpp:1127
C10Wc2MW
double C10Wc2MW(double x)
Wilson coefficient C10Wc at M_W scale.
Definition special_SM.cpp:1116
C7t2mt
double C7t2mt(double x)
Wilson coefficient C7t at 2m_t scale.
Definition special_SM.cpp:1042
E1t
double E1t(double x, double l)
Wilson coefficient E1t depending on x and scale l.
Definition special_SM.cpp:855
LhaID
Represents an identifier of a LHA element, possibly containing several sub-ids.
Definition LhaID.h:56
MatchingInfo
Matching metadata for a Wilson coefficient at a specific QCD order.
Definition Wilson.h:62