a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
NPEffectiveBS.cpp
Go to the documentation of this file.
1 /*
2  * Copyright (C) 2013 HEPfit Collaboration
3  *
4  *
5  * For the licensing terms see doc/COPYING.
6  */
7 
8 #include "NPEffectiveBS.h"
9 #include <stdexcept>
10 
11 const std::string NPEffectiveBS::NPEffectiveBSVars[NNPEffectiveBSVars]
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"};
16 
17 const std::string NPEffectiveBS::NPEffectiveBSVars_LFU[NNPEffectiveBSVars_LFU]
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"};
21 
22 const std::string NPEffectiveBS::NPEffectiveBSVars_QFU[NNPEffectiveBSVars_QFU]
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"};
26 
27 const std::string NPEffectiveBS::NPEffectiveBSVars_LFU_QFU[NNPEffectiveBSVars_LFU_QFU]
28  = {"cWB_NP", "cH_NP", "cLL_NP", "cHLp_NP", "cHQp_NP",
29  "cHL_NP", "cHQ_NP", "cHE_NP", "cHU_NP", "cHD_NP", "Lambda_NP"};
30 
31 NPEffectiveBS::NPEffectiveBS(const bool FlagLeptonUniversal_in, const bool FlagQuarkUniversal_in)
32 : NPbase(), FlagLeptonUniversal(FlagLeptonUniversal_in), FlagQuarkUniversal(FlagQuarkUniversal_in)
33 {
34  ModelParamMap.insert(std::make_pair("cWB_NP", std::cref(cWB)));
35  ModelParamMap.insert(std::make_pair("cH_NP", std::cref(cH)));
36  if (FlagLeptonUniversal) {
37  ModelParamMap.insert(std::make_pair("cLL_NP", std::cref(cL1L1)));
38  ModelParamMap.insert(std::make_pair("cHLp_NP", std::cref(cHL1p)));
39  ModelParamMap.insert(std::make_pair("cHL_NP", std::cref(cHL1)));
40  ModelParamMap.insert(std::make_pair("cHE_NP", std::cref(cHE1)));
41  } else {
42  ModelParamMap.insert(std::make_pair("cL1L1_NP", std::cref(cL1L1)));
43  ModelParamMap.insert(std::make_pair("cL1L2_NP", std::cref(cL1L2)));
44  ModelParamMap.insert(std::make_pair("cL1L3_NP", std::cref(cL1L3)));
45  ModelParamMap.insert(std::make_pair("cL2L2_NP", std::cref(cL2L2)));
46  ModelParamMap.insert(std::make_pair("cL2L3_NP", std::cref(cL2L3)));
47  ModelParamMap.insert(std::make_pair("cL3L3_NP", std::cref(cL3L3)));
48  ModelParamMap.insert(std::make_pair("cHL1p_NP", std::cref(cHL1p)));
49  ModelParamMap.insert(std::make_pair("cHL2p_NP", std::cref(cHL2p)));
50  ModelParamMap.insert(std::make_pair("cHL3p_NP", std::cref(cHL3p)));
51  ModelParamMap.insert(std::make_pair("cHL1_NP", std::cref(cHL1)));
52  ModelParamMap.insert(std::make_pair("cHL2_NP", std::cref(cHL2)));
53  ModelParamMap.insert(std::make_pair("cHL3_NP", std::cref(cHL3)));
54  ModelParamMap.insert(std::make_pair("cHE1_NP", std::cref(cHE1)));
55  ModelParamMap.insert(std::make_pair("cHE2_NP", std::cref(cHE2)));
56  ModelParamMap.insert(std::make_pair("cHE3_NP", std::cref(cHE3)));
57  }
58  if (FlagQuarkUniversal) {
59  ModelParamMap.insert(std::make_pair("cHQp_NP", std::cref(cHQ1p)));
60  ModelParamMap.insert(std::make_pair("cHQ_NP", std::cref(cHQ1)));
61  ModelParamMap.insert(std::make_pair("cHU_NP", std::cref(cHU1)));
62  ModelParamMap.insert(std::make_pair("cHD_NP", std::cref(cHD1)));
63  } else {
64  ModelParamMap.insert(std::make_pair("cHQ1p_NP", std::cref(cHQ1p)));
65  ModelParamMap.insert(std::make_pair("cHQ2p_NP", std::cref(cHQ2p)));
66  ModelParamMap.insert(std::make_pair("cHQ3p_NP", std::cref(cHQ3p)));
67  ModelParamMap.insert(std::make_pair("cHQ1_NP", std::cref(cHQ1)));
68  ModelParamMap.insert(std::make_pair("cHQ2_NP", std::cref(cHQ2)));
69  ModelParamMap.insert(std::make_pair("cHQ3_NP", std::cref(cHQ3)));
70  ModelParamMap.insert(std::make_pair("cHU1_NP", std::cref(cHU1)));
71  ModelParamMap.insert(std::make_pair("cHU2_NP", std::cref(cHU2)));
72  ModelParamMap.insert(std::make_pair("cHU3_NP", std::cref(cHU3)));
73  ModelParamMap.insert(std::make_pair("cHD1_NP", std::cref(cHD1)));
74  ModelParamMap.insert(std::make_pair("cHD2_NP", std::cref(cHD2)));
75  ModelParamMap.insert(std::make_pair("cHD3_NP", std::cref(cHD3)));
76  }
77  ModelParamMap.insert(std::make_pair("Lambda_NP", std::cref(Lambda_NP)));
78 }
79 
80 
82 
83 void NPEffectiveBS::setParameter(const std::string name, const double& value)
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 }
183 
184 bool NPEffectiveBS::CheckParameters(const std::map<std::string, double>& DPars)
185 {
186  if (FlagLeptonUniversal && FlagQuarkUniversal) {
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  raiseMissingModelParameterCount();
192  addMissingModelParameter(NPEffectiveBSVars_LFU_QFU[i]);
193  }
194  }
195  } else if (FlagLeptonUniversal && !FlagQuarkUniversal) {
196  for (int i = 0; i < NNPEffectiveBSVars_LFU; i++) {
197  if (DPars.find(NPEffectiveBSVars_LFU[i]) == DPars.end()) {
198  std::cout << "ERROR: Missing mandatory NPEffectiveBS_LFU parameter "
199  << NPEffectiveBSVars_LFU[i] << std::endl;
200  raiseMissingModelParameterCount();
201  addMissingModelParameter(NPEffectiveBSVars_LFU[i]);
202  }
203  }
204  } else if (!FlagLeptonUniversal && FlagQuarkUniversal) {
205  for (int i = 0; i < NNPEffectiveBSVars_QFU; i++) {
206  if (DPars.find(NPEffectiveBSVars_QFU[i]) == DPars.end()) {
207  std::cout << "ERROR: Missing mandatory NPEffectiveBS_QFU parameter "
208  << NPEffectiveBSVars_QFU[i] << std::endl;
209  raiseMissingModelParameterCount();
210  addMissingModelParameter(NPEffectiveBSVars_QFU[i]);
211  }
212  }
213  } else if (!FlagLeptonUniversal && !FlagQuarkUniversal) {
214  for (int i = 0; i < NNPEffectiveBSVars; i++) {
215  if (DPars.find(NPEffectiveBSVars[i]) == DPars.end()) {
216  std::cout << "ERROR: Missing mandatory NPEffectiveBS parameter "
217  << NPEffectiveBSVars[i] << std::endl;
218  raiseMissingModelParameterCount();
219  addMissingModelParameter(NPEffectiveBSVars[i]);
220  }
221  }
222  } else
223  throw std::runtime_error("Error in NPEffectiveBS::CheckParameters()");
224 
225  return (NPbase::CheckParameters(DPars));
226 }
227 
228 
229 
231 
232 double NPEffectiveBS::v() const
233 {
234  //return ( sqrt( (1.0 - (cL1L2 - cHL1p - cHL2p)/sqrt(2.0)/GF/LambdaNP/LambdaNP)
235  // /sqrt(2.0)/GF ) );
236 
237  /* use the tree-level relation */
238  return trueSM.v();
239 }
240 
241 double NPEffectiveBS::Mw_tree() const
242 {
243  double GF0 = GF * (1.0 - DeltaGF());
244  double tmp = 4.0 * M_PI * ale / sqrt(2.0) / GF0 / Mz / Mz;
245  return ( Mz / sqrt(2.0) * sqrt(1.0 + sqrt(1.0 - tmp)));
246 }
247 
248 double NPEffectiveBS::DeltaGF() const
249 {
250  double ratio = v() * v() / Lambda_NP / Lambda_NP;
251 
252  return ( -(cL1L2 - cHL1p - cHL2p) * ratio);
253 }
254 
255 
257 
258 double NPEffectiveBS::GammaW() const
259 {
260  double Gamma_W = trueSM.GammaW();
261 
262  double alpha = alphaMz();
263  double c2 = trueSM.cW2();
264  double s2 = trueSM.sW2();
265  double ratio = v() * v() / Lambda_NP / Lambda_NP;
266 
267  Gamma_W *= 1.0 - 3.0 * alpha / 4.0 / (c2 - s2)
268  *(obliqueS() - 2.0 * c2 * obliqueT() - (c2 - s2) * obliqueU() / 2.0 / s2)
269  - (1.0 + c2) / 2.0 / (c2 - s2) * DeltaGF()
270  //+ (cHL1p + cHL2p + cHL3p + cHQ1p + cHQ2p) * ratio; // incorrect
271  + 2.0 / 9.0 * (cHL1p + cHL2p + cHL3p + 3.0 * cHQ1p + 3.0 * cHQ2p) * ratio;
272 
273  return Gamma_W;
274 }
275 
276 
278 
279 double NPEffectiveBS::obliqueS() const
280 {
281  double ratio = v() * v() / Lambda_NP / Lambda_NP;
282  double sW_SM = sqrt(trueSM.sW2()); /* This has to be the SM value. */
283  double cW_SM = sqrt(trueSM.cW2()); /* This has to be the SM value. */
284 
285  return ( 4.0 * sW_SM * cW_SM * cWB / alphaMz() * ratio);
286 }
287 
288 double NPEffectiveBS::obliqueT() const
289 {
290  double ratio = v() * v() / Lambda_NP / Lambda_NP;
291 
292  return ( -cH / 2.0 / alphaMz() * ratio);
293 }
294 
295 double NPEffectiveBS::obliqueU() const
296 {
297  return 0.0;
298 }
299 
300 
302 
303 double NPEffectiveBS::deltaGV_f(const Particle f) const
304 {
305  return ( deltaGL_f_tmp(f) + deltaGR_f_tmp(f) + NPbase::deltaGV_f(f));
306 }
307 
308 double NPEffectiveBS::deltaGA_f(const Particle f) const
309 {
310  return ( deltaGL_f_tmp(f) - deltaGR_f_tmp(f) + NPbase::deltaGA_f(f));
311 }
312 
313 
315 
316 double NPEffectiveBS::deltaGL_f_tmp(const Particle f) const
317 {
318  double ratio = v() * v() / Lambda_NP / Lambda_NP;
319  if (f.is("NEUTRINO_1"))
320  return ( (cHL1p - cHL1) / 2.0 * ratio);
321  if (f.is("NEUTRINO_2"))
322  return ( (cHL2p - cHL2) / 2.0 * ratio);
323  if (f.is("NEUTRINO_3"))
324  return ( (cHL3p - cHL3) / 2.0 * ratio);
325  if (f.is("ELECTRON"))
326  return ( -(cHL1p + cHL1) / 2.0 * ratio);
327  if (f.is("MU"))
328  return ( -(cHL2p + cHL2) / 2.0 * ratio);
329  if (f.is("TAU"))
330  return ( -(cHL3p + cHL3) / 2.0 * ratio);
331  if (f.is("UP"))
332  return ( (cHQ1p - cHQ1) / 2.0 * ratio);
333  if (f.is("CHARM"))
334  return ( (cHQ2p - cHQ2) / 2.0 * ratio);
335  if (f.is("TOP"))
336  return 0.0;
337  if (f.is("DOWN"))
338  return ( -(cHQ1p + cHQ1) / 2.0 * ratio);
339  if (f.is("STRANGE"))
340  return ( -(cHQ2p + cHQ2) / 2.0 * ratio);
341  if (f.is("BOTTOM"))
342  return ( -(cHQ3p + cHQ3) / 2.0 * ratio);
343  throw std::runtime_error("Error in NPEffectiveBS::deltaGL_f()");
344 }
345 
346 double NPEffectiveBS::deltaGR_f_tmp(const Particle f) const
347 {
348  double ratio = v() * v() / Lambda_NP / Lambda_NP;
349  if (f.is("NEUTRINO_1"))
350  return ( 0.);
351  if (f.is("NEUTRINO_2"))
352  return ( 0.);
353  if (f.is("NEUTRINO_3"))
354  return ( 0.);
355  if (f.is("ELECTRON"))
356  return ( -cHE1 / 2.0 * ratio);
357  if (f.is("MU"))
358  return ( -cHE2 / 2.0 * ratio);
359  if (f.is("TAU"))
360  return ( -cHE3 / 2.0 * ratio);
361  if (f.is("UP"))
362  return ( -cHU1 / 2.0 * ratio);
363  if (f.is("CHARM"))
364  return ( -cHU2 / 2.0 * ratio);
365  if (f.is("TOP"))
366  return 0.0;
367  if (f.is("DOWN"))
368  return ( -cHD1 / 2.0 * ratio);
369  if (f.is("STRANGE"))
370  return ( -cHD2 / 2.0 * ratio);
371  if (f.is("BOTTOM"))
372  return ( -cHD3 / 2.0 * ratio);
373  throw std::runtime_error("Error in NPEffectiveBS::deltaGL_f()");
374 }
375 
376 
NPEffectiveBS::NPEffectiveBSVars_LFU_QFU
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...
Definition: NPEffectiveBS.h:309
NPEffectiveBS::cHL1
double cHL1
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:480
StandardModel::cW2
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 .
Definition: StandardModel.cpp:1015
StandardModel::setParameter
virtual void setParameter(const std::string name, const double &value)
A method to set the value of a parameter of StandardModel.
Definition: StandardModel.cpp:257
StandardModel::v
virtual double v() const
The Higgs vacuum expectation value.
Definition: StandardModel.cpp:943
NPEffectiveBS::cL1L1
double cL1L1
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:468
NPEffectiveBS::cL3L3
double cL3L3
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:473
NPEffectiveBS::cHQ3
double cHQ3
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:485
NPEffectiveBS::NPEffectiveBSVars_QFU
static const std::string NPEffectiveBSVars_QFU[NNPEffectiveBSVars_QFU]
A string array containing the labels of the model parameters in NPEffectiveBS with quark flavour univ...
Definition: NPEffectiveBS.h:303
NPEffectiveBS::cHD1
double cHD1
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:492
Particle::is
bool is(std::string name_i) const
Definition: Particle.cpp:23
Particle
A class for particles.
Definition: Particle.h:26
NPEffectiveBS::v
virtual double v() const
The SM Higgs vacuum expectation value .
Definition: NPEffectiveBS.cpp:232
NPEffectiveBS::cHQ1p
double cHQ1p
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:477
NPEffectiveBS::obliqueT
virtual double obliqueT() const
The oblique parameter .
Definition: NPEffectiveBS.cpp:288
NPEffectiveBS::FlagQuarkUniversal
const bool FlagQuarkUniversal
An internal boolean flag that is true if assuming quark flavour universality.
Definition: NPEffectiveBS.h:510
NPEffectiveBS::cHE2
double cHE2
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:487
StandardModel::GF
double GF
The Fermi constant in .
Definition: StandardModel.h:2555
Model::addMissingModelParameter
void addMissingModelParameter(const std::string &missingParameterName)
Definition: Model.h:240
StandardModel::CheckParameters
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...
Definition: StandardModel.cpp:339
StandardModel::alphaMz
double alphaMz() const
The electromagnetic coupling at the -mass scale, .
Definition: StandardModel.cpp:893
NPEffectiveBS::deltaGA_f
virtual double deltaGA_f(const Particle f) const
New physics contribution to the neutral-current axial-vector coupling .
Definition: NPEffectiveBS.cpp:308
NPbase::deltaGA_f
virtual double deltaGA_f(const Particle f) const
New physics contribution to the neutral-current axial-vector coupling .
Definition: NPbase.cpp:107
NPEffectiveBS::NPEffectiveBS
NPEffectiveBS(const bool FlagLeptonUniversal_in=false, const bool FlagQuarkUniversal_in=false)
Constructor.
Definition: NPEffectiveBS.cpp:31
NPEffectiveBS::NNPEffectiveBSVars_QFU
static const int NNPEffectiveBSVars_QFU
The number of the model parameters in NPEffectiveBS with quark flavour universality.
Definition: NPEffectiveBS.h:279
NPEffectiveBS::cHD3
double cHD3
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:494
NPEffectiveBS::cH
double cH
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:467
NPEffectiveBS::cHL1p
double cHL1p
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:474
StandardModel::sW2
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 .
Definition: StandardModel.cpp:1026
StandardModel::ale
double ale
The fine-structure constant .
Definition: StandardModel.h:2556
Model::ModelParamMap
std::map< std::string, std::reference_wrapper< const double > > ModelParamMap
Definition: Model.h:270
NPEffectiveBS::cHL2
double cHL2
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:481
NPEffectiveBS::cL2L2
double cL2L2
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:471
NPEffectiveBS::cHL3
double cHL3
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:482
NPEffectiveBS::cHQ2
double cHQ2
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:484
NPbase::trueSM
StandardModel trueSM
Definition: NPbase.h:2787
NPEffectiveBS.h
NPEffectiveBS::obliqueU
virtual double obliqueU() const
The oblique parameter .
Definition: NPEffectiveBS.cpp:295
NPEffectiveBS::deltaGR_f_tmp
double deltaGR_f_tmp(const Particle f) const
New physics contribution to .
Definition: NPEffectiveBS.cpp:346
NPbase
The auxiliary base model class for other model classes.
Definition: NPbase.h:66
NPEffectiveBS::cHU1
double cHU1
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:489
Model::raiseMissingModelParameterCount
void raiseMissingModelParameterCount()
Definition: Model.h:250
gslpp::sqrt
complex sqrt(const complex &z)
Definition: gslpp_complex.cpp:385
NPEffectiveBS::cHQ2p
double cHQ2p
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:478
NPEffectiveBS::cHL2p
double cHL2p
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:475
NPEffectiveBS::cL1L3
double cL1L3
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:470
NPEffectiveBS::Mw_tree
virtual double Mw_tree() const
The tree-level mass of the boson, .
Definition: NPEffectiveBS.cpp:241
NPEffectiveBS::cHL3p
double cHL3p
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:476
StandardModel::GammaW
virtual double GammaW(const Particle fi, const Particle fj) const
A partial decay width of the boson decay into a SM fermion pair.
Definition: StandardModel.cpp:1166
NPbase::deltaGV_f
virtual double deltaGV_f(const Particle f) const
New physics contribution to the neutral-current vector coupling .
Definition: NPbase.cpp:85
NPEffectiveBS::FlagLeptonUniversal
const bool FlagLeptonUniversal
An internal boolean flag that is true if assuming lepton flavour universality.
Definition: NPEffectiveBS.h:504
NPEffectiveBS::NPEffectiveBSVars
static const std::string NPEffectiveBSVars[NNPEffectiveBSVars]
A string array containing the labels of the model parameters in NPEffectiveBS.
Definition: NPEffectiveBS.h:291
NPEffectiveBS::obliqueS
virtual double obliqueS() const
The oblique parameter .
Definition: NPEffectiveBS.cpp:279
NPEffectiveBS::cWB
double cWB
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:466
NPEffectiveBS::Lambda_NP
double Lambda_NP
The new physics scale .
Definition: NPEffectiveBS.h:495
NPEffectiveBS::NNPEffectiveBSVars_LFU
static const int NNPEffectiveBSVars_LFU
The number of the model parameters in NPEffectiveBS with lepton flavour universality.
Definition: NPEffectiveBS.h:273
NPEffectiveBS::cHD2
double cHD2
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:493
NPEffectiveBS::GammaW
virtual double GammaW() const
A partial decay width of the boson decay into a SM fermion pair.
Definition: NPEffectiveBS.cpp:258
Model::name
std::string name
The name of the model.
Definition: Model.h:275
StandardModel::Mz
double Mz
The mass of the boson in GeV.
Definition: StandardModel.h:2554
NPEffectiveBS::NPEffectiveBSVars_LFU
static const std::string NPEffectiveBSVars_LFU[NNPEffectiveBSVars_LFU]
A string array containing the labels of the model parameters in NPEffectiveBS with lepton flavour uni...
Definition: NPEffectiveBS.h:297
NPEffectiveBS::deltaGV_f
virtual double deltaGV_f(const Particle f) const
New physics contribution to the neutral-current vector coupling .
Definition: NPEffectiveBS.cpp:303
NPEffectiveBS::cHQ1
double cHQ1
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:483
NPEffectiveBS::setParameter
virtual void setParameter(const std::string name, const double &value)
A method to set the value of a parameter of the model.
Definition: NPEffectiveBS.cpp:83
NPEffectiveBS::cHE3
double cHE3
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:488
NPEffectiveBS::cL1L2
double cL1L2
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:469
NPEffectiveBS::cHU2
double cHU2
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:490
NPEffectiveBS::cHU3
double cHU3
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:491
NPEffectiveBS::NNPEffectiveBSVars_LFU_QFU
static const int NNPEffectiveBSVars_LFU_QFU
The number of the model parameters in NPEffectiveBS with lepton and quark flavour universalities.
Definition: NPEffectiveBS.h:285
NPEffectiveBS::DeltaGF
virtual double DeltaGF() const
New physics contribution to the Fermi constant.
Definition: NPEffectiveBS.cpp:248
NPEffectiveBS::deltaGL_f_tmp
double deltaGL_f_tmp(const Particle f) const
New physics contribution to .
Definition: NPEffectiveBS.cpp:316
NPEffectiveBS::CheckParameters
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...
Definition: NPEffectiveBS.cpp:184
NPEffectiveBS::NNPEffectiveBSVars
static const int NNPEffectiveBSVars
The number of the model parameters in NPEffectiveBS.
Definition: NPEffectiveBS.h:267
NPEffectiveBS::cL2L3
double cL2L3
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:472
NPEffectiveBS::cHQ3p
double cHQ3p
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:479
NPEffectiveBS::cHE1
double cHE1
The dimension-6 operator coefficient .
Definition: NPEffectiveBS.h:486