54 double delta_vT =
mySMEFT.getDelta_v();
55 double vTosq2 = vT / sqrt(2.);
103 double cbar =
mySMEFT.getXWZ_tree();
104 double sbar = -
mySMEFT.getXBZ_tree();
105 double sbar2 = sbar * sbar;
107 double delta_sbar =
mySMEFT.getDelta_xBZ();
108 double g1bar =
mySMEFT.getG1_tree();
110 double g2bar =
mySMEFT.getG2_tree();
112 double delta_MZ2 =
mySMEFT.getDelta_Mz2();
113 double ebar =
mySMEFT.getEeMz();
116 double gZbar = ebar / sbar / cbar;
117 double delta_gZbar = (g1bar * g1bar + g2bar * g2bar) / (2. * g1bar * g2bar) *
v2 *
mySMEFT.getSMEFTCoeffEW(
"CHWB");
119 delta_gZbar += -(g1bar * g1bar * g1bar * g1bar + g2bar * g2bar * g2bar * g2bar) / (2. * (g1bar * g1bar + g2bar * g2bar) * g1bar * g2bar) *
v2 *
mySMEFT.getSMEFTCoeffEW(
"CHWB");
122 double delta_gZ2oMZ2 = 2. * delta_gZbar - delta_MZ2;
123 double g22oMW2 = 4. /
v2;
124 double delta_g22oMW2 = -2. * delta_vT;
129 double deltagZbar2oMZ4 = 2. * delta_gZbar - 2. * delta_MZ2;
133 double deltag2bar2oMW4 = -4. * delta_vT;
135 double lambda =
mySMEFT.getLambdaH_tree() * 2.;
136 double oneoMh2 =
mySMEFT.FlagHiggsExch * 1. / (lambda *
v2);
137 double deltaoneoMh2 = (- 2. *
mySMEFT.getSMEFTCoeffEW(
"CHbox") + 0.5 *
mySMEFT.getSMEFTCoeffEW(
"CHD") + 3 *
mySMEFT.getSMEFTCoeffEW(
"CH") / lambda) *
v2;
139 std::array<double, 3> Me = {
144 std::array<double, 3> Mu = {
148 std::array<double, 3> Md = {
163 gslpp::matrix<complex> VCKMd =
VCKM.hconjugate();
167 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.)
168 throw(
"Compiler is not putting to zero correctly the 2-d arrays of Wilson coefficients");
231#ifdef NOLEPTONFLAVOURVIOLATION
235 for (
int i = 0; i < 3; i++)
236 for (
int j = 0; j < 3; j++)
238 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);
240 CnunuVLL.at(i).at(i).at(j).at(j) += -0.0625 * (delta_gZ2oMZ2 * gZ2oMZ2);
241 CnunuVLL.at(i).at(j).at(j).at(i) += -0.0625 * (delta_gZ2oMZ2 * gZ2oMZ2);
242 CeeVLL.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * (-1 + 2 * sbar2) * (delta_gZ2oMZ2 - 2 * delta_gZ2oMZ2 * sbar2 - 8 * delta_sbar * sbar2)) / 16.;
243 CeeVLL.at(i).at(j).at(j).at(i) += (gZ2oMZ2 * (-1 + 2 * sbar2) * (delta_gZ2oMZ2 - 2 * delta_gZ2oMZ2 * sbar2 - 8 * delta_sbar * sbar2)) / 16.;
244 CnueVLL.at(i).at(i).at(j).at(j) += -0.25 * (gZ2oMZ2 * (4 * delta_sbar * sbar2 + delta_gZ2oMZ2 * (-1 + 2 * sbar2)));
245 CnueVLL.at(i).at(j).at(j).at(i) += -0.5 * (delta_g22oMW2 * g22oMW2);
246 CnudVLL.at(i).at(i).at(j).at(j) += -0.08333333333333333 * (gZ2oMZ2 * (4 * delta_sbar * sbar2 + delta_gZ2oMZ2 * (-3 + 2 * sbar2)));
247 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)));
248 CddVLL.at(i).at(i).at(j).at(j) += -0.013888888888888888 * (gZ2oMZ2 * (-3 + 2 * sbar2) * (8 * delta_sbar * sbar2 + delta_gZ2oMZ2 * (-3 + 2 * sbar2)));
249 CeeVRR.at(i).at(i).at(j).at(j) += -0.25 * ((delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2);
250 CeeVRR.at(i).at(j).at(j).at(i) += -0.25 * ((delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2);
251 CedVRR.at(i).at(i).at(j).at(j) += -0.3333333333333333 * ((delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2);
252 CddVRR.at(i).at(i).at(j).at(j) += -0.05555555555555555 * ((delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2);
253 CnueVLR.at(i).at(i).at(j).at(j) += -0.5 * ((delta_gZ2oMZ2 + 2 * delta_sbar) * gZ2oMZ2 * sbar2);
254 CeeVLR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * sbar2 * (delta_gZ2oMZ2 + delta_sbar * (2 - 8 * sbar2) - 2 * delta_gZ2oMZ2 * sbar2)) / 2.;
255 CnudVLR.at(i).at(i).at(j).at(j) += -0.16666666666666666 * ((delta_gZ2oMZ2 + 2 * delta_sbar) * gZ2oMZ2 * sbar2);
256 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));
257 CdeVLR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * (delta_sbar * (6 - 8 * sbar2) + delta_gZ2oMZ2 * (3 - 2 * sbar2)) * sbar2) / 6.;
258 CddV1LR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * (delta_sbar * (6 - 8 * sbar2) + delta_gZ2oMZ2 * (3 - 2 * sbar2)) * sbar2) / 18.;
260 CnueVLR.at(i).at(j).at(j).at(i) += -0.25*((1 + deltag2bar2oMW4)*g2bar2oMW4*Me[i]*Me[j]);
261 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]);
262 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.;
263 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.;
264 CeeVLR.at(i).at(j).at(j).at(i) += -0.5*((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Me[i]*Me[j]);
265 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]);
266 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.;
267 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.;
268 CddV1LR.at(i).at(j).at(j).at(i) += -0.018518518518518517*((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[i]*Md[j]);
269 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]);
270 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.;
271 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.;
272 CddV8LR.at(i).at(j).at(j).at(i) += -0.1111111111111111*((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[i]*Md[j]);
273 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.;
274 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.;
275 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.;
276 CedSRL.at(i).at(i).at(j).at(j) += ((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[j]*Me[i])/3.;
277 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.;
278 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]);
279 CeeSRR.at(i).at(i).at(j).at(j) += -0.5*((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Me[i]*Me[j]);
280 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.;
281 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.;
282 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]);
283 CedSRR.at(i).at(i).at(j).at(j) += -0.3333333333333333*((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[j]*Me[i]);
284 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.;
285 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.;
286 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]);
287 CddS1RR.at(i).at(i).at(j).at(j) += -0.1111111111111111*((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[i]*Md[j]);
288 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.;
290 for (
int k = 0; k < 3; k++)
292 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.;
293 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.;
294 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.;
295 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.;
296 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.;
297 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.;
298 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.;
299 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.;
300 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.;
301 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.;
302 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)));
303 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)));
304 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.;
305 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.;
306 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.;
307 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.);
308 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.;
309 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.;
310 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.;
311 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.;
312 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.;
313 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.;
314 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.;
315 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.;
316 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.;
317 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.;
319 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));
320 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);
321 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);
322 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);
323 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);
324 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);
325 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));
326 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);
327 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);
328 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);
329 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);
330 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);
331 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);
332 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]);
333 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]);
334 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);
335 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.;
336 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.;
337 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);
338 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));
339 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);
340 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);
341 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.;
342 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.;
343 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.;
344 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.;
346 for (
int p = 0; p < 3; p++)
347 for (
int r = 0; r < 3; r++)
349 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);
350 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);
351 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);
352 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);
353 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);
354 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);
355 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);
356 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);
357 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);
358 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);
359 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);
360 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);
361 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);
363 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);
364 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);
365 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);
366 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);
367 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);
368 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);
369 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);
370 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);
371 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);
372 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);
373 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);
374 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);
375 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);
376 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);
377 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);
378 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);
379 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);
380 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);
381 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);
382 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);
383 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);
384 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);
385 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);
386 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);
387 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);
388 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);
389 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);
390 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);
391 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);
392 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);
393 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);
394 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);
395 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);
396 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);
397 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);
398 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);
399 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);
400 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);
401 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);
402 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);
406 for (
int l = 0; l < 3; l++)
408 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));
409 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));
410 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));
411 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));
412 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));
413 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));
415 CeeVLR.at(i).at(j).at(k).at(l) += (oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", i, l) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", i, l))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", k, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", k, j)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", j, k) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", j, k)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YeR", k, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", k, j)))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", l, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", l, i))))/4.;
416 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.);
417 CeeSRR.at(i).at(j).at(k).at(l) += -0.25*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", k, l) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", k, l))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", i, j)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i)))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", l, k) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", l, k))));
419 for (
int p = 0; p < 3; p++)
420 for (
int r = 0; r < 3; r++)
422 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);
423 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);
424 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);
425 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);
426 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);
427 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);
428 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);
430 CedSRL.at(i).at(j).at(k).at(l) +=
VdRd(k,p) * (-0.5*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", i, j))*(
mySMEFT.getSMEFTCoeffEW(
"YdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI", p, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", r, p)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI", p, r)))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i)))))*
VdL(r,l);
431 CedSRR.at(i).at(j).at(k).at(l) +=
VdLd(k,p) * (-0.5*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", i, j))*(
mySMEFT.getSMEFTCoeffEW(
"YdR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI", r, p)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", p, r)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YdR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI", r, p)))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i))))) *
VdR(r,l);
434 for (
int s = 0;
s < 3;
s++)
435 for (
int t = 0;
t < 3;
t++)
437 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);
438 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);
439 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);
440 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);
441 CddVLL.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * (-0.125*(gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t))*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))) - (gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))/8. - (gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t))))/8.) *
VdL(r,j) *
VdL(
t,l);
442 CddVRR.at(i).at(j).at(k).at(l) +=
VdRd(i,p) *
VdRd(k,
s) * ((gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s,
t)))/8.) *
VdR(r,j) *
VdR(
t,l);
443 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * (-0.25*(gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s,
t))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))))) *
VdL(r,j) *
VdR(
t,l);
445 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", r,
s)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
t, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
t, p))))/12.) *
VdL(r,j) *
VdR(
t,l);
446 CddV8LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", r,
s)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
t, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
t, p))))/2.) *
VdL(r,j) *
VdR(
t,l);
447 CddS1RR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * (-0.25*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI",
s,
t))*(
mySMEFT.getSMEFTCoeffEW(
"YdR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI", r, p)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", p, r)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YdR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI", r, p)))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
t,
s))))) *
VdR(r,j) *
VdR(
t,l);
448 CddVLL.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI",
t,
s)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI",
t,
s)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[k])/(8.*sqrt(2))) *
VdL(r,j) *
VdL(
t,l);
449 CddVLL.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI",
s,
t)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI",
s,
t)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[l])/(8.*sqrt(2))) *
VdL(r,j) *
VdR(
t,l);
450 CddVLL.at(i).at(j).at(k).at(l) +=
VdRd(i,p) *
VdLd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))*Md[i])/(8.*sqrt(2))) *
VdL(r,j) *
VdL(
t,l);
451 CddVLL.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))*Md[j])/(8.*sqrt(2))) *
VdR(r,j) *
VdL(
t,l);
452 CddVRR.at(i).at(j).at(k).at(l) +=
VdRd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI",
t,
s)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI",
t,
s)))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*Md[l])/(4.*sqrt(2)))*
VdR(r,j) *
VdL(
t,l);
453 CddVRR.at(i).at(j).at(k).at(l) +=
VdRd(i,p) *
VdLd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI",
s,
t)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI",
s,
t)))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p, r))*Md[k])/(4.*sqrt(2))) *
VdR(r,j) *
VdR(
t,l);
454 CddVRR.at(i).at(j).at(k).at(l) +=
VdRd(i,p) *
VdRd(k,
s) *((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s,
t))*Md[j])/(4.*sqrt(2))) *
VdL(r,j) *
VdR(
t,l);
455 CddVRR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) *((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p ,r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p ,r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p ,r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p ,r)))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s,
t))*Md[i])/(4.*sqrt(2))) *
VdR(r,j) *
VdR(
t,l);
456 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * (-0.041666666666666664*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", p,
t))*(3*
mySMEFT.FlagQuadraticWC*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", r,
s)))))) *
VdL(r,j) *
VdR(
t,l);
457 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI",
t,
s)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI",
t,
s)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[l])/(2.*sqrt(2))) *
VdL(r,j) *
VdR(
t,l);
458 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI",
s,
t)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI",
s,
t)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[k])/(2.*sqrt(2))) *
VdL(r,j) *
VdR(
t,l);
459 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", r, p)))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s,
t))*Md[i])/(2.*sqrt(2))) *
VdL(r,j) *
VdR(
t,l);
460 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CdBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CdWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdWI", p, r)))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s,
t))*Md[j])/(2.*sqrt(2))) *
VdL(r,j) *
VdR(
t,l);
461 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * (-0.041666666666666664*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s, r)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p,
t)))*Md[k]*Md[l])) *
VdL(r,j) *
VdR(
t,l);
462 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s, r))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p,
t)))*Md[j]*Md[l])/24.) *
VdL(r,j) *
VdR(
t,l);
463 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p,
t))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s, r)))*Md[i]*Md[k])/24.) *
VdL(r,j) *
VdR(
t,l);
464 CddV1LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * (-0.041666666666666664*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s, r))*Md[i]*Md[j])) *
VdL(r,j) *
VdR(
t,l);
465 CddV8LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * (-0.25*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", p,
t))*(3*
mySMEFT.FlagQuadraticWC*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", r,
s)))))) *
VdL(r,j) *
VdR(
t,l);
466 CddV8LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * (-0.25*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s, r)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p,
t)))*Md[k]*Md[l])) *
VdL(r,j) *
VdR(
t,l);
467 CddV8LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s, r))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p,
t)))*Md[j]*Md[l])/4.) *
VdL(r,j) *
VdR(
t,l);
468 CddV8LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p,
t))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s, r)))*Md[i]*Md[k])/4.) *
VdL(r,j) *
VdR(
t,l);
469 CddV8LR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdRd(k,
s) * (-0.25*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s, r))*Md[i]*Md[j])) *
VdL(r,j) *
VdR(
t,l);
470 CddS1RR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", p, r))*(3*
mySMEFT.FlagQuadraticWC*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI",
s,
t)))))/8.) *
VdR(r,j) *
VdR(
t,l);
471 CddS1RR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s,
t))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*Md[j]*Md[k])/4.) *
VdR(r,j) *
VdR(
t,l);
472 CddS1RR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * (-0.25*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))*Md[j]*Md[l])) *
VdR(r,j) *
VdR(
t,l);
473 CddS1RR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * (-0.25*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p ,r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p ,r))*(
mySMEFT.getSMEFTCoeffEW(
"CHdR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI",
s,
t))*Md[i]*Md[k])) *
VdR(r,j) *
VdR(
t,l);
474 CddS1RR.at(i).at(j).at(k).at(l) +=
VdLd(i,p) *
VdLd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHdR", p ,r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHdI", p ,r))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))*Md[i]*Md[l])/4.) *
VdR(r,j) *
VdR(
t,l);
479 for (
int l = 0; l < 2; l++)
481 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))));
483 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));
484 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]));
486 for (
int p = 0; p < 3; p++)
487 for (
int r = 0; r < 3; r++)
489 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);
490 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);
491 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);
492 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);
496 for (
int k = 0; k < 2; k++)
498 CnueduVLL.at(i).at(i).at(j).at(k) += VCKMd(j, k) * (-0.5 * (delta_g22oMW2 * g22oMW2));
500 CnueduSRR.at(i).at(i).at(j).at(k) += VCKMd(j,k) * (-0.5*((1 + deltag2bar2oMW4)*g2bar2oMW4*Me[i]*Mu[k]));
502 for (
int p = 0; p < 3; p++)
503 for (
int r = 0; r < 3; r++)
505 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);
506 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);
508 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);
509 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);
510 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);
511 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);
512 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);
513 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);
514 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);
515 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);
518 for (
int l = 0; l < 2; l++)
521 CduV1LR.at(i).at(j).at(k).at(l) += VCKMd(i,l) *
VCKM(k,j) * (-0.08333333333333333*((1 + deltag2bar2oMW4)*g2bar2oMW4*Mu[k]*Mu[l]));
522 CduV8LR.at(i).at(j).at(k).at(l) += VCKMd(i,l) *
VCKM(k,j) * (-0.5*((1 + deltag2bar2oMW4)*g2bar2oMW4*Mu[k]*Mu[l]));
524 for (
int p = 0; p < 3; p++)
525 for (
int r = 0; r < 3; r++)
527 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);
528 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);
529 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);
530 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);
531 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);
532 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);
533 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);
535 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);
536 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);
537 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);
538 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);
539 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);
540 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);
541 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);
542 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);
543 CeuSRL.at(i).at(j).at(k).at(l) +=
VuRd(k,p) * (-0.5*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", r, p))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", i, j)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i)))*(
mySMEFT.getSMEFTCoeffEW(
"YuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI", p, r)))))*
VuL(r,l);
544 CeuSRR.at(i).at(j).at(k).at(l) +=
VuLd(k,p) * (-0.5*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p, r))*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CeHR", i, j) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CeHI", i, j)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YeR", j, i) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YeI", j, i)))*(
mySMEFT.getSMEFTCoeffEW(
"YuR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI", r, p))))) *
VuR(r,l);
546 for (
int s = 0;
s < 3;
s++)
547 for (
int t = 0;
t < 3;
t++)
549 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);
550 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);
555 for (
int p = 0; p < 3; p++)
556 for (
int r = 0; r < 3; r++)
558 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);
559 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);
563 for (
int i = 0; i < 3; i++)
564 for (
int j = 0; j < 2; j++)
566 CnuuVLL.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * (8 * delta_sbar * sbar2 + delta_gZ2oMZ2 * (-3 + 4 * sbar2))) / 12.;
567 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.;
568 CeuVRR.at(i).at(i).at(j).at(j) += (2 * (delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2) / 3.;
569 CnuuVLR.at(i).at(i).at(j).at(j) += ((delta_gZ2oMZ2 + 2 * delta_sbar) * gZ2oMZ2 * sbar2) / 3.;
570 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.;
571 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.;
573 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]);
574 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.;
575 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.;
576 CeuSRL.at(i).at(i).at(j).at(j) += (-2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Me[i]*Mu[j])/3.;
577 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.;
578 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.;
579 CeuSRR.at(i).at(i).at(j).at(j) += (2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Me[i]*Mu[j])/3.;
580 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.;
582 for (
int k = 0; k < 2; k++)
583 for (
int p = 0; p < 3; p++)
584 for (
int r = 0; r < 3; r++)
586 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);
587 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);
588 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);
589 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);
590 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);
591 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);
592 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);
593 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);
594 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);
595 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);
597 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);
598 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);
599 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);
600 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);
601 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);
602 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);
603 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);
604 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);
605 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);
606 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);
607 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);
608 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);
609 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);
610 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);
611 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);
612 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);
613 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);
614 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);
615 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);
616 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);
617 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);
618 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);
619 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);
620 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);
621 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);
622 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);
623 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);
624 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);
625 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);
626 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);
631 for (
int i = 0; i < 2; i++)
632 for (
int j = 0; j < 3; j++)
634 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.;
635 CudV1RR.at(i).at(i).at(j).at(j) += (2 * (delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2) / 9.;
636 CudV1LR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * sbar2 * (delta_gZ2oMZ2 * (-3 + 4 * sbar2) + 2 * delta_sbar * (-3 + 8 * sbar2))) / 18.;
637 CueVLR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * sbar2 * (delta_gZ2oMZ2 * (-3 + 4 * sbar2) + 2 * delta_sbar * (-3 + 8 * sbar2))) / 6.;
639 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.;
640 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]);
641 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]);
642 CudduV1LR.at(i).at(j).at(j).at(i) += ((1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[j]*Mu[i])/27.;
643 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.;
644 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]);
645 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]);
646 CudduV8LR.at(i).at(j).at(j).at(i) += (2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[j]*Mu[i])/9.;
647 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]);
648 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.;
649 CudS1RR.at(i).at(i).at(j).at(j) += (2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Md[j]*Mu[i])/9.;
650 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]);
652 for (
int k = 0; k < 3; k++)
654 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))));
655 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))));
656 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)));
657 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.;
658 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))));
659 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.;
661 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));
662 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.;
663 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.;
664 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));
665 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.;
666 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.;
668 for (
int p = 0; p < 3; p++)
669 for (
int r = 0; r < 3; r++)
671 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);
672 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);
673 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);
674 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);
676 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);
677 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);
678 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);
679 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);
680 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);
681 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);
682 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);
683 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);
684 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);
685 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);
686 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);
687 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);
688 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);
689 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);
692 for (
int l = 0; l < 2; l++)
695 CudduS1RR.at(i).at(j).at(k).at(l) +=
VCKM(i,j) * VCKMd(k,l) * (-0.5*((1 + deltag2bar2oMW4)*g2bar2oMW4*Md[j]*Mu[l]));
697 for (
int p = 0; p < 3; p++)
698 for (
int r = 0; r < 3; r++)
700 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);
702 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);
703 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);
704 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);
705 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);
706 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);
707 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);
709 for (
int s = 0;
s < 3;
s++)
710 for (
int t = 0;
t < 3;
t++)
712 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);
713 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);
715 CudduV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VdRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", r,
s)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)))*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
t, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
t, p))))/12.) *
VdL(r,j) *
VuR(
t,l);
716 CudduV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VdRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI", r,
s)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
s, r)))*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
t, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
t, p))))/2.) *
VdL(r,j) *
VuR(
t,l);
724 for (
int i = 0; i < 2; i++)
725 for (
int j = 0; j < 2; j++)
727 CuuVLL.at(i).at(i).at(j).at(j) += -0.013888888888888888 * (gZ2oMZ2 * (-3 + 4 * sbar2) * (16 * delta_sbar * sbar2 + delta_gZ2oMZ2 * (-3 + 4 * sbar2)));
728 CuuVRR.at(i).at(i).at(j).at(j) += (-2 * (delta_gZ2oMZ2 + 4 * delta_sbar) * gZ2oMZ2 * sbar2 * sbar2) / 9.;
729 CuuV1LR.at(i).at(i).at(j).at(j) += (gZ2oMZ2 * (2 * delta_sbar * (3 - 8 * sbar2) + delta_gZ2oMZ2 * (3 - 4 * sbar2)) * sbar2) / 9.;
731 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]);
732 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.;
733 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.;
734 CuuV1LR.at(i).at(j).at(j).at(i) += (-2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Mu[i]*Mu[j])/27.;
735 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]);
736 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.;
737 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.;
738 CuuV8LR.at(i).at(j).at(j).at(i) += (-4*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Mu[i]*Mu[j])/9.;
739 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.;
740 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]);
741 CuuS1RR.at(i).at(i).at(j).at(j) += (-2*(1 + deltagZbar2oMZ4 + 4*delta_sbar)*gZbar2oMZ4*pow(sbar2,2)*Mu[i]*Mu[j])/9.;
742 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.;
744 for (
int k = 0; k < 3; k++)
746 for (
int l = 0; l < 3; l++)
748 CudV1LL.at(i).at(j).at(k).at(l) +=
VCKM(i, l) * VCKMd(k, j) * (-0.16666666666666666 * (delta_g22oMW2 * g22oMW2));
749 CudV8LL.at(i).at(j).at(k).at(l) +=
VCKM(i, l) * VCKMd(k, j) * (-(delta_g22oMW2 * g22oMW2));
751 CudV1LR.at(i).at(j).at(k).at(l) +=
VCKM(i,l) * VCKMd(k,j) * (-0.08333333333333333*((1 + deltag2bar2oMW4)*g2bar2oMW4*Md[k]*Md[l]));
752 CudV8LR.at(i).at(j).at(k).at(l) +=
VCKM(i,l) * VCKMd(k,j) * (-0.5*((1 + deltag2bar2oMW4)*g2bar2oMW4*Md[k]*Md[l]));
754 for (
int p = 0; p < 3; p++)
755 for (
int r = 0; r < 3; r++)
757 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);
758 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);
759 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);
760 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);
761 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);
763 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);
764 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);
765 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);
766 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);
767 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);
768 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);
769 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);
770 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);
771 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);
772 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);
773 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);
774 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);
775 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);
776 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);
777 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);
778 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);
780 for (
int s = 0;
s < 3;
s++)
781 for (
int t = 0;
t < 3;
t++)
783 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);
784 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);
785 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);
786 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);
787 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);
788 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);
789 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);
790 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);
792 CudS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VdLd(k,
s) * ((oneoMh2*(-3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p, r))*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
t,
s)) + (-3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CdHR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CdHI",
s,
t)) + 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YdR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YdI",
t,
s)))*(
mySMEFT.getSMEFTCoeffEW(
"YuR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI", r, p))))/2.) *
VuR(r,j) *
VdR(
t,l);
798 for (
int k = 0; k < 2; k++)
800 for (
int p = 0; p < 3; p++)
801 for (
int r = 0; r < 3; r++)
803 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);
804 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);
805 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);
806 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);
807 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);
808 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);
810 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);
811 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);
812 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);
813 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);
814 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);
815 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);
816 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);
817 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);
818 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);
819 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);
820 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);
821 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);
822 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);
823 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);
824 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);
825 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);
826 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);
827 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);
828 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);
829 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);
830 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);
831 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);
832 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);
833 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);
836 for (
int l = 0; l < 2; l++)
837 for (
int p = 0; p < 3; p++)
838 for (
int r = 0; r < 3; r++)
839 for (
int s = 0;
s < 3;
s++)
840 for (
int t = 0;
t < 3;
t++)
842 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);
843 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);
844 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);
845 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);
846 CuuVLL.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * ((gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t))*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)) - gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)) + gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r))*(-(
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) + (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t))))/8.) *
VuL(r,j) *
VuL(
t,l);
847 CuuVRR.at(i).at(j).at(k).at(l) +=
VuRd(i,p) *
VuRd(k,
s) * (-0.125*(gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t)))) *
VuR(r,j) *
VuR(
t,l);
848 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * (-0.25*(gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t))) + (gZ2oMZ2*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t)))/4.) *
VuL(r,j) *
VuR(
t,l);
850 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", r,
s)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
s, r)))*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
t, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
t, p))))/12.) *
VuL(r,j) *
VuR(
t,l);
851 CuuV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
s, r)) + (3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", r,
s)) - 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
s, r)))*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
t, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
t, p))))/2.) *
VuL(r,j) *
VuR(
t,l);
852 CuuS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * ((oneoMh2*(-3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI",
s,
t))*(
mySMEFT.getSMEFTCoeffEW(
"YuR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI", r, p)) + (-3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p, r)) + 2*(1 + 2*(
mySMEFT.getSMEFTCoeffEW(
"CHbox") - 0.25 *
mySMEFT.getSMEFTCoeffEW(
"CHD")) *
v2 + deltaoneoMh2)*(
mySMEFT.getSMEFTCoeffEW(
"YuR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI", r, p)))*(
mySMEFT.getSMEFTCoeffEW(
"YuR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"YuI",
t,
s))))/4.) *
VuR(r,j) *
VuR(
t,l);
853 CuuVLL.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI",
t,
s)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI",
t,
s)))*Mu[k])/(4.*sqrt(2))) *
VuL(r,j) *
VuL(
t,l);
854 CuuVLL.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI",
s,
t)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI",
s,
t)))*Mu[l])/(4.*sqrt(2))) *
VuL(r,j) *
VuR(
t,l);
855 CuuVLL.at(i).at(j).at(k).at(l) +=
VuRd(i,p) *
VuLd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[i])/(4.*sqrt(2))) *
VuL(r,j) *
VuL(
t,l);
856 CuuVLL.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[j])/(4.*sqrt(2))) *
VuR(r,j) *
VuL(
t,l);
857 CuuVRR.at(i).at(j).at(k).at(l) +=
VuRd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI",
t,
s)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI",
t,
s)))*Mu[l])/(4.*sqrt(2)))*
VuR(r,j) *
VuL(
t,l);
858 CuuVRR.at(i).at(j).at(k).at(l) +=
VuRd(i,p) *
VuLd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p, r))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI",
s,
t)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI",
s,
t)))*Mu[k])/(4.*sqrt(2))) *
VuR(r,j) *
VuR(
t,l);
859 CuuVRR.at(i).at(j).at(k).at(l) +=
VuRd(i,p) *
VuRd(k,
s) *((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[j])/(4.*sqrt(2))) *
VuL(r,j) *
VuR(
t,l);
860 CuuVRR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) *((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p ,r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p ,r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p ,r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p ,r)))*Mu[i])/(4.*sqrt(2))) *
VuR(r,j) *
VuR(
t,l);
861 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * (-0.041666666666666664*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p,
t))*(3*
mySMEFT.FlagQuadraticWC*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", r,
s)))))) *
VuL(r,j) *
VuR(
t,l);
862 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI",
t,
s)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR",
t,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI",
t,
s)))*Mu[l])/(2.*sqrt(2))) *
VuL(r,j) *
VuR(
t,l);
863 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI",
s,
t)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI",
s,
t)))*Mu[k])/(2.*sqrt(2))) *
VuL(r,j) *
VuR(
t,l);
864 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", r, p)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", r, p) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", r, p)))*Mu[i])/(2.*sqrt(2))) *
VuL(r,j) *
VuR(
t,l);
865 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ3*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t))*(sbar*(
mySMEFT.getSMEFTCoeffEW(
"CuBR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuBI", p, r)) + cbar*(
mySMEFT.getSMEFTCoeffEW(
"CuWR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuWI", p, r)))*Mu[j])/(2.*sqrt(2))) *
VuL(r,j) *
VuR(
t,l);
866 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * (-0.041666666666666664*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s, r)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p,
t)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p,
t)))*Mu[k]*Mu[l])) *
VuL(r,j) *
VuR(
t,l);
867 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p,
t)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p,
t)))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s, r))*Mu[j]*Mu[l])/24.) *
VuL(r,j) *
VuR(
t,l);
868 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s, r)))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p,
t))*Mu[i]*Mu[k])/24.) *
VuL(r,j) *
VuR(
t,l);
869 CuuV1LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * (-0.041666666666666664*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s, r))*Mu[i]*Mu[j])) *
VuL(r,j) *
VuR(
t,l);
870 CuuV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * (-0.25*(oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p,
t))*(3*
mySMEFT.FlagQuadraticWC*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", r,
s) - gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", r,
s)))))) *
VuL(r,j) *
VuR(
t,l);
871 CuuV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * (-0.25*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s, r)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p,
t)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p,
t)))*Mu[k]*Mu[l])) *
VuL(r,j) *
VuR(
t,l);
872 CuuV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p,
t)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p,
t)))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s, r))*Mu[j]*Mu[l])/4.) *
VuL(r,j) *
VuR(
t,l);
873 CuuV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s, r)))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p,
t))*Mu[i]*Mu[k])/4.) *
VuL(r,j) *
VuR(
t,l);
874 CuuV8LR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuRd(k,
s) * (-0.25*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p,
t))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s, r))*Mu[i]*Mu[j])) *
VuL(r,j) *
VuR(
t,l);
875 CuuS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * ((oneoMh2*(3*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI", p, r))*(3*
mySMEFT.FlagQuadraticWC*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CuHR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CuHI",
s,
t)))))/8.) *
VuR(r,j) *
VuR(
t,l);
876 CuuS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t))*Mu[j]*Mu[k])/8.) *
VuR(r,j) *
VuR(
t,l);
877 CuuS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * (-0.125*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I", p, r)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R", p, r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I", p, r)))*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))*Mu[j]*Mu[l])) *
VuR(r,j) *
VuR(
t,l);
878 CuuS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * (-0.125*(gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p ,r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p ,r))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI",
s,
t))*Mu[i]*Mu[k])) *
VuR(r,j) *
VuR(
t,l);
879 CuuS1RR.at(i).at(j).at(k).at(l) +=
VuLd(i,p) *
VuLd(k,
s) * ((gZbar2oMZ4*
mySMEFT.FlagQuadraticWC*
v2*
v2*((
mySMEFT.getSMEFTCoeffEW(
"CHq1R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq1I",
s,
t)) - (
mySMEFT.getSMEFTCoeffEW(
"CHq3R",
s,
t) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHq3I",
s,
t)))*(
mySMEFT.getSMEFTCoeffEW(
"CHuR", p ,r) + gslpp::complex::i() *
mySMEFT.getSMEFTCoeffEW(
"CHuI", p ,r))*Mu[i]*Mu[l])/8.) *
VuR(r,j) *
VuR(
t,l);
884 for (
int p = 0; p < 3; p++)
885 for (
int r = 0; r < 3; r++)
887 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);
888 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 Particle & getLeptons(const QCD::lepton p) const
A get method to retrieve the member object of a lepton.
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 