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