11 #include <gsl/gsl_integration.h>
23 mcmueconv(8,
NDR,
LO),
24 mcgminus2mu(2,
NDR,
LO),
153 for (
int a=0;a<4;a++) {
154 Mdiag.assign(a,a,
MChi0(a));
156 MN_tmp =
myN.transpose() * Mdiag *
myN;
158 MN_tmp.eigensystem(cmplxONT,
MNeig);
159 for (
int a=0;a<4;a++) {
160 for (
int b=0;b<4;b++) {
161 ON.assign(a,b,cmplxONT(b,a).real());
171 double MW =
mySUSY.Mw_tree();
173 double piconst = 1.0/(32.0 * pi * pi);
174 double sw2 =
mySUSY.StandardModel::sW2(MW);
175 double stw =
sqrt(sw2);
176 double ctw =
sqrt(1.0 - sw2);
177 double ttw = stw/ctw;
183 double cdenn = MW*
cosb;
185 double g2t = g2/
sqrt(2.0);
190 for (
int a=0;a<4;a++) {
191 for (
int x=0;x<6;x++) {
193 NRlE.assign(a, x, - (g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 0) + (mE/cdenn)*
ON(a, 2)*
myRl(x, 3)));
194 NRlMU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 1) + (mMU/cdenn)*
ON(a, 2)*
myRl(x, 4)));
195 NRlTAU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 2) + (mTAU/cdenn)*
ON(a, 2)*
myRl(x, 5)));
197 NLlE.assign(a, x, -(g2t)*((mE/cdenn)*
ON(a, 2)*
myRl(x, 0) + 2.0*
ON(a, 0)*ttw*
myRl(x, 3)));
198 NLlMU.assign(a, x, -(g2t)*((mMU/cdenn)*
ON(a, 2)*
myRl(x, 1) + 2.0*
ON(a, 0)*ttw*
myRl(x, 4)));
199 NLlTAU.assign(a, x, -(g2t)*((mTAU/cdenn)*
ON(a, 2)*
myRl(x, 2) + 2.0*
ON(a, 0)*ttw*
myRl(x, 5)));
208 for (
int a=0;a<2;a++) {
209 for (
int x=0;x<3;x++) {
215 CLlE.assign(a, x, g2*mE/cdenc*
myU(a, 1).conjugate()*
myRn(x, 0));
216 CLlMU.assign(a, x, g2*mMU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 1));
217 CLlTAU.assign(a, x, g2*mTAU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 2));
222 for (
int a=0;a<4;a++) {
223 for (
int x=0;x<6;x++) {
228 for (
int a=0;a<2;a++) {
229 for (
int x=0;x<3;x++) {
234 for (
int a=0;a<4;a++) {
235 for (
int x=0;x<6;x++) {
236 if (fabs(1.0 -
Lepty(a, x)) > 0.01) {
242 (6.0 *
pow((1.0 -
Lepty(a,x)),4.0)) );
254 for (
int a=0;a<2;a++) {
255 for (
int x=0;x<3;x++) {
256 if(fabs(1.0-
Leptz(a, x)) > 0.01) {
262 (6.0*
pow((1.0 -
Leptz(a, x)),4.0))) );
278 for (
int a=0;a<4;a++) {
279 for (
int x=0;x<6;x++) {
298 for (
int a=0;a<4;a++) {
299 for (
int x=0;x<6;x++) {
307 for (
int a=0;a<2;a++) {
308 for (
int x=0;x<3;x++) {
327 for (
int a=0;a<2;a++) {
328 for (
int x=0;x<3;x++) {
346 for (
int a=0;a<4;a++) {
347 for (
int x=0;x<6;x++) {
366 for (
int a=0;a<4;a++) {
367 for (
int x=0;x<6;x++) {
375 for (
int a=0;a<2;a++) {
376 for (
int x=0;x<3;x++) {
395 for (
int a=0;a<2;a++) {
396 for (
int x=0;x<3;x++) {
413 for (
int a=0;a<4;a++) {
414 for (
int x=0;x<6;x++) {
433 for (
int a=0;a<4;a++) {
434 for (
int x=0;x<6;x++) {
442 for (
int a=0;a<2;a++) {
443 for (
int x=0;x<3;x++) {
462 for (
int a=0;a<2;a++) {
463 for (
int x=0;x<3;x++) {
486 double MW =
mySUSY.Mw_tree();
488 double sw2 =
mySUSY.StandardModel::sW2(MW);
489 double stw =
sqrt(sw2);
490 double ctw =
sqrt(1.0 - sw2);
491 double ttw = stw/ctw;
497 double cdenn = MW*
cosb;
499 double g2t = g2/
sqrt(2.0);
500 double alph =
mySUSY.getAle();
505 for (
int a=0;a<4;a++) {
506 for (
int x=0;x<6;x++) {
508 NRlE.assign(a, x, - (g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 0) + (mE/cdenn)*
ON(a, 2)*
myRl(x, 3)));
509 NRlMU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 1) + (mMU/cdenn)*
ON(a, 2)*
myRl(x, 4)));
510 NRlTAU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 2) + (mTAU/cdenn)*
ON(a, 2)*
myRl(x, 5)));
512 NLlE.assign(a, x, -(g2t)*((mE/cdenn)*
ON(a, 2)*
myRl(x, 0) + 2.0*
ON(a, 0)*ttw*
myRl(x, 3)));
513 NLlMU.assign(a, x, -(g2t)*((mMU/cdenn)*
ON(a, 2)*
myRl(x, 1) + 2.0*
ON(a, 0)*ttw*
myRl(x, 4)));
514 NLlTAU.assign(a, x, -(g2t)*((mTAU/cdenn)*
ON(a, 2)*
myRl(x, 2) + 2.0*
ON(a, 0)*ttw*
myRl(x, 5)));
523 for (
int a=0;a<2;a++) {
524 for (
int x=0;x<3;x++) {
530 CLlE.assign(a, x, g2*mE/cdenc*
myU(a, 1).conjugate()*
myRn(x, 0));
531 CLlMU.assign(a, x, g2*mMU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 1));
532 CLlTAU.assign(a, x, g2*mTAU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 2));
551 for (
int a=0;a<4;a++) {
552 for (
int b=0;b<4;b++) {
553 for (
int x=0;x<6;x++) {
554 for (
int t=0;t<6;t++) {
555 B1nRMu3E = B1nRMu3E + (1.0/(4.0*pi*alph))*(0.5*
NLlMU(a,x)*
NLlE(a,t)*
NLlE(b,t)*
NLlE(b,x)
559 B2nRMu3E = B2nRMu3E + (1.0/(4.0*pi*alph))*(0.25*(
NLlMU(a,x)*
NLlE(a,t)*
NRlE(b,t)*
NRlE(b,x)
570 B1nLMu3E = B1nLMu3E + (1.0/(4.0*pi*alph))*(0.5*
NRlMU(a,x)*
NRlE(a,t)*
NRlE(b,t)*
NRlE(b,x)
574 B2nLMu3E = B2nLMu3E + (1.0/(4.0*pi*alph))*(0.25*(
NRlMU(a,x)*
NRlE(a,t)*
NLlE(b,t)*
NLlE(b,x)
600 for (
int a=0;a<2;a++) {
601 for (
int b=0;b<2;b++) {
602 for (
int x=0;x<3;x++) {
603 for (
int t=0;t<3;t++) {
604 B1cRMu3E = B1cRMu3E + (1.0/(8.0*pi*alph))*
CLlMU(a,x)*
CLlE(a,t)*
CLlE(b,t)*
CLlE(b,x)
606 B2cRMu3E = B2cRMu3E + (1.0/(4.0*pi*alph))*(0.25*
CLlMU(a,x)*
CLlE(a,t)*
CRlE(b,t)*
CRlE(b,x)
612 B1cLMu3E = B1cLMu3E + (1.0/(8.0*pi*alph))*(
CRlMU(a,x)*
CRlE(a,t)*
CRlE(b,t)*
CRlE(b,x)
614 B2cLMu3E = B2cLMu3E + (1.0/(4.0*pi*alph))*(0.25*
CRlMU(a,x)*
CRlE(a,t)*
CLlE(b,t)*
CLlE(b,x)
646 for (
int a=0;a<4;a++) {
647 for (
int b=0;b<4;b++) {
648 for (
int x=0;x<6;x++) {
649 for (
int t=0;t<6;t++) {
691 for (
int a=0;a<2;a++) {
692 for (
int b=0;b<2;b++) {
693 for (
int x=0;x<3;x++) {
694 for (
int t=0;t<3;t++) {
716 BFunctions.assign(0, B1nRTau3Mu + B1cRTau3Mu );
717 BFunctions.assign(1, B1nLTau3Mu + B1cLTau3Mu );
718 BFunctions.assign(2, B2nRTau3Mu + B2cRTau3Mu );
719 BFunctions.assign(3, B2nLTau3Mu + B2cLTau3Mu );
720 BFunctions.assign(4, B3nRTau3Mu + B3cRTau3Mu );
721 BFunctions.assign(5, B3nLTau3Mu + B3cLTau3Mu );
736 for (
int a=0;a<4;a++) {
737 for (
int b=0;b<4;b++) {
738 for (
int x=0;x<6;x++) {
739 for (
int t=0;t<6;t++) {
740 B1nRTau3E = B1nRTau3E + (1.0/(4.0*pi*alph))*(0.5*
NLlTAU(a,x)*
NLlE(a,t)*
NLlE(b,t)*
NLlE(b,x)
744 B2nRTau3E = B2nRTau3E + (1.0/(4.0*pi*alph))*(0.25*(
NLlTAU(a,x)*
NLlE(a,t)*
NRlE(b,t)*
NRlE(b,x)
755 B1nLTau3E = B1nLTau3E + (1.0/(4.0*pi*alph))*(0.5*
NRlTAU(a,x)*
NRlE(a,t)*
NRlE(b,t)*
NRlE(b,x)
759 B2nLTau3E = B2nLTau3E + (1.0/(4.0*pi*alph))*(0.25*(
NRlTAU(a,x)*
NRlE(a,t)*
NLlE(b,t)*
NLlE(b,x)
781 for (
int a=0;a<2;a++) {
782 for (
int b=0;b<2;b++) {
783 for (
int x=0;x<3;x++) {
784 for (
int t=0;t<3;t++) {
787 B2cRTau3E = B2cRTau3E + (1.0/(4.0*pi*alph))*(0.25*
CLlTAU(a,x)*
CLlE(a,t)*
CRlE(b,t)*
CRlE(b,x)
795 B2cLTau3E = B2cLTau3E + (1.0/(4.0*pi*alph))*(0.25*
CRlTAU(a,x)*
CRlE(a,t)*
CLlE(b,t)*
CLlE(b,x)
834 double MZ =
mySUSY.getMz();
835 double MW =
mySUSY.Mw_tree();
837 double piconst = 1.0/(32.0 * pi * pi);
838 double sw2 =
mySUSY.StandardModel::sW2(MW);
839 double stw =
sqrt(sw2);
840 double ctw =
sqrt(1.0 - sw2);
841 double ttw = stw/ctw;
852 double mh =
mySUSY.getMHl();
853 double mH =
mySUSY.getMHh();
854 double mA =
mySUSY.getMHa();
861 double cdenn = MW*
cosb;
863 double g2t = g2/
sqrt(2.0);
864 double alph =
mySUSY.getAle();
870 for (
int a=0;a<4;a++) {
871 for (
int x=0;x<6;x++) {
873 NRlE.assign(a, x, - (g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 0) + (mE/cdenn)*
ON(a, 2)*
myRl(x, 3)));
874 NRlMU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 1) + (mMU/cdenn)*
ON(a, 2)*
myRl(x, 4)));
875 NRlTAU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 2) + (mTAU/cdenn)*
ON(a, 2)*
myRl(x, 5)));
877 NLlE.assign(a, x, -(g2t)*((mE/cdenn)*
ON(a, 2)*
myRl(x, 0) + 2.0*
ON(a, 0)*ttw*
myRl(x, 3)));
878 NLlMU.assign(a, x, -(g2t)*((mMU/cdenn)*
ON(a, 2)*
myRl(x, 1) + 2.0*
ON(a, 0)*ttw*
myRl(x, 4)));
879 NLlTAU.assign(a, x, -(g2t)*((mTAU/cdenn)*
ON(a, 2)*
myRl(x, 2) + 2.0*
ON(a, 0)*ttw*
myRl(x, 5)));
888 for (
int a=0;a<2;a++) {
889 for (
int x=0;x<3;x++) {
895 CLlE.assign(a, x, g2*mE/cdenc*
myU(a, 1).conjugate()*
myRn(x, 0));
896 CLlMU.assign(a, x, g2*mMU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 1));
897 CLlTAU.assign(a, x, g2*mTAU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 2));
906 sigma1.assign(0, sina);
907 sigma1.assign(1, -cosa);
909 sigma2.assign(0, cosa);
910 sigma2.assign(1, sina);
912 sigma3.assign(0, sinapb);
913 sigma3.assign(1, -cosapb);
914 sigma3.assign(2, 0.);
915 sigma4.assign(0, -sina);
916 sigma4.assign(1, cosa);
917 sigma4.assign(2, 0.);
918 sigma5.assign(0, -cosbma);
919 sigma5.assign(1, sinbma);
923 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.);
924 for (
int a=0;a<4;a++) {
925 for (
int b=0;b<4;b++) {
926 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))) );
927 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) );
928 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)) );
929 DR0.assign(b, a, DL0(b,a).conjugate() );
930 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)) );
931 DR1.assign(b, a, DL1(b,a).conjugate() );
932 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)) );
933 DR2.assign(b, a, DL2(b,a).conjugate() );
938 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.);
939 for (
int a=0;a<2;a++) {
940 for (
int b=0;b<2;b++) {
941 Qch.assign(a, b,
myU(a,1)*
myV(b,0)/
sqrt(2.0) );
943 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)) );
944 WL0.assign(b, a, WR0(a,b).conjugate() );
945 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)) );
946 WL1.assign(b, a, WR1(a,b).conjugate() );
947 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)) );
948 WL2.assign(b, a, WR2(a,b).conjugate() );
956 for (
int p=0;p<3;p++) {
957 gLLE.assign(p, MZ/ctw*sigma3(p)*(0.5-sw2) + mE*mE/(MW*
cosb)*sigma4(p));
958 gLLMU.assign(p, MZ/ctw*sigma3(p)*(0.5-sw2) + mMU*mMU/(MW*
cosb)*sigma4(p));
959 gLLTAU.assign(p, MZ/ctw*sigma3(p)*(0.5-sw2) + mTAU*mTAU/(MW*
cosb)*sigma4(p));
960 gRRE.assign(p, MZ/ctw*sigma3(p)*sw2 + mE*mE/(MW*
cosb)*sigma4(p));
961 gRRMU.assign(p, MZ/ctw*sigma3(p)*sw2 + mMU*mMU/(MW*
cosb)*sigma4(p));
962 gRRTAU.assign(p, MZ/ctw*sigma3(p)*sw2 + mTAU*mTAU/(MW*
cosb)*sigma4(p));
963 gLRE.assign(p, (-sigma1(p)*
TEhat(0,0)/mE*
v1/
sqrt(2.0)-sigma2(p).conjugate()*muH)*mE/(2.0*MW*
cosb));
964 gLRMU.assign(p, (-sigma1(p)*
TEhat(1,1)/mMU*
v1/
sqrt(2.0)-sigma2(p).conjugate()*muH)*mMU/(2.0*MW*
cosb));
965 gLRTAU.assign(p, (-sigma1(p)*
TEhat(2,2)/mTAU*
v1/
sqrt(2.0)-sigma2(p).conjugate()*muH)*mTAU/(2.0*MW*
cosb));
966 gRLE.assign(p, gLRE(p).conjugate());
967 gRLMU.assign(p, gLRMU(p).conjugate());
968 gRLTAU.assign(p, gLRTAU(p).conjugate());
969 gLLNU.assign(p, -0.5*MZ/ctw*sigma3(p));
986 for (
int x=0;x<6;x++) {
987 for (
int y=0;y<6;y++) {
988 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)
989 +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)
990 +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)));
991 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)
992 +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)
993 +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)));
994 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)
995 +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)
996 +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)));
1001 for (
int x=0;x<3;x++) {
1002 Gnu0.assign(x, x, -g2*gLLNU(0) );
1003 Gnu1.assign(x, x, -g2*gLLNU(1) );
1004 Gnu2.assign(x, x, -g2*gLLNU(2) );
1007 gslpp::vector<gslpp::complex> SRE(3, 0.), SLE(3, 0.), SRMU(3, 0.), SLMU(3, 0.), SRTAU(3, 0.), SLTAU(3, 0.);
1008 for (
int p=0;p<3;p++) {
1009 SRE.assign(p, g2*mE/(2.0*MW*
cosb) * sigma1(p));
1010 SLE.assign(p, g2*mE/(2.0*MW*
cosb) * sigma1(p).conjugate());
1011 SRMU.assign(p, g2*mMU/(2.0*MW*
cosb) * sigma1(p));
1012 SLMU.assign(p, g2*mMU/(2.0*MW*
cosb) * sigma1(p).conjugate());
1013 SRTAU.assign(p, g2*mTAU/(2.0*MW*
cosb) * sigma1(p));
1014 SLTAU.assign(p, g2*mTAU/(2.0*MW*
cosb) * sigma1(p).conjugate());
1026 for (
int x=0;x<6;x++) {
1027 for (
int a=0;a<4;a++) {
1028 for (
int b=0;b<4;b++) {
1132 for (
int y=0;y<6;y++) {
1222 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);
1223 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);
1224 gslpp::complex B3HiggsnR = (HpengMuEEENR0*SRE(0)/(mh*mh)+HpengMuEEENR1*SRE(1)/(mH*mH)+HpengMuEEENR2*SRE(2)/(mA*mA))/(4.0*pi*alph);
1225 gslpp::complex B3HiggsnL = (HpengMuEEENL0*SLE(0)/(mh*mh)+HpengMuEEENL1*SLE(1)/(mH*mH)+HpengMuEEENL2*SLE(2)/(mA*mA))/(4.0*pi*alph);
1234 for (
int x=0;x<3;x++) {
1235 for (
int a=0;a<2;a++) {
1236 for (
int b=0;b<2;b++) {
1340 for (
int y=0;y<3;y++) {
1373 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))
1382 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))
1391 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))
1400 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))
1409 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))
1418 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))
1430 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);
1431 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);
1432 gslpp::complex B3HiggscR = (HpengMuEEECR0*SRE(0)/(mh*mh)+HpengMuEEECR1*SRE(1)/(mH*mH)+HpengMuEEECR2*SRE(2)/(mA*mA))/(4.0*pi*alph);
1433 gslpp::complex B3HiggscL = (HpengMuEEECL0*SLE(0)/(mh*mh)+HpengMuEEECL1*SLE(1)/(mH*mH)+HpengMuEEECL2*SLE(2)/(mA*mA))/(4.0*pi*alph);
1451 for (
int x=0;x<6;x++) {
1452 for (
int a=0;a<4;a++) {
1453 for (
int b=0;b<4;b++) {
1557 for (
int y=0;y<6;y++) {
1590 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))
1599 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))
1608 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))
1617 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))
1626 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))
1635 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))
1647 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);
1648 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);
1649 gslpp::complex B3HiggsnRtm = (HpengTauMUMUMUNR0*SRMU(0)/(mh*mh)+HpengTauMUMUMUNR1*SRMU(1)/(mH*mH)+HpengTauMUMUMUNR2*SRMU(2)/(mA*mA))/(4.0*pi*alph);
1650 gslpp::complex B3HiggsnLtm = (HpengTauMUMUMUNL0*SLMU(0)/(mh*mh)+HpengTauMUMUMUNL1*SLMU(1)/(mH*mH)+HpengTauMUMUMUNL2*SLMU(2)/(mA*mA))/(4.0*pi*alph);
1659 for (
int x=0;x<3;x++) {
1660 for (
int a=0;a<2;a++) {
1661 for (
int b=0;b<2;b++) {
1765 for (
int y=0;y<3;y++) {
1798 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))
1807 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))
1816 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))
1825 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))
1834 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))
1843 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))
1855 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);
1856 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);
1857 gslpp::complex B3HiggscRtm = (HpengTauMUMUMUCR0*SRMU(0)/(mh*mh)+HpengTauMUMUMUCR1*SRMU(1)/(mH*mH)+HpengTauMUMUMUCR2*SRMU(2)/(mA*mA))/(4.0*pi*alph);
1858 gslpp::complex B3HiggscLtm = (HpengTauMUMUMUCL0*SLMU(0)/(mh*mh)+HpengTauMUMUMUCL1*SLMU(1)/(mH*mH)+HpengTauMUMUMUCL2*SLMU(2)/(mA*mA))/(4.0*pi*alph);
1876 for (
int x=0;x<6;x++) {
1877 for (
int a=0;a<4;a++) {
1878 for (
int b=0;b<4;b++) {
1982 for (
int y=0;y<6;y++) {
2072 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);
2073 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);
2074 gslpp::complex B3HiggsnRte = (HpengTauEEENR0*SRE(0)/(mh*mh)+HpengTauEEENR1*SRE(1)/(mH*mH)+HpengTauEEENR2*SRE(2)/(mA*mA))/(4.0*pi*alph);
2075 gslpp::complex B3HiggsnLte = (HpengTauEEENL0*SLE(0)/(mh*mh)+HpengTauEEENL1*SLE(1)/(mH*mH)+HpengTauEEENL2*SLE(2)/(mA*mA))/(4.0*pi*alph);
2084 for (
int x=0;x<3;x++) {
2085 for (
int a=0;a<2;a++) {
2086 for (
int b=0;b<2;b++) {
2190 for (
int y=0;y<3;y++) {
2223 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))
2232 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))
2241 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))
2250 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))
2259 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))
2268 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))
2280 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);
2281 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);
2282 gslpp::complex B3HiggscRte = (HpengTauEEECR0*SRE(0)/(mh*mh)+HpengTauEEECR1*SRE(1)/(mH*mH)+HpengTauEEECR2*SRE(2)/(mA*mA))/(4.0*pi*alph);
2283 gslpp::complex B3HiggscLte = (HpengTauEEECL0*SLE(0)/(mh*mh)+HpengTauEEECL1*SLE(1)/(mH*mH)+HpengTauEEECL2*SLE(2)/(mA*mA))/(4.0*pi*alph);
2300 double MW =
mySUSY.Mw_tree();
2302 double sw2 =
mySUSY.StandardModel::sW2(MW);
2303 double stw =
sqrt(sw2);
2304 double ctw =
sqrt(1.0 - sw2);
2305 double ttw = stw/ctw;
2311 double cdenn = MW*
cosb;
2313 double g2t = g2/
sqrt(2.0);
2314 double alph =
mySUSY.getAle();
2318 for (
int a=0;a<4;a++) {
2319 for (
int x=0;x<6;x++) {
2321 NRlE.assign(a, x, - (g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 0) + (mE/cdenn)*
ON(a, 2)*
myRl(x, 3)));
2322 NRlMU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 1) + (mMU/cdenn)*
ON(a, 2)*
myRl(x, 4)));
2326 NLlE.assign(a, x, -(g2t)*((mE/cdenn)*
ON(a, 2)*
myRl(x, 0) + 2.0*
ON(a, 0)*ttw*
myRl(x, 3)));
2327 NLlMU.assign(a, x, -(g2t)*((mMU/cdenn)*
ON(a, 2)*
myRl(x, 1) + 2.0*
ON(a, 0)*ttw*
myRl(x, 4)));
2333 for (
int a=0;a<2;a++) {
2334 for (
int x=0;x<3;x++) {
2339 CLlE.assign(a, x, g2*mE/cdenc*
myU(a, 1).conjugate()*
myRn(x, 0));
2340 CLlMU.assign(a, x, g2*mMU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 1));
2343 for (
int a=0;a<2;a++) {
2344 for (
int x=0;x<6;x++) {
2359 for (
int a=0;a<4;a++) {
2360 for (
int b=0;b<4;b++) {
2361 for (
int x=0;x<6;x++) {
2362 for (
int y=0;y<6;y++) {
2397 for (
int a=0;a<2;a++) {
2398 for (
int b=0;b<2;b++) {
2399 for (
int x=0;x<3;x++) {
2400 for (
int y=0;y<6;y++) {
2436 double MZ =
mySUSY.getMz();
2437 double MW =
mySUSY.Mw_tree();
2439 double piconst = 1.0/(32.0 * pi * pi);
2440 double sw2 =
mySUSY.StandardModel::sW2(MW);
2441 double stw =
sqrt(sw2);
2442 double ctw =
sqrt(1.0 - sw2);
2443 double ttw = stw/ctw;
2449 double cdenn = MW*
cosb;
2451 double g2t = g2/
sqrt(2.0);
2458 for (
int a=0;a<4;a++) {
2459 for (
int x=0;x<6;x++) {
2461 NRlE.assign(a, x, - (g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 0) + (mE/cdenn)*
ON(a, 2)*
myRl(x, 3)));
2462 NRlMU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 1) + (mMU/cdenn)*
ON(a, 2)*
myRl(x, 4)));
2463 NRlTAU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 2) + (mTAU/cdenn)*
ON(a, 2)*
myRl(x, 5)));
2465 NLlE.assign(a, x, -(g2t)*((mE/cdenn)*
ON(a, 2)*
myRl(x, 0) + 2.0*
ON(a, 0)*ttw*
myRl(x, 3)));
2466 NLlMU.assign(a, x, -(g2t)*((mMU/cdenn)*
ON(a, 2)*
myRl(x, 1) + 2.0*
ON(a, 0)*ttw*
myRl(x, 4)));
2467 NLlTAU.assign(a, x, -(g2t)*((mTAU/cdenn)*
ON(a, 2)*
myRl(x, 2) + 2.0*
ON(a, 0)*ttw*
myRl(x, 5)));
2476 for (
int a=0;a<2;a++) {
2477 for (
int x=0;x<3;x++) {
2483 CLlE.assign(a, x, g2*mE/cdenc*
myU(a, 1).conjugate()*
myRn(x, 0));
2484 CLlMU.assign(a, x, g2*mMU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 1));
2485 CLlTAU.assign(a, x, g2*mTAU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 2));
2491 for (
int a=0;a<4;a++) {
2492 for (
int x=0;x<6;x++) {
2497 for (
int x=0;x<6;x++) {
2498 for (
int a=0;a<4;a++) {
2499 for (
int b=0;b<4;b++) {
2500 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) {
2506 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) {
2510 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) {
2514 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) {
2527 for (
int a=0;a<2;a++) {
2528 for (
int x=0;x<3;x++) {
2533 for (
int x=0;x<3;x++) {
2534 for (
int a=0;a<2;a++) {
2535 for (
int b=0;b<2;b++) {
2536 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) {
2542 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) {
2546 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) {
2550 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) {
2567 for (
int x=0;x<6;x++) {
2568 for (
int a=0;a<4;a++) {
2569 for (
int b=0;b<4;b++) {
2570 ZpengMuEEENR = ZpengMuEEENR -
NLlE(a,x)*
NLlMU(b,x)*piconst*
2572 ZpengMuEEENL = ZpengMuEEENL +
NRlE(a,x)*
NRlMU(b,x)*piconst*
2579 for (
int x=0;x<3;x++) {
2580 for (
int a=0;a<2;a++) {
2581 for (
int b=0;b<2;b++) {
2582 ZpengMuEEEC = ZpengMuEEEC +
CRlE(a,x)*
CRlMU(b,x)*piconst*
2589 FFunctions.assign(0, ZpengMuEEENR/(MZ*MZ*ctw*ctw) );
2590 FFunctions.assign(1, ZpengMuEEENR*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
2591 FFunctions.assign(2, (ZpengMuEEENL + ZpengMuEEEC)/(MZ*MZ*ctw*ctw) );
2592 FFunctions.assign(3, (ZpengMuEEENL + ZpengMuEEEC)*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
2599 for (
int x=0;x<6;x++) {
2600 for (
int a=0;a<4;a++) {
2601 for (
int b=0;b<4;b++) {
2602 ZpengTauMuMuMuNR = ZpengTauMuMuMuNR -
NLlMU(a,x)*
NLlTAU(b,x)*piconst*
2604 ZpengTauMuMuMuNL = ZpengTauMuMuMuNL +
NRlMU(a,x)*
NRlTAU(b,x)*piconst*
2611 for (
int x=0;x<3;x++) {
2612 for (
int a=0;a<2;a++) {
2613 for (
int b=0;b<2;b++) {
2614 ZpengTauMuMuMuC = ZpengTauMuMuMuC +
CRlMU(a,x)*
CRlTAU(b,x)*piconst*
2621 FFunctions.assign(0, ZpengTauMuMuMuNR/(MZ*MZ*ctw*ctw) );
2622 FFunctions.assign(1, ZpengTauMuMuMuNR*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
2623 FFunctions.assign(2, (ZpengTauMuMuMuNL + ZpengTauMuMuMuC)/(MZ*MZ*ctw*ctw) );
2624 FFunctions.assign(3, (ZpengTauMuMuMuNL + ZpengTauMuMuMuC)*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
2631 for (
int x=0;x<6;x++) {
2632 for (
int a=0;a<4;a++) {
2633 for (
int b=0;b<4;b++) {
2634 ZpengTauEEENR = ZpengTauEEENR -
NLlE(a,x)*
NLlTAU(b,x)*piconst*
2636 ZpengTauEEENL = ZpengTauEEENL +
NRlE(a,x)*
NRlTAU(b,x)*piconst*
2643 for (
int x=0;x<3;x++) {
2644 for (
int a=0;a<2;a++) {
2645 for (
int b=0;b<2;b++) {
2646 ZpengTauEEEC = ZpengTauEEEC +
CRlE(a,x)*
CRlTAU(b,x)*piconst*
2653 FFunctions.assign(0, ZpengTauEEENR/(MZ*MZ*ctw*ctw) );
2654 FFunctions.assign(1, ZpengTauEEENR*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
2655 FFunctions.assign(2, (ZpengTauEEENL + ZpengTauEEEC)/(MZ*MZ*ctw*ctw) );
2656 FFunctions.assign(3, (ZpengTauEEENL + ZpengTauEEEC)*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
2663 for (
int x=0;x<6;x++) {
2664 for (
int a=0;a<4;a++) {
2665 for (
int b=0;b<4;b++) {
2666 ZpengTauMuEENR = ZpengTauMuEENR -
NLlMU(a,x)*
NLlTAU(b,x)*piconst*
2668 ZpengTauMuEENL = ZpengTauMuEENL +
NRlMU(a,x)*
NRlTAU(b,x)*piconst*
2675 for (
int x=0;x<3;x++) {
2676 for (
int a=0;a<2;a++) {
2677 for (
int b=0;b<2;b++) {
2678 ZpengTauMuEEC = ZpengTauMuEEC +
CRlMU(a,x)*
CRlTAU(b,x)*piconst*
2685 FFunctions.assign(0, ZpengTauMuEENR/(MZ*MZ*ctw*ctw) );
2686 FFunctions.assign(1, ZpengTauMuEENR*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
2687 FFunctions.assign(2, (ZpengTauMuEENL + ZpengTauMuEEC)/(MZ*MZ*ctw*ctw) );
2688 FFunctions.assign(3, (ZpengTauMuEENL + ZpengTauMuEEC)*(-0.5+sw2)/(MZ*MZ*sw2*ctw*ctw) );
2698 double MW =
mySUSY.Mw_tree();
2700 double piconst = 1.0/(32.0 * pi * pi);
2701 double sw2 =
mySUSY.StandardModel::sW2(MW);
2702 double stw =
sqrt(sw2);
2703 double ctw =
sqrt(1.0 - sw2);
2704 double ttw = stw/ctw;
2708 double cdenn = MW*
cosb;
2710 double g2t = g2/
sqrt(2.0);
2715 for (
int a=0;a<4;a++) {
2716 for (
int x=0;x<6;x++) {
2718 NRlMU.assign(a, x, -(g2t)*((-
ON(a, 1) -
ON(a, 0)*ttw)*
myRl(x, 1) + (mMU/cdenn)*
ON(a, 2)*
myRl(x, 4)));
2720 NLlMU.assign(a, x, -(g2t)*((mMU/cdenn)*
ON(a, 2)*
myRl(x, 1) + 2.0*
ON(a, 0)*ttw*
myRl(x, 4)));
2727 for (
int a=0;a<2;a++) {
2728 for (
int x=0;x<3;x++) {
2732 CLlMU.assign(a, x, g2*mMU/cdenc*
myU(a, 1).conjugate()*
myRn(x, 1));
2736 for (
int a=0;a<4;a++) {
2737 for (
int x=0;x<6;x++) {
2742 for (
int a=0;a<2;a++) {
2743 for (
int x=0;x<3;x++) {
2748 for (
int a=0;a<4;a++) {
2749 for (
int x=0;x<6;x++) {
2750 if (fabs(1.0 -
Lepty(a, x)) > 0.01) {
2753 (6.0 *
pow((1.0 -
Lepty(a,x)),4.0)) );
2764 for (
int a=0;a<2;a++) {
2765 for (
int x=0;x<3;x++) {
2766 if(fabs(1.0-
Leptz(a, x)) > 0.01) {
2769 (6.0*
pow((1.0 -
Leptz(a, x)),4.0))) );
2784 for (
int a=0;a<4;a++) {
2785 for (
int x=0;x<6;x++) {
2786 g2ARN = g2ARN -mMU*mMU*piconst*(4.0*
NRlMU(a,x)*
NRlMU(a,x).conjugate()*
Leptf1(a,x)
2788 g2ALN = g2ALN -mMU*mMU*piconst*(4.0*
NLlMU(a,x)*
NLlMU(a,x).conjugate()*
Leptf1(a,x)
2796 for (
int a=0;a<2;a++) {
2797 for (
int x=0;x<3;x++) {
2798 g2ARC = g2ARC +mMU*mMU*piconst*(4.0*
CRlMU(a,x)*
CRlMU(a,x).conjugate()*
Leptf3(a,x)
2800 g2ALC = g2ALC +mMU*mMU*piconst*(4.0*
CLlMU(a,x)*
CLlMU(a,x).conjugate()*
Leptf3(a,x)
2820 __fPS_params ¶ms= *reinterpret_cast<__fPS_params *>(p);
2821 double r = params.
a*
log(x*(1.-x)/params.
a) / (x*(1.-x)-params.
a);
2831 gsl_integration_glfixed_table * w
2832 = gsl_integration_glfixed_table_alloc(100);
2836 F.params = reinterpret_cast<void *>(¶ms);
2838 result = gsl_integration_glfixed (&F, 0, 1, w);
2840 gsl_integration_glfixed_table_free (w);
2848 r=(2.*x-1.)*
fPS(x)-2.*x*(2.+
log(x));
2855 r=x*(2.+
log(x)-
fPS(x))/2.0;
2862 if (std::fabs(a-b) < 0.0005 && std::fabs(a-c) < 0.0005 && std::fabs(b-c) < 0.0005)
2864 r=9./(2.*(a+b+c)*(a+b+c));
2868 if (std::fabs(a-b) < 0.0005)
2870 if (std::fabs(a-c) < 0.0005)
2872 r=(-4.*(4.*b*b*c*c*
log(b*b/(c*c))
2873 +((b+c)*(b+c))*(c*c*
log((4.*c*c)/((b+c)*(b+c)))
2874 + b*b*
log(((b+c)*(b+c))/(4.*b*b)))))
2875 /(
pow(b-c,3)*(b+c)*(3.*b+c)*(b+3.*c));
2879 if (std::fabs(b-c) < 0.0005)
2881 r=(4.*(
pow(a-2.*b+c,2)*(-3.*
pow(a,8) - 46.*
pow(a,7)*c - 42.*
pow(a,6)*c*c + 234.*
pow(a,5)*c*c*c
2882 - 234.*a*a*a*
pow(c,5) + 42.*a*a*
pow(c,6) + 46.*a*
pow(c,7) + 3.*
pow(c,8)
2883 + 8.*a*a*c*c*(11.*
pow(a,4) + 28.*a*a*a*c + 50.*a*a*c*c + 28.*a*c*c*c + 11.*
pow(c,4))
2885 + a*a*(7.*
pow(a,6) + 62.*
pow(a,5)*c + 265.*
pow(a,4)*c*c
2886 + 356.*a*a*a*c*c*c + 265.*a*a*
pow(c,4) + 62.*a*
pow(c,5) + 7.*
pow(c,6))
2887 *
log((4.*a*a)/((a+c)*(a+c)))
2888 + 7.*
pow(a,6)*c*c*
log(((a+c)*(a+c))/(4.*c*c))
2889 + 62.*
pow(a,5)*c*c*c*
log(((a+c)*(a+c))/(4.*c*c))
2890 + 265.*
pow(a,4)*
pow(c,4)*
log(((a+c)*(a+c))/(4.*c*c))
2891 + 356.*a*a*a*
pow(c,5)*
log(((a+c)*(a+c))/(4.*c*c))
2892 + 265.*a*a*
pow(c,6)*
log(((a+c)*(a+c))/(4.*c*c))
2893 + 62.*a*
pow(c,7)*
log(((a+c)*(a+c))/(4.*c*c))
2894 + 7.*
pow(c,8)*
log(((a+c)*(a+c))/(4.*c*c)))
2895 +((a-c)*(a-c))*
pow(3.*a+c,2)*
pow(a+3.*c,2)
2896 *(4.*a*a*c*c*
log((c*c)/(a*a))
2897 +((a+c)*(a+c))*(a*a*
log((4.*a*a)/((a+c)*(a+c))) + c*c*
log(((a+c)*(a+c))/(4.*c*c))))))
2898 /(
pow(a-c,5)*(a+c)*
pow(3.*a+c,3)*
pow(a+3.*c,3));
2902 r=((b*b) - (c*c) - (c*c)*
log((b*b)/(c*c)))/
pow((b*b) - (c*c),2) +
2903 ((a - b)*(-
pow(b,4) +
pow(c,4) + 2*(b*b)*(c*c)*
log((b*b)/(c*c))))/(b*
pow((b*b) - (c*c),3));
2909 if (std::fabs(b-c) < 0.0005)
2911 if (std::fabs(a-c) < 0.0005)
2913 r=(-4*(4*(a*a)*(b*b)*
log((a*a)/(b*b)) + ((a+b)*(a+b))*((b*b)*
log((4*(b*b))/((a+b)*(a+b))) + (a*a)*
log(((a+b)*(a+b))/(4.*(a*a))))))/
2914 (
pow(a - b,3)*(a + b)*(3*a + b)*(a + 3*b));
2918 r=(-(a*a) + (c*c) - (a*a)*
log((c*c)/(a*a)))/
pow((a*a) - (c*c),2) +
2919 ((b - c)*(
pow(a,4) -
pow(c,4) + 2*(a*a)*(c*c)*
log((c*c)/(a*a))))/(c*
pow(-(a*a) + (c*c),3));
2924 if (std::fabs(a-c) < 0.0005)
2926 r=(-(b*b) + (c*c) + (b*b)*
log((b*b)/(c*c)))/
pow((b*b) - (c*c),2) +
2927 ((a - c)*(
pow(b,4) -
pow(c,4) - 2*(b*b)*(c*c)*
log((b*b)/(c*c))))/(c*
pow(-(b*b) + (c*c),3));
2931 r=( a*a*b*b*
log(a*a/(b*b))+b*b*c*c*
log(b*b/(c*c))+c*c*a*a*
log(c*c/(a*a)) )
2932 /( (a*a-b*b)*(b*b-c*c)*(a*a-c*c) );
2949 if (std::fabs(x-1.) < 0.005)
2951 r=1. + (1059.*(x-1.))/1175.;
2955 r=(1. - x) * (151. * x*x - 335. * x + 592.)
2956 + 6. * (21. * x*x*x - 108. * x*x - 93. * x + 50.) *
log(x)
2957 - 54. * x * (x*x - 2. * x - 2.) *
pow(
log(x),2)
2959 r*=4./(141.*
pow(1. - x,4));
2967 if (std::fabs(x-1.) < 0.005)
2969 r=1. - (45.*(x-1.))/122.;
2973 r=8.*(x*x-3.*x+2.)+(11.*x*x-40.*x+5.)*
log(x)
2974 -2.*(x*x-2.*x-2.)*
log(x)*
log(x)
2976 r=r*(-9.)/(122.*
pow(1-x,3));
2984 if (std::fabs(x-1.) < 0.005)
2986 r=1. + (76.*(x-1.))/875.;
2990 r=(1-x)*(-97.*x*x-529.*x+2.)+6.*x*x*(13.*x+81.)*
log(x)
2992 r=r*4./(105.*
pow(1-x,4));
3000 if (std::fabs(x-1.) < 0.005)
3002 r=1. - (111.*(x-1.)*(x-1.))/800.;
3007 r=r*(-9.)/(4.*
pow(1-x,3));
3015 if (std::fabs(x-y) < 0.005)
3017 if (std::fabs(y-1.) < 0.005)
3019 r=(49. - 22.*x - 22.*y + 10.*x*y)/60.;
3023 r=(2. + 3.*x - 6.*x*x + x*x*x + 6.*x*
log(x))/(4.*
pow(x-1.,4)*x)
3024 +(2. + 3.*y - 6.*y*y + y*y*y + 6.*y*
log(y))/(4.*
pow(y-1.,4)*y);
3029 if (std::fabs(x-1.) < 0.005)
3031 r=((x-1.)*(-25. + 48.*y - 36.*y*y + 16.*y*y*y - 3.*y*y*y*y - 12.*
log(y)))/(12.*
pow(y-1.,5))
3032 + (-11. + 18.*y - 9.*y*y + 2.*y*y*y - 6.*
log(y))/(6.*
pow(y-1.,4));
3034 else if (std::fabs(y-1.) < 0.005)
3036 r=((y-1.)*(-25. + 48.*x - 36.*x*x + 16.*x*x*x - 3.*x*x*x*x - 12.*
log(x)))/(12.*
pow(x-1.,5))
3037 + (-11. + 18.*x - 9.*x*x + 2.*x*x*x - 6.*
log(x))/(6.*
pow(x-1.,4));
3042 G3x=(1./(2.*
pow(x-1,3))) *( (x-1)*(x-3)+2.*
log(x) );
3043 G3y=(1./(2.*
pow(y-1,3))) *( (y-1)*(y-3)+2.*
log(y) );
3053 if (std::fabs(x-y) < 0.005)
3055 if (std::fabs(y-1.) < 0.005)
3057 r=(13. - 5.*x - 5.*y + 2.*x*y)/60.;
3061 r=((2. + 3.*x - 6.*x*x + x*x*x + 6.*x*
log(x))/
pow(x-1.,4)
3062 + (2. + 3.*y - 6.*y*y + y*y*y + 6.*y*
log(y))/
pow(y-1.,4))/12.;
3067 if (std::fabs(x-1.) < 0.005)
3069 r=(2. + 3.*y - 6.*y*y + y*y*y + 6.*y*
log(y))/(6.*
pow(y-1.,4))
3070 + ((x-1.)*(3. + 10.*y - 18.*y*y + 6.*y*y*y - y*y*y*y + 12.*y*
log(y)))/(12.*
pow(y-1.,5));
3072 else if (std::fabs(y-1.) < 0.005)
3074 r=(2. + 3.*x - 6.*x*x + x*x*x + 6.*x*
log(x))/(6.*
pow(x-1.,4))
3075 + ((y-1.)*(3. + 10.*x - 18.*x*x + 6.*x*x*x - x*x*x*x + 12.*x*
log(x)))/(12.*
pow(x-1.,5));
3080 G4x=(1./(2.*
pow(x-1,3))) *( (x-1)*(x+1) -2.*x*
log(x) );
3081 G4y=(1./(2.*
pow(y-1,3))) *( (y-1)*(y+1) -2.*y*
log(y) );
3093 double MZ =
mySUSY.getMz();
3096 double g1atMZ = 0.357456;
3097 double g2atMZ = 0.651721;
3098 double g3atMZ =
sqrt(4.0*pi*
mySUSY.getAlsMz());
3099 double mh =
mySUSY.getMHl();
3100 double mhh =
mySUSY.getMHh();
3101 double ma =
mySUSY.getMHa();
3106 double vew =
v/
sqrt(2.);
3128 MsQhat2 =
mySUSY.getMsQhat2();
3129 MsUhat2 =
mySUSY.getMsUhat2();
3130 MsDhat2 =
mySUSY.getMsDhat2();
3131 MsLhat2 =
mySUSY.getMsLhat2();
3132 MsEhat2 =
mySUSY.getMsEhat2();
3133 double msq1L = MsQhat2(0,0).real();
3134 double msq2L = MsQhat2(1,1).real();
3135 double msq3L = MsQhat2(2,2).real();
3136 double msuR = MsUhat2(0,0).real();
3137 double mscR = MsUhat2(1,1).real();
3138 double mstR = MsUhat2(2,2).real();
3139 double msdR = MsDhat2(0,0).real();
3140 double mssR = MsDhat2(1,1).real();
3141 double msbR = MsDhat2(2,2).real();
3142 double mseL = MsLhat2(0,0).real();
3143 double msmuL = MsLhat2(1,1).real();
3144 double mstauL = MsLhat2(2,2).real();
3145 double mseR = MsEhat2(0,0).real();
3146 double msmuR = MsEhat2(1,1).real();
3147 double mstauR = MsEhat2(2,2).real();
3154 b1[0]=41./6., b1[1]=-19./6., b1[2]=-7.;
3155 gi[0]=g1atMZ, gi[1]=g2atMZ, gi[2]=g3atMZ;
3157 res_g=1./gi[k]/gi[k] -b1[k]/(8.*pi*pi)*
log(
sqrt(msmuL)/MZ);
3158 gi[k]=
sqrt(1/res_g);
3160 g1=gi[0], g2=gi[1], g3=gi[2];
3162 double alp = (g1*g1*g2*g2/(g1*g1+g2*g2))/(4.0*pi);
3165 double mzq =
sqrt( 0.5*(g1*g1+g2*g2)*vew*vew );
3166 double mwq =
sqrt( 0.5*(g2*g2)*vew*vew );
3167 double sw2 = g1*g1/(g1*g1+g2*g2);
3171 double msneu2 = msmuL+0.5*mzq*mzq*c2b;
3173 TEhat.assign(1,1, 0.);
3181 Rsmu.assign(0,0, msmuL+mmu*mmu+c2b*mzq*mzq*(-1.0/2.0+sw2) );
3184 Rsmu.assign(1,1, msmuR+mmu*mmu-c2b*mzq*mzq*sw2);
3185 Rsmu.eigensystem(Xm,msmu2);
3193 double gm21loop = 1.06978e-9;
3200 double gm2cor, res2, res3;
3201 double res01, res02;
3215 for(
int i=0; i<2; i++){
3217 lxh1.assign(i, tmp1*(
myU(i,0)*
myV(i,1)*ca +
myU(i,1)*
myV(i,0)*(-sa) ) );
3218 lxh2.assign(i, tmp1*(
myU(i,0)*
myV(i,1)*sa +
myU(i,1)*
myV(i,0)*( ca) ) );
3223 double xps, xh1, xh2;
3224 for(
int i=0; i<2; i++){
3229 res += (lA*lxA(i)).real()*
fPS(xps) + (lh1*lxh1(i)).real()*
fS(xh1)
3230 +(lh2*lxh2(i)).real()*
fS(xh2);
3232 res *= alp*alp * mmu*mmu / (8.*pi*pi*mwq*mwq*sw2);
3240 stauM.assign(0,0, mstauL+mTAU*mTAU + mzq*mzq*c2b*(0.5-(-1.)*sw2) );
3241 stauM.assign(1,1, mstauR+mTAU*mTAU + mzq*mzq*c2b*(-1.)*sw2);
3242 stauM.assign(0,1, mTAU*(a3tau-muH*
tanb));
3243 stauM.assign(1,0, mTAU*(a3tau-muH*
tanb));
3245 sbottomM.assign(0,0, msq3L+mb*mb + mzq*mzq*c2b*(0.5-(-1./3.)*sw2) );
3246 sbottomM.assign(1,1, msbR+mb*mb + mzq*mzq*c2b*(-1./3.)*sw2);
3247 sbottomM.assign(0,1, mb*(a3b-muH*
tanb));
3248 sbottomM.assign(1,0, mb*(a3b-muH*
tanb));
3250 stopM.assign(0,0, msq3L+mt*mt + mzq*mzq*c2b*(0.5-(2./3.)*sw2) );
3251 stopM.assign(1,1, mstR+mt*mt + mzq*mzq*c2b*(2./3.)*sw2);
3252 stopM.assign(0,1, mt*(a3t-muH/
tanb));
3253 stopM.assign(1,0, mt*(a3t-muH/
tanb));
3255 stauM.eigensystem(Ustau,mstau2);
3256 sbottomM.eigensystem(Usbottom,msbottom2);
3257 stopM.eigensystem(Ustop,mstop2);
3259 Ustau=Ustau.hconjugate();
3260 Usbottom=Usbottom.hconjugate();
3261 Ustop=Ustop.hconjugate();
3271 for(
int i=0; i<2; i++){
3273 rr=2.*mTAU/(
cosb*mstau2(i)) * Ustau(i,0).
conjugate() * Ustau(i,1);
3274 lstauh1.assign(i,rr*(-muH*ca + a3tau*(-sa)) );
3275 lstauh2.assign(i,rr*(-muH*sa + a3tau*(ca)) );
3277 rr=2.*mb/(
cosb*msbottom2(i)) * Usbottom(i,0).
conjugate() * Usbottom(i,1);
3278 lsbottomh1.assign(i,rr*(-muH*ca + a3b*(-sa)) );
3279 lsbottomh2.assign(i,rr*(-muH*sa + a3b*(ca)) );
3281 rr=2.*mt/(
cosb*mstop2(i)) * Ustop(i,0).conjugate() * Ustop(i,1);
3282 lstoph1.assign(i,rr*(-muH*sa + a3t*(ca)) );
3283 lstoph2.assign(i,rr*(-muH*(-ca) + a3t*(sa)) );
3292 for(
int i=0;i<2;i++){
3293 xx1=mstau2(i)/(mh*mh);
3294 xx2=mstau2(i)/(mhh*mhh);
3295 res1 += (lh1*lstauh1(i)).real()*
fft(xx1)
3296 +(lh2*lstauh2(i)).real()*
fft(xx2);
3300 xx1=msbottom2(i)/(mh*mh);
3301 xx2=msbottom2(i)/(mhh*mhh);
3302 res1 += (lh1*lsbottomh1(i)).real()*qe2*
fft(xx1)
3303 +(lh2*lsbottomh2(i)).real()*qe2*
fft(xx2);
3307 xx1=mstop2(i)/(mh*mh);
3308 xx2=mstop2(i)/(mhh*mhh);
3309 res1 += (lh1*lstoph1(i)).real()*qe2*
fft(xx1)
3310 +(lh2*lstoph2(i)).real()*qe2*
fft(xx2);
3313 res1 *= alp*alp * mmu*mmu / (8.*pi*pi*mwq*mwq*sw2);
3325 gm2cor=0, res2=0, res3=0;
3336 double x0, x1a, x2a, xL, xR;
3339 x0=
sqrt(std::fabs(msneu2));
3340 xL=
sqrt( msmuL-mzq*mzq*(sw2-0.5) );
3341 xR=
sqrt( msmuR +mzq*mzq*sw2 );
3343 tmp2=M2.
abs2()+muH.
abs2() +2.*mwq*mwq;
3344 tmp3= tmp2*tmp2 -4.*M2.
abs2()*muH.
abs2();
3345 x1a=
sqrt( 0.5*( tmp2-
sqrt(tmp3) ) );
3346 x2a=
sqrt( 0.5*( tmp2+
sqrt(tmp3) ) );
3348 dmu=-muH.
real()*
tanb*g2*g2*M2.
real()/(16.*pi*pi)*(
It(x1a,x2a,x0) + 0.5*
It(x1a,x2a,xL))
3352 gm2cor=gm21loop/(1+dmu);
3361 res2=gm2cor*alp/(4.*pi)*16.*
log(mmu/
sqrt(msmuL));
3370 if (sub_leading==1){
3375 double ymu=mmu/(vew*
cosb);
3377 ckR.assign(0, ymu*
myU(0,1));
3378 ckR.assign(1, ymu*
myU(1,1));
3381 ckL.assign(0, -g2*
myV(0,0));
3382 ckL.assign(1, -g2*
myV(1,0));
3385 for(
int i=0; i<2; i++){
3387 amch=amch -(47.*mmu/(12.*msneu2))*( (ckL(i).abs()*ckL(i).abs() + ckR(i).abs()*ckR(i).abs())*
F3C(xk) );
3388 amch=amch -(122.*
MChi(i)/(9.*msneu2))*
F4C(xk)*( (ckL(i)*ckR(i)).real() );
3390 amch=(1./(1+dmu))*amch*mmu/(16.*pi*pi)*alp/(4.*pi);
3397 for(
int i=0; i<4; i++){
3398 nR1.assign(i,
sqrt(2.)*g1*
myN(i,0)*Xm(0,1) + ymu*
myN(i,2)*Xm(0,0));
3399 nR2.assign(i,
sqrt(2.)*g1*
myN(i,0)*Xm(1,1) + ymu*
myN(i,2)*Xm(1,0));
3400 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());
3401 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());
3407 for(
int i=0; i<4; i++){
3410 tmp4=nL1(i).abs2() + nR1(i).abs2();
3411 tmp5=nL2(i).abs2() + nR2(i).abs2();
3412 r1=35.*mmu/(72.*msmu2(0))*
F3N(xi1)*tmp4 - 16.*
MChi0(i)/(9.*msmu2(0))*
F4N(xi1)*( (nL1(i)*nR1(i)).real() );
3413 r2=35.*mmu/(72.*msmu2(1))*
F3N(xi2)*tmp5 - 16.*
MChi0(i)/(9.*msmu2(1))*
F4N(xi2)*( (nL2(i)*nR2(i)).real() );
3417 amne*=(1./(1.+dmu))*mmu/(16.*pi*pi)*alp/(4.*pi);
3421 res2=res2+amch+amne;
3434 double x1,y1,x2,y2,x3,y3,x4,y4,x5,y5;
3435 double awhn=mmu*mmu*muH.
real()*
tanb/(1+dmu);
3442 x1=M2.
abs2()/msneu2;
3443 y1=muH.
abs2()/msneu2;
3444 awhn *= g2*g2*M2.
real()/(8.*pi*pi*msneu2*msneu2)*
Fa(x1,y1);
3448 y2=muH.
abs2()/msmuL;
3449 awhl *= -g2*g2*M2.
real()/(16.*pi*pi*
pow(msmuL,2))*
Fb(x2,y2);
3453 y3=muH.
abs2()/msmuL;
3454 abhl *= g1*g1*M1.
real()/(16.*pi*pi*
pow(msmuL,2))*
Fb(x3,y3);
3458 y4=muH.
abs2()/msmuR;
3459 abhr *= -g1*g1*M1.
real()/(8.*pi*pi*
pow(msmuR,2))*
Fb(x4,y4);
3462 x5=msmuL/(M1.
abs2());
3463 y5=msmuR/(M1.
abs2());
3464 ablr *= g1*g1/(8.*pi*pi*
pow(M1.
real(),3))*
Fb(x5,y5);
3466 double dg2, dg1, dh, dwh, dbh, dtb;
3467 double msusy=
sqrt(msmuL);
3469 dg1=g1*g1/(16.*pi*pi)*(4./3.)
3476 dg2=g2*g2/(16.*pi*pi)*(4./3.)
3481 double yb, ytau, yt;
3482 double as_mt, delta_mt;
3487 as_mt=as_mt - (-7.)*
log(mt/
sqrt(msmuL))/(8.*pi*pi);
3488 as_mt=(1./as_mt)/(4.*pi);
3490 delta_mt=-4./3.*(as_mt/pi)-9.1*
pow( (as_mt/pi), 2)-80.*
pow((as_mt/pi),3);
3493 yt=mt/(vew*
sinb)*(1+delta_mt);
3494 ytau=mTAU/(vew*
cosb);
3497 yt=yt*(1+as_mt/(8.*pi)*(-4./3.));
3498 yb=yb*(1+as_mt/(8.*pi)*(-4./3.));
3500 dh=0.5/(16.*pi*pi)*(3.*yt*yt*
log(
sqrt(mstR)/msusy) +3.*yb*yb*
log(
sqrt(msbR)/msusy)
3501 +3.*(yt*yt+yb*yb)*
log(
sqrt(msq3L)/msusy)
3502 +ytau*ytau*
log(
sqrt(mstauR)/msusy) +ytau*ytau*
log(
sqrt(mstauL)/msusy));
3504 dwh=yt*yt/(16.*pi*pi)*(-6.*
log(
sqrt(msq3L)/msusy));
3506 dbh=yt*yt/(16.*pi*pi)*( 2.*
log(
sqrt(msq3L)/msusy)-8.*
log(
sqrt(mstR)/msusy) );
3510 dtb=1./(16.*pi*pi)*( 3.*yb*yb -3.*yt*yt +ytau*ytau)*
log(
Q_S/msusy);
3513 res3= awhn*(dg2 + dh + dwh + dtb)
3514 + awhl*(dg2 + dh + dwh + dtb)
3515 + abhl*(dg1 + dh + dbh + dtb)
3516 + abhr*(dg1 + dh + dbh + dtb)
3526 res01=gm2cor+res2+res3-gm21loop;
3549 double MZ =
mySUSY.getMz();
3550 double MW =
mySUSY.Mw_tree();
3551 double sw2 =
mySUSY.StandardModel::sW2(MW);
3552 double cw2 = 1.0 - sw2;
3558 + (1.0-1.0/(4.0*sw2))/(cw2*MZ*MZ)*
FFunctions(li_to_lj)(1)
3561 + (1.0-1.0/(4.0*sw2))/(cw2*MZ*MZ)*
FFunctions(li_to_lj)(0)
3569 double MZ =
mySUSY.getMz();
3570 double MW =
mySUSY.Mw_tree();
3571 double sw2 =
mySUSY.StandardModel::sW2(MW);
3572 double cw2 = 1.0 - sw2;
3577 C10.assign(0, 1.0/(4.0*sw2*cw2*MZ*MZ)*
FFunctions(li_to_lj)(1)
3579 C10.assign(1, 1.0/(4.0*sw2*cw2*MZ*MZ)*
FFunctions(li_to_lj)(0)
3620 vmcDLij = StandardModelMatching::CMDLij(li_lj);
3631 std::stringstream out;
3633 throw std::runtime_error(
"SUSYMatching::CMDLij(): order " + out.str() +
" not implemented.\nFor lepton flavour violating observables only Leading Order (LO) necessary.");
3636 vmcDLij.push_back(
mcDLij);
3643 vmcDLi3j = StandardModelMatching::CMDLi3j(li_lj);
3687 std::stringstream out;
3689 throw std::runtime_error(
"SUSYMatching::CMDLi3j(): order " + out.str() +
" not implemented.\nFor lepton flavour violating observables only Leading Order (LO) necessary.");
3699 double MW =
mySUSY.Mw_tree();
3700 double sw2 =
mySUSY.StandardModel::sW2(MW);
3702 vmcmueconv = StandardModelMatching::CMmueconv();
3721 std::stringstream out;
3723 throw std::runtime_error(
"SUSYMatching::CMmueconv(): order " + out.str() +
" not implemented.\nFor lepton flavour violating observables only Leading Order (LO) necessary.");
3733 vmcgminus2mu = StandardModelMatching::CMgminus2mu();
3749 std::stringstream out;
3751 throw std::runtime_error(
"SUSYMatching::CMgminus2mu(): order " + out.str() +
" not implemented.\nFor lepton flavour violating observables only Leading Order (LO) necessary.");
3755 return(vmcgminus2mu);