a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
NPd6SILH.cpp
Go to the documentation of this file.
1/* Copyright (C) 2025 HEPfit Collaboration
2*
3* For the licensing terms see doc/COPYING.
4*
5* Created on Tue 29 Jul 2025 16:05:01
6*
7*/
8
9#include "NPd6SILH.h"
10
11
12std::string NPd6SILH::NPd6SILHVars[NNPd6SILHVars] = {
13 "cH_LNP", "cT_LNP", "c6_LNP", "cB_LNP", "cW_LNP", "c2B_LNP", "c2W_LNP", "c2G_LNP", "c3W_LNP", "c3G_LNP", "cHW_LNP", "cHB_LNP", "cgam_LNP", "cg_LNP", "cHq1_LNP", "cHq3_LNP", "cHt_LNP", "ctD_LNP", "cqD1_LNP", "cqD3_LNP", "cqq1_LNP", "cqq3_LNP", "cqt1_LNP", "cqt8_LNP", "ctt_LNP", "ctG_LNP", "ctB_LNP", "ctW_LNP", "cu_LNP", "cd_LNP", "ce_LNP", "Lambda_NP"
14};
15
16
17NPd6SILH::NPd6SILH()
18: NPSMEFTd6General(), YuUV(3,3,0.), YuUVhc(3,3,0.), YdUV(3,3,0.), YdUVhc(3,3,0.), YeUV(3,3,0.), YeUVhc(3,3,0.) {
19 setModelName("NPd6SILH");
20
21 FlagRGEci = true;
22
23 ModelParamMap.insert(std::make_pair("cH_LNP", std::cref(cH_LNP)));
24 ModelParamMap.insert(std::make_pair("cT_LNP", std::cref(cT_LNP)));
25 ModelParamMap.insert(std::make_pair("c6_LNP", std::cref(c6_LNP)));
26 ModelParamMap.insert(std::make_pair("cB_LNP", std::cref(cB_LNP)));
27 ModelParamMap.insert(std::make_pair("cW_LNP", std::cref(cW_LNP)));
28 ModelParamMap.insert(std::make_pair("c2B_LNP", std::cref(c2B_LNP)));
29 ModelParamMap.insert(std::make_pair("c2W_LNP", std::cref(c2W_LNP)));
30 ModelParamMap.insert(std::make_pair("c2G_LNP", std::cref(c2G_LNP)));
31 ModelParamMap.insert(std::make_pair("c3W_LNP", std::cref(c3W_LNP)));
32 ModelParamMap.insert(std::make_pair("c3G_LNP", std::cref(c3G_LNP)));
33 ModelParamMap.insert(std::make_pair("cHW_LNP", std::cref(cHW_LNP)));
34 ModelParamMap.insert(std::make_pair("cHB_LNP", std::cref(cHB_LNP)));
35 ModelParamMap.insert(std::make_pair("cgam_LNP", std::cref(cgam_LNP)));
36 ModelParamMap.insert(std::make_pair("cg_LNP", std::cref(cg_LNP)));
37 ModelParamMap.insert(std::make_pair("cHq1_LNP", std::cref(cHq1_LNP)));
38 ModelParamMap.insert(std::make_pair("cHq3_LNP", std::cref(cHq3_LNP)));
39 ModelParamMap.insert(std::make_pair("cHt_LNP", std::cref(cHt_LNP)));
40 ModelParamMap.insert(std::make_pair("ctD_LNP", std::cref(ctD_LNP)));
41 ModelParamMap.insert(std::make_pair("cqD1_LNP", std::cref(cqD1_LNP)));
42 ModelParamMap.insert(std::make_pair("cqD3_LNP", std::cref(cqD3_LNP)));
43 ModelParamMap.insert(std::make_pair("cqq1_LNP", std::cref(cqq1_LNP)));
44 ModelParamMap.insert(std::make_pair("cqq3_LNP", std::cref(cqq3_LNP)));
45 ModelParamMap.insert(std::make_pair("cqt1_LNP", std::cref(cqt1_LNP)));
46 ModelParamMap.insert(std::make_pair("cqt8_LNP", std::cref(cqt8_LNP)));
47 ModelParamMap.insert(std::make_pair("ctt_LNP", std::cref(ctt_LNP)));
48 ModelParamMap.insert(std::make_pair("ctG_LNP", std::cref(ctG_LNP)));
49 ModelParamMap.insert(std::make_pair("ctB_LNP", std::cref(ctB_LNP)));
50 ModelParamMap.insert(std::make_pair("ctW_LNP", std::cref(ctW_LNP)));
51 ModelParamMap.insert(std::make_pair("cu_LNP", std::cref(cu_LNP)));
52 ModelParamMap.insert(std::make_pair("cd_LNP", std::cref(cd_LNP)));
53 ModelParamMap.insert(std::make_pair("ce_LNP", std::cref(ce_LNP)));
54
55}
56
57
58void NPd6SILH::setParameter(const std::string name, const double& value)
59{
60
61 if(name.compare("cH_LNP") == 0 ) {
62
63 cH_LNP = value;
64
65 } else if(name.compare("cT_LNP")==0) {
66
67 cT_LNP = value;
68
69 } else if(name.compare("c6_LNP")==0) {
70
71 c6_LNP = value;
72
73 } else if(name.compare("cB_LNP")==0) {
74
75 cB_LNP = value;
76
77 } else if(name.compare("cW_LNP")==0) {
78
79 cW_LNP = value;
80
81 } else if(name.compare("c2B_LNP")==0) {
82
83 c2B_LNP = value;
84
85 } else if(name.compare("c2W_LNP")==0) {
86
87 c2W_LNP = value;
88
89 } else if(name.compare("c2G_LNP")==0) {
90
91 c2G_LNP = value;
92
93 } else if(name.compare("c3W_LNP")==0) {
94
95 c3W_LNP = value;
96
97 } else if(name.compare("c3G_LNP")==0) {
98
99 c3G_LNP = value;
100
101 } else if(name.compare("cHW_LNP")==0) {
102
103 cHW_LNP = value;
104
105 } else if(name.compare("cHB_LNP")==0) {
106
107 cHB_LNP = value;
108
109 } else if(name.compare("cgam_LNP")==0) {
110
111 cgam_LNP = value;
112
113 } else if(name.compare("cg_LNP")==0) {
114
115 cg_LNP = value;
116
117 } else if(name.compare("cHq1_LNP")==0) {
118
119 cHq1_LNP = value;
120
121 } else if(name.compare("cHq3_LNP")==0) {
122
123 cHq3_LNP = value;
124
125 } else if(name.compare("cHt_LNP")==0) {
126
127 cHt_LNP = value;
128
129 } else if(name.compare("ctD_LNP")==0) {
130
131 ctD_LNP = value;
132
133 } else if(name.compare("cqD1_LNP")==0) {
134
135 cqD1_LNP = value;
136
137 } else if(name.compare("cqD3_LNP")==0) {
138
139 cqD3_LNP = value;
140
141 } else if(name.compare("cqq1_LNP")==0) {
142
143 cqq1_LNP = value;
144
145 } else if(name.compare("cqq3_LNP")==0) {
146
147 cqq3_LNP = value;
148
149 } else if(name.compare("cqt1_LNP")==0) {
150
151 cqt1_LNP = value;
152
153 } else if(name.compare("cqt8_LNP")==0) {
154
155 cqt8_LNP = value;
156
157 } else if(name.compare("ctt_LNP")==0) {
158
159 ctt_LNP = value;
160
161 } else if(name.compare("ctG_LNP")==0) {
162
163 ctG_LNP = value;
164
165 } else if(name.compare("ctB_LNP")==0) {
166
167 ctB_LNP = value;
168
169 } else if(name.compare("ctW_LNP")==0) {
170
171 ctW_LNP = value;
172
173 } else if(name.compare("cu_LNP")==0) {
174
175 cu_LNP = value;
176
177 } else if(name.compare("cd_LNP")==0) {
178
179 cd_LNP = value;
180
181 } else if(name.compare("ce_LNP")==0) {
182
183 ce_LNP = value;
184
185 } else if(name.compare("Lambda_NP")==0) {
186
187 Lambda_NP = value;
188
189 } else {
190 NPSMEFTd6General::setParameter(name, value);
191 }
192}
193
194
195bool NPd6SILH::setFlag(const std::string name, const bool value) {
196 bool res = false;
197 if (name.compare("RGEci") == 0) {
198 FlagRGEci = value;
199 //We need to fix FlagRGEci also in the NPSMEFTd6General
200 res = NPSMEFTd6General::setFlag(name, value);
201
202 } else
203 res = NPSMEFTd6General::setFlag(name, value);
204
205 return (res);
206}
207
208
209void NPd6SILH::setNPSMEFTd6GeneralParameters()
210{
211 // SM parameters at the UV scale
212 g1UV = getSMEFTCoeffEW("g1");
213 g2UV = getSMEFTCoeffEW("g2");
214 g3UV = getSMEFTCoeffEW("g3");
215 lambdaHUV = getSMEFTCoeffEW("lambda");
216
217 g1UV2 = g1UV * g1UV;
218 g2UV2 = g2UV * g2UV;
219 g3UV2 = g3UV * g3UV;
220
221 for (int i = 0; i < 3; i++) {
222 for (int j = 0; j < 3; j++) {
223 YuUV.assignre(i, j, getSMEFTCoeffEW("YuR", i, j) );
224 YuUV.assignim(i, j, getSMEFTCoeffEW("YuI", i, j) );
225 YdUV.assignre(i, j, getSMEFTCoeffEW("YdR", i, j) );
226 YdUV.assignim(i, j, getSMEFTCoeffEW("YdI", i, j) );
227 YeUV.assignre(i, j, getSMEFTCoeffEW("YeR", i, j) );
228 YeUV.assignim(i, j, getSMEFTCoeffEW("YeI", i, j) );
229 }
230 }
231
232 YuUVhc = YuUV.hconjugate();
233 YdUVhc = YdUV.hconjugate();
234 YeUVhc = YeUV.hconjugate();
235
236
237 // Matching
238
239 CG_LNP = c3G_LNP;
240 CW_LNP = c3W_LNP;
241 CH_LNP = -c6_LNP - 2*cHW_LNP*g2UV*lambdaHUV - 4*cW_LNP*g2UV*lambdaHUV - 2*c2W_LNP*g2UV2*lambdaHUV;
242 CHbox_LNP = -cH_LNP - cT_LNP - (cB_LNP*g1UV)/2. - (cHB_LNP*g1UV)/4. - (c2B_LNP*g1UV2)/4. - (3*cHW_LNP*g2UV)/4. - (3*cW_LNP*g2UV)/2. - (3*c2W_LNP*g2UV2)/4.;
243 CHD_LNP = -4*cT_LNP - 2*cB_LNP*g1UV - cHB_LNP*g1UV - c2B_LNP*g1UV2;
244 CHG_LNP = cg_LNP;
245 CHB_LNP = cgam_LNP + (cHB_LNP*g1UV)/4.;
246 CHW_LNP = (cHW_LNP*g2UV)/4.;
247 CHWB_LNP = (cHW_LNP*g1UV)/4. + (cHB_LNP*g2UV)/4.;
248 CuH_11r_LNP = ( cu_LNP*YuUVhc(0,0) - (cHW_LNP*g2UV*YuUVhc(0,0))/2. - cW_LNP*g2UV*YuUVhc(0,0) - (c2W_LNP*g2UV2*YuUVhc(0,0))/2. ).real();
249 CuH_11i_LNP = ( cu_LNP*YuUVhc(0,0) - (cHW_LNP*g2UV*YuUVhc(0,0))/2. - cW_LNP*g2UV*YuUVhc(0,0) - (c2W_LNP*g2UV2*YuUVhc(0,0))/2. ).imag();
250 CuH_12r_LNP = ( cu_LNP*YuUVhc(0,1) - (cHW_LNP*g2UV*YuUVhc(0,1))/2. - cW_LNP*g2UV*YuUVhc(0,1) - (c2W_LNP*g2UV2*YuUVhc(0,1))/2. ).real();
251 CuH_12i_LNP = ( cu_LNP*YuUVhc(0,1) - (cHW_LNP*g2UV*YuUVhc(0,1))/2. - cW_LNP*g2UV*YuUVhc(0,1) - (c2W_LNP*g2UV2*YuUVhc(0,1))/2. ).imag();
252 CuH_13r_LNP = ( cu_LNP*YuUVhc(0,2) - (cHW_LNP*g2UV*YuUVhc(0,2))/2. - cW_LNP*g2UV*YuUVhc(0,2) - (c2W_LNP*g2UV2*YuUVhc(0,2))/2. ).real();
253 CuH_13i_LNP = ( cu_LNP*YuUVhc(0,2) - (cHW_LNP*g2UV*YuUVhc(0,2))/2. - cW_LNP*g2UV*YuUVhc(0,2) - (c2W_LNP*g2UV2*YuUVhc(0,2))/2. ).imag();
254 CuH_21r_LNP = ( cu_LNP*YuUVhc(1,0) - (cHW_LNP*g2UV*YuUVhc(1,0))/2. - cW_LNP*g2UV*YuUVhc(1,0) - (c2W_LNP*g2UV2*YuUVhc(1,0))/2. ).real();
255 CuH_21i_LNP = ( cu_LNP*YuUVhc(1,0) - (cHW_LNP*g2UV*YuUVhc(1,0))/2. - cW_LNP*g2UV*YuUVhc(1,0) - (c2W_LNP*g2UV2*YuUVhc(1,0))/2. ).imag();
256 CuH_22r_LNP = ( cu_LNP*YuUVhc(1,1) - (cHW_LNP*g2UV*YuUVhc(1,1))/2. - cW_LNP*g2UV*YuUVhc(1,1) - (c2W_LNP*g2UV2*YuUVhc(1,1))/2. ).real();
257 CuH_22i_LNP = ( cu_LNP*YuUVhc(1,1) - (cHW_LNP*g2UV*YuUVhc(1,1))/2. - cW_LNP*g2UV*YuUVhc(1,1) - (c2W_LNP*g2UV2*YuUVhc(1,1))/2. ).imag();
258 CuH_23r_LNP = ( cu_LNP*YuUVhc(1,2) - (cHW_LNP*g2UV*YuUVhc(1,2))/2. - cW_LNP*g2UV*YuUVhc(1,2) - (c2W_LNP*g2UV2*YuUVhc(1,2))/2. ).real();
259 CuH_23i_LNP = ( cu_LNP*YuUVhc(1,2) - (cHW_LNP*g2UV*YuUVhc(1,2))/2. - cW_LNP*g2UV*YuUVhc(1,2) - (c2W_LNP*g2UV2*YuUVhc(1,2))/2. ).imag();
260 CuH_31r_LNP = ( cu_LNP*YuUVhc(2,0) - (cHW_LNP*g2UV*YuUVhc(2,0))/2. - cW_LNP*g2UV*YuUVhc(2,0) - (c2W_LNP*g2UV2*YuUVhc(2,0))/2. ).real();
261 CuH_31i_LNP = ( cu_LNP*YuUVhc(2,0) - (cHW_LNP*g2UV*YuUVhc(2,0))/2. - cW_LNP*g2UV*YuUVhc(2,0) - (c2W_LNP*g2UV2*YuUVhc(2,0))/2. ).imag();
262 CuH_32r_LNP = ( cu_LNP*YuUVhc(2,1) - (cHW_LNP*g2UV*YuUVhc(2,1))/2. - cW_LNP*g2UV*YuUVhc(2,1) - (c2W_LNP*g2UV2*YuUVhc(2,1))/2. ).real();
263 CuH_32i_LNP = ( cu_LNP*YuUVhc(2,1) - (cHW_LNP*g2UV*YuUVhc(2,1))/2. - cW_LNP*g2UV*YuUVhc(2,1) - (c2W_LNP*g2UV2*YuUVhc(2,1))/2. ).imag();
264 CuH_33r_LNP = ( cu_LNP*YuUVhc(2,2) - (cHW_LNP*g2UV*YuUVhc(2,2))/2. - cW_LNP*g2UV*YuUVhc(2,2) - (c2W_LNP*g2UV2*YuUVhc(2,2))/2. ).real();
265 CuH_33i_LNP = ( cu_LNP*YuUVhc(2,2) - (cHW_LNP*g2UV*YuUVhc(2,2))/2. - cW_LNP*g2UV*YuUVhc(2,2) - (c2W_LNP*g2UV2*YuUVhc(2,2))/2. ).imag();
266 CdH_11r_LNP = ( cd_LNP*YdUVhc(0,0) - (cHW_LNP*g2UV*YdUVhc(0,0))/2. - cW_LNP*g2UV*YdUVhc(0,0) - (c2W_LNP*g2UV2*YdUVhc(0,0))/2. ).real();
267 CdH_11i_LNP = ( cd_LNP*YdUVhc(0,0) - (cHW_LNP*g2UV*YdUVhc(0,0))/2. - cW_LNP*g2UV*YdUVhc(0,0) - (c2W_LNP*g2UV2*YdUVhc(0,0))/2. ).imag();
268 CdH_12r_LNP = ( cd_LNP*YdUVhc(0,1) - (cHW_LNP*g2UV*YdUVhc(0,1))/2. - cW_LNP*g2UV*YdUVhc(0,1) - (c2W_LNP*g2UV2*YdUVhc(0,1))/2. ).real();
269 CdH_12i_LNP = ( cd_LNP*YdUVhc(0,1) - (cHW_LNP*g2UV*YdUVhc(0,1))/2. - cW_LNP*g2UV*YdUVhc(0,1) - (c2W_LNP*g2UV2*YdUVhc(0,1))/2. ).imag();
270 CdH_13r_LNP = ( cd_LNP*YdUVhc(0,2) - (cHW_LNP*g2UV*YdUVhc(0,2))/2. - cW_LNP*g2UV*YdUVhc(0,2) - (c2W_LNP*g2UV2*YdUVhc(0,2))/2. ).real();
271 CdH_13i_LNP = ( cd_LNP*YdUVhc(0,2) - (cHW_LNP*g2UV*YdUVhc(0,2))/2. - cW_LNP*g2UV*YdUVhc(0,2) - (c2W_LNP*g2UV2*YdUVhc(0,2))/2. ).imag();
272 CdH_21r_LNP = ( cd_LNP*YdUVhc(1,0) - (cHW_LNP*g2UV*YdUVhc(1,0))/2. - cW_LNP*g2UV*YdUVhc(1,0) - (c2W_LNP*g2UV2*YdUVhc(1,0))/2. ).real();
273 CdH_21i_LNP = ( cd_LNP*YdUVhc(1,0) - (cHW_LNP*g2UV*YdUVhc(1,0))/2. - cW_LNP*g2UV*YdUVhc(1,0) - (c2W_LNP*g2UV2*YdUVhc(1,0))/2. ).imag();
274 CdH_22r_LNP = ( cd_LNP*YdUVhc(1,1) - (cHW_LNP*g2UV*YdUVhc(1,1))/2. - cW_LNP*g2UV*YdUVhc(1,1) - (c2W_LNP*g2UV2*YdUVhc(1,1))/2. ).real();
275 CdH_22i_LNP = ( cd_LNP*YdUVhc(1,1) - (cHW_LNP*g2UV*YdUVhc(1,1))/2. - cW_LNP*g2UV*YdUVhc(1,1) - (c2W_LNP*g2UV2*YdUVhc(1,1))/2. ).imag();
276 CdH_23r_LNP = ( cd_LNP*YdUVhc(1,2) - (cHW_LNP*g2UV*YdUVhc(1,2))/2. - cW_LNP*g2UV*YdUVhc(1,2) - (c2W_LNP*g2UV2*YdUVhc(1,2))/2. ).real();
277 CdH_23i_LNP = ( cd_LNP*YdUVhc(1,2) - (cHW_LNP*g2UV*YdUVhc(1,2))/2. - cW_LNP*g2UV*YdUVhc(1,2) - (c2W_LNP*g2UV2*YdUVhc(1,2))/2. ).imag();
278 CdH_31r_LNP = ( cd_LNP*YdUVhc(2,0) - (cHW_LNP*g2UV*YdUVhc(2,0))/2. - cW_LNP*g2UV*YdUVhc(2,0) - (c2W_LNP*g2UV2*YdUVhc(2,0))/2. ).real();
279 CdH_31i_LNP = ( cd_LNP*YdUVhc(2,0) - (cHW_LNP*g2UV*YdUVhc(2,0))/2. - cW_LNP*g2UV*YdUVhc(2,0) - (c2W_LNP*g2UV2*YdUVhc(2,0))/2. ).imag();
280 CdH_32r_LNP = ( cd_LNP*YdUVhc(2,1) - (cHW_LNP*g2UV*YdUVhc(2,1))/2. - cW_LNP*g2UV*YdUVhc(2,1) - (c2W_LNP*g2UV2*YdUVhc(2,1))/2. ).real();
281 CdH_32i_LNP = ( cd_LNP*YdUVhc(2,1) - (cHW_LNP*g2UV*YdUVhc(2,1))/2. - cW_LNP*g2UV*YdUVhc(2,1) - (c2W_LNP*g2UV2*YdUVhc(2,1))/2. ).imag();
282 CdH_33r_LNP = ( cd_LNP*YdUVhc(2,2) - (cHW_LNP*g2UV*YdUVhc(2,2))/2. - cW_LNP*g2UV*YdUVhc(2,2) - (c2W_LNP*g2UV2*YdUVhc(2,2))/2. ).real();
283 CdH_33i_LNP = ( cd_LNP*YdUVhc(2,2) - (cHW_LNP*g2UV*YdUVhc(2,2))/2. - cW_LNP*g2UV*YdUVhc(2,2) - (c2W_LNP*g2UV2*YdUVhc(2,2))/2. ).imag();
284 CeH_11r_LNP = ( ce_LNP*YeUVhc(0,0) - (cHW_LNP*g2UV*YeUVhc(0,0))/2. - cW_LNP*g2UV*YeUVhc(0,0) - (c2W_LNP*g2UV2*YeUVhc(0,0))/2. ).real();
285 CeH_11i_LNP = ( ce_LNP*YeUVhc(0,0) - (cHW_LNP*g2UV*YeUVhc(0,0))/2. - cW_LNP*g2UV*YeUVhc(0,0) - (c2W_LNP*g2UV2*YeUVhc(0,0))/2. ).imag();
286 CeH_12r_LNP = ( ce_LNP*YeUVhc(0,1) - (cHW_LNP*g2UV*YeUVhc(0,1))/2. - cW_LNP*g2UV*YeUVhc(0,1) - (c2W_LNP*g2UV2*YeUVhc(0,1))/2. ).real();
287 CeH_12i_LNP = ( ce_LNP*YeUVhc(0,1) - (cHW_LNP*g2UV*YeUVhc(0,1))/2. - cW_LNP*g2UV*YeUVhc(0,1) - (c2W_LNP*g2UV2*YeUVhc(0,1))/2. ).imag();
288 CeH_13r_LNP = ( ce_LNP*YeUVhc(0,2) - (cHW_LNP*g2UV*YeUVhc(0,2))/2. - cW_LNP*g2UV*YeUVhc(0,2) - (c2W_LNP*g2UV2*YeUVhc(0,2))/2. ).real();
289 CeH_13i_LNP = ( ce_LNP*YeUVhc(0,2) - (cHW_LNP*g2UV*YeUVhc(0,2))/2. - cW_LNP*g2UV*YeUVhc(0,2) - (c2W_LNP*g2UV2*YeUVhc(0,2))/2. ).imag();
290 CeH_21r_LNP = ( ce_LNP*YeUVhc(1,0) - (cHW_LNP*g2UV*YeUVhc(1,0))/2. - cW_LNP*g2UV*YeUVhc(1,0) - (c2W_LNP*g2UV2*YeUVhc(1,0))/2. ).real();
291 CeH_21i_LNP = ( ce_LNP*YeUVhc(1,0) - (cHW_LNP*g2UV*YeUVhc(1,0))/2. - cW_LNP*g2UV*YeUVhc(1,0) - (c2W_LNP*g2UV2*YeUVhc(1,0))/2. ).imag();
292 CeH_22r_LNP = ( ce_LNP*YeUVhc(1,1) - (cHW_LNP*g2UV*YeUVhc(1,1))/2. - cW_LNP*g2UV*YeUVhc(1,1) - (c2W_LNP*g2UV2*YeUVhc(1,1))/2. ).real();
293 CeH_22i_LNP = ( ce_LNP*YeUVhc(1,1) - (cHW_LNP*g2UV*YeUVhc(1,1))/2. - cW_LNP*g2UV*YeUVhc(1,1) - (c2W_LNP*g2UV2*YeUVhc(1,1))/2. ).imag();
294 CeH_23r_LNP = ( ce_LNP*YeUVhc(1,2) - (cHW_LNP*g2UV*YeUVhc(1,2))/2. - cW_LNP*g2UV*YeUVhc(1,2) - (c2W_LNP*g2UV2*YeUVhc(1,2))/2. ).real();
295 CeH_23i_LNP = ( ce_LNP*YeUVhc(1,2) - (cHW_LNP*g2UV*YeUVhc(1,2))/2. - cW_LNP*g2UV*YeUVhc(1,2) - (c2W_LNP*g2UV2*YeUVhc(1,2))/2. ).imag();
296 CeH_31r_LNP = ( ce_LNP*YeUVhc(2,0) - (cHW_LNP*g2UV*YeUVhc(2,0))/2. - cW_LNP*g2UV*YeUVhc(2,0) - (c2W_LNP*g2UV2*YeUVhc(2,0))/2. ).real();
297 CeH_31i_LNP = ( ce_LNP*YeUVhc(2,0) - (cHW_LNP*g2UV*YeUVhc(2,0))/2. - cW_LNP*g2UV*YeUVhc(2,0) - (c2W_LNP*g2UV2*YeUVhc(2,0))/2. ).imag();
298 CeH_32r_LNP = ( ce_LNP*YeUVhc(2,1) - (cHW_LNP*g2UV*YeUVhc(2,1))/2. - cW_LNP*g2UV*YeUVhc(2,1) - (c2W_LNP*g2UV2*YeUVhc(2,1))/2. ).real();
299 CeH_32i_LNP = ( ce_LNP*YeUVhc(2,1) - (cHW_LNP*g2UV*YeUVhc(2,1))/2. - cW_LNP*g2UV*YeUVhc(2,1) - (c2W_LNP*g2UV2*YeUVhc(2,1))/2. ).imag();
300 CeH_33r_LNP = ( ce_LNP*YeUVhc(2,2) - (cHW_LNP*g2UV*YeUVhc(2,2))/2. - cW_LNP*g2UV*YeUVhc(2,2) - (c2W_LNP*g2UV2*YeUVhc(2,2))/2. ).real();
301 CeH_33i_LNP = ( ce_LNP*YeUVhc(2,2) - (cHW_LNP*g2UV*YeUVhc(2,2))/2. - cW_LNP*g2UV*YeUVhc(2,2) - (c2W_LNP*g2UV2*YeUVhc(2,2))/2. ).imag();
302 CuG_33r_LNP = ctG_LNP;
303 CuG_33i_LNP = 0.0;
304 CuW_33r_LNP = ctW_LNP;
305 CuW_33i_LNP = 0.0;
306 CuB_33r_LNP = ctB_LNP;
307 CuB_33i_LNP = 0.0;
308 CHl1_11r_LNP = (cB_LNP*g1UV)/2. + (cHB_LNP*g1UV)/4. + (c2B_LNP*g1UV2)/2.;
309 CHl1_22r_LNP = (cB_LNP*g1UV)/2. + (cHB_LNP*g1UV)/4. + (c2B_LNP*g1UV2)/2.;
310 CHl1_33r_LNP = (cB_LNP*g1UV)/2. + (cHB_LNP*g1UV)/4. + (c2B_LNP*g1UV2)/2.;
311 CHl3_11r_LNP = -0.25*(cHW_LNP*g2UV) - (cW_LNP*g2UV)/2. - (c2W_LNP*g2UV2)/2.;
312 CHl3_22r_LNP = -0.25*(cHW_LNP*g2UV) - (cW_LNP*g2UV)/2. - (c2W_LNP*g2UV2)/2.;
313 CHl3_33r_LNP = -0.25*(cHW_LNP*g2UV) - (cW_LNP*g2UV)/2. - (c2W_LNP*g2UV2)/2.;
314 CHe_11r_LNP = cB_LNP*g1UV + (cHB_LNP*g1UV)/2. + c2B_LNP*g1UV2;
315 CHe_22r_LNP = cB_LNP*g1UV + (cHB_LNP*g1UV)/2. + c2B_LNP*g1UV2;
316 CHe_33r_LNP = cB_LNP*g1UV + (cHB_LNP*g1UV)/2. + c2B_LNP*g1UV2;
317 CHq1_11r_LNP = -0.16666666666666666*(cB_LNP*g1UV) - (cHB_LNP*g1UV)/12. - (c2B_LNP*g1UV2)/6.;
318 CHq1_22r_LNP = -0.16666666666666666*(cB_LNP*g1UV) - (cHB_LNP*g1UV)/12. - (c2B_LNP*g1UV2)/6.;
319 CHq1_33r_LNP = cHq1_LNP - (cB_LNP*g1UV)/6. - (cHB_LNP*g1UV)/12. - (cqD1_LNP*g1UV)/2. - (c2B_LNP*g1UV2)/6.;
320 CHq3_11r_LNP = -0.25*(cHW_LNP*g2UV) - (cW_LNP*g2UV)/2. - (c2W_LNP*g2UV2)/2.;
321 CHq3_22r_LNP = -0.25*(cHW_LNP*g2UV) - (cW_LNP*g2UV)/2. - (c2W_LNP*g2UV2)/2.;
322 CHq3_33r_LNP = cHq3_LNP - (cHW_LNP*g2UV)/4. - (cqD3_LNP*g2UV)/2. - (cW_LNP*g2UV)/2. - (c2W_LNP*g2UV2)/2.;
323 CHu_11r_LNP = (-2*cB_LNP*g1UV)/3. - (cHB_LNP*g1UV)/3. - (2*c2B_LNP*g1UV2)/3.;
324 CHu_22r_LNP = (-2*cB_LNP*g1UV)/3. - (cHB_LNP*g1UV)/3. - (2*c2B_LNP*g1UV2)/3.;
325 CHu_33r_LNP = cHt_LNP - (2*cB_LNP*g1UV)/3. - (cHB_LNP*g1UV)/3. - (ctD_LNP*g1UV)/2. - (2*c2B_LNP*g1UV2)/3.;
326 CHd_11r_LNP = (cB_LNP*g1UV)/3. + (cHB_LNP*g1UV)/6. + (c2B_LNP*g1UV2)/3.;
327 CHd_22r_LNP = (cB_LNP*g1UV)/3. + (cHB_LNP*g1UV)/6. + (c2B_LNP*g1UV2)/3.;
328 CHd_33r_LNP = (cB_LNP*g1UV)/3. + (cHB_LNP*g1UV)/6. + (c2B_LNP*g1UV2)/3.;
329 Cll_1111r_LNP = -0.25*(c2B_LNP*g1UV2) - (c2W_LNP*g2UV2)/4.;
330 Cll_1122r_LNP = -0.25*(c2B_LNP*g1UV2) + (c2W_LNP*g2UV2)/4.;
331 Cll_1133r_LNP = -0.25*(c2B_LNP*g1UV2) + (c2W_LNP*g2UV2)/4.;
332 Cll_1221r_LNP = -0.5*(c2W_LNP*g2UV2);
333 Cll_1331r_LNP = -0.5*(c2W_LNP*g2UV2);
334 Cll_2222r_LNP = -0.25*(c2B_LNP*g1UV2) - (c2W_LNP*g2UV2)/4.;
335 Cll_2233r_LNP = -0.25*(c2B_LNP*g1UV2) + (c2W_LNP*g2UV2)/4.;
336 Cll_2332r_LNP = -0.5*(c2W_LNP*g2UV2);
337 Cll_3333r_LNP = -0.25*(c2B_LNP*g1UV2) - (c2W_LNP*g2UV2)/4.;
338 Cqq1_1111r_LNP = -0.027777777777777776*(c2B_LNP*g1UV2) - (c2G_LNP*g3UV2)/12.;
339 Cqq1_1122r_LNP = -0.027777777777777776*(c2B_LNP*g1UV2) + (c2G_LNP*g3UV2)/6.;
340 Cqq1_1133r_LNP = -0.16666666666666666*(cqD1_LNP*g1UV) - (c2B_LNP*g1UV2)/36. + (c2G_LNP*g3UV2)/6.;
341 Cqq1_1221r_LNP = -0.25*(c2G_LNP*g3UV2);
342 Cqq1_1331r_LNP = -0.25*(c2G_LNP*g3UV2);
343 Cqq1_2222r_LNP = -0.027777777777777776*(c2B_LNP*g1UV2) - (c2G_LNP*g3UV2)/12.;
344 Cqq1_2233r_LNP = -0.16666666666666666*(cqD1_LNP*g1UV) - (c2B_LNP*g1UV2)/36. + (c2G_LNP*g3UV2)/6.;
345 Cqq1_2332r_LNP = -0.25*(c2G_LNP*g3UV2);
346 Cqq1_3333r_LNP = cqq1_LNP - (cqD1_LNP*g1UV)/6. - (c2B_LNP*g1UV2)/36. - (c2G_LNP*g3UV2)/12.;
347 Cqq3_1111r_LNP = -(c2W_LNP*g2UV2) - (c2G_LNP*g3UV2)/4.;
348 Cqq3_1122r_LNP = -(c2W_LNP*g2UV2);
349 Cqq3_1133r_LNP = -0.5*(cqD3_LNP*g2UV) - c2W_LNP*g2UV2;
350 Cqq3_1221r_LNP = -0.25*(c2G_LNP*g3UV2);
351 Cqq3_1331r_LNP = -0.25*(c2G_LNP*g3UV2);
352 Cqq3_2222r_LNP = -(c2W_LNP*g2UV2) - (c2G_LNP*g3UV2)/4.;
353 Cqq3_2233r_LNP = -0.5*(cqD3_LNP*g2UV) - c2W_LNP*g2UV2;
354 Cqq3_2332r_LNP = -0.25*(c2G_LNP*g3UV2);
355 Cqq3_3333r_LNP = cqq3_LNP - (cqD3_LNP*g2UV)/2. - c2W_LNP*g2UV2 - (c2G_LNP*g3UV2)/4.;
356 Clq1_1111r_LNP = (c2B_LNP*g1UV2)/6.;
357 Clq1_1122r_LNP = (c2B_LNP*g1UV2)/6.;
358 Clq1_1133r_LNP = (c2B_LNP*g1UV2)/6.;
359 Clq1_2211r_LNP = (c2B_LNP*g1UV2)/6.;
360 Clq1_2222r_LNP = (c2B_LNP*g1UV2)/6.;
361 Clq1_2233r_LNP = (c2B_LNP*g1UV2)/6.;
362 Clq1_3311r_LNP = (c2B_LNP*g1UV2)/6.;
363 Clq1_3322r_LNP = (c2B_LNP*g1UV2)/6.;
364 Clq1_3333r_LNP = (c2B_LNP*g1UV2)/6.;
365 Clq3_1111r_LNP = -0.5*(c2W_LNP*g2UV2);
366 Clq3_1122r_LNP = -0.5*(c2W_LNP*g2UV2);
367 Clq3_1133r_LNP = -0.5*(c2W_LNP*g2UV2);
368 Clq3_2211r_LNP = -0.5*(c2W_LNP*g2UV2);
369 Clq3_2222r_LNP = -0.5*(c2W_LNP*g2UV2);
370 Clq3_2233r_LNP = -0.5*(c2W_LNP*g2UV2);
371 Clq3_3311r_LNP = -0.5*(c2W_LNP*g2UV2);
372 Clq3_3322r_LNP = -0.5*(c2W_LNP*g2UV2);
373 Clq3_3333r_LNP = -0.5*(c2W_LNP*g2UV2);
374 Cee_1111r_LNP = -(c2B_LNP*g1UV2);
375 Cee_1122r_LNP = -0.5*(c2B_LNP*g1UV2);
376 Cee_1133r_LNP = -0.5*(c2B_LNP*g1UV2);
377 Cee_2222r_LNP = -(c2B_LNP*g1UV2);
378 Cee_2233r_LNP = -0.5*(c2B_LNP*g1UV2);
379 Cee_3333r_LNP = -(c2B_LNP*g1UV2);
380 Cuu_1111r_LNP = (-4*c2B_LNP*g1UV2)/9. - (c2G_LNP*g3UV2)/3.;
381 Cuu_1122r_LNP = (-4*c2B_LNP*g1UV2)/9. + (c2G_LNP*g3UV2)/6.;
382 Cuu_1133r_LNP = (-2*ctD_LNP*g1UV)/3. - (4*c2B_LNP*g1UV2)/9. + (c2G_LNP*g3UV2)/6.;
383 Cuu_1221r_LNP = -0.5*(c2G_LNP*g3UV2);
384 Cuu_1331r_LNP = -0.5*(c2G_LNP*g3UV2);
385 Cuu_2222r_LNP = (-4*c2B_LNP*g1UV2)/9. - (c2G_LNP*g3UV2)/3.;
386 Cuu_2233r_LNP = (-2*ctD_LNP*g1UV)/3. - (4*c2B_LNP*g1UV2)/9. + (c2G_LNP*g3UV2)/6.;
387 Cuu_2332r_LNP = -0.5*(c2G_LNP*g3UV2);
388 Cuu_3333r_LNP = ctt_LNP - (2*ctD_LNP*g1UV)/3. - (4*c2B_LNP*g1UV2)/9. - (c2G_LNP*g3UV2)/3.;
389 Cdd_1111r_LNP = -0.1111111111111111*(c2B_LNP*g1UV2) - (c2G_LNP*g3UV2)/3.;
390 Cdd_1122r_LNP = -0.1111111111111111*(c2B_LNP*g1UV2) + (c2G_LNP*g3UV2)/6.;
391 Cdd_1133r_LNP = -0.1111111111111111*(c2B_LNP*g1UV2) + (c2G_LNP*g3UV2)/6.;
392 Cdd_1221r_LNP = -0.5*(c2G_LNP*g3UV2);
393 Cdd_1331r_LNP = -0.5*(c2G_LNP*g3UV2);
394 Cdd_2222r_LNP = -0.1111111111111111*(c2B_LNP*g1UV2) - (c2G_LNP*g3UV2)/3.;
395 Cdd_2233r_LNP = -0.1111111111111111*(c2B_LNP*g1UV2) + (c2G_LNP*g3UV2)/6.;
396 Cdd_2332r_LNP = -0.5*(c2G_LNP*g3UV2);
397 Cdd_3333r_LNP = -0.1111111111111111*(c2B_LNP*g1UV2) - (c2G_LNP*g3UV2)/3.;
398 Ceu_1111r_LNP = (4*c2B_LNP*g1UV2)/3.;
399 Ceu_1122r_LNP = (4*c2B_LNP*g1UV2)/3.;
400 Ceu_1133r_LNP = ctD_LNP*g1UV + (4*c2B_LNP*g1UV2)/3.;
401 Ceu_2211r_LNP = (4*c2B_LNP*g1UV2)/3.;
402 Ceu_2222r_LNP = (4*c2B_LNP*g1UV2)/3.;
403 Ceu_2233r_LNP = ctD_LNP*g1UV + (4*c2B_LNP*g1UV2)/3.;
404 Ceu_3311r_LNP = (4*c2B_LNP*g1UV2)/3.;
405 Ceu_3322r_LNP = (4*c2B_LNP*g1UV2)/3.;
406 Ceu_3333r_LNP = ctD_LNP*g1UV + (4*c2B_LNP*g1UV2)/3.;
407 Ced_1111r_LNP = (-2*c2B_LNP*g1UV2)/3.;
408 Ced_1122r_LNP = (-2*c2B_LNP*g1UV2)/3.;
409 Ced_1133r_LNP = (-2*c2B_LNP*g1UV2)/3.;
410 Ced_2211r_LNP = (-2*c2B_LNP*g1UV2)/3.;
411 Ced_2222r_LNP = (-2*c2B_LNP*g1UV2)/3.;
412 Ced_2233r_LNP = (-2*c2B_LNP*g1UV2)/3.;
413 Ced_3311r_LNP = (-2*c2B_LNP*g1UV2)/3.;
414 Ced_3322r_LNP = (-2*c2B_LNP*g1UV2)/3.;
415 Ced_3333r_LNP = (-2*c2B_LNP*g1UV2)/3.;
416 Cud1_1111r_LNP = (4*c2B_LNP*g1UV2)/9.;
417 Cud1_1122r_LNP = (4*c2B_LNP*g1UV2)/9.;
418 Cud1_1133r_LNP = (4*c2B_LNP*g1UV2)/9.;
419 Cud1_2211r_LNP = (4*c2B_LNP*g1UV2)/9.;
420 Cud1_2222r_LNP = (4*c2B_LNP*g1UV2)/9.;
421 Cud1_2233r_LNP = (4*c2B_LNP*g1UV2)/9.;
422 Cud1_3311r_LNP = (ctD_LNP*g1UV)/3. + (4*c2B_LNP*g1UV2)/9.;
423 Cud1_3322r_LNP = (ctD_LNP*g1UV)/3. + (4*c2B_LNP*g1UV2)/9.;
424 Cud1_3333r_LNP = (ctD_LNP*g1UV)/3. + (4*c2B_LNP*g1UV2)/9.;
425 Cud8_1111r_LNP = -2*c2G_LNP*g3UV2;
426 Cud8_1122r_LNP = -2*c2G_LNP*g3UV2;
427 Cud8_1133r_LNP = -2*c2G_LNP*g3UV2;
428 Cud8_2211r_LNP = -2*c2G_LNP*g3UV2;
429 Cud8_2222r_LNP = -2*c2G_LNP*g3UV2;
430 Cud8_2233r_LNP = -2*c2G_LNP*g3UV2;
431 Cud8_3311r_LNP = -2*c2G_LNP*g3UV2;
432 Cud8_3322r_LNP = -2*c2G_LNP*g3UV2;
433 Cud8_3333r_LNP = -2*c2G_LNP*g3UV2;
434 Cle_1111r_LNP = -(c2B_LNP*g1UV2);
435 Cle_1122r_LNP = -(c2B_LNP*g1UV2);
436 Cle_1133r_LNP = -(c2B_LNP*g1UV2);
437 Cle_2211r_LNP = -(c2B_LNP*g1UV2);
438 Cle_2222r_LNP = -(c2B_LNP*g1UV2);
439 Cle_2233r_LNP = -(c2B_LNP*g1UV2);
440 Cle_3311r_LNP = -(c2B_LNP*g1UV2);
441 Cle_3322r_LNP = -(c2B_LNP*g1UV2);
442 Cle_3333r_LNP = -(c2B_LNP*g1UV2);
443 Clu_1111r_LNP = (2*c2B_LNP*g1UV2)/3.;
444 Clu_1122r_LNP = (2*c2B_LNP*g1UV2)/3.;
445 Clu_1133r_LNP = (ctD_LNP*g1UV)/2. + (2*c2B_LNP*g1UV2)/3.;
446 Clu_2211r_LNP = (2*c2B_LNP*g1UV2)/3.;
447 Clu_2222r_LNP = (2*c2B_LNP*g1UV2)/3.;
448 Clu_2233r_LNP = (ctD_LNP*g1UV)/2. + (2*c2B_LNP*g1UV2)/3.;
449 Clu_3311r_LNP = (2*c2B_LNP*g1UV2)/3.;
450 Clu_3322r_LNP = (2*c2B_LNP*g1UV2)/3.;
451 Clu_3333r_LNP = (ctD_LNP*g1UV)/2. + (2*c2B_LNP*g1UV2)/3.;
452 Cld_1111r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
453 Cld_1122r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
454 Cld_1133r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
455 Cld_2211r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
456 Cld_2222r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
457 Cld_2233r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
458 Cld_3311r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
459 Cld_3322r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
460 Cld_3333r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
461 Cqe_1111r_LNP = (c2B_LNP*g1UV2)/3.;
462 Cqe_1122r_LNP = (c2B_LNP*g1UV2)/3.;
463 Cqe_1133r_LNP = (c2B_LNP*g1UV2)/3.;
464 Cqe_2211r_LNP = (c2B_LNP*g1UV2)/3.;
465 Cqe_2222r_LNP = (c2B_LNP*g1UV2)/3.;
466 Cqe_2233r_LNP = (c2B_LNP*g1UV2)/3.;
467 Cqe_3311r_LNP = cqD1_LNP*g1UV + (c2B_LNP*g1UV2)/3.;
468 Cqe_3322r_LNP = cqD1_LNP*g1UV + (c2B_LNP*g1UV2)/3.;
469 Cqe_3333r_LNP = cqD1_LNP*g1UV + (c2B_LNP*g1UV2)/3.;
470 Cqu1_1111r_LNP = (-2*c2B_LNP*g1UV2)/9.;
471 Cqu1_1122r_LNP = (-2*c2B_LNP*g1UV2)/9.;
472 Cqu1_1133r_LNP = -0.16666666666666666*(ctD_LNP*g1UV) - (2*c2B_LNP*g1UV2)/9.;
473 Cqu1_2211r_LNP = (-2*c2B_LNP*g1UV2)/9.;
474 Cqu1_2222r_LNP = (-2*c2B_LNP*g1UV2)/9.;
475 Cqu1_2233r_LNP = -0.16666666666666666*(ctD_LNP*g1UV) - (2*c2B_LNP*g1UV2)/9.;
476 Cqu1_3311r_LNP = (-2*cqD1_LNP*g1UV)/3. - (2*c2B_LNP*g1UV2)/9.;
477 Cqu1_3322r_LNP = (-2*cqD1_LNP*g1UV)/3. - (2*c2B_LNP*g1UV2)/9.;
478 Cqu1_3333r_LNP = cqt1_LNP - (2*cqD1_LNP*g1UV)/3. - (ctD_LNP*g1UV)/6. - (2*c2B_LNP*g1UV2)/9.;
479 Cqu8_1111r_LNP = -2*c2G_LNP*g3UV2;
480 Cqu8_1122r_LNP = -2*c2G_LNP*g3UV2;
481 Cqu8_1133r_LNP = -2*c2G_LNP*g3UV2;
482 Cqu8_2211r_LNP = -2*c2G_LNP*g3UV2;
483 Cqu8_2222r_LNP = -2*c2G_LNP*g3UV2;
484 Cqu8_2233r_LNP = -2*c2G_LNP*g3UV2;
485 Cqu8_3311r_LNP = -2*c2G_LNP*g3UV2;
486 Cqu8_3322r_LNP = -2*c2G_LNP*g3UV2;
487 Cqu8_3333r_LNP = cqt8_LNP - 2*c2G_LNP*g3UV2;
488 Cqd1_1111r_LNP = (c2B_LNP*g1UV2)/9.;
489 Cqd1_1122r_LNP = (c2B_LNP*g1UV2)/9.;
490 Cqd1_1133r_LNP = (c2B_LNP*g1UV2)/9.;
491 Cqd1_2211r_LNP = (c2B_LNP*g1UV2)/9.;
492 Cqd1_2222r_LNP = (c2B_LNP*g1UV2)/9.;
493 Cqd1_2233r_LNP = (c2B_LNP*g1UV2)/9.;
494 Cqd1_3311r_LNP = (cqD1_LNP*g1UV)/3. + (c2B_LNP*g1UV2)/9.;
495 Cqd1_3322r_LNP = (cqD1_LNP*g1UV)/3. + (c2B_LNP*g1UV2)/9.;
496 Cqd1_3333r_LNP = (cqD1_LNP*g1UV)/3. + (c2B_LNP*g1UV2)/9.;
497 Cqd8_1111r_LNP = -2*c2G_LNP*g3UV2;
498 Cqd8_1122r_LNP = -2*c2G_LNP*g3UV2;
499 Cqd8_1133r_LNP = -2*c2G_LNP*g3UV2;
500 Cqd8_2211r_LNP = -2*c2G_LNP*g3UV2;
501 Cqd8_2222r_LNP = -2*c2G_LNP*g3UV2;
502 Cqd8_2233r_LNP = -2*c2G_LNP*g3UV2;
503 Cqd8_3311r_LNP = -2*c2G_LNP*g3UV2;
504 Cqd8_3322r_LNP = -2*c2G_LNP*g3UV2;
505 Cqd8_3333r_LNP = -2*c2G_LNP*g3UV2;
506}
507
508
509bool NPd6SILH::PostUpdate()
510{
511 // 1) Post-update operations involving SM parameters only
512
513 v2LambdaNP2 = v() * v() / Lambda_NP/Lambda_NP;
514
515// Obtain the SM parameters at the UV scale
516 GenerateSMInitialConditions();
517
518// Perform the matching
519 setNPSMEFTd6GeneralParameters();
520
521 if (!NPSMEFTd6General::PostUpdate()) return (false);
522
523 return (true);
524
525}
526
527const double NPd6SILH::obliqueS() const
528{
529 return 0.0;
530}
531
532const double NPd6SILH::obliqueT() const
533{
534 return 0.0;
535}
536
537const double NPd6SILH::obliqueW() const
538{
539 return ( g2UV2 * c2W_LNP * v2LambdaNP2 / 2.0);
540}
541
542const double NPd6SILH::obliqueY() const
543{
544 return ( g2UV2 * c2B_LNP * v2LambdaNP2 / 2.0);
545}
546
547
548const double NPd6SILH::AuxObs_NP30() const
549{
550 // To be used for some temporary observable
551
552 // WY analysis at 13 TeV for HL-LHC 3/ab
553 double Wpar, Ypar, Wpar2, Ypar2;
554 double Chi2NC13, Chi2CC13, Chi2Tot;
555
556 Wpar = 10000.0 * obliqueW();
557 Ypar = 10000.0 * obliqueY();
558
559 Wpar2 = Wpar*Wpar;
560 Ypar2 = Ypar*Ypar;
561
562 Chi2CC13 = Wpar2 * (18.365037149441695 + 2.422904241798858 * Wpar + 0.12120594308623695 * Wpar2);
563
564 Chi2NC13 = 0.032772034538390675 * Wpar2 * Wpar2 + 2.815243944990361 * Ypar2 - 0.36522061776278516 * Ypar2 * Ypar
565 + 0.017375258924241194 * Ypar2 * Ypar2 + Wpar2 * Wpar * (-0.7059117582389635 + 0.006816297425306027 * Ypar)
566 + Wpar * Ypar * (7.988302197022343 + Ypar * (-0.5450119819316416 + 0.0050292149953719766 * Ypar))
567 + Wpar2 * (5.68581760491364 + Ypar * (-0.5794111075840261 + 0.048026245835369625 * Ypar));
568
569 Chi2Tot = Chi2CC13 + Chi2NC13;
570
571 // To be used as Gaussian observable with mean=0, var=1 I must return the sqrt.
572 return sqrt(Chi2Tot);
573}
Model::ModelParamMap
std::map< std::string, std::reference_wrapper< const double > > ModelParamMap
Definition: Model.h:280
Model::name
std::string name
The name of the model.
Definition: Model.h:285
Model::setModelName
void setModelName(const std::string name)
A method to set the name of the model.
Definition: Model.h:50
NPSMEFTd6General
A model class for new physics in the form of the dimension-six effective Lagrangian.
Definition: NPSMEFTd6General.h:593
NPSMEFTd6General::PostUpdate
virtual bool PostUpdate()
The post-update method for NPSMEFTd6General.
Definition: NPSMEFTd6General.cpp:8450
NPSMEFTd6General::setFlag
virtual bool setFlag(const std::string name, const bool value)
A method to check if all the mandatory parameters for NPSMEFTd6General have been provided in model in...
Definition: NPSMEFTd6General.cpp:14082
NPSMEFTd6General::GenerateSMInitialConditions
void GenerateSMInitialConditions()
Generates the initial condition for the Standard Model parameters.
Definition: NPSMEFTd6General.cpp:8355
NPd6SILH::cqq3_LNP
double cqq3_LNP
Definition: NPd6SILH.h:100
NPd6SILH::c3G_LNP
double c3G_LNP
Definition: NPd6SILH.h:88
NPd6SILH::c3W_LNP
double c3W_LNP
Definition: NPd6SILH.h:87
NPd6SILH::cgam_LNP
double cgam_LNP
Definition: NPd6SILH.h:91
NPd6SILH::obliqueT
virtual const double obliqueT() const
The oblique parameter . (Simplified implementation. Contribution only from .)
Definition: NPd6SILH.cpp:532
NPd6SILH::YuUV
gslpp::matrix< gslpp::complex > YuUV
Definition: NPd6SILH.h:126
NPd6SILH::ctt_LNP
double ctt_LNP
Definition: NPd6SILH.h:103
NPd6SILH::g1UV2
double g1UV2
Definition: NPd6SILH.h:125
NPd6SILH::YeUVhc
gslpp::matrix< gslpp::complex > YeUVhc
Yukawa matrices and h.c. at the UV scale.
Definition: NPd6SILH.h:126
NPd6SILH::NPd6SILHVars
static std::string NPd6SILHVars[NNPd6SILHVars]
Definition: NPd6SILH.h:19
NPd6SILH::cHq1_LNP
double cHq1_LNP
Definition: NPd6SILH.h:93
NPd6SILH::g2UV
double g2UV
Definition: NPd6SILH.h:124
NPd6SILH::g3UV
double g3UV
Definition: NPd6SILH.h:124
NPd6SILH::cg_LNP
double cg_LNP
Definition: NPd6SILH.h:92
NPd6SILH::ctB_LNP
double ctB_LNP
Definition: NPd6SILH.h:105
NPd6SILH::obliqueY
virtual const double obliqueY() const
The oblique parameter . (Simplified implementation. Contribution only from .)
Definition: NPd6SILH.cpp:542
NPd6SILH::cqD1_LNP
double cqD1_LNP
Definition: NPd6SILH.h:97
NPd6SILH::g3UV2
double g3UV2
SM gauge couplings at the UV scale (squared).
Definition: NPd6SILH.h:125
NPd6SILH::lambdaHUV
double lambdaHUV
SM couplings at the UV scale.
Definition: NPd6SILH.h:124
NPd6SILH::cu_LNP
double cu_LNP
Definition: NPd6SILH.h:107
NPd6SILH::v2LambdaNP2
double v2LambdaNP2
The ratio between the EW vev and the new physics scale, squared .
Definition: NPd6SILH.h:111
NPd6SILH::c2G_LNP
double c2G_LNP
Definition: NPd6SILH.h:86
NPd6SILH::obliqueW
virtual const double obliqueW() const
The oblique parameter . (Simplified implementation. Contribution only from .)
Definition: NPd6SILH.cpp:537
NPd6SILH::YuUVhc
gslpp::matrix< gslpp::complex > YuUVhc
Definition: NPd6SILH.h:126
NPd6SILH::cqt8_LNP
double cqt8_LNP
Definition: NPd6SILH.h:102
NPd6SILH::AuxObs_NP30
virtual const double AuxObs_NP30() const
Auxiliary observable AuxObs_NP30.
Definition: NPd6SILH.cpp:548
NPd6SILH::ctG_LNP
double ctG_LNP
Definition: NPd6SILH.h:104
NPd6SILH::obliqueS
virtual const double obliqueS() const
The oblique parameter . (Simplified implementation. Contribution only from .)
Definition: NPd6SILH.cpp:527
NPd6SILH::FlagRGEci
bool FlagRGEci
A boolean for the model flag RGEci , to include RGE effects. (Overwrittes NPSMEFTGeneral....
Definition: NPd6SILH.h:128
NPd6SILH::YdUV
gslpp::matrix< gslpp::complex > YdUV
Definition: NPd6SILH.h:126
NPd6SILH::cT_LNP
double cT_LNP
Definition: NPd6SILH.h:80
NPd6SILH::cHt_LNP
double cHt_LNP
Definition: NPd6SILH.h:95
NPd6SILH::cqt1_LNP
double cqt1_LNP
Definition: NPd6SILH.h:101
NPd6SILH::g1UV
double g1UV
Definition: NPd6SILH.h:124
NPd6SILH::cqq1_LNP
double cqq1_LNP
Definition: NPd6SILH.h:99
NPd6SILH::cqD3_LNP
double cqD3_LNP
Definition: NPd6SILH.h:98
NPd6SILH::setParameter
virtual void setParameter(const std::string name, const double &value)
Definition: NPd6SILH.cpp:58
NPd6SILH::cd_LNP
double cd_LNP
Definition: NPd6SILH.h:108
NPd6SILH::cW_LNP
double cW_LNP
Definition: NPd6SILH.h:83
NPd6SILH::ctW_LNP
double ctW_LNP
Definition: NPd6SILH.h:106
NPd6SILH::ce_LNP
double ce_LNP
Definition: NPd6SILH.h:109
NPd6SILH::cB_LNP
double cB_LNP
Definition: NPd6SILH.h:82
NPd6SILH::c6_LNP
double c6_LNP
Definition: NPd6SILH.h:81
NPd6SILH::cH_LNP
double cH_LNP
Definition: NPd6SILH.h:79
NPd6SILH::g2UV2
double g2UV2
Definition: NPd6SILH.h:125
NPd6SILH::YeUV
gslpp::matrix< gslpp::complex > YeUV
Definition: NPd6SILH.h:126
NPd6SILH::setNPSMEFTd6GeneralParameters
void setNPSMEFTd6GeneralParameters()
Definition: NPd6SILH.cpp:209
NPd6SILH::cHB_LNP
double cHB_LNP
Definition: NPd6SILH.h:90
NPd6SILH::NPd6SILH
NPd6SILH()
Definition: NPd6SILH.cpp:17
NPd6SILH::YdUVhc
gslpp::matrix< gslpp::complex > YdUVhc
Definition: NPd6SILH.h:126
NPd6SILH::c2B_LNP
double c2B_LNP
Definition: NPd6SILH.h:84
NPd6SILH::setFlag
virtual bool setFlag(const std::string name, const bool value)
Definition: NPd6SILH.cpp:195
NPd6SILH::cHW_LNP
double cHW_LNP
Definition: NPd6SILH.h:89
NPd6SILH::PostUpdate
virtual bool PostUpdate()
Definition: NPd6SILH.cpp:509
NPd6SILH::cHq3_LNP
double cHq3_LNP
Definition: NPd6SILH.h:94
NPd6SILH::ctD_LNP
double ctD_LNP
Definition: NPd6SILH.h:96
NPd6SILH::c2W_LNP
double c2W_LNP
Definition: NPd6SILH.h:85
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/src/StandardModel.cpp:280
StandardModel::v
const double v() const
The Higgs vacuum expectation value.
Definition: StandardModel/src/StandardModel.cpp:995