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_i a reference to an object of type StandardModel

Definition at line 15 of file EWSMcache.cpp.

: SM(SM_i), PV(true)
{
    FlagDebug = false;
    FlagCacheInEWSMcache = true; // use caches in the current class
    //FlagCacheInEWSMcache = false;// do not use caches in the current class (for test)
 
    log2 = log(2.0);
 
    /* zeta functions */
    zeta2 = gsl_sf_zeta_int(2);
    zeta3 = gsl_sf_zeta_int(3);
    zeta4 = gsl_sf_zeta_int(4);
    zeta5 = gsl_sf_zeta_int(5);
 
    /* Constants for three-loop contribution */
    double Cl2_Pi_3 = Clausen.Cl2(M_PI / 3.0);
    S2 = 4.0 / 9.0 / sqrt(3.0) * Cl2_Pi_3;
    D3 = 6.0 * zeta3 - 15.0 / 4.0 * zeta4 - 6.0 * Cl2_Pi_3*Cl2_Pi_3;
    B4 = -1.76280008707377;
    //double Li4_1_2 = ??;
    //B4 = 16.0*Li4_1_2 - 4.0*zeta2*log2*log2 + 2.0/3.0*pow(log2,4.0) - 13.0/2.0*zeta4;
 
    // Initializations of the cache
    for (int i = 0; i < 12; ++i) {
        mf_atMz_cache[i] = 0.0;
        for (int j = 0; j < StandardModel::NumSMParamsForEWPO; ++j)
            mf_atMz_params_cache[i][j] = 0.0;
    }
}

Member Function Documentation

◆ A0_Mw2_mh2()

double EWSMcache::A0_Mw2_mh2 ( const double Mw_i ) const

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

Returns: \(A_0(m_h^2)\) with \(\mu=M_W\)

Definition at line 247 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {Mw_i, SM.getMHl()};
 
    if (CacheCheck(A0_Mw2_mh2_cache, NumPar, params))
        return A0_Mw2_mh2_cache[NumPar];
    else {
        double newResult = PV.A0(Mw_i*Mw_i, SM.getMHl() * SM.getMHl());
        newCacheForDouble(A0_Mw2_mh2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ A0_Mw2_Mw2()

double EWSMcache::A0_Mw2_Mw2 ( const double Mw_i ) const

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

Returns: \(A_0(M_W^2)\) with \(\mu=M_W\)

Definition at line 275 of file EWSMcache.cpp.

{
    int NumPar = 1;
    double params[] = {Mw_i};
 
    if (CacheCheck(A0_Mw2_Mw2_cache, NumPar, params))
        return A0_Mw2_Mw2_cache[NumPar];
    else {
        double newResult = PV.A0(Mw_i*Mw_i, Mw_i * Mw_i);
        newCacheForDouble(A0_Mw2_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ A0_Mw2_Mz2()

double EWSMcache::A0_Mw2_Mz2 ( const double Mw_i ) const

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 233 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {Mw_i, SM.getMz()};
 
    if (CacheCheck(A0_Mw2_Mz2_cache, NumPar, params))
        return A0_Mw2_Mz2_cache[NumPar];
    else {
        double newResult = PV.A0(Mw_i*Mw_i, SM.getMz() * SM.getMz());
        newCacheForDouble(A0_Mw2_Mz2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ 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 219 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), SM.getMHl()};
 
    if (CacheCheck(A0_Mz2_mh2_cache, NumPar, params))
        return A0_Mz2_mh2_cache[NumPar];
    else {
        double newResult = PV.A0(SM.getMz() * SM.getMz(), SM.getMHl() * SM.getMHl());
        newCacheForDouble(A0_Mz2_mh2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ A0_Mz2_Mw2()

double EWSMcache::A0_Mz2_Mw2 ( const double Mw_i ) const

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 205 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), Mw_i};
 
    if (CacheCheck(A0_Mz2_Mw2_cache, NumPar, params))
        return A0_Mz2_Mw2_cache[NumPar];
    else {
        double newResult = PV.A0(SM.getMz() * SM.getMz(), Mw_i * Mw_i);
        newCacheForDouble(A0_Mz2_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ 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 261 of file EWSMcache.cpp.

{
    int NumPar = 1;
    double params[] = {SM.getMz()};
 
    if (CacheCheck(A0_Mz2_Mz2_cache, NumPar, params))
        return A0_Mz2_Mz2_cache[NumPar];
    else {
        double newResult = PV.A0(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz());
        newCacheForDouble(A0_Mz2_Mz2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ 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] f a lepton or quark

Returns: \(a_f\)

Definition at line 301 of file EWSMcache.h.

    {
        return ( f.getIsospin());
    }

◆ Als()

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

inline

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

Parameters

[in]	mu	renormalization scale \(\mu\) in GeV
[in]	order	order in QCD (see orders)

Returns: \(\alpha_s(\mu)\) at the given order

Definition at line 366 of file EWSMcache.h.

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

◆ 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.

    {
        if (FlagDebug)
            return ( 0.1074432788); // for debug
        else
            return ( SM.Als(SM.getMtpole(), FULLNNLO));
    }

◆ B0_Mw2_Mw2_0_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

Returns: \(B_0(M_W^2;0,M_W^2)\) with \(\mu=M_W\)

Definition at line 469 of file EWSMcache.cpp.

{
    int NumPar = 1;
    double params[] = {Mw_i};
 
    if (CacheCheck(B0_Mw2_Mw2_0_Mw2_cache, NumPar, params))
        return gslpp::complex(B0_Mw2_Mw2_0_Mw2_cache[NumPar],
            B0_Mw2_Mw2_0_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0(Mw_i*Mw_i, Mw_i*Mw_i, 0.0, Mw_i * Mw_i);
        newCacheForComplex(B0_Mw2_Mw2_0_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0_Mw2_Mw2_mh2_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 454 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {Mw_i, SM.getMHl()};
 
    if (CacheCheck(B0_Mw2_Mw2_mh2_Mw2_cache, NumPar, params))
        return gslpp::complex(B0_Mw2_Mw2_mh2_Mw2_cache[NumPar],
            B0_Mw2_Mw2_mh2_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0(Mw_i*Mw_i, Mw_i*Mw_i, SM.getMHl() * SM.getMHl(), Mw_i * Mw_i);
        newCacheForComplex(B0_Mw2_Mw2_mh2_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0_Mw2_Mw2_Mz2_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 439 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {Mw_i, SM.getMz()};
 
    if (CacheCheck(B0_Mw2_Mw2_Mz2_Mw2_cache, NumPar, params))
        return gslpp::complex(B0_Mw2_Mw2_Mz2_Mw2_cache[NumPar],
            B0_Mw2_Mw2_Mz2_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0(Mw_i*Mw_i, Mw_i*Mw_i, SM.getMz() * SM.getMz(), Mw_i * Mw_i);
        newCacheForComplex(B0_Mw2_Mw2_Mz2_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0_Mw2_Mz2_Mt2_Mt2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

Returns: \(B_0(M_Z^2;m_t^2,m_t^2)\) with \(\mu=M_W\)

Definition at line 319 of file EWSMcache.cpp.

{
    int NumPar = 3;
    double params[] = {Mw_i, SM.getMz(), SM.getMtpole()};
 
    if (CacheCheck(B0_Mw2_Mz2_Mt2_Mt2_cache, NumPar, params))
        return gslpp::complex(B0_Mw2_Mz2_Mt2_Mt2_cache[NumPar],
            B0_Mw2_Mz2_Mt2_Mt2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0(Mw_i*Mw_i, SM.getMz() * SM.getMz(), SM.getMtpole() * SM.getMtpole(), SM.getMtpole() * SM.getMtpole());
        newCacheForComplex(B0_Mw2_Mz2_Mt2_Mt2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0_Mw2_Mz2_Mw2_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 424 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {Mw_i, SM.getMz()};
 
    if (CacheCheck(B0_Mw2_Mz2_Mw2_Mw2_cache, NumPar, params))
        return gslpp::complex(B0_Mw2_Mz2_Mw2_Mw2_cache[NumPar],
            B0_Mw2_Mz2_Mw2_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0(Mw_i*Mw_i, SM.getMz() * SM.getMz(), Mw_i*Mw_i, Mw_i * Mw_i);
        newCacheForComplex(B0_Mw2_Mz2_Mw2_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0_Mz2_0_0_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

Returns: \(B_0(0;0,M_W^2)\) with \(\mu=M_Z\)

Definition at line 409 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), Mw_i};
 
    if (CacheCheck(B0_Mz2_0_0_Mw2_cache, NumPar, params))
        return gslpp::complex(B0_Mz2_0_0_Mw2_cache[NumPar],
            B0_Mz2_0_0_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), 0.0, 0.0, Mw_i * Mw_i);
        newCacheForComplex(B0_Mz2_0_0_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0_Mz2_0_mh2_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 304 of file EWSMcache.cpp.

{
    int NumPar = 3;
    double params[] = {SM.getMz(), SM.getMHl(), Mw_i};
 
    if (CacheCheck(B0_Mz2_0_mh2_Mw2_cache, NumPar, params))
        return gslpp::complex(B0_Mz2_0_mh2_Mw2_cache[NumPar],
            B0_Mz2_0_mh2_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), 0.0, SM.getMHl() * SM.getMHl(), Mw_i * Mw_i);
        newCacheForComplex(B0_Mz2_0_mh2_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0_Mz2_0_Mz2_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

Returns: \(B_0(0;M_Z^2,M_W^2)\) with \(\mu=M_Z\)

Definition at line 394 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), Mw_i};
 
    if (CacheCheck(B0_Mz2_0_Mz2_Mw2_cache, NumPar, params))
        return gslpp::complex(B0_Mz2_0_Mz2_Mw2_cache[NumPar],
            B0_Mz2_0_Mz2_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), 0.0, SM.getMz() * SM.getMz(), Mw_i * Mw_i);
        newCacheForComplex(B0_Mz2_0_Mz2_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0_Mz2_Mw2_0_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

Returns: \(B_0(M_W^2;0,M_W^2)\) with \(\mu=M_Z\)

Definition at line 379 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), Mw_i};
 
    if (CacheCheck(B0_Mz2_Mw2_0_Mw2_cache, NumPar, params))
        return gslpp::complex(B0_Mz2_Mw2_0_Mw2_cache[NumPar],
            B0_Mz2_Mw2_0_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), Mw_i*Mw_i, 0.0, Mw_i * Mw_i);
        newCacheForComplex(B0_Mz2_Mw2_0_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0_Mz2_Mw2_mh2_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 289 of file EWSMcache.cpp.

{
    int NumPar = 3;
    double params[] = {SM.getMz(), Mw_i, SM.getMHl()};
 
    if (CacheCheck(B0_Mz2_Mw2_mh2_Mw2_cache, NumPar, params))
        return gslpp::complex(B0_Mz2_Mw2_mh2_Mw2_cache[NumPar],
            B0_Mz2_Mw2_mh2_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), Mw_i*Mw_i, SM.getMHl() * SM.getMHl(), Mw_i * Mw_i);
        newCacheForComplex(B0_Mz2_Mw2_mh2_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0_Mz2_Mw2_Mz2_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 364 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), Mw_i};
 
    if (CacheCheck(B0_Mz2_Mw2_Mz2_Mw2_cache, NumPar, params))
        return gslpp::complex(B0_Mz2_Mw2_Mz2_Mw2_cache[NumPar],
            B0_Mz2_Mw2_Mz2_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), Mw_i*Mw_i, SM.getMz() * SM.getMz(), Mw_i * Mw_i);
        newCacheForComplex(B0_Mz2_Mw2_Mz2_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0_Mz2_Mz2_mf2_mf2()

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

A cache method.

Parameters

[in] f a lepton or quark

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

Definition at line 484 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), mf(f, SM.getMz())};
    int ind = f.getIndex();
    if (CacheCheck(B0_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params))
        return gslpp::complex(B0_Mz2_Mz2_mf2_mf2_cache[ind][NumPar],
            B0_Mz2_Mz2_mf2_mf2_cache[ind][NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), mf2(f, SM.getMz()), mf2(f, SM.getMz()));
        newCacheForComplex(B0_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params, newResult);
        return newResult;
    }
}

◆ 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 349 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), SM.getMHl()};
 
    if (CacheCheck(B0_Mz2_Mz2_mh2_Mz2_cache, NumPar, params))
        return gslpp::complex(B0_Mz2_Mz2_mh2_Mz2_cache[NumPar],
            B0_Mz2_Mz2_mh2_Mz2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), SM.getMHl() * SM.getMHl(), SM.getMz() * SM.getMz());
        newCacheForComplex(B0_Mz2_Mz2_mh2_Mz2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0_Mz2_Mz2_Mw2_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 334 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), Mw_i};
 
    if (CacheCheck(B0_Mz2_Mz2_Mw2_Mw2_cache, NumPar, params))
        return gslpp::complex(B0_Mz2_Mz2_Mw2_Mw2_cache[NumPar],
            B0_Mz2_Mz2_Mw2_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), Mw_i*Mw_i, Mw_i * Mw_i);
        newCacheForComplex(B0_Mz2_Mz2_Mw2_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0p_Mw2_Mw2_0_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 589 of file EWSMcache.cpp.

{
    int NumPar = 1;
    double params[] = {Mw_i};
 
    if (CacheCheck(B0p_Mw2_Mw2_0_Mw2_cache, NumPar, params))
        return gslpp::complex(B0p_Mw2_Mw2_0_Mw2_cache[NumPar],
            B0p_Mw2_Mw2_0_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0p(Mw_i*Mw_i, Mw_i*Mw_i, 0.0, Mw_i * Mw_i);
        newCacheForComplex(B0p_Mw2_Mw2_0_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0p_Mw2_Mw2_mh2_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 574 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {Mw_i, SM.getMHl()};
 
    if (CacheCheck(B0p_Mw2_Mw2_mh2_Mw2_cache, NumPar, params))
        return gslpp::complex(B0p_Mw2_Mw2_mh2_Mw2_cache[NumPar],
            B0p_Mw2_Mw2_mh2_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0p(Mw_i*Mw_i, Mw_i*Mw_i, SM.getMHl() * SM.getMHl(), Mw_i * Mw_i);
        newCacheForComplex(B0p_Mw2_Mw2_mh2_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0p_Mw2_Mw2_Mz2_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 559 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {Mw_i, SM.getMz()};
 
    if (CacheCheck(B0p_Mw2_Mw2_Mz2_Mw2_cache, NumPar, params))
        return gslpp::complex(B0p_Mw2_Mw2_Mz2_Mw2_cache[NumPar],
            B0p_Mw2_Mw2_Mz2_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0p(Mw_i*Mw_i, Mw_i*Mw_i, SM.getMz() * SM.getMz(), Mw_i * Mw_i);
        newCacheForComplex(B0p_Mw2_Mw2_Mz2_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0p_Mz2_0_mh2_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 499 of file EWSMcache.cpp.

{
    int NumPar = 3;
    double params[] = {SM.getMz(), SM.getMHl(), Mw_i};
 
    if (CacheCheck(B0p_Mz2_0_mh2_Mw2_cache, NumPar, params))
        return gslpp::complex(B0p_Mz2_0_mh2_Mw2_cache[NumPar],
            B0p_Mz2_0_mh2_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0p(SM.getMz() * SM.getMz(), 0.0, SM.getMHl() * SM.getMHl(), Mw_i * Mw_i);
        newCacheForComplex(B0p_Mz2_0_mh2_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0p_Mz2_0_Mz2_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 529 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), Mw_i};
 
    if (CacheCheck(B0p_Mz2_0_Mz2_Mw2_cache, NumPar, params))
        return gslpp::complex(B0p_Mz2_0_Mz2_Mw2_cache[NumPar],
            B0p_Mz2_0_Mz2_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0p(SM.getMz() * SM.getMz(), 0.0, SM.getMz() * SM.getMz(), Mw_i * Mw_i);
        newCacheForComplex(B0p_Mz2_0_Mz2_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0p_Mz2_Mz2_mf2_mf2()

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

A cache method.

Parameters

[in] f a lepton or quark

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

Definition at line 604 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), mf(f, SM.getMz())};
    int ind = f.getIndex();
    if (CacheCheck(B0p_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params))
        return gslpp::complex(B0p_Mz2_Mz2_mf2_mf2_cache[ind][NumPar],
            B0p_Mz2_Mz2_mf2_mf2_cache[ind][NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0p(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), mf2(f, SM.getMz()), mf2(f, SM.getMz()));
        newCacheForComplex(B0p_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params, newResult);
        return newResult;
    }
}

◆ 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 514 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), SM.getMHl()};
 
    if (CacheCheck(B0p_Mz2_Mz2_mh2_Mz2_cache, NumPar, params))
        return gslpp::complex(B0p_Mz2_Mz2_mh2_Mz2_cache[NumPar],
            B0p_Mz2_Mz2_mh2_Mz2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0p(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), SM.getMHl() * SM.getMHl(), SM.getMz() * SM.getMz());
        newCacheForComplex(B0p_Mz2_Mz2_mh2_Mz2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B0p_Mz2_Mz2_Mw2_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 544 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), Mw_i};
 
    if (CacheCheck(B0p_Mz2_Mz2_Mw2_Mw2_cache, NumPar, params))
        return gslpp::complex(B0p_Mz2_Mz2_Mw2_Mw2_cache[NumPar],
            B0p_Mz2_Mz2_Mw2_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.B0p(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), Mw_i*Mw_i, Mw_i * Mw_i);
        newCacheForComplex(B0p_Mz2_Mz2_Mw2_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ B1_Mz2_0_mf2_mfprime2()

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

A cache method.

Parameters

[in] gen the generation index of a lepton doublet

Returns: \(B_1(0;m_l^2,m_{l'}^2)\) with \(\mu=M_Z\)

Definition at line 619 of file EWSMcache.cpp.

{
    int NumPar = 3;
    double mymf, mymfprime;
    double Mz = SM.getMz();
    if (gen < 3) {
        mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
        mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
    } else {
        int genq = gen - 3;
        mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
        mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
    }
    double params[] = {Mz, mymf, mymfprime};
 
    if (CacheCheck(B1_Mz2_0_mf2_mfprime2_cache[gen], NumPar, params))
        return gslpp::complex(B1_Mz2_0_mf2_mfprime2_cache[gen][NumPar],
            B1_Mz2_0_mf2_mfprime2_cache[gen][NumPar + 1], false);
    else {
        double mf2 = mymf*mymf;
        double mfprime2 = mymfprime*mymfprime;
        gslpp::complex newResult = PV.B1(Mz * Mz, 0.0, mf2, mfprime2);
        newCacheForComplex(B1_Mz2_0_mf2_mfprime2_cache[gen], NumPar, params, newResult);
        return newResult;
    }
}

◆ B1_Mz2_0_mfprime2_mf2()

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

A cache method.

Parameters

[in] gen the generation index of a lepton doublet

Returns: \(B_1(0;m_{l'}^2,m_l^2)\) with \(\mu=M_Z\)

Definition at line 646 of file EWSMcache.cpp.

{
    int NumPar = 3;
    double mymf, mymfprime;
    double Mz = SM.getMz();
    if (gen < 3) {
        mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
        mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
    } else {
        int genq = gen - 3;
        mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
        mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
    }
    double params[] = {Mz, mymf, mymfprime};
 
    if (CacheCheck(B1_Mz2_0_mfprime2_mf2_cache[gen], NumPar, params))
        return gslpp::complex(B1_Mz2_0_mfprime2_mf2_cache[gen][NumPar],
            B1_Mz2_0_mfprime2_mf2_cache[gen][NumPar + 1], false);
    else {
        double mf2 = mymf*mymf;
        double mfprime2 = mymfprime*mymfprime;
        gslpp::complex newResult = PV.B1(Mz * Mz, 0.0, mfprime2, mf2);
        newCacheForComplex(B1_Mz2_0_mfprime2_mf2_cache[gen], NumPar, params, newResult);
        return newResult;
    }
}

◆ 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]	gen	the generation index of a lepton doublet
[in]	Mw_i	the \(W\)-boson mass

Returns: \(B_1(M_W^2;m_l^2,m_{l'}^2)\) with \(\mu=M_Z\)

Definition at line 673 of file EWSMcache.cpp.

{
    int NumPar = 4;
    double mymf, mymfprime;
    double Mz = SM.getMz();
    if (gen < 3) {
        mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
        mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
    } else {
        int genq = gen - 3;
        mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
        mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
    }
    double params[] = {Mz, Mw_i, mymf, mymfprime};
 
    if (CacheCheck(B1_Mz2_Mw2_mf2_mfprime2_cache[gen], NumPar, params))
        return gslpp::complex(B1_Mz2_Mw2_mf2_mfprime2_cache[gen][NumPar],
            B1_Mz2_Mw2_mf2_mfprime2_cache[gen][NumPar + 1], false);
    else {
        double mf2 = mymf*mymf;
        double mfprime2 = mymfprime*mymfprime;
        gslpp::complex newResult = PV.B1(Mz * Mz, Mw_i*Mw_i, mf2, mfprime2);
        newCacheForComplex(B1_Mz2_Mw2_mf2_mfprime2_cache[gen], NumPar, params, newResult);
        return newResult;
    }
}

◆ 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]	gen	the generation index of a lepton doublet
[in]	Mw_i	the \(W\)-boson mass

Returns: \(B_1(M_W^2;m_{l'}^2,m_l^2)\) with \(\mu=M_Z\)

Definition at line 700 of file EWSMcache.cpp.

{
    int NumPar = 4;
    double mymf, mymfprime;
    double Mz = SM.getMz();
    if (gen < 3) {
        mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
        mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
    } else {
        int genq = gen - 3;
        mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
        mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
    }
    double params[] = {Mz, Mw_i, mymf, mymfprime};
 
    if (CacheCheck(B1_Mz2_Mw2_mfprime2_mf2_cache[gen], NumPar, params))
        return gslpp::complex(B1_Mz2_Mw2_mfprime2_mf2_cache[gen][NumPar],
            B1_Mz2_Mw2_mfprime2_mf2_cache[gen][NumPar + 1], false);
    else {
        double mf2 = mymf*mymf;
        double mfprime2 = mymfprime*mymfprime;
        gslpp::complex newResult = PV.B1(Mz * Mz, Mw_i*Mw_i, mfprime2, mf2);
        newCacheForComplex(B1_Mz2_Mw2_mfprime2_mf2_cache[gen], NumPar, params, newResult);
        return newResult;
    }
}

◆ 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]	gen	the generation index of a lepton doublet
[in]	Mw_i	the \(W\)-boson mass

Returns: \(B_{1p}(M_W^2;m_l^2,m_{l'}^2)\) with \(\mu=M_W\)

Definition at line 727 of file EWSMcache.cpp.

{
    int NumPar = 3;
    double mymf, mymfprime;
    double Mz = SM.getMz();
    if (gen < 3) {
        mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
        mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
    } else {
        int genq = gen - 3;
        mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
        mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
    }
    double params[] = {Mw_i, mymf, mymfprime};
 
    if (CacheCheck(B1p_Mw2_Mw2_mf2_mfprime2_cache[gen], NumPar, params))
        return gslpp::complex(B1p_Mw2_Mw2_mf2_mfprime2_cache[gen][NumPar],
            B1p_Mw2_Mw2_mf2_mfprime2_cache[gen][NumPar + 1], false);
    else {
        double mf2 = mymf*mymf;
        double mfprime2 = mymfprime*mymfprime;
        gslpp::complex newResult = PV.B1p(Mw_i*Mw_i, Mw_i*Mw_i, mf2, mfprime2);
        newCacheForComplex(B1p_Mw2_Mw2_mf2_mfprime2_cache[gen], NumPar, params, newResult);
        return newResult;
    }
}

◆ 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]	gen	the generation index of a lepton doublet
[in]	Mw_i	the \(W\)-boson mass

Returns: \(B_{1p}(M_W^2;m_{l'}^2,m_l^2)\) with \(\mu=M_W\)

Definition at line 754 of file EWSMcache.cpp.

{
    int NumPar = 3;
    double mymf, mymfprime;
    double Mz = SM.getMz();
    if (gen < 3) {
        mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
        mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
    } else {
        int genq = gen - 3;
        mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
        mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
    }
    double params[] = {Mw_i, mymf, mymfprime};
 
    if (CacheCheck(B1p_Mw2_Mw2_mfprime2_mf2_cache[gen], NumPar, params))
        return gslpp::complex(B1p_Mw2_Mw2_mfprime2_mf2_cache[gen][NumPar],
            B1p_Mw2_Mw2_mfprime2_mf2_cache[gen][NumPar + 1], false);
    else {
        double mf2 = mymf*mymf;
        double mfprime2 = mymfprime*mymfprime;
        gslpp::complex newResult = PV.B1p(Mw_i*Mw_i, Mw_i*Mw_i, mfprime2, mf2);
        newCacheForComplex(B1p_Mw2_Mw2_mfprime2_mf2_cache[gen], NumPar, params, newResult);
        return newResult;
    }
}

◆ 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]	gen	the generation index of a lepton doublet
[in]	Mw_i	the \(W\)-boson mass

Returns: \(B_f(M_W^2;m_{l'}^2,m_l^2)\) with \(\mu=M_W\)

Definition at line 866 of file EWSMcache.cpp.

{
    int NumPar = 3;
    double mymf, mymfprime;
    double Mz = SM.getMz();
    if (gen < 3) {
        mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
        mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
    } else {
        int genq = gen - 3;
        mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
        mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
    }
    double params[] = {Mw_i, mymf, mymfprime};
 
    if (CacheCheck(Bf_Mw2_Mw2_mfprime2_mf2_cache[gen], NumPar, params))
        return gslpp::complex(Bf_Mw2_Mw2_mfprime2_mf2_cache[gen][NumPar],
            Bf_Mw2_Mw2_mfprime2_mf2_cache[gen][NumPar + 1], false);
    else {
        double mf2 = mymf*mymf;
        double mfprime2 = mymfprime*mymfprime;
        gslpp::complex newResult = PV.Bf(Mw_i*Mw_i, Mw_i*Mw_i, mfprime2, mf2);
        newCacheForComplex(Bf_Mw2_Mw2_mfprime2_mf2_cache[gen], NumPar, params, newResult);
        return newResult;
    }
}

◆ Bf_Mz2_0_mf2_mf2()

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

A cache method.

Parameters

[in] f a lepton or quark

Returns: \(B_f(0;m_f^2,m_f^2)\) with \(\mu=M_Z\)

Definition at line 796 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), mf(f, SM.getMz())};
    if (params[1] == 0.0)
        throw std::runtime_error("Error in EWSMcache::Bf_Mz_0_mf_mf()");
    int ind = f.getIndex();
    if (CacheCheck(Bf_Mz2_0_mf2_mf2_cache[ind], NumPar, params))
        return gslpp::complex(Bf_Mz2_0_mf2_mf2_cache[ind][NumPar],
            Bf_Mz2_0_mf2_mf2_cache[ind][NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.Bf(SM.getMz() * SM.getMz(), 0.0, mf2(f, SM.getMz()), mf2(f, SM.getMz()));
        newCacheForComplex(Bf_Mz2_0_mf2_mf2_cache[ind], NumPar, params, newResult);
        return newResult;
    }
}

◆ Bf_Mz2_0_mfprime2_mf2()

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

A cache method.

Parameters

[in] gen the generation index of a lepton doublet

Returns: \(B_f(0;m_{l'}^2,m_l^2)\) with \(\mu=M_Z\)

Definition at line 839 of file EWSMcache.cpp.

{
    int NumPar = 3;
    double mymf, mymfprime;
    double Mz = SM.getMz();
    if (gen < 3) {
        mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
        mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
    } else {
        int genq = gen - 3;
        mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
        mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
    }
    double params[] = {Mz, mymf, mymfprime};
 
    if (CacheCheck(Bf_Mz2_0_mfprime2_mf2_cache[gen], NumPar, params))
        return gslpp::complex(Bf_Mz2_0_mfprime2_mf2_cache[gen][NumPar],
            Bf_Mz2_0_mfprime2_mf2_cache[gen][NumPar + 1], false);
    else {
        double mf2 = mymf*mymf;
        double mfprime2 = mymfprime*mymfprime;
        gslpp::complex newResult = PV.Bf(Mz * Mz, 0.0, mfprime2, mf2);
        newCacheForComplex(Bf_Mz2_0_mfprime2_mf2_cache[gen], NumPar, params, newResult);
        return newResult;
    }
}

◆ 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]	gen	the generation index of a lepton doublet
[in]	Mw_i	the \(W\)-boson mass

Returns: \(B_f(M_W^2;m_{l'}^2,m_l^2)\) with \(\mu=M_Z\)

Definition at line 813 of file EWSMcache.cpp.

{
    int NumPar = 4;
    double mymf, mymfprime;
    double Mz = SM.getMz();
    if (gen < 3) {
        mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
        mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
    } else {
        int genq = gen - 3;
        mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
        mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
    }
    double params[] = {Mz, Mw_i, mymf, mymfprime};
    if (CacheCheck(Bf_Mz2_Mw2_mfprime2_mf2_cache[gen], NumPar, params))
        return gslpp::complex(Bf_Mz2_Mw2_mfprime2_mf2_cache[gen][NumPar],
            Bf_Mz2_Mw2_mfprime2_mf2_cache[gen][NumPar + 1], false);
    else {
        double mf2 = mymf*mymf;
        double mfprime2 = mymfprime*mymfprime;
        gslpp::complex newResult = PV.Bf(Mz * Mz, Mw_i*Mw_i, mfprime2, mf2);
        newCacheForComplex(Bf_Mz2_Mw2_mfprime2_mf2_cache[gen], NumPar, params, newResult);
        return newResult;
    }
}

◆ Bf_Mz2_Mz2_mf2_mf2()

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

A cache method.

Parameters

[in] f a lepton or quark

Returns: \(B_f(M_Z^2;m_f^2,m_f^2)\) with \(\mu=M_Z\)

Definition at line 781 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), mf(f, SM.getMz())};
    int ind = f.getIndex();
    if (CacheCheck(Bf_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params))
        return gslpp::complex(Bf_Mz2_Mz2_mf2_mf2_cache[ind][NumPar],
            Bf_Mz2_Mz2_mf2_mf2_cache[ind][NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.Bf(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), mf2(f, SM.getMz()), mf2(f, SM.getMz()));
        newCacheForComplex(Bf_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params, newResult);
        return newResult;
    }
}

◆ 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]	gen	the generation index of a lepton doublet
[in]	Mw_i	the \(W\)-boson mass

Returns: \(B_{fp}(M_W^2;m_{l'}^2,m_l^2)\) with \(\mu=M_W\)

Definition at line 908 of file EWSMcache.cpp.

{
    int NumPar = 3;
    double mymf, mymfprime;
    double Mz = SM.getMz();
    if (gen < 3) {
        mymf = mf(SM.getLeptons((QCD::lepton) (2 * gen)), Mz);
        mymfprime = mf(SM.getLeptons((QCD::lepton) (2 * gen + 1)), Mz);
    } else {
        int genq = gen - 3;
        mymf = mf(SM.getQuarks((QCD::quark) (2 * genq)), Mz);
        mymfprime = mf(SM.getQuarks((QCD::quark) (2 * genq + 1)), Mz);
    }
    double params[] = {Mw_i, mymf, mymfprime};
 
    if (CacheCheck(Bfp_Mw2_Mw2_mfprime2_mf2_cache[gen], NumPar, params))
        return gslpp::complex(Bfp_Mw2_Mw2_mfprime2_mf2_cache[gen][NumPar],
            Bfp_Mw2_Mw2_mfprime2_mf2_cache[gen][NumPar + 1], false);
    else {
        double mf2 = mymf*mymf;
        double mfprime2 = mymfprime*mymfprime;
        gslpp::complex newResult = PV.Bfp(Mw_i*Mw_i, Mw_i*Mw_i, mfprime2, mf2);
        newCacheForComplex(Bfp_Mw2_Mw2_mfprime2_mf2_cache[gen], NumPar, params, newResult);
        return newResult;
    }
}

◆ Bfp_Mz2_Mz2_mf2_mf2()

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

A cache method.

Parameters

[in] f a lepton or quark

Returns: \(B_{fp}(M_Z^2;m_f^2,m_f^2)\) with \(\mu=M_Z\)

Definition at line 893 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), mf(f, SM.getMz())};
    int ind = f.getIndex();
    if (CacheCheck(Bfp_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params))
        return gslpp::complex(Bfp_Mz2_Mz2_mf2_mf2_cache[ind][NumPar],
            Bfp_Mz2_Mz2_mf2_mf2_cache[ind][NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.Bfp(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), mf2(f, SM.getMz()), mf2(f, SM.getMz()));
        newCacheForComplex(Bfp_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params, newResult);
        return newResult;
    }
}

◆ C0_Mw2_0_Mz2_0()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 1010 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {Mw_i, SM.getMz()};
 
    if (CacheCheck(C0_Mw2_0_Mz2_0_cache, NumPar, params))
        return gslpp::complex(C0_Mw2_0_Mz2_0_cache[NumPar],
            C0_Mw2_0_Mz2_0_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.C0(Mw_i*Mw_i, 0.0, SM.getMz() * SM.getMz(), 0.0);
        newCacheForComplex(C0_Mw2_0_Mz2_0_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ C0_Mw2_Mw2_0_Mz2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 995 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {Mw_i, SM.getMz()};
 
    if (CacheCheck(C0_Mw2_Mw2_0_Mz2_cache, NumPar, params))
        return gslpp::complex(C0_Mw2_Mw2_0_Mz2_cache[NumPar],
            C0_Mw2_Mw2_0_Mz2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.C0(Mw_i*Mw_i, Mw_i*Mw_i, 0.0, SM.getMz() * SM.getMz());
        newCacheForComplex(C0_Mw2_Mw2_0_Mz2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ C0_Mz2_0_Mw2_0()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 965 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), Mw_i};
 
    if (CacheCheck(C0_Mz2_0_Mw2_0_cache, NumPar, params))
        return gslpp::complex(C0_Mz2_0_Mw2_0_cache[NumPar],
            C0_Mz2_0_Mw2_0_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.C0(SM.getMz() * SM.getMz(), 0.0, Mw_i*Mw_i, 0.0);
        newCacheForComplex(C0_Mz2_0_Mw2_0_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ 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 1025 of file EWSMcache.cpp.

{
    int NumPar = 1;
    double params[] = {SM.getMz()};
 
    if (CacheCheck(C0_Mz2_0_Mz2_0_cache, NumPar, params))
        return gslpp::complex(C0_Mz2_0_Mz2_0_cache[NumPar],
            C0_Mz2_0_Mz2_0_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.C0(SM.getMz() * SM.getMz(), 0.0, SM.getMz() * SM.getMz(), 0.0);
        newCacheForComplex(C0_Mz2_0_Mz2_0_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ C0_Mz2_Mt2_Mw2_Mt2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

Returns: \(C_0(M_Z^2;M_t^2,M_W^2,m_t^2)\)

Definition at line 950 of file EWSMcache.cpp.

{
    int NumPar = 3;
    double params[] = {SM.getMz(), SM.getMtpole(), Mw_i};
 
    if (CacheCheck(C0_Mz2_Mt2_Mw2_Mt2_cache, NumPar, params))
        return gslpp::complex(C0_Mz2_Mt2_Mw2_Mt2_cache[NumPar],
            C0_Mz2_Mt2_Mw2_Mt2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.C0(SM.getMz() * SM.getMz(), SM.getMtpole() * SM.getMtpole(), Mw_i*Mw_i, SM.getMtpole() * SM.getMtpole());
        newCacheForComplex(C0_Mz2_Mt2_Mw2_Mt2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ C0_Mz2_Mw2_0_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

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

Definition at line 980 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), Mw_i};
 
    if (CacheCheck(C0_Mz2_Mw2_0_Mw2_cache, NumPar, params))
        return gslpp::complex(C0_Mz2_Mw2_0_Mw2_cache[NumPar],
            C0_Mz2_Mw2_0_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.C0(SM.getMz() * SM.getMz(), Mw_i*Mw_i, 0.0, Mw_i * Mw_i);
        newCacheForComplex(C0_Mz2_Mw2_0_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ C0_Mz2_Mw2_Mt2_Mw2()

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

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

Returns: \(C_0(M_Z^2;M_W^2,m_t^2,M_W^2)\)

Definition at line 935 of file EWSMcache.cpp.

{
    int NumPar = 3;
    double params[] = {SM.getMz(), Mw_i, SM.getMtpole()};
 
    if (CacheCheck(C0_Mz2_Mw2_Mt2_Mw2_cache, NumPar, params))
        return gslpp::complex(C0_Mz2_Mw2_Mt2_Mw2_cache[NumPar],
            C0_Mz2_Mw2_Mt2_Mw2_cache[NumPar + 1], false);
    else {
        gslpp::complex newResult = PV.C0(SM.getMz() * SM.getMz(), Mw_i*Mw_i, SM.getMtpole() * SM.getMtpole(), Mw_i * Mw_i);
        newCacheForComplex(C0_Mz2_Mw2_Mt2_Mw2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ 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]	cache	a cache of the parameters to be checked
[in]	NumPar	the number of the parameters to be checked
[in]	params	the 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.

    {
        if (!FlagCacheInEWSMcache) return false;
        bool bCache = true;
        for (int i = 0; i < NumPar; ++i)
            bCache &= (params[i] == cache[i]);
        return bCache;
    }

◆ 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.

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

◆ 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.

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

◆ 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.

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

◆ 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.

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

◆ delta_f()

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

inline

\(\delta_f = v_f-a_f\).

Parameters

[in]	f	a lepton or quark
[in]	Mw_i	the \(W\)-boson mass

Returns: \(\delta_f\)

Definition at line 323 of file EWSMcache.h.

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

◆ 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_i the \(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.

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

◆ 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.

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

◆ 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.

    {
        return B4;
    }

◆ 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.

    {
        return Clausen;
    }

◆ 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.

    {
        return D3;
    }

◆ 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.

    {
        return log2;
    }

◆ 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.

    {
        return PolyLog;
    }

◆ 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.

    {
        return PV;
    }

◆ 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.

    {
        return S2;
    }

◆ getSM()

const StandardModel & EWSMcache::getSM ( ) const

inline

Returns

Definition at line 56 of file EWSMcache.h.

    {
        return SM;
    }

◆ 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.

    {
        return zeta2;
    }

◆ 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.

    {
        return zeta3;
    }

◆ 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.

    {
        return zeta4;
    }

◆ 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.

    {
        return zeta5;
    }

◆ 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.

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

◆ I3_f()

double EWSMcache::I3_f ( const Particle f ) const

inline

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

Parameters

[in] f a lepton or quark

Returns: \(I_3^f\)

Definition at line 278 of file EWSMcache.h.

    {
        return f.getIsospin();
    }

◆ Li2_MW2toMTOP2()

double EWSMcache::Li2_MW2toMTOP2 ( const double Mw_i ) const

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

Returns: \(\mathrm{Li}_2(M_W^2/m_t^2)\)

Definition at line 162 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {Mw_i, SM.getMtpole()};
 
    if (CacheCheck(Li2_MW2toMTOP2_cache, NumPar, params))
        return Li2_MW2toMTOP2_cache[NumPar];
    else {
        double newResult = PolyLog.Li2(Mw_i * Mw_i / SM.getMtpole() / SM.getMtpole()).real();
        newCacheForDouble(Li2_MW2toMTOP2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ Li3_for_F1()

double EWSMcache::Li3_for_F1 ( const double Mw_i ) const

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

Returns: \(\mathrm{Li}_3(-M_W^2/m_t^2/(1-M_W^2/m_t^2))\)

Definition at line 190 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {Mw_i, SM.getMtpole()};
 
    if (CacheCheck(Li3_for_F1_cache, NumPar, params))
        return Li3_for_F1_cache[NumPar];
    else {
        double tmp = Mw_i * Mw_i / SM.getMtpole() / SM.getMtpole();
        double newResult = PolyLog.Li3(-tmp / (1.0 - tmp));
        newCacheForDouble(Li3_for_F1_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ Li3_MW2toMTOP2()

double EWSMcache::Li3_MW2toMTOP2 ( const double Mw_i ) const

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

Returns: \(\mathrm{Li}_3(M_W^2/m_t^2)\)

Definition at line 176 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {Mw_i, SM.getMtpole()};
 
    if (CacheCheck(Li3_MW2toMTOP2_cache, NumPar, params))
        return Li3_MW2toMTOP2_cache[NumPar];
    else {
        double newResult = PolyLog.Li3(Mw_i * Mw_i / SM.getMtpole() / SM.getMtpole());
        newCacheForDouble(Li3_MW2toMTOP2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ log_cW2()

double EWSMcache::log_cW2 ( const double Mw_i ) const

A cache method.

Parameters

[in] Mw_i the \(W\)-boson mass

Returns: \(\ln c_W^2\)

Definition at line 148 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), Mw_i};
 
    if (CacheCheck(log_cW2_cache, NumPar, params))
        return log_cW2_cache[NumPar];
    else {
        double newResult = log(SM.cW2(Mw_i));
        newCacheForDouble(log_cW2_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ logMTOPtoMH()

double EWSMcache::logMTOPtoMH ( ) const

A cache method.

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

Definition at line 134 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMtpole(), SM.getMHl()};
 
    if (CacheCheck(logMTOPtoMH_cache, NumPar, params))
        return logMTOPtoMH_cache[NumPar];
    else {
        double newResult = log(SM.getMtpole() / SM.getMHl());
        newCacheForDouble(logMTOPtoMH_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ logMZtoME()

double EWSMcache::logMZtoME ( ) const

A cache method.

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

Definition at line 78 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), mf(SM.getLeptons(SM.ELECTRON))};
 
    if (CacheCheck(logMZtoME_cache, NumPar, params))
        return logMZtoME_cache[NumPar];
    else {
        double newResult = log(SM.getMz() / mf(SM.getLeptons(SM.ELECTRON)));
        newCacheForDouble(logMZtoME_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ logMZtoMMU()

double EWSMcache::logMZtoMMU ( ) const

A cache method.

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

Definition at line 92 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), mf(SM.getLeptons(SM.MU))};
 
    if (CacheCheck(logMZtoMMU_cache, NumPar, params))
        return logMZtoMMU_cache[NumPar];
    else {
        double newResult = log(SM.getMz() / mf(SM.getLeptons(SM.MU)));
        newCacheForDouble(logMZtoMMU_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ logMZtoMTAU()

double EWSMcache::logMZtoMTAU ( ) const

A cache method.

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

Definition at line 106 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), mf(SM.getLeptons(SM.TAU))};
 
    if (CacheCheck(logMZtoMTAU_cache, NumPar, params))
        return logMZtoMTAU_cache[NumPar];
    else {
        double newResult = log(SM.getMz() / mf(SM.getLeptons(SM.TAU)));
        newCacheForDouble(logMZtoMTAU_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ logMZtoMTOP()

double EWSMcache::logMZtoMTOP ( ) const

A cache method.

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

Definition at line 120 of file EWSMcache.cpp.

{
    int NumPar = 2;
    double params[] = {SM.getMz(), SM.getMtpole()};
 
    if (CacheCheck(logMZtoMTOP_cache, NumPar, params))
        return logMZtoMTOP_cache[NumPar];
    else {
        double newResult = log(SM.getMz() / SM.getMtpole());
        newCacheForDouble(logMZtoMTOP_cache, NumPar, params, newResult);
        return newResult;
    }
}

◆ 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.

    {
        gsl_complex OneMinusE2Iphi = gsl_complex_rect(1.0 - cos(2.0 * Phi_QCD2()),
                -sin(2.0 * Phi_QCD2()));
        gsl_complex OneMinusE4Iphi = gsl_complex_rect(1.0 - cos(4.0 * Phi_QCD2()),
                -sin(4.0 * Phi_QCD2()));
        return ( GSL_REAL(gsl_complex_log(OneMinusE2Iphi))
                - 2.0 * GSL_REAL(gsl_complex_log(OneMinusE4Iphi)));
    }

◆ 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]	f	a lepton or quark
[in]	mu	renormalization scale
[in]	order	order 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.

{
    if (f.is("TOP"))
        return SM.getMtpole(); // the pole mass
    else if (f.is("BOTTOM") && !FlagDebug) {
        /* These codes are slow and not effective. */
        //if (mu == SM.getMz()) {
        //    if (FlagCacheInEWSMcache && order == FULLNNLO)
        //        if (SM.checkSMparams(mf_atMz_params_cache[f.getIndex()]))
        //            return mf_atMz_cache[f.getIndex()];
        //    mf_atMz_cache[f.getIndex()] = SM.Mrun(mu, f.getMass_scale(), f.getMass(), order);
        //    return mf_atMz_cache[f.getIndex()];
        //}
        return SM.Mrun(mu, f.getMass_scale(), f.getMass(), QCD::BOTTOM, order);
    } else if (f.is("CHARM") && !FlagDebug) {
        return SM.Mrun(mu, f.getMass_scale(), f.getMass(), QCD::CHARM, order);
    } else if (f.is("STRANGE") && !FlagDebug) {
        return SM.Mrun(mu, f.getMass_scale(), f.getMass(), QCD::STRANGE, order);
    } else if (f.is("UP") && !FlagDebug) {
        return SM.Mrun(mu, f.getMass_scale(), f.getMass(), QCD::UP, order);
    } else if (f.is("DOWN") && !FlagDebug) {
        return SM.Mrun(mu, f.getMass_scale(), f.getMass(), QCD::DOWN, order);    
    } else
        return f.getMass();
}

◆ 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]	f	a lepton or quark
[in]	mu	renormalization scale
[in]	order	order 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.

    {
        double mf1 = mf(f, mu, order);
        return ( mf1 * mf1);
    }

◆ 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]	cache	a cache of the parameters and the quantity
[in]	NumPar	the number of the parameters
[in]	params	an array of the parameters
[in]	newResult	the 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.

    {
        if (!FlagCacheInEWSMcache) return;
        for (int i = 0; i < NumPar; ++i)
            cache[i] = params[i];
        cache[NumPar] = newResult.real();
        cache[NumPar + 1] = newResult.imag();
    }

◆ 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]	cache	a cache of the parameters and the quantity
[in]	NumPar	the number of the parameters
[in]	params	an array of the parameters
[in]	newResult	the 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.

    {
        if (!FlagCacheInEWSMcache) return;
        for (int i = 0; i < NumPar; ++i)
            cache[i] = params[i];
        cache[NumPar] = newResult;
    }

◆ 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.

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

◆ Q_f()

double EWSMcache::Q_f ( const Particle f ) const

inline

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

Parameters

[in] f a lepton or quark

Returns: \(Q_f\)

Definition at line 268 of file EWSMcache.h.

    {
        return f.getCharge();
    }

◆ 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] FlagCacheInEWSMcache a boolean flag for caching

See also: the description of the StandardModel flags

Definition at line 83 of file EWSMcache.h.

    {
        this->FlagCacheInEWSMcache = FlagCacheInEWSMcache;
    }

◆ 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] FlagDebug a boolean flag for debugging

Definition at line 68 of file EWSMcache.h.

    {
        this->FlagDebug = FlagDebug;
    }

◆ sigma_f()

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

inline

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

Parameters

[in]	f	a lepton or quark
[in]	Mw_i	the \(W\)-boson mass

Returns: \(\sigma_f\)

Definition at line 312 of file EWSMcache.h.

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

◆ 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]	f	a lepton or quark
[in]	Mw_i	the \(W\)-boson mass

Returns: \(v_f\)

Definition at line 290 of file EWSMcache.h.

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

◆ Xt_alpha()

double EWSMcache::Xt_alpha ( const double Mw_i ) const

inline

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

Parameters

[in] Mw_i the \(W\)-boson mass

Returns: \(X_t^{\alpha}\)

See also: Xt_GF() and f_AlphaToGF()

Definition at line 355 of file EWSMcache.h.

    {
        return ( Xt_GF() / f_AlphaToGF(Mw_i));
    }

◆ 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.

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

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:

Detailed Description

Public Member Functions

Private Member Functions

Private Attributes

Constructor & Destructor Documentation

◆ EWSMcache()

Member Function Documentation

◆ A0_Mw2_mh2()

◆ A0_Mw2_Mw2()

◆ A0_Mw2_Mz2()

◆ A0_Mz2_mh2()

◆ A0_Mz2_Mw2()

◆ A0_Mz2_Mz2()

◆ a_f()

◆ Als()

◆ alsMt()

◆ B0_Mw2_Mw2_0_Mw2()

◆ B0_Mw2_Mw2_mh2_Mw2()

◆ B0_Mw2_Mw2_Mz2_Mw2()

◆ B0_Mw2_Mz2_Mt2_Mt2()

◆ B0_Mw2_Mz2_Mw2_Mw2()

◆ B0_Mz2_0_0_Mw2()

◆ B0_Mz2_0_mh2_Mw2()

◆ B0_Mz2_0_Mz2_Mw2()

◆ B0_Mz2_Mw2_0_Mw2()

◆ B0_Mz2_Mw2_mh2_Mw2()

◆ B0_Mz2_Mw2_Mz2_Mw2()

◆ B0_Mz2_Mz2_mf2_mf2()

◆ B0_Mz2_Mz2_mh2_Mz2()

◆ B0_Mz2_Mz2_Mw2_Mw2()

◆ B0p_Mw2_Mw2_0_Mw2()

◆ B0p_Mw2_Mw2_mh2_Mw2()

◆ B0p_Mw2_Mw2_Mz2_Mw2()

◆ B0p_Mz2_0_mh2_Mw2()

◆ B0p_Mz2_0_Mz2_Mw2()

◆ B0p_Mz2_Mz2_mf2_mf2()

◆ B0p_Mz2_Mz2_mh2_Mz2()

◆ B0p_Mz2_Mz2_Mw2_Mw2()

◆ B1_Mz2_0_mf2_mfprime2()

◆ B1_Mz2_0_mfprime2_mf2()

◆ B1_Mz2_Mw2_mf2_mfprime2()

◆ B1_Mz2_Mw2_mfprime2_mf2()

◆ B1p_Mw2_Mw2_mf2_mfprime2()

◆ B1p_Mw2_Mw2_mfprime2_mf2()

◆ Bf_Mw2_Mw2_mfprime2_mf2()

◆ Bf_Mz2_0_mf2_mf2()

◆ Bf_Mz2_0_mfprime2_mf2()

◆ Bf_Mz2_Mw2_mfprime2_mf2()

◆ Bf_Mz2_Mz2_mf2_mf2()

◆ Bfp_Mw2_Mw2_mfprime2_mf2()

◆ Bfp_Mz2_Mz2_mf2_mf2()

◆ C0_Mw2_0_Mz2_0()

◆ C0_Mw2_Mw2_0_Mz2()

◆ C0_Mz2_0_Mw2_0()

◆ C0_Mz2_0_Mz2_0()

◆ C0_Mz2_Mt2_Mw2_Mt2()

◆ C0_Mz2_Mw2_0_Mw2()

◆ C0_Mz2_Mw2_Mt2_Mw2()

◆ CacheCheck()

◆ Cl2_2Phi()

◆ Cl2_4Phi()

◆ Cl3_2Phi()

◆ Cl3_4Phi()

◆ delta_f()

◆ f_AlphaToGF()

◆ gamma_QCD2()

◆ getB4()

◆ getClausen()

◆ getD3()

◆ getLog2()

◆ getPolyLog()

◆ getPV()

◆ getS2()

◆ getSM()

◆ getZeta2()

◆ getZeta3()

◆ getZeta4()

◆ getZeta5()

◆ h_QCD2()

◆ I3_f()