Calculates Higgs penguin amplitudes for the process \( \ell_j \to \ell_i \ell_i \ell_i \) from [19].
Calculates Higgs penguin amplitudes for m->3e (1), t->3m (2) and t->3e (3)**/ /**
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);