SUSYMassInsertionMatching Class Reference

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

#include <SUSYMassInsertionMatching.h>

Inheritance diagram for SUSYMassInsertionMatching:
[legend]
Collaboration diagram for SUSYMassInsertionMatching:
[legend]

Detailed Description

A class for the matching of SUSY MIA.

Author
HEPfit Collaboration

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: 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
double A0t (double x) const
 loop function which appear in the Wilson coefficient for the magnetic operator in the non-effective Misiak basis, Bobeth et al hep-ph/9910220 More...
 
double A1t (double x, double mu) const
 loop function which appear in the Wilson coefficient for the semileptonic operator in the non-effective Misiak basis, Misiak and Urban hep-ph/9901278v1 More...
 
double B0b (double x) const
 loop functions for non-leptonic B decays, Buiras Basis Buras et al, hep-ph/9512380v1 More...
 
double B0t (double x) const
 loop function which appear in the Wilson coefficient for the semileptonic operator in the non-effective Misiak basis, Bobeth et al hep-ph/9910220 More...
 
double B1t (double x, double mu) const
 loop function which appear in the Wilson coefficient for the semileptonic operator in the non-effective Misiak basis, Misiak and Urban hep-ph/9901278v1 More...
 
double C0b (double x) const
 loop functions for non-leptonic B decays, Buiras Basis Buras et al, hep-ph/9512380v1 More...
 
double C0t (double x) const
 loop function which appear in the Wilson coefficient for the magnetic operator in the non-effective Misiak basis, Bobeth et al hep-ph/9910220 More...
 
double C1t (double x, double mu) const
 loop function which appear in the Wilson coefficient for the magnetic operator in the non-effective Misiak basis, Misiak and Urban hep-ph/9901278v1 More...
 
double C7LOeff (double x) const
 loop function which appear in the Wilson coefficient for the magnetic operator in the effective Misiak basis, LO term, Chetyrkin et al hep-ph/9612313 More...
 
double C7NLOeff (double x) const
 loop function which appear in the Wilson coefficient for the magnetic operator in the effective Misiak basis, NLO term, Chetyrkin et al hep-ph/9612313 More...
 
double C8LOeff (double x) const
 loop function which appear in the Wilson coefficient for the chromomagnetic operator in the effective Misiak basis, LO term, Chetyrkin et al hep-ph/9612313 More...
 
double C8NLOeff (double x) const
 loop function which appear in the Wilson coefficient for the chromomagnetic operator in the effective Misiak basis, LO term, Chetyrkin et al hep-ph/9612313 More...
 
virtual std::vector< WilsonCoefficient > & CMbdmm ()
 
virtual std::vector< WilsonCoefficient > & CMBMll ()
 operator basis: current current; qcd penguins; magnetic and chromomagnetic penguins; semileptonic More...
 
virtual std::vector< WilsonCoefficient > & CMbnlep (int a)
 operator basis: More...
 
virtual std::vector< WilsonCoefficient > & CMbnlepCC (int a)
 operator basis: - current current opertors More...
 
virtual std::vector< WilsonCoefficient > & CMbsg ()
 operator basis: current current; qcd penguins; magnetic and chromomagnetic penguins; semileptonic More...
 
virtual std::vector< WilsonCoefficient > & CMbsmm ()
 
virtual std::vector< WilsonCoefficient > & CMbtaunu ()
 
virtual std::vector< WilsonCoefficient > & CMBXdnn ()
 
virtual std::vector< WilsonCoefficient > & CMBXsnn ()
 
virtual std::vector< WilsonCoefficient > & CMd1Buras ()
 current-current oerators, Buras basis More...
 
virtual std::vector< WilsonCoefficient > & CMDLi3j (int li_lj)
 
virtual std::vector< WilsonCoefficient > & CMDLij (int li_lj)
 
virtual std::vector< WilsonCoefficient > & CMgminus2mu ()
 
virtual std::vector< WilsonCoefficient > & CMK ()
 operator basis: More...
 
virtual std::vector< WilsonCoefficient > & CMKCC ()
 operator basis: More...
 
virtual std::vector< WilsonCoefficient > & CMkmm ()
 
virtual std::vector< WilsonCoefficient > & CMkpnn ()
 
virtual std::vector< WilsonCoefficient > & CMmueconv ()
 
virtual std::vector< WilsonCoefficient > & CMprimeBMll ()
 operator basis: current current; qcd penguins; magnetic and chromomagnetic penguins; semileptonic More...
 
virtual std::vector< WilsonCoefficient > & CMprimebsg ()
 operator basis: current current; qcd penguins; magnetic and chromomagnetic penguins; semileptonic More...
 
double D0b (double x) const
 loop functions for non-leptonic B decays, Buiras Basis Buras et al, hep-ph/9512380v1 More...
 
double D0t (double x) const
 loop function which appear in the Wilson coefficient for the magnetic operator in the non-effective Misiak basis, Bobeth et al hep-ph/9910220 More...
 
double D1t (double x, double mu) const
 loop function which appear in the Wilson coefficient for the magnetic operator in the non-effective Misiak basis, Misiak and Urban hep-ph/9901278v1 More...
 
double E0b (double x) const
 loop functions for non-leptonic B decays, Buiras Basis Buras et al, hep-ph/9512380v1 More...
 
double E0t (double x) const
 loop function which appear in the Wilson coefficient for the chromomagnetic operator in the Misiak basis, Chetyrkin et al hep-ph/9612313 More...
 
double E1t (double x, double mu) const
 loop function which appear in the Wilson coefficient in the non-effective Misiak basis, Misiak and Urban hep-ph/9910220 More...
 
double Eet (double x) const
 loop function which appear in the Wilson coefficient in the non-effective Misiak basis, Misiak and Urban hep-ph/0512066 More...
 
double F0t (double x) const
 loop function which appear in the Wilson coefficient for the chromomagnetic operator in the non-effective Misiak basis, Bobeth et al hep-ph/9910220 More...
 
double F1t (double x, double mu) const
 loop function which appear in the Wilson coefficient for the semileptonic operator in the non-effective Misiak basis, Misiak and Urban hep-ph/9901278v1 More...
 
double G1t (double x, double mu) const
 loop function which appear in the Wilson coefficient in the non-effective Misiak basis, Misiak and Urban hep-ph/9910220 More...
 
double mt2omh2 (const double mu, const orders order=FULLNNLO) const
 
double Rest (double x, double mu) const
 approximation of two-loops EW correction for Q_10 operator in the non-effective Misiak basis, Misiak and Urban hep-ph/1311.1348 More...
 
double S0 (double, double) const
 
gslpp::complex S0c () const
 hep-ph/9512380 More...
 
gslpp::complex S0ct () const
 hep-ph/9512380 More...
 
gslpp::complex S0tt () const
 hep-ph/9512380v1 More...
 
double S1 (double x) const
 
 StandardModelMatching (const StandardModel &SM_i)
 
double Tt (double x) const
 loop function which appear in the Wilson coefficient in the non-effective Misiak basis, Misiak and Urban hep-ph/9910220 More...
 
void updateSMParameters ()
 Updates to new Standard Model parameter sets. More...
 
double Wt (double x) const
 loop function which appear in the Wilson coefficient in the non-effective Misiak basis, Misiak and Urban hep-ph/0512066 More...
 
double X0t (double x) const
 hep-ph/9512380v1 More...
 
double X1t (double x) const
 hep-ph/1009.0947v2 More...
 
double x_c (const double mu, const orders order=FULLNNLO) const
 
double x_t (const double mu, const orders order=FULLNNLO) const
 
double Xewt (double x, double a, double mu) const
 hep-ph/1009.0947v2 More...
 
double Y0 (double x) const
 
double Y1 (double x, double mu) const
 
- Public Member Functions inherited from ModelMatching
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 SUSYMassInsertionSusyMI
 

Additional Inherited Members

- Protected Attributes inherited from StandardModelMatching
std::vector< WilsonCoefficientvmcbdmm
 
std::vector< WilsonCoefficientvmcbdnn
 
std::vector< WilsonCoefficientvmcBMll
 
std::vector< WilsonCoefficientvmcbnlep
 
std::vector< WilsonCoefficientvmcbnlepCC
 
std::vector< WilsonCoefficientvmcbsg
 
std::vector< WilsonCoefficientvmcbsmm
 
std::vector< WilsonCoefficientvmcbsnn
 
std::vector< WilsonCoefficientvmcbtaunu
 
std::vector< WilsonCoefficientvmcd1
 
std::vector< WilsonCoefficientvmcd1Buras
 
std::vector< WilsonCoefficientvmcd2
 
std::vector< WilsonCoefficientvmcdb
 
std::vector< WilsonCoefficientvmcDLi3j
 
std::vector< WilsonCoefficientvmcDLij
 
std::vector< WilsonCoefficientvmcds
 
std::vector< WilsonCoefficientvmcgminus2mu
 
std::vector< WilsonCoefficientvmck
 
std::vector< WilsonCoefficientvmck2
 
std::vector< WilsonCoefficientvmckcc
 
std::vector< WilsonCoefficientvmckmm
 
std::vector< WilsonCoefficientvmckpnn
 
std::vector< WilsonCoefficientvmcmueconv
 
std::vector< WilsonCoefficientvmcprimeBMll
 
std::vector< WilsonCoefficientvmcprimebsg
 

Constructor & Destructor Documentation

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 }
const SUSYMassInsertion & SusyMI
double getNc() const
A get method to access the number of colours .
Definition: QCD.h:840
StandardModelMatching(const StandardModel &SM_i)
complex log(const complex &z)
gslpp::matrix< double > drNDRLRI

Member Function Documentation

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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
std::vector< WilsonCoefficient > & SUSYMassInsertionMatching::CMd1 ( )
virtual

\( \Delta F = 1 \)

Returns
the vector of 8 wilson coefficients: SM + SusyMI

Reimplemented from StandardModelMatching.

Definition at line 849 of file SUSYMassInsertionMatching.cpp.

849  {
850 
851  vmcd1.clear();
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 }
virtual void setMu(double mu)
const SUSYMassInsertion & SusyMI
double B2(double x) const
loop functions, LO term
double P2(double x) const
loop functions, LO term
complex pow(const complex &z1, const complex &z2)
double M4(double x) const
loop functions, LO term
double Als(const double mu, const orders order=FULLNLO) const
Computes the running strong coupling in the scheme. In the cases of LO, NLO and FULLNNLO...
Definition: QCD.cpp:1004
double getMuM() const
double getMsq() const
schemes getScheme() const
void setCoeff(const gslpp::vector< gslpp::complex > &z, orders order_i)
gslpp::matrix< gslpp::complex > getDu_LR() const
virtual std::vector< WilsonCoefficient > & CMd1()
current-current oerators, Misiak basis
double P1(double x) const
loop functions, LO term
Definition: OrderScheme.h:33
gslpp::matrix< gslpp::complex > getDu_LL() const
gslpp::matrix< double > EffectiveBase() const
the basis is defined, for example, in Chetyrkin et al hep-ph/9612313
gslpp::matrix< double > RtoMisiak() const
Misiak et al, hep-ph/0005183.
double M2(double x) const
loop functions, LO term
double M1(double x) const
loop functions, LO term
double getMuc() const
A get method to access the threshold between four- and three-flavour theory in GeV.
Definition: QCD.h:914
std::vector< WilsonCoefficient > vmcd1
double B1(double x) const
loop functions, LO term
orders getOrder() const
double M3(double x) const
loop functions, LO term
double getM3() const
gslpp::matrix< gslpp::complex > getDu_RL() const
gslpp::matrix< gslpp::complex > getDu_RR() const
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 
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 }
double C4LOA(double x)
loop functions RI-MOM scheme, LO term
double C3NLOA(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
double C4NLOA(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
virtual void setMu(double mu)
double C1NLO(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
double C0NLO(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
const SUSYMassInsertion & SusyMI
gslpp::matrix< gslpp::complex > getDd_LL() const
double C4NLOB(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
complex pow(const complex &z1, const complex &z2)
gslpp::matrix< gslpp::complex > getDd_LR() const
double C4LOB(double x)
loop functions RI-MOM scheme, LO term
double Als(const double mu, const orders order=FULLNLO) const
Computes the running strong coupling in the scheme. In the cases of LO, NLO and FULLNNLO...
Definition: QCD.cpp:1004
double getMuM() const
virtual std::vector< WilsonCoefficient > & CMdbd2()
,
double getMsq() const
gslpp::matrix< gslpp::complex > getDd_RR() const
schemes getScheme() const
void setCoeff(const gslpp::vector< gslpp::complex > &z, orders order_i)
double C2NLO(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
double C3LOA(double x)
loop functions RI-MOM scheme, LO term
gslpp::matrix< gslpp::complex > getDd_RL() const
std::vector< WilsonCoefficient > vmcdb
void LRItoNDR(int i)
it changes renormalization scheme from LRI to NDR
Definition: OrderScheme.h:33
double C3NLOB(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
double C3LOB(double x)
loop functions RI-MOM scheme, LO term
double C0LO(double x)
loop functions RI-MOM scheme, LO term
double C2LO(double x)
loop functions RI-MOM scheme, LO term
orders getOrder() const
double getM3() const
double C1LO(double x)
loop functions RI-MOM scheme, LO term
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 
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 }
double C4LOA(double x)
loop functions RI-MOM scheme, LO term
double C3NLOA(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
double C4NLOA(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
virtual void setMu(double mu)
double C1NLO(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
double C0NLO(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
const SUSYMassInsertion & SusyMI
gslpp::matrix< gslpp::complex > getDd_LL() const
double C4NLOB(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
complex pow(const complex &z1, const complex &z2)
gslpp::matrix< gslpp::complex > getDd_LR() const
double C4LOB(double x)
loop functions RI-MOM scheme, LO term
double Als(const double mu, const orders order=FULLNLO) const
Computes the running strong coupling in the scheme. In the cases of LO, NLO and FULLNNLO...
Definition: QCD.cpp:1004
double getMuM() const
double getMsq() const
gslpp::matrix< gslpp::complex > getDd_RR() const
schemes getScheme() const
void setCoeff(const gslpp::vector< gslpp::complex > &z, orders order_i)
double C2NLO(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
double C3LOA(double x)
loop functions RI-MOM scheme, LO term
gslpp::matrix< gslpp::complex > getDd_RL() const
void LRItoNDR(int i)
it changes renormalization scheme from LRI to NDR
Definition: OrderScheme.h:33
double C3NLOB(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
virtual std::vector< WilsonCoefficient > & CMdbs2()
,
std::vector< WilsonCoefficient > vmcds
double C3LOB(double x)
loop functions RI-MOM scheme, LO term
double C0LO(double x)
loop functions RI-MOM scheme, LO term
double C2LO(double x)
loop functions RI-MOM scheme, LO term
orders getOrder() const
double getM3() const
double C1LO(double x)
loop functions RI-MOM scheme, LO term
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();
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;
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 }
double C4LOA(double x)
loop functions RI-MOM scheme, LO term
std::vector< WilsonCoefficient > vmcd2
double C3NLOA(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
double C4NLOA(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
virtual void setMu(double mu)
double C1NLO(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
double C0NLO(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
const SUSYMassInsertion & SusyMI
double C4NLOB(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
complex pow(const complex &z1, const complex &z2)
double C4LOB(double x)
loop functions RI-MOM scheme, LO term
double Als(const double mu, const orders order=FULLNLO) const
Computes the running strong coupling in the scheme. In the cases of LO, NLO and FULLNNLO...
Definition: QCD.cpp:1004
double getMuM() const
double getMsq() const
schemes getScheme() const
void setCoeff(const gslpp::vector< gslpp::complex > &z, orders order_i)
gslpp::matrix< gslpp::complex > getDu_LR() const
double C2NLO(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
double C3LOA(double x)
loop functions RI-MOM scheme, LO term
void LRItoNDR(int i)
it changes renormalization scheme from LRI to NDR
Definition: OrderScheme.h:33
gslpp::matrix< gslpp::complex > getDu_LL() const
double C3NLOB(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
double C3LOB(double x)
loop functions RI-MOM scheme, LO term
double C0LO(double x)
loop functions RI-MOM scheme, LO term
double C2LO(double x)
loop functions RI-MOM scheme, LO term
orders getOrder() const
double getM3() const
gslpp::matrix< gslpp::complex > getDu_RL() const
double C1LO(double x)
loop functions RI-MOM scheme, LO term
virtual std::vector< WilsonCoefficient > & CMdd2()
,
gslpp::matrix< gslpp::complex > getDu_RR() const
std::vector< WilsonCoefficient > & SUSYMassInsertionMatching::CMdk2 ( )
virtual

\( \Delta S = 2 \)

Returns
the vector of 8 wilson coefficients: SM + SusyMI

Reimplemented from StandardModelMatching.

Definition at line 1214 of file SUSYMassInsertionMatching.cpp.

1214  {
1215 
1216  vmck2.clear();
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 }
double C4LOA(double x)
loop functions RI-MOM scheme, LO term
double C3NLOA(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
double C4NLOA(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
virtual void setMu(double mu)
double C1NLO(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
double C0NLO(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
std::vector< WilsonCoefficient > vmck2
const SUSYMassInsertion & SusyMI
gslpp::matrix< gslpp::complex > getDd_LL() const
double C4NLOB(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
complex pow(const complex &z1, const complex &z2)
gslpp::matrix< gslpp::complex > getDd_LR() const
double C4LOB(double x)
loop functions RI-MOM scheme, LO term
double Als(const double mu, const orders order=FULLNLO) const
Computes the running strong coupling in the scheme. In the cases of LO, NLO and FULLNNLO...
Definition: QCD.cpp:1004
double getMuM() const
double getMsq() const
gslpp::matrix< gslpp::complex > getDd_RR() const
schemes getScheme() const
void setCoeff(const gslpp::vector< gslpp::complex > &z, orders order_i)
double C2NLO(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
double C3LOA(double x)
loop functions RI-MOM scheme, LO term
gslpp::matrix< gslpp::complex > getDd_RL() const
void LRItoNDR(int i)
it changes renormalization scheme from LRI to NDR
Definition: OrderScheme.h:33
double C3NLOB(double x, double mumatch2, double Ms2)
loop functions RI-MOM scheme, LO term
double C3LOB(double x)
loop functions RI-MOM scheme, LO term
double C0LO(double x)
loop functions RI-MOM scheme, LO term
double C2LO(double x)
loop functions RI-MOM scheme, LO term
virtual std::vector< WilsonCoefficient > & CMdk2()
orders getOrder() const
double getM3() const
double C1LO(double x)
loop functions RI-MOM scheme, LO term
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 }
A class for constructing and defining operations on real matrices.
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
const SUSYMassInsertion & SusyMI
double Als(const double mu, const orders order=FULLNLO) const
Computes the running strong coupling in the scheme. In the cases of LO, NLO and FULLNNLO...
Definition: QCD.cpp:1004
gslpp::vector< gslpp::complex > ** getCoeff() const
void setCoeff(const gslpp::vector< gslpp::complex > &z, orders order_i)
Definition: OrderScheme.h:33
double getMu() const
gslpp::matrix< double > drNDRLRI
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
complex pow(const complex &z1, const complex &z2)
complex log(const complex &z)
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 }
A class for constructing and defining operations on real matrices.

Member Data Documentation

gslpp::complex SUSYMassInsertionMatching::DLL
private

Definition at line 277 of file SUSYMassInsertionMatching.h.

gslpp::complex SUSYMassInsertionMatching::DLR
private

Definition at line 277 of file SUSYMassInsertionMatching.h.

gslpp::complex SUSYMassInsertionMatching::DRL
private

Definition at line 277 of file SUSYMassInsertionMatching.h.

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

Definition at line 278 of file SUSYMassInsertionMatching.h.

gslpp::complex SUSYMassInsertionMatching::DRR
private

Definition at line 277 of file SUSYMassInsertionMatching.h.

WilsonCoefficient SUSYMassInsertionMatching::mcbd
private

Definition at line 280 of file SUSYMassInsertionMatching.h.

WilsonCoefficient SUSYMassInsertionMatching::mcbs
private

Definition at line 280 of file SUSYMassInsertionMatching.h.

WilsonCoefficient SUSYMassInsertionMatching::mcd1
private

Definition at line 280 of file SUSYMassInsertionMatching.h.

WilsonCoefficient SUSYMassInsertionMatching::mcd2
private

Definition at line 280 of file SUSYMassInsertionMatching.h.

WilsonCoefficient SUSYMassInsertionMatching::mck2
private

Definition at line 280 of file SUSYMassInsertionMatching.h.

double SUSYMassInsertionMatching::Ms2
private

Definition at line 275 of file SUSYMassInsertionMatching.h.

double SUSYMassInsertionMatching::MuM2
private

Definition at line 275 of file SUSYMassInsertionMatching.h.

unsigned int SUSYMassInsertionMatching::Nf
private

Definition at line 276 of file SUSYMassInsertionMatching.h.

const SUSYMassInsertion& SUSYMassInsertionMatching::SusyMI
private

Definition at line 274 of file SUSYMassInsertionMatching.h.


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