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 = 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", l, i) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", l, i))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", k, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", k, j)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", k, j) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", k, j)) - 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", l, k) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", l, k))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", j, i) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", j, i)) - 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", j, i) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", j, i))*(
mySMEFT.getSMEFTCoeffEW(
"YdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI", p, r)) + (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", 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", j, i) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", j, i))*(
mySMEFT.getSMEFTCoeffEW(
"YdR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI", r, p)) + (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", 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);
430 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR",
t, p) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI",
t, p))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR",
s, r) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI",
s, r)) - 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);
431 CddV8LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR",
t, p) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI",
t, p))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR",
s, r) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI",
s, r)) - 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);
432 CddS1RR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * (-0.25*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR",
t,
s) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI",
t,
s))*(
mySMEFT.getSMEFTCoeffEW(
"YdR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI", r, p)) + (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", 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);
437 for (
int l = 0; l < 2; l++)
439 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))));
441 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));
442 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]));
444 for (
int p = 0; p < 3; p++)
445 for (
int r = 0; r < 3; r++)
447 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);
448 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);
449 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);
450 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);
454 for (
int k = 0; k < 2; k++)
456 CnueduVLL.at(i).at(i).at(j).at(k) += VCKMd(j, k) * (-0.5 * (delta_g22oMW2 * g22oMW2));
458 CnueduSRR.at(i).at(i).at(j).at(k) += VCKMd(j,k) * (-0.5*((1 + deltag2bar2oMW4)*g2bar2oMW4*Me[i]*Mu[k]));
460 for (
int p = 0; p < 3; p++)
461 for (
int r = 0; r < 3; r++)
463 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);
464 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);
466 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);
467 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);
468 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);
469 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);
470 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);
471 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);
472 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);
473 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);
476 for (
int l = 0; l < 2; l++)
479 CduV1LR.at(i).at(j).at(k).at(l) += VCKMd(i,l) *
VCKM(k,j) * (-0.08333333333333333*((1 + deltag2bar2oMW4)*g2bar2oMW4*Mu[k]*Mu[l]));
480 CduV8LR.at(i).at(j).at(k).at(l) += VCKMd(i,l) *
VCKM(k,j) * (-0.5*((1 + deltag2bar2oMW4)*g2bar2oMW4*Mu[k]*Mu[l]));
482 for (
int p = 0; p < 3; p++)
483 for (
int r = 0; r < 3; r++)
485 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);
486 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);
487 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);
488 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);
489 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);
490 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);
491 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);
493 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);
494 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);
495 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);
496 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);
497 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);
498 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);
499 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);
500 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);
501 CeuSRL.at(i).at(j).at(k).at(l) +=
VuRd(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", j, i) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", j, i)) - 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);
502 CeuSRR.at(i).at(j).at(k).at(l) +=
VuLd(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", j, i) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", j, i)) - 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);
504 for (
int s = 0;
s < 3;
s++)
505 for (
int t = 0;
t < 3;
t++)
507 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);
508 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);
513 for (
int p = 0; p < 3; p++)
514 for (
int r = 0; r < 3; r++)
516 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);
517 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);
521 for (
int i = 0; i < 3; i++)
522 for (
int j = 0; j < 2; j++)
524 CnuuVLL.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * (8 * delta_sbar * sbar2 + delta_gZ2oMZ2 * (-3 + 4 * sbar2))) / 12.;
525 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.;
526 CeuVRR.at(i).at(i).at(j).at(j) += (2 * (delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2) / 3.;
527 CnuuVLR.at(i).at(i).at(j).at(j) += ((delta_gZ2oMZ2 + 2 * delta_sbar) * gZ2oMZ2 * sbar2) / 3.;
528 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.;
529 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.;
531 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]);
532 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.;
533 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.;
534 CeuSRL.at(i).at(i).at(j).at(j) += (-2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Me[i]*Mu[j])/3.;
535 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.;
536 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.;
537 CeuSRR.at(i).at(i).at(j).at(j) += (2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Me[i]*Mu[j])/3.;
538 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.;
540 for (
int k = 0; k < 2; k++)
541 for (
int p = 0; p < 3; p++)
542 for (
int r = 0; r < 3; r++)
544 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);
545 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);
546 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);
547 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);
548 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);
549 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);
550 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);
551 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);
552 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);
553 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);
555 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);
556 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);
557 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);
558 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);
559 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);
560 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);
561 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);
562 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);
563 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);
564 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);
565 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);
566 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);
567 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);
568 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);
569 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);
570 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);
571 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);
572 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);
573 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);
574 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);
575 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);
576 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);
577 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);
578 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);
579 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);
580 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);
581 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);
582 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);
583 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);
584 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);
589 for (
int i = 0; i < 2; i++)
590 for (
int j = 0; j < 3; j++)
592 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.;
593 CudV1RR.at(i).at(i).at(j).at(j) += (2 * (delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2) / 9.;
594 CudV1LR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * sbar2 * (delta_gZ2oMZ2 * (-3 + 4 * sbar2) + 2 * delta_sbar * (-3 + 8 * sbar2))) / 18.;
595 CueVLR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * sbar2 * (delta_gZ2oMZ2 * (-3 + 4 * sbar2) + 2 * delta_sbar * (-3 + 8 * sbar2))) / 6.;
597 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.;
598 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]);
599 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]);
600 CudduV1LR.at(i).at(j).at(j).at(i) += ((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[j]*Mu[i])/27.;
601 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.;
602 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]);
603 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]);
604 CudduV8LR.at(i).at(j).at(j).at(i) += (2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[j]*Mu[i])/9.;
605 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]);
606 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.;
607 CudS1RR.at(i).at(i).at(j).at(j) += (2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[j]*Mu[i])/9.;
608 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]);
610 for (
int k = 0; k < 3; k++)
612 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))));
613 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))));
614 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)));
615 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.;
616 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))));
617 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.;
619 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));
620 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.;
621 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.;
622 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));
623 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.;
624 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.;
626 for (
int p = 0; p < 3; p++)
627 for (
int r = 0; r < 3; r++)
629 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);
630 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);
631 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);
632 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);
634 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);
635 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);
636 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);
637 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);
638 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);
639 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);
640 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);
641 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);
642 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);
643 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);
644 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);
645 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);
646 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);
647 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);
650 for (
int l = 0; l < 2; l++)
653 CudduS1RR.at(i).at(j).at(k).at(l) +=
VCKM(i,j) * VCKMd(k,l) * (-0.5*((1 + deltag2bar2oMW4)*g2bar2oMW4*Md[j]*Mu[l]));
655 for (
int p = 0; p < 3; p++)
656 for (
int r = 0; r < 3; r++)
658 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);
660 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);
661 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);
662 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);
663 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);
664 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);
665 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);
667 for (
int s = 0;
s < 3;
s++)
668 for (
int t = 0;
t < 3;
t++)
670 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);
671 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);
673 CudduV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VdRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR",
t, p) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI",
t, p))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR",
s, r) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI",
s, r)) - 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);
674 CudduV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VdRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR",
t, p) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI",
t, p))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR",
s, r) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI",
s, r)) - 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);
675 CudduS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VdLd(k,
s) * (-0.25*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR",
t,
s) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI",
t,
s))*(
mySMEFT.getSMEFTCoeffEW(
"YdR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI", r, p)) + (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", 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) *
VuR(
t,l);
683 for (
int i = 0; i < 2; i++)
684 for (
int j = 0; j < 2; j++)
686 CuuVLL.at(i).at(i).at(j).at(j) += -0.013888888888888888 * (gZ2oMZ2 * (-3 + 4 * sbar2) * (16 * delta_sbar * sbar2 + delta_gZ2oMZ2 * (-3 + 4 * sbar2)));
687 CuuVRR.at(i).at(i).at(j).at(j) += (-2 * (delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2) / 9.;
688 CuuV1LR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * (2 * delta_sbar * (3 - 8 * sbar2) + delta_gZ2oMZ2 * (3 - 4 * sbar2)) * sbar2) / 9.;
690 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]);
691 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.;
692 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.;
693 CuuV1LR.at(i).at(j).at(j).at(i) += (-2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Mu[i]*Mu[j])/27.;
694 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]);
695 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.;
696 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.;
697 CuuV8LR.at(i).at(j).at(j).at(i) += (-4*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Mu[i]*Mu[j])/9.;
698 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.;
699 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]);
700 CuuS1RR.at(i).at(i).at(j).at(j) += (-2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Mu[i]*Mu[j])/9.;
701 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.;
703 for (
int k = 0; k < 3; k++)
705 for (
int l = 0; l < 3; l++)
707 CudV1LL.at(i).at(j).at(k).at(l) +=
VCKM(i, l) * VCKMd(k, j) * (-0.16666666666666666 * (delta_g22oMW2 * g22oMW2));
708 CudV8LL.at(i).at(j).at(k).at(l) +=
VCKM(i, l) * VCKMd(k, j) * (-(delta_g22oMW2 * g22oMW2));
710 CudV1LR.at(i).at(j).at(k).at(l) +=
VCKM(i,l) * VCKMd(k,j) * (-0.08333333333333333*((1 + deltag2bar2oMW4)*g2bar2oMW4*Md[k]*Md[l]));
711 CudV8LR.at(i).at(j).at(k).at(l) +=
VCKM(i,l) * VCKMd(k,j) * (-0.5*((1 + deltag2bar2oMW4)*g2bar2oMW4*Md[k]*Md[l]));
713 for (
int p = 0; p < 3; p++)
714 for (
int r = 0; r < 3; r++)
716 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);
717 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);
718 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);
719 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);
720 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);
722 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);
723 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);
724 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);
725 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);
726 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);
727 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);
728 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);
729 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);
730 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);
731 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);
732 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);
733 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);
734 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);
735 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);
736 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);
737 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);
739 for (
int s = 0;
s < 3;
s++)
740 for (
int t = 0;
t < 3;
t++)
742 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);
743 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);
744 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);
745 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);
746 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);
747 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);
748 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);
749 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);
751 CudS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VdLd(k,
s) * ((oneoMh2*(-3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", r, p) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", r, p))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
t,
s)) + (-3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR",
t,
s) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI",
t,
s)) + 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);
757 for (
int k = 0; k < 2; k++)
759 for (
int p = 0; p < 3; p++)
760 for (
int r = 0; r < 3; r++)
762 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);
763 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);
764 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);
765 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);
766 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);
767 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);
769 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);
770 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);
771 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);
772 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);
773 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);
774 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);
775 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);
776 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);
777 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);
778 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);
779 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);
780 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);
781 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);
782 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);
783 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);
784 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);
785 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);
786 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);
787 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);
788 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);
789 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);
790 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);
791 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);
792 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);
795 for (
int l = 0; l < 2; l++)
796 for (
int p = 0; p < 3; p++)
797 for (
int r = 0; r < 3; r++)
798 for (
int s = 0;
s < 3;
s++)
799 for (
int t = 0;
t < 3;
t++)
801 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);
802 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);
803 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);
804 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);
806 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR",
t, p) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI",
t, p))*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR",
s, r) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI",
s, r)) - 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);
807 CuuV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR",
t, p) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI",
t, p))*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR",
s, r) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI",
s, r)) - 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);
808 CuuS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * ((oneoMh2*(-3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR",
t,
s) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI",
t,
s))*(
mySMEFT.getSMEFTCoeffEW(
"YuR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI", r, p)) + (-3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", r, p) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", r, p)) + 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);
813 for (
int p = 0; p < 3; p++)
814 for (
int r = 0; r < 3; r++)
816 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);
817 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 