19 : br_tt(19961, 2, 0.),
21 br_tautau(19961, 2, 0.),
23 br_mumu(19961, 2, 0.),
26 GammaHtot_SM(19961, 2, 0.),
27 log_cs_ggH_8(199, 2, 0.),
28 log_cs_VBF_8(199, 2, 0.),
29 log_cs_WH_8(199, 2, 0.),
30 log_cs_ZH_8(199, 2, 0.),
31 log_cs_ggH_13(199, 2, 0.),
32 log_cs_VBF_13(199, 2, 0.),
33 log_cs_WH_13(199, 2, 0.),
34 log_cs_ZH_13(199, 2, 0.),
35 log_cs_ttH_8(199, 2, 0.),
36 log_cs_ttH_13(199, 2, 0.),
37 log_cs_bbH_8(199, 2, 0.),
38 log_cs_bbH_13(199, 2, 0.),
39 log_cs_ggA_8(199, 2, 0.),
40 log_cs_ttA_8(199, 2, 0.),
41 log_cs_bbA_8(199, 2, 0.),
42 log_cs_ggA_13(199, 2, 0.),
43 log_cs_ttA_13(199, 2, 0.),
44 log_cs_bbA_13(199, 2, 0.),
45 log_cs_ggHp_8(744, 3, 0.),
46 log_cs_ggHp_13(1104, 3, 0.),
47 csrH_top_8(199, 2, 0.),
48 csrH_bottom_8(199, 2, 0.),
49 csrA_top_8(199, 2, 0.),
50 csrA_bottom_8(199, 2, 0.),
51 csrH_top_13(199, 2, 0.),
52 csrH_bottom_13(199, 2, 0.),
53 csrA_top_13(199, 2, 0.),
54 csrA_bottom_13(199, 2, 0.),
56 ATLAS8_gg_phi_tautau(92, 2, 0.),
57 ATLAS8_bb_phi_tautau(92, 2, 0.),
58 ATLAS8_gg_phi_gaga(108, 2, 0.),
59 ATLAS8_pp_phi_Zga_llga(141, 2, 0.),
60 ATLAS8_gg_phi_ZZ(173, 2, 0.),
61 ATLAS8_VV_phi_ZZ(173, 2, 0.),
62 ATLAS8_gg_phi_WW(13, 2, 0.),
63 ATLAS8_VV_phi_WW(13, 2, 0.),
64 ATLAS8_gg_phi_phi1phi1(75, 2, 0.),
65 ATLAS8_gg_phi_phi1Z_bbZ(79, 2, 0.),
66 ATLAS8_gg_phi_phi1Z_tautauZ(79, 2, 0.),
67 CMS8_mu_pp_phi_VV(172, 2, 0.),
68 CMS8_bb_phi_bb(81, 2, 0.),
69 CMS8_gg_phi_bb(88, 2, 0.),
70 CMS8_gg_phi_tautau(92, 2, 0.),
71 CMS8_bb_phi_tautau(92, 2, 0.),
72 CMS8_pp_phi_Zga_llga(101, 2, 0.),
73 CMS8_pp_phi_phi1phi1_bbbb(167, 2, 0.),
74 CMS8_pp_phi_phi1phi1_bbgaga(85, 2, 0.),
75 CMS8_gg_phi_phi1phi1_bbtautau(10, 2, 0.),
76 CMS8_pp_phi_phi1phi1_bbtautau(71, 2, 0.),
77 CMS8_gg_phi_phi1Z_bbll(16, 2, 0.),
78 CMS8_gg_phi_phi1Z_tautaull(14, 2, 0.),
79 CMS8_pp_phii_phijZ_bbll_1(28718, 3, 0.),
80 CMS8_pp_phii_phijZ_bbll_2(29050, 3, 0.),
81 CMS8_pp_phii_phijZ_tautaull_1(400, 3, 0.),
82 CMS8_pp_phii_phijZ_tautaull_2(400, 3, 0.),
83 ATLAS13_tt_phi_tt(61, 2, 0.),
84 ATLAS13_bb_phi_tt(61, 2, 0.),
85 ATLAS13_gg_phi_tautau(206, 2, 0.),
86 ATLAS13_bb_phi_tautau(206, 2, 0.),
87 ATLAS13_pp_phi_gaga(251, 2, 0.),
88 ATLAS13_gg_phi_Zga_llga(216, 2, 0.),
89 ATLAS13_gg_phi_Zga_qqga(581, 2, 0.),
90 ATLAS13_gg_phi_ZZ_llllnunu(101, 2, 0.),
91 ATLAS13_VV_phi_ZZ_llllnunu(101, 2, 0.),
92 ATLAS13_gg_phi_ZZ_qqllnunu(271, 2, 0.),
93 ATLAS13_VV_phi_ZZ_qqllnunu(271, 2, 0.),
94 ATLAS13_gg_phi_WW_enumunu(381, 2, 0.),
95 ATLAS13_VV_phi_WW_enumunu(281, 2, 0.),
96 ATLAS13_gg_phi_WW_lnuqq(271, 2, 0.),
97 ATLAS13_VV_phi_WW_lnuqq(271, 2, 0.),
98 ATLAS13_pp_phi_VV_qqqq(181, 2, 0.),
99 ATLAS13_pp_phi_phi1phi1_bbbb(275, 2, 0.),
100 ATLAS13_pp_phi_phi1phi1_bbgaga(75, 2, 0.),
101 ATLAS13_pp_phi_phi1phi1_bbtautau(75, 2, 0.),
102 ATLAS13_pp_phi_phi1phi1_bbWW(51, 2, 0.),
103 ATLAS13_gg_phi_phi1phi1_gagaWW(25, 2, 0.),
104 ATLAS13_gg_phi_phi1Z_bbZ(181, 2, 0.),
105 ATLAS13_bb_phi_phi1Z_bbZ(181, 2, 0.),
106 ATLAS13_gg_phii_phijZ_bbZ(3364, 3, 0.),
107 ATLAS13_bb_phii_phijZ_bbZ(3364, 3, 0.),
108 CMS13_pp_phi_bb(66, 2, 0.),
109 CMS13_bb_phi_bb(101, 2, 0.),
110 CMS13_gg_phi_tautau(312, 2, 0.),
111 CMS13_bb_phi_tautau(312, 2, 0.),
112 CMS13_gg_phi_gaga(351, 2, 0.),
113 CMS13_gg_phi_Zga(366, 2, 0.),
114 CMS13_pp_phi_ZZ_llqqnunull(288, 2, 0.),
115 CMS13_pp_phi_ZZ_qqnunu(301, 2, 0.),
116 CMS13_ggVV_phi_WW_lnulnu(81, 2, 0.),
117 CMS13_pp_phi_WW_lnuqq(341, 2, 0.),
118 CMS13_pp_phi_phi1phi1_bbbb_1(95, 2, 0.),
119 CMS13_pp_phi_phi1phi1_bbbb_2(181, 2, 0.),
120 CMS13_pp_phi_phi1phi1_bbgaga(66, 2, 0.),
121 CMS13_pp_phi_phi1phi1_bbtautau_1(66, 2, 0.),
122 CMS13_pp_phi_phi1phi1_bbtautau_2(311, 2, 0.),
123 CMS13_pp_phi_phi1phi1_bbVV(65, 2, 0.),
124 CMS13_gg_phi_phi1Z_bbZ_1(79, 2, 0.),
125 CMS13_gg_phi_phi1Z_bbZ_2(121, 2, 0.),
126 CMS13_bb_phi_phi1Z_bbZ_1(79, 2, 0.),
127 CMS13_bb_phi_phi1Z_bbZ_2(121, 2, 0.),
128 ATLAS8_pp_Hpm_taunu(83, 2, 0.),
129 ATLAS8_pp_Hpm_tb(41, 2, 0.),
130 CMS8_pp_Hp_taunu(43, 2, 0.),
131 CMS8_pp_Hp_tb(43, 2, 0.),
132 ATLAS13_pp_Hpm_taunu(192, 2, 0.),
133 ATLAS13_pp_Hpm_tb(181, 2, 0.),
134 CMS13_pp_Hpm_taunu(283, 2, 0.),
135 arraybsgamma(1111, 3, 0.),
136 Mu_GTHDM(3,3,0.), Md_GTHDM(3,3,0.), Ml_GTHDM(3,3,0.),
137 Nu_GTHDM(3,3,0.), Nd_GTHDM(3,3,0.), Nl_GTHDM(3,3,0.),
138 Yu1_GTHDM(3,3,0.), Yu2_GTHDM(3,3,0.), Yd1_GTHDM(3,3,0.), Yd2_GTHDM(3,3,0.),
139 Yl1_GTHDM(3,3,0.), Yl2_GTHDM(3,3,0.),
152 const int NumPar,
const double params[])
const {
156 for(
int j=0; j<NumPar; j++)
157 bCache &= (params[j] == cache[j][i].real());
158 if (bCache)
return i;
164 const int NumPar,
const double params[])
const {
168 for(
int j=0; j<NumPar; j++)
169 bCache &= (params[j] == cache[j][i]);
170 if (bCache)
return i;
180 for(
int j=0; j<NumPar+1; j++)
181 cache[j][i] = cache[j][i-1];
184 for(
int j=0; j<NumPar; j++) {
186 cache[NumPar][0] = newResult;
191 const double params[],
const double newResult)
const {
194 for(
int j=0; j<NumPar+1; j++)
195 cache[j][i] = cache[j][i-1];
198 for(
int j=0; j<NumPar; j++) {
199 cache[j][0] = params[j];
200 cache[NumPar][0] = newResult;
210 double params[] = {MZ2, MW2, mHh2};
212 int i =
CacheCheck(B0_MZ2_0_MW2_mHh2_cache, NumPar, params);
214 return ( B0_MZ2_0_MW2_mHh2_cache[NumPar][i] );
217 CacheShift(B0_MZ2_0_MW2_mHh2_cache, NumPar, params, newResult);
224 double params[] = {MZ2, MW2, mHl2};
226 int i =
CacheCheck(B0_MZ2_0_MW2_mHl2_cache, NumPar, params);
228 return ( B0_MZ2_0_MW2_mHl2_cache[NumPar][i] );
231 CacheShift(B0_MZ2_0_MW2_mHl2_cache, NumPar, params, newResult);
238 double params[] = {MZ2, mHh2};
240 int i =
CacheCheck(B0_MZ2_0_MZ2_mHh2_cache, NumPar, params);
242 return ( B0_MZ2_0_MZ2_mHh2_cache[NumPar][i] );
245 CacheShift(B0_MZ2_0_MZ2_mHh2_cache, NumPar, params, newResult);
252 double params[] = {MZ2, mHl2};
254 int i =
CacheCheck(B0_MZ2_0_MZ2_mHl2_cache, NumPar, params);
256 return ( B0_MZ2_0_MZ2_mHl2_cache[NumPar][i] );
259 CacheShift(B0_MZ2_0_MZ2_mHl2_cache, NumPar, params, newResult);
266 double params[] = {MZ2, MW2, mHh2};
268 int i =
CacheCheck(B0_MZ2_MW2_MW2_mHh2_cache, NumPar, params);
270 return ( B0_MZ2_MW2_MW2_mHh2_cache[NumPar][i] );
273 CacheShift(B0_MZ2_MW2_MW2_mHh2_cache, NumPar, params, newResult);
280 double params[] = {MZ2, MW2, mHl2};
282 int i =
CacheCheck(B0_MZ2_MW2_MW2_mHl2_cache, NumPar, params);
284 return ( B0_MZ2_MW2_MW2_mHl2_cache[NumPar][i] );
287 CacheShift(B0_MZ2_MW2_MW2_mHl2_cache, NumPar, params, newResult);
294 double params[] = {MZ2, mHh2};
296 int i =
CacheCheck(B0_MZ2_MZ2_MZ2_mHh2_cache, NumPar, params);
298 return ( B0_MZ2_MZ2_MZ2_mHh2_cache[NumPar][i] );
301 CacheShift(B0_MZ2_MZ2_MZ2_mHh2_cache, NumPar, params, newResult);
308 double params[] = {MZ2, mHl2};
310 int i =
CacheCheck(B0_MZ2_MZ2_MZ2_mHl2_cache, NumPar, params);
312 return ( B0_MZ2_MZ2_MZ2_mHl2_cache[NumPar][i] );
315 CacheShift(B0_MZ2_MZ2_MZ2_mHl2_cache, NumPar, params, newResult);
322 double params[] = {MZ2, mHl2};
324 int i =
CacheCheck(B0_MZ2_0_0_mHl2_cache, NumPar, params);
326 return ( B0_MZ2_0_0_mHl2_cache[NumPar][i] );
329 CacheShift(B0_MZ2_0_0_mHl2_cache, NumPar, params, newResult);
336 double params[] = {MZ2, mHh2};
338 int i =
CacheCheck(B0_MZ2_0_0_mHh2_cache, NumPar, params);
340 return ( B0_MZ2_0_0_mHh2_cache[NumPar][i] );
343 CacheShift(B0_MZ2_0_0_mHh2_cache, NumPar, params, newResult);
350 double params[] = {MZ2, mHp2, mHl2};
352 int i =
CacheCheck(B0_MZ2_0_mHp2_mHl2_cache, NumPar, params);
354 return ( B0_MZ2_0_mHp2_mHl2_cache[NumPar][i] );
357 CacheShift(B0_MZ2_0_mHp2_mHl2_cache, NumPar, params, newResult);
364 double params[] = {MZ2, mHp2, mHh2};
366 int i =
CacheCheck(B0_MZ2_0_mHp2_mHh2_cache, NumPar, params);
368 return ( B0_MZ2_0_mHp2_mHh2_cache[NumPar][i] );
371 CacheShift(B0_MZ2_0_mHp2_mHh2_cache, NumPar, params, newResult);
378 double params[] = {MZ2, mA2, mHl2};
380 int i =
CacheCheck(B0_MZ2_0_mA2_mHl2_cache, NumPar, params);
382 return ( B0_MZ2_0_mA2_mHl2_cache[NumPar][i] );
385 CacheShift(B0_MZ2_0_mA2_mHl2_cache, NumPar, params, newResult);
392 double params[] = {MZ2, mA2, mHh2};
394 int i =
CacheCheck(B0_MZ2_0_mA2_mHh2_cache, NumPar, params);
396 return ( B0_MZ2_0_mA2_mHh2_cache[NumPar][i] );
399 CacheShift(B0_MZ2_0_mA2_mHh2_cache, NumPar, params, newResult);
406 double params[] = {MZ2, mHl2};
408 int i =
CacheCheck(B0_MZ2_mHl2_0_0_cache, NumPar, params);
410 return ( B0_MZ2_mHl2_0_0_cache[NumPar][i] );
413 CacheShift(B0_MZ2_mHl2_0_0_cache, NumPar, params, newResult);
420 double params[] = {MZ2, mHl2, mHp2};
422 int i =
CacheCheck(B0_MZ2_mHl2_0_mHp2_cache, NumPar, params);
424 return ( B0_MZ2_mHl2_0_mHp2_cache[NumPar][i] );
427 CacheShift(B0_MZ2_mHl2_0_mHp2_cache, NumPar, params, newResult);
434 double params[] = {MZ2, mHl2, mA2};
436 int i =
CacheCheck(B0_MZ2_mHl2_0_mA2_cache, NumPar, params);
438 return ( B0_MZ2_mHl2_0_mA2_cache[NumPar][i] );
441 CacheShift(B0_MZ2_mHl2_0_mA2_cache, NumPar, params, newResult);
448 double params[] = {MZ2, mHl2};
450 int i =
CacheCheck(B0_MZ2_mHl2_mHl2_mHl2_cache, NumPar, params);
452 return ( B0_MZ2_mHl2_mHl2_mHl2_cache[NumPar][i] );
455 CacheShift(B0_MZ2_mHl2_mHl2_mHl2_cache, NumPar, params, newResult);
462 double params[] = {MZ2, mHl2, mHh2};
464 int i =
CacheCheck(B0_MZ2_mHl2_mHh2_mHl2_cache, NumPar, params);
466 return ( B0_MZ2_mHl2_mHh2_mHl2_cache[NumPar][i] );
469 CacheShift(B0_MZ2_mHl2_mHh2_mHl2_cache, NumPar, params, newResult);
476 double params[] = {MZ2, mHl2, mHh2};
478 int i =
CacheCheck(B0_MZ2_mHl2_mHh2_mHh2_cache, NumPar, params);
480 return ( B0_MZ2_mHl2_mHh2_mHh2_cache[NumPar][i] );
483 CacheShift(B0_MZ2_mHl2_mHh2_mHh2_cache, NumPar, params, newResult);
490 double params[] = {MZ2, mHl2, mHp2};
492 int i =
CacheCheck(B0_MZ2_mHl2_mHp2_mHp2_cache, NumPar, params);
494 return ( B0_MZ2_mHl2_mHp2_mHp2_cache[NumPar][i] );
497 CacheShift(B0_MZ2_mHl2_mHp2_mHp2_cache, NumPar, params, newResult);
504 double params[] = {MZ2, mHl2, mA2};
506 int i =
CacheCheck(B0_MZ2_mHl2_mA2_mA2_cache, NumPar, params);
508 return ( B0_MZ2_mHl2_mA2_mA2_cache[NumPar][i] );
511 CacheShift(B0_MZ2_mHl2_mA2_mA2_cache, NumPar, params, newResult);
518 double params[] = {MZ2, mHh2};
520 int i =
CacheCheck(B0_MZ2_mHh2_0_0_cache, NumPar, params);
522 return ( B0_MZ2_mHh2_0_0_cache[NumPar][i] );
525 CacheShift(B0_MZ2_mHh2_0_0_cache, NumPar, params, newResult);
532 double params[] = {MZ2, mHh2, mHp2};
534 int i =
CacheCheck(B0_MZ2_mHh2_0_mHp2_cache, NumPar, params);
536 return ( B0_MZ2_mHh2_0_mHp2_cache[NumPar][i] );
539 CacheShift(B0_MZ2_mHh2_0_mHp2_cache, NumPar, params, newResult);
546 double params[] = {MZ2, mHh2, mA2};
548 int i =
CacheCheck(B0_MZ2_mHh2_0_mA2_cache, NumPar, params);
550 return ( B0_MZ2_mHh2_0_mA2_cache[NumPar][i] );
553 CacheShift(B0_MZ2_mHh2_0_mA2_cache, NumPar, params, newResult);
560 double params[] = {MZ2, mHh2, mHl2};
562 int i =
CacheCheck(B0_MZ2_mHh2_mHl2_mHl2_cache, NumPar, params);
564 return ( B0_MZ2_mHh2_mHl2_mHl2_cache[NumPar][i] );
567 CacheShift(B0_MZ2_mHh2_mHl2_mHl2_cache, NumPar, params, newResult);
574 double params[] = {MZ2, mHh2, mHl2};
576 int i =
CacheCheck(B0_MZ2_mHh2_mHh2_mHl2_cache, NumPar, params);
578 return ( B0_MZ2_mHh2_mHh2_mHl2_cache[NumPar][i] );
581 CacheShift(B0_MZ2_mHh2_mHh2_mHl2_cache, NumPar, params, newResult);
588 double params[] = {MZ2, mHh2};
590 int i =
CacheCheck(B0_MZ2_mHh2_mHh2_mHh2_cache, NumPar, params);
592 return ( B0_MZ2_mHh2_mHh2_mHh2_cache[NumPar][i] );
595 CacheShift(B0_MZ2_mHh2_mHh2_mHh2_cache, NumPar, params, newResult);
602 double params[] = {MZ2, mHh2, mHp2};
604 int i =
CacheCheck(B0_MZ2_mHh2_mHp2_mHp2_cache, NumPar, params);
606 return ( B0_MZ2_mHh2_mHp2_mHp2_cache[NumPar][i] );
609 CacheShift(B0_MZ2_mHh2_mHp2_mHp2_cache, NumPar, params, newResult);
616 double params[] = {MZ2, mHh2, mA2};
618 int i =
CacheCheck(B0_MZ2_mHh2_mA2_mA2_cache, NumPar, params);
620 return ( B0_MZ2_mHh2_mA2_mA2_cache[NumPar][i] );
623 CacheShift(B0_MZ2_mHh2_mA2_mA2_cache, NumPar, params, newResult);
630 double params[] = {MZ2, mHp2, mHl2};
632 int i =
CacheCheck(B0_MZ2_mHp2_0_mHl2_cache, NumPar, params);
634 return ( B0_MZ2_mHp2_0_mHl2_cache[NumPar][i] );
637 CacheShift(B0_MZ2_mHp2_0_mHl2_cache, NumPar, params, newResult);
644 double params[] = {MZ2, mHp2, mHh2};
646 int i =
CacheCheck(B0_MZ2_mHp2_0_mHh2_cache, NumPar, params);
648 return ( B0_MZ2_mHp2_0_mHh2_cache[NumPar][i] );
651 CacheShift(B0_MZ2_mHp2_0_mHh2_cache, NumPar, params, newResult);
658 double params[] = {MZ2, mHp2, mHl2};
660 int i =
CacheCheck(B0_MZ2_mHp2_mHp2_mHl2_cache, NumPar, params);
662 return ( B0_MZ2_mHp2_mHp2_mHl2_cache[NumPar][i] );
665 CacheShift(B0_MZ2_mHp2_mHp2_mHl2_cache, NumPar, params, newResult);
672 double params[] = {MZ2, mHp2, mHh2};
674 int i =
CacheCheck(B0_MZ2_mHp2_mHp2_mHh2_cache, NumPar, params);
676 return ( B0_MZ2_mHp2_mHp2_mHh2_cache[NumPar][i] );
679 CacheShift(B0_MZ2_mHp2_mHp2_mHh2_cache, NumPar, params, newResult);
686 double params[] = {MZ2, mA2, mHl2};
688 int i =
CacheCheck(B0_MZ2_mA2_0_mHl2_cache, NumPar, params);
690 return ( B0_MZ2_mA2_0_mHl2_cache[NumPar][i] );
693 CacheShift(B0_MZ2_mA2_0_mHl2_cache, NumPar, params, newResult);
700 double params[] = {MZ2, mA2, mHh2};
702 int i =
CacheCheck(B0_MZ2_mA2_0_mHh2_cache, NumPar, params);
704 return ( B0_MZ2_mA2_0_mHh2_cache[NumPar][i] );
707 CacheShift(B0_MZ2_mA2_0_mHh2_cache, NumPar, params, newResult);
714 double params[] = {MZ2, mA2, mHl2};
716 int i =
CacheCheck(B0_MZ2_mA2_mA2_mHl2_cache, NumPar, params);
718 return ( B0_MZ2_mA2_mA2_mHl2_cache[NumPar][i] );
721 CacheShift(B0_MZ2_mA2_mA2_mHl2_cache, NumPar, params, newResult);
728 double params[] = {MZ2, mA2, mHh2};
730 int i =
CacheCheck(B0_MZ2_mA2_mA2_mHh2_cache, NumPar, params);
732 return ( B0_MZ2_mA2_mA2_mHh2_cache[NumPar][i] );
735 CacheShift(B0_MZ2_mA2_mA2_mHh2_cache, NumPar, params, newResult);
746 double params[] = {MZ2, mHl2};
748 int i =
CacheCheck(B0p_MZ2_0_0_mHl2_cache, NumPar, params);
750 return ( B0p_MZ2_0_0_mHl2_cache[NumPar][i] );
753 CacheShift(B0p_MZ2_0_0_mHl2_cache, NumPar, params, newResult);
760 double params[] = {MZ2, mHh2};
762 int i =
CacheCheck(B0p_MZ2_0_0_mHh2_cache, NumPar, params);
764 return ( B0p_MZ2_0_0_mHh2_cache[NumPar][i] );
767 CacheShift(B0p_MZ2_0_0_mHh2_cache, NumPar, params, newResult);
774 double params[] = {MZ2, mHp2, mHl2};
776 int i =
CacheCheck(B0p_MZ2_0_mHp2_mHl2_cache, NumPar, params);
778 return ( B0p_MZ2_0_mHp2_mHl2_cache[NumPar][i] );
781 CacheShift(B0p_MZ2_0_mHp2_mHl2_cache, NumPar, params, newResult);
788 double params[] = {MZ2, mHp2, mHh2};
790 int i =
CacheCheck(B0p_MZ2_0_mHp2_mHh2_cache, NumPar, params);
792 return ( B0p_MZ2_0_mHp2_mHh2_cache[NumPar][i] );
795 CacheShift(B0p_MZ2_0_mHp2_mHh2_cache, NumPar, params, newResult);
802 double params[] = {MZ2, mHp2, mA2};
804 int i =
CacheCheck(B0p_MZ2_0_mHp2_mA2_cache, NumPar, params);
806 return ( B0p_MZ2_0_mHp2_mA2_cache[NumPar][i] );
809 CacheShift(B0p_MZ2_0_mHp2_mA2_cache, NumPar, params, newResult);
816 double params[] = {MZ2, mA2, mHl2};
818 int i =
CacheCheck(B0p_MZ2_0_mA2_mHl2_cache, NumPar, params);
820 return ( B0p_MZ2_0_mA2_mHl2_cache[NumPar][i] );
823 CacheShift(B0p_MZ2_0_mA2_mHl2_cache, NumPar, params, newResult);
830 double params[] = {MZ2, mA2, mHh2};
832 int i =
CacheCheck(B0p_MZ2_0_mA2_mHh2_cache, NumPar, params);
834 return ( B0p_MZ2_0_mA2_mHh2_cache[NumPar][i] );
837 CacheShift(B0p_MZ2_0_mA2_mHh2_cache, NumPar, params, newResult);
844 double params[] = {MZ2, mHl2};
846 int i =
CacheCheck(B0p_MZ2_0_0_mHl2_cache, NumPar, params);
848 return ( B0p_MZ2_0_0_mHl2_cache[NumPar][i] );
851 CacheShift(B0p_MZ2_0_0_mHl2_cache, NumPar, params, newResult);
858 double params[] = {MZ2, mHl2, mHp2};
860 int i =
CacheCheck(B0p_MZ2_mHl2_0_mHp2_cache, NumPar, params);
862 return ( B0p_MZ2_mHl2_0_mHp2_cache[NumPar][i] );
865 CacheShift(B0p_MZ2_mHl2_0_mHp2_cache, NumPar, params, newResult);
872 double params[] = {MZ2, mHl2, mA2};
874 int i =
CacheCheck(B0p_MZ2_mHl2_0_mA2_cache, NumPar, params);
876 return ( B0p_MZ2_mHl2_0_mA2_cache[NumPar][i] );
879 CacheShift(B0p_MZ2_mHl2_0_mA2_cache, NumPar, params, newResult);
886 double params[] = {MZ2, mHl2};
888 int i =
CacheCheck(B0p_MZ2_mHl2_mHl2_mHl2_cache, NumPar, params);
890 return ( B0p_MZ2_mHl2_mHl2_mHl2_cache[NumPar][i] );
893 CacheShift(B0p_MZ2_mHl2_mHl2_mHl2_cache, NumPar, params, newResult);
900 double params[] = {MZ2, mHl2, mHh2};
902 int i =
CacheCheck(B0p_MZ2_mHl2_mHh2_mHl2_cache, NumPar, params);
904 return ( B0p_MZ2_mHl2_mHh2_mHl2_cache[NumPar][i] );
907 CacheShift(B0p_MZ2_mHl2_mHh2_mHl2_cache, NumPar, params, newResult);
914 double params[] = {MZ2, mHl2, mHh2};
916 int i =
CacheCheck(B0p_MZ2_mHl2_mHh2_mHh2_cache, NumPar, params);
918 return ( B0p_MZ2_mHl2_mHh2_mHh2_cache[NumPar][i] );
921 CacheShift(B0p_MZ2_mHl2_mHh2_mHh2_cache, NumPar, params, newResult);
928 double params[] = {MZ2, mHl2, mHp2};
930 int i =
CacheCheck(B0p_MZ2_mHl2_mHp2_mHp2_cache, NumPar, params);
932 return ( B0p_MZ2_mHl2_mHp2_mHp2_cache[NumPar][i] );
935 CacheShift(B0p_MZ2_mHl2_mHp2_mHp2_cache, NumPar, params, newResult);
942 double params[] = {MZ2, mHl2, mA2};
944 int i =
CacheCheck(B0p_MZ2_mHl2_mA2_mA2_cache, NumPar, params);
946 return ( B0p_MZ2_mHl2_mA2_mA2_cache[NumPar][i] );
949 CacheShift(B0p_MZ2_mHl2_mA2_mA2_cache, NumPar, params, newResult);
956 double params[] = {MZ2, mHh2};
958 int i =
CacheCheck(B0p_MZ2_0_0_mHh2_cache, NumPar, params);
960 return ( B0p_MZ2_0_0_mHh2_cache[NumPar][i] );
963 CacheShift(B0p_MZ2_0_0_mHh2_cache, NumPar, params, newResult);
970 double params[] = {MZ2, mHh2, mHp2};
972 int i =
CacheCheck(B0p_MZ2_mHh2_0_mHp2_cache, NumPar, params);
974 return ( B0p_MZ2_mHh2_0_mHp2_cache[NumPar][i] );
977 CacheShift(B0p_MZ2_mHh2_0_mHp2_cache, NumPar, params, newResult);
984 double params[] = {MZ2, mHh2, mA2};
986 int i =
CacheCheck(B0p_MZ2_mHh2_0_mA2_cache, NumPar, params);
988 return ( B0p_MZ2_mHh2_0_mA2_cache[NumPar][i] );
991 CacheShift(B0p_MZ2_mHh2_0_mA2_cache, NumPar, params, newResult);
998 double params[] = {MZ2, mHh2, mHl2};
1000 int i =
CacheCheck(B0p_MZ2_mHh2_mHl2_mHl2_cache, NumPar, params);
1002 return ( B0p_MZ2_mHh2_mHl2_mHl2_cache[NumPar][i] );
1005 CacheShift(B0p_MZ2_mHh2_mHl2_mHl2_cache, NumPar, params, newResult);
1012 double params[] = {MZ2, mHh2, mHl2};
1014 int i =
CacheCheck(B0p_MZ2_mHh2_mHh2_mHl2_cache, NumPar, params);
1016 return ( B0p_MZ2_mHh2_mHh2_mHl2_cache[NumPar][i] );
1019 CacheShift(B0p_MZ2_mHh2_mHh2_mHl2_cache, NumPar, params, newResult);
1026 double params[] = {MZ2, mHh2};
1028 int i =
CacheCheck(B0p_MZ2_mHh2_mHh2_mHh2_cache, NumPar, params);
1030 return ( B0p_MZ2_mHh2_mHh2_mHh2_cache[NumPar][i] );
1033 CacheShift(B0p_MZ2_mHh2_mHh2_mHh2_cache, NumPar, params, newResult);
1040 double params[] = {MZ2, mHh2, mHp2};
1042 int i =
CacheCheck(B0p_MZ2_mHh2_mHp2_mHp2_cache, NumPar, params);
1044 return ( B0p_MZ2_mHh2_mHp2_mHp2_cache[NumPar][i] );
1047 CacheShift(B0p_MZ2_mHh2_mHp2_mHp2_cache, NumPar, params, newResult);
1054 double params[] = {MZ2, mHh2, mA2};
1056 int i =
CacheCheck(B0p_MZ2_mHh2_mA2_mA2_cache, NumPar, params);
1058 return ( B0p_MZ2_mHh2_mA2_mA2_cache[NumPar][i] );
1061 CacheShift(B0p_MZ2_mHh2_mA2_mA2_cache, NumPar, params, newResult);
1068 double params[] = {MZ2, mHp2, mHl2};
1070 int i =
CacheCheck(B0p_MZ2_mHp2_0_mHl2_cache, NumPar, params);
1072 return ( B0p_MZ2_mHp2_0_mHl2_cache[NumPar][i] );
1075 CacheShift(B0p_MZ2_mHp2_0_mHl2_cache, NumPar, params, newResult);
1082 double params[] = {MZ2, mHp2, mHh2};
1084 int i =
CacheCheck(B0p_MZ2_mHp2_0_mHh2_cache, NumPar, params);
1086 return ( B0p_MZ2_mHp2_0_mHh2_cache[NumPar][i] );
1089 CacheShift(B0p_MZ2_mHp2_0_mHh2_cache, NumPar, params, newResult);
1096 double params[] = {MZ2, mHp2, mA2};
1098 int i =
CacheCheck(B0p_MZ2_mHp2_0_mA2_cache, NumPar, params);
1100 return ( B0p_MZ2_mHp2_0_mA2_cache[NumPar][i] );
1103 CacheShift(B0p_MZ2_mHp2_0_mA2_cache, NumPar, params, newResult);
1110 double params[] = {MZ2, mHp2, mHl2};
1112 int i =
CacheCheck(B0p_MZ2_mHp2_mHp2_mHl2_cache, NumPar, params);
1114 return ( B0p_MZ2_mHp2_mHp2_mHl2_cache[NumPar][i] );
1117 CacheShift(B0p_MZ2_mHp2_mHp2_mHl2_cache, NumPar, params, newResult);
1124 double params[] = {MZ2, mHp2, mHh2};
1126 int i =
CacheCheck(B0p_MZ2_mHp2_mHp2_mHh2_cache, NumPar, params);
1128 return ( B0p_MZ2_mHp2_mHp2_mHh2_cache[NumPar][i] );
1131 CacheShift(B0p_MZ2_mHp2_mHp2_mHh2_cache, NumPar, params, newResult);
1138 double params[] = {MZ2, mA2, mHl2};
1140 int i =
CacheCheck(B0p_MZ2_mA2_0_mHl2_cache, NumPar, params);
1142 return ( B0p_MZ2_mA2_0_mHl2_cache[NumPar][i] );
1145 CacheShift(B0p_MZ2_mA2_0_mHl2_cache, NumPar, params, newResult);
1152 double params[] = {MZ2, mA2, mHh2};
1154 int i =
CacheCheck(B0p_MZ2_mA2_0_mHh2_cache, NumPar, params);
1156 return ( B0p_MZ2_mA2_0_mHh2_cache[NumPar][i] );
1159 CacheShift(B0p_MZ2_mA2_0_mHh2_cache, NumPar, params, newResult);
1166 double params[] = {MZ2, mA2, mHp2};
1168 int i =
CacheCheck(B0p_MZ2_mA2_0_mHp2_cache, NumPar, params);
1170 return ( B0p_MZ2_mA2_0_mHp2_cache[NumPar][i] );
1173 CacheShift(B0p_MZ2_mA2_0_mHp2_cache, NumPar, params, newResult);
1180 double params[] = {MZ2, mA2, mHl2};
1182 int i =
CacheCheck(B0p_MZ2_mA2_mA2_mHl2_cache, NumPar, params);
1184 return ( B0p_MZ2_mA2_mA2_mHl2_cache[NumPar][i] );
1187 CacheShift(B0p_MZ2_mA2_mA2_mHl2_cache, NumPar, params, newResult);
1194 double params[] = {MZ2, mA2, mHh2};
1196 int i =
CacheCheck(B0p_MZ2_mA2_mA2_mHh2_cache, NumPar, params);
1198 return ( B0p_MZ2_mA2_mA2_mHh2_cache[NumPar][i] );
1201 CacheShift(B0p_MZ2_mA2_mA2_mHh2_cache, NumPar, params, newResult);
1210 double params[] = {MZ2, mA2, mHp2};
1212 int i =
CacheCheck(B00_MZ2_0_mA2_mHp2_cache, NumPar, params);
1214 return ( B00_MZ2_0_mA2_mHp2_cache[NumPar][i] );
1217 CacheShift(B00_MZ2_0_mA2_mHp2_cache, NumPar, params, newResult);
1224 double params[] = {MZ2, mHh2, mA2};
1226 int i =
CacheCheck(B00_MZ2_0_mHh2_mA2_cache, NumPar, params);
1228 return ( B00_MZ2_0_mHh2_mA2_cache[NumPar][i] );
1231 CacheShift(B00_MZ2_0_mHh2_mA2_cache, NumPar, params, newResult);
1238 double params[] = {MZ2, mHh2, mHp2};
1240 int i =
CacheCheck(B00_MZ2_0_mHh2_mHp2_cache, NumPar, params);
1242 return ( B00_MZ2_0_mHh2_mHp2_cache[NumPar][i] );
1245 CacheShift(B00_MZ2_0_mHh2_mHp2_cache, NumPar, params, newResult);
1252 double params[] = {MZ2, mHl2, mA2};
1254 int i =
CacheCheck(B00_MZ2_0_mHl2_mA2_cache, NumPar, params);
1256 return ( B00_MZ2_0_mHl2_mA2_cache[NumPar][i] );
1259 CacheShift(B00_MZ2_0_mHl2_mA2_cache, NumPar, params, newResult);
1266 double params[] = {MZ2, mHl2, mHp2};
1268 int i =
CacheCheck(B00_MZ2_0_mHl2_mHp2_cache, NumPar, params);
1270 return ( B00_MZ2_0_mHl2_mHp2_cache[NumPar][i] );
1273 CacheShift(B00_MZ2_0_mHl2_mHp2_cache, NumPar, params, newResult);
1280 double params[] = {MZ2, mHp2};
1282 int i =
CacheCheck(B00_MZ2_0_mHp2_mHp2_cache, NumPar, params);
1284 return ( B00_MZ2_0_mHp2_mHp2_cache[NumPar][i] );
1287 CacheShift(B00_MZ2_0_mHp2_mHp2_cache, NumPar, params, newResult);
1294 double params[] = {MZ2, MW2, mHh2};
1296 int i =
CacheCheck(B00_MZ2_0_MW2_mHh2_cache, NumPar, params);
1298 return ( B00_MZ2_0_MW2_mHh2_cache[NumPar][i] );
1301 CacheShift(B00_MZ2_0_MW2_mHh2_cache, NumPar, params, newResult);
1308 double params[] = {MZ2, MW2, mHl2};
1310 int i =
CacheCheck(B00_MZ2_0_MW2_mHl2_cache, NumPar, params);
1312 return ( B00_MZ2_0_MW2_mHl2_cache[NumPar][i] );
1315 CacheShift(B00_MZ2_0_MW2_mHl2_cache, NumPar, params, newResult);
1322 double params[] = {MZ2, mHh2};
1324 int i =
CacheCheck(B00_MZ2_0_MZ2_mHh2_cache, NumPar, params);
1326 return ( B00_MZ2_0_MZ2_mHh2_cache[NumPar][i] );
1329 CacheShift(B00_MZ2_0_MZ2_mHh2_cache, NumPar, params, newResult);
1336 double params[] = {MZ2, mHl2};
1338 int i =
CacheCheck(B00_MZ2_0_MZ2_mHl2_cache, NumPar, params);
1340 return ( B00_MZ2_0_MZ2_mHl2_cache[NumPar][i] );
1343 CacheShift(B00_MZ2_0_MZ2_mHl2_cache, NumPar, params, newResult);
1350 double params[] = {MZ2, MW2, mA2, mHp2};
1352 int i =
CacheCheck(B00_MZ2_MW2_mA2_mHp2_cache, NumPar, params);
1354 return ( B00_MZ2_MW2_mA2_mHp2_cache[NumPar][i] );
1357 CacheShift(B00_MZ2_MW2_mA2_mHp2_cache, NumPar, params, newResult);
1364 double params[] = {MZ2, MW2, mHh2, mHp2};
1366 int i =
CacheCheck(B00_MZ2_MW2_mHh2_mHp2_cache, NumPar, params);
1368 return ( B00_MZ2_MW2_mHh2_mHp2_cache[NumPar][i] );
1371 CacheShift(B00_MZ2_MW2_mHh2_mHp2_cache, NumPar, params, newResult);
1378 double params[] = {MZ2, MW2, mHl2, mHp2};
1380 int i =
CacheCheck(B00_MZ2_MW2_mHl2_mHp2_cache, NumPar, params);
1382 return ( B00_MZ2_MW2_mHl2_mHp2_cache[NumPar][i] );
1385 CacheShift(B00_MZ2_MW2_mHl2_mHp2_cache, NumPar, params, newResult);
1392 double params[] = {MZ2, MW2, mHp2};
1394 int i =
CacheCheck(B00_MZ2_MW2_mHp2_mHp2_cache, NumPar, params);
1396 return ( B00_MZ2_MW2_mHp2_mHp2_cache[NumPar][i] );
1399 CacheShift(B00_MZ2_MW2_mHp2_mHp2_cache, NumPar, params, newResult);
1406 double params[] = {MZ2, MW2, mHh2};
1408 int i =
CacheCheck(B00_MZ2_MW2_MW2_mHh2_cache, NumPar, params);
1410 return ( B00_MZ2_MW2_MW2_mHh2_cache[NumPar][i] );
1413 CacheShift(B00_MZ2_MW2_MW2_mHh2_cache, NumPar, params, newResult);
1420 double params[] = {MZ2, MW2, mHl2};
1422 int i =
CacheCheck(B00_MZ2_MW2_MW2_mHl2_cache, NumPar, params);
1424 return ( B00_MZ2_MW2_MW2_mHl2_cache[NumPar][i] );
1427 CacheShift(B00_MZ2_MW2_MW2_mHl2_cache, NumPar, params, newResult);
1434 double params[] = {MZ2, mHh2, mA2};
1436 int i =
CacheCheck(B00_MZ2_MZ2_mHh2_mA2_cache, NumPar, params);
1438 return ( B00_MZ2_MZ2_mHh2_mA2_cache[NumPar][i] );
1441 CacheShift(B00_MZ2_MZ2_mHh2_mA2_cache, NumPar, params, newResult);
1448 double params[] = {MZ2, mHl2, mA2};
1450 int i =
CacheCheck(B00_MZ2_MZ2_mHl2_mA2_cache, NumPar, params);
1452 return ( B00_MZ2_MZ2_mHl2_mA2_cache[NumPar][i] );
1455 CacheShift(B00_MZ2_MZ2_mHl2_mA2_cache, NumPar, params, newResult);
1462 double params[] = {MZ2, mHp2};
1464 int i =
CacheCheck(B00_MZ2_MZ2_mHp2_mHp2_cache, NumPar, params);
1466 return ( B00_MZ2_MZ2_mHp2_mHp2_cache[NumPar][i] );
1469 CacheShift(B00_MZ2_MZ2_mHp2_mHp2_cache, NumPar, params, newResult);
1476 double params[] = {MZ2, mHh2};
1478 int i =
CacheCheck(B00_MZ2_MZ2_MZ2_mHh2_cache, NumPar, params);
1480 return ( B00_MZ2_MZ2_MZ2_mHh2_cache[NumPar][i] );
1483 CacheShift(B00_MZ2_MZ2_MZ2_mHh2_cache, NumPar, params, newResult);
1490 double params[] = {MZ2, mHl2};
1492 int i =
CacheCheck(B00_MZ2_MZ2_MZ2_mHl2_cache, NumPar, params);
1494 return ( B00_MZ2_MZ2_MZ2_mHl2_cache[NumPar][i] );
1497 CacheShift(B00_MZ2_MZ2_MZ2_mHl2_cache, NumPar, params, newResult);
1507 std::ifstream INfile;
1508 std::string lineTab;
1509 INfile.open( filename.c_str() );
1511 std::cout<<
"error: in GeneralTHDMcache, table doesn't exist!"<<std::endl;
1519 while(INfile.good()){
1520 while(getline(INfile, lineTab)){
1521 if( lineTab[0]==
'#' )
continue;
1523 std::istringstream streamTab(lineTab);
1525 while(streamTab >>v){
1543 int rowN=arrayTab.
size_i();
1545 double xmin = arrayTab(0,0);
1546 double xmax = arrayTab(rowN-1,0);
1547 double interval = arrayTab(1,0)-arrayTab(0,0);
1548 int Nintervals = (x-xmin)/interval;
1560 y =(arrayTab(Nintervals+1,1)-arrayTab(Nintervals,1))/(arrayTab(Nintervals+1,0)
1561 -arrayTab(Nintervals,0))*(x-arrayTab(Nintervals,0))+arrayTab(Nintervals,1);
1570 int rowN=arrayTab.
size_i();
1572 double xmin = arrayTab(0,0);
1573 double xmax = arrayTab(rowN-1,0);
1574 double ymin = arrayTab(0,1);
1575 double ymax = arrayTab(rowN-1,1);
1576 double intervalx = arrayTab(1,0)-arrayTab(0,0);
1579 while(arrayTab(i,1)-arrayTab(i-1,1)==0&&i<30000);
1580 double intervaly = arrayTab(i,1)-arrayTab(i-1,1);
1581 int Nintervalsx = (x-xmin)/intervalx;
1582 int Nintervalsy = (y-ymin)/intervaly;
1583 if(x<xmin||x>xmax||y<ymin||y>ymax){
1588 double x1=arrayTab(i*Nintervalsy+Nintervalsx,0);
1589 double x2=arrayTab(i*Nintervalsy+Nintervalsx+1,0);
1590 double y1=arrayTab(i*Nintervalsy+Nintervalsx,1);
1591 double y2=arrayTab(i*(Nintervalsy+1)+Nintervalsx,1);
1592 return (arrayTab(i*Nintervalsy+Nintervalsx,2) * (x2-x) * (y2-y)
1593 +arrayTab(i*Nintervalsy+Nintervalsx+1,2) * (x-x1) * (y2-y)
1594 +arrayTab(i*(Nintervalsy+1)+Nintervalsx,2) * (x2-x) * (y-y1)
1595 +arrayTab(i*(Nintervalsy+1)+Nintervalsx+1,2) * (x-x1) * (y-y1))
1601 std::stringstream br1,br2,br3,br4,br5,br6,br7;
1602 std::stringstream dw1;
1603 std::stringstream cs1,cs2,cs3,cs4,cs5,cs6,cs7,cs8,cs9;
1604 std::stringstream cs11,cs12,cs13,cs14,cs15,cs16,cs17,cs18,cs19;
1605 std::stringstream cs20,cs21;
1606 std::stringstream csr1,csr2,csr3,csr4;
1607 std::stringstream csr11,csr12,csr13,csr14;
1608 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;
1609 std::stringstream ex24,ex25,ex26,ex27,ex28,ex29,ex30,ex31,ex32,ex33,ex34,ex35,ex36,ex37,ex38,ex39,ex40,ex41,ex42,ex43,ex44,\
1610 ex45,ex46,ex46a,ex47,ex48,ex49,ex50,ex51,ex52,ex53,ex54,ex55,ex56;
1611 std::stringstream ex57,ex58,ex59,ex60,ex61,ex62,ex63,ex64,ex65,ex66,ex67,ex68,ex69,ex70,ex71,ex72,ex73,ex74,ex75,ex76,ex77,\
1613 std::stringstream bsg1;
1615 std::cout<<
"reading tables"<<std::endl;
1618 std::stringstream path;
1619 path << getenv(
"HEPFITTABS") <<
"/THDM/tabs/";
1620 std::string tablepath=path.str();
1622 br1 << tablepath <<
"br1.dat";
1624 br2 << tablepath <<
"br2.dat";
1626 br3 << tablepath <<
"br3.dat";
1628 br4 << tablepath <<
"br4.dat";
1630 br5 << tablepath <<
"br5.dat";
1632 br6 << tablepath <<
"br6.dat";
1634 br7 << tablepath <<
"br7.dat";
1636 dw1 << tablepath <<
"dw1.dat";
1638 cs1 << tablepath <<
"log_cs_ggH_8.dat";
1640 cs11 << tablepath <<
"log_cs_ggH_13.dat";
1642 cs2 << tablepath <<
"log_cs_VBF_8.dat";
1644 cs12 << tablepath <<
"log_cs_VBF_13.dat";
1646 cs3 << tablepath <<
"log_cs_WH_8.dat";
1648 cs13 << tablepath <<
"log_cs_WH_13.dat";
1650 cs4 << tablepath <<
"log_cs_ZH_8.dat";
1652 cs14 << tablepath <<
"log_cs_ZH_13.dat";
1654 cs5 << tablepath <<
"log_cs_ttH_8.dat";
1656 cs15 << tablepath <<
"log_cs_ttH_13.dat";
1658 cs6 << tablepath <<
"log_cs_bbH_8.dat";
1660 cs16 << tablepath <<
"log_cs_bbH_13.dat";
1662 cs7 << tablepath <<
"log_cs_ggA_8.dat";
1664 cs17 << tablepath <<
"log_cs_ggA_13.dat";
1666 cs8 << tablepath <<
"log_cs_ttA_8.dat";
1668 cs18 << tablepath <<
"log_cs_ttA_13.dat";
1670 cs9 << tablepath <<
"log_cs_bbA_8.dat";
1672 cs19 << tablepath <<
"log_cs_bbA_13.dat";
1674 cs20 << tablepath <<
"log_cs_ggHp_8.dat";
1676 cs21 << tablepath <<
"log_cs_ggHp_13.dat";
1678 csr1 << tablepath <<
"csrH_top_8.dat";
1680 csr11 << tablepath <<
"csrH_top_13.dat";
1682 csr2 << tablepath <<
"csrH_bottom_8.dat";
1684 csr12 << tablepath <<
"csrH_bottom_13.dat";
1686 csr3 << tablepath <<
"csrA_top_8.dat";
1688 csr13 << tablepath <<
"csrA_top_13.dat";
1690 csr4 << tablepath <<
"csrA_bottom_8.dat";
1692 csr14 << tablepath <<
"csrA_bottom_13.dat";
1695 ex1 << tablepath <<
"14096064_a.dat";
1697 ex2 << tablepath <<
"14096064_b.dat";
1699 ex3 << tablepath <<
"14076583.dat";
1701 ex4 << tablepath <<
"14078150.dat";
1703 ex5 << tablepath <<
"150705930_a.dat";
1705 ex6 << tablepath <<
"150705930_b.dat";
1707 ex7 << tablepath <<
"150900389_a.dat";
1709 ex8 << tablepath <<
"150900389_b.dat";
1711 ex9 << tablepath <<
"150904670.dat";
1713 ex10 << tablepath <<
"150204478_b.dat";
1715 ex11 << tablepath <<
"150204478_a.dat";
1717 ex12 << tablepath <<
"150400936.dat";
1719 ex13 << tablepath <<
"150608329.dat";
1721 ex14 << tablepath <<
"180206149.dat";
1723 ex15 << tablepath <<
"CMS-PAS-HIG-14-029_a.dat";
1725 ex16 << tablepath <<
"CMS-PAS-HIG-14-029_b.dat";
1727 ex17 << tablepath <<
"CMS-PAS-HIG-16-014.dat";
1729 ex18 << tablepath <<
"150304114.dat";
1731 ex19 << tablepath <<
"160306896.dat";
1733 ex20 << tablepath <<
"151001181_a.dat";
1735 ex21 << tablepath <<
"170700350.dat";
1737 ex22 << tablepath <<
"150404710.dat";
1739 ex23 << tablepath <<
"151001181_b.dat";
1742 ex24 << tablepath <<
"160302991_a.dat";
1744 ex25 << tablepath <<
"160302991_b.dat";
1746 ex26 << tablepath <<
"160302991_c.dat";
1748 ex27 << tablepath <<
"160302991_d.dat";
1751 ex28 << tablepath <<
"180711883.dat";
1753 ex29 << tablepath <<
"ATLAS-CONF-2016-104_b.dat";
1755 ex30 << tablepath <<
"170907242_a.dat";
1757 ex31 << tablepath <<
"170907242_b.dat";
1759 ex32 << tablepath <<
"170704147.dat";
1761 ex33 << tablepath <<
"170800212.dat";
1763 ex34 << tablepath <<
"180501908.dat";
1765 ex35 << tablepath <<
"171206386_a.dat";
1767 ex36 << tablepath <<
"171206386_b.dat";
1769 ex37 << tablepath <<
"170809638_a.dat";
1771 ex38 << tablepath <<
"170809638_b.dat";
1773 ex39 << tablepath <<
"171001123_a.dat";
1775 ex40 << tablepath <<
"171001123_b.dat";
1777 ex41 << tablepath <<
"171007235_a.dat";
1779 ex42 << tablepath <<
"171007235_b.dat";
1781 ex43 << tablepath <<
"170804445.dat";
1783 ex44 << tablepath <<
"180406174.dat";
1785 ex45 << tablepath <<
"180704873.dat";
1787 ex46 << tablepath <<
"180800336.dat";
1789 ex46a << tablepath <<
"181104671.dat";
1791 ex47 << tablepath <<
"180708567.dat";
1793 ex48 << tablepath <<
"171206518_a.dat";
1795 ex49 << tablepath <<
"171206518_b.dat";
1798 ex50 << tablepath <<
"180401126_a.dat";
1800 ex51 << tablepath <<
"180401126_b.dat";
1803 ex52 << tablepath <<
"CMS-PAS-HIG-16-025.dat";
1805 ex53 << tablepath <<
"180512191.dat";
1807 ex54 << tablepath <<
"180306553_a.dat";
1809 ex55 << tablepath <<
"180306553_b.dat";
1811 ex56 << tablepath <<
"160902507.dat";
1813 ex57 << tablepath <<
"171203143.dat";
1815 ex58 << tablepath <<
"180401939_a.dat";
1817 ex59 << tablepath <<
"180303838.dat";
1819 ex60 << tablepath <<
"CMS-PAS-HIG-16-023.dat";
1821 ex61 << tablepath <<
"180209407.dat";
1823 ex62 << tablepath <<
"180603548.dat";
1825 ex63 << tablepath <<
"180801473.dat";
1827 ex64 << tablepath <<
"180600408.dat";
1829 ex65 << tablepath <<
"170702909.dat";
1831 ex66 << tablepath <<
"180801365.dat";
1833 ex67 << tablepath <<
"170804188.dat";
1835 ex68 << tablepath <<
"CMS-PAS-HIG-18-005_a.dat";
1837 ex69 << tablepath <<
"180702826_a.dat";
1839 ex70 << tablepath <<
"CMS-PAS-HIG-18-005_b.dat";
1841 ex71 << tablepath <<
"180702826_b.dat";
1844 ex72 << tablepath <<
"14126663.dat";
1846 ex73 << tablepath <<
"151203704.dat";
1848 ex74 << tablepath <<
"150807774_a.dat";
1850 ex75 << tablepath <<
"150807774_b.dat";
1852 ex76 << tablepath <<
"180707915.dat";
1854 ex77 << tablepath <<
"180803599.dat";
1856 ex78 << tablepath <<
"CMS-PAS-HIG-16-031.dat";
1859 bsg1 << tablepath <<
"bsgammatable.dat";
1867 double params[] = {mass};
1871 return ( ip_Br_HPtott_cache[NumPar][i] );
1883 double params[] = {mass};
1887 return ( ip_Br_HPtobb_cache[NumPar][i] );
1899 double params[] = {mass};
1903 return ( ip_Br_HPtotautau_cache[NumPar][i] );
1906 CacheShiftReal(ip_Br_HPtotautau_cache, NumPar, params, newResult);
1915 double params[] = {mass};
1919 return ( ip_Br_HPtocc_cache[NumPar][i] );
1931 double params[] = {mass};
1935 return ( ip_Br_HPtomumu_cache[NumPar][i] );
1947 double params[] = {mass};
1951 return ( ip_Br_HPtoZZ_cache[NumPar][i] );
1963 double params[] = {mass};
1967 return ( ip_Br_HPtoWW_cache[NumPar][i] );
1979 double params[] = {mass};
1983 return ( ip_GammaHPtotSM_cache[NumPar][i] );
1986 CacheShiftReal(ip_GammaHPtotSM_cache, NumPar, params, newResult);
1995 double params[] = {mass};
1999 return ( ip_cs_ggtoH_8_cache[NumPar][i] );
2001 double newResult = 0.0;
2002 if (mass>=20. && mass <=2000.) {
2014 double params[] = {mass};
2018 return ( ip_cs_ggtoH_13_cache[NumPar][i] );
2020 double newResult = 0.0;
2021 if (mass>=20. && mass <=2000.) {
2033 double params[] = {mass};
2037 return ( ip_cs_VBFtoH_8_cache[NumPar][i] );
2039 double newResult = 0.0;
2040 if (mass>=20. && mass <=2000.) {
2052 double params[] = {mass};
2056 return ( ip_cs_VBFtoH_13_cache[NumPar][i] );
2058 double newResult = 0.0;
2059 if (mass>=20. && mass <=2000.) {
2062 CacheShiftReal(ip_cs_VBFtoH_13_cache, NumPar, params, newResult);
2071 double params[] = {mass};
2075 return ( ip_cs_WtoWH_8_cache[NumPar][i] );
2077 double newResult = 0.0;
2078 if (mass>=20. && mass <=2000.) {
2090 double params[] = {mass};
2094 return ( ip_cs_WtoWH_13_cache[NumPar][i] );
2096 double newResult = 0.0;
2097 if (mass>=20. && mass <=2000.) {
2109 double params[] = {mass};
2113 return ( ip_cs_ZtoZH_8_cache[NumPar][i] );
2115 double newResult = 0.0;
2116 if (mass>=20. && mass <=2000.) {
2128 double params[] = {mass};
2132 return ( ip_cs_ZtoZH_13_cache[NumPar][i] );
2134 double newResult = 0.0;
2135 if (mass>=20. && mass <=2000.) {
2147 double params[] = {mass};
2151 return ( ip_cs_pptottH_8_cache[NumPar][i] );
2153 double newResult = 0.0;
2154 if (mass>=20. && mass <=2000.) {
2157 CacheShiftReal(ip_cs_pptottH_8_cache, NumPar, params, newResult);
2166 double params[] = {mass};
2170 return ( ip_cs_pptottH_13_cache[NumPar][i] );
2172 double newResult = 0.0;
2173 if (mass>=20. && mass <=2000.) {
2176 CacheShiftReal(ip_cs_pptottH_13_cache, NumPar, params, newResult);
2185 double params[] = {mass};
2189 return ( ip_cs_pptobbH_8_cache[NumPar][i] );
2191 double newResult = 0.0;
2192 if (mass>=20. && mass <=2000.) {
2195 CacheShiftReal(ip_cs_pptobbH_8_cache, NumPar, params, newResult);
2204 double params[] = {mass};
2208 return ( ip_cs_pptobbH_13_cache[NumPar][i] );
2210 double newResult = 0.0;
2211 if (mass>=20. && mass <=2000.) {
2214 CacheShiftReal(ip_cs_pptobbH_13_cache, NumPar, params, newResult);
2223 double params[] = {mass};
2227 return ( ip_cs_ggtoA_8_cache[NumPar][i] );
2229 double newResult = 0.0;
2230 if (mass>=20. && mass <=2000.) {
2242 double params[] = {mass};
2246 return ( ip_cs_ggtoA_13_cache[NumPar][i] );
2248 double newResult = 0.0;
2249 if (mass>=20. && mass <=2000.) {
2261 double params[] = {mass};
2265 return ( ip_cs_pptottA_8_cache[NumPar][i] );
2267 double newResult = 0.0;
2268 if (mass>=20. && mass <=2000.) {
2271 CacheShiftReal(ip_cs_pptottA_8_cache, NumPar, params, newResult);
2280 double params[] = {mass};
2284 return ( ip_cs_pptottA_13_cache[NumPar][i] );
2286 double newResult = 0.0;
2287 if (mass>=20. && mass <=2000.) {
2290 CacheShiftReal(ip_cs_pptottA_13_cache, NumPar, params, newResult);
2299 double params[] = {mass};
2303 return ( ip_cs_pptobbA_8_cache[NumPar][i] );
2305 double newResult = 0.0;
2306 if (mass>=20. && mass <=2000.) {
2309 CacheShiftReal(ip_cs_pptobbA_8_cache, NumPar, params, newResult);
2318 double params[] = {mass};
2322 return ( ip_cs_pptobbA_13_cache[NumPar][i] );
2324 double newResult = 0.0;
2325 if (mass>=20. && mass <=2000.) {
2328 CacheShiftReal(ip_cs_pptobbA_13_cache, NumPar, params, newResult);
2337 double params[] = {mHp, logtb};
2341 return ( ip_cs_ggtoHp_8_cache[NumPar][i] );
2343 double newResult = 0.0;
2344 if (mHp>=180. && mHp <=1400. && logtb>=-1. && logtb<=1.75) {
2356 double params[] = {mHp, logtb};
2360 return ( ip_cs_ggtoHp_13_cache[NumPar][i] );
2362 double newResult = 0.0;
2363 if (mHp>=180. && mHp <=2000. && logtb>=-1. && logtb<=1.75) {
2366 CacheShiftReal(ip_cs_ggtoHp_13_cache, NumPar, params, newResult);
2375 double params[] = {mass};
2379 return ( ip_csr_ggH_t_8_cache[NumPar][i] );
2391 double params[] = {mass};
2395 return ( ip_csr_ggH_t_13_cache[NumPar][i] );
2398 CacheShiftReal(ip_csr_ggH_t_13_cache, NumPar, params, newResult);
2407 double params[] = {mass};
2411 return ( ip_csr_ggH_b_8_cache[NumPar][i] );
2423 double params[] = {mass};
2427 return ( ip_csr_ggH_b_13_cache[NumPar][i] );
2430 CacheShiftReal(ip_csr_ggH_b_13_cache, NumPar, params, newResult);
2439 double params[] = {mass};
2443 return ( ip_csr_ggA_t_8_cache[NumPar][i] );
2455 double params[] = {mass};
2459 return ( ip_csr_ggA_t_13_cache[NumPar][i] );
2462 CacheShiftReal(ip_csr_ggA_t_13_cache, NumPar, params, newResult);
2471 double params[] = {mass};
2475 return ( ip_csr_ggA_b_8_cache[NumPar][i] );
2487 double params[] = {mass};
2491 return ( ip_csr_ggA_b_13_cache[NumPar][i] );
2494 CacheShiftReal(ip_csr_ggA_b_13_cache, NumPar, params, newResult);
2501 double params[] = {mass};
2503 int i =
CacheCheckReal(ip_ex_tt_phi_tt_ATLAS13_cache, NumPar, params);
2505 return ( ip_ex_tt_phi_tt_ATLAS13_cache[NumPar][i] );
2508 CacheShiftReal(ip_ex_tt_phi_tt_ATLAS13_cache, NumPar, params, newResult);
2515 double params[] = {mass};
2517 int i =
CacheCheckReal(ip_ex_bb_phi_tt_ATLAS13_cache, NumPar, params);
2519 return ( ip_ex_bb_phi_tt_ATLAS13_cache[NumPar][i] );
2522 CacheShiftReal(ip_ex_bb_phi_tt_ATLAS13_cache, NumPar, params, newResult);
2529 double params[] = {mass};
2531 int i =
CacheCheckReal(ip_ex_bb_phi_bb_CMS8_cache, NumPar, params);
2533 return ( ip_ex_bb_phi_bb_CMS8_cache[NumPar][i] );
2536 CacheShiftReal(ip_ex_bb_phi_bb_CMS8_cache, NumPar, params, newResult);
2543 double params[] = {mass};
2545 int i =
CacheCheckReal(ip_ex_gg_phi_bb_CMS8_cache, NumPar, params);
2547 return ( ip_ex_gg_phi_bb_CMS8_cache[NumPar][i] );
2550 CacheShiftReal(ip_ex_gg_phi_bb_CMS8_cache, NumPar, params, newResult);
2557 double params[] = {mass};
2559 int i =
CacheCheckReal(ip_ex_pp_phi_bb_CMS13_cache, NumPar, params);
2561 return ( ip_ex_pp_phi_bb_CMS13_cache[NumPar][i] );
2564 CacheShiftReal(ip_ex_pp_phi_bb_CMS13_cache, NumPar, params, newResult);
2571 double params[] = {mass};
2573 int i =
CacheCheckReal(ip_ex_bb_phi_bb_CMS13_cache, NumPar, params);
2575 return ( ip_ex_bb_phi_bb_CMS13_cache[NumPar][i] );
2578 CacheShiftReal(ip_ex_bb_phi_bb_CMS13_cache, NumPar, params, newResult);
2585 double params[] = {mass};
2587 int i =
CacheCheckReal(ip_ex_gg_phi_tautau_ATLAS8_cache, NumPar, params);
2589 return ( ip_ex_gg_phi_tautau_ATLAS8_cache[NumPar][i] );
2592 CacheShiftReal(ip_ex_gg_phi_tautau_ATLAS8_cache, NumPar, params, newResult);
2599 double params[] = {mass};
2601 int i =
CacheCheckReal(ip_ex_gg_phi_tautau_CMS8_cache, NumPar, params);
2603 return ( ip_ex_gg_phi_tautau_CMS8_cache[NumPar][i] );
2606 CacheShiftReal(ip_ex_gg_phi_tautau_CMS8_cache, NumPar, params, newResult);
2613 double params[] = {mass};
2615 int i =
CacheCheckReal(ip_ex_bb_phi_tautau_ATLAS8_cache, NumPar, params);
2617 return ( ip_ex_bb_phi_tautau_ATLAS8_cache[NumPar][i] );
2620 CacheShiftReal(ip_ex_bb_phi_tautau_ATLAS8_cache, NumPar, params, newResult);
2627 double params[] = {mass};
2629 int i =
CacheCheckReal(ip_ex_bb_phi_tautau_CMS8_cache, NumPar, params);
2631 return ( ip_ex_bb_phi_tautau_CMS8_cache[NumPar][i] );
2634 CacheShiftReal(ip_ex_bb_phi_tautau_CMS8_cache, NumPar, params, newResult);
2641 double params[] = {mass};
2643 int i =
CacheCheckReal(ip_ex_gg_phi_tautau_ATLAS13_cache, NumPar, params);
2645 return ( ip_ex_gg_phi_tautau_ATLAS13_cache[NumPar][i] );
2648 CacheShiftReal(ip_ex_gg_phi_tautau_ATLAS13_cache, NumPar, params, newResult);
2655 double params[] = {mass};
2657 int i =
CacheCheckReal(ip_ex_gg_phi_tautau_CMS13_cache, NumPar, params);
2659 return ( ip_ex_gg_phi_tautau_CMS13_cache[NumPar][i] );
2662 CacheShiftReal(ip_ex_gg_phi_tautau_CMS13_cache, NumPar, params, newResult);
2669 double params[] = {mass};
2671 int i =
CacheCheckReal(ip_ex_bb_phi_tautau_ATLAS13_cache, NumPar, params);
2673 return ( ip_ex_bb_phi_tautau_ATLAS13_cache[NumPar][i] );
2676 CacheShiftReal(ip_ex_bb_phi_tautau_ATLAS13_cache, NumPar, params, newResult);
2683 double params[] = {mass};
2685 int i =
CacheCheckReal(ip_ex_bb_phi_tautau_CMS13_cache, NumPar, params);
2687 return ( ip_ex_bb_phi_tautau_CMS13_cache[NumPar][i] );
2690 CacheShiftReal(ip_ex_bb_phi_tautau_CMS13_cache, NumPar, params, newResult);
2697 double params[] = {mass};
2699 int i =
CacheCheckReal(ip_ex_gg_phi_gaga_ATLAS8_cache, NumPar, params);
2701 return ( ip_ex_gg_phi_gaga_ATLAS8_cache[NumPar][i] );
2704 CacheShiftReal(ip_ex_gg_phi_gaga_ATLAS8_cache, NumPar, params, newResult);
2711 double params[] = {mass};
2713 int i =
CacheCheckReal(ip_ex_pp_phi_gaga_ATLAS13_cache, NumPar, params);
2715 return ( ip_ex_pp_phi_gaga_ATLAS13_cache[NumPar][i] );
2718 CacheShiftReal(ip_ex_pp_phi_gaga_ATLAS13_cache, NumPar, params, newResult);
2725 double params[] = {mass};
2727 int i =
CacheCheckReal(ip_ex_gg_phi_gaga_CMS13_cache, NumPar, params);
2729 return ( ip_ex_gg_phi_gaga_CMS13_cache[NumPar][i] );
2732 CacheShiftReal(ip_ex_gg_phi_gaga_CMS13_cache, NumPar, params, newResult);
2739 double params[] = {mass};
2741 int i =
CacheCheckReal(ip_ex_pp_phi_Zga_llga_ATLAS8_cache, NumPar, params);
2743 return ( ip_ex_pp_phi_Zga_llga_ATLAS8_cache[NumPar][i] );
2746 CacheShiftReal(ip_ex_pp_phi_Zga_llga_ATLAS8_cache, NumPar, params, newResult);
2753 double params[] = {mass};
2755 int i =
CacheCheckReal(ip_ex_pp_phi_Zga_llga_CMS8_cache, NumPar, params);
2757 return ( ip_ex_pp_phi_Zga_llga_CMS8_cache[NumPar][i] );
2760 CacheShiftReal(ip_ex_pp_phi_Zga_llga_CMS8_cache, NumPar, params, newResult);
2767 double params[] = {mass};
2769 int i =
CacheCheckReal(ip_ex_gg_phi_Zga_llga_ATLAS13_cache, NumPar, params);
2771 return ( ip_ex_gg_phi_Zga_llga_ATLAS13_cache[NumPar][i] );
2774 CacheShiftReal(ip_ex_gg_phi_Zga_llga_ATLAS13_cache, NumPar, params, newResult);
2781 double params[] = {mass};
2783 int i =
CacheCheckReal(ip_ex_gg_phi_Zga_qqga_ATLAS13_cache, NumPar, params);
2785 return ( ip_ex_gg_phi_Zga_qqga_ATLAS13_cache[NumPar][i] );
2788 CacheShiftReal(ip_ex_gg_phi_Zga_qqga_ATLAS13_cache, NumPar, params, newResult);
2795 double params[] = {mass};
2797 int i =
CacheCheckReal(ip_ex_gg_phi_Zga_CMS13_cache, NumPar, params);
2799 return ( ip_ex_gg_phi_Zga_CMS13_cache[NumPar][i] );
2802 CacheShiftReal(ip_ex_gg_phi_Zga_CMS13_cache, NumPar, params, newResult);
2809 double params[] = {mass};
2811 int i =
CacheCheckReal(ip_ex_gg_phi_ZZ_ATLAS8_cache, NumPar, params);
2813 return ( ip_ex_gg_phi_ZZ_ATLAS8_cache[NumPar][i] );
2816 CacheShiftReal(ip_ex_gg_phi_ZZ_ATLAS8_cache, NumPar, params, newResult);
2823 double params[] = {mass};
2825 int i =
CacheCheckReal(ip_ex_VV_phi_ZZ_ATLAS8_cache, NumPar, params);
2827 return ( ip_ex_VV_phi_ZZ_ATLAS8_cache[NumPar][i] );
2830 CacheShiftReal(ip_ex_VV_phi_ZZ_ATLAS8_cache, NumPar, params, newResult);
2837 double params[] = {mass};
2839 int i =
CacheCheckReal(ip_ex_gg_phi_ZZ_llllnunu_ATLAS13_cache, NumPar, params);
2841 return ( ip_ex_gg_phi_ZZ_llllnunu_ATLAS13_cache[NumPar][i] );
2844 CacheShiftReal(ip_ex_gg_phi_ZZ_llllnunu_ATLAS13_cache, NumPar, params, newResult);
2851 double params[] = {mass};
2853 int i =
CacheCheckReal(ip_ex_VV_phi_ZZ_llllnunu_ATLAS13_cache, NumPar, params);
2855 return ( ip_ex_VV_phi_ZZ_llllnunu_ATLAS13_cache[NumPar][i] );
2858 CacheShiftReal(ip_ex_VV_phi_ZZ_llllnunu_ATLAS13_cache, NumPar, params, newResult);
2865 double params[] = {mass};
2867 int i =
CacheCheckReal(ip_ex_gg_phi_ZZ_qqllnunu_ATLAS13_cache, NumPar, params);
2869 return ( ip_ex_gg_phi_ZZ_qqllnunu_ATLAS13_cache[NumPar][i] );
2872 CacheShiftReal(ip_ex_gg_phi_ZZ_qqllnunu_ATLAS13_cache, NumPar, params, newResult);
2879 double params[] = {mass};
2881 int i =
CacheCheckReal(ip_ex_VV_phi_ZZ_qqllnunu_ATLAS13_cache, NumPar, params);
2883 return ( ip_ex_VV_phi_ZZ_qqllnunu_ATLAS13_cache[NumPar][i] );
2886 CacheShiftReal(ip_ex_VV_phi_ZZ_qqllnunu_ATLAS13_cache, NumPar, params, newResult);
2893 double params[] = {mass};
2895 int i =
CacheCheckReal(ip_ex_pp_phi_ZZ_llqqnunull_CMS13_cache, NumPar, params);
2897 return ( ip_ex_pp_phi_ZZ_llqqnunull_CMS13_cache[NumPar][i] );
2900 CacheShiftReal(ip_ex_pp_phi_ZZ_llqqnunull_CMS13_cache, NumPar, params, newResult);
2907 double params[] = {mass};
2909 int i =
CacheCheckReal(ip_ex_pp_phi_ZZ_qqnunu_CMS13_cache, NumPar, params);
2911 return ( ip_ex_pp_phi_ZZ_qqnunu_CMS13_cache[NumPar][i] );
2914 CacheShiftReal(ip_ex_pp_phi_ZZ_qqnunu_CMS13_cache, NumPar, params, newResult);
2921 double params[] = {mass};
2923 int i =
CacheCheckReal(ip_ex_gg_phi_WW_ATLAS8_cache, NumPar, params);
2925 return ( ip_ex_gg_phi_WW_ATLAS8_cache[NumPar][i] );
2928 CacheShiftReal(ip_ex_gg_phi_WW_ATLAS8_cache, NumPar, params, newResult);
2935 double params[] = {mass};
2937 int i =
CacheCheckReal(ip_ex_VV_phi_WW_ATLAS8_cache, NumPar, params);
2939 return ( ip_ex_VV_phi_WW_ATLAS8_cache[NumPar][i] );
2942 CacheShiftReal(ip_ex_VV_phi_WW_ATLAS8_cache, NumPar, params, newResult);
2949 double params[] = {mass};
2951 int i =
CacheCheckReal(ip_ex_gg_phi_WW_enumunu_ATLAS13_cache, NumPar, params);
2953 return ( ip_ex_gg_phi_WW_enumunu_ATLAS13_cache[NumPar][i] );
2956 CacheShiftReal(ip_ex_gg_phi_WW_enumunu_ATLAS13_cache, NumPar, params, newResult);
2963 double params[] = {mass};
2965 int i =
CacheCheckReal(ip_ex_VV_phi_WW_enumunu_ATLAS13_cache, NumPar, params);
2967 return ( ip_ex_VV_phi_WW_enumunu_ATLAS13_cache[NumPar][i] );
2970 CacheShiftReal(ip_ex_VV_phi_WW_enumunu_ATLAS13_cache, NumPar, params, newResult);
2977 double params[] = {mass};
2979 int i =
CacheCheckReal(ip_ex_ggVV_phi_WW_lnulnu_CMS13_cache, NumPar, params);
2981 return ( ip_ex_ggVV_phi_WW_lnulnu_CMS13_cache[NumPar][i] );
2984 CacheShiftReal(ip_ex_ggVV_phi_WW_lnulnu_CMS13_cache, NumPar, params, newResult);
2991 double params[] = {mass};
2993 int i =
CacheCheckReal(ip_ex_gg_phi_WW_lnuqq_ATLAS13_cache, NumPar, params);
2995 return ( ip_ex_gg_phi_WW_lnuqq_ATLAS13_cache[NumPar][i] );
2998 CacheShiftReal(ip_ex_gg_phi_WW_lnuqq_ATLAS13_cache, NumPar, params, newResult);
3005 double params[] = {mass};
3007 int i =
CacheCheckReal(ip_ex_VV_phi_WW_lnuqq_ATLAS13_cache, NumPar, params);
3009 return ( ip_ex_VV_phi_WW_lnuqq_ATLAS13_cache[NumPar][i] );
3012 CacheShiftReal(ip_ex_VV_phi_WW_lnuqq_ATLAS13_cache, NumPar, params, newResult);
3019 double params[] = {mass};
3021 int i =
CacheCheckReal(ip_ex_pp_phi_WW_lnuqq_CMS13_cache, NumPar, params);
3023 return ( ip_ex_pp_phi_WW_lnuqq_CMS13_cache[NumPar][i] );
3026 CacheShiftReal(ip_ex_pp_phi_WW_lnuqq_CMS13_cache, NumPar, params, newResult);
3033 double params[] = {mass};
3035 int i =
CacheCheckReal(ip_ex_mu_pp_phi_VV_CMS8_cache, NumPar, params);
3037 return ( ip_ex_mu_pp_phi_VV_CMS8_cache[NumPar][i] );
3040 CacheShiftReal(ip_ex_mu_pp_phi_VV_CMS8_cache, NumPar, params, newResult);
3047 double params[] = {mass};
3049 int i =
CacheCheckReal(ip_ex_pp_phi_VV_qqqq_ATLAS13_cache, NumPar, params);
3051 return ( ip_ex_pp_phi_VV_qqqq_ATLAS13_cache[NumPar][i] );
3054 CacheShiftReal(ip_ex_pp_phi_VV_qqqq_ATLAS13_cache, NumPar, params, newResult);
3061 double params[] = {mass};
3063 int i =
CacheCheckReal(ip_ex_gg_phi_phi1phi1_ATLAS8_cache, NumPar, params);
3065 return(ip_ex_gg_phi_phi1phi1_ATLAS8_cache[NumPar][i] );
3068 CacheShiftReal(ip_ex_gg_phi_phi1phi1_ATLAS8_cache, NumPar, params, newResult);
3075 double params[] = {mass};
3077 int i =
CacheCheckReal(ip_ex_pp_phi_phi1phi1_bbbb_CMS8_cache, NumPar, params);
3079 return(ip_ex_pp_phi_phi1phi1_bbbb_CMS8_cache[NumPar][i] );
3082 CacheShiftReal(ip_ex_pp_phi_phi1phi1_bbbb_CMS8_cache, NumPar, params, newResult);
3089 double params[] = {mass};
3091 int i =
CacheCheckReal(ip_ex_pp_phi_phi1phi1_bbgaga_CMS8_cache, NumPar, params);
3093 return(ip_ex_pp_phi_phi1phi1_bbgaga_CMS8_cache[NumPar][i] );
3096 CacheShiftReal(ip_ex_pp_phi_phi1phi1_bbgaga_CMS8_cache, NumPar, params, newResult);
3103 double params[] = {mass};
3105 int i =
CacheCheckReal(ip_ex_gg_phi_phi1phi1_bbtautau_CMS8_cache, NumPar, params);
3107 return(ip_ex_gg_phi_phi1phi1_bbtautau_CMS8_cache[NumPar][i] );
3110 CacheShiftReal(ip_ex_gg_phi_phi1phi1_bbtautau_CMS8_cache, NumPar, params, newResult);
3117 double params[] = {mass};
3119 int i =
CacheCheckReal(ip_ex_pp_phi_phi1phi1_bbtautau_CMS8_cache, NumPar, params);
3121 return(ip_ex_pp_phi_phi1phi1_bbtautau_CMS8_cache[NumPar][i] );
3124 CacheShiftReal(ip_ex_pp_phi_phi1phi1_bbtautau_CMS8_cache, NumPar, params, newResult);
3131 double params[] = {mass};
3133 int i =
CacheCheckReal(ip_ex_pp_phi_phi1phi1_bbbb_ATLAS13_cache, NumPar, params);
3135 return(ip_ex_pp_phi_phi1phi1_bbbb_ATLAS13_cache[NumPar][i] );
3138 CacheShiftReal(ip_ex_pp_phi_phi1phi1_bbbb_ATLAS13_cache, NumPar, params, newResult);
3145 double params[] = {mass};
3147 int i =
CacheCheckReal(ip_ex_pp_phi_phi1phi1_bbbb_1_CMS13_cache, NumPar, params);
3149 return(ip_ex_pp_phi_phi1phi1_bbbb_1_CMS13_cache[NumPar][i] );
3152 CacheShiftReal(ip_ex_pp_phi_phi1phi1_bbbb_1_CMS13_cache, NumPar, params, newResult);
3159 double params[] = {mass};
3161 int i =
CacheCheckReal(ip_ex_pp_phi_phi1phi1_bbbb_2_CMS13_cache, NumPar, params);
3163 return(ip_ex_pp_phi_phi1phi1_bbbb_2_CMS13_cache[NumPar][i] );
3166 CacheShiftReal(ip_ex_pp_phi_phi1phi1_bbbb_2_CMS13_cache, NumPar, params, newResult);
3173 double params[] = {mass};
3175 int i =
CacheCheckReal(ip_ex_pp_phi_phi1phi1_bbgaga_ATLAS13_cache, NumPar, params);
3177 return(ip_ex_pp_phi_phi1phi1_bbgaga_ATLAS13_cache[NumPar][i] );
3180 CacheShiftReal(ip_ex_pp_phi_phi1phi1_bbgaga_ATLAS13_cache, NumPar, params, newResult);
3187 double params[] = {mass};
3189 int i =
CacheCheckReal(ip_ex_pp_phi_phi1phi1_bbgaga_CMS13_cache, NumPar, params);
3191 return(ip_ex_pp_phi_phi1phi1_bbgaga_CMS13_cache[NumPar][i] );
3194 CacheShiftReal(ip_ex_pp_phi_phi1phi1_bbgaga_CMS13_cache, NumPar, params, newResult);
3201 double params[] = {mass};
3203 int i =
CacheCheckReal(ip_ex_pp_phi_phi1phi1_bbtautau_ATLAS13_cache, NumPar, params);
3205 return(ip_ex_pp_phi_phi1phi1_bbtautau_ATLAS13_cache[NumPar][i] );
3208 CacheShiftReal(ip_ex_pp_phi_phi1phi1_bbtautau_ATLAS13_cache, NumPar, params, newResult);
3215 double params[] = {mass};
3217 int i =
CacheCheckReal(ip_ex_pp_phi_phi1phi1_bbtautau_1_CMS13_cache, NumPar, params);
3219 return(ip_ex_pp_phi_phi1phi1_bbtautau_1_CMS13_cache[NumPar][i] );
3222 CacheShiftReal(ip_ex_pp_phi_phi1phi1_bbtautau_1_CMS13_cache, NumPar, params, newResult);
3229 double params[] = {mass};
3231 int i =
CacheCheckReal(ip_ex_pp_phi_phi1phi1_bbtautau_2_CMS13_cache, NumPar, params);
3233 return(ip_ex_pp_phi_phi1phi1_bbtautau_2_CMS13_cache[NumPar][i] );
3236 CacheShiftReal(ip_ex_pp_phi_phi1phi1_bbtautau_2_CMS13_cache, NumPar, params, newResult);
3243 double params[] = {mass};
3245 int i =
CacheCheckReal(ip_ex_pp_phi_phi1phi1_bbVV_CMS13_cache, NumPar, params);
3247 return(ip_ex_pp_phi_phi1phi1_bbVV_CMS13_cache[NumPar][i] );
3250 CacheShiftReal(ip_ex_pp_phi_phi1phi1_bbVV_CMS13_cache, NumPar, params, newResult);
3258 double params[] = {mass};
3260 int i =
CacheCheckReal(ip_ex_pp_phi_phi1phi1_bbWW_ATLAS13_cache, NumPar, params);
3262 return(ip_ex_pp_phi_phi1phi1_bbWW_ATLAS13_cache[NumPar][i] );
3265 CacheShiftReal(ip_ex_pp_phi_phi1phi1_bbWW_ATLAS13_cache, NumPar, params, newResult);
3272 double params[] = {mass};
3274 int i =
CacheCheckReal(ip_ex_gg_phi_phi1phi1_gagaWW_ATLAS13_cache, NumPar, params);
3276 return(ip_ex_gg_phi_phi1phi1_gagaWW_ATLAS13_cache[NumPar][i] );
3279 CacheShiftReal(ip_ex_gg_phi_phi1phi1_gagaWW_ATLAS13_cache, NumPar, params, newResult);
3286 double params[] = {mass};
3288 int i =
CacheCheckReal(ip_ex_gg_phi_phi1Z_bbZ_ATLAS8_cache, NumPar, params);
3290 return ( ip_ex_gg_phi_phi1Z_bbZ_ATLAS8_cache[NumPar][i] );
3293 CacheShiftReal(ip_ex_gg_phi_phi1Z_bbZ_ATLAS8_cache, NumPar, params, newResult);
3300 double params[] = {mass};
3302 int i =
CacheCheckReal(ip_ex_gg_phi_phi1Z_bbll_CMS8_cache, NumPar, params);
3304 return ( ip_ex_gg_phi_phi1Z_bbll_CMS8_cache[NumPar][i] );
3307 CacheShiftReal(ip_ex_gg_phi_phi1Z_bbll_CMS8_cache, NumPar, params, newResult);
3314 double params[] = {mass};
3316 int i =
CacheCheckReal(ip_ex_gg_phi_phi1Z_tautauZ_ATLAS8_cache, NumPar, params);
3318 return ( ip_ex_gg_phi_phi1Z_tautauZ_ATLAS8_cache[NumPar][i] );
3321 CacheShiftReal(ip_ex_gg_phi_phi1Z_tautauZ_ATLAS8_cache, NumPar, params, newResult);
3328 double params[] = {mass};
3330 int i =
CacheCheckReal(ip_ex_gg_phi_phi1Z_tautaull_CMS8_cache, NumPar, params);
3332 return ( ip_ex_gg_phi_phi1Z_tautaull_CMS8_cache[NumPar][i] );
3335 CacheShiftReal(ip_ex_gg_phi_phi1Z_tautaull_CMS8_cache, NumPar, params, newResult);
3342 double params[] = {mass};
3344 int i =
CacheCheckReal(ip_ex_gg_phi_phi1Z_bbZ_ATLAS13_cache, NumPar, params);
3346 return ( ip_ex_gg_phi_phi1Z_bbZ_ATLAS13_cache[NumPar][i] );
3349 CacheShiftReal(ip_ex_gg_phi_phi1Z_bbZ_ATLAS13_cache, NumPar, params, newResult);
3356 double params[] = {mass};
3358 int i =
CacheCheckReal(ip_ex_gg_phi_phi1Z_bbZ_1_CMS13_cache, NumPar, params);
3360 return ( ip_ex_gg_phi_phi1Z_bbZ_1_CMS13_cache[NumPar][i] );
3363 CacheShiftReal(ip_ex_gg_phi_phi1Z_bbZ_1_CMS13_cache, NumPar, params, newResult);
3370 double params[] = {mass};
3372 int i =
CacheCheckReal(ip_ex_gg_phi_phi1Z_bbZ_2_CMS13_cache, NumPar, params);
3374 return ( ip_ex_gg_phi_phi1Z_bbZ_2_CMS13_cache[NumPar][i] );
3377 CacheShiftReal(ip_ex_gg_phi_phi1Z_bbZ_2_CMS13_cache, NumPar, params, newResult);
3384 double params[] = {mass};
3386 int i =
CacheCheckReal(ip_ex_bb_phi_phi1Z_bbZ_ATLAS13_cache, NumPar, params);
3388 return ( ip_ex_bb_phi_phi1Z_bbZ_ATLAS13_cache[NumPar][i] );
3391 CacheShiftReal(ip_ex_bb_phi_phi1Z_bbZ_ATLAS13_cache, NumPar, params, newResult);
3398 double params[] = {mass};
3400 int i =
CacheCheckReal(ip_ex_bb_phi_phi1Z_bbZ_1_CMS13_cache, NumPar, params);
3402 return ( ip_ex_bb_phi_phi1Z_bbZ_1_CMS13_cache[NumPar][i] );
3405 CacheShiftReal(ip_ex_bb_phi_phi1Z_bbZ_1_CMS13_cache, NumPar, params, newResult);
3412 double params[] = {mass};
3414 int i =
CacheCheckReal(ip_ex_bb_phi_phi1Z_bbZ_2_CMS13_cache, NumPar, params);
3416 return ( ip_ex_bb_phi_phi1Z_bbZ_2_CMS13_cache[NumPar][i] );
3419 CacheShiftReal(ip_ex_bb_phi_phi1Z_bbZ_2_CMS13_cache, NumPar, params, newResult);
3426 double params[] = {mi, mj};
3428 int i =
CacheCheckReal(ip_ex_pp_phii_phijZ_bbll_1_CMS8_cache, NumPar, params);
3430 return ( ip_ex_pp_phii_phijZ_bbll_1_CMS8_cache[NumPar][i] );
3433 CacheShiftReal(ip_ex_pp_phii_phijZ_bbll_1_CMS8_cache, NumPar, params, newResult);
3442 double params[] = {mi, mj};
3444 int i =
CacheCheckReal(ip_ex_pp_phii_phijZ_bbll_2_CMS8_cache, NumPar, params);
3446 return ( ip_ex_pp_phii_phijZ_bbll_2_CMS8_cache[NumPar][i] );
3449 CacheShiftReal(ip_ex_pp_phii_phijZ_bbll_2_CMS8_cache, NumPar, params, newResult);
3457 double params[] = {mi, mj};
3459 int i =
CacheCheckReal(ip_ex_pp_phii_phijZ_tautaull_1_CMS8_cache, NumPar, params);
3461 return ( ip_ex_pp_phii_phijZ_tautaull_1_CMS8_cache[NumPar][i] );
3464 CacheShiftReal(ip_ex_pp_phii_phijZ_tautaull_1_CMS8_cache, NumPar, params, newResult);
3473 double params[] = {mi, mj};
3475 int i =
CacheCheckReal(ip_ex_pp_phii_phijZ_tautaull_2_CMS8_cache, NumPar, params);
3477 return ( ip_ex_pp_phii_phijZ_tautaull_2_CMS8_cache[NumPar][i] );
3480 CacheShiftReal(ip_ex_pp_phii_phijZ_tautaull_2_CMS8_cache, NumPar, params, newResult);
3489 double params[] = {mj, mi};
3491 int i =
CacheCheckReal(ip_ex_gg_phii_phijZ_bbZ_ATLAS13_cache, NumPar, params);
3493 return ( ip_ex_gg_phii_phijZ_bbZ_ATLAS13_cache[NumPar][i] );
3496 CacheShiftReal(ip_ex_gg_phii_phijZ_bbZ_ATLAS13_cache, NumPar, params, newResult);
3504 double params[] = {mj, mi};
3506 int i =
CacheCheckReal(ip_ex_bb_phii_phijZ_bbZ_ATLAS13_cache, NumPar, params);
3508 return ( ip_ex_bb_phii_phijZ_bbZ_ATLAS13_cache[NumPar][i] );
3511 CacheShiftReal(ip_ex_bb_phii_phijZ_bbZ_ATLAS13_cache, NumPar, params, newResult);
5120 double params[] = {mass};
5122 int i =
CacheCheckReal(ip_ex_pp_Hpm_taunu_ATLAS8_cache, NumPar, params);
5124 return(ip_ex_pp_Hpm_taunu_ATLAS8_cache[NumPar][i] );
5127 CacheShiftReal(ip_ex_pp_Hpm_taunu_ATLAS8_cache, NumPar, params, newResult);
5136 double params[] = {mass};
5138 int i =
CacheCheckReal(ip_ex_pp_Hp_taunu_CMS8_cache, NumPar, params);
5140 return(ip_ex_pp_Hp_taunu_CMS8_cache[NumPar][i] );
5143 CacheShiftReal(ip_ex_pp_Hp_taunu_CMS8_cache, NumPar, params, newResult);
5152 double params[] = {mass};
5154 int i =
CacheCheckReal(ip_ex_pp_Hpm_taunu_ATLAS13_cache, NumPar, params);
5156 return(ip_ex_pp_Hpm_taunu_ATLAS13_cache[NumPar][i] );
5159 CacheShiftReal(ip_ex_pp_Hpm_taunu_ATLAS13_cache, NumPar, params, newResult);
5168 double params[] = {mass};
5170 int i =
CacheCheckReal(ip_ex_pp_Hpm_taunu_CMS13_cache, NumPar, params);
5172 return(ip_ex_pp_Hpm_taunu_CMS13_cache[NumPar][i] );
5175 CacheShiftReal(ip_ex_pp_Hpm_taunu_CMS13_cache, NumPar, params, newResult);
5184 double params[] = {mass};
5186 int i =
CacheCheckReal(ip_ex_pp_Hpm_tb_ATLAS8_cache, NumPar, params);
5188 return(ip_ex_pp_Hpm_tb_ATLAS8_cache[NumPar][i] );
5191 CacheShiftReal(ip_ex_pp_Hpm_tb_ATLAS8_cache, NumPar, params, newResult);
5200 double params[] = {mass};
5202 int i =
CacheCheckReal(ip_ex_pp_Hp_tb_CMS8_cache, NumPar, params);
5204 return(ip_ex_pp_Hp_tb_CMS8_cache[NumPar][i] );
5207 CacheShiftReal(ip_ex_pp_Hp_tb_CMS8_cache, NumPar, params, newResult);
5214 double params[] = {mass};
5216 int i =
CacheCheckReal(ip_ex_pp_Hpm_tb_ATLAS13_cache, NumPar, params);
5218 return(ip_ex_pp_Hpm_tb_ATLAS13_cache[NumPar][i] );
5221 CacheShiftReal(ip_ex_pp_Hpm_tb_ATLAS13_cache, NumPar, params, newResult);
5228 double params[] = {logtb, logmHp};
5232 return ( ip_ex_bsgamma_cache[NumPar][i] );
5248 return pow(asin(
sqrt(1.0/x)),2);
5264 gslpp::complex GeneralTHDMcache::Int1(
const double tau,
const double lambda)
const{
5265 return tau*lambda/(tau-lambda)/2.0+tau*tau*lambda*lambda/((tau-lambda)
5266 *(tau-lambda))/2.0*(f_func(tau)-f_func(lambda))+tau*tau*lambda/((tau-lambda)
5267 *(tau-lambda))*(g_func(tau)-g_func(lambda));
5270 gslpp::complex GeneralTHDMcache::Int2(
const double tau,
const double lambda)
const{
5271 return -tau*lambda/(tau-lambda)/2.0*(f_func(tau)-f_func(lambda));
5279 double params[] = {mHl2, Mu, Mc, Mt};
5281 int i =
CacheCheck(I_h_U_cache, NumPar, params);
5283 return ( I_h_U_cache[NumPar][i] );
5285 double TAUu=4.0*Mu*Mu/mHl2;
5286 double TAUc=4.0*Mc*Mc/mHl2;
5287 double TAUt=4.0*Mt*Mt/mHl2;
5288 gslpp::complex newResult = -(8./3.)*(TAUu*(1.0+(1.0-TAUu)*f_func(TAUu))
5289 +TAUc*(1.0+(1.0-TAUc)*f_func(TAUc))+TAUt*(1.0+(1.0-TAUt)*f_func(TAUt)));
5290 CacheShift(I_h_U_cache, NumPar, params, newResult);
5297 double params[] = {mHh2, Mc, Mt};
5299 int i =
CacheCheck(I_HH_U_cache, NumPar, params);
5301 return ( I_HH_U_cache[NumPar][i] );
5303 double TAUc=4.0*Mc*Mc/mHh2;
5304 double TAUt=4.0*Mt*Mt/mHh2;
5305 gslpp::complex newResult = -(8./3.)*(TAUc*(1.0+(1.0-TAUc)*f_func(TAUc))
5306 +TAUt*(1.0+(1.0-TAUt)*f_func(TAUt)));
5307 CacheShift(I_HH_U_cache, NumPar, params, newResult);
5314 double params[] = {mA2, Mc, Mt};
5316 int i =
CacheCheck(I_A_U_cache, NumPar, params);
5318 return ( I_A_U_cache[NumPar][i] );
5320 double TAUc=4.0*Mc*Mc/mA2;
5321 double TAUt=4.0*Mt*Mt/mA2;
5322 gslpp::complex newResult = -(8./3.)*(TAUc*f_func(TAUc)+TAUt*f_func(TAUt));
5323 CacheShift(I_A_U_cache, NumPar, params, newResult);
5330 double params[] = {mHl2, Md, Ms, Mb};
5332 int i =
CacheCheck(I_h_D_cache, NumPar, params);
5334 return ( I_h_D_cache[NumPar][i] );
5336 double TAUd=4.0*Md*Md/mHl2;
5337 double TAUs=4.0*Ms*Ms/mHl2;
5338 double TAUb=4.0*Mb*Mb/mHl2;
5339 gslpp::complex newResult = -(2./3.)*(TAUd*(1.0+(1.0-TAUd)*f_func(TAUd))
5340 +TAUs*(1.0+(1.0-TAUs)*f_func(TAUs))+TAUb*(1.0+(1.0-TAUb)*f_func(TAUb)));
5341 CacheShift(I_h_D_cache, NumPar, params, newResult);
5348 double params[] = {mHh2, Ms, Mb};
5350 int i =
CacheCheck(I_HH_D_cache, NumPar, params);
5352 return ( I_HH_D_cache[NumPar][i] );
5354 double TAUs=4.0*Ms*Ms/mHh2;
5355 double TAUb=4.0*Mb*Mb/mHh2;
5356 gslpp::complex newResult = -(2./3.)*(TAUs*(1.0+(1.0-TAUs)*f_func(TAUs))
5357 +TAUb*(1.0+(1.0-TAUb)*f_func(TAUb)));
5358 CacheShift(I_HH_D_cache, NumPar, params, newResult);
5365 double params[] = {mA2, Ms, Mb};
5367 int i =
CacheCheck(I_A_D_cache, NumPar, params);
5369 return ( I_A_D_cache[NumPar][i] );
5371 double TAUs=4.0*Ms*Ms/mA2;
5372 double TAUb=4.0*Mb*Mb/mA2;
5373 gslpp::complex newResult = -(2./3.)*(TAUs*f_func(TAUs)+TAUb*f_func(TAUb));
5374 CacheShift(I_A_D_cache, NumPar, params, newResult);
5381 double params[] = {mHl2, Me, Mmu, Mtau};
5383 int i =
CacheCheck(I_h_L_cache, NumPar, params);
5385 return ( I_h_L_cache[NumPar][i] );
5387 double TAUe=4.0*Me*Me/mHl2;
5388 double TAUmu=4.0*Mmu*Mmu/mHl2;
5389 double TAUtau=4.0*Mtau*Mtau/mHl2;
5390 gslpp::complex newResult = -2.0*(TAUe*(1.0+(1.0-TAUe)*f_func(TAUe))
5391 +TAUmu*(1.0+(1.0-TAUmu)*f_func(TAUmu))
5392 +TAUtau*(1.0+(1.0-TAUtau)*f_func(TAUtau)));
5393 CacheShift(I_h_L_cache, NumPar, params, newResult);
5400 double params[] = {mHh2, Mmu, Mtau};
5402 int i =
CacheCheck(I_HH_L_cache, NumPar, params);
5404 return ( I_HH_L_cache[NumPar][i] );
5406 double TAUmu=4.0*Mmu*Mmu/mHh2;
5407 double TAUtau=4.0*Mtau*Mtau/mHh2;
5408 gslpp::complex newResult = -2.0*(TAUmu*(1.0+(1.0-TAUmu)*f_func(TAUmu))+
5409 TAUtau*(1.0+(1.0-TAUtau)*f_func(TAUtau)));
5410 CacheShift(I_HH_L_cache, NumPar, params, newResult);
5417 double params[] = {mA2, Mmu, Mtau};
5419 int i =
CacheCheck(I_A_L_cache, NumPar, params);
5421 return ( I_A_L_cache[NumPar][i] );
5423 double TAUmu=4.0*Mmu*Mmu/mA2;
5424 double TAUtau=4.0*Mtau*Mtau/mA2;
5425 gslpp::complex newResult = -2.0*(TAUmu*f_func(TAUmu)+TAUtau*f_func(TAUtau));
5426 CacheShift(I_A_L_cache, NumPar, params, newResult);
5433 double params[] = {mH, MW};
5435 int i =
CacheCheck(I_H_W_cache, NumPar, params);
5437 return ( I_H_W_cache[NumPar][i] );
5439 double TAUw=4.0*MW*MW/(mH*mH);
5440 gslpp::complex newResult = 2.0 + 3.0*TAUw + 3.0*TAUw*(2.0-TAUw)*f_func(TAUw);
5441 CacheShift(I_H_W_cache, NumPar, params, newResult);
5448 double params[] = {mHp2, mH};
5450 int i =
CacheCheck(I_H_Hp_cache, NumPar, params);
5452 return ( I_H_Hp_cache[NumPar][i] );
5454 double TAUhp=4.0*mHp2/(mH*mH);
5456 CacheShift(I_H_Hp_cache, NumPar, params, newResult);
5463 double params[] = {mHl2, cW2, Mu, Mc, Mt, MZ};
5465 int i =
CacheCheck(A_h_U_cache, NumPar, params);
5467 return ( A_h_U_cache[NumPar][i] );
5469 double TAUu=4.0*Mu*Mu/mHl2;
5470 double TAUc=4.0*Mc*Mc/mHl2;
5471 double TAUt=4.0*Mt*Mt/mHl2;
5472 double LAMu=4.0*Mu*Mu/(MZ*MZ);
5473 double LAMc=4.0*Mc*Mc/(MZ*MZ);
5474 double LAMt=4.0*Mt*Mt/(MZ*MZ);
5476 gslpp::complex newResult = -4.0*(1.0/2.0-4.0/3.0*sW2)*(Int1(TAUu,LAMu)+Int1(TAUc,LAMc)
5477 +Int1(TAUt,LAMt)-Int2(TAUu,LAMu)-Int2(TAUc,LAMc)-Int2(TAUt,LAMt));
5478 CacheShift(A_h_U_cache, NumPar, params, newResult);
5486 double params[] = {mHh2, cW2, Mc, Mt, MZ};
5488 int i =
CacheCheck(A_HH_U_cache, NumPar, params);
5490 return ( A_HH_U_cache[NumPar][i] );
5492 double TAUc=4.0*Mc*Mc/mHh2;
5493 double TAUt=4.0*Mt*Mt/mHh2;
5494 double LAMc=4.0*Mc*Mc/(MZ*MZ);
5495 double LAMt=4.0*Mt*Mt/(MZ*MZ);
5497 gslpp::complex newResult = -4.0*(1.0/2.0-4.0/3.0*sW2)*(Int1(TAUc,LAMc)-Int2(TAUc,LAMc)
5498 +Int1(TAUt,LAMt)-Int2(TAUt,LAMt));
5499 CacheShift(A_HH_U_cache, NumPar, params, newResult);
5506 double params[] = {mA2, cW2, Mc, Mt, MZ};
5508 int i =
CacheCheck(A_A_U_cache, NumPar, params);
5510 return ( A_A_U_cache[NumPar][i] );
5512 double TAUc=4.0*Mc*Mc/mA2;
5513 double TAUt=4.0*Mt*Mt/mA2;
5514 double LAMc=4.0*Mc*Mc/(MZ*MZ);
5515 double LAMt=4.0*Mt*Mt/(MZ*MZ);
5517 gslpp::complex newResult = -4.0*(1.0/2.0-4.0/3.0*sW2)*(-Int2(TAUc,LAMc)-Int2(TAUt,LAMt))/
sqrt(sW2*cW2);
5518 CacheShift(A_A_U_cache, NumPar, params, newResult);
5525 double params[] = {mHl2, cW2, Md, Ms, Mb, MZ};
5527 int i =
CacheCheck(A_h_D_cache, NumPar, params);
5529 return ( A_h_D_cache[NumPar][i] );
5531 double TAUd=4.0*Md*Md/mHl2;
5532 double TAUs=4.0*Ms*Ms/mHl2;
5533 double TAUb=4.0*Mb*Mb/mHl2;
5534 double LAMd=4.0*Md*Md/(MZ*MZ);
5535 double LAMs=4.0*Ms*Ms/(MZ*MZ);
5536 double LAMb=4.0*Mb*Mb/(MZ*MZ);
5538 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0/3.0*sW2)*(Int1(TAUd,LAMd)+Int1(TAUs,LAMs)
5539 +Int1(TAUb,LAMb)-Int2(TAUd,LAMd)-Int2(TAUs,LAMs)-Int2(TAUb,LAMb));
5540 CacheShift(A_h_D_cache, NumPar, params, newResult);
5547 double params[] = {mHh2, cW2, Ms, Mb, MZ};
5549 int i =
CacheCheck(A_HH_D_cache, NumPar, params);
5551 return ( A_HH_D_cache[NumPar][i] );
5553 double TAUs=4.0*Ms*Ms/mHh2;
5554 double TAUb=4.0*Mb*Mb/mHh2;
5555 double LAMs=4.0*Ms*Ms/(MZ*MZ);
5556 double LAMb=4.0*Mb*Mb/(MZ*MZ);
5558 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0/3.0*sW2)*(Int1(TAUs,LAMs)-Int2(TAUs,LAMs)
5559 +Int1(TAUb,LAMb)-Int2(TAUb,LAMb));
5560 CacheShift(A_HH_D_cache, NumPar, params, newResult);
5567 double params[] = {mA2, cW2, Ms, Mb, MZ};
5569 int i =
CacheCheck(A_A_D_cache, NumPar, params);
5571 return ( A_A_D_cache[NumPar][i] );
5573 double TAUs=4.0*Ms*Ms/mA2;
5574 double TAUb=4.0*Mb*Mb/mA2;
5575 double LAMs=4.0*Ms*Ms/(MZ*MZ);
5576 double LAMb=4.0*Mb*Mb/(MZ*MZ);
5578 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0/3.0*sW2)*(-Int2(TAUs,LAMs)-Int2(TAUb,LAMb))/
sqrt(sW2*cW2);
5579 CacheShift(A_A_D_cache, NumPar, params, newResult);
5586 double params[] = {mHl2, cW2, Me, Mmu, Mtau, MZ};
5588 int i =
CacheCheck(A_h_L_cache, NumPar, params);
5590 return ( A_h_L_cache[NumPar][i] );
5592 double TAUe=4.0*Me*Me/mHl2;
5593 double TAUmu=4.0*Mmu*Mmu/mHl2;
5594 double TAUtau=4.0*Mtau*Mtau/mHl2;
5595 double LAMe=4.0*Me*Me/(MZ*MZ);
5596 double LAMmu=4.0*Mmu*Mmu/(MZ*MZ);
5597 double LAMtau=4.0*Mtau*Mtau/(MZ*MZ);
5599 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0*sW2)*(Int1(TAUe,LAMe)+Int1(TAUmu,LAMmu)
5600 +Int1(TAUtau,LAMtau)-Int2(TAUe,LAMe)-Int2(TAUmu,LAMmu)
5601 -Int2(TAUtau,LAMtau));
5602 CacheShift(A_h_L_cache, NumPar, params, newResult);
5609 double params[] = {mHh2, cW2, Mmu, Mtau, MZ};
5611 int i =
CacheCheck(A_HH_L_cache, NumPar, params);
5613 return ( A_HH_L_cache[NumPar][i] );
5615 double TAUmu=4.0*Mmu*Mmu/mHh2;
5616 double TAUtau=4.0*Mtau*Mtau/mHh2;
5617 double LAMmu=4.0*Mmu*Mmu/(MZ*MZ);
5618 double LAMtau=4.0*Mtau*Mtau/(MZ*MZ);
5620 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0*sW2)*(Int1(TAUmu,LAMmu)-Int2(TAUmu,LAMmu)
5621 +Int1(TAUtau,LAMtau)-Int2(TAUtau,LAMtau));
5622 CacheShift(A_HH_L_cache, NumPar, params, newResult);
5631 double params[] = {mA2, cW2, Mmu, Mtau, MZ};
5633 int i =
CacheCheck(A_A_L_cache, NumPar, params);
5635 return ( A_A_L_cache[NumPar][i] );
5637 double TAUmu=4.0*Mmu*Mmu/mA2;
5638 double TAUtau=4.0*Mtau*Mtau/mA2;
5639 double LAMmu=4.0*Mmu*Mmu/(MZ*MZ);
5640 double LAMtau=4.0*Mtau*Mtau/(MZ*MZ);
5642 gslpp::complex newResult = 2.0*(-1.0/2.0+2.0*sW2)*(-Int2(TAUmu,LAMmu)-Int2(TAUtau,LAMtau))/
sqrt(sW2*cW2);
5643 CacheShift(A_A_L_cache, NumPar, params, newResult);
5650 double params[] = {mH, cW2, MW, MZ};
5652 int i =
CacheCheck(A_H_W_cache, NumPar, params);
5654 return ( A_H_W_cache[NumPar][i] );
5656 double TAUw=4.0*MW*MW/(mH*mH);
5657 double LAMw=4.0*MW*MW/(MZ*MZ);
5660 +((1.0+2.0/TAUw)*sW2/cW2-(5.0+2.0/TAUw))*Int1(TAUw,LAMw));
5661 CacheShift(A_H_W_cache, NumPar, params, newResult);
5668 double params[] = {mHp2, mH, cW2, MZ};
5670 int i =
CacheCheck(A_H_Hp_cache, NumPar, params);
5672 return ( A_H_Hp_cache[NumPar][i] );
5674 double TAUhp=4.0*mHp2/(mH*mH);
5675 double LAMhp=4.0*mHp2/(MZ*MZ);
5678 CacheShift(A_H_Hp_cache, NumPar, params, newResult);
5683 int GeneralTHDMcache::HSTheta (
const double x)
const{
5691 double params[] = {a2, b2, c2};
5695 return ( KaellenFunction_cache[NumPar][i] );
5698 double newResult = 0.0;
5699 double x = (a2-b2-c2)*(a2-b2-c2)-4.0*b2*c2;
5700 if(x>0) newResult =
sqrt(std::fabs(x/a2))/2.0;
5701 CacheShiftReal(KaellenFunction_cache, NumPar, params, newResult);
5721 return sqrt(1.0-4.0*mf*mf/m_2);
5727 if (4.0*mt*mt/m_2 < 1 )
5729 return 1.0-4.0*mt*mt/m_2;
5733 return -(1.0-4.0*mt*mt/m_2);
5738 double GeneralTHDMcache::lambdaijk(
const double Ri1,
const double Ri2,
const double Ri3,
const double Rj1,
const double Rj2,
const double Rj3,
const double Rk1,
const double Rk2,
const double Rk3,
const double lambda1H,
const double lambda3H,
const double lambda4H,
const double Relambda5H,
const double Imlambda5H,
const double Relambda6H,
const double Imlambda6H,
const double Relambda7H,
const double Imlambda7H)
const
5740 return (1.0/2.0)*vev*(Imlambda7H*(-Ri3*Rj3*Rk3 - Ri2*Rj2*Rk3) - 3.0*Imlambda6H*Ri1*Rj1*Rk3
5741 + lambda1H*Ri1*Rj1*Rk1 + Relambda7H*Ri2*Rj2*Rk2 + 3.0*Relambda6H*Ri1*Rj1*Rk2
5742 +(Relambda5H + lambda3H + lambda4H)*Ri1*Rj2*Rk2 - (2.0*Relambda5H - lambda3H - lambda4H)*Ri1*Rj3*Rk3
5743 + Relambda7H*Ri2*Rj3*Rk3 - Imlambda5H*Ri1*Rj2*Rk3);
5748 return -vev*(lambda3H*Ri1 + Relambda7H*Ri2 - Imlambda7H*Ri3);
5760 double GF = 1/(
sqrt(2.0)*vev*vev);
5765 if (!myGTHDM->getATHDMflag())
5767 throw std::runtime_error(
"Signal strengths are only available in the A2HDM.");
5773 Mt=myGTHDM->getQuarks(
QCD::TOP).getMass();
5779 Mu=myGTHDM->getQuarks(
QCD::UP).getMass();
5780 Md=myGTHDM->getQuarks(
QCD::DOWN).getMass();
5782 MW=
MWGTHDM(myGTHDM->Mw_tree());
5789 double m1_2 = mH1sq;
5798 yu1R = myGTHDM->getyu1R();
5799 yd1R = myGTHDM->getyd1R();
5800 yl1R = myGTHDM->getyl1R();
5808 double BrSM_htoWW = 2.15e-1;
5809 double BrSM_htoZZ = 2.64e-2;
5810 double BrSM_htogg = 8.57e-2;
5811 double BrSM_htoZga = 1.54e-3;
5812 double BrSM_htocc = 2.91e-2;
5815 double SigmaggF8 = myGTHDM->computeSigmaggH(8.0);
5817 double SigmaggF13 = myGTHDM->computeSigmaggH(13.0);
5819 double Sigmaggh_tt8 = myGTHDM->computeSigmaggH_tt(8.0);
5823 double Sigmaggh_bb8 = myGTHDM->computeSigmaggH_bb(8.0);
5827 double Sigmatth8 = myGTHDM->computeSigmattH(8.0);
5829 double Sigmatth13 = myGTHDM->computeSigmattH(13.0);
5833 double SigmaVBFVh13 = (myGTHDM->computeSigmaVBF(13.0)+myGTHDM->computeSigmaWH(13.0)+myGTHDM->computeSigmaZH(13.0));
5852 beta_h_t =
sqrt(-1.0+4.0*Mt*Mt/
m1_2);
5871 rh_gg = rh_ggE+rh_ggO*(SigmagghO_8/SigmaggF8);
5903 rh_gaga = ((I_hE_F+I_h_W+I_h_Hp).abs2()+ (I_hO_F).abs2())/(I_hSM_F +I_hSM_W).
abs2();
5927 rh_Zga = ((A_hE_F+A_h_W+A_h_Hp).abs2()+ A_hO_F.
abs2())/(A_hSM_F +A_hSM_W ).
abs2();
5931 double Gamma_hggSM=GF*Als*Als*m1*m1*m1/(
sqrt(2.0)*16.0*M_PI*M_PI*M_PI)*(9.0/4.0)*(fermU/4.0+fermD).abs2();
5933 double Gamma_hgg=rh_gg*GF*Als*Als*m1*m1*m1/(
sqrt(2.0)*16.0*M_PI*M_PI*M_PI)*(9.0/4.0)*(fermU/4.0+fermD).abs2();
5934 double lambda122 = (2.0)*(
lambdaijk(R11, R21, R31, R12, R22, R32, R11, R21, R31, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R12, R22, R32, R11, R21, R31, R11, R21, R31, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R12, R22, R32, R11, R21, R31, R11, R21, R31, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) );
5937 double lambda133 = (vev)*(Relambda7H*R21 - (2.0*Relambda5H - lambda3H - lambda4H)*R11);
5951 ggF_tth8 = (SigmaggF8*rh_gg + Sigmatth8*(rh_QuQuE + rh_QuQuO/(
beta(Mc,
m1_2)*
beta(Mc,
m1_2))))/(SigmaggF8 + Sigmatth8);
5953 ggF_tth13 = (SigmaggF13*rh_gg + Sigmatth8*(rh_QuQuE + rh_QuQuO/(
beta(Mc,
m1_2)*
beta(Mc,
m1_2))))/(SigmaggF13 + Sigmatth13);
5955 pph13 = (SigmaggF13*rh_gg+ SigmaVBFVh13*rh_VV + Sigmatth13*(rh_QuQuE + rh_QuQuO/(
beta(Mc,
m1_2)*
beta(Mc,
m1_2))) + Sigmabbh13*(rh_QdQdE + rh_QdQdE/(
beta(Mb,
m1_2)*
beta(Mb,
m1_2))))/(SigmaggF13 + SigmaVBFVh13 + Sigmatth13 + Sigmabbh13);
5960 + (BrSM_htoWW+BrSM_htoZZ)*rh_VV
5964 + BrSM_htoZga*rh_Zga
5966 Gamma_h = sumModBRs*myGTHDM->computeGammaHTotal() + Gamma_hHH + Gamma_hAA;
5969 GTHDM_BR_h_WW = rh_VV*BrSM_htoWW/sumModBRs;
5970 GTHDM_BR_h_ZZ = rh_VV*BrSM_htoZZ/sumModBRs;
5972 GTHDM_BR_h_cc =(rh_QuQuE + rh_QuQuO/(
beta(Mc,
m1_2)*
beta(Mc,
m1_2)))*BrSM_htocc/sumModBRs;
5974 GTHDM_BR_h_gg = (Gamma_hgg/Gamma_hggSM)*BrSM_htogg/sumModBRs;
5986 double GF=1/(
sqrt(2.0)*vev*vev);
5990 if (!myGTHDM->getATHDMflag())
5992 throw std::runtime_error(
"Direct Searches are only aviable in the A2HDM.");
6036 double rphi2_QuQuE= yu2.
real()*yu2.
real();
6037 double rphi2_QuQuO= yu2.
imag()*yu2.
imag();
6038 double rphi2_QdQdE= yd2.
real()*yd2.
real();
6039 double rphi2_QdQdO= yd2.
imag()*yd2.
imag();
6040 double rphi2_QlQlE= yl2.real()*yl2.real();
6041 double rphi2_QlQlO= yl2.imag()*yl2.imag();
6062 double lambdaphi2HpHm =
lambdaipm(R21, R22, R32);
6074 double Gamma_phi2gaga=(GF*Ale*Ale*m2*m2*m2/(
sqrt(2.0)*128.0*M_PI*M_PI*M_PI))*((I_phi2E_F+I_phi2_W+I_phi2_Hp).abs2()
6075 + (I_phi2O_F).abs2());
6100 double Gamma_phi2Zga=HSTheta(m2-MZ)*GF*Ale*Ale*m2*m2*m2/(
sqrt(2.0)*64.0*M_PI*M_PI*M_PI)*(1.0-MZ*MZ/(m2*m2))*(1.0-MZ*MZ/(m2*m2))*(1.0-MZ*MZ/(m2*m2))*((A_phi2E_F+A_phi2_W+A_phi2_Hp).abs2()
6101 + A_phi2O_F.
abs2());
6106 double Gamma_phi2gg=(rphi2_ggE)*GF*Als*Als*m2*m2*m2/(
sqrt(2.0)*16.0*M_PI*M_PI*M_PI)*(9.0/4.0)*(I_HH2_Ux/4.0+I_HH2_Dx).abs2()
6107 +rphi2_ggO*GF*Als*Als*m2*m2*m2/(
sqrt(2.0)*16.0*M_PI*M_PI*M_PI)*(9.0/4.0)*(I_A2_Ux/4.0+I_A2_Dx).abs2();
6128 SigmaTotSM_phi2_8 = 1.0e-15;
6130 if (m2>=20. && m2 <=2000.) {
6133 SigmaSumphi2_8 = SigmaggF_phi2_8 + SigmaVBF_phi2_8 + SigmaVH_phi2_8 + SigmattF_phi2_8 + SigmabbF_phi2_8;
6155 SigmaSumphi2_13 = SigmaggF_phi2_13 + SigmaVBF_phi2_13 + SigmaVH_phi2_13 + SigmattF_phi2_13 + SigmabbF_phi2_13;
6168 double lambda132 =
lambdaijk(R11, R21, R31, R13, R23, R33, R12, R22, R32, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R12, R22, R32, R11, R22, R32, R12, R22, R33, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R12, R22, R32, R13, R23, R33, R11, R21, R33, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R11, R21, R31, R13, R23, R33, R12, R22, R32, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R13, R23,R33, R11, R21, R31, R12, R22, R32,lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R13, R23,R33, R12, R22, R32, R11, R21, R31,lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) ;
6169 double lambda332 =
lambdaijk(R13, R23, R33, R13, R23, R33, R12, R22, R32, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R12, R22, R32, R13, R23, R33, R12, R22, R32, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R13, R23, R33, R12, R22, R32, R12, R22, R32, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) ;
6172 double Gammaphi2_phi3phi3=HSTheta(m2 - 2.0*m3)*
sqrt(std::fabs(1.0 - (4.0*
m3_2)/
m2_2))*lambda332*lambda332/(32.0*m2*M_PI);
6177 double lambda112 = (2.0)*(
lambdaijk(R11, R21, R31, R12, R22, R32, R11, R21, R31, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R11, R21, R31, R11, R21, R31, R12, R22, R32, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R12, R22, R32, R11, R21, R31, R11, R21, R31, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) );
6181 double Gammaphi2_phi1phi1=HSTheta(m2 - 2.0*m1)*
sqrt(std::fabs(1.0 - (4.0*
m1_2)/
m2_2))*lambda112*lambda112/(32.0*m2*M_PI);
6184 double Gammaphi2_HpHm=HSTheta(m2 - 2.0*
sqrt(mHp2))*
sqrt(std::fabs(1.0 - (4.0*mHp2)/
m2_2))*lambdaphi2HpHm*lambdaphi2HpHm/(8.0*m2*M_PI);
6186 double Gammaphi2_phi1Z=HSTheta(m2-(m1+MZ))*
pow(
KaellenFunction(
m2_2,MZ*MZ,
m1_2),3)*(R23*R12 + R22*R13)*(R23*R12 + R22*R13)/(2.0*M_PI*vev*vev);
6189 double Gammaphi2_phi3Z=HSTheta(m2-(m3+MZ))*
pow(
KaellenFunction(
m2_2,MZ*MZ,
m3_2),3)*(R33*R22 + R32*R23)*(R33*R22 + R32*R23)/(2.0*M_PI*vev*vev);
6195 Gammaphi2tot = 1.e-10;
6197 Gammaphi2tot= Gammaphi2tot + (BrSM_phi2tott*(rphi2_QuQuE + rphi2_QuQuO/(
beta_mt_sq(Mt,
m2_2)))
6200 +BrSM_phi2totautau*(rphi2_QlQlE + rphi2_QlQlO/(
beta(Mtau,
m2_2)*
beta(Mtau,
m2_2)))
6201 +(BrSM_phi2toWW+BrSM_phi2toZZ)*rphi2_VV)*Gammaphi2totSM+Gamma_phi2gaga+Gamma_phi2Zga+Gamma_phi2gg
6202 +Gammaphi2_phi1phi3+ Gammaphi2_phi1phi1+Br_phi2tophi3phi3+Gammaphi2_HpHm+Gammaphi2_phi1Z+Gammaphi2_phi3Z
6206 Br_phi2tott=BrSM_phi2tott*(rphi2_QuQuE + rphi2_QuQuO/(
beta_mt_sq(Mt,
m2_2)))*Gammaphi2totSM/Gammaphi2tot;
6207 Br_phi2tobb=BrSM_phi2tobb*(rphi2_QdQdE + rphi2_QdQdO/(
beta(Mb,
m2_2)*
beta(Mb,
m2_2)))*Gammaphi2totSM/Gammaphi2tot;
6208 Br_phi2totautau=BrSM_phi2totautau*(rphi2_QlQlE + rphi2_QlQlO/(
beta(Mtau,
m2_2)*
beta(Mtau,
m2_2)))*Gammaphi2totSM/Gammaphi2tot;
6210 Br_phi2toWW=BrSM_phi2toWW*rphi2_VV*Gammaphi2totSM/Gammaphi2tot;
6211 Br_phi2toZZ=BrSM_phi2toZZ*rphi2_VV*Gammaphi2totSM/Gammaphi2tot;
6212 Br_phi2togaga=Gamma_phi2gaga/Gammaphi2tot;
6213 Br_phi2toZga=Gamma_phi2Zga/Gammaphi2tot;
6214 Br_phi2tophi1phi1=Gammaphi2_phi1phi1/Gammaphi2tot;
6215 Br_phi2tophi3phi3=Gammaphi2_phi3phi3/Gammaphi2tot;
6216 Br_phi2tophi1phi3=Gammaphi2_phi1phi3/Gammaphi2tot;
6217 Br_phi2toHpHm=Gammaphi2_HpHm/Gammaphi2tot;
6218 Br_phi2tophi1Z=Gammaphi2_phi1Z/Gammaphi2tot;
6219 Br_phi2tophi3Z=Gammaphi2_phi3Z/Gammaphi2tot;
6220 Br_phi2toHpW=Gammaphi2_HpW/Gammaphi2tot;
6237 double GF=1/(
sqrt(2.0)*vev*vev);
6242 if (!myGTHDM->getATHDMflag())
6244 throw std::runtime_error(
"Direct Searches are only aviable in the A2HDM.");
6289 double rphi3_QuQuE= yu3.
real()*yu3.
real();
6290 double rphi3_QuQuO= yu3.
imag()*yu3.
imag();
6291 double rphi3_QdQdE= yd3.
real()*yd3.
real();
6292 double rphi3_QdQdO= yd3.
imag()*yd3.
imag();
6293 double rphi3_QlQlE= yl3.real()*yl3.real();
6294 double rphi3_QlQlO= yl3.imag()*yl3.imag();
6314 double lambdaphi3HpHm =
lambdaipm(R31, R32, R33);
6326 double Gamma_phi3gaga=(GF*Ale*Ale*m3*m3*m3/(
sqrt(2.0)*128.0*M_PI*M_PI*M_PI))*((I_phi3E_F+I_phi3_W+I_phi3_Hp).abs2()
6327 + (I_phi3O_F).abs2());
6351 double Gamma_phi3Zga=HSTheta(m3-MZ)*GF*Ale*Ale*m3*m3*m3/(
sqrt(2.0)*64.0*M_PI*M_PI*M_PI)*(1.0-MZ*MZ/(m3*m3))*(1.0-MZ*MZ/(m3*m3))*(1.0-MZ*MZ/(m3*m3))*((A_phi3E_F+A_phi3_W+A_phi3_Hp).abs2()
6352 + A_phi3O_F.
abs2());
6357 double Gamma_phi3gg=(rphi3_ggE)*GF*Als*Als*m3*m3*m3/(
sqrt(2.0)*16.0*M_PI*M_PI*M_PI)*(9.0/4.0)*(I_HH_Ux/4.0+I_HH_Dx).abs2()
6358 +rphi3_ggO*GF*Als*Als*m3*m3*m3/(
sqrt(2.0)*16.0*M_PI*M_PI*M_PI)*(9.0/4.0)*(I_A_Ux/4.0+I_A_Dx).abs2();
6375 SigmaTotSM_phi3_8 = 1.0e-15;
6377 if (m3>=20. && m3 <=2000.) {
6380 SigmaSumphi3_8 = SigmaggF_phi3_8 + SigmaVBF_phi3_8 + SigmaVH_phi3_8 + SigmattF_phi3_8 + SigmabbF_phi3_8;
6397 SigmaSumphi3_13 = SigmaggF_phi3_13 + SigmaVBF_phi3_13 + SigmaVH_phi3_13 + SigmattF_phi3_13 + SigmabbF_phi3_13;
6413 double lambda123 =
lambdaijk(R11, R21, R13, R12, R22, R32, R13, R23, R33, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R12, R22, R32, R11, R22, R32, R12, R22, R33, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R12, R22, R32, R13, R23, R33, R11, R21, R33, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R11, R21, R31, R13, R23, R33, R12, R22, R32, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R13, R23,R33, R11, R21, R31, R12, R22, R32,lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R13, R23,R33, R12, R22, R32, R11, R21, R31,lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) ;
6414 double lambda223 =
lambdaijk(R11, R22, R32, R12, R22, R32, R13, R23, R33, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R12, R22, R32, R13, R23, R33, R12, R22, R32, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R13, R23, R33, R12, R22, R32, R12, R22, R32, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) ;
6415 double lambda113 =
lambdaijk(R11, R21, R31, R11, R21, R31, R13, R23, R33, lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R11, R21, R31, R13, R23, R33, R11, R21, R31,lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) +
lambdaijk(R13, R23, R33, R11, R21, R31, R11, R21, R31,lambda1H, lambda3H, lambda4H, Relambda5H, Imlambda5H, Relambda6H, Imlambda6H, Relambda7H, Imlambda7H) ;
6419 double Gammaphi3_phi2phi2=HSTheta(m3 - 2.0*m2)*
sqrt(std::fabs(1.0 - (4.0*
m2_2)/
m3_2))*lambda223*lambda223/(32.0*m3*M_PI);
6421 double Gammaphi3_phi1phi1=HSTheta(m3 - 2.0*m1)*
sqrt(std::fabs(1.0 - (4.0*
m1_2)/
m3_2))*lambda113*lambda113/(32.0*m3*M_PI);
6423 double Gammaphi3_HpHm=HSTheta(m3 - 2.0*
sqrt(mHp2))*
sqrt(std::fabs(1.0 - (4.0*mHp2)/
m3_2))*lambdaphi3HpHm*lambdaphi3HpHm/(32.0*m3*M_PI);
6425 double Gammaphi3_phi1Z=HSTheta(m3-(m1+MZ))*
pow(
KaellenFunction(
m3_2,MZ*MZ,
m1_2),3)*(R33*R12 + R32*R13)*(R33*R12 + R32*R13)/(2.0*M_PI*vev*vev);
6427 double Gammaphi3_phi2Z=HSTheta(m3-(m2+MZ))*
pow(
KaellenFunction(
m3_2,MZ*MZ,
m2_2),3)*(R33*R22 + R32*R23)*(R33*R22 + R32*R23)/(2.0*M_PI*vev*vev);
6431 Gammaphi3tot = 1.e-10;
6433 Gammaphi3tot= Gammaphi3tot + (BrSM_phi3tott*(rphi3_QuQuE + rphi3_QuQuO/(
beta_mt_sq(Mt,
m3_2)))
6436 +BrSM_phi3totautau*(rphi3_QlQlE + rphi3_QlQlO/(
beta(Mtau,
m3_2)*
beta(Mtau,
m3_2)))
6437 +(BrSM_phi3toWW+BrSM_phi3toZZ)*rphi3_VV)*Gammaphi3totSM+Gamma_phi3gaga+Gamma_phi3Zga+Gamma_phi3gg
6438 + Gammaphi3_phi1phi1+Gammaphi3_phi2phi2+Gammaphi3_phi1phi2+Gammaphi3_HpHm+Gammaphi3_phi1Z+Gammaphi3_phi2Z
6441 Br_phi3tott=BrSM_phi3tott*(rphi3_QuQuE + rphi3_QuQuO/(
beta_mt_sq(Mt,
m3_2)))*Gammaphi3totSM/Gammaphi3tot;
6442 Br_phi3tobb=BrSM_phi3tobb*(rphi3_QdQdE + rphi3_QdQdO/(
beta(Mb,
m3_2)*
beta(Mb,
m3_2)))*Gammaphi3totSM/Gammaphi3tot;
6445 Br_phi3totautau=BrSM_phi3totautau*(rphi3_QlQlE + rphi3_QlQlO/(
beta(Mtau,
m3_2)*
beta(Mtau,
m3_2)))*Gammaphi3totSM/Gammaphi3tot;
6446 Br_phi3toWW=BrSM_phi3toWW*rphi3_VV*Gammaphi3totSM/Gammaphi3tot;
6447 Br_phi3toZZ=BrSM_phi3toZZ*rphi3_VV*Gammaphi3totSM/Gammaphi3tot;
6448 Br_phi3togaga=Gamma_phi3gaga/Gammaphi3tot;
6449 Br_phi3toZga=Gamma_phi3Zga/Gammaphi3tot;
6450 Br_phi3tophi1phi1=Gammaphi3_phi1phi1/Gammaphi3tot;
6451 Br_phi3tophi2phi2=Gammaphi3_phi2phi2/Gammaphi3tot;
6453 Br_phi3tophi1phi2=Gammaphi3_phi1phi2/Gammaphi3tot;
6454 Br_phi3toHpHm=Gammaphi3_HpHm/Gammaphi3tot;
6455 Br_phi3tophi1Z=Gammaphi3_phi1Z/Gammaphi3tot;
6456 Br_phi3tophi2Z=Gammaphi3_phi2Z/Gammaphi3tot;
6457 Br_phi3toHpW=Gammaphi3_HpW/Gammaphi3tot;
6478 double Mtau2=Mtau*Mtau;
6480 double Vtb=myGTHDM->getCKM().getV_tb().abs();
6487 if (!myGTHDM->getATHDMflag())
6489 throw std::runtime_error(
"Direct Searches are only aviable in the A2HDM.");
6499 double GammaHptaunu=HSTheta(mHp-Mtau)*(Mtau2*(mHp2-Mtau2)*(mHp2-Mtau2)*sl.abs2())/(8.0*mHp*mHp2*M_PI*vev*vev);
6503 double GammaHptb=HSTheta(mHp-Mt-Mb)*(Vtb*Vtb/(8.0*mHp*M_PI*vev*vev))*3.0*(-4.0*(su*sd).real()*Mb2*Mt2
6504 -sd.abs2()*Mb2*(Mb2-mHp2+Mt2)-su.abs2()*Mt2*(Mb2-mHp2+Mt2))
6505 *
sqrt((Mb2*Mb2+(mHp2-Mt2)*(mHp2-Mt2)-2.0*Mb2*(mHp2+Mt2))/(mHp2*mHp2));
6510 *MW2*MW2/mHp*(R21-R31)*(R21-R31)/(2.0*M_PI*vev*vev);
6514 *MW2*MW2/mHp*(R22-R32)*(R22-R32)/(2.0*M_PI*vev*vev);
6518 *MW2*MW2/mHp*(R23-R33)*(R23-R33)/(2.0*M_PI*vev*vev);
6520 GammaHptot = 1.e-10;
6522 GammaHptot= GammaHptot + GammaHptaunu + GammaHptb + GammaHpHlW + GammaHpphi2W + GammaHpphi3W;
6525 Br_Hptotaunu=GammaHptaunu/GammaHptot;
6526 Br_Hptotb=GammaHptb/GammaHptot;
6544 if (!myGTHDM->getATHDMflag())
6546 throw std::runtime_error(
"Direct Searches are only aviable in the A2HDM.");
6552 double Br_Ztoee=0.03363;
6553 double Br_Ztomumu=0.03366;
6554 double Br_Ztotautau=0.0337;
6555 double Br_Ztoinv=0.2;
6556 double Br_Wtoenu=0.1071;
6557 double Br_Wtomunu=0.1063;
6558 double Br_Wtotaunu=0.1138;
6560 THoEX_bb_phi2_tt_ATLAS13=0.0;
6561 THoEX_bb_phi3_tt_ATLAS13=0.0;
6562 THoEX_bb_phi2_bb_CMS8=0.0;
6563 THoEX_bb_phi3_bb_CMS8=0.0;
6564 THoEX_gg_phi2_bb_CMS8=0.0;
6565 THoEX_gg_phi3_bb_CMS8=0.0;
6566 THoEX_pp_phi2_bb_CMS13=0.0;
6567 THoEX_pp_phi3_bb_CMS13=0.0;
6568 THoEX_bb_phi2_bb_CMS13=0.0;
6569 THoEX_bb_phi3_bb_CMS13=0.0;
6570 THoEX_gg_phi2_tautau_ATLAS8=0.0;
6571 THoEX_gg_phi3_tautau_ATLAS8=0.0;
6572 THoEX_gg_phi2_tautau_CMS8=0.0;
6573 THoEX_gg_phi3_tautau_CMS8=0.0;
6574 THoEX_bb_phi2_tautau_ATLAS8=0.0;
6575 THoEX_bb_phi3_tautau_ATLAS8=0.0;
6576 THoEX_bb_phi2_tautau_CMS8=0.0;
6577 THoEX_bb_phi3_tautau_CMS8=0.0;
6578 THoEX_gg_phi2_tautau_ATLAS13=0.0;
6579 THoEX_gg_phi3_tautau_ATLAS13=0.0;
6580 THoEX_gg_phi2_tautau_CMS13=0.0;
6581 THoEX_gg_phi3_tautau_CMS13=0.0;
6582 THoEX_bb_phi2_tautau_ATLAS13=0.0;
6583 THoEX_bb_phi3_tautau_ATLAS13=0.0;
6584 THoEX_bb_phi2_tautau_CMS13=0.0;
6585 THoEX_bb_phi3_tautau_CMS13=0.0;
6586 THoEX_gg_phi2_gaga_ATLAS8=0.0;
6587 THoEX_gg_phi3_gaga_ATLAS8=0.0;
6588 THoEX_pp_phi2_gaga_ATLAS13=0.0;
6589 THoEX_pp_phi3_gaga_ATLAS13=0.0;
6590 THoEX_gg_phi2_gaga_CMS13=0.0;
6591 THoEX_gg_phi3_gaga_CMS13=0.0;
6592 THoEX_pp_phi2_Zga_llga_ATLAS8=0.0;
6593 THoEX_pp_phi3_Zga_llga_ATLAS8=0.0;
6594 THoEX_pp_phi2_Zga_llga_CMS8=0.0;
6595 THoEX_pp_phi3_Zga_llga_CMS8=0.0;
6596 THoEX_gg_phi2_Zga_llga_ATLAS13=0.0;
6597 THoEX_gg_phi3_Zga_llga_ATLAS13=0.0;
6598 THoEX_gg_phi2_Zga_qqga_ATLAS13=0.0;
6599 THoEX_gg_phi3_Zga_qqga_ATLAS13=0.0;
6600 THoEX_gg_phi2_Zga_CMS13=0.0;
6601 THoEX_gg_phi3_Zga_CMS13=0.0;
6602 THoEX_gg_phi2_ZZ_ATLAS8=0.0;
6603 THoEX_gg_phi3_ZZ_ATLAS8=0.0;
6604 THoEX_VV_phi2_ZZ_ATLAS8=0.0;
6605 THoEX_VV_phi3_ZZ_ATLAS8=0.0;
6606 THoEX_gg_phi2_ZZ_llllnunu_ATLAS13=0.0;
6607 THoEX_gg_phi3_ZZ_llllnunu_ATLAS13=0.0;
6608 THoEX_VV_phi2_ZZ_llllnunu_ATLAS13=0.0;
6609 THoEX_VV_phi3_ZZ_llllnunu_ATLAS13=0.0;
6610 THoEX_gg_phi2_ZZ_qqllnunu_ATLAS13=0.0;
6611 THoEX_gg_phi3_ZZ_qqllnunu_ATLAS13=0.0;
6612 THoEX_VV_phi2_ZZ_qqllnunu_ATLAS13=0.0;
6613 THoEX_VV_phi3_ZZ_qqllnunu_ATLAS13=0.0;
6614 THoEX_pp_phi2_ZZ_llqqnunull_CMS13=0.0;
6615 THoEX_pp_phi3_ZZ_llqqnunull_CMS13=0.0;
6616 THoEX_pp_phi2_ZZ_qqnunu_CMS13=0.0;
6617 THoEX_pp_phi3_ZZ_qqnunu_CMS13=0.0;
6618 THoEX_gg_phi2_WW_ATLAS8=0.0;
6619 THoEX_gg_phi3_WW_ATLAS8=0.0;
6620 THoEX_VV_phi2_WW_ATLAS8=0.0;
6621 THoEX_VV_phi3_WW_ATLAS8=0.0;
6622 THoEX_gg_phi2_WW_enumunu_ATLAS13=0.0;
6623 THoEX_gg_phi3_WW_enumunu_ATLAS13=0.0;
6624 THoEX_VV_phi2_WW_enumunu_ATLAS13=0.0;
6625 THoEX_VV_phi3_WW_enumunu_ATLAS13=0.0;
6626 THoEX_ggVV_phi2_WW_lnulnu_CMS13=0.0;
6627 THoEX_ggVV_phi3_WW_lnulnu_CMS13=0.0;
6628 THoEX_gg_phi2_WW_lnuqq_ATLAS13=0.0;
6629 THoEX_gg_phi3_WW_lnuqq_ATLAS13=0.0;
6630 THoEX_VV_phi2_WW_lnuqq_ATLAS13=0.0;
6631 THoEX_VV_phi3_WW_lnuqq_ATLAS13=0.0;
6632 THoEX_pp_phi2_WW_lnuqq_CMS13=0.0;
6633 THoEX_pp_phi3_WW_lnuqq_CMS13=0.0;
6634 THoEX_mu_pp_phi2_VV_CMS8=0.0;
6635 THoEX_mu_pp_phi3_VV_CMS8=0.0;
6636 THoEX_pp_phi2_VV_qqqq_ATLAS13=0.0;
6637 THoEX_pp_phi3_VV_qqqq_ATLAS13=0.0;
6638 THoEX_gg_phi2_phi1phi1_ATLAS8=0.0;
6639 THoEX_gg_phi3_phi1phi1_ATLAS8=0.0;
6640 THoEX_pp_phi2_phi1phi1_bbbb_CMS8=0.0;
6641 THoEX_pp_phi3_phi1phi1_bbbb_CMS8=0.0;
6642 THoEX_pp_phi2_phi1phi1_bbgaga_CMS8=0.0;
6643 THoEX_pp_phi3_phi1phi1_bbgaga_CMS8=0.0;
6644 THoEX_gg_phi2_phi1phi1_bbtautau_CMS8=0.0;
6645 THoEX_gg_phi3_phi1phi1_bbtautau_CMS8=0.0;
6646 THoEX_pp_phi2_phi1phi1_bbtautau_CMS8=0.0;
6647 THoEX_pp_phi3_phi1phi1_bbtautau_CMS8=0.0;
6648 THoEX_pp_phi2_phi1phi1_bbbb_ATLAS13=0.0;
6649 THoEX_pp_phi3_phi1phi1_bbbb_ATLAS13=0.0;
6650 THoEX_pp_phi2_phi1phi1_bbbb_1_CMS13=0.0;
6651 THoEX_pp_phi3_phi1phi1_bbbb_1_CMS13=0.0;
6652 THoEX_pp_phi2_phi1phi1_bbbb_2_CMS13=0.0;
6653 THoEX_pp_phi3_phi1phi1_bbbb_2_CMS13=0.0;
6654 THoEX_pp_phi2_phi1phi1_bbgaga_ATLAS13=0.0;
6655 THoEX_pp_phi3_phi1phi1_bbgaga_ATLAS13=0.0;
6656 THoEX_pp_phi2_phi1phi1_bbgaga_CMS13=0.0;
6657 THoEX_pp_phi3_phi1phi1_bbgaga_CMS13=0.0;
6658 THoEX_pp_phi2_phi1phi1_bbtautau_ATLAS13=0.0;
6659 THoEX_pp_phi3_phi1phi1_bbtautau_ATLAS13=0.0;
6660 THoEX_pp_phi2_phi1phi1_bbtautau_1_CMS13=0.0;
6661 THoEX_pp_phi3_phi1phi1_bbtautau_1_CMS13=0.0;
6662 THoEX_pp_phi2_phi1phi1_bbtautau_2_CMS13=0.0;
6663 THoEX_pp_phi3_phi1phi1_bbtautau_2_CMS13=0.0;
6664 THoEX_pp_phi2_phi1phi1_bbVV_CMS13=0.0;
6665 THoEX_pp_phi3_phi1phi1_bbVV_CMS13=0.0;
6666 THoEX_pp_phi2_phi1phi1_bbWW_ATLAS13=0.0;
6667 THoEX_pp_phi3_phi1phi1_bbWW_ATLAS13=0.0;
6668 THoEX_gg_phi2_phi1phi1_gagaWW_ATLAS13=0.0;
6669 THoEX_gg_phi3_phi1phi1_gagaWW_ATLAS13=0.0;
6670 THoEX_gg_phi2_phi1Z_bbZ_ATLAS8=0.0;
6671 THoEX_gg_phi3_phi1Z_bbZ_ATLAS8=0.0;
6672 THoEX_gg_phi2_phi1Z_bbll_CMS8=0.0;
6673 THoEX_gg_phi3_phi1Z_bbll_CMS8=0.0;
6674 THoEX_gg_phi2_phi1Z_tautauZ_ATLAS8=0.0;
6675 THoEX_gg_phi3_phi1Z_tautauZ_ATLAS8=0.0;
6676 THoEX_gg_phi2_phi1Z_tautaull_CMS8=0.0;
6677 THoEX_gg_phi3_phi1Z_tautaull_CMS8=0.0;
6678 THoEX_gg_phi2_phi1Z_bbZ_ATLAS13=0.0;
6679 THoEX_gg_phi3_phi1Z_bbZ_ATLAS13=0.0;
6680 THoEX_gg_phi2_phi1Z_bbZ_1_CMS13=0.0;
6681 THoEX_gg_phi3_phi1Z_bbZ_1_CMS13=0.0;
6682 THoEX_gg_phi2_phi1Z_bbZ_2_CMS13=0.0;
6683 THoEX_gg_phi3_phi1Z_bbZ_2_CMS13=0.0;
6684 THoEX_bb_phi2_phi1Z_bbZ_ATLAS13=0.0;
6685 THoEX_bb_phi3_phi1Z_bbZ_ATLAS13=0.0;
6686 THoEX_bb_phi2_phi1Z_bbZ_1_CMS13=0.0;
6687 THoEX_bb_phi3_phi1Z_bbZ_1_CMS13=0.0;
6688 THoEX_bb_phi2_phi1Z_bbZ_2_CMS13=0.0;
6689 THoEX_bb_phi3_phi1Z_bbZ_2_CMS13=0.0;
6691 THoEX_pp_phi3_phi2Z_bbll_1_CMS8=0.0;
6692 THoEX_pp_phi2_phi3Z_bbll_1_CMS8=0.0;
6693 THoEX_pp_phi3_phi2Z_bbll_2_CMS8=0.0;
6694 THoEX_pp_phi2_phi3Z_bbll_2_CMS8=0.0;
6695 THoEX_pp_phi3_phi2Z_tautaull_1_CMS8=0.0;
6696 THoEX_pp_phi2_phi3Z_tautaull_1_CMS8=0.0;
6697 THoEX_pp_phi3_phi2Z_tautaull_2_CMS8=0.0;
6698 THoEX_pp_phi2_phi3Z_tautaull_2_CMS8=0.0;
6699 THoEX_gg_phi3_phi2Z_bbZ_ATLAS13=0.0;
6700 THoEX_gg_phi2_phi3Z_bbZ_ATLAS13=0.0;
6701 THoEX_bb_phi3_phi2Z_bbZ_ATLAS13=0.0;
6702 THoEX_bb_phi2_phi3Z_bbZ_ATLAS13=0.0;
6704 THoEX_pp_Hpm_taunu_ATLAS8=0.0;
6705 THoEX_pp_Hp_taunu_CMS8=0.0;
6706 THoEX_pp_Hpm_taunu_ATLAS13=0.0;
6707 THoEX_pp_Hpm_taunu_CMS13=0.0;
6708 THoEX_pp_Hpm_tb_ATLAS8=0.0;
6709 THoEX_pp_Hp_tb_CMS8=0.0;
6710 THoEX_pp_Hpm_tb_ATLAS13=0.0;
6714 tt_phi2_tt_TH13=SigmattF_phi2_13*Br_phi2tott;
6715 tt_phi3_tt_TH13=SigmattF_phi3_13*Br_phi3tott;
6716 bb_phi2_tt_TH13=SigmabbF_phi2_13*Br_phi2tott;
6717 bb_phi3_tt_TH13=SigmabbF_phi3_13*Br_phi3tott;
6719 bb_phi2_bb_TH8=SigmabbF_phi2_8*Br_phi2tobb;
6720 bb_phi3_bb_TH8=SigmabbF_phi3_8*Br_phi3tobb;
6721 gg_phi2_bb_TH8=SigmaggF_phi2_8*Br_phi2tobb;
6722 gg_phi3_bb_TH8=SigmaggF_phi3_8*Br_phi3tobb;
6723 pp_phi2_bb_TH13=SigmaSumphi2_13*Br_phi2tobb;
6724 pp_phi3_bb_TH13=SigmaSumphi3_13*Br_phi3tobb;
6725 bb_phi2_bb_TH13=SigmabbF_phi2_13*Br_phi2tobb;
6726 bb_phi3_bb_TH13=SigmabbF_phi3_13*Br_phi3tobb;
6728 gg_phi2_tautau_TH8=SigmaggF_phi2_8*Br_phi2totautau;
6729 gg_phi3_tautau_TH8=SigmaggF_phi3_8*Br_phi3totautau;
6730 bb_phi2_tautau_TH8=SigmabbF_phi2_8*Br_phi2totautau;
6731 bb_phi3_tautau_TH8=SigmabbF_phi3_8*Br_phi3totautau;
6732 gg_phi2_tautau_TH13=SigmaggF_phi2_13*Br_phi2totautau;
6733 gg_phi3_tautau_TH13=SigmaggF_phi3_13*Br_phi3totautau;
6734 bb_phi2_tautau_TH13=SigmabbF_phi2_13*Br_phi2totautau;
6735 bb_phi3_tautau_TH13=SigmabbF_phi3_13*Br_phi3totautau;
6737 gg_phi2_gaga_TH8=SigmaggF_phi2_8*Br_phi2togaga;
6738 gg_phi3_gaga_TH8=SigmaggF_phi3_8*Br_phi3togaga;
6739 pp_phi2_gaga_TH13=SigmaSumphi2_13*Br_phi2togaga;
6740 pp_phi3_gaga_TH13=SigmaSumphi3_13*Br_phi3togaga;
6741 gg_phi2_gaga_TH13=SigmaggF_phi2_13*Br_phi2togaga;
6742 gg_phi3_gaga_TH13=SigmaggF_phi3_13*Br_phi3togaga;
6744 pp_phi2_Zga_llga_TH8=SigmaSumphi2_8*Br_phi2toZga*(Br_Ztoee+Br_Ztomumu);
6745 pp_phi3_Zga_llga_TH8=SigmaSumphi3_8*Br_phi3toZga*(Br_Ztoee+Br_Ztomumu);
6746 gg_phi2_Zga_TH13=SigmaggF_phi2_13*Br_phi2toZga;
6747 gg_phi3_Zga_TH13=SigmaggF_phi3_13*Br_phi3toZga;
6749 gg_phi2_ZZ_TH8=SigmaggF_phi2_8*Br_phi2toZZ;
6750 gg_phi3_ZZ_TH8=SigmaggF_phi3_8*Br_phi3toZZ;
6751 VV_phi2_ZZ_TH8=SigmaVBF_phi2_8*Br_phi2toZZ;
6752 VV_phi3_ZZ_TH8=SigmaVBF_phi3_8*Br_phi3toZZ;
6753 gg_phi2_ZZ_TH13=SigmaggF_phi2_13*Br_phi2toZZ;
6754 gg_phi3_ZZ_TH13=SigmaggF_phi3_13*Br_phi3toZZ;
6755 VV_phi2_ZZ_TH13=SigmaVBF_phi2_13*Br_phi2toZZ;
6756 VV_phi3_ZZ_TH13=SigmaVBF_phi3_13*Br_phi3toZZ;
6757 pp_phi2_ZZ_TH13=SigmaSumphi2_13*Br_phi2toZZ;
6758 pp_phi3_ZZ_TH13=SigmaSumphi3_13*Br_phi3toZZ;
6760 gg_phi2_WW_TH8=SigmaggF_phi2_8*Br_phi2toWW;
6761 gg_phi3_WW_TH8=SigmaggF_phi3_8*Br_phi3toWW;
6762 VV_phi2_WW_TH8=SigmaVBF_phi2_8*Br_phi2toWW;
6763 VV_phi3_WW_TH8=SigmaVBF_phi3_8*Br_phi3toWW;
6764 gg_phi2_WW_TH13=SigmaggF_phi2_13*Br_phi2toWW;
6765 gg_phi3_WW_TH13=SigmaggF_phi3_13*Br_phi3toWW;
6766 VV_phi2_WW_TH13=SigmaVBF_phi2_13*Br_phi2toWW;
6767 VV_phi3_WW_TH13=SigmaVBF_phi3_13*Br_phi3toWW;
6768 ggVV_phi2_WW_lnulnu_TH13=(SigmaggF_phi2_13+SigmaVBF_phi2_13)*Br_phi2toWW*(Br_Wtoenu+Br_Wtomunu)*(Br_Wtoenu+Br_Wtomunu);
6769 ggVV_phi3_WW_lnulnu_TH13=(SigmaggF_phi3_13+SigmaVBF_phi3_13)*Br_phi3toWW*(Br_Wtoenu+Br_Wtomunu)*(Br_Wtoenu+Br_Wtomunu);
6770 pp_phi2_WW_TH13=SigmaSumphi2_13*Br_phi2toWW;
6771 pp_phi3_WW_TH13=SigmaSumphi3_13*Br_phi3toWW;
6773 mu_pp_phi2_VV_TH8=SigmaSumphi2_8/SigmaTotSM_phi2_8*rphi2_VV*Gammaphi2totSM/Gammaphi2tot;
6774 mu_pp_phi3_VV_TH8=SigmaSumphi3_8/SigmaTotSM_phi3_8*rphi3_VV*Gammaphi3totSM/Gammaphi3tot;
6775 pp_phi2_VV_TH13=SigmaSumphi2_13*(Br_phi2toZZ+Br_phi2toWW);
6776 pp_phi3_VV_TH13=SigmaSumphi3_13*(Br_phi3toZZ+Br_phi3toWW);
6778 gg_phi2_phi1phi1_TH8=SigmaggF_phi2_8*Br_phi2tophi1phi1;
6779 gg_phi3_phi1phi1_TH8=SigmaggF_phi3_8*Br_phi3tophi1phi1;
6786 pp_phi2_phi1phi1_TH8=SigmaSumphi2_8*Br_phi2tophi1phi1;
6787 pp_phi3_phi1phi1_TH8=SigmaSumphi3_8*Br_phi3tophi1phi1;
6790 pp_phi2_phi1phi1_TH13=SigmaSumphi2_13*Br_phi2tophi1phi1;
6791 pp_phi3_phi1phi1_TH13=SigmaSumphi3_13*Br_phi3tophi1phi1;
6796 pp_phi2_phi1phi1_bbVV_TH13=SigmaSumphi2_13*Br_phi2tophi1phi1*2.0*
GTHDM_BR_h_bb*(GTHDM_BR_h_WW*
pow(Br_Wtoenu+Br_Wtomunu+Br_Wtotaunu*0.352,2)
6797 +GTHDM_BR_h_ZZ*2.0*Br_Ztoinv*(Br_Ztoee+Br_Ztomumu+Br_Ztotautau*0.124));
6798 pp_phi3_phi1phi1_bbVV_TH13=SigmaSumphi3_13*Br_phi3tophi1phi1*2.0*
GTHDM_BR_h_bb*(GTHDM_BR_h_WW*
pow(Br_Wtoenu+Br_Wtomunu+Br_Wtotaunu*0.352,2)
6799 +GTHDM_BR_h_ZZ*2.0*Br_Ztoinv*(Br_Ztoee+Br_Ztomumu+Br_Ztotautau*0.124));
6800 pp_phi2_phi1phi1_bbWW_TH13=SigmaSumphi2_13*Br_phi2tophi1phi1*2.0*5.77e-1*2.15e-1;
6801 pp_phi3_phi1phi1_bbWW_TH13=SigmaSumphi3_13*Br_phi3tophi1phi1*2.0*5.77e-1*2.15e-1;
6802 gg_phi2_phi1phi1_gagaWW_TH13=SigmaggF_phi2_13*Br_phi2tophi1phi1*
GTHDM_BR_h_gaga*GTHDM_BR_h_WW*2.0;
6803 gg_phi3_phi1phi1_gagaWW_TH13=SigmaggF_phi3_13*Br_phi3tophi1phi1*
GTHDM_BR_h_gaga*GTHDM_BR_h_WW*2.0;
6805 gg_phi2_phi1Z_bbZ_TH8=SigmaggF_phi2_8*Br_phi2tophi1Z*
GTHDM_BR_h_bb;
6806 gg_phi3_phi1Z_bbZ_TH8=SigmaggF_phi3_8*Br_phi3tophi1Z*
GTHDM_BR_h_bb;
6807 gg_phi2_phi1Z_bbll_TH8=SigmaggF_phi2_8*Br_phi2tophi1Z*
GTHDM_BR_h_bb*(Br_Ztoee+Br_Ztomumu);
6808 gg_phi3_phi1Z_bbll_TH8=SigmaggF_phi3_8*Br_phi3tophi1Z*
GTHDM_BR_h_bb*(Br_Ztoee+Br_Ztomumu);
6811 gg_phi2_phi1Z_tautaull_TH8=SigmaggF_phi2_8*Br_phi2tophi1Z*
GTHDM_BR_h_tautau*(Br_Ztoee+Br_Ztomumu);
6812 gg_phi3_phi1Z_tautaull_TH8=SigmaggF_phi3_8*Br_phi3tophi1Z*
GTHDM_BR_h_tautau*(Br_Ztoee+Br_Ztomumu);
6813 gg_phi2_phi1Z_bbZ_TH13=SigmaggF_phi2_13*Br_phi2tophi1Z*
GTHDM_BR_h_bb;
6814 gg_phi3_phi1Z_bbZ_TH13=SigmaggF_phi3_13*Br_phi3tophi1Z*
GTHDM_BR_h_bb;
6815 bb_phi2_phi1Z_bbZ_TH13=SigmabbF_phi2_13*Br_phi2tophi1Z*
GTHDM_BR_h_bb;
6816 bb_phi3_phi1Z_bbZ_TH13=SigmabbF_phi3_13*Br_phi3tophi1Z*
GTHDM_BR_h_bb;
6818 pp_phi3_phi2Z_bbll_TH8=SigmaSumphi3_8*Br_phi3tophi2Z*Br_phi2tobb*(Br_Ztoee+Br_Ztomumu);
6819 pp_phi2_phi3Z_bbll_TH8=SigmaSumphi2_8*Br_phi2tophi3Z*Br_phi3tobb*(Br_Ztoee+Br_Ztomumu);
6820 pp_phi3_phi2Z_tautaull_TH8=SigmaSumphi3_8*Br_phi3tophi2Z*Br_phi2totautau*(Br_Ztoee+Br_Ztomumu);
6821 pp_phi2_phi3Z_tautaull_TH8=SigmaSumphi2_8*Br_phi2tophi3Z*Br_phi3totautau*(Br_Ztoee+Br_Ztomumu);
6822 gg_phi3_phi2Z_bbZ_TH13=SigmaggF_phi3_13*Br_phi3tophi2Z*Br_phi2tobb;
6823 gg_phi2_phi3Z_bbZ_TH13=SigmaggF_phi2_13*Br_phi2tophi3Z*Br_phi3tobb;
6824 bb_phi3_phi2Z_bbZ_TH13=SigmabbF_phi3_13*Br_phi3tophi2Z*Br_phi2tobb;
6825 bb_phi2_phi3Z_bbZ_TH13=SigmabbF_phi2_13*Br_phi2tophi3Z*Br_phi3tobb;
6827 pp_Hpm_taunu_TH8=2.0*SigmaHp8*Br_Hptotaunu;
6828 pp_Hp_taunu_TH8=SigmaHp8*Br_Hptotaunu;
6829 pp_Hpm_taunu_TH13=2.0*SigmaHpm13*Br_Hptotaunu;
6830 pp_Hpm_tb_TH8=2.0*SigmaHp8*Br_Hptotb;
6831 pp_Hp_tb_TH8=SigmaHp8*Br_Hptotb;
6832 pp_Hpm_tb_TH13=2.0*SigmaHpm13*Br_Hptotb;
6993 void GeneralTHDMcache::runGeneralTHDMparameters()
6996 std::string RGEorder=myGTHDM->getRGEorderflag();
7000 if( RGEorder ==
"LO" ) flag=0;
7004 throw std::runtime_error(
"RGEorder can be only \"LO\" at the moment");
7007 if(myGTHDM->getATHDMflag() && myGTHDM->getCPconservationflag())
7009 double g1_at_MZ=
sqrt(4.0*M_PI*Ale/cW2);
7010 double g2_at_MZ=
sqrt(4.0*M_PI*Ale/(1-cW2));
7011 double g3_at_MZ=
sqrt(4.0*M_PI*Als);
7014 double etaU1_at_MZ=0.;
7015 double etaU2_at_MZ=0.;
7016 double etaD1_at_MZ=0.;
7017 double etaD2_at_MZ=0.;
7018 double etaL1_at_MZ=0.;
7019 double etaL2_at_MZ=0.;
7020 double m11sq_at_MZ=0.;
7021 double m12sq_at_MZ=0.;
7022 double m22sq_at_MZ=0.;
7027 double lambda5_at_MZ=Relambda5;
7028 double lambda6_at_MZ=Relambda6;
7029 double lambda7_at_MZ=Relambda7;
7030 double RpepsGTHDM=myGTHDM->getRpepsGTHDM();
7031 double NLOuniscale=myGTHDM->getNLOuniscaleGTHDM();
7033 if(fabs(Q_GTHDM-
log10(MZ))<0.005)
7042 etaU1_at_Q = etaU1_at_MZ;
7043 etaU2_at_Q = etaU2_at_MZ;
7044 etaD1_at_Q = etaD1_at_MZ;
7045 etaD2_at_Q = etaD2_at_MZ;
7046 etaL1_at_Q = etaL1_at_MZ;
7047 etaL2_at_Q = etaL2_at_MZ;
7048 m11sq_at_Q = m11sq_at_MZ;
7049 m12sq_at_Q = m12sq_at_MZ;
7050 m22sq_at_Q = m22sq_at_MZ;
7051 lambda1_at_Q = lambda1_at_MZ;
7052 lambda2_at_Q = lambda2_at_MZ;
7053 lambda3_at_Q = lambda3_at_MZ;
7054 lambda4_at_Q = lambda4_at_MZ;
7055 Relambda5_at_Q = lambda5_at_MZ;
7056 Relambda6_at_Q = lambda6_at_MZ;
7057 Relambda7_at_Q = lambda7_at_MZ;
7061 double InitVals[21];
7062 InitVals[0]=g1_at_MZ;
7063 InitVals[1]=g2_at_MZ;
7064 InitVals[2]=g3_at_MZ;
7065 InitVals[3]=v1_at_MZ;
7066 InitVals[4]=v2_at_MZ;
7067 InitVals[5]=etaU1_at_MZ;
7068 InitVals[6]=etaU2_at_MZ;
7069 InitVals[7]=etaD1_at_MZ;
7070 InitVals[8]=etaD2_at_MZ;
7071 InitVals[9]=etaL1_at_MZ;
7072 InitVals[10]=etaL2_at_MZ;
7073 InitVals[11]=m11sq_at_MZ;
7074 InitVals[12]=m12sq_at_MZ;
7075 InitVals[13]=m22sq_at_MZ;
7076 InitVals[14]=lambda1_at_MZ;
7077 InitVals[15]=lambda2_at_MZ;
7078 InitVals[16]=lambda3_at_MZ;
7079 InitVals[17]=lambda4_at_MZ;
7080 InitVals[18]=lambda5_at_MZ;
7081 InitVals[19]=lambda6_at_MZ;
7082 InitVals[20]=lambda7_at_MZ;
7084 Q_cutoff=myRunnerGTHDM->RGEGeneralTHDMRunner(InitVals, 21,
log10(MZ), Q_GTHDM, flag, RpepsGTHDM, NLOuniscale);
7086 g1_at_Q = InitVals[0];
7087 g2_at_Q = InitVals[1];
7088 g3_at_Q = InitVals[2];
7089 v1_at_Q = InitVals[3];
7090 v2_at_Q = InitVals[4];
7091 etaU1_at_Q = InitVals[5];
7092 etaU2_at_Q = InitVals[6];
7093 etaD1_at_Q = InitVals[7];
7094 etaD2_at_Q = InitVals[8];
7095 etaL1_at_Q = InitVals[9];
7096 etaL2_at_Q = InitVals[10];
7097 m11sq_at_Q = InitVals[11];
7098 m12sq_at_Q = InitVals[12];
7099 m22sq_at_Q = InitVals[13];
7100 lambda1_at_Q = InitVals[14];
7101 lambda2_at_Q = InitVals[15];
7102 lambda3_at_Q = InitVals[16];
7103 lambda4_at_Q = InitVals[17];
7104 Relambda5_at_Q = InitVals[18];
7105 Relambda6_at_Q = InitVals[19];
7106 Relambda7_at_Q = InitVals[20];
7110 throw std::runtime_error(
"RGE are only defined for the CP conserving ATHDM at the moment");
7116 mHl=myGTHDM->getMHl();
7118 mH2sq=myGTHDM->getmH2sq();
7119 mH3sq=myGTHDM->getmH3sq();
7120 mHp2=myGTHDM->getmHp2();
7122 tanb=myGTHDM->gettanb();
7123 cosb=myGTHDM->getcosb();
7124 sinb=myGTHDM->getsinb();
7125 cosa1=myGTHDM->getcosalpha1();
7126 sina1=myGTHDM->getsinalpha1();
7127 cosa2=myGTHDM->getcosalpha2();
7128 sina2=myGTHDM->getsinalpha2();
7129 cosa3=myGTHDM->getcosalpha3();
7130 sina3=myGTHDM->getsinalpha3();
7131 Relambda5=myGTHDM->getRelambda5();
7132 Imlambda5=myGTHDM->getImlambda5();
7133 Relambda6=myGTHDM->getRelambda6();
7134 Relambda7=myGTHDM->getRelambda7();
7139 R21_GTHDM = cosa1*sina2*sina3 - sina1*cosa3;
7140 R22_GTHDM = sina1*sina2*sina3 + cosa1*cosa3;
7142 R31_GTHDM = cosa2*sina2*cosa3 + sina1*sina3;
7143 R32_GTHDM = sina1*sina2*cosa3 - cosa1*sina3;
7149 R21 = cosa1*sina2*sina3 - sina1*cosa3;
7150 R22 = sina1*sina2*sina3 + cosa1*cosa3;
7152 R31 = cosa2*sina2*cosa3 + sina1*sina3;
7153 R32 = sina1*sina2*cosa3 - cosa1*sina3;
7159 M11_2 = -0.5*(mH1sq*cosa1*cosa1*cosa2*cosa2
7160 +mH2sq*sina1*sina1*cosa2*cosa2 + mH3sq*sina2*sina2);
7161 M12_2 = 0.5*cosa2*((mH1sq-mH2sq)*cosa1*cosa3*sina1
7162 +(-mH3sq+mH1sq*cosa1*cosa1+mH2sq*sina1*sina1)*sina2*sina3);
7163 M13_2 = 0.5*cosa2*(cosa3*(-mH3sq+mH1sq*cosa1*cosa1+mH2sq*sina1*sina1)*sina2
7164 +(mH2sq-mH1sq)*cosa1*sina1*sina3);
7165 M22_2 = -0.5*(mH3sq*cosa2*cosa2*sina3*sina3
7166 +mH1sq*(cosa3*sina1+cosa1*sina2*sina3)*(cosa3*sina1+cosa1*sina2*sina3)
7167 +mH2sq*(cosa1*cosa3-sina1*sina2*sina3)*(cosa1*cosa3-sina1*sina2*sina3));
7168 M23_2 = 0.5*((mH2sq-mH1sq)*cosa1*(cosa3*cosa3-sina3*sina3)*sina1*sina2
7169 +cosa1*cosa1*cosa3*(mH2sq-mH1sq*sina2*sina2)*sina3
7170 -cosa3*sina3*(mH3sq*cosa2*cosa2+sina1*sina1*(-mH1sq+mH2sq*sina2*sina2)));
7171 M33_2 = -0.5*(mH3sq*cosa2*cosa2*cosa3*cosa3
7172 +mH2sq*(cosa3*sina1*sina2+cosa1*sina3)*(cosa3*sina1*sina2+cosa1*sina3)
7173 +mH1sq*(cosa1*cosa3*sina2-sina1*sina3)*(cosa1*cosa3*sina2-sina1*sina3));
7176 m11sq = M11_2 - M33_2 - M12_2*tanb + 0.5*Relambda5*vev*vev
7177 + (M33_2-0.5*Relambda5*vev*vev)*(cosb*cosb-sinb*sinb)
7178 + 0.5*vev*vev*((Relambda6-Relambda7)*sinb*cosb+Relambda7*tanb);
7180 m22sq = M11_2 - M33_2 + M12_2/tanb + 0.5*Relambda5*vev*vev
7181 - (M33_2-0.5*Relambda5*vev*vev)*(cosb*cosb-sinb*sinb)
7182 + 0.25*vev*vev*(Relambda6+Relambda7+(Relambda6-Relambda7)*(cosb*cosb-sinb*sinb))/tanb;
7184 Rem12sq = 0.25*vev*vev*(Relambda6+Relambda7+(Relambda6-Relambda7)*(cosb*cosb-sinb*sinb))
7185 - (2.0*M33_2-Relambda5*vev*vev)*sinb*cosb;
7189 lambda1 = (-2.0*(M11_2-M22_2+M33_2) + Relambda5*vev*vev
7190 - (2.0*M22_2-2.0*M33_2+Relambda5*vev*vev)/(cosb*cosb)
7191 + (4.0*M12_2-2.0*Relambda6*vev*vev)*tanb)/(vev*vev);
7193 lambda2 = (-2.0*(M11_2-M22_2+M33_2) + Relambda5*vev*vev
7194 - (2.0*M22_2-2.0*M33_2+Relambda5*vev*vev)/(sinb*sinb)
7195 - (4.0*M12_2+2.0*Relambda7*vev*vev)/tanb)/(vev*vev);
7197 lambda3 = -(2.0*(M11_2-M22_2-M33_2-mHp2) + Relambda5*vev*vev
7198 + (2.0*M12_2+Relambda6*vev*vev)/tanb
7199 - (2.0*M12_2-Relambda7*vev*vev)*tanb)/(vev*vev);
7201 lambda4 = Relambda5 - (2.0*mHp2+4.0*M33_2)/(vev*vev);
7203 Imlambda6 = (2.0*M13_2-(2.0*M23_2+0.5*Imlambda5*vev*vev)*tanb)/(vev*vev);
7205 Imlambda7 = 2.0*M13_2/(vev*vev) + (-0.5*Imlambda5+(2.0*M23_2)/(vev*vev))/tanb;
7210 m22sqH = M11_2-2.0*M33_2+Relambda5*vev*vev
7211 +(M12_2+0.5*(Relambda6*vev*vev))/tanb
7212 -(M12_2-0.5*(Relambda7*vev*vev))*tanb;
7215 lambda1H = -2.0*M11_2/(vev*vev);
7216 lambda2H = -2.0*((2.0*M12_2+Relambda7*vev*vev)/tanb
7217 +M11_2-4.0*M22_2+4.0*M33_2-2.0*Relambda5*vev*vev
7218 +(M22_2-M33_2+0.5*Relambda5*vev*vev)/(sinb*sinb*cosb*cosb)
7219 -2.0*M12_2*tanb+Relambda6*vev*vev*tanb)/(vev*vev);
7220 lambda3H = -((2.0*(M11_2-2.0*M33_2-mHp2+Relambda5*vev*vev)
7221 +(2.0*M12_2+Relambda6*vev*vev)/tanb
7222 -(2.0*M12_2-Relambda7*vev*vev)*tanb)/(vev*vev));
7223 lambda4H = -2.0*(M22_2+M33_2+mHp2)/(vev*vev);
7224 Relambda5H = -2.0*(M22_2-M33_2)/(vev*vev);
7225 Imlambda5H = 4.0*M23_2/(vev*vev);
7226 Relambda6H = -2.0*M12_2/(vev*vev);
7227 Imlambda6H = 2.0*M13_2/(vev*vev);
7228 Relambda7H = (-2.0*M12_2+(Relambda6-Relambda7)*vev*vev
7229 +(2.0*M22_2-2.0*M33_2+Relambda5*vev*vev)*(tanb-1.0/tanb))/(vev*vev);
7230 Imlambda7H = 2.0*(M13_2-M23_2*(tanb-1.0/tanb))/(vev*vev)-0.5*Imlambda5/(sinb*cosb);
7233 M2 = Rem12sq/(sinb*cosb);
7253 if(mH1sq<mH3sq && mH3sq<mH2sq)
7260 else if(mH3sq<mH2sq && mH2sq<mH1sq)
7267 else if(mH2sq<mH1sq && mH1sq<mH3sq)
7274 else if(mH2sq<mH3sq && mH3sq<mH1sq)
7281 else if(mH3sq<mH1sq && mH1sq<mH2sq)
7289 else if(mH1sq<mH2sq && mH2sq<mH3sq)
7326 Mu_GTHDM.assign(0,0, myGTHDM->getQuarks(
QCD::UP).getMass());
7327 Mu_GTHDM.assign(1,1, myGTHDM->getQuarks(
QCD::CHARM).getMass());
7328 Mu_GTHDM.assign(2,2, myGTHDM->getQuarks(
QCD::TOP).getMass());
7330 Md_GTHDM.assign(0,0, myGTHDM->getQuarks(
QCD::DOWN).getMass());
7331 Md_GTHDM.assign(1,1, myGTHDM->getQuarks(
QCD::STRANGE).getMass());
7332 Md_GTHDM.assign(2,2, myGTHDM->getQuarks(
QCD::BOTTOM).getMass());
7338 if(myGTHDM->getATHDMflag() ==
true)
7342 sigmal_ATHDM = myGTHDM->getNl_11()/myGTHDM->getLeptons(
StandardModel::TAU).getMass();
7344 Nu_GTHDM.assign(0,0, sigmau_ATHDM*Mu_GTHDM(0,0));
7345 Nu_GTHDM.assign(1,1, sigmau_ATHDM*Mu_GTHDM(1,1));
7346 Nu_GTHDM.assign(2,2, sigmau_ATHDM*Mu_GTHDM(2,2));
7348 Nd_GTHDM.assign(0,0, sigmad_ATHDM*Mu_GTHDM(0,0));
7349 Nd_GTHDM.assign(1,1, sigmad_ATHDM*Md_GTHDM(1,1));
7350 Nd_GTHDM.assign(2,2, sigmad_ATHDM*Mu_GTHDM(2,2));
7352 Nl_GTHDM.assign(0,0, sigmal_ATHDM*Ml_GTHDM(0,0));
7353 Nl_GTHDM.assign(1,1, sigmal_ATHDM*Ml_GTHDM(1,1));
7354 Nl_GTHDM.assign(2,2, sigmal_ATHDM*Ml_GTHDM(2,2));
7358 Nu_GTHDM.assign(0,0, myGTHDM->getNu_11());
7359 Nu_GTHDM.assign(0,1, myGTHDM->getNu_12());
7360 Nu_GTHDM.assign(0,2, myGTHDM->getNu_13());
7361 Nu_GTHDM.assign(1,0, myGTHDM->getNu_21());
7362 Nu_GTHDM.assign(1,1, myGTHDM->getNu_22());
7363 Nu_GTHDM.assign(1,2, myGTHDM->getNu_23());
7364 Nu_GTHDM.assign(2,0, myGTHDM->getNu_31());
7365 Nu_GTHDM.assign(2,1, myGTHDM->getNu_32());
7366 Nu_GTHDM.assign(2,2, myGTHDM->getNu_33());
7368 Nd_GTHDM.assign(0,0, myGTHDM->getNd_11());
7369 Nd_GTHDM.assign(0,1, myGTHDM->getNd_12());
7370 Nd_GTHDM.assign(0,2, myGTHDM->getNd_13());
7371 Nd_GTHDM.assign(1,0, myGTHDM->getNd_21());
7372 Nd_GTHDM.assign(1,1, myGTHDM->getNd_22());
7373 Nd_GTHDM.assign(1,2, myGTHDM->getNd_23());
7374 Nd_GTHDM.assign(2,0, myGTHDM->getNd_31());
7375 Nd_GTHDM.assign(2,1, myGTHDM->getNd_32());
7376 Nd_GTHDM.assign(2,2, myGTHDM->getNd_33());
7378 Nl_GTHDM.assign(0,0, myGTHDM->getNl_11());
7379 Nl_GTHDM.assign(0,1, myGTHDM->getNl_12());
7380 Nl_GTHDM.assign(0,2, myGTHDM->getNl_13());
7381 Nl_GTHDM.assign(1,0, myGTHDM->getNl_21());
7382 Nl_GTHDM.assign(1,1, myGTHDM->getNl_22());
7383 Nl_GTHDM.assign(1,2, myGTHDM->getNl_23());
7384 Nl_GTHDM.assign(2,0, myGTHDM->getNl_31());
7385 Nl_GTHDM.assign(2,1, myGTHDM->getNl_32());
7386 Nl_GTHDM.assign(2,2, myGTHDM->getNl_33());
7390 Yu1_GTHDM = (cosb*Mu_GTHDM - sinb*Nu_GTHDM)*
sqrt(2.)/vev;
7391 Yu2_GTHDM = (cosb*Nu_GTHDM + sinb*Mu_GTHDM)*
sqrt(2.)/vev;
7392 Yd1_GTHDM = (cosb*Md_GTHDM - sinb*Nd_GTHDM)*
sqrt(2.)/vev;
7393 Yd2_GTHDM = (cosb*Nd_GTHDM + sinb*Md_GTHDM)*
sqrt(2.)/vev;
7394 Yl1_GTHDM = (cosb*Ml_GTHDM - sinb*Nl_GTHDM)*
sqrt(2.)/vev;
7395 Yl2_GTHDM = (cosb*Nl_GTHDM + sinb*Ml_GTHDM)*
sqrt(2.)/vev;
7398 su = myGTHDM->getNu_11();
7399 sd = myGTHDM->getNd_11();
7400 sl = myGTHDM->getNl_11();
7550 return std::numeric_limits<double>::quiet_NaN();
7554 Q_GTHDM=myGTHDM->getQ_GTHDM();
7555 Ale=myGTHDM->getAle();
7556 Als=myGTHDM->getAlsMz();
7566 MZ=myGTHDM->getMz();
7570 runGeneralTHDMparameters();