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
212 CG_LNP = c3G_LNP;
213 CW_LNP = c3W_LNP;
214 CH_LNP = -c6_LNP - 2*cHW_LNP*g2UV*lambdaHUV - 4*cW_LNP*g2UV*lambdaHUV - 2*c2W_LNP*g2UV2*lambdaHUV;
215 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.;
216 CHD_LNP = -4*cT_LNP - 2*cB_LNP*g1UV - cHB_LNP*g1UV - c2B_LNP*g1UV2;
217 CHG_LNP = cg_LNP;
218 CHB_LNP = cgam_LNP + (cHB_LNP*g1UV)/4.;
219 CHW_LNP = (cHW_LNP*g2UV)/4.;
220 CHWB_LNP = (cHW_LNP*g1UV)/4. + (cHB_LNP*g2UV)/4.;
221 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();
222 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();
223 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();
224 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();
225 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();
226 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();
227 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();
228 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();
229 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();
230 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();
231 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();
232 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();
233 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();
234 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();
235 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();
236 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();
237 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();
238 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();
239 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();
240 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();
241 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();
242 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();
243 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();
244 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();
245 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();
246 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();
247 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();
248 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();
249 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();
250 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();
251 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();
252 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();
253 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();
254 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();
255 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();
256 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();
257 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();
258 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();
259 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();
260 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();
261 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();
262 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();
263 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();
264 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();
265 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();
266 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();
267 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();
268 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();
269 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();
270 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();
271 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();
272 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();
273 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();
274 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();
275 CuG_33r_LNP = ctG_LNP;
276 CuG_33i_LNP = 0.0;
277 CuW_33r_LNP = ctW_LNP;
278 CuW_33i_LNP = 0.0;
279 CuB_33r_LNP = ctB_LNP;
280 CuB_33i_LNP = 0.0;
281 CHl1_11r_LNP = (cB_LNP*g1UV)/2. + (cHB_LNP*g1UV)/4. + (c2B_LNP*g1UV2)/2.;
282 CHl1_22r_LNP = (cB_LNP*g1UV)/2. + (cHB_LNP*g1UV)/4. + (c2B_LNP*g1UV2)/2.;
283 CHl1_33r_LNP = (cB_LNP*g1UV)/2. + (cHB_LNP*g1UV)/4. + (c2B_LNP*g1UV2)/2.;
284 CHl3_11r_LNP = -0.25*(cHW_LNP*g2UV) - (cW_LNP*g2UV)/2. - (c2W_LNP*g2UV2)/2.;
285 CHl3_22r_LNP = -0.25*(cHW_LNP*g2UV) - (cW_LNP*g2UV)/2. - (c2W_LNP*g2UV2)/2.;
286 CHl3_33r_LNP = -0.25*(cHW_LNP*g2UV) - (cW_LNP*g2UV)/2. - (c2W_LNP*g2UV2)/2.;
287 CHe_11r_LNP = cB_LNP*g1UV + (cHB_LNP*g1UV)/2. + c2B_LNP*g1UV2;
288 CHe_22r_LNP = cB_LNP*g1UV + (cHB_LNP*g1UV)/2. + c2B_LNP*g1UV2;
289 CHe_33r_LNP = cB_LNP*g1UV + (cHB_LNP*g1UV)/2. + c2B_LNP*g1UV2;
290 CHq1_11r_LNP = -0.16666666666666666*(cB_LNP*g1UV) - (cHB_LNP*g1UV)/12. - (c2B_LNP*g1UV2)/6.;
291 CHq1_22r_LNP = -0.16666666666666666*(cB_LNP*g1UV) - (cHB_LNP*g1UV)/12. - (c2B_LNP*g1UV2)/6.;
292 CHq1_33r_LNP = cHq1_LNP - (cB_LNP*g1UV)/6. - (cHB_LNP*g1UV)/12. - (cqD1_LNP*g1UV)/2. - (c2B_LNP*g1UV2)/6.;
293 CHq3_11r_LNP = -0.25*(cHW_LNP*g2UV) - (cW_LNP*g2UV)/2. - (c2W_LNP*g2UV2)/2.;
294 CHq3_22r_LNP = -0.25*(cHW_LNP*g2UV) - (cW_LNP*g2UV)/2. - (c2W_LNP*g2UV2)/2.;
295 CHq3_33r_LNP = cHq3_LNP - (cHW_LNP*g2UV)/4. - (cqD3_LNP*g2UV)/2. - (cW_LNP*g2UV)/2. - (c2W_LNP*g2UV2)/2.;
296 CHu_11r_LNP = (-2*cB_LNP*g1UV)/3. - (cHB_LNP*g1UV)/3. - (2*c2B_LNP*g1UV2)/3.;
297 CHu_22r_LNP = (-2*cB_LNP*g1UV)/3. - (cHB_LNP*g1UV)/3. - (2*c2B_LNP*g1UV2)/3.;
298 CHu_33r_LNP = cHt_LNP - (2*cB_LNP*g1UV)/3. - (cHB_LNP*g1UV)/3. - (ctD_LNP*g1UV)/2. - (2*c2B_LNP*g1UV2)/3.;
299 CHd_11r_LNP = (cB_LNP*g1UV)/3. + (cHB_LNP*g1UV)/6. + (c2B_LNP*g1UV2)/3.;
300 CHd_22r_LNP = (cB_LNP*g1UV)/3. + (cHB_LNP*g1UV)/6. + (c2B_LNP*g1UV2)/3.;
301 CHd_33r_LNP = (cB_LNP*g1UV)/3. + (cHB_LNP*g1UV)/6. + (c2B_LNP*g1UV2)/3.;
302 Cll_1111r_LNP = -0.25*(c2B_LNP*g1UV2) - (c2W_LNP*g2UV2)/4.;
303 Cll_1122r_LNP = -0.25*(c2B_LNP*g1UV2) + (c2W_LNP*g2UV2)/4.;
304 Cll_1133r_LNP = -0.25*(c2B_LNP*g1UV2) + (c2W_LNP*g2UV2)/4.;
305 Cll_1221r_LNP = -0.5*(c2W_LNP*g2UV2);
306 Cll_1331r_LNP = -0.5*(c2W_LNP*g2UV2);
307 Cll_2222r_LNP = -0.25*(c2B_LNP*g1UV2) - (c2W_LNP*g2UV2)/4.;
308 Cll_2233r_LNP = -0.25*(c2B_LNP*g1UV2) + (c2W_LNP*g2UV2)/4.;
309 Cll_2332r_LNP = -0.5*(c2W_LNP*g2UV2);
310 Cll_3333r_LNP = -0.25*(c2B_LNP*g1UV2) - (c2W_LNP*g2UV2)/4.;
311 Cqq1_1111r_LNP = -0.027777777777777776*(c2B_LNP*g1UV2) - (c2G_LNP*g3UV2)/12.;
312 Cqq1_1122r_LNP = -0.027777777777777776*(c2B_LNP*g1UV2) + (c2G_LNP*g3UV2)/6.;
313 Cqq1_1133r_LNP = -0.16666666666666666*(cqD1_LNP*g1UV) - (c2B_LNP*g1UV2)/36. + (c2G_LNP*g3UV2)/6.;
314 Cqq1_1221r_LNP = -0.25*(c2G_LNP*g3UV2);
315 Cqq1_1331r_LNP = -0.25*(c2G_LNP*g3UV2);
316 Cqq1_2222r_LNP = -0.027777777777777776*(c2B_LNP*g1UV2) - (c2G_LNP*g3UV2)/12.;
317 Cqq1_2233r_LNP = -0.16666666666666666*(cqD1_LNP*g1UV) - (c2B_LNP*g1UV2)/36. + (c2G_LNP*g3UV2)/6.;
318 Cqq1_2332r_LNP = -0.25*(c2G_LNP*g3UV2);
319 Cqq1_3333r_LNP = cqq1_LNP - (cqD1_LNP*g1UV)/6. - (c2B_LNP*g1UV2)/36. - (c2G_LNP*g3UV2)/12.;
320 Cqq3_1111r_LNP = -(c2W_LNP*g2UV2) - (c2G_LNP*g3UV2)/4.;
321 Cqq3_1122r_LNP = -(c2W_LNP*g2UV2);
322 Cqq3_1133r_LNP = -0.5*(cqD3_LNP*g2UV) - c2W_LNP*g2UV2;
323 Cqq3_1221r_LNP = -0.25*(c2G_LNP*g3UV2);
324 Cqq3_1331r_LNP = -0.25*(c2G_LNP*g3UV2);
325 Cqq3_2222r_LNP = -(c2W_LNP*g2UV2) - (c2G_LNP*g3UV2)/4.;
326 Cqq3_2233r_LNP = -0.5*(cqD3_LNP*g2UV) - c2W_LNP*g2UV2;
327 Cqq3_2332r_LNP = -0.25*(c2G_LNP*g3UV2);
328 Cqq3_3333r_LNP = cqq3_LNP - (cqD3_LNP*g2UV)/2. - c2W_LNP*g2UV2 - (c2G_LNP*g3UV2)/4.;
329 Clq1_1111r_LNP = (c2B_LNP*g1UV2)/6.;
330 Clq1_1122r_LNP = (c2B_LNP*g1UV2)/6.;
331 Clq1_1133r_LNP = (c2B_LNP*g1UV2)/6.;
332 Clq1_2211r_LNP = (c2B_LNP*g1UV2)/6.;
333 Clq1_2222r_LNP = (c2B_LNP*g1UV2)/6.;
334 Clq1_2233r_LNP = (c2B_LNP*g1UV2)/6.;
335 Clq1_3311r_LNP = (c2B_LNP*g1UV2)/6.;
336 Clq1_3322r_LNP = (c2B_LNP*g1UV2)/6.;
337 Clq1_3333r_LNP = (c2B_LNP*g1UV2)/6.;
338 Clq3_1111r_LNP = -0.5*(c2W_LNP*g2UV2);
339 Clq3_1122r_LNP = -0.5*(c2W_LNP*g2UV2);
340 Clq3_1133r_LNP = -0.5*(c2W_LNP*g2UV2);
341 Clq3_2211r_LNP = -0.5*(c2W_LNP*g2UV2);
342 Clq3_2222r_LNP = -0.5*(c2W_LNP*g2UV2);
343 Clq3_2233r_LNP = -0.5*(c2W_LNP*g2UV2);
344 Clq3_3311r_LNP = -0.5*(c2W_LNP*g2UV2);
345 Clq3_3322r_LNP = -0.5*(c2W_LNP*g2UV2);
346 Clq3_3333r_LNP = -0.5*(c2W_LNP*g2UV2);
347 Cee_1111r_LNP = -(c2B_LNP*g1UV2);
348 Cee_1122r_LNP = -(c2B_LNP*g1UV2);
349 Cee_1133r_LNP = -(c2B_LNP*g1UV2);
350 Cee_2222r_LNP = -(c2B_LNP*g1UV2);
351 Cee_2233r_LNP = -(c2B_LNP*g1UV2);
352 Cee_3333r_LNP = -(c2B_LNP*g1UV2);
353 Cuu_1111r_LNP = (-4*c2B_LNP*g1UV2)/9. - (c2G_LNP*g3UV2)/3.;
354 Cuu_1122r_LNP = (-4*c2B_LNP*g1UV2)/9. + (c2G_LNP*g3UV2)/6.;
355 Cuu_1133r_LNP = (-2*ctD_LNP*g1UV)/3. - (4*c2B_LNP*g1UV2)/9. + (c2G_LNP*g3UV2)/6.;
356 Cuu_1221r_LNP = -0.5*(c2G_LNP*g3UV2);
357 Cuu_1331r_LNP = -0.5*(c2G_LNP*g3UV2);
358 Cuu_2222r_LNP = (-4*c2B_LNP*g1UV2)/9. - (c2G_LNP*g3UV2)/3.;
359 Cuu_2233r_LNP = (-2*ctD_LNP*g1UV)/3. - (4*c2B_LNP*g1UV2)/9. + (c2G_LNP*g3UV2)/6.;
360 Cuu_2332r_LNP = -0.5*(c2G_LNP*g3UV2);
361 Cuu_3333r_LNP = ctt_LNP - (2*ctD_LNP*g1UV)/3. - (4*c2B_LNP*g1UV2)/9. - (c2G_LNP*g3UV2)/3.;
362 Cdd_1111r_LNP = -0.1111111111111111*(c2B_LNP*g1UV2) - (c2G_LNP*g3UV2)/3.;
363 Cdd_1122r_LNP = -0.1111111111111111*(c2B_LNP*g1UV2) + (c2G_LNP*g3UV2)/6.;
364 Cdd_1133r_LNP = -0.1111111111111111*(c2B_LNP*g1UV2) + (c2G_LNP*g3UV2)/6.;
365 Cdd_1221r_LNP = -0.5*(c2G_LNP*g3UV2);
366 Cdd_1331r_LNP = -0.5*(c2G_LNP*g3UV2);
367 Cdd_2222r_LNP = -0.1111111111111111*(c2B_LNP*g1UV2) - (c2G_LNP*g3UV2)/3.;
368 Cdd_2233r_LNP = -0.1111111111111111*(c2B_LNP*g1UV2) + (c2G_LNP*g3UV2)/6.;
369 Cdd_2332r_LNP = -0.5*(c2G_LNP*g3UV2);
370 Cdd_3333r_LNP = -0.1111111111111111*(c2B_LNP*g1UV2) - (c2G_LNP*g3UV2)/3.;
371 Ceu_1111r_LNP = (4*c2B_LNP*g1UV2)/3.;
372 Ceu_1122r_LNP = (4*c2B_LNP*g1UV2)/3.;
373 Ceu_1133r_LNP = ctD_LNP*g1UV + (4*c2B_LNP*g1UV2)/3.;
374 Ceu_2211r_LNP = (4*c2B_LNP*g1UV2)/3.;
375 Ceu_2222r_LNP = (4*c2B_LNP*g1UV2)/3.;
376 Ceu_2233r_LNP = ctD_LNP*g1UV + (4*c2B_LNP*g1UV2)/3.;
377 Ceu_3311r_LNP = (4*c2B_LNP*g1UV2)/3.;
378 Ceu_3322r_LNP = (4*c2B_LNP*g1UV2)/3.;
379 Ceu_3333r_LNP = ctD_LNP*g1UV + (4*c2B_LNP*g1UV2)/3.;
380 Ced_1111r_LNP = (-2*c2B_LNP*g1UV2)/3.;
381 Ced_1122r_LNP = (-2*c2B_LNP*g1UV2)/3.;
382 Ced_1133r_LNP = (-2*c2B_LNP*g1UV2)/3.;
383 Ced_2211r_LNP = (-2*c2B_LNP*g1UV2)/3.;
384 Ced_2222r_LNP = (-2*c2B_LNP*g1UV2)/3.;
385 Ced_2233r_LNP = (-2*c2B_LNP*g1UV2)/3.;
386 Ced_3311r_LNP = (-2*c2B_LNP*g1UV2)/3.;
387 Ced_3322r_LNP = (-2*c2B_LNP*g1UV2)/3.;
388 Ced_3333r_LNP = (-2*c2B_LNP*g1UV2)/3.;
389 Cud1_1111r_LNP = (4*c2B_LNP*g1UV2)/9.;
390 Cud1_1122r_LNP = (4*c2B_LNP*g1UV2)/9.;
391 Cud1_1133r_LNP = (4*c2B_LNP*g1UV2)/9.;
392 Cud1_2211r_LNP = (4*c2B_LNP*g1UV2)/9.;
393 Cud1_2222r_LNP = (4*c2B_LNP*g1UV2)/9.;
394 Cud1_2233r_LNP = (4*c2B_LNP*g1UV2)/9.;
395 Cud1_3311r_LNP = (ctD_LNP*g1UV)/3. + (4*c2B_LNP*g1UV2)/9.;
396 Cud1_3322r_LNP = (ctD_LNP*g1UV)/3. + (4*c2B_LNP*g1UV2)/9.;
397 Cud1_3333r_LNP = (ctD_LNP*g1UV)/3. + (4*c2B_LNP*g1UV2)/9.;
398 Cud8_1111r_LNP = -2*c2G_LNP*g3UV2;
399 Cud8_1122r_LNP = -2*c2G_LNP*g3UV2;
400 Cud8_1133r_LNP = -2*c2G_LNP*g3UV2;
401 Cud8_2211r_LNP = -2*c2G_LNP*g3UV2;
402 Cud8_2222r_LNP = -2*c2G_LNP*g3UV2;
403 Cud8_2233r_LNP = -2*c2G_LNP*g3UV2;
404 Cud8_3311r_LNP = -2*c2G_LNP*g3UV2;
405 Cud8_3322r_LNP = -2*c2G_LNP*g3UV2;
406 Cud8_3333r_LNP = -2*c2G_LNP*g3UV2;
407 Cle_1111r_LNP = -(c2B_LNP*g1UV2);
408 Cle_1122r_LNP = -(c2B_LNP*g1UV2);
409 Cle_1133r_LNP = -(c2B_LNP*g1UV2);
410 Cle_2211r_LNP = -(c2B_LNP*g1UV2);
411 Cle_2222r_LNP = -(c2B_LNP*g1UV2);
412 Cle_2233r_LNP = -(c2B_LNP*g1UV2);
413 Cle_3311r_LNP = -(c2B_LNP*g1UV2);
414 Cle_3322r_LNP = -(c2B_LNP*g1UV2);
415 Cle_3333r_LNP = -(c2B_LNP*g1UV2);
416 Clu_1111r_LNP = (2*c2B_LNP*g1UV2)/3.;
417 Clu_1122r_LNP = (2*c2B_LNP*g1UV2)/3.;
418 Clu_1133r_LNP = (ctD_LNP*g1UV)/2. + (2*c2B_LNP*g1UV2)/3.;
419 Clu_2211r_LNP = (2*c2B_LNP*g1UV2)/3.;
420 Clu_2222r_LNP = (2*c2B_LNP*g1UV2)/3.;
421 Clu_2233r_LNP = (ctD_LNP*g1UV)/2. + (2*c2B_LNP*g1UV2)/3.;
422 Clu_3311r_LNP = (2*c2B_LNP*g1UV2)/3.;
423 Clu_3322r_LNP = (2*c2B_LNP*g1UV2)/3.;
424 Clu_3333r_LNP = (ctD_LNP*g1UV)/2. + (2*c2B_LNP*g1UV2)/3.;
425 Cld_1111r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
426 Cld_1122r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
427 Cld_1133r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
428 Cld_2211r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
429 Cld_2222r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
430 Cld_2233r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
431 Cld_3311r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
432 Cld_3322r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
433 Cld_3333r_LNP = -0.3333333333333333*(c2B_LNP*g1UV2);
434 Cqe_1111r_LNP = (c2B_LNP*g1UV2)/3.;
435 Cqe_1122r_LNP = (c2B_LNP*g1UV2)/3.;
436 Cqe_1133r_LNP = (c2B_LNP*g1UV2)/3.;
437 Cqe_2211r_LNP = (c2B_LNP*g1UV2)/3.;
438 Cqe_2222r_LNP = (c2B_LNP*g1UV2)/3.;
439 Cqe_2233r_LNP = (c2B_LNP*g1UV2)/3.;
440 Cqe_3311r_LNP = cqD1_LNP*g1UV + (c2B_LNP*g1UV2)/3.;
441 Cqe_3322r_LNP = cqD1_LNP*g1UV + (c2B_LNP*g1UV2)/3.;
442 Cqe_3333r_LNP = cqD1_LNP*g1UV + (c2B_LNP*g1UV2)/3.;
443 Cqu1_1111r_LNP = (-2*c2B_LNP*g1UV2)/9.;
444 Cqu1_1122r_LNP = (-2*c2B_LNP*g1UV2)/9.;
445 Cqu1_1133r_LNP = -0.16666666666666666*(ctD_LNP*g1UV) - (2*c2B_LNP*g1UV2)/9.;
446 Cqu1_2211r_LNP = (-2*c2B_LNP*g1UV2)/9.;
447 Cqu1_2222r_LNP = (-2*c2B_LNP*g1UV2)/9.;
448 Cqu1_2233r_LNP = -0.16666666666666666*(ctD_LNP*g1UV) - (2*c2B_LNP*g1UV2)/9.;
449 Cqu1_3311r_LNP = (-2*cqD1_LNP*g1UV)/3. - (2*c2B_LNP*g1UV2)/9.;
450 Cqu1_3322r_LNP = (-2*cqD1_LNP*g1UV)/3. - (2*c2B_LNP*g1UV2)/9.;
451 Cqu1_3333r_LNP = cqt1_LNP - (2*cqD1_LNP*g1UV)/3. - (ctD_LNP*g1UV)/6. - (2*c2B_LNP*g1UV2)/9.;
452 Cqu8_1111r_LNP = -2*c2G_LNP*g3UV2;
453 Cqu8_1122r_LNP = -2*c2G_LNP*g3UV2;
454 Cqu8_1133r_LNP = -2*c2G_LNP*g3UV2;
455 Cqu8_2211r_LNP = -2*c2G_LNP*g3UV2;
456 Cqu8_2222r_LNP = -2*c2G_LNP*g3UV2;
457 Cqu8_2233r_LNP = -2*c2G_LNP*g3UV2;
458 Cqu8_3311r_LNP = -2*c2G_LNP*g3UV2;
459 Cqu8_3322r_LNP = -2*c2G_LNP*g3UV2;
460 Cqu8_3333r_LNP = cqt8_LNP - 2*c2G_LNP*g3UV2;
461 Cqd1_1111r_LNP = (c2B_LNP*g1UV2)/9.;
462 Cqd1_1122r_LNP = (c2B_LNP*g1UV2)/9.;
463 Cqd1_1133r_LNP = (c2B_LNP*g1UV2)/9.;
464 Cqd1_2211r_LNP = (c2B_LNP*g1UV2)/9.;
465 Cqd1_2222r_LNP = (c2B_LNP*g1UV2)/9.;
466 Cqd1_2233r_LNP = (c2B_LNP*g1UV2)/9.;
467 Cqd1_3311r_LNP = (cqD1_LNP*g1UV)/3. + (c2B_LNP*g1UV2)/9.;
468 Cqd1_3322r_LNP = (cqD1_LNP*g1UV)/3. + (c2B_LNP*g1UV2)/9.;
469 Cqd1_3333r_LNP = (cqD1_LNP*g1UV)/3. + (c2B_LNP*g1UV2)/9.;
470 Cqd8_1111r_LNP = -2*c2G_LNP*g3UV2;
471 Cqd8_1122r_LNP = -2*c2G_LNP*g3UV2;
472 Cqd8_1133r_LNP = -2*c2G_LNP*g3UV2;
473 Cqd8_2211r_LNP = -2*c2G_LNP*g3UV2;
474 Cqd8_2222r_LNP = -2*c2G_LNP*g3UV2;
475 Cqd8_2233r_LNP = -2*c2G_LNP*g3UV2;
476 Cqd8_3311r_LNP = -2*c2G_LNP*g3UV2;
477 Cqd8_3322r_LNP = -2*c2G_LNP*g3UV2;
478 Cqd8_3333r_LNP = -2*c2G_LNP*g3UV2;
479}
480
481
482bool NPd6SILH::PostUpdate()
483{
484
485// Obtain the values of the SM parameters at the UV scale, to do the matching
486 ChangeToEvolutorsBasisPureSM();
487 double Mu_LEW[3] = {mu_LEW, mc_LEW, mt_LEW};
488 double Md_LEW[3] = {md_LEW, ms_LEW, mb_LEW};
489 double Me_LEW[3] = {me_LEW, mmu_LEW, mtau_LEW};
490
491 if (FlagRGEci) {
492 // SM initial conditions at the UV scale. Use RGEsolver SMEFTEvolEW
493 SMEFTEvolEW.GenerateSMInitialConditions(muw, Lambda_NP, SMEFTBasisFlag, "Numeric",
494 g1_LEW, g2_LEW, g3_LEW, lambdaH_LEW, mH2_LEW,
495 Mu_LEW, Md_LEW, Me_LEW, s12CKM_LEW, s13CKM_LEW, s23CKM_LEW, dCKM_LEW);
496
497 } else {
498 // SM initial conditions at the UV scale. Use RGEsolver SMEFTEvolEW
499 // Skip RGE by setting the two scales at Lambda_NP for the EFT and to muw for the SM pars
500 SMEFTEvolEW.GenerateSMInitialConditions(muw, muw, SMEFTBasisFlag, "Numeric",
501 g1_LEW, g2_LEW, g3_LEW, lambdaH_LEW, mH2_LEW,
502 Mu_LEW, Md_LEW, Me_LEW, s12CKM_LEW, s13CKM_LEW, s23CKM_LEW, dCKM_LEW);
503 }
504
505// SM parameters at the UV scale
506 g1UV = getSMEFTCoeffEW("g1");
507 g2UV = getSMEFTCoeffEW("g2");
508 g3UV = getSMEFTCoeffEW("g3");
509 lambdaHUV = getSMEFTCoeffEW("lambda");
510
511 g1UV2 = g1UV * g1UV;
512 g2UV2 = g2UV * g2UV;
513 g3UV2 = g3UV * g3UV;
514
515 for (int i = 0; i < 3; i++) {
516 for (int j = 0; j < 3; j++) {
517 YuUV.assignre(i, j, getSMEFTCoeffEW("YuR", i, j) );
518 YuUV.assignim(i, j, getSMEFTCoeffEW("YuI", i, j) );
519 YdUV.assignre(i, j, getSMEFTCoeffEW("YdR", i, j) );
520 YdUV.assignim(i, j, getSMEFTCoeffEW("YdI", i, j) );
521 YeUV.assignre(i, j, getSMEFTCoeffEW("YeR", i, j) );
522 YeUV.assignim(i, j, getSMEFTCoeffEW("YeI", i, j) );
523 }
524 }
525
526 YuUVhc = YuUV.hconjugate();
527 YdUVhc = YdUV.hconjugate();
528 YeUVhc = YeUV.hconjugate();
529
530
531 setNPSMEFTd6GeneralParameters();
532
533 if (!NPSMEFTd6General::PostUpdate()) return (false);
534
535 return (true);
536
537}
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::Cqu8_3333r_LNP
double Cqu8_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7362
NPSMEFTd6General::CW_LNP
double CW_LNP
The dimension-6 operator coefficient .
Definition: NPSMEFTd6General.h:7157
NPSMEFTd6General::CHbox_LNP
double CHbox_LNP
The dimension-6 operator coefficient .
Definition: NPSMEFTd6General.h:7165
NPSMEFTd6General::CHu_33r_LNP
double CHu_33r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7183
NPSMEFTd6General::CeH_33i_LNP
double CeH_33i_LNP
The dimension-6 operator coefficient (Imaginary part).
Definition: NPSMEFTd6General.h:7190
NPSMEFTd6General::PostUpdate
virtual bool PostUpdate()
The post-update method for NPSMEFTd6General.
Definition: NPSMEFTd6General.cpp:8354
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:14054
NPSMEFTd6General::Cqe_3333r_LNP
double Cqe_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7288
NPSMEFTd6General::Cqq1_3333r_LNP
double Cqq1_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7313
NPSMEFTd6General::CdH_33i_LNP
double CdH_33i_LNP
The dimension-6 operator coefficient (Imaginary part).
Definition: NPSMEFTd6General.h:7194
NPSMEFTd6General::Cqd8_3333r_LNP
double Cqd8_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7380
NPSMEFTd6General::setParameter
virtual void setParameter(const std::string name, const double &value)
A method to set the value of a parameter of the model.
Definition: NPSMEFTd6General.cpp:8870
NPSMEFTd6General::CHl1_33r_LNP
double CHl1_33r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7173
NPSMEFTd6General::Cld_3333r_LNP
double Cld_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7279
NPSMEFTd6General::Ced_3333r_LNP
double Ced_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7252
NPSMEFTd6General::Ceu_3333r_LNP
double Ceu_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7243
NPSMEFTd6General::CuW_33r_LNP
double CuW_33r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7197
NPSMEFTd6General::ChangeToEvolutorsBasisPureSM
void ChangeToEvolutorsBasisPureSM()
This function computes the SM parameters needed for the evolutor, neglecting any NP contribution.
Definition: NPSMEFTd6General.cpp:3024
NPSMEFTd6General::Clq3_3333r_LNP
double Clq3_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7229
NPSMEFTd6General::Cqd1_3333r_LNP
double Cqd1_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7371
NPSMEFTd6General::CdH_33r_LNP
double CdH_33r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7193
NPSMEFTd6General::CHq1_33r_LNP
double CHq1_33r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7179
NPSMEFTd6General::CHG_LNP
double CHG_LNP
The dimension-6 operator coefficient .
Definition: NPSMEFTd6General.h:7158
NPSMEFTd6General::Cll_3333r_LNP
double Cll_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7213
NPSMEFTd6General::Cqu1_3333r_LNP
double Cqu1_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7353
NPSMEFTd6General::CHD_LNP
double CHD_LNP
The dimension-6 operator coefficient .
Definition: NPSMEFTd6General.h:7164
NPSMEFTd6General::CHd_33r_LNP
double CHd_33r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7185
NPSMEFTd6General::Cqq3_3333r_LNP
double Cqq3_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7318
NPSMEFTd6General::CHB_LNP
double CHB_LNP
The dimension-6 operator coefficient .
Definition: NPSMEFTd6General.h:7160
NPSMEFTd6General::CHe_33r_LNP
double CHe_33r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7177
NPSMEFTd6General::CHWB_LNP
double CHWB_LNP
The dimension-6 operator coefficient .
Definition: NPSMEFTd6General.h:7163
NPSMEFTd6General::Cle_3333r_LNP
double Cle_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7261
NPSMEFTd6General::CuG_33i_LNP
double CuG_33i_LNP
The dimension-6 operator coefficient (Imaginary part).
Definition: NPSMEFTd6General.h:7196
NPSMEFTd6General::Clq1_3333r_LNP
double Clq1_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7220
NPSMEFTd6General::CHq3_33r_LNP
double CHq3_33r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7181
NPSMEFTd6General::SMEFTBasisFlag
std::string SMEFTBasisFlag
Definition: NPSMEFTd6General.h:7114
NPSMEFTd6General::CuG_33r_LNP
double CuG_33r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7195
NPSMEFTd6General::CuH_33i_LNP
double CuH_33i_LNP
The dimension-6 operator coefficient (Imaginary part).
Definition: NPSMEFTd6General.h:7192
NPSMEFTd6General::CeH_33r_LNP
double CeH_33r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7189
NPSMEFTd6General::CG_LNP
double CG_LNP
The dimension-6 operator coefficient .
Definition: NPSMEFTd6General.h:7156
NPSMEFTd6General::CuB_33r_LNP
double CuB_33r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7199
NPSMEFTd6General::Clu_3333r_LNP
double Clu_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7270
NPSMEFTd6General::Cud1_3333r_LNP
double Cud1_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7335
NPSMEFTd6General::Cuu_3333r_LNP
double Cuu_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7323
NPSMEFTd6General::CHW_LNP
double CHW_LNP
The dimension-6 operator coefficient .
Definition: NPSMEFTd6General.h:7159
NPSMEFTd6General::CuB_33i_LNP
double CuB_33i_LNP
The dimension-6 operator coefficient (Imaginary part).
Definition: NPSMEFTd6General.h:7200
NPSMEFTd6General::Cud8_3333r_LNP
double Cud8_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7344
NPSMEFTd6General::CH_LNP
double CH_LNP
The dimension-6 operator coefficient .
Definition: NPSMEFTd6General.h:7166
NPSMEFTd6General::Cee_3333r_LNP
double Cee_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7236
NPSMEFTd6General::CuW_33i_LNP
double CuW_33i_LNP
The dimension-6 operator coefficient (Imaginary part).
Definition: NPSMEFTd6General.h:7198
NPSMEFTd6General::getSMEFTCoeffEW
double getSMEFTCoeffEW(const std::string name) const
Definition: NPSMEFTd6General.h:6790
NPSMEFTd6General::CuH_33r_LNP
double CuH_33r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7191
NPSMEFTd6General::CHl3_33r_LNP
double CHl3_33r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7175
NPSMEFTd6General::Lambda_NP
double Lambda_NP
The new physics scale [GeV].
Definition: NPSMEFTd6General.h:7458
NPSMEFTd6General::Cdd_3333r_LNP
double Cdd_3333r_LNP
The dimension-6 operator coefficient (Real part and pure real operator).
Definition: NPSMEFTd6General.h:7328
NPd6SILH::cqq3_LNP
double cqq3_LNP
Definition: NPd6SILH.h:65
NPd6SILH::c3G_LNP
double c3G_LNP
Definition: NPd6SILH.h:53
NPd6SILH::c3W_LNP
double c3W_LNP
Definition: NPd6SILH.h:52
NPd6SILH::cgam_LNP
double cgam_LNP
Definition: NPd6SILH.h:56
NPd6SILH::YuUV
gslpp::matrix< gslpp::complex > YuUV
Definition: NPd6SILH.h:90
NPd6SILH::ctt_LNP
double ctt_LNP
Definition: NPd6SILH.h:68
NPd6SILH::g1UV2
double g1UV2
Definition: NPd6SILH.h:89
NPd6SILH::YeUVhc
gslpp::matrix< gslpp::complex > YeUVhc
Yukawa matrices and h.c. at the UV scale.
Definition: NPd6SILH.h:90
NPd6SILH::NPd6SILHVars
static std::string NPd6SILHVars[NNPd6SILHVars]
Definition: NPd6SILH.h:19
NPd6SILH::cHq1_LNP
double cHq1_LNP
Definition: NPd6SILH.h:58
NPd6SILH::g2UV
double g2UV
Definition: NPd6SILH.h:88
NPd6SILH::g3UV
double g3UV
Definition: NPd6SILH.h:88
NPd6SILH::cg_LNP
double cg_LNP
Definition: NPd6SILH.h:57
NPd6SILH::ctB_LNP
double ctB_LNP
Definition: NPd6SILH.h:70
NPd6SILH::cqD1_LNP
double cqD1_LNP
Definition: NPd6SILH.h:62
NPd6SILH::g3UV2
double g3UV2
SM gauge couplings at the UV scale (squared).
Definition: NPd6SILH.h:89
NPd6SILH::lambdaHUV
double lambdaHUV
SM couplings at the UV scale.
Definition: NPd6SILH.h:88
NPd6SILH::cu_LNP
double cu_LNP
Definition: NPd6SILH.h:72
NPd6SILH::c2G_LNP
double c2G_LNP
Definition: NPd6SILH.h:51
NPd6SILH::YuUVhc
gslpp::matrix< gslpp::complex > YuUVhc
Definition: NPd6SILH.h:90
NPd6SILH::cqt8_LNP
double cqt8_LNP
Definition: NPd6SILH.h:67
NPd6SILH::ctG_LNP
double ctG_LNP
Definition: NPd6SILH.h:69
NPd6SILH::FlagRGEci
bool FlagRGEci
A boolean for the model flag RGEci , to include RGE effects. (Overwrittes NPSMEFTGeneral....
Definition: NPd6SILH.h:92
NPd6SILH::YdUV
gslpp::matrix< gslpp::complex > YdUV
Definition: NPd6SILH.h:90
NPd6SILH::cT_LNP
double cT_LNP
Definition: NPd6SILH.h:45
NPd6SILH::cHt_LNP
double cHt_LNP
Definition: NPd6SILH.h:60
NPd6SILH::cqt1_LNP
double cqt1_LNP
Definition: NPd6SILH.h:66
NPd6SILH::g1UV
double g1UV
Definition: NPd6SILH.h:88
NPd6SILH::cqq1_LNP
double cqq1_LNP
Definition: NPd6SILH.h:64
NPd6SILH::cqD3_LNP
double cqD3_LNP
Definition: NPd6SILH.h:63
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:73
NPd6SILH::cW_LNP
double cW_LNP
Definition: NPd6SILH.h:48
NPd6SILH::ctW_LNP
double ctW_LNP
Definition: NPd6SILH.h:71
NPd6SILH::ce_LNP
double ce_LNP
Definition: NPd6SILH.h:74
NPd6SILH::cB_LNP
double cB_LNP
Definition: NPd6SILH.h:47
NPd6SILH::c6_LNP
double c6_LNP
Definition: NPd6SILH.h:46
NPd6SILH::cH_LNP
double cH_LNP
Definition: NPd6SILH.h:44
NPd6SILH::g2UV2
double g2UV2
Definition: NPd6SILH.h:89
NPd6SILH::YeUV
gslpp::matrix< gslpp::complex > YeUV
Definition: NPd6SILH.h:90
NPd6SILH::setNPSMEFTd6GeneralParameters
void setNPSMEFTd6GeneralParameters()
Definition: NPd6SILH.cpp:209
NPd6SILH::cHB_LNP
double cHB_LNP
Definition: NPd6SILH.h:55
NPd6SILH::NPd6SILH
NPd6SILH()
Definition: NPd6SILH.cpp:17
NPd6SILH::YdUVhc
gslpp::matrix< gslpp::complex > YdUVhc
Definition: NPd6SILH.h:90
NPd6SILH::c2B_LNP
double c2B_LNP
Definition: NPd6SILH.h:49
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:54
NPd6SILH::PostUpdate
virtual bool PostUpdate()
Definition: NPd6SILH.cpp:482
NPd6SILH::cHq3_LNP
double cHq3_LNP
Definition: NPd6SILH.h:59
NPd6SILH::ctD_LNP
double ctD_LNP
Definition: NPd6SILH.h:61
NPd6SILH::c2W_LNP
double c2W_LNP
Definition: NPd6SILH.h:50
StandardModel::muw
double muw
A matching scale around the weak scale in GeV.
Definition: StandardModel.h:3470