a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
GeneralTHDM.cpp
Go to the documentation of this file.
1 /*
2  * Copyright (C) 2016 HEPfit Collaboration
3  *
4  *
5  * For the licensing terms see doc/COPYING.
6  */
7 
8 #include "GeneralTHDM.h"
9 #include "GeneralTHDMcache.h"
10 
11 std::string GeneralTHDM::GeneralTHDMvars[NGeneralTHDMvars] = {"logtb", "mHp2", "mH2sq", "mH3sq", "alpha1", "alpha2", "alpha3", "Relambda5", "Imlambda5", "Relambda6", "Relambda7",
12 "Nu_11r", "Nu_11i", "Nu_12r", "Nu_12i", "Nu_13r", "Nu_13i",
13 "Nu_21r", "Nu_21i", "Nu_22r", "Nu_22i", "Nu_23r", "Nu_23i",
14 "Nu_31r", "Nu_31i", "Nu_32r", "Nu_32i", "Nu_33r", "Nu_33i",
15 "Nd_11r", "Nd_11i", "Nd_12r", "Nd_12i", "Nd_13r", "Nd_13i",
16 "Nd_21r", "Nd_21i", "Nd_22r", "Nd_22i", "Nd_23r", "Nd_23i",
17 "Nd_31r", "Nd_31i", "Nd_32r", "Nd_32i", "Nd_33r", "Nd_33i",
18 "Nl_11r", "Nl_11i", "Nl_12r", "Nl_12i", "Nl_13r", "Nl_13i",
19 "Nl_21r", "Nl_21i", "Nl_22r", "Nl_22i", "Nl_23r", "Nl_23i",
20 "Nl_31r", "Nl_31i", "Nl_32r", "Nl_32i", "Nl_33r", "Nl_33i",
21 "Q_GTHDM","RpepsGTHDM","NLOuniscaleGTHDM"};
22 
23 GeneralTHDM::GeneralTHDM() : NPbase(), GTHDMM(*this) {
24 
25  SMM.setObj((StandardModelMatching&) GTHDMM.getObj());
26  ModelParamMap.insert(std::make_pair("logtb", std::cref(logtb)));
27  ModelParamMap.insert(std::make_pair("mHp1", std::cref(mHp1)));
28  ModelParamMap.insert(std::make_pair("mH21", std::cref(mH21)));
29  ModelParamMap.insert(std::make_pair("mH31", std::cref(mH31)));
30  ModelParamMap.insert(std::make_pair("mHp2", std::cref(mHp2)));
31  ModelParamMap.insert(std::make_pair("mH2sq", std::cref(mH2sq)));
32  ModelParamMap.insert(std::make_pair("mH3sq", std::cref(mH3sq)));
33  ModelParamMap.insert(std::make_pair("alpha1", std::cref(alpha1)));
34  ModelParamMap.insert(std::make_pair("alpha2", std::cref(alpha2)));
35  ModelParamMap.insert(std::make_pair("alpha3", std::cref(alpha3)));
36  ModelParamMap.insert(std::make_pair("Relambda5", std::cref(Relambda5)));
37  ModelParamMap.insert(std::make_pair("Imlambda5", std::cref(Imlambda5)));
38  ModelParamMap.insert(std::make_pair("Relambda6", std::cref(Relambda6)));
39  ModelParamMap.insert(std::make_pair("Relambda7", std::cref(Relambda7)));
40 
41  ModelParamMap.insert(std::make_pair("yu1R_GTHDM", std::cref(yu1R_GTHDM)));
42  ModelParamMap.insert(std::make_pair("yd1R_GTHDM", std::cref(yd1R_GTHDM)));
43  ModelParamMap.insert(std::make_pair("yl1R_GTHDM", std::cref(yl1R_GTHDM)));
44 
45 
46 
47  ModelParamMap.insert(std::make_pair("Nu_11r", std::cref(Nu_11r)));
48  ModelParamMap.insert(std::make_pair("Nu_11i", std::cref(Nu_11i)));
49  ModelParamMap.insert(std::make_pair("Nu_12r", std::cref(Nu_12r)));
50  ModelParamMap.insert(std::make_pair("Nu_12i", std::cref(Nu_12i)));
51  ModelParamMap.insert(std::make_pair("Nu_13r", std::cref(Nu_13r)));
52  ModelParamMap.insert(std::make_pair("Nu_13i", std::cref(Nu_13i)));
53  ModelParamMap.insert(std::make_pair("Nu_21r", std::cref(Nu_21r)));
54  ModelParamMap.insert(std::make_pair("Nu_21i", std::cref(Nu_21i)));
55  ModelParamMap.insert(std::make_pair("Nu_22r", std::cref(Nu_22r)));
56  ModelParamMap.insert(std::make_pair("Nu_22i", std::cref(Nu_22i)));
57  ModelParamMap.insert(std::make_pair("Nu_23r", std::cref(Nu_23r)));
58  ModelParamMap.insert(std::make_pair("Nu_23i", std::cref(Nu_23i)));
59  ModelParamMap.insert(std::make_pair("Nu_31r", std::cref(Nu_31r)));
60  ModelParamMap.insert(std::make_pair("Nu_31i", std::cref(Nu_31i)));
61  ModelParamMap.insert(std::make_pair("Nu_32r", std::cref(Nu_32r)));
62  ModelParamMap.insert(std::make_pair("Nu_32i", std::cref(Nu_32i)));
63  ModelParamMap.insert(std::make_pair("Nu_33r", std::cref(Nu_33r)));
64  ModelParamMap.insert(std::make_pair("Nu_33i", std::cref(Nu_33i)));
65  ModelParamMap.insert(std::make_pair("Nd_11r", std::cref(Nd_11r)));
66  ModelParamMap.insert(std::make_pair("Nd_11i", std::cref(Nd_11i)));
67  ModelParamMap.insert(std::make_pair("Nd_12r", std::cref(Nd_12r)));
68  ModelParamMap.insert(std::make_pair("Nd_12i", std::cref(Nd_12i)));
69  ModelParamMap.insert(std::make_pair("Nd_13r", std::cref(Nd_13r)));
70  ModelParamMap.insert(std::make_pair("Nd_13i", std::cref(Nd_13i)));
71  ModelParamMap.insert(std::make_pair("Nd_21r", std::cref(Nd_21r)));
72  ModelParamMap.insert(std::make_pair("Nd_21i", std::cref(Nd_21i)));
73  ModelParamMap.insert(std::make_pair("Nd_22r", std::cref(Nd_22r)));
74  ModelParamMap.insert(std::make_pair("Nd_22i", std::cref(Nd_22i)));
75  ModelParamMap.insert(std::make_pair("Nd_23r", std::cref(Nd_23r)));
76  ModelParamMap.insert(std::make_pair("Nd_23i", std::cref(Nd_23i)));
77  ModelParamMap.insert(std::make_pair("Nd_31r", std::cref(Nd_31r)));
78  ModelParamMap.insert(std::make_pair("Nd_31i", std::cref(Nd_31i)));
79  ModelParamMap.insert(std::make_pair("Nd_32r", std::cref(Nd_32r)));
80  ModelParamMap.insert(std::make_pair("Nd_32i", std::cref(Nd_32i)));
81  ModelParamMap.insert(std::make_pair("Nd_33r", std::cref(Nd_33r)));
82  ModelParamMap.insert(std::make_pair("Nd_33i", std::cref(Nd_33i)));
83  ModelParamMap.insert(std::make_pair("Nl_11r", std::cref(Nl_11r)));
84  ModelParamMap.insert(std::make_pair("Nl_11i", std::cref(Nl_11i)));
85  ModelParamMap.insert(std::make_pair("Nl_12r", std::cref(Nl_12r)));
86  ModelParamMap.insert(std::make_pair("Nl_12i", std::cref(Nl_12i)));
87  ModelParamMap.insert(std::make_pair("Nl_13r", std::cref(Nl_13r)));
88  ModelParamMap.insert(std::make_pair("Nl_13i", std::cref(Nl_13i)));
89  ModelParamMap.insert(std::make_pair("Nl_21r", std::cref(Nl_21r)));
90  ModelParamMap.insert(std::make_pair("Nl_21i", std::cref(Nl_21i)));
91  ModelParamMap.insert(std::make_pair("Nl_22r", std::cref(Nl_22r)));
92  ModelParamMap.insert(std::make_pair("Nl_22i", std::cref(Nl_22i)));
93  ModelParamMap.insert(std::make_pair("Nl_23r", std::cref(Nl_23r)));
94  ModelParamMap.insert(std::make_pair("Nl_23i", std::cref(Nl_23i)));
95  ModelParamMap.insert(std::make_pair("Nl_31r", std::cref(Nl_31r)));
96  ModelParamMap.insert(std::make_pair("Nl_31i", std::cref(Nl_31i)));
97  ModelParamMap.insert(std::make_pair("Nl_32r", std::cref(Nl_32r)));
98  ModelParamMap.insert(std::make_pair("Nl_32i", std::cref(Nl_32i)));
99  ModelParamMap.insert(std::make_pair("Nl_33r", std::cref(Nl_33r)));
100  ModelParamMap.insert(std::make_pair("Nl_33i", std::cref(Nl_33i)));
101  ModelParamMap.insert(std::make_pair("Q_GTHDM", std::cref(Q_GTHDM)));
102  ModelParamMap.insert(std::make_pair("RpepsGTHDM", std::cref(RpepsGTHDM)));
103  ModelParamMap.insert(std::make_pair("NLOuniscaleGTHDM", std::cref(NLOuniscaleGTHDM)));
104  flag_use_sq_masses=true;
105  flag_sigma=true;
106 
107 
108 }
109 
110 GeneralTHDM::~GeneralTHDM(){
111  if (IsModelInitialized()) {
112  if (myGTHDMcache != NULL) delete(myGTHDMcache);
113  }
114 }
115 
117 // Initialization
118 
119 bool GeneralTHDM::InitializeModel()
120 {
121  myGTHDMcache = new GeneralTHDMcache(*this);
122  setModelInitialized(StandardModel::InitializeModel());
123  setModelGeneralTHDM();
124  return(true);
125 }
126 
127 bool GeneralTHDM::Init(const std::map<std::string, double>& DPars) {
128  return(StandardModel::Init(DPars));
129 }
130 
131 bool GeneralTHDM::PreUpdate()
132 {
133  if(!StandardModel::PreUpdate()) return (false);
134 
135  return (true);
136 }
137 
138 bool GeneralTHDM::Update(const std::map<std::string, double>& DPars) {
139 
140  if(!PreUpdate()) return (false);
141 
142  UpdateError = false;
143 
144  for (std::map<std::string, double>::const_iterator it = DPars.begin(); it != DPars.end(); it++)
145  setParameter(it->first, it->second);
146 
147  if (UpdateError) return (false);
148 
149  if(!PostUpdate()) return (false);
150 
151  return (true);
152 }
153 
154 bool GeneralTHDM::PostUpdate()
155 {
156  if(!StandardModel::PostUpdate()) return (false);
157 
158  myGTHDMcache->updateCache();
159 
160  return (true);
161 }
162 
163 void GeneralTHDM::setParameter(const std::string name, const double& value){
164 
165  if(name.compare("logtb") == 0) {
166  logtb = value;
167  tanb = pow(10.,logtb);
168  if(tanb > 0.) {
169  sinb = tanb / sqrt(1. + tanb*tanb);
170  cosb = 1. / sqrt(1. + tanb*tanb);
171  }
172  else {
173  throw std::runtime_error("error in GeneralTHDM::SetParameter, tanb < 0!");
174  }
175  }
176  else if(name.compare("mH21") == 0 && !flag_use_sq_masses)
177  mH21 = value;
178  else if(name.compare("mH31") == 0 && !flag_use_sq_masses)
179  mH31 = value;
180  else if(name.compare("mHp1") == 0 && !flag_use_sq_masses)
181  mHp1 = value;
182  else if(name.compare("mH2sq") == 0 && flag_use_sq_masses)
183  mH2sq = value;
184  else if(name.compare("mH3sq") == 0 && flag_use_sq_masses)
185  mH3sq = value;
186  else if(name.compare("mHp2") == 0 && flag_use_sq_masses)
187  mHp2 = value;
188  else if(name.compare("alpha1") == 0) {
189  alpha1 = value;
190  cosalpha1 = cos(alpha1);
191  sinalpha1 = sin(alpha1);
192  // std::cout<<"alpha1 before ordering = "<<alpha1<<std::endl;
193 
194  }
195  else if(name.compare("alpha2") == 0) {
196  alpha2 = value;
197  cosalpha2 = cos(alpha2);
198  sinalpha2 = sin(alpha2);
199  //std::cout<<"alpha2 before ordering = "<<alpha2<<std::endl;
200  }
201  else if(name.compare("alpha3") == 0) {
202  alpha3 = value;
203  cosalpha3 = cos(alpha3);
204  sinalpha3 = sin(alpha3);
205  //std::cout<<"alpha3 before ordering = "<<alpha3<<std::endl;
206  }
207  else if(name.compare("Relambda5") == 0)
208  Relambda5 = value;
209  else if(name.compare("Imlambda5") == 0)
210  Imlambda5 = value;
211  else if(name.compare("Relambda6") == 0)
212  Relambda6 = value;
213  else if(name.compare("Relambda7") == 0)
214  Relambda7 = value;
215  else if(name.compare("Nu_11r") == 0 && flag_sigma)
216  Nu_11r = value;
217  else if(name.compare("yu1R_GTHDM") == 0 && !flag_sigma){
218  if(alpha1==0)
219  yu1R_GTHDM=1;
220  else
221  yu1R_GTHDM = value;
222  }
223  else if(name.compare("Nu_11i") == 0)
224  Nu_11i = value;
225  else if(name.compare("Nu_12r") == 0)
226  Nu_12r = value;
227  else if(name.compare("Nu_12i") == 0)
228  Nu_12i = value;
229  else if(name.compare("Nu_13r") == 0)
230  Nu_13r = value;
231  else if(name.compare("Nu_13i") == 0)
232  Nu_13i = value;
233  else if(name.compare("Nu_21r") == 0)
234  Nu_21r = value;
235  else if(name.compare("Nu_21i") == 0)
236  Nu_21i = value;
237  else if(name.compare("Nu_22r") == 0)
238  Nu_22r = value;
239  else if(name.compare("Nu_22i") == 0)
240  Nu_22i = value;
241  else if(name.compare("Nu_23r") == 0)
242  Nu_23r = value;
243  else if(name.compare("Nu_23i") == 0)
244  Nu_23i = value;
245  else if(name.compare("Nu_31r") == 0)
246  Nu_31r = value;
247  else if(name.compare("Nu_31i") == 0)
248  Nu_31i = value;
249  else if(name.compare("Nu_32r") == 0)
250  Nu_32r = value;
251  else if(name.compare("Nu_32i") == 0)
252  Nu_32i = value;
253  else if(name.compare("Nu_33r") == 0)
254  Nu_33r = value;
255  else if(name.compare("Nu_33i") == 0)
256  Nu_33i = value;
257  else if(name.compare("Nd_11r") == 0 && flag_sigma)
258  Nd_11r = value;
259  else if(name.compare("yd1R_GTHDM") == 0 && !flag_sigma)
260  if(alpha1==0)
261  yd1R_GTHDM=1;
262  else
263  yd1R_GTHDM = value;
264  else if(name.compare("Nd_11i") == 0)
265  Nd_11i = value;
266  else if(name.compare("Nd_12r") == 0)
267  Nd_12r = value;
268  else if(name.compare("Nd_12i") == 0)
269  Nd_12i = value;
270  else if(name.compare("Nd_13r") == 0)
271  Nd_13r = value;
272  else if(name.compare("Nd_13i") == 0)
273  Nd_13i = value;
274  else if(name.compare("Nd_21r") == 0)
275  Nd_21r = value;
276  else if(name.compare("Nd_21i") == 0)
277  Nd_21i = value;
278  else if(name.compare("Nd_22r") == 0)
279  Nd_22r = value;
280  else if(name.compare("Nd_22i") == 0)
281  Nd_22i = value;
282  else if(name.compare("Nd_23r") == 0)
283  Nd_23r = value;
284  else if(name.compare("Nd_23i") == 0)
285  Nd_23i = value;
286  else if(name.compare("Nd_31r") == 0)
287  Nd_31r = value;
288  else if(name.compare("Nd_31i") == 0)
289  Nd_31i = value;
290  else if(name.compare("Nd_32r") == 0)
291  Nd_32r = value;
292  else if(name.compare("Nd_32i") == 0)
293  Nd_32i = value;
294  else if(name.compare("Nd_33r") == 0)
295  Nd_33r = value;
296  else if(name.compare("Nd_33i") == 0)
297  Nd_33i = value;
298  else if(name.compare("Nl_11r") == 0 && flag_sigma)
299  Nl_11r = value;
300  else if(name.compare("yl1R_GTHDM") == 0 && !flag_sigma)
301  if(alpha1==0)
302  yl1R_GTHDM=1;
303  else
304  yl1R_GTHDM = value;
305  else if(name.compare("Nl_11i") == 0)
306  Nl_11i = value;
307  else if(name.compare("Nl_12r") == 0)
308  Nl_12r = value;
309  else if(name.compare("Nl_12i") == 0)
310  Nl_12i = value;
311  else if(name.compare("Nl_13r") == 0)
312  Nl_13r = value;
313  else if(name.compare("Nl_13i") == 0)
314  Nl_13i = value;
315  else if(name.compare("Nl_21r") == 0)
316  Nl_21r = value;
317  else if(name.compare("Nl_21i") == 0)
318  Nl_21i = value;
319  else if(name.compare("Nl_22r") == 0)
320  Nl_22r = value;
321  else if(name.compare("Nl_22i") == 0)
322  Nl_22i = value;
323  else if(name.compare("Nl_23r") == 0)
324  Nl_23r = value;
325  else if(name.compare("Nl_23i") == 0)
326  Nl_23i = value;
327  else if(name.compare("Nl_31r") == 0)
328  Nl_31r = value;
329  else if(name.compare("Nl_31i") == 0)
330  Nl_31i = value;
331  else if(name.compare("Nl_32r") == 0)
332  Nl_32r = value;
333  else if(name.compare("Nl_32i") == 0)
334  Nl_32i = value;
335  else if(name.compare("Nl_33r") == 0)
336  Nl_33r = value;
337  else if(name.compare("Nl_33i") == 0)
338  Nl_33i = value;
339  else if(name.compare("Q_GTHDM") == 0)
340  Q_GTHDM = value;
341  else if(name.compare("RpepsGTHDM") == 0)
342  RpepsGTHDM = value;
343  else if(name.compare("NLOuniscaleGTHDM") == 0)
344  NLOuniscaleGTHDM = value;
345  else
346  StandardModel::setParameter(name,value);
347 }
348 
349 bool GeneralTHDM::CheckParameters(const std::map<std::string, double>& DPars) {
350  for (int i = 0; i < NGeneralTHDMvars; i++) {
351  if (DPars.find(GeneralTHDMvars[i]) == DPars.end()) {
352  std::cout << "ERROR: missing mandatory GeneralTHDM parameter " << GeneralTHDMvars[i] << std::endl;
353  raiseMissingModelParameterCount();
354  addMissingModelParameter(GeneralTHDMvars[i]);
355  }
356  }
357  return(StandardModel::CheckParameters(DPars));
358 }
359 
361 // Flags
362 
363 bool GeneralTHDM::setFlagStr(const std::string name, const std::string value)
364 {
365  bool res = false;
366  if(name.compare("RGEorder") == 0)
367  {
368  if (checkRGEorder(value))
369  {
370  flag_RGEorder = value;
371  res = true;
372  }
373  else
374  {
375  throw std::runtime_error("GeneralTHDM::setFlagStr(): Invalid flag "
376  + name + "=" + value);
377  }
378  }
379  else
380  {
381  res = StandardModel::setFlagStr(name,value);
382  }
383 
384  return(res);
385 }
386 
387 bool GeneralTHDM::setFlag(const std::string name, const bool value)
388 {
389  bool res = false;
390  if(name.compare("use_sq_masses") == 0) {
391  flag_use_sq_masses = value;
392  res = true;
393  if (!flag_use_sq_masses) {
394  GeneralTHDMvars[std::distance(GeneralTHDMvars,std::find(GeneralTHDMvars,GeneralTHDMvars+NGeneralTHDMvars,"mH2sq"))] = "mH21";
395  GeneralTHDMvars[std::distance(GeneralTHDMvars,std::find(GeneralTHDMvars,GeneralTHDMvars+NGeneralTHDMvars,"mH3sq"))] = "mH31";
396  GeneralTHDMvars[std::distance(GeneralTHDMvars,std::find(GeneralTHDMvars,GeneralTHDMvars+NGeneralTHDMvars,"mHp2"))] = "mHp1";
397 
398  ModelParamMap.insert(std::make_pair("mH21", std::cref(mH21)));
399  ModelParamMap.insert(std::make_pair("mH31", std::cref(mH31)));
400  ModelParamMap.insert(std::make_pair("mHp1", std::cref(mHp1)));
401 
402  }
403  }
404  else if(name.compare("use_sigma") == 0) {
405  std::cout<<"use_sigma = "<< value << std::endl;
406  flag_sigma = value;
407  res = true;
408  if (!flag_sigma && flag_CPconservation) {
409  GeneralTHDMvars[std::distance(GeneralTHDMvars,std::find(GeneralTHDMvars,GeneralTHDMvars+NGeneralTHDMvars,"Nu_11r"))] = "yu1R_GTHDM";
410  GeneralTHDMvars[std::distance(GeneralTHDMvars,std::find(GeneralTHDMvars,GeneralTHDMvars+NGeneralTHDMvars,"Nd_11r"))] = "yd1R_GTHDM";
411  GeneralTHDMvars[std::distance(GeneralTHDMvars,std::find(GeneralTHDMvars,GeneralTHDMvars+NGeneralTHDMvars,"Nl_11r"))] = "yl1R_GTHDM";
412 
413  ModelParamMap.insert(std::make_pair("yu1R_GTHDM", std::cref(yu1R_GTHDM)));
414  ModelParamMap.insert(std::make_pair("yd1R_GTHDM", std::cref(yd1R_GTHDM)));
415  ModelParamMap.insert(std::make_pair("yl1R_GTHDM", std::cref(yl1R_GTHDM)));
416  }
417  if(!flag_sigma && !flag_CPconservation){
418  throw std::runtime_error("sigma flag different from true only in the CP-conserving");
419  }
420  }
421  else if(name.compare("ATHDMflag") == 0) {
422  std::cout<<"ATHDMflag = "<< value<<std::endl;
423  flag_ATHDM = value;
424  res = true;
425  }
426  else if(name.compare("CPconservation") == 0) {
427  std::cout<<"CPconservation = "<< value<<std::endl;
428  flag_CPconservation = value;
429  res = true;
430  }
431  else
432  res = StandardModel::setFlag(name,value);
433 
434  return(res);
435 }
436 
437 double GeneralTHDM::muggH(const double sqrt_s) const
438 {
439  return getMyGTHDMCache()->rh_gg;
440 }
441 
442 double GeneralTHDM::muVBF(const double sqrt_s) const
443 {
444  return getMyGTHDMCache()->rh_VV;
445 }
446 
447 double GeneralTHDM::mueeWBF(const double sqrt_s) const
448 {
449  return getMyGTHDMCache()->rh_VV;
450 }
451 
452 double GeneralTHDM::muWH(const double sqrt_s) const
453 {
454  return getMyGTHDMCache()->rh_VV;
455 }
456 
457 double GeneralTHDM::muZH(const double sqrt_s) const
458 {
459  return getMyGTHDMCache()->rh_VV;
460 }
461 
462 double GeneralTHDM::mueeZH(const double sqrt_s) const
463 {
464  return getMyGTHDMCache()->rh_VV;
465 }
466 
467 double GeneralTHDM::muVH(const double sqrt_s) const
468 {
469  return getMyGTHDMCache()->rh_VV;
470 }
471 
472 double GeneralTHDM::muVBFpVH(const double sqrt_s) const
473 {
474  return getMyGTHDMCache()->rh_VV;
475 }
476 
477 double GeneralTHDM::muttH(const double sqrt_s) const
478 {
479  return getMyGTHDMCache()->rh_QuQuE + getMyGTHDMCache()->rh_QuQuO/(getMyGTHDMCache()->beta_h_t*getMyGTHDMCache()->beta_h_t);
480 }
481 
482 double GeneralTHDM::computeGammaTotalRatio() const
483 {
484  return getMyGTHDMCache()->sumModBRs;
485 }
486 
487 double GeneralTHDM::GammaTotal() const
488 {
489  return getMyGTHDMCache()->Gamma_h;
490 }
491 
492 double GeneralTHDM::BrHggRatio() const
493 {
494  return getMyGTHDMCache()->rh_gg / computeGammaTotalRatio();
495 }
496 
497 double GeneralTHDM::BrHWWRatio() const
498 {
499  return getMyGTHDMCache()->rh_VV / computeGammaTotalRatio();
500 }
501 
502 double GeneralTHDM::BrHZZRatio() const
503 {
504  return getMyGTHDMCache()->rh_VV / computeGammaTotalRatio();
505 }
506 
507 double GeneralTHDM::BrHZgaRatio() const
508 {
509  return getMyGTHDMCache()->rh_Zga / computeGammaTotalRatio();
510 }
511 
512 double GeneralTHDM::BrHgagaRatio() const
513 {
514  return getMyGTHDMCache()->rh_gaga / computeGammaTotalRatio();
515 }
516 
517 double GeneralTHDM::BrHmumuRatio() const
518 {
519  return getMyGTHDMCache()->rh_QlQlE + getMyGTHDMCache()->rh_QlQlO/(getMyGTHDMCache()->beta_h_mu*getMyGTHDMCache()->beta_h_mu) / computeGammaTotalRatio();
520 }
521 
522 double GeneralTHDM::BrHtautauRatio() const
523 {
524  return getMyGTHDMCache()->rh_QlQlE + getMyGTHDMCache()->rh_QlQlO/(getMyGTHDMCache()->beta_h_tau*getMyGTHDMCache()->beta_h_tau) / computeGammaTotalRatio();
525 }
526 
527 double GeneralTHDM::BrHccRatio() const
528 {
529  return getMyGTHDMCache()->rh_QuQuE + getMyGTHDMCache()->rh_QuQuO/(getMyGTHDMCache()->beta_h_c*getMyGTHDMCache()->beta_h_c) / computeGammaTotalRatio();
530 }
531 
532 double GeneralTHDM::BrHbbRatio() const
533 {
534  return getMyGTHDMCache()->rh_QdQdE + getMyGTHDMCache()->rh_QdQdO/(getMyGTHDMCache()->beta_h_b*getMyGTHDMCache()->beta_h_b) / computeGammaTotalRatio();
535 }
536 
537 double GeneralTHDM::muggHgaga(const double sqrt_s) const
538 {
539  return getMyGTHDMCache()->rh_gg * getMyGTHDMCache()->rh_gaga / computeGammaTotalRatio();
540 }
541 
542 double GeneralTHDM::muVBFHgaga(const double sqrt_s) const
543 {
544  return getMyGTHDMCache()->rh_VV * getMyGTHDMCache()->rh_gaga / computeGammaTotalRatio();
545 }
546 
547 double GeneralTHDM::muVHgaga(const double sqrt_s) const
548 {
549  return getMyGTHDMCache()->rh_VV * getMyGTHDMCache()->rh_gaga / computeGammaTotalRatio();
550 }
551 
552 double GeneralTHDM::muttHgaga(const double sqrt_s) const
553 {
554  return (getMyGTHDMCache()->rh_QuQuE + getMyGTHDMCache()->rh_QuQuO/(getMyGTHDMCache()->beta_h_t*getMyGTHDMCache()->beta_h_t)) * getMyGTHDMCache()->rh_gaga / computeGammaTotalRatio();
555 }
556 
557 double GeneralTHDM::muggHZZ(const double sqrt_s) const
558 {
559  return getMyGTHDMCache()->rh_gg * getMyGTHDMCache()->rh_VV / computeGammaTotalRatio();
560 }
561 
562 double GeneralTHDM::muVBFHZZ(const double sqrt_s) const
563 {
564  return getMyGTHDMCache()->rh_VV * getMyGTHDMCache()->rh_VV / computeGammaTotalRatio();
565 }
566 
567 double GeneralTHDM::muVHZZ(const double sqrt_s) const
568 {
569  return getMyGTHDMCache()->rh_VV * getMyGTHDMCache()->rh_VV / computeGammaTotalRatio();
570 }
571 
572 double GeneralTHDM::muttHZZ(const double sqrt_s) const
573 {
574  return (getMyGTHDMCache()->rh_QuQuE + getMyGTHDMCache()->rh_QuQuO/(getMyGTHDMCache()->beta_h_t*getMyGTHDMCache()->beta_h_t)) * getMyGTHDMCache()->rh_VV / computeGammaTotalRatio();
575 }
576 
577 double GeneralTHDM::muggHWW(const double sqrt_s) const
578 {
579  return getMyGTHDMCache()->rh_gg * getMyGTHDMCache()->rh_VV / computeGammaTotalRatio();
580 }
581 
582 double GeneralTHDM::muVBFHWW(const double sqrt_s) const
583 {
584  return getMyGTHDMCache()->rh_VV * getMyGTHDMCache()->rh_VV / computeGammaTotalRatio();
585 }
586 
587 double GeneralTHDM::muVHWW(const double sqrt_s) const
588 {
589  return getMyGTHDMCache()->rh_VV * getMyGTHDMCache()->rh_VV / computeGammaTotalRatio();
590 }
591 
592 double GeneralTHDM::muttHWW(const double sqrt_s) const
593 {
594  return (getMyGTHDMCache()->rh_QuQuE + getMyGTHDMCache()->rh_QuQuO/(getMyGTHDMCache()->beta_h_t*getMyGTHDMCache()->beta_h_t)) * getMyGTHDMCache()->rh_VV / computeGammaTotalRatio();
595 }
596 
597 double GeneralTHDM::muggHtautau(const double sqrt_s) const
598 {
599  return getMyGTHDMCache()->rh_gg * (getMyGTHDMCache()->rh_QlQlE + getMyGTHDMCache()->rh_QlQlO/(getMyGTHDMCache()->beta_h_tau*getMyGTHDMCache()->beta_h_tau)) / computeGammaTotalRatio();
600 }
601 
602 double GeneralTHDM::muVBFHtautau(const double sqrt_s) const
603 {
604  return getMyGTHDMCache()->rh_VV * (getMyGTHDMCache()->rh_QlQlE + getMyGTHDMCache()->rh_QlQlO/(getMyGTHDMCache()->beta_h_tau*getMyGTHDMCache()->beta_h_tau)) / computeGammaTotalRatio();
605 }
606 
607 double GeneralTHDM::muVHtautau(const double sqrt_s) const
608 {
609  return getMyGTHDMCache()->rh_VV * (getMyGTHDMCache()->rh_QlQlE + getMyGTHDMCache()->rh_QlQlO/(getMyGTHDMCache()->beta_h_tau*getMyGTHDMCache()->beta_h_tau)) / computeGammaTotalRatio();
610 }
611 
612 double GeneralTHDM::muttHtautau(const double sqrt_s) const
613 {
614  return (getMyGTHDMCache()->rh_QuQuE + getMyGTHDMCache()->rh_QuQuO/(getMyGTHDMCache()->beta_h_t*getMyGTHDMCache()->beta_h_t))
615  * (getMyGTHDMCache()->rh_QlQlE + getMyGTHDMCache()->rh_QlQlO/(getMyGTHDMCache()->beta_h_tau*getMyGTHDMCache()->beta_h_tau)) / computeGammaTotalRatio();
616 }
617 
618 double GeneralTHDM::muggHbb(const double sqrt_s) const
619 {
620  return getMyGTHDMCache()->rh_gg * (getMyGTHDMCache()->rh_QdQdE + getMyGTHDMCache()->rh_QdQdO/(getMyGTHDMCache()->beta_h_b*getMyGTHDMCache()->beta_h_b)) / computeGammaTotalRatio();
621 }
622 
623 double GeneralTHDM::muVBFHbb(const double sqrt_s) const
624 {
625  return getMyGTHDMCache()->rh_VV * (getMyGTHDMCache()->rh_QdQdE + getMyGTHDMCache()->rh_QdQdO/(getMyGTHDMCache()->beta_h_b*getMyGTHDMCache()->beta_h_b)) / computeGammaTotalRatio();
626 }
627 
628 double GeneralTHDM::muVHbb(const double sqrt_s) const
629 {
630  return getMyGTHDMCache()->rh_VV * (getMyGTHDMCache()->rh_QdQdE + getMyGTHDMCache()->rh_QdQdO/(getMyGTHDMCache()->beta_h_b*getMyGTHDMCache()->beta_h_b)) / computeGammaTotalRatio();
631 }
632 
633 double GeneralTHDM::muttHbb(const double sqrt_s) const
634 {
635  return (getMyGTHDMCache()->rh_QuQuE + getMyGTHDMCache()->rh_QuQuO/(getMyGTHDMCache()->beta_h_t*getMyGTHDMCache()->beta_h_t))
636  * (getMyGTHDMCache()->rh_QdQdE + getMyGTHDMCache()->rh_QdQdO/(getMyGTHDMCache()->beta_h_b*getMyGTHDMCache()->beta_h_b)) / computeGammaTotalRatio();
637 }
638 
639 double GeneralTHDM::muppHmumu(const double sqrt_s) const
640 {
641  if(sqrt_s==8)
642  {
643  return (0.872 * getMyGTHDMCache()->rh_gg + 0.122 * getMyGTHDMCache()->rh_VV + 0.006 * (getMyGTHDMCache()->rh_QuQuE + getMyGTHDMCache()->rh_QuQuO/(getMyGTHDMCache()->beta_h_t*getMyGTHDMCache()->beta_h_t)))
644  * (getMyGTHDMCache()->rh_QlQlE + getMyGTHDMCache()->rh_QlQlO/(getMyGTHDMCache()->beta_h_mu*getMyGTHDMCache()->beta_h_mu)) / computeGammaTotalRatio();
645  }
646  if(sqrt_s==13)
647  {
648  return (0.871 * getMyGTHDMCache()->rh_gg + 0.119 * getMyGTHDMCache()->rh_VV + 0.010 * (getMyGTHDMCache()->rh_QuQuE + getMyGTHDMCache()->rh_QuQuO/(getMyGTHDMCache()->beta_h_t*getMyGTHDMCache()->beta_h_t)))
649  * (getMyGTHDMCache()->rh_QlQlE + getMyGTHDMCache()->rh_QlQlO/(getMyGTHDMCache()->beta_h_mu*getMyGTHDMCache()->beta_h_mu)) / computeGammaTotalRatio();
650  }
651  else
652  {
653  throw std::runtime_error("The observable muppHmumu is only defined for 8 or 13 TeV.");
654  }
655 }
656 
657 double GeneralTHDM::muppHZga(const double sqrt_s) const
658 {
659  if(sqrt_s==8)
660  {
661  return (0.872 * getMyGTHDMCache()->rh_gg + 0.122 * getMyGTHDMCache()->rh_VV + 0.006 * (getMyGTHDMCache()->rh_QuQuE + getMyGTHDMCache()->rh_QuQuO/(getMyGTHDMCache()->beta_h_t*getMyGTHDMCache()->beta_h_t))) * getMyGTHDMCache()->rh_Zga / computeGammaTotalRatio();
662  }
663  if(sqrt_s==13)
664  {
665  return (0.871 * getMyGTHDMCache()->rh_gg + 0.119 * getMyGTHDMCache()->rh_VV + 0.010 * (getMyGTHDMCache()->rh_QuQuE + getMyGTHDMCache()->rh_QuQuO/(getMyGTHDMCache()->beta_h_t*getMyGTHDMCache()->beta_h_t))) * getMyGTHDMCache()->rh_Zga / computeGammaTotalRatio();
666  }
667  else
668  {
669  throw std::runtime_error("The observable muppHZga is only defined for 8 or 13 TeV.");
670  }
671 }
672 
673 double GeneralTHDM::Mw() const{
674  double MZ = StandardModel::Mz;
675  return ( MZ / sqrt(2.0) * sqrt(1.0 + sqrt(1.0 - 4.0 * M_PI * StandardModel::ale / (sqrt(2.0) * StandardModel::GF * MZ* MZ))));
676 }
GeneralTHDM::BrHZgaRatio
virtual double BrHZgaRatio() const
Definition: GeneralTHDM.cpp:507
GeneralTHDM::checkRGEorder
bool checkRGEorder(const std::string RGEorder) const
A method to check if the RGE order name in string form is valid.
Definition: GeneralTHDM.h:1295
GeneralTHDM::mH2sq
double mH2sq
Definition: GeneralTHDM.h:1310
gslpp::cos
complex cos(const complex &z)
Definition: gslpp_complex.cpp:429
GeneralTHDM::Nu_12r
double Nu_12r
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nu_21r
double Nu_21r
Definition: GeneralTHDM.h:1310
GeneralTHDM::muVBFHbb
virtual double muVBFHbb(const double sqrt_s) const
Definition: GeneralTHDM.cpp:623
GeneralTHDM::muttHWW
virtual double muttHWW(const double sqrt_s) const
Definition: GeneralTHDM.cpp:592
GeneralTHDM::flag_use_sq_masses
bool flag_use_sq_masses
Definition: GeneralTHDM.h:1323
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:231
GeneralTHDM::Nd_11i
double Nd_11i
Definition: GeneralTHDM.h:1310
GeneralTHDM::Update
virtual bool Update(const std::map< std::string, double > &DPars)
The update method for GeneralTHDM.
Definition: GeneralTHDM.cpp:138
GeneralTHDM::Nd_12i
double Nd_12i
Definition: GeneralTHDM.h:1310
yu1R_GTHDM
Coupling of the SM-Higgs to up quarks real part.
Definition: GeneralTHDMLightHiggs.h:738
GeneralTHDM::Nl_22i
double Nl_22i
Definition: GeneralTHDM.h:1310
GeneralTHDM::mueeWBF
virtual double mueeWBF(const double sqrt_s) const
Definition: GeneralTHDM.cpp:447
GeneralTHDM::BrHccRatio
virtual double BrHccRatio() const
Definition: GeneralTHDM.cpp:527
GeneralTHDM::GeneralTHDM
GeneralTHDM()
GeneralTHDM constructor.
Definition: GeneralTHDM.cpp:23
GeneralTHDM::Nl_21r
double Nl_21r
Definition: GeneralTHDM.h:1310
GeneralTHDM::muppHmumu
virtual double muppHmumu(const double sqrt_s) const
Definition: GeneralTHDM.cpp:639
GeneralTHDM::muttHtautau
virtual double muttHtautau(const double sqrt_s) const
Definition: GeneralTHDM.cpp:612
GeneralTHDM::cosalpha2
double cosalpha2
Definition: GeneralTHDM.h:1310
GeneralTHDM::muggHtautau
virtual double muggHtautau(const double sqrt_s) const
Definition: GeneralTHDM.cpp:597
GeneralTHDM::BrHgagaRatio
virtual double BrHgagaRatio() const
Definition: GeneralTHDM.cpp:512
GeneralTHDM::muVH
virtual double muVH(const double sqrt_s) const
Definition: GeneralTHDM.cpp:467
GeneralTHDM::muVBF
virtual double muVBF(const double sqrt_s) const
Definition: GeneralTHDM.cpp:442
GeneralTHDM::Nd_21r
double Nd_21r
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nl_33r
double Nl_33r
Definition: GeneralTHDM.h:1310
GeneralTHDM::GeneralTHDMvars
static std::string GeneralTHDMvars[NGeneralTHDMvars]
Definition: GeneralTHDM.h:467
gslpp::sin
complex sin(const complex &z)
Definition: gslpp_complex.cpp:420
GeneralTHDM::alpha1
double alpha1
Definition: GeneralTHDM.h:1310
Matching::getObj
T & getObj()
Definition: Matching.h:14
yd1R_GTHDM
Coupling of the SM-Higgs to down quarks real part.
Definition: GeneralTHDMLightHiggs.h:759
GeneralTHDM::flag_RGEorder
std::string flag_RGEorder
Definition: GeneralTHDM.h:1324
GeneralTHDM::Nl_11i
double Nl_11i
Definition: GeneralTHDM.h:1310
GeneralTHDM::Relambda5
double Relambda5
Definition: GeneralTHDM.h:1310
GeneralTHDM::getMyGTHDMCache
GeneralTHDMcache * getMyGTHDMCache() const
Definition: GeneralTHDM.h:564
GeneralTHDM::setParameter
virtual void setParameter(const std::string, const double &)
A method to set the value of a parameter of GeneralTHDM.
Definition: GeneralTHDM.cpp:163
StandardModel::GF
double GF
The Fermi constant in .
Definition: StandardModel.h:2511
GeneralTHDM::Mw
virtual double Mw() const
Definition: GeneralTHDM.cpp:673
GeneralTHDM::Nd_32i
double Nd_32i
Definition: GeneralTHDM.h:1310
Model::addMissingModelParameter
void addMissingModelParameter(const std::string &missingParameterName)
Definition: Model.h:232
GeneralTHDM::Nd_22i
double Nd_22i
Definition: GeneralTHDM.h:1310
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:313
GeneralTHDM::muVBFHtautau
virtual double muVBFHtautau(const double sqrt_s) const
Definition: GeneralTHDM.cpp:602
GeneralTHDM::Nu_21i
double Nu_21i
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nl_12i
double Nl_12i
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nd_11r
double Nd_11r
Definition: GeneralTHDM.h:1310
GeneralTHDM::Relambda7
double Relambda7
Definition: GeneralTHDM.h:1310
GeneralTHDM::muVHZZ
virtual double muVHZZ(const double sqrt_s) const
Definition: GeneralTHDM.cpp:567
GeneralTHDM::alpha3
double alpha3
Definition: GeneralTHDM.h:1310
GeneralTHDM::mHp2
double mHp2
Definition: GeneralTHDM.h:1310
Matching::setObj
void setObj(T &obji)
Definition: Matching.h:15
Model::setModelGeneralTHDM
void setModelGeneralTHDM()
Definition: Model.h:194
GeneralTHDM::muppHZga
virtual double muppHZga(const double sqrt_s) const
Definition: GeneralTHDM.cpp:657
GeneralTHDM::Nd_12r
double Nd_12r
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nl_21i
double Nl_21i
Definition: GeneralTHDM.h:1310
GeneralTHDM::BrHmumuRatio
virtual double BrHmumuRatio() const
Definition: GeneralTHDM.cpp:517
Model::UpdateError
bool UpdateError
A boolean set to false if update is successful.
Definition: Model.h:254
StandardModel::SMM
Matching< StandardModelMatching, StandardModel > SMM
An object of type Matching.
Definition: StandardModel.h:2506
GeneralTHDM::mH21
double mH21
Definition: GeneralTHDM.h:1310
GeneralTHDM::muggHWW
virtual double muggHWW(const double sqrt_s) const
Definition: GeneralTHDM.cpp:577
GeneralTHDM::setFlag
virtual bool setFlag(const std::string, const bool)
A method to set a flag of GeneralTHDM.
Definition: GeneralTHDM.cpp:387
GeneralTHDM::muVHgaga
virtual double muVHgaga(const double sqrt_s) const
Definition: GeneralTHDM.cpp:547
GeneralTHDM::Nu_33i
double Nu_33i
Definition: GeneralTHDM.h:1310
GeneralTHDM::CheckParameters
virtual bool CheckParameters(const std::map< std::string, double > &DPars)
A method to check if all the mandatory parameters for GeneralTHDM have been provided in model initial...
Definition: GeneralTHDM.cpp:349
GeneralTHDM::cosalpha1
double cosalpha1
Definition: GeneralTHDM.h:1310
GeneralTHDM::~GeneralTHDM
~GeneralTHDM()
GeneralTHDM destructor.
Definition: GeneralTHDM.cpp:110
GeneralTHDM::PreUpdate
virtual bool PreUpdate()
The pre-update method for GeneralTHDM.
Definition: GeneralTHDM.cpp:131
GeneralTHDM::Nl_31i
double Nl_31i
Definition: GeneralTHDM.h:1310
StandardModel::ale
double ale
The fine-structure constant .
Definition: StandardModel.h:2512
StandardModel::setFlag
virtual bool setFlag(const std::string name, const bool value)
A method to set a flag of StandardModel.
Definition: StandardModel.cpp:378
Model::ModelParamMap
std::map< std::string, std::reference_wrapper< const double > > ModelParamMap
Definition: Model.h:262
GeneralTHDM::Nl_33i
double Nl_33i
Definition: GeneralTHDM.h:1310
GeneralTHDM::sinalpha2
double sinalpha2
Definition: GeneralTHDM.h:1310
yl1R_GTHDM
Coupling of the SM-Higgs to leptons real part.
Definition: GeneralTHDMLightHiggs.h:780
GeneralTHDM::Nu_31r
double Nu_31r
Definition: GeneralTHDM.h:1310
GeneralTHDM::muVBFpVH
virtual double muVBFpVH(const double sqrt_s) const
Definition: GeneralTHDM.cpp:472
GeneralTHDM::myGTHDMcache
GeneralTHDMcache * myGTHDMcache
Definition: GeneralTHDM.h:1308
GeneralTHDM::BrHtautauRatio
virtual double BrHtautauRatio() const
Definition: GeneralTHDM.cpp:522
GeneralTHDM::Nd_13i
double Nd_13i
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nl_22r
double Nl_22r
Definition: GeneralTHDM.h:1310
StandardModel::Init
virtual bool Init(const std::map< std::string, double > &DPars)
A method to initialize the model parameters.
Definition: StandardModel.cpp:159
GeneralTHDM::GTHDMM
Matching< GeneralTHDMMatching, GeneralTHDM > GTHDMM
An object of type Matching.
Definition: GeneralTHDM.h:1304
StandardModelMatching
A class for the matching in the Standard Model.
Definition: StandardModelMatching.h:26
GeneralTHDM::Nl_13r
double Nl_13r
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nd_22r
double Nd_22r
Definition: GeneralTHDM.h:1310
NPbase
The auxiliary base model class for other model classes.
Definition: NPbase.h:66
GeneralTHDM::Nu_22i
double Nu_22i
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nl_23r
double Nl_23r
Definition: GeneralTHDM.h:1310
GeneralTHDM::flag_ATHDM
bool flag_ATHDM
Definition: GeneralTHDM.h:1323
StandardModel::PreUpdate
virtual bool PreUpdate()
The pre-update method for StandardModel.
Definition: StandardModel.cpp:172
GeneralTHDM::Nu_32i
double Nu_32i
Definition: GeneralTHDM.h:1310
GeneralTHDM::flag_CPconservation
bool flag_CPconservation
Definition: GeneralTHDM.h:1323
GeneralTHDM::Nu_33r
double Nu_33r
Definition: GeneralTHDM.h:1310
GeneralTHDM::logtb
double logtb
Definition: GeneralTHDM.h:1310
GeneralTHDM::setFlagStr
virtual bool setFlagStr(const std::string name, const std::string value)
A method to set a string flag of GeneralTHDM.
Definition: GeneralTHDM.cpp:363
gslpp::pow
complex pow(const complex &z1, const complex &z2)
Definition: gslpp_complex.cpp:395
GeneralTHDM::mueeZH
virtual double mueeZH(const double sqrt_s) const
Definition: GeneralTHDM.cpp:462
Model::raiseMissingModelParameterCount
void raiseMissingModelParameterCount()
Definition: Model.h:242
gslpp::sqrt
complex sqrt(const complex &z)
Definition: gslpp_complex.cpp:385
GeneralTHDM.h
GeneralTHDM::RpepsGTHDM
double RpepsGTHDM
Definition: GeneralTHDM.h:1310
GeneralTHDM::muggHZZ
virtual double muggHZZ(const double sqrt_s) const
Definition: GeneralTHDM.cpp:557
GeneralTHDM::flag_sigma
bool flag_sigma
Definition: GeneralTHDM.h:1323
GeneralTHDM::Nu_23r
double Nu_23r
Definition: GeneralTHDM.h:1310
GeneralTHDM::computeGammaTotalRatio
virtual double computeGammaTotalRatio() const
Definition: GeneralTHDM.cpp:482
StandardModel::PostUpdate
virtual bool PostUpdate()
The post-update method for StandardModel.
Definition: StandardModel.cpp:199
GeneralTHDM::Imlambda5
double Imlambda5
Definition: GeneralTHDM.h:1310
GeneralTHDMcache
Definition: GeneralTHDMcache.h:22
GeneralTHDM::BrHbbRatio
virtual double BrHbbRatio() const
Definition: GeneralTHDM.cpp:532
GeneralTHDM::sinalpha1
double sinalpha1
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nd_31i
double Nd_31i
Definition: GeneralTHDM.h:1310
StandardModel::setFlagStr
virtual bool setFlagStr(const std::string name, const std::string value)
A method to set a flag of StandardModel.
Definition: StandardModel.cpp:418
Model::IsModelInitialized
bool IsModelInitialized() const
A method to check if the model is initialized.
Definition: Model.h:136
GeneralTHDM::alpha2
double alpha2
Definition: GeneralTHDM.h:1310
GeneralTHDM::muttHgaga
virtual double muttHgaga(const double sqrt_s) const
Definition: GeneralTHDM.cpp:552
GeneralTHDM::NLOuniscaleGTHDM
double NLOuniscaleGTHDM
Definition: GeneralTHDM.h:1310
GeneralTHDM::sinalpha3
double sinalpha3
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nu_11i
double Nu_11i
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nl_32i
double Nl_32i
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nu_13i
double Nu_13i
Definition: GeneralTHDM.h:1310
GeneralTHDM::sinb
double sinb
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nl_12r
double Nl_12r
Definition: GeneralTHDM.h:1310
GeneralTHDMcache.h
GeneralTHDM::muVBFHZZ
virtual double muVBFHZZ(const double sqrt_s) const
Definition: GeneralTHDM.cpp:562
GeneralTHDM::muVBFHgaga
virtual double muVBFHgaga(const double sqrt_s) const
Definition: GeneralTHDM.cpp:542
GeneralTHDM::Nu_23i
double Nu_23i
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nu_22r
double Nu_22r
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nu_31i
double Nu_31i
Definition: GeneralTHDM.h:1310
GeneralTHDM::muggH
virtual double muggH(const double sqrt_s) const
Definition: GeneralTHDM.cpp:437
GeneralTHDM::muVHbb
virtual double muVHbb(const double sqrt_s) const
Definition: GeneralTHDM.cpp:628
GeneralTHDM::Nd_32r
double Nd_32r
Definition: GeneralTHDM.h:1310
GeneralTHDM::muVHtautau
virtual double muVHtautau(const double sqrt_s) const
Definition: GeneralTHDM.cpp:607
StandardModel::InitializeModel
virtual bool InitializeModel()
A method to initialize the model.
Definition: StandardModel.cpp:140
GeneralTHDM::Nl_11r
double Nl_11r
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nu_11r
double Nu_11r
Definition: GeneralTHDM.h:1310
GeneralTHDM::BrHWWRatio
virtual double BrHWWRatio() const
Definition: GeneralTHDM.cpp:497
GeneralTHDM::Nu_32r
double Nu_32r
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nd_33r
double Nd_33r
Definition: GeneralTHDM.h:1310
GeneralTHDM::BrHggRatio
virtual double BrHggRatio() const
Definition: GeneralTHDM.cpp:492
GeneralTHDM::Nl_23i
double Nl_23i
Definition: GeneralTHDM.h:1310
GeneralTHDM::muggHgaga
virtual double muggHgaga(const double sqrt_s) const
Definition: GeneralTHDM.cpp:537
GeneralTHDM::tanb
double tanb
Definition: GeneralTHDM.h:1310
GeneralTHDM::mH31
double mH31
Definition: GeneralTHDM.h:1310
GeneralTHDM::muWH
virtual double muWH(const double sqrt_s) const
Definition: GeneralTHDM.cpp:452
GeneralTHDM::BrHZZRatio
virtual double BrHZZRatio() const
Definition: GeneralTHDM.cpp:502
GeneralTHDM::Nd_23r
double Nd_23r
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nl_13i
double Nl_13i
Definition: GeneralTHDM.h:1310
GeneralTHDM::Q_GTHDM
double Q_GTHDM
Definition: GeneralTHDM.h:1310
GeneralTHDM::mH3sq
double mH3sq
Definition: GeneralTHDM.h:1310
GeneralTHDM::InitializeModel
virtual bool InitializeModel()
A method to initialize the model.
Definition: GeneralTHDM.cpp:119
GeneralTHDM::Nd_31r
double Nd_31r
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nd_13r
double Nd_13r
Definition: GeneralTHDM.h:1310
Model::name
std::string name
The name of the model.
Definition: Model.h:267
StandardModel::Mz
double Mz
The mass of the boson in GeV.
Definition: StandardModel.h:2510
GeneralTHDM::muVBFHWW
virtual double muVBFHWW(const double sqrt_s) const
Definition: GeneralTHDM.cpp:582
GeneralTHDM::Nd_23i
double Nd_23i
Definition: GeneralTHDM.h:1310
GeneralTHDM::muttH
virtual double muttH(const double sqrt_s) const
Definition: GeneralTHDM.cpp:477
GeneralTHDM::Nl_31r
double Nl_31r
Definition: GeneralTHDM.h:1310
GeneralTHDM::GammaTotal
virtual double GammaTotal() const
Definition: GeneralTHDM.cpp:487
GeneralTHDM::muggHbb
virtual double muggHbb(const double sqrt_s) const
Definition: GeneralTHDM.cpp:618
GeneralTHDMcache::updateCache
double updateCache()
Definition: GeneralTHDMcache.cpp:7114
GeneralTHDM::muVHWW
virtual double muVHWW(const double sqrt_s) const
Definition: GeneralTHDM.cpp:587
GeneralTHDM::Relambda6
double Relambda6
Definition: GeneralTHDM.h:1310
GeneralTHDM::muZH
virtual double muZH(const double sqrt_s) const
Definition: GeneralTHDM.cpp:457
GeneralTHDM::muttHZZ
virtual double muttHZZ(const double sqrt_s) const
Definition: GeneralTHDM.cpp:572
GeneralTHDM::NGeneralTHDMvars
static const int NGeneralTHDMvars
Definition: GeneralTHDM.h:466
GeneralTHDM::muttHbb
virtual double muttHbb(const double sqrt_s) const
Definition: GeneralTHDM.cpp:633
GeneralTHDM::Nd_33i
double Nd_33i
Definition: GeneralTHDM.h:1310
GeneralTHDM::cosb
double cosb
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nu_12i
double Nu_12i
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nu_13r
double Nu_13r
Definition: GeneralTHDM.h:1310
GeneralTHDM::Init
virtual bool Init(const std::map< std::string, double > &DPars)
Initializes the GeneralTHDM parameters found in the argument.
Definition: GeneralTHDM.cpp:127
GeneralTHDM::mHp1
double mHp1
Definition: GeneralTHDM.h:1310
GeneralTHDM::cosalpha3
double cosalpha3
Definition: GeneralTHDM.h:1310
GeneralTHDM::Nd_21i
double Nd_21i
Definition: GeneralTHDM.h:1310
Model::setModelInitialized
void setModelInitialized(bool ModelInitialized)
A set method to fix the failure or success of the initialization of the model.
Definition: Model.h:145
GeneralTHDM::Nl_32r
double Nl_32r
Definition: GeneralTHDM.h:1310
GeneralTHDM::PostUpdate
virtual bool PostUpdate()
The post-update method for GeneralTHDM.
Definition: GeneralTHDM.cpp:154