54 double delta_vT =
mySMEFT.getDelta_v();
55 double vTosq2 = vT / sqrt(2.);
103 double cbar =
mySMEFT.getXWZ_tree();
104 double sbar = -
mySMEFT.getXBZ_tree();
105 double sbar2 = sbar * sbar;
107 double delta_sbar =
mySMEFT.getDelta_xBZ();
108 double g1bar =
mySMEFT.getG1_tree();
110 double g2bar =
mySMEFT.getG2_tree();
112 double delta_MZ2 =
mySMEFT.getDelta_Mz2();
113 double ebar =
mySMEFT.getEeMz();
116 double gZbar = ebar / sbar / cbar;
117 double delta_gZbar = (g1bar * g1bar + g2bar * g2bar) / (2. * g1bar * g2bar) *
v2 *
mySMEFT.getSMEFTCoeffEW(
"CHWB");
119 delta_gZbar += -(g1bar * g1bar * g1bar * g1bar + g2bar * g2bar * g2bar * g2bar) / (2. * (g1bar * g1bar + g2bar * g2bar) * g1bar * g2bar) *
v2 *
mySMEFT.getSMEFTCoeffEW(
"CHWB");
122 double delta_gZ2oMZ2 = 2. * delta_gZbar - delta_MZ2;
123 double g22oMW2 = 4. /
v2;
124 double delta_g22oMW2 = -2. * delta_vT;
129 double deltagZbar2oMZ4 = 2. * delta_gZbar - 2. * delta_MZ2;
133 double deltag2bar2oMW4 = -4. * delta_vT;
135 double lambda =
mySMEFT.getLambdaH_tree() * 2.;
136 double oneoMh2 =
mySMEFT.FlagHiggsExch * 1. / (lambda *
v2);
137 double deltaoneoMh2 = (- 2. *
mySMEFT.getSMEFTCoeffEW(
"CHbox") + 0.5 *
mySMEFT.getSMEFTCoeffEW(
"CHD") + 3 *
mySMEFT.getSMEFTCoeffEW(
"CH") / lambda) *
v2;
138 std::array<double, 3> Me =
mySMEFT.getMe_LEW();
139 std::array<double, 3> Mu =
mySMEFT.getMu_LEW();
140 std::array<double, 3> Md =
mySMEFT.getMd_LEW();
151 gslpp::matrix<complex> VCKMd =
VCKM.hconjugate();
155 if (
Ceg.at(0).at(1) != 0. ||
Ceg.at(0).at(2) != 0. ||
Ceg.at(1).at(0) != 0. ||
Ceg.at(1).at(2) != 0. ||
Ceg.at(2).at(0) != 0. ||
Ceg.at(2).at(1) != 0.)
156 throw(
"Compiler is not putting to zero correctly the 2-d arrays of Wilson coefficients");
219#ifdef NOLEPTONFLAVOURVIOLATION
223 for (
int i = 0; i < 3; i++)
224 for (
int j = 0; j < 3; j++)
226 Ceg.at(i).at(j) += vTosq2 * (-(
mySMEFT.getSMEFTCoeffEW(
"CeWR", i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", i, j)) * sbar + (
mySMEFT.getSMEFTCoeffEW(
"CeBR", i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", i, j)) * cbar);
228 CnunuVLL.at(i).at(i).at(j).at(j) += -0.0625 * (delta_gZ2oMZ2 * gZ2oMZ2);
229 CnunuVLL.at(i).at(j).at(j).at(i) += -0.0625 * (delta_gZ2oMZ2 * gZ2oMZ2);
230 CeeVLL.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * (-1 + 2 * sbar2) * (delta_gZ2oMZ2 - 2 * delta_gZ2oMZ2 * sbar2 - 8 * delta_sbar * sbar2)) / 16.;
231 CeeVLL.at(i).at(j).at(j).at(i) += (gZ2oMZ2 * (-1 + 2 * sbar2) * (delta_gZ2oMZ2 - 2 * delta_gZ2oMZ2 * sbar2 - 8 * delta_sbar * sbar2)) / 16.;
232 CnueVLL.at(i).at(i).at(j).at(j) += -0.25 * (gZ2oMZ2 * (4 * delta_sbar * sbar2 + delta_gZ2oMZ2 * (-1 + 2 * sbar2)));
233 CnueVLL.at(i).at(j).at(j).at(i) += -0.5 * (delta_g22oMW2 * g22oMW2);
234 CnudVLL.at(i).at(i).at(j).at(j) += -0.08333333333333333 * (gZ2oMZ2 * (4 * delta_sbar * sbar2 + delta_gZ2oMZ2 * (-3 + 2 * sbar2)));
235 CedVLL.at(i).at(i).at(j).at(j) += -0.08333333333333333 * (gZ2oMZ2 * (16 * delta_sbar * (-1 + sbar2) * sbar2 + delta_gZ2oMZ2 * (3 - 8 * sbar2 + 4 * sbar2 * sbar2)));
236 CddVLL.at(i).at(i).at(j).at(j) += -0.013888888888888888 * (gZ2oMZ2 * (-3 + 2 * sbar2) * (8 * delta_sbar * sbar2 + delta_gZ2oMZ2 * (-3 + 2 * sbar2)));
237 CeeVRR.at(i).at(i).at(j).at(j) += -0.25 * ((delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2);
238 CeeVRR.at(i).at(j).at(j).at(i) += -0.25 * ((delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2);
239 CedVRR.at(i).at(i).at(j).at(j) += -0.3333333333333333 * ((delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2);
240 CddVRR.at(i).at(i).at(j).at(j) += -0.05555555555555555 * ((delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2);
241 CnueVLR.at(i).at(i).at(j).at(j) += -0.5 * ((delta_gZ2oMZ2 + 2 * delta_sbar) * gZ2oMZ2 * sbar2);
242 CeeVLR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * sbar2 * (delta_gZ2oMZ2 + delta_sbar * (2 - 8 * sbar2) - 2 * delta_gZ2oMZ2 * sbar2)) / 2.;
243 CnudVLR.at(i).at(i).at(j).at(j) += -0.16666666666666666 * ((delta_gZ2oMZ2 + 2 * delta_sbar) * gZ2oMZ2 * sbar2);
244 CedVLR.at(i).at(i).at(j).at(j) += -0.16666666666666666 * (gZ2oMZ2 * sbar2 * (-delta_gZ2oMZ2 - 2 * delta_sbar + 2 * delta_gZ2oMZ2 * sbar2 + 8 * delta_sbar * sbar2));
245 CdeVLR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * (delta_sbar * (6 - 8 * sbar2) + delta_gZ2oMZ2 * (3 - 2 * sbar2)) * sbar2) / 6.;
246 CddV1LR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * (delta_sbar * (6 - 8 * sbar2) + delta_gZ2oMZ2 * (3 - 2 * sbar2)) * sbar2) / 18.;
248 CnueVLR.at(i).at(j).at(j).at(i) += -0.25*((1 + deltag2bar2oMW4)*g2bar2oMW4*Me[i]*Me[j]);
249 CeeVLR.at(i).at(j).at(j).at(i) += -0.125*(gZbar2oMZ4*(-1 + 2*sbar2)*(-1 + (2 + 8*delta_sbar)*sbar2 + deltagZbar2oMZ4*(-1 + 2*sbar2))*Me[i]*Me[j]);
250 CeeVLR.at(i).at(j).at(j).at(i) += (gZbar2oMZ4*sbar2*(-1 + 2*sbar2 + deltagZbar2oMZ4*(-1 + 2*sbar2) + delta_sbar*(-2 + 8*sbar2))*Me[i]*Me[j])/4.;
251 CeeVLR.at(i).at(j).at(j).at(i) += (gZbar2oMZ4*sbar2*(-1 + 2*sbar2 + deltagZbar2oMZ4*(-1 + 2*sbar2) + delta_sbar*(-2 + 8*sbar2))*Me[i]*Me[j])/4.;
252 CeeVLR.at(i).at(j).at(j).at(i) += -0.5*((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Me[i]*Me[j]);
253 CddV1LR.at(i).at(j).at(j).at(i) += -0.004629629629629629*(gZbar2oMZ4*(-3 + 2*sbar2)*(-3 + (2 + 8*delta_sbar)*sbar2 + deltagZbar2oMZ4*(-3 + 2*sbar2))*Md[i]*Md[j]);
254 CddV1LR.at(i).at(j).at(j).at(i) += (gZbar2oMZ4*sbar2*(-3 + 2*sbar2 + deltagZbar2oMZ4*(-3 + 2*sbar2) + delta_sbar*(-6 + 8*sbar2))*Md[i]*Md[j])/108.;
255 CddV1LR.at(i).at(j).at(j).at(i) += (gZbar2oMZ4*sbar2*(-3 + 2*sbar2 + deltagZbar2oMZ4*(-3 + 2*sbar2) + delta_sbar*(-6 + 8*sbar2))*Md[i]*Md[j])/108.;
256 CddV1LR.at(i).at(j).at(j).at(i) += -0.018518518518518517*((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[i]*Md[j]);
257 CddV8LR.at(i).at(j).at(j).at(i) += -0.027777777777777776*(gZbar2oMZ4*(-3 + 2*sbar2)*(-3 + (2 + 8*delta_sbar)*sbar2 + deltagZbar2oMZ4*(-3 + 2*sbar2))*Md[i]*Md[j]);
258 CddV8LR.at(i).at(j).at(j).at(i) += (gZbar2oMZ4*sbar2*(-3 + 2*sbar2 + deltagZbar2oMZ4*(-3 + 2*sbar2) + delta_sbar*(-6 + 8*sbar2))*Md[i]*Md[j])/18.;
259 CddV8LR.at(i).at(j).at(j).at(i) += (gZbar2oMZ4*sbar2*(-3 + 2*sbar2 + deltagZbar2oMZ4*(-3 + 2*sbar2) + delta_sbar*(-6 + 8*sbar2))*Md[i]*Md[j])/18.;
260 CddV8LR.at(i).at(j).at(j).at(i) += -0.1111111111111111*((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[i]*Md[j]);
261 CedSRL.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*(3 - 8*(1 + 2*delta_sbar)*sbar2 + 4*(1 + 4*delta_sbar)*pow(sbar2,2) + deltagZbar2oMZ4*(3 - 8*sbar2 + 4*pow(sbar2,2)))*Md[j]*Me[i])/12.;
262 CedSRL.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*sbar2*(1 + deltagZbar2oMZ4 + delta_sbar*(2 - 8*sbar2) - 2*sbar2 - 2*deltagZbar2oMZ4*sbar2)*Md[j]*Me[i])/6.;
263 CedSRL.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*(3 + delta_sbar*(6 - 8*sbar2) + deltagZbar2oMZ4*(3 - 2*sbar2) - 2*sbar2)*sbar2*Md[j]*Me[i])/6.;
264 CedSRL.at(i).at(i).at(j).at(j) += ((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[j]*Me[i])/3.;
265 CeeSRR.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*sbar2*(-1 + 2*sbar2 + deltagZbar2oMZ4*(-1 + 2*sbar2) + delta_sbar*(-2 + 8*sbar2))*Me[i]*Me[j])/4.;
266 CeeSRR.at(i).at(i).at(j).at(j) += -0.125*(gZbar2oMZ4*(-1 + 2*sbar2)*(-1 + (2 + 8*delta_sbar)*sbar2 + deltagZbar2oMZ4*(-1 + 2*sbar2))*Me[i]*Me[j]);
267 CeeSRR.at(i).at(i).at(j).at(j) += -0.5*((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Me[i]*Me[j]);
268 CeeSRR.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*sbar2*(-1 + 2*sbar2 + deltagZbar2oMZ4*(-1 + 2*sbar2) + delta_sbar*(-2 + 8*sbar2))*Me[i]*Me[j])/4.;
269 CedSRR.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*sbar2*(-1 + 2*sbar2 + deltagZbar2oMZ4*(-1 + 2*sbar2) + delta_sbar*(-2 + 8*sbar2))*Md[j]*Me[i])/6.;
270 CedSRR.at(i).at(i).at(j).at(j) += -0.08333333333333333*(gZbar2oMZ4*(3 - 8*(1 + 2*delta_sbar)*sbar2 + 4*(1 + 4*delta_sbar)*pow(sbar2,2) + deltagZbar2oMZ4*(3 - 8*sbar2 + 4*pow(sbar2,2)))*Md[j]*Me[i]);
271 CedSRR.at(i).at(i).at(j).at(j) += -0.3333333333333333*((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[j]*Me[i]);
272 CedSRR.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*sbar2*(-3 + 2*sbar2 + deltagZbar2oMZ4*(-3 + 2*sbar2) + delta_sbar*(-6 + 8*sbar2))*Md[j]*Me[i])/6.;
273 CddS1RR.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*sbar2*(-3 + 2*sbar2 + deltagZbar2oMZ4*(-3 + 2*sbar2) + delta_sbar*(-6 + 8*sbar2))*Md[i]*Md[j])/18.;
274 CddS1RR.at(i).at(i).at(j).at(j) += -0.027777777777777776*(gZbar2oMZ4*(-3 + 2*sbar2)*(-3 + (2 + 8*delta_sbar)*sbar2 + deltagZbar2oMZ4*(-3 + 2*sbar2))*Md[i]*Md[j]);
275 CddS1RR.at(i).at(i).at(j).at(j) += -0.1111111111111111*((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[i]*Md[j]);
276 CddS1RR.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*sbar2*(-3 + 2*sbar2 + deltagZbar2oMZ4*(-3 + 2*sbar2) + delta_sbar*(-6 + 8*sbar2))*Md[i]*Md[j])/18.;
278 for (
int k = 0; k < 3; k++)
280 CnunuVLL.at(i).at(i).at(j).at(k) += (gZ2oMZ2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) - (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))) / 16.;
281 CnunuVLL.at(j).at(k).at(i).at(i) += (gZ2oMZ2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) - (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))) / 16.;
282 CnunuVLL.at(j).at(i).at(i).at(k) += (gZ2oMZ2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) - (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))) / 16.;
283 CnunuVLL.at(i).at(j).at(k).at(i) += (gZ2oMZ2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", k, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", k, j)) - (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", k, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", k, j)))) / 16.;
284 CeeVLL.at(i).at(i).at(j).at(k) += (gZ2oMZ2 * (-1 + 2 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))) / 16.;
285 CeeVLL.at(j).at(k).at(i).at(i) += (gZ2oMZ2 * (-1 + 2 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))) / 16.;
286 CeeVLL.at(j).at(i).at(i).at(k) += (gZ2oMZ2 * (-1 + 2 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))) / 16.;
287 CeeVLL.at(i).at(j).at(k).at(i) += (gZ2oMZ2 * (-1 + 2 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", k, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", k, j)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", k, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", k, j)))) / 16.;
288 CnueVLL.at(i).at(i).at(j).at(k) += (gZ2oMZ2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))) / 4.;
289 CnueVLL.at(j).at(k).at(i).at(i) += (gZ2oMZ2 * (-1 + 2 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) - (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))) / 4.;
290 CnueVLL.at(j).at(i).at(i).at(k) += -0.5 * (g22oMW2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)));
291 CnueVLL.at(i).at(j).at(k).at(i) += -0.5 * (g22oMW2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)));
292 CnudVLL.at(j).at(k).at(i).at(i) += (gZ2oMZ2 * (-3 + 2 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) - (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))) / 12.;
293 CedVLL.at(j).at(k).at(i).at(i) += (gZ2oMZ2 * (-3 + 2 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))) / 12.;
294 CeeVRR.at(i).at(i).at(j).at(k) += (gZ2oMZ2 * sbar2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHeR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", j, k))) / 8.;
295 CeeVRR.at(j).at(k).at(i).at(i) += ((gZ2oMZ2 * sbar2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHeR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", j, k))) / 8.);
296 CeeVRR.at(j).at(i).at(i).at(k) += (gZ2oMZ2 * sbar2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHeR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", j, k))) / 8.;
297 CeeVRR.at(i).at(j).at(k).at(i) += (gZ2oMZ2 * sbar2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHeR", k, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", k, j))) / 8.;
298 CedVRR.at(j).at(k).at(i).at(i) += (gZ2oMZ2 * sbar2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHeR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", j, k))) / 6.;
299 CnueVLR.at(i).at(i).at(j).at(k) += (gZ2oMZ2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHeR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", j, k))) / 4.;
300 CnueVLR.at(j).at(k).at(i).at(i) += (gZ2oMZ2 * sbar2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) - (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))) / 2.;
301 CeeVLR.at(i).at(i).at(j).at(k) += (gZ2oMZ2 * (-1 + 2 * sbar2) *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHeR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", j, k))) / 4.;
302 CeeVLR.at(j).at(k).at(i).at(i) += (gZ2oMZ2 * sbar2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))) / 2.;
303 CnudVLR.at(j).at(k).at(i).at(i) += (gZ2oMZ2 * sbar2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) - (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))) / 6.;
304 CedVLR.at(j).at(k).at(i).at(i) += (gZ2oMZ2 * sbar2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))) / 6.;
305 CdeVLR.at(i).at(i).at(j).at(k) += (gZ2oMZ2 * (-3 + 2 * sbar2) *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHeR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", j, k))) / 12.;
307 CeeVLL.at(i).at(i).at(j).at(k) += (gZbar2oMZ3*(0.5 - sbar2)*
v2*((sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j)))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k)))*Me[k]))/(4.*sqrt(2));
308 CeeVLL.at(j).at(k).at(i).at(i) += -0.125*(gZbar2oMZ3*(-1 + 2*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j))*Me[j] + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k)))*Me[k]))/sqrt(2);
309 CeeVLL.at(j).at(i).at(i).at(k) += -0.125*(gZbar2oMZ3*(-1 + 2*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j))*Me[j] + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k)))*Me[k]))/sqrt(2);
310 CeeVLL.at(i).at(j).at(k).at(i) += -0.125*(gZbar2oMZ3*(-1 + 2*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j))*Me[j] + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k)))*Me[k]))/sqrt(2);
311 CnueVLL.at(i).at(i).at(j).at(k) += -0.5*(gZbar2oMZ3*
v2*((sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j)))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k)))*Me[k]))/sqrt(2);
312 CnueVLL.at(j).at(i).at(i).at(k) += (g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k))*Me[k])/sqrt(2);
313 CnueVLL.at(i).at(j).at(k).at(i) += -((g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k))*Me[k])/sqrt(2));
314 CedVLL.at(j).at(k).at(i).at(i) += -0.16666666666666666*(gZbar2oMZ3*(-3 + 2*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j))*Me[j] + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k)))*Me[k]))/sqrt(2);
315 CeeVRR.at(i).at(i).at(j).at(k) += -0.25*(gZbar2oMZ3*sbar2*
v2*((sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k)))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j)))*Me[k]))/sqrt(2);
316 CeeVRR.at(j).at(k).at(i).at(i) += -0.25*(gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k))*Me[j] + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j)))*Me[k]))/sqrt(2);
317 CeeVRR.at(j).at(i).at(i).at(k) += -0.25*(gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k))*Me[j] + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j)))*Me[k]))/sqrt(2);
318 CeeVRR.at(i).at(j).at(k).at(i) += -0.25*(gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j))*Me[j] + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k)))*Me[k]))/sqrt(2);
319 CedVRR.at(j).at(k).at(i).at(i) += -0.3333333333333333*(gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k))*Me[j] + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j)))*Me[k]))/sqrt(2);
320 CnueVLR.at(i).at(j).at(k).at(i) += -0.25*(g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k))*Me[i]*Me[k]);
321 CnueVLR.at(j).at(i).at(i).at(k) += -0.25*(g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k))*Me[i]*Me[k]);
322 CnueVLR.at(i).at(i).at(j).at(k) += -0.5*(gZbar2oMZ3*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k))*Me[j] + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j)))*Me[k]))/sqrt(2);
323 CeeVLR.at(i).at(j).at(k).at(i) += (gZbar2oMZ4*
v2*Me[i]*((
mySMEFT.getSMEFTCoeffEW(
"CHeR", k, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", k, j))*Me[j] - ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", k, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", k, j)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", k, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", k, j)))*Me[k]))/8.;
324 CeeVLR.at(j).at(i).at(i).at(k) += (gZbar2oMZ4*
v2*Me[i]*((
mySMEFT.getSMEFTCoeffEW(
"CHeR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", j, k))*Me[j] - ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))*Me[k]))/8.;
325 CeeVLR.at(i).at(i).at(j).at(k) += -0.5*(gZbar2oMZ3*(-1 + 2*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k))*Me[j] + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j)))*Me[k]))/sqrt(2);
326 CeeVLR.at(j).at(k).at(i).at(i) += -((gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j))*Me[j] + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k)))*Me[k]))/sqrt(2));
327 CedVLR.at(j).at(k).at(i).at(i) += -0.3333333333333333*(gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j))*Me[j] + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k)))*Me[k]))/sqrt(2);
328 CdeVLR.at(i).at(i).at(j).at(k) += -0.16666666666666666*(gZbar2oMZ3*(-3 + 2*sbar2)*
v2*((sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k)))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j)))*Me[k]))/sqrt(2);
329 CedSRL.at(j).at(k).at(i).at(i) += (gZbar2oMZ4*
v2*Md[i]*(-((
mySMEFT.getSMEFTCoeffEW(
"CHeR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", j, k))*Me[j]) + ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))*Me[k]))/4.;
330 CeeSRR.at(i).at(i).at(j).at(k) += (gZbar2oMZ4*
v2*Me[i]*((
mySMEFT.getSMEFTCoeffEW(
"CHeR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", j, k))*Me[j] - ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))*Me[k]))/8.;
331 CeeSRR.at(j).at(k).at(i).at(i) += (gZbar2oMZ4*
v2*Me[i]*((
mySMEFT.getSMEFTCoeffEW(
"CHeR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", j, k))*Me[j] - ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))*Me[k]))/8.;
332 CedSRR.at(j).at(k).at(i).at(i) += (gZbar2oMZ4*
v2*Md[i]*((
mySMEFT.getSMEFTCoeffEW(
"CHeR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", j, k))*Me[j] - ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))*Me[k]))/4.;
334 for (
int p = 0; p < 3; p++)
335 for (
int r = 0; r < 3; r++)
337 CnudVLL.at(i).at(i).at(j).at(k) +=
VdLd(j, p) * ((gZ2oMZ2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))) / 4.) *
VdL(r, k);
338 CedVLL.at(i).at(i).at(j).at(k) +=
VdLd(j, p) * ((gZ2oMZ2 * (-1 + 2 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))) / 4.) *
VdL(r, k);
339 CddVLL.at(i).at(i).at(j).at(k) +=
VdLd(j, p) * ((gZ2oMZ2 * (-3 + 2 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))) / 24.) *
VdL(r, k);
340 CddVLL.at(j).at(k).at(i).at(i) +=
VdLd(j, p) * ((gZ2oMZ2 * (-3 + 2 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))) / 24.) *
VdL(r, k);
341 CedVRR.at(i).at(i).at(j).at(k) +=
VdRd(j, p) * ((gZ2oMZ2 * sbar2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))) / 2.) *
VdR(r, k);
342 CddVRR.at(i).at(i).at(j).at(k) +=
VdRd(j, p) * ((gZ2oMZ2 * sbar2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))) / 12.) *
VdR(r, k);
343 CddVRR.at(j).at(k).at(i).at(i) +=
VdRd(j, p) * ((gZ2oMZ2 * sbar2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))) / 12.) *
VdR(r, k);
344 CnudVLR.at(i).at(i).at(j).at(k) +=
VdRd(j, p) * ((gZ2oMZ2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))) / 4.) *
VdR(r, k);
345 CnudVLR.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * (-0.5*(gZbar2oMZ3*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[j])/sqrt(2)) *
VdR(r,k);
346 CedVLR.at(i).at(i).at(j).at(k) +=
VdRd(j, p) * ((gZ2oMZ2 * (-1 + 2 * sbar2) *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))) / 4.) *
VdR(r, k);
347 CdeVLR.at(j).at(k).at(i).at(i) +=
VdLd(j, p) * ((gZ2oMZ2 * sbar2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))) / 2.) *
VdL(r, k);
348 CddV1LR.at(i).at(i).at(j).at(k) +=
VdRd(j, p) * ((gZ2oMZ2 * (-3 + 2 * sbar2) *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))) / 12.) *
VdR(r, k);
349 CddV1LR.at(j).at(k).at(i).at(i) +=
VdLd(j, p) * ((gZ2oMZ2 * sbar2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))) / 6.) *
VdL(r, k);
351 CnudVLL.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * (-0.5*(gZbar2oMZ3*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*Md[j])/sqrt(2)) *
VdL(r,k);
352 CnudVLL.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * (-0.5*(gZbar2oMZ3*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[k])/sqrt(2)) *
VdR(r,k);
353 CedVLL.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * (-((gZbar2oMZ3*(-0.5 + sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*Md[j])/sqrt(2))) *
VdL(r,k);
354 CedVLL.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * (-((gZbar2oMZ3*(-0.5 + sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[k])/sqrt(2))) *
VdR(r,k);
355 CddVLL.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * (-0.041666666666666664*(gZbar2oMZ3*(-3 + 2*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*Md[j])/sqrt(2)) *
VdL(r,k);
356 CddVLL.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * (-0.041666666666666664*(gZbar2oMZ3*(-3 + 2*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[k])/sqrt(2)) *
VdR(r,k);
357 CddVLL.at(j).at(k).at(i).at(i) +=
VdRd(j,p) * (-0.041666666666666664*(gZbar2oMZ3*(-3 + 2*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*Md[j])/sqrt(2)) *
VdL(r,k);
358 CddVLL.at(j).at(k).at(i).at(i) +=
VdLd(j,p) * (-0.041666666666666664*(gZbar2oMZ3*(-3 + 2*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[k])/sqrt(2)) *
VdR(r,k);
359 CedVRR.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * (-((gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*Md[k])/sqrt(2))) *
VdL(r,k);
360 CedVRR.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * (-((gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[j])/sqrt(2))) *
VdR(r,k);
361 CddVRR.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * (-0.16666666666666666*(gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*Md[k])/sqrt(2)) *
VdL(r,k);
362 CddVRR.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * (-0.16666666666666666*(gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[j])/sqrt(2)) *
VdR(r,k);
363 CddVRR.at(j).at(k).at(i).at(i) +=
VdRd(j,p) * (-0.16666666666666666*(gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*Md[k])/sqrt(2)) *
VdL(r,k);
364 CddVRR.at(j).at(k).at(i).at(i) +=
VdLd(j,p) * (-0.16666666666666666*(gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[j])/sqrt(2)) *
VdR(r,k);
365 CnudVLR.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * (-0.5*(gZbar2oMZ3*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*Md[k])/sqrt(2)) *
VdL(r,k);
366 CnudVLR.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * (-0.5*(gZbar2oMZ3*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[j])/sqrt(2)) *
VdR(r,k);
367 CedVLR.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * (-((gZbar2oMZ3*(-0.5 + sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*Md[k])/sqrt(2))) *
VdL(r,k);
368 CedVLR.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * (-((gZbar2oMZ3*(-0.5 + sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[j])/sqrt(2))) *
VdR(r,k);
369 CdeVLR.at(j).at(k).at(i).at(i) +=
VdRd(j,p) * (-((gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*Md[j])/sqrt(2))) *
VdL(r,k);
370 CdeVLR.at(j).at(k).at(i).at(i) +=
VdLd(j,p) * (-((gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[k])/sqrt(2))) *
VdR(r,k);
371 CddV1LR.at(i).at(j).at(k).at(i) +=
VdLd(k,p) * (-0.041666666666666664*(gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[i]*Md[k])) *
VdL(r,j);
372 CddV1LR.at(i).at(j).at(k).at(i) +=
VdRd(k,p) * ((gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*Md[i]*Md[j])/24.) *
VdR(r,j);
373 CddV1LR.at(j).at(i).at(i).at(k) +=
VdLd(j,p) * (-0.041666666666666664*(gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[i]*Md[k])) *
VdL(r,k);
374 CddV1LR.at(j).at(i).at(i).at(k) +=
VdRd(j,p) * ((gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*Md[i]*Md[j])/24.) *
VdR(r,k);
375 CddV1LR.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * (-0.16666666666666666*(gZbar2oMZ3*(-3 + 2*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*Md[k])/sqrt(2)) *
VdL(r,k);
376 CddV1LR.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * (-0.16666666666666666*(gZbar2oMZ3*(-3 + 2*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[j])/sqrt(2)) *
VdR(r,k);
377 CddV1LR.at(j).at(k).at(i).at(i) +=
VdRd(j,p) * (-0.3333333333333333*(gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*Md[j])/sqrt(2)) *
VdL(r,k);
378 CddV1LR.at(j).at(k).at(i).at(i) +=
VdLd(j,p) * (-0.3333333333333333*(gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[k])/sqrt(2)) *
VdR(r,k);
379 CddV8LR.at(i).at(j).at(k).at(i) +=
VdLd(k,p) * (-0.25*(gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[i]*Md[k])) *
VdL(r,j);
380 CddV8LR.at(i).at(j).at(k).at(i) +=
VdRd(k,p) * ((gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*Md[i]*Md[j])/4.) *
VdR(r,j);
381 CddV8LR.at(j).at(i).at(i).at(k) +=
VdLd(j,p) * (-0.25*(gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[i]*Md[k])) *
VdL(r,k);
382 CddV8LR.at(j).at(i).at(i).at(k) +=
VdRd(j,p) * ((gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*Md[i]*Md[j])/4.) *
VdR(r,k);
383 CedSRL.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * ((gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[j]*Me[i])/4.) *
VdL(r,k);
384 CedSRL.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * (-0.25*(gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*Md[k]*Me[i])) *
VdR(r,k);
385 CedSRR.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * ((gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*Md[j]*Me[i])/4.) *
VdR(r,k);
386 CedSRR.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * (-0.25*(gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[k]*Me[i])) *
VdL(r,k);
387 CddS1RR.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * ((gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*Md[i]*Md[j])/4.) *
VdR(r,k);
388 CddS1RR.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * (-0.25*(gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[i]*Md[k])) *
VdL(r,k);
389 CddS1RR.at(j).at(k).at(i).at(i) +=
VdLd(j,p) * (-0.25*(gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[i]*Md[k])) *
VdL(r,k);
390 CddS1RR.at(j).at(k).at(i).at(i) +=
VdRd(j,p) * ((gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*Md[i]*Md[j])/4.) *
VdR(r,k);
394 for (
int l = 0; l < 3; l++)
396 CnunuVLL.at(i).at(j).at(k).at(l) += (
mySMEFT.getSMEFTCoeffEW(
"CllR", i, j, k, l) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CllI", i, j, k, l));
397 CeeVLL.at(i).at(j).at(k).at(l) += (
mySMEFT.getSMEFTCoeffEW(
"CllR", i, j, k, l) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CllI", i, j, k, l));
398 CnueVLL.at(i).at(j).at(k).at(l) += (
mySMEFT.getSMEFTCoeffEW(
"CllR", i, j, k, l) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CllI", i, j, k, l)) + (
mySMEFT.getSMEFTCoeffEW(
"CllR", k, l, i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CllI", k, l, i, j));
399 CeeVRR.at(i).at(j).at(k).at(l) += (
mySMEFT.getSMEFTCoeffEW(
"CeeR", i, j, k, l) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeeI", i, j, k, l));
400 CnueVLR.at(i).at(j).at(k).at(l) += (
mySMEFT.getSMEFTCoeffEW(
"CleR", i, j, k, l) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CleI", i, j, k, l));
401 CeeVLR.at(i).at(j).at(k).at(l) += (
mySMEFT.getSMEFTCoeffEW(
"CleR", i, j, k, l) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CleI", i, j, k, l));
403 CeeVLR.at(i).at(j).at(k).at(l) += (oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", i, l) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", i, l))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", k, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", k, j)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", j, k) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", j, k)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YeR", k, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", k, j)))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", l, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", l, i))))/4.;
404 CedSRL.at(i).at(j).at(k).at(l) += VCKMd(k,l) * ((g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHl3R", i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", i, j))*Md[k]*Me[j])/2.);
405 CeeSRR.at(i).at(j).at(k).at(l) += -0.25*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", k, l) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", k, l))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", i, j)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i)))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", l, k) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", l, k))));
407 for (
int p = 0; p < 3; p++)
408 for (
int r = 0; r < 3; r++)
410 CnudVLL.at(i).at(j).at(k).at(l) +=
VdLd(k, p) * ((
mySMEFT.getSMEFTCoeffEW(
"Clq1R", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Clq1I", i, j, p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"Clq3R", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Clq3I", i, j, p, r))) *
VdL(r, l);
411 CedVLL.at(i).at(j).at(k).at(l) +=
VdLd(k, p) * ((
mySMEFT.getSMEFTCoeffEW(
"Clq1R", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Clq1I", i, j, p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"Clq3R", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Clq3I", i, j, p, r))) *
VdL(r, l);
412 CedVRR.at(i).at(j).at(k).at(l) +=
VdRd(k, p) * ((
mySMEFT.getSMEFTCoeffEW(
"CedR", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CedI", i, j, p, r))) *
VdR(r, l);
413 CnudVLR.at(i).at(j).at(k).at(l) +=
VdRd(k, p) * ((
mySMEFT.getSMEFTCoeffEW(
"CldR", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CldI", i, j, p, r))) *
VdR(r, l);
414 CedVLR.at(i).at(j).at(k).at(l) +=
VdRd(k, p) * ((
mySMEFT.getSMEFTCoeffEW(
"CldR", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CldI", i, j, p, r))) *
VdR(r, l);
415 CdeVLR.at(i).at(j).at(k).at(l) +=
VdLd(i, p) * ((
mySMEFT.getSMEFTCoeffEW(
"CqeR", p, r, k, l) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CqeI", p, r, k, l))) *
VdL(r, j);
416 CedSRL.at(i).at(j).at(k).at(l) +=
VdRd(k, p) * ((
mySMEFT.getSMEFTCoeffEW(
"CledqR", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CledqI", i, j, p, r))) *
VdL(r, l);
418 CedSRL.at(i).at(j).at(k).at(l) +=
VdRd(k,p) * (-0.5*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", i, j))*(
mySMEFT.getSMEFTCoeffEW(
"YdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI", p, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", r, p)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI", p, r)))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i)))))*
VdL(r,l);
419 CedSRR.at(i).at(j).at(k).at(l) +=
VdLd(k,p) * (-0.5*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", i, j))*(
mySMEFT.getSMEFTCoeffEW(
"YdR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI", r, p)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", p, r)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YdR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI", r, p)))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i))))) *
VdR(r,l);
422 for (
int s = 0;
s < 3;
s++)
423 for (
int t = 0;
t < 3;
t++)
425 CddVLL.at(i).at(j).at(k).at(l) +=
VdLd(i, p) *
VdLd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"Cqq1R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqq1I", p, r,
s,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"Cqq3R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqq3I", p, r,
s,
t))) *
VdL(r, j) *
VdL(
t, l);
426 CddVRR.at(i).at(j).at(k).at(l) +=
VdRd(i, p) *
VdRd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"CddR", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CddI", p, r,
s,
t))) *
VdR(r, j) *
VdR(
t, l);
427 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i, p) *
VdRd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"Cqd1R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqd1I", p, r,
s,
t))) *
VdL(r, j) *
VdR(
t, l);
428 CddV8LR.at(i).at(j).at(k).at(l) +=
VdLd(i, p) *
VdRd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"Cqd8R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqd8I", p, r,
s,
t))) *
VdL(r, j) *
VdR(
t, l);
429 CddVLL.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * (-0.125*(gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t))*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))) - (gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))/8. - (gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t))))/8.) *
VdL(r,j) *
VdL(
t,l);
430 CddVRR.at(i).at(j).at(k).at(l) +=
VdRd(i,p) *
VdRd(k,
s) * ((gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s,
t)))/8.) *
VdR(r,j) *
VdR(
t,l);
431 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * (-0.25*(gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s,
t))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))))) *
VdL(r,j) *
VdR(
t,l);
433 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", r,
s)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
t, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
t, p))))/12.) *
VdL(r,j) *
VdR(
t,l);
434 CddV8LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", r,
s)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
t, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
t, p))))/2.) *
VdL(r,j) *
VdR(
t,l);
435 CddS1RR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * (-0.25*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI",
s,
t))*(
mySMEFT.getSMEFTCoeffEW(
"YdR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI", r, p)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", p, r)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YdR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI", r, p)))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
t,
s))))) *
VdR(r,j) *
VdR(
t,l);
436 CddVLL.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI",
t,
s)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI",
t,
s)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[k])/(8.*sqrt(2))) *
VdL(r,j) *
VdL(
t,l);
437 CddVLL.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI",
s,
t)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI",
s,
t)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[l])/(8.*sqrt(2))) *
VdL(r,j) *
VdR(
t,l);
438 CddVLL.at(i).at(j).at(k).at(l) +=
VdRd(i,p) *
VdLd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))*Md[i])/(8.*sqrt(2))) *
VdL(r,j) *
VdL(
t,l);
439 CddVLL.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))*Md[j])/(8.*sqrt(2))) *
VdR(r,j) *
VdL(
t,l);
440 CddVRR.at(i).at(j).at(k).at(l) +=
VdRd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI",
t,
s)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI",
t,
s)))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*Md[l])/(4.*sqrt(2)))*
VdR(r,j) *
VdL(
t,l);
441 CddVRR.at(i).at(j).at(k).at(l) +=
VdRd(i,p) *
VdLd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI",
s,
t)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI",
s,
t)))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*Md[k])/(4.*sqrt(2))) *
VdR(r,j) *
VdR(
t,l);
442 CddVRR.at(i).at(j).at(k).at(l) +=
VdRd(i,p) *
VdRd(k,
s) *((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s,
t))*Md[j])/(4.*sqrt(2))) *
VdL(r,j) *
VdR(
t,l);
443 CddVRR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) *((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p ,r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p ,r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p ,r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p ,r)))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s,
t))*Md[i])/(4.*sqrt(2))) *
VdR(r,j) *
VdR(
t,l);
444 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * (-0.041666666666666664*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", p,
t))*(3*
mySMEFT.FlagQuadraticWC*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", r,
s)))))) *
VdL(r,j) *
VdR(
t,l);
445 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI",
t,
s)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI",
t,
s)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[l])/(2.*sqrt(2))) *
VdL(r,j) *
VdR(
t,l);
446 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI",
s,
t)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI",
s,
t)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[k])/(2.*sqrt(2))) *
VdL(r,j) *
VdR(
t,l);
447 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s,
t))*Md[i])/(2.*sqrt(2))) *
VdL(r,j) *
VdR(
t,l);
448 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s,
t))*Md[j])/(2.*sqrt(2))) *
VdL(r,j) *
VdR(
t,l);
449 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * (-0.041666666666666664*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s, r)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p,
t)))*Md[k]*Md[l])) *
VdL(r,j) *
VdR(
t,l);
450 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s, r))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p,
t)))*Md[j]*Md[l])/24.) *
VdL(r,j) *
VdR(
t,l);
451 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p,
t))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s, r)))*Md[i]*Md[k])/24.) *
VdL(r,j) *
VdR(
t,l);
452 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * (-0.041666666666666664*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s, r))*Md[i]*Md[j])) *
VdL(r,j) *
VdR(
t,l);
453 CddV8LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * (-0.25*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", p,
t))*(3*
mySMEFT.FlagQuadraticWC*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", r,
s)))))) *
VdL(r,j) *
VdR(
t,l);
454 CddV8LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * (-0.25*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s, r)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p,
t)))*Md[k]*Md[l])) *
VdL(r,j) *
VdR(
t,l);
455 CddV8LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s, r))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p,
t)))*Md[j]*Md[l])/4.) *
VdL(r,j) *
VdR(
t,l);
456 CddV8LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p,
t))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s, r)))*Md[i]*Md[k])/4.) *
VdL(r,j) *
VdR(
t,l);
457 CddV8LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * (-0.25*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s, r))*Md[i]*Md[j])) *
VdL(r,j) *
VdR(
t,l);
458 CddS1RR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", p, r))*(3*
mySMEFT.FlagQuadraticWC*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI",
s,
t)))))/8.) *
VdR(r,j) *
VdR(
t,l);
459 CddS1RR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s,
t))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[j]*Md[k])/4.) *
VdR(r,j) *
VdR(
t,l);
460 CddS1RR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * (-0.25*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))*Md[j]*Md[l])) *
VdR(r,j) *
VdR(
t,l);
461 CddS1RR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * (-0.25*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p ,r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p ,r))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s,
t))*Md[i]*Md[k])) *
VdR(r,j) *
VdR(
t,l);
462 CddS1RR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p ,r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p ,r))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))*Md[i]*Md[l])/4.) *
VdR(r,j) *
VdR(
t,l);
467 for (
int l = 0; l < 2; l++)
469 CnueduVLL.at(i).at(j).at(k).at(l) += VCKMd(k, l) * (-0.5 * (g22oMW2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", i, j))));
471 CnueduVLL.at(i).at(j).at(k).at(l) += VCKMd(k,l) *((g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", i, j))*Me[j])/sqrt(2));
472 CnueduSRR.at(i).at(j).at(k).at(l) += VCKMd(k,l) * (-0.5*(g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHl3R", i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", i, j))*Me[j]*Mu[l]));
474 for (
int p = 0; p < 3; p++)
475 for (
int r = 0; r < 3; r++)
477 CnueduVLL.at(i).at(j).at(k).at(l) +=
VdLd(k, p) * (2 * (
mySMEFT.getSMEFTCoeffEW(
"Clq3R", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Clq3I", i, j, p, r))) *
VuL(r, l);
478 CnueduSRL.at(i).at(j).at(k).at(l) +=
VdRd(k, p) * ((
mySMEFT.getSMEFTCoeffEW(
"CledqR", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CledqI", i, j, p, r))) *
VuL(r, l);
479 CnueduSRR.at(i).at(j).at(k).at(l) +=
VdLd(k, p) * ((
mySMEFT.getSMEFTCoeffEW(
"Clequ1R", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Clequ1I", i, j, p, r))) *
VuR(r, l);
480 CnueduTRR.at(i).at(j).at(k).at(l) +=
VdLd(k, p) * ((
mySMEFT.getSMEFTCoeffEW(
"Clequ3R", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Clequ3I", i, j, p, r))) *
VuR(r, l);
484 for (
int k = 0; k < 2; k++)
486 CnueduVLL.at(i).at(i).at(j).at(k) += VCKMd(j, k) * (-0.5 * (delta_g22oMW2 * g22oMW2));
488 CnueduSRR.at(i).at(i).at(j).at(k) += VCKMd(j,k) * (-0.5*((1 + deltag2bar2oMW4)*g2bar2oMW4*Me[i]*Mu[k]));
490 for (
int p = 0; p < 3; p++)
491 for (
int r = 0; r < 3; r++)
493 CnueduVLL.at(i).at(i).at(j).at(k) +=
VdLd(j, p) * (-0.5 * (g22oMW2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", r, p)))) *
VuL(r, k);
494 CnueduVLR.at(i).at(i).at(j).at(k) +=
VdRd(j, p) * (-0.25 * (g22oMW2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHudR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHudI", r, p)))) *
VuR(r, k);
496 CnueduVLL.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * (-((g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p))*Md[j])/sqrt(2))) *
VuL(r,k);
497 CnueduVLL.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * (-((g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r))*Mu[k])/sqrt(2))) *
VuR(r,k);
498 CnueduVLR.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * (-((g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p))*Mu[k])/sqrt(2))) *
VuL(r,k);
499 CnueduVLR.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * (-((g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r))*Md[j])/sqrt(2))) *
VuR(r,k);
500 CnueduSRL.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * ((g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", r, p))*Md[j]*Me[i])/2.) *
VuL(r,k);
501 CnueduSRL.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * (-0.25*(g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHudR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHudI", r, p))*Me[i]*Mu[k])) *
VuR(r,k);
502 CnueduSRR.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * ((g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHudR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHudI", r, p))*Md[j]*Me[i])/4.) *
VuR(r,k);
503 CnueduSRR.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * (-0.5*(g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", r, p))*Me[i]*Mu[k])) *
VuL(r,k);
506 for (
int l = 0; l < 2; l++)
509 CduV1LR.at(i).at(j).at(k).at(l) += VCKMd(i,l) *
VCKM(k,j) * (-0.08333333333333333*((1 + deltag2bar2oMW4)*g2bar2oMW4*Mu[k]*Mu[l]));
510 CduV8LR.at(i).at(j).at(k).at(l) += VCKMd(i,l) *
VCKM(k,j) * (-0.5*((1 + deltag2bar2oMW4)*g2bar2oMW4*Mu[k]*Mu[l]));
512 for (
int p = 0; p < 3; p++)
513 for (
int r = 0; r < 3; r++)
515 CnuuVLL.at(i).at(j).at(k).at(l) +=
VuLd(k, p) * ((
mySMEFT.getSMEFTCoeffEW(
"Clq1R", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Clq1I", i, j, p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"Clq3R", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Clq3I", i, j, p, r))) *
VuL(r, l);
516 CeuVLL.at(i).at(j).at(k).at(l) +=
VuLd(k, p) * ((
mySMEFT.getSMEFTCoeffEW(
"Clq1R", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Clq1I", i, j, p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"Clq3R", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Clq3I", i, j, p, r))) *
VuL(r, l);
517 CeuVRR.at(i).at(j).at(k).at(l) +=
VuRd(k, p) * ((
mySMEFT.getSMEFTCoeffEW(
"CeuR", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeuI", i, j, p, r))) *
VuR(r, l);
518 CnuuVLR.at(i).at(j).at(k).at(l) +=
VuRd(k, p) * ((
mySMEFT.getSMEFTCoeffEW(
"CluR", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CluI", i, j, p, r))) *
VuR(r, l);
519 CeuVLR.at(i).at(j).at(k).at(l) +=
VuRd(k, p) * ((
mySMEFT.getSMEFTCoeffEW(
"CluR", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CluI", i, j, p, r))) *
VuR(r, l);
520 CeuSRR.at(i).at(j).at(k).at(l) +=
VuLd(k, p) * (-(
mySMEFT.getSMEFTCoeffEW(
"Clequ1R", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Clequ1I", i, j, p, r))) *
VuR(r, l);
521 CeuTRR.at(i).at(j).at(k).at(l) +=
VuLd(k, p) * (-(
mySMEFT.getSMEFTCoeffEW(
"Clequ3R", i, j, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Clequ3I", i, j, p, r))) *
VuR(r, l);
523 CduV1LR.at(i).at(j).at(k).at(l) += VCKMd(i,l) *
VuLd(k,p) * (-0.08333333333333333*(g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))*Mu[k]*Mu[l])) *
VdL(r,j);
524 CduV1LR.at(i).at(j).at(k).at(l) += VCKMd(i,l) *
VuRd(k,p) * ((g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHudR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHudI", p, r))*Md[j]*Mu[l])/24.) *
VdR(r,j);
525 CduV1LR.at(i).at(j).at(k).at(l) +=
VCKM(k,j) *
VdLd(i,p) * (-0.08333333333333333*(g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", r, p))*Mu[k]*Mu[l])) *
VuL(r,l);
526 CduV1LR.at(i).at(j).at(k).at(l) +=
VCKM(k,j) *
VdRd(i,p) * ((g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHudR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHudI", r, p))*Md[i]*Mu[k])/24.) *
VuR(r,l);
527 CduV8LR.at(i).at(j).at(k).at(l) += VCKMd(i,l) *
VuLd(k,p) * (-0.5*(g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))*Mu[k]*Mu[l])) *
VdL(r,j);
528 CduV8LR.at(i).at(j).at(k).at(l) += VCKMd(i,l) *
VuRd(k,p) * ((g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHudR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHudI", p, r))*Md[j]*Mu[l])/4.) *
VdR(r,j);
529 CduV8LR.at(i).at(j).at(k).at(l) +=
VCKM(k,j) *
VdLd(i,p) * (-0.5*(g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", r, p))*Mu[k]*Mu[l])) *
VuL(r,l);
530 CduV8LR.at(i).at(j).at(k).at(l) +=
VCKM(k,j) *
VdRd(i,p) * ((g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHudR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHudI", r, p))*Md[i]*Mu[k])/4.) *
VuR(r,l);
531 CeuSRL.at(i).at(j).at(k).at(l) +=
VuRd(k,p) * (-0.5*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", r, p))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", i, j)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i)))*(
mySMEFT.getSMEFTCoeffEW(
"YuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI", p, r)))))*
VuL(r,l);
532 CeuSRR.at(i).at(j).at(k).at(l) +=
VuLd(k,p) * (-0.5*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p, r))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", i, j)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i)))*(
mySMEFT.getSMEFTCoeffEW(
"YuR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI", r, p))))) *
VuR(r,l);
534 for (
int s = 0;
s < 3;
s++)
535 for (
int t = 0;
t < 3;
t++)
537 CduV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i, p) *
VuRd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"Cqu1R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqu1I", p, r,
s,
t))) *
VdL(r, j) *
VuR(
t, l);
538 CduV8LR.at(i).at(j).at(k).at(l) +=
VdLd(i, p) *
VuRd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"Cqu8R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqu8I", p, r,
s,
t))) *
VdL(r, j) *
VuR(
t, l);
543 for (
int p = 0; p < 3; p++)
544 for (
int r = 0; r < 3; r++)
546 Cdg.at(i).at(j) += vTosq2 *
VdLd(i, p) * (-(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)) * sbar + (
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) * cbar) *
VdR(r, j);
547 CdG.at(i).at(j) += vTosq2 *
VdLd(i, p) * (
mySMEFT.getSMEFTCoeffEW(
"CdGR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdGI", p, r)) *
VdR(r, j);
551 for (
int i = 0; i < 3; i++)
552 for (
int j = 0; j < 2; j++)
554 CnuuVLL.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * (8 * delta_sbar * sbar2 + delta_gZ2oMZ2 * (-3 + 4 * sbar2))) / 12.;
555 CeuVLL.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * (4 * delta_sbar * sbar2 * (-5 + 8 * sbar2) + delta_gZ2oMZ2 * (3 - 10 * sbar2 + 8 * sbar2 * sbar2))) / 12.;
556 CeuVRR.at(i).at(i).at(j).at(j) += (2 * (delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2) / 3.;
557 CnuuVLR.at(i).at(i).at(j).at(j) += ((delta_gZ2oMZ2 + 2 * delta_sbar) * gZ2oMZ2 * sbar2) / 3.;
558 CeuVLR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * sbar2 * (-delta_gZ2oMZ2 - 2 * delta_sbar + 2 * delta_gZ2oMZ2 * sbar2 + 8 * delta_sbar * sbar2)) / 3.;
559 CduV1LR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * sbar2 * (-3 * delta_gZ2oMZ2 - 6 * delta_sbar + 2 * delta_gZ2oMZ2 * sbar2 + 8 * delta_sbar * sbar2)) / 9.;
561 CeuSRL.at(i).at(i).at(j).at(j) += -0.08333333333333333*(gZbar2oMZ4*(3 - 10*(1 + 2*delta_sbar)*sbar2 + 8*(1 + 4*delta_sbar)*pow(sbar2,2) + deltagZbar2oMZ4*(3 - 10*sbar2 + 8*pow(sbar2,2)))*Me[i]*Mu[j]);
562 CeuSRL.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*sbar2*(-1 + 2*sbar2 + deltagZbar2oMZ4*(-1 + 2*sbar2) + delta_sbar*(-2 + 8*sbar2))*Me[i]*Mu[j])/3.;
563 CeuSRL.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*sbar2*(-3 + 4*sbar2 + deltagZbar2oMZ4*(-3 + 4*sbar2) + 2*delta_sbar*(-3 + 8*sbar2))*Me[i]*Mu[j])/6.;
564 CeuSRL.at(i).at(i).at(j).at(j) += (-2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Me[i]*Mu[j])/3.;
565 CeuSRR.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*sbar2*(1 + deltagZbar2oMZ4 + delta_sbar*(2 - 8*sbar2) - 2*sbar2 - 2*deltagZbar2oMZ4*sbar2)*Me[i]*Mu[j])/3.;
566 CeuSRR.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*(3 - 10*(1 + 2*delta_sbar)*sbar2 + 8*(1 + 4*delta_sbar)*pow(sbar2,2) + deltagZbar2oMZ4*(3 - 10*sbar2 + 8*pow(sbar2,2)))*Me[i]*Mu[j])/12.;
567 CeuSRR.at(i).at(i).at(j).at(j) += (2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Me[i]*Mu[j])/3.;
568 CeuSRR.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*(3 + delta_sbar*(6 - 16*sbar2) + deltagZbar2oMZ4*(3 - 4*sbar2) - 4*sbar2)*sbar2*Me[i]*Mu[j])/6.;
570 for (
int k = 0; k < 2; k++)
571 for (
int p = 0; p < 3; p++)
572 for (
int r = 0; r < 3; r++)
574 CudV1LL.at(j).at(k).at(i).at(i) +=
VuLd(j, p) * ((gZ2oMZ2 * (-3 + 2 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))) / 12.) *
VuL(r, k);
575 CnuuVLL.at(i).at(i).at(j).at(k) +=
VuLd(j, p) * ((gZ2oMZ2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))) / 4.) *
VuL(r, k);
576 CeuVLL.at(i).at(i).at(j).at(k) +=
VuLd(j, p) * ((gZ2oMZ2 * (-1 + 2 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))) / 4.) *
VuL(r, k);
577 CudV1RR.at(j).at(k).at(i).at(i) +=
VuRd(j, p) * ((gZ2oMZ2 * sbar2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))) / 6.) *
VuR(r, k);
578 CeuVRR.at(i).at(i).at(j).at(k) +=
VuRd(j, p) * ((gZ2oMZ2 * sbar2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))) / 2.) *
VuR(r, k);
579 CnuuVLR.at(i).at(i).at(j).at(k) +=
VuRd(j, p) * ((gZ2oMZ2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))) / 4.) *
VuR(r, k);
580 CeuVLR.at(i).at(i).at(j).at(k) +=
VuRd(j, p) * ((gZ2oMZ2 * (-1 + 2 * sbar2) *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))) / 4.) *
VuR(r, k);
581 CduV1LR.at(i).at(i).at(j).at(k) +=
VuRd(j, p) * ((gZ2oMZ2 * (-3 + 2 * sbar2) *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))) / 12.) *
VuR(r, k);
582 CudV1LR.at(j).at(k).at(i).at(i) +=
VuLd(j, p) * ((gZ2oMZ2 * sbar2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))) / 6.) *
VuL(r, k);
583 CueVLR.at(j).at(k).at(i).at(i) +=
VuLd(j, p) * ((gZ2oMZ2 * sbar2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))) / 2.) *
VuL(r, k);
585 CnuuVLL.at(i).at(i).at(j).at(k) +=
VuRd(j,p) * (-0.5*(gZbar2oMZ3*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[j])/sqrt(2)) *
VuL(r,k);
586 CnuuVLL.at(i).at(i).at(j).at(k) +=
VuLd(j,p) * (-0.5*(gZbar2oMZ3*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[k])/sqrt(2)) *
VuR(r,k);
587 CeuVLL.at(i).at(i).at(j).at(k) +=
VuRd(j,p) * (-((gZbar2oMZ3*(-0.5 + sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[j])/sqrt(2))) *
VuL(r,k);
588 CeuVLL.at(i).at(i).at(j).at(k) +=
VuLd(j,p) * (-((gZbar2oMZ3*(-0.5 + sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[k])/sqrt(2))) *
VuR(r,k);
589 CudV1LL.at(j).at(k).at(i).at(i) +=
VuRd(j,p) * (-0.16666666666666666*(gZbar2oMZ3*(-3 + 2*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[j])/sqrt(2)) *
VuL(r,k);
590 CudV1LL.at(j).at(k).at(i).at(i) +=
VuLd(j,p) * (-0.16666666666666666*(gZbar2oMZ3*(-3 + 2*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[k])/sqrt(2)) *
VuR(r,k);
591 CeuVRR.at(i).at(i).at(j).at(k) +=
VuRd(j,p) * (-((gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[k])/sqrt(2))) *
VuL(r,k);
592 CeuVRR.at(i).at(i).at(j).at(k) +=
VuLd(j,p) * (-((gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[j])/sqrt(2))) *
VuR(r,k);
593 CudV1RR.at(j).at(k).at(i).at(i) +=
VuRd(j,p) * (-0.3333333333333333*(gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[k])/sqrt(2)) *
VuL(r,k);
594 CudV1RR.at(j).at(k).at(i).at(i) +=
VuLd(j,p) * (-0.3333333333333333*(gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[j])/sqrt(2)) *
VuR(r,k);
595 CnuuVLR.at(i).at(i).at(j).at(k) +=
VuRd(j,p) * (-0.5*(gZbar2oMZ3*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[k])/sqrt(2)) *
VuL(r,k);
596 CnuuVLR.at(i).at(i).at(j).at(k) +=
VuLd(j,p) * (-0.5*(gZbar2oMZ3*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[j])/sqrt(2)) *
VuR(r,k);
597 CeuVLR.at(i).at(i).at(j).at(k) +=
VuRd(j,p) * (-((gZbar2oMZ3*(-0.5 + sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[k])/sqrt(2))) *
VuL(r,k);
598 CeuVLR.at(i).at(i).at(j).at(k) +=
VuLd(j,p) * (-((gZbar2oMZ3*(-0.5 + sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[j])/sqrt(2))) *
VuR(r,k);
599 CueVLR.at(j).at(k).at(i).at(i) +=
VuRd(j,p) * (-((gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[j])/sqrt(2))) *
VuL(r,k);
600 CueVLR.at(j).at(k).at(i).at(i) +=
VuLd(j,p) * (-((gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[k])/sqrt(2))) *
VuR(r,k);
601 CudV1LR.at(j).at(k).at(i).at(i) +=
VuRd(j,p) * (-0.3333333333333333*(gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[j])/sqrt(2)) *
VuL(r,k);
602 CudV1LR.at(j).at(k).at(i).at(i) +=
VuLd(j,p) * (-0.3333333333333333*(gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[k])/sqrt(2)) *
VuR(r,k);
603 CduV1LR.at(i).at(i).at(j).at(k) +=
VuRd(j,p) * (-0.16666666666666666*(gZbar2oMZ3*(-3 + 2*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[k])/sqrt(2)) *
VuL(r,k);
604 CduV1LR.at(i).at(i).at(j).at(k) +=
VuLd(j,p) * (-0.16666666666666666*(gZbar2oMZ3*(-3 + 2*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[j])/sqrt(2)) *
VuR(r,k);
605 CudduV1LR.at(j).at(i).at(i).at(k) +=
VuLd(j,p) * (-0.041666666666666664*(gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[i]*Mu[k])) *
VuL(r,k);
606 CudduV1LR.at(j).at(i).at(i).at(k) +=
VuRd(j,p) * ((gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*Md[i]*Mu[j])/24.) *
VuR(r,k);
607 CudduV8LR.at(j).at(i).at(i).at(k) +=
VuLd(j,p) * (-0.25*(gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[i]*Mu[k])) *
VuL(r,k);
608 CudduV8LR.at(j).at(i).at(i).at(k) +=
VuRd(j,p) * ((gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*Md[i]*Mu[j])/4.) *
VuR(r,k);
609 CeuSRL.at(i).at(i).at(j).at(k) +=
VuLd(j,p) * ((gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Me[i]*Mu[j])/4.) *
VuL(r,k);
610 CeuSRL.at(i).at(i).at(j).at(k) +=
VuRd(j,p) * (-0.25*(gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*Me[i]*Mu[k])) *
VuR(r,k);
611 CeuSRR.at(i).at(i).at(j).at(k) +=
VuRd(j,p) * ((gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*Me[i]*Mu[j])/4.) *
VuR(r,k);
612 CeuSRR.at(i).at(i).at(j).at(k) +=
VuLd(j,p) * (-0.25*(gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Me[i]*Mu[k])) *
VuL(r,k);
613 CudS1RR.at(j).at(k).at(i).at(i) +=
VuLd(j,p) * (-0.25*(gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[i]*Mu[k])) *
VuL(r,k);
614 CudS1RR.at(j).at(k).at(i).at(i) +=
VuRd(j,p) * ((gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*Md[i]*Mu[j])/4.) *
VuR(r,k);
619 for (
int i = 0; i < 2; i++)
620 for (
int j = 0; j < 3; j++)
622 CudV1LL.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * (4 * delta_sbar * sbar2 * (-9 + 8 * sbar2) + delta_gZ2oMZ2 * (9 - 18 * sbar2 + 8 * sbar2 * sbar2))) / 36.;
623 CudV1RR.at(i).at(i).at(j).at(j) += (2 * (delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2) / 9.;
624 CudV1LR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * sbar2 * (delta_gZ2oMZ2 * (-3 + 4 * sbar2) + 2 * delta_sbar * (-3 + 8 * sbar2))) / 18.;
625 CueVLR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * sbar2 * (delta_gZ2oMZ2 * (-3 + 4 * sbar2) + 2 * delta_sbar * (-3 + 8 * sbar2))) / 6.;
627 CudduV1LR.at(i).at(j).at(j).at(i) += (gZbar2oMZ4*(9 - 18*(1 + 2*delta_sbar)*sbar2 + 8*(1 + 4*delta_sbar)*pow(sbar2,2) + deltagZbar2oMZ4*(9 - 18*sbar2 + 8*pow(sbar2,2)))*Md[j]*Mu[i])/216.;
628 CudduV1LR.at(i).at(j).at(j).at(i) += -0.009259259259259259*(gZbar2oMZ4*sbar2*(-3 + 4*sbar2 + deltagZbar2oMZ4*(-3 + 4*sbar2) + 2*delta_sbar*(-3 + 8*sbar2))*Md[j]*Mu[i]);
629 CudduV1LR.at(i).at(j).at(j).at(i) += -0.018518518518518517*(gZbar2oMZ4*sbar2*(-3 + 2*sbar2 + deltagZbar2oMZ4*(-3 + 2*sbar2) + delta_sbar*(-6 + 8*sbar2))*Md[j]*Mu[i]);
630 CudduV1LR.at(i).at(j).at(j).at(i) += ((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[j]*Mu[i])/27.;
631 CudduV8LR.at(i).at(j).at(j).at(i) += (gZbar2oMZ4*(9 - 18*(1 + 2*delta_sbar)*sbar2 + 8*(1 + 4*delta_sbar)*pow(sbar2,2) + deltagZbar2oMZ4*(9 - 18*sbar2 + 8*pow(sbar2,2)))*Md[j]*Mu[i])/36.;
632 CudduV8LR.at(i).at(j).at(j).at(i) += -0.05555555555555555*(gZbar2oMZ4*sbar2*(-3 + 4*sbar2 + deltagZbar2oMZ4*(-3 + 4*sbar2) + 2*delta_sbar*(-3 + 8*sbar2))*Md[j]*Mu[i]);
633 CudduV8LR.at(i).at(j).at(j).at(i) += -0.1111111111111111*(gZbar2oMZ4*sbar2*(-3 + 2*sbar2 + deltagZbar2oMZ4*(-3 + 2*sbar2) + delta_sbar*(-6 + 8*sbar2))*Md[j]*Mu[i]);
634 CudduV8LR.at(i).at(j).at(j).at(i) += (2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[j]*Mu[i])/9.;
635 CudS1RR.at(i).at(i).at(j).at(j) += -0.05555555555555555*(gZbar2oMZ4*sbar2*(-3 + 4*sbar2 + deltagZbar2oMZ4*(-3 + 4*sbar2) + 2*delta_sbar*(-3 + 8*sbar2))*Md[j]*Mu[i]);
636 CudS1RR.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*(9 - 18*(1 + 2*delta_sbar)*sbar2 + 8*(1 + 4*delta_sbar)*pow(sbar2,2) + deltagZbar2oMZ4*(9 - 18*sbar2 + 8*pow(sbar2,2)))*Md[j]*Mu[i])/36.;
637 CudS1RR.at(i).at(i).at(j).at(j) += (2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[j]*Mu[i])/9.;
638 CudS1RR.at(i).at(i).at(j).at(j) += -0.1111111111111111*(gZbar2oMZ4*sbar2*(-3 + 2*sbar2 + deltagZbar2oMZ4*(-3 + 2*sbar2) + delta_sbar*(-6 + 8*sbar2))*Md[j]*Mu[i]);
640 for (
int k = 0; k < 3; k++)
642 CeuVLL.at(j).at(k).at(i).at(i) += -0.08333333333333333 * (gZ2oMZ2 * (-3 + 4 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k))));
643 CnuuVLL.at(j).at(k).at(i).at(i) += -0.08333333333333333 * (gZ2oMZ2 * (-3 + 4 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) - (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k))));
644 CeuVRR.at(j).at(k).at(i).at(i) += -0.3333333333333333 * (gZ2oMZ2 * sbar2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHeR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", j, k)));
645 CnuuVLR.at(j).at(k).at(i).at(i) += (gZ2oMZ2 * sbar2 *
v2 * (-(
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))) / 3.;
646 CeuVLR.at(j).at(k).at(i).at(i) += -0.3333333333333333 * (gZ2oMZ2 * sbar2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k))));
647 CueVLR.at(i).at(i).at(j).at(k) += (gZ2oMZ2 * (3 - 4 * sbar2) *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHeR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", j, k))) / 12.;
649 CeuVLL.at(j).at(k).at(i).at(i) += (gZbar2oMZ3*(-3 + 4*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j))*Me[j] + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k)))*Me[k]))/(6.*sqrt(2));
650 CeuVRR.at(j).at(k).at(i).at(i) += (sqrt(2)*gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k))*Me[j] + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j)))*Me[k]))/3.;
651 CeuVLR.at(j).at(k).at(i).at(i) += (sqrt(2)*gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j))*Me[j] + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k)))*Me[k]))/3.;
652 CueVLR.at(i).at(i).at(j).at(k) += (gZbar2oMZ3*(-3 + 4*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", j, k))*Me[j] + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", j, k))*Me[j] + (sbar*(
mySMEFT.getSMEFTCoeffEW(
"CeBR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeBI", k, j)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CeWR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeWI", k, j)))*Me[k]))/(6.*sqrt(2));
653 CeuSRL.at(j).at(k).at(i).at(i) += (gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHeR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", j, k))*Me[j] - ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))*Me[k])*Mu[i])/4.;
654 CeuSRR.at(j).at(k).at(i).at(i) += (gZbar2oMZ4*
v2*(-((
mySMEFT.getSMEFTCoeffEW(
"CHeR", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHeI", j, k))*Me[j]) + ((
mySMEFT.getSMEFTCoeffEW(
"CHl1R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl1I", j, k)) + (
mySMEFT.getSMEFTCoeffEW(
"CHl3R", j, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHl3I", j, k)))*Me[k])*Mu[i])/4.;
656 for (
int p = 0; p < 3; p++)
657 for (
int r = 0; r < 3; r++)
659 CudV1LL.at(i).at(i).at(j).at(k) +=
VdLd(j, p) * (-0.08333333333333333 * (gZ2oMZ2 * (-3 + 4 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))))) *
VdL(r, k);
660 CudV1RR.at(i).at(i).at(j).at(k) +=
VdRd(j, p) * (-0.3333333333333333 * (gZ2oMZ2 * sbar2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r)))) *
VdR(r, k);
661 CudV1LR.at(i).at(i).at(j).at(k) +=
VdRd(j, p) * ((gZ2oMZ2 * (3 - 4 * sbar2) *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))) / 12.) *
VdR(r, k);
662 CduV1LR.at(j).at(k).at(i).at(i) +=
VdLd(j, p) * (-0.3333333333333333 * (gZ2oMZ2 * sbar2 *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))))) *
VdL(r, k);
664 CudV1LL.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * ((gZbar2oMZ3*(-3 + 4*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*Md[j])/(6.*sqrt(2))) *
VdL(r,k);
665 CudV1LL.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * ((gZbar2oMZ3*(-3 + 4*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[k])/(6.*sqrt(2))) *
VdR(r,k);
666 CudV1RR.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * ((sqrt(2)*gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*Md[k])/3.) *
VdL(r,k);
667 CudV1RR.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * ((sqrt(2)*gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[j])/3.) *
VdR(r,k);
668 CudV1LR.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * ((gZbar2oMZ3*(-3 + 4*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*Md[k])/(6.*sqrt(2))) *
VdL(r,k);
669 CudV1LR.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * ((gZbar2oMZ3*(-3 + 4*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[j])/(6.*sqrt(2))) *
VdR(r,k);
670 CduV1LR.at(j).at(k).at(i).at(i) +=
VdRd(j,p) * ((sqrt(2)*gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*Md[j])/3.) *
VdL(r,k);
671 CduV1LR.at(j).at(k).at(i).at(i) +=
VdLd(j,p) * ((sqrt(2)*gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*Md[k])/3.) *
VdR(r,k);
672 CudduV1LR.at(i).at(j).at(k).at(i) +=
VdLd(k,p) * ((gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[k]*Mu[i])/24.) *
VdL(r,j);
673 CudduV1LR.at(i).at(j).at(k).at(i) +=
VdRd(k,p) * (-0.041666666666666664*(gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*Md[j]*Mu[i])) *
VdR(r,j);
674 CudduV8LR.at(i).at(j).at(k).at(i) +=
VdLd(k,p) * ((gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[k]*Mu[i])/4.) *
VdL(r,j);
675 CudduV8LR.at(i).at(j).at(k).at(i) +=
VdRd(k,p) * (-0.25*(gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*Md[j]*Mu[i])) *
VdR(r,j);
676 CudS1RR.at(i).at(i).at(j).at(k) +=
VdRd(j,p) * (-0.25*(gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*Md[j]*Mu[i])) *
VdR(r,k);
677 CudS1RR.at(i).at(i).at(j).at(k) +=
VdLd(j,p) * ((gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[k]*Mu[i])/4.) *
VdL(r,k);
680 for (
int l = 0; l < 2; l++)
683 CudduS1RR.at(i).at(j).at(k).at(l) +=
VCKM(i,j) * VCKMd(k,l) * (-0.5*((1 + deltag2bar2oMW4)*g2bar2oMW4*Md[j]*Mu[l]));
685 for (
int p = 0; p < 3; p++)
686 for (
int r = 0; r < 3; r++)
688 CudduV1LR.at(i).at(j).at(k).at(l) +=
VCKM(i, j) *
VdRd(k, p) * (-0.25 * (g22oMW2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHudR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHudI", r, p)))) *
VuR(r, l);
690 CudduV1LR.at(i).at(j).at(k).at(l) +=
VCKM(i,j) *
VdRd(k,p) * (-((g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p))*Mu[l])/sqrt(2))) *
VuL(r,l);
691 CudduV1LR.at(i).at(j).at(k).at(l) +=
VCKM(i,j) *
VdLd(k,p) * (-((g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r))*Md[k])/sqrt(2))) *
VuR(r,l);
692 CudduS1RR.at(i).at(j).at(k).at(l) +=
VCKM(i,j) *
VdRd(k,p) * ((g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHudR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHudI", r, p))*Md[j]*Md[k])/4.) *
VuR(r,l);
693 CudduS1RR.at(i).at(j).at(k).at(l) +=
VCKM(i,j) *
VdLd(k,p) * (-0.5*(g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", r, p))*Md[j]*Mu[l])) *
VuL(r,l);
694 CudduS1RR.at(i).at(j).at(k).at(l) += VCKMd(k,l) *
VuLd(i,p) * (-0.5*(g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))*Md[j]*Mu[l])) *
VdL(r,j);
695 CudduS1RR.at(i).at(j).at(k).at(l) += VCKMd(k,l) *
VuRd(i,p) * ((g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHudR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHudI", p, r))*Mu[i]*Mu[l])/4.) *
VdR(r,j);
697 for (
int s = 0;
s < 3;
s++)
698 for (
int t = 0;
t < 3;
t++)
700 CudduS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i, p) *
VdLd(k,
s) * (-(
mySMEFT.getSMEFTCoeffEW(
"Cquqd1R",
s,
t, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cquqd1I",
s,
t, p, r))) *
VdR(r, j) *
VuR(
t, l);
701 CudduS8RR.at(i).at(j).at(k).at(l) +=
VuLd(i, p) *
VdLd(k,
s) * (-(
mySMEFT.getSMEFTCoeffEW(
"Cquqd8R",
s,
t, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cquqd8I",
s,
t, p, r))) *
VdR(r, j) *
VuR(
t, l);
703 CudduV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VdRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", r,
s)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)))*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
t, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
t, p))))/12.) *
VdL(r,j) *
VuR(
t,l);
704 CudduV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VdRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", r,
s)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)))*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
t, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
t, p))))/2.) *
VdL(r,j) *
VuR(
t,l);
712 for (
int i = 0; i < 2; i++)
713 for (
int j = 0; j < 2; j++)
715 CuuVLL.at(i).at(i).at(j).at(j) += -0.013888888888888888 * (gZ2oMZ2 * (-3 + 4 * sbar2) * (16 * delta_sbar * sbar2 + delta_gZ2oMZ2 * (-3 + 4 * sbar2)));
716 CuuVRR.at(i).at(i).at(j).at(j) += (-2 * (delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2) / 9.;
717 CuuV1LR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * (2 * delta_sbar * (3 - 8 * sbar2) + delta_gZ2oMZ2 * (3 - 4 * sbar2)) * sbar2) / 9.;
719 CuuV1LR.at(i).at(j).at(j).at(i) += -0.004629629629629629*(gZbar2oMZ4*(-3 + 4*sbar2)*(-3 + 4*(1 + 4*delta_sbar)*sbar2 + deltagZbar2oMZ4*(-3 + 4*sbar2))*Mu[i]*Mu[j]);
720 CuuV1LR.at(i).at(j).at(j).at(i) += (gZbar2oMZ4*sbar2*(-3 + 4*sbar2 + deltagZbar2oMZ4*(-3 + 4*sbar2) + 2*delta_sbar*(-3 + 8*sbar2))*Mu[i]*Mu[j])/54.;
721 CuuV1LR.at(i).at(j).at(j).at(i) += (gZbar2oMZ4*sbar2*(-3 + 4*sbar2 + deltagZbar2oMZ4*(-3 + 4*sbar2) + 2*delta_sbar*(-3 + 8*sbar2))*Mu[i]*Mu[j])/54.;
722 CuuV1LR.at(i).at(j).at(j).at(i) += (-2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Mu[i]*Mu[j])/27.;
723 CuuV8LR.at(i).at(j).at(j).at(i) += -0.027777777777777776*(gZbar2oMZ4*(-3 + 4*sbar2)*(-3 + 4*(1 + 4*delta_sbar)*sbar2 + deltagZbar2oMZ4*(-3 + 4*sbar2))*Mu[i]*Mu[j]);
724 CuuV8LR.at(i).at(j).at(j).at(i) += (gZbar2oMZ4*sbar2*(-3 + 4*sbar2 + deltagZbar2oMZ4*(-3 + 4*sbar2) + 2*delta_sbar*(-3 + 8*sbar2))*Mu[i]*Mu[j])/9.;
725 CuuV8LR.at(i).at(j).at(j).at(i) += (gZbar2oMZ4*sbar2*(-3 + 4*sbar2 + deltagZbar2oMZ4*(-3 + 4*sbar2) + 2*delta_sbar*(-3 + 8*sbar2))*Mu[i]*Mu[j])/9.;
726 CuuV8LR.at(i).at(j).at(j).at(i) += (-4*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Mu[i]*Mu[j])/9.;
727 CuuS1RR.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*sbar2*(-3 + 4*sbar2 + deltagZbar2oMZ4*(-3 + 4*sbar2) + 2*delta_sbar*(-3 + 8*sbar2))*Mu[i]*Mu[j])/18.;
728 CuuS1RR.at(i).at(i).at(j).at(j) += -0.013888888888888888*(gZbar2oMZ4*(-3 + 4*sbar2)*(-3 + 4*(1 + 4*delta_sbar)*sbar2 + deltagZbar2oMZ4*(-3 + 4*sbar2))*Mu[i]*Mu[j]);
729 CuuS1RR.at(i).at(i).at(j).at(j) += (-2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Mu[i]*Mu[j])/9.;
730 CuuS1RR.at(i).at(i).at(j).at(j) += (gZbar2oMZ4*sbar2*(-3 + 4*sbar2 + deltagZbar2oMZ4*(-3 + 4*sbar2) + 2*delta_sbar*(-3 + 8*sbar2))*Mu[i]*Mu[j])/18.;
732 for (
int k = 0; k < 3; k++)
734 for (
int l = 0; l < 3; l++)
736 CudV1LL.at(i).at(j).at(k).at(l) +=
VCKM(i, l) * VCKMd(k, j) * (-0.16666666666666666 * (delta_g22oMW2 * g22oMW2));
737 CudV8LL.at(i).at(j).at(k).at(l) +=
VCKM(i, l) * VCKMd(k, j) * (-(delta_g22oMW2 * g22oMW2));
739 CudV1LR.at(i).at(j).at(k).at(l) +=
VCKM(i,l) * VCKMd(k,j) * (-0.08333333333333333*((1 + deltag2bar2oMW4)*g2bar2oMW4*Md[k]*Md[l]));
740 CudV8LR.at(i).at(j).at(k).at(l) +=
VCKM(i,l) * VCKMd(k,j) * (-0.5*((1 + deltag2bar2oMW4)*g2bar2oMW4*Md[k]*Md[l]));
742 for (
int p = 0; p < 3; p++)
743 for (
int r = 0; r < 3; r++)
745 CudV1LL.at(i).at(j).at(k).at(l) += VCKMd(k, j) *
VuLd(i, p) * (-0.16666666666666666 * (g22oMW2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))) *
VdL(r, l);
746 CudV1LL.at(i).at(j).at(k).at(l) +=
VCKM(i, l) *
VdLd(k, p) * (-0.16666666666666666 * (g22oMW2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", r, p)))) *
VuL(r, j);
747 CudV8LL.at(i).at(j).at(k).at(l) += VCKMd(k, j) *
VuLd(i, p) * (-(g22oMW2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))) *
VdL(r, l);
748 CudV8LL.at(i).at(j).at(k).at(l) +=
VCKM(i, l) *
VdLd(k, p) * (-(g22oMW2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", r, p)))) *
VuL(r, j);
749 CueVLR.at(i).at(j).at(k).at(l) +=
VuLd(i, p) * ((
mySMEFT.getSMEFTCoeffEW(
"CqeR", p, r, k, l) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CqeI", p, r, k, l))) *
VuL(r, j);
751 CudV1LL.at(i).at(j).at(k).at(l) += VCKMd(k,j) *
VuRd(i,p) * ((g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p))*Mu[i])/(3.*sqrt(2))) *
VdL(r,l);
752 CudV1LL.at(i).at(j).at(k).at(l) += VCKMd(k,j) *
VuLd(i,p) * ((g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r))*Md[l])/(3.*sqrt(2))) *
VdR(r,l);
753 CudV1LL.at(i).at(j).at(k).at(l) +=
VCKM(i,l) *
VdRd(k,p) * (-0.3333333333333333*(g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p))*Md[k])/sqrt(2)) *
VuL(r,j);
754 CudV1LL.at(i).at(j).at(k).at(l) +=
VCKM(i,l) *
VdLd(k,p) * (-0.3333333333333333*(g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r))*Mu[j])/sqrt(2)) *
VuR(r,j);
755 CudV8LL.at(i).at(j).at(k).at(l) += VCKMd(k,j) *
VuRd(i,p) * (sqrt(2)*g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p))*Mu[i]) *
VdL(r,l);
756 CudV8LL.at(i).at(j).at(k).at(l) += VCKMd(k,j) *
VuLd(i,p) * (sqrt(2)*g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r))*Md[l]) *
VdR(r,l);
757 CudV8LL.at(i).at(j).at(k).at(l) +=
VCKM(i,l) *
VdRd(k,p) * (-(sqrt(2)*g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p))*Md[k])) *
VuL(r,j);
758 CudV8LL.at(i).at(j).at(k).at(l) +=
VCKM(i,l) *
VdLd(k,p) * (-(sqrt(2)*g2bar2oMW3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r))*Mu[j])) *
VuR(r,j);
759 CudV1LR.at(i).at(j).at(k).at(l) +=
VCKM(i,l) *
VdLd(k,p) * (-0.08333333333333333*(g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", r, p))*Md[k]*Md[l])) *
VuL(r,j);
760 CudV1LR.at(i).at(j).at(k).at(l) +=
VCKM(i,l) *
VdRd(k,p) * ((g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHudR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHudI", r, p))*Md[l]*Mu[j])/24.) *
VuR(r,j);
761 CudV1LR.at(i).at(j).at(k).at(l) += VCKMd(k,j) *
VuLd(i,p) * (-0.08333333333333333*(g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))*Md[k]*Md[l])) *
VdL(r,l);
762 CudV1LR.at(i).at(j).at(k).at(l) += VCKMd(k,j) *
VuRd(i,p) * ((g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHudR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHudI", p, r))*Md[k]*Mu[i])/24.) *
VdR(r,l);
763 CudV8LR.at(i).at(j).at(k).at(l) +=
VCKM(i,l) *
VdLd(k,p) * (-0.5*(g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", r, p))*Md[k]*Md[l])) *
VuL(r,j);
764 CudV8LR.at(i).at(j).at(k).at(l) +=
VCKM(i,l) *
VdRd(k,p) * ((g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHudR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHudI", r, p))*Md[l]*Mu[j])/4.) *
VuR(r,j);
765 CudV8LR.at(i).at(j).at(k).at(l) += VCKMd(k,j) *
VuLd(i,p) * (-0.5*(g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))*Md[k]*Md[l])) *
VdL(r,l);
766 CudV8LR.at(i).at(j).at(k).at(l) += VCKMd(k,j) *
VuRd(i,p) * ((g2bar2oMW4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHudR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHudI", p, r))*Md[k]*Mu[i])/4.) *
VdR(r,l);
768 for (
int s = 0;
s < 3;
s++)
769 for (
int t = 0;
t < 3;
t++)
771 CudV1LL.at(i).at(j).at(k).at(l) +=
VuLd(i, p) *
VdLd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"Cqq1R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqq1I", p, r,
s,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"Cqq1R",
s,
t, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqq1I",
s,
t, p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"Cqq3R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqq3I", p, r,
s,
t)) + (2 * (
mySMEFT.getSMEFTCoeffEW(
"Cqq3R", p,
t,
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqq3I", p,
t,
s, r))) / 3. + (2 * (
mySMEFT.getSMEFTCoeffEW(
"Cqq3R",
s, r, p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqq3I",
s, r, p,
t))) / 3. - (
mySMEFT.getSMEFTCoeffEW(
"Cqq3R",
s,
t, p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqq3I",
s,
t, p, r))) *
VuL(r, j) *
VdL(
t, l);
772 CudV8LL.at(i).at(j).at(k).at(l) +=
VuLd(i, p) *
VdLd(k,
s) * (4 * ((
mySMEFT.getSMEFTCoeffEW(
"Cqq3R", p,
t,
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqq3I", p,
t,
s, r)) + (
mySMEFT.getSMEFTCoeffEW(
"Cqq3R",
s, r, p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqq3I",
s, r, p,
t)))) *
VuL(r, j) *
VdL(
t, l);
773 CudV1RR.at(i).at(j).at(k).at(l) +=
VuRd(i, p) *
VdRd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"Cud1R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cud1I", p, r,
s,
t))) *
VuR(r, j) *
VdR(
t, l);
774 CudV8RR.at(i).at(j).at(k).at(l) +=
VuRd(i, p) *
VdRd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"Cud8R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cud8I", p, r,
s,
t))) *
VuR(r, j) *
VdR(
t, l);
775 CudV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i, p) *
VdRd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"Cqd1R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqd1I", p, r,
s,
t))) *
VuL(r, j) *
VdR(
t, l);
776 CudV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i, p) *
VdRd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"Cqd8R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqd8I", p, r,
s,
t))) *
VuL(r, j) *
VdR(
t, l);
777 CudS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i, p) *
VdLd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"Cquqd1R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cquqd1I", p, r,
s,
t))) *
VuR(r, j) *
VdR(
t, l);
778 CudS8RR.at(i).at(j).at(k).at(l) +=
VuLd(i, p) *
VdLd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"Cquqd8R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cquqd8I", p, r,
s,
t))) *
VuR(r, j) *
VdR(
t, l);
780 CudS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VdLd(k,
s) * ((oneoMh2*(-3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p, r))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
t,
s)) + (-3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI",
s,
t)) + 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
t,
s)))*(
mySMEFT.getSMEFTCoeffEW(
"YuR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI", r, p))))/2.) *
VuR(r,j) *
VdR(
t,l);
786 for (
int k = 0; k < 2; k++)
788 for (
int p = 0; p < 3; p++)
789 for (
int r = 0; r < 3; r++)
791 CuuVLL.at(i).at(i).at(j).at(k) +=
VuLd(j, p) * (-0.041666666666666664 * (gZ2oMZ2 * (-3 + 4 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))))) *
VuL(r, k);
792 CuuVLL.at(j).at(k).at(i).at(i) +=
VuLd(j, p) * (-0.041666666666666664 * (gZ2oMZ2 * (-3 + 4 * sbar2) *
v2 * ((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))))) *
VuL(r, k);
793 CuuVRR.at(i).at(i).at(j).at(k) +=
VuRd(j, p) * (-0.16666666666666666 * (gZ2oMZ2 * sbar2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r)))) *
VuR(r, k);
794 CuuVRR.at(j).at(k).at(i).at(i) +=
VuRd(j, p) * (-0.16666666666666666 * (gZ2oMZ2 * sbar2 *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r)))) *
VuR(r, k);
795 CuuV1LR.at(i).at(i).at(j).at(k) +=
VuRd(j, p) * ((gZ2oMZ2 * (3 - 4 * sbar2) *
v2 * (
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))) / 12.) *
VuR(r, k);
796 CuuV1LR.at(j).at(k).at(i).at(i) +=
VuLd(j, p) * ((gZ2oMZ2 * sbar2 *
v2 * (-(
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))) / 3.) *
VuL(r, k);
798 CuuVLL.at(i).at(i).at(j).at(k) +=
VuRd(j,p) * ((gZbar2oMZ3*(-3 + 4*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[j])/(12.*sqrt(2))) *
VuL(r,k);
799 CuuVLL.at(i).at(i).at(j).at(k) +=
VuLd(j,p) * ((gZbar2oMZ3*(-3 + 4*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[k])/(12.*sqrt(2))) *
VuR(r,k);
800 CuuVLL.at(j).at(k).at(i).at(i) +=
VuRd(j,p) * ((gZbar2oMZ3*(-3 + 4*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[j])/(12.*sqrt(2))) *
VuL(r,k);
801 CuuVLL.at(j).at(k).at(i).at(i) +=
VuLd(j,p) * ((gZbar2oMZ3*(-3 + 4*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[k])/(12.*sqrt(2))) *
VuR(r,k);
802 CuuVRR.at(i).at(i).at(j).at(k) +=
VuRd(j,p) * ((gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[k])/(3.*sqrt(2))) *
VuL(r,k);
803 CuuVRR.at(i).at(i).at(j).at(k) +=
VuLd(j,p) * ((gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[j])/(3.*sqrt(2))) *
VuR(r,k);
804 CuuVRR.at(j).at(k).at(i).at(i) +=
VuRd(j,p) * ((gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[k])/(3.*sqrt(2))) *
VuL(r,k);
805 CuuVRR.at(j).at(k).at(i).at(i) +=
VuLd(j,p) * ((gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[j])/(3.*sqrt(2))) *
VuR(r,k);
806 CuuV1LR.at(i).at(j).at(k).at(i) +=
VuLd(k,p) * ((gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Mu[i]*Mu[k])/24.) *
VuL(r,j);
807 CuuV1LR.at(i).at(j).at(k).at(i) +=
VuRd(k,p) * (-0.041666666666666664*(gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*Mu[i]*Mu[j])) *
VuR(r,j);
808 CuuV1LR.at(j).at(i).at(i).at(k) +=
VuLd(j,p) * ((gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Mu[i]*Mu[k])/24.) *
VuL(r,k);
809 CuuV1LR.at(j).at(i).at(i).at(k) +=
VuRd(j,p) * (-0.041666666666666664*(gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*Mu[i]*Mu[j])) *
VuR(r,k);
810 CuuV1LR.at(i).at(i).at(j).at(k) +=
VuRd(j,p) * ((gZbar2oMZ3*(-3 + 4*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[k])/(6.*sqrt(2))) *
VuL(r,k);
811 CuuV1LR.at(i).at(i).at(j).at(k) +=
VuLd(j,p) * ((gZbar2oMZ3*(-3 + 4*sbar2)*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[j])/(6.*sqrt(2))) *
VuR(r,k);
812 CuuV1LR.at(j).at(k).at(i).at(i) +=
VuRd(j,p) * ((sqrt(2)*gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[j])/3.) *
VuL(r,k);
813 CuuV1LR.at(j).at(k).at(i).at(i) +=
VuLd(j,p) * ((sqrt(2)*gZbar2oMZ3*sbar2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[k])/3.) *
VuR(r,k);
814 CuuV8LR.at(i).at(j).at(k).at(i) +=
VuLd(k,p) * ((gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Mu[i]*Mu[k])/4.) *
VuL(r,j);
815 CuuV8LR.at(i).at(j).at(k).at(i) +=
VuRd(k,p) * (-0.25*(gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*Mu[i]*Mu[j])) *
VuR(r,j);
816 CuuV8LR.at(j).at(i).at(i).at(k) +=
VuLd(j,p) * ((gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Mu[i]*Mu[k])/4.) *
VuL(r,k);
817 CuuV8LR.at(j).at(i).at(i).at(k) +=
VuRd(j,p) * (-0.25*(gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*Mu[i]*Mu[j])) *
VuR(r,k);
818 CuuS1RR.at(i).at(i).at(j).at(k) +=
VuRd(j,p) * (-0.125*(gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*Mu[i]*Mu[j])) *
VuR(r,k);
819 CuuS1RR.at(i).at(i).at(j).at(k) +=
VuLd(j,p) * ((gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Mu[i]*Mu[k])/8.) *
VuL(r,k);
820 CuuS1RR.at(j).at(k).at(i).at(i) +=
VuLd(j,p) * ((gZbar2oMZ4*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Mu[i]*Mu[k])/8.) *
VuL(r,k);
821 CuuS1RR.at(j).at(k).at(i).at(i) +=
VuRd(j,p) * (-0.125*(gZbar2oMZ4*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*Mu[i]*Mu[j])) *
VuR(r,k);
824 for (
int l = 0; l < 2; l++)
825 for (
int p = 0; p < 3; p++)
826 for (
int r = 0; r < 3; r++)
827 for (
int s = 0;
s < 3;
s++)
828 for (
int t = 0;
t < 3;
t++)
830 CuuVLL.at(i).at(j).at(k).at(l) +=
VuLd(i, p) *
VuLd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"Cqq1R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqq1I", p, r,
s,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"Cqq3R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqq3I", p, r,
s,
t))) *
VuL(r, j) *
VuL(
t, l);
831 CuuVRR.at(i).at(j).at(k).at(l) +=
VuRd(i, p) *
VuRd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"CuuR", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuuI", p, r,
s,
t))) *
VuR(r, j) *
VuR(
t, l);
832 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i, p) *
VuRd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"Cqu1R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqu1I", p, r,
s,
t))) *
VuL(r, j) *
VuR(
t, l);
833 CuuV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i, p) *
VuRd(k,
s) * ((
mySMEFT.getSMEFTCoeffEW(
"Cqu8R", p, r,
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"Cqu8I", p, r,
s,
t))) *
VuL(r, j) *
VuR(
t, l);
834 CuuVLL.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * ((gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t))*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)) - gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)) + gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r))*(-(
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t))))/8.) *
VuL(r,j) *
VuL(
t,l);
835 CuuVRR.at(i).at(j).at(k).at(l) +=
VuRd(i,p) *
VuRd(k,
s) * (-0.125*(gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t)))) *
VuR(r,j) *
VuR(
t,l);
836 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * (-0.25*(gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t))) + (gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t)))/4.) *
VuL(r,j) *
VuR(
t,l);
838 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", r,
s)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
s, r)))*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
t, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
t, p))))/12.) *
VuL(r,j) *
VuR(
t,l);
839 CuuV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", r,
s)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
s, r)))*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
t, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
t, p))))/2.) *
VuL(r,j) *
VuR(
t,l);
840 CuuS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * ((oneoMh2*(-3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI",
s,
t))*(
mySMEFT.getSMEFTCoeffEW(
"YuR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI", r, p)) + (-3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p, r)) + 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YuR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI", r, p)))*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
t,
s))))/4.) *
VuR(r,j) *
VuR(
t,l);
841 CuuVLL.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI",
t,
s)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI",
t,
s)))*Mu[k])/(4.*sqrt(2))) *
VuL(r,j) *
VuL(
t,l);
842 CuuVLL.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI",
s,
t)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI",
s,
t)))*Mu[l])/(4.*sqrt(2))) *
VuL(r,j) *
VuR(
t,l);
843 CuuVLL.at(i).at(j).at(k).at(l) +=
VuRd(i,p) *
VuLd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[i])/(4.*sqrt(2))) *
VuL(r,j) *
VuL(
t,l);
844 CuuVLL.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[j])/(4.*sqrt(2))) *
VuR(r,j) *
VuL(
t,l);
845 CuuVRR.at(i).at(j).at(k).at(l) +=
VuRd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI",
t,
s)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI",
t,
s)))*Mu[l])/(4.*sqrt(2)))*
VuR(r,j) *
VuL(
t,l);
846 CuuVRR.at(i).at(j).at(k).at(l) +=
VuRd(i,p) *
VuLd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI",
s,
t)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI",
s,
t)))*Mu[k])/(4.*sqrt(2))) *
VuR(r,j) *
VuR(
t,l);
847 CuuVRR.at(i).at(j).at(k).at(l) +=
VuRd(i,p) *
VuRd(k,
s) *((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[j])/(4.*sqrt(2))) *
VuL(r,j) *
VuR(
t,l);
848 CuuVRR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) *((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p ,r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p ,r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p ,r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p ,r)))*Mu[i])/(4.*sqrt(2))) *
VuR(r,j) *
VuR(
t,l);
849 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * (-0.041666666666666664*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p,
t))*(3*
mySMEFT.FlagQuadraticWC*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", r,
s)))))) *
VuL(r,j) *
VuR(
t,l);
850 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI",
t,
s)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI",
t,
s)))*Mu[l])/(2.*sqrt(2))) *
VuL(r,j) *
VuR(
t,l);
851 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI",
s,
t)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI",
s,
t)))*Mu[k])/(2.*sqrt(2))) *
VuL(r,j) *
VuR(
t,l);
852 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[i])/(2.*sqrt(2))) *
VuL(r,j) *
VuR(
t,l);
853 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[j])/(2.*sqrt(2))) *
VuL(r,j) *
VuR(
t,l);
854 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * (-0.041666666666666664*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s, r)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p,
t)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p,
t)))*Mu[k]*Mu[l])) *
VuL(r,j) *
VuR(
t,l);
855 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p,
t)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p,
t)))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s, r))*Mu[j]*Mu[l])/24.) *
VuL(r,j) *
VuR(
t,l);
856 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s, r)))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p,
t))*Mu[i]*Mu[k])/24.) *
VuL(r,j) *
VuR(
t,l);
857 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * (-0.041666666666666664*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s, r))*Mu[i]*Mu[j])) *
VuL(r,j) *
VuR(
t,l);
858 CuuV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * (-0.25*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p,
t))*(3*
mySMEFT.FlagQuadraticWC*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", r,
s)))))) *
VuL(r,j) *
VuR(
t,l);
859 CuuV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * (-0.25*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s, r)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p,
t)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p,
t)))*Mu[k]*Mu[l])) *
VuL(r,j) *
VuR(
t,l);
860 CuuV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p,
t)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p,
t)))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s, r))*Mu[j]*Mu[l])/4.) *
VuL(r,j) *
VuR(
t,l);
861 CuuV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s, r)))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p,
t))*Mu[i]*Mu[k])/4.) *
VuL(r,j) *
VuR(
t,l);
862 CuuV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * (-0.25*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s, r))*Mu[i]*Mu[j])) *
VuL(r,j) *
VuR(
t,l);
863 CuuS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p, r))*(3*
mySMEFT.FlagQuadraticWC*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI",
s,
t)))))/8.) *
VuR(r,j) *
VuR(
t,l);
864 CuuS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t))*Mu[j]*Mu[k])/8.) *
VuR(r,j) *
VuR(
t,l);
865 CuuS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * (-0.125*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))*Mu[j]*Mu[l])) *
VuR(r,j) *
VuR(
t,l);
866 CuuS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * (-0.125*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p ,r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p ,r))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t))*Mu[i]*Mu[k])) *
VuR(r,j) *
VuR(
t,l);
867 CuuS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p ,r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p ,r))*Mu[i]*Mu[l])/8.) *
VuR(r,j) *
VuR(
t,l);
872 for (
int p = 0; p < 3; p++)
873 for (
int r = 0; r < 3; r++)
875 Cug.at(i).at(j) += vTosq2 *
VuLd(i, p) * ((
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)) * sbar + (
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) * cbar) *
VuR(r, j);
876 CuG.at(i).at(j) += vTosq2 *
VuLd(i, p) * (
mySMEFT.getSMEFTCoeffEW(
"CuGR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuGI", p, r)) *
VuR(r, j);
const gslpp::matrix< gslpp::complex > getCKM() const
A member for returning the CKM matrix.
double getLambda_NP() const
Return Lambda_NP.
const std::array< std::array< std::array< std::array< gslpp::complex, 2 >, 2 >, 2 >, 2 > zero2222
std::array< std::array< gslpp::complex, 3 >, 3 > Ceg
The real part of the dimension-5 operator coefficient .
const std::array< std::array< std::array< std::array< gslpp::complex, 2 >, 2 >, 3 >, 3 > zero3322
std::array< std::array< gslpp::complex, 2 >, 2 > Cug
The real part of the dimension-5 operator coefficient .
const std::array< std::array< std::array< std::array< gslpp::complex, 3 >, 3 >, 2 >, 2 > zero2233
std::array< std::array< gslpp::complex, 2 >, 2 > CuG
The real part of the dimension-5 operator coefficient .
const std::array< std::array< std::array< std::array< gslpp::complex, 2 >, 3 >, 3 >, 3 > zero3332
const std::array< std::array< std::array< std::array< gslpp::complex, 3 >, 3 >, 3 >, 3 > zero3333
const std::array< std::array< gslpp::complex, 2 >, 2 > zero22
const std::array< std::array< gslpp::complex, 3 >, 3 > zero33
const std::array< std::array< std::array< std::array< gslpp::complex, 2 >, 3 >, 3 >, 2 > zero2332
const double getMz() const
A get method to access the mass of the boson .
const double getMw() const
A get method to access the input value of the mass of the boson .
const double v() const
The Higgs vacuum expectation value.
void updateSMParameters()
Updates to new Standard Model parameter sets.
A class for the CKM elements .
master
Inheritance diagram for NPSMEFTd6GeneralMatching:
1.9.2 