24 br_tautau(19981, 2, 0.),
26 br_mumu(19981, 2, 0.),
29 GammaHtot_SM(19981, 2, 0.),
30 log_cs_ggH_8(200, 2, 0.),
31 log_cs_VBF_8(200, 2, 0.),
32 log_cs_WH_8(200, 2, 0.),
33 log_cs_ZH_8(200, 2, 0.),
34 log_cs_ggH_13(200, 2, 0.),
35 log_cs_VBF_13(200, 2, 0.),
36 log_cs_WH_13(200, 2, 0.),
37 log_cs_ZH_13(200, 2, 0.),
38 log_cs_ttH_8(200, 2, 0.),
39 log_cs_ttH_13(200, 2, 0.),
40 log_cs_bbH_8(200, 2, 0.),
41 log_cs_bbH_13(200, 2, 0.),
42 log_cs_ggA_8(200, 2, 0.),
43 log_cs_ttA_8(200, 2, 0.),
44 log_cs_bbA_8(200, 2, 0.),
45 log_cs_ggA_13(200, 2, 0.),
46 log_cs_ttA_13(200, 2, 0.),
47 log_cs_bbA_13(200, 2, 0.),
48 log_cs_ggHp_8(744, 3, 0.),
49 log_cs_ggHp_13(1104, 3, 0.),
50 csrH_top_charm_8(200, 2, 0.),
51 csrH_bottom_8(200, 2, 0.),
52 csrA_top_charm_8(200, 2, 0.),
53 csrA_bottom_8(200, 2, 0.),
54 csrH_top_charm_13(200, 2, 0.),
55 csrH_bottom_13(200, 2, 0.),
56 csrA_top_charm_13(200, 2, 0.),
57 csrA_bottom_13(200, 2, 0.),
59 CMS8_gg_phi_mumu(78, 2, 0.),
60 CMS8_bb_phi_mumu(78, 2, 0.),
61 CMS13_gg_phi_mumu(175, 2, 0.),
62 CMS13_bb_phi_mumu(175, 2, 0.),
63 ATLAS13_gg_phi_mumu(81, 2, 0.),
64 ATLAS13_bb_phi_mumu(81, 2, 0.),
65 ATLAS8_gg_phi_tautau(92, 2, 0.),
66 ATLAS8_bb_phi_tautau(92, 2, 0.),
67 ATLAS8_gg_phi_gaga(108, 2, 0.),
68 ATLAS8_pp_phi_Zga_llga(141, 2, 0.),
69 ATLAS8_gg_phi_ZZ(173, 2, 0.),
70 ATLAS8_VV_phi_ZZ(173, 2, 0.),
71 ATLAS8_gg_phi_WW(13, 2, 0.),
72 ATLAS8_VV_phi_WW(13, 2, 0.),
73 ATLAS8_gg_phi_phi1phi1(75, 2, 0.),
74 ATLAS8_gg_phi_phi1Z_bbZ(79, 2, 0.),
75 ATLAS8_gg_phi_phi1Z_tautauZ(79, 2, 0.),
76 CMS8_pp_phi_VV(172, 2, 0.),
77 CMS8_bb_phi_bb(81, 2, 0.),
78 CMS8_gg_phi_bb(88, 2, 0.),
79 CMS8_gg_phi_tautau(92, 2, 0.),
80 CMS8_bb_phi_tautau(92, 2, 0.),
81 CMS8_pp_phi_Zga_llga(101, 2, 0.),
82 CMS8_pp_phi_phi1phi1_bbbb(167, 2, 0.),
83 CMS8_pp_phi_phi1phi1_bbgaga(85, 2, 0.),
84 CMS8_gg_phi_phi1phi1_bbtautau(10, 2, 0.),
85 CMS8_pp_phi_phi1phi1_bbtautau(71, 2, 0.),
86 CMS8_gg_phi_phi1Z_bbll(16, 2, 0.),
87 CMS8_gg_phi_phi1Z_tautaull(14, 2, 0.),
88 CMS8_pp_phii_phijZ_bbll_1(28718, 3, 0.),
89 CMS8_pp_phii_phijZ_bbll_2(29050, 3, 0.),
90 CMS8_pp_phii_phijZ_tautaull_1(400, 3, 0.),
91 CMS8_pp_phii_phijZ_tautaull_2(400, 3, 0.),
92 ATLAS13_bb_phi_bb(96, 2, 0.),
93 CMS13_gg_phi_WW_heavy(71, 2, 0.),
94 CMS13_VV_phi_WW_heavy(71, 2, 0.),
95 CMS13_gg_phi_WW(281, 2, 0.),
96 CMS13_VV_phi_WW(561, 2, 0.),
97 ATLAS13_gg_phi_VV_llqq(95, 2, 0.),
98 ATLAS13_VV_phi_VV_llqq(95, 2, 0.),
100 ATLAS13_tt_phi_tt(13, 2, 0.),
101 ATLAS13_bb_phi_tt(61, 2, 0.),
103 ATLAS13_gg_phi_tautau(47, 2, 0.),
105 ATLAS13_bb_phi_tautau(47, 2, 0.),
107 ATLAS13_pp_phi_gaga(285, 2, 0.),
108 ATLAS13_gg_phi_Zga_llga(1590, 2, 0.),
109 ATLAS13_gg_phi_Zga_qqga(291, 2, 0.),
112 ATLAS13_gg_phi_ZZ_llllnunu(359, 2, 0.),
113 ATLAS13_VV_phi_ZZ_llllnunu(359, 2, 0.),
114 ATLAS13_gg_phi_ZZ_qqllnunu(271, 2, 0.),
115 ATLAS13_VV_phi_ZZ_qqllnunu(271, 2, 0.),
118 ATLAS13_gg_phi_WW_enumunu(72, 2, 0.),
119 ATLAS13_VV_phi_WW_enumunu(72, 2, 0.),
120 ATLAS13_gg_phi_WW_lnuqq(271, 2, 0.),
121 ATLAS13_VV_phi_WW_lnuqq(271, 2, 0.),
122 ATLAS13_pp_phi_VV_qqqq(181, 2, 0.),
124 ATLAS13_pp_phi_phi1phi1_bbbb(476, 2, 0.),
126 ATLAS13_pp_phi_phi1phi1_bbgaga(76, 2, 0.),
128 ATLAS13_pp_phi_phi1phi1_bbtautau_1(136, 2, 0.),
129 ATLAS13_pp_phi_phi1phi1_bbtautau_2(41, 2, 0.),
130 ATLAS13_pp_phi_phi1phi1_bbWW(51, 2, 0.),
131 ATLAS13_gg_phi_phi1phi1_gagaWW(25, 2, 0.),
132 ATLAS13_gg_phi_phi1Z_bbZ(45, 2, 0.),
133 ATLAS13_bb_phi_phi1Z_bbZ(45, 2, 0.),
136 ATLAS13_gg_phii_phijZ_bbZ(1711, 3, 0.),
137 ATLAS13_bb_phii_phijZ_bbZ(1711, 3, 0.),
138 ATLAS13_gg_phii_phijZ_WWZ(1326, 3, 0.),
139 CMS13_gg_phi_phi1Z_tautaull(19, 2, 0.),
140 CMS13_pp_phi2_bb_light(61, 2, 0.),
141 CMS13_pp_phi3_bb_light(61, 2, 0.),
142 CMS13_tt_phi2_tt(31, 2, 0.),
143 CMS13_tt_phi3_tt(31, 2, 0.),
144 CMS13_pp_phi_bb(66, 2, 0.),
145 CMS13_bb_phi_bb(101, 2, 0.),
150 CMS13_gg_phi_tautau(172, 2, 0.),
151 CMS13_bb_phi_tautau(172, 2, 0.),
152 CMS13_gg_phi_gaga(176, 2, 0.),
153 CMS13_gg_phi_Zga(366, 2, 0.),
154 CMS13_pp_phi_ZZ_llqqnunull(288, 2, 0.),
155 CMS13_pp_phi_ZZ_qqnunu(301, 2, 0.),
156 CMS13_ggVV_phi_WW_lnulnu(81, 2, 0.),
157 CMS13_pp_phi_WW_lnuqq(341, 2, 0.),
158 CMS13_pp_phi_phi1phi1_bbbb_1(95, 2, 0.),
160 CMS13_pp_phi_phi1phi1_bbbb_2(41, 2, 0.),
162 CMS13_pp_phi_phi1phi1_bbgaga(74, 2, 0.),
163 CMS13_pp_phi_phi1phi1_bbtautau_1(66, 2, 0.),
164 CMS13_pp_phi_phi1phi1_bbtautau_2(311, 2, 0.),
165 CMS13_pp_phi_phi1phi1_bbVV(65, 2, 0.),
166 CMS13_gg_phi_phi1Z_bbZ_1(79, 2, 0.),
167 CMS13_gg_phi_phi1Z_bbZ_2(121, 2, 0.),
168 CMS13_bb_phi_phi1Z_bbZ_1(79, 2, 0.),
169 CMS13_bb_phi_phi1Z_bbZ_2(121, 2, 0.),
170 CMS13_pp_phi_phi1phi1_4WOr2W2tauOr4tau(76, 2, 0.),
171 CMS13_pp_phi_phi1phi1_bbWW_qqlnu(55, 2, 0.),
172 CMS13_pp_phi_phi1phi1_bbZZ_lljj(149, 2, 0.),
173 CMS13_pp_phi_phi1phi1_bbZZ_llnunu(151, 2, 0.),
174 CMS13_pp_phi_phi1phi1_bbWWorbbtautau(75, 2, 0.),
175 CMS13_pp_phi_phi1phi1_bbWWorbbtautau_low_masses(65, 2, 0.),
176 ATLAS8_pp_Hpm_taunu(83, 2, 0.),
177 ATLAS8_pp_Hpm_tb(41, 2, 0.),
178 CMS8_pp_Hp_taunu(43, 2, 0.),
179 CMS8_pp_Hp_tb(43, 2, 0.),
180 ATLAS13_pp_Hpm_taunu(192, 2, 0.),
182 ATLAS13_pp_Hpm_tb(181, 2, 0.),
184 CMS13_pp_Hpm_taunu(585, 2, 0.),
185 CMS13_pp_Hpm_tb(281, 2, 0.),
187 CMS13_pp_h_phi3phi3_mumutautau(48, 2, 0.),
190 CMS13_pp_h_phi3phi3_bbtautau(11, 2, 0.),
193 CMS13_pp_h_phi3phi3_bbmumu(95, 2, 0.),
195 CMS13_pp_h_phi23Z_mumull(352, 2, 0.),
197 CMS13_pp_h_phi23phi23_mumumumu(462, 2, 0.),
199 CMS13_pp_h_phi3phi3_gagagaga(95, 2, 0.),
201 CMS13_pp_h_phi3phi3_tautautautau(12, 2, 0.),
203 CMS13_pp_phi2_gaga(81, 2, 0.),
205 CMS13_pp_bbphi3_bbtautau(10, 2, 0.),
207 CMS13_pp_ttphi3_ttmumu(220, 2, 0.),
209 CMS13_pp_ttphi3_tttautau(18, 2, 0.),
211 CMS13_Vh_h_phi3phi3_bbbb(10, 2, 0.),
213 ATLAS13_pp_h_phi3phi3_bbmumu(47, 2, 0.),
215 ATLAS13_gg_h_phi23phi23_mumumumu(110, 2, 0.),
217 ATLAS13_gg_h_phi23Z_mumull(16, 2, 0.),
219 ATLAS13_Vh_h_phi23phi23_bbbb(41, 2, 0.),
221 ATLAS13_Zh_h_phi23phi23_bbbb(7, 2, 0.),
223 ATLAS13_pp_h_phi23phi23_bbmumu_old(101, 2, 0.),
225 ATLAS13_pp_h_phi23phi23_gagagg(13, 2, 0.),
227 ATLAS13_pp_phi2_gaga_low(89, 2, 0.),
229 ATLAS13_pp_ttphi3_ttmumu(37, 2, 0.),
231 ATLAS13_pp_h_phi3phi3_gagagaga(53, 2, 0.),
233 ATLAS13_pp_h_phi3phi3_bbtautau(8, 2, 0.),
235 ATLAS13_gg_phi3_tautau_low(12, 2, 0.),
237 ATLAS8_pp_h_phi3phi3_gagagaga(105, 2, 0.),
239 ATLAS8_gg_h_phi3phi3_tautautautau(101, 2, 0.),
241 CMS8_pp_h_phi3phi3_tautautautau(6, 2, 0.),
243 CMS8_pp_h_phi3phi3_bbmumu(76, 2, 0.),
245 CMS8_pp_h_phi3phi3_mumutautau(96, 2, 0.),
247 CMS8_pp_phi2_gaga(300, 2, 0.),
249 CMS8_pp_bbphi3_bbtautau(12, 2, 0.),
251 CMS8_pp_bbphi3_bbmumu(36, 2, 0.),
253 LEP209_phi2Z_gagaZ(49, 2, 0.),
254 LEP209_phi2Z_bbZ(217, 2, 0.),
255 LEP209_phi2Z_tautauZ(233, 2, 0.),
256 LEP209_phiEphi3_bbbb(324, 3, 0.),
257 LEP209_phiEphi3_tautautautau(365, 3, 0.),
259 CMS8_t_Hpb_csb(7, 2, 0.),
260 CMS8_t_Hpb_taunub(8, 2, 0.),
261 CMS8_t_Hpb_cbb(7, 2, 0.),
262 CMS13_t_Hpb_WAb_Wmumub(95, 2, 0.),
263 CMS13_t_Hpb_csb(8, 2, 0.),
264 ATLAS8_t_Hpb_taunub(9, 2, 0.),
265 ATLAS13_t_Hpb_cbb(11, 2, 0.),
266 ATLAS13_t_Hpb_WAb_Wmumub(31, 2, 0.),
267 ATLAS13_t_Hpb_csb(12, 2, 0.),
270 LEP209_HpHm_taunutaunu(105, 2, 0.),
271 LEP209_HpHm_qqqq(105, 2, 0.),
272 OPAL209_HpHm_qqtaunu(87, 2, 0.),
273 OPAL172_HpHm_qqtaunu(41, 2, 0.),
274 OPAL209_HpHm_AWstAWst_bbbb(2835, 3, 0.),
275 OPAL209_HpHm_taunuAWst_taunubb(2835, 3, 0.),
277 ATLAS13_pp_HpHm_taunutaunu(19, 2, 0.),
278 ATLAS13_pp_HpHm_munumunu(11, 2, 0.),
279 CMS13_pp_HpHm_taunutaunu(16, 2, 0.),
280 CMS13_pp_HpHm_munumunu(22, 2, 0.),
282 csrHpHm_pp_13_LO(12, 2, 0.),
283 csrHpHm_pp_13_NLO(31, 2, 0.),
285 integral_x2_1mx_G_log(62500, 3, 0.),
286 integral_x2_1px_G_log(62500, 3, 0.),
287 integral_x2_G_log(62500, 3, 0.),
288 integral_x_1mx2_G_log(62500, 3, 0.),
289 integral_x_1mx_1px_G_log(62500, 3, 0.),
291 integral_x2_1mx_G_variable_set_1_log(2000, 2, 0.),
292 integral_x2_G_variable_set_1_log(2000, 2, 0.),
293 integral_x_1mx2_G_variable_set_0_log(2000, 2, 0.),
294 integral_perturb_ATHDM_values_log(2000, 2, 0.),
296 integral_Hp_A_Wstar(2835, 3, 0.),
299 arraybsgamma(1111, 3, 0.),
303 Mu_GTHDM(3,3,0.), Md_GTHDM(3,3,0.), Ml_GTHDM(3,3,0.),
304 Nu_GTHDM(3,3,0.), Nd_GTHDM(3,3,0.), Nl_GTHDM(3,3,0.),
309 myGTHDM(static_cast<const GeneralTHDM*> (&SM_i)),
320 const int NumPar,
const double params[])
const {
324 for(
int j=0; j<NumPar; j++)
325 bCache &= (params[j] == cache[j][i].real());
326 if (bCache)
return i;
332 const int NumPar,
const double params[])
const {
336 for(
int j=0; j<NumPar; j++)
337 bCache &= (params[j] == cache[j][i]);
338 if (bCache)
return i;
345 const double params[],
const gslpp::complex newResult)
const {
348 for(
int j=0; j<NumPar+1; j++)
349 cache[j][i] = cache[j][i-1];
352 for(
int j=0; j<NumPar; j++) {
353 cache[j][0] = gslpp::complex(params[j], 0.0,
false);
354 cache[NumPar][0] = newResult;
359 const double params[],
const double newResult)
const {
362 for(
int j=0; j<NumPar+1; j++)
363 cache[j][i] = cache[j][i-1];
366 for(
int j=0; j<NumPar; j++) {
367 cache[j][0] = params[j];
368 cache[NumPar][0] = newResult;
378 double params[] = {MZ2, MW2, mHh2};
384 gslpp::complex newResult =
PV.
B0(MZ2, 0., MW2, mHh2);
392 double params[] = {MZ2, MW2, mHl2};
398 gslpp::complex newResult =
PV.
B0(MZ2, 0., MW2, mHl2);
406 double params[] = {MZ2, mHh2};
412 gslpp::complex newResult =
PV.
B0(MZ2, 0., MZ2, mHh2);
420 double params[] = {MZ2, mHl2};
426 gslpp::complex newResult =
PV.
B0(MZ2, 0., MZ2, mHl2);
434 double params[] = {MZ2, MW2, mHh2};
440 gslpp::complex newResult =
PV.
B0(MZ2, MW2, MW2, mHh2);
448 double params[] = {MZ2, MW2, mHl2};
454 gslpp::complex newResult =
PV.
B0(MZ2, MW2, MW2, mHl2);
462 double params[] = {MZ2, mHh2};
468 gslpp::complex newResult =
PV.
B0(MZ2, MZ2, MZ2, mHh2);
476 double params[] = {MZ2, mHl2};
482 gslpp::complex newResult =
PV.
B0(MZ2, MZ2, MZ2, mHl2);
490 double params[] = {MZ2, mHl2};
496 gslpp::complex newResult =
PV.
B0(MZ2, 0.0, 0.0, mHl2);
504 double params[] = {MZ2, mHh2};
510 gslpp::complex newResult =
PV.
B0(MZ2, 0.0, 0.0, mHh2);
518 double params[] = {MZ2,
mHp2, mHl2};
524 gslpp::complex newResult =
PV.
B0(MZ2, 0.0,
mHp2, mHl2);
532 double params[] = {MZ2,
mHp2, mHh2};
538 gslpp::complex newResult =
PV.
B0(MZ2, 0.0,
mHp2, mHh2);
546 double params[] = {MZ2, mA2, mHl2};
552 gslpp::complex newResult =
PV.
B0(MZ2, 0.0, mA2, mHl2);
560 double params[] = {MZ2, mA2, mHh2};
566 gslpp::complex newResult =
PV.
B0(MZ2, 0.0, mA2, mHh2);
574 double params[] = {MZ2, mHl2};
580 gslpp::complex newResult =
PV.
B0(MZ2, mHl2, 0.0, 0.0);
588 double params[] = {MZ2, mHl2,
mHp2};
594 gslpp::complex newResult =
PV.
B0(MZ2, mHl2, 0.0,
mHp2);
602 double params[] = {MZ2, mHl2, mA2};
608 gslpp::complex newResult =
PV.
B0(MZ2, mHl2, 0.0, mA2);
616 double params[] = {MZ2, mHl2};
622 gslpp::complex newResult =
PV.
B0(MZ2, mHl2, mHl2, mHl2);
630 double params[] = {MZ2, mHl2, mHh2};
636 gslpp::complex newResult =
PV.
B0(MZ2, mHl2, mHh2, mHl2);
644 double params[] = {MZ2, mHl2, mHh2};
650 gslpp::complex newResult =
PV.
B0(MZ2, mHl2, mHh2, mHh2);
658 double params[] = {MZ2, mHl2,
mHp2};
664 gslpp::complex newResult =
PV.
B0(MZ2, mHl2,
mHp2,
mHp2);
672 double params[] = {MZ2, mHl2, mA2};
678 gslpp::complex newResult =
PV.
B0(MZ2, mHl2, mA2, mA2);
686 double params[] = {MZ2, mHh2};
692 gslpp::complex newResult =
PV.
B0(MZ2, mHh2, 0.0, 0.0);
700 double params[] = {MZ2, mHh2,
mHp2};
706 gslpp::complex newResult =
PV.
B0(MZ2, mHh2, 0.0,
mHp2);
714 double params[] = {MZ2, mHh2, mA2};
720 gslpp::complex newResult =
PV.
B0(MZ2, mHh2, 0.0, mA2);
728 double params[] = {MZ2, mHh2, mHl2};
734 gslpp::complex newResult =
PV.
B0(MZ2, mHh2, mHl2, mHl2);
742 double params[] = {MZ2, mHh2, mHl2};
748 gslpp::complex newResult =
PV.
B0(MZ2, mHh2, mHh2, mHl2);
756 double params[] = {MZ2, mHh2};
762 gslpp::complex newResult =
PV.
B0(MZ2, mHh2, mHh2, mHh2);
770 double params[] = {MZ2, mHh2,
mHp2};
776 gslpp::complex newResult =
PV.
B0(MZ2, mHh2,
mHp2,
mHp2);
784 double params[] = {MZ2, mHh2, mA2};
790 gslpp::complex newResult =
PV.
B0(MZ2, mHh2, mA2, mA2);
798 double params[] = {MZ2,
mHp2, mHl2};
804 gslpp::complex newResult =
PV.
B0(MZ2,
mHp2, 0.0, mHl2);
812 double params[] = {MZ2,
mHp2, mHh2};
818 gslpp::complex newResult =
PV.
B0(MZ2,
mHp2, 0.0, mHh2);
826 double params[] = {MZ2,
mHp2, mHl2};
832 gslpp::complex newResult =
PV.
B0(MZ2,
mHp2,
mHp2, mHl2);
840 double params[] = {MZ2,
mHp2, mHh2};
846 gslpp::complex newResult =
PV.
B0(MZ2,
mHp2,
mHp2, mHh2);
854 double params[] = {MZ2, mA2, mHl2};
860 gslpp::complex newResult =
PV.
B0(MZ2, mA2, 0.0, mHl2);
868 double params[] = {MZ2, mA2, mHh2};
874 gslpp::complex newResult =
PV.
B0(MZ2, mA2, 0.0, mHh2);
882 double params[] = {MZ2, mA2, mHl2};
888 gslpp::complex newResult =
PV.
B0(MZ2, mA2, mA2, mHl2);
896 double params[] = {MZ2, mA2, mHh2};
902 gslpp::complex newResult =
PV.
B0(MZ2, mA2, mA2, mHh2);
914 double params[] = {MZ2, mHl2};
920 gslpp::complex newResult =
PV.
B0p(MZ2, 0.0, 0.0, mHl2);
928 double params[] = {MZ2, mHh2};
934 gslpp::complex newResult =
PV.
B0p(MZ2, 0.0, 0.0, mHh2);
942 double params[] = {MZ2,
mHp2, mHl2};
948 gslpp::complex newResult =
PV.
B0p(MZ2, 0.0,
mHp2, mHl2);
956 double params[] = {MZ2,
mHp2, mHh2};
962 gslpp::complex newResult =
PV.
B0p(MZ2, 0.0,
mHp2, mHh2);
970 double params[] = {MZ2,
mHp2, mA2};
976 gslpp::complex newResult =
PV.
B0p(MZ2, 0.0,
mHp2, mA2);
984 double params[] = {MZ2, mA2, mHl2};
990 gslpp::complex newResult =
PV.
B0p(MZ2, 0.0, mA2, mHl2);
998 double params[] = {MZ2, mA2, mHh2};
1004 gslpp::complex newResult =
PV.
B0p(MZ2, 0.0, mA2, mHh2);
1012 double params[] = {MZ2, mHl2};
1018 gslpp::complex newResult =
PV.
B0p(MZ2, mHl2, 0.0, 0.0);
1026 double params[] = {MZ2, mHl2,
mHp2};
1032 gslpp::complex newResult =
PV.
B0p(MZ2, mHl2, 0.0,
mHp2);
1040 double params[] = {MZ2, mHl2, mA2};
1046 gslpp::complex newResult =
PV.
B0p(MZ2, mHl2, 0.0, mA2);
1054 double params[] = {MZ2, mHl2};
1060 gslpp::complex newResult =
PV.
B0p(MZ2, mHl2, mHl2, mHl2);
1068 double params[] = {MZ2, mHl2, mHh2};
1074 gslpp::complex newResult =
PV.
B0p(MZ2, mHl2, mHh2, mHl2);
1082 double params[] = {MZ2, mHl2, mHh2};
1088 gslpp::complex newResult =
PV.
B0p(MZ2, mHl2, mHh2, mHh2);
1096 double params[] = {MZ2, mHl2,
mHp2};
1110 double params[] = {MZ2, mHl2, mA2};
1116 gslpp::complex newResult =
PV.
B0p(MZ2, mHl2, mA2, mA2);
1124 double params[] = {MZ2, mHh2};
1130 gslpp::complex newResult =
PV.
B0p(MZ2, mHh2, 0.0, 0.0);
1138 double params[] = {MZ2, mHh2,
mHp2};
1144 gslpp::complex newResult =
PV.
B0p(MZ2, mHh2, 0.0,
mHp2);
1152 double params[] = {MZ2, mHh2, mA2};
1158 gslpp::complex newResult =
PV.
B0p(MZ2, mHh2, 0.0, mA2);
1166 double params[] = {MZ2, mHh2, mHl2};
1172 gslpp::complex newResult =
PV.
B0p(MZ2, mHh2, mHl2, mHl2);
1180 double params[] = {MZ2, mHh2, mHl2};
1186 gslpp::complex newResult =
PV.
B0p(MZ2, mHh2, mHh2, mHl2);
1194 double params[] = {MZ2, mHh2};
1200 gslpp::complex newResult =
PV.
B0p(MZ2, mHh2, mHh2, mHh2);
1208 double params[] = {MZ2, mHh2,
mHp2};
1222 double params[] = {MZ2, mHh2, mA2};
1228 gslpp::complex newResult =
PV.
B0p(MZ2, mHh2, mA2, mA2);
1236 double params[] = {MZ2,
mHp2, mHl2};
1242 gslpp::complex newResult =
PV.
B0p(MZ2,
mHp2, 0.0, mHl2);
1250 double params[] = {MZ2,
mHp2, mHh2};
1256 gslpp::complex newResult =
PV.
B0p(MZ2,
mHp2, 0.0, mHh2);
1264 double params[] = {MZ2,
mHp2, mA2};
1270 gslpp::complex newResult =
PV.
B0p(MZ2,
mHp2, 0.0, mA2);
1278 double params[] = {MZ2,
mHp2, mHl2};
1292 double params[] = {MZ2,
mHp2, mHh2};
1306 double params[] = {MZ2, mA2, mHl2};
1312 gslpp::complex newResult =
PV.
B0p(MZ2, mA2, 0.0, mHl2);
1320 double params[] = {MZ2, mA2, mHh2};
1326 gslpp::complex newResult =
PV.
B0p(MZ2, mA2, 0.0, mHh2);
1334 double params[] = {MZ2, mA2,
mHp2};
1340 gslpp::complex newResult =
PV.
B0p(MZ2, mA2, 0.0,
mHp2);
1348 double params[] = {MZ2, mA2, mHl2};
1354 gslpp::complex newResult =
PV.
B0p(MZ2, mA2, mA2, mHl2);
1362 double params[] = {MZ2, mA2, mHh2};
1368 gslpp::complex newResult =
PV.
B0p(MZ2, mA2, mA2, mHh2);
1378 double params[] = {MZ2, mA2,
mHp2};
1384 gslpp::complex newResult =
PV.
B00(MZ2, 0., mA2,
mHp2);
1392 double params[] = {MZ2, mHh2, mA2};
1398 gslpp::complex newResult =
PV.
B00(MZ2, 0., mHh2, mA2);
1406 double params[] = {MZ2, mHh2,
mHp2};
1412 gslpp::complex newResult =
PV.
B00(MZ2, 0., mHh2,
mHp2);
1420 double params[] = {MZ2, mHl2, mA2};
1426 gslpp::complex newResult =
PV.
B00(MZ2, 0., mHl2, mA2);
1434 double params[] = {MZ2, mHl2,
mHp2};
1440 gslpp::complex newResult =
PV.
B00(MZ2, 0., mHl2,
mHp2);
1448 double params[] = {MZ2,
mHp2};
1462 double params[] = {MZ2, MW2, mHh2};
1468 gslpp::complex newResult =
PV.
B00(MZ2, MW2, MW2, mHh2);
1476 double params[] = {MZ2, MW2, mHl2};
1482 gslpp::complex newResult =
PV.
B00(MZ2, 0., MW2, mHl2);
1490 double params[] = {MZ2, mHh2};
1496 gslpp::complex newResult =
PV.
B00(MZ2, 0., MZ2, mHh2);
1504 double params[] = {MZ2, mHl2};
1510 gslpp::complex newResult =
PV.
B00(MZ2, 0., MZ2, mHl2);
1518 double params[] = {MZ2, MW2, mA2,
mHp2};
1524 gslpp::complex newResult =
PV.
B00(MZ2, MW2, mA2,
mHp2);
1532 double params[] = {MZ2, MW2, mHh2,
mHp2};
1538 gslpp::complex newResult =
PV.
B00(MZ2, MW2, mHh2,
mHp2);
1546 double params[] = {MZ2, MW2, mHl2,
mHp2};
1552 gslpp::complex newResult =
PV.
B00(MZ2, MW2, mHl2,
mHp2);
1560 double params[] = {MZ2, MW2,
mHp2};
1574 double params[] = {MZ2, MW2, mHh2};
1580 gslpp::complex newResult =
PV.
B00(MZ2, MW2, MW2, mHh2);
1588 double params[] = {MZ2, MW2, mHl2};
1594 gslpp::complex newResult =
PV.
B00(MZ2, MW2, MW2, mHl2);
1602 double params[] = {MZ2, mHh2, mA2};
1608 gslpp::complex newResult =
PV.
B00(MZ2, MZ2, mHh2, mA2);
1616 double params[] = {MZ2, mHl2, mA2};
1622 gslpp::complex newResult =
PV.
B00(MZ2, MZ2, mHl2, mA2);
1630 double params[] = {MZ2,
mHp2};
1644 double params[] = {MZ2, mHh2};
1650 gslpp::complex newResult =
PV.
B00(MZ2, MZ2, MZ2, mHh2);
1658 double params[] = {MZ2, mHl2};
1664 gslpp::complex newResult =
PV.
B00(MZ2, MZ2, MZ2, mHl2);
1675 std::ifstream INfile;
1676 std::string lineTab;
1677 INfile.open( filename.c_str() );
1679 std::cout<<
"error: in GeneralTHDMcache, table doesn't exist!"<<std::endl;
1682 gslpp::matrix<double> arrayTab(rowN, colN, 0.);
1687 while(INfile.good()){
1688 while(getline(INfile, lineTab)){
1689 if( lineTab[0]==
'#' )
continue;
1691 std::istringstream streamTab(lineTab);
1693 while(streamTab >>v){
1694 arrayTab.assign(a,b,v);
1711 int rowN=arrayTab.size_i();
1713 double xmin = arrayTab(0,0);
1714 double xmax = arrayTab(rowN-1,0);
1715 double interval = arrayTab(1,0)-arrayTab(0,0);
1716 int Nintervals = (x-xmin)/interval;
1719 std::cout<<
"warning: your table parameter value is smaller than the minimum allowed value for m = "<<x<<
"\n while xmin = "<<xmin<<std::endl;
1728 std::cout<<
"warning: your table parameter value is greater than the maximum allowed value for m = "<<x<<
"\n while xmax = "<<xmax<<std::endl;
1735 return arrayTab(rowN-1,1);
1738 return ((arrayTab(Nintervals+1,1)-arrayTab(Nintervals,1))/(arrayTab(Nintervals+1,0)
1739 -arrayTab(Nintervals,0))*(x-arrayTab(Nintervals,0))+arrayTab(Nintervals,1));
1747 int rowN=arrayTab.size_i();
1749 double xmin = arrayTab(0,0);
1750 double xmax = arrayTab(rowN-1,0);
1751 double test = x - (xmax + xmin) / 2;
1756 for(row = 1; row <= rowN-1; row++)
1758 if(arrayTab(row,0) - x > 0)
1764 for(row = rowN-2; row >= 0; row--)
1766 if(x - arrayTab(row,0) > 0)
1778 return arrayTab(rowN-1,1);
1780 return ((arrayTab(row,1)-arrayTab(row-1,1))/(arrayTab(row,0)
1781 -arrayTab(row-1,0))*(x-arrayTab(row-1,0))+arrayTab(row-1,1));
1788 int rowN=arrayTab.size_i();
1790 double xmin = arrayTab(0,0);
1791 double xmax = arrayTab(rowN-1,0);
1792 double ymin = arrayTab(0,1);
1793 double ymax = arrayTab(rowN-1,1);
1794 double intervalx = arrayTab(1,0)-arrayTab(0,0);
1797 while(arrayTab(i,1)-arrayTab(i-1,1)==0&&i<30000);
1798 double intervaly = arrayTab(i,1)-arrayTab(i-1,1);
1799 int Nintervalsx = (x-xmin)/intervalx;
1800 int Nintervalsy = (y-ymin)/intervaly;
1801 if(x<xmin||x>xmax||y<ymin||y>ymax){
1802 std::cout<<
"warning: the parameter point lies outside the table on interpolate2D"<<std::endl;
1806 double x1=arrayTab(i*Nintervalsy+Nintervalsx,0);
1807 double x2=arrayTab(i*Nintervalsy+Nintervalsx+1,0);
1808 double y1=arrayTab(i*Nintervalsy+Nintervalsx,1);
1809 double y2=arrayTab(i*(Nintervalsy+1)+Nintervalsx,1);
1810 return (arrayTab(i*Nintervalsy+Nintervalsx,2) * (x2-x) * (y2-y)
1811 +arrayTab(i*Nintervalsy+Nintervalsx+1,2) * (x-x1) * (y2-y)
1812 +arrayTab(i*(Nintervalsy+1)+Nintervalsx,2) * (x2-x) * (y-y1)
1813 +arrayTab(i*(Nintervalsy+1)+Nintervalsx+1,2) * (x-x1) * (y-y1))
1821 int rowN = arrayTab.size_i();
1823 double xmin = arrayTab(0,0);
1824 double xmax = arrayTab(rowN-1,0);
1825 double ymin = arrayTab(0,1);
1826 double ymax = arrayTab(rowN-1,1);
1827 double intervaly = arrayTab(1,1) - arrayTab(0,1);
1830 while (arrayTab(i,0) - arrayTab(i-1,0) == 0 && i < 30000);
1831 double intervalx = arrayTab(i,0) - arrayTab(i-1,0);
1833 int Nintervalsy = (y - ymin) / intervaly;
1834 int Nintervalsx = (x - xmin) / intervalx;
1835 if (x < xmin || x > xmax || y < ymin || y > ymax) {
1836 std::cout <<
"Warning: the parameter point lies outside the table on interpolate2Dv2" << std::endl;
1840 double x1 = arrayTab(Nintervalsy + Nintervalsx*i, 0);
1841 double x2 = arrayTab(Nintervalsy + (Nintervalsx+1)*i, 0);
1842 double y1 = arrayTab(Nintervalsy + Nintervalsx*i, 1);
1843 double y2 = arrayTab(Nintervalsy + 1 + Nintervalsx*i, 1);
1845 return ((arrayTab(Nintervalsy + Nintervalsx*i, 2) * (x2 - x) * (y2 - y)
1846 + arrayTab(Nintervalsy + (Nintervalsx + 1)*i, 2) * (x - x1) * (y2 - y)
1847 + arrayTab(Nintervalsy + 1 + Nintervalsx*i, 2) * (x2 - x) * (y - y1)
1848 + arrayTab((Nintervalsy + 1) + (Nintervalsx + 1)*i, 2) * (x - x1) * (y - y1))
1849 / ((x2 - x1) * (y2 - y1)));
1859 int rowN=arrayTab.size_i();
1861 double xmin = arrayTab(0,0);
1862 double xmax = arrayTab(rowN-1,0);
1863 double ymin = arrayTab(0,1);
1864 double ymax = arrayTab(rowN-1,1);
1865 double intervaly = arrayTab(1,1)-arrayTab(0,1);
1868 while(arrayTab(i,0)-arrayTab(i-1,0)==0&&i<30000);
1869 double intervalx = arrayTab(i,0)-arrayTab(i-1,0);
1870 int Nintervalsx = (x-xmin)/intervalx;
1871 int Nintervalsy = (y-ymin)/intervaly;
1872 if(x<xmin||x>xmax||y<ymin||y>ymax){
1873 std::cout<<
"warning: the parameter point lies outside the table on interpolate2DtriangularData"<<std::endl;
1878 int Nx1y1=i*Nintervalsx-(Nintervalsx*(Nintervalsx+1)/2) +Nintervalsy;
1879 int Nx1y2=i*Nintervalsx-(Nintervalsx*(Nintervalsx+1)/2) +Nintervalsy+1;
1880 int Nx2y1=i*(Nintervalsx+1)-((Nintervalsx+1)*((Nintervalsx+1)+1)/2) +Nintervalsy;
1881 int Nx2y2=i*(Nintervalsx+1)-((Nintervalsx+1)*((Nintervalsx+1)+1)/2) +Nintervalsy+1;
1916 double y1 = arrayTab(Nx1y1,1);
1917 double y2 = arrayTab(Nx1y2,1);
1918 double x1 = arrayTab(Nx1y1,0);
1919 double x2 = arrayTab(Nx2y1,0);
1922 return (arrayTab(Nx2y2,2) * (y2-y) * (x2-x)
1923 +arrayTab(Nx2y1,2) * (y-y1) * (x2-x)
1924 +arrayTab(Nx1y2,2) * (y2-y) * (x-x1)
1925 +arrayTab(Nx1y1,2) * (y-y1) * (x-x1))
1934 unsigned int interpol, Npoints, Nbreak;
1935 std::vector<double> row, aux;
1936 double xdx = x / dx;
1937 double ydy = y / dy;
1939 if(xdx - ((
int) xdx) == 0. && ydy - ((
int) ydy) == 0.){
1943 else if(xdx - ((
int) xdx) == 0.){
1947 else if(ydy - ((
int) ydy) == 0.){
1958 for(
unsigned int i = 0; Nbreak < Npoints; i++)
1960 if(i == arrayTab.size_i())
1963 if(std::fabs(x-arrayTab(i,0)) < dx && std::fabs(y-arrayTab(i,1)) < dy)
1966 Nbreak = Nbreak + 1;
1976 return (arrayTab(row[0],2));
1978 return (arrayTab(row[0],2) + (arrayTab(row[1],2)-arrayTab(row[0],2)) * (y-arrayTab(row[0],1)) / dy);
1980 return (arrayTab(row[0],2) + (arrayTab(row[1],2)-arrayTab(row[0],2)) * (x-arrayTab(row[0],0)) / dx);
1982 aux.push_back(arrayTab(row[0],2) + (arrayTab(row[1],2)-arrayTab(row[0],2)) * (y-arrayTab(row[0],1)) / dy);
1983 aux.push_back(arrayTab(row[2],2) + (arrayTab(row[3],2)-arrayTab(row[2],2)) * (y-arrayTab(row[2],1)) / dy);
1984 return (aux[0] + (aux[1]-aux[0]) * (x-arrayTab(row[0],0)) / dx);
1986 aux.push_back(arrayTab(row[0],2) + (arrayTab(row[1],2)-arrayTab(row[0],2)) * (y-arrayTab(row[0],1)) / dy);
1987 aux.push_back(arrayTab(row[1],2) + (arrayTab(row[2],2)-arrayTab(row[1],2)) * (y-arrayTab(row[0],1)) / dy);
1988 return (aux[0] + (aux[1]-aux[0]) * (x-arrayTab(row[0],0)) / dx);
1990 throw std::runtime_error(
"Review applicability of the mass ranges provided to GeneralTHDMcache::interpolate2Dsteps");
2000 std::stringstream br1,br2,br3,br4,br5,br6,br7;
2001 std::stringstream dw1;
2002 std::stringstream cs1,cs2,cs3,cs4,cs5,cs6,cs7,cs8,cs9;
2003 std::stringstream cs11,cs12,cs13,cs14,cs15,cs16,cs17,cs18,cs19;
2004 std::stringstream cs20,cs21;
2005 std::stringstream csr1,csr2,csr3,csr4;
2006 std::stringstream csr11,csr12,csr13,csr14;
2007 std::stringstream ex1m6,ex1m5,ex1m4,ex1m3,ex1m2,ex1m1;
2008 std::stringstream ex1,ex2,ex3,ex4,ex5,ex6,ex7,ex8,ex9,ex10,ex11,ex12,ex13,ex14,ex15,ex16,ex17,ex18,ex19,ex20,ex21,ex22,ex23;
2009 std::stringstream ex24,ex25,ex26,ex27,ex28m1,ex28,ex29,ex30,ex31,ex32,ex33,ex34,ex35,ex36,ex37,ex38,ex39m4,ex39m3,ex39m2,ex39m1,ex39,ex40,ex41,ex42,ex43,ex43p1,ex43p2,ex44,\
2010 ex45,ex46n1,ex46n2,ex46a,ex47,ex48,ex49,ex49p2,ex50a,ex50b,ex51,ex52m2,ex52m1,ex52,ex52p1,ex52p2,ex53,ex54,ex55,ex56;
2011 std::stringstream ex57,ex58,ex59,ex60,ex61,ex62,ex63,ex64,ex65,ex66,ex67,ex67p1,ex67p2,ex67p3,ex67p4,ex67p5,ex67p6,ex68,ex69,ex70,ex71,ex72,ex73,ex74,ex75,ex76,ex77,\
2013 std::stringstream lowC01,lowC02,lowC03,lowC04,lowC05,lowC06,lowC07,lowC08, lowC09, lowC10, lowC11, lowC12;
2014 std::stringstream lowA01, lowA02, lowA03, lowA04, lowA05, lowA06, lowA07, lowA08, lowA09, lowA10, lowA11, lowA12;
2015 std::stringstream lowA801, lowA802;
2016 std::stringstream lowC801, lowC802, lowC803, lowC804, lowC805, lowC806;
2017 std::stringstream low209a, low209b, low209c, low209d, low209e;
2018 std::stringstream lowHpC801, lowHpC802, lowHpC803, lowHpC1301, lowHpC1302;
2019 std::stringstream lowHpA801, lowHpA1301, lowHpA1302, lowHpA1303;
2020 std::stringstream lowHp209a, lowHp209b, lowHp209c, lowHp172a, lowHp209d, lowHp209e;
2021 std::stringstream thint01, thint02, thint03, thint04, thint05, thint06, thint07, thint08, thint09, thint10;
2022 std::stringstream csrslepLO, csrslepNLO, susyHpA01, susyHpA02, susyHpC01, susyHpC02;
2023 std::stringstream bsg1;
2025 std::cout<<
"reading tables"<<std::endl;
2028 std::stringstream path;
2029 path << getenv(
"HEPFITTABS") <<
"/THDM/tabs/";
2030 std::string tablepath=path.str();
2032 br1 << tablepath <<
"br1.dat";
2034 br2 << tablepath <<
"br2.dat";
2036 br3 << tablepath <<
"br3.dat";
2038 br4 << tablepath <<
"br4.dat";
2040 br5 << tablepath <<
"br5.dat";
2042 br6 << tablepath <<
"br6.dat";
2044 br7 << tablepath <<
"br7.dat";
2046 dw1 << tablepath <<
"dw1.dat";
2048 cs1 << tablepath <<
"log_cs_ggH_8.dat";
2050 cs11 << tablepath <<
"log_cs_ggH_13.dat";
2052 cs2 << tablepath <<
"log_cs_VBF_8.dat";
2054 cs12 << tablepath <<
"log_cs_VBF_13.dat";
2056 cs3 << tablepath <<
"log_cs_WH_8.dat";
2058 cs13 << tablepath <<
"log_cs_WH_13.dat";
2060 cs4 << tablepath <<
"log_cs_ZH_8.dat";
2062 cs14 << tablepath <<
"log_cs_ZH_13.dat";
2064 cs5 << tablepath <<
"log_cs_ttH_8.dat";
2066 cs15 << tablepath <<
"log_cs_ttH_13.dat";
2068 cs6 << tablepath <<
"log_cs_bbH_8.dat";
2070 cs16 << tablepath <<
"log_cs_bbH_13.dat";
2072 cs7 << tablepath <<
"log_cs_ggA_8.dat";
2074 cs17 << tablepath <<
"log_cs_ggA_13.dat";
2076 cs8 << tablepath <<
"log_cs_ttA_8.dat";
2078 cs18 << tablepath <<
"log_cs_ttA_13.dat";
2080 cs9 << tablepath <<
"log_cs_bbA_8.dat";
2082 cs19 << tablepath <<
"log_cs_bbA_13.dat";
2084 cs20 << tablepath <<
"log_cs_ggHp_8.dat";
2086 cs21 << tablepath <<
"log_cs_ggHp_13.dat";
2088 csr1 << tablepath <<
"csrH_top_charm_8.dat";
2090 csr11 << tablepath <<
"csrH_top_charm_13.dat";
2092 csr2 << tablepath <<
"csrH_bottom_8.dat";
2094 csr12 << tablepath <<
"csrH_bottom_13.dat";
2096 csr3 << tablepath <<
"csrA_top_charm_8.dat";
2098 csr13 << tablepath <<
"csrA_top_charm_13.dat";
2100 csr4 << tablepath <<
"csrA_bottom_8.dat";
2102 csr14 << tablepath <<
"csrA_bottom_13.dat";
2104 csrslepLO << tablepath <<
"csrHpHm_pp_slepton_13_LO.dat";
2106 csrslepNLO << tablepath <<
"csrHpHm_pp_slepton_13_NLONLL.dat";
2110 ex1m6<< tablepath <<
"150801437_9b.dat";
2112 ex1m5<< tablepath <<
"150801437_9a.dat";
2114 ex1m4<< tablepath <<
"190703152_6b.dat";
2116 ex1m3<< tablepath <<
"190703152_6a.dat";
2118 ex1m2<< tablepath <<
"190108144_4b.dat";
2120 ex1m1<< tablepath <<
"190108144_4a.dat";
2127 ex1 << tablepath <<
"14096064_a.dat";
2129 ex2 << tablepath <<
"14096064_b.dat";
2131 ex3 << tablepath <<
"14076583.dat";
2133 ex4 << tablepath <<
"14078150.dat";
2135 ex5 << tablepath <<
"150705930_a.dat";
2137 ex6 << tablepath <<
"150705930_b.dat";
2139 ex7 << tablepath <<
"150900389_a.dat";
2141 ex8 << tablepath <<
"150900389_b.dat";
2143 ex9 << tablepath <<
"150904670.dat";
2145 ex10 << tablepath <<
"150204478_b.dat";
2147 ex11 << tablepath <<
"150204478_a.dat";
2149 ex12 << tablepath <<
"150400936.dat";
2151 ex13 << tablepath <<
"150608329.dat";
2153 ex14 << tablepath <<
"180206149.dat";
2155 ex15 << tablepath <<
"CMS-PAS-HIG-14-029_a.dat";
2157 ex16 << tablepath <<
"CMS-PAS-HIG-14-029_b.dat";
2159 ex17 << tablepath <<
"CMS-PAS-HIG-16-014.dat";
2161 ex18 << tablepath <<
"150304114.dat";
2163 ex19 << tablepath <<
"160306896.dat";
2165 ex20 << tablepath <<
"151001181_a.dat";
2167 ex21 << tablepath <<
"170700350.dat";
2169 ex22 << tablepath <<
"150404710.dat";
2171 ex23 << tablepath <<
"151001181_b.dat";
2174 ex24 << tablepath <<
"160302991_a.dat";
2176 ex25 << tablepath <<
"160302991_b.dat";
2178 ex26 << tablepath <<
"160302991_c.dat";
2180 ex27 << tablepath <<
"160302991_d.dat";
2185 ex28m1 << tablepath <<
"190702749.dat";
2189 ex28 << tablepath <<
"ATLAS_CONF_2022_008.dat";
2191 ex29 << tablepath <<
"ATLAS-CONF-2016-104_b.dat";
2203 ex30 << tablepath <<
"200212223_2a.dat";
2205 ex31 << tablepath <<
"200212223_2b.dat";
2212 ex32 << tablepath <<
"210213405.dat";
2219 ex33 << tablepath <<
"230904364.dat";
2227 ex34 << tablepath <<
"230411962.dat";
2236 ex35 << tablepath <<
"200914791_4a.dat";
2238 ex36 << tablepath <<
"200914791_4b.dat";
2243 ex37 << tablepath <<
"170809638_a.dat";
2245 ex38 << tablepath <<
"170809638_b.dat";
2249 ex39m4 << tablepath <<
"210906055_7a.dat";
2251 ex39m3 << tablepath <<
"210906055_7b.dat";
2255 ex39m2 << tablepath <<
"191201594_6c.dat";
2257 ex39m1 << tablepath <<
"191201594_6d.dat";
2266 ex39 << tablepath <<
"ATLAS_CONF_2022_066_ggF.dat";
2268 ex40 << tablepath <<
"ATLAS_CONF_2022_066_VBF.dat";
2273 ex41 << tablepath <<
"171007235_a.dat";
2275 ex42 << tablepath <<
"171007235_b.dat";
2277 ex43 << tablepath <<
"170804445.dat";
2281 ex43p1 << tablepath <<
"200414636_12a.dat";
2283 ex43p2 << tablepath <<
"200414636_12b.dat";
2288 ex44 << tablepath <<
"220207288.dat";
2292 ex45 << tablepath <<
"211211876.dat";
2296 ex46n1 << tablepath <<
"ATLAS-CONF-2021-030.dat";
2298 ex46n2 << tablepath <<
"200714811.dat";
2300 ex46a << tablepath <<
"181104671.dat";
2302 ex47 << tablepath <<
"180708567.dat";
2308 ex48 << tablepath <<
"220700230_ggF.dat";
2310 ex49 << tablepath <<
"220700230_bb.dat";
2313 ex49p2 << tablepath <<
"191011634.dat";
2321 ex50a << tablepath <<
"201105639_ggF_decay_bb.dat";
2325 ex51 << tablepath <<
"201105639_bb_decay_bb.dat";
2328 ex50b << tablepath <<
"201105639_ggF_decay_WW.dat";
2332 ex52m2 << tablepath <<
"190806463_7a.dat";
2334 ex52m1 << tablepath <<
"190806463_7b.dat";
2338 ex52 << tablepath <<
"CMS-PAS-HIG-16-025.dat";
2342 ex52p1 << tablepath <<
"181011822_7a.dat";
2344 ex52p2 << tablepath <<
"181011822_8a.dat";
2349 ex53 << tablepath <<
"180512191.dat";
2366 ex54 << tablepath <<
"220802717_ggF.dat";
2368 ex55 << tablepath <<
"220802717_bb.dat";
2380 ex56 << tablepath <<
"240509320.dat";
2384 ex57 << tablepath <<
"171203143.dat";
2386 ex58 << tablepath <<
"180401939_a.dat";
2388 ex59 << tablepath <<
"180303838.dat";
2390 ex60 << tablepath <<
"CMS-PAS-HIG-16-023.dat";
2392 ex61 << tablepath <<
"180209407.dat";
2394 ex62 << tablepath <<
"180603548.dat";
2400 ex63 << tablepath <<
"CMS-PAS-B2G-20-004.dat";
2405 ex64 << tablepath <<
"231001643.dat";
2409 ex65 << tablepath <<
"170702909.dat";
2411 ex66 << tablepath <<
"180801365.dat";
2413 ex67 << tablepath <<
"170804188.dat";
2415 ex67p1 << tablepath <<
"220610268.dat";
2417 ex67p2 << tablepath <<
"190404193.dat";
2421 ex67p3 << tablepath <<
"200606391_bblljj.dat";
2424 ex67p4 << tablepath <<
"200606391_bbllnunu.dat";
2428 ex67p5 << tablepath <<
"211203161.dat";
2432 ex67p6 << tablepath <<
"240309430.dat";
2436 ex68 << tablepath <<
"CMS-PAS-HIG-18-005_a.dat";
2438 ex69 << tablepath <<
"180702826_a.dat";
2440 ex70 << tablepath <<
"CMS-PAS-HIG-18-005_b.dat";
2442 ex71 << tablepath <<
"180702826_b.dat";
2445 ex72 << tablepath <<
"14126663.dat";
2447 ex73 << tablepath <<
"151203704.dat";
2449 ex74 << tablepath <<
"150807774_a.dat";
2451 ex75 << tablepath <<
"150807774_b.dat";
2453 ex76 << tablepath <<
"180707915.dat";
2457 ex77 << tablepath <<
"210210076.dat";
2459 ex78 << tablepath <<
"200107763.dat";
2463 ex79 << tablepath <<
"190304560.dat";
2466 lowC01 << tablepath <<
"CMS-HIG-17-029_4e.dat";
2472 lowC02 << tablepath <<
"CMS-HIG-22-007_11d.dat";
2478 lowC03 << tablepath <<
"CMS-HIG-22-007_10.dat";
2481 lowC04 << tablepath <<
"CMS-HIG-19-007_4a.dat";
2484 lowC05 << tablepath <<
"CMS-HIG-19-007_5a.dat";
2487 lowC06 << tablepath <<
"CMS-HIG-21-003_6.dat";
2490 lowC07 << tablepath <<
"CMS-HIG-18-006_10.dat";
2493 lowC08 << tablepath <<
"CMS-PAS-HIG-20-002_5b.dat";
2496 lowC09 << tablepath <<
"CMS-HIG-17-014_4.dat";
2499 lowC10 << tablepath <<
"CMS-EXO-21-018_14d.dat";
2502 lowC11 << tablepath <<
"CMS-EXO-21-018_14f.dat";
2505 lowC12 << tablepath <<
"CMS-HIG-18-026_3c.dat";
2508 lowA01 << tablepath <<
"ATLAS_CERN-EP-2021-157_9.dat";
2511 lowA02 << tablepath <<
"ATLAS_CERN-EP-2021-193_14b.dat";
2514 lowA03 << tablepath <<
"ATLAS_CERN-EP-2021-193_17b.dat";
2517 lowA04 << tablepath <<
"ATLAS_CERN-EP-2018-128_9c.dat";
2520 lowA05 << tablepath <<
"ATLAS_CERN-EP-2020-067_10.dat";
2523 lowA06 << tablepath <<
"ATLAS_CERN-EP-2018-153_6a.dat";
2526 lowA07 << tablepath <<
"ATLAS_CERN-EP-2017-295_2.dat";
2529 lowA08 << tablepath <<
"ATLAS-CONF-2023-035_7b.dat";
2532 lowA09 << tablepath <<
"ATLAS_CERN-EP-2023-070_5a.dat";
2535 lowA10 << tablepath <<
"ATLAS_CERN-EP-2024-235_7.dat";
2538 lowA11 << tablepath <<
"ATLAS_CERN-EP-2023-202_5d.dat";
2541 lowA12 << tablepath <<
"ATLAS_CERN-EP-2024-164_10a.dat";
2544 lowA801 << tablepath <<
"ATLAS_CERN-PH-EP-2015-187_4b.dat";
2547 lowA802 << tablepath <<
"ATLAS_CERN-PH-EP-2015-057_6a.dat";
2550 lowC801 << tablepath <<
"CMS-HIG-16-015_4.dat";
2553 lowC802 << tablepath <<
"CMS-HIG-16-015_5.dat";
2556 lowC803 << tablepath <<
"CMS-HIG-16-015_6.dat";
2559 lowC804 << tablepath <<
"CMS-HIG-17-013_4a.dat";
2562 lowC805 << tablepath <<
"CMS-HIG-14-033_3.dat";
2565 lowC806 << tablepath <<
"CMS-HIG-15-009_5b.dat";
2568 low209a << tablepath <<
"LEP_LHWG-NOTE-2001-08_4.dat";
2571 low209b << tablepath <<
"LEP_CERN-PH-EP-2006-001_t14b.dat";
2574 low209c << tablepath <<
"LEP_CERN-PH-EP-2006-001_t14c.dat";
2577 low209d << tablepath <<
"LEP_CERN-PH-EP-2006-001_t18.dat";
2580 low209e << tablepath <<
"LEP_CERN-PH-EP-2006-001_t19.dat";
2583 lowHpC801 << tablepath <<
"CMS-HIG-13-035_t4.dat";
2586 lowHpC802 << tablepath <<
"CMS-HIG-14-023_t10.dat";
2589 lowHpC803 << tablepath <<
"CMS-HIG-16-030_3.dat";
2592 lowHpC1301 << tablepath <<
"CMS-HIG-18-020_2a.dat";
2595 lowHpC1302 << tablepath <<
"CMS-HIG-18-021_4c.dat";
2598 lowHpA801 << tablepath <<
"ATLAS_CERN-PH-EP-2014-274_7a.dat";
2601 lowHpA1301 << tablepath <<
"ATLAS_CERN-EP-2022-207_8.dat";
2604 lowHpA1302 << tablepath <<
"ATLAS_CERN-EP-2023-070_5b.dat";
2607 lowHpA1303 << tablepath <<
"ATLAS_CERN-EP-2024-185_6.dat";
2610 lowHp209a << tablepath <<
"LEP_CERN-PH-EP-2012-369_4a.dat";
2613 lowHp209b << tablepath <<
"LEP_CERN-PH-EP-2012-369_4d.dat";
2616 lowHp209c << tablepath <<
"LEP_CERN-PH-EP-2008-016_7b.dat";
2619 lowHp172a << tablepath <<
"LEP_CERN-PPE-97-168_3b_tnqq.dat";
2622 lowHp209d << tablepath <<
"LEP_CERN-PH-EP-2008-016_10a.dat";
2625 lowHp209e << tablepath <<
"LEP_CERN-PH-EP-2008-016_10b.dat";
2628 susyHpA01 << tablepath <<
"ATLAS_CERN-EP-2019-191_7b.dat";
2631 susyHpA02 << tablepath <<
"ATLAS_CERN-EP-2022-132_6a.dat";
2634 susyHpC01 << tablepath <<
"CMS-SUS-21-001_5a.dat";
2637 susyHpC02 << tablepath <<
"CMS-SUS-21-008_16a.dat";
2640 thint01 << tablepath <<
"integral_x2_1mx_G_values_log.dat";
2643 thint02 << tablepath <<
"integral_x2_1px_G_values_log.dat";
2646 thint03 << tablepath <<
"integral_x2_G_values_log.dat";
2649 thint04 << tablepath <<
"integral_x_1mx2_G_values_log.dat";
2652 thint05 << tablepath <<
"integral_x_1mx_1px_G_values_log.dat";
2655 thint06 << tablepath <<
"integral_x2_1mx_G_variable_set_1_values_log.dat";
2658 thint07 << tablepath <<
"integral_x2_G_variable_set_1_values_log.dat";
2661 thint08 << tablepath <<
"integral_x_1mx2_G_variable_set_0_values_log.dat";
2664 thint09 << tablepath <<
"integral_perturb_ATHDM_values_log.dat";
2667 thint10 << tablepath <<
"integral_Hp_A_Wstar.dat";
2672 bsg1 << tablepath <<
"bsgammatable.dat";
2680 double params[] = {mass};
2696 double params[] = {mass};
2712 double params[] = {mass};
2728 double params[] = {mass};
2744 double params[] = {mass};
2760 double params[] = {mass};
2776 double params[] = {mass};
2792 double params[] = {mass};
2808 double params[] = {mass};
2814 double newResult = 0.0;
2815 if (mass>=10. && mass <=2000.) {
2827 double params[] = {mass};
2833 double newResult = 0.0;
2834 if (mass>=10. && mass <=2000.) {
2846 double params[] = {mass};
2852 double newResult = 0.0;
2853 if (mass>=10. && mass <=2000.) {
2865 double params[] = {mass};
2871 double newResult = 0.0;
2872 if (mass>=10. && mass <=2000.) {
2884 double params[] = {mass};
2890 double newResult = 0.0;
2891 if (mass>=10. && mass <=2000.) {
2903 double params[] = {mass};
2909 double newResult = 0.0;
2910 if (mass>=10. && mass <=2000.) {
2922 double params[] = {mass};
2928 double newResult = 0.0;
2929 if (mass>=10. && mass <=2000.) {
2941 double params[] = {mass};
2947 double newResult = 0.0;
2948 if (mass>=10. && mass <=2000.) {
2960 double params[] = {mass};
2966 double newResult = 0.0;
2967 if (mass>=10. && mass <=2000.) {
2979 double params[] = {mass};
2985 double newResult = 0.0;
2986 if (mass>=10. && mass <=2000.) {
2998 double params[] = {mass};
3004 double newResult = 0.0;
3005 if (mass>=10. && mass <=2000.) {
3017 double params[] = {mass};
3023 double newResult = 0.0;
3024 if (mass>=10. && mass <=2000.) {
3036 double params[] = {mass};
3042 double newResult = 0.0;
3043 if (mass>=20. && mass <=2000.) {
3055 double params[] = {mass};
3061 double newResult = 0.0;
3062 if (mass>=20. && mass <=2000.) {
3074 double params[] = {mass};
3080 double newResult = 0.0;
3081 if (mass>=20. && mass <=2000.) {
3093 double params[] = {mass};
3099 double newResult = 0.0;
3100 if (mass>=20. && mass <=2000.) {
3112 double params[] = {mass};
3118 double newResult = 0.0;
3119 if (mass>=20. && mass <=2000.) {
3131 double params[] = {mass};
3137 double newResult = 0.0;
3138 if (mass>=20. && mass <=2000.) {
3150 double params[] = {
mHp, logtb};
3156 double newResult = 0.0;
3157 if (
mHp>=180. && mHp <=1400. && logtb>=-1. && logtb<=1.75) {
3169 double params[] = {
mHp, logtb};
3175 double newResult = 0.0;
3176 if (
mHp>=180. && mHp <=2000. && logtb>=-1. && logtb<=1.75) {
3189 double params[] = {mass};
3205 double params[] = {mass};
3221 double params[] = {mass};
3237 double params[] = {mass};
3253 double params[] = {mass};
3269 double params[] = {mass};
3285 double params[] = {mass};
3301 double params[] = {mass};
3316 double params[] = {mass};
3331 double params[] = {mass};
3347 double params[] = {mass};
3364 double params[] = {mass};
3378 double params[] = {mass};
3392 double params[] = {mass};
3406 double params[] = {mass};
3422 double params[] = {mass};
3437 double params[] = {mass};
3457 double params[] = {mass};
3478 double params[] = {mass};
3495 double params[] = {mass};
3526 double params[] = {mass};
3541 double params[] = {mass};
3556 double params[] = {mass};
3572 double params[] = {mass};
3587 double params[] = {mass};
3603 double params[] = {mass};
3631 double params[] = {mass};
3645 double params[] = {mass};
3659 double params[] = {mass};
3673 double params[] = {mass};
3687 double params[] = {mass};
3701 double params[] = {mass};
3715 double params[] = {mass};
3729 double params[] = {mass};
3743 double params[] = {mass};
3757 double params[] = {mass};
3771 double params[] = {mass};
3785 double params[] = {mass};
3799 double params[] = {mass};
3813 double params[] = {mass};
3827 double params[] = {mass};
3841 double params[] = {mass};
3855 double params[] = {mass};
3869 double params[] = {mass};
3883 double params[] = {mass};
3897 double params[] = {mass};
3911 double params[] = {mass};
3925 double params[] = {mass};
3939 double params[] = {mass};
3953 double params[] = {mass};
3967 double params[] = {mass};
3981 double params[] = {mass};
3995 double params[] = {mass};
4014 double params[] = {mass};
4032 double params[] = {mass};
4056 double params[] = {mass};
4074 double params[] = {mass};
4105 double params[] = {mass};
4119 double params[] = {mass};
4133 double params[] = {mass};
4147 double params[] = {mass};
4161 double params[] = {mass};
4175 double params[] = {mass};
4189 double params[] = {mass};
4203 double params[] = {mass};
4223 double params[] = {mass};
4238 double params[] = {mass};
4259 double params[] = {mass};
4273 double params[] = {mass};
4287 double params[] = {mass};
4301 double params[] = {mass};
4315 double params[] = {mass};
4329 double params[] = {mass};
4343 double params[] = {mass};
4357 double params[] = {mass};
4371 double params[] = {mass};
4385 double params[] = {mass};
4418 double params[] = {mass};
4433 double params[] = {mass};
4447 double params[] = {mass};
4461 double params[] = {mass};
4475 double params[] = {mass};
4491 double params[] = {mass};
4509 double params[] = {mass};
4539 double params[] = {mass};
4557 double params[] = {mass};
4575 double params[] = {mass};
4591 double params[] = {mass};
4608 double params[] = {mass};
4629 double params[] = {mass};
4643 double params[] = {mass};
4657 double params[] = {mass};
4671 double params[] = {mass};
4685 double params[] = {mass};
4699 double params[] = {mass};
4713 double params[] = {mass};
4727 double params[] = {mass};
4741 double params[] = {mass};
4758 double params[] = {mass};
4778 double params[] = {mass};
4792 double params[] = {mass};
4806 double params[] = {mi, mj};
4822 double params[] = {mi, mj};
4837 double params[] = {mi, mj};
4853 double params[] = {mi, mj};
4869 double params[] = {mj, mi};
4885 double params[] = {mj, mi};
4902 double params[] = {mj, mi};
4919 double params[] = {mass};
4935 double params[] = {mass};
4951 double params[] = {mass};
4967 double params[] = {mass};
4983 double params[] = {mass};
4999 double params[] = {mass};
5013 double params[] = {mass};
5029 double params[] = {mass};
5044 double params[] = {mass};
5058 double params[] = {mass};
5072 double params[] = {mass};
5086 double params[] = {mass};
5100 double params[] = {mass};
5114 double params[] = {mass};
5128 double params[] = {mass};
5142 double params[] = {mass};
5156 double params[] = {mass};
5170 double params[] = {mass};
5184 double params[] = {mass};
5197 double params[] = {mass};
5210 double params[] = {mass};
5223 double params[] = {mass};
5236 double params[] = {mass};
5249 double params[] = {mass};
5262 double params[] = {mass};
5275 double params[] = {mass};
5288 double params[] = {mass};
5301 double params[] = {mass};
5314 double params[] = {mass};
5327 double params[] = {mass};
5340 double params[] = {mass};
5353 double params[] = {mass};
5366 double params[] = {mass};
5380 double params[] = {mass};
5393 double params[] = {mass};
5406 double params[] = {mass};
5419 double params[] = {mass};
5432 double params[] = {mass};
5445 double params[] = {mass};
5458 double params[] = {mass};
5471 double params[] = {mass};
5484 double params[] = {mass};
5497 double params[] = {mass};
5510 double params[] = {mi, mj};
5523 double params[] = {mi, mj};
5536 double params[] = {mass};
5549 double params[] = {mass};
5562 double params[] = {mass};
5575 double params[] = {mass};
5588 double params[] = {mass};
5601 double params[] = {mass};
5614 double params[] = {mass};
5627 double params[] = {mass};
5640 double params[] = {mass};
5653 double params[] = {mass};
5666 double params[] = {mass};
5679 double params[] = {mass};
5692 double params[] = {mass};
5705 double params[] = {mi, mj};
5718 double params[] = {mi, mj};
5731 double params[] = {mass};
5744 double params[] = {mass};
5757 double params[] = {mass};
5770 double params[] = {mass};
5787 if(mass >= 45. && mass <= 65.)
5789 else if(mass > 65. && mass <= 76.)
5791 else if(mass > 76. && mass <= 86.)
5793 else if(mass > 86. && mass <= 89.)
5795 else if(mass > 89. && mass <= 96.)
5797 else if(mass > 96. && mass <= 97.)
5799 else if(mass > 97. && mass <= 98.)
5801 else if(mass > 98. && mass <= 99.)
5803 else if(mass > 99. && mass <= 101.)
5805 else if(mass > 101. && mass <= 102.)
5818 double params[] = {wa, wb};
5823 double log10wa = log10(wa);
5824 double log10wb = log10(wb);
5834 double params[] = {wa, wb};
5839 double log10wa = log10(wa);
5840 double log10wb = log10(wb);
5851 double params[] = {wa, wb};
5856 double log10wa = log10(wa);
5857 double log10wb = log10(wb);
5866 double params[] = {wa, wb};
5871 double log10wa = log10(wa);
5872 double log10wb = log10(wb);
5882 double params[] = {wa, wb};
5887 double log10wa = log10(wa);
5888 double log10wb = log10(wb);
5898 double params[] = {wb};
5903 double log10wb = log10(wb);
5913 double params[] = {wb};
5918 double log10wb = log10(wb);
5928 double params[] = {wb};
5933 double log10wb = log10(wb);
5943 double params[] = {ratio};
5948 double log10ratio = log10(ratio);
5958 double params[] = {mi, mj};
5973 double params[] = {logtb, logmHp};
5992 gslpp::complex z = -gslpp::complex::i()*M_PI;
5993 return -pow(log((1.0+sqrt(1.0-x))/(1.0-sqrt(1.0-x)))+z,2)/4.0;
5996 return pow(asin(sqrt(1.0/x)),2);
6002 gslpp::complex z = -gslpp::complex::i()*M_PI;
6003 gslpp::complex gs1 = sqrt(1.0-x)*(log((1.0+sqrt(1.0-x))/(1.0-sqrt(1.0-x)))+z)/2.0;
6007 gslpp::complex gg1 = sqrt(x-1.0)*asin(sqrt(1.0/x));
6013 return tau*lambda/(tau-lambda)/2.0+tau*tau*lambda*lambda/((tau-lambda)
6014 *(tau-lambda))/2.0*(
f_func(tau)-
f_func(lambda))+tau*tau*lambda/((tau-lambda)
6019 return -tau*lambda/(tau-lambda)/2.0*(
f_func(tau)-
f_func(lambda));
6027 double params[] = {mHl2,
Mu,
Mc,
Mt};
6033 double TAUu=4.0*
Mu*
Mu/mHl2;
6034 double TAUc=4.0*
Mc*
Mc/mHl2;
6035 double TAUt=4.0*
Mt*
Mt/mHl2;
6036 gslpp::complex newResult = -(8./3.)*(TAUu*(1.0+(1.0-TAUu)*
f_func(TAUu))
6037 +TAUc*(1.0+(1.0-TAUc)*
f_func(TAUc))+TAUt*(1.0+(1.0-TAUt)*
f_func(TAUt)));
6045 double params[] = {mHh2,
Mc,
Mt};
6051 double TAUc=4.0*
Mc*
Mc/mHh2;
6052 double TAUt=4.0*
Mt*
Mt/mHh2;
6053 gslpp::complex newResult = -(8./3.)*(TAUc*(1.0+(1.0-TAUc)*
f_func(TAUc))
6054 +TAUt*(1.0+(1.0-TAUt)*
f_func(TAUt)));
6062 double params[] = {mA2,
Mc,
Mt};
6068 double TAUc=4.0*
Mc*
Mc/mA2;
6069 double TAUt=4.0*
Mt*
Mt/mA2;
6070 gslpp::complex newResult = -(8./3.)*(TAUc*
f_func(TAUc)+TAUt*
f_func(TAUt));
6078 double params[] = {mHl2,
Md,
Ms,
Mb};
6084 double TAUd=4.0*
Md*
Md/mHl2;
6085 double TAUs=4.0*
Ms*
Ms/mHl2;
6086 double TAUb=4.0*
Mb*
Mb/mHl2;
6087 gslpp::complex newResult = -(2./3.)*(TAUd*(1.0+(1.0-TAUd)*
f_func(TAUd))
6088 +TAUs*(1.0+(1.0-TAUs)*
f_func(TAUs))+TAUb*(1.0+(1.0-TAUb)*
f_func(TAUb)));
6096 double params[] = {mHh2,
Ms,
Mb};
6102 double TAUs=4.0*
Ms*
Ms/mHh2;
6103 double TAUb=4.0*
Mb*
Mb/mHh2;
6104 gslpp::complex newResult = -(2./3.)*(TAUs*(1.0+(1.0-TAUs)*
f_func(TAUs))
6105 +TAUb*(1.0+(1.0-TAUb)*
f_func(TAUb)));
6113 double params[] = {mA2,
Ms,
Mb};
6119 double TAUs=4.0*
Ms*
Ms/mA2;
6120 double TAUb=4.0*
Mb*
Mb/mA2;
6121 gslpp::complex newResult = -(2./3.)*(TAUs*
f_func(TAUs)+TAUb*
f_func(TAUb));
6129 double params[] = {mHl2,
Me,
Mmu,
Mtau};
6135 double TAUe=4.0*
Me*
Me/mHl2;
6136 double TAUmu=4.0*
Mmu*
Mmu/mHl2;
6138 gslpp::complex newResult = -2.0*(TAUe*(1.0+(1.0-TAUe)*
f_func(TAUe))
6139 +TAUmu*(1.0+(1.0-TAUmu)*
f_func(TAUmu))
6140 +TAUtau*(1.0+(1.0-TAUtau)*
f_func(TAUtau)));
6148 double params[] = {mHh2,
Mmu,
Mtau};
6154 double TAUmu=4.0*
Mmu*
Mmu/mHh2;
6156 gslpp::complex newResult = -2.0*(TAUmu*(1.0+(1.0-TAUmu)*
f_func(TAUmu))+
6157 TAUtau*(1.0+(1.0-TAUtau)*
f_func(TAUtau)));
6165 double params[] = {mA2,
Mmu,
Mtau};
6171 double TAUmu=4.0*
Mmu*
Mmu/mA2;
6173 gslpp::complex newResult = -2.0*(TAUmu*
f_func(TAUmu)+TAUtau*
f_func(TAUtau));
6181 double params[] = {mH,
MW};
6187 double TAUw=4.0*
MW*
MW/(mH*mH);
6188 gslpp::complex newResult = 2.0 + 3.0*TAUw + 3.0*TAUw*(2.0-TAUw)*
f_func(TAUw);
6196 double params[] = {
mHp2, mH};
6202 double TAUhp=4.0*
mHp2/(mH*mH);
6203 gslpp::complex newResult = -TAUhp*(1.0-TAUhp*
f_func(TAUhp));
6209gslpp::complex
GeneralTHDMcache::A_h_U(
const double mHl2,
const double cW2,
const double Mu,
const double Mc,
const double Mt,
const double MZ)
const {
6217 double TAUu=4.0*
Mu*
Mu/mHl2;
6218 double TAUc=4.0*
Mc*
Mc/mHl2;
6219 double TAUt=4.0*
Mt*
Mt/mHl2;
6224 gslpp::complex newResult = -4.0*(1.0/2.0-4.0/3.0*
sW2)*(
Int1(TAUu,LAMu)+
Int1(TAUc,LAMc)
6234 double params[] = {mHh2,
cW2,
Mc,
Mt,
MZ};
6240 double TAUc=4.0*
Mc*
Mc/mHh2;
6241 double TAUt=4.0*
Mt*
Mt/mHh2;
6252 gslpp::complex newResult = -4.0*(1.0/2.0-4.0/3.0*
sW2)*(
Int1(TAUc,LAMc)-
Int2(TAUc,LAMc)
6259gslpp::complex
GeneralTHDMcache::A_A_U(
const double mA2,
const double cW2,
const double Mc,
const double Mt,
const double MZ)
const {
6261 double params[] = {mA2,
cW2,
Mc,
Mt,
MZ};
6267 double TAUc=4.0*
Mc*
Mc/mA2;
6268 double TAUt=4.0*
Mt*
Mt/mA2;
6272 gslpp::complex newResult = -4.0*(1.0/2.0-4.0/3.0*
sW2)*(-
Int2(TAUc,LAMc)-
Int2(TAUt,LAMt))/sqrt(
sW2*
cW2);
6278gslpp::complex
GeneralTHDMcache::A_h_D(
const double mHl2,
const double cW2,
const double Md,
const double Ms,
const double Mb,
const double MZ)
const {
6286 double TAUd=4.0*
Md*
Md/mHl2;
6287 double TAUs=4.0*
Ms*
Ms/mHl2;
6288 double TAUb=4.0*
Mb*
Mb/mHl2;
6293 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0/3.0*
sW2)*(
Int1(TAUd,LAMd)+
Int1(TAUs,LAMs)
6302 double params[] = {mHh2,
cW2,
Ms,
Mb,
MZ};
6308 double TAUs=4.0*
Ms*
Ms/mHh2;
6309 double TAUb=4.0*
Mb*
Mb/mHh2;
6313 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0/3.0*
sW2)*(
Int1(TAUs,LAMs)-
Int2(TAUs,LAMs)
6320gslpp::complex
GeneralTHDMcache::A_A_D(
const double mA2,
const double cW2,
const double Ms,
const double Mb,
const double MZ)
const {
6322 double params[] = {mA2,
cW2,
Ms,
Mb,
MZ};
6328 double TAUs=4.0*
Ms*
Ms/mA2;
6329 double TAUb=4.0*
Mb*
Mb/mA2;
6333 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0/3.0*
sW2)*(-
Int2(TAUs,LAMs)-
Int2(TAUb,LAMb))/sqrt(
sW2*
cW2);
6339gslpp::complex
GeneralTHDMcache::A_h_L(
const double mHl2,
const double cW2,
const double Me,
const double Mmu,
const double Mtau,
const double MZ)
const {
6347 double TAUe=4.0*
Me*
Me/mHl2;
6348 double TAUmu=4.0*
Mmu*
Mmu/mHl2;
6354 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0*
sW2)*(
Int1(TAUe,LAMe)+
Int1(TAUmu,LAMmu)
6356 -
Int2(TAUtau,LAMtau));
6362gslpp::complex
GeneralTHDMcache::A_HH_L(
const double mHh2,
const double cW2,
const double Mmu,
const double Mtau,
const double MZ)
const {
6370 double TAUmu=4.0*
Mmu*
Mmu/mHh2;
6375 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0*
sW2)*(
Int1(TAUmu,LAMmu)-
Int2(TAUmu,LAMmu)
6376 +
Int1(TAUtau,LAMtau)-
Int2(TAUtau,LAMtau));
6384gslpp::complex
GeneralTHDMcache::A_A_L(
const double mA2,
const double cW2,
const double Mmu,
const double Mtau,
const double MZ)
const {
6392 double TAUmu=4.0*
Mmu*
Mmu/mA2;
6397 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0*
sW2)*(-
Int2(TAUmu,LAMmu)-
Int2(TAUtau,LAMtau))/sqrt(
sW2*
cW2);
6405 double params[] = {mH,
cW2,
MW,
MZ};
6411 double TAUw=4.0*
MW*
MW/(mH*mH);
6414 gslpp::complex newResult = -sqrt(
cW2/
sW2)*(4.0*(3.0-
sW2/
cW2)*
Int2(TAUw,LAMw)
6415 +((1.0+2.0/TAUw)*
sW2/
cW2-(5.0+2.0/TAUw))*
Int1(TAUw,LAMw));
6429 double TAUhp=4.0*
mHp2/(mH*mH);
6432 gslpp::complex newResult = (1.0-2.0*
sW2)/sqrt(
cW2*
sW2)*
Int1(TAUhp,LAMhp);
6441 double s_LEP = sqrts * sqrts;
6442 double one_MZ2_s = 1. -
MZ*
MZ / s_LEP;
6444 double cWsW = sqrt(
cW2) * sqrt(
sW2);
6445 double ahat_e = -1. / 4. / cWsW;
6446 double vhat_e = (-1. + 4. *
sW2) / 4. / cWsW;
6447 double vhat_H = (-1. + 2. *
sW2) / 2. / cWsW;
6448 double beta_H = sqrt(std::fabs(1. - 4. * mass2 / s_LEP));
6450 return (2. *
GF*
GF *
MW*
MW*
MW*
MW *
sW2*
sW2 * (1. + 2. * vhat_e * vhat_H / one_MZ2_s +
6451 (ahat_e*ahat_e + vhat_e*vhat_e) * vhat_H*vhat_H / one_MZ2_s / one_MZ2_s) *
6452 beta_H*beta_H*beta_H / (3. * M_PI * s_LEP));
6464 double params[] = {a2, b2, c2};
6471 double newResult = 0.0;
6472 double x = (a2-b2-c2)*(a2-b2-c2)-4.0*b2*c2;
6473 if(x>0) newResult = sqrt(std::fabs(x/a2))/2.0;
6496 return sqrt(std::fabs(1.0 - 4.0*mf*mf / m_2));
6501 double GeneralTHDMcache::lambdaijk(
const double Ri1,
const double Ri2,
const double Ri3,
const double Rj1,
const double Rj2,
const double Rj3,
const double Rk1,
const double Rk2,
const double Rk3,
const double lambda1H,
const double lambda3H,
const double lambda4H,
const double Relambda5H,
const double Imlambda5H,
const double Relambda6H,
const double Imlambda6H,
const double Relambda7H,
const double Imlambda7H)
const
6503 return (1.0/2.0)*
vev*(Imlambda7H*(-Ri3*Rj3*Rk3 - Ri2*Rj2*Rk3) - 3.0*Imlambda6H*Ri1*Rj1*Rk3
6504 + lambda1H*Ri1*Rj1*Rk1 + Relambda7H*Ri2*Rj2*Rk2 + 3.0*Relambda6H*Ri1*Rj1*Rk2
6505 +(Relambda5H + lambda3H + lambda4H)*Ri1*Rj2*Rk2 - (2.0*Relambda5H - lambda3H - lambda4H)*Ri1*Rj3*Rk3
6506 + Relambda7H*Ri2*Rj3*Rk3 - Imlambda5H*Ri1*Rj2*Rk3);
6519 throw std::runtime_error(
"Signal strengths are only available in the A2HDM.");
6523 gslpp::complex i = gslpp::complex::i();
6544 double SigmaggF8 =
myGTHDM->computeSigmaggH(8.0);
6546 double Sigmaggh_tt8 =
myGTHDM->computeSigmaggH_tt(8.0);
6550 double Sigmaggh_bb8 =
myGTHDM->computeSigmaggH_bb(8.0);
6554 double Sigmatth8 =
myGTHDM->computeSigmattH(8.0);
6569 double Sigmatth13 =
myGTHDM->computeSigmattH(13.0);
6573 double SigmaVBFVh13 = (
myGTHDM->computeSigmaVBF(13.0)+
myGTHDM->computeSigmaWH(13.0)+
myGTHDM->computeSigmaZH(13.0));
6610 rh_ggE = yu1.real()*
yd1.real() + (yu1.real()*yu1.real() - yu1.real()*
yd1.real())*(Sigmaggh_tt8/SigmaggF8) + (
yd1.real()*
yd1.real() - yu1.real()*
yd1.real())*(Sigmaggh_bb8/SigmaggF8);
6611 rh_ggO = yu1.imag()*yu1.imag() + (yu1.imag()*yu1.imag() - yu1.imag()*
yd1.imag())*rSigmagghO_b8 + (
yd1.imag()*
yd1.imag() - yu1.imag()*
yd1.imag())*rSigmagghO_b8;
6624 gslpp::complex I_h_W=
Rij_GTHDM(0,0)*I_hSM_W;
6627 gslpp::complex I_hSM_F= fermU+fermD+fermL;
6628 gslpp::complex I_hE_F= yu1.real()*fermU+
yd1.real()*fermD+yl1.real()*fermL;
6644 gslpp::complex I_hO_F = yu1.imag()*I_h_Ux +
yd1.imag()*I_h_Dx + yl1.imag()*I_h_Lx;
6647 rh_gaga = ((I_hE_F+I_h_W+I_h_Hp).abs2()+ (I_hO_F).abs2())/(I_hSM_F +I_hSM_W).abs2();
6669 gslpp::complex A_h_W =
Rij_GTHDM(0,0)*A_hSM_W;
6671 gslpp::complex A_hSM_F = (A_hE_Ux+ A_hE_Dx+ A_hE_Lx)/sqrt(
sW2*
cW2);
6672 gslpp::complex A_hE_F = (yu1.real()*A_hE_Ux+
yd1.real()*A_hE_Dx+ yl1.real()*A_hE_Lx)/sqrt(
sW2*
cW2);
6682 gslpp::complex A_hO_F=yu1.imag()*A_hO_Ux +
yd1.imag()*A_hO_Dx + yl1.imag()*A_hO_Lx;
6685 rh_Zga = ((A_hE_F+A_h_W+A_h_Hp).abs2()+ A_hO_F.abs2())/(A_hSM_F +A_hSM_W ).abs2();
6710 double Gamma_hHH =
HSTheta(
m1 - 2.0*
m2)*sqrt(std::fabs(1 - 4*(
m2*
m2)/(
m1*
m1)))*(lambda122*lambda122)/(32.0*M_PI*
m1);
6718 double Gamma_hAA =
HSTheta(
m1 - 2.0*
m3)*sqrt(std::fabs(1 - 4*(
m3*
m3)/(
m1*
m1)))*(lambda133*lambda133)/(32.0*M_PI*
m1);
6720 double Gamma_hHpHm =
HSTheta(
m1-2.0*
mHp)*sqrt(std::fabs(1 - 4*(
mHp*
mHp)/(
m1*
m1)))*(lambda1pm*lambda1pm)/(16.0*M_PI*
m1);
6762 double Gamma_h_exp =
myGTHDM->computeGammaHTotal();
6776 Gamma_h_inv = Gamma_hHpHm + Gamma_hHH + Gamma_hAA + Gamma_hHZ + Gamma_hAZ;
6778 Gamma_h = (Gamma_htobb + Gamma_htoWW + Gamma_htoZZ + Gamma_htotautau + Gamma_htogaga +
6779 Gamma_htogg + Gamma_htoZga + Gamma_htocc +
Gamma_h_inv);
6804 throw std::runtime_error(
"Direct Searches are only aviable in the A2HDM.");
6808 gslpp::complex i = gslpp::complex::i();
6843 double rphi2_QuQuE= yu2.real()*yu2.real();
6844 double rphi2_QuQuO= yu2.imag()*yu2.imag();
6845 double rphi2_QdQdE=
yd2.real()*
yd2.real();
6846 double rphi2_QdQdO=
yd2.imag()*
yd2.imag();
6847 double rphi2_QlQlE= yl2.real()*yl2.real();
6848 double rphi2_QlQlO= yl2.imag()*yl2.imag();
6849 rphi2_ggE = yu2.real()*
yd2.real() + (yu2.real()*yu2.real() - yu2.real()*
yd2.real())*rSigmaggphi2E_tc8 + (
yd2.real()*
yd2.real() - yu2.real()*
yd2.real())*rSigmaggphi2E_b8;
6850 rphi2_ggO = yu2.imag()*yu2.imag() + (yu2.imag()*yu2.imag() - yu2.imag()*
yd2.imag())*rSigmaggphi2O_tc8 + (
yd2.imag()*
yd2.imag() - yu2.imag()*
yd2.imag())*rSigmaggphi2O_b8;
6864 gslpp::complex I_phi2E_F= yu2.real()*I_HH2_Ux+
yd2.real()*I_HH2_Dx+yl2.real()*I_HH2_Lx;
6880 gslpp::complex I_phi2O_F = yu2.imag()*I_A2_Ux +
yd2.imag()*I_A2_Dx + yl2.imag()*I_A2_Lx;
6882 double Gamma_phi2gaga=(
GF*
Ale*
Ale*
m2*
m2*
m2/(sqrt(2.0)*128.0*M_PI*M_PI*M_PI))*((I_phi2E_F+I_phi2_W+I_phi2_Hp).abs2()
6883 + (I_phi2O_F).abs2());
6894 gslpp::complex A_phi2E_F = (yu2.real()*A_HH2_Ux+
yd2.real()*A_HH2_Dx+ yl2.real()*A_HH2_Lx)/sqrt(
sW2*
cW2);
6908 gslpp::complex A_phi2O_F=yu2.imag()*A_A2_Ux +
yd2.imag()*A_A2_Dx + yl2.imag()*A_A2_Lx;
6913 + A_phi2O_F.abs2());
6923 double Gamma_phi2gg=(
rphi2_ggE)*
GF*
Als*
Als*
m2*
m2*
m2/(sqrt(2.0)*16.0*M_PI*M_PI*M_PI)*(9.0/4.0)*(I_HH2_Ux/4.0+I_HH2_Dx).abs2()
6947 if (
m2>=20. &&
m2 <=2000.) {
7002double lambda132 =
lambdaijk(
R11,
R12,
R13,
R31,
R32,
R33,
R21,
R22,
R23,
lambda1,
lambda3,
lambda4,
Relambda5,
7004 lambdaijk(
R21,
R22,
R23,
R11,
R12,
R13,
R31,
R32,
R33,
lambda1,
lambda3,
lambda4,
Relambda5,
7006 lambdaijk(
R21,
R22,
R23,
R31,
R32,
R33,
R11,
R12,
R13,
lambda1,
lambda3,
lambda4,
Relambda5,
7008 lambdaijk(
R11,
R12,
R13,
R21,
R22,
R23,
R32,
R32,
R32,
lambda1,
lambda3,
lambda4,
Relambda5,
7010 lambdaijk(
R31,
R32,
R33,
R11,
R12,
R13,
R21,
R22,
R23,
lambda1,
lambda3,
lambda4,
Relambda5,
7012 lambdaijk(
R31,
R32,
R33,
R21,
R22,
R23,
R11,
R12,
R13,
lambda1,
lambda3,
lambda4,
Relambda5,
7015double lambda332 =
lambdaijk(
R31,
R32,
R33,
R31,
R32,
R33,
R21,
R22,
R23,
lambda1,
lambda3,
lambda4,
Relambda5,
7017 lambdaijk(
R21,
R22,
R23,
R31,
R32,
R33,
R31,
R32,
R33,
lambda1,
lambda3,
lambda4,
Relambda5,
7019 lambdaijk(
R31,
R32,
R33,
R21,
R22,
R23,
R31,
R32,
R33,
lambda1,
lambda3,
lambda4,
Relambda5,
7024double Gammaphi2_phi3phi3=
HSTheta(
m2 - 2.0*
m3)*sqrt(std::fabs(1.0 - (4.0*
m3_2)/
m2_2))*lambda332*lambda332/(32.0*
m2*M_PI);
7027double lambda112 = (2.0)*(
lambdaijk(
R11,
R12,
R13,
R21,
R22,
R23,
R11,
R12,
R13,
lambda1,
lambda3,
lambda4,
Relambda5,
7029 lambdaijk(
R11,
R12,
R13,
R11,
R12,
R13,
R21,
R22,
R23,
lambda1,
lambda3,
lambda4,
Relambda5,
7031 lambdaijk(
R21,
R22,
R23,
R11,
R12,
R13,
R11,
R12,
R13,
lambda1,
lambda3,
lambda4,
Relambda5,
7035double Gammaphi2_phi1phi1=
HSTheta(
m2 - 2.0*
m1)*sqrt(std::fabs(1.0 - (4.0*
m1_2)/
m2_2))*lambda112*lambda112/(32.0*
m2*M_PI);
7052double Gammaphi2_HpHm=
HSTheta(
m2 - 2.0*sqrt(
mHp2))*sqrt(std::fabs(1.0 - (4.0*
mHp2)/
m2_2))*lambdaphi2HpHm*lambdaphi2HpHm/(8.0*
m2*M_PI);
7083 +Gammaphi2_bb+Gammaphi2_tautau+Gammaphi2_mumu
7084 +Gammaphi2_WW+Gammaphi2_ZZ+Gamma_phi2gaga
7085 +Gamma_phi2Zga+Gamma_phi2gg +Gammaphi2_phi1phi3
7086 +Gammaphi2_phi1phi1+Gammaphi2_phi3phi3+Gammaphi2_HpHm
7087 +Gammaphi2_phi1Z+Gammaphi2_phi3Z+Gammaphi2_HpW;
7151 throw std::runtime_error(
"Direct Searches are only aviable in the A2HDM.");
7155 gslpp::complex i = gslpp::complex::i();
7192 double rphi3_QuQuE= yu3.real()*yu3.real();
7193 double rphi3_QuQuO= yu3.imag()*yu3.imag();
7194 double rphi3_QdQdE=
yd3.real()*
yd3.real();
7195 double rphi3_QdQdO=
yd3.imag()*
yd3.imag();
7196 double rphi3_QlQlE= yl3.real()*yl3.real();
7197 double rphi3_QlQlO= yl3.imag()*yl3.imag();
7198 rphi3_ggE = yu3.real()*
yd3.real() + (yu3.real()*yu3.real() - yu3.real()*
yd3.real())*rSigmaggphi3E_tc8 + (
yd3.real()*
yd3.real() - yu3.real()*
yd3.real())*rSigmaggphi3E_b8;
7199 rphi3_ggO = yu3.imag()*
yd3.imag() + (yu3.imag()*yu3.imag() - yu3.imag()*
yd3.imag())*rSigmaggphi3O_tc8 + (
yd3.imag()*
yd3.imag() - yu3.imag()*
yd3.imag())*rSigmaggphi3O_b8;
7211 gslpp::complex I_phi3E_F= yu3.real()*I_HH_Ux+
yd3.real()*I_HH_Dx+yl3.real()*I_HH_Lx;
7227 gslpp::complex I_phi3O_F = yu3.imag()*I_A_Ux +
yd3.imag()*I_A_Dx + yl3.imag()*I_A_Lx;
7229 double Gamma_phi3gaga=(
GF*
Ale*
Ale*
m3*
m3*
m3/(sqrt(2.0)*128.0*M_PI*M_PI*M_PI))*((I_phi3E_F+I_phi3_W+I_phi3_Hp).abs2() + (I_phi3O_F).abs2());
7239 gslpp::complex A_phi3E_F = (yu3.real()*A_HH_Ux+
yd3.real()*A_HH_Dx+ yl3.real()*A_HH_Lx)/sqrt(
sW2*
cW2);
7250 gslpp::complex A_phi3O_F=yu3.imag()*A_A_Ux +
yd3.imag()*A_A_Dx + yl3.imag()*A_A_Lx;
7253 + A_phi3O_F.abs2());
7258 double Gamma_phi3gg=(
rphi3_ggE)*
GF*
Als*
Als*
m3*
m3*
m3/(sqrt(2.0)*16.0*M_PI*M_PI*M_PI)*(9.0/4.0)*(I_HH_Ux/4.0+I_HH_Dx).abs2()
7282 if (
m3>=20. &&
m3 <=2000.) {
7318double lambda123 =
lambdaijk(
R11,
R12,
R13,
R21,
R22,
R23,
R31,
R32,
R33,
lambda1,
lambda3,
lambda4,
Relambda5,
7320 lambdaijk(
R21,
R22,
R23,
R11,
R12,
R13,
R31,
R32,
R33,
lambda1,
lambda3,
lambda4,
Relambda5,
7322 lambdaijk(
R21,
R22,
R23,
R31,
R32,
R33,
R11,
R12,
R13,
lambda1,
lambda3,
lambda4,
Relambda5,
7324 lambdaijk(
R11,
R12,
R13,
R31,
R32,
R33,
R21,
R22,
R23,
lambda1,
lambda3,
lambda4,
Relambda5,
7326 lambdaijk(
R31,
R32,
R33,
R11,
R12,
R13,
R21,
R22,
R23,
lambda1,
lambda3,
lambda4,
Relambda5,
7328 lambdaijk(
R31,
R32,
R33,
R21,
R22,
R23,
R11,
R12,
R13,
lambda1,
lambda3,
lambda4,
Relambda5,
7330double lambda223 =
lambdaijk(
R21,
R22,
R23,
R21,
R22,
R23,
R31,
R32,
R33,
lambda1,
lambda3,
lambda4,
Relambda5,
7332 lambdaijk(
R21,
R22,
R23,
R31,
R32,
R33,
R21,
R22,
R23,
lambda1,
lambda3,
lambda4,
Relambda5,
7334 lambdaijk(
R31,
R32,
R33,
R21,
R22,
R23,
R21,
R22,
R23,
lambda1,
lambda3,
lambda4,
Relambda5,
7336double lambda113 =
lambdaijk(
R11,
R12,
R13,
R11,
R13,
R13,
R31,
R32,
R33,
lambda1,
lambda3,
lambda4,
Relambda5,
7338 lambdaijk(
R11,
R11,
R13,
R31,
R32,
R33,
R11,
R12,
R13,
lambda1,
lambda3,
lambda4,
Relambda5,
7340 lambdaijk(
R31,
R32,
R33,
R11,
R12,
R13,
R11,
R12,
R13,
lambda1,
lambda3,
lambda4,
Relambda5,
7346double Gammaphi3_phi2phi2=
HSTheta(
m3 - 2.0*
m2)*sqrt(std::fabs(1.0 - (4.0*
m2_2)/
m3_2))*lambda223*lambda223/(32.0*
m3*M_PI);
7348double Gammaphi3_phi1phi1=
HSTheta(
m3 - 2.0*
m1)*sqrt(std::fabs(1.0 - (4.0*
m1_2)/
m3_2))*lambda113*lambda113/(32.0*
m3*M_PI);
7350double Gammaphi3_HpHm=
HSTheta(
m3 - 2.0*sqrt(
mHp2))*sqrt(std::fabs(1.0 - (4.0*
mHp2)/
m3_2))*lambdaphi3HpHm*lambdaphi3HpHm/(32.0*
m3*M_PI);
7376 + Gammaphi3_tautau + Gammaphi3_mumu + Gammaphi3_WW
7377 + Gammaphi3_ZZ + Gamma_phi3gaga + Gamma_phi3Zga + Gamma_phi3gg
7378 + Gammaphi3_phi1phi1 + Gammaphi3_phi2phi2 + Gammaphi3_phi1phi2
7379 + Gammaphi3_HpHm + Gammaphi3_phi1Z + Gammaphi3_phi2Z
7427 throw std::runtime_error(
"Direct Searches are only available in the A2HDM.");
7431 gslpp::complex i = gslpp::complex::i();
7436 double Mcp =
Mc*(1.+4.*alsmuc/3./M_PI+alsmuc*alsmuc/M_PI/M_PI*(-1.0414*(1.-4.*
Ms/3.*
Mc)+13.4434));
7438 double Mtp =
myGTHDM->getMtpole();
7442 double Mc2 = Mcp*Mcp;
7443 double Mb2 = Mbp*Mbp;
7444 double Mt2 = Mtp*Mtp;
7445 double r2hot =
mHp*
mHp/Mtp/Mtp;
7446 double r2bot = Mbp*Mbp/Mtp/Mtp;
7447 double Vtb =
myGTHDM->getCKM().getV_tb().abs();
7448 double Vcb =
myGTHDM->getCKM().getV_cb().abs();
7449 double Vcs =
myGTHDM->getCKM().getV_cs().abs();
7466 double GammaHpcs = 0.;
7467 double GammaHpcb = 0.;
7468 double GammaHptb = 0.;
7469 double GammaHpHlW = 0.;
7470 double GammaHpphi2W = 0.;
7471 double GammaHpphi3W = 0.;
7472 double GammaHpphi3Wst = 0.;
7477 -
sd.abs2()*Mb2*(Mb2-
mHp2+Mt2)-
su.abs2()*Mt2*(Mb2-
mHp2+Mt2))
7484 -
sd.abs2()*Mb2*(Mb2-
mHp2+Mc2)-
su.abs2()*Mc2*(Mb2-
mHp2+Mc2))
7491 -
sd.abs2()*Ms2*(Ms2-
mHp2+Mc2)-
su.abs2()*Mc2*(Ms2-
mHp2+Mc2))
7524 if(
mHp >= 40.0 && mHp <= 93.0 && mH3 >= 12.0 &&
mH3 <= (
mHp - 3.))
7530 GammaHpcs + GammaHpcb + GammaHptb +
7531 GammaHpHlW + GammaHpphi2W + GammaHpphi3W + GammaHpphi3Wst;
7550 double GammatHpb = 0.;
7551 double pre_tHpb =
GF * Mt2 * Mtp * Vtb * Vtb / 4. / sqrt(2.) / M_PI;
7552 gslpp::complex aCD =
su -
sd * Mbp / Mtp;
7553 gslpp::complex bCD =
su +
sd * Mbp / Mtp;
7554 double aCD2 = (aCD * aCD.conjugate()).real();
7555 double bCD2 = (bCD * bCD.conjugate()).real();
7559 GammatHpb = pre_tHpb * ((1. - r2hot + r2bot) * (aCD2 + bCD2) / 2. + r2bot * (aCD2 - bCD2)) *
7560 sqrt(1. + r2hot*r2hot + r2bot*r2bot - 2. * (r2hot + r2bot + r2hot*r2bot)) / 2.;
7563 double GammatWb = 1.35;
7573 double MZ2 =
MZ *
MZ;
7574 double MW2 =
MW *
MW;
7575 double preZ =
GF *
MZ * MZ2 / 24. / sqrt(2.) / M_PI;
7576 double preW =
GF *
MW * MW2 / 24. / sqrt(2.) / M_PI;
7578 double gZAH_2 =
R22 *
R22;
7580 double gWHpH_2 =
R22 *
R22;
7581 double gWHpA_2 =
R33 *
R33;
7584 double func_ZHpHm = 1. - 4. *
mHp2 / MZ2;
7588 double Gamma_ZAH =
HSTheta(
MZ-
mH3-
mH2) * preZ * gZAH_2 * sqrt(std::fabs(func_ZAH)) * func_ZAH;
7589 double Gamma_ZHpHm =
HSTheta(
MZ-2.*
mHp) * preZ * gZHpHm_2 * sqrt(std::fabs(func_ZHpHm)) * func_ZHpHm;
7590 double Gamma_WHpH =
HSTheta(
MW-
mHp-
mH2) * preW * gWHpH_2 * sqrt(std::fabs(func_WHpH)) * func_WHpH;
7591 double Gamma_WHpA =
HSTheta(
MW-
mHp-
mH3) * preW * gWHpA_2 * sqrt(std::fabs(func_WHpA)) * func_WHpA;
7606 throw std::runtime_error(
"Direct Searches are only aviable in the A2HDM.");
7610 gslpp::complex i = gslpp::complex::i();
7612 double Br_Ztoee=0.03363;
7613 double Br_Ztomumu=0.03366;
7614 double Br_Ztotautau=0.0337;
7615 double Br_Ztoinv=0.2;
7617 double Br_Wtoenu=0.1071;
7618 double Br_Wtomunu=0.1063;
7619 double Br_Wtotaunu=0.1138;
7620 double Br_Wtohadrons=0.6741;
7621 double Br_tautoleptons=0.3521;
7622 double Br_tautohadrons=1-Br_tautoleptons;
7998 +
GTHDM_BR_h_ZZ*2.0*Br_Ztoinv*(Br_Ztoee+Br_Ztomumu+Br_Ztotautau*Br_tautoleptons*Br_tautoleptons));
8000 +
GTHDM_BR_h_ZZ*2.0*Br_Ztoinv*(Br_Ztoee+Br_Ztomumu+Br_Ztotautau*Br_tautoleptons*Br_tautoleptons));
8007 +2.0*
GTHDM_BR_h_WW*
GTHDM_BR_h_tautau*Br_tautohadrons*Br_tautohadrons*(Br_Wtoenu+Br_Wtomunu+Br_Wtotaunu*Br_tautohadrons)*(Br_Wtoenu+Br_Wtomunu+Br_Wtotaunu*Br_tautohadrons)
8010 +2.0*
BrSM_htoWW*
BrSM_htotautau*Br_tautohadrons*Br_tautohadrons*(Br_Wtoenu+Br_Wtomunu+Br_Wtotaunu*Br_tautohadrons)*(Br_Wtoenu+Br_Wtomunu+Br_Wtotaunu*Br_tautohadrons)
8014 +2.0*
GTHDM_BR_h_WW*
GTHDM_BR_h_tautau*Br_tautohadrons*Br_tautohadrons*(Br_Wtoenu+Br_Wtomunu+Br_Wtotaunu*Br_tautohadrons)*(Br_Wtoenu+Br_Wtomunu+Br_Wtotaunu*Br_tautohadrons)
8017 +2.0*
BrSM_htoWW*
BrSM_htotautau*Br_tautohadrons*Br_tautohadrons*(Br_Wtoenu+Br_Wtomunu+Br_Wtotaunu*Br_tautohadrons)*(Br_Wtoenu+Br_Wtomunu+Br_Wtotaunu*Br_tautohadrons)
8058 (
GTHDM_BR_h_WW*((Br_Wtoenu+Br_Wtomunu)*(Br_Wtoenu+Br_Wtomunu)
8059 +2.0*(Br_Wtoenu+Br_Wtomunu)*Br_Wtohadrons)
8062 (
BrSM_htoWW*((Br_Wtoenu+Br_Wtomunu)*(Br_Wtoenu+Br_Wtomunu)
8063 +2.0*(Br_Wtoenu+Br_Wtomunu)*Br_Wtohadrons)
8068 (
GTHDM_BR_h_WW*((Br_Wtoenu+Br_Wtomunu)*(Br_Wtoenu+Br_Wtomunu)
8069 +2.0*(Br_Wtoenu+Br_Wtomunu)*Br_Wtohadrons)
8072 (
BrSM_htoWW*((Br_Wtoenu+Br_Wtomunu)*(Br_Wtoenu+Br_Wtomunu)
8073 +2.0*(Br_Wtoenu+Br_Wtomunu)*Br_Wtohadrons)
8208 if(
m2>= 450.0 &&
m2<1400.0) {
8215 if(
m3>= 450.0 &&
m3<1400.0){
8225 if(
m2>= 400.0 &&
m2<1000.0) {
8232 if(
m3>= 400.0 &&
m3<1000.0){
8240 if(
m2>= 400.0 &&
m2<1000.0) {
8248 if(
m3>= 400.0 &&
m3<1000.0){
8257 if(
m2>= 100.0 &&
m2< 900.0) {
8265 if(
m3>= 100.0 &&
m3< 900.0) {
8272 if(
m2>= 330.0 &&
m2<1200.0) {
8279 if(
m3>= 330.0 &&
m3<1200.0) {
8293 if(
m2>= 350.0 &&
m2<650.0) {
8300 if(
m3>= 350.0 &&
m3<650.0) {
8315 if(
m2>= 550.0 &&
m2<1200.0) {
8322 if(
m3>= 550.0 &&
m3<1200.0) {
8336 if(
m2>= 50.0 &&
m2<350.0) {
8343 if(
m3>= 50.0 &&
m3<350.0) {
8359 if(
m2>= 300.0 &&
m2<1300.0) {
8366 if(
m3>= 300.0 &&
m3<1300.0) {
8383 if(
m2>= 120.0 &&
m2<500.0) {
8390 if(
m3>= 120.0 &&
m3<500.0) {
8399 if(
m2>= 120.0 &&
m2<500.0) {
8406 if(
m3>= 120.0 &&
m3<500.0) {
8426 if(
m2>= 140.0 &&
m2<1000.0) {
8439 if(
m3>= 140.0 &&
m3<1000.0) {
8448 if(
m2>= 140.0 &&
m2<1000.0) {
8462 if(
m3>= 140.0 &&
m3<1000.0) {
8482 if(
m2>= 200.0 &&
m2<1000.0) {
8492 if(
m3>= 200.0 &&
m3<1000.0) {
8501 if(
m2>= 200.0 &&
m2<1000.0) {
8508 if(
m3>= 200.0 &&
m3<1000.0) {
8522 if(
m2>= 90.0 &&
m2<1000.0) {
8530 if(
m3>= 90.0 &&
m3<1000.0) {
8541 if(
m2>= 90.0 &&
m2<1000.0) {
8548 if(
m3>= 90.0 &&
m3<1000.0) {
8555 if(
m2>= 90.0 &&
m2<1000.0) {
8562 if(
m3>= 90.0 &&
m3<1000.0) {
8569 if(
m2>= 90.0 &&
m2<1000.0) {
8576 if(
m3>= 90.0 &&
m3<1000.0) {
8583 if(
m2>= 200.0 &&
m2<2250.0) {
8594 if(
m3>= 200.0 &&
m3<2250.0) {
8604 if(
m2>= 90.0 &&
m2<3200.0) {
8614 if(
m3>= 90.0 &&
m3<3200.0) {
8621 if(
m2>= 200.0 &&
m2<2250.0) {
8630 if(
m3>= 200.0 &&
m3<2250.0) {
8637 if(
m2>= 60.0 &&
m2<3500.0) {
8644 if(
m3>= 60.0 &&
m3<3500.0) {
8653 if(
m2>= 65.0 &&
m2< 600.0) {
8660 if(
m3>= 65.0 &&
m3< 600.0) {
8667 if(
m2>= 160.0 &&
m2<3000.0) {
8679 if(
m3>= 160.0 &&
m3<3000.0)
8685 if(
m2>= 600.0 &&
m2<5000.0)
8691 if(
m3>= 600.0 &&
m3<5000.0)
8697 if(
m2>= 200.0 &&
m2<1600.0)
8704 if(
m3>= 200.0 &&
m3<1600.0)
8712 if(
m2>= 200.0 &&
m2<1200.0)
8720 if(
m3>= 200.0 &&
m3<1200.0)
8728 if(
m2>= 250.0 &&
m2<2400.0)
8736 if(
m3>= 250.0 &&
m3<2400.0)
8744 if(
m2>=1000.0 &&
m2<6800.0)
8752 if(
m3>=1000.0 &&
m3<6800.0)
8761 if(
m2>= 350.0 &&
m2<4000.0)
8769 if(
m3>= 350.0 &&
m3<4000.0)
8777 if(
m2>= 140.0 &&
m2<1000.0)
8785 if(
m3>= 140.0 &&
m3<1000.0)
8793 if(
m2>= 140.0 &&
m2<1000.0)
8802 if(
m3>= 140.0 &&
m3<1000.0)
8810 if(
m2>= 210.0 &&
m2<1200.0)
8818 if(
m3>= 210.0 &&
m3<1200.0)
8826 if(
m2>= 210.0 &&
m2<2000.0)
8834 if(
m3>= 210.0 &&
m3<2000.0)
8842 if(
m2>= 300.0 &&
m2<3000.0)
8850 if(
m3>= 300.0 &&
m3<3000.0)
8858 if(
m2>= 300.0 &&
m2<3000.0)
8866 if(
m3>= 300.0 &&
m3<3000.0)
8875 if(
m2>= 130.0 &&
m2<3000.0)
8883 if(
m3>= 130.0 &&
m3<3000.0)
8891 if(
m2>=1000.0 &&
m2<4000.0)
8899 if(
m3>=1000.0 &&
m3<4000.0)
8907 if(
m2>= 300.0 &&
m2<1500.0)
8915 if(
m3>= 300.0 &&
m3<1500.0)
8923 if(
m2>= 300.0 &&
m2<1500.0)
8931 if(
m3>= 300.0 &&
m3<1500.0)
8942 if(
m2>= 200.0 &&
m2<3000.0)
8950 if(
m3>= 200.0 &&
m3<3000.0)
8958 if(
m2>= 200.0 &&
m2<3000.0)
8966 if(
m3>= 200.0 &&
m3<3000.0)
8984 if(
m2>= 1000.0 &&
m2<4500.0)
8992 if(
m3>= 1000.0 &&
m3<4500.0)
9000 if(
m2>= 1000.0 &&
m2<4500.0)
9008 if(
m3>= 1000.0 &&
m3<4500.0)
9028 if(
m2>= 250.0 &&
m2<4000.0)
9036 if(
m3>= 250.0 &&
m3<4000.0)
9044 if(
m2>= 250.0 &&
m2<3000.0)
9052 if(
m3>= 250.0 &&
m3<3000.0)
9060 if(
m2>= 200.0 &&
m2<1000.0)
9068 if(
m3>= 200.0 &&
m3<1000.0)
9076 if(
m2>= 300.0 &&
m2<3000.0)
9084 if(
m3>= 300.0 &&
m3<3000.0)
9092 if(
m2>= 300.0 &&
m2<3000.0)
9100 if(
m3>= 300.0 &&
m3<3000.0)
9108 if(
m2>=1000.0 &&
m2<4400.0)
9116 if(
m3>=1000.0 &&
m3<4400.0)
9124 if(
m2>= 145.0 &&
m2<1000.0)
9130 if(
m3>= 145.0 &&
m3<1000.0)
9136 if(
m2>=1200.0 &&
m2<3000.0)
9143 if(
m3>=1200.0 &&
m3<3000.0)
9156 if(
m2>=300.0 &&
m2<5000.0)
9164 if(
m3>=300.0 &&
m3<5000.0)
9174 if(
m2>=300.0 &&
m2<5000.0)
9182 if(
m3>=300.0 &&
m3<5000.0)
9195 if(
m2>= 260.0 &&
m2<1000.0)
9202 if(
m3>= 260.0 &&
m3<1000.0)
9210 if(
m2>= 270.0 &&
m2<1100.0)
9218 if(
m3>= 270.0 &&
m3<1100.0)
9226 if(
m2>= 260.0 &&
m2<1100.0)
9234 if(
m3>= 260.0 &&
m3<1100.0)
9242 if(
m2>= 260.0 &&
m2< 350.0)
9250 if(
m3>= 260.0 &&
m3< 350.0)
9258 if(
m2>= 350.0 &&
m2<1000.0)
9266 if(
m3>= 350.0 &&
m3<1000.0)
9274 if(
m2>= 250.0 &&
m2<3000.0)
9282 if(
m3>= 250.0 &&
m3<3000.0)
9290 if(
m2>= 260.0 &&
m2<1200.0)
9298 if(
m3>= 260.0 &&
m3<1200.0)
9306 if(
m2>=1000.0 &&
m2<3000.0)
9314 if(
m3>=1000.0 &&
m3<3000.0)
9322 if(
m2>= 251.0 &&
m2<1000.0)
9330 if(
m3>= 251.0 &&
m3<1000.0)
9337 if(
m2>= 260.0 &&
m2< 1000.0)
9345 if(
m3>= 260.0 &&
m3< 1000.0)
9363 if(
m2>= 251.0 &&
m2<1600.0)
9371 if(
m3>= 251.0 &&
m3<1600.0)
9379 if(
m2>= 1000.0 &&
m2<3000.0)
9387 if(
m3>= 1000.0 &&
m3<3000.0)
9396 if(
m2>= 250.0 &&
m2< 900.0)
9404 if(
m3>= 250.0 &&
m3< 900.0)
9412 if(
m2>= 900.0 &&
m2<4000.0)
9420 if(
m3>= 900.0 &&
m3<4000.0)
9428 if(
m2>= 260.0 &&
m2< 900.0)
9436 if(
m3>= 260.0 &&
m3< 900.0)
9446 if(
m2>= 250.0 &&
m2< 1000.0)
9454 if(
m3>= 250.0 &&
m3< 1000.0)
9465 if(
m2>= 800.0 &&
m2< 3500.0)
9473 if(
m3>= 800.0 &&
m3< 3500.0)
9487 if(
m2>= 260.0 &&
m2< 1000.0)
9495 if(
m3>= 260.0 &&
m3< 1000.0)
9503 if(
m2>= 250.0 &&
m2< 1000.0)
9511 if(
m3>= 250.0 &&
m3< 1000.0)
9523 if(
m2>= 800.0 &&
m2< 4500.0)
9529 if(
m3>= 800.0 &&
m3< 4500.0)
9538 if(
m2>= 250.0 &&
m2 < 899.0)
9544 if(
m3>= 250.0 &&
m3< 899.0)
9552 if(
m2>= 500.0 &&
m2< 3000.0)
9560 if(
m3>= 500.0 &&
m3< 3000.0)
9568 if(
m2>= 260.0 &&
m2< 500.0)
9576 if(
m3>= 260.0 &&
m3< 500.0)
9584 if(
m2>= 220.0 &&
m2<1000.0)
9592 if(
m3>= 220.0 &&
m3<1000.0)
9600 if(
m2>= 225.0 &&
m2< 600.0)
9608 if(
m3>= 225.0 &&
m3< 600.0)
9616 if(
m2>= 220.0 &&
m2<1000.0)
9624 if(
m3>= 220.0 &&
m3<1000.0)
9632 if(
m2>= 220.0 &&
m2< 350.0)
9640 if(
m3>= 220.0 &&
m3< 350.0)
9646 if(
m2>= 220.0 &&
m2<2000.0)
9656 if(
m3>= 220.0 &&
m3<2000.0)
9664 if(
m2>= 220.0 &&
m2< 800.0)
9671 if(
m3>= 220.0 &&
m3< 800.0)
9677 if(
m2>= 800.0 &&
m2<2000.0)
9685 if(
m3>= 800.0 &&
m3<2000.0)
9694 if(
m2>= 220.0 &&
m2<2000.0)
9702 if(
m3>= 220.0 &&
m3<2000.0)
9711 if(
m2>= 220.0 &&
m2<400.0)
9719 if(
m3>= 220.0 &&
m3<400.0)
9733 if(
m2>= 220.0 &&
m2< 800.0)
9741 if(
m3>= 220.0 &&
m3< 800.0)
9749 if(
m2>= 800.0 &&
m2<2000.0)
9757 if(
m3>= 800.0 &&
m3<2000.0)
9765 if(
m3>= 175.0 && m3<1000.0 && m2 >=50.0 &&
m2 <910.0)
9773 if(
m2>= 175.0 && m2<1000.0 && m3 >=50.0 &&
m3 <910.0)
9781 if(
m3>= 50.0 && m3<1000.0 && m2 >=50.0 &&
m2 <1000.0)
9788 if(
m2>= 50.0 && m2<1000.0 && m3 >=50.0 &&
m3 <1000.0)
9795 if(
m3>= 50.0 && m3<1000.0 && m2 >=50.0 &&
m2 <1000.0)
9802 if(
m2>= 50.0 && m2<1000.0 && m3 >=50.0 &&
m3 <1000.0)
9809 if(
m2 >= 230.0 && m2 <800.0 && m3>=130.0 && m3<700.0 && m2-m3>=100)
9818 if(
m3 >= 230.0 && m3 <800.0 && m2>=130.0 && m2<700.0 && m3-m2>=100)
9825 if(
m2 >= 230.0 && m2 <800.0 && m3>=130.0 && m3<700.0 && m2-m3>=100)
9831 if(
m3 >= 230.0 && m3 <800.0 && m2>=130.0 && m2<700.0 && m3-m2>=100)
9840 if(
m2 >= 300.0 && m2 <800.0 && m3>=200.0 && m3<700.0 && m2-m3>=100)
9846 if(
m3 >= 300.0 && m3 <800.0 && m2>=200.0 && m2<700.0 && m3-m2>=100)
9855 if(
mHp>= 180.0 &&
mHp<1000.0)
9875 if(
mHp>= 180.0 &&
mHp< 600.0)
9890 if(
mHp>= 180.0 &&
mHp<2000.0)
9905 if(
mHp>= 180.0 &&
mHp<3000.0)
9919 if(
mHp>= 200.0 &&
mHp< 600.0)
9933 if(
mHp>= 180.0 &&
mHp< 600.0)
9946 if(
mHp>= 200.0 &&
mHp<2000.0)
9963 if(
mHp>= 200.0 &&
mHp<3000.0)
9998 if(
mH3 >= 15.0 &&
mH3 <= 62.0)
10012 if(
mH3 >= 12.0 &&
mH3 <= 60.0)
10026 if(
mH3 >= 15.0 &&
mH3 <= 62.0)
10033 if(
mH3 >= 4.2045 &&
mH3 <= 34.843)
10040 if(
mH3 >= 4.2 &&
mH3 <= 59.946)
10047 if(
mH3 >= 15.0 &&
mH3 <= 62.0)
10055 if(
mH3 >= 4.0 &&
mH3 <= 15.0)
10062 if(
mH3 >= 15.0 &&
mH3 <= 350.0)
10069 if(
mH3 >= 15.0 &&
mH3 <= 350.0)
10076 if(
mH3 >= 15.0 &&
mH3 <= 60.0)
10083 if(
mH3 >= 16.0 &&
mH3 <= 62.0)
10090 if((
mH3 >= 1.5 &&
mH3 <= 2.0)||(
mH3 >= 4.5 &&
mH3 <= 8.5)||(
mH3 >= 11.5 &&
mH3 <= 14.999)||(
mH3 >= 15.0 &&
mH3 <= 60.0))
10098 if(
mH3 >= 15.0 &&
mH3 <= 30.0)
10106 if(
mH3 >= 20.0 &&
mH3 <= 60.0)
10113 if(
mH3 >= 15.0 &&
mH3 <= 30.0)
10120 if(
mH3 >= 20.0 &&
mH3 <= 60.0)
10127 if((
mH3 >= 20.0 &&
mH3 <= 24.999)||(
mH3 >= 25.0 &&
mH3 <= 34.999)||(
mH3 >= 35.0 &&
mH3 <= 44.999)||(
mH3 >= 45.0 &&
mH3 <= 54.999)||(
mH3 >= 55.0 &&
mH3 <= 60.0))
10134 if(
mH3 >= 20.0 &&
mH3 <= 90.0)
10141 if(
mH3 >= 10.0 &&
mH3 <= 62.0)
10148 if(
mH3 >= 3.7 &&
mH3 <= 50.0)
10155 if(
mH3 >= 5.0 &&
mH3 <= 15.0)
10162 if(
mH3 >= 25.0 &&
mH3 <= 62.5)
10169 if(
mH3 >= 15.0 &&
mH3 <= 62.5)
10176 if(
mH3 >= 25.0 &&
mH3 <= 70.0)
10183 if(
mH3 >= 25.0 &&
mH3 <= 80.0)
10190 if(
mH3 >= 15.1 &&
mH3 <= 73.3)
10197 if(
mH3 >= 25.0 &&
mH3 <= 60.0)
10204 if(
mH3 >= 10.0 &&
mH3 <= 62.0)
10211 if(
mH3 >= 12.0 &&
mH3 <= 60.0)
10222 if(
mH2 >= 4.2045 &&
mH2 <= 34.843)
10229 if(
mH2 >= 4.2 &&
mH2 <= 59.946)
10236 if((
mH2 >= 1.5 &&
mH2 <= 2.0)||(
mH2 >= 4.5 &&
mH2 <= 8.5)||(
mH2 >= 11.5 &&
mH2 <= 14.999)||(
mH2 >= 15.0 &&
mH2 <= 60.0))
10244 if(
mH2 >= 15.0 &&
mH2 <= 30.0)
10252 if(
mH2 >= 20.0 &&
mH2 <= 60.0)
10259 if(
mH2 >= 15.0 &&
mH2 <= 30.0)
10266 if(
mH2 >= 20.0 &&
mH2 <= 60.0)
10273 if((
mH2 >= 20.0 &&
mH2 <= 24.999)||(
mH2 >= 25.0 &&
mH2 <= 34.999)||(
mH2 >= 35.0 &&
mH2 <= 44.999)||(
mH2 >= 45.0 &&
mH2 <= 54.999)||(
mH2 >= 55.0 &&
mH2 <= 60.0))
10280 if(
mH2 >= 80.0 &&
mH2 <= 109.9)
10287 if(
mH2 >= 70.0 &&
mH2 <= 110.0)
10294 if(
mH2 >= 66.0 &&
mH2 <= 110.0)
10301 if(
mH2 >= 20.0 &&
mH2 <= 116.0)
10308 if(
mH2 >= 12.0 &&
mH2 <= 120.0)
10315 if(
mH2 >= 4.0 &&
mH2 <= 120.0)
10366 if(
mHp >= 90.0 &&
mHp <= 160.0)
10373 if(
mHp >= 80.0 &&
mHp <= 160.0)
10380 if(
mHp >= 90.0 &&
mHp <= 150.0)
10387 if(
mH3 >= 15.0 && mH3 <= 75.0 && mHp >= (
mH3+85.0) &&
mHp <= 160.)
10394 if(
mHp >= 80.0 &&
mHp <= 160.0)
10401 if(
mHp >= 80.0 &&
mHp <= 160.0)
10408 if(
mHp >= 90.0 &&
mHp <= 150.0)
10415 if(
mH3 >= 15.1 && mH3 <= 70.2 && mHp >= 120.0 &&
mHp <= 160.0)
10422 if(
mHp >= 60.0 &&
mHp <= 168.0)
10428 double GeV2pb = 0.389379e9;
10432 if(
mHp >= 43.0 &&
mHp <= 95.0)
10439 if(
mHp >= 43.0 &&
mHp <= 95.0)
10446 if(
mHp >= 50.0 &&
mHp <= 93.0)
10453 if(
mHp >= 40.0 &&
mHp < 50.0)
10461 if(
mHp >= 40.0 && mHp <= 93.0 && mH3 >= 12.0 &&
mH3 <= (
mHp - 3.))
10470 if(
mHp >= 80.0 &&
mHp < 440.0)
10478 if(
mHp >= 90.0 &&
mHp < 500.0)
10486 if(
mHp >= 90.0 &&
mHp < 350.0)
10494 if(
mHp >= 100.0 &&
mHp < 800.0)
10502 if(
mHp >= 45.0 &&
mHp <= 102.0)
10513 if( RGEorder ==
"LO" ) flag=0;
10517 throw std::runtime_error(
"RGEorder can be only \"LO\" at the moment");
10522 double g1_at_MZ=sqrt(4.0*M_PI*
Ale/
cW2);
10523 double g2_at_MZ=sqrt(4.0*M_PI*
Ale/(1-
cW2));
10524 double g3_at_MZ=sqrt(4.0*M_PI*
Als);
10525 double v1_at_MZ=0.;
10526 double v2_at_MZ=0.;
10527 double etaU1_at_MZ=0.;
10528 double etaU2_at_MZ=0.;
10529 double etaD1_at_MZ=0.;
10530 double etaD2_at_MZ=0.;
10531 double etaL1_at_MZ=0.;
10532 double etaL2_at_MZ=0.;
10533 double m11sq_at_MZ=0.;
10534 double m12sq_at_MZ=0.;
10535 double m22sq_at_MZ=0.;
10536 double lambda1_at_MZ=
lambda1;
10537 double lambda2_at_MZ=
lambda2;
10538 double lambda3_at_MZ=
lambda3;
10539 double lambda4_at_MZ=
lambda4;
10543 double RpepsGTHDM=
myGTHDM->getRpepsGTHDM();
10544 double NLOuniscale=
myGTHDM->getNLOuniscaleGTHDM();
10548 double Q_A2HDM = log10(
MZ);
10550 if(fabs(Q_A2HDM-log10(
MZ))<0.005)
10578 double InitVals[21];
10579 InitVals[0]=g1_at_MZ;
10580 InitVals[1]=g2_at_MZ;
10581 InitVals[2]=g3_at_MZ;
10582 InitVals[3]=v1_at_MZ;
10583 InitVals[4]=v2_at_MZ;
10584 InitVals[5]=etaU1_at_MZ;
10585 InitVals[6]=etaU2_at_MZ;
10586 InitVals[7]=etaD1_at_MZ;
10587 InitVals[8]=etaD2_at_MZ;
10588 InitVals[9]=etaL1_at_MZ;
10589 InitVals[10]=etaL2_at_MZ;
10590 InitVals[11]=m11sq_at_MZ;
10591 InitVals[12]=m12sq_at_MZ;
10592 InitVals[13]=m22sq_at_MZ;
10593 InitVals[14]=lambda1_at_MZ;
10594 InitVals[15]=lambda2_at_MZ;
10595 InitVals[16]=lambda3_at_MZ;
10596 InitVals[17]=lambda4_at_MZ;
10597 InitVals[18]=lambda5_at_MZ;
10598 InitVals[19]=lambda6_at_MZ;
10599 InitVals[20]=lambda7_at_MZ;
10627 throw std::runtime_error(
"RGE are only defined for the CP conserving ATHDM at the moment");
10672 return std::numeric_limits<double>::quiet_NaN();
10891 if(
myGTHDM->getATHDMflag() ==
true)
10983 if (
myGTHDM->getModelName() ==
"GeneralTHDMZ2") {
const GeneralTHDM & myGTHDM
virtual double RGEGeneralTHDMRunner(double InitialValues[], unsigned long int NumberOfRGEs, double Q1, double Q2, int order, double Rpeps, double tNLOuni)
double ip_low_pp_h_phi3phi3_mumutautau_CMS8(double mass)
double THoEX_gg_phi3_phi1Z_bbll_CMS8
double ip_ex_pp_phi_phi1phi1_bbWW_ATLAS13(double mass)
double THoEX_gg_h_phi3Z_mumull_ATLAS13
double gg_phi2_ZZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_csr_ggA_tc_13(double mass)
Interpolating function for the gluon-gluon fusion A cross section ratio of the charm and top-loop and...
gslpp::matrix< double > CMS13_t_Hpb_WAb_Wmumub
CMS observed upper branching fraction limits at 13 TeV, depending on charged and pseudoscalar masses...
double ip_low_pp_h_phi3phi3_bbmumu_CMS13_cache[2][CacheSize]
double ip_Br_HPtoWW_cache[2][CacheSize]
double ip_low_HpHm_qqqq_LEP209_cache[2][CacheSize]
gslpp::complex B00_MZ2_MW2_mHp2_mHp2(const double MZ2, const double MW2, const double mHp2) const
double THoEX_pp_phi3_phi1phi1_bbtautau_1_ATLAS13
double THoEX_gg_phi2_VV_llqq_ATLAS13
double pp_phi2_phi1phi1_bbZZ_TH13
double ip_ex_bb_phi_bb_CMS13(double mass)
double ip_low_pp_bbphi3_bbmumu_CMS8(double mass)
gslpp::complex A_A_U(const double mA2, const double cW2, const double Mc, const double Mt, const double MZ) const
double THoEX_VV_phi2_VV_llqq_ATLAS13
gslpp::matrix< double > CMS13_pp_h_phi3phi3_tautautautau
double gg_phi2_phi3Z_WWZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex A_A_D_cache[6][CacheSize]
double ip_ex_pp_phi_ZZ_llqqnunull_CMS13_cache[2][CacheSize]
gslpp::matrix< double > CMS13_gg_phi_WW
gslpp::complex I_A_D(const double mA2, const double Ms, const double Mb) const
double beta_h_c
beta function evaluated at Mc and m1_2
double THoEX_bb_phi2_tautau_CMS8
void runGeneralTHDMparameters(std::string RGEorder)
double ip_ex_pp_Hpm_taunu_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a singly charged scalar resonance decaying...
double ip_integral_x_1mx_1px_G(double wa, double wb)
gslpp::complex B0p_MZ2_mA2_0_mHp2(const double MZ2, const double mA2, const double mHp2) const
gslpp::complex B00_MZ2_0_MW2_mHl2(const double MZ2, const double MW2, const double mHl2) const
double ip_ex_bb_phi_tautau_CMS13_cache[2][CacheSize]
double THoEX_bb_phi2_phi3Z_bbZ_ATLAS13
gslpp::complex B0_MZ2_mHh2_mHp2_mHp2(const double MZ2, const double mHh2, const double mHp2) const
double ip_cs_WtoWH_8(double mass)
Interpolating function for the W associated H production cross section at 8 TeV.
gslpp::complex A_A_L(const double mA2, const double cW2, const double Mmu, const double Mtau, const double MZ) const
gslpp::complex B0_MZ2_0_mHp2_mHl2(const double MZ2, const double mHp2, const double mHl2) const
double VV_phi3_VV_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_low_pp_bbphi3_bbtautau_CMS8(double mass)
double ip_ex_gg_phi_tautau_CMS13(double mass)
double ip_low_pp_h_phi3phi3_bbmumu_CMS8_cache[2][CacheSize]
gslpp::matrix< double > CMS13_pp_h_phi23phi23_mumumumu
double THoEX_pp_phi2_phi1phi1_bbVV_CMS13
gslpp::complex B0_MZ2_mHh2_mHh2_mHl2(const double MZ2, const double mHh2, const double mHl2) const
gslpp::matrix< double > ATLAS13_bb_phii_phijZ_bbZ
double THoEX_gg_phi2_WW_CMS13
gslpp::matrix< double > CMS13_gg_phi_phi1Z_tautaull
double ip_low_pp_h_phi3phi3_gagagaga_ATLAS8_cache[2][CacheSize]
double ip_ex_pp_Hp_taunu_CMS8_cache[2][CacheSize]
double ip_csr_ggA_tc_8(double mass)
Interpolating function for the gluon-gluon fusion A cross section ratio of the charm and top-loop and...
gslpp::complex B0p_MZ2_mHp2_0_mA2(const double MZ2, const double mHp2, const double mA2) const
double THoEX_gg_phi3_VV_llqq_ATLAS13
gslpp::complex B0_MZ2_mHl2_mHp2_mHp2(const double MZ2, const double mHl2, const double mHp2) const
gslpp::complex B0_MZ2_MW2_MW2_mHh2_cache[4][CacheSize]
gslpp::matrix< double > CMS13_pp_ttphi3_tttautau
gslpp::matrix< double > CMS13_pp_HpHm_taunutaunu
double interpolate2Dv2(gslpp::matrix< double > arrayTab, double x, double y)
Linearly interpolates a table with two parameter dimensions. In this case the y variable changes firs...
gslpp::complex B0_MZ2_mHl2_mA2_mA2_cache[4][CacheSize]
gslpp::matrix< double > ATLAS13_pp_phi_phi1phi1_bbtautau_2
gslpp::complex B0_MZ2_0_0_mHl2_cache[3][CacheSize]
double ip_low_phi2Z_gagaZ_LEP209_cache[2][CacheSize]
double ip_ex_gg_phi_tautau_ATLAS13_cache[2][CacheSize]
void computephi2quantities()
gslpp::matrix< double > readTable(std::string filename, int rowN, int colN)
This function reads values from a table and returns them as an array.
double THoEX_pp_phi3_phi1phi1_bbbb_ATLAS13
double gg_phi3_bb_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::complex sigmau_ATHDM
gslpp::matrix< double > ATLAS13_pp_h_phi23phi23_gagagg
double pp_phi2_phi1phi1_bbbb_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double THoEX_pp_phi2_phi1phi1_bbWW_qqlnu_CMS13
gslpp::complex B00_MZ2_MZ2_MZ2_mHl2(const double MZ2, const double mHl2) const
double pph13
Ratio of GTHDM and SM cross sections for the production of h at 13 TeV.
gslpp::complex B0_MZ2_mHh2_0_mHp2(const double MZ2, const double mHh2, const double mHp2) const
double ip_Br_HPtocc(double mass)
Interpolating function for the SM branching ratio to two charm quarks.
gslpp::matrix< double > ATLAS8_gg_phi_phi1phi1
double ip_Br_HPtoWW(double mass)
Interpolating function for the SM branching ratio to two bosons.
double ip_ex_pp_phii_phijZ_bbll_1_CMS8_cache[3][CacheSize]
double THoEX_gg_phi3_bb_CMS8
gslpp::complex B0_MZ2_mHl2_mHh2_mHh2(const double MZ2, const double mHl2, const double mHh2) const
GeneralTHDMRunner * myRunnerGTHDM
double THoEX_pp_h_phi3phi3_tautautautau_CMS13
double ip_ex_VV_phi_WW_CMS13(double mass)
double ip_ex_pp_phi_phi1phi1_bbZZ_llnunu_CMS13_cache[2][CacheSize]
double THoEX_pp_phi2_phi1phi1_bbbb_1_CMS13
double THoEX_gg_phi2_Zga_llga_ATLAS13
gslpp::complex A_h_U(const double mHl2, const double cW2, const double Mu, const double Mc, const double Mt, const double MZ) const
gslpp::complex I_HH_D(const double mHh2, const double Ms, const double Mb) const
double THoEX_t_Hpb_WAb_Wmumub_ATLAS13
double THoEX_gg_h_phi3phi3_mumumumu_ATLAS13
gslpp::complex B0_MZ2_mHh2_0_0_cache[3][CacheSize]
gslpp::complex yd1
Coupling of the SM-Higgs to up quarks.
gslpp::matrix< double > ATLAS13_pp_phi_phi1phi1_bbWW
gslpp::complex B0p_MZ2_mHl2_0_0(const double MZ2, const double mHl2) const
gslpp::complex B0p_MZ2_mHl2_mHh2_mHh2(const double MZ2, const double mHl2, const double mHh2) const
double ip_ex_tt_phi_tt_ATLAS13(double mass)
double ip_ex_gg_phi_phi1Z_bbZ_1_CMS13(double mass)
double pp_phi3_phi1phi1_bbbb_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double THoEX_VV_phi3_WW_enumunu_ATLAS13
gslpp::complex B0p_MZ2_mHh2_mHl2_mHl2(const double MZ2, const double mHh2, const double mHl2) const
double ip_integral_x2_G(double wa, double wb)
double ip_csr_ggH_tc_13(double mass)
Interpolating function for the gluon-gluon fusion H cross section ratio of the top-loop and the total...
double GTHDM_BR_h_WW
branching ratio to two bosons in the GTHDM.
double ip_ex_gg_phi_phi1Z_bbll_CMS8(double mass)
double gg_phi2_VV_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_pp_phi_phi1phi1_bbgaga_CMS8_cache[2][CacheSize]
gslpp::matrix< double > CMS13_pp_h_phi23Z_mumull
gslpp::complex B00_MZ2_MW2_mA2_mHp2(const double MZ2, const double MW2, const double mA2, const double mHp2) const
double THoEX_gg_phi3_phi1phi1_ATLAS8
double THoEX_gg_phi3_tautau_ATLAS13_low
double ip_low_pp_h_phi3phi3_tautautautau_CMS8_cache[2][CacheSize]
gslpp::complex A_h_D(const double mHl2, const double cW2, const double Md, const double Ms, const double Mb, const double MZ) const
gslpp::complex B0_MZ2_mHh2_mA2_mA2_cache[4][CacheSize]
double ip_cs_pptobbH_8(double mass)
Interpolating function for the bottom associated H production cross section at 8 TeV.
double ip_low_pp_phi2_gaga_ATLAS13_cache[2][CacheSize]
double THoEX_pp_phi2_phi3Z_tautaull_2_CMS8
void computeWZquantities()
gslpp::complex B0p_MZ2_0_mHp2_mHh2(const double MZ2, const double mHp2, const double mHh2) const
double gg_phi3_phi1Z_tautaull_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double THoEX_gg_phi3_phi1Z_tautauZ_ATLAS8
double THoEX_pp_Hpm_tb_ATLAS13
double ip_ex_VV_phi_VV_llqq_ATLAS13(double mass)
gslpp::complex B0_MZ2_mHp2_mHp2_mHl2(const double MZ2, const double mHp2, const double mHl2) const
double ip_ex_gg_phi_WW_CMS13_cache[2][CacheSize]
double ip_integral_x2_1px_G_cache[3][CacheSize]
double pp_phi3_gaga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex B0_MZ2_0_mHp2_mHh2_cache[4][CacheSize]
double ip_csr_ggA_b_13(double mass)
Interpolating function for the gluon-gluon fusion A cross section ratio of the bottom-loop and the to...
double THoEX_gg_phi3_phi2Z_WWZ_ATLAS13
double ip_susy_HpHm_munumunu_LEP208(double mass)
double pp_phi2_phi1phi1_bbtautau_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double pp_phi2_VV_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > integral_x2_1mx_G_variable_set_1_log
double ip_Br_HPtotautau(double mass)
Interpolating function for the SM branching ratio to two tau leptons.
double THoEX_gg_phi3_Zga_llga_ATLAS13
gslpp::matrix< double > ATLAS13_pp_ttphi3_ttmumu
double THoEX_VV_phi3_VV_llqq_ATLAS13
double ip_ex_gg_phi_phi1Z_tautaull_CMS13(double mass)
double THoEX_t_Hpb_WAb_Wmumub_CMS13
double ip_ex_bb_phi_tt_ATLAS13_cache[2][CacheSize]
double THoEX_pp_HpHm_munumunu_ATLAS13
double interpolate2DtriangularData(gslpp::matrix< double > arrayTab, double x, double y)
Linearly interpolates a table with two parameter dimensions. In this case the y variable changes firs...
double ip_low_pp_ttphi3_ttmumu_CMS13(double mass)
double ip_low_pp_h_phi3phi3_bbmumu_ATLAS13_cache[2][CacheSize]
double THoEX_Vh_h_phi3phi3_bbbb_ATLAS13
gslpp::matrix< double > CMS13_pp_phi2_gaga
double THoEX_phi2Z_tautauZ_LEP209
double THoEX_VV_phi3_WW_ATLAS8
double pp_phi3_phi1phi1_bbgaga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex I_h_D_cache[5][CacheSize]
double ip_ex_pp_phi_phi1phi1_bbtautau_2_ATLAS13(double mass)
double ip_ex_pp_phi_phi1phi1_bbgaga_CMS13(double mass)
double THoEX_phi2phi3_tautautautau_LEP209
gslpp::complex B0_MZ2_mA2_0_mHl2(const double MZ2, const double mA2, const double mHl2) const
double ggVV_phi3_WW_lnulnu_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex B0p_MZ2_mHl2_mHp2_mHp2_cache[4][CacheSize]
double ip_ex_gg_phi_phi1Z_bbZ_ATLAS13(double mass)
double THoEX_gg_phi3_WW_CMS13
gslpp::complex B0_MZ2_mHp2_mHp2_mHh2(const double MZ2, const double mHp2, const double mHh2) const
double ip_cs_ggtoHp_8_cache[3][CacheSize]
gslpp::complex I_A_D_cache[4][CacheSize]
gslpp::complex I_A_L(const double mA2, const double Mmu, const double Mtau) const
gslpp::complex A_HH_D(const double mHh2, const double cW2, const double Ms, const double Mb, const double MZ) const
double pp_phi3_phi1phi1_bbbb_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_bb_phi_phi1Z_bbZ_1_CMS13_cache[2][CacheSize]
double ip_integral_Hp_A_Wstar(double mi, double mj, double di, double dj)
double pp_phi3_phi1phi1_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double THoEX_pp_phi3_ZZ_qqnunu_CMS13
gslpp::matrix< double > integral_x_1mx2_G_log
gslpp::complex B00_MZ2_0_mHl2_mA2_cache[4][CacheSize]
gslpp::matrix< double > ATLAS13_gg_phii_phijZ_bbZ
double ip_integral_x2_1mx_G_variable_set_1_cache[2][CacheSize]
double SigmaZh13
SM cross sections for the associated production of h with a Z boson at 13 TeV at LHC.
gslpp::matrix< gslpp::complex > Nl_GTHDM
double THoEX_gg_phi2_mumu_CMS8
gslpp::complex A_A_D(const double mA2, const double cW2, const double Ms, const double Mb, const double MZ) const
double ip_integral_x2_1mx_G_variable_set_1(double wb)
double ip_ex_bsgamma(double logtb, double logmHp)
Interpolating function for the NNLO value for the branching ratio of decays in the GTHDM.
double SigmaggF13
SM cross sections for the production of h at 13 TeV at LHC through gluon fusion.
double ip_ex_pp_phi_phi1phi1_bbbb_2_CMS13(double mass)
gslpp::complex B0p_MZ2_mHh2_mA2_mA2(const double MZ2, const double mHh2, const double mA2) const
double THoEX_pp_phi2_phi1phi1_bbtautau_1_CMS13
double ip_integral_x_1mx_1px_G_cache[3][CacheSize]
double ip_ex_gg_phi_phi1Z_bbZ_1_CMS13_cache[2][CacheSize]
double ip_ex_bb_phi_mumu_CMS13(double mass)
double ip_ex_pp_phi_phi1phi1_4WOr2W2tauOr4tau_CMS13_cache[2][CacheSize]
gslpp::complex B0_MZ2_0_MW2_mHl2(const double MZ2, const double MW2, const double mHl2) const
double pp_phi3_phi1phi1_bbtautau_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_csr_ggA_tc_13_cache[2][CacheSize]
double ip_Br_HPtoZZ(double mass)
Interpolating function for the SM branching ratio to two bosons.
double THoEX_bb_phi3_phi1Z_bbZ_ATLAS13
gslpp::matrix< double > ATLAS13_pp_h_phi3phi3_gagagaga
void computephi3quantities()
double THoEX_gg_phi3_WW_lnuqq_ATLAS13
double ip_cs_pptottA_8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_bb_phi_tt
double ip_low_pp_h_phi23phi23_gagagg_ATLAS13_cache[2][CacheSize]
double ip_low_t_Hpb_csb_CMS8_cache[2][CacheSize]
double ip_Br_HPtoZZ_cache[2][CacheSize]
double THoEX_t_Hpb_taunub_ATLAS8
double ip_ex_gg_phii_phijZ_bbZ_ATLAS13(double m3, double m2)
double ip_ex_pp_phi_phi1phi1_bbgaga_ATLAS13(double mass)
double THoEX_Vh_h_phi2phi2_bbbb_ATLAS13
gslpp::matrix< double > CMS13_pp_phi_phi1phi1_bbtautau_2
gslpp::complex B0p_MZ2_mA2_0_mHl2_cache[4][CacheSize]
gslpp::complex B0_MZ2_mHl2_mHh2_mHl2_cache[4][CacheSize]
gslpp::complex B0p_MZ2_mHh2_mHp2_mHp2_cache[4][CacheSize]
double GTHDM_BR_h_AZ
branching ratio to bosons in the GTHDM.
gslpp::matrix< double > CMS13_pp_phi2_bb_light
gslpp::complex B0p_MZ2_mHp2_0_mHl2_cache[4][CacheSize]
gslpp::matrix< double > CMS13_pp_h_phi3phi3_bbtautau
double THoEX_pp_phi3_gaga_ATLAS13
double THoEX_pp_h_phi3phi3_mumutautau_CMS8
gslpp::complex B0p_MZ2_mHl2_mHl2_mHl2_cache[3][CacheSize]
gslpp::matrix< double > ATLAS13_pp_phi_VV_qqqq
double beta_h_b
beta function evaluated at Mb and m1_2
gslpp::matrix< double > CMS8_pp_Hp_taunu
CMS observed upper cross section limits at 8 TeV, depending on the charged Higgs mass.
gslpp::complex A_HH_D_cache[6][CacheSize]
gslpp::complex B0_MZ2_mA2_mA2_mHl2_cache[4][CacheSize]
double THoEX_gg_phi2_phi1phi1_ATLAS8
double THoEX_gg_phi2_gaga_ATLAS8
double VV_phi2_WW_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_cs_pptobbA_13(double mass)
Interpolating function for the bottom associated A production cross section at 13 TeV.
double ip_cs_ZtoZH_13(double mass)
Interpolating function for the Z associated H production cross section at 13 TeV.
double THoEX_pp_phi2_gaga_CMS8
gslpp::matrix< double > CMS8_pp_h_phi3phi3_tautautautau
double THoEX_pp_phi2_phi1phi1_bbtautau_CMS8
gslpp::complex A_HH_L_cache[6][CacheSize]
gslpp::complex B0_MZ2_0_mA2_mHl2_cache[4][CacheSize]
double ip_low_pp_h_phi3phi3_bbmumu_ATLAS13(double mass)
double ip_ex_pp_phii_phijZ_bbll_2_CMS8_cache[3][CacheSize]
gslpp::matrix< double > ATLAS8_gg_phi_tautau
ATLAS observed upper cross section limits at 8 TeV, depending on the Higgs mass.
double ip_ex_pp_phi_phi1phi1_bbZZ_lljj_CMS13(double mass)
gslpp::matrix< double > CMS13_pp_phi_phi1phi1_4WOr2W2tauOr4tau
double ip_low_HpHm_AWstAWst_bbbb_OPAL209_cache[3][CacheSize]
void CacheShiftReal(double cache[][CacheSize], const int NumPar, const double params[], const double newResult) const
Adds a new result and its parameters into the cache.
double ip_ex_gg_phi_tautau_ATLAS8_cache[2][CacheSize]
double ip_low_HpHm_qqqq_LEP209(double mass)
double ip_ex_pp_Hpm_tb_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > br_tt
SM Higgs branching ratio tables (obtained with HDECAY 6.10), depending on the Higgs mass.
double beta_h_mu
beta function evaluated at Mmu and m1_2
double KaellenFunction(const double a2, const double b2, const double c2) const
Kaellen function.
double THoEX_VV_phi3_ZZ_qqllnunu_ATLAS13
gslpp::complex B0_MZ2_0_MW2_mHh2(const double MZ2, const double MW2, const double mHh2) const
gslpp::matrix< double > LEP209_phi2Z_gagaZ
LEP observed upper branching fraction limits for GeV, depending on the scalar mass.
double yu3R
Coupling of A to leptons.
gslpp::complex B0_MZ2_MZ2_MZ2_mHl2(const double MZ2, const double mHl2) const
gslpp::complex B0p_MZ2_mHh2_mHp2_mHp2(const double MZ2, const double mHh2, const double mHp2) const
gslpp::matrix< double > ATLAS13_gg_phi_phi1Z_bbZ
gslpp::complex B0p_MZ2_mHp2_0_mHh2_cache[4][CacheSize]
double bb_phi3_bb_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex B0p_MZ2_mHh2_0_mA2(const double MZ2, const double mHh2, const double mA2) const
double ip_ex_pp_phi_gaga_ATLAS13(double mass)
double THoEX_pp_phi2_ZZ_qqnunu_CMS13
gslpp::complex B0p_MZ2_mHh2_mHh2_mHh2(const double MZ2, const double mHh2) const
double rh_gaga
Squared relative coupling of to two photons.
double ip_ex_VV_phi_ZZ_llllnunu_ATLAS13(double mass)
double ip_low_gg_h_phi3phi3_tautautautau_ATLAS8(double mass)
double ip_cs_ggtoH_13_cache[2][CacheSize]
double THoEX_phi1phi3_tautautautau_LEP209
double ip_ex_gg_phi_Zga_llga_ATLAS13(double mass)
gslpp::matrix< double > ATLAS8_VV_phi_ZZ
double ip_ex_gg_phi_WW_ATLAS8(double mass)
double ip_ex_gg_phi_gaga_CMS13(double mass)
gslpp::complex I_A_U_cache[4][CacheSize]
gslpp::complex I_HH_D_cache[4][CacheSize]
gslpp::complex B00_MZ2_MZ2_mHp2_mHp2_cache[3][CacheSize]
double ip_ex_bb_phii_phijZ_bbZ_ATLAS13(double m3, double m2)
gslpp::matrix< double > ATLAS13_pp_phi_phi1phi1_bbtautau_1
double THoEX_bb_phi3_tautau_ATLAS13
gslpp::complex I_A_L_cache[4][CacheSize]
double ip_low_HpHm_taunuAWst_taunubb_OPAL209(double mi, double mj, double di, double dj)
gslpp::complex B0p_MZ2_mHl2_mHh2_mHh2_cache[4][CacheSize]
double GTHDM_BR_h_HZ
branching ratio to bosons in the GTHDM.
gslpp::complex B0_MZ2_mHl2_0_mA2(const double MZ2, const double mHl2, const double mA2) const
double BrSM_htogaga
SM branching ratio of .
double THoEX_gg_h_phi3phi3_tautautautau_ATLAS8
double THoEX_pp_h_phi3phi3_gagagaga_ATLAS13
gslpp::complex B0p_MZ2_mHl2_0_mA2_cache[4][CacheSize]
double VV_phi3_ZZ_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_low_t_Hpb_taunub_CMS8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_pp_h_phi3phi3_bbmumu
ATLAS observed upper cross section (or branching fraction) limits at 13 TeV, depending on the pseudo...
gslpp::matrix< double > ATLAS8_pp_Hpm_tb
gslpp::matrix< double > CMS13_pp_phi_phi1phi1_bbVV
gslpp::complex yd2
Coupling of H to up quarks.
gslpp::complex B0_MZ2_0_MW2_mHh2_cache[4][CacheSize]
double ip_cs_pptobbH_8_cache[2][CacheSize]
double ip_ex_gg_phi_WW_enumunu_ATLAS13_cache[2][CacheSize]
double pp_phi3_phi2Z_tautaull_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::matrix< double > ATLAS8_pp_h_phi3phi3_gagagaga
ATLAS & CMS observed upper cross section (or branching fraction) limits at 8 TeV,...
double THoEX_gg_phi2_phi1phi1_gagaWW_ATLAS13
double ip_ex_gg_phi_mumu_CMS8(double mass)
gslpp::matrix< double > ATLAS13_t_Hpb_csb
ATLAS observed upper branching fraction limits at 13 TeV, depending on the charged scalar mass.
double ip_ex_pp_phi_Zga_llga_ATLAS8_cache[2][CacheSize]
double pp_phi2_phi1phi1_with_channel_bbWW_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double THoEX_pp_phi2_WW_lnuqq_CMS13
gslpp::matrix< double > CMS8_pp_phii_phijZ_tautaull_2
double GTHDM_BR_h_gaga
branching ratio to two photons in the GTHDM.
double THoEX_ggVV_phi3_WW_lnulnu_CMS13
gslpp::matrix< double > ATLAS8_gg_phi_WW
double THoEX_gg_phi2_phi1Z_tautaull_CMS8
double ip_cs_pptobbA_8_cache[2][CacheSize]
double ip_low_gg_phi3_tautau_ATLAS13(double mass)
double THoEX_pp_phi2_gaga_CMS13
double ip_ex_pp_phi_phi1phi1_bbZZ_llnunu_CMS13(double mass)
gslpp::matrix< double > csrA_top_charm_13
double ip_ex_bb_phi_phi1Z_bbZ_1_CMS13(double mass)
double ip_csr_ggH_tc_8(double mass)
Interpolating function for the gluon-gluon fusion H cross section ratio of the top-loop and the total...
double rh_gg
Squared relative coupling of to two gluons.
gslpp::complex B0p_MZ2_0_mA2_mHl2_cache[4][CacheSize]
gslpp::complex B00_MZ2_MZ2_mHh2_mA2_cache[4][CacheSize]
double ip_ex_pp_phi_phi1phi1_bbbb_1_CMS13(double mass)
double ip_integral_Hp_A_Wstar_cache[3][CacheSize]
double ip_ex_VV_phi_WW_lnuqq_ATLAS13(double mass)
double rh_ggO
Squared relative coupling of to two gluons, odd part.
double THoEX_pp_phi2_phi3Z_tautaull_1_CMS8
double GTHDM_BR_h_cc
branching ratio to two quarks in the GTHDM.
gslpp::complex B00_MZ2_0_mA2_mHp2_cache[4][CacheSize]
double THoEX_pp_h_phi3phi3_gagagaga_ATLAS8
gslpp::matrix< double > CMS8_bb_phi_tautau
double THoEX_pp_phi3_ZZ_llqqnunull_CMS13
double Gammaphi2tot
Total decay width of the heavy CP-even Higgs .
gslpp::complex B0_MZ2_mHl2_0_0(const double MZ2, const double mHl2) const
double THoEX_pp_phi3_phi1phi1_bbtautau_2_CMS13
gslpp::complex I_HH_U(const double mHh2, const double Mc, const double Mt) const
gslpp::matrix< double > CMS8_pp_phii_phijZ_bbll_1
double ip_low_Vh_h_phi3phi3_bbbb_CMS13(double mass)
gslpp::complex B0p_MZ2_mHp2_0_mHl2(const double MZ2, const double mHp2, const double mHl2) const
double THoEX_bb_phi2_mumu_ATLAS13
double ip_csr_HpHm_pp_13_LO_cache[2][CacheSize]
gslpp::complex B0_MZ2_MW2_MW2_mHl2(const double MZ2, const double MW2, const double mHl2) const
double ip_Br_HPtobb_cache[2][CacheSize]
~GeneralTHDMcache()
GeneralTHDMcache destructor.
gslpp::matrix< double > br_ZZ
double ip_low_pp_h_phi3phi3_mumutautau_CMS8_cache[2][CacheSize]
double ip_ex_gg_phi_gaga_ATLAS8_cache[2][CacheSize]
gslpp::complex B0p_MZ2_0_mHp2_mA2_cache[4][CacheSize]
double gg_phi2_tautau_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_low_pp_h_phi23phi23_bbmumu_ATLAS13_old(double mass)
double ip_Br_HPtomumu(double mass)
Interpolating function for the SM branching ratio to two muons.
double ip_low_t_Hpb_WAb_Wmumub_CMS13(double mass)
double ip_integral_x_1mx2_G(double wa, double wb)
gslpp::complex B0_MZ2_mHh2_mHl2_mHl2_cache[4][CacheSize]
gslpp::complex B0p_MZ2_mA2_mA2_mHl2_cache[4][CacheSize]
gslpp::matrix< double > CMS13_gg_phi_phi1Z_bbZ_2
double beta_h_tau
beta function evaluated at Mtau and m1_2
double ip_low_HpHm_taunuAWst_taunubb_OPAL209_cache[3][CacheSize]
double gg_phi3_phi2Z_bbZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > br_tautau
double THoEX_gg_phi2_phi1Z_bbZ_2_CMS13
double pp_phi3_phi2Z_bbll_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_low_Vh_h_phi23phi23_bbbb_ATLAS13(double mass)
double ip_ex_pp_Hpm_tb_ATLAS8_cache[2][CacheSize]
double gg_phi2_WW_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double gg_phi3_phi1Z_tautauZ_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_gg_phi_ZZ_ATLAS8_cache[2][CacheSize]
double THoEX_gg_phi3_tautau_CMS13
double gg_phi3_gaga_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::matrix< double > ATLAS13_bb_phi_bb
double ip_low_pp_h_phi23phi23_gagagg_ATLAS13(double mass)
double gg_phi3_Zga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_cs_ZtoZH_13_cache[2][CacheSize]
double BrSM_Ztoll
branching ratio to in the SM.
double THoEX_pp_h_phi2phi2_gagagg_ATLAS13
double ip_cs_pptobbA_8(double mass)
Interpolating function for the bottom associated A production cross section at 8 TeV.
gslpp::complex A_h_L(const double mHl2, const double cW2, const double Me, const double Mmu, const double Mtau, const double MZ) const
double bb_phi3_phi2Z_bbZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double SigSM_pph13
SM cross sections for the production of h at 13 TeV at LHC.
double THoEX_bb_phi3_phi1Z_bbZ_1_CMS13
double ip_cs_ggtoA_8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_pp_phi_gaga
gslpp::matrix< double > CMS13_VV_phi_WW_heavy
double THoEX_VV_phi2_ZZ_llllnunu_ATLAS13
double THoEX_gg_phi3_phi1Z_tautaull_CMS13
int CacheCheck(const gslpp::complex cache[][CacheSize], const int NumPar, const double params[]) const
Check whether for the latest set of parameters a value is in the cache.
gslpp::complex B0p_MZ2_0_mHp2_mHl2_cache[4][CacheSize]
gslpp::complex B0_MZ2_mHp2_0_mHl2(const double MZ2, const double mHp2, const double mHl2) const
gslpp::complex B0p_MZ2_mHp2_0_mHh2(const double MZ2, const double mHp2, const double mHh2) const
double gg_phi3_mumu_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double pp_phi2_VV_TH8
Signal strength for the process with $VV=WW,ZZ$ at the LHC with 8 TeV.
double ip_ex_pp_Hpm_taunu_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a singly charged scalar resonance decayi...
gslpp::matrix< double > CMS8_pp_bbphi3_bbmumu
gslpp::matrix< double > CMS8_pp_phi_Zga_llga
double BrSM_htotautau
SM branching ratio of .
gslpp::matrix< double > CMS13_bb_phi_phi1Z_bbZ_1
gslpp::complex B0_MZ2_MZ2_MZ2_mHh2_cache[3][CacheSize]
double ip_ex_pp_phi_phi1phi1_bbbb_ATLAS13(double mass)
double ip_ex_pp_phii_phijZ_tautaull_1_CMS8(double m2, double m3)
gslpp::matrix< double > log_cs_bbH_13
SM Higgs production cross section table at 13 TeV obtained with SusHi 1.5, depending on the Higgs mas...
double ip_cs_pptottA_13_cache[2][CacheSize]
gslpp::complex B0_MZ2_mA2_mA2_mHl2(const double MZ2, const double mA2, const double mHl2) const
double THoEX_pp_phi3_VV_CMS8
double ip_ex_pp_phi_phi1phi1_bbtautau_2_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_bb_phi_phi1Z_bbZ
gslpp::complex I_h_D(const double mHl2, const double Md, const double Ms, const double Mb) const
double SigmaVBF8
SM cross sections for the production of h through vector boson fusion at 8 TeV at LHC.
gslpp::complex A_h_D_cache[7][CacheSize]
double pp_phi3_phi1phi1_bbZZ_TH13
double ip_ex_gg_phi_tautau_ATLAS13(double mass)
double THoEX_pp_Hpm_tb_ATLAS8
gslpp::complex B0_MZ2_mHh2_0_mA2(const double MZ2, const double mHh2, const double mA2) const
double ip_ex_gg_phii_phijZ_bbZ_ATLAS13_cache[3][CacheSize]
double THoEX_bb_phi2_phi1Z_bbZ_2_CMS13
double ip_ex_pp_phii_phijZ_tautaull_2_CMS8_cache[3][CacheSize]
gslpp::matrix< double > CMS13_VV_phi_WW
gslpp::matrix< double > integral_x2_G_log
gslpp::complex B00_MZ2_0_MZ2_mHh2(const double MZ2, const double mHh2) const
double ip_low_phiEphi3_tautautautau_LEP209(double mi, double mj, double di, double dj)
gslpp::matrix< double > ATLAS13_gg_phi_mumu
double THoEX_gg_phi2_phi1Z_bbll_CMS8
double ip_low_pp_h_phi3phi3_tautautautau_CMS13(double mass)
double ip_ex_pp_phi_phi1phi1_bbVV_CMS13_cache[2][CacheSize]
double ip_ex_pp_phi_phi1phi1_bbgaga_CMS13_cache[2][CacheSize]
gslpp::matrix< double > CMS13_gg_phi_mumu
double ip_low_pp_phi2_gaga_CMS8(double mass)
gslpp::matrix< double > CMS13_pp_Hpm_tb
gslpp::complex A_H_W(const double mH, const double cW2, const double MW, const double MZ) const
double ip_Br_HPtocc_cache[2][CacheSize]
double gg_phi3_phi2Z_WWZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double pp_Hpm_taunu_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex B0p_MZ2_mA2_0_mHl2(const double MZ2, const double mA2, const double mHl2) const
gslpp::complex B0p_MZ2_mHh2_mHl2_mHl2_cache[4][CacheSize]
double THoEX_pp_bbphi3_bbtautau_CMS8
double THoEX_pp_h_phi3phi3_bbmumu_CMS8
gslpp::matrix< double > CMS8_gg_phi_phi1phi1_bbtautau
gslpp::complex B0_MZ2_mHh2_0_0(const double MZ2, const double mHh2) const
double VBF_Vh
Ratio of GTHDM and SM cross sections for VBF and Vh production of h.
double THoEX_pp_phi2_phi1phi1_bbtautau_2_CMS13
double THoEX_gg_phi2_WW_lnuqq_ATLAS13
double ip_low_t_Hpb_csb_ATLAS13(double mass)
double ip_low_gg_phi3_tautau_ATLAS13_cache[2][CacheSize]
double rh_QuQuO
Squared relative coupling of to two up quarks, odd part.
double THoEX_pp_phi3_phi1phi1_bbZZ_lljj_CMS13
double ip_ex_pp_phi2_bb_light_CMS13(double mass)
double pp_phi3_WW_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_pp_Hp_tb_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a singly charged scalar resonance decaying...
double ip_low_pp_h_phi3phi3_gagagaga_ATLAS13_cache[2][CacheSize]
double GTHDM_BR_h_ZZ
branching ratio to two bosons in the GTHDM.
double THoEX_gg_phi2_ZZ_llllnunu_ATLAS13
gslpp::matrix< double > CMS13_pp_phi_phi1phi1_bbWWorbbtautau
double ip_ex_gg_phii_phijZ_WWZ_ATLAS13(double m3, double m2)
gslpp::complex I_A_U(const double mA2, const double Mc, const double Mt) const
gslpp::complex B0_MZ2_mHh2_mHh2_mHl2_cache[4][CacheSize]
gslpp::complex yd3
Coupling of A to up quarks.
gslpp::complex A_HH_U_cache[6][CacheSize]
double lambdaipm(const double R1i, const double R2i, const double R3i) const
double THoEX_gg_phi3_phi1phi1_bbtautau_CMS8
gslpp::complex B0_MZ2_mHp2_0_mHh2_cache[4][CacheSize]
double bb_phi3_mumu_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::complex B0_MZ2_mA2_mA2_mHh2(const double MZ2, const double mA2, const double mHh2) const
double ip_low_pp_phi2_gaga_CMS13(double mass)
double THoEX_pp_HpHm_munumunu_CMS13
double bb_phi2_bb_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::matrix< double > CMS13_gg_phi_Zga
double ip_ex_VV_phi_ZZ_ATLAS8(double mass)
double THoEX_gg_phi2_tautau_CMS8
double ip_integral_x2_G_variable_set_1(double wb)
double ip_low_pp_bbphi3_bbtautau_CMS13(double mass)
double ip_ex_pp_Hp_taunu_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a singly charged scalar resonance decaying...
double beta(const double mf, const double m_2) const
beta function
double ip_low_HpHm_taunutaunu_LEP209(double mass)
double pp_phi3_phi1phi1_bbVV_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_low_HpHm_qqtaunu_OPAL209(double mass)
double ip_cs_ggtoH_8(double mass)
Interpolating function for the H production cross section via gluon-gluon fusion at 8 TeV.
double ip_low_Zh_h_phi23phi23_bbbb_ATLAS13(double mass)
double THoEX_pp_h_phi3phi3_mumutautau_CMS13
double ip_ex_gg_phi_phi1phi1_ATLAS8(double mass)
double ip_ex_gg_phi_phi1phi1_bbtautau_CMS8_cache[2][CacheSize]
gslpp::matrix< double > log_cs_ttH_8
SM Higgs production cross section table at 8 TeV obtained with MadGraph 5, depending on the Higgs mas...
double pp_phi3_phi1phi1_with_channel_bbtautau_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV where the decay channel o...
double ip_ex_bb_phi_mumu_CMS8_cache[2][CacheSize]
gslpp::matrix< double > LEP209_phi2Z_tautauZ
double pp_phi2_phi1phi1_bbbb_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double THoEX_pp_Hpm_taunu_ATLAS8
gslpp::matrix< double > CMS13_pp_h_phi3phi3_mumutautau
CMS observed upper cross section (or branching fraction) limits at 13 TeV, depending on the pseudosc...
double pp_phi3_phi1phi1_with_channel_bbgaga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV where the decay channel o...
double ip_csr_ggH_tc_13_cache[2][CacheSize]
gslpp::complex B0p_MZ2_0_0_mHl2(const double MZ2, const double mHl2) const
double ip_ex_pp_Hpm_tb_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a singly charged scalar resonance decayi...
double ip_low_pp_h_phi3phi3_bbtautau_ATLAS13(double mass)
double THoEX_gg_phi3_phi1Z_tautaull_CMS8
double VV_phi2_VV_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double THoEX_gg_phi3_tautau_ATLAS8
double ip_csr_ggA_b_8(double mass)
Interpolating function for the gluon-gluon fusion A cross section ratio of the bottom-loop and the to...
double ip_ex_gg_phi_mumu_CMS13(double mass)
double THoEX_gg_phi3_WW_enumunu_ATLAS13
gslpp::matrix< double > ATLAS13_bb_phi_mumu
gslpp::complex A_H_W_cache[5][CacheSize]
gslpp::complex B00_MZ2_MW2_mHl2_mHp2_cache[5][CacheSize]
double THoEX_Zh_h_phi3phi3_bbbb_ATLAS13
double ip_ex_bb_phi_tautau_ATLAS13_cache[2][CacheSize]
double pp_phi2_phi1phi1_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > CMS8_gg_phi_phi1Z_tautaull
double THoEX_gg_phi3_tautau_CMS8
double THoEX_pp_phi2_gaga_ATLAS13_low
gslpp::matrix< double > CMS13_tt_phi2_tt
gslpp::complex B0p_MZ2_mHp2_mHp2_mHh2_cache[4][CacheSize]
gslpp::matrix< double > ATLAS13_pp_HpHm_munumunu
double cW2GTHDM(const double c02) const
double THoEX_HpHm_taunuAWst_taunubb_OPAL209
gslpp::complex A_HH_U(const double mHh2, const double cW2, const double Mc, const double Mt, const double MZ) const
gslpp::matrix< double > ATLAS13_gg_phi_WW_enumunu
gslpp::complex B0p_MZ2_mHl2_0_mA2(const double MZ2, const double mHl2, const double mA2) const
double ip_ex_gg_phi_phi1Z_bbZ_2_CMS13_cache[2][CacheSize]
double ip_low_t_Hpb_cbb_ATLAS13(double mass)
double SigSM_pph8
SM cross sections for the production of h at 8 TeV at LHC.
gslpp::complex B0p_MZ2_mHh2_mHh2_mHl2_cache[4][CacheSize]
gslpp::matrix< double > CMS13_gg_phi_WW_heavy
double ip_ex_pp_phi_Zga_llga_CMS8(double mass)
gslpp::complex B0_MZ2_0_mA2_mHh2(const double MZ2, const double mA2, const double mHh2) const
double ip_ex_VV_phi_WW_ATLAS8_cache[2][CacheSize]
double ip_ex_bb_phi_mumu_ATLAS13(double mass)
gslpp::matrix< double > br_WW
double gg_phi3_tautau_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex B0p_MZ2_mHh2_0_0(const double MZ2, const double mHh2) const
double THoEX_pp_ttphi3_ttmumu_ATLAS13
double pp_phi2_phi1phi1_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::matrix< double > CMS13_pp_phi_phi1phi1_bbZZ_lljj
double bb_phi3_tautau_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_gg_phi_phi1Z_bbZ_2_CMS13(double mass)
double THoEX_pp_h_phi3phi3_bbtautau_CMS13
gslpp::matrix< double > log_cs_WH_13
double THoEX_bb_phi3_phi2Z_bbZ_ATLAS13
double ip_ex_bb_phi_tautau_CMS8_cache[2][CacheSize]
double ip_cs_ggtoHp_8(double mHp, double logtb)
Interpolating function for the H+ production cross section from two gluons at 8 TeV.
double ip_cs_WtoWH_13_cache[2][CacheSize]
double gg_phi3_ZZ_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_VV_phi_WW_CMS13_cache[2][CacheSize]
double ip_ex_pp_phi_phi1phi1_bbWWorbbtautau_CMS13(double mass)
gslpp::matrix< double > log_cs_ttA_8
double ip_ex_VV_phi_WW_enumunu_ATLAS13(double mass)
double THoEX_Vh_h_phi3phi3_bbbb_CMS13
double ip_ex_bb_phi_bb_CMS13_cache[2][CacheSize]
double THoEX_VV_phi3_WW_CMS13
gslpp::matrix< double > ATLAS13_Vh_h_phi23phi23_bbbb
gslpp::matrix< double > CMS8_gg_phi_tautau
gslpp::matrix< gslpp::complex > Nu_GTHDM
gslpp::matrix< double > CMS13_Vh_h_phi3phi3_bbbb
double gg_phi3_VV_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_cs_pptobbH_13(double mass)
Interpolating function for the bottom associated H production cross section at 13 TeV.
gslpp::complex B0_MZ2_mHl2_mHp2_mHp2_cache[4][CacheSize]
double bb_phi3_phi1Z_bbZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_pp_phi_phi1phi1_bbtautau_CMS8_cache[2][CacheSize]
double ip_low_pp_h_phi3phi3_gagagaga_ATLAS8(double mass)
double BrSM_htogg
SM branching ratio of .
gslpp::matrix< double > ATLAS13_gg_h_phi23Z_mumull
gslpp::matrix< double > OPAL209_HpHm_qqtaunu
OPAL observed upper branching fraction limits for GeV, depending on the charged masses.
double pp_Hp_taunu_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double VV_phi3_ZZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_bb_phi_mumu_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > CMS8_bb_phi_bb
CMS observed upper cross section limits at 8 TeV, depending on the Higgs mass.
double THoEX_gg_phi3_WW_ATLAS8
gslpp::complex B0_MZ2_mA2_0_mHh2_cache[4][CacheSize]
double THoEX_bb_phi3_phi1Z_bbZ_2_CMS13
double ip_ex_gg_phi_Zga_llga_ATLAS13_cache[2][CacheSize]
double ip_cs_ggtoHp_13_cache[3][CacheSize]
gslpp::matrix< double > log_cs_VBF_8
double ip_ex_gg_phi_bb_CMS8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_pp_phi_phi1phi1_bbgaga
gslpp::complex B0_MZ2_mHl2_mHl2_mHl2_cache[3][CacheSize]
double ip_ex_pp_phi_phi1phi1_bbWW_ATLAS13_cache[2][CacheSize]
gslpp::complex B0p_MZ2_mA2_mA2_mHh2_cache[4][CacheSize]
double ip_ex_gg_phi_phi1phi1_gagaWW_ATLAS13_cache[2][CacheSize]
double THoEX_t_Hpb_csb_ATLAS13
double THoEX_pp_phi3_phi1phi1_bbgaga_CMS13
double THoEX_pp_phi3_phi1phi1_bbbb_2_CMS13
double ip_susy_pp_HpHm_taunutaunu_CMS13(double mass)
double THoEX_HpHm_munumunu_LEP208
double ip_ex_pp_phi_phi1phi1_bbbb_ATLAS13_cache[2][CacheSize]
double bb_phi2_tautau_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_bb_phi_bb_ATLAS13_cache[2][CacheSize]
gslpp::complex B0_MZ2_MZ2_MZ2_mHl2_cache[3][CacheSize]
double pp_phi2_phi1phi1_with_channel_4WOr2W2tauOr4tau_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV where the SM-like Higgs d...
gslpp::matrix< double > LEP209_phiEphi3_tautautautau
double pp_Hpm_tb_TH13
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_gg_phi_VV_llqq_ATLAS13(double mass)
double THoEX_tt_phi2_tt_ATLAS13
gslpp::matrix< double > integral_x2_G_variable_set_1_log
double ip_ex_gg_phi_phi1Z_bbZ_ATLAS8_cache[2][CacheSize]
int CacheCheckReal(const double cache[][CacheSize], const int NumPar, const double params[]) const
Check whether for the latest set of parameters a value is in the cache.
double THoEX_pp_phi3_phi1phi1_bbWW_qqlnu_CMS13
gslpp::complex A_h_L_cache[7][CacheSize]
double ip_low_pp_bbphi3_bbtautau_CMS8_cache[2][CacheSize]
gslpp::complex B0_MZ2_0_mA2_mHh2_cache[4][CacheSize]
double ip_low_HpHm_taunutaunu_LEP209_cache[2][CacheSize]
gslpp::complex I_H_Hp(const double mHp2, const double mH) const
gslpp::complex B00_MZ2_0_mHh2_mA2(const double MZ2, const double mHh2, const double mA2) const
double THoEX_pp_phi2_VV_qqqq_ATLAS13
double ip_ex_VV_phi_ZZ_qqllnunu_ATLAS13_cache[2][CacheSize]
double THoEX_pp_phi3_phi1phi1_4WOr2W2tauOr4tau_CMS13
double THoEX_bb_phi3_tautau_ATLAS8
double ip_integral_perturb_ATHDM(double ratio)
double THoEX_gg_phi3_WW_heavy_CMS13
gslpp::complex B0p_MZ2_0_mHp2_mHh2_cache[4][CacheSize]
gslpp::complex B00_MZ2_MZ2_MZ2_mHh2_cache[3][CacheSize]
gslpp::matrix< double > CMS8_pp_Hp_tb
double THoEX_t_Hpb_cbb_CMS8
gslpp::matrix< double > log_cs_WH_8
double gg_phi3_WW_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > CMS8_t_Hpb_csb
CMS observed upper branching fraction limits at 8 TeV, depending on the charged scalar mass.
gslpp::complex B0p_MZ2_mHl2_mHp2_mHp2(const double MZ2, const double mHl2, const double mHp2) const
gslpp::matrix< double > ATLAS8_gg_h_phi3phi3_tautautautau
double ip_csr_ggH_b_13(double mass)
Interpolating function for the gluon-gluon fusion H cross section ratio of the bottom-loop and the to...
double ip_low_HpHm_qqtaunu_OPAL172(double mass)
gslpp::matrix< double > CMS13_pp_h_phi3phi3_gagagaga
gslpp::matrix< double > CMS13_ggVV_phi_WW_lnulnu
gslpp::matrix< double > CMS13_pp_phi_phi1phi1_bbWW_qqlnu
double GTHDM_BR_h_tautau
branching ratio to two leptons in the GTHDM.
double THoEX_gg_phi3_ZZ_llllnunu_ATLAS13
double gg_phi2_bb_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double THoEX_bb_phi3_mumu_CMS13
double ip_ex_gg_phi_Zga_CMS13(double mass)
gslpp::complex B0_MZ2_mHl2_mHh2_mHh2_cache[4][CacheSize]
double THoEX_tt_phi2_tt_CMS13
gslpp::complex B0p_MZ2_0_mA2_mHh2(const double MZ2, const double mA2, const double mHh2) const
double tt_phi2_tt_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double KaellenFunction_cache[4][CacheSize]
gslpp::matrix< double > CMS13_pp_phi_ZZ_qqnunu
gslpp::matrix< double > CMS13_gg_phi_tautau
double THoEX_VV_phi2_ZZ_qqllnunu_ATLAS13
double THoEX_gg_phi3_phi2Z_bbZ_ATLAS13
double pp_phi2_phi1phi1_with_channel_bbWWorbbtautau_TH13
double ip_ex_pp_phi_VV_CMS8_cache[2][CacheSize]
double ip_ex_tt_phi3_tt_CMS13(double mass)
double THoEX_bb_phi2_mumu_CMS8
gslpp::matrix< double > OPAL172_HpHm_qqtaunu
OPAL observed upper branching fraction limits for GeV, depending on the charged masses.
gslpp::complex B0p_MZ2_mHh2_mHh2_mHh2_cache[3][CacheSize]
double THoEX_HpHm_qqtaunu_OPAL209
gslpp::matrix< double > CMS13_pp_bbphi3_bbtautau
gslpp::matrix< double > ATLAS13_pp_phi2_gaga_low
double THoEX_pp_phi2_phi1phi1_bbbb_2_CMS13
double GTHDM_BR_h_bb
branching ratio to two quarks in the GTHDM.
double ip_ex_gg_phi_VV_llqq_ATLAS13_cache[2][CacheSize]
double ip_ex_tt_phi_tt_ATLAS13_cache[2][CacheSize]
double rh_QdQdO
Squared relative coupling of to two down quarks, odd part.
gslpp::matrix< double > CMS8_gg_phi_phi1Z_bbll
double ip_integral_x_1mx2_G_variable_set_0_cache[2][CacheSize]
int HSTheta(const double x) const
Heaviside function.
double ip_low_HpHm_qqtaunu_OPAL209_cache[2][CacheSize]
double ip_low_pp_ttphi3_ttmumu_ATLAS13(double mass)
double pp_phi2_phi3Z_bbll_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_pp_phi_bb_CMS13(double mass)
double THoEX_pp_phi3_Zga_llga_ATLAS8
double ip_low_gg_h_phi23phi23_mumumumu_ATLAS13_cache[2][CacheSize]
double bb_phi2_mumu_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::matrix< double > log_cs_ggH_13
SM Higgs production cross section tables at 13 TeV from the LHC Higgs Cross Section Working Group,...
double THoEX_bb_phi3_mumu_CMS8
double THoEX_pp_phi2_phi3Z_bbll_1_CMS8
gslpp::complex B00_MZ2_0_mHh2_mHp2_cache[4][CacheSize]
double sumModBRs
Sum of the modified branching ratios.
double ip_cs_pptottA_8(double mass)
Interpolating function for the top associated A production cross section at 8 TeV.
double ip_ex_pp_Hpm_taunu_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a singly charged scalar resonance decayi...
double THoEX_pp_phi3_phi1phi1_bbbb_CMS8
double THoEX_pp_HpHm_taunutaunu_ATLAS13
double ip_cs_ZtoZH_8(double mass)
Interpolating function for the Z associated H production cross section at 8 TeV.
gslpp::matrix< double > ATLAS13_VV_phi_ZZ_llllnunu
double bb_phi3_bb_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_VV_phi_VV_llqq_ATLAS13_cache[2][CacheSize]
double ip_low_pp_h_phi3phi3_tautautautau_CMS13_cache[2][CacheSize]
gslpp::matrix< double > CMS13_pp_phi_WW_lnuqq
gslpp::matrix< double > ATLAS8_pp_phi_Zga_llga
gslpp::matrix< double > ATLAS13_pp_HpHm_taunutaunu
SUSY observed upper branching fraction limits at 13 TeV, depending on charged scalar mass.
double ip_low_t_Hpb_cbb_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_pp_h_phi3phi3_bbtautau
gslpp::matrix< double > CMS13_bb_phi_mumu
double THoEX_pp_h_phi3phi3_bbtautau_ATLAS13
double ip_cs_VBFtoH_13(double mass)
Interpolating function for the H production cross section via vector boson fusion at 13 TeV.
double THoEX_bb_phi2_phi1Z_bbZ_ATLAS13
gslpp::complex B00_MZ2_MW2_mA2_mHp2_cache[5][CacheSize]
gslpp::complex B0_MZ2_mHh2_mHh2_mHh2_cache[3][CacheSize]
gslpp::matrix< double > ATLAS8_t_Hpb_taunub
ATLAS observed upper branching fraction limits at 8 TeV, depending on the charged scalar mass.
gslpp::matrix< double > LEP209_HpHm_taunutaunu
LEP observed upper cross-section times branching fraction limits for GeV, depending on the charged ...
gslpp::complex B0_MZ2_mHl2_0_0_cache[3][CacheSize]
double ip_susy_pp_HpHm_munumunu_CMS13(double mass)
gslpp::complex B0p_MZ2_mHh2_mHh2_mHl2(const double MZ2, const double mHh2, const double mHl2) const
gslpp::matrix< double > ATLAS13_gg_phi_WW_lnuqq
double THoEX_pp_Hp_taunu_CMS8
double THoEX_pp_ttphi3_ttmumu_CMS13
double gg_phi2_ZZ_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_cs_ggtoA_13_cache[2][CacheSize]
double ip_csr_ggA_b_13_cache[2][CacheSize]
double THoEX_gg_phi3_phi1Z_bbZ_1_CMS13
double THoEX_gg_phi2_phi3Z_WWZ_ATLAS13
double ip_low_pp_h_phi23Z_mumull_CMS13_cache[2][CacheSize]
gslpp::matrix< double > OPAL209_HpHm_AWstAWst_bbbb
OPAL observed upper branching fraction limits for GeV, depending on the charged and pseudoscalar ma...
double gg_phi2_WW_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double THoEX_bb_phi3_tt_ATLAS13
double THoEX_HpHm_taunutaunu_LEP209
double THoEX_pp_h_phi3phi3_gagagaga_CMS13
gslpp::complex B0p_MZ2_mHh2_0_mA2_cache[4][CacheSize]
double pp_phi2_WW_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_bsgamma_cache[3][CacheSize]
gslpp::complex B0p_MZ2_0_mA2_mHl2(const double MZ2, const double mA2, const double mHl2) const
double THoEX_pp_ttphi3_tttautau_CMS13
gslpp::complex I_h_L(const double mHl2, const double Me, const double Mmu, const double Mtau) const
double ip_ex_pp_phii_phijZ_tautaull_2_CMS8(double m2, double m3)
double ip_cs_pptottA_13(double mass)
Interpolating function for the top associated A production cross section at 13 TeV.
gslpp::matrix< double > ATLAS13_gg_phi_ZZ_qqllnunu
double ip_integral_x_1mx2_G_cache[3][CacheSize]
double gg_phi2_Zga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double THoEX_pp_phi2_bb_light_CMS13
gslpp::complex B0_MZ2_mHh2_0_mA2_cache[4][CacheSize]
gslpp::complex B0_MZ2_0_MZ2_mHl2(const double MZ2, const double mHl2) const
double ip_ex_VV_phi_ZZ_ATLAS8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_gg_phii_phijZ_WWZ
gslpp::complex I_HH_U_cache[4][CacheSize]
double ip_ex_gg_phi_WW_heavy_CMS13_cache[2][CacheSize]
double gg_phi3_mumu_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double gg_phi2_phi1Z_tautaull_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_pp_phii_phijZ_bbll_2_CMS8(double m2, double m3)
gslpp::matrix< double > CMS8_pp_phi_phi1phi1_bbgaga
double VV_phi2_WW_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex B0p_MZ2_0_0_mHh2_cache[3][CacheSize]
static const int CacheSize
Cache size.
double THoEX_gg_phi2_phi1Z_tautauZ_ATLAS8
double ip_ex_VV_phi_WW_heavy_CMS13(double mass)
gslpp::matrix< double > integral_x2_1mx_G_log
gslpp::matrix< double > log_cs_bbA_13
double THoEX_VV_phi2_WW_enumunu_ATLAS13
gslpp::matrix< double > CMS13_gg_phi_gaga
gslpp::matrix< double > CMS13_gg_phi_phi1Z_bbZ_1
double THoEX_HpHm_qqtaunu_OPAL172
double yl3R
Coupling of A to down quarks real part.
double THoEX_pp_phi3_phi1phi1_bbWWorbbtautau_CMS13
double ip_ex_gg_phi_ZZ_llllnunu_ATLAS13_cache[2][CacheSize]
gslpp::complex B0_MZ2_MW2_MW2_mHl2_cache[4][CacheSize]
gslpp::complex B00_MZ2_MZ2_MZ2_mHl2_cache[3][CacheSize]
double Gammaphi3tot
Total decay width of the heavy CP-even Higgs .
double ip_ex_bb_phi_phi1Z_bbZ_ATLAS13(double mass)
double THoEX_pp_phi2_bb_CMS13
double THoEX_pp_h_phi3phi3_tautautautau_CMS8
double ip_ex_pp_Hpm_taunu_CMS13_cache[2][CacheSize]
double ip_susy_pp_HpHm_munumunu_ATLAS13(double mass)
gslpp::complex B0p_MZ2_mA2_0_mHp2_cache[4][CacheSize]
double ip_ex_gg_phii_phijZ_WWZ_ATLAS13_cache[3][CacheSize]
double SigmaVh8
SM cross sections for the associated production of h with a vector boson at 8 TeV at LHC.
double ip_susy_pp_HpHm_munumunu_CMS13_cache[2][CacheSize]
double ip_ex_pp_phi_phi1phi1_bbbb_CMS8_cache[2][CacheSize]
gslpp::complex B0_MZ2_mHh2_mA2_mA2(const double MZ2, const double mHh2, const double mA2) const
gslpp::matrix< double > CMS13_pp_h_phi3phi3_bbmumu
double THoEX_pp_phi3_phi1phi1_bbWW_ATLAS13
double pp_phi3_bb_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > csrH_bottom_8
double THoEX_pp_phi2_phi1phi1_bbtautau_2_ATLAS13
gslpp::matrix< double > integral_x2_1px_G_log
double ip_ex_bb_phi_bb_CMS8(double mass)
gslpp::matrix< double > CMS13_pp_ttphi3_ttmumu
double ip_low_Zh_h_phi23phi23_bbbb_ATLAS13_cache[2][CacheSize]
gslpp::complex A_HH_L(const double mHh2, const double cW2, const double Mmu, const double Mtau, const double MZ) const
gslpp::complex B0_MZ2_mHl2_mHh2_mHl2(const double MZ2, const double mHl2, const double mHh2) const
double THoEX_gg_phi3_ZZ_qqllnunu_ATLAS13
double THoEX_gg_phi2_phi1Z_bbZ_1_CMS13
double pp_phi3_VV_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_bb_phi_bb_ATLAS13(double mass)
double ip_csr_ggH_b_8(double mass)
Interpolating function for the gluon-gluon fusion H cross section ratio of the bottom-loop and the to...
double ip_ex_bb_phi_phi1Z_bbZ_2_CMS13(double mass)
double THoEX_VV_phi2_ZZ_ATLAS8
double pp_phi2_Zga_llga_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double THoEX_gg_phi3_ZZ_ATLAS8
double bb_phi2_mumu_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_gg_phi_phi1Z_tautauZ_ATLAS8_cache[2][CacheSize]
double THoEX_pp_HpHm_taunutaunu_CMS13
double ip_low_pp_h_phi23phi23_mumumumu_CMS13(double mass)
double gg_phi3_phi1phi1_bbtautau_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double rphi3_ggE
Squared relative coupling of to two gluons, even part.
double ip_cs_ggtoA_8(double mass)
Interpolating function for the A production cross section via gluon-gluon fusion at 8 TeV.
double VV_phi2_ZZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double THoEX_t_Hpb_cbb_ATLAS13
double pp_phi2_phi1phi1_bbVV_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double THoEX_pp_h_phi3phi3_gagagg_ATLAS13
double ip_ex_gg_phi_phi1Z_tautauZ_ATLAS8(double mass)
double ip_ex_gg_phi_bb_CMS8(double mass)
gslpp::matrix< double > csrHpHm_pp_13_LO
charged higgs pair production cross-sections at 13 TeV, depending on charged scalar mass (taken from ...
double gg_phi2_tautau_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double bb_phi2_phi3Z_bbZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double pp_phi3_phi1phi1_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double THoEX_phi1phi3_bbbb_LEP209
double gg_phi2_phi1Z_bbll_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double gg_phi2_mumu_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double THoEX_bb_phi3_bb_CMS13
gslpp::matrix< double > CMS13_pp_phi_phi1phi1_bbbb_1
gslpp::matrix< gslpp::complex > Md_GTHDM
double ip_ex_gg_phi_gaga_CMS13_cache[2][CacheSize]
gslpp::complex I_HH_L(const double mHh2, const double Mmu, const double Mtau) const
double THoEX_pp_phi3_phi1phi1_bbbb_1_CMS13
double ip_ex_gg_phi_gaga_ATLAS8(double mass)
double THoEX_VV_phi3_WW_lnuqq_ATLAS13
double THoEX_VV_phi3_ZZ_ATLAS8
double ip_low_pp_phi2_gaga_ATLAS13(double mass)
double ip_ex_pp_phi_ZZ_qqnunu_CMS13_cache[2][CacheSize]
gslpp::matrix< double > CMS8_pp_h_phi3phi3_mumutautau
gslpp::matrix< double > CMS8_t_Hpb_taunub
gslpp::matrix< double > CMS8_pp_phi_phi1phi1_bbtautau
double gg_phi3_phi1phi1_gagaWW_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double THoEX_gg_phi3_gaga_CMS13
double ip_susy_pp_HpHm_munumunu_ATLAS13_cache[2][CacheSize]
double pp_phi3_phi1phi1_with_channel_bbWW_qqlnu_TH13
gslpp::matrix< gslpp::complex > Nd_GTHDM
double gg_phi2_phi1Z_tautauZ_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_gg_phi_WW_CMS13(double mass)
gslpp::complex B0p_MZ2_mHl2_mA2_mA2(const double MZ2, const double mHl2, const double mA2) const
double THoEX_bb_phi2_tautau_ATLAS13
double ip_low_t_Hpb_csb_CMS13_cache[2][CacheSize]
double ip_ex_pp_phi_phi1phi1_bbgaga_ATLAS13_cache[2][CacheSize]
gslpp::complex B00_MZ2_MW2_mHp2_mHp2_cache[4][CacheSize]
gslpp::matrix< double > CMS13_pp_phi_phi1phi1_bbtautau_1
double ip_ex_pp_phi_ZZ_llqqnunull_CMS13(double mass)
double pp_phi3_phi1phi1_with_channel_4WOr2W2tauOr4tau_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV where the SM-like Higgs d...
double ip_cs_pptottH_8(double mass)
Interpolating function for the top associated H production cross section at 8 TeV.
gslpp::matrix< double > CMS13_pp_phi_bb
CMS observed upper cross section limits at 13 TeV, depending on the Higgs mass.
double ip_low_gg_h_phi23Z_mumull_ATLAS13_cache[2][CacheSize]
double ip_susy_pp_HpHm_taunutaunu_ATLAS13(double mass)
double ip_integral_x2_G_variable_set_1_cache[2][CacheSize]
double THoEX_pp_phi2_phi1phi1_bbWWorbbtautau_CMS13
gslpp::complex I_HH_L_cache[4][CacheSize]
gslpp::complex B0p_MZ2_0_mHp2_mHl2(const double MZ2, const double mHp2, const double mHl2) const
double THoEX_gg_phi2_mumu_CMS13
double pp_Hpm_taunu_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double gg_phi2_mumu_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double gg_phi3_ZZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex B0p_MZ2_mA2_0_mHh2(const double MZ2, const double mA2, const double mHh2) const
double ip_low_gg_h_phi23phi23_mumumumu_ATLAS13(double mass)
double ip_ex_gg_phi_phi1Z_tautaull_CMS13_cache[2][CacheSize]
double ip_cs_pptobbH_13_cache[2][CacheSize]
gslpp::matrix< double > arraybsgamma
table, depending on logtb and the logarithm of the charged Higgs mass.
double ip_integral_x_1mx2_G_variable_set_0(double wb)
double ip_ex_gg_phi_phi1phi1_bbtautau_CMS8(double mass)
double THoEX_gg_phi2_tautau_CMS13
double ip_ex_pp_phii_phijZ_bbll_1_CMS8(double m2, double m3)
double THoEX_pp_bbphi3_bbtautau_CMS13
double THoEX_gg_phi2_phi1Z_bbZ_ATLAS13
double ip_ex_tt_phi3_tt_CMS13_cache[2][CacheSize]
double ip_ex_pp_phi_phi1phi1_bbZZ_lljj_CMS13_cache[2][CacheSize]
double THoEX_pp_phi2_phi1phi1_bbZZ_lljj_CMS13
gslpp::complex B0_MZ2_0_MW2_mHl2_cache[4][CacheSize]
GeneralTHDMcache(const StandardModel &SM_i)
GeneralTHDMcache constructor.
gslpp::matrix< double > ATLAS8_gg_phi_phi1Z_tautauZ
gslpp::complex I_h_U_cache[5][CacheSize]
gslpp::matrix< double > log_cs_bbA_8
double ip_low_pp_h_phi3phi3_mumutautau_CMS13_cache[2][CacheSize]
double ip_low_phiEphi3_bbbb_LEP209(double mi, double mj, double di, double dj)
double ip_ex_pp_phi_phi1phi1_bbbb_2_CMS13_cache[2][CacheSize]
gslpp::complex B0_MZ2_mHl2_0_mHp2(const double MZ2, const double mHl2, const double mHp2) const
gslpp::matrix< double > ATLAS13_pp_phi_phi1phi1_bbbb
gslpp::complex I_H_W(const double mH, const double MW) const
gslpp::complex B0_MZ2_mHl2_0_mA2_cache[4][CacheSize]
double ip_low_pp_h_phi3phi3_bbtautau_ATLAS13_cache[2][CacheSize]
double Sigma_HpHm_LEP(const double mass2, const double sqrts) const
Cross section for charged-Higgs pair production at LEP, from hep-ph/9602250.
gslpp::matrix< double > ATLAS13_gg_phi3_tautau_low
double ip_ex_pp_phi_phi1phi1_bbtautau_1_CMS13(double mass)
double ip_low_pp_phi2_gaga_CMS8_cache[2][CacheSize]
double ip_low_pp_h_phi3phi3_tautautautau_CMS8(double mass)
gslpp::matrix< double > log_cs_ttH_13
SM Higgs production cross section table at 13 TeV obtained with MadGraph 5, depending on the Higgs ma...
double THoEX_VV_phi2_WW_heavy_CMS13
gslpp::complex B00_MZ2_0_MZ2_mHh2_cache[3][CacheSize]
gslpp::complex B00_MZ2_MZ2_mHp2_mHp2(const double MZ2, const double mHp2) const
double ip_ex_gg_phi_tautau_CMS13_cache[2][CacheSize]
double ip_low_pp_h_phi3phi3_gagagaga_CMS13(double mass)
double ip_ex_gg_phi_ZZ_ATLAS8(double mass)
double gg_phi3_WW_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::matrix< double > ATLAS13_pp_Hpm_taunu
ATLAS observed upper cross section limits at 13 TeV, depending on the charged Higgs mass.
double THoEX_pp_h_phi2phi2_mumumumu_CMS13
double ip_ex_pp_phi_bb_CMS13_cache[2][CacheSize]
double ip_ex_gg_phi_mumu_ATLAS13_cache[2][CacheSize]
double gg_phi3_gaga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double BrSM_htobb
SM branching ratio of .
gslpp::matrix< double > CMS8_gg_phi_bb
double ip_low_pp_h_phi3phi3_gagagaga_CMS13_cache[2][CacheSize]
gslpp::complex A_A_U_cache[6][CacheSize]
gslpp::matrix< gslpp::complex > Ml_GTHDM
double bb_phi2_tt_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > CMS13_pp_phi_ZZ_llqqnunull
double THoEX_pp_h_phi3phi3_bbmumu_ATLAS13
gslpp::matrix< double > ATLAS8_bb_phi_tautau
double THoEX_pp_phi3_bb_CMS13
double Gamma_Z_inv
decay rate in the GTHDM.
double ip_ex_pp_phi_phi1phi1_bbbb_1_CMS13_cache[2][CacheSize]
double THoEX_phi2Z_gagaZ_LEP209
double THoEX_pp_h_phi3phi3_bbmumu_ATLAS13_old
double ip_Br_HPtobb(double mass)
Interpolating function for the SM branching ratio to two bottom quarks.
gslpp::matrix< double > log_cs_ggA_13
CP-odd Higgs production cross section tables at 13 TeV obtained with HIGLU 4.34, depending on the Hig...
double ip_ex_bb_phi_bb_CMS8_cache[2][CacheSize]
double ip_csr_HpHm_pp_13_LO(double mass)
double ip_ex_pp_phi_gaga_ATLAS13_cache[2][CacheSize]
double THoEX_bb_phi3_mumu_ATLAS13
double ip_cs_VBFtoH_8(double mass)
Interpolating function for the H production cross section via vector boson fusion at 8 TeV.
double rh_VV
Squared relative coupling of to two massive vector bosons.
gslpp::complex sigmal_ATHDM
double ip_ex_pp_phi_phi1phi1_bbbb_CMS8(double mass)
gslpp::matrix< double > CMS13_pp_Hpm_taunu
CMS observed upper cross section limits at 13 TeV, depending on the charged Higgs mass.
double ip_low_t_Hpb_taunub_ATLAS8(double mass)
double ip_integral_perturb_ATHDM_cache[2][CacheSize]
double pp_phi3_phi1phi1_bbgaga_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double gg_phi2_phi1phi1_bbtautau_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_gg_phi_mumu_CMS8_cache[2][CacheSize]
gslpp::complex B00_MZ2_MW2_MW2_mHl2_cache[4][CacheSize]
gslpp::complex B0_MZ2_0_mA2_mHl2(const double MZ2, const double mA2, const double mHl2) const
double THoEX_bb_phi2_mumu_CMS13
double ip_ex_gg_phi_phi1Z_bbZ_ATLAS13_cache[2][CacheSize]
gslpp::complex B0p_MZ2_mHl2_mHl2_mHl2(const double MZ2, const double mHl2) const
double THoEX_bb_phi2_tt_ATLAS13
double THoEX_pp_phi3_WW_lnuqq_CMS13
gslpp::complex A_h_U_cache[7][CacheSize]
double pp_phi3_phi1phi1_with_channel_bbbb_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double THoEX_bb_phi2_tautau_CMS13
double ip_low_t_Hpb_taunub_CMS8(double mass)
gslpp::complex B0_MZ2_mHl2_mA2_mA2(const double MZ2, const double mHl2, const double mA2) const
double ip_cs_ggtoA_13(double mass)
Interpolating function for the A production cross section via gluon-gluon fusion at 13 TeV.
double ip_Br_HPtotautau_cache[2][CacheSize]
double ip_low_pp_h_phi3phi3_bbmumu_CMS8(double mass)
gslpp::complex B00_MZ2_0_mHl2_mHp2_cache[4][CacheSize]
double THoEX_pp_Hpm_tb_CMS13
double THoEX_pp_phi3_bb_light_CMS13
double ip_ex_pp_phi_VV_qqqq_ATLAS13(double mass)
double THoEX_pp_phi2_phi1phi1_4WOr2W2tauOr4tau_CMS13
double GTHDM_BR_h_AA
branching ratio to two bosons in the GTHDM.
gslpp::matrix< double > CMS13_pp_phi_phi1phi1_bbWWorbbtautau_low_masses
double THoEX_pp_Hpm_taunu_CMS13
gslpp::matrix< double > ATLAS13_VV_phi_WW_lnuqq
double THoEX_pp_phi2_phi1phi1_bbtautau_1_ATLAS13
gslpp::complex B0p_MZ2_0_0_mHl2_cache[3][CacheSize]
gslpp::complex B00_MZ2_MZ2_mHl2_mA2_cache[4][CacheSize]
double THoEX_pp_h_phi3Z_mumull_CMS13
double ip_ex_pp_phi_phi1phi1_bbtautau_CMS8(double mass)
double THoEX_gg_phi3_phi1Z_bbZ_ATLAS13
double THoEX_pp_phi2_phi1phi1_bbgaga_CMS8
gslpp::matrix< double > CMS13_pp_phi_phi1phi1_bbgaga
double THoEX_pp_phi3_phi1phi1_bbtautau_2_ATLAS13
gslpp::matrix< double > ATLAS8_gg_phi_gaga
gslpp::matrix< double > csrH_bottom_13
double ip_ex_pp_phi_Zga_llga_ATLAS8(double mass)
gslpp::matrix< double > Rij_GTHDM
double pp_phi2_phi1phi1_with_channel_bbtautau_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV where the decay channel o...
double THoEX_pp_phi2_Zga_llga_ATLAS8
double rh_QdQdE
Squared relative coupling of to two down quarks, even part.
gslpp::complex A_H_Hp(const double mHp2, const double mH, const double cW2, const double MZ) const
double THoEX_pp_phi3_phi1phi1_bbgaga_ATLAS13
double THoEX_gg_h_phi2Z_mumull_ATLAS13
void computeSignalStrengths()
double pp_phi2_gaga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double rphi3_VV
Squared relative coupling of to two massive vector bosons.
double THoEX_bb_phi3_bb_ATLAS13
double THoEX_bb_phi2_phi1Z_bbZ_1_CMS13
gslpp::matrix< double > CMS8_pp_h_phi3phi3_bbmumu
double ip_low_t_Hpb_csb_ATLAS13_cache[2][CacheSize]
double THoEX_gg_phi2_phi3Z_bbZ_ATLAS13
double ip_integral_x2_1px_G(double wa, double wb)
gslpp::matrix< double > CMS13_tt_phi3_tt
double THoEX_gg_phi2_WW_heavy_CMS13
double ip_ex_pp_phi_Zga_llga_CMS8_cache[2][CacheSize]
double Gamma_h_inv
decay rate in the GTHDM.
gslpp::matrix< double > log_cs_ggHp_8
Charged Higgs production cross section table at 8 TeV from LHCHXSWGMSSMCharged, depending on the char...
gslpp::complex B00_MZ2_0_MZ2_mHl2(const double MZ2, const double mHl2) const
double ip_cs_WtoWH_13(double mass)
Interpolating function for the W associated H production cross section at 13 TeV.
double ip_ex_pp_phi_VV_qqqq_ATLAS13_cache[2][CacheSize]
double ip_low_HpHm_AWstAWst_bbbb_OPAL209(double mi, double mj, double di, double dj)
double ip_ex_pp_phi_phi1phi1_bbWWorbbtautau_CMS13_cache[2][CacheSize]
double pp_phi3_phi1phi1_with_channel_bbWW_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double THoEX_gg_phi2_Zga_CMS13
double bb_phi3_mumu_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > ATLAS13_VV_phi_ZZ_qqllnunu
gslpp::complex f_func(const double x) const
f function for the gamma gamma coupling to h, H and A
double ip_ex_gg_phi_phi1Z_bbll_CMS8_cache[2][CacheSize]
double SigmaVh13
SM cross sections for the associated production of h with a vector boson at 13 TeV at LHC.
double gg_phi2_phi1Z_bbZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > integral_x_1mx_1px_G_log
double ip_ex_bb_phi_tt_ATLAS13(double mass)
double ip_low_t_Hpb_taunub_ATLAS8_cache[2][CacheSize]
double ip_ex_pp_Hpm_tb_CMS13_cache[2][CacheSize]
gslpp::matrix< double > CMS8_pp_phi_phi1phi1_bbbb
double MWGTHDM(const double MW) const
double ip_low_pp_ttphi3_tttautau_CMS13(double mass)
double pp_Hp_tb_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double THoEX_tt_phi3_tt_ATLAS13
gslpp::matrix< double > ATLAS13_Zh_h_phi23phi23_bbbb
double THoEX_VV_phi3_WW_heavy_CMS13
double THoEX_phi2phi3_bbbb_LEP209
double ip_ex_gg_phi_phi1phi1_gagaWW_ATLAS13(double mass)
double ip_low_pp_bbphi3_bbtautau_CMS13_cache[2][CacheSize]
gslpp::matrix< double > CMS13_bb_phi_phi1Z_bbZ_2
double ip_ex_gg_phi_phi1Z_bbZ_ATLAS8(double mass)
double interpolate2D(gslpp::matrix< double > arrayTab, double x, double y)
Linearly interpolates a table with two parameter dimensions. In this case the x variable changes firs...
double ggF_tth8
Ratio of GTHDM and SM cross sections for ggF and tth production of h at 8 TeV.
double BrSM_htocc
SM branching ratio of .
double ip_ex_gg_phi_phi1Z_tautaull_CMS8_cache[2][CacheSize]
gslpp::complex B0_MZ2_0_0_mHh2(const double MZ2, const double mHh2) const
double ip_Br_HPtott(double mass)
Interpolating function for the SM branching ratio to two top quarks.
double THoEX_pp_phi2_gaga_ATLAS13
double ip_ex_pp_phii_phijZ_tautaull_1_CMS8_cache[3][CacheSize]
double ip_ex_pp_phi_VV_CMS8(double mass)
double ip_ex_pp_phi_phi1phi1_bbVV_CMS13(double mass)
double ip_low_Vh_h_phi23phi23_bbbb_ATLAS13_cache[2][CacheSize]
double ip_low_gg_h_phi23Z_mumull_ATLAS13(double mass)
double ip_GammaHPtotSM_cache[2][CacheSize]
double THoEX_pp_phi3_phi2Z_bbll_2_CMS8
double THoEX_t_Hpb_csb_CMS13
gslpp::matrix< double > csrH_top_charm_8
Production cross section ratio tables at 8 TeV obtained with HIGLU 4.34, depending on the Higgs mass.
gslpp::complex Int1(const double tau, const double lambda) const
function for Z gamma coupling to h, H and A
double THoEX_VV_phi2_WW_lnuqq_ATLAS13
double ip_low_t_Hpb_WAb_Wmumub_CMS13_cache[2][CacheSize]
double ip_ex_pp_Hpm_tb_CMS13(double mass)
gslpp::complex B0p_MZ2_mHl2_mHh2_mHl2_cache[4][CacheSize]
gslpp::complex A_H_Hp_cache[5][CacheSize]
gslpp::complex B00_MZ2_MW2_mHl2_mHp2(const double MZ2, const double MW2, const double mHl2, const double mHp2) const
double ip_integral_x2_G_cache[3][CacheSize]
gslpp::matrix< double > csrA_bottom_13
gslpp::matrix< double > LEP209_phiEphi3_bbbb
double ip_ex_gg_phi_WW_lnuqq_ATLAS13(double mass)
double VV_phi3_WW_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > csrA_bottom_8
gslpp::complex B0p_MZ2_mA2_mA2_mHl2(const double MZ2, const double mA2, const double mHl2) const
double ip_cs_ggtoH_8_cache[2][CacheSize]
double ip_ex_pp_phi_phi1phi1_bbtautau_1_CMS13_cache[2][CacheSize]
double gg_phi3_phi1Z_bbZ_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_pp_phi_phi1phi1_bbWWorbbtautau_low_masses_CMS13(double mass)
double ip_ex_pp_phi3_bb_light_CMS13(double mass)
gslpp::complex B0_MZ2_0_MZ2_mHl2_cache[3][CacheSize]
double ip_integral_x2_1mx_G(double wa, double wb)
double ip_susy_pp_HpHm_taunutaunu_CMS13_cache[2][CacheSize]
double ip_integral_x2_1mx_G_cache[3][CacheSize]
gslpp::complex B00_MZ2_0_mA2_mHp2(const double MZ2, const double mA2, const double mHp2) const
gslpp::matrix< double > CMS13_t_Hpb_csb
CMS observed upper branching fraction limits at 13 TeV, depending on the charged scalar mass.
double ip_ex_pp_phi_phi1phi1_bbtautau_1_ATLAS13(double mass)
gslpp::matrix< double > log_cs_ggHp_13
Charged Higgs production cross section table at 13 TeV from LHCHXSWGMSSMCharged, depending on the cha...
double THoEX_pp_phi2_phi1phi1_bbgaga_CMS13
gslpp::matrix< double > ATLAS13_gg_h_phi23phi23_mumumumu
gslpp::complex B00_MZ2_MZ2_mHh2_mA2(const double MZ2, const double mHh2, const double mA2) const
double gg_phi2_phi1Z_bbZ_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::complex B00_MZ2_MW2_MW2_mHh2_cache[4][CacheSize]
double ip_ex_gg_phi_phi1Z_tautaull_CMS8(double mass)
gslpp::complex B0_MZ2_0_0_mHh2_cache[3][CacheSize]
double ip_low_pp_bbphi3_bbmumu_CMS8_cache[2][CacheSize]
double yl1R
Coupling of the SM-Higgs to down quarks real part.
double THoEX_VV_phi3_ZZ_llllnunu_ATLAS13
gslpp::complex B0p_MZ2_mA2_mA2_mHh2(const double MZ2, const double mA2, const double mHh2) const
void read()
Fills all required arrays with the values read from the tables.
double THoEX_gg_phi3_mumu_CMS13
double ip_ex_gg_phi_tautau_CMS8(double mass)
gslpp::complex B0_MZ2_mA2_0_mHh2(const double MZ2, const double mA2, const double mHh2) const
gslpp::matrix< double > csrHpHm_pp_13_NLO
double ip_ex_pp_phi_phi1phi1_bbWW_qqlnu_CMS13_cache[2][CacheSize]
double ip_csr_ggH_b_13_cache[2][CacheSize]
double ip_low_pp_phi2_gaga_CMS13_cache[2][CacheSize]
double ip_ex_pp_phi_phi1phi1_bbtautau_2_CMS13(double mass)
double beta_h_t
beta function evaluated at Mt and m1_2
double gg_phi2_phi1phi1_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double pp_phi2_phi1phi1_with_channel_bbWW_qqlnu_TH13
double THoEX_tt_phi3_tt_CMS13
double ip_low_t_Hpb_cbb_CMS8_cache[2][CacheSize]
double THoEX_phi2Z_bbZ_LEP209
gslpp::matrix< double > br_cc
double ip_ex_gg_phi_Zga_CMS13_cache[2][CacheSize]
double THoEX_gg_h_phi2phi2_mumumumu_ATLAS13
gslpp::matrix< double > ATLAS13_t_Hpb_cbb
ATLAS observed upper branching fraction limits at 13 TeV, depending on the charged scalar mass.
gslpp::matrix< double > log_cs_ZH_8
double ip_ex_VV_phi_ZZ_qqllnunu_ATLAS13(double mass)
gslpp::complex B0_MZ2_0_mHp2_mHl2_cache[4][CacheSize]
double ip_ex_bb_phi_phi1Z_bbZ_ATLAS13_cache[2][CacheSize]
double ip_low_pp_h_phi23phi23_bbmumu_ATLAS13_old_cache[2][CacheSize]
double THoEX_pp_phi2_ZZ_llqqnunull_CMS13
double ggVV_phi2_WW_lnulnu_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double gg_phi3_tautau_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double THoEX_pp_phi3_Zga_llga_CMS8
double bb_phi3_tautau_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::complex B00_MZ2_0_MW2_mHh2(const double MZ2, const double MW2, const double mHh2) const
double ip_ex_bb_phi_tautau_CMS13(double mass)
double ip_ex_VV_phi_WW_ATLAS8(double mass)
double yu2R
Coupling of H to leptons.
double THoEX_gg_phi3_phi1Z_bbZ_2_CMS13
double pp_phi3_VV_TH8
Signal strength for the process with $VV=WW,ZZ$ at the LHC with 8 TeV.
double ip_ex_gg_phi_WW_lnuqq_ATLAS13_cache[2][CacheSize]
double ip_ex_bb_phi_mumu_CMS13_cache[2][CacheSize]
double THoEX_pp_h_phi3phi3_mumumumu_CMS13
double ip_csr_ggH_b_8_cache[2][CacheSize]
double THoEX_pp_Hp_tb_CMS8
double gg_phi3_phi1Z_bbll_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_VV_phi_WW_lnuqq_ATLAS13_cache[2][CacheSize]
double ip_low_pp_h_phi3phi3_bbtautau_CMS13_cache[2][CacheSize]
double Gamma_h
Total decay rate in the GTHDM.
gslpp::complex B0p_MZ2_mHl2_0_mHp2_cache[4][CacheSize]
double THoEX_gg_phi2_tautau_ATLAS8
gslpp::complex B00_MZ2_MW2_MW2_mHl2(const double MZ2, const double MW2, const double mHl2) const
gslpp::complex B0p_MZ2_0_mHp2_mA2(const double MZ2, const double mHp2, const double mA2) const
double ip_low_phi2Z_tautauZ_LEP209_cache[2][CacheSize]
gslpp::complex B0p_MZ2_0_mA2_mHh2_cache[4][CacheSize]
double THoEX_pp_h_phi2Z_mumull_CMS13
double THoEX_pp_phi3_phi1phi1_bbgaga_CMS8
double ip_low_phiEphi3_bbbb_LEP209_cache[3][CacheSize]
double gg_phi2_phi1Z_tautaull_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > OPAL209_HpHm_taunuAWst_taunubb
double THoEX_pp_phi3_phi1phi1_bbtautau_1_CMS13
gslpp::complex B0p_MZ2_mHh2_0_mHp2_cache[4][CacheSize]
double gg_phi3_phi1Z_tautaull_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::complex B0p_MZ2_mHh2_mA2_mA2_cache[4][CacheSize]
double THoEX_gg_phi3_tautau_ATLAS13
double ip_ex_pp_Hpm_tb_ATLAS13(double mass)
gslpp::matrix< double > CMS13_bb_phi_tautau
double ip_ex_gg_phi_ZZ_llllnunu_ATLAS13(double mass)
gslpp::complex B00_MZ2_0_mHh2_mHp2(const double MZ2, const double mHh2, const double mHp2) const
double ip_susy_pp_HpHm_taunutaunu_ATLAS13_cache[2][CacheSize]
double pp_phi2_phi1phi1_bbgaga_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double THoEX_gg_phi3_Zga_qqga_ATLAS13
double bb_phi2_tautau_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double THoEX_gg_phi2_ZZ_ATLAS8
double THoEX_HpHm_qqqq_LEP209
void CacheShift(gslpp::complex cache[][CacheSize], const int NumPar, const double params[], const gslpp::complex newResult) const
Adds a new result and its parameters into the cache.
double THoEX_bb_phi2_bb_ATLAS13
double THoEX_Zh_h_phi2phi2_bbbb_ATLAS13
double THoEX_gg_phi2_phi1Z_bbZ_ATLAS8
gslpp::matrix< double > CMS8_pp_phi2_gaga
double THoEX_pp_Hpm_taunu_ATLAS13
double pp_phi2_phi1phi1_with_channel_bbbb_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double bb_phi3_phi1Z_with_channel_bbZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV. The channel used is .
gslpp::complex B0_MZ2_mHp2_0_mHh2(const double MZ2, const double mHp2, const double mHh2) const
gslpp::matrix< double > br_mumu
double ip_ex_pp_Hp_tb_CMS8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_t_Hpb_WAb_Wmumub
ATLAS observed upper branching fraction limits at 13 TeV, depending on charged and pseudoscalar mass...
gslpp::matrix< double > ATLAS13_tt_phi_tt
ATLAS observed upper cross section limits at 13 TeV, depending on the Higgs mass.
double THoEX_bb_phi2_bb_CMS8
double interpolateNU(gslpp::matrix< double > arrayTab, double x)
Linearly interpolates a table with one parameter dimension and non-uniform intervals between bin valu...
double BrSM_htoZga
SM branching ratio of .
gslpp::complex B00_MZ2_0_mHp2_mHp2_cache[3][CacheSize]
double THoEX_pp_phi3_phi2Z_tautaull_1_CMS8
double gg_phi2_phi1phi1_gagaWW_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double yl2R
Coupling of H to down quarks real part.
double gg_phi3_phi1phi1_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double BrSM_htoZZ
SM branching ratio of .
gslpp::complex I_H_W_cache[3][CacheSize]
double ggF_VBF_Vh13
Ratio of GTHDM and SM cross sections for the production of h through ggF, VBF and Vh at 13 TeV.
double ggF_tth13
Ratio of GTHDM and SM cross sections for ggF and tth production of h at 13 TeV.
double ip_low_pp_h_phi3phi3_gagagaga_ATLAS13(double mass)
gslpp::complex B00_MZ2_MW2_MW2_mHh2(const double MZ2, const double MW2, const double mHh2) const
double lambdaijk(const double R1i, const double R2i, const double R3i, const double R1j, const double R2j, const double R3j, const double R1k, const double R2k, const double R3k, const double lambda1H, const double lambda3H, const double lambda4H, const double Relambda5H, const double Imlambda5H, const double Relambda6H, const double Imlambda6H, const double Relambda7H, const double Imlambda7H) const
lambdaijk function
double ip_ex_tt_phi2_tt_CMS13_cache[2][CacheSize]
gslpp::matrix< double > CMS8_bb_phi_mumu
double ip_low_t_Hpb_csb_CMS8(double mass)
double ip_ex_gg_phi_Zga_qqga_ATLAS13_cache[2][CacheSize]
double ip_ex_gg_phi_tautau_CMS8_cache[2][CacheSize]
double ip_cs_pptobbA_13_cache[2][CacheSize]
double THoEX_pp_phi2_phi1phi1_bbgaga_ATLAS13
double THoEX_gg_phi2_Zga_qqga_ATLAS13
double ip_ex_gg_phi_WW_heavy_CMS13(double mass)
void computeHpquantities()
double rh_QlQlE
Squared relative coupling of to two leptons, even part.
gslpp::matrix< double > csrH_top_charm_13
Production cross section ratio tables at 13 TeV obtained with HIGLU 4.34, depending on the Higgs mass...
gslpp::matrix< double > ATLAS13_bb_phi_tautau
double ip_cs_ZtoZH_8_cache[2][CacheSize]
double ip_ex_pp_Hpm_taunu_ATLAS13_cache[2][CacheSize]
double THoEX_bb_phi2_bb_CMS13
double ip_low_t_Hpb_WAb_Wmumub_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > log_cs_ggA_8
CP-odd Higgs production cross section tables at 8 TeV obtained with HIGLU 4.34, depending on the Higg...
gslpp::complex B0_MZ2_mHl2_0_mHp2_cache[4][CacheSize]
double ip_ex_pp_phi_phi1phi1_bbWW_qqlnu_CMS13(double mass)
gslpp::matrix< double > ATLAS13_gg_phi_Zga_llga
double ip_GammaHPtotSM(double mass)
Interpolating function for the total SM Higgs decay width.
double THoEX_gg_phi2_mumu_ATLAS13
double Gamma_W_inv
decay rate in the GTHDM.
double THoEX_pp_h_phi3phi3_bbmumu_CMS13
gslpp::complex I_H_Hp_cache[3][CacheSize]
double ip_ex_tt_phi2_tt_CMS13(double mass)
gslpp::complex B00_MZ2_0_mHl2_mHp2(const double MZ2, const double mHl2, const double mHp2) const
double ip_cs_pptottH_13(double mass)
Interpolating function for the top associated H production cross section at 13 TeV.
double pp_phi2_phi3Z_tautaull_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double THoEX_bb_phi3_tautau_CMS13
gslpp::complex B0_MZ2_0_MZ2_mHh2_cache[3][CacheSize]
double pp_phi3_phi1phi1_with_channel_bbWWorbbtautau_TH13
gslpp::complex sigmad_ATHDM
double ip_ex_gg_phi_WW_enumunu_ATLAS13(double mass)
gslpp::complex B00_MZ2_0_MW2_mHh2_cache[4][CacheSize]
double ip_ex_gg_phi_Zga_qqga_ATLAS13(double mass)
gslpp::matrix< double > integral_x_1mx2_G_variable_set_0_log
double THoEX_HpHm_AWstAWst_bbbb_OPAL209
double THoEX_gg_phi3_phi1Z_bbZ_ATLAS8
double rh_QlQlO
Squared relative coupling of to two leptons, odd part.
gslpp::matrix< double > CMS8_pp_phi_VV
CMS observed upper signal strength limits at 8 TeV, depending on the Higgs mass.
gslpp::complex B0_MZ2_mHp2_mHp2_mHh2_cache[4][CacheSize]
double THoEX_gg_phi2_tautau_ATLAS13
gslpp::matrix< double > ATLAS13_gg_phi_tautau
gslpp::matrix< double > ATLAS13_gg_phi_ZZ_llllnunu
gslpp::complex B00_MZ2_MZ2_mHl2_mA2(const double MZ2, const double mHl2, const double mA2) const
gslpp::matrix< double > ATLAS13_gg_phi_Zga_qqga
gslpp::complex B0p_MZ2_mHl2_mHh2_mHl2(const double MZ2, const double mHl2, const double mHh2) const
double ip_ex_ggVV_phi_WW_lnulnu_CMS13(double mass)
double THoEX_pp_phi3_VV_qqqq_ATLAS13
double interpolate2Dsteps(gslpp::matrix< double > arrayTab, double x, double y, double dx, double dy)
Linearly interpolates a table with two parameter dimensions, where the number of rows with the same x...
gslpp::matrix< double > ATLAS13_VV_phi_VV_llqq
gslpp::matrix< double > log_cs_bbH_8
SM Higgs production cross section table at 8 TeV obtained with SusHi 1.5, depending on the Higgs mass...
double GTHDM_BR_h_gg
branching ratio to two gluons in the GTHDM.
double THoEX_pp_phi2_phi3Z_bbll_2_CMS8
gslpp::matrix< double > CMS8_pp_phii_phijZ_bbll_2
double THoEX_pp_phi3_phi1phi1_bbZZ_llnunu_CMS13
double ip_low_pp_ttphi3_ttmumu_CMS13_cache[2][CacheSize]
gslpp::complex B0_MZ2_mHp2_0_mHl2_cache[4][CacheSize]
double pph8
Ratio of GTHDM and SM cross sections for the production of h at 8 TeV.
gslpp::complex B00_MZ2_0_mHp2_mHp2(const double MZ2, const double mHp2) const
double ip_ex_pp_phi_WW_lnuqq_CMS13_cache[2][CacheSize]
double bb_phi2_phi1Z_bbZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_csr_HpHm_pp_13_NLO(double mass)
double ip_low_pp_ttphi3_ttmumu_ATLAS13_cache[2][CacheSize]
double ip_ex_bb_phi_mumu_CMS8(double mass)
double pp_phi3_Zga_llga_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_bb_phi_tautau_CMS8(double mass)
double ip_ex_gg_phi_phi1phi1_ATLAS8_cache[2][CacheSize]
double ip_Br_HPtott_cache[2][CacheSize]
gslpp::complex B0p_MZ2_mHp2_0_mA2_cache[4][CacheSize]
double THoEX_VV_phi2_WW_ATLAS8
double THoEX_ggVV_phi2_WW_lnulnu_CMS13
gslpp::matrix< double > CMS13_bb_phi_bb
double THoEX_VV_phi2_WW_CMS13
double THoEX_bb_phi3_tautau_CMS8
gslpp::complex B0p_MZ2_mHh2_0_mHp2(const double MZ2, const double mHh2, const double mHp2) const
gslpp::matrix< double > log_cs_VBF_13
double THoEX_gg_phi3_gaga_ATLAS8
double THoEX_pp_phi2_phi1phi1_bbbb_ATLAS13
double BrSM_htoWW
SM branching ratio of .
double ip_cs_VBFtoH_13_cache[2][CacheSize]
gslpp::complex B0_MZ2_MZ2_MZ2_mHh2(const double MZ2, const double mHh2) const
double rphi3_ggO
Squared relative coupling of to two gluons, odd part.
gslpp::matrix< double > ATLAS8_pp_Hpm_taunu
ATLAS observed upper cross section limits at 8 TeV, depending on the charged Higgs mass.
double gg_phi3_phi1Z_bbZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_gg_phi_mumu_CMS13_cache[2][CacheSize]
double ip_ex_gg_phi_WW_ATLAS8_cache[2][CacheSize]
double THoEX_pp_h_phi2phi2_bbmumu_ATLAS13_old
double ip_low_pp_h_phi23phi23_mumumumu_CMS13_cache[2][CacheSize]
gslpp::complex B0_MZ2_mHh2_mHp2_mHp2_cache[4][CacheSize]
double GTHDM_BR_h_HH
branching ratio to two bosons in the GTHDM.
double rh_QuQuE
Squared relative coupling of to two up quarks, even part.
double ip_low_t_Hpb_cbb_CMS8(double mass)
double THoEX_bb_phi3_bb_CMS8
double THoEX_pp_phi2_phi1phi1_bbWW_ATLAS13
gslpp::matrix< double > CMS13_pp_phi_phi1phi1_bbbb_2
double ip_ex_gg_phi_ZZ_qqllnunu_ATLAS13(double mass)
gslpp::matrix< double > ATLAS13_pp_h_phi23phi23_bbmumu_old
double ip_low_t_Hpb_WAb_Wmumub_ATLAS13(double mass)
double ip_ex_pp_phi_phi1phi1_4WOr2W2tauOr4tau_CMS13(double mass)
double rh_Zga
Squared relative coupling of to a Z boson and a photon.
double ip_ex_gg_phi_mumu_ATLAS13(double mass)
double ip_ex_pp_phi_ZZ_qqnunu_CMS13(double mass)
double ip_ex_VV_phi_ZZ_llllnunu_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > CMS13_pp_phi3_bb_light
gslpp::complex B0p_MZ2_mHp2_mHp2_mHl2(const double MZ2, const double mHp2, const double mHl2) const
double ip_ex_gg_phi_tautau_ATLAS8(double mass)
double ip_Br_HPtomumu_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_gg_phi_phi1phi1_gagaWW
double ip_cs_WtoWH_8_cache[2][CacheSize]
double THoEX_pp_phi2_phi1phi1_bbbb_CMS8
double ip_ex_bb_phi_tautau_ATLAS13(double mass)
gslpp::matrix< double > CMS13_pp_phi_phi1phi1_bbZZ_llnunu
gslpp::matrix< double > log_cs_ggH_8
SM Higgs production cross section tables at 8 TeV from the LHC Higgs Cross Section Working Group,...
double ip_low_pp_h_phi3phi3_bbtautau_CMS13(double mass)
gslpp::matrix< double > ATLAS13_pp_Hpm_tb
double THoEX_gg_phi3_mumu_ATLAS13
double ip_low_pp_ttphi3_tttautau_CMS13_cache[2][CacheSize]
gslpp::complex B0_MZ2_0_mHp2_mHh2(const double MZ2, const double mHp2, const double mHh2) const
gslpp::complex B0p_MZ2_0_0_mHh2(const double MZ2, const double mHh2) const
double ip_ex_pp_phi_WW_lnuqq_CMS13(double mass)
double rphi2_ggO
Squared relative coupling of to two gluons, odd part.
double THoEX_gg_phi2_phi1Z_tautaull_CMS13
double ip_low_Vh_h_phi3phi3_bbbb_CMS13_cache[2][CacheSize]
double ip_low_pp_h_phi3phi3_bbmumu_CMS13(double mass)
double ip_csr_HpHm_pp_13_NLO_cache[2][CacheSize]
double ip_low_phiEphi3_tautautautau_LEP209_cache[3][CacheSize]
double THoEX_pp_phi2_Zga_llga_CMS8
double ip_low_t_Hpb_csb_CMS13(double mass)
gslpp::complex B0p_MZ2_mHl2_0_mHp2(const double MZ2, const double mHl2, const double mHp2) const
gslpp::matrix< double > GammaHtot_SM
Total SM decay width (obtained with HDECAY 6.10), depending on the Higgs mass.
gslpp::complex B0_MZ2_0_MZ2_mHh2(const double MZ2, const double mHh2) const
double ip_cs_VBFtoH_8_cache[2][CacheSize]
gslpp::complex B00_MZ2_0_MW2_mHl2_cache[4][CacheSize]
double THoEX_pp_bbphi3_bbmumu_CMS8
double THoEX_bb_phi2_tautau_ATLAS8
gslpp::matrix< double > csrA_top_charm_8
gslpp::matrix< gslpp::complex > Mu_GTHDM
gslpp::matrix< double > ATLAS8_gg_phi_phi1Z_bbZ
double THoEX_pp_phi3_phi2Z_bbll_1_CMS8
double pp_phi2_phi1phi1_bbgaga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex B0p_MZ2_mHp2_mHp2_mHl2_cache[4][CacheSize]
double ip_ex_ggVV_phi_WW_lnulnu_CMS13_cache[2][CacheSize]
double THoEX_pp_phi2_phi1phi1_bbZZ_llnunu_CMS13
double ip_low_pp_h_phi23Z_mumull_CMS13(double mass)
double ip_ex_bb_phii_phijZ_bbZ_ATLAS13_cache[3][CacheSize]
double rphi2_VV
Squared relative coupling of to two massive vector bosons.
double gg_phi2_gaga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_csr_ggH_tc_8_cache[2][CacheSize]
double gg_phi2_phi3Z_bbZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double THoEX_pp_phi3_phi1phi1_bbtautau_CMS8
double THoEX_gg_phi3_phi1phi1_gagaWW_ATLAS13
double ip_low_phi2Z_bbZ_LEP209_cache[2][CacheSize]
gslpp::matrix< double > log_cs_ttA_13
double bb_phi2_bb_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > integral_perturb_ATHDM_values_log
gslpp::matrix< double > CMS13_pp_HpHm_munumunu
double gg_phi3_phi1Z_with_channel_bbZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV. The SM Higgs is decaying...
gslpp::complex B00_MZ2_MW2_mHh2_mHp2(const double MZ2, const double MW2, const double mHh2, const double mHp2) const
gslpp::complex g_func(const double x) const
g function for the Int1 function
gslpp::matrix< double > ATLAS13_gg_phi_VV_llqq
double ip_low_phi2Z_gagaZ_LEP209(double mass)
double ip_low_gg_h_phi3phi3_tautautautau_ATLAS8_cache[2][CacheSize]
double ip_cs_ggtoH_13(double mass)
Interpolating function for the H production cross section via gluon-gluon fusion at 13 TeV.
double interpolate(gslpp::matrix< double > arrayTab, double x)
Linearly interpolates a table with one parameter dimension.
double THoEX_gg_phi2_phi1phi1_bbtautau_CMS8
double gg_phi2_phi1Z_with_channel_bbZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV. The SM Higgs is decaying...
double rh_ggE
Squared relative coupling of to two gluons, even part.
double THoEX_gg_phi2_WW_enumunu_ATLAS13
double THoEX_pp_phi3_phi2Z_tautaull_2_CMS8
double ip_ex_pp_phi_phi1phi1_bbtautau_1_ATLAS13_cache[2][CacheSize]
double ip_ex_bb_phi_tautau_ATLAS8(double mass)
double VV_phi3_WW_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::complex B0_MZ2_mHh2_mHh2_mHh2(const double MZ2, const double mHh2) const
double ip_low_phi2Z_tautauZ_LEP209(double mass)
gslpp::complex B0_MZ2_mHp2_mHp2_mHl2_cache[4][CacheSize]
gslpp::matrix< double > LEP209_phi2Z_bbZ
LEP observed upper cross-section times branching fraction limits for GeV, depending on the scalar m...
double ip_low_phi2Z_bbZ_LEP209(double mass)
double THoEX_gg_phi2_WW_ATLAS8
gslpp::matrix< double > CMS8_pp_bbphi3_bbtautau
double VV_phi2_ZZ_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_pp_phi3_bb_light_CMS13_cache[2][CacheSize]
gslpp::complex B00_MZ2_MW2_mHh2_mHp2_cache[5][CacheSize]
gslpp::complex B00_MZ2_0_mHl2_mA2(const double MZ2, const double mHl2, const double mA2) const
double THoEX_gg_phi3_mumu_CMS8
gslpp::complex B0_MZ2_MW2_MW2_mHh2(const double MZ2, const double MW2, const double mHh2) const
double ip_ex_gg_phi_ZZ_qqllnunu_ATLAS13_cache[2][CacheSize]
gslpp::complex B00_MZ2_0_MZ2_mHl2_cache[3][CacheSize]
double THoEX_gg_phi2_gaga_CMS13
gslpp::complex I_h_U(const double mHl2, const double Mu, const double Mc, const double Mt) const
double rphi2_ggE
Squared relative coupling of to two gluons, even part.
double ip_ex_bb_phi_phi1Z_bbZ_2_CMS13_cache[2][CacheSize]
double THoEX_pp_phi3_phi1phi1_bbVV_CMS13
double bb_phi3_tt_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex Int2(const double tau, const double lambda) const
function for Z gamma coupling to h, H and A
double ip_ex_pp_phi_phi1phi1_bbtautau_2_CMS13_cache[2][CacheSize]
double yu1R
Coupling of the SM-Higgs to leptons.
gslpp::complex B0_MZ2_mHh2_mHl2_mHl2(const double MZ2, const double mHh2, const double mHl2) const
double pp_phi2_bb_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double tt_phi3_tt_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double pp_phi3_ZZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double bb_phi2_phi1Z_with_channel_bbZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV. The channel used is .
double gg_phi2_gaga_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_cs_pptottH_8_cache[2][CacheSize]
double ip_ex_pp_phi2_bb_light_CMS13_cache[2][CacheSize]
gslpp::complex B0_MZ2_mA2_mA2_mHh2_cache[4][CacheSize]
gslpp::complex B00_MZ2_MZ2_MZ2_mHh2(const double MZ2, const double mHh2) const
gslpp::matrix< double > LEP209_HpHm_qqqq
gslpp::complex B0_MZ2_mA2_0_mHl2_cache[4][CacheSize]
gslpp::matrix< double > ATLAS8_gg_phi_ZZ
double THoEX_t_Hpb_csb_CMS8
gslpp::complex B0p_MZ2_mHl2_mA2_mA2_cache[4][CacheSize]
const GeneralTHDM * myGTHDM
double pp_Hpm_tb_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_VV_phi_WW_enumunu_ATLAS13_cache[2][CacheSize]
double ip_cs_ggtoHp_13(double mHp, double logtb)
Interpolating function for the H+ production cross section from two gluons at 13 TeV.
double ip_ex_pp_phi_phi1phi1_bbWWorbbtautau_low_masses_CMS13_cache[2][CacheSize]
gslpp::complex B00_MZ2_0_mHh2_mA2_cache[4][CacheSize]
gslpp::matrix< double > integral_Hp_A_Wstar
double ip_ex_pp_Hpm_taunu_ATLAS8_cache[2][CacheSize]
gslpp::matrix< double > log_cs_ZH_13
gslpp::complex I_h_L_cache[5][CacheSize]
gslpp::complex B0_MZ2_mHh2_0_mHp2_cache[4][CacheSize]
gslpp::matrix< double > br_bb
double ip_csr_ggA_tc_8_cache[2][CacheSize]
double ip_cs_pptottH_13_cache[2][CacheSize]
double ip_low_HpHm_qqtaunu_OPAL172_cache[2][CacheSize]
double ip_csr_ggA_b_8_cache[2][CacheSize]
double ip_ex_bb_phi_tautau_ATLAS8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_VV_phi_WW_enumunu
double ip_ex_pp_phi_phi1phi1_bbgaga_CMS8(double mass)
double THoEX_t_Hpb_taunub_CMS8
gslpp::matrix< double > CMS8_pp_phii_phijZ_tautaull_1
double THoEX_gg_phi2_ZZ_qqllnunu_ATLAS13
gslpp::matrix< double > ATLAS8_VV_phi_WW
double THoEX_pp_phi2_VV_CMS8
double ip_ex_VV_phi_WW_heavy_CMS13_cache[2][CacheSize]
double pp_phi2_ZZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex B0p_MZ2_mA2_0_mHh2_cache[4][CacheSize]
gslpp::complex B0p_MZ2_mHp2_mHp2_mHh2(const double MZ2, const double mHp2, const double mHh2) const
double GTHDM_BR_h_HpHm
branching ratio to bosons in the GTHDM.
gslpp::matrix< double > CMS8_t_Hpb_cbb
double THoEX_gg_phi2_bb_CMS8
gslpp::complex B0_MZ2_0_0_mHl2(const double MZ2, const double mHl2) const
gslpp::matrix< double > CMS8_gg_phi_mumu
double THoEX_gg_phi3_Zga_CMS13
gslpp::complex B0_MZ2_mHl2_mHl2_mHl2(const double MZ2, const double mHl2) const
double ip_low_pp_h_phi3phi3_mumutautau_CMS13(double mass)
double pp_phi2_phi1phi1_with_channel_bbgaga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV. where the decay channel ...
gslpp::complex A_A_L_cache[6][CacheSize]
gslpp::complex B00(const double mu2, const double p2, const double m02, const double m12) const
.
gslpp::complex B0(const double mu2, const double p2, const double m02, const double m12) const
.
gslpp::complex B0p(const double muIR2, const double p2, const double m02, const double m12) const
.
A model class for the Standard Model.
parameter of the Higgs potential
An observable class for the quartic Higgs potential coupling .
An observable class for the quartic Higgs potential coupling combination .
An observable class for the quartic Higgs potential coupling .
An observable class for the quartic Higgs potential coupling .
master
1.9.2