20 br_tautau(19981, 2, 0.),
22 br_mumu(19981, 2, 0.),
25 GammaHtot_SM(19981, 2, 0.),
26 log_cs_ggH_8(200, 2, 0.),
27 log_cs_VBF_8(200, 2, 0.),
28 log_cs_WH_8(200, 2, 0.),
29 log_cs_ZH_8(200, 2, 0.),
30 log_cs_ggH_13(200, 2, 0.),
31 log_cs_VBF_13(200, 2, 0.),
32 log_cs_WH_13(200, 2, 0.),
33 log_cs_ZH_13(200, 2, 0.),
34 log_cs_ttH_8(200, 2, 0.),
35 log_cs_ttH_13(200, 2, 0.),
36 log_cs_bbH_8(200, 2, 0.),
37 log_cs_bbH_13(200, 2, 0.),
38 log_cs_ggA_8(200, 2, 0.),
39 log_cs_ttA_8(200, 2, 0.),
40 log_cs_bbA_8(200, 2, 0.),
41 log_cs_ggA_13(200, 2, 0.),
42 log_cs_ttA_13(200, 2, 0.),
43 log_cs_bbA_13(200, 2, 0.),
44 log_cs_ggHp_8(744, 3, 0.),
45 log_cs_ggHp_13(1104, 3, 0.),
46 csrH_top_8(200, 2, 0.),
47 csrH_bottom_8(200, 2, 0.),
48 csrA_top_8(200, 2, 0.),
49 csrA_bottom_8(200, 2, 0.),
50 csrH_top_13(200, 2, 0.),
51 csrH_bottom_13(200, 2, 0.),
52 csrA_top_13(200, 2, 0.),
53 csrA_bottom_13(200, 2, 0.),
54 ATLAS8_pp_phi_gaga(108, 2, 0.),
55 ATLAS8_pp_phi_Zga_llga(141, 2, 0.),
56 ATLAS8_gg_phi_tautau(92, 2, 0.),
57 ATLAS8_bb_phi_tautau(92, 2, 0.),
58 ATLAS8_gg_A_hZ_tautauZ(79, 2, 0.),
59 ATLAS8_gg_A_hZ_bbZ(79, 2, 0.),
60 ATLAS8_gg_phi_tt(53, 2, 0.),
61 ATLAS8_gg_H_WW(13,2,0.),
62 ATLAS8_VBF_H_WW(13,2,0.),
63 ATLAS8_gg_H_ZZ(173,2,0.),
64 ATLAS8_VBF_H_ZZ(173,2,0.),
65 ATLAS8_gg_H_hh(75,2,0.),
66 ATLAS8_pp_phi_gaga_e(108, 2, 0.),
67 ATLAS8_pp_phi_Zga_llga_e(141, 2, 0.),
68 ATLAS8_gg_phi_tautau_e(92, 2, 0.),
69 ATLAS8_bb_phi_tautau_e(92, 2, 0.),
70 ATLAS8_gg_A_hZ_tautauZ_e(79, 2, 0.),
71 ATLAS8_gg_A_hZ_bbZ_e(79, 2, 0.),
72 ATLAS8_gg_phi_tt_e(53, 2, 0.),
73 ATLAS8_gg_H_WW_e(13,2,0.),
74 ATLAS8_VBF_H_WW_e(13,2,0.),
75 ATLAS8_gg_H_ZZ_e(173,2,0.),
76 ATLAS8_VBF_H_ZZ_e(173,2,0.),
77 ATLAS8_gg_H_hh_e(75,2,0.),
78 CMS8_mu_pp_H_VV(172, 2, 0.),
79 CMS8_mu_pp_H_VV_e(172, 2, 0.),
80 CMS8_gg_A_hZ_bbll(16, 2, 0.),
81 CMS8_pp_H_hh(71, 2, 0.),
82 CMS8_pp_H_hh_gagabb(85, 2, 0.),
83 CMS8_pp_H_hh_bbbb(167, 2, 0.),
84 CMS8_bb_phi_bb(81, 2, 0.),
85 CMS8_gg_phi_tautau(92,2,0.),
86 CMS8_bb_phi_tautau(92,2,0.),
87 CMS8_gg_phi_gaga(141,2,0.),
88 CMS8_pp_A_Zga_llga(101,2,0.),
89 CMS8_pp_phi_Zga(101,2,0.),
90 CMS8_gg_H_hh_bbtautau(10,2,0.),
91 CMS8_gg_A_hZ_tautaull(14,2,0.),
92 CMS8_pp_A_HZ_bbll(28718, 3, 0.),
93 CMS8_pp_H_AZ_bbll(29050, 3, 0.),
94 CMS8_pp_A_HZ_tautaull(400, 3, 0.),
95 CMS8_pp_H_AZ_tautaull(400, 3, 0.),
96 CMS8_gg_A_hZ_bbll_e(16, 2, 0.),
97 CMS8_pp_H_hh_e(71, 2, 0.),
98 CMS8_pp_H_hh_gagabb_e(85, 2, 0.),
99 CMS8_pp_H_hh_bbbb_e(167, 2, 0.),
100 CMS8_bb_phi_bb_e(81, 2, 0.),
101 CMS8_gg_phi_tautau_e(92,2,0.),
102 CMS8_bb_phi_tautau_e(92,2,0.),
103 CMS8_gg_phi_gaga_e(141,2,0.),
104 CMS8_pp_A_Zga_llga_e(101,2,0.),
105 CMS8_gg_H_hh_bbtautau_e(10,2,0.),
106 CMS8_gg_A_hZ_tautaull_e(14,2,0.),
111 ATLAS13_bb_phi_tt(61,2,0.),
112 ATLAS13_tt_phi_tt(61,2,0.),
113 ATLAS13_gg_phi_tautau(206,2,0.),
114 ATLAS13_bb_phi_tautau(206,2,0.),
115 ATLAS13_pp_phi_gaga(251,2,0.),
116 ATLAS13_pp_phi_Zga(216,2,0.),
117 ATLAS13_gg_phi_Zga_llga(216,2,0.),
118 ATLAS13_gg_H_ZZ_llllnunu(101,2,0.),
119 ATLAS13_VBF_H_ZZ_llllnunu(101,2,0.),
120 ATLAS13_gg_H_ZZ_llnunu(71,2,0.),
121 ATLAS13_gg_H_ZZ_llll(81,2,0.),
122 ATLAS13_VBF_H_ZZ_llll(81,2,0.),
123 ATLAS13_gg_H_ZZ_qqllnunu(271,2,0.),
124 ATLAS13_VBF_H_ZZ_qqllnunu(271,2,0.),
125 ATLAS13_gg_H_ZZ_llqq(271,2,0.),
126 ATLAS13_VBF_H_ZZ_llqq(271,2,0.),
127 ATLAS13_gg_H_ZZ_nunuqq(251,2,0.),
128 ATLAS13_gg_H_WW_enumumu(381,2,0.),
129 ATLAS13_VBF_H_WW_enumumu(281,2,0.),
130 ATLAS13_gg_H_WW_lnuqq(271,2,0.),
131 ATLAS13_VBF_H_WW_lnuqq(271,2,0.),
132 ATLAS13_pp_H_VV_qqqq(181,2,0.),
133 ATLAS13_pp_H_hh_bbbb(271,2,0.),
134 ATLAS13_pp_H_hh_gagabb(26,2,0.),
135 ATLAS13_pp_H_hh_gagaWW(25,2,0.),
136 ATLAS13_gg_A_Zh_Zbb(181,2,0.),
137 ATLAS13_bb_A_Zh_Zbb(181,2,0.),
138 ATLAS13_bb_phi_tt_e(61,2,0.),
139 ATLAS13_tt_phi_tt_e(61,2,0.),
140 ATLAS13_gg_phi_tautau_e(206,2,0.),
141 ATLAS13_bb_phi_tautau_e(206,2,0.),
142 ATLAS13_pp_phi_gaga_e(251,2,0.),
143 ATLAS13_pp_phi_Zga_e(216,2,0.),
144 ATLAS13_gg_phi_Zga_llga_e(216,2,0.),
145 ATLAS13_gg_H_ZZ_llllnunu_e(101,2,0.),
146 ATLAS13_VBF_H_ZZ_llllnunu_e(101,2,0.),
147 ATLAS13_gg_H_ZZ_llnunu_e(71,2,0.),
148 ATLAS13_gg_H_ZZ_llll_e(81,2,0.),
149 ATLAS13_VBF_H_ZZ_llll_e(81,2,0.),
150 ATLAS13_gg_H_ZZ_qqllnunu_e(271,2,0.),
151 ATLAS13_VBF_H_ZZ_qqllnunu_e(271,2,0.),
152 ATLAS13_gg_H_ZZ_llqq_e(271,2,0.),
153 ATLAS13_VBF_H_ZZ_llqq_e(271,2,0.),
154 ATLAS13_gg_H_ZZ_nunuqq_e(251,2,0.),
155 ATLAS13_gg_H_WW_enumumu_e(381,2,0.),
156 ATLAS13_VBF_H_WW_enumumu_e(281,2,0.),
157 ATLAS13_gg_H_WW_lnuqq_e(271,2,0.),
158 ATLAS13_VBF_H_WW_lnuqq_e(271,2,0.),
159 ATLAS13_pp_H_VV_qqqq_e(181,2,0.),
160 ATLAS13_pp_H_hh_bbbb_e(271,2,0.),
161 ATLAS13_pp_H_hh_gagabb_e(26,2,0.),
162 ATLAS13_pp_H_hh_gagaWW_e(25,2,0.),
163 ATLAS13_gg_A_Zh_Zbb_e(181,2,0.),
164 ATLAS13_bb_A_Zh_Zbb_e(181,2,0.),
165 CMS13_pp_phi_bb(66,2,0.),
166 CMS13_gg_phi_tautau(312,2,0.),
167 CMS13_bb_phi_tautau(312,2,0.),
168 CMS13_gg_phi_gaga(351,2,0.),
169 CMS13_pp_phi_Zga_llga(171,2,0.),
170 CMS13_pp_phi_Zga_qqga(236,2,0.),
171 CMS13_ggF_phi_Zga(366,2,0.),
172 CMS13_pp_H_ZZ_llnunu(191,2,0.),
173 CMS13_gg_H_ZZ_llnunu(131,2,0.),
174 CMS13_VBF_H_ZZ_llnunu(131,2,0.),
175 CMS13_pp_H_ZZ_llll(241,2,0.),
176 CMS13_VBFVH_H_ZZ_llll(241,2,0.),
177 CMS13_pp_H_ZZ_llqq(151,2,0.),
178 CMS13_ggFVBF_H_WW_lnulnu(81,2,0.),
179 CMS13_pp_H_hh_bbbb(95,2,0.),
180 CMS13_ggF_H_hh_bbbb(226,2,0.),
181 CMS13_pp_H_hh_gagabb(66,2,0.),
182 CMS13_pp_H_hh_bbtautau(66,2,0.),
183 CMS13_pp_H_hh_bbtautau1(66,2,0.),
184 CMS13_pp_H_hh_bblnulnu(65,2,0.),
185 CMS13_pp_H_hh_bbVV(65,2,0.),
186 CMS13_pp_phi_bb_e(66,2,0.),
187 CMS13_gg_phi_tautau_e(312,2,0.),
188 CMS13_bb_phi_tautau_e(312,2,0.),
189 CMS13_gg_phi_gaga_e(351,2,0.),
190 CMS13_pp_phi_Zga_llga_e(171,2,0.),
191 CMS13_pp_phi_Zga_qqga_e(236,2,0.),
192 CMS13_ggF_phi_Zga_e(366,2,0.),
193 CMS13_pp_H_ZZ_llnunu_e(191,2,0.),
194 CMS13_gg_H_ZZ_llnunu_e(131,2,0.),
195 CMS13_VBF_H_ZZ_llnunu_e(131,2,0.),
196 CMS13_pp_H_ZZ_llll_e(241,2,0.),
197 CMS13_VBFVH_H_ZZ_llll_e(241,2,0.),
198 CMS13_pp_H_ZZ_llqq_e(151,2,0.),
199 CMS13_ggFVBF_H_WW_lnulnu_e(81,2,0.),
200 CMS13_pp_H_hh_bbbb_e(95,2,0.),
201 CMS13_ggF_H_hh_bbbb_e(226,2,0.),
202 CMS13_pp_H_hh_gagabb_e(66,2,0.),
203 CMS13_pp_H_hh_bbtautau_e(66,2,0.),
204 CMS13_pp_H_hh_bbtautau1_e(66,2,0.),
205 CMS13_pp_H_hh_bblnulnu_e(65,2,0.),
206 CMS13_pp_H_hh_bbVV_e(65,2,0.),
207 temp1(1,2,0.), temp2(1,2,0.), temp3(1,2,0.), temp4(1,2,0.), temp5(1,2,0.), temp6(1,2,0.), temp7(1,2,0.), temp8(1,2,0.), temp9(1,2,0.), temp10(1,2,0.),
208 temp11(1,2,0.), temp12(1,2,0.), temp13(1,2,0.), temp14(1,2,0.), temp15(1,2,0.), temp16(1,2,0.), temp17(1,2,0.), temp18(1,2,0.), temp19(1,2,0.), temp20(1,2,0.),
209 temp21(1,2,0.), temp22(1,2,0.), temp23(1,2,0.), temp24(1,2,0.), temp25(1,2,0.), temp26(1,2,0.), temp27(1,2,0.), temp28(1,2,0.), temp29(1,2,0.), temp30(1,2,0.),
210 temp31(1,2,0.), temp32(1,2,0.), temp33(1,2,0.), temp34(1,2,0.), temp35(1,2,0.), temp36(1,2,0.), temp37(1,2,0.), temp38(1,2,0.), temp39(1,2,0.), temp40(1,2,0.),
211 temp1e(1,2,0.), temp2e(1,2,0.), temp3e(1,2,0.), temp4e(1,2,0.), temp5e(1,2,0.), temp6e(1,2,0.), temp7e(1,2,0.), temp8e(1,2,0.), temp9e(1,2,0.), temp10e(1,2,0.),
212 temp11e(1,2,0.), temp12e(1,2,0.), temp13e(1,2,0.), temp14e(1,2,0.), temp15e(1,2,0.), temp16e(1,2,0.), temp17e(1,2,0.), temp18e(1,2,0.), temp19e(1,2,0.), temp20e(1,2,0.),
213 temp21e(1,2,0.), temp22e(1,2,0.), temp23e(1,2,0.), temp24e(1,2,0.), temp25e(1,2,0.), temp26e(1,2,0.), temp27e(1,2,0.), temp28e(1,2,0.), temp29e(1,2,0.), temp30e(1,2,0.),
214 temp31e(1,2,0.), temp32e(1,2,0.), temp33e(1,2,0.), temp34e(1,2,0.), temp35e(1,2,0.), temp36e(1,2,0.), temp37e(1,2,0.), temp38e(1,2,0.), temp39e(1,2,0.), temp40e(1,2,0.),
215 ATLAS8_pp_Hpm_taunu(83,2,0.),
216 ATLAS8_pp_Hpm_tb(41,2,0.),
217 ATLAS8_pp_Hpm_taunu_e(83,2,0.),
218 ATLAS8_pp_Hpm_tb_e(41,2,0.),
219 CMS8_pp_Hp_taunu(43,2,0.),
220 CMS8_pp_Hp_tb(43,2,0.),
221 CMS8_pp_Hp_taunu_e(43,2,0.),
222 CMS8_pp_Hp_tb_e(43,2,0.),
223 ATLAS13_pp_Hpm_taunu(181,2,0.),
224 ATLAS13_pp_Hp_tb1(71,2,0.),
225 ATLAS13_pp_Hp_tb2(181,2,0.),
226 ATLAS13_pp_Hpm_taunu_e(181,2,0.),
227 ATLAS13_pp_Hp_tb1_e(71,2,0.),
228 ATLAS13_pp_Hp_tb2_e(181,2,0.),
229 CMS13_pp_Hpm_taunu(283,2,0.),
230 CMS13_pp_Hpm_taunu_e(283,2,0.),
231 arraybsgamma(1111, 3, 0.),
232 myTHDM(static_cast<const
THDM*> (&SM_i)),
263 const int NumPar,
const double params[])
const {
267 for(
int j=0; j<NumPar; j++)
268 bCache &= (params[j] == cache[j][i].real());
269 if (bCache)
return i;
275 const int NumPar,
const double params[])
const {
279 for(
int j=0; j<NumPar; j++)
280 bCache &= (params[j] == cache[j][i]);
281 if (bCache)
return i;
287 const double params[],
const gslpp::complex newResult)
const {
290 for(
int j=0; j<NumPar+1; j++)
291 cache[j][i] = cache[j][i-1];
294 for(
int j=0; j<NumPar; j++) {
295 cache[j][0] = gslpp::complex(params[j], 0.0,
false);
296 cache[NumPar][0] = newResult;
301 const double params[],
const double newResult)
const {
304 for(
int j=0; j<NumPar+1; j++)
305 cache[j][i] = cache[j][i-1];
308 for(
int j=0; j<NumPar; j++) {
309 cache[j][0] = params[j];
310 cache[NumPar][0] = newResult;
320 double params[] = {MZ2, MW2,
mHh2};
326 gslpp::complex newResult =
PV.
B0(MZ2, 0., MW2,
mHh2);
334 double params[] = {MZ2, MW2,
mHl2};
340 gslpp::complex newResult =
PV.
B0(MZ2, 0., MW2,
mHl2);
348 double params[] = {MZ2,
mHh2};
354 gslpp::complex newResult =
PV.
B0(MZ2, 0., MZ2,
mHh2);
362 double params[] = {MZ2,
mHl2};
368 gslpp::complex newResult =
PV.
B0(MZ2, 0., MZ2,
mHl2);
376 double params[] = {MZ2, MW2,
mHh2};
382 gslpp::complex newResult =
PV.
B0(MZ2, MW2, MW2,
mHh2);
390 double params[] = {MZ2, MW2,
mHl2};
396 gslpp::complex newResult =
PV.
B0(MZ2, MW2, MW2,
mHl2);
404 double params[] = {MZ2,
mHh2};
410 gslpp::complex newResult =
PV.
B0(MZ2, MZ2, MZ2,
mHh2);
418 double params[] = {MZ2,
mHl2};
424 gslpp::complex newResult =
PV.
B0(MZ2, MZ2, MZ2,
mHl2);
432 double params[] = {MZ2,
mHl2};
438 gslpp::complex newResult =
PV.
B0(MZ2, 0.0, 0.0,
mHl2);
446 double params[] = {MZ2,
mHh2};
452 gslpp::complex newResult =
PV.
B0(MZ2, 0.0, 0.0,
mHh2);
460 double params[] = {MZ2,
mHp2,
mHl2};
474 double params[] = {MZ2,
mHp2,
mHh2};
488 double params[] = {MZ2,
mA2,
mHl2};
494 gslpp::complex newResult =
PV.
B0(MZ2, 0.0,
mA2,
mHl2);
502 double params[] = {MZ2,
mA2,
mHh2};
508 gslpp::complex newResult =
PV.
B0(MZ2, 0.0,
mA2,
mHh2);
516 double params[] = {MZ2,
mHl2};
522 gslpp::complex newResult =
PV.
B0(MZ2,
mHl2, 0.0, 0.0);
530 double params[] = {MZ2,
mHl2,
mHp2};
544 double params[] = {MZ2,
mHl2,
mA2};
550 gslpp::complex newResult =
PV.
B0(MZ2,
mHl2, 0.0,
mA2);
558 double params[] = {MZ2,
mHl2};
572 double params[] = {MZ2,
mHl2,
mHh2};
586 double params[] = {MZ2,
mHl2,
mHh2};
600 double params[] = {MZ2,
mHl2,
mHp2};
614 double params[] = {MZ2,
mHl2,
mA2};
628 double params[] = {MZ2,
mHh2};
634 gslpp::complex newResult =
PV.
B0(MZ2,
mHh2, 0.0, 0.0);
642 double params[] = {MZ2,
mHh2,
mHp2};
656 double params[] = {MZ2,
mHh2,
mA2};
662 gslpp::complex newResult =
PV.
B0(MZ2,
mHh2, 0.0,
mA2);
670 double params[] = {MZ2,
mHh2,
mHl2};
684 double params[] = {MZ2,
mHh2,
mHl2};
698 double params[] = {MZ2,
mHh2};
712 double params[] = {MZ2,
mHh2,
mHp2};
726 double params[] = {MZ2,
mHh2,
mA2};
740 double params[] = {MZ2,
mHp2,
mHl2};
754 double params[] = {MZ2,
mHp2,
mHh2};
768 double params[] = {MZ2,
mHp2,
mHl2};
782 double params[] = {MZ2,
mHp2,
mHh2};
796 double params[] = {MZ2,
mA2,
mHl2};
802 gslpp::complex newResult =
PV.
B0(MZ2,
mA2, 0.0,
mHl2);
810 double params[] = {MZ2,
mA2,
mHh2};
816 gslpp::complex newResult =
PV.
B0(MZ2,
mA2, 0.0,
mHh2);
824 double params[] = {MZ2,
mA2,
mHl2};
838 double params[] = {MZ2,
mA2,
mHh2};
854 double params[] = {MZ2,
mHl2};
860 gslpp::complex newResult =
PV.
B0p(MZ2, 0.0, 0.0,
mHl2);
868 double params[] = {MZ2,
mHh2};
874 gslpp::complex newResult =
PV.
B0p(MZ2, 0.0, 0.0,
mHh2);
882 double params[] = {MZ2,
mHp2,
mHl2};
896 double params[] = {MZ2,
mHp2,
mHh2};
910 double params[] = {MZ2,
mHp2,
mA2};
924 double params[] = {MZ2,
mA2,
mHl2};
938 double params[] = {MZ2,
mA2,
mHh2};
952 double params[] = {MZ2,
mHl2};
958 gslpp::complex newResult =
PV.
B0p(MZ2,
mHl2, 0.0, 0.0);
966 double params[] = {MZ2,
mHl2,
mHp2};
980 double params[] = {MZ2,
mHl2,
mA2};
994 double params[] = {MZ2,
mHl2};
1008 double params[] = {MZ2,
mHl2,
mHh2};
1022 double params[] = {MZ2,
mHl2,
mHh2};
1036 double params[] = {MZ2,
mHl2,
mHp2};
1050 double params[] = {MZ2,
mHl2,
mA2};
1064 double params[] = {MZ2,
mHh2};
1070 gslpp::complex newResult =
PV.
B0p(MZ2,
mHh2, 0.0, 0.0);
1078 double params[] = {MZ2,
mHh2,
mHp2};
1092 double params[] = {MZ2,
mHh2,
mA2};
1098 gslpp::complex newResult =
PV.
B0p(MZ2,
mHh2, 0.0,
mA2);
1106 double params[] = {MZ2,
mHh2,
mHl2};
1120 double params[] = {MZ2,
mHh2,
mHl2};
1134 double params[] = {MZ2,
mHh2};
1148 double params[] = {MZ2,
mHh2,
mHp2};
1162 double params[] = {MZ2,
mHh2,
mA2};
1176 double params[] = {MZ2,
mHp2,
mHl2};
1190 double params[] = {MZ2,
mHp2,
mHh2};
1204 double params[] = {MZ2,
mHp2,
mA2};
1210 gslpp::complex newResult =
PV.
B0p(MZ2,
mHp2, 0.0,
mA2);
1218 double params[] = {MZ2,
mHp2,
mHl2};
1232 double params[] = {MZ2,
mHp2,
mHh2};
1246 double params[] = {MZ2,
mA2,
mHl2};
1252 gslpp::complex newResult =
PV.
B0p(MZ2,
mA2, 0.0,
mHl2);
1260 double params[] = {MZ2,
mA2,
mHh2};
1266 gslpp::complex newResult =
PV.
B0p(MZ2,
mA2, 0.0,
mHh2);
1274 double params[] = {MZ2,
mA2,
mHp2};
1280 gslpp::complex newResult =
PV.
B0p(MZ2,
mA2, 0.0,
mHp2);
1288 double params[] = {MZ2,
mA2,
mHl2};
1302 double params[] = {MZ2,
mA2,
mHh2};
1318 double params[] = {MZ2,
mA2,
mHp2};
1332 double params[] = {MZ2,
mHh2,
mA2};
1346 double params[] = {MZ2,
mHh2,
mHp2};
1360 double params[] = {MZ2,
mHl2,
mA2};
1374 double params[] = {MZ2,
mHl2,
mHp2};
1388 double params[] = {MZ2,
mHp2};
1402 double params[] = {MZ2, MW2,
mHh2};
1408 gslpp::complex newResult =
PV.
B00(MZ2, MW2, MW2,
mHh2);
1416 double params[] = {MZ2, MW2,
mHl2};
1422 gslpp::complex newResult =
PV.
B00(MZ2, 0., MW2,
mHl2);
1430 double params[] = {MZ2,
mHh2};
1436 gslpp::complex newResult =
PV.
B00(MZ2, 0., MZ2,
mHh2);
1444 double params[] = {MZ2,
mHl2};
1450 gslpp::complex newResult =
PV.
B00(MZ2, 0., MZ2,
mHl2);
1458 double params[] = {MZ2, MW2,
mA2,
mHp2};
1464 gslpp::complex newResult =
PV.
B00(MZ2, MW2,
mA2,
mHp2);
1472 double params[] = {MZ2, MW2,
mHh2,
mHp2};
1486 double params[] = {MZ2, MW2,
mHl2,
mHp2};
1500 double params[] = {MZ2, MW2,
mHp2};
1514 double params[] = {MZ2, MW2,
mHh2};
1520 gslpp::complex newResult =
PV.
B00(MZ2, MW2, MW2,
mHh2);
1528 double params[] = {MZ2, MW2,
mHl2};
1534 gslpp::complex newResult =
PV.
B00(MZ2, MW2, MW2,
mHl2);
1542 double params[] = {MZ2,
mHh2,
mA2};
1548 gslpp::complex newResult =
PV.
B00(MZ2, MZ2,
mHh2,
mA2);
1556 double params[] = {MZ2,
mHl2,
mA2};
1562 gslpp::complex newResult =
PV.
B00(MZ2, MZ2,
mHl2,
mA2);
1570 double params[] = {MZ2,
mHp2};
1584 double params[] = {MZ2,
mHh2};
1590 gslpp::complex newResult =
PV.
B00(MZ2, MZ2, MZ2,
mHh2);
1598 double params[] = {MZ2,
mHl2};
1604 gslpp::complex newResult =
PV.
B00(MZ2, MZ2, MZ2,
mHl2);
1612 std::stringstream br1,br2,br3,br4,br5,br6,br7;
1613 std::stringstream dw1;
1614 std::stringstream cs1,cs2,cs3,cs4,cs5,cs6,cs7,cs8,cs9;
1615 std::stringstream cs11,cs12,cs13,cs14,cs15,cs16,cs17,cs18,cs19;
1616 std::stringstream cs20,cs21;
1617 std::stringstream csr1,csr2,csr3,csr4;
1618 std::stringstream csr11,csr12,csr13,csr14;
1619 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;
1620 std::stringstream ex1e,ex2e,ex3e,ex4e,ex5e,ex6e,ex7e,ex8e,ex9e,ex10e,ex11e,ex12e,ex13e,ex14e,ex15e,ex16e,ex17e,ex18e,ex19e,ex20e,ex21e,ex22e,ex23e;
1622 std::stringstream ex24,ex25,ex26,ex27,ex28,ex29,ex30,ex31,ex32,ex33,ex34,ex35,ex36,ex37,ex38,ex39,ex40,ex41,ex42,ex43,ex44,\
1623 ex45,ex46,ex47,ex48,ex49,ex50,ex51,ex52,ex53,ex54,ex55,ex56;
1624 std::stringstream ex24e,ex25e,ex26e,ex27e,ex28e,ex29e,ex30e,ex31e,ex32e,ex33e,ex34e,ex35e,ex36e,ex37e,ex38e,ex39e,ex40e,ex41e,ex42e,ex43e,ex44e,\
1625 ex45e,ex46e,ex47e,ex48e,ex49e,ex50e,ex51e,ex52e,ex53e,ex54e,ex55e,ex56e;
1626 std::stringstream ex57,ex58,ex59,ex60,ex61,ex62,ex63,ex64,ex65,ex66,ex67,ex68,ex69,ex70,ex71,ex72,ex73,ex74,ex75,ex76,ex77,\
1627 ex78,ex79,ex80,ex81,ex82,ex83,ex84,ex85,ex86,ex87,ex88,ex89,ex90,ex91,ex92,ex93,ex94,ex95,ex96,ex97,ex98;
1628 std::stringstream ex57e,ex58e,ex59e,ex60e,ex61e,ex62e,ex63e,ex64e,ex65e,ex66e,ex67e,ex68e,ex69e,ex70e,ex71e,ex72e,ex73e,ex74e,ex75e,ex76e,ex77e,\
1629 ex78e,ex79e,ex80e,ex81e,ex82e,ex83e,ex84e,ex85e,ex86e,ex87e,ex88e,ex89e,ex90e,ex91e,ex92e,ex93e,ex94e,ex95e,ex96e,ex97e,ex98e;
1630 std::stringstream ex99,ex100,ex101,ex102,ex103,ex104,ex105,ex106,ex107,ex108,ex109,ex110,ex111,ex112,ex113,ex114,ex115,ex116,ex117,ex118,ex119,ex120;
1631 std::stringstream ex99e,ex100e,ex101e,ex102e,ex103e,ex104e,ex105e,ex106e,ex107e,ex108e,ex109e,ex110e,ex111e,ex112e,ex113e,ex114e,ex115e,ex116e,ex117e,ex118e,ex119e,ex120e;
1632 std::stringstream ex121,ex122,ex123,ex124;
1633 std::stringstream bsg1;
1635 std::cout<<
"reading tables"<<std::endl;
1638 std::stringstream path;
1639 path << getenv(
"HEPFITTABS") <<
"/THDM/tabs/";
1640 std::string tablepath=path.str();
1642 br1 << tablepath <<
"br1.dat";
1657 br2 << tablepath <<
"br2.dat";
1659 br3 << tablepath <<
"br3.dat";
1661 br4 << tablepath <<
"br4.dat";
1663 br5 << tablepath <<
"br5.dat";
1665 br6 << tablepath <<
"br6.dat";
1667 br7 << tablepath <<
"br7.dat";
1669 dw1 << tablepath <<
"dw1.dat";
1671 cs1 << tablepath <<
"log_cs_ggH_8.dat";
1673 cs11 << tablepath <<
"log_cs_ggH_13.dat";
1675 cs2 << tablepath <<
"log_cs_VBF_8.dat";
1677 cs12 << tablepath <<
"log_cs_VBF_13.dat";
1679 cs3 << tablepath <<
"log_cs_WH_8.dat";
1681 cs13 << tablepath <<
"log_cs_WH_13.dat";
1683 cs4 << tablepath <<
"log_cs_ZH_8.dat";
1685 cs14 << tablepath <<
"log_cs_ZH_13.dat";
1687 cs5 << tablepath <<
"log_cs_ttH_8.dat";
1689 cs15 << tablepath <<
"log_cs_ttH_13.dat";
1691 cs6 << tablepath <<
"log_cs_bbH_8.dat";
1693 cs16 << tablepath <<
"log_cs_bbH_13.dat";
1695 cs7 << tablepath <<
"log_cs_ggA_8.dat";
1697 cs17 << tablepath <<
"log_cs_ggA_13.dat";
1699 cs8 << tablepath <<
"log_cs_ttA_8.dat";
1701 cs18 << tablepath <<
"log_cs_ttA_13.dat";
1703 cs9 << tablepath <<
"log_cs_bbA_8.dat";
1705 cs19 << tablepath <<
"log_cs_bbA_13.dat";
1707 cs20 << tablepath <<
"log_cs_ggHp_8.dat";
1709 cs21 << tablepath <<
"log_cs_ggHp_13.dat";
1711 csr1 << tablepath <<
"csrH_top_8.dat";
1713 csr11 << tablepath <<
"csrH_top_13.dat";
1715 csr2 << tablepath <<
"csrH_bottom_8.dat";
1717 csr12 << tablepath <<
"csrH_bottom_13.dat";
1719 csr3 << tablepath <<
"csrA_top_8.dat";
1721 csr13 << tablepath <<
"csrA_top_13.dat";
1723 csr4 << tablepath <<
"csrA_bottom_8.dat";
1725 csr14 << tablepath <<
"csrA_bottom_13.dat";
1727 ex1 << tablepath <<
"150400936.dat";
1729 ex1e << tablepath <<
"150400936_e.dat";
1731 ex2 << tablepath <<
"150404710.dat";
1733 ex2e << tablepath <<
"150404710_e.dat";
1735 ex3 << tablepath <<
"160306896.dat";
1737 ex3e << tablepath <<
"160306896_e.dat";
1739 ex4 << tablepath <<
"150304114.dat";
1741 ex4e << tablepath <<
"150304114_e.dat";
1743 ex5 << tablepath <<
"14076583.dat";
1745 ex5e << tablepath <<
"14076583_e.dat";
1747 ex6 << tablepath <<
"14096064_a.dat";
1749 ex6e << tablepath <<
"14096064_a_e.dat";
1751 ex7 << tablepath <<
"14096064_b.dat";
1753 ex7e << tablepath <<
"14096064_b_e.dat";
1755 ex8 << tablepath <<
"150204478_a.dat";
1757 ex8e << tablepath <<
"150204478_a_e.dat";
1759 ex9 << tablepath <<
"150204478_b.dat";
1761 ex9e << tablepath <<
"150204478_b_e.dat";
1763 ex10 << tablepath <<
"150608329.dat";
1765 ex10e << tablepath <<
"150608329_e.dat";
1767 ex11 << tablepath <<
"150507018.dat";
1769 ex11e << tablepath <<
"150507018_e.dat";
1771 ex12 << tablepath <<
"CMS-PAS-HIG-14-029_a.dat";
1773 ex12e << tablepath <<
"CMS-PAS-HIG-14-029_a_e.dat";
1775 ex13 << tablepath <<
"CMS-PAS-HIG-14-029_b.dat";
1777 ex13e << tablepath <<
"CMS-PAS-HIG-14-029_b_e.dat";
1779 ex14 << tablepath <<
"150602301.dat";
1781 ex14e << tablepath <<
"150602301_e.dat";
1795 ex15 << tablepath <<
"150900389_a.dat";
1797 ex15e << tablepath <<
"150900389_a_e.dat";
1799 ex16 << tablepath <<
"150900389_b.dat";
1801 ex16e << tablepath <<
"150900389_b_e.dat";
1803 ex17 << tablepath <<
"150904670.dat";
1805 ex17e << tablepath <<
"150904670_e.dat";
1807 ex18 << tablepath <<
"151001181_a.dat";
1809 ex18e << tablepath <<
"151001181_a_e.dat";
1811 ex19 << tablepath <<
"151001181_b.dat";
1813 ex19e << tablepath <<
"151001181_b_e.dat";
1815 ex20 << tablepath <<
"150705930_a.dat";
1817 ex20e << tablepath <<
"150705930_a_e.dat";
1819 ex21 << tablepath <<
"150705930_b.dat";
1821 ex21e << tablepath <<
"150705930_b_e.dat";
1823 ex22 << tablepath <<
"CMS-PAS-HIG-15-013.dat";
1825 ex22e << tablepath <<
"CMS-PAS-HIG-15-013_e.dat";
1827 ex23 << tablepath <<
"CMS-PAS-HIG-16-014.dat";
1829 ex23e << tablepath <<
"CMS-PAS-HIG-16-014_e.dat";
1832 ex24 << tablepath <<
"ATLAS-CONF-2016-104_b.dat";
1834 ex24e << tablepath <<
"ATLAS-CONF-2016-104_b_e.dat";
1836 ex25 << tablepath <<
"ATLAS-CONF-2016-104_a.dat";
1838 ex25e << tablepath <<
"ATLAS-CONF-2016-104_a_e.dat";
1840 ex26 << tablepath <<
"170907242_a.dat";
1842 ex26e << tablepath <<
"170907242_a_e.dat";
1844 ex27 << tablepath <<
"170907242_b.dat";
1846 ex27e << tablepath <<
"170907242_b_e.dat";
1848 ex28 << tablepath <<
"170704147.dat";
1850 ex28e << tablepath <<
"170704147_e.dat";
1852 ex29 << tablepath <<
"ATLAS-CONF-2016-044.dat";
1854 ex29e << tablepath <<
"ATLAS-CONF-2016-044_e.dat";
1856 ex30 << tablepath <<
"ATLAS-CONF-2016-056.dat";
1858 ex30e << tablepath <<
"ATLAS-CONF-2016-056_e.dat";
1860 ex31 << tablepath <<
"ATLAS-CONF-2016-079_a.dat";
1862 ex31e << tablepath <<
"ATLAS-CONF-2016-079_a_e.dat";
1864 ex32 << tablepath <<
"ATLAS-CONF-2016-079_b.dat";
1866 ex32e << tablepath <<
"ATLAS-CONF-2016-079_b_e.dat";
1868 ex33 << tablepath <<
"ATLAS-CONF-2016-082_a.dat";
1870 ex33e << tablepath <<
"ATLAS-CONF-2016-082_a_e.dat";
1872 ex34 << tablepath <<
"ATLAS-CONF-2016-082_b.dat";
1874 ex34e << tablepath <<
"ATLAS-CONF-2016-082_b_e.dat";
1876 ex35 << tablepath <<
"ATLAS-CONF-2016-082_c.dat";
1878 ex35e << tablepath <<
"ATLAS-CONF-2016-082_c_e.dat";
1880 ex36 << tablepath <<
"171001123_a.dat";
1882 ex36e << tablepath <<
"171001123_a_e.dat";
1884 ex37 << tablepath <<
"171001123_b.dat";
1886 ex37e << tablepath <<
"171001123_b_e.dat";
1888 ex38 << tablepath <<
"171007235_a.dat";
1890 ex38e << tablepath <<
"171007235_a_e.dat";
1892 ex39 << tablepath <<
"ATLAS-CONF-2016-049.dat";
1894 ex39e << tablepath <<
"ATLAS-CONF-2016-049_e.dat";
1896 ex40 << tablepath <<
"ATLAS-CONF-2016-004.dat";
1898 ex40e << tablepath <<
"ATLAS-CONF-2016-004_e.dat";
1900 ex41 << tablepath <<
"ATLAS-CONF-2016-071.dat";
1902 ex41e << tablepath <<
"ATLAS-CONF-2016-071_e.dat";
1904 ex42 << tablepath <<
"171206518_a.dat";
1906 ex42e << tablepath <<
"171206518_a_e.dat";
1908 ex43 << tablepath <<
"171206518_b.dat";
1910 ex43e << tablepath <<
"171206518_b_e.dat";
1912 ex44 << tablepath <<
"CMS-PAS-HIG-16-025.dat";
1914 ex44e << tablepath <<
"CMS-PAS-HIG-16-025_e.dat";
1916 ex45 << tablepath <<
"CMS-PAS-HIG-16-037_a.dat";
1918 ex45e << tablepath <<
"CMS-PAS-HIG-16-037_a_e.dat";
1920 ex46 << tablepath <<
"CMS-PAS-HIG-16-037_b.dat";
1922 ex46e << tablepath <<
"CMS-PAS-HIG-16-037_b_e.dat";
1924 ex47 << tablepath <<
"160902507.dat";
1926 ex47e << tablepath <<
"160902507_e.dat";
1928 ex48 << tablepath <<
"CMS-PAS-EXO-16-034.dat";
1930 ex48e << tablepath <<
"CMS-PAS-EXO-16-034_e.dat";
1932 ex49 << tablepath <<
"CMS-PAS-EXO-16-035.dat";
1934 ex49e << tablepath <<
"CMS-PAS-EXO-16-035_e.dat";
1936 ex50 << tablepath <<
"CMS-PAS-HIG-16-033_a.dat";
1938 ex50e << tablepath <<
"CMS-PAS-HIG-16-033_a_e.dat";
1940 ex51 << tablepath <<
"CMS-PAS-HIG-16-033_b.dat";
1942 ex51e << tablepath <<
"CMS-PAS-HIG-16-033_b_e.dat";
1944 ex52 << tablepath <<
"CMS-PAS-HIG-16-023.dat";
1946 ex52e << tablepath <<
"CMS-PAS-HIG-16-023_e.dat";
1948 ex53 << tablepath <<
"180603548.dat";
1950 ex53e << tablepath <<
"180603548_e.dat";
1952 ex54 << tablepath <<
"180600408.dat";
1954 ex54e << tablepath <<
"180600408_e.dat";
1956 ex55 << tablepath <<
"CMS-PAS-HIG-16-029.dat";
1958 ex55e << tablepath <<
"CMS-PAS-HIG-16-029_e.dat";
1960 ex56 << tablepath <<
"CMS-PAS-HIG-16-011.dat";
1962 ex56e << tablepath <<
"CMS-PAS-HIG-16-011_e.dat";
1965 ex57 << tablepath <<
"t1.dat";
1967 ex57e << tablepath <<
"t1_e.dat";
1969 ex58 << tablepath <<
"t2.dat";
1971 ex58e << tablepath <<
"t2_e.dat";
1973 ex59 << tablepath <<
"t3.dat";
1975 ex59e << tablepath <<
"t3_e.dat";
1977 ex60 << tablepath <<
"t4.dat";
1979 ex60e << tablepath <<
"t4_e.dat";
1981 ex61 << tablepath <<
"t5.dat";
1983 ex61e << tablepath <<
"t5_e.dat";
1985 ex62 << tablepath <<
"t6.dat";
1987 ex62e << tablepath <<
"t6_e.dat";
1989 ex63 << tablepath <<
"t7.dat";
1991 ex63e << tablepath <<
"t7_e.dat";
1993 ex64 << tablepath <<
"t8.dat";
1995 ex64e << tablepath <<
"t8_e.dat";
1997 ex65 << tablepath <<
"t9.dat";
1999 ex65e << tablepath <<
"t9_e.dat";
2001 ex66 << tablepath <<
"t10.dat";
2003 ex66e << tablepath <<
"t10_e.dat";
2005 ex67 << tablepath <<
"t11.dat";
2007 ex67e << tablepath <<
"t11_e.dat";
2009 ex68 << tablepath <<
"t12.dat";
2011 ex68e << tablepath <<
"t12_e.dat";
2013 ex69 << tablepath <<
"t13.dat";
2015 ex69e << tablepath <<
"t13_e.dat";
2017 ex70 << tablepath <<
"t14.dat";
2019 ex70e << tablepath <<
"t14_e.dat";
2021 ex71 << tablepath <<
"t15.dat";
2023 ex71e << tablepath <<
"t15_e.dat";
2025 ex72 << tablepath <<
"t16.dat";
2027 ex72e << tablepath <<
"t16_e.dat";
2029 ex73 << tablepath <<
"t17.dat";
2031 ex73e << tablepath <<
"t17_e.dat";
2033 ex74 << tablepath <<
"t18.dat";
2035 ex74e << tablepath <<
"t18_e.dat";
2037 ex75 << tablepath <<
"t19.dat";
2039 ex75e << tablepath <<
"t19_e.dat";
2041 ex76 << tablepath <<
"t20.dat";
2043 ex76e << tablepath <<
"t20_e.dat";
2045 ex77 << tablepath <<
"t21.dat";
2047 ex77e << tablepath <<
"t21_e.dat";
2049 ex78 << tablepath <<
"t22.dat";
2051 ex78e << tablepath <<
"t22_e.dat";
2053 ex79 << tablepath <<
"t23.dat";
2055 ex79e << tablepath <<
"t23_e.dat";
2057 ex80 << tablepath <<
"t24.dat";
2059 ex80e << tablepath <<
"t24_e.dat";
2061 ex81 << tablepath <<
"t25.dat";
2063 ex81e << tablepath <<
"t25_e.dat";
2065 ex82 << tablepath <<
"t26.dat";
2067 ex82e << tablepath <<
"t26_e.dat";
2069 ex83 << tablepath <<
"t27.dat";
2071 ex83e << tablepath <<
"t27_e.dat";
2073 ex84 << tablepath <<
"t28.dat";
2075 ex84e << tablepath <<
"t28_e.dat";
2077 ex85 << tablepath <<
"t29.dat";
2079 ex85e << tablepath <<
"t29_e.dat";
2081 ex86 << tablepath <<
"t30.dat";
2083 ex86e << tablepath <<
"t30_e.dat";
2085 ex87 << tablepath <<
"t31.dat";
2087 ex87e << tablepath <<
"t31_e.dat";
2089 ex88 << tablepath <<
"t32.dat";
2091 ex88e << tablepath <<
"t32_e.dat";
2093 ex89 << tablepath <<
"t33.dat";
2095 ex89e << tablepath <<
"t33_e.dat";
2097 ex90 << tablepath <<
"t34.dat";
2099 ex90e << tablepath <<
"t34_e.dat";
2101 ex91 << tablepath <<
"t35.dat";
2103 ex91e << tablepath <<
"t35_e.dat";
2105 ex92 << tablepath <<
"t36.dat";
2107 ex92e << tablepath <<
"t36_e.dat";
2109 ex93 << tablepath <<
"t37.dat";
2111 ex93e << tablepath <<
"t37_e.dat";
2113 ex94 << tablepath <<
"t38.dat";
2115 ex94e << tablepath <<
"t38_e.dat";
2117 ex95 << tablepath <<
"t39.dat";
2119 ex95e << tablepath <<
"t39_e.dat";
2121 ex96 << tablepath <<
"t40.dat";
2123 ex96e << tablepath <<
"t40_e.dat";
2126 ex97 << tablepath <<
"CMS-PAS-HIG-16-034.dat";
2128 ex97e << tablepath <<
"CMS-PAS-HIG-16-034_e.dat";
2130 ex98 << tablepath <<
"14078150.dat";
2132 ex98e << tablepath <<
"14078150.dat";
2135 ex99 << tablepath <<
"14126663.dat";
2137 ex99e << tablepath <<
"14126663_e.dat";
2139 ex100 << tablepath <<
"151203704.dat";
2141 ex100e << tablepath <<
"151203704_e.dat";
2143 ex101 << tablepath <<
"150807774_a.dat";
2145 ex101e << tablepath <<
"150807774_a_e.dat";
2147 ex102 << tablepath <<
"150807774_b.dat";
2149 ex102e << tablepath <<
"150807774_b_e.dat";
2151 ex103 << tablepath <<
"ATLAS-CONF-2016-088.dat";
2153 ex103e << tablepath <<
"ATLAS-CONF-2016-088_e.dat";
2155 ex104 << tablepath <<
"ATLAS-CONF-2016-089.dat";
2157 ex104e << tablepath <<
"ATLAS-CONF-2016-089_e.dat";
2159 ex105 << tablepath <<
"ATLAS-CONF-2016-104_c.dat";
2161 ex105e << tablepath <<
"ATLAS-CONF-2016-104_c_e.dat";
2163 ex106 << tablepath <<
"CMS-PAS-HIG-16-031.dat";
2165 ex106e << tablepath <<
"CMS-PAS-HIG-16-031_e.dat";
2168 ex107 << tablepath <<
"CMS-PAS-HIG-17-002.dat";
2170 ex107e << tablepath <<
"CMS-PAS-HIG-17-002_e.dat";
2172 ex108 << tablepath <<
"170804188.dat";
2174 ex108e << tablepath <<
"170804188_e.dat";
2176 ex109 << tablepath <<
"171203143.dat";
2178 ex109e << tablepath <<
"171203143_e.dat";
2180 ex110 << tablepath <<
"171004960.dat";
2182 ex110e << tablepath <<
"171004960_e.dat";
2185 ex111 << tablepath <<
"171007235_b.dat";
2187 ex111e << tablepath <<
"171007235_b_e.dat";
2189 ex112 << tablepath <<
"170800212.dat";
2191 ex112e << tablepath <<
"170800212_e.dat";
2193 ex113 << tablepath <<
"171206386_a.dat";
2195 ex113e << tablepath <<
"171206386_a_e.dat";
2197 ex114 << tablepath <<
"171206386_b.dat";
2199 ex114e << tablepath <<
"171206386_b_e.dat";
2201 ex115 << tablepath <<
"170809638_a.dat";
2203 ex115e << tablepath <<
"170809638_a_e.dat";
2205 ex116 << tablepath <<
"170809638_b.dat";
2207 ex116e << tablepath <<
"170809638_b_e.dat";
2209 ex117 << tablepath <<
"170804445.dat";
2211 ex117e << tablepath <<
"170804445_e.dat";
2213 ex118 << tablepath <<
"CMS-PAS-B2G-16-023.dat";
2215 ex118e << tablepath <<
"CMS-PAS-B2G-16-023_e.dat";
2217 ex119 << tablepath <<
"CMS-PAS-HIG-16-001_a.dat";
2219 ex119e << tablepath <<
"CMS-PAS-HIG-16-001_a_e.dat";
2221 ex120 << tablepath <<
"CMS-PAS-HIG-16-001_b.dat";
2223 ex120e << tablepath <<
"CMS-PAS-HIG-16-001_b_e.dat";
2226 ex121 << tablepath <<
"160302991_a.dat";
2228 ex122 << tablepath <<
"160302991_b.dat";
2230 ex123 << tablepath <<
"160302991_c.dat";
2232 ex124 << tablepath <<
"160302991_d.dat";
2235 bsg1 << tablepath <<
"bsgammatable.dat";
2243 double params[] = {mass};
2259 double params[] = {mass};
2275 double params[] = {mass};
2291 double params[] = {mass};
2307 double params[] = {mass};
2323 double params[] = {mass};
2339 double params[] = {mass};
2355 double params[] = {mass};
2371 double params[] = {mass};
2377 double newResult = 0.0;
2378 if (mass>=20. && mass <=2000.) {
2390 double params[] = {mass};
2396 double newResult = 0.0;
2397 if (mass>=20. && mass <=2000.) {
2409 double params[] = {mass};
2415 double newResult = 0.0;
2416 if (mass>=20. && mass <=2000.) {
2428 double params[] = {mass};
2434 double newResult = 0.0;
2435 if (mass>=20. && mass <=2000.) {
2447 double params[] = {mass};
2453 double newResult = 0.0;
2454 if (mass>=20. && mass <=2000.) {
2466 double params[] = {mass};
2472 double newResult = 0.0;
2473 if (mass>=20. && mass <=2000.) {
2485 double params[] = {mass};
2491 double newResult = 0.0;
2492 if (mass>=20. && mass <=2000.) {
2504 double params[] = {mass};
2510 double newResult = 0.0;
2511 if (mass>=20. && mass <=2000.) {
2523 double params[] = {mass};
2529 double newResult = 0.0;
2530 if (mass>=20. && mass <=2000.) {
2542 double params[] = {mass};
2548 double newResult = 0.0;
2549 if (mass>=20. && mass <=2000.) {
2561 double params[] = {mass};
2567 double newResult = 0.0;
2568 if (mass>=20. && mass <=2000.) {
2580 double params[] = {mass};
2586 double newResult = 0.0;
2587 if (mass>=20. && mass <=2000.) {
2599 double params[] = {mass};
2605 double newResult = 0.0;
2606 if (mass>=20. && mass <=2000.) {
2618 double params[] = {mass};
2624 double newResult = 0.0;
2625 if (mass>=20. && mass <=2000.) {
2637 double params[] = {mass};
2643 double newResult = 0.0;
2644 if (mass>=20. && mass <=2000.) {
2656 double params[] = {mass};
2662 double newResult = 0.0;
2663 if (mass>=20. && mass <=2000.) {
2675 double params[] = {mass};
2681 double newResult = 0.0;
2682 if (mass>=20. && mass <=2000.) {
2694 double params[] = {mass};
2700 double newResult = 0.0;
2701 if (mass>=20. && mass <=2000.) {
2713 double params[] = {mHp,
logtb};
2719 double newResult = 0.0;
2720 if (mHp>=180. && mHp <=1400. && logtb>=-1. &&
logtb<=1.75) {
2732 double params[] = {mHp,
logtb};
2738 double newResult = 0.0;
2739 if (mHp>=180. && mHp <=2000. && logtb>=-1. &&
logtb<=1.75) {
2751 double params[] = {mass};
2767 double params[] = {mass};
2783 double params[] = {mass};
2799 double params[] = {mass};
2815 double params[] = {mass};
2831 double params[] = {mass};
2847 double params[] = {mass};
2863 double params[] = {mass};
2879 double params[] = {mass};
2895 double params[] = {mass};
2911 double params[] = {mass};
2927 double params[] = {mass};
2943 double params[] = {mass};
2959 double params[] = {mass};
2975 double params[] = {mass};
2991 double params[] = {mass};
3007 double params[] = {mass};
3023 double params[] = {mass};
3039 double params[] = {mass};
3057 double params[] = {mass};
3073 double params[] = {mass};
3089 double params[] = {mass};
3105 double params[] = {mass};
3121 double params[] = {mass};
3137 double params[] = {mass};
3153 double params[] = {mass};
3169 double params[] = {mass};
3185 double params[] = {mass};
3201 double params[] = {mass};
3217 double params[] = {mass};
3233 double params[] = {mass};
3249 double params[] = {mass};
3265 double params[] = {mass};
3281 double params[] = {mass};
3297 double params[] = {mass};
3313 double params[] = {mass};
3329 double params[] = {mass};
3345 double params[] = {mass};
3361 double params[] = {mass};
3377 double params[] = {mass};
3393 double params[] = {mass};
3409 double params[] = {mass};
3425 double params[] = {mass};
3441 double params[] = {mass};
3457 double params[] = {mass};
3473 double params[] = {mass};
3489 double params[] = {mass};
3505 double params[] = {mass};
3521 double params[] = {mass};
3537 double params[] = {mass};
3612 double params[] = {mass};
3628 double params[] = {mass};
3644 double params[] = {mass};
3660 double params[] = {mass};
3676 double params[] = {mass};
3692 double params[] = {mass};
3708 double params[] = {mA, mH};
3724 double params[] = {mA, mH};
3740 double params[] = {mA, mH};
3756 double params[] = {mA, mH};
3772 double params[] = {mass};
3788 double params[] = {mass};
3804 double params[] = {mass};
3820 double params[] = {mass};
3836 double params[] = {mass};
3852 double params[] = {mass};
3868 double params[] = {mass};
3884 double params[] = {mass};
3900 double params[] = {mass};
3916 double params[] = {mass};
3932 double params[] = {mass};
3948 double params[] = {mass};
3964 double params[] = {mass};
3980 double params[] = {mass};
3996 double params[] = {mass};
4012 double params[] = {mass};
4028 double params[] = {mass};
4044 double params[] = {mass};
4060 double params[] = {mass};
4076 double params[] = {mass};
4092 double params[] = {mass};
4108 double params[] = {mass};
4124 double params[] = {mass};
4140 double params[] = {mass};
4156 double params[] = {mass};
4172 double params[] = {mass};
4188 double params[] = {mass};
4204 double params[] = {mass};
4220 double params[] = {mass};
4236 double params[] = {mass};
4252 double params[] = {mass};
4268 double params[] = {mass};
4284 double params[] = {mass};
4300 double params[] = {mass};
4316 double params[] = {mass};
4332 double params[] = {mass};
4348 double params[] = {mass};
4364 double params[] = {mass};
4380 double params[] = {mass};
4396 double params[] = {mass};
4412 double params[] = {mass};
4428 double params[] = {mass};
4444 double params[] = {mass};
4460 double params[] = {mass};
4476 double params[] = {mass};
4492 double params[] = {mass};
4508 double params[] = {mass};
4524 double params[] = {mass};
4540 double params[] = {mass};
4556 double params[] = {mass};
4572 double params[] = {mass};
4588 double params[] = {mass};
4604 double params[] = {mass};
4620 double params[] = {mass};
4636 double params[] = {mass};
4652 double params[] = {mass};
4668 double params[] = {mass};
4684 double params[] = {mass};
4700 double params[] = {mass};
4716 double params[] = {mass};
4732 double params[] = {mass};
4748 double params[] = {mass};
4764 double params[] = {mass};
4780 double params[] = {mass};
4796 double params[] = {mass};
4812 double params[] = {mass};
4828 double params[] = {mass};
4844 double params[] = {mass};
4860 double params[] = {mass};
4876 double params[] = {mass};
4892 double params[] = {mass};
4908 double params[] = {mass};
4924 double params[] = {mass};
4940 double params[] = {mass};
4956 double params[] = {mass};
4972 double params[] = {mass};
4988 double params[] = {mass};
5004 double params[] = {mass};
5020 double params[] = {mass};
5036 double params[] = {mass};
5052 double params[] = {mass};
5068 double params[] = {mass};
5084 double params[] = {mass};
5100 double params[] = {mass};
5116 double params[] = {mass};
5132 double params[] = {mass};
5148 double params[] = {mass};
5164 double params[] = {mass};
5180 double params[] = {mass};
5196 double params[] = {mass};
5212 double params[] = {mass};
5228 double params[] = {mass};
5244 double params[] = {mass};
5260 double params[] = {mass};
5276 double params[] = {mass};
5292 double params[] = {mass};
5308 double params[] = {mass};
5312 return(
ie1c[NumPar][i] );
5324 double params[] = {mass};
5328 return(
ie1ec[NumPar][i] );
5340 double params[] = {mass};
5344 return(
ie2c[NumPar][i] );
5356 double params[] = {mass};
5360 return(
ie2ec[NumPar][i] );
5372 double params[] = {mass};
5376 return(
ie3c[NumPar][i] );
5388 double params[] = {mass};
5392 return(
ie3ec[NumPar][i] );
5404 double params[] = {mass};
5408 return(
ie4c[NumPar][i] );
5420 double params[] = {mass};
5424 return(
ie4ec[NumPar][i] );
5436 double params[] = {mass};
5440 return(
ie5c[NumPar][i] );
5452 double params[] = {mass};
5456 return(
ie5ec[NumPar][i] );
5468 double params[] = {mass};
5472 return(
ie6c[NumPar][i] );
5484 double params[] = {mass};
5488 return(
ie6ec[NumPar][i] );
5500 double params[] = {mass};
5504 return(
ie7c[NumPar][i] );
5516 double params[] = {mass};
5520 return(
ie7ec[NumPar][i] );
5532 double params[] = {mass};
5536 return(
ie8c[NumPar][i] );
5548 double params[] = {mass};
5552 return(
ie8ec[NumPar][i] );
5564 double params[] = {mass};
5568 return(
ie9c[NumPar][i] );
5580 double params[] = {mass};
5584 return(
ie9ec[NumPar][i] );
5596 double params[] = {mass};
5600 return(
ie10c[NumPar][i] );
5612 double params[] = {mass};
5616 return(
ie10ec[NumPar][i] );
5628 double params[] = {mass};
5632 return(
ie11c[NumPar][i] );
5644 double params[] = {mass};
5648 return(
ie11ec[NumPar][i] );
5660 double params[] = {mass};
5664 return(
ie12c[NumPar][i] );
5676 double params[] = {mass};
5680 return(
ie12ec[NumPar][i] );
5692 double params[] = {mass};
5696 return(
ie13c[NumPar][i] );
5708 double params[] = {mass};
5712 return(
ie13ec[NumPar][i] );
5724 double params[] = {mass};
5728 return(
ie14c[NumPar][i] );
5740 double params[] = {mass};
5744 return(
ie14ec[NumPar][i] );
5756 double params[] = {mass};
5760 return(
ie15c[NumPar][i] );
5772 double params[] = {mass};
5776 return(
ie15ec[NumPar][i] );
5788 double params[] = {mass};
5792 return(
ie16c[NumPar][i] );
5804 double params[] = {mass};
5808 return(
ie16ec[NumPar][i] );
5820 double params[] = {mass};
5824 return(
ie17c[NumPar][i] );
5836 double params[] = {mass};
5840 return(
ie17ec[NumPar][i] );
5852 double params[] = {mass};
5856 return(
ie18c[NumPar][i] );
5868 double params[] = {mass};
5872 return(
ie18ec[NumPar][i] );
5884 double params[] = {mass};
5888 return(
ie19c[NumPar][i] );
5900 double params[] = {mass};
5904 return(
ie19ec[NumPar][i] );
5916 double params[] = {mass};
5920 return(
ie20c[NumPar][i] );
5932 double params[] = {mass};
5936 return(
ie20ec[NumPar][i] );
5948 double params[] = {mass};
5952 return(
ie21c[NumPar][i] );
5964 double params[] = {mass};
5968 return(
ie21ec[NumPar][i] );
5980 double params[] = {mass};
5984 return(
ie22c[NumPar][i] );
5996 double params[] = {mass};
6000 return(
ie22ec[NumPar][i] );
6012 double params[] = {mass};
6016 return(
ie23c[NumPar][i] );
6028 double params[] = {mass};
6032 return(
ie23ec[NumPar][i] );
6044 double params[] = {mass};
6048 return(
ie24c[NumPar][i] );
6060 double params[] = {mass};
6064 return(
ie24ec[NumPar][i] );
6076 double params[] = {mass};
6080 return(
ie25c[NumPar][i] );
6092 double params[] = {mass};
6096 return(
ie25ec[NumPar][i] );
6108 double params[] = {mass};
6112 return(
ie26c[NumPar][i] );
6124 double params[] = {mass};
6128 return(
ie26ec[NumPar][i] );
6140 double params[] = {mass};
6144 return(
ie27c[NumPar][i] );
6156 double params[] = {mass};
6160 return(
ie27ec[NumPar][i] );
6172 double params[] = {mass};
6176 return(
ie28c[NumPar][i] );
6188 double params[] = {mass};
6192 return(
ie28ec[NumPar][i] );
6204 double params[] = {mass};
6208 return(
ie29c[NumPar][i] );
6220 double params[] = {mass};
6224 return(
ie29ec[NumPar][i] );
6236 double params[] = {mass};
6240 return(
ie30c[NumPar][i] );
6252 double params[] = {mass};
6256 return(
ie30ec[NumPar][i] );
6268 double params[] = {mass};
6272 return(
ie31c[NumPar][i] );
6284 double params[] = {mass};
6288 return(
ie31ec[NumPar][i] );
6300 double params[] = {mass};
6304 return(
ie32c[NumPar][i] );
6316 double params[] = {mass};
6320 return(
ie32ec[NumPar][i] );
6332 double params[] = {mass};
6336 return(
ie33c[NumPar][i] );
6348 double params[] = {mass};
6352 return(
ie33ec[NumPar][i] );
6364 double params[] = {mass};
6368 return(
ie34c[NumPar][i] );
6380 double params[] = {mass};
6384 return(
ie34ec[NumPar][i] );
6396 double params[] = {mass};
6400 return(
ie35c[NumPar][i] );
6412 double params[] = {mass};
6416 return(
ie35ec[NumPar][i] );
6428 double params[] = {mass};
6432 return(
ie36c[NumPar][i] );
6444 double params[] = {mass};
6448 return(
ie36ec[NumPar][i] );
6460 double params[] = {mass};
6464 return(
ie37c[NumPar][i] );
6476 double params[] = {mass};
6480 return(
ie37ec[NumPar][i] );
6492 double params[] = {mass};
6496 return(
ie38c[NumPar][i] );
6508 double params[] = {mass};
6512 return(
ie38ec[NumPar][i] );
6524 double params[] = {mass};
6528 return(
ie39c[NumPar][i] );
6540 double params[] = {mass};
6544 return(
ie39ec[NumPar][i] );
6556 double params[] = {mass};
6560 return(
ie40c[NumPar][i] );
6572 double params[] = {mass};
6576 return(
ie40ec[NumPar][i] );
6588 double params[] = {mass};
6604 double params[] = {mass};
6620 double params[] = {mass};
6636 double params[] = {mass};
6652 double params[] = {mass};
6668 double params[] = {mass};
6684 double params[] = {mass};
6700 double params[] = {mass};
6716 double params[] = {mass};
6732 double params[] = {mass};
6748 double params[] = {mass};
6764 double params[] = {mass};
6780 double params[] = {mass};
6796 double params[] = {mass};
6812 double params[] = {mass};
6828 double params[] = {mass};
6844 double params[] = {
logtb, logmHp};
6860 std::ifstream INfile;
6861 std::string lineTab;
6862 INfile.open( filename.c_str() );
6864 std::cout<<
"error: in THDMcache, table doesn't exist!"<<std::endl;
6867 gslpp::matrix<double> arrayTab(rowN, colN, 0.);
6872 while(INfile.good()){
6873 while(getline(INfile, lineTab)){
6874 if( lineTab[0]==
'#' )
continue;
6876 std::istringstream streamTab(lineTab);
6878 while(streamTab >>v){
6879 arrayTab.assign(a,b,v);
6896 int rowN=arrayTab.size_i();
6898 double xmin = arrayTab(0,0);
6899 double xmax = arrayTab(rowN-1,0);
6900 double interval = arrayTab(1,0)-arrayTab(0,0);
6901 int Nintervals = (x-xmin)/interval;
6913 y =(arrayTab(Nintervals+1,1)-arrayTab(Nintervals,1))/(arrayTab(Nintervals+1,0)
6914 -arrayTab(Nintervals,0))*(x-arrayTab(Nintervals,0))+arrayTab(Nintervals,1);
6923 int rowN=arrayTab.size_i();
6925 double xmin = arrayTab(0,0);
6926 double xmax = arrayTab(rowN-1,0);
6927 double ymin = arrayTab(0,1);
6928 double ymax = arrayTab(rowN-1,1);
6929 double intervalx = arrayTab(1,0)-arrayTab(0,0);
6932 while(arrayTab(i,1)-arrayTab(i-1,1)==0&&i<30000);
6933 double intervaly = arrayTab(i,1)-arrayTab(i-1,1);
6934 int Nintervalsx = (x-xmin)/intervalx;
6935 int Nintervalsy = (y-ymin)/intervaly;
6936 if(x<xmin||x>xmax||y<ymin||y>ymax){
6941 double x1=arrayTab(i*Nintervalsy+Nintervalsx,0);
6942 double x2=arrayTab(i*Nintervalsy+Nintervalsx+1,0);
6943 double y1=arrayTab(i*Nintervalsy+Nintervalsx,1);
6944 double y2=arrayTab(i*(Nintervalsy+1)+Nintervalsx,1);
6945 return (arrayTab(i*Nintervalsy+Nintervalsx,2) * (x2-x) * (y2-y)
6946 +arrayTab(i*Nintervalsy+Nintervalsx+1,2) * (x-x1) * (y2-y)
6947 +arrayTab(i*(Nintervalsy+1)+Nintervalsx,2) * (x2-x) * (y-y1)
6948 +arrayTab(i*(Nintervalsy+1)+Nintervalsx+1,2) * (x-x1) * (y-y1))
6953double THDMcache::ghHpHm(
const double mHp2,
const double tanb,
const double m12_2,
const double bma,
const double mHl2,
const double vev)
const {
6969double THDMcache::g_HH_HpHm(
const double mHp2,
const double mHh2,
const double tanb,
const double m12_2,
const double bma,
const double vev)
const {
6985gslpp::complex
THDMcache::I_h_U(
const double mHl2,
const double Mu,
const double Mc,
const double Mt)
const {
6996 gslpp::complex newResult = -(8./3.)*(TAUu*(1.0+(1.0-TAUu)*
f_func(TAUu))
6997 +TAUc*(1.0+(1.0-TAUc)*
f_func(TAUc))+TAUt*(1.0+(1.0-TAUt)*
f_func(TAUt)));
7013 gslpp::complex newResult = -(8./3.)*(TAUc*(1.0+(1.0-TAUc)*
f_func(TAUc))
7014 +TAUt*(1.0+(1.0-TAUt)*
f_func(TAUt)));
7022 double params[] = {
mA2,
Mc,
Mt};
7030 gslpp::complex newResult = -(8./3.)*(TAUc*
f_func(TAUc)+TAUt*
f_func(TAUt));
7036gslpp::complex
THDMcache::I_h_D(
const double mHl2,
const double Md,
const double Ms,
const double Mb)
const {
7047 gslpp::complex newResult = -(2./3.)*(TAUd*(1.0+(1.0-TAUd)*
f_func(TAUd))
7048 +TAUs*(1.0+(1.0-TAUs)*
f_func(TAUs))+TAUb*(1.0+(1.0-TAUb)*
f_func(TAUb)));
7064 gslpp::complex newResult = -(2./3.)*(TAUs*(1.0+(1.0-TAUs)*
f_func(TAUs))
7065 +TAUb*(1.0+(1.0-TAUb)*
f_func(TAUb)));
7073 double params[] = {
mA2,
Ms,
Mb};
7081 gslpp::complex newResult = -(2./3.)*(TAUs*
f_func(TAUs)+TAUb*
f_func(TAUb));
7087gslpp::complex
THDMcache::I_h_L(
const double mHl2,
const double Me,
const double Mmu,
const double Mtau)
const {
7098 gslpp::complex newResult = -2.0*(TAUe*(1.0+(1.0-TAUe)*
f_func(TAUe))
7099 +TAUmu*(1.0+(1.0-TAUmu)*
f_func(TAUmu))
7100 +TAUtau*(1.0+(1.0-TAUtau)*
f_func(TAUtau)));
7116 gslpp::complex newResult = -2.0*(TAUmu*(1.0+(1.0-TAUmu)*
f_func(TAUmu))+
7117 TAUtau*(1.0+(1.0-TAUtau)*
f_func(TAUtau)));
7133 gslpp::complex newResult = -2.0*(TAUmu*
f_func(TAUmu)+TAUtau*
f_func(TAUtau));
7141 double params[] = {mH,
MW};
7147 double TAUw=4.0*
MW*
MW/(mH*mH);
7148 gslpp::complex newResult = 2.0 + 3.0*TAUw + 3.0*TAUw*(2.0-TAUw)*
f_func(TAUw);
7156 double params[] = {
mHp2, mH};
7162 double TAUhp=4.0*
mHp2/(mH*mH);
7163 gslpp::complex newResult = -TAUhp*(1.0-TAUhp*
f_func(TAUhp));
7169gslpp::complex
THDMcache::A_h_U(
const double mHl2,
const double cW2,
const double Mu,
const double Mc,
const double Mt,
const double MZ)
const {
7184 gslpp::complex newResult = -4.0*(1.0/2.0-4.0/3.0*sW2)*(
Int1(TAUu,LAMu)+
Int1(TAUc,LAMc)
7191gslpp::complex
THDMcache::A_HH_U(
const double mHh2,
const double cW2,
const double Mc,
const double Mt,
const double MZ)
const {
7204 gslpp::complex newResult = -4.0*(1.0/2.0-4.0/3.0*sW2)*(
Int1(TAUc,LAMc)-
Int2(TAUc,LAMc)
7211gslpp::complex
THDMcache::A_A_U(
const double mA2,
const double cW2,
const double Mc,
const double Mt,
const double MZ)
const {
7224 gslpp::complex newResult = -4.0*(1.0/2.0-4.0/3.0*sW2)*(-
Int2(TAUc,LAMc)-
Int2(TAUt,LAMt))/sqrt(sW2*
cW2);
7230gslpp::complex
THDMcache::A_h_D(
const double mHl2,
const double cW2,
const double Md,
const double Ms,
const double Mb,
const double MZ)
const {
7245 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0/3.0*sW2)*(
Int1(TAUd,LAMd)+
Int1(TAUs,LAMs)
7252gslpp::complex
THDMcache::A_HH_D(
const double mHh2,
const double cW2,
const double Ms,
const double Mb,
const double MZ)
const {
7265 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0/3.0*sW2)*(
Int1(TAUs,LAMs)-
Int2(TAUs,LAMs)
7272gslpp::complex
THDMcache::A_A_D(
const double mA2,
const double cW2,
const double Ms,
const double Mb,
const double MZ)
const {
7285 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0/3.0*sW2)*(-
Int2(TAUs,LAMs)-
Int2(TAUb,LAMb))/sqrt(sW2*
cW2);
7291gslpp::complex
THDMcache::A_h_L(
const double mHl2,
const double cW2,
const double Me,
const double Mmu,
const double Mtau,
const double MZ)
const {
7306 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0*sW2)*(
Int1(TAUe,LAMe)+
Int1(TAUmu,LAMmu)
7308 -
Int2(TAUtau,LAMtau));
7314gslpp::complex
THDMcache::A_HH_L(
const double mHh2,
const double cW2,
const double Mmu,
const double Mtau,
const double MZ)
const {
7327 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0*sW2)*(
Int1(TAUmu,LAMmu)-
Int2(TAUmu,LAMmu)
7328 +
Int1(TAUtau,LAMtau)-
Int2(TAUtau,LAMtau));
7334gslpp::complex
THDMcache::A_A_L(
const double mA2,
const double cW2,
const double Mmu,
const double Mtau,
const double MZ)
const {
7347 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0*sW2)*(-
Int2(TAUmu,LAMmu)-
Int2(TAUtau,LAMtau))/sqrt(sW2*
cW2);
7353gslpp::complex
THDMcache::A_H_W(
const double mH,
const double cW2,
const double MW,
const double MZ)
const {
7355 double params[] = {mH,
cW2,
MW,
MZ};
7361 double TAUw=4.0*
MW*
MW/(mH*mH);
7364 gslpp::complex newResult = -sqrt(
cW2/sW2)*(4.0*(3.0-sW2/
cW2)*
Int2(TAUw,LAMw)
7365 +((1.0+2.0/TAUw)*sW2/
cW2-(5.0+2.0/TAUw))*
Int1(TAUw,LAMw));
7371gslpp::complex
THDMcache::A_H_Hp(
const double mHp2,
const double mH,
const double cW2,
const double MZ)
const {
7379 double TAUhp=4.0*
mHp2/(mH*mH);
7382 gslpp::complex newResult = (1.0-2.0*sW2)/sqrt(
cW2*sW2)*
Int1(TAUhp,LAMhp);
7390 gslpp::complex z = -gslpp::complex::i()*M_PI;
7391 return -pow(log((1.0+sqrt(1.0-x))/(1.0-sqrt(1.0-x)))+z,2)/4.0;
7394 return pow(asin(sqrt(1.0/x)),2);
7402 gslpp::complex z = -gslpp::complex::i()*M_PI;
7403 gslpp::complex gs1 = sqrt(1.0-x)*(log((1.0+sqrt(1.0-x))/(1.0-sqrt(1.0-x)))+z)/2.0;
7407 gslpp::complex gg1 = sqrt(x-1.0)*asin(sqrt(1.0/x));
7415 return tau*lambda/(tau-lambda)/2.0+tau*tau*lambda*lambda/((tau-lambda)
7416 *(tau-lambda))/2.0*(
f_func(tau)-
f_func(lambda))+tau*tau*lambda/((tau-lambda)
7423 return -tau*lambda/(tau-lambda)/2.0*(
f_func(tau)-
f_func(lambda));
7435 double params[] = {a2, b2, c2};
7442 double newResult = 0.0;
7443 double x = (a2-b2-c2)*(a2-b2-c2)-4.0*b2*c2;
7444 if(x>0) newResult = sqrt(std::fabs(x/a2))/2.0;
7466 double sin_ba=sin(
bma);
7468 double cosb=1.0/sqrt(1.0+
tanb*
tanb);
7469 double sina=sinb*cos(
bma)-cosb*sin(
bma);
7470 double cosa=cosb*cos(
bma)+sinb*sin(
bma);
7477 double BrSM_htoWW = 2.15e-1;
7478 double BrSM_htoZZ = 2.64e-2;
7479 double BrSM_htogg = 8.57e-2;
7480 double BrSM_htoZga = 1.54e-3;
7481 double BrSM_htocc = 2.91e-2;
7506 rh_QuQu=cosa*cosa/(sinb*sinb);
7507 rh_VV=sin_ba*sin_ba;
7516 gslpp::complex I_h_F=0.0;
7521 gslpp::complex I_h_W=sin_ba*I_hSM_W;
7524 double ABSgagaTHDM=0.0;
7525 double ABSgagaSM=0.0;
7529 gslpp::complex A_h_F = 0.0;
7534 gslpp::complex A_h_W = sin_ba*A_hSM_W;
7537 double ABSZgaTHDM=0.0;
7538 double ABSZgaSM=0.0;
7541 rh_gg=cosa/sinb*cosa/sinb;
7543 rh_ll=cosa/sinb*cosa/sinb;
7544 I_h_F=cosa/sinb*(fermU+fermD+fermL);
7545 A_h_F = cosa/sinb*(A_h_Ux+A_h_Dx+A_h_Lx)/sqrt(sW2*
cW2);
7548 rh_gg=-cosa/sinb*sina/cosb+(cosa/sinb+sina/cosb)
7549 *(Sigmaggh_tt8*cosa/sinb+Sigmaggh_bb8*sina/cosb)/SigmaggF8;
7551 rh_ll=sina/cosb*sina/cosb;
7552 I_h_F=cosa/sinb*fermU -sina/cosb*(fermD+fermL);
7553 A_h_F = (cosa/sinb*A_h_Ux-sina/cosb*(A_h_Dx+A_h_Lx))/sqrt(sW2*
cW2);
7556 rh_gg=cosa/sinb*cosa/sinb;
7558 rh_ll=sina/cosb*sina/cosb;
7559 I_h_F = cosa/sinb*(fermU+fermD) -sina/cosb*fermL;
7560 A_h_F = (cosa/sinb*(A_h_Ux+A_h_Dx)-sina/cosb*A_h_Lx)/sqrt(sW2*
cW2);
7563 rh_gg=-cosa/sinb*sina/cosb+(cosa/sinb+sina/cosb)
7564 *(Sigmaggh_tt8*cosa/sinb+Sigmaggh_bb8*sina/cosb)/SigmaggF8;
7566 rh_ll=cosa/sinb*cosa/sinb;
7567 I_h_F = cosa/sinb*(fermU+fermL) -sina/cosb*fermD;
7568 A_h_F = (cosa/sinb*(A_h_Ux+A_h_Lx)-sina/cosb*A_h_Dx)/sqrt(sW2*
cW2);
7571 throw std::runtime_error(
"modelflag can be only any of \"type1\", \"type2\", \"typeX\" or \"typeY\"");
7574 ABSgagaTHDM=(I_h_F+I_h_W+I_h_Hp).abs2();
7575 ABSgagaSM=(fermU+fermL+fermD+I_hSM_W).abs2();
7576 rh_gaga=ABSgagaTHDM/ABSgagaSM;
7578 ABSZgaTHDM=(A_h_F+A_h_W+A_h_Hp).abs2();
7579 ABSZgaSM=((A_h_Ux+A_h_Lx+A_h_Dx)/sqrt(sW2*
cW2)+A_hSM_W).abs2();
7580 rh_Zga=ABSZgaTHDM/ABSZgaSM;
7588 pph13 = (SigmaggF13*
rh_gg + SigmaVBFVh13*
rh_VV + Sigmatth13*
rh_QuQu + Sigmabbh13*
rh_QdQd)/(SigmaggF13 + SigmaVBFVh13 + Sigmatth13 + Sigmabbh13);
7608 double GF=1/(sqrt(2.0)*
vev*
vev);
7610 double mHh=sqrt(
mHh2);
7611 double sin_ba=sin(
bma);
7613 double cosb=1.0/sqrt(1.0+
tanb*
tanb);
7614 double sina=sinb*cos(
bma)-cosb*sin(
bma);
7615 double cosa=cosb*cos(
bma)+sinb*sin(
bma);
7616 double cos_2b=cosb*cosb-sinb*sinb;
7617 double cos_ba=cos(
bma);
7631 double rHH_QuQu=sina/sinb*sina/sinb;
7632 double rHH_QdQd=0.0;
7646 gslpp::complex I_HH_F=0.0;
7650 gslpp::complex I_HH_W=cos_ba*
I_H_W(mHh,
MW);
7657 gslpp::complex A_HH_F = 0.0;
7667 rHH_gg=sina/sinb*sina/sinb;
7668 rHH_QdQd=sina/sinb*sina/sinb;
7669 rHH_ll=sina/sinb*sina/sinb;
7670 I_HH_F=sina/sinb*(I_HH_Ux+I_HH_Dx+I_HH_Lx);
7671 A_HH_F = sina/sinb*(A_HH_Ux+A_HH_Dx+A_HH_Lx)/sqrt(sW2*
cW2);
7674 rHH_gg=sina/sinb*cosa/cosb+rSigmaggH_t8*sina/sinb*(sina/sinb-cosa/cosb)
7675 +rSigmaggH_b8*cosa/cosb*(cosa/cosb-sina/sinb);
7676 rHH_QdQd=cosa/cosb*cosa/cosb;
7677 rHH_ll=cosa/cosb*cosa/cosb;
7678 I_HH_F=sina/sinb*I_HH_Ux+cosa/cosb*(I_HH_Dx+I_HH_Lx);
7679 A_HH_F = (sina/sinb*A_HH_Ux+cosa/cosb*(A_HH_Dx+A_HH_Lx))/sqrt(sW2*
cW2);
7682 rHH_gg=sina/sinb*sina/sinb;
7683 rHH_QdQd=sina/sinb*sina/sinb;
7684 rHH_ll=cosa/cosb*cosa/cosb;
7685 I_HH_F=sina/sinb*(I_HH_Ux+I_HH_Dx)+cosa/cosb*I_HH_Lx;
7686 A_HH_F = (sina/sinb*(A_HH_Ux+A_HH_Dx)+cosa/cosb*A_HH_Lx)/sqrt(sW2*
cW2);
7689 rHH_gg=sina/sinb*cosa/cosb+rSigmaggH_t8*sina/sinb*(sina/sinb-cosa/cosb)
7690 +rSigmaggH_b8*cosa/cosb*(cosa/cosb-sina/sinb);
7691 rHH_QdQd=cosa/cosb*cosa/cosb;
7692 rHH_ll=sina/sinb*sina/sinb;
7693 I_HH_F=sina/sinb*(I_HH_Ux+I_HH_Lx)+cosa/cosb*I_HH_Dx;
7694 A_HH_F = (sina/sinb*(A_HH_Ux+A_HH_Lx)+cosa/cosb*A_HH_Dx)/sqrt(sW2*
cW2);
7697 throw std::runtime_error(
"modelflag can be only any of \"type1\", \"type2\", \"typeX\" or \"typeY\"");
7701 double Gamma_Hgaga=GF*
Ale*
Ale*mHh*mHh*mHh/(sqrt(2.0)*128.0*M_PI*M_PI*M_PI)
7702 *(I_HH_F+I_HH_W+I_HH_Hp).abs2();
7705 double Gamma_HZga=
HSTheta(mHh-
MZ)*GF*
Ale*
Ale*mHh*mHh*mHh/(sqrt(2.0)*64.0*M_PI*M_PI*M_PI)
7706 *(1.0-
MZ*
MZ/(mHh*mHh))*(1.0-
MZ*
MZ/(mHh*mHh))*(1.0-
MZ*
MZ/(mHh*mHh))
7707 *(A_HH_F+A_HH_W+A_HH_Hp).abs2();
7710 double Gamma_Hgg=
rHH_gg*GF*
Als*
Als*mHh*mHh*mHh/(sqrt(2.0)*16.0*M_PI*M_PI*M_PI)
7711 *(9.0/4.0)*(I_HH_Ux/4.0+I_HH_Dx).abs2();
7721 if (mHh>=20. && mHh <=2000.) {
7748 *std::fabs((cos_ba*cos_ba/(4.0*sinb*cosb*sinb*cosb)
7750 *sina*cosa,2))/(
vev*
vev))/(8.0*mHh*M_PI);
7753 *std::fabs(pow(cos_ba*(
mHh2 + 2.0*
mHp2 - 2.0*
m12_2/sinb/cosb)
7757 double GammaHAA=
HSTheta(mHh-2.0*sqrt(
mA2))*sqrt(std::fabs(1.0 - (4.0*
mA2)/
mHh2))
7758 *std::fabs(pow(cos_ba*(2.0*
mA2 +
mHh2 - 2.0*
m12_2/sinb/cosb)
7763 *sin_ba*sin_ba/(2.0*M_PI*
vev*
vev);
7767 GammaHtot= ((BrSM_Htott+BrSM_Htocc)*rHH_QuQu
7768 +BrSM_Htobb*rHH_QdQd
7769 +(BrSM_Htotautau+BrSM_Htomumu)*rHH_ll
7771 +Gamma_Hgaga+Gamma_HZga+Gamma_Hgg
7772 +GammaHhh+GammaHHpHm+GammaHAA+GammaHAZ+GammaHHpW;
7808 double mHh=sqrt(
mHh2);
7809 double mA=sqrt(
mA2);
7811 double Br_Ztoee=0.03363;
7812 double Br_Ztomumu=0.03366;
7813 double Br_Ztotautau=0.0337;
7814 double Br_Ztoinv=0.2;
7815 double Br_Wtoenu=0.1071;
7816 double Br_Wtomunu=0.1063;
7817 double Br_Wtotaunu=0.1138;
7818 double Br_tautoleptons=0.3521;
7919 +
THDM_BR_h_ZZ*2.0*Br_Ztoinv*(Br_Ztoee+Br_Ztomumu+Br_Ztotautau*Br_tautoleptons*Br_tautoleptons));
8020 double nftos=1.95996398454;
8022 if(mHh>=65.0 && mHh<90.0)
8026 if(mA>=50.0 && mA<=1000.0)
8032 else if(mHh>=90.0 && mHh<100.0)
8044 if(mA>=50.0 && mA<=1000.0)
8055 else if(mHh>=100.0 && mHh<130.0)
8067 if(mA>=50.0 && mA<=1000.0)
8078 else if(mHh>=130.0 && mHh<140.0)
8092 if(mA>=50.0 && mA<=1000.0)
8107 else if(mHh>=140.0 && mHh<145.0)
8125 if(mA>=50.0 && mA<=1000.0)
8140 else if(mHh>=145.0 && mHh<150.0)
8160 if(mA>=50.0 && mA<=1000.0)
8175 else if(mHh>=150.0 && mHh<175.0)
8195 if(mA>=50.0 && mA<=1000.0)
8210 else if(mHh>=175.0 && mHh<200.0)
8230 if(mA>=40.0 && mA<=910.0)
8235 if(mA>=50.0 && mA<=1000.0)
8250 else if(mHh>=200.0 && mHh<220.0)
8274 if(mA>=40.0 && mA<=910.0)
8279 if(mA>=50.0 && mA<=1000.0)
8314 else if(mHh>=220.0 && mHh<250.0)
8338 if(mA>=40.0 && mA<=910.0)
8343 if(mA>=50.0 && mA<=1000.0)
8378 else if(mHh>=250.0 && mHh<260.0)
8408 if(mA>=40.0 && mA<=910.0)
8413 if(mA>=50.0 && mA<=1000.0)
8458 else if(mHh>=260.0 && mHh<270.0)
8494 if(mA>=40.0 && mA<=910.0)
8499 if(mA>=50.0 && mA<=1000.0)
8552 else if(mHh>=270.0 && mHh<275.0)
8588 if(mA>=40.0 && mA<=910.0)
8593 if(mA>=50.0 && mA<=1000.0)
8646 else if(mHh>=275.0 && mHh<300.0)
8684 if(mA>=40.0 && mA<=910.0)
8689 if(mA>=50.0 && mA<=1000.0)
8744 else if(mHh>=300.0 && mHh<350.0)
8788 if(mA>=40.0 && mA<=910.0)
8793 if(mA>=50.0 && mA<=1000.0)
8870 else if(mHh>=350.0 && mHh<400.0)
8912 if(mA>=40.0 && mA<=910.0)
8917 if(mA>=50.0 && mA<=1000.0)
8996 else if(mHh>=400.0 && mHh<500.0)
9038 if(mA>=40.0 && mA<=910.0)
9043 if(mA>=50.0 && mA<=1000.0)
9124 else if(mHh>=500.0 && mHh<550.0)
9166 if(mA>=40.0 && mA<=910.0)
9171 if(mA>=50.0 && mA<=1000.0)
9256 else if(mHh>=550.0 && mHh<600.0)
9300 if(mA>=40.0 && mA<=910.0)
9305 if(mA>=50.0 && mA<=1000.0)
9392 else if(mHh>=600.0 && mHh<650.0)
9432 if(mA>=40.0 && mA<=910.0)
9437 if(mA>=50.0 && mA<=1000.0)
9522 else if(mHh>=650.0 && mHh<760.0)
9564 if(mA>=40.0 && mA<=910.0)
9569 if(mA>=50.0 && mA<=1000.0)
9656 else if(mHh>=760.0 && mHh<850.0)
9698 if(mA>=40.0 && mA<=910.0)
9703 if(mA>=50.0 && mA<=1000.0)
9790 else if(mHh>=850.0 && mHh<900.0)
9826 if(mA>=40.0 && mA<=910.0)
9831 if(mA>=50.0 && mA<=1000.0)
9918 else if(mHh>=900.0 && mHh<1000.0)
9952 if(mA>=40.0 && mA<=910.0)
9957 if(mA>=50.0 && mA<=1000.0)
10034 else if(mHh>=1000.0 && mHh<1100.0)
10110 else if(mHh>=1100.0 && mHh<1200.0)
10182 else if(mHh>=1200.0 && mHh<1500.0)
10248 else if(mHh>=1500.0 && mHh<1600.0)
10310 else if(mHh>=1600.0 && mHh<2000.0)
10370 else if(mHh>=2000.0 && mHh<2250.0)
10426 else if(mHh>=2250.0 && mHh<2400.0)
10478 else if(mHh>=2400.0 && mHh<2500.0)
10528 else if(mHh>=2500.0 && mHh<2530.0)
10574 else if(mHh>=2530.0 && mHh<2700.0)
10616 else if(mHh>=2700.0 && mHh<3000.0)
10656 else if(mHh>=3000.0 && mHh<3200.0)
10669 else if(mHh>=3200.0 && mHh<4000.0)
10682 double GF=1/(sqrt(2.0)*
vev*
vev);
10683 double mHh=sqrt(
mHh2);
10684 double mA=sqrt(
mA2);
10686 double cosb=1.0/sqrt(1.0+
tanb*
tanb);
10687 double mHp=sqrt(
mHp2);
10688 double sin_ba=sin(
bma);
10689 double cos_ba=cos(
bma);
10699 double rA_QuQu=(cosb*cosb)/(sinb*sinb);
10700 double rA_QdQd=0.0;
10713 gslpp::complex I_A_F=0.0;
10720 gslpp::complex A_A_F = 0.0;
10726 rA_gg=-cosb/sinb*cosb/sinb+2.0*cosb/sinb*cosb/sinb*(rSigmaggA_t8+rSigmaggA_b8);
10727 rA_QdQd=cosb/sinb*cosb/sinb;
10728 rA_ll=cosb/sinb*cosb/sinb;
10729 I_A_F=cosb/sinb*(I_A_Ux-I_A_Dx-I_A_Lx);
10730 A_A_F=cosb/sinb*(A_A_Ux-A_A_Dx-A_A_Lx);
10733 rA_gg= 1.0+(cosb/sinb-sinb/cosb)*(rSigmaggA_t8*cosb/sinb-rSigmaggA_b8*sinb/cosb);
10734 rA_QdQd=sinb/cosb*sinb/cosb;
10735 rA_ll=sinb/cosb*sinb/cosb;
10736 I_A_F=cosb/sinb*I_A_Ux+sinb/cosb*(I_A_Dx+I_A_Lx);
10737 A_A_F=cosb/sinb*A_A_Ux+sinb/cosb*(A_A_Dx+A_A_Lx);
10740 rA_gg=-cosb/sinb*cosb/sinb+2.0*cosb/sinb*cosb/sinb*(rSigmaggA_t8+rSigmaggA_b8);
10741 rA_QdQd=cosb/sinb*cosb/sinb;
10742 rA_ll=sinb/cosb*sinb/cosb;
10743 I_A_F=cosb/sinb*(I_A_Ux-I_A_Dx)+sinb/cosb*I_A_Lx;
10744 A_A_F=cosb/sinb*(A_A_Ux-A_A_Dx)+sinb/cosb*A_A_Lx;
10747 rA_gg=1.0+(cosb/sinb-sinb/cosb)*(rSigmaggA_t8*cosb/sinb-rSigmaggA_b8*sinb/cosb);
10748 rA_QdQd=sinb/cosb*sinb/cosb;
10749 rA_ll=cosb/sinb*cosb/sinb;
10750 I_A_F=cosb/sinb*(I_A_Ux-I_A_Lx)+sinb/cosb*I_A_Dx;
10751 A_A_F=cosb/sinb*(A_A_Ux-A_A_Lx)+sinb/cosb*A_A_Dx;
10754 throw std::runtime_error(
"modelflag can be only any of \"type1\", \"type2\", \"typeX\" or \"typeY\"");
10758 double Gamma_Agaga=GF*
Ale*
Ale*mA*mA*mA/(sqrt(2.0)*128.0*M_PI*M_PI*M_PI)
10762 *(1.0-
MZ*
MZ/(mA*mA))*(1.0-
MZ*
MZ/(mA*mA))*(1.0-
MZ*
MZ/(mA*mA))
10765 double Gamma_Agg=
rA_gg*GF*
Als*
Als*mA*mA*mA/(sqrt(2.0)*16.0*M_PI*M_PI*M_PI)
10766 *(9.0/4.0)*(I_A_Ux/4.0+I_A_Dx).abs2();
10789 *sin_ba*sin_ba/(2.0*M_PI*
vev*
vev);
10792 *cos_ba*cos_ba/(2.0*M_PI*
vev*
vev);
10797 GammaAtot= ((BrSM_Atott+BrSM_Atocc)*rA_QuQu
10798 +BrSM_Atobb*rA_QdQd
10799 +(BrSM_Atotautau+BrSM_Atomumu)*rA_ll)*GammaAtotSM
10800 +Gamma_Agaga+
Gamma_AZga+Gamma_Agg+GammaAHZ+GammaAhZ+GammaAHpW;
10813 double mHh=sqrt(
mHh2);
10814 double mA=sqrt(
mA2);
10816 double Br_Ztoee=0.03363;
10817 double Br_Ztomumu=0.03366;
10918 double nftos=1.95996398454;
10920 if(mA>=65.0 && mA<90.0)
10924 if(mHh>=50.0 && mHh<=1000.0)
10930 else if(mA>=90.0 && mA<100.0)
10942 if(mHh>=50.0 && mHh<=1000.0)
10953 else if(mA>=100.0 && mA<150.0)
10967 if(mHh>=50.0 && mHh<=1000.0)
10978 else if(mA>=150.0 && mA<175.0)
10994 if(mHh>=50.0 && mHh<=1000.0)
11005 else if(mA>=175.0 && mA<200.0)
11021 if(mHh>=50.0 && mHh<=910.0)
11026 if(mHh>=50.0 && mHh<=1000.0)
11037 else if(mA>=200.0 && mA<220.0)
11057 if(mHh>=50.0 && mHh<=910.0)
11062 if(mHh>=50.0 && mHh<=1000.0)
11083 else if(mA>=220.0 && mA<225.0)
11109 if(mHh>=50.0 && mHh<=910.0)
11114 if(mHh>=50.0 && mHh<=1000.0)
11135 else if(mA>=225.0 && mA<250.0)
11163 if(mHh>=50.0 && mHh<=910.0)
11168 if(mHh>=50.0 && mHh<=1000.0)
11189 else if(mA>=250.0 && mA<300.0)
11217 if(mHh>=50.0 && mHh<=910.0)
11222 if(mHh>=50.0 && mHh<=1000.0)
11247 else if(mA>=300.0 && mA<350.0)
11275 if(mHh>=50.0 && mHh<=910.0)
11280 if(mHh>=50.0 && mHh<=1000.0)
11307 else if(mA>=350.0 && mA<400.0)
11333 if(mHh>=50.0 && mHh<=910.0)
11338 if(mHh>=50.0 && mHh<=1000.0)
11367 else if(mA>=400.0 && mA<500.0)
11395 if(mHh>=50.0 && mHh<=910.0)
11400 if(mHh>=50.0 && mHh<=1000.0)
11433 else if(mA>=500.0 && mA<550.0)
11461 if(mHh>=50.0 && mHh<=910.0)
11466 if(mHh>=50.0 && mHh<=1000.0)
11501 else if(mA>=550.0 && mA<600.0)
11529 if(mHh>=50.0 && mHh<=910.0)
11534 if(mHh>=50.0 && mHh<=1000.0)
11571 else if(mA>=600.0 && mA<650.0)
11595 if(mHh>=50.0 && mHh<=910.0)
11600 if(mHh>=50.0 && mHh<=1000.0)
11637 else if(mA>=650.0 && mA<850.0)
11661 if(mHh>=50.0 && mHh<=910.0)
11666 if(mHh>=50.0 && mHh<=1000.0)
11705 else if(mA>=850.0 && mA<900.0)
11727 if(mHh>=50.0 && mHh<=910.0)
11732 if(mHh>=50.0 && mHh<=1000.0)
11771 else if(mA>=900.0 && mA<1000.0)
11791 if(mHh>=50.0 && mHh<=910.0)
11796 if(mHh>=50.0 && mHh<=1000.0)
11835 else if(mA>=1000.0 && mA<1200.0)
11873 else if(mA>=1200.0 && mA<1600.0)
11907 else if(mA>=1600.0 && mA<2000.0)
11939 else if(mA>=2000.0 && mA<2250.0)
11965 else if(mA>=2250.0 && mA<2400.0)
11987 else if(mA>=2400.0 && mA<2500.0)
12007 else if(mA>=2500.0 && mA<2700.0)
12025 else if(mA>=2700.0 && mA<3000.0)
12041 else if(mA>=3000.0 && mA<3200.0)
12052 else if(mA>=3200.0 && mA<4000.0)
12064 double mHh=sqrt(
mHh2);
12065 double mA=sqrt(
mA2);
12066 double mHp=sqrt(
mHp2);
12071 double gHpt=-1.0/
tanb;
12076 double SigmaHp8=0.0;
12077 double SigmaHp13=0.0;
12110 throw std::runtime_error(
"modelflag can be only any of \"type1\", \"type2\", \"typeX\" or \"typeY\"");
12113 double GammaHptaunu=
HSTheta(mHp-
Mtau)*(Mtau2*(
mHp2-Mtau2)*(
mHp2-Mtau2)*gHptau*gHptau)/(8.0*mHp*
mHp2*M_PI*vev2);
12114 double GammaHptb=
HSTheta(mHp-
Mt-
Mb)*(Vtb*Vtb/(8.0*mHp*M_PI*vev2))*3.0*(-4.0*gHpb*gHpt*Mb2*Mt2
12115 -gHpb*gHpb*Mb2*(Mb2-
mHp2+Mt2)-gHpt*gHpt*Mt2*(Mb2-
mHp2+Mt2))
12118 *MW2*MW2/mHp*cos(
bma)*cos(
bma)/(2.0*M_PI*vev2);
12120 *MW2*MW2/mHp*sin(
bma)*sin(
bma)/(2.0*M_PI*vev2);
12122 *MW2*MW2/mHp/(2.0*M_PI*vev2);
12123 GammaHptot= GammaHptaunu + GammaHptb + GammaHpHlW + GammaHpHhW + GammaHpAW;
12125 double Br_Hptotaunu=GammaHptaunu/
GammaHptot;
12159 double nftos=1.95996398454;
12161 if(mHp>=180.0 && mHp<200.0)
12173 else if(mHp>=200.0 && mHp<300.0)
12193 else if(mHp>=300.0 && mHp<600.0)
12215 else if(mHp>=600.0 && mHp<1000.0)
12231 else if(mHp>=1000.0 && mHp<1400.0)
12242 else if(mHp>=1400.0 && mHp<2000.0)
12265 if( RGEorder ==
"LO" ) flag=0;
12266 else if( RGEorder ==
"approxNLO" ) flag=1;
12267 else if( RGEorder ==
"NLO" ) flag=2;
12269 throw std::runtime_error(
"RGEorder can be only any of \"LO\", \"approxNLO\" or \"NLO\"");
12272 double g1_at_MZ=sqrt(4.0*M_PI*
Ale/
cW2);
12273 double g2_at_MZ=sqrt(4.0*M_PI*
Ale/(1-
cW2));
12274 double g3_at_MZ=sqrt(4.0*M_PI*
Als);
12277 double Ybottom1_at_MZ=0.0;
12278 double Ybottom2_at_MZ=0.0;
12279 double Ytau1_at_MZ=0.0;
12280 double Ytau2_at_MZ=0.0;
12283 double Mb_at_MZ=2.96;
12284 double Mtau_at_MZ=1.75;
12287 Ybottom2_at_MZ=(sqrt(2.0)*Mb_at_MZ)/(
vev*sinb);
12288 Ytau2_at_MZ=(sqrt(2.0)*Mtau_at_MZ)/(
vev*sinb);
12291 Ybottom1_at_MZ=(sqrt(2.0)*Mb_at_MZ)/(
vev*cosb);
12292 Ytau1_at_MZ=(sqrt(2.0)*Mtau_at_MZ)/(
vev*cosb);
12296 Ybottom2_at_MZ=(sqrt(2.0)*Mb_at_MZ)/(
vev*sinb);
12297 Ytau1_at_MZ=(sqrt(2.0)*Mtau_at_MZ)/(
vev*cosb);
12301 Ybottom1_at_MZ=(sqrt(2.0)*Mb_at_MZ)/(
vev*cosb);
12302 Ytau2_at_MZ=(sqrt(2.0)*Mtau_at_MZ)/(
vev*sinb);
12306 throw std::runtime_error(
"modelflag can be only any of \"type1\", \"type2\", \"typeX\" or \"typeY\"");
12344 double InitVals[14];
12345 InitVals[0]=g1_at_MZ;
12346 InitVals[1]=g2_at_MZ;
12347 InitVals[2]=g3_at_MZ;
12348 InitVals[3]=Ytop_at_MZ;
12349 InitVals[4]=Ybottom1_at_MZ+Ybottom2_at_MZ;
12350 InitVals[5]=Ytau1_at_MZ+Ytau2_at_MZ;
12351 InitVals[6]=m11_2_at_MZ;
12352 InitVals[7]=m22_2_at_MZ;
12353 InitVals[8]=m12_2_at_MZ;
12354 InitVals[9]=lambda1_at_MZ;
12355 InitVals[10]=lambda2_at_MZ;
12356 InitVals[11]=lambda3_at_MZ;
12357 InitVals[12]=lambda4_at_MZ;
12358 InitVals[13]=lambda5_at_MZ;
12387 throw std::runtime_error(
"modelflag can be only any of \"type1\", \"type2\", \"typeX\" or \"typeY\"");
12401 std::cout <<
"\n \t Warning: UnitarityTestflag is switched on.\n \t All lambda's are hard coded in this mode!\n" << std::endl;
12424 double WFRcomb1a = 0.0;
12425 double WFRcomb1b = 0.0;
12426 double WFRcomb2a = 0.0;
12427 double WFRcomb3a = 0.0;
12428 double WFRcomb3b = 0.0;
12429 double WFRcomb4a = 0.0;
12435 double beta=atan(
tanb);
12436 double alpha=beta-
bma;
12504 WFRcomb1a = 3.0*
mHl2*
mHl2*cosb*cosb*sin(
bma)*sin(
bma) * ddpB000mh
12511 + 1.5*
mHl2*
mHl2*sin(alpha)*sin(alpha)*sin(
bma)*sin(
bma) * ddpB0mh00
12514 + 9.0*sin(alpha)*sin(alpha)*pow(-
mHl2*(3.0*sin(
bma)+sin(3.0*
bma)+sin(3.0*alpha+beta)+3.0*sin(alpha+3.0*beta))
12515 +16.0*
m12_2*cos(
bma)*cos(
bma)*cos(alpha+beta),2)/(512.0*pow(cosb*sinb,4)) * ddpB0mhmhmh
12516 + sin(alpha)*sin(alpha)*pow((cos(alpha)/sinb+sin(alpha)/cosb)*(
m12_2+cos(alpha)*sin(alpha)*(
mHh2+2.0*
mHl2-(3.0*
m12_2)/(cosb*sinb))),2) * ddpB0mhmHmh
12517 + sin(alpha)*sin(alpha)*sin(
bma)*sin(
bma)*pow(-2.0*
m12_2+(2.0*
mHh2+
mHl2-(3.0*
m12_2)/(cosb*sinb))*sin(2.0*alpha),2)/(8.0*cosb*cosb*sinb*sinb) * ddpB0mhmHmH
12518 + (sin(alpha)*sin(alpha)*pow((
mHl2-2.0*
mHp2)*cos(alpha-3.0*beta)*2.0*sinb*cosb
12519 +cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*2.0*sinb*cosb),2))/(64.0*pow(cosb*sinb,4)) * ddpB0mhmHpmHp
12520 + (sin(alpha)*sin(alpha)*pow((2.0*
mA2-
mHl2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(8.0*
m12_2-(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb),2))/(128.0*pow(cosb*sinb,4)) * ddpB0mhmAmA
12521 + 1.5*
mHh2*
mHh2*cos(alpha)*cos(alpha)*cos(
bma)*cos(
bma) * ddpB0mH00
12524 + cos(alpha)*cos(alpha)*cos(
bma)*cos(
bma)*pow(
m12_2+cos(alpha)*sin(alpha)*(
mHh2-3.0*
m12_2/(cosb*sinb))+
mHl2*sin(2.0*alpha),2)/(2.0*cosb*cosb*sinb*sinb) * ddpB0mHmhmh
12525 + cos(alpha)*cos(alpha)*sin(
bma)*sin(
bma)*pow(
m12_2*cosb*sinb+0.5*sin(2.0*alpha)*(3.0*
m12_2-(2.0*
mHh2+
mHl2)*cosb*sinb),2)/pow(cosb*sinb,4) * ddpB0mHmHmh
12526 + 9.0*cos(alpha)*cos(alpha)*pow(
mHh2*(-3.0*cos(
bma)+cos(3.0*
bma)-cos(3.0*alpha+beta)+3.0*cos(alpha+3.0*beta))
12527 +16.0*
m12_2*sin(
bma)*sin(
bma)*sin(alpha+beta),2)/(512.0*pow(cosb*sinb,4)) * ddpB0mHmHmH
12528 + (cos(alpha)*cos(alpha)*pow((
mHh2-2.0*
mHp2)*cos(alpha-5.0*beta)
12529 +2.0*(
mHh2+2.0*
mHp2)*cos(
bma)-(3.0*
mHh2+2.0*
mHp2)*cos(alpha+3.0*beta)-16.0*
m12_2*sin(alpha+beta),2))/(256.0*pow(cosb*sinb,4)) * ddpB0mHmHpmHp
12530 + (cos(alpha)*cos(alpha)*pow((2.0*
mA2-
mHh2)*cos(alpha-5.0*beta)
12531 -2.0*(2.0*
mA2+
mHh2)*cos(
bma)+(2.0*
mA2+3.0*
mHh2)*cos(alpha+3.0*beta)+16.0*
m12_2*sin(alpha+beta),2))/(512.0*pow(cosb*sinb,4)) * ddpB0mHmAmA
12535 + 2.0*pow((
m12_2*cos(alpha+beta))/(sinb*sinb*cosb*cosb)-(
mHl2*cos(
bma)*cos(2.0*beta))/(cosb*sinb)-(
mHl2+2.0*
mHp2)*sin(
bma),2)*sinb*sinb * ddpB0mHpmHpmh
12536 + 2.0*pow(sinb*((
mHh2+2.0*
mHp2)*cos(
bma)-
mHh2*cos(2.0*beta)*sin(
bma)/(cosb*sinb)-
m12_2*sin(alpha+beta)/(sinb*sinb*cosb*cosb)),2) * ddpB0mHpmHpmH
12540 + pow((2.0*
mA2-
mHl2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(8.0*
m12_2-(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb),2)/(64.0*pow(cosb,4)*sinb*sinb) * ddpB0mAmAmh
12541 + pow((2.0*
mA2-
mHh2)*cos(alpha-5.0*beta)-2.0*(2.0*
mA2+
mHh2)*cos(
bma)+2.0*
mA2*cos(alpha+3.0*beta)
12542 + 3.0*
mHh2*cos(alpha+3.0*beta)+16.0*
m12_2*sin(alpha+beta),2)/(256.0*pow(cosb,4)*sinb*sinb) * ddpB0mAmAmH;
12547 +cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*2.0*sinb*cosb)))/(2.0*
mHp2*sinb*cosb) * B00mHpmh
12549 -2.0*
mHp2*cos(alpha+3.0*beta)-16.0*
m12_2*sin(alpha+beta)))/(4.0*
mHp2*sinb*cosb) * B00mHpmH
12551 +cos(alpha+beta)*(-8.0*
m12_2+(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb)))/(4.0*
mA2*sinb*cosb) * B00mAmh
12553 +3.0*
mHh2*cos(alpha+3.0*beta)+16.0*
m12_2*sin(alpha+beta)))/(8.0*
mA2*sinb*cosb) * B00mAmH
12557 + 3.0*cos(
bma)*sin(2.0*alpha)*(
m12_2+cos(alpha)*sin(alpha)*(
mHh2-(3.0*
m12_2)/(sinb*cosb))+
mHl2*sin(2.0*alpha))
12558 *(-
mHl2*(3.0*sin(
bma)+sin(3.0*
bma)+sin(3.0*alpha+beta)+3.0*sin(alpha+3.0*beta))
12559 +16.0*
m12_2*cos(
bma)*cos(
bma)*cos(alpha+beta))/(32.0*(
mHl2-
mHh2)*pow(sinb*cosb,3)) * B0mhmhmh
12560 + (sin(2.0*
bma)*sin(2.0*alpha)
12563 +(9.0*
m12_2*
m12_2)/(sinb*sinb*cosb*cosb))*sin(2.0*alpha)*sin(2.0*alpha)))/(8.0*(
mHh2-
mHl2)*sinb*sinb*cosb*cosb) * B0mhmHmh
12564 - 3.0*sin(2.0*alpha)*sin(
bma)*(
m12_2-cos(alpha)*sin(alpha)*(2.0*
mHh2+
mHl2-(3.0*
m12_2)/(sinb*cosb)))
12565 *(
mHh2*(-3.0*cos(
bma)+cos(3.0*
bma)-cos(3.0*alpha+beta)+3.0*cos(alpha+3.0*beta))
12566 +16.0*
m12_2*sin(
bma)*sin(
bma)*sin(alpha+beta))/(32.0*(
mHh2-
mHl2)*pow(sinb*cosb,3)) * B0mhmHmH
12567 + sin(2.0*alpha)*((
mHl2-2.0*
mHp2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*2.0*sinb*cosb))
12569 -2.0*
mHp2*cos(alpha+3.0*beta)-16.0*
m12_2*sin(alpha+beta))/(128.0*(
mHh2-
mHl2)*pow(cosb*sinb,4)) * B0mhmHpmHp
12570 + sin(2.0*alpha)*((
mHl2-2.0*
mA2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(-8.0*
m12_2+(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb))
12572 -3.0*
mHh2*cos(alpha+3.0*beta)-16.0*
m12_2*sin(alpha+beta))/(256.0*(
mHh2-
mHl2)*pow(cosb*sinb,4)) * B0mhmAmA
12576 + 3.0*cos(
bma)*sin(2.0*alpha)*(
m12_2+cos(alpha)*sin(alpha)*(
mHh2-(3.0*
m12_2)/(sinb*cosb))+
mHl2*sin(2.0*alpha))
12577 *(-
mHl2*(3.0*sin(
bma)+sin(3.0*
bma)+sin(3.0*alpha+beta)+3.0*sin(alpha+3.0*beta))
12578 +16.0*
m12_2*cos(
bma)*cos(
bma)*cos(alpha+beta))/(32.0*(
mHh2-
mHl2)*pow(sinb*cosb,3)) * B0mHmhmh
12582 +(9.0*
m12_2*
m12_2)/(sinb*sinb*cosb*cosb))*sin(2.0*alpha)*sin(2.0*alpha)))
12583 /(8.0*(
mHh2-
mHl2)*sinb*sinb*cosb*cosb) * B0mHmHmh
12584 + 3.0*sin(
bma)*sin(2.0*alpha)*(
m12_2-cos(alpha)*sin(alpha)*(2.0*
mHh2+
mHl2-(3.0*
m12_2)/(sinb*cosb)))
12585 *(
mHh2*(-3.0*cos(
bma)+cos(3.0*
bma)-cos(3.0*alpha+beta)+3.0*cos(alpha+3.0*beta))
12586 +16.0*
m12_2*sin(
bma)*sin(
bma)*sin(alpha+beta))/(32.0*(
mHh2-
mHl2)*pow(sinb*cosb,3)) * B0mHmHmH
12587 - sin(2.0*alpha)*((
mHl2-2.0*
mHp2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*2.0*sinb*cosb))
12589 -2.0*
mHp2*cos(alpha+3.0*beta)-16.0*
m12_2*sin(alpha+beta))/(128.0*(
mHh2-
mHl2)*pow(cosb*sinb,4)) * B0mHmHpmHp
12590 - sin(2.0*alpha)*((
mHl2-2.0*
mA2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mA2)*2.0*sinb*cosb))
12592 -2.0*
mA2*cos(alpha+3.0*beta)-16.0*
m12_2*sin(alpha+beta))/(256.0*(
mHh2-
mHl2)*pow(cosb*sinb,4)) * B0mHmAmA
12596 +cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*2.0*sinb*cosb)))/(2.0*
mHp2*sinb*cosb) * B0mHpmHpmh
12598 -2.0*
mHp2*cos(alpha+3.0*beta)-16.0*
m12_2*sin(alpha+beta)))/(4.0*
mHp2*sinb*cosb) * B0mHpmHpmH
12602 +cos(alpha+beta)*(-8.0*
m12_2+(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb)))/(4.0*
mA2*sinb*cosb) * B0mAmAmh
12604 +3.0*
mHh2*cos(alpha+3.0*beta)+16.0*
m12_2*sin(alpha+beta)))/(8.0*
mA2*sinb*cosb) * B0mAmAmH;
12616 + 9.0*pow(16.0*
m12_2*cos(
bma)*cos(
bma)*cos(alpha+beta)
12617 -
mHl2*(3.0*sin(
bma)+sin(3.0*
bma)+sin(3.0*alpha+beta)+3.0*sin(alpha+3.0*beta)),2)/(1024.0*pow(cosb*sinb,4)) * ddpB0mhmhmh
12618 + 0.5*pow(cos(alpha)/sinb + sin(alpha)/cosb,2)
12619 *pow(
m12_2+cos(alpha)*sin(alpha)*(
mHh2+2.0*
mHl2-3.0*
m12_2/(cosb*sinb)),2) * ddpB0mhmHmh
12621 *sin(
bma)*sin(
bma))/(16.0*cosb*cosb*sinb*sinb) * ddpB0mhmHmH
12622 + pow((
mHl2-2.0*
mHp2)*cos(alpha-3.0*beta)*cosb*sinb
12623 +cos(alpha+beta)*(-4.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*cosb*sinb),2)/(32.0*pow(cosb*sinb,4)) * ddpB0mhmHpmHp
12624 + pow((2.0*
mA2-
mHl2)*cos(alpha-3.0*beta)*cosb*sinb
12625 +cos(alpha+beta)*(4.0*
m12_2-(2.0*
mA2+3.0*
mHl2)*cosb*sinb),2)/(64.0*pow(cosb*sinb,4)) * ddpB0mhmAmA
12629 + 0.25*pow(cos(alpha)/sinb + sin(alpha)/cosb,2)
12630 *pow(
m12_2+cos(alpha)*sin(alpha)*(
mHh2+2.0*
mHl2-3.0*
m12_2/(cosb*sinb)),2) * ddpB0mHmhmh
12632 *sin(
bma)*sin(
bma))/(8.0*cosb*cosb*sinb*sinb) * ddpB0mHmHmh
12633 + 9.0*pow(
mHh2*(-3.0*cos(
bma)+cos(3.0*
bma)-cos(3.0*alpha+beta)+3.0*cos(alpha+3.0*beta))
12634 +16.0*
m12_2*sin(
bma)*sin(
bma)*sin(alpha+beta),2)/(1024.0*pow(cosb*sinb,4)) * ddpB0mHmHmH
12636 -16.0*
m12_2*sin(alpha+beta),2)/(512*pow(cosb*sinb,4)) * ddpB0mHmHpmHp
12638 +16.0*
m12_2*sin(alpha+beta),2)/(1024*pow(cosb*sinb,4)) * ddpB0mHmAmA
12642 + pow((
m12_2*cos(alpha+beta))/(cosb*cosb*sinb*sinb)-(
mHl2*cos(
bma)*cos(2.0*beta))/(cosb*sinb)
12645 -(
m12_2*sin(alpha+beta))/(cosb*cosb*sinb*sinb),2) * ddpB0mHpmHpmH
12649 + pow((2.0*
mA2-
mHl2)*cos(alpha-3.0*beta)*2.0*cosb*sinb
12650 +cos(alpha+beta)*(8.0*
m12_2-(2.0*
mA2+3.0*
mHl2)*2.0*cosb*sinb),2)/(128.0*pow(cosb*sinb,4)) * ddpB0mAmAmh
12652 +16.0*
m12_2*sin(alpha+beta),2)/(512.0*pow(cosb*sinb,4)) * ddpB0mAmAmH;
12654 WFRcomb3a = 0.5*
mHl2*
mHl2*(3.0-cos(2.0*beta))*sin(
bma)*sin(
bma) * ddpB000mh
12655 + 0.5*
mHh2*
mHh2*(3.0-cos(2.0*beta))*cos(
bma)*cos(
bma) * ddpB000mH
12661 + 1.5*
mHl2*
mHl2*sin(alpha)*sin(alpha)*sin(
bma)*sin(
bma) * ddpB0mh00
12664 + 9.0*sin(alpha)*sin(alpha)*pow(-
mHl2*(3.0*sin(
bma)+sin(3.0*
bma)+sin(3.0*alpha+beta)+3.0*sin(alpha+3.0*beta))
12665 +16.0*
m12_2*cos(
bma)*cos(
bma)*cos(alpha+beta),2)/(512.0*pow(cosb*sinb,4)) * ddpB0mhmhmh
12666 + sin(alpha)*sin(alpha)*pow((cos(alpha)/sinb+sin(alpha)/cosb)*(
m12_2+cos(alpha)*sin(alpha)*(
mHh2+2.0*
mHl2-(3.0*
m12_2)/(cosb*sinb))),2) * ddpB0mhmHmh
12667 + sin(alpha)*sin(alpha)*sin(
bma)*sin(
bma)*pow(-2.0*
m12_2+(2.0*
mHh2+
mHl2-(3.0*
m12_2)/(cosb*sinb))*sin(2.0*alpha),2)/(8.0*cosb*cosb*sinb*sinb) * ddpB0mhmHmH
12668 + (sin(alpha)*sin(alpha)*pow((
mHl2-2.0*
mHp2)*cos(alpha-3.0*beta)*2.0*sinb*cosb
12669 +cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*2.0*sinb*cosb),2))/(64.0*pow(cosb*sinb,4)) * ddpB0mhmHpmHp
12670 + (sin(alpha)*sin(alpha)*pow((2.0*
mA2-
mHl2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(8.0*
m12_2-(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb),2))/(128.0*pow(cosb*sinb,4)) * ddpB0mhmAmA
12671 + 1.5*
mHh2*
mHh2*cos(alpha)*cos(alpha)*cos(
bma)*cos(
bma) * ddpB0mH00
12674 + cos(alpha)*cos(alpha)*cos(
bma)*cos(
bma)*pow(
m12_2+cos(alpha)*sin(alpha)*(
mHh2-3.0*
m12_2/(cosb*sinb))+
mHl2*sin(2.0*alpha),2)/(2.0*cosb*cosb*sinb*sinb) * ddpB0mHmhmh
12675 + cos(alpha)*cos(alpha)*sin(
bma)*sin(
bma)*pow(
m12_2*cosb*sinb+0.5*sin(2.0*alpha)*(3.0*
m12_2-(2.0*
mHh2+
mHl2)*cosb*sinb),2)/pow(cosb*sinb,4) * ddpB0mHmHmh
12676 + 9.0*cos(alpha)*cos(alpha)*pow(
mHh2*(-3.0*cos(
bma)+cos(3.0*
bma)-cos(3.0*alpha+beta)+3.0*cos(alpha+3.0*beta))
12677 +16.0*
m12_2*sin(
bma)*sin(
bma)*sin(alpha+beta),2)/(512.0*pow(cosb*sinb,4)) * ddpB0mHmHmH
12678 + (cos(alpha)*cos(alpha)*pow((
mHh2-2.0*
mHp2)*cos(alpha-5.0*beta)
12679 +2.0*(
mHh2+2.0*
mHp2)*cos(
bma)-(3.0*
mHh2+2.0*
mHp2)*cos(alpha+3.0*beta)-16.0*
m12_2*sin(alpha+beta),2))/(256.0*pow(cosb*sinb,4)) * ddpB0mHmHpmHp
12680 + (cos(alpha)*cos(alpha)*pow((2.0*
mA2-
mHh2)*cos(alpha-5.0*beta)
12681 -2.0*(2.0*
mA2+
mHh2)*cos(
bma)+(2.0*
mA2+3.0*
mHh2)*cos(alpha+3.0*beta)+16.0*
m12_2*sin(alpha+beta),2))/(512.0*pow(cosb*sinb,4)) * ddpB0mHmAmA
12685 + 2.0*cosb*cosb*pow((
m12_2*cos(alpha+beta))/(cosb*cosb*sinb*sinb)-(
mHl2*cos(
bma)*cos(2.0*beta))/(cosb*sinb)-(
mHl2+2.0*
mHp2)*sin(
bma),2) * ddpB0mHpmHpmh
12686 + 2.0*cosb*cosb*pow((
mHh2+2.0*
mHp2)*cos(
bma)-(
mHh2*cos(2.0*beta)*sin(
bma))/(cosb*sinb)-(
m12_2*sin(alpha+beta))/(cosb*cosb*sinb*sinb),2) * ddpB0mHpmHpmH
12690 + pow((2.0*
mA2-
mHl2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(8.0*
m12_2-(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb),2)/(64.0*pow(cosb,4)*sinb*sinb) * ddpB0mAmAmh
12691 + pow((2.0*
mA2-
mHh2)*cos(alpha-5.0*beta)-2.0*(2.0*
mA2+
mHh2)*cos(
bma)+2.0*
mA2*cos(alpha+3.0*beta)
12692 + 3.0*
mHh2*cos(alpha+3.0*beta)+16.0*
m12_2*sin(alpha+beta),2)/(256.0*pow(cosb,4)*sinb*sinb) * ddpB0mAmAmH;
12697 +cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*2.0*sinb*cosb)))/(2.0*
mHp2*sinb*cosb) * B00mHpmh
12699 -2.0*
mHp2*cos(alpha+3.0*beta)-16.0*
m12_2*sin(alpha+beta)))/(4.0*
mHp2*sinb*cosb) * B00mHpmH
12701 +cos(alpha+beta)*(-8.0*
m12_2+(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb)))/(4.0*
mA2*sinb*cosb) * B00mAmh
12703 +3.0*
mHh2*cos(alpha+3.0*beta)+16.0*
m12_2*sin(alpha+beta)))/(8.0*
mA2*sinb*cosb) * B00mAmH
12707 + 3.0*cos(
bma)*sin(2.0*alpha)*(
m12_2+cos(alpha)*sin(alpha)*(
mHh2-(3.0*
m12_2)/(sinb*cosb))+
mHl2*sin(2.0*alpha))
12708 *(-
mHl2*(3.0*sin(
bma)+sin(3.0*
bma)+sin(3.0*alpha+beta)+3.0*sin(alpha+3.0*beta))
12709 +16.0*
m12_2*cos(
bma)*cos(
bma)*cos(alpha+beta))/(32.0*(
mHl2-
mHh2)*pow(sinb*cosb,3)) * B0mhmhmh
12710 + (sin(2.0*
bma)*sin(2.0*alpha)
12713 +(9.0*
m12_2*
m12_2)/(sinb*sinb*cosb*cosb))*sin(2.0*alpha)*sin(2.0*alpha)))/(8.0*(
mHh2-
mHl2)*sinb*sinb*cosb*cosb) * B0mhmHmh
12714 - 3.0*sin(2.0*alpha)*sin(
bma)*(
m12_2-cos(alpha)*sin(alpha)*(2.0*
mHh2+
mHl2-(3.0*
m12_2)/(sinb*cosb)))
12715 *(
mHh2*(-3.0*cos(
bma)+cos(3.0*
bma)-cos(3.0*alpha+beta)+3.0*cos(alpha+3.0*beta))
12716 +16.0*
m12_2*sin(
bma)*sin(
bma)*sin(alpha+beta))/(32.0*(
mHh2-
mHl2)*pow(sinb*cosb,3)) * B0mhmHmH
12717 + sin(2.0*alpha)*((
mHl2-2.0*
mHp2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*2.0*sinb*cosb))
12719 -2.0*
mHp2*cos(alpha+3.0*beta)-16.0*
m12_2*sin(alpha+beta))/(128.0*(
mHh2-
mHl2)*pow(cosb*sinb,4)) * B0mhmHpmHp
12720 + sin(2.0*alpha)*((
mHl2-2.0*
mA2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(-8.0*
m12_2+(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb))
12722 -3.0*
mHh2*cos(alpha+3.0*beta)-16.0*
m12_2*sin(alpha+beta))/(256.0*(
mHh2-
mHl2)*pow(cosb*sinb,4)) * B0mhmAmA
12726 + 3.0*cos(
bma)*sin(2.0*alpha)*(
m12_2+cos(alpha)*sin(alpha)*(
mHh2-(3.0*
m12_2)/(sinb*cosb))+
mHl2*sin(2.0*alpha))
12727 *(-
mHl2*(3.0*sin(
bma)+sin(3.0*
bma)+sin(3.0*alpha+beta)+3.0*sin(alpha+3.0*beta))
12728 +16.0*
m12_2*cos(
bma)*cos(
bma)*cos(alpha+beta))/(32.0*(
mHh2-
mHl2)*pow(sinb*cosb,3)) * B0mHmhmh
12732 +(9.0*
m12_2*
m12_2)/(sinb*sinb*cosb*cosb))*sin(2.0*alpha)*sin(2.0*alpha)))
12733 /(8.0*(
mHh2-
mHl2)*sinb*sinb*cosb*cosb) * B0mHmHmh
12734 + 3.0*sin(
bma)*sin(2.0*alpha)*(
m12_2-cos(alpha)*sin(alpha)*(2.0*
mHh2+
mHl2-(3.0*
m12_2)/(sinb*cosb)))
12735 *(
mHh2*(-3.0*cos(
bma)+cos(3.0*
bma)-cos(3.0*alpha+beta)+3.0*cos(alpha+3.0*beta))
12736 +16.0*
m12_2*sin(
bma)*sin(
bma)*sin(alpha+beta))/(32.0*(
mHh2-
mHl2)*pow(sinb*cosb,3)) * B0mHmHmH
12737 - sin(2.0*alpha)*((
mHl2-2.0*
mHp2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*2.0*sinb*cosb))
12739 -2.0*
mHp2*cos(alpha+3.0*beta)-16.0*
m12_2*sin(alpha+beta))/(128.0*(
mHh2-
mHl2)*pow(cosb*sinb,4)) * B0mHmHpmHp
12740 - sin(2.0*alpha)*((
mHl2-2.0*
mA2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mA2)*2.0*sinb*cosb))
12742 -2.0*
mA2*cos(alpha+3.0*beta)-16.0*
m12_2*sin(alpha+beta))/(256.0*(
mHh2-
mHl2)*pow(cosb*sinb,4)) * B0mHmAmA
12746 +cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*2.0*sinb*cosb))/(2.0*
mHp2*sinb*cosb) * B0mHpmHpmh
12748 -3.0*
mHh2*cos(alpha+3.0*beta)-2.0*
mHp2*cos(alpha+3.0*beta)-16.0*
m12_2*sin(alpha+beta))/(4.0*
mHp2*sinb*cosb) * B0mHpmHpmH
12752 +cos(alpha+beta)*(-8.0*
m12_2+(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb)))/(4.0*
mA2*sinb*cosb) * B0mAmAmh
12754 +3.0*
mHh2*cos(alpha+3.0*beta)+16.0*
m12_2*sin(alpha+beta)))/(8.0*
mA2*sinb*cosb) * B0mAmAmH;
12763 + 9.0*pow(16.0*
m12_2*cos(
bma)*cos(
bma)*cos(alpha+beta)
12764 -
mHl2*(3.0*sin(
bma)+sin(3.0*
bma)+sin(3.0*alpha+beta)+3.0*sin(alpha+3.0*beta)),2)/(1024.0*pow(cosb*sinb,4)) * ddpB0mhmhmh
12765 + 0.5*pow(cos(alpha)/sinb + sin(alpha)/cosb,2)
12766 *pow(
m12_2+cos(alpha)*sin(alpha)*(
mHh2+2.0*
mHl2-3.0*
m12_2/(cosb*sinb)),2) * ddpB0mhmHmh
12768 *sin(
bma)*sin(
bma))/(16.0*cosb*cosb*sinb*sinb) * ddpB0mhmHmH
12769 + pow((
mHl2-2.0*
mHp2)*cos(alpha-3.0*beta)*cosb*sinb
12770 +cos(alpha+beta)*(-4.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*cosb*sinb),2)/(32.0*pow(cosb*sinb,4)) * ddpB0mhmHpmHp
12771 + pow((2.0*
mA2-
mHl2)*cos(alpha-3.0*beta)*cosb*sinb
12772 +cos(alpha+beta)*(4.0*
m12_2-(2.0*
mA2+3.0*
mHl2)*cosb*sinb),2)/(64.0*pow(cosb*sinb,4)) * ddpB0mhmAmA
12776 + 0.25*pow(cos(alpha)/sinb + sin(alpha)/cosb,2)
12777 *pow(
m12_2+cos(alpha)*sin(alpha)*(
mHh2+2.0*
mHl2-3.0*
m12_2/(cosb*sinb)),2) * ddpB0mHmhmh
12779 *sin(
bma)*sin(
bma))/(8.0*cosb*cosb*sinb*sinb) * ddpB0mHmHmh
12780 + 9.0*pow(
mHh2*(-3.0*cos(
bma)+cos(3.0*
bma)-cos(3.0*alpha+beta)+3.0*cos(alpha+3.0*beta))
12781 +16.0*
m12_2*sin(
bma)*sin(
bma)*sin(alpha+beta),2)/(1024.0*pow(cosb*sinb,4)) * ddpB0mHmHmH
12783 -16.0*
m12_2*sin(alpha+beta),2)/(512*pow(cosb*sinb,4)) * ddpB0mHmHpmHp
12785 +16.0*
m12_2*sin(alpha+beta),2)/(1024*pow(cosb*sinb,4)) * ddpB0mHmAmA
12789 + pow((2.0*
mA2-
mHl2)*cos(alpha-3.0*beta)*2.0*cosb*sinb
12790 +cos(alpha+beta)*(8.0*
m12_2-(2.0*
mA2+3.0*
mHl2)*2.0*cosb*sinb),2)/(128.0*pow(cosb*sinb,4)) * ddpB0mAmAmh
12792 +16.0*
m12_2*sin(alpha+beta),2)/(512.0*pow(cosb*sinb,4)) * ddpB0mAmAmH;
const gslpp::complex getV_tb() const
A member for returning the value of the CKM element .
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
.
const double & getMass() const
A get method to access the particle mass.
const Particle & getQuarks(const QCD::quark q) const
A get method to access a quark as an object of the type Particle.
An RGE running algorithm.
virtual double RGERunner(double InitialValues[], unsigned long int NumberOfRGEs, double Q1, double Q2, int order, double Rpeps, double tNLOuni)
A model class for the Standard Model.
const double computeGammaHTotal() const
The Higgs total width in the Standard Model.
const Particle & getLeptons(const QCD::lepton p) const
A get method to retrieve the member object of a lepton.
const double computeSigmattH(const double sqrt_s) const
The ttH production cross section in the Standard Model.
const double computeSigmaggH(const double sqrt_s) const
The ggH cross section in the Standard Model.
const double getMz() const
A get method to access the mass of the boson .
const double computeSigmaVBF(const double sqrt_s) const
The VBF cross section in the Standard Model.
const CKM & getCKM() const
A get method to retrieve the member object of type CKM.
const double getAlsMz() const
A get method to access the value of .
const double computeSigmaWH(const double sqrt_s) const
The WH production cross section in the Standard Model.
virtual const double getMHl() const
A get method to retrieve the Higgs mass .
const double computeSigmaggH_tt(const double sqrt_s) const
The square of the top-quark contribution to the ggH cross section in the Standard Model.
const double computeSigmaggH_bb(const double sqrt_s) const
The square of the bottom-quark contribution to the ggH cross section in the Standard Model.
const double computeSigmaZH(const double sqrt_s) const
The ZH production cross section in the Standard Model.
const double c02() const
The square of the cosine of the weak mixing angle defined without weak radiative corrections.
const double Mw_tree() const
The tree-level mass of the boson, .
const double v() const
The Higgs vacuum expectation value.
const double getAle() const
A get method to retrieve the fine-structure constant .
THDM branching ratio of .
THDM branching ratio of .
THDM branching ratio of .
A base class for symmetric Two-Higgs-Doublet models.
double getBDtaunu_SM() const
A method get the SM expectation.
std::string getModelTypeflag() const
A method get the THDM model type.
double gettanb() const
A method get .
double getBDstartaunu_B() const
A method get the coefficient B.
double getNLOuniscale() const
A method get the minimal NLO unitarity check scale.
double getBDstartaunu_A() const
A method get the coefficient A.
double getBDtaunu_A() const
A method get the coefficient A.
double getBDtaunu_B() const
A method get the coefficient B.
double getbma() const
A method get .
double getQ_THDM() const
A method get the THDM scale.
double getmHh2() const
A method get the squared mass of the "non-125 GeV" neutral scalar Higgs.
double getsinb() const
A method get .
double getmHp2() const
A method get the squared charged Higgs mass.
double getRpeps() const
A method get the minimal R' value.
std::string getRGEorderflag() const
A method get switch for NLO RGE and approximate NLO RGE.
double getm12_2() const
A method get the parameter of the Higgs potential .
double getcosb() const
A method get .
double getlogtb() const
A method get .
bool getWFRflag() const
A method get the flag to switch on wavefunction renormalization for the NLO unitarity conditions.
double getmA2() const
A method get the squared mass of the pseudoscalar Higgs A.
double getBDstartaunu_SM() const
A method get the SM expectation.
gslpp::matrix< double > CMS13_ggF_phi_Zga_e
double ip_ex_pp_phi_bb_CMS13_cache_e[2][CacheSize]
gslpp::complex B0_MZ2_mHl2_mHl2_mHl2_cache[3][CacheSize]
double ip_ex_pp_A_HZ_tautaull_CMS8(double mA, double mH)
Interpolating function for the observed CMS upper limit on a pseudoscalar resonance decaying to whic...
gslpp::matrix< double > CMS13_ggF_H_hh_bbbb_e
double ip_ex_bb_phi_tautau_CMS8_cache_e[2][CacheSize]
gslpp::matrix< double > CMS13_pp_H_ZZ_llnunu
gslpp::matrix< double > CMS13_VBFVH_H_ZZ_llll
double ip_ex_pp_A_HZ_tautaull_CMS8_cache[3][CacheSize]
gslpp::matrix< double > temp1e
double ip_ex_VBF_H_ZZ_llllnunu_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a VBF produced scalar resonance decaying...
gslpp::complex B0p_MZ2_mHl2_mHl2_mHl2_cache[3][CacheSize]
double ip_cs_pptottA_8(double mass)
Interpolating function for the top associated A production cross section at 8 TeV.
gslpp::matrix< double > ATLAS13_pp_phi_Zga_e
gslpp::complex B0_MZ2_mHl2_mHh2_mHh2_cache[4][CacheSize]
gslpp::complex B00_MZ2_0_mHl2_mHp2(const double MZ2, const double mHl2, const double mHp2) const
.
gslpp::matrix< double > temp40
gslpp::matrix< double > temp9e
gslpp::matrix< double > log_cs_ggH_8
double ie16e(double mass)
gslpp::complex B00_MZ2_0_mA2_mHp2_cache[4][CacheSize]
double ip_ex_gg_A_hZ_tautaull_CMS8_cache_e[2][CacheSize]
gslpp::matrix< double > ATLAS13_bb_A_Zh_Zbb_e
gslpp::matrix< double > temp8
gslpp::matrix< double > log_cs_WH_13
gslpp::complex B0p_MZ2_mHp2_0_mA2_cache[4][CacheSize]
double ip_ex_pp_phi_hh_gagabb_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double R_pp_A_gaga_ATLAS8
gslpp::matrix< double > temp18e
gslpp::complex B0_MZ2_mHl2_0_mA2_cache[4][CacheSize]
double ie8c[2][CacheSize]
double ip_ex_VBF_H_ZZ_qqllnunu_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a VBF produced scalar resonance decaying...
gslpp::complex A_HH_L_cache[6][CacheSize]
gslpp::matrix< double > csrA_bottom_13
gslpp::complex B0_MZ2_mHl2_0_0_cache[3][CacheSize]
gslpp::complex B00_MZ2_MZ2_MZ2_mHl2(const double MZ2, const double mHl2) const
.
double ip_ex_VBF_H_ZZ_llqq_ATLAS13_cache_e[2][CacheSize]
double ip_ex_bb_phi_tautau_CMS8_cache[2][CacheSize]
void computeHpquantities()
gslpp::complex I_A_U(const double mA2, const double Mc, const double Mt) const
Amplitude for a CP-odd Higgs boson decay to diphotons including the charm and top quarks in the loop.
gslpp::matrix< double > ATLAS13_VBF_H_ZZ_llll
double ip_ex_bb_phi_tt_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a bb associated scalar resonance decayin...
double THoEX_ggF_A_hZ_bbll_CMS8
double R_pp_H_hh_bbgaga_ATLAS13
double ip_ex_ggF_phi_Zga_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced scalar resonance de...
double KaellenFunction(const double a2, const double b2, const double c2) const
Kaellen function.
double R_pp_H_Zga_ATLAS13
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 THoEX_ggF_H_ZZ_llll_ATLAS13
double ie2c[2][CacheSize]
double ip_ex_pp_H_hh_bbbb_ATLAS13_cache_e[2][CacheSize]
gslpp::complex B0p_MZ2_mHl2_mA2_mA2_cache[4][CacheSize]
gslpp::complex B0p_MZ2_mHp2_mHp2_mHl2_cache[4][CacheSize]
double THoEX_VBF_H_ZZ_qqllnunu_ATLAS13
double ip_ex_gg_A_hZ_bbZ_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced pseudoscalar reso...
gslpp::matrix< double > ATLAS13_gg_phi_tautau_e
double ip_ex_pp_phi_hh_gagabb_CMS8_cache[2][CacheSize]
double R_VBF_H_ZZ_llqq_ATLAS13
gslpp::matrix< double > CMS8_gg_A_hZ_bbll_e
double THoEX_ggF_A_tautau_ATLAS8
gslpp::complex B0p_MZ2_mHh2_0_mHp2(const double MZ2, const double mHh2, const double mHp2) const
.
double ip_ex_VBF_H_ZZ_llll_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a vector boson fusion produced scalar re...
double R_pp_H_hh_bbVV_CMS13
gslpp::matrix< double > ATLAS13_bb_phi_tautau_e
double ghHpHm(const double mHp2, const double tanb, const double m12_2, const double bma, const double mHl2, const double vev) const
double ip_ex_pp_phi_gaga_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a scalar resonance decaying to two photo...
double ip_ex_pp_H_ZZ_llqq_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double R_pp_A_Zga_llga_ATLAS8
gslpp::complex A_H_Hp(const double mHp2, const double mH, const double cW2, const double MZ) const
Amplitude for a CP-even Higgs boson decay to a photon and a Z boson including the charged Higgs boson...
double ie39e(double mass)
double ip_cs_pptobbA_8_cache[2][CacheSize]
double ip_ex_VBF_H_ZZ_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a vector boson fusion produced scalar re...
double ie33c[2][CacheSize]
double ie26c[2][CacheSize]
gslpp::complex B00_MZ2_MZ2_MZ2_mHh2_cache[3][CacheSize]
double ip_cs_ggtoA_8_cache[2][CacheSize]
double cW2THDM(const double c02) const
gslpp::complex I_h_L(const double mHl2, const double Me, const double Mmu, const double Mtau) const
Amplitude for the SM Higgs boson decay to diphotons including the leptons in the loop.
gslpp::complex B0p_MZ2_mHl2_0_0(const double MZ2, const double mHl2) const
.
double ip_GammaHPtotSM(double mass)
Interpolating function for the total SM Higgs decay width.
gslpp::matrix< double > CMS13_pp_H_ZZ_llqq_e
double ip_ex_gg_H_ZZ_llnunu_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced scalar resonance de...
gslpp::matrix< double > log_cs_ttH_8
gslpp::complex B0p_MZ2_mA2_0_mHl2(const double MZ2, const double mA2, const double mHl2) const
.
double THoEX_ggF_H_ZZ_llnunu_ATLAS13
gslpp::matrix< double > CMS8_pp_H_hh_gagabb_e
double ip_ex_pp_H_VV_qqqq_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a scalar resonance decaying to two or ...
double ip_ex_pp_H_hh_bblnulnu_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::complex B0_MZ2_0_MZ2_mHl2_cache[3][CacheSize]
double ie37e(double mass)
double THoEX_ggF_A_tautau_CMS8
double THoEX_pp_H_Zga_ATLAS13
gslpp::matrix< double > temp2
double R_pp_H_Zga_llga_ATLAS8
gslpp::complex A_H_Hp_cache[5][CacheSize]
gslpp::matrix< double > ATLAS13_gg_H_ZZ_llnunu
double ip_ex_gg_A_hZ_bbZ_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced pseudoscalar reso...
gslpp::matrix< double > temp20e
double ip_cs_VBFtoH_8(double mass)
Interpolating function for the H production cross section via vector boson fusion at 8 TeV.
gslpp::matrix< double > ATLAS8_bb_phi_tautau
double ip_cs_pptottH_8(double mass)
Interpolating function for the top associated H production cross section at 8 TeV.
gslpp::matrix< double > br_WW
gslpp::complex I_HH_D_cache[4][CacheSize]
gslpp::complex B0_MZ2_mHp2_mHp2_mHh2(const double MZ2, const double mHp2, const double mHh2) const
.
gslpp::matrix< double > ATLAS8_pp_Hpm_taunu_e
gslpp::complex B0_MZ2_mHl2_mHh2_mHh2(const double MZ2, const double mHl2, const double mHh2) const
.
gslpp::complex B0p_MZ2_mHh2_0_mHp2_cache[4][CacheSize]
double THoEX_pp_H_Zga_llga_CMS8
gslpp::matrix< double > br_ZZ
gslpp::matrix< double > br_tautau
double ip_ex_pp_phi_Zga_llga_CMS13_cache_e[2][CacheSize]
double ip_Br_HPtoWW_cache[2][CacheSize]
double R_ggF_A_tautau_CMS8
gslpp::matrix< double > CMS13_pp_phi_Zga_llga
gslpp::matrix< double > temp37e
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_ex_mu_pp_H_VV_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two massive...
gslpp::complex B0p_MZ2_mHp2_mHp2_mHh2_cache[4][CacheSize]
double THoEX_VBF_H_ZZ_ATLAS8
double ip_ex_pp_phi_gaga_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a scalar resonance decaying to two photo...
gslpp::matrix< double > temp15e
gslpp::complex B0_MZ2_mHh2_0_mA2(const double MZ2, const double mHh2, const double mA2) const
.
gslpp::matrix< double > br_bb
double ip_ex_pp_Hpm_tb_ATLAS8_cache_e[2][CacheSize]
double ip_ex_VBF_VH_H_ZZ_llll_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a vector boson fusion produced scalar reso...
double R_ggF_A_Zga_llga_ATLAS13
double ie6ec[2][CacheSize]
double R_bbF_A_tautau_ATLAS8
gslpp::matrix< double > ATLAS13_gg_H_ZZ_nunuqq_e
gslpp::matrix< double > temp19
double ip_ex_gg_A_hZ_tautaull_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced pseudoscalar resona...
double THoEX_pp_H_ZZ_llnunu_CMS13
gslpp::matrix< double > temp39
double ip_ex_gg_phi_tt_ATLAS8_cache_e[2][CacheSize]
gslpp::matrix< double > ATLAS13_pp_H_hh_gagabb_e
gslpp::matrix< double > ATLAS13_bb_phi_tautau
double R_pp_H_VV_qqqq_ATLAS13
gslpp::matrix< double > ATLAS8_pp_phi_gaga
gslpp::complex B0_MZ2_0_MW2_mHl2(const double MZ2, const double MW2, const double mHl2) const
.
double ip_ex_pp_H_AZ_tautaull_CMS8(double mA, double mH)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to which furt...
gslpp::complex B0p_MZ2_mHl2_mHl2_mHl2(const double MZ2, const double mHl2) const
.
double R_ggF_H_tautau_CMS13
gslpp::matrix< double > CMS13_pp_phi_bb
gslpp::complex B0_MZ2_mHl2_mHp2_mHp2_cache[4][CacheSize]
double R_ggF_H_ZZ_qqllnunu_ATLAS13
double THoEX_ggF_H_ZZ_ATLAS8
double pp_H_hh_gagabb_TH13
double ip_ex_VBF_H_ZZ_llnunu_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a VBF produced scalar resonance decaying t...
double THoEX_bbF_A_tautau_CMS8
double ie15e(double mass)
double ip_ex_pp_phi_Zga_llga_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to a boson an...
gslpp::matrix< double > CMS13_pp_H_hh_bbtautau1_e
double R_ggF_H_WW_lnuqq_ATLAS13
gslpp::matrix< double > temp14
double R_VBF_VH_H_ZZ_llll_CMS13
double R_ggF_VBF_H_WW_lnulnu_CMS13
double THoEX_ggF_H_ZZ_nunuqq_ATLAS13
gslpp::complex B00_MZ2_MZ2_mHh2_mA2(const double MZ2, const double mHh2, const double mA2) const
.
double ie7c[2][CacheSize]
double ip_ex_pp_phi_Zga_llga_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a scalar resonance decaying to a boson ...
gslpp::matrix< double > CMS8_pp_Hp_taunu_e
double ip_ex_pp_phi_Zga_llga_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a pseudoscalar resonance decaying to a Z...
double THoEX_pp_H_hh_bbbb_CMS13
gslpp::complex B0_MZ2_mHh2_mA2_mA2(const double MZ2, const double mHh2, const double mA2) const
.
double ie40c[2][CacheSize]
double ip_Br_HPtomumu(double mass)
Interpolating function for the SM branching ratio to two muons.
gslpp::complex B0p_MZ2_0_mA2_mHl2_cache[4][CacheSize]
gslpp::matrix< double > CMS13_pp_H_hh_bbVV
double ie20c[2][CacheSize]
double ip_ex_pp_H_hh_bbtautau_CMS13_cache[2][CacheSize]
gslpp::complex B00_MZ2_0_mHp2_mHp2(const double MZ2, const double mHp2) const
.
gslpp::matrix< double > ATLAS13_VBF_H_WW_lnuqq_e
double ip_ex_VBF_VH_H_ZZ_llll_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a vector boson fusion produced scalar reso...
double ip_ex_gg_H_WW_ATLAS8_cache[2][CacheSize]
gslpp::complex B0_MZ2_0_MZ2_mHh2_cache[3][CacheSize]
double THoEX_pp_H_gaga_ATLAS13
double ip_Br_HPtoZZ_cache[2][CacheSize]
double ie30e(double mass)
gslpp::complex A_HH_D_cache[6][CacheSize]
double ie31e(double mass)
double ip_ex_gg_phi_tautau_ATLAS8_cache_e[2][CacheSize]
gslpp::matrix< double > ATLAS13_bb_phi_tt_e
double ip_ex_pp_Hp_tb_ATLAS13_2_e(double mass)
Interpolating function for the expected ATLAS upper limit on a singly charged scalar resonance decayi...
double R_VBF_H_ZZ_qqllnunu_ATLAS13
double ie13c[2][CacheSize]
double ip_ex_bb_A_Zh_Zbb_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a bb associated pseudoscalar resonance d...
double ip_Br_HPtoZZ(double mass)
Interpolating function for the SM branching ratio to two bosons.
double ip_ex_bb_phi_tautau_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a b associated scalar resonance decaying t...
gslpp::complex B0p_MZ2_0_mA2_mHh2_cache[4][CacheSize]
gslpp::complex B00_MZ2_0_MZ2_mHh2_cache[3][CacheSize]
double THoEX_bbF_A_tautau_ATLAS13
double R_bbF_H_tautau_ATLAS13
double ggF_A_hZ_tautaull_TH8
double ip_ex_pp_Hp_tb_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a singly charged scalar resonance decaying...
gslpp::complex B0p_MZ2_mA2_0_mHp2_cache[4][CacheSize]
gslpp::complex B0p_MZ2_0_0_mHh2(const double MZ2, const double mHh2) const
.
gslpp::matrix< double > ATLAS13_pp_phi_Zga
double ip_cs_WtoWH_13(double mass)
Interpolating function for the W associated H production cross section at 13 TeV.
gslpp::matrix< double > temp40e
double R_bbF_A_tautau_ATLAS13
double ip_ex_gg_A_hZ_bbZ_ATLAS8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS8_gg_H_ZZ
double VBF_H_ZZ_llll_TH13
double ip_ex_pp_Hp_tb_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a singly charged scalar resonance decaying...
gslpp::complex B0_MZ2_mHh2_mHh2_mHh2(const double MZ2, const double mHh2) const
.
double THoEX_ggF_H_gaga_CMS13
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 > temp5e
gslpp::matrix< double > CMS8_pp_Hp_taunu
double ie29ec[2][CacheSize]
gslpp::complex B0_MZ2_MZ2_MZ2_mHh2(const double MZ2, const double mHh2) const
.
double ip_ex_ggVV_H_WW_lnulnu_CMS13_cache_e[2][CacheSize]
double ip_csr_ggA_t_8_cache[2][CacheSize]
double R_bbF_A_hZ_bbZ_ATLAS13
double ip_ex_VBF_H_ZZ_llllnunu_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a VBF produced scalar resonance decaying...
double R_VBF_H_ZZ_llnunu_CMS13
gslpp::matrix< double > CMS8_gg_phi_gaga_e
double ip_ex_gg_A_Zh_Zbb_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > ATLAS8_pp_phi_gaga_e
gslpp::matrix< double > CMS8_pp_H_hh_e
gslpp::complex B0p_MZ2_0_0_mHl2(const double MZ2, const double mHl2) const
.
double R_pp_Hp_tb_ATLAS13_2
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::matrix< double > CMS13_gg_phi_tautau
double ip_ex_pp_H_hh_gagabb_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double ip_ex_gg_A_hZ_tautauZ_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced pseudoscalar reso...
gslpp::matrix< double > ATLAS13_tt_phi_tt
gslpp::matrix< double > temp25
gslpp::matrix< double > log_cs_WH_8
gslpp::complex B0_MZ2_mHl2_mHp2_mHp2(const double MZ2, const double mHl2, const double mHp2) const
.
double ie28c[2][CacheSize]
double THoEX_bbF_H_tautau_CMS13
gslpp::matrix< double > temp10
double ip_ex_gg_H_hh_ATLAS8_cache_e[2][CacheSize]
gslpp::complex B0p_MZ2_mHp2_mHp2_mHh2(const double MZ2, const double mHp2, const double mHh2) const
.
gslpp::complex B00_MZ2_0_mHh2_mA2(const double MZ2, const double mHh2, const double mA2) const
.
gslpp::matrix< double > ATLAS13_tt_phi_tt_e
double ip_ex_gg_phi_gaga_CMS8_cache[2][CacheSize]
double R_ggF_H_tautau_ATLAS13
double ip_ex_gg_H_WW_lnuqq_ATLAS13_cache_e[2][CacheSize]
double ie25c[2][CacheSize]
double ip_ex_gg_A_hZ_tautaull_CMS8_cache[2][CacheSize]
double ip_ex_pp_phi_Zga_qqga_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to a boson an...
gslpp::matrix< double > ATLAS13_gg_H_WW_lnuqq_e
double R_ggF_A_hZ_tautaull_CMS8
double THoEX_ggF_H_WW_lnuqq_ATLAS13
gslpp::matrix< double > ATLAS8_gg_H_ZZ_e
double ip_ex_bb_phi_tautau_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a bb associated scalar resonance decayin...
gslpp::matrix< double > temp5
gslpp::matrix< double > CMS8_pp_H_hh_gagabb
gslpp::matrix< double > ATLAS13_pp_Hp_tb1_e
double R_ggF_A_hZ_bbZ_ATLAS8
double R_ggF_A_hZ_bbll_CMS8
gslpp::complex B0_MZ2_mHp2_0_mHl2(const double MZ2, const double mHp2, const double mHl2) const
.
double ip_ex_gg_H_ZZ_llnunu_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_pp_phi_gaga_ATLAS13_cache[2][CacheSize]
gslpp::complex B0_MZ2_mA2_0_mHl2_cache[4][CacheSize]
double ip_ex_pp_H_AZ_bbll_CMS8_cache[3][CacheSize]
gslpp::complex A_A_U_cache[6][CacheSize]
gslpp::matrix< double > ATLAS13_gg_H_ZZ_llllnunu
gslpp::matrix< double > temp25e
double ip_ex_pp_Hp_tb_CMS8_cache_e[2][CacheSize]
double ie40ec[2][CacheSize]
double ip_ex_pp_H_hh_gagaWW_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::matrix< double > csrA_top_8
gslpp::matrix< double > ATLAS13_gg_H_ZZ_llll
double ip_ex_gg_phi_tautau_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced scalar resonance de...
gslpp::matrix< double > CMS13_bb_phi_tautau_e
double ip_ex_VBF_H_WW_enumunu_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a vector boson fusion produced scalar re...
double ie11ec[2][CacheSize]
double R_pp_H_ZZ_llll_CMS13
gslpp::complex Int1(const double tau, const double lambda) const
function for Z gamma coupling to h, H and A
double ie23c[2][CacheSize]
double ie33e(double mass)
gslpp::matrix< double > arraybsgamma
double R_ggF_H_gaga_CMS13
gslpp::matrix< double > ATLAS8_pp_Hpm_tb
double ie10e(double mass)
double ie19c[2][CacheSize]
gslpp::complex B0_MZ2_mHh2_0_0(const double MZ2, const double mHh2) const
.
double ip_ex_pp_H_hh_bbbb_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a spin-2 resonance decaying to two boso...
gslpp::matrix< double > ATLAS13_VBF_H_ZZ_llllnunu
gslpp::complex B0_MZ2_0_mHp2_mHh2_cache[4][CacheSize]
double ip_ex_pp_Hpm_taunu_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a singly charged scalar resonance decaying...
gslpp::complex B0p_MZ2_0_0_mHh2_cache[3][CacheSize]
double ip_csr_ggA_t_13_cache[2][CacheSize]
gslpp::matrix< double > temp1
gslpp::complex I_H_Hp(const double mHp2, const double mH) const
Amplitude for a CP-even Higgs boson decay to diphotons including the charged Higgs boson in the loop.
gslpp::matrix< double > br_tt
double ie21e(double mass)
gslpp::complex B0p_MZ2_mA2_0_mHh2(const double MZ2, const double mA2, const double mHh2) const
.
gslpp::matrix< double > temp34
double THoEX_ttF_A_tt_ATLAS13
double ip_ex_bsgamma_cache[3][CacheSize]
double ip_ex_VBF_H_ZZ_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a vector boson fusion produced scalar re...
double THoEX_bbF_A_tautau_ATLAS8
gslpp::complex B0p_MZ2_mHl2_mA2_mA2(const double MZ2, const double mHl2, const double mA2) const
.
double ie4c[2][CacheSize]
gslpp::complex B0_MZ2_MW2_MW2_mHl2_cache[4][CacheSize]
double R_ggF_A_hZ_tautauZ_ATLAS8
double THoEX_pp_Hp_tb_ATLAS13_2
double interpolate2D(gslpp::matrix< double > arrayTab, double x, double y)
Linearly interpolates a table with two parameter dimensions.
double ip_ex_pp_H_ZZ_llqq_CMS13_cache_e[2][CacheSize]
gslpp::complex B0p_MZ2_mHl2_0_mA2(const double MZ2, const double mHl2, const double mA2) const
.
double ip_ex_pp_H_hh_CMS8_cache_e[2][CacheSize]
gslpp::complex B00_MZ2_MW2_mHp2_mHp2_cache[4][CacheSize]
gslpp::matrix< double > temp27e
double ie18ec[2][CacheSize]
gslpp::complex B0_MZ2_0_0_mHl2(const double MZ2, const double mHl2) const
.
double R_pp_H_hh_bbtautau_CMS13
void computeAquantities()
gslpp::complex B0p_MZ2_mHh2_0_mA2_cache[4][CacheSize]
double VBF_VH_H_ZZ_llll_TH13
double THoEX_pp_A_Zga_llga_CMS13
gslpp::matrix< double > temp37
double THoEX_pp_Hpm_taunu_ATLAS13
gslpp::matrix< double > GammaHtot_SM
gslpp::complex B0_MZ2_0_mHp2_mHl2_cache[4][CacheSize]
double ip_ex_pp_Hpm_taunu_ATLAS13_cache_e[2][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 ggF_H_hh_bbtautau_TH8
gslpp::matrix< double > ATLAS8_gg_A_hZ_bbZ_e
double ie34e(double mass)
double ip_ex_pp_Hpm_taunu_ATLAS8_cache[2][CacheSize]
double ie30ec[2][CacheSize]
double ip_cs_pptobbH_13_cache[2][CacheSize]
double ip_ex_gg_H_ZZ_llnunu_CMS13_cache[2][CacheSize]
double ie9c[2][CacheSize]
double ie11c[2][CacheSize]
gslpp::complex B0_MZ2_MZ2_MZ2_mHl2_cache[3][CacheSize]
double ie37ec[2][CacheSize]
gslpp::matrix< double > ATLAS8_VBF_H_ZZ
double ip_ex_gg_A_Zh_Zbb_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced pseudoscalar reso...
double THoEX_pp_H_ZZ_llll_CMS13
double ip_ex_VBF_H_ZZ_llll_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a vector boson fusion produced scalar re...
gslpp::matrix< double > CMS13_VBF_H_ZZ_llnunu
double ip_ex_VBF_H_WW_enumunu_ATLAS13_cache_e[2][CacheSize]
double THoEX_ggF_A_hZ_bbZ_ATLAS8
gslpp::matrix< double > ATLAS13_gg_phi_Zga_llga_e
double ip_ex_pp_H_hh_gagaWW_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > temp10e
gslpp::matrix< double > ATLAS13_bb_phi_tt
double ip_ex_gg_A_hZ_tautauZ_ATLAS8_cache[2][CacheSize]
double ie32c[2][CacheSize]
double THoEX_pp_H_AZ_bbll_CMS8
gslpp::matrix< double > ATLAS8_gg_A_hZ_tautauZ_e
double ip_ex_ggF_H_hh_bbbb_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double ip_ex_pp_H_ZZ_llll_CMS13_cache_e[2][CacheSize]
gslpp::matrix< double > temp17
double ip_ex_pp_Hpm_taunu_ATLAS8_cache_e[2][CacheSize]
double THoEX_bbF_H_tt_ATLAS13
double THoEX_ggF_H_tautau_ATLAS8
double R_ggF_A_hZ_bbZ_ATLAS13
double ip_ex_pp_phi_hh_gagabb_CMS8_cache_e[2][CacheSize]
double ip_ex_gg_phi_gaga_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced scalar resonance de...
gslpp::complex A_h_U_cache[7][CacheSize]
double ip_ex_pp_phi_gaga_ATLAS8_cache[2][CacheSize]
double ie36e(double mass)
double ggF_H_ZZ_llll_TH13
double ip_ex_gg_phi_tautau_CMS13_cache_e[2][CacheSize]
gslpp::matrix< double > ATLAS13_gg_H_ZZ_llnunu_e
double ip_ex_bb_phi_tt_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a bb associated scalar resonance decayin...
gslpp::matrix< double > CMS8_pp_H_hh
double ip_ex_pp_Hpm_taunu_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a singly charged scalar resonance decayi...
double ip_ex_gg_phi_Zga_llga_ATLAS13_cache_e[2][CacheSize]
double THoEX_pp_A_Zga_qqga_CMS13
gslpp::matrix< double > log_cs_bbH_13
double THoEX_bbF_H_tautau_ATLAS8
gslpp::complex B0p_MZ2_mHh2_mHh2_mHh2(const double MZ2, const double mHh2) const
.
double ie3ec[2][CacheSize]
double THoEX_ttF_H_tt_ATLAS13
gslpp::matrix< double > ATLAS13_pp_H_hh_gagaWW_e
double ip_ex_pp_H_hh_gagabb_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
double ip_ex_bb_phi_tautau_ATLAS8_cache_e[2][CacheSize]
double THoEX_bbF_A_bb_CMS8
gslpp::matrix< double > temp14e
double R_ggF_H_hh_bbbb_CMS13
gslpp::matrix< double > CMS8_mu_pp_H_VV_e
double ip_ex_VBF_H_ZZ_ATLAS8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS8_gg_A_hZ_tautauZ
gslpp::matrix< double > CMS13_pp_H_hh_gagabb
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_H_hh_CMS8
double ie20ec[2][CacheSize]
double THoEX_pp_A_Zga_llga_CMS8
double THoEX_bbF_H_bb_CMS8
double ip_Br_HPtott_cache[2][CacheSize]
double THoEX_pp_H_hh_gagabb_CMS8
gslpp::matrix< double > log_cs_ttH_13
double ie13e(double mass)
double ip_ex_gg_phi_gaga_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced scalar resonance de...
double ip_ex_pp_Hp_tb_ATLAS13_1(double mass)
Interpolating function for the observed ATLAS upper limit on a singly charged scalar resonance decayi...
gslpp::complex B0_MZ2_mHp2_mHp2_mHl2_cache[4][CacheSize]
double THoEX_ggF_A_gaga_CMS13
double ip_ex_pp_H_hh_bblnulnu_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::matrix< double > temp11e
gslpp::matrix< double > CMS13_pp_H_hh_bbtautau
gslpp::complex A_A_L_cache[6][CacheSize]
double ip_ex_tt_phi_tt_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a tt associated scalar resonance decayin...
gslpp::matrix< double > CMS13_pp_H_hh_bbVV_e
double ip_ex_gg_A_hZ_bbll_CMS8_cache_e[2][CacheSize]
double ip_ex_pp_phi_Zga_qqga_CMS13_cache_e[2][CacheSize]
gslpp::complex B0p_MZ2_mHl2_0_mA2_cache[4][CacheSize]
double ie8ec[2][CacheSize]
gslpp::complex A_A_D(const double mA2, const double cW2, const double Ms, const double Mb, const double MZ) const
Amplitude for a CP-odd Higgs boson decay to a photon and a Z boson including the strange and bottom q...
gslpp::matrix< double > CMS13_pp_H_hh_bblnulnu
double ip_ex_gg_H_ZZ_llll_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double R_pp_H_hh_gagabb_CMS8
gslpp::complex B0p_MZ2_mHh2_mA2_mA2(const double MZ2, const double mHh2, const double mA2) const
.
double ip_ex_pp_Hp_tb_ATLAS13_2_cache_e[2][CacheSize]
double THoEX_ggF_H_Zga_CMS13
double ip_ex_VBF_H_WW_lnuqq_ATLAS13_cache[2][CacheSize]
gslpp::complex A_h_D_cache[7][CacheSize]
void computeWFRcombinations()
double ip_ex_pp_A_Zga_llga_CMS8_cache_e[2][CacheSize]
gslpp::complex B0_MZ2_0_0_mHh2_cache[3][CacheSize]
double ip_Br_HPtobb_cache[2][CacheSize]
double THoEX_pp_A_HZ_bbll_CMS8
double R_ggF_H_ZZ_llqq_ATLAS13
double ip_csr_ggH_b_13_cache[2][CacheSize]
double ip_cs_VBFtoH_8_cache[2][CacheSize]
gslpp::complex A_A_D_cache[6][CacheSize]
gslpp::matrix< double > ATLAS13_gg_H_ZZ_llqq
double ip_ex_pp_H_hh_gagabb_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a scalar resonance decaying to two boso...
double ie17e(double mass)
gslpp::complex I_A_L(const double mA2, const double Mmu, const double Mtau) const
Amplitude for a CP-odd Higgs boson decay to diphotons including muons and taus in the loop.
double ip_ex_pp_H_hh_bbtautau1_CMS13_cache[2][CacheSize]
double ip_ex_pp_H_hh_bbVV_CMS13_cache[2][CacheSize]
gslpp::matrix< double > ATLAS8_gg_phi_tautau
gslpp::complex A_HH_L(const double mHh2, const double cW2, const double Mmu, const double Mtau, const double MZ) const
Amplitude for a heavy CP-even Higgs boson decay to a photon and a Z boson including muons and taus in...
gslpp::matrix< double > ATLAS8_gg_H_hh
double ie32ec[2][CacheSize]
double R_ggF_H_WW_enumunu_ATLAS13
gslpp::matrix< double > temp11
gslpp::matrix< double > temp3e
gslpp::complex B0p_MZ2_0_mHp2_mA2_cache[4][CacheSize]
double ip_cs_pptobbA_8(double mass)
Interpolating function for the bottom associated A production cross section at 8 TeV.
gslpp::complex B0_MZ2_mHp2_mHp2_mHh2_cache[4][CacheSize]
double THoEX_ggF_H_tautau_CMS8
double ip_ex_gg_phi_gaga_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced scalar resonance de...
gslpp::matrix< double > ATLAS8_pp_Hpm_taunu
static const int CacheSize
Cache size.
double ip_ex_pp_H_hh_gagabb_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a scalar resonance decaying to two boso...
gslpp::complex I_h_U(const double mHl2, const double Mu, const double Mc, const double Mt) const
Amplitude for the SM Higgs boson decay to diphotons including the up-type quarks in the loop.
gslpp::matrix< double > ATLAS8_VBF_H_ZZ_e
double ip_cs_ggtoH_8_cache[2][CacheSize]
double THoEX_ggF_H_hh_ATLAS8
double ip_ex_pp_H_ZZ_llnunu_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double ip_ex_bb_phi_tautau_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a bottom quark produced scalar resonance d...
double ip_ex_pp_H_hh_gagabb_ATLAS13_cache_e[2][CacheSize]
gslpp::matrix< double > temp3
double ie12c[2][CacheSize]
gslpp::complex B0p_MZ2_mA2_mA2_mHh2_cache[4][CacheSize]
double ie1c[2][CacheSize]
gslpp::matrix< double > temp13
gslpp::matrix< double > temp29
double ip_ex_pp_phi_gaga_ATLAS8_cache_e[2][CacheSize]
double ip_ex_gg_H_ZZ_llqq_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::complex B0_MZ2_MW2_MW2_mHl2(const double MZ2, const double MW2, const double mHl2) const
.
gslpp::complex B0p_MZ2_mHp2_0_mA2(const double MZ2, const double mHp2, const double mA2) const
.
void read()
Fills all required arrays with the values read from the tables.
double ip_ex_pp_A_Zga_llga_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a pseudoscalar resonance decaying to a Z b...
double ip_ex_pp_H_hh_bbVV_CMS13_cache_e[2][CacheSize]
gslpp::complex B0_MZ2_0_0_mHh2(const double MZ2, const double mHh2) const
.
gslpp::complex B0p_MZ2_mA2_0_mHl2_cache[4][CacheSize]
double ip_ex_pp_H_ZZ_llnunu_CMS13_cache_e[2][CacheSize]
double ie38c[2][CacheSize]
double R_pp_H_gaga_ATLAS8
double R_VBF_H_ZZ_llllnunu_ATLAS13
gslpp::complex A_H_W(const double mH, const double cW2, const double MW, const double MZ) const
Amplitude for a CP-even Higgs boson decay to a photon and a Z boson including the W boson in the loop...
double ip_ex_ggVV_H_WW_lnulnu_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon or vector boson fusion produ...
double ip_ex_tt_phi_tt_ATLAS13_cache[2][CacheSize]
gslpp::complex B00_MZ2_0_mHh2_mHp2(const double MZ2, const double mHh2, const double mHp2) const
.
double ggF_H_hh_bbbb_TH13
gslpp::complex B00_MZ2_0_MW2_mHl2_cache[4][CacheSize]
double ip_ex_pp_H_hh_bbtautau_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::matrix< double > temp39e
gslpp::matrix< double > ATLAS13_pp_Hpm_taunu_e
double THoEX_pp_Hpm_taunu_ATLAS8
gslpp::matrix< double > CMS13_ggFVBF_H_WW_lnulnu_e
double R_pp_A_Zga_llga_CMS13
double ip_ex_VBF_H_ZZ_qqllnunu_ATLAS13_cache_e[2][CacheSize]
double R_pp_H_AZ_bbll_CMS8
double THoEX_pp_Hp_tb_CMS8
gslpp::complex B0_MZ2_mHl2_mHl2_mHl2(const double MZ2, const double mHl2) const
.
gslpp::matrix< double > CMS13_pp_H_hh_bbtautau1
gslpp::complex B00_MZ2_0_mHl2_mHp2_cache[4][CacheSize]
double ip_ex_pp_Hpm_tb_ATLAS8_cache[2][CacheSize]
double ie16ec[2][CacheSize]
gslpp::complex B0p_MZ2_mHh2_0_0(const double MZ2, const double mHh2) const
.
gslpp::matrix< double > ATLAS8_gg_phi_tautau_e
double ip_ex_pp_H_ZZ_llqq_CMS13_cache[2][CacheSize]
gslpp::complex B0p_MZ2_mHp2_0_mHh2(const double MZ2, const double mHp2, const double mHh2) const
.
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_gaga_CMS13_cache[2][CacheSize]
double R_pp_H_hh_bbbb_ATLAS13
double ip_ex_gg_H_ZZ_qqllnunu_ATLAS13_cache[2][CacheSize]
double ip_ex_VBF_VH_H_ZZ_llll_CMS13_cache_e[2][CacheSize]
double THoEX_ggF_H_hh_bbbb_CMS13
gslpp::matrix< double > ATLAS13_VBF_H_WW_enumumu
double ip_ex_gg_phi_tautau_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced scalar resonance de...
double ip_GammaHPtotSM_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_VBF_H_ZZ_qqllnunu
double THoEX_bbF_A_tt_ATLAS13
double ip_ex_bb_A_Zh_Zbb_ATLAS13_cache_e[2][CacheSize]
gslpp::matrix< double > ATLAS8_pp_phi_Zga_llga_e
double ip_ex_pp_H_hh_gagaWW_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::matrix< double > ATLAS13_VBF_H_ZZ_llqq_e
gslpp::matrix< double > br_mumu
double ie34c[2][CacheSize]
double ie5c[2][CacheSize]
double ip_ex_bb_phi_tautau_ATLAS13_cache_e[2][CacheSize]
double pp_H_hh_bblnulnu_TH13
gslpp::complex B0p_MZ2_mHh2_mHh2_mHl2(const double MZ2, const double mHh2, const double mHl2) const
.
gslpp::complex B00_MZ2_MW2_mHl2_mHp2_cache[5][CacheSize]
gslpp::complex B0_MZ2_mHl2_mA2_mA2(const double MZ2, const double mHl2, const double mA2) const
.
gslpp::complex B00_MZ2_MW2_MW2_mHl2(const double MZ2, const double MW2, const double mHl2) const
.
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_pptottA_13(double mass)
Interpolating function for the top associated A production cross section at 13 TeV.
double R_pp_Hp_tb_ATLAS13
double ip_ex_pp_phi_hh_bbbb_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
int HSTheta(const double x) const
Heaviside function.
double THoEX_pp_H_bb_CMS13
double ip_ex_bb_phi_tautau_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a bottom quark produced scalar resonance...
double ip_ex_pp_A_HZ_bbll_CMS8(double mA, double mH)
Interpolating function for the observed CMS upper limit on a pseudoscalar resonance decaying to whic...
gslpp::complex I_HH_D(const double mHh2, const double Ms, const double Mb) const
Amplitude for a heavy CP-even Higgs boson decay to diphotons including the strange and bottom quarks ...
double ie12ec[2][CacheSize]
gslpp::matrix< double > temp36
double ip_cs_ggtoH_13_cache[2][CacheSize]
double ip_ex_pp_A_HZ_bbll_CMS8_cache[3][CacheSize]
double ip_ex_pp_phi_hh_bbbb_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::complex B00_MZ2_MW2_mA2_mHp2(const double MZ2, const double MW2, const double mA2, const double mHp2) const
.
double R_pp_Hp_taunu_CMS8
gslpp::matrix< double > ATLAS8_gg_H_WW_e
double R_ggF_A_tautau_ATLAS8
double ie10ec[2][CacheSize]
gslpp::complex B0_MZ2_mHh2_0_0_cache[3][CacheSize]
double THoEX_pp_H_hh_bblnulnu_CMS13
gslpp::matrix< double > ATLAS13_pp_Hp_tb2
gslpp::complex I_HH_U(const double mHh2, const double Mc, const double Mt) const
Amplitude for a heavy CP-even Higgs boson decay to diphotons including the charm and top quarks in th...
gslpp::matrix< double > csrA_top_13
gslpp::complex B0p_MZ2_mHh2_mHl2_mHl2(const double MZ2, const double mHh2, const double mHl2) const
.
double ip_csr_ggH_t_8_cache[2][CacheSize]
double ie36c[2][CacheSize]
gslpp::matrix< double > ATLAS13_VBF_H_ZZ_llll_e
void computeHHquantities()
double ip_cs_ggtoA_8(double mass)
Interpolating function for the A production cross section via gluon-gluon fusion at 8 TeV.
double ip_ex_VBF_H_ZZ_llqq_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a vector boson fusion produced scalar re...
double ip_ex_pp_H_hh_bbtautau1_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double THoEX_pp_H_hh_bbbb_CMS8
gslpp::complex B0_MZ2_mHl2_mHh2_mHl2_cache[4][CacheSize]
double ip_ex_gg_H_ZZ_ATLAS8_cache[2][CacheSize]
gslpp::matrix< double > log_cs_bbA_13
gslpp::complex B00_MZ2_0_MW2_mHh2(const double MZ2, const double MW2, const double mHh2) const
.
double ip_ex_gg_phi_tt_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_pp_H_hh_bbtautau1_CMS13_cache_e[2][CacheSize]
double ip_ex_VBF_H_WW_enumunu_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a vector boson fusion produced scalar re...
double ip_csr_ggA_t_8(double mass)
Interpolating function for the gluon-gluon fusion A cross section ratio of the top-loop and the total...
gslpp::matrix< double > ATLAS8_VBF_H_WW_e
gslpp::complex B00_MZ2_MW2_MW2_mHl2_cache[4][CacheSize]
double ip_ex_pp_Hpm_taunu_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a singly charged scalar resonance decayi...
gslpp::complex B0p_MZ2_mA2_mA2_mHh2(const double MZ2, const double mA2, const double mHh2) const
.
double ie19e(double mass)
double ie37c[2][CacheSize]
gslpp::complex g_func(const double x) const
g function for the Int1 function
~THDMcache()
THDMcache destructor.
gslpp::complex B0_MZ2_MZ2_MZ2_mHh2_cache[3][CacheSize]
double ip_ex_pp_Hp_tb_ATLAS13_2_cache[2][CacheSize]
gslpp::matrix< double > log_cs_ZH_8
gslpp::matrix< double > temp31
gslpp::matrix< double > temp30
double ip_ex_gg_phi_gaga_CMS13_cache_e[2][CacheSize]
double THoEX_ggF_H_ZZ_llnunu_CMS13
double ip_ex_pp_H_hh_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::matrix< double > ATLAS13_pp_H_hh_gagabb
double ip_ex_pp_H_ZZ_llll_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::complex B00_MZ2_MZ2_MZ2_mHl2_cache[3][CacheSize]
double ip_ex_gg_phi_tautau_CMS8_cache_e[2][CacheSize]
double R_bbF_H_tautau_CMS13
double Br_AtoHpW
branching ratio to a charged Higgs and a boson.
gslpp::matrix< double > CMS8_pp_A_HZ_tautaull
double ip_ex_VBF_H_ZZ_llnunu_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a VBF produced scalar resonance decaying t...
double ie26e(double mass)
double ip_ex_gg_A_Zh_Zbb_ATLAS13_cache_e[2][CacheSize]
gslpp::matrix< double > temp24
double THoEX_pp_H_hh_bbbb_ATLAS13
double ip_ex_gg_A_hZ_tautauZ_ATLAS8_cache_e[2][CacheSize]
gslpp::complex B0_MZ2_mHl2_0_0(const double MZ2, const double mHl2) const
.
double R_pp_Hpm_taunu_ATLAS13
gslpp::matrix< double > temp4
gslpp::complex B00_MZ2_0_MW2_mHh2_cache[4][CacheSize]
double ip_ex_ggF_phi_Zga_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced scalar resonance de...
gslpp::matrix< double > ATLAS13_pp_H_hh_gagaWW
gslpp::matrix< double > CMS13_pp_phi_Zga_qqga_e
double ie4ec[2][CacheSize]
double R_ttF_H_tt_ATLAS13
gslpp::matrix< double > ATLAS13_pp_H_hh_bbbb
double ip_ex_bb_phi_tautau_CMS13_cache[2][CacheSize]
gslpp::matrix< double > CMS8_pp_A_Zga_llga
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 ip_ex_gg_H_ZZ_llllnunu_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double THoEX_ggF_H_WW_enumunu_ATLAS13
double ip_ex_VBF_VH_H_ZZ_llll_CMS13_cache[2][CacheSize]
gslpp::complex B0p_MZ2_mA2_0_mHp2(const double MZ2, const double mA2, const double mHp2) const
.
gslpp::complex I_A_D_cache[4][CacheSize]
double ip_ex_VBF_H_ZZ_llllnunu_ATLAS13_cache[2][CacheSize]
double ip_ex_pp_H_AZ_bbll_CMS8(double mA, double mH)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to which furt...
double ip_cs_VBFtoH_13_cache[2][CacheSize]
gslpp::matrix< double > temp36e
double ip_ex_VBF_H_WW_ATLAS8_cache_e[2][CacheSize]
gslpp::matrix< double > CMS8_pp_H_hh_bbbb_e
double ip_csr_ggA_t_13(double mass)
Interpolating function for the gluon-gluon fusion A cross section ratio of the top-loop and the total...
double ip_ex_mu_pp_H_VV_CMS8_cache[2][CacheSize]
double ip_ex_gg_A_hZ_bbll_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced pseudoscalar resona...
gslpp::complex B0_MZ2_mHp2_0_mHh2_cache[4][CacheSize]
double ip_ex_ggF_H_hh_bbbb_CMS13_cache_e[2][CacheSize]
double ip_ex_gg_H_WW_enumunu_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::matrix< double > CMS13_VBF_H_ZZ_llnunu_e
double THoEX_bbF_H_tautau_CMS8
double THoEX_ggF_H_tautau_CMS13
double R_pp_Hpm_taunu_ATLAS8
double ip_ex_pp_Hp_taunu_CMS8_cache[2][CacheSize]
gslpp::matrix< double > temp23
double R_pp_A_Zga_llga_CMS8
gslpp::matrix< double > temp33
double ie22ec[2][CacheSize]
double THoEX_ggF_H_ZZ_qqllnunu_ATLAS13
gslpp::matrix< double > ATLAS13_pp_phi_gaga
gslpp::matrix< double > CMS13_pp_Hpm_taunu
double R_ggF_H_ZZ_llll_ATLAS13
gslpp::complex A_h_D(const double mHl2, const double cW2, const double Md, const double Ms, const double Mb, const double MZ) const
Amplitude for the SM Higgs boson decay to a photon and a Z boson including the down-type quarks in th...
gslpp::matrix< double > CMS13_gg_phi_gaga
double ip_cs_pptobbA_13(double mass)
Interpolating function for the bottom associated A production cross section at 13 TeV.
gslpp::matrix< double > CMS8_pp_Hp_tb_e
gslpp::complex Int2(const double tau, const double lambda) const
function for Z gamma coupling to h, H and A
double THoEX_pp_H_hh_bbgaga_ATLAS13
double ip_ex_pp_H_VV_qqqq_ATLAS13_cache_e[2][CacheSize]
double ip_ex_gg_H_ZZ_qqllnunu_ATLAS13_cache_e[2][CacheSize]
gslpp::matrix< double > CMS8_gg_A_hZ_tautaull
double ip_Br_HPtocc_cache[2][CacheSize]
gslpp::complex B0_MZ2_0_MW2_mHl2_cache[4][CacheSize]
gslpp::matrix< double > temp32
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 ip_ex_gg_phi_tt_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_pp_H_hh_gagabb_CMS13_cache_e[2][CacheSize]
double ip_ex_pp_Hp_taunu_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a singly charged scalar resonance decaying...
double THoEX_ggF_A_Zga_llga_ATLAS13
gslpp::complex B0_MZ2_mHp2_mHp2_mHl2(const double MZ2, const double mHp2, const double mHl2) const
.
double ie20e(double mass)
double R_pp_H_hh_bbgaga_CMS13
double ie17ec[2][CacheSize]
double ip_ex_bb_A_Zh_Zbb_ATLAS13_cache[2][CacheSize]
gslpp::complex f_func(const double x) const
f function for the gamma gamma coupling to h, H and A
gslpp::matrix< double > ATLAS13_gg_H_WW_enumumu_e
double R_ggF_H_ZZ_nunuqq_ATLAS13
double ip_ex_gg_A_hZ_bbll_CMS8_cache[2][CacheSize]
gslpp::complex B0_MZ2_mA2_0_mHl2(const double MZ2, const double mA2, const double mHl2) const
.
gslpp::matrix< double > ATLAS8_gg_phi_tt
gslpp::complex B0p_MZ2_mHl2_mHh2_mHh2_cache[4][CacheSize]
gslpp::complex B0p_MZ2_0_mHp2_mHh2_cache[4][CacheSize]
double ip_ex_pp_H_hh_bbtautau_CMS13_cache_e[2][CacheSize]
gslpp::matrix< double > log_cs_ggHp_8
double ip_ex_pp_phi_Zga_llga_CMS13_cache[2][CacheSize]
gslpp::complex B00_MZ2_MZ2_mHh2_mA2_cache[4][CacheSize]
gslpp::complex B0_MZ2_0_mA2_mHl2_cache[4][CacheSize]
gslpp::matrix< double > ATLAS13_gg_phi_Zga_llga
double ip_ex_pp_H_hh_bbVV_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::complex B0_MZ2_mHl2_mA2_mA2_cache[4][CacheSize]
gslpp::matrix< double > CMS13_pp_H_ZZ_llll
gslpp::matrix< double > CMS8_pp_H_hh_bbbb
gslpp::complex B0_MZ2_mHh2_mHp2_mHp2(const double MZ2, const double mHh2, const double mHp2) const
.
gslpp::complex B0p_MZ2_0_mA2_mHh2(const double MZ2, const double mA2, const double mHh2) const
.
gslpp::complex I_h_D_cache[5][CacheSize]
gslpp::complex B00_MZ2_0_mHh2_mHp2_cache[4][CacheSize]
gslpp::matrix< double > ATLAS13_gg_A_Zh_Zbb
double ie30c[2][CacheSize]
gslpp::complex B0p_MZ2_mHh2_mHp2_mHp2_cache[4][CacheSize]
double THoEX_VBF_H_ZZ_llllnunu_ATLAS13
gslpp::matrix< double > temp27
double R_pp_A_Zga_llga_ATLAS13
double ip_ex_pp_phi_Zga_qqga_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to a boson an...
double THoEX_pp_H_hh_bbtautau1_CMS13
gslpp::complex B0_MZ2_0_mA2_mHh2_cache[4][CacheSize]
gslpp::matrix< double > CMS13_pp_H_ZZ_llqq
double ip_cs_pptottH_8_cache[2][CacheSize]
gslpp::matrix< double > temp13e
double GammaHptot
Total decay width of the charged Higgs .
gslpp::complex B00_MZ2_MW2_mHh2_mHp2(const double MZ2, const double MW2, const double mHh2, const double mHp2) const
.
gslpp::matrix< double > log_cs_ggHp_13
gslpp::complex B0p_MZ2_mA2_mA2_mHl2(const double MZ2, const double mA2, const double mHl2) const
.
double ip_ex_bsgamma(double logtb, double logmHp)
Interpolating function for the NNLO value for the branching ratio of decays in the THDM.
double ie31c[2][CacheSize]
double ip_ex_bb_phi_tt_ATLAS13_cache_e[2][CacheSize]
double ip_ex_gg_H_WW_lnuqq_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double R_ggF_A_gaga_CMS13
double ip_ex_gg_phi_gaga_CMS8_cache_e[2][CacheSize]
gslpp::complex I_HH_L_cache[4][CacheSize]
gslpp::complex A_HH_U_cache[6][CacheSize]
double THoEX_bbF_A_hZ_bbZ_ATLAS13
gslpp::matrix< double > temp21e
double ip_ex_pp_H_hh_bbbb_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > CMS8_bb_phi_tautau_e
double ip_ex_VBF_H_WW_ATLAS8_cache[2][CacheSize]
gslpp::matrix< double > CMS8_pp_Hp_tb
double ip_ex_gg_H_hh_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_VBF_H_ZZ_llqq_ATLAS13_cache[2][CacheSize]
double ie25e(double mass)
double ip_ex_gg_H_hh_bbtautau_CMS8_cache[2][CacheSize]
gslpp::complex B0_MZ2_0_mHp2_mHh2(const double MZ2, const double mHp2, const double mHh2) const
.
double R_pp_H_Zga_qqga_CMS13
gslpp::matrix< double > temp6e
double ip_ex_gg_H_ZZ_llqq_ATLAS13_cache_e[2][CacheSize]
gslpp::complex B0p_MZ2_mHl2_mHh2_mHl2_cache[4][CacheSize]
gslpp::complex B0p_MZ2_0_mHp2_mA2(const double MZ2, const double mHp2, const double mA2) const
.
double ip_ex_gg_H_ZZ_nunuqq_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ie22e(double mass)
double ip_Br_HPtotautau_cache[2][CacheSize]
gslpp::complex B0p_MZ2_mHh2_mA2_mA2_cache[4][CacheSize]
gslpp::matrix< double > temp24e
double ip_ex_pp_phi_hh_bbbb_CMS8_cache_e[2][CacheSize]
double ie10c[2][CacheSize]
double ip_ex_VBF_H_ZZ_llll_ATLAS13_cache_e[2][CacheSize]
gslpp::complex B0_MZ2_mHh2_mHh2_mHl2_cache[4][CacheSize]
gslpp::matrix< double > temp30e
double ie36ec[2][CacheSize]
gslpp::matrix< double > temp7e
double pp_A_HZ_tautaull_TH8
double THoEX_pp_H_hh_bbtautau_CMS13
gslpp::matrix< double > temp19e
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_ex_gg_phi_gaga_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced scalar resonance de...
THDMcache(const StandardModel &SM_i)
THDMcache constructor.
double pp_H_hh_gagabb_TH8
gslpp::matrix< double > CMS8_gg_A_hZ_tautaull_e
double ie34ec[2][CacheSize]
double R_bbF_H_tt_ATLAS13
double R_ggF_H_tautau_ATLAS8
double R_ggF_H_ZZ_llllnunu_ATLAS13
double ip_ex_gg_H_ZZ_llqq_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::complex A_H_W_cache[5][CacheSize]
double ip_ex_gg_H_hh_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::matrix< double > log_cs_ggH_13
double ip_cs_ZtoZH_13_cache[2][CacheSize]
gslpp::matrix< double > log_cs_ggA_8
gslpp::complex B0p_MZ2_0_0_mHl2_cache[3][CacheSize]
gslpp::complex B00_MZ2_MZ2_mHl2_mA2_cache[4][CacheSize]
double ip_ex_pp_H_hh_bbbb_CMS13_cache_e[2][CacheSize]
gslpp::matrix< double > ATLAS13_gg_H_WW_lnuqq
gslpp::complex A_HH_D(const double mHh2, const double cW2, const double Ms, const double Mb, const double MZ) const
Amplitude for a heavy CP-even Higgs boson decay to a photon and a Z boson including the strange and b...
gslpp::matrix< double > log_cs_VBF_8
gslpp::matrix< double > log_cs_ttA_8
gslpp::matrix< double > CMS13_pp_H_ZZ_llll_e
double ie3c[2][CacheSize]
double R_bbF_H_tautau_CMS8
gslpp::complex B0_MZ2_mHh2_mHl2_mHl2(const double MZ2, const double mHh2, const double mHl2) const
.
double R_pp_H_hh_bbtautau1_CMS13
double ie24c[2][CacheSize]
double pp_H_hh_bbtautau_TH13
double g_HH_HpHm_cache[7][CacheSize]
double ip_cs_pptobbH_13(double mass)
Interpolating function for the bottom associated H production cross section at 13 TeV.
gslpp::complex B00_MZ2_0_MZ2_mHh2(const double MZ2, const double mHh2) const
.
double THoEX_pp_H_Zga_llga_ATLAS8
gslpp::matrix< double > ATLAS8_gg_H_hh_e
double THoEX_pp_Hp_taunu_CMS8
double ip_ex_pp_H_hh_gagabb_CMS13_cache[2][CacheSize]
double ip_ex_pp_H_VV_qqqq_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a scalar resonance decaying to two or ...
double ip_ex_gg_H_WW_lnuqq_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_ggVV_H_WW_lnulnu_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon or vector boson fusion produ...
gslpp::matrix< double > ATLAS13_pp_phi_gaga_e
gslpp::matrix< double > temp17e
gslpp::complex B0p_MZ2_0_mHp2_mHh2(const double MZ2, const double mHp2, const double mHh2) const
.
gslpp::matrix< double > ATLAS13_gg_H_WW_enumumu
double ip_ex_pp_H_hh_bbtautau_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double ie24e(double mass)
double ip_ex_VBF_H_WW_lnuqq_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a VBF produced scalar resonance decaying...
double R_pp_Hpm_taunu_CMS13
gslpp::matrix< double > ATLAS13_VBF_H_ZZ_qqllnunu_e
gslpp::matrix< double > CMS8_bb_phi_tautau
double ip_ex_gg_H_ZZ_llllnunu_ATLAS13_cache[2][CacheSize]
double ip_ex_bb_phi_tt_ATLAS13_cache[2][CacheSize]
double ip_ex_gg_H_WW_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_gg_phi_Zga_llga_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double THoEX_ggF_H_tautau_ATLAS13
gslpp::matrix< double > CMS13_ggF_phi_Zga
double ip_ex_gg_H_ZZ_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double THoEX_ggF_A_gaga_CMS8
double ip_ex_gg_H_WW_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
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 ip_cs_pptobbH_8(double mass)
Interpolating function for the bottom associated H production cross section at 8 TeV.
gslpp::matrix< double > log_cs_ZH_13
gslpp::matrix< double > CMS8_pp_H_AZ_tautaull
gslpp::matrix< double > ATLAS13_gg_H_ZZ_llllnunu_e
double THoEX_ggF_A_tautau_CMS13
double ip_ex_bb_phi_tautau_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a bottom quark produced scalar resonance...
double R_pp_H_gaga_ATLAS13
double ip_ex_pp_phi_gaga_ATLAS13_cache_e[2][CacheSize]
double THoEX_ggF_A_hZ_tautauZ_ATLAS8
double THoEX_VBF_H_WW_lnuqq_ATLAS13
double ip_ex_pp_Hp_tb_CMS8_cache[2][CacheSize]
double ip_csr_ggH_t_8(double mass)
Interpolating function for the gluon-gluon fusion H cross section ratio of the top-loop and the total...
double ie26ec[2][CacheSize]
double THoEX_pp_A_Zga_llga_ATLAS8
double ip_csr_ggH_t_13(double mass)
Interpolating function for the gluon-gluon fusion H cross section ratio of the top-loop and the total...
double THoEX_ggF_A_tt_ATLAS8
gslpp::complex B0_MZ2_0_mA2_mHh2(const double MZ2, const double mA2, const double mHh2) const
.
gslpp::complex B0_MZ2_0_MW2_mHh2_cache[4][CacheSize]
double R_bbF_A_tt_ATLAS13
gslpp::matrix< double > CMS13_ggFVBF_H_WW_lnulnu
double ip_ex_gg_phi_tautau_CMS13_cache[2][CacheSize]
double ip_ex_gg_H_ZZ_llll_ATLAS13_cache_e[2][CacheSize]
double R_ggF_A_tautau_CMS13
gslpp::matrix< double > ATLAS13_pp_H_hh_bbbb_e
double THoEX_pp_H_Zga_llga_CMS13
double ip_ex_pp_H_hh_gagaWW_ATLAS13_cache_e[2][CacheSize]
double R_pp_H_Zga_llga_CMS8
gslpp::complex B0_MZ2_0_MZ2_mHh2(const double MZ2, const double mHh2) const
.
gslpp::complex B0_MZ2_0_0_mHl2_cache[3][CacheSize]
double R_bbF_H_tautau_ATLAS8
double ie14c[2][CacheSize]
double ip_csr_ggH_b_8_cache[2][CacheSize]
gslpp::complex B0p_MZ2_0_mHp2_mHl2(const double MZ2, const double mHp2, const double mHl2) const
.
gslpp::complex B0_MZ2_mHp2_0_mHh2(const double MZ2, const double mHp2, const double mHh2) const
.
gslpp::complex B0_MZ2_mHh2_mHh2_mHh2_cache[3][CacheSize]
double ip_ex_gg_H_ZZ_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::complex A_h_L_cache[7][CacheSize]
gslpp::matrix< double > CMS8_bb_phi_bb_e
double ie15c[2][CacheSize]
gslpp::matrix< double > ATLAS8_gg_H_WW
double ggF_VBF_H_WW_lnulnu_TH13
void computeSignalStrengthQuantities()
gslpp::matrix< double > temp21
gslpp::matrix< double > ATLAS13_pp_Hp_tb1
double ip_ex_gg_H_ZZ_llnunu_CMS13_cache_e[2][CacheSize]
double ip_ex_tt_phi_tt_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a tt associated scalar resonance decayin...
double ip_ex_VBF_H_ZZ_llnunu_CMS13_cache_e[2][CacheSize]
gslpp::matrix< double > ATLAS13_gg_H_ZZ_qqllnunu_e
double ip_cs_VBFtoH_13(double mass)
Interpolating function for the H production cross section via vector boson fusion at 13 TeV.
double ip_ex_pp_H_hh_CMS8_cache[2][CacheSize]
gslpp::complex B0p_MZ2_mHp2_0_mHl2(const double MZ2, const double mHp2, const double mHl2) const
.
double R_pp_A_HZ_bbll_CMS8
double ggF_A_hZ_tautauZ_TH8
gslpp::matrix< double > ATLAS13_pp_H_VV_qqqq_e
double ip_ex_pp_H_ZZ_llqq_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::matrix< double > ATLAS13_gg_H_ZZ_qqllnunu
double THoEX_VBF_VH_H_ZZ_llll_CMS13
double ip_ex_pp_phi_Zga_llga_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a scalar resonance decaying to a boson ...
gslpp::matrix< double > temp7
double THoEX_ggF_H_hh_bbtautau_CMS8
double ip_ex_pp_H_hh_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double ip_ex_bb_A_Zh_Zbb_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a bb associated pseudoscalar resonance d...
double g_HH_HpHm(const double mHp2, const double mHh2, const double tanb, const double m12_2, const double bma, const double vev) const
double ie23e(double mass)
double ip_ex_pp_phi_Zga_llga_ATLAS8_cache[2][CacheSize]
double ip_ex_ggF_H_hh_bbbb_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::complex B00_MZ2_0_mHp2_mHp2_cache[3][CacheSize]
gslpp::complex A_h_U(const double mHl2, const double cW2, const double Mu, const double Mc, const double Mt, const double MZ) const
Amplitude for the SM Higgs boson decay to a photon and a Z boson including the up-type quarks in the ...
double ip_ex_gg_phi_tautau_ATLAS13_cache_e[2][CacheSize]
gslpp::matrix< double > ATLAS8_pp_Hpm_tb_e
gslpp::complex B0p_MZ2_mHh2_0_mA2(const double MZ2, const double mHh2, const double mA2) const
.
double ip_ex_pp_Hpm_taunu_CMS13_cache_e[2][CacheSize]
gslpp::complex B0p_MZ2_mHh2_mHh2_mHl2_cache[4][CacheSize]
double ip_ex_bb_phi_tautau_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a bb associated scalar resonance decayin...
gslpp::complex B0_MZ2_mA2_0_mHh2_cache[4][CacheSize]
double pp_H_AZ_tautaull_TH8
double ip_ex_gg_phi_tautau_ATLAS13_cache[2][CacheSize]
gslpp::complex B0_MZ2_MW2_MW2_mHh2_cache[4][CacheSize]
gslpp::complex B00_MZ2_0_MZ2_mHl2_cache[3][CacheSize]
gslpp::complex B00_MZ2_MZ2_mHp2_mHp2_cache[3][CacheSize]
gslpp::matrix< double > CMS8_gg_phi_gaga
gslpp::matrix< double > temp28
double ip_ex_gg_H_WW_ATLAS8_cache_e[2][CacheSize]
double ie2ec[2][CacheSize]
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 ip_ex_pp_Hp_taunu_CMS8_cache_e[2][CacheSize]
double THoEX_ggF_H_ZZ_llllnunu_ATLAS13
double ip_cs_ZtoZH_8(double mass)
Interpolating function for the Z associated H production cross section at 8 TeV.
gslpp::complex A_HH_U(const double mHh2, const double cW2, const double Mc, const double Mt, const double MZ) const
Amplitude for a heavy CP-even Higgs boson decay to a photon and a Z boson including the charm and top...
double ip_ex_VBF_H_WW_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a vector boson fusion produced scalar re...
double ie11e(double mass)
double ip_ex_pp_phi_hh_gagabb_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::matrix< double > temp15
gslpp::complex B0_MZ2_mHh2_mA2_mA2_cache[4][CacheSize]
double THoEX_pp_Hp_tb_ATLAS13
double THoEX_bbF_A_tautau_CMS13
double ip_ex_pp_H_ZZ_llnunu_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::complex I_h_L_cache[5][CacheSize]
double ip_Br_HPtotautau(double mass)
Interpolating function for the SM branching ratio to two tau leptons.
gslpp::matrix< double > ATLAS13_VBF_H_ZZ_llllnunu_e
double ip_ex_VBF_H_ZZ_llnunu_CMS13_cache[2][CacheSize]
double ie23ec[2][CacheSize]
gslpp::complex I_H_W(const double mH, const double MW) const
Amplitude for a CP-even Higgs boson decay to diphotons including the W boson in the loop.
gslpp::matrix< double > CMS13_pp_H_hh_bbbb
gslpp::matrix< double > temp33e
double ip_cs_WtoWH_13_cache[2][CacheSize]
double THoEX_ggF_A_Zga_CMS13
gslpp::matrix< double > CMS13_gg_phi_tautau_e
gslpp::matrix< double > csrH_bottom_8
gslpp::complex B0_MZ2_mHl2_0_mA2(const double MZ2, const double mHl2, const double mA2) const
.
gslpp::complex I_h_U_cache[5][CacheSize]
double ip_ex_pp_phi_Zga_llga_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to a boson an...
gslpp::matrix< double > CMS8_gg_A_hZ_bbll
double THoEX_mu_pp_H_VV_CMS8
gslpp::matrix< double > CMS13_pp_phi_Zga_qqga
double ip_ex_gg_H_WW_enumunu_ATLAS13_cache_e[2][CacheSize]
gslpp::matrix< double > CMS13_ggF_H_hh_bbbb
double ip_ex_pp_Hp_tb_ATLAS13_1_cache_e[2][CacheSize]
double MWTHDM(const double MW) const
double ip_ex_gg_H_ZZ_nunuqq_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::complex B0p_MZ2_mHl2_mHh2_mHh2(const double MZ2, const double mHl2, const double mHh2) const
.
gslpp::complex B00_MZ2_MW2_MW2_mHh2_cache[4][CacheSize]
double ip_ex_gg_phi_tautau_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double THoEX_pp_A_bb_CMS13
double R_ggF_H_tautau_CMS8
gslpp::matrix< double > ATLAS13_pp_Hp_tb2_e
double ip_ex_pp_H_VV_qqqq_ATLAS13_cache[2][CacheSize]
double ip_ex_gg_A_Zh_Zbb_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced pseudoscalar reso...
double ie39ec[2][CacheSize]
double ip_ex_pp_Hp_tb_ATLAS13_2(double mass)
Interpolating function for the observed ATLAS upper limit on a singly charged scalar resonance decayi...
double THoEX_pp_Hpm_tb_ATLAS8
double ie6c[2][CacheSize]
double ie28ec[2][CacheSize]
double ip_ex_pp_Hp_tb_ATLAS13_1_cache[2][CacheSize]
double ip_ex_gg_H_hh_ATLAS8_cache[2][CacheSize]
double ip_Br_HPtott(double mass)
Interpolating function for the SM branching ratio to two top quarks.
gslpp::matrix< double > CMS8_gg_phi_tautau_e
double ie14ec[2][CacheSize]
double ip_cs_ggtoA_13_cache[2][CacheSize]
double R_pp_Hpm_tb_ATLAS8
gslpp::complex B00_MZ2_MZ2_MZ2_mHh2(const double MZ2, const double mHh2) const
.
gslpp::matrix< double > CMS13_gg_H_ZZ_llnunu
gslpp::complex I_H_W_cache[3][CacheSize]
gslpp::complex I_H_Hp_cache[3][CacheSize]
gslpp::matrix< double > CMS8_mu_pp_H_VV
double ip_ex_gg_H_hh_bbtautau_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced scalar resonance de...
double THoEX_VBF_H_ZZ_llnunu_CMS13
double ie22c[2][CacheSize]
gslpp::matrix< double > temp18
double ie18e(double mass)
double ip_ex_gg_H_WW_lnuqq_ATLAS13_cache[2][CacheSize]
double THoEX_ggF_A_tautau_ATLAS13
double THoEX_ggF_VBF_H_WW_lnulnu_CMS13
gslpp::matrix< double > CMS13_pp_phi_Zga_llga_e
double R_VBF_H_WW_enumunu_ATLAS13
double ip_ex_VBF_H_ZZ_ATLAS8_cache_e[2][CacheSize]
gslpp::matrix< double > ATLAS13_pp_Hpm_taunu
gslpp::complex B00_MZ2_0_MZ2_mHl2(const double MZ2, const double mHl2) const
.
double ip_ex_ggF_H_hh_bbbb_CMS13_cache[2][CacheSize]
gslpp::matrix< double > log_cs_ggA_13
double R_bbF_A_tautau_CMS13
gslpp::complex B0_MZ2_MW2_MW2_mHh2(const double MZ2, const double MW2, const double mHh2) const
.
double ie14e(double mass)
double THoEX_pp_H_hh_bbVV_CMS13
gslpp::matrix< double > temp35e
gslpp::complex B00_MZ2_MW2_MW2_mHh2(const double MZ2, const double MW2, const double mHh2) const
.
double R_ggF_H_Zga_llga_ATLAS13
double ie16c[2][CacheSize]
gslpp::matrix< double > CMS13_bb_phi_tautau
gslpp::complex I_A_L_cache[4][CacheSize]
gslpp::complex I_HH_U_cache[4][CacheSize]
double ip_cs_pptobbA_13_cache[2][CacheSize]
double ie27e(double mass)
gslpp::complex B00_MZ2_0_mHh2_mA2_cache[4][CacheSize]
gslpp::matrix< double > temp12
double THoEX_ggF_H_Zga_llga_ATLAS13
double ip_ex_VBF_H_WW_lnuqq_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a VBF produced scalar resonance decaying...
gslpp::complex B0_MZ2_mA2_mA2_mHl2(const double MZ2, const double mA2, const double mHl2) const
.
double ip_ex_pp_phi_hh_bbbb_CMS8_cache[2][CacheSize]
gslpp::complex B00_MZ2_0_mHl2_mA2_cache[4][CacheSize]
double THoEX_ggF_A_hZ_bbZ_ATLAS13
gslpp::matrix< double > ATLAS13_pp_H_VV_qqqq
double ip_cs_ggtoHp_13_cache[3][CacheSize]
gslpp::complex B0_MZ2_mHl2_0_mHp2_cache[4][CacheSize]
gslpp::matrix< double > temp4e
double R_pp_H_ZZ_llnunu_CMS13
double ie25ec[2][CacheSize]
double ip_ex_gg_A_hZ_bbll_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced pseudoscalar resona...
double ip_cs_WtoWH_8_cache[2][CacheSize]
gslpp::complex A_h_L(const double mHl2, const double cW2, const double Me, const double Mmu, const double Mtau, const double MZ) const
Amplitude for the SM Higgs boson decay to a photon and a Z boson including the leptons in the loop.
double ip_ex_VBF_H_WW_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a vector boson fusion produced scalar re...
double ip_cs_pptottH_13_cache[2][CacheSize]
double ip_ex_pp_H_hh_bbbb_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
double ip_ex_pp_phi_Zga_llga_ATLAS8_cache_e[2][CacheSize]
double THoEX_VBF_H_ZZ_llqq_ATLAS13
double ie7ec[2][CacheSize]
gslpp::matrix< double > temp26
gslpp::matrix< double > log_cs_bbA_8
gslpp::matrix< double > temp34e
double ip_ex_gg_H_ZZ_qqllnunu_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::matrix< double > csrH_top_13
double ip_ex_VBF_H_ZZ_qqllnunu_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a VBF produced scalar resonance decaying...
double ip_ex_bb_phi_tautau_ATLAS13_cache[2][CacheSize]
double ip_ex_pp_H_ZZ_llll_CMS13_cache[2][CacheSize]
double ip_ex_bb_phi_tautau_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a bottom quark produced scalar resonance d...
double R_pp_H_ZZ_llqq_CMS13
double ip_ex_gg_H_ZZ_llllnunu_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::matrix< double > log_cs_VBF_13
double THoEX_ggF_H_gaga_CMS8
gslpp::matrix< double > CMS8_pp_H_AZ_bbll
double R_pp_A_Zga_qqga_CMS13
double ip_cs_pptobbH_8_cache[2][CacheSize]
gslpp::complex B0_MZ2_0_mHp2_mHl2(const double MZ2, const double mHp2, const double mHl2) const
.
double ie9ec[2][CacheSize]
double ip_ex_gg_H_hh_bbtautau_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced scalar resonance de...
double ip_ex_bb_phi_bb_CMS8_cache_e[2][CacheSize]
double ip_csr_ggH_t_13_cache[2][CacheSize]
gslpp::complex B00_MZ2_MZ2_mHl2_mA2(const double MZ2, const double mHl2, const double mA2) const
.
double ip_ex_pp_phi_bb_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to a b quark p...
double THoEX_pp_A_Zga_llga_ATLAS13
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 ie17c[2][CacheSize]
double ip_cs_ggtoA_13(double mass)
Interpolating function for the A production cross section via gluon-gluon fusion at 13 TeV.
double ip_ex_gg_phi_tautau_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::complex B0_MZ2_mA2_mA2_mHl2_cache[4][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_Br_HPtocc(double mass)
Interpolating function for the SM branching ratio to two charm quarks.
double ip_ex_pp_phi_Zga_qqga_CMS13_cache[2][CacheSize]
double ie35ec[2][CacheSize]
gslpp::complex B0_MZ2_mA2_mA2_mHh2(const double MZ2, const double mA2, const double mHh2) const
.
double ip_ex_gg_phi_tautau_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced scalar resonance de...
gslpp::matrix< double > temp22
gslpp::matrix< double > ATLAS13_VBF_H_WW_enumumu_e
double ie13ec[2][CacheSize]
gslpp::complex B0p_MZ2_mHh2_mHh2_mHh2_cache[3][CacheSize]
gslpp::complex B0p_MZ2_mHl2_0_mHp2_cache[4][CacheSize]
gslpp::complex B00_MZ2_MW2_mA2_mHp2_cache[5][CacheSize]
double ip_ex_pp_A_Zga_llga_CMS8_cache[2][CacheSize]
gslpp::complex B00_MZ2_MW2_mHp2_mHp2(const double MZ2, const double MW2, const double mHp2) const
.
gslpp::matrix< double > ATLAS8_gg_A_hZ_bbZ
double interpolate(gslpp::matrix< double > arrayTab, double x)
Linearly interpolates a table with one parameter dimension.
gslpp::matrix< double > temp38
gslpp::complex B00_MZ2_0_MW2_mHl2(const double MZ2, const double MW2, const double mHl2) const
.
double THoEX_pp_H_VV_qqqq_ATLAS13
gslpp::matrix< double > log_cs_bbH_8
gslpp::complex B0_MZ2_0_mA2_mHl2(const double MZ2, const double mA2, const double mHl2) const
.
double ip_Br_HPtomumu_cache[2][CacheSize]
double ip_ex_pp_H_hh_bbVV_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double ie1ec[2][CacheSize]
double R_ggF_H_ZZ_llnunu_ATLAS13
double ie27c[2][CacheSize]
gslpp::complex I_A_U_cache[4][CacheSize]
double R_ggF_A_tautau_ATLAS13
double ip_ex_gg_H_ZZ_llllnunu_ATLAS13_cache_e[2][CacheSize]
gslpp::matrix< double > ATLAS13_gg_H_ZZ_llqq_e
gslpp::complex B0_MZ2_mHh2_0_mHp2_cache[4][CacheSize]
double ip_ex_pp_H_ZZ_llll_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double ie28e(double mass)
double R_pp_H_hh_bblnulnu_CMS13
double ie15ec[2][CacheSize]
double R_bbF_A_tautau_CMS8
gslpp::matrix< double > ATLAS13_bb_A_Zh_Zbb
gslpp::matrix< double > ATLAS8_bb_phi_tautau_e
gslpp::matrix< double > CMS8_gg_H_hh_bbtautau
gslpp::matrix< double > ATLAS13_gg_H_ZZ_nunuqq
double KaellenFunction_cache[4][CacheSize]
double ip_ex_gg_A_hZ_tautaull_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced pseudoscalar resona...
double ip_cs_pptottA_8_cache[2][CacheSize]
gslpp::matrix< double > temp9
gslpp::matrix< double > temp16
double ie38e(double mass)
double ip_ex_ggF_phi_Zga_CMS13_cache_e[2][CacheSize]
double ip_ex_gg_H_ZZ_ATLAS8_cache_e[2][CacheSize]
gslpp::complex B0_MZ2_mHh2_0_mHp2(const double MZ2, const double mHh2, const double mHp2) const
.
double ie38ec[2][CacheSize]
double ip_ex_ggVV_H_WW_lnulnu_CMS13_cache[2][CacheSize]
double ie5ec[2][CacheSize]
gslpp::matrix< double > CMS13_VBFVH_H_ZZ_llll_e
double ie29c[2][CacheSize]
double ie33ec[2][CacheSize]
double ip_ex_mu_pp_H_VV_CMS8_cache_e[2][CacheSize]
gslpp::matrix< double > temp8e
double ip_ex_gg_H_ZZ_llnunu_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::matrix< double > CMS8_bb_phi_bb
gslpp::complex B0p_MZ2_mHl2_0_mHp2(const double MZ2, const double mHl2, const double mHp2) const
.
double ie31ec[2][CacheSize]
double ip_ex_pp_A_Zga_llga_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a pseudoscalar resonance decaying to a Z b...
double ip_csr_ggA_b_8_cache[2][CacheSize]
double ip_ex_bb_phi_tautau_CMS13_cache_e[2][CacheSize]
gslpp::complex B0_MZ2_mHh2_mHp2_mHp2_cache[4][CacheSize]
double ip_ex_gg_H_ZZ_nunuqq_ATLAS13_cache_e[2][CacheSize]
double THoEX_ggF_H_hh_gagaWW_ATLAS13
double ip_ex_VBF_H_ZZ_llqq_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a vector boson fusion produced scalar re...
gslpp::complex B0_MZ2_0_MW2_mHh2(const double MZ2, const double MW2, const double mHh2) const
.
gslpp::complex B0p_MZ2_mHh2_mHl2_mHl2_cache[4][CacheSize]
gslpp::matrix< double > CMS8_gg_H_hh_bbtautau_e
gslpp::complex B00_MZ2_MW2_mHl2_mHp2(const double MZ2, const double MW2, const double mHl2, const double mHp2) const
.
double R_ggF_H_hh_gagaWW_ATLAS13
gslpp::complex B0_MZ2_mHp2_0_mHl2_cache[4][CacheSize]
gslpp::matrix< double > temp20
gslpp::complex I_A_D(const double mA2, const double Ms, const double Mb) const
Amplitude for a CP-odd Higgs boson decay to diphotons including the strange and bottom quarks in the ...
double ie27ec[2][CacheSize]
gslpp::matrix< double > ATLAS13_gg_phi_tautau
double ie19ec[2][CacheSize]
double ip_ex_pp_H_hh_bbbb_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a spin-2 resonance decaying to two boso...
gslpp::matrix< double > CMS13_pp_H_hh_bbbb_e
double ip_ex_pp_H_hh_bblnulnu_CMS13_cache_e[2][CacheSize]
gslpp::complex A_A_L(const double mA2, const double cW2, const double Mmu, const double Mtau, const double MZ) const
Amplitude for a CP-odd Higgs boson decay to a photon and a Z boson including muons and taus in the lo...
double THoEX_VBF_H_ZZ_llll_ATLAS13
double ip_ex_gg_phi_tt_ATLAS8_cache[2][CacheSize]
double ip_ex_pp_phi_bb_CMS13_cache[2][CacheSize]
gslpp::complex B0p_MZ2_mHl2_mHp2_mHp2(const double MZ2, const double mHl2, const double mHp2) const
.
double THoEX_pp_H_hh_bbgaga_CMS13
double THoEX_pp_Hpm_taunu_CMS13
gslpp::matrix< double > CMS13_pp_phi_bb_e
double ie12e(double mass)
double R_pp_H_Zga_llga_CMS13
gslpp::matrix< double > CMS13_pp_H_hh_bbtautau_e
gslpp::complex B00_MZ2_0_mHl2_mA2(const double MZ2, const double mHl2, const double mA2) const
.
double ip_ex_gg_H_ZZ_llnunu_ATLAS13_cache_e[2][CacheSize]
double ie39c[2][CacheSize]
double R_pp_H_AZ_tautaull_CMS8
gslpp::matrix< double > csrA_bottom_8
double R_pp_Hp_tb_ATLAS13_1
double THoEX_ggF_H_tt_ATLAS8
double ip_ex_gg_H_WW_enumunu_ATLAS13_cache[2][CacheSize]
gslpp::complex B0_MZ2_0_MZ2_mHl2(const double MZ2, const double mHl2) const
.
double ip_ex_gg_H_WW_enumunu_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_pp_H_hh_bbtautau1_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
double ip_ex_bb_phi_bb_CMS8_cache[2][CacheSize]
gslpp::matrix< double > csrH_bottom_13
double ip_ex_gg_H_ZZ_llnunu_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced scalar resonance de...
gslpp::complex B00_MZ2_MW2_mHh2_mHp2_cache[5][CacheSize]
double ip_Br_HPtoWW(double mass)
Interpolating function for the SM branching ratio to two bosons.
gslpp::matrix< double > CMS8_pp_A_Zga_llga_e
double ip_ex_pp_H_hh_bblnulnu_CMS13_cache[2][CacheSize]
double ip_cs_ggtoHp_8_cache[3][CacheSize]
double THoEX_pp_Hp_tb_ATLAS13_1
gslpp::matrix< double > csrH_top_8
gslpp::complex B0_MZ2_mA2_mA2_mHh2_cache[4][CacheSize]
double THoEX_ggF_H_ZZ_llqq_ATLAS13
double ip_csr_ggA_b_13_cache[2][CacheSize]
double ip_ex_pp_Hpm_taunu_CMS13_cache[2][CacheSize]
double ip_cs_pptottA_13_cache[2][CacheSize]
double ip_ex_gg_phi_Zga_llga_ATLAS13_cache[2][CacheSize]
gslpp::complex B0p_MZ2_0_mA2_mHl2(const double MZ2, const double mA2, const double mHl2) const
.
double ip_ex_pp_phi_Zga_llga_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a pseudoscalar resonance decaying to a Z...
gslpp::matrix< double > temp22e
double ip_ex_pp_phi_bb_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to a b quark p...
double ip_ex_VBF_H_ZZ_llllnunu_ATLAS13_cache_e[2][CacheSize]
gslpp::complex B0_MZ2_mHh2_0_mA2_cache[4][CacheSize]
gslpp::matrix< double > br_cc
double ip_cs_ggtoH_8(double mass)
Interpolating function for the H production cross section via gluon-gluon fusion at 8 TeV.
double THoEX_bbF_H_tautau_ATLAS13
gslpp::matrix< double > temp12e
double THoEX_VBF_H_WW_enumunu_ATLAS13
gslpp::matrix< double > temp23e
gslpp::complex B0p_MZ2_mHp2_mHp2_mHl2(const double MZ2, const double mHp2, const double mHl2) const
.
gslpp::matrix< double > ATLAS13_gg_H_ZZ_llll_e
double ip_ex_VBF_H_ZZ_qqllnunu_ATLAS13_cache[2][CacheSize]
double ip_ex_VBF_H_ZZ_llll_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > ATLAS8_gg_phi_tt_e
double ip_ex_gg_phi_Zga_llga_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::complex B0p_MZ2_mHl2_mHh2_mHl2(const double MZ2, const double mHl2, const double mHh2) const
.
gslpp::matrix< double > ATLAS13_VBF_H_ZZ_llqq
gslpp::matrix< double > temp28e
double ip_cs_ZtoZH_13(double mass)
Interpolating function for the Z associated H production cross section at 13 TeV.
gslpp::matrix< double > ATLAS8_pp_phi_Zga_llga
double ip_ex_gg_phi_tautau_CMS8_cache[2][CacheSize]
double ip_ex_pp_phi_gaga_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a scalar resonance decaying to two photo...
double ip_ex_gg_A_hZ_bbZ_ATLAS8_cache_e[2][CacheSize]
double ip_ex_bb_phi_tautau_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a b associated scalar resonance decaying t...
gslpp::complex B0_MZ2_mHh2_mHl2_mHl2_cache[4][CacheSize]
gslpp::complex B0p_MZ2_0_mHp2_mHl2_cache[4][CacheSize]
gslpp::complex I_h_D(const double mHl2, const double Md, const double Ms, const double Mb) const
Amplitude for the SM Higgs boson decay to diphotons including the down-type quarks in the loop.
gslpp::matrix< double > CMS13_pp_H_hh_gagabb_e
double R_ttF_A_tt_ATLAS13
gslpp::complex B00_MZ2_0_mA2_mHp2(const double MZ2, const double mA2, const double mHp2) const
.
gslpp::complex B0p_MZ2_mHl2_mHp2_mHp2_cache[4][CacheSize]
gslpp::matrix< double > CMS8_gg_phi_tautau
double ip_ex_gg_phi_tautau_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::matrix< double > temp6
double ip_ex_VBF_H_WW_lnuqq_ATLAS13_cache_e[2][CacheSize]
double ie40e(double mass)
double THoEX_ggF_A_hZ_tautaull_CMS8
gslpp::matrix< double > temp32e
double ip_ex_pp_H_ZZ_llnunu_CMS13_cache[2][CacheSize]
double ip_ex_tt_phi_tt_ATLAS13_cache_e[2][CacheSize]
gslpp::matrix< double > ATLAS8_VBF_H_WW
double ip_cs_ggtoH_13(double mass)
Interpolating function for the H production cross section via gluon-gluon fusion at 13 TeV.
gslpp::matrix< double > CMS13_pp_H_hh_bblnulnu_e
double ip_ex_pp_phi_Zga_llga_ATLAS13_cache[2][CacheSize]
double ip_ex_pp_phi_gaga_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a scalar resonance decaying to two photo...
gslpp::matrix< double > temp26e
double ip_ex_bb_phi_bb_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a bottom quark produced scalar resonance d...
double ip_ex_ggF_phi_Zga_CMS13_cache[2][CacheSize]
gslpp::complex B0_MZ2_MZ2_MZ2_mHl2(const double MZ2, const double mHl2) const
.
double ip_ex_mu_pp_H_VV_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two massive...
double ip_ex_VBF_H_WW_enumunu_ATLAS13_cache[2][CacheSize]
double ie21c[2][CacheSize]
double R_pp_A_gaga_ATLAS13
gslpp::matrix< double > temp29e
double THoEX_pp_H_gaga_ATLAS8
double ip_ex_bb_phi_tautau_ATLAS8_cache[2][CacheSize]
gslpp::complex B0_MZ2_mHl2_mHh2_mHl2(const double MZ2, const double mHl2, const double mHh2) const
.
gslpp::matrix< double > CMS13_pp_Hpm_taunu_e
double ie29e(double mass)
double ip_ex_gg_phi_tautau_ATLAS8_cache[2][CacheSize]
gslpp::matrix< double > CMS13_pp_H_ZZ_llnunu_e
double ip_ex_gg_H_ZZ_llqq_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > temp35
gslpp::matrix< double > temp16e
gslpp::complex I_HH_L(const double mHh2, const double Mmu, const double Mtau) const
Amplitude for a heavy CP-even Higgs boson decay to diphotons including muons and taus in the loop.
double ie35e(double mass)
gslpp::matrix< double > CMS8_pp_A_HZ_bbll
double ip_ex_pp_Hpm_taunu_ATLAS13_cache[2][CacheSize]
double R_ggF_H_hh_bbtautau_CMS8
gslpp::complex B0_MZ2_mHh2_mHh2_mHl2(const double MZ2, const double mHh2, const double mHl2) const
.
gslpp::complex B0p_MZ2_mA2_mA2_mHl2_cache[4][CacheSize]
double R_pp_A_HZ_tautaull_CMS8
gslpp::complex B0p_MZ2_mHp2_0_mHh2_cache[4][CacheSize]
double THoEX_pp_A_HZ_tautaull_CMS8
double R_pp_H_hh_bbbb_CMS8
double ip_ex_gg_H_ZZ_llnunu_ATLAS13_cache[2][CacheSize]
double THoEX_pp_H_ZZ_llqq_CMS13
gslpp::matrix< double > temp31e
gslpp::complex B0p_MZ2_mA2_0_mHh2_cache[4][CacheSize]
double R_pp_H_hh_bbbb_CMS13
double THoEX_pp_A_gaga_ATLAS13
gslpp::complex B0_MZ2_mHl2_0_mHp2(const double MZ2, const double mHl2, const double mHp2) const
.
gslpp::matrix< double > ATLAS13_gg_A_Zh_Zbb_e
double ip_ex_gg_A_hZ_tautauZ_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced pseudoscalar reso...
double ip_cs_pptottH_13(double mass)
Interpolating function for the top associated H production cross section at 13 TeV.
double ip_ex_pp_Hp_tb_ATLAS13_1_e(double mass)
Interpolating function for the expected ATLAS upper limit on a singly charged scalar resonance decayi...
gslpp::matrix< double > temp2e
double R_VBF_H_ZZ_llll_ATLAS13
double R_VBF_H_WW_lnuqq_ATLAS13
double THoEX_pp_H_AZ_tautaull_CMS8
double ie21ec[2][CacheSize]
double ie18c[2][CacheSize]
double ie35c[2][CacheSize]
double THoEX_pp_H_Zga_qqga_CMS13
double ip_ex_gg_H_ZZ_llll_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::complex B00_MZ2_MZ2_mHp2_mHp2(const double MZ2, const double mHp2) const
.
gslpp::complex B0p_MZ2_mHp2_0_mHl2_cache[4][CacheSize]
gslpp::matrix< double > CMS13_gg_H_ZZ_llnunu_e
double THoEX_ggF_H_WW_ATLAS8
double ip_ex_pp_H_AZ_tautaull_CMS8_cache[3][CacheSize]
double ip_ex_gg_H_ZZ_qqllnunu_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double R_ggF_H_ZZ_llnunu_CMS13
double ie24ec[2][CacheSize]
gslpp::matrix< double > ATLAS13_VBF_H_WW_lnuqq
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 ip_ex_pp_H_hh_gagabb_ATLAS13_cache[2][CacheSize]
gslpp::complex A_A_U(const double mA2, const double cW2, const double Mc, const double Mt, const double MZ) const
Amplitude for a CP-odd Higgs boson decay to a photon and a Z boson including the charm and top quarks...
double ip_ex_bb_phi_bb_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a bottom quark produced scalar resonance d...
double ip_ex_gg_H_ZZ_llll_ATLAS13_cache[2][CacheSize]
gslpp::complex B0p_MZ2_mHh2_mHp2_mHp2(const double MZ2, const double mHh2, const double mHp2) const
.
gslpp::complex B0_MZ2_mA2_0_mHh2(const double MZ2, const double mA2, const double mHh2) const
.
double THoEX_VBF_H_WW_ATLAS8
double ip_cs_ZtoZH_8_cache[2][CacheSize]
double ie32e(double mass)
double ip_ex_gg_phi_tautau_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_gg_H_hh_bbtautau_CMS8_cache_e[2][CacheSize]
double ip_ex_pp_phi_Zga_llga_ATLAS13_cache_e[2][CacheSize]
double ip_ex_gg_H_ZZ_nunuqq_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > temp38e
double ip_ex_pp_H_hh_bbbb_CMS13_cache[2][CacheSize]
double THoEX_pp_A_gaga_ATLAS8
double ip_ex_pp_Hpm_tb_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a singly charged scalar resonance decayi...
gslpp::matrix< double > log_cs_ttA_13
double ip_ex_gg_phi_tautau_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced scalar resonance de...
gslpp::matrix< double > CMS13_gg_phi_gaga_e
double ip_ex_pp_H_hh_bbbb_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double ghHpHm_cache[7][CacheSize]
double ip_cs_WtoWH_8(double mass)
Interpolating function for the W associated H production cross section at 8 TeV.
double ip_Br_HPtobb(double mass)
Interpolating function for the SM branching ratio to two bottom quarks.
parameter of the Higgs potential
double computeThValue()
The quartic coupling .
An observable class for the quartic Higgs potential coupling .
double computeThValue()
The quartic coupling .
An observable class for the quartic Higgs potential coupling .
double computeThValue()
The quartic coupling .
An observable class for the quartic Higgs potential coupling .
double computeThValue()
The quartic coupling .
An observable class for the quartic Higgs potential coupling .
double computeThValue()
The quartic coupling .
An observable class for the quadratic Higgs potential coupling .
double computeThValue()
The quartic coupling .
An observable class for the quadratic Higgs potential coupling .
double computeThValue()
The quartic coupling .