623 complex_t C7SMeps_0 = (epsilonb - epsilonbp) / (1. + epsilonb * tanb) * tanb *
F7_2(xt);
624 complex_t C7Heps_0 = (-epsilon0p - epsilonb) / (1. + epsilonb * tanb) * tanb *
F7_2(yt);
642 if ((m_A02 == 0.) && (m_H03 == 0.)) {
643 C7Heps2_0 = -epsilon2 * epsilon1p *
pow(tanb, 2.) / (1. + epsilonb * tanb) / (1. + epsilon0 * tanb) *
F7_2(yt);
644 C7Heps2_0 += epsilon2 /
pow(1. + epsilonb * tanb, 2.) * (1. +
pow(tanb, 2.)) / (1. + epsilon0 * tanb) / 72. *
645 ((cos(alpha) + sin(alpha) * tanb) * (-sin(alpha) + epsilonb * cos(alpha)) *
pow(mass_b_muW / mh, 2.) +
646 (sin(alpha) - cos(alpha) * tanb) * (cos(alpha) + epsilonb * sin(alpha)) *
pow(mass_b_muW / mH0, 2.) +
647 (-cos(atan(tanb)) - sin(atan(tanb)) * tanb) * (sin(atan(tanb)) - epsilonb * cos(atan(tanb))) *
pow(mass_b_muW / mA0, 2.));
649 C7Heps2_0 = -epsilon2 * epsilon1p *
pow(tanb, 2.) / (1. + epsilonb * tanb) / (1. + epsilon0 * tanb) *
F7_2(yt);
650 C7Heps2_0 += epsilon2 /
pow(1. + epsilonb * tanb, 2.) * (1. +
pow(tanb, 2.)) / (1. + epsilon0 * tanb) / 72. *
666 auto calculateContribution = [src, mW, tanb, epsilonb](
auto hFunc,
int X,
int X2,
int ie,
int ae,
bool isChargeps) ->
double {
668 double ratio = std::pow(mW / Mch, 2);
670 double factor = isChargeps ? (-epsilonb / (1.0 + epsilonb * tanb) * tanb) : 1.0;
680 for (
int ie = 0; ie < 2; ++ie) {
681 for (
int ae = 0; ae < 6; ++ae) {
683 C7charg_0 += calculateContribution(
h10, 3, 3, ie, ae,
false);
684 C7charg_0 += Mch / mass_b_muW * calculateContribution(
h20, 3, 4, ie, ae,
false);
685 C7_chargeps_0 += Mch / mass_b_muW * calculateContribution(
h20, 3, 4, ie, ae,
true);
691 return scalar_t(C7SMeps_0 + C7Heps_0 + C7Heps2_0 + C7charg_0 + C7_chargeps_0 + C7H_0);
1004 complex_t C8SMeps_0 = (epsb - epsbp) / (1. + epsb * tanb) * tanb *
F8_2(xt);
1005 complex_t C8Heps_0 = (-eps0p - epsb) / (1. + epsb * tanb) * tanb *
F8_2(yt);
1009 if ((m_A02 == 0.) && (m_H03 == 0.)) {
1010 C8Heps2_0 = -eps2 * eps1p *
pow(tanb, 2.) / ((1. + epsb * tanb) * (1. + eps0 * tanb)) *
F8_2(yt);
1011 C8Heps2_0 += eps2 /
pow(1. + epsb * tanb, 2.) * (1. +
pow(tanb, 2.)) / (1. + eps0 * tanb) / 72. *
1012 ((cos(alpha) + sin(alpha) * tanb) * (-sin(alpha) + epsb * cos(alpha)) *
pow(mb_muW / mh, 2.) +
1013 (sin(alpha) - cos(alpha) * tanb) * (cos(alpha) + epsb * sin(alpha)) *
pow(mb_muW / mH0, 2.) +
1014 (-cos(atan(tanb)) - sin(atan(tanb)) * tanb) * (sin(atan(tanb)) - epsb * cos(atan(tanb))) *
pow(mb_muW / mA0, 2.));
1016 C8Heps2_0 = -eps2 * eps1p *
pow(tanb, 2.) / ((1. + epsb * tanb) * (1. + eps0 * tanb)) *
F8_2(yt);
1017 C8Heps2_0 += -3. * eps2 /
pow(1. + epsb * tanb, 2.) * (1. +
pow(tanb, 2.)) / (1. + eps0 * tanb) / 72. *
1030 auto calculateContribution = [src, mW, tanb, epsb](
auto hFunc,
int X,
int X2,
int ie,
int ae,
bool isChargeps) ->
double {
1033 double ratio = std::pow(mW / Mch, 2.0);
1034 double factor = isChargeps ? (-epsb / (1. + epsb * tanb) * tanb) : 1.0;
1040 hFunc(x) * kappa * factor;
1046 for (
int ie = 0; ie < 2; ++ie) {
1048 for (
int ae = 0; ae < 6; ++ae) {
1049 C8charg_0 += calculateContribution(
h50, 3, 3, ie, ae,
false) + Mch / mb_muW * calculateContribution(
h60, 3, 4, ie, ae,
false);
1050 C8_chargeps_0 += Mch / mb_muW * calculateContribution(
h60, 3, 4, ie, ae,
true);
1056 return scalar_t(C8SMeps_0 + C8Heps_0 + C8Heps2_0 + C8charg_0 + C8_chargeps_0 + C8H_0);
1267 for (
int ie = 0; ie < 2; ie++) {
1269 for (
int ae = 0; ae < 6; ae++) {
1274 D91c += std::pow(mW / src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.0) * src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 1}) * src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) *
1275 h31(ratio_MsqU_ae_Mch_ie, log_mu_W_MsqU_ae);
1279 for (
int je = 0; je < 2; je++) {
1281 double factor_abs = 2.0 * std::fabs(ratio_Mch_je_ie);
1282 double factor_f31_f30 = (
f31(std::pow(ratio_Mch_je_ie, 2), ratio_MsqU_ae_Mch_ie) +
1283 4.0 * (
f30(std::pow(ratio_Mch_je_ie, 2), ratio_MsqU_ae_Mch_ie) +
1284 (
f30(std::pow(ratio_Mch_je_ie, 2), ratio_MsqU_ae_Mch_ie * 1.0001) -
1285 f30(std::pow(ratio_Mch_je_ie, 2), ratio_MsqU_ae_Mch_ie * 0.9999)) / 0.0002) * log_mu_W_MsqU_ae);
1286 double factor_f41_f40 = (
f41(std::pow(ratio_Mch_je_ie, 2), ratio_MsqU_ae_Mch_ie) +
1287 4.0 * (
f40(std::pow(ratio_Mch_je_ie, 2), ratio_MsqU_ae_Mch_ie) +
1288 (
f40(std::pow(ratio_Mch_je_ie, 2), ratio_MsqU_ae_Mch_ie * 1.0001) -
1289 f40(std::pow(ratio_Mch_je_ie, 2), ratio_MsqU_ae_Mch_ie * 0.9999)) / 0.0002) * log_mu_W_MsqU_ae);
1291 C91c += src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,je, ae, 1}) * src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) *
1295 for (
int de=0; de<6; de++) {
1296 for (
int ke=0; ke<6; ke++) {
1297 C91f+= src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,de, ke}) *
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,ke, ae}) * (1+
log(
pow(Q_match/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke}),2.)))*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,je, de, 1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) * (
1298 2.*fabs(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})) *
f60(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.))*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 0+1})*src.get_val(
ParameterType::BSM,
"UMIX", {ie+1, 0+1})
1299 -
f50(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.)) * src.get_val(
ParameterType::BSM,
"VMIX", {je+1, 0+1})*src.get_val(
ParameterType::BSM,
"VMIX", {ie+1, 0+1}));
1304 for (
int be=0; be<3; be++) {
1309 B1c1 += src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,je, ae, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) /
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2) *
1310 (0.5 * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {5, ie, be, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {5, je, be, 1}) *
1311 (
f81(
pow(ratio_Mch, 2),
pow(ratio_MsqU, 2),
pow(ratio_Msn, 2)) +
1312 4 * (
f50(
pow(ratio_Mch, 2),
pow(ratio_MsqU, 2),
pow(ratio_Msn, 2)) +
1313 (
f50(
pow(ratio_Mch, 2),
pow(ratio_MsqU, 2) * 1.0001,
pow(ratio_Msn, 2)) -
1314 f50(
pow(ratio_Mch, 2),
pow(ratio_MsqU, 2) * 0.9999,
pow(ratio_Msn, 2))) / 0.0002) *
1317 B1c2 += src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,je, ae, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) /
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2) *
1318 (src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {6, ie, be, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {6, je, be, 1}) * fabs(ratio_Mch) *
1319 (
f91(
pow(ratio_Mch, 2),
pow(ratio_MsqU, 2),
pow(ratio_Msn, 2)) +
1320 4 * (
f60(
pow(ratio_Mch, 2),
pow(ratio_MsqU, 2),
pow(ratio_Msn, 2)) +
1321 (
f60(
pow(ratio_Mch, 2),
pow(ratio_MsqU, 2) * 1.0001,
pow(ratio_Msn, 2)) -
1322 f60(
pow(ratio_Mch, 2),
pow(ratio_MsqU, 2) * 0.9999,
pow(ratio_Msn, 2))) / 0.0002) *
1326 for (
int ce = 0; ce < 6; ce++) {
1327 for (
int fe = 0; fe < 3; fe++) {
1328 for (
int de = 0; de < 6; de++) {
1329 for (
int ke = 0; ke < 6; ke++) {
1333 (1.0 + log_scale_MsqU_ke) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ce, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) *
1340 (1.0 + log_scale_MsqU_ke) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, de, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) *
1348 for (
int je = 0; je < 2; je++) {
1349 double factor_common = mass24_Mch_ie_squared * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,ae, de}) *
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2)*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,de, ce}) *
1350 (1+
log(
pow(Q_match/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de}),2))) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,je, ae, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ce, 2});
1356 src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {5, ie, fe, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {5, je, fe, 1});
1363 src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {6, ie, fe, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {6, je, fe, 1});
1366 C91c += src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ce, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) *
1380 for (
int de = 0; de < 6; de++) {
1385 src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, de, 2}) *
1395 C91c *= -kappa / 8.0;
1397 complex_t B91c = -(B1c1 - B1c2) * kappa * mW * mW / 2.0 /
pow(g2, 2);
1399 C91f *= kappa / 6.0;
1402 complex_t B91f = (B1f1 - B1f2) * 2.0 / 3.0 * kappa /
pow(g2, 2);
1404 complex_t C9four_1 = (1. - 4. * sw2) / sw2 * C91f - B91f / sw2 - D91f;
1406 complex_t C9charg_1=(1.-4.*sw2)/sw2*C91c-B91c/sw2-D91c;
1410 return scalar_t(C9four_1+C9charg_1 + C9H_1);
1601 for (
int ie = 0; ie < 2; ie++) {
1602 for (
int ae = 0; ae < 6; ae++) {
1607 for (
int je = 0; je < 2; je++) {
1609 double factor_abs = 2.0 * std::fabs(ratio_Mch_je_ie);
1610 double factor_f31_f30 = (
f31(std::pow(ratio_Mch_je_ie, 2), ratio_MsqU_ae_Mch_ie) +
1611 4.0 * (
f30(std::pow(ratio_Mch_je_ie, 2), ratio_MsqU_ae_Mch_ie) +
1612 (
f30(std::pow(ratio_Mch_je_ie, 2), ratio_MsqU_ae_Mch_ie * 1.0001) -
1613 f30(std::pow(ratio_Mch_je_ie, 2), ratio_MsqU_ae_Mch_ie * 0.9999)) / 0.0002) * log_mu_W_MsqU_ae);
1614 double factor_f41_f40 = (
f41(std::pow(ratio_Mch_je_ie, 2), ratio_MsqU_ae_Mch_ie) +
1615 4.0 * (
f40(std::pow(ratio_Mch_je_ie, 2), ratio_MsqU_ae_Mch_ie) +
1616 (
f40(std::pow(ratio_Mch_je_ie, 2), ratio_MsqU_ae_Mch_ie * 1.0001) -
1617 f40(std::pow(ratio_Mch_je_ie, 2), ratio_MsqU_ae_Mch_ie * 0.9999)) / 0.0002) * log_mu_W_MsqU_ae);
1619 C91c += src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,je, ae, 1}) * src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) *
1624 for (
int de=0; de<6; de++) {
1625 for (
int ke=0; ke<6; ke++) {
1626 C91f+=src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,de, ke}) *
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,ke, ae}) * (1.+
log(
pow(Q_match/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke}),2.)))*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,je, de, 1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) * (
1627 2.*fabs(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})) *
f60(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.))*src.get_val(
ParameterType::BSM,
"UMIX", {je+1,0+1})*src.get_val(
ParameterType::BSM,
"UMIX", {ie+1,0+1})
1628 -
f50(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.)) * src.get_val(
ParameterType::BSM,
"VMIX", {je+1,0+1})*src.get_val(
ParameterType::BSM,
"VMIX", {ie+1,0+1}));
1630 C91f_b +=src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,de, ke}) *
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,ke, ae}) * (1+
log(
pow(Q_match/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke}),2.)))*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,je, de, 1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) * (
1631 2.*fabs(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})) *
f60(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.))*src.get_val(
ParameterType::BSM,
"UMIX", {je+1,0+1})*src.get_val(
ParameterType::BSM,
"UMIX", {ie+1,0+1})
1632 -
f50(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.)) * src.get_val(
ParameterType::BSM,
"VMIX", {je+1,0+1})*src.get_val(
ParameterType::BSM,
"VMIX", {ie+1,0+1}));
1634 C91f_60 += 2.*fabs(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})) *
f60(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.))*src.get_val(
ParameterType::BSM,
"UMIX", {je+1,0+1})*src.get_val(
ParameterType::BSM,
"UMIX", {ie+1,0+1});
1635 C91f_50 +=
f50(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2.)) * src.get_val(
ParameterType::BSM,
"VMIX", {je+1,0+1})*src.get_val(
ParameterType::BSM,
"VMIX", {ie+1,0+1});
1636 C91f_beg += src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,de, ke}) *
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,ke, ae}) * (1+
log(
pow(Q_match/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke}),2.)))*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,je, de, 1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2});
1638 C91f_beg_1 += (1.+
log(
pow(Q_match/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke}),2.)))*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,je, de, 1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2});
1642 C91f_beg_11 += src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,je, de, 1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2});
1648 for (
int be=0; be<3; be++) {
1653 B1c1 += src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,je, ae, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) /
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2) *
1654 (0.5 * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {5, ie, be, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {5, je, be, 1}) *
1655 (
f81(
pow(ratio_Mch, 2),
pow(ratio_MsqU, 2),
pow(ratio_Msn, 2)) +
1656 4 * (
f50(
pow(ratio_Mch, 2),
pow(ratio_MsqU, 2),
pow(ratio_Msn, 2)) +
1657 (
f50(
pow(ratio_Mch, 2),
pow(ratio_MsqU, 2) * 1.0001,
pow(ratio_Msn, 2)) -
1658 f50(
pow(ratio_Mch, 2),
pow(ratio_MsqU, 2) * 0.9999,
pow(ratio_Msn, 2))) / 0.0002) *
1661 B1c2 += src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,je, ae, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) /
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}), 2) *
1662 (src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {6, ie, be, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {6, je, be, 1}) * fabs(ratio_Mch) *
1663 (
f91(
pow(ratio_Mch, 2),
pow(ratio_MsqU, 2),
pow(ratio_Msn, 2)) +
1664 4 * (
f60(
pow(ratio_Mch, 2),
pow(ratio_MsqU, 2),
pow(ratio_Msn, 2)) +
1665 (
f60(
pow(ratio_Mch, 2),
pow(ratio_MsqU, 2) * 1.0001,
pow(ratio_Msn, 2)) -
1666 f60(
pow(ratio_Mch, 2),
pow(ratio_MsqU, 2) * 0.9999,
pow(ratio_Msn, 2))) / 0.0002) *
1671 for (
int ce = 0; ce < 6; ce++) {
1672 for (
int fe = 0; fe < 3; fe++) {
1673 for (
int de = 0; de < 6; de++) {
1674 for (
int ke = 0; ke < 6; ke++) {
1678 (1.0 + log_scale_MsqU_ke) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ce, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) *
1685 (1.0 + log_scale_MsqU_ke) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, de, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) *
1692 (1.0 + log_scale_MsqU_ke) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ce, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) *
1699 (1.0 + log_scale_MsqU_ke) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, de, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) *
1707 for (
int je = 0; je < 2; je++) {
1708 double factor_common = mass24_Mch_ie_squared * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,ae, de}) *
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2)*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,de, ce}) *
1709 (1+
log(
pow(Q_match/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de}),2))) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,je, ae, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ce, 2});
1715 src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {5, ie, fe, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {5, je, fe, 1});
1722 src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {6, ie, fe, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {6, je, fe, 1});
1728 C91c += src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ce, 1}) * src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3,ie, ae, 2}) *
1764 C91c *= -kappa / 8.0;
1766 complex_t B101c = (B1c1 + B1c2) * kappa * mW * mW / 2.0 /
pow(g2, 2);
1768 C91f *= kappa / 6.0;
1770 complex_t B101f = -(B1f1 + B1f2) * 2.0 / 3.0 * kappa /
pow(g2, 2);
1776 complex_t C10four_1 = (B101f - C91f) / sw2;
1779 double C10H_1 = -
C9llH1(xt,yt,lu,
log(
pow(Q_match/mH,2.)))/sw2;
1785 return C10four_1 + C10charg_1 + C10H_1;
3052 complex_t CQ1H_1=(NQ11H+BQ11H)*mass_b_muW/sw2;
3065 complex_t Dp{0}, Dm{0}, temp{0}, temp2{0};
3066 complex_t a0a{0}, a0b{0}, a0c{0}, a0Q1{0}, a0Q2{0}, a0p{0}, a1{0},a2p{0};
3068 for(
int ie=0;ie<2;ie++) {
3069 for(
int ae=0;ae<6;ae++){
3070 for(
int je=0;je<2;je++) {
3071 for(
int be=0;be<3;be++) {
3072 BQ11c1+=src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3, je, ae, 1})*src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {4, ie, ae, 2})/src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*(src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {6, ie, be, 1})*src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {5, je, be, 1})*(
f121(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {16, be})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))+4.*(
f50(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*1.0001,
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {16, be})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))-
f50(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*0.9999,
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {16, be})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)))/0.0002*log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae}),2.))));
3073 BQ11c2+=src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3, je, ae, 1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {4, ie, ae, 2})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*(src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {5, ie, be, 1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {6, je, be, 1})*fabs(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}))*(
f131(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {16, be})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))+4.*(
f60(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*1.0001,
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {16, be})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))-
f60(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*0.9999,
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {16, be})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)))/0.0002*
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae}),2.))));
3075 for(
int me=0;me<6;me++){
3076 for(
int ne=0;ne<3;ne++) {
3079 for(
int fe=0;fe<3;fe++) {
3080 Dp+=src.get_val(
ParameterType::WILSON,
"WPARAM_MATCH_BSM", {2, fe})/
sqrt(2.)/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*muQ*(src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {2, ae, fe})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, me, fe})+src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, ae, fe})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {2, me, fe}));
3081 Dm+=src.get_val(
ParameterType::WILSON,
"WPARAM_MATCH_BSM", {2, fe})/
sqrt(2.)/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*muQ*(src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {2, ae, fe})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, me, fe})-src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, ae, fe})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {2, me, fe}));
3083 a0a=-(fabs(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}))*(
f181(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.))+4.*(
f30(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.)*1.0001)-
f30(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.)*0.9999))/0.0002*
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae}),2.)))*src.get_val(
ParameterType::BSM,
"UMIX", {ie+1, 1+1})*src.get_val(
ParameterType::BSM,
"VMIX", {je+1, 0+1}))*
kron(ae,me);
3084 a0b=-((
f191(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.))+4.*(
f40(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.)*1.0001)-
f40(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.)*0.9999))/0.0002*
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae}),2.)))*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})*src.get_val(
ParameterType::BSM,
"UMIX", {ie+1, 0+1}))*
kron(ae,me);
3085 a0c=1./mW*(
f171(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, me})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))+4.*(
f30(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, me})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))+(
f30(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*1.0001,
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, me})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))-
f30(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*0.9999,
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, me})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)))/0.0002+(
f30(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, me})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*1.0001)-
f30(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, me})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*0.9999))/0.0002)*
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae}),2.)))*
kron(ie,je);
3087 a0Q1=a0a+a0b+Dp*a0c;
3088 a0Q2=-a0a+a0b+Dm*a0c;
3089 a0p=4.*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {12,ae, ie, be, ne})/mW/(src.get_val(
ParameterType::SM,
"VCKM", {be, 2})*src.get_val(
ParameterType::SM,
"VCKM", {ne, 1})/src.get_val(
ParameterType::SM,
"VCKM", {2, 2})/src.get_val(
ParameterType::SM,
"VCKM", {2, 1}))/src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})*
f151(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))*
kron(ie,je)*
kron(ae,me)*
kron(be,ne);
3090 a1=src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/
sqrt(2.)/mW*(
f161(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))+4.*(
f80(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*1.0001)-
f80(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*0.9999))/0.0002*
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae}),2.)))*
kron(ie,je)*
kron(ae,me);
3091 a2p=src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, me, be})*(src.get_val(
ParameterType::SM,
"VCKM", {be, 2})*src.get_val(
ParameterType::SM,
"VCKM", {ne, 1})/src.get_val(
ParameterType::SM,
"VCKM", {2, 2})/src.get_val(
ParameterType::SM,
"VCKM", {2, 1}))*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})/2./mW*
f151(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))*
kron(ie,je)*
kron(ae,me)*
kron(be,ne);
3093 NQ11c+=src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {12,ae, ie, be, ne})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, me, be})*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})*(a0Q1+a1*tanb)
3095 +src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, me, be})*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})*a2p*
pow(src.get_val(
ParameterType::SM,
"MASS", 3)*tanb,2.);
3096 NQ21c+=src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {12,ae, ie, be, ne})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, me, be})*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})*(a0Q2+a1*tanb)
3098 +src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, me, be})*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})*a2p*
pow(src.get_val(
ParameterType::SM,
"MASS", 3)*tanb,2.);
3104 for(
int be=0;be<6;be++) {
3105 for(
int ce=0;ce<6;ce++) {
3106 for(
int fe=0;fe<3;fe++) {
3107 BQ11f1+=-src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3, je, be, 1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {4, ie, ae, 2})*
pow(mW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9, ae, ce})*src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ce})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9, ce, be})*(1.+
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ce}),2.))) *(
f90(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, be})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {16, fe})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {6, ie, fe, 2})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {5, je, fe, 2}));
3108 BQ11f2+=-src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3, je, be, 1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {4, ie, ae, 2})*
pow(mW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9, ae, ce})*src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ce})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9, ce, be})*(1.+
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ce}),2.))) *(fabs(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}))*
f100(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, be})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {16, fe})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {5, ie, fe, 1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {6, je, fe, 1}));
3117 for(
int me=0;me<3;me++) {
3118 for(
int ne=0;ne<3;ne++) {
3119 for(
int de=0;de<6;de++) {
3120 for(
int ke=0;ke<6;ke++) {
3122 temp2 = src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {12,ae, ie, me, ne})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, de, me})*src.get_val(
ParameterType::BSM,
"UMIX", {ie+1, 1+1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9, ae, ke})*src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9, ke, de})*(1.+
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke}),2.)))*tanb*src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/
sqrt(2.)*
f30(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.));
3125 for(
int ce=0;ce<6;ce++) {
3128 for(
int fe=0;fe<3;fe++)
3130 Dp+=src.get_val(
ParameterType::WILSON,
"WPARAM_MATCH_BSM", {2, fe})/
sqrt(2.)/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*muQ*(src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {2, de, fe})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, ce, fe})+src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, de, fe})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {2, ce, fe}));
3131 Dm+=src.get_val(
ParameterType::WILSON,
"WPARAM_MATCH_BSM", {2, fe})/
sqrt(2.)/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*muQ*(src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {2, de, fe})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, ce, fe})-src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, de, fe})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {2, ce, fe}));
3133 temp=src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {12,ae, ie, me, ne})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, ce, me})*src.get_val(
ParameterType::BSM,
"UMIX", {ie+1, 1+1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9, ae, ke})*src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9, ke, de})*
3134 (1.+
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke}),2.)))*
f60(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ce})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.));
3138 for(
int je=0;je<2;je++) {
3139 temp=-src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {12,ae, ie, me, ne})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, de, me})*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9, ae, ke})*src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9, ke, de})*
3140 (1.+
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke}),2.)))*mW*(fabs(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}))*
f60(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.))*src.get_val(
ParameterType::BSM,
"UMIX", {ie+1, 1+1})*src.get_val(
ParameterType::BSM,
"VMIX", {je+1, 0+1})+
3141 f50(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.))*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})*src.get_val(
ParameterType::BSM,
"VMIX", {ie+1, 0+1}));
3156 complex_t BQ11c=(BQ11c1+BQ11c2)*kappa*mW*mW/2./g2/g2/sw2;
3160 NQ11c*=src.get_val(
ParameterType::WILSON,
"WPARAM_SI_SM", lepton_mass_slot)*(tanb)*tanb/mW/(mH*mH-mW*mW)*aY*mass_b_muW/sw2;
3161 NQ21c*=-src.get_val(
ParameterType::WILSON,
"WPARAM_SI_SM", lepton_mass_slot)*(tanb)*tanb/mW/(mH*mH-mW*mW)*aY*mass_b_muW/sw2;
3170 complex_t BQ11f=(BQ11f1+BQ11f2)*2./3.*kappa/g2/g2/sw2;
3174 NQ11f*=-4./3.*src.get_val(
ParameterType::WILSON,
"WPARAM_SI_SM", lepton_mass_slot)*(tanb)*tanb/mW/mW/(mH*mH-mW*mW)*aY*mass_b_muW/sw2;
3176 NQ21f*=4./3.*src.get_val(
ParameterType::WILSON,
"WPARAM_SI_SM", lepton_mass_slot)*(tanb)*tanb/mW/mW/(mH*mH-mW*mW)*aY*mass_b_muW/sw2;
3180 complex_t coeff_temp = (CQ1H_1+CQ1charg_1)/epsfac+CQ1four_1;
3578 complex_t CQ1H_1=(NQ11H+BQ11H)*mass_b_muW/sw2;
3592 complex_t Dp{0}, Dm{0}, temp{0}, temp2{0};
3593 complex_t a0a{0}, a0b{0}, a0c{0}, a0Q1{0}, a0Q2{0}, a0p{0}, a1{0},a2p{0};
3595 for(
int ie=0;ie<2;ie++) {
3596 for(
int ae=0;ae<6;ae++){
3597 for(
int je=0;je<2;je++) {
3598 for(
int be=0;be<3;be++) {
3599 BQ11c1+=src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3, je, ae, 1})*src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {4, ie, ae, 2})/src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*(src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {6, ie, be, 1})*src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {5, je, be, 1})*(
f121(
pow(src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {16, be})/src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))+4.*(
f50(
pow(src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*1.0001,
pow(src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {16, be})/src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))-
f50(
pow(src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*0.9999,
pow(src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {16, be})/src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)))/0.0002*log(
pow(mass_top_muW/src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae}),2.))));
3600 BQ11c2+=src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3, je, ae, 1})*src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {4, ie, ae, 2})/src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.
get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*(src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {5, ie, be, 1})*src.
get_val(
ParameterType::WILSON,
"MATRIX_BSM", {6, je, be, 1})*fabs(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}))*(
f131(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {16, be})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))+4.*(
f60(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*1.0001,
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {16, be})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))-
f60(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*0.9999,
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {16, be})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)))/0.0002*
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae}),2.))));
3602 for(
int me=0;me<6;me++){
3603 for(
int ne=0;ne<3;ne++) {
3606 for(
int fe=1;fe<3;fe++) {
3607 Dp+=src.get_val(
ParameterType::WILSON,
"WPARAM_MATCH_BSM", {2, fe})/
sqrt(2.)/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*muQ*(src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {2, ae, fe})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, me, fe})+src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, ae, fe})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {2, me, fe}));
3608 Dm+=src.get_val(
ParameterType::WILSON,
"WPARAM_MATCH_BSM", {2, fe})/
sqrt(2.)/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*muQ*(src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {2, ae, fe})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, me, fe})-src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, ae, fe})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {2, me, fe}));
3610 a0a=-(fabs(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}))*(
f181(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.))+4.*(
f30(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.)*1.0001)-
f30(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.)*0.9999))/0.0002*
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae}),2.)))*src.get_val(
ParameterType::BSM,
"UMIX", {ie+1, 1+1})*src.get_val(
ParameterType::BSM,
"VMIX", {je+1, 0+1}))*
kron(ae,me);
3611 a0b=-((
f191(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.))+4.*(
f40(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.)*1.0001)-
f40(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.)*0.9999))/0.0002*
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae}),2.)))*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})*src.get_val(
ParameterType::BSM,
"UMIX", {ie+1, 0+1}))*
kron(ae,me);
3612 a0c=1./mW*(
f171(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, me})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))+4.*(
f30(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, me})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))+(
f30(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*1.0001,
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, me})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))-
f30(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*0.9999,
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, me})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)))/0.0002+(
f30(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, me})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*1.0001)-
f30(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, me})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*0.9999))/0.0002)*
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae}),2.)))*
kron(ie,je);
3614 a0Q1=a0a+a0b+Dp*a0c;
3615 a0Q2=-a0a+a0b+Dm*a0c;
3616 a0p=4.*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {12,ae, ie, be, ne})/mW/(src.get_val(
ParameterType::SM,
"VCKM", {be, 2})*src.get_val(
ParameterType::SM,
"VCKM", {ne, 1})/src.get_val(
ParameterType::SM,
"VCKM", {2, 2})/src.get_val(
ParameterType::SM,
"VCKM", {2, 1}))/src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})*
f151(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))*
kron(ie,je)*
kron(ae,me)*
kron(be,ne);
3617 a1=src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/
sqrt(2.)/mW*(
f161(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))+4.*(
f80(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*1.0001)-
f80(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*0.9999))/0.0002*
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae}),2.)))*
kron(ie,je)*
kron(ae,me);
3618 a2p=src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, me, be})*(src.get_val(
ParameterType::SM,
"VCKM", {be, 2})*src.get_val(
ParameterType::SM,
"VCKM", {ne, 1})/src.get_val(
ParameterType::SM,
"VCKM", {2, 2})/src.get_val(
ParameterType::SM,
"VCKM", {2, 1}))*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})/2./mW*
f151(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))*
kron(ie,je)*
kron(ae,me)*
kron(be,ne);
3620 NQ11c+=src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {12,ae, ie, be, ne})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, me, be})*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})*(a0Q1+a1*tanb)
3622 +src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, me, be})*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})*a2p*
pow(src.get_val(
ParameterType::SM,
"MASS", 3)*tanb,2.);
3623 NQ21c+=src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {12,ae, ie, be, ne})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, me, be})*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})*(a0Q2+a1*tanb)
3625 +src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, me, be})*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})*a2p*
pow(src.get_val(
ParameterType::SM,
"MASS", 3)*tanb,2.);
3630 for(
int be=0;be<6;be++) {
3631 for(
int ce=0;ce<6;ce++) {
3632 for(
int fe=0;fe<3;fe++) {
3633 BQ11f1+=-src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3, je, be, 1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {4, ie, ae, 2})*
pow(mW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,ae, ce})*src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ce})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9, ce, be})*(1.+
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ce}),2.))) *(
f90(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, be})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {16, fe})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {6,ie, fe, 2})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {5, je, fe, 2}));
3634 BQ11f2+=-src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {3, je, be, 1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {4, ie, ae, 2})*
pow(mW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.)*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,ae, ce})*src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ce})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9, ce, be})*(1.+
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ce}),2.))) *(fabs(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}))*
f100(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, be})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {16, fe})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.))*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {5, ie, fe, 1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {6, je, fe, 1}));
3643 for(
int me=0;me<3;me++) {
3644 for(
int ne=0;ne<3;ne++) {
3645 for(
int de=0;de<6;de++) {
3646 for(
int ke=0;ke<6;ke++) {
3648 temp2 = src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {12,ae, ie, me, ne})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, de, me})*src.get_val(
ParameterType::BSM,
"UMIX", {ie+1, 1+1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,ae, ke})*src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9, ke, de})*(1.+
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke}),2.)))*tanb*src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/
sqrt(2.)*
f30(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.));
3651 for(
int ce=0;ce<6;ce++) {
3654 for(
int fe=0;fe<3;fe++)
3656 Dp+=src.get_val(
ParameterType::WILSON,
"WPARAM_MATCH_BSM", {2, fe})/
sqrt(2.)/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*muQ*(src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {2, de, fe})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, ce, fe})+src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, de, fe})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {2, ce, fe}));
3657 Dm+=src.get_val(
ParameterType::WILSON,
"WPARAM_MATCH_BSM", {2, fe})/
sqrt(2.)/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*muQ*(src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {2, de, fe})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, ce, fe})-src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, de, fe})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {2, ce, fe}));
3659 temp=src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {12,ae, ie, me, ne})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, ce, me})*src.get_val(
ParameterType::BSM,
"UMIX", {ie+1, 1+1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,ae, ke})*src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9, ke, de})*
3660 (1.+
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke}),2.)))*
f60(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ce})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie}),2.));
3664 for(
int je=0;je<2;je++) {
3665 temp=-src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {12,ae, ie, me, ne})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {1, de, me})*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9,ae, ke})*src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je})*src.get_val(
ParameterType::WILSON,
"MATRIX_BSM", {9, ke, de})*
3666 (1.+
log(
pow(mass_top_muW/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ke}),2.)))*mW*(fabs(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}))*
f60(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.))*src.get_val(
ParameterType::BSM,
"UMIX", {ie+1, 1+1})*src.get_val(
ParameterType::BSM,
"VMIX", {je+1, 0+1})+
3667 f50(
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, ie})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, ae})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.),
pow(src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {14, de})/src.get_val(
ParameterType::WILSON,
"WPARAM_SI_BSM", {13, je}),2.))*src.get_val(
ParameterType::BSM,
"UMIX", {je+1, 1+1})*src.get_val(
ParameterType::BSM,
"VMIX", {ie+1, 0+1}));
3684 complex_t BQ21c=-(BQ11c1-BQ11c2)*kappa*mW*mW/2./g2/g2/sw2;
3687 NQ11c*=src.get_val(
ParameterType::WILSON,
"WPARAM_SI_SM", lepton_mass_slot)*(tanb)*tanb/mW/(mH*mH-mW*mW)*aY*mass_b_muW/sw2;
3688 NQ21c*=-src.get_val(
ParameterType::WILSON,
"WPARAM_SI_SM", lepton_mass_slot)*(tanb)*tanb/mW/(mH*mH-mW*mW)*aY*mass_b_muW/sw2;
3693 complex_t BQ21f=-(BQ11f1-BQ11f2)*2./3.*kappa/g2/g2/sw2;
3696 NQ11f*=-4./3.*src.get_val(
ParameterType::WILSON,
"WPARAM_SI_SM", lepton_mass_slot)*(tanb)*tanb/mW/mW/(mH*mH-mW*mW)*aY*mass_b_muW/sw2;
3698 NQ21f*=4./3.*src.get_val(
ParameterType::WILSON,
"WPARAM_SI_SM", lepton_mass_slot)*(tanb)*tanb/mW/mW/(mH*mH-mW*mW)*aY*mass_b_muW/sw2;
3702 complex_t coeff_temp = (CQ2H_1+CQ2charg_1)/epsfac+CQ2four_1;