a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
HiggsChiral.cpp
Go to the documentation of this file.
1 /*
2  * Copyright (C) 2017 HEPfit Collaboration
3  *
4  *
5  * For the licensing terms see doc/COPYING.
6  */
7 
8 #include "HiggsChiral.h"
9 
10 const std::string HiggsChiral::HChiralvars[NHChiralvars] = {
11  "cv", "ct", "cb", "cc", "ctau", "cmu", "cg", "cga", "cZga", "obsZgaLimitATLAS13", "obsZgaLimitCMS13", "obsZgaLimitATLAS", "obsZgaLimitCMS", "expZgaLimitATLAS13", "expZgaLimitCMS13", "expZgaLimitATLAS", "expZgaLimitCMS"
12 };
13 
15 : NPbase()
16 {
17  FlagUniversalcf = false;
18 
19  ModelParamMap.insert(std::make_pair("cv", std::cref(cv)));
20  ModelParamMap.insert(std::make_pair("ct", std::cref(ct)));
21  ModelParamMap.insert(std::make_pair("cb", std::cref(cb)));
22  ModelParamMap.insert(std::make_pair("cc", std::cref(cc)));
23  ModelParamMap.insert(std::make_pair("ctau", std::cref(ctau)));
24  ModelParamMap.insert(std::make_pair("cmu", std::cref(cmu)));
25  ModelParamMap.insert(std::make_pair("cg", std::cref(cg)));
26  ModelParamMap.insert(std::make_pair("cga", std::cref(cga)));
27  ModelParamMap.insert(std::make_pair("cZga", std::cref(cZga)));
28  ModelParamMap.insert(std::make_pair("obsZgaLimitATLAS13", std::cref(obsZgaLimitATLAS13)));
29  ModelParamMap.insert(std::make_pair("obsZgaLimitCMS13", std::cref(obsZgaLimitCMS13)));
30  ModelParamMap.insert(std::make_pair("obsZgaLimitATLAS", std::cref(obsZgaLimitATLAS)));
31  ModelParamMap.insert(std::make_pair("obsZgaLimitCMS", std::cref(obsZgaLimitCMS)));
32  ModelParamMap.insert(std::make_pair("expZgaLimitATLAS13", std::cref(expZgaLimitATLAS13)));
33  ModelParamMap.insert(std::make_pair("expZgaLimitCMS13", std::cref(expZgaLimitCMS13)));
34  ModelParamMap.insert(std::make_pair("expZgaLimitATLAS", std::cref(expZgaLimitATLAS)));
35  ModelParamMap.insert(std::make_pair("expZgaLimitCMS", std::cref(expZgaLimitCMS)));
36 
37  loopComputed = false;
38  cg_loop = 0.;
39  cga_loop = 0.;
40  cZga_loop = 0.;
41 }
42 
44 {
45  if (!NPbase::PostUpdate()) return (false);
46 
47 // Check first the flags that control the values of the parameters
48 
49  if (FlagUniversalcf) {
50 // Assign to all cf the value of ct
51  cb = ct;
52  cc = ct;
53  ctau = ct;
54  cmu = ct;
55  }
56 
57 // Then the add the values of the loops to the cache
58 
59  loopComputed = false;
60 
61  cg_loop = computecg();
64 
65  loopComputed = true;
66 
67  return (true);
68 }
69 
70 void HiggsChiral::setParameter(const std::string name, const double& value)
71 {
72  if (name.compare("cv") == 0)
73  cv = value;
74  else if (name.compare("ct") == 0)
75  ct = value;
76  else if (name.compare("cb") == 0)
77  cb = value;
78  else if (name.compare("cc") == 0)
79  cc = value;
80  else if (name.compare("ctau") == 0)
81  ctau = value;
82  else if (name.compare("cmu") == 0)
83  cmu = value;
84  else if (name.compare("cg") == 0)
85  cg = value;
86  else if (name.compare("cga") == 0)
87  cga = value;
88  else if (name.compare("cZga") == 0)
89  cZga = value;
90  else if (name.compare("obsZgaLimitATLAS13") == 0)
91  obsZgaLimitATLAS13 = value;
92  else if (name.compare("obsZgaLimitCMS13") == 0)
93  obsZgaLimitCMS13 = value;
94  else if (name.compare("obsZgaLimitATLAS") == 0)
95  obsZgaLimitATLAS = value;
96  else if (name.compare("obsZgaLimitCMS") == 0)
97  obsZgaLimitCMS = value;
98  else if (name.compare("expZgaLimitATLAS13") == 0)
99  expZgaLimitATLAS13 = value;
100  else if (name.compare("expZgaLimitCMS13") == 0)
101  expZgaLimitCMS13 = value;
102  else if (name.compare("expZgaLimitATLAS") == 0)
103  expZgaLimitATLAS = value;
104  else if (name.compare("expZgaLimitCMS") == 0)
105  expZgaLimitCMS = value;
106  else
107  NPbase::setParameter(name, value);
108 }
109 
110 bool HiggsChiral::CheckParameters(const std::map<std::string, double>& DPars)
111 {
112  for (int i = 0; i < NHChiralvars; i++) {
113  if (DPars.find(HChiralvars[i]) == DPars.end()) {
114  std::cout << "ERROR: missing mandatory HiggsChiral parameter " << HChiralvars[i] << std::endl;
117  }
118  }
119  return (NPbase::CheckParameters(DPars));
120 }
121 
122 bool HiggsChiral::setFlag(const std::string name, const bool value)
123 {
124  bool res = false;
125  if (name.compare("Universalcf") == 0) {
126  FlagUniversalcf = value;
127  res = true;
128  } else
129  res = NPbase::setFlag(name, value);
130 
131  return (res);
132 }
133 
135 
136 double HiggsChiral::obliqueS() const
137 {
138  double Lambda;
139  if (fabs(1.0 - cv * cv) < pow(10.0, -32.0))
140  Lambda = pow(10.0, 19.0);
141  else
142  Lambda = 4.0 * M_PI * v() / sqrt(fabs(1.0 - cv * cv));
143 
144  return ( 1.0 / 12.0 / M_PI * (1.0 - cv * cv) * log(Lambda * Lambda / mHl / mHl));
145 }
146 
147 double HiggsChiral::obliqueT() const
148 {
149  double Lambda;
150  double cW2_SM = trueSM.cW2();
151  if (fabs(1.0 - cv * cv) < pow(10.0, -32.0))
152  Lambda = pow(10.0, 19.0);
153  else
154  Lambda = 4.0 * M_PI * v() / sqrt(fabs(1.0 - cv * cv));
155 
156  return ( -3.0 / 16.0 / M_PI / cW2_SM * (1.0 - cv * cv) * log(Lambda * Lambda / mHl / mHl));
157 }
158 
159 double HiggsChiral::obliqueU() const
160 {
161  return 0.0;
162 }
163 
165 
166 double HiggsChiral::muggH(const double sqrt_s) const
167 {
168 // return Gammagg() / (trueSM.computeBrHtogg()*trueSM.computeGammaHTotal());
169 
170  return computecg() * computecg();
171 }
172 
173 double HiggsChiral::muVBF(const double sqrt_s) const
174 {
175  return computecV() * computecV();
176 }
177 
178 double HiggsChiral::muVBFgamma(const double sqrt_s) const
179 {
180  return computecV() * computecV();
181 }
182 
183 double HiggsChiral::mueeWBF(const double sqrt_s) const
184 {
185  return computecV() * computecV();
186 }
187 
188 double HiggsChiral::mueeWBFPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const
189 {
190  return computecV() * computecV();
191 }
192 
193 double HiggsChiral::mueeHvv(const double sqrt_s) const
194 {
195  return computecV() * computecV();
196 }
197 
198 double HiggsChiral::mueeHvvPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const
199 {
200  return computecV() * computecV();
201 }
202 
203 double HiggsChiral::mueeZBF(const double sqrt_s) const
204 {
205  return computecV() * computecV();
206 }
207 
208 double HiggsChiral::mueeZBFPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const
209 {
210  return computecV() * computecV();
211 }
212 
213 double HiggsChiral::muWH(const double sqrt_s) const
214 {
215  return computecV() * computecV();
216 }
217 
218 double HiggsChiral::muZH(const double sqrt_s) const
219 {
220  return computecV() * computecV();
221 }
222 
223 double HiggsChiral::mueeZH(const double sqrt_s) const
224 {
225  return computecV() * computecV();
226 }
227 
228 double HiggsChiral::mueeZllH(const double sqrt_s) const
229 {
230  return computecV() * computecV();
231 }
232 
233 double HiggsChiral::mueeZqqH(const double sqrt_s) const
234 {
235  return computecV() * computecV();
236 }
237 
238 double HiggsChiral::mueeZHPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const
239 {
240  return computecV() * computecV();
241 }
242 
243 double HiggsChiral::mueeZllHPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const
244 {
245  return computecV() * computecV();
246 }
247 
248 double HiggsChiral::mueeZqqHPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const
249 {
250  return computecV() * computecV();
251 }
252 
253 double HiggsChiral::muVH(const double sqrt_s) const
254 {
255  return computecV() * computecV();
256 }
257 
258 double HiggsChiral::muVBFpVH(const double sqrt_s) const
259 {
260  return computecV() * computecV();
261 }
262 
263 double HiggsChiral::muttH(const double sqrt_s) const
264 {
265  return computect() * computect();
266 }
267 
268 double HiggsChiral::mutHq(const double sqrt_s) const
269 {
270  double mu = 1.0;
271  double Kt = computect();
272  double Kw = computecV();
273 
274 // From ATLAS parametrization of the K-framework (at 14 TeV?)
275  mu = 2.63 * Kt * Kt + 3.58 * Kw * Kw - 5.21 * Kt * Kw;
276 
277  return mu;
278 }
279 
280 double HiggsChiral::muggHpttH(const double sqrt_s) const
281 {
282  double sigmaggH_SM = trueSM.computeSigmaggH(sqrt_s);
283  double sigmattH_SM = trueSM.computeSigmattH(sqrt_s);
284 
285  double sigmaggH = muggH(sqrt_s) * sigmaggH_SM;
286  double sigmattH = muttH(sqrt_s) * sigmattH_SM;
287 
288  return (sigmaggH + sigmattH) / (sigmaggH_SM + sigmattH_SM);
289 }
290 
291 double HiggsChiral::mueettH(const double sqrt_s) const
292 {
293  return computect() * computect();
294 }
295 
296 double HiggsChiral::mueettHPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const
297 {
298  return computect() * computect();
299 }
300 
301 double HiggsChiral::UpperLimitZgammaA13(const double sqrt_s) const
302 {
303  double sigmaggH_SM = trueSM.computeSigmaggH(sqrt_s);
304  double sigmattH_SM = trueSM.computeSigmattH(sqrt_s);
305  double sigmaWH_SM = trueSM.computeSigmaWH(sqrt_s);
306  double sigmaZH_SM = trueSM.computeSigmaZH(sqrt_s);
307  double sigmaWF_SM = trueSM.computeSigmaWF(sqrt_s);
308  double sigmaZF_SM = trueSM.computeSigmaZF(sqrt_s);
309  double sigmaZWF_SM = trueSM.computeSigmaZWF(sqrt_s);
310  double sigmaVBF_SM = sigmaWF_SM + sigmaZF_SM + sigmaZWF_SM;
311 
312  double sigmaggH = muggH(sqrt_s) * sigmaggH_SM;
313  double sigmattH = muttH(sqrt_s) * sigmattH_SM;
314  double sigmaWH = muWH(sqrt_s) * sigmaWH_SM;
315  double sigmaZH = muZH(sqrt_s) * sigmaZH_SM;
316  double sigmaVBF = muVBF(sqrt_s) * sigmaVBF_SM;
317 
318  double muppH=((sigmaggH + sigmattH + sigmaWH + sigmaZH + sigmaVBF) / (sigmaggH_SM + sigmattH_SM + sigmaWH_SM + sigmaZH_SM + sigmaVBF_SM));
319 // double BrHZgaRatio=;
320  double nftos=1.95996398454;
321 
322  return (1.0+(muppH*BrHZgaRatio()-getobsZgaLimitATLAS13())/getexpZgaLimitATLAS13() ) * nftos;
323 }
324 
325 double HiggsChiral::UpperLimitZgammaC13(const double sqrt_s) const
326 {
327  double sigmaggH_SM = trueSM.computeSigmaggH(sqrt_s);
328  double sigmattH_SM = trueSM.computeSigmattH(sqrt_s);
329  double sigmaWH_SM = trueSM.computeSigmaWH(sqrt_s);
330  double sigmaZH_SM = trueSM.computeSigmaZH(sqrt_s);
331  double sigmaWF_SM = trueSM.computeSigmaWF(sqrt_s);
332  double sigmaZF_SM = trueSM.computeSigmaZF(sqrt_s);
333  double sigmaZWF_SM = trueSM.computeSigmaZWF(sqrt_s);
334  double sigmaVBF_SM = sigmaWF_SM + sigmaZF_SM + sigmaZWF_SM;
335 
336  double sigmaggH = muggH(sqrt_s) * sigmaggH_SM;
337  double sigmattH = muttH(sqrt_s) * sigmattH_SM;
338  double sigmaWH = muWH(sqrt_s) * sigmaWH_SM;
339  double sigmaZH = muZH(sqrt_s) * sigmaZH_SM;
340  double sigmaVBF = muVBF(sqrt_s) * sigmaVBF_SM;
341 
342  double muppH=((sigmaggH + sigmattH + sigmaWH + sigmaZH + sigmaVBF) / (sigmaggH_SM + sigmattH_SM + sigmaWH_SM + sigmaZH_SM + sigmaVBF_SM));
343 // double BrHZgaRatio=;
344  double nftos=1.95996398454;
345 
346  return (1.0+(muppH*BrHZgaRatio()-getobsZgaLimitCMS13())/getexpZgaLimitCMS13() ) * nftos;
347 }
348 
349 double HiggsChiral::UpperLimitZgammaA(const double sqrt_s) const
350 {
351  double sigmaggH_SM = trueSM.computeSigmaggH(sqrt_s);
352  double sigmattH_SM = trueSM.computeSigmattH(sqrt_s);
353  double sigmaWH_SM = trueSM.computeSigmaWH(sqrt_s);
354  double sigmaZH_SM = trueSM.computeSigmaZH(sqrt_s);
355  double sigmaWF_SM = trueSM.computeSigmaWF(sqrt_s);
356  double sigmaZF_SM = trueSM.computeSigmaZF(sqrt_s);
357  double sigmaZWF_SM = trueSM.computeSigmaZWF(sqrt_s);
358  double sigmaVBF_SM = sigmaWF_SM + sigmaZF_SM + sigmaZWF_SM;
359 
360  double sigmaggH = muggH(sqrt_s) * sigmaggH_SM;
361  double sigmattH = muttH(sqrt_s) * sigmattH_SM;
362  double sigmaWH = muWH(sqrt_s) * sigmaWH_SM;
363  double sigmaZH = muZH(sqrt_s) * sigmaZH_SM;
364  double sigmaVBF = muVBF(sqrt_s) * sigmaVBF_SM;
365 
366  double muppH=((sigmaggH + sigmattH + sigmaWH + sigmaZH + sigmaVBF) / (sigmaggH_SM + sigmattH_SM + sigmaWH_SM + sigmaZH_SM + sigmaVBF_SM));
367 // double BrHZgaRatio=;
368  double nftos=1.95996398454;
369 
370  return (1.0+(muppH*BrHZgaRatio()-getobsZgaLimitATLAS())/getexpZgaLimitATLAS() ) * nftos;
371 }
372 
373 double HiggsChiral::UpperLimitZgammaC(const double sqrt_s) const
374 {
375  double sigmaggH_SM = trueSM.computeSigmaggH(sqrt_s);
376  double sigmattH_SM = trueSM.computeSigmattH(sqrt_s);
377  double sigmaWH_SM = trueSM.computeSigmaWH(sqrt_s);
378  double sigmaZH_SM = trueSM.computeSigmaZH(sqrt_s);
379  double sigmaWF_SM = trueSM.computeSigmaWF(sqrt_s);
380  double sigmaZF_SM = trueSM.computeSigmaZF(sqrt_s);
381  double sigmaZWF_SM = trueSM.computeSigmaZWF(sqrt_s);
382  double sigmaVBF_SM = sigmaWF_SM + sigmaZF_SM + sigmaZWF_SM;
383 
384  double sigmaggH = muggH(sqrt_s) * sigmaggH_SM;
385  double sigmattH = muttH(sqrt_s) * sigmattH_SM;
386  double sigmaWH = muWH(sqrt_s) * sigmaWH_SM;
387  double sigmaZH = muZH(sqrt_s) * sigmaZH_SM;
388  double sigmaVBF = muVBF(sqrt_s) * sigmaVBF_SM;
389 
390  double muppH=((sigmaggH + sigmattH + sigmaWH + sigmaZH + sigmaVBF) / (sigmaggH_SM + sigmattH_SM + sigmaWH_SM + sigmaZH_SM + sigmaVBF_SM));
391 // double BrHZgaRatio=;
392  double nftos=1.95996398454;
393 
394  return (1.0+(muppH*BrHZgaRatio()-getobsZgaLimitCMS())/getexpZgaLimitCMS() ) * nftos;
395 }
396 
397 double HiggsChiral::cgplusct() const
398 {
399  return cg + ct;
400 }
401 
403 {
404  return cga + ct;
405 }
406 
408 {
409  return cg - cga;
410 }
411 
412 double HiggsChiral::cVpluscb() const
413 {
414  return cv + cb;
415 }
416 
418 {
419  return cv + ctau;
420 }
421 
423 {
424  return cb - cc;
425 }
426 
428 {
429  return cb - ctau;
430 }
431 
433 {
434  return cc - ctau;
435 }
436 
437 double HiggsChiral::Gammagg() const
438 {
439  double Mt=trueSM.getQuarks(QCD::TOP).getMass();
440  double Mb=trueSM.getQuarks(QCD::BOTTOM).getMass();
441  double TAUt=4.0*Mt*Mt/(mHl*mHl);
442  double TAUb=4.0*Mb*Mb/(mHl*mHl);
443 
444  double factor = GF*AlsMz*AlsMz*mHl*mHl*mHl/(sqrt(2.0)*4.0*M_PI*M_PI*M_PI);
445 
446  double ksoftNLO = 1.0 + (AlsMz/M_PI) * ( 73.0/4.0 -7.0*5.0/6.0);
447 
448  return factor * ksoftNLO * ( ( ct * 0.5*TAUt*(1.0+(1.0-TAUt)*f_func(TAUt))
449  +cb * 0.5*TAUb*(1.0+(1.0-TAUb)*f_func(TAUb)) ) * (1.0+11.0*AlsMz/(4.0*M_PI))
450  +cg * 0.5 ).abs2();
451 }
452 
453 double HiggsChiral::GammaWW() const
454 {
456 }
457 
458 double HiggsChiral::GammaZZ() const
459 {
461 }
462 
463 double HiggsChiral::GammaZga() const
464 {
465  double Mt=trueSM.getQuarks(QCD::TOP).getMass();
466  double Mb=trueSM.getQuarks(QCD::BOTTOM).getMass();
468  double MW=trueSM.Mw_tree();
469  double cW2=trueSM.c02();
470  double sW2=1.0-cW2;
471  double TAUt=4.0*Mt*Mt/(mHl*mHl);
472  double TAUb=4.0*Mb*Mb/(mHl*mHl);
473  double TAUtau=4.0*Mtau*Mtau/(mHl*mHl);
474  double TAUw=4.0*MW*MW/(mHl*mHl);
475  double LAMt=4.0*Mt*Mt/(Mz*Mz);
476  double LAMb=4.0*Mb*Mb/(Mz*Mz);
477  double LAMtau=4.0*Mtau*Mtau/(Mz*Mz);
478  double LAMw=4.0*MW*MW/(Mz*Mz);
479 
480  double factor = GF*ale*ale/(sqrt(2.0)*64.0*M_PI*M_PI*M_PI)*pow(mHl*(1.0-Mz*Mz/(mHl*mHl)),3);
481 
482  return factor * ((-ct * 4.0*(0.5-4.0/3.0*sW2)*(Int1(TAUt,LAMt)-Int2(TAUt,LAMt)) * (1.0-AlsMz/M_PI)
483  +cb * 2.0*(-0.5+2.0/3.0*sW2)*(Int1(TAUb,LAMb)-Int2(TAUb,LAMb))
484  +ctau * 2.0*(-0.5+2.0*sW2)*(Int1(TAUtau,LAMtau)-Int2(TAUtau,LAMtau)) )/sqrt(sW2*cW2)
485  -cv * sqrt(cW2/sW2)*(4.0*(3.0-sW2/cW2)*Int2(TAUw,LAMw)
486  +((1.0+2.0/TAUw)*sW2/cW2-(5.0+2.0/TAUw))*Int1(TAUw,LAMw))
487  +cZga).abs2(); /*CHANGE THE cZga COEFFICIENT!*/
488 }
489 
491 {
492  double Mt=trueSM.getQuarks(QCD::TOP).getMass();
493  double Mb=trueSM.getQuarks(QCD::BOTTOM).getMass();
495  double MW=trueSM.Mw_tree();
496  double TAUt=4.0*Mt*Mt/(mHl*mHl);
497  double TAUb=4.0*Mb*Mb/(mHl*mHl);
498  double TAUtau=4.0*Mtau*Mtau/(mHl*mHl);
499  double TAUw=4.0*MW*MW/(mHl*mHl);
500 
501  double factor = GF*ale*ale*mHl*mHl*mHl/(sqrt(2.0)*128.0*M_PI*M_PI*M_PI);
502 
503  return factor * ( ct * (8./3.)*TAUt*(1.+(1.-TAUt)*f_func(TAUt)) * (1.0-AlsMz/M_PI)
504  +cb * (2./3.)*TAUb*(1.+(1.-TAUb)*f_func(TAUb))
505  +ctau * 2.0*TAUtau*(1.+(1.-TAUtau)*f_func(TAUtau))
506  -cv * (2.0+3.0*TAUw+3.0*TAUw*(2.0-TAUw)*f_func(TAUw))
507  +cga * 2.0 ).abs2();
508 }
509 
511 {
513 }
514 
516 {
518 }
519 
520 double HiggsChiral::Gammacc() const
521 {
523 }
524 
525 double HiggsChiral::Gammabb() const
526 {
528 }
529 
531 {
532  return Gammagg() + GammaWW() + GammaZZ() + GammaZga() + Gammagaga()
533  + Gammamumu() + Gammatautau() + Gammacc() + Gammabb();
534 }
535 
537 {
538 // return Gammagg() / GammaTotal() / trueSM.computeBrHtogg();
539 
540  return (computecg() * computecg() / computeGammaTotalRatio());
541 }
542 
544 {
545 // return GammaWW() / GammaTotal() / trueSM.computeBrHtoWW();
546 
547  return (computecV() * computecV() / computeGammaTotalRatio());
548 }
549 
551 {
552 // return GammaWW() / GammaTotal() / trueSM.computeBrHtoWW();
553 
554  return BrHWWRatio();
555 }
556 
558 {
559 // return GammaZZ() / GammaTotal() / trueSM.computeBrHtoZZ();
560 
561  return (computecV() * computecV() / computeGammaTotalRatio());
562 }
563 
565 {
566 // return GammaZZ() / GammaTotal() / trueSM.computeBrHtoZZ();
567 
568  return BrHZZRatio();
569 }
570 
572 {
573 // return GammaZZ() / GammaTotal() / trueSM.computeBrHtoZZ();
574 
575  return BrHZZRatio();
576 }
577 
579 {
580 // return GammaZZ() / GammaTotal() / trueSM.computeBrHtoZZ();
581 
582  return BrHZZRatio();
583 }
584 
586 {
587 // return GammaZZ() / GammaTotal() / trueSM.computeBrHtoZZ();
588 
589  return BrHZZRatio();
590 }
591 
593 {
594 // return GammaZga() / GammaTotal() / trueSM.computeBrHtoZga();
595 
597 }
598 
600 {
601 // return GammaZga() / GammaTotal() / trueSM.computeBrHtoZga();
602 
603  return BrHZgaRatio();
604 }
605 
607 {
608 // return GammaZga() / GammaTotal() / trueSM.computeBrHtoZga();
609 
610  return BrHZgaRatio();
611 }
612 
614 {
615 // return GammaZga() / GammaTotal() / trueSM.computeBrHtoZga();
616 
617  return BrHZgaRatio();
618 }
619 
621 {
622 // return Gammagaga() / GammaTotal() / trueSM.computeBrHtogaga();
623 
625 }
626 
628 {
629 // return Gammamumu() / GammaTotal() / trueSM.computeBrHtomumu();
630 
631  return (computecmu() * computecmu() / computeGammaTotalRatio());
632 }
633 
635 {
636 // return Gammatautau() / GammaTotal() / trueSM.computeBrHtotautau();
637 
639 }
640 
642 {
643 // return Gammacc() / GammaTotal() / trueSM.computeBrHtocc();
644 
645  return (computecc() * computecc() / computeGammaTotalRatio());
646 }
647 
649 {
650 // return Gammabb() / GammaTotal() / trueSM.computeBrHtobb();
651 
652  return (computecb() * computecb() / computeGammaTotalRatio());
653 }
654 
656 {
657 // No extra decays so it scales as the H->ZZ
658 
659  return BrHZZRatio();
660 }
661 
662 double HiggsChiral::muggHgaga(const double sqrt_s) const
663 {
665 }
666 
667 double HiggsChiral::muggHgagaInt(const double sqrt_s) const
668 {
669  double muNWA, GammaRatio, fki;
670  double kt,kb,kc,ks,ku,kd,kta,km,ke;
671  double kw;
672  double kg,kgamma;
673 
674  muNWA = muggH(sqrt_s)*BrHgagaRatio();
675 
676  kt = computect();
677  kb = computecb();
678  kc = computecc();
679  ks = 1.0;
680  ku = 1.0;
681  kd = 1.0;
682  kta = computectau();
683  km = computecmu();
684  ke = 1.0;
685 
686  kw = computecV();
687 
688  kg = (3.0/2.0) * computecg();
689  kgamma = (3.0/2.0) * computecgaga();
690 
691  GammaRatio = computeGammaTotalRatio();
692 
693  fki = 0.000802422 *kb*kb + 0.000312884 *kb*kc + 0.0000182107 *kc*kc +
694  5.94769e-9 *kb*kd + 9.62554e-10 *kc*kd + 6.2785e-15 *kd*kd +
695  5.53251e-10 *kb*ke + 3.51863e-11 *kc*ke + 1.09243e-15 *kd*ke -
696  0.00905016 *kb*kg - 0.00190706 *kc*kg - 5.97591e-9 *kd*kg -
697  6.72288e-10 *ke*kg - 0.0271505 *kb*kgamma - 0.00143029 *kc*kgamma -
698  1.79277e-8 *kd*kgamma - 0.174392 *kg*kgamma + 8.97565e-6 *kb*km +
699  6.21013e-7 *kc*km + 2.37616e-11 *kd*km - 0.0000460022 *kg*km +
700  4.93348e-6 *kb*ks + 8.51176e-7 *kc*ks + 1.29558e-11 *kd*ks +
701  1.16267e-12 *ke*ks - 0.0000123381 *kg*ks - 0.0000370143 *kgamma*ks +
702  2.22544e-8 *km*ks + 6.08665e-9 *ks*ks - 0.0467672 *kb*kt -
703  0.00394193 *kc*kt - 3.08808e-8 *kd*kt - 6.94817e-10 *ke*kt -
704  0.240315 *kg*kt - 0.180236 *kgamma*kt - 0.0000475437 *km*kt -
705  0.0000637578 *ks*kt - 0.248368 *kt*kt + 0.00100168 *kb*kta +
706  0.0000759092 *kc*kta + 3.44671e-9 *kd*kta - 0.00975386 *kg*kta +
707  2.93009e-6 *ks*kta - 0.0100807 *kt*kta + 5.30126e-8 *kb*ku +
708  5.54256e-9 *kc*ku + 1.15815e-13 *kd*ku + 4.05206e-15 *ke*ku -
709  1.03323e-7 *kg*ku - 7.74926e-8 *kgamma*ku +
710  8.62762e-11 *km*ku + 1.17664e-10 *ks*ku - 2.13572e-7 *kt*ku +
711  1.2332e-8 *kta*ku + 3.40922e-13 *ku*ku + 0.169912 *kb*kw +
712  0.00895098 *kc*kw + 1.12194e-7 *kd*kw + 1.09137 *kg*kw +
713  0.000231641 *ks*kw + 1.12795 *kt*kw + 4.8496e-7 *ku*kw;
714 
715  return (muNWA - 0.022 *GammaRatio * fki)/0.978;
716 
717 }
718 
719 double HiggsChiral::muVBFHgaga(const double sqrt_s) const
720 {
722 }
723 
724 double HiggsChiral::muZHgaga(const double sqrt_s) const
725 {
727 }
728 
729 double HiggsChiral::muWHgaga(const double sqrt_s) const
730 {
732 }
733 
734 double HiggsChiral::muVHgaga(const double sqrt_s) const
735 {
737 }
738 
739 double HiggsChiral::muttHgaga(const double sqrt_s) const
740 {
742 }
743 
744 double HiggsChiral::muggHZga(const double sqrt_s) const
745 {
747 }
748 
749 double HiggsChiral::muVBFHZga(const double sqrt_s) const
750 {
752 }
753 
754 double HiggsChiral::muZHZga(const double sqrt_s) const
755 {
757 }
758 
759 double HiggsChiral::muWHZga(const double sqrt_s) const
760 {
762 }
763 
764 double HiggsChiral::muVHZga(const double sqrt_s) const
765 {
767 }
768 
769 double HiggsChiral::muttHZga(const double sqrt_s) const
770 {
772 }
773 
774 double HiggsChiral::muggHZZ(const double sqrt_s) const
775 {
777 }
778 
779 double HiggsChiral::muVBFHZZ(const double sqrt_s) const
780 {
782 }
783 
784 double HiggsChiral::muZHZZ(const double sqrt_s) const
785 {
787 }
788 
789 double HiggsChiral::muWHZZ(const double sqrt_s) const
790 {
792 }
793 
794 double HiggsChiral::muVHZZ(const double sqrt_s) const
795 {
797 }
798 
799 double HiggsChiral::muttHZZ(const double sqrt_s) const
800 {
802 }
803 
804 double HiggsChiral::muggHZZ4l(const double sqrt_s) const
805 {
807 }
808 
809 double HiggsChiral::muVBFHZZ4l(const double sqrt_s) const
810 {
812 }
813 
814 double HiggsChiral::muZHZZ4l(const double sqrt_s) const
815 {
817 }
818 
819 double HiggsChiral::muWHZZ4l(const double sqrt_s) const
820 {
822 }
823 
824 double HiggsChiral::muVHZZ4l(const double sqrt_s) const
825 {
827 }
828 
829 double HiggsChiral::muttHZZ4l(const double sqrt_s) const
830 {
832 }
833 
834 double HiggsChiral::muggHWW(const double sqrt_s) const
835 {
837 }
838 
839 double HiggsChiral::muVBFHWW(const double sqrt_s) const
840 {
842 }
843 
844 double HiggsChiral::muZHWW(const double sqrt_s) const
845 {
847 }
848 
849 double HiggsChiral::muWHWW(const double sqrt_s) const
850 {
852 }
853 
854 double HiggsChiral::muVHWW(const double sqrt_s) const
855 {
857 }
858 
859 double HiggsChiral::muttHWW(const double sqrt_s) const
860 {
862 }
863 
864 double HiggsChiral::muggHWW2l2v(const double sqrt_s) const
865 {
867 }
868 
869 double HiggsChiral::muVBFHWW2l2v(const double sqrt_s) const
870 {
872 }
873 
874 double HiggsChiral::muZHWW2l2v(const double sqrt_s) const
875 {
877 }
878 
879 double HiggsChiral::muWHWW2l2v(const double sqrt_s) const
880 {
882 }
883 
884 double HiggsChiral::muVHWW2l2v(const double sqrt_s) const
885 {
887 }
888 
889 double HiggsChiral::muttHWW2l2v(const double sqrt_s) const
890 {
892 }
893 
894 double HiggsChiral::muggHmumu(const double sqrt_s) const
895 {
897 }
898 
899 double HiggsChiral::muVBFHmumu(const double sqrt_s) const
900 {
902 }
903 
904 double HiggsChiral::muZHmumu(const double sqrt_s) const
905 {
907 }
908 
909 double HiggsChiral::muWHmumu(const double sqrt_s) const
910 {
912 }
913 
914 double HiggsChiral::muVHmumu(const double sqrt_s) const
915 {
917 }
918 
919 double HiggsChiral::muttHmumu(const double sqrt_s) const
920 {
922 }
923 
924 double HiggsChiral::muggHtautau(const double sqrt_s) const
925 {
927 }
928 
929 double HiggsChiral::muVBFHtautau(const double sqrt_s) const
930 {
932 }
933 
934 double HiggsChiral::muZHtautau(const double sqrt_s) const
935 {
937 }
938 
939 double HiggsChiral::muWHtautau(const double sqrt_s) const
940 {
942 }
943 
944 double HiggsChiral::muVHtautau(const double sqrt_s) const
945 {
947 }
948 
949 double HiggsChiral::muttHtautau(const double sqrt_s) const
950 {
952 }
953 
954 double HiggsChiral::muggHbb(const double sqrt_s) const
955 {
957 }
958 
959 double HiggsChiral::muVBFHbb(const double sqrt_s) const
960 {
962 }
963 
964 double HiggsChiral::muZHbb(const double sqrt_s) const
965 {
967 }
968 
969 double HiggsChiral::muWHbb(const double sqrt_s) const
970 {
972 }
973 
974 double HiggsChiral::muVHbb(const double sqrt_s) const
975 {
977 }
978 
979 double HiggsChiral::muttHbb(const double sqrt_s) const
980 {
982 }
983 
984 double HiggsChiral::muppHmumu(const double sqrt_s) const
985 {
986  if(sqrt_s==8)
987  {
988  return (0.872 * computecg() * computecg() + 0.122 * computecV() * computecV() + 0.006 * computect() * computect()) * computecmu() * computecmu() / computeGammaTotalRatio();
989  }
990  if(sqrt_s==13)
991  {
992  return (0.871 * computecg() * computecg() + 0.119 * computecV() * computecV() + 0.010 * computect() * computect()) * computecmu() * computecmu() / computeGammaTotalRatio();
993  }
994  else
995  {
996  throw std::runtime_error("The observable muppHmumu is only defined for 8 or 13 TeV.");
997  }
998 }
999 
1000 
1002 //-----------------------------------------------------------------------------------------
1003 //-- Special Hadron collider signal strengths with separate full TH unc U(prod x decay) ---
1004 //-----------------------------------------------------------------------------------------
1006 
1007 double HiggsChiral::muTHUggHgaga(const double sqrt_s) const
1008 {
1009  return muggH(sqrt_s)*BrHgagaRatio();
1010 }
1011 
1012 double HiggsChiral::muTHUVBFHgaga(const double sqrt_s) const
1013 {
1014  return muVBF(sqrt_s)*BrHgagaRatio();
1015 }
1016 
1017 double HiggsChiral::muTHUZHgaga(const double sqrt_s) const
1018 {
1019  return muZH(sqrt_s)*BrHgagaRatio();
1020 }
1021 
1022 double HiggsChiral::muTHUWHgaga(const double sqrt_s) const
1023 {
1024  return muWH(sqrt_s)*BrHgagaRatio();
1025 }
1026 
1027 double HiggsChiral::muTHUVHgaga(const double sqrt_s) const
1028 {
1029  return muVH(sqrt_s)*BrHgagaRatio();
1030 }
1031 
1032 double HiggsChiral::muTHUttHgaga(const double sqrt_s) const
1033 {
1034  return muttH(sqrt_s)*BrHgagaRatio();
1035 }
1036 
1037 double HiggsChiral::muTHUggHZga(const double sqrt_s) const
1038 {
1039  return muggH(sqrt_s)*BrHZgaRatio();
1040 }
1041 
1042 double HiggsChiral::muTHUVBFHZga(const double sqrt_s) const
1043 {
1044  return muVBF(sqrt_s)*BrHZgaRatio();
1045 }
1046 
1047 double HiggsChiral::muTHUZHZga(const double sqrt_s) const
1048 {
1049  return muZH(sqrt_s)*BrHZgaRatio();
1050 }
1051 
1052 double HiggsChiral::muTHUWHZga(const double sqrt_s) const
1053 {
1054  return muWH(sqrt_s)*BrHZgaRatio();
1055 }
1056 
1057 double HiggsChiral::muTHUVHZga(const double sqrt_s) const
1058 {
1059  return muVH(sqrt_s)*BrHZgaRatio();
1060 }
1061 
1062 double HiggsChiral::muTHUttHZga(const double sqrt_s) const
1063 {
1064  return muttH(sqrt_s)*BrHZgaRatio();
1065 }
1066 
1067 double HiggsChiral::muTHUggHZZ(const double sqrt_s) const
1068 {
1069  return muggH(sqrt_s)*BrHZZRatio();
1070 }
1071 
1072 double HiggsChiral::muTHUVBFHZZ(const double sqrt_s) const
1073 {
1074  return muVBF(sqrt_s)*BrHZZRatio();
1075 }
1076 
1077 double HiggsChiral::muTHUZHZZ(const double sqrt_s) const
1078 {
1079  return muZH(sqrt_s)*BrHZZRatio();
1080 }
1081 
1082 double HiggsChiral::muTHUWHZZ(const double sqrt_s) const
1083 {
1084  return muWH(sqrt_s)*BrHZZRatio();
1085 }
1086 
1087 double HiggsChiral::muTHUVHZZ(const double sqrt_s) const
1088 {
1089  return muVH(sqrt_s)*BrHZZRatio();
1090 }
1091 
1092 double HiggsChiral::muTHUttHZZ(const double sqrt_s) const
1093 {
1094  return muttH(sqrt_s)*BrHZZRatio();
1095 }
1096 
1097 double HiggsChiral::muTHUggHZZ4l(const double sqrt_s) const
1098 {
1099  return muggH(sqrt_s)*BrHZZ4lRatio();
1100 }
1101 
1102 double HiggsChiral::muTHUVBFHZZ4l(const double sqrt_s) const
1103 {
1104  return muVBF(sqrt_s)*BrHZZ4lRatio();
1105 }
1106 
1107 double HiggsChiral::muTHUZHZZ4l(const double sqrt_s) const
1108 {
1109  return muZH(sqrt_s)*BrHZZ4lRatio();
1110 }
1111 
1112 double HiggsChiral::muTHUWHZZ4l(const double sqrt_s) const
1113 {
1114  return muWH(sqrt_s)*BrHZZ4lRatio();
1115 }
1116 
1117 double HiggsChiral::muTHUVHZZ4l(const double sqrt_s) const
1118 {
1119  return muVH(sqrt_s)*BrHZZ4lRatio();
1120 }
1121 
1122 double HiggsChiral::muTHUttHZZ4l(const double sqrt_s) const
1123 {
1124  return muttH(sqrt_s)*BrHZZ4lRatio();
1125 }
1126 
1127 double HiggsChiral::muTHUggHWW(const double sqrt_s) const
1128 {
1129  return muggH(sqrt_s)*BrHWWRatio();
1130 }
1131 
1132 double HiggsChiral::muTHUVBFHWW(const double sqrt_s) const
1133 {
1134  return muVBF(sqrt_s)*BrHWWRatio();
1135 }
1136 
1137 double HiggsChiral::muTHUZHWW(const double sqrt_s) const
1138 {
1139  return muZH(sqrt_s)*BrHWWRatio();
1140 }
1141 
1142 double HiggsChiral::muTHUWHWW(const double sqrt_s) const
1143 {
1144  return muWH(sqrt_s)*BrHWWRatio();
1145 }
1146 
1147 double HiggsChiral::muTHUVHWW(const double sqrt_s) const
1148 {
1149  return muVH(sqrt_s)*BrHWWRatio();
1150 }
1151 
1152 double HiggsChiral::muTHUttHWW(const double sqrt_s) const
1153 {
1154  return muttH(sqrt_s)*BrHWWRatio();
1155 }
1156 
1157 double HiggsChiral::muTHUggHWW2l2v(const double sqrt_s) const
1158 {
1159  return muggH(sqrt_s)*BrHWW2l2vRatio();
1160 }
1161 
1162 double HiggsChiral::muTHUVBFHWW2l2v(const double sqrt_s) const
1163 {
1164  return muVBF(sqrt_s)*BrHWW2l2vRatio();
1165 }
1166 
1167 double HiggsChiral::muTHUZHWW2l2v(const double sqrt_s) const
1168 {
1169  return muZH(sqrt_s)*BrHWW2l2vRatio();
1170 }
1171 
1172 double HiggsChiral::muTHUWHWW2l2v(const double sqrt_s) const
1173 {
1174  return muWH(sqrt_s)*BrHWW2l2vRatio();
1175 }
1176 
1177 double HiggsChiral::muTHUVHWW2l2v(const double sqrt_s) const
1178 {
1179  return muVH(sqrt_s)*BrHWW2l2vRatio();
1180 }
1181 
1182 double HiggsChiral::muTHUttHWW2l2v(const double sqrt_s) const
1183 {
1184  return muttH(sqrt_s)*BrHWW2l2vRatio();
1185 }
1186 
1187 double HiggsChiral::muTHUggHmumu(const double sqrt_s) const
1188 {
1189  return muggH(sqrt_s)*BrHmumuRatio();
1190 }
1191 
1192 double HiggsChiral::muTHUVBFHmumu(const double sqrt_s) const
1193 {
1194  return muVBF(sqrt_s)*BrHmumuRatio();
1195 }
1196 
1197 double HiggsChiral::muTHUZHmumu(const double sqrt_s) const
1198 {
1199  return muZH(sqrt_s)*BrHmumuRatio();
1200 }
1201 
1202 double HiggsChiral::muTHUWHmumu(const double sqrt_s) const
1203 {
1204  return muWH(sqrt_s)*BrHmumuRatio();
1205 }
1206 
1207 double HiggsChiral::muTHUVHmumu(const double sqrt_s) const
1208 {
1209  return muVH(sqrt_s)*BrHmumuRatio();
1210 }
1211 
1212 double HiggsChiral::muTHUttHmumu(const double sqrt_s) const
1213 {
1214  return muttH(sqrt_s)*BrHmumuRatio();
1215 }
1216 
1217 double HiggsChiral::muTHUggHtautau(const double sqrt_s) const
1218 {
1219  return muggH(sqrt_s)*BrHtautauRatio();
1220 }
1221 
1222 double HiggsChiral::muTHUVBFHtautau(const double sqrt_s) const
1223 {
1224  return muVBF(sqrt_s)*BrHtautauRatio();
1225 }
1226 
1227 double HiggsChiral::muTHUZHtautau(const double sqrt_s) const
1228 {
1229  return muZH(sqrt_s)*BrHtautauRatio();
1230 }
1231 
1232 double HiggsChiral::muTHUWHtautau(const double sqrt_s) const
1233 {
1234  return muWH(sqrt_s)*BrHtautauRatio();
1235 }
1236 
1237 double HiggsChiral::muTHUVHtautau(const double sqrt_s) const
1238 {
1239  return muVH(sqrt_s)*BrHtautauRatio();
1240 }
1241 
1242 double HiggsChiral::muTHUttHtautau(const double sqrt_s) const
1243 {
1244  return muttH(sqrt_s)*BrHtautauRatio();
1245 }
1246 
1247 double HiggsChiral::muTHUggHbb(const double sqrt_s) const
1248 {
1249  return muggH(sqrt_s)*BrHbbRatio();
1250 }
1251 
1252 double HiggsChiral::muTHUVBFHbb(const double sqrt_s) const
1253 {
1254  return muVBF(sqrt_s)*BrHbbRatio();
1255 }
1256 
1257 double HiggsChiral::muTHUZHbb(const double sqrt_s) const
1258 {
1259  return muZH(sqrt_s)*BrHbbRatio();
1260 }
1261 
1262 double HiggsChiral::muTHUWHbb(const double sqrt_s) const
1263 {
1264  return muWH(sqrt_s)*BrHbbRatio();
1265 }
1266 
1267 double HiggsChiral::muTHUVHbb(const double sqrt_s) const
1268 {
1269  return muVH(sqrt_s)*BrHbbRatio();
1270 }
1271 
1272 double HiggsChiral::muTHUttHbb(const double sqrt_s) const
1273 {
1274  return muttH(sqrt_s)*BrHbbRatio();
1275 }
1276 
1277 double HiggsChiral::muTHUVBFBRinv(const double sqrt_s) const
1278 {
1279  return muVBF(sqrt_s)*Br_H_inv();
1280 }
1281 
1282 double HiggsChiral::muTHUVBFHinv(const double sqrt_s) const
1283 {
1284  return muVBF(sqrt_s)*BrHtoinvRatio();
1285 }
1286 
1287 double HiggsChiral::muTHUVHBRinv(const double sqrt_s) const
1288 {
1289  return muVH(sqrt_s)*Br_H_inv();
1290 }
1291 
1292 double HiggsChiral::muTHUVHinv(const double sqrt_s) const
1293 {
1294  return muVH(sqrt_s)*BrHtoinvRatio();
1295 }
1296 
1297 double HiggsChiral::muTHUggHZZ4mu(const double sqrt_s) const
1298 {
1299  return muggH(sqrt_s)*BrHZZ4muRatio();
1300 }
1301 
1302 double HiggsChiral::muTHUggHZgamumu(const double sqrt_s) const
1303 {
1304  return muggH(sqrt_s)*BrHZgamumuRatio();
1305 }
1306 
1308 //-----------------------------------------------------------------------------------------
1309 //-- Decays ---
1310 //-----------------------------------------------------------------------------------------
1312 
1313 
1315 {
1316 
1317  return ((computecg() * computecg() * trueSM.computeBrHtogg()
1326  / (trueSM.computeBrHtogg()
1334  + trueSM.computeBrHtobb()));
1335 
1336 // return (GammaTotal() / trueSM.computeGammaHTotal());
1337 }
1338 
1340 
1342 {
1343  if (loopComputed) return cg_loop;
1344 
1345  double Mt=trueSM.getQuarks(QCD::TOP).getMass();
1346  double Mb=trueSM.getQuarks(QCD::BOTTOM).getMass();
1347  double Mc=trueSM.getQuarks(QCD::CHARM).getMass();
1348  double TAUt=4.0*Mt*Mt/(mHl*mHl);
1349  double TAUb=4.0*Mb*Mb/(mHl*mHl);
1350  double TAUc=4.0*Mc*Mc/(mHl*mHl);
1351  double cgEff;
1352 
1353  cgEff = ( ( ct * 0.5*TAUt*(1.0+(1.0-TAUt)*f_func(TAUt))
1354  +cb * 0.5*TAUb*(1.0+(1.0-TAUb)*f_func(TAUb))
1355  +cc * 0.5*TAUc*(1.0+(1.0-TAUc)*f_func(TAUc)) ) * (1.0+11.0*AlsMz/(4.0*M_PI))
1356  +cg * 0.5 ).abs2();
1357 
1358  cgEff = cgEff / ( ( 0.5*TAUt*(1.0+(1.0-TAUt)*f_func(TAUt))
1359  + 0.5*TAUb*(1.0+(1.0-TAUb)*f_func(TAUb))
1360  + 0.5*TAUc*(1.0+(1.0-TAUc)*f_func(TAUc)) ) * (1.0+11.0*AlsMz/(4.0*M_PI))).abs2();
1361 
1362  return (sqrt(cgEff));
1363 }
1364 
1366 {
1367  return cv;
1368 }
1369 
1371 {
1372  if (loopComputed) return cZga_loop;
1373 
1374  double Mt=trueSM.getQuarks(QCD::TOP).getMass();
1375  double Mb=trueSM.getQuarks(QCD::BOTTOM).getMass();
1376  double Mc=trueSM.getQuarks(QCD::CHARM).getMass();
1379  double MW=trueSM.Mw_tree();
1380  double cW2=trueSM.c02();
1381  double sW2=1.0-cW2;
1382  double TAUt=4.0*Mt*Mt/(mHl*mHl);
1383  double TAUb=4.0*Mb*Mb/(mHl*mHl);
1384  double TAUc=4.0*Mc*Mc/(mHl*mHl);
1385  double TAUtau=4.0*Mtau*Mtau/(mHl*mHl);
1386  double TAUmu=4.0*Mmu*Mmu/(mHl*mHl);
1387  double TAUw=4.0*MW*MW/(mHl*mHl);
1388  double LAMt=4.0*Mt*Mt/(Mz*Mz);
1389  double LAMb=4.0*Mb*Mb/(Mz*Mz);
1390  double LAMc=4.0*Mc*Mc/(Mz*Mz);
1391  double LAMtau=4.0*Mtau*Mtau/(Mz*Mz);
1392  double LAMmu=4.0*Mmu*Mmu/(Mz*Mz);
1393  double LAMw=4.0*MW*MW/(Mz*Mz);
1394  double cZgaEff;
1395 
1396  cZgaEff = ((-ct * 4.0*(0.5-4.0/3.0*sW2)*(Int1(TAUt,LAMt)-Int2(TAUt,LAMt)) * (1.0-AlsMz/M_PI)
1397  +cb * 2.0*(-0.5+2.0/3.0*sW2)*(Int1(TAUb,LAMb)-Int2(TAUb,LAMb))
1398  -cc * 4.0*(0.5-4.0/3.0*sW2)*(Int1(TAUc,LAMc)-Int2(TAUc,LAMc))
1399  +ctau * 2.0*(-0.5+2.0*sW2)*(Int1(TAUtau,LAMtau)-Int2(TAUtau,LAMtau))
1400  +cmu * 2.0*(-0.5+2.0*sW2)*(Int1(TAUmu,LAMmu)-Int2(TAUmu,LAMmu)) )/sqrt(sW2*cW2)
1401  -cv * sqrt(cW2/sW2)*(4.0*(3.0-sW2/cW2)*Int2(TAUw,LAMw)
1402  +((1.0+2.0/TAUw)*sW2/cW2-(5.0+2.0/TAUw))*Int1(TAUw,LAMw))
1403  +cZga).abs2();
1404 
1405  cZgaEff = cZgaEff / ((-4.0*(0.5-4.0/3.0*sW2)*(Int1(TAUt,LAMt)-Int2(TAUt,LAMt)) * (1.0-AlsMz/M_PI)
1406  + 2.0*(-0.5+2.0/3.0*sW2)*(Int1(TAUb,LAMb)-Int2(TAUb,LAMb))
1407  - 4.0*(0.5-4.0/3.0*sW2)*(Int1(TAUc,LAMc)-Int2(TAUc,LAMc))
1408  + 2.0*(-0.5+2.0*sW2)*(Int1(TAUtau,LAMtau)-Int2(TAUtau,LAMtau))
1409  + 2.0*(-0.5+2.0*sW2)*(Int1(TAUmu,LAMmu)-Int2(TAUmu,LAMmu)) )/sqrt(sW2*cW2)
1410  - sqrt(cW2/sW2)*(4.0*(3.0-sW2/cW2)*Int2(TAUw,LAMw)
1411  +((1.0+2.0/TAUw)*sW2/cW2-(5.0+2.0/TAUw))*Int1(TAUw,LAMw))).abs2();
1412 
1413  return (sqrt(cZgaEff));
1414 }
1415 
1417 {
1418  if (loopComputed) return cga_loop;
1419 
1420  double Mt=trueSM.getQuarks(QCD::TOP).getMass();
1421  double Mb=trueSM.getQuarks(QCD::BOTTOM).getMass();
1422  double Mc=trueSM.getQuarks(QCD::CHARM).getMass();
1425  double MW=trueSM.Mw_tree();
1426  double TAUt=4.0*Mt*Mt/(mHl*mHl);
1427  double TAUb=4.0*Mb*Mb/(mHl*mHl);
1428  double TAUc=4.0*Mc*Mc/(mHl*mHl);
1429  double TAUtau=4.0*Mtau*Mtau/(mHl*mHl);
1430  double TAUmu=4.0*Mmu*Mmu/(mHl*mHl);
1431  double TAUw=4.0*MW*MW/(mHl*mHl);
1432  double cgagaEff;
1433 
1434  cgagaEff = ( ct * (8./3.)*TAUt*(1.+(1.-TAUt)*f_func(TAUt)) * (1.0-AlsMz/M_PI)
1435  +cb * (2./3.)*TAUb*(1.+(1.-TAUb)*f_func(TAUb))
1436  +cc * (8./3.)*TAUc*(1.+(1.-TAUc)*f_func(TAUc))
1437  +ctau * 2.0*TAUtau*(1.+(1.-TAUtau)*f_func(TAUtau))
1438  +cmu * 2.0*TAUmu*(1.+(1.-TAUmu)*f_func(TAUmu))
1439  -cv * (2.0+3.0*TAUw+3.0*TAUw*(2.0-TAUw)*f_func(TAUw))
1440  +cga * 2.0 ).abs2();
1441 
1442  cgagaEff = cgagaEff / ( (8./3.)*TAUt*(1.+(1.-TAUt)*f_func(TAUt)) * (1.0-AlsMz/M_PI)
1443  + (2./3.)*TAUb*(1.+(1.-TAUb)*f_func(TAUb))
1444  + (8./3.)*TAUc*(1.+(1.-TAUc)*f_func(TAUc))
1445  + 2.0*TAUtau*(1.+(1.-TAUtau)*f_func(TAUtau))
1446  + 2.0*TAUmu*(1.+(1.-TAUmu)*f_func(TAUmu))
1447  - (2.0+3.0*TAUw+3.0*TAUw*(2.0-TAUw)*f_func(TAUw)) ).abs2();
1448 
1449  return (sqrt(cgagaEff));
1450 }
1451 
1453 {
1454  return cmu;
1455 }
1456 
1458 {
1459  return ctau;
1460 }
1461 
1463 {
1464  return cc;
1465 }
1466 
1468 {
1469  return ct;
1470 }
1471 
1473 {
1474  return cb;
1475 }
1476 
1477 
1479 
1481 {
1482  // Not needed/implemented yet
1483  return 0.0;
1484 }
1485 
1487  return (computecg() - 1.0);
1488 }
1489 
1491  return 0.0;
1492 }
1493 
1495  return 0.0;
1496 }
1497 
1499 
1500  double gSM = 2.0 * (trueSM.Mw_tree())* (trueSM.Mw_tree()) / (trueSM.v());
1501 
1502  return ( gSM * (computecV() - 1.0 ) );
1503 }
1504 
1506  return 0.0;
1507 }
1508 
1510  return 0.0;
1511 }
1512 
1514 
1515  double gSM = (trueSM.getMz()) * (trueSM.getMz()) / (trueSM.v());
1516 
1517  return ( gSM * (computecV() - 1.0 ) );
1518 }
1519 
1521  // Not needed/implemented yet
1522  return 0.0;
1523 }
1524 
1526  return (computecZga() - 1.0);
1527 }
1528 
1530  return 0.0;
1531 }
1532 
1534  // Not needed/implemented yet
1535  return 0.0;
1536 }
1537 
1539  return (computecgaga() - 1.0);
1540 }
1541 
1543 
1544  double gSM = -(p.getMass()) / (trueSM.v());
1545 
1546  if ( p.is("ELECTRON") ) {
1547  return ( 0.0 );
1548  } else if ( p.is("MU") ) {
1549  return ( gSM * (computecmu() - 1.0 ) );
1550  } else if ( p.is("TAU") ) {
1551  return ( gSM * (computectau() - 1.0 ) );
1552  } else if ( p.is("UP") ) {
1553  return ( 0.0 );
1554  } else if ( p.is("CHARM") ) {
1555  return ( gSM * (computecc() - 1.0 ) );
1556  } else if ( p.is("TOP") ) {
1557  return ( gSM * (computect() - 1.0 ) );
1558  } else if ( p.is("DOWN") ) {
1559  return ( 0.0 );
1560  } else if ( p.is("STRANGE") ) {
1561  return ( 0.0 );
1562  } else if ( p.is("BOTTOM") ) {
1563  return ( gSM * (computecb() - 1.0 ) );
1564  } else {
1565  return 0.0;
1566  }
1567 }
1568 
1570 
1572 {
1573  return computecmu();
1574 }
1575 
1577 {
1578  return computectau();
1579 }
1580 
1582 {
1583  return computecc();
1584 }
1585 
1587 {
1588  return computecb();
1589 }
1590 
1592 {
1593  return computecg();
1594 }
1595 
1597 {
1598  return computecV();
1599 }
1600 
1602 {
1603  return computecV();
1604 }
1605 
1607 {
1608  return computecgaga();
1609 }
1610 
1612 {
1613  return computecZga();
1614 }
1615 
1617 
1618 gslpp::complex HiggsChiral::f_func(const double x) const{
1619  if(x<1) {
1620  gslpp::complex z = -gslpp::complex::i()*M_PI;
1621  return -pow(log((1+sqrt(1-x))/(1-sqrt(1-x)))+z,2)/4.0;
1622  }
1623  else {
1624  return pow(asin(sqrt(1.0/x)),2);
1625  }
1626 }
1627 
1628 gslpp::complex HiggsChiral::g_func(const double x) const{
1629  if(x<1) {
1630  gslpp::complex z = -gslpp::complex::i()*M_PI;
1631  gslpp::complex gs1 = sqrt(1.0-x)*(log((1.0+sqrt(1.0-x))/(1.0-sqrt(1.0-x)))+z)/2.0;
1632  return gs1;
1633  }
1634  else {
1635  gslpp::complex gg1 = sqrt(x-1.0)*asin(sqrt(1.0/x));
1636  return gg1;
1637  }
1638 }
1639 
1640 gslpp::complex HiggsChiral::Int1(const double tau, const double lambda) const{
1641  return tau*lambda/(tau-lambda)/2.0+tau*tau*lambda*lambda/((tau-lambda)
1642  *(tau-lambda))/2.0*(f_func(tau)-f_func(lambda))+tau*tau*lambda/((tau-lambda)
1643  *(tau-lambda))*(g_func(tau)-g_func(lambda));
1644 }
1645 
1646 gslpp::complex HiggsChiral::Int2(const double tau, const double lambda) const{
1647  return -tau*lambda/(tau-lambda)/2.0*(f_func(tau)-f_func(lambda));
1648 }
QCD::TAU
Definition: QCD.h:316
HiggsChiral::muTHUttHbb
virtual double muTHUttHbb(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into in the current model ...
Definition: HiggsChiral.cpp:1272
sigmattH
Definition: NPSMEFT6dtopquark.h:606
HiggsChiral::muVBFHWW2l2v
virtual double muVBFHWW2l2v(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into in the current ...
Definition: HiggsChiral.cpp:869
HiggsChiral::muttHtautau
virtual double muttHtautau(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into in the current model ...
Definition: HiggsChiral.cpp:949
HiggsChiral::cbminuscc
virtual double cbminuscc() const
The value of .
Definition: HiggsChiral.cpp:422
HiggsChiral::deltaG1_hWW
virtual double deltaG1_hWW() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1490
HiggsChiral::BrHZZRatio
virtual double BrHZZRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:557
HiggsChiral::muggHpttH
virtual double muggHpttH(const double sqrt_s) const
The ratio between the sum of gluon-gluon fusion and t-tbar-Higgs associated production cross-section...
Definition: HiggsChiral.cpp:280
HiggsChiral::muggHbb
virtual double muggHbb(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:954
HiggsChiral::muVHWW2l2v
virtual double muVHWW2l2v(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:884
HiggsChiral::mueeHvvPol
virtual double mueeHvvPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const
The ratio between the associated production cross-section in the current model and in the Standard ...
Definition: HiggsChiral.cpp:198
HiggsChiral::HiggsChiral
HiggsChiral()
The default constructor.
Definition: HiggsChiral.cpp:14
HiggsChiral::getobsZgaLimitCMS
double getobsZgaLimitCMS() const
a getter for the value of the observed upper limit in from CMS
Definition: HiggsChiral.h:246
HiggsChiral::UpperLimitZgammaC13
virtual double UpperLimitZgammaC13(const double sqrt_s) const
Observable implementing the contribution to the likelihood from the upper limit in from CMS at 13 Te...
Definition: HiggsChiral.cpp:325
HiggsChiral::muTHUZHtautau
virtual double muTHUZHtautau(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1227
HiggsChiral::muTHUWHmumu
virtual double muTHUWHmumu(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1202
HiggsChiral::cZga_loop
double cZga_loop
Definition: HiggsChiral.h:1923
HiggsChiral::kappaWeff
virtual double kappaWeff() const
The effective coupling .
Definition: HiggsChiral.cpp:1601
HiggsChiral::deltaG_hff
virtual gslpp::complex deltaG_hff(const Particle p) const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1542
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:989
HiggsChiral::muggHmumu
virtual double muggHmumu(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:894
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
StandardModel::v
virtual double v() const
The Higgs vacuum expectation value.
Definition: StandardModel.cpp:917
HiggsChiral::deltaG2_hZA
virtual double deltaG2_hZA() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1529
HiggsChiral::muTHUWHtautau
virtual double muTHUWHtautau(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1232
HiggsChiral::cZga
double cZga
Definition: HiggsChiral.h:1911
HiggsChiral::GammaZZ
virtual double GammaZZ() const
The decay width in the current model.
Definition: HiggsChiral.cpp:458
HiggsChiral::muTHUttHmumu
virtual double muTHUttHmumu(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into in the current model ...
Definition: HiggsChiral.cpp:1212
HiggsChiral::expZgaLimitATLAS
double expZgaLimitATLAS
Definition: HiggsChiral.h:1918
QCD::BOTTOM
Definition: QCD.h:329
HiggsChiral::kappaZeff
virtual double kappaZeff() const
The effective coupling .
Definition: HiggsChiral.cpp:1596
HiggsChiral::muTHUVBFHWW
virtual double muTHUVBFHWW(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into in the current ...
Definition: HiggsChiral.cpp:1132
Particle::is
bool is(std::string name_i) const
Definition: Particle.cpp:23
Particle
A class for particles.
Definition: Particle.h:26
StandardModel::computeSigmaWH
double computeSigmaWH(const double sqrt_s) const
The WH production cross section in the Standard Model.
Definition: StandardModel.h:2088
HiggsChiral::muTHUVHWW
virtual double muTHUVHWW(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1147
HiggsChiral::getexpZgaLimitCMS
double getexpZgaLimitCMS() const
a getter for the experimental value of the expected upper limit in from CMS
Definition: HiggsChiral.h:282
HiggsChiral::deltaG_hAARatio
virtual double deltaG_hAARatio() const
The full new physics contribution to the coupling of the effective interaction , including new local ...
Definition: HiggsChiral.cpp:1538
HiggsChiral::muTHUVHZZ4l
virtual double muTHUVHZZ4l(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1117
HiggsChiral::muWHZga
virtual double muWHZga(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:759
HiggsChiral::muWHbb
virtual double muWHbb(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:969
HiggsChiral::deltaG3_hWW
virtual double deltaG3_hWW() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1498
HiggsChiral::kappaceff
virtual double kappaceff() const
The effective coupling .
Definition: HiggsChiral.cpp:1581
HiggsChiral::muWHZZ
virtual double muWHZZ(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:789
HiggsChiral::muVBFgamma
virtual double muVBFgamma(const double sqrt_s) const
The ratio between the vector-boson fusion Higgs production cross-section in association with a hard ...
Definition: HiggsChiral.cpp:178
HiggsChiral::muTHUVBFHZZ4l
virtual double muTHUVBFHZZ4l(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into in the current ...
Definition: HiggsChiral.cpp:1102
StandardModel::computeBrHtotautau
double computeBrHtotautau() const
The Br in the Standard Model.
Definition: StandardModel.h:2264
HiggsChiral::ct
double ct
Definition: HiggsChiral.h:1904
HiggsChiral::muggHZga
virtual double muggHZga(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:744
HiggsChiral::BrHggRatio
virtual double BrHggRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:536
HiggsChiral::muTHUVHmumu
virtual double muTHUVHmumu(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1207
HiggsChiral::muTHUZHgaga
virtual double muTHUZHgaga(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into 2 photons in the curren...
Definition: HiggsChiral.cpp:1017
HiggsChiral::deltaG1_hZZ
virtual double deltaG1_hZZ() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1505
HiggsChiral::BrHmumuRatio
virtual double BrHmumuRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:627
HiggsChiral::computecmu
virtual double computecmu() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1452
HiggsChiral::muVHtautau
virtual double muVHtautau(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:944
HiggsChiral::CheckParameters
virtual bool CheckParameters(const std::map< std::string, double > &DPars)
A method to check if all the mandatory parameters for HiggsChiral have been provided in model initial...
Definition: HiggsChiral.cpp:110
HiggsChiral::muVHWW
virtual double muVHWW(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:854
HiggsChiral::obsZgaLimitCMS13
double obsZgaLimitCMS13
Definition: HiggsChiral.h:1913
HiggsChiral::computectau
virtual double computectau() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1457
HiggsChiral::FlagUniversalcf
bool FlagUniversalcf
A boolean flag that is true if all cf take the same universal value.
Definition: HiggsChiral.h:1926
HiggsChiral::muTHUttHWW2l2v
virtual double muTHUttHWW2l2v(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into in the current model ...
Definition: HiggsChiral.cpp:1182
HiggsChiral::muVH
virtual double muVH(const double sqrt_s) const
The ratio between the WH+ZH associated production cross-section in the current model and in the Stan...
Definition: HiggsChiral.cpp:253
HiggsChiral::muTHUVHWW2l2v
virtual double muTHUVHWW2l2v(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1177
HiggsChiral::muZHZga
virtual double muZHZga(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:754
HiggsChiral::ctau
double ctau
Definition: HiggsChiral.h:1907
HiggsChiral::muTHUZHWW
virtual double muTHUZHWW(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1137
HiggsChiral::muVBFHZZ4l
virtual double muVBFHZZ4l(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into in the current ...
Definition: HiggsChiral.cpp:809
StandardModel::computeSigmaWF
double computeSigmaWF(const double sqrt_s) const
The W fusion contribution to higgs-production cross section in the Standard Model.
Definition: StandardModel.h:2017
HiggsChiral::muWH
virtual double muWH(const double sqrt_s) const
The ratio between the W-Higgs associated production cross-section in the current model and in the St...
Definition: HiggsChiral.cpp:213
HiggsChiral::muZHWW
virtual double muZHWW(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:844
HiggsChiral::muTHUZHZga
virtual double muTHUZHZga(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1047
HiggsChiral::muTHUZHZZ4l
virtual double muTHUZHZZ4l(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1107
StandardModel::GF
double GF
The Fermi constant in .
Definition: StandardModel.h:2511
HiggsChiral::deltaG2_hZZ
virtual double deltaG2_hZZ() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1509
HiggsChiral::deltaG_hgg
virtual double deltaG_hgg() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1480
HiggsChiral::muTHUggHtautau
virtual double muTHUggHtautau(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:1217
HiggsChiral::muTHUttHZga
virtual double muTHUttHZga(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into in the current model ...
Definition: HiggsChiral.cpp:1062
Model::addMissingModelParameter
void addMissingModelParameter(const std::string &missingParameterName)
Definition: Model.h:232
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
HiggsChiral::muZHmumu
virtual double muZHmumu(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:904
HiggsChiral::muZHtautau
virtual double muZHtautau(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:934
QCD::CHARM
Definition: QCD.h:326
HiggsChiral::muTHUWHbb
virtual double muTHUWHbb(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1262
HiggsChiral::muWHZZ4l
virtual double muWHZZ4l(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:819
HiggsChiral::muTHUggHZZ4mu
virtual double muTHUggHZZ4mu(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:1297
HiggsChiral::muVBFHZZ
virtual double muVBFHZZ(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into in the current ...
Definition: HiggsChiral.cpp:779
StandardModel::computeBrHtobb
double computeBrHtobb() const
The Br in the Standard Model.
Definition: StandardModel.h:2299
HiggsChiral::mueeWBF
virtual double mueeWBF(const double sqrt_s) const
The ratio between the production cross-section in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:183
gslpp::complex
A class for defining operations on and functions of complex numbers.
Definition: gslpp_complex.h:35
HiggsChiral::Gammacc
virtual double Gammacc() const
The decay width in the current model.
Definition: HiggsChiral.cpp:520
HiggsChiral::muTHUttHtautau
virtual double muTHUttHtautau(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into in the current model ...
Definition: HiggsChiral.cpp:1242
StandardModel::mHl
double mHl
The Higgs mass in GeV.
Definition: StandardModel.h:2514
gslpp::log
complex log(const complex &z)
Definition: gslpp_complex.cpp:342
HiggsChiral::muTHUggHgaga
virtual double muTHUggHgaga(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into 2...
Definition: HiggsChiral.cpp:1007
HiggsChiral::muTHUggHmumu
virtual double muTHUggHmumu(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:1187
HiggsChiral::obliqueU
virtual double obliqueU() const
The oblique parameter .
Definition: HiggsChiral.cpp:159
gslpp::complex::abs2
double abs2() const
Definition: gslpp_complex.cpp:86
HiggsChiral::muTHUZHWW2l2v
virtual double muTHUZHWW2l2v(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1167
HiggsChiral::cb
double cb
Definition: HiggsChiral.h:1905
HiggsChiral::GammaZga
virtual double GammaZga() const
The decay width in the current model.
Definition: HiggsChiral.cpp:463
HiggsChiral::deltaG1_hZARatio
virtual double deltaG1_hZARatio() const
The full new physics contribution to the coupling of the effective interaction , including new local ...
Definition: HiggsChiral.cpp:1525
HiggsChiral::computect
virtual double computect() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1467
HiggsChiral::obliqueT
virtual double obliqueT() const
The oblique parameter .
Definition: HiggsChiral.cpp:147
HiggsChiral::getobsZgaLimitATLAS13
double getobsZgaLimitATLAS13() const
a getter for the value of the observed upper limit in from ATLAS at 13 TeV
Definition: HiggsChiral.h:219
HiggsChiral::muWHgaga
virtual double muWHgaga(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into 2 photons in the curren...
Definition: HiggsChiral.cpp:729
HiggsChiral::muTHUVHZZ
virtual double muTHUVHZZ(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1087
HiggsChiral::mueettHPol
virtual double mueettHPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const
The ratio between the production cross-section in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:296
HiggsChiral::expZgaLimitCMS
double expZgaLimitCMS
Definition: HiggsChiral.h:1919
HiggsChiral::muTHUggHZZ4l
virtual double muTHUggHZZ4l(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:1097
HiggsChiral::muTHUggHZZ
virtual double muTHUggHZZ(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:1067
HiggsChiral::muTHUWHZZ4l
virtual double muTHUWHZZ4l(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1112
StandardModel::cW2
virtual double cW2() const
Definition: StandardModel.cpp:994
HiggsChiral::muTHUZHmumu
virtual double muTHUZHmumu(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1197
HiggsChiral::muttHbb
virtual double muttHbb(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into in the current model ...
Definition: HiggsChiral.cpp:979
StandardModel::ale
double ale
The fine-structure constant .
Definition: StandardModel.h:2512
HiggsChiral::BrHZZ2e2muRatio
virtual double BrHZZ2e2muRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:578
HiggsChiral::BrHZgamumuRatio
virtual double BrHZgamumuRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:613
HiggsChiral::muTHUVHBRinv
virtual double muTHUVHBRinv(const double sqrt_s) const
The ratio between the VH production cross-section in the current model and in the Standard Model,...
Definition: HiggsChiral.cpp:1287
HiggsChiral::UpperLimitZgammaC
virtual double UpperLimitZgammaC(const double sqrt_s) const
Observable implementing the contribution to the likelihood from the upper limit in from CMS.
Definition: HiggsChiral.cpp:373
HiggsChiral::muTHUggHWW2l2v
virtual double muTHUggHWW2l2v(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:1157
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
HiggsChiral::deltaG_hAA
virtual double deltaG_hAA() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1533
HiggsChiral::muttHZZ4l
virtual double muttHZZ4l(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into in the current model ...
Definition: HiggsChiral.cpp:829
HiggsChiral::muTHUggHbb
virtual double muTHUggHbb(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:1247
HiggsChiral::obsZgaLimitATLAS
double obsZgaLimitATLAS
Definition: HiggsChiral.h:1914
HiggsChiral::PostUpdate
virtual bool PostUpdate()
The post-update method for HiggsChiral.
Definition: HiggsChiral.cpp:43
HiggsChiral::muTHUggHWW
virtual double muTHUggHWW(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:1127
HiggsChiral::muVBF
virtual double muVBF(const double sqrt_s) const
The ratio between the vector-boson fusion Higgs production cross-section in the current model and in...
Definition: HiggsChiral.cpp:173
HiggsChiral::muTHUVBFHtautau
virtual double muTHUVBFHtautau(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into in the current ...
Definition: HiggsChiral.cpp:1222
HiggsChiral::expZgaLimitCMS13
double expZgaLimitCMS13
Definition: HiggsChiral.h:1917
HiggsChiral::NHChiralvars
static const int NHChiralvars
The number of the model parameters.
Definition: HiggsChiral.h:100
HiggsChiral::mueeZHPol
virtual double mueeZHPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const
The ratio between the associated production cross-section in the current model and in the Standard ...
Definition: HiggsChiral.cpp:238
HiggsChiral::muttHZga
virtual double muttHZga(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into in the current model ...
Definition: HiggsChiral.cpp:769
HiggsChiral::kappaAeff
virtual double kappaAeff() const
The effective coupling .
Definition: HiggsChiral.cpp:1606
NPbase::Br_H_inv
virtual double Br_H_inv() const
The branching ratio of the of the Higgs into invisible particles.
Definition: NPbase.h:1474
HiggsChiral::BrHZgallRatio
virtual double BrHZgallRatio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:599
StandardModel::computeSigmaZF
double computeSigmaZF(const double sqrt_s) const
The Z fusion contribution to higgs-production cross section in the Standard Model.
Definition: StandardModel.h:2045
HiggsChiral::muTHUttHZZ4l
virtual double muTHUttHZZ4l(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into in the current model ...
Definition: HiggsChiral.cpp:1122
HiggsChiral::muVBFHgaga
virtual double muVBFHgaga(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into 2 photons in the...
Definition: HiggsChiral.cpp:719
HiggsChiral::getexpZgaLimitATLAS
double getexpZgaLimitATLAS() const
a getter for the experimental value of the expected upper limit in from ATLAS
Definition: HiggsChiral.h:273
HiggsChiral::muggHWW
virtual double muggHWW(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:834
NPbase::trueSM
StandardModel trueSM
Definition: NPbase.h:2787
StandardModel::computeBrHtoZZ
double computeBrHtoZZ() const
The Br in the Standard Model.
Definition: StandardModel.h:2208
StandardModel::c02
double c02() const
The square of the cosine of the weak mixing angle defined without weak radiative corrections.
Definition: StandardModel.cpp:939
HiggsChiral::muZHZZ4l
virtual double muZHZZ4l(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:814
HiggsChiral::muTHUVBFHZga
virtual double muTHUVBFHZga(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into in the current ...
Definition: HiggsChiral.cpp:1042
HiggsChiral::muTHUVBFHZZ
virtual double muTHUVBFHZZ(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into in the current ...
Definition: HiggsChiral.cpp:1072
HiggsChiral::UpperLimitZgammaA
virtual double UpperLimitZgammaA(const double sqrt_s) const
Observable implementing the contribution to the likelihood from the upper limit in from ATLAS.
Definition: HiggsChiral.cpp:349
HiggsChiral::deltaG3_hZZ
virtual double deltaG3_hZZ() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1513
HiggsChiral::deltaG_hggRatio
virtual double deltaG_hggRatio() const
The full new physics contribution to the coupling of the effective interaction , including new local ...
Definition: HiggsChiral.cpp:1486
HiggsChiral::BrHZZ4eRatio
virtual double BrHZZ4eRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:571
HiggsChiral::Int1
gslpp::complex Int1(const double tau, const double lambda) const
Definition: HiggsChiral.cpp:1640
HiggsChiral::BrHZZ4lRatio
virtual double BrHZZ4lRatio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:564
HiggsChiral::f_func
gslpp::complex f_func(const double x) const
Definition: HiggsChiral.cpp:1618
Particle::getMass
const double & getMass() const
A get method to access the particle mass.
Definition: Particle.h:61
HiggsChiral::cga
double cga
Definition: HiggsChiral.h:1910
HiggsChiral::muggHgagaInt
virtual double muggHgagaInt(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into 2...
Definition: HiggsChiral.cpp:667
HiggsChiral::muTHUVHtautau
virtual double muTHUVHtautau(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1237
HiggsChiral::cga_loop
double cga_loop
Definition: HiggsChiral.h:1922
HiggsChiral::Gammabb
virtual double Gammabb() const
The decay width in the current model.
Definition: HiggsChiral.cpp:525
StandardModel::AlsMz
double AlsMz
The strong coupling constant at the Z-boson mass, .
Definition: StandardModel.h:2509
NPbase
The auxiliary base model class for other model classes.
Definition: NPbase.h:66
HiggsChiral::muTHUVHinv
virtual double muTHUVHinv(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into invisible states in the...
Definition: HiggsChiral.cpp:1292
StandardModel::sW2
double sW2() const
Definition: StandardModel.cpp:1005
HiggsChiral::muppHmumu
virtual double muppHmumu(const double sqrt_s) const
The ratio between the Higgs production cross-section with subsequent decay into in the current mode...
Definition: HiggsChiral.cpp:984
HiggsChiral::Gammagg
virtual double Gammagg() const
The decay width in the current model.
Definition: HiggsChiral.cpp:437
HiggsChiral::obsZgaLimitATLAS13
double obsZgaLimitATLAS13
Definition: HiggsChiral.h:1912
HiggsChiral::UpperLimitZgammaA13
virtual double UpperLimitZgammaA13(const double sqrt_s) const
Observable implementing the contribution to the likelihood from the upper limit in from ATLAS at 13 ...
Definition: HiggsChiral.cpp:301
HiggsChiral::muWHWW2l2v
virtual double muWHWW2l2v(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:879
QCD::TOP
Definition: QCD.h:328
HiggsChiral::muggHWW2l2v
virtual double muggHWW2l2v(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:864
HiggsChiral::muTHUttHWW
virtual double muTHUttHWW(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into in the current model ...
Definition: HiggsChiral.cpp:1152
HiggsChiral::getobsZgaLimitATLAS
double getobsZgaLimitATLAS() const
a getter for the value of the observed upper limit in from ATLAS
Definition: HiggsChiral.h:237
HiggsChiral::GammaTotal
virtual double GammaTotal() const
The total decay width of the Higgs boson in the current model.
Definition: HiggsChiral.cpp:530
HiggsChiral::muVBFHZga
virtual double muVBFHZga(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into in the current ...
Definition: HiggsChiral.cpp:749
StandardModel::computeGammaHTotal
double computeGammaHTotal() const
The Higgs total width in the Standard Model.
Definition: StandardModel.h:2311
gslpp::pow
complex pow(const complex &z1, const complex &z2)
Definition: gslpp_complex.cpp:395
HiggsChiral::muVBFHtautau
virtual double muVBFHtautau(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into in the current ...
Definition: HiggsChiral.cpp:929
HiggsChiral::muTHUVBFHgaga
virtual double muTHUVBFHgaga(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into 2 photons in the...
Definition: HiggsChiral.cpp:1012
HiggsChiral::computeGammaTotalRatio
virtual double computeGammaTotalRatio() const
The ratio of the in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:1314
Model::raiseMissingModelParameterCount
void raiseMissingModelParameterCount()
Definition: Model.h:242
gslpp::sqrt
complex sqrt(const complex &z)
Definition: gslpp_complex.cpp:385
HiggsChiral::cbminusctau
virtual double cbminusctau() const
The value of .
Definition: HiggsChiral.cpp:427
gslpp::complex::i
static const complex & i()
Definition: gslpp_complex.cpp:154
HiggsChiral::mueeZllH
virtual double mueeZllH(const double sqrt_s) const
The ratio between the associated production cross-section in the current model and in the Standard ...
Definition: HiggsChiral.cpp:228
StandardModel::computeBrHtoZga
double computeBrHtoZga() const
The Br in the Standard Model.
Definition: StandardModel.h:2230
HiggsChiral::kappaGeff
virtual double kappaGeff() const
The effective coupling .
Definition: HiggsChiral.cpp:1591
StandardModel::lambda
double lambda
The CKM parameter in the Wolfenstein parameterization.
Definition: StandardModel.h:2524
HiggsChiral::muVBFHWW
virtual double muVBFHWW(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into in the current ...
Definition: HiggsChiral.cpp:839
HiggsChiral::cgminuscga
virtual double cgminuscga() const
The value of .
Definition: HiggsChiral.cpp:407
HiggsChiral::muTHUttHgaga
virtual double muTHUttHgaga(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into 2 photons in the curre...
Definition: HiggsChiral.cpp:1032
StandardModel::computeBrHtogaga
double computeBrHtogaga() const
The Br in the Standard Model.
Definition: StandardModel.h:2242
HiggsChiral::deltaG2_hWW
virtual double deltaG2_hWW() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1494
StandardModel::computeSigmaggH
double computeSigmaggH(const double sqrt_s) const
The ggH cross section in the Standard Model.
Definition: StandardModel.h:1883
HiggsChiral::cg_loop
double cg_loop
Definition: HiggsChiral.h:1921
HiggsChiral::mueeZllHPol
virtual double mueeZllHPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const
The ratio between the associated production cross-section in the current model and in the Standard ...
Definition: HiggsChiral.cpp:243
HiggsChiral::cVplusctau
virtual double cVplusctau() const
The value of .
Definition: HiggsChiral.cpp:417
HiggsChiral::kappamueff
virtual double kappamueff() const
The effective coupling .
Definition: HiggsChiral.cpp:1571
HiggsChiral::mueeZH
virtual double mueeZH(const double sqrt_s) const
The ratio between the associated production cross-section in the current model and in the Standard ...
Definition: HiggsChiral.cpp:223
HiggsChiral::muZHWW2l2v
virtual double muZHWW2l2v(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:874
HiggsChiral::getexpZgaLimitATLAS13
double getexpZgaLimitATLAS13() const
a getter for the experimental value of the expected upper limit in from ATLAS at 13 TeV
Definition: HiggsChiral.h:255
HiggsChiral::muVHZZ
virtual double muVHZZ(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:794
HiggsChiral::kappataueff
virtual double kappataueff() const
The effective coupling .
Definition: HiggsChiral.cpp:1576
HiggsChiral::getexpZgaLimitCMS13
double getexpZgaLimitCMS13() const
a getter for the experimental value of the expected upper limit in from CMS at 13 TeV
Definition: HiggsChiral.h:264
HiggsChiral::muTHUWHZZ
virtual double muTHUWHZZ(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1082
HiggsChiral::muTHUVHbb
virtual double muTHUVHbb(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1267
StandardModel::computeSigmattH
double computeSigmattH(const double sqrt_s) const
The ttH production cross section in the Standard Model.
Definition: StandardModel.h:2157
QCD::getQuarks
Particle getQuarks(const QCD::quark q) const
A get method to access a quark as an object of the type Particle.
Definition: QCD.h:534
HiggsChiral::cv
double cv
Definition: HiggsChiral.h:1903
HiggsChiral::Int2
gslpp::complex Int2(const double tau, const double lambda) const
Definition: HiggsChiral.cpp:1646
HiggsChiral::mueeZqqH
virtual double mueeZqqH(const double sqrt_s) const
The ratio between the associated production cross-section in the current model and in the Standard ...
Definition: HiggsChiral.cpp:233
HiggsChiral::muttHgaga
virtual double muttHgaga(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into 2 photons in the curre...
Definition: HiggsChiral.cpp:739
HiggsChiral::GammaWW
virtual double GammaWW() const
The decay width in the current model.
Definition: HiggsChiral.cpp:453
HiggsChiral::expZgaLimitATLAS13
double expZgaLimitATLAS13
Definition: HiggsChiral.h:1916
HiggsChiral::muVHmumu
virtual double muVHmumu(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:914
HiggsChiral::obliqueS
virtual double obliqueS() const
The oblique parameter .
Definition: HiggsChiral.cpp:136
HiggsChiral::muggHZZ
virtual double muggHZZ(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:774
HiggsChiral::BrHgagaRatio
virtual double BrHgagaRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:620
HiggsChiral::muTHUVBFBRinv
virtual double muTHUVBFBRinv(const double sqrt_s) const
The ratio between the VBF production cross-section in the current model and in the Standard Model,...
Definition: HiggsChiral.cpp:1277
HiggsChiral::loopComputed
bool loopComputed
Definition: HiggsChiral.h:1924
HiggsChiral::cVpluscb
virtual double cVpluscb() const
The value of .
Definition: HiggsChiral.cpp:412
HiggsChiral::mutHq
virtual double mutHq(const double sqrt_s) const
The ratio between the t-q-Higgs associated production cross-section in the current model and in the ...
Definition: HiggsChiral.cpp:268
HiggsChiral::setFlag
virtual bool setFlag(const std::string name, const bool value)
A method to set a flag of HiggsChiral.
Definition: HiggsChiral.cpp:122
HiggsChiral::muWHmumu
virtual double muWHmumu(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:909
HiggsChiral.h
HiggsChiral::muTHUVBFHWW2l2v
virtual double muTHUVBFHWW2l2v(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into in the current ...
Definition: HiggsChiral.cpp:1162
HiggsChiral::muZHbb
virtual double muZHbb(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:964
HiggsChiral::muggHtautau
virtual double muggHtautau(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:924
HiggsChiral::mueeZBF
virtual double mueeZBF(const double sqrt_s) const
The ratio between the production cross-section in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:203
HiggsChiral::computecc
virtual double computecc() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1462
HiggsChiral::BrHZgaRatio
virtual double BrHZgaRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:592
HiggsChiral::muttHWW2l2v
virtual double muttHWW2l2v(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into in the current model ...
Definition: HiggsChiral.cpp:889
HiggsChiral::muttHWW
virtual double muttHWW(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into in the current model ...
Definition: HiggsChiral.cpp:859
HiggsChiral::muZH
virtual double muZH(const double sqrt_s) const
The ratio between the Z-Higgs associated production cross-section in the current model and in the St...
Definition: HiggsChiral.cpp:218
StandardModel::computeSigmaZH
double computeSigmaZH(const double sqrt_s) const
The ZH production cross section in the Standard Model.
Definition: StandardModel.h:2121
StandardModel::getMz
double getMz() const
A get method to access the mass of the boson .
Definition: StandardModel.h:718
HiggsChiral::setParameter
virtual void setParameter(const std::string name, const double &value)
A method to set the value of a parameter of HiggsChiral.
Definition: HiggsChiral.cpp:70
HiggsChiral::kappabeff
virtual double kappabeff() const
The effective coupling .
Definition: HiggsChiral.cpp:1586
StandardModel::computeBrHtocc
double computeBrHtocc() const
The Br in the Standard Model.
Definition: StandardModel.h:2276
HiggsChiral::Gammagaga
virtual double Gammagaga() const
The decay width in the current model.
Definition: HiggsChiral.cpp:490
HiggsChiral::muTHUWHWW2l2v
virtual double muTHUWHWW2l2v(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1172
HiggsChiral::muTHUVBFHinv
virtual double muTHUVBFHinv(const double sqrt_s) const
The ratio between the VBF production cross-section with subsequent decay into invisible states in th...
Definition: HiggsChiral.cpp:1282
HiggsChiral::muTHUZHbb
virtual double muTHUZHbb(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1257
HiggsChiral::ccminusctau
virtual double ccminusctau() const
The value of .
Definition: HiggsChiral.cpp:432
HiggsChiral::BrHccRatio
virtual double BrHccRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:641
HiggsChiral::g_func
gslpp::complex g_func(const double x) const
Definition: HiggsChiral.cpp:1628
StandardModel::Mw_tree
virtual double Mw_tree() const
The tree-level mass of the boson, .
Definition: StandardModel.cpp:925
HiggsChiral::muWHWW
virtual double muWHWW(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:849
HiggsChiral::kappaZAeff
virtual double kappaZAeff() const
The effective coupling .
Definition: HiggsChiral.cpp:1611
StandardModel::computeBrHtomumu
double computeBrHtomumu() const
The Br in the Standard Model.
Definition: StandardModel.h:2253
HiggsChiral::muTHUWHWW
virtual double muTHUWHWW(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1142
HiggsChiral::computecV
virtual double computecV() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1365
StandardModel::computeBrHtoWW
double computeBrHtoWW() const
The Br in the Standard Model.
Definition: StandardModel.h:2196
HiggsChiral::muZHZZ
virtual double muZHZZ(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:784
HiggsChiral::muggH
virtual double muggH(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section in the current model and in ...
Definition: HiggsChiral.cpp:166
HiggsChiral::HChiralvars
static const std::string HChiralvars[NHChiralvars]
A string array containing the labels of the model parameters in HiggsKvKf.
Definition: HiggsChiral.h:105
HiggsChiral::muVHgaga
virtual double muVHgaga(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into 2 photons in the curren...
Definition: HiggsChiral.cpp:734
HiggsChiral::muttHmumu
virtual double muttHmumu(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into in the current model ...
Definition: HiggsChiral.cpp:919
HiggsChiral::mueeWBFPol
virtual double mueeWBFPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const
The ratio between the production cross-section in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:188
HiggsChiral::muTHUggHZgamumu
virtual double muTHUggHZgamumu(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:1302
HiggsChiral::BrHtoinvRatio
virtual double BrHtoinvRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:655
HiggsChiral::Gammamumu
virtual double Gammamumu() const
The decay width in the current model.
Definition: HiggsChiral.cpp:510
HiggsChiral::muTHUVBFHmumu
virtual double muTHUVBFHmumu(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into in the current ...
Definition: HiggsChiral.cpp:1192
NPbase::PostUpdate
virtual bool PostUpdate()
The postupdate method for NPbase.
Definition: NPbase.cpp:23
HiggsChiral::mueeHvv
virtual double mueeHvv(const double sqrt_s) const
The ratio between the associated production cross-section in the current model and in the Standard ...
Definition: HiggsChiral.cpp:193
StandardModel::computeBrHtogg
double computeBrHtogg() const
The Br in the Standard Model.
Definition: StandardModel.h:2185
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
HiggsChiral::muTHUZHZZ
virtual double muTHUZHZZ(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1077
StandardModel::computeSigmaZWF
double computeSigmaZWF(const double sqrt_s) const
The Z W interference fusion contribution to higgs-production cross section in the Standard Model.
Definition: StandardModel.h:2072
HiggsChiral::computecg
virtual double computecg() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1341
HiggsChiral::muWHtautau
virtual double muWHtautau(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:939
HiggsChiral::muVHZZ4l
virtual double muVHZZ4l(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:824
HiggsChiral::getobsZgaLimitCMS13
double getobsZgaLimitCMS13() const
a getter for the value of the observed upper limit in from CMS at 13 TeV
Definition: HiggsChiral.h:228
HiggsChiral::deltaG1_hZA
virtual double deltaG1_hZA() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1520
HiggsChiral::cgplusct
virtual double cgplusct() const
The value of .
Definition: HiggsChiral.cpp:397
HiggsChiral::muggHZZ4l
virtual double muggHZZ4l(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:804
HiggsChiral::muttHZZ
virtual double muttHZZ(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into in the current model ...
Definition: HiggsChiral.cpp:799
HiggsChiral::muggHgaga
virtual double muggHgaga(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into 2...
Definition: HiggsChiral.cpp:662
HiggsChiral::cgaplusct
virtual double cgaplusct() const
The value of .
Definition: HiggsChiral.cpp:402
HiggsChiral::obsZgaLimitCMS
double obsZgaLimitCMS
Definition: HiggsChiral.h:1915
HiggsChiral::muTHUggHZga
virtual double muTHUggHZga(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section with subsequent decay into ...
Definition: HiggsChiral.cpp:1037
HiggsChiral::mueeZBFPol
virtual double mueeZBFPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const
The ratio between the production cross-section in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:208
HiggsChiral::computecgaga
virtual double computecgaga() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1416
HiggsChiral::BrHWWRatio
virtual double BrHWWRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:543
HiggsChiral::cc
double cc
Definition: HiggsChiral.h:1906
HiggsChiral::mueettH
virtual double mueettH(const double sqrt_s) const
The ratio between the production cross-section in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:291
HiggsChiral::muVHbb
virtual double muVHbb(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:974
HiggsChiral::muTHUVBFHbb
virtual double muTHUVBFHbb(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into in the current ...
Definition: HiggsChiral.cpp:1252
HiggsChiral::muVBFHmumu
virtual double muVBFHmumu(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into in the current ...
Definition: HiggsChiral.cpp:899
HiggsChiral::computecb
virtual double computecb() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1472
HiggsChiral::cg
double cg
Definition: HiggsChiral.h:1909
HiggsChiral::muVBFHbb
virtual double muVBFHbb(const double sqrt_s) const
The ratio between the VBF Higgs production cross-section with subsequent decay into in the current ...
Definition: HiggsChiral.cpp:959
HiggsChiral::muTHUWHZga
virtual double muTHUWHZga(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1052
HiggsChiral::BrHtautauRatio
virtual double BrHtautauRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:634
HiggsChiral::BrHZgaeeRatio
virtual double BrHZgaeeRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:606
HiggsChiral::muttH
virtual double muttH(const double sqrt_s) const
The ratio between the t-tbar-Higgs associated production cross-section in the current model and in t...
Definition: HiggsChiral.cpp:263
HiggsChiral::muTHUttHZZ
virtual double muTHUttHZZ(const double sqrt_s) const
The ratio between the ttH production cross-section with subsequent decay into in the current model ...
Definition: HiggsChiral.cpp:1092
HiggsChiral::muTHUWHgaga
virtual double muTHUWHgaga(const double sqrt_s) const
The ratio between the WH production cross-section with subsequent decay into 2 photons in the curren...
Definition: HiggsChiral.cpp:1022
HiggsChiral::muVBFpVH
virtual double muVBFpVH(const double sqrt_s) const
The ratio between the sum of VBF and WH+ZH associated production cross-section in the current model ...
Definition: HiggsChiral.cpp:258
HiggsChiral::muZHgaga
virtual double muZHgaga(const double sqrt_s) const
The ratio between the ZH production cross-section with subsequent decay into 2 photons in the curren...
Definition: HiggsChiral.cpp:724
HiggsChiral::BrHWW2l2vRatio
virtual double BrHWW2l2vRatio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:550
HiggsChiral::cmu
double cmu
Definition: HiggsChiral.h:1908
HiggsChiral::BrHbbRatio
virtual double BrHbbRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:648
QCD::MU
Definition: QCD.h:314
HiggsChiral::computecZga
virtual double computecZga() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1370
HiggsChiral::BrHZZ4muRatio
virtual double BrHZZ4muRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:585
HiggsChiral::muTHUVHgaga
virtual double muTHUVHgaga(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into 2 photons in the curren...
Definition: HiggsChiral.cpp:1027
HiggsChiral::muVHZga
virtual double muVHZga(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:764
HiggsChiral::Gammatautau
virtual double Gammatautau() const
The decay width in the current model.
Definition: HiggsChiral.cpp:515
HiggsChiral::mueeZqqHPol
virtual double mueeZqqHPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const
The ratio between the associated production cross-section in the current model and in the Standard ...
Definition: HiggsChiral.cpp:248
StandardModel::getLeptons
Particle getLeptons(const QCD::lepton p) const
A get method to retrieve the member object of a lepton.
Definition: StandardModel.h:709
HiggsChiral::muTHUVHZga
virtual double muTHUVHZga(const double sqrt_s) const
The ratio between the VH production cross-section with subsequent decay into in the current model a...
Definition: HiggsChiral.cpp:1057