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
10const std::string HiggsChiral::HChiralvars[NHChiralvars] = {
11 "cv", "ct", "cb", "cc", "cs", "ctau", "cmu", "cg", "cga", "cZga", "obsZgaLimitATLAS13", "obsZgaLimitCMS13", "obsZgaLimitATLAS", "obsZgaLimitCMS", "expZgaLimitATLAS13", "expZgaLimitCMS13", "expZgaLimitATLAS", "expZgaLimitCMS"
12};
13
14HiggsChiral::HiggsChiral()
15: NPbase()
16{
17 FlagUniversalcf = false;
18 FlagUniversalcvcf = false;
19
20 ModelParamMap.insert(std::make_pair("cv", std::cref(cv)));
21 ModelParamMap.insert(std::make_pair("ct", std::cref(ct)));
22 ModelParamMap.insert(std::make_pair("cb", std::cref(cb)));
23 ModelParamMap.insert(std::make_pair("cc", std::cref(cc)));
24 ModelParamMap.insert(std::make_pair("cs", std::cref(cs)));
25 ModelParamMap.insert(std::make_pair("ctau", std::cref(ctau)));
26 ModelParamMap.insert(std::make_pair("cmu", std::cref(cmu)));
27 ModelParamMap.insert(std::make_pair("cg", std::cref(cg)));
28 ModelParamMap.insert(std::make_pair("cga", std::cref(cga)));
29 ModelParamMap.insert(std::make_pair("cZga", std::cref(cZga)));
30 ModelParamMap.insert(std::make_pair("obsZgaLimitATLAS13", std::cref(obsZgaLimitATLAS13)));
31 ModelParamMap.insert(std::make_pair("obsZgaLimitCMS13", std::cref(obsZgaLimitCMS13)));
32 ModelParamMap.insert(std::make_pair("obsZgaLimitATLAS", std::cref(obsZgaLimitATLAS)));
33 ModelParamMap.insert(std::make_pair("obsZgaLimitCMS", std::cref(obsZgaLimitCMS)));
34 ModelParamMap.insert(std::make_pair("expZgaLimitATLAS13", std::cref(expZgaLimitATLAS13)));
35 ModelParamMap.insert(std::make_pair("expZgaLimitCMS13", std::cref(expZgaLimitCMS13)));
36 ModelParamMap.insert(std::make_pair("expZgaLimitATLAS", std::cref(expZgaLimitATLAS)));
37 ModelParamMap.insert(std::make_pair("expZgaLimitCMS", std::cref(expZgaLimitCMS)));
38
39 loopComputed = false;
40 cg_loop = 0.;
41 cga_loop = 0.;
42 cZga_loop = 0.;
43}
44
45bool HiggsChiral::PostUpdate()
46{
47 if (!NPbase::PostUpdate()) return (false);
48
49// Check first the flags that control the values of the parameters
50
51 if (FlagUniversalcf) {
52// Assign to all cf the value of ct
53 cb = ct;
54 cc = ct;
55 cs = ct;
56 ctau = ct;
57 cmu = ct;
58 }
59
60 if (FlagUniversalcvcf) {
61// Assign to all cv and cf the value of ct
62 cv = ct;
63 cb = ct;
64 cc = ct;
65 cs = ct;
66 ctau = ct;
67 cmu = ct;
68 }
69
70// Then the add the values of the loops to the cache
71
72 loopComputed = false;
73
74 cg_loop = computecg();
75 cga_loop = computecgaga();
76 cZga_loop = computecZga();
77
78 loopComputed = true;
79
80 return (true);
81}
82
83void HiggsChiral::setParameter(const std::string name, const double& value)
84{
85 if (name.compare("cv") == 0)
86 cv = value;
87 else if (name.compare("ct") == 0)
88 ct = value;
89 else if (name.compare("cb") == 0)
90 cb = value;
91 else if (name.compare("cc") == 0)
92 cc = value;
93 else if (name.compare("cs") == 0)
94 cs = value;
95 else if (name.compare("ctau") == 0)
96 ctau = value;
97 else if (name.compare("cmu") == 0)
98 cmu = value;
99 else if (name.compare("cg") == 0)
100 cg = value;
101 else if (name.compare("cga") == 0)
102 cga = value;
103 else if (name.compare("cZga") == 0)
104 cZga = value;
105 else if (name.compare("obsZgaLimitATLAS13") == 0)
106 obsZgaLimitATLAS13 = value;
107 else if (name.compare("obsZgaLimitCMS13") == 0)
108 obsZgaLimitCMS13 = value;
109 else if (name.compare("obsZgaLimitATLAS") == 0)
110 obsZgaLimitATLAS = value;
111 else if (name.compare("obsZgaLimitCMS") == 0)
112 obsZgaLimitCMS = value;
113 else if (name.compare("expZgaLimitATLAS13") == 0)
114 expZgaLimitATLAS13 = value;
115 else if (name.compare("expZgaLimitCMS13") == 0)
116 expZgaLimitCMS13 = value;
117 else if (name.compare("expZgaLimitATLAS") == 0)
118 expZgaLimitATLAS = value;
119 else if (name.compare("expZgaLimitCMS") == 0)
120 expZgaLimitCMS = value;
121 else
122 NPbase::setParameter(name, value);
123}
124
125bool HiggsChiral::CheckParameters(const std::map<std::string, double>& DPars)
126{
127 for (int i = 0; i < NHChiralvars; i++) {
128 if (DPars.find(HChiralvars[i]) == DPars.end()) {
129 std::cout << "ERROR: missing mandatory HiggsChiral parameter " << HChiralvars[i] << std::endl;
130 raiseMissingModelParameterCount();
131 addMissingModelParameter(HChiralvars[i]);
132 }
133 }
134 return (NPbase::CheckParameters(DPars));
135}
136
137bool HiggsChiral::setFlag(const std::string name, const bool value)
138{
139 bool res = false;
140 if (name.compare("Universalcf") == 0) {
141 FlagUniversalcf = value;
142 res = true;
143 } else if (name.compare("Universalcvcf") == 0) {
144 FlagUniversalcvcf = value;
145 res = true;
146 } else
147 res = NPbase::setFlag(name, value);
148
149 return (res);
150}
151
153
154const double HiggsChiral::obliqueS() const
155{
156 double Lambda;
157 if (fabs(1.0 - cv * cv) < pow(10.0, -32.0))
158 Lambda = pow(10.0, 19.0);
159 else
160 Lambda = 4.0 * M_PI * v() / sqrt(fabs(1.0 - cv * cv));
161
162 return ( 1.0 / 12.0 / M_PI * (1.0 - cv * cv) * log(Lambda * Lambda / mHl / mHl));
163}
164
165const double HiggsChiral::obliqueT() const
166{
167 double Lambda;
168 double cW2_SM = trueSM.cW2();
169 if (fabs(1.0 - cv * cv) < pow(10.0, -32.0))
170 Lambda = pow(10.0, 19.0);
171 else
172 Lambda = 4.0 * M_PI * v() / sqrt(fabs(1.0 - cv * cv));
173
174 return ( -3.0 / 16.0 / M_PI / cW2_SM * (1.0 - cv * cv) * log(Lambda * Lambda / mHl / mHl));
175}
176
177const double HiggsChiral::obliqueU() const
178{
179 return 0.0;
180}
181
183
184const double HiggsChiral::muggH(const double sqrt_s) const
185{
186// return Gammagg() / (trueSM.computeBrHtogg()*trueSM.computeGammaHTotal());
187
188 return computecg() * computecg();
189}
190
191const double HiggsChiral::muVBF(const double sqrt_s) const
192{
193 return computecV() * computecV();
194}
195
196const double HiggsChiral::muVBFgamma(const double sqrt_s) const
197{
198 return computecV() * computecV();
199}
200
201const double HiggsChiral::mueeWBF(const double sqrt_s, const double Pol_em, const double Pol_ep) const
202{
203 return computecV() * computecV();
204}
205
206const double HiggsChiral::mueeHvv(const double sqrt_s, const double Pol_em, const double Pol_ep) const
207{
208 return computecV() * computecV();
209}
210
211const double HiggsChiral::mueeZBF(const double sqrt_s, const double Pol_em, const double Pol_ep) const
212{
213 return computecV() * computecV();
214}
215
216const double HiggsChiral::muWH(const double sqrt_s) const
217{
218 return computecV() * computecV();
219}
220
221const double HiggsChiral::muZH(const double sqrt_s) const
222{
223 return computecV() * computecV();
224}
225
226const double HiggsChiral::mueeZH(const double sqrt_s, const double Pol_em, const double Pol_ep) const
227{
228 return computecV() * computecV();
229}
230
231const double HiggsChiral::mueeZllH(const double sqrt_s, const double Pol_em, const double Pol_ep) const
232{
233 return computecV() * computecV();
234}
235
236const double HiggsChiral::mueeZqqH(const double sqrt_s, const double Pol_em, const double Pol_ep) const
237{
238 return computecV() * computecV();
239}
240
241const double HiggsChiral::muVH(const double sqrt_s) const
242{
243 return computecV() * computecV();
244}
245
246const double HiggsChiral::muVBFpVH(const double sqrt_s) const
247{
248 return computecV() * computecV();
249}
250
251const double HiggsChiral::muttH(const double sqrt_s) const
252{
253 return computect() * computect();
254}
255
256const double HiggsChiral::mutHq(const double sqrt_s) const
257{
258 double mu = 1.0;
259 double Kt = computect();
260 double Kw = computecV();
261
262// From ATLAS parametrization of the K-framework (at 14 TeV?)
263 mu = 2.63 * Kt * Kt + 3.58 * Kw * Kw - 5.21 * Kt * Kw;
264
265 return mu;
266}
267
268const double HiggsChiral::muggHpttH(const double sqrt_s) const
269{
270 double sigmaggH_SM = trueSM.computeSigmaggH(sqrt_s);
271 double sigmattH_SM = trueSM.computeSigmattH(sqrt_s);
272
273 double sigmaggH = muggH(sqrt_s) * sigmaggH_SM;
274 double sigmattH = muttH(sqrt_s) * sigmattH_SM;
275
276 return (sigmaggH + sigmattH) / (sigmaggH_SM + sigmattH_SM);
277}
278
279const double HiggsChiral::mueettH(const double sqrt_s, const double Pol_em, const double Pol_ep) const
280{
281 return computect() * computect();
282}
283
284
285const double HiggsChiral::mummH(const double sqrt_s) const
286{
287 double mu = 1.0;
288// mu = (computecmu() * computecmu());
289
290 if (sqrt_s == 0.125) {
291
292// Peak production cross section mu mu -> H -> X = 4 pi/mH^2 * BR(H->mu mu) * BR(H-> X)
293// Use mu mu -> H = 4 pi/mH^2 * BR(H->mu mu), so the xs BR formulae still applies
294 mu = BrHmumuRatio();
295
296 } else
297 throw std::runtime_error("Bad argument in HiggsChiral::mummH()");
298
299 return mu;
300}
301
302const double HiggsChiral::mummHNWA(const double sqrt_s) const
303{
304 double mu = 1.0;
305 mu = (computecmu() * computecmu());
306
307 return mu;
308}
309
310
311const double HiggsChiral::mummZH(const double sqrt_s) const
312{
313 double mu = 1.0;
314 mu = (computecV() * computecV());
315 return mu;
316}
317
318
319const double HiggsChiral::mummHvv(const double sqrt_s) const
320{
321// In this model we neglect the contributions from HZ->vv with off-shell Z
322// Dependence is exact if Custodial symmetry option is on
323 double mu = 1.0;
324 mu = (computecV() * computecV());
325 return mu;
326}
327
328
329const double HiggsChiral::mummHmm(const double sqrt_s) const
330{
331 double mu = 1.0;
332 mu = (computecV() * computecV());
333 return mu;
334}
335
336
337const double HiggsChiral::mummttH(const double sqrt_s) const
338{
339 double mu = 1.0;
340 mu = (computect() * computect());
341 return mu;
342}
343
344
345const double HiggsChiral::UpperLimitZgammaA13(const double sqrt_s) const
346{
347 double sigmaggH_SM = trueSM.computeSigmaggH(sqrt_s);
348 double sigmattH_SM = trueSM.computeSigmattH(sqrt_s);
349 double sigmaWH_SM = trueSM.computeSigmaWH(sqrt_s);
350 double sigmaZH_SM = trueSM.computeSigmaZH(sqrt_s);
351 double sigmaWF_SM = trueSM.computeSigmaWF(sqrt_s);
352 double sigmaZF_SM = trueSM.computeSigmaZF(sqrt_s);
353 double sigmaZWF_SM = trueSM.computeSigmaZWF(sqrt_s);
354 double sigmaVBF_SM = sigmaWF_SM + sigmaZF_SM + sigmaZWF_SM;
355
356 double sigmaggH = muggH(sqrt_s) * sigmaggH_SM;
357 double sigmattH = muttH(sqrt_s) * sigmattH_SM;
358 double sigmaWH = muWH(sqrt_s) * sigmaWH_SM;
359 double sigmaZH = muZH(sqrt_s) * sigmaZH_SM;
360 double sigmaVBF = muVBF(sqrt_s) * sigmaVBF_SM;
361
362 double muppH=((sigmaggH + sigmattH + sigmaWH + sigmaZH + sigmaVBF) / (sigmaggH_SM + sigmattH_SM + sigmaWH_SM + sigmaZH_SM + sigmaVBF_SM));
363// double BrHZgaRatio=;
364 double nftos=1.95996398454;
365
366 return (1.0+(muppH*BrHZgaRatio()-getobsZgaLimitATLAS13())/getexpZgaLimitATLAS13() ) * nftos;
367}
368
369const double HiggsChiral::UpperLimitZgammaC13(const double sqrt_s) const
370{
371 double sigmaggH_SM = trueSM.computeSigmaggH(sqrt_s);
372 double sigmattH_SM = trueSM.computeSigmattH(sqrt_s);
373 double sigmaWH_SM = trueSM.computeSigmaWH(sqrt_s);
374 double sigmaZH_SM = trueSM.computeSigmaZH(sqrt_s);
375 double sigmaWF_SM = trueSM.computeSigmaWF(sqrt_s);
376 double sigmaZF_SM = trueSM.computeSigmaZF(sqrt_s);
377 double sigmaZWF_SM = trueSM.computeSigmaZWF(sqrt_s);
378 double sigmaVBF_SM = sigmaWF_SM + sigmaZF_SM + sigmaZWF_SM;
379
380 double sigmaggH = muggH(sqrt_s) * sigmaggH_SM;
381 double sigmattH = muttH(sqrt_s) * sigmattH_SM;
382 double sigmaWH = muWH(sqrt_s) * sigmaWH_SM;
383 double sigmaZH = muZH(sqrt_s) * sigmaZH_SM;
384 double sigmaVBF = muVBF(sqrt_s) * sigmaVBF_SM;
385
386 double muppH=((sigmaggH + sigmattH + sigmaWH + sigmaZH + sigmaVBF) / (sigmaggH_SM + sigmattH_SM + sigmaWH_SM + sigmaZH_SM + sigmaVBF_SM));
387// double BrHZgaRatio=;
388 double nftos=1.95996398454;
389
390 return (1.0+(muppH*BrHZgaRatio()-getobsZgaLimitCMS13())/getexpZgaLimitCMS13() ) * nftos;
391}
392
393const double HiggsChiral::UpperLimitZgammaA(const double sqrt_s) const
394{
395 double sigmaggH_SM = trueSM.computeSigmaggH(sqrt_s);
396 double sigmattH_SM = trueSM.computeSigmattH(sqrt_s);
397 double sigmaWH_SM = trueSM.computeSigmaWH(sqrt_s);
398 double sigmaZH_SM = trueSM.computeSigmaZH(sqrt_s);
399 double sigmaWF_SM = trueSM.computeSigmaWF(sqrt_s);
400 double sigmaZF_SM = trueSM.computeSigmaZF(sqrt_s);
401 double sigmaZWF_SM = trueSM.computeSigmaZWF(sqrt_s);
402 double sigmaVBF_SM = sigmaWF_SM + sigmaZF_SM + sigmaZWF_SM;
403
404 double sigmaggH = muggH(sqrt_s) * sigmaggH_SM;
405 double sigmattH = muttH(sqrt_s) * sigmattH_SM;
406 double sigmaWH = muWH(sqrt_s) * sigmaWH_SM;
407 double sigmaZH = muZH(sqrt_s) * sigmaZH_SM;
408 double sigmaVBF = muVBF(sqrt_s) * sigmaVBF_SM;
409
410 double muppH=((sigmaggH + sigmattH + sigmaWH + sigmaZH + sigmaVBF) / (sigmaggH_SM + sigmattH_SM + sigmaWH_SM + sigmaZH_SM + sigmaVBF_SM));
411// double BrHZgaRatio=;
412 double nftos=1.95996398454;
413
414 return (1.0+(muppH*BrHZgaRatio()-getobsZgaLimitATLAS())/getexpZgaLimitATLAS() ) * nftos;
415}
416
417const double HiggsChiral::UpperLimitZgammaC(const double sqrt_s) const
418{
419 double sigmaggH_SM = trueSM.computeSigmaggH(sqrt_s);
420 double sigmattH_SM = trueSM.computeSigmattH(sqrt_s);
421 double sigmaWH_SM = trueSM.computeSigmaWH(sqrt_s);
422 double sigmaZH_SM = trueSM.computeSigmaZH(sqrt_s);
423 double sigmaWF_SM = trueSM.computeSigmaWF(sqrt_s);
424 double sigmaZF_SM = trueSM.computeSigmaZF(sqrt_s);
425 double sigmaZWF_SM = trueSM.computeSigmaZWF(sqrt_s);
426 double sigmaVBF_SM = sigmaWF_SM + sigmaZF_SM + sigmaZWF_SM;
427
428 double sigmaggH = muggH(sqrt_s) * sigmaggH_SM;
429 double sigmattH = muttH(sqrt_s) * sigmattH_SM;
430 double sigmaWH = muWH(sqrt_s) * sigmaWH_SM;
431 double sigmaZH = muZH(sqrt_s) * sigmaZH_SM;
432 double sigmaVBF = muVBF(sqrt_s) * sigmaVBF_SM;
433
434 double muppH=((sigmaggH + sigmattH + sigmaWH + sigmaZH + sigmaVBF) / (sigmaggH_SM + sigmattH_SM + sigmaWH_SM + sigmaZH_SM + sigmaVBF_SM));
435// double BrHZgaRatio=;
436 double nftos=1.95996398454;
437
438 return (1.0+(muppH*BrHZgaRatio()-getobsZgaLimitCMS())/getexpZgaLimitCMS() ) * nftos;
439}
440
441const double HiggsChiral::cgplusct() const
442{
443 return cg + ct;
444}
445
446const double HiggsChiral::cgaplusct() const
447{
448 return cga + ct;
449}
450
451const double HiggsChiral::cgminuscga() const
452{
453 return cg - cga;
454}
455
456const double HiggsChiral::cVpluscb() const
457{
458 return cv + cb;
459}
460
461const double HiggsChiral::cVplusctau() const
462{
463 return cv + ctau;
464}
465
466const double HiggsChiral::cbminuscc() const
467{
468 return cb - cc;
469}
470
471const double HiggsChiral::cbminusctau() const
472{
473 return cb - ctau;
474}
475
476const double HiggsChiral::ccminusctau() const
477{
478 return cc - ctau;
479}
480
481const double HiggsChiral::Gammagg() const
482{
483 double Mt=trueSM.getQuarks(QCD::TOP).getMass();
484 double Mb=trueSM.getQuarks(QCD::BOTTOM).getMass();
485 double Mc=trueSM.getQuarks(QCD::CHARM).getMass();
486 double Ms=trueSM.getQuarks(QCD::STRANGE).getMass();
487 double TAUt=4.0*Mt*Mt/(mHl*mHl);
488 double TAUb=4.0*Mb*Mb/(mHl*mHl);
489 double TAUc=4.0*Mc*Mc/(mHl*mHl);
490 double TAUs=4.0*Ms*Ms/(mHl*mHl);
491
492 double factor = GF*AlsMz*AlsMz*mHl*mHl*mHl/(sqrt(2.0)*4.0*M_PI*M_PI*M_PI);
493
494 double ksoftNLO = 1.0 + (AlsMz/M_PI) * ( 73.0/4.0 -7.0*5.0/6.0);
495
496 return factor * ksoftNLO * ( ( ct * 0.5*TAUt*(1.0+(1.0-TAUt)*f_func(TAUt))
497 +cb * 0.5*TAUb*(1.0+(1.0-TAUb)*f_func(TAUb))
498 +cc * 0.5*TAUc*(1.0+(1.0-TAUc)*f_func(TAUc))
499 +cs * 0.5*TAUs*(1.0+(1.0-TAUs)*f_func(TAUs))) * (1.0+11.0*AlsMz/(4.0*M_PI))
500 +cg * 0.5 ).abs2();
501}
502
503const double HiggsChiral::GammaWW() const
504{
505 return computecV() * computecV() * trueSM.computeBrHtoWW() * trueSM.computeGammaHTotal();
506}
507
508const double HiggsChiral::GammaZZ() const
509{
510 return computecV() * computecV() * trueSM.computeBrHtoZZ() * trueSM.computeGammaHTotal();
511}
512
513const double HiggsChiral::GammaZga() const
514{
515 double Mt=trueSM.getQuarks(QCD::TOP).getMass();
516 double Mb=trueSM.getQuarks(QCD::BOTTOM).getMass();
517 double Mc=trueSM.getQuarks(QCD::CHARM).getMass();
518 double Ms=trueSM.getQuarks(QCD::STRANGE).getMass();
519 double Mtau=trueSM.getLeptons(StandardModel::TAU).getMass();
520 double MW=trueSM.Mw_tree();
521 double cW2=trueSM.c02();
522 double sW2=1.0-cW2;
523 double TAUt=4.0*Mt*Mt/(mHl*mHl);
524 double TAUb=4.0*Mb*Mb/(mHl*mHl);
525 double TAUc=4.0*Mc*Mc/(mHl*mHl);
526 double TAUs=4.0*Ms*Ms/(mHl*mHl);
527 double TAUtau=4.0*Mtau*Mtau/(mHl*mHl);
528 double TAUw=4.0*MW*MW/(mHl*mHl);
529 double LAMt=4.0*Mt*Mt/(Mz*Mz);
530 double LAMb=4.0*Mb*Mb/(Mz*Mz);
531 double LAMc=4.0*Mc*Mc/(Mz*Mz);
532 double LAMs=4.0*Ms*Ms/(Mz*Mz);
533 double LAMtau=4.0*Mtau*Mtau/(Mz*Mz);
534 double LAMw=4.0*MW*MW/(Mz*Mz);
535
536 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);
537
538 return factor * ((-ct * 4.0*(0.5-4.0/3.0*sW2)*(Int1(TAUt,LAMt)-Int2(TAUt,LAMt)) * (1.0-AlsMz/M_PI)
539 +cb * 2.0*(-0.5+2.0/3.0*sW2)*(Int1(TAUb,LAMb)-Int2(TAUb,LAMb))
540 -cc * 4.0*(0.5-4.0/3.0*sW2)*(Int1(TAUc,LAMc)-Int2(TAUc,LAMc))
541 +cs * 2.0*(-0.5+2.0/3.0*sW2)*(Int1(TAUs,LAMs)-Int2(TAUs,LAMs))
542 +ctau * 2.0*(-0.5+2.0*sW2)*(Int1(TAUtau,LAMtau)-Int2(TAUtau,LAMtau)) )/sqrt(sW2*cW2)
543 -cv * sqrt(cW2/sW2)*(4.0*(3.0-sW2/cW2)*Int2(TAUw,LAMw)
544 +((1.0+2.0/TAUw)*sW2/cW2-(5.0+2.0/TAUw))*Int1(TAUw,LAMw))
545 +cZga).abs2(); /*CHANGE THE cZga COEFFICIENT!*/
546}
547
548const double HiggsChiral::Gammagaga() const
549{
550 double Mt=trueSM.getQuarks(QCD::TOP).getMass();
551 double Mb=trueSM.getQuarks(QCD::BOTTOM).getMass();
552 double Mc=trueSM.getQuarks(QCD::CHARM).getMass();
553 double Ms=trueSM.getQuarks(QCD::STRANGE).getMass();
554 double Mtau=trueSM.getLeptons(StandardModel::TAU).getMass();
555 double MW=trueSM.Mw_tree();
556 double TAUt=4.0*Mt*Mt/(mHl*mHl);
557 double TAUb=4.0*Mb*Mb/(mHl*mHl);
558 double TAUc=4.0*Mc*Mc/(mHl*mHl);
559 double TAUs=4.0*Ms*Ms/(mHl*mHl);
560 double TAUtau=4.0*Mtau*Mtau/(mHl*mHl);
561 double TAUw=4.0*MW*MW/(mHl*mHl);
562
563 double factor = GF*ale*ale*mHl*mHl*mHl/(sqrt(2.0)*128.0*M_PI*M_PI*M_PI);
564
565 return factor * ( ct * (8./3.)*TAUt*(1.+(1.-TAUt)*f_func(TAUt)) * (1.0-AlsMz/M_PI)
566 +cb * (2./3.)*TAUb*(1.+(1.-TAUb)*f_func(TAUb))
567 +cc * (8./3.)*TAUc*(1.+(1.-TAUc)*f_func(TAUc))
568 +cs * (2./3.)*TAUs*(1.+(1.-TAUs)*f_func(TAUs))
569 +ctau * 2.0*TAUtau*(1.+(1.-TAUtau)*f_func(TAUtau))
570 -cv * (2.0+3.0*TAUw+3.0*TAUw*(2.0-TAUw)*f_func(TAUw))
571 +cga * 2.0 ).abs2();
572}
573
574const double HiggsChiral::Gammamumu() const
575{
576 return computecmu() * computecmu() * trueSM.computeBrHtomumu() * trueSM.computeGammaHTotal();
577}
578
579const double HiggsChiral::Gammatautau() const
580{
581 return computectau() * computectau() * trueSM.computeBrHtotautau() * trueSM.computeGammaHTotal();
582}
583
584const double HiggsChiral::Gammacc() const
585{
586 return computecc() * computecc() * trueSM.computeBrHtocc() * trueSM.computeGammaHTotal();
587}
588
589const double HiggsChiral::Gammass() const
590{
591 return computecs() * computecs() * trueSM.computeBrHtoss() * trueSM.computeGammaHTotal();
592}
593
594const double HiggsChiral::Gammabb() const
595{
596 return computecb() * computecb() * trueSM.computeBrHtobb() * trueSM.computeGammaHTotal();
597}
598
599const double HiggsChiral::GammaTotal() const
600{
601 return Gammagg() + GammaWW() + GammaZZ() + GammaZga() + Gammagaga()
602 + Gammamumu() + Gammatautau() + Gammacc() + Gammass() + Gammabb();
603}
604
605const double HiggsChiral::BrHggRatio() const
606{
607// return Gammagg() / GammaTotal() / trueSM.computeBrHtogg();
608
609 return (computecg() * computecg() / computeGammaTotalRatio());
610}
611
612const double HiggsChiral::BrHWWRatio() const
613{
614// return GammaWW() / GammaTotal() / trueSM.computeBrHtoWW();
615
616 return (computecV() * computecV() / computeGammaTotalRatio());
617}
618
619const double HiggsChiral::BrHZZRatio() const
620{
621// return GammaZZ() / GammaTotal() / trueSM.computeBrHtoZZ();
622
623 return (computecV() * computecV() / computeGammaTotalRatio());
624}
625
626const double HiggsChiral::BrHVVRatio() const
627{
628 double BrZZSM = trueSM.computeBrHtoZZ(), BrWWSM = trueSM.computeBrHtoWW();
629
630 return (BrZZSM * BrHZZRatio() + BrWWSM * BrHWWRatio() )/(BrZZSM + BrWWSM);
631}
632
633const double HiggsChiral::BrHZgaRatio() const
634{
635// return GammaZga() / GammaTotal() / trueSM.computeBrHtoZga();
636
637 return (computecZga() * computecZga() / computeGammaTotalRatio());
638}
639
640const double HiggsChiral::BrHZgallRatio() const
641{
642// return GammaZga() / GammaTotal() / trueSM.computeBrHtoZga();
643
644 return BrHZgaRatio();
645}
646
647const double HiggsChiral::BrHZgaeeRatio() const
648{
649// return GammaZga() / GammaTotal() / trueSM.computeBrHtoZga();
650
651 return BrHZgaRatio();
652}
653
654const double HiggsChiral::BrHZgamumuRatio() const
655{
656// return GammaZga() / GammaTotal() / trueSM.computeBrHtoZga();
657
658 return BrHZgaRatio();
659}
660
661const double HiggsChiral::BrHgagaRatio() const
662{
663// return Gammagaga() / GammaTotal() / trueSM.computeBrHtogaga();
664
665 return (computecgaga() * computecgaga() / computeGammaTotalRatio());
666}
667
668const double HiggsChiral::BrHmumuRatio() const
669{
670// return Gammamumu() / GammaTotal() / trueSM.computeBrHtomumu();
671
672 return (computecmu() * computecmu() / computeGammaTotalRatio());
673}
674
675const double HiggsChiral::BrHtautauRatio() const
676{
677// return Gammatautau() / GammaTotal() / trueSM.computeBrHtotautau();
678
679 return (computectau() * computectau() / computeGammaTotalRatio());
680}
681
682const double HiggsChiral::BrHccRatio() const
683{
684// return Gammacc() / GammaTotal() / trueSM.computeBrHtocc();
685
686 return (computecc() * computecc() / computeGammaTotalRatio());
687}
688
689const double HiggsChiral::BrHssRatio() const
690{
691// return Gammass() / GammaTotal() / trueSM.computeBrHtocc();
692
693 return (computecs() * computecs() / computeGammaTotalRatio());
694}
695
696const double HiggsChiral::BrHbbRatio() const
697{
698// return Gammabb() / GammaTotal() / trueSM.computeBrHtobb();
699
700 return (computecb() * computecb() / computeGammaTotalRatio());
701}
702
704
705const double HiggsChiral::BrH2L2LRatio() const
706{
707 return BrHZZRatio();
708}
709
710const double HiggsChiral::BrH2e2muRatio() const
711{
712 return BrHZZRatio();
713}
714
715const double HiggsChiral::BrH2v2vRatio() const
716{
717 return BrHZZRatio();
718}
719
720const double HiggsChiral::BrH2L2vRatio() const
721{
722 return BrHZZRatio();
723}
724
725const double HiggsChiral::BrH2L2v2Ratio() const
726{
727 return BrHZZRatio();
728}
729
730const double HiggsChiral::BrH2e2vRatio() const
731{
732 return BrHZZRatio();
733}
734
735const double HiggsChiral::BrH2mu2vRatio() const
736{
737 return BrHZZRatio();
738}
739
740const double HiggsChiral::BrH2u2uRatio() const
741{
742 return BrHZZRatio();
743}
744
745const double HiggsChiral::BrH2d2dRatio() const
746{
747 return BrHZZRatio();
748}
749
750const double HiggsChiral::BrH2u2dRatio() const
751{
752 return BrHZZRatio();
753}
754
755const double HiggsChiral::BrH2L2uRatio() const
756{
757 return BrHZZRatio();
758}
759
760const double HiggsChiral::BrH2L2dRatio() const
761{
762 return BrHZZRatio();
763}
764
765const double HiggsChiral::BrH2v2uRatio() const
766{
767 return BrHZZRatio();
768}
769
770const double HiggsChiral::BrH2v2dRatio() const
771{
772 return BrHZZRatio();
773}
774
775const double HiggsChiral::BrH4LRatio() const
776{
777 return BrHZZRatio();
778}
779
780const double HiggsChiral::BrH4L2Ratio() const
781{
782 return BrHZZRatio();
783}
784
785const double HiggsChiral::BrH4eRatio() const
786{
787 return BrHZZRatio();
788}
789
790const double HiggsChiral::BrH4muRatio() const
791{
792 return BrHZZRatio();
793}
794
795const double HiggsChiral::BrH4vRatio() const
796{
797 return BrHZZRatio();
798}
799
800const double HiggsChiral::BrH4uRatio() const
801{
802 return BrHZZRatio();
803}
804
805const double HiggsChiral::BrH4dRatio() const
806{
807 return BrHZZRatio();
808}
809
810const double HiggsChiral::BrHLvvLRatio() const
811{
812 return BrHWWRatio();
813}
814
815const double HiggsChiral::BrHevmuvRatio() const
816{
817 return BrHWWRatio();
818}
819
820const double HiggsChiral::BrHudduRatio() const
821{
822 return BrHWWRatio();
823}
824
825const double HiggsChiral::BrHLvudRatio() const
826{
827 return BrHWWRatio();
828}
829
830const double HiggsChiral::BrH2udRatio() const
831{
832 // Approximate by pure CC
833 return BrHWWRatio();
834}
835
836const double HiggsChiral::BrH2LvRatio() const
837{
838 // Approximate by pure CC
839 return BrHWWRatio();
840}
841
842const double HiggsChiral::BrH2Lv2Ratio() const
843{
844 // Approximate by pure CC
845 return BrHWWRatio();
846}
847
848const double HiggsChiral::BrH2evRatio() const
849{
850 // Approximate by pure CC
851 return BrHWWRatio();
852}
853
854const double HiggsChiral::BrH2muvRatio() const
855{
856 // Approximate by pure CC
857 return BrHWWRatio();
858}
859
860const double HiggsChiral::BrH4fRatio() const
861{
862 double Br = 0.0;
863
864 // SM decay widths (from MG simulations)
865 double wH2L2LSM=0.65682e-06, wH2v2vSM=0.28126e-05, wH2L2vSM=0.27224e-05;
866 double wH2u2uSM=0.22500e-05, wH2d2dSM=0.11906e-04, wH2u2dSM=0.12361e-04;
867 double wH2L2uSM=0.45029e-05, wH2L2dSM=0.85830e-05, wH2v2uSM=0.93233e-05;
868 double wH2v2dSM=0.17794e-04, wH4LSM=0.33973e-06, wH4vSM=0.16884e-05;
869 double wH4uSM=0.23669e-05, wH4dSM=0.60254e-05;
870 double wHLvvLSM=0.58098e-04, wHudduSM=0.13384e-03, wHLvudSM=0.34149e-03 ;
871 double wH2udSM=0.13711e-03, wH2LvSM=0.27557e-04;
872
873 // Sum: SM
874 double wH4fSM=wH2L2LSM+wH2v2vSM+wH2L2vSM+wH2u2uSM+wH2d2dSM+wH2u2dSM+
875 wH2L2uSM+wH2L2dSM+wH2v2uSM+wH2v2dSM+wH4LSM+wH4vSM+wH4uSM+wH4dSM+wHLvvLSM+wHudduSM+
876 wHLvudSM+wH2udSM+wH2LvSM;
877
878 // Sum: SM+NP
879 double wH4f= (computecV() * computecV()) * (wH2L2LSM+wH2v2vSM+wH2L2vSM+wH2u2uSM+wH2d2dSM+wH2u2dSM
880 + wH2L2uSM+wH2L2dSM+wH2v2uSM+wH2v2dSM+wH4LSM+wH4vSM+wH4uSM+wH4dSM)
881 + (computecV() * computecV()) * (wHLvvLSM+wHudduSM
882 + wHLvudSM+wH2udSM+wH2LvSM);
883
884 Br = (wH4f/wH4fSM) / computeGammaTotalRatio();
885
886 return Br;
887}
888
889 // DECAYS INVOLVING ONLY ELECTRONS, MUONS OR NEUTRINOS IN THE FINAL STATES
890
891const double HiggsChiral::BrH4lRatio() const
892{
893 return BrHZZRatio();
894}
895
896const double HiggsChiral::BrH2l2vRatio() const
897{
898 // Approximate by pure CC
899 return BrHWWRatio();
900}
901
902const double HiggsChiral::BrHtoinvRatio() const
903{
904// No extra decays so it scales as the H->ZZ
905
906 return BrHZZRatio();
907}
908
909const double HiggsChiral::muggHgaga(const double sqrt_s) const
910{
911 return computecg() * computecg() * computecgaga() * computecgaga() / computeGammaTotalRatio();
912}
913
914const double HiggsChiral::muggHgagaInt(const double sqrt_s) const
915{
916 double muNWA, GammaRatio, fki;
917 double kt,kb,kc,ks,ku,kd,kta,km,ke;
918 double kw;
919 double kg,kgamma;
920
921 muNWA = muggH(sqrt_s)*BrHgagaRatio();
922
923 kt = computect();
924 kb = computecb();
925 kc = computecc();
926 ks = 1.0;
927 ku = 1.0;
928 kd = 1.0;
929 kta = computectau();
930 km = computecmu();
931 ke = 1.0;
932
933 kw = computecV();
934
935 kg = (3.0/2.0) * computecg();
936 kgamma = (3.0/2.0) * computecgaga();
937
938 GammaRatio = computeGammaTotalRatio();
939
940 fki = 0.000802422 *kb*kb + 0.000312884 *kb*kc + 0.0000182107 *kc*kc +
941 5.94769e-9 *kb*kd + 9.62554e-10 *kc*kd + 6.2785e-15 *kd*kd +
942 5.53251e-10 *kb*ke + 3.51863e-11 *kc*ke + 1.09243e-15 *kd*ke -
943 0.00905016 *kb*kg - 0.00190706 *kc*kg - 5.97591e-9 *kd*kg -
944 6.72288e-10 *ke*kg - 0.0271505 *kb*kgamma - 0.00143029 *kc*kgamma -
945 1.79277e-8 *kd*kgamma - 0.174392 *kg*kgamma + 8.97565e-6 *kb*km +
946 6.21013e-7 *kc*km + 2.37616e-11 *kd*km - 0.0000460022 *kg*km +
947 4.93348e-6 *kb*ks + 8.51176e-7 *kc*ks + 1.29558e-11 *kd*ks +
948 1.16267e-12 *ke*ks - 0.0000123381 *kg*ks - 0.0000370143 *kgamma*ks +
949 2.22544e-8 *km*ks + 6.08665e-9 *ks*ks - 0.0467672 *kb*kt -
950 0.00394193 *kc*kt - 3.08808e-8 *kd*kt - 6.94817e-10 *ke*kt -
951 0.240315 *kg*kt - 0.180236 *kgamma*kt - 0.0000475437 *km*kt -
952 0.0000637578 *ks*kt - 0.248368 *kt*kt + 0.00100168 *kb*kta +
953 0.0000759092 *kc*kta + 3.44671e-9 *kd*kta - 0.00975386 *kg*kta +
954 2.93009e-6 *ks*kta - 0.0100807 *kt*kta + 5.30126e-8 *kb*ku +
955 5.54256e-9 *kc*ku + 1.15815e-13 *kd*ku + 4.05206e-15 *ke*ku -
956 1.03323e-7 *kg*ku - 7.74926e-8 *kgamma*ku +
957 8.62762e-11 *km*ku + 1.17664e-10 *ks*ku - 2.13572e-7 *kt*ku +
958 1.2332e-8 *kta*ku + 3.40922e-13 *ku*ku + 0.169912 *kb*kw +
959 0.00895098 *kc*kw + 1.12194e-7 *kd*kw + 1.09137 *kg*kw +
960 0.000231641 *ks*kw + 1.12795 *kt*kw + 4.8496e-7 *ku*kw;
961
962 return (muNWA - 0.022 *GammaRatio * fki)/0.978;
963
964}
965
966const double HiggsChiral::muVBFHgaga(const double sqrt_s) const
967{
968 return computecV() * computecV() * computecgaga() * computecgaga() / computeGammaTotalRatio();
969}
970
971const double HiggsChiral::muZHgaga(const double sqrt_s) const
972{
973 return computecV() * computecV() * computecgaga() * computecgaga() / computeGammaTotalRatio();
974}
975
976const double HiggsChiral::muWHgaga(const double sqrt_s) const
977{
978 return computecV() * computecV() * computecgaga() * computecgaga() / computeGammaTotalRatio();
979}
980
981const double HiggsChiral::muVHgaga(const double sqrt_s) const
982{
983 return computecV() * computecV() * computecgaga() * computecgaga() / computeGammaTotalRatio();
984}
985
986const double HiggsChiral::muttHgaga(const double sqrt_s) const
987{
988 return computect() * computect() * computecgaga() * computecgaga() / computeGammaTotalRatio();
989}
990
991const double HiggsChiral::muggHZga(const double sqrt_s) const
992{
993 return computecg() * computecg() * computecZga() * computecZga() / computeGammaTotalRatio();
994}
995
996const double HiggsChiral::muVBFHZga(const double sqrt_s) const
997{
998 return computecV() * computecV() * computecZga() * computecZga() / computeGammaTotalRatio();
999}
1000
1001const double HiggsChiral::muZHZga(const double sqrt_s) const
1002{
1003 return computecV() * computecV() * computecZga() * computecZga() / computeGammaTotalRatio();
1004}
1005
1006const double HiggsChiral::muWHZga(const double sqrt_s) const
1007{
1008 return computecV() * computecV() * computecZga() * computecZga() / computeGammaTotalRatio();
1009}
1010
1011const double HiggsChiral::muVHZga(const double sqrt_s) const
1012{
1013 return computecV() * computecV() * computecZga() * computecZga() / computeGammaTotalRatio();
1014}
1015
1016const double HiggsChiral::muttHZga(const double sqrt_s) const
1017{
1018 return computect() * computect() * computecZga() * computecZga() / computeGammaTotalRatio();
1019}
1020
1021const double HiggsChiral::muggHZZ(const double sqrt_s) const
1022{
1023 return computecg() * computecg() * computecV() * computecV() / computeGammaTotalRatio();
1024}
1025
1026const double HiggsChiral::muVBFHZZ(const double sqrt_s) const
1027{
1028 return computecV() * computecV() * computecV() * computecV() / computeGammaTotalRatio();
1029}
1030
1031const double HiggsChiral::muZHZZ(const double sqrt_s) const
1032{
1033 return computecV() * computecV() * computecV() * computecV() / computeGammaTotalRatio();
1034}
1035
1036const double HiggsChiral::muWHZZ(const double sqrt_s) const
1037{
1038 return computecV() * computecV() * computecV() * computecV() / computeGammaTotalRatio();
1039}
1040
1041const double HiggsChiral::muVHZZ(const double sqrt_s) const
1042{
1043 return computecV() * computecV() * computecV() * computecV() / computeGammaTotalRatio();
1044}
1045
1046const double HiggsChiral::muttHZZ(const double sqrt_s) const
1047{
1048 return computect() * computect() * computecV() * computecV() / computeGammaTotalRatio();
1049}
1050
1051const double HiggsChiral::muggHZZ4l(const double sqrt_s) const
1052{
1053 return computecg() * computecg() * computecV() * computecV() / computeGammaTotalRatio();
1054}
1055
1056const double HiggsChiral::muVBFHZZ4l(const double sqrt_s) const
1057{
1058 return computecV() * computecV() * computecV() * computecV() / computeGammaTotalRatio();
1059}
1060
1061const double HiggsChiral::muZHZZ4l(const double sqrt_s) const
1062{
1063 return computecV() * computecV() * computecV() * computecV() / computeGammaTotalRatio();
1064}
1065
1066const double HiggsChiral::muWHZZ4l(const double sqrt_s) const
1067{
1068 return computecV() * computecV() * computecV() * computecV() / computeGammaTotalRatio();
1069}
1070
1071const double HiggsChiral::muVHZZ4l(const double sqrt_s) const
1072{
1073 return computecV() * computecV() * computecV() * computecV() / computeGammaTotalRatio();
1074}
1075
1076const double HiggsChiral::muttHZZ4l(const double sqrt_s) const
1077{
1078 return computect() * computect() * computecV() * computecV() / computeGammaTotalRatio();
1079}
1080
1081const double HiggsChiral::muggHWW(const double sqrt_s) const
1082{
1083 return computecg() * computecg() * computecV() * computecV() / computeGammaTotalRatio();
1084}
1085
1086const double HiggsChiral::muVBFHWW(const double sqrt_s) const
1087{
1088 return computecV() * computecV() * computecV() * computecV() / computeGammaTotalRatio();
1089}
1090
1091const double HiggsChiral::muZHWW(const double sqrt_s) const
1092{
1093 return computecV() * computecV() * computecV() * computecV() / computeGammaTotalRatio();
1094}
1095
1096const double HiggsChiral::muWHWW(const double sqrt_s) const
1097{
1098 return computecV() * computecV() * computecV() * computecV() / computeGammaTotalRatio();
1099}
1100
1101const double HiggsChiral::muVHWW(const double sqrt_s) const
1102{
1103 return computecV() * computecV() * computecV() * computecV() / computeGammaTotalRatio();
1104}
1105
1106const double HiggsChiral::muttHWW(const double sqrt_s) const
1107{
1108 return computect() * computect() * computecV() * computecV() / computeGammaTotalRatio();
1109}
1110
1111const double HiggsChiral::muggHWW2l2v(const double sqrt_s) const
1112{
1113 return computecg() * computecg() * computecV() * computecV() / computeGammaTotalRatio();
1114}
1115
1116const double HiggsChiral::muVBFHWW2l2v(const double sqrt_s) const
1117{
1118 return computecV() * computecV() * computecV() * computecV() / computeGammaTotalRatio();
1119}
1120
1121const double HiggsChiral::muZHWW2l2v(const double sqrt_s) const
1122{
1123 return computecV() * computecV() * computecV() * computecV() / computeGammaTotalRatio();
1124}
1125
1126const double HiggsChiral::muWHWW2l2v(const double sqrt_s) const
1127{
1128 return computecV() * computecV() * computecV() * computecV() / computeGammaTotalRatio();
1129}
1130
1131const double HiggsChiral::muVHWW2l2v(const double sqrt_s) const
1132{
1133 return computecV() * computecV() * computecV() * computecV() / computeGammaTotalRatio();
1134}
1135
1136const double HiggsChiral::muttHWW2l2v(const double sqrt_s) const
1137{
1138 return computect() * computect() * computecV() * computecV() / computeGammaTotalRatio();
1139}
1140
1141const double HiggsChiral::muggHmumu(const double sqrt_s) const
1142{
1143 return computecg() * computecg() * computecmu() * computecmu() / computeGammaTotalRatio();
1144}
1145
1146const double HiggsChiral::muVBFHmumu(const double sqrt_s) const
1147{
1148 return computecV() * computecV() * computecmu() * computecmu() / computeGammaTotalRatio();
1149}
1150
1151const double HiggsChiral::muZHmumu(const double sqrt_s) const
1152{
1153 return computecV() * computecV() * computecmu() * computecmu() / computeGammaTotalRatio();
1154}
1155
1156const double HiggsChiral::muWHmumu(const double sqrt_s) const
1157{
1158 return computecV() * computecV() * computecmu() * computecmu() / computeGammaTotalRatio();
1159}
1160
1161const double HiggsChiral::muVHmumu(const double sqrt_s) const
1162{
1163 return computecV() * computecV() * computecmu() * computecmu() / computeGammaTotalRatio();
1164}
1165
1166const double HiggsChiral::muttHmumu(const double sqrt_s) const
1167{
1168 return computect() * computect() * computecmu() * computecmu() / computeGammaTotalRatio();
1169}
1170
1171const double HiggsChiral::muggHtautau(const double sqrt_s) const
1172{
1173 return computecg() * computecg() * computectau() * computectau() / computeGammaTotalRatio();
1174}
1175
1176const double HiggsChiral::muVBFHtautau(const double sqrt_s) const
1177{
1178 return computecV() * computecV() * computectau() * computectau() / computeGammaTotalRatio();
1179}
1180
1181const double HiggsChiral::muZHtautau(const double sqrt_s) const
1182{
1183 return computecV() * computecV() * computectau() * computectau() / computeGammaTotalRatio();
1184}
1185
1186const double HiggsChiral::muWHtautau(const double sqrt_s) const
1187{
1188 return computecV() * computecV() * computectau() * computectau() / computeGammaTotalRatio();
1189}
1190
1191const double HiggsChiral::muVHtautau(const double sqrt_s) const
1192{
1193 return computecV() * computecV() * computectau() * computectau() / computeGammaTotalRatio();
1194}
1195
1196const double HiggsChiral::muttHtautau(const double sqrt_s) const
1197{
1198 return computect() * computect() * computectau() * computectau() / computeGammaTotalRatio();
1199}
1200
1201const double HiggsChiral::muggHbb(const double sqrt_s) const
1202{
1203 return computecg() * computecg() * computecb() * computecb() / computeGammaTotalRatio();
1204}
1205
1206const double HiggsChiral::muVBFHbb(const double sqrt_s) const
1207{
1208 return computecV() * computecV() * computecb() * computecb() / computeGammaTotalRatio();
1209}
1210
1211const double HiggsChiral::muZHbb(const double sqrt_s) const
1212{
1213 return computecV() * computecV() * computecb() * computecb() / computeGammaTotalRatio();
1214}
1215
1216const double HiggsChiral::muWHbb(const double sqrt_s) const
1217{
1218 return computecV() * computecV() * computecb() * computecb() / computeGammaTotalRatio();
1219}
1220
1221const double HiggsChiral::muVHbb(const double sqrt_s) const
1222{
1223 return computecV() * computecV() * computecb() * computecb() / computeGammaTotalRatio();
1224}
1225
1226const double HiggsChiral::muttHbb(const double sqrt_s) const
1227{
1228 return computect() * computect() * computecb() * computecb() / computeGammaTotalRatio();
1229}
1230
1231const double HiggsChiral::muppHmumu(const double sqrt_s) const
1232{
1233 if(sqrt_s==8)
1234 {
1235 return (0.872 * computecg() * computecg() + 0.122 * computecV() * computecV() + 0.006 * computect() * computect()) * computecmu() * computecmu() / computeGammaTotalRatio();
1236 }
1237 if(sqrt_s==13)
1238 {
1239 return (0.871 * computecg() * computecg() + 0.119 * computecV() * computecV() + 0.010 * computect() * computect()) * computecmu() * computecmu() / computeGammaTotalRatio();
1240 }
1241 else
1242 {
1243 throw std::runtime_error("The observable muppHmumu is only defined for 8 or 13 TeV.");
1244 }
1245}
1246
1247
1249//-----------------------------------------------------------------------------------------
1250//-- Special Hadron collider signal strengths with separate full TH unc U(prod x decay) ---
1251//-----------------------------------------------------------------------------------------
1253
1254const double HiggsChiral::muTHUggHgaga(const double sqrt_s) const
1255{
1256 return muggH(sqrt_s)*BrHgagaRatio();
1257}
1258
1259const double HiggsChiral::muTHUVBFHgaga(const double sqrt_s) const
1260{
1261 return muVBF(sqrt_s)*BrHgagaRatio();
1262}
1263
1264const double HiggsChiral::muTHUZHgaga(const double sqrt_s) const
1265{
1266 return muZH(sqrt_s)*BrHgagaRatio();
1267}
1268
1269const double HiggsChiral::muTHUWHgaga(const double sqrt_s) const
1270{
1271 return muWH(sqrt_s)*BrHgagaRatio();
1272}
1273
1274const double HiggsChiral::muTHUVHgaga(const double sqrt_s) const
1275{
1276 return muVH(sqrt_s)*BrHgagaRatio();
1277}
1278
1279const double HiggsChiral::muTHUttHgaga(const double sqrt_s) const
1280{
1281 return muttH(sqrt_s)*BrHgagaRatio();
1282}
1283
1284const double HiggsChiral::muTHUggHZga(const double sqrt_s) const
1285{
1286 return muggH(sqrt_s)*BrHZgaRatio();
1287}
1288
1289const double HiggsChiral::muTHUVBFHZga(const double sqrt_s) const
1290{
1291 return muVBF(sqrt_s)*BrHZgaRatio();
1292}
1293
1294const double HiggsChiral::muTHUZHZga(const double sqrt_s) const
1295{
1296 return muZH(sqrt_s)*BrHZgaRatio();
1297}
1298
1299const double HiggsChiral::muTHUWHZga(const double sqrt_s) const
1300{
1301 return muWH(sqrt_s)*BrHZgaRatio();
1302}
1303
1304const double HiggsChiral::muTHUVHZga(const double sqrt_s) const
1305{
1306 return muVH(sqrt_s)*BrHZgaRatio();
1307}
1308
1309const double HiggsChiral::muTHUttHZga(const double sqrt_s) const
1310{
1311 return muttH(sqrt_s)*BrHZgaRatio();
1312}
1313
1314const double HiggsChiral::muTHUggHZZ(const double sqrt_s) const
1315{
1316 return muggH(sqrt_s)*BrHZZRatio();
1317}
1318
1319const double HiggsChiral::muTHUVBFHZZ(const double sqrt_s) const
1320{
1321 return muVBF(sqrt_s)*BrHZZRatio();
1322}
1323
1324const double HiggsChiral::muTHUZHZZ(const double sqrt_s) const
1325{
1326 return muZH(sqrt_s)*BrHZZRatio();
1327}
1328
1329const double HiggsChiral::muTHUWHZZ(const double sqrt_s) const
1330{
1331 return muWH(sqrt_s)*BrHZZRatio();
1332}
1333
1334const double HiggsChiral::muTHUVHZZ(const double sqrt_s) const
1335{
1336 return muVH(sqrt_s)*BrHZZRatio();
1337}
1338
1339const double HiggsChiral::muTHUttHZZ(const double sqrt_s) const
1340{
1341 return muttH(sqrt_s)*BrHZZRatio();
1342}
1343
1344const double HiggsChiral::muTHUggHZZ4l(const double sqrt_s) const
1345{
1346 return muggH(sqrt_s)*BrH4lRatio();
1347}
1348
1349const double HiggsChiral::muTHUVBFHZZ4l(const double sqrt_s) const
1350{
1351 return muVBF(sqrt_s)*BrH4lRatio();
1352}
1353
1354const double HiggsChiral::muTHUZHZZ4l(const double sqrt_s) const
1355{
1356 return muZH(sqrt_s)*BrH4lRatio();
1357}
1358
1359const double HiggsChiral::muTHUWHZZ4l(const double sqrt_s) const
1360{
1361 return muWH(sqrt_s)*BrH4lRatio();
1362}
1363
1364const double HiggsChiral::muTHUVHZZ4l(const double sqrt_s) const
1365{
1366 return muVH(sqrt_s)*BrH4lRatio();
1367}
1368
1369const double HiggsChiral::muTHUttHZZ4l(const double sqrt_s) const
1370{
1371 return muttH(sqrt_s)*BrH4lRatio();
1372}
1373
1374const double HiggsChiral::muTHUggHWW(const double sqrt_s) const
1375{
1376 return muggH(sqrt_s)*BrHWWRatio();
1377}
1378
1379const double HiggsChiral::muTHUVBFHWW(const double sqrt_s) const
1380{
1381 return muVBF(sqrt_s)*BrHWWRatio();
1382}
1383
1384const double HiggsChiral::muTHUZHWW(const double sqrt_s) const
1385{
1386 return muZH(sqrt_s)*BrHWWRatio();
1387}
1388
1389const double HiggsChiral::muTHUWHWW(const double sqrt_s) const
1390{
1391 return muWH(sqrt_s)*BrHWWRatio();
1392}
1393
1394const double HiggsChiral::muTHUVHWW(const double sqrt_s) const
1395{
1396 return muVH(sqrt_s)*BrHWWRatio();
1397}
1398
1399const double HiggsChiral::muTHUttHWW(const double sqrt_s) const
1400{
1401 return muttH(sqrt_s)*BrHWWRatio();
1402}
1403
1404const double HiggsChiral::muTHUggHWW2l2v(const double sqrt_s) const
1405{
1406 return muggH(sqrt_s)*BrH2l2vRatio();
1407}
1408
1409const double HiggsChiral::muTHUVBFHWW2l2v(const double sqrt_s) const
1410{
1411 return muVBF(sqrt_s)*BrH2l2vRatio();
1412}
1413
1414const double HiggsChiral::muTHUZHWW2l2v(const double sqrt_s) const
1415{
1416 return muZH(sqrt_s)*BrH2l2vRatio();
1417}
1418
1419const double HiggsChiral::muTHUWHWW2l2v(const double sqrt_s) const
1420{
1421 return muWH(sqrt_s)*BrH2l2vRatio();
1422}
1423
1424const double HiggsChiral::muTHUVHWW2l2v(const double sqrt_s) const
1425{
1426 return muVH(sqrt_s)*BrH2l2vRatio();
1427}
1428
1429const double HiggsChiral::muTHUttHWW2l2v(const double sqrt_s) const
1430{
1431 return muttH(sqrt_s)*BrH2l2vRatio();
1432}
1433
1434const double HiggsChiral::muTHUggHmumu(const double sqrt_s) const
1435{
1436 return muggH(sqrt_s)*BrHmumuRatio();
1437}
1438
1439const double HiggsChiral::muTHUVBFHmumu(const double sqrt_s) const
1440{
1441 return muVBF(sqrt_s)*BrHmumuRatio();
1442}
1443
1444const double HiggsChiral::muTHUZHmumu(const double sqrt_s) const
1445{
1446 return muZH(sqrt_s)*BrHmumuRatio();
1447}
1448
1449const double HiggsChiral::muTHUWHmumu(const double sqrt_s) const
1450{
1451 return muWH(sqrt_s)*BrHmumuRatio();
1452}
1453
1454const double HiggsChiral::muTHUVHmumu(const double sqrt_s) const
1455{
1456 return muVH(sqrt_s)*BrHmumuRatio();
1457}
1458
1459const double HiggsChiral::muTHUttHmumu(const double sqrt_s) const
1460{
1461 return muttH(sqrt_s)*BrHmumuRatio();
1462}
1463
1464const double HiggsChiral::muTHUggHtautau(const double sqrt_s) const
1465{
1466 return muggH(sqrt_s)*BrHtautauRatio();
1467}
1468
1469const double HiggsChiral::muTHUVBFHtautau(const double sqrt_s) const
1470{
1471 return muVBF(sqrt_s)*BrHtautauRatio();
1472}
1473
1474const double HiggsChiral::muTHUZHtautau(const double sqrt_s) const
1475{
1476 return muZH(sqrt_s)*BrHtautauRatio();
1477}
1478
1479const double HiggsChiral::muTHUWHtautau(const double sqrt_s) const
1480{
1481 return muWH(sqrt_s)*BrHtautauRatio();
1482}
1483
1484const double HiggsChiral::muTHUVHtautau(const double sqrt_s) const
1485{
1486 return muVH(sqrt_s)*BrHtautauRatio();
1487}
1488
1489const double HiggsChiral::muTHUttHtautau(const double sqrt_s) const
1490{
1491 return muttH(sqrt_s)*BrHtautauRatio();
1492}
1493
1494const double HiggsChiral::muTHUggHbb(const double sqrt_s) const
1495{
1496 return muggH(sqrt_s)*BrHbbRatio();
1497}
1498
1499const double HiggsChiral::muTHUVBFHbb(const double sqrt_s) const
1500{
1501 return muVBF(sqrt_s)*BrHbbRatio();
1502}
1503
1504const double HiggsChiral::muTHUZHbb(const double sqrt_s) const
1505{
1506 return muZH(sqrt_s)*BrHbbRatio();
1507}
1508
1509const double HiggsChiral::muTHUWHbb(const double sqrt_s) const
1510{
1511 return muWH(sqrt_s)*BrHbbRatio();
1512}
1513
1514const double HiggsChiral::muTHUVHbb(const double sqrt_s) const
1515{
1516 return muVH(sqrt_s)*BrHbbRatio();
1517}
1518
1519const double HiggsChiral::muTHUttHbb(const double sqrt_s) const
1520{
1521 return muttH(sqrt_s)*BrHbbRatio();
1522}
1523
1524const double HiggsChiral::muTHUVBFBRinv(const double sqrt_s) const
1525{
1526 return muVBF(sqrt_s)*Br_H_inv();
1527}
1528
1529const double HiggsChiral::muTHUVBFHinv(const double sqrt_s) const
1530{
1531 return muVBF(sqrt_s)*BrHtoinvRatio();
1532}
1533
1534const double HiggsChiral::muTHUVHBRinv(const double sqrt_s) const
1535{
1536 return muVH(sqrt_s)*Br_H_inv();
1537}
1538
1539const double HiggsChiral::muTHUVHinv(const double sqrt_s) const
1540{
1541 return muVH(sqrt_s)*BrHtoinvRatio();
1542}
1543
1544const double HiggsChiral::muTHUggHZZ4mu(const double sqrt_s) const
1545{
1546 return muggH(sqrt_s)*BrH4muRatio();
1547}
1548
1549const double HiggsChiral::muTHUggHZgamumu(const double sqrt_s) const
1550{
1551 return muggH(sqrt_s)*BrHZgamumuRatio();
1552}
1553
1555//-----------------------------------------------------------------------------------------
1556//-- Decays ---
1557//-----------------------------------------------------------------------------------------
1559
1560
1561const double HiggsChiral::computeGammaTotalRatio() const
1562{
1563
1564 return ((computecg() * computecg() * trueSM.computeBrHtogg()
1565 + computecV() * computecV() * trueSM.computeBrHtoWW()
1566 + computecV() * computecV() * trueSM.computeBrHtoZZ()
1567 + computecZga() * computecZga() * trueSM.computeBrHtoZga()
1568 + computecgaga() * computecgaga() * trueSM.computeBrHtogaga()
1569 + computecmu() * computecmu() * trueSM.computeBrHtomumu()
1570 + computectau() * computectau() * trueSM.computeBrHtotautau()
1571 + computecc() * computecc() * trueSM.computeBrHtocc()
1572 + computecs() * computecs() * trueSM.computeBrHtoss()
1573 + computecb() * computecb() * trueSM.computeBrHtobb())
1574 / (trueSM.computeBrHtogg()
1575 + trueSM.computeBrHtoWW()
1576 + trueSM.computeBrHtoZZ()
1577 + trueSM.computeBrHtoZga()
1578 + trueSM.computeBrHtogaga()
1579 + trueSM.computeBrHtomumu()
1580 + trueSM.computeBrHtotautau()
1581 + trueSM.computeBrHtocc()
1582 + trueSM.computeBrHtoss()
1583 + trueSM.computeBrHtobb()));
1584
1585// return (GammaTotal() / trueSM.computeGammaHTotal());
1586}
1587
1589
1590const double HiggsChiral::computecg() const
1591{
1592 if (loopComputed) return cg_loop;
1593
1594 double Mt=trueSM.getQuarks(QCD::TOP).getMass();
1595 double Mb=trueSM.getQuarks(QCD::BOTTOM).getMass();
1596 double Mc=trueSM.getQuarks(QCD::CHARM).getMass();
1597 double Ms=trueSM.getQuarks(QCD::STRANGE).getMass();
1598 double TAUt=4.0*Mt*Mt/(mHl*mHl);
1599 double TAUb=4.0*Mb*Mb/(mHl*mHl);
1600 double TAUc=4.0*Mc*Mc/(mHl*mHl);
1601 double TAUs=4.0*Ms*Ms/(mHl*mHl);
1602 double cgEff;
1603
1604 cgEff = ( ( ct * 0.5*TAUt*(1.0+(1.0-TAUt)*f_func(TAUt))
1605 +cb * 0.5*TAUb*(1.0+(1.0-TAUb)*f_func(TAUb))
1606 +cc * 0.5*TAUc*(1.0+(1.0-TAUc)*f_func(TAUc))
1607 +cs * 0.5*TAUs*(1.0+(1.0-TAUs)*f_func(TAUs)) ) * (1.0+11.0*AlsMz/(4.0*M_PI))
1608 +cg * 0.5 ).abs2();
1609
1610 cgEff = cgEff / ( ( 0.5*TAUt*(1.0+(1.0-TAUt)*f_func(TAUt))
1611 + 0.5*TAUb*(1.0+(1.0-TAUb)*f_func(TAUb))
1612 + 0.5*TAUc*(1.0+(1.0-TAUc)*f_func(TAUc))
1613 + 0.5*TAUs*(1.0+(1.0-TAUs)*f_func(TAUs)) ) * (1.0+11.0*AlsMz/(4.0*M_PI))).abs2();
1614
1615 return (sqrt(cgEff));
1616}
1617
1618const double HiggsChiral::computecV() const
1619{
1620 return cv;
1621}
1622
1623const double HiggsChiral::computecZga() const
1624{
1625 if (loopComputed) return cZga_loop;
1626
1627 double Mt=trueSM.getQuarks(QCD::TOP).getMass();
1628 double Mb=trueSM.getQuarks(QCD::BOTTOM).getMass();
1629 double Mc=trueSM.getQuarks(QCD::CHARM).getMass();
1630 double Ms=trueSM.getQuarks(QCD::STRANGE).getMass();
1631 double Mtau=trueSM.getLeptons(StandardModel::TAU).getMass();
1632 double Mmu=trueSM.getLeptons(StandardModel::MU).getMass();
1633 double MW=trueSM.Mw_tree();
1634 double cW2=trueSM.c02();
1635 double sW2=1.0-cW2;
1636 double TAUt=4.0*Mt*Mt/(mHl*mHl);
1637 double TAUb=4.0*Mb*Mb/(mHl*mHl);
1638 double TAUc=4.0*Mc*Mc/(mHl*mHl);
1639 double TAUs=4.0*Ms*Ms/(mHl*mHl);
1640 double TAUtau=4.0*Mtau*Mtau/(mHl*mHl);
1641 double TAUmu=4.0*Mmu*Mmu/(mHl*mHl);
1642 double TAUw=4.0*MW*MW/(mHl*mHl);
1643 double LAMt=4.0*Mt*Mt/(Mz*Mz);
1644 double LAMb=4.0*Mb*Mb/(Mz*Mz);
1645 double LAMc=4.0*Mc*Mc/(Mz*Mz);
1646 double LAMs=4.0*Ms*Ms/(Mz*Mz);
1647 double LAMtau=4.0*Mtau*Mtau/(Mz*Mz);
1648 double LAMmu=4.0*Mmu*Mmu/(Mz*Mz);
1649 double LAMw=4.0*MW*MW/(Mz*Mz);
1650 double cZgaEff;
1651
1652 cZgaEff = ((-ct * 4.0*(0.5-4.0/3.0*sW2)*(Int1(TAUt,LAMt)-Int2(TAUt,LAMt)) * (1.0-AlsMz/M_PI)
1653 +cb * 2.0*(-0.5+2.0/3.0*sW2)*(Int1(TAUb,LAMb)-Int2(TAUb,LAMb))
1654 -cc * 4.0*(0.5-4.0/3.0*sW2)*(Int1(TAUc,LAMc)-Int2(TAUc,LAMc))
1655 +cs * 2.0*(-0.5+2.0/3.0*sW2)*(Int1(TAUs,LAMs)-Int2(TAUs,LAMs))
1656 +ctau * 2.0*(-0.5+2.0*sW2)*(Int1(TAUtau,LAMtau)-Int2(TAUtau,LAMtau))
1657 +cmu * 2.0*(-0.5+2.0*sW2)*(Int1(TAUmu,LAMmu)-Int2(TAUmu,LAMmu)) )/sqrt(sW2*cW2)
1658 -cv * sqrt(cW2/sW2)*(4.0*(3.0-sW2/cW2)*Int2(TAUw,LAMw)
1659 +((1.0+2.0/TAUw)*sW2/cW2-(5.0+2.0/TAUw))*Int1(TAUw,LAMw))
1660 +cZga).abs2();
1661
1662 cZgaEff = cZgaEff / ((-4.0*(0.5-4.0/3.0*sW2)*(Int1(TAUt,LAMt)-Int2(TAUt,LAMt)) * (1.0-AlsMz/M_PI)
1663 + 2.0*(-0.5+2.0/3.0*sW2)*(Int1(TAUb,LAMb)-Int2(TAUb,LAMb))
1664 - 4.0*(0.5-4.0/3.0*sW2)*(Int1(TAUc,LAMc)-Int2(TAUc,LAMc))
1665 + 2.0*(-0.5+2.0/3.0*sW2)*(Int1(TAUs,LAMs)-Int2(TAUs,LAMs))
1666 + 2.0*(-0.5+2.0*sW2)*(Int1(TAUtau,LAMtau)-Int2(TAUtau,LAMtau))
1667 + 2.0*(-0.5+2.0*sW2)*(Int1(TAUmu,LAMmu)-Int2(TAUmu,LAMmu)) )/sqrt(sW2*cW2)
1668 - sqrt(cW2/sW2)*(4.0*(3.0-sW2/cW2)*Int2(TAUw,LAMw)
1669 +((1.0+2.0/TAUw)*sW2/cW2-(5.0+2.0/TAUw))*Int1(TAUw,LAMw))).abs2();
1670
1671 return (sqrt(cZgaEff));
1672}
1673
1674const double HiggsChiral::computecgaga() const
1675{
1676 if (loopComputed) return cga_loop;
1677
1678 double Mt=trueSM.getQuarks(QCD::TOP).getMass();
1679 double Mb=trueSM.getQuarks(QCD::BOTTOM).getMass();
1680 double Mc=trueSM.getQuarks(QCD::CHARM).getMass();
1681 double Ms=trueSM.getQuarks(QCD::STRANGE).getMass();
1682 double Mtau=trueSM.getLeptons(StandardModel::TAU).getMass();
1683 double Mmu=trueSM.getLeptons(StandardModel::MU).getMass();
1684 double MW=trueSM.Mw_tree();
1685 double TAUt=4.0*Mt*Mt/(mHl*mHl);
1686 double TAUb=4.0*Mb*Mb/(mHl*mHl);
1687 double TAUc=4.0*Mc*Mc/(mHl*mHl);
1688 double TAUs=4.0*Ms*Ms/(mHl*mHl);
1689 double TAUtau=4.0*Mtau*Mtau/(mHl*mHl);
1690 double TAUmu=4.0*Mmu*Mmu/(mHl*mHl);
1691 double TAUw=4.0*MW*MW/(mHl*mHl);
1692 double cgagaEff;
1693
1694 cgagaEff = ( ct * (8./3.)*TAUt*(1.+(1.-TAUt)*f_func(TAUt)) * (1.0-AlsMz/M_PI)
1695 +cb * (2./3.)*TAUb*(1.+(1.-TAUb)*f_func(TAUb))
1696 +cc * (8./3.)*TAUc*(1.+(1.-TAUc)*f_func(TAUc))
1697 +cs * (2./3.)*TAUs*(1.+(1.-TAUs)*f_func(TAUs))
1698 +ctau * 2.0*TAUtau*(1.+(1.-TAUtau)*f_func(TAUtau))
1699 +cmu * 2.0*TAUmu*(1.+(1.-TAUmu)*f_func(TAUmu))
1700 -cv * (2.0+3.0*TAUw+3.0*TAUw*(2.0-TAUw)*f_func(TAUw))
1701 +cga * 2.0 ).abs2();
1702
1703 cgagaEff = cgagaEff / ( (8./3.)*TAUt*(1.+(1.-TAUt)*f_func(TAUt)) * (1.0-AlsMz/M_PI)
1704 + (2./3.)*TAUb*(1.+(1.-TAUb)*f_func(TAUb))
1705 + (8./3.)*TAUc*(1.+(1.-TAUc)*f_func(TAUc))
1706 + (2./3.)*TAUs*(1.+(1.-TAUs)*f_func(TAUs))
1707 + 2.0*TAUtau*(1.+(1.-TAUtau)*f_func(TAUtau))
1708 + 2.0*TAUmu*(1.+(1.-TAUmu)*f_func(TAUmu))
1709 - (2.0+3.0*TAUw+3.0*TAUw*(2.0-TAUw)*f_func(TAUw)) ).abs2();
1710
1711 return (sqrt(cgagaEff));
1712}
1713
1714const double HiggsChiral::computecmu() const
1715{
1716 return cmu;
1717}
1718
1719const double HiggsChiral::computectau() const
1720{
1721 return ctau;
1722}
1723
1724const double HiggsChiral::computecc() const
1725{
1726 return cc;
1727}
1728
1729const double HiggsChiral::computecs() const
1730{
1731 return cs;
1732}
1733
1734const double HiggsChiral::computect() const
1735{
1736 return ct;
1737}
1738
1739const double HiggsChiral::computecb() const
1740{
1741 return cb;
1742}
1743
1744
1746
1747const double HiggsChiral::deltaG_hgg() const
1748{
1749 // Not needed/implemented yet
1750 return 0.0;
1751}
1752
1753const double HiggsChiral::deltaG_hggRatio() const{
1754 return (computecg() - 1.0);
1755}
1756
1757const double HiggsChiral::deltaG1_hWW() const{
1758 return 0.0;
1759}
1760
1761const double HiggsChiral::deltaG2_hWW() const{
1762 return 0.0;
1763}
1764
1765const double HiggsChiral::deltaG3_hWW() const{
1766
1767 double gSM = 2.0 * (trueSM.Mw_tree())* (trueSM.Mw_tree()) / (trueSM.v());
1768
1769 return ( gSM * (computecV() - 1.0 ) );
1770}
1771
1772const double HiggsChiral::deltaG1_hZZ() const{
1773 return 0.0;
1774}
1775
1776const double HiggsChiral::deltaG2_hZZ() const{
1777 return 0.0;
1778}
1779
1780const double HiggsChiral::deltaG3_hZZ() const{
1781
1782 double gSM = (trueSM.getMz()) * (trueSM.getMz()) / (trueSM.v());
1783
1784 return ( gSM * (computecV() - 1.0 ) );
1785}
1786
1787const double HiggsChiral::deltaG1_hZA() const{
1788 // Not needed/implemented yet
1789 return 0.0;
1790}
1791
1792const double HiggsChiral::deltaG1_hZARatio() const{
1793 return (computecZga() - 1.0);
1794}
1795
1796const double HiggsChiral::deltaG2_hZA() const{
1797 return 0.0;
1798}
1799
1800const double HiggsChiral::deltaG_hAA() const{
1801 // Not needed/implemented yet
1802 return 0.0;
1803}
1804
1805const double HiggsChiral::deltaG_hAARatio() const{
1806 return (computecgaga() - 1.0);
1807}
1808
1809gslpp::complex HiggsChiral::deltaG_hff(const Particle p) const{
1810
1811 double gSM = -(p.getMass()) / (trueSM.v());
1812
1813 if ( p.is("ELECTRON") ) {
1814 return ( 0.0 );
1815 } else if ( p.is("MU") ) {
1816 return ( gSM * (computecmu() - 1.0 ) );
1817 } else if ( p.is("TAU") ) {
1818 return ( gSM * (computectau() - 1.0 ) );
1819 } else if ( p.is("UP") ) {
1820 return ( 0.0 );
1821 } else if ( p.is("CHARM") ) {
1822 return ( gSM * (computecc() - 1.0 ) );
1823 } else if ( p.is("TOP") ) {
1824 return ( gSM * (computect() - 1.0 ) );
1825 } else if ( p.is("DOWN") ) {
1826 return ( 0.0 );
1827 } else if ( p.is("STRANGE") ) {
1828 return ( gSM * (computecs() - 1.0 ) );
1829 } else if ( p.is("BOTTOM") ) {
1830 return ( gSM * (computecb() - 1.0 ) );
1831 } else {
1832 return 0.0;
1833 }
1834}
1835
1837
1838const double HiggsChiral::kappamueff() const
1839{
1840 return computecmu();
1841}
1842
1843const double HiggsChiral::kappataueff() const
1844{
1845 return computectau();
1846}
1847
1848const double HiggsChiral::kappaceff() const
1849{
1850 return computecc();
1851}
1852
1853const double HiggsChiral::kappaseff() const
1854{
1855 return computecs();
1856}
1857
1858const double HiggsChiral::kappabeff() const
1859{
1860 return computecb();
1861}
1862
1863const double HiggsChiral::kappaGeff() const
1864{
1865 return computecg();
1866}
1867
1868const double HiggsChiral::kappaZeff() const
1869{
1870 return computecV();
1871}
1872
1873const double HiggsChiral::kappaWeff() const
1874{
1875 return computecV();
1876}
1877
1878const double HiggsChiral::kappaAeff() const
1879{
1880 return computecgaga();
1881}
1882
1883const double HiggsChiral::kappaZAeff() const
1884{
1885 return computecZga();
1886}
1887
1889
1890gslpp::complex HiggsChiral::f_func(const double x) const{
1891 if(x<1) {
1892 gslpp::complex z = -gslpp::complex::i()*M_PI;
1893 return -pow(log((1+sqrt(1-x))/(1-sqrt(1-x)))+z,2)/4.0;
1894 }
1895 else {
1896 return pow(asin(sqrt(1.0/x)),2);
1897 }
1898}
1899
1900gslpp::complex HiggsChiral::g_func(const double x) const{
1901 if(x<1) {
1902 gslpp::complex z = -gslpp::complex::i()*M_PI;
1903 gslpp::complex gs1 = sqrt(1.0-x)*(log((1.0+sqrt(1.0-x))/(1.0-sqrt(1.0-x)))+z)/2.0;
1904 return gs1;
1905 }
1906 else {
1907 gslpp::complex gg1 = sqrt(x-1.0)*asin(sqrt(1.0/x));
1908 return gg1;
1909 }
1910}
1911
1912gslpp::complex HiggsChiral::Int1(const double tau, const double lambda) const{
1913 return tau*lambda/(tau-lambda)/2.0+tau*tau*lambda*lambda/((tau-lambda)
1914 *(tau-lambda))/2.0*(f_func(tau)-f_func(lambda))+tau*tau*lambda/((tau-lambda)
1915 *(tau-lambda))*(g_func(tau)-g_func(lambda));
1916}
1917
1918gslpp::complex HiggsChiral::Int2(const double tau, const double lambda) const{
1919 return -tau*lambda/(tau-lambda)/2.0*(f_func(tau)-f_func(lambda));
1920}
DPars
std::map< std::string, double > DPars
Definition: Minimal.cpp:11
HiggsChiral::kappataueff
virtual const double kappataueff() const
The effective coupling .
Definition: HiggsChiral.cpp:1843
HiggsChiral::muTHUWHZZ
virtual const 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:1329
HiggsChiral::BrHmumuRatio
virtual const double BrHmumuRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:668
HiggsChiral::deltaG2_hZA
virtual const double deltaG2_hZA() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1796
HiggsChiral::f_func
gslpp::complex f_func(const double x) const
Definition: HiggsChiral.cpp:1890
HiggsChiral::cgminuscga
virtual const double cgminuscga() const
The value of .
Definition: HiggsChiral.cpp:451
HiggsChiral::BrHssRatio
virtual const double BrHssRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:689
HiggsChiral::muTHUWHgaga
virtual const 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:1269
HiggsChiral::muTHUWHtautau
virtual const 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:1479
HiggsChiral::muggHbb
virtual const 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:1201
HiggsChiral::computecZga
virtual const double computecZga() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1623
HiggsChiral::BrH2L2vRatio
virtual const double BrH2L2vRatio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:720
HiggsChiral::computecg
virtual const double computecg() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1590
HiggsChiral::BrH4eRatio
virtual const double BrH4eRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:785
HiggsChiral::muWH
virtual const 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:216
HiggsChiral::muTHUggHWW
virtual const 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:1374
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:1809
HiggsChiral::muTHUttHWW2l2v
virtual const 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:1429
HiggsChiral::kappabeff
virtual const double kappabeff() const
The effective coupling .
Definition: HiggsChiral.cpp:1858
HiggsChiral::computectau
virtual const double computectau() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1719
HiggsChiral::ccminusctau
virtual const double ccminusctau() const
The value of .
Definition: HiggsChiral.cpp:476
HiggsChiral::muppHmumu
virtual const 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:1231
HiggsChiral::muTHUZHbb
virtual const 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:1504
HiggsChiral::HChiralvars
static const std::string HChiralvars[NHChiralvars]
A string array containing the labels of the model parameters in HiggsKvKf.
Definition: HiggsChiral.h:117
HiggsChiral::expZgaLimitCMS
double expZgaLimitCMS
Definition: HiggsChiral.h:2238
HiggsChiral::obsZgaLimitATLAS13
double obsZgaLimitATLAS13
Definition: HiggsChiral.h:2231
HiggsChiral::kappaceff
virtual const double kappaceff() const
The effective coupling .
Definition: HiggsChiral.cpp:1848
HiggsChiral::muTHUggHtautau
virtual const 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:1464
HiggsChiral::BrH2v2uRatio
virtual const double BrH2v2uRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:765
HiggsChiral::muZH
virtual const 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:221
HiggsChiral::mueeZH
virtual const double mueeZH(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:226
HiggsChiral::muggH
virtual const 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:184
HiggsChiral::muTHUWHZga
virtual const 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:1299
HiggsChiral::muTHUWHZZ4l
virtual const 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:1359
HiggsChiral::muVHbb
virtual const 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:1221
HiggsChiral::kappaseff
virtual const double kappaseff() const
The effective coupling .
Definition: HiggsChiral.cpp:1853
HiggsChiral::BrH2u2dRatio
virtual const double BrH2u2dRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:750
HiggsChiral::muTHUVBFHZZ
virtual const 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:1319
HiggsChiral::kappaWeff
virtual const double kappaWeff() const
The effective coupling .
Definition: HiggsChiral.cpp:1873
HiggsChiral::muVBFHWW2l2v
virtual const 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:1116
HiggsChiral::mueettH
virtual const double mueettH(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:279
HiggsChiral::BrHZgaRatio
virtual const double BrHZgaRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:633
HiggsChiral::muttHZZ4l
virtual const 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:1076
HiggsChiral::muTHUggHZZ4mu
virtual const 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:1544
HiggsChiral::muVBFgamma
virtual const 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:196
HiggsChiral::expZgaLimitATLAS13
double expZgaLimitATLAS13
Definition: HiggsChiral.h:2235
HiggsChiral::muTHUVHZZ4l
virtual const 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:1364
HiggsChiral::muZHgaga
virtual const 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:971
HiggsChiral::BrH2v2vRatio
virtual const double BrH2v2vRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:715
HiggsChiral::muZHWW2l2v
virtual const 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:1121
HiggsChiral::muTHUZHmumu
virtual const 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:1444
HiggsChiral::muTHUttHbb
virtual const 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:1519
HiggsChiral::obsZgaLimitCMS13
double obsZgaLimitCMS13
Definition: HiggsChiral.h:2232
HiggsChiral::muTHUZHZga
virtual const 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:1294
HiggsChiral::HiggsChiral
HiggsChiral()
The default constructor.
Definition: HiggsChiral.cpp:14
HiggsChiral::muggHZZ4l
virtual const 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:1051
HiggsChiral::muTHUttHZZ4l
virtual const 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:1369
HiggsChiral::g_func
gslpp::complex g_func(const double x) const
Definition: HiggsChiral.cpp:1900
HiggsChiral::muTHUVHZga
virtual const 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:1304
HiggsChiral::expZgaLimitATLAS
double expZgaLimitATLAS
Definition: HiggsChiral.h:2237
HiggsChiral::BrH4lRatio
virtual const double BrH4lRatio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:891
HiggsChiral::BrH2L2dRatio
virtual const double BrH2L2dRatio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:760
HiggsChiral::muVBFHZga
virtual const 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:996
HiggsChiral::BrHccRatio
virtual const double BrHccRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:682
HiggsChiral::muVHmumu
virtual const 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:1161
HiggsChiral::BrH4vRatio
virtual const double BrH4vRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:795
HiggsChiral::BrH2L2uRatio
virtual const double BrH2L2uRatio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:755
HiggsChiral::muWHbb
virtual const 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:1216
HiggsChiral::getobsZgaLimitCMS
double getobsZgaLimitCMS() const
a getter for the value of the observed upper limit in from CMS
Definition: HiggsChiral.h:268
HiggsChiral::computecV
virtual const double computecV() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1618
HiggsChiral::muTHUttHtautau
virtual const 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:1489
HiggsChiral::muZHmumu
virtual const 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:1151
HiggsChiral::muVBFHWW
virtual const 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:1086
HiggsChiral::BrH2u2uRatio
virtual const double BrH2u2uRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:740
HiggsChiral::BrH2d2dRatio
virtual const double BrH2d2dRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:745
HiggsChiral::muTHUttHgaga
virtual const 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:1279
HiggsChiral::BrHLvvLRatio
virtual const double BrHLvvLRatio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:810
HiggsChiral::UpperLimitZgammaA
virtual const double UpperLimitZgammaA(const double sqrt_s) const
Observable implementing the contribution to the likelihood from the upper limit in from ATLAS.
Definition: HiggsChiral.cpp:393
HiggsChiral::muTHUggHZga
virtual const 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:1284
HiggsChiral::BrHLvudRatio
virtual const double BrHLvudRatio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:825
HiggsChiral::obliqueT
virtual const double obliqueT() const
The oblique parameter .
Definition: HiggsChiral.cpp:165
HiggsChiral::muWHmumu
virtual const 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:1156
HiggsChiral::muZHWW
virtual const 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:1091
HiggsChiral::kappaGeff
virtual const double kappaGeff() const
The effective coupling .
Definition: HiggsChiral.cpp:1863
HiggsChiral::muTHUZHWW
virtual const 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:1384
HiggsChiral::muTHUVBFHWW2l2v
virtual const 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:1409
HiggsChiral::BrH2l2vRatio
virtual const double BrH2l2vRatio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:896
HiggsChiral::cVplusctau
virtual const double cVplusctau() const
The value of .
Definition: HiggsChiral.cpp:461
HiggsChiral::muTHUVHtautau
virtual const 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:1484
HiggsChiral::cga_loop
double cga_loop
Definition: HiggsChiral.h:2241
HiggsChiral::muTHUVHinv
virtual const 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:1539
HiggsChiral::Int1
gslpp::complex Int1(const double tau, const double lambda) const
Definition: HiggsChiral.cpp:1912
HiggsChiral::muTHUVBFHmumu
virtual const 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:1439
HiggsChiral::muTHUVBFHZZ4l
virtual const 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:1349
HiggsChiral::muTHUggHWW2l2v
virtual const 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:1404
HiggsChiral::Gammagg
virtual const double Gammagg() const
The decay width in the current model.
Definition: HiggsChiral.cpp:481
HiggsChiral::deltaG_hAA
virtual const double deltaG_hAA() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1800
HiggsChiral::BrH2e2muRatio
virtual const double BrH2e2muRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:710
HiggsChiral::muttHmumu
virtual const 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:1166
HiggsChiral::mummHNWA
virtual const double mummHNWA(const double sqrt_s) const
The ratio between the production cross-section in the current model and in the Standard Model,...
Definition: HiggsChiral.cpp:302
HiggsChiral::mummHvv
virtual const double mummHvv(const double sqrt_s) const
The ratio between the production cross-section in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:319
HiggsChiral::deltaG3_hWW
virtual const double deltaG3_hWW() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1765
HiggsChiral::deltaG2_hWW
virtual const double deltaG2_hWW() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1761
HiggsChiral::muggHWW
virtual const 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:1081
HiggsChiral::muWHWW
virtual const 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:1096
HiggsChiral::kappaZAeff
virtual const double kappaZAeff() const
The effective coupling .
Definition: HiggsChiral.cpp:1883
HiggsChiral::muVHZga
virtual const 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:1011
HiggsChiral::muttHZga
virtual const 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:1016
HiggsChiral::muTHUVHgaga
virtual const 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:1274
HiggsChiral::muTHUZHZZ
virtual const 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:1324
HiggsChiral::BrHZgamumuRatio
virtual const double BrHZgamumuRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:654
HiggsChiral::deltaG1_hZARatio
virtual const double deltaG1_hZARatio() const
The full new physics contribution to the coupling of the effective interaction , including new local ...
Definition: HiggsChiral.cpp:1792
HiggsChiral::muTHUggHZZ
virtual const 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:1314
HiggsChiral::muTHUVHmumu
virtual const 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:1454
HiggsChiral::deltaG1_hZA
virtual const double deltaG1_hZA() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1787
HiggsChiral::muTHUVBFHgaga
virtual const 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:1259
HiggsChiral::muWHZZ
virtual const 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:1036
HiggsChiral::computecgaga
virtual const double computecgaga() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1674
HiggsChiral::GammaZga
virtual const double GammaZga() const
The decay width in the current model.
Definition: HiggsChiral.cpp:513
HiggsChiral::cg_loop
double cg_loop
Definition: HiggsChiral.h:2240
HiggsChiral::GammaTotal
virtual const double GammaTotal() const
The total decay width of the Higgs boson in the current model.
Definition: HiggsChiral.cpp:599
HiggsChiral::expZgaLimitCMS13
double expZgaLimitCMS13
Definition: HiggsChiral.h:2236
HiggsChiral::muVHgaga
virtual const 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:981
HiggsChiral::cVpluscb
virtual const double cVpluscb() const
The value of .
Definition: HiggsChiral.cpp:456
HiggsChiral::muttH
virtual const 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:251
HiggsChiral::Int2
gslpp::complex Int2(const double tau, const double lambda) const
Definition: HiggsChiral.cpp:1918
HiggsChiral::muTHUWHWW2l2v
virtual const 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:1419
HiggsChiral::deltaG1_hWW
virtual const double deltaG1_hWW() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1757
HiggsChiral::BrHggRatio
virtual const double BrHggRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:605
HiggsChiral::BrH2udRatio
virtual const double BrH2udRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:830
HiggsChiral::muTHUVBFHbb
virtual const 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:1499
HiggsChiral::muTHUWHWW
virtual const 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:1389
HiggsChiral::muTHUZHZZ4l
virtual const 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:1354
HiggsChiral::muTHUttHZZ
virtual const 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:1339
HiggsChiral::muVHWW
virtual const 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:1101
HiggsChiral::BrHZgallRatio
virtual const double BrHZgallRatio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:640
HiggsChiral::muTHUVHWW2l2v
virtual const 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:1424
HiggsChiral::muTHUggHbb
virtual const 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:1494
HiggsChiral::kappamueff
virtual const double kappamueff() const
The effective coupling .
Definition: HiggsChiral.cpp:1838
HiggsChiral::kappaAeff
virtual const double kappaAeff() const
The effective coupling .
Definition: HiggsChiral.cpp:1878
HiggsChiral::muTHUggHgaga
virtual const 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:1254
HiggsChiral::BrH2evRatio
virtual const double BrH2evRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:848
HiggsChiral::muTHUVBFHWW
virtual const 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:1379
HiggsChiral::getexpZgaLimitATLAS
double getexpZgaLimitATLAS() const
a getter for the experimental value of the expected upper limit in from ATLAS
Definition: HiggsChiral.h:295
HiggsChiral::muTHUVBFHtautau
virtual const 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:1469
HiggsChiral::muTHUttHZga
virtual const 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:1309
HiggsChiral::computeGammaTotalRatio
virtual const double computeGammaTotalRatio() const
The ratio of the in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:1561
HiggsChiral::BrHbbRatio
virtual const double BrHbbRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:696
HiggsChiral::BrH4muRatio
virtual const double BrH4muRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:790
HiggsChiral::mueeWBF
virtual const double mueeWBF(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:201
HiggsChiral::mummHmm
virtual const double mummHmm(const double sqrt_s) const
The ratio between the production cross-section in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:329
HiggsChiral::deltaG3_hZZ
virtual const double deltaG3_hZZ() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1780
HiggsChiral::muggHZga
virtual const 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:991
HiggsChiral::BrHZgaeeRatio
virtual const double BrHZgaeeRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:647
HiggsChiral::UpperLimitZgammaA13
virtual const 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:345
HiggsChiral::muttHZZ
virtual const 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:1046
HiggsChiral::getexpZgaLimitCMS
double getexpZgaLimitCMS() const
a getter for the experimental value of the expected upper limit in from CMS
Definition: HiggsChiral.h:304
HiggsChiral::muTHUggHZZ4l
virtual const 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:1344
HiggsChiral::getobsZgaLimitCMS13
double getobsZgaLimitCMS13() const
a getter for the value of the observed upper limit in from CMS at 13 TeV
Definition: HiggsChiral.h:250
HiggsChiral::obliqueS
virtual const double obliqueS() const
The oblique parameter .
Definition: HiggsChiral.cpp:154
HiggsChiral::BrH2mu2vRatio
virtual const double BrH2mu2vRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:735
HiggsChiral::muZHtautau
virtual const 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:1181
HiggsChiral::muTHUVHBRinv
virtual const 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:1534
HiggsChiral::muZHbb
virtual const 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:1211
HiggsChiral::obsZgaLimitATLAS
double obsZgaLimitATLAS
Definition: HiggsChiral.h:2233
HiggsChiral::Gammabb
virtual const double Gammabb() const
The decay width in the current model.
Definition: HiggsChiral.cpp:594
HiggsChiral::muVBFHZZ
virtual const 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:1026
HiggsChiral::getobsZgaLimitATLAS
double getobsZgaLimitATLAS() const
a getter for the value of the observed upper limit in from ATLAS
Definition: HiggsChiral.h:259
HiggsChiral::mummZH
virtual const double mummZH(const double sqrt_s) const
The ratio between the production cross-section in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:311
HiggsChiral::BrHevmuvRatio
virtual const double BrHevmuvRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:815
HiggsChiral::muttHWW2l2v
virtual const 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:1136
HiggsChiral::computect
virtual const double computect() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1734
HiggsChiral::mueeZllH
virtual const double mueeZllH(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:231
HiggsChiral::muVHWW2l2v
virtual const 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:1131
HiggsChiral::cgplusct
virtual const double cgplusct() const
The value of .
Definition: HiggsChiral.cpp:441
HiggsChiral::muTHUttHmumu
virtual const 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:1459
HiggsChiral::BrH2L2v2Ratio
virtual const double BrH2L2v2Ratio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:725
HiggsChiral::muVBFHtautau
virtual const 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:1176
HiggsChiral::muTHUggHmumu
virtual const 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:1434
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:125
HiggsChiral::muTHUVHZZ
virtual const 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:1334
HiggsChiral::kappaZeff
virtual const double kappaZeff() const
The effective coupling .
Definition: HiggsChiral.cpp:1868
HiggsChiral::BrH4uRatio
virtual const double BrH4uRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:800
HiggsChiral::muTHUVBFHinv
virtual const 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:1529
HiggsChiral::getobsZgaLimitATLAS13
double getobsZgaLimitATLAS13() const
a getter for the value of the observed upper limit in from ATLAS at 13 TeV
Definition: HiggsChiral.h:241
HiggsChiral::BrH4dRatio
virtual const double BrH4dRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:805
HiggsChiral::PostUpdate
virtual bool PostUpdate()
The post-update method for HiggsChiral.
Definition: HiggsChiral.cpp:45
HiggsChiral::muWHZZ4l
virtual const 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:1066
HiggsChiral::BrH4fRatio
virtual const double BrH4fRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:860
HiggsChiral::BrH2muvRatio
virtual const double BrH2muvRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:854
HiggsChiral::BrHZZRatio
virtual const double BrHZZRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:619
HiggsChiral::muZHZga
virtual const 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:1001
HiggsChiral::mummH
virtual const double mummH(const double sqrt_s) const
The ratio between the production cross-section in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:285
HiggsChiral::computecmu
virtual const double computecmu() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1714
HiggsChiral::muTHUttHWW
virtual const 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:1399
HiggsChiral::BrHWWRatio
virtual const double BrHWWRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:612
HiggsChiral::muVBFHZZ4l
virtual const 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:1056
HiggsChiral::Gammatautau
virtual const double Gammatautau() const
The decay width in the current model.
Definition: HiggsChiral.cpp:579
HiggsChiral::cgaplusct
virtual const double cgaplusct() const
The value of .
Definition: HiggsChiral.cpp:446
HiggsChiral::UpperLimitZgammaC13
virtual const 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:369
HiggsChiral::muWHtautau
virtual const 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:1186
HiggsChiral::setFlag
virtual bool setFlag(const std::string name, const bool value)
A method to set a flag of HiggsChiral.
Definition: HiggsChiral.cpp:137
HiggsChiral::muggHpttH
virtual const 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:268
HiggsChiral::Gammass
virtual const double Gammass() const
The decay width in the current model.
Definition: HiggsChiral.cpp:589
HiggsChiral::BrH4L2Ratio
virtual const double BrH4L2Ratio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:780
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:277
HiggsChiral::deltaG2_hZZ
virtual const double deltaG2_hZZ() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1776
HiggsChiral::BrHgagaRatio
virtual const double BrHgagaRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:661
HiggsChiral::muWHZga
virtual const 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:1006
HiggsChiral::BrHtoinvRatio
virtual const double BrHtoinvRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:902
HiggsChiral::mueeZqqH
virtual const double mueeZqqH(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:236
HiggsChiral::BrHVVRatio
virtual const double BrHVVRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:626
HiggsChiral::muTHUVBFHZga
virtual const 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:1289
HiggsChiral::muTHUggHZgamumu
virtual const 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:1549
HiggsChiral::deltaG_hggRatio
virtual const double deltaG_hggRatio() const
The full new physics contribution to the coupling of the effective interaction , including new local ...
Definition: HiggsChiral.cpp:1753
HiggsChiral::muWHgaga
virtual const 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:976
HiggsChiral::mutHq
virtual const 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:256
HiggsChiral::computecb
virtual const double computecb() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1739
HiggsChiral::muZHZZ4l
virtual const 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:1061
HiggsChiral::muTHUVBFBRinv
virtual const 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:1524
HiggsChiral::cZga_loop
double cZga_loop
Definition: HiggsChiral.h:2242
HiggsChiral::muggHmumu
virtual const 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:1141
HiggsChiral::mueeZBF
virtual const double mueeZBF(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:211
HiggsChiral::cbminuscc
virtual const double cbminuscc() const
The value of .
Definition: HiggsChiral.cpp:466
HiggsChiral::muttHtautau
virtual const 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:1196
HiggsChiral::muVBFHbb
virtual const 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:1206
HiggsChiral::muTHUVHWW
virtual const 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:1394
HiggsChiral::muWHWW2l2v
virtual const 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:1126
HiggsChiral::muVHtautau
virtual const 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:1191
HiggsChiral::loopComputed
bool loopComputed
Definition: HiggsChiral.h:2243
HiggsChiral::muVH
virtual const 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:241
HiggsChiral::Gammacc
virtual const double Gammacc() const
The decay width in the current model.
Definition: HiggsChiral.cpp:584
HiggsChiral::muggHtautau
virtual const 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:1171
HiggsChiral::computecs
virtual const double computecs() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1729
HiggsChiral::muTHUZHWW2l2v
virtual const 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:1414
HiggsChiral::BrH2LvRatio
virtual const double BrH2LvRatio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:836
HiggsChiral::muttHgaga
virtual const 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:986
HiggsChiral::deltaG1_hZZ
virtual const double deltaG1_hZZ() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1772
HiggsChiral::FlagUniversalcf
bool FlagUniversalcf
A boolean flag that is true if all cf take the same universal value.
Definition: HiggsChiral.h:2245
HiggsChiral::deltaG_hgg
virtual const double deltaG_hgg() const
The new physics contribution to the coupling of the effective interaction .
Definition: HiggsChiral.cpp:1747
HiggsChiral::BrHudduRatio
virtual const double BrHudduRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:820
HiggsChiral::Gammamumu
virtual const double Gammamumu() const
The decay width in the current model.
Definition: HiggsChiral.cpp:574
HiggsChiral::GammaWW
virtual const double GammaWW() const
The decay width in the current model.
Definition: HiggsChiral.cpp:503
HiggsChiral::muVHZZ4l
virtual const 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:1071
HiggsChiral::FlagUniversalcvcf
bool FlagUniversalcvcf
A boolean flag that is true if all cv and cf take the same universal value.
Definition: HiggsChiral.h:2246
HiggsChiral::muTHUVHbb
virtual const 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:1514
HiggsChiral::computecc
virtual const double computecc() const
A method to compute the ratio of the coupling in the current model and in the SM.
Definition: HiggsChiral.cpp:1724
HiggsChiral::muVBFpVH
virtual const 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:246
HiggsChiral::mummttH
virtual const double mummttH(const double sqrt_s) const
The ratio between the production cross-section in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:337
HiggsChiral::BrHtautauRatio
virtual const double BrHtautauRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:675
HiggsChiral::obliqueU
virtual const double obliqueU() const
The oblique parameter .
Definition: HiggsChiral.cpp:177
HiggsChiral::BrH4LRatio
virtual const double BrH4LRatio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:775
HiggsChiral::muVBFHmumu
virtual const 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:1146
HiggsChiral::BrH2L2LRatio
virtual const double BrH2L2LRatio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:705
HiggsChiral::muVHZZ
virtual const 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:1041
HiggsChiral::muZHZZ
virtual const 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:1031
HiggsChiral::GammaZZ
virtual const double GammaZZ() const
The decay width in the current model.
Definition: HiggsChiral.cpp:508
HiggsChiral::BrH2v2dRatio
virtual const double BrH2v2dRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:770
HiggsChiral::muTHUZHtautau
virtual const 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:1474
HiggsChiral::muggHZZ
virtual const 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:1021
HiggsChiral::muggHgaga
virtual const 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:909
HiggsChiral::muTHUWHmumu
virtual const 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:1449
HiggsChiral::cbminusctau
virtual const double cbminusctau() const
The value of .
Definition: HiggsChiral.cpp:471
HiggsChiral::deltaG_hAARatio
virtual const double deltaG_hAARatio() const
The full new physics contribution to the coupling of the effective interaction , including new local ...
Definition: HiggsChiral.cpp:1805
HiggsChiral::Gammagaga
virtual const double Gammagaga() const
The decay width in the current model.
Definition: HiggsChiral.cpp:548
HiggsChiral::muTHUWHbb
virtual const 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:1509
HiggsChiral::obsZgaLimitCMS
double obsZgaLimitCMS
Definition: HiggsChiral.h:2234
HiggsChiral::muggHWW2l2v
virtual const 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:1111
HiggsChiral::mueeHvv
virtual const double mueeHvv(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:206
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:83
HiggsChiral::muVBF
virtual const 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:191
HiggsChiral::muTHUZHgaga
virtual const 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:1264
HiggsChiral::muggHgagaInt
virtual const 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:914
HiggsChiral::BrH2e2vRatio
virtual const double BrH2e2vRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:730
HiggsChiral::NHChiralvars
static const int NHChiralvars
The number of the model parameters.
Definition: HiggsChiral.h:112
HiggsChiral::muttHWW
virtual const 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:1106
HiggsChiral::UpperLimitZgammaC
virtual const double UpperLimitZgammaC(const double sqrt_s) const
Observable implementing the contribution to the likelihood from the upper limit in from CMS.
Definition: HiggsChiral.cpp:417
HiggsChiral::muVBFHgaga
virtual const 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:966
HiggsChiral::muttHbb
virtual const 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:1226
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:286
HiggsChiral::BrH2Lv2Ratio
virtual const double BrH2Lv2Ratio() const
The ratio of the Br ( ) in the current model and in the Standard Model.
Definition: HiggsChiral.cpp:842
Model::addMissingModelParameter
void addMissingModelParameter(const std::string &missingParameterName)
Definition: Model.h:250
Model::ModelParamMap
std::map< std::string, std::reference_wrapper< const double > > ModelParamMap
Definition: Model.h:280
Model::name
std::string name
The name of the model.
Definition: Model.h:285
Model::raiseMissingModelParameterCount
void raiseMissingModelParameterCount()
Definition: Model.h:260
NPbase
The auxiliary base model class for other model classes.
Definition: NPbase.h:66
NPbase::trueSM
StandardModel trueSM
Definition: NPbase.h:5731
NPbase::Br_H_inv
virtual const double Br_H_inv() const
The branching ratio of the of the Higgs into invisible particles.
Definition: NPbase.h:2645
NPbase::setFlag
virtual bool setFlag(const std::string name, const bool value)
A method to set a flag of NPbase.
Definition: NPbase.h:97
Particle
A class for particles.
Definition: Particle.h:26
Particle::is
bool is(std::string name_i) const
Definition: Particle.cpp:23
Particle::getMass
const double & getMass() const
A get method to access the particle mass.
Definition: Particle.h:61
QCD::BOTTOM
@ BOTTOM
Definition: QCD.h:329
QCD::TOP
@ TOP
Definition: QCD.h:328
QCD::STRANGE
@ STRANGE
Definition: QCD.h:327
QCD::CHARM
@ CHARM
Definition: QCD.h:326
QCD::MU
@ MU
Definition: QCD.h:314
QCD::TAU
@ TAU
Definition: QCD.h:316
QCD::getQuarks
const Particle & getQuarks(const QCD::quark q) const
A get method to access a quark as an object of the type Particle.
Definition: QCD.h:536
StandardModel::computeBrHtomumu
const double computeBrHtomumu() const
The Br in the Standard Model.
Definition: StandardModel.h:2810
StandardModel::computeGammaHTotal
const double computeGammaHTotal() const
The Higgs total width in the Standard Model.
Definition: StandardModel.h:2957
StandardModel::computeBrHtoZZ
const double computeBrHtoZZ() const
The Br in the Standard Model.
Definition: StandardModel.h:2775
StandardModel::getLeptons
const Particle & getLeptons(const QCD::lepton p) const
A get method to retrieve the member object of a lepton.
Definition: StandardModel.h:744
StandardModel::computeSigmattH
const double computeSigmattH(const double sqrt_s) const
The ttH production cross section in the Standard Model.
Definition: StandardModel.h:2686
StandardModel::computeSigmaggH
const double computeSigmaggH(const double sqrt_s) const
The ggH cross section in the Standard Model.
Definition: StandardModel.h:2371
StandardModel::Mz
double Mz
The mass of the boson in GeV.
Definition: StandardModel.h:3446
StandardModel::computeBrHtocc
const double computeBrHtocc() const
The Br in the Standard Model.
Definition: StandardModel.h:2833
StandardModel::getMz
const double getMz() const
A get method to access the mass of the boson .
Definition: StandardModel.h:753
StandardModel::computeSigmaZWF
const 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:2591
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/src/StandardModel.cpp:364
StandardModel::computeSigmaWH
const double computeSigmaWH(const double sqrt_s) const
The WH production cross section in the Standard Model.
Definition: StandardModel.h:2609
StandardModel::computeBrHtotautau
const double computeBrHtotautau() const
The Br in the Standard Model.
Definition: StandardModel.h:2821
StandardModel::computeBrHtobb
const double computeBrHtobb() const
The Br in the Standard Model.
Definition: StandardModel.h:2856
StandardModel::computeBrHtogg
const double computeBrHtogg() const
The Br in the Standard Model.
Definition: StandardModel.h:2752
StandardModel::GF
double GF
The Fermi constant in .
Definition: StandardModel.h:3448
StandardModel::computeBrHtoZga
const double computeBrHtoZga() const
The Br in the Standard Model.
Definition: StandardModel.h:2787
StandardModel::computeSigmaZH
const double computeSigmaZH(const double sqrt_s) const
The ZH production cross section in the Standard Model.
Definition: StandardModel.h:2646
StandardModel::computeSigmaWF
const double computeSigmaWF(const double sqrt_s) const
The W fusion contribution to higgs-production cross section in the Standard Model.
Definition: StandardModel.h:2529
StandardModel::computeBrHtogaga
const double computeBrHtogaga() const
The Br in the Standard Model.
Definition: StandardModel.h:2799
StandardModel::lambda
double lambda
The CKM parameter in the Wolfenstein parameterization.
Definition: StandardModel.h:3461
StandardModel::computeSigmaZF
const double computeSigmaZF(const double sqrt_s) const
The Z fusion contribution to higgs-production cross section in the Standard Model.
Definition: StandardModel.h:2562
StandardModel::c02
const double c02() const
The square of the cosine of the weak mixing angle defined without weak radiative corrections.
Definition: StandardModel/src/StandardModel.cpp:1014
StandardModel::Mw_tree
const double Mw_tree() const
The tree-level mass of the boson, .
Definition: StandardModel/src/StandardModel.cpp:997
StandardModel::cW2
virtual const 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/src/StandardModel.cpp:1076
StandardModel::computeBrHtoss
const double computeBrHtoss() const
The Br in the Standard Model.
Definition: StandardModel.h:2845
StandardModel::mHl
double mHl
The Higgs mass in GeV.
Definition: StandardModel.h:3451
StandardModel::ale
double ale
The fine-structure constant .
Definition: StandardModel.h:3449
StandardModel::AlsMz
double AlsMz
The strong coupling constant at the Z-boson mass, .
Definition: StandardModel.h:3445
StandardModel::PostUpdate
virtual bool PostUpdate()
The post-update method for StandardModel.
Definition: StandardModel/src/StandardModel.cpp:241
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/src/StandardModel.cpp:278
StandardModel::v
const double v() const
The Higgs vacuum expectation value.
Definition: StandardModel/src/StandardModel.cpp:989
StandardModel::cW2
virtual const double cW2() const
Definition: StandardModel/src/StandardModel.cpp:1081
StandardModel::computeBrHtoWW
const double computeBrHtoWW() const
The Br in the Standard Model.
Definition: StandardModel.h:2763