a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
EWSM_Output Class Reference

A class for testing SM radiative corrections to the EW precision obsrvables. More...

#include <EWSM_Output.h>

+ Inheritance diagram for EWSM_Output:

Detailed Description

A class for testing SM radiative corrections to the EW precision obsrvables.

Author
HEPfit Collaboration

Definition at line 22 of file EWSM_Output.h.

Public Member Functions

 EWSM_Output (const StandardModel &SM_in)
 Constructor. More...
 
void outputEachDeltaKappaZ (const double f_AlphaToGF, const double cW2overSW2, const double DeltaRho[StandardModel::orders_EW_size], const double deltaKappa_rem[StandardModel::orders_EW_size], const double DeltaRbar_rem, const bool bool_Zbb, const double taub[StandardModel::orders_EW_size], const double ZbbSubtract, const double Zgamma_EW2) const
 
void outputEachDeltaKappaZ_l (const QCD::lepton l, const double Mw_i) const
 
void outputEachDeltaKappaZ_q (const QCD::quark q, const double Mw_i) const
 
void outputEachDeltaR (const double Mw_i) const
 
void outputEachDeltaRhoZ (const double f_AlphaToGF, const double DeltaRho[StandardModel::orders_EW_size], const double deltaRho_rem[StandardModel::orders_EW_size], const double DeltaRbar_rem, const bool bool_Zbb, const double taub[StandardModel::orders_EW_size], const double ZbbSubtract) const
 
void outputEachDeltaRhoZ_l (const QCD::lepton l, const double Mw_i) const
 
void outputEachDeltaRhoZ_q (const QCD::quark q, const double Mw_i) const
 
- Public Member Functions inherited from StandardModel
virtual double A_f (const Particle f) const
 The left-right asymmetry in \(e^+e^-\to Z\to \ell \bar{\ell}\) at the \(Z\)-pole, \(\mathcal{A}_\ell\). More...
 
virtual double AFB (const Particle f) const
 
double Ale (double mu, orders order, bool Nf_thr=true) const
 The running electromagnetic coupling \(\alpha_e(\mu)\) in the \(\overline{MS}\) scheme. More...
 
double ale_OS (const double mu, orders order=FULLNLO) const
 The running electromagnetic coupling \(\alpha(\mu)\) in the on-shell scheme. More...
 
double alphaMz () const
 The electromagnetic coupling at the \(Z\)-mass scale, \(\alpha(M_Z^2)=\alpha/(1-\Delta\alpha(M_Z^2))\). More...
 
double Als (double mu, orders order=FULLNLO, bool qed_flag=false, bool Nf_thr=true) const
 The running QCD coupling \(\alpha(\mu)\) in the \(\overline{MS}\) scheme including QED corrections. More...
 
double AlsByOrder (double mu, orders order=FULLNLO, bool qed_flag=false, bool Nf_thr=true) const
 
double Alstilde5 (const double mu) const
 The value of \(\frac{\alpha_s^{\mathrm{FULLNLO}}}{4\pi}\) at any scale \(\mu\) with the number of flavours \(n_f = 4\) and full EW corrections. More...
 
double Beta_e (int nm, unsigned int nf) const
 QED beta function coefficients - eq. (36) hep-ph/0512066. More...
 
double Beta_s (int nm, unsigned int nf) const
 QCD beta function coefficients including QED corrections - eq. (36) hep-ph/0512066. More...
 
double c02 () const
 The square of the cosine of the weak mixing angle \(c_0^2\) defined without weak radiative corrections. More...
 
virtual bool CheckFlags () const
 A method to check the sanity of the set of model flags. More...
 
virtual bool CheckParameters (const std::map< std::string, double > &DPars)
 A method to check if all the mandatory parameters for StandardModel have been provided in model initialization. More...
 
bool checkSMparamsForEWPO ()
 A method to check whether the parameters relevant to the EWPO are updated. More...
 
double computeBrHtobb () const
 The Br \((H\to bb)\) in the Standard Model. More...
 
double computeBrHtocc () const
 The Br \((H\to cc)\) in the Standard Model. More...
 
double computeBrHtogaga () const
 The Br \((H\to\gamma\gamma)\) in the Standard Model. More...
 
double computeBrHtogg () const
 The Br \((H\to gg)\) in the Standard Model. More...
 
double computeBrHtomumu () const
 The Br \((H\to \mu\mu)\) in the Standard Model. More...
 
double computeBrHtoss () const
 The Br \((H\to ss)\) in the Standard Model. More...
 
double computeBrHtotautau () const
 The Br \((H\to \tau\tau)\) in the Standard Model. More...
 
double computeBrHtoWW () const
 The Br \((H\to WW)\) in the Standard Model. More...
 
double computeBrHtoZga () const
 The Br \((H\to Z\gamma)\) in the Standard Model. More...
 
double computeBrHtoZZ () const
 The Br \((H\to ZZ)\) in the Standard Model. More...
 
double computeBrHtoZZinv () const
 The Br \((H\to ZZ \to inv)\) in the Standard Model. More...
 
void ComputeDeltaR_rem (const double Mw_i, double DeltaR_rem[orders_EW_size]) const
 A method to collect \(\Delta r_{\mathrm{rem}}\) computed via subclasses. More...
 
void ComputeDeltaRho (const double Mw_i, double DeltaRho[orders_EW_size]) const
 A method to collect \(\Delta\rho\) computed via subclasses. More...
 
double computeGammaHgaga_tt () const
 The top loop contribution to \(H\to\gamma\gamma\) in the Standard Model. More...
 
double computeGammaHgaga_tW () const
 The mixed \(t-W\) loop contribution to \(H\to\gamma\gamma\) in the Standard Model. More...
 
double computeGammaHgaga_WW () const
 The \(W\) loop contribution to \(H\to\gamma\gamma\) in the Standard Model. More...
 
double computeGammaHgg_bb () const
 The bottom loop contribution to \(H\to gg\) in the Standard Model. More...
 
double computeGammaHgg_tb () const
 The top-bottom interference contribution to \(H\to gg\) in the Standard Model. More...
 
double computeGammaHgg_tt () const
 The top loop contribution to \(H\to gg\) in the Standard Model. More...
 
double computeGammaHTotal () const
 The Higgs total width in the Standard Model. More...
 
double computeGammaHZga_tt () const
 The top loop contribution to \(H\to Z\gamma\) in the Standard Model. More...
 
double computeGammaHZga_tW () const
 The mixed \(t-W\) loop contribution to \(H\to Z\gamma\) in the Standard Model. More...
 
double computeGammaHZga_WW () const
 The \(W\) loop contribution to \(H\to Z\gamma\) in the Standard Model. Currently it returns the value of tab 41 in ref. [138]. More...
 
double computeSigmaggH (const double sqrt_s) const
 The ggH cross section in the Standard Model. More...
 
double computeSigmaggH_bb (const double sqrt_s) const
 The square of the bottom-quark contribution to the ggH cross section in the Standard Model. More...
 
double computeSigmaggH_tb (const double sqrt_s) const
 The top-bottom interference contribution to the ggH cross section in the Standard Model. More...
 
double computeSigmaggH_tt (const double sqrt_s) const
 The square of the top-quark contribution to the ggH cross section in the Standard Model. More...
 
double computeSigmattH (const double sqrt_s) const
 The ttH production cross section in the Standard Model. More...
 
double computeSigmaVBF (const double sqrt_s) const
 The VBF cross section in the Standard Model. More...
 
double computeSigmaWF (const double sqrt_s) const
 The W fusion contribution \(\sigma_{WF}\) to higgs-production cross section in the Standard Model. More...
 
double computeSigmaWH (const double sqrt_s) const
 The WH production cross section in the Standard Model. More...
 
double computeSigmaZF (const double sqrt_s) const
 The Z fusion contribution \(\sigma_{ZF}\) to higgs-production cross section in the Standard Model. More...
 
double computeSigmaZH (const double sqrt_s) const
 The ZH production cross section in the Standard Model. More...
 
double computeSigmaZWF (const double sqrt_s) const
 The Z W interference fusion contribution \(\sigma_{ZWF}\) to higgs-production cross section in the Standard Model. More...
 
virtual double cW2 () const
 
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 \(c_W^2\). More...
 
double DeltaAlpha () const
 The total corrections to the electromagnetic coupling \(\alpha\) at the \(Z\)-mass scale, denoted as \(\Delta\alpha(M_Z^2)\). More...
 
double DeltaAlphaL5q () const
 The sum of the leptonic and the five-flavour hadronic corrections to the electromagnetic coupling \(\alpha\) at the \(Z\)-mass scale, denoted as \(\Delta\alpha^{\ell+5q}(M_Z^2)\). More...
 
double DeltaAlphaLepton (const double s) const
 Leptonic contribution to the electromagnetic coupling \(\alpha\), denoted as \(\Delta\alpha_{\mathrm{lept}}(s)\). More...
 
double DeltaAlphaTop (const double s) const
 Top-quark contribution to the electromagnetic coupling \(\alpha\), denoted as \(\Delta\alpha_{\mathrm{top}}(s)\). More...
 
virtual gslpp::complex deltaKappaZ_f (const Particle f) const
 Flavour non-universal vertex corrections to \(\kappa_Z^l\), denoted by \(\Delta\kappa_Z^l\). More...
 
virtual double DeltaR () const
 The SM prediction for \(\Delta r\) derived from that for the \(W\) boson mass. More...
 
virtual double DeltaRbar () const
 The SM prediction for \(\Delta \overline{r}\) derived from that for the \(W\)-boson mass. More...
 
virtual gslpp::complex deltaRhoZ_f (const Particle f) const
 Flavour non-universal vertex corrections to \(\rho_Z^l\), denoted by \(\Delta\rho_Z^l\). More...
 
virtual double epsilon1 () const
 The SM contribution to the epsilon parameter \(\varepsilon_1\). More...
 
virtual double epsilon2 () const
 The SM contribution to the epsilon parameter \(\varepsilon_2\). More...
 
virtual double epsilon3 () const
 The SM contribution to the epsilon parameter \(\varepsilon_3\). More...
 
virtual double epsilonb () const
 The SM contribution to the epsilon parameter \(\varepsilon_b\). More...
 
virtual gslpp::complex gA_f (const Particle f) const
 The effective leptonic neutral-current axial-vector coupling \(g_A^l\) in the SM. More...
 
virtual double Gamma_had () const
 The hadronic decay width of the \(Z\) boson, \(\Gamma_{h}\). More...
 
virtual double Gamma_inv () const
 The invisible partial decay width of the \(Z\) boson, \(\Gamma_{\mathrm{inv}}\). More...
 
virtual double Gamma_Z () const
 The total decay width of the \(Z\) boson, \(\Gamma_Z\). More...
 
virtual double GammaW () const
 The total width of the \(W\) boson, \(\Gamma_W\). More...
 
virtual double GammaW (const Particle fi, const Particle fj) const
 A partial decay width of the \(W\) boson decay into a SM fermion pair. More...
 
virtual double GammaZ (const Particle f) const
 The \(Z\to \ell\bar{\ell}\) partial decay width, \(\Gamma_\ell\). More...
 
double getAle () const
 A get method to retrieve the fine-structure constant \(\alpha\). More...
 
double getAlsMz () const
 A get method to access the value of \(\alpha_s(M_Z)\). More...
 
virtual double getCBd () const
 The ratio of the absolute value of the $B_d$ mixing amplitude over the Standard Model value. More...
 
virtual double getCBs () const
 The ratio of the absolute value of the $B_s$ mixing amplitude over the Standard Model value. More...
 
virtual double getCCC1 () const
 A virtual implementation for the RealWeakEFTCC class. More...
 
virtual double getCCC2 () const
 A virtual implementation for the RealWeakEFTCC class. More...
 
virtual double getCCC3 () const
 A virtual implementation for the RealWeakEFTCC class. More...
 
virtual double getCCC4 () const
 A virtual implementation for the RealWeakEFTCC class. More...
 
virtual double getCCC5 () const
 A virtual implementation for the RealWeakEFTCC class. More...
 
virtual double getCDMK () const
 The ratio of the real part of the $K$ mixing amplitude over the Standard Model value. More...
 
virtual double getCepsK () const
 The ratio of the imaginary part of the $K$ mixing amplitude over the Standard Model value. More...
 
CKM getCKM () const
 A get method to retrieve the member object of type CKM. More...
 
double getDAle5Mz () const
 A get method to retrieve the five-flavour hadronic contribution to the electromagnetic coupling, \(\Delta\alpha_{\mathrm{had}}^{(5)}(M_Z^2)\). More...
 
double getDelGammaZ () const
 A get method to retrieve the theoretical uncertainty in \(\Gamma_Z\), denoted as \(\delta\,\Gamma_Z\). More...
 
double getDelMw () const
 A get method to retrieve the theoretical uncertainty in \(M_W\), denoted as \(\delta\,M_W\). More...
 
double getDelR0b () const
 A get method to retrieve the theoretical uncertainty in \(R_b^0\), denoted as \(\delta\,R_b^0\). More...
 
double getDelR0c () const
 A get method to retrieve the theoretical uncertainty in \(R_c^0\), denoted as \(\delta\,R_c^0\). More...
 
double getDelR0l () const
 A get method to retrieve the theoretical uncertainty in \(R_l^0\), denoted as \(\delta\,R_l^0\). More...
 
double getDelSigma0H () const
 A get method to retrieve the theoretical uncertainty in \(\sigma_{Hadron}^0\), denoted as \(\delta\,\sigma_{Hadron}^0\). More...
 
double getDelSin2th_b () const
 A get method to retrieve the theoretical uncertainty in \(\sin^2\theta_{\rm eff}^{b}\), denoted as \(\delta\sin^2\theta_{\rm eff}^{b}\). More...
 
double getDelSin2th_l () const
 A get method to retrieve the theoretical uncertainty in \(\sin^2\theta_{\rm eff}^{\rm lept}\), denoted as \(\delta\sin^2\theta_{\rm eff}^{\rm lept}\). More...
 
double getDelSin2th_q () const
 A get method to retrieve the theoretical uncertainty in \(\sin^2\theta_{\rm eff}^{q\not = b,t}\), denoted as \(\delta\sin^2\theta_{\rm eff}^{q\not = b,t}\). More...
 
std::string getFlagKappaZ () const
 A method to retrieve the model flag KappaZ. More...
 
std::string getFlagMw () const
 A method to retrieve the model flag Mw. More...
 
std::string getFlagRhoZ () const
 A method to retrieve the model flag RhoZ. More...
 
const FlavourgetFlavour () const
 
double getGF () const
 A get method to retrieve the Fermi constant \(G_\mu\). More...
 
int getIterationNo () const
 
Particle getLeptons (const QCD::lepton p) const
 A get method to retrieve the member object of a lepton. More...
 
virtual StandardModelMatchinggetMatching () const
 A get method to access the member reference of type StandardModelMatching. More...
 
virtual double getMHl () const
 A get method to retrieve the Higgs mass \(m_h\). More...
 
virtual double getmq (const QCD::quark q, const double mu) const
 
double getMuw () const
 A get method to retrieve the matching scale \(\mu_W\) around the weak scale. More...
 
EWSMApproximateFormulaegetMyApproximateFormulae () const
 A get method to retrieve the member pointer of type EWSMApproximateFormulae. More...
 
EWSMcachegetMyEWSMcache () const
 A get method to retrieve the member pointer of type EWSMcache. More...
 
LeptonFlavourgetMyLeptonFlavour () const
 
EWSMOneLoopEWgetMyOneLoopEW () const
 A get method to retrieve the member pointer of type EWSMOneLoopEW,. More...
 
EWSMThreeLoopEWgetMyThreeLoopEW () const
 
EWSMThreeLoopEW2QCDgetMyThreeLoopEW2QCD () const
 
EWSMThreeLoopQCDgetMyThreeLoopQCD () const
 
EWSMTwoFermionsLEP2getMyTwoFermionsLEP2 () const
 A get method to retrieve the member pointer of type EWSMTwoFermionsLEP2. More...
 
EWSMTwoLoopEWgetMyTwoLoopEW () const
 
EWSMTwoLoopQCDgetMyTwoLoopQCD () const
 
double getMz () const
 A get method to access the mass of the \(Z\) boson \(M_Z\). More...
 
virtual double getPhiBd () const
 Half the relative phase of the $B_d$ mixing amplitude w.r.t. the Standard Model one. More...
 
virtual double getPhiBs () const
 Half the relative phase of the $B_s$ mixing amplitude w.r.t. the Standard Model one. More...
 
virtual StandardModel getTrueSM () const
 
gslpp::matrix< gslpp::complexgetUPMNS () const
 A get method to retrieve the object of the PMNS matrix. More...
 
gslpp::matrix< gslpp::complexgetVCKM () const
 A get method to retrieve the CKM matrix. More...
 
gslpp::matrix< gslpp::complexgetYd () const
 A get method to retrieve the Yukawa matrix of the down-type quarks, \(Y_d\). More...
 
gslpp::matrix< gslpp::complexgetYe () const
 A get method to retrieve the Yukawa matrix of the charged leptons, \(Y_e\). More...
 
gslpp::matrix< gslpp::complexgetYn () const
 A get method to retrieve the Yukawa matrix of the neutrinos, \(Y_\nu\). More...
 
gslpp::matrix< gslpp::complexgetYu () const
 A get method to retrieve the Yukawa matrix of the up-type quarks, \(Y_u\). More...
 
virtual gslpp::complex gV_f (const Particle f) const
 The effective leptonic neutral-current vector coupling \(g_V^l\) in the SM. More...
 
virtual bool Init (const std::map< std::string, double > &DPars)
 A method to initialize the model parameters. More...
 
virtual bool InitializeModel ()
 A method to initialize the model. More...
 
bool IsFlagNoApproximateGammaZ () const
 A method to retrieve the model flag NoApproximateGammaZ. More...
 
bool IsFlagWithoutNonUniversalVC () const
 A method to retrieve the model flag WithoutNonUniversalVC. More...
 
virtual gslpp::complex kappaZ_f (const Particle f) const
 The effective leptonic neutral-current coupling \(\kappa_Z^l\) in the SM. More...
 
virtual double LEP2AFBbottom (const double s) const
 
virtual double LEP2AFBcharm (const double s) const
 
virtual double LEP2AFBmu (const double s) const
 
virtual double LEP2AFBtau (const double s) const
 
virtual double LEP2Rbottom (const double s) const
 
virtual double LEP2Rcharm (const double s) const
 
virtual double LEP2sigmaBottom (const double s) const
 
virtual double LEP2sigmaCharm (const double s) const
 
virtual double LEP2sigmaHadron (const double s) const
 
virtual double LEP2sigmaMu (const double s) const
 
virtual double LEP2sigmaTau (const double s) const
 
virtual double Mw () const
 The SM prediction for the \(W\)-boson mass in the on-shell scheme, \(M_{W,\mathrm{SM}}\). More...
 
virtual double Mw_tree () const
 The tree-level mass of the \(W\) boson, \(M_W^{\mathrm{tree}}\). More...
 
double MwbarFromMw (const double Mw) const
 A method to convert the \(W\)-boson mass in the experimental/running-width scheme to that in the complex-pole/fixed-width scheme. More...
 
double MwFromMwbar (const double Mwbar) const
 A method to convert the \(W\)-boson mass in the complex-pole/fixed-width scheme to that in the experimental/running-width scheme. More...
 
double Mzbar () const
 The \(Z\)-boson mass \(\overline{M}_Z\) in the complex-pole/fixed-width scheme. More...
 
virtual double N_nu () const
 The number of neutrinos obtained indirectly from the measurements at the Z pole, \(N_{\nu}\). More...
 
virtual bool PostUpdate ()
 The post-update method for StandardModel. More...
 
virtual bool PreUpdate ()
 The pre-update method for StandardModel. More...
 
virtual double R0_f (const Particle f) const
 The ratio \(R_\ell^0=\Gamma(Z\to {\rm hadrons})/\Gamma(Z\to \ell^+ \ell^-)\). More...
 
virtual double R_inv () const
 The ratio of the invisible and leptonic (electron) decay widths of the \(Z\) boson, \(R_{inv}\). More...
 
virtual double rho_GammaW (const Particle fi, const Particle fj) const
 EW radiative corrections to the width of \(W \to f_i \bar{f}_j\), denoted as \(\rho^W_{ij}\). More...
 
virtual gslpp::complex rhoZ_f (const Particle f) const
 The effective leptonic neutral-current coupling \(\rho_Z^l\) in the SM. More...
 
double s02 () const
 The square of the sine of the weak mixing angle \(s_0^2\) defined without weak radiative corrections. More...
 
virtual bool setFlag (const std::string name, const bool value)
 A method to set a flag of StandardModel. More...
 
void setFlagCacheInStandardModel (bool FlagCacheInStandardModel)
 A set method to change the model flag CacheInStandardModel of StandardModel. More...
 
void setFlagNoApproximateGammaZ (bool FlagNoApproximateGammaZ)
 
bool setFlagSigmaForAFB (const bool flagSigmaForAFB_i)
 
bool setFlagSigmaForR (const bool flagSigmaForR_i)
 
virtual bool setFlagStr (const std::string name, const std::string value)
 A method to set a flag of StandardModel. More...
 
virtual double sigma0_had () const
 The hadronic cross section for \(e^+e^- \to Z \to \mathrm{hadrons}\) at the \(Z\)-pole, \(\sigma_h^0\). More...
 
virtual double sin2thetaEff (const Particle f) const
 The effective weak mixing angle \(\sin^2\theta_{\rm eff}^{\,\ell}\) for \(Z\ell\bar{\ell}\) at the the \(Z\)-mass scale. More...
 
 StandardModel ()
 The default constructor. More...
 
double sW2 () const
 
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 \(s_W^2\). More...
 
virtual bool Update (const std::map< std::string, double > &DPars)
 The update method for StandardModel. More...
 
virtual double v () const
 The Higgs vacuum expectation value. More...
 
virtual ~StandardModel ()
 The default destructor. More...
 
- Public Member Functions inherited from QCD
double AboveTh (const double mu) const
 The active flavour threshold above the scale \(\mu\) as defined in QCD::Thresholds(). More...
 
void addParameters (std::vector< std::string > params_i)
 A method to add parameters that are specific to only one set of observables. More...
 
virtual double Als (const double mu, const orders order=FULLNLO, bool Nf_thr=true) const
 
double Als4 (const double mu) const
 The value of \(\alpha_s^{\mathrm{FULLNLO}}\) at any scale \(\mu\) with the number of flavours \(n_f = 4\). More...
 
virtual double AlsByOrder (const double mu, const orders order=FULLNLO, bool Nf_thr=true) const
 
double AlsOLD (const double mu, const orders order=FULLNLO) const
 Computes the running strong coupling \(\alpha_s(\mu)\) in the \(\overline{\mathrm{MS}}\) scheme. In the cases of LO, NLO and FULLNNLO, the coupling is computed with AlsWithInit(). On the other hand, in the cases of NNLO and FULLNNLO, the coupling is computed with AlsWithLambda(). More...
 
double AlsWithInit (const double mu, const double alsi, const double mu_i, const orders order) const
 Computes the running strong coupling \(\alpha_s(\mu)\) from \(\alpha_s(\mu_i)\) in the \(\overline{\mathrm{MS}}\) scheme, where it is forbidden to across a flavour threshold in the RG running from \(\mu_i\) to \(\mu\). More...
 
double AlsWithLambda (const double mu, const orders order) const
 Computes the running strong coupling \(\alpha_s(\mu)\) in the \(\overline{\mathrm{MS}}\) scheme with the use of \(\Lambda_{\rm QCD}\). More...
 
double BelowTh (const double mu) const
 The active flavour threshold below the scale \(\mu\) as defined in QCD::Thresholds(). More...
 
double Beta0 (const double nf) const
 The \(\beta_0(n_f)\) coefficient for a certain number of flavours \(n_f\). More...
 
double Beta1 (const double nf) const
 The \(\beta_1(n_f)\) coefficient for a certain number of flavours \(n_f\). More...
 
double Beta2 (const double nf) const
 The \(\beta_2(n_f)\) coefficient for a certain number of flavours \(n_f\). More...
 
double Beta3 (const double nf) const
 The \(\beta_3(n_f)\) coefficient for a certain number of flavours \(n_f\). More...
 
void CacheShift (double cache[][5], int n) const
 A member used to manage the caching for this class. More...
 
void CacheShift (int cache[][5], int n) const
 
orders FullOrder (orders order) const
 Return the FULLORDER enum corresponding to order. More...
 
double Gamma0 (const double nf) const
 The \(\gamma_0\) coefficient used to compute the running of a mass. More...
 
double Gamma1 (const double nf) const
 The \(\gamma_1\) coefficient used to compute the running of a mass. More...
 
double Gamma2 (const double nf) const
 The \(\gamma_2\) coefficient used to compute the running of a mass. More...
 
double getAlsM () const
 A get method to access the value of \(\alpha_s(M_{\alpha_s})\). More...
 
BParameter getBBd () const
 For getting the bag parameters corresponding to the operator basis \(O_1 -O_5\) in \(\Delta b = 2\) process in the \(B_d\) meson system. More...
 
BParameter getBBs () const
 For getting the bag parameters corresponding to the operator basis \(O_1 -O_5\) in \(\Delta b = 2\) process in the \(B_s\) meson system. More...
 
BParameter getBD () const
 For getting the bag parameters corresponding to the operator basis \(O_1 -O_5\) in \(\Delta c = 2\) process in the \(D^0\) meson system. More...
 
BParameter getBK () const
 For getting the bag parameters corresponding to the operator basis \(O_1 -O_5\) in \(\Delta s = 2\) process in the \(K^0\) meson system. More...
 
BParameter getBKd1 () const
 
BParameter getBKd3 () const
 
double getCF () const
 A get method to access the Casimir factor of QCD. More...
 
double getMAls () const
 A get method to access the mass scale \(M_{\alpha_s}\) at which the strong coupling constant measurement is provided. More...
 
Meson getMesons (const QCD::meson m) const
 A get method to access a meson as an object of the type Meson. More...
 
double getMtpole () const
 A get method to access the pole mass of the top quark. More...
 
double getMub () const
 A get method to access the threshold between five- and four-flavour theory in GeV. More...
 
double getMuc () const
 A get method to access the threshold between four- and three-flavour theory in GeV. More...
 
double getMut () const
 A get method to access the threshold between six- and five-flavour theory in GeV. More...
 
double getNc () const
 A get method to access the number of colours \(N_c\). More...
 
double getOptionalParameter (std::string name) const
 A method to get parameters that are specific to only one set of observables. More...
 
Particle getQuarks (const QCD::quark q) const
 A get method to access a quark as an object of the type Particle. More...
 
std::vector< std::string > getUnknownParameters ()
 A method to get the vector of the parameters that have been specified in the configuration file but not being used. More...
 
void initializeBParameter (std::string name_i) const
 A method to initialize B Parameter and the corresponding meson. More...
 
void initializeMeson (QCD::meson meson_i) const
 A method to initialize a meson. More...
 
double logLambda (const double nf, orders order) const
 Computes \(\ln\Lambda_\mathrm{QCD}\) with nf flavours in GeV. More...
 
double Mbar2Mp (const double mbar, const orders order=FULLNNLO) const
 Converts the \(\overline{\mathrm{MS}}\) mass \(m(m)\) to the pole mass. More...
 
double Mp2Mbar (const double mp, const orders order=FULLNNLO) const
 Converts a quark pole mass to the corresponding \(\overline{\mathrm{MS}}\) mass \(m(m)\). More...
 
double Mrun (const double mu, const double m, const orders order=FULLNNLO) const
 Computes a running quark mass \(m(\mu)\) from \(m(m)\). More...
 
double Mrun (const double mu_f, const double mu_i, const double m, const orders order=FULLNNLO) const
 Runs a quark mass from \(\mu_i\) to \(\mu_f\). More...
 
double Mrun4 (const double mu_f, const double mu_i, const double m) const
 The running of a mass with the number of flavours \(n_f = 4\). More...
 
double MS2DRqmass (const double MSbar) const
 Converts a quark mass from the \(\overline{\mathrm{MS}}\) scheme to the \(\overline{\mathrm{DR}}\) scheme. More...
 
double MS2DRqmass (const double MSscale, const double MSbar) const
 Converts a quark mass from the \(\overline{\mathrm{MS}}\) scheme to the \(\overline{\mathrm{DR}}\) scheme. More...
 
double Nf (const double mu) const
 The number of active flavour at scale \(\mu\). More...
 
double NfThresholdCorrections (double mu, double M, double als, int nf, orders order) const
 Threshold corrections in matching \(\alpha_s(n_f+1)\) with \(\alpha_s(n_f)\) from eq. (34) of hep-ph/0512060. More...
 
std::string orderToString (const orders order) const
 Converts an object of the enum type "orders" to the corresponding string. More...
 
 QCD ()
 Constructor. More...
 
void setNc (double Nc)
 A set method to change the number of colours \(N_c\). More...
 
void setOptionalParameter (std::string name, double value)
 A method to set the parameter value for the parameters that are specific to only one set of observables. More...
 
double Thresholds (const int i) const
 For accessing the active flavour threshold scales. More...
 
- Public Member Functions inherited from Model
void addMissingModelParameter (const std::string &missingParameterName)
 
std::vector< std::string > getmissingModelParameters ()
 
unsigned int getMissingModelParametersCount ()
 
std::string getModelName () const
 A method to fetch the name of the model. More...
 
const double & getModelParam (std::string name) const
 
bool isModelFWC_DF2 () const
 
bool isModelGeneralTHDM () const
 
bool isModelGeorgiMachacek () const
 
bool IsModelInitialized () const
 A method to check if the model is initialized. More...
 
bool isModelLinearized () const
 
bool isModelParam (std::string name) const
 
bool isModelSUSY () const
 
bool isModelTHDM () const
 
bool isModelTHDMW () const
 
bool IsUpdateError () const
 A method to check if there was any error in the model update process. More...
 
 Model ()
 The default constructor. More...
 
void raiseMissingModelParameterCount ()
 
void setModelFWC_DF2 ()
 
void setModelGeneralTHDM ()
 
void setModelGeorgiMachacek ()
 
void setModelInitialized (bool ModelInitialized)
 A set method to fix the failure or success of the initialization of the model. More...
 
void setModelLinearized (bool linearized=true)
 
void setModelName (const std::string name)
 A method to set the name of the model. More...
 
void setModelSUSY ()
 
void setModelTHDM ()
 
void setModelTHDMW ()
 
void setSliced (bool Sliced)
 
void setUpdateError (bool UpdateError)
 A set method to fix the update status as success or failure. More...
 
virtual ~Model ()
 The default destructor. More...
 

Additional Inherited Members

- Public Types inherited from StandardModel
enum  LEP2RCs { Weak = 0, WeakBox, ISR, QEDFSR, QCDFSR, NUMofLEP2RCs }
 
enum  orders_EW { EW1 = 0, EW1QCD1, EW1QCD2, EW2, EW2QCD1, EW3, orders_EW_size }
 An enumerated type representing perturbative orders of radiative corrections to EW precision observables. More...
 
- Public Types inherited from QCD
enum  lepton { NEUTRINO_1, ELECTRON, NEUTRINO_2, MU, NEUTRINO_3, TAU, NOLEPTON }
 An enum type for leptons. More...
 
enum  meson { P_0, P_P, K_0, K_P, D_0, D_P, B_D, B_P, B_S, B_C, PHI, K_star, K_star_P, D_star_P, RHO, RHO_P, OMEGA, MESON_END }
 An enum type for mesons. More...
 
enum  quark { UP, DOWN, CHARM, STRANGE, TOP, BOTTOM }
 An enum type for quarks. More...
 
- Static Public Attributes inherited from StandardModel
static const double GeVminus2_to_nb = 389379.338
 
static const double Mw_error = 0.00001
 The target accuracy of the iterative calculation of the \(W\)-boson mass in units of GeV. More...
 
static const int NSMvars = 26
 The number of the model parameters in StandardModel. More...
 
static const int NumSMParamsForEWPO = 33
 The number of the SM parameters that are relevant to the EW precision observables. More...
 
static std::string SMvars [NSMvars]
 A string array containing the labels of the model parameters in StandardModel. More...
 
- Static Public Attributes inherited from QCD
static const int NQCDvars = 11
 The number of model parameters in QCD. More...
 
static std::string QCDvars [NQCDvars]
 An array containing the labels under which all QCD parameters are stored in a vector of ModelParameter via InputParser::ReadParameters(). More...
 
- Protected Member Functions inherited from StandardModel
double AFB_NoISR_l (const QCD::lepton l_flavor, const double s) const
 
double AFB_NoISR_q (const QCD::quark q_flavor, const double s) const
 
bool checkEWPOscheme (const std::string scheme) const
 A method to check if a given scheme name in string form is valid. More...
 
virtual void computeCKM ()
 The method to compute the CKM matrix. More...
 
virtual void computeYukawas ()
 The method to compute the Yukawa matrices. More...
 
double Delta_EWQCD (const QCD::quark q) const
 The non-factorizable EW-QCD corrections to the partial widths for \(Z\to q\bar{q}\), denoted as \(\Delta_{\mathrm{EW/QCD}}\). More...
 
double getIntegrand_AFBnumeratorWithISR_bottom133 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_bottom167 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_bottom172 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_bottom183 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_bottom189 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_bottom192 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_bottom196 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_bottom200 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_bottom202 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_bottom205 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_bottom207 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_charm133 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_charm167 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_charm172 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_charm183 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_charm189 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_charm192 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_charm196 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_charm200 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_charm202 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_charm205 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_charm207 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_mu130 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_mu136 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_mu161 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_mu172 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_mu183 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_mu189 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_mu192 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_mu196 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_mu200 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_mu202 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_mu205 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_mu207 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_tau130 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_tau136 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_tau161 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_tau172 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_tau183 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_tau189 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_tau192 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_tau196 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_tau200 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_tau202 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_tau205 (double x) const
 
double getIntegrand_AFBnumeratorWithISR_tau207 (double x) const
 
double getIntegrand_dsigmaBox_bottom130 (double x) const
 
double getIntegrand_dsigmaBox_bottom133 (double x) const
 
double getIntegrand_dsigmaBox_bottom136 (double x) const
 
double getIntegrand_dsigmaBox_bottom161 (double x) const
 
double getIntegrand_dsigmaBox_bottom167 (double x) const
 
double getIntegrand_dsigmaBox_bottom172 (double x) const
 
double getIntegrand_dsigmaBox_bottom183 (double x) const
 
double getIntegrand_dsigmaBox_bottom189 (double x) const
 
double getIntegrand_dsigmaBox_bottom192 (double x) const
 
double getIntegrand_dsigmaBox_bottom196 (double x) const
 
double getIntegrand_dsigmaBox_bottom200 (double x) const
 
double getIntegrand_dsigmaBox_bottom202 (double x) const
 
double getIntegrand_dsigmaBox_bottom205 (double x) const
 
double getIntegrand_dsigmaBox_bottom207 (double x) const
 
double getIntegrand_dsigmaBox_charm130 (double x) const
 
double getIntegrand_dsigmaBox_charm133 (double x) const
 
double getIntegrand_dsigmaBox_charm136 (double x) const
 
double getIntegrand_dsigmaBox_charm161 (double x) const
 
double getIntegrand_dsigmaBox_charm167 (double x) const
 
double getIntegrand_dsigmaBox_charm172 (double x) const
 
double getIntegrand_dsigmaBox_charm183 (double x) const
 
double getIntegrand_dsigmaBox_charm189 (double x) const
 
double getIntegrand_dsigmaBox_charm192 (double x) const
 
double getIntegrand_dsigmaBox_charm196 (double x) const
 
double getIntegrand_dsigmaBox_charm200 (double x) const
 
double getIntegrand_dsigmaBox_charm202 (double x) const
 
double getIntegrand_dsigmaBox_charm205 (double x) const
 
double getIntegrand_dsigmaBox_charm207 (double x) const
 
double getIntegrand_dsigmaBox_down130 (double x) const
 
double getIntegrand_dsigmaBox_down133 (double x) const
 
double getIntegrand_dsigmaBox_down136 (double x) const
 
double getIntegrand_dsigmaBox_down161 (double x) const
 
double getIntegrand_dsigmaBox_down167 (double x) const
 
double getIntegrand_dsigmaBox_down172 (double x) const
 
double getIntegrand_dsigmaBox_down183 (double x) const
 
double getIntegrand_dsigmaBox_down189 (double x) const
 
double getIntegrand_dsigmaBox_down192 (double x) const
 
double getIntegrand_dsigmaBox_down196 (double x) const
 
double getIntegrand_dsigmaBox_down200 (double x) const
 
double getIntegrand_dsigmaBox_down202 (double x) const
 
double getIntegrand_dsigmaBox_down205 (double x) const
 
double getIntegrand_dsigmaBox_down207 (double x) const
 
double getIntegrand_dsigmaBox_mu130 (double x) const
 
double getIntegrand_dsigmaBox_mu133 (double x) const
 
double getIntegrand_dsigmaBox_mu136 (double x) const
 
double getIntegrand_dsigmaBox_mu161 (double x) const
 
double getIntegrand_dsigmaBox_mu167 (double x) const
 
double getIntegrand_dsigmaBox_mu172 (double x) const
 
double getIntegrand_dsigmaBox_mu183 (double x) const
 
double getIntegrand_dsigmaBox_mu189 (double x) const
 
double getIntegrand_dsigmaBox_mu192 (double x) const
 
double getIntegrand_dsigmaBox_mu196 (double x) const
 
double getIntegrand_dsigmaBox_mu200 (double x) const
 
double getIntegrand_dsigmaBox_mu202 (double x) const
 
double getIntegrand_dsigmaBox_mu205 (double x) const
 
double getIntegrand_dsigmaBox_mu207 (double x) const
 
double getIntegrand_dsigmaBox_strange130 (double x) const
 
double getIntegrand_dsigmaBox_strange133 (double x) const
 
double getIntegrand_dsigmaBox_strange136 (double x) const
 
double getIntegrand_dsigmaBox_strange161 (double x) const
 
double getIntegrand_dsigmaBox_strange167 (double x) const
 
double getIntegrand_dsigmaBox_strange172 (double x) const
 
double getIntegrand_dsigmaBox_strange183 (double x) const
 
double getIntegrand_dsigmaBox_strange189 (double x) const
 
double getIntegrand_dsigmaBox_strange192 (double x) const
 
double getIntegrand_dsigmaBox_strange196 (double x) const
 
double getIntegrand_dsigmaBox_strange200 (double x) const
 
double getIntegrand_dsigmaBox_strange202 (double x) const
 
double getIntegrand_dsigmaBox_strange205 (double x) const
 
double getIntegrand_dsigmaBox_strange207 (double x) const
 
double getIntegrand_dsigmaBox_tau130 (double x) const
 
double getIntegrand_dsigmaBox_tau133 (double x) const
 
double getIntegrand_dsigmaBox_tau136 (double x) const
 
double getIntegrand_dsigmaBox_tau161 (double x) const
 
double getIntegrand_dsigmaBox_tau167 (double x) const
 
double getIntegrand_dsigmaBox_tau172 (double x) const
 
double getIntegrand_dsigmaBox_tau183 (double x) const
 
double getIntegrand_dsigmaBox_tau189 (double x) const
 
double getIntegrand_dsigmaBox_tau192 (double x) const
 
double getIntegrand_dsigmaBox_tau196 (double x) const
 
double getIntegrand_dsigmaBox_tau200 (double x) const
 
double getIntegrand_dsigmaBox_tau202 (double x) const
 
double getIntegrand_dsigmaBox_tau205 (double x) const
 
double getIntegrand_dsigmaBox_tau207 (double x) const
 
double getIntegrand_dsigmaBox_up130 (double x) const
 
double getIntegrand_dsigmaBox_up133 (double x) const
 
double getIntegrand_dsigmaBox_up136 (double x) const
 
double getIntegrand_dsigmaBox_up161 (double x) const
 
double getIntegrand_dsigmaBox_up167 (double x) const
 
double getIntegrand_dsigmaBox_up172 (double x) const
 
double getIntegrand_dsigmaBox_up183 (double x) const
 
double getIntegrand_dsigmaBox_up189 (double x) const
 
double getIntegrand_dsigmaBox_up192 (double x) const
 
double getIntegrand_dsigmaBox_up196 (double x) const
 
double getIntegrand_dsigmaBox_up200 (double x) const
 
double getIntegrand_dsigmaBox_up202 (double x) const
 
double getIntegrand_dsigmaBox_up205 (double x) const
 
double getIntegrand_dsigmaBox_up207 (double x) const
 
double getIntegrand_sigmaWithISR_bottom130 (double x) const
 
double getIntegrand_sigmaWithISR_bottom133 (double x) const
 
double getIntegrand_sigmaWithISR_bottom136 (double x) const
 
double getIntegrand_sigmaWithISR_bottom161 (double x) const
 
double getIntegrand_sigmaWithISR_bottom167 (double x) const
 
double getIntegrand_sigmaWithISR_bottom172 (double x) const
 
double getIntegrand_sigmaWithISR_bottom183 (double x) const
 
double getIntegrand_sigmaWithISR_bottom189 (double x) const
 
double getIntegrand_sigmaWithISR_bottom192 (double x) const
 
double getIntegrand_sigmaWithISR_bottom196 (double x) const
 
double getIntegrand_sigmaWithISR_bottom200 (double x) const
 
double getIntegrand_sigmaWithISR_bottom202 (double x) const
 
double getIntegrand_sigmaWithISR_bottom205 (double x) const
 
double getIntegrand_sigmaWithISR_bottom207 (double x) const
 
double getIntegrand_sigmaWithISR_charm130 (double x) const
 
double getIntegrand_sigmaWithISR_charm133 (double x) const
 
double getIntegrand_sigmaWithISR_charm136 (double x) const
 
double getIntegrand_sigmaWithISR_charm161 (double x) const
 
double getIntegrand_sigmaWithISR_charm167 (double x) const
 
double getIntegrand_sigmaWithISR_charm172 (double x) const
 
double getIntegrand_sigmaWithISR_charm183 (double x) const
 
double getIntegrand_sigmaWithISR_charm189 (double x) const
 
double getIntegrand_sigmaWithISR_charm192 (double x) const
 
double getIntegrand_sigmaWithISR_charm196 (double x) const
 
double getIntegrand_sigmaWithISR_charm200 (double x) const
 
double getIntegrand_sigmaWithISR_charm202 (double x) const
 
double getIntegrand_sigmaWithISR_charm205 (double x) const
 
double getIntegrand_sigmaWithISR_charm207 (double x) const
 
double getIntegrand_sigmaWithISR_down130 (double x) const
 
double getIntegrand_sigmaWithISR_down133 (double x) const
 
double getIntegrand_sigmaWithISR_down136 (double x) const
 
double getIntegrand_sigmaWithISR_down161 (double x) const
 
double getIntegrand_sigmaWithISR_down167 (double x) const
 
double getIntegrand_sigmaWithISR_down172 (double x) const
 
double getIntegrand_sigmaWithISR_down183 (double x) const
 
double getIntegrand_sigmaWithISR_down189 (double x) const
 
double getIntegrand_sigmaWithISR_down192 (double x) const
 
double getIntegrand_sigmaWithISR_down196 (double x) const
 
double getIntegrand_sigmaWithISR_down200 (double x) const
 
double getIntegrand_sigmaWithISR_down202 (double x) const
 
double getIntegrand_sigmaWithISR_down205 (double x) const
 
double getIntegrand_sigmaWithISR_down207 (double x) const
 
double getIntegrand_sigmaWithISR_mu130 (double x) const
 
double getIntegrand_sigmaWithISR_mu136 (double x) const
 
double getIntegrand_sigmaWithISR_mu161 (double x) const
 
double getIntegrand_sigmaWithISR_mu172 (double x) const
 
double getIntegrand_sigmaWithISR_mu183 (double x) const
 
double getIntegrand_sigmaWithISR_mu189 (double x) const
 
double getIntegrand_sigmaWithISR_mu192 (double x) const
 
double getIntegrand_sigmaWithISR_mu196 (double x) const
 
double getIntegrand_sigmaWithISR_mu200 (double x) const
 
double getIntegrand_sigmaWithISR_mu202 (double x) const
 
double getIntegrand_sigmaWithISR_mu205 (double x) const
 
double getIntegrand_sigmaWithISR_mu207 (double x) const
 
double getIntegrand_sigmaWithISR_strange130 (double x) const
 
double getIntegrand_sigmaWithISR_strange133 (double x) const
 
double getIntegrand_sigmaWithISR_strange136 (double x) const
 
double getIntegrand_sigmaWithISR_strange161 (double x) const
 
double getIntegrand_sigmaWithISR_strange167 (double x) const
 
double getIntegrand_sigmaWithISR_strange172 (double x) const
 
double getIntegrand_sigmaWithISR_strange183 (double x) const
 
double getIntegrand_sigmaWithISR_strange189 (double x) const
 
double getIntegrand_sigmaWithISR_strange192 (double x) const
 
double getIntegrand_sigmaWithISR_strange196 (double x) const
 
double getIntegrand_sigmaWithISR_strange200 (double x) const
 
double getIntegrand_sigmaWithISR_strange202 (double x) const
 
double getIntegrand_sigmaWithISR_strange205 (double x) const
 
double getIntegrand_sigmaWithISR_strange207 (double x) const
 
double getIntegrand_sigmaWithISR_tau130 (double x) const
 
double getIntegrand_sigmaWithISR_tau136 (double x) const
 
double getIntegrand_sigmaWithISR_tau161 (double x) const
 
double getIntegrand_sigmaWithISR_tau172 (double x) const
 
double getIntegrand_sigmaWithISR_tau183 (double x) const
 
double getIntegrand_sigmaWithISR_tau189 (double x) const
 
double getIntegrand_sigmaWithISR_tau192 (double x) const
 
double getIntegrand_sigmaWithISR_tau196 (double x) const
 
double getIntegrand_sigmaWithISR_tau200 (double x) const
 
double getIntegrand_sigmaWithISR_tau202 (double x) const
 
double getIntegrand_sigmaWithISR_tau205 (double x) const
 
double getIntegrand_sigmaWithISR_tau207 (double x) const
 
double getIntegrand_sigmaWithISR_up130 (double x) const
 
double getIntegrand_sigmaWithISR_up133 (double x) const
 
double getIntegrand_sigmaWithISR_up136 (double x) const
 
double getIntegrand_sigmaWithISR_up161 (double x) const
 
double getIntegrand_sigmaWithISR_up167 (double x) const
 
double getIntegrand_sigmaWithISR_up172 (double x) const
 
double getIntegrand_sigmaWithISR_up183 (double x) const
 
double getIntegrand_sigmaWithISR_up189 (double x) const
 
double getIntegrand_sigmaWithISR_up192 (double x) const
 
double getIntegrand_sigmaWithISR_up196 (double x) const
 
double getIntegrand_sigmaWithISR_up200 (double x) const
 
double getIntegrand_sigmaWithISR_up202 (double x) const
 
double getIntegrand_sigmaWithISR_up205 (double x) const
 
double getIntegrand_sigmaWithISR_up207 (double x) const
 
double Integrand_AFBnumeratorWithISR_l (double x, const QCD::lepton l_flavor, const double s) const
 
double Integrand_AFBnumeratorWithISR_q (double x, const QCD::quark q_flavor, const double s) const
 
double Integrand_dsigmaBox_l (double cosTheta, const QCD::lepton l_flavor, const double s) const
 
double Integrand_dsigmaBox_q (double cosTheta, const QCD::quark q_flavor, const double s) const
 
double Integrand_sigmaWithISR_l (double x, const QCD::lepton l_flavor, const double s) const
 
double Integrand_sigmaWithISR_q (double x, const QCD::quark q_flavor, const double s) const
 
double m_q (const QCD::quark q, const double mu, const orders order=FULLNLO) const
 
double RAq (const QCD::quark q) const
 The radiator factor associated with the final-state QED and QCD corrections to the the axial-vector-current interactions, \(R_A^q(M_Z^2)\). More...
 
double resumKappaZ (const double DeltaRho[orders_EW_size], const double deltaKappa_rem[orders_EW_size], const double DeltaRbar_rem, const bool bool_Zbb) const
 A method to compute the real part of the effetvive coupling \(\kappa_Z^f\) from \(\Delta\rho\), \(\delta\rho_{\rm rem}^{f}\) and \(\Delta r_{\mathrm{rem}}\). More...
 
double resumMw (const double Mw_i, const double DeltaRho[orders_EW_size], const double DeltaR_rem[orders_EW_size]) const
 A method to compute the \(W\)-boson mass from \(\Delta\rho\) and \(\Delta r_{\mathrm{rem}}\). More...
 
double resumRhoZ (const double DeltaRho[orders_EW_size], const double deltaRho_rem[orders_EW_size], const double DeltaRbar_rem, const bool bool_Zbb) const
 A method to compute the real part of the effective coupling \(\rho_Z^f\) from \(\Delta\rho\), \(\delta\rho_{\rm rem}^{f}\) and \(\Delta r_{\mathrm{rem}}\). More...
 
double RVh () const
 The singlet vector corrections to the hadronic \(Z\)-boson width, denoted as \(R_V^h\). More...
 
double RVq (const QCD::quark q) const
 The radiator factor associated with the final-state QED and QCD corrections to the the vector-current interactions, \(R_V^q(M_Z^2)\). More...
 
double SchemeToDouble (const std::string scheme) const
 A method to convert a given scheme name in string form into a floating-point number with double precision. More...
 
virtual void setParameter (const std::string name, const double &value)
 A method to set the value of a parameter of StandardModel. More...
 
double sigma_NoISR_l (const QCD::lepton l_flavor, const double s) const
 
double sigma_NoISR_q (const QCD::quark q_flavor, const double s) const
 
double taub () const
 Top-mass corrections to the \(Zb\bar{b}\) vertex, denoted by \(\tau_b\). More...
 
- Protected Member Functions inherited from QCD
double MassOfNf (int nf) const
 The Mbar mass of the heaviest quark in the theory with Nf active flavour. More...
 
- Protected Attributes inherited from StandardModel
double A
 The CKM parameter \(A\) in the Wolfenstein parameterization. More...
 
double ale
 The fine-structure constant \(\alpha\). More...
 
double alpha21
 
double alpha31
 
double AlsMz
 The strong coupling constant at the Z-boson mass, \(\alpha_s(M_Z)\). More...
 
bool bSigmaForAFB
 
bool bSigmaForR
 
double dAle5Mz
 The five-flavour hadronic contribution to the electromagnetic coupling, \(\Delta\alpha_{\mathrm{had}}^{(5)}(M_Z^2)\). More...
 
double delGammaZ
 The theoretical uncertainty in \(\Gamma_Z\), denoted as \(\delta\,\Gamma_Z\), in GeV. More...
 
double delMw
 The theoretical uncertainty in \(M_W\), denoted as \(\delta\,M_W\), in GeV. More...
 
double delR0b
 The theoretical uncertainty in \(R_b^0\), denoted as \(\delta\,R_b^0\). More...
 
double delR0c
 The theoretical uncertainty in \(R_c^0\), denoted as \(\delta\,R_c^0\). More...
 
double delR0l
 The theoretical uncertainty in \(R_l^0\), denoted as \(\delta\,R_l^0\). More...
 
double delsigma0H
 The theoretical uncertainty in \(\sigma_{Hadron}^0\), denoted as \(\delta\,\sigma_{Hadron}^0\) in nb. More...
 
double delSin2th_b
 The theoretical uncertainty in \(\sin^2\theta_{\rm eff}^{b}\), denoted as \(\delta\sin^2\theta_{\rm eff}^{b}\). More...
 
double delSin2th_l
 The theoretical uncertainty in \(\sin^2\theta_{\rm eff}^{\rm lept}\), denoted as \(\delta\sin^2\theta_{\rm eff}^{\rm lept}\). More...
 
double delSin2th_q
 The theoretical uncertainty in \(\sin^2\theta_{\rm eff}^{q\not = b,t}\), denoted as \(\delta\sin^2\theta_{\rm eff}^{q\not = b,t}\). More...
 
double delta
 
double etab
 The CKM parameter \(\bar{\eta}\) in the Wolfenstein parameterization. More...
 
bool flag_order [orders_EW_size]
 An array of internal flags controlling the inclusions of higher-order corrections. More...
 
bool flagLEP2 [NUMofLEP2RCs]
 
double gamma
 \(\gamma \) used as an input for FlagWolfenstein = FALSE More...
 
double GF
 The Fermi constant \(G_\mu\) in \({\rm GeV}^{-2}\). More...
 
double lambda
 The CKM parameter \(\lambda\) in the Wolfenstein parameterization. More...
 
Particle leptons [6]
 An array of Particle objects for the leptons. More...
 
double mHl
 The Higgs mass \(m_h\) in GeV. More...
 
double muw
 A matching scale \(\mu_W\) around the weak scale in GeV. More...
 
CKM myCKM
 An object of type CKM. More...
 
PMNS myPMNS
 
double Mz
 The mass of the \(Z\) boson in GeV. More...
 
bool requireCKM
 An internal flag to control whether the CKM matrix has to be recomputed. More...
 
bool requireYe
 An internal flag to control whether the charged-lepton Yukawa matrix has to be recomputed. More...
 
bool requireYn
 An internal flag to control whether the neutrino Yukawa matrix has to be recomputed. More...
 
double rhob
 The CKM parameter \(\bar{\rho}\) in the Wolfenstein parameterization. More...
 
double s12
 
double s13
 
double s23
 
Flavour SMFlavour
 An object of type Flavour. More...
 
Matching< StandardModelMatching, StandardModelSMM
 An object of type Matching. More...
 
double Vcb
 \(\vert V_{cb} \vert \) used as an input for FlagWolfenstein = FALSE More...
 
double Vub
 \(\vert V_{ub} \vert \) used as an input for FlagWolfenstein = FALSE More...
 
double Vus
 \(\vert V_{us} \vert \) used as an input for FlagWolfenstein = FALSE More...
 
gslpp::matrix< gslpp::complexYd
 The Yukawa matrix of the down-type quarks. More...
 
gslpp::matrix< gslpp::complexYe
 The Yukawa matrix of the charged leptons. More...
 
gslpp::matrix< gslpp::complexYn
 The Yukawa matrix of the neutrinos. More...
 
gslpp::matrix< gslpp::complexYu
 The Yukawa matrix of the up-type quarks. More...
 
- Protected Attributes inherited from QCD
double AlsM
 The strong coupling constant at the mass scale MAls, \(\alpha_s(M_{\alpha_s})\). More...
 
double CA
 
double CF
 
bool computemt
 Switch for computing the \(\overline{\mathrm{MS}}\) mass of the top quark. More...
 
double dAdA_NA
 
double dFdA_NA
 
double dFdF_NA
 
double MAls
 The mass scale in GeV at which the strong coupling measurement is provided. More...
 
double mtpole
 The pole mass of the top quark. More...
 
double mub
 The threshold between five- and four-flavour theory in GeV. More...
 
double muc
 The threshold between four- and three-flavour theory in GeV. More...
 
double mut
 The threshold between six- and five-flavour theory in GeV. More...
 
double NA
 
double Nc
 The number of colours. More...
 
Particle quarks [6]
 The vector of all SM quarks. More...
 
bool requireYd
 Switch for generating the Yukawa couplings to the down-type quarks. More...
 
bool requireYu
 Switch for generating the Yukawa couplings to the up-type quarks. More...
 
double TF
 
- Protected Attributes inherited from Model
bool isSliced
 A boolean set to true if the current istance is a slice of an extended object. More...
 
std::map< std::string, std::reference_wrapper< const double > > ModelParamMap
 
bool UpdateError
 A boolean set to false if update is successful. More...
 

Constructor & Destructor Documentation

◆ EWSM_Output()

EWSM_Output::EWSM_Output ( const StandardModel SM_in)

Constructor.

Parameters
[in]EWSM_ina reference to an object of type EWSM

Definition at line 20 of file EWSM_Output.cpp.

21 : StandardModel(SM_in)
22 {
23 }

Member Function Documentation

◆ outputEachDeltaKappaZ()

void EWSM_Output::outputEachDeltaKappaZ ( const double  f_AlphaToGF,
const double  cW2overSW2,
const double  DeltaRho[StandardModel::orders_EW_size],
const double  deltaKappa_rem[StandardModel::orders_EW_size],
const double  DeltaRbar_rem,
const bool  bool_Zbb,
const double  taub[StandardModel::orders_EW_size],
const double  ZbbSubtract,
const double  Zgamma_EW2 
) const

Definition at line 583 of file EWSM_Output.cpp.

592 {
593  /* rescale */
594  double DeltaRho_new[StandardModel::orders_EW_size];
595  for (int j = 0; j < StandardModel::orders_EW_size; ++j)
596  DeltaRho_new[j] = cW2overSW2 * DeltaRho[j];
597 
598  if (getFlagKappaZ().compare("APPROXIMATEFORMULA") == 0) {
599  std::cout << "Delta kappaZb (from the approximate formula of sin2thb) = "
600  << kappaZ_f(getQuarks(QCD::BOTTOM)) - 1.0 << std::endl;
601  } else if (getFlagKappaZ().compare("NORESUM") == 0) {
602  std::cout << "Leading contributions: alpha or Gmu" << std::endl;
603  std::cout << " DeltaRho[EW1]=" << DeltaRho_new[StandardModel::EW1] << " "
604  << f_AlphaToGF * DeltaRho_new[StandardModel::EW1] << std::endl;
605  std::cout << " DeltaRho[EW1QCD1]=" << DeltaRho_new[StandardModel::EW1QCD1] << " "
606  << f_AlphaToGF * DeltaRho_new[StandardModel::EW1QCD1] << std::endl;
607  std::cout << " DeltaRho[EW1QCD2]=" << DeltaRho_new[StandardModel::EW1QCD2] << " "
608  << f_AlphaToGF * DeltaRho_new[StandardModel::EW1QCD2] << std::endl;
609  std::cout << " DeltaRho[EW2]=" << DeltaRho_new[StandardModel::EW2] << " "
610  << pow(f_AlphaToGF, 2.0) * DeltaRho_new[StandardModel::EW2] << std::endl;
611  std::cout << " DeltaRho[EW2QCD1]=" << DeltaRho_new[StandardModel::EW2QCD1] << " "
612  << pow(f_AlphaToGF, 2.0) * DeltaRho_new[StandardModel::EW2QCD1] << std::endl;
613  std::cout << " DeltaRho[EW3]=" << DeltaRho_new[StandardModel::EW3] << " "
614  << pow(f_AlphaToGF, 3.0) * DeltaRho_new[StandardModel::EW3] << std::endl;
615  std::cout << "EW2 from Z-gamma = " << Zgamma_EW2 << std::endl;
616  std::cout << "Remainder contributions: alpha or Gmu" << std::endl;
617  std::cout << " DeltaRbar_rem[EW1]=" << DeltaRbar_rem << " "
618  << f_AlphaToGF * DeltaRbar_rem << std::endl;
619  std::cout << " deltaKappa_rem[EW1]=" << deltaKappa_rem[StandardModel::EW1] - ZbbSubtract << " "
620  << f_AlphaToGF * (deltaKappa_rem[StandardModel::EW1] - ZbbSubtract) << std::endl;
621  std::cout << " deltaKappa_rem[EW1QCD1]=" << deltaKappa_rem[StandardModel::EW1QCD1] << " "
622  << f_AlphaToGF * deltaKappa_rem[StandardModel::EW1QCD1] << std::endl;
623  std::cout << " deltaKappa_rem[EW1QCD2]=" << deltaKappa_rem[StandardModel::EW1QCD2] << " "
624  << f_AlphaToGF * deltaKappa_rem[StandardModel::EW1QCD2] << std::endl;
625  std::cout << " deltaKappa_rem[EW2]=" << deltaKappa_rem[StandardModel::EW2] << " "
626  << pow(f_AlphaToGF, 2.0) * deltaKappa_rem[StandardModel::EW2] << std::endl;
627  if (bool_Zbb) {
628  std::cout << "Taub: alpha or Gmu" << std::endl;
629  std::cout << " taub[EW1]=" << taub[StandardModel::EW1] << " "
630  << f_AlphaToGF * taub[StandardModel::EW1] << std::endl;
631  std::cout << " taub[EW1QCD1]=" << taub[StandardModel::EW1QCD1] << " "
632  << f_AlphaToGF * taub[StandardModel::EW1QCD1] << std::endl;
633  std::cout << " taub[EW2]=" << taub[StandardModel::EW2] << " "
634  << pow(f_AlphaToGF, 2.0) * taub[StandardModel::EW2] << std::endl;
635  }
636  std::cout << "Each order: alpha or Gmu" << std::endl;
637  double dKappa_EW1 = DeltaRho_new[StandardModel::EW1] + deltaKappa_rem[StandardModel::EW1] - ZbbSubtract;
638  double dKappa_EW1QCD1 = DeltaRho_new[StandardModel::EW1QCD1] + deltaKappa_rem[StandardModel::EW1QCD1];
639  double dKappa_EW2 = (deltaKappa_rem[StandardModel::EW1] - ZbbSubtract) * DeltaRho_new[StandardModel::EW1]
640  - DeltaRho_new[StandardModel::EW1] * DeltaRbar_rem
641  + DeltaRho_new[StandardModel::EW2]
642  + deltaKappa_rem[StandardModel::EW2];
643  double dKappa_EW1QCD2 = DeltaRho_new[StandardModel::EW1QCD2] + deltaKappa_rem[StandardModel::EW1QCD2];
644  double dKappa_EW2QCD1 = DeltaRho_new[StandardModel::EW2QCD1] + deltaKappa_rem[StandardModel::EW1QCD1] * DeltaRho_new[StandardModel::EW1];
645  double dKappa_EW3 = DeltaRho_new[StandardModel::EW3];
646  double dKappa_EW2QCD2 = deltaKappa_rem[StandardModel::EW1QCD2] * DeltaRho_new[StandardModel::EW1];
647  //
648  double dKappa_EW1_TMP = dKappa_EW1;
649  double dKappa_EW1QCD1_TMP = dKappa_EW1QCD1;
650  double dKappa_EW2_TMP = dKappa_EW2;
651  double dKappa_EW1QCD2_TMP = dKappa_EW1QCD2;
652  double dKappa_EW2QCD1_TMP = dKappa_EW2QCD1;
653  double dKappa_EW3_TMP = dKappa_EW3;
654  double dKappa_EW2QCD2_TMP = dKappa_EW2QCD2;
655  //
656  if (bool_Zbb) {
657  dKappa_EW1 = dKappa_EW1_TMP - taub[StandardModel::EW1];
658  dKappa_EW1QCD1 = dKappa_EW1QCD1_TMP - taub[StandardModel::EW1QCD1];
659  dKappa_EW2 = dKappa_EW2_TMP - taub[StandardModel::EW2] + taub[StandardModel::EW1] * taub[StandardModel::EW1]
660  - dKappa_EW1_TMP * taub[StandardModel::EW1];
661  dKappa_EW1QCD2 = dKappa_EW1QCD2_TMP;
662  dKappa_EW2QCD1 = dKappa_EW2QCD1_TMP - dKappa_EW1_TMP * taub[StandardModel::EW1QCD1]
663  - dKappa_EW1QCD1_TMP * taub[StandardModel::EW1] + 2.0 * taub[StandardModel::EW1] * taub[StandardModel::EW1QCD1];
664  dKappa_EW3 = dKappa_EW3_TMP - dKappa_EW1_TMP * taub[StandardModel::EW2]
665  - dKappa_EW2_TMP * taub[StandardModel::EW1] + 2.0 * taub[StandardModel::EW2] * taub[StandardModel::EW1];
666  dKappa_EW2QCD2 = dKappa_EW2QCD2_TMP - dKappa_EW1QCD2_TMP * taub[StandardModel::EW1]
667  - dKappa_EW1QCD1_TMP * taub[StandardModel::EW1QCD1]
669  }
670  std::cout << " EW1: " << dKappa_EW1 << " " << f_AlphaToGF * dKappa_EW1 << std::endl;
671  std::cout << " EW1QCD1: " << dKappa_EW1QCD1 << " " << f_AlphaToGF * dKappa_EW1QCD1 << std::endl;
672  std::cout << " EW2: " << dKappa_EW2 + Zgamma_EW2
673  << " " << pow(f_AlphaToGF, 2.0) * dKappa_EW2 + Zgamma_EW2 << std::endl;
674  std::cout << " EW1QCD2: " << dKappa_EW1QCD2 << " " << f_AlphaToGF * dKappa_EW1QCD2 << std::endl;
675  std::cout << " EW2QCD1: " << dKappa_EW2QCD1 << " " << pow(f_AlphaToGF, 2.0) * dKappa_EW2QCD1 << std::endl;
676  std::cout << " EW3: " << dKappa_EW3 << " " << pow(f_AlphaToGF, 3.0) * dKappa_EW3 << std::endl;
677  std::cout << " EW2QCD2: " << dKappa_EW2QCD2 << " " << pow(f_AlphaToGF, 2.0) * dKappa_EW2QCD2 << std::endl;
678  std::cout << "Total contribution: alpha or Gmu" << std::endl;
679  std::cout << " kappaZ="
680  << 1.0 + dKappa_EW1 + dKappa_EW1QCD1 + dKappa_EW2
681  + dKappa_EW1QCD2 + dKappa_EW2QCD1 + dKappa_EW3
682  + dKappa_EW2QCD2 + Zgamma_EW2
683  << " "
684  << 1.0 + f_AlphaToGF * dKappa_EW1 + f_AlphaToGF * dKappa_EW1QCD1
685  + pow(f_AlphaToGF, 2.0) * dKappa_EW2 + f_AlphaToGF * dKappa_EW1QCD2
686  + pow(f_AlphaToGF, 2.0) * dKappa_EW2QCD1
687  + pow(f_AlphaToGF, 3.0) * dKappa_EW3
688  + pow(f_AlphaToGF, 2.0) * dKappa_EW2QCD2
689  + Zgamma_EW2 << std::endl;
690  if (bool_Zbb) {
691  std::cout << " kappaZ(taub resummed)="
692  << (1.0 + dKappa_EW1_TMP + dKappa_EW1QCD1_TMP
693  + dKappa_EW2_TMP + dKappa_EW1QCD2_TMP
694  + dKappa_EW2QCD1_TMP + dKappa_EW3_TMP
695  + dKappa_EW2QCD2_TMP)
697  + Zgamma_EW2
698  << " "
699  << (1.0 + f_AlphaToGF * dKappa_EW1_TMP
700  + f_AlphaToGF * dKappa_EW1QCD1_TMP
701  + pow(f_AlphaToGF, 2.0) * dKappa_EW2_TMP
702  + f_AlphaToGF * dKappa_EW1QCD2_TMP
703  + pow(f_AlphaToGF, 2.0) * dKappa_EW2QCD1_TMP
704  + pow(f_AlphaToGF, 3.0) * dKappa_EW3_TMP
705  + pow(f_AlphaToGF, 2.0) * dKappa_EW2QCD2_TMP)
706  / (1.0 + f_AlphaToGF * taub[StandardModel::EW1]
707  + f_AlphaToGF * taub[StandardModel::EW1QCD1]
708  + pow(f_AlphaToGF, 2.0) * taub[StandardModel::EW2])
709  + Zgamma_EW2
710  << std::endl;
711  }
712  } else
713  std::cout << "EWSM_Output::outputEachDeltaKappaZ(): Not implemented for schemeKappaZ="
714  << getFlagKappaZ() << std::endl;
715 }

◆ outputEachDeltaKappaZ_l()

void EWSM_Output::outputEachDeltaKappaZ_l ( const QCD::lepton  l,
const double  Mw_i 
) const

Definition at line 443 of file EWSM_Output.cpp.

444 {
445  std::cout << "================================================" << std::endl;
446  std::cout << "kappaZ_l[(QCD::lepton)" << l << "]" << std::endl;
447  std::cout << "Mw(input) = " << Mw_i << std::endl;
448 
449  double cW2_TMP = Mw_i * Mw_i / getMz() / getMz();
450  double sW2_TMP = 1.0 - cW2_TMP;
451 
452  double DeltaRho[StandardModel::orders_EW_size];
453  DeltaRho[StandardModel::EW1] = getMyOneLoopEW()->DeltaRho(Mw_i);
454  DeltaRho[StandardModel::EW1QCD1] = getMyTwoLoopQCD()->DeltaRho(Mw_i);
456  DeltaRho[StandardModel::EW2] = getMyTwoLoopEW()->DeltaRho(Mw_i);
458  DeltaRho[StandardModel::EW3] = getMyThreeLoopEW()->DeltaRho(Mw_i);
459 
460  /* compute delta kappa_rem^f */
462  deltaKappa_rem_f[StandardModel::EW1] = getMyOneLoopEW()->deltaKappa_rem_f(getLeptons(l), Mw_i);
463 #ifdef WITHIMTWOLOOPQCD
464  deltaKappa_rem_f[StandardModel::EW1QCD1] = gslpp::complex(getMyTwoLoopQCD()->deltaKappa_rem_f(getLeptons(l), Mw_i).real(),
465  getMyTwoLoopQCD()->deltaKappa_rem_f(getLeptons(l), Mw_i).imag(), false);
466 #else
467  deltaKappa_rem_f[StandardModel::EW1QCD1] = gslpp::complex(getMyTwoLoopQCD()->deltaKappa_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
468 #endif
469  deltaKappa_rem_f[StandardModel::EW1QCD2] = gslpp::complex(getMyThreeLoopQCD()->deltaKappa_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
470  deltaKappa_rem_f[StandardModel::EW2] = gslpp::complex(getMyTwoLoopEW()->deltaKappa_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
471  deltaKappa_rem_f[StandardModel::EW2QCD1] = gslpp::complex(getMyThreeLoopEW2QCD()->deltaKappa_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
472  deltaKappa_rem_f[StandardModel::EW3] = gslpp::complex(getMyThreeLoopEW()->deltaKappa_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
473 
474  /* compute Delta rbar_rem */
475  double DeltaRbar_rem = getMyOneLoopEW()->DeltaRbar_rem(Mw_i);
476 
477  /* conversion factor */
478  double f_AlphaToGF = sqrt(2.0) * getGF() * pow(getMz(), 2.0)
479  * sW2_TMP * cW2_TMP / M_PI / getAle();
480 
481  /* Re[Kappa_Z^f] with or without resummation */
482  double deltaKappa_rem_f_real[StandardModel::orders_EW_size];
483  for (int j = 0; j < StandardModel::orders_EW_size; ++j)
484  deltaKappa_rem_f_real[j] = deltaKappa_rem_f[j].real();
485 
486  /* O(alpha^2) correction to Re[kappa_Z^f] from the Z-gamma mixing */
487  double ReKappaZf = resumKappaZ(DeltaRho, deltaKappa_rem_f_real,
488  DeltaRbar_rem, false);
489  double Zgamma_EW2 = 35.0 * alphaMz() * alphaMz() / 18.0 / sW2_TMP
490  * (1.0 - 8.0 / 3.0 * ReKappaZf * sW2_TMP);
491 
492  double dummy[StandardModel::orders_EW_size];
493  outputEachDeltaKappaZ(f_AlphaToGF, cW2_TMP / sW2_TMP,
494  DeltaRho, deltaKappa_rem_f_real,
495  DeltaRbar_rem, false, dummy, 0.0, Zgamma_EW2);
496 
497  /* Im[kappa_Z^f] without resummation */
498  double ImKappaZf = 0.0;
499  for (int j = 0; j < StandardModel::orders_EW_size; ++j)
500  ImKappaZf += deltaKappa_rem_f[j].imag();
501  std::cout << "ImKappaZf(with alpha)=" << ImKappaZf << std::endl;
502 
503  std::cout << "================================================" << std::endl;
504 }

◆ outputEachDeltaKappaZ_q()

void EWSM_Output::outputEachDeltaKappaZ_q ( const QCD::quark  q,
const double  Mw_i 
) const

Definition at line 506 of file EWSM_Output.cpp.

507 {
508  std::cout << "================================================" << std::endl;
509  std::cout << "kappaZ_q[(QCD::quark)" << q << "]" << std::endl;
510  std::cout << "Mw(input) = " << Mw_i << std::endl;
511 
512  double cW2_TMP = Mw_i * Mw_i / getMz() / getMz();
513  double sW2_TMP = 1.0 - cW2_TMP;
514 
515  double DeltaRho[StandardModel::orders_EW_size];
516  DeltaRho[StandardModel::EW1] = getMyOneLoopEW()->DeltaRho(Mw_i);
517  DeltaRho[StandardModel::EW1QCD1] = getMyTwoLoopQCD()->DeltaRho(Mw_i);
519  DeltaRho[StandardModel::EW2] = getMyTwoLoopEW()->DeltaRho(Mw_i);
521  DeltaRho[StandardModel::EW3] = getMyThreeLoopEW()->DeltaRho(Mw_i);
522 
523  /* compute delta kappa_rem^f */
525  deltaKappa_rem_f[StandardModel::EW1] = getMyOneLoopEW()->deltaKappa_rem_f(getQuarks(q), Mw_i);
526 #ifdef WITHIMTWOLOOPQCD
527  deltaKappa_rem_f[StandardModel::EW1QCD1] = gslpp::complex(getMyTwoLoopQCD()->deltaKappa_rem_f(getQuarks(q), Mw_i).real(),
528  getMyTwoLoopQCD()->deltaKappa_rem_f(getQuarks(q), Mw_i).imag(), false);
529 #else
530  deltaKappa_rem_f[StandardModel::EW1QCD1] = gslpp::complex(getMyTwoLoopQCD()->deltaKappa_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
531 #endif
532  deltaKappa_rem_f[StandardModel::EW1QCD2] = gslpp::complex(getMyThreeLoopQCD()->deltaKappa_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
533  deltaKappa_rem_f[StandardModel::EW2] = gslpp::complex(getMyTwoLoopEW()->deltaKappa_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
534  deltaKappa_rem_f[StandardModel::EW2QCD1] = gslpp::complex(getMyThreeLoopEW2QCD()->deltaKappa_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
535  deltaKappa_rem_f[StandardModel::EW3] = gslpp::complex(getMyThreeLoopEW()->deltaKappa_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
536 
537  /* compute Delta rbar_rem */
538  double DeltaRbar_rem = getMyOneLoopEW()->DeltaRbar_rem(Mw_i);
539 
540  /* conversion factor */
541  double f_AlphaToGF = sqrt(2.0) * getGF() * pow(getMz(), 2.0)
542  * sW2_TMP * cW2_TMP / M_PI / getAle();
543 
544  /* Zbb */
545  bool bool_Zbb = false;
546  if (q == QCD::BOTTOM) bool_Zbb = true;
547  double ZbbSubtract = 0.0;
548  if (bool_Zbb)
549  ZbbSubtract = getAle() / 8.0 / M_PI / sW2_TMP
550  * pow(getMtpole() / Mw_i, 2.0);
552  double Xt = getMyEWSMcache()->Xt_alpha(Mw_i);
553  if (bool_Zbb) {
554  taub[StandardModel::EW1] = -2.0 * Xt;
555  taub[StandardModel::EW1QCD1] = 2.0 / 3.0 * M_PI * Xt * getMyEWSMcache()->alsMt();
556  taub[StandardModel::EW2] = -2.0 * Xt * Xt * getMyTwoLoopEW()->tau_2();
557  }
558 
559  /* Re[Kappa_Z^f] with or without resummation */
560  double deltaKappa_rem_f_real[StandardModel::orders_EW_size];
561  for (int j = 0; j < StandardModel::orders_EW_size; ++j)
562  deltaKappa_rem_f_real[j] = deltaKappa_rem_f[j].real();
563 
564  /* O(alpha^2) correction to Re[kappa_Z^f] from the Z-gamma mixing */
565  double ReKappaZf = resumKappaZ(DeltaRho, deltaKappa_rem_f_real,
566  DeltaRbar_rem, bool_Zbb);
567  double Zgamma_EW2 = 35.0 * alphaMz() * alphaMz() / 18.0 / sW2_TMP
568  * (1.0 - 8.0 / 3.0 * ReKappaZf * sW2_TMP);
569 
570  outputEachDeltaKappaZ(f_AlphaToGF, cW2_TMP / sW2_TMP,
571  DeltaRho, deltaKappa_rem_f_real,
572  DeltaRbar_rem, bool_Zbb, taub, ZbbSubtract, Zgamma_EW2);
573 
574  /* Im[kappa_Z^f] without resummation */
575  double ImKappaZf = 0.0;
576  for (int j = 0; j < StandardModel::orders_EW_size; ++j)
577  ImKappaZf += deltaKappa_rem_f[j].imag();
578  std::cout << "ImKappaZf(with alpha)=" << ImKappaZf << std::endl;
579 
580  std::cout << "================================================" << std::endl;
581 }

◆ outputEachDeltaR()

void EWSM_Output::outputEachDeltaR ( const double  Mw_i) const

Definition at line 25 of file EWSM_Output.cpp.

26 {
27  std::cout << "Mw_SM = " << Mw() << std::endl;
28  std::cout << "DeltaR_SM() = " << DeltaR() << std::endl;
29  std::cout << "DeltaRbar_SM() = " << DeltaRbar() << std::endl;
30  std::cout << "Mw(input) = " << Mw_i << std::endl;
31 
32  double cW2_TMP = Mw_i * Mw_i / getMz() / getMz();
33  double sW2_TMP = 1.0 - cW2_TMP;
34 
35  double DeltaRho[StandardModel::orders_EW_size];
36  DeltaRho[StandardModel::EW1] = getMyOneLoopEW()->DeltaRho(Mw_i);
37  DeltaRho[StandardModel::EW1QCD1] = getMyTwoLoopQCD()->DeltaRho(Mw_i);
39  DeltaRho[StandardModel::EW2] = getMyTwoLoopEW()->DeltaRho(Mw_i);
41  DeltaRho[StandardModel::EW3] = getMyThreeLoopEW()->DeltaRho(Mw_i);
42 
43  double DeltaR_rem[StandardModel::orders_EW_size];
44  DeltaR_rem[StandardModel::EW1] = getMyOneLoopEW()->DeltaR_rem(Mw_i);
45  DeltaR_rem[StandardModel::EW1QCD1] = getMyTwoLoopQCD()->DeltaR_rem(Mw_i);
47  DeltaR_rem[StandardModel::EW2] = getMyTwoLoopEW()->DeltaR_rem(Mw_i);
49  DeltaR_rem[StandardModel::EW3] = getMyThreeLoopEW()->DeltaR_rem(Mw_i);
50 
51  double f_AlphaToGF = sqrt(2.0) * getGF() * pow(getMz(), 2.0) * sW2_TMP * cW2_TMP / M_PI / getAle();
52  //f_AlphaToGF = 1.0; /* for test */
53  double DeltaRho_sum = f_AlphaToGF * DeltaRho[StandardModel::EW1]
54  + f_AlphaToGF * DeltaRho[StandardModel::EW1QCD1]
55  + f_AlphaToGF * DeltaRho[StandardModel::EW1QCD2]
56  + pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2]
57  + pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2QCD1]
58  + pow(f_AlphaToGF, 3.0) * DeltaRho[StandardModel::EW3];
59  double DeltaRho_G = f_AlphaToGF * DeltaRho[StandardModel::EW1];
60 
61  if (getFlagMw().compare("NORESUM") == 0) {
62 
63  f_AlphaToGF = 1.0;
64  DeltaRho[StandardModel::EW1QCD2] *= f_AlphaToGF;
65  DeltaRho[StandardModel::EW2QCD1] *= pow(f_AlphaToGF, 2.0);
66  DeltaRho[StandardModel::EW3] *= pow(f_AlphaToGF, 3.0);
67 
68  // Full EW one-loop contribution (without the full DeltaAlphaL5q)
69  double DeltaR_EW1 = -cW2_TMP / sW2_TMP * DeltaRho[StandardModel::EW1] + DeltaR_rem[StandardModel::EW1];
70 
71  // Full EW two-loop contribution with reducible corrections
72  double DeltaR_EW2_rem = getMyApproximateFormulae()->DeltaR_TwoLoopEW_rem(Mw_i);
73 
74  // EW two-loop irreducible contributions with large-mt expansion
75  double DeltaR_EW2_old_red = DeltaAlphaL5q() * DeltaAlphaL5q()
76  - 2.0 * cW2_TMP / sW2_TMP * DeltaAlphaL5q() * DeltaRho[StandardModel::EW1]
77  + pow(cW2_TMP / sW2_TMP * DeltaRho[StandardModel::EW1], 2.0);
78  double DeltaR_EW2_old_irred = -cW2_TMP / sW2_TMP * DeltaRho[StandardModel::EW2] + DeltaR_rem[StandardModel::EW2];
79 
80  // Delta r, including the full EW two-loop contribution
81  double deltaR = DeltaAlphaL5q();
82  for (int j = 0; j < StandardModel::orders_EW_size; ++j) {
83  deltaR += -cW2_TMP / sW2_TMP * DeltaRho[(StandardModel::orders_EW)j];
84  deltaR += DeltaR_rem[(StandardModel::orders_EW)j];
85  }
86  deltaR -= -cW2_TMP / sW2_TMP * DeltaRho[StandardModel::EW2];
87  deltaR -= DeltaR_rem[StandardModel::EW2];
88  deltaR += DeltaAlphaL5q() * DeltaAlphaL5q() + 2.0 * DeltaAlphaL5q() * DeltaR_EW1 + DeltaR_EW2_rem;
89 
90  std::cout << "(1+dr) - 1 = " << deltaR << std::endl;
91  std::cout << " EW1 = " << DeltaAlphaL5q() + DeltaR_EW1 << std::endl;
92  std::cout << " DeltaAlphaL5q = " << DeltaAlphaL5q() << std::endl;
93  std::cout << " dR = " << DeltaR_EW1 << std::endl;
94  std::cout << " EW1QCD1 = " << -cW2_TMP / sW2_TMP * DeltaRho[StandardModel::EW1QCD1] + DeltaR_rem[StandardModel::EW1QCD1] << std::endl;
95  std::cout << " EW2(full) = " << DeltaR_EW2_rem + DeltaAlphaL5q() * DeltaAlphaL5q() + 2.0 * DeltaAlphaL5q() * DeltaR_EW1 << std::endl;
96  std::cout << " dAle*dAle = " << DeltaAlphaL5q() * DeltaAlphaL5q() << std::endl;
97  std::cout << " 2*dAle*dR = " << 2.0 * DeltaAlphaL5q() * DeltaR_EW1 << std::endl;
98  std::cout << " others = " << DeltaR_EW2_rem << std::endl;
99  std::cout << " EW1QCD2 = " << -cW2_TMP / sW2_TMP * DeltaRho[StandardModel::EW1QCD2] + DeltaR_rem[StandardModel::EW1QCD2] << std::endl;
100  std::cout << " EW2QCD1 = " << -cW2_TMP / sW2_TMP * DeltaRho[StandardModel::EW2QCD1] + DeltaR_rem[StandardModel::EW2QCD1] << std::endl;
101  std::cout << " EW3 = " << -cW2_TMP / sW2_TMP * DeltaRho[StandardModel::EW3] + DeltaR_rem[StandardModel::EW3] << std::endl;
102  std::cout << " EW2(old,irreducible) = " << DeltaR_EW2_old_irred << std::endl;
103  std::cout << " EW2(old,red+irred) = " << DeltaR_EW2_old_red + DeltaR_EW2_old_irred << std::endl;
104  std::cout << " EW2(old,red+irred-dAle*dAle-2*dAle*dR) = "
105  << DeltaR_EW2_old_red + DeltaR_EW2_old_irred
107  - 2.0 * DeltaAlphaL5q() * DeltaR_EW1 << std::endl;
108 
109  } else if (getFlagMw().compare("OMSI") == 0) {
110 
111  // R = 1/(1 - Delta r)
112  double R = 1.0 / (1.0 + cW2_TMP / sW2_TMP * DeltaRho_sum)
113  / (1.0 - DeltaAlphaL5q() - DeltaR_rem[StandardModel::EW1] - DeltaR_rem[StandardModel::EW1QCD1] - DeltaR_rem[StandardModel::EW2]);
114 
115  std::cout << "1/(1-dr) - 1 (exact) = " << R - 1.0 << std::endl;
116  std::cout << " --> dr = " << 1.0 - 1.0 / R << std::endl;
117 
118  // each contribution
119  double DeltaR_EW1 = DeltaAlphaL5q() - cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1] + DeltaR_rem[StandardModel::EW1];
120  double DeltaR_EW1QCD1 = -cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1QCD1] + DeltaR_rem[StandardModel::EW1QCD1];
121  double DeltaR_EW2 = -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2]
122  + DeltaR_rem[StandardModel::EW2]
123  + cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1]
124  *(DeltaAlphaL5q() + DeltaR_rem[StandardModel::EW1])
125  + DeltaR_EW1*DeltaR_EW1;
126  double DeltaR_EW1QCD2 = -cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1QCD2];
127  double DeltaR_EW2QCD1 = -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2QCD1]
128  + cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1QCD1]
129  *(DeltaAlphaL5q() + DeltaR_rem[StandardModel::EW1])
130  + 2.0 * DeltaR_EW1*DeltaR_EW1QCD1;
131  double DeltaR_EW3 = -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 3.0) * DeltaRho[StandardModel::EW3]
132  + cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2]
133  *(DeltaAlphaL5q() + DeltaR_rem[StandardModel::EW1])
134  + pow(DeltaR_EW1, 3.0)
135  + 2.0 * DeltaR_EW1 * (DeltaR_EW2 - DeltaR_EW1 * DeltaR_EW1);
136 
137  std::cout << " EW1 = " << DeltaR_EW1 << std::endl;
138  std::cout << " DeltaAlphaL5q = " << DeltaAlphaL5q() << std::endl;
139  std::cout << " -cW2/sW2*dRho1= " << -cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1] << std::endl;
140  std::cout << " DeltaR1_rem = " << DeltaR_rem[StandardModel::EW1] << std::endl;
141  std::cout << " EW1QCD1 = " << DeltaR_EW1QCD1 << std::endl;
142  std::cout << " EW2(full) = " << DeltaR_EW2 << std::endl;
143  std::cout << " EW1*EW1 = " << DeltaR_EW1 * DeltaR_EW1 << std::endl;
144  std::cout << " dAle*dAle = " << DeltaAlphaL5q() * DeltaAlphaL5q() << std::endl;
145  std::cout << " others = " << DeltaR_EW1 * DeltaR_EW1 - DeltaAlphaL5q() * DeltaAlphaL5q() << std::endl;
146  std::cout << " -cW2/sW2*dRho2= " << -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2] << std::endl;
147  std::cout << " DeltaR2_rem = " << DeltaR_rem[StandardModel::EW2] << std::endl;
148  std::cout << " others = " << cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1]*(DeltaAlphaL5q() + DeltaR_rem[StandardModel::EW1]) << std::endl;
149  std::cout << " EW1QCD2 = " << DeltaR_EW1QCD2 << std::endl;
150  std::cout << " EW2QCD1 = " << DeltaR_EW2QCD1 << std::endl;
151  std::cout << " EW3 = " << DeltaR_EW3 << std::endl;
152  std::cout << " -cW2/sW2*dRho3= " << -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 3.0) * DeltaRho[StandardModel::EW3] << std::endl;
153  std::cout << " EW1^3 = " << pow(DeltaR_EW1, 3.0) << std::endl;
154  std::cout << " 2*EW1*(EW2-EW1^2)=" << 2.0 * DeltaR_EW1 * (DeltaR_EW2 - DeltaR_EW1 * DeltaR_EW1) << std::endl;
155  std::cout << " others = " << cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2]*(DeltaAlphaL5q() + DeltaR_rem[StandardModel::EW1]) << std::endl;
156 
157  } else if (getFlagMw().compare("OMSII") == 0) {
158 
159  // R = 1/(1 - Delta r)
160  double R = 1.0 / ((1.0 + cW2_TMP / sW2_TMP * DeltaRho_sum)*(1.0 - DeltaAlphaL5q())
161  - (1.0 + cW2_TMP / sW2_TMP * DeltaRho_G) * DeltaR_rem[StandardModel::EW1]
162  - DeltaR_rem[StandardModel::EW1QCD1] - DeltaR_rem[StandardModel::EW2]);
163 
164  // each contribution
165  double DeltaR_EW1 = DeltaAlphaL5q() - cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1] + DeltaR_rem[StandardModel::EW1];
166  double DeltaR_EW1QCD1 = -cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1QCD1] + DeltaR_rem[StandardModel::EW1QCD1];
167  double DeltaR_EW2 = -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2]
168  + DeltaR_rem[StandardModel::EW2]
169  + cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1]
170  *(DeltaAlphaL5q() + DeltaR_rem[StandardModel::EW1])
171  + DeltaR_EW1*DeltaR_EW1;
172  double DeltaR_EW1QCD2 = -cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1QCD2];
173  double DeltaR_EW2QCD1 = -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2QCD1]
174  + cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1QCD1] * DeltaAlphaL5q()
175  + 2.0 * DeltaR_EW1*DeltaR_EW1QCD1;
176  double DeltaR_EW3 = -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 3.0) * DeltaRho[StandardModel::EW3]
177  + cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2] * DeltaAlphaL5q()
178  + pow(DeltaR_EW1, 3.0)
179  + 2.0 * DeltaR_EW1 * (DeltaR_EW2 - DeltaR_EW1 * DeltaR_EW1);
180 
181  std::cout << "1/(1-dr) - 1 (exact) = " << R - 1.0 << std::endl;
182  std::cout << " --> dr = " << 1.0 - 1.0 / R << std::endl;
183  std::cout << "1/(1-dr) - 1 (sum of expanded terms) = "
184  << DeltaR_EW1 + DeltaR_EW1QCD1 + DeltaR_EW2 + DeltaR_EW1QCD2
185  + DeltaR_EW2QCD1 + DeltaR_EW3 << std::endl;
186  std::cout << " EW1 = " << DeltaR_EW1 << std::endl;
187  std::cout << " DeltaAlphaL5q = " << DeltaAlphaL5q() << std::endl;
188  std::cout << " -cW2/sW2*dRho1= " << -cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1] << std::endl;
189  std::cout << " DeltaR1_rem = " << DeltaR_rem[StandardModel::EW1] << std::endl;
190  std::cout << " EW1QCD1 = " << DeltaR_EW1QCD1 << std::endl;
191  std::cout << " EW2(full) = " << DeltaR_EW2 << std::endl;
192  std::cout << " EW1*EW1 = " << DeltaR_EW1 * DeltaR_EW1 << std::endl;
193  std::cout << " dAle*dAle = " << DeltaAlphaL5q() * DeltaAlphaL5q() << std::endl;
194  std::cout << " others = " << DeltaR_EW1 * DeltaR_EW1 - DeltaAlphaL5q() * DeltaAlphaL5q() << std::endl;
195  std::cout << " -cW2/sW2*dRho2= " << -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2] << std::endl;
196  std::cout << " DeltaR2_rem = " << DeltaR_rem[StandardModel::EW2] << std::endl;
197  std::cout << " others = " << cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1]*(DeltaAlphaL5q() + DeltaR_rem[StandardModel::EW1]) << std::endl;
198  std::cout << " EW1QCD2 = " << DeltaR_EW1QCD2 << std::endl;
199  std::cout << " EW2QCD1 = " << DeltaR_EW2QCD1 << std::endl;
200  std::cout << " EW3 = " << DeltaR_EW3 << std::endl;
201  std::cout << " -cW2/sW2*dRho3= " << -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 3.0) * DeltaRho[StandardModel::EW3] << std::endl;
202  std::cout << " EW1^3 = " << pow(DeltaR_EW1, 3.0) << std::endl;
203  std::cout << " 2*EW1*(EW2-EW1^2)=" << 2.0 * DeltaR_EW1 * (DeltaR_EW2 - DeltaR_EW1 * DeltaR_EW1) << std::endl;
204  std::cout << " others = " << cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2] * DeltaAlphaL5q() << std::endl;
205 
206  } else
207  std::cout << "EWSM_Output::outputEachDeltaR(): Not implemented for schemeMw="
208  << getFlagMw() << std::endl;
209 }

◆ outputEachDeltaRhoZ()

void EWSM_Output::outputEachDeltaRhoZ ( const double  f_AlphaToGF,
const double  DeltaRho[StandardModel::orders_EW_size],
const double  deltaRho_rem[StandardModel::orders_EW_size],
const double  DeltaRbar_rem,
const bool  bool_Zbb,
const double  taub[StandardModel::orders_EW_size],
const double  ZbbSubtract 
) const

Definition at line 336 of file EWSM_Output.cpp.

343 {
344  if (getFlagRhoZ().compare("APPROXIMATEFORMULA") == 0) {
345 
346  } else if (getFlagRhoZ().compare("NORESUM") == 0) {
347  std::cout << "Leading contributions: alpha or Gmu" << std::endl;
348  std::cout << " DeltaRho[EW1]=" << DeltaRho[StandardModel::EW1] << " "
349  << f_AlphaToGF * DeltaRho[StandardModel::EW1] << std::endl;
350  std::cout << " DeltaRho[EW1QCD1]=" << DeltaRho[StandardModel::EW1QCD1] << " "
351  << f_AlphaToGF * DeltaRho[StandardModel::EW1QCD1] << std::endl;
352  std::cout << " DeltaRho[EW1QCD2]=" << DeltaRho[StandardModel::EW1QCD2] << " "
353  << f_AlphaToGF * DeltaRho[StandardModel::EW1QCD2] << std::endl;
354  std::cout << " DeltaRho[EW2]=" << DeltaRho[StandardModel::EW2] << " "
355  << pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2] << std::endl;
356  std::cout << " DeltaRho[EW2QCD1]=" << DeltaRho[StandardModel::EW2QCD1] << " "
357  << pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2QCD1] << std::endl;
358  std::cout << " DeltaRho[EW3]=" << DeltaRho[StandardModel::EW3] << " "
359  << pow(f_AlphaToGF, 3.0) * DeltaRho[StandardModel::EW3] << std::endl;
360  std::cout << "Remainder contributions: alpha or Gmu" << std::endl;
361  std::cout << " DeltaRbar_rem[EW1]=" << DeltaRbar_rem << " "
362  << f_AlphaToGF * DeltaRbar_rem << std::endl;
363  std::cout << " deltaRho_rem[EW1]=" << deltaRho_rem[StandardModel::EW1] - ZbbSubtract << " "
364  << f_AlphaToGF * (deltaRho_rem[StandardModel::EW1] - ZbbSubtract) << std::endl;
365  std::cout << " deltaRho_rem[EW1QCD1]=" << deltaRho_rem[StandardModel::EW1QCD1] << " "
366  << f_AlphaToGF * deltaRho_rem[StandardModel::EW1QCD1] << std::endl;
367  std::cout << " deltaRho_rem[EW2]=" << deltaRho_rem[StandardModel::EW2] << " "
368  << pow(f_AlphaToGF, 2.0) * deltaRho_rem[StandardModel::EW2] << std::endl;
369  if (bool_Zbb) {
370  std::cout << "Taub: alpha or Gmu" << std::endl;
371  std::cout << " taub[EW1]=" << taub[StandardModel::EW1] << " "
372  << f_AlphaToGF * taub[StandardModel::EW1] << std::endl;
373  std::cout << " taub[EW1QCD1]=" << taub[StandardModel::EW1QCD1] << " "
374  << f_AlphaToGF * taub[StandardModel::EW1QCD1] << std::endl;
375  std::cout << " taub[EW2]=" << taub[StandardModel::EW2] << " "
376  << pow(f_AlphaToGF, 2.0) * taub[StandardModel::EW2] << std::endl;
377  }
378  std::cout << "Each order: alpha or Gmu" << std::endl;
379  double dRho_EW1 = DeltaRho[StandardModel::EW1] + deltaRho_rem[StandardModel::EW1] - ZbbSubtract;
380  double dRho_EW1QCD1 = DeltaRho[StandardModel::EW1QCD1] + deltaRho_rem[StandardModel::EW1QCD1];
381  double dRho_EW2 = DeltaRho[StandardModel::EW1] * DeltaRho[StandardModel::EW1]
382  + (deltaRho_rem[StandardModel::EW1] - ZbbSubtract) * DeltaRho[StandardModel::EW1]
383  - DeltaRho[StandardModel::EW1] * DeltaRbar_rem
384  + DeltaRho[StandardModel::EW2]
385  + deltaRho_rem[StandardModel::EW2];
386  double dRho_EW1QCD2 = DeltaRho[StandardModel::EW1QCD2];
387  double dRho_EW2QCD1 = DeltaRho[StandardModel::EW2QCD1] + deltaRho_rem[StandardModel::EW1QCD1] * DeltaRho[StandardModel::EW1];
388  double dRho_EW3 = DeltaRho[StandardModel::EW3];
389  //
390  double dRho_EW1_TMP = dRho_EW1;
391  double dRho_EW1QCD1_TMP = dRho_EW1QCD1;
392  double dRho_EW2_TMP = dRho_EW2;
393  double dRho_EW1QCD2_TMP = dRho_EW1QCD2;
394  double dRho_EW2QCD1_TMP = dRho_EW2QCD1;
395  double dRho_EW3_TMP = dRho_EW3;
396  //
397  if (bool_Zbb) {
398  dRho_EW1 = dRho_EW1_TMP + 2.0 * taub[StandardModel::EW1];
399  dRho_EW1QCD1 = dRho_EW1QCD1_TMP + 2.0 * taub[StandardModel::EW1QCD1];
400  dRho_EW2 = dRho_EW2_TMP + 2.0 * taub[StandardModel::EW2] + taub[StandardModel::EW1] * taub[StandardModel::EW1]
401  + dRho_EW1_TMP * 2.0 * taub[StandardModel::EW1];
402  dRho_EW1QCD2 = dRho_EW1QCD2_TMP;
403  dRho_EW2QCD1 = dRho_EW2QCD1_TMP + dRho_EW1_TMP * 2.0 * taub[StandardModel::EW1QCD1]
404  + dRho_EW1QCD1_TMP * 2.0 * taub[StandardModel::EW1] + 2.0 * taub[StandardModel::EW1] * taub[StandardModel::EW1QCD1];
405  dRho_EW3 = dRho_EW3_TMP + dRho_EW1_TMP * 2.0 * taub[StandardModel::EW2]
406  + dRho_EW2_TMP * 2.0 * taub[StandardModel::EW1] + 2.0 * taub[StandardModel::EW2] * taub[StandardModel::EW1];
407  }
408  std::cout << " EW1: " << dRho_EW1 << " " << f_AlphaToGF * dRho_EW1 << std::endl;
409  std::cout << " EW1QCD1: " << dRho_EW1QCD1 << " " << f_AlphaToGF * dRho_EW1QCD1 << std::endl;
410  std::cout << " EW2: " << dRho_EW2 << " " << pow(f_AlphaToGF, 2.0) * dRho_EW2 << std::endl;
411  std::cout << " EW1QCD2: " << dRho_EW1QCD2 << " " << f_AlphaToGF * dRho_EW1QCD2 << std::endl;
412  std::cout << " EW2QCD1: " << dRho_EW2QCD1 << " " << pow(f_AlphaToGF, 2.0) * dRho_EW2QCD1 << std::endl;
413  std::cout << " EW3: " << dRho_EW3 << " " << pow(f_AlphaToGF, 3.0) * dRho_EW3 << std::endl;
414  std::cout << "Total contribution: alpha or Gmu" << std::endl;
415  std::cout << " rhoZ="
416  << 1.0 + dRho_EW1 + dRho_EW1QCD1 + dRho_EW2
417  + dRho_EW1QCD2 + dRho_EW2QCD1 + dRho_EW3
418  << " "
419  << 1.0 + f_AlphaToGF * dRho_EW1 + f_AlphaToGF * dRho_EW1QCD1
420  + pow(f_AlphaToGF, 2.0) * dRho_EW2
421  + f_AlphaToGF * dRho_EW1QCD2 + pow(f_AlphaToGF, 2.0) * dRho_EW2QCD1
422  + pow(f_AlphaToGF, 3.0) * dRho_EW3
423  << std::endl;
424  if (bool_Zbb) {
425  std::cout << " rhoZ(taub resummed)="
426  << (1.0 + dRho_EW1_TMP + dRho_EW1QCD1_TMP + dRho_EW2_TMP
427  + dRho_EW1QCD2_TMP + dRho_EW2QCD1_TMP + dRho_EW3_TMP)
429  << " "
430  << (1.0 + f_AlphaToGF * dRho_EW1_TMP + f_AlphaToGF * dRho_EW1QCD1_TMP
431  + pow(f_AlphaToGF, 2.0) * dRho_EW2_TMP
432  + f_AlphaToGF * dRho_EW1QCD2_TMP + pow(f_AlphaToGF, 2.0) * dRho_EW2QCD1_TMP
433  + pow(f_AlphaToGF, 3.0) * dRho_EW3_TMP)
434  * pow(1.0 + f_AlphaToGF * taub[StandardModel::EW1] + f_AlphaToGF * taub[StandardModel::EW1QCD1]
435  + pow(f_AlphaToGF, 2.0) * taub[StandardModel::EW2], 2.0)
436  << std::endl;
437  }
438  } else
439  std::cout << "EWSM_Output::outputEachDeltaRhoZ(): Not implemented for schemeRhoZ="
440  << getFlagRhoZ() << std::endl;
441 }

◆ outputEachDeltaRhoZ_l()

void EWSM_Output::outputEachDeltaRhoZ_l ( const QCD::lepton  l,
const double  Mw_i 
) const

Definition at line 211 of file EWSM_Output.cpp.

212 {
213  std::cout << "================================================" << std::endl;
214  std::cout << "rhoZ_l[(QCD::lepton)" << l << "]" << std::endl;
215  std::cout << "Mw(input) = " << Mw_i << std::endl;
216 
217  double cW2_TMP = Mw_i * Mw_i / getMz() / getMz();
218  double sW2_TMP = 1.0 - cW2_TMP;
219 
220  double DeltaRho[StandardModel::orders_EW_size];
221  DeltaRho[StandardModel::EW1] = getMyOneLoopEW()->DeltaRho(Mw_i);
222  DeltaRho[StandardModel::EW1QCD1] = getMyTwoLoopQCD()->DeltaRho(Mw_i);
224  DeltaRho[StandardModel::EW2] = getMyTwoLoopEW()->DeltaRho(Mw_i);
226  DeltaRho[StandardModel::EW3] = getMyThreeLoopEW()->DeltaRho(Mw_i);
227 
228  /* compute delta rho_rem^f */
230  deltaRho_rem_f[StandardModel::EW1] = getMyOneLoopEW()->deltaRho_rem_f(getLeptons(l), Mw_i);
231 #ifdef WITHIMTWOLOOPQCD
232  deltaRho_rem_f[StandardModel::EW1QCD1] = gslpp::complex(getMyTwoLoopQCD()->deltaRho_rem_f(getLeptons(l), Mw_i).real(),
233  getMyTwoLoopQCD()->deltaRho_rem_f(getLeptons(l), Mw_i).imag(), false);
234 #else
235  deltaRho_rem_f[StandardModel::EW1QCD1] = gslpp::complex(getMyTwoLoopQCD()->deltaRho_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
236 #endif
237  deltaRho_rem_f[StandardModel::EW1QCD2] = gslpp::complex(getMyThreeLoopQCD()->deltaRho_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
238  deltaRho_rem_f[StandardModel::EW2] = gslpp::complex(getMyTwoLoopEW()->deltaRho_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
239  deltaRho_rem_f[StandardModel::EW2QCD1] = gslpp::complex(getMyThreeLoopEW2QCD()->deltaRho_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
240  deltaRho_rem_f[StandardModel::EW3] = gslpp::complex(getMyThreeLoopEW()->deltaRho_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
241 
242  /* compute Delta rbar_rem */
243  double DeltaRbar_rem = getMyOneLoopEW()->DeltaRbar_rem(Mw_i);
244 
245  double f_AlphaToGF = sqrt(2.0) * getGF() * pow(getMz(), 2.0)
246  * sW2_TMP * cW2_TMP / M_PI / getAle();
247 
248  /* Re[rho_Z^f] with or without resummation */
249  double deltaRho_rem_f_real[StandardModel::orders_EW_size];
250  for (int j = 0; j < StandardModel::orders_EW_size; ++j)
251  deltaRho_rem_f_real[j] = deltaRho_rem_f[j].real();
252 
253  double dummy[StandardModel::orders_EW_size];
254  outputEachDeltaRhoZ(f_AlphaToGF, DeltaRho, deltaRho_rem_f_real,
255  DeltaRbar_rem, false, dummy, 0.0);
256 
257  /* Im[rho_Z^f] without resummation */
258  double ImRhoZf = 0.0;
259  for (int j = 0; j < StandardModel::orders_EW_size; ++j)
260  ImRhoZf += deltaRho_rem_f[j].imag();
261  std::cout << "ImRhoZf(with alpha)=" << ImRhoZf << std::endl;
262 
263  std::cout << "================================================" << std::endl;
264 }

◆ outputEachDeltaRhoZ_q()

void EWSM_Output::outputEachDeltaRhoZ_q ( const QCD::quark  q,
const double  Mw_i 
) const

Definition at line 266 of file EWSM_Output.cpp.

267 {
268  std::cout << "================================================" << std::endl;
269  std::cout << "rhoZ_q[(QCD::quark)" << q << "]" << std::endl;
270  std::cout << "Mw(input) = " << Mw_i << std::endl;
271 
272  double cW2_TMP = Mw_i * Mw_i / getMz() / getMz();
273  double sW2_TMP = 1.0 - cW2_TMP;
274 
275  double DeltaRho[StandardModel::orders_EW_size];
276  DeltaRho[StandardModel::EW1] = getMyOneLoopEW()->DeltaRho(Mw_i);
277  DeltaRho[StandardModel::EW1QCD1] = getMyTwoLoopQCD()->DeltaRho(Mw_i);
279  DeltaRho[StandardModel::EW2] = getMyTwoLoopEW()->DeltaRho(Mw_i);
281  DeltaRho[StandardModel::EW3] = getMyThreeLoopEW()->DeltaRho(Mw_i);
282 
283  /* compute delta rho_rem^f */
285  deltaRho_rem_f[StandardModel::EW1] = getMyOneLoopEW()->deltaRho_rem_f(getQuarks(q), Mw_i);
286 #ifdef WITHIMTWOLOOPQCD
287  deltaRho_rem_f[StandardModel::EW1QCD1] = gslpp::complex(getMyTwoLoopQCD()->deltaRho_rem_f(getQuarks(q), Mw_i).real(),
288  getMyTwoLoopQCD()->deltaRho_rem_f(getQuarks(q), Mw_i).imag(), false);
289 #else
290  deltaRho_rem_f[StandardModel::EW1QCD1] = gslpp::complex(getMyTwoLoopQCD()->deltaRho_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
291 #endif
292  deltaRho_rem_f[StandardModel::EW1QCD2] = gslpp::complex(getMyThreeLoopQCD()->deltaRho_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
293  deltaRho_rem_f[StandardModel::EW2] = gslpp::complex(getMyTwoLoopEW()->deltaRho_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
294  deltaRho_rem_f[StandardModel::EW2QCD1] = gslpp::complex(getMyThreeLoopEW2QCD()->deltaRho_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
295  deltaRho_rem_f[StandardModel::EW3] = gslpp::complex(getMyThreeLoopEW()->deltaRho_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
296 
297  /* compute Delta rbar_rem */
298  double DeltaRbar_rem = getMyOneLoopEW()->DeltaRbar_rem(Mw_i);
299 
300  /* conversion factor */
301  double f_AlphaToGF = sqrt(2.0) * getGF() * pow(getMz(), 2.0)
302  * sW2_TMP * cW2_TMP / M_PI / getAle();
303 
304  /* Zbb */
305  bool bool_Zbb = false;
306  if (q == QCD::BOTTOM) bool_Zbb = true;
307  double ZbbSubtract = 0.0;
308  if (bool_Zbb)
309  ZbbSubtract = -getAle() / 4.0 / M_PI / sW2_TMP
310  * pow(getMtpole() / Mw_i, 2.0);
312  double Xt = getMyEWSMcache()->Xt_alpha(Mw_i);
313  if (bool_Zbb) {
314  taub[StandardModel::EW1] = -2.0 * Xt;
315  taub[StandardModel::EW1QCD1] = 2.0 / 3.0 * M_PI * Xt * getMyEWSMcache()->alsMt();
316  taub[StandardModel::EW2] = -2.0 * Xt * Xt * getMyTwoLoopEW()->tau_2();
317  }
318 
319  /* Re[rho_Z^f] with or without resummation */
320  double deltaRho_rem_f_real[StandardModel::orders_EW_size];
321  for (int j = 0; j < StandardModel::orders_EW_size; ++j)
322  deltaRho_rem_f_real[j] = deltaRho_rem_f[j].real();
323 
324  outputEachDeltaRhoZ(f_AlphaToGF, DeltaRho, deltaRho_rem_f_real,
325  DeltaRbar_rem, bool_Zbb, taub, ZbbSubtract);
326 
327  /* Im[rho_Z^f] without resummation */
328  double ImRhoZf = 0.0;
329  for (int j = 0; j < StandardModel::orders_EW_size; ++j)
330  ImRhoZf += deltaRho_rem_f[j].imag();
331  std::cout << "ImRhoZf(with alpha)=" << ImRhoZf << std::endl;
332 
333  std::cout << "================================================" << std::endl;
334 }

The documentation for this class was generated from the following files:
EWSMThreeLoopEW::DeltaR_rem
double DeltaR_rem(const double Mw_i) const
Remainder contribution of to , denoted as .
Definition: EWSMThreeLoopEW.cpp:109
StandardModel::getMyOneLoopEW
EWSMOneLoopEW * getMyOneLoopEW() const
A get method to retrieve the member pointer of type EWSMOneLoopEW,.
Definition: StandardModel.h:970
StandardModel::EW2
Two-loop of .
Definition: StandardModel.h:499
StandardModel::taub
double taub() const
Top-mass corrections to the vertex, denoted by .
Definition: StandardModel.cpp:2104
QCD::BOTTOM
Definition: QCD.h:329
StandardModel::DeltaRbar
virtual double DeltaRbar() const
The SM prediction for derived from that for the -boson mass.
Definition: StandardModel.cpp:1146
EWSMThreeLoopEW::DeltaRho
double DeltaRho(const double Mw_i) const
Leading three-loop contribution of to , denoted as .
Definition: EWSMThreeLoopEW.cpp:70
EWSMOneLoopEW::DeltaRbar_rem
double DeltaRbar_rem(const double Mw_i) const
.
Definition: EWSMOneLoopEW.cpp:73
EWSM_Output::outputEachDeltaRhoZ
void outputEachDeltaRhoZ(const double f_AlphaToGF, const double DeltaRho[StandardModel::orders_EW_size], const double deltaRho_rem[StandardModel::orders_EW_size], const double DeltaRbar_rem, const bool bool_Zbb, const double taub[StandardModel::orders_EW_size], const double ZbbSubtract) const
Definition: EWSM_Output.cpp:336
EWSMOneLoopEW::DeltaRho
double DeltaRho(const double Mw_i) const
Leading one-loop contribution of to , denoted as .
Definition: EWSMOneLoopEW.cpp:43
StandardModel::getFlagRhoZ
std::string getFlagRhoZ() const
A method to retrieve the model flag RhoZ.
Definition: StandardModel.h:672
StandardModel::EW1
One-loop of .
Definition: StandardModel.h:496
StandardModel::getMyEWSMcache
EWSMcache * getMyEWSMcache() const
A get method to retrieve the member pointer of type EWSMcache.
Definition: StandardModel.h:961
StandardModel::getMyTwoLoopQCD
EWSMTwoLoopQCD * getMyTwoLoopQCD() const
Definition: StandardModel.h:1015
StandardModel::getMyThreeLoopQCD
EWSMThreeLoopQCD * getMyThreeLoopQCD() const
Definition: StandardModel.h:1005
StandardModel::alphaMz
double alphaMz() const
The electromagnetic coupling at the -mass scale, .
Definition: StandardModel.cpp:893
StandardModel::DeltaAlphaL5q
double DeltaAlphaL5q() const
The sum of the leptonic and the five-flavour hadronic corrections to the electromagnetic coupling at...
Definition: StandardModel.cpp:856
gslpp::complex
A class for defining operations on and functions of complex numbers.
Definition: gslpp_complex.h:35
EWSMOneLoopEW::deltaKappa_rem_f
gslpp::complex deltaKappa_rem_f(const Particle f, const double Mw_i) const
Remainder contribution of to the effective couplings , denoted as .
Definition: EWSMOneLoopEW.cpp:143
StandardModel::getMyThreeLoopEW
EWSMThreeLoopEW * getMyThreeLoopEW() const
Definition: StandardModel.h:995
EWSMOneLoopEW::deltaRho_rem_f
gslpp::complex deltaRho_rem_f(const Particle f, const double Mw_i) const
Remainder contribution of to the effective couplings , denoted as .
Definition: EWSMOneLoopEW.cpp:113
EWSMThreeLoopEW2QCD::DeltaRho
double DeltaRho(const double Mw_i) const
Leading three-loop contribution of to , denoted as .
Definition: EWSMThreeLoopEW2QCD.cpp:29
StandardModel::DeltaR
virtual double DeltaR() const
The SM prediction for derived from that for the boson mass.
Definition: StandardModel.cpp:1036
EWSMTwoLoopQCD::DeltaR_rem
double DeltaR_rem(const double Mw_i) const
Remainder contribution of to , denoted as .
Definition: EWSMTwoLoopQCD.cpp:44
StandardModel::StandardModel
StandardModel()
The default constructor.
Definition: StandardModel.cpp:40
StandardModel::EW1QCD1
Two-loop of .
Definition: StandardModel.h:497
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:739
gslpp::sqrt
complex sqrt(const complex &z)
Definition: gslpp_complex.cpp:385
EWSMThreeLoopQCD::DeltaRho
double DeltaRho(const double Mw_i) const
Leading three-loop QCD contribution of to , denoted as .
Definition: EWSMThreeLoopQCD.cpp:42
EWSM_Output::outputEachDeltaKappaZ
void outputEachDeltaKappaZ(const double f_AlphaToGF, const double cW2overSW2, const double DeltaRho[StandardModel::orders_EW_size], const double deltaKappa_rem[StandardModel::orders_EW_size], const double DeltaRbar_rem, const bool bool_Zbb, const double taub[StandardModel::orders_EW_size], const double ZbbSubtract, const double Zgamma_EW2) const
Definition: EWSM_Output.cpp:583
StandardModel::EW3
Three-loop of .
Definition: StandardModel.h:501
QCD::getMtpole
double getMtpole() const
A get method to access the pole mass of the top quark.
Definition: QCD.h:588
EWSMTwoLoopQCD::DeltaRho
double DeltaRho(const double Mw_i) const
Leading two-loop QCD contribution of to , denoted as .
Definition: EWSMTwoLoopQCD.cpp:38
EWSMThreeLoopEW2QCD::DeltaR_rem
double DeltaR_rem(const double Mw_i) const
Remainder contribution of to , denoted as .
Definition: EWSMThreeLoopEW2QCD.cpp:63
StandardModel::getMyTwoLoopEW
EWSMTwoLoopEW * getMyTwoLoopEW() const
Definition: StandardModel.h:1010
EWSMTwoLoopEW::tau_2
double tau_2() const
The function .
Definition: EWSMTwoLoopEW.cpp:151
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
StandardModel::getMyApproximateFormulae
EWSMApproximateFormulae * getMyApproximateFormulae() const
A get method to retrieve the member pointer of type EWSMApproximateFormulae.
Definition: StandardModel.h:979
EWSMApproximateFormulae::DeltaR_TwoLoopEW_rem
double DeltaR_TwoLoopEW_rem(const double Mw_i) const
.
Definition: EWSMApproximateFormulae.cpp:252
EWSMTwoLoopEW::DeltaRho
double DeltaRho(const double Mw_i) const
Leading two-loop contribution of to , denoted as .
Definition: EWSMTwoLoopEW.cpp:54
StandardModel::resumKappaZ
double resumKappaZ(const double DeltaRho[orders_EW_size], const double deltaKappa_rem[orders_EW_size], const double DeltaRbar_rem, const bool bool_Zbb) const
A method to compute the real part of the effetvive coupling from , and .
Definition: StandardModel.cpp:2027
StandardModel::getMz
double getMz() const
A get method to access the mass of the boson .
Definition: StandardModel.h:721
EWSMOneLoopEW::DeltaR_rem
double DeltaR_rem(const double Mw_i) const
Remainder contribution of to , denoted as .
Definition: EWSMOneLoopEW.cpp:49
StandardModel::getFlagKappaZ
std::string getFlagKappaZ() const
A method to retrieve the model flag KappaZ.
Definition: StandardModel.h:682
EWSMThreeLoopQCD::DeltaR_rem
double DeltaR_rem(const double Mw_i) const
Remainder contribution of to , denoted as .
Definition: EWSMThreeLoopQCD.cpp:48
StandardModel::getFlagMw
std::string getFlagMw() const
A method to retrieve the model flag Mw.
Definition: StandardModel.h:662
StandardModel::EW2QCD1
Three-loop of .
Definition: StandardModel.h:500
StandardModel::Mw
virtual double Mw() const
The SM prediction for the -boson mass in the on-shell scheme, .
Definition: StandardModel.cpp:970
EWSMcache::alsMt
double alsMt() const
The strong coupling at NNLO.
Definition: EWSMcache.h:378
StandardModel::kappaZ_f
virtual gslpp::complex kappaZ_f(const Particle f) const
The effective leptonic neutral-current coupling in the SM.
Definition: StandardModel.cpp:1644
StandardModel::EW1QCD2
Three-loop of .
Definition: StandardModel.h:498
EWSMcache::Xt_alpha
double Xt_alpha(const double Mw_i) const
The quantity with the coupling .
Definition: EWSMcache.h:355
StandardModel::getAle
double getAle() const
A get method to retrieve the fine-structure constant .
Definition: StandardModel.h:748
StandardModel::orders_EW_size
The size of this enum.
Definition: StandardModel.h:502
EWSMTwoLoopEW::DeltaR_rem
double DeltaR_rem(const double Mw_i) const
Remainder contribution of to , denoted as .
Definition: EWSMTwoLoopEW.cpp:81
StandardModel::getMyThreeLoopEW2QCD
EWSMThreeLoopEW2QCD * getMyThreeLoopEW2QCD() const
Definition: StandardModel.h:1000
StandardModel::getLeptons
Particle getLeptons(const QCD::lepton p) const
A get method to retrieve the member object of a lepton.
Definition: StandardModel.h:712
StandardModel::orders_EW
orders_EW
An enumerated type representing perturbative orders of radiative corrections to EW precision observab...
Definition: StandardModel.h:495