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

Copyright

GNU General Public License

Definition at line 22 of file EWSM_Output.h.

Public Member Functions
	EWSM_Output (const StandardModel &SM_in)
	Constructor.
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 const 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\).
virtual const double	AFB (const Particle f) const
gslpp::complex	AH_f (const double tau) const
	Fermionic loop function entering in the calculation of the effective \(Hgg\) and \(H\gamma\gamma\) couplings.
gslpp::complex	AH_W (const double tau) const
	W loop function entering in the calculation of the effective \(H\gamma\gamma\) coupling.
gslpp::complex	AHZga_f (const double tau, const double lambda) const
	Fermionic loop function entering in the calculation of the effective \(HZ\gamma\) coupling.
gslpp::complex	AHZga_W (const double tau, const double lambda) const
	W loop function entering in the calculation of the effective \(HZ\gamma\) coupling.
const double	Ale (double mu, orders order, bool Nf_thr=true) const
	The running electromagnetic coupling \(\alpha_e(\mu)\) in the \(\overline{MS}\) scheme.
const double	ale_OS (const double mu, orders order=FULLNLO) const
	The running electromagnetic coupling \(\alpha(\mu)\) in the on-shell scheme.
virtual const double	alphaMz () const
	The electromagnetic coupling at the \(Z\)-mass scale, \(\alpha(M_Z^2)=\alpha/(1-\Delta\alpha(M_Z^2))\).
virtual const double	alrmoller (const double q2, const double y) const
	The computation of the parity violating asymmetry in Moller scattering.
const double	Als (const double mu, const int Nf_in, const orders order=FULLNLO) const
	Computes the running strong coupling \(\alpha_s(\mu)\) with \(N_f\) active flavours in the \(\overline{\mathrm{MS}}\) scheme. In the cases of LO, NLO and FULLNLO, the coupling is computed with AlsWithInit(). On the other hand, in the cases of NNLO and FULLNNLO, the coupling is computed with AlsWithLambda().
const double	Als (const double mu, const orders order, const bool Nf_thr, const bool qed_flag) const
	The running QCD coupling \(\alpha(\mu)\) in the \(\overline{MS}\) scheme including QED corrections.
const double	Als (const double mu, const orders order=FULLNLO, const bool Nf_thr=true) const
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.
virtual const double	amuon () const
	The computation of the anomalous magnetic moment of the muon \(a_\mu=(g_\mu-2)/2\).
const double	Beta_e (int nm, unsigned int nf) const
	QED beta function coefficients - eq. (36) hep-ph/0512066.
const double	Beta_s (int nm, unsigned int nf) const
	QCD beta function coefficients including QED corrections - eq. (36) hep-ph/0512066.
virtual const double	BrHtobb () const
	The Br \((H\to b \bar{b})\) in the Standard Model.
virtual const double	BrHtocc () const
	The Br \((H\to c \bar{c})\) in the Standard Model.
virtual const double	BrHtogaga () const
	The Br \((H\to \gamma \gamma)\) in the Standard Model.
virtual const double	BrHtogg () const
	The Br \(\(H\to gg)\) in the Standard Model.
virtual const double	BrHtomumu () const
	The Br \((H\to \mu^+ \mu^-)\) in the Standard Model.
virtual const double	BrHtoss () const
	The Br \((H\to s \bar{s})\) in the Standard Model.
virtual const double	BrHtotautau () const
	The Br \((H\to \tau^+ \tau^-)\) in the Standard Model.
virtual const double	BrHtoWWstar () const
	The Br \((H\to W W^*)\) in the Standard Model.
virtual const double	BrHtoZga () const
	The Br \((H\to Z \gamma)\) in the Standard Model.
virtual const double	BrHtoZZstar () const
	The Br \((H\to Z Z^*)\) in the Standard Model.
virtual const double	BrW (const Particle fi, const Particle fj) const
	The branching ratio of the \(W\) boson decaying into a SM fermion pair, \(Br(W\to f_i f_j)\).
const double	c02 () const
	The square of the cosine of the weak mixing angle \(c_0^2\) defined without weak radiative corrections.
virtual bool	CheckFlags () const
	A method to check the sanity of the set of model flags.
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.
bool	checkSMparamsForEWPO ()
	A method to check whether the parameters relevant to the EWPO are updated.
const double	computeBrHto4f () const
	The Br \((H\to 4f)\) in the Standard Model.
const double	computeBrHto4l2 () const
	The Br \((H\to 4l)\) \(l=e,\mu\) in the Standard Model.
const double	computeBrHto4l3 () const
	The Br \((H\to 4l)\) \(l=e,\mu,\tau\) in the Standard Model.
const double	computeBrHto4q () const
	The Br \((H\to 4q)\) in the Standard Model.
const double	computeBrHto4v () const
	The Br \((H\to 4\nu)\) in the Standard Model.
const double	computeBrHtobb () const
	The Br \((H\to bb)\) in the Standard Model.
const double	computeBrHtocc () const
	The Br \((H\to cc)\) in the Standard Model.
const double	computeBrHtoevmuv () const
	The Br \((H\to e \nu \mu \nu)\) in the Standard Model.
const double	computeBrHtogaga () const
	The Br \((H\to\gamma\gamma)\) in the Standard Model.
const double	computeBrHtogg () const
	The Br \((H\to gg)\) in the Standard Model.
const double	computeBrHtollvv2 () const
	The Br \((H\to l^+ l^- \nu \nu)\) \(l=e,\mu\) in the Standard Model.
const double	computeBrHtollvv3 () const
	The Br \((H\to l^+ l^- \nu \nu)\) \(l=e,\mu,\tau\) in the Standard Model.
const double	computeBrHtomumu () const
	The Br \((H\to \mu\mu)\) in the Standard Model.
const double	computeBrHtoss () const
	The Br \((H\to ss)\) in the Standard Model.
const double	computeBrHtotautau () const
	The Br \((H\to \tau\tau)\) in the Standard Model.
const double	computeBrHtoWW () const
	The Br \((H\to WW)\) in the Standard Model.
const double	computeBrHtoZga () const
	The Br \((H\to Z\gamma)\) in the Standard Model.
const double	computeBrHtoZZ () const
	The Br \((H\to ZZ)\) in the Standard Model.
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.
void	ComputeDeltaRho (const double Mw_i, double DeltaRho[orders_EW_size]) const
	A method to collect \(\Delta\rho\) computed via subclasses.
const double	computeGammaHgaga_tt () const
	The top loop contribution to \(H\to\gamma\gamma\) in the Standard Model.
const double	computeGammaHgaga_tW () const
	The mixed \(t-W\) loop contribution to \(H\to\gamma\gamma\) in the Standard Model.
const double	computeGammaHgaga_WW () const
	The \(W\) loop contribution to \(H\to\gamma\gamma\) in the Standard Model.
const double	computeGammaHgg_bb () const
	The bottom loop contribution to \(H\to gg\) in the Standard Model.
const double	computeGammaHgg_tb () const
	The top-bottom interference contribution to \(H\to gg\) in the Standard Model.
const double	computeGammaHgg_tt () const
	The top loop contribution to \(H\to gg\) in the Standard Model.
const double	computeGammaHTotal () const
	The Higgs total width in the Standard Model.
const double	computeGammaHZga_tt () const
	The top loop contribution to \(H\to Z\gamma\) in the Standard Model.
const double	computeGammaHZga_tW () const
	The mixed \(t-W\) loop contribution to \(H\to Z\gamma\) in the Standard Model.
const 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. [Heinemeyer:2013tqa].
const double	computeSigmabbH (const double sqrt_s) const
	The bbH production cross section in the Standard Model.
const double	computeSigmaggH (const double sqrt_s) const
	The ggH cross section in the Standard Model.
const double	computeSigmaggH_bb (const double sqrt_s) const
	The square of the bottom-quark contribution to the ggH cross section in the Standard Model.
const double	computeSigmaggH_tb (const double sqrt_s) const
	The top-bottom interference contribution to the ggH cross section in the Standard Model.
const double	computeSigmaggH_tt (const double sqrt_s) const
	The square of the top-quark contribution to the ggH cross section in the Standard Model.
const double	computeSigmatHq (const double sqrt_s) const
	The tHq production cross section in the Standard Model.
const double	computeSigmattH (const double sqrt_s) const
	The ttH production cross section in the Standard Model.
const double	computeSigmaVBF (const double sqrt_s) const
	The VBF cross section in the Standard Model.
const double	computeSigmaWF (const double sqrt_s) const
	The W fusion contribution \(\sigma_{WF}\) to higgs-production cross section in the Standard Model.
const double	computeSigmaWH (const double sqrt_s) const
	The WH production cross section in the Standard Model.
const double	computeSigmaZF (const double sqrt_s) const
	The Z fusion contribution \(\sigma_{ZF}\) to higgs-production cross section in the Standard Model.
const double	computeSigmaZH (const double sqrt_s) const
	The ZH production cross section in the Standard Model.
const double	computeSigmaZWF (const double sqrt_s) const
	The Z W interference fusion contribution \(\sigma_{ZWF}\) to higgs-production cross section in the Standard Model.
virtual const double	cW2 () const
virtual const 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\).
virtual const double	Dalpha5hMz () const
	The 5-quark contribution to the running of the em constant to the \(Z\) pole. \(\Delta\alpha_{had}^{(5)}(M_Z)\).
const double	DeltaAlpha () const
	The total corrections to the electromagnetic coupling \(\alpha\) at the \(Z\)-mass scale, denoted as \(\Delta\alpha(M_Z^2)\).
const 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)\).
const double	DeltaAlphaLepton (const double s) const
	Leptonic contribution to the electromagnetic coupling \(\alpha\), denoted as \(\Delta\alpha_{\mathrm{lept}}(s)\).
const double	DeltaAlphaTop (const double s) const
	Top-quark contribution to the electromagnetic coupling \(\alpha\), denoted as \(\Delta\alpha_{\mathrm{top}}(s)\).
virtual const gslpp::complex	deltaKappaZ_f (const Particle f) const
	Flavour non-universal vertex corrections to \(\kappa_Z^l\), denoted by \(\Delta\kappa_Z^l\).
virtual const double	DeltaR () const
	The SM prediction for \(\Delta r\) derived from that for the \(W\) boson mass.
virtual const double	DeltaRbar () const
	The SM prediction for \(\Delta \overline{r}\) derived from that for the \(W\)-boson mass.
virtual const gslpp::complex	deltaRhoZ_f (const Particle f) const
	Flavour non-universal vertex corrections to \(\rho_Z^l\), denoted by \(\Delta\rho_Z^l\).
virtual const double	eeffAFBbottom (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffAFBcharm (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffAFBe (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffAFBetsub (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffAFBmu (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffAFBstrange (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffAFBtau (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffRbottom (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffRcharm (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffRelectron (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffRelectrontsub (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffRmuon (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffRstrange (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffRtau (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffsigma (const Particle f, const double pol_e, const double pol_p, const double s, const double cosmin, const double cosmax) const
virtual const double	eeffsigmaBottom (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffsigmaCharm (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffsigmaE (const double pol_e, const double pol_p, const double s) const
const double	eeffsigmaEbin (const double pol_e, const double pol_p, const double s, const double cosmin, const double cosmax) const
virtual const double	eeffsigmaEtsub (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffsigmaHadron (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffsigmaMu (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffsigmaStrange (const double pol_e, const double pol_p, const double s) const
virtual const double	eeffsigmaTau (const double pol_e, const double pol_p, const double s) const
virtual const double	epsilon1 () const
	The SM contribution to the epsilon parameter \(\varepsilon_1\).
virtual const double	epsilon2 () const
	The SM contribution to the epsilon parameter \(\varepsilon_2\).
virtual const double	epsilon3 () const
	The SM contribution to the epsilon parameter \(\varepsilon_3\).
virtual const double	epsilonb () const
	The SM contribution to the epsilon parameter \(\varepsilon_b\).
gslpp::complex	f_triangle (const double tau) const
	Loop function entering in the calculation of the effective \(Hgg\) and \(H\gamma\gamma\) couplings.
gslpp::complex	g_triangle (const double tau) const
	Loop function entering in the calculation of the effective \(HZ\gamma\) coupling.
virtual const gslpp::complex	gA_f (const Particle f) const
	The effective leptonic neutral-current axial-vector coupling \(g_A^l\) in the SM.
virtual const double	Gamma_had () const
	The hadronic decay width of the \(Z\) boson, \(\Gamma_{h}\).
virtual const double	Gamma_inv () const
	The invisible partial decay width of the \(Z\) boson, \(\Gamma_{\mathrm{inv}}\).
virtual const double	Gamma_muon () const
	The computation of the muon decay.
virtual const double	Gamma_tau_l_nunu (const Particle l) const
	The computation of the leptonic tau decays.
virtual const double	Gamma_Z () const
	The total decay width of the \(Z\) boson, \(\Gamma_Z\).
virtual const double	GammaHtobb () const
	The \(\Gamma(H\to b \bar{b})\) in the Standard Model.
virtual const double	GammaHtocc () const
	The \(\Gamma(H\to c \bar{c})\) in the Standard Model.
virtual const double	GammaHtogaga () const
	The \(\Gamma(H\to \gamma \gamma)\) in the Standard Model.
virtual const double	GammaHtogg () const
	The \(\Gamma(H\to gg)\) in the Standard Model.
virtual const double	GammaHtomumu () const
	The \(\Gamma(H\to \mu^+ \mu^-)\) in the Standard Model.
virtual const double	GammaHtoss () const
	The \(\Gamma(H\to s \bar{s})\) in the Standard Model.
virtual const double	GammaHTot () const
	The total Higgs width \(\Gamma(H)\) in the Standard Model.
virtual const double	GammaHtotautau () const
	The \(\Gamma(H\to \tau^+ \tau^-)\) in the Standard Model.
virtual const double	GammaHtoWWstar () const
	The \(\Gamma(H\to W W^*)\) in the Standard Model.
virtual const double	GammaHtoZga () const
	The \(\Gamma(H\to Z \gamma)\) in the Standard Model.
virtual const double	GammaHtoZZstar () const
	The \(\Gamma(H\to Z Z^*)\) in the Standard Model.
virtual const double	GammaW () const
	The total width of the \(W\) boson, \(\Gamma_W\).
virtual const double	GammaW (const Particle fi, const Particle fj) const
	A partial decay width of the \(W\) boson decay into a SM fermion pair.
virtual const double	GammaZ (const Particle f) const
	The \(Z\to \ell\bar{\ell}\) partial decay width, \(\Gamma_\ell\).
virtual const double	gAnue () const
	The effective (muon) neutrino-electron axial-vector coupling: gAnue.
const double	getAle () const
	A get method to retrieve the fine-structure constant \(\alpha\).
const double	getAlsMz () const
	A get method to access the value of \(\alpha_s(M_Z)\).
virtual const double	getCBd () const
	The ratio of the absolute value of the $B_d$ mixing amplitude over the Standard Model value.
virtual const double	getCBs () const
	The ratio of the absolute value of the $B_s$ mixing amplitude over the Standard Model value.
virtual const double	getCCC1 () const
	A virtual implementation for the RealWeakEFTCC class.
virtual const double	getCCC2 () const
	A virtual implementation for the RealWeakEFTCC class.
virtual const double	getCCC3 () const
	A virtual implementation for the RealWeakEFTCC class.
virtual const double	getCCC4 () const
	A virtual implementation for the RealWeakEFTCC class.
virtual const double	getCCC5 () const
	A virtual implementation for the RealWeakEFTCC class.
virtual const double	getCDMK () const
	The ratio of the real part of the $K$ mixing amplitude over the Standard Model value.
virtual const double	getCepsK () const
	The ratio of the imaginary part of the $K$ mixing amplitude over the Standard Model value.
const CKM &	getCKM () const
	A get method to retrieve the member object of type CKM.
const 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)\).
const double	getDelGammaWlv () const
	A get method to retrieve the theoretical uncertainty in \(\Gamma_W_{l\nu}\), denoted as \(\delta\,\Gamma_W_{l\nu}\).
const double	getDelGammaWqq () const
	A get method to retrieve the theoretical uncertainty in \(\Gamma_W_{qq}\), denoted as \(\delta\,\Gamma_W_{qq}\).
const double	getDelGammaZ () const
	A get method to retrieve the theoretical uncertainty in \(\Gamma_Z\), denoted as \(\delta\,\Gamma_Z\).
const double	getDelMw () const
	A get method to retrieve the theoretical uncertainty in \(M_W\), denoted as \(\delta\,M_W\).
const double	getDelR0b () const
	A get method to retrieve the theoretical uncertainty in \(R_b^0\), denoted as \(\delta\,R_b^0\).
const double	getDelR0c () const
	A get method to retrieve the theoretical uncertainty in \(R_c^0\), denoted as \(\delta\,R_c^0\).
const double	getDelR0l () const
	A get method to retrieve the theoretical uncertainty in \(R_l^0\), denoted as \(\delta\,R_l^0\).
const double	getDelSigma0H () const
	A get method to retrieve the theoretical uncertainty in \(\sigma_{Hadron}^0\), denoted as \(\delta\,\sigma_{Hadron}^0\).
const 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}\).
const 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}\).
const 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}\).
const std::string	getFlagKappaZ () const
	A method to retrieve the model flag KappaZ.
const std::string	getFlagMw () const
	A method to retrieve the model flag Mw.
const std::string	getFlagRhoZ () const
	A method to retrieve the model flag RhoZ.
const Flavour &	getFlavour () const
const double	getGF () const
	A get method to retrieve the Fermi constant \(G_\mu\).
const int	getIterationNo () const
const Particle &	getLeptons (const QCD::lepton p) const
	A get method to retrieve the member object of a lepton.
virtual StandardModelMatching &	getMatching () const
	A get method to access the member reference of type StandardModelMatching.
virtual const double	getMHl () const
	A get method to retrieve the Higgs mass \(m_h\).
virtual const double	getmq (const QCD::quark q, const double mu) const
	The MSbar running quark mass computed at NLO.
const double	getMuw () const
	A get method to retrieve the matching scale \(\mu_W\) around the weak scale.
const double	getMw () const
	A get method to access the input value of the mass of the \(W\) boson \(M_W\).
EWSMApproximateFormulae *	getMyApproximateFormulae () const
	A get method to retrieve the member pointer of type EWSMApproximateFormulae.
EWSMcache *	getMyEWSMcache () const
	A get method to retrieve the member pointer of type EWSMcache.
LeptonFlavour *	getMyLeptonFlavour () const
EWSMOneLoopEW *	getMyOneLoopEW () const
	A get method to retrieve the member pointer of type EWSMOneLoopEW,.
EWSMThreeLoopEW *	getMyThreeLoopEW () const
EWSMThreeLoopEW2QCD *	getMyThreeLoopEW2QCD () const
EWSMThreeLoopQCD *	getMyThreeLoopQCD () const
EWSMTwoFermionsLEP2 *	getMyTwoFermionsLEP2 () const
	A get method to retrieve the member pointer of type EWSMTwoFermionsLEP2.
EWSMTwoLoopEW *	getMyTwoLoopEW () const
EWSMTwoLoopQCD *	getMyTwoLoopQCD () const
const double	getMz () const
	A get method to access the mass of the \(Z\) boson \(M_Z\).
virtual const double	getPhiBd () const
	Half the relative phase of the $B_d$ mixing amplitude w.r.t. the Standard Model one.
virtual const double	getPhiBs () const
	Half the relative phase of the $B_s$ mixing amplitude w.r.t. the Standard Model one.
virtual const StandardModel &	getTrueSM () const
const gslpp::matrix< gslpp::complex >	getUPMNS () const
	A get method to retrieve the object of the PMNS matrix.
const gslpp::matrix< gslpp::complex >	getVCKM () const
	A get method to retrieve the CKM matrix.
const gslpp::matrix< gslpp::complex > &	getYd () const
	A get method to retrieve the Yukawa matrix of the down-type quarks, \(Y_d\).
const gslpp::matrix< gslpp::complex > &	getYe () const
	A get method to retrieve the Yukawa matrix of the charged leptons, \(Y_e\).
const gslpp::matrix< gslpp::complex > &	getYn () const
	A get method to retrieve the Yukawa matrix of the neutrinos, \(Y_\nu\).
const gslpp::matrix< gslpp::complex > &	getYu () const
	A get method to retrieve the Yukawa matrix of the up-type quarks, \(Y_u\).
virtual const double	gLnuN2 () const
	The effective neutrino nucleon LH coupling: gLnuN2.
virtual const double	gRnuN2 () const
	The effective neutrino nucleon RH coupling: gRnuN2.
virtual const gslpp::complex	gV_f (const Particle f) const
	The effective leptonic neutral-current vector coupling \(g_V^l\) in the SM.
virtual const double	gVnue () const
	The effective (muon) neutrino-electron vector coupling: gVnue.
gslpp::complex	I_triangle_1 (const double tau, const double lambda) const
	Loop function entering in the calculation of the effective \(HZ\gamma\) coupling.
gslpp::complex	I_triangle_2 (const double tau, const double lambda) const
	Loop function entering in the calculation of the effective \(HZ\gamma\) coupling.
virtual bool	Init (const std::map< std::string, double > &DPars)
	A method to initialize the model parameters.
virtual bool	InitializeModel ()
	A method to initialize the model.
const double	intMLL2eeeeus2 (const double s, const double t0, const double t1) const
const double	intMLR2eeeets2 (const double s, const double t0, const double t1) const
const double	intMLRtilde2eeeest2 (const double s, const double t0, const double t1) const
const double	intMRR2eeeeus2 (const double s, const double t0, const double t1) const
const bool	IsFlagNoApproximateGammaZ () const
	A method to retrieve the model flag NoApproximateGammaZ.
const bool	IsFlagWithoutNonUniversalVC () const
	A method to retrieve the model flag WithoutNonUniversalVC.
const bool	isSMSuccess () const
	A get method to retrieve the success status of the Standard Model update and matching.
virtual const gslpp::complex	kappaZ_f (const Particle f) const
	The effective leptonic neutral-current coupling \(\kappa_Z^l\) in the SM.
virtual const double	LEP2AFBbottom (const double s) const
virtual const double	LEP2AFBcharm (const double s) const
virtual const double	LEP2AFBe (const double s) const
virtual const double	LEP2AFBmu (const double s) const
virtual const double	LEP2AFBtau (const double s) const
virtual const double	LEP2dsigmadcosBinE (const double s, const double cos, const double cosmin, const double cosmax) const
virtual const double	LEP2dsigmadcosBinMu (const double s, const double cos, const double cosmin, const double cosmax) const
virtual const double	LEP2dsigmadcosBinTau (const double s, const double cos, const double cosmin, const double cosmax) const
virtual const double	LEP2dsigmadcosE (const double s, const double cos) const
virtual const double	LEP2dsigmadcosMu (const double s, const double cos) const
virtual const double	LEP2dsigmadcosTau (const double s, const double cos) const
virtual const double	LEP2Rbottom (const double s) const
virtual const double	LEP2Rcharm (const double s) const
virtual const double	LEP2sigmaBottom (const double s) const
virtual const double	LEP2sigmaCharm (const double s) const
virtual const double	LEP2sigmaE (const double s) const
virtual const double	LEP2sigmaHadron (const double s) const
virtual const double	LEP2sigmaMu (const double s) const
virtual const double	LEP2sigmaTau (const double s) const
const double	MLL2eeff (const Particle f, const double s, const double t) const
const double	MLR2eeff (const Particle f, const double s) const
const double	MRL2eeff (const Particle f, const double s) const
const double	MRR2eeff (const Particle f, const double s, const double t) const
virtual const double	Mw () const
	The SM prediction for the \(W\)-boson mass in the on-shell scheme, \(M_{W,\mathrm{SM}}\).
const double	Mw_tree () const
	The tree-level mass of the \(W\) boson, \(M_W^{\mathrm{tree}}\).
const 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.
const 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.
double	Mzbar () const
	The \(Z\)-boson mass \(\overline{M}_Z\) in the complex-pole/fixed-width scheme.
virtual const double	N_nu () const
	The number of neutrinos obtained indirectly from the measurements at the Z pole, \(N_{\nu}\).
virtual bool	PostUpdate ()
	The post-update method for StandardModel.
virtual bool	PreUpdate ()
	The pre-update method for StandardModel.
virtual const double	Qwemoller (const double q2, const double y) const
	The computation of the electron's weak charge.
virtual const double	Qwn () const
	The computation of the neutron weak charge: Qwn.
virtual const double	Qwp () const
	The computation of the proton weak charge: Qwp.
virtual const double	R0_f (const Particle f) const
	The ratio \(R_\ell^0=\Gamma(Z\to {\rm hadrons})/\Gamma(Z\to \ell^+ \ell^-)\).
virtual const double	R_inv () const
	The ratio of the invisible and leptonic (electron) decay widths of the \(Z\) boson, \(R_{inv}\).
virtual const 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}\).
virtual const gslpp::complex	rhoZ_f (const Particle f) const
	The effective leptonic neutral-current coupling \(\rho_Z^l\) in the SM.
virtual const double	Ruc () const
virtual const double	RWc () const
	The ratio \(R_{W,c)=\Gamma(W\to c + X)/\Gamma(W\to had)\).
virtual const double	RWlilj (const Particle li, const Particle lj) const
	The lepton universality ratio \(R_{W,l_i/l_j)=\Gamma(W\to l_i \nu_i)/\Gamma(W\to l_j \nu_j)\).
virtual const double	RZlilj (const Particle li, const Particle lj) const
	The lepton universality ratio \(R_{Z,l_i/l_j)=\Gamma(Z\to l_i^+ l_i^-)/\Gamma(Z\to l_j^+ l_j^-)\).
const double	s02 () const
	The square of the sine of the weak mixing angle \(s_0^2\) defined without weak radiative corrections.
void	setCKM (const CKM &CKMMatrix)
	A set method to change the CKM matrix.
virtual bool	setFlag (const std::string name, const bool value)
	A method to set a flag of StandardModel.
void	setFlagCacheInStandardModel (bool FlagCacheInStandardModel)
	A set method to change the model flag CacheInStandardModel of StandardModel.
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.
void	setRequireCKM (bool requireCKM)
	A set method to change the value of requireCKM.
void	setSMSuccess (bool success) const
	A set method to change the success status of the Standard Model update and matching.
void	setYd (const gslpp::matrix< gslpp::complex > &Yd)
	A set method to set the Yukawa matrix of the down-type quarks, \(Y_d\).
void	setYe (const gslpp::matrix< gslpp::complex > &Ye)
	A set method to set the Yukawa matrix of the charged leptons, \(Y_e\).
void	setYu (const gslpp::matrix< gslpp::complex > &Yu)
	A set method to set the Yukawa matrix of the up-type quarks, \(Y_u\).
virtual const double	sigma0_had () const
	The hadronic cross section for \(e^+e^- \to Z \to \mathrm{hadrons}\) at the \(Z\)-pole, \(\sigma_h^0\).
virtual const double	SigmaeeHee (const double sqrt_s, const double Pe, const double Pp) const
	The \(\sigma(e^+ e^- \to e^+ e^- H)\) in the Standard Model.
virtual const double	SigmaeeHvv (const double sqrt_s, const double Pe, const double Pp) const
	The \(\sigma(e^+ e^- \to \nu \bar{\nu} H)\) in the Standard Model.
virtual const double	SigmaeeZH (const double sqrt_s, const double Pe, const double Pp) const
	The \(\sigma(e^+ e^- \to Z H)\) in the Standard Model.
virtual const 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.
	StandardModel ()
	The default constructor.
const double	sW2 () const
virtual const 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\).
const double	sW2_MSbar_Approx () const
	The (approximated formula for the) square of the sine of the weak mixing angle in the MSbar scheme, denoted as \(\hat{s}_{W}^2\). See: PDG 22, R.L. Workman et al. (Particle Data Group), Prog. Theor. Exp. Phys. 2022, 083C01 (2022)
const double	sW2_ND () const
	The square of the sine of the weak mixing angle in the MSbar-ND scheme (w/o decoupling $\alpha\ln(m_t/M_Z)$ terms), denoted as \(\hat{s}_{ND}^2\). See: PDG 22, R.L. Workman et al. (Particle Data Group), Prog. Theor. Exp. Phys. 2022, 083C01 (2022) (eq. 10.13a/10.13b)
virtual const double	TauLFU_gmuge () const
	The computation of the LFU ratio \(g_\mu/ g_e \).
virtual const double	TauLFU_gtauge () const
	The computation of the LFU ratio \(g_\tau/ g_e \).
virtual const double	TauLFU_gtaugmu () const
	The computation of the LFU ratio \(g_\tau/ g_\mu \).
virtual const double	TauLFU_gtaugmuK () const
	The computation of the LFU ratio \(\left(g_\tau/ g_\mu\right)_K \).
virtual const double	TauLFU_gtaugmuPi () const
	The computation of the LFU ratio \(\left(g_\tau/ g_\mu\right)_\pi \).
virtual const double	ThetaLnuN () const
	The effective neutrino nucleon LH parameter: ThetaLnuN.
virtual const double	ThetaRnuN () const
	The effective neutrino nucleon RH parameter: ThetaRnuN.
const double	tovers2 (const double cosmin, const double cosmax) const
const double	uovers2 (const double cosmin, const double cosmax) const
virtual bool	Update (const std::map< std::string, double > &DPars)
	The update method for StandardModel.
const double	v () const
	The Higgs vacuum expectation value.
virtual	~StandardModel ()
	The default destructor.
Public Member Functions inherited from QCD
const double	AboveTh (const double mu) const
	The active flavour threshold above the scale \(\mu\) as defined in QCD::Thresholds().
void	addParameters (std::vector< std::string > params_i)
	A method to add parameters that are specific to only one set of observables.
const double	Als (const double mu, const int Nf_in, const orders order=FULLNLO) const
	Computes the running strong coupling \(\alpha_s(\mu)\) with \(N_f\) active flavours in the \(\overline{\mathrm{MS}}\) scheme. In the cases of LO, NLO and FULLNLO, the coupling is computed with AlsWithInit(). On the other hand, in the cases of NNLO and FULLNNLO, the coupling is computed with AlsWithLambda().
const double	Als (const double mu, const orders order=FULLNLO, const bool Nf_thr=true) const
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\).
const double	AlsByOrder (const double mu, const int Nf_in, const orders order=FULLNLO) const
const double	AlsByOrder (const double mu, const orders order=FULLNLO, bool Nf_thr=true) const
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().
const double	AlsWithInit (const double mu, const double alsi, const double mu_i, const int nf, 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\).
const 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}\).
const double	BelowTh (const double mu) const
	The active flavour threshold below the scale \(\mu\) as defined in QCD::Thresholds().
const double	Beta0 (const double nf) const
	The \(\beta_0(n_f)\) coefficient for a certain number of flavours \(n_f\).
const double	Beta1 (const double nf) const
	The \(\beta_1(n_f)\) coefficient for a certain number of flavours \(n_f\).
const double	Beta2 (const double nf) const
	The \(\beta_2(n_f)\) coefficient for a certain number of flavours \(n_f\).
const double	Beta3 (const double nf) const
	The \(\beta_3(n_f)\) coefficient for a certain number of flavours \(n_f\).
void	CacheShift (double cache[][5], int n) const
	A member used to manage the caching for this class.
void	CacheShift (int cache[][5], int n) const
const orders	FullOrder (orders order) const
	Return the FULLORDER enum corresponding to order.
const double	Gamma0 (const double nf) const
	The \(\gamma_0\) coefficient used to compute the running of a mass.
const double	Gamma1 (const double nf) const
	The \(\gamma_1\) coefficient used to compute the running of a mass.
const double	Gamma2 (const double nf) const
	The \(\gamma_2\) coefficient used to compute the running of a mass.
const double	getAlsM () const
	A get method to access the value of \(\alpha_s(M_{\alpha_s})\).
const 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.
const BParameter &	getBBd_subleading () const
	For getting the subleading bag parameters \(R_2 - R_3\) in \(\Delta b = 2\) process in the \(B_d\) meson system.
const 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.
const BParameter &	getBBs_subleading () const
	For getting the subleading bag parameters \(R_2 - R_3\) in \(\Delta b = 2\) process in the \(B_s\) meson system.
const 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.
const 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.
const BParameter &	getBKd1 () const
const BParameter &	getBKd3 () const
const double	getCF () const
	A get method to access the Casimir factor of QCD.
const double	getMAls () const
	A get method to access the mass scale \(M_{\alpha_s}\) at which the strong coupling constant measurement is provided.
const Meson &	getMesons (const QCD::meson m) const
	A get method to access a meson as an object of the type Meson.
const double	getMtpole () const
	A get method to access the pole mass of the top quark.
const double	getMub () const
	A get method to access the threshold between five- and four-flavour theory in GeV.
const double	getMuc () const
	A get method to access the threshold between four- and three-flavour theory in GeV.
const double	getMut () const
	A get method to access the threshold between six- and five-flavour theory in GeV.
const double	getNc () const
	A get method to access the number of colours \(N_c\).
const double	getOptionalParameter (std::string name) const
	A method to get parameters that are specific to only one set of observables.
const Particle &	getQuarks (const QCD::quark q) const
	A get method to access a quark as an object of the type Particle.
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.
void	initializeBParameter (std::string name_i) const
	A method to initialize B Parameter and the corresponding meson.
void	initializeMeson (QCD::meson meson_i) const
	A method to initialize a meson.
bool	isQCDsuccess () const
	A getter for the QCDsuccess flag.
const double	logLambda (const double nf, orders order) const
	Computes \(\ln\Lambda_\mathrm{QCD}\) with nf flavours in GeV.
const double	Mbar2Mp (const double mbar, const quark q, const orders order=FULLNNLO) const
	Converts the \(\overline{\mathrm{MS}}\) mass \(m(m)\) to the pole mass.
const double	Mofmu2Mbar (const double m, const double mu, const quark q) const
	Converts a quark running mass at an arbitrary scale to the corresponding \(\overline{\mathrm{MS}}\) mass \(m(m)\).
const double	Mp2Mbar (const double mp, const quark q, orders order=FULLNNLO) const
	Converts a quark pole mass to the corresponding \(\overline{\mathrm{MS}}\) mass \(m(m)\).
const double	Mrun (const double mu, const double m, const quark q, const orders order=FULLNNLO) const
	Computes a running quark mass \(m(\mu)\) from \(m(m)\).
const double	Mrun (const double mu_f, const double mu_i, const double m, const quark q, const orders order=FULLNNLO) const
	Runs a quark mass from \(\mu_i\) to \(\mu_f\).
const 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\).
const double	MS2DRqmass (const double MSbar) const
	Converts a quark mass from the \(\overline{\mathrm{MS}}\) scheme to the \(\overline{\mathrm{DR}}\) scheme.
const 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.
const double	Nf (const double mu) const
	The number of active flavour at scale \(\mu\).
const 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.
const std::string	orderToString (const orders order) const
	Converts an object of the enum type "orders" to the corresponding string.
	QCD ()
	Constructor.
void	setComputemt (bool computemt)
	A set method to change the value of computemt.
void	setMtpole (double mtpole_in)
	A method to set the pole mass of the top quark.
void	setNc (double Nc)
	A set method to change the number of colours \(N_c\).
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.
void	setQuarkMass (const quark q, const double mass)
	A set method to change the mass of a quark.
const double	Thresholds (const int i) const
	For accessing the active flavour threshold scales.
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.
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.
bool	isModelLinearized () const
bool	isModelNPquadratic () 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.
	Model ()
	The default constructor.
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.
void	setModelLinearized (bool linearized=true)
void	setModelName (const std::string name)
	A method to set the name of the model.
void	setModelNPquadratic (bool NPquadratic=true)
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.
virtual	~Model ()
	The default destructor.

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 , D_S , B_D , B_P , B_S , B_C , PHI , K_star , K_star_P , K_S , 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.
static const int	NSMvars = 28
	The number of the model parameters in StandardModel.
static const int	NumSMParamsForEWPO = 35
	The number of the SM parameters that are relevant to the EW precision observables.
static std::string	SMvars [NSMvars]
	A string array containing the labels of the model parameters in StandardModel.
Static Public Attributes inherited from QCD
static const int	NQCDvars = 11
	The number of model parameters in QCD.
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().
Protected Member Functions inherited from StandardModel
const double	AFB_NoISR_l (const QCD::lepton l_flavor, const double s) const
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.
virtual void	computeCKM ()
	The method to compute the CKM matrix.
virtual void	computeYukawas ()
	The method to compute the Yukawas matrix.
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}}\).
const double	getIntegrand_AFBnumeratorWithISR_bottom133 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_bottom167 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_bottom172 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_bottom183 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_bottom189 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_bottom192 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_bottom196 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_bottom200 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_bottom202 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_bottom205 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_bottom207 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_charm133 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_charm167 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_charm172 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_charm183 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_charm189 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_charm192 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_charm196 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_charm200 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_charm202 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_charm205 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_charm207 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_mu130 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_mu136 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_mu161 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_mu172 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_mu183 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_mu189 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_mu192 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_mu196 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_mu200 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_mu202 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_mu205 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_mu207 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_tau130 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_tau136 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_tau161 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_tau172 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_tau183 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_tau189 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_tau192 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_tau196 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_tau200 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_tau202 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_tau205 (double x) const
const double	getIntegrand_AFBnumeratorWithISR_tau207 (double x) const
const double	getIntegrand_dsigmaBox_bottom130 (double x) const
const double	getIntegrand_dsigmaBox_bottom133 (double x) const
const double	getIntegrand_dsigmaBox_bottom136 (double x) const
const double	getIntegrand_dsigmaBox_bottom161 (double x) const
const double	getIntegrand_dsigmaBox_bottom167 (double x) const
const double	getIntegrand_dsigmaBox_bottom172 (double x) const
const double	getIntegrand_dsigmaBox_bottom183 (double x) const
const double	getIntegrand_dsigmaBox_bottom189 (double x) const
const double	getIntegrand_dsigmaBox_bottom192 (double x) const
const double	getIntegrand_dsigmaBox_bottom196 (double x) const
const double	getIntegrand_dsigmaBox_bottom200 (double x) const
const double	getIntegrand_dsigmaBox_bottom202 (double x) const
const double	getIntegrand_dsigmaBox_bottom205 (double x) const
const double	getIntegrand_dsigmaBox_bottom207 (double x) const
const double	getIntegrand_dsigmaBox_charm130 (double x) const
const double	getIntegrand_dsigmaBox_charm133 (double x) const
const double	getIntegrand_dsigmaBox_charm136 (double x) const
const double	getIntegrand_dsigmaBox_charm161 (double x) const
const double	getIntegrand_dsigmaBox_charm167 (double x) const
const double	getIntegrand_dsigmaBox_charm172 (double x) const
const double	getIntegrand_dsigmaBox_charm183 (double x) const
const double	getIntegrand_dsigmaBox_charm189 (double x) const
const double	getIntegrand_dsigmaBox_charm192 (double x) const
const double	getIntegrand_dsigmaBox_charm196 (double x) const
const double	getIntegrand_dsigmaBox_charm200 (double x) const
const double	getIntegrand_dsigmaBox_charm202 (double x) const
const double	getIntegrand_dsigmaBox_charm205 (double x) const
const double	getIntegrand_dsigmaBox_charm207 (double x) const
const double	getIntegrand_dsigmaBox_down130 (double x) const
const double	getIntegrand_dsigmaBox_down133 (double x) const
const double	getIntegrand_dsigmaBox_down136 (double x) const
const double	getIntegrand_dsigmaBox_down161 (double x) const
const double	getIntegrand_dsigmaBox_down167 (double x) const
const double	getIntegrand_dsigmaBox_down172 (double x) const
const double	getIntegrand_dsigmaBox_down183 (double x) const
const double	getIntegrand_dsigmaBox_down189 (double x) const
const double	getIntegrand_dsigmaBox_down192 (double x) const
const double	getIntegrand_dsigmaBox_down196 (double x) const
const double	getIntegrand_dsigmaBox_down200 (double x) const
const double	getIntegrand_dsigmaBox_down202 (double x) const
const double	getIntegrand_dsigmaBox_down205 (double x) const
const double	getIntegrand_dsigmaBox_down207 (double x) const
const double	getIntegrand_dsigmaBox_mu130 (double x) const
const double	getIntegrand_dsigmaBox_mu133 (double x) const
const double	getIntegrand_dsigmaBox_mu136 (double x) const
const double	getIntegrand_dsigmaBox_mu161 (double x) const
const double	getIntegrand_dsigmaBox_mu167 (double x) const
const double	getIntegrand_dsigmaBox_mu172 (double x) const
const double	getIntegrand_dsigmaBox_mu183 (double x) const
const double	getIntegrand_dsigmaBox_mu189 (double x) const
const double	getIntegrand_dsigmaBox_mu192 (double x) const
const double	getIntegrand_dsigmaBox_mu196 (double x) const
const double	getIntegrand_dsigmaBox_mu200 (double x) const
const double	getIntegrand_dsigmaBox_mu202 (double x) const
const double	getIntegrand_dsigmaBox_mu205 (double x) const
const double	getIntegrand_dsigmaBox_mu207 (double x) const
const double	getIntegrand_dsigmaBox_strange130 (double x) const
const double	getIntegrand_dsigmaBox_strange133 (double x) const
const double	getIntegrand_dsigmaBox_strange136 (double x) const
const double	getIntegrand_dsigmaBox_strange161 (double x) const
const double	getIntegrand_dsigmaBox_strange167 (double x) const
const double	getIntegrand_dsigmaBox_strange172 (double x) const
const double	getIntegrand_dsigmaBox_strange183 (double x) const
const double	getIntegrand_dsigmaBox_strange189 (double x) const
const double	getIntegrand_dsigmaBox_strange192 (double x) const
const double	getIntegrand_dsigmaBox_strange196 (double x) const
const double	getIntegrand_dsigmaBox_strange200 (double x) const
const double	getIntegrand_dsigmaBox_strange202 (double x) const
const double	getIntegrand_dsigmaBox_strange205 (double x) const
const double	getIntegrand_dsigmaBox_strange207 (double x) const
const double	getIntegrand_dsigmaBox_tau130 (double x) const
const double	getIntegrand_dsigmaBox_tau133 (double x) const
const double	getIntegrand_dsigmaBox_tau136 (double x) const
const double	getIntegrand_dsigmaBox_tau161 (double x) const
const double	getIntegrand_dsigmaBox_tau167 (double x) const
const double	getIntegrand_dsigmaBox_tau172 (double x) const
const double	getIntegrand_dsigmaBox_tau183 (double x) const
const double	getIntegrand_dsigmaBox_tau189 (double x) const
const double	getIntegrand_dsigmaBox_tau192 (double x) const
const double	getIntegrand_dsigmaBox_tau196 (double x) const
const double	getIntegrand_dsigmaBox_tau200 (double x) const
const double	getIntegrand_dsigmaBox_tau202 (double x) const
const double	getIntegrand_dsigmaBox_tau205 (double x) const
const double	getIntegrand_dsigmaBox_tau207 (double x) const
const double	getIntegrand_dsigmaBox_up130 (double x) const
const double	getIntegrand_dsigmaBox_up133 (double x) const
const double	getIntegrand_dsigmaBox_up136 (double x) const
const double	getIntegrand_dsigmaBox_up161 (double x) const
const double	getIntegrand_dsigmaBox_up167 (double x) const
const double	getIntegrand_dsigmaBox_up172 (double x) const
const double	getIntegrand_dsigmaBox_up183 (double x) const
const double	getIntegrand_dsigmaBox_up189 (double x) const
const double	getIntegrand_dsigmaBox_up192 (double x) const
const double	getIntegrand_dsigmaBox_up196 (double x) const
const double	getIntegrand_dsigmaBox_up200 (double x) const
const double	getIntegrand_dsigmaBox_up202 (double x) const
const double	getIntegrand_dsigmaBox_up205 (double x) const
const double	getIntegrand_dsigmaBox_up207 (double x) const
const double	getIntegrand_sigmaWithISR_bottom130 (double x) const
const double	getIntegrand_sigmaWithISR_bottom133 (double x) const
const double	getIntegrand_sigmaWithISR_bottom136 (double x) const
const double	getIntegrand_sigmaWithISR_bottom161 (double x) const
const double	getIntegrand_sigmaWithISR_bottom167 (double x) const
const double	getIntegrand_sigmaWithISR_bottom172 (double x) const
const double	getIntegrand_sigmaWithISR_bottom183 (double x) const
const double	getIntegrand_sigmaWithISR_bottom189 (double x) const
const double	getIntegrand_sigmaWithISR_bottom192 (double x) const
const double	getIntegrand_sigmaWithISR_bottom196 (double x) const
const double	getIntegrand_sigmaWithISR_bottom200 (double x) const
const double	getIntegrand_sigmaWithISR_bottom202 (double x) const
const double	getIntegrand_sigmaWithISR_bottom205 (double x) const
const double	getIntegrand_sigmaWithISR_bottom207 (double x) const
const double	getIntegrand_sigmaWithISR_charm130 (double x) const
const double	getIntegrand_sigmaWithISR_charm133 (double x) const
const double	getIntegrand_sigmaWithISR_charm136 (double x) const
const double	getIntegrand_sigmaWithISR_charm161 (double x) const
const double	getIntegrand_sigmaWithISR_charm167 (double x) const
const double	getIntegrand_sigmaWithISR_charm172 (double x) const
const double	getIntegrand_sigmaWithISR_charm183 (double x) const
const double	getIntegrand_sigmaWithISR_charm189 (double x) const
const double	getIntegrand_sigmaWithISR_charm192 (double x) const
const double	getIntegrand_sigmaWithISR_charm196 (double x) const
const double	getIntegrand_sigmaWithISR_charm200 (double x) const
const double	getIntegrand_sigmaWithISR_charm202 (double x) const
const double	getIntegrand_sigmaWithISR_charm205 (double x) const
const double	getIntegrand_sigmaWithISR_charm207 (double x) const
const double	getIntegrand_sigmaWithISR_down130 (double x) const
const double	getIntegrand_sigmaWithISR_down133 (double x) const
const double	getIntegrand_sigmaWithISR_down136 (double x) const
const double	getIntegrand_sigmaWithISR_down161 (double x) const
const double	getIntegrand_sigmaWithISR_down167 (double x) const
const double	getIntegrand_sigmaWithISR_down172 (double x) const
const double	getIntegrand_sigmaWithISR_down183 (double x) const
const double	getIntegrand_sigmaWithISR_down189 (double x) const
const double	getIntegrand_sigmaWithISR_down192 (double x) const
const double	getIntegrand_sigmaWithISR_down196 (double x) const
const double	getIntegrand_sigmaWithISR_down200 (double x) const
const double	getIntegrand_sigmaWithISR_down202 (double x) const
const double	getIntegrand_sigmaWithISR_down205 (double x) const
const double	getIntegrand_sigmaWithISR_down207 (double x) const
const double	getIntegrand_sigmaWithISR_mu130 (double x) const
const double	getIntegrand_sigmaWithISR_mu136 (double x) const
const double	getIntegrand_sigmaWithISR_mu161 (double x) const
const double	getIntegrand_sigmaWithISR_mu172 (double x) const
const double	getIntegrand_sigmaWithISR_mu183 (double x) const
const double	getIntegrand_sigmaWithISR_mu189 (double x) const
const double	getIntegrand_sigmaWithISR_mu192 (double x) const
const double	getIntegrand_sigmaWithISR_mu196 (double x) const
const double	getIntegrand_sigmaWithISR_mu200 (double x) const
const double	getIntegrand_sigmaWithISR_mu202 (double x) const
const double	getIntegrand_sigmaWithISR_mu205 (double x) const
const double	getIntegrand_sigmaWithISR_mu207 (double x) const
const double	getIntegrand_sigmaWithISR_strange130 (double x) const
const double	getIntegrand_sigmaWithISR_strange133 (double x) const
const double	getIntegrand_sigmaWithISR_strange136 (double x) const
const double	getIntegrand_sigmaWithISR_strange161 (double x) const
const double	getIntegrand_sigmaWithISR_strange167 (double x) const
const double	getIntegrand_sigmaWithISR_strange172 (double x) const
const double	getIntegrand_sigmaWithISR_strange183 (double x) const
const double	getIntegrand_sigmaWithISR_strange189 (double x) const
const double	getIntegrand_sigmaWithISR_strange192 (double x) const
const double	getIntegrand_sigmaWithISR_strange196 (double x) const
const double	getIntegrand_sigmaWithISR_strange200 (double x) const
const double	getIntegrand_sigmaWithISR_strange202 (double x) const
const double	getIntegrand_sigmaWithISR_strange205 (double x) const
const double	getIntegrand_sigmaWithISR_strange207 (double x) const
const double	getIntegrand_sigmaWithISR_tau130 (double x) const
const double	getIntegrand_sigmaWithISR_tau136 (double x) const
const double	getIntegrand_sigmaWithISR_tau161 (double x) const
const double	getIntegrand_sigmaWithISR_tau172 (double x) const
const double	getIntegrand_sigmaWithISR_tau183 (double x) const
const double	getIntegrand_sigmaWithISR_tau189 (double x) const
const double	getIntegrand_sigmaWithISR_tau192 (double x) const
const double	getIntegrand_sigmaWithISR_tau196 (double x) const
const double	getIntegrand_sigmaWithISR_tau200 (double x) const
const double	getIntegrand_sigmaWithISR_tau202 (double x) const
const double	getIntegrand_sigmaWithISR_tau205 (double x) const
const double	getIntegrand_sigmaWithISR_tau207 (double x) const
const double	getIntegrand_sigmaWithISR_up130 (double x) const
const double	getIntegrand_sigmaWithISR_up133 (double x) const
const double	getIntegrand_sigmaWithISR_up136 (double x) const
const double	getIntegrand_sigmaWithISR_up161 (double x) const
const double	getIntegrand_sigmaWithISR_up167 (double x) const
const double	getIntegrand_sigmaWithISR_up172 (double x) const
const double	getIntegrand_sigmaWithISR_up183 (double x) const
const double	getIntegrand_sigmaWithISR_up189 (double x) const
const double	getIntegrand_sigmaWithISR_up192 (double x) const
const double	getIntegrand_sigmaWithISR_up196 (double x) const
const double	getIntegrand_sigmaWithISR_up200 (double x) const
const double	getIntegrand_sigmaWithISR_up202 (double x) const
const double	getIntegrand_sigmaWithISR_up205 (double x) const
const double	getIntegrand_sigmaWithISR_up207 (double x) const
const double	Integrand_AFBnumeratorWithISR_l (double x, const QCD::lepton l_flavor, const double s) const
const double	Integrand_AFBnumeratorWithISR_q (double x, const QCD::quark q_flavor, const double s) const
const double	Integrand_dsigmaBox_l (double cosTheta, const QCD::lepton l_flavor, const double s) const
const double	Integrand_dsigmaBox_q (double cosTheta, const QCD::quark q_flavor, const double s) const
const double	Integrand_sigmaWithISR_l (double x, const QCD::lepton l_flavor, const double s) const
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)\).
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}}\).
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}}\).
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}}\).
double	RVh () const
	The singlet vector corrections to the hadronic \(Z\)-boson width, denoted as \(R_V^h\).
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)\).
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.
virtual void	setParameter (const std::string name, const double &value)
	A method to set the value of a parameter of StandardModel.
const double	sigma_NoISR_l (const QCD::lepton l_flavor, const double s) const
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\).
Protected Member Functions inherited from QCD
const double	MassOfNf (int nf) const
	The Mbar mass of the heaviest quark in the theory with Nf active flavour.
Protected Attributes inherited from StandardModel
double	A
	The CKM parameter \(A\) in the Wolfenstein parameterization.
double	ale
	The fine-structure constant \(\alpha\).
double	alpha21
double	alpha31
double	AlsMz
	The strong coupling constant at the Z-boson mass, \(\alpha_s(M_Z)\).
bool	bSigmaForAFB
bool	bSigmaForR
double	dAl5hMz
	The five-flavour hadronic contribution to the electromagnetic coupling, \(\Delta\alpha_{\mathrm{had}}^{(5)}(M_Z^2)\). (Non-input parameter)
double	dAle5Mz
	The five-flavour hadronic contribution to the electromagnetic coupling, \(\Delta\alpha_{\mathrm{had}}^{(5)}(M_Z^2)\), used as input for FlagMWinput = FALSE.
double	delGammaWlv
	The theoretical uncertainty in \(\Gamma_W_{l\nu}\), denoted as \(\delta\,\Gamma_W_{l\nu}\).
double	delGammaWqq
	The theoretical uncertainty in \(\Gamma_W_{qq}\), denoted as \(\delta\,\Gamma_W_{qq}\).
double	delGammaZ
	The theoretical uncertainty in \(\Gamma_Z\), denoted as \(\delta\,\Gamma_Z\), in GeV.
double	delMw
	The theoretical uncertainty in \(M_W\), denoted as \(\delta\,M_W\), in GeV.
double	delR0b
	The theoretical uncertainty in \(R_b^0\), denoted as \(\delta\,R_b^0\).
double	delR0c
	The theoretical uncertainty in \(R_c^0\), denoted as \(\delta\,R_c^0\).
double	delR0l
	The theoretical uncertainty in \(R_l^0\), denoted as \(\delta\,R_l^0\).
double	delsigma0H
	The theoretical uncertainty in \(\sigma_{Hadron}^0\), denoted as \(\delta\,\sigma_{Hadron}^0\) in nb.
double	delSin2th_b
	The theoretical uncertainty in \(\sin^2\theta_{\rm eff}^{b}\), denoted as \(\delta\sin^2\theta_{\rm eff}^{b}\).
double	delSin2th_l
	The theoretical uncertainty in \(\sin^2\theta_{\rm eff}^{\rm lept}\), denoted as \(\delta\sin^2\theta_{\rm eff}^{\rm lept}\).
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}\).
double	delta
double	etab
	The CKM parameter \(\bar{\eta}\) in the Wolfenstein parameterization.
bool	flag_order [orders_EW_size]
	An array of internal flags controlling the inclusions of higher-order corrections.
bool	FlagFixMuwMut
	A boolean for the model flag FixMuwMut.
bool	flagLEP2 [NUMofLEP2RCs]
double	gamma
	\(\gamma \) used as an input for FlagWolfenstein = FALSE
double	GF
	The Fermi constant \(G_\mu\) in \({\rm GeV}^{-2}\).
double	lambda
	The CKM parameter \(\lambda\) in the Wolfenstein parameterization.
Particle	leptons [6]
	An array of Particle objects for the leptons.
double	mHl
	The Higgs mass \(m_h\) in GeV.
double	muw
	A matching scale \(\mu_W\) around the weak scale in GeV.
double	Mw_inp
	The mass of the \(W\) boson in GeV used as input for FlagMWinput = TRUE.
CKM	myCKM
	An object of type CKM.
PMNS	myPMNS
double	Mz
	The mass of the \(Z\) boson in GeV.
bool	requireCKM
	An internal flag to control whether the CKM matrix has to be recomputed.
bool	requireYe
	An internal flag to control whether the charged-lepton Yukawa matrix has to be recomputed.
bool	requireYn
	An internal flag to control whether the neutrino Yukawa matrix has to be recomputed.
double	rhob
	The CKM parameter \(\bar{\rho}\) in the Wolfenstein parameterization.
double	s12
double	s13
double	s23
Flavour	SMFlavour
	An object of type Flavour.
Matching< StandardModelMatching, StandardModel >	SMM
	An object of type Matching.
double	Vcb
	\(\vert V_{cb} \vert \) used as an input for FlagWolfenstein = FALSE
double	Vub
	\(\vert V_{ub} \vert \) used as an input for FlagWolfenstein = FALSE
double	Vud
	\(\vert V_{ud} \vert \) used as an input for FlagWolfenstein = FALSE and FlagUseVud = TRUE
double	Vus
	\(\vert V_{us} \vert \) used as an input for FlagWolfenstein = FALSE
gslpp::matrix< gslpp::complex >	Yd
	The Yukawa matrix of the down-type quarks.
gslpp::matrix< gslpp::complex >	Ye
	The Yukawa matrix of the charged leptons.
gslpp::matrix< gslpp::complex >	Yn
	The Yukawa matrix of the neutrinos.
gslpp::matrix< gslpp::complex >	Yu
	The Yukawa matrix of the up-type quarks.
Protected Attributes inherited from QCD
double	AlsM
	The strong coupling constant at the mass scale MAls, \(\alpha_s(M_{\alpha_s})\).
double	CA
double	CF
bool	computemt
	Switch for computing the \(\overline{\mathrm{MS}}\) mass of the top quark.
double	dAdA_NA
double	dFdA_NA
double	dFdF_NA
bool	FlagMpole2MbarNumeric
	A flag to determine whether the pole mass to \(\over \mathrm{MS}\) mass conversion is done numerically.
bool	FlagMtPole
	A flag to determine whether the pole mass of the top quark is used as input.
double	MAls
	The mass scale in GeV at which the strong coupling measurement is provided.
double	mtpole
	The pole mass of the top quark.
double	mub
	The threshold between five- and four-flavour theory in GeV.
double	muc
	The threshold between four- and three-flavour theory in GeV.
double	mut
	The threshold between six- and five-flavour theory in GeV.
double	NA
double	Nc
	The number of colours.
bool	QCDsuccess =true
Particle	quarks [6]
	The vector of all SM quarks.
bool	requireYd
	Switch for generating the Yukawa couplings to the down-type quarks.
bool	requireYu
	Switch for generating the Yukawa couplings to the up-type quarks.
double	TF
Protected Attributes inherited from Model
bool	isSliced = false
	A boolean set to true if the current istance is a slice of an extended object.
std::map< std::string, std::reference_wrapper< const double > >	ModelParamMap
bool	UpdateError = false
	A boolean set to false if update is successful.

Constructor & Destructor Documentation

EWSM_Output()

EWSM_Output::EWSM_Output

(

const StandardModel &

SM_in

)

Constructor.

Parameters

[in]

EWSM_in

a reference to an object of type EWSM

Definition at line 20 of file EWSM_Output.cpp.

: StandardModel(SM_in)
{
}

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.

{
    /* rescale */
    double DeltaRho_new[StandardModel::orders_EW_size];
    for (int j = 0; j < StandardModel::orders_EW_size; ++j)
        DeltaRho_new[j] = cW2overSW2 * DeltaRho[j];
 
    if (getFlagKappaZ().compare("APPROXIMATEFORMULA") == 0) {
        std::cout << "Delta kappaZb (from the approximate formula of sin2thb) = "
                << kappaZ_f(getQuarks(QCD::BOTTOM)) - 1.0 << std::endl;
    } else if (getFlagKappaZ().compare("NORESUM") == 0) {
        std::cout << "Leading contributions: alpha or Gmu" << std::endl;
        std::cout << "  DeltaRho[EW1]=" << DeltaRho_new[StandardModel::EW1] << " "
                << f_AlphaToGF * DeltaRho_new[StandardModel::EW1] << std::endl;
        std::cout << "  DeltaRho[EW1QCD1]=" << DeltaRho_new[StandardModel::EW1QCD1] << " "
                << f_AlphaToGF * DeltaRho_new[StandardModel::EW1QCD1] << std::endl;
        std::cout << "  DeltaRho[EW1QCD2]=" << DeltaRho_new[StandardModel::EW1QCD2] << " "
                << f_AlphaToGF * DeltaRho_new[StandardModel::EW1QCD2] << std::endl;
        std::cout << "  DeltaRho[EW2]=" << DeltaRho_new[StandardModel::EW2] << " "
                << pow(f_AlphaToGF, 2.0) * DeltaRho_new[StandardModel::EW2] << std::endl;
        std::cout << "  DeltaRho[EW2QCD1]=" << DeltaRho_new[StandardModel::EW2QCD1] << " "
                << pow(f_AlphaToGF, 2.0) * DeltaRho_new[StandardModel::EW2QCD1] << std::endl;
        std::cout << "  DeltaRho[EW3]=" << DeltaRho_new[StandardModel::EW3] << " "
                << pow(f_AlphaToGF, 3.0) * DeltaRho_new[StandardModel::EW3] << std::endl;
        std::cout << "EW2 from Z-gamma = " << Zgamma_EW2 << std::endl;
        std::cout << "Remainder contributions: alpha or Gmu" << std::endl;
        std::cout << "  DeltaRbar_rem[EW1]=" << DeltaRbar_rem << " "
                << f_AlphaToGF * DeltaRbar_rem << std::endl;
        std::cout << "  deltaKappa_rem[EW1]=" << deltaKappa_rem[StandardModel::EW1] - ZbbSubtract << " "
                << f_AlphaToGF * (deltaKappa_rem[StandardModel::EW1] - ZbbSubtract) << std::endl;
        std::cout << "  deltaKappa_rem[EW1QCD1]=" << deltaKappa_rem[StandardModel::EW1QCD1] << " "
                << f_AlphaToGF * deltaKappa_rem[StandardModel::EW1QCD1] << std::endl;
        std::cout << "  deltaKappa_rem[EW1QCD2]=" << deltaKappa_rem[StandardModel::EW1QCD2] << " "
                << f_AlphaToGF * deltaKappa_rem[StandardModel::EW1QCD2] << std::endl;
        std::cout << "  deltaKappa_rem[EW2]=" << deltaKappa_rem[StandardModel::EW2] << " "
                << pow(f_AlphaToGF, 2.0) * deltaKappa_rem[StandardModel::EW2] << std::endl;
        if (bool_Zbb) {
            std::cout << "Taub: alpha or Gmu" << std::endl;
            std::cout << "  taub[EW1]=" << taub[StandardModel::EW1] << " "
                    << f_AlphaToGF * taub[StandardModel::EW1] << std::endl;
            std::cout << "  taub[EW1QCD1]=" << taub[StandardModel::EW1QCD1] << " "
                    << f_AlphaToGF * taub[StandardModel::EW1QCD1] << std::endl;
            std::cout << "  taub[EW2]=" << taub[StandardModel::EW2] << " "
                    << pow(f_AlphaToGF, 2.0) * taub[StandardModel::EW2] << std::endl;
        }
        std::cout << "Each order: alpha or Gmu" << std::endl;
        double dKappa_EW1 = DeltaRho_new[StandardModel::EW1] + deltaKappa_rem[StandardModel::EW1] - ZbbSubtract;
        double dKappa_EW1QCD1 = DeltaRho_new[StandardModel::EW1QCD1] + deltaKappa_rem[StandardModel::EW1QCD1];
        double dKappa_EW2 = (deltaKappa_rem[StandardModel::EW1] - ZbbSubtract) * DeltaRho_new[StandardModel::EW1]
                - DeltaRho_new[StandardModel::EW1] * DeltaRbar_rem
                + DeltaRho_new[StandardModel::EW2]
                + deltaKappa_rem[StandardModel::EW2];
        double dKappa_EW1QCD2 = DeltaRho_new[StandardModel::EW1QCD2] + deltaKappa_rem[StandardModel::EW1QCD2];
        double dKappa_EW2QCD1 = DeltaRho_new[StandardModel::EW2QCD1] + deltaKappa_rem[StandardModel::EW1QCD1] * DeltaRho_new[StandardModel::EW1];
        double dKappa_EW3 = DeltaRho_new[StandardModel::EW3];
        double dKappa_EW2QCD2 = deltaKappa_rem[StandardModel::EW1QCD2] * DeltaRho_new[StandardModel::EW1];
        //
        double dKappa_EW1_TMP = dKappa_EW1;
        double dKappa_EW1QCD1_TMP = dKappa_EW1QCD1;
        double dKappa_EW2_TMP = dKappa_EW2;
        double dKappa_EW1QCD2_TMP = dKappa_EW1QCD2;
        double dKappa_EW2QCD1_TMP = dKappa_EW2QCD1;
        double dKappa_EW3_TMP = dKappa_EW3;
        double dKappa_EW2QCD2_TMP = dKappa_EW2QCD2;
        //
        if (bool_Zbb) {
            dKappa_EW1 = dKappa_EW1_TMP - taub[StandardModel::EW1];
            dKappa_EW1QCD1 = dKappa_EW1QCD1_TMP - taub[StandardModel::EW1QCD1];
            dKappa_EW2 = dKappa_EW2_TMP - taub[StandardModel::EW2] + taub[StandardModel::EW1] * taub[StandardModel::EW1]
                    - dKappa_EW1_TMP * taub[StandardModel::EW1];
            dKappa_EW1QCD2 = dKappa_EW1QCD2_TMP;
            dKappa_EW2QCD1 = dKappa_EW2QCD1_TMP - dKappa_EW1_TMP * taub[StandardModel::EW1QCD1]
                    - dKappa_EW1QCD1_TMP * taub[StandardModel::EW1] + 2.0 * taub[StandardModel::EW1] * taub[StandardModel::EW1QCD1];
            dKappa_EW3 = dKappa_EW3_TMP - dKappa_EW1_TMP * taub[StandardModel::EW2]
                    - dKappa_EW2_TMP * taub[StandardModel::EW1] + 2.0 * taub[StandardModel::EW2] * taub[StandardModel::EW1];
            dKappa_EW2QCD2 = dKappa_EW2QCD2_TMP - dKappa_EW1QCD2_TMP * taub[StandardModel::EW1]
                    - dKappa_EW1QCD1_TMP * taub[StandardModel::EW1QCD1]
                    + taub[StandardModel::EW1QCD1] * taub[StandardModel::EW1QCD1];
        }
        std::cout << "  EW1: " << dKappa_EW1 << " " << f_AlphaToGF * dKappa_EW1 << std::endl;
        std::cout << "  EW1QCD1: " << dKappa_EW1QCD1 << " " << f_AlphaToGF * dKappa_EW1QCD1 << std::endl;
        std::cout << "  EW2: " << dKappa_EW2 + Zgamma_EW2
                << " " << pow(f_AlphaToGF, 2.0) * dKappa_EW2 + Zgamma_EW2 << std::endl;
        std::cout << "  EW1QCD2: " << dKappa_EW1QCD2 << " " << f_AlphaToGF * dKappa_EW1QCD2 << std::endl;
        std::cout << "  EW2QCD1: " << dKappa_EW2QCD1 << " " << pow(f_AlphaToGF, 2.0) * dKappa_EW2QCD1 << std::endl;
        std::cout << "  EW3: " << dKappa_EW3 << " " << pow(f_AlphaToGF, 3.0) * dKappa_EW3 << std::endl;
        std::cout << "  EW2QCD2: " << dKappa_EW2QCD2 << " " << pow(f_AlphaToGF, 2.0) * dKappa_EW2QCD2 << std::endl;
        std::cout << "Total contribution: alpha or Gmu" << std::endl;
        std::cout << "  kappaZ="
                << 1.0 + dKappa_EW1 + dKappa_EW1QCD1 + dKappa_EW2
                + dKappa_EW1QCD2 + dKappa_EW2QCD1 + dKappa_EW3
                + dKappa_EW2QCD2 + Zgamma_EW2
                << " "
                << 1.0 + f_AlphaToGF * dKappa_EW1 + f_AlphaToGF * dKappa_EW1QCD1
                + pow(f_AlphaToGF, 2.0) * dKappa_EW2 + f_AlphaToGF * dKappa_EW1QCD2
                + pow(f_AlphaToGF, 2.0) * dKappa_EW2QCD1
                + pow(f_AlphaToGF, 3.0) * dKappa_EW3
                + pow(f_AlphaToGF, 2.0) * dKappa_EW2QCD2
                + Zgamma_EW2 << std::endl;
        if (bool_Zbb) {
            std::cout << "  kappaZ(taub resummed)="
                    << (1.0 + dKappa_EW1_TMP + dKappa_EW1QCD1_TMP
                    + dKappa_EW2_TMP + dKappa_EW1QCD2_TMP
                    + dKappa_EW2QCD1_TMP + dKappa_EW3_TMP
                    + dKappa_EW2QCD2_TMP)
                    / (1.0 + taub[StandardModel::EW1] + taub[StandardModel::EW1QCD1] + taub[StandardModel::EW2])
                    + Zgamma_EW2
                    << " "
                    << (1.0 + f_AlphaToGF * dKappa_EW1_TMP
                    + f_AlphaToGF * dKappa_EW1QCD1_TMP
                    + pow(f_AlphaToGF, 2.0) * dKappa_EW2_TMP
                    + f_AlphaToGF * dKappa_EW1QCD2_TMP
                    + pow(f_AlphaToGF, 2.0) * dKappa_EW2QCD1_TMP
                    + pow(f_AlphaToGF, 3.0) * dKappa_EW3_TMP
                    + pow(f_AlphaToGF, 2.0) * dKappa_EW2QCD2_TMP)
                    / (1.0 + f_AlphaToGF * taub[StandardModel::EW1]
                    + f_AlphaToGF * taub[StandardModel::EW1QCD1]
                    + pow(f_AlphaToGF, 2.0) * taub[StandardModel::EW2])
                    + Zgamma_EW2
                    << std::endl;
        }
    } else
        std::cout << "EWSM_Output::outputEachDeltaKappaZ(): Not implemented for schemeKappaZ="
            << getFlagKappaZ() << std::endl;
}

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.

{
    std::cout << "================================================" << std::endl;
    std::cout << "kappaZ_l[(QCD::lepton)" << l << "]" << std::endl;
    std::cout << "Mw(input)   = " << Mw_i << std::endl;
 
    double cW2_TMP = Mw_i * Mw_i / getMz() / getMz();
    double sW2_TMP = 1.0 - cW2_TMP;
 
    double DeltaRho[StandardModel::orders_EW_size];
    DeltaRho[StandardModel::EW1] = getMyOneLoopEW()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW1QCD1] = getMyTwoLoopQCD()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW1QCD2] = getMyThreeLoopQCD()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW2] = getMyTwoLoopEW()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW2QCD1] = getMyThreeLoopEW2QCD()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW3] = getMyThreeLoopEW()->DeltaRho(Mw_i);
 
    /* compute delta kappa_rem^f */
    gslpp::complex deltaKappa_rem_f[StandardModel::orders_EW_size];
    deltaKappa_rem_f[StandardModel::EW1] = getMyOneLoopEW()->deltaKappa_rem_f(getLeptons(l), Mw_i);
#ifdef WITHIMTWOLOOPQCD
    deltaKappa_rem_f[StandardModel::EW1QCD1] = gslpp::complex(getMyTwoLoopQCD()->deltaKappa_rem_f(getLeptons(l), Mw_i).real(),
            getMyTwoLoopQCD()->deltaKappa_rem_f(getLeptons(l), Mw_i).imag(), false);
#else
    deltaKappa_rem_f[StandardModel::EW1QCD1] = gslpp::complex(getMyTwoLoopQCD()->deltaKappa_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
#endif
    deltaKappa_rem_f[StandardModel::EW1QCD2] = gslpp::complex(getMyThreeLoopQCD()->deltaKappa_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
    deltaKappa_rem_f[StandardModel::EW2] = gslpp::complex(getMyTwoLoopEW()->deltaKappa_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
    deltaKappa_rem_f[StandardModel::EW2QCD1] = gslpp::complex(getMyThreeLoopEW2QCD()->deltaKappa_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
    deltaKappa_rem_f[StandardModel::EW3] = gslpp::complex(getMyThreeLoopEW()->deltaKappa_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
 
    /* compute Delta rbar_rem */
    double DeltaRbar_rem = getMyOneLoopEW()->DeltaRbar_rem(Mw_i);
 
    /* conversion factor */
    double f_AlphaToGF = sqrt(2.0) * getGF() * pow(getMz(), 2.0)
            * sW2_TMP * cW2_TMP / M_PI / getAle();
 
    /* Re[Kappa_Z^f] with or without resummation */
    double deltaKappa_rem_f_real[StandardModel::orders_EW_size];
    for (int j = 0; j < StandardModel::orders_EW_size; ++j)
        deltaKappa_rem_f_real[j] = deltaKappa_rem_f[j].real();
 
    /* O(alpha^2) correction to Re[kappa_Z^f] from the Z-gamma mixing */
    double ReKappaZf = resumKappaZ(DeltaRho, deltaKappa_rem_f_real,
            DeltaRbar_rem, false);
    double Zgamma_EW2 = 35.0 * alphaMz() * alphaMz() / 18.0 / sW2_TMP
            * (1.0 - 8.0 / 3.0 * ReKappaZf * sW2_TMP);
 
    double dummy[StandardModel::orders_EW_size];
    outputEachDeltaKappaZ(f_AlphaToGF, cW2_TMP / sW2_TMP,
            DeltaRho, deltaKappa_rem_f_real,
            DeltaRbar_rem, false, dummy, 0.0, Zgamma_EW2);
 
    /* Im[kappa_Z^f] without resummation */
    double ImKappaZf = 0.0;
    for (int j = 0; j < StandardModel::orders_EW_size; ++j)
        ImKappaZf += deltaKappa_rem_f[j].imag();
    std::cout << "ImKappaZf(with alpha)=" << ImKappaZf << std::endl;
 
    std::cout << "================================================" << std::endl;
}

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.

{
    std::cout << "================================================" << std::endl;
    std::cout << "kappaZ_q[(QCD::quark)" << q << "]" << std::endl;
    std::cout << "Mw(input)   = " << Mw_i << std::endl;
 
    double cW2_TMP = Mw_i * Mw_i / getMz() / getMz();
    double sW2_TMP = 1.0 - cW2_TMP;
 
    double DeltaRho[StandardModel::orders_EW_size];
    DeltaRho[StandardModel::EW1] = getMyOneLoopEW()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW1QCD1] = getMyTwoLoopQCD()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW1QCD2] = getMyThreeLoopQCD()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW2] = getMyTwoLoopEW()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW2QCD1] = getMyThreeLoopEW2QCD()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW3] = getMyThreeLoopEW()->DeltaRho(Mw_i);
 
    /* compute delta kappa_rem^f */
    gslpp::complex deltaKappa_rem_f[StandardModel::orders_EW_size];
    deltaKappa_rem_f[StandardModel::EW1] = getMyOneLoopEW()->deltaKappa_rem_f(getQuarks(q), Mw_i);
#ifdef WITHIMTWOLOOPQCD
    deltaKappa_rem_f[StandardModel::EW1QCD1] = gslpp::complex(getMyTwoLoopQCD()->deltaKappa_rem_f(getQuarks(q), Mw_i).real(),
            getMyTwoLoopQCD()->deltaKappa_rem_f(getQuarks(q), Mw_i).imag(), false);
#else
    deltaKappa_rem_f[StandardModel::EW1QCD1] = gslpp::complex(getMyTwoLoopQCD()->deltaKappa_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
#endif
    deltaKappa_rem_f[StandardModel::EW1QCD2] = gslpp::complex(getMyThreeLoopQCD()->deltaKappa_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
    deltaKappa_rem_f[StandardModel::EW2] = gslpp::complex(getMyTwoLoopEW()->deltaKappa_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
    deltaKappa_rem_f[StandardModel::EW2QCD1] = gslpp::complex(getMyThreeLoopEW2QCD()->deltaKappa_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
    deltaKappa_rem_f[StandardModel::EW3] = gslpp::complex(getMyThreeLoopEW()->deltaKappa_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
 
    /* compute Delta rbar_rem */
    double DeltaRbar_rem = getMyOneLoopEW()->DeltaRbar_rem(Mw_i);
 
    /* conversion factor */
    double f_AlphaToGF = sqrt(2.0) * getGF() * pow(getMz(), 2.0)
            * sW2_TMP * cW2_TMP / M_PI / getAle();
 
    /* Zbb */
    bool bool_Zbb = false;
    if (q == QCD::BOTTOM) bool_Zbb = true;
    double ZbbSubtract = 0.0;
    if (bool_Zbb)
        ZbbSubtract = getAle() / 8.0 / M_PI / sW2_TMP
            * pow(getMtpole() / Mw_i, 2.0);
    double taub[StandardModel::orders_EW_size];
    double Xt = getMyEWSMcache()->Xt_alpha(Mw_i);
    if (bool_Zbb) {
        taub[StandardModel::EW1] = -2.0 * Xt;
        taub[StandardModel::EW1QCD1] = 2.0 / 3.0 * M_PI * Xt * getMyEWSMcache()->alsMt();
        taub[StandardModel::EW2] = -2.0 * Xt * Xt * getMyTwoLoopEW()->tau_2();
    }
 
    /* Re[Kappa_Z^f] with or without resummation */
    double deltaKappa_rem_f_real[StandardModel::orders_EW_size];
    for (int j = 0; j < StandardModel::orders_EW_size; ++j)
        deltaKappa_rem_f_real[j] = deltaKappa_rem_f[j].real();
 
    /* O(alpha^2) correction to Re[kappa_Z^f] from the Z-gamma mixing */
    double ReKappaZf = resumKappaZ(DeltaRho, deltaKappa_rem_f_real,
            DeltaRbar_rem, bool_Zbb);
    double Zgamma_EW2 = 35.0 * alphaMz() * alphaMz() / 18.0 / sW2_TMP
            * (1.0 - 8.0 / 3.0 * ReKappaZf * sW2_TMP);
 
    outputEachDeltaKappaZ(f_AlphaToGF, cW2_TMP / sW2_TMP,
            DeltaRho, deltaKappa_rem_f_real,
            DeltaRbar_rem, bool_Zbb, taub, ZbbSubtract, Zgamma_EW2);
 
    /* Im[kappa_Z^f] without resummation */
    double ImKappaZf = 0.0;
    for (int j = 0; j < StandardModel::orders_EW_size; ++j)
        ImKappaZf += deltaKappa_rem_f[j].imag();
    std::cout << "ImKappaZf(with alpha)=" << ImKappaZf << std::endl;
 
    std::cout << "================================================" << std::endl;
}

outputEachDeltaR()

void EWSM_Output::outputEachDeltaR

(

const double

Mw_i

)

const

Definition at line 25 of file EWSM_Output.cpp.

{
    std::cout << "Mw_SM          = " << Mw() << std::endl;
    std::cout << "DeltaR_SM()    = " << DeltaR() << std::endl;
    std::cout << "DeltaRbar_SM() = " << DeltaRbar() << std::endl;
    std::cout << "Mw(input)      = " << Mw_i << std::endl;
 
    double cW2_TMP = Mw_i * Mw_i / getMz() / getMz();
    double sW2_TMP = 1.0 - cW2_TMP;
 
    double DeltaRho[StandardModel::orders_EW_size];
    DeltaRho[StandardModel::EW1] = getMyOneLoopEW()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW1QCD1] = getMyTwoLoopQCD()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW1QCD2] = getMyThreeLoopQCD()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW2] = getMyTwoLoopEW()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW2QCD1] = getMyThreeLoopEW2QCD()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW3] = getMyThreeLoopEW()->DeltaRho(Mw_i);
 
    double DeltaR_rem[StandardModel::orders_EW_size];
    DeltaR_rem[StandardModel::EW1] = getMyOneLoopEW()->DeltaR_rem(Mw_i);
    DeltaR_rem[StandardModel::EW1QCD1] = getMyTwoLoopQCD()->DeltaR_rem(Mw_i);
    DeltaR_rem[StandardModel::EW1QCD2] = getMyThreeLoopQCD()->DeltaR_rem(Mw_i);
    DeltaR_rem[StandardModel::EW2] = getMyTwoLoopEW()->DeltaR_rem(Mw_i);
    DeltaR_rem[StandardModel::EW2QCD1] = getMyThreeLoopEW2QCD()->DeltaR_rem(Mw_i);
    DeltaR_rem[StandardModel::EW3] = getMyThreeLoopEW()->DeltaR_rem(Mw_i);
 
    double f_AlphaToGF = sqrt(2.0) * getGF() * pow(getMz(), 2.0) * sW2_TMP * cW2_TMP / M_PI / getAle();
    //f_AlphaToGF = 1.0; /* for test */
    double DeltaRho_sum = f_AlphaToGF * DeltaRho[StandardModel::EW1]
            + f_AlphaToGF * DeltaRho[StandardModel::EW1QCD1]
            + f_AlphaToGF * DeltaRho[StandardModel::EW1QCD2]
            + pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2]
            + pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2QCD1]
            + pow(f_AlphaToGF, 3.0) * DeltaRho[StandardModel::EW3];
    double DeltaRho_G = f_AlphaToGF * DeltaRho[StandardModel::EW1];
 
    if (getFlagMw().compare("NORESUM") == 0) {
 
        f_AlphaToGF = 1.0;
        DeltaRho[StandardModel::EW1QCD2] *= f_AlphaToGF;
        DeltaRho[StandardModel::EW2QCD1] *= pow(f_AlphaToGF, 2.0);
        DeltaRho[StandardModel::EW3] *= pow(f_AlphaToGF, 3.0);
 
        // Full EW one-loop contribution (without the full DeltaAlphaL5q)
        double DeltaR_EW1 = -cW2_TMP / sW2_TMP * DeltaRho[StandardModel::EW1] + DeltaR_rem[StandardModel::EW1];
 
        // Full EW two-loop contribution with reducible corrections
        double DeltaR_EW2_rem = getMyApproximateFormulae()->DeltaR_TwoLoopEW_rem(Mw_i);
 
        // EW two-loop irreducible contributions with large-mt expansion
        double DeltaR_EW2_old_red = DeltaAlphaL5q() * DeltaAlphaL5q()
                - 2.0 * cW2_TMP / sW2_TMP * DeltaAlphaL5q() * DeltaRho[StandardModel::EW1]
                + pow(cW2_TMP / sW2_TMP * DeltaRho[StandardModel::EW1], 2.0);
        double DeltaR_EW2_old_irred = -cW2_TMP / sW2_TMP * DeltaRho[StandardModel::EW2] + DeltaR_rem[StandardModel::EW2];
 
        // Delta r, including the full EW two-loop contribution
        double deltaR = DeltaAlphaL5q();
        for (int j = 0; j < StandardModel::orders_EW_size; ++j) {
            deltaR += -cW2_TMP / sW2_TMP * DeltaRho[(StandardModel::orders_EW)j];
            deltaR += DeltaR_rem[(StandardModel::orders_EW)j];
        }
        deltaR -= -cW2_TMP / sW2_TMP * DeltaRho[StandardModel::EW2];
        deltaR -= DeltaR_rem[StandardModel::EW2];
        deltaR += DeltaAlphaL5q() * DeltaAlphaL5q() + 2.0 * DeltaAlphaL5q() * DeltaR_EW1 + DeltaR_EW2_rem;
 
        std::cout << "(1+dr) - 1        =  " << deltaR << std::endl;
        std::cout << "  EW1             =  " << DeltaAlphaL5q() + DeltaR_EW1 << std::endl;
        std::cout << "    DeltaAlphaL5q =  " << DeltaAlphaL5q() << std::endl;
        std::cout << "    dR            = " << DeltaR_EW1 << std::endl;
        std::cout << "  EW1QCD1         =  " << -cW2_TMP / sW2_TMP * DeltaRho[StandardModel::EW1QCD1] + DeltaR_rem[StandardModel::EW1QCD1] << std::endl;
        std::cout << "  EW2(full)       =  " << DeltaR_EW2_rem + DeltaAlphaL5q() * DeltaAlphaL5q() + 2.0 * DeltaAlphaL5q() * DeltaR_EW1 << std::endl;
        std::cout << "    dAle*dAle     =  " << DeltaAlphaL5q() * DeltaAlphaL5q() << std::endl;
        std::cout << "    2*dAle*dR     = " << 2.0 * DeltaAlphaL5q() * DeltaR_EW1 << std::endl;
        std::cout << "    others        =  " << DeltaR_EW2_rem << std::endl;
        std::cout << "  EW1QCD2         =  " << -cW2_TMP / sW2_TMP * DeltaRho[StandardModel::EW1QCD2] + DeltaR_rem[StandardModel::EW1QCD2] << std::endl;
        std::cout << "  EW2QCD1         = " << -cW2_TMP / sW2_TMP * DeltaRho[StandardModel::EW2QCD1] + DeltaR_rem[StandardModel::EW2QCD1] << std::endl;
        std::cout << "  EW3             = " << -cW2_TMP / sW2_TMP * DeltaRho[StandardModel::EW3] + DeltaR_rem[StandardModel::EW3] << std::endl;
        std::cout << "  EW2(old,irreducible) = " << DeltaR_EW2_old_irred << std::endl;
        std::cout << "  EW2(old,red+irred)   = " << DeltaR_EW2_old_red + DeltaR_EW2_old_irred << std::endl;
        std::cout << "  EW2(old,red+irred-dAle*dAle-2*dAle*dR) = "
                << DeltaR_EW2_old_red + DeltaR_EW2_old_irred
                - DeltaAlphaL5q() * DeltaAlphaL5q()
                - 2.0 * DeltaAlphaL5q() * DeltaR_EW1 << std::endl;
 
    } else if (getFlagMw().compare("OMSI") == 0) {
 
        // R = 1/(1 - Delta r)
        double R = 1.0 / (1.0 + cW2_TMP / sW2_TMP * DeltaRho_sum)
                / (1.0 - DeltaAlphaL5q() - DeltaR_rem[StandardModel::EW1] - DeltaR_rem[StandardModel::EW1QCD1] - DeltaR_rem[StandardModel::EW2]);
 
        std::cout << "1/(1-dr) - 1 (exact)                 = " << R - 1.0 << std::endl;
        std::cout << "     --> dr = " << 1.0 - 1.0 / R << std::endl;
 
        // each contribution
        double DeltaR_EW1 = DeltaAlphaL5q() - cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1] + DeltaR_rem[StandardModel::EW1];
        double DeltaR_EW1QCD1 = -cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1QCD1] + DeltaR_rem[StandardModel::EW1QCD1];
        double DeltaR_EW2 = -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2]
                + DeltaR_rem[StandardModel::EW2]
                + cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1]
                *(DeltaAlphaL5q() + DeltaR_rem[StandardModel::EW1])
                + DeltaR_EW1*DeltaR_EW1;
        double DeltaR_EW1QCD2 = -cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1QCD2];
        double DeltaR_EW2QCD1 = -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2QCD1]
                + cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1QCD1]
                *(DeltaAlphaL5q() + DeltaR_rem[StandardModel::EW1])
                + 2.0 * DeltaR_EW1*DeltaR_EW1QCD1;
        double DeltaR_EW3 = -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 3.0) * DeltaRho[StandardModel::EW3]
                + cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2]
                *(DeltaAlphaL5q() + DeltaR_rem[StandardModel::EW1])
                + pow(DeltaR_EW1, 3.0)
                + 2.0 * DeltaR_EW1 * (DeltaR_EW2 - DeltaR_EW1 * DeltaR_EW1);
 
        std::cout << "  EW1             =  " << DeltaR_EW1 << std::endl;
        std::cout << "    DeltaAlphaL5q =  " << DeltaAlphaL5q() << std::endl;
        std::cout << "    -cW2/sW2*dRho1= " << -cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1] << std::endl;
        std::cout << "    DeltaR1_rem   =  " << DeltaR_rem[StandardModel::EW1] << std::endl;
        std::cout << "  EW1QCD1         =  " << DeltaR_EW1QCD1 << std::endl;
        std::cout << "  EW2(full)       =  " << DeltaR_EW2 << std::endl;
        std::cout << "    EW1*EW1       =  " << DeltaR_EW1 * DeltaR_EW1 << std::endl;
        std::cout << "      dAle*dAle   =  " << DeltaAlphaL5q() * DeltaAlphaL5q() << std::endl;
        std::cout << "      others      = " << DeltaR_EW1 * DeltaR_EW1 - DeltaAlphaL5q() * DeltaAlphaL5q() << std::endl;
        std::cout << "    -cW2/sW2*dRho2=  " << -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2] << std::endl;
        std::cout << "    DeltaR2_rem   =  " << DeltaR_rem[StandardModel::EW2] << std::endl;
        std::cout << "    others        =  " << cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1]*(DeltaAlphaL5q() + DeltaR_rem[StandardModel::EW1]) << std::endl;
        std::cout << "  EW1QCD2         =  " << DeltaR_EW1QCD2 << std::endl;
        std::cout << "  EW2QCD1         = " << DeltaR_EW2QCD1 << std::endl;
        std::cout << "  EW3             =  " << DeltaR_EW3 << std::endl;
        std::cout << "    -cW2/sW2*dRho3= " << -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 3.0) * DeltaRho[StandardModel::EW3] << std::endl;
        std::cout << "    EW1^3         =  " << pow(DeltaR_EW1, 3.0) << std::endl;
        std::cout << "    2*EW1*(EW2-EW1^2)=" << 2.0 * DeltaR_EW1 * (DeltaR_EW2 - DeltaR_EW1 * DeltaR_EW1) << std::endl;
        std::cout << "    others        = " << cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2]*(DeltaAlphaL5q() + DeltaR_rem[StandardModel::EW1]) << std::endl;
 
    } else if (getFlagMw().compare("OMSII") == 0) {
 
        // R = 1/(1 - Delta r)
        double R = 1.0 / ((1.0 + cW2_TMP / sW2_TMP * DeltaRho_sum)*(1.0 - DeltaAlphaL5q())
                - (1.0 + cW2_TMP / sW2_TMP * DeltaRho_G) * DeltaR_rem[StandardModel::EW1]
                - DeltaR_rem[StandardModel::EW1QCD1] - DeltaR_rem[StandardModel::EW2]);
 
        // each contribution
        double DeltaR_EW1 = DeltaAlphaL5q() - cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1] + DeltaR_rem[StandardModel::EW1];
        double DeltaR_EW1QCD1 = -cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1QCD1] + DeltaR_rem[StandardModel::EW1QCD1];
        double DeltaR_EW2 = -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2]
                + DeltaR_rem[StandardModel::EW2]
                + cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1]
                *(DeltaAlphaL5q() + DeltaR_rem[StandardModel::EW1])
                + DeltaR_EW1*DeltaR_EW1;
        double DeltaR_EW1QCD2 = -cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1QCD2];
        double DeltaR_EW2QCD1 = -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2QCD1]
                + cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1QCD1] * DeltaAlphaL5q()
                + 2.0 * DeltaR_EW1*DeltaR_EW1QCD1;
        double DeltaR_EW3 = -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 3.0) * DeltaRho[StandardModel::EW3]
                + cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2] * DeltaAlphaL5q()
                + pow(DeltaR_EW1, 3.0)
                + 2.0 * DeltaR_EW1 * (DeltaR_EW2 - DeltaR_EW1 * DeltaR_EW1);
 
        std::cout << "1/(1-dr) - 1 (exact)                 = " << R - 1.0 << std::endl;
        std::cout << "     --> dr = " << 1.0 - 1.0 / R << std::endl;
        std::cout << "1/(1-dr) - 1 (sum of expanded terms) = "
                << DeltaR_EW1 + DeltaR_EW1QCD1 + DeltaR_EW2 + DeltaR_EW1QCD2
                + DeltaR_EW2QCD1 + DeltaR_EW3 << std::endl;
        std::cout << "  EW1             =  " << DeltaR_EW1 << std::endl;
        std::cout << "    DeltaAlphaL5q =  " << DeltaAlphaL5q() << std::endl;
        std::cout << "    -cW2/sW2*dRho1= " << -cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1] << std::endl;
        std::cout << "    DeltaR1_rem   =  " << DeltaR_rem[StandardModel::EW1] << std::endl;
        std::cout << "  EW1QCD1         =  " << DeltaR_EW1QCD1 << std::endl;
        std::cout << "  EW2(full)       =  " << DeltaR_EW2 << std::endl;
        std::cout << "    EW1*EW1       =  " << DeltaR_EW1 * DeltaR_EW1 << std::endl;
        std::cout << "      dAle*dAle   =  " << DeltaAlphaL5q() * DeltaAlphaL5q() << std::endl;
        std::cout << "      others      = " << DeltaR_EW1 * DeltaR_EW1 - DeltaAlphaL5q() * DeltaAlphaL5q() << std::endl;
        std::cout << "    -cW2/sW2*dRho2=  " << -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2] << std::endl;
        std::cout << "    DeltaR2_rem   =  " << DeltaR_rem[StandardModel::EW2] << std::endl;
        std::cout << "    others        =  " << cW2_TMP / sW2_TMP * f_AlphaToGF * DeltaRho[StandardModel::EW1]*(DeltaAlphaL5q() + DeltaR_rem[StandardModel::EW1]) << std::endl;
        std::cout << "  EW1QCD2         =  " << DeltaR_EW1QCD2 << std::endl;
        std::cout << "  EW2QCD1         = " << DeltaR_EW2QCD1 << std::endl;
        std::cout << "  EW3             =  " << DeltaR_EW3 << std::endl;
        std::cout << "    -cW2/sW2*dRho3= " << -cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 3.0) * DeltaRho[StandardModel::EW3] << std::endl;
        std::cout << "    EW1^3         =  " << pow(DeltaR_EW1, 3.0) << std::endl;
        std::cout << "    2*EW1*(EW2-EW1^2)=" << 2.0 * DeltaR_EW1 * (DeltaR_EW2 - DeltaR_EW1 * DeltaR_EW1) << std::endl;
        std::cout << "    others        = " << cW2_TMP / sW2_TMP * pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2] * DeltaAlphaL5q() << std::endl;
 
    } else
        std::cout << "EWSM_Output::outputEachDeltaR(): Not implemented for schemeMw="
            << getFlagMw() << std::endl;
}

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.

{
    if (getFlagRhoZ().compare("APPROXIMATEFORMULA") == 0) {
 
    } else if (getFlagRhoZ().compare("NORESUM") == 0) {
        std::cout << "Leading contributions: alpha or Gmu" << std::endl;
        std::cout << "  DeltaRho[EW1]=" << DeltaRho[StandardModel::EW1] << " "
                << f_AlphaToGF * DeltaRho[StandardModel::EW1] << std::endl;
        std::cout << "  DeltaRho[EW1QCD1]=" << DeltaRho[StandardModel::EW1QCD1] << " "
                << f_AlphaToGF * DeltaRho[StandardModel::EW1QCD1] << std::endl;
        std::cout << "  DeltaRho[EW1QCD2]=" << DeltaRho[StandardModel::EW1QCD2] << " "
                << f_AlphaToGF * DeltaRho[StandardModel::EW1QCD2] << std::endl;
        std::cout << "  DeltaRho[EW2]=" << DeltaRho[StandardModel::EW2] << " "
                << pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2] << std::endl;
        std::cout << "  DeltaRho[EW2QCD1]=" << DeltaRho[StandardModel::EW2QCD1] << " "
                << pow(f_AlphaToGF, 2.0) * DeltaRho[StandardModel::EW2QCD1] << std::endl;
        std::cout << "  DeltaRho[EW3]=" << DeltaRho[StandardModel::EW3] << " "
                << pow(f_AlphaToGF, 3.0) * DeltaRho[StandardModel::EW3] << std::endl;
        std::cout << "Remainder contributions: alpha or Gmu" << std::endl;
        std::cout << "  DeltaRbar_rem[EW1]=" << DeltaRbar_rem << " "
                << f_AlphaToGF * DeltaRbar_rem << std::endl;
        std::cout << "  deltaRho_rem[EW1]=" << deltaRho_rem[StandardModel::EW1] - ZbbSubtract << " "
                << f_AlphaToGF * (deltaRho_rem[StandardModel::EW1] - ZbbSubtract) << std::endl;
        std::cout << "  deltaRho_rem[EW1QCD1]=" << deltaRho_rem[StandardModel::EW1QCD1] << " "
                << f_AlphaToGF * deltaRho_rem[StandardModel::EW1QCD1] << std::endl;
        std::cout << "  deltaRho_rem[EW2]=" << deltaRho_rem[StandardModel::EW2] << " "
                << pow(f_AlphaToGF, 2.0) * deltaRho_rem[StandardModel::EW2] << std::endl;
        if (bool_Zbb) {
            std::cout << "Taub: alpha or Gmu" << std::endl;
            std::cout << "  taub[EW1]=" << taub[StandardModel::EW1] << " "
                    << f_AlphaToGF * taub[StandardModel::EW1] << std::endl;
            std::cout << "  taub[EW1QCD1]=" << taub[StandardModel::EW1QCD1] << " "
                    << f_AlphaToGF * taub[StandardModel::EW1QCD1] << std::endl;
            std::cout << "  taub[EW2]=" << taub[StandardModel::EW2] << " "
                    << pow(f_AlphaToGF, 2.0) * taub[StandardModel::EW2] << std::endl;
        }
        std::cout << "Each order: alpha or Gmu" << std::endl;
        double dRho_EW1 = DeltaRho[StandardModel::EW1] + deltaRho_rem[StandardModel::EW1] - ZbbSubtract;
        double dRho_EW1QCD1 = DeltaRho[StandardModel::EW1QCD1] + deltaRho_rem[StandardModel::EW1QCD1];
        double dRho_EW2 = DeltaRho[StandardModel::EW1] * DeltaRho[StandardModel::EW1]
                + (deltaRho_rem[StandardModel::EW1] - ZbbSubtract) * DeltaRho[StandardModel::EW1]
                - DeltaRho[StandardModel::EW1] * DeltaRbar_rem
                + DeltaRho[StandardModel::EW2]
                + deltaRho_rem[StandardModel::EW2];
        double dRho_EW1QCD2 = DeltaRho[StandardModel::EW1QCD2];
        double dRho_EW2QCD1 = DeltaRho[StandardModel::EW2QCD1] + deltaRho_rem[StandardModel::EW1QCD1] * DeltaRho[StandardModel::EW1];
        double dRho_EW3 = DeltaRho[StandardModel::EW3];
        //
        double dRho_EW1_TMP = dRho_EW1;
        double dRho_EW1QCD1_TMP = dRho_EW1QCD1;
        double dRho_EW2_TMP = dRho_EW2;
        double dRho_EW1QCD2_TMP = dRho_EW1QCD2;
        double dRho_EW2QCD1_TMP = dRho_EW2QCD1;
        double dRho_EW3_TMP = dRho_EW3;
        //
        if (bool_Zbb) {
            dRho_EW1 = dRho_EW1_TMP + 2.0 * taub[StandardModel::EW1];
            dRho_EW1QCD1 = dRho_EW1QCD1_TMP + 2.0 * taub[StandardModel::EW1QCD1];
            dRho_EW2 = dRho_EW2_TMP + 2.0 * taub[StandardModel::EW2] + taub[StandardModel::EW1] * taub[StandardModel::EW1]
                    + dRho_EW1_TMP * 2.0 * taub[StandardModel::EW1];
            dRho_EW1QCD2 = dRho_EW1QCD2_TMP;
            dRho_EW2QCD1 = dRho_EW2QCD1_TMP + dRho_EW1_TMP * 2.0 * taub[StandardModel::EW1QCD1]
                    + dRho_EW1QCD1_TMP * 2.0 * taub[StandardModel::EW1] + 2.0 * taub[StandardModel::EW1] * taub[StandardModel::EW1QCD1];
            dRho_EW3 = dRho_EW3_TMP + dRho_EW1_TMP * 2.0 * taub[StandardModel::EW2]
                    + dRho_EW2_TMP * 2.0 * taub[StandardModel::EW1] + 2.0 * taub[StandardModel::EW2] * taub[StandardModel::EW1];
        }
        std::cout << "  EW1: " << dRho_EW1 << " " << f_AlphaToGF * dRho_EW1 << std::endl;
        std::cout << "  EW1QCD1: " << dRho_EW1QCD1 << " " << f_AlphaToGF * dRho_EW1QCD1 << std::endl;
        std::cout << "  EW2: " << dRho_EW2 << " " << pow(f_AlphaToGF, 2.0) * dRho_EW2 << std::endl;
        std::cout << "  EW1QCD2: " << dRho_EW1QCD2 << " " << f_AlphaToGF * dRho_EW1QCD2 << std::endl;
        std::cout << "  EW2QCD1: " << dRho_EW2QCD1 << " " << pow(f_AlphaToGF, 2.0) * dRho_EW2QCD1 << std::endl;
        std::cout << "  EW3: " << dRho_EW3 << " " << pow(f_AlphaToGF, 3.0) * dRho_EW3 << std::endl;
        std::cout << "Total contribution: alpha or Gmu" << std::endl;
        std::cout << "  rhoZ="
                << 1.0 + dRho_EW1 + dRho_EW1QCD1 + dRho_EW2
                + dRho_EW1QCD2 + dRho_EW2QCD1 + dRho_EW3
                << " "
                << 1.0 + f_AlphaToGF * dRho_EW1 + f_AlphaToGF * dRho_EW1QCD1
                + pow(f_AlphaToGF, 2.0) * dRho_EW2
                + f_AlphaToGF * dRho_EW1QCD2 + pow(f_AlphaToGF, 2.0) * dRho_EW2QCD1
                + pow(f_AlphaToGF, 3.0) * dRho_EW3
                << std::endl;
        if (bool_Zbb) {
            std::cout << "  rhoZ(taub resummed)="
                    << (1.0 + dRho_EW1_TMP + dRho_EW1QCD1_TMP + dRho_EW2_TMP
                    + dRho_EW1QCD2_TMP + dRho_EW2QCD1_TMP + dRho_EW3_TMP)
                    * pow(1.0 + taub[StandardModel::EW1] + taub[StandardModel::EW1QCD1] + taub[StandardModel::EW2], 2.0)
                    << " "
                    << (1.0 + f_AlphaToGF * dRho_EW1_TMP + f_AlphaToGF * dRho_EW1QCD1_TMP
                    + pow(f_AlphaToGF, 2.0) * dRho_EW2_TMP
                    + f_AlphaToGF * dRho_EW1QCD2_TMP + pow(f_AlphaToGF, 2.0) * dRho_EW2QCD1_TMP
                    + pow(f_AlphaToGF, 3.0) * dRho_EW3_TMP)
                    * pow(1.0 + f_AlphaToGF * taub[StandardModel::EW1] + f_AlphaToGF * taub[StandardModel::EW1QCD1]
                    + pow(f_AlphaToGF, 2.0) * taub[StandardModel::EW2], 2.0)
                    << std::endl;
        }
    } else
        std::cout << "EWSM_Output::outputEachDeltaRhoZ(): Not implemented for schemeRhoZ="
            << getFlagRhoZ() << std::endl;
}

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.

{
    std::cout << "================================================" << std::endl;
    std::cout << "rhoZ_l[(QCD::lepton)" << l << "]" << std::endl;
    std::cout << "Mw(input)   = " << Mw_i << std::endl;
 
    double cW2_TMP = Mw_i * Mw_i / getMz() / getMz();
    double sW2_TMP = 1.0 - cW2_TMP;
 
    double DeltaRho[StandardModel::orders_EW_size];
    DeltaRho[StandardModel::EW1] = getMyOneLoopEW()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW1QCD1] = getMyTwoLoopQCD()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW1QCD2] = getMyThreeLoopQCD()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW2] = getMyTwoLoopEW()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW2QCD1] = getMyThreeLoopEW2QCD()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW3] = getMyThreeLoopEW()->DeltaRho(Mw_i);
 
    /* compute delta rho_rem^f */
    gslpp::complex deltaRho_rem_f[StandardModel::orders_EW_size];
    deltaRho_rem_f[StandardModel::EW1] = getMyOneLoopEW()->deltaRho_rem_f(getLeptons(l), Mw_i);
#ifdef WITHIMTWOLOOPQCD
    deltaRho_rem_f[StandardModel::EW1QCD1] = gslpp::complex(getMyTwoLoopQCD()->deltaRho_rem_f(getLeptons(l), Mw_i).real(),
            getMyTwoLoopQCD()->deltaRho_rem_f(getLeptons(l), Mw_i).imag(), false);
#else
    deltaRho_rem_f[StandardModel::EW1QCD1] = gslpp::complex(getMyTwoLoopQCD()->deltaRho_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
#endif
    deltaRho_rem_f[StandardModel::EW1QCD2] = gslpp::complex(getMyThreeLoopQCD()->deltaRho_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
    deltaRho_rem_f[StandardModel::EW2] = gslpp::complex(getMyTwoLoopEW()->deltaRho_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
    deltaRho_rem_f[StandardModel::EW2QCD1] = gslpp::complex(getMyThreeLoopEW2QCD()->deltaRho_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
    deltaRho_rem_f[StandardModel::EW3] = gslpp::complex(getMyThreeLoopEW()->deltaRho_rem_f(getLeptons(l), Mw_i).real(), 0.0, false);
 
    /* compute Delta rbar_rem */
    double DeltaRbar_rem = getMyOneLoopEW()->DeltaRbar_rem(Mw_i);
 
    double f_AlphaToGF = sqrt(2.0) * getGF() * pow(getMz(), 2.0)
            * sW2_TMP * cW2_TMP / M_PI / getAle();
 
    /* Re[rho_Z^f] with or without resummation */
    double deltaRho_rem_f_real[StandardModel::orders_EW_size];
    for (int j = 0; j < StandardModel::orders_EW_size; ++j)
        deltaRho_rem_f_real[j] = deltaRho_rem_f[j].real();
 
    double dummy[StandardModel::orders_EW_size];
    outputEachDeltaRhoZ(f_AlphaToGF, DeltaRho, deltaRho_rem_f_real,
            DeltaRbar_rem, false, dummy, 0.0);
 
    /* Im[rho_Z^f] without resummation */
    double ImRhoZf = 0.0;
    for (int j = 0; j < StandardModel::orders_EW_size; ++j)
        ImRhoZf += deltaRho_rem_f[j].imag();
    std::cout << "ImRhoZf(with alpha)=" << ImRhoZf << std::endl;
 
    std::cout << "================================================" << std::endl;
}

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.

{
    std::cout << "================================================" << std::endl;
    std::cout << "rhoZ_q[(QCD::quark)" << q << "]" << std::endl;
    std::cout << "Mw(input)   = " << Mw_i << std::endl;
 
    double cW2_TMP = Mw_i * Mw_i / getMz() / getMz();
    double sW2_TMP = 1.0 - cW2_TMP;
 
    double DeltaRho[StandardModel::orders_EW_size];
    DeltaRho[StandardModel::EW1] = getMyOneLoopEW()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW1QCD1] = getMyTwoLoopQCD()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW1QCD2] = getMyThreeLoopQCD()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW2] = getMyTwoLoopEW()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW2QCD1] = getMyThreeLoopEW2QCD()->DeltaRho(Mw_i);
    DeltaRho[StandardModel::EW3] = getMyThreeLoopEW()->DeltaRho(Mw_i);
 
    /* compute delta rho_rem^f */
    gslpp::complex deltaRho_rem_f[StandardModel::orders_EW_size];
    deltaRho_rem_f[StandardModel::EW1] = getMyOneLoopEW()->deltaRho_rem_f(getQuarks(q), Mw_i);
#ifdef WITHIMTWOLOOPQCD
    deltaRho_rem_f[StandardModel::EW1QCD1] = gslpp::complex(getMyTwoLoopQCD()->deltaRho_rem_f(getQuarks(q), Mw_i).real(),
            getMyTwoLoopQCD()->deltaRho_rem_f(getQuarks(q), Mw_i).imag(), false);
#else
    deltaRho_rem_f[StandardModel::EW1QCD1] = gslpp::complex(getMyTwoLoopQCD()->deltaRho_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
#endif
    deltaRho_rem_f[StandardModel::EW1QCD2] = gslpp::complex(getMyThreeLoopQCD()->deltaRho_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
    deltaRho_rem_f[StandardModel::EW2] = gslpp::complex(getMyTwoLoopEW()->deltaRho_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
    deltaRho_rem_f[StandardModel::EW2QCD1] = gslpp::complex(getMyThreeLoopEW2QCD()->deltaRho_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
    deltaRho_rem_f[StandardModel::EW3] = gslpp::complex(getMyThreeLoopEW()->deltaRho_rem_f(getQuarks(q), Mw_i).real(), 0.0, false);
 
    /* compute Delta rbar_rem */
    double DeltaRbar_rem = getMyOneLoopEW()->DeltaRbar_rem(Mw_i);
 
    /* conversion factor */
    double f_AlphaToGF = sqrt(2.0) * getGF() * pow(getMz(), 2.0)
            * sW2_TMP * cW2_TMP / M_PI / getAle();
 
    /* Zbb */
    bool bool_Zbb = false;
    if (q == QCD::BOTTOM) bool_Zbb = true;
    double ZbbSubtract = 0.0;
    if (bool_Zbb)
        ZbbSubtract = -getAle() / 4.0 / M_PI / sW2_TMP
            * pow(getMtpole() / Mw_i, 2.0);
    double taub[StandardModel::orders_EW_size];
    double Xt = getMyEWSMcache()->Xt_alpha(Mw_i);
    if (bool_Zbb) {
        taub[StandardModel::EW1] = -2.0 * Xt;
        taub[StandardModel::EW1QCD1] = 2.0 / 3.0 * M_PI * Xt * getMyEWSMcache()->alsMt();
        taub[StandardModel::EW2] = -2.0 * Xt * Xt * getMyTwoLoopEW()->tau_2();
    }
 
    /* Re[rho_Z^f] with or without resummation */
    double deltaRho_rem_f_real[StandardModel::orders_EW_size];
    for (int j = 0; j < StandardModel::orders_EW_size; ++j)
        deltaRho_rem_f_real[j] = deltaRho_rem_f[j].real();
 
    outputEachDeltaRhoZ(f_AlphaToGF, DeltaRho, deltaRho_rem_f_real,
            DeltaRbar_rem, bool_Zbb, taub, ZbbSubtract);
 
    /* Im[rho_Z^f] without resummation */
    double ImRhoZf = 0.0;
    for (int j = 0; j < StandardModel::orders_EW_size; ++j)
        ImRhoZf += deltaRho_rem_f[j].imag();
    std::cout << "ImRhoZf(with alpha)=" << ImRhoZf << std::endl;
 
    std::cout << "================================================" << std::endl;
}

---

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

• EWSM_Output.h
• EWSM_Output.cpp