20 mcmueconv(8,
NDR,
LO),
21 mcgminus2mu(2,
NDR,
LO),
130 for (
int a=0;a<4;a++) {
131 Mdiag.assign(a,a,
MChi0(a));
135 MN_tmp =
myN.transpose() * Mdiag *
myN;
137 MN_tmp.eigensystem(cmplxONT,
MNeig);
138 for (
int a=0;a<4;a++) {
139 for (
int b=0;b<4;b++) {
140 ON.assign(a,b,cmplxONT(b,a).real());
154 double piconst = 1.0/(32.0 * pi * pi);
155 double sw2 =
mySUSY.StandardModel::sW2(MW);
156 double stw =
sqrt(sw2);
157 double ctw =
sqrt(1.0 - sw2);
158 double ttw = stw/ctw;
164 double cdenn = MW*
cosb;
166 double g2t = g2/
sqrt(2.0);
171 for (
int a=0;a<4;a++) {
172 for (
int x=0;x<6;x++) {
174 NRlE.assign(a, x, - (g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 0) + (mE/cdenn)*
ON(a, 2)*
myRl(x, 3)));
175 NRlMU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 1) + (mMU/cdenn)*
ON(a, 2)*
myRl(x, 4)));
176 NRlTAU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 2) + (mTAU/cdenn)*
ON(a, 2)*
myRl(x, 5)));
178 NLlE.assign(a, x, -(g2t)*((mE/cdenn)*
ON(a, 2)*
myRl(x, 0) + 2.0*
ON(a, 0)*ttw*
myRl(x, 3)));
179 NLlMU.assign(a, x, -(g2t)*((mMU/cdenn)*
ON(a, 2)*
myRl(x, 1) + 2.0*
ON(a, 0)*ttw*
myRl(x, 4)));
180 NLlTAU.assign(a, x, -(g2t)*((mTAU/cdenn)*
ON(a, 2)*
myRl(x, 2) + 2.0*
ON(a, 0)*ttw*
myRl(x, 5)));
189 for (
int a=0;a<2;a++) {
190 for (
int x=0;x<3;x++) {
196 CLlE.assign(a, x, g2*mE/cdenc*
myU(a, 1).conjugate()*
myRn(x, 0));
197 CLlMU.assign(a, x, g2*mMU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 1));
198 CLlTAU.assign(a, x, g2*mTAU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 2));
203 for (
int a=0;a<4;a++) {
204 for (
int x=0;x<6;x++) {
209 for (
int a=0;a<2;a++) {
210 for (
int x=0;x<3;x++) {
215 for (
int a=0;a<4;a++) {
216 for (
int x=0;x<6;x++) {
217 if (fabs(1.0 -
Lepty(a, x)) > 0.01) {
223 (6.0 *
pow((1.0 -
Lepty(a,x)),4.0)) );
235 for (
int a=0;a<2;a++) {
236 for (
int x=0;x<3;x++) {
237 if(fabs(1.0-
Leptz(a, x)) > 0.01) {
243 (6.0*
pow((1.0 -
Leptz(a, x)),4.0))) );
259 for (
int a=0;a<4;a++) {
260 for (
int x=0;x<6;x++) {
279 for (
int a=0;a<4;a++) {
280 for (
int x=0;x<6;x++) {
288 for (
int a=0;a<2;a++) {
289 for (
int x=0;x<3;x++) {
308 for (
int a=0;a<2;a++) {
309 for (
int x=0;x<3;x++) {
317 AFunctions.assign(0, A1RN + A1RC );
318 AFunctions.assign(1, A1LN + A1LC );
319 AFunctions.assign(2, ARN + ARC );
320 AFunctions.assign(3, ALN + ALC );
327 for (
int a=0;a<4;a++) {
328 for (
int x=0;x<6;x++) {
347 for (
int a=0;a<4;a++) {
348 for (
int x=0;x<6;x++) {
356 for (
int a=0;a<2;a++) {
357 for (
int x=0;x<3;x++) {
376 for (
int a=0;a<2;a++) {
377 for (
int x=0;x<3;x++) {
385 AFunctions.assign(0, TauA1RC + TauA1RN );
386 AFunctions.assign(1, TauA1LC + TauA1LN );
387 AFunctions.assign(2, TauARC + TauARN );
388 AFunctions.assign(3, TauALC + TauALN );
394 for (
int a=0;a<4;a++) {
395 for (
int x=0;x<6;x++) {
414 for (
int a=0;a<4;a++) {
415 for (
int x=0;x<6;x++) {
423 for (
int a=0;a<2;a++) {
424 for (
int x=0;x<3;x++) {
443 for (
int a=0;a<2;a++) {
444 for (
int x=0;x<3;x++) {
452 AFunctions.assign(0, TEA1RC + TEA1RN );
453 AFunctions.assign(1, TEA1LC + TEA1LN );
454 AFunctions.assign(2, TEARC + TEARN );
455 AFunctions.assign(3, TEALC + TEALN );
469 double sw2 =
mySUSY.StandardModel::sW2(MW);
470 double stw =
sqrt(sw2);
471 double ctw =
sqrt(1.0 - sw2);
472 double ttw = stw/ctw;
478 double cdenn = MW*
cosb;
480 double g2t = g2/
sqrt(2.0);
486 for (
int a=0;a<4;a++) {
487 for (
int x=0;x<6;x++) {
489 NRlE.assign(a, x, - (g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 0) + (mE/cdenn)*
ON(a, 2)*
myRl(x, 3)));
490 NRlMU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 1) + (mMU/cdenn)*
ON(a, 2)*
myRl(x, 4)));
491 NRlTAU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 2) + (mTAU/cdenn)*
ON(a, 2)*
myRl(x, 5)));
493 NLlE.assign(a, x, -(g2t)*((mE/cdenn)*
ON(a, 2)*
myRl(x, 0) + 2.0*
ON(a, 0)*ttw*
myRl(x, 3)));
494 NLlMU.assign(a, x, -(g2t)*((mMU/cdenn)*
ON(a, 2)*
myRl(x, 1) + 2.0*
ON(a, 0)*ttw*
myRl(x, 4)));
495 NLlTAU.assign(a, x, -(g2t)*((mTAU/cdenn)*
ON(a, 2)*
myRl(x, 2) + 2.0*
ON(a, 0)*ttw*
myRl(x, 5)));
504 for (
int a=0;a<2;a++) {
505 for (
int x=0;x<3;x++) {
511 CLlE.assign(a, x, g2*mE/cdenc*
myU(a, 1).conjugate()*
myRn(x, 0));
512 CLlMU.assign(a, x, g2*mMU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 1));
513 CLlTAU.assign(a, x, g2*mTAU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 2));
532 for (
int a=0;a<4;a++) {
533 for (
int b=0;b<4;b++) {
534 for (
int x=0;x<6;x++) {
535 for (
int t=0;t<6;t++) {
536 B1nRMu3E = B1nRMu3E + (1.0/(4.0*pi*alph))*(0.5*
NLlMU(a,x)*
NLlE(a,t)*
NLlE(b,t)*
NLlE(b,x)
540 B2nRMu3E = B2nRMu3E + (1.0/(4.0*pi*alph))*(0.25*(
NLlMU(a,x)*
NLlE(a,t)*
NRlE(b,t)*
NRlE(b,x)
551 B1nLMu3E = B1nLMu3E + (1.0/(4.0*pi*alph))*(0.5*
NRlMU(a,x)*
NRlE(a,t)*
NRlE(b,t)*
NRlE(b,x)
555 B2nLMu3E = B2nLMu3E + (1.0/(4.0*pi*alph))*(0.25*(
NRlMU(a,x)*
NRlE(a,t)*
NLlE(b,t)*
NLlE(b,x)
581 for (
int a=0;a<2;a++) {
582 for (
int b=0;b<2;b++) {
583 for (
int x=0;x<3;x++) {
584 for (
int t=0;t<3;t++) {
585 B1cRMu3E = B1cRMu3E + (1.0/(8.0*pi*alph))*
CLlMU(a,x)*
CLlE(a,t)*
CLlE(b,t)*
CLlE(b,x)
587 B2cRMu3E = B2cRMu3E + (1.0/(4.0*pi*alph))*(0.25*
CLlMU(a,x)*
CLlE(a,t)*
CRlE(b,t)*
CRlE(b,x)
593 B1cLMu3E = B1cLMu3E + (1.0/(8.0*pi*alph))*(
CRlMU(a,x)*
CRlE(a,t)*
CRlE(b,t)*
CRlE(b,x)
595 B2cLMu3E = B2cLMu3E + (1.0/(4.0*pi*alph))*(0.25*
CRlMU(a,x)*
CRlE(a,t)*
CLlE(b,t)*
CLlE(b,x)
607 BFunctions.assign(0, B1nRMu3E + B1cRMu3E );
608 BFunctions.assign(1, B1nLMu3E + B1cLMu3E );
609 BFunctions.assign(2, B2nRMu3E + B2cRMu3E );
610 BFunctions.assign(3, B2nLMu3E + B2cLMu3E );
611 BFunctions.assign(4, B3nRMu3E + B3cRMu3E );
612 BFunctions.assign(5, B3nLMu3E + B3cLMu3E );
613 BFunctions.assign(6, B4nRMu3E );
614 BFunctions.assign(7, B4nLMu3E );
627 for (
int a=0;a<4;a++) {
628 for (
int b=0;b<4;b++) {
629 for (
int x=0;x<6;x++) {
630 for (
int t=0;t<6;t++) {
672 for (
int a=0;a<2;a++) {
673 for (
int b=0;b<2;b++) {
674 for (
int x=0;x<3;x++) {
675 for (
int t=0;t<3;t++) {
697 BFunctions.assign(0, B1nRTau3Mu + B1cRTau3Mu );
698 BFunctions.assign(1, B1nLTau3Mu + B1cLTau3Mu );
699 BFunctions.assign(2, B2nRTau3Mu + B2cRTau3Mu );
700 BFunctions.assign(3, B2nLTau3Mu + B2cLTau3Mu );
701 BFunctions.assign(4, B3nRTau3Mu + B3cRTau3Mu );
702 BFunctions.assign(5, B3nLTau3Mu + B3cLTau3Mu );
703 BFunctions.assign(6, B4nRTau3Mu );
704 BFunctions.assign(7, B4nLTau3Mu );
717 for (
int a=0;a<4;a++) {
718 for (
int b=0;b<4;b++) {
719 for (
int x=0;x<6;x++) {
720 for (
int t=0;t<6;t++) {
721 B1nRTau3E = B1nRTau3E + (1.0/(4.0*pi*alph))*(0.5*
NLlTAU(a,x)*
NLlE(a,t)*
NLlE(b,t)*
NLlE(b,x)
725 B2nRTau3E = B2nRTau3E + (1.0/(4.0*pi*alph))*(0.25*(
NLlTAU(a,x)*
NLlE(a,t)*
NRlE(b,t)*
NRlE(b,x)
736 B1nLTau3E = B1nLTau3E + (1.0/(4.0*pi*alph))*(0.5*
NRlTAU(a,x)*
NRlE(a,t)*
NRlE(b,t)*
NRlE(b,x)
740 B2nLTau3E = B2nLTau3E + (1.0/(4.0*pi*alph))*(0.25*(
NRlTAU(a,x)*
NRlE(a,t)*
NLlE(b,t)*
NLlE(b,x)
762 for (
int a=0;a<2;a++) {
763 for (
int b=0;b<2;b++) {
764 for (
int x=0;x<3;x++) {
765 for (
int t=0;t<3;t++) {
768 B2cRTau3E = B2cRTau3E + (1.0/(4.0*pi*alph))*(0.25*
CLlTAU(a,x)*
CLlE(a,t)*
CRlE(b,t)*
CRlE(b,x)
776 B2cLTau3E = B2cLTau3E + (1.0/(4.0*pi*alph))*(0.25*
CRlTAU(a,x)*
CRlE(a,t)*
CLlE(b,t)*
CLlE(b,x)
787 BFunctions.assign(0, B1nRTau3E + B1cRTau3E );
788 BFunctions.assign(1, B1nLTau3E + B1cLTau3E );
789 BFunctions.assign(2, B2nRTau3E + B2cRTau3E );
790 BFunctions.assign(3, B2nLTau3E + B2cLTau3E );
791 BFunctions.assign(4, B3nRTau3E + B3cRTau3E );
792 BFunctions.assign(5, B3nLTau3E + B3cLTau3E );
793 BFunctions.assign(6, B4nRTau3E );
794 BFunctions.assign(7, B4nLTau3E );
808 std::cout<<
"BFunctions("<<li_to_lj<<
") = "<<BFunctions<<std::endl;
831 double piconst = 1.0/(32.0 * pi * pi);
832 double sw2 =
mySUSY.StandardModel::sW2(MW);
833 double stw =
sqrt(sw2);
834 double ctw =
sqrt(1.0 - sw2);
835 double ttw = stw/ctw;
855 double cdenn = MW*
cosb;
857 double g2t = g2/
sqrt(2.0);
893 for (
int a=0;a<4;a++) {
894 for (
int x=0;x<6;x++) {
896 NRlE.assign(a, x, - (g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 0) + (mE/cdenn)*
ON(a, 2)*
myRl(x, 3)));
897 NRlMU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 1) + (mMU/cdenn)*
ON(a, 2)*
myRl(x, 4)));
898 NRlTAU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 2) + (mTAU/cdenn)*
ON(a, 2)*
myRl(x, 5)));
900 NLlE.assign(a, x, -(g2t)*((mE/cdenn)*
ON(a, 2)*
myRl(x, 0) + 2.0*
ON(a, 0)*ttw*
myRl(x, 3)));
901 NLlMU.assign(a, x, -(g2t)*((mMU/cdenn)*
ON(a, 2)*
myRl(x, 1) + 2.0*
ON(a, 0)*ttw*
myRl(x, 4)));
902 NLlTAU.assign(a, x, -(g2t)*((mTAU/cdenn)*
ON(a, 2)*
myRl(x, 2) + 2.0*
ON(a, 0)*ttw*
myRl(x, 5)));
911 for (
int a=0;a<2;a++) {
912 for (
int x=0;x<3;x++) {
918 CLlE.assign(a, x, g2*mE/cdenc*
myU(a, 1).conjugate()*
myRn(x, 0));
919 CLlMU.assign(a, x, g2*mMU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 1));
920 CLlTAU.assign(a, x, g2*mTAU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 2));
929 sigma1.assign(0, sina);
930 sigma1.assign(1, -cosa);
932 sigma2.assign(0, cosa);
933 sigma2.assign(1, sina);
935 sigma3.assign(0, sinapb);
936 sigma3.assign(1, -cosapb);
937 sigma3.assign(2, 0.);
938 sigma4.assign(0, -sina);
939 sigma4.assign(1, cosa);
940 sigma4.assign(2, 0.);
941 sigma5.assign(0, -cosbma);
942 sigma5.assign(1, sinbma);
946 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.);
947 for (
int a=0;a<4;a++) {
948 for (
int b=0;b<4;b++) {
949 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))) );
950 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) );
951 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)) );
952 DR0.assign(b, a, DL0(b,a).conjugate() );
953 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)) );
954 DR1.assign(b, a, DL1(b,a).conjugate() );
955 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)) );
956 DR2.assign(b, a, DL2(b,a).conjugate() );
961 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.);
962 for (
int a=0;a<2;a++) {
963 for (
int b=0;b<2;b++) {
964 Qch.assign(a, b,
myU(a,1)*
myV(b,0)/
sqrt(2.0) );
966 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)) );
967 WL0.assign(b, a, WR0(a,b).conjugate() );
968 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)) );
969 WL1.assign(b, a, WR1(a,b).conjugate() );
970 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)) );
971 WL2.assign(b, a, WR2(a,b).conjugate() );
979 for (
int p=0;p<3;p++) {
980 gLLE.assign(p, MZ/ctw*sigma3(p)*(0.5-sw2) + mE*mE/(MW*cosb)*sigma4(p));
981 gLLMU.assign(p, MZ/ctw*sigma3(p)*(0.5-sw2) + mMU*mMU/(MW*cosb)*sigma4(p));
982 gLLTAU.assign(p, MZ/ctw*sigma3(p)*(0.5-sw2) + mTAU*mTAU/(MW*cosb)*sigma4(p));
983 gRRE.assign(p, MZ/ctw*sigma3(p)*sw2 + mE*mE/(MW*cosb)*sigma4(p));
984 gRRMU.assign(p, MZ/ctw*sigma3(p)*sw2 + mMU*mMU/(MW*cosb)*sigma4(p));
985 gRRTAU.assign(p, MZ/ctw*sigma3(p)*sw2 + mTAU*mTAU/(MW*cosb)*sigma4(p));
986 gLRE.assign(p, (-sigma1(p)*
TEhat(0,0)/mE*
v1/
sqrt(2.0)-sigma2(p).conjugate()*muH)*mE/(2.0*MW*cosb));
987 gLRMU.assign(p, (-sigma1(p)*
TEhat(1,1)/mMU*
v1/
sqrt(2.0)-sigma2(p).conjugate()*muH)*mMU/(2.0*MW*cosb));
988 gLRTAU.assign(p, (-sigma1(p)*
TEhat(2,2)/mTAU*
v1/
sqrt(2.0)-sigma2(p).conjugate()*muH)*mTAU/(2.0*MW*cosb));
989 gRLE.assign(p, gLRE(p).conjugate());
990 gRLMU.assign(p, gLRMU(p).conjugate());
991 gRLTAU.assign(p, gLRTAU(p).conjugate());
992 gLLNU.assign(p, -0.5*MZ/ctw*sigma3(p));
1009 for (
int x=0;x<6;x++) {
1010 for (
int y=0;y<6;y++) {
1011 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)
1012 +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)
1013 +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)));
1014 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)
1015 +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)
1016 +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)));
1017 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)
1018 +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)
1019 +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)));
1024 for (
int x=0;x<3;x++) {
1025 Gnu0.assign(x, x, -g2*gLLNU(0) );
1026 Gnu1.assign(x, x, -g2*gLLNU(1) );
1027 Gnu2.assign(x, x, -g2*gLLNU(2) );
1030 gslpp::vector<gslpp::complex> SRE(3, 0.), SLE(3, 0.), SRMU(3, 0.), SLMU(3, 0.), SRTAU(3, 0.), SLTAU(3, 0.);
1031 for (
int p=0;p<3;p++) {
1032 SRE.assign(p, g2*mE/(2.0*MW*cosb) * sigma1(p));
1033 SLE.assign(p, g2*mE/(2.0*MW*cosb) * sigma1(p).conjugate());
1034 SRMU.assign(p, g2*mMU/(2.0*MW*cosb) * sigma1(p));
1035 SLMU.assign(p, g2*mMU/(2.0*MW*cosb) * sigma1(p).conjugate());
1036 SRTAU.assign(p, g2*mTAU/(2.0*MW*cosb) * sigma1(p));
1037 SLTAU.assign(p, g2*mTAU/(2.0*MW*cosb) * sigma1(p).conjugate());
1049 for (
int x=0;x<6;x++) {
1050 for (
int a=0;a<4;a++) {
1051 for (
int b=0;b<4;b++) {
1155 for (
int y=0;y<6;y++) {
1245 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);
1246 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);
1247 gslpp::complex B3HiggsnR = (HpengMuEEENR0*SRE(0)/(mh*mh)+HpengMuEEENR1*SRE(1)/(mH*mH)+HpengMuEEENR2*SRE(2)/(mA*mA))/(4.0*pi*alph);
1248 gslpp::complex B3HiggsnL = (HpengMuEEENL0*SLE(0)/(mh*mh)+HpengMuEEENL1*SLE(1)/(mH*mH)+HpengMuEEENL2*SLE(2)/(mA*mA))/(4.0*pi*alph);
1257 for (
int x=0;x<3;x++) {
1258 for (
int a=0;a<2;a++) {
1259 for (
int b=0;b<2;b++) {
1363 for (
int y=0;y<3;y++) {
1396 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))
1405 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))
1414 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))
1423 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))
1432 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))
1441 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))
1453 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);
1454 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);
1455 gslpp::complex B3HiggscR = (HpengMuEEECR0*SRE(0)/(mh*mh)+HpengMuEEECR1*SRE(1)/(mH*mH)+HpengMuEEECR2*SRE(2)/(mA*mA))/(4.0*pi*alph);
1456 gslpp::complex B3HiggscL = (HpengMuEEECL0*SLE(0)/(mh*mh)+HpengMuEEECL1*SLE(1)/(mH*mH)+HpengMuEEECL2*SLE(2)/(mA*mA))/(4.0*pi*alph);
1459 BHFunctions.assign(0, B2HiggsnR+B2HiggscR );
1460 BHFunctions.assign(1, B2HiggsnL+B2HiggscL );
1461 BHFunctions.assign(2, B3HiggsnR+B3HiggscR );
1462 BHFunctions.assign(3, B3HiggsnL+B3HiggscL );
1474 for (
int x=0;x<6;x++) {
1475 for (
int a=0;a<4;a++) {
1476 for (
int b=0;b<4;b++) {
1580 for (
int y=0;y<6;y++) {
1613 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))
1622 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))
1631 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))
1640 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))
1649 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))
1658 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))
1670 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);
1671 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);
1672 gslpp::complex B3HiggsnRtm = (HpengTauMUMUMUNR0*SRMU(0)/(mh*mh)+HpengTauMUMUMUNR1*SRMU(1)/(mH*mH)+HpengTauMUMUMUNR2*SRMU(2)/(mA*mA))/(4.0*pi*alph);
1673 gslpp::complex B3HiggsnLtm = (HpengTauMUMUMUNL0*SLMU(0)/(mh*mh)+HpengTauMUMUMUNL1*SLMU(1)/(mH*mH)+HpengTauMUMUMUNL2*SLMU(2)/(mA*mA))/(4.0*pi*alph);
1682 for (
int x=0;x<3;x++) {
1683 for (
int a=0;a<2;a++) {
1684 for (
int b=0;b<2;b++) {
1788 for (
int y=0;y<3;y++) {
1821 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))
1830 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))
1839 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))
1848 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))
1857 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))
1866 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))
1878 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);
1879 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);
1880 gslpp::complex B3HiggscRtm = (HpengTauMUMUMUCR0*SRMU(0)/(mh*mh)+HpengTauMUMUMUCR1*SRMU(1)/(mH*mH)+HpengTauMUMUMUCR2*SRMU(2)/(mA*mA))/(4.0*pi*alph);
1881 gslpp::complex B3HiggscLtm = (HpengTauMUMUMUCL0*SLMU(0)/(mh*mh)+HpengTauMUMUMUCL1*SLMU(1)/(mH*mH)+HpengTauMUMUMUCL2*SLMU(2)/(mA*mA))/(4.0*pi*alph);
1884 BHFunctions.assign(0, B2HiggsnRtm+B2HiggscRtm );
1885 BHFunctions.assign(1, B2HiggsnLtm+B2HiggscLtm );
1886 BHFunctions.assign(2, B3HiggsnRtm+B3HiggscRtm );
1887 BHFunctions.assign(3, B3HiggsnLtm+B3HiggscLtm );
1899 for (
int x=0;x<6;x++) {
1900 for (
int a=0;a<4;a++) {
1901 for (
int b=0;b<4;b++) {
2005 for (
int y=0;y<6;y++) {
2095 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);
2096 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);
2097 gslpp::complex B3HiggsnRte = (HpengTauEEENR0*SRE(0)/(mh*mh)+HpengTauEEENR1*SRE(1)/(mH*mH)+HpengTauEEENR2*SRE(2)/(mA*mA))/(4.0*pi*alph);
2098 gslpp::complex B3HiggsnLte = (HpengTauEEENL0*SLE(0)/(mh*mh)+HpengTauEEENL1*SLE(1)/(mH*mH)+HpengTauEEENL2*SLE(2)/(mA*mA))/(4.0*pi*alph);
2107 for (
int x=0;x<3;x++) {
2108 for (
int a=0;a<2;a++) {
2109 for (
int b=0;b<2;b++) {
2213 for (
int y=0;y<3;y++) {
2246 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))
2255 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))
2264 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))
2273 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))
2282 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))
2291 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))
2303 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);
2304 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);
2305 gslpp::complex B3HiggscRte = (HpengTauEEECR0*SRE(0)/(mh*mh)+HpengTauEEECR1*SRE(1)/(mH*mH)+HpengTauEEECR2*SRE(2)/(mA*mA))/(4.0*pi*alph);
2306 gslpp::complex B3HiggscLte = (HpengTauEEECL0*SLE(0)/(mh*mh)+HpengTauEEECL1*SLE(1)/(mH*mH)+HpengTauEEECL2*SLE(2)/(mA*mA))/(4.0*pi*alph);
2309 BHFunctions.assign(0, B2HiggsnRte+B2HiggscRte );
2310 BHFunctions.assign(1, B2HiggsnLte+B2HiggscLte );
2311 BHFunctions.assign(2, B3HiggsnRte+B3HiggscRte );
2312 BHFunctions.assign(3, B3HiggsnLte+B3HiggscLte );
2315 return(BHFunctions);
2325 double sw2 =
mySUSY.StandardModel::sW2(MW);
2326 double stw =
sqrt(sw2);
2327 double ctw =
sqrt(1.0 - sw2);
2328 double ttw = stw/ctw;
2334 double cdenn = MW*
cosb;
2336 double g2t = g2/
sqrt(2.0);
2341 for (
int a=0;a<4;a++) {
2342 for (
int x=0;x<6;x++) {
2344 NRlE.assign(a, x, - (g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 0) + (mE/cdenn)*
ON(a, 2)*
myRl(x, 3)));
2345 NRlMU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 1) + (mMU/cdenn)*
ON(a, 2)*
myRl(x, 4)));
2347 NRqDOWN.assign(a, x, -(g2t)*((-
ON(a, 1)+
ON(a, 0)*ttw/3.0)*
myRd(x, 0)+mDOWN/(MW*cosb)*
ON(a, 2)*
myRd(x, 3)));
2349 NLlE.assign(a, x, -(g2t)*((mE/cdenn)*
ON(a, 2)*
myRl(x, 0) + 2.0*
ON(a, 0)*ttw*
myRl(x, 3)));
2350 NLlMU.assign(a, x, -(g2t)*((mMU/cdenn)*
ON(a, 2)*
myRl(x, 1) + 2.0*
ON(a, 0)*ttw*
myRl(x, 4)));
2351 NLqUP.assign(a, x, -(g2t)*(mUP/(MW*sinb)*
ON(a, 3)*
myRu(x, 0)-4.0*
ON(a, 0)*
myRu(x, 3)*ttw/3.0));
2352 NLqDOWN.assign(a, x, -(g2t)*(mDOWN/(MW*cosb)*
ON(a, 2)*
myRd(x, 0)+2.0*
ON(a, 0)*
myRd(x, 3)*ttw/3.0));
2356 for (
int a=0;a<2;a++) {
2357 for (
int x=0;x<3;x++) {
2362 CLlE.assign(a, x, g2*mE/cdenc*
myU(a, 1).conjugate()*
myRn(x, 0));
2363 CLlMU.assign(a, x, g2*mMU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 1));
2366 for (
int a=0;a<2;a++) {
2367 for (
int x=0;x<6;x++) {
2382 for (
int a=0;a<4;a++) {
2383 for (
int b=0;b<4;b++) {
2384 for (
int x=0;x<6;x++) {
2385 for (
int y=0;y<6;y++) {
2420 for (
int a=0;a<2;a++) {
2421 for (
int b=0;b<2;b++) {
2422 for (
int x=0;x<3;x++) {
2423 for (
int y=0;y<6;y++) {
2446 DFunctions.assign(0, DunR+DucR );
2447 DFunctions.assign(1, DunL+DucL );
2448 DFunctions.assign(2, DdnR+DdcR );
2449 DFunctions.assign(3, DdnL+DdcL );
2462 double piconst = 1.0/(32.0 * pi * pi);
2463 double sw2 =
mySUSY.StandardModel::sW2(MW);
2464 double stw =
sqrt(sw2);
2465 double ctw =
sqrt(1.0 - sw2);
2466 double ttw = stw/ctw;
2472 double cdenn = MW*
cosb;
2474 double g2t = g2/
sqrt(2.0);
2481 for (
int a=0;a<4;a++) {
2482 for (
int x=0;x<6;x++) {
2484 NRlE.assign(a, x, - (g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 0) + (mE/cdenn)*
ON(a, 2)*
myRl(x, 3)));
2485 NRlMU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 1) + (mMU/cdenn)*
ON(a, 2)*
myRl(x, 4)));
2486 NRlTAU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 2) + (mTAU/cdenn)*
ON(a, 2)*
myRl(x, 5)));
2488 NLlE.assign(a, x, -(g2t)*((mE/cdenn)*
ON(a, 2)*
myRl(x, 0) + 2.0*
ON(a, 0)*ttw*
myRl(x, 3)));
2489 NLlMU.assign(a, x, -(g2t)*((mMU/cdenn)*
ON(a, 2)*
myRl(x, 1) + 2.0*
ON(a, 0)*ttw*
myRl(x, 4)));
2490 NLlTAU.assign(a, x, -(g2t)*((mTAU/cdenn)*
ON(a, 2)*
myRl(x, 2) + 2.0*
ON(a, 0)*ttw*
myRl(x, 5)));
2499 for (
int a=0;a<2;a++) {
2500 for (
int x=0;x<3;x++) {
2506 CLlE.assign(a, x, g2*mE/cdenc*
myU(a, 1).conjugate()*
myRn(x, 0));
2507 CLlMU.assign(a, x, g2*mMU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 1));
2508 CLlTAU.assign(a, x, g2*mTAU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 2));
2514 for (
int a=0;a<4;a++) {
2515 for (
int x=0;x<6;x++) {
2520 for (
int x=0;x<6;x++) {
2521 for (
int a=0;a<4;a++) {
2522 for (
int b=0;b<4;b++) {
2523 if (a != b && std::abs(
Lepty(a,x)-
Lepty(b,x)) > 0.01 && std::abs(1.0-
Lepty(a,x)) > 0.01 && std::abs(1.0-
Lepty(b,x)) > 0.01) {
2529 else if (a != b && std::abs(
Lepty(a,x)-
Lepty(b,x)) > 0.01 && std::abs(1.0-
Lepty(a,x)) > 0.01 && std::abs(1.0-
Lepty(b,x)) <= 0.01) {
2533 else if (a != b && std::abs(
Lepty(a,x)-
Lepty(b,x)) > 0.01 && std::abs(1.0-
Lepty(b,x)) > 0.01 && std::abs(1.0-
Lepty(a,x)) <= 0.01) {
2537 else if ((a == b || std::abs(
Lepty(a,x)-
Lepty(b,x)) <= 0.01) && std::abs(1.0-
Lepty(a,x)) > 0.01 && std::abs(1.0-
Lepty(b,x)) > 0.01) {
2550 for (
int a=0;a<2;a++) {
2551 for (
int x=0;x<3;x++) {
2556 for (
int x=0;x<3;x++) {
2557 for (
int a=0;a<2;a++) {
2558 for (
int b=0;b<2;b++) {
2559 if (a != b && std::abs(
Leptz(a,x)-
Leptz(b,x)) > 0.01 && std::abs(1.0-
Leptz(a,x)) > 0.01 && std::abs(1.0-
Leptz(b,x)) > 0.01) {
2565 else if (a != b && std::abs(
Leptz(a,x)-
Leptz(b,x)) > 0.01 && std::abs(1.0-
Leptz(a,x)) > 0.01 && std::abs(1.0-
Leptz(b,x)) <= 0.01) {
2569 else if (a != b && std::abs(
Leptz(a,x)-
Leptz(b,x)) > 0.01 && std::abs(1.0-
Leptz(b,x)) > 0.01 && std::abs(1.0-
Leptz(a,x)) <= 0.01) {
2573 else if ((a == b || std::abs(
Leptz(a,x)-
Leptz(b,x)) <= 0.01) && std::abs(1.0-
Leptz(a,x)) > 0.01 && std::abs(1.0-
Leptz(b,x)) > 0.01) {
2590 for (
int x=0;x<6;x++) {
2591 for (
int a=0;a<4;a++) {
2592 for (
int b=0;b<4;b++) {
2593 ZpengMuEEENR = ZpengMuEEENR -
NLlE(a,x)*
NLlMU(b,x)*piconst*
2595 ZpengMuEEENL = ZpengMuEEENL +
NRlE(a,x)*
NRlMU(b,x)*piconst*
2602 for (
int x=0;x<3;x++) {
2603 for (
int a=0;a<2;a++) {
2604 for (
int b=0;b<2;b++) {
2605 ZpengMuEEEC = ZpengMuEEEC +
CRlE(a,x)*
CRlMU(b,x)*piconst*
2612 FFunctions.assign(0, ZpengMuEEENR/(MZ*MZ*ctw*ctw) );
2613 FFunctions.assign(1, ZpengMuEEENR*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
2614 FFunctions.assign(2, (ZpengMuEEENL + ZpengMuEEEC)/(MZ*MZ*ctw*ctw) );
2615 FFunctions.assign(3, (ZpengMuEEENL + ZpengMuEEEC)*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
2622 for (
int x=0;x<6;x++) {
2623 for (
int a=0;a<4;a++) {
2624 for (
int b=0;b<4;b++) {
2625 ZpengTauMuMuMuNR = ZpengTauMuMuMuNR -
NLlMU(a,x)*
NLlTAU(b,x)*piconst*
2627 ZpengTauMuMuMuNL = ZpengTauMuMuMuNL +
NRlMU(a,x)*
NRlTAU(b,x)*piconst*
2634 for (
int x=0;x<3;x++) {
2635 for (
int a=0;a<2;a++) {
2636 for (
int b=0;b<2;b++) {
2637 ZpengTauMuMuMuC = ZpengTauMuMuMuC +
CRlMU(a,x)*
CRlTAU(b,x)*piconst*
2644 FFunctions.assign(0, ZpengTauMuMuMuNR/(MZ*MZ*ctw*ctw) );
2645 FFunctions.assign(1, ZpengTauMuMuMuNR*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
2646 FFunctions.assign(2, (ZpengTauMuMuMuNL + ZpengTauMuMuMuC)/(MZ*MZ*ctw*ctw) );
2647 FFunctions.assign(3, (ZpengTauMuMuMuNL + ZpengTauMuMuMuC)*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
2654 for (
int x=0;x<6;x++) {
2655 for (
int a=0;a<4;a++) {
2656 for (
int b=0;b<4;b++) {
2657 ZpengTauEEENR = ZpengTauEEENR -
NLlE(a,x)*
NLlTAU(b,x)*piconst*
2659 ZpengTauEEENL = ZpengTauEEENL +
NRlE(a,x)*
NRlTAU(b,x)*piconst*
2666 for (
int x=0;x<3;x++) {
2667 for (
int a=0;a<2;a++) {
2668 for (
int b=0;b<2;b++) {
2669 ZpengTauEEEC = ZpengTauEEEC +
CRlE(a,x)*
CRlTAU(b,x)*piconst*
2676 FFunctions.assign(0, ZpengTauEEENR/(MZ*MZ*ctw*ctw) );
2677 FFunctions.assign(1, ZpengTauEEENR*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
2678 FFunctions.assign(2, (ZpengTauEEENL + ZpengTauEEEC)/(MZ*MZ*ctw*ctw) );
2679 FFunctions.assign(3, (ZpengTauEEENL + ZpengTauEEEC)*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
2686 for (
int x=0;x<6;x++) {
2687 for (
int a=0;a<4;a++) {
2688 for (
int b=0;b<4;b++) {
2689 ZpengTauMuEENR = ZpengTauMuEENR -
NLlMU(a,x)*
NLlTAU(b,x)*piconst*
2691 ZpengTauMuEENL = ZpengTauMuEENL +
NRlMU(a,x)*
NRlTAU(b,x)*piconst*
2698 for (
int x=0;x<3;x++) {
2699 for (
int a=0;a<2;a++) {
2700 for (
int b=0;b<2;b++) {
2701 ZpengTauMuEEC = ZpengTauMuEEC +
CRlMU(a,x)*
CRlTAU(b,x)*piconst*
2708 FFunctions.assign(0, ZpengTauMuEENR/(MZ*MZ*ctw*ctw) );
2709 FFunctions.assign(1, ZpengTauMuEENR*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
2710 FFunctions.assign(2, (ZpengTauMuEENL + ZpengTauMuEEC)/(MZ*MZ*ctw*ctw) );
2711 FFunctions.assign(3, (ZpengTauMuEENL + ZpengTauMuEEC)*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
2723 double piconst = 1.0/(32.0 * pi * pi);
2724 double sw2 =
mySUSY.StandardModel::sW2(MW);
2725 double stw =
sqrt(sw2);
2726 double ctw =
sqrt(1.0 - sw2);
2727 double ttw = stw/ctw;
2731 double cdenn = MW*
cosb;
2733 double g2t = g2/
sqrt(2.0);
2738 for (
int a=0;a<4;a++) {
2739 for (
int x=0;x<6;x++) {
2741 NRlMU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 1) + (mMU/cdenn)*
ON(a, 2)*
myRl(x, 4)));
2743 NLlMU.assign(a, x, -(g2t)*((mMU/cdenn)*
ON(a, 2)*
myRl(x, 1) + 2.0*
ON(a, 0)*ttw*
myRl(x, 4)));
2750 for (
int a=0;a<2;a++) {
2751 for (
int x=0;x<3;x++) {
2755 CLlMU.assign(a, x, g2*mMU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 1));
2759 for (
int a=0;a<4;a++) {
2760 for (
int x=0;x<6;x++) {
2765 for (
int a=0;a<2;a++) {
2766 for (
int x=0;x<3;x++) {
2771 for (
int a=0;a<4;a++) {
2772 for (
int x=0;x<6;x++) {
2773 if (fabs(1.0 -
Lepty(a, x)) > 0.01) {
2776 (6.0 *
pow((1.0 -
Lepty(a,x)),4.0)) );
2787 for (
int a=0;a<2;a++) {
2788 for (
int x=0;x<3;x++) {
2789 if(fabs(1.0-
Leptz(a, x)) > 0.01) {
2792 (6.0*
pow((1.0 -
Leptz(a, x)),4.0))) );
2807 for (
int a=0;a<4;a++) {
2808 for (
int x=0;x<6;x++) {
2809 g2ARN = g2ARN -mMU*mMU*piconst*(4.0*
NRlMU(a,x)*
NRlMU(a,x).conjugate()*
Leptf1(a,x)
2811 g2ALN = g2ALN -mMU*mMU*piconst*(4.0*
NLlMU(a,x)*
NLlMU(a,x).conjugate()*
Leptf1(a,x)
2819 for (
int a=0;a<2;a++) {
2820 for (
int x=0;x<3;x++) {
2821 g2ARC = g2ARC +mMU*mMU*piconst*(4.0*
CRlMU(a,x)*
CRlMU(a,x).conjugate()*
Leptf3(a,x)
2823 g2ALC = g2ALC +mMU*mMU*piconst*(4.0*
CLlMU(a,x)*
CLlMU(a,x).conjugate()*
Leptf3(a,x)
2829 gminus2mu.assign(0, g2ARN + g2ARC );
2830 gminus2mu.assign(1, g2ALN + g2ALC );
2849 double sw2 =
mySUSY.StandardModel::sW2(MW);
2850 double cw2 = 1.0 - sw2;
2856 + (1.0-1.0/(4.0*sw2))/(cw2*MZ*MZ)*
FFunctions(li_to_lj)(1)
2859 + (1.0-1.0/(4.0*sw2))/(cw2*MZ*MZ)*
FFunctions(li_to_lj)(0)
2869 double sw2 =
mySUSY.StandardModel::sW2(MW);
2870 double cw2 = 1.0 - sw2;
2875 C10.assign(0, 1.0/(4.0*sw2*cw2*MZ*MZ)*
FFunctions(li_to_lj)(1)
2877 C10.assign(1, 1.0/(4.0*sw2*cw2*MZ*MZ)*
FFunctions(li_to_lj)(0)
2929 std::stringstream out;
2931 throw std::runtime_error(
"SUSYMatching::CMDLij(): order " + out.str() +
" not implemented.\nFor lepton flavour violating observables only Leading Order (LO) necessary.");
2985 std::stringstream out;
2987 throw std::runtime_error(
"SUSYMatching::CMDLi3j(): order " + out.str() +
" not implemented.\nFor lepton flavour violating observables only Leading Order (LO) necessary.");
2998 double sw2 =
mySUSY.StandardModel::sW2(MW);
3019 std::stringstream out;
3021 throw std::runtime_error(
"SUSYMatching::CMmueconv(): order " + out.str() +
" not implemented.\nFor lepton flavour violating observables only Leading Order (LO) necessary.");
3042 std::stringstream out;
3044 throw std::runtime_error(
"SUSYMatching::CMgminus2mu(): order " + out.str() +
" not implemented.\nFor lepton flavour violating observables only Leading Order (LO) necessary.");
virtual gslpp::vector< gslpp::complex > FFunctions(int n)
Calculates Z penguin amplitudes for the process from .
WilsonCoefficient mcgminus2mu
virtual gslpp::vector< gslpp::complex > gminus2mu()
Calculates amplitudes for at one-loop from .
std::complex< double > Leptfzn[6][4][4]
gslpp::matrix< gslpp::complex > AmpTauA1LN
gslpp::matrix< gslpp::complex > myRl
gslpp::matrix< gslpp::complex > AmpTauALN
WilsonCoefficient mcmueconv
gslpp::complex C0(const double p2, const double m02, const double m12, const double m22) const
.
virtual std::vector< WilsonCoefficient > & CMDLi3j(int li_lj)
gslpp::matrix< gslpp::complex > AmpTEARN
virtual std::vector< WilsonCoefficient > & CMmueconv()
Wilson coefficient for the process conversion in Nuclei.
gslpp::vector< double > mym_se_sq
gslpp::matrix< gslpp::complex > NLlE
gslpp::vector< double > MChi
Chargino mass-eigenvalue.
gslpp::vector< double > mym_sn_sq
gslpp::matrix< gslpp::complex > getN() const
Gets the rotation matrix for neutralinos.
virtual gslpp::vector< gslpp::complex > C7_Lepton(int n)
Calculates and Wilson coefficients for the process .
Particle getLeptons(const StandardModel::lepton p) const
A get method to retrieve the member object of a lepton.
gslpp::vector< double > mym_sd_sq
virtual double Mw_tree() const
The tree-level mass of the boson, .
gslpp::matrix< gslpp::complex > NLqUP
std::complex< double > Leptfzc[3][2][2]
gslpp::matrix< double > Lepty
Functions needed to calculate various LFV observables.
gslpp::matrix< gslpp::complex > ON
gslpp::matrix< double > Leptf4
complex conjugate() const
virtual std::vector< WilsonCoefficient > & CMDLi3j(int li_lj)
Wilson coefficient for the process .
gslpp::matrix< gslpp::complex > CRlTAU
virtual gslpp::vector< gslpp::complex > C10_Lepton(int n)
Calculates and Wilson coefficients for the process .
gslpp::matrix< double > Leptfa1
complex pow(const complex &z1, const complex &z2)
virtual double getMHl() const
Gets the light Higgs mass.
gslpp::matrix< double > Leptf2
gslpp::matrix< gslpp::complex > AmpTEA1LC
virtual gslpp::vector< gslpp::complex > DFunctions()
Calculates box diagram amplitudes for the process conversion in Nuclei from from ...
double getMHa() const
Gets the pseudo-scalar Higgs mass.
virtual gslpp::vector< gslpp::complex > BHFunctions(int n)
Calculates Higgs penguin amplitudes for the process from .
gslpp::matrix< gslpp::complex > CLlTAU
gslpp::matrix< gslpp::complex > AmpA1LC
virtual std::vector< WilsonCoefficient > & CMmueconv()
gslpp::matrix< gslpp::complex > CLqDOWN
gslpp::matrix< gslpp::complex > AmpTauALC
gslpp::matrix< gslpp::complex > NRlMU
std::vector< WilsonCoefficient > vmcgminus2mu
static const complex & i()
int delta_ab(int a, int b)
Kronecker delta.
gslpp::matrix< gslpp::complex > AmpA1LN
Amplitudes of Chargino and Neutralino contribution to various LFV observables.
gslpp::vector< double > MNeig
gslpp::matrix< gslpp::complex > myRd
gslpp::matrix< gslpp::complex > CLqUP
gslpp::matrix< gslpp::complex > AmpTauA1RN
gslpp::vector< double > getMch() const
Gets the chargino masses.
virtual gslpp::vector< gslpp::complex > BFunctions(int n)
Calculates box diagram amplitudes for the process from from .
gslpp::matrix< double > Leptz
gslpp::matrix< gslpp::complex > AmpTEALN
gslpp::matrix< gslpp::complex > AmpARN
virtual gslpp::vector< gslpp::complex > CS_Lepton(int n)
Calculates and Wilson coefficients for the process .
gslpp::complex getSaeff() const
Gets the sine of the effective mixing angle for the CP-even neutral Higgs bosons. ...
void setCoeff(const gslpp::vector< gslpp::complex > &z, orders order_i)
void assign(const size_t &i, const size_t &j, const double &a)
gslpp::matrix< gslpp::complex > AmpA1RN
A class for the matching in the Standard Model.
virtual std::vector< WilsonCoefficient > & CMgminus2mu()
Wilson coefficient for at one-loop.
gslpp::matrix< gslpp::complex > CRlMU
double getMHh() const
Gets the heavy Higgs mass.
gslpp::matrix< gslpp::complex > myRn
virtual std::vector< WilsonCoefficient > & CMDLij(int li_lj)
Wilson coefficient for the process .
gslpp::matrix< gslpp::complex > AmpTEARC
gslpp::matrix< gslpp::complex > getTEhat() const
Gets the trilinear-coupling matrix for charged sleptons.
gslpp::matrix< gslpp::complex > NRlE
gslpp::matrix< gslpp::complex > NRqDOWN
gslpp::matrix< gslpp::complex > myRu
Sfermion mixing matrices.
gslpp::matrix< gslpp::complex > AmpTEA1RC
double getTanb() const
Gets .
WilsonCoefficient mcDLi3j
gslpp::matrix< gslpp::complex > myV
Chargino mixing matrices.
gslpp::matrix< double > Leptf3
gslpp::matrix< gslpp::complex > AmpTauARN
virtual double v() const
The Higgs vacuum expectation value. where is the Fermi constant, measured through muon decays...
std::vector< WilsonCoefficient > vmcmueconv
gslpp::vector< double > getMsn2() const
Gets the sneutrino mass squared.
gslpp::matrix< gslpp::complex > AmpARC
double getGF() const
A get method to retrieve the Fermi constant .
gslpp::matrix< gslpp::complex > NRqUP
gslpp::matrix< gslpp::complex > CRlE
Chargino and Neutralino couplings to sfermions.
void NeutralinoRemixing()
gslpp::matrix< gslpp::complex > AmpTEA1RN
gslpp::vector< double > mym_su_sq
Sfermion mass-eigenvalue squared.
gslpp::complex B1(const double mu2, const double p2, const double m02, const double m12) const
.
gslpp::matrix< gslpp::complex > CRqDOWN
gslpp::vector< double > getMse2() const
Gets the charged slepton mass squared.
gslpp::matrix< gslpp::complex > CLlMU
gslpp::complex D0(const double s, const double t, const double m02, const double m12, const double m22, const double m32) const
.
gslpp::matrix< gslpp::complex > getRu() const
Gets the rotation matrix for up-type squarks.
gslpp::matrix< gslpp::complex > CRqUP
double C11(const double m12, const double m22, const double m32) const
.
gslpp::complex B0(const double mu2, const double p2, const double m02, const double m12) const
.
Particle getQuarks(const quark q) const
A get method to access a quark as an object of the type Particle.
gslpp::matrix< gslpp::complex > NLlTAU
double C12(const double m12, const double m22, const double m32) const
.
gslpp::matrix< gslpp::complex > NRlTAU
gslpp::matrix< gslpp::complex > TEhat
Slepton tri-linear coupling matrix.
virtual gslpp::vector< gslpp::complex > C9_Lepton(int n)
Calculates and Wilson coefficients for the process .
virtual std::vector< WilsonCoefficient > & CMgminus2mu()
double getSinb() const
Gets .
double getAle() const
A get method to retrieve the fine-structure constant .
virtual gslpp::vector< gslpp::complex > CP_Lepton(int n)
Calculates and Wilson coefficients for the process .
gslpp::matrix< double > Leptf1
double getCosb() const
Gets .
gslpp::matrix< double > Leptfa2
virtual gslpp::vector< gslpp::complex > AFunctions(int n)
Calculates gamma penguin amplitudes for the process from .
const double & getMass() const
A get method to access the particle mass.
gslpp::matrix< gslpp::complex > getRn() const
Gets the rotation matrix for sneutrinos.
gslpp::complex D00(const double s, const double t, const double m02, const double m12, const double m22, const double m32) const
.
gslpp::matrix< gslpp::complex > AmpTEALC
gslpp::matrix< gslpp::complex > AmpA1RC
gslpp::matrix< gslpp::complex > AmpTauA1RC
virtual gslpp::vector< gslpp::complex > CT_Lepton(int n)
Calculates and Wilson coefficients for the process .
gslpp::matrix< gslpp::complex > CLlE
std::vector< WilsonCoefficient > vmcDLij
gslpp::complex getM2() const
Gets the wino mass.
complex log(const complex &z)
gslpp::matrix< gslpp::complex > getRd() const
Gets the rotation matrix for down-type squarks.
gslpp::vector< double > MChi0
Neutralino mass-eigenvalue.
gslpp::matrix< gslpp::complex > NLqDOWN
gslpp::matrix< gslpp::complex > getV() const
Gets the rotation matrix for positive charginos.
std::vector< WilsonCoefficient > vmcDLi3j
gslpp::matrix< gslpp::complex > getRl() const
Gets the rotation matrix for charged sleptons.
gslpp::complex getMuH() const
Gets the parameter in the superpotential.
gslpp::matrix< gslpp::complex > getU() const
Gets the rotation matrix for negative charginos.
gslpp::complex getM1() const
Gets the bino mass.
A class for defining operations on and functions of complex numbers.
gslpp::matrix< gslpp::complex > NLlMU
gslpp::matrix< gslpp::complex > AmpALN
A base class for SUSY models.
gslpp::matrix< gslpp::complex > AmpALC
gslpp::matrix< gslpp::complex > AmpTauARC
gslpp::vector< double > getMneu() const
Gets the neutralino masses.
std::complex< double > Leptgzn[6][4][4]
gslpp::matrix< gslpp::complex > myU
SUSYMatching(const SUSY &SUSY_i)
void updateSUSYParameters()
Updates to new SUSY parameter sets.
double getMz() const
A get method to access the mass of the boson .
virtual std::vector< WilsonCoefficient > & CMDLij(int li_lj)
gslpp::matrix< gslpp::complex > myN
Neutralino mixing matrix.
gslpp::matrix< gslpp::complex > AmpTauA1LC
std::complex< double > Leptgzc[3][2][2]
complex sqrt(const complex &z)
gslpp::matrix< gslpp::complex > AmpTEA1LN