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;
1979 return (arrayTab(row[0],2));
1981 return (arrayTab(row[0],2) + (arrayTab(row[1],2)-arrayTab(row[0],2)) * (y-arrayTab(row[0],1)) / dy);
1983 return (arrayTab(row[0],2) + (arrayTab(row[1],2)-arrayTab(row[0],2)) * (x-arrayTab(row[0],0)) / dx);
1985 aux.push_back(arrayTab(row[0],2) + (arrayTab(row[1],2)-arrayTab(row[0],2)) * (y-arrayTab(row[0],1)) / dy);
1986 aux.push_back(arrayTab(row[2],2) + (arrayTab(row[3],2)-arrayTab(row[2],2)) * (y-arrayTab(row[2],1)) / dy);
1987 return (aux[0] + (aux[1]-aux[0]) * (x-arrayTab(row[0],0)) / dx);
1989 aux.push_back(arrayTab(row[0],2) + (arrayTab(row[1],2)-arrayTab(row[0],2)) * (y-arrayTab(row[0],1)) / dy);
1990 aux.push_back(arrayTab(row[1],2) + (arrayTab(row[2],2)-arrayTab(row[1],2)) * (y-arrayTab(row[0],1)) / dy);
1991 return (aux[0] + (aux[1]-aux[0]) * (x-arrayTab(row[0],0)) / dx);
1993 return (arrayTab(row[0],2) + arrayTab(row[1],2)) / 2.;
1995 throw std::runtime_error(
"Review applicability of the mass ranges provided to GeneralTHDMcache::interpolate2Dsteps");
2005 std::stringstream br1,br2,br3,br4,br5,br6,br7;
2006 std::stringstream dw1;
2007 std::stringstream cs1,cs2,cs3,cs4,cs5,cs6,cs7,cs8,cs9;
2008 std::stringstream cs11,cs12,cs13,cs14,cs15,cs16,cs17,cs18,cs19;
2009 std::stringstream cs20,cs21;
2010 std::stringstream csr1,csr2,csr3,csr4;
2011 std::stringstream csr11,csr12,csr13,csr14;
2012 std::stringstream ex1m6,ex1m5,ex1m4,ex1m3,ex1m2,ex1m1;
2013 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;
2014 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,\
2015 ex45,ex46n1,ex46n2,ex46a,ex47,ex48,ex49,ex49p2,ex50a,ex50b,ex51,ex52m2,ex52m1,ex52,ex52p1,ex52p2,ex53,ex54,ex55,ex56;
2016 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,\
2018 std::stringstream lowC01,lowC02,lowC03,lowC04,lowC05,lowC06,lowC07,lowC08, lowC09, lowC10, lowC11, lowC12;
2019 std::stringstream lowA01, lowA02, lowA03, lowA04, lowA05, lowA06, lowA07, lowA08, lowA09, lowA10, lowA11, lowA12;
2020 std::stringstream lowA801, lowA802;
2021 std::stringstream lowC801, lowC802, lowC803, lowC804, lowC805, lowC806;
2022 std::stringstream low209a, low209b, low209c, low209d, low209e;
2023 std::stringstream lowHpC801, lowHpC802, lowHpC803, lowHpC1301, lowHpC1302;
2024 std::stringstream lowHpA801, lowHpA1301, lowHpA1302, lowHpA1303;
2025 std::stringstream lowHp209a, lowHp209b, lowHp209c, lowHp172a, lowHp209d, lowHp209e;
2026 std::stringstream thint01, thint02, thint03, thint04, thint05, thint06, thint07, thint08, thint09, thint10;
2027 std::stringstream csrslepLO, csrslepNLO, susyHpA01, susyHpA02, susyHpC01, susyHpC02;
2028 std::stringstream bsg1;
2030 std::cout<<
"reading tables"<<std::endl;
2033 std::stringstream path;
2034 path << getenv(
"HEPFITTABS") <<
"/THDM/tabs/";
2035 std::string tablepath=path.str();
2037 br1 << tablepath <<
"br1.dat";
2039 br2 << tablepath <<
"br2.dat";
2041 br3 << tablepath <<
"br3.dat";
2043 br4 << tablepath <<
"br4.dat";
2045 br5 << tablepath <<
"br5.dat";
2047 br6 << tablepath <<
"br6.dat";
2049 br7 << tablepath <<
"br7.dat";
2051 dw1 << tablepath <<
"dw1.dat";
2053 cs1 << tablepath <<
"log_cs_ggH_8.dat";
2055 cs11 << tablepath <<
"log_cs_ggH_13.dat";
2057 cs2 << tablepath <<
"log_cs_VBF_8.dat";
2059 cs12 << tablepath <<
"log_cs_VBF_13.dat";
2061 cs3 << tablepath <<
"log_cs_WH_8.dat";
2063 cs13 << tablepath <<
"log_cs_WH_13.dat";
2065 cs4 << tablepath <<
"log_cs_ZH_8.dat";
2067 cs14 << tablepath <<
"log_cs_ZH_13.dat";
2069 cs5 << tablepath <<
"log_cs_ttH_8.dat";
2071 cs15 << tablepath <<
"log_cs_ttH_13.dat";
2073 cs6 << tablepath <<
"log_cs_bbH_8.dat";
2075 cs16 << tablepath <<
"log_cs_bbH_13.dat";
2077 cs7 << tablepath <<
"log_cs_ggA_8.dat";
2079 cs17 << tablepath <<
"log_cs_ggA_13.dat";
2081 cs8 << tablepath <<
"log_cs_ttA_8.dat";
2083 cs18 << tablepath <<
"log_cs_ttA_13.dat";
2085 cs9 << tablepath <<
"log_cs_bbA_8.dat";
2087 cs19 << tablepath <<
"log_cs_bbA_13.dat";
2089 cs20 << tablepath <<
"log_cs_ggHp_8.dat";
2091 cs21 << tablepath <<
"log_cs_ggHp_13.dat";
2093 csr1 << tablepath <<
"csrH_top_charm_8.dat";
2095 csr11 << tablepath <<
"csrH_top_charm_13.dat";
2097 csr2 << tablepath <<
"csrH_bottom_8.dat";
2099 csr12 << tablepath <<
"csrH_bottom_13.dat";
2101 csr3 << tablepath <<
"csrA_top_charm_8.dat";
2103 csr13 << tablepath <<
"csrA_top_charm_13.dat";
2105 csr4 << tablepath <<
"csrA_bottom_8.dat";
2107 csr14 << tablepath <<
"csrA_bottom_13.dat";
2109 csrslepLO << tablepath <<
"csrHpHm_pp_slepton_13_LO.dat";
2111 csrslepNLO << tablepath <<
"csrHpHm_pp_slepton_13_NLONLL.dat";
2115 ex1m6<< tablepath <<
"150801437_9b.dat";
2117 ex1m5<< tablepath <<
"150801437_9a.dat";
2119 ex1m4<< tablepath <<
"190703152_6b.dat";
2121 ex1m3<< tablepath <<
"190703152_6a.dat";
2123 ex1m2<< tablepath <<
"190108144_4b.dat";
2125 ex1m1<< tablepath <<
"190108144_4a.dat";
2132 ex1 << tablepath <<
"14096064_a.dat";
2134 ex2 << tablepath <<
"14096064_b.dat";
2136 ex3 << tablepath <<
"14076583.dat";
2138 ex4 << tablepath <<
"14078150.dat";
2140 ex5 << tablepath <<
"150705930_a.dat";
2142 ex6 << tablepath <<
"150705930_b.dat";
2144 ex7 << tablepath <<
"150900389_a.dat";
2146 ex8 << tablepath <<
"150900389_b.dat";
2148 ex9 << tablepath <<
"150904670.dat";
2150 ex10 << tablepath <<
"150204478_b.dat";
2152 ex11 << tablepath <<
"150204478_a.dat";
2154 ex12 << tablepath <<
"150400936.dat";
2156 ex13 << tablepath <<
"150608329.dat";
2158 ex14 << tablepath <<
"180206149.dat";
2160 ex15 << tablepath <<
"CMS-PAS-HIG-14-029_a.dat";
2162 ex16 << tablepath <<
"CMS-PAS-HIG-14-029_b.dat";
2164 ex17 << tablepath <<
"CMS-PAS-HIG-16-014.dat";
2166 ex18 << tablepath <<
"150304114.dat";
2168 ex19 << tablepath <<
"160306896.dat";
2170 ex20 << tablepath <<
"151001181_a.dat";
2172 ex21 << tablepath <<
"170700350.dat";
2174 ex22 << tablepath <<
"150404710.dat";
2176 ex23 << tablepath <<
"151001181_b.dat";
2179 ex24 << tablepath <<
"160302991_a.dat";
2181 ex25 << tablepath <<
"160302991_b.dat";
2183 ex26 << tablepath <<
"160302991_c.dat";
2185 ex27 << tablepath <<
"160302991_d.dat";
2190 ex28m1 << tablepath <<
"190702749.dat";
2194 ex28 << tablepath <<
"ATLAS_CONF_2022_008.dat";
2196 ex29 << tablepath <<
"ATLAS-CONF-2016-104_b.dat";
2208 ex30 << tablepath <<
"200212223_2a.dat";
2210 ex31 << tablepath <<
"200212223_2b.dat";
2217 ex32 << tablepath <<
"210213405.dat";
2224 ex33 << tablepath <<
"230904364.dat";
2232 ex34 << tablepath <<
"230411962.dat";
2241 ex35 << tablepath <<
"200914791_4a.dat";
2243 ex36 << tablepath <<
"200914791_4b.dat";
2248 ex37 << tablepath <<
"170809638_a.dat";
2250 ex38 << tablepath <<
"170809638_b.dat";
2254 ex39m4 << tablepath <<
"210906055_7a.dat";
2256 ex39m3 << tablepath <<
"210906055_7b.dat";
2260 ex39m2 << tablepath <<
"191201594_6c.dat";
2262 ex39m1 << tablepath <<
"191201594_6d.dat";
2271 ex39 << tablepath <<
"ATLAS_CONF_2022_066_ggF.dat";
2273 ex40 << tablepath <<
"ATLAS_CONF_2022_066_VBF.dat";
2278 ex41 << tablepath <<
"171007235_a.dat";
2280 ex42 << tablepath <<
"171007235_b.dat";
2282 ex43 << tablepath <<
"170804445.dat";
2286 ex43p1 << tablepath <<
"200414636_12a.dat";
2288 ex43p2 << tablepath <<
"200414636_12b.dat";
2293 ex44 << tablepath <<
"220207288.dat";
2297 ex45 << tablepath <<
"211211876.dat";
2301 ex46n1 << tablepath <<
"ATLAS-CONF-2021-030.dat";
2303 ex46n2 << tablepath <<
"200714811.dat";
2305 ex46a << tablepath <<
"181104671.dat";
2307 ex47 << tablepath <<
"180708567.dat";
2313 ex48 << tablepath <<
"220700230_ggF.dat";
2315 ex49 << tablepath <<
"220700230_bb.dat";
2318 ex49p2 << tablepath <<
"191011634.dat";
2326 ex50a << tablepath <<
"201105639_ggF_decay_bb.dat";
2330 ex51 << tablepath <<
"201105639_bb_decay_bb.dat";
2333 ex50b << tablepath <<
"201105639_ggF_decay_WW.dat";
2337 ex52m2 << tablepath <<
"190806463_7a.dat";
2339 ex52m1 << tablepath <<
"190806463_7b.dat";
2343 ex52 << tablepath <<
"CMS-PAS-HIG-16-025.dat";
2347 ex52p1 << tablepath <<
"181011822_7a.dat";
2349 ex52p2 << tablepath <<
"181011822_8a.dat";
2354 ex53 << tablepath <<
"180512191.dat";
2371 ex54 << tablepath <<
"220802717_ggF.dat";
2373 ex55 << tablepath <<
"220802717_bb.dat";
2385 ex56 << tablepath <<
"240509320.dat";
2389 ex57 << tablepath <<
"171203143.dat";
2391 ex58 << tablepath <<
"180401939_a.dat";
2393 ex59 << tablepath <<
"180303838.dat";
2395 ex60 << tablepath <<
"CMS-PAS-HIG-16-023.dat";
2397 ex61 << tablepath <<
"180209407.dat";
2399 ex62 << tablepath <<
"180603548.dat";
2405 ex63 << tablepath <<
"CMS-PAS-B2G-20-004.dat";
2410 ex64 << tablepath <<
"231001643.dat";
2414 ex65 << tablepath <<
"170702909.dat";
2416 ex66 << tablepath <<
"180801365.dat";
2418 ex67 << tablepath <<
"170804188.dat";
2420 ex67p1 << tablepath <<
"220610268.dat";
2422 ex67p2 << tablepath <<
"190404193.dat";
2426 ex67p3 << tablepath <<
"200606391_bblljj.dat";
2429 ex67p4 << tablepath <<
"200606391_bbllnunu.dat";
2433 ex67p5 << tablepath <<
"211203161.dat";
2437 ex67p6 << tablepath <<
"240309430.dat";
2441 ex68 << tablepath <<
"CMS-PAS-HIG-18-005_a.dat";
2443 ex69 << tablepath <<
"180702826_a.dat";
2445 ex70 << tablepath <<
"CMS-PAS-HIG-18-005_b.dat";
2447 ex71 << tablepath <<
"180702826_b.dat";
2450 ex72 << tablepath <<
"14126663.dat";
2452 ex73 << tablepath <<
"151203704.dat";
2454 ex74 << tablepath <<
"150807774_a.dat";
2456 ex75 << tablepath <<
"150807774_b.dat";
2458 ex76 << tablepath <<
"180707915.dat";
2462 ex77 << tablepath <<
"210210076.dat";
2464 ex78 << tablepath <<
"200107763.dat";
2468 ex79 << tablepath <<
"190304560.dat";
2471 lowC01 << tablepath <<
"CMS-HIG-17-029_4e.dat";
2477 lowC02 << tablepath <<
"CMS-HIG-22-007_11d.dat";
2483 lowC03 << tablepath <<
"CMS-HIG-22-007_10.dat";
2486 lowC04 << tablepath <<
"CMS-HIG-19-007_4a.dat";
2489 lowC05 << tablepath <<
"CMS-HIG-19-007_5a.dat";
2492 lowC06 << tablepath <<
"CMS-HIG-21-003_6.dat";
2495 lowC07 << tablepath <<
"CMS-HIG-18-006_10.dat";
2498 lowC08 << tablepath <<
"CMS-PAS-HIG-20-002_5b.dat";
2501 lowC09 << tablepath <<
"CMS-HIG-17-014_4.dat";
2504 lowC10 << tablepath <<
"CMS-EXO-21-018_14d.dat";
2507 lowC11 << tablepath <<
"CMS-EXO-21-018_14f.dat";
2510 lowC12 << tablepath <<
"CMS-HIG-18-026_3c.dat";
2513 lowA01 << tablepath <<
"ATLAS_CERN-EP-2021-157_9.dat";
2516 lowA02 << tablepath <<
"ATLAS_CERN-EP-2021-193_14b.dat";
2519 lowA03 << tablepath <<
"ATLAS_CERN-EP-2021-193_17b.dat";
2522 lowA04 << tablepath <<
"ATLAS_CERN-EP-2018-128_9c.dat";
2525 lowA05 << tablepath <<
"ATLAS_CERN-EP-2020-067_10.dat";
2528 lowA06 << tablepath <<
"ATLAS_CERN-EP-2018-153_6a.dat";
2531 lowA07 << tablepath <<
"ATLAS_CERN-EP-2017-295_2.dat";
2534 lowA08 << tablepath <<
"ATLAS-CONF-2023-035_7b.dat";
2537 lowA09 << tablepath <<
"ATLAS_CERN-EP-2023-070_5a.dat";
2540 lowA10 << tablepath <<
"ATLAS_CERN-EP-2024-235_7.dat";
2543 lowA11 << tablepath <<
"ATLAS_CERN-EP-2023-202_5d.dat";
2546 lowA12 << tablepath <<
"ATLAS_CERN-EP-2024-164_10a.dat";
2549 lowA801 << tablepath <<
"ATLAS_CERN-PH-EP-2015-187_4b.dat";
2552 lowA802 << tablepath <<
"ATLAS_CERN-PH-EP-2015-057_6a.dat";
2555 lowC801 << tablepath <<
"CMS-HIG-16-015_4.dat";
2558 lowC802 << tablepath <<
"CMS-HIG-16-015_5.dat";
2561 lowC803 << tablepath <<
"CMS-HIG-16-015_6.dat";
2564 lowC804 << tablepath <<
"CMS-HIG-17-013_4a.dat";
2567 lowC805 << tablepath <<
"CMS-HIG-14-033_3.dat";
2570 lowC806 << tablepath <<
"CMS-HIG-15-009_5b.dat";
2573 low209a << tablepath <<
"LEP_LHWG-NOTE-2001-08_4.dat";
2576 low209b << tablepath <<
"LEP_CERN-PH-EP-2006-001_t14b.dat";
2579 low209c << tablepath <<
"LEP_CERN-PH-EP-2006-001_t14c.dat";
2582 low209d << tablepath <<
"LEP_CERN-PH-EP-2006-001_t18.dat";
2585 low209e << tablepath <<
"LEP_CERN-PH-EP-2006-001_t19.dat";
2588 lowHpC801 << tablepath <<
"CMS-HIG-13-035_t4.dat";
2591 lowHpC802 << tablepath <<
"CMS-HIG-14-023_t10.dat";
2594 lowHpC803 << tablepath <<
"CMS-HIG-16-030_3.dat";
2597 lowHpC1301 << tablepath <<
"CMS-HIG-18-020_2a.dat";
2600 lowHpC1302 << tablepath <<
"CMS-HIG-18-021_4c.dat";
2603 lowHpA801 << tablepath <<
"ATLAS_CERN-PH-EP-2014-274_7a.dat";
2606 lowHpA1301 << tablepath <<
"ATLAS_CERN-EP-2022-207_8.dat";
2609 lowHpA1302 << tablepath <<
"ATLAS_CERN-EP-2023-070_5b.dat";
2612 lowHpA1303 << tablepath <<
"ATLAS_CERN-EP-2024-185_6.dat";
2615 lowHp209a << tablepath <<
"LEP_CERN-PH-EP-2012-369_4a.dat";
2618 lowHp209b << tablepath <<
"LEP_CERN-PH-EP-2012-369_4d.dat";
2621 lowHp209c << tablepath <<
"LEP_CERN-PH-EP-2008-016_7b.dat";
2624 lowHp172a << tablepath <<
"LEP_CERN-PPE-97-168_3b_tnqq.dat";
2627 lowHp209d << tablepath <<
"LEP_CERN-PH-EP-2008-016_10a.dat";
2630 lowHp209e << tablepath <<
"LEP_CERN-PH-EP-2008-016_10b.dat";
2633 susyHpA01 << tablepath <<
"ATLAS_CERN-EP-2019-191_7b.dat";
2636 susyHpA02 << tablepath <<
"ATLAS_CERN-EP-2022-132_6a.dat";
2639 susyHpC01 << tablepath <<
"CMS-SUS-21-001_5a.dat";
2642 susyHpC02 << tablepath <<
"CMS-SUS-21-008_16a.dat";
2645 thint01 << tablepath <<
"integral_x2_1mx_G_values_log.dat";
2648 thint02 << tablepath <<
"integral_x2_1px_G_values_log.dat";
2651 thint03 << tablepath <<
"integral_x2_G_values_log.dat";
2654 thint04 << tablepath <<
"integral_x_1mx2_G_values_log.dat";
2657 thint05 << tablepath <<
"integral_x_1mx_1px_G_values_log.dat";
2660 thint06 << tablepath <<
"integral_x2_1mx_G_variable_set_1_values_log.dat";
2663 thint07 << tablepath <<
"integral_x2_G_variable_set_1_values_log.dat";
2666 thint08 << tablepath <<
"integral_x_1mx2_G_variable_set_0_values_log.dat";
2669 thint09 << tablepath <<
"integral_perturb_ATHDM_values_log.dat";
2672 thint10 << tablepath <<
"integral_Hp_A_Wstar.dat";
2677 bsg1 << tablepath <<
"bsgammatable.dat";
2685 double params[] = {mass};
2701 double params[] = {mass};
2717 double params[] = {mass};
2733 double params[] = {mass};
2749 double params[] = {mass};
2765 double params[] = {mass};
2781 double params[] = {mass};
2797 double params[] = {mass};
2813 double params[] = {mass};
2819 double newResult = 0.0;
2820 if (mass>=10. && mass <=2000.) {
2832 double params[] = {mass};
2838 double newResult = 0.0;
2839 if (mass>=10. && mass <=2000.) {
2851 double params[] = {mass};
2857 double newResult = 0.0;
2858 if (mass>=10. && mass <=2000.) {
2870 double params[] = {mass};
2876 double newResult = 0.0;
2877 if (mass>=10. && mass <=2000.) {
2889 double params[] = {mass};
2895 double newResult = 0.0;
2896 if (mass>=10. && mass <=2000.) {
2908 double params[] = {mass};
2914 double newResult = 0.0;
2915 if (mass>=10. && mass <=2000.) {
2927 double params[] = {mass};
2933 double newResult = 0.0;
2934 if (mass>=10. && mass <=2000.) {
2946 double params[] = {mass};
2952 double newResult = 0.0;
2953 if (mass>=10. && mass <=2000.) {
2965 double params[] = {mass};
2971 double newResult = 0.0;
2972 if (mass>=10. && mass <=2000.) {
2984 double params[] = {mass};
2990 double newResult = 0.0;
2991 if (mass>=10. && mass <=2000.) {
3003 double params[] = {mass};
3009 double newResult = 0.0;
3010 if (mass>=10. && mass <=2000.) {
3022 double params[] = {mass};
3028 double newResult = 0.0;
3029 if (mass>=10. && mass <=2000.) {
3041 double params[] = {mass};
3047 double newResult = 0.0;
3048 if (mass>=20. && mass <=2000.) {
3060 double params[] = {mass};
3066 double newResult = 0.0;
3067 if (mass>=20. && mass <=2000.) {
3079 double params[] = {mass};
3085 double newResult = 0.0;
3086 if (mass>=20. && mass <=2000.) {
3098 double params[] = {mass};
3104 double newResult = 0.0;
3105 if (mass>=20. && mass <=2000.) {
3117 double params[] = {mass};
3123 double newResult = 0.0;
3124 if (mass>=20. && mass <=2000.) {
3136 double params[] = {mass};
3142 double newResult = 0.0;
3143 if (mass>=20. && mass <=2000.) {
3155 double params[] = {
mHp, logtb};
3161 double newResult = 0.0;
3162 if (
mHp>=180. && mHp <=1400. && logtb>=-1. && logtb<=1.75) {
3174 double params[] = {
mHp, logtb};
3180 double newResult = 0.0;
3181 if (
mHp>=180. && mHp <=2000. && logtb>=-1. && logtb<=1.75) {
3194 double params[] = {mass};
3210 double params[] = {mass};
3226 double params[] = {mass};
3242 double params[] = {mass};
3258 double params[] = {mass};
3274 double params[] = {mass};
3290 double params[] = {mass};
3306 double params[] = {mass};
3321 double params[] = {mass};
3336 double params[] = {mass};
3352 double params[] = {mass};
3369 double params[] = {mass};
3383 double params[] = {mass};
3397 double params[] = {mass};
3411 double params[] = {mass};
3427 double params[] = {mass};
3442 double params[] = {mass};
3462 double params[] = {mass};
3483 double params[] = {mass};
3500 double params[] = {mass};
3531 double params[] = {mass};
3546 double params[] = {mass};
3561 double params[] = {mass};
3577 double params[] = {mass};
3592 double params[] = {mass};
3608 double params[] = {mass};
3636 double params[] = {mass};
3650 double params[] = {mass};
3664 double params[] = {mass};
3678 double params[] = {mass};
3692 double params[] = {mass};
3706 double params[] = {mass};
3720 double params[] = {mass};
3734 double params[] = {mass};
3748 double params[] = {mass};
3762 double params[] = {mass};
3776 double params[] = {mass};
3790 double params[] = {mass};
3804 double params[] = {mass};
3818 double params[] = {mass};
3832 double params[] = {mass};
3846 double params[] = {mass};
3860 double params[] = {mass};
3874 double params[] = {mass};
3888 double params[] = {mass};
3902 double params[] = {mass};
3916 double params[] = {mass};
3930 double params[] = {mass};
3944 double params[] = {mass};
3958 double params[] = {mass};
3972 double params[] = {mass};
3986 double params[] = {mass};
4000 double params[] = {mass};
4019 double params[] = {mass};
4037 double params[] = {mass};
4061 double params[] = {mass};
4079 double params[] = {mass};
4110 double params[] = {mass};
4124 double params[] = {mass};
4138 double params[] = {mass};
4152 double params[] = {mass};
4166 double params[] = {mass};
4180 double params[] = {mass};
4194 double params[] = {mass};
4208 double params[] = {mass};
4228 double params[] = {mass};
4243 double params[] = {mass};
4264 double params[] = {mass};
4278 double params[] = {mass};
4292 double params[] = {mass};
4306 double params[] = {mass};
4320 double params[] = {mass};
4334 double params[] = {mass};
4348 double params[] = {mass};
4362 double params[] = {mass};
4376 double params[] = {mass};
4390 double params[] = {mass};
4423 double params[] = {mass};
4438 double params[] = {mass};
4452 double params[] = {mass};
4466 double params[] = {mass};
4480 double params[] = {mass};
4496 double params[] = {mass};
4514 double params[] = {mass};
4544 double params[] = {mass};
4562 double params[] = {mass};
4580 double params[] = {mass};
4596 double params[] = {mass};
4613 double params[] = {mass};
4634 double params[] = {mass};
4648 double params[] = {mass};
4662 double params[] = {mass};
4676 double params[] = {mass};
4690 double params[] = {mass};
4704 double params[] = {mass};
4718 double params[] = {mass};
4732 double params[] = {mass};
4746 double params[] = {mass};
4763 double params[] = {mass};
4783 double params[] = {mass};
4797 double params[] = {mass};
4811 double params[] = {mi, mj};
4827 double params[] = {mi, mj};
4842 double params[] = {mi, mj};
4858 double params[] = {mi, mj};
4874 double params[] = {mj, mi};
4890 double params[] = {mj, mi};
4907 double params[] = {mj, mi};
4924 double params[] = {mass};
4940 double params[] = {mass};
4956 double params[] = {mass};
4972 double params[] = {mass};
4988 double params[] = {mass};
5004 double params[] = {mass};
5018 double params[] = {mass};
5034 double params[] = {mass};
5049 double params[] = {mass};
5063 double params[] = {mass};
5077 double params[] = {mass};
5091 double params[] = {mass};
5105 double params[] = {mass};
5119 double params[] = {mass};
5133 double params[] = {mass};
5147 double params[] = {mass};
5161 double params[] = {mass};
5175 double params[] = {mass};
5189 double params[] = {mass};
5202 double params[] = {mass};
5215 double params[] = {mass};
5228 double params[] = {mass};
5241 double params[] = {mass};
5254 double params[] = {mass};
5267 double params[] = {mass};
5280 double params[] = {mass};
5293 double params[] = {mass};
5306 double params[] = {mass};
5319 double params[] = {mass};
5332 double params[] = {mass};
5345 double params[] = {mass};
5358 double params[] = {mass};
5371 double params[] = {mass};
5385 double params[] = {mass};
5398 double params[] = {mass};
5411 double params[] = {mass};
5424 double params[] = {mass};
5437 double params[] = {mass};
5450 double params[] = {mass};
5463 double params[] = {mass};
5476 double params[] = {mass};
5489 double params[] = {mass};
5502 double params[] = {mass};
5515 double params[] = {mi, mj};
5528 double params[] = {mi, mj};
5541 double params[] = {mass};
5554 double params[] = {mass};
5567 double params[] = {mass};
5580 double params[] = {mass};
5593 double params[] = {mass};
5606 double params[] = {mass};
5619 double params[] = {mass};
5632 double params[] = {mass};
5645 double params[] = {mass};
5658 double params[] = {mass};
5671 double params[] = {mass};
5684 double params[] = {mass};
5697 double params[] = {mass};
5710 double params[] = {mi, mj};
5723 double params[] = {mi, mj};
5736 double params[] = {mass};
5749 double params[] = {mass};
5762 double params[] = {mass};
5775 double params[] = {mass};
5792 if(mass >= 45. && mass <= 65.)
5794 else if(mass > 65. && mass <= 76.)
5796 else if(mass > 76. && mass <= 86.)
5798 else if(mass > 86. && mass <= 89.)
5800 else if(mass > 89. && mass <= 96.)
5802 else if(mass > 96. && mass <= 97.)
5804 else if(mass > 97. && mass <= 98.)
5806 else if(mass > 98. && mass <= 99.)
5808 else if(mass > 99. && mass <= 101.)
5810 else if(mass > 101. && mass <= 102.)
5823 double params[] = {wa, wb};
5828 double log10wa = log10(wa);
5829 double log10wb = log10(wb);
5839 double params[] = {wa, wb};
5844 double log10wa = log10(wa);
5845 double log10wb = log10(wb);
5856 double params[] = {wa, wb};
5861 double log10wa = log10(wa);
5862 double log10wb = log10(wb);
5871 double params[] = {wa, wb};
5876 double log10wa = log10(wa);
5877 double log10wb = log10(wb);
5887 double params[] = {wa, wb};
5892 double log10wa = log10(wa);
5893 double log10wb = log10(wb);
5903 double params[] = {wb};
5908 double log10wb = log10(wb);
5918 double params[] = {wb};
5923 double log10wb = log10(wb);
5933 double params[] = {wb};
5938 double log10wb = log10(wb);
5948 double params[] = {ratio};
5953 double log10ratio = log10(ratio);
5963 double params[] = {mi, mj};
5978 double params[] = {logtb, logmHp};
5997 gslpp::complex z = -gslpp::complex::i()*M_PI;
5998 return -pow(log((1.0+sqrt(1.0-x))/(1.0-sqrt(1.0-x)))+z,2)/4.0;
6001 return pow(asin(sqrt(1.0/x)),2);
6007 gslpp::complex z = -gslpp::complex::i()*M_PI;
6008 gslpp::complex gs1 = sqrt(1.0-x)*(log((1.0+sqrt(1.0-x))/(1.0-sqrt(1.0-x)))+z)/2.0;
6012 gslpp::complex gg1 = sqrt(x-1.0)*asin(sqrt(1.0/x));
6018 return tau*lambda/(tau-lambda)/2.0+tau*tau*lambda*lambda/((tau-lambda)
6019 *(tau-lambda))/2.0*(
f_func(tau)-
f_func(lambda))+tau*tau*lambda/((tau-lambda)
6024 return -tau*lambda/(tau-lambda)/2.0*(
f_func(tau)-
f_func(lambda));
6032 double params[] = {mHl2,
Mu,
Mc,
Mt};
6038 double TAUu=4.0*
Mu*
Mu/mHl2;
6039 double TAUc=4.0*
Mc*
Mc/mHl2;
6040 double TAUt=4.0*
Mt*
Mt/mHl2;
6041 gslpp::complex newResult = -(8./3.)*(TAUu*(1.0+(1.0-TAUu)*
f_func(TAUu))
6042 +TAUc*(1.0+(1.0-TAUc)*
f_func(TAUc))+TAUt*(1.0+(1.0-TAUt)*
f_func(TAUt)));
6050 double params[] = {mHh2,
Mc,
Mt};
6056 double TAUc=4.0*
Mc*
Mc/mHh2;
6057 double TAUt=4.0*
Mt*
Mt/mHh2;
6058 gslpp::complex newResult = -(8./3.)*(TAUc*(1.0+(1.0-TAUc)*
f_func(TAUc))
6059 +TAUt*(1.0+(1.0-TAUt)*
f_func(TAUt)));
6067 double params[] = {mA2,
Mc,
Mt};
6073 double TAUc=4.0*
Mc*
Mc/mA2;
6074 double TAUt=4.0*
Mt*
Mt/mA2;
6075 gslpp::complex newResult = -(8./3.)*(TAUc*
f_func(TAUc)+TAUt*
f_func(TAUt));
6083 double params[] = {mHl2,
Md,
Ms,
Mb};
6089 double TAUd=4.0*
Md*
Md/mHl2;
6090 double TAUs=4.0*
Ms*
Ms/mHl2;
6091 double TAUb=4.0*
Mb*
Mb/mHl2;
6092 gslpp::complex newResult = -(2./3.)*(TAUd*(1.0+(1.0-TAUd)*
f_func(TAUd))
6093 +TAUs*(1.0+(1.0-TAUs)*
f_func(TAUs))+TAUb*(1.0+(1.0-TAUb)*
f_func(TAUb)));
6101 double params[] = {mHh2,
Ms,
Mb};
6107 double TAUs=4.0*
Ms*
Ms/mHh2;
6108 double TAUb=4.0*
Mb*
Mb/mHh2;
6109 gslpp::complex newResult = -(2./3.)*(TAUs*(1.0+(1.0-TAUs)*
f_func(TAUs))
6110 +TAUb*(1.0+(1.0-TAUb)*
f_func(TAUb)));
6118 double params[] = {mA2,
Ms,
Mb};
6124 double TAUs=4.0*
Ms*
Ms/mA2;
6125 double TAUb=4.0*
Mb*
Mb/mA2;
6126 gslpp::complex newResult = -(2./3.)*(TAUs*
f_func(TAUs)+TAUb*
f_func(TAUb));
6134 double params[] = {mHl2,
Me,
Mmu,
Mtau};
6140 double TAUe=4.0*
Me*
Me/mHl2;
6141 double TAUmu=4.0*
Mmu*
Mmu/mHl2;
6143 gslpp::complex newResult = -2.0*(TAUe*(1.0+(1.0-TAUe)*
f_func(TAUe))
6144 +TAUmu*(1.0+(1.0-TAUmu)*
f_func(TAUmu))
6145 +TAUtau*(1.0+(1.0-TAUtau)*
f_func(TAUtau)));
6153 double params[] = {mHh2,
Mmu,
Mtau};
6159 double TAUmu=4.0*
Mmu*
Mmu/mHh2;
6161 gslpp::complex newResult = -2.0*(TAUmu*(1.0+(1.0-TAUmu)*
f_func(TAUmu))+
6162 TAUtau*(1.0+(1.0-TAUtau)*
f_func(TAUtau)));
6170 double params[] = {mA2,
Mmu,
Mtau};
6176 double TAUmu=4.0*
Mmu*
Mmu/mA2;
6178 gslpp::complex newResult = -2.0*(TAUmu*
f_func(TAUmu)+TAUtau*
f_func(TAUtau));
6186 double params[] = {mH,
MW};
6192 double TAUw=4.0*
MW*
MW/(mH*mH);
6193 gslpp::complex newResult = 2.0 + 3.0*TAUw + 3.0*TAUw*(2.0-TAUw)*
f_func(TAUw);
6201 double params[] = {
mHp2, mH};
6207 double TAUhp=4.0*
mHp2/(mH*mH);
6208 gslpp::complex newResult = -TAUhp*(1.0-TAUhp*
f_func(TAUhp));
6214gslpp::complex
GeneralTHDMcache::A_h_U(
const double mHl2,
const double cW2,
const double Mu,
const double Mc,
const double Mt,
const double MZ)
const {
6222 double TAUu=4.0*
Mu*
Mu/mHl2;
6223 double TAUc=4.0*
Mc*
Mc/mHl2;
6224 double TAUt=4.0*
Mt*
Mt/mHl2;
6229 gslpp::complex newResult = -4.0*(1.0/2.0-4.0/3.0*
sW2)*(
Int1(TAUu,LAMu)+
Int1(TAUc,LAMc)
6239 double params[] = {mHh2,
cW2,
Mc,
Mt,
MZ};
6245 double TAUc=4.0*
Mc*
Mc/mHh2;
6246 double TAUt=4.0*
Mt*
Mt/mHh2;
6257 gslpp::complex newResult = -4.0*(1.0/2.0-4.0/3.0*
sW2)*(
Int1(TAUc,LAMc)-
Int2(TAUc,LAMc)
6264gslpp::complex
GeneralTHDMcache::A_A_U(
const double mA2,
const double cW2,
const double Mc,
const double Mt,
const double MZ)
const {
6266 double params[] = {mA2,
cW2,
Mc,
Mt,
MZ};
6272 double TAUc=4.0*
Mc*
Mc/mA2;
6273 double TAUt=4.0*
Mt*
Mt/mA2;
6277 gslpp::complex newResult = -4.0*(1.0/2.0-4.0/3.0*
sW2)*(-
Int2(TAUc,LAMc)-
Int2(TAUt,LAMt))/sqrt(
sW2*
cW2);
6283gslpp::complex
GeneralTHDMcache::A_h_D(
const double mHl2,
const double cW2,
const double Md,
const double Ms,
const double Mb,
const double MZ)
const {
6291 double TAUd=4.0*
Md*
Md/mHl2;
6292 double TAUs=4.0*
Ms*
Ms/mHl2;
6293 double TAUb=4.0*
Mb*
Mb/mHl2;
6298 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0/3.0*
sW2)*(
Int1(TAUd,LAMd)+
Int1(TAUs,LAMs)
6307 double params[] = {mHh2,
cW2,
Ms,
Mb,
MZ};
6313 double TAUs=4.0*
Ms*
Ms/mHh2;
6314 double TAUb=4.0*
Mb*
Mb/mHh2;
6318 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0/3.0*
sW2)*(
Int1(TAUs,LAMs)-
Int2(TAUs,LAMs)
6325gslpp::complex
GeneralTHDMcache::A_A_D(
const double mA2,
const double cW2,
const double Ms,
const double Mb,
const double MZ)
const {
6327 double params[] = {mA2,
cW2,
Ms,
Mb,
MZ};
6333 double TAUs=4.0*
Ms*
Ms/mA2;
6334 double TAUb=4.0*
Mb*
Mb/mA2;
6338 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0/3.0*
sW2)*(-
Int2(TAUs,LAMs)-
Int2(TAUb,LAMb))/sqrt(
sW2*
cW2);
6344gslpp::complex
GeneralTHDMcache::A_h_L(
const double mHl2,
const double cW2,
const double Me,
const double Mmu,
const double Mtau,
const double MZ)
const {
6352 double TAUe=4.0*
Me*
Me/mHl2;
6353 double TAUmu=4.0*
Mmu*
Mmu/mHl2;
6359 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0*
sW2)*(
Int1(TAUe,LAMe)+
Int1(TAUmu,LAMmu)
6361 -
Int2(TAUtau,LAMtau));
6367gslpp::complex
GeneralTHDMcache::A_HH_L(
const double mHh2,
const double cW2,
const double Mmu,
const double Mtau,
const double MZ)
const {
6375 double TAUmu=4.0*
Mmu*
Mmu/mHh2;
6380 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0*
sW2)*(
Int1(TAUmu,LAMmu)-
Int2(TAUmu,LAMmu)
6381 +
Int1(TAUtau,LAMtau)-
Int2(TAUtau,LAMtau));
6389gslpp::complex
GeneralTHDMcache::A_A_L(
const double mA2,
const double cW2,
const double Mmu,
const double Mtau,
const double MZ)
const {
6397 double TAUmu=4.0*
Mmu*
Mmu/mA2;
6402 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0*
sW2)*(-
Int2(TAUmu,LAMmu)-
Int2(TAUtau,LAMtau))/sqrt(
sW2*
cW2);
6410 double params[] = {mH,
cW2,
MW,
MZ};
6416 double TAUw=4.0*
MW*
MW/(mH*mH);
6419 gslpp::complex newResult = -sqrt(
cW2/
sW2)*(4.0*(3.0-
sW2/
cW2)*
Int2(TAUw,LAMw)
6420 +((1.0+2.0/TAUw)*
sW2/
cW2-(5.0+2.0/TAUw))*
Int1(TAUw,LAMw));
6434 double TAUhp=4.0*
mHp2/(mH*mH);
6437 gslpp::complex newResult = (1.0-2.0*
sW2)/sqrt(
cW2*
sW2)*
Int1(TAUhp,LAMhp);
6446 double s_LEP = sqrts * sqrts;
6447 double one_MZ2_s = 1. -
MZ*
MZ / s_LEP;
6449 double cWsW = sqrt(
cW2) * sqrt(
sW2);
6450 double ahat_e = -1. / 4. / cWsW;
6451 double vhat_e = (-1. + 4. *
sW2) / 4. / cWsW;
6452 double vhat_H = (-1. + 2. *
sW2) / 2. / cWsW;
6453 double beta_H = sqrt(std::fabs(1. - 4. * mass2 / s_LEP));
6455 return (2. *
GF*
GF *
MW*
MW*
MW*
MW *
sW2*
sW2 * (1. + 2. * vhat_e * vhat_H / one_MZ2_s +
6456 (ahat_e*ahat_e + vhat_e*vhat_e) * vhat_H*vhat_H / one_MZ2_s / one_MZ2_s) *
6457 beta_H*beta_H*beta_H / (3. * M_PI * s_LEP));
6469 double params[] = {a2, b2, c2};
6476 double newResult = 0.0;
6477 double x = (a2-b2-c2)*(a2-b2-c2)-4.0*b2*c2;
6478 if(x>0) newResult = sqrt(std::fabs(x/a2))/2.0;
6501 return sqrt(std::fabs(1.0 - 4.0*mf*mf / m_2));
6506 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
6508 return (1.0/2.0)*
vev*(Imlambda7H*(-Ri3*Rj3*Rk3 - Ri2*Rj2*Rk3) - 3.0*Imlambda6H*Ri1*Rj1*Rk3
6509 + lambda1H*Ri1*Rj1*Rk1 + Relambda7H*Ri2*Rj2*Rk2 + 3.0*Relambda6H*Ri1*Rj1*Rk2
6510 +(Relambda5H + lambda3H + lambda4H)*Ri1*Rj2*Rk2 - (2.0*Relambda5H - lambda3H - lambda4H)*Ri1*Rj3*Rk3
6511 + Relambda7H*Ri2*Rj3*Rk3 - Imlambda5H*Ri1*Rj2*Rk3);
6524 throw std::runtime_error(
"Signal strengths are only available in the A2HDM.");
6528 gslpp::complex i = gslpp::complex::i();
6549 double SigmaggF8 =
myGTHDM->computeSigmaggH(8.0);
6551 double Sigmaggh_tt8 =
myGTHDM->computeSigmaggH_tt(8.0);
6555 double Sigmaggh_bb8 =
myGTHDM->computeSigmaggH_bb(8.0);
6559 double Sigmatth8 =
myGTHDM->computeSigmattH(8.0);
6574 double Sigmatth13 =
myGTHDM->computeSigmattH(13.0);
6578 double SigmaVBFVh13 = (
myGTHDM->computeSigmaVBF(13.0)+
myGTHDM->computeSigmaWH(13.0)+
myGTHDM->computeSigmaZH(13.0));
6615 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);
6616 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;
6629 gslpp::complex I_h_W=
Rij_GTHDM(0,0)*I_hSM_W;
6632 gslpp::complex I_hSM_F= fermU+fermD+fermL;
6633 gslpp::complex I_hE_F= yu1.real()*fermU+
yd1.real()*fermD+yl1.real()*fermL;
6649 gslpp::complex I_hO_F = yu1.imag()*I_h_Ux +
yd1.imag()*I_h_Dx + yl1.imag()*I_h_Lx;
6652 rh_gaga = ((I_hE_F+I_h_W+I_h_Hp).abs2()+ (I_hO_F).abs2())/(I_hSM_F +I_hSM_W).abs2();
6674 gslpp::complex A_h_W =
Rij_GTHDM(0,0)*A_hSM_W;
6676 gslpp::complex A_hSM_F = (A_hE_Ux+ A_hE_Dx+ A_hE_Lx)/sqrt(
sW2*
cW2);
6677 gslpp::complex A_hE_F = (yu1.real()*A_hE_Ux+
yd1.real()*A_hE_Dx+ yl1.real()*A_hE_Lx)/sqrt(
sW2*
cW2);
6687 gslpp::complex A_hO_F=yu1.imag()*A_hO_Ux +
yd1.imag()*A_hO_Dx + yl1.imag()*A_hO_Lx;
6690 rh_Zga = ((A_hE_F+A_h_W+A_h_Hp).abs2()+ A_hO_F.abs2())/(A_hSM_F +A_hSM_W ).abs2();
6715 double Gamma_hHH =
HSTheta(
m1 - 2.0*
m2)*sqrt(std::fabs(1 - 4*(
m2*
m2)/(
m1*
m1)))*(lambda122*lambda122)/(32.0*M_PI*
m1);
6723 double Gamma_hAA =
HSTheta(
m1 - 2.0*
m3)*sqrt(std::fabs(1 - 4*(
m3*
m3)/(
m1*
m1)))*(lambda133*lambda133)/(32.0*M_PI*
m1);
6725 double Gamma_hHpHm =
HSTheta(
m1-2.0*
mHp)*sqrt(std::fabs(1 - 4*(
mHp*
mHp)/(
m1*
m1)))*(lambda1pm*lambda1pm)/(16.0*M_PI*
m1);
6767 double Gamma_h_exp =
myGTHDM->computeGammaHTotal();
6781 Gamma_h_inv = Gamma_hHpHm + Gamma_hHH + Gamma_hAA + Gamma_hHZ + Gamma_hAZ;
6783 Gamma_h = (Gamma_htobb + Gamma_htoWW + Gamma_htoZZ + Gamma_htotautau + Gamma_htogaga +
6784 Gamma_htogg + Gamma_htoZga + Gamma_htocc +
Gamma_h_inv);
6809 throw std::runtime_error(
"Direct Searches are only aviable in the A2HDM.");
6813 gslpp::complex i = gslpp::complex::i();
6848 double rphi2_QuQuE= yu2.real()*yu2.real();
6849 double rphi2_QuQuO= yu2.imag()*yu2.imag();
6850 double rphi2_QdQdE=
yd2.real()*
yd2.real();
6851 double rphi2_QdQdO=
yd2.imag()*
yd2.imag();
6852 double rphi2_QlQlE= yl2.real()*yl2.real();
6853 double rphi2_QlQlO= yl2.imag()*yl2.imag();
6854 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;
6855 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;
6869 gslpp::complex I_phi2E_F= yu2.real()*I_HH2_Ux+
yd2.real()*I_HH2_Dx+yl2.real()*I_HH2_Lx;
6885 gslpp::complex I_phi2O_F = yu2.imag()*I_A2_Ux +
yd2.imag()*I_A2_Dx + yl2.imag()*I_A2_Lx;
6887 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()
6888 + (I_phi2O_F).abs2());
6899 gslpp::complex A_phi2E_F = (yu2.real()*A_HH2_Ux+
yd2.real()*A_HH2_Dx+ yl2.real()*A_HH2_Lx)/sqrt(
sW2*
cW2);
6913 gslpp::complex A_phi2O_F=yu2.imag()*A_A2_Ux +
yd2.imag()*A_A2_Dx + yl2.imag()*A_A2_Lx;
6918 + A_phi2O_F.abs2());
6928 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()
6952 if (
m2>=20. &&
m2 <=2000.) {
7007double lambda132 =
lambdaijk(
R11,
R12,
R13,
R31,
R32,
R33,
R21,
R22,
R23,
lambda1,
lambda3,
lambda4,
Relambda5,
7009 lambdaijk(
R21,
R22,
R23,
R11,
R12,
R13,
R31,
R32,
R33,
lambda1,
lambda3,
lambda4,
Relambda5,
7011 lambdaijk(
R21,
R22,
R23,
R31,
R32,
R33,
R11,
R12,
R13,
lambda1,
lambda3,
lambda4,
Relambda5,
7013 lambdaijk(
R11,
R12,
R13,
R21,
R22,
R23,
R32,
R32,
R32,
lambda1,
lambda3,
lambda4,
Relambda5,
7015 lambdaijk(
R31,
R32,
R33,
R11,
R12,
R13,
R21,
R22,
R23,
lambda1,
lambda3,
lambda4,
Relambda5,
7017 lambdaijk(
R31,
R32,
R33,
R21,
R22,
R23,
R11,
R12,
R13,
lambda1,
lambda3,
lambda4,
Relambda5,
7020double lambda332 =
lambdaijk(
R31,
R32,
R33,
R31,
R32,
R33,
R21,
R22,
R23,
lambda1,
lambda3,
lambda4,
Relambda5,
7022 lambdaijk(
R21,
R22,
R23,
R31,
R32,
R33,
R31,
R32,
R33,
lambda1,
lambda3,
lambda4,
Relambda5,
7024 lambdaijk(
R31,
R32,
R33,
R21,
R22,
R23,
R31,
R32,
R33,
lambda1,
lambda3,
lambda4,
Relambda5,
7029double 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);
7032double lambda112 = (2.0)*(
lambdaijk(
R11,
R12,
R13,
R21,
R22,
R23,
R11,
R12,
R13,
lambda1,
lambda3,
lambda4,
Relambda5,
7034 lambdaijk(
R11,
R12,
R13,
R11,
R12,
R13,
R21,
R22,
R23,
lambda1,
lambda3,
lambda4,
Relambda5,
7036 lambdaijk(
R21,
R22,
R23,
R11,
R12,
R13,
R11,
R12,
R13,
lambda1,
lambda3,
lambda4,
Relambda5,
7040double 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);
7057double 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);
7088 +Gammaphi2_bb+Gammaphi2_tautau+Gammaphi2_mumu
7089 +Gammaphi2_WW+Gammaphi2_ZZ+Gamma_phi2gaga
7090 +Gamma_phi2Zga+Gamma_phi2gg +Gammaphi2_phi1phi3
7091 +Gammaphi2_phi1phi1+Gammaphi2_phi3phi3+Gammaphi2_HpHm
7092 +Gammaphi2_phi1Z+Gammaphi2_phi3Z+Gammaphi2_HpW;
7156 throw std::runtime_error(
"Direct Searches are only aviable in the A2HDM.");
7160 gslpp::complex i = gslpp::complex::i();
7197 double rphi3_QuQuE= yu3.real()*yu3.real();
7198 double rphi3_QuQuO= yu3.imag()*yu3.imag();
7199 double rphi3_QdQdE=
yd3.real()*
yd3.real();
7200 double rphi3_QdQdO=
yd3.imag()*
yd3.imag();
7201 double rphi3_QlQlE= yl3.real()*yl3.real();
7202 double rphi3_QlQlO= yl3.imag()*yl3.imag();
7203 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;
7204 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;
7216 gslpp::complex I_phi3E_F= yu3.real()*I_HH_Ux+
yd3.real()*I_HH_Dx+yl3.real()*I_HH_Lx;
7232 gslpp::complex I_phi3O_F = yu3.imag()*I_A_Ux +
yd3.imag()*I_A_Dx + yl3.imag()*I_A_Lx;
7234 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());
7244 gslpp::complex A_phi3E_F = (yu3.real()*A_HH_Ux+
yd3.real()*A_HH_Dx+ yl3.real()*A_HH_Lx)/sqrt(
sW2*
cW2);
7255 gslpp::complex A_phi3O_F=yu3.imag()*A_A_Ux +
yd3.imag()*A_A_Dx + yl3.imag()*A_A_Lx;
7258 + A_phi3O_F.abs2());
7263 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()
7287 if (
m3>=20. &&
m3 <=2000.) {
7323double lambda123 =
lambdaijk(
R11,
R12,
R13,
R21,
R22,
R23,
R31,
R32,
R33,
lambda1,
lambda3,
lambda4,
Relambda5,
7325 lambdaijk(
R21,
R22,
R23,
R11,
R12,
R13,
R31,
R32,
R33,
lambda1,
lambda3,
lambda4,
Relambda5,
7327 lambdaijk(
R21,
R22,
R23,
R31,
R32,
R33,
R11,
R12,
R13,
lambda1,
lambda3,
lambda4,
Relambda5,
7329 lambdaijk(
R11,
R12,
R13,
R31,
R32,
R33,
R21,
R22,
R23,
lambda1,
lambda3,
lambda4,
Relambda5,
7331 lambdaijk(
R31,
R32,
R33,
R11,
R12,
R13,
R21,
R22,
R23,
lambda1,
lambda3,
lambda4,
Relambda5,
7333 lambdaijk(
R31,
R32,
R33,
R21,
R22,
R23,
R11,
R12,
R13,
lambda1,
lambda3,
lambda4,
Relambda5,
7335double lambda223 =
lambdaijk(
R21,
R22,
R23,
R21,
R22,
R23,
R31,
R32,
R33,
lambda1,
lambda3,
lambda4,
Relambda5,
7337 lambdaijk(
R21,
R22,
R23,
R31,
R32,
R33,
R21,
R22,
R23,
lambda1,
lambda3,
lambda4,
Relambda5,
7339 lambdaijk(
R31,
R32,
R33,
R21,
R22,
R23,
R21,
R22,
R23,
lambda1,
lambda3,
lambda4,
Relambda5,
7341double lambda113 =
lambdaijk(
R11,
R12,
R13,
R11,
R13,
R13,
R31,
R32,
R33,
lambda1,
lambda3,
lambda4,
Relambda5,
7343 lambdaijk(
R11,
R11,
R13,
R31,
R32,
R33,
R11,
R12,
R13,
lambda1,
lambda3,
lambda4,
Relambda5,
7345 lambdaijk(
R31,
R32,
R33,
R11,
R12,
R13,
R11,
R12,
R13,
lambda1,
lambda3,
lambda4,
Relambda5,
7351double 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);
7353double 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);
7355double 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);
7381 + Gammaphi3_tautau + Gammaphi3_mumu + Gammaphi3_WW
7382 + Gammaphi3_ZZ + Gamma_phi3gaga + Gamma_phi3Zga + Gamma_phi3gg
7383 + Gammaphi3_phi1phi1 + Gammaphi3_phi2phi2 + Gammaphi3_phi1phi2
7384 + Gammaphi3_HpHm + Gammaphi3_phi1Z + Gammaphi3_phi2Z
7432 throw std::runtime_error(
"Direct Searches are only available in the A2HDM.");
7436 gslpp::complex i = gslpp::complex::i();
7441 double Mcp =
Mc*(1.+4.*alsmuc/3./M_PI+alsmuc*alsmuc/M_PI/M_PI*(-1.0414*(1.-4.*
Ms/3.*
Mc)+13.4434));
7443 double Mtp =
myGTHDM->getMtpole();
7447 double Mc2 = Mcp*Mcp;
7448 double Mb2 = Mbp*Mbp;
7449 double Mt2 = Mtp*Mtp;
7450 double r2hot =
mHp*
mHp/Mtp/Mtp;
7451 double r2bot = Mbp*Mbp/Mtp/Mtp;
7452 double Vtb =
myGTHDM->getCKM().getV_tb().abs();
7453 double Vcb =
myGTHDM->getCKM().getV_cb().abs();
7454 double Vcs =
myGTHDM->getCKM().getV_cs().abs();
7471 double GammaHpcs = 0.;
7472 double GammaHpcb = 0.;
7473 double GammaHptb = 0.;
7474 double GammaHpHlW = 0.;
7475 double GammaHpphi2W = 0.;
7476 double GammaHpphi3W = 0.;
7477 double GammaHpphi3Wst = 0.;
7482 -
sd.abs2()*Mb2*(Mb2-
mHp2+Mt2)-
su.abs2()*Mt2*(Mb2-
mHp2+Mt2))
7489 -
sd.abs2()*Mb2*(Mb2-
mHp2+Mc2)-
su.abs2()*Mc2*(Mb2-
mHp2+Mc2))
7496 -
sd.abs2()*Ms2*(Ms2-
mHp2+Mc2)-
su.abs2()*Mc2*(Ms2-
mHp2+Mc2))
7529 if(
mHp >= 40.0 && mHp <= 93.0 && mH3 >= 12.0 &&
mH3 <= (
mHp - 3.))
7535 GammaHpcs + GammaHpcb + GammaHptb +
7536 GammaHpHlW + GammaHpphi2W + GammaHpphi3W + GammaHpphi3Wst;
7555 double GammatHpb = 0.;
7556 double pre_tHpb =
GF * Mt2 * Mtp * Vtb * Vtb / 4. / sqrt(2.) / M_PI;
7557 gslpp::complex aCD =
su -
sd * Mbp / Mtp;
7558 gslpp::complex bCD =
su +
sd * Mbp / Mtp;
7559 double aCD2 = (aCD * aCD.conjugate()).real();
7560 double bCD2 = (bCD * bCD.conjugate()).real();
7564 GammatHpb = pre_tHpb * ((1. - r2hot + r2bot) * (aCD2 + bCD2) / 2. + r2bot * (aCD2 - bCD2)) *
7565 sqrt(1. + r2hot*r2hot + r2bot*r2bot - 2. * (r2hot + r2bot + r2hot*r2bot)) / 2.;
7568 double GammatWb = 1.35;
7578 double MZ2 =
MZ *
MZ;
7579 double MW2 =
MW *
MW;
7580 double preZ =
GF *
MZ * MZ2 / 24. / sqrt(2.) / M_PI;
7581 double preW =
GF *
MW * MW2 / 24. / sqrt(2.) / M_PI;
7583 double gZAH_2 =
R22 *
R22;
7585 double gWHpH_2 =
R22 *
R22;
7586 double gWHpA_2 =
R33 *
R33;
7589 double func_ZHpHm = 1. - 4. *
mHp2 / MZ2;
7593 double Gamma_ZAH =
HSTheta(
MZ-
mH3-
mH2) * preZ * gZAH_2 * sqrt(std::fabs(func_ZAH)) * func_ZAH;
7594 double Gamma_ZHpHm =
HSTheta(
MZ-2.*
mHp) * preZ * gZHpHm_2 * sqrt(std::fabs(func_ZHpHm)) * func_ZHpHm;
7595 double Gamma_WHpH =
HSTheta(
MW-
mHp-
mH2) * preW * gWHpH_2 * sqrt(std::fabs(func_WHpH)) * func_WHpH;
7596 double Gamma_WHpA =
HSTheta(
MW-
mHp-
mH3) * preW * gWHpA_2 * sqrt(std::fabs(func_WHpA)) * func_WHpA;
7611 throw std::runtime_error(
"Direct Searches are only aviable in the A2HDM.");
7615 gslpp::complex i = gslpp::complex::i();
7617 double Br_Ztoee=0.03363;
7618 double Br_Ztomumu=0.03366;
7619 double Br_Ztotautau=0.0337;
7620 double Br_Ztoinv=0.2;
7622 double Br_Wtoenu=0.1071;
7623 double Br_Wtomunu=0.1063;
7624 double Br_Wtotaunu=0.1138;
7625 double Br_Wtohadrons=0.6741;
7626 double Br_tautoleptons=0.3521;
7627 double Br_tautohadrons=1-Br_tautoleptons;
8003 +
GTHDM_BR_h_ZZ*2.0*Br_Ztoinv*(Br_Ztoee+Br_Ztomumu+Br_Ztotautau*Br_tautoleptons*Br_tautoleptons));
8005 +
GTHDM_BR_h_ZZ*2.0*Br_Ztoinv*(Br_Ztoee+Br_Ztomumu+Br_Ztotautau*Br_tautoleptons*Br_tautoleptons));
8012 +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)
8015 +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)
8019 +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)
8022 +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)
8063 (
GTHDM_BR_h_WW*((Br_Wtoenu+Br_Wtomunu)*(Br_Wtoenu+Br_Wtomunu)
8064 +2.0*(Br_Wtoenu+Br_Wtomunu)*Br_Wtohadrons)
8067 (
BrSM_htoWW*((Br_Wtoenu+Br_Wtomunu)*(Br_Wtoenu+Br_Wtomunu)
8068 +2.0*(Br_Wtoenu+Br_Wtomunu)*Br_Wtohadrons)
8073 (
GTHDM_BR_h_WW*((Br_Wtoenu+Br_Wtomunu)*(Br_Wtoenu+Br_Wtomunu)
8074 +2.0*(Br_Wtoenu+Br_Wtomunu)*Br_Wtohadrons)
8077 (
BrSM_htoWW*((Br_Wtoenu+Br_Wtomunu)*(Br_Wtoenu+Br_Wtomunu)
8078 +2.0*(Br_Wtoenu+Br_Wtomunu)*Br_Wtohadrons)
8213 if(
m2>= 450.0 &&
m2<1400.0) {
8220 if(
m3>= 450.0 &&
m3<1400.0){
8230 if(
m2>= 400.0 &&
m2<1000.0) {
8237 if(
m3>= 400.0 &&
m3<1000.0){
8245 if(
m2>= 400.0 &&
m2<1000.0) {
8253 if(
m3>= 400.0 &&
m3<1000.0){
8262 if(
m2>= 100.0 &&
m2< 900.0) {
8270 if(
m3>= 100.0 &&
m3< 900.0) {
8277 if(
m2>= 330.0 &&
m2<1200.0) {
8284 if(
m3>= 330.0 &&
m3<1200.0) {
8298 if(
m2>= 350.0 &&
m2<650.0) {
8305 if(
m3>= 350.0 &&
m3<650.0) {
8320 if(
m2>= 550.0 &&
m2<1200.0) {
8327 if(
m3>= 550.0 &&
m3<1200.0) {
8341 if(
m2>= 50.0 &&
m2<350.0) {
8348 if(
m3>= 50.0 &&
m3<350.0) {
8364 if(
m2>= 300.0 &&
m2<1300.0) {
8371 if(
m3>= 300.0 &&
m3<1300.0) {
8388 if(
m2>= 120.0 &&
m2<500.0) {
8395 if(
m3>= 120.0 &&
m3<500.0) {
8404 if(
m2>= 120.0 &&
m2<500.0) {
8411 if(
m3>= 120.0 &&
m3<500.0) {
8431 if(
m2>= 140.0 &&
m2<1000.0) {
8444 if(
m3>= 140.0 &&
m3<1000.0) {
8453 if(
m2>= 140.0 &&
m2<1000.0) {
8467 if(
m3>= 140.0 &&
m3<1000.0) {
8487 if(
m2>= 200.0 &&
m2<1000.0) {
8497 if(
m3>= 200.0 &&
m3<1000.0) {
8506 if(
m2>= 200.0 &&
m2<1000.0) {
8513 if(
m3>= 200.0 &&
m3<1000.0) {
8527 if(
m2>= 90.0 &&
m2<1000.0) {
8535 if(
m3>= 90.0 &&
m3<1000.0) {
8546 if(
m2>= 90.0 &&
m2<1000.0) {
8553 if(
m3>= 90.0 &&
m3<1000.0) {
8560 if(
m2>= 90.0 &&
m2<1000.0) {
8567 if(
m3>= 90.0 &&
m3<1000.0) {
8574 if(
m2>= 90.0 &&
m2<1000.0) {
8581 if(
m3>= 90.0 &&
m3<1000.0) {
8588 if(
m2>= 200.0 &&
m2<2250.0) {
8599 if(
m3>= 200.0 &&
m3<2250.0) {
8609 if(
m2>= 90.0 &&
m2<3200.0) {
8619 if(
m3>= 90.0 &&
m3<3200.0) {
8626 if(
m2>= 200.0 &&
m2<2250.0) {
8635 if(
m3>= 200.0 &&
m3<2250.0) {
8642 if(
m2>= 60.0 &&
m2<3500.0) {
8649 if(
m3>= 60.0 &&
m3<3500.0) {
8658 if(
m2>= 65.0 &&
m2< 600.0) {
8665 if(
m3>= 65.0 &&
m3< 600.0) {
8672 if(
m2>= 160.0 &&
m2<3000.0) {
8684 if(
m3>= 160.0 &&
m3<3000.0)
8690 if(
m2>= 600.0 &&
m2<5000.0)
8696 if(
m3>= 600.0 &&
m3<5000.0)
8702 if(
m2>= 200.0 &&
m2<1600.0)
8709 if(
m3>= 200.0 &&
m3<1600.0)
8717 if(
m2>= 200.0 &&
m2<1200.0)
8725 if(
m3>= 200.0 &&
m3<1200.0)
8733 if(
m2>= 250.0 &&
m2<2400.0)
8741 if(
m3>= 250.0 &&
m3<2400.0)
8749 if(
m2>=1000.0 &&
m2<6800.0)
8757 if(
m3>=1000.0 &&
m3<6800.0)
8766 if(
m2>= 350.0 &&
m2<4000.0)
8774 if(
m3>= 350.0 &&
m3<4000.0)
8782 if(
m2>= 140.0 &&
m2<1000.0)
8790 if(
m3>= 140.0 &&
m3<1000.0)
8798 if(
m2>= 140.0 &&
m2<1000.0)
8807 if(
m3>= 140.0 &&
m3<1000.0)
8815 if(
m2>= 210.0 &&
m2<1200.0)
8823 if(
m3>= 210.0 &&
m3<1200.0)
8831 if(
m2>= 210.0 &&
m2<2000.0)
8839 if(
m3>= 210.0 &&
m3<2000.0)
8847 if(
m2>= 300.0 &&
m2<3000.0)
8855 if(
m3>= 300.0 &&
m3<3000.0)
8863 if(
m2>= 300.0 &&
m2<3000.0)
8871 if(
m3>= 300.0 &&
m3<3000.0)
8880 if(
m2>= 130.0 &&
m2<3000.0)
8888 if(
m3>= 130.0 &&
m3<3000.0)
8896 if(
m2>=1000.0 &&
m2<4000.0)
8904 if(
m3>=1000.0 &&
m3<4000.0)
8912 if(
m2>= 300.0 &&
m2<1500.0)
8920 if(
m3>= 300.0 &&
m3<1500.0)
8928 if(
m2>= 300.0 &&
m2<1500.0)
8936 if(
m3>= 300.0 &&
m3<1500.0)
8947 if(
m2>= 200.0 &&
m2<3000.0)
8955 if(
m3>= 200.0 &&
m3<3000.0)
8963 if(
m2>= 200.0 &&
m2<3000.0)
8971 if(
m3>= 200.0 &&
m3<3000.0)
8989 if(
m2>= 1000.0 &&
m2<4500.0)
8997 if(
m3>= 1000.0 &&
m3<4500.0)
9005 if(
m2>= 1000.0 &&
m2<4500.0)
9013 if(
m3>= 1000.0 &&
m3<4500.0)
9033 if(
m2>= 250.0 &&
m2<4000.0)
9041 if(
m3>= 250.0 &&
m3<4000.0)
9049 if(
m2>= 250.0 &&
m2<3000.0)
9057 if(
m3>= 250.0 &&
m3<3000.0)
9065 if(
m2>= 200.0 &&
m2<1000.0)
9073 if(
m3>= 200.0 &&
m3<1000.0)
9081 if(
m2>= 300.0 &&
m2<3000.0)
9089 if(
m3>= 300.0 &&
m3<3000.0)
9097 if(
m2>= 300.0 &&
m2<3000.0)
9105 if(
m3>= 300.0 &&
m3<3000.0)
9113 if(
m2>=1000.0 &&
m2<4400.0)
9121 if(
m3>=1000.0 &&
m3<4400.0)
9129 if(
m2>= 145.0 &&
m2<1000.0)
9135 if(
m3>= 145.0 &&
m3<1000.0)
9141 if(
m2>=1200.0 &&
m2<3000.0)
9148 if(
m3>=1200.0 &&
m3<3000.0)
9161 if(
m2>=300.0 &&
m2<5000.0)
9169 if(
m3>=300.0 &&
m3<5000.0)
9179 if(
m2>=300.0 &&
m2<5000.0)
9187 if(
m3>=300.0 &&
m3<5000.0)
9200 if(
m2>= 260.0 &&
m2<1000.0)
9207 if(
m3>= 260.0 &&
m3<1000.0)
9215 if(
m2>= 270.0 &&
m2<1100.0)
9223 if(
m3>= 270.0 &&
m3<1100.0)
9231 if(
m2>= 260.0 &&
m2<1100.0)
9239 if(
m3>= 260.0 &&
m3<1100.0)
9247 if(
m2>= 260.0 &&
m2< 350.0)
9255 if(
m3>= 260.0 &&
m3< 350.0)
9263 if(
m2>= 350.0 &&
m2<1000.0)
9271 if(
m3>= 350.0 &&
m3<1000.0)
9279 if(
m2>= 250.0 &&
m2<3000.0)
9287 if(
m3>= 250.0 &&
m3<3000.0)
9295 if(
m2>= 260.0 &&
m2<1200.0)
9303 if(
m3>= 260.0 &&
m3<1200.0)
9311 if(
m2>=1000.0 &&
m2<3000.0)
9319 if(
m3>=1000.0 &&
m3<3000.0)
9327 if(
m2>= 251.0 &&
m2<1000.0)
9335 if(
m3>= 251.0 &&
m3<1000.0)
9342 if(
m2>= 260.0 &&
m2< 1000.0)
9350 if(
m3>= 260.0 &&
m3< 1000.0)
9368 if(
m2>= 251.0 &&
m2<1600.0)
9376 if(
m3>= 251.0 &&
m3<1600.0)
9384 if(
m2>= 1000.0 &&
m2<3000.0)
9392 if(
m3>= 1000.0 &&
m3<3000.0)
9401 if(
m2>= 250.0 &&
m2< 900.0)
9409 if(
m3>= 250.0 &&
m3< 900.0)
9417 if(
m2>= 900.0 &&
m2<4000.0)
9425 if(
m3>= 900.0 &&
m3<4000.0)
9433 if(
m2>= 260.0 &&
m2< 900.0)
9441 if(
m3>= 260.0 &&
m3< 900.0)
9451 if(
m2>= 250.0 &&
m2< 1000.0)
9459 if(
m3>= 250.0 &&
m3< 1000.0)
9470 if(
m2>= 800.0 &&
m2< 3500.0)
9478 if(
m3>= 800.0 &&
m3< 3500.0)
9492 if(
m2>= 260.0 &&
m2< 1000.0)
9500 if(
m3>= 260.0 &&
m3< 1000.0)
9508 if(
m2>= 250.0 &&
m2< 1000.0)
9516 if(
m3>= 250.0 &&
m3< 1000.0)
9528 if(
m2>= 800.0 &&
m2< 4500.0)
9534 if(
m3>= 800.0 &&
m3< 4500.0)
9543 if(
m2>= 250.0 &&
m2 < 899.0)
9549 if(
m3>= 250.0 &&
m3< 899.0)
9557 if(
m2>= 500.0 &&
m2< 3000.0)
9565 if(
m3>= 500.0 &&
m3< 3000.0)
9573 if(
m2>= 260.0 &&
m2< 500.0)
9581 if(
m3>= 260.0 &&
m3< 500.0)
9589 if(
m2>= 220.0 &&
m2<1000.0)
9597 if(
m3>= 220.0 &&
m3<1000.0)
9605 if(
m2>= 225.0 &&
m2< 600.0)
9613 if(
m3>= 225.0 &&
m3< 600.0)
9621 if(
m2>= 220.0 &&
m2<1000.0)
9629 if(
m3>= 220.0 &&
m3<1000.0)
9637 if(
m2>= 220.0 &&
m2< 350.0)
9645 if(
m3>= 220.0 &&
m3< 350.0)
9651 if(
m2>= 220.0 &&
m2<2000.0)
9661 if(
m3>= 220.0 &&
m3<2000.0)
9669 if(
m2>= 220.0 &&
m2< 800.0)
9676 if(
m3>= 220.0 &&
m3< 800.0)
9682 if(
m2>= 800.0 &&
m2<2000.0)
9690 if(
m3>= 800.0 &&
m3<2000.0)
9699 if(
m2>= 220.0 &&
m2<2000.0)
9707 if(
m3>= 220.0 &&
m3<2000.0)
9716 if(
m2>= 220.0 &&
m2<400.0)
9724 if(
m3>= 220.0 &&
m3<400.0)
9738 if(
m2>= 220.0 &&
m2< 800.0)
9746 if(
m3>= 220.0 &&
m3< 800.0)
9754 if(
m2>= 800.0 &&
m2<2000.0)
9762 if(
m3>= 800.0 &&
m3<2000.0)
9770 if(
m3>= 175.0 && m3<1000.0 && m2 >=50.0 &&
m2 <910.0)
9778 if(
m2>= 175.0 && m2<1000.0 && m3 >=50.0 &&
m3 <910.0)
9786 if(
m3>= 50.0 && m3<1000.0 && m2 >=50.0 &&
m2 <1000.0)
9793 if(
m2>= 50.0 && m2<1000.0 && m3 >=50.0 &&
m3 <1000.0)
9800 if(
m3>= 50.0 && m3<1000.0 && m2 >=50.0 &&
m2 <1000.0)
9807 if(
m2>= 50.0 && m2<1000.0 && m3 >=50.0 &&
m3 <1000.0)
9814 if(
m2 >= 230.0 && m2 <800.0 && m3>=130.0 && m3<700.0 && m2-m3>=100)
9823 if(
m3 >= 230.0 && m3 <800.0 && m2>=130.0 && m2<700.0 && m3-m2>=100)
9830 if(
m2 >= 230.0 && m2 <800.0 && m3>=130.0 && m3<700.0 && m2-m3>=100)
9836 if(
m3 >= 230.0 && m3 <800.0 && m2>=130.0 && m2<700.0 && m3-m2>=100)
9845 if(
m2 >= 300.0 && m2 <800.0 && m3>=200.0 && m3<700.0 && m2-m3>=100)
9851 if(
m3 >= 300.0 && m3 <800.0 && m2>=200.0 && m2<700.0 && m3-m2>=100)
9860 if(
mHp>= 180.0 &&
mHp<1000.0)
9880 if(
mHp>= 180.0 &&
mHp< 600.0)
9895 if(
mHp>= 180.0 &&
mHp<2000.0)
9910 if(
mHp>= 180.0 &&
mHp<3000.0)
9924 if(
mHp>= 200.0 &&
mHp< 600.0)
9938 if(
mHp>= 180.0 &&
mHp< 600.0)
9951 if(
mHp>= 200.0 &&
mHp<2000.0)
9968 if(
mHp>= 200.0 &&
mHp<3000.0)
10003 if(
mH3 >= 15.0 &&
mH3 <= 62.0)
10017 if(
mH3 >= 12.0 &&
mH3 <= 60.0)
10031 if(
mH3 >= 15.0 &&
mH3 <= 62.0)
10038 if(
mH3 >= 4.2045 &&
mH3 <= 34.843)
10045 if(
mH3 >= 4.2 &&
mH3 <= 59.946)
10052 if(
mH3 >= 15.0 &&
mH3 <= 62.0)
10060 if(
mH3 >= 4.0 &&
mH3 <= 15.0)
10067 if(
mH3 >= 15.0 &&
mH3 <= 350.0)
10074 if(
mH3 >= 15.0 &&
mH3 <= 350.0)
10081 if(
mH3 >= 15.0 &&
mH3 <= 60.0)
10088 if(
mH3 >= 16.0 &&
mH3 <= 62.0)
10095 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))
10103 if(
mH3 >= 15.0 &&
mH3 <= 30.0)
10111 if(
mH3 >= 20.0 &&
mH3 <= 60.0)
10118 if(
mH3 >= 15.0 &&
mH3 <= 30.0)
10125 if(
mH3 >= 20.0 &&
mH3 <= 60.0)
10132 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))
10139 if(
mH3 >= 20.0 &&
mH3 <= 90.0)
10146 if(
mH3 >= 10.0 &&
mH3 <= 62.0)
10153 if(
mH3 >= 3.7 &&
mH3 <= 50.0)
10160 if(
mH3 >= 5.0 &&
mH3 <= 15.0)
10167 if(
mH3 >= 25.0 &&
mH3 <= 62.5)
10174 if(
mH3 >= 15.0 &&
mH3 <= 62.5)
10181 if(
mH3 >= 25.0 &&
mH3 <= 70.0)
10188 if(
mH3 >= 25.0 &&
mH3 <= 80.0)
10195 if(
mH3 >= 15.1 &&
mH3 <= 73.3)
10202 if(
mH3 >= 25.0 &&
mH3 <= 60.0)
10209 if(
mH3 >= 10.0 &&
mH3 <= 62.0)
10216 if(
mH3 >= 12.0 &&
mH3 <= 60.0)
10227 if(
mH2 >= 4.2045 &&
mH2 <= 34.843)
10234 if(
mH2 >= 4.2 &&
mH2 <= 59.946)
10241 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))
10249 if(
mH2 >= 15.0 &&
mH2 <= 30.0)
10257 if(
mH2 >= 20.0 &&
mH2 <= 60.0)
10264 if(
mH2 >= 15.0 &&
mH2 <= 30.0)
10271 if(
mH2 >= 20.0 &&
mH2 <= 60.0)
10278 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))
10285 if(
mH2 >= 80.0 &&
mH2 <= 109.9)
10292 if(
mH2 >= 70.0 &&
mH2 <= 110.0)
10299 if(
mH2 >= 66.0 &&
mH2 <= 110.0)
10306 if(
mH2 >= 20.0 &&
mH2 <= 116.0)
10313 if(
mH2 >= 12.0 &&
mH2 <= 120.0)
10320 if(
mH2 >= 4.0 &&
mH2 <= 120.0)
10371 if(
mHp >= 90.0 &&
mHp <= 160.0)
10378 if(
mHp >= 80.0 &&
mHp <= 160.0)
10385 if(
mHp >= 90.0 &&
mHp <= 150.0)
10392 if(
mH3 >= 15.0 && mH3 <= 75.0 && mHp >= (
mH3+85.0) &&
mHp <= 160.)
10399 if(
mHp >= 80.0 &&
mHp <= 160.0)
10406 if(
mHp >= 80.0 &&
mHp <= 160.0)
10413 if(
mHp >= 90.0 &&
mHp <= 150.0)
10420 if(
mH3 >= 15.1 && mH3 <= 70.2 && mHp >= 120.0 &&
mHp <= 160.0)
10427 if(
mHp >= 60.0 &&
mHp <= 168.0)
10433 double GeV2pb = 0.389379e9;
10437 if(
mHp >= 43.0 &&
mHp <= 95.0)
10444 if(
mHp >= 43.0 &&
mHp <= 95.0)
10451 if(
mHp >= 50.0 &&
mHp <= 93.0)
10458 if(
mHp >= 40.0 &&
mHp < 50.0)
10466 if(
mHp >= 40.0 && mHp <= 93.0 && mH3 >= 12.0 &&
mH3 <= (
mHp - 3.))
10475 if(
mHp >= 80.0 &&
mHp < 440.0)
10483 if(
mHp >= 90.0 &&
mHp < 500.0)
10491 if(
mHp >= 90.0 &&
mHp < 350.0)
10499 if(
mHp >= 100.0 &&
mHp < 800.0)
10507 if(
mHp >= 45.0 &&
mHp <= 102.0)
10518 if( RGEorder ==
"LO" ) flag=0;
10522 throw std::runtime_error(
"RGEorder can be only \"LO\" at the moment");
10527 double g1_at_MZ=sqrt(4.0*M_PI*
Ale/
cW2);
10528 double g2_at_MZ=sqrt(4.0*M_PI*
Ale/(1-
cW2));
10529 double g3_at_MZ=sqrt(4.0*M_PI*
Als);
10530 double v1_at_MZ=0.;
10531 double v2_at_MZ=0.;
10532 double etaU1_at_MZ=0.;
10533 double etaU2_at_MZ=0.;
10534 double etaD1_at_MZ=0.;
10535 double etaD2_at_MZ=0.;
10536 double etaL1_at_MZ=0.;
10537 double etaL2_at_MZ=0.;
10538 double m11sq_at_MZ=0.;
10539 double m12sq_at_MZ=0.;
10540 double m22sq_at_MZ=0.;
10541 double lambda1_at_MZ=
lambda1;
10542 double lambda2_at_MZ=
lambda2;
10543 double lambda3_at_MZ=
lambda3;
10544 double lambda4_at_MZ=
lambda4;
10548 double RpepsGTHDM=
myGTHDM->getRpepsGTHDM();
10549 double NLOuniscale=
myGTHDM->getNLOuniscaleGTHDM();
10553 double Q_A2HDM = log10(
MZ);
10555 if(fabs(Q_A2HDM-log10(
MZ))<0.005)
10583 double InitVals[21];
10584 InitVals[0]=g1_at_MZ;
10585 InitVals[1]=g2_at_MZ;
10586 InitVals[2]=g3_at_MZ;
10587 InitVals[3]=v1_at_MZ;
10588 InitVals[4]=v2_at_MZ;
10589 InitVals[5]=etaU1_at_MZ;
10590 InitVals[6]=etaU2_at_MZ;
10591 InitVals[7]=etaD1_at_MZ;
10592 InitVals[8]=etaD2_at_MZ;
10593 InitVals[9]=etaL1_at_MZ;
10594 InitVals[10]=etaL2_at_MZ;
10595 InitVals[11]=m11sq_at_MZ;
10596 InitVals[12]=m12sq_at_MZ;
10597 InitVals[13]=m22sq_at_MZ;
10598 InitVals[14]=lambda1_at_MZ;
10599 InitVals[15]=lambda2_at_MZ;
10600 InitVals[16]=lambda3_at_MZ;
10601 InitVals[17]=lambda4_at_MZ;
10602 InitVals[18]=lambda5_at_MZ;
10603 InitVals[19]=lambda6_at_MZ;
10604 InitVals[20]=lambda7_at_MZ;
10632 throw std::runtime_error(
"RGE are only defined for the CP conserving ATHDM at the moment");
10677 return std::numeric_limits<double>::quiet_NaN();
10896 if(
myGTHDM->getATHDMflag() ==
true)
10988 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