NPEffectiveBS Class Reference

A model class for new physics in the form the dimension-six effective Lagrangian in the basis of Barbieri and Strumia [27]. More...

#include <NPEffectiveBS.h>

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

Detailed Description

A model class for new physics in the form the dimension-six effective Lagrangian in the basis of Barbieri and Strumia [27].

Author
HEPfit Collaboration

This is a Model class containing parameters and functions associated with the dimension-six effective Lagrangian:

\[ \mathcal{L}_\mathrm{eff} = \mathcal{L}_\mathrm{SM} + \sum_i \frac{C_i}{\Lambda^2} \mathcal{O}_i \]

where only the SM gauge-invariant and flavour-diagonal dimension-six operators that are relevant to the electroweak precision observables are considered:

\begin{align} \mathcal{O}_{WB} &= (H^\dagger\tau^a H) W^a_{\mu\nu}B^{\mu\nu}, & \mathcal{O}_{H} &= | H^\dagger D_\mu H|^2\,, \nonumber\\ \mathcal{O}_{LL} &= \frac{1}{2}(\overline{L}\gamma_\mu\tau^a L)^2\,, & \mathcal{O}^\prime_{HL} &= i(H^\dagger D_\mu\tau^a H)(\overline{L}\gamma^\mu\tau^a L)\,, \nonumber\\ \mathcal{O}^\prime_{HQ} &= i(H^\dagger D_\mu\tau^a H)(\overline{Q}\gamma^\mu\tau^a Q)\,, & \mathcal{O}_{HL} &= i(H^\dagger D_\mu H)(\overline{L}\gamma^\mu L)\,, \nonumber\\ \mathcal{O}_{HQ} &= i(H^\dagger D_\mu H)(\overline{Q}\gamma^\mu Q)\,, & \mathcal{O}_{HE} &= i(H^\dagger D_\mu H)(\overline{E}\gamma^\mu E)\,, \nonumber\\ \mathcal{O}_{HU} &= i(H^\dagger D_\mu H)(\overline{U}\gamma^\mu U)\,, & \mathcal{O}_{HD} &= i(H^\dagger D_\mu H)(\overline{D}\gamma^\mu D)\,, \end{align}

as well as the Hermitian conjugate for operators \(\mathcal{O}^\prime_{HL}\) to \(\mathcal{O}_{HD}\).

Initialization

After creating an instance of the current class with the constructor NPEffectiveBS(), it is required to call the initialization method InitializeModel(). In the Monte Carlo run, the constructor as well as the initialization method are called in InputParser::ReadParameters().

Model parameters

The model parameters of NPEffectiveBS are summarized below:

Label LaTeX symbol Description
cWB_NP \(C_{WB} \) The coefficient of the operator \(\mathcal{O}_{WB}= \big(H^\dagger\tau^aH\big)W_{\mu\nu}^a B^{\mu\nu}\).
cH_NP \(C_H\) The coefficient of the operator \(\mathcal{O}_{H}= \big|H^\dagger D_\mu H\big|^2\).
cLL_NP \(C_{LL} \) The coefficient of the operator \(\mathcal{O}_{LL}=\frac 12 \big(\overline{L}\,\gamma^\mu\tau^a L\big) \big(\overline{L}\,\gamma_\mu\tau^a L\big)\).
cHLp_NP \(C_{HL}^\prime \) The coefficient of the operator \(\mathcal{O}_{HL}^\prime=i \big(H^\dagger D_\mu\tau^a H\big) \big(\overline{L}\,\gamma^\mu\tau^a L\big)\).
cHQp_NP \(C_{HQ}^\prime \) The coefficient of the operator \(\mathcal{O}_{HQ}^\prime=i \big(H^\dagger D_\mu\tau^a H\big) \big(\overline{Q}\,\gamma^\mu\tau^a Q\big)\).
cHL_NP \(C_{HL} \) The coefficient of the operator \(\mathcal{O}_{HL}=i \big(H^\dagger D_\mu H\big) \big(\overline{L}\,\gamma^\mu L\big) \).
cHQ_NP \(C_{HQ} \) The coefficient of the operator \(\mathcal{O}_{HQ}=i \big(H^\dagger D_\mu H\big) \big(\overline{Q}\,\gamma^\mu Q\big) \).
cHE_NP \(C_{HE} \) The coefficient of the operator \(\mathcal{O}_{HE}=i \big(H^\dagger D_\mu H\big) \big(\overline{E}\,\gamma^\mu E\big) \).
cHU_NP \(C_{HU} \) The coefficient of the operator \(\mathcal{O}_{HU}=i \big(H^\dagger D_\mu H\big) \big(\overline{U}\,\gamma^\mu U\big) \).
cHD_NP \(C_{HD} \) The coefficient of the operator \(\mathcal{O}_{HD}=i \big(H^\dagger D_\mu H\big) \big(\overline{D}\,\gamma^\mu D\big) \).
Lambda_NP \(\Lambda \) The new physics scale.

The parameters to be used as inputs when lepton or quark universality are not assumed are the following:

Label LaTeX symbol Description
cWB_NP \(C_{WB} \) The coefficient of the operator \({\cal O}_{WB}=\big(H^\dagger\tau^a H\big)W_{\mu\nu}^a B^{\mu\nu}\).
cH_NP \(C_H\) The coefficient of the operator \({\cal O}_{H}=\big|H^\dagger D_\mu H\big|^2\).
cLL_NP \(C_{LL} \) The coefficient of the operator \({\cal O}_{LL}=\frac 12\big(\overline{L}\,\gamma^\mu\tau^a L\big) \big(\overline{L}\,\gamma_\mu\tau^a L\big)\).
cHL1p_NP, cHL2p_NP, cHL3p_NP \( C_{HL_i}^\prime \) The coefficient of the operator \({\cal O}_{HL_i}^\prime =i\big(H^\dagger D_\mu\tau^a H\big) \big(\overline{L_i}\,\gamma^\mu\tau^a L_i\big)\) for \(i=1,2,3\).
cHQ1p_NP, cHQ2p_NP, cHQ3p_NP \( C_{HQ_i}^\prime \) The coefficient of the operator \({\cal O}_{HQ_i}^\prime =i\big(H^\dagger D_\mu\tau^a H\big) \big(\overline{Q_i}\,\gamma^\mu\tau^a Q_i\big)\) for \(i=1,2,3\).
cHL1_NP, cHL2_NP, cHL3_NP \( C_{HL_i} \) The coefficient of the operator \({\cal O}_{HL_i} =i\big(H^\dagger D_\mu H\big) \big(\overline{L_i}\,\gamma^\mu L_i\big)\) for \(i=1,2,3\).
cHQ1_NP, cHQ2_NP, cHQ3_NP \( C_{HQ_i} \) The coefficient of the operator \({\cal O}_{HQ_i} =i\big(H^\dagger D_\mu H\big) \big(\overline{Q_i}\,\gamma^\mu Q_i\big)\) for \(i=1,2,3\).
cHE1_NP, cHE2_NP, cHE3_NP \( C_{HE_i} \) The coefficient of the operator \({\cal O}_{HE_i} =i\big(H^\dagger D_\mu H\big) \big(\overline{E_i}\,\gamma^\mu E_i\big)\) for \(i=1,2,3\).
cHU1_NP, cHU2_NP, cHU3_NP \( C_{HU_i} \) The coefficient of the operator \({\cal O}_{HU_i} =i\big(H^\dagger D_\mu H\big) \big(\overline{U_i}\,\gamma^\mu U_i\big)\) for \(i=1,2,3\).
cHD1_NP, cHD2_NP, cHD3_NP \( C_{HD_i} \) The coefficient of the operator \({\cal O}_{HD_i} =i\big(H^\dagger D_\mu H\big) \big(\overline{D_i}\,\gamma^\mu D_i\big)\) for \(i=1,2,3\).
Lambda_NP \(\Lambda \) The new physics scale.

Model flags

There is no model flag in the current class.

Important member functions

Compared to the base class NPbase, the functions for the following quantities are reimplemented in the current class:

  • \(v\)  (with v()),
  • \(M_{W}^{\mathrm{tree}}\)  (with Mw_tree()),
  • \(\Delta G\)  (with DeltaGF()),
  • \(S\), \(T\) and \(U\)   (with obliqueS(), obliqueT() and obliqueU()),
  • \(\Gamma_W\)   (with GammaW()),
  • \(\delta g_V^f\)  (with deltaGVl() and deltaGVq()),
  • \(\delta g_A^f\)  (with deltaGAl() and deltaGAq()),

and also the get methods getCoeff() to retrieve the value of each coefficient.

Definition at line 268 of file NPEffectiveBS.h.

Public Member Functions

virtual bool CheckParameters (const std::map< std::string, double > &DPars)
 A method to check if all the mandatory parameters for the model have been provided in model initialization. More...
 
virtual double deltaGA_f (const Particle f) const
 New physics contribution to the neutral-current axial-vector coupling \(g_A^f\). More...
 
virtual double DeltaGF () const
 New physics contribution to the Fermi constant. More...
 
virtual double deltaGV_f (const Particle f) const
 New physics contribution to the neutral-current vector coupling \(g_V^f\). More...
 
virtual double GammaW () const
 The total width of the \(W\) boson, \(\Gamma_W\). More...
 
virtual double Mw_tree () const
 The tree-level mass of the \(W\) boson, \(M_W^{\mathrm{tree}}\). More...
 
 NPEffectiveBS (const bool FlagLeptonUniversal_in=false, const bool FlagQuarkUniversal_in=false)
 Constructor. More...
 
virtual double obliqueS () const
 The oblique parameter \(S\). More...
 
virtual double obliqueT () const
 The oblique parameter \(T\). More...
 
virtual double obliqueU () const
 The oblique parameter \(U\). More...
 
virtual double v () const
 The SM Higgs vacuum expectation value \(v\). More...
 
- Public Member Functions inherited from NPbase
virtual double A_f (const Particle f) const
 The left-right asymmetry in \(e^+e^-\to Z\to f \bar{f}\) at the \(Z\)-pole, \(\mathcal{A}_f\). More...
 
virtual double AFB (const Particle f) const
 The forward-backward asymmetry in \(e^+e^-\to Z\to f \bar{f}\) at the \(Z\)-pole, \(A^f_{FB}\). More...
 
virtual double BrHbbRatio () const
 The ratio of the Br \((H\to b\bar{b})\) in the current model and in the Standard Model. More...
 
virtual double BrHccRatio () const
 The ratio of the Br \((H\to c\bar{c})\) in the current model and in the Standard Model. More...
 
virtual double BrHgagaRatio () const
 The ratio of the Br \((H\to \gamma\gamma)\) in the current model and in the Standard Model. More...
 
virtual double BrHggRatio () const
 The ratio of the Br \((H\to gg)\) in the current model and in the Standard Model. More...
 
virtual double BrHtautauRatio () const
 The ratio of the Br \((H\to \tau^+\tau^-)\) in the current model and in the Standard Model. More...
 
virtual double BrHWWRatio () const
 The ratio of the Br \((H\to WW)\) in the current model and in the Standard Model. More...
 
virtual double BrHZgaRatio () const
 The ratio of the Br \((H\to Z\gamma)\) in the current model and in the Standard Model. More...
 
virtual double BrHZZRatio () const
 The ratio of the Br \((H\to ZZ)\) in the current model and in the Standard Model. More...
 
virtual double computeGammaTotalRatio () const
 The ratio of the \(\Gamma(H)\) in the current model and in the Standard Model. More...
 
virtual double deltaA_f (const Particle f) const
 The new physics contribution to the left-right asymmetry in \(e^+e^-\to Z\to f \bar{f}\) at the \(Z\)-pole, \(\delta \mathcal{A}_f\). More...
 
virtual double deltaAFB (const Particle f) const
 The new physics contribution to the forward-backward asymmetry in \(e^+e^-\to Z\to f \bar{f}\) at the \(Z\)-pole, \(\delta A^f_{FB}\). More...
 
virtual double deltaGamma_Z () const
 The new physics contribution to the total decay width of the \(Z\) boson, \(\delta \Gamma_Z\). More...
 
virtual double deltaR0_f (const Particle f) const
 The new physics contribution to the ratio \(R_\ell^0=\Gamma_{\mathrm{had}}/\Gamma_\ell\) or \(R_q^0=\Gamma_q/\Gamma_{\mathrm{had}}\), for leptons or quarks, respectively. More...
 
virtual double deltaSigmaHadron () const
 The new physics contribution to the cross section for the process \(e^+ e^-\to Z\to \mathrm{hadrons}\) at the \(Z\) pole, \(\delta \sigma_h^0\). More...
 
virtual double deltaSin2thetaEff_e () const
 The new physics contribution to the effective leptonic weak angle \(\delta \sin^2\theta_{\rm eff}^{\rm lept}\) at the \(Z\) pole. More...
 
virtual gslpp::complex gA_f (const Particle f) const
 The total (SM+NP) contribution to the neutral-current axial-vector coupling \(g_A^f\). More...
 
virtual double Gamma_Z () const
 The total decay width of the \(Z\) boson, \(\Gamma_Z\). More...
 
virtual StandardModel getTrueSM () const
 A method to return a StandardModel object from NPbase. More...
 
virtual gslpp::complex gV_f (const Particle f) const
 The total (SM+NP) contribution to the neutral-current vector coupling \(g_V^f\). More...
 
virtual gslpp::complex kappaZ_f (const Particle f) const
 The effective neutral-current coupling \(\kappa_Z^f\) including SM plus NP contributions. More...
 
virtual double mueeZH (const double sqrt_s) const
 The ratio \(\mu_{eeZH}\) between the \( e^{+}e^{-}\to ZH \) associated production cross-section in the current model and in the Standard Model. More...
 
virtual double muggH (const double sqrt_s) const
 The ratio \(\mu_{ggH}\) between the gluon-gluon fusion Higgs production cross-section in the current model and in the Standard Model. More...
 
virtual double muggHpttH (const double sqrt_s) const
 The ratio \(\mu_{ggH+ttH}\) between the sum of gluon-gluon fusion and t-tbar-Higgs associated production cross-section in the current model and in the Standard Model. More...
 
virtual double muttH (const double sqrt_s) const
 The ratio \(\mu_{ttH}\) between the t-tbar-Higgs associated production cross-section in the current model and in the Standard Model. More...
 
virtual double muVBF (const double sqrt_s) const
 The ratio \(\mu_{VBF}\) between the vector-boson fusion Higgs production cross-section in the current model and in the Standard Model. More...
 
virtual double muVBFpVH (const double sqrt_s) const
 The ratio \(\mu_{VBF+VH}\) between the sum of VBF and WH+ZH associated production cross-section in the current model and in the Standard Model. More...
 
virtual double muVH (const double sqrt_s) const
 The ratio \(\mu_{VH}\) between the WH+ZH associated production cross-section in the current model and in the Standard Model. More...
 
virtual double muWH (const double sqrt_s) const
 The ratio \(\mu_{WH}\) between the W-Higgs associated production cross-section in the current model and in the Standard Model. More...
 
virtual double muZH (const double sqrt_s) const
 The ratio \(\mu_{ZH}\) between the Z-Higgs associated production cross-section in the current model and in the Standard Model. More...
 
virtual double Mw () const
 The mass of the \(W\) boson, \(M_W\). More...
 
 NPbase ()
 The default constructor. More...
 
virtual bool PostUpdate ()
 The post-update method for NPbase. More...
 
virtual double R0_f (const Particle f) const
 The ratio \(R_\ell^0=\Gamma_{\mathrm{had}}/\Gamma_\ell\) or \(R_q^0=\Gamma_q/\Gamma_{\mathrm{had}}\), for leptons or quarks, respectively. More...
 
virtual gslpp::complex rhoZ_f (const Particle f) const
 The effective neutral-current coupling \(\rho_Z^f\) including SM plus NP contributions. More...
 
virtual double sigma0_had () const
 The cross section for the process \(e^+ e^-\to Z\to \mathrm{hadrons}\) at the \(Z\) pole, \(\sigma_h^0\). More...
 
virtual double sin2thetaEff (const Particle f) const
 The leptonic effective weak mixing angle \(\sin^2\theta_{\rm eff}^{\rm lept}\) at the the \(Z\) pole. More...
 
- Public Member Functions inherited from StandardModel
double ale_OS (const double mu, orders order=FULLNLO) const
 The running electromagnetic coupling \(\alpha(\mu)\) in the on-shell schem. More...
 
double alphaMz () const
 The electromagnetic coupling at the \(Z\)-mass scale, \(\alpha(M_Z^2)=\alpha/(1-\Delta\alpha(M_Z^2))\). More...
 
double c02 () const
 The square of the cosine of the weak mixing angle \(c_0^2\) defined without weak radiative corrections. More...
 
virtual bool CheckFlags () const
 A method to check the sanity of the set of model flags. More...
 
bool checkSMparamsForEWPO ()
 A method to check whether the parameters relevant to the EWPO are updated. More...
 
double computeAlpha () const
 The CKM angle \(\alpha\). More...
 
double computeBeta () const
 The CKM angle \(\beta\). More...
 
double computeBetas () const
 The CKM angle \(\beta_s\). More...
 
double computeBrHtobb () const
 The Br \((H\to bb)\) in the Standard Model. More...
 
double computeBrHtocc () const
 The Br \((H\to cc)\) in the Standard Model. More...
 
double computeBrHtogaga () const
 The Br \((H\to\gamma\gamma)\) in the Standard Model. More...
 
double computeBrHtogg () const
 The Br \((H\to gg)\) in the Standard Model. More...
 
double computeBrHtotautau () const
 The Br \((H\to \tau\tau)\) in the Standard Model. More...
 
double computeBrHtoWW () const
 The Br \((H\to WW)\) in the Standard Model. More...
 
double computeBrHtoZga () const
 The Br \((H\to Z\gamma)\) in the Standard Model. More...
 
double computeBrHtoZZ () const
 The Br \((H\to ZZ)\) in the Standard Model. More...
 
void ComputeDeltaR_rem (const double Mw_i, double DeltaR_rem[orders_EW_size]) const
 A method to collect \(\Delta r_{\mathrm{rem}}\) computed via subclasses. More...
 
void ComputeDeltaRho (const double Mw_i, double DeltaRho[orders_EW_size]) const
 A method to collect \(\Delta\rho\) computed via subclasses. More...
 
double computeGamma () const
 The CKM angle \(\gamma\). More...
 
double computeGammaHgaga_tt () const
 The top loop contribution to \(H\to\gamma\gamma\) in the Standard Model. More...
 
double computeGammaHgaga_tW () const
 The mixed \(t-W\) loop contribution to \(H\to\gamma\gamma\) in the Standard Model. More...
 
double computeGammaHgaga_WW () const
 The \(W\) loop contribution to \(H\to\gamma\gamma\) in the Standard Model. More...
 
double computeGammaHgg_bb () const
 The bottom loop contribution to \(H\to gg\) in the Standard Model. More...
 
double computeGammaHgg_tb () const
 The top-bottom interference contribution to \(H\to gg\) in the Standard Model. More...
 
double computeGammaHgg_tt () const
 The top loop contribution to \(H\to gg\) in the Standard Model. More...
 
double computeGammaHTotal () const
 The Higgs total width in the Standard Model. More...
 
double computeGammaHZga_tt () const
 The top loop contribution to \(H\to Z\gamma\) in the Standard Model. More...
 
double computeGammaHZga_tW () const
 The mixed \(t-W\) loop contribution to \(H\to Z\gamma\) in the Standard Model. More...
 
double computeGammaHZga_WW () const
 The \(W\) loop contribution to \(H\to Z\gamma\) in the Standard Model. Currently it returns the value of tab 41 in ref. [89]. More...
 
gslpp::complex computelamc () const
 The product of the CKM elements \(V_{cd} V_{cs}^*\). More...
 
gslpp::complex computelamc_d () const
 The product of the CKM elements \(V_{cd} V_{cb}^*\). More...
 
gslpp::complex computelamc_s () const
 The product of the CKM elements \(V_{cs} V_{cb}^*\). More...
 
gslpp::complex computelamt () const
 The product of the CKM elements \(V_{td} V_{ts}^*\). More...
 
gslpp::complex computelamt_d () const
 The product of the CKM elements \(V_{td} V_{tb}^*\). More...
 
gslpp::complex computelamt_s () const
 The product of the CKM elements \(V_{ts} V_{tb}^*\). More...
 
gslpp::complex computelamu () const
 The product of the CKM elements \(V_{ud} V_{us}^*\). More...
 
gslpp::complex computelamu_d () const
 The product of the CKM elements \(V_{ud} V_{ub}^*\). More...
 
gslpp::complex computelamu_s () const
 The product of the CKM elements \(V_{us} V_{ub}^*\). More...
 
double computeRb () const
 \(R_b=|(V_{ud}V_{ub}^*)/(V_{ud}V_{ub}^*)|\). More...
 
double computeRt () const
 \(R_t=|(V_{td} V_{tb}^*)/(V_{cd}V_{cb}^*)|\). More...
 
double computeRts () const
 \(R_{ts}=|(V_{ts}V_{tb}^*)/(V_{cs}V_{cb}^*)|\). More...
 
double computeSigmaggH (const double sqrt_s) const
 The ggH cross section in the Standard Model. More...
 
double computeSigmaggH_bb (const double sqrt_s) const
 The square of the bottom-quark contribution to the ggH cross section in the Standard Model. More...
 
double computeSigmaggH_tb (const double sqrt_s) const
 The top-bottom interference contribution to the ggH cross section in the Standard Model. More...
 
double computeSigmaggH_tt (const double sqrt_s) const
 The square of the top-quark contribution to the ggH cross section in the Standard Model. More...
 
double computeSigmattH (const double sqrt_s) const
 The ttH production cross section in the Standard Model. More...
 
double computeSigmaVBF (const double sqrt_s) const
 The VBF cross section in the Standard Model. More...
 
double computeSigmaWF (const double sqrt_s) const
 The W fusion contribution \(\sigma_{WF}\) to higgs-production cross section in the Standard Model. More...
 
double computeSigmaWH (const double sqrt_s) const
 The WH production cross section in the Standard Model. More...
 
double computeSigmaZF (const double sqrt_s) const
 The Z fusion contribution \(\sigma_{ZF}\) to higgs-production cross section in the Standard Model. More...
 
double computeSigmaZH (const double sqrt_s) const
 The ZH production cross section in the Standard Model. More...
 
double computeSigmaZWF (const double sqrt_s) const
 The Z W interference fusion contribution \(\sigma_{ZWF}\) to higgs-production cross section in the Standard Model. More...
 
virtual double cW2 (const double Mw_i) const
 The square of the cosine of the weak mixing angle in the on-shell scheme, denoted as \(c_W^2\). More...
 
virtual double cW2 () const
 
double DeltaAlpha () const
 The total corrections to the electromagnetic coupling \(\alpha\) at the \(Z\)-mass scale, denoted as \(\Delta\alpha(M_Z^2)\). More...
 
double DeltaAlphaL5q () const
 The sum of the leptonic and the five-flavour hadronic corrections to the electromagnetic coupling \(\alpha\) at the \(Z\)-mass scale, denoted as \(\Delta\alpha^{\ell+5q}(M_Z^2)\). More...
 
double DeltaAlphaLepton (const double s) const
 Leptonic contribution to the electromagnetic coupling \(\alpha\), denoted as \(\Delta\alpha_{\mathrm{lept}}(s)\). More...
 
double DeltaAlphaTop (const double s) const
 Top-quark contribution to the electromagnetic coupling \(\alpha\), denoted as \(\Delta\alpha_{\mathrm{top}}(s)\). More...
 
virtual gslpp::complex deltaKappaZ_f (const Particle f) const
 Flavour non-universal vertex corrections to \(\kappa_Z^l\), denoted by \(\Delta\kappa_Z^l\). More...
 
virtual double DeltaR () const
 The SM prediction for \(\Delta r\) derived from that for the \(W\) boson mass. More...
 
virtual double DeltaRbar () const
 The SM prediction for \(\Delta \overline{r}\) derived from that for the \(W\)-boson mass. More...
 
virtual gslpp::complex deltaRhoZ_f (const Particle f) const
 Flavour non-universal vertex corrections to \(\rho_Z^l\), denoted by \(\Delta\rho_Z^l\). More...
 
virtual double epsilon1 () const
 The SM contribution to the epsilon parameter \(\varepsilon_1\). More...
 
virtual double epsilon2 () const
 The SM contribution to the epsilon parameter \(\varepsilon_2\). More...
 
virtual double epsilon3 () const
 The SM contribution to the epsilon parameter \(\varepsilon_3\). More...
 
virtual double epsilonb () const
 The SM contribution to the epsilon parameter \(\varepsilon_b\). More...
 
virtual double Gamma_had () const
 The hadronic decay width of the \(Z\) boson, \(\Gamma_{h}\). More...
 
virtual double Gamma_inv () const
 The invisible partial decay width of the \(Z\) boson, \(\Gamma_{\mathrm{inv}}\). More...
 
virtual double GammaW (const Particle fi, const Particle fj) const
 A partial decay width of the \(W\) boson decay into a SM fermion pair. More...
 
virtual double GammaZ (const Particle f) const
 The \(Z\to \ell\bar{\ell}\) partial decay width, \(\Gamma_\ell\). More...
 
double getA () const
 A get method to retrieve the CKM element \(A\). More...
 
double getAle () const
 A get method to retrieve the fine-structure constant \(\alpha\). More...
 
double getAlsMz () const
 A get method to access the value of \(\alpha_s(M_Z)\). More...
 
CKM getCKM () const
 A get method to retrieve the member object of type CKM. More...
 
double getDAle5Mz () const
 A get method to retrieve the five-flavour hadronic contribution to the electromagnetic coupling, \(\Delta\alpha_{\mathrm{had}}^{(5)}(M_Z^2)\). More...
 
double getDelGammaZ () const
 A get method to retrieve the theoretical uncertainty in \(\Gamma_Z\), denoted as \(\delta\,\Gamma_Z\). More...
 
double getDelMw () const
 A get method to retrieve the theoretical uncertainty in \(M_W\), denoted as \(\delta\,M_W\). More...
 
double getDelSin2th_l () const
 A get method to retrieve the theoretical uncertainty in \(\sin^2\theta_{\rm eff}^{\rm lept}\), denoted as \(\delta\sin^2\theta_{\rm eff}^{\rm lept}\). More...
 
double getDeltaMK () const
 
double getDmk () const
 
double getEpsK () const
 
double getEtab () const
 A get method to retrieve the CKM element \(\bar{\eta}\). More...
 
std::string getFlagKappaZ () const
 A method to retrieve the model flag KappaZ. More...
 
std::string getFlagMw () const
 A method to retrieve the model flag Mw. More...
 
std::string getFlagRhoZ () const
 A method to retrieve the model flag RhoZ. More...
 
double getGF () const
 A get method to retrieve the Fermi constant \(G_\mu\). More...
 
int getIterationNo () const
 
double getKbarEpsK () const
 
double getLambda () const
 A get method to retrieve the CKM element \(\lambda\). More...
 
Particle getLeptons (const StandardModel::lepton p) const
 A get method to retrieve the member object of a lepton. More...
 
double getMHl () const
 A get method to retrieve the Higgs mass \(m_h\). More...
 
double getMuw () const
 A get method to retrieve the matching scale \(\mu_W\) around the weak scale. More...
 
EWSMApproximateFormulaegetMyApproximateFormulae () const
 A get method to retrieve the member pointer of type EWSMApproximateFormulae. More...
 
EWSMcachegetMyEWSMcache () const
 A get method to retrieve the member pointer of type EWSMcache. More...
 
FlavourgetMyFlavour () const
 
LeptonFlavourgetMyLeptonFlavour () const
 
virtual StandardModelMatchinggetMyMatching () const
 A get method to access the member pointer of type StandardModelMatching. More...
 
EWSMOneLoopEWgetMyOneLoopEW () const
 A get method to retrieve the member pointer of type EWSMOneLoopEW,. More...
 
EWSMThreeLoopEWgetMyThreeLoopEW () const
 
EWSMThreeLoopEW2QCDgetMyThreeLoopEW2QCD () const
 
EWSMThreeLoopQCDgetMyThreeLoopQCD () const
 
EWSMTwoFermionsLEP2 * getMyTwoFermionsLEP2 () const
 A get method to retrieve the member pointer of type EWSMTwoFermionsLEP2. More...
 
EWSMTwoLoopEWgetMyTwoLoopEW () const
 
EWSMTwoLoopQCDgetMyTwoLoopQCD () const
 
double getMz () const
 A get method to access the mass of the \(Z\) boson \(M_Z\). More...
 
double getphiEpsK () const
 
double getRhob () const
 A get method to retrieve the CKM element \(\bar{\rho}\). More...
 
double getSM_M12D () const
 
gslpp::matrix< gslpp::complexgetUPMNS () const
 A get method to retrieve the object of the PMNS matrix. More...
 
gslpp::matrix< gslpp::complexgetVCKM () const
 A get method to retrieve the CKM matrix. More...
 
gslpp::matrix< gslpp::complexgetYd () const
 A get method to retrieve the Yukawa matrix of the down-type quarks, \(Y_d\). More...
 
gslpp::matrix< gslpp::complexgetYe () const
 A get method to retrieve the Yukawa matrix of the charged leptons, \(Y_e\). More...
 
gslpp::matrix< gslpp::complexgetYn () const
 A get method to retrieve the Yukawa matrix of the neutrinos, \(Y_\nu\). More...
 
gslpp::matrix< gslpp::complexgetYu () const
 A get method to retrieve the Yukawa matrix of the up-type quarks, \(Y_u\). More...
 
virtual bool Init (const std::map< std::string, double > &DPars)
 A method to initialize the model parameters. More...
 
virtual bool InitializeModel ()
 A method to initialize the model. More...
 
bool IsFlagNoApproximateGammaZ () const
 A method to retrieve the model flag NoApproximateGammaZ. More...
 
bool IsFlagWithoutNonUniversalVC () const
 A method to retrieve the model flag WithoutNonUniversalVC. More...
 
double MwbarFromMw (const double Mw) const
 A method to convert the \(W\)-boson mass in the experimental/running-width scheme to that in the complex-pole/fixed-width scheme. More...
 
double MwFromMwbar (const double Mwbar) const
 A method to convert the \(W\)-boson mass in the complex-pole/fixed-width scheme to that in the experimental/running-width scheme. More...
 
double Mzbar () const
 The \(Z\)-boson mass \(\overline{M}_Z\) in the complex-pole/fixed-width scheme. More...
 
virtual bool PreUpdate ()
 The pre-update method for StandardModel. More...
 
virtual double rho_GammaW (const Particle fi, const Particle fj) const
 EW radiative corrections to the width of \(W \to f_i \bar{f}_j\), denoted as \(\rho^W_{ij}\). More...
 
double s02 () const
 The square of the sine of the weak mixing angle \(s_0^2\) defined without weak radiative corrections. More...
 
virtual bool setFlag (const std::string name, const bool value)
 A method to set a flag of StandardModel. More...
 
void setFlagCacheInStandardModel (bool FlagCacheInStandardModel)
 A set method to change the model flag CacheInStandardModel of StandardModel. More...
 
void setFlagNoApproximateGammaZ (bool FlagNoApproximateGammaZ)
 
virtual bool setFlagStr (const std::string name, const std::string value)
 A method to set a flag of StandardModel. More...
 
 StandardModel ()
 The default constructor. More...
 
virtual double sW2 (const double Mw_i) const
 The square of the sine of the weak mixing angle in the on-shell scheme, denoted as \(s_W^2\). More...
 
double sW2 () const
 
virtual bool Update (const std::map< std::string, double > &DPars)
 The update method for StandardModel. More...
 
virtual ~StandardModel ()
 The default destructor. More...
 
- Public Member Functions inherited from QCD
double AboveTh (const double mu) const
 The active flavour threshold above the scale \(\mu\) as defined in QCD::Thresholds(). More...
 
double Als (const double mu, const orders order=FULLNLO) const
 Computes the running strong coupling \(\alpha_s(\mu)\) in the \(\overline{\mathrm{MS}}\) scheme. In the cases of LO, NLO and FULLNNLO, the coupling is computed with AlsWithInit(). On the other hand, in the cases of NNLO and FULLNNLO, the coupling is computed with AlsWithLambda(). More...
 
double Als4 (const double mu) const
 The value of \(\alpha_s^{\mathrm{FULLNLO}}\) at any scale \(\mu\) with the number of flavours \(n_f = 4\). More...
 
double AlsWithInit (const double mu, const double alsi, const double mu_i, const orders order) const
 Computes the running strong coupling \(\alpha_s(\mu)\) from \(\alpha_s(\mu_i)\) in the \(\overline{\mathrm{MS}}\) scheme, where it is forbidden to across a flavour threshold in the RG running from \(\mu_i\) to \(\mu\). More...
 
double AlsWithLambda (const double mu, const orders order) const
 Computes the running strong coupling \(\alpha_s(\mu)\) in the \(\overline{\mathrm{MS}}\) scheme with the use of \(\Lambda_{\rm QCD}\). More...
 
double BelowTh (const double mu) const
 The active flavour threshold below the scale \(\mu\) as defined in QCD::Thresholds(). More...
 
double Beta0 (const double nf) const
 The \(\beta_0(n_f)\) coefficient for a certain number of flavours \(n_f\). More...
 
double Beta1 (const double nf) const
 The \(\beta_1(n_f)\) coefficient for a certain number of flavours \(n_f\). More...
 
double Beta2 (const double nf) const
 The \(\beta_2(n_f)\) coefficient for a certain number of flavours \(n_f\). More...
 
double geta_0A0 () const
 
double geta_0A0phi () const
 
double geta_0A1 () const
 
double geta_0A12 () const
 
double geta_0A12phi () const
 
double geta_0A1phi () const
 
double geta_0T1 () const
 
double geta_0T1phi () const
 
double geta_0T2 () const
 
double geta_0T23 () const
 
double geta_0T23phi () const
 
double geta_0T2phi () const
 
double geta_0V () const
 
double geta_0Vphi () const
 
double geta_1A0 () const
 
double geta_1A0phi () const
 
double geta_1A1 () const
 
double geta_1A12 () const
 
double geta_1A12phi () const
 
double geta_1A1phi () const
 
double geta_1T1 () const
 
double geta_1T1phi () const
 
double geta_1T2 () const
 
double geta_1T23 () const
 
double geta_1T23phi () const
 
double geta_1T2phi () const
 
double geta_1V () const
 
double geta_1Vphi () const
 
double geta_2A0 () const
 
double geta_2A0phi () const
 
double geta_2A1 () const
 
double geta_2A12 () const
 
double geta_2A12phi () const
 
double geta_2A1phi () const
 
double geta_2T1 () const
 
double geta_2T1phi () const
 
double geta_2T2 () const
 
double geta_2T23 () const
 
double geta_2T23phi () const
 
double geta_2T2phi () const
 
double geta_2V () const
 
double geta_2Vphi () const
 
double getAlsM () const
 A get method to access the value of \(\alpha_s(M_{\alpha_s})\). More...
 
BParameter getBBd () const
 For getting the bag parameters corresponding to the operator basis \(O_1 -O_5\) in \(\Delta b = 2\) process in the \(B_d\) meson system. More...
 
BParameter getBBs () const
 For getting the bag parameters corresponding to the operator basis \(O_1 -O_5\) in \(\Delta b = 2\) process in the \(B_s\) meson system. More...
 
BParameter getBD () const
 For getting the bag parameters corresponding to the operator basis \(O_1 -O_5\) in \(\Delta c = 2\) process in the \(D^0\) meson system. More...
 
BParameter getBK () const
 For getting the bag parameters corresponding to the operator basis \(O_1 -O_5\) in \(\Delta s = 2\) process in the \(K^0\) meson system. More...
 
BParameter getBKd1 () const
 
BParameter getBKd3 () const
 
double getBLNPcorr () const
 
double getBr_B_Xcenu () const
 
double getBr_Kp_munu () const
 
double getBr_Kp_P0enu () const
 
double getbsgamma_E0 () const
 
double getCF () const
 A get method to access the Casimir factor of QCD. More...
 
double getDeltaP_cu () const
 
double getFKstarp () const
 
double getGambino_BRsem () const
 
double getGambino_Mbkin () const
 
double getGambino_Mcatmuc () const
 
double getGambino_muG2 () const
 
double getGambino_mukin () const
 
double getGambino_mupi2 () const
 
double getGambino_rhoD3 () const
 
double getGambino_rhoLS3 () const
 
gslpp::complex geth_0 () const
 
gslpp::complex geth_0_1 () const
 
gslpp::complex geth_0_1_MP () const
 
gslpp::complex geth_0_2 () const
 
gslpp::complex geth_0_MP () const
 
gslpp::complex geth_m () const
 
gslpp::complex geth_m_1 () const
 
gslpp::complex geth_m_2 () const
 
gslpp::complex geth_p () const
 
gslpp::complex geth_p_1 () const
 
gslpp::complex geth_p_2 () const
 
double getIB_Kl () const
 
double getIB_Kp () const
 
double getm_fit2_f0 () const
 
double getm_fit2_fplus () const
 
double getm_fit2_fT () const
 
double getMAls () const
 A get method to access the mass scale \(M_{\alpha_s}\) at which the strong coupling constant measurement is provided. More...
 
Meson getMesons (const meson m) const
 A get method to access a meson as an object of the type Meson. More...
 
double getMRA0 () const
 
double getMRA0phi () const
 
double getMRA1 () const
 
double getMRA12 () const
 
double getMRA12phi () const
 
double getMRA1phi () const
 
double getMRT1 () const
 
double getMRT1phi () const
 
double getMRT2 () const
 
double getMRT23 () const
 
double getMRT23phi () const
 
double getMRT2phi () const
 
double getMRV () const
 
double getMRVphi () const
 
double getMtpole () const
 A get method to access the pole mass of the top quark. More...
 
double getMub () const
 A get method to access the threshold between five- and four-flavour theory in GeV. More...
 
double getMuc () const
 A get method to access the threshold between four- and three-flavour theory in GeV. More...
 
double getMut () const
 A get method to access the threshold between six- and five-flavour theory in GeV. More...
 
double getNc () const
 A get method to access the number of colours \(N_c\). More...
 
double getOmega_eta_etap () const
 
Particle getQuarks (const quark q) const
 A get method to access a quark as an object of the type Particle. More...
 
double getr_1_fplus () const
 
double getr_1_fT () const
 
double getr_2_f0 () const
 
double getr_2_fplus () const
 
double getr_2_fT () const
 
double getReA0_Kd () const
 
double getReA2_Kd () const
 
double logLambda (const double nf, orders order) const
 Computes \(\ln\Lambda_\mathrm{QCD}\) with nf flavours in GeV. More...
 
double Nf (const double mu) const
 The number of active flavour at scale \(\mu\). More...
 
std::string orderToString (const orders order) const
 Converts an object of the enum type "orders" to the corresponding string. More...
 
 QCD ()
 Constructor. More...
 
void setNc (double Nc)
 A set method to change the number of colours \(N_c\). More...
 
double Thresholds (const int i) const
 For accessing the active flavour threshold scales. More...
 
- Public Member Functions inherited from Model
const double & getModelParam (std::string name) const
 
bool IsModelInitialized () const
 A method to check if the model is initialized. More...
 
bool isModelParam (std::string name) const
 
bool isModelSUSY () const
 
bool isModelTHDM () const
 
bool IsUpdateError () const
 A method to check if there was any error in the model update process. More...
 
 Model ()
 The default constructor. More...
 
std::string ModelName () const
 A method to fetch the name of the model. More...
 
void setModelInitialized (bool ModelInitialized)
 A set method to fix the failure or success of the initialization of the model. More...
 
void setModelName (const std::string name)
 A method to set the name of the model. More...
 
void setModelSUSY ()
 
void setModelTHDM ()
 
void setUpdateError (bool UpdateError)
 A set method to fix the update status as success or failure. More...
 
virtual ~Model ()
 The default destructor. More...
 

Static Public Attributes

static const int NNPEffectiveBSVars = 30
 The number of the model parameters in NPEffectiveBS. More...
 
static const int NNPEffectiveBSVars_LFU = 19
 The number of the model parameters in NPEffectiveBS with lepton flavour universality. More...
 
static const int NNPEffectiveBSVars_LFU_QFU = 11
 The number of the model parameters in NPEffectiveBS with lepton and quark flavour universalities. More...
 
static const int NNPEffectiveBSVars_QFU = 22
 The number of the model parameters in NPEffectiveBS with quark flavour universality. More...
 
static const std::string NPEffectiveBSVars [NNPEffectiveBSVars]
 A string array containing the labels of the model parameters in NPEffectiveBS. More...
 
static const std::string NPEffectiveBSVars_LFU [NNPEffectiveBSVars_LFU]
 A string array containing the labels of the model parameters in NPEffectiveBS with lepton flavour universality. More...
 
static const std::string NPEffectiveBSVars_LFU_QFU [NNPEffectiveBSVars_LFU_QFU]
 A string array containing the labels of the model parameters in NPEffectiveBS with lepton and quark flavour universalities. More...
 
static const std::string NPEffectiveBSVars_QFU [NNPEffectiveBSVars_QFU]
 A string array containing the labels of the model parameters in NPEffectiveBS with quark flavour universality. More...
 
- Static Public Attributes inherited from StandardModel
static const double GeVminus2_to_nb = 389379.338
 
static const double Mw_error = 0.00001
 The target accuracy of the iterative calculation of the \(W\)-boson mass in units of GeV. More...
 
static const int NSMvars = 26
 The number of the model parameters in StandardModel. More...
 
static const int NumSMParamsForEWPO = 27
 The number of the SM parameters that are relevant to the EW precision observables. More...
 
static const std::string SMvars [NSMvars]
 A string array containing the labels of the model parameters in StandardModel. More...
 
- Static Public Attributes inherited from QCD
static const int NQCDvars = 186
 The number of model parameters in QCD. More...
 
static const std::string QCDvars [NQCDvars]
 An array containing the labels under which all QCD parameters are stored in a vector of ModelParameter via InputParser::ReadParameters(). More...
 

Protected Member Functions

virtual void setParameter (const std::string name, const double &value)
 A method to set the value of a parameter of the model. More...
 
- Protected Member Functions inherited from StandardModel
bool checkEWPOscheme (const std::string scheme) const
 A method to check if a given scheme name in string form is valid. More...
 
virtual void computeCKM ()
 The method to compute the CKM matrix. More...
 
virtual void computeYukawas ()
 The method to compute the Yukawa matrices. More...
 
double Delta_EWQCD (const QCD::quark q) const
 The non-factorizable EW-QCD corrections to the partial widths for \(Z\to q\bar{q}\), denoted as \(\Delta_{\mathrm{EW/QCD}}\). More...
 
double RAq (const QCD::quark q) const
 The radiator factor associated with the final-state QED and QCD corrections to the the axial-vector-current interactions, \(R_A^q(M_Z^2)\). More...
 
double resumKappaZ (const double DeltaRho[orders_EW_size], const double deltaKappa_rem[orders_EW_size], const double DeltaRbar_rem, const bool bool_Zbb) const
 A method to compute the real part of the effetvive coupling \(\kappa_Z^f\) from \(\Delta\rho\), \(\delta\rho_{\rm rem}^{f}\) and \(\Delta r_{\mathrm{rem}}\). More...
 
double resumMw (const double Mw_i, const double DeltaRho[orders_EW_size], const double DeltaR_rem[orders_EW_size]) const
 A method to compute the \(W\)-boson mass from \(\Delta\rho\) and \(\Delta r_{\mathrm{rem}}\). More...
 
double resumRhoZ (const double DeltaRho[orders_EW_size], const double deltaRho_rem[orders_EW_size], const double DeltaRbar_rem, const bool bool_Zbb) const
 A method to compute the real part of the effective coupling \(\rho_Z^f\) from \(\Delta\rho\), \(\delta\rho_{\rm rem}^{f}\) and \(\Delta r_{\mathrm{rem}}\). More...
 
double RVh () const
 The singlet vector corrections to the hadronic \(Z\)-boson width, denoted as \(R_V^h\). More...
 
double RVq (const QCD::quark q) const
 The radiator factor associated with the final-state QED and QCD corrections to the the vector-current interactions, \(R_V^q(M_Z^2)\). More...
 
double SchemeToDouble (const std::string scheme) const
 A method to convert a given scheme name in string form into a floating-point number with double precision. More...
 
double taub () const
 Top-mass corrections to the \(Zb\bar{b}\) vertex, denoted by \(\tau_b\). More...
 

Protected Attributes

double cH
 The dimension-6 operator coefficient \(C_{H}\). More...
 
double cHD1
 The dimension-6 operator coefficient \(C_{HD_1}\). More...
 
double cHD2
 The dimension-6 operator coefficient \(C_{HD_2}\). More...
 
double cHD3
 The dimension-6 operator coefficient \(C_{HD_3}\). More...
 
double cHE1
 The dimension-6 operator coefficient \(C_{HE_1}\). More...
 
double cHE2
 The dimension-6 operator coefficient \(C_{HE_2}\). More...
 
double cHE3
 The dimension-6 operator coefficient \(C_{HE_3}\). More...
 
double cHL1
 The dimension-6 operator coefficient \(C_{HL_1}\). More...
 
double cHL1p
 The dimension-6 operator coefficient \(C_{HL_1}^\prime\). More...
 
double cHL2
 The dimension-6 operator coefficient \(C_{HL_2}\). More...
 
double cHL2p
 The dimension-6 operator coefficient \(C_{HL_2}^\prime\). More...
 
double cHL3
 The dimension-6 operator coefficient \(C_{HL_3}\). More...
 
double cHL3p
 The dimension-6 operator coefficient \(C_{HL_3}^\prime\). More...
 
double cHQ1
 The dimension-6 operator coefficient \(C_{HQ_1}\). More...
 
double cHQ1p
 The dimension-6 operator coefficient \(C_{HQ_1}^\prime\). More...
 
double cHQ2
 The dimension-6 operator coefficient \(C_{HQ_2}\). More...
 
double cHQ2p
 The dimension-6 operator coefficient \(C_{HQ_2}^\prime\). More...
 
double cHQ3
 The dimension-6 operator coefficient \(C_{HQ_3}\). More...
 
double cHQ3p
 The dimension-6 operator coefficient \(C_{HQ_3}^\prime\). More...
 
double cHU1
 The dimension-6 operator coefficient \(C_{HU_1}\). More...
 
double cHU2
 The dimension-6 operator coefficient \(C_{HU_2}\). More...
 
double cHU3
 The dimension-6 operator coefficient \(C_{HU_3}\). More...
 
double cL1L1
 The dimension-6 operator coefficient \(C_{L_1L_1}\). More...
 
double cL1L2
 The dimension-6 operator coefficient \(C_{L_1L_2}\). More...
 
double cL1L3
 The dimension-6 operator coefficient \(C_{L_1L_3}\). More...
 
double cL2L2
 The dimension-6 operator coefficient \(C_{L_2L_2}\). More...
 
double cL2L3
 The dimension-6 operator coefficient \(C_{L_2L_3}\). More...
 
double cL3L3
 The dimension-6 operator coefficient \(C_{L_3L_3}\). More...
 
double cWB
 The dimension-6 operator coefficient \(C_{WB}\). More...
 
double Lambda_NP
 The new physics scale \(\Lambda\). More...
 
- Protected Attributes inherited from NPbase
StandardModel trueSM
 
- Protected Attributes inherited from StandardModel
double A
 The CKM parameter \(A\) in the Wolfenstein parameterization. More...
 
double ale
 The fine-structure constant \(\alpha\). More...
 
double AlsMz
 The strong coupling constant at the Z-boson mass, \(\alpha_s(M_Z)\). More...
 
double dAle5Mz
 The five-flavour hadronic contribution to the electromagnetic coupling, \(\Delta\alpha_{\mathrm{had}}^{(5)}(M_Z^2)\). More...
 
double delGammaZ
 The theoretical uncertainty in \(\Gamma_Z\), denoted as \(\delta\,\Gamma_Z\), in GeV. More...
 
double delMw
 The theoretical uncertainty in \(M_W\), denoted as \(\delta\,M_W\), in GeV. More...
 
double delSin2th_l
 The theoretical uncertainty in \(\sin^2\theta_{\rm eff}^{\rm lept}\), denoted as \(\delta\sin^2\theta_{\rm eff}^{\rm lept}\). More...
 
double DeltaMK
 
double Dmk
 
double EpsK
 
double etab
 The CKM parameter \(\bar{\eta}\) in the Wolfenstein parameterization. More...
 
bool flag_order [orders_EW_size]
 An array of internal flags controlling the inclusions of higher-order corrections. More...
 
double GF
 The Fermi constant \(G_\mu\) in \({\rm GeV}^{-2}\). More...
 
double KbarEpsK
 
double lambda
 The CKM parameter \(\lambda\) in the Wolfenstein parameterization. More...
 
Particle leptons [6]
 An array of Particle objects for the leptons. More...
 
double mHl
 The Higgs mass \(m_h\) in GeV. More...
 
double muw
 A matching scale \(\mu_W\) around the weak scale in GeV. More...
 
CKM myCKM
 An object of type CKM. More...
 
double Mz
 The mass of the \(Z\) boson in GeV. More...
 
double phiEpsK
 
bool requireCKM
 An internal flag to control whether the CKM matrix has to be recomputed. More...
 
bool requireYe
 An internal flag to control whether the charged-lepton Yukawa matrix has to be recomputed. More...
 
bool requireYn
 An internal flag to control whether the neutrino Yukawa matrix has to be recomputed. More...
 
double rhob
 The CKM parameter \(\bar{\rho}\) in the Wolfenstein parameterization. More...
 
double SM_M12D
 
gslpp::matrix< gslpp::complexUPMNS
 The PMNS matrix. More...
 
gslpp::matrix< gslpp::complexVCKM
 The CKM matrix. More...
 
gslpp::matrix< gslpp::complexYd
 The Yukawa matrix of the down-type quarks. More...
 
gslpp::matrix< gslpp::complexYe
 The Yukawa matrix of the charged leptons. More...
 
gslpp::matrix< gslpp::complexYn
 The Yukawa matrix of the neutrinos. More...
 
gslpp::matrix< gslpp::complexYu
 The Yukawa matrix of the up-type quarks. More...
 
- Protected Attributes inherited from Model
std::map< std::string, boost::reference_wrapper< const double > > ModelParamMap
 
bool UpdateError
 A boolean set to false if update is successful. More...
 

Private Member Functions

double deltaGL_f_tmp (const Particle f) const
 New physics contribution to \(g_L^f\). More...
 
double deltaGR_f_tmp (const Particle f) const
 New physics contribution to \(g_R^f\). More...
 

Private Attributes

const bool FlagLeptonUniversal
 An internal boolean flag that is true if assuming lepton flavour universality. More...
 
const bool FlagQuarkUniversal
 An internal boolean flag that is true if assuming quark flavour universality. More...
 

Additional Inherited Members

- Public Types inherited from StandardModel
enum  lepton { NEUTRINO_1, ELECTRON, NEUTRINO_2, MU, NEUTRINO_3, TAU }
 An enum type for leptons. More...
 
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  meson { P_0, P_P, K_0, K_P, D_0, B_D, B_P, B_S, PHI, K_star, MESON_END }
 An enum type for mesons. More...
 
enum  quark { UP, DOWN, CHARM, STRANGE, TOP, BOTTOM }
 An enum type for quarks. More...
 

Constructor & Destructor Documentation

NPEffectiveBS::NPEffectiveBS ( const bool  FlagLeptonUniversal_in = false,
const bool  FlagQuarkUniversal_in = false 
)

Constructor.

Parameters
[in]FlagLeptonUniversal_inan internal boolean flag that is true if assuming lepton flavour universality
[in]FlagQuarkUniversal_inan internal boolean flag that is true if assuming quark flavour universality

Definition at line 31 of file NPEffectiveBS.cpp.

32 : NPbase(), FlagLeptonUniversal(FlagLeptonUniversal_in), FlagQuarkUniversal(FlagQuarkUniversal_in)
33 {
34  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cWB_NP", boost::cref(cWB)));
35  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cH_NP", boost::cref(cH)));
36  if (FlagLeptonUniversal) {
37  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cLL_NP", boost::cref(cL1L1)));
38  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHLp_NP", boost::cref(cHL1p)));
39  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHL_NP", boost::cref(cHL1)));
40  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHE_NP", boost::cref(cHE1)));
41  } else {
42  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cL1L1_NP", boost::cref(cL1L1)));
43  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cL1L2_NP", boost::cref(cL1L2)));
44  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cL1L3_NP", boost::cref(cL1L3)));
45  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cL2L2_NP", boost::cref(cL2L2)));
46  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cL2L3_NP", boost::cref(cL2L3)));
47  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cL3L3_NP", boost::cref(cL3L3)));
48  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHL1p_NP", boost::cref(cHL1p)));
49  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHL2p_NP", boost::cref(cHL2p)));
50  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHL3p_NP", boost::cref(cHL3p)));
51  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHL1_NP", boost::cref(cHL1)));
52  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHL2_NP", boost::cref(cHL2)));
53  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHL3_NP", boost::cref(cHL3)));
54  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHE1_NP", boost::cref(cHE1)));
55  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHE2_NP", boost::cref(cHE2)));
56  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHE3_NP", boost::cref(cHE3)));
57  }
58  if (FlagQuarkUniversal) {
59  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHQp_NP", boost::cref(cHQ1p)));
60  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHQ_NP", boost::cref(cHQ1)));
61  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHU_NP", boost::cref(cHU1)));
62  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHD_NP", boost::cref(cHD1)));
63  } else {
64  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHQ1p_NP", boost::cref(cHQ1p)));
65  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHQ2p_NP", boost::cref(cHQ2p)));
66  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHQ3p_NP", boost::cref(cHQ3p)));
67  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHQ1_NP", boost::cref(cHQ1)));
68  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHQ2_NP", boost::cref(cHQ2)));
69  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHQ3_NP", boost::cref(cHQ3)));
70  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHU1_NP", boost::cref(cHU1)));
71  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHU2_NP", boost::cref(cHU2)));
72  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHU3_NP", boost::cref(cHU3)));
73  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHD1_NP", boost::cref(cHD1)));
74  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHD2_NP", boost::cref(cHD2)));
75  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHD3_NP", boost::cref(cHD3)));
76  }
77  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("Lambda_NP", boost::cref(Lambda_NP)));
78 }
NPbase()
The default constructor.
Definition: NPbase.cpp:10
double cHL3
The dimension-6 operator coefficient .
const bool FlagQuarkUniversal
An internal boolean flag that is true if assuming quark flavour universality.
double cHQ3p
The dimension-6 operator coefficient .
double cHL3p
The dimension-6 operator coefficient .
double cHQ2p
The dimension-6 operator coefficient .
double cL2L2
The dimension-6 operator coefficient .
double cHQ1
The dimension-6 operator coefficient .
const bool FlagLeptonUniversal
An internal boolean flag that is true if assuming lepton flavour universality.
double cWB
The dimension-6 operator coefficient .
double cHD1
The dimension-6 operator coefficient .
double cL2L3
The dimension-6 operator coefficient .
double cHU2
The dimension-6 operator coefficient .
double cL1L2
The dimension-6 operator coefficient .
double cHQ2
The dimension-6 operator coefficient .
double cL3L3
The dimension-6 operator coefficient .
double cHQ3
The dimension-6 operator coefficient .
double cHE1
The dimension-6 operator coefficient .
double cL1L1
The dimension-6 operator coefficient .
double Lambda_NP
The new physics scale .
double cHQ1p
The dimension-6 operator coefficient .
double cHU1
The dimension-6 operator coefficient .
double cL1L3
The dimension-6 operator coefficient .
double cHE3
The dimension-6 operator coefficient .
double cHE2
The dimension-6 operator coefficient .
double cHD3
The dimension-6 operator coefficient .
double cHL1
The dimension-6 operator coefficient .
double cHD2
The dimension-6 operator coefficient .
double cHL1p
The dimension-6 operator coefficient .
double cHU3
The dimension-6 operator coefficient .
std::map< std::string, boost::reference_wrapper< const double > > ModelParamMap
Definition: Model.h:200
double cHL2p
The dimension-6 operator coefficient .
double cHL2
The dimension-6 operator coefficient .
double cH
The dimension-6 operator coefficient .

Member Function Documentation

bool NPEffectiveBS::CheckParameters ( const std::map< std::string, double > &  DPars)
virtual

A method to check if all the mandatory parameters for the model have been provided in model initialization.

Parameters
[in]DParsa map of the parameters that are being updated in the Monte Carlo run (including parameters that are varied and those that are held constant)
Returns
a boolean that is true if the execution is successful

Reimplemented from StandardModel.

Definition at line 184 of file NPEffectiveBS.cpp.

185 {
187  for (int i = 0; i < NNPEffectiveBSVars_LFU_QFU; i++) {
188  if (DPars.find(NPEffectiveBSVars_LFU_QFU[i]) == DPars.end()) {
189  std::cout << "ERROR: Missing mandatory NPEffectiveBS_LFU_QFU parameter "
190  << NPEffectiveBSVars_LFU_QFU[i] << std::endl;
191  return false;
192  }
193  }
194  } else if (FlagLeptonUniversal && !FlagQuarkUniversal) {
195  for (int i = 0; i < NNPEffectiveBSVars_LFU; i++) {
196  if (DPars.find(NPEffectiveBSVars_LFU[i]) == DPars.end()) {
197  std::cout << "ERROR: Missing mandatory NPEffectiveBS_LFU parameter "
198  << NPEffectiveBSVars_LFU[i] << std::endl;
199  return false;
200  }
201  }
202  } else if (!FlagLeptonUniversal && FlagQuarkUniversal) {
203  for (int i = 0; i < NNPEffectiveBSVars_QFU; i++) {
204  if (DPars.find(NPEffectiveBSVars_QFU[i]) == DPars.end()) {
205  std::cout << "ERROR: Missing mandatory NPEffectiveBS_QFU parameter "
206  << NPEffectiveBSVars_QFU[i] << std::endl;
207  return false;
208  }
209  }
210  } else if (!FlagLeptonUniversal && !FlagQuarkUniversal) {
211  for (int i = 0; i < NNPEffectiveBSVars; i++) {
212  if (DPars.find(NPEffectiveBSVars[i]) == DPars.end()) {
213  std::cout << "ERROR: Missing mandatory NPEffectiveBS parameter "
214  << NPEffectiveBSVars[i] << std::endl;
215  return false;
216  }
217  }
218  } else
219  throw std::runtime_error("Error in NPEffectiveBS::CheckParameters()");
220 
221  return (NPbase::CheckParameters(DPars));
222 }
static const int NNPEffectiveBSVars
The number of the model parameters in NPEffectiveBS.
const bool FlagQuarkUniversal
An internal boolean flag that is true if assuming quark flavour universality.
static const int NNPEffectiveBSVars_QFU
The number of the model parameters in NPEffectiveBS with quark flavour universality.
static const std::string NPEffectiveBSVars[NNPEffectiveBSVars]
A string array containing the labels of the model parameters in NPEffectiveBS.
static const std::string NPEffectiveBSVars_QFU[NNPEffectiveBSVars_QFU]
A string array containing the labels of the model parameters in NPEffectiveBS with quark flavour univ...
const bool FlagLeptonUniversal
An internal boolean flag that is true if assuming lepton flavour universality.
static const std::string NPEffectiveBSVars_LFU[NNPEffectiveBSVars_LFU]
A string array containing the labels of the model parameters in NPEffectiveBS with lepton flavour uni...
static const std::string NPEffectiveBSVars_LFU_QFU[NNPEffectiveBSVars_LFU_QFU]
A string array containing the labels of the model parameters in NPEffectiveBS with lepton and quark f...
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 initi...
static const int NNPEffectiveBSVars_LFU_QFU
The number of the model parameters in NPEffectiveBS with lepton and quark flavour universalities...
static const int NNPEffectiveBSVars_LFU
The number of the model parameters in NPEffectiveBS with lepton flavour universality.
double NPEffectiveBS::deltaGA_f ( const Particle  f) const
virtual

New physics contribution to the neutral-current axial-vector coupling \(g_A^f\).

New physics contribution to the neutral-current axial-vector coupling \(g_A^f\) is given by

\[ \delta g_A^f = \delta g_L^f - \delta g_R^f. \]

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

Reimplemented from NPbase.

Definition at line 304 of file NPEffectiveBS.cpp.

305 {
306  return ( deltaGL_f_tmp(f) - deltaGR_f_tmp(f) + NPbase::deltaGA_f(f));
307 }
double deltaGR_f_tmp(const Particle f) const
New physics contribution to .
double deltaGL_f_tmp(const Particle f) const
New physics contribution to .
virtual double deltaGA_f(const Particle f) const
New physics contribution to the neutral-current axial-vector coupling .
Definition: NPbase.cpp:82
double NPEffectiveBS::DeltaGF ( ) const
virtual

New physics contribution to the Fermi constant.

The new physics contribution \(\Delta G\) is defined as

\[ G_\mu = G_{\mu,\mathrm{SM}}(1+\Delta G)\,, \]

where \(G_\mu\) is the experimentl value measured through muon decays, and \(G_{\mu,\mathrm{SM}}\) is the Fermi constant in the SM. The dimension-six operators yield the following corrections:

\[ \Delta G = - C_{LL}\left(\frac{v}{\Lambda}\right)^2 + C^\prime_{HL_1}\left(\frac{v}{\Lambda}\right)^2 + C^\prime_{HL_2}\left(\frac{v}{\Lambda}\right)^2. \]

See [58] and references therein.

Returns
\(\Delta G\)

Reimplemented from NPbase.

Definition at line 244 of file NPEffectiveBS.cpp.

245 {
246  double ratio = v() * v() / Lambda_NP / Lambda_NP;
247 
248  return ( -(cL1L2 - cHL1p - cHL2p) * ratio);
249 }
virtual double v() const
The SM Higgs vacuum expectation value .
double cL1L2
The dimension-6 operator coefficient .
double Lambda_NP
The new physics scale .
double cHL1p
The dimension-6 operator coefficient .
double cHL2p
The dimension-6 operator coefficient .
double NPEffectiveBS::deltaGL_f_tmp ( const Particle  f) const
private

New physics contribution to \(g_L^f\).

New physics contributions to the neutral-current left-handed coupling \(g_L^f\): e.g.,

\[ \delta g_L^{\nu_i} = \frac{C_{HL_i}^\prime-C_{HL_i}}{2} \left(\frac{v}{\Lambda}\right)^2, \qquad \delta g_L^{e_i} = -\frac{C_{HL_i}^\prime+C_{HL_i}}{2} \left(\frac{v}{\Lambda}\right)^2, \]

where \(L\to Q\), \(\nu_i\to u_i\) and \(e_i\to d_i\) for quarks. See [58] and references therein.

Parameters
[in]fa lepton or quark
Returns
\(\delta g_L^f\)
Attention
The new physics contribution via \(S\), \(T\), \(U\) and \(\Delta G\) are not included in this function.

Definition at line 312 of file NPEffectiveBS.cpp.

313 {
314  double ratio = v() * v() / Lambda_NP / Lambda_NP;
315  if (f.is("NEUTRINO_1"))
316  return ( (cHL1p - cHL1) / 2.0 * ratio);
317  if (f.is("NEUTRINO_2"))
318  return ( (cHL2p - cHL2) / 2.0 * ratio);
319  if (f.is("NEUTRINO_3"))
320  return ( (cHL3p - cHL3) / 2.0 * ratio);
321  if (f.is("ELECTRON"))
322  return ( -(cHL1p + cHL1) / 2.0 * ratio);
323  if (f.is("MU"))
324  return ( -(cHL2p + cHL2) / 2.0 * ratio);
325  if (f.is("TAU"))
326  return ( -(cHL3p + cHL3) / 2.0 * ratio);
327  if (f.is("UP"))
328  return ( (cHQ1p - cHQ1) / 2.0 * ratio);
329  if (f.is("CHARM"))
330  return ( (cHQ2p - cHQ2) / 2.0 * ratio);
331  if (f.is("TOP"))
332  return 0.0;
333  if (f.is("DOWN"))
334  return ( -(cHQ1p + cHQ1) / 2.0 * ratio);
335  if (f.is("STRANGE"))
336  return ( -(cHQ2p + cHQ2) / 2.0 * ratio);
337  if (f.is("BOTTOM"))
338  return ( -(cHQ3p + cHQ3) / 2.0 * ratio);
339  throw std::runtime_error("Error in NPEffectiveBS::deltaGL_f()");
340 }
double cHL3
The dimension-6 operator coefficient .
double cHQ3p
The dimension-6 operator coefficient .
double cHL3p
The dimension-6 operator coefficient .
double cHQ2p
The dimension-6 operator coefficient .
double cHQ1
The dimension-6 operator coefficient .
bool is(std::string name_i) const
Definition: Particle.cpp:23
virtual double v() const
The SM Higgs vacuum expectation value .
double cHQ2
The dimension-6 operator coefficient .
double cHQ3
The dimension-6 operator coefficient .
double Lambda_NP
The new physics scale .
double cHQ1p
The dimension-6 operator coefficient .
double cHL1
The dimension-6 operator coefficient .
double cHL1p
The dimension-6 operator coefficient .
double cHL2p
The dimension-6 operator coefficient .
double cHL2
The dimension-6 operator coefficient .
double NPEffectiveBS::deltaGR_f_tmp ( const Particle  f) const
private

New physics contribution to \(g_R^f\).

New physics contributions to the neutral-current right-handed coupling \(g_R^f\): e.g.,

\[ \delta g_R^{e_i} = -\frac{C_{HE_i}}{2} \left(\frac{v}{\Lambda}\right)^2. \]

See [58] and references therein.

Parameters
[in]fa lepton or quark
Returns
\(\delta g_R^f\)
Attention
The new physics contribution via \(S\), \(T\), \(U\) and \(\Delta G\) are not included in this function.

Definition at line 342 of file NPEffectiveBS.cpp.

343 {
344  double ratio = v() * v() / Lambda_NP / Lambda_NP;
345  if (f.is("NEUTRINO_1"))
346  return ( 0.);
347  if (f.is("NEUTRINO_2"))
348  return ( 0.);
349  if (f.is("NEUTRINO_3"))
350  return ( 0.);
351  if (f.is("ELECTRON"))
352  return ( -cHE1 / 2.0 * ratio);
353  if (f.is("MU"))
354  return ( -cHE2 / 2.0 * ratio);
355  if (f.is("TAU"))
356  return ( -cHE3 / 2.0 * ratio);
357  if (f.is("UP"))
358  return ( -cHU1 / 2.0 * ratio);
359  if (f.is("CHARM"))
360  return ( -cHU2 / 2.0 * ratio);
361  if (f.is("TOP"))
362  return 0.0;
363  if (f.is("DOWN"))
364  return ( -cHD1 / 2.0 * ratio);
365  if (f.is("STRANGE"))
366  return ( -cHD2 / 2.0 * ratio);
367  if (f.is("BOTTOM"))
368  return ( -cHD3 / 2.0 * ratio);
369  throw std::runtime_error("Error in NPEffectiveBS::deltaGL_f()");
370 }
bool is(std::string name_i) const
Definition: Particle.cpp:23
double cHD1
The dimension-6 operator coefficient .
double cHU2
The dimension-6 operator coefficient .
virtual double v() const
The SM Higgs vacuum expectation value .
double cHE1
The dimension-6 operator coefficient .
double Lambda_NP
The new physics scale .
double cHU1
The dimension-6 operator coefficient .
double cHE3
The dimension-6 operator coefficient .
double cHE2
The dimension-6 operator coefficient .
double cHD3
The dimension-6 operator coefficient .
double cHD2
The dimension-6 operator coefficient .
double NPEffectiveBS::deltaGV_f ( const Particle  f) const
virtual

New physics contribution to the neutral-current vector coupling \(g_V^f\).

New physics contribution to the neutral-current vector coupling \(g_V^f\) is given by

\[ \delta g_V^f = \delta g_L^f + \delta g_R^f. \]

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

Reimplemented from NPbase.

Definition at line 299 of file NPEffectiveBS.cpp.

300 {
301  return ( deltaGL_f_tmp(f) + deltaGR_f_tmp(f) + NPbase::deltaGV_f(f));
302 }
virtual double deltaGV_f(const Particle f) const
New physics contribution to the neutral-current vector coupling .
Definition: NPbase.cpp:60
double deltaGR_f_tmp(const Particle f) const
New physics contribution to .
double deltaGL_f_tmp(const Particle f) const
New physics contribution to .
double NPEffectiveBS::GammaW ( ) const
virtual

The total width of the \(W\) boson, \(\Gamma_W\).

The width of the \(W\) boson receives the corrections from the operators \(\mathcal{O}^\prime_{HL}\) and \(\mathcal{O}^\prime_{HQ}\), in addition to the corrections via \(S\), \(T\), \(U\) and \(\Delta G\):

\[ \Gamma_W = \Gamma_{W,\mathrm{SM}} \left[ 1 - \frac{3\alpha(M_Z^2)}{4(c_W^2-s_W^2)} \left( S - 2c_W^2\,T - \frac{c_W^2-s_W^2}{2s_W^2}\,U \right) - \frac{1+c_W^2}{2(c_W^2-s_W^2)}\, \Delta G + \left( C_{HL_1}^\prime + C_{HL_2}^\prime + C_{HL_3}^\prime + C_{HQ_1}^\prime + C_{HQ_2}^\prime\right) \left(\frac{v}{\Lambda}\right)^2 \right]. \]

See [58] and references therein.

Returns
the total width of the \(W\) boson in GeV

Reimplemented from NPbase.

Definition at line 254 of file NPEffectiveBS.cpp.

255 {
256  double Gamma_W = trueSM.GammaW();
257 
258  double alpha = alphaMz();
259  double c2 = trueSM.cW2();
260  double s2 = trueSM.sW2();
261  double ratio = v() * v() / Lambda_NP / Lambda_NP;
262 
263  Gamma_W *= 1.0 - 3.0 * alpha / 4.0 / (c2 - s2)
264  *(obliqueS() - 2.0 * c2 * obliqueT() - (c2 - s2) * obliqueU() / 2.0 / s2)
265  - (1.0 + c2) / 2.0 / (c2 - s2) * DeltaGF()
266  //+ (cHL1p + cHL2p + cHL3p + cHQ1p + cHQ2p) * ratio; // incorrect
267  + 2.0 / 9.0 * (cHL1p + cHL2p + cHL3p + 3.0 * cHQ1p + 3.0 * cHQ2p) * ratio;
268 
269  return Gamma_W;
270 }
virtual double sW2(const double Mw_i) const
The square of the sine of the weak mixing angle in the on-shell scheme, denoted as ...
virtual double DeltaGF() const
New physics contribution to the Fermi constant.
double cHL3p
The dimension-6 operator coefficient .
StandardModel trueSM
Definition: NPbase.h:543
virtual double obliqueU() const
The oblique parameter .
double cHQ2p
The dimension-6 operator coefficient .
virtual double obliqueT() const
The oblique parameter .
virtual double GammaW(const Particle fi, const Particle fj) const
A partial decay width of the boson decay into a SM fermion pair.
virtual double v() const
The SM Higgs vacuum expectation value .
double Lambda_NP
The new physics scale .
double cHQ1p
The dimension-6 operator coefficient .
virtual double obliqueS() const
The oblique parameter .
virtual double cW2(const double Mw_i) const
The square of the cosine of the weak mixing angle in the on-shell scheme, denoted as ...
double alphaMz() const
The electromagnetic coupling at the -mass scale, .
double cHL1p
The dimension-6 operator coefficient .
double cHL2p
The dimension-6 operator coefficient .
double NPEffectiveBS::Mw_tree ( ) const
virtual

The tree-level mass of the \(W\) boson, \(M_W^{\mathrm{tree}}\).

The tree-level mass is obtained by subtracting the effect of \(\Delta G\).

Returns
\(M_W^{\mathrm{tree}}\) in GeV.

Reimplemented from StandardModel.

Definition at line 237 of file NPEffectiveBS.cpp.

238 {
239  double GF0 = GF * (1.0 - DeltaGF());
240  double tmp = 4.0 * M_PI * ale / sqrt(2.0) / GF0 / Mz / Mz;
241  return ( Mz / sqrt(2.0) * sqrt(1.0 + sqrt(1.0 - tmp)));
242 }
virtual double DeltaGF() const
New physics contribution to the Fermi constant.
double ale
The fine-structure constant .
double GF
The Fermi constant in .
double Mz
The mass of the boson in GeV.
complex sqrt(const complex &z)
double NPEffectiveBS::obliqueS ( ) const
virtual

The oblique parameter \(S\).

The operator \(\mathcal{O}_{WB}\) contributes to the \(S\) parameter:

\[ S = \frac{4s_Wc_W\, C_{WB} }{\alpha(M_Z^2)} \left(\frac{v}{\Lambda}\right)^2. \]

See [58] and references therein.

Returns
the value of \(S\)

Reimplemented from NPbase.

Definition at line 275 of file NPEffectiveBS.cpp.

276 {
277  double ratio = v() * v() / Lambda_NP / Lambda_NP;
278  double sW_SM = sqrt(trueSM.sW2()); /* This has to be the SM value. */
279  double cW_SM = sqrt(trueSM.cW2()); /* This has to be the SM value. */
280 
281  return ( 4.0 * sW_SM * cW_SM * cWB / alphaMz() * ratio);
282 }
virtual double sW2(const double Mw_i) const
The square of the sine of the weak mixing angle in the on-shell scheme, denoted as ...
StandardModel trueSM
Definition: NPbase.h:543
double cWB
The dimension-6 operator coefficient .
virtual double v() const
The SM Higgs vacuum expectation value .
double Lambda_NP
The new physics scale .
virtual double cW2(const double Mw_i) const
The square of the cosine of the weak mixing angle in the on-shell scheme, denoted as ...
double alphaMz() const
The electromagnetic coupling at the -mass scale, .
complex sqrt(const complex &z)
double NPEffectiveBS::obliqueT ( ) const
virtual

The oblique parameter \(T\).

The operator \(\mathcal{O}_{H}\) contributes to the \(T\) parameter:

\[ T = - \frac{C_H}{2\alpha(M_Z^2)}\left(\frac{v}{\Lambda}\right)^2. \]

See [58] and references therein.

Returns
the value of \(T\)

Reimplemented from NPbase.

Definition at line 284 of file NPEffectiveBS.cpp.

285 {
286  double ratio = v() * v() / Lambda_NP / Lambda_NP;
287 
288  return ( -cH / 2.0 / alphaMz() * ratio);
289 }
virtual double v() const
The SM Higgs vacuum expectation value .
double Lambda_NP
The new physics scale .
double alphaMz() const
The electromagnetic coupling at the -mass scale, .
double cH
The dimension-6 operator coefficient .
double NPEffectiveBS::obliqueU ( ) const
virtual

The oblique parameter \(U\).

Returns
\(U=0\)

Reimplemented from NPbase.

Definition at line 291 of file NPEffectiveBS.cpp.

292 {
293  return 0.0;
294 }
void NPEffectiveBS::setParameter ( const std::string  name,
const double &  value 
)
protectedvirtual

A method to set the value of a parameter of the model.

Parameters
[in]namename of a model parameter
[in]valuethe value to be assigned to the parameter specified by name

Reimplemented from StandardModel.

Definition at line 83 of file NPEffectiveBS.cpp.

84 {
85  if (name.compare("cWB_NP") == 0)
86  cWB = value;
87  else if (name.compare("cH_NP") == 0)
88  cH = value;
89  else if (name.compare("cL1L1_NP") == 0)
90  cL1L1 = value;
91  else if (name.compare("cL1L2_NP") == 0)
92  cL1L2 = value;
93  else if (name.compare("cL1L3_NP") == 0)
94  cL1L3 = value;
95  else if (name.compare("cL2L2_NP") == 0)
96  cL2L2 = value;
97  else if (name.compare("cL2L3_NP") == 0)
98  cL2L3 = value;
99  else if (name.compare("cL3L3_NP") == 0)
100  cL3L3 = value;
101  else if (name.compare("cLL_NP") == 0) {
102  cL1L1 = value;
103  cL1L2 = value;
104  cL1L3 = value;
105  cL2L2 = value;
106  cL2L3 = value;
107  cL3L3 = value;
108  } else if (name.compare("cHL1p_NP") == 0)
109  cHL1p = value;
110  else if (name.compare("cHL2p_NP") == 0)
111  cHL2p = value;
112  else if (name.compare("cHL3p_NP") == 0)
113  cHL3p = value;
114  else if (name.compare("cHLp_NP") == 0) {
115  cHL1p = value;
116  cHL2p = value;
117  cHL3p = value;
118  } else if (name.compare("cHQ1p_NP") == 0)
119  cHQ1p = value;
120  else if (name.compare("cHQ2p_NP") == 0)
121  cHQ2p = value;
122  else if (name.compare("cHQ3p_NP") == 0)
123  cHQ3p = value;
124  else if (name.compare("cHQp_NP") == 0) {
125  cHQ1p = value;
126  cHQ2p = value;
127  cHQ3p = value;
128  } else if (name.compare("cHL1_NP") == 0)
129  cHL1 = value;
130  else if (name.compare("cHL2_NP") == 0)
131  cHL2 = value;
132  else if (name.compare("cHL3_NP") == 0)
133  cHL3 = value;
134  else if (name.compare("cHL_NP") == 0) {
135  cHL1 = value;
136  cHL2 = value;
137  cHL3 = value;
138  } else if (name.compare("cHQ1_NP") == 0)
139  cHQ1 = value;
140  else if (name.compare("cHQ2_NP") == 0)
141  cHQ2 = value;
142  else if (name.compare("cHQ3_NP") == 0)
143  cHQ3 = value;
144  else if (name.compare("cHQ_NP") == 0) {
145  cHQ1 = value;
146  cHQ2 = value;
147  cHQ3 = value;
148  } else if (name.compare("cHE1_NP") == 0)
149  cHE1 = value;
150  else if (name.compare("cHE2_NP") == 0)
151  cHE2 = value;
152  else if (name.compare("cHE3_NP") == 0)
153  cHE3 = value;
154  else if (name.compare("cHE_NP") == 0) {
155  cHE1 = value;
156  cHE2 = value;
157  cHE3 = value;
158  } else if (name.compare("cHU1_NP") == 0)
159  cHU1 = value;
160  else if (name.compare("cHU2_NP") == 0)
161  cHU2 = value;
162  else if (name.compare("cHU3_NP") == 0)
163  cHU3 = value;
164  else if (name.compare("cHU_NP") == 0) {
165  cHU1 = value;
166  cHU2 = value;
167  cHU3 = value;
168  } else if (name.compare("cHD1_NP") == 0)
169  cHD1 = value;
170  else if (name.compare("cHD2_NP") == 0)
171  cHD2 = value;
172  else if (name.compare("cHD3_NP") == 0)
173  cHD3 = value;
174  else if (name.compare("cHD_NP") == 0) {
175  cHD1 = value;
176  cHD2 = value;
177  cHD3 = value;
178  } else if (name.compare("Lambda_NP") == 0)
179  Lambda_NP = value;
180  else
181  NPbase::setParameter(name, value);
182 }
double cHL3
The dimension-6 operator coefficient .
double cHQ3p
The dimension-6 operator coefficient .
double cHL3p
The dimension-6 operator coefficient .
double cHQ2p
The dimension-6 operator coefficient .
double cL2L2
The dimension-6 operator coefficient .
std::string name
The name of the model.
Definition: Model.h:203
double cHQ1
The dimension-6 operator coefficient .
double cWB
The dimension-6 operator coefficient .
virtual void setParameter(const std::string name, const double &value)
A method to set the value of a parameter of StandardModel.
double cHD1
The dimension-6 operator coefficient .
double cL2L3
The dimension-6 operator coefficient .
double cHU2
The dimension-6 operator coefficient .
double cL1L2
The dimension-6 operator coefficient .
double cHQ2
The dimension-6 operator coefficient .
double cL3L3
The dimension-6 operator coefficient .
double cHQ3
The dimension-6 operator coefficient .
double cHE1
The dimension-6 operator coefficient .
double cL1L1
The dimension-6 operator coefficient .
double Lambda_NP
The new physics scale .
double cHQ1p
The dimension-6 operator coefficient .
double cHU1
The dimension-6 operator coefficient .
double cL1L3
The dimension-6 operator coefficient .
double cHE3
The dimension-6 operator coefficient .
double cHE2
The dimension-6 operator coefficient .
double cHD3
The dimension-6 operator coefficient .
double cHL1
The dimension-6 operator coefficient .
double cHD2
The dimension-6 operator coefficient .
double cHL1p
The dimension-6 operator coefficient .
double cHU3
The dimension-6 operator coefficient .
double cHL2p
The dimension-6 operator coefficient .
double cHL2
The dimension-6 operator coefficient .
double cH
The dimension-6 operator coefficient .
double NPEffectiveBS::v ( ) const
virtual

The SM Higgs vacuum expectation value \(v\).

\[ v = \left(\frac{1}{\sqrt{2} G_F}\right)^{1/2}. \]

Returns
the SM value of the Higgs VEV in GeV

Reimplemented from StandardModel.

Definition at line 228 of file NPEffectiveBS.cpp.

229 {
230  //return ( sqrt( (1.0 - (cL1L2 - cHL1p - cHL2p)/sqrt(2.0)/GF/LambdaNP/LambdaNP)
231  // /sqrt(2.0)/GF ) );
232 
233  /* use the tree-level relation */
234  return trueSM.v();
235 }
StandardModel trueSM
Definition: NPbase.h:543
virtual double v() const
The Higgs vacuum expectation value. where is the Fermi constant, measured through muon decays...

Member Data Documentation

double NPEffectiveBS::cH
protected

The dimension-6 operator coefficient \(C_{H}\).

Definition at line 474 of file NPEffectiveBS.h.

double NPEffectiveBS::cHD1
protected

The dimension-6 operator coefficient \(C_{HD_1}\).

Definition at line 499 of file NPEffectiveBS.h.

double NPEffectiveBS::cHD2
protected

The dimension-6 operator coefficient \(C_{HD_2}\).

Definition at line 500 of file NPEffectiveBS.h.

double NPEffectiveBS::cHD3
protected

The dimension-6 operator coefficient \(C_{HD_3}\).

Definition at line 501 of file NPEffectiveBS.h.

double NPEffectiveBS::cHE1
protected

The dimension-6 operator coefficient \(C_{HE_1}\).

Definition at line 493 of file NPEffectiveBS.h.

double NPEffectiveBS::cHE2
protected

The dimension-6 operator coefficient \(C_{HE_2}\).

Definition at line 494 of file NPEffectiveBS.h.

double NPEffectiveBS::cHE3
protected

The dimension-6 operator coefficient \(C_{HE_3}\).

Definition at line 495 of file NPEffectiveBS.h.

double NPEffectiveBS::cHL1
protected

The dimension-6 operator coefficient \(C_{HL_1}\).

Definition at line 487 of file NPEffectiveBS.h.

double NPEffectiveBS::cHL1p
protected

The dimension-6 operator coefficient \(C_{HL_1}^\prime\).

Definition at line 481 of file NPEffectiveBS.h.

double NPEffectiveBS::cHL2
protected

The dimension-6 operator coefficient \(C_{HL_2}\).

Definition at line 488 of file NPEffectiveBS.h.

double NPEffectiveBS::cHL2p
protected

The dimension-6 operator coefficient \(C_{HL_2}^\prime\).

Definition at line 482 of file NPEffectiveBS.h.

double NPEffectiveBS::cHL3
protected

The dimension-6 operator coefficient \(C_{HL_3}\).

Definition at line 489 of file NPEffectiveBS.h.

double NPEffectiveBS::cHL3p
protected

The dimension-6 operator coefficient \(C_{HL_3}^\prime\).

Definition at line 483 of file NPEffectiveBS.h.

double NPEffectiveBS::cHQ1
protected

The dimension-6 operator coefficient \(C_{HQ_1}\).

Definition at line 490 of file NPEffectiveBS.h.

double NPEffectiveBS::cHQ1p
protected

The dimension-6 operator coefficient \(C_{HQ_1}^\prime\).

Definition at line 484 of file NPEffectiveBS.h.

double NPEffectiveBS::cHQ2
protected

The dimension-6 operator coefficient \(C_{HQ_2}\).

Definition at line 491 of file NPEffectiveBS.h.

double NPEffectiveBS::cHQ2p
protected

The dimension-6 operator coefficient \(C_{HQ_2}^\prime\).

Definition at line 485 of file NPEffectiveBS.h.

double NPEffectiveBS::cHQ3
protected

The dimension-6 operator coefficient \(C_{HQ_3}\).

Definition at line 492 of file NPEffectiveBS.h.

double NPEffectiveBS::cHQ3p
protected

The dimension-6 operator coefficient \(C_{HQ_3}^\prime\).

Definition at line 486 of file NPEffectiveBS.h.

double NPEffectiveBS::cHU1
protected

The dimension-6 operator coefficient \(C_{HU_1}\).

Definition at line 496 of file NPEffectiveBS.h.

double NPEffectiveBS::cHU2
protected

The dimension-6 operator coefficient \(C_{HU_2}\).

Definition at line 497 of file NPEffectiveBS.h.

double NPEffectiveBS::cHU3
protected

The dimension-6 operator coefficient \(C_{HU_3}\).

Definition at line 498 of file NPEffectiveBS.h.

double NPEffectiveBS::cL1L1
protected

The dimension-6 operator coefficient \(C_{L_1L_1}\).

Definition at line 475 of file NPEffectiveBS.h.

double NPEffectiveBS::cL1L2
protected

The dimension-6 operator coefficient \(C_{L_1L_2}\).

Definition at line 476 of file NPEffectiveBS.h.

double NPEffectiveBS::cL1L3
protected

The dimension-6 operator coefficient \(C_{L_1L_3}\).

Definition at line 477 of file NPEffectiveBS.h.

double NPEffectiveBS::cL2L2
protected

The dimension-6 operator coefficient \(C_{L_2L_2}\).

Definition at line 478 of file NPEffectiveBS.h.

double NPEffectiveBS::cL2L3
protected

The dimension-6 operator coefficient \(C_{L_2L_3}\).

Definition at line 479 of file NPEffectiveBS.h.

double NPEffectiveBS::cL3L3
protected

The dimension-6 operator coefficient \(C_{L_3L_3}\).

Definition at line 480 of file NPEffectiveBS.h.

double NPEffectiveBS::cWB
protected

The dimension-6 operator coefficient \(C_{WB}\).

Definition at line 473 of file NPEffectiveBS.h.

const bool NPEffectiveBS::FlagLeptonUniversal
private

An internal boolean flag that is true if assuming lepton flavour universality.

Definition at line 511 of file NPEffectiveBS.h.

const bool NPEffectiveBS::FlagQuarkUniversal
private

An internal boolean flag that is true if assuming quark flavour universality.

Definition at line 517 of file NPEffectiveBS.h.

double NPEffectiveBS::Lambda_NP
protected

The new physics scale \(\Lambda\).

Definition at line 502 of file NPEffectiveBS.h.

const int NPEffectiveBS::NNPEffectiveBSVars = 30
static

The number of the model parameters in NPEffectiveBS.

 

Definition at line 274 of file NPEffectiveBS.h.

const int NPEffectiveBS::NNPEffectiveBSVars_LFU = 19
static

The number of the model parameters in NPEffectiveBS with lepton flavour universality.

 

Definition at line 280 of file NPEffectiveBS.h.

const int NPEffectiveBS::NNPEffectiveBSVars_LFU_QFU = 11
static

The number of the model parameters in NPEffectiveBS with lepton and quark flavour universalities.

 

Definition at line 292 of file NPEffectiveBS.h.

const int NPEffectiveBS::NNPEffectiveBSVars_QFU = 22
static

The number of the model parameters in NPEffectiveBS with quark flavour universality.

 

Definition at line 286 of file NPEffectiveBS.h.

const std::string NPEffectiveBS::NPEffectiveBSVars
static
Initial value:
= {"cWB", "cH", "cL1L1", "cL1L2", "cL1L3", "cL2L2", "cL2L3", "cL3L3",
"cHL1p", "cHL2p", "cHL3p", "cHQ1p", "cHQ2p", "cHQ3p",
"cHL1", "cHL2", "cHL3", "cHQ1", "cHQ2", "cHQ3", "cHE1", "cHE2", "cHE3",
"cHU1", "cHU2", "cHU3", "cHD1", "cHD2", "cHD3", "Lambda_NP"}

A string array containing the labels of the model parameters in NPEffectiveBS.

Definition at line 298 of file NPEffectiveBS.h.

const std::string NPEffectiveBS::NPEffectiveBSVars_LFU
static
Initial value:
= {"cWB_NP", "cH_NP", "cLL_NP", "cHLp_NP", "cHL_NP",
"cHQ1p_NP", "cHQ2p_NP", "cHQ3p_NP", "cHQ1_NP", "cHQ2_NP", "cHQ3_NP",
"cHU1_NP", "cHU2_NP", "cHU3_NP", "cHD1_NP", "cHD2_NP", "cHD3_NP", "cHE_NP", "Lambda_NP"}

A string array containing the labels of the model parameters in NPEffectiveBS with lepton flavour universality.

Definition at line 304 of file NPEffectiveBS.h.

const std::string NPEffectiveBS::NPEffectiveBSVars_LFU_QFU
static
Initial value:
= {"cWB_NP", "cH_NP", "cLL_NP", "cHLp_NP", "cHQp_NP",
"cHL_NP", "cHQ_NP", "cHE_NP", "cHU_NP", "cHD_NP", "Lambda_NP"}

A string array containing the labels of the model parameters in NPEffectiveBS with lepton and quark flavour universalities.

Definition at line 316 of file NPEffectiveBS.h.

const std::string NPEffectiveBS::NPEffectiveBSVars_QFU
static
Initial value:
= {"cWB_NP", "cH_NP", "cL1L1", "cL1L2", "cL1L3", "cL2L2", "cL2L3", "cL3L3",
"cHL1p", "cHL2p", "cHL3p", "cHQp_NP", "cHL1", "cHL2", "cHL3",
"cHQ_NP", "cHE1", "cHE2", "cHE3", "cHU_NP", "cHD_NP", "Lambda_NP"}

A string array containing the labels of the model parameters in NPEffectiveBS with quark flavour universality.

Definition at line 310 of file NPEffectiveBS.h.


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