1881 gslpp::matrix<gslpp::complex> myR(6, 6, 0.);
1882 gslpp::vector<double> myM2Squarks(6, 0.);
1888 if (Dmixingflag == 0) {
1892 else if (Dmixingflag == 1) {
1901 gslpp::complex CLO(0., 0.,
false);
1902 gslpp::vector<gslpp::complex> VCLO(8, 0.);
1909 for (O = 1; O < 9; O++) {
1911 CLO.assign(0., 0.,0);
1914 for (h = 0; h < 6; h++) {
1915 for (k = 0; k < 6; k++) {
1918 myR(h, b) * myR(k, b)
1919 * myR(h, q).conjugate() * myR(k, q).conjugate() *
1921 Dk(myM2Squarks(h), myM2Squarks(k), M2g, M2g, 0) +
1922 11. / 9. *
Dk(myM2Squarks(h), myM2Squarks(k), M2g, M2g, 2));
1927 for (h = 0; h < 6; h++) {
1928 for (k = 0; k < 6; k++) {
1930 CLO += -
Als *
Als * 17. / 18. * M2g * myR(h, b)
1931 * myR(k, b) * myR(h, q + 3).conjugate()
1932 * myR(k, q + 3).conjugate() *
1933 Dk(myM2Squarks(h), myM2Squarks(k), M2g, M2g, 0);
1938 for (h = 0; h < 6; h++) {
1939 for (k = 0; k < 6; k++) {
1941 CLO +=
Als *
Als * 1. / 6. * M2g * myR(h, b)
1942 * myR(k, b) * myR(h, q + 3).conjugate()
1943 * myR(k, q + 3).conjugate() *
1944 Dk(myM2Squarks(h), myM2Squarks(k), M2g, M2g, 0);
1949 for (h = 0; h < 6; h++) {
1950 for (k = 0; k < 6; k++) {
1952 CLO += -
Als *
Als * 7. / 3. * M2g * myR(h, b)
1953 * myR(k, b + 3) * myR(h, q).conjugate()
1954 * myR(k, q + 3).conjugate() *
1955 Dk(myM2Squarks(h), myM2Squarks(k), M2g, M2g, 0) +
1957 Dk(myM2Squarks(h), myM2Squarks(k), M2g, M2g, 2)
1958 * myR(h, b) * myR(k, b + 3) *
1959 (6. * myR(h, q).conjugate() * myR(k, q + 3).conjugate()
1960 + 11. * myR(k, q).conjugate() * myR(h, q + 3).conjugate())
1966 for (h = 0; h < 6; h++) {
1967 for (k = 0; k < 6; k++) {
1969 CLO += -
Als *
Als * 1. / 9. * M2g * myR(h, b)
1970 * myR(k, b + 3) * myR(h, q).conjugate()
1971 * myR(k, q + 3).conjugate() *
1972 Dk(myM2Squarks(h), myM2Squarks(k), M2g, M2g, 0) +
1974 Dk(myM2Squarks(h), myM2Squarks(k), M2g, M2g, 2)
1975 * myR(h, b) * myR(k, b + 3) *
1976 (3. * myR(k, q).conjugate() * myR(h, q + 3).conjugate()
1977 - 2. * myR(h, q).conjugate() * myR(k, q + 3).conjugate())
1983 for (h = 0; h < 6; h++) {
1984 for (k = 0; k < 6; k++) {
1986 CLO += -
Als *
Als * myR(h, b + 3) * myR(k, b + 3)
1987 * myR(h, q + 3).conjugate() * myR(k, q + 3).conjugate() *
1988 (1. / 9. * M2g *
Dk(myM2Squarks(h), myM2Squarks(k), M2g, M2g, 0) -
1989 11. / 9. *
Dk(myM2Squarks(h), myM2Squarks(k), M2g, M2g, 2));
1994 for (h = 0; h < 6; h++) {
1995 for (k = 0; k < 6; k++) {
1997 CLO += -
Als *
Als * 17. / 18. * M2g * myR(h, b + 3)
1998 * myR(k, b + 3) * myR(h, q).conjugate()
1999 * myR(k, q).conjugate() *
2000 Dk(myM2Squarks(h), myM2Squarks(k), M2g, M2g, 0);
2005 for (h = 0; h < 6; h++) {
2006 for (k = 0; k < 6; k++) {
2008 CLO +=
Als *
Als * 1. / 6. * M2g * myR(h, b + 3)
2009 * myR(k, b + 3) * myR(h, q).conjugate()
2010 * myR(k, q).conjugate() *
2011 Dk(myM2Squarks(h), myM2Squarks(k), M2g, M2g, 0);
2017 VCLO.assign(O - 1, CLO);
2350 gslpp::matrix<gslpp::complex> myR(6, 6, 0.);
2351 gslpp::vector<double> myM2Squarks(6, 0.);
2352 gslpp::vector<gslpp::complex> VCLO(8, 0.);
2353 gslpp::complex CLO(0., 0.,
false);
2362 if (Dmixingflag == 0) {
2366 else if (Dmixingflag == 1) {
2374 int D = Dmixingflag;
2377 for (O = 1; O < 9; O++) {
2379 CLO.assign(0., 0., 0);
2382 for (i = 0; i < 4; i++) {
2383 for (h = 0; h < 6; h++) {
2384 for (k = 0; k < 6; k++) {
2386 CLO += -
Als * 2. / 3. / 4. / M_PI *
Dk(myM2Squarks(k), myM2Squarks(h),
2387 MChi0(i) *
MChi0(i), M2g, 2) * myR(k, q).conjugate()
2388 * myR(h, b) *
VdDNL(q, h, i, 1, D).conjugate() *
2391 , M2g, 0) * (myR(h, q).conjugate() * myR(k, q).conjugate()
2392 *
VdDNL(b, h, i, 1, D) *
VdDNL(b, k, i, 1, D)
2393 + myR(h, b) * myR(k, b) *
VdDNL(q, h, i, 1, D).conjugate() *
2394 VdDNL(q, k, i, 1, D).conjugate());
2401 for (i = 0; i < 4; i++) {
2402 for (h = 0; h < 6; h++) {
2403 for (k = 0; k < 6; k++) {
2406 Dk(myM2Squarks(k), myM2Squarks(h),
MChi0(i)
2407 *
MChi0(i), M2g, 0) * (3. * myR(h, b) *
2408 myR(k, q + 3).conjugate() *
VdDNL(b, k, i, 1, D) *
2409 VdDNR(q, h, i, 1, D).conjugate()
2410 + myR(k, b) * myR(h, b) *
2411 VdDNR(q, k, i, 1, D).conjugate() *
VdDNR(q, h, i, 1, D).conjugate() +
2412 myR(k, q + 3).conjugate() * myR(h, q + 3).conjugate() *
2413 VdDNL(b, k, i, 1, D) *
VdDNL(b, h, i, 1, D));
2420 for (i = 0; i < 4; i++) {
2421 for (h = 0; h < 6; h++) {
2422 for (k = 0; k < 6; k++) {
2424 CLO += -
Als / 4. / M_PI *
MChi0(i) *
Mg / 3. *
2425 Dk(myM2Squarks(k), myM2Squarks(h),
MChi0(i)
2426 *
MChi0(i), M2g, 0) * (
2427 myR(h, b) * myR(k, q + 3).conjugate() *
VdDNL(b, k, i, 1, D) *
2428 VdDNR(q, h, i, 1, D).conjugate() - myR(k, b) *
2429 myR(h, b) *
VdDNR(q, k, i, 1, D).conjugate() *
2430 VdDNR(q, h, i, 1, D).conjugate() - myR(k, q + 3).conjugate() *
2431 myR(h, q + 3).conjugate() *
VdDNL(b, k, i, 1, D) *
2432 VdDNL(b, h, i, 1, D));
2439 for (i = 0; i < 4; i++) {
2440 for (h = 0; h < 6; h++) {
2441 for (k = 0; k < 6; k++) {
2443 CLO += -
Als / 4. / M_PI * 2. / 3. *
2444 Dk(myM2Squarks(k), myM2Squarks(h),
MChi0(i)
2445 *
MChi0(i), M2g, 2) * (
2446 myR(h, b) * myR(k, q).conjugate() *
VdDNR(b, k, i, 1, D) *
2447 VdDNR(q, h, i, 1, D).conjugate() +
2449 myR(k, q + 3).conjugate() *
VdDNL(b, k, i, 1, D) *
2450 VdDNL(q, h, i, 1, D).conjugate()
2451 - myR(h, b) * myR(k, b + 3) *
2452 VdDNL(q, k, i, 1, D).conjugate() *
VdDNR(q, h, i, 1, D).conjugate()
2453 - myR(h, q).conjugate() * myR(k, q + 3).conjugate() *
2455 - 3. * myR(k, b + 3) * myR(h, b) *
2456 VdDNL(q, h, i, 1, D).conjugate() *
VdDNR(q, k, i, 1, D).conjugate()
2457 - 3. * myR(k, q + 3).conjugate() *
2458 myR(h, q).conjugate() *
VdDNL(b, h, i, 1, D) *
VdDNR(b, k, i, 1, D)) +
2460 Dk(myM2Squarks(k), myM2Squarks(h),
MChi0(i)*
MChi0(i), M2g, 0) *
2461 (myR(h, b) * myR(k, q + 3).conjugate() *
2462 VdDNR(b, k, i, 1, D) *
VdDNL(q, h, i, 1, D).conjugate() +
2463 myR(h, b + 3) * myR(k, q).conjugate() *
2464 VdDNL(b, k, i, 1, D) *
VdDNR(q, h, i, 1, D).conjugate());
2471 for (i = 0; i < 4; i++) {
2472 for (h = 0; h < 6; h++) {
2473 for (k = 0; k < 6; k++) {
2475 CLO +=
Als / 4. / M_PI * 2. / 3. *
2476 Dk(myM2Squarks(k), myM2Squarks(h),
MChi0(i)
2477 *
MChi0(i), M2g, 2)*(
2478 3. * myR(h, b) * myR(k, q).conjugate() *
VdDNR(b, k, i, 1, D) *
2479 VdDNR(q, h, i, 1, D).conjugate() + 3. * myR(h, b + 3) *
2480 myR(k, q + 3).conjugate() *
VdDNL(b, k, i, 1, D) *
2481 VdDNL(q, h, i, 1, D).conjugate()
2482 - 3. * myR(h, b) * myR(k, b + 3) *
2483 VdDNL(q, k, i, 1, D).conjugate() *
VdDNR(q, h, i, 1, D).conjugate()
2484 - 3. * myR(h, q).conjugate() * myR(k, q + 3).conjugate() *
2485 VdDNL(b, k, i, 1, D) *
VdDNR(b, h, i, 1, D) - myR(k, b + 3) *
2486 myR(h, b) *
VdDNL(q, h, i, 1, D).conjugate() *
2487 VdDNR(q, k, i, 1, D).conjugate() - myR(k, q + 3).conjugate() *
2488 myR(h, q).conjugate() *
VdDNL(b, h, i, 1, D) *
VdDNR(b, k, i, 1, D)) -
2490 Dk(myM2Squarks(k), myM2Squarks(h),
MChi0(i)
2491 *
MChi0(i), M2g, 0)* (myR(h, b) * myR(k, q + 3).conjugate() *
2492 VdDNR(b, k, i, 1, D) *
VdDNL(q, h, i, 1, D).conjugate() +
2493 myR(h, b + 3) * myR(k, q).conjugate() *
2494 VdDNL(b, k, i, 1, D) *
VdDNR(q, h, i, 1, D).conjugate());
2501 for (i = 0; i < 4; i++) {
2502 for (h = 0; h < 6; h++) {
2503 for (k = 0; k < 6; k++) {
2505 CLO += -
Als * 2. / 3. / 4. / M_PI *
Dk(myM2Squarks(k), myM2Squarks(h),
2506 MChi0(i) *
MChi0(i), M2g, 2) * myR(k, q + 3).conjugate()
2507 * myR(h, b + 3) *
VdDNR(q, h, i, 1, D).conjugate() *
2510 , M2g, 0) * (myR(h, q + 3).conjugate() * myR(k, q + 3).conjugate()
2511 *
VdDNR(b, h, i, 1, D) *
VdDNR(b, k, i, 1, D)
2512 + myR(h, b + 3) * myR(k, b + 3) *
VdDNR(q, h, i, 1, D).conjugate() *
2513 VdDNR(q, k, i, 1, D).conjugate());
2520 for (i = 0; i < 4; i++) {
2521 for (h = 0; h < 6; h++) {
2522 for (k = 0; k < 6; k++) {
2525 Dk(myM2Squarks(k), myM2Squarks(h),
MChi0(i)
2526 *
MChi0(i), M2g, 0) * (3. * myR(h, b + 3) *
2527 myR(k, q).conjugate() *
VdDNR(b, k, i, 1, D) *
2528 VdDNL(q, h, i, 1, D).conjugate()
2529 + myR(k, b + 3) * myR(h, b + 3) *
2530 VdDNL(q, k, i, 1, D).conjugate() *
VdDNL(q, h, i, 1, D).conjugate() +
2531 myR(k, q).conjugate() * myR(h, q).conjugate() *
2532 VdDNR(b, k, i, 1, D) *
VdDNR(b, h, i, 1, D));
2539 for (i = 0; i < 4; i++) {
2540 for (h = 0; h < 6; h++) {
2541 for (k = 0; k < 6; k++) {
2543 CLO += -
Als / 4. / M_PI *
MChi0(i) *
Mg / 3. *
2544 Dk(myM2Squarks(k), myM2Squarks(h),
MChi0(i)
2545 *
MChi0(i), M2g, 0) * (
2546 myR(h, b + 3) * myR(k, q).conjugate() *
VdDNR(b, k, i, 1, D) *
2547 VdDNL(q, h, i, 1, D).conjugate() - myR(k, b + 3) *
2548 myR(h, b + 3) *
VdDNL(q, k, i, 1, D).conjugate() *
2549 VdDNL(q, h, i, 1, D).conjugate() - myR(k, q).conjugate() *
2550 myR(h, q).conjugate() *
VdDNR(b, k, i, 1, D) *
2551 VdDNR(b, h, i, 1, D));
2558 VCLO.assign(O - 1, CLO);
4057 double piconst = 1.0/(32.0 * pi * pi);
4058 double sw2 =
mySUSY.StandardModel::sW2(MW);
4059 double stw = sqrt(sw2);
4060 double ctw = sqrt(1.0 - sw2);
4061 double ttw = stw/ctw;
4066 double cdenc = sqrt(2.0)*MW*
cosb;
4067 double cdenn = MW*
cosb;
4069 double g2t = g2/sqrt(2.0);
4071 gslpp::vector<gslpp::complex>
AFunctions(4, 0.);
4074 for (
int a=0;a<4;a++) {
4075 for (
int x=0;x<6;x++) {
4077 NRlE.assign(a, x, - (g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 0) + (mE/cdenn)*
ON(a, 2)*
myRl(x, 3)));
4078 NRlMU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 1) + (mMU/cdenn)*
ON(a, 2)*
myRl(x, 4)));
4079 NRlTAU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 2) + (mTAU/cdenn)*
ON(a, 2)*
myRl(x, 5)));
4081 NLlE.assign(a, x, -(g2t)*((mE/cdenn)*
ON(a, 2)*
myRl(x, 0) + 2.0*
ON(a, 0)*ttw*
myRl(x, 3)));
4082 NLlMU.assign(a, x, -(g2t)*((mMU/cdenn)*
ON(a, 2)*
myRl(x, 1) + 2.0*
ON(a, 0)*ttw*
myRl(x, 4)));
4083 NLlTAU.assign(a, x, -(g2t)*((mTAU/cdenn)*
ON(a, 2)*
myRl(x, 2) + 2.0*
ON(a, 0)*ttw*
myRl(x, 5)));
4092 for (
int a=0;a<2;a++) {
4093 for (
int x=0;x<3;x++) {
4099 CLlE.assign(a, x, g2*mE/cdenc*
myU(a, 1).conjugate()*
myRn(x, 0));
4100 CLlMU.assign(a, x, g2*mMU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 1));
4101 CLlTAU.assign(a, x, g2*mTAU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 2));
4106 for (
int a=0;a<4;a++) {
4107 for (
int x=0;x<6;x++) {
4112 for (
int a=0;a<2;a++) {
4113 for (
int x=0;x<3;x++) {
4118 for (
int a=0;a<4;a++) {
4119 for (
int x=0;x<6;x++) {
4120 if (fabs(1.0 -
Lepty(a, x)) > 0.01) {
4122 (2.0 - 9.0*
Lepty(a, x) + 18.0*pow(
Lepty(a, x),2.0) - 11.0*pow(
Lepty(a, x),3.0)
4123 + 6.0*pow(
Lepty(a, x),3.0)*log(
Lepty(a, x))) );
4125 2.0*pow(
Lepty(a, x),3.0) - 6.0*pow(
Lepty(a,x),2.0)*log(
Lepty(a, x))))/
4126 (6.0 * pow((1.0 -
Lepty(a,x)),4.0)) );
4128 (pow((1.0-
Lepty(a, x)),3.0)));
4132 Leptf1.assign(a, x, 1.0/12.0 - (
Lepty(a, x) - 1.0)/30.0);
4133 Leptf2.assign(a, x, 1.0/3.0 - (
Lepty(a, x) - 1.0)/6.0);
4138 for (
int a=0;a<2;a++) {
4139 for (
int x=0;x<3;x++) {
4140 if(fabs(1.0-
Leptz(a, x)) > 0.01) {
4142 (16.0 - 45.0*
Leptz(a,x) + 36.0*pow(
Leptz(a,x),2.0) - 7.0*pow(
Leptz(a,x),3.0)
4143 + 6.0*(2.0 - 3.0*
Leptz(a,x))*log(
Leptz(a,x))) );
4146 (6.0*pow((1.0 -
Leptz(a, x)),4.0))) );
4148 - 2.0*log(
Leptz(a, x)))/
4149 pow((1.0 -
Leptz(a, x)),3.0)) );
4153 Leptf3.assign(a, x, 1.0/12.0 - (
Leptz(a, x) - 1.0)/20.0 );
4154 Leptf4.assign(a, x, 2.0/3.0 - (
Leptz(a, x) - 1.0)/2.0 );
4162 for (
int a=0;a<4;a++) {
4163 for (
int x=0;x<6;x++) {
4178 gslpp::complex A1RN = 0.0;
4179 gslpp::complex A1LN = 0.0;
4180 gslpp::complex ARN = 0.0;
4181 gslpp::complex ALN = 0.0;
4182 for (
int a=0;a<4;a++) {
4183 for (
int x=0;x<6;x++) {
4191 for (
int a=0;a<2;a++) {
4192 for (
int x=0;x<3;x++) {
4207 gslpp::complex A1RC = 0.0;
4208 gslpp::complex A1LC = 0.0;
4209 gslpp::complex ARC = 0.0;
4210 gslpp::complex ALC = 0.0;
4211 for (
int a=0;a<2;a++) {
4212 for (
int x=0;x<3;x++) {
4230 for (
int a=0;a<4;a++) {
4231 for (
int x=0;x<6;x++) {
4246 gslpp::complex TauA1RN = 0.0;
4247 gslpp::complex TauA1LN = 0.0;
4248 gslpp::complex TauARN = 0.0;
4249 gslpp::complex TauALN = 0.0;
4250 for (
int a=0;a<4;a++) {
4251 for (
int x=0;x<6;x++) {
4259 for (
int a=0;a<2;a++) {
4260 for (
int x=0;x<3;x++) {
4275 gslpp::complex TauA1RC = 0.0;
4276 gslpp::complex TauA1LC = 0.0;
4277 gslpp::complex TauARC = 0.0;
4278 gslpp::complex TauALC = 0.0;
4279 for (
int a=0;a<2;a++) {
4280 for (
int x=0;x<3;x++) {
4297 for (
int a=0;a<4;a++) {
4298 for (
int x=0;x<6;x++) {
4313 gslpp::complex TEA1RN = 0.0;
4314 gslpp::complex TEA1LN = 0.0;
4315 gslpp::complex TEARN = 0.0;
4316 gslpp::complex TEALN = 0.0;
4317 for (
int a=0;a<4;a++) {
4318 for (
int x=0;x<6;x++) {
4326 for (
int a=0;a<2;a++) {
4327 for (
int x=0;x<3;x++) {
4342 gslpp::complex TEA1RC = 0.0;
4343 gslpp::complex TEA1LC = 0.0;
4344 gslpp::complex TEARC = 0.0;
4345 gslpp::complex TEALC = 0.0;
4346 for (
int a=0;a<2;a++) {
4347 for (
int x=0;x<3;x++) {
4372 double sw2 =
mySUSY.StandardModel::sW2(MW);
4373 double stw = sqrt(sw2);
4374 double ctw = sqrt(1.0 - sw2);
4375 double ttw = stw/ctw;
4380 double cdenc = sqrt(2.0)*MW*
cosb;
4381 double cdenn = MW*
cosb;
4383 double g2t = g2/sqrt(2.0);
4386 gslpp::vector<gslpp::complex>
BFunctions(8, 0.);
4389 for (
int a=0;a<4;a++) {
4390 for (
int x=0;x<6;x++) {
4392 NRlE.assign(a, x, - (g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 0) + (mE/cdenn)*
ON(a, 2)*
myRl(x, 3)));
4393 NRlMU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 1) + (mMU/cdenn)*
ON(a, 2)*
myRl(x, 4)));
4394 NRlTAU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 2) + (mTAU/cdenn)*
ON(a, 2)*
myRl(x, 5)));
4396 NLlE.assign(a, x, -(g2t)*((mE/cdenn)*
ON(a, 2)*
myRl(x, 0) + 2.0*
ON(a, 0)*ttw*
myRl(x, 3)));
4397 NLlMU.assign(a, x, -(g2t)*((mMU/cdenn)*
ON(a, 2)*
myRl(x, 1) + 2.0*
ON(a, 0)*ttw*
myRl(x, 4)));
4398 NLlTAU.assign(a, x, -(g2t)*((mTAU/cdenn)*
ON(a, 2)*
myRl(x, 2) + 2.0*
ON(a, 0)*ttw*
myRl(x, 5)));
4407 for (
int a=0;a<2;a++) {
4408 for (
int x=0;x<3;x++) {
4414 CLlE.assign(a, x, g2*mE/cdenc*
myU(a, 1).conjugate()*
myRn(x, 0));
4415 CLlMU.assign(a, x, g2*mMU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 1));
4416 CLlTAU.assign(a, x, g2*mTAU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 2));
4427 gslpp::complex B1nRMu3E = 0.0;
4428 gslpp::complex B2nRMu3E = 0.0;
4429 gslpp::complex B3nRMu3E = 0.0;
4430 gslpp::complex B4nRMu3E = 0.0;
4431 gslpp::complex B1nLMu3E = 0.0;
4432 gslpp::complex B2nLMu3E = 0.0;
4433 gslpp::complex B3nLMu3E = 0.0;
4434 gslpp::complex B4nLMu3E = 0.0;
4435 for (
int a=0;a<4;a++) {
4436 for (
int b=0;b<4;b++) {
4437 for (
int x=0;x<6;x++) {
4438 for (
int t=0;
t<6;
t++) {
4478 gslpp::complex B1cRMu3E = 0.0;
4479 gslpp::complex B2cRMu3E = 0.0;
4480 gslpp::complex B3cRMu3E = 0.0;
4481 gslpp::complex B1cLMu3E = 0.0;
4482 gslpp::complex B2cLMu3E = 0.0;
4483 gslpp::complex B3cLMu3E = 0.0;
4484 for (
int a=0;a<2;a++) {
4485 for (
int b=0;b<2;b++) {
4486 for (
int x=0;x<3;x++) {
4487 for (
int t=0;
t<3;
t++) {
4522 gslpp::complex B1nRTau3Mu = 0.0;
4523 gslpp::complex B2nRTau3Mu = 0.0;
4524 gslpp::complex B3nRTau3Mu = 0.0;
4525 gslpp::complex B4nRTau3Mu = 0.0;
4526 gslpp::complex B1nLTau3Mu = 0.0;
4527 gslpp::complex B2nLTau3Mu = 0.0;
4528 gslpp::complex B3nLTau3Mu = 0.0;
4529 gslpp::complex B4nLTau3Mu = 0.0;
4530 for (
int a=0;a<4;a++) {
4531 for (
int b=0;b<4;b++) {
4532 for (
int x=0;x<6;x++) {
4533 for (
int t=0;
t<6;
t++) {
4569 gslpp::complex B1cRTau3Mu = 0.0;
4570 gslpp::complex B2cRTau3Mu = 0.0;
4571 gslpp::complex B3cRTau3Mu = 0.0;
4572 gslpp::complex B1cLTau3Mu = 0.0;
4573 gslpp::complex B2cLTau3Mu = 0.0;
4574 gslpp::complex B3cLTau3Mu = 0.0;
4575 for (
int a=0;a<2;a++) {
4576 for (
int b=0;b<2;b++) {
4577 for (
int x=0;x<3;x++) {
4578 for (
int t=0;
t<3;
t++) {
4600 BFunctions.assign(0, B1nRTau3Mu + B1cRTau3Mu );
4601 BFunctions.assign(1, B1nLTau3Mu + B1cLTau3Mu );
4602 BFunctions.assign(2, B2nRTau3Mu + B2cRTau3Mu );
4603 BFunctions.assign(3, B2nLTau3Mu + B2cLTau3Mu );
4604 BFunctions.assign(4, B3nRTau3Mu + B3cRTau3Mu );
4605 BFunctions.assign(5, B3nLTau3Mu + B3cLTau3Mu );
4612 gslpp::complex B1nRTau3E = 0.0;
4613 gslpp::complex B2nRTau3E = 0.0;
4614 gslpp::complex B3nRTau3E = 0.0;
4615 gslpp::complex B4nRTau3E = 0.0;
4616 gslpp::complex B1nLTau3E = 0.0;
4617 gslpp::complex B2nLTau3E = 0.0;
4618 gslpp::complex B3nLTau3E = 0.0;
4619 gslpp::complex B4nLTau3E = 0.0;
4620 for (
int a=0;a<4;a++) {
4621 for (
int b=0;b<4;b++) {
4622 for (
int x=0;x<6;x++) {
4623 for (
int t=0;
t<6;
t++) {
4659 gslpp::complex B1cRTau3E = 0.0;
4660 gslpp::complex B2cRTau3E = 0.0;
4661 gslpp::complex B3cRTau3E = 0.0;
4662 gslpp::complex B1cLTau3E = 0.0;
4663 gslpp::complex B2cLTau3E = 0.0;
4664 gslpp::complex B3cLTau3E = 0.0;
4665 for (
int a=0;a<2;a++) {
4666 for (
int b=0;b<2;b++) {
4667 for (
int x=0;x<3;x++) {
4668 for (
int t=0;
t<3;
t++) {
4690 BFunctions.assign(0, B1nRTau3E + B1cRTau3E );
4691 BFunctions.assign(1, B1nLTau3E + B1cLTau3E );
4692 BFunctions.assign(2, B2nRTau3E + B2cRTau3E );
4693 BFunctions.assign(3, B2nLTau3E + B2cLTau3E );
4694 BFunctions.assign(4, B3nRTau3E + B3cRTau3E );
4695 BFunctions.assign(5, B3nLTau3E + B3cLTau3E );
4734 double piconst = 1.0/(32.0 * pi * pi);
4735 double sw2 =
mySUSY.StandardModel::sW2(MW);
4736 double stw = sqrt(sw2);
4737 double ctw = sqrt(1.0 - sw2);
4738 double ttw = stw/ctw;
4744 gslpp::complex cosa = sqrt(1.0-sina*sina);
4745 gslpp::complex sinapb = sina*
cosb+cosa*
sinb;
4746 gslpp::complex cosapb = cosa*
cosb-sina*
sinb;
4747 gslpp::complex cosbma =
cosb*cosa+
sinb*sina;
4748 gslpp::complex sinbma =
sinb*cosa-
cosb*sina;
4757 double cdenc = sqrt(2.0)*MW*
cosb;
4758 double cdenn = MW*
cosb;
4760 double g2t = g2/sqrt(2.0);
4796 for (
int a=0;a<4;a++) {
4797 for (
int x=0;x<6;x++) {
4799 NRlE.assign(a, x, - (g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 0) + (mE/cdenn)*
ON(a, 2)*
myRl(x, 3)));
4800 NRlMU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 1) + (mMU/cdenn)*
ON(a, 2)*
myRl(x, 4)));
4801 NRlTAU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 2) + (mTAU/cdenn)*
ON(a, 2)*
myRl(x, 5)));
4803 NLlE.assign(a, x, -(g2t)*((mE/cdenn)*
ON(a, 2)*
myRl(x, 0) + 2.0*
ON(a, 0)*ttw*
myRl(x, 3)));
4804 NLlMU.assign(a, x, -(g2t)*((mMU/cdenn)*
ON(a, 2)*
myRl(x, 1) + 2.0*
ON(a, 0)*ttw*
myRl(x, 4)));
4805 NLlTAU.assign(a, x, -(g2t)*((mTAU/cdenn)*
ON(a, 2)*
myRl(x, 2) + 2.0*
ON(a, 0)*ttw*
myRl(x, 5)));
4814 for (
int a=0;a<2;a++) {
4815 for (
int x=0;x<3;x++) {
4821 CLlE.assign(a, x, g2*mE/cdenc*
myU(a, 1).conjugate()*
myRn(x, 0));
4822 CLlMU.assign(a, x, g2*mMU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 1));
4823 CLlTAU.assign(a, x, g2*mTAU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 2));
4827 gslpp::vector<gslpp::complex> sigma1(3, 0.);
4828 gslpp::vector<gslpp::complex> sigma2(3, 0.);
4829 gslpp::vector<gslpp::complex> sigma3(3, 0.);
4830 gslpp::vector<gslpp::complex> sigma4(3, 0.);
4831 gslpp::vector<gslpp::complex> sigma5(3, 0.);
4832 sigma1.assign(0, sina);
4833 sigma1.assign(1, -cosa);
4834 sigma1.assign(2,
sinb*gslpp::complex::i());
4835 sigma2.assign(0, cosa);
4836 sigma2.assign(1, sina);
4837 sigma2.assign(2, -
cosb*gslpp::complex::i());
4838 sigma3.assign(0, sinapb);
4839 sigma3.assign(1, -cosapb);
4840 sigma3.assign(2, 0.);
4841 sigma4.assign(0, -sina);
4842 sigma4.assign(1, cosa);
4843 sigma4.assign(2, 0.);
4844 sigma5.assign(0, -cosbma);
4845 sigma5.assign(1, sinbma);
4846 sigma5.assign(2, cos2b*gslpp::complex::i());
4848 gslpp::matrix<gslpp::complex> Qpp(4, 4, 0.), Rpp(4, 4, 0.);
4849 gslpp::matrix<gslpp::complex>
DL0(4, 4, 0.), DR0(4, 4, 0.), DL1(4, 4, 0.), DR1(4, 4, 0.), DL2(4, 4, 0.), DR2(4, 4, 0.);
4850 for (
int a=0;a<4;a++) {
4851 for (
int b=0;b<4;b++) {
4852 Qpp.assign(a, b, 0.5*(
ON(a,2)*(
ON(b,1)-ttw*
ON(b,0))+
ON(b,2)*(
ON(a,1)-ttw*
ON(a,0))) );
4853 Rpp.assign(a, b, (M2.conjugate()*
ON(a,1)*
ON(b,1) +M1.conjugate()*
ON(a,0)*
ON(b,0) -muH.conjugate()*(
ON(a,2)*
ON(b,3)+
ON(a,3)*
ON(b,2)))/(2.0*MW) );
4854 DL0.assign(b, a, -g2/
sinb * (Qpp(a,b).conjugate()*sigma5(0) -Rpp(a,b).conjugate()*sigma2(0) +
MNeig(a)/(2.0*MW)*sigma2(0)*
delta_ab(a,b)) );
4855 DR0.assign(b, a,
DL0(b,a).conjugate() );
4856 DL1.assign(b, a, -g2/
sinb * (Qpp(a,b).conjugate()*sigma5(1) -Rpp(a,b).conjugate()*sigma2(1) +
MNeig(a)/(2.0*MW)*sigma2(1)*
delta_ab(a,b)) );
4857 DR1.assign(b, a, DL1(b,a).conjugate() );
4858 DL2.assign(b, a, -g2/
sinb * (Qpp(a,b).conjugate()*sigma5(2) -Rpp(a,b).conjugate()*sigma2(2) +
MNeig(a)/(2.0*MW)*sigma2(2)*
delta_ab(a,b)) );
4859 DR2.assign(b, a, DL2(b,a).conjugate() );
4863 gslpp::matrix<gslpp::complex> Qch(2, 2, 0.), Rch(2, 2, 0.);
4864 gslpp::matrix<gslpp::complex> WL0(2, 2, 0.), WR0(2, 2, 0.), WL1(2, 2, 0.), WR1(2, 2, 0.), WL2(2, 2, 0.), WR2(2, 2, 0.);
4865 for (
int a=0;a<2;a++) {
4866 for (
int b=0;b<2;b++) {
4867 Qch.assign(a, b,
myU(a,1)*
myV(b,0)/sqrt(2.0) );
4868 Rch.assign(a, b, (M2.conjugate()*
myU(a,0)*
myV(b,0) +muH.conjugate()*
myU(a,1)*
myV(b,1))/(2.0*MW) );
4869 WR0.assign(a, b, -g2/
sinb * (Qch(a,b)*sigma5(0).conjugate() -Rch(a,b)*sigma2(0).conjugate() +
MChi(a)/(2.0*MW)*sigma2(0).conjugate()*
delta_ab(a,b)) );
4870 WL0.assign(b, a, WR0(a,b).conjugate() );
4871 WR1.assign(a, b, -g2/
sinb * (Qch(a,b)*sigma5(1).conjugate() -Rch(a,b)*sigma2(1).conjugate() +
MChi(a)/(2.0*MW)*sigma2(1).conjugate()*
delta_ab(a,b)) );
4872 WL1.assign(b, a, WR1(a,b).conjugate() );
4873 WR2.assign(a, b, -g2/
sinb * (Qch(a,b)*sigma5(2).conjugate() -Rch(a,b)*sigma2(2).conjugate() +
MChi(a)/(2.0*MW)*sigma2(2).conjugate()*
delta_ab(a,b)) );
4874 WL2.assign(b, a, WR2(a,b).conjugate() );
4878 gslpp::vector<gslpp::complex> gLLE(3, 0.), gRRE(3, 0.), gLRE(3, 0.), gRLE(3, 0.);
4879 gslpp::vector<gslpp::complex> gLLMU(3, 0.), gRRMU(3, 0.), gLRMU(3, 0.), gRLMU(3, 0.);
4880 gslpp::vector<gslpp::complex> gLLTAU(3, 0.), gRRTAU(3, 0.), gLRTAU(3, 0.), gRLTAU(3, 0.);
4881 gslpp::vector<gslpp::complex> gLLNU(3, 0.);
4882 for (
int p=0;p<3;p++) {
4883 gLLE.assign(p, MZ/ctw*sigma3(p)*(0.5-sw2) + mE*mE/(MW*
cosb)*sigma4(p));
4884 gLLMU.assign(p, MZ/ctw*sigma3(p)*(0.5-sw2) + mMU*mMU/(MW*
cosb)*sigma4(p));
4885 gLLTAU.assign(p, MZ/ctw*sigma3(p)*(0.5-sw2) + mTAU*mTAU/(MW*
cosb)*sigma4(p));
4886 gRRE.assign(p, MZ/ctw*sigma3(p)*sw2 + mE*mE/(MW*
cosb)*sigma4(p));
4887 gRRMU.assign(p, MZ/ctw*sigma3(p)*sw2 + mMU*mMU/(MW*
cosb)*sigma4(p));
4888 gRRTAU.assign(p, MZ/ctw*sigma3(p)*sw2 + mTAU*mTAU/(MW*
cosb)*sigma4(p));
4889 gLRE.assign(p, (-sigma1(p)*
TEhat(0,0)/mE*
v1/sqrt(2.0)-sigma2(p).conjugate()*muH)*mE/(2.0*MW*
cosb));
4890 gLRMU.assign(p, (-sigma1(p)*
TEhat(1,1)/mMU*
v1/sqrt(2.0)-sigma2(p).conjugate()*muH)*mMU/(2.0*MW*
cosb));
4891 gLRTAU.assign(p, (-sigma1(p)*
TEhat(2,2)/mTAU*
v1/sqrt(2.0)-sigma2(p).conjugate()*muH)*mTAU/(2.0*MW*
cosb));
4892 gRLE.assign(p, gLRE(p).conjugate());
4893 gRLMU.assign(p, gLRMU(p).conjugate());
4894 gRLTAU.assign(p, gLRTAU(p).conjugate());
4895 gLLNU.assign(p, -0.5*MZ/ctw*sigma3(p));
4911 gslpp::matrix<gslpp::complex> Gl0(6, 6, 0.), Gl1(6, 6, 0.), Gl2(6, 6, 0.);
4912 for (
int x=0;x<6;x++) {
4913 for (
int y=0;y<6;y++) {
4914 Gl0.assign(x, y, -g2*( gLLE(0)*
myRl(x,0).conjugate()*
myRl(y,0) +gRRE(0)*
myRl(x,3).conjugate()*
myRl(y,3) +gLRE(0)*
myRl(x,0).conjugate()*
myRl(y,3) +gRLE(0)*
myRl(x,3).conjugate()*
myRl(y,0)
4915 +gLLMU(0)*
myRl(x,1).conjugate()*
myRl(y,1) +gRRMU(0)*
myRl(x,4).conjugate()*
myRl(y,4) +gLRMU(0)*
myRl(x,1).conjugate()*
myRl(y,4) +gRLMU(0)*
myRl(x,4).conjugate()*
myRl(y,1)
4916 +gLLTAU(0)*
myRl(x,2).conjugate()*
myRl(y,2) +gRRTAU(0)*
myRl(x,5).conjugate()*
myRl(y,5) +gLRTAU(0)*
myRl(x,2).conjugate()*
myRl(y,5) +gRLTAU(0)*
myRl(x,5).conjugate()*
myRl(y,2)));
4917 Gl1.assign(x, y, -g2*( gLLE(1)*
myRl(x,0).conjugate()*
myRl(y,0) +gRRE(1)*
myRl(x,3).conjugate()*
myRl(y,3) +gLRE(1)*
myRl(x,0).conjugate()*
myRl(y,3) +gRLE(1)*
myRl(x,3).conjugate()*
myRl(y,0)
4918 +gLLMU(1)*
myRl(x,1).conjugate()*
myRl(y,1) +gRRMU(1)*
myRl(x,4).conjugate()*
myRl(y,4) +gLRMU(1)*
myRl(x,1).conjugate()*
myRl(y,4) +gRLMU(1)*
myRl(x,4).conjugate()*
myRl(y,1)
4919 +gLLTAU(1)*
myRl(x,2).conjugate()*
myRl(y,2) +gRRTAU(1)*
myRl(x,5).conjugate()*
myRl(y,5) +gLRTAU(1)*
myRl(x,2).conjugate()*
myRl(y,5) +gRLTAU(1)*
myRl(x,5).conjugate()*
myRl(y,2)));
4920 Gl2.assign(x, y, -g2*( gLLE(2)*
myRl(x,0).conjugate()*
myRl(y,0) +gRRE(2)*
myRl(x,3).conjugate()*
myRl(y,3) +gLRE(2)*
myRl(x,0).conjugate()*
myRl(y,3) +gRLE(2)*
myRl(x,3).conjugate()*
myRl(y,0)
4921 +gLLMU(2)*
myRl(x,1).conjugate()*
myRl(y,1) +gRRMU(2)*
myRl(x,4).conjugate()*
myRl(y,4) +gLRMU(2)*
myRl(x,1).conjugate()*
myRl(y,4) +gRLMU(2)*
myRl(x,4).conjugate()*
myRl(y,1)
4922 +gLLTAU(2)*
myRl(x,2).conjugate()*
myRl(y,2) +gRRTAU(2)*
myRl(x,5).conjugate()*
myRl(y,5) +gLRTAU(2)*
myRl(x,2).conjugate()*
myRl(y,5) +gRLTAU(2)*
myRl(x,5).conjugate()*
myRl(y,2)));
4926 gslpp::matrix<gslpp::complex> Gnu0(3, 3, 0.), Gnu1(3, 3, 0.), Gnu2(3, 3, 0.);
4927 for (
int x=0;x<3;x++) {
4928 Gnu0.assign(x, x, -g2*gLLNU(0) );
4929 Gnu1.assign(x, x, -g2*gLLNU(1) );
4930 Gnu2.assign(x, x, -g2*gLLNU(2) );
4933 gslpp::vector<gslpp::complex> SRE(3, 0.), SLE(3, 0.), SRMU(3, 0.), SLMU(3, 0.), SRTAU(3, 0.), SLTAU(3, 0.);
4934 for (
int p=0;p<3;p++) {
4935 SRE.assign(p, g2*mE/(2.0*MW*
cosb) * sigma1(p));
4936 SLE.assign(p, g2*mE/(2.0*MW*
cosb) * sigma1(p).conjugate());
4937 SRMU.assign(p, g2*mMU/(2.0*MW*
cosb) * sigma1(p));
4938 SLMU.assign(p, g2*mMU/(2.0*MW*
cosb) * sigma1(p).conjugate());
4939 SRTAU.assign(p, g2*mTAU/(2.0*MW*
cosb) * sigma1(p));
4940 SLTAU.assign(p, g2*mTAU/(2.0*MW*
cosb) * sigma1(p).conjugate());
4946 gslpp::complex HpengMuEEENR0 = 0.0;
4947 gslpp::complex HpengMuEEENL0 = 0.0;
4948 gslpp::complex HpengMuEEENR1 = 0.0;
4949 gslpp::complex HpengMuEEENL1 = 0.0;
4950 gslpp::complex HpengMuEEENR2 = 0.0;
4951 gslpp::complex HpengMuEEENL2 = 0.0;
4952 for (
int x=0;x<6;x++) {
4953 for (
int a=0;a<4;a++) {
4954 for (
int b=0;b<4;b++) {
5058 for (
int y=0;y<6;y++) {
5148 gslpp::complex B2HiggsnR = (-0.5*HpengMuEEENR0*SLE(0)/(mh*mh)-0.5*HpengMuEEENR1*SLE(1)/(mH*mH)-0.5*HpengMuEEENR2*SLE(2)/(mA*mA))/(4.0*pi*alph);
5149 gslpp::complex B2HiggsnL = (-0.5*HpengMuEEENL0*SRE(0)/(mh*mh)-0.5*HpengMuEEENL1*SRE(1)/(mH*mH)-0.5*HpengMuEEENL2*SRE(2)/(mA*mA))/(4.0*pi*alph);
5150 gslpp::complex B3HiggsnR = (HpengMuEEENR0*SRE(0)/(mh*mh)+HpengMuEEENR1*SRE(1)/(mH*mH)+HpengMuEEENR2*SRE(2)/(mA*mA))/(4.0*pi*alph);
5151 gslpp::complex B3HiggsnL = (HpengMuEEENL0*SLE(0)/(mh*mh)+HpengMuEEENL1*SLE(1)/(mH*mH)+HpengMuEEENL2*SLE(2)/(mA*mA))/(4.0*pi*alph);
5154 gslpp::complex HpengMuEEECR0 = 0.0;
5155 gslpp::complex HpengMuEEECL0 = 0.0;
5156 gslpp::complex HpengMuEEECR1 = 0.0;
5157 gslpp::complex HpengMuEEECL1 = 0.0;
5158 gslpp::complex HpengMuEEECR2 = 0.0;
5159 gslpp::complex HpengMuEEECL2 = 0.0;
5160 for (
int x=0;x<3;x++) {
5161 for (
int a=0;a<2;a++) {
5162 for (
int b=0;b<2;b++) {
5266 for (
int y=0;y<3;y++) {
5299 HpengMuEEECR0 = HpengMuEEECR0 - 2.0*piconst*(SRMU(0)/(mE*mE-mMU*mMU)*(-
CRlE(a,x)*
CRlMU(a,x).conjugate()*mE*mE*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
5308 HpengMuEEECL0 = HpengMuEEECL0 - 2.0*piconst*(SLMU(0)/(mE*mE-mMU*mMU)*(-
CLlE(a,x)*
CLlMU(a,x).conjugate()*mE*mE*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
5317 HpengMuEEECR1 = HpengMuEEECR1 - 2.0*piconst*(SRMU(1)/(mE*mE-mMU*mMU)*(-
CRlE(a,x)*
CRlMU(a,x).conjugate()*mE*mE*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
5326 HpengMuEEECL1 = HpengMuEEECL1 - 2.0*piconst*(SLMU(1)/(mE*mE-mMU*mMU)*(-
CLlE(a,x)*
CLlMU(a,x).conjugate()*mE*mE*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
5335 HpengMuEEECR2 = HpengMuEEECR2 - 2.0*piconst*(SRMU(2)/(mE*mE-mMU*mMU)*(-
CRlE(a,x)*
CRlMU(a,x).conjugate()*mE*mE*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
5344 HpengMuEEECL2 = HpengMuEEECL2 - 2.0*piconst*(SLMU(2)/(mE*mE-mMU*mMU)*(-
CLlE(a,x)*
CLlMU(a,x).conjugate()*mE*mE*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
5356 gslpp::complex B2HiggscR = (-0.5*HpengMuEEECR0*SLE(0)/(mh*mh)-0.5*HpengMuEEECR1*SLE(1)/(mH*mH)-0.5*HpengMuEEECR2*SLE(2)/(mA*mA))/(4.0*pi*alph);
5357 gslpp::complex B2HiggscL = (-0.5*HpengMuEEECL0*SRE(0)/(mh*mh)-0.5*HpengMuEEECL1*SRE(1)/(mH*mH)-0.5*HpengMuEEECL2*SRE(2)/(mA*mA))/(4.0*pi*alph);
5358 gslpp::complex B3HiggscR = (HpengMuEEECR0*SRE(0)/(mh*mh)+HpengMuEEECR1*SRE(1)/(mH*mH)+HpengMuEEECR2*SRE(2)/(mA*mA))/(4.0*pi*alph);
5359 gslpp::complex B3HiggscL = (HpengMuEEECL0*SLE(0)/(mh*mh)+HpengMuEEECL1*SLE(1)/(mH*mH)+HpengMuEEECL2*SLE(2)/(mA*mA))/(4.0*pi*alph);
5371 gslpp::complex HpengTauMUMUMUNR0 = 0.0;
5372 gslpp::complex HpengTauMUMUMUNL0 = 0.0;
5373 gslpp::complex HpengTauMUMUMUNR1 = 0.0;
5374 gslpp::complex HpengTauMUMUMUNL1 = 0.0;
5375 gslpp::complex HpengTauMUMUMUNR2 = 0.0;
5376 gslpp::complex HpengTauMUMUMUNL2 = 0.0;
5377 for (
int x=0;x<6;x++) {
5378 for (
int a=0;a<4;a++) {
5379 for (
int b=0;b<4;b++) {
5483 for (
int y=0;y<6;y++) {
5516 HpengTauMUMUMUNR0 = HpengTauMUMUMUNR0 - 2.0*piconst*(SRTAU(0)/(mMU*mMU-mTAU*mTAU)*(-
NRlMU(a,x)*
NRlTAU(a,x).conjugate()*mMU*mMU*
PV.
B0(1.,0.,
MNeig(a)*
MNeig(a),
mym_se_sq(x))
5525 HpengTauMUMUMUNL0 = HpengTauMUMUMUNL0 - 2.0*piconst*(SLTAU(0)/(mMU*mMU-mTAU*mTAU)*(-
NLlMU(a,x)*
NLlTAU(a,x).conjugate()*mMU*mMU*
PV.
B0(1.,0.,
MNeig(a)*
MNeig(a),
mym_se_sq(x))
5534 HpengTauMUMUMUNR1 = HpengTauMUMUMUNR1 - 2.0*piconst*(SRTAU(1)/(mMU*mMU-mTAU*mTAU)*(-
NRlMU(a,x)*
NRlTAU(a,x).conjugate()*mMU*mMU*
PV.
B0(1.,0.,
MNeig(a)*
MNeig(a),
mym_se_sq(x))
5543 HpengTauMUMUMUNL1 = HpengTauMUMUMUNL1 - 2.0*piconst*(SLTAU(1)/(mMU*mMU-mTAU*mTAU)*(-
NLlMU(a,x)*
NLlTAU(a,x).conjugate()*mMU*mMU*
PV.
B0(1.,0.,
MNeig(a)*
MNeig(a),
mym_se_sq(x))
5552 HpengTauMUMUMUNR2 = HpengTauMUMUMUNR2 - 2.0*piconst*(SRTAU(2)/(mMU*mMU-mTAU*mTAU)*(-
NRlMU(a,x)*
NRlTAU(a,x).conjugate()*mMU*mMU*
PV.
B0(1.,0.,
MNeig(a)*
MNeig(a),
mym_se_sq(x))
5561 HpengTauMUMUMUNL2 = HpengTauMUMUMUNL2 - 2.0*piconst*(SLTAU(2)/(mMU*mMU-mTAU*mTAU)*(-
NLlMU(a,x)*
NLlTAU(a,x).conjugate()*mMU*mMU*
PV.
B0(1.,0.,
MNeig(a)*
MNeig(a),
mym_se_sq(x))
5573 gslpp::complex B2HiggsnRtm = (-0.5*HpengTauMUMUMUNR0*SLMU(0)/(mh*mh)-0.5*HpengTauMUMUMUNR1*SLMU(1)/(mH*mH)-0.5*HpengTauMUMUMUNR2*SLMU(2)/(mA*mA))/(4.0*pi*alph);
5574 gslpp::complex B2HiggsnLtm = (-0.5*HpengTauMUMUMUNL0*SRMU(0)/(mh*mh)-0.5*HpengTauMUMUMUNL1*SRMU(1)/(mH*mH)-0.5*HpengTauMUMUMUNL2*SRMU(2)/(mA*mA))/(4.0*pi*alph);
5575 gslpp::complex B3HiggsnRtm = (HpengTauMUMUMUNR0*SRMU(0)/(mh*mh)+HpengTauMUMUMUNR1*SRMU(1)/(mH*mH)+HpengTauMUMUMUNR2*SRMU(2)/(mA*mA))/(4.0*pi*alph);
5576 gslpp::complex B3HiggsnLtm = (HpengTauMUMUMUNL0*SLMU(0)/(mh*mh)+HpengTauMUMUMUNL1*SLMU(1)/(mH*mH)+HpengTauMUMUMUNL2*SLMU(2)/(mA*mA))/(4.0*pi*alph);
5579 gslpp::complex HpengTauMUMUMUCR0 = 0.0;
5580 gslpp::complex HpengTauMUMUMUCL0 = 0.0;
5581 gslpp::complex HpengTauMUMUMUCR1 = 0.0;
5582 gslpp::complex HpengTauMUMUMUCL1 = 0.0;
5583 gslpp::complex HpengTauMUMUMUCR2 = 0.0;
5584 gslpp::complex HpengTauMUMUMUCL2 = 0.0;
5585 for (
int x=0;x<3;x++) {
5586 for (
int a=0;a<2;a++) {
5587 for (
int b=0;b<2;b++) {
5691 for (
int y=0;y<3;y++) {
5724 HpengTauMUMUMUCR0 = HpengTauMUMUMUCR0 - 2.0*piconst*(SRTAU(0)/(mMU*mMU-mTAU*mTAU)*(-
CRlMU(a,x)*
CRlTAU(a,x).conjugate()*mMU*mMU*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
5733 HpengTauMUMUMUCL0 = HpengTauMUMUMUCL0 - 2.0*piconst*(SLTAU(0)/(mMU*mMU-mTAU*mTAU)*(-
CLlMU(a,x)*
CLlTAU(a,x).conjugate()*mMU*mMU*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
5742 HpengTauMUMUMUCR1 = HpengTauMUMUMUCR1 - 2.0*piconst*(SRTAU(1)/(mMU*mMU-mTAU*mTAU)*(-
CRlMU(a,x)*
CRlTAU(a,x).conjugate()*mMU*mMU*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
5751 HpengTauMUMUMUCL1 = HpengTauMUMUMUCL1 - 2.0*piconst*(SLTAU(1)/(mMU*mMU-mTAU*mTAU)*(-
CLlMU(a,x)*
CLlTAU(a,x).conjugate()*mMU*mMU*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
5760 HpengTauMUMUMUCR2 = HpengTauMUMUMUCR2 - 2.0*piconst*(SRTAU(2)/(mMU*mMU-mTAU*mTAU)*(-
CRlMU(a,x)*
CRlTAU(a,x).conjugate()*mMU*mMU*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
5769 HpengTauMUMUMUCL2 = HpengTauMUMUMUCL2 - 2.0*piconst*(SLTAU(2)/(mMU*mMU-mTAU*mTAU)*(-
CLlMU(a,x)*
CLlTAU(a,x).conjugate()*mMU*mMU*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
5781 gslpp::complex B2HiggscRtm = (-0.5*HpengTauMUMUMUCR0*SLMU(0)/(mh*mh)-0.5*HpengTauMUMUMUCR1*SLMU(1)/(mH*mH)-0.5*HpengTauMUMUMUCR2*SLMU(2)/(mA*mA))/(4.0*pi*alph);
5782 gslpp::complex B2HiggscLtm = (-0.5*HpengTauMUMUMUCL0*SRMU(0)/(mh*mh)-0.5*HpengTauMUMUMUCL1*SRMU(1)/(mH*mH)-0.5*HpengTauMUMUMUCL2*SRMU(2)/(mA*mA))/(4.0*pi*alph);
5783 gslpp::complex B3HiggscRtm = (HpengTauMUMUMUCR0*SRMU(0)/(mh*mh)+HpengTauMUMUMUCR1*SRMU(1)/(mH*mH)+HpengTauMUMUMUCR2*SRMU(2)/(mA*mA))/(4.0*pi*alph);
5784 gslpp::complex B3HiggscLtm = (HpengTauMUMUMUCL0*SLMU(0)/(mh*mh)+HpengTauMUMUMUCL1*SLMU(1)/(mH*mH)+HpengTauMUMUMUCL2*SLMU(2)/(mA*mA))/(4.0*pi*alph);
5796 gslpp::complex HpengTauEEENR0 = 0.0;
5797 gslpp::complex HpengTauEEENL0 = 0.0;
5798 gslpp::complex HpengTauEEENR1 = 0.0;
5799 gslpp::complex HpengTauEEENL1 = 0.0;
5800 gslpp::complex HpengTauEEENR2 = 0.0;
5801 gslpp::complex HpengTauEEENL2 = 0.0;
5802 for (
int x=0;x<6;x++) {
5803 for (
int a=0;a<4;a++) {
5804 for (
int b=0;b<4;b++) {
5908 for (
int y=0;y<6;y++) {
5998 gslpp::complex B2HiggsnRte = (-0.5*HpengTauEEENR0*SLE(0)/(mh*mh)-0.5*HpengTauEEENR1*SLE(1)/(mH*mH)-0.5*HpengTauEEENR2*SLE(2)/(mA*mA))/(4.0*pi*alph);
5999 gslpp::complex B2HiggsnLte = (-0.5*HpengTauEEENL0*SRE(0)/(mh*mh)-0.5*HpengTauEEENL1*SRE(1)/(mH*mH)-0.5*HpengTauEEENL2*SRE(2)/(mA*mA))/(4.0*pi*alph);
6000 gslpp::complex B3HiggsnRte = (HpengTauEEENR0*SRE(0)/(mh*mh)+HpengTauEEENR1*SRE(1)/(mH*mH)+HpengTauEEENR2*SRE(2)/(mA*mA))/(4.0*pi*alph);
6001 gslpp::complex B3HiggsnLte = (HpengTauEEENL0*SLE(0)/(mh*mh)+HpengTauEEENL1*SLE(1)/(mH*mH)+HpengTauEEENL2*SLE(2)/(mA*mA))/(4.0*pi*alph);
6004 gslpp::complex HpengTauEEECR0 = 0.0;
6005 gslpp::complex HpengTauEEECL0 = 0.0;
6006 gslpp::complex HpengTauEEECR1 = 0.0;
6007 gslpp::complex HpengTauEEECL1 = 0.0;
6008 gslpp::complex HpengTauEEECR2 = 0.0;
6009 gslpp::complex HpengTauEEECL2 = 0.0;
6010 for (
int x=0;x<3;x++) {
6011 for (
int a=0;a<2;a++) {
6012 for (
int b=0;b<2;b++) {
6116 for (
int y=0;y<3;y++) {
6149 HpengTauEEECR0 = HpengTauEEECR0 - 2.0*piconst*(SRTAU(0)/(mE*mE-mTAU*mTAU)*(-
CRlE(a,x)*
CRlTAU(a,x).conjugate()*mE*mE*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
6158 HpengTauEEECL0 = HpengTauEEECL0 - 2.0*piconst*(SLTAU(0)/(mE*mE-mTAU*mTAU)*(-
CLlE(a,x)*
CLlTAU(a,x).conjugate()*mE*mE*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
6167 HpengTauEEECR1 = HpengTauEEECR1 - 2.0*piconst*(SRTAU(1)/(mE*mE-mTAU*mTAU)*(-
CRlE(a,x)*
CRlTAU(a,x).conjugate()*mE*mE*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
6176 HpengTauEEECL1 = HpengTauEEECL1 - 2.0*piconst*(SLTAU(1)/(mE*mE-mTAU*mTAU)*(-
CLlE(a,x)*
CLlTAU(a,x).conjugate()*mE*mE*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
6185 HpengTauEEECR2 = HpengTauEEECR2 - 2.0*piconst*(SRTAU(2)/(mE*mE-mTAU*mTAU)*(-
CRlE(a,x)*
CRlTAU(a,x).conjugate()*mE*mE*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
6194 HpengTauEEECL2 = HpengTauEEECL2 - 2.0*piconst*(SLTAU(2)/(mE*mE-mTAU*mTAU)*(-
CLlE(a,x)*
CLlTAU(a,x).conjugate()*mE*mE*
PV.
B0(1.,0.,
MChi(a)*
MChi(a),
mym_sn_sq(x))
6206 gslpp::complex B2HiggscRte = (-0.5*HpengTauEEECR0*SLE(0)/(mh*mh)-0.5*HpengTauEEECR1*SLE(1)/(mH*mH)-0.5*HpengTauEEECR2*SLE(2)/(mA*mA))/(4.0*pi*alph);
6207 gslpp::complex B2HiggscLte = (-0.5*HpengTauEEECL0*SRE(0)/(mh*mh)-0.5*HpengTauEEECL1*SRE(1)/(mH*mH)-0.5*HpengTauEEECL2*SRE(2)/(mA*mA))/(4.0*pi*alph);
6208 gslpp::complex B3HiggscRte = (HpengTauEEECR0*SRE(0)/(mh*mh)+HpengTauEEECR1*SRE(1)/(mH*mH)+HpengTauEEECR2*SRE(2)/(mA*mA))/(4.0*pi*alph);
6209 gslpp::complex B3HiggscLte = (HpengTauEEECL0*SLE(0)/(mh*mh)+HpengTauEEECL1*SLE(1)/(mH*mH)+HpengTauEEECL2*SLE(2)/(mA*mA))/(4.0*pi*alph);
6365 double piconst = 1.0/(32.0 * pi * pi);
6366 double sw2 =
mySUSY.StandardModel::sW2(MW);
6367 double stw = sqrt(sw2);
6368 double ctw = sqrt(1.0 - sw2);
6369 double ttw = stw/ctw;
6374 double cdenc = sqrt(2.0)*MW*
cosb;
6375 double cdenn = MW*
cosb;
6377 double g2t = g2/sqrt(2.0);
6379 gslpp::vector<gslpp::complex>
FFunctions(4, 0.);
6384 for (
int a=0;a<4;a++) {
6385 for (
int x=0;x<6;x++) {
6387 NRlE.assign(a, x, - (g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 0) + (mE/cdenn)*
ON(a, 2)*
myRl(x, 3)));
6388 NRlMU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 1) + (mMU/cdenn)*
ON(a, 2)*
myRl(x, 4)));
6389 NRlTAU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 2) + (mTAU/cdenn)*
ON(a, 2)*
myRl(x, 5)));
6391 NLlE.assign(a, x, -(g2t)*((mE/cdenn)*
ON(a, 2)*
myRl(x, 0) + 2.0*
ON(a, 0)*ttw*
myRl(x, 3)));
6392 NLlMU.assign(a, x, -(g2t)*((mMU/cdenn)*
ON(a, 2)*
myRl(x, 1) + 2.0*
ON(a, 0)*ttw*
myRl(x, 4)));
6393 NLlTAU.assign(a, x, -(g2t)*((mTAU/cdenn)*
ON(a, 2)*
myRl(x, 2) + 2.0*
ON(a, 0)*ttw*
myRl(x, 5)));
6402 for (
int a=0;a<2;a++) {
6403 for (
int x=0;x<3;x++) {
6409 CLlE.assign(a, x, g2*mE/cdenc*
myU(a, 1).conjugate()*
myRn(x, 0));
6410 CLlMU.assign(a, x, g2*mMU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 1));
6411 CLlTAU.assign(a, x, g2*mTAU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 2));
6417 for (
int a=0;a<4;a++) {
6418 for (
int x=0;x<6;x++) {
6423 for (
int x=0;x<6;x++) {
6424 for (
int a=0;a<4;a++) {
6425 for (
int b=0;b<4;b++) {
6426 if (a != b && std::fabs(
Lepty(a,x)-
Lepty(b,x)) > 0.01 && std::fabs(1.0-
Lepty(a,x)) > 0.01 && std::fabs(1.0-
Lepty(b,x)) > 0.01) {
6432 else if (a != b && std::fabs(
Lepty(a,x)-
Lepty(b,x)) > 0.01 && std::fabs(1.0-
Lepty(a,x)) > 0.01 && std::fabs(1.0-
Lepty(b,x)) <= 0.01) {
6436 else if (a != b && std::fabs(
Lepty(a,x)-
Lepty(b,x)) > 0.01 && std::fabs(1.0-
Lepty(b,x)) > 0.01 && std::fabs(1.0-
Lepty(a,x)) <= 0.01) {
6440 else if ((a == b || std::fabs(
Lepty(a,x)-
Lepty(b,x)) <= 0.01) && std::fabs(1.0-
Lepty(a,x)) > 0.01 && std::fabs(1.0-
Lepty(b,x)) > 0.01) {
6453 for (
int a=0;a<2;a++) {
6454 for (
int x=0;x<3;x++) {
6459 for (
int x=0;x<3;x++) {
6460 for (
int a=0;a<2;a++) {
6461 for (
int b=0;b<2;b++) {
6462 if (a != b && std::fabs(
Leptz(a,x)-
Leptz(b,x)) > 0.01 && std::fabs(1.0-
Leptz(a,x)) > 0.01 && std::fabs(1.0-
Leptz(b,x)) > 0.01) {
6468 else if (a != b && std::fabs(
Leptz(a,x)-
Leptz(b,x)) > 0.01 && std::fabs(1.0-
Leptz(a,x)) > 0.01 && std::fabs(1.0-
Leptz(b,x)) <= 0.01) {
6472 else if (a != b && std::fabs(
Leptz(a,x)-
Leptz(b,x)) > 0.01 && std::fabs(1.0-
Leptz(b,x)) > 0.01 && std::fabs(1.0-
Leptz(a,x)) <= 0.01) {
6476 else if ((a == b || std::fabs(
Leptz(a,x)-
Leptz(b,x)) <= 0.01) && std::fabs(1.0-
Leptz(a,x)) > 0.01 && std::fabs(1.0-
Leptz(b,x)) > 0.01) {
6491 gslpp::complex ZpengMuEEENR = 0.0;
6492 gslpp::complex ZpengMuEEENL = 0.0;
6493 for (
int x=0;x<6;x++) {
6494 for (
int a=0;a<4;a++) {
6495 for (
int b=0;b<4;b++) {
6496 ZpengMuEEENR = ZpengMuEEENR -
NLlE(a,x)*
NLlMU(b,x)*piconst*
6498 ZpengMuEEENL = ZpengMuEEENL +
NRlE(a,x)*
NRlMU(b,x)*piconst*
6504 gslpp::complex ZpengMuEEEC = 0.0;
6505 for (
int x=0;x<3;x++) {
6506 for (
int a=0;a<2;a++) {
6507 for (
int b=0;b<2;b++) {
6508 ZpengMuEEEC = ZpengMuEEEC +
CRlE(a,x)*
CRlMU(b,x)*piconst*
6510 -
myV(a,1)*
myV(b,1)*gslpp::complex(std::real(
Leptfzc[x][a][b]),std::imag(
Leptfzc[x][a][b]))/2.0);
6515 FFunctions.assign(0, ZpengMuEEENR/(MZ*MZ*ctw*ctw) );
6516 FFunctions.assign(1, ZpengMuEEENR*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
6517 FFunctions.assign(2, (ZpengMuEEENL + ZpengMuEEEC)/(MZ*MZ*ctw*ctw) );
6518 FFunctions.assign(3, (ZpengMuEEENL + ZpengMuEEEC)*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
6523 gslpp::complex ZpengTauMuMuMuNR = 0.0;
6524 gslpp::complex ZpengTauMuMuMuNL = 0.0;
6525 for (
int x=0;x<6;x++) {
6526 for (
int a=0;a<4;a++) {
6527 for (
int b=0;b<4;b++) {
6528 ZpengTauMuMuMuNR = ZpengTauMuMuMuNR -
NLlMU(a,x)*
NLlTAU(b,x)*piconst*
6530 ZpengTauMuMuMuNL = ZpengTauMuMuMuNL +
NRlMU(a,x)*
NRlTAU(b,x)*piconst*
6536 gslpp::complex ZpengTauMuMuMuC = 0.0;
6537 for (
int x=0;x<3;x++) {
6538 for (
int a=0;a<2;a++) {
6539 for (
int b=0;b<2;b++) {
6540 ZpengTauMuMuMuC = ZpengTauMuMuMuC +
CRlMU(a,x)*
CRlTAU(b,x)*piconst*
6542 -
myV(a,1)*
myV(b,1)*gslpp::complex(std::real(
Leptfzc[x][a][b]),std::imag(
Leptfzc[x][a][b]))/2.0);
6547 FFunctions.assign(0, ZpengTauMuMuMuNR/(MZ*MZ*ctw*ctw) );
6548 FFunctions.assign(1, ZpengTauMuMuMuNR*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
6549 FFunctions.assign(2, (ZpengTauMuMuMuNL + ZpengTauMuMuMuC)/(MZ*MZ*ctw*ctw) );
6550 FFunctions.assign(3, (ZpengTauMuMuMuNL + ZpengTauMuMuMuC)*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
6555 gslpp::complex ZpengTauEEENR = 0.0;
6556 gslpp::complex ZpengTauEEENL = 0.0;
6557 for (
int x=0;x<6;x++) {
6558 for (
int a=0;a<4;a++) {
6559 for (
int b=0;b<4;b++) {
6560 ZpengTauEEENR = ZpengTauEEENR -
NLlE(a,x)*
NLlTAU(b,x)*piconst*
6562 ZpengTauEEENL = ZpengTauEEENL +
NRlE(a,x)*
NRlTAU(b,x)*piconst*
6568 gslpp::complex ZpengTauEEEC = 0.0;
6569 for (
int x=0;x<3;x++) {
6570 for (
int a=0;a<2;a++) {
6571 for (
int b=0;b<2;b++) {
6572 ZpengTauEEEC = ZpengTauEEEC +
CRlE(a,x)*
CRlTAU(b,x)*piconst*
6574 -
myV(a,1)*
myV(b,1)*gslpp::complex(std::real(
Leptfzc[x][a][b]),std::imag(
Leptfzc[x][a][b]))/2.0);
6579 FFunctions.assign(0, ZpengTauEEENR/(MZ*MZ*ctw*ctw) );
6580 FFunctions.assign(1, ZpengTauEEENR*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
6581 FFunctions.assign(2, (ZpengTauEEENL + ZpengTauEEEC)/(MZ*MZ*ctw*ctw) );
6582 FFunctions.assign(3, (ZpengTauEEENL + ZpengTauEEEC)*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
6587 gslpp::complex ZpengTauMuEENR = 0.0;
6588 gslpp::complex ZpengTauMuEENL = 0.0;
6589 for (
int x=0;x<6;x++) {
6590 for (
int a=0;a<4;a++) {
6591 for (
int b=0;b<4;b++) {
6592 ZpengTauMuEENR = ZpengTauMuEENR -
NLlMU(a,x)*
NLlTAU(b,x)*piconst*
6594 ZpengTauMuEENL = ZpengTauMuEENL +
NRlMU(a,x)*
NRlTAU(b,x)*piconst*
6600 gslpp::complex ZpengTauMuEEC = 0.0;
6601 for (
int x=0;x<3;x++) {
6602 for (
int a=0;a<2;a++) {
6603 for (
int b=0;b<2;b++) {
6604 ZpengTauMuEEC = ZpengTauMuEEC +
CRlMU(a,x)*
CRlTAU(b,x)*piconst*
6606 -
myV(a,1)*
myV(b,1)*gslpp::complex(std::real(
Leptfzc[x][a][b]),std::imag(
Leptfzc[x][a][b]))/2.0);
6611 FFunctions.assign(0, ZpengTauMuEENR/(MZ*MZ*ctw*ctw) );
6612 FFunctions.assign(1, ZpengTauMuEENR*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
6613 FFunctions.assign(2, (ZpengTauMuEENL + ZpengTauMuEEC)/(MZ*MZ*ctw*ctw) );
6614 FFunctions.assign(3, (ZpengTauMuEENL + ZpengTauMuEEC)*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
7022 double g1atMZ = 0.357456;
7023 double g2atMZ = 0.651721;
7029 gslpp::complex ca = sqrt(1.0-sa*sa);
7030 gslpp::complex s2a = 2.0*ca*sa;
7031 gslpp::complex c2a = ca*ca-sa*sa;
7032 double vew =
v/sqrt(2.);
7047 gslpp::complex a3tau =
TEhat(2,2)*
mySUSY.
v1()/sqrt(2.0)/mTAU;
7049 gslpp::matrix<gslpp::complex> MsQhat2(3,3,0);
7050 gslpp::matrix<gslpp::complex> MsUhat2(3,3,0);
7051 gslpp::matrix<gslpp::complex> MsDhat2(3,3,0);
7052 gslpp::matrix<gslpp::complex> MsLhat2(3,3,0);
7053 gslpp::matrix<gslpp::complex> MsEhat2(3,3,0);
7059 double msq1L = MsQhat2(0,0).real();
7060 double msq2L = MsQhat2(1,1).real();
7061 double msq3L = MsQhat2(2,2).real();
7062 double msuR = MsUhat2(0,0).real();
7063 double mscR = MsUhat2(1,1).real();
7064 double mstR = MsUhat2(2,2).real();
7065 double msdR = MsDhat2(0,0).real();
7066 double mssR = MsDhat2(1,1).real();
7067 double msbR = MsDhat2(2,2).real();
7068 double mseL = MsLhat2(0,0).real();
7069 double msmuL = MsLhat2(1,1).real();
7070 double mstauL = MsLhat2(2,2).real();
7071 double mseR = MsEhat2(0,0).real();
7072 double msmuR = MsEhat2(1,1).real();
7073 double mstauR = MsEhat2(2,2).real();
7080 b1[0]=41./6., b1[1]=-19./6., b1[2]=-7.;
7081 gi[0]=g1atMZ, gi[1]=g2atMZ, gi[2]=g3atMZ;
7083 res_g=1./gi[k]/gi[k] -b1[k]/(8.*pi*pi)*log(sqrt(msmuL)/MZ);
7084 gi[k]=sqrt(1/res_g);
7086 g1=gi[0], g2=gi[1], g3=gi[2];
7088 double alp = (g1*g1*g2*g2/(g1*g1+g2*g2))/(4.0*pi);
7091 double mzq = sqrt( 0.5*(g1*g1+g2*g2)*vew*vew );
7092 double mwq = sqrt( 0.5*(g2*g2)*vew*vew );
7093 double sw2 = g1*g1/(g1*g1+g2*g2);
7094 gslpp::complex lh1 = -sa/
cosb;
7095 gslpp::complex lh2 = ca/
cosb;
7097 double msneu2 = msmuL+0.5*mzq*mzq*c2b;
7099 TEhat.assign(1,1, 0.);
7105 gslpp::matrix<gslpp::complex> Rsmu(2,2,0.), Xm(2,2,0.);
7106 gslpp::vector<double> msmu2(2,0.);
7107 Rsmu.assign(0,0, msmuL+mmu*mmu+c2b*mzq*mzq*(-1.0/2.0+sw2) );
7110 Rsmu.assign(1,1, msmuR+mmu*mmu-c2b*mzq*mzq*sw2);
7111 Rsmu.eigensystem(Xm,msmu2);
7113 gslpp::vector<gslpp::complex> gminus2muvector=
gminus2mu();
7119 double gm21loop = 1.06978e-9;
7126 double gm2cor, res2, res3;
7127 double res01, res02;
7136 gslpp::vector<gslpp::complex> lxh1(2,0);
7137 gslpp::vector<gslpp::complex> lxh2(2,0);
7138 gslpp::vector<gslpp::complex> lxA(2,0);
7141 for(
int i=0; i<2; i++){
7142 tmp1=sqrt(2.0)*mwq/
MChi(i);
7143 lxh1.assign(i, tmp1*(
myU(i,0)*
myV(i,1)*ca +
myU(i,1)*
myV(i,0)*(-sa) ) );
7144 lxh2.assign(i, tmp1*(
myU(i,0)*
myV(i,1)*sa +
myU(i,1)*
myV(i,0)*( ca) ) );
7149 double xps, xh1, xh2;
7150 for(
int i=0; i<2; i++){
7155 res += (lA*lxA(i)).real()*
fPS(xps) + (lh1*lxh1(i)).real()*
fS(xh1)
7156 +(lh2*lxh2(i)).real()*
fS(xh2);
7158 res *= alp*alp * mmu*mmu / (8.*pi*pi*mwq*mwq*sw2);
7162 gslpp::matrix<gslpp::complex> stauM(2,2,0), sbottomM(2,2,0), stopM(2,2,0);
7163 gslpp::vector<double> mstau2(2,0), msbottom2(2,0), mstop2(2,0);
7164 gslpp::matrix<gslpp::complex> Ustau(2,2,0), Usbottom(2,2,0), Ustop(2,2,0);
7166 stauM.assign(0,0, mstauL+mTAU*mTAU + mzq*mzq*c2b*(0.5-(-1.)*sw2) );
7167 stauM.assign(1,1, mstauR+mTAU*mTAU + mzq*mzq*c2b*(-1.)*sw2);
7168 stauM.assign(0,1, mTAU*(a3tau-muH*
tanb));
7169 stauM.assign(1,0, mTAU*(a3tau-muH*
tanb));
7171 sbottomM.assign(0,0, msq3L+mb*mb + mzq*mzq*c2b*(0.5-(-1./3.)*sw2) );
7172 sbottomM.assign(1,1, msbR+mb*mb + mzq*mzq*c2b*(-1./3.)*sw2);
7173 sbottomM.assign(0,1, mb*(a3b-muH*
tanb));
7174 sbottomM.assign(1,0, mb*(a3b-muH*
tanb));
7176 stopM.assign(0,0, msq3L+mt*mt + mzq*mzq*c2b*(0.5-(2./3.)*sw2) );
7177 stopM.assign(1,1, mstR+mt*mt + mzq*mzq*c2b*(2./3.)*sw2);
7178 stopM.assign(0,1, mt*(a3t-muH/
tanb));
7179 stopM.assign(1,0, mt*(a3t-muH/
tanb));
7181 stauM.eigensystem(Ustau,mstau2);
7182 sbottomM.eigensystem(Usbottom,msbottom2);
7183 stopM.eigensystem(Ustop,mstop2);
7185 Ustau=Ustau.hconjugate();
7186 Usbottom=Usbottom.hconjugate();
7187 Ustop=Ustop.hconjugate();
7189 gslpp::vector<gslpp::complex> lstauh1(2,0);
7190 gslpp::vector<gslpp::complex> lstauh2(2,0);
7191 gslpp::vector<gslpp::complex> lsbottomh1(2,0);
7192 gslpp::vector<gslpp::complex> lsbottomh2(2,0);
7193 gslpp::vector<gslpp::complex> lstoph1(2,0);
7194 gslpp::vector<gslpp::complex> lstoph2(2,0);
7197 for(
int i=0; i<2; i++){
7199 rr=2.*mTAU/(
cosb*mstau2(i)) * Ustau(i,0).conjugate() * Ustau(i,1);
7200 lstauh1.assign(i,rr*(-muH*ca + a3tau*(-sa)) );
7201 lstauh2.assign(i,rr*(-muH*sa + a3tau*(ca)) );
7203 rr=2.*mb/(
cosb*msbottom2(i)) * Usbottom(i,0).conjugate() * Usbottom(i,1);
7204 lsbottomh1.assign(i,rr*(-muH*ca + a3b*(-sa)) );
7205 lsbottomh2.assign(i,rr*(-muH*sa + a3b*(ca)) );
7207 rr=2.*mt/(
cosb*mstop2(i)) * Ustop(i,0).conjugate() * Ustop(i,1);
7208 lstoph1.assign(i,rr*(-muH*sa + a3t*(ca)) );
7209 lstoph2.assign(i,rr*(-muH*(-ca) + a3t*(sa)) );
7218 for(
int i=0;i<2;i++){
7219 xx1=mstau2(i)/(mh*mh);
7220 xx2=mstau2(i)/(mhh*mhh);
7221 res1 += (lh1*lstauh1(i)).real()*
fft(xx1)
7222 +(lh2*lstauh2(i)).real()*
fft(xx2);
7226 xx1=msbottom2(i)/(mh*mh);
7227 xx2=msbottom2(i)/(mhh*mhh);
7228 res1 += (lh1*lsbottomh1(i)).real()*qe2*
fft(xx1)
7229 +(lh2*lsbottomh2(i)).real()*qe2*
fft(xx2);
7233 xx1=mstop2(i)/(mh*mh);
7234 xx2=mstop2(i)/(mhh*mhh);
7235 res1 += (lh1*lstoph1(i)).real()*qe2*
fft(xx1)
7236 +(lh2*lstoph2(i)).real()*qe2*
fft(xx2);
7239 res1 *= alp*alp * mmu*mmu / (8.*pi*pi*mwq*mwq*sw2);
7251 gm2cor=0, res2=0, res3=0;
7262 double x0, x1a, x2a, xL, xR;
7265 x0=sqrt(std::fabs(msneu2));
7266 xL=sqrt( msmuL-mzq*mzq*(sw2-0.5) );
7267 xR=sqrt( msmuR +mzq*mzq*sw2 );
7269 tmp2=M2.abs2()+muH.abs2() +2.*mwq*mwq;
7270 tmp3= tmp2*tmp2 -4.*M2.abs2()*muH.abs2();
7271 x1a=sqrt( 0.5*( tmp2-sqrt(tmp3) ) );
7272 x2a=sqrt( 0.5*( tmp2+sqrt(tmp3) ) );
7274 dmu=-muH.real()*
tanb*g2*g2*M2.real()/(16.*pi*pi)*(
It(x1a,x2a,x0) + 0.5*
It(x1a,x2a,xL))
7275 -muH.real()*
tanb*g1*g1*M1.real()/(16.*pi*pi)*(
It(muH.real(),M1.real(),xR) - 0.5*
It(muH.real(),M1.real(),xL) -
It(M1.real(),xL,xR));
7278 gm2cor=gm21loop/(1+dmu);
7287 res2=gm2cor*alp/(4.*pi)*16.*log(mmu/sqrt(msmuL));
7296 if (sub_leading==1){
7300 gslpp::vector<gslpp::complex> ckL(2,0), ckR(2,0);
7301 double ymu=mmu/(vew*
cosb);
7303 ckR.assign(0, ymu*
myU(0,1));
7304 ckR.assign(1, ymu*
myU(1,1));
7307 ckL.assign(0, -g2*
myV(0,0));
7308 ckL.assign(1, -g2*
myV(1,0));
7311 for(
int i=0; i<2; i++){
7313 amch=amch -(47.*mmu/(12.*msneu2))*( (ckL(i).abs()*ckL(i).abs() + ckR(i).abs()*ckR(i).abs())*
F3C(xk) );
7314 amch=amch -(122.*
MChi(i)/(9.*msneu2))*
F4C(xk)*( (ckL(i)*ckR(i)).real() );
7316 amch=(1./(1+dmu))*amch*mmu/(16.*pi*pi)*alp/(4.*pi);
7321 gslpp::vector<gslpp::complex> nR1(4,0), nR2(4,0), nL1(4,0), nL2(4,0);
7323 for(
int i=0; i<4; i++){
7324 nR1.assign(i, sqrt(2.)*g1*
myN(i,0)*Xm(0,1) + ymu*
myN(i,2)*Xm(0,0));
7325 nR2.assign(i, sqrt(2.)*g1*
myN(i,0)*Xm(1,1) + ymu*
myN(i,2)*Xm(1,0));
7326 nL1.assign(i, (1./sqrt(2.))*(g2*
myN(i,1) + g1*
myN(i,0))*Xm(0,0).conjugate() -ymu*
myN(i,2)*Xm(0,1).conjugate());
7327 nL2.assign(i, (1./sqrt(2.))*(g2*
myN(i,1) + g1*
myN(i,0))*Xm(1,0).conjugate() -ymu*
myN(i,2)*Xm(1,1).conjugate());
7333 for(
int i=0; i<4; i++){
7336 tmp4=nL1(i).abs2() + nR1(i).abs2();
7337 tmp5=nL2(i).abs2() + nR2(i).abs2();
7338 r1=35.*mmu/(72.*msmu2(0))*
F3N(xi1)*tmp4 - 16.*
MChi0(i)/(9.*msmu2(0))*
F4N(xi1)*( (nL1(i)*nR1(i)).real() );
7339 r2=35.*mmu/(72.*msmu2(1))*
F3N(xi2)*tmp5 - 16.*
MChi0(i)/(9.*msmu2(1))*
F4N(xi2)*( (nL2(i)*nR2(i)).real() );
7343 amne*=(1./(1.+dmu))*mmu/(16.*pi*pi)*alp/(4.*pi);
7347 res2=res2+amch+amne;
7360 double x1,y1,x2,y2,x3,y3,x4,y4,x5,y5;
7361 double awhn=mmu*mmu*muH.real()*
tanb/(1+dmu);
7368 x1=M2.abs2()/msneu2;
7369 y1=muH.abs2()/msneu2;
7370 awhn *= g2*g2*M2.real()/(8.*pi*pi*msneu2*msneu2)*
Fa(x1,y1);
7374 y2=muH.abs2()/msmuL;
7375 awhl *= -g2*g2*M2.real()/(16.*pi*pi*pow(msmuL,2))*
Fb(x2,y2);
7379 y3=muH.abs2()/msmuL;
7380 abhl *= g1*g1*M1.real()/(16.*pi*pi*pow(msmuL,2))*
Fb(x3,y3);
7384 y4=muH.abs2()/msmuR;
7385 abhr *= -g1*g1*M1.real()/(8.*pi*pi*pow(msmuR,2))*
Fb(x4,y4);
7388 x5=msmuL/(M1.abs2());
7389 y5=msmuR/(M1.abs2());
7390 ablr *= g1*g1/(8.*pi*pi*pow(M1.real(),3))*
Fb(x5,y5);
7392 double dg2, dg1, dh, dwh, dbh, dtb;
7393 double msusy=sqrt(msmuL);
7395 dg1=g1*g1/(16.*pi*pi)*(4./3.)
7396 *(4./3.*log(sqrt(msuR)/msusy) + 4./3.*log(sqrt(mscR)/msusy) + 4./3.*log(sqrt(mstR)/msusy)
7397 + 2./3.*log(sqrt(msdR)/msusy) + 2./3.*log(sqrt(mssR)/msusy) + 2./3.*log(sqrt(msbR)/msusy)
7398 + 1./3.*log(sqrt(msq1L)/msusy) + 1./3.*log(sqrt(msq2L)/msusy) + 1./3.*log(sqrt(msq3L)/msusy)
7399 + log(sqrt(mseR)/msusy) + log(sqrt(mstauR)/msusy)
7400 + 1./2.*log(sqrt(mseL)/msusy) + 1./2.*log(sqrt(mstauL)/msusy) );
7402 dg2=g2*g2/(16.*pi*pi)*(4./3.)
7403 *(3./2.*log(sqrt(msq1L)/msusy) + 3./2.*log(sqrt(msq2L)/msusy) + 3./2.*log(sqrt(msq3L)/msusy)
7404 + 1./2.*log(sqrt(mseL)/msusy) + 1./2.*log(sqrt(mstauL)/msusy) );
7407 double yb, ytau, yt;
7408 double as_mt, delta_mt;
7413 as_mt=as_mt - (-7.)*log(mt/sqrt(msmuL))/(8.*pi*pi);
7414 as_mt=(1./as_mt)/(4.*pi);
7416 delta_mt=-4./3.*(as_mt/pi)-9.1*pow( (as_mt/pi), 2)-80.*pow((as_mt/pi),3);
7419 yt=mt/(vew*
sinb)*(1+delta_mt);
7420 ytau=mTAU/(vew*
cosb);
7423 yt=yt*(1+as_mt/(8.*pi)*(-4./3.));
7424 yb=yb*(1+as_mt/(8.*pi)*(-4./3.));
7426 dh=0.5/(16.*pi*pi)*(3.*yt*yt*log(sqrt(mstR)/msusy) +3.*yb*yb*log(sqrt(msbR)/msusy)
7427 +3.*(yt*yt+yb*yb)*log(sqrt(msq3L)/msusy)
7428 +ytau*ytau*log(sqrt(mstauR)/msusy) +ytau*ytau*log(sqrt(mstauL)/msusy));
7430 dwh=yt*yt/(16.*pi*pi)*(-6.*log(sqrt(msq3L)/msusy));
7432 dbh=yt*yt/(16.*pi*pi)*( 2.*log(sqrt(msq3L)/msusy)-8.*log(sqrt(mstR)/msusy) );
7436 dtb=1./(16.*pi*pi)*( 3.*yb*yb -3.*yt*yt +ytau*ytau)*log(
Q_S/msusy);
7439 res3= awhn*(dg2 + dh + dwh + dtb)
7440 + awhl*(dg2 + dh + dwh + dtb)
7441 + abhl*(dg1 + dh + dbh + dtb)
7442 + abhr*(dg1 + dh + dbh + dtb)
7452 res01=gm2cor+res2+res3-gm21loop;