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

A class for cache variables used in computing radiative corrections to the EW precision observables. More...

#include <EWSMcache.h>

Detailed Description

A class for cache variables used in computing radiative corrections to the EW precision observables.

Author
HEPfit Collaboration
Copyright
GNU General Public License

This class provides caching methods for a bunch of the functions appearing in EWSMOneLoopEW, EWSMTwoLoopQCD, EWSMTwoLoopEW, EWSMThreeLoopQCD, EWSMThreeLoopEW2QCD and EWSMThreeLoopEW classes. Each caching method calls the private member function CacheCheck() and newCacheForDouble() (or newCacheForComplex()).

Moreover, this class contains methods to access model parameters and functions in QCD and StandardModel class, such that the classes listed above do not call directly the functions in QCD nor StandardModel.

The internal flags FlagDebug and FlagCacheInEWSMcache, which can be changed with setFlagDebug() and setFlagCacheInEWSMcache(), respectively, are designed for debugging. The latter flag can be controlled with the model flag CacheInEWSMcache of StandardModel.

Definition at line 40 of file EWSMcache.h.

Public Member Functions

double A0_Mw2_mh2 (const double Mw_i) const
 A cache method. More...
 
double A0_Mw2_Mw2 (const double Mw_i) const
 A cache method. More...
 
double A0_Mw2_Mz2 (const double Mw_i) const
 A cache method. More...
 
double A0_Mz2_mh2 () const
 A cache method. More...
 
double A0_Mz2_Mw2 (const double Mw_i) const
 A cache method. More...
 
double A0_Mz2_Mz2 () const
 A cache method. More...
 
double a_f (const Particle f) const
 The tree-level axial-vector coupling for \(Z\to f\bar{f}\), denoted as \(a_f\). More...
 
double Als (const double mu, const orders order) const
 The strong coupling \(\alpha_s(\mu)\). More...
 
double alsMt () const
 The strong coupling \(\alpha_s(m_t^2)\) at NNLO. More...
 
gslpp::complex B0_Mw2_Mw2_0_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0_Mw2_Mw2_mh2_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0_Mw2_Mw2_Mz2_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0_Mw2_Mz2_Mt2_Mt2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0_Mw2_Mz2_Mw2_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0_Mz2_0_0_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0_Mz2_0_mh2_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0_Mz2_0_Mz2_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0_Mz2_Mw2_0_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0_Mz2_Mw2_mh2_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0_Mz2_Mw2_Mz2_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0_Mz2_Mz2_mf2_mf2 (const Particle f) const
 A cache method. More...
 
gslpp::complex B0_Mz2_Mz2_mh2_Mz2 () const
 A cache method. More...
 
gslpp::complex B0_Mz2_Mz2_Mw2_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0p_Mw2_Mw2_0_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0p_Mw2_Mw2_mh2_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0p_Mw2_Mw2_Mz2_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0p_Mz2_0_mh2_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0p_Mz2_0_Mz2_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B0p_Mz2_Mz2_mf2_mf2 (const Particle f) const
 A cache method. More...
 
gslpp::complex B0p_Mz2_Mz2_mh2_Mz2 () const
 A cache method. More...
 
gslpp::complex B0p_Mz2_Mz2_Mw2_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex B1_Mz2_0_mf2_mfprime2 (const int gen) const
 A cache method. More...
 
gslpp::complex B1_Mz2_0_mfprime2_mf2 (const int gen) const
 A cache method. More...
 
gslpp::complex B1_Mz2_Mw2_mf2_mfprime2 (const int gen, const double Mw_i) const
 A cache method. More...
 
gslpp::complex B1_Mz2_Mw2_mfprime2_mf2 (const int gen, const double Mw_i) const
 A cache method. More...
 
gslpp::complex B1p_Mw2_Mw2_mf2_mfprime2 (const int gen, const double Mw_i) const
 A cache method. More...
 
gslpp::complex B1p_Mw2_Mw2_mfprime2_mf2 (const int gen, const double Mw_i) const
 A cache method. More...
 
gslpp::complex Bf_Mw2_Mw2_mfprime2_mf2 (const int gen, const double Mw_i) const
 A cache method. More...
 
gslpp::complex Bf_Mz2_0_mf2_mf2 (const Particle f) const
 A cache method. More...
 
gslpp::complex Bf_Mz2_0_mfprime2_mf2 (const int gen) const
 A cache method. More...
 
gslpp::complex Bf_Mz2_Mw2_mfprime2_mf2 (const int gen, const double Mw_i) const
 A cache method. More...
 
gslpp::complex Bf_Mz2_Mz2_mf2_mf2 (const Particle f) const
 A cache method. More...
 
gslpp::complex Bfp_Mw2_Mw2_mfprime2_mf2 (const int gen, const double Mw_i) const
 A cache method. More...
 
gslpp::complex Bfp_Mz2_Mz2_mf2_mf2 (const Particle f) const
 A cache method. More...
 
gslpp::complex C0_Mw2_0_Mz2_0 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex C0_Mw2_Mw2_0_Mz2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex C0_Mz2_0_Mw2_0 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex C0_Mz2_0_Mz2_0 () const
 A cache method. More...
 
gslpp::complex C0_Mz2_Mt2_Mw2_Mt2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex C0_Mz2_Mw2_0_Mw2 (const double Mw_i) const
 A cache method. More...
 
gslpp::complex C0_Mz2_Mw2_Mt2_Mw2 (const double Mw_i) const
 A cache method. More...
 
double Cl2_2Phi () const
 The constant \({\rm Cl}_2(2 \Phi)\). More...
 
double Cl2_4Phi () const
 The constant \({\rm Cl}_2(4 \Phi)\). More...
 
double Cl3_2Phi () const
 The constant \({\rm Cl}_3(2 \Phi)\). More...
 
double Cl3_4Phi () const
 The constant \({\rm Cl}_3(4 \Phi)\). More...
 
double delta_f (const Particle f, const double Mw_i) const
 \(\delta_f = v_f-a_f\). More...
 
 EWSMcache (const StandardModel &SM_i)
 Constructor. More...
 
double f_AlphaToGF (const double Mw_i) const
 The conversion factor from \(\alpha\) to \(G_\mu\), denoted as \(f\). More...
 
double gamma_QCD2 () const
 The constant \(\gamma\) for two-loop QCD contribution. More...
 
double getB4 () const
 A get method to access the constant \(B_4\). More...
 
const ClausenFunctions getClausen () const
 A get method to access the member object of type ClausenFunctions. More...
 
double getD3 () const
 A get method to access the constant \(D_3\). More...
 
double getLog2 () const
 A get method to access the constant \(\ln 2\). More...
 
const Polylogarithms getPolyLog () const
 A get method to access the member object of type Polylogarithms. More...
 
const PVfunctions getPV () const
 A get method to access the member reference to the object of type StandardModel passed to the constructor. More...
 
double getS2 () const
 A get method to access the constant \(S_2\). More...
 
const StandardModel & getSM () const
 
double getZeta2 () const
 A get method to access the value of the zeta function \(\zeta(2)\). More...
 
double getZeta3 () const
 A get method to access the value of the zeta function \(\zeta(3)\). More...
 
double getZeta4 () const
 A get method to access the value of the zeta function \(\zeta(4)\). More...
 
double getZeta5 () const
 A get method to access the value of the zeta function \(\zeta(5)\). More...
 
double h_QCD2 () const
 The constant \(h\) for two-loop QCD contribution. More...
 
double I3_f (const Particle f) const
 The isospin of an SM fermion \(I_3^f\). More...
 
double Li2_MW2toMTOP2 (const double Mw_i) const
 A cache method. More...
 
double Li3_for_F1 (const double Mw_i) const
 A cache method. More...
 
double Li3_MW2toMTOP2 (const double Mw_i) const
 A cache method. More...
 
double log_cW2 (const double Mw_i) const
 A cache method. More...
 
double logMTOPtoMH () const
 A cache method. More...
 
double logMZtoME () const
 A cache method. More...
 
double logMZtoMMU () const
 A cache method. More...
 
double logMZtoMTAU () const
 A cache method. More...
 
double logMZtoMTOP () const
 A cache method. More...
 
double logV1primeAndA1prime () const
 A logarithm appearing in the functions \(V_1'\) and \(A_1'\) for two-loop QCD contribution. More...
 
double mf (const Particle f, const double mu=0.0, const orders order=FULLNNLO) const
 The mass of an SM fermion. More...
 
double mf2 (const Particle f, const double mu=0.0, const orders order=FULLNNLO) const
 The mass squared of an SM fermion. More...
 
double Phi_QCD2 () const
 The constant \(\Phi\) for two-loop QCD contribution. More...
 
double Q_f (const Particle f) const
 The charge of an SM fermion \(Q_f\). More...
 
void setFlagCacheInEWSMcache (bool FlagCacheInEWSMcache)
 A set method to change the model flag CacheInEWSMcache in StandardModel. More...
 
void setFlagDebug (bool FlagDebug)
 A set method to change the internal boolean flag FlagDebug. More...
 
double sigma_f (const Particle f, const double Mw_i) const
 \(\sigma_f = |v_f+a_f|\). More...
 
double v_f (const Particle f, const double Mw_i) const
 The tree-level vector coupling for \(Z\to f\bar{f}\), denoted as \(v_f\). More...
 
double Xt_alpha (const double Mw_i) const
 The quantity \(X_t\) with the coupling \(\alpha\). More...
 
double Xt_GF () const
 The quantity \(X_t\) with the coupling \(G_\mu\). More...
 

Private Member Functions

bool CacheCheck (const double cache[], const int NumPar, const double params[]) const
 A check method for caching. More...
 
void newCacheForComplex (double cache[], const int NumPar, const double params[], const gslpp::complex newResult) const
 A method to update a cache of the parameters and the quantity under consideration. More...
 
void newCacheForDouble (double cache[], const int NumPar, const double params[], const double newResult) const
 A method to update a cache of the parameters and the quantity under consideration. More...
 

Private Attributes

double A0_Mw2_mh2_cache [3]
 A cache of a PV function. More...
 
double A0_Mw2_Mw2_cache [2]
 A cache of a PV function. More...
 
double A0_Mw2_Mz2_cache [3]
 A cache of a PV function. More...
 
double A0_Mz2_mh2_cache [3]
 A cache of a PV function. More...
 
double A0_Mz2_Mw2_cache [3]
 A cache of a PV function. More...
 
double A0_Mz2_Mz2_cache [2]
 A cache of a PV function. More...
 
double B0_Mw2_Mw2_0_Mw2_cache [3]
 A cache of a PV function. More...
 
double B0_Mw2_Mw2_mh2_Mw2_cache [4]
 A cache of a PV function. More...
 
double B0_Mw2_Mw2_Mz2_Mw2_cache [4]
 A cache of a PV function. More...
 
double B0_Mw2_Mz2_Mt2_Mt2_cache [5]
 A cache of a PV function. More...
 
double B0_Mw2_Mz2_Mw2_Mw2_cache [4]
 A cache of a PV function. More...
 
double B0_Mz2_0_0_Mw2_cache [4]
 A cache of a PV function. More...
 
double B0_Mz2_0_mh2_Mw2_cache [5]
 A cache of a PV function. More...
 
double B0_Mz2_0_Mz2_Mw2_cache [4]
 A cache of a PV function. More...
 
double B0_Mz2_Mw2_0_Mw2_cache [4]
 A cache of a PV function. More...
 
double B0_Mz2_Mw2_mh2_Mw2_cache [5]
 A cache of a PV function. More...
 
double B0_Mz2_Mw2_Mz2_Mw2_cache [4]
 A cache of a PV function. More...
 
double B0_Mz2_Mz2_mf2_mf2_cache [12][4]
 A cache of a PV function. More...
 
double B0_Mz2_Mz2_mh2_Mz2_cache [4]
 A cache of a PV function. More...
 
double B0_Mz2_Mz2_Mw2_Mw2_cache [4]
 A cache of a PV function. More...
 
double B0p_Mw2_Mw2_0_Mw2_cache [3]
 A cache of a PV function. More...
 
double B0p_Mw2_Mw2_mh2_Mw2_cache [4]
 A cache of a PV function. More...
 
double B0p_Mw2_Mw2_Mz2_Mw2_cache [4]
 A cache of a PV function. More...
 
double B0p_Mz2_0_mh2_Mw2_cache [5]
 A cache of a PV function. More...
 
double B0p_Mz2_0_Mz2_Mw2_cache [4]
 A cache of a PV function. More...
 
double B0p_Mz2_Mz2_mf2_mf2_cache [12][4]
 A cache of a PV function. More...
 
double B0p_Mz2_Mz2_mh2_Mz2_cache [4]
 A cache of a PV function. More...
 
double B0p_Mz2_Mz2_Mw2_Mw2_cache [4]
 A cache of a PV function. More...
 
double B1_Mz2_0_mf2_mfprime2_cache [6][5]
 A cache of a PV function. More...
 
double B1_Mz2_0_mfprime2_mf2_cache [6][5]
 A cache of a PV function. More...
 
double B1_Mz2_Mw2_mf2_mfprime2_cache [6][6]
 A cache of a PV function. More...
 
double B1_Mz2_Mw2_mfprime2_mf2_cache [6][6]
 A cache of a PV function. More...
 
double B1p_Mw2_Mw2_mf2_mfprime2_cache [6][5]
 A cache of a PV function. More...
 
double B1p_Mw2_Mw2_mfprime2_mf2_cache [6][5]
 A cache of a PV function. More...
 
double B4
 The constant \(B_4=16\mathrm{Li}_4(1/2) - 4\zeta(2)\ln^2(2) + (2/3)\ln^4(2) - (13/2)\zeta(4)\). More...
 
double Bf_Mw2_Mw2_mfprime2_mf2_cache [6][5]
 A cache of a PV function. More...
 
double Bf_Mz2_0_mf2_mf2_cache [12][4]
 A cache of a PV function. More...
 
double Bf_Mz2_0_mfprime2_mf2_cache [6][5]
 A cache of a PV function. More...
 
double Bf_Mz2_Mw2_mfprime2_mf2_cache [6][6]
 A cache of a PV function. More...
 
double Bf_Mz2_Mz2_mf2_mf2_cache [12][4]
 A cache of a PV function. More...
 
double Bfp_Mw2_Mw2_mfprime2_mf2_cache [6][5]
 A cache of a PV function. More...
 
double Bfp_Mz2_Mz2_mf2_mf2_cache [12][4]
 A cache of a PV function. More...
 
double C0_Mw2_0_Mz2_0_cache [4]
 A cache of a PV function. More...
 
double C0_Mw2_Mw2_0_Mz2_cache [4]
 A cache of a PV function. More...
 
double C0_Mz2_0_Mw2_0_cache [4]
 A cache of a PV function. More...
 
double C0_Mz2_0_Mz2_0_cache [3]
 A cache of a PV function. More...
 
double C0_Mz2_Mt2_Mw2_Mt2_cache [5]
 A cache of a PV function. More...
 
double C0_Mz2_Mw2_0_Mw2_cache [4]
 A cache of a PV function. More...
 
double C0_Mz2_Mw2_Mt2_Mw2_cache [5]
 A cache of a PV function. More...
 
const ClausenFunctions Clausen
 An object of type ClausenFunctions. More...
 
double D3
 The constant \(D_3=6\zeta(3) - (15/4)\zeta(4) - 6[\mathrm{Cl}_2(\pi/3)]^2\). More...
 
bool FlagCacheInEWSMcache
 A flag for caching (true by default). More...
 
bool FlagDebug
 A flag for debugging (false by default). More...
 
double Li2_MW2toMTOP2_cache [3]
 A cache of \(\mathrm{Li}_2(M_W^2/m_t^2)\). More...
 
double Li3_for_F1_cache [3]
 A cache of \(\mathrm{Li}_3(-M_W^2/m_t^2/(1-M_W^2/m_t^2))\). More...
 
double Li3_MW2toMTOP2_cache [3]
 A cache of \(\mathrm{Li}_3(M_W^2/m_t^2)\). More...
 
double log2
 The constant \(\ln 2\). More...
 
double log_cW2_cache [3]
 A cache of \(\ln c_W^2\). More...
 
double logMTOPtoMH_cache [3]
 A cache of \(\ln (m_t/m_h)\). More...
 
double logMZtoME_cache [3]
 A cache of \(\ln (M_Z/m_e)\). More...
 
double logMZtoMMU_cache [3]
 A cache of \(\ln (M_Z/m_\mu)\). More...
 
double logMZtoMTAU_cache [3]
 A cache of \(\ln (M_Z/m_\tau)\). More...
 
double logMZtoMTOP_cache [3]
 A cache of \(\ln (M_Z/m_t)\). More...
 
double mf_atMz_cache [12]
 A cache of the fermion masses at \(\mu=M_Z\). More...
 
double mf_atMz_params_cache [12][StandardModel::NumSMParamsForEWPO]
 A cache array of a set of SM parameters, used together with mf_atMz_cache. More...
 
const Polylogarithms PolyLog
 An object of type Polylogarithms. More...
 
const PVfunctions PV
 An object of type PVfunctions. More...
 
double S2
 The constant \(S_2=(4/9/\sqrt{3})\mathrm{Cl}_2(\pi/3)\). More...
 
const StandardModel & SM
 A reference to an object of type StandardModel. More...
 
double zeta2
 The constant \(\zeta(2)\). More...
 
double zeta3
 The constant \(\zeta(3)\). More...
 
double zeta4
 The constant \(\zeta(4)\). More...
 
double zeta5
 The constant \(\zeta(5)\). More...
 

Constructor & Destructor Documentation

◆ EWSMcache()

EWSMcache::EWSMcache ( const StandardModel &  SM_i)

Constructor.

Parameters
[in]SM_ia reference to an object of type StandardModel

Definition at line 15 of file EWSMcache.cpp.

16 : SM(SM_i), PV(true)
17 {
18  FlagDebug = false;
19  FlagCacheInEWSMcache = true; // use caches in the current class
20  //FlagCacheInEWSMcache = false;// do not use caches in the current class (for test)
21 
22  log2 = log(2.0);
23 
24  /* zeta functions */
25  zeta2 = gsl_sf_zeta_int(2);
26  zeta3 = gsl_sf_zeta_int(3);
27  zeta4 = gsl_sf_zeta_int(4);
28  zeta5 = gsl_sf_zeta_int(5);
29 
30  /* Constants for three-loop contribution */
31  double Cl2_Pi_3 = Clausen.Cl2(M_PI / 3.0);
32  S2 = 4.0 / 9.0 / sqrt(3.0) * Cl2_Pi_3;
33  D3 = 6.0 * zeta3 - 15.0 / 4.0 * zeta4 - 6.0 * Cl2_Pi_3*Cl2_Pi_3;
34  B4 = -1.76280008707377;
35  //double Li4_1_2 = ??;
36  //B4 = 16.0*Li4_1_2 - 4.0*zeta2*log2*log2 + 2.0/3.0*pow(log2,4.0) - 13.0/2.0*zeta4;
37 
38  // Initializations of the cache
39  for (int i = 0; i < 12; ++i) {
40  mf_atMz_cache[i] = 0.0;
41  for (int j = 0; j < StandardModel::NumSMParamsForEWPO; ++j)
42  mf_atMz_params_cache[i][j] = 0.0;
43  }
44 }

Member Function Documentation

◆ A0_Mw2_mh2()

double EWSMcache::A0_Mw2_mh2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(A_0(m_h^2)\) with \(\mu=M_W\)

Definition at line 239 of file EWSMcache.cpp.

240 {
241  int NumPar = 2;
242  double params[] = {Mw_i, SM.getMHl()};
243 
244  if (CacheCheck(A0_Mw2_mh2_cache, NumPar, params))
245  return A0_Mw2_mh2_cache[NumPar];
246  else {
247  double newResult = PV.A0(Mw_i*Mw_i, SM.getMHl() * SM.getMHl());
248  newCacheForDouble(A0_Mw2_mh2_cache, NumPar, params, newResult);
249  return newResult;
250  }
251 }

◆ A0_Mw2_Mw2()

double EWSMcache::A0_Mw2_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(A_0(M_W^2)\) with \(\mu=M_W\)

Definition at line 267 of file EWSMcache.cpp.

268 {
269  int NumPar = 1;
270  double params[] = {Mw_i};
271 
272  if (CacheCheck(A0_Mw2_Mw2_cache, NumPar, params))
273  return A0_Mw2_Mw2_cache[NumPar];
274  else {
275  double newResult = PV.A0(Mw_i*Mw_i, Mw_i * Mw_i);
276  newCacheForDouble(A0_Mw2_Mw2_cache, NumPar, params, newResult);
277  return newResult;
278  }
279 }

◆ A0_Mw2_Mz2()

double EWSMcache::A0_Mw2_Mz2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(A_0(M_Z^2)\) with \(\mu=M_W\)

Definition at line 225 of file EWSMcache.cpp.

226 {
227  int NumPar = 2;
228  double params[] = {Mw_i, SM.getMz()};
229 
230  if (CacheCheck(A0_Mw2_Mz2_cache, NumPar, params))
231  return A0_Mw2_Mz2_cache[NumPar];
232  else {
233  double newResult = PV.A0(Mw_i*Mw_i, SM.getMz() * SM.getMz());
234  newCacheForDouble(A0_Mw2_Mz2_cache, NumPar, params, newResult);
235  return newResult;
236  }
237 }

◆ A0_Mz2_mh2()

double EWSMcache::A0_Mz2_mh2 ( ) const

A cache method.

Returns
\(A_0(m_h^2)\) with \(\mu=M_Z\)

Definition at line 211 of file EWSMcache.cpp.

212 {
213  int NumPar = 2;
214  double params[] = {SM.getMz(), SM.getMHl()};
215 
216  if (CacheCheck(A0_Mz2_mh2_cache, NumPar, params))
217  return A0_Mz2_mh2_cache[NumPar];
218  else {
219  double newResult = PV.A0(SM.getMz() * SM.getMz(), SM.getMHl() * SM.getMHl());
220  newCacheForDouble(A0_Mz2_mh2_cache, NumPar, params, newResult);
221  return newResult;
222  }
223 }

◆ A0_Mz2_Mw2()

double EWSMcache::A0_Mz2_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(A_0(M_W^2)\) with \(\mu=M_Z\)

Definition at line 197 of file EWSMcache.cpp.

198 {
199  int NumPar = 2;
200  double params[] = {SM.getMz(), Mw_i};
201 
202  if (CacheCheck(A0_Mz2_Mw2_cache, NumPar, params))
203  return A0_Mz2_Mw2_cache[NumPar];
204  else {
205  double newResult = PV.A0(SM.getMz() * SM.getMz(), Mw_i * Mw_i);
206  newCacheForDouble(A0_Mz2_Mw2_cache, NumPar, params, newResult);
207  return newResult;
208  }
209 }

◆ A0_Mz2_Mz2()

double EWSMcache::A0_Mz2_Mz2 ( ) const

A cache method.

Returns
\(A_0(M_Z^2)\) with \(\mu=M_Z\)

Definition at line 253 of file EWSMcache.cpp.

254 {
255  int NumPar = 1;
256  double params[] = {SM.getMz()};
257 
258  if (CacheCheck(A0_Mz2_Mz2_cache, NumPar, params))
259  return A0_Mz2_Mz2_cache[NumPar];
260  else {
261  double newResult = PV.A0(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz());
262  newCacheForDouble(A0_Mz2_Mz2_cache, NumPar, params, newResult);
263  return newResult;
264  }
265 }

◆ a_f()

double EWSMcache::a_f ( const Particle  f) const
inline

The tree-level axial-vector coupling for \(Z\to f\bar{f}\), denoted as \(a_f\).

Parameters
[in]fa lepton or quark
Returns
\(a_f\)

Definition at line 301 of file EWSMcache.h.

302  {
303  return ( f.getIsospin());
304  }

◆ Als()

double EWSMcache::Als ( const double  mu,
const orders  order 
) const
inline

The strong coupling \(\alpha_s(\mu)\).

Parameters
[in]murenormalization scale \(\mu\) in GeV
[in]orderorder in QCD (see orders)
Returns
\(\alpha_s(\mu)\) at the given order

Definition at line 366 of file EWSMcache.h.

367  {
368  return ( SM.Als(mu, order));
369  }

◆ alsMt()

double EWSMcache::alsMt ( ) const
inline

The strong coupling \(\alpha_s(m_t^2)\) at NNLO.

Returns
\(\alpha_s(m_t^2)\) at NNLO
Attention
The constant value \(\alpha_s(m_t^2)=0.1074432788\) is returned when the flag FlagDebug is set to true.

Definition at line 378 of file EWSMcache.h.

379  {
380  if (FlagDebug)
381  return ( 0.1074432788); // for debug
382  else
383  return ( SM.Als(SM.getMtpole(), FULLNNLO));
384  }

◆ B0_Mw2_Mw2_0_Mw2()

gslpp::complex EWSMcache::B0_Mw2_Mw2_0_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_0(M_W^2;0,M_W^2)\) with \(\mu=M_W\)

Definition at line 461 of file EWSMcache.cpp.

462 {
463  int NumPar = 1;
464  double params[] = {Mw_i};
465 
466  if (CacheCheck(B0_Mw2_Mw2_0_Mw2_cache, NumPar, params))
467  return gslpp::complex(B0_Mw2_Mw2_0_Mw2_cache[NumPar],
468  B0_Mw2_Mw2_0_Mw2_cache[NumPar + 1], false);
469  else {
470  gslpp::complex newResult = PV.B0(Mw_i*Mw_i, Mw_i*Mw_i, 0.0, Mw_i * Mw_i);
471  newCacheForComplex(B0_Mw2_Mw2_0_Mw2_cache, NumPar, params, newResult);
472  return newResult;
473  }
474 }

◆ B0_Mw2_Mw2_mh2_Mw2()

gslpp::complex EWSMcache::B0_Mw2_Mw2_mh2_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_0(M_W^2;m_h^2,M_W^2)\) with \(\mu=M_W\)

Definition at line 446 of file EWSMcache.cpp.

447 {
448  int NumPar = 2;
449  double params[] = {Mw_i, SM.getMHl()};
450 
451  if (CacheCheck(B0_Mw2_Mw2_mh2_Mw2_cache, NumPar, params))
452  return gslpp::complex(B0_Mw2_Mw2_mh2_Mw2_cache[NumPar],
453  B0_Mw2_Mw2_mh2_Mw2_cache[NumPar + 1], false);
454  else {
455  gslpp::complex newResult = PV.B0(Mw_i*Mw_i, Mw_i*Mw_i, SM.getMHl() * SM.getMHl(), Mw_i * Mw_i);
456  newCacheForComplex(B0_Mw2_Mw2_mh2_Mw2_cache, NumPar, params, newResult);
457  return newResult;
458  }
459 }

◆ B0_Mw2_Mw2_Mz2_Mw2()

gslpp::complex EWSMcache::B0_Mw2_Mw2_Mz2_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_0(M_W^2;M_Z^2,M_W^2)\) with \(\mu=M_W\)

Definition at line 431 of file EWSMcache.cpp.

432 {
433  int NumPar = 2;
434  double params[] = {Mw_i, SM.getMz()};
435 
436  if (CacheCheck(B0_Mw2_Mw2_Mz2_Mw2_cache, NumPar, params))
437  return gslpp::complex(B0_Mw2_Mw2_Mz2_Mw2_cache[NumPar],
438  B0_Mw2_Mw2_Mz2_Mw2_cache[NumPar + 1], false);
439  else {
440  gslpp::complex newResult = PV.B0(Mw_i*Mw_i, Mw_i*Mw_i, SM.getMz() * SM.getMz(), Mw_i * Mw_i);
441  newCacheForComplex(B0_Mw2_Mw2_Mz2_Mw2_cache, NumPar, params, newResult);
442  return newResult;
443  }
444 }

◆ B0_Mw2_Mz2_Mt2_Mt2()

gslpp::complex EWSMcache::B0_Mw2_Mz2_Mt2_Mt2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_0(M_Z^2;m_t^2,m_t^2)\) with \(\mu=M_W\)

Definition at line 311 of file EWSMcache.cpp.

312 {
313  int NumPar = 3;
314  double params[] = {Mw_i, SM.getMz(), SM.getMtpole()};
315 
316  if (CacheCheck(B0_Mw2_Mz2_Mt2_Mt2_cache, NumPar, params))
317  return gslpp::complex(B0_Mw2_Mz2_Mt2_Mt2_cache[NumPar],
318  B0_Mw2_Mz2_Mt2_Mt2_cache[NumPar + 1], false);
319  else {
320  gslpp::complex newResult = PV.B0(Mw_i*Mw_i, SM.getMz() * SM.getMz(), SM.getMtpole() * SM.getMtpole(), SM.getMtpole() * SM.getMtpole());
321  newCacheForComplex(B0_Mw2_Mz2_Mt2_Mt2_cache, NumPar, params, newResult);
322  return newResult;
323  }
324 }

◆ B0_Mw2_Mz2_Mw2_Mw2()

gslpp::complex EWSMcache::B0_Mw2_Mz2_Mw2_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_0(M_Z^2;M_W^2,M_W^2)\) with \(\mu=M_W\)

Definition at line 416 of file EWSMcache.cpp.

417 {
418  int NumPar = 2;
419  double params[] = {Mw_i, SM.getMz()};
420 
421  if (CacheCheck(B0_Mw2_Mz2_Mw2_Mw2_cache, NumPar, params))
422  return gslpp::complex(B0_Mw2_Mz2_Mw2_Mw2_cache[NumPar],
423  B0_Mw2_Mz2_Mw2_Mw2_cache[NumPar + 1], false);
424  else {
425  gslpp::complex newResult = PV.B0(Mw_i*Mw_i, SM.getMz() * SM.getMz(), Mw_i*Mw_i, Mw_i * Mw_i);
426  newCacheForComplex(B0_Mw2_Mz2_Mw2_Mw2_cache, NumPar, params, newResult);
427  return newResult;
428  }
429 }

◆ B0_Mz2_0_0_Mw2()

gslpp::complex EWSMcache::B0_Mz2_0_0_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_0(0;0,M_W^2)\) with \(\mu=M_Z\)

Definition at line 401 of file EWSMcache.cpp.

402 {
403  int NumPar = 2;
404  double params[] = {SM.getMz(), Mw_i};
405 
406  if (CacheCheck(B0_Mz2_0_0_Mw2_cache, NumPar, params))
407  return gslpp::complex(B0_Mz2_0_0_Mw2_cache[NumPar],
408  B0_Mz2_0_0_Mw2_cache[NumPar + 1], false);
409  else {
410  gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), 0.0, 0.0, Mw_i * Mw_i);
411  newCacheForComplex(B0_Mz2_0_0_Mw2_cache, NumPar, params, newResult);
412  return newResult;
413  }
414 }

◆ B0_Mz2_0_mh2_Mw2()

gslpp::complex EWSMcache::B0_Mz2_0_mh2_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_0(0;m_h^2,M_W^2)\) with \(\mu=M_Z\)

Definition at line 296 of file EWSMcache.cpp.

297 {
298  int NumPar = 3;
299  double params[] = {SM.getMz(), SM.getMHl(), Mw_i};
300 
301  if (CacheCheck(B0_Mz2_0_mh2_Mw2_cache, NumPar, params))
302  return gslpp::complex(B0_Mz2_0_mh2_Mw2_cache[NumPar],
303  B0_Mz2_0_mh2_Mw2_cache[NumPar + 1], false);
304  else {
305  gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), 0.0, SM.getMHl() * SM.getMHl(), Mw_i * Mw_i);
306  newCacheForComplex(B0_Mz2_0_mh2_Mw2_cache, NumPar, params, newResult);
307  return newResult;
308  }
309 }

◆ B0_Mz2_0_Mz2_Mw2()

gslpp::complex EWSMcache::B0_Mz2_0_Mz2_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_0(0;M_Z^2,M_W^2)\) with \(\mu=M_Z\)

Definition at line 386 of file EWSMcache.cpp.

387 {
388  int NumPar = 2;
389  double params[] = {SM.getMz(), Mw_i};
390 
391  if (CacheCheck(B0_Mz2_0_Mz2_Mw2_cache, NumPar, params))
392  return gslpp::complex(B0_Mz2_0_Mz2_Mw2_cache[NumPar],
393  B0_Mz2_0_Mz2_Mw2_cache[NumPar + 1], false);
394  else {
395  gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), 0.0, SM.getMz() * SM.getMz(), Mw_i * Mw_i);
396  newCacheForComplex(B0_Mz2_0_Mz2_Mw2_cache, NumPar, params, newResult);
397  return newResult;
398  }
399 }

◆ B0_Mz2_Mw2_0_Mw2()

gslpp::complex EWSMcache::B0_Mz2_Mw2_0_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_0(M_W^2;0,M_W^2)\) with \(\mu=M_Z\)

Definition at line 371 of file EWSMcache.cpp.

372 {
373  int NumPar = 2;
374  double params[] = {SM.getMz(), Mw_i};
375 
376  if (CacheCheck(B0_Mz2_Mw2_0_Mw2_cache, NumPar, params))
377  return gslpp::complex(B0_Mz2_Mw2_0_Mw2_cache[NumPar],
378  B0_Mz2_Mw2_0_Mw2_cache[NumPar + 1], false);
379  else {
380  gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), Mw_i*Mw_i, 0.0, Mw_i * Mw_i);
381  newCacheForComplex(B0_Mz2_Mw2_0_Mw2_cache, NumPar, params, newResult);
382  return newResult;
383  }
384 }

◆ B0_Mz2_Mw2_mh2_Mw2()

gslpp::complex EWSMcache::B0_Mz2_Mw2_mh2_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_0(M_W^2;m_h^2,M_W^2)\) with \(\mu=M_Z\)

Definition at line 281 of file EWSMcache.cpp.

282 {
283  int NumPar = 3;
284  double params[] = {SM.getMz(), Mw_i, SM.getMHl()};
285 
286  if (CacheCheck(B0_Mz2_Mw2_mh2_Mw2_cache, NumPar, params))
287  return gslpp::complex(B0_Mz2_Mw2_mh2_Mw2_cache[NumPar],
288  B0_Mz2_Mw2_mh2_Mw2_cache[NumPar + 1], false);
289  else {
290  gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), Mw_i*Mw_i, SM.getMHl() * SM.getMHl(), Mw_i * Mw_i);
291  newCacheForComplex(B0_Mz2_Mw2_mh2_Mw2_cache, NumPar, params, newResult);
292  return newResult;
293  }
294 }

◆ B0_Mz2_Mw2_Mz2_Mw2()

gslpp::complex EWSMcache::B0_Mz2_Mw2_Mz2_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_0(M_W^2;M_Z^2,M_W^2)\) with \(\mu=M_Z\)

Definition at line 356 of file EWSMcache.cpp.

357 {
358  int NumPar = 2;
359  double params[] = {SM.getMz(), Mw_i};
360 
361  if (CacheCheck(B0_Mz2_Mw2_Mz2_Mw2_cache, NumPar, params))
362  return gslpp::complex(B0_Mz2_Mw2_Mz2_Mw2_cache[NumPar],
363  B0_Mz2_Mw2_Mz2_Mw2_cache[NumPar + 1], false);
364  else {
365  gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), Mw_i*Mw_i, SM.getMz() * SM.getMz(), Mw_i * Mw_i);
366  newCacheForComplex(B0_Mz2_Mw2_Mz2_Mw2_cache, NumPar, params, newResult);
367  return newResult;
368  }
369 }

◆ B0_Mz2_Mz2_mf2_mf2()

gslpp::complex EWSMcache::B0_Mz2_Mz2_mf2_mf2 ( const Particle  f) const

A cache method.

Parameters
[in]fa lepton or quark
Returns
\(B_0(M_Z^2;m_f^2,m_f^2)\) with \(\mu=M_Z\)

Definition at line 476 of file EWSMcache.cpp.

477 {
478  int NumPar = 2;
479  double params[] = {SM.getMz(), mf(f, SM.getMz())};
480  int ind = f.getIndex();
481  if (CacheCheck(B0_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params))
482  return gslpp::complex(B0_Mz2_Mz2_mf2_mf2_cache[ind][NumPar],
483  B0_Mz2_Mz2_mf2_mf2_cache[ind][NumPar + 1], false);
484  else {
485  gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), mf2(f, SM.getMz()), mf2(f, SM.getMz()));
486  newCacheForComplex(B0_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params, newResult);
487  return newResult;
488  }
489 }

◆ B0_Mz2_Mz2_mh2_Mz2()

gslpp::complex EWSMcache::B0_Mz2_Mz2_mh2_Mz2 ( ) const

A cache method.

Returns
\(B_0(M_Z^2;m_h^2,M_Z^2)\) with \(\mu=M_Z\)

Definition at line 341 of file EWSMcache.cpp.

342 {
343  int NumPar = 2;
344  double params[] = {SM.getMz(), SM.getMHl()};
345 
346  if (CacheCheck(B0_Mz2_Mz2_mh2_Mz2_cache, NumPar, params))
347  return gslpp::complex(B0_Mz2_Mz2_mh2_Mz2_cache[NumPar],
348  B0_Mz2_Mz2_mh2_Mz2_cache[NumPar + 1], false);
349  else {
350  gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), SM.getMHl() * SM.getMHl(), SM.getMz() * SM.getMz());
351  newCacheForComplex(B0_Mz2_Mz2_mh2_Mz2_cache, NumPar, params, newResult);
352  return newResult;
353  }
354 }

◆ B0_Mz2_Mz2_Mw2_Mw2()

gslpp::complex EWSMcache::B0_Mz2_Mz2_Mw2_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_0(M_Z^2;M_W^2,M_W^2)\) with \(\mu=M_Z\)

Definition at line 326 of file EWSMcache.cpp.

327 {
328  int NumPar = 2;
329  double params[] = {SM.getMz(), Mw_i};
330 
331  if (CacheCheck(B0_Mz2_Mz2_Mw2_Mw2_cache, NumPar, params))
332  return gslpp::complex(B0_Mz2_Mz2_Mw2_Mw2_cache[NumPar],
333  B0_Mz2_Mz2_Mw2_Mw2_cache[NumPar + 1], false);
334  else {
335  gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), Mw_i*Mw_i, Mw_i * Mw_i);
336  newCacheForComplex(B0_Mz2_Mz2_Mw2_Mw2_cache, NumPar, params, newResult);
337  return newResult;
338  }
339 }

◆ B0p_Mw2_Mw2_0_Mw2()

gslpp::complex EWSMcache::B0p_Mw2_Mw2_0_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_{0p}(M_W^2;0,M_W^2)\) with \(\mu=M_W\)

Definition at line 581 of file EWSMcache.cpp.

582 {
583  int NumPar = 1;
584  double params[] = {Mw_i};
585 
586  if (CacheCheck(B0p_Mw2_Mw2_0_Mw2_cache, NumPar, params))
587  return gslpp::complex(B0p_Mw2_Mw2_0_Mw2_cache[NumPar],
588  B0p_Mw2_Mw2_0_Mw2_cache[NumPar + 1], false);
589  else {
590  gslpp::complex newResult = PV.B0p(Mw_i*Mw_i, Mw_i*Mw_i, 0.0, Mw_i * Mw_i);
591  newCacheForComplex(B0p_Mw2_Mw2_0_Mw2_cache, NumPar, params, newResult);
592  return newResult;
593  }
594 }

◆ B0p_Mw2_Mw2_mh2_Mw2()

gslpp::complex EWSMcache::B0p_Mw2_Mw2_mh2_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_{0p}(M_W^2;m_h^2,M_W^2)\) with \(\mu=M_W\)

Definition at line 566 of file EWSMcache.cpp.

567 {
568  int NumPar = 2;
569  double params[] = {Mw_i, SM.getMHl()};
570 
571  if (CacheCheck(B0p_Mw2_Mw2_mh2_Mw2_cache, NumPar, params))
572  return gslpp::complex(B0p_Mw2_Mw2_mh2_Mw2_cache[NumPar],
573  B0p_Mw2_Mw2_mh2_Mw2_cache[NumPar + 1], false);
574  else {
575  gslpp::complex newResult = PV.B0p(Mw_i*Mw_i, Mw_i*Mw_i, SM.getMHl() * SM.getMHl(), Mw_i * Mw_i);
576  newCacheForComplex(B0p_Mw2_Mw2_mh2_Mw2_cache, NumPar, params, newResult);
577  return newResult;
578  }
579 }

◆ B0p_Mw2_Mw2_Mz2_Mw2()

gslpp::complex EWSMcache::B0p_Mw2_Mw2_Mz2_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_{0p}(M_W^2;M_Z^2,M_W^2)\) with \(\mu=M_W\)

Definition at line 551 of file EWSMcache.cpp.

552 {
553  int NumPar = 2;
554  double params[] = {Mw_i, SM.getMz()};
555 
556  if (CacheCheck(B0p_Mw2_Mw2_Mz2_Mw2_cache, NumPar, params))
557  return gslpp::complex(B0p_Mw2_Mw2_Mz2_Mw2_cache[NumPar],
558  B0p_Mw2_Mw2_Mz2_Mw2_cache[NumPar + 1], false);
559  else {
560  gslpp::complex newResult = PV.B0p(Mw_i*Mw_i, Mw_i*Mw_i, SM.getMz() * SM.getMz(), Mw_i * Mw_i);
561  newCacheForComplex(B0p_Mw2_Mw2_Mz2_Mw2_cache, NumPar, params, newResult);
562  return newResult;
563  }
564 }

◆ B0p_Mz2_0_mh2_Mw2()

gslpp::complex EWSMcache::B0p_Mz2_0_mh2_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_{0p}(0;m_h^2,M_W^2)\) with \(\mu=M_Z\)

Definition at line 491 of file EWSMcache.cpp.

492 {
493  int NumPar = 3;
494  double params[] = {SM.getMz(), SM.getMHl(), Mw_i};
495 
496  if (CacheCheck(B0p_Mz2_0_mh2_Mw2_cache, NumPar, params))
497  return gslpp::complex(B0p_Mz2_0_mh2_Mw2_cache[NumPar],
498  B0p_Mz2_0_mh2_Mw2_cache[NumPar + 1], false);
499  else {
500  gslpp::complex newResult = PV.B0p(SM.getMz() * SM.getMz(), 0.0, SM.getMHl() * SM.getMHl(), Mw_i * Mw_i);
501  newCacheForComplex(B0p_Mz2_0_mh2_Mw2_cache, NumPar, params, newResult);
502  return newResult;
503  }
504 }

◆ B0p_Mz2_0_Mz2_Mw2()

gslpp::complex EWSMcache::B0p_Mz2_0_Mz2_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_{0p}(0;M_Z^2,M_W^2)\) with \(\mu=M_Z\)

Definition at line 521 of file EWSMcache.cpp.

522 {
523  int NumPar = 2;
524  double params[] = {SM.getMz(), Mw_i};
525 
526  if (CacheCheck(B0p_Mz2_0_Mz2_Mw2_cache, NumPar, params))
527  return gslpp::complex(B0p_Mz2_0_Mz2_Mw2_cache[NumPar],
528  B0p_Mz2_0_Mz2_Mw2_cache[NumPar + 1], false);
529  else {
530  gslpp::complex newResult = PV.B0p(SM.getMz() * SM.getMz(), 0.0, SM.getMz() * SM.getMz(), Mw_i * Mw_i);
531  newCacheForComplex(B0p_Mz2_0_Mz2_Mw2_cache, NumPar, params, newResult);
532  return newResult;
533  }
534 }

◆ B0p_Mz2_Mz2_mf2_mf2()

gslpp::complex EWSMcache::B0p_Mz2_Mz2_mf2_mf2 ( const Particle  f) const

A cache method.

Parameters
[in]fa lepton or quark
Returns
\(B_{0p}(M_Z^2;m_f^2,m_f^2)\) with \(\mu=M_Z\)

Definition at line 596 of file EWSMcache.cpp.

597 {
598  int NumPar = 2;
599  double params[] = {SM.getMz(), mf(f, SM.getMz())};
600  int ind = f.getIndex();
601  if (CacheCheck(B0p_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params))
602  return gslpp::complex(B0p_Mz2_Mz2_mf2_mf2_cache[ind][NumPar],
603  B0p_Mz2_Mz2_mf2_mf2_cache[ind][NumPar + 1], false);
604  else {
605  gslpp::complex newResult = PV.B0p(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), mf2(f, SM.getMz()), mf2(f, SM.getMz()));
606  newCacheForComplex(B0p_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params, newResult);
607  return newResult;
608  }
609 }

◆ B0p_Mz2_Mz2_mh2_Mz2()

gslpp::complex EWSMcache::B0p_Mz2_Mz2_mh2_Mz2 ( ) const

A cache method.

Returns
\(B_{0p}(M_Z^2;m_h^2,M_Z^2)\) with \(\mu=M_Z\)

Definition at line 506 of file EWSMcache.cpp.

507 {
508  int NumPar = 2;
509  double params[] = {SM.getMz(), SM.getMHl()};
510 
511  if (CacheCheck(B0p_Mz2_Mz2_mh2_Mz2_cache, NumPar, params))
512  return gslpp::complex(B0p_Mz2_Mz2_mh2_Mz2_cache[NumPar],
513  B0p_Mz2_Mz2_mh2_Mz2_cache[NumPar + 1], false);
514  else {
515  gslpp::complex newResult = PV.B0p(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), SM.getMHl() * SM.getMHl(), SM.getMz() * SM.getMz());
516  newCacheForComplex(B0p_Mz2_Mz2_mh2_Mz2_cache, NumPar, params, newResult);
517  return newResult;
518  }
519 }

◆ B0p_Mz2_Mz2_Mw2_Mw2()

gslpp::complex EWSMcache::B0p_Mz2_Mz2_Mw2_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_{0p}(M_Z^2;M_W^2,M_W^2)\) with \(\mu=M_Z\)

Definition at line 536 of file EWSMcache.cpp.

537 {
538  int NumPar = 2;
539  double params[] = {SM.getMz(), Mw_i};
540 
541  if (CacheCheck(B0p_Mz2_Mz2_Mw2_Mw2_cache, NumPar, params))
542  return gslpp::complex(B0p_Mz2_Mz2_Mw2_Mw2_cache[NumPar],
543  B0p_Mz2_Mz2_Mw2_Mw2_cache[NumPar + 1], false);
544  else {
545  gslpp::complex newResult = PV.B0p(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), Mw_i*Mw_i, Mw_i * Mw_i);
546  newCacheForComplex(B0p_Mz2_Mz2_Mw2_Mw2_cache, NumPar, params, newResult);
547  return newResult;
548  }
549 }

◆ B1_Mz2_0_mf2_mfprime2()

gslpp::complex EWSMcache::B1_Mz2_0_mf2_mfprime2 ( const int  gen) const

A cache method.

Parameters
[in]genthe generation index of a lepton doublet
Returns
\(B_1(0;m_l^2,m_{l'}^2)\) with \(\mu=M_Z\)

Definition at line 611 of file EWSMcache.cpp.

612 {
613  int NumPar = 3;
614  double mymf, mymfprime;
615  double Mz = SM.getMz();
616  if (gen < 3) {
617  mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
618  mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
619  } else {
620  int genq = gen - 3;
621  mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
622  mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
623  }
624  double params[] = {Mz, mymf, mymfprime};
625 
626  if (CacheCheck(B1_Mz2_0_mf2_mfprime2_cache[gen], NumPar, params))
627  return gslpp::complex(B1_Mz2_0_mf2_mfprime2_cache[gen][NumPar],
628  B1_Mz2_0_mf2_mfprime2_cache[gen][NumPar + 1], false);
629  else {
630  double mf2 = mymf*mymf;
631  double mfprime2 = mymfprime*mymfprime;
632  gslpp::complex newResult = PV.B1(Mz * Mz, 0.0, mf2, mfprime2);
633  newCacheForComplex(B1_Mz2_0_mf2_mfprime2_cache[gen], NumPar, params, newResult);
634  return newResult;
635  }
636 }

◆ B1_Mz2_0_mfprime2_mf2()

gslpp::complex EWSMcache::B1_Mz2_0_mfprime2_mf2 ( const int  gen) const

A cache method.

Parameters
[in]genthe generation index of a lepton doublet
Returns
\(B_1(0;m_{l'}^2,m_l^2)\) with \(\mu=M_Z\)

Definition at line 638 of file EWSMcache.cpp.

639 {
640  int NumPar = 3;
641  double mymf, mymfprime;
642  double Mz = SM.getMz();
643  if (gen < 3) {
644  mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
645  mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
646  } else {
647  int genq = gen - 3;
648  mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
649  mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
650  }
651  double params[] = {Mz, mymf, mymfprime};
652 
653  if (CacheCheck(B1_Mz2_0_mfprime2_mf2_cache[gen], NumPar, params))
654  return gslpp::complex(B1_Mz2_0_mfprime2_mf2_cache[gen][NumPar],
655  B1_Mz2_0_mfprime2_mf2_cache[gen][NumPar + 1], false);
656  else {
657  double mf2 = mymf*mymf;
658  double mfprime2 = mymfprime*mymfprime;
659  gslpp::complex newResult = PV.B1(Mz * Mz, 0.0, mfprime2, mf2);
660  newCacheForComplex(B1_Mz2_0_mfprime2_mf2_cache[gen], NumPar, params, newResult);
661  return newResult;
662  }
663 }

◆ B1_Mz2_Mw2_mf2_mfprime2()

gslpp::complex EWSMcache::B1_Mz2_Mw2_mf2_mfprime2 ( const int  gen,
const double  Mw_i 
) const

A cache method.

Parameters
[in]genthe generation index of a lepton doublet
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_1(M_W^2;m_l^2,m_{l'}^2)\) with \(\mu=M_Z\)

Definition at line 665 of file EWSMcache.cpp.

666 {
667  int NumPar = 4;
668  double mymf, mymfprime;
669  double Mz = SM.getMz();
670  if (gen < 3) {
671  mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
672  mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
673  } else {
674  int genq = gen - 3;
675  mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
676  mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
677  }
678  double params[] = {Mz, Mw_i, mymf, mymfprime};
679 
680  if (CacheCheck(B1_Mz2_Mw2_mf2_mfprime2_cache[gen], NumPar, params))
681  return gslpp::complex(B1_Mz2_Mw2_mf2_mfprime2_cache[gen][NumPar],
682  B1_Mz2_Mw2_mf2_mfprime2_cache[gen][NumPar + 1], false);
683  else {
684  double mf2 = mymf*mymf;
685  double mfprime2 = mymfprime*mymfprime;
686  gslpp::complex newResult = PV.B1(Mz * Mz, Mw_i*Mw_i, mf2, mfprime2);
687  newCacheForComplex(B1_Mz2_Mw2_mf2_mfprime2_cache[gen], NumPar, params, newResult);
688  return newResult;
689  }
690 }

◆ B1_Mz2_Mw2_mfprime2_mf2()

gslpp::complex EWSMcache::B1_Mz2_Mw2_mfprime2_mf2 ( const int  gen,
const double  Mw_i 
) const

A cache method.

Parameters
[in]genthe generation index of a lepton doublet
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_1(M_W^2;m_{l'}^2,m_l^2)\) with \(\mu=M_Z\)

Definition at line 692 of file EWSMcache.cpp.

693 {
694  int NumPar = 4;
695  double mymf, mymfprime;
696  double Mz = SM.getMz();
697  if (gen < 3) {
698  mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
699  mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
700  } else {
701  int genq = gen - 3;
702  mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
703  mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
704  }
705  double params[] = {Mz, Mw_i, mymf, mymfprime};
706 
707  if (CacheCheck(B1_Mz2_Mw2_mfprime2_mf2_cache[gen], NumPar, params))
708  return gslpp::complex(B1_Mz2_Mw2_mfprime2_mf2_cache[gen][NumPar],
709  B1_Mz2_Mw2_mfprime2_mf2_cache[gen][NumPar + 1], false);
710  else {
711  double mf2 = mymf*mymf;
712  double mfprime2 = mymfprime*mymfprime;
713  gslpp::complex newResult = PV.B1(Mz * Mz, Mw_i*Mw_i, mfprime2, mf2);
714  newCacheForComplex(B1_Mz2_Mw2_mfprime2_mf2_cache[gen], NumPar, params, newResult);
715  return newResult;
716  }
717 }

◆ B1p_Mw2_Mw2_mf2_mfprime2()

gslpp::complex EWSMcache::B1p_Mw2_Mw2_mf2_mfprime2 ( const int  gen,
const double  Mw_i 
) const

A cache method.

Parameters
[in]genthe generation index of a lepton doublet
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_{1p}(M_W^2;m_l^2,m_{l'}^2)\) with \(\mu=M_W\)

Definition at line 719 of file EWSMcache.cpp.

720 {
721  int NumPar = 3;
722  double mymf, mymfprime;
723  double Mz = SM.getMz();
724  if (gen < 3) {
725  mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
726  mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
727  } else {
728  int genq = gen - 3;
729  mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
730  mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
731  }
732  double params[] = {Mw_i, mymf, mymfprime};
733 
734  if (CacheCheck(B1p_Mw2_Mw2_mf2_mfprime2_cache[gen], NumPar, params))
735  return gslpp::complex(B1p_Mw2_Mw2_mf2_mfprime2_cache[gen][NumPar],
736  B1p_Mw2_Mw2_mf2_mfprime2_cache[gen][NumPar + 1], false);
737  else {
738  double mf2 = mymf*mymf;
739  double mfprime2 = mymfprime*mymfprime;
740  gslpp::complex newResult = PV.B1p(Mw_i*Mw_i, Mw_i*Mw_i, mf2, mfprime2);
741  newCacheForComplex(B1p_Mw2_Mw2_mf2_mfprime2_cache[gen], NumPar, params, newResult);
742  return newResult;
743  }
744 }

◆ B1p_Mw2_Mw2_mfprime2_mf2()

gslpp::complex EWSMcache::B1p_Mw2_Mw2_mfprime2_mf2 ( const int  gen,
const double  Mw_i 
) const

A cache method.

Parameters
[in]genthe generation index of a lepton doublet
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_{1p}(M_W^2;m_{l'}^2,m_l^2)\) with \(\mu=M_W\)

Definition at line 746 of file EWSMcache.cpp.

747 {
748  int NumPar = 3;
749  double mymf, mymfprime;
750  double Mz = SM.getMz();
751  if (gen < 3) {
752  mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
753  mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
754  } else {
755  int genq = gen - 3;
756  mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
757  mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
758  }
759  double params[] = {Mw_i, mymf, mymfprime};
760 
761  if (CacheCheck(B1p_Mw2_Mw2_mfprime2_mf2_cache[gen], NumPar, params))
762  return gslpp::complex(B1p_Mw2_Mw2_mfprime2_mf2_cache[gen][NumPar],
763  B1p_Mw2_Mw2_mfprime2_mf2_cache[gen][NumPar + 1], false);
764  else {
765  double mf2 = mymf*mymf;
766  double mfprime2 = mymfprime*mymfprime;
767  gslpp::complex newResult = PV.B1p(Mw_i*Mw_i, Mw_i*Mw_i, mfprime2, mf2);
768  newCacheForComplex(B1p_Mw2_Mw2_mfprime2_mf2_cache[gen], NumPar, params, newResult);
769  return newResult;
770  }
771 }

◆ Bf_Mw2_Mw2_mfprime2_mf2()

gslpp::complex EWSMcache::Bf_Mw2_Mw2_mfprime2_mf2 ( const int  gen,
const double  Mw_i 
) const

A cache method.

Parameters
[in]genthe generation index of a lepton doublet
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_f(M_W^2;m_{l'}^2,m_l^2)\) with \(\mu=M_W\)

Definition at line 858 of file EWSMcache.cpp.

859 {
860  int NumPar = 3;
861  double mymf, mymfprime;
862  double Mz = SM.getMz();
863  if (gen < 3) {
864  mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
865  mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
866  } else {
867  int genq = gen - 3;
868  mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
869  mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
870  }
871  double params[] = {Mw_i, mymf, mymfprime};
872 
873  if (CacheCheck(Bf_Mw2_Mw2_mfprime2_mf2_cache[gen], NumPar, params))
874  return gslpp::complex(Bf_Mw2_Mw2_mfprime2_mf2_cache[gen][NumPar],
875  Bf_Mw2_Mw2_mfprime2_mf2_cache[gen][NumPar + 1], false);
876  else {
877  double mf2 = mymf*mymf;
878  double mfprime2 = mymfprime*mymfprime;
879  gslpp::complex newResult = PV.Bf(Mw_i*Mw_i, Mw_i*Mw_i, mfprime2, mf2);
880  newCacheForComplex(Bf_Mw2_Mw2_mfprime2_mf2_cache[gen], NumPar, params, newResult);
881  return newResult;
882  }
883 }

◆ Bf_Mz2_0_mf2_mf2()

gslpp::complex EWSMcache::Bf_Mz2_0_mf2_mf2 ( const Particle  f) const

A cache method.

Parameters
[in]fa lepton or quark
Returns
\(B_f(0;m_f^2,m_f^2)\) with \(\mu=M_Z\)

Definition at line 788 of file EWSMcache.cpp.

789 {
790  int NumPar = 2;
791  double params[] = {SM.getMz(), mf(f, SM.getMz())};
792  if (params[1] == 0.0)
793  throw std::runtime_error("Error in EWSMcache::Bf_Mz_0_mf_mf()");
794  int ind = f.getIndex();
795  if (CacheCheck(Bf_Mz2_0_mf2_mf2_cache[ind], NumPar, params))
796  return gslpp::complex(Bf_Mz2_0_mf2_mf2_cache[ind][NumPar],
797  Bf_Mz2_0_mf2_mf2_cache[ind][NumPar + 1], false);
798  else {
799  gslpp::complex newResult = PV.Bf(SM.getMz() * SM.getMz(), 0.0, mf2(f, SM.getMz()), mf2(f, SM.getMz()));
800  newCacheForComplex(Bf_Mz2_0_mf2_mf2_cache[ind], NumPar, params, newResult);
801  return newResult;
802  }
803 }

◆ Bf_Mz2_0_mfprime2_mf2()

gslpp::complex EWSMcache::Bf_Mz2_0_mfprime2_mf2 ( const int  gen) const

A cache method.

Parameters
[in]genthe generation index of a lepton doublet
Returns
\(B_f(0;m_{l'}^2,m_l^2)\) with \(\mu=M_Z\)

Definition at line 831 of file EWSMcache.cpp.

832 {
833  int NumPar = 3;
834  double mymf, mymfprime;
835  double Mz = SM.getMz();
836  if (gen < 3) {
837  mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
838  mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
839  } else {
840  int genq = gen - 3;
841  mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
842  mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
843  }
844  double params[] = {Mz, mymf, mymfprime};
845 
846  if (CacheCheck(Bf_Mz2_0_mfprime2_mf2_cache[gen], NumPar, params))
847  return gslpp::complex(Bf_Mz2_0_mfprime2_mf2_cache[gen][NumPar],
848  Bf_Mz2_0_mfprime2_mf2_cache[gen][NumPar + 1], false);
849  else {
850  double mf2 = mymf*mymf;
851  double mfprime2 = mymfprime*mymfprime;
852  gslpp::complex newResult = PV.Bf(Mz * Mz, 0.0, mfprime2, mf2);
853  newCacheForComplex(Bf_Mz2_0_mfprime2_mf2_cache[gen], NumPar, params, newResult);
854  return newResult;
855  }
856 }

◆ Bf_Mz2_Mw2_mfprime2_mf2()

gslpp::complex EWSMcache::Bf_Mz2_Mw2_mfprime2_mf2 ( const int  gen,
const double  Mw_i 
) const

A cache method.

Parameters
[in]genthe generation index of a lepton doublet
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_f(M_W^2;m_{l'}^2,m_l^2)\) with \(\mu=M_Z\)

Definition at line 805 of file EWSMcache.cpp.

806 {
807  int NumPar = 4;
808  double mymf, mymfprime;
809  double Mz = SM.getMz();
810  if (gen < 3) {
811  mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
812  mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
813  } else {
814  int genq = gen - 3;
815  mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
816  mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
817  }
818  double params[] = {Mz, Mw_i, mymf, mymfprime};
819  if (CacheCheck(Bf_Mz2_Mw2_mfprime2_mf2_cache[gen], NumPar, params))
820  return gslpp::complex(Bf_Mz2_Mw2_mfprime2_mf2_cache[gen][NumPar],
821  Bf_Mz2_Mw2_mfprime2_mf2_cache[gen][NumPar + 1], false);
822  else {
823  double mf2 = mymf*mymf;
824  double mfprime2 = mymfprime*mymfprime;
825  gslpp::complex newResult = PV.Bf(Mz * Mz, Mw_i*Mw_i, mfprime2, mf2);
826  newCacheForComplex(Bf_Mz2_Mw2_mfprime2_mf2_cache[gen], NumPar, params, newResult);
827  return newResult;
828  }
829 }

◆ Bf_Mz2_Mz2_mf2_mf2()

gslpp::complex EWSMcache::Bf_Mz2_Mz2_mf2_mf2 ( const Particle  f) const

A cache method.

Parameters
[in]fa lepton or quark
Returns
\(B_f(M_Z^2;m_f^2,m_f^2)\) with \(\mu=M_Z\)

Definition at line 773 of file EWSMcache.cpp.

774 {
775  int NumPar = 2;
776  double params[] = {SM.getMz(), mf(f, SM.getMz())};
777  int ind = f.getIndex();
778  if (CacheCheck(Bf_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params))
779  return gslpp::complex(Bf_Mz2_Mz2_mf2_mf2_cache[ind][NumPar],
780  Bf_Mz2_Mz2_mf2_mf2_cache[ind][NumPar + 1], false);
781  else {
782  gslpp::complex newResult = PV.Bf(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), mf2(f, SM.getMz()), mf2(f, SM.getMz()));
783  newCacheForComplex(Bf_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params, newResult);
784  return newResult;
785  }
786 }

◆ Bfp_Mw2_Mw2_mfprime2_mf2()

gslpp::complex EWSMcache::Bfp_Mw2_Mw2_mfprime2_mf2 ( const int  gen,
const double  Mw_i 
) const

A cache method.

Parameters
[in]genthe generation index of a lepton doublet
[in]Mw_ithe \(W\)-boson mass
Returns
\(B_{fp}(M_W^2;m_{l'}^2,m_l^2)\) with \(\mu=M_W\)

Definition at line 900 of file EWSMcache.cpp.

901 {
902  int NumPar = 3;
903  double mymf, mymfprime;
904  double Mz = SM.getMz();
905  if (gen < 3) {
906  mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
907  mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
908  } else {
909  int genq = gen - 3;
910  mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
911  mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
912  }
913  double params[] = {Mw_i, mymf, mymfprime};
914 
915  if (CacheCheck(Bfp_Mw2_Mw2_mfprime2_mf2_cache[gen], NumPar, params))
916  return gslpp::complex(Bfp_Mw2_Mw2_mfprime2_mf2_cache[gen][NumPar],
917  Bfp_Mw2_Mw2_mfprime2_mf2_cache[gen][NumPar + 1], false);
918  else {
919  double mf2 = mymf*mymf;
920  double mfprime2 = mymfprime*mymfprime;
921  gslpp::complex newResult = PV.Bfp(Mw_i*Mw_i, Mw_i*Mw_i, mfprime2, mf2);
922  newCacheForComplex(Bfp_Mw2_Mw2_mfprime2_mf2_cache[gen], NumPar, params, newResult);
923  return newResult;
924  }
925 }

◆ Bfp_Mz2_Mz2_mf2_mf2()

gslpp::complex EWSMcache::Bfp_Mz2_Mz2_mf2_mf2 ( const Particle  f) const

A cache method.

Parameters
[in]fa lepton or quark
Returns
\(B_{fp}(M_Z^2;m_f^2,m_f^2)\) with \(\mu=M_Z\)

Definition at line 885 of file EWSMcache.cpp.

886 {
887  int NumPar = 2;
888  double params[] = {SM.getMz(), mf(f, SM.getMz())};
889  int ind = f.getIndex();
890  if (CacheCheck(Bfp_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params))
891  return gslpp::complex(Bfp_Mz2_Mz2_mf2_mf2_cache[ind][NumPar],
892  Bfp_Mz2_Mz2_mf2_mf2_cache[ind][NumPar + 1], false);
893  else {
894  gslpp::complex newResult = PV.Bfp(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), mf2(f, SM.getMz()), mf2(f, SM.getMz()));
895  newCacheForComplex(Bfp_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params, newResult);
896  return newResult;
897  }
898 }

◆ C0_Mw2_0_Mz2_0()

gslpp::complex EWSMcache::C0_Mw2_0_Mz2_0 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(C_0(M_W^2;0,M_Z^2,0)\)

Definition at line 1002 of file EWSMcache.cpp.

1003 {
1004  int NumPar = 2;
1005  double params[] = {Mw_i, SM.getMz()};
1006 
1007  if (CacheCheck(C0_Mw2_0_Mz2_0_cache, NumPar, params))
1008  return gslpp::complex(C0_Mw2_0_Mz2_0_cache[NumPar],
1009  C0_Mw2_0_Mz2_0_cache[NumPar + 1], false);
1010  else {
1011  gslpp::complex newResult = PV.C0(Mw_i*Mw_i, 0.0, SM.getMz() * SM.getMz(), 0.0);
1012  newCacheForComplex(C0_Mw2_0_Mz2_0_cache, NumPar, params, newResult);
1013  return newResult;
1014  }
1015 }

◆ C0_Mw2_Mw2_0_Mz2()

gslpp::complex EWSMcache::C0_Mw2_Mw2_0_Mz2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(C_0(M_W^2;M_W^2,0,M_Z^2)\)

Definition at line 987 of file EWSMcache.cpp.

988 {
989  int NumPar = 2;
990  double params[] = {Mw_i, SM.getMz()};
991 
992  if (CacheCheck(C0_Mw2_Mw2_0_Mz2_cache, NumPar, params))
993  return gslpp::complex(C0_Mw2_Mw2_0_Mz2_cache[NumPar],
994  C0_Mw2_Mw2_0_Mz2_cache[NumPar + 1], false);
995  else {
996  gslpp::complex newResult = PV.C0(Mw_i*Mw_i, Mw_i*Mw_i, 0.0, SM.getMz() * SM.getMz());
997  newCacheForComplex(C0_Mw2_Mw2_0_Mz2_cache, NumPar, params, newResult);
998  return newResult;
999  }
1000 }

◆ C0_Mz2_0_Mw2_0()

gslpp::complex EWSMcache::C0_Mz2_0_Mw2_0 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(C_0(M_Z^2;0,M_W^2,0)\)

Definition at line 957 of file EWSMcache.cpp.

958 {
959  int NumPar = 2;
960  double params[] = {SM.getMz(), Mw_i};
961 
962  if (CacheCheck(C0_Mz2_0_Mw2_0_cache, NumPar, params))
963  return gslpp::complex(C0_Mz2_0_Mw2_0_cache[NumPar],
964  C0_Mz2_0_Mw2_0_cache[NumPar + 1], false);
965  else {
966  gslpp::complex newResult = PV.C0(SM.getMz() * SM.getMz(), 0.0, Mw_i*Mw_i, 0.0);
967  newCacheForComplex(C0_Mz2_0_Mw2_0_cache, NumPar, params, newResult);
968  return newResult;
969  }
970 }

◆ C0_Mz2_0_Mz2_0()

gslpp::complex EWSMcache::C0_Mz2_0_Mz2_0 ( ) const

A cache method.

Returns
\(C_0(M_Z^2;0,M_Z^2,0)\)

Definition at line 1017 of file EWSMcache.cpp.

1018 {
1019  int NumPar = 1;
1020  double params[] = {SM.getMz()};
1021 
1022  if (CacheCheck(C0_Mz2_0_Mz2_0_cache, NumPar, params))
1023  return gslpp::complex(C0_Mz2_0_Mz2_0_cache[NumPar],
1024  C0_Mz2_0_Mz2_0_cache[NumPar + 1], false);
1025  else {
1026  gslpp::complex newResult = PV.C0(SM.getMz() * SM.getMz(), 0.0, SM.getMz() * SM.getMz(), 0.0);
1027  newCacheForComplex(C0_Mz2_0_Mz2_0_cache, NumPar, params, newResult);
1028  return newResult;
1029  }
1030 }

◆ C0_Mz2_Mt2_Mw2_Mt2()

gslpp::complex EWSMcache::C0_Mz2_Mt2_Mw2_Mt2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(C_0(M_Z^2;M_t^2,M_W^2,m_t^2)\)

Definition at line 942 of file EWSMcache.cpp.

943 {
944  int NumPar = 3;
945  double params[] = {SM.getMz(), SM.getMtpole(), Mw_i};
946 
947  if (CacheCheck(C0_Mz2_Mt2_Mw2_Mt2_cache, NumPar, params))
948  return gslpp::complex(C0_Mz2_Mt2_Mw2_Mt2_cache[NumPar],
949  C0_Mz2_Mt2_Mw2_Mt2_cache[NumPar + 1], false);
950  else {
951  gslpp::complex newResult = PV.C0(SM.getMz() * SM.getMz(), SM.getMtpole() * SM.getMtpole(), Mw_i*Mw_i, SM.getMtpole() * SM.getMtpole());
952  newCacheForComplex(C0_Mz2_Mt2_Mw2_Mt2_cache, NumPar, params, newResult);
953  return newResult;
954  }
955 }

◆ C0_Mz2_Mw2_0_Mw2()

gslpp::complex EWSMcache::C0_Mz2_Mw2_0_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(C_0(M_Z^2;M_W^2,0,M_W^2)\)

Definition at line 972 of file EWSMcache.cpp.

973 {
974  int NumPar = 2;
975  double params[] = {SM.getMz(), Mw_i};
976 
977  if (CacheCheck(C0_Mz2_Mw2_0_Mw2_cache, NumPar, params))
978  return gslpp::complex(C0_Mz2_Mw2_0_Mw2_cache[NumPar],
979  C0_Mz2_Mw2_0_Mw2_cache[NumPar + 1], false);
980  else {
981  gslpp::complex newResult = PV.C0(SM.getMz() * SM.getMz(), Mw_i*Mw_i, 0.0, Mw_i * Mw_i);
982  newCacheForComplex(C0_Mz2_Mw2_0_Mw2_cache, NumPar, params, newResult);
983  return newResult;
984  }
985 }

◆ C0_Mz2_Mw2_Mt2_Mw2()

gslpp::complex EWSMcache::C0_Mz2_Mw2_Mt2_Mw2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(C_0(M_Z^2;M_W^2,m_t^2,M_W^2)\)

Definition at line 927 of file EWSMcache.cpp.

928 {
929  int NumPar = 3;
930  double params[] = {SM.getMz(), Mw_i, SM.getMtpole()};
931 
932  if (CacheCheck(C0_Mz2_Mw2_Mt2_Mw2_cache, NumPar, params))
933  return gslpp::complex(C0_Mz2_Mw2_Mt2_Mw2_cache[NumPar],
934  C0_Mz2_Mw2_Mt2_Mw2_cache[NumPar + 1], false);
935  else {
936  gslpp::complex newResult = PV.C0(SM.getMz() * SM.getMz(), Mw_i*Mw_i, SM.getMtpole() * SM.getMtpole(), Mw_i * Mw_i);
937  newCacheForComplex(C0_Mz2_Mw2_Mt2_Mw2_cache, NumPar, params, newResult);
938  return newResult;
939  }
940 }

◆ CacheCheck()

bool EWSMcache::CacheCheck ( const double  cache[],
const int  NumPar,
const double  params[] 
) const
inlineprivate

A check method for caching.

This function checks if the values of the parameters in the array params[] are all identical to those stored in the array cache[]. When they are identical to each other, the current function returns true. Otherwise, the function returns false.

Parameters
[in]cachea cache of the parameters to be checked
[in]NumParthe number of the parameters to be checked
[in]paramsthe parameters to be checked
Returns
true (false) if the parameters in params[] are (not) identical to those in cache[].
Attention
If the flag FlagCacheInEWSMcache is set to false, the current function always returns false.

Definition at line 1008 of file EWSMcache.h.

1010  {
1011  if (!FlagCacheInEWSMcache) return false;
1012  bool bCache = true;
1013  for (int i = 0; i < NumPar; ++i)
1014  bCache &= (params[i] == cache[i]);
1015  return bCache;
1016  }

◆ Cl2_2Phi()

double EWSMcache::Cl2_2Phi ( ) const
inline

The constant \({\rm Cl}_2(2 \Phi)\).

This constant appears in two-loop QCD contribution, where \(\Phi=\mathrm{arcsin}\big(M_Z/(2m_t)\big)\).

Returns
\({\rm Cl}_2(2 \Phi)\)

Definition at line 438 of file EWSMcache.h.

439  {
440  double Phi = asin(SM.getMz() / 2.0 / SM.getMtpole());
441  return ( Clausen.Cl2(2.0 * Phi));
442  }

◆ Cl2_4Phi()

double EWSMcache::Cl2_4Phi ( ) const
inline

The constant \({\rm Cl}_2(4 \Phi)\).

This constant appears in two-loop QCD contribution, where \(\Phi=\mathrm{arcsin}\big(M_Z/(2m_t)\big)\).

Returns
\({\rm Cl}_2(4 \Phi)\)

Definition at line 450 of file EWSMcache.h.

451  {
452  double Phi = asin(SM.getMz() / 2.0 / SM.getMtpole());
453  return ( Clausen.Cl2(4.0 * Phi));
454  }

◆ Cl3_2Phi()

double EWSMcache::Cl3_2Phi ( ) const
inline

The constant \({\rm Cl}_3(2 \Phi)\).

This constant appears in two-loop QCD contribution, where \(\Phi=\mathrm{arcsin}\big(M_Z/(2m_t)\big)\).

Returns
\({\rm Cl}_3(2 \Phi)\)

Definition at line 462 of file EWSMcache.h.

463  {
464  double Phi = asin(SM.getMz() / 2.0 / SM.getMtpole());
465  return ( Clausen.Cl3(2.0 * Phi));
466  }

◆ Cl3_4Phi()

double EWSMcache::Cl3_4Phi ( ) const
inline

The constant \({\rm Cl}_3(4 \Phi)\).

This constant appears in two-loop QCD contribution, where \(\Phi=\mathrm{arcsin}\big(\sqrt{M_Z/(2m_t)}\big)\).

Returns
\({\rm Cl}_3(4 \Phi)\)

Definition at line 474 of file EWSMcache.h.

475  {
476  double Phi = asin(SM.getMz() / 2.0 / SM.getMtpole());
477  return ( Clausen.Cl3(4.0 * Phi));
478  }

◆ delta_f()

double EWSMcache::delta_f ( const Particle  f,
const double  Mw_i 
) const
inline

\(\delta_f = v_f-a_f\).

Parameters
[in]fa lepton or quark
[in]Mw_ithe \(W\)-boson mass
Returns
\(\delta_f\)

Definition at line 323 of file EWSMcache.h.

324  {
325  return ( -2.0 * Q_f(f) * SM.sW2(Mw_i));
326  }

◆ f_AlphaToGF()

double EWSMcache::f_AlphaToGF ( const double  Mw_i) const
inline

The conversion factor from \(\alpha\) to \(G_\mu\), denoted as \(f\).

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\( f = \sqrt{2}G_\mu M_Z^2 s_W^2 c_W^2/(\pi\alpha)\).

Definition at line 334 of file EWSMcache.h.

335  {
336  return ( sqrt(2.0) * SM.getGF() * pow(SM.getMz(), 2.0) * SM.sW2(Mw_i) * SM.cW2(Mw_i) / M_PI / SM.getAle());
337  }

◆ gamma_QCD2()

double EWSMcache::gamma_QCD2 ( ) const
inline

The constant \(\gamma\) for two-loop QCD contribution.

Returns
\(\gamma = \ln (2\sqrt{r})\) with \(r = M_Z^2/(4m_t^2)\)

Definition at line 400 of file EWSMcache.h.

401  {
402  double r_QCD2 = SM.getMz() * SM.getMz() / 4.0 / SM.getMtpole() / SM.getMtpole();
403  return ( log(2.0 * sqrt(r_QCD2)));
404  }

◆ getB4()

double EWSMcache::getB4 ( ) const
inline

A get method to access the constant \(B_4\).

The constant \(B_4\) is defined as

\[ B_4 = 16\, {\rm Li}_4 \left( \frac{1}{2} \right) - 4\, \zeta(2) \ln^2 2 + \frac{2}{3} \ln^4 2 - \frac{13}{2}\, \zeta(4) \approx -1.76280008707377\,, \]

which appears in three-loop amplitudes.

Returns
\(B_4\)

Definition at line 214 of file EWSMcache.h.

215  {
216  return B4;
217  }

◆ getClausen()

const ClausenFunctions EWSMcache::getClausen ( ) const
inline

A get method to access the member object of type ClausenFunctions.

Returns
the object of type ClausenFunctions in the current class

Definition at line 124 of file EWSMcache.h.

125  {
126  return Clausen;
127  }

◆ getD3()

double EWSMcache::getD3 ( ) const
inline

A get method to access the constant \(D_3\).

The constant \(D_3\) is defined as

\[ D_3 = 6\, \zeta(3) - \frac{15}{4}\, \zeta(4) - 6 \left[ {\rm Cl}_2 \left( \frac{\pi}{3} \right) \right]^2 \approx -3.02700949398765\,, \]

which appears in three-loop amplitudes.

Returns
\(D_3\)

Definition at line 197 of file EWSMcache.h.

198  {
199  return D3;
200  }

◆ getLog2()

double EWSMcache::getLog2 ( ) const
inline

A get method to access the constant \(\ln 2\).

Returns
\(\ln 2\)

Definition at line 223 of file EWSMcache.h.

224  {
225  return log2;
226  }

◆ getPolyLog()

const Polylogarithms EWSMcache::getPolyLog ( ) const
inline

A get method to access the member object of type Polylogarithms.

Returns
the object of type Polylogarithms in the current class

Definition at line 115 of file EWSMcache.h.

116  {
117  return PolyLog;
118  }

◆ getPV()

const PVfunctions EWSMcache::getPV ( ) const
inline

A get method to access the member reference to the object of type StandardModel passed to the constructor.

Returns
the reference to the object of type StandardModel passed to the constructor.

A get method to access the member object of type PVfunctions.

Returns
the object of type PVfunctions in the current class

Definition at line 106 of file EWSMcache.h.

107  {
108  return PV;
109  }

◆ getS2()

double EWSMcache::getS2 ( ) const
inline

A get method to access the constant \(S_2\).

The constant \(S_2\) is defined as

\[ S_2 = \frac{4}{9 \sqrt{3}}\, {\rm Cl}_2 \left( \frac{\pi}{3} \right) \approx 0.260434137632162\,, \]

which appears in three-loop amplitudes.

Returns
\(S_2\)

Definition at line 180 of file EWSMcache.h.

181  {
182  return S2;
183  }

◆ getSM()

const StandardModel& EWSMcache::getSM ( ) const
inline
Returns

Definition at line 56 of file EWSMcache.h.

57  {
58  return SM;
59  }

◆ getZeta2()

double EWSMcache::getZeta2 ( ) const
inline

A get method to access the value of the zeta function \(\zeta(2)\).

Returns
\(\zeta(2)\)

Definition at line 137 of file EWSMcache.h.

138  {
139  return zeta2;
140  }

◆ getZeta3()

double EWSMcache::getZeta3 ( ) const
inline

A get method to access the value of the zeta function \(\zeta(3)\).

Returns
\(\zeta(3)\)

Definition at line 146 of file EWSMcache.h.

147  {
148  return zeta3;
149  }

◆ getZeta4()

double EWSMcache::getZeta4 ( ) const
inline

A get method to access the value of the zeta function \(\zeta(4)\).

Returns
\(\zeta(4)\)

Definition at line 155 of file EWSMcache.h.

156  {
157  return zeta4;
158  }

◆ getZeta5()

double EWSMcache::getZeta5 ( ) const
inline

A get method to access the value of the zeta function \(\zeta(5)\).

Returns
\(\zeta(5)\)

Definition at line 164 of file EWSMcache.h.

165  {
166  return zeta5;
167  }

◆ h_QCD2()

double EWSMcache::h_QCD2 ( ) const
inline

The constant \(h\) for two-loop QCD contribution.

Returns
\(h= \ln(2\sqrt{1-r})\) with \(r = M_Z^2/(4m_t^2)\)

Definition at line 410 of file EWSMcache.h.

411  {
412  double r_QCD2 = SM.getMz() * SM.getMz() / 4.0 / SM.getMtpole() / SM.getMtpole();
413  return ( log(2.0 * sqrt(1.0 - r_QCD2)));
414  }

◆ I3_f()

double EWSMcache::I3_f ( const Particle  f) const
inline

The isospin of an SM fermion \(I_3^f\).

Parameters
[in]fa lepton or quark
Returns
\(I_3^f\)

Definition at line 278 of file EWSMcache.h.

279  {
280  return f.getIsospin();
281  }

◆ Li2_MW2toMTOP2()

double EWSMcache::Li2_MW2toMTOP2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(\mathrm{Li}_2(M_W^2/m_t^2)\)

Definition at line 154 of file EWSMcache.cpp.

155 {
156  int NumPar = 2;
157  double params[] = {Mw_i, SM.getMtpole()};
158 
159  if (CacheCheck(Li2_MW2toMTOP2_cache, NumPar, params))
160  return Li2_MW2toMTOP2_cache[NumPar];
161  else {
162  double newResult = PolyLog.Li2(Mw_i * Mw_i / SM.getMtpole() / SM.getMtpole()).real();
163  newCacheForDouble(Li2_MW2toMTOP2_cache, NumPar, params, newResult);
164  return newResult;
165  }
166 }

◆ Li3_for_F1()

double EWSMcache::Li3_for_F1 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(\mathrm{Li}_3(-M_W^2/m_t^2/(1-M_W^2/m_t^2))\)

Definition at line 182 of file EWSMcache.cpp.

183 {
184  int NumPar = 2;
185  double params[] = {Mw_i, SM.getMtpole()};
186 
187  if (CacheCheck(Li3_for_F1_cache, NumPar, params))
188  return Li3_for_F1_cache[NumPar];
189  else {
190  double tmp = Mw_i * Mw_i / SM.getMtpole() / SM.getMtpole();
191  double newResult = PolyLog.Li3(-tmp / (1.0 - tmp));
192  newCacheForDouble(Li3_for_F1_cache, NumPar, params, newResult);
193  return newResult;
194  }
195 }

◆ Li3_MW2toMTOP2()

double EWSMcache::Li3_MW2toMTOP2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(\mathrm{Li}_3(M_W^2/m_t^2)\)

Definition at line 168 of file EWSMcache.cpp.

169 {
170  int NumPar = 2;
171  double params[] = {Mw_i, SM.getMtpole()};
172 
173  if (CacheCheck(Li3_MW2toMTOP2_cache, NumPar, params))
174  return Li3_MW2toMTOP2_cache[NumPar];
175  else {
176  double newResult = PolyLog.Li3(Mw_i * Mw_i / SM.getMtpole() / SM.getMtpole());
177  newCacheForDouble(Li3_MW2toMTOP2_cache, NumPar, params, newResult);
178  return newResult;
179  }
180 }

◆ log_cW2()

double EWSMcache::log_cW2 ( const double  Mw_i) const

A cache method.

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(\ln c_W^2\)

Definition at line 140 of file EWSMcache.cpp.

141 {
142  int NumPar = 2;
143  double params[] = {SM.getMz(), Mw_i};
144 
145  if (CacheCheck(log_cW2_cache, NumPar, params))
146  return log_cW2_cache[NumPar];
147  else {
148  double newResult = log(SM.cW2(Mw_i));
149  newCacheForDouble(log_cW2_cache, NumPar, params, newResult);
150  return newResult;
151  }
152 }

◆ logMTOPtoMH()

double EWSMcache::logMTOPtoMH ( ) const

A cache method.

Returns
\(\ln (m_t/m_h)\)

Definition at line 126 of file EWSMcache.cpp.

127 {
128  int NumPar = 2;
129  double params[] = {SM.getMtpole(), SM.getMHl()};
130 
131  if (CacheCheck(logMTOPtoMH_cache, NumPar, params))
132  return logMTOPtoMH_cache[NumPar];
133  else {
134  double newResult = log(SM.getMtpole() / SM.getMHl());
135  newCacheForDouble(logMTOPtoMH_cache, NumPar, params, newResult);
136  return newResult;
137  }
138 }

◆ logMZtoME()

double EWSMcache::logMZtoME ( ) const

A cache method.

Returns
\(\ln (M_Z/m_e)\)

Definition at line 70 of file EWSMcache.cpp.

71 {
72  int NumPar = 2;
73  double params[] = {SM.getMz(), mf(SM.getLeptons(SM.ELECTRON))};
74 
75  if (CacheCheck(logMZtoME_cache, NumPar, params))
76  return logMZtoME_cache[NumPar];
77  else {
78  double newResult = log(SM.getMz() / mf(SM.getLeptons(SM.ELECTRON)));
79  newCacheForDouble(logMZtoME_cache, NumPar, params, newResult);
80  return newResult;
81  }
82 }

◆ logMZtoMMU()

double EWSMcache::logMZtoMMU ( ) const

A cache method.

Returns
\(\ln (M_Z/m_\mu)\)

Definition at line 84 of file EWSMcache.cpp.

85 {
86  int NumPar = 2;
87  double params[] = {SM.getMz(), mf(SM.getLeptons(SM.MU))};
88 
89  if (CacheCheck(logMZtoMMU_cache, NumPar, params))
90  return logMZtoMMU_cache[NumPar];
91  else {
92  double newResult = log(SM.getMz() / mf(SM.getLeptons(SM.MU)));
93  newCacheForDouble(logMZtoMMU_cache, NumPar, params, newResult);
94  return newResult;
95  }
96 }

◆ logMZtoMTAU()

double EWSMcache::logMZtoMTAU ( ) const

A cache method.

Returns
\(\ln (M_Z/m_\tau)\)

Definition at line 98 of file EWSMcache.cpp.

99 {
100  int NumPar = 2;
101  double params[] = {SM.getMz(), mf(SM.getLeptons(SM.TAU))};
102 
103  if (CacheCheck(logMZtoMTAU_cache, NumPar, params))
104  return logMZtoMTAU_cache[NumPar];
105  else {
106  double newResult = log(SM.getMz() / mf(SM.getLeptons(SM.TAU)));
107  newCacheForDouble(logMZtoMTAU_cache, NumPar, params, newResult);
108  return newResult;
109  }
110 }

◆ logMZtoMTOP()

double EWSMcache::logMZtoMTOP ( ) const

A cache method.

Returns
\(\ln (M_Z/m_t)\)

Definition at line 112 of file EWSMcache.cpp.

113 {
114  int NumPar = 2;
115  double params[] = {SM.getMz(), SM.getMtpole()};
116 
117  if (CacheCheck(logMZtoMTOP_cache, NumPar, params))
118  return logMZtoMTOP_cache[NumPar];
119  else {
120  double newResult = log(SM.getMz() / SM.getMtpole());
121  newCacheForDouble(logMZtoMTOP_cache, NumPar, params, newResult);
122  return newResult;
123  }
124 }

◆ logV1primeAndA1prime()

double EWSMcache::logV1primeAndA1prime ( ) const
inline

A logarithm appearing in the functions \(V_1'\) and \(A_1'\) for two-loop QCD contribution.

Returns
\(\mathrm{Re}[\ln (1-e^{2i\Phi}) - 2\ln (1-e^{4i\Phi})]\) with \(\Phi=\)Phi_QCD2()

Definition at line 422 of file EWSMcache.h.

423  {
424  gsl_complex OneMinusE2Iphi = gsl_complex_rect(1.0 - cos(2.0 * Phi_QCD2()),
425  -sin(2.0 * Phi_QCD2()));
426  gsl_complex OneMinusE4Iphi = gsl_complex_rect(1.0 - cos(4.0 * Phi_QCD2()),
427  -sin(4.0 * Phi_QCD2()));
428  return ( GSL_REAL(gsl_complex_log(OneMinusE2Iphi))
429  - 2.0 * GSL_REAL(gsl_complex_log(OneMinusE4Iphi)));
430  }

◆ mf()

double EWSMcache::mf ( const Particle  f,
const double  mu = 0.0,
const orders  order = FULLNNLO 
) const

The mass of an SM fermion.

Parameters
[in]fa lepton or quark
[in]murenormalization scale
[in]orderorder in QCD (= LO, FULLNLO, FULLNNLO[defalut])
Returns
the MSbar mass of u, d, s, c, b at the scale mu or the pole mass of t and leptons
Attention
If the flag FlagDebug is set to true, \(m_{u,d,s}(2\,\mathrm{GeV})\), \(m_c(m_c)\), \(m_b(m_b)\) or \(m_t^{\mathrm{pole}}\) is returned.

Definition at line 49 of file EWSMcache.cpp.

50 {
51  if (f.is("TOP"))
52  return SM.getMtpole(); // the pole mass
53  else if (f.is("QUARK") && !FlagDebug) {
54  /* These codes are slow and not effective. */
55  //if (mu == SM.getMz()) {
56  // if (FlagCacheInEWSMcache && order == FULLNNLO)
57  // if (SM.checkSMparams(mf_atMz_params_cache[f.getIndex()]))
58  // return mf_atMz_cache[f.getIndex()];
59  // mf_atMz_cache[f.getIndex()] = SM.Mrun(mu, f.getMass_scale(), f.getMass(), order);
60  // return mf_atMz_cache[f.getIndex()];
61  //}
62  return SM.Mrun(mu, f.getMass_scale(), f.getMass(), order);
63  } else
64  return f.getMass();
65 }

◆ mf2()

double EWSMcache::mf2 ( const Particle  f,
const double  mu = 0.0,
const orders  order = FULLNNLO 
) const
inline

The mass squared of an SM fermion.

Parameters
[in]fa lepton or quark
[in]murenormalization scale
[in]orderorder in QCD (= LO, FULLNLO, FULLNNLO[defalut])
Returns
the MSbar mass squared of u, d, s, c, b at the scale mu or the pole mass squared of t and leptons
Attention
If the flag FlagDebug is set to true, \((m_{u,d,s}(2\,\mathrm{GeV}))^2\), \((m_c(m_c))^2\), \((m_b(m_b))^2\) or \((m_t^{\mathrm{pole}})^2\) is returned.

Definition at line 257 of file EWSMcache.h.

258  {
259  double mf1 = mf(f, mu, order);
260  return ( mf1 * mf1);
261  }

◆ newCacheForComplex()

void EWSMcache::newCacheForComplex ( double  cache[],
const int  NumPar,
const double  params[],
const gslpp::complex  newResult 
) const
inlineprivate

A method to update a cache of the parameters and the quantity under consideration.

This function updates cache[] with params[] and newResult, where newResult for the quantity under consideration depends on the parameters in params[]. Both the values of params[] and that of newResult are stored into cache[], whose last two elements correspond to the real and imaginary parts of the latter.

Parameters
[out]cachea cache of the parameters and the quantity
[in]NumParthe number of the parameters
[in]paramsan array of the parameters
[in]newResultthe new result of the quantity
Attention
If the flag FlagCacheInEWSMcache is set to false, the current function does not modify cache[].

Definition at line 1058 of file EWSMcache.h.

1060  {
1061  if (!FlagCacheInEWSMcache) return;
1062  for (int i = 0; i < NumPar; ++i)
1063  cache[i] = params[i];
1064  cache[NumPar] = newResult.real();
1065  cache[NumPar + 1] = newResult.imag();
1066  }

◆ newCacheForDouble()

void EWSMcache::newCacheForDouble ( double  cache[],
const int  NumPar,
const double  params[],
const double  newResult 
) const
inlineprivate

A method to update a cache of the parameters and the quantity under consideration.

This function updates cache[] with params[] and newResult, where newResult for the quantity under consideration depends on the parameters in params[]. Both the values of params[] and that of newResult are stored into cache[], whose last element corresponds to the latter.

Parameters
[out]cachea cache of the parameters and the quantity
[in]NumParthe number of the parameters
[in]paramsan array of the parameters
[in]newResultthe new result of the quantity
Attention
If the flag FlagCacheInEWSMcache is set to false, the current function does not modify cache[].

Definition at line 1033 of file EWSMcache.h.

1035  {
1036  if (!FlagCacheInEWSMcache) return;
1037  for (int i = 0; i < NumPar; ++i)
1038  cache[i] = params[i];
1039  cache[NumPar] = newResult;
1040  }

◆ Phi_QCD2()

double EWSMcache::Phi_QCD2 ( ) const
inline

The constant \(\Phi\) for two-loop QCD contribution.

Returns
\(\Phi = {\rm arcsin}(\sqrt{r})\) with \(r = M_Z^2/(4m_t^2)\)

Definition at line 390 of file EWSMcache.h.

391  {
392  double r_QCD2 = SM.getMz() * SM.getMz() / 4.0 / SM.getMtpole() / SM.getMtpole();
393  return ( asin(sqrt(r_QCD2)));
394  }

◆ Q_f()

double EWSMcache::Q_f ( const Particle  f) const
inline

The charge of an SM fermion \(Q_f\).

Parameters
[in]fa lepton or quark
Returns
\(Q_f\)

Definition at line 268 of file EWSMcache.h.

269  {
270  return f.getCharge();
271  }

◆ setFlagCacheInEWSMcache()

void EWSMcache::setFlagCacheInEWSMcache ( bool  FlagCacheInEWSMcache)
inline

A set method to change the model flag CacheInEWSMcache in StandardModel.

Setting CacheInEWSMcache to false, the caching methods defined in the current class are not employed in numerical computations. The flag is set to true in the constructor EWSMcache() by default.

Parameters
[in]FlagCacheInEWSMcachea boolean flag for caching
See also
the description of the StandardModel flags

Definition at line 83 of file EWSMcache.h.

84  {
85  this->FlagCacheInEWSMcache = FlagCacheInEWSMcache;
86  }

◆ setFlagDebug()

void EWSMcache::setFlagDebug ( bool  FlagDebug)
inline

A set method to change the internal boolean flag FlagDebug.

The flag FlagDebug=true is used for testing and debugging the codes associated with the current class. The flag FlagDebug is set to false in the constructor EWSMcache() by default.

Parameters
[in]FlagDebuga boolean flag for debugging

Definition at line 68 of file EWSMcache.h.

69  {
70  this->FlagDebug = FlagDebug;
71  }

◆ sigma_f()

double EWSMcache::sigma_f ( const Particle  f,
const double  Mw_i 
) const
inline

\(\sigma_f = |v_f+a_f|\).

Parameters
[in]fa lepton or quark
[in]Mw_ithe \(W\)-boson mass
Returns
\(\sigma_f\)

Definition at line 312 of file EWSMcache.h.

313  {
314  return ( 1.0 - 2.0 * fabs(Q_f(f)) * SM.sW2(Mw_i));
315  }

◆ v_f()

double EWSMcache::v_f ( const Particle  f,
const double  Mw_i 
) const
inline

The tree-level vector coupling for \(Z\to f\bar{f}\), denoted as \(v_f\).

Parameters
[in]fa lepton or quark
[in]Mw_ithe \(W\)-boson mass
Returns
\(v_f\)

Definition at line 290 of file EWSMcache.h.

291  {
292  return ( a_f(f) - 2.0 * Q_f(f) * SM.sW2(Mw_i));
293  }

◆ Xt_alpha()

double EWSMcache::Xt_alpha ( const double  Mw_i) const
inline

The quantity \(X_t\) with the coupling \(\alpha\).

Parameters
[in]Mw_ithe \(W\)-boson mass
Returns
\(X_t^{\alpha}\)
See also
Xt_GF() and f_AlphaToGF()

Definition at line 355 of file EWSMcache.h.

356  {
357  return ( Xt_GF() / f_AlphaToGF(Mw_i));
358  }

◆ Xt_GF()

double EWSMcache::Xt_GF ( ) const
inline

The quantity \(X_t\) with the coupling \(G_\mu\).

Returns
\(X_t^{G_\mu}=G_\mu m_t^2/(8\sqrt{2}\pi^2)\)

Definition at line 343 of file EWSMcache.h.

344  {
345  return ( SM.getGF() * SM.getMtpole() * SM.getMtpole() / 8.0 / sqrt(2.0) / M_PI / M_PI);
346  }

Member Data Documentation

◆ A0_Mw2_mh2_cache

double EWSMcache::A0_Mw2_mh2_cache[3]
mutableprivate

A cache of a PV function.

Definition at line 925 of file EWSMcache.h.

◆ A0_Mw2_Mw2_cache

double EWSMcache::A0_Mw2_Mw2_cache[2]
mutableprivate

A cache of a PV function.

Definition at line 927 of file EWSMcache.h.

◆ A0_Mw2_Mz2_cache

double EWSMcache::A0_Mw2_Mz2_cache[3]
mutableprivate

A cache of a PV function.

Definition at line 924 of file EWSMcache.h.

◆ A0_Mz2_mh2_cache

double EWSMcache::A0_Mz2_mh2_cache[3]
mutableprivate

A cache of a PV function.

Definition at line 923 of file EWSMcache.h.

◆ A0_Mz2_Mw2_cache

double EWSMcache::A0_Mz2_Mw2_cache[3]
mutableprivate

A cache of a PV function.

Definition at line 922 of file EWSMcache.h.

◆ A0_Mz2_Mz2_cache

double EWSMcache::A0_Mz2_Mz2_cache[2]
mutableprivate

A cache of a PV function.

Definition at line 926 of file EWSMcache.h.

◆ B0_Mw2_Mw2_0_Mw2_cache

double EWSMcache::B0_Mw2_Mw2_0_Mw2_cache[3]
mutableprivate

A cache of a PV function.

Definition at line 952 of file EWSMcache.h.

◆ B0_Mw2_Mw2_mh2_Mw2_cache

double EWSMcache::B0_Mw2_Mw2_mh2_Mw2_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 951 of file EWSMcache.h.

◆ B0_Mw2_Mw2_Mz2_Mw2_cache

double EWSMcache::B0_Mw2_Mw2_Mz2_Mw2_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 950 of file EWSMcache.h.

◆ B0_Mw2_Mz2_Mt2_Mt2_cache

double EWSMcache::B0_Mw2_Mz2_Mt2_Mt2_cache[5]
mutableprivate

A cache of a PV function.

Definition at line 942 of file EWSMcache.h.

◆ B0_Mw2_Mz2_Mw2_Mw2_cache

double EWSMcache::B0_Mw2_Mz2_Mw2_Mw2_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 949 of file EWSMcache.h.

◆ B0_Mz2_0_0_Mw2_cache

double EWSMcache::B0_Mz2_0_0_Mw2_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 948 of file EWSMcache.h.

◆ B0_Mz2_0_mh2_Mw2_cache

double EWSMcache::B0_Mz2_0_mh2_Mw2_cache[5]
mutableprivate

A cache of a PV function.

Definition at line 941 of file EWSMcache.h.

◆ B0_Mz2_0_Mz2_Mw2_cache

double EWSMcache::B0_Mz2_0_Mz2_Mw2_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 947 of file EWSMcache.h.

◆ B0_Mz2_Mw2_0_Mw2_cache

double EWSMcache::B0_Mz2_Mw2_0_Mw2_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 946 of file EWSMcache.h.

◆ B0_Mz2_Mw2_mh2_Mw2_cache

double EWSMcache::B0_Mz2_Mw2_mh2_Mw2_cache[5]
mutableprivate

A cache of a PV function.

Definition at line 940 of file EWSMcache.h.

◆ B0_Mz2_Mw2_Mz2_Mw2_cache

double EWSMcache::B0_Mz2_Mw2_Mz2_Mw2_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 945 of file EWSMcache.h.

◆ B0_Mz2_Mz2_mf2_mf2_cache

double EWSMcache::B0_Mz2_Mz2_mf2_mf2_cache[12][4]
mutableprivate

A cache of a PV function.

Definition at line 953 of file EWSMcache.h.

◆ B0_Mz2_Mz2_mh2_Mz2_cache

double EWSMcache::B0_Mz2_Mz2_mh2_Mz2_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 944 of file EWSMcache.h.

◆ B0_Mz2_Mz2_Mw2_Mw2_cache

double EWSMcache::B0_Mz2_Mz2_Mw2_Mw2_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 943 of file EWSMcache.h.

◆ B0p_Mw2_Mw2_0_Mw2_cache

double EWSMcache::B0p_Mw2_Mw2_0_Mw2_cache[3]
mutableprivate

A cache of a PV function.

Definition at line 961 of file EWSMcache.h.

◆ B0p_Mw2_Mw2_mh2_Mw2_cache

double EWSMcache::B0p_Mw2_Mw2_mh2_Mw2_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 960 of file EWSMcache.h.

◆ B0p_Mw2_Mw2_Mz2_Mw2_cache

double EWSMcache::B0p_Mw2_Mw2_Mz2_Mw2_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 959 of file EWSMcache.h.

◆ B0p_Mz2_0_mh2_Mw2_cache

double EWSMcache::B0p_Mz2_0_mh2_Mw2_cache[5]
mutableprivate

A cache of a PV function.

Definition at line 955 of file EWSMcache.h.

◆ B0p_Mz2_0_Mz2_Mw2_cache

double EWSMcache::B0p_Mz2_0_Mz2_Mw2_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 957 of file EWSMcache.h.

◆ B0p_Mz2_Mz2_mf2_mf2_cache

double EWSMcache::B0p_Mz2_Mz2_mf2_mf2_cache[12][4]
mutableprivate

A cache of a PV function.

Definition at line 962 of file EWSMcache.h.

◆ B0p_Mz2_Mz2_mh2_Mz2_cache

double EWSMcache::B0p_Mz2_Mz2_mh2_Mz2_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 956 of file EWSMcache.h.

◆ B0p_Mz2_Mz2_Mw2_Mw2_cache

double EWSMcache::B0p_Mz2_Mz2_Mw2_Mw2_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 958 of file EWSMcache.h.

◆ B1_Mz2_0_mf2_mfprime2_cache

double EWSMcache::B1_Mz2_0_mf2_mfprime2_cache[6][5]
mutableprivate

A cache of a PV function.

Definition at line 964 of file EWSMcache.h.

◆ B1_Mz2_0_mfprime2_mf2_cache

double EWSMcache::B1_Mz2_0_mfprime2_mf2_cache[6][5]
mutableprivate

A cache of a PV function.

Definition at line 965 of file EWSMcache.h.

◆ B1_Mz2_Mw2_mf2_mfprime2_cache

double EWSMcache::B1_Mz2_Mw2_mf2_mfprime2_cache[6][6]
mutableprivate

A cache of a PV function.

Definition at line 966 of file EWSMcache.h.

◆ B1_Mz2_Mw2_mfprime2_mf2_cache

double EWSMcache::B1_Mz2_Mw2_mfprime2_mf2_cache[6][6]
mutableprivate

A cache of a PV function.

Definition at line 967 of file EWSMcache.h.

◆ B1p_Mw2_Mw2_mf2_mfprime2_cache

double EWSMcache::B1p_Mw2_Mw2_mf2_mfprime2_cache[6][5]
mutableprivate

A cache of a PV function.

Definition at line 969 of file EWSMcache.h.

◆ B1p_Mw2_Mw2_mfprime2_mf2_cache

double EWSMcache::B1p_Mw2_Mw2_mfprime2_mf2_cache[6][5]
mutableprivate

A cache of a PV function.

Definition at line 970 of file EWSMcache.h.

◆ B4

double EWSMcache::B4
private

The constant \(B_4=16\mathrm{Li}_4(1/2) - 4\zeta(2)\ln^2(2) + (2/3)\ln^4(2) - (13/2)\zeta(4)\).

Definition at line 899 of file EWSMcache.h.

◆ Bf_Mw2_Mw2_mfprime2_mf2_cache

double EWSMcache::Bf_Mw2_Mw2_mfprime2_mf2_cache[6][5]
mutableprivate

A cache of a PV function.

Definition at line 976 of file EWSMcache.h.

◆ Bf_Mz2_0_mf2_mf2_cache

double EWSMcache::Bf_Mz2_0_mf2_mf2_cache[12][4]
mutableprivate

A cache of a PV function.

Definition at line 973 of file EWSMcache.h.

◆ Bf_Mz2_0_mfprime2_mf2_cache

double EWSMcache::Bf_Mz2_0_mfprime2_mf2_cache[6][5]
mutableprivate

A cache of a PV function.

Definition at line 975 of file EWSMcache.h.

◆ Bf_Mz2_Mw2_mfprime2_mf2_cache

double EWSMcache::Bf_Mz2_Mw2_mfprime2_mf2_cache[6][6]
mutableprivate

A cache of a PV function.

Definition at line 974 of file EWSMcache.h.

◆ Bf_Mz2_Mz2_mf2_mf2_cache

double EWSMcache::Bf_Mz2_Mz2_mf2_mf2_cache[12][4]
mutableprivate

A cache of a PV function.

Definition at line 972 of file EWSMcache.h.

◆ Bfp_Mw2_Mw2_mfprime2_mf2_cache

double EWSMcache::Bfp_Mw2_Mw2_mfprime2_mf2_cache[6][5]
mutableprivate

A cache of a PV function.

Definition at line 979 of file EWSMcache.h.

◆ Bfp_Mz2_Mz2_mf2_mf2_cache

double EWSMcache::Bfp_Mz2_Mz2_mf2_mf2_cache[12][4]
mutableprivate

A cache of a PV function.

Definition at line 978 of file EWSMcache.h.

◆ C0_Mw2_0_Mz2_0_cache

double EWSMcache::C0_Mw2_0_Mz2_0_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 986 of file EWSMcache.h.

◆ C0_Mw2_Mw2_0_Mz2_cache

double EWSMcache::C0_Mw2_Mw2_0_Mz2_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 985 of file EWSMcache.h.

◆ C0_Mz2_0_Mw2_0_cache

double EWSMcache::C0_Mz2_0_Mw2_0_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 983 of file EWSMcache.h.

◆ C0_Mz2_0_Mz2_0_cache

double EWSMcache::C0_Mz2_0_Mz2_0_cache[3]
mutableprivate

A cache of a PV function.

Definition at line 987 of file EWSMcache.h.

◆ C0_Mz2_Mt2_Mw2_Mt2_cache

double EWSMcache::C0_Mz2_Mt2_Mw2_Mt2_cache[5]
mutableprivate

A cache of a PV function.

Definition at line 982 of file EWSMcache.h.

◆ C0_Mz2_Mw2_0_Mw2_cache

double EWSMcache::C0_Mz2_Mw2_0_Mw2_cache[4]
mutableprivate

A cache of a PV function.

Definition at line 984 of file EWSMcache.h.

◆ C0_Mz2_Mw2_Mt2_Mw2_cache

double EWSMcache::C0_Mz2_Mw2_Mt2_Mw2_cache[5]
mutableprivate

A cache of a PV function.

Definition at line 981 of file EWSMcache.h.

◆ Clausen

const ClausenFunctions EWSMcache::Clausen
private

An object of type ClausenFunctions.

Definition at line 889 of file EWSMcache.h.

◆ D3

double EWSMcache::D3
private

The constant \(D_3=6\zeta(3) - (15/4)\zeta(4) - 6[\mathrm{Cl}_2(\pi/3)]^2\).

Definition at line 898 of file EWSMcache.h.

◆ FlagCacheInEWSMcache

bool EWSMcache::FlagCacheInEWSMcache
private

A flag for caching (true by default).

Definition at line 885 of file EWSMcache.h.

◆ FlagDebug

bool EWSMcache::FlagDebug
private

A flag for debugging (false by default).

Definition at line 884 of file EWSMcache.h.

◆ Li2_MW2toMTOP2_cache

double EWSMcache::Li2_MW2toMTOP2_cache[3]
mutableprivate

A cache of \(\mathrm{Li}_2(M_W^2/m_t^2)\).

Definition at line 917 of file EWSMcache.h.

◆ Li3_for_F1_cache

double EWSMcache::Li3_for_F1_cache[3]
mutableprivate

A cache of \(\mathrm{Li}_3(-M_W^2/m_t^2/(1-M_W^2/m_t^2))\).

Definition at line 919 of file EWSMcache.h.

◆ Li3_MW2toMTOP2_cache

double EWSMcache::Li3_MW2toMTOP2_cache[3]
mutableprivate

A cache of \(\mathrm{Li}_3(M_W^2/m_t^2)\).

Definition at line 918 of file EWSMcache.h.

◆ log2

double EWSMcache::log2
private

The constant \(\ln 2\).

Definition at line 900 of file EWSMcache.h.

◆ log_cW2_cache

double EWSMcache::log_cW2_cache[3]
mutableprivate

A cache of \(\ln c_W^2\).

Definition at line 914 of file EWSMcache.h.

◆ logMTOPtoMH_cache

double EWSMcache::logMTOPtoMH_cache[3]
mutableprivate

A cache of \(\ln (m_t/m_h)\).

Definition at line 913 of file EWSMcache.h.

◆ logMZtoME_cache

double EWSMcache::logMZtoME_cache[3]
mutableprivate

A cache of \(\ln (M_Z/m_e)\).

Definition at line 909 of file EWSMcache.h.

◆ logMZtoMMU_cache

double EWSMcache::logMZtoMMU_cache[3]
mutableprivate

A cache of \(\ln (M_Z/m_\mu)\).

Definition at line 910 of file EWSMcache.h.

◆ logMZtoMTAU_cache

double EWSMcache::logMZtoMTAU_cache[3]
mutableprivate

A cache of \(\ln (M_Z/m_\tau)\).

Definition at line 911 of file EWSMcache.h.

◆ logMZtoMTOP_cache

double EWSMcache::logMZtoMTOP_cache[3]
mutableprivate

A cache of \(\ln (M_Z/m_t)\).

Definition at line 912 of file EWSMcache.h.

◆ mf_atMz_cache

double EWSMcache::mf_atMz_cache[12]
mutableprivate

A cache of the fermion masses at \(\mu=M_Z\).

Definition at line 933 of file EWSMcache.h.

◆ mf_atMz_params_cache

double EWSMcache::mf_atMz_params_cache[12][StandardModel::NumSMParamsForEWPO]
mutableprivate

A cache array of a set of SM parameters, used together with mf_atMz_cache.

Definition at line 932 of file EWSMcache.h.

◆ PolyLog

const Polylogarithms EWSMcache::PolyLog
private

An object of type Polylogarithms.

Definition at line 890 of file EWSMcache.h.

◆ PV

const PVfunctions EWSMcache::PV
private

An object of type PVfunctions.

Definition at line 888 of file EWSMcache.h.

◆ S2

double EWSMcache::S2
private

The constant \(S_2=(4/9/\sqrt{3})\mathrm{Cl}_2(\pi/3)\).

Definition at line 897 of file EWSMcache.h.

◆ SM

const StandardModel& EWSMcache::SM
private

A reference to an object of type StandardModel.

Definition at line 887 of file EWSMcache.h.

◆ zeta2

double EWSMcache::zeta2
private

The constant \(\zeta(2)\).

Definition at line 893 of file EWSMcache.h.

◆ zeta3

double EWSMcache::zeta3
private

The constant \(\zeta(3)\).

Definition at line 894 of file EWSMcache.h.

◆ zeta4

double EWSMcache::zeta4
private

The constant \(\zeta(4)\).

Definition at line 895 of file EWSMcache.h.

◆ zeta5

double EWSMcache::zeta5
private

The constant \(\zeta(5)\).

Definition at line 896 of file EWSMcache.h.

The documentation for this class was generated from the following files:
QCD::TAU
Definition: QCD.h:316
EWSMcache::mf2
double mf2(const Particle f, const double mu=0.0, const orders order=FULLNNLO) const
The mass squared of an SM fermion.
Definition: EWSMcache.h:257
EWSMcache::B0p_Mz2_Mz2_Mw2_Mw2_cache
double B0p_Mz2_Mz2_Mw2_Mw2_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:958
EWSMcache::B0p_Mw2_Mw2_0_Mw2_cache
double B0p_Mw2_Mw2_0_Mw2_cache[3]
A cache of a PV function.
Definition: EWSMcache.h:961
gslpp::cos
complex cos(const complex &z)
Definition: gslpp_complex.cpp:429
EWSMcache::log2
double log2
The constant .
Definition: EWSMcache.h:900
EWSMcache::Bf_Mz2_0_mf2_mf2_cache
double Bf_Mz2_0_mf2_mf2_cache[12][4]
A cache of a PV function.
Definition: EWSMcache.h:973
StandardModel::cW2
virtual double cW2(const double Mw_i) const
The square of the cosine of the weak mixing angle in the on-shell scheme, denoted as .
Definition: StandardModel.cpp:989
EWSMcache::A0_Mw2_mh2_cache
double A0_Mw2_mh2_cache[3]
A cache of a PV function.
Definition: EWSMcache.h:925
PVfunctions::A0
double A0(const double mu2, const double m2) const
.
Definition: PVfunctions.cpp:23
EWSMcache::logMZtoMTOP_cache
double logMZtoMTOP_cache[3]
A cache of .
Definition: EWSMcache.h:912
EWSMcache::Li2_MW2toMTOP2_cache
double Li2_MW2toMTOP2_cache[3]
A cache of .
Definition: EWSMcache.h:917
PVfunctions::B1p
gslpp::complex B1p(const double mu2, const double p2, const double m02, const double m12) const
.
Definition: PVfunctions.cpp:356
EWSMcache::zeta2
double zeta2
The constant .
Definition: EWSMcache.h:893
EWSMcache::B0p_Mz2_0_mh2_Mw2_cache
double B0p_Mz2_0_mh2_Mw2_cache[5]
A cache of a PV function.
Definition: EWSMcache.h:955
Particle::is
bool is(std::string name_i) const
Definition: Particle.cpp:23
EWSMcache::B0_Mz2_Mw2_mh2_Mw2_cache
double B0_Mz2_Mw2_mh2_Mw2_cache[5]
A cache of a PV function.
Definition: EWSMcache.h:940
EWSMcache::B0_Mz2_0_Mz2_Mw2_cache
double B0_Mz2_0_Mz2_Mw2_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:947
PVfunctions::Bf
gslpp::complex Bf(const double mu2, const double p2, const double m02, const double m12) const
.
Definition: PVfunctions.cpp:265
EWSMcache::zeta5
double zeta5
The constant .
Definition: EWSMcache.h:896
EWSMcache::f_AlphaToGF
double f_AlphaToGF(const double Mw_i) const
The conversion factor from to , denoted as .
Definition: EWSMcache.h:334
EWSMcache::Li3_for_F1_cache
double Li3_for_F1_cache[3]
A cache of .
Definition: EWSMcache.h:919
StandardModel::NumSMParamsForEWPO
static const int NumSMParamsForEWPO
The number of the SM parameters that are relevant to the EW precision observables.
Definition: StandardModel.h:1849
EWSMcache::logMZtoMMU_cache
double logMZtoMMU_cache[3]
A cache of .
Definition: EWSMcache.h:910
EWSMcache::B0_Mz2_Mz2_mf2_mf2_cache
double B0_Mz2_Mz2_mf2_mf2_cache[12][4]
A cache of a PV function.
Definition: EWSMcache.h:953
EWSMcache::B0_Mz2_0_0_Mw2_cache
double B0_Mz2_0_0_Mw2_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:948
gslpp::sin
complex sin(const complex &z)
Definition: gslpp_complex.cpp:420
EWSMcache::C0_Mz2_Mw2_Mt2_Mw2_cache
double C0_Mz2_Mw2_Mt2_Mw2_cache[5]
A cache of a PV function.
Definition: EWSMcache.h:981
EWSMcache::C0_Mz2_Mw2_0_Mw2_cache
double C0_Mz2_Mw2_0_Mw2_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:984
EWSMcache::A0_Mz2_Mw2_cache
double A0_Mz2_Mw2_cache[3]
A cache of a PV function.
Definition: EWSMcache.h:922
Polylogarithms::Li3
double Li3(const double x) const
The trilogarithm .
Definition: Polylogarithms.cpp:36
EWSMcache::newCacheForDouble
void newCacheForDouble(double cache[], const int NumPar, const double params[], const double newResult) const
A method to update a cache of the parameters and the quantity under consideration.
Definition: EWSMcache.h:1033
EWSMcache::logMZtoME_cache
double logMZtoME_cache[3]
A cache of .
Definition: EWSMcache.h:909
Polylogarithms::Li2
gslpp::complex Li2(const double x) const
The dilogarithm with a real argument, .
Definition: Polylogarithms.cpp:22
EWSMcache::D3
double D3
The constant .
Definition: EWSMcache.h:898
EWSMcache::Bfp_Mz2_Mz2_mf2_mf2_cache
double Bfp_Mz2_Mz2_mf2_mf2_cache[12][4]
A cache of a PV function.
Definition: EWSMcache.h:978
gslpp::complex
A class for defining operations on and functions of complex numbers.
Definition: gslpp_complex.h:35
gslpp::log
complex log(const complex &z)
Definition: gslpp_complex.cpp:342
EWSMcache::B0_Mz2_Mw2_0_Mw2_cache
double B0_Mz2_Mw2_0_Mw2_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:946
Particle::getMass_scale
double getMass_scale() const
A get method to access the scale at which the particle mass is defined.
Definition: Particle.h:133
EWSMcache::B0_Mz2_Mz2_mh2_Mz2_cache
double B0_Mz2_Mz2_mh2_Mz2_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:944
EWSMcache::B1p_Mw2_Mw2_mf2_mfprime2_cache
double B1p_Mw2_Mw2_mf2_mfprime2_cache[6][5]
A cache of a PV function.
Definition: EWSMcache.h:969
QCD::ELECTRON
Definition: QCD.h:312
Particle::getIsospin
double getIsospin() const
A get method to access the particle isospin.
Definition: Particle.h:115
EWSMcache::PV
const PVfunctions PV
An object of type PVfunctions.
Definition: EWSMcache.h:888
EWSMcache::Bf_Mw2_Mw2_mfprime2_mf2_cache
double Bf_Mw2_Mw2_mfprime2_mf2_cache[6][5]
A cache of a PV function.
Definition: EWSMcache.h:976
EWSMcache::B0p_Mz2_0_Mz2_Mw2_cache
double B0p_Mz2_0_Mz2_Mw2_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:957
StandardModel::sW2
virtual double sW2(const double Mw_i) const
The square of the sine of the weak mixing angle in the on-shell scheme, denoted as .
Definition: StandardModel.cpp:1000
EWSMcache::zeta4
double zeta4
The constant .
Definition: EWSMcache.h:895
EWSMcache::log_cW2_cache
double log_cW2_cache[3]
A cache of .
Definition: EWSMcache.h:914
EWSMcache::B0_Mw2_Mw2_0_Mw2_cache
double B0_Mw2_Mw2_0_Mw2_cache[3]
A cache of a PV function.
Definition: EWSMcache.h:952
EWSMcache::mf_atMz_cache
double mf_atMz_cache[12]
A cache of the fermion masses at .
Definition: EWSMcache.h:933
ClausenFunctions::Cl3
double Cl3(const double phi) const
The Clausen function of index 3, .
Definition: ClausenFunctions.cpp:26
EWSMcache::B0p_Mw2_Mw2_Mz2_Mw2_cache
double B0p_Mw2_Mw2_Mz2_Mw2_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:959
gslpp::complex::imag
const double & imag() const
Definition: gslpp_complex.cpp:59
EWSMcache::PolyLog
const Polylogarithms PolyLog
An object of type Polylogarithms.
Definition: EWSMcache.h:890
EWSMcache::B0_Mw2_Mz2_Mw2_Mw2_cache
double B0_Mw2_Mz2_Mw2_Mw2_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:949
Particle::getMass
const double & getMass() const
A get method to access the particle mass.
Definition: Particle.h:61
EWSMcache::Phi_QCD2
double Phi_QCD2() const
The constant for two-loop QCD contribution.
Definition: EWSMcache.h:390
EWSMcache::C0_Mz2_0_Mz2_0_cache
double C0_Mz2_0_Mz2_0_cache[3]
A cache of a PV function.
Definition: EWSMcache.h:987
EWSMcache::B1p_Mw2_Mw2_mfprime2_mf2_cache
double B1p_Mw2_Mw2_mfprime2_mf2_cache[6][5]
A cache of a PV function.
Definition: EWSMcache.h:970
EWSMcache::C0_Mz2_Mt2_Mw2_Mt2_cache
double C0_Mz2_Mt2_Mw2_Mt2_cache[5]
A cache of a PV function.
Definition: EWSMcache.h:982
ClausenFunctions::Cl2
double Cl2(const double phi) const
The Clausen function of index 2, .
Definition: ClausenFunctions.cpp:21
PVfunctions::C0
gslpp::complex C0(const double p2, const double m02, const double m12, const double m22) const
.
Definition: PVfunctions.cpp:469
EWSMcache::mf
double mf(const Particle f, const double mu=0.0, const orders order=FULLNNLO) const
The mass of an SM fermion.
Definition: EWSMcache.cpp:49
EWSMcache::B4
double B4
The constant .
Definition: EWSMcache.h:899
EWSMcache::B1_Mz2_Mw2_mf2_mfprime2_cache
double B1_Mz2_Mw2_mf2_mfprime2_cache[6][6]
A cache of a PV function.
Definition: EWSMcache.h:966
gslpp::pow
complex pow(const complex &z1, const complex &z2)
Definition: gslpp_complex.cpp:395
StandardModel::getGF
double getGF() const
A get method to retrieve the Fermi constant .
Definition: StandardModel.h:736
gslpp::sqrt
complex sqrt(const complex &z)
Definition: gslpp_complex.cpp:385
EWSMcache::C0_Mw2_Mw2_0_Mz2_cache
double C0_Mw2_Mw2_0_Mz2_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:985
EWSMcache::A0_Mw2_Mz2_cache
double A0_Mw2_Mz2_cache[3]
A cache of a PV function.
Definition: EWSMcache.h:924
QCD::getMtpole
double getMtpole() const
A get method to access the pole mass of the top quark.
Definition: QCD.h:588
Particle::getCharge
double getCharge() const
A get method to access the particle charge.
Definition: Particle.h:97
EWSMcache::B1_Mz2_0_mf2_mfprime2_cache
double B1_Mz2_0_mf2_mfprime2_cache[6][5]
A cache of a PV function.
Definition: EWSMcache.h:964
EWSMcache::newCacheForComplex
void newCacheForComplex(double cache[], const int NumPar, const double params[], const gslpp::complex newResult) const
A method to update a cache of the parameters and the quantity under consideration.
Definition: EWSMcache.h:1058
EWSMcache::B0_Mz2_Mz2_Mw2_Mw2_cache
double B0_Mz2_Mz2_Mw2_Mw2_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:943
EWSMcache::B1_Mz2_0_mfprime2_mf2_cache
double B1_Mz2_0_mfprime2_mf2_cache[6][5]
A cache of a PV function.
Definition: EWSMcache.h:965
EWSMcache::S2
double S2
The constant .
Definition: EWSMcache.h:897
EWSMcache::B0p_Mz2_Mz2_mf2_mf2_cache
double B0p_Mz2_Mz2_mf2_mf2_cache[12][4]
A cache of a PV function.
Definition: EWSMcache.h:962
EWSMcache::a_f
double a_f(const Particle f) const
The tree-level axial-vector coupling for , denoted as .
Definition: EWSMcache.h:301
EWSMcache::logMZtoMTAU_cache
double logMZtoMTAU_cache[3]
A cache of .
Definition: EWSMcache.h:911
PVfunctions::B0p
gslpp::complex B0p(const double muIR2, const double p2, const double m02, const double m12) const
.
Definition: PVfunctions.cpp:276
EWSMcache::B0_Mz2_0_mh2_Mw2_cache
double B0_Mz2_0_mh2_Mw2_cache[5]
A cache of a PV function.
Definition: EWSMcache.h:941
StandardModel::Als
double Als(double mu, orders order=FULLNLO, bool qed_flag=false, bool Nf_thr=true) const
The running QCD coupling in the scheme including QED corrections.
Definition: StandardModel.cpp:576
EWSMcache::B1_Mz2_Mw2_mfprime2_mf2_cache
double B1_Mz2_Mw2_mfprime2_mf2_cache[6][6]
A cache of a PV function.
Definition: EWSMcache.h:967
QCD::getQuarks
Particle getQuarks(const QCD::quark q) const
A get method to access a quark as an object of the type Particle.
Definition: QCD.h:534
EWSMcache::Bf_Mz2_Mz2_mf2_mf2_cache
double Bf_Mz2_Mz2_mf2_mf2_cache[12][4]
A cache of a PV function.
Definition: EWSMcache.h:972
EWSMcache::Bfp_Mw2_Mw2_mfprime2_mf2_cache
double Bfp_Mw2_Mw2_mfprime2_mf2_cache[6][5]
A cache of a PV function.
Definition: EWSMcache.h:979
EWSMcache::SM
const StandardModel & SM
A reference to an object of type StandardModel.
Definition: EWSMcache.h:887
QCD::quark
quark
An enum type for quarks.
Definition: QCD.h:323
EWSMcache::C0_Mw2_0_Mz2_0_cache
double C0_Mw2_0_Mz2_0_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:986
EWSMcache::B0_Mz2_Mw2_Mz2_Mw2_cache
double B0_Mz2_Mw2_Mz2_Mw2_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:945
EWSMcache::A0_Mz2_Mz2_cache
double A0_Mz2_Mz2_cache[2]
A cache of a PV function.
Definition: EWSMcache.h:926
EWSMcache::B0_Mw2_Mz2_Mt2_Mt2_cache
double B0_Mw2_Mz2_Mt2_Mt2_cache[5]
A cache of a PV function.
Definition: EWSMcache.h:942
EWSMcache::Bf_Mz2_Mw2_mfprime2_mf2_cache
double Bf_Mz2_Mw2_mfprime2_mf2_cache[6][6]
A cache of a PV function.
Definition: EWSMcache.h:974
EWSMcache::zeta3
double zeta3
The constant .
Definition: EWSMcache.h:894
EWSMcache::B0_Mw2_Mw2_mh2_Mw2_cache
double B0_Mw2_Mw2_mh2_Mw2_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:951
EWSMcache::A0_Mw2_Mw2_cache
double A0_Mw2_Mw2_cache[2]
A cache of a PV function.
Definition: EWSMcache.h:927
EWSMcache::Xt_GF
double Xt_GF() const
The quantity with the coupling .
Definition: EWSMcache.h:343
StandardModel::getMz
double getMz() const
A get method to access the mass of the boson .
Definition: StandardModel.h:718
EWSMcache::FlagDebug
bool FlagDebug
A flag for debugging (false by default).
Definition: EWSMcache.h:884
EWSMcache::CacheCheck
bool CacheCheck(const double cache[], const int NumPar, const double params[]) const
A check method for caching.
Definition: EWSMcache.h:1008
EWSMcache::Li3_MW2toMTOP2_cache
double Li3_MW2toMTOP2_cache[3]
A cache of .
Definition: EWSMcache.h:918
EWSMcache::B0p_Mz2_Mz2_mh2_Mz2_cache
double B0p_Mz2_Mz2_mh2_Mz2_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:956
EWSMcache::A0_Mz2_mh2_cache
double A0_Mz2_mh2_cache[3]
A cache of a PV function.
Definition: EWSMcache.h:923
QCD::Mrun
double Mrun(const double mu, const double m, const orders order=FULLNNLO) const
Computes a running quark mass from .
Definition: QCD.cpp:1064
EWSMcache::mf_atMz_params_cache
double mf_atMz_params_cache[12][StandardModel::NumSMParamsForEWPO]
A cache array of a set of SM parameters, used together with mf_atMz_cache.
Definition: EWSMcache.h:932
gslpp::complex::real
const double & real() const
Definition: gslpp_complex.cpp:53
Particle::getIndex
int getIndex() const
Definition: Particle.h:160
PVfunctions::Bfp
gslpp::complex Bfp(const double mu2, const double p2, const double m02, const double m12) const
.
Definition: PVfunctions.cpp:458
StandardModel::getMHl
virtual double getMHl() const
A get method to retrieve the Higgs mass .
Definition: StandardModel.h:765
EWSMcache::B0p_Mw2_Mw2_mh2_Mw2_cache
double B0p_Mw2_Mw2_mh2_Mw2_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:960
PVfunctions::B0
gslpp::complex B0(const double mu2, const double p2, const double m02, const double m12) const
.
Definition: PVfunctions.cpp:41
EWSMcache::Q_f
double Q_f(const Particle f) const
The charge of an SM fermion .
Definition: EWSMcache.h:268
FULLNNLO
Definition: OrderScheme.h:38
PVfunctions::B1
gslpp::complex B1(const double mu2, const double p2, const double m02, const double m12) const
.
Definition: PVfunctions.cpp:120
EWSMcache::logMTOPtoMH_cache
double logMTOPtoMH_cache[3]
A cache of .
Definition: EWSMcache.h:913
EWSMcache::FlagCacheInEWSMcache
bool FlagCacheInEWSMcache
A flag for caching (true by default).
Definition: EWSMcache.h:885
EWSMcache::B0_Mw2_Mw2_Mz2_Mw2_cache
double B0_Mw2_Mw2_Mz2_Mw2_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:950
StandardModel::getAle
double getAle() const
A get method to retrieve the fine-structure constant .
Definition: StandardModel.h:745
EWSMcache::C0_Mz2_0_Mw2_0_cache
double C0_Mz2_0_Mw2_0_cache[4]
A cache of a PV function.
Definition: EWSMcache.h:983
EWSMcache::Bf_Mz2_0_mfprime2_mf2_cache
double Bf_Mz2_0_mfprime2_mf2_cache[6][5]
A cache of a PV function.
Definition: EWSMcache.h:975
QCD::MU
Definition: QCD.h:314
QCD::lepton
lepton
An enum type for leptons.
Definition: QCD.h:310
EWSMcache::Clausen
const ClausenFunctions Clausen
An object of type ClausenFunctions.
Definition: EWSMcache.h:889
StandardModel::getLeptons
Particle getLeptons(const QCD::lepton p) const
A get method to retrieve the member object of a lepton.
Definition: StandardModel.h:709