A base class for $Z_2$ symmetric Two-Higgs-Doublet models. More...

#include <THDM.h>

Inheritance diagram for THDM:

Detailed Description

A base class for $Z_2$ symmetric Two-Higgs-Doublet models.

Author: HEPfit Collaboration

Copyright: GNU General Public License

Model parameters

The model parameters of THDM are summarized below. The current implementation only allows for a Two-Higgs-Doublet model with a softly broken $Z_2$ symmetry and without CP violation in the Higgs potential. The scalar 125 GeV resonance is assumed to be the one of the CP-even Higgs states of the model; the other is attributed to mHh in the THDM configuration.

Label	LaTeX symbol	Description
logtb	$\log_{10} \tan \beta$	The decadic logarithm of the tangent of the mixing angle $\beta \equiv \arctan \frac{v_2}{v_1}$.
bma	$\beta - \alpha$	The difference of the mixing angles $\beta$ and $\alpha$.
mHh2	$m_H^2$	The mass square of the "non-125 GeV" CP-even Higgs state.
mA2	$m_A^2$	The squared masses of the CP-odd Higgs $A$.
mHp2	$m_{H^+}^2$	The squared masses of the charged Higgs.
m12_2	$m_{12}^2$	The soft $Z_2$ symmetry breaking parameter.
bsgamma_theoryerror	$\theta_{b\to s\gamma}^{\text{theo}}$	A nuisance parameter between -1 and +1 for the theoretical error in the determination of $BR(B\to X_s \gamma)$.
Q_THDM	$Q_{\text{THDM}}$	The THDM scale.
Rpeps	$\varepsilon_{\text{R'}}$	Minimal value for which the R' perturbativity criterion should be applied for the unitarity bounds.
NLOuniscale	$Q_{\text{min}}^{\text{NLOuni}}$	Minimal scale at which the NLO unitarity conditions are checked.

Model flags

Label	Value	Description
modelTypeflag	typeI / typeII / typeX / typeY	This flag determines the type of $Z_2$ symmetry.
RGEorder	LO / approxNLO / NLO	This flag determines the order in perturbation theory of the renormalization group equations.
wavefunctionrenormalization	true / false	Whether to use wavefunction renormalization for NLO unitarity constraints.

Definition at line 120 of file THDM.h.

Public Member Functions
virtual bool	CheckParameters (const std::map< std::string, double > &DPars)
	A method to check if all the mandatory parameters for THDM have been provided in model initialization. More...

double	getBDstartaunu_A () const
	A method get the $ B\to D^*\tau\nu $ coefficient A. More...

double	getBDstartaunu_B () const
	A method get the $ B\to D \tau\nu$ coefficient B. More...

double	getBDstartaunu_SM () const
	A method get the $ B\to D^*\tau\nu $ SM expectation. More...

double	getBDtaunu_A () const
	A method get the $ B\to D \tau\nu$ coefficient A. More...

double	getBDtaunu_B () const
	A method get the $ B\to D \tau\nu$ coefficient B. More...

double	getBDtaunu_SM () const
	A method get the $ B\to D \tau\nu$ SM expectation. More...

double	getbma () const
	A method get $\beta-\alpha$. More...

double	getbsgamma_theoryerror () const
	A method get the nuisance parameter for the theoretical error on $ b\to s\gamma$. More...

double	getcosa () const
	A method get $\cos \alpha$. More...

double	getcosb () const
	A method get $\cos \beta$. More...

double	getlogtb () const
	A method get $\log_{10}(\tan \beta)$. More...

double	getm12_2 () const
	A method get the parameter of the Higgs potential $m_{12}^2$. More...

double	getmA () const
	A method get the mass of the pseudoscalar Higgs A. More...

double	getmA2 () const
	A method get the squared mass of the pseudoscalar Higgs A. More...

virtual THDMMatching &	getMatching () const
	A get method to access the member reference of type THDMMatching. More...

double	getmHh () const
	A method get the mass of the "non-125 GeV" neutral scalar Higgs. More...

double	getmHh2 () const
	A method get the squared mass of the "non-125 GeV" neutral scalar Higgs. More...

double	getmHl2 () const
	A method get the squared mass of the lighter neutral scalar Higgs. More...

double	getmHp () const
	A method get the charged Higgs mass. More...

double	getmHp2 () const
	A method get the squared charged Higgs mass. More...

std::string	getModelTypeflag () const
	A method get the THDM model type. More...

THDMcache *	getMyTHDMCache () const
	A method get the THDMCache. More...

double	getNLOuniscale () const
	A method get the minimal NLO unitarity check scale. More...

double	getQ_THDM () const
	A method get the THDM scale. More...

std::string	getRGEorderflag () const
	A method get switch for NLO RGE and approximate NLO RGE. More...

double	getRpeps () const
	A method get the minimal R' value. More...

double	getsin_ba () const
	A method get $\sin(\beta-\alpha)$. More...

double	getsina () const
	A method get $\sin \alpha$. More...

double	getsinb () const
	A method get $\sin \beta$. More...

bool	getsqmassesflag () const
	A method choose if you want to use the THDM masses or rather their squares. More...

double	gettanb () const
	A method get $\tan \beta$. More...

double	getv1 () const
	A method get the VEV $v_1$. More...

double	getv2 () const
	A method get the VEV $v_2$. More...

bool	getWFRflag () const
	A method get the flag to switch on wavefunction renormalization for the NLO unitarity conditions. More...

virtual bool	Init (const std::map< std::string, double > &DPars)
	Initializes the THDM parameters found in the argument. More...

virtual bool	InitializeModel ()
	A method to initialize the model. More...

virtual bool	PostUpdate ()
	The post-update method for THDM. More...

virtual bool	PreUpdate ()
	The pre-update method for THDM. More...

virtual bool	setFlag (const std::string name, const bool value)
	A method to set a flag of THDM. More...

virtual bool	setFlagStr (const std::string name, const std::string value)
	A method to set a string flag of THDM. More...

	THDM ()
	THDM constructor. More...

virtual bool	Update (const std::map< std::string, double > &DPars)
	The update method for THDM. More...

	~THDM ()
	THDM destructor. More...

Public Member Functions inherited from StandardModel
virtual double	A_f (const Particle f) const
	The left-right asymmetry in $e^+e^-\to Z\to \ell \bar{\ell}$ at the $Z$-pole, $\mathcal{A}_\ell$. More...

virtual double	AFB (const Particle f) const

double	Ale (double mu, orders order, bool Nf_thr=true) const
	The running electromagnetic coupling $\alpha_e(\mu)$ in the $\overline{MS}$ scheme. More...

double	ale_OS (const double mu, orders order=FULLNLO) const
	The running electromagnetic coupling $\alpha(\mu)$ in the on-shell scheme. More...

double	alphaMz () const
	The electromagnetic coupling at the $Z$-mass scale, $\alpha(M_Z^2)=\alpha/(1-\Delta\alpha(M_Z^2))$. More...

double	Als (double mu, orders order=FULLNLO, bool qed_flag=false, bool Nf_thr=true) const
	The running QCD coupling $\alpha(\mu)$ in the $\overline{MS}$ scheme including QED corrections. More...

double	AlsByOrder (double mu, orders order=FULLNLO, bool qed_flag=false, bool Nf_thr=true) const

double	Alstilde5 (const double mu) const
	The value of $\frac{\alpha_s^{\mathrm{FULLNLO}}}{4\pi}$ at any scale $\mu$ with the number of flavours $n_f = 4$ and full EW corrections. More...

double	Beta_e (int nm, unsigned int nf) const
	QED beta function coefficients - eq. (36) hep-ph/0512066. More...

double	Beta_s (int nm, unsigned int nf) const
	QCD beta function coefficients including QED corrections - eq. (36) hep-ph/0512066. More...

double	c02 () const
	The square of the cosine of the weak mixing angle $c_0^2$ defined without weak radiative corrections. More...

virtual bool	CheckFlags () const
	A method to check the sanity of the set of model flags. More...

bool	checkSMparamsForEWPO ()
	A method to check whether the parameters relevant to the EWPO are updated. More...

double	computeBrHtobb () const
	The Br $(H\to bb)$ in the Standard Model. More...

double	computeBrHtocc () const
	The Br $(H\to cc)$ in the Standard Model. More...

double	computeBrHtogaga () const
	The Br $(H\to\gamma\gamma)$ in the Standard Model. More...

double	computeBrHtogg () const
	The Br $(H\to gg)$ in the Standard Model. More...

double	computeBrHtomumu () const
	The Br $(H\to \mu\mu)$ in the Standard Model. More...

double	computeBrHtoss () const
	The Br $(H\to ss)$ in the Standard Model. More...

double	computeBrHtotautau () const
	The Br $(H\to \tau\tau)$ in the Standard Model. More...

double	computeBrHtoWW () const
	The Br $(H\to WW)$ in the Standard Model. More...

double	computeBrHtoZga () const
	The Br $(H\to Z\gamma)$ in the Standard Model. More...

double	computeBrHtoZZ () const
	The Br $(H\to ZZ)$ in the Standard Model. More...

double	computeBrHtoZZinv () const
	The Br $(H\to ZZ \to inv)$ in the Standard Model. More...

void	ComputeDeltaR_rem (const double Mw_i, double DeltaR_rem[orders_EW_size]) const
	A method to collect $\Delta r_{\mathrm{rem}}$ computed via subclasses. More...

void	ComputeDeltaRho (const double Mw_i, double DeltaRho[orders_EW_size]) const
	A method to collect $\Delta\rho$ computed via subclasses. More...

double	computeGammaHgaga_tt () const
	The top loop contribution to $H\to\gamma\gamma$ in the Standard Model. More...

double	computeGammaHgaga_tW () const
	The mixed $t-W$ loop contribution to $H\to\gamma\gamma$ in the Standard Model. More...

double	computeGammaHgaga_WW () const
	The $W$ loop contribution to $H\to\gamma\gamma$ in the Standard Model. More...

double	computeGammaHgg_bb () const
	The bottom loop contribution to $H\to gg$ in the Standard Model. More...

double	computeGammaHgg_tb () const
	The top-bottom interference contribution to $H\to gg$ in the Standard Model. More...

double	computeGammaHgg_tt () const
	The top loop contribution to $H\to gg$ in the Standard Model. More...

double	computeGammaHTotal () const
	The Higgs total width in the Standard Model. More...

double	computeGammaHZga_tt () const
	The top loop contribution to $H\to Z\gamma$ in the Standard Model. More...

double	computeGammaHZga_tW () const
	The mixed $t-W$ loop contribution to $H\to Z\gamma$ in the Standard Model. More...

double	computeGammaHZga_WW () const
	The $W$ loop contribution to $H\to Z\gamma$ in the Standard Model. Currently it returns the value of tab 41 in ref. [137]. More...

double	computeSigmaggH (const double sqrt_s) const
	The ggH cross section in the Standard Model. More...

double	computeSigmaggH_bb (const double sqrt_s) const
	The square of the bottom-quark contribution to the ggH cross section in the Standard Model. More...

double	computeSigmaggH_tb (const double sqrt_s) const
	The top-bottom interference contribution to the ggH cross section in the Standard Model. More...

double	computeSigmaggH_tt (const double sqrt_s) const
	The square of the top-quark contribution to the ggH cross section in the Standard Model. More...

double	computeSigmattH (const double sqrt_s) const
	The ttH production cross section in the Standard Model. More...

double	computeSigmaVBF (const double sqrt_s) const
	The VBF cross section in the Standard Model. More...

double	computeSigmaWF (const double sqrt_s) const
	The W fusion contribution $\sigma_{WF}$ to higgs-production cross section in the Standard Model. More...

double	computeSigmaWH (const double sqrt_s) const
	The WH production cross section in the Standard Model. More...

double	computeSigmaZF (const double sqrt_s) const
	The Z fusion contribution $\sigma_{ZF}$ to higgs-production cross section in the Standard Model. More...

double	computeSigmaZH (const double sqrt_s) const
	The ZH production cross section in the Standard Model. More...

double	computeSigmaZWF (const double sqrt_s) const
	The Z W interference fusion contribution $\sigma_{ZWF}$ to higgs-production cross section in the Standard Model. More...

virtual double	cW2 () const

virtual double	cW2 (const double Mw_i) const
	The square of the cosine of the weak mixing angle in the on-shell scheme, denoted as $c_W^2$. More...

double	DeltaAlpha () const
	The total corrections to the electromagnetic coupling $\alpha$ at the $Z$-mass scale, denoted as $\Delta\alpha(M_Z^2)$. More...

double	DeltaAlphaL5q () const
	The sum of the leptonic and the five-flavour hadronic corrections to the electromagnetic coupling $\alpha$ at the $Z$-mass scale, denoted as $\Delta\alpha^{\ell+5q}(M_Z^2)$. More...

double	DeltaAlphaLepton (const double s) const
	Leptonic contribution to the electromagnetic coupling $\alpha$, denoted as $\Delta\alpha_{\mathrm{lept}}(s)$. More...

double	DeltaAlphaTop (const double s) const
	Top-quark contribution to the electromagnetic coupling $\alpha$, denoted as $\Delta\alpha_{\mathrm{top}}(s)$. More...

virtual gslpp::complex	deltaKappaZ_f (const Particle f) const
	Flavour non-universal vertex corrections to $\kappa_Z^l$, denoted by $\Delta\kappa_Z^l$. More...

virtual double	DeltaR () const
	The SM prediction for $\Delta r$ derived from that for the $W$ boson mass. More...

virtual double	DeltaRbar () const
	The SM prediction for $\Delta \overline{r}$ derived from that for the $W$-boson mass. More...

virtual gslpp::complex	deltaRhoZ_f (const Particle f) const
	Flavour non-universal vertex corrections to $\rho_Z^l$, denoted by $\Delta\rho_Z^l$. More...

virtual double	epsilon1 () const
	The SM contribution to the epsilon parameter $\varepsilon_1$. More...

virtual double	epsilon2 () const
	The SM contribution to the epsilon parameter $\varepsilon_2$. More...

virtual double	epsilon3 () const
	The SM contribution to the epsilon parameter $\varepsilon_3$. More...

virtual double	epsilonb () const
	The SM contribution to the epsilon parameter $\varepsilon_b$. More...

virtual gslpp::complex	gA_f (const Particle f) const
	The effective leptonic neutral-current axial-vector coupling $g_A^l$ in the SM. More...

virtual double	Gamma_had () const
	The hadronic decay width of the $Z$ boson, $\Gamma_{h}$. More...

virtual double	Gamma_inv () const
	The invisible partial decay width of the $Z$ boson, $\Gamma_{\mathrm{inv}}$. More...

virtual double	Gamma_Z () const
	The total decay width of the $Z$ boson, $\Gamma_Z$. More...

virtual double	GammaW () const
	The total width of the $W$ boson, $\Gamma_W$. More...

virtual double	GammaW (const Particle fi, const Particle fj) const
	A partial decay width of the $W$ boson decay into a SM fermion pair. More...

virtual double	GammaZ (const Particle f) const
	The $Z\to \ell\bar{\ell}$ partial decay width, $\Gamma_\ell$. More...

double	getAle () const
	A get method to retrieve the fine-structure constant $\alpha$. More...

double	getAlsMz () const
	A get method to access the value of $\alpha_s(M_Z)$. More...

virtual double	getCBd () const
	The ratio of the absolute value of the $B_d$ mixing amplitude over the Standard Model value. More...

virtual double	getCBs () const
	The ratio of the absolute value of the $B_s$ mixing amplitude over the Standard Model value. More...

virtual double	getCCC1 () const
	A virtual implementation for the RealWeakEFTCC class. More...

virtual double	getCCC2 () const
	A virtual implementation for the RealWeakEFTCC class. More...

virtual double	getCCC3 () const
	A virtual implementation for the RealWeakEFTCC class. More...

virtual double	getCCC4 () const
	A virtual implementation for the RealWeakEFTCC class. More...

virtual double	getCCC5 () const
	A virtual implementation for the RealWeakEFTCC class. More...

virtual double	getCDMK () const
	The ratio of the real part of the $K$ mixing amplitude over the Standard Model value. More...

virtual double	getCepsK () const
	The ratio of the imaginary part of the $K$ mixing amplitude over the Standard Model value. More...

CKM	getCKM () const
	A get method to retrieve the member object of type CKM. More...

double	getDAle5Mz () const
	A get method to retrieve the five-flavour hadronic contribution to the electromagnetic coupling, $\Delta\alpha_{\mathrm{had}}^{(5)}(M_Z^2)$. More...

double	getDelGammaZ () const
	A get method to retrieve the theoretical uncertainty in $\Gamma_Z$, denoted as $\delta\,\Gamma_Z$. More...

double	getDelMw () const
	A get method to retrieve the theoretical uncertainty in $M_W$, denoted as $\delta\,M_W$. More...

double	getDelR0b () const
	A get method to retrieve the theoretical uncertainty in $R_b^0$, denoted as $\delta\,R_b^0$. More...

double	getDelR0c () const
	A get method to retrieve the theoretical uncertainty in $R_c^0$, denoted as $\delta\,R_c^0$. More...

double	getDelR0l () const
	A get method to retrieve the theoretical uncertainty in $R_l^0$, denoted as $\delta\,R_l^0$. More...

double	getDelSigma0H () const
	A get method to retrieve the theoretical uncertainty in $\sigma_{Hadron}^0$, denoted as $\delta\,\sigma_{Hadron}^0$. More...

double	getDelSin2th_b () const
	A get method to retrieve the theoretical uncertainty in $\sin^2\theta_{\rm eff}^{b}$, denoted as $\delta\sin^2\theta_{\rm eff}^{b}$. More...

double	getDelSin2th_l () const
	A get method to retrieve the theoretical uncertainty in $\sin^2\theta_{\rm eff}^{\rm lept}$, denoted as $\delta\sin^2\theta_{\rm eff}^{\rm lept}$. More...

double	getDelSin2th_q () const
	A get method to retrieve the theoretical uncertainty in $\sin^2\theta_{\rm eff}^{q\not = b,t}$, denoted as $\delta\sin^2\theta_{\rm eff}^{q\not = b,t}$. More...

std::string	getFlagKappaZ () const
	A method to retrieve the model flag KappaZ. More...

std::string	getFlagMw () const
	A method to retrieve the model flag Mw. More...

std::string	getFlagRhoZ () const
	A method to retrieve the model flag RhoZ. More...

const Flavour &	getFlavour () const

double	getGF () const
	A get method to retrieve the Fermi constant $G_\mu$. More...

int	getIterationNo () const

Particle	getLeptons (const QCD::lepton p) const
	A get method to retrieve the member object of a lepton. More...

virtual double	getMHl () const
	A get method to retrieve the Higgs mass $m_h$. More...

virtual double	getmq (const QCD::quark q, const double mu) const

double	getMuw () const
	A get method to retrieve the matching scale $\mu_W$ around the weak scale. More...

EWSMApproximateFormulae *	getMyApproximateFormulae () const
	A get method to retrieve the member pointer of type EWSMApproximateFormulae. More...

EWSMcache *	getMyEWSMcache () const
	A get method to retrieve the member pointer of type EWSMcache. More...

LeptonFlavour *	getMyLeptonFlavour () const

EWSMOneLoopEW *	getMyOneLoopEW () const
	A get method to retrieve the member pointer of type EWSMOneLoopEW,. More...

EWSMThreeLoopEW *	getMyThreeLoopEW () const

EWSMThreeLoopEW2QCD *	getMyThreeLoopEW2QCD () const

EWSMThreeLoopQCD *	getMyThreeLoopQCD () const

EWSMTwoLoopEW *	getMyTwoLoopEW () const

EWSMTwoLoopQCD *	getMyTwoLoopQCD () const

double	getMz () const
	A get method to access the mass of the $Z$ boson $M_Z$. More...

virtual double	getPhiBd () const
	Half the relative phase of the $B_d$ mixing amplitude w.r.t. the Standard Model one. More...

virtual double	getPhiBs () const
	Half the relative phase of the $B_s$ mixing amplitude w.r.t. the Standard Model one. More...

virtual StandardModel	getTrueSM () const

gslpp::matrix< gslpp::complex >	getUPMNS () const
	A get method to retrieve the object of the PMNS matrix. More...

gslpp::matrix< gslpp::complex >	getVCKM () const
	A get method to retrieve the CKM matrix. More...

gslpp::matrix< gslpp::complex >	getYd () const
	A get method to retrieve the Yukawa matrix of the down-type quarks, $Y_d$. More...

gslpp::matrix< gslpp::complex >	getYe () const
	A get method to retrieve the Yukawa matrix of the charged leptons, $Y_e$. More...

gslpp::matrix< gslpp::complex >	getYn () const
	A get method to retrieve the Yukawa matrix of the neutrinos, $Y_\nu$. More...

gslpp::matrix< gslpp::complex >	getYu () const
	A get method to retrieve the Yukawa matrix of the up-type quarks, $Y_u$. More...

virtual gslpp::complex	gV_f (const Particle f) const
	The effective leptonic neutral-current vector coupling $g_V^l$ in the SM. More...

bool	IsFlagNoApproximateGammaZ () const
	A method to retrieve the model flag NoApproximateGammaZ. More...

bool	IsFlagWithoutNonUniversalVC () const
	A method to retrieve the model flag WithoutNonUniversalVC. More...

virtual gslpp::complex	kappaZ_f (const Particle f) const
	The effective leptonic neutral-current coupling $\kappa_Z^l$ in the SM. More...

virtual double	Mw () const
	The SM prediction for the $W$-boson mass in the on-shell scheme, $M_{W,\mathrm{SM}}$. More...

virtual double	Mw_tree () const
	The tree-level mass of the $W$ boson, $M_W^{\mathrm{tree}}$. More...

double	MwbarFromMw (const double Mw) const
	A method to convert the $W$-boson mass in the experimental/running-width scheme to that in the complex-pole/fixed-width scheme. More...

double	MwFromMwbar (const double Mwbar) const
	A method to convert the $W$-boson mass in the complex-pole/fixed-width scheme to that in the experimental/running-width scheme. More...

double	Mzbar () const
	The $Z$-boson mass $\overline{M}_Z$ in the complex-pole/fixed-width scheme. More...

virtual double	N_nu () const
	The number of neutrinos obtained indirectly from the measurements at the Z pole, $N_{\nu}$. More...

virtual double	R0_f (const Particle f) const
	The ratio $R_\ell^0=\Gamma(Z\to {\rm hadrons})/\Gamma(Z\to \ell^+ \ell^-)$. More...

virtual double	R_inv () const
	The ratio of the invisible and leptonic (electron) decay widths of the $Z$ boson, $R_{inv}$. More...

virtual double	rho_GammaW (const Particle fi, const Particle fj) const
	EW radiative corrections to the width of $W \to f_i \bar{f}_j$, denoted as $\rho^W_{ij}$. More...

virtual gslpp::complex	rhoZ_f (const Particle f) const
	The effective leptonic neutral-current coupling $\rho_Z^l$ in the SM. More...

double	s02 () const
	The square of the sine of the weak mixing angle $s_0^2$ defined without weak radiative corrections. More...

void	setFlagCacheInStandardModel (bool FlagCacheInStandardModel)
	A set method to change the model flag CacheInStandardModel of StandardModel. More...

void	setFlagNoApproximateGammaZ (bool FlagNoApproximateGammaZ)

virtual double	sigma0_had () const
	The hadronic cross section for $e^+e^- \to Z \to \mathrm{hadrons}$ at the $Z$-pole, $\sigma_h^0$. More...

virtual double	sin2thetaEff (const Particle f) const
	The effective weak mixing angle $\sin^2\theta_{\rm eff}^{\,\ell}$ for $Z\ell\bar{\ell}$ at the the $Z$-mass scale. More...

	StandardModel ()
	The default constructor. More...

double	sW2 () const

virtual double	sW2 (const double Mw_i) const
	The square of the sine of the weak mixing angle in the on-shell scheme, denoted as $s_W^2$. More...

virtual double	v () const
	The Higgs vacuum expectation value. More...

virtual	~StandardModel ()
	The default destructor. More...

Public Member Functions inherited from QCD
double	AboveTh (const double mu) const
	The active flavour threshold above the scale $\mu$ as defined in QCD::Thresholds(). More...

void	addParameters (std::vector< std::string > params_i)
	A method to add parameters that are specific to only one set of observables. More...

virtual double	Als (const double mu, const orders order=FULLNLO, bool Nf_thr=true) const

double	Als4 (const double mu) const
	The value of $\alpha_s^{\mathrm{FULLNLO}}$ at any scale $\mu$ with the number of flavours $n_f = 4$. More...

virtual double	AlsByOrder (const double mu, const orders order=FULLNLO, bool Nf_thr=true) const

double	AlsOLD (const double mu, const orders order=FULLNLO) const
	Computes the running strong coupling $\alpha_s(\mu)$ in the $\overline{\mathrm{MS}}$ scheme. In the cases of LO, NLO and FULLNNLO, the coupling is computed with AlsWithInit(). On the other hand, in the cases of NNLO and FULLNNLO, the coupling is computed with AlsWithLambda(). More...

double	AlsWithInit (const double mu, const double alsi, const double mu_i, const orders order) const
	Computes the running strong coupling $\alpha_s(\mu)$ from $\alpha_s(\mu_i)$ in the $\overline{\mathrm{MS}}$ scheme, where it is forbidden to across a flavour threshold in the RG running from $\mu_i$ to $\mu$. More...

double	AlsWithLambda (const double mu, const orders order) const
	Computes the running strong coupling $\alpha_s(\mu)$ in the $\overline{\mathrm{MS}}$ scheme with the use of $\Lambda_{\rm QCD}$. More...

double	BelowTh (const double mu) const
	The active flavour threshold below the scale $\mu$ as defined in QCD::Thresholds(). More...

double	Beta0 (const double nf) const
	The $\beta_0(n_f)$ coefficient for a certain number of flavours $n_f$. More...

double	Beta1 (const double nf) const
	The $\beta_1(n_f)$ coefficient for a certain number of flavours $n_f$. More...

double	Beta2 (const double nf) const
	The $\beta_2(n_f)$ coefficient for a certain number of flavours $n_f$. More...

double	Beta3 (const double nf) const
	The $\beta_3(n_f)$ coefficient for a certain number of flavours $n_f$. More...

void	CacheShift (double cache[][5], int n) const
	A member used to manage the caching for this class. More...

void	CacheShift (int cache[][5], int n) const

orders	FullOrder (orders order) const
	Return the FULLORDER enum corresponding to order. More...

double	Gamma0 (const double nf) const
	The $\gamma_0$ coefficient used to compute the running of a mass. More...

double	Gamma1 (const double nf) const
	The $\gamma_1$ coefficient used to compute the running of a mass. More...

double	Gamma2 (const double nf) const
	The $\gamma_2$ coefficient used to compute the running of a mass. More...

double	getAlsM () const
	A get method to access the value of $\alpha_s(M_{\alpha_s})$. More...

BParameter	getBBd () const
	For getting the bag parameters corresponding to the operator basis $O_1 -O_5$ in $\Delta b = 2$ process in the $B_d$ meson system. More...

BParameter	getBBs () const
	For getting the bag parameters corresponding to the operator basis $O_1 -O_5$ in $\Delta b = 2$ process in the $B_s$ meson system. More...

BParameter	getBD () const
	For getting the bag parameters corresponding to the operator basis $O_1 -O_5$ in $\Delta c = 2$ process in the $D^0$ meson system. More...

BParameter	getBK () const
	For getting the bag parameters corresponding to the operator basis $O_1 -O_5$ in $\Delta s = 2$ process in the $K^0$ meson system. More...

BParameter	getBKd1 () const

BParameter	getBKd3 () const

double	getCF () const
	A get method to access the Casimir factor of QCD. More...

double	getMAls () const
	A get method to access the mass scale $M_{\alpha_s}$ at which the strong coupling constant measurement is provided. More...

Meson	getMesons (const QCD::meson m) const
	A get method to access a meson as an object of the type Meson. More...

double	getMtpole () const
	A get method to access the pole mass of the top quark. More...

double	getMub () const
	A get method to access the threshold between five- and four-flavour theory in GeV. More...

double	getMuc () const
	A get method to access the threshold between four- and three-flavour theory in GeV. More...

double	getMut () const
	A get method to access the threshold between six- and five-flavour theory in GeV. More...

double	getNc () const
	A get method to access the number of colours $N_c$. More...

double	getOptionalParameter (std::string name) const
	A method to get parameters that are specific to only one set of observables. More...

Particle	getQuarks (const QCD::quark q) const
	A get method to access a quark as an object of the type Particle. More...

std::vector< std::string >	getUnknownParameters ()
	A method to get the vector of the parameters that have been specified in the configuration file but not being used. More...

void	initializeBParameter (std::string name_i) const
	A method to initialize B Parameter and the corresponding meson. More...

void	initializeMeson (QCD::meson meson_i) const
	A method to initialize a meson. More...

double	logLambda (const double nf, orders order) const
	Computes $\ln\Lambda_\mathrm{QCD}$ with nf flavours in GeV. More...

double	Mbar2Mp (const double mbar, const orders order=FULLNNLO) const
	Converts the $\overline{\mathrm{MS}}$ mass $m(m)$ to the pole mass. More...

double	Mp2Mbar (const double mp, const orders order=FULLNNLO) const
	Converts a quark pole mass to the corresponding $\overline{\mathrm{MS}}$ mass $m(m)$. More...

double	Mrun (const double mu, const double m, const orders order=FULLNNLO) const
	Computes a running quark mass $m(\mu)$ from $m(m)$. More...

double	Mrun (const double mu_f, const double mu_i, const double m, const orders order=FULLNNLO) const
	Runs a quark mass from $\mu_i$ to $\mu_f$. More...

double	Mrun4 (const double mu_f, const double mu_i, const double m) const
	The running of a mass with the number of flavours $n_f = 4$. More...

double	MS2DRqmass (const double MSbar) const
	Converts a quark mass from the $\overline{\mathrm{MS}}$ scheme to the $\overline{\mathrm{DR}}$ scheme. More...

double	MS2DRqmass (const double MSscale, const double MSbar) const
	Converts a quark mass from the $\overline{\mathrm{MS}}$ scheme to the $\overline{\mathrm{DR}}$ scheme. More...

double	Nf (const double mu) const
	The number of active flavour at scale $\mu$. More...

double	NfThresholdCorrections (double mu, double M, double als, int nf, orders order) const
	Threshold corrections in matching $\alpha_s(n_f+1)$ with $\alpha_s(n_f)$ from eq. (34) of hep-ph/0512060. More...

std::string	orderToString (const orders order) const
	Converts an object of the enum type "orders" to the corresponding string. More...

	QCD ()
	Constructor. More...

void	setNc (double Nc)
	A set method to change the number of colours $N_c$. More...

void	setOptionalParameter (std::string name, double value)
	A method to set the parameter value for the parameters that are specific to only one set of observables. More...

double	Thresholds (const int i) const
	For accessing the active flavour threshold scales. More...

Public Member Functions inherited from Model
void	addMissingModelParameter (const std::string &missingParameterName)

std::vector< std::string >	getmissingModelParameters ()

unsigned int	getMissingModelParametersCount ()

std::string	getModelName () const
	A method to fetch the name of the model. More...

const double &	getModelParam (std::string name) const

bool	isModelGeneralTHDM () const

bool	isModelGeorgiMachacek () const

bool	IsModelInitialized () const
	A method to check if the model is initialized. More...

bool	isModelLinearized () const

bool	isModelParam (std::string name) const

bool	isModelSUSY () const

bool	isModelTHDM () const

bool	isModelTHDMW () const

bool	IsUpdateError () const
	A method to check if there was any error in the model update process. More...

	Model ()
	The default constructor. More...

void	raiseMissingModelParameterCount ()

void	setModelGeneralTHDM ()

void	setModelGeorgiMachacek ()

void	setModelInitialized (bool ModelInitialized)
	A set method to fix the failure or success of the initialization of the model. More...

void	setModelLinearized (bool linearized=true)

void	setModelName (const std::string name)
	A method to set the name of the model. More...

void	setModelSUSY ()

void	setModelTHDM ()

void	setModelTHDMW ()

void	setSliced (bool Sliced)

void	setUpdateError (bool UpdateError)
	A set method to fix the update status as success or failure. More...

virtual	~Model ()
	The default destructor. More...

Static Public Attributes
static const int	NTHDMvars = 16

static std::string	THDMvars [NTHDMvars]

Static Public Attributes inherited from StandardModel
static const double	GeVminus2_to_nb = 389379.338

static const double	Mw_error = 0.00001
	The target accuracy of the iterative calculation of the $W$-boson mass in units of GeV. More...

static const int	NSMvars = 26
	The number of the model parameters in StandardModel. More...

static const int	NumSMParamsForEWPO = 33
	The number of the SM parameters that are relevant to the EW precision observables. More...

static std::string	SMvars [NSMvars]
	A string array containing the labels of the model parameters in StandardModel. More...

Static Public Attributes inherited from QCD
static const int	NQCDvars = 11
	The number of model parameters in QCD. More...

static std::string	QCDvars [NQCDvars]
	An array containing the labels under which all QCD parameters are stored in a vector of ModelParameter via InputParser::ReadParameters(). More...

Protected Member Functions
bool	checkmodelType (const std::string modeltype) const
	A method to check if the model type name in string form is valid. More...

bool	checkRGEorder (const std::string RGEorder) const
	A method to check if the RGE order name in string form is valid. More...

virtual void	setParameter (const std::string name, const double &value)
	A method to set the value of a parameter of THDM. 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	m_q (const QCD::quark q, const double mu, const orders order=FULLNLO) const

double	RAq (const QCD::quark q) const
	The radiator factor associated with the final-state QED and QCD corrections to the the axial-vector-current interactions, $R_A^q(M_Z^2)$. More...

double	resumKappaZ (const double DeltaRho[orders_EW_size], const double deltaKappa_rem[orders_EW_size], const double DeltaRbar_rem, const bool bool_Zbb) const
	A method to compute the real part of the effetvive coupling $\kappa_Z^f$ from $\Delta\rho$, $\delta\rho_{\rm rem}^{f}$ and $\Delta r_{\mathrm{rem}}$. More...

double	resumMw (const double Mw_i, const double DeltaRho[orders_EW_size], const double DeltaR_rem[orders_EW_size]) const
	A method to compute the $W$-boson mass from $\Delta\rho$ and $\Delta r_{\mathrm{rem}}$. More...

double	resumRhoZ (const double DeltaRho[orders_EW_size], const double deltaRho_rem[orders_EW_size], const double DeltaRbar_rem, const bool bool_Zbb) const
	A method to compute the real part of the effective coupling $\rho_Z^f$ from $\Delta\rho$, $\delta\rho_{\rm rem}^{f}$ and $\Delta r_{\mathrm{rem}}$. More...

double	RVh () const
	The singlet vector corrections to the hadronic $Z$-boson width, denoted as $R_V^h$. More...

double	RVq (const QCD::quark q) const
	The radiator factor associated with the final-state QED and QCD corrections to the the vector-current interactions, $R_V^q(M_Z^2)$. More...

double	SchemeToDouble (const std::string scheme) const
	A method to convert a given scheme name in string form into a floating-point number with double precision. More...

double	taub () const
	Top-mass corrections to the $Zb\bar{b}$ vertex, denoted by $\tau_b$. More...

Protected Member Functions inherited from QCD
double	MassOfNf (int nf) const
	The Mbar mass of the heaviest quark in the theory with Nf active flavour. More...

Protected Attributes
Matching< THDMMatching, THDM >	THDMM
	An object of type THDMMatching. More...

Protected Attributes inherited from StandardModel
double	A
	The CKM parameter $A$ in the Wolfenstein parameterization. More...

double	ale
	The fine-structure constant $\alpha$. More...

double	alpha21

double	alpha31

double	AlsMz
	The strong coupling constant at the Z-boson mass, $\alpha_s(M_Z)$. More...

double	dAle5Mz
	The five-flavour hadronic contribution to the electromagnetic coupling, $\Delta\alpha_{\mathrm{had}}^{(5)}(M_Z^2)$. More...

double	delGammaZ
	The theoretical uncertainty in $\Gamma_Z$, denoted as $\delta\,\Gamma_Z$, in GeV. More...

double	delMw
	The theoretical uncertainty in $M_W$, denoted as $\delta\,M_W$, in GeV. More...

double	delR0b
	The theoretical uncertainty in $R_b^0$, denoted as $\delta\,R_b^0$. More...

double	delR0c
	The theoretical uncertainty in $R_c^0$, denoted as $\delta\,R_c^0$. More...

double	delR0l
	The theoretical uncertainty in $R_l^0$, denoted as $\delta\,R_l^0$. More...

double	delsigma0H
	The theoretical uncertainty in $\sigma_{Hadron}^0$, denoted as $\delta\,\sigma_{Hadron}^0$ in nb. More...

double	delSin2th_b
	The theoretical uncertainty in $\sin^2\theta_{\rm eff}^{b}$, denoted as $\delta\sin^2\theta_{\rm eff}^{b}$. More...

double	delSin2th_l
	The theoretical uncertainty in $\sin^2\theta_{\rm eff}^{\rm lept}$, denoted as $\delta\sin^2\theta_{\rm eff}^{\rm lept}$. More...

double	delSin2th_q
	The theoretical uncertainty in $\sin^2\theta_{\rm eff}^{q\not = b,t}$, denoted as $\delta\sin^2\theta_{\rm eff}^{q\not = b,t}$. More...

double	delta

double	etab
	The CKM parameter $\bar{\eta}$ in the Wolfenstein parameterization. More...

bool	flag_order [orders_EW_size]
	An array of internal flags controlling the inclusions of higher-order corrections. More...

double	gamma
	$\gamma $ used as an input for FlagWolfenstein = FALSE More...

double	GF
	The Fermi constant $G_\mu$ in ${\rm GeV}^{-2}$. More...

double	lambda
	The CKM parameter $\lambda$ in the Wolfenstein parameterization. More...

Particle	leptons [6]
	An array of Particle objects for the leptons. More...

double	mHl
	The Higgs mass $m_h$ in GeV. More...

double	muw
	A matching scale $\mu_W$ around the weak scale in GeV. More...

CKM	myCKM
	An object of type CKM. More...

PMNS	myPMNS

double	Mz
	The mass of the $Z$ boson in GeV. More...

bool	requireCKM
	An internal flag to control whether the CKM matrix has to be recomputed. More...

bool	requireYe
	An internal flag to control whether the charged-lepton Yukawa matrix has to be recomputed. More...

bool	requireYn
	An internal flag to control whether the neutrino Yukawa matrix has to be recomputed. More...

double	rhob
	The CKM parameter $\bar{\rho}$ in the Wolfenstein parameterization. More...

double	s12

double	s13

double	s23

Flavour	SMFlavour
	An object of type Flavour. More...

Matching< StandardModelMatching, StandardModel >	SMM
	An object of type Matching. More...

double	Vcb
	$\vert V_{cb} \vert $ used as an input for FlagWolfenstein = FALSE More...

double	Vub
	$\vert V_{ub} \vert $ used as an input for FlagWolfenstein = FALSE More...

double	Vus
	$\vert V_{us} \vert $ used as an input for FlagWolfenstein = FALSE More...

gslpp::matrix< gslpp::complex >	Yd
	The Yukawa matrix of the down-type quarks. More...

gslpp::matrix< gslpp::complex >	Ye
	The Yukawa matrix of the charged leptons. More...

gslpp::matrix< gslpp::complex >	Yn
	The Yukawa matrix of the neutrinos. More...

gslpp::matrix< gslpp::complex >	Yu
	The Yukawa matrix of the up-type quarks. More...

Protected Attributes inherited from QCD
double	AlsM
	The strong coupling constant at the mass scale MAls, $\alpha_s(M_{\alpha_s})$. More...

double	CA

double	CF

bool	computemt
	Switch for computing the $\overline{\mathrm{MS}}$ mass of the top quark. More...

double	dAdA_NA

double	dFdA_NA

double	dFdF_NA

double	MAls
	The mass scale in GeV at which the strong coupling measurement is provided. More...

double	mtpole
	The pole mass of the top quark. More...

double	mub
	The threshold between five- and four-flavour theory in GeV. More...

double	muc
	The threshold between four- and three-flavour theory in GeV. More...

double	mut
	The threshold between six- and five-flavour theory in GeV. More...

double	NA

double	Nc
	The number of colours. More...

Particle	quarks [6]
	The vector of all SM quarks. More...

bool	requireYd
	Switch for generating the Yukawa couplings to the down-type quarks. More...

bool	requireYu
	Switch for generating the Yukawa couplings to the up-type quarks. More...

double	TF

Protected Attributes inherited from Model
bool	isSliced
	A boolean set to true if the current istance is a slice of an extended object. More...

std::map< std::string, std::reference_wrapper< const double > >	ModelParamMap

bool	UpdateError
	A boolean set to false if update is successful. More...

Private Attributes
double	BDstartaunu_A

double	BDstartaunu_B

double	BDstartaunu_SM

double	BDtaunu_A

double	BDtaunu_B

double	BDtaunu_SM

double	bma

double	bsgamma_theoryerror

double	cosb

std::string	flag_model

std::string	flag_RGEorder

bool	flag_use_sq_masses

bool	flag_wfr

double	logtb

double	m12_2

double	mA1

double	mA2

double	mHh1

double	mHh2

double	mHl2

double	mHp1

double	mHp2

THDMcache *	myTHDMcache
	An object of type THDMCache. More...

double	NLOuniscale

double	Q_THDM

double	Rpeps

double	sin_ba

double	sinb

double	tanb

Additional Inherited Members
Public Types inherited from StandardModel
enum	LEP2RCs { Weak = 0, WeakBox, ISR, QEDFSR, QCDFSR, NUMofLEP2RCs }

enum	orders_EW { EW1 = 0, EW1QCD1, EW1QCD2, EW2, EW2QCD1, EW3, orders_EW_size }
	An enumerated type representing perturbative orders of radiative corrections to EW precision observables. More...

Public Types inherited from QCD
enum	lepton { NEUTRINO_1, ELECTRON, NEUTRINO_2, MU, NEUTRINO_3, TAU, NOLEPTON }
	An enum type for leptons. More...

enum	meson { P_0, P_P, K_0, K_P, D_0, D_P, B_D, B_P, B_S, B_C, PHI, K_star, K_star_P, D_star_P, RHO, RHO_P, OMEGA, MESON_END }
	An enum type for mesons. More...

enum	quark { UP, DOWN, CHARM, STRANGE, TOP, BOTTOM }
	An enum type for quarks. More...

Constructor & Destructor Documentation

◆ THDM()

THDM::THDM ( )

THDM constructor.

Definition at line 15 of file THDM.cpp.

            : StandardModel(), THDMM(*this) {
     SMM.setObj((StandardModelMatching&) THDMM.getObj());
     ModelParamMap.insert(std::make_pair("logtb", std::cref(logtb)));
     ModelParamMap.insert(std::make_pair("bma", std::cref(bma)));
     ModelParamMap.insert(std::make_pair("mHh1", std::cref(mHh1)));
     ModelParamMap.insert(std::make_pair("mA1", std::cref(mA1)));
     ModelParamMap.insert(std::make_pair("mHp1", std::cref(mHp1)));
     ModelParamMap.insert(std::make_pair("mHh2", std::cref(mHh2)));
     ModelParamMap.insert(std::make_pair("mA2", std::cref(mA2)));
     ModelParamMap.insert(std::make_pair("mHp2", std::cref(mHp2)));
     ModelParamMap.insert(std::make_pair("m12_2", std::cref(m12_2)));
     ModelParamMap.insert(std::make_pair("BDtaunu_SM", std::cref(BDtaunu_SM)));
     ModelParamMap.insert(std::make_pair("BDtaunu_A", std::cref(BDtaunu_A)));
     ModelParamMap.insert(std::make_pair("BDtaunu_B", std::cref(BDtaunu_B)));
     ModelParamMap.insert(std::make_pair("BDstartaunu_SM", std::cref(BDstartaunu_SM)));
     ModelParamMap.insert(std::make_pair("BDstartaunu_A", std::cref(BDstartaunu_A)));
     ModelParamMap.insert(std::make_pair("BDstartaunu_B", std::cref(BDstartaunu_B)));
     ModelParamMap.insert(std::make_pair("bsgamma_theoryerror", std::cref(bsgamma_theoryerror)));
     ModelParamMap.insert(std::make_pair("Q_THDM", std::cref(Q_THDM)));
     ModelParamMap.insert(std::make_pair("Rpeps", std::cref(Rpeps)));
     ModelParamMap.insert(std::make_pair("NLOuniscale", std::cref(NLOuniscale)));
     flag_use_sq_masses=true;
 }

◆ ~THDM()

THDM::~THDM ( )

THDM destructor.

Definition at line 39 of file THDM.cpp.

            {
     if (IsModelInitialized()) {
             if (myTHDMcache != NULL) delete(myTHDMcache);
         }
 }

Member Function Documentation

◆ checkmodelType()

bool THDM::checkmodelType ( const std::string modeltype ) const

inlineprotected

A method to check if the model type name in string form is valid.

Parameters

[in] modeltype THDM model type name

Returns: a boolean that is true if the model type name is valid

Definition at line 590 of file THDM.h.

     {
         if (modeltype.compare("type1") == 0
                 || modeltype.compare("type2") == 0
                 || modeltype.compare("typeX") == 0
                 || modeltype.compare("typeY") == 0)
             return true;
         else
             return false;
     }

◆ CheckParameters()

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

virtual

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

Parameters

[in] DPars a 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 147 of file THDM.cpp.

                                                                  {
     for (int i = 0; i < NTHDMvars; i++) {
         if (DPars.find(THDMvars[i]) == DPars.end()) {
             std::cout << "ERROR: missing mandatory THDM parameter " << THDMvars[i] << std::endl;
             raiseMissingModelParameterCount();
             addMissingModelParameter(THDMvars[i]);
         }
     }
     return(StandardModel::CheckParameters(DPars));
 }

◆ checkRGEorder()

bool THDM::checkRGEorder ( const std::string RGEorder ) const

inlineprotected

A method to check if the RGE order name in string form is valid.

Parameters

[in] RGEorder THDM RGE order

Returns: a boolean that is true if the RGE order string is valid

Definition at line 606 of file THDM.h.

     {
         if (RGEorder.compare("LO") == 0
                 || RGEorder.compare("approxNLO") == 0
                 || RGEorder.compare("NLO") == 0)
             return true;
         else
             return false;
     }

◆ getBDstartaunu_A()

double THDM::getBDstartaunu_A ( ) const

inline

A method get the $ B\to D^*\tau\nu $ coefficient A.

Returns: $ B\to D^*\tau\nu $ coefficient A

Definition at line 531 of file THDM.h.

                                     {
         return BDstartaunu_A;
     }

◆ getBDstartaunu_B()

double THDM::getBDstartaunu_B ( ) const

inline

A method get the $ B\to D \tau\nu$ coefficient B.

Returns: $ B\to D \tau\nu$ coefficient B

Definition at line 539 of file THDM.h.

                                     {
         return BDstartaunu_B;
     }

◆ getBDstartaunu_SM()

double THDM::getBDstartaunu_SM ( ) const

inline

A method get the $ B\to D^*\tau\nu $ SM expectation.

Returns: $ B\to D^*\tau\nu $ SM expectation

Definition at line 523 of file THDM.h.

                                      {
         return BDstartaunu_SM;
     }

◆ getBDtaunu_A()

double THDM::getBDtaunu_A ( ) const

inline

A method get the $ B\to D \tau\nu$ coefficient A.

Returns: $ B\to D \tau\nu$ coefficient A

Definition at line 507 of file THDM.h.

                                 {
         return BDtaunu_A;
     }

◆ getBDtaunu_B()

double THDM::getBDtaunu_B ( ) const

inline

A method get the $ B\to D \tau\nu$ coefficient B.

Returns: $ B\to D \tau\nu$ coefficient B

Definition at line 515 of file THDM.h.

                                 {
         return BDtaunu_B;
     }

◆ getBDtaunu_SM()

double THDM::getBDtaunu_SM ( ) const

inline

A method get the $ B\to D \tau\nu$ SM expectation.

Returns: $ B\to D \tau\nu$ SM expectation

Definition at line 499 of file THDM.h.

                                  {
         return BDtaunu_SM;
     }

◆ getbma()

double THDM::getbma ( ) const

inline

A method get $\beta-\alpha$.

Returns: $\beta-\alpha$

Definition at line 307 of file THDM.h.

                           {
         return bma;
     }

◆ getbsgamma_theoryerror()

double THDM::getbsgamma_theoryerror ( ) const

inline

A method get the nuisance parameter for the theoretical error on $ b\to s\gamma$.

Returns: nuisance parameter for the theoretical error on $ b\to s\gamma$

Definition at line 547 of file THDM.h.

                                           {
         return bsgamma_theoryerror;
     }

◆ getcosa()

double THDM::getcosa ( ) const

inline

A method get $\cos \alpha$.

Returns: $\cos \alpha$

Definition at line 323 of file THDM.h.

                           {
         return cos(atan(pow(10.,logtb))-bma);
     }

◆ getcosb()

double THDM::getcosb ( ) const

inline

A method get $\cos \beta$.

Returns: $\cos \beta$

Definition at line 299 of file THDM.h.

                            {
         return cosb;
     }

◆ getlogtb()

double THDM::getlogtb ( ) const

inline

A method get $\log_{10}(\tan \beta)$.

Returns: $\log_{10}(\tan \beta)$

Definition at line 275 of file THDM.h.

                             {
         return logtb;
     }

◆ getm12_2()

double THDM::getm12_2 ( ) const

inline

A method get the parameter of the Higgs potential $m_{12}^2$.

Returns: parameter of the Higgs potential $m_{12}^2$

Definition at line 491 of file THDM.h.

                             {
         return m12_2;
     }

◆ getmA()

double THDM::getmA ( ) const

inline

A method get the mass of the pseudoscalar Higgs A.

Returns: mass of the pseudoscalar Higgs A

Definition at line 437 of file THDM.h.

                          {
         if(flag_use_sq_masses) {
             if(mA2 < 0.) {
                 return 0.;
             }
             else
             {
                 return sqrt(mA2);
             }
         }
         else
         {
                 return mA1;
         }
     }

◆ getmA2()

double THDM::getmA2 ( ) const

inline

A method get the squared mass of the pseudoscalar Higgs A.

Returns: squared mass of the pseudoscalar Higgs A

Definition at line 423 of file THDM.h.

                           {
         if(flag_use_sq_masses) {
             return mA2;
         }
         else
         {
             return mA1*mA1;
         }
     }

◆ getMatching()

virtual THDMMatching& THDM::getMatching ( ) const

inlinevirtual

A get method to access the member reference of type THDMMatching.

Returns: a reference to a THDMMatching object

Reimplemented from StandardModel.

Definition at line 185 of file THDM.h.

     {
         return THDMM.getObj();
     }

◆ getmHh()

double THDM::getmHh ( ) const

inline

A method get the mass of the "non-125 GeV" neutral scalar Higgs.

Returns: mass of the "non-125 GeV" neutral scalar Higgs

Definition at line 394 of file THDM.h.

                           {
         if(flag_use_sq_masses) {
             if(mHh2 < 0.) {
                 return 0.;
             }
             else if(mHh2 < mHl2) {
                 return sqrt(mHl2);
             }
             else
             {
                 return sqrt(mHh2);
             }
         }
         else
         {
             if(mHh1*mHh1 < mHl2) {
                 return sqrt(mHl2);
             }
             else
             {
                 return mHh1;
             }
         }
     }

◆ getmHh2()

double THDM::getmHh2 ( ) const

inline

A method get the squared mass of the "non-125 GeV" neutral scalar Higgs.

Returns: squared mass of the "non-125 GeV" neutral scalar Higgs

Definition at line 365 of file THDM.h.

                            {
         if(flag_use_sq_masses) {
             if(mHh2 < 0.) {
                 return 0.;
             }
             else if(mHh2 < mHl2) {
                 return mHl2;
             }
             else
             {
                 return mHh2;
             }
         }
         else
         {
             if(mHh1*mHh1 < mHl2) {
                 return mHl2;
             }
             else
             {
                 return mHh1*mHh1;
             }
         }
     }

◆ getmHl2()

double THDM::getmHl2 ( ) const

inline

A method get the squared mass of the lighter neutral scalar Higgs.

Returns: squared mass of the lighter neutral scalar Higgs

Definition at line 339 of file THDM.h.

                            {
         if(flag_use_sq_masses) {
             if(mHh2 < mHl2) {
                 return mHh2;
             }
             else
             {
                 return mHl2;
             }
         }
         else
         {
             if(mHh1*mHh1 < mHl2) {
                 return mHh1*mHh1;
             }
             else
             {
                 return mHl2;
             }
         }
     }

◆ getmHp()

double THDM::getmHp ( ) const

inline

A method get the charged Higgs mass.

Returns: charged Higgs mass

Definition at line 471 of file THDM.h.

                           {
         if(flag_use_sq_masses) {
             if(mHp2 < 0.) {
                 return 0.;
             }
             else
             {
                 return sqrt(mHp2);
             }
         }
         else
         {
                 return mHp1;
         }
     }

◆ getmHp2()

double THDM::getmHp2 ( ) const

inline

A method get the squared charged Higgs mass.

Returns: squared charged Higgs mass

Definition at line 457 of file THDM.h.

                            {
         if(flag_use_sq_masses) {
             return mHp2;
         }
         else
         {
             return mHp1*mHp1;
         }
     }

◆ getModelTypeflag()

std::string THDM::getModelTypeflag ( ) const

inline

A method get the THDM model type.

Returns: THDM model type

Definition at line 242 of file THDM.h.

                                        {
         return flag_model;
     }

◆ getMyTHDMCache()

THDMcache* THDM::getMyTHDMCache ( ) const

inline

A method get the THDMCache.

Returns: a object of the type THDMCache

Definition at line 213 of file THDM.h.

     {
         return myTHDMcache;
     }

◆ getNLOuniscale()

double THDM::getNLOuniscale ( ) const

inline

A method get the minimal NLO unitarity check scale.

Returns: Minimal NLO unitarity check scale

Definition at line 571 of file THDM.h.

                                   {
         return NLOuniscale;
     }

◆ getQ_THDM()

double THDM::getQ_THDM ( ) const

inline

A method get the THDM scale.

Returns: THDM scale

Definition at line 555 of file THDM.h.

                              {
         return Q_THDM;
     }

◆ getRGEorderflag()

std::string THDM::getRGEorderflag ( ) const

inline

A method get switch for NLO RGE and approximate NLO RGE.

Returns: Switch for NLO RGE and approximate NLO RGE

Definition at line 250 of file THDM.h.

                                       {
         return flag_RGEorder;
     }

◆ getRpeps()

double THDM::getRpeps ( ) const

inline

A method get the minimal R' value.

Returns: Minimal R' value

Definition at line 563 of file THDM.h.

                             {
         return Rpeps;
     }

◆ getsin_ba()

double THDM::getsin_ba ( ) const

inline

A method get $\sin(\beta-\alpha)$.

Returns: $\sin(\beta-\alpha)$

Definition at line 315 of file THDM.h.

                              {
         return sin_ba;
     }

◆ getsina()

double THDM::getsina ( ) const

inline

A method get $\sin \alpha$.

Returns: $\sin \alpha$

Definition at line 331 of file THDM.h.

                           {
         return sin(atan(pow(10.,logtb))-bma);
     }

◆ getsinb()

double THDM::getsinb ( ) const

inline

A method get $\sin \beta$.

Returns: $\sin \beta$

Definition at line 291 of file THDM.h.

                            {
         return sinb;
     }

◆ getsqmassesflag()

bool THDM::getsqmassesflag ( ) const

inline

A method choose if you want to use the THDM masses or rather their squares.

Returns: Choose if you want to use the THDM masses or rather their squares

Definition at line 258 of file THDM.h.

                                  {
         return flag_use_sq_masses;
     }

◆ gettanb()

double THDM::gettanb ( ) const

inline

A method get $\tan \beta$.

Returns: $\tan \beta$

Definition at line 283 of file THDM.h.

                            {
         return tanb;
     }

◆ getv1()

double THDM::getv1 ( ) const

inline

A method get the VEV $v_1$.

Returns: the VEV $v_1$

Definition at line 224 of file THDM.h.

                          {
         return v() * cosb;
     }

◆ getv2()

double THDM::getv2 ( ) const

inline

A method get the VEV $v_2$.

Returns: the VEV $v_2$

Definition at line 232 of file THDM.h.

                          {
         return v() * sinb;
     }

◆ getWFRflag()

bool THDM::getWFRflag ( ) const

inline

A method get the flag to switch on wavefunction renormalization for the NLO unitarity conditions.

Returns: Flag to switch on wavefunction renormalization for the NLO unitarity conditions

Definition at line 267 of file THDM.h.

                             {
         return flag_wfr;
     }

◆ Init()

bool THDM::Init ( const std::map< std::string, double > & DPars )

virtual

Initializes the THDM parameters found in the argument.

Parameters

[in] DPars a map containing the parameters (all as double) to be used in Monte Carlo

Reimplemented from StandardModel.

Definition at line 56 of file THDM.cpp.

                                                       {
     return(StandardModel::Init(DPars));
 }

◆ InitializeModel()

bool THDM::InitializeModel ( )

virtual

A method to initialize the model.

This method, called via InputParser::ReadParameters(), allocates memory to the pointers defined in the current class.

Returns: a boolean that is true if model initialization is successful

Reimplemented from StandardModel.

Definition at line 48 of file THDM.cpp.

 {
     myTHDMcache = new THDMcache(*this);
     setModelInitialized(StandardModel::InitializeModel());
     setModelTHDM();
     return(true);
 }

◆ PostUpdate()

bool THDM::PostUpdate ( )

virtual

The post-update method for THDM.

This method runs all the procedures that are need to be executed after the model is successfully updated.

Returns: a boolean that is true if the execution is successful

Reimplemented from StandardModel.

Definition at line 83 of file THDM.cpp.

 {
     if(!StandardModel::PostUpdate()) return (false);
  
     mHl2=myTHDMcache->updateCache();
  
     return (true);
 }

◆ PreUpdate()

bool THDM::PreUpdate ( )

virtual

The pre-update method for THDM.

Returns: a boolean that is true if the execution is successful

Reimplemented from StandardModel.

Definition at line 60 of file THDM.cpp.

 {
     if(!StandardModel::PreUpdate()) return (false);
  
     return (true);
 }

◆ setFlag()

bool THDM::setFlag	(	const std::string	name,
		const bool	value
	)

virtual

A method to set a flag of THDM.

Parameters

[in]	name	name of a model flag
[in]	value	the boolean to be assigned to the flag specified by name

Returns: a boolean that is true if the execution is successful

Reimplemented from StandardModel.

Definition at line 197 of file THDM.cpp.

 {
     bool res = false;
     if(name.compare("use_sq_masses") == 0) {
         flag_use_sq_masses = value;
         res = true;
         if (!flag_use_sq_masses) {
             THDMvars[std::distance(THDMvars,std::find(THDMvars,THDMvars+NTHDMvars,"mHh2"))] = "mHh1";
             THDMvars[std::distance(THDMvars,std::find(THDMvars,THDMvars+NTHDMvars,"mA2"))] = "mA1";
             THDMvars[std::distance(THDMvars,std::find(THDMvars,THDMvars+NTHDMvars,"mHp2"))] = "mHp1";
             
             ModelParamMap.insert(std::make_pair("mHh1", std::cref(mHh1)));
             ModelParamMap.insert(std::make_pair("mA1", std::cref(mA1)));
             ModelParamMap.insert(std::make_pair("mHp1", std::cref(mHp1)));
         }
     }
     else if(name.compare("wavefunctionrenormalization") == 0) {
         flag_wfr = value;
         res = true;
     }
     else
         res = StandardModel::setFlag(name,value);
  
     return(res);
 }

◆ setFlagStr()

bool THDM::setFlagStr	(	const std::string	name,
		const std::string	value
	)

virtual

A method to set a string flag of THDM.

Parameters

[in]	name	name of a model flag
[in]	value	the string to be assigned to the flag specified by name

Returns: a boolean that is true if the execution is successful

Reimplemented from StandardModel.

Definition at line 161 of file THDM.cpp.

 {
     bool res = false;
     if(name.compare("modelTypeflag") == 0)
     {
         if (checkmodelType(value))
         {
             flag_model = value;
             res = true;
         }
         else
         {
             throw std::runtime_error("THDM::setFlagStr(): Invalid flag "
                 + name + "=" + value);
         }
     }
     else if(name.compare("RGEorder") == 0) {
         if (checkRGEorder(value))
         {
             flag_RGEorder = value;
             res = true;
         }
         else
         {
             throw std::runtime_error("THDM::setFlagStr(): Invalid flag "
                 + name + "=" + value);
         }
     }
     else
     {
         res = StandardModel::setFlagStr(name,value);
     }
  
     return(res);
 }

◆ setParameter()

void THDM::setParameter	(	const std::string	name,
		const double &	value
	)

protectedvirtual

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

Parameters

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

Reimplemented from StandardModel.

Definition at line 92 of file THDM.cpp.

                                                                 {
  
     if(name.compare("logtb") == 0) {
         logtb = value;
         tanb = pow(10.,logtb);
         if(tanb > 0.) {
             sinb = tanb / sqrt(1. + tanb*tanb);
             cosb = 1. / sqrt(1. + tanb*tanb);
         }
         else {
             throw std::runtime_error("error in THDM::SetParameter, tanb < 0!"); 
           }
         }
     else if(name.compare("bma") == 0) {
         bma = value;
         sin_ba = sin(bma);
     }
     else if(name.compare("mHh1") == 0 && !flag_use_sq_masses)
         mHh1 = value;
     else if(name.compare("mA1") == 0 && !flag_use_sq_masses)
         mA1 = value;
     else if(name.compare("mHp1") == 0 && !flag_use_sq_masses)
         mHp1 = value;
     else if(name.compare("mHh2") == 0 && flag_use_sq_masses)
         mHh2 = value;
     else if(name.compare("mA2") == 0 && flag_use_sq_masses)
         mA2 = value;
     else if(name.compare("mHp2") == 0 && flag_use_sq_masses)
         mHp2 = value;
     else if(name.compare("m12_2") == 0)
         m12_2 = value;
     else if(name.compare("BDtaunu_SM") == 0)
         BDtaunu_SM = value;
     else if(name.compare("BDtaunu_A") == 0)
         BDtaunu_A = value;
     else if(name.compare("BDtaunu_B") == 0)
         BDtaunu_B = value;
     else if(name.compare("BDstartaunu_SM") == 0)
         BDstartaunu_SM = value;
     else if(name.compare("BDstartaunu_A") == 0)
         BDstartaunu_A = value;
     else if(name.compare("BDstartaunu_B") == 0)
         BDstartaunu_B = value;
     else if(name.compare("bsgamma_theoryerror") == 0)
         bsgamma_theoryerror = value;
     else if(name.compare("Q_THDM") == 0)
         Q_THDM = value;
     else if(name.compare("Rpeps") == 0)
         Rpeps = value;
     else if(name.compare("NLOuniscale") == 0)
         NLOuniscale = value;
     else
         StandardModel::setParameter(name,value);
 }

◆ Update()

bool THDM::Update ( const std::map< std::string, double > & DPars )

virtual

The update method for THDM.

This method updates all the model parameters with given DPars.

Parameters

[in] DPars a map of the parameters that are being updated in the Monte Carlo run

Returns: a boolean that is true if the execution is successful

Reimplemented from StandardModel.

Definition at line 67 of file THDM.cpp.

                                                         {
  
     if(!PreUpdate()) return (false);
  
     UpdateError = false;
  
     for (std::map<std::string, double>::const_iterator it = DPars.begin(); it != DPars.end(); it++)
         setParameter(it->first, it->second);
  
     if (UpdateError) return (false);
  
     if(!PostUpdate()) return (false);
  
     return (true);
 }

Member Data Documentation

◆ BDstartaunu_A

double THDM::BDstartaunu_A

private

Definition at line 624 of file THDM.h.

◆ BDstartaunu_B

double THDM::BDstartaunu_B

private

Definition at line 624 of file THDM.h.

◆ BDstartaunu_SM

double THDM::BDstartaunu_SM

private

Definition at line 624 of file THDM.h.

◆ BDtaunu_A

double THDM::BDtaunu_A

private

Definition at line 624 of file THDM.h.

◆ BDtaunu_B

double THDM::BDtaunu_B

private

Definition at line 624 of file THDM.h.

◆ BDtaunu_SM

double THDM::BDtaunu_SM

private

Definition at line 624 of file THDM.h.

◆ bma

double THDM::bma

private

Definition at line 622 of file THDM.h.

◆ bsgamma_theoryerror

double THDM::bsgamma_theoryerror

private

Definition at line 622 of file THDM.h.

◆ cosb

double THDM::cosb

private

Definition at line 622 of file THDM.h.

◆ flag_model

std::string THDM::flag_model

private

Definition at line 625 of file THDM.h.

◆ flag_RGEorder

std::string THDM::flag_RGEorder

private

Definition at line 625 of file THDM.h.

◆ flag_use_sq_masses

bool THDM::flag_use_sq_masses

private

Definition at line 626 of file THDM.h.

◆ flag_wfr

bool THDM::flag_wfr

private

Definition at line 626 of file THDM.h.

◆ logtb

double THDM::logtb

private

Definition at line 622 of file THDM.h.

◆ m12_2

double THDM::m12_2

private

Definition at line 622 of file THDM.h.

◆ mA1

double THDM::mA1

private

Definition at line 622 of file THDM.h.

◆ mA2

double THDM::mA2

private

Definition at line 622 of file THDM.h.

◆ mHh1

double THDM::mHh1

private

Definition at line 622 of file THDM.h.

◆ mHh2

double THDM::mHh2

private

Definition at line 622 of file THDM.h.

◆ mHl2

double THDM::mHl2

private

Definition at line 623 of file THDM.h.

◆ mHp1

double THDM::mHp1

private

Definition at line 622 of file THDM.h.

◆ mHp2

double THDM::mHp2

private

Definition at line 622 of file THDM.h.

◆ myTHDMcache

THDMcache* THDM::myTHDMcache

private

An object of type THDMCache.

Definition at line 620 of file THDM.h.

◆ NLOuniscale

double THDM::NLOuniscale

private

Definition at line 622 of file THDM.h.

◆ NTHDMvars

const int THDM::NTHDMvars = 16

static

Definition at line 123 of file THDM.h.

◆ Q_THDM

double THDM::Q_THDM

private

Definition at line 622 of file THDM.h.

◆ Rpeps

double THDM::Rpeps

private

Definition at line 622 of file THDM.h.

◆ sin_ba

double THDM::sin_ba

private

Definition at line 622 of file THDM.h.

◆ sinb

double THDM::sinb

private

Definition at line 622 of file THDM.h.

◆ tanb

double THDM::tanb

private

Definition at line 622 of file THDM.h.

◆ THDMM

Matching<THDMMatching,THDM> THDM::THDMM

mutableprotected

An object of type THDMMatching.

Definition at line 616 of file THDM.h.

◆ THDMvars

std::string THDM::THDMvars

static

Initial value:

= {"logtb","bma","mHh2","mA2","mHp2","m12_2",
                                               "BDtaunu_SM","BDtaunu_A","BDtaunu_B","BDstartaunu_SM","BDstartaunu_A","BDstartaunu_B",
                                               "bsgamma_theoryerror","Q_THDM","Rpeps","NLOuniscale"}

Definition at line 124 of file THDM.h.

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

Label	LaTeX symbol	Description
logtb	\(\log_{10} \tan \beta\)	The decadic logarithm of the tangent of the mixing angle \(\beta \equiv \arctan \frac{v_2}{v_1}\).
bma	\(\beta - \alpha\)	The difference of the mixing angles \(\beta\) and \(\alpha\).
mHh2	\(m_H^2\)	The mass square of the "non-125 GeV" CP-even Higgs state.
mA2	\(m_A^2\)	The squared masses of the CP-odd Higgs \(A\).
mHp2	\(m_{H^+}^2\)	The squared masses of the charged Higgs.
m12_2	\(m_{12}^2\)	The soft \(Z_2\) symmetry breaking parameter.
bsgamma_theoryerror	\(\theta_{b\to s\gamma}^{\text{theo}}\)	A nuisance parameter between -1 and +1 for the theoretical error in the determination of \(BR(B\to X_s \gamma)\).
Q_THDM	\(Q_{\text{THDM}}\)	The THDM scale.
Rpeps	\(\varepsilon_{\text{R'}}\)	Minimal value for which the R' perturbativity criterion should be applied for the unitarity bounds.
NLOuniscale	\(Q_{\text{min}}^{\text{NLOuni}}\)	Minimal scale at which the NLO unitarity conditions are checked.

Label	Value	Description
modelTypeflag	typeI / typeII / typeX / typeY	This flag determines the type of \(Z_2\) symmetry.
RGEorder	LO / approxNLO / NLO	This flag determines the order in perturbation theory of the renormalization group equations.
wavefunctionrenormalization	true / false	Whether to use wavefunction renormalization for NLO unitarity constraints.

Detailed Description

Model parameters

Model flags

Public Member Functions

Static Public Attributes

Protected Member Functions

Protected Attributes

Private Attributes

Additional Inherited Members

Constructor & Destructor Documentation

◆ THDM()

◆ ~THDM()

Member Function Documentation

◆ checkmodelType()

◆ CheckParameters()

◆ checkRGEorder()

◆ getBDstartaunu_A()

◆ getBDstartaunu_B()

◆ getBDstartaunu_SM()

◆ getBDtaunu_A()

◆ getBDtaunu_B()

◆ getBDtaunu_SM()

◆ getbma()

◆ getbsgamma_theoryerror()

◆ getcosa()

◆ getcosb()

◆ getlogtb()

◆ getm12_2()

◆ getmA()

◆ getmA2()

◆ getMatching()

◆ getmHh()

◆ getmHh2()

◆ getmHl2()

◆ getmHp()

◆ getmHp2()

◆ getModelTypeflag()

◆ getMyTHDMCache()

◆ getNLOuniscale()

◆ getQ_THDM()

◆ getRGEorderflag()

◆ getRpeps()

◆ getsin_ba()

◆ getsina()

◆ getsinb()

◆ getsqmassesflag()

◆ gettanb()

◆ getv1()

◆ getv2()

◆ getWFRflag()

◆ Init()

◆ InitializeModel()

◆ PostUpdate()

◆ PreUpdate()

◆ setFlag()

◆ setFlagStr()

◆ setParameter()

◆ Update()

Member Data Documentation

◆ BDstartaunu_A

◆ BDstartaunu_B

◆ BDstartaunu_SM

◆ BDtaunu_A

◆ BDtaunu_B

◆ BDtaunu_SM

◆ bma

◆ bsgamma_theoryerror

◆ cosb

◆ flag_model

◆ flag_RGEorder

◆ flag_use_sq_masses

◆ flag_wfr

◆ logtb

◆ m12_2

◆ mA1

◆ mA2

◆ mHh1

◆ mHh2

◆ mHl2

◆ mHp1