NPEffectiveBS.cpp
Go to the documentation of this file.
1 /*
2  * Copyright (C) 2013 HEPfit Collaboration
3  * All rights reserved.
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::pair<std::string, boost::reference_wrapper<const double> >("cWB_NP", boost::cref(cWB)));
35  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cH_NP", boost::cref(cH)));
36  if (FlagLeptonUniversal) {
37  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cLL_NP", boost::cref(cL1L1)));
38  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHLp_NP", boost::cref(cHL1p)));
39  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHL_NP", boost::cref(cHL1)));
40  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHE_NP", boost::cref(cHE1)));
41  } else {
42  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cL1L1_NP", boost::cref(cL1L1)));
43  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cL1L2_NP", boost::cref(cL1L2)));
44  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cL1L3_NP", boost::cref(cL1L3)));
45  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cL2L2_NP", boost::cref(cL2L2)));
46  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cL2L3_NP", boost::cref(cL2L3)));
47  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cL3L3_NP", boost::cref(cL3L3)));
48  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHL1p_NP", boost::cref(cHL1p)));
49  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHL2p_NP", boost::cref(cHL2p)));
50  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHL3p_NP", boost::cref(cHL3p)));
51  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHL1_NP", boost::cref(cHL1)));
52  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHL2_NP", boost::cref(cHL2)));
53  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHL3_NP", boost::cref(cHL3)));
54  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHE1_NP", boost::cref(cHE1)));
55  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHE2_NP", boost::cref(cHE2)));
56  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHE3_NP", boost::cref(cHE3)));
57  }
58  if (FlagQuarkUniversal) {
59  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHQp_NP", boost::cref(cHQ1p)));
60  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHQ_NP", boost::cref(cHQ1)));
61  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHU_NP", boost::cref(cHU1)));
62  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHD_NP", boost::cref(cHD1)));
63  } else {
64  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHQ1p_NP", boost::cref(cHQ1p)));
65  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHQ2p_NP", boost::cref(cHQ2p)));
66  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHQ3p_NP", boost::cref(cHQ3p)));
67  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHQ1_NP", boost::cref(cHQ1)));
68  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHQ2_NP", boost::cref(cHQ2)));
69  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHQ3_NP", boost::cref(cHQ3)));
70  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHU1_NP", boost::cref(cHU1)));
71  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHU2_NP", boost::cref(cHU2)));
72  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHU3_NP", boost::cref(cHU3)));
73  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHD1_NP", boost::cref(cHD1)));
74  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHD2_NP", boost::cref(cHD2)));
75  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("cHD3_NP", boost::cref(cHD3)));
76  }
77  ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("Lambda_NP", boost::cref(Lambda_NP)));
78 }
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 {
187  for (int i = 0; i < NNPEffectiveBSVars_LFU_QFU; i++) {
188  if (DPars.find(NPEffectiveBSVars_LFU_QFU[i]) == DPars.end()) {
189  std::cout << "ERROR: Missing mandatory NPEffectiveBS_LFU_QFU parameter "
190  << NPEffectiveBSVars_LFU_QFU[i] << std::endl;
191  return false;
192  }
193  }
194  } else if (FlagLeptonUniversal && !FlagQuarkUniversal) {
195  for (int i = 0; i < NNPEffectiveBSVars_LFU; i++) {
196  if (DPars.find(NPEffectiveBSVars_LFU[i]) == DPars.end()) {
197  std::cout << "ERROR: Missing mandatory NPEffectiveBS_LFU parameter "
198  << NPEffectiveBSVars_LFU[i] << std::endl;
199  return false;
200  }
201  }
202  } else if (!FlagLeptonUniversal && FlagQuarkUniversal) {
203  for (int i = 0; i < NNPEffectiveBSVars_QFU; i++) {
204  if (DPars.find(NPEffectiveBSVars_QFU[i]) == DPars.end()) {
205  std::cout << "ERROR: Missing mandatory NPEffectiveBS_QFU parameter "
206  << NPEffectiveBSVars_QFU[i] << std::endl;
207  return false;
208  }
209  }
210  } else if (!FlagLeptonUniversal && !FlagQuarkUniversal) {
211  for (int i = 0; i < NNPEffectiveBSVars; i++) {
212  if (DPars.find(NPEffectiveBSVars[i]) == DPars.end()) {
213  std::cout << "ERROR: Missing mandatory NPEffectiveBS parameter "
214  << NPEffectiveBSVars[i] << std::endl;
215  return false;
216  }
217  }
218  } else
219  throw std::runtime_error("Error in NPEffectiveBS::CheckParameters()");
220 
221  return (NPbase::CheckParameters(DPars));
222 }
223 
224 
225 
227 
228 double NPEffectiveBS::v() const
229 {
230  //return ( sqrt( (1.0 - (cL1L2 - cHL1p - cHL2p)/sqrt(2.0)/GF/LambdaNP/LambdaNP)
231  // /sqrt(2.0)/GF ) );
232 
233  /* use the tree-level relation */
234  return trueSM.v();
235 }
236 
238 {
239  double GF0 = GF * (1.0 - DeltaGF());
240  double tmp = 4.0 * M_PI * ale / sqrt(2.0) / GF0 / Mz / Mz;
241  return ( Mz / sqrt(2.0) * sqrt(1.0 + sqrt(1.0 - tmp)));
242 }
243 
245 {
246  double ratio = v() * v() / Lambda_NP / Lambda_NP;
247 
248  return ( -(cL1L2 - cHL1p - cHL2p) * ratio);
249 }
250 
251 
253 
254 double NPEffectiveBS::GammaW() const
255 {
256  double Gamma_W = trueSM.GammaW();
257 
258  double alpha = alphaMz();
259  double c2 = trueSM.cW2();
260  double s2 = trueSM.sW2();
261  double ratio = v() * v() / Lambda_NP / Lambda_NP;
262 
263  Gamma_W *= 1.0 - 3.0 * alpha / 4.0 / (c2 - s2)
264  *(obliqueS() - 2.0 * c2 * obliqueT() - (c2 - s2) * obliqueU() / 2.0 / s2)
265  - (1.0 + c2) / 2.0 / (c2 - s2) * DeltaGF()
266  //+ (cHL1p + cHL2p + cHL3p + cHQ1p + cHQ2p) * ratio; // incorrect
267  + 2.0 / 9.0 * (cHL1p + cHL2p + cHL3p + 3.0 * cHQ1p + 3.0 * cHQ2p) * ratio;
268 
269  return Gamma_W;
270 }
271 
272 
274 
276 {
277  double ratio = v() * v() / Lambda_NP / Lambda_NP;
278  double sW_SM = sqrt(trueSM.sW2()); /* This has to be the SM value. */
279  double cW_SM = sqrt(trueSM.cW2()); /* This has to be the SM value. */
280 
281  return ( 4.0 * sW_SM * cW_SM * cWB / alphaMz() * ratio);
282 }
283 
285 {
286  double ratio = v() * v() / Lambda_NP / Lambda_NP;
287 
288  return ( -cH / 2.0 / alphaMz() * ratio);
289 }
290 
292 {
293  return 0.0;
294 }
295 
296 
298 
299 double NPEffectiveBS::deltaGV_f(const Particle f) const
300 {
301  return ( deltaGL_f_tmp(f) + deltaGR_f_tmp(f) + NPbase::deltaGV_f(f));
302 }
303 
304 double NPEffectiveBS::deltaGA_f(const Particle f) const
305 {
306  return ( deltaGL_f_tmp(f) - deltaGR_f_tmp(f) + NPbase::deltaGA_f(f));
307 }
308 
309 
311 
313 {
314  double ratio = v() * v() / Lambda_NP / Lambda_NP;
315  if (f.is("NEUTRINO_1"))
316  return ( (cHL1p - cHL1) / 2.0 * ratio);
317  if (f.is("NEUTRINO_2"))
318  return ( (cHL2p - cHL2) / 2.0 * ratio);
319  if (f.is("NEUTRINO_3"))
320  return ( (cHL3p - cHL3) / 2.0 * ratio);
321  if (f.is("ELECTRON"))
322  return ( -(cHL1p + cHL1) / 2.0 * ratio);
323  if (f.is("MU"))
324  return ( -(cHL2p + cHL2) / 2.0 * ratio);
325  if (f.is("TAU"))
326  return ( -(cHL3p + cHL3) / 2.0 * ratio);
327  if (f.is("UP"))
328  return ( (cHQ1p - cHQ1) / 2.0 * ratio);
329  if (f.is("CHARM"))
330  return ( (cHQ2p - cHQ2) / 2.0 * ratio);
331  if (f.is("TOP"))
332  return 0.0;
333  if (f.is("DOWN"))
334  return ( -(cHQ1p + cHQ1) / 2.0 * ratio);
335  if (f.is("STRANGE"))
336  return ( -(cHQ2p + cHQ2) / 2.0 * ratio);
337  if (f.is("BOTTOM"))
338  return ( -(cHQ3p + cHQ3) / 2.0 * ratio);
339  throw std::runtime_error("Error in NPEffectiveBS::deltaGL_f()");
340 }
341 
343 {
344  double ratio = v() * v() / Lambda_NP / Lambda_NP;
345  if (f.is("NEUTRINO_1"))
346  return ( 0.);
347  if (f.is("NEUTRINO_2"))
348  return ( 0.);
349  if (f.is("NEUTRINO_3"))
350  return ( 0.);
351  if (f.is("ELECTRON"))
352  return ( -cHE1 / 2.0 * ratio);
353  if (f.is("MU"))
354  return ( -cHE2 / 2.0 * ratio);
355  if (f.is("TAU"))
356  return ( -cHE3 / 2.0 * ratio);
357  if (f.is("UP"))
358  return ( -cHU1 / 2.0 * ratio);
359  if (f.is("CHARM"))
360  return ( -cHU2 / 2.0 * ratio);
361  if (f.is("TOP"))
362  return 0.0;
363  if (f.is("DOWN"))
364  return ( -cHD1 / 2.0 * ratio);
365  if (f.is("STRANGE"))
366  return ( -cHD2 / 2.0 * ratio);
367  if (f.is("BOTTOM"))
368  return ( -cHD3 / 2.0 * ratio);
369  throw std::runtime_error("Error in NPEffectiveBS::deltaGL_f()");
370 }
371 
372 
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.
Definition: NPbase.h:66
static const int NNPEffectiveBSVars
The number of the model parameters in NPEffectiveBS.
A class for particles.
Definition: Particle.h:26
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 .
Definition: NPbase.cpp:60
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 .
StandardModel trueSM
Definition: NPbase.h:543
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
Definition: Particle.cpp:23
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
Definition: Model.h:200
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 .
Definition: NPbase.cpp:82
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 .