12 = {
"cWB",
"cH",
"cL1L1",
"cL1L2",
"cL1L3",
"cL2L2",
"cL2L3",
"cL3L3",
13 "cHL1p",
"cHL2p",
"cHL3p",
"cHQ1p",
"cHQ2p",
"cHQ3p",
14 "cHL1",
"cHL2",
"cHL3",
"cHQ1",
"cHQ2",
"cHQ3",
"cHE1",
"cHE2",
"cHE3",
15 "cHU1",
"cHU2",
"cHU3",
"cHD1",
"cHD2",
"cHD3",
"Lambda_NP"};
18 = {
"cWB_NP",
"cH_NP",
"cLL_NP",
"cHLp_NP",
"cHL_NP",
19 "cHQ1p_NP",
"cHQ2p_NP",
"cHQ3p_NP",
"cHQ1_NP",
"cHQ2_NP",
"cHQ3_NP",
20 "cHU1_NP",
"cHU2_NP",
"cHU3_NP",
"cHD1_NP",
"cHD2_NP",
"cHD3_NP",
"cHE_NP",
"Lambda_NP"};
23 = {
"cWB_NP",
"cH_NP",
"cL1L1",
"cL1L2",
"cL1L3",
"cL2L2",
"cL2L3",
"cL3L3",
24 "cHL1p",
"cHL2p",
"cHL3p",
"cHQp_NP",
"cHL1",
"cHL2",
"cHL3",
25 "cHQ_NP",
"cHE1",
"cHE2",
"cHE3",
"cHU_NP",
"cHD_NP",
"Lambda_NP"};
28 = {
"cWB_NP",
"cH_NP",
"cLL_NP",
"cHLp_NP",
"cHQp_NP",
29 "cHL_NP",
"cHQ_NP",
"cHE_NP",
"cHU_NP",
"cHD_NP",
"Lambda_NP"};
32 :
NPbase(), FlagLeptonUniversal(FlagLeptonUniversal_in), FlagQuarkUniversal(FlagQuarkUniversal_in)
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)));
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)));
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)));
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)));
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)));
77 ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >(
"Lambda_NP", boost::cref(
Lambda_NP)));
85 if (name.compare(
"cWB_NP") == 0)
87 else if (name.compare(
"cH_NP") == 0)
89 else if (name.compare(
"cL1L1_NP") == 0)
91 else if (name.compare(
"cL1L2_NP") == 0)
93 else if (name.compare(
"cL1L3_NP") == 0)
95 else if (name.compare(
"cL2L2_NP") == 0)
97 else if (name.compare(
"cL2L3_NP") == 0)
99 else if (name.compare(
"cL3L3_NP") == 0)
101 else if (name.compare(
"cLL_NP") == 0) {
108 }
else if (name.compare(
"cHL1p_NP") == 0)
110 else if (name.compare(
"cHL2p_NP") == 0)
112 else if (name.compare(
"cHL3p_NP") == 0)
114 else if (name.compare(
"cHLp_NP") == 0) {
118 }
else if (name.compare(
"cHQ1p_NP") == 0)
120 else if (name.compare(
"cHQ2p_NP") == 0)
122 else if (name.compare(
"cHQ3p_NP") == 0)
124 else if (name.compare(
"cHQp_NP") == 0) {
128 }
else if (name.compare(
"cHL1_NP") == 0)
130 else if (name.compare(
"cHL2_NP") == 0)
132 else if (name.compare(
"cHL3_NP") == 0)
134 else if (name.compare(
"cHL_NP") == 0) {
138 }
else if (name.compare(
"cHQ1_NP") == 0)
140 else if (name.compare(
"cHQ2_NP") == 0)
142 else if (name.compare(
"cHQ3_NP") == 0)
144 else if (name.compare(
"cHQ_NP") == 0) {
148 }
else if (name.compare(
"cHE1_NP") == 0)
150 else if (name.compare(
"cHE2_NP") == 0)
152 else if (name.compare(
"cHE3_NP") == 0)
154 else if (name.compare(
"cHE_NP") == 0) {
158 }
else if (name.compare(
"cHU1_NP") == 0)
160 else if (name.compare(
"cHU2_NP") == 0)
162 else if (name.compare(
"cHU3_NP") == 0)
164 else if (name.compare(
"cHU_NP") == 0) {
168 }
else if (name.compare(
"cHD1_NP") == 0)
170 else if (name.compare(
"cHD2_NP") == 0)
172 else if (name.compare(
"cHD3_NP") == 0)
174 else if (name.compare(
"cHD_NP") == 0) {
178 }
else if (name.compare(
"Lambda_NP") == 0)
189 std::cout <<
"ERROR: Missing mandatory NPEffectiveBS_LFU_QFU parameter "
197 std::cout <<
"ERROR: Missing mandatory NPEffectiveBS_LFU parameter "
205 std::cout <<
"ERROR: Missing mandatory NPEffectiveBS_QFU parameter "
213 std::cout <<
"ERROR: Missing mandatory NPEffectiveBS parameter "
219 throw std::runtime_error(
"Error in NPEffectiveBS::CheckParameters()");
240 double tmp = 4.0 * M_PI *
ale /
sqrt(2.0) / GF0 /
Mz /
Mz;
263 Gamma_W *= 1.0 - 3.0 * alpha / 4.0 / (c2 - s2)
265 - (1.0 + c2) / 2.0 / (c2 - s2) *
DeltaGF()
281 return ( 4.0 * sW_SM * cW_SM *
cWB /
alphaMz() * ratio);
315 if (f.
is(
"NEUTRINO_1"))
317 if (f.
is(
"NEUTRINO_2"))
319 if (f.
is(
"NEUTRINO_3"))
321 if (f.
is(
"ELECTRON"))
339 throw std::runtime_error(
"Error in NPEffectiveBS::deltaGL_f()");
345 if (f.
is(
"NEUTRINO_1"))
347 if (f.
is(
"NEUTRINO_2"))
349 if (f.
is(
"NEUTRINO_3"))
351 if (f.
is(
"ELECTRON"))
352 return ( -
cHE1 / 2.0 * ratio);
354 return ( -
cHE2 / 2.0 * ratio);
356 return ( -
cHE3 / 2.0 * ratio);
358 return ( -
cHU1 / 2.0 * ratio);
360 return ( -
cHU2 / 2.0 * ratio);
364 return ( -
cHD1 / 2.0 * ratio);
366 return ( -
cHD2 / 2.0 * ratio);
368 return ( -
cHD3 / 2.0 * ratio);
369 throw std::runtime_error(
"Error in NPEffectiveBS::deltaGL_f()");
double cHL3
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 the model.
virtual double deltaGA_f(const Particle f) const
New physics contribution to the neutral-current axial-vector coupling .
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 ...
The auxiliary base model class for other model classes.
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.
virtual double deltaGV_f(const Particle f) const
New physics contribution to the neutral-current vector coupling .
double cHQ3p
The dimension-6 operator coefficient .
virtual double DeltaGF() const
New physics contribution to the Fermi constant.
double cHL3p
The dimension-6 operator coefficient .
virtual double obliqueU() const
The oblique parameter .
double cHQ2p
The dimension-6 operator coefficient .
double cL2L2
The dimension-6 operator coefficient .
double deltaGR_f_tmp(const Particle f) const
New physics contribution to .
virtual double obliqueT() const
The oblique parameter .
double ale
The fine-structure constant .
static const int NNPEffectiveBSVars_QFU
The number of the model parameters in NPEffectiveBS with quark flavour universality.
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 initializ...
double GF
The Fermi constant in .
static const std::string NPEffectiveBSVars[NNPEffectiveBSVars]
A string array containing the labels of the model parameters in NPEffectiveBS.
double cHQ1
The dimension-6 operator coefficient .
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.
bool is(std::string name_i) const
double Mz
The mass of the boson in GeV.
virtual double v() const
The Higgs vacuum expectation value. where is the Fermi constant, measured through muon decays...
double cWB
The dimension-6 operator coefficient .
double deltaGL_f_tmp(const Particle f) const
New physics contribution to .
static const std::string NPEffectiveBSVars_LFU[NNPEffectiveBSVars_LFU]
A string array containing the labels of the model parameters in NPEffectiveBS with lepton flavour uni...
NPEffectiveBS(const bool FlagLeptonUniversal_in=false, const bool FlagQuarkUniversal_in=false)
Constructor.
virtual double GammaW(const Particle fi, const Particle fj) const
A partial decay width of the boson decay into a SM fermion pair.
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 .
virtual double v() const
The SM Higgs vacuum expectation value .
double cL1L2
The dimension-6 operator coefficient .
double cHQ2
The dimension-6 operator coefficient .
double cL3L3
The dimension-6 operator coefficient .
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...
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 .
virtual double GammaW() const
The total width of the boson, .
double cHE2
The dimension-6 operator coefficient .
virtual double obliqueS() const
The oblique parameter .
double cHD3
The dimension-6 operator coefficient .
virtual double Mw_tree() const
The tree-level mass of the boson, .
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 cHL1
The dimension-6 operator coefficient .
double cHD2
The dimension-6 operator coefficient .
double alphaMz() const
The electromagnetic coupling at the -mass scale, .
double cHL1p
The dimension-6 operator coefficient .
double cHU3
The dimension-6 operator coefficient .
static const int NNPEffectiveBSVars_LFU
The number of the model parameters in NPEffectiveBS with lepton flavour universality.
std::map< std::string, boost::reference_wrapper< const double > > ModelParamMap
double cHL2p
The dimension-6 operator coefficient .
virtual double deltaGA_f(const Particle f) const
New physics contribution to the neutral-current axial-vector coupling .
double cHL2
The dimension-6 operator coefficient .
virtual double deltaGV_f(const Particle f) const
New physics contribution to the neutral-current vector coupling .
complex sqrt(const complex &z)
double cH
The dimension-6 operator coefficient .