a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
SUSYMassInsertionMatching.cpp
Go to the documentation of this file.
1/*
2 * Copyright (C) 2012 HEPfit Collaboration
3 *
4 *
5 * For the licensing terms see doc/COPYING.
6 */
7
8#include "SUSYMassInsertionMatching.h"
9#include "SUSYMassInsertion.h"
10
11SUSYMassInsertionMatching::SUSYMassInsertionMatching(const SUSYMassInsertion & SUSYMassInsertion_i):
12
13 StandardModelMatching(SUSYMassInsertion_i),
14 SusyMI(SUSYMassInsertion_i),
15
16 drNDRLRI(5, 5, 0.),
17 mcd2(5, NDR, NLO),
18 mcd1(10, NDR, NLO),
19 mcbd(5, NDR, NLO), mcbs(5, NDR, NLO), mck2(5, NDR, NLO) {
20
21 Nf = 6;
22
23 double Nc = SusyMI.getNc();
24
25 drNDRLRI(0,0) = -(((-1. + Nc)*(-7. + log(4096.))) / Nc);
26 drNDRLRI(1,1) = (-2.*(-1. + 6.*Nc*Nc - 8.*log(2.) + Nc*(-13. + log(1024.))))/(3.*Nc);
27 drNDRLRI(1,2) = (-2.*(13. - 10.*log(2.) + Nc*(-5. + log(256.))))/(3.*Nc);
28 drNDRLRI(2,1) = (-8. + 6.*Nc*Nc + 20.*log(2.) - 8.*Nc*(1. + log(4.)))/(3.*Nc);
29 drNDRLRI(2,2) = (2.*(4. + Nc - 10.*Nc*log(2.) + log(256.)))/(3.*Nc);
30 drNDRLRI(3,3) = (2. - 4.*Nc*Nc + log(4.))/Nc;
31 drNDRLRI(3,4) = 2. - log(4.);
32 drNDRLRI(4,3) = -2.*(1. + log(2.));
33 drNDRLRI(4,4) = (2. + log(4.))/Nc;
34
35}
36
37// delta F = 1 loop functions
38
39double SUSYMassInsertionMatching::B1(double x) const {
40 return ((1. + 4.*x - 5.*x*x + 4.*x*log(x) + 2.*x*x*log(x))/(8.*pow(x-1.,4.)));
41}
42
43double SUSYMassInsertionMatching::B2(double x) const {
44 return (x*(5. - 4.*x - x*x + 2.*log(x) + 4.*x*log(x))/(2.*pow(x-1.,4.)));
45}
46
47double SUSYMassInsertionMatching::P1(double x) const {
48 return ((1. - 6.*x + 18.*x*x - 10.*x*x*x - 3.*x*x*x*x +
49 12.*x*x*x*log(x))/(18.*pow(x-1.,5.)));
50}
51
52double SUSYMassInsertionMatching::P2(double x) const {
53 return ((7. - 18.*x + 9.*x*x + 2.*x*x*x +
54 3.*log(x) - 9.*x*x*log(x))/(9.*pow(x-1.,5.)));
55}
56
57double SUSYMassInsertionMatching::M1(double x) const{
58 return((1. + 4. * x - 5. * x*x + 4.*x*log(x) + 2.*x*x*log(x))/(2.*pow(x-1.,4.)));
59}
60
61double SUSYMassInsertionMatching::M2(double x) const{
62 return((-5. + 4.*x + x*x - 2.*log(x) - 4.*x*log(x))/(2.*pow(x-1.,4.)));
63}
64
65double SUSYMassInsertionMatching::M3(double x) const{
66 return((-1. + 9.*x + 9.*x*x - 17.*x*x*x + 18.*x*x*log(x) + 6.*x*x*x*log(x))/(12.*pow(x-1.,5.)));
67}
68
69double SUSYMassInsertionMatching::M4(double x) const {
70 return((-1.-9.*x + 9.*x*x + x*x*x - 6.*x*log(x) - 6.*x*x* log(x))/(6.*pow(x-1.,5.)));
71}
72
73// delta F = 2 loop functions
74
75double SUSYMassInsertionMatching::f6(double x){
76 return(x*(17. - 9.*x - 9.*x*x + x*x*x + 18.*x*log(x) +
77 6.*log(x))/(6.*pow(x-1.,5.)));
78}
79
80double SUSYMassInsertionMatching::f6t(double x){
81 return((1. + 9.*x - 9.*x*x - x*x*x + 6.*x*x*log(x) +
82 6.*x*log(x))/(3.*pow(x-1.,5.)));
83}
84
85//the LO terms have to be multiplyed by als^2(mumatch)/Msq^2
91double SUSYMassInsertionMatching::C0LO(double x){
92 //return(16./pow(1.-x,5.)*(11./1728. + 133.*x/1728. - 13.*x*x/192. - 29.*x*x*x/1728. +
93 // x*x*x*x/864 + 13.*x*log(x)/288. + 17.*x*x*log(x)/288.));
94 double pow5 = pow(1.-x,5.);
95 return( (1./pow5) *
96 (11./108.+133./108.*x-13./12.*x*x-29./108.*x*x*x+1./54.*x*x*x*x+(13./18.*x+17./18.*x*x)*log(x)
97 ) );
98}
99
100
101double SUSYMassInsertionMatching::C1LO(double x){
102 return(16.*(289.*x/1728. - 17.*x*x/192. - 17.*x*x*x/192. + 17.*x*x*x*x/1728. +
103 17.*x*log(x)/288. + 17.*x*x*log(x)/96. )/pow(1.-x,5.));
104}
105
106double SUSYMassInsertionMatching::C2LO(double x){
107 return(16.*(-17.*x/576. + x*x/64. + x*x*x/64. - x*x*x*x/576. - x*log(x)/96. -
108 x*x*log(x)/32.)/(pow(1.-x,5)));
109}
110
111double SUSYMassInsertionMatching::C3LOA(double x){
112 return(16.*(-11./864. - 11.*x/96. + 11.*x*x/96. + 11.*x*x*x/864. - 11.*x*log(x)/144. -
113 11.*x*x*log(x)/144.)/pow(1.-x,5.) );
114}
115
116double SUSYMassInsertionMatching::C3LOB(double x){
117 return(16.*(-1./144. + 101.*x/288. - 5.*x*x/32. - 61.*x*x*x/288. + 7.*x*x*x*x/288.
118 + 5.*x*log(x)/48. + 19.*x*x*log(x)/48.)/(pow(1.-x,5.)));
119}
120
121double SUSYMassInsertionMatching::C4LOA(double x){
122 return (16.*(-5./288. - 5.*x/32. + 5.*x*x/32. + 5.*x*x*x/288.
123 - 5.*x*log(x)/48. - 5.*x*x*log(x)/48.)/(pow(1.-x,5.))
124 );
125}
126
127double SUSYMassInsertionMatching::C4LOB(double x){
128 return(16.*(5./432. + 107.*x/864. - 11.*x*x/96. - 19.*x*x*x/864. + x*x*x*x/864.
129 + 11.*x*log(x)/144. + 13.*x*x*log(x)/144.)/(pow(1.-x,5.)));
130}
131
132//to multiplyed by als^3(mumatch)/M_PI/Msq^2
138double SUSYMassInsertionMatching::C0NLO(double x, double mumatch2, double Ms2){
139 double Li2 = gsl_sf_dilog(1-x);
140 double a = log(2.);
141 double pow7 = pow(1.-x,7.);
142 double pow6 = pow(1.-x,6.);
143 double pow5 = pow(1.-x,5.);
144 double pow4 = pow(1.-x,4.);
145
146 return ( 64.*
147 ((-11./(20736.*pow7) +
148 4459./(165888.*pow6) + (209.*Nf)/(27648.*
149 pow6) - 1./(1728.*pow5) +
150 151./(41472.*pow4) - (89.*x)/(20736.*
151 pow7) + (13235.*x)/(62208.*pow6) - (83.*
152 Nf*x)/(6912.*pow6) + (3853.*x)/(15552.*
153 pow5) - (173.*x)/(27648.*pow4) + (77.*
154 x*x)/(2592.*pow7) + (33361.*
155 x*x)/(248832.*pow6) - (33.*Nf*
156 x*x)/(512.*pow6) - (9355.*
157 x*x)/(41472.*pow5) + (217.*
158 x*x)/(82944.*pow4) - (419.*
159 x*x*x)/(10368.*pow7) - (7.*
160 x*x*x)/(18.*pow6) + (481.*Nf*
161 x*x*x)/(6912.*pow6) - (503.*
162 x*x*x)/(20736.*pow5) + (233.*
163 x*x*x*x)/(20736.*pow7) + (10085.*
164 x*x*x*x)/(497664.*pow6) - (19.*Nf*
165 x*x*x*x)/(27648.*pow6) + (331.*
166 x*x*x*x)/(124416.*pow5) + (95.*
167 x*x*x*x*x)/(20736.*pow7) - (91.*
168 x*x*x*x*x)/(15552.*pow6) -
169 x*x*x*x*x*x/(3456.*pow7) +
170 x*x*x*x*x*x/(1296.*pow6) - (203.*Li2)/(13824.*
171 pow6) - (11.*Nf*Li2)/(2304.*
172 pow6) + (47.*Li2)/(4608.*
173 pow5) + (53.*x*Li2)/(864.*
174 pow6) - (Nf*x*Li2)/(192.*
175 pow6) + (367.*x*Li2)/(6912.*
176 pow5) - (85.*x*Li2)/(13824.*
177 pow4) - (127.*x*x*Li2)/(13824.*
178 pow6) + (19.*Nf*x*x*Li2)/(768.*
179 pow6) - (35.*x*x*Li2)/(13824.*
180 pow5) - (259.*x*x*x*Li2)/(6912.*
181 pow6) - (17.*Nf*x*x*x*Li2)/(1152.*
182 pow6) + (11.*a)/(3456.*
183 pow5) + (133.*x*a)/(3456.*
184 pow5) - (13.*x*x*a)/(384.*
185 pow5) - (29.*x*x*x*a)/(3456.*
186 pow5) + (x*x*x*x*a)/(1728.*
187 pow5) - (13.*x*log(x))/(3456.*
188 pow7) + (7.*x*log(x))/(256.*
189 pow6) - (Nf*x*log(x))/(1152.*
190 pow6) + (6431.*x*log(x))/(82944.*
191 pow5) - (425.*x*log(x))/(82944.*
192 pow4) + (29.*x*x*log(x))/(2592.*
193 pow7) + (10535.*x*x*log(x))/(41472.*
194 pow6) - (Nf*x*x*log(x))/(48.*
195 pow6) + (1589.*x*x*log(x))/(6912.*
196 pow5) + (55.*x*x*x*log(x))/(2592.*
197 pow7) + (3751.*x*x*x*log(x))/(10368.*
198 pow6) - (25.*Nf*x*x*x*log(x))/(576.*
199 pow6) + (1681.*x*x*x*log(x))/(82944.*
200 pow5) - (5.*x*x*x*x*log(x))/(192.*
201 pow7) - (91.*x*x*x*x*log(x))/(4608.*
202 pow6) - (29.*x*x*x*x*log(x))/(20736.*
203 pow5) - (29.*x*x*x*x*x*log(x))/(10368.*
204 pow7) + (139.*x*x*x*x*x*log(x))/(20736.*
205 pow6) + (x*x*x*x*x*x*log(x))/(5184.*
206 pow7) - (x*x*x*x*x*x*log(x))/(1296.*
207 pow6) + (13.*x*a*log(x))/(576.*
208 pow5) + (17.*x*x*a*log(x))/(576.*
209 pow5) - (73.*x*pow(log(x), 2.))/(3072.*
210 pow6) - (29.*x*pow(log(x), 2.))/(1152.*
211 pow5) - (43.*x*pow(log(x), 2.))/(9216.*
212 pow4) - (1075.*x*x*
213 pow(log(x), 2.))/(4608.*pow6) - (205.*
214 x*x*pow(log(x), 2.))/(2304.*pow5) + (13.*
215 x*x*x*pow(log(x), 2.))/(1728.*
216 pow7) - (3053.*x*x*x*
217 pow(log(x), 2.))/(27648.*pow6) + (17.*
218 x*x*x*x*pow(log(x), 2.))/(1728.*pow7) +
219 log(Ms2/mumatch2)*((11.*Nf)/(13824.*pow6) -
220 55./(6912.*pow5) - (1153.*x)/(6912.*
221 pow6) + (97.*Nf*x)/(3456.*
222 pow6) - (665.*x)/(6912.*
223 pow5) - (1819.*x*x)/(20736.*
224 pow6) + (Nf*x*x)/(192.*
225 pow6) + (65.*x*x)/(768.*
226 pow5) + (2789.*x*x*x)/(10368.*
227 pow6) - (113.*Nf*x*x*x)/(3456.*
228 pow6) + (145.*x*x*x)/(6912.*
229 pow5) - (143.*x*x*x*x)/(6912.*
230 pow6) - (19.*Nf*x*x*x*x)/(13824.*
231 pow6) - (5.*x*x*x*x)/(3456.*
232 pow5) + (145.*x*x*x*x*x)/(20736.*
233 pow6) -
234 x*x*x*x*x*x/(1296.*pow6) - (13.*x*log(x))/(256.*
235 pow6) + (13.*Nf*x*log(x))/(1152.*
236 pow6) - (65.*x*log(x))/(1152.*
237 pow5) - (953.*x*x*log(x))/(3456.*
238 pow6) + (5.*Nf*x*x*log(x))/(128.*
239 pow6) - (85.*x*x*log(x))/(1152.*
240 pow5) - (593.*x*x*x*log(x))/(6912.*
241 pow6) + (17.*Nf*x*x*x*log(x))/(1152.*
242 pow6))))
243
244 );
245}
246
247
248double SUSYMassInsertionMatching::C1NLO(double x, double mumatch2, double Ms2){
249 double Li2 = gsl_sf_dilog(1-x);
250 double a = log(2.);
251 double pow7 = pow(1.-x,7.);
252 double pow6 = pow(1.-x,6.);
253 double pow5 = pow(1.-x,5.);
254
255 return(64.*(((-289.*x)/(20736.*pow7) + (134447.*x)/(248832.*
256 pow6) + (323.*Nf*x)/(6912.*
257 pow6) + (43957.*x)/(82944.*
258 pow5) + (1309.*x*x)/(20736.*
259 pow7) + (35261.*x*x)/(124416.*
260 pow6) - (187.*Nf*x*x)/(768.*
261 pow6) - (30427.*x*x)/(82944.*
262 pow5) - (221.*x*x*x)/(3456.*
263 pow7) - (8107.*x*x*x)/(10368.*
264 pow6) + (51.*Nf*x*x*x)/(256.*
265 pow6) - (15493.*x*x*x)/(82944.*
266 pow5) - (85.*x*x*x*x)/(10368.*
267 pow7) - (1813.*x*x*x*x)/(124416.*
268 pow6) - (17.*Nf*x*x*x*x)/(6912.*
269 pow6) + (1963.*x*x*x*x)/(82944.*
270 pow5) + (527.*x*x*x*x*x)/(20736.*
271 pow7) - (8407.*x*x*x*x*x)/(248832.*
272 pow6) - (17.*x*x*x*x*x*x)/(6912.*
273 pow7) + (17.*x*x*x*x*x*x)/(2592.*
274 pow6) + (89.*x*Li2)/(3456.*
275 pow6) - (17.*Nf*x*Li2)/(384.*
276 pow6) + (2.*x*Li2)/(27.*
277 pow5) + (61.*x*x*Li2)/(576.*
278 pow6) + (17.*Nf*x*x*Li2)/(192.*
279 pow6) - (127.*x*x*Li2)/(864.*
280 pow5) - (455.*x*x*x*Li2)/(3456.*
281 pow6) - (17.*Nf*x*x*x*Li2)/(384.*
282 pow6) + (1411.*x*a)/(31104.*
283 pow5) - (83.*x*x*a)/(3456.*
284 pow5) - (83.*x*x*x*a)/(3456.*
285 pow5) + (83.*x*x*x*x*a)/(31104.*
286 pow5) - (17.*x*log(x))/(3456.*
287 pow7) + (377.*x*log(x))/(3456.*
288 pow6) - (17.*Nf*x*log(x))/(2304.*
289 pow6) + (6443.*x*log(x))/(41472.*
290 pow5) + (17.*x*x*log(x))/(3456.*
291 pow7) + (35.*x*x*log(x))/(81.*
292 pow6) - (17.*Nf*x*x*log(x))/(1152.*
293 pow6) + (16055.*x*x*log(x))/(41472.*
294 pow5) + (187.*x*x*x*log(x))/(2592.*
295 pow7) + (251.*x*x*x*log(x))/(288.*
296 pow6) - (289.*Nf*x*x*x*log(x))/(2304.*
297 pow6) + (35.*x*x*x*log(x))/(576.*
298 pow5) - (17.*x*x*x*x*log(x))/(288.*
299 pow7) - (31.*x*x*x*x*log(x))/(2304.*
300 pow6) - (35.*x*x*x*x*log(x))/(5184.*
301 pow5) - (17.*x*x*x*x*x*log(x))/(1152.*
302 pow7) + (701.*x*x*x*x*x*log(x))/(20736.*
303 pow6) + (17.*x*x*x*x*x*x*log(x))/(10368.*
304 pow7) - (17.*x*x*x*x*x*x*log(x))/(2592.*
305 pow6) + (83.*x*a*log(x))/(5184.*
306 pow5) + (83.*x*x*a*log(x))/(1728.*
307 pow5) - (7.*x*pow(log(x), 2.))/(1536.*
308 pow6) - (17.*x*pow(log(x), 2.))/(1728.*
309 pow5) - (2599.*x*x*
310 pow(log(x), 2.))/(6912.*pow6) - (167.*
311 x*x*pow(log(x), 2.))/(864.*pow5) + (17.*
312 x*x*x*pow(log(x), 2.))/(1728.*
313 pow7) - (4673.*x*x*x*
314 pow(log(x), 2.))/(13824.*pow6) + (17.*
315 x*x*x*x*pow(log(x), 2.))/(576.*pow7) +
316 log(Ms2/mumatch2)*((-12869.*x)/(41472.*pow6) + (85.*Nf*
317 x)/(1728.*pow6) - (697.*x)/(5184.*
318 pow5) - (7429.*x*x)/(20736.*
319 pow6) + (17.*Nf*x*x)/(384.*
320 pow6) + (41.*x*x)/(576.*
321 pow5) + (3859.*x*x*x)/(5184.*
322 pow6) - (17.*Nf*x*x*x)/(192.*
323 pow6) + (41.*x*x*x)/(576.*
324 pow5) - (2363.*x*x*x*x)/(20736.*
325 pow6) - (17.*Nf*x*x*x*x)/(3456.*
326 pow6) - (41.*x*x*x*x)/(5184.*
327 pow5) + (1853.*x*x*x*x*x)/(41472.*
328 pow6) - (17.*x*x*x*x*x*x)/(2592.*
329 pow6) - (595.*x*log(x))/(6912.*
330 pow6) + (17.*Nf*x*log(x))/(1152.*
331 pow6) - (41.*x*log(x))/(864.*
332 pow5) - (2159.*x*x*log(x))/(3456.*
333 pow6) + (17.*Nf*x*x*log(x))/(192.*
334 pow6) - (41.*x*x*log(x))/(288.*
335 pow5) - (1547.*x*x*x*log(x))/(6912.*
336 pow6) + (17.*Nf*x*x*x*log(x))/(384.*
337 pow6))))
338 );
339}
340
341double SUSYMassInsertionMatching::C2NLO(double x, double mumatch2, double Ms2){
342 double Li2 = gsl_sf_dilog(1-x);
343 double a = log(2.);
344 double pow7 = pow(1.-x,7.);
345 double pow6 = pow(1.-x,6.);
346 double pow5 = pow(1.-x,5.);
347
348 return (64.*(((17.*x)/(6912.*pow7) - (10861.*x)/(248832.*
349 pow6) - (19.*Nf*x)/(2304.*
350 pow6) - (18253.*x)/(248832.*
351 pow5) - (77.*x*x)/(6912.*
352 pow7) - (19159.*x*x)/(124416.*
353 pow6) + (11.*Nf*x*x)/(256.*
354 pow6) + (4289.*x*x)/(82944.*
355 pow5) + (13.*x*x*x)/(1152.*
356 pow7) + (659.*x*x*x)/(3456.*
357 pow6) - (9.*Nf*x*x*x)/(256.*
358 pow6) + (77.*x*x*x)/(3072.*
359 pow5) + (5.*x*x*x*x)/(3456.*
360 pow7) + (287.*x*x*x*x)/(124416.*
361 pow6) + (Nf*x*x*x*x)/(2304.*
362 pow6) - (851.*x*x*x*x)/(248832.*
363 pow5) - (31.*x*x*x*x*x)/(6912.*
364 pow7) + (1445.*x*x*x*x*x)/(248832.*
365 pow6) + x*x*x*x*x*x/(2304.*pow7) -
366 x*x*x*x*x*x/(864.*pow6) - (35.*x*Li2)/(3456.*
367 pow6) + (Nf*x*Li2)/(128.*
368 pow6) - (7.*x*Li2)/(432.*
369 pow5) - (13.*x*x*Li2)/(1728.*
370 pow6) - (Nf*x*x*Li2)/(64.*
371 pow6) - (5.*x*x*Li2)/(864.*
372 pow5) + (61.*x*x*x*Li2)/(3456.*
373 pow6) + (Nf*x*x*x*Li2)/(128.*
374 pow6) + (731.*x*a)/(31104.*
375 pow5) - (43.*x*x*a)/(3456.*
376 pow5) - (43.*x*x*x*a)/(3456.*
377 pow5) + (43.*x*x*x*x*a)/(31104.*
378 pow5) + (x*log(x))/(1152.*pow7) + (35.*
379 x*log(x))/(1152.*pow6) + (Nf*x*log(x))/(768.*
380 pow6) - (401.*x*log(x))/(41472.*
381 pow5) - (x*x*log(x))/(1152.*
382 pow7) - (53.*x*x*log(x))/(648.*
383 pow6) + (Nf*x*x*log(x))/(384.*
384 pow6) - (3605.*x*x*log(x))/(41472.*
385 pow5) - (11.*x*x*x*log(x))/(864.*
386 pow7) - (533.*x*x*x*log(x))/(2592.*
387 pow6) + (17.*Nf*x*x*x*log(x))/(768.*
388 pow6) - (11.*x*x*x*log(x))/(576.*
389 pow5) + (x*x*x*x*log(x))/(96.*
390 pow7) + (229.*x*x*x*x*log(x))/(20736.*
391 pow6) + (11.*x*x*x*x*log(x))/(5184.*
392 pow5) + (x*x*x*x*x*log(x))/(384.*
393 pow7) - (143.*x*x*x*x*x*log(x))/(20736.*
394 pow6) - (x*x*x*x*x*x*log(x))/(3456.*
395 pow7) + (x*x*x*x*x*x*log(x))/(864.*
396 pow6) + (43.*x*a*log(x))/(5184.*
397 pow5) + (43.*x*x*a*log(x))/(1728.*
398 pow5) + (7.*x*pow(log(x), 2.))/(512.*
399 pow6) + (11.*x*pow(log(x), 2.))/(1728.*
400 pow5) + (215.*x*x*
401 pow(log(x), 2.))/(2304.*pow6) + (35.*x*x*
402 pow(log(x), 2.))/(864.*pow5) - (x*x*x*
403 pow(log(x), 2.))/(576.*pow7) + (257.*x*x*x*
404 pow(log(x), 2.))/(4608.*pow6) - (x*x*x*x*
405 pow(log(x), 2.))/(192.*pow7) +
406 log(Ms2/mumatch2)*((757.*x)/(13824.*pow6) - (5.*Nf*x)/(576.*
407 pow6) + (187.*x)/(5184.*
408 pow5) + (437.*x*x)/(6912.*
409 pow6) - (Nf*x*x)/(128.*
410 pow6) - (11.*x*x)/(576.*
411 pow5) - (227.*x*x*x)/(1728.*
412 pow6) + (Nf*x*x*x)/(64.*
413 pow6) - (11.*x*x*x)/(576.*
414 pow5) + (139.*x*x*x*x)/(6912.*
415 pow6) + (Nf*x*x*x*x)/(1152.*
416 pow6) + (11.*x*x*x*x)/(5184.*
417 pow5) - (109.*x*x*x*x*x)/(13824.*
418 pow6) +
419 x*x*x*x*x*x/(864.*pow6) + (35.*x*log(x))/(2304.*
420 pow6) - (Nf*x*log(x))/(384.*
421 pow6) + (11.*x*log(x))/(864.*
422 pow5) + (127.*x*x*log(x))/(1152.*
423 pow6) - (Nf*x*x*log(x))/(64.*
424 pow6) + (11.*x*x*log(x))/(288.*
425 pow5) + (91.*x*x*x*log(x))/(2304.*
426 pow6) - (Nf*x*x*x*log(x))/(128.*
427 pow6))))
428
429
430 );
431}
432
433double SUSYMassInsertionMatching::C3NLOA(double x, double mumatch2, double Ms2){
434 double Li2 = gsl_sf_dilog(1-x);
435 double a = log(2.);
436 double pow7 = pow(1.-x,7.);
437 double pow6 = pow(1.-x,6.);
438 double pow5 = pow(1.-x,5.);
439 double pow4 = pow(1.-x,4.);
440
441 return (64.*(11./(10368.*pow7) -
442 8047./(248832.*pow6) - (209.*Nf)/(13824.*
443 pow6) - 15881./(248832.*pow5) -
444 3./(256.*pow4) + (55.*x)/(10368.*
445 pow7) - (16187.*x)/(62208.*
446 pow6) + (121.*Nf*x)/(3456.*
447 pow6) - (38227.*x)/(82944.*
448 pow5) + (35.*x)/(1536.*
449 pow4) - (77.*x*x)/(1728.*
450 pow7) - (6361.*x*x)/(13824.*
451 pow6) + (55.*Nf*x*x)/(768.*
452 pow6) + (43003.*x*x)/(82944.*
453 pow5) - (17.*x*x)/(1536.*
454 pow4) + (341.*x*x*x)/(5184.*
455 pow7) + (56701.*x*x*x)/(62208.*
456 pow6) - (319.*Nf*x*x*x)/(3456.*
457 pow6) + (1553.*x*x*x)/(248832.*
458 pow5) - (253.*x*x*x*x)/(10368.*
459 pow7) - (40567.*x*x*x*x)/(248832.*
460 pow6) + (11.*Nf*x*x*x*x)/(13824.*
461 pow6) - (11.*x*x*x*x*x)/(3456.*
462 pow7) + (11.*x*x*x*x*x)/(2592.*
463 pow6) - (11.*Li2)/(6912.*
464 pow6) + (11.*Nf*Li2)/(1152.*
465 pow6) - (5.*x*Li2)/(108.*
466 pow6) + (61.*x*Li2)/(1728.*
467 pow5) + (x*Li2)/(256.*
468 pow4) - (31.*x*x*Li2)/(6912.*
469 pow6) - (11.*Nf*x*x*
470 Li2)/(384.*pow6) + (125.*x*x*
471 Li2)/(864.*pow5) + (181.*x*x*x*
472 Li2)/(3456.*pow6) + (11.*Nf*
473 x*x*x*Li2)/(576.*pow6) - (17.*
474 a)/(10368.*pow5) - (17.*x*
475 a)/(1152.*pow5) + (17.*x*x*
476 a)/(1152.*pow5) + (17.*x*x*x*
477 a)/(10368.*pow5) + (11.*x*
478 log(x))/(1728.*pow7) + (83.*x*
479 log(x))/(1152.*pow6) - (47.*x*log(x))/(432.*
480 pow5) + (5.*x*log(x))/(1536.*
481 pow4) - (55.*x*x*log(x))/(2592.*
482 pow7) - (1037.*x*x*log(x))/(1728.*
483 pow6) + (11.*Nf*x*x*log(x))/(288.*
484 pow6) - (455.*x*x*log(x))/(1536.*
485 pow5) - (11.*x*x*x*log(x))/(432.*
486 pow7) - (293.*x*x*x*log(x))/(648.*
487 pow6) + (11.*Nf*x*x*x*log(x))/(192.*
488 pow6) - (5.*x*x*x*log(x))/(512.*
489 pow5) + (11.*x*x*x*x*log(x))/(288.*
490 pow7) + (637.*x*x*x*x*log(x))/(10368.*
491 pow6) + (11.*x*x*x*x*x*log(x))/(5184.*
492 pow7) - (11.*x*x*x*x*x*log(x))/(2592.*
493 pow6) - (17.*x*a*log(x))/(1728.*
494 pow5) - (17.*x*x*a*
495 log(x))/(1728.*pow5) + (19.*x*
496 pow(log(x), 2.))/(192.*pow6) + (23.*x*
497 pow(log(x), 2.))/(256.*pow5) + (19.*x*
498 pow(log(x), 2.))/(1536.*pow4) + (337.*
499 x*x*pow(log(x), 2.))/(1152.*
500 pow6) + (45.*x*x*
501 pow(log(x), 2.))/(256.*pow5) - (11.*
502 x*x*x*pow(log(x), 2.))/(864.*
503 pow7) + (589.*x*x*x*
504 pow(log(x), 2.))/(3456.*pow6) - (11.*
505 x*x*x*x*pow(log(x), 2.))/(864.*
506 pow7) +
507 log(Ms2/mumatch2)*(11./(2592.*pow6) - (11.*Nf)/(6912.*
508 pow6) +
509 65./(13824.*pow5) + (1595.*x)/(5184.*
510 pow6) - (77.*Nf*x)/(1728.*
511 pow6) + (65.*x)/(1536.*
512 pow5) - (55.*x*x)/(10368.*
513 pow6) - (65.*x*x)/(1536.*
514 pow5) - (1705.*x*x*x)/(5184.*
515 pow6) + (77.*Nf*x*x*x)/(1728.*
516 pow6) - (65.*x*x*x)/(13824.*
517 pow5) + (275.*x*x*x*x)/(10368.*
518 pow6) + (11.*Nf*x*x*x*x)/(6912.*
519 pow6) - (11.*x*x*x*x*x)/(2592.*
520 pow6) + (385.*x*log(x))/(3456.*
521 pow6) - (11.*Nf*x*log(x))/(576.*
522 pow6) + (65.*x*log(x))/(2304.*
523 pow5) + (715.*x*x*
524 log(x))/(1728.*pow6) - (11.*Nf*
525 x*x*log(x))/(192.*
526 pow6) + (65.*x*x*
527 log(x))/(2304.*pow5) + (275.*
528 x*x*x*log(x))/(3456.*
529 pow6) - (11.*Nf*x*x*x*
530 log(x))/(576.*pow6)))
531 );
532}
533
534double SUSYMassInsertionMatching::C3NLOB(double x, double mumatch2, double Ms2){
535 double Li2 = gsl_sf_dilog(1-x);
536 double a = log(2.);
537 double pow7 = pow(1.-x,7.);
538 double pow6 = pow(1.-x,6.);
539 double pow5 = pow(1.-x,5.);
540 double pow4 = pow(1.-x,4.);
541
542 return (64.*
543
544 (((5./(3456.*pow7) -
545 1465./(82944.*pow6) - (95.*Nf)/(4608.*
546 pow6) - 4151./(82944.*pow5) -
547 5./(768.*pow4) + (25.*x)/(3456.*
548 pow7) - (11465.*x)/(20736.*
549 pow6) + (55.*Nf*x)/(1152.*
550 pow6) - (12493.*x)/(27648.*
551 pow5) + (5.*x)/(512.*pow4) - (35.*
552 pow(x, 2.))/(576.*pow7) - (535.*
553 pow(x, 2.))/(4608.*pow6) + (25.*Nf*
554 pow(x, 2.))/(256.*pow6) + (13093.*
555 pow(x, 2.))/(27648.*pow5) - (5.*
556 pow(x, 2.))/(1536.*pow4) + (155.*
557 pow(x, 3.))/(1728.*pow7) + (14695.*
558 pow(x, 3.))/(20736.*pow6) - (145.*Nf*
559 pow(x, 3.))/(1152.*pow6) + (2351.*
560 pow(x, 3.))/(82944.*pow5) - (115.*
561 pow(x, 4.))/(3456.*pow7) - (2305.*
562 pow(x, 4.))/(82944.*pow6) + (5.*Nf*
563 pow(x, 4.))/(4608.*pow6) - (5.*
564 pow(x, 5.))/(1152.*pow7) + (5.*
565 pow(x, 5.))/(864.*pow6) - (5.*
566 Li2)/(2304.*pow6) + (5.*Nf*
567 Li2)/(384.*pow6) - (5.*x*
568 Li2)/(144.*pow6) - (95.*x*
569 Li2)/(576.*pow5) + (5.*x*
570 Li2)/(256.*pow4) - (145.*pow(x, 2.)*
571 Li2)/(2304.*pow6) - (5.*Nf*
572 pow(x, 2.)*Li2)/(128.*pow6) + (5.*
573 pow(x, 2.)*Li2)/(72.*pow5) + (115.*
574 pow(x, 3.)*Li2)/(1152.*pow6) + (5.*
575 Nf*pow(x, 3.)*Li2)/(192.*pow6) -
576 a/(1152.*pow5) - (x*a)/(128.*
577 pow5) + (pow(x, 2.)*a)/(128.*
578 pow5) + (pow(x, 3.)*a)/(1152.*
579 pow5) + (5.*x*log(x))/(576.*
580 pow7) + (5.*x*log(x))/(384.*
581 pow6) - (139.*x*log(x))/(576.*
582 pow5) + (25.*x*log(x))/(1536.*
583 pow4) - (25.*pow(x, 2.)*log(x))/(864.*
584 pow7) - (95.*pow(x, 2.)*log(x))/(144.*
585 pow6) + (5.*Nf*pow(x, 2.)*log(x))/(96.*
586 pow6) - (1897.*pow(x, 2.)*log(x))/(4608.*
587 pow5) - (5.*pow(x, 3.)*log(x))/(144.*
588 pow7) - (545.*pow(x, 3.)*log(x))/(864.*
589 pow6) + (5.*Nf*pow(x, 3.)*log(x))/(64.*
590 pow6) - (35.*pow(x, 3.)*log(x))/(1536.*
591 pow5) + (5.*pow(x, 4.)*log(x))/(96.*
592 pow7) + (85.*pow(x, 4.)*log(x))/(3456.*
593 pow6) + (5.*pow(x, 5.)*log(x))/(1728.*
594 pow7) - (5.*pow(x, 5.)*log(x))/(864.*
595 pow6) - (x*a*log(x))/(192.*
596 pow5) - (pow(x, 2.)*a*log(x))/(192.*
597 pow5) + (5.*x*pow(log(x), 2.))/(64.*
598 pow6) + (15.*x*pow(log(x), 2.))/(256.*
599 pow5) + (5.*x*pow(log(x), 2.))/(512.*
600 pow4) + (175.*pow(x, 2.)*
601 pow(log(x), 2.))/(384.*pow6) + (45.*
602 pow(x, 2.)*pow(log(x), 2.))/(256.*
603 pow5) - (5.*pow(x, 3.)*
604 pow(log(x), 2.))/(288.*pow7) + (235.*
605 pow(x, 3.)*pow(log(x), 2.))/(1152.*
606 pow6) - (5.*pow(x, 4.)*
607 pow(log(x), 2.))/(288.*pow7) +
608 log(Ms2/mumatch2)*(5./(864.*pow6) - (5.*Nf)/(2304.*
609 pow6) +
610 95./(4608.*pow5) + (725.*x)/(1728.*
611 pow6) - (35.*Nf*x)/(576.*
612 pow6) + (95.*x)/(512.*
613 pow5) - (25.*pow(x, 2.))/(3456.*
614 pow6) - (95.*pow(x, 2.))/(512.*
615 pow5) - (775.*pow(x, 3.))/(1728.*
616 pow6) + (35.*Nf*pow(x, 3.))/(576.*
617 pow6) - (95.*pow(x, 3.))/(4608.*
618 pow5) + (125.*pow(x, 4.))/(3456.*
619 pow6) + (5.*Nf*pow(x, 4.))/(2304.*
620 pow6) - (5.*pow(x, 5.))/(864.*
621 pow6) + (175.*x*log(x))/(1152.*
622 pow6) - (5.*Nf*x*log(x))/(192.*
623 pow6) + (95.*x*log(x))/(768.*
624 pow5) + (325.*pow(x, 2.)*
625 log(x))/(576.*pow6) - (5.*Nf*
626 pow(x, 2.)*log(x))/(64.*pow6) + (95.*
627 pow(x, 2.)*log(x))/(768.*
628 pow5) + (125.*pow(x, 3.)*
629 log(x))/(1152.*pow6) - (5.*Nf*
630 pow(x, 3.)*log(x))/(192.*pow6)))))
631
632 );
633}
634
635double SUSYMassInsertionMatching::C4NLOA(double x, double mumatch2, double Ms2){
636 double Li2 = gsl_sf_dilog(1.-x);
637 double a = log(2.);
638 double pow7 = pow(1.-x,7.);
639 double pow6 = pow(1.-x,6.);
640 double pow5 = pow(1.-x,5.);
641 double pow4 = pow(1.-x,4.);
642
643 return (
644 (5./(3456.*pow7) -
645 1465./(82944.*pow6) - (95.*Nf)/(4608.*
646 pow6) - 4151./(82944.*pow5) -
647 5./(768.*pow4) + (25.*x)/(3456.*
648 pow7) - (11465.*x)/(20736.*
649 pow6) + (55.*Nf*x)/(1152.*
650 pow6) - (12493.*x)/(27648.*
651 pow5) + (5.*x)/(512.*pow4) - (35.*
652 x*x)/(576.*pow7) - (535.*
653 x*x)/(4608.*pow6) + (25.*Nf*
654 x*x)/(256.*pow6) + (13093.*
655 x*x)/(27648.*pow5) - (5.*
656 x*x)/(1536.*pow4) + (155.*
657 x*x*x)/(1728.*pow7) + (14695.*
658 x*x*x)/(20736.*pow6) - (145.*Nf*
659 x*x*x)/(1152.*pow6) + (2351.*
660 x*x*x)/(82944.*pow5) - (115.*
661 x*x*x*x)/(3456.*pow7) - (2305.*
662 x*x*x*x)/(82944.*pow6) + (5.*Nf*
663 x*x*x*x)/(4608.*pow6) - (5.*
664 x*x*x*x*x)/(1152.*pow7) + (5.*
665 x*x*x*x*x)/(864.*pow6) - (5.*
666 Li2)/(2304.*pow6) + (5.*Nf*
667 Li2)/(384.*pow6) - (5.*x*
668 Li2)/(144.*pow6) - (95.*x*
669 Li2)/(576.*pow5) + (5.*x*
670 Li2)/(256.*pow4) - (145.*x*x*
671 Li2)/(2304.*pow6) - (5.*Nf*
672 x*x*Li2)/(128.*pow6) + (5.*
673 x*x*Li2)/(72.*pow5) + (115.*
674 x*x*x*Li2)/(1152.*pow6) + (5.*
675 Nf*x*x*x*Li2)/(192.*pow6) -
676 a/(1152.*pow5) - (x*a)/(128.*
677 pow5) + (x*x*a)/(128.*
678 pow5) + (x*x*x*a)/(1152.*
679 pow5) + (5.*x*log(x))/(576.*
680 pow7) + (5.*x*log(x))/(384.*
681 pow6) - (139.*x*log(x))/(576.*
682 pow5) + (25.*x*log(x))/(1536.*
683 pow4) - (25.*x*x*log(x))/(864.*
684 pow7) - (95.*x*x*log(x))/(144.*
685 pow6) + (5.*Nf*x*x*log(x))/(96.*
686 pow6) - (1897.*x*x*log(x))/(4608.*
687 pow5) - (5.*x*x*x*log(x))/(144.*
688 pow7) - (545.*x*x*x*log(x))/(864.*
689 pow6) + (5.*Nf*x*x*x*log(x))/(64.*
690 pow6) - (35.*x*x*x*log(x))/(1536.*
691 pow5) + (5.*x*x*x*x*log(x))/(96.*
692 pow7) + (85.*x*x*x*x*log(x))/(3456.*
693 pow6) + (5.*x*x*x*x*x*log(x))/(1728.*
694 pow7) - (5.*x*x*x*x*x*log(x))/(864.*
695 pow6) - (x*a*log(x))/(192.*
696 pow5) - (x*x*a*log(x))/(192.*
697 pow5) + (5.*x*pow(log(x), 2.))/(64.*
698 pow6) + (15.*x*pow(log(x), 2.))/(256.*
699 pow5) + (5.*x*pow(log(x), 2.))/(512.*
700 pow4) + (175.*x*x*
701 pow(log(x), 2.))/(384.*pow6) + (45.*
702 x*x*pow(log(x), 2.))/(256.*
703 pow5) - (5.*x*x*x*
704 pow(log(x), 2.))/(288.*pow7) + (235.*
705 x*x*x*pow(log(x), 2.))/(1152.*
706 pow6) - (5.*x*x*x*x*
707 pow(log(x), 2.))/(288.*pow7) +
708 log(Ms2/mumatch2)*(5./(864.*pow6) - (5.*Nf)/(2304.*
709 pow6) +
710 95./(4608.*pow5) + (725.*x)/(1728.*
711 pow6) - (35.*Nf*x)/(576.*
712 pow6) + (95.*x)/(512.*
713 pow5) - (25.*x*x)/(3456.*
714 pow6) - (95.*x*x)/(512.*
715 pow5) - (775.*x*x*x)/(1728.*
716 pow6) + (35.*Nf*x*x*x)/(576.*
717 pow6) - (95.*x*x*x)/(4608.*
718 pow5) + (125.*x*x*x*x)/(3456.*
719 pow6) + (5.*Nf*x*x*x*x)/(2304.*
720 pow6) - (5.*x*x*x*x*x)/(864.*
721 pow6) + (175.*x*log(x))/(1152.*
722 pow6) - (5.*Nf*x*log(x))/(192.*
723 pow6) + (95.*x*log(x))/(768.*
724 pow5) + (325.*x*x*
725 log(x))/(576.*pow6) - (5.*Nf*
726 x*x*log(x))/(64.*pow6) + (95.*
727 x*x*log(x))/(768.*
728 pow5) + (125.*x*x*x*
729 log(x))/(1152.*pow6) - (5.*Nf*
730 x*x*x*log(x))/(192.*pow6)))
731 );
732}
733
734double SUSYMassInsertionMatching::C4NLOB(double x, double mumatch2, double Ms2){
735 double Li2 = gsl_sf_dilog(1.-x);
736 double a = log(2.);
737 double pow7 = pow(1.-x,7.);
738 double pow6 = pow(1.-x,6.);
739 double pow5 = pow(1.-x,5.);
740 double pow4 = pow(1.-x,4.);
741
742 return (64.*(-5./(5184.*pow7) +
743 149./(1944.*pow6) + (95.*Nf)/(6912.*
744 pow6) - 3653./(124416.*pow5) +
745 29./(10368.*pow4) - (67.*x)/(10368.*
746 pow7) + (3971.*x)/(15552.*
747 pow6) - (91.*Nf*x)/(3456.*
748 pow6) + (42647.*x)/(82944.*
749 pow5) - (17.*x)/(3456.*
750 pow4) + (497.*x*x)/(10368.*
751 pow7) + (2059.*x*x)/(6912.*
752 pow6) - (3.*Nf*x*x)/(32.*
753 pow6) - (1397.*x*x)/(3072.*
754 pow5) + (11.*x*x)/(5184.*
755 pow4) - (349.*x*x*x)/(5184.*
756 pow7) - (41527.*x*x*x)/(62208.*
757 pow6) + (371.*Nf*x*x*x)/(3456.*
758 pow6) - (8525.*x*x*x)/(248832.*
759 pow5) + (55.*x*x*x*x)/(2592.*
760 pow7) + (2773.*x*x*x*x)/(62208.*
761 pow6) - (7.*Nf*x*x*x*x)/(6912.*
762 pow6) + (349.*x*x*x*x)/(82944.*
763 pow5) + (61.*x*x*x*x*x)/(10368.*
764 pow7) - (53.*x*x*x*x*x)/(6912.*
765 pow6) -
766 x*x*x*x*x*x/(3456.*pow7) +
767 x*x*x*x*x*x/(1296.*pow6) - (19.*Li2)/(432.*
768 pow6) - (5.*Nf*Li2)/(576.*
769 pow6) + (41.*Li2)/(1152.*
770 pow5) + (599.*x*Li2)/(3456.*
771 pow6) - (Nf*x*Li2)/(192.*
772 pow6) + (13.*x*Li2)/(576.*
773 pow5) - (7.*x*Li2)/(1728.*
774 pow4) - (179.*x*x*
775 Li2)/(1728.*pow6) + (7.*Nf*
776 x*x*Li2)/(192.*pow6) + (7.*
777 x*x*Li2)/(384.*pow5) - (89.*
778 x*x*x*Li2)/(3456.*pow6) - (13.*
779 Nf*x*x*x*Li2)/(576.*
780 pow6) - (7.*a)/(5184.*
781 pow5) + (401.*x*a)/(10368.*
782 pow5) - (17.*x*x*a)/(1152.*
783 pow5) - (265.*x*x*x*a)/(10368.*
784 pow5) + (31.*x*x*x*x*a)/(10368.*
785 pow5) - (11.*x*log(x))/(1728.*
786 pow7) + (41.*x*log(x))/(768.*
787 pow6) - (Nf*x*log(x))/(1152.*
788 pow6) + (5075.*x*log(x))/(41472.*
789 pow5) - (35.*x*log(x))/(10368.*
790 pow4) + (103.*x*x*log(x))/(5184.*
791 pow7) + (6985.*x*x*log(x))/(20736.*
792 pow6) - (7.*Nf*x*x*log(x))/(192.*
793 pow6) + (5741.*x*x*log(x))/(13824.*
794 pow5) + (41.*x*x*x*log(x))/(1296.*
795 pow7) + (4411.*x*x*x*log(x))/(6912.*
796 pow6) - (77.*Nf*x*x*x*log(x))/(1152.*
797 pow6) + (443.*x*x*x*log(x))/(20736.*
798 pow5) - (x*x*x*x*log(x))/(24.*
799 pow7) - (67.*x*x*x*x*log(x))/(2304.*
800 pow6) - (7.*x*x*x*x*log(x))/(5184.*
801 pow5) - (19.*x*x*x*x*x*log(x))/(5184.*
802 pow7) + (29.*x*x*x*x*x*log(x))/(3456.*
803 pow6) + (x*x*x*x*x*x*log(x))/(5184.*
804 pow7) - (x*x*x*x*x*x*log(x))/(1296.*
805 pow6) + (17.*x*a*log(x))/(1728.*
806 pow5) + (79.*x*x*a*
807 log(x))/(1728.*pow5) - (71.*x*
808 pow(log(x), 2.))/(1536.*pow6) - (89.*x*
809 pow(log(x), 2.))/(2304.*pow5) - (x*
810 pow(log(x), 2.))/(288.*pow4) - (1007.*
811 x*x*pow(log(x), 2.))/(2304.*
812 pow6) - (281.*x*x*
813 pow(log(x), 2.))/(2304.*pow5) + (11.*
814 x*x*x*pow(log(x), 2.))/(864.*
815 pow7) - (2551.*x*x*x*
816 pow(log(x), 2.))/(13824.*pow6) + (13.*
817 x*x*x*x*pow(log(x), 2.))/(864.*
818 pow7) +
819 log(Ms2/mumatch2)*((5.*Nf)/(3456.*pow6) -
820 95./(6912.*pow5) - (1927.*x)/(6912.*
821 pow6) + (5.*Nf*x)/(108.*
822 pow6) - (2033.*x)/(13824.*
823 pow5) - (1211.*x*x)/(10368.*
824 pow6) + (Nf*x*x)/(192.*
825 pow6) + (209.*x*x)/(1536.*
826 pow5) + (541.*x*x*x)/(1296.*
827 pow6) - (11.*Nf*x*x*x)/(216.*
828 pow6) + (361.*x*x*x)/(13824.*
829 pow5) - (103.*x*x*x*x)/(3456.*
830 pow6) - (7.*Nf*x*x*x*x)/(3456.*
831 pow6) - (19.*x*x*x*x)/(13824.*
832 pow5) + (181.*x*x*x*x*x)/(20736.*
833 pow6) -
834 x*x*x*x*x*x/(1296.*pow6) - (11.*x*
835 log(x))/(128.*pow6) + (11.*Nf*x*
836 log(x))/(576.*pow6) - (209.*x*
837 log(x))/(2304.*pow5) - (769.*
838 x*x*log(x))/(1728.*
839 pow6) + (Nf*x*x*log(x))/(16.*
840 pow6) - (247.*x*x*
841 log(x))/(2304.*pow5) - (445.*
842 x*x*x*log(x))/(3456.*
843 pow6) + (13.*Nf*x*x*x*
844 log(x))/(576.*pow6)))
845
846 );
847}
848
849 std::vector<WilsonCoefficient>& SUSYMassInsertionMatching::CMd1() {
850
851 vmcd1.clear();
852 vmcd1 = StandardModelMatching::CMd1();
853
854 double x = pow(SusyMI.getM3() / SusyMI.getMsq(), 2.);
855 double constLO = SusyMI.Als(SusyMI.getMuM()) / SusyMI.getMsq();
856 double mcharm = SusyMI.getMuc();
857 double mgluino = SusyMI.getM3();
858
859 DLL = SusyMI.getDu_LL()(0,1);
860 DLR = SusyMI.getDu_LR()(0,1);
861 DRL = SusyMI.getDu_RL()(0,1);
862 DRR = SusyMI.getDu_RR()(0,1);
863 /*
864 std::cout << " DLL --> " << DLL << std::endl;
865 std::cout << " DLR --> " << DLR << std::endl;
866 std::cout << " DRL --> " << DRL << std::endl;
867 std::cout << " DRR --> " << DRR << std::endl;
868 */
869
870 switch (mcd1.getScheme()) {
871 case NDR:
872 //case HV:
873 //case LRI:
874 break;
875 default:
876 std::stringstream out;
877 out << mcd1.getScheme();
878 throw std::runtime_error("SUSYMassInsertionMatching::CMD1(): scheme " + out.str() + "not implemented");
879 }
880
881 mcd1.setMu(SusyMI.getMuM());
882
883 gslpp::vector<gslpp::complex> C0_B(10,0.), C0_M(10,0.), C0_M_eff(10,0.);
884
885 C0_B.assign(2, constLO * constLO / 4. * (-1./9.*B1(x) - 5./9.*B2(x)
886 - 1./18.*P1(x) - 1./2.*P2(x))* DLL);
887 C0_B.assign(3, constLO * constLO / 4. * (-7./3.*B1(x) + 1./3.*B2(x)
888 + 1./6.*P1(x) + 3./2.*P2(x))* DLL);
889 C0_B.assign(4, constLO * constLO / 4. * (10./9.*B1(x) + 1./18.*B2(x)
890 - 1./18.*P1(x) - 1./2.*P2(x))* DLL);
891 C0_B.assign(5, constLO * constLO / 4. * (-2./3.*B1(x) + 7./6.*B2(x)
892 + 1./6.*P1(x) + 3./2.*P2(x)) * DLL);
893 C0_B.assign(6, constLO/SusyMI.getMsq() / 4. * M_PI * ((8./3.*M3(x)) * DLL
894 + mgluino/mcharm * (8./3.*M1(x)) * DLR));
895 C0_B.assign(7, constLO/SusyMI.getMsq() / 4. * M_PI * ((-1./3.*M3(x)-3.*M4(x)) * DLL
896 + mgluino/mcharm * (-1./3.*M1(x)-3.*M2(x)) * DLR));
897
898 C0_M = RtoMisiak().transpose() * C0_B;
899 C0_M_eff = EffectiveBase() * C0_M ;
900
901 /* for(int i= 0; i< 10; i ++){
902 std::cout << "i" << " "<< "C0_B(i)" << " " << "C0_M" << " " << "C0_M_eff" << std::endl;
903 std::cout << i << " "<< C0_B(i) << " " << C0_M(i) << " " << C0_M_eff(i) << std::endl;
904 } */
905
906 switch (mcd1.getOrder()) {
907 case NLO:
908 for(int k = 0; k<10; k++){
909 mcd1.setCoeff(k, 0., NLO);
910 }
911 case LO:
912 for(int k = 0; k<10; k++){
913 mcd1.setCoeff(k, C0_M_eff(k), LO);
914 }
915 break;
916 default:
917 std::stringstream out;
918 out << mcd1.getOrder();
919 throw std::runtime_error("SUSYMassInsertionMatching::CMd1(): order " + out.str() + "not implemented");
920 }
921
922 vmcd1.push_back(mcd1);
923
924 C0_B.assign(2, constLO * constLO / 4. * (-1./9.*B1(x) - 5./9.*B2(x)
925 - 1./18.*P1(x) - 1./2.*P2(x))* DRR);
926 C0_B.assign(3, constLO * constLO / 4. * (-7./3.*B1(x) + 1./3.*B2(x)
927 + 1./6.*P1(x) + 3./2.*P2(x))* DRR);
928 C0_B.assign(4, constLO * constLO / 4. * (10./9.*B1(x) + 1./18.*B2(x)
929 - 1./18.*P1(x) - 1./2.*P2(x))* DRR);
930 C0_B.assign(5, constLO * constLO / 4. * (-2./3.*B1(x) + 7./6.*B2(x)
931 + 1./6.*P1(x) + 3./2.*P2(x)) * DRR);
932 C0_B.assign(6, constLO/SusyMI.getMsq() / 4. * M_PI * ((8./3.*M3(x)) * DRR
933 + mgluino/mcharm * (8./3.*M1(x) * DRL)));
934 C0_B.assign(7, (constLO/SusyMI.getMsq() / 4. * M_PI * ((-1./3.*M3(x)-3.*M4(x)) * DRR
935 + mgluino/mcharm * (-1./3.*M1(x)-3.*M2(x))* DRL)));
936 C0_M = RtoMisiak().transpose() * C0_B;
937 C0_M_eff = EffectiveBase() * C0_M ;
938
939 /* for(int i= 0; i< 10; i ++){
940 std::cout << "i" << " "<< "C0_B(i)" << " " << "C0_M" << " " << "C0_M_eff" << std::endl;
941 std::cout << i << " "<< C0_B(i) << " " << C0_M(i) << " " << C0_M_eff(i) << std::endl;
942 } */
943
944 switch (mcd1.getOrder()) {
945 case NLO:
946 for(int k = 0; k<10; k++){
947 mcd1.setCoeff(k, 0., NLO);
948 }
949 case LO:
950 for(int k=0; k<10; k++){
951 mcd1.setCoeff(k, C0_M_eff(k), LO);
952 }
953 break;
954 default:
955 std::stringstream out;
956 out << mcd1.getOrder();
957 throw std::runtime_error("SUSYMassInsertionMatching::CMd1(): order " + out.str() + "not implemented");
958 }
959
960 vmcd1.push_back(mcd1);
961
962 return(vmcd1);
963}
964
965 std::vector<WilsonCoefficient>& SUSYMassInsertionMatching::CMdd2 () {
966
967 vmcd2.clear();
968 vmcd2 = StandardModelMatching::CMdd2();
969
970 double als = SusyMI.Als(SusyMI.getMuM())*(1.+SusyMI.Als(SusyMI.getM3())/4./M_PI);
971 double x = pow(SusyMI.getM3() / SusyMI.getMsq(), 2.);
972 double coLO = pow(als/ SusyMI.getMsq(), 2.);
973 double coNLO = coLO* als / M_PI;
974 MuM2 = SusyMI.getMuM()*SusyMI.getMuM();
975 Ms2 = SusyMI.getMsq()*SusyMI.getMsq();
976 DLL = SusyMI.getDu_LL()(0,1);
977 DLR = SusyMI.getDu_LR()(0,1);
978 DRL = SusyMI.getDu_RL()(0,1);
979 DRR = SusyMI.getDu_RR()(0,1);
980
981 switch (mcd2.getScheme()) {
982 case NDR:
983 //case HV:
984 //case LRI:
985 break;
986 default:
987 std::stringstream out;
988 out << mcd2.getScheme();
989 throw std::runtime_error("StandardModel::CMbd2(): scheme " + out.str() + "not implemented");
990 }
991 mcd2.setMu(SusyMI.getMuM());
992 switch (mcd2.getOrder()) {
993 case NLO:
994 mcd2.setCoeff(0, coNLO * C0NLO(x,MuM2,Ms2) * DLL * DLL, NLO);
995 mcd2.setCoeff(1, coNLO * C1NLO(x,MuM2,Ms2) * DRL * DRL, NLO);
996 mcd2.setCoeff(2, coNLO * C2NLO(x,MuM2,Ms2) * DRL * DRL, NLO);
997 mcd2.setCoeff(3, coNLO * (C3NLOA(x,MuM2,Ms2) * DLR * DRL +
998 C3NLOB(x,MuM2,Ms2) * DLL * DRR), NLO);
999 mcd2.setCoeff(4, coNLO * (C4NLOA(x,MuM2,Ms2) * DLR * DRL +
1000 C4NLOB(x,MuM2,Ms2) * DLL * DRR), NLO);
1001 case LO:
1002 mcd2.setCoeff(0, coLO * C0LO(x) * DLL * DLL, LO);
1003 mcd2.setCoeff(1, coLO * C1LO(x) * DRL * DRL, LO);
1004 mcd2.setCoeff(2, coLO * C2LO(x) * DRL * DRL, LO);
1005 mcd2.setCoeff(3, coLO * (C3LOA(x) * DLR * DRL + C3LOB(x) * DLL * DRR), LO);
1006 mcd2.setCoeff(4, coLO * (C4LOA(x) * DLR * DRL + C4LOB(x) * DLL * DRR), LO);
1007 break;
1008 default:
1009 std::stringstream out;
1010 out << mcd2.getOrder();
1011 throw std::runtime_error("StandardModelMatching::CMd2(): order " + out.str() + "not implemented");
1012 }
1013
1014 LRItoNDR(1);
1015
1016 vmcd2.push_back(mcd2);
1017
1018 switch (mcd2.getOrder()) {
1019 case NLO:
1020 mcd2.setCoeff(0, coNLO * C0NLO(x,MuM2,Ms2) * DRR * DRR, NLO);
1021 mcd2.setCoeff(1, coNLO * C1NLO(x,MuM2,Ms2) * DLR * DLR, NLO);
1022 mcd2.setCoeff(2, coNLO * C2NLO(x,MuM2,Ms2) * DLR * DLR, NLO);
1023 mcd2.setCoeff(3, 0., NLO);
1024 mcd2.setCoeff(4, 0., NLO);
1025 case LO:
1026 mcd2.setCoeff(0, coLO * C0LO(x) * DRR * DRR, LO);
1027 mcd2.setCoeff(1, coLO * C1LO(x) * DLR * DLR, LO);
1028 mcd2.setCoeff(2, coLO * C2LO(x) * DLR * DLR, LO);
1029 mcd2.setCoeff(3, 0., LO);
1030 mcd2.setCoeff(4, 0., LO);
1031 break;
1032 default:
1033 std::stringstream out;
1034 out << mcd2.getOrder();
1035 throw std::runtime_error("StandardModelMatching::CMd2(): order " + out.str() + "not implemented");
1036 }
1037
1038 LRItoNDR(1);
1039
1040 vmcd2.push_back(mcd2);
1041
1042 return(vmcd2);
1043}
1044
1045 std::vector<WilsonCoefficient>& SUSYMassInsertionMatching::CMdbd2() {
1046
1047 vmcdb.clear();
1048
1049 vmcdb = StandardModelMatching::CMdbd2();
1050
1051 double als = SusyMI.Als(SusyMI.getMuM())*(1.+SusyMI.Als(SusyMI.getM3())/4./M_PI);
1052 double x = pow(SusyMI.getM3() / SusyMI.getMsq(), 2.);
1053 double coLO = pow(als/ SusyMI.getMsq(), 2.);
1054 double coNLO = coLO * als / M_PI;
1055 MuM2 = SusyMI.getMuM()*SusyMI.getMuM();
1056 Ms2 = SusyMI.getMsq()*SusyMI.getMsq();
1057
1058 DLL = SusyMI.getDd_LL()(0,2);
1059 DLR = SusyMI.getDd_LR()(0,2);
1060 DRL = SusyMI.getDd_RL()(0,2);
1061 DRR = SusyMI.getDd_RR()(0,2);
1062
1063 switch (mcbd.getScheme()) {
1064 case NDR:
1065 //case HV:
1066 //case LRI:
1067 break;
1068 default:
1069 std::stringstream out;
1070 out << mcbd.getScheme();
1071 throw std::runtime_error("StandardModel::CMbbd(): scheme " + out.str() + "not implemented");
1072 }
1073
1074 mcbd.setMu(SusyMI.getMuM());
1075
1076 switch (mcbd.getOrder()) {
1077 case NLO:
1078 mcbd.setCoeff(0, coNLO * C0NLO(x,MuM2,Ms2) * DLL * DLL, NLO);
1079 mcbd.setCoeff(1, coNLO * C1NLO(x,MuM2,Ms2) * DRL * DRL, NLO);
1080 mcbd.setCoeff(2, coNLO * C2NLO(x,MuM2,Ms2) * DRL * DRL, NLO);
1081 mcbd.setCoeff(3, coNLO * (C3NLOA(x,MuM2,Ms2) * DLR * DRL +
1082 C3NLOB(x,MuM2,Ms2) * DLL * DRR), NLO);
1083 mcbd.setCoeff(4, coNLO * (C4NLOA(x,MuM2,Ms2) * DLR * DRL +
1084 C4NLOB(x,MuM2,Ms2) * DLL * DRR), NLO);
1085 case LO:
1086 mcbd.setCoeff(0, coLO * C0LO(x) * DLL * DLL, LO);
1087 mcbd.setCoeff(1, coLO * C1LO(x) * DRL * DRL, LO);
1088 mcbd.setCoeff(2, coLO * C2LO(x) * DRL * DRL, LO);
1089 mcbd.setCoeff(3, coLO * (C3LOA(x) * DLR * DRL + C3LOB(x) * DLL * DRR), LO);
1090 mcbd.setCoeff(4, coLO * (C4LOA(x) * DLR * DRL + C4LOB(x) * DLL * DRR), LO);
1091 break;
1092 default:
1093 std::stringstream out;
1094 out << mcbd.getOrder();
1095 throw std::runtime_error("StandardModelMatching::CMbd(): order " + out.str() + "not implemented");
1096 }
1097
1098 LRItoNDR(2);
1099
1100 vmcdb.push_back(mcbd);
1101
1102 switch (mcbd.getOrder()) {
1103 case NLO:
1104 mcbd.setCoeff(0, coNLO * C0NLO(x,MuM2,Ms2) * DRR * DRR, NLO);
1105 mcbd.setCoeff(1, coNLO * C1NLO(x,MuM2,Ms2) * DLR * DLR, NLO);
1106 mcbd.setCoeff(2, coNLO * C2NLO(x,MuM2,Ms2) * DLR * DLR, NLO);
1107 mcbd.setCoeff(3, 0., NLO);
1108 mcbd.setCoeff(4, 0., NLO);
1109 case LO:
1110 mcbd.setCoeff(0, coLO * C0LO(x) * DRR * DRR, LO);
1111 mcbd.setCoeff(1, coLO * C1LO(x) * DLR * DLR, LO);
1112 mcbd.setCoeff(2, coLO * C2LO(x) * DLR * DLR, LO);
1113 mcbd.setCoeff(3, 0., LO);
1114 mcbd.setCoeff(4, 0., LO);
1115 break;
1116 default:
1117 std::stringstream out;
1118 out << mcbd.getOrder();
1119 throw std::runtime_error("StandardModelMatching::CMbd(): order " + out.str() + "not implemented");
1120 }
1121
1122 LRItoNDR(2);
1123
1124 vmcdb.push_back(mcbd);
1125
1126 return(vmcdb);
1127}
1128
1129 std::vector<WilsonCoefficient>& SUSYMassInsertionMatching::CMdbs2() {
1130
1131
1132 vmcds.clear();
1133
1134 vmcds = StandardModelMatching::CMdbs2();
1135
1136 double als = SusyMI.Als(SusyMI.getMuM())*(1.+SusyMI.Als(SusyMI.getM3())/4./M_PI);
1137 double x = pow(SusyMI.getM3() / SusyMI.getMsq(), 2.);
1138 double coLO = pow(als / SusyMI.getMsq(), 2.);
1139 double coNLO = coLO * als / M_PI;
1140 MuM2 = SusyMI.getMuM()*SusyMI.getMuM();
1141 Ms2 = SusyMI.getMsq()*SusyMI.getMsq();
1142
1143 DLL = SusyMI.getDd_LL()(1,2);
1144 DLR = SusyMI.getDd_LR()(1,2);
1145 DRL = SusyMI.getDd_RL()(1,2);
1146 DRR = SusyMI.getDd_RR()(1,2);
1147
1148 switch (mcbs.getScheme()) {
1149 case NDR:
1150 //case HV:
1151 //case LRI:
1152 break;
1153 default:
1154 std::stringstream out;
1155 out << mcbs.getScheme();
1156 throw std::runtime_error("StandardModel::CMbbs(): scheme " + out.str() + "not implemented");
1157 }
1158
1159 mcbs.setMu(SusyMI.getMuM());
1160
1161 switch (mcbs.getOrder()) {
1162 case NLO:
1163 mcbs.setCoeff(0, coNLO * C0NLO(x,MuM2,Ms2) * DLL * DLL, NLO);
1164 mcbs.setCoeff(1, coNLO * C1NLO(x,MuM2,Ms2) * DRL * DRL, NLO);
1165 mcbs.setCoeff(2, coNLO * C2NLO(x,MuM2,Ms2) * DRL * DRL, NLO);
1166 mcbs.setCoeff(3, coNLO * (C3NLOA(x,MuM2,Ms2) * DLR * DRL
1167 + C3NLOB(x,MuM2,Ms2) * DLL * DRR), NLO);
1168 mcbs.setCoeff(4, coNLO * (C4NLOA(x,MuM2,Ms2) * DLR * DRL +
1169 C4NLOB(x,MuM2,Ms2) * DLL * DRR), NLO);
1170 case LO:
1171 mcbs.setCoeff(0, coLO * C0LO(x) * DLL * DLL, LO);
1172 mcbs.setCoeff(1, coLO * C1LO(x) * DRL * DRL, LO);
1173 mcbs.setCoeff(2, coLO * C2LO(x) * DRL * DRL, LO);
1174 mcbs.setCoeff(3, coLO * (C3LOA(x) * DLR * DRL + C3LOB(x) * DLL * DRR), LO);
1175 mcbs.setCoeff(4, coLO * (C4LOA(x) * DLR * DRL + C4LOB(x) * DLL * DRR), LO);
1176 break;
1177 default:
1178 std::stringstream out;
1179 out << mcbs.getOrder();
1180 throw std::runtime_error("StandardModelMatching::CMbs(): order " + out.str() + "not implemented");
1181 }
1182
1183 LRItoNDR(3);
1184
1185 vmcds.push_back(mcbs);
1186
1187 switch (mcbs.getOrder()) {
1188 case NLO:
1189 mcbs.setCoeff(0, coNLO * C0NLO(x,MuM2,Ms2) * DRR * DRR, NLO);
1190 mcbs.setCoeff(1, coNLO * C1NLO(x,MuM2,Ms2) * DLR * DLR, NLO);
1191 mcbs.setCoeff(2, coNLO * C2NLO(x,MuM2,Ms2) * DLR * DRL, NLO);
1192 mcbs.setCoeff(3, 0., NLO);
1193 mcbs.setCoeff(4, 0., NLO);
1194 case LO:
1195 mcbs.setCoeff(0, coLO * C0LO(x) * DRR * DRR, LO);
1196 mcbs.setCoeff(1, coLO * C1LO(x) * DLR * DLR, LO);
1197 mcbs.setCoeff(2, coLO * C2LO(x) * DLR * DLR, LO);
1198 mcbs.setCoeff(3, 0., LO);
1199 mcbs.setCoeff(4, 0., LO);
1200 break;
1201 default:
1202 std::stringstream out;
1203 out << mcbs.getOrder();
1204 throw std::runtime_error("StandardModelMatching::CMbs(): order " + out.str() + "not implemented");
1205 }
1206
1207 LRItoNDR(3);
1208
1209 vmcds.push_back(mcbs);
1210
1211 return(vmcds);
1212}
1213
1214 std::vector<WilsonCoefficient>& SUSYMassInsertionMatching::CMdk2() {
1215
1216 vmck2.clear();
1217 vmck2 = StandardModelMatching::CMdk2();
1218
1219 double als = SusyMI.Als(SusyMI.getMuM())*(1.+SusyMI.Als(SusyMI.getM3())/4./M_PI);
1220 double x = pow(SusyMI.getM3() / SusyMI.getMsq(), 2.);
1221 double coLO = pow(als/ SusyMI.getMsq(), 2.);
1222 double coNLO = coLO * als / M_PI;
1223 MuM2 = SusyMI.getMuM()*SusyMI.getMuM();
1224 Ms2 = SusyMI.getMsq()*SusyMI.getMsq();
1225
1226 DLL = SusyMI.getDd_LL()(0,1);
1227 DLR = SusyMI.getDd_LR()(0,1);
1228 DRL = SusyMI.getDd_RL()(0,1);
1229 DRR = SusyMI.getDd_RR()(0,1);
1230
1231 switch (mck2.getScheme()) {
1232 case NDR:
1233 //case HV:
1234 //case LRI:
1235 break;
1236 default:
1237 std::stringstream out;
1238 out << mck2.getScheme();
1239 throw std::runtime_error("StandardModel::CMbk2(): scheme " + out.str() + "not implemented");
1240 }
1241
1242 mck2.setMu(SusyMI.getMuM());
1243
1244 switch (mck2.getOrder()) {
1245 case NLO:
1246 mck2.setCoeff(0, coNLO * C0NLO(x,MuM2,Ms2) * DLL * DLL, NLO);
1247 mck2.setCoeff(1, coNLO * C1NLO(x,MuM2,Ms2) * DRL * DRL, NLO);
1248 mck2.setCoeff(2, coNLO * C2NLO(x,MuM2,Ms2) * DRL * DRL, NLO);
1249 mck2.setCoeff(3, coNLO * C3NLOA(x,MuM2,Ms2) * DLR * DRL +
1250 coNLO * C3NLOB(x,MuM2,Ms2) * DLL * DRR, NLO);
1251 mck2.setCoeff(4, coNLO * C4NLOA(x,MuM2,Ms2) * DLR * DRL +
1252 coNLO * C4NLOB(x,MuM2,Ms2) * DLL * DRR, NLO);
1253 case LO:
1254 mck2.setCoeff(0, coLO * C0LO(x) * DLL * DLL, LO);
1255 mck2.setCoeff(1, coLO * C1LO(x) * DRL * DRL, LO);
1256 mck2.setCoeff(2, coLO * C2LO(x) * DRL * DRL, LO);
1257 mck2.setCoeff(3, coLO * C3LOA(x) * DLR * DRL + coLO * coLO * C3LOB(x) * DLL * DRR, LO);
1258 mck2.setCoeff(4, coLO * C4LOA(x) * DLR * DRL + coLO * coLO * C4LOB(x) * DLL * DRR, LO);
1259 break;
1260 default:
1261 std::stringstream out;
1262 out << mck2.getOrder();
1263 throw std::runtime_error("StandardModelMatching::CMk2(): order " + out.str() + "not implemented");
1264 }
1265
1266 LRItoNDR(4);
1267
1268 vmck2.push_back(mck2);
1269
1270 switch (mck2.getOrder()) {
1271 case NLO:
1272 mck2.setCoeff(0, coNLO * C0NLO(x,MuM2,Ms2) * DRR * DRR, NLO);
1273 mck2.setCoeff(1, coNLO * C1NLO(x,MuM2,Ms2) * DLR * DLR, NLO);
1274 mck2.setCoeff(2, coNLO * C2NLO(x,MuM2,Ms2) * DLR * DLR, NLO);
1275 mck2.setCoeff(3, 0., NLO);
1276 mck2.setCoeff(4, 0., NLO);
1277 case LO:
1278 mck2.setCoeff(0, coLO * C0LO(x) * DRR * DRR, LO);
1279 mck2.setCoeff(1, coLO * C1LO(x) * DLR * DLR, LO);
1280 mck2.setCoeff(2, coLO * C2LO(x) * DLR * DLR, LO);
1281 mck2.setCoeff(3, 0., LO);
1282 mck2.setCoeff(4, 0., LO);
1283 break;
1284 default:
1285 std::stringstream out;
1286 out << mck2.getOrder();
1287 throw std::runtime_error("StandardModelMatching::CMk2(): order " + out.str() + "not implemented");
1288 }
1289
1290 LRItoNDR(4);
1291
1292 vmck2.push_back(mck2);
1293
1294 return(vmck2);
1295}
1296
1297gslpp::matrix<double> SUSYMassInsertionMatching::RtoMisiak() const{
1298
1299 gslpp::matrix<double> R(10,0.);
1300
1301 R(0,0) = 2.;
1302 R(0,1) = 1./3.;
1303 R(1,1) = 1.;
1304 R(2,2) = -1./3.;
1305 R(2,4) = 1./12.;
1306 R(3,2) = -1./9.;
1307 R(3,3) = -2./3.;
1308 R(3,4) = 1./36.;
1309 R(3,5) = 1./6.;
1310 R(4,2) = 4./3.;
1311 R(4,4) = -1./12.;
1312 R(5,2) = 4./9.;
1313 R(5,3) = 8./3.;
1314 R(5,4) = -1./36.;
1315 R(5,5) = -1./6.;
1316 R(6,6) = 1.;
1317 R(7,7) = 1.;
1318 R(8,8) = 1.;
1319 R(9,9) = 1.;
1320 return(R);
1321
1322}
1323
1324gslpp::matrix<double> SUSYMassInsertionMatching::EffectiveBase() const{
1325
1326 gslpp::matrix<double> y(10, 0.);
1327
1328 y(0,0) = 1.;
1329 y(1,1) = 1.;
1330 y(2,2) = 1.;
1331 y(3,3) = 1.;
1332 y(4,4) = 1.;
1333 y(5,5) = 1.;
1334 y(6,6) = 1.;
1335 y(7,7) = 1.;
1336 y(8,8) = 1.;
1337 y(9,9) = 1.;
1338
1339 y(6,2) = -1./3.;
1340 y(6,3) = -4./9.;
1341 y(6,4) = -20./3.;
1342 y(6,5) = -80./9.;
1343
1344 y(7,2) = 1.;
1345 y(7,3) = -1./6.;
1346 y(7,4) = 20.;
1347 y(7,5) = -10./3.;
1348
1349 return(y);
1350
1351}
1352
1353void SUSYMassInsertionMatching::LRItoNDR(int i){
1354 switch (i){
1355 case 1:
1356 mcd2.setCoeff(*mcd2.getCoeff(NLO) + SusyMI.Als(mcd2.getMu()) / 4. / M_PI *
1357 drNDRLRI.transpose() * (*mcd2.getCoeff(LO)), NLO);
1358 break;
1359 case 2:
1360 mcbd.setCoeff(*mcbd.getCoeff(NLO) + SusyMI.Als(mcbd.getMu()) / 4. / M_PI *
1361 drNDRLRI.transpose() * (*mcbd.getCoeff(LO)), NLO);
1362 break;
1363 case 3:
1364 mcbs.setCoeff(*mcbs.getCoeff(NLO) + SusyMI.Als(mcbs.getMu()) / 4. / M_PI *
1365 drNDRLRI.transpose() * (*mcbs.getCoeff(LO)), NLO);
1366 break;
1367 case 4:
1368 mck2.setCoeff(*mck2.getCoeff(NLO) + SusyMI.Als(mck2.getMu()) / 4. / M_PI *
1369 drNDRLRI.transpose() * (*mck2.getCoeff(LO)), NLO);
1370 break;
1371 default:
1372 throw std::runtime_error("SUSYMassInsertionMatching::LRItoNDR : change of scheme not implemented");
1373 }
1374}
LO
@ LO
Definition: OrderScheme.h:34
NLO
@ NLO
Definition: OrderScheme.h:35
NDR
@ NDR
Definition: OrderScheme.h:21
ModelMatching::CMd1
virtual std::vector< WilsonCoefficient > & CMd1()=0
QCD::getMuc
const double getMuc() const
A get method to access the threshold between four- and three-flavour theory in GeV.
Definition: QCD.h:591
QCD::getNc
const double getNc() const
A get method to access the number of colours .
Definition: QCD.h:507
SUSYMassInsertion::getMsq
double getMsq() const
Definition: SUSYMassInsertion.h:93
SUSYMassInsertion::getDu_LL
gslpp::matrix< gslpp::complex > getDu_LL() const
Definition: SUSYMassInsertion.h:109
SUSYMassInsertion::getDu_RR
gslpp::matrix< gslpp::complex > getDu_RR() const
Definition: SUSYMassInsertion.h:157
SUSYMassInsertion::getDd_RL
gslpp::matrix< gslpp::complex > getDd_RL() const
Definition: SUSYMassInsertion.h:205
SUSYMassInsertion::getMuM
double getMuM() const
Definition: SUSYMassInsertion.h:236
SUSYMassInsertion::getM3
double getM3() const
Definition: SUSYMassInsertion.h:77
SUSYMassInsertion::getDd_LL
gslpp::matrix< gslpp::complex > getDd_LL() const
Definition: SUSYMassInsertion.h:173
SUSYMassInsertion::getDu_RL
gslpp::matrix< gslpp::complex > getDu_RL() const
Definition: SUSYMassInsertion.h:141
SUSYMassInsertion::getDu_LR
gslpp::matrix< gslpp::complex > getDu_LR() const
Definition: SUSYMassInsertion.h:125
SUSYMassInsertion::getDd_LR
gslpp::matrix< gslpp::complex > getDd_LR() const
Definition: SUSYMassInsertion.h:189
SUSYMassInsertion::getDd_RR
gslpp::matrix< gslpp::complex > getDd_RR() const
Definition: SUSYMassInsertion.h:221
SUSYMassInsertionMatching::C0NLO
double C0NLO(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
Definition: SUSYMassInsertionMatching.cpp:138
SUSYMassInsertionMatching::B2
double B2(double x) const
loop functions, LO term
Definition: SUSYMassInsertionMatching.cpp:43
SUSYMassInsertionMatching::P1
double P1(double x) const
loop functions, LO term
Definition: SUSYMassInsertionMatching.cpp:47
SUSYMassInsertionMatching::M2
double M2(double x) const
loop functions, LO term
Definition: SUSYMassInsertionMatching.cpp:61
SUSYMassInsertionMatching::M4
double M4(double x) const
loop functions, LO term
Definition: SUSYMassInsertionMatching.cpp:69
SUSYMassInsertionMatching::C2LO
double C2LO(double x)
loop functions RI-MOM scheme, LO term
Definition: SUSYMassInsertionMatching.cpp:106
SUSYMassInsertionMatching::C1LO
double C1LO(double x)
loop functions RI-MOM scheme, LO term
Definition: SUSYMassInsertionMatching.cpp:101
SUSYMassInsertionMatching::CMdbd2
virtual std::vector< WilsonCoefficient > & CMdbd2()
,
Definition: SUSYMassInsertionMatching.cpp:1045
SUSYMassInsertionMatching::CMd1
virtual std::vector< WilsonCoefficient > & CMd1()
Definition: SUSYMassInsertionMatching.cpp:849
SUSYMassInsertionMatching::C3LOB
double C3LOB(double x)
loop functions RI-MOM scheme, LO term
Definition: SUSYMassInsertionMatching.cpp:116
SUSYMassInsertionMatching::C4LOB
double C4LOB(double x)
loop functions RI-MOM scheme, LO term
Definition: SUSYMassInsertionMatching.cpp:127
SUSYMassInsertionMatching::C3LOA
double C3LOA(double x)
loop functions RI-MOM scheme, LO term
Definition: SUSYMassInsertionMatching.cpp:111
SUSYMassInsertionMatching::CMdbs2
virtual std::vector< WilsonCoefficient > & CMdbs2()
,
Definition: SUSYMassInsertionMatching.cpp:1129
SUSYMassInsertionMatching::P2
double P2(double x) const
loop functions, LO term
Definition: SUSYMassInsertionMatching.cpp:52
SUSYMassInsertionMatching::C0LO
double C0LO(double x)
loop functions RI-MOM scheme, LO term
Definition: SUSYMassInsertionMatching.cpp:91
SUSYMassInsertionMatching::C3NLOA
double C3NLOA(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
Definition: SUSYMassInsertionMatching.cpp:433
SUSYMassInsertionMatching::C2NLO
double C2NLO(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
Definition: SUSYMassInsertionMatching.cpp:341
SUSYMassInsertionMatching::C1NLO
double C1NLO(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
Definition: SUSYMassInsertionMatching.cpp:248
SUSYMassInsertionMatching::M3
double M3(double x) const
loop functions, LO term
Definition: SUSYMassInsertionMatching.cpp:65
SUSYMassInsertionMatching::B1
double B1(double x) const
loop functions, LO term
Definition: SUSYMassInsertionMatching.cpp:39
SUSYMassInsertionMatching::C4NLOA
double C4NLOA(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
Definition: SUSYMassInsertionMatching.cpp:635
SUSYMassInsertionMatching::C3NLOB
double C3NLOB(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
Definition: SUSYMassInsertionMatching.cpp:534
SUSYMassInsertionMatching::SUSYMassInsertionMatching
SUSYMassInsertionMatching(const SUSYMassInsertion &SUSYMassInsertion_i)
constructor
Definition: SUSYMassInsertionMatching.cpp:11
SUSYMassInsertionMatching::EffectiveBase
gslpp::matrix< double > EffectiveBase() const
the basis is defined, for example, in Chetyrkin et al hep-ph/9612313
Definition: SUSYMassInsertionMatching.cpp:1324
SUSYMassInsertionMatching::drNDRLRI
gslpp::matrix< double > drNDRLRI
Definition: SUSYMassInsertionMatching.h:278
SUSYMassInsertionMatching::LRItoNDR
void LRItoNDR(int i)
it changes renormalization scheme from LRI to NDR
Definition: SUSYMassInsertionMatching.cpp:1353
SUSYMassInsertionMatching::CMdk2
virtual std::vector< WilsonCoefficient > & CMdk2()
Definition: SUSYMassInsertionMatching.cpp:1214
SUSYMassInsertionMatching::C4LOA
double C4LOA(double x)
loop functions RI-MOM scheme, LO term
Definition: SUSYMassInsertionMatching.cpp:121
SUSYMassInsertionMatching::CMdd2
virtual std::vector< WilsonCoefficient > & CMdd2()
Definition: SUSYMassInsertionMatching.cpp:965
SUSYMassInsertionMatching::M1
double M1(double x) const
loop functions, LO term
Definition: SUSYMassInsertionMatching.cpp:57
SUSYMassInsertionMatching::SusyMI
const SUSYMassInsertion & SusyMI
Definition: SUSYMassInsertionMatching.h:274
SUSYMassInsertionMatching::C4NLOB
double C4NLOB(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
Definition: SUSYMassInsertionMatching.cpp:734
SUSYMassInsertionMatching::RtoMisiak
gslpp::matrix< double > RtoMisiak() const
Misiak et al, hep-ph/0005183.
Definition: SUSYMassInsertionMatching.cpp:1297
StandardModel::Als
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.
Definition: StandardModel.h:1088
StandardModelMatching
A class for the matching in the Standard Model.
Definition: StandardModelMatching.h:26
StandardModelMatching::CMdbs2
virtual std::vector< WilsonCoefficient > & CMdbs2()
,
Definition: StandardModelMatching.cpp:1016
StandardModelMatching::CMdbd2
virtual std::vector< WilsonCoefficient > & CMdbd2()
,
Definition: StandardModelMatching.cpp:966
StandardModelMatching::CMdd2
virtual std::vector< WilsonCoefficient > & CMdd2()
,
Definition: StandardModelMatching.cpp:1169
WilsonCoefficient::getCoeff
gslpp::vector< gslpp::complex > ** getCoeff() const
Definition: WilsonCoefficient.h:29
WilsonCoefficient::setCoeff
void setCoeff(const gslpp::vector< gslpp::complex > &z, orders order_i)
Definition: WilsonCoefficient.h:34
WilsonTemplate::getMu
double getMu() const
Definition: WilsonTemplate.h:75
WilsonTemplate::getScheme
schemes getScheme() const
Definition: WilsonTemplate.h:98
WilsonTemplate::setMu
virtual void setMu(double mu)
Definition: WilsonTemplate.h:92
WilsonTemplate::getOrder
orders getOrder() const
Definition: WilsonTemplate.h:65
Li2
cd Li2(cd x)
Definition: hpl.h:1011