a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
NPEffectiveGIMR.cpp
Go to the documentation of this file.
1 /*
2  * Copyright (C) 2014 HEPfit Collaboration
3  *
4  *
5  * For the licensing terms see doc/COPYING.
6  */
7 
8 #include "NPEffectiveGIMR.h"
9 #include <limits>
10 
11 const std::string NPEffectiveGIMR::NPEffectiveGIMRVars[NNPEffectiveGIMRVars]
12  = {"CG", "CW", "CHG", "CHW", "CHB", "CHWB", "CHD", "CHbox", "CH",
13  "CHL1_11", "CHL1_12r", "CHL1_13r", "CHL1_22", "CHL1_23r", "CHL1_33",
14  "CHL1_12i", "CHL1_13i", "CHL1_23i",
15  "CHL3_11", "CHL3_12r", "CHL3_13r", "CHL3_22", "CHL3_23r", "CHL3_33",
16  "CHL3_12i", "CHL3_13i", "CHL3_23i",
17  "CHe_11", "CHe_12r", "CHe_13r", "CHe_22", "CHe_23r", "CHe_33",
18  "CHe_12i", "CHe_13i", "CHe_23i",
19  "CHQ1_11", "CHQ1_12r", "CHQ1_13r", "CHQ1_22", "CHQ1_23r", "CHQ1_33",
20  "CHQ1_12i", "CHQ1_13i", "CHQ1_23i",
21  "CHQ3_11", "CHQ3_12r", "CHQ3_13r", "CHQ3_22", "CHQ3_23r", "CHQ3_33",
22  "CHQ3_12i", "CHQ3_13i", "CHQ3_23i",
23  "CHu_11", "CHu_12r", "CHu_13r", "CHu_22", "CHu_23r", "CHu_33",
24  "CHu_12i", "CHu_13i", "CHu_23i",
25  "CHd_11", "CHd_12r", "CHd_13r", "CHd_22", "CHd_23r", "CHd_33",
26  "CHd_12i", "CHd_13i", "CHd_23i",
27  "CHud_11r", "CHud_12r", "CHud_13r", "CHud_22r", "CHud_23r", "CHud_33r",
28  "CHud_11i", "CHud_12i", "CHud_13i", "CHud_22i", "CHud_23i", "CHud_33i",
29  "CeH_11r", "CeH_12r", "CeH_13r", "CeH_22r", "CeH_23r", "CeH_33r",
30  "CeH_11i", "CeH_12i", "CeH_13i", "CeH_22i", "CeH_23i", "CeH_33i",
31  "CuH_11r", "CuH_12r", "CuH_13r", "CuH_22r", "CuH_23r", "CuH_33r",
32  "CuH_11i", "CuH_12i", "CuH_13i", "CuH_22i", "CuH_23i", "CuH_33i",
33  "CdH_11r", "CdH_12r", "CdH_13r", "CdH_22r", "CdH_23r", "CdH_33r",
34  "CdH_11i", "CdH_12i", "CdH_13i", "CdH_22i", "CdH_23i", "CdH_33i",
35  "CuG_11r", "CuG_12r", "CuG_13r", "CuG_22r", "CuG_23r", "CuG_33r",
36  "CuG_11i", "CuG_12i", "CuG_13i", "CuG_22i", "CuG_23i", "CuG_33i",
37  "CuW_11r", "CuW_12r", "CuW_13r", "CuW_22r", "CuW_23r", "CuW_33r",
38  "CuW_11i", "CuW_12i", "CuW_13i", "CuW_22i", "CuW_23i", "CuW_33i",
39  "CuB_11r", "CuB_12r", "CuB_13r", "CuB_22r", "CuB_23r", "CuB_33r",
40  "CuB_11i", "CuB_12i", "CuB_13i", "CuB_22i", "CuB_23i", "CuB_33i",
41  "CLL_1111","CLL_1221","CLL_2112","CLL_2211","CLL_1122",
42  "CLL_1133","CLL_3311","CLL_1331","CLL_3113",
43  "CLQ1_1111","CLQ1_1122","CLQ1_2211","CLQ1_1221","CLQ1_2112",
44  "CLQ1_1133","CLQ1_3311","CLQ1_1331","CLQ1_3113",
45  "CLQ1_1123","CLQ1_2223","CLQ1_3323",
46  "CLQ1_1132","CLQ1_2232","CLQ1_3332",
47  "CLQ3_1111","CLQ3_1122","CLQ3_2211","CLQ3_1221","CLQ3_2112",
48  "CLQ3_1133","CLQ3_3311","CLQ3_1331","CLQ3_3113",
49  "CLQ3_1123","CLQ3_2223","CLQ3_3323",
50  "CLQ3_1132","CLQ3_2232","CLQ3_3332",
51  "Cee_1111","Cee_1122","Cee_2211","Cee_1133","Cee_3311",
52  "Ceu_1111","Ceu_1122","Ceu_2211","Ceu_1133","Ceu_3311",
53  "Ced_1111","Ced_1122","Ced_2211","Ced_1133","Ced_3311",
54  "Ced_1123","Ced_2223","Ced_3323",
55  "Ced_1132","Ced_2232","Ced_3332",
56  "CLe_1111","CLe_1122","CLe_2211","CLe_1133","CLe_3311",
57  "CLu_1111","CLu_1122","CLu_2211","CLu_1133","CLu_3311",
58  "CLd_1111","CLd_1122","CLd_2211","CLd_1133","CLd_3311",
59  "CLd_1123","CLd_2223","CLd_3323",
60  "CLd_1132","CLd_2232","CLd_3332",
61  "CQe_1111","CQe_1122","CQe_2211","CQe_1133","CQe_3311",
62  "CQe_2311","CQe_2322","CQe_2333",
63  "CQe_3211","CQe_3222","CQe_3233",
64  "Lambda_NP",
65  "eVBF2_HZZ1", "eVBF2_HZZ2", "eVBF2_HZZ3", "eVBF2_HZA1", "eVBF2_HZA2", "eVBF2_HAA",
66  "eVBF2_HWW1", "eVBF2_HWW2", "eVBF2_HWW3", "eVBF2_Hgg", "eVBF2_HZuL", "eVBF2_HZuR",
67  "eVBF2_HZdL", "eVBF2_HZdR", "eVBF2_HWud", "eVBF2_ZuL", "eVBF2_ZuR", "eVBF2_ZdL",
68  "eVBF2_ZdR", "eVBF2_Wud",
69  "eVBF78_HZZ1", "eVBF78_HZZ2", "eVBF78_HZZ3", "eVBF78_HZA1", "eVBF78_HZA2", "eVBF78_HAA",
70  "eVBF78_HWW1", "eVBF78_HWW2", "eVBF78_HWW3", "eVBF78_Hgg", "eVBF78_HZuL", "eVBF78_HZuR",
71  "eVBF78_HZdL", "eVBF78_HZdR", "eVBF78_HWud", "eVBF78_ZuL", "eVBF78_ZuR", "eVBF78_ZdL",
72  "eVBF78_ZdR", "eVBF78_Wud",
73  "eWH2_HWW1", "eWH2_HWW2", "eWH2_HWW3", "eWH2_HWud", "eWH2_Wud",
74  "eWH78_HWW1", "eWH78_HWW2", "eWH78_HWW3", "eWH78_HWud", "eWH78_Wud",
75  "eZH2_HZZ1", "eZH2_HZZ2", "eZH2_HZZ3", "eZH2_HZA1", "eZH2_HZA2", "eZH2_HZuL", "eZH2_HZuR",
76  "eZH2_HZdL", "eZH2_HZdR", "eZH2_ZuL", "eZH2_ZuR", "eZH2_ZdL", "eZH2_ZdR",
77  "eZH78_HZZ1", "eZH78_HZZ2", "eZH78_HZZ3", "eZH78_HZA1", "eZH78_HZA2", "eZH78_HZuL", "eZH78_HZuR",
78  "eZH78_HZdL", "eZH78_HZdR", "eZH78_ZuL", "eZH78_ZuR", "eZH78_ZdL", "eZH78_ZdR",
79  "ettH2_Htt", "ettH2_Hgg",
80  "ettH78_Htt", "ettH78_Hgg"};
81 
82 const std::string NPEffectiveGIMR::NPEffectiveGIMRVars_LFU_QFU[NNPEffectiveGIMRVars_LFU_QFU]
83  = {"CG", "CW", "CHG", "CHW", "CHB", "CHWB", "CHD", "CHbox", "CH",
84  "CHL1", "CHL3", "CHe", "CHQ1", "CHQ3", "CHu", "CHd", "CHud_r", "CHud_i",
85  "CeH_r", "CeH_i", "CuH_r", "CuH_i", "CdH_r", "CdH_i",
86  "CuG_r", "CuG_i", "CuW_r", "CuW_i", "CuB_r", "CuB_i",
87  "CLL", "CLQ1", "CLQ3",
88  "Cee", "Ceu", "Ced", "CLe", "CLu", "CLd", "CQe","Lambda_NP",
89  "eVBF2_HZZ1", "eVBF2_HZZ2", "eVBF2_HZZ3", "eVBF2_HZA1", "eVBF2_HZA2", "eVBF2_HAA",
90  "eVBF2_HWW1", "eVBF2_HWW2", "eVBF2_HWW3", "eVBF2_Hgg", "eVBF2_HZuL", "eVBF2_HZuR",
91  "eVBF2_HZdL", "eVBF2_HZdR", "eVBF2_HWud", "eVBF2_ZuL", "eVBF2_ZuR", "eVBF2_ZdL",
92  "eVBF2_ZdR", "eVBF2_Wud",
93  "eVBF78_HZZ1", "eVBF78_HZZ2", "eVBF78_HZZ3", "eVBF78_HZA1", "eVBF78_HZA2", "eVBF78_HAA",
94  "eVBF78_HWW1", "eVBF78_HWW2", "eVBF78_HWW3", "eVBF78_Hgg", "eVBF78_HZuL", "eVBF78_HZuR",
95  "eVBF78_HZdL", "eVBF78_HZdR", "eVBF78_HWud", "eVBF78_ZuL", "eVBF78_ZuR", "eVBF78_ZdL",
96  "eVBF78_ZdR", "eVBF78_Wud",
97  "eWH2_HWW1", "eWH2_HWW2", "eWH2_HWW3", "eWH2_HWud", "eWH2_Wud",
98  "eWH78_HWW1", "eWH78_HWW2", "eWH78_HWW3", "eWH78_HWud", "eWH78_Wud",
99  "eZH2_HZZ1", "eZH2_HZZ2", "eZH2_HZZ3", "eZH2_HZA1", "eZH2_HZA2", "eZH2_HZuL", "eZH2_HZuR",
100  "eZH2_HZdL", "eZH2_HZdR", "eZH2_ZuL", "eZH2_ZuR", "eZH2_ZdL", "eZH2_ZdR",
101  "eZH78_HZZ1", "eZH78_HZZ2", "eZH78_HZZ3", "eZH78_HZA1", "eZH78_HZA2", "eZH78_HZuL", "eZH78_HZuR",
102  "eZH78_HZdL", "eZH78_HZdR", "eZH78_ZuL", "eZH78_ZuR", "eZH78_ZdL", "eZH78_ZdR",
103  "ettH2_Htt", "ettH2_Hgg",
104  "ettH78_Htt", "ettH78_Hgg"};
105 
106 NPEffectiveGIMR::NPEffectiveGIMR(const bool FlagLeptonUniversal_in, const bool FlagQuarkUniversal_in)
107 : NPbase(), myLEP2GIMR(trueSM), FlagLeptonUniversal(FlagLeptonUniversal_in), FlagQuarkUniversal(FlagQuarkUniversal_in)
108 {
111  throw std::runtime_error("Invalid arguments for NPEffectiveGIMR::NPEffectiveGIMR()");
112 
113  FlagMwInput = false;
114  FlagQuadraticTerms = false;
116 
117  ModelParamMap.insert(std::make_pair("CG", std::cref(CG)));
118  ModelParamMap.insert(std::make_pair("CW", std::cref(CW)));
119  ModelParamMap.insert(std::make_pair("CHG", std::cref(CHG)));
120  ModelParamMap.insert(std::make_pair("CHW", std::cref(CHW)));
121  ModelParamMap.insert(std::make_pair("CHB", std::cref(CHB)));
122  ModelParamMap.insert(std::make_pair("CHWB", std::cref(CHWB)));
123  ModelParamMap.insert(std::make_pair("CHD", std::cref(CHD)));
124  ModelParamMap.insert(std::make_pair("CHbox", std::cref(CHbox)));
125  ModelParamMap.insert(std::make_pair("CH", std::cref(CH)));
126  if (FlagLeptonUniversal) {
127  ModelParamMap.insert(std::make_pair("CHL1", std::cref(CHL1_11)));
128  ModelParamMap.insert(std::make_pair("CHL3", std::cref(CHL3_11)));
129  ModelParamMap.insert(std::make_pair("CHe", std::cref(CHe_11)));
130  ModelParamMap.insert(std::make_pair("CeH_r", std::cref(CeH_11r)));
131  ModelParamMap.insert(std::make_pair("CeH_i", std::cref(CeH_11i)));
132  ModelParamMap.insert(std::make_pair("CLL", std::cref(CLL_1221)));
133  ModelParamMap.insert(std::make_pair("Cee", std::cref(Cee_1111)));
134  ModelParamMap.insert(std::make_pair("CLe", std::cref(CLe_1111)));
135  } else {
136  ModelParamMap.insert(std::make_pair("CHL1_11", std::cref(CHL1_11)));
137  ModelParamMap.insert(std::make_pair("CHL1_12r", std::cref(CHL1_12r)));
138  ModelParamMap.insert(std::make_pair("CHL1_13r", std::cref(CHL1_13r)));
139  ModelParamMap.insert(std::make_pair("CHL1_22", std::cref(CHL1_22)));
140  ModelParamMap.insert(std::make_pair("CHL1_23r", std::cref(CHL1_23r)));
141  ModelParamMap.insert(std::make_pair("CHL1_33", std::cref(CHL1_33)));
142  ModelParamMap.insert(std::make_pair("CHL1_12i", std::cref(CHL1_12i)));
143  ModelParamMap.insert(std::make_pair("CHL1_13i", std::cref(CHL1_13i)));
144  ModelParamMap.insert(std::make_pair("CHL1_23i", std::cref(CHL1_23i)));
145  ModelParamMap.insert(std::make_pair("CHL3_11", std::cref(CHL3_11)));
146  ModelParamMap.insert(std::make_pair("CHL3_12r", std::cref(CHL3_12r)));
147  ModelParamMap.insert(std::make_pair("CHL3_13r", std::cref(CHL3_13r)));
148  ModelParamMap.insert(std::make_pair("CHL3_22", std::cref(CHL3_22)));
149  ModelParamMap.insert(std::make_pair("CHL3_23r", std::cref(CHL3_23r)));
150  ModelParamMap.insert(std::make_pair("CHL3_33", std::cref(CHL3_33)));
151  ModelParamMap.insert(std::make_pair("CHL3_12i", std::cref(CHL3_12i)));
152  ModelParamMap.insert(std::make_pair("CHL3_13i", std::cref(CHL3_13i)));
153  ModelParamMap.insert(std::make_pair("CHL3_23i", std::cref(CHL3_23i)));
154  ModelParamMap.insert(std::make_pair("CHe_11", std::cref(CHe_11)));
155  ModelParamMap.insert(std::make_pair("CHe_12r", std::cref(CHe_12r)));
156  ModelParamMap.insert(std::make_pair("CHe_13r", std::cref(CHe_13r)));
157  ModelParamMap.insert(std::make_pair("CHe_22", std::cref(CHe_22)));
158  ModelParamMap.insert(std::make_pair("CHe_23r", std::cref(CHe_23r)));
159  ModelParamMap.insert(std::make_pair("CHe_33", std::cref(CHe_33)));
160  ModelParamMap.insert(std::make_pair("CHe_12i", std::cref(CHe_12i)));
161  ModelParamMap.insert(std::make_pair("CHe_13i", std::cref(CHe_13i)));
162  ModelParamMap.insert(std::make_pair("CHe_23i", std::cref(CHe_23i)));
163  ModelParamMap.insert(std::make_pair("CeH_11r", std::cref(CeH_11r)));
164  ModelParamMap.insert(std::make_pair("CeH_12r", std::cref(CeH_12r)));
165  ModelParamMap.insert(std::make_pair("CeH_13r", std::cref(CeH_13r)));
166  ModelParamMap.insert(std::make_pair("CeH_22r", std::cref(CeH_22r)));
167  ModelParamMap.insert(std::make_pair("CeH_23r", std::cref(CeH_23r)));
168  ModelParamMap.insert(std::make_pair("CeH_33r", std::cref(CeH_33r)));
169  ModelParamMap.insert(std::make_pair("CeH_11i", std::cref(CeH_11i)));
170  ModelParamMap.insert(std::make_pair("CeH_12i", std::cref(CeH_12i)));
171  ModelParamMap.insert(std::make_pair("CeH_13i", std::cref(CeH_13i)));
172  ModelParamMap.insert(std::make_pair("CeH_22i", std::cref(CeH_22i)));
173  ModelParamMap.insert(std::make_pair("CeH_23i", std::cref(CeH_23i)));
174  ModelParamMap.insert(std::make_pair("CeH_33i", std::cref(CeH_33i)));
175  ModelParamMap.insert(std::make_pair("CLL_1111", std::cref(CLL_1111)));
176  ModelParamMap.insert(std::make_pair("CLL_1221", std::cref(CLL_1221)));
177  ModelParamMap.insert(std::make_pair("CLL_2112", std::cref(CLL_2112)));
178  ModelParamMap.insert(std::make_pair("CLL_2211", std::cref(CLL_2211)));
179  ModelParamMap.insert(std::make_pair("CLL_1122", std::cref(CLL_1122)));
180  ModelParamMap.insert(std::make_pair("CLL_1331", std::cref(CLL_1331)));
181  ModelParamMap.insert(std::make_pair("CLL_3113", std::cref(CLL_3113)));
182  ModelParamMap.insert(std::make_pair("CLL_3311", std::cref(CLL_3311)));
183  ModelParamMap.insert(std::make_pair("CLL_1133", std::cref(CLL_1133)));
184  ModelParamMap.insert(std::make_pair("Cee_1111", std::cref(Cee_1111)));
185  ModelParamMap.insert(std::make_pair("Cee_1122", std::cref(Cee_1122)));
186  ModelParamMap.insert(std::make_pair("Cee_2211", std::cref(Cee_2211)));
187  ModelParamMap.insert(std::make_pair("Cee_1133", std::cref(Cee_1133)));
188  ModelParamMap.insert(std::make_pair("Cee_3311", std::cref(Cee_3311)));
189  ModelParamMap.insert(std::make_pair("CLe_1111", std::cref(CLe_1111)));
190  ModelParamMap.insert(std::make_pair("CLe_1122", std::cref(CLe_1122)));
191  ModelParamMap.insert(std::make_pair("CLe_2211", std::cref(CLe_2211)));
192  ModelParamMap.insert(std::make_pair("CLe_1133", std::cref(CLe_1133)));
193  ModelParamMap.insert(std::make_pair("CLe_3311", std::cref(CLe_3311)));
194  ModelParamMap.insert(std::make_pair("CLQ1_1123", std::cref(CLQ1_1123)));
195  ModelParamMap.insert(std::make_pair("CLQ1_2223", std::cref(CLQ1_2223)));
196  ModelParamMap.insert(std::make_pair("CLQ1_3323", std::cref(CLQ1_3323)));
197  ModelParamMap.insert(std::make_pair("CLQ1_1132", std::cref(CLQ1_1132)));
198  ModelParamMap.insert(std::make_pair("CLQ1_2232", std::cref(CLQ1_2232)));
199  ModelParamMap.insert(std::make_pair("CLQ1_3332", std::cref(CLQ1_3332)));
200  ModelParamMap.insert(std::make_pair("CLQ3_1123", std::cref(CLQ3_1123)));
201  ModelParamMap.insert(std::make_pair("CLQ3_2223", std::cref(CLQ3_2223)));
202  ModelParamMap.insert(std::make_pair("CLQ3_3323", std::cref(CLQ3_3323)));
203  ModelParamMap.insert(std::make_pair("CLQ3_1132", std::cref(CLQ3_1132)));
204  ModelParamMap.insert(std::make_pair("CLQ3_2232", std::cref(CLQ3_2232)));
205  ModelParamMap.insert(std::make_pair("CLQ3_3332", std::cref(CLQ3_3332)));
206  ModelParamMap.insert(std::make_pair("Ced_1123", std::cref(Ced_1123)));
207  ModelParamMap.insert(std::make_pair("Ced_2223", std::cref(Ced_2223)));
208  ModelParamMap.insert(std::make_pair("Ced_3323", std::cref(Ced_3323)));
209  ModelParamMap.insert(std::make_pair("Ced_1132", std::cref(Ced_1132)));
210  ModelParamMap.insert(std::make_pair("Ced_2232", std::cref(Ced_2232)));
211  ModelParamMap.insert(std::make_pair("Ced_3332", std::cref(Ced_3332)));
212  ModelParamMap.insert(std::make_pair("CLd_1123", std::cref(CLd_1123)));
213  ModelParamMap.insert(std::make_pair("CLd_2223", std::cref(CLd_2223)));
214  ModelParamMap.insert(std::make_pair("CLd_3323", std::cref(CLd_3323)));
215  ModelParamMap.insert(std::make_pair("CLd_1132", std::cref(CLd_1132)));
216  ModelParamMap.insert(std::make_pair("CLd_2232", std::cref(CLd_2232)));
217  ModelParamMap.insert(std::make_pair("CLd_3332", std::cref(CLd_3332)));
218  ModelParamMap.insert(std::make_pair("CQe_2311", std::cref(CQe_2311)));
219  ModelParamMap.insert(std::make_pair("CQe_2322", std::cref(CQe_2322)));
220  ModelParamMap.insert(std::make_pair("CQe_2333", std::cref(CQe_2333)));
221  ModelParamMap.insert(std::make_pair("CQe_3211", std::cref(CQe_3211)));
222  ModelParamMap.insert(std::make_pair("CQe_3222", std::cref(CQe_3222)));
223  ModelParamMap.insert(std::make_pair("CQe_3233", std::cref(CQe_3233)));
224  }
225  if (FlagQuarkUniversal) {
226  ModelParamMap.insert(std::make_pair("CHQ1", std::cref(CHQ1_11)));
227  ModelParamMap.insert(std::make_pair("CHQ3", std::cref(CHQ3_11)));
228  ModelParamMap.insert(std::make_pair("CHu", std::cref(CHu_11)));
229  ModelParamMap.insert(std::make_pair("CHd", std::cref(CHd_11)));
230  ModelParamMap.insert(std::make_pair("CHud_r", std::cref(CHud_11r)));
231  ModelParamMap.insert(std::make_pair("CHud_i", std::cref(CHud_11i)));
232  ModelParamMap.insert(std::make_pair("CuH_r", std::cref(CuH_11r)));
233  ModelParamMap.insert(std::make_pair("CuH_i", std::cref(CuH_11i)));
234  ModelParamMap.insert(std::make_pair("CdH_r", std::cref(CdH_11r)));
235  ModelParamMap.insert(std::make_pair("CdH_i", std::cref(CdH_11i)));
236  ModelParamMap.insert(std::make_pair("CuG_r", std::cref(CuG_11r)));
237  ModelParamMap.insert(std::make_pair("CuG_i", std::cref(CuG_11i)));
238  ModelParamMap.insert(std::make_pair("CuW_r", std::cref(CuW_11r)));
239  ModelParamMap.insert(std::make_pair("CuW_i", std::cref(CuW_11i)));
240  ModelParamMap.insert(std::make_pair("CuB_r", std::cref(CuB_11r)));
241  ModelParamMap.insert(std::make_pair("CuB_i", std::cref(CuB_11i)));
242  } else {
243  ModelParamMap.insert(std::make_pair("CHQ1_11", std::cref(CHQ1_11)));
244  ModelParamMap.insert(std::make_pair("CHQ1_12r", std::cref(CHQ1_12r)));
245  ModelParamMap.insert(std::make_pair("CHQ1_13r", std::cref(CHQ1_13r)));
246  ModelParamMap.insert(std::make_pair("CHQ1_22", std::cref(CHQ1_22)));
247  ModelParamMap.insert(std::make_pair("CHQ1_23r", std::cref(CHQ1_23r)));
248  ModelParamMap.insert(std::make_pair("CHQ1_33", std::cref(CHQ1_33)));
249  ModelParamMap.insert(std::make_pair("CHQ1_12i", std::cref(CHQ1_12i)));
250  ModelParamMap.insert(std::make_pair("CHQ1_13i", std::cref(CHQ1_13i)));
251  ModelParamMap.insert(std::make_pair("CHQ1_23i", std::cref(CHQ1_23i)));
252  ModelParamMap.insert(std::make_pair("CHQ3_11", std::cref(CHQ3_11)));
253  ModelParamMap.insert(std::make_pair("CHQ3_12r", std::cref(CHQ3_12r)));
254  ModelParamMap.insert(std::make_pair("CHQ3_13r", std::cref(CHQ3_13r)));
255  ModelParamMap.insert(std::make_pair("CHQ3_22", std::cref(CHQ3_22)));
256  ModelParamMap.insert(std::make_pair("CHQ3_23r", std::cref(CHQ3_23r)));
257  ModelParamMap.insert(std::make_pair("CHQ3_33", std::cref(CHQ3_33)));
258  ModelParamMap.insert(std::make_pair("CHQ3_12i", std::cref(CHQ3_12i)));
259  ModelParamMap.insert(std::make_pair("CHQ3_13i", std::cref(CHQ3_13i)));
260  ModelParamMap.insert(std::make_pair("CHQ3_23i", std::cref(CHQ3_23i)));
261  ModelParamMap.insert(std::make_pair("CHu_11", std::cref(CHu_11)));
262  ModelParamMap.insert(std::make_pair("CHu_12r", std::cref(CHu_12r)));
263  ModelParamMap.insert(std::make_pair("CHu_13r", std::cref(CHu_13r)));
264  ModelParamMap.insert(std::make_pair("CHu_22", std::cref(CHu_22)));
265  ModelParamMap.insert(std::make_pair("CHu_23r", std::cref(CHu_23r)));
266  ModelParamMap.insert(std::make_pair("CHu_33", std::cref(CHu_33)));
267  ModelParamMap.insert(std::make_pair("CHu_12i", std::cref(CHu_12i)));
268  ModelParamMap.insert(std::make_pair("CHu_13i", std::cref(CHu_13i)));
269  ModelParamMap.insert(std::make_pair("CHu_23i", std::cref(CHu_23i)));
270  ModelParamMap.insert(std::make_pair("CHd_11", std::cref(CHd_11)));
271  ModelParamMap.insert(std::make_pair("CHd_12r", std::cref(CHd_12r)));
272  ModelParamMap.insert(std::make_pair("CHd_13r", std::cref(CHd_13r)));
273  ModelParamMap.insert(std::make_pair("CHd_22", std::cref(CHd_22)));
274  ModelParamMap.insert(std::make_pair("CHd_23r", std::cref(CHd_23r)));
275  ModelParamMap.insert(std::make_pair("CHd_33", std::cref(CHd_33)));
276  ModelParamMap.insert(std::make_pair("CHd_12i", std::cref(CHd_12i)));
277  ModelParamMap.insert(std::make_pair("CHd_13i", std::cref(CHd_13i)));
278  ModelParamMap.insert(std::make_pair("CHd_23i", std::cref(CHd_23i)));
279  ModelParamMap.insert(std::make_pair("CHud_11r", std::cref(CHud_11r)));
280  ModelParamMap.insert(std::make_pair("CHud_12r", std::cref(CHud_12r)));
281  ModelParamMap.insert(std::make_pair("CHud_13r", std::cref(CHud_13r)));
282  ModelParamMap.insert(std::make_pair("CHud_22r", std::cref(CHud_22r)));
283  ModelParamMap.insert(std::make_pair("CHud_23r", std::cref(CHud_23r)));
284  ModelParamMap.insert(std::make_pair("CHud_33r", std::cref(CHud_33r)));
285  ModelParamMap.insert(std::make_pair("CHud_11i", std::cref(CHud_11i)));
286  ModelParamMap.insert(std::make_pair("CHud_12i", std::cref(CHud_12i)));
287  ModelParamMap.insert(std::make_pair("CHud_13i", std::cref(CHud_13i)));
288  ModelParamMap.insert(std::make_pair("CHud_22i", std::cref(CHud_22i)));
289  ModelParamMap.insert(std::make_pair("CHud_23i", std::cref(CHud_23i)));
290  ModelParamMap.insert(std::make_pair("CHud_33i", std::cref(CHud_33i)));
291  ModelParamMap.insert(std::make_pair("CuH_11r", std::cref(CuH_11r)));
292  ModelParamMap.insert(std::make_pair("CuH_12r", std::cref(CuH_12r)));
293  ModelParamMap.insert(std::make_pair("CuH_13r", std::cref(CuH_13r)));
294  ModelParamMap.insert(std::make_pair("CuH_22r", std::cref(CuH_22r)));
295  ModelParamMap.insert(std::make_pair("CuH_23r", std::cref(CuH_23r)));
296  ModelParamMap.insert(std::make_pair("CuH_33r", std::cref(CuH_33r)));
297  ModelParamMap.insert(std::make_pair("CuH_11i", std::cref(CuH_11i)));
298  ModelParamMap.insert(std::make_pair("CuH_12i", std::cref(CuH_12i)));
299  ModelParamMap.insert(std::make_pair("CuH_13i", std::cref(CuH_13i)));
300  ModelParamMap.insert(std::make_pair("CuH_22i", std::cref(CuH_22i)));
301  ModelParamMap.insert(std::make_pair("CuH_23i", std::cref(CuH_23i)));
302  ModelParamMap.insert(std::make_pair("CuH_33i", std::cref(CuH_33i)));
303  ModelParamMap.insert(std::make_pair("CdH_11r", std::cref(CdH_11r)));
304  ModelParamMap.insert(std::make_pair("CdH_12r", std::cref(CdH_12r)));
305  ModelParamMap.insert(std::make_pair("CdH_13r", std::cref(CdH_13r)));
306  ModelParamMap.insert(std::make_pair("CdH_22r", std::cref(CdH_22r)));
307  ModelParamMap.insert(std::make_pair("CdH_23r", std::cref(CdH_23r)));
308  ModelParamMap.insert(std::make_pair("CdH_33r", std::cref(CdH_33r)));
309  ModelParamMap.insert(std::make_pair("CdH_11i", std::cref(CdH_11i)));
310  ModelParamMap.insert(std::make_pair("CdH_12i", std::cref(CdH_12i)));
311  ModelParamMap.insert(std::make_pair("CdH_13i", std::cref(CdH_13i)));
312  ModelParamMap.insert(std::make_pair("CdH_22i", std::cref(CdH_22i)));
313  ModelParamMap.insert(std::make_pair("CdH_23i", std::cref(CdH_23i)));
314  ModelParamMap.insert(std::make_pair("CdH_33i", std::cref(CdH_33i)));
315  ModelParamMap.insert(std::make_pair("CuG_11r", std::cref(CuG_11r)));
316  ModelParamMap.insert(std::make_pair("CuG_12r", std::cref(CuG_12r)));
317  ModelParamMap.insert(std::make_pair("CuG_13r", std::cref(CuG_13r)));
318  ModelParamMap.insert(std::make_pair("CuG_22r", std::cref(CuG_22r)));
319  ModelParamMap.insert(std::make_pair("CuG_23r", std::cref(CuG_23r)));
320  ModelParamMap.insert(std::make_pair("CuG_33r", std::cref(CuG_33r)));
321  ModelParamMap.insert(std::make_pair("CuG_11i", std::cref(CuG_11i)));
322  ModelParamMap.insert(std::make_pair("CuG_12i", std::cref(CuG_12i)));
323  ModelParamMap.insert(std::make_pair("CuG_13i", std::cref(CuG_13i)));
324  ModelParamMap.insert(std::make_pair("CuG_22i", std::cref(CuG_22i)));
325  ModelParamMap.insert(std::make_pair("CuG_23i", std::cref(CuG_23i)));
326  ModelParamMap.insert(std::make_pair("CuG_33i", std::cref(CuG_33i)));
327  ModelParamMap.insert(std::make_pair("CuW_11r", std::cref(CuW_11r)));
328  ModelParamMap.insert(std::make_pair("CuW_12r", std::cref(CuW_12r)));
329  ModelParamMap.insert(std::make_pair("CuW_13r", std::cref(CuW_13r)));
330  ModelParamMap.insert(std::make_pair("CuW_22r", std::cref(CuW_22r)));
331  ModelParamMap.insert(std::make_pair("CuW_23r", std::cref(CuW_23r)));
332  ModelParamMap.insert(std::make_pair("CuW_33r", std::cref(CuW_33r)));
333  ModelParamMap.insert(std::make_pair("CuW_11i", std::cref(CuW_11i)));
334  ModelParamMap.insert(std::make_pair("CuW_12i", std::cref(CuW_12i)));
335  ModelParamMap.insert(std::make_pair("CuW_13i", std::cref(CuW_13i)));
336  ModelParamMap.insert(std::make_pair("CuW_22i", std::cref(CuW_22i)));
337  ModelParamMap.insert(std::make_pair("CuW_23i", std::cref(CuW_23i)));
338  ModelParamMap.insert(std::make_pair("CuW_33i", std::cref(CuW_33i)));
339  ModelParamMap.insert(std::make_pair("CuB_11r", std::cref(CuB_11r)));
340  ModelParamMap.insert(std::make_pair("CuB_12r", std::cref(CuB_12r)));
341  ModelParamMap.insert(std::make_pair("CuB_13r", std::cref(CuB_13r)));
342  ModelParamMap.insert(std::make_pair("CuB_22r", std::cref(CuB_22r)));
343  ModelParamMap.insert(std::make_pair("CuB_23r", std::cref(CuB_23r)));
344  ModelParamMap.insert(std::make_pair("CuB_33r", std::cref(CuB_33r)));
345  ModelParamMap.insert(std::make_pair("CuB_11i", std::cref(CuB_11i)));
346  ModelParamMap.insert(std::make_pair("CuB_12i", std::cref(CuB_12i)));
347  ModelParamMap.insert(std::make_pair("CuB_13i", std::cref(CuB_13i)));
348  ModelParamMap.insert(std::make_pair("CuB_22i", std::cref(CuB_22i)));
349  ModelParamMap.insert(std::make_pair("CuB_23i", std::cref(CuB_23i)));
350  ModelParamMap.insert(std::make_pair("CuB_33i", std::cref(CuB_33i)));
351  }
353  ModelParamMap.insert(std::make_pair("CLQ1", std::cref(CLQ1_1111)));
354  ModelParamMap.insert(std::make_pair("CLQ3", std::cref(CLQ3_1111)));
355  ModelParamMap.insert(std::make_pair("Ceu", std::cref(Ceu_1111)));
356  ModelParamMap.insert(std::make_pair("Ced", std::cref(Ced_1111)));
357  ModelParamMap.insert(std::make_pair("CLu", std::cref(CLu_1111)));
358  ModelParamMap.insert(std::make_pair("CLd", std::cref(CLd_1111)));
359  ModelParamMap.insert(std::make_pair("CQe", std::cref(CQe_1111)));
360  } else {
361  ModelParamMap.insert(std::make_pair("CLQ1_1111", std::cref(CLQ1_1111)));
362  ModelParamMap.insert(std::make_pair("CLQ1_1122", std::cref(CLQ1_1122)));
363  ModelParamMap.insert(std::make_pair("CLQ1_2211", std::cref(CLQ1_2211)));
364  ModelParamMap.insert(std::make_pair("CLQ1_1221", std::cref(CLQ1_1221)));
365  ModelParamMap.insert(std::make_pair("CLQ1_2112", std::cref(CLQ1_2112)));
366  ModelParamMap.insert(std::make_pair("CLQ1_1133", std::cref(CLQ1_1133)));
367  ModelParamMap.insert(std::make_pair("CLQ1_3311", std::cref(CLQ1_3311)));
368  ModelParamMap.insert(std::make_pair("CLQ1_1331", std::cref(CLQ1_1331)));
369  ModelParamMap.insert(std::make_pair("CLQ1_3113", std::cref(CLQ1_3113)));
370  ModelParamMap.insert(std::make_pair("CLQ3_1111", std::cref(CLQ3_1111)));
371  ModelParamMap.insert(std::make_pair("CLQ3_1122", std::cref(CLQ3_1122)));
372  ModelParamMap.insert(std::make_pair("CLQ3_2211", std::cref(CLQ3_2211)));
373  ModelParamMap.insert(std::make_pair("CLQ3_1221", std::cref(CLQ3_1221)));
374  ModelParamMap.insert(std::make_pair("CLQ3_2112", std::cref(CLQ3_2112)));
375  ModelParamMap.insert(std::make_pair("CLQ3_1133", std::cref(CLQ3_1133)));
376  ModelParamMap.insert(std::make_pair("CLQ3_3311", std::cref(CLQ3_3311)));
377  ModelParamMap.insert(std::make_pair("CLQ3_1331", std::cref(CLQ3_1331)));
378  ModelParamMap.insert(std::make_pair("CLQ3_3113", std::cref(CLQ3_3113)));
379  ModelParamMap.insert(std::make_pair("Ceu_1111", std::cref(Ceu_1111)));
380  ModelParamMap.insert(std::make_pair("Ceu_1122", std::cref(Ceu_1122)));
381  ModelParamMap.insert(std::make_pair("Ceu_2211", std::cref(Ceu_2211)));
382  ModelParamMap.insert(std::make_pair("Ceu_1133", std::cref(Ceu_1133)));
383  ModelParamMap.insert(std::make_pair("Ceu_3311", std::cref(Ceu_3311)));
384  ModelParamMap.insert(std::make_pair("Ced_1111", std::cref(Ced_1111)));
385  ModelParamMap.insert(std::make_pair("Ced_1122", std::cref(Ced_1122)));
386  ModelParamMap.insert(std::make_pair("Ced_2211", std::cref(Ced_2211)));
387  ModelParamMap.insert(std::make_pair("Ced_1133", std::cref(Ced_1133)));
388  ModelParamMap.insert(std::make_pair("Ced_3311", std::cref(Ced_3311)));
389  ModelParamMap.insert(std::make_pair("CLe_1111", std::cref(CLe_1111)));
390  ModelParamMap.insert(std::make_pair("CLe_1122", std::cref(CLe_1122)));
391  ModelParamMap.insert(std::make_pair("CLe_2211", std::cref(CLe_2211)));
392  ModelParamMap.insert(std::make_pair("CLe_1133", std::cref(CLe_1133)));
393  ModelParamMap.insert(std::make_pair("CLe_3311", std::cref(CLe_3311)));
394  ModelParamMap.insert(std::make_pair("CLu_1111", std::cref(CLu_1111)));
395  ModelParamMap.insert(std::make_pair("CLu_1122", std::cref(CLu_1122)));
396  ModelParamMap.insert(std::make_pair("CLu_2211", std::cref(CLu_2211)));
397  ModelParamMap.insert(std::make_pair("CLu_1133", std::cref(CLu_1133)));
398  ModelParamMap.insert(std::make_pair("CLu_3311", std::cref(CLu_3311)));
399  ModelParamMap.insert(std::make_pair("CLd_1111", std::cref(CLd_1111)));
400  ModelParamMap.insert(std::make_pair("CLd_1122", std::cref(CLd_1122)));
401  ModelParamMap.insert(std::make_pair("CLd_2211", std::cref(CLd_2211)));
402  ModelParamMap.insert(std::make_pair("CLd_1133", std::cref(CLd_1133)));
403  ModelParamMap.insert(std::make_pair("CLd_3311", std::cref(CLd_3311)));
404  ModelParamMap.insert(std::make_pair("CQe_1111", std::cref(CQe_1111)));
405  ModelParamMap.insert(std::make_pair("CQe_1122", std::cref(CQe_1122)));
406  ModelParamMap.insert(std::make_pair("CQe_2211", std::cref(CQe_2211)));
407  ModelParamMap.insert(std::make_pair("CQe_1133", std::cref(CQe_1133)));
408  ModelParamMap.insert(std::make_pair("CQe_3311", std::cref(CQe_3311)));
409  }
410  ModelParamMap.insert(std::make_pair("Lambda_NP", std::cref(Lambda_NP)));
411  ModelParamMap.insert(std::make_pair("eVBF2_HZZ1", std::cref(eVBF2_HZZ1)));
412  ModelParamMap.insert(std::make_pair("eVBF2_HZZ2", std::cref(eVBF2_HZZ2)));
413  ModelParamMap.insert(std::make_pair("eVBF2_HZZ3", std::cref(eVBF2_HZZ3)));
414  ModelParamMap.insert(std::make_pair("eVBF2_HZA1", std::cref(eVBF2_HZA1)));
415  ModelParamMap.insert(std::make_pair("eVBF2_HZA2", std::cref(eVBF2_HZA2)));
416  ModelParamMap.insert(std::make_pair("eVBF2_HAA", std::cref(eVBF2_HAA)));
417  ModelParamMap.insert(std::make_pair("eVBF2_HWW1", std::cref(eVBF2_HWW1)));
418  ModelParamMap.insert(std::make_pair("eVBF2_HWW2", std::cref(eVBF2_HWW2)));
419  ModelParamMap.insert(std::make_pair("eVBF2_HWW3", std::cref(eVBF2_HWW3)));
420  ModelParamMap.insert(std::make_pair("eVBF2_Hgg", std::cref(eVBF2_Hgg)));
421  ModelParamMap.insert(std::make_pair("eVBF2_HZuL", std::cref(eVBF2_HZuL)));
422  ModelParamMap.insert(std::make_pair("eVBF2_HZuR", std::cref(eVBF2_HZuR)));
423  ModelParamMap.insert(std::make_pair("eVBF2_HZdL", std::cref(eVBF2_HZdL)));
424  ModelParamMap.insert(std::make_pair("eVBF2_HZdR", std::cref(eVBF2_HZdR)));
425  ModelParamMap.insert(std::make_pair("eVBF2_HWud", std::cref(eVBF2_HWud)));
426  ModelParamMap.insert(std::make_pair("eVBF2_ZuL", std::cref(eVBF2_ZuL)));
427  ModelParamMap.insert(std::make_pair("eVBF2_ZuR", std::cref(eVBF2_ZuR)));
428  ModelParamMap.insert(std::make_pair("eVBF2_ZdL", std::cref(eVBF2_ZdL)));
429  ModelParamMap.insert(std::make_pair("eVBF2_ZdR", std::cref(eVBF2_ZdR)));
430  ModelParamMap.insert(std::make_pair("eVBF2_Wud", std::cref(eVBF2_Wud)));
431  ModelParamMap.insert(std::make_pair("eVBF78_HZZ1", std::cref(eVBF78_HZZ1)));
432  ModelParamMap.insert(std::make_pair("eVBF78_HZZ2", std::cref(eVBF78_HZZ2)));
433  ModelParamMap.insert(std::make_pair("eVBF78_HZZ3", std::cref(eVBF78_HZZ3)));
434  ModelParamMap.insert(std::make_pair("eVBF78_HZA1", std::cref(eVBF78_HZA1)));
435  ModelParamMap.insert(std::make_pair("eVBF78_HZA2", std::cref(eVBF78_HZA2)));
436  ModelParamMap.insert(std::make_pair("eVBF78_HAA", std::cref(eVBF78_HAA)));
437  ModelParamMap.insert(std::make_pair("eVBF78_HWW1", std::cref(eVBF78_HWW1)));
438  ModelParamMap.insert(std::make_pair("eVBF78_HWW2", std::cref(eVBF78_HWW2)));
439  ModelParamMap.insert(std::make_pair("eVBF78_HWW3", std::cref(eVBF78_HWW3)));
440  ModelParamMap.insert(std::make_pair("eVBF78_Hgg", std::cref(eVBF78_Hgg)));
441  ModelParamMap.insert(std::make_pair("eVBF78_HZuL", std::cref(eVBF78_HZuL)));
442  ModelParamMap.insert(std::make_pair("eVBF78_HZuR", std::cref(eVBF78_HZuR)));
443  ModelParamMap.insert(std::make_pair("eVBF78_HZdL", std::cref(eVBF78_HZdL)));
444  ModelParamMap.insert(std::make_pair("eVBF78_HZdR", std::cref(eVBF78_HZdR)));
445  ModelParamMap.insert(std::make_pair("eVBF78_HWud", std::cref(eVBF78_HWud)));
446  ModelParamMap.insert(std::make_pair("eVBF78_ZuL", std::cref(eVBF78_ZuL)));
447  ModelParamMap.insert(std::make_pair("eVBF78_ZuR", std::cref(eVBF78_ZuR)));
448  ModelParamMap.insert(std::make_pair("eVBF78_ZdL", std::cref(eVBF78_ZdL)));
449  ModelParamMap.insert(std::make_pair("eVBF78_ZdR", std::cref(eVBF78_ZdR)));
450  ModelParamMap.insert(std::make_pair("eVBF78_Wud", std::cref(eVBF78_Wud)));
451  ModelParamMap.insert(std::make_pair("eWH2_HWW1", std::cref(eWH2_HWW1)));
452  ModelParamMap.insert(std::make_pair("eWH2_HWW2", std::cref(eWH2_HWW2)));
453  ModelParamMap.insert(std::make_pair("eWH2_HWW3", std::cref(eWH2_HWW3)));
454  ModelParamMap.insert(std::make_pair("eWH2_HWud", std::cref(eWH2_HWud)));
455  ModelParamMap.insert(std::make_pair("eWH2_Wud", std::cref(eWH2_Wud)));
456  ModelParamMap.insert(std::make_pair("eWH78_HWW1", std::cref(eWH78_HWW1)));
457  ModelParamMap.insert(std::make_pair("eWH78_HWW2", std::cref(eWH78_HWW2)));
458  ModelParamMap.insert(std::make_pair("eWH78_HWW3", std::cref(eWH78_HWW3)));
459  ModelParamMap.insert(std::make_pair("eWH78_HWud", std::cref(eWH78_HWud)));
460  ModelParamMap.insert(std::make_pair("eWH78_Wud", std::cref(eWH78_Wud)));
461  ModelParamMap.insert(std::make_pair("eZH2_HZZ1", std::cref(eZH2_HZZ1)));
462  ModelParamMap.insert(std::make_pair("eZH2_HZZ2", std::cref(eZH2_HZZ2)));
463  ModelParamMap.insert(std::make_pair("eZH2_HZZ3", std::cref(eZH2_HZZ3)));
464  ModelParamMap.insert(std::make_pair("eZH2_HZA1", std::cref(eZH2_HZA1)));
465  ModelParamMap.insert(std::make_pair("eZH2_HZA2", std::cref(eZH2_HZA2)));
466  ModelParamMap.insert(std::make_pair("eZH2_HZuL", std::cref(eZH2_HZuL)));
467  ModelParamMap.insert(std::make_pair("eZH2_HZuR", std::cref(eZH2_HZuR)));
468  ModelParamMap.insert(std::make_pair("eZH2_HZdL", std::cref(eZH2_HZdL)));
469  ModelParamMap.insert(std::make_pair("eZH2_HZdR", std::cref(eZH2_HZdR)));
470  ModelParamMap.insert(std::make_pair("eZH2_ZuL", std::cref(eZH2_ZuL)));
471  ModelParamMap.insert(std::make_pair("eZH2_ZuR", std::cref(eZH2_ZuR)));
472  ModelParamMap.insert(std::make_pair("eZH2_ZdL", std::cref(eZH2_ZdL)));
473  ModelParamMap.insert(std::make_pair("eZH2_ZdR", std::cref(eZH2_ZdR)));
474  ModelParamMap.insert(std::make_pair("eZH78_HZZ1", std::cref(eZH78_HZZ1)));
475  ModelParamMap.insert(std::make_pair("eZH78_HZZ2", std::cref(eZH78_HZZ2)));
476  ModelParamMap.insert(std::make_pair("eZH78_HZZ3", std::cref(eZH78_HZZ3)));
477  ModelParamMap.insert(std::make_pair("eZH78_HZA1", std::cref(eZH78_HZA1)));
478  ModelParamMap.insert(std::make_pair("eZH78_HZA2", std::cref(eZH78_HZA2)));
479  ModelParamMap.insert(std::make_pair("eZH78_HZuL", std::cref(eZH78_HZuL)));
480  ModelParamMap.insert(std::make_pair("eZH78_HZuR", std::cref(eZH78_HZuR)));
481  ModelParamMap.insert(std::make_pair("eZH78_HZdL", std::cref(eZH78_HZdL)));
482  ModelParamMap.insert(std::make_pair("eZH78_HZdR", std::cref(eZH78_HZdR)));
483  ModelParamMap.insert(std::make_pair("eZH78_ZuL", std::cref(eZH78_ZuL)));
484  ModelParamMap.insert(std::make_pair("eZH78_ZuR", std::cref(eZH78_ZuR)));
485  ModelParamMap.insert(std::make_pair("eZH78_ZdL", std::cref(eZH78_ZdL)));
486  ModelParamMap.insert(std::make_pair("eZH78_ZdR", std::cref(eZH78_ZdR)));
487  ModelParamMap.insert(std::make_pair("ettH2_Htt", std::cref(ettH2_Htt)));
488  ModelParamMap.insert(std::make_pair("ettH2_Hgg", std::cref(ettH2_Hgg)));
489  ModelParamMap.insert(std::make_pair("ettH78_Htt", std::cref(ettH78_Htt)));
490  ModelParamMap.insert(std::make_pair("ettH78_Hgg", std::cref(ettH78_Hgg)));
491  if (FlagMwInput) {
492  ModelParamMap.insert(std::make_pair("MwInput", std::cref(MwInput)));
493  }
494 
495  if (FlagLeptonUniversal) {
496 // bsll/sbll entries only interesting (for the moment) if non-lepton universal. Set to 0 otherwise
497  CLQ1_1123 = 0.0;
498  CLQ1_2223 = 0.0;
499  CLQ1_3323 = 0.0;
500  CLQ1_1132 = 0.0;
501  CLQ1_2232 = 0.0;
502  CLQ1_3332 = 0.0;
503 
504  CLQ3_1123 = 0.0;
505  CLQ3_2223 = 0.0;
506  CLQ3_3323 = 0.0;
507  CLQ3_1132 = 0.0;
508  CLQ3_2232 = 0.0;
509  CLQ3_3332 = 0.0;
510 
511  Ced_1123 = 0.0;
512  Ced_2223 = 0.0;
513  Ced_3323 = 0.0;
514  Ced_1132 = 0.0;
515  Ced_2232 = 0.0;
516  Ced_3332 = 0.0;
517 
518  CLd_1123 = 0.0;
519  CLd_2223 = 0.0;
520  CLd_3323 = 0.0;
521  CLd_1132 = 0.0;
522  CLd_2232 = 0.0;
523  CLd_3332 = 0.0;
524 
525  CQe_2311 = 0.0;
526  CQe_2322 = 0.0;
527  CQe_2333 = 0.0;
528  CQe_3211 = 0.0;
529  CQe_3222 = 0.0;
530  CQe_3233 = 0.0;
531  }
532 
533 }
534 
536 {
537  if (!NPbase::PostUpdate()) return (false);
538 
540  v2_over_LambdaNP2 = v() * v() / LambdaNP2;
541  if (FlagMwInput)
542  cW_tree = MwInput / Mz;
543  else
544  cW_tree = Mw_tree() / Mz;
546  sW2_tree = 1.0 - cW2_tree;
547  sW_tree = sqrt(sW2_tree);
548 
551  delta_AZ = 2.0 * sW_tree * cW_tree * (CHW - CHB) * v2_over_LambdaNP2
553  delta_h = (-CHD / 4.0 + CHbox) * v2_over_LambdaNP2;
554 
555  return (true);
556 }
557 
558 void NPEffectiveGIMR::setParameter(const std::string name, const double& value)
559 {
560  if (name.compare("CG") == 0)
561  CG = value;
562  else if (name.compare("CW") == 0)
563  CW = value;
564  else if (name.compare("CHG") == 0)
565  CHG = value;
566  else if (name.compare("CHW") == 0)
567  CHW = value;
568  else if (name.compare("CHB") == 0)
569  CHB = value;
570  else if (name.compare("CHWB") == 0)
571  CHWB = value;
572  else if (name.compare("CHD") == 0)
573  CHD = value;
574  else if (name.compare("CHbox") == 0)
575  CHbox = value;
576  else if (name.compare("CH") == 0)
577  CH = value;
578  else if (name.compare("CHL1_11") == 0)
579  CHL1_11 = value;
580  else if (name.compare("CHL1_12r") == 0)
581  CHL1_12r = value;
582  else if (name.compare("CHL1_13r") == 0)
583  CHL1_13r = value;
584  else if (name.compare("CHL1_22") == 0)
585  CHL1_22 = value;
586  else if (name.compare("CHL1_23r") == 0)
587  CHL1_23r = value;
588  else if (name.compare("CHL1_33") == 0)
589  CHL1_33 = value;
590  else if (name.compare("CHL1_12i") == 0)
591  CHL1_12i = value;
592  else if (name.compare("CHL1_13i") == 0)
593  CHL1_13i = value;
594  else if (name.compare("CHL1_23i") == 0)
595  CHL1_23i = value;
596  else if (name.compare("CHL1") == 0) {
597  CHL1_11 = value;
598  CHL1_12r = 0.0;
599  CHL1_13r = 0.0;
600  CHL1_22 = value;
601  CHL1_23r = 0.0;
602  CHL1_33 = value;
603  CHL1_12i = 0.0;
604  CHL1_13i = 0.0;
605  CHL1_23i = 0.0;
606  } else if (name.compare("CHL3_11") == 0)
607  CHL3_11 = value;
608  else if (name.compare("CHL3_12r") == 0)
609  CHL3_12r = value;
610  else if (name.compare("CHL3_13r") == 0)
611  CHL3_13r = value;
612  else if (name.compare("CHL3_22") == 0)
613  CHL3_22 = value;
614  else if (name.compare("CHL3_23r") == 0)
615  CHL3_23r = value;
616  else if (name.compare("CHL3_33") == 0)
617  CHL3_33 = value;
618  else if (name.compare("CHL3_12i") == 0)
619  CHL3_12i = value;
620  else if (name.compare("CHL3_13i") == 0)
621  CHL3_13i = value;
622  else if (name.compare("CHL3_23i") == 0)
623  CHL3_23i = value;
624  else if (name.compare("CHL3") == 0) {
625  CHL3_11 = value;
626  CHL3_12r = 0.0;
627  CHL3_13r = 0.0;
628  CHL3_22 = value;
629  CHL3_23r = 0.0;
630  CHL3_33 = value;
631  CHL3_12i = 0.0;
632  CHL3_13i = 0.0;
633  CHL3_23i = 0.0;
634  } else if (name.compare("CHe_11") == 0)
635  CHe_11 = value;
636  else if (name.compare("CHe_12r") == 0)
637  CHe_12r = value;
638  else if (name.compare("CHe_13r") == 0)
639  CHe_13r = value;
640  else if (name.compare("CHe_22") == 0)
641  CHe_22 = value;
642  else if (name.compare("CHe_23r") == 0)
643  CHe_23r = value;
644  else if (name.compare("CHe_33") == 0)
645  CHe_33 = value;
646  else if (name.compare("CHe_12i") == 0)
647  CHe_12i = value;
648  else if (name.compare("CHe_13i") == 0)
649  CHe_13i = value;
650  else if (name.compare("CHe_23i") == 0)
651  CHe_23i = value;
652  else if (name.compare("CHe") == 0) {
653  CHe_11 = value;
654  CHe_12r = 0.0;
655  CHe_13r = 0.0;
656  CHe_22 = value;
657  CHe_23r = 0.0;
658  CHe_33 = value;
659  CHe_12i = 0.0;
660  CHe_13i = 0.0;
661  CHe_23i = 0.0;
662  } else if (name.compare("CHQ1_11") == 0)
663  CHQ1_11 = value;
664  else if (name.compare("CHQ1_12r") == 0)
665  CHQ1_12r = value;
666  else if (name.compare("CHQ1_13r") == 0)
667  CHQ1_13r = value;
668  else if (name.compare("CHQ1_22") == 0)
669  CHQ1_22 = value;
670  else if (name.compare("CHQ1_23r") == 0)
671  CHQ1_23r = value;
672  else if (name.compare("CHQ1_33") == 0)
673  CHQ1_33 = value;
674  else if (name.compare("CHQ1_12i") == 0)
675  CHQ1_12i = value;
676  else if (name.compare("CHQ1_13i") == 0)
677  CHQ1_13i = value;
678  else if (name.compare("CHQ1_23i") == 0)
679  CHQ1_23i = value;
680  else if (name.compare("CHQ1") == 0) {
681  CHQ1_11 = value;
682  CHQ1_12r = 0.0;
683  CHQ1_13r = 0.0;
684  CHQ1_22 = value;
685  CHQ1_23r = 0.0;
686  CHQ1_33 = value;
687  CHQ1_12i = 0.0;
688  CHQ1_13i = 0.0;
689  CHQ1_23i = 0.0;
690  } else if (name.compare("CHQ3_11") == 0)
691  CHQ3_11 = value;
692  else if (name.compare("CHQ3_12r") == 0)
693  CHQ3_12r = value;
694  else if (name.compare("CHQ3_13r") == 0)
695  CHQ3_13r = value;
696  else if (name.compare("CHQ3_22") == 0)
697  CHQ3_22 = value;
698  else if (name.compare("CHQ3_23r") == 0)
699  CHQ3_23r = value;
700  else if (name.compare("CHQ3_33") == 0)
701  CHQ3_33 = value;
702  else if (name.compare("CHQ3_12i") == 0)
703  CHQ3_12i = value;
704  else if (name.compare("CHQ3_13i") == 0)
705  CHQ3_13i = value;
706  else if (name.compare("CHQ3_23i") == 0)
707  CHQ3_23i = value;
708  else if (name.compare("CHQ3") == 0) {
709  CHQ3_11 = value;
710  CHQ3_12r = 0.0;
711  CHQ3_13r = 0.0;
712  CHQ3_22 = value;
713  CHQ3_23r = 0.0;
714  CHQ3_33 = value;
715  CHQ3_12i = 0.0;
716  CHQ3_13i = 0.0;
717  CHQ3_23i = 0.0;
718  } else if (name.compare("CHu_11") == 0)
719  CHu_11 = value;
720  else if (name.compare("CHu_12r") == 0)
721  CHu_12r = value;
722  else if (name.compare("CHu_13r") == 0)
723  CHu_13r = value;
724  else if (name.compare("CHu_22") == 0)
725  CHu_22 = value;
726  else if (name.compare("CHu_23r") == 0)
727  CHu_23r = value;
728  else if (name.compare("CHu_33") == 0)
729  CHu_33 = value;
730  else if (name.compare("CHu_12i") == 0)
731  CHu_12i = value;
732  else if (name.compare("CHu_13i") == 0)
733  CHu_13i = value;
734  else if (name.compare("CHu_23i") == 0)
735  CHu_23i = value;
736  else if (name.compare("CHu") == 0) {
737  CHu_11 = value;
738  CHu_12r = 0.0;
739  CHu_13r = 0.0;
740  CHu_22 = value;
741  CHu_23r = 0.0;
742  CHu_33 = value;
743  CHu_12i = 0.0;
744  CHu_13i = 0.0;
745  CHu_23i = 0.0;
746  } else if (name.compare("CHd_11") == 0)
747  CHd_11 = value;
748  else if (name.compare("CHd_12r") == 0)
749  CHd_12r = value;
750  else if (name.compare("CHd_13r") == 0)
751  CHd_13r = value;
752  else if (name.compare("CHd_22") == 0)
753  CHd_22 = value;
754  else if (name.compare("CHd_23r") == 0)
755  CHd_23r = value;
756  else if (name.compare("CHd_33") == 0)
757  CHd_33 = value;
758  else if (name.compare("CHd_12i") == 0)
759  CHd_12i = value;
760  else if (name.compare("CHd_13i") == 0)
761  CHd_13i = value;
762  else if (name.compare("CHd_23i") == 0)
763  CHd_23i = value;
764  else if (name.compare("CHd") == 0) {
765  CHd_11 = value;
766  CHd_12r = 0.0;
767  CHd_13r = 0.0;
768  CHd_22 = value;
769  CHd_23r = 0.0;
770  CHd_33 = value;
771  CHd_12i = 0.0;
772  CHd_13i = 0.0;
773  CHd_23i = 0.0;
774  } else if (name.compare("CHud_11r") == 0)
775  CHud_11r = value;
776  else if (name.compare("CHud_12r") == 0)
777  CHud_12r = value;
778  else if (name.compare("CHud_13r") == 0)
779  CHud_13r = value;
780  else if (name.compare("CHud_22r") == 0)
781  CHud_22r = value;
782  else if (name.compare("CHud_23r") == 0)
783  CHud_23r = value;
784  else if (name.compare("CHud_33r") == 0)
785  CHud_33r = value;
786  else if (name.compare("CHud_r") == 0) {
787  CHud_11r = value;
788  CHud_12r = 0.0;
789  CHud_13r = 0.0;
790  CHud_22r = value;
791  CHud_23r = 0.0;
792  CHud_33r = value;
793  } else if (name.compare("CHud_11i") == 0)
794  CHud_11i = value;
795  else if (name.compare("CHud_12i") == 0)
796  CHud_12i = value;
797  else if (name.compare("CHud_13i") == 0)
798  CHud_13i = value;
799  else if (name.compare("CHud_22i") == 0)
800  CHud_22i = value;
801  else if (name.compare("CHud_23i") == 0)
802  CHud_23i = value;
803  else if (name.compare("CHud_33i") == 0)
804  CHud_33i = value;
805  else if (name.compare("CHud_i") == 0) {
806  CHud_11i = value;
807  CHud_12i = 0.0;
808  CHud_13i = 0.0;
809  CHud_22i = value;
810  CHud_23i = 0.0;
811  CHud_33i = value;
812  } else if (name.compare("CeH_11r") == 0)
813  CeH_11r = value;
814  else if (name.compare("CeH_12r") == 0)
815  CeH_12r = value;
816  else if (name.compare("CeH_13r") == 0)
817  CeH_13r = value;
818  else if (name.compare("CeH_22r") == 0)
819  CeH_22r = value;
820  else if (name.compare("CeH_23r") == 0)
821  CeH_23r = value;
822  else if (name.compare("CeH_33r") == 0)
823  CeH_33r = value;
824  else if (name.compare("CeH_r") == 0) {
825  CeH_11r = value;
826  CeH_12r = 0.0;
827  CeH_13r = 0.0;
828  CeH_22r = value;
829  CeH_23r = 0.0;
830  CeH_33r = value;
831  } else if (name.compare("CeH_11i") == 0)
832  CeH_11i = value;
833  else if (name.compare("CeH_12i") == 0)
834  CeH_12i = value;
835  else if (name.compare("CeH_13i") == 0)
836  CeH_13i = value;
837  else if (name.compare("CeH_22i") == 0)
838  CeH_22i = value;
839  else if (name.compare("CeH_23i") == 0)
840  CeH_23i = value;
841  else if (name.compare("CeH_33i") == 0)
842  CeH_33i = value;
843  else if (name.compare("CeH_i") == 0) {
844  CeH_11i = value;
845  CeH_12i = 0.0;
846  CeH_13i = 0.0;
847  CeH_22i = value;
848  CeH_23i = 0.0;
849  CeH_33i = value;
850  } else if (name.compare("CuH_11r") == 0)
851  CuH_11r = value;
852  else if (name.compare("CuH_12r") == 0)
853  CuH_12r = value;
854  else if (name.compare("CuH_13r") == 0)
855  CuH_13r = value;
856  else if (name.compare("CuH_22r") == 0)
857  CuH_22r = value;
858  else if (name.compare("CuH_23r") == 0)
859  CuH_23r = value;
860  else if (name.compare("CuH_33r") == 0)
861  CuH_33r = value;
862  else if (name.compare("CuH_r") == 0) {
863  CuH_11r = value;
864  CuH_12r = 0.0;
865  CuH_13r = 0.0;
866  CuH_22r = value;
867  CuH_23r = 0.0;
868  CuH_33r = value;
869  } else if (name.compare("CuH_11i") == 0)
870  CuH_11i = value;
871  else if (name.compare("CuH_12i") == 0)
872  CuH_12i = value;
873  else if (name.compare("CuH_13i") == 0)
874  CuH_13i = value;
875  else if (name.compare("CuH_22i") == 0)
876  CuH_22i = value;
877  else if (name.compare("CuH_23i") == 0)
878  CuH_23i = value;
879  else if (name.compare("CuH_33i") == 0)
880  CuH_33i = value;
881  else if (name.compare("CuH_i") == 0) {
882  CuH_11i = value;
883  CuH_12i = 0.0;
884  CuH_13i = 0.0;
885  CuH_22i = value;
886  CuH_23i = 0.0;
887  CuH_33i = value;
888  } else if (name.compare("CdH_11r") == 0)
889  CdH_11r = value;
890  else if (name.compare("CdH_12r") == 0)
891  CdH_12r = value;
892  else if (name.compare("CdH_13r") == 0)
893  CdH_13r = value;
894  else if (name.compare("CdH_22r") == 0)
895  CdH_22r = value;
896  else if (name.compare("CdH_23r") == 0)
897  CdH_23r = value;
898  else if (name.compare("CdH_33r") == 0)
899  CdH_33r = value;
900  else if (name.compare("CdH_r") == 0) {
901  CdH_11r = value;
902  CdH_12r = 0.0;
903  CdH_13r = 0.0;
904  CdH_22r = value;
905  CdH_23r = 0.0;
906  CdH_33r = value;
907  } else if (name.compare("CdH_11i") == 0)
908  CdH_11i = value;
909  else if (name.compare("CdH_12i") == 0)
910  CdH_12i = value;
911  else if (name.compare("CdH_13i") == 0)
912  CdH_13i = value;
913  else if (name.compare("CdH_22i") == 0)
914  CdH_22i = value;
915  else if (name.compare("CdH_23i") == 0)
916  CdH_23i = value;
917  else if (name.compare("CdH_33i") == 0)
918  CdH_33i = value;
919  else if (name.compare("CdH_i") == 0) {
920  CdH_11i = value;
921  CdH_12i = 0.0;
922  CdH_13i = 0.0;
923  CdH_22i = value;
924  CdH_23i = 0.0;
925  CdH_33i = value;
926  } else if (name.compare("CuG_11r") == 0)
927  CuG_11r = value;
928  else if (name.compare("CuG_12r") == 0)
929  CuG_12r = value;
930  else if (name.compare("CuG_13r") == 0)
931  CuG_13r = value;
932  else if (name.compare("CuG_22r") == 0)
933  CuG_22r = value;
934  else if (name.compare("CuG_23r") == 0)
935  CuG_23r = value;
936  else if (name.compare("CuG_33r") == 0)
937  CuG_33r = value;
938  else if (name.compare("CuG_r") == 0) {
939  CuG_11r = value;
940  CuG_12r = 0.0;
941  CuG_13r = 0.0;
942  CuG_22r = value;
943  CuG_23r = 0.0;
944  CuG_33r = value;
945  } else if (name.compare("CuG_11i") == 0)
946  CuG_11i = value;
947  else if (name.compare("CuG_12i") == 0)
948  CuG_12i = value;
949  else if (name.compare("CuG_13i") == 0)
950  CuG_13i = value;
951  else if (name.compare("CuG_22i") == 0)
952  CuG_22i = value;
953  else if (name.compare("CuG_23i") == 0)
954  CuG_23i = value;
955  else if (name.compare("CuG_33i") == 0)
956  CuG_33i = value;
957  else if (name.compare("CuG_i") == 0) {
958  CuG_11i = value;
959  CuG_12i = 0.0;
960  CuG_13i = 0.0;
961  CuG_22i = value;
962  CuG_23i = 0.0;
963  CuG_33i = value;
964  } else if (name.compare("CuW_11r") == 0)
965  CuW_11r = value;
966  else if (name.compare("CuW_12r") == 0)
967  CuW_12r = value;
968  else if (name.compare("CuW_13r") == 0)
969  CuW_13r = value;
970  else if (name.compare("CuW_22r") == 0)
971  CuW_22r = value;
972  else if (name.compare("CuW_23r") == 0)
973  CuW_23r = value;
974  else if (name.compare("CuW_33r") == 0)
975  CuW_33r = value;
976  else if (name.compare("CuW_r") == 0) {
977  CuW_11r = value;
978  CuW_12r = 0.0;
979  CuW_13r = 0.0;
980  CuW_22r = value;
981  CuW_23r = 0.0;
982  CuW_33r = value;
983  } else if (name.compare("CuW_11i") == 0)
984  CuW_11i = value;
985  else if (name.compare("CuW_12i") == 0)
986  CuW_12i = value;
987  else if (name.compare("CuW_13i") == 0)
988  CuW_13i = value;
989  else if (name.compare("CuW_22i") == 0)
990  CuW_22i = value;
991  else if (name.compare("CuW_23i") == 0)
992  CuW_23i = value;
993  else if (name.compare("CuW_33i") == 0)
994  CuW_33i = value;
995  else if (name.compare("CuW_i") == 0) {
996  CuW_11i = value;
997  CuW_12i = 0.0;
998  CuW_13i = 0.0;
999  CuW_22i = value;
1000  CuW_23i = 0.0;
1001  CuW_33i = value;
1002  } else if (name.compare("CuB_11r") == 0)
1003  CuB_11r = value;
1004  else if (name.compare("CuB_12r") == 0)
1005  CuB_12r = value;
1006  else if (name.compare("CuB_13r") == 0)
1007  CuB_13r = value;
1008  else if (name.compare("CuB_22r") == 0)
1009  CuB_22r = value;
1010  else if (name.compare("CuB_23r") == 0)
1011  CuB_23r = value;
1012  else if (name.compare("CuB_33r") == 0)
1013  CuB_33r = value;
1014  else if (name.compare("CuB_r") == 0) {
1015  CuB_11r = value;
1016  CuB_12r = 0.0;
1017  CuB_13r = 0.0;
1018  CuB_22r = value;
1019  CuB_23r = 0.0;
1020  CuB_33r = value;
1021  } else if (name.compare("CuB_11i") == 0)
1022  CuB_11i = value;
1023  else if (name.compare("CuB_12i") == 0)
1024  CuB_12i = value;
1025  else if (name.compare("CuB_13i") == 0)
1026  CuB_13i = value;
1027  else if (name.compare("CuB_22i") == 0)
1028  CuB_22i = value;
1029  else if (name.compare("CuB_23i") == 0)
1030  CuB_23i = value;
1031  else if (name.compare("CuB_33i") == 0)
1032  CuB_33i = value;
1033  else if (name.compare("CuB_i") == 0) {
1034  CuB_11i = value;
1035  CuB_12i = 0.0;
1036  CuB_13i = 0.0;
1037  CuB_22i = value;
1038  CuB_23i = 0.0;
1039  CuB_33i = value;
1040  } else if (name.compare("CLL_1111") == 0) {
1041  CLL_1111 = value;
1042  } else if (name.compare("CLL_1221") == 0) {
1043  CLL_1221 = value;
1044  CLL_2112 = value;
1045  CLL_2211 = value;
1046  CLL_1122 = value;
1047  } else if (name.compare("CLL_1133") == 0) {
1048  CLL_1331 = value;
1049  CLL_3113 = value;
1050  CLL_3311 = value;
1051  CLL_1133 = value;
1052  } else if (name.compare("CLL") == 0) {
1053  CLL_1111 = value;
1054  CLL_1221 = value;
1055  CLL_2112 = value;
1056  CLL_2211 = value;
1057  CLL_1122 = value;
1058  CLL_3311 = value;
1059  CLL_1133 = value;
1060  CLL_1331 = value;
1061  CLL_3113 = value;
1062  } else if (name.compare("CLQ1_1111") == 0) {
1063  CLQ1_1111 = value;
1064  } else if (name.compare("CLQ1_1122") == 0) {
1065  CLQ1_1122 = value;
1066  CLQ1_2211 = value;
1067  CLQ1_1221 = value;
1068  CLQ1_2112 = value;
1069  } else if (name.compare("CLQ1_1133") == 0) {
1070  CLQ1_1133 = value;
1071  CLQ1_3311 = value;
1072  CLQ1_1331 = value;
1073  CLQ1_3113 = value;
1074  } else if (name.compare("CLQ1") == 0) {
1075  CLQ1_1111 = value;
1076  CLQ1_1122 = value;
1077  CLQ1_2211 = value;
1078  CLQ1_1221 = value;
1079  CLQ1_2112 = value;
1080  CLQ1_1133 = value;
1081  CLQ1_3311 = value;
1082  CLQ1_1331 = value;
1083  CLQ1_3113 = value;
1084  } else if (name.compare("CLQ1_1111") == 0) {
1085  CLQ1_1111 = value;
1086  } else if (name.compare("CLQ1_1122") == 0) {
1087  CLQ1_1122 = value;
1088  CLQ1_2211 = value;
1089  CLQ1_2112 = value;
1090  CLQ1_1221 = value;
1091  } else if (name.compare("CLQ1_1133") == 0) {
1092  CLQ1_1133 = value;
1093  CLQ1_3311 = value;
1094  CLQ1_3113 = value;
1095  CLQ1_1331 = value;
1096  } else if (name.compare("CLQ1_1123") == 0) {
1097  CLQ1_1123 = value;
1098  } else if (name.compare("CLQ1_2223") == 0) {
1099  CLQ1_2223 = value;
1100  } else if (name.compare("CLQ1_3323") == 0) {
1101  CLQ1_3323 = value;
1102  } else if (name.compare("CLQ1_1132") == 0) {
1103  CLQ1_1132 = value;
1104  } else if (name.compare("CLQ1_2232") == 0) {
1105  CLQ1_2232 = value;
1106  } else if (name.compare("CLQ1_3332") == 0) {
1107  CLQ1_3332 = value;
1108  } else if (name.compare("CLQ3") == 0) {
1109  CLQ3_1111 = value;
1110  CLQ3_1122 = value;
1111  CLQ3_2211 = value;
1112  CLQ3_1221 = value;
1113  CLQ3_2112 = value;
1114  CLQ3_1133 = value;
1115  CLQ3_3311 = value;
1116  CLQ3_1331 = value;
1117  CLQ3_3113 = value;
1118  } else if (name.compare("CLQ3_1111") == 0) {
1119  CLQ3_1111 = value;
1120  } else if (name.compare("CLQ3_1122") == 0) {
1121  CLQ3_1122 = value;
1122  CLQ3_2211 = value;
1123  CLQ3_2112 = value;
1124  CLQ3_1221 = value;
1125  } else if (name.compare("CLQ3_1133") == 0) {
1126  CLQ3_1133 = value;
1127  CLQ3_3311 = value;
1128  CLQ3_3113 = value;
1129  CLQ3_1331 = value;
1130  } else if (name.compare("CLQ3_1123") == 0) {
1131  CLQ3_1123 = value;
1132  } else if (name.compare("CLQ3_2223") == 0) {
1133  CLQ3_2223 = value;
1134  } else if (name.compare("CLQ3_3323") == 0) {
1135  CLQ3_3323 = value;
1136  } else if (name.compare("CLQ3_1132") == 0) {
1137  CLQ3_1132 = value;
1138  } else if (name.compare("CLQ3_2232") == 0) {
1139  CLQ3_2232 = value;
1140  } else if (name.compare("CLQ3_3332") == 0) {
1141  CLQ3_3332 = value;
1142  } else if (name.compare("Cee") == 0) {
1143  Cee_1111 = value;
1144  Cee_1122 = value;
1145  Cee_2211 = value;
1146  Cee_1133 = value;
1147  Cee_3311 = value;
1148  } else if (name.compare("Cee_1111") == 0) {
1149  Cee_1111 = value;
1150  } else if (name.compare("Cee_1122") == 0) {
1151  Cee_1122 = value;
1152  Cee_2211 = value;
1153  } else if (name.compare("Cee_1133") == 0) {
1154  Cee_1133 = value;
1155  Cee_3311 = value;
1156  } else if (name.compare("Ceu") == 0) {
1157  Ceu_1111 = value;
1158  Ceu_1122 = value;
1159  Ceu_2211 = value;
1160  Ceu_1133 = value;
1161  Ceu_3311 = value;
1162  } else if (name.compare("Ceu_1111") == 0) {
1163  Ceu_1111 = value;
1164  } else if (name.compare("Ceu_1122") == 0) {
1165  Ceu_1122 = value;
1166  Ceu_2211 = value;
1167  } else if (name.compare("Ceu_1133") == 0) {
1168  Ceu_1133 = value;
1169  Ceu_3311 = value;
1170  } else if (name.compare("Ced") == 0) {
1171  Ced_1111 = value;
1172  Ced_1122 = value;
1173  Ced_2211 = value;
1174  Ced_1133 = value;
1175  Ced_3311 = value;
1176  } else if (name.compare("Ced_1111") == 0) {
1177  Ced_1111 = value;
1178  } else if (name.compare("Ced_1122") == 0) {
1179  Ced_1122 = value;
1180  Ced_2211 = value;
1181  } else if (name.compare("Ced_1133") == 0) {
1182  Ced_1133 = value;
1183  Ced_3311 = value;
1184  } else if (name.compare("Ced_1123") == 0) {
1185  Ced_1123 = value;
1186  } else if (name.compare("Ced_2223") == 0) {
1187  Ced_2223 = value;
1188  } else if (name.compare("Ced_3323") == 0) {
1189  Ced_3323 = value;
1190  } else if (name.compare("Ced_1132") == 0) {
1191  Ced_1132 = value;
1192  } else if (name.compare("Ced_2232") == 0) {
1193  Ced_2232 = value;
1194  } else if (name.compare("Ced_3332") == 0) {
1195  Ced_3332 = value;
1196  } else if (name.compare("CLe") == 0) {
1197  CLe_1111 = value;
1198  CLe_1122 = value;
1199  CLe_2211 = value;
1200  CLe_1133 = value;
1201  CLe_3311 = value;
1202  } else if (name.compare("CLe_1111") == 0) {
1203  CLe_1111 = value;
1204  } else if (name.compare("CLe_1122") == 0) {
1205  CLe_1122 = value;
1206  CLe_2211 = value;
1207  } else if (name.compare("CLe_1133") == 0) {
1208  CLe_1133 = value;
1209  CLe_3311 = value;
1210  } else if (name.compare("CLu") == 0) {
1211  CLu_1111 = value;
1212  CLu_1122 = value;
1213  CLu_2211 = value;
1214  CLu_1133 = value;
1215  CLu_3311 = value;
1216  } else if (name.compare("CLu_1111") == 0) {
1217  CLu_1111 = value;
1218  } else if (name.compare("CLu_1122") == 0) {
1219  CLu_1122 = value;
1220  CLu_2211 = value;
1221  } else if (name.compare("CLu_1133") == 0) {
1222  CLu_1133 = value;
1223  CLu_3311 = value;
1224  } else if (name.compare("CLd") == 0) {
1225  CLd_1111 = value;
1226  CLd_1122 = value;
1227  CLd_2211 = value;
1228  CLd_1133 = value;
1229  CLd_3311 = value;
1230  } else if (name.compare("CLd_1111") == 0) {
1231  CLd_1111 = value;
1232  } else if (name.compare("CLd_1122") == 0) {
1233  CLd_1122 = value;
1234  CLd_2211 = value;
1235  } else if (name.compare("CLd_1133") == 0) {
1236  CLd_1133 = value;
1237  CLd_3311 = value;
1238  } else if (name.compare("CLd_1123") == 0) {
1239  CLd_1123 = value;
1240  } else if (name.compare("CLd_2223") == 0) {
1241  CLd_2223 = value;
1242  } else if (name.compare("CLd_3323") == 0) {
1243  CLd_3323 = value;
1244  } else if (name.compare("CLd_1132") == 0) {
1245  CLd_1132 = value;
1246  } else if (name.compare("CLd_2232") == 0) {
1247  CLd_2232 = value;
1248  } else if (name.compare("CLd_3332") == 0) {
1249  CLd_3332 = value;
1250  } else if (name.compare("CQe") == 0) {
1251  CQe_1111 = value;
1252  CQe_1122 = value;
1253  CQe_2211 = value;
1254  CQe_1133 = value;
1255  CQe_3311 = value;
1256  } else if (name.compare("CQe_1111") == 0) {
1257  CQe_1111 = value;
1258  } else if (name.compare("CQe_1122") == 0) {
1259  CQe_1122 = value;
1260  CQe_2211 = value;
1261  } else if (name.compare("CQe_1133") == 0) {
1262  CQe_1133 = value;
1263  CQe_3311 = value;
1264  } else if (name.compare("CQe_2311") == 0) {
1265  CQe_2311 = value;
1266  } else if (name.compare("CQe_2322") == 0) {
1267  CQe_2322 = value;
1268  } else if (name.compare("CQe_2333") == 0) {
1269  CQe_2333 = value;
1270  } else if (name.compare("CQe_3211") == 0) {
1271  CQe_3211 = value;
1272  } else if (name.compare("CQe_3222") == 0) {
1273  CQe_3222 = value;
1274  } else if (name.compare("CQe_3233") == 0) {
1275  CQe_3233 = value;
1276  } else if (name.compare("Lambda_NP") == 0) {
1277  Lambda_NP = value;
1278  } else if (name.compare("eVBF2_HZZ1") == 0) {
1279  eVBF2_HZZ1 = value;
1280  } else if (name.compare("eVBF2_HZZ2") == 0) {
1281  eVBF2_HZZ2 = value;
1282  } else if (name.compare("eVBF2_HZZ3") == 0) {
1283  eVBF2_HZZ3 = value;
1284  } else if (name.compare("eVBF2_HZA1") == 0) {
1285  eVBF2_HZA1 = value;
1286  } else if (name.compare("eVBF2_HZA2") == 0) {
1287  eVBF2_HZA2 = value;
1288  } else if (name.compare("eVBF2_HAA") == 0) {
1289  eVBF2_HAA = value;
1290  } else if (name.compare("eVBF2_HWW1") == 0) {
1291  eVBF2_HWW1 = value;
1292  } else if (name.compare("eVBF2_HWW2") == 0) {
1293  eVBF2_HWW2 = value;
1294  } else if (name.compare("eVBF2_HWW3") == 0) {
1295  eVBF2_HWW3 = value;
1296  } else if (name.compare("eVBF2_Hgg") == 0) {
1297  eVBF2_Hgg = value;
1298  } else if (name.compare("eVBF2_HZuL") == 0) {
1299  eVBF2_HZuL = value;
1300  } else if (name.compare("eVBF2_HZuR") == 0) {
1301  eVBF2_HZuR = value;
1302  } else if (name.compare("eVBF2_HZdL") == 0) {
1303  eVBF2_HZdL = value;
1304  } else if (name.compare("eVBF2_HZdR") == 0) {
1305  eVBF2_HZdR = value;
1306  } else if (name.compare("eVBF2_HWud") == 0) {
1307  eVBF2_HWud = value;
1308  } else if (name.compare("eVBF2_ZuL") == 0) {
1309  eVBF2_ZuL = value;
1310  } else if (name.compare("eVBF2_ZuR") == 0) {
1311  eVBF2_ZuR = value;
1312  } else if (name.compare("eVBF2_ZdL") == 0) {
1313  eVBF2_ZdL = value;
1314  } else if (name.compare("eVBF2_ZdR") == 0) {
1315  eVBF2_ZdR = value;
1316  } else if (name.compare("eVBF2_Wud") == 0) {
1317  eVBF2_Wud = value;
1318  } else if (name.compare("eVBF78_HZZ1") == 0) {
1319  eVBF78_HZZ1 = value;
1320  } else if (name.compare("eVBF78_HZZ2") == 0) {
1321  eVBF78_HZZ2 = value;
1322  } else if (name.compare("eVBF78_HZZ3") == 0) {
1323  eVBF78_HZZ3 = value;
1324  } else if (name.compare("eVBF78_HZA1") == 0) {
1325  eVBF78_HZA1 = value;
1326  } else if (name.compare("eVBF78_HZA2") == 0) {
1327  eVBF78_HZA2 = value;
1328  } else if (name.compare("eVBF78_HAA") == 0) {
1329  eVBF78_HAA = value;
1330  } else if (name.compare("eVBF78_HWW1") == 0) {
1331  eVBF78_HWW1 = value;
1332  } else if (name.compare("eVBF78_HWW2") == 0) {
1333  eVBF78_HWW2 = value;
1334  } else if (name.compare("eVBF78_HWW3") == 0) {
1335  eVBF78_HWW3 = value;
1336  } else if (name.compare("eVBF78_Hgg") == 0) {
1337  eVBF78_Hgg = value;
1338  } else if (name.compare("eVBF78_HZuL") == 0) {
1339  eVBF78_HZuL = value;
1340  } else if (name.compare("eVBF78_HZuR") == 0) {
1341  eVBF78_HZuR = value;
1342  } else if (name.compare("eVBF78_HZdL") == 0) {
1343  eVBF78_HZdL = value;
1344  } else if (name.compare("eVBF78_HZdR") == 0) {
1345  eVBF78_HZdR = value;
1346  } else if (name.compare("eVBF78_HWud") == 0) {
1347  eVBF78_HWud = value;
1348  } else if (name.compare("eVBF78_ZuL") == 0) {
1349  eVBF78_ZuL = value;
1350  } else if (name.compare("eVBF78_ZuR") == 0) {
1351  eVBF78_ZuR = value;
1352  } else if (name.compare("eVBF78_ZdL") == 0) {
1353  eVBF78_ZdL = value;
1354  } else if (name.compare("eVBF78_ZdR") == 0) {
1355  eVBF78_ZdR = value;
1356  } else if (name.compare("eVBF78_Wud") == 0) {
1357  eVBF78_Wud = value;
1358  } else if (name.compare("eWH2_HWW1") == 0) {
1359  eWH2_HWW1 = value;
1360  } else if (name.compare("eWH2_HWW2") == 0) {
1361  eWH2_HWW2 = value;
1362  } else if (name.compare("eWH2_HWW3") == 0) {
1363  eWH2_HWW3 = value;
1364  } else if (name.compare("eWH2_HWud") == 0) {
1365  eWH2_HWud = value;
1366  } else if (name.compare("eWH2_Wud") == 0) {
1367  eWH2_Wud = value;
1368  } else if (name.compare("eWH78_HWW1") == 0) {
1369  eWH78_HWW1 = value;
1370  } else if (name.compare("eWH78_HWW2") == 0) {
1371  eWH78_HWW2 = value;
1372  } else if (name.compare("eWH78_HWW3") == 0) {
1373  eWH78_HWW3 = value;
1374  } else if (name.compare("eWH78_HWud") == 0) {
1375  eWH78_HWud = value;
1376  } else if (name.compare("eWH78_Wud") == 0) {
1377  eWH78_Wud = value;
1378  } else if (name.compare("eZH2_HZZ1") == 0) {
1379  eZH2_HZZ1 = value;
1380  } else if (name.compare("eZH2_HZZ2") == 0) {
1381  eZH2_HZZ2 = value;
1382  } else if (name.compare("eZH2_HZZ3") == 0) {
1383  eZH2_HZZ3 = value;
1384  } else if (name.compare("eZH2_HZA1") == 0) {
1385  eZH2_HZA1 = value;
1386  } else if (name.compare("eZH2_HZA2") == 0) {
1387  eZH2_HZA2 = value;
1388  } else if (name.compare("eZH2_HZuL") == 0) {
1389  eZH2_HZuL = value;
1390  } else if (name.compare("eZH2_HZuR") == 0) {
1391  eZH2_HZuR = value;
1392  } else if (name.compare("eZH2_HZdL") == 0) {
1393  eZH2_HZdL = value;
1394  } else if (name.compare("eZH2_HZdR") == 0) {
1395  eZH2_HZdR = value;
1396  } else if (name.compare("eZH2_ZuL") == 0) {
1397  eZH2_ZuL = value;
1398  } else if (name.compare("eZH2_ZuR") == 0) {
1399  eZH2_ZuR = value;
1400  } else if (name.compare("eZH2_ZdL") == 0) {
1401  eZH2_ZdL = value;
1402  } else if (name.compare("eZH2_ZdR") == 0) {
1403  eZH2_ZdR = value;
1404  } else if (name.compare("eZH78_HZZ1") == 0) {
1405  eZH78_HZZ1 = value;
1406  } else if (name.compare("eZH78_HZZ2") == 0) {
1407  eZH78_HZZ2 = value;
1408  } else if (name.compare("eZH78_HZZ3") == 0) {
1409  eZH78_HZZ3 = value;
1410  } else if (name.compare("eZH78_HZA1") == 0) {
1411  eZH78_HZA1 = value;
1412  } else if (name.compare("eZH78_HZA2") == 0) {
1413  eZH78_HZA2 = value;
1414  } else if (name.compare("eZH78_HZuL") == 0) {
1415  eZH78_HZuL = value;
1416  } else if (name.compare("eZH78_HZuR") == 0) {
1417  eZH78_HZuR = value;
1418  } else if (name.compare("eZH78_HZdL") == 0) {
1419  eZH78_HZdL = value;
1420  } else if (name.compare("eZH78_HZdR") == 0) {
1421  eZH78_HZdR = value;
1422  } else if (name.compare("eZH78_ZuL") == 0) {
1423  eZH78_ZuL = value;
1424  } else if (name.compare("eZH78_ZuR") == 0) {
1425  eZH78_ZuR = value;
1426  } else if (name.compare("eZH78_ZdL") == 0) {
1427  eZH78_ZdL = value;
1428  } else if (name.compare("eZH78_ZdR") == 0) {
1429  eZH78_ZdR = value;
1430  } else if (name.compare("ettH2_Htt") == 0) {
1431  ettH2_Htt = value;
1432  } else if (name.compare("ettH2_Hgg") == 0) {
1433  ettH2_Hgg = value;
1434  } else if (name.compare("ettH78_Htt") == 0) {
1435  ettH78_Htt = value;
1436  } else if (name.compare("ettH78_Hgg") == 0) {
1437  ettH78_Hgg = value;
1438  } else if (name.compare("MwInput") == 0)
1439  MwInput = value;
1440  else
1441  NPbase::setParameter(name, value);
1442 }
1443 
1444 bool NPEffectiveGIMR::CheckParameters(const std::map<std::string, double>& DPars)
1445 {
1447  if (FlagMwInput) {
1448  if (DPars.find("MwInput") == DPars.end()) {
1449  std::cout << "ERROR: Missing mandatory NPEffectiveGIMR_LFU_QFU parameter MwInput" << std::endl;
1451  addMissingModelParameter("MwInput");
1452  }
1453  }
1454  for (int i = 0; i < NNPEffectiveGIMRVars_LFU_QFU; i++) {
1455  if (DPars.find(NPEffectiveGIMRVars_LFU_QFU[i]) == DPars.end()) {
1456  std::cout << "ERROR: Missing mandatory NPEffectiveGIMR_LFU_QFU parameter "
1457  << NPEffectiveGIMRVars_LFU_QFU[i] << std::endl;
1460  }
1461  }
1462  //} else if (FlagLeptonUniversal && !FlagQuarkUniversal) {
1463  //} else if (!FlagLeptonUniversal && FlagQuarkUniversal) {
1464  } else if (!FlagLeptonUniversal && !FlagQuarkUniversal) {
1465  if (FlagMwInput) {
1466  if (DPars.find("MwInput") == DPars.end()) {
1467  std::cout << "ERROR: Missing mandatory NPEffectiveGIMR parameter MwInput" << std::endl;
1469  addMissingModelParameter("MwInput");
1470  }
1471  }
1472  for (int i = 0; i < NNPEffectiveGIMRVars; i++) {
1473  if (DPars.find(NPEffectiveGIMRVars[i]) == DPars.end()) {
1474  std::cout << "ERROR: Missing mandatory NPEffectiveGIMR parameter"
1475  << NPEffectiveGIMRVars[i] << std::endl;
1478  }
1479  }
1480  } else
1481  throw std::runtime_error("Error in NPEffectiveGIMR::CheckParameters()");
1482 
1483  return (NPbase::CheckParameters(DPars));
1484 }
1485 
1486 bool NPEffectiveGIMR::setFlag(const std::string name, const bool value)
1487 {
1488  bool res = false;
1489  if (name.compare("MwInput") == 0) {
1490  FlagMwInput = value;
1491  res = true;
1492  } else if (name.compare("QuadraticTerms") == 0) {
1493  FlagQuadraticTerms = value;
1494  if(value) setModelLinearized(false);
1495  res = true;
1496  } else
1497  res = NPbase::setFlag(name, value);
1498 
1499  return (res);
1500 }
1501 
1502 
1504 
1506 {
1507  if (F.is("NEUTRINO_1") || F.is("ELECTRON"))
1508  return CHL1_11;
1509  else if (F.is("NEUTRINO_2") || F.is("MU"))
1510  return CHL1_22;
1511  else if (F.is("NEUTRINO_3") || F.is("TAU"))
1512  return CHL1_33;
1513  else if (F.is("UP") || F.is("DOWN"))
1514  return CHQ1_11;
1515  else if (F.is("CHARM") || F.is("STRANGE"))
1516  return CHQ1_22;
1517  else if (F.is("TOP") || F.is("BOTTOM"))
1518  return CHQ1_33;
1519  else
1520  throw std::runtime_error("NPEffectiveGIMR::CHF1_diag(): wrong argument");
1521 }
1522 
1524 {
1525  if (F.is("NEUTRINO_1") || F.is("ELECTRON"))
1526  return CHL3_11;
1527  else if (F.is("NEUTRINO_2") || F.is("MU"))
1528  return CHL3_22;
1529  else if (F.is("NEUTRINO_3") || F.is("TAU"))
1530  return CHL3_33;
1531  else if (F.is("UP") || F.is("DOWN"))
1532  return CHQ3_11;
1533  else if (F.is("CHARM") || F.is("STRANGE"))
1534  return CHQ3_22;
1535  else if (F.is("TOP") || F.is("BOTTOM"))
1536  return CHQ3_33;
1537  else
1538  throw std::runtime_error("NPEffectiveGIMR::CHF3_diag(): wrong argument");
1539 }
1540 
1542 {
1543  if (f.is("NEUTRINO_1") || f.is("NEUTRINO_2") || f.is("NEUTRINO_3"))
1544  return 0.0;
1545  else if (f.is("ELECTRON"))
1546  return CHe_11;
1547  else if (f.is("MU"))
1548  return CHe_22;
1549  else if (f.is("TAU"))
1550  return CHe_33;
1551  else if (f.is("UP"))
1552  return CHu_11;
1553  else if (f.is("CHARM"))
1554  return CHu_22;
1555  else if (f.is("TOP"))
1556  return CHu_33;
1557  else if (f.is("DOWN"))
1558  return CHd_11;
1559  else if (f.is("STRANGE"))
1560  return CHd_22;
1561  else if (f.is("BOTTOM"))
1562  return CHd_33;
1563  else
1564  throw std::runtime_error("NPEffectiveGIMR::CHf_diag(): wrong argument");
1565 }
1566 
1568 {
1569  if (!u.is("QUARK") || u.getIndex() % 2 != 0)
1570  throw std::runtime_error("NPEffectiveGIMR::CHud_diag(): wrong argument");
1571 
1572  if (u.is("UP"))
1573  return gslpp::complex(CHud_11r, CHud_11i, false);
1574  else if (u.is("CHARM"))
1575  return gslpp::complex(CHud_22r, CHud_22i, false);
1576  else if (u.is("TOP"))
1577  return gslpp::complex(CHud_22r, CHud_33i, false);
1578  else
1579  throw std::runtime_error("NPEffectiveGIMR::CHud_diag(): wrong argument");
1580 }
1581 
1583 {
1584  if (f.is("NEUTRINO_1") || f.is("NEUTRINO_2") || f.is("NEUTRINO_3"))
1585  return 0.0;
1586  else if (f.is("ELECTRON"))
1587  return gslpp::complex(CeH_11r, CeH_11i, false);
1588  else if (f.is("MU"))
1589  return gslpp::complex(CeH_22r, CeH_22i, false);
1590  else if (f.is("TAU"))
1591  return gslpp::complex(CeH_33r, CeH_33i, false);
1592  else if (f.is("UP"))
1593  return gslpp::complex(CuH_11r, CuH_11i, false);
1594  else if (f.is("CHARM"))
1595  return gslpp::complex(CuH_22r, CuH_22i, false);
1596  else if (f.is("TOP"))
1597  return gslpp::complex(CuH_33r, CuH_33i, false);
1598  else if (f.is("DOWN"))
1599  return gslpp::complex(CdH_11r, CdH_11i, false);
1600  else if (f.is("STRANGE"))
1601  return gslpp::complex(CdH_22r, CdH_22i, false);
1602  else if (f.is("BOTTOM"))
1603  return gslpp::complex(CdH_33r, CdH_33i, false);
1604  else
1605  throw std::runtime_error("NPEffectiveGIMR::CfH_diag(): wrong argument");
1606 }
1607 
1609 {
1610  if (f.is("NEUTRINO_1") || f.is("NEUTRINO_2") || f.is("NEUTRINO_3"))
1611  return 0.0;
1612  else if (f.is("ELECTRON"))
1613  return 0.0;
1614  else if (f.is("MU"))
1615  return 0.0;
1616  else if (f.is("TAU"))
1617  return 0.0;
1618  else if (f.is("UP"))
1619  return gslpp::complex(CuG_11r, CuG_11i, false);
1620  else if (f.is("CHARM"))
1621  return gslpp::complex(CuG_22r, CuG_22i, false);
1622  else if (f.is("TOP"))
1623  return gslpp::complex(CuG_33r, CuG_33i, false);
1624  else if (f.is("DOWN"))
1625  return 0.0;
1626  else if (f.is("STRANGE"))
1627  return 0.0;
1628  else if (f.is("BOTTOM"))
1629  return 0.0;
1630  else
1631  throw std::runtime_error("NPEffectiveGIMR::CfG_diag(): wrong argument");
1632 }
1633 
1635 {
1636  if (f.is("NEUTRINO_1") || f.is("NEUTRINO_2") || f.is("NEUTRINO_3"))
1637  return 0.0;
1638  else if (f.is("ELECTRON"))
1639  return 0.0;
1640  else if (f.is("MU"))
1641  return 0.0;
1642  else if (f.is("TAU"))
1643  return 0.0;
1644  else if (f.is("UP"))
1645  return gslpp::complex(CuW_11r, CuW_11i, false);
1646  else if (f.is("CHARM"))
1647  return gslpp::complex(CuW_22r, CuW_22i, false);
1648  else if (f.is("TOP"))
1649  return gslpp::complex(CuW_33r, CuW_33i, false);
1650  else if (f.is("DOWN"))
1651  return 0.0;
1652  else if (f.is("STRANGE"))
1653  return 0.0;
1654  else if (f.is("BOTTOM"))
1655  return 0.0;
1656  else
1657  throw std::runtime_error("NPEffectiveGIMR::CfW_diag(): wrong argument");
1658 }
1659 
1661 {
1662  if (f.is("NEUTRINO_1") || f.is("NEUTRINO_2") || f.is("NEUTRINO_3"))
1663  return 0.0;
1664  else if (f.is("ELECTRON"))
1665  return 0.0;
1666  else if (f.is("MU"))
1667  return 0.0;
1668  else if (f.is("TAU"))
1669  return 0.0;
1670  else if (f.is("UP"))
1671  return gslpp::complex(CuB_11r, CuB_11i, false);
1672  else if (f.is("CHARM"))
1673  return gslpp::complex(CuB_22r, CuB_22i, false);
1674  else if (f.is("TOP"))
1675  return gslpp::complex(CuB_33r, CuB_33i, false);
1676  else if (f.is("DOWN"))
1677  return 0.0;
1678  else if (f.is("STRANGE"))
1679  return 0.0;
1680  else if (f.is("BOTTOM"))
1681  return 0.0;
1682  else
1683  throw std::runtime_error("NPEffectiveGIMR::CfB_diag(): wrong argument");
1684 }
1685 
1686 
1688 
1690 {
1691  return ((CHL3_11 + CHL3_22 - 0.5 * (CLL_1221 + CLL_2112)) * v2_over_LambdaNP2);
1692 }
1693 
1695 {
1696  return (4.0 * sW_tree * cW_tree * CHWB / alphaMz() * v2_over_LambdaNP2);
1697 }
1698 
1700 {
1701  return (-CHD / 2.0 / alphaMz() * v2_over_LambdaNP2);
1702 }
1703 
1705 {
1706  return 0.0;
1707 }
1708 
1710 {
1712 }
1713 
1714 double NPEffectiveGIMR::Mw() const
1715 {
1716  if (FlagMwInput)
1717  return MwInput;
1718  else
1719  return (trueSM.Mw() - Mw_tree() / 4.0 / (cW2_tree - sW2_tree)
1720  *(4.0 * sW_tree * cW_tree * CHWB * v2_over_LambdaNP2
1722  + 2.0 * sW2_tree * DeltaGF()));
1723 }
1724 
1726 {
1727  double G0 = GF * pow(Mw(), 3.0) / 6.0 / sqrt(2.0) / M_PI;
1728  double GammaW_tree = (3.0 + 2.0 * Nc) * G0;
1729 
1730  if (FlagMwInput)
1731  throw std::runtime_error("Write codes in NPEffectiveGIMR::GammaW()!");
1732  else
1733  return (trueSM.GammaW()
1734  - 3.0 * GammaW_tree / 4.0 / (cW2_tree - sW2_tree)
1735  *(4.0 * sW_tree * cW_tree * CHWB * v2_over_LambdaNP2
1737  + 2.0 * (1.0 + cW2_tree) / 3.0 * DeltaGF())
1738  + 2.0 * GammaW_tree / 3.0 * (CHL3_11 + CHQ3_11 + CHQ3_22) * v2_over_LambdaNP2);
1739 }
1740 
1742 {
1743  return (deltaGL_f(p) + deltaGR_f(p));
1744 }
1745 
1747 {
1748  return (deltaGL_f(p) - deltaGR_f(p));
1749 }
1750 
1752 {
1753  double I3p = p.getIsospin(), Qp = p.getCharge();
1754  double CHF1 = CHF1_diag(p);
1755  double CHF3 = CHF3_diag(p);
1756  double NPindirect;
1757  if (FlagMwInput) {
1758  NPindirect = -I3p / 4.0 * (CHD * v2_over_LambdaNP2 + 2.0 * DeltaGF())
1759  + Qp * sW2_tree
1760  * ((cW_tree / sW_tree * CHWB + (1.0 + cW2_tree) / 4.0 / sW2_tree * CHD) * v2_over_LambdaNP2 + 0.5 * DeltaGF());
1761  } else {
1762  NPindirect = -I3p / 4.0 * (CHD * v2_over_LambdaNP2 + 2.0 * DeltaGF())
1763  - Qp * sW2_tree / 4.0 / (cW2_tree - sW2_tree)
1764  *((4.0 * cW_tree / sW_tree * CHWB + CHD) * v2_over_LambdaNP2 + 2.0 * DeltaGF());
1765  }
1766  double NPdirect = -0.5 * (CHF1 - 2.0 * I3p * CHF3) * v2_over_LambdaNP2;
1767  return (NPindirect + NPdirect);
1768 }
1769 
1771 {
1772  double Qp = p.getCharge();
1773  double CHf = CHf_diag(p);
1774  double NPindirect;
1775  if (FlagMwInput) {
1776  NPindirect = Qp * sW2_tree
1777  * ((cW_tree / sW_tree * CHWB + (1.0 + cW2_tree) / 4.0 / sW2_tree * CHD) * v2_over_LambdaNP2 + 0.5 * DeltaGF());
1778  } else {
1779  NPindirect = -Qp * sW2_tree / 4.0 / (cW2_tree - sW2_tree)
1780  *((4.0 * cW_tree / sW_tree * CHWB + CHD) * v2_over_LambdaNP2 + 2.0 * DeltaGF());
1781  }
1782  double NPdirect = -0.5 * CHf*v2_over_LambdaNP2;
1783  return (NPindirect + NPdirect);
1784 }
1785 
1786 
1788 
1790 {
1791  if (pbar.getIndex() + 1 != p.getIndex() || pbar.getIndex() % 2 != 0)
1792  throw std::runtime_error("NPEffectiveGIMR::deltaGL_Wff(): Not implemented");
1793 
1794  double CHF3 = CHF3_diag(pbar);
1795  double NPindirect;
1796  if (FlagMwInput) {
1797  NPindirect = -0.5 * DeltaGF();
1798  } else {
1799  NPindirect = -cW2_tree / 4.0 / (cW2_tree - sW2_tree)
1800  * ((4.0 * sW_tree / cW_tree * CHWB + CHD) * v2_over_LambdaNP2 + 2.0 * DeltaGF());
1801  }
1802  double NPdirect = CHF3 * v2_over_LambdaNP2;
1803  return (NPindirect + NPdirect);
1804 }
1805 
1807 {
1808  if (pbar.getIndex() + 1 != p.getIndex() || pbar.getIndex() % 2 != 0)
1809  throw std::runtime_error("NPEffectiveGIMR::deltaGR_Wff(): Not implemented");
1810 
1811  gslpp::complex CHud = CHud_diag(pbar);
1812  return (0.5 * CHud * v2_over_LambdaNP2);
1813 }
1814 
1816 {
1817  return (CHG * v2_over_LambdaNP2 / v());
1818 }
1819 
1821 {
1822  return (2.0 * CHW * v2_over_LambdaNP2 / v());
1823 }
1824 
1826 {
1827  return 0.0;
1828 }
1829 
1831 {
1832  double NPindirect;
1833  if (FlagMwInput) {
1834  NPindirect = 2.0 * MwInput * MwInput / v() * (delta_h - 0.5 * DeltaGF());
1835  } else {
1836  NPindirect = 2.0 * cW2_tree * Mz * Mz / v()
1837  * (delta_h - 1.0 / 2.0 / (cW2_tree - sW2_tree)
1838  * ((4.0 * sW_tree * cW_tree * CHWB + cW2_tree * CHD) * v2_over_LambdaNP2 + DeltaGF()));
1839  }
1840  return NPindirect;
1841 }
1842 
1844 {
1845  return (delta_ZZ / v());
1846 }
1847 
1849 {
1850  return 0.0;
1851 }
1852 
1854 {
1855  double NPindirect = Mz * Mz / v() * (-0.5 * CHD * v2_over_LambdaNP2 + delta_h - 0.5 * DeltaGF());
1856  double NPdirect = Mz * Mz / v() * CHD * v2_over_LambdaNP2;
1857  return (NPindirect + NPdirect);
1858 }
1859 
1861 {
1862  return (delta_AZ / v());
1863 }
1864 
1866 {
1867  return 0.0;
1868 }
1869 
1871 {
1872  return (delta_AA / v());
1873 }
1874 
1876 {
1877  /* The effects of the RG running are neglected. */
1878  double mf;
1879  if (p.is("TOP"))
1880  //mf = p.getMass(); // m_t(m_t)
1881  mf = mtpole; // pole mass
1882  else
1883  mf = p.getMass();
1884  gslpp::complex CfH = CfH_diag(p);
1885  return (-mf / v() * (delta_h - 0.5 * DeltaGF())
1886  + CfH * v2_over_LambdaNP2 / sqrt(2.0));
1887 }
1888 
1890 {
1891  if (pbar.getIndex() + 1 != p.getIndex() || pbar.getIndex() % 2 != 0)
1892  throw std::runtime_error("NPEffectiveGIMR::deltaGL_Wffh(): Not implemented");
1893 
1894  double CHF3 = CHF3_diag(pbar);
1895  return (2.0 * sqrt(2.0) * Mz * cW_tree / v() / v() * CHF3 * v2_over_LambdaNP2);
1896 }
1897 
1899 {
1900  if (pbar.getIndex() + 1 != p.getIndex() || pbar.getIndex() % 2 != 0)
1901  throw std::runtime_error("NPEffectiveGIMR::deltaGR_Wffh(): Not implemented");
1902 
1903  gslpp::complex CHud = CHud_diag(pbar);
1904  return (sqrt(2.0) * Mz * cW_tree / v() / v() * CHud * v2_over_LambdaNP2);
1905 }
1906 
1908 {
1909  double I3p = p.getIsospin();
1910  double CHF1 = CHF1_diag(p);
1911  double CHF3 = CHF3_diag(p);
1912  return (-2.0 * Mz / v() / v() * (CHF1 - 2.0 * I3p * CHF3) * v2_over_LambdaNP2);
1913 }
1914 
1916 {
1917  double CHf = CHf_diag(p);
1918  return (-2.0 * Mz / v() / v() * CHf * v2_over_LambdaNP2);
1919 }
1920 
1922 {
1923  /* Set to 0. for the moment */
1924 
1925  return 0.;
1926 }
1927 
1929 {
1930  /* Set to 0. for the moment */
1931 
1932  return 0.;
1933 }
1934 
1936 {
1937  /* Set to 0. for the moment */
1938 
1939  return 0.;
1940 }
1941 
1943 {
1944  /* Set to 0. for the moment */
1945 
1946  return 0.;
1947 }
1948 
1950 {
1951  /* Set to 0. for the moment */
1952 
1953  return 0.;
1954 }
1955 
1957 {
1958  /* Set to 0. for the moment */
1959 
1960  return 0.;
1961 }
1962 
1964 {
1965  /* Set to 0. for the moment */
1966 
1967  return 0.;
1968 }
1969 
1970 
1972 
1974 {
1975  gslpp::complex tmp;
1976  if (tau >= 1.0) {
1977  tmp = asin(1.0 / sqrt(tau));
1978  return (tmp * tmp);
1979  } else {
1980  tmp = log((1.0 + sqrt(1.0 - tau)) / (1.0 - sqrt(1.0 - tau))) - M_PI * gslpp::complex::i();
1981  return (-0.25 * tmp * tmp);
1982  }
1983 }
1984 
1985 gslpp::complex NPEffectiveGIMR::AH_f(const double tau) const
1986 {
1987  return (2.0 * tau * (1.0 + (1.0 - tau) * f_triangle(tau)));
1988 }
1989 
1990 double NPEffectiveGIMR::muggH(const double sqrt_s) const
1991 {
1992  double m_t = mtpole;
1993  //doulbe m_t = quarks[TOP].getMass();
1994  //double m_b = quarks[BOTTOM].getMass();
1995 
1996  gslpp::complex dKappa_t = deltaG_hff(quarks[TOP]) / (-m_t / v());
1997  //gslpp::complex dKappa_b = deltaG_hff(quarks[BOTTOM]) / (-m_b / v());
1998 
1999  /* L_eff = G_eff_t_SM*hGG */
2000  gslpp::complex G_eff_t_SM = AlsMz / 16.0 / M_PI / v() * AH_f(4.0 * m_t * m_t / mHl / mHl);
2001 
2002  //double sigma_tt_SM = trueSM.computeSigmaggH_tt(sqrt_s);
2003  //double sigma_bb_SM = trueSM.computeSigmaggH_bb(sqrt_s);
2004  //double sigma_tb_SM = trueSM.computeSigmaggH_tb(sqrt_s);
2005  //gslpp::complex tmp = (2.0 * dKappa_t * sigma_tt_SM
2006  // + 2.0 * dKappa_b * sigma_bb_SM
2007  // + (dKappa_t + dKappa_b) * sigma_tb_SM)
2008  // / (sigma_tt_SM + sigma_bb_SM + sigma_tb_SM);
2009  gslpp::complex tmp = CHG / v() * v2_over_LambdaNP2 / G_eff_t_SM;
2010 
2011  double mu = (1.0 + 2.0 * ( dKappa_t.real() + tmp.real() ) );
2012 
2013  if (FlagQuadraticTerms) {
2014  //Add contributions that are quadratic in the effective coefficients
2015  gslpp::complex tmp2 = dKappa_t + tmp;
2016 
2017  mu += tmp2.abs2();
2018 
2019  }
2020 
2021  if (mu < 0) return std::numeric_limits<double>::quiet_NaN();
2022 
2023  return mu;
2024 }
2025 
2026 double NPEffectiveGIMR::muVBF(const double sqrt_s) const
2027 {
2028  double mu = 1.0;
2029  if (sqrt_s == 1.96) {
2030  mu += +1.123 * (1. + eVBF2_ZuL ) * deltaGL_f(quarks[UP])
2031  -0.531 * (1. + eVBF2_ZuR ) * deltaGR_f(quarks[UP])
2032  -0.705 * (1. + eVBF2_ZdL ) * deltaGL_f(quarks[DOWN])
2033  +0.136 * (1. + eVBF2_ZdR ) * deltaGR_f(quarks[DOWN])
2034  +2.662 * (1. + eVBF2_Wud ) * deltaGL_Wff(quarks[UP],quarks[DOWN]).real()
2035  -1407.72 * (1. + eVBF2_HWud ) * deltaGL_Wffh(quarks[UP], quarks[DOWN]).real()
2036  +14928.1 * (1. + eVBF2_Hgg ) * deltaG_hgg()
2037  -12.451 * (1. + eVBF2_HAA ) * deltaG_hAA()
2038  -21.274 * (1. + eVBF2_HZA1 ) * deltaG1_hZA()
2039  +45.617 * (1. + eVBF2_HZA2 ) * deltaG2_hZA()
2040  -84.016 * (1. + eVBF2_HWW1 ) * deltaG1_hWW()
2041  +390.524 * (1. + eVBF2_HWW2 ) * deltaG2_hWW()
2042  +0.026 * (1. + eVBF2_HWW3 ) * deltaG3_hWW()
2043  -45.832 * (1. + eVBF2_HZZ1 ) * deltaG1_hZZ()
2044  +88.358 * (1. + eVBF2_HZZ2 ) * deltaG2_hZZ()
2045  +0.012 * (1. + eVBF2_HZZ3 ) * deltaG3_hZZ()
2046  -129.338 * (1. + eVBF2_HZuL ) * deltaGL_Zffh(quarks[UP])
2047  +84.325 * (1. + eVBF2_HZuR ) * deltaGR_Zffh(quarks[UP])
2048  +164.195 * (1. + eVBF2_HZdL ) * deltaGL_Zffh(quarks[DOWN])
2049  -32.751 * (1. + eVBF2_HZdR ) * deltaGR_Zffh(quarks[DOWN]);
2050 
2051  if (FlagQuadraticTerms) {
2052  //Add contributions that are quadratic in the effective coefficients
2053  //(Only valid under the assumptions of one dim 6 operator at a time)
2054  mu += +2.478 * pow(deltaGL_f(quarks[UP]),2.0)
2055  +1.878 * pow(deltaGR_f(quarks[UP]),2.0)
2056  +1.214 * pow(deltaGL_f(quarks[DOWN]),2.0)
2057  +0.898 * pow(deltaGR_f(quarks[DOWN]),2.0)
2058  +2.659 * pow(deltaGL_Wff(quarks[UP],quarks[DOWN]).real(),2.0)
2059  +1917816. * pow(deltaGL_Wffh(quarks[UP], quarks[DOWN]).real(),2.0)
2060  +524312994. * pow(deltaG_hgg(),2.0)
2061  +831253. * pow(deltaG_hAA(),2.0)
2062  +151140. * pow(deltaG1_hZA(),2.0)
2063  +58067.7 * pow(deltaG2_hZA(),2.0)
2064  +106835. * pow(deltaG1_hWW(),2.0)
2065  +219369. * pow(deltaG2_hWW(),2.0)
2066  +145840. * pow(deltaG1_hZZ(),2.0)
2067  +66461.2 * pow(deltaG2_hZZ(),2.0)
2068  +1608277. * pow(deltaGL_Zffh(quarks[UP]),2.0)
2069  +1449825. * pow(deltaGR_Zffh(quarks[UP]),2.0)
2070  +409700. * pow(deltaGL_Zffh(quarks[DOWN]),2.0)
2071  +385965. * pow(deltaGR_Zffh(quarks[DOWN]),2.0);
2072  }
2073 
2074  } else if (sqrt_s == 7.0) {
2075  mu += +1.188 * (1. + eVBF78_ZuL ) * deltaGL_f(quarks[UP])
2076  -0.536 * (1. + eVBF78_ZuR ) * deltaGR_f(quarks[UP])
2077  -0.976 * (1. + eVBF78_ZdL ) * deltaGL_f(quarks[DOWN])
2078  +0.179 * (1. + eVBF78_ZdR ) * deltaGR_f(quarks[DOWN])
2079  +2.592 * (1. + eVBF78_Wud ) * deltaGL_Wff(quarks[UP],quarks[DOWN]).real()
2080  -1826.63 * (1. + eVBF78_HWud ) * deltaGL_Wffh(quarks[UP], quarks[DOWN]).real()
2081  +14265.8 * (1. + eVBF78_Hgg ) * deltaG_hgg()
2082  -40.051 * (1. + eVBF78_HAA ) * deltaG_hAA()
2083  -42.43 * (1. + eVBF78_HZA1 ) * deltaG1_hZA()
2084  +88.972 * (1. + eVBF78_HZA2 ) * deltaG2_hZA()
2085  -108.107 * (1. + eVBF78_HWW1 ) * deltaG1_hWW()
2086  +547.508 * (1. + eVBF78_HWW2 ) * deltaG2_hWW()
2087  +0.026 * (1. + eVBF78_HWW3 ) * deltaG3_hWW()
2088  -67.672 * (1. + eVBF78_HZZ1 ) * deltaG1_hZZ()
2089  +168.86 * (1. + eVBF78_HZZ2 ) * deltaG2_hZZ()
2090  +0.014 * (1. + eVBF78_HZZ3 ) * deltaG3_hZZ()
2091  -466.198 * (1. + eVBF78_HZuL ) * deltaGL_Zffh(quarks[UP])
2092  +211.308 * (1. + eVBF78_HZuR ) * deltaGR_Zffh(quarks[UP])
2093  +374.597 * (1. + eVBF78_HZdL ) * deltaGL_Zffh(quarks[DOWN])
2094  -69.916 * (1. + eVBF78_HZdR ) * deltaGR_Zffh(quarks[DOWN]);
2095 
2096  if (FlagQuadraticTerms) {
2097  //Add contributions that are quadratic in the effective coefficients
2098  //(Only valid under the assumptions of one dim 6 operator at a time)
2099  mu += +2.534 * pow(deltaGL_f(quarks[UP]),2.0)
2100  +1.9 * pow(deltaGR_f(quarks[UP]),2.0)
2101  +1.695 * pow(deltaGL_f(quarks[DOWN]),2.0)
2102  +1.177 * pow(deltaGR_f(quarks[DOWN]),2.0)
2103  +2.608 * pow(deltaGL_Wff(quarks[UP],quarks[DOWN]).real(),2.0)
2104  +2862580. * pow(deltaGL_Wffh(quarks[UP], quarks[DOWN]).real(),2.0)
2105  +519301209. * pow(deltaG_hgg(),2.0)
2106  +777159. * pow(deltaG_hAA(),2.0)
2107  +206157. * pow(deltaG1_hZA(),2.0)
2108  +94511.2 * pow(deltaG2_hZA(),2.0)
2109  +174828. * pow(deltaG1_hWW(),2.0)
2110  +414624. * pow(deltaG2_hWW(),2.0)
2111  +209132. * pow(deltaG1_hZZ(),2.0)
2112  +120250. * pow(deltaG2_hZZ(),2.0)
2113  +1311032. * pow(deltaGL_Zffh(quarks[UP]),2.0)
2114  +1130789. * pow(deltaGR_Zffh(quarks[UP]),2.0)
2115  +757088. * pow(deltaGL_Zffh(quarks[DOWN]),2.0)
2116  +651756. * pow(deltaGR_Zffh(quarks[DOWN]),2.0);
2117  }
2118 
2119  } else if (sqrt_s == 8.0) {
2120  mu += +1.179 * (1. + eVBF78_ZuL ) * deltaGL_f(quarks[UP])
2121  -0.532 * (1. + eVBF78_ZuR ) * deltaGR_f(quarks[UP])
2122  -0.984 * (1. + eVBF78_ZdL ) * deltaGL_f(quarks[DOWN])
2123  +0.181 * (1. + eVBF78_ZdR ) * deltaGR_f(quarks[DOWN])
2124  +2.591 * (1. + eVBF78_Wud ) * deltaGL_Wff(quarks[UP],quarks[DOWN]).real()
2125  -1858.03 * (1. + eVBF78_HWud ) * deltaGL_Wffh(quarks[UP], quarks[DOWN]).real()
2126  +14247.4 * (1. + eVBF78_Hgg ) * deltaG_hgg()
2127  -40.46 * (1. + eVBF78_HAA ) * deltaG_hAA()
2128  -41.713 * (1. + eVBF78_HZA1 ) * deltaG1_hZA()
2129  +90.462 * (1. + eVBF78_HZA2 ) * deltaG2_hZA()
2130  -106.576 * (1. + eVBF78_HWW1 ) * deltaG1_hWW()
2131  +562.98 * (1. + eVBF78_HWW2 ) * deltaG2_hWW()
2132  +0.026 * (1. + eVBF78_HWW3 ) * deltaG3_hWW()
2133  -67.57 * (1. + eVBF78_HZZ1 ) * deltaG1_hZZ()
2134  +174.474 * (1. + eVBF78_HZZ2 ) * deltaG2_hZZ()
2135  +0.014 * (1. + eVBF78_HZZ3 ) * deltaG3_hZZ()
2136  -472.887 * (1. + eVBF78_HZuL ) * deltaGL_Zffh(quarks[UP])
2137  +214.739 * (1. + eVBF78_HZuR ) * deltaGR_Zffh(quarks[UP])
2138  +386.582 * (1. + eVBF78_HZdL ) * deltaGL_Zffh(quarks[DOWN])
2139  -72.228 * (1. + eVBF78_HZdR ) * deltaGR_Zffh(quarks[DOWN]);
2140 
2141  if (FlagQuadraticTerms) {
2142  //Add contributions that are quadratic in the effective coefficients
2143  //(Only valid under the assumptions of one dim 6 operator at a time)
2144  mu += +2.503 * pow(deltaGL_f(quarks[UP]),2.0)
2145  +1.877 * pow(deltaGR_f(quarks[UP]),2.0)
2146  +1.712 * pow(deltaGL_f(quarks[DOWN]),2.0)
2147  +1.191 * pow(deltaGR_f(quarks[DOWN]),2.0)
2148  +2.606 * pow(deltaGL_Wff(quarks[UP],quarks[DOWN]).real(),2.0)
2149  +3057041. * pow(deltaGL_Wffh(quarks[UP], quarks[DOWN]).real(),2.0)
2150  +517064803. * pow(deltaG_hgg(),2.0)
2151  +766750. * pow(deltaG_hAA(),2.0)
2152  +207500. * pow(deltaG1_hZA(),2.0)
2153  +101779. * pow(deltaG2_hZA(),2.0)
2154  +177714. * pow(deltaG1_hWW(),2.0)
2155  +454117. * pow(deltaG2_hWW(),2.0)
2156  +210212. * pow(deltaG1_hZZ(),2.0)
2157  +131594. * pow(deltaG2_hZZ(),2.0)
2158  +1399281. * pow(deltaGL_Zffh(quarks[UP]),2.0)
2159  +1231240. * pow(deltaGR_Zffh(quarks[UP]),2.0)
2160  +820259. * pow(deltaGL_Zffh(quarks[DOWN]),2.0)
2161  +713820. * pow(deltaGR_Zffh(quarks[DOWN]),2.0);
2162  }
2163 
2164  } else if (sqrt_s == 14.0) {
2165  mu += +1.129 * deltaGL_f(quarks[UP])
2166  -0.505 * deltaGR_f(quarks[UP])
2167  -1.05 * deltaGL_f(quarks[DOWN])
2168  +0.191 * deltaGR_f(quarks[DOWN])
2169  +2.586 * deltaGL_Wff(quarks[UP],quarks[DOWN]).real()
2170  -1989.34 * deltaGL_Wffh(quarks[UP], quarks[DOWN]).real()
2171  +14228.8 * deltaG_hgg()
2172  -35.554 * deltaG_hAA()
2173  -39.847 * deltaG1_hZA()
2174  +98.522 * deltaG2_hZA()
2175  -99.287 * deltaG1_hWW()
2176  +622.352 * deltaG2_hWW()
2177  +0.026 * deltaG3_hWW()
2178  -66.196 * deltaG1_hZZ()
2179  +196.676 * deltaG2_hZZ()
2180  +0.014 * deltaG3_hZZ()
2181  -493.198 * deltaGL_Zffh(quarks[UP])
2182  +217.017 * deltaGR_Zffh(quarks[UP])
2183  +447.396 * deltaGL_Zffh(quarks[DOWN])
2184  -82.396 * deltaGR_Zffh(quarks[DOWN]);
2185 
2186  if (FlagQuadraticTerms) {
2187  //Add contributions that are quadratic in the effective coefficients
2188  //(Only valid under the assumptions of one dim 6 operator at a time)
2189  mu += +2.319 * pow(deltaGL_f(quarks[UP]),2.0)
2190  +1.783 * pow(deltaGR_f(quarks[UP]),2.0)
2191  +1.849 * pow(deltaGL_f(quarks[DOWN]),2.0)
2192  +1.263 * pow(deltaGR_f(quarks[DOWN]),2.0)
2193  +2.592 * pow(deltaGL_Wff(quarks[UP],quarks[DOWN]).real(),2.0)
2194  +4077238. * pow(deltaGL_Wffh(quarks[UP], quarks[DOWN]).real(),2.0)
2195  +507787376. * pow(deltaG_hgg(),2.0)
2196  +702353. * pow(deltaG_hAA(),2.0)
2197  +212082. * pow(deltaG1_hZA(),2.0)
2198  +141422. * pow(deltaG2_hZA(),2.0)
2199  +195770. * pow(deltaG1_hWW(),2.0)
2200  +655804. * pow(deltaG2_hWW(),2.0)
2201  +0. * pow(deltaG3_hWW(),2.0)
2202  +240333. * pow(deltaG1_hZZ(),2.0)
2203  +192371. * pow(deltaG2_hZZ(),2.0)
2204  +0. * pow(deltaG3_hZZ(),2.0)
2205  +1904757. * pow(deltaGL_Zffh(quarks[UP]),2.0)
2206  +1743849. * pow(deltaGR_Zffh(quarks[UP]),2.0)
2207  +1185212. * pow(deltaGL_Zffh(quarks[DOWN]),2.0)
2208  +1061826. * pow(deltaGR_Zffh(quarks[DOWN]),2.0);
2209  }
2210 
2211  } else
2212  throw std::runtime_error("Bad argument in NPEffectiveGIMR::muVBF()");
2213 
2214  if (mu < 0) return std::numeric_limits<double>::quiet_NaN();
2215 
2216  return mu;
2217 }
2218 
2219 double NPEffectiveGIMR::mueeWBF(const double sqrt_s) const
2220 {
2221  double mu = 1.0;
2222  if (sqrt_s == 0.24) {
2223  mu += 985.974 * deltaG1_hZZ()
2224  +23.622 * deltaG2_hZZ()
2225  +0. * deltaG3_hZZ()
2226  +377.441 * deltaG1_hZA()
2227  -394.144 * deltaG2_hZA()
2228  -30.997 * deltaG1_hWW()
2229  +273.526 * deltaG2_hWW()
2230  +0.038 * deltaG3_hWW()
2231  -62.505 * deltaGL_Zffh(leptons[NEUTRINO_1])
2232  -1781.02 * deltaGL_Zffh(leptons[ELECTRON])
2233  +183.495 * deltaGR_Zffh(leptons[ELECTRON])
2235  -0.062 * deltaGL_f(leptons[NEUTRINO_1])
2236  -1.436 * deltaGL_f(leptons[ELECTRON])
2237  +0.004 * deltaGR_f(leptons[ELECTRON])
2239 
2240  if (FlagQuadraticTerms) {
2241  //Add contributions that are quadratic in the effective coefficients
2242  //(Only valid under the assumptions of one dim 6 operator at a time)
2243  mu += +0.0;
2244  }
2245 
2246  } else if (sqrt_s == 0.25) {
2247  mu += 903.947 * deltaG1_hZZ()
2248  +93.416 * deltaG2_hZZ()
2249  -0.002 * deltaG3_hZZ()
2250  +277.754 * deltaG1_hZA()
2251  -271.135 * deltaG2_hZA()
2252  -29.647 * deltaG1_hWW()
2253  +298.034 * deltaG2_hWW()
2254  +0.039 * deltaG3_hWW()
2255  -61.73 * deltaGL_Zffh(leptons[NEUTRINO_1])
2256  -1409.1 * deltaGL_Zffh(leptons[ELECTRON])
2257  +182.068 * deltaGR_Zffh(leptons[ELECTRON])
2259  -0.174 * deltaGL_f(leptons[NEUTRINO_1])
2260  -1.003 * deltaGL_f(leptons[ELECTRON])
2261  +0. * deltaGR_f(leptons[ELECTRON])
2263 
2264  if (FlagQuadraticTerms) {
2265  //Add contributions that are quadratic in the effective coefficients
2266  //(Only valid under the assumptions of one dim 6 operator at a time)
2267  mu += +0.0;
2268  }
2269 
2270  } else if (sqrt_s == 0.35) {
2271  mu += -63.056 * deltaG1_hZZ()
2272  +99.834 * deltaG2_hZZ()
2273  -0.001 * deltaG3_hZZ()
2274  -100.006 * deltaG1_hZA()
2275  +120.858 * deltaG2_hZA()
2276  -77.085 * deltaG1_hWW()
2277  +424.21 * deltaG2_hWW()
2278  +0.039 * deltaG3_hWW()
2279  -36.176 * deltaGL_Zffh(leptons[NEUTRINO_1])
2280  +108.441 * deltaGL_Zffh(leptons[ELECTRON])
2281  +137.702 * deltaGR_Zffh(leptons[ELECTRON])
2283  -0.091 * deltaGL_f(leptons[NEUTRINO_1])
2284  +0.074 * deltaGL_f(leptons[ELECTRON])
2285  +0.033 * deltaGR_f(leptons[ELECTRON])
2287 
2288  if (FlagQuadraticTerms) {
2289  //Add contributions that are quadratic in the effective coefficients
2290  //(Only valid under the assumptions of one dim 6 operator at a time)
2291  mu += +0.0;
2292  }
2293 
2294  } else if (sqrt_s == 0.5) {
2295  mu += -82.771 * deltaG1_hZZ()
2296  +48.73 * deltaG2_hZZ()
2297  +0. * deltaG3_hZZ()
2298  -78.056 * deltaG1_hZA()
2299  +78.264 * deltaG2_hZA()
2300  -98.794 * deltaG1_hWW()
2301  +579.5 * deltaG2_hWW()
2302  +0.039 * deltaG3_hWW()
2303  -26.448 * deltaGL_Zffh(leptons[NEUTRINO_1])
2304  +163.236 * deltaGL_Zffh(leptons[ELECTRON])
2305  +56.583 * deltaGR_Zffh(leptons[ELECTRON])
2307  -0.02 * deltaGL_f(leptons[NEUTRINO_1])
2308  +0.037 * deltaGL_f(leptons[ELECTRON])
2309  +0.009 * deltaGR_f(leptons[ELECTRON])
2311 
2312  if (FlagQuadraticTerms) {
2313  //Add contributions that are quadratic in the effective coefficients
2314  //(Only valid under the assumptions of one dim 6 operator at a time)
2315  mu += +0.0;
2316  }
2317 
2318  } else if (sqrt_s == 1.0) {
2319  mu += -32.198 * deltaG1_hZZ()
2320  +13.389 * deltaG2_hZZ()
2321  +0. * deltaG3_hZZ()
2322  -27.018 * deltaG1_hZA()
2323  +18.957 * deltaG2_hZA()
2324  -100.42 * deltaG1_hWW()
2325  +884.402 * deltaG2_hWW()
2326  +0.039 * deltaG3_hWW()
2327  -11.556 * deltaGL_Zffh(leptons[NEUTRINO_1])
2328  +56.187 * deltaGL_Zffh(leptons[ELECTRON])
2329  +16.258 * deltaGR_Zffh(leptons[ELECTRON])
2331  -0.001 * deltaGL_f(leptons[NEUTRINO_1])
2332  +0.004 * deltaGL_f(leptons[ELECTRON])
2333  +0.002 * deltaGR_f(leptons[ELECTRON])
2335 
2336  if (FlagQuadraticTerms) {
2337  //Add contributions that are quadratic in the effective coefficients
2338  //(Only valid under the assumptions of one dim 6 operator at a time)
2339  mu += +0.0;
2340  }
2341 
2342  } else
2343  throw std::runtime_error("Bad argument in NPEffectiveGIMR::mueeWBF()");
2344 
2345  if (mu < 0) return std::numeric_limits<double>::quiet_NaN();
2346 
2347  return mu;
2348 }
2349 
2350 double NPEffectiveGIMR::muWH(const double sqrt_s) const
2351 {
2352  double mu = 1.0;
2353  if (sqrt_s == 1.96) {
2354  mu += +2.032 * (1. + eWH2_Wud ) * deltaGL_Wff(quarks[UP], quarks[DOWN]).real()
2355  +1738.87 * (1. + eWH2_HWW1 ) * deltaG1_hWW()
2356  -3432.64 * (1. + eWH2_HWW2 ) * deltaG2_hWW()
2357  +0.039 * (1. + eWH2_HWW3 ) * deltaG3_hWW()
2358  +6523.35 * (1. + eWH2_HWud ) * deltaGL_Wffh(quarks[UP], quarks[DOWN]).real();
2359 
2360  if (FlagQuadraticTerms) {
2361  //Add contributions that are quadratic in the effective coefficients
2362  //(Only valid under the assumptions of one dim 6 operator at a time)
2363  mu += +1.042 * pow(deltaGL_Wff(quarks[UP], quarks[DOWN]).real(),2.0)
2364  +1075949. * pow(deltaG1_hWW(),2.0)
2365  +3978950. * pow(deltaG2_hWW(),2.0)
2366  +15887131. * pow(deltaGL_Wffh(quarks[UP], quarks[DOWN]).real(),2.0);
2367  }
2368 
2369  } else if (sqrt_s == 7.0) {
2370  mu += +1.979 * (1. + eWH78_Wud ) * deltaGL_Wff(quarks[UP], quarks[DOWN]).real()
2371  +1777.77 * (1. + eWH78_HWW1 ) * deltaG1_hWW()
2372  -3890.65 * (1. + eWH78_HWW2 ) * deltaG2_hWW()
2373  +0.039 * (1. + eWH78_HWW3 ) * deltaG3_hWW()
2374  +7344.73 * (1. + eWH78_HWud ) * deltaGL_Wffh(quarks[UP], quarks[DOWN]).real();
2375 
2376  if (FlagQuadraticTerms) {
2377  //Add contributions that are quadratic in the effective coefficients
2378  //(Only valid under the assumptions of one dim 6 operator at a time)
2379  mu += +1.015 * pow(deltaGL_Wff(quarks[UP], quarks[DOWN]).real(),2.0)
2380  +1294405. * pow(deltaG1_hWW(),2.0)
2381  +7356224. * pow(deltaG2_hWW(),2.0)
2382  +31355627. * pow(deltaGL_Wffh(quarks[UP], quarks[DOWN]).real(),2.0);
2383  }
2384 
2385  } else if (sqrt_s == 8.0) {
2386  mu += +1.978 * (1. + eWH78_Wud ) * deltaGL_Wff(quarks[UP], quarks[DOWN]).real()
2387  +1784.47 * (1. + eWH78_HWW1 ) * deltaG1_hWW()
2388  -3967.38 * (1. + eWH78_HWW2 ) * deltaG2_hWW()
2389  +0.039 * (1. + eWH78_HWW3 ) * deltaG3_hWW()
2390  +7507.02 * (1. + eWH78_HWud ) * deltaGL_Wffh(quarks[UP], quarks[DOWN]).real();
2391 
2392  if (FlagQuadraticTerms) {
2393  //Add contributions that are quadratic in the effective coefficients
2394  //(Only valid under the assumptions of one dim 6 operator at a time)
2395  mu += +1.016 * pow(deltaGL_Wff(quarks[UP], quarks[DOWN]).real(),2.0)
2396  +1331512. * pow(deltaG1_hWW(),2.0)
2397  +8168916. * pow(deltaG2_hWW(),2.0)
2398  +35201222. * pow(deltaGL_Wffh(quarks[UP], quarks[DOWN]).real(),2.0);
2399  }
2400 
2401  } else if (sqrt_s == 14.0) {
2402  mu += +1.963 * deltaGL_Wff(quarks[UP], quarks[DOWN]).real()
2403  +1799.45 * deltaG1_hWW()
2404  -4252.03 * deltaG2_hWW()
2405  +0.039 * deltaG3_hWW()
2406  +8047.59 * deltaGL_Wffh(quarks[UP], quarks[DOWN]).real();
2407 
2408  if (FlagQuadraticTerms) {
2409  //Add contributions that are quadratic in the effective coefficients
2410  //(Only valid under the assumptions of one dim 6 operator at a time)
2411  mu += +1.007 * pow(deltaGL_Wff(quarks[UP], quarks[DOWN]).real(),2.0)
2412  +1467903. * pow(deltaG1_hWW(),2.0)
2413  +13173439. * pow(deltaG2_hWW(),2.0)
2414  +58780336. * pow(deltaGL_Wffh(quarks[UP], quarks[DOWN]).real(),2.0);
2415  }
2416 
2417  } else
2418  throw std::runtime_error("Bad argument in NPEffectiveGIMR::muWH()");
2419 
2420  if (mu < 0) return std::numeric_limits<double>::quiet_NaN();
2421 
2422  return mu;
2423 }
2424 
2425 double NPEffectiveGIMR::muZH(const double sqrt_s) const
2426 {
2427  double mu = 1.0;
2428  if (sqrt_s == 1.96) {
2429  mu += +3.529 * (1. + eZH2_ZuL ) * deltaGL_f(quarks[UP])
2430  -1.598 * (1. + eZH2_ZuR ) * deltaGR_f(quarks[UP])
2431  -1.229 * (1. + eZH2_ZdL ) * deltaGL_f(quarks[DOWN])
2432  +0.227 * (1. + eZH2_ZdR ) * deltaGR_f(quarks[DOWN])
2433  +3215.38 * (1. + eZH2_HZZ1 ) * deltaG1_hZZ()
2434  -2922.42 * (1. + eZH2_HZZ2 ) * deltaG2_hZZ()
2435  +0.059 * (1. + eZH2_HZZ3 ) * deltaG3_hZZ()
2436  +495.399 * (1. + eZH2_HZA1 ) * deltaG1_hZA()
2437  -838.743 * (1. + eZH2_HZA2 ) * deltaG2_hZA()
2438  +5931.99 * (1. + eZH2_HZuL ) * deltaGL_Zffh(quarks[UP])
2439  -2684.23 * (1. + eZH2_HZuR ) * deltaGR_Zffh(quarks[UP])
2440  -1878.46 * (1. + eZH2_HZdL ) * deltaGL_Zffh(quarks[DOWN])
2441  +346.694 * (1. + eZH2_HZdR ) * deltaGR_Zffh(quarks[DOWN]);
2442 
2443  if (FlagQuadraticTerms) {
2444  //Add contributions that are quadratic in the effective coefficients
2445  //(Only valid under the assumptions of one dim 6 operator at a time)
2446  mu += +5.126 * pow(deltaGL_f(quarks[UP]),2.0)
2447  +5.126 * pow(deltaGR_f(quarks[UP]),2.0)
2448  +1.456 * pow(deltaGL_f(quarks[DOWN]),2.0)
2449  +1.454 * pow(deltaGR_f(quarks[DOWN]),2.0)
2450  +3525123. * pow(deltaG1_hZZ(),2.0)
2451  +2844179. * pow(deltaG2_hZZ(),2.0)
2452  +0.001 * pow(deltaG3_hZZ(),2.0)
2453  +662397. * pow(deltaG1_hZA(),2.0)
2454  +2006248. * pow(deltaG2_hZA(),2.0)
2455  +21799545. * pow(deltaGL_Zffh(quarks[UP]),2.0)
2456  +21795795. * pow(deltaGR_Zffh(quarks[UP]),2.0)
2457  +4723149. * pow(deltaGL_Zffh(quarks[DOWN]),2.0)
2458  +4725123. * pow(deltaGR_Zffh(quarks[DOWN]),2.0);
2459  }
2460 
2461  } else if (sqrt_s == 7.0) {
2462  mu += +2.583 * (1. + eZH78_ZuL ) * deltaGL_f(quarks[UP])
2463  -1.17 * (1. + eZH78_ZuR ) * deltaGR_f(quarks[UP])
2464  -2.127 * (1. + eZH78_ZdL ) * deltaGL_f(quarks[DOWN])
2465  +0.392 * (1. + eZH78_ZdR ) * deltaGR_f(quarks[DOWN])
2466  +3269.53 * (1. + eZH78_HZZ1 ) * deltaG1_hZZ()
2467  -3201.65 * (1. + eZH78_HZZ2 ) * deltaG2_hZZ()
2468  +0.059 * (1. + eZH78_HZZ3 ) * deltaG3_hZZ()
2469  +473.267 * (1. + eZH78_HZA1 ) * deltaG1_hZA()
2470  -873.421 * (1. + eZH78_HZA2 ) * deltaG2_hZA()
2471  +4763.44 * (1. + eZH78_HZuL ) * deltaGL_Zffh(quarks[UP])
2472  -2156.99 * (1. + eZH78_HZuR ) * deltaGR_Zffh(quarks[UP])
2473  -3853.2 * (1. + eZH78_HZdL ) * deltaGL_Zffh(quarks[DOWN])
2474  +712.124 * (1. + eZH78_HZdR ) * deltaGR_Zffh(quarks[DOWN]);
2475 
2476  if (FlagQuadraticTerms) {
2477  //Add contributions that are quadratic in the effective coefficients
2478  //(Only valid under the assumptions of one dim 6 operator at a time)
2479  mu += +3.752 * pow(deltaGL_f(quarks[UP]),2.0)
2480  +3.753 * pow(deltaGR_f(quarks[UP]),2.0)
2481  +2.519 * pow(deltaGL_f(quarks[DOWN]),2.0)
2482  +2.517 * pow(deltaGR_f(quarks[DOWN]),2.0)
2483  +4051505. * pow(deltaG1_hZZ(),2.0)
2484  +4597749. * pow(deltaG2_hZZ(),2.0)
2485  +0.001 * pow(deltaG3_hZZ(),2.0)
2486  +610510. * pow(deltaG1_hZA(),2.0)
2487  +2766996. * pow(deltaG2_hZA(),2.0)
2488  +27425400. * pow(deltaGL_Zffh(quarks[UP]),2.0)
2489  +27416894. * pow(deltaGR_Zffh(quarks[UP]),2.0)
2490  +17043782. * pow(deltaGL_Zffh(quarks[DOWN]),2.0)
2491  +17039528. * pow(deltaGR_Zffh(quarks[DOWN]),2.0);
2492  }
2493 
2494  } else if (sqrt_s == 8.0) {
2495  mu += +2.569 * (1. + eZH78_ZuL ) * deltaGL_f(quarks[UP])
2496  -1.163 * (1. + eZH78_ZuR ) * deltaGR_f(quarks[UP])
2497  -2.14 * (1. + eZH78_ZdL ) * deltaGL_f(quarks[DOWN])
2498  +0.395 * (1. + eZH78_ZdR ) * deltaGR_f(quarks[DOWN])
2499  +3282.79 * (1. + eZH78_HZZ1 ) * deltaG1_hZZ()
2500  -3262.46 * (1. + eZH78_HZZ2 ) * deltaG2_hZZ()
2501  +0.059 * (1. + eZH78_HZZ3 ) * deltaG3_hZZ()
2502  +475.044 * (1. + eZH78_HZA1 ) * deltaG1_hZA()
2503  -892.243 * (1. + eZH78_HZA2 ) * deltaG2_hZA()
2504  +4847.78 * (1. + eZH78_HZuL ) * deltaGL_Zffh(quarks[UP])
2505  -2193.61 * (1. + eZH78_HZuR ) * deltaGR_Zffh(quarks[UP])
2506  -3960.46 * (1. + eZH78_HZdL ) * deltaGL_Zffh(quarks[DOWN])
2507  +731.438 * (1. + eZH78_HZdR ) * deltaGR_Zffh(quarks[DOWN]);
2508 
2509  if (FlagQuadraticTerms) {
2510  //Add contributions that are quadratic in the effective coefficients
2511  //(Only valid under the assumptions of one dim 6 operator at a time)
2512  mu += +3.732 * pow(deltaGL_f(quarks[UP]),2.0)
2513  +3.736 * pow(deltaGR_f(quarks[UP]),2.0)
2514  +2.535 * pow(deltaGL_f(quarks[DOWN]),2.0)
2515  +2.536 * pow(deltaGR_f(quarks[DOWN]),2.0)
2516  +4164701. * pow(deltaG1_hZZ(),2.0)
2517  +5067698. * pow(deltaG2_hZZ(),2.0)
2518  +0.001 * pow(deltaG3_hZZ(),2.0)
2519  +627966. * pow(deltaG1_hZA(),2.0)
2520  +3087745. * pow(deltaG2_hZA(),2.0)
2521  +30566228. * pow(deltaGL_Zffh(quarks[UP]),2.0)
2522  +30559313. * pow(deltaGR_Zffh(quarks[UP]),2.0)
2523  +19107837. * pow(deltaGL_Zffh(quarks[DOWN]),2.0)
2524  +19109134. * pow(deltaGR_Zffh(quarks[DOWN]),2.0);
2525  }
2526 
2527  } else if (sqrt_s == 14.0) {
2528  mu += +2.477 * deltaGL_f(quarks[UP])
2529  -1.103 * deltaGR_f(quarks[UP])
2530  -2.226 * deltaGL_f(quarks[DOWN])
2531  +0.405 * deltaGR_f(quarks[DOWN])
2532  +3321.75 * deltaG1_hZZ()
2533  -3494.38 * deltaG2_hZZ()
2534  +0.059 * deltaG3_hZZ()
2535  +481.727 * deltaG1_hZA()
2536  -967.239 * deltaG2_hZA()
2537  +5106.92 * deltaGL_Zffh(quarks[UP])
2538  -2270.81 * deltaGR_Zffh(quarks[UP])
2539  -4434.64 * deltaGL_Zffh(quarks[DOWN])
2540  +807.186 * deltaGR_Zffh(quarks[DOWN]);
2541  if (FlagQuadraticTerms) {
2542  //Add contributions that are quadratic in the effective coefficients
2543  //(Only valid under the assumptions of one dim 6 operator at a time)
2544  mu += +3.579 * pow(deltaGL_f(quarks[UP]),2.0)
2545  +3.58 * pow(deltaGR_f(quarks[UP]),2.0)
2546  +2.631 * pow(deltaGL_f(quarks[DOWN]),2.0)
2547  +2.629 * pow(deltaGR_f(quarks[DOWN]),2.0)
2548  +4609160. * pow(deltaG1_hZZ(),2.0)
2549  +7946470. * pow(deltaG2_hZZ(),2.0)
2550  +0.001 * pow(deltaG3_hZZ(),2.0)
2551  +683466. * pow(deltaG1_hZA(),2.0)
2552  +5019397. * pow(deltaG2_hZA(),2.0)
2553  +50036976. * pow(deltaGL_Zffh(quarks[UP]),2.0)
2554  +50008570. * pow(deltaGR_Zffh(quarks[UP]),2.0)
2555  +31660707. * pow(deltaGL_Zffh(quarks[DOWN]),2.0)
2556  +31666009. * pow(deltaGR_Zffh(quarks[DOWN]),2.0);
2557  }
2558 
2559  } else
2560  throw std::runtime_error("Bad argument in NPEffectiveGIMR::muZH()");
2561 
2562  if (mu < 0) return std::numeric_limits<double>::quiet_NaN();
2563 
2564  return mu;
2565 }
2566 
2567 double NPEffectiveGIMR::mueeZH(const double sqrt_s) const
2568 {
2569  double mu = 1.0;
2570 
2571  if (sqrt_s == 0.24) {
2572  mu += 2690.84 * deltaG1_hZZ()
2573  -1951.93 * deltaG2_hZZ()
2574  +0.059 * deltaG3_hZZ()
2575  +147.761 * deltaG1_hZA()
2576  -185.735 * deltaG2_hZA()
2577  -4217.73 * deltaGL_Zffh(leptons[ELECTRON])
2578  +3619.82 * deltaGR_Zffh(leptons[ELECTRON])
2579  -4.282 * deltaGL_f(leptons[ELECTRON])
2580  +3.674 * deltaGR_f(leptons[ELECTRON]);
2581 
2582  if (FlagQuadraticTerms) {
2583  //Add contributions that are quadratic in the effective coefficients
2584  //(Only valid under the assumptions of one dim 6 operator at a time)
2585  mu += +7.966 * pow(deltaGL_f(leptons[DOWN]),2.0)
2586  +7.966 * pow(deltaGR_f(leptons[DOWN]),2.0)
2587  +1841343. * pow(deltaG1_hZZ(),2.0)
2588  +952412. * pow(deltaG2_hZZ(),2.0)
2589  +0.001 * pow(deltaG3_hZZ(),2.0)
2590  +961714. * pow(deltaG1_hZA(),2.0)
2591  +1520521. * pow(deltaG2_hZA(),2.0)
2592  +7731703. * pow(deltaGL_Zffh(leptons[DOWN]),2.0)
2593  +7731703. * pow(deltaGR_Zffh(leptons[DOWN]),2.0);
2594  }
2595 
2596  } else if (sqrt_s == 0.25) {
2597  mu += 2829.45 * deltaG1_hZZ()
2598  -2097.01 * deltaG2_hZZ()
2599  +0.059 * deltaG3_hZZ()
2600  +156.787 * deltaG1_hZA()
2601  -204.357 * deltaG2_hZA()
2602  -4635.23 * deltaGL_Zffh(leptons[ELECTRON])
2603  +3979.84 * deltaGR_Zffh(leptons[ELECTRON])
2604  -4.282 * deltaGL_f(leptons[ELECTRON])
2605  +3.674 * deltaGR_f(leptons[ELECTRON]);
2606 
2607  if (FlagQuadraticTerms) {
2608  //Add contributions that are quadratic in the effective coefficients
2609  //(Only valid under the assumptions of one dim 6 operator at a time)
2610  mu += +0.0;
2611  }
2612 
2613  } else if (sqrt_s == 0.35) {
2614  mu += 3893.41 * deltaG1_hZZ()
2615  -3873.83 * deltaG2_hZZ()
2616  +0.059 * deltaG3_hZZ()
2617  +231.963 * deltaG1_hZA()
2618  -424.266 * deltaG2_hZA()
2619  -9763.89 * deltaGL_Zffh(leptons[ELECTRON])
2620  +8387.72 * deltaGR_Zffh(leptons[ELECTRON])
2621  -4.282 * deltaGL_f(leptons[ELECTRON])
2622  +3.674 * deltaGR_f(leptons[ELECTRON]);
2623 
2624  if (FlagQuadraticTerms) {
2625  //Add contributions that are quadratic in the effective coefficients
2626  //(Only valid under the assumptions of one dim 6 operator at a time)
2627  mu += +0.0;
2628  }
2629 
2630  } else if (sqrt_s == 0.5) {
2631  mu += 4747.11 * deltaG1_hZZ()
2632  -7649.28 * deltaG2_hZZ()
2633  +0.059 * deltaG3_hZZ()
2634  +291.854 * deltaG1_hZA()
2635  -902.663 * deltaG2_hZA()
2636  -20668.9 * deltaGL_Zffh(leptons[ELECTRON])
2637  +17754.2 * deltaGR_Zffh(leptons[ELECTRON])
2638  -4.282 * deltaGL_f(leptons[ELECTRON])
2639  +3.674 * deltaGR_f(leptons[ELECTRON]);
2640 
2641  if (FlagQuadraticTerms) {
2642  //Add contributions that are quadratic in the effective coefficients
2643  //(Only valid under the assumptions of one dim 6 operator at a time)
2644  mu += +0.0;
2645  }
2646 
2647  } else if (sqrt_s == 1.0) {
2648  mu += 5576.18 * deltaG1_hZZ()
2649  -29856.9 * deltaG2_hZZ()
2650  +0.059 * deltaG3_hZZ()
2651  +351.186 * deltaG1_hZA()
2652  -3727.98 * deltaG2_hZA()
2653  -84814.3 * deltaGL_Zffh(leptons[ELECTRON])
2654  +72844.9 * deltaGR_Zffh(leptons[ELECTRON])
2655  -4.282 * deltaGL_f(leptons[ELECTRON])
2656  +3.673 * deltaGR_f(leptons[ELECTRON]);
2657 
2658  if (FlagQuadraticTerms) {
2659  //Add contributions that are quadratic in the effective coefficients
2660  //(Only valid under the assumptions of one dim 6 operator at a time)
2661  mu += +0.0;
2662  }
2663 
2664  } else
2665  throw std::runtime_error("Bad argument in NPEffectiveGIMR::mueeZH()");
2666 
2667  if (mu < 0) return std::numeric_limits<double>::quiet_NaN();
2668 
2669  return mu;
2670 }
2671 
2672 double NPEffectiveGIMR::muVH(const double sqrt_s) const
2673 {
2674  double sigmaWH_SM = computeSigmaWH(sqrt_s);
2675  double sigmaZH_SM = computeSigmaZH(sqrt_s);
2676  double sigmaWH = muWH(sqrt_s) * sigmaWH_SM;
2677  double sigmaZH = muZH(sqrt_s) * sigmaZH_SM;
2678  double mu = ((sigmaWH + sigmaZH) / (sigmaWH_SM + sigmaZH_SM));
2679 
2680  if (mu < 0) return std::numeric_limits<double>::quiet_NaN();
2681 
2682  return mu;
2683 }
2684 
2685 double NPEffectiveGIMR::muVBFpVH(const double sqrt_s) const
2686 {
2687  double sigmaWH_SM = computeSigmaWH(sqrt_s);
2688  double sigmaZH_SM = computeSigmaZH(sqrt_s);
2689  double sigmaVBF_SM = computeSigmaVBF(sqrt_s);
2690  double sigmaWH = muWH(sqrt_s) * sigmaWH_SM;
2691  double sigmaZH = muZH(sqrt_s) * sigmaZH_SM;
2692  double sigmaVBF = muVBF(sqrt_s) * sigmaVBF_SM;
2693  double mu = ((sigmaWH + sigmaZH + sigmaVBF) / (sigmaWH_SM + sigmaZH_SM + sigmaVBF_SM));
2694 
2695  if (mu < 0) return std::numeric_limits<double>::quiet_NaN();
2696 
2697  return mu;
2698 }
2699 
2700 double NPEffectiveGIMR::muttH(const double sqrt_s) const
2701 {
2702  double mu = 1.0;
2703  if (sqrt_s == 1.96) {
2704  mu += -2.863 * (1. + ettH2_Htt ) * deltaG_hff(quarks[TOP]).real()
2705  +1737.35 * (1. + ettH2_Hgg ) * deltaG_hgg();
2706 
2707  if (FlagQuadraticTerms) {
2708  //Add contributions that are quadratic in the effective coefficients
2709  //(Only valid under the assumptions of one dim 6 operator at a time)
2710  mu += +2.036 * pow(deltaG_hff(quarks[TOP]).real(),2.0)
2711  +885586. * pow(deltaG_hgg(),2.0);
2712  }
2713 
2714  } else if (sqrt_s == 7.0) {
2715  mu += -2.861 * (1. + ettH78_Htt ) * deltaG_hff(quarks[TOP]).real()
2716  +2583.3 * (1. + ettH78_Hgg ) * deltaG_hgg();
2717 
2718  if (FlagQuadraticTerms) {
2719  //Add contributions that are quadratic in the effective coefficients
2720  //(Only valid under the assumptions of one dim 6 operator at a time)
2721  mu += +2.073 * pow(deltaG_hff(quarks[TOP]).real(),2.0)
2722  +3909554. * pow(deltaG_hgg(),2.0);
2723  }
2724 
2725  } else if (sqrt_s == 8.0) {
2726  mu += -2.861 * (1. + ettH78_Htt ) * deltaG_hff(quarks[TOP]).real()
2727  +2636.88 * (1. + ettH78_Hgg ) * deltaG_hgg();
2728 
2729  if (FlagQuadraticTerms) {
2730  //Add contributions that are quadratic in the effective coefficients
2731  //(Only valid under the assumptions of one dim 6 operator at a time)
2732  mu += +1.963 * pow(deltaG_hff(quarks[TOP]).real(),2.0)
2733  +4367338. * pow(deltaG_hgg(),2.0);
2734  }
2735 
2736  } else if (sqrt_s == 14.0) {
2737  mu += -2.861 * deltaG_hff(quarks[TOP]).real()
2738  +2769.79 * deltaG_hgg();
2739 
2740  if (FlagQuadraticTerms) {
2741  //Add contributions that are quadratic in the effective coefficients
2742  //(Only valid under the assumptions of one dim 6 operator at a time)
2743  mu += +2.012 * pow(deltaG_hff(quarks[TOP]).real(),2.0)
2744  +5689423. * pow(deltaG_hgg(),2.0);
2745  }
2746 
2747  } else
2748  throw std::runtime_error("Bad argument in NPEffectiveGIMR::muttH()");
2749 
2750  if (mu < 0) return std::numeric_limits<double>::quiet_NaN();
2751 
2752  return mu;
2753 }
2754 
2755 double NPEffectiveGIMR::muggHpttH(const double sqrt_s) const
2756 {
2757  double sigmaggH_SM = computeSigmaggH(sqrt_s);
2758  double sigmattH_SM = computeSigmattH(sqrt_s);
2759  double sigmaggH = muggH(sqrt_s) * sigmaggH_SM;
2760  double sigmattH = muttH(sqrt_s) * sigmattH_SM;
2761 
2762  double mu = ((sigmaggH + sigmattH) / (sigmaggH_SM + sigmattH_SM));
2763 
2764  if (mu < 0) return std::numeric_limits<double>::quiet_NaN();
2765 
2766  return mu;
2767 }
2768 
2769 double NPEffectiveGIMR::mueettH(const double sqrt_s) const
2770 {
2771  double mu = 1.0;
2772  if (sqrt_s == 0.5) {
2773  mu += 85.139 * deltaG1_hZZ()
2774  -51.41 * deltaG2_hZZ()
2775  +0. * deltaG3_hZZ()
2776  +276.673 * deltaG1_hZA()
2777  -159.708 * deltaG2_hZA()
2778  +1017.44 * deltaG_hAA()
2779  -2.833 * deltaG_hff(quarks[TOP]).real()
2780  -287.92 * deltaGL_Zffh(leptons[ELECTRON])
2781  -132.849 * deltaGR_Zffh(leptons[ELECTRON])
2782  +84.883 * deltaGL_Zffh(quarks[TOP])
2783  +72.935 * deltaGR_Zffh(quarks[TOP])
2784  +863716. * deltaG_hAff(quarks[TOP]).real()
2785  +154393. * deltaG_hZff(quarks[TOP]).real()
2786  -1.046 * deltaGL_f(leptons[ELECTRON])
2787  -0.608 * deltaGR_f(leptons[ELECTRON])
2788  +0.663 * deltaGL_f(quarks[TOP])
2789  +0.585 * deltaGR_f(quarks[TOP])
2790  +8833.35 * deltaG_Aff(quarks[TOP]).real()
2791  +1650.94 * deltaG_Zff(quarks[TOP]).real();
2792 
2793  if (FlagQuadraticTerms) {
2794  //Add contributions that are quadratic in the effective coefficients
2795  //(Only valid under the assumptions of one dim 6 operator at a time)
2796  mu += +0.0;
2797  }
2798 
2799  } else if (sqrt_s == 1.0) {
2800  mu += 446.758 * deltaG1_hZZ()
2801  -1500.92 * deltaG2_hZZ()
2802  +0.003 * deltaG3_hZZ()
2803  +657.283 * deltaG1_hZA()
2804  -572.102 * deltaG2_hZA()
2805  +2443.18 * deltaG_hAA()
2806  -2.701 * deltaG_hff(quarks[TOP]).real()
2807  -4591.53 * deltaGL_Zffh(leptons[ELECTRON])
2808  +2945.96 * deltaGR_Zffh(leptons[ELECTRON])
2809  +251.003 * deltaGL_Zffh(quarks[TOP])
2810  +49.581 * deltaGR_Zffh(quarks[TOP])
2811  +3025550. * deltaG_hAff(quarks[TOP]).real()
2812  +519896. * deltaG_hZff(quarks[TOP]).real()
2813  -1.426 * deltaGL_f(leptons[ELECTRON])
2814  -0.041 * deltaGR_f(leptons[ELECTRON])
2815  +1.066 * deltaGL_f(quarks[TOP])
2816  -0.038 * deltaGR_f(quarks[TOP])
2817  +12745.4 * deltaG_Aff(quarks[TOP]).real()
2818  +2238.61 * deltaG_Zff(quarks[TOP]).real();
2819 
2820  if (FlagQuadraticTerms) {
2821  //Add contributions that are quadratic in the effective coefficients
2822  //(Only valid under the assumptions of one dim 6 operator at a time)
2823  mu += +0.0;
2824  }
2825 
2826  } else
2827  throw std::runtime_error("Bad argument in NPEffectiveGIMR::mueettH()");
2828 
2829  if (mu < 0) return std::numeric_limits<double>::quiet_NaN();
2830 
2831  return mu;
2832 }
2833 
2835 {
2836  double Br = 1.0;
2837 
2839 
2840  if (FlagQuadraticTerms) {
2841  //Add contributions that are quadratic in the effective coefficients
2842  //(Only valid under the assumptions of one dim 6 operator at a time)
2845  + pow(deltaGammaTotalRatio1(),2.0);
2846  }
2847 
2848  if (Br < 0) return std::numeric_limits<double>::quiet_NaN();
2849 
2850  return Br;
2851 
2852 }
2853 
2855 {
2856  double Br = 1.0;
2857 
2859 
2860  if (FlagQuadraticTerms) {
2861  //Add contributions that are quadratic in the effective coefficients
2862  //(Only valid under the assumptions of one dim 6 operator at a time)
2865  + pow(deltaGammaTotalRatio1(),2.0);
2866  }
2867 
2868  if (Br < 0) return std::numeric_limits<double>::quiet_NaN();
2869 
2870  return Br;
2871 
2872 }
2873 
2875 {
2876  double Br = 1.0;
2877 
2879 
2880  if (FlagQuadraticTerms) {
2881  //Add contributions that are quadratic in the effective coefficients
2882  //(Only valid under the assumptions of one dim 6 operator at a time)
2885  + pow(deltaGammaTotalRatio1(),2.0);
2886  }
2887 
2888  if (Br < 0) return std::numeric_limits<double>::quiet_NaN();
2889 
2890  return Br;
2891 
2892 }
2893 
2895 {
2896  double Br = 1.0;
2897 
2899 
2900  if (FlagQuadraticTerms) {
2901  //Add contributions that are quadratic in the effective coefficients
2902  //(Only valid under the assumptions of one dim 6 operator at a time)
2905  + pow(deltaGammaTotalRatio1(),2.0);
2906  }
2907 
2908  if (Br < 0) return std::numeric_limits<double>::quiet_NaN();
2909 
2910  return Br;
2911 
2912 }
2913 
2915 {
2916  double Br = 1.0;
2917 
2919 
2920  if (FlagQuadraticTerms) {
2921  //Add contributions that are quadratic in the effective coefficients
2922  //(Only valid under the assumptions of one dim 6 operator at a time)
2925  + pow(deltaGammaTotalRatio1(),2.0);
2926  }
2927 
2928  if (Br < 0) return std::numeric_limits<double>::quiet_NaN();
2929 
2930  return Br;
2931 
2932 }
2933 
2935 {
2936  double Br = 1.0;
2937 
2939 
2940  if (FlagQuadraticTerms) {
2941  //Add contributions that are quadratic in the effective coefficients
2942  //(Only valid under the assumptions of one dim 6 operator at a time)
2945  + pow(deltaGammaTotalRatio1(),2.0);
2946  }
2947 
2948  if (Br < 0) return std::numeric_limits<double>::quiet_NaN();
2949 
2950  return Br;
2951 
2952 }
2953 
2955 {
2956  double Br = 1.0;
2957 
2959 
2960  if (FlagQuadraticTerms) {
2961  //Add contributions that are quadratic in the effective coefficients
2962  //(Only valid under the assumptions of one dim 6 operator at a time)
2965  + pow(deltaGammaTotalRatio1(),2.0);
2966  }
2967 
2968  if (Br < 0) return std::numeric_limits<double>::quiet_NaN();
2969 
2970  return Br;
2971 
2972 }
2973 
2975 {
2976  double Br = 1.0;
2977 
2979 
2980  if (FlagQuadraticTerms) {
2981  //Add contributions that are quadratic in the effective coefficients
2982  //(Only valid under the assumptions of one dim 6 operator at a time)
2985  + pow(deltaGammaTotalRatio1(),2.0);
2986  }
2987 
2988  if (Br < 0) return std::numeric_limits<double>::quiet_NaN();
2989 
2990  return Br;
2991 
2992 }
2993 
2995 {
2996  double Br = 1.0;
2997 
2999 
3000  if (FlagQuadraticTerms) {
3001  //Add contributions that are quadratic in the effective coefficients
3002  //(Only valid under the assumptions of one dim 6 operator at a time)
3005  + pow(deltaGammaTotalRatio1(),2.0);
3006  }
3007 
3008  if (Br < 0) return std::numeric_limits<double>::quiet_NaN();
3009 
3010  return Br;
3011 
3012 }
3013 
3015 {
3016  return (trueSM.computeBrHtogg() * GammaHggRatio()
3025 }
3026 
3028 {
3038 }
3039 
3041 {
3051 }
3052 
3054 {
3055  double width = 1.0;
3056 
3057  width += deltaGammaHggRatio1();
3058 
3059  if (FlagQuadraticTerms) {
3060  //Add contributions that are quadratic in the effective coefficients
3061  //(Only valid under the assumptions of one dim 6 operator at a time)
3062  width += deltaGammaHggRatio2();
3063  }
3064 
3065  return width;
3066 
3067 }
3068 
3070 {
3071  return ( +151669. * deltaG_hgg()
3072  -3.006 * deltaG_hff(quarks[TOP]).real()
3073  +5.853 * deltaG_hff(quarks[BOTTOM]).real()
3074  +4.71 * deltaG_hff(quarks[CHARM]).real() );
3075 }
3076 
3078 {
3079  //Contributions that are quadratic in the effective coefficients
3080  //(Only valid under the assumptions of one dim 6 operator at a time)
3081  return ( +5879800851. * pow(deltaG_hgg(),2.0)
3082  +2.284 * pow(deltaG_hff(quarks[TOP]).real(),2.0)
3083  +40.881 * pow(deltaG_hff(quarks[BOTTOM]).real(),2.0)
3084  +2.17 * pow(deltaG_hff(quarks[CHARM]).real(),2.0) );
3085 
3086 }
3087 
3089 {
3090  double width = 1.0;
3091 
3092  width += deltaGammaHWWRatio1();
3093 
3094  if (FlagQuadraticTerms) {
3095  //Add contributions that are quadratic in the effective coefficients
3096  //(Only valid under the assumptions of one dim 6 operator at a time)
3097  width += deltaGammaHWWRatio2();
3098  }
3099 
3100  return width;
3101 
3102 }
3103 
3105 {
3106 
3107  return ( -183.404 * deltaG1_hWW()
3108  -274.568 * deltaG2_hWW()
3109  +0.039 * deltaG3_hWW() );
3110 
3111 }
3112 
3114 {
3115  //Contributions that are quadratic in the effective coefficients
3116  //(Only valid under the assumptions of one dim 6 operator at a time)
3117  return ( +1267. * pow(deltaG1_hWW(),2.0)
3118  +868.393 * pow(deltaG2_hWW(),2.0) );
3119 
3120 }
3121 
3123 {
3124  double width = 1.0;
3125 
3126  width += deltaGammaHZZRatio1();
3127 
3128  if (FlagQuadraticTerms) {
3129  //Add contributions that are quadratic in the effective coefficients
3130  //(Only valid under the assumptions of one dim 6 operator at a time)
3131  width += deltaGammaHZZRatio2();
3132  }
3133 
3134  return width;
3135 
3136 }
3137 
3139 {
3140 
3141  return ( -246.654 * deltaG1_hZZ()
3142  -240.846 * deltaG2_hZZ()
3143  +0.059 * deltaG3_hZZ() );
3144 
3145 }
3146 
3148 {
3149  //Contributions that are quadratic in the effective coefficients
3150  //(Only valid under the assumptions of one dim 6 operator at a time)
3151  return ( +6391.57 * pow(deltaG1_hZZ(),2.0)
3152  +2088.67 * pow(deltaG2_hZZ(),2.0)
3153  +0.001 * pow(deltaG3_hZZ(),2.0) );
3154 
3155 }
3156 
3158 {
3159  double width = 1.0;
3160 
3161  width += deltaGammaHZgaRatio1();
3162 
3163  if (FlagQuadraticTerms) {
3164  //Add contributions that are quadratic in the effective coefficients
3165  //(Only valid under the assumptions of one dim 6 operator at a time)
3166  width += deltaGammaHZgaRatio2();
3167  }
3168 
3169  return width;
3170 
3171 }
3172 
3174 {
3175 
3176  return ( -71321.5 * deltaG1_hZA()
3177  +0.041 * deltaG3_hWW()
3178  +0.172 * deltaG_hff(quarks[TOP]).real()
3179  -0.301 * deltaG_hff(quarks[BOTTOM]).real()
3180  +0.196 * deltaG_hff(leptons[TAU]).real()
3181  +0.232 * deltaG_hff(quarks[CHARM]).real() );
3182 
3183 }
3184 
3186 {
3187  //Contributions that are quadratic in the effective coefficients
3188  //(Only valid under the assumptions of one dim 6 operator at a time)
3189  return ( +1271853409. * pow(deltaG1_hZA(),2.0)
3190  +0.003 * pow(deltaG_hff(quarks[TOP]).real(),2.0)
3191  +3.539 * pow(deltaG_hff(quarks[BOTTOM]).real(),2.0)
3192  -14.568 * pow(deltaG_hff(leptons[TAU]).real(),2.0)
3193  -31.197 * pow(deltaG_hff(quarks[CHARM]).real(),2.0) );
3194 
3195 }
3196 
3198 {
3199  double width = 1.0;
3200 
3201  width += deltaGammaHgagaRatio1();
3202 
3203  if (FlagQuadraticTerms) {
3204  //Add contributions that are quadratic in the effective coefficients
3205  //(Only valid under the assumptions of one dim 6 operator at a time)
3206  width += deltaGammaHgagaRatio2();
3207  }
3208 
3209  return width;
3210 
3211 }
3212 
3214 {
3215  return ( -257366. * deltaG_hAA()
3216  +0.049 * deltaG3_hWW()
3217  +0.761 * deltaG_hff(quarks[TOP]).real()
3218  -0.441 * deltaG_hff(quarks[BOTTOM]).real()
3219  -1.087 * deltaG_hff(leptons[TAU]).real()
3220  -0.646 * deltaG_hff(quarks[CHARM]).real() );
3221 
3222 }
3223 
3225 {
3226  //Contributions that are quadratic in the effective coefficients
3227  //(Only valid under the assumptions of one dim 6 operator at a time)
3228  return ( +16479108529. * pow(deltaG_hAA(),2.0)
3229  +0.001 * pow(deltaG3_hWW(),2.0)
3230  +0.146 * pow(deltaG_hff(quarks[TOP]).real(),2.0)
3231  +1.828 * pow(deltaG_hff(quarks[BOTTOM]).real(),2.0)
3232  +6.672 * pow(deltaG_hff(leptons[TAU]).real(),2.0)
3233  +9.962 * pow(deltaG_hff(quarks[CHARM]).real(),2.0) );
3234 
3235 }
3236 
3238 {
3239  double width = 1.0;
3240 
3241  width += deltaGammaHmumuRatio1();
3242 
3243  if (FlagQuadraticTerms) {
3244  //Add contributions that are quadratic in the effective coefficients
3245  //(Only valid under the assumptions of one dim 6 operator at a time)
3246  width += deltaGammaHmumuRatio2();
3247  }
3248 
3249  return width;
3250 
3251 }
3252 
3254 {
3255  return ( -4653.43 * deltaG_hff(leptons[MU]).real() );
3256 
3257 }
3258 
3260 {
3261  //Contributions that are quadratic in the effective coefficients
3262  //(Only valid under the assumptions of one dim 6 operator at a time)
3263  return 0.0;
3264 
3265 }
3266 
3268 {
3269  double width = 1.0;
3270 
3271  width += deltaGammaHtautauRatio1();
3272 
3273  if (FlagQuadraticTerms) {
3274  //Add contributions that are quadratic in the effective coefficients
3275  //(Only valid under the assumptions of one dim 6 operator at a time)
3276  width += deltaGammaHtautauRatio2();
3277  }
3278 
3279  return width;
3280 
3281 }
3282 
3284 {
3285  return ( -277.458 * deltaG_hff(leptons[TAU]).real() );
3286 
3287 }
3288 
3290 {
3291  //Contributions that are quadratic in the effective coefficients
3292  //(Only valid under the assumptions of one dim 6 operator at a time)
3293  return ( +19223. * pow(deltaG_hff(leptons[TAU]).real(),2.0) );
3294 
3295 }
3296 
3298 {
3299  double width = 1.0;
3300 
3301  width += deltaGammaHccRatio1();
3302 
3303  if (FlagQuadraticTerms) {
3304  //Add contributions that are quadratic in the effective coefficients
3305  //(Only valid under the assumptions of one dim 6 operator at a time)
3306  width += deltaGammaHccRatio2();
3307  }
3308 
3309  return width;
3310 
3311 }
3312 
3314 {
3315  return ( -383.036 * deltaG_hff(quarks[CHARM]).real() );
3316 }
3317 
3319 {
3320  //Contributions that are quadratic in the effective coefficients
3321  //(Only valid under the assumptions of one dim 6 operator at a time)
3322  return ( +36709.1 * pow(deltaG_hff(quarks[CHARM]).real(),2.0) );
3323 
3324 }
3325 
3327 {
3328  double width = 1.0;
3329 
3330  width += deltaGammaHbbRatio1();
3331 
3332  if (FlagQuadraticTerms) {
3333  //Add contributions that are quadratic in the effective coefficients
3334  //(Only valid under the assumptions of one dim 6 operator at a time)
3335  width += deltaGammaHbbRatio2();
3336  }
3337 
3338  return width;
3339 }
3340 
3342 {
3343  return ( -0.013 * deltaG_hff(quarks[TOP]).real()
3344  -117.431 * deltaG_hff(quarks[BOTTOM]).real() );
3345 }
3346 
3348 {
3349  //Contributions that are quadratic in the effective coefficients
3350  //(Only valid under the assumptions of one dim 6 operator at a time)
3351  return ( +3443.96 * pow(deltaG_hff(quarks[BOTTOM]).real(),2.0) );
3352 
3353 }
3354 
3356 
3358 {
3359  return ((CLL_1122 + CLL_1221) / LambdaNP2);
3360 }
3361 
3363 {
3364  return ((CLL_1133 + CLL_1331) / LambdaNP2);
3365 }
3366 
3368 {
3369  return ((CLQ1_1111 - CLQ3_1111) / LambdaNP2);
3370 }
3371 
3373 {
3374  return ((CLQ1_1111 + CLQ3_1111) / LambdaNP2);
3375 }
3376 
3378 {
3379  return ((CLQ1_1122 - CLQ3_1122) / LambdaNP2);
3380 }
3381 
3383 {
3384  return ((CLQ1_1122 + CLQ3_1122) / LambdaNP2);
3385 }
3386 
3388 {
3389  return ((CLQ1_1133 + CLQ3_1133) / LambdaNP2);
3390 }
3391 
3393 {
3394  return ((CLe_1122 + CLe_2211) / LambdaNP2);
3395 }
3396 
3398 {
3399  return ((CLe_1133 + CLe_3311) / LambdaNP2);
3400 }
3401 
3403 {
3404  return (CLu_1111 / LambdaNP2);
3405 }
3406 
3408 {
3409  return (CLd_1111 / LambdaNP2);
3410 }
3411 
3413 {
3414  return (CLu_1122 / LambdaNP2);
3415 }
3416 
3418 {
3419  return (CLd_1122 / LambdaNP2);
3420 }
3421 
3423 {
3424  return (CLd_1133 / LambdaNP2);
3425 }
3426 
3428 {
3429  return (CLe_1122 / LambdaNP2);
3430 }
3431 
3433 {
3434  return (CLe_1133 / LambdaNP2);
3435 }
3436 
3438 {
3439  return (CQe_1111 / LambdaNP2);
3440 }
3441 
3443 {
3444  return (CQe_1111 / LambdaNP2);
3445 }
3446 
3448 {
3449  return (CQe_1122 / LambdaNP2);
3450 }
3451 
3453 {
3454  return (CQe_1122 / LambdaNP2);
3455 }
3456 
3458 {
3459  return (CQe_1133 / LambdaNP2);
3460 }
3461 
3463 {
3464  return (Cee_1122 / LambdaNP2);
3465 }
3466 
3468 {
3469  return (Cee_1133 / LambdaNP2);
3470 }
3471 
3472 
3474 {
3475  return (Ceu_1111 / LambdaNP2);
3476 }
3477 
3479 {
3480  return (Ced_1111 / LambdaNP2);
3481 }
3482 
3484 {
3485  return (Ceu_1122 / LambdaNP2);
3486 }
3487 
3489 {
3490  return (Ced_1122 / LambdaNP2);
3491 }
3492 
3494 {
3495  return (Ced_1133 / LambdaNP2);
3496 }
3497 
3498 
3499 
3500 double NPEffectiveGIMR::LEP2sigmaMu(const double s) const
3501 {
3502  double sigma_mu;
3503  double GIMRParam[10] = {CLL_mu(), CLR_mu(), CRL_mu(), CRR_mu(), deltaGamma_Z(),
3504  deltaGR_f(leptons[MU]),
3505  deltaGL_f(leptons[MU]),
3508  deltaMz2()};
3509 
3510  sigma_mu = trueSM.LEP2sigmaMu(s) +
3511  myLEP2GIMR.sigma_l_LEP2_GIMR(QCD::lepton(MU), s, GIMRParam);
3512 
3513  return sigma_mu;
3514 }
3515 
3516 double NPEffectiveGIMR::LEP2sigmaTau(const double s) const
3517 {
3518  double sigma_tau;
3519  double GIMRParam[10] = {CLL_tau(), CLR_tau(), CRL_tau(), CRR_tau(), deltaGamma_Z(),
3520  deltaGR_f(leptons[TAU]),
3521  deltaGL_f(leptons[TAU]),
3524  deltaMz2()};
3525 
3526  sigma_tau = trueSM.LEP2sigmaTau(s) +
3527  myLEP2GIMR.sigma_l_LEP2_GIMR(QCD::lepton(TAU), s, GIMRParam);
3528 
3529  return sigma_tau;
3530 }
3531 
3532 double NPEffectiveGIMR::LEP2sigmaHadron(const double s) const
3533 {
3534  double sigma_had;
3535  double GIMRParamUP[10] = {CLL_up(), CLR_up(), CRL_up(), CRR_up(), deltaGamma_Z(),
3536  deltaGR_f(quarks[UP]),
3537  deltaGL_f(quarks[UP]),
3540  deltaMz2()};
3541 
3542  double GIMRParamDOWN[10] = {CLL_down(), CLR_down(), CRL_down(), CRR_down(), deltaGamma_Z(),
3543  deltaGR_f(quarks[DOWN]),
3544  deltaGL_f(quarks[DOWN]),
3547  deltaMz2()};
3548 
3549  double GIMRParamSTRANGE[10] = {CLL_strange(), CLR_strange(), CRL_strange(), CRR_strange(), deltaGamma_Z(),
3554  deltaMz2()};
3555 
3556  double GIMRParamCHARM[10] = {CLL_charm(), CLR_charm(), CRL_charm(), CRR_charm(), deltaGamma_Z(),
3557  deltaGR_f(quarks[CHARM]),
3561  deltaMz2()};
3562 
3563  double GIMRParamBOTTOM[10] = {CLL_bottom(), CLR_bottom(), CRL_bottom(), CRR_bottom(), deltaGamma_Z(),
3564  deltaGR_f(quarks[BOTTOM]),
3568  deltaMz2()};
3569 
3570 
3571  sigma_had = trueSM.LEP2sigmaHadron(s) +
3572  myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::UP, s, GIMRParamUP) +
3573  myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::DOWN, s, GIMRParamDOWN) +
3574  myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::CHARM, s, GIMRParamCHARM) +
3575  myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::STRANGE, s, GIMRParamSTRANGE) +
3576  myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::BOTTOM, s, GIMRParamBOTTOM);
3577 
3578  return sigma_had;
3579 }
3580 
3581 
3582 double NPEffectiveGIMR::LEP2AFBmu(const double s) const
3583 {
3584  double AFB_mu = 0., sigma = 0.0,dFminusB =0.0;
3585  double GIMRParam[10] = {CLL_mu(), CLR_mu(), CRL_mu(), CRR_mu(), deltaGamma_Z(),
3586  deltaGR_f(leptons[MU]),
3587  deltaGL_f(leptons[MU]),
3590  deltaMz2()};
3591 
3592  sigma = trueSM.LEP2sigmaMu(s);
3593  dFminusB = myLEP2GIMR.sigmaFminusB_l_LEP2_GIMR(QCD::lepton(MU), s, GIMRParam);
3594 
3595  AFB_mu = trueSM.LEP2AFBmu(s) * (1. - myLEP2GIMR.sigma_l_LEP2_GIMR(QCD::lepton(MU), s, GIMRParam) / sigma)
3596  + dFminusB / sigma;
3597 
3598  return AFB_mu;
3599 }
3600 
3601 double NPEffectiveGIMR::LEP2AFBtau(const double s) const
3602 {
3603  double AFB_tau, sigma = 0.0,dFminusB =0.0;
3604  double GIMRParam[10] = {CLL_tau(), CLR_tau(), CRL_tau(), CRR_tau(), deltaGamma_Z(),
3605  deltaGR_f(leptons[TAU]),
3606  deltaGL_f(leptons[TAU]),
3609  deltaMz2()};
3610  sigma = trueSM.LEP2sigmaTau(s);
3611  dFminusB = myLEP2GIMR.sigmaFminusB_l_LEP2_GIMR(QCD::lepton(TAU), s, GIMRParam);
3612 
3613  AFB_tau = trueSM.LEP2AFBtau(s) * (1. - myLEP2GIMR.sigma_l_LEP2_GIMR(QCD::lepton(TAU), s, GIMRParam) / sigma)
3614  + dFminusB / sigma;
3615 
3616  return AFB_tau;
3617 }
3618 
3619 double NPEffectiveGIMR::LEP2AFBcharm(const double s) const
3620 {
3621  double AFB_charm, sigma = 0.0,dFminusB =0.0;
3622  double GIMRParam[10] = {CLL_charm(), CLR_charm(), CRL_charm(), CRR_charm(), deltaGamma_Z(),
3623  deltaGR_f(quarks[CHARM]),
3627  deltaMz2()};
3628  sigma = trueSM.LEP2sigmaCharm(s);
3629  dFminusB = myLEP2GIMR.sigmaFminusB_q_LEP2_GIMR(QCD::quark(CHARM), s, GIMRParam);
3630 
3631  AFB_charm = trueSM.LEP2AFBcharm(s) * (1. - myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::quark(CHARM), s, GIMRParam) / sigma)
3632  + dFminusB / sigma;
3633 
3634  return AFB_charm;
3635 }
3636 
3637 double NPEffectiveGIMR::LEP2AFBbottom(const double s) const
3638 {
3639  double AFB_bottom, sigma = 0.0,dFminusB =0.0;;
3640  double GIMRParam[10] = {CLL_bottom(), CLR_bottom(), CRL_bottom(), CRR_bottom(), deltaGamma_Z(),
3641  deltaGR_f(quarks[BOTTOM]),
3645  deltaMz2()};
3646  sigma = trueSM.LEP2sigmaBottom(s);
3647  dFminusB = myLEP2GIMR.sigmaFminusB_q_LEP2_GIMR(QCD::quark(BOTTOM), s, GIMRParam);
3648 
3649  AFB_bottom = trueSM.LEP2AFBbottom(s) * (1. - myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::quark(BOTTOM), s, GIMRParam) / sigma)
3650  + dFminusB / sigma;
3651 
3652  return AFB_bottom;
3653 }
3654 
3655 
3656 double NPEffectiveGIMR::LEP2Rcharm(const double s) const
3657 {
3658  double R_charm = 0., dsigma_c = 0.0, dsigma_had =0.0, sigma_hSM = 0.0;
3659  double GIMRParamUP[10] = {CLL_up(), CLR_up(), CRL_up(), CRR_up(), deltaGamma_Z(),
3660  deltaGR_f(quarks[UP]),
3661  deltaGL_f(quarks[UP]),
3664  deltaMz2()};
3665 
3666  double GIMRParamDOWN[10] = {CLL_down(), CLR_down(), CRL_down(), CRR_down(), deltaGamma_Z(),
3667  deltaGR_f(quarks[DOWN]),
3668  deltaGL_f(quarks[DOWN]),
3671  deltaMz2()};
3672 
3673  double GIMRParamSTRANGE[10] = {CLL_strange(), CLR_strange(), CRL_strange(), CRR_strange(), deltaGamma_Z(),
3678  deltaMz2()};
3679 
3680  double GIMRParamCHARM[10] = {CLL_charm(), CLR_charm(), CRL_charm(), CRR_charm(), deltaGamma_Z(),
3681  deltaGR_f(quarks[CHARM]),
3685  deltaMz2()};
3686 
3687  double GIMRParamBOTTOM[10] = {CLL_bottom(), CLR_bottom(), CRL_bottom(), CRR_bottom(), deltaGamma_Z(),
3688  deltaGR_f(quarks[BOTTOM]),
3692  deltaMz2()};
3693 
3694  sigma_hSM = trueSM.LEP2sigmaHadron(s);
3695 
3696  dsigma_had = myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::UP, s, GIMRParamUP) +
3697  myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::DOWN, s, GIMRParamDOWN) +
3698  myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::CHARM, s, GIMRParamCHARM) +
3699  myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::STRANGE, s, GIMRParamSTRANGE) +
3700  myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::BOTTOM, s, GIMRParamBOTTOM);
3701 
3702  dsigma_c = myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::CHARM, s, GIMRParamCHARM);
3703 
3704  R_charm = trueSM.LEP2Rcharm(s) * (1 - dsigma_had / sigma_hSM) + dsigma_c / sigma_hSM;
3705 
3706  return R_charm;
3707 }
3708 
3709 double NPEffectiveGIMR::LEP2Rbottom(const double s) const
3710 {
3711  double R_bottom = 0., dsigma_b = 0.0, dsigma_had =0.0, sigma_hSM = 0.0;
3712  double GIMRParamUP[10] = {CLL_up(), CLR_up(), CRL_up(), CRR_up(), deltaGamma_Z(),
3713  deltaGR_f(quarks[UP]),
3714  deltaGL_f(quarks[UP]),
3717  deltaMz2()};
3718 
3719  double GIMRParamDOWN[10] = {CLL_down(), CLR_down(), CRL_down(), CRR_down(), deltaGamma_Z(),
3720  deltaGR_f(quarks[DOWN]),
3721  deltaGL_f(quarks[DOWN]),
3724  deltaMz2()};
3725 
3726  double GIMRParamSTRANGE[10] = {CLL_strange(), CLR_strange(), CRL_strange(), CRR_strange(), deltaGamma_Z(),
3731  deltaMz2()};
3732 
3733  double GIMRParamCHARM[10] = {CLL_charm(), CLR_charm(), CRL_charm(), CRR_charm(), deltaGamma_Z(),
3734  deltaGR_f(quarks[CHARM]),
3738  deltaMz2()};
3739 
3740  double GIMRParamBOTTOM[10] = {CLL_bottom(), CLR_bottom(), CRL_bottom(), CRR_bottom(), deltaGamma_Z(),
3741  deltaGR_f(quarks[BOTTOM]),
3745  deltaMz2()};
3746 
3747  sigma_hSM = trueSM.LEP2sigmaHadron(s);
3748 
3749  dsigma_had = myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::UP, s, GIMRParamUP) +
3750  myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::DOWN, s, GIMRParamDOWN) +
3751  myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::CHARM, s, GIMRParamCHARM) +
3752  myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::STRANGE, s, GIMRParamSTRANGE) +
3753  myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::BOTTOM, s, GIMRParamBOTTOM);
3754 
3755  dsigma_b = myLEP2GIMR.sigma_q_LEP2_GIMR(QCD::BOTTOM, s, GIMRParamCHARM);
3756 
3757  R_bottom = trueSM.LEP2Rbottom(s) * (1 - dsigma_had / sigma_hSM) + dsigma_b / sigma_hSM;
3758 
3759  return R_bottom;
3760 }
3761 
3762 
QCD::TAU
Definition: QCD.h:316
NPEffectiveGIMR::CuG_11i
double CuG_11i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1448
NPEffectiveGIMR::eWH2_HWW1
double eWH2_HWW1
Theoretical uncertainty in the (linear) new physics contribution from to WH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1561
NPEffectiveGIMR::AH_f
gslpp::complex AH_f(const double tau) const
Fermionic loop function entering in the calculation of the effective and couplings.
Definition: NPEffectiveGIMR.cpp:1985
NPEffectiveGIMR::CeH_11r
double CeH_11r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1406
sigmattH
Definition: NPSMEFT6dtopquark.h:659
NPEffectiveGIMR::deltaGammaHbbRatio2
double deltaGammaHbbRatio2() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3347
NPEffectiveGIMR::CuB_23i
double CuB_23i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1476
NPEffectiveGIMR::CHL1_13i
double CHL1_13i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1338
NPEffectiveGIMR::deltaGammaHccRatio1
double deltaGammaHccRatio1() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3313
NPEffectiveGIMR::eZH78_HZuL
double eZH78_HZuL
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1590
NPEffectiveGIMR::GammaHZZRatio
double GammaHZZRatio() const
The ratio of the in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:3122
NPEffectiveGIMR::delta_ZZ
double delta_ZZ
Combination of dimension 6 coefficients modifying the canonical field definition.
Definition: NPEffectiveGIMR.h:1612
NPEffectiveGIMR::CuB_23r
double CuB_23r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1470
NPEffectiveGIMR::deltaMz2
virtual double deltaMz2() const
Definition: NPEffectiveGIMR.cpp:1709
NPEffectiveGIMR::muVBFpVH
virtual double muVBFpVH(const double sqrt_s) const
The ratio between the sum of VBF and WH+ZH associated production cross-section in the current model ...
Definition: NPEffectiveGIMR.cpp:2685
NPEffectiveGIMR::CLe_1111
double CLe_1111
Definition: NPEffectiveGIMR.h:1502
NPEffectiveGIMR::MwInput
double MwInput
The input value for the -boson mass if FlagMwInput is true.
Definition: NPEffectiveGIMR.h:1604
NPEffectiveGIMR::LEP2AFBbottom
virtual double LEP2AFBbottom(const double s) const
Definition: NPEffectiveGIMR.cpp:3637
NPEffectiveGIMR::CHf_diag
double CHf_diag(const Particle f) const
The diagonal entry of the dimension-6 operator coefficient corresponding to particle f.
Definition: NPEffectiveGIMR.cpp:1541
NPEffectiveGIMR::deltaG_hff
virtual gslpp::complex deltaG_hff(const Particle p) const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1875
NPEffectiveGIMR::CHQ1_12i
double CHQ1_12i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1364
StandardModel::setParameter
virtual void setParameter(const std::string name, const double &value)
A method to set the value of a parameter of StandardModel.
Definition: StandardModel.cpp:257
NPEffectiveGIMR::CLQ3_1111
double CLQ3_1111
Definition: NPEffectiveGIMR.h:1486
StandardModel::v
virtual double v() const
The Higgs vacuum expectation value.
Definition: StandardModel.cpp:943
NPEffectiveGIMR::CLd_2211
double CLd_2211
Definition: NPEffectiveGIMR.h:1509
NPEffectiveGIMR::eZH2_HZZ3
double eZH2_HZZ3
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1574
NPEffectiveGIMR::CHL1_12i
double CHL1_12i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1337
NPEffectiveGIMR::CuW_23i
double CuW_23i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1464
NPEffectiveGIMR::eVBF2_HZdL
double eVBF2_HZdL
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1532
NPEffectiveGIMR::CHud_22r
double CHud_22r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1397
NPEffectiveGIMR::CeH_12i
double CeH_12i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1413
NPEffectiveGIMR::eVBF2_HAA
double eVBF2_HAA
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1525
NPEffectiveGIMR::CLe_3311
double CLe_3311
Definition: NPEffectiveGIMR.h:1504
NPEffectiveGIMR::LEP2Rcharm
virtual double LEP2Rcharm(const double s) const
Definition: NPEffectiveGIMR.cpp:3656
NPEffectiveGIMR::Ced_1111
double Ced_1111
Definition: NPEffectiveGIMR.h:1497
QCD::BOTTOM
Definition: QCD.h:329
NPEffectiveGIMR::CHQ1_13r
double CHQ1_13r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1360
NPEffectiveGIMR::eVBF2_HZA2
double eVBF2_HZA2
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1524
NPEffectiveGIMR::eWH78_Wud
double eWH78_Wud
Theoretical uncertainty in the (linear) new physics contribution from to WH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1570
NPEffectiveGIMR::sW_tree
double sW_tree
The tree level values for the sine of the weak angle.
Definition: NPEffectiveGIMR.h:1609
NPEffectiveGIMR::CRL_bottom
double CRL_bottom() const
Definition: NPEffectiveGIMR.cpp:3457
NPEffectiveGIMR::deltaGammaHZgaRatio2
double deltaGammaHZgaRatio2() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3185
NPEffectiveGIMR::eVBF2_HZdR
double eVBF2_HZdR
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1533
NPEffectiveGIMR::CHQ3_13i
double CHQ3_13i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1374
NPEffectiveGIMR::deltaG_Gff
gslpp::complex deltaG_Gff(const Particle p) const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1942
Particle::is
bool is(std::string name_i) const
Definition: Particle.cpp:23
Particle
A class for particles.
Definition: Particle.h:26
NPEffectiveGIMR::deltaGammaHggRatio1
double deltaGammaHggRatio1() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3069
NPEffectiveGIMR::eVBF78_ZuR
double eVBF78_ZuR
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1556
StandardModel::computeSigmaWH
double computeSigmaWH(const double sqrt_s) const
The WH production cross section in the Standard Model.
Definition: StandardModel.h:2102
NPEffectiveGIMR::eZH2_HZA2
double eZH2_HZA2
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1576
NPEffectiveGIMR::CuB_13r
double CuB_13r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1468
StandardModel::LEP2sigmaTau
virtual double LEP2sigmaTau(const double s) const
Definition: StandardModel.cpp:2674
NPEffectiveGIMR::delta_AZ
double delta_AZ
Combination of dimension 6 coefficients modifying the canonical field definition.
Definition: NPEffectiveGIMR.h:1614
NPEffectiveGIMR::deltaG_hZff
gslpp::complex deltaG_hZff(const Particle p) const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1928
NPEffectiveGIMR::CLQ1_1123
double CLQ1_1123
Definition: NPEffectiveGIMR.h:1484
NPEffectiveGIMR::eVBF2_HZZ1
double eVBF2_HZZ1
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1520
NPEffectiveGIMR::eVBF2_Wud
double eVBF2_Wud
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1539
NPEffectiveGIMR::CuH_11i
double CuH_11i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1424
NPEffectiveGIMR::CheckParameters
virtual bool CheckParameters(const std::map< std::string, double > &DPars)
A method to check if all the mandatory parameters for NPEffectiveGIMR have been provided in model ini...
Definition: NPEffectiveGIMR.cpp:1444
NPEffectiveGIMR::deltaGammaHZZRatio2
double deltaGammaHZZRatio2() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3147
NPEffectiveGIMR::CuH_22r
double CuH_22r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1421
NPEffectiveGIMR::eVBF78_HWud
double eVBF78_HWud
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1554
NPEffectiveGIMR::CLQ3_1221
double CLQ3_1221
Definition: NPEffectiveGIMR.h:1487
NPEffectiveGIMR::CHL1_23r
double CHL1_23r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1335
NPEffectiveGIMR::deltaGV_f
virtual double deltaGV_f(const Particle p) const
New physics contribution to the neutral-current vector coupling .
Definition: NPEffectiveGIMR.cpp:1741
StandardModel::computeBrHtotautau
double computeBrHtotautau() const
The Br in the Standard Model.
Definition: StandardModel.h:2278
NPEffectiveGIMR::setFlag
virtual bool setFlag(const std::string name, const bool value)
A method to set a flag of NPEffectiveGIMR.
Definition: NPEffectiveGIMR.cpp:1486
NPEffectiveGIMR::deltaGammaHmumuRatio2
double deltaGammaHmumuRatio2() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3259
NPEffectiveGIMR::CHL1_11
double CHL1_11
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1331
NPEffectiveGIMR::eWH78_HWW3
double eWH78_HWW3
Theoretical uncertainty in the (linear) new physics contribution from to WH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1568
NPEffectiveGIMR::deltaGammaHWWRatio2
double deltaGammaHWWRatio2() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3113
NPEffectiveGIMR::CLL_1111
double CLL_1111
Definition: NPEffectiveGIMR.h:1478
NPEffectiveGIMR::CW
double CW
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1323
NPEffectiveGIMR::CLQ1_2112
double CLQ1_2112
Definition: NPEffectiveGIMR.h:1482
NPEffectiveGIMR::CHu_23r
double CHu_23r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1380
NPEffectiveGIMR::eWH2_HWud
double eWH2_HWud
Theoretical uncertainty in the (linear) new physics contribution from to WH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1564
NPEffectiveGIMR::eWH2_HWW2
double eWH2_HWW2
Theoretical uncertainty in the (linear) new physics contribution from to WH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1562
NPEffectiveGIMR::GammaHtautauRatio
double GammaHtautauRatio() const
The ratio of the in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:3267
NPEffectiveGIMR::GammaHmumuRatio
double GammaHmumuRatio() const
The ratio of the in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:3237
NPEffectiveGIMR::CRL_up
double CRL_up() const
Definition: NPEffectiveGIMR.cpp:3437
NPEffectiveGIMR::CHQ3_12i
double CHQ3_12i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1373
NPEffectiveGIMR::CHu_13i
double CHu_13i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1383
NPEffectiveGIMR::CLL_1331
double CLL_1331
Definition: NPEffectiveGIMR.h:1480
NPEffectiveGIMR::eVBF2_HZZ3
double eVBF2_HZZ3
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1522
NPEffectiveGIMR::CdH_23i
double CdH_23i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1440
NPEffectiveGIMR::GammaHbbRatio
double GammaHbbRatio() const
The ratio of the in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:3326
NPEffectiveGIMR::CLQ1_2223
double CLQ1_2223
Definition: NPEffectiveGIMR.h:1484
NPEffectiveGIMR::CfW_diag
gslpp::complex CfW_diag(const Particle f) const
The diagonal entry of the dimension-6 operator coefficient corresponding to particle f.
Definition: NPEffectiveGIMR.cpp:1634
NPEffectiveGIMR::CLu_1133
double CLu_1133
Definition: NPEffectiveGIMR.h:1507
NPEffectiveGIMR::BrHgagaRatio
virtual double BrHgagaRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:2914
NPEffectiveGIMR::CQe_1122
double CQe_1122
Definition: NPEffectiveGIMR.h:1514
NPEffectiveGIMR::CLL_1221
double CLL_1221
Definition: NPEffectiveGIMR.h:1479
NPEffectiveGIMR::CHud_33i
double CHud_33i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1405
NPEffectiveGIMR::Ced_3332
double Ced_3332
Definition: NPEffectiveGIMR.h:1501
NPEffectiveGIMR::CHe_33
double CHe_33
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1354
NPEffectiveGIMR::CuW_12i
double CuW_12i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1461
NPEffectiveGIMR::deltaGammaHgagaRatio1
double deltaGammaHgagaRatio1() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3213
NPEffectiveGIMR::eVBF2_HWW1
double eVBF2_HWW1
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1526
NPEffectiveGIMR::CLd_3311
double CLd_3311
Definition: NPEffectiveGIMR.h:1510
NPEffectiveGIMR::eWH2_Wud
double eWH2_Wud
Theoretical uncertainty in the (linear) new physics contribution from to WH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1565
NPEffectiveGIMR::CLQ1_3113
double CLQ1_3113
Definition: NPEffectiveGIMR.h:1483
NPEffectiveGIMR::CLL_1122
double CLL_1122
Definition: NPEffectiveGIMR.h:1479
NPEffectiveGIMR::CHD
double CHD
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1328
NPEffectiveGIMR::eWH2_HWW3
double eWH2_HWW3
Theoretical uncertainty in the (linear) new physics contribution from to WH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1563
NPEffectiveGIMR::CRR_up
double CRR_up() const
Definition: NPEffectiveGIMR.cpp:3473
NPEffectiveGIMR::CLQ3_2112
double CLQ3_2112
Definition: NPEffectiveGIMR.h:1487
NPEffectiveGIMR::CHu_12r
double CHu_12r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1377
QCD::UP
Definition: QCD.h:324
NPEffectiveGIMR::CeH_23i
double CeH_23i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1416
NPEffectiveGIMR::eVBF2_HWW3
double eVBF2_HWW3
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1528
StandardModel::GF
double GF
The Fermi constant in .
Definition: StandardModel.h:2555
NPEffectiveGIMR::eWH78_HWud
double eWH78_HWud
Theoretical uncertainty in the (linear) new physics contribution from to WH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1569
NPEffectiveGIMR::LambdaNP2
double LambdaNP2
The square of the new physics scale [GeV ].
Definition: NPEffectiveGIMR.h:1606
NPEffectiveGIMR::eZH2_ZdR
double eZH2_ZdR
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1584
NPEffectiveGIMR::CHL3_22
double CHL3_22
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1343
NPEffectiveGIMR::Cee_2211
double Cee_2211
Definition: NPEffectiveGIMR.h:1492
NPEffectiveGIMR::LEP2AFBtau
virtual double LEP2AFBtau(const double s) const
Definition: NPEffectiveGIMR.cpp:3601
NPEffectiveGIMR::CLQ3_1122
double CLQ3_1122
Definition: NPEffectiveGIMR.h:1487
NPEffectiveGIMR::muWH
virtual double muWH(const double sqrt_s) const
The ratio between the W-Higgs associated production cross-section in the current model and in the St...
Definition: NPEffectiveGIMR.cpp:2350
Model::addMissingModelParameter
void addMissingModelParameter(const std::string &missingParameterName)
Definition: Model.h:240
NPEffectiveGIMR::CRL_mu
double CRL_mu() const
Definition: NPEffectiveGIMR.cpp:3427
NPEffectiveGIMR::CLd_1123
double CLd_1123
Definition: NPEffectiveGIMR.h:1511
StandardModel::CheckParameters
virtual bool CheckParameters(const std::map< std::string, double > &DPars)
A method to check if all the mandatory parameters for StandardModel have been provided in model initi...
Definition: StandardModel.cpp:339
NPEffectiveGIMR::CHu_11
double CHu_11
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1376
NPEffectiveGIMR::deltaGR_Wffh
gslpp::complex deltaGR_Wffh(const Particle pbar, const Particle p) const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1898
NPEffectiveGIMR::CQe_3311
double CQe_3311
Definition: NPEffectiveGIMR.h:1515
StandardModel::alphaMz
double alphaMz() const
The electromagnetic coupling at the -mass scale, .
Definition: StandardModel.cpp:893
NPEffectiveGIMR::CeH_11i
double CeH_11i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1412
NPEffectiveGIMR::CLQ1_1122
double CLQ1_1122
Definition: NPEffectiveGIMR.h:1482
NPEffectiveGIMR::eVBF2_HZZ2
double eVBF2_HZZ2
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1521
NPEffectiveGIMR::Mw
virtual double Mw() const
The mass of the boson, .
Definition: NPEffectiveGIMR.cpp:1714
StandardModel::LEP2AFBcharm
virtual double LEP2AFBcharm(const double s) const
Definition: StandardModel.cpp:4882
NPEffectiveGIMR::v2_over_LambdaNP2
double v2_over_LambdaNP2
The ratio between the EW vev and the new physics scale, squared .
Definition: NPEffectiveGIMR.h:1607
NPEffectiveGIMR::deltaGammaTotalRatio1
virtual double deltaGammaTotalRatio1() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3027
NPEffectiveGIMR::CfG_diag
gslpp::complex CfG_diag(const Particle f) const
The diagonal entry of the dimension-6 operator coefficient corresponding to particle f.
Definition: NPEffectiveGIMR.cpp:1608
QCD::CHARM
Definition: QCD.h:326
NPEffectiveGIMR::eZH2_ZuR
double eZH2_ZuR
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1582
NPEffectiveGIMR::eVBF78_Hgg
double eVBF78_Hgg
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1549
NPEffectiveGIMR::CuW_33i
double CuW_33i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1465
NPEffectiveGIMR::CRR_bottom
double CRR_bottom() const
Definition: NPEffectiveGIMR.cpp:3493
NPEffectiveGIMR::BrHtautauRatio
virtual double BrHtautauRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:2954
NPEffectiveGIMR::CHud_13r
double CHud_13r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1396
NPbase::deltaGamma_Z
virtual double deltaGamma_Z() const
The new physics contribution to the total decay width of the boson, .
Definition: NPbase.cpp:176
NPEffectiveGIMR::CeH_22i
double CeH_22i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1415
NPEffectiveGIMR::CHL1_23i
double CHL1_23i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1339
NPEffectiveGIMR::eVBF78_HWW1
double eVBF78_HWW1
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1546
StandardModel::computeBrHtobb
double computeBrHtobb() const
The Br in the Standard Model.
Definition: StandardModel.h:2313
NPEffectiveGIMR::deltaG_Zff
gslpp::complex deltaG_Zff(const Particle p) const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1949
NPEffectiveGIMR::CuW_11i
double CuW_11i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1460
NPEffectiveGIMR::Ceu_1122
double Ceu_1122
Definition: NPEffectiveGIMR.h:1495
NPEffectiveGIMR::BrHccRatio
virtual double BrHccRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:2974
gslpp::complex
A class for defining operations on and functions of complex numbers.
Definition: gslpp_complex.h:35
StandardModel::LEP2AFBtau
virtual double LEP2AFBtau(const double s) const
Definition: StandardModel.cpp:5784
NPEffectiveGIMR::CHL3_23i
double CHL3_23i
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1348
NPEffectiveGIMR::muttH
virtual double muttH(const double sqrt_s) const
The ratio between the t-tbar-Higgs associated production cross-section in the current model and in t...
Definition: NPEffectiveGIMR.cpp:2700
NPEffectiveGIMR::CHud_12i
double CHud_12i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1401
StandardModel::mHl
double mHl
The Higgs mass in GeV.
Definition: StandardModel.h:2558
NPEffectiveGIMR::eZH78_HZdL
double eZH78_HZdL
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1592
gslpp::log
complex log(const complex &z)
Definition: gslpp_complex.cpp:342
NPEffectiveGIMR::CHQ1_33
double CHQ1_33
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1363
NPEffectiveGIMR::ettH2_Htt
double ettH2_Htt
Theoretical uncertainty in the (linear) new physics contribution from to ttH production at Tevatron ...
Definition: NPEffectiveGIMR.h:1599
StandardModel::LEP2Rbottom
virtual double LEP2Rbottom(const double s) const
Definition: StandardModel.cpp:6274
NPEffectiveGIMR::CHud_33r
double CHud_33r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1399
NPEffectiveGIMR::eVBF2_HZuR
double eVBF2_HZuR
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1531
NPEffectiveGIMR::LEP2sigmaMu
virtual double LEP2sigmaMu(const double s) const
Definition: NPEffectiveGIMR.cpp:3500
NPEffectiveGIMR::CHd_33
double CHd_33
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1390
NPEffectiveGIMR::CQe_1111
double CQe_1111
Definition: NPEffectiveGIMR.h:1513
NPEffectiveGIMR::CLQ3_3311
double CLQ3_3311
Definition: NPEffectiveGIMR.h:1488
QCD::ELECTRON
Definition: QCD.h:312
NPEffectiveGIMR::CLQ3_2211
double CLQ3_2211
Definition: NPEffectiveGIMR.h:1487
NPEffectiveGIMR::CHQ1_23r
double CHQ1_23r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1362
Particle::getIsospin
double getIsospin() const
A get method to access the particle isospin.
Definition: Particle.h:115
NPEffectiveGIMR::CuH_22i
double CuH_22i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1427
NPEffectiveGIMR::CLL_bottom
double CLL_bottom() const
Definition: NPEffectiveGIMR.cpp:3387
gslpp::complex::abs2
double abs2() const
Definition: gslpp_complex.cpp:86
NPEffectiveGIMR::CHL3_11
double CHL3_11
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1340
NPEffectiveGIMR::CHQ3_23i
double CHQ3_23i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1375
NPEffectiveGIMR::deltaGammaHtautauRatio1
double deltaGammaHtautauRatio1() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3283
NPEffectiveGIMR::Cee_1122
double Cee_1122
Definition: NPEffectiveGIMR.h:1492
NPEffectiveGIMR::delta_AA
double delta_AA
Combination of dimension 6 coefficients modifying the canonical field definition.
Definition: NPEffectiveGIMR.h:1613
NPEffectiveGIMR::Ceu_1133
double Ceu_1133
Definition: NPEffectiveGIMR.h:1496
NPEffectiveGIMR::CLQ1_1133
double CLQ1_1133
Definition: NPEffectiveGIMR.h:1483
NPEffectiveGIMR::Ced_3323
double Ced_3323
Definition: NPEffectiveGIMR.h:1500
NPEffectiveGIMR::Ceu_1111
double Ceu_1111
Definition: NPEffectiveGIMR.h:1494
NPEffectiveGIMR::CQe_1133
double CQe_1133
Definition: NPEffectiveGIMR.h:1515
NPEffectiveGIMR::CfB_diag
gslpp::complex CfB_diag(const Particle f) const
The diagonal entry of the dimension-6 operator coefficient corresponding to particle f.
Definition: NPEffectiveGIMR.cpp:1660
NPEffectiveGIMR::Ceu_2211
double Ceu_2211
Definition: NPEffectiveGIMR.h:1495
NPEffectiveGIMR::CHd_13i
double CHd_13i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1392
NPEffectiveGIMR::mueettH
virtual double mueettH(const double sqrt_s) const
The ratio between the production cross-section in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:2769
NPEffectiveGIMR::muggH
virtual double muggH(const double sqrt_s) const
The ratio between the gluon-gluon fusion Higgs production cross-section in the current model and in ...
Definition: NPEffectiveGIMR.cpp:1990
NPEffectiveGIMR::deltaG_hgg
virtual double deltaG_hgg() const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1815
NPEffectiveGIMR::CLR_mu
double CLR_mu() const
Definition: NPEffectiveGIMR.cpp:3392
NPEffectiveGIMR::CuW_13r
double CuW_13r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1456
NPEffectiveGIMR::deltaGammaHWWRatio1
double deltaGammaHWWRatio1() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3104
NPEffectiveGIMR::NNPEffectiveGIMRVars
static const int NNPEffectiveGIMRVars
The number of the model parameters in NPEffectiveGIMR.
Definition: NPEffectiveGIMR.h:616
NPEffectiveGIMR::CuB_12r
double CuB_12r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1467
NPEffectiveGIMR::ettH2_Hgg
double ettH2_Hgg
Theoretical uncertainty in the (linear) new physics contribution from to ttH production at Tevatron ...
Definition: NPEffectiveGIMR.h:1600
NPEffectiveGIMR::CuG_13r
double CuG_13r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1444
StandardModel::LEP2AFBbottom
virtual double LEP2AFBbottom(const double s) const
Definition: StandardModel.cpp:4470
NPEffectiveGIMR::CeH_22r
double CeH_22r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1409
QCD::mtpole
double mtpole
The pole mass of the top quark.
Definition: QCD.h:927
NPEffectiveGIMR::CLQ3_2232
double CLQ3_2232
Definition: NPEffectiveGIMR.h:1490
NPEffectiveGIMR::CfH_diag
gslpp::complex CfH_diag(const Particle f) const
The diagonal entry of the dimension-6 operator coefficient corresponding to particle f.
Definition: NPEffectiveGIMR.cpp:1582
StandardModel::setFlag
virtual bool setFlag(const std::string name, const bool value)
A method to set a flag of StandardModel.
Definition: StandardModel.cpp:404
Model::ModelParamMap
std::map< std::string, std::reference_wrapper< const double > > ModelParamMap
Definition: Model.h:270
NPEffectiveGIMR::CuH_13r
double CuH_13r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1420
NPEffectiveGIMR::deltaGA_f
virtual double deltaGA_f(const Particle p) const
New physics contribution to the neutral-current axial-vector coupling .
Definition: NPEffectiveGIMR.cpp:1746
NPEffectiveGIMR::CLQ1_2232
double CLQ1_2232
Definition: NPEffectiveGIMR.h:1485
NPEffectiveGIMR::CLL_charm
double CLL_charm() const
Definition: NPEffectiveGIMR.cpp:3377
NPEffectiveGIMR::CuG_11r
double CuG_11r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1442
NPEffectiveGIMR::deltaGR_Zffh
double deltaGR_Zffh(const Particle p) const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1915
NPEffectiveGIMR::CuW_22i
double CuW_22i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1463
NPEffectiveGIMR::LEP2AFBcharm
virtual double LEP2AFBcharm(const double s) const
Definition: NPEffectiveGIMR.cpp:3619
NPEffectiveGIMR::eVBF2_ZuR
double eVBF2_ZuR
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1536
NPEffectiveGIMR::cW2_tree
double cW2_tree
The sqaure of the tree level values for the cosine of the weak angle.
Definition: NPEffectiveGIMR.h:1610
NPEffectiveGIMR::CHG
double CHG
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1324
NPEffectiveGIMR::muZH
virtual double muZH(const double sqrt_s) const
The ratio between the Z-Higgs associated production cross-section in the current model and in the St...
Definition: NPEffectiveGIMR.cpp:2425
NPEffectiveGIMR::NNPEffectiveGIMRVars_LFU_QFU
static const int NNPEffectiveGIMRVars_LFU_QFU
The number of the model parameters in NPEffectiveGIMR with lepton and quark flavour universalities.
Definition: NPEffectiveGIMR.h:628
NPEffectiveGIMR::CLQ1_3323
double CLQ1_3323
Definition: NPEffectiveGIMR.h:1484
NPEffectiveGIMR::CHud_23r
double CHud_23r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1398
NPEffectiveGIMR::deltaGL_f
double deltaGL_f(const Particle p) const
New physics contribution to the neutral-current left-handed coupling .
Definition: NPEffectiveGIMR.cpp:1751
NPEffectiveGIMR::BrHZZRatio
virtual double BrHZZRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:2874
NPEffectiveGIMR::CRR_charm
double CRR_charm() const
Definition: NPEffectiveGIMR.cpp:3483
NPEffectiveGIMR::CLQ1_1221
double CLQ1_1221
Definition: NPEffectiveGIMR.h:1482
NPEffectiveGIMR::BrHmumuRatio
virtual double BrHmumuRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:2934
NPEffectiveGIMR::CRR_strange
double CRR_strange() const
Definition: NPEffectiveGIMR.cpp:3488
NPEffectiveGIMR::CLL_3311
double CLL_3311
Definition: NPEffectiveGIMR.h:1480
NPbase::trueSM
StandardModel trueSM
Definition: NPbase.h:2787
NPEffectiveGIMR::CuG_23i
double CuG_23i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1452
NPEffectiveGIMR::Ced_1132
double Ced_1132
Definition: NPEffectiveGIMR.h:1501
StandardModel::computeBrHtoZZ
double computeBrHtoZZ() const
The Br in the Standard Model.
Definition: StandardModel.h:2222
NPEffectiveGIMR::CHL3_33
double CHL3_33
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1345
NPEffectiveGIMR::CHe_12r
double CHe_12r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1350
NPEffectiveGIMR::CuH_12i
double CuH_12i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1425
NPEffectiveGIMR::CQe_2311
double CQe_2311
Definition: NPEffectiveGIMR.h:1516
NPEffectiveGIMR::CHe_13r
double CHe_13r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1351
NPEffectiveGIMR::eZH78_HZA1
double eZH78_HZA1
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1588
NPEffectiveGIMR::CuG_22r
double CuG_22r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1445
NPEffectiveGIMR::CLR_tau
double CLR_tau() const
Definition: NPEffectiveGIMR.cpp:3397
NPEffectiveGIMR::CG
double CG
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1322
NPEffectiveGIMR::CLL_down
double CLL_down() const
Definition: NPEffectiveGIMR.cpp:3372
NPEffectiveGIMR::deltaGammaHgagaRatio2
double deltaGammaHgagaRatio2() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3224
NPEffectiveGIMR::eVBF78_Wud
double eVBF78_Wud
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1559
NPEffectiveGIMR::eVBF78_HZZ3
double eVBF78_HZZ3
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1542
NPEffectiveGIMR::CLR_charm
double CLR_charm() const
Definition: NPEffectiveGIMR.cpp:3412
Particle::getMass
const double & getMass() const
A get method to access the particle mass.
Definition: Particle.h:61
NPEffectiveGIMR::CHd_22
double CHd_22
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1388
NPEffectiveGIMR::deltaGammaHccRatio2
double deltaGammaHccRatio2() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3318
StandardModel::AlsMz
double AlsMz
The strong coupling constant at the Z-boson mass, .
Definition: StandardModel.h:2553
NPbase
The auxiliary base model class for other model classes.
Definition: NPbase.h:66
NPEffectiveGIMR::eZH78_ZuL
double eZH78_ZuL
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1594
NPEffectiveGIMR::deltaG1_hWW
virtual double deltaG1_hWW() const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1820
NPEffectiveGIMR::deltag3G
double deltag3G() const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1963
NPEffectiveGIMR::CHu_22
double CHu_22
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1379
LEP2GIMR::sigmaFminusB_q_LEP2_GIMR
double sigmaFminusB_q_LEP2_GIMR(const QCD::quark q, const double s, const double GIMRParam_i[]) const
Definition: LEP2GIMR.cpp:155
NPEffectiveGIMR::obliqueT
virtual double obliqueT() const
The oblique parameter .
Definition: NPEffectiveGIMR.cpp:1699
NPEffectiveGIMR::eVBF78_HZdL
double eVBF78_HZdL
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1552
NPEffectiveGIMR::eVBF78_HZuR
double eVBF78_HZuR
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1551
NPEffectiveGIMR::CHud_13i
double CHud_13i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1402
NPEffectiveGIMR::NPEffectiveGIMRVars
static const std::string NPEffectiveGIMRVars[NNPEffectiveGIMRVars]
A string array containing the labels of the model parameters in NPEffectiveGIMR.
Definition: NPEffectiveGIMR.h:622
NPEffectiveGIMR::CuG_33r
double CuG_33r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1447
NPEffectiveGIMR::computeGammaTotalRatio
virtual double computeGammaTotalRatio() const
The ratio of the in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:3014
NPEffectiveGIMR::CHL3_23r
double CHL3_23r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1344
NPEffectiveGIMR::CuB_22r
double CuB_22r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1469
NPEffectiveGIMR::CLu_3311
double CLu_3311
Definition: NPEffectiveGIMR.h:1507
NPEffectiveGIMR::eZH2_HZZ1
double eZH2_HZZ1
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1572
QCD::TOP
Definition: QCD.h:328
NPEffectiveGIMR::eVBF2_HZuL
double eVBF2_HZuL
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1530
NPEffectiveGIMR::CHbox
double CHbox
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1329
NPEffectiveGIMR::CLe_1122
double CLe_1122
Definition: NPEffectiveGIMR.h:1503
NPEffectiveGIMR::eZH78_HZA2
double eZH78_HZA2
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1589
NPEffectiveGIMR::eVBF2_HWud
double eVBF2_HWud
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1534
NPEffectiveGIMR::eZH78_ZdR
double eZH78_ZdR
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1597
NPEffectiveGIMR::CuH_13i
double CuH_13i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1426
NPEffectiveGIMR::CdH_11r
double CdH_11r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1430
NPEffectiveGIMR::CHF1_diag
double CHF1_diag(const Particle F) const
The diagonal entry of the dimension-6 operator coefficient corresponding to particle F.
Definition: NPEffectiveGIMR.cpp:1505
NPEffectiveGIMR::CLd_2232
double CLd_2232
Definition: NPEffectiveGIMR.h:1512
LEP2GIMR::sigma_q_LEP2_GIMR
double sigma_q_LEP2_GIMR(const QCD::quark q, const double s, const double GIMRParam_i[]) const
Definition: LEP2GIMR.cpp:111
NPEffectiveGIMR::eVBF78_ZdR
double eVBF78_ZdR
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1558
NPEffectiveGIMR::CHud_diag
gslpp::complex CHud_diag(const Particle u) const
The diagonal entry of the dimension-6 operator coefficient corresponding to particle f.
Definition: NPEffectiveGIMR.cpp:1567
NPEffectiveGIMR::CLQ3_2223
double CLQ3_2223
Definition: NPEffectiveGIMR.h:1489
gslpp::pow
complex pow(const complex &z1, const complex &z2)
Definition: gslpp_complex.cpp:395
NPEffectiveGIMR::CLQ1_2211
double CLQ1_2211
Definition: NPEffectiveGIMR.h:1482
NPEffectiveGIMR::CdH_33i
double CdH_33i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1441
NPEffectiveGIMR::CLQ3_1133
double CLQ3_1133
Definition: NPEffectiveGIMR.h:1488
NPEffectiveGIMR::eVBF78_HWW3
double eVBF78_HWW3
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1548
NPEffectiveGIMR::CdH_22r
double CdH_22r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1433
NPEffectiveGIMR::FlagQuarkUniversal
const bool FlagQuarkUniversal
An internal boolean flag that is true if assuming quark flavour universality.
Definition: NPEffectiveGIMR.h:1691
Model::raiseMissingModelParameterCount
void raiseMissingModelParameterCount()
Definition: Model.h:250
NPEffectiveGIMR::deltaG_hAff
gslpp::complex deltaG_hAff(const Particle p) const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1935
gslpp::sqrt
complex sqrt(const complex &z)
Definition: gslpp_complex.cpp:385
NPEffectiveGIMR::CLL_3113
double CLL_3113
Definition: NPEffectiveGIMR.h:1480
NPEffectiveGIMR::CuW_33r
double CuW_33r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1459
NPEffectiveGIMR::CLQ3_1123
double CLQ3_1123
Definition: NPEffectiveGIMR.h:1489
NPEffectiveGIMR::CHud_22i
double CHud_22i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1403
gslpp::complex::i
static const complex & i()
Definition: gslpp_complex.cpp:154
NPEffectiveGIMR::deltaG2_hWW
virtual double deltaG2_hWW() const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1825
NPEffectiveGIMR::muggHpttH
virtual double muggHpttH(const double sqrt_s) const
The ratio between the sum of gluon-gluon fusion and t-tbar-Higgs associated production cross-section...
Definition: NPEffectiveGIMR.cpp:2755
NPEffectiveGIMR::GammaW
virtual double GammaW() const
The total width of the boson, .
Definition: NPEffectiveGIMR.cpp:1725
NPEffectiveGIMR::CHQ1_22
double CHQ1_22
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1361
NPEffectiveGIMR::eZH78_ZuR
double eZH78_ZuR
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1595
StandardModel::computeBrHtoZga
double computeBrHtoZga() const
The Br in the Standard Model.
Definition: StandardModel.h:2244
NPEffectiveGIMR::eVBF2_Hgg
double eVBF2_Hgg
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1529
NPEffectiveGIMR::Cee_3311
double Cee_3311
Definition: NPEffectiveGIMR.h:1493
NPEffectiveGIMR::CHB
double CHB
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1326
Particle::getCharge
double getCharge() const
A get method to access the particle charge.
Definition: Particle.h:97
NPEffectiveGIMR::eVBF78_HZuL
double eVBF78_HZuL
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1550
NPEffectiveGIMR::NPEffectiveGIMRVars_LFU_QFU
static const std::string NPEffectiveGIMRVars_LFU_QFU[NNPEffectiveGIMRVars_LFU_QFU]
A string array containing the labels of the model parameters in NPEffectiveGIMR with lepton and quark...
Definition: NPEffectiveGIMR.h:634
NPEffectiveGIMR::DeltaGF
virtual double DeltaGF() const
New physics contribution to the Fermi constant.
Definition: NPEffectiveGIMR.cpp:1689
StandardModel::computeBrHtogaga
double computeBrHtogaga() const
The Br in the Standard Model.
Definition: StandardModel.h:2256
NPEffectiveGIMR::eVBF78_HZA1
double eVBF78_HZA1
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1543
StandardModel::computeSigmaggH
double computeSigmaggH(const double sqrt_s) const
The ggH cross section in the Standard Model.
Definition: StandardModel.h:1897
NPEffectiveGIMR::CLe_1133
double CLe_1133
Definition: NPEffectiveGIMR.h:1504
NPEffectiveGIMR::CHL3_13r
double CHL3_13r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1342
NPEffectiveGIMR::CLR_strange
double CLR_strange() const
Definition: NPEffectiveGIMR.cpp:3417
NPEffectiveGIMR::deltaG1_hZA
virtual double deltaG1_hZA() const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1860
NPEffectiveGIMR::CHd_11
double CHd_11
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1385
NPEffectiveGIMR::eZH78_ZdL
double eZH78_ZdL
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1596
NPEffectiveGIMR::CLR_down
double CLR_down() const
Definition: NPEffectiveGIMR.cpp:3407
NPEffectiveGIMR::CdH_12r
double CdH_12r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1431
NPEffectiveGIMR::CHQ3_13r
double CHQ3_13r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1369
NPEffectiveGIMR::CLL_1133
double CLL_1133
Definition: NPEffectiveGIMR.h:1480
NPEffectiveGIMR::CHL1_22
double CHL1_22
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1334
StandardModel::LEP2sigmaCharm
virtual double LEP2sigmaCharm(const double s) const
Definition: StandardModel.cpp:2924
NPEffectiveGIMR::Ced_2211
double Ced_2211
Definition: NPEffectiveGIMR.h:1498
NPEffectiveGIMR::CuW_12r
double CuW_12r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1455
NPEffectiveGIMR::eZH2_ZdL
double eZH2_ZdL
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1583
NPEffectiveGIMR::CHWB
double CHWB
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1327
NPEffectiveGIMR::CdH_13i
double CdH_13i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1438
NPEffectiveGIMR::CRL_down
double CRL_down() const
Definition: NPEffectiveGIMR.cpp:3442
NPEffectiveGIMR.h
StandardModel::computeSigmattH
double computeSigmattH(const double sqrt_s) const
The ttH production cross section in the Standard Model.
Definition: StandardModel.h:2171
NPEffectiveGIMR::eVBF78_ZuL
double eVBF78_ZuL
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1555
NPEffectiveGIMR::CHL1_12r
double CHL1_12r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1332
NPEffectiveGIMR::eZH2_HZuR
double eZH2_HZuR
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1578
NPEffectiveGIMR::CLQ1_3311
double CLQ1_3311
Definition: NPEffectiveGIMR.h:1483
NPEffectiveGIMR::LEP2Rbottom
virtual double LEP2Rbottom(const double s) const
Definition: NPEffectiveGIMR.cpp:3709
NPEffectiveGIMR::muVH
virtual double muVH(const double sqrt_s) const
The ratio between the WH+ZH associated production cross-section in the current model and in the Stan...
Definition: NPEffectiveGIMR.cpp:2672
NPEffectiveGIMR::setParameter
virtual void setParameter(const std::string name, const double &value)
A method to set the value of a parameter of the model.
Definition: NPEffectiveGIMR.cpp:558
NPEffectiveGIMR::CHL1_33
double CHL1_33
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1336
NPEffectiveGIMR::eZH78_HZuR
double eZH78_HZuR
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1591
NPEffectiveGIMR::GammaHggRatio
double GammaHggRatio() const
The ratio of the in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:3053
NPEffectiveGIMR::Ced_2223
double Ced_2223
Definition: NPEffectiveGIMR.h:1500
NPEffectiveGIMR::muVBF
virtual double muVBF(const double sqrt_s) const
The ratio between the vector-boson fusion Higgs production cross-section in the current model and in...
Definition: NPEffectiveGIMR.cpp:2026
NPEffectiveGIMR::eWH78_HWW1
double eWH78_HWW1
Theoretical uncertainty in the (linear) new physics contribution from to WH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1566
NPEffectiveGIMR::CHQ3_12r
double CHQ3_12r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1368
NPEffectiveGIMR::CeH_33r
double CeH_33r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1411
NPEffectiveGIMR::f_triangle
gslpp::complex f_triangle(const double tau) const
Loop function entering in the calculation of the effective and couplings.
Definition: NPEffectiveGIMR.cpp:1973
QCD::quark
quark
An enum type for quarks.
Definition: QCD.h:323
NPEffectiveGIMR::deltaG_hGff
gslpp::complex deltaG_hGff(const Particle p) const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1921
NPEffectiveGIMR::eVBF2_ZuL
double eVBF2_ZuL
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1535
NPEffectiveGIMR::Cee_1133
double Cee_1133
Definition: NPEffectiveGIMR.h:1493
NPEffectiveGIMR::CuH_11r
double CuH_11r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1418
StandardModel::GammaW
virtual double GammaW(const Particle fi, const Particle fj) const
A partial decay width of the boson decay into a SM fermion pair.
Definition: StandardModel.cpp:1166
NPEffectiveGIMR::CHud_11r
double CHud_11r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1394
NPEffectiveGIMR::CLQ1_1132
double CLQ1_1132
Definition: NPEffectiveGIMR.h:1485
NPEffectiveGIMR::CRR_tau
double CRR_tau() const
Definition: NPEffectiveGIMR.cpp:3467
NPEffectiveGIMR::CHQ1_12r
double CHQ1_12r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1359
NPEffectiveGIMR::CuH_33r
double CuH_33r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1423
NPEffectiveGIMR::CdH_23r
double CdH_23r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1434
NPEffectiveGIMR::eZH2_HZdL
double eZH2_HZdL
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1579
NPEffectiveGIMR::FlagMwInput
bool FlagMwInput
A boolean flag that is true if the W mass is taken as an input parameter. (Warning: The W width is no...
Definition: NPEffectiveGIMR.h:1678
NPEffectiveGIMR::CHQ3_33
double CHQ3_33
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1372
NPEffectiveGIMR::CHW
double CHW
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1325
NPEffectiveGIMR::GammaHWWRatio
double GammaHWWRatio() const
The ratio of the in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:3088
NPEffectiveGIMR::eVBF78_HAA
double eVBF78_HAA
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1545
NPEffectiveGIMR::deltaGammaTotalRatio2
virtual double deltaGammaTotalRatio2() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3040
NPEffectiveGIMR::CHQ1_23i
double CHQ1_23i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1366
NPEffectiveGIMR::CdH_11i
double CdH_11i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1436
NPEffectiveGIMR::CHu_12i
double CHu_12i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1382
NPEffectiveGIMR::CLL_tau
double CLL_tau() const
Definition: NPEffectiveGIMR.cpp:3362
NPEffectiveGIMR::CdH_12i
double CdH_12i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1437
NPEffectiveGIMR::Ced_3311
double Ced_3311
Definition: NPEffectiveGIMR.h:1499
NPEffectiveGIMR::FlagLeptonUniversal
const bool FlagLeptonUniversal
An internal boolean flag that is true if assuming lepton flavour universality.
Definition: NPEffectiveGIMR.h:1685
NPEffectiveGIMR::eVBF78_HZdR
double eVBF78_HZdR
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1553
NPEffectiveGIMR::CRL_tau
double CRL_tau() const
Definition: NPEffectiveGIMR.cpp:3432
NPEffectiveGIMR::LEP2sigmaTau
virtual double LEP2sigmaTau(const double s) const
Definition: NPEffectiveGIMR.cpp:3516
NPEffectiveGIMR::CLL_2112
double CLL_2112
Definition: NPEffectiveGIMR.h:1479
NPEffectiveGIMR::mueeWBF
virtual double mueeWBF(const double sqrt_s) const
The ratio between the production cross-section in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:2219
NPEffectiveGIMR::CRR_down
double CRR_down() const
Definition: NPEffectiveGIMR.cpp:3478
NPEffectiveGIMR::eVBF78_HZA2
double eVBF78_HZA2
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1544
NPEffectiveGIMR::eZH2_HZuL
double eZH2_HZuL
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1577
NPEffectiveGIMR::CLQ1_1111
double CLQ1_1111
Definition: NPEffectiveGIMR.h:1481
NPEffectiveGIMR::eVBF2_ZdL
double eVBF2_ZdL
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1537
NPEffectiveGIMR::CuG_22i
double CuG_22i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1451
NPEffectiveGIMR::eVBF78_HWW2
double eVBF78_HWW2
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1547
NPEffectiveGIMR::CuG_12i
double CuG_12i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1449
NPEffectiveGIMR::CHQ3_23r
double CHQ3_23r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1371
NPEffectiveGIMR::Ced_1133
double Ced_1133
Definition: NPEffectiveGIMR.h:1499
NPEffectiveGIMR::CHe_22
double CHe_22
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1352
NPEffectiveGIMR::CHL3_12i
double CHL3_12i
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1346
NPEffectiveGIMR::eVBF2_HWW2
double eVBF2_HWW2
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1527
NPEffectiveGIMR::CuW_22r
double CuW_22r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1457
NPEffectiveGIMR::deltaGR_Wff
virtual gslpp::complex deltaGR_Wff(const Particle pbar, const Particle p) const
New physics contribution to the charged current coupling .
Definition: NPEffectiveGIMR.cpp:1806
NPEffectiveGIMR::CLL_strange
double CLL_strange() const
Definition: NPEffectiveGIMR.cpp:3382
NPEffectiveGIMR::CHQ3_22
double CHQ3_22
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1370
NPEffectiveGIMR::CdH_22i
double CdH_22i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1439
StandardModel::computeSigmaZH
double computeSigmaZH(const double sqrt_s) const
The ZH production cross section in the Standard Model.
Definition: StandardModel.h:2135
NPEffectiveGIMR::eVBF78_ZdL
double eVBF78_ZdL
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1557
NPEffectiveGIMR::LEP2sigmaHadron
virtual double LEP2sigmaHadron(const double s) const
Definition: NPEffectiveGIMR.cpp:3532
NPEffectiveGIMR::CHud_12r
double CHud_12r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1395
NPEffectiveGIMR::deltaGammaHmumuRatio1
double deltaGammaHmumuRatio1() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3253
StandardModel::computeBrHtocc
double computeBrHtocc() const
The Br in the Standard Model.
Definition: StandardModel.h:2290
NPEffectiveGIMR::CHQ1_13i
double CHQ1_13i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1365
NPEffectiveGIMR::CQe_3211
double CQe_3211
Definition: NPEffectiveGIMR.h:1517
NPEffectiveGIMR::CHu_23i
double CHu_23i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1384
NPEffectiveGIMR::obliqueU
virtual double obliqueU() const
The oblique parameter .
Definition: NPEffectiveGIMR.cpp:1704
NPEffectiveGIMR::CLQ3_3332
double CLQ3_3332
Definition: NPEffectiveGIMR.h:1490
NPEffectiveGIMR::BrHZgaRatio
virtual double BrHZgaRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:2894
NPEffectiveGIMR::deltaG3_hWW
virtual double deltaG3_hWW() const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1830
NPEffectiveGIMR::PostUpdate
virtual bool PostUpdate()
The post-update method for NPEffectiveGIMR.
Definition: NPEffectiveGIMR.cpp:535
NPEffectiveGIMR::CuH_12r
double CuH_12r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1419
NPEffectiveGIMR::CHd_12i
double CHd_12i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1391
NPEffectiveGIMR::CHL3_12r
double CHL3_12r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1341
StandardModel::Mw_tree
virtual double Mw_tree() const
The tree-level mass of the boson, .
Definition: StandardModel.cpp:951
NPEffectiveGIMR::CHd_23r
double CHd_23r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1389
NPEffectiveGIMR::CHd_12r
double CHd_12r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1386
NPEffectiveGIMR::Ced_1123
double Ced_1123
Definition: NPEffectiveGIMR.h:1500
NPEffectiveGIMR::deltaGammaHZgaRatio1
double deltaGammaHZgaRatio1() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3173
StandardModel::computeBrHtomumu
double computeBrHtomumu() const
The Br in the Standard Model.
Definition: StandardModel.h:2267
StandardModel::computeSigmaVBF
double computeSigmaVBF(const double sqrt_s) const
The VBF cross section in the Standard Model.
Definition: StandardModel.h:2003
NPEffectiveGIMR::GammaHccRatio
double GammaHccRatio() const
The ratio of the in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:3297
StandardModel::computeBrHtoWW
double computeBrHtoWW() const
The Br in the Standard Model.
Definition: StandardModel.h:2210
NPEffectiveGIMR::deltaG_Aff
gslpp::complex deltaG_Aff(const Particle p) const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1956
NPEffectiveGIMR::deltaG_hAA
virtual double deltaG_hAA() const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1870
NPEffectiveGIMR::deltaGL_Wff
virtual gslpp::complex deltaGL_Wff(const Particle pbar, const Particle p) const
New physics contribution to the charged current coupling .
Definition: NPEffectiveGIMR.cpp:1789
NPEffectiveGIMR::CLQ1_3332
double CLQ1_3332
Definition: NPEffectiveGIMR.h:1485
NPEffectiveGIMR::CuW_23r
double CuW_23r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1458
NPEffectiveGIMR::CQe_3233
double CQe_3233
Definition: NPEffectiveGIMR.h:1517
StandardModel::LEP2sigmaHadron
virtual double LEP2sigmaHadron(const double s) const
Definition: StandardModel.cpp:3346
NPEffectiveGIMR::CuH_23r
double CuH_23r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1422
QCD::STRANGE
Definition: QCD.h:327
NPEffectiveGIMR::CLd_1122
double CLd_1122
Definition: NPEffectiveGIMR.h:1509
NPEffectiveGIMR::ettH78_Hgg
double ettH78_Hgg
Theoretical uncertainty in the (linear) new physics contribution from to ttH production at the LHC (...
Definition: NPEffectiveGIMR.h:1602
NPEffectiveGIMR::deltaGammaHggRatio2
double deltaGammaHggRatio2() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3077
NPEffectiveGIMR::CuB_22i
double CuB_22i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1475
NPEffectiveGIMR::CLL_2211
double CLL_2211
Definition: NPEffectiveGIMR.h:1479
NPEffectiveGIMR::FlagQuadraticTerms
bool FlagQuadraticTerms
A boolean flag that is true if the quadratic terms in cross sections and widths are switched on.
Definition: NPEffectiveGIMR.h:1679
NPEffectiveGIMR::CuG_33i
double CuG_33i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1453
NPEffectiveGIMR::Ceu_3311
double Ceu_3311
Definition: NPEffectiveGIMR.h:1496
NPEffectiveGIMR::eZH78_HZZ2
double eZH78_HZZ2
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1586
NPEffectiveGIMR::CLQ3_3113
double CLQ3_3113
Definition: NPEffectiveGIMR.h:1488
NPEffectiveGIMR::deltaG2_hZA
virtual double deltaG2_hZA() const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1865
gslpp::complex::real
const double & real() const
Definition: gslpp_complex.cpp:53
NPEffectiveGIMR::myLEP2GIMR
const LEP2GIMR myLEP2GIMR
Definition: NPEffectiveGIMR.h:1673
NPEffectiveGIMR::CLQ3_1331
double CLQ3_1331
Definition: NPEffectiveGIMR.h:1488
NPEffectiveGIMR::CQe_2322
double CQe_2322
Definition: NPEffectiveGIMR.h:1516
NPEffectiveGIMR::deltaGL_Wffh
gslpp::complex deltaGL_Wffh(const Particle pbar, const Particle p) const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1889
NPEffectiveGIMR::GammaHgagaRatio
double GammaHgagaRatio() const
The ratio of the in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:3197
NPEffectiveGIMR::CLd_3323
double CLd_3323
Definition: NPEffectiveGIMR.h:1511
StandardModel::leptons
Particle leptons[6]
An array of Particle objects for the leptons.
Definition: StandardModel.h:2540
NPEffectiveGIMR::CHe_13i
double CHe_13i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1356
Particle::getIndex
int getIndex() const
Definition: Particle.h:160
NPEffectiveGIMR::eVBF78_HZZ2
double eVBF78_HZZ2
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1541
NPEffectiveGIMR::CLR_bottom
double CLR_bottom() const
Definition: NPEffectiveGIMR.cpp:3422
NPEffectiveGIMR::BrHggRatio
virtual double BrHggRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:2834
NPEffectiveGIMR::CLd_1133
double CLd_1133
Definition: NPEffectiveGIMR.h:1510
NPEffectiveGIMR::eZH2_ZuL
double eZH2_ZuL
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1581
NPEffectiveGIMR::CeH_13r
double CeH_13r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1408
StandardModel::LEP2AFBmu
virtual double LEP2AFBmu(const double s) const
Definition: StandardModel.cpp:5294
NPEffectiveGIMR::BrHWWRatio
virtual double BrHWWRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:2854
NPEffectiveGIMR::CuB_33i
double CuB_33i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1477
NPEffectiveGIMR::CLu_1122
double CLu_1122
Definition: NPEffectiveGIMR.h:1506
NPEffectiveGIMR::Lambda_NP
double Lambda_NP
The new physics scale [GeV].
Definition: NPEffectiveGIMR.h:1518
NPbase::PostUpdate
virtual bool PostUpdate()
The postupdate method for NPbase.
Definition: NPbase.cpp:23
StandardModel::computeBrHtogg
double computeBrHtogg() const
The Br in the Standard Model.
Definition: StandardModel.h:2199
Model::name
std::string name
The name of the model.
Definition: Model.h:275
NPEffectiveGIMR::CLQ3_1132
double CLQ3_1132
Definition: NPEffectiveGIMR.h:1490
StandardModel::Mz
double Mz
The mass of the boson in GeV.
Definition: StandardModel.h:2554
Model::setModelLinearized
void setModelLinearized(bool linearized=true)
Definition: Model.h:231
NPEffectiveGIMR::CRL_strange
double CRL_strange() const
Definition: NPEffectiveGIMR.cpp:3452
NPEffectiveGIMR::CHd_23i
double CHd_23i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1393
NPEffectiveGIMR::CHd_13r
double CHd_13r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1387
NPEffectiveGIMR::CuB_11r
double CuB_11r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1466
QCD::Nc
double Nc
The number of colours.
Definition: QCD.h:932
NPEffectiveGIMR::CuH_23i
double CuH_23i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1428
NPEffectiveGIMR::obliqueS
virtual double obliqueS() const
The oblique parameter .
Definition: NPEffectiveGIMR.cpp:1694
NPEffectiveGIMR::CQe_2211
double CQe_2211
Definition: NPEffectiveGIMR.h:1514
NPEffectiveGIMR::CQe_2333
double CQe_2333
Definition: NPEffectiveGIMR.h:1516
NPEffectiveGIMR::deltaGL_Zffh
double deltaGL_Zffh(const Particle p) const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1907
NPEffectiveGIMR::CLd_2223
double CLd_2223
Definition: NPEffectiveGIMR.h:1511
NPEffectiveGIMR::Cee_1111
double Cee_1111
Definition: NPEffectiveGIMR.h:1491
NPEffectiveGIMR::mueeZH
virtual double mueeZH(const double sqrt_s) const
The ratio between the associated production cross-section in the current model and in the Standard ...
Definition: NPEffectiveGIMR.cpp:2567
NPEffectiveGIMR::deltaG3_hZZ
virtual double deltaG3_hZZ() const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1853
StandardModel::Mw
virtual double Mw() const
The SM prediction for the -boson mass in the on-shell scheme, .
Definition: StandardModel.cpp:970
NPEffectiveGIMR::eVBF2_HZA1
double eVBF2_HZA1
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1523
NPEffectiveGIMR::CuB_12i
double CuB_12i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1473
NPEffectiveGIMR::CuG_12r
double CuG_12r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1443
NPEffectiveGIMR::CHe_11
double CHe_11
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1349
NPEffectiveGIMR::CHL3_13i
double CHL3_13i
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1347
NPEffectiveGIMR::CRL_charm
double CRL_charm() const
Definition: NPEffectiveGIMR.cpp:3447
StandardModel::LEP2sigmaBottom
virtual double LEP2sigmaBottom(const double s) const
Definition: StandardModel.cpp:3135
NPEffectiveGIMR::eVBF2_ZdR
double eVBF2_ZdR
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at Tevatron ...
Definition: NPEffectiveGIMR.h:1538
NPEffectiveGIMR::eZH78_HZdR
double eZH78_HZdR
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1593
NPEffectiveGIMR::CuH_33i
double CuH_33i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1429
NPEffectiveGIMR::eZH2_HZZ2
double eZH2_HZZ2
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1573
NPEffectiveGIMR::sW2_tree
double sW2_tree
The sqaure of the tree level values for the sine of the weak angle.
Definition: NPEffectiveGIMR.h:1611
NPEffectiveGIMR::eZH78_HZZ1
double eZH78_HZZ1
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1585
NPEffectiveGIMR::CeH_12r
double CeH_12r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1407
StandardModel::LEP2sigmaMu
virtual double LEP2sigmaMu(const double s) const
Definition: StandardModel.cpp:2425
NPEffectiveGIMR::deltaGammaHZZRatio1
double deltaGammaHZZRatio1() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3138
StandardModel::LEP2Rcharm
virtual double LEP2Rcharm(const double s) const
Definition: StandardModel.cpp:6286
NPEffectiveGIMR::CQe_3222
double CQe_3222
Definition: NPEffectiveGIMR.h:1517
NPEffectiveGIMR::CHe_23r
double CHe_23r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1353
NPEffectiveGIMR::eZH2_HZdR
double eZH2_HZdR
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1580
NPEffectiveGIMR::CeH_33i
double CeH_33i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1417
NPEffectiveGIMR::CLQ3_3323
double CLQ3_3323
Definition: NPEffectiveGIMR.h:1489
NPEffectiveGIMR::CdH_33r
double CdH_33r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1435
NPEffectiveGIMR::Ced_2232
double Ced_2232
Definition: NPEffectiveGIMR.h:1501
NPEffectiveGIMR::CHu_13r
double CHu_13r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1378
NPEffectiveGIMR::CuB_33r
double CuB_33r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1471
NPEffectiveGIMR::NPEffectiveGIMR
NPEffectiveGIMR(const bool FlagLeptonUniversal_in=false, const bool FlagQuarkUniversal_in=false)
Constructor.
Definition: NPEffectiveGIMR.cpp:106
NPEffectiveGIMR::eWH78_HWW2
double eWH78_HWW2
Theoretical uncertainty in the (linear) new physics contribution from to WH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1567
NPEffectiveGIMR::CuG_23r
double CuG_23r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1446
QCD::DOWN
Definition: QCD.h:325
NPEffectiveGIMR::CHQ3_11
double CHQ3_11
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1367
NPEffectiveGIMR::CLL_mu
double CLL_mu() const
Definition: NPEffectiveGIMR.cpp:3357
NPEffectiveGIMR::CLd_3332
double CLd_3332
Definition: NPEffectiveGIMR.h:1512
NPEffectiveGIMR::CLe_2211
double CLe_2211
Definition: NPEffectiveGIMR.h:1503
NPEffectiveGIMR::Ced_1122
double Ced_1122
Definition: NPEffectiveGIMR.h:1498
NPEffectiveGIMR::deltaGammaHtautauRatio2
double deltaGammaHtautauRatio2() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3289
NPEffectiveGIMR::CHud_23i
double CHud_23i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1404
NPEffectiveGIMR::CuG_13i
double CuG_13i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1450
NPEffectiveGIMR::cW_tree
double cW_tree
The tree level values for the cosine of the weak angle.
Definition: NPEffectiveGIMR.h:1608
NPEffectiveGIMR::CHud_11i
double CHud_11i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1400
NPEffectiveGIMR::CeH_13i
double CeH_13i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1414
NPEffectiveGIMR::CuB_13i
double CuB_13i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1474
NPEffectiveGIMR::CRR_mu
double CRR_mu() const
Definition: NPEffectiveGIMR.cpp:3462
NPEffectiveGIMR::CLu_1111
double CLu_1111
Definition: NPEffectiveGIMR.h:1505
NPEffectiveGIMR::CuW_11r
double CuW_11r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1454
NPEffectiveGIMR::deltaGammaHbbRatio1
double deltaGammaHbbRatio1() const
The new physics contribution to the ratio of the in the current model and in the Standard Model....
Definition: NPEffectiveGIMR.cpp:3341
NPEffectiveGIMR::CLu_2211
double CLu_2211
Definition: NPEffectiveGIMR.h:1506
NPEffectiveGIMR::CHL1_13r
double CHL1_13r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1333
NPEffectiveGIMR::CLQ1_1331
double CLQ1_1331
Definition: NPEffectiveGIMR.h:1483
NPEffectiveGIMR::deltaG1_hZZ
virtual double deltaG1_hZZ() const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1843
NPEffectiveGIMR::deltaG2_hZZ
virtual double deltaG2_hZZ() const
The new physics contribution to the coupling of the effective interaction .
Definition: NPEffectiveGIMR.cpp:1848
NPEffectiveGIMR::CLL_up
double CLL_up() const
Definition: NPEffectiveGIMR.cpp:3367
NPEffectiveGIMR::deltaGR_f
double deltaGR_f(const Particle p) const
New physics contribution to the neutral-current right-handed coupling .
Definition: NPEffectiveGIMR.cpp:1770
NPEffectiveGIMR::CLR_up
double CLR_up() const
Definition: NPEffectiveGIMR.cpp:3402
NPEffectiveGIMR::CeH_23r
double CeH_23r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1410
NPEffectiveGIMR::CLd_1111
double CLd_1111
Definition: NPEffectiveGIMR.h:1508
NPEffectiveGIMR::LEP2AFBmu
virtual double LEP2AFBmu(const double s) const
Definition: NPEffectiveGIMR.cpp:3582
QCD::NEUTRINO_1
Definition: QCD.h:311
NPEffectiveGIMR::CHu_33
double CHu_33
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1381
QCD::quarks
Particle quarks[6]
The vector of all SM quarks.
Definition: QCD.h:934
NPEffectiveGIMR::CLd_1132
double CLd_1132
Definition: NPEffectiveGIMR.h:1512
NPEffectiveGIMR::BrHbbRatio
virtual double BrHbbRatio() const
The ratio of the Br in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:2994
NPEffectiveGIMR::eZH78_HZZ3
double eZH78_HZZ3
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at the LHC (7...
Definition: NPEffectiveGIMR.h:1587
NPEffectiveGIMR::ettH78_Htt
double ettH78_Htt
Theoretical uncertainty in the (linear) new physics contribution from to ttH production at the LHC (...
Definition: NPEffectiveGIMR.h:1601
NPEffectiveGIMR::CdH_13r
double CdH_13r
The dimension-6 operator coefficient (real part).
Definition: NPEffectiveGIMR.h:1432
NPEffectiveGIMR::CuW_13i
double CuW_13i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1462
QCD::MU
Definition: QCD.h:314
NPEffectiveGIMR::CHF3_diag
double CHF3_diag(const Particle F) const
The diagonal entry of the dimension-6 operator coefficient corresponding to particle F.
Definition: NPEffectiveGIMR.cpp:1523
NPEffectiveGIMR::CHQ1_11
double CHQ1_11
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1358
NPEffectiveGIMR::CHe_12i
double CHe_12i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1355
QCD::lepton
lepton
An enum type for leptons.
Definition: QCD.h:310
LEP2GIMR::sigma_l_LEP2_GIMR
double sigma_l_LEP2_GIMR(const QCD::lepton l, const double s, const double GIMRParam_i[]) const
Definition: LEP2GIMR.cpp:18
NPEffectiveGIMR::CHe_23i
double CHe_23i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1357
LEP2GIMR::sigmaFminusB_l_LEP2_GIMR
double sigmaFminusB_l_LEP2_GIMR(const QCD::lepton l, const double s, const double GIMRParam_i[]) const
Definition: LEP2GIMR.cpp:64
NPEffectiveGIMR::CH
double CH
The dimension-6 operator coefficient .
Definition: NPEffectiveGIMR.h:1330
NPEffectiveGIMR::delta_h
double delta_h
Combinations of dimension 6 coefficients modifying the canonical field definition.
Definition: NPEffectiveGIMR.h:1615
NPEffectiveGIMR::eZH2_HZA1
double eZH2_HZA1
Theoretical uncertainty in the (linear) new physics contribution from to ZH production at Tevatron (...
Definition: NPEffectiveGIMR.h:1575
NPEffectiveGIMR::CuB_11i
double CuB_11i
The dimension-6 operator coefficient (imaginary part).
Definition: NPEffectiveGIMR.h:1472
NPEffectiveGIMR::eVBF78_HZZ1
double eVBF78_HZZ1
Theoretical uncertainty in the (linear) new physics contribution from to VBF production at the LHC (...
Definition: NPEffectiveGIMR.h:1540
NPEffectiveGIMR::GammaHZgaRatio
double GammaHZgaRatio() const
The ratio of the in the current model and in the Standard Model.
Definition: NPEffectiveGIMR.cpp:3157