a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
Public Member Functions | Private Attributes | List of all members
SUSYMassInsertionMatching Class Reference

A class for the matching of SUSY MIA. More...

#include <SUSYMassInsertionMatching.h>

+ Inheritance diagram for SUSYMassInsertionMatching:

Detailed Description

A class for the matching of SUSY MIA.

Author
HEPfit Collaboration
Copyright
GNU General Public License

this class contains the Wilson Coefficients at the SUSY matching scale for \( \Delta F = 1 \) and \( \Delta F = 2 \) processes in the mass insertion approximation; references: \Delta F = 2 \( Wilson coefficients: M. Ciuchini et al, hep-ph/0606197 \Delta F = 1 \) Wilson coefficients: F.Gabbiani et al, hep-ph/9604387

Parameters
SUSYMassInsertion_iobject of type SUSYMassInsertion

Definition at line 31 of file SUSYMassInsertionMatching.h.

Public Member Functions

double B1 (double x) const
 \( \DeltaF = 1 \) loop functions, LO term More...
 
double B2 (double x) const
 \( \DeltaF = 1 \) loop functions, LO term More...
 
double C0LO (double x)
 \( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term More...
 
double C0NLO (double x, double mumatch2, double Ms2)
 \( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term More...
 
double C1LO (double x)
 \( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term More...
 
double C1NLO (double x, double mumatch2, double Ms2)
 \( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term More...
 
double C2LO (double x)
 \( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term More...
 
double C2NLO (double x, double mumatch2, double Ms2)
 \( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term More...
 
double C3LOA (double x)
 \( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term More...
 
double C3LOB (double x)
 \( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term More...
 
double C3NLOA (double x, double mumatch2, double Ms2)
 \( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term More...
 
double C3NLOB (double x, double mumatch2, double Ms2)
 \( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term More...
 
double C4LOA (double x)
 \( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term More...
 
double C4LOB (double x)
 \( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term More...
 
double C4NLOA (double x, double mumatch2, double Ms2)
 \( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term More...
 
double C4NLOB (double x, double mumatch2, double Ms2)
 \( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term More...
 
virtual std::vector< WilsonCoefficient > & CMd1 ()
 \( \Delta F = 1 \) More...
 
virtual std::vector< WilsonCoefficient > & CMdbd2 ()
 \( \Delta B = 2 \), \( B_{d} \) More...
 
virtual std::vector< WilsonCoefficient > & CMdbs2 ()
 \( \Delta B = 2 \), \( B_{s} \) More...
 
virtual std::vector< WilsonCoefficient > & CMdd2 ()
 \( \Delta C = 2 \) More...
 
virtual std::vector< WilsonCoefficient > & CMdk2 ()
 \( \Delta S = 2 \) More...
 
gslpp::matrix< double > EffectiveBase () const
 the basis is defined, for example, in Chetyrkin et al hep-ph/9612313 More...
 
double f6 (double x)
 
double f6t (double x)
 
void LRItoNDR (int i)
 it changes renormalization scheme from LRI to NDR More...
 
double M1 (double x) const
 \( \DeltaF = 1 \) loop functions, LO term More...
 
double M2 (double x) const
 \( \DeltaF = 1 \) loop functions, LO term More...
 
double M3 (double x) const
 \( \DeltaF = 1 \) loop functions, LO term More...
 
double M4 (double x) const
 \( \DeltaF = 1 \) loop functions, LO term More...
 
double P1 (double x) const
 \( \DeltaF = 1 \) loop functions, LO term More...
 
double P2 (double x) const
 \( \DeltaF = 1 \) loop functions, LO term More...
 
gslpp::matrix< double > RtoMisiak () const
 Misiak et al, hep-ph/0005183. More...
 
 SUSYMassInsertionMatching (const SUSYMassInsertion &SUSYMassInsertion_i)
 constructor More...
 
- Public Member Functions inherited from StandardModelMatching
 StandardModelMatching (const StandardModel &SM_i)
 
void updateSMParameters ()
 Updates to new Standard Model parameter sets. More...
 
virtual ~StandardModelMatching ()
 
- Public Member Functions inherited from ModelMatching
virtual std::vector< WilsonCoefficient > & CMBMll (QCD::lepton lepton)=0
 
virtual std::vector< WilsonCoefficient > & CMbnlep (int a)=0
 
virtual std::vector< WilsonCoefficient > & CMbnlepCC (const int a)=0
 
virtual std::vector< WilsonCoefficient > & CMbsg ()=0
 
virtual std::vector< WilsonCoefficient > & CMBXsnn (QCD::lepton lepton)=0
 
virtual std::vector< WilsonCoefficient > & CMd1Buras ()=0
 
virtual std::vector< WilsonCoefficientNew > & CMDF1 (std::string blocks, unsigned int nops)=0
 
virtual std::vector< WilsonCoefficient > & CMprimeBMll (QCD::lepton lepton)=0
 
virtual std::vector< WilsonCoefficient > & CMprimebsg ()=0
 
virtual ~ModelMatching ()
 

Private Attributes

gslpp::complex DLL
 
gslpp::complex DLR
 
gslpp::complex DRL
 
gslpp::matrix< double > drNDRLRI
 
gslpp::complex DRR
 
WilsonCoefficient mcbd
 
WilsonCoefficient mcbs
 
WilsonCoefficient mcd1
 
WilsonCoefficient mcd2
 
WilsonCoefficient mck2
 
double Ms2
 
double MuM2
 
unsigned int Nf
 
const SUSYMassInsertion & SusyMI
 

Constructor & Destructor Documentation

◆ SUSYMassInsertionMatching()

SUSYMassInsertionMatching::SUSYMassInsertionMatching ( const SUSYMassInsertion &  SUSYMassInsertion_i)

constructor

Parameters
SUSYMassInsertion_iobject of type SUSYMassInsertion

Definition at line 11 of file SUSYMassInsertionMatching.cpp.

11 :
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}
NLO
@ NLO
Definition: OrderScheme.h:35
NDR
@ NDR
Definition: OrderScheme.h:21
QCD::getNc
const double getNc() const
A get method to access the number of colours .
Definition: QCD.h:507
SUSYMassInsertionMatching::drNDRLRI
gslpp::matrix< double > drNDRLRI
Definition: SUSYMassInsertionMatching.h:278
SUSYMassInsertionMatching::SusyMI
const SUSYMassInsertion & SusyMI
Definition: SUSYMassInsertionMatching.h:274
StandardModelMatching::StandardModelMatching
StandardModelMatching(const StandardModel &SM_i)
Definition: StandardModelMatching.cpp:16

Member Function Documentation

◆ B1()

double SUSYMassInsertionMatching::B1 ( double  x) const

\( \DeltaF = 1 \) loop functions, LO term

Parameters
xthe square ratio between squark mass and gluino mass

Definition at line 39 of file SUSYMassInsertionMatching.cpp.

39 {
40 return ((1. + 4.*x - 5.*x*x + 4.*x*log(x) + 2.*x*x*log(x))/(8.*pow(x-1.,4.)));
41}

◆ B2()

double SUSYMassInsertionMatching::B2 ( double  x) const

\( \DeltaF = 1 \) loop functions, LO term

Parameters
xthe square ratio between squark mass and gluino mass

Definition at line 43 of file SUSYMassInsertionMatching.cpp.

43 {
44 return (x*(5. - 4.*x - x*x + 2.*log(x) + 4.*x*log(x))/(2.*pow(x-1.,4.)));
45}

◆ C0LO()

double SUSYMassInsertionMatching::C0LO ( double  x)

\( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term

Parameters
xthe square ratio between squark mass and gluino mass

Wilson Coefficients in RI-MOM scheme

Definition at line 91 of file SUSYMassInsertionMatching.cpp.

91 {
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}

◆ C0NLO()

double SUSYMassInsertionMatching::C0NLO ( double  x,
double  mumatch2,
double  Ms2 
)

\( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term

Parameters
xthe square ratio between squark mass and gluino mass
mumatch2double for the square of the SUSY matching scale
Ms2double for the square of the average squark mass

Wilson Coefficients in RI-MOM scheme

Definition at line 138 of file SUSYMassInsertionMatching.cpp.

138 {
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}
Li2
cd Li2(cd x)
Definition: hpl.h:1011

◆ C1LO()

double SUSYMassInsertionMatching::C1LO ( double  x)

\( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term

Parameters
xthe square ratio between squark mass and gluino mass

Definition at line 101 of file SUSYMassInsertionMatching.cpp.

101 {
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}

◆ C1NLO()

double SUSYMassInsertionMatching::C1NLO ( double  x,
double  mumatch2,
double  Ms2 
)

\( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term

Parameters
xthe square ratio between squark mass and gluino mass
mumatch2double for the square of the SUSY matching scale
Ms2double for the square of the average squark mass

Definition at line 248 of file SUSYMassInsertionMatching.cpp.

248 {
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}

◆ C2LO()

double SUSYMassInsertionMatching::C2LO ( double  x)

\( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term

Parameters
xthe square ratio between squark mass and gluino mass

Definition at line 106 of file SUSYMassInsertionMatching.cpp.

106 {
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}

◆ C2NLO()

double SUSYMassInsertionMatching::C2NLO ( double  x,
double  mumatch2,
double  Ms2 
)

\( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term

Parameters
xthe square ratio between squark mass and gluino mass
mumatch2double for the square of the SUSY matching scale
Ms2double for the square of the average squark mass

Definition at line 341 of file SUSYMassInsertionMatching.cpp.

341 {
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}

◆ C3LOA()

double SUSYMassInsertionMatching::C3LOA ( double  x)

\( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term

Parameters
xthe square ratio between squark mass and gluino mass

Definition at line 111 of file SUSYMassInsertionMatching.cpp.

111 {
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}

◆ C3LOB()

double SUSYMassInsertionMatching::C3LOB ( double  x)

\( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term

Parameters
xthe square ratio between squark mass and gluino mass

Definition at line 116 of file SUSYMassInsertionMatching.cpp.

116 {
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}

◆ C3NLOA()

double SUSYMassInsertionMatching::C3NLOA ( double  x,
double  mumatch2,
double  Ms2 
)

\( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term

Parameters
xthe square ratio between squark mass and gluino mass
mumatch2double for the square of the SUSY matching scale
Ms2double for the square of the average squark mass

Definition at line 433 of file SUSYMassInsertionMatching.cpp.

433 {
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}

◆ C3NLOB()

double SUSYMassInsertionMatching::C3NLOB ( double  x,
double  mumatch2,
double  Ms2 
)

\( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term

Parameters
xthe square ratio between squark mass and gluino mass
mumatch2double for the square of the SUSY matching scale
Ms2double for the square of the average squark mass

Definition at line 534 of file SUSYMassInsertionMatching.cpp.

534 {
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}

◆ C4LOA()

double SUSYMassInsertionMatching::C4LOA ( double  x)

\( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term

Parameters
xthe square ratio between squark mass and gluino mass

Definition at line 121 of file SUSYMassInsertionMatching.cpp.

121 {
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}

◆ C4LOB()

double SUSYMassInsertionMatching::C4LOB ( double  x)

\( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term

Parameters
xthe square ratio between squark mass and gluino mass

Definition at line 127 of file SUSYMassInsertionMatching.cpp.

127 {
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}

◆ C4NLOA()

double SUSYMassInsertionMatching::C4NLOA ( double  x,
double  mumatch2,
double  Ms2 
)

\( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term

Parameters
xthe square ratio between squark mass and gluino mass
mumatch2double for the square of the SUSY matching scale
Ms2double for the square of the average squark mass

Definition at line 635 of file SUSYMassInsertionMatching.cpp.

635 {
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}

◆ C4NLOB()

double SUSYMassInsertionMatching::C4NLOB ( double  x,
double  mumatch2,
double  Ms2 
)

\( \DeltaF = 2 \) loop functions RI-MOM scheme, LO term

Parameters
xthe square ratio between squark mass and gluino mass
mumatch2double for the square of the SUSY matching scale
Ms2double for the square of the average squark mass

Definition at line 734 of file SUSYMassInsertionMatching.cpp.

734 {
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}

◆ CMd1()

std::vector< WilsonCoefficient > & SUSYMassInsertionMatching::CMd1 ( )
virtual

\( \Delta F = 1 \)

Returns
the vector of 8 wilson coefficients: SM + SusyMI

Implements ModelMatching.

Definition at line 849 of file SUSYMassInsertionMatching.cpp.

849 {
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}
LO
@ LO
Definition: OrderScheme.h:34
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
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::getMuM
double getMuM() const
Definition: SUSYMassInsertion.h:236
SUSYMassInsertion::getM3
double getM3() const
Definition: SUSYMassInsertion.h:77
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
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::P2
double P2(double x) const
loop functions, LO term
Definition: SUSYMassInsertionMatching.cpp:52
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::EffectiveBase
gslpp::matrix< double > EffectiveBase() const
the basis is defined, for example, in Chetyrkin et al hep-ph/9612313
Definition: SUSYMassInsertionMatching.cpp:1324
SUSYMassInsertionMatching::M1
double M1(double x) const
loop functions, LO term
Definition: SUSYMassInsertionMatching.cpp:57
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
WilsonCoefficient::setCoeff
void setCoeff(const gslpp::vector< gslpp::complex > &z, orders order_i)
Definition: WilsonCoefficient.h:34
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

◆ CMdbd2()

std::vector< WilsonCoefficient > & SUSYMassInsertionMatching::CMdbd2 ( )
virtual

\( \Delta B = 2 \), \( B_{d} \)

Returns
the vector of 8 wilson coefficients: SM + SusyMI

Reimplemented from StandardModelMatching.

Definition at line 1045 of file SUSYMassInsertionMatching.cpp.

1045 {
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}
SUSYMassInsertion::getDd_RL
gslpp::matrix< gslpp::complex > getDd_RL() const
Definition: SUSYMassInsertion.h:205
SUSYMassInsertion::getDd_LL
gslpp::matrix< gslpp::complex > getDd_LL() const
Definition: SUSYMassInsertion.h:173
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::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::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::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::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::LRItoNDR
void LRItoNDR(int i)
it changes renormalization scheme from LRI to NDR
Definition: SUSYMassInsertionMatching.cpp:1353
SUSYMassInsertionMatching::C4LOA
double C4LOA(double x)
loop functions RI-MOM scheme, LO term
Definition: SUSYMassInsertionMatching.cpp:121
SUSYMassInsertionMatching::C4NLOB
double C4NLOB(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
Definition: SUSYMassInsertionMatching.cpp:734
StandardModelMatching::CMdbd2
virtual std::vector< WilsonCoefficient > & CMdbd2()
,
Definition: StandardModelMatching.cpp:966

◆ CMdbs2()

std::vector< WilsonCoefficient > & SUSYMassInsertionMatching::CMdbs2 ( )
virtual

\( \Delta B = 2 \), \( B_{s} \)

Returns
the vector of 8 wilson coefficients: SM + SusyMI

Reimplemented from StandardModelMatching.

Definition at line 1129 of file SUSYMassInsertionMatching.cpp.

1129 {
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}
StandardModelMatching::CMdbs2
virtual std::vector< WilsonCoefficient > & CMdbs2()
,
Definition: StandardModelMatching.cpp:1016

◆ CMdd2()

std::vector< WilsonCoefficient > & SUSYMassInsertionMatching::CMdd2 ( )
virtual

\( \Delta C = 2 \)

Returns
the vector of 8 wilson coefficients: SM + SusyMI

Reimplemented from StandardModelMatching.

Definition at line 965 of file SUSYMassInsertionMatching.cpp.

965 {
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}
StandardModelMatching::CMdd2
virtual std::vector< WilsonCoefficient > & CMdd2()
,
Definition: StandardModelMatching.cpp:1169

◆ CMdk2()

std::vector< WilsonCoefficient > & SUSYMassInsertionMatching::CMdk2 ( )
virtual

\( \Delta S = 2 \)

Returns
the vector of 8 wilson coefficients: SM + SusyMI

Definition at line 1214 of file SUSYMassInsertionMatching.cpp.

1214 {
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}

◆ EffectiveBase()

gslpp::matrix< double > SUSYMassInsertionMatching::EffectiveBase ( ) const

the basis is defined, for example, in Chetyrkin et al hep-ph/9612313

Returns
matrix for the change to effective basis

Definition at line 1324 of file SUSYMassInsertionMatching.cpp.

1324 {
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}

◆ f6()

double SUSYMassInsertionMatching::f6 ( double  x)

Definition at line 75 of file SUSYMassInsertionMatching.cpp.

75 {
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}

◆ f6t()

double SUSYMassInsertionMatching::f6t ( double  x)

Definition at line 80 of file SUSYMassInsertionMatching.cpp.

80 {
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}

◆ LRItoNDR()

void SUSYMassInsertionMatching::LRItoNDR ( int  i)

it changes renormalization scheme from LRI to NDR

Parameters
iint flag: 1 for D, 2 for B_d, 3 for B_s, 4 for K

Definition at line 1353 of file SUSYMassInsertionMatching.cpp.

1353 {
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}
WilsonCoefficient::getCoeff
gslpp::vector< gslpp::complex > ** getCoeff() const
Definition: WilsonCoefficient.h:29
WilsonTemplate::getMu
double getMu() const
Definition: WilsonTemplate.h:75

◆ M1()

double SUSYMassInsertionMatching::M1 ( double  x) const

\( \DeltaF = 1 \) loop functions, LO term

Parameters
xthe square ratio between squark mass and gluino mass

Definition at line 57 of file SUSYMassInsertionMatching.cpp.

57 {
58 return((1. + 4. * x - 5. * x*x + 4.*x*log(x) + 2.*x*x*log(x))/(2.*pow(x-1.,4.)));
59}

◆ M2()

double SUSYMassInsertionMatching::M2 ( double  x) const

\( \DeltaF = 1 \) loop functions, LO term

Parameters
xthe square ratio between squark mass and gluino mass

Definition at line 61 of file SUSYMassInsertionMatching.cpp.

61 {
62 return((-5. + 4.*x + x*x - 2.*log(x) - 4.*x*log(x))/(2.*pow(x-1.,4.)));
63}

◆ M3()

double SUSYMassInsertionMatching::M3 ( double  x) const

\( \DeltaF = 1 \) loop functions, LO term

Parameters
xthe square ratio between squark mass and gluino mass

Definition at line 65 of file SUSYMassInsertionMatching.cpp.

65 {
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}

◆ M4()

double SUSYMassInsertionMatching::M4 ( double  x) const

\( \DeltaF = 1 \) loop functions, LO term

Parameters
xthe square ratio between squark mass and gluino mass

Definition at line 69 of file SUSYMassInsertionMatching.cpp.

69 {
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}

◆ P1()

double SUSYMassInsertionMatching::P1 ( double  x) const

\( \DeltaF = 1 \) loop functions, LO term

Parameters
xthe square ratio between squark mass and gluino mass

Definition at line 47 of file SUSYMassInsertionMatching.cpp.

47 {
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}

◆ P2()

double SUSYMassInsertionMatching::P2 ( double  x) const

\( \DeltaF = 1 \) loop functions, LO term

Parameters
xthe square ratio between squark mass and gluino mass

Definition at line 52 of file SUSYMassInsertionMatching.cpp.

52 {
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}

◆ RtoMisiak()

gslpp::matrix< double > SUSYMassInsertionMatching::RtoMisiak ( ) const

Misiak et al, hep-ph/0005183.

Returns
matrix for the change of basis from the standard to Misiak one

Definition at line 1297 of file SUSYMassInsertionMatching.cpp.

1297 {
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}

Member Data Documentation

◆ DLL

gslpp::complex SUSYMassInsertionMatching::DLL
private

Definition at line 277 of file SUSYMassInsertionMatching.h.

◆ DLR

gslpp::complex SUSYMassInsertionMatching::DLR
private

Definition at line 277 of file SUSYMassInsertionMatching.h.

◆ DRL

gslpp::complex SUSYMassInsertionMatching::DRL
private

Definition at line 277 of file SUSYMassInsertionMatching.h.

◆ drNDRLRI

gslpp::matrix<double> SUSYMassInsertionMatching::drNDRLRI
private

Definition at line 278 of file SUSYMassInsertionMatching.h.

◆ DRR

gslpp::complex SUSYMassInsertionMatching::DRR
private

Definition at line 277 of file SUSYMassInsertionMatching.h.

◆ mcbd

WilsonCoefficient SUSYMassInsertionMatching::mcbd
private

Definition at line 280 of file SUSYMassInsertionMatching.h.

◆ mcbs

WilsonCoefficient SUSYMassInsertionMatching::mcbs
private

Definition at line 280 of file SUSYMassInsertionMatching.h.

◆ mcd1

WilsonCoefficient SUSYMassInsertionMatching::mcd1
private

Definition at line 280 of file SUSYMassInsertionMatching.h.

◆ mcd2

WilsonCoefficient SUSYMassInsertionMatching::mcd2
private

Definition at line 280 of file SUSYMassInsertionMatching.h.

◆ mck2

WilsonCoefficient SUSYMassInsertionMatching::mck2
private

Definition at line 280 of file SUSYMassInsertionMatching.h.

◆ Ms2

double SUSYMassInsertionMatching::Ms2
private

Definition at line 275 of file SUSYMassInsertionMatching.h.

◆ MuM2

double SUSYMassInsertionMatching::MuM2
private

Definition at line 275 of file SUSYMassInsertionMatching.h.

◆ Nf

unsigned int SUSYMassInsertionMatching::Nf
private

Definition at line 276 of file SUSYMassInsertionMatching.h.

◆ SusyMI

const SUSYMassInsertion& SUSYMassInsertionMatching::SusyMI
private

Definition at line 274 of file SUSYMassInsertionMatching.h.

The documentation for this class was generated from the following files: