--- /dev/null
+#include "isajet/pilot.h"
+ SUBROUTINE SSGLBF
+C-----------------------------------------------------------------------
+C
+C This subroutine gives gluino branching fractions to gauginos
+C according to Baer,Barger,Karatas,Tata (Phys.Rev.D36,96(1987)
+C (Now includes gluino->gluon+zino1,2,3,4 loop decays. Jan 1990)
+C Also includes contribution due to non-degenerate t1-t2 stops
+C Also includes contribution due to non-degenerate b_L and b_R
+C Updated to include mixed sbottom states b1 and b2: 10/9/96
+C
+C Auxiliary functions are called SSGxyi, where normally x
+C indicates the SUSY particle, y the SM particle(s), and i is
+C a counter.
+C
+C Baer's GLUBF
+C
+C-----------------------------------------------------------------------
+#if defined(CERNLIB_IMPNONE)
+ IMPLICIT NONE
+#endif
+#include "isajet/sslun.inc"
+#include "isajet/ssmode.inc"
+#include "isajet/sssm.inc"
+#include "isajet/sspar.inc"
+#include "isajet/sstmp.inc"
+#include "isajet/sstype.inc"
+C
+ EXTERNAL SSGWQ1,SSGWQ2,SSGZG1,SSGZG2,SSGZG3,SSGZT
+ EXTERNAL SSGWT1,SSGWT2,SSGWT3,SSGWT4,SSGWT5,SSGWT6,
+ $SSGWT7,SSGWT8
+ EXTERNAL SSGX1,SSGX2,SSGX3,SSGX4,SSGX5,SSGX6,SSGX7,SSGX8
+ $,SSGX9
+ REAL WIDU,WIDD,WIDI,WIDL,WIDR
+ REAL SSXINT,SSXLAM,XUPPER,XTCON,SUM,GMQK,C,G,FT,GP
+ $,UPPER,GS,FB,TANB,PI,SR2,UPPR,ALF3
+ REAL MW1,MW2,SNW1,SNW2,XM,YM,THX,THY,COST,SINT,COSB,SINB
+ REAL MZIZ,AUIZ,ADIZ,BUIZ,BDIZ
+ REAL WID,SNIZ,FACT
+ REAL XT1,XT2,XT3,XT4,XT5,XT6,XT7,XT8,XT9,AL1,AL2,BE1,BE2
+ REAL XLL,XRR,XL1R1,XL2R2,XL1R2,XL2R1,XL1L2,XR1R2,XLR1
+ REAL XMST1,XMST2,XMST12,XLR2,XMSB1,XMSB2,XMSB12
+ REAL BW1,BW2,GT1,GT2,GT1T2,GB1,GB2,GT1B1,GT1B2,GT2B1,GT2B2
+ REAL KUL,KUR,KDL,KDR,KCL,KCR,KSL,KSR,KBL,KBR,KT1,KT2
+ REAL XKUL,XKUR,XKDL,XKDR,XKSL,XKSR,XKCL,XKCR,XKBL,XKBR,
+ $XKT1,XKT2
+ REAL XI1UL,XI1UR,XI1DL,XI1DR,XI1SL,XI1SR,XI1CL,XI1CR,
+ $XI1BL,XI1BR,XI1T1,XI1T2,XIT1,XIT2
+ REAL ALT1,ART1,ALT2,ART2,ALB1,ARB1,ALB2,ARB2
+ REAL BETA,ASMB,MBMB,MBQ,ASMT,MTMT,MTQ,SUALFS,AMPL
+ INTEGER IZ,ISZI(4),THIZ
+ COMPLEX ZONE,ZI,ZAT1(2),ZAT2(2),ZADW1,ZADW2,ZAUW1,ZAUW2
+ COMPLEX ZAUIZ,ZADIZ,ZBUIZ,ZBDIZ,Z1(2),Z2(2)
+ DOUBLE PRECISION SSALFS,SSMQCD
+ SAVE ZONE,ZI
+ DATA ZONE,ZI/(1.,0.),(0.,1.)/
+
+C
+C Partly duplicated from SSMASS.
+C
+ AMPL=2.4E18
+ PI=4.*ATAN(1.)
+ SR2=SQRT(2.)
+ G=SQRT(4*PI*ALFAEM/SN2THW)
+ GP=G*SQRT(SN2THW/(1.-SN2THW))
+ ALF3=SSALFS(DBLE(AMGLSS**2))
+ GS=SQRT(4.*PI*ALF3)
+ TANB=1./RV2V1
+ BETA=ATAN(TANB)
+ ASMB=SUALFS(AMBT**2,.36,AMTP,3)
+ MBMB=AMBT*(1.-4*ASMB/3./PI)
+ MBQ=SSMQCD(DBLE(MBMB),DBLE(AMGLSS))
+ ASMT=SUALFS(AMTP**2,.36,AMTP,3)
+ MTMT=AMTP/(1.+4*ASMT/3./PI+(16.11-1.04*(5.-6.63/AMTP))*
+ $(ASMT/PI)**2)
+ MTQ=SSMQCD(DBLE(MTMT),DBLE(AMGLSS))
+ FB=G*MBQ/SR2/AMW/COS(BETA)
+ FT=G*MTQ/SR2/AMW/SIN(BETA)
+ MW1=ABS(AMW1SS)
+ MW2=ABS(AMW2SS)
+ SNW1=SIGN(1.,AMW1SS)
+ SNW2=SIGN(1.,AMW2SS)
+ XM=1./TAN(GAMMAL)
+ YM=1./TAN(GAMMAR)
+ THX=SIGN(1.,XM)
+ THY=SIGN(1.,YM)
+ COST=COS(THETAT)
+ SINT=SIN(THETAT)
+ COSB=COS(THETAB)
+ SINB=SIN(THETAB)
+C
+ ZADW1=ZI*G*SNW1*SIN(GAMMAR)
+ ZAUW1=ZI*G*SIN(GAMMAL)
+ ZADW2=ZI*G*SNW2*COS(GAMMAR)*THY
+ ZAUW2=ZI*G*COS(GAMMAL)*THX
+ BW1=-FT*SNW1*COS(GAMMAR)
+ BW2=FT*SNW2*SIN(GAMMAR)*THY
+ ZAT1(1)=ZADW1*COST+ZI*BW1*SINT
+ ZAT2(1)=ZADW1*SINT-ZI*BW1*COST
+ ZAT1(2)=ZADW2*COST+ZI*BW2*SINT
+ ZAT2(2)=ZADW2*SINT-ZI*BW2*COST
+ C=SSALFS(DBLE(AMGLSS**2))*AMGLSS/8./PI**2
+C
+C gluino --> w1 + qk + qb
+C
+C Now includes sbottom as well as stop mixing/Yukawa effects
+C as of 3/31/97, thanks to M. Drees
+ TMP(1)=MW1
+ UPPR=(AMGLSS**2-MW1**2)/2./AMGLSS
+ IF (AMGLSS.GT.(MW1+AMUP+AMDN)) THEN
+ IF (AMGLSS.LT.AMULSS.AND.AMGLSS.LT.AMDLSS) THEN
+ TMP(2)=AMULSS
+ WIDU=ZADW1*CONJG(ZADW1)*SSXINT(0.,SSGWQ1,UPPR)
+ TMP(2)=AMDLSS
+ WIDD=ZAUW1*CONJG(ZAUW1)*SSXINT(0.,SSGWQ1,UPPR)
+ WIDI=2*REAL(ZAUW1*ZADW1)*SSXINT(0.,SSGWQ2,UPPR)
+ WID=GS**2/2./AMGLSS/(2*PI)**5*(WIDU+WIDD+WIDI)
+ ELSE IF (AMGLSS.LT.AMULSS.AND.AMGLSS.GE.AMDLSS) THEN
+ TMP(2)=AMULSS
+ WIDU=ZADW1*CONJG(ZADW1)*SSXINT(0.,SSGWQ1,UPPR)
+ WID=GS**2/2./AMGLSS/(2*PI)**5*WIDU
+ ELSE IF (AMGLSS.LT.AMDLSS.AND.AMGLSS.GE.AMULSS) THEN
+ TMP(2)=AMDLSS
+ WIDD=ZAUW1*CONJG(ZAUW1)*SSXINT(0.,SSGWQ1,UPPR)
+ WID=GS**2/2./AMGLSS/(2*PI)**5*WIDD
+ ELSE
+ WID=0.
+ END IF
+ CALL SSSAVE(ISGL,WID,+ISW1,+IDDN,-IDUP,0,0)
+ CALL SSSAVE(ISGL,WID,-ISW1,+IDUP,-IDDN,0,0)
+ Z1(1)=1.
+ Z1(2)=-Z1(1)
+ Z2(1)=G*SIN(GAMMAL)
+ Z2(2)=Z2(1)
+ CALL SSME3(2,AMDLSS,Z1,Z2)
+ Z1(1)=G*SIN(GAMMAR)
+ Z1(2)=-Z1(1)
+ Z2(1)=1.
+ Z2(2)=Z2(1)
+ CALL SSME3(3,AMULSS,Z1,Z2)
+ END IF
+C
+ IF (AMGLSS.GT.(MW1+AMST+AMCH)) THEN
+ IF (AMGLSS.LT.AMCLSS.AND.AMGLSS.LT.AMSLSS) THEN
+ TMP(2)=AMCLSS
+ WIDU=ZADW1*CONJG(ZADW1)*SSXINT(0.,SSGWQ1,UPPR)
+ TMP(2)=AMSLSS
+ WIDD=ZAUW1*CONJG(ZAUW1)*SSXINT(0.,SSGWQ1,UPPR)
+ WIDI=2*REAL(ZAUW1*ZADW1)*SSXINT(0.,SSGWQ2,UPPR)
+ WID=GS**2/2./AMGLSS/(2*PI)**5*(WIDU+WIDD+WIDI)
+ ELSE IF (AMGLSS.LT.AMCLSS.AND.AMGLSS.GE.AMSLSS) THEN
+ TMP(2)=AMCLSS
+ WIDU=ZADW1*CONJG(ZADW1)*SSXINT(0.,SSGWQ1,UPPR)
+ WID=GS**2/2./AMGLSS/(2*PI)**5*WIDU
+ ELSE IF (AMGLSS.LT.AMSLSS.AND.AMGLSS.GE.AMCLSS) THEN
+ TMP(2)=AMSLSS
+ WIDD=ZAUW1*CONJG(ZAUW1)*SSXINT(0.,SSGWQ1,UPPR)
+ WID=GS**2/2./AMGLSS/(2*PI)**5*WIDD
+ ELSE
+ WID=0.
+ END IF
+ CALL SSSAVE(ISGL,WID,+ISW1,+IDST,-IDCH,0,0)
+ CALL SSSAVE(ISGL,WID,-ISW1,+IDCH,-IDST,0,0)
+ Z1(1)=1.
+ Z1(2)=-Z1(1)
+ Z2(1)=G*SIN(GAMMAL)
+ Z2(2)=Z2(1)
+ CALL SSME3(2,AMSLSS,Z1,Z2)
+ Z1(1)=G*SIN(GAMMAR)
+ Z1(2)=-Z1(1)
+ Z2(1)=1.
+ Z2(2)=Z2(1)
+ CALL SSME3(3,AMCLSS,Z1,Z2)
+ END IF
+C
+ IF (AMGLSS.GT.(MW1+AMTP+AMBT)) THEN
+ ALT1=-G*SIN(GAMMAR)*COST+FT*COS(GAMMAR)*SINT
+ ART1=-FB*COS(GAMMAL)*COST
+ ALT2=-G*SIN(GAMMAR)*SINT-FT*COS(GAMMAR)*COST
+ ART2=-FB*COS(GAMMAL)*SINT
+ ALB1=-G*SIN(GAMMAL)*COSB+FB*COS(GAMMAL)*SINB
+ ARB1=-FT*COS(GAMMAR)*COSB
+ ALB2=-G*SIN(GAMMAL)*SINB-FB*COS(GAMMAL)*COSB
+ ARB2=-FT*COS(GAMMAR)*SINB
+ UPPER=(AMGLSS**2+AMTP**2-(MW1+AMBT)**2)/2./AMGLSS
+ FACT=GS**2*PI**2/(2*PI)**5/2./AMGLSS
+ TMP(1)=MW1
+ TMP(2)=AMGLSS
+ TMP(3)=AMTP
+ IF (AMGLSS.LT.(AMTP+AMT1SS)) THEN
+ TMP(6)=AMT1SS
+ TMP(7)=AMT1SS
+ GT1=(ALT1**2+ART1**2)*(SSXINT(AMTP,SSGWT1,UPPER)+
+ , SIN(2*THETAT)*SSXINT(AMTP,SSGWT8,UPPER))
+ TMP(7)=AMT2SS
+ GT1T2=2*(ALT1*ALT2+ART1*ART2)*(SINT**2-COST**2)*
+ , SSXINT(AMTP,SSGWT8,UPPER)
+ ELSE
+ GT1=0.
+ GT1T2=0.
+ END IF
+ IF (AMGLSS.LT.(AMTP+AMT2SS)) THEN
+ TMP(6)=AMT2SS
+ TMP(7)=AMT2SS
+ GT2=(ALT2**2+ART2**2)*(SSXINT(AMTP,SSGWT1,UPPER)-
+ , SIN(2*THETAT)*SSXINT(AMTP,SSGWT8,UPPER))
+ ELSE
+ GT2=0.
+ END IF
+ IF (AMGLSS.LT.(AMB1SS+AMBT)) THEN
+ TMP(4)=AMB1SS
+ TMP(8)=SNW1
+C Rewrite UPPR=(AMGLSS**2+AMBT**2-(AMTP+MW1)**2)/2./AMGLSS
+ UPPR=((AMGLSS-AMTP-MW1)*(AMGLSS+AMTP+MW1)+AMBT**2)/(2*AMGLSS)
+ GB1=(ALB1**2+ARB1**2)*SSXINT(AMBT,SSGWT2,UPPR)-
+ , ALB1*ARB1*SSXINT(AMBT,SSGWT3,UPPR)
+ ELSE
+ GB1=0.
+ END IF
+ IF (AMGLSS.LT.(AMB2SS+AMBT)) THEN
+ TMP(4)=AMB2SS
+ TMP(8)=SNW1
+C Rewrite UPPR=(AMGLSS**2+AMBT**2-(AMTP+MW1)**2)/2./AMGLSS
+ UPPR=((AMGLSS-AMTP-MW1)*(AMGLSS+AMTP+MW1)+AMBT**2)/(2*AMGLSS)
+ GB2=(ALB2**2+ARB2**2)*SSXINT(AMBT,SSGWT2,UPPR)-
+ , ALB2*ARB2*SSXINT(AMBT,SSGWT3,UPPR)
+ ELSE
+ GB2=0.
+ END IF
+ IF (AMGLSS.LT.(AMTP+AMT1SS).AND.AMGLSS.LT.
+ $ (AMB1SS+AMBT)) THEN
+ TMP(4)=AMB1SS
+ TMP(6)=AMT1SS
+ TMP(8)=SNW1
+ GT1B1=(COST*SINB*ART1*ALB1+SINT*COSB*ALT1*ARB1)*
+ , SSXINT(AMTP,SSGWT6,UPPER)-
+ , (COST*COSB*ALT1*ALB1+SINT*SINB*ART1*ARB1)*
+ , SSXINT(AMTP,SSGWT4,UPPER)+
+ , (COST*COSB*ALT1*ARB1+SINT*SINB*ART1*ALB1)*
+ , SSXINT(AMTP,SSGWT5,UPPER)-
+ , (COST*SINB*ART1*ARB1+SINT*COSB*ALT1*ALB1)*
+ , SSXINT(AMTP,SSGWT7,UPPER)
+ ELSE
+ GT1B1=0.
+ END IF
+ IF (AMGLSS.LT.(AMTP+AMT1SS).AND.AMGLSS.LT.
+ $ (AMB2SS+AMBT)) THEN
+ TMP(4)=AMB2SS
+ TMP(6)=AMT1SS
+ TMP(8)=SNW1
+ GT1B2=(-COST*COSB*ART1*ALB2+SINT*SINB*ALT1*ARB2)*
+ , SSXINT(AMTP,SSGWT6,UPPER)-
+ , (COST*COSB*ALT1*ALB2-SINT*COSB*ART1*ARB2)*
+ , SSXINT(AMTP,SSGWT4,UPPER)+
+ , (COST*SINB*ALT1*ARB2-SINT*COSB*ART1*ALB2)*
+ , SSXINT(AMTP,SSGWT5,UPPER)-
+ , (-COST*COSB*ART1*ARB2+SINT*SINB*ALT1*ALB2)*
+ , SSXINT(AMTP,SSGWT7,UPPER)
+ ELSE
+ GT1B2=0.
+ END IF
+ IF (AMGLSS.LT.(AMTP+AMT2SS).AND.AMGLSS.LT.
+ $ (AMB1SS+AMBT)) THEN
+ TMP(4)=AMB1SS
+ TMP(6)=AMT2SS
+ TMP(8)=SNW1
+ GT2B1=(SINT*SINB*ART2*ALB1-COST*COSB*ALT2*ARB1)*
+ , SSXINT(AMTP,SSGWT6,UPPER)-
+ , (SINT*COSB*ALT2*ALB1-COST*SINB*ART2*ARB1)*
+ , SSXINT(AMTP,SSGWT4,UPPER)+
+ , (SINT*COSB*ALT2*ARB1-COST*SINB*ART2*ALB1)*
+ , SSXINT(AMTP,SSGWT5,UPPER)-
+ , (SINT*SINB*ART2*ARB1-COST*COSB*ALT2*ALB1)*
+ , SSXINT(AMTP,SSGWT7,UPPER)
+ ELSE
+ GT2B1=0.
+ END IF
+ IF (AMGLSS.LT.(AMTP+AMT2SS).AND.AMGLSS.LT.
+ $ (AMB2SS+AMBT)) THEN
+ TMP(4)=AMB2SS
+ TMP(6)=AMT2SS
+ TMP(8)=SNW1
+ GT2B2=(-SINT*COSB*ART2*ALB2-COST*SINB*ALT2*ARB2)*
+ , SSXINT(AMTP,SSGWT6,UPPER)-
+ , (SINT*SINB*ALT2*ALB2+COST*COSB*ART2*ARB2)*
+ , SSXINT(AMTP,SSGWT4,UPPER)+
+ , (SINT*SINB*ALT2*ARB2+COST*COSB*ART2*ALB2)*
+ , SSXINT(AMTP,SSGWT5,UPPER)-
+ , (-SINT*COSB*ART2*ARB2-COST*SINB*ALT2*ALB2)*
+ , SSXINT(AMTP,SSGWT7,UPPER)
+ ELSE
+ GT2B2=0.
+ END IF
+ WID=GT1+GT2+GT1T2+GB1+GB2+GT1B1+GT1B2+GT2B1+GT2B2
+ WID=FACT*WID
+ IF (WID.GT.0.) THEN
+ CALL SSSAVE(ISGL,WID,+ISW1,+IDBT,-IDTP,0,0)
+ Z1(1)=SINT-COST
+ Z1(2)=SINT+COST
+ Z2(1)=ALT1+SNW1*ART1
+ Z2(2)=ALT1-SNW1*ART1
+ CALL SSME3(2,AMT1SS,Z1,Z2)
+ Z1(1)=-COST-SINT
+ Z1(2)=-COST+SINT
+ Z2(1)=ALT2+SNW1*ART2
+ Z2(2)=ALT2-SNW1*ART2
+ CALL SSME3(2,AMT2SS,Z1,Z2)
+ Z1(1)=SNW1*ALB1+ARB1
+ Z1(2)=-SNW1*ALB1+ARB1
+ Z2(1)=-COSB+SINB
+ Z2(2)=-COSB-SINB
+ CALL SSME3(3,AMB1SS,Z1,Z2)
+ Z1(1)=SNW1*ALB2+ARB2
+ Z1(2)=-SNW1*ALB2+ARB2
+ Z2(1)=-SINB-COSB
+ Z2(2)=-SINB+COSB
+ CALL SSME3(3,AMB2SS,Z1,Z2)
+ CALL SSSAVE(ISGL,WID,-ISW1,+IDTP,-IDBT,0,0)
+ Z1(1)=SINB-COSB
+ Z1(2)=SINB+COSB
+ Z2(1)=ALB1+SNW1*ARB1
+ Z2(2)=ALB1-SNW1*ARB1
+ CALL SSME3(2,AMB1SS,Z1,Z2)
+ Z1(1)=-COSB-SINB
+ Z1(2)=-COSB+SINB
+ Z2(1)=ALB2+SNW1*ARB2
+ Z2(2)=ALB2-SNW1*ARB2
+ CALL SSME3(2,AMB2SS,Z1,Z2)
+ Z1(1)=SNW1*ALT1+ART1
+ Z1(2)=-SNW1*ALT1+ART1
+ Z2(1)=-COST+SINT
+ Z2(2)=-COST-SINT
+ CALL SSME3(3,AMT1SS,Z1,Z2)
+ Z1(1)=SNW1*ALT2+ART2
+ Z1(2)=-SNW1*ALT2+ART2
+ Z2(1)=-SINT-COST
+ Z2(2)=-SINT+COST
+ CALL SSME3(3,AMT2SS,Z1,Z2)
+ END IF
+ END IF
+C
+C gluino --> w2 + qk + qb
+C
+ TMP(1)=MW2
+ UPPR=(AMGLSS**2-MW2**2)/2./AMGLSS
+ IF (AMGLSS.GT.(MW2+AMUP+AMDN)) THEN
+ IF (AMGLSS.LT.AMULSS.AND.AMGLSS.LT.AMDLSS) THEN
+ TMP(2)=AMULSS
+ WIDU=ZADW2*CONJG(ZADW2)*SSXINT(0.,SSGWQ1,UPPR)
+ TMP(2)=AMDLSS
+ WIDD=ZAUW2*CONJG(ZAUW2)*SSXINT(0.,SSGWQ1,UPPR)
+ WIDI=2*REAL(ZAUW2*ZADW2)*SSXINT(0.,SSGWQ2,UPPR)
+ WID=GS**2/2./AMGLSS/(2*PI)**5*(WIDU+WIDD+WIDI)
+ ELSE IF (AMGLSS.LT.AMULSS.AND.AMGLSS.GE.AMDLSS) THEN
+ TMP(2)=AMULSS
+ WIDU=ZADW2*CONJG(ZADW2)*SSXINT(0.,SSGWQ1,UPPR)
+ WID=GS**2/2./AMGLSS/(2*PI)**5*WIDU
+ ELSE IF (AMGLSS.LT.AMDLSS.AND.AMGLSS.GE.AMULSS) THEN
+ TMP(2)=AMDLSS
+ WIDD=ZAUW2*CONJG(ZAUW2)*SSXINT(0.,SSGWQ1,UPPR)
+ WID=GS**2/2./AMGLSS/(2*PI)**5*WIDD
+ ELSE
+ WID=0.
+ END IF
+ CALL SSSAVE(ISGL,WID,+ISW2,+IDDN,-IDUP,0,0)
+ CALL SSSAVE(ISGL,WID,-ISW2,+IDUP,-IDDN,0,0)
+ Z1(1)=1.
+ Z1(2)=-Z1(1)
+ Z2(1)=G*THX*COS(GAMMAL)
+ Z2(2)=Z2(1)
+ CALL SSME3(2,AMDLSS,Z1,Z2)
+ Z1(1)=G*THY*COS(GAMMAR)
+ Z1(2)=-Z1(1)
+ Z2(1)=1.
+ Z2(2)=Z2(1)
+ CALL SSME3(3,AMULSS,Z1,Z2)
+ END IF
+C
+ IF (AMGLSS.GT.(MW2+AMST+AMCH)) THEN
+ IF (AMGLSS.LT.AMCLSS.AND.AMGLSS.LT.AMSLSS) THEN
+ TMP(2)=AMCLSS
+ WIDU=ZADW2*CONJG(ZADW2)*SSXINT(0.,SSGWQ1,UPPR)
+ TMP(2)=AMSLSS
+ WIDD=ZAUW2*CONJG(ZAUW2)*SSXINT(0.,SSGWQ1,UPPR)
+ WIDI=2*REAL(ZAUW2*ZADW2)*SSXINT(0.,SSGWQ2,UPPR)
+ WID=GS**2/2./AMGLSS/(2*PI)**5*(WIDU+WIDD+WIDI)
+ ELSE IF (AMGLSS.LT.AMCLSS.AND.AMGLSS.GE.AMSLSS) THEN
+ TMP(2)=AMCLSS
+ WIDU=ZADW2*CONJG(ZADW2)*SSXINT(0.,SSGWQ1,UPPR)
+ WID=GS**2/2./AMGLSS/(2*PI)**5*WIDU
+ ELSE IF (AMGLSS.LT.AMSLSS.AND.AMGLSS.GE.AMCLSS) THEN
+ TMP(2)=AMSLSS
+ WIDD=ZAUW2*CONJG(ZAUW2)*SSXINT(0.,SSGWQ1,UPPR)
+ WID=GS**2/2./AMGLSS/(2*PI)**5*WIDD
+ ELSE
+ WID=0.
+ END IF
+ CALL SSSAVE(ISGL,WID,+ISW2,+IDST,-IDCH,0,0)
+ CALL SSSAVE(ISGL,WID,-ISW2,+IDCH,-IDST,0,0)
+ Z1(1)=1.
+ Z1(2)=-Z1(1)
+ Z2(1)=G*THX*COS(GAMMAL)
+ Z2(2)=Z2(1)
+ CALL SSME3(2,AMSLSS,Z1,Z2)
+ Z1(1)=G*THY*COS(GAMMAR)
+ Z1(2)=-Z1(1)
+ Z2(1)=1.
+ Z2(2)=Z2(1)
+ CALL SSME3(3,AMCLSS,Z1,Z2)
+ END IF
+C
+ IF (AMGLSS.GT.(MW2+AMTP+AMBT)) THEN
+ ALT1=-G*THY*COS(GAMMAR)*COST-FT*THY*SIN(GAMMAR)*SINT
+ ART1=FB*THX*SIN(GAMMAL)*COST
+ ALT2=-G*THY*COS(GAMMAR)*SINT+FT*THY*SIN(GAMMAR)*COST
+ ART2=FB*THX*SIN(GAMMAL)*SINT
+ ALB1=-G*THX*COS(GAMMAL)*COSB-FB*THX*SIN(GAMMAL)*SINB
+ ARB1=FT*THY*SIN(GAMMAR)*COSB
+ ALB2=-G*THX*COS(GAMMAL)*SINB+FB*THX*SIN(GAMMAL)*COSB
+ ARB2=FT*THY*SIN(GAMMAR)*SINB
+ UPPER=(AMGLSS**2+AMTP**2-(MW2+AMBT)**2)/2./AMGLSS
+ FACT=GS**2*PI**2/(2*PI)**5/2./AMGLSS
+ TMP(1)=MW2
+ TMP(2)=AMGLSS
+ TMP(3)=AMTP
+ IF (AMGLSS.LT.(AMTP+AMT1SS)) THEN
+ TMP(6)=AMT1SS
+ TMP(7)=AMT1SS
+ GT1=(ALT1**2+ART1**2)*(SSXINT(AMTP,SSGWT1,UPPER)+
+ , SIN(2*THETAT)*SSXINT(AMTP,SSGWT8,UPPER))
+ TMP(7)=AMT2SS
+ GT1T2=2*(ALT1*ALT2+ART1*ART2)*(SINT**2-COST**2)*
+ , SSXINT(AMTP,SSGWT8,UPPER)
+ ELSE
+ GT1=0.
+ GT1T2=0.
+ END IF
+ IF (AMGLSS.LT.(AMTP+AMT2SS)) THEN
+ TMP(6)=AMT2SS
+ TMP(7)=AMT2SS
+ GT2=(ALT2**2+ART2**2)*(SSXINT(AMTP,SSGWT1,UPPER)-
+ , SIN(2*THETAT)*SSXINT(AMTP,SSGWT8,UPPER))
+ ELSE
+ GT2=0.
+ END IF
+ IF (AMGLSS.LT.(AMB1SS+AMBT)) THEN
+ TMP(4)=AMB1SS
+ TMP(8)=SNW2
+C Rewrite UPPR=(AMGLSS**2+AMBT**2-(AMTP+MW2)**2)/2./AMGLSS
+ UPPR=((AMGLSS-AMTP-MW2)*(AMGLSS+AMTP+MW2)+AMBT**2)/(2*AMGLSS)
+ GB1=(ALB1**2+ARB1**2)*SSXINT(AMBT,SSGWT2,UPPR)-
+ , ALB1*ARB1*SSXINT(AMBT,SSGWT3,UPPR)
+ ELSE
+ GB1=0.
+ END IF
+ IF (AMGLSS.LT.(AMB2SS+AMBT)) THEN
+ TMP(4)=AMB2SS
+ TMP(8)=SNW2
+ UPPR=((AMGLSS-AMTP-MW2)*(AMGLSS+AMTP+MW2)+AMBT**2)/(2*AMGLSS)
+ GB2=(ALB2**2+ARB2**2)*SSXINT(AMBT,SSGWT2,UPPR)-
+ , ALB2*ARB2*SSXINT(AMBT,SSGWT3,UPPR)
+ ELSE
+ GB2=0.
+ END IF
+ IF (AMGLSS.LT.(AMTP+AMT1SS).AND.AMGLSS.LT.
+ $ (AMB1SS+AMBT)) THEN
+ TMP(4)=AMB1SS
+ TMP(6)=AMT1SS
+ TMP(8)=SNW2
+ GT1B1=(COST*SINB*ART1*ALB1+SINT*COSB*ALT1*ARB1)*
+ , SSXINT(AMTP,SSGWT6,UPPER)-
+ , (COST*COSB*ALT1*ALB1+SINT*SINB*ART1*ARB1)*
+ , SSXINT(AMTP,SSGWT4,UPPER)+
+ , (COST*COSB*ALT1*ARB1+SINT*SINB*ART1*ALB1)*
+ , SSXINT(AMTP,SSGWT5,UPPER)-
+ , (COST*SINB*ART1*ARB1+SINT*COSB*ALT1*ALB1)*
+ , SSXINT(AMTP,SSGWT7,UPPER)
+ ELSE
+ GT1B1=0.
+ END IF
+ IF (AMGLSS.LT.(AMTP+AMT1SS).AND.AMGLSS.LT.
+ $ (AMB2SS+AMBT)) THEN
+ TMP(4)=AMB2SS
+ TMP(6)=AMT1SS
+ TMP(8)=SNW2
+ GT1B2=(-COST*COSB*ART1*ALB2+SINT*SINB*ALT1*ARB2)*
+ , SSXINT(AMTP,SSGWT6,UPPER)-
+ , (COST*COSB*ALT1*ALB2-SINT*COSB*ART1*ARB2)*
+ , SSXINT(AMTP,SSGWT4,UPPER)+
+ , (COST*SINB*ALT1*ARB2-SINT*COSB*ART1*ALB2)*
+ , SSXINT(AMTP,SSGWT5,UPPER)-
+ , (-COST*COSB*ART1*ARB2+SINT*SINB*ALT1*ALB2)*
+ , SSXINT(AMTP,SSGWT7,UPPER)
+ ELSE
+ GT1B2=0.
+ END IF
+ IF (AMGLSS.LT.(AMTP+AMT2SS).AND.AMGLSS.LT.
+ $ (AMB1SS+AMBT)) THEN
+ TMP(4)=AMB1SS
+ TMP(6)=AMT2SS
+ TMP(8)=SNW2
+ GT2B1=(SINT*SINB*ART2*ALB1-COST*COSB*ALT2*ARB1)*
+ , SSXINT(AMTP,SSGWT6,UPPER)-
+ , (SINT*COSB*ALT2*ALB1-COST*SINB*ART2*ARB1)*
+ , SSXINT(AMTP,SSGWT4,UPPER)+
+ , (SINT*COSB*ALT2*ARB1-COST*SINB*ART2*ALB1)*
+ , SSXINT(AMTP,SSGWT5,UPPER)-
+ , (SINT*SINB*ART2*ARB1-COST*COSB*ALT2*ALB1)*
+ , SSXINT(AMTP,SSGWT7,UPPER)
+ ELSE
+ GT2B1=0.
+ END IF
+ IF (AMGLSS.LT.(AMTP+AMT2SS).AND.AMGLSS.LT.
+ $ (AMB2SS+AMBT)) THEN
+ TMP(4)=AMB2SS
+ TMP(6)=AMT2SS
+ TMP(8)=SNW2
+ GT2B2=(-SINT*COSB*ART2*ALB2-COST*SINB*ALT2*ARB2)*
+ , SSXINT(AMTP,SSGWT6,UPPER)-
+ , (SINT*SINB*ALT2*ALB2+COST*COSB*ART2*ARB2)*
+ , SSXINT(AMTP,SSGWT4,UPPER)+
+ , (SINT*SINB*ALT2*ARB2+COST*COSB*ART2*ALB2)*
+ , SSXINT(AMTP,SSGWT5,UPPER)-
+ , (-SINT*COSB*ART2*ARB2-COST*SINB*ALT2*ALB2)*
+ , SSXINT(AMTP,SSGWT7,UPPER)
+ ELSE
+ GT2B2=0.
+ END IF
+ WID=GT1+GT2+GT1T2+GB1+GB2+GT1B1+GT1B2+GT2B1+GT2B2
+ WID=FACT*WID
+ IF (WID.GT.0.) THEN
+ CALL SSSAVE(ISGL,WID,+ISW2,+IDBT,-IDTP,0,0)
+ Z1(1)=SINT-COST
+ Z1(2)=SINT+COST
+ Z2(1)=ALT1+SNW2*ART1
+ Z2(2)=ALT1-SNW2*ART1
+ CALL SSME3(2,AMT1SS,Z1,Z2)
+ Z1(1)=-COST-SINT
+ Z1(2)=-COST+SINT
+ Z2(1)=ALT2+SNW2*ART2
+ Z2(2)=ALT2-SNW2*ART2
+ CALL SSME3(2,AMT2SS,Z1,Z2)
+ Z1(1)=SNW2*ALB1+ARB1
+ Z1(2)=-SNW2*ALB1+ARB1
+ Z2(1)=-COSB+SINB
+ Z2(2)=-COSB-SINB
+ CALL SSME3(3,AMB1SS,Z1,Z2)
+ Z1(1)=SNW2*ALB2+ARB2
+ Z1(2)=-SNW2*ALB2+ARB2
+ Z2(1)=-SINB-COSB
+ Z2(2)=-SINB+COSB
+ CALL SSME3(3,AMB2SS,Z1,Z2)
+ CALL SSSAVE(ISGL,WID,-ISW2,+IDTP,-IDBT,0,0)
+ Z1(1)=SINB-COSB
+ Z1(2)=SINB+COSB
+ Z2(1)=ALB1+SNW2*ARB1
+ Z2(2)=ALB1-SNW2*ARB1
+ CALL SSME3(2,AMB1SS,Z1,Z2)
+ Z1(1)=-COSB-SINB
+ Z1(2)=-COSB+SINB
+ Z2(1)=ALB2+SNW2*ARB2
+ Z2(2)=ALB2-SNW2*ARB2
+ CALL SSME3(2,AMB2SS,Z1,Z2)
+ Z1(1)=SNW2*ALT1+ART1
+ Z1(2)=-SNW2*ALT1+ART1
+ Z2(1)=-COST+SINT
+ Z2(2)=-COST-SINT
+ CALL SSME3(3,AMT1SS,Z1,Z2)
+ Z1(1)=SNW2*ALT2+ART2
+ Z1(2)=-SNW2*ALT2+ART2
+ Z2(1)=-SINT-COST
+ Z2(2)=-SINT+COST
+ CALL SSME3(3,AMT2SS,Z1,Z2)
+ END IF
+ END IF
+C
+C gluino --> zi decays, zi = z1, z2, z3, z4
+C the auiz etc, below are Atilde's etc. of PRD 42,1568 (1990)
+C
+ ISZI(1)=ISZ1
+ ISZI(2)=ISZ2
+ ISZI(3)=ISZ3
+ ISZI(4)=ISZ4
+ DO 100 IZ=1,4
+ MZIZ=ABS(AMZISS(IZ))
+ AUIZ=G/SR2*ZMIXSS(3,IZ)+GP/3./SR2*ZMIXSS(4,IZ)
+ ADIZ=-G/SR2*ZMIXSS(3,IZ)+GP/3./SR2*ZMIXSS(4,IZ)
+ BUIZ=4*GP*ZMIXSS(4,IZ)/3./SR2
+ BDIZ=-2*GP/3./SR2*ZMIXSS(4,IZ)
+ SNIZ=SIGN(1.,AMZISS(IZ))
+ THIZ=0
+ IF (AMZISS(IZ).LT.0.) THIZ=1
+ ZAUIZ=ZI**(THIZ-1)*SNIZ
+ $ *(-G/SR2*ZMIXSS(3,IZ)-GP/3./SR2*ZMIXSS(4,IZ))
+ ZBUIZ=ZI**(THIZ-1)*4*GP*ZMIXSS(4,IZ)/3./SR2
+ ZADIZ=ZI**(THIZ-1)*SNIZ
+ $ *(G/SR2*ZMIXSS(3,IZ)-GP/3./SR2*ZMIXSS(4,IZ))
+ ZBDIZ=-2*ZI**(THIZ-1)*GP*ZMIXSS(4,IZ)/3./SR2
+C Radiative gluino --> gluon + zi loop decay
+ IF (AMGLSS.GT.MZIZ) THEN
+ IF (AMGLSS.LT.(AMT1SS+AMTP)) THEN
+ TMP(1)=AMTP
+ TMP(2)=MZIZ
+ TMP(3)=AMT1SS
+ XIT1=SSXINT(0.,SSGZG1,1.)
+ XI1T1=SSXINT(0.,SSGZG2,1.)
+ XKT1=SSXINT(0.,SSGZG3,1.)
+ ELSE
+ XIT1=0.
+ XI1T1=0.
+ XKT1=0.
+ END IF
+ IF (AMGLSS.LT.(AMT2SS+AMTP)) THEN
+ TMP(1)=AMTP
+ TMP(2)=MZIZ
+ TMP(3)=AMT2SS
+ XIT2=SSXINT(0.,SSGZG1,1.)
+ XI1T2=SSXINT(0.,SSGZG2,1.)
+ XKT2=SSXINT(0.,SSGZG3,1.)
+ ELSE
+ XIT2=0.
+ XI1T2=0.
+ XKT2=0.
+ END IF
+C !!! NEEDS UPDATE FOR SBOTTOM MIXING !!!
+ IF (AMGLSS.LT.(AMB1SS+AMBT)) THEN
+ TMP(1)=AMBT
+ TMP(2)=MZIZ
+ TMP(3)=AMB1SS
+ XI1BL=SSXINT(0.,SSGZG2,1.)
+ XKBL=SSXINT(0.,SSGZG3,1.)
+ ELSE
+ XI1BL=0.
+ XKBL=0.
+ END IF
+ IF (AMGLSS.LT.(AMB2SS+AMBT)) THEN
+ TMP(1)=AMBT
+ TMP(2)=MZIZ
+ TMP(3)=AMB2SS
+ XI1BR=SSXINT(0.,SSGZG2,1.)
+ XKBR=SSXINT(0.,SSGZG3,1.)
+ ELSE
+ XI1BR=0.
+ XKBR=0.
+ END IF
+ IF (AMGLSS.LT.(AMULSS+AMUP)) THEN
+ TMP(1)=AMUP
+ TMP(2)=MZIZ
+ TMP(3)=AMULSS
+ XI1UL=SSXINT(0.,SSGZG2,1.)
+ XKUL=SSXINT(0.,SSGZG3,1.)
+ ELSE
+ XI1UL=0.
+ XKUL=0.
+ END IF
+ IF (AMGLSS.LT.(AMURSS+AMUP)) THEN
+ TMP(1)=AMUP
+ TMP(2)=MZIZ
+ TMP(3)=AMURSS
+ XI1UR=SSXINT(0.,SSGZG2,1.)
+ XKUR=SSXINT(0.,SSGZG3,1.)
+ ELSE
+ XI1UR=0.
+ XKUR=0.
+ END IF
+ IF (AMGLSS.LT.(AMDLSS+AMDN)) THEN
+ TMP(1)=AMDN
+ TMP(2)=MZIZ
+ TMP(3)=AMDLSS
+ XI1DL=SSXINT(0.,SSGZG2,1.)
+ XKDL=SSXINT(0.,SSGZG3,1.)
+ ELSE
+ XI1DL=0.
+ XKDL=0.
+ END IF
+ IF (AMGLSS.LT.(AMDRSS+AMDN)) THEN
+ TMP(1)=AMDN
+ TMP(2)=MZIZ
+ TMP(3)=AMDRSS
+ XI1DR=SSXINT(0.,SSGZG2,1.)
+ XKDR=SSXINT(0.,SSGZG3,1.)
+ ELSE
+ XI1DR=0.
+ XKDR=0.
+ END IF
+ IF (AMGLSS.LT.(AMSLSS+AMST)) THEN
+ TMP(1)=AMST
+ TMP(2)=MZIZ
+ TMP(3)=AMSLSS
+ XI1SL=SSXINT(0.,SSGZG2,1.)
+ XKSL=SSXINT(0.,SSGZG3,1.)
+ ELSE
+ XI1SL=0.
+ XKSL=0.
+ END IF
+ IF (AMGLSS.LT.(AMSRSS+AMST)) THEN
+ TMP(1)=AMST
+ TMP(2)=MZIZ
+ TMP(3)=AMSRSS
+ XI1SR=SSXINT(0.,SSGZG2,1.)
+ XKSR=SSXINT(0.,SSGZG3,1.)
+ ELSE
+ XI1SR=0.
+ XKSR=0.
+ END IF
+ IF (AMGLSS.LT.(AMCLSS+AMCH)) THEN
+ TMP(1)=AMCH
+ TMP(2)=MZIZ
+ TMP(3)=AMCLSS
+ XI1CL=SSXINT(0.,SSGZG2,1.)
+ XKCL=SSXINT(0.,SSGZG3,1.)
+ ELSE
+ XI1CL=0.
+ XKCL=0.
+ END IF
+ IF (AMGLSS.LT.(AMCRSS+AMCH)) THEN
+ TMP(1)=AMCH
+ TMP(2)=MZIZ
+ TMP(3)=AMCRSS
+ XI1CR=SSXINT(0.,SSGZG2,1.)
+ XKCR=SSXINT(0.,SSGZG3,1.)
+ ELSE
+ XI1CR=0.
+ XKCR=0.
+ END IF
+ KUL=AUIZ*(XKUL*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1UL)
+ KUR=-BUIZ*(XKUR*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1UR)
+ KDL=ADIZ*(XKDL*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1DL)
+ KDR=-BDIZ*(XKDR*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1DR)
+ KCL=AUIZ*(XKCL*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1CL)
+ KCR=-BUIZ*(XKCR*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1CR)
+ KSL=ADIZ*(XKSL*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1SL)
+ KSR=-BDIZ*(XKSR*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1SR)
+ KBL=ADIZ*(XKBL*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1BL)
+ KBR=-BDIZ*(XKBR*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1BR)
+ KT1=(MZIZ*(XKT1+XI1T1)*(AUIZ*COST-FT*ZMIXSS(1,IZ)*SINT)
+ $ -SNIZ*AMGLSS*XKT1*(AUIZ*COST-FT*ZMIXSS(1,IZ)*SINT)+
+ $ SNIZ*AMTP*XIT1*(BUIZ*SINT+FT*ZMIXSS(1,IZ)*COST))*COST+
+ $ (MZIZ*(XKT1+XI1T1)*(-BUIZ*SINT-FT*ZMIXSS(1,IZ)*COST)
+ $ -SNIZ*AMGLSS*XKT1*(-BUIZ*SINT-FT*ZMIXSS(1,IZ)*COST)-
+ $ SNIZ*AMTP*XIT1*(AUIZ*COST-FT*ZMIXSS(1,IZ)*SINT))*SINT
+ KT2=(MZIZ*(XKT2+XI1T2)*(AUIZ*SINT+FT*ZMIXSS(1,IZ)*COST)
+ $ -SNIZ*AMGLSS*XKT2*(AUIZ*SINT+FT*ZMIXSS(1,IZ)*COST)+
+ $ SNIZ*AMTP*XIT2*(-BUIZ*COST+FT*ZMIXSS(1,IZ)*SINT))*SINT+
+ $ (-MZIZ*(XKT2+XI1T2)*(BUIZ*COST-FT*ZMIXSS(1,IZ)*SINT)
+ $ +SNIZ*AMGLSS*XKT2*(BUIZ*COST-FT*ZMIXSS(1,IZ)*SINT)+
+ $ SNIZ*AMTP*XIT2*(AUIZ*SINT+FT*ZMIXSS(1,IZ)*COST))*COST
+ SUM=(KUL+KUR+KDL+KDR+KSL+KSR+KCL+KCR+KBL+KBR+KT1+
+ $ KT2)**2/AMGLSS**2
+ WID=ALF3**2*AMGLSS*(1.-MZIZ**2/AMGLSS**2)/256./PI**3*SUM
+ CALL SSSAVE(ISGL,WID,ISZI(IZ),IDGL,0,0,0)
+ END IF
+C 3 body gluino --> q + qb + zi decay, q=u,d
+ UPPR=(AMGLSS**2-MZIZ**2)/2./AMGLSS
+ IF (AMGLSS.GT.(MZIZ+2*AMUP)) THEN
+ IF (AMGLSS.LT.AMULSS.AND.AMGLSS.LT.AMURSS) THEN
+ TMP(1)=MZIZ
+ TMP(2)=AMULSS
+ WIDL=2*AUIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
+ $ SSXINT(0.,SSGWQ2,UPPR))
+ TMP(2)=AMURSS
+ WIDR=2*BUIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
+ $ SSXINT(0.,SSGWQ2,UPPR))
+ WID=WIDL+WIDR
+ ELSE IF (AMGLSS.LT.AMULSS.AND.AMGLSS.GE.AMURSS) THEN
+ TMP(2)=AMULSS
+ WID=2*AUIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
+ $ SSXINT(0.,SSGWQ2,UPPR))
+ ELSE IF (AMGLSS.LT.AMURSS.AND.AMGLSS.GE.AMULSS) THEN
+ TMP(2)=AMURSS
+ WID=2*BUIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
+ $ SSXINT(0.,SSGWQ2,UPPR))
+ ELSE
+ WID=0.
+ END IF
+ WID=GS**2/AMGLSS/2./(2*PI)**5*WID
+ IF (WID.GT.0.) THEN
+ CALL SSSAVE(ISGL,WID,ISZI(IZ),IDUP,-IDUP,0,0)
+C Enter decay matrix element info
+ Z1(1)=1.
+ Z1(2)=-Z1(1)
+ Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)*(THIZ+1)*AUIZ)
+ Z2(2)=Z2(1)
+ CALL SSME3(2,AMULSS,Z1,Z2)
+ Z1(1)=1.
+ Z1(2)=Z1(1)
+ Z2(1)=-CONJG(ZI**(THIZ-1)*BUIZ)
+ Z2(2)=-Z2(1)
+ CALL SSME3(2,AMURSS,Z1,Z2)
+ Z1(1)=ZI**(THIZ-1)*(-1.)*(THIZ+1)*AUIZ
+ Z1(2)=-Z1(1)
+ Z2(1)=1.
+ Z2(2)=Z2(1)
+ CALL SSME3(3,AMULSS,Z1,Z2)
+ Z1(1)=ZI**(THIZ-1)*BUIZ
+ Z1(2)=Z1(1)
+ Z2(1)=1.
+ Z2(2)=-Z2(1)
+ CALL SSME3(3,AMURSS,Z1,Z2)
+ END IF
+ END IF
+C
+ IF (AMGLSS.GT.(MZIZ+2*AMDN)) THEN
+ IF (AMGLSS.LT.AMDLSS.AND.AMGLSS.LT.AMDRSS) THEN
+ TMP(1)=MZIZ
+ TMP(2)=AMDLSS
+ WIDL=2*ADIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
+ $ SSXINT(0.,SSGWQ2,UPPR))
+ TMP(2)=AMDRSS
+ WIDR=2*BDIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
+ $ SSXINT(0.,SSGWQ2,UPPR))
+ WID=WIDL+WIDR
+ ELSE IF (AMGLSS.LT.AMDLSS.AND.AMGLSS.GE.AMDRSS) THEN
+ TMP(2)=AMDLSS
+ WID=2*ADIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
+ $ SSXINT(0.,SSGWQ2,UPPR))
+ ELSE IF (AMGLSS.LT.AMDRSS.AND.AMGLSS.GE.AMDLSS) THEN
+ TMP(2)=AMDRSS
+ WID=2*BDIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
+ $ SSXINT(0.,SSGWQ2,UPPR))
+ ELSE
+ WID=0.
+ END IF
+ WID=GS**2/AMGLSS/2./(2*PI)**5*WID
+ IF (WID.GT.0.) THEN
+ CALL SSSAVE(ISGL,WID,ISZI(IZ),IDDN,-IDDN,0,0)
+C Enter decay matrix element info
+ Z1(1)=1.
+ Z1(2)=-Z1(1)
+ Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)*(THIZ+1)*ADIZ)
+ Z2(2)=Z2(1)
+ CALL SSME3(2,AMDLSS,Z1,Z2)
+ Z1(1)=1.
+ Z1(2)=Z1(1)
+ Z2(1)=-CONJG(ZI**(THIZ-1)*BDIZ)
+ Z2(2)=-Z2(1)
+ CALL SSME3(2,AMDRSS,Z1,Z2)
+ Z1(1)=ZI**(THIZ-1)*(-1.)*(THIZ+1)*ADIZ
+ Z1(2)=-Z1(1)
+ Z2(1)=1.
+ Z2(2)=Z2(1)
+ CALL SSME3(3,AMDLSS,Z1,Z2)
+ Z1(1)=ZI**(THIZ-1)*BDIZ
+ Z1(2)=Z1(1)
+ Z2(1)=1.
+ Z2(2)=-Z2(1)
+ CALL SSME3(3,AMDRSS,Z1,Z2)
+ END IF
+ END IF
+C 3 body gluino --> q + qb + zi decay, q=s
+ IF (AMGLSS.GT.(MZIZ+2*AMST)) THEN
+ IF (AMGLSS.LT.AMSLSS.AND.AMGLSS.LT.AMSRSS) THEN
+ TMP(1)=MZIZ
+ TMP(2)=AMSLSS
+ WIDL=2*ADIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
+ $ SSXINT(0.,SSGWQ2,UPPR))
+ TMP(2)=AMSRSS
+ WIDR=2*BDIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
+ $ SSXINT(0.,SSGWQ2,UPPR))
+ WID=WIDL+WIDR
+ ELSE IF (AMGLSS.LT.AMSLSS.AND.AMGLSS.GE.AMSRSS) THEN
+ TMP(2)=AMSLSS
+ WID=2*ADIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
+ $ SSXINT(0.,SSGWQ2,UPPR))
+ ELSE IF (AMGLSS.LT.AMSRSS.AND.AMGLSS.GE.AMSLSS) THEN
+ TMP(2)=AMSRSS
+ WID=2*BDIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
+ $ SSXINT(0.,SSGWQ2,UPPR))
+ ELSE
+ WID=0.
+ END IF
+ WID=GS**2/AMGLSS/2./(2*PI)**5*WID
+ IF (WID.GT.0.) THEN
+ CALL SSSAVE(ISGL,WID,ISZI(IZ),IDST,-IDST,0,0)
+C Enter decay matrix element info
+ Z1(1)=1.
+ Z1(2)=-Z1(1)
+ Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)*(THIZ+1)*ADIZ)
+ Z2(2)=Z2(1)
+ CALL SSME3(2,AMDLSS,Z1,Z2)
+ Z1(1)=1.
+ Z1(2)=Z1(1)
+ Z2(1)=-CONJG(ZI**(THIZ-1)*BDIZ)
+ Z2(2)=-Z2(1)
+ CALL SSME3(2,AMDRSS,Z1,Z2)
+ Z1(1)=ZI**(THIZ-1)*(-1.)*(THIZ+1)*ADIZ
+ Z1(2)=-Z1(1)
+ Z2(1)=1.
+ Z2(2)=Z2(1)
+ CALL SSME3(3,AMDLSS,Z1,Z2)
+ Z1(1)=ZI**(THIZ-1)*BDIZ
+ Z1(2)=Z1(1)
+ Z2(1)=1.
+ Z2(2)=-Z2(1)
+ CALL SSME3(3,AMDRSS,Z1,Z2)
+ END IF
+ END IF
+C 3 body gluino --> q + qb + zi decay, q=c
+ IF (AMGLSS.GT.(MZIZ+2*AMCH)) THEN
+ IF (AMGLSS.LT.AMCLSS.AND.AMGLSS.LT.AMCRSS) THEN
+ TMP(1)=MZIZ
+ TMP(2)=AMCLSS
+ WIDL=2*AUIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
+ $ SSXINT(0.,SSGWQ2,UPPR))
+ TMP(2)=AMCRSS
+ WIDR=2*BUIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
+ $ SSXINT(0.,SSGWQ2,UPPR))
+ WID=WIDL+WIDR
+ ELSE IF (AMGLSS.LT.AMCLSS.AND.AMGLSS.GE.AMCRSS) THEN
+ TMP(2)=AMCLSS
+ WID=2*AUIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
+ $ SSXINT(0.,SSGWQ2,UPPR))
+ ELSE IF (AMGLSS.LT.AMCRSS.AND.AMGLSS.GE.AMCLSS) THEN
+ TMP(2)=AMCRSS
+ WID=2*BUIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
+ $ SSXINT(0.,SSGWQ2,UPPR))
+ ELSE
+ WID=0.
+ END IF
+ WID=GS**2/AMGLSS/2./(2*PI)**5*WID
+ IF (WID.GT.0.) THEN
+ CALL SSSAVE(ISGL,WID,ISZI(IZ),IDCH,-IDCH,0,0)
+C Enter decay matrix element info
+ Z1(1)=1.
+ Z1(2)=-Z1(1)
+ Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)*(THIZ+1)*AUIZ)
+ Z2(2)=Z2(1)
+ CALL SSME3(2,AMULSS,Z1,Z2)
+ Z1(1)=1.
+ Z1(2)=Z1(1)
+ Z2(1)=-CONJG(ZI**(THIZ-1)*BUIZ)
+ Z2(2)=-Z2(1)
+ CALL SSME3(2,AMURSS,Z1,Z2)
+ Z1(1)=ZI**(THIZ-1)*(-1.)*(THIZ+1)*AUIZ
+ Z1(2)=-Z1(1)
+ Z2(1)=1.
+ Z2(2)=Z2(1)
+ CALL SSME3(3,AMULSS,Z1,Z2)
+ Z1(1)=ZI**(THIZ-1)*BUIZ
+ Z1(2)=Z1(1)
+ Z2(1)=1.
+ Z2(2)=-Z2(1)
+ CALL SSME3(3,AMURSS,Z1,Z2)
+ END IF
+ END IF
+C 3 body gluino --> q + qb + zi decay, q=b
+ XTCON=ALF3/8./PI**4/AMGLSS
+ IF (AMGLSS.GT.(MZIZ+2*AMBT).AND.AMGLSS.LT.
+ $ (AMB1SS+AMBT)) THEN
+ TMP(1)=AMGLSS
+ TMP(2)=AMBT
+ TMP(3)=MZIZ
+ TMP(4)=AMB1SS
+ TMP(5)=AMB1SS
+ XUPPER=(AMGLSS**2+AMBT**2-(AMBT+MZIZ)**2)/2./AMGLSS
+ XT1=SSXINT(AMBT,SSGX1,XUPPER)
+ XT2=SSXINT(AMBT,SSGX2,XUPPER)
+ XT3=SSXINT(AMBT,SSGX3,XUPPER)
+ XT4=SSXINT(AMBT,SSGX4,XUPPER)
+ XT5=SSXINT(AMBT,SSGX5,XUPPER)
+ XT6=SSXINT(AMBT,SSGX6,XUPPER)
+ XT7=SSXINT(AMBT,SSGX7,XUPPER)
+ XT8=SSXINT(AMBT,SSGX8,XUPPER)
+ XT9=SSXINT(AMBT,SSGX9,XUPPER)
+ AL1=ADIZ*COSB-FB*ZMIXSS(2,IZ)*SINB
+ BE1=FB*ZMIXSS(2,IZ)*COSB+BDIZ*SINB
+C ---- here, al2 is (-) al2 of tata notes-----
+ AL2=BDIZ*SINB+FB*ZMIXSS(2,IZ)*COSB
+ BE2=-FB*ZMIXSS(2,IZ)*SINB+ADIZ*COSB
+ XLL=(AL1**2+BE1**2)*XT1-4*AMBT*MZIZ*SNIZ*AL1*
+ $ BE1*XT3-AMGLSS*(SNIZ*MZIZ*(AL1**2*XT2/AMGLSS/
+ $ MZIZ+BE1**2*AMBT**2*XT5)-AL1*BE1*AMBT*(XT4-
+ $ MZIZ**2*XT5))
+ XRR=(AL2**2+BE2**2)*XT1-4*AMBT*MZIZ*SNIZ*AL2*
+ $ BE2*XT3-AMGLSS*(SNIZ*MZIZ*(AL2**2*XT2/AMGLSS/
+ $ MZIZ+BE2**2*AMBT**2*XT5)-AL2*BE2*AMBT*(XT4-
+ $ MZIZ**2*XT5))
+ XL1R1=-2*AMGLSS*AMBT*((AL1*AL2+BE1*BE2)*SNIZ*AMBT*
+ $ MZIZ*XT6-(AL2*BE1+AL1*BE2)*XT7)
+ XL2R2=XL1R1
+ XL1R2=BE1*BE2*XT8+AL1*AL2*AMBT**2*XT4-AMBT*MZIZ*
+ $ SNIZ*(AL1*BE2+AL2*BE1)*XT9
+ XL2R1=XL1R2
+ XMSB1=COSB**2*XLL+SINB**2*XRR-SINB*COSB*(XL1R1+XL1R2+
+ $ XL2R1+XL2R2)
+ ELSE
+ XMSB1=0.
+ END IF
+ IF (AMGLSS.GT.(MZIZ+2*AMBT).AND.AMGLSS.LT.
+ $ (AMB2SS+AMBT)) THEN
+ TMP(1)=AMGLSS
+ TMP(2)=AMBT
+ TMP(3)=MZIZ
+ TMP(4)=AMB2SS
+ TMP(5)=AMB2SS
+ XUPPER=(AMGLSS**2+AMBT**2-(AMBT+MZIZ)**2)/2./AMGLSS
+ XT1=SSXINT(AMBT,SSGX1,XUPPER)
+ XT2=SSXINT(AMBT,SSGX2,XUPPER)
+ XT3=SSXINT(AMBT,SSGX3,XUPPER)
+ XT4=SSXINT(AMBT,SSGX4,XUPPER)
+ XT5=SSXINT(AMBT,SSGX5,XUPPER)
+ XT6=SSXINT(AMBT,SSGX6,XUPPER)
+ XT7=SSXINT(AMBT,SSGX7,XUPPER)
+ XT8=SSXINT(AMBT,SSGX8,XUPPER)
+ XT9=SSXINT(AMBT,SSGX9,XUPPER)
+ AL1=ADIZ*SINB+FB*ZMIXSS(2,IZ)*COSB
+ BE1=FB*ZMIXSS(2,IZ)*SINB-BDIZ*COSB
+C ---- here, al2 is (-) al2 of tata notes-----
+ AL2=-BDIZ*COSB+FB*ZMIXSS(2,IZ)*SINB
+ BE2=FB*ZMIXSS(2,IZ)*COSB+ADIZ*SINB
+ XLL=(AL1**2+BE1**2)*XT1-4*AMBT*MZIZ*SNIZ*AL1*
+ $ BE1*XT3-AMGLSS*(SNIZ*MZIZ*(AL1**2*XT2/AMGLSS/
+ $ MZIZ+BE1**2*AMBT**2*XT5)-AL1*BE1*AMBT*(XT4-
+ $ MZIZ**2*XT5))
+ XRR=(AL2**2+BE2**2)*XT1-4*AMBT*MZIZ*SNIZ*AL2*
+ $ BE2*XT3-AMGLSS*(SNIZ*MZIZ*(AL2**2*XT2/AMGLSS/
+ $ MZIZ+BE2**2*AMBT**2*XT5)-AL2*BE2*AMBT*(XT4-
+ $ MZIZ**2*XT5))
+ XL1R1=-2*AMGLSS*AMBT*((AL1*AL2+BE1*BE2)*SNIZ*AMBT*
+ $ MZIZ*XT6-(AL2*BE1+AL1*BE2)*XT7)
+ XL2R2=XL1R1
+ XL1R2=BE1*BE2*XT8+AL1*AL2*AMBT**2*XT4-AMBT*MZIZ*
+ $ SNIZ*(AL1*BE2+AL2*BE1)*XT9
+ XL2R1=XL1R2
+ XMSB2=SINB**2*XLL+COSB**2*XRR+SINB*COSB*(XL1R1+XL1R2+
+ $ XL2R1+XL2R2)
+ ELSE
+ XMSB2=0.
+ END IF
+C ----cross term between b_1 and b_2 graphs -----------
+ IF (AMGLSS.GT.(MZIZ+2*AMBT).AND.AMGLSS.LT.
+ $ (AMB1SS+AMBT)) THEN
+ TMP(1)=AMGLSS
+ TMP(2)=AMBT
+ TMP(3)=MZIZ
+ TMP(4)=AMB1SS
+ TMP(5)=AMB2SS
+ XUPPER=(AMGLSS**2+AMBT**2-(AMBT+MZIZ)**2)/2./AMGLSS
+ XT1=SSXINT(AMBT,SSGX1,XUPPER)
+ XT2=SSXINT(AMBT,SSGX2,XUPPER)
+ XT3=SSXINT(AMBT,SSGX3,XUPPER)
+ XT4=SSXINT(AMBT,SSGX4,XUPPER)
+ XT5=SSXINT(AMBT,SSGX5,XUPPER)
+ XT6=SSXINT(AMBT,SSGX6,XUPPER)
+ XT7=SSXINT(AMBT,SSGX7,XUPPER)
+ XT8=SSXINT(AMBT,SSGX8,XUPPER)
+ XT9=SSXINT(AMBT,SSGX9,XUPPER)
+ AL1=ADIZ*COSB-FB*ZMIXSS(2,IZ)*SINB
+ AL2=ADIZ*SINB+FB*ZMIXSS(2,IZ)*COSB
+ BE1=FB*ZMIXSS(2,IZ)*COSB+BDIZ*SINB
+ BE2=FB*ZMIXSS(2,IZ)*SINB-BDIZ*COSB
+ XL1L2=COSB*SINB*(2*(AL1*AL2+BE1*BE2)*XT1-4*SNIZ*AMBT*
+ $ MZIZ*(AL1*BE2+AL2*BE1)*XT3-AMGLSS*(2*MZIZ*
+ $ SNIZ*(AL1*AL2*XT2/AMGLSS/MZIZ+BE1*BE2*AMBT**2*
+ $ XT5)-(AL1*BE2+AL2*BE1)*AMBT*(XT4-MZIZ**2*XT5)))
+ AL1=-BDIZ*SINB-FB*ZMIXSS(2,IZ)*COSB
+ AL2=BDIZ*COSB-FB*ZMIXSS(2,IZ)*SINB
+ BE1=-FB*ZMIXSS(2,IZ)*SINB+ADIZ*COSB
+ BE2=FB*ZMIXSS(2,IZ)*COSB+ADIZ*SINB
+ XR1R2=-COSB*SINB*(2*(AL1*AL2+BE1*BE2)*XT1+4*SNIZ*AMBT*
+ $ MZIZ*(AL1*BE2+AL2*BE1)*XT3-AMGLSS*(2*MZIZ*
+ $ SNIZ*(AL1*AL2*XT2/AMGLSS/MZIZ+BE1*BE2*AMBT**2*
+ $ XT5)+(AL1*BE2+AL2*BE1)*AMBT*(XT4-MZIZ**2*XT5)))
+ AL1=ADIZ*COSB-FB*ZMIXSS(2,IZ)*SINB
+ AL2=BDIZ*COSB-FB*ZMIXSS(2,IZ)*SINB
+ BE1=FB*ZMIXSS(2,IZ)*COSB+BDIZ*SINB
+ BE2=FB*ZMIXSS(2,IZ)*COSB+ADIZ*SINB
+ XL1R1=2*AMGLSS*AMBT*COSB**2*(SNIZ*(AL1*AL2-BE1*BE2)*
+ $ AMBT*MZIZ*XT6-(AL2*BE1-AL1*BE2)*XT7)
+ XL1R2=COSB**2*(BE1*BE2*XT8-AL1*AL2*AMBT**2*XT4+AMBT*
+ $ MZIZ*SNIZ*XT9*(-AL1*BE2+BE1*AL2))
+ XLR1=2*(XL1R1+XL1R2)
+ AL1=ADIZ*SINB+FB*ZMIXSS(2,IZ)*COSB
+ AL2=-BDIZ*SINB-FB*ZMIXSS(2,IZ)*COSB
+ BE1=FB*ZMIXSS(2,IZ)*SINB-BDIZ*COSB
+ BE2=-FB*ZMIXSS(2,IZ)*SINB+ADIZ*COSB
+ TMP(4)=AMB2SS
+ TMP(5)=AMB1SS
+ XT8=SSXINT(AMBT,SSGX8,XUPPER)
+ XT9=SSXINT(AMBT,SSGX9,XUPPER)
+ XL1R1=2*AMGLSS*AMBT*SINB**2*(SNIZ*(-AL1*AL2+BE1*BE2)*
+ $ AMBT*MZIZ*XT6+(AL2*BE1-AL1*BE2)*XT7)
+ XL1R2=-SINB**2*(BE1*BE2*XT8-AL1*AL2*AMBT**2*XT4+AMBT*
+ $ MZIZ*SNIZ*XT9*(-AL1*BE2+BE1*AL2))
+ XLR2=2*(XL1R1+XL1R2)
+ XMSB12=XL1L2+XR1R2+XLR1+XLR2
+ ELSE
+ XMSB12=0.
+ END IF
+ WID=XTCON*(XMSB1+XMSB2+XMSB12)
+ IF (WID.GT.0.) THEN
+ CALL SSSAVE(ISGL,WID,ISZI(IZ),IDBT,-IDBT,0,0)
+ Z1(1)=((ZI*ZADIZ-FB*ZMIXSS(2,IZ)*ZI**THIZ)*COSB-
+ $ (ZI*ZBDIZ-FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINB)/2.
+ Z1(2)=((-ZI*ZADIZ-FB*ZMIXSS(2,IZ)*ZI**THIZ)*COSB-
+ $ (ZI*ZBDIZ+FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINB)/2.
+ Z2(1)=(COSB-SINB)/2.
+ Z2(2)=-(COSB+SINB)/2.
+ CALL SSME3(3,AMB1SS,Z1,Z2)
+ Z1(1)=(COSB-SINB)/2.
+ Z1(2)=-(COSB+SINB)/2.
+ Z2(1)=CONJG((ZI*ZADIZ-FB*ZMIXSS(2,IZ)*ZI**THIZ)*COSB-
+ $ (ZI*ZBDIZ-FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINB)/2.
+ Z2(2)=-CONJG((-ZI*ZADIZ-FB*ZMIXSS(2,IZ)*ZI**THIZ)*COSB-
+ $ (ZI*ZBDIZ+FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINB)/2.
+ CALL SSME3(2,AMB1SS,Z1,Z2)
+ Z1(1)=((ZI*ZADIZ-FB*ZMIXSS(2,IZ)*ZI**THIZ)*SINB+
+ $ (ZI*ZBDIZ-FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSB)/2.
+ Z1(2)=((-ZI*ZADIZ-FB*ZMIXSS(2,IZ)*ZI**THIZ)*SINB+
+ $ (ZI*ZBDIZ+FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSB)/2.
+ Z2(1)=(COSB+SINB)/2.
+ Z2(2)=(COSB-SINB)/2.
+ CALL SSME3(3,AMB2SS,Z1,Z2)
+ Z1(1)=(COSB+SINB)/2.
+ Z1(2)=(COSB-SINB)/2.
+ Z2(1)=CONJG((ZI*ZADIZ-FB*ZMIXSS(2,IZ)*ZI**THIZ)*SINB+
+ $ (ZI*ZBDIZ-FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSB)/2.
+ Z2(2)=-CONJG((-ZI*ZADIZ-FB*ZMIXSS(2,IZ)*ZI**THIZ)*SINB+
+ $ (ZI*ZBDIZ+FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSB)/2.
+ CALL SSME3(2,AMB2SS,Z1,Z2)
+ END IF
+C 3 body gluino --> q + qb + zi decay, q=t
+ IF (AMGLSS.GT.(MZIZ+2*AMTP).AND.AMGLSS.LT.
+ $ (AMT1SS+AMTP)) THEN
+ TMP(1)=AMGLSS
+ TMP(2)=AMTP
+ TMP(3)=MZIZ
+ TMP(4)=AMT1SS
+ TMP(5)=AMT1SS
+ XUPPER=(AMGLSS**2+AMTP**2-(AMTP+MZIZ)**2)/2./AMGLSS
+ XT1=SSXINT(AMTP,SSGX1,XUPPER)
+ XT2=SSXINT(AMTP,SSGX2,XUPPER)
+ XT3=SSXINT(AMTP,SSGX3,XUPPER)
+ XT4=SSXINT(AMTP,SSGX4,XUPPER)
+ XT5=SSXINT(AMTP,SSGX5,XUPPER)
+ XT6=SSXINT(AMTP,SSGX6,XUPPER)
+ XT7=SSXINT(AMTP,SSGX7,XUPPER)
+ XT8=SSXINT(AMTP,SSGX8,XUPPER)
+ XT9=SSXINT(AMTP,SSGX9,XUPPER)
+ AL1=AUIZ*COST-FT*ZMIXSS(1,IZ)*SINT
+ BE1=FT*ZMIXSS(1,IZ)*COST+BUIZ*SINT
+C ---- here, al2 is (-) al2 of tata notes-----
+ AL2=BUIZ*SINT+FT*ZMIXSS(1,IZ)*COST
+ BE2=-FT*ZMIXSS(1,IZ)*SINT+AUIZ*COST
+ XLL=(AL1**2+BE1**2)*XT1-4*AMTP*MZIZ*SNIZ*AL1*
+ $ BE1*XT3-AMGLSS*(SNIZ*MZIZ*(AL1**2*XT2/AMGLSS/
+ $ MZIZ+BE1**2*AMTP**2*XT5)-AL1*BE1*AMTP*(XT4-
+ $ MZIZ**2*XT5))
+ XRR=(AL2**2+BE2**2)*XT1-4*AMTP*MZIZ*SNIZ*AL2*
+ $ BE2*XT3-AMGLSS*(SNIZ*MZIZ*(AL2**2*XT2/AMGLSS/
+ $ MZIZ+BE2**2*AMTP**2*XT5)-AL2*BE2*AMTP*(XT4-
+ $ MZIZ**2*XT5))
+ XL1R1=-2*AMGLSS*AMTP*((AL1*AL2+BE1*BE2)*SNIZ*AMTP*
+ $ MZIZ*XT6-(AL2*BE1+AL1*BE2)*XT7)
+ XL2R2=XL1R1
+ XL1R2=BE1*BE2*XT8+AL1*AL2*AMTP**2*XT4-AMTP*MZIZ*
+ $ SNIZ*(AL1*BE2+AL2*BE1)*XT9
+ XL2R1=XL1R2
+ XMST1=COST**2*XLL+SINT**2*XRR-SINT*COST*(XL1R1+XL1R2+
+ $ XL2R1+XL2R2)
+ ELSE
+ XMST1=0.
+ END IF
+ IF (AMGLSS.GT.(MZIZ+2*AMTP).AND.AMGLSS.LT.
+ $ (AMT2SS+AMTP)) THEN
+ TMP(1)=AMGLSS
+ TMP(2)=AMTP
+ TMP(3)=MZIZ
+ TMP(4)=AMT2SS
+ TMP(5)=AMT2SS
+ XUPPER=(AMGLSS**2+AMTP**2-(AMTP+MZIZ)**2)/2./AMGLSS
+ XT1=SSXINT(AMTP,SSGX1,XUPPER)
+ XT2=SSXINT(AMTP,SSGX2,XUPPER)
+ XT3=SSXINT(AMTP,SSGX3,XUPPER)
+ XT4=SSXINT(AMTP,SSGX4,XUPPER)
+ XT5=SSXINT(AMTP,SSGX5,XUPPER)
+ XT6=SSXINT(AMTP,SSGX6,XUPPER)
+ XT7=SSXINT(AMTP,SSGX7,XUPPER)
+ XT8=SSXINT(AMTP,SSGX8,XUPPER)
+ XT9=SSXINT(AMTP,SSGX9,XUPPER)
+ AL1=AUIZ*SINT+FT*ZMIXSS(1,IZ)*COST
+ BE1=FT*ZMIXSS(1,IZ)*SINT-BUIZ*COST
+C ---- here, al2 is (-) al2 of tata notes-----
+ AL2=-BUIZ*COST+FT*ZMIXSS(1,IZ)*SINT
+ BE2=FT*ZMIXSS(1,IZ)*COST+AUIZ*SINT
+ XLL=(AL1**2+BE1**2)*XT1-4*AMTP*MZIZ*SNIZ*AL1*
+ $ BE1*XT3-AMGLSS*(SNIZ*MZIZ*(AL1**2*XT2/AMGLSS/
+ $ MZIZ+BE1**2*AMTP**2*XT5)-AL1*BE1*AMTP*(XT4-
+ $ MZIZ**2*XT5))
+ XRR=(AL2**2+BE2**2)*XT1-4*AMTP*MZIZ*SNIZ*AL2*
+ $ BE2*XT3-AMGLSS*(SNIZ*MZIZ*(AL2**2*XT2/AMGLSS/
+ $ MZIZ+BE2**2*AMTP**2*XT5)-AL2*BE2*AMTP*(XT4-
+ $ MZIZ**2*XT5))
+ XL1R1=-2*AMGLSS*AMTP*((AL1*AL2+BE1*BE2)*SNIZ*AMTP*
+ $ MZIZ*XT6-(AL2*BE1+AL1*BE2)*XT7)
+ XL2R2=XL1R1
+ XL1R2=BE1*BE2*XT8+AL1*AL2*AMTP**2*XT4-AMTP*MZIZ*
+ $ SNIZ*(AL1*BE2+AL2*BE1)*XT9
+ XL2R1=XL1R2
+ XMST2=SINT**2*XLL+COST**2*XRR+SINT*COST*(XL1R1+XL1R2+
+ $ XL2R1+XL2R2)
+ ELSE
+ XMST2=0.
+ END IF
+C ----cross term between t_1 and t_2 graphs -----------
+ IF (AMGLSS.GT.(MZIZ+2*AMTP).AND.AMGLSS.LT.
+ $ (AMT1SS+AMTP)) THEN
+ TMP(1)=AMGLSS
+ TMP(2)=AMTP
+ TMP(3)=MZIZ
+ TMP(4)=AMT1SS
+ TMP(5)=AMT2SS
+ XUPPER=(AMGLSS**2+AMTP**2-(AMTP+MZIZ)**2)/2./AMGLSS
+ XT1=SSXINT(AMTP,SSGX1,XUPPER)
+ XT2=SSXINT(AMTP,SSGX2,XUPPER)
+ XT3=SSXINT(AMTP,SSGX3,XUPPER)
+ XT4=SSXINT(AMTP,SSGX4,XUPPER)
+ XT5=SSXINT(AMTP,SSGX5,XUPPER)
+ XT6=SSXINT(AMTP,SSGX6,XUPPER)
+ XT7=SSXINT(AMTP,SSGX7,XUPPER)
+ XT8=SSXINT(AMTP,SSGX8,XUPPER)
+ XT9=SSXINT(AMTP,SSGX9,XUPPER)
+ AL1=AUIZ*COST-FT*ZMIXSS(1,IZ)*SINT
+ AL2=AUIZ*SINT+FT*ZMIXSS(1,IZ)*COST
+ BE1=FT*ZMIXSS(1,IZ)*COST+BUIZ*SINT
+ BE2=FT*ZMIXSS(1,IZ)*SINT-BUIZ*COST
+ XL1L2=COST*SINT*(2*(AL1*AL2+BE1*BE2)*XT1-4*SNIZ*AMTP*
+ $ MZIZ*(AL1*BE2+AL2*BE1)*XT3-AMGLSS*(2*MZIZ*
+ $ SNIZ*(AL1*AL2*XT2/AMGLSS/MZIZ+BE1*BE2*AMTP**2*
+ $ XT5)-(AL1*BE2+AL2*BE1)*AMTP*(XT4-MZIZ**2*XT5)))
+ AL1=-BUIZ*SINT-FT*ZMIXSS(1,IZ)*COST
+ AL2=BUIZ*COST-FT*ZMIXSS(1,IZ)*SINT
+ BE1=-FT*ZMIXSS(1,IZ)*SINT+AUIZ*COST
+ BE2=FT*ZMIXSS(1,IZ)*COST+AUIZ*SINT
+ XR1R2=-COST*SINT*(2*(AL1*AL2+BE1*BE2)*XT1+4*SNIZ*AMTP*
+ $ MZIZ*(AL1*BE2+AL2*BE1)*XT3-AMGLSS*(2*MZIZ*
+ $ SNIZ*(AL1*AL2*XT2/AMGLSS/MZIZ+BE1*BE2*AMTP**2*
+ $ XT5)+(AL1*BE2+AL2*BE1)*AMTP*(XT4-MZIZ**2*XT5)))
+ AL1=AUIZ*COST-FT*ZMIXSS(1,IZ)*SINT
+ AL2=BUIZ*COST-FT*ZMIXSS(1,IZ)*SINT
+ BE1=FT*ZMIXSS(1,IZ)*COST+BUIZ*SINT
+ BE2=FT*ZMIXSS(1,IZ)*COST+AUIZ*SINT
+ XL1R1=2*AMGLSS*AMTP*COST**2*(SNIZ*(AL1*AL2-BE1*BE2)*
+ $ AMTP*MZIZ*XT6-(AL2*BE1-AL1*BE2)*XT7)
+ XL1R2=COST**2*(BE1*BE2*XT8-AL1*AL2*AMTP**2*XT4+AMTP*
+ $ MZIZ*SNIZ*XT9*(-AL1*BE2+BE1*AL2))
+ XLR1=2*(XL1R1+XL1R2)
+ AL1=AUIZ*SINT+FT*ZMIXSS(1,IZ)*COST
+ AL2=-BUIZ*SINT-FT*ZMIXSS(1,IZ)*COST
+ BE1=FT*ZMIXSS(1,IZ)*SINT-BUIZ*COST
+ BE2=-FT*ZMIXSS(1,IZ)*SINT+AUIZ*COST
+ TMP(4)=AMT2SS
+ TMP(5)=AMT1SS
+ XT8=SSXINT(AMTP,SSGX8,XUPPER)
+ XT9=SSXINT(AMTP,SSGX9,XUPPER)
+ XL1R1=2*AMGLSS*AMTP*SINT**2*(SNIZ*(-AL1*AL2+BE1*BE2)*
+ $ AMTP*MZIZ*XT6+(AL2*BE1-AL1*BE2)*XT7)
+ XL1R2=-SINT**2*(BE1*BE2*XT8-AL1*AL2*AMTP**2*XT4+AMTP*
+ $ MZIZ*SNIZ*XT9*(-AL1*BE2+BE1*AL2))
+ XLR2=2*(XL1R1+XL1R2)
+ XMST12=XL1L2+XR1R2+XLR1+XLR2
+ ELSE
+ XMST12=0.
+ END IF
+ WID=XTCON*(XMST1+XMST2+XMST12)
+ IF (WID.GT.0.) THEN
+ CALL SSSAVE(ISGL,WID,ISZI(IZ),IDTP,-IDTP,0,0)
+ Z1(1)=((ZI*ZAUIZ-FT*ZMIXSS(1,IZ)*ZI**THIZ)*COST-
+ $ (ZI*ZBUIZ-FT*ZMIXSS(1,IZ)*(-ZI)**THIZ)*SINT)/2.
+ Z1(2)=((-ZI*ZAUIZ-FT*ZMIXSS(1,IZ)*ZI**THIZ)*COST-
+ $ (ZI*ZBUIZ+FT*ZMIXSS(1,IZ)*(-ZI)**THIZ)*SINT)/2.
+ Z2(1)=(COST-SINT)/2.
+ Z2(2)=-(COST+SINT)/2.
+ CALL SSME3(3,AMT1SS,Z1,Z2)
+ Z1(1)=(COST-SINT)/2.
+ Z1(2)=-(COST+SINT)/2.
+ Z2(1)=CONJG((ZI*ZAUIZ-FT*ZMIXSS(1,IZ)*ZI**THIZ)*COST-
+ $ (ZI*ZBUIZ-FT*ZMIXSS(1,IZ)*(-ZI)**THIZ)*SINT)/2.
+ Z2(2)=-CONJG((-ZI*ZAUIZ-FT*ZMIXSS(1,IZ)*ZI**THIZ)*COST-
+ $ (ZI*ZBUIZ+FT*ZMIXSS(1,IZ)*(-ZI)**THIZ)*SINT)/2.
+ CALL SSME3(2,AMT1SS,Z1,Z2)
+ Z1(1)=((ZI*ZAUIZ-FT*ZMIXSS(1,IZ)*ZI**THIZ)*SINT+
+ $ (ZI*ZBUIZ-FT*ZMIXSS(1,IZ)*(-ZI)**THIZ)*COST)/2.
+ Z1(2)=((-ZI*ZAUIZ-FT*ZMIXSS(1,IZ)*ZI**THIZ)*SINT+
+ $ (ZI*ZBUIZ+FT*ZMIXSS(1,IZ)*(-ZI)**THIZ)*COST)/2.
+ Z2(1)=(COST+SINT)/2.
+ Z2(2)=(COST-SINT)/2.
+ CALL SSME3(3,AMT2SS,Z1,Z2)
+ Z1(1)=(COST+SINT)/2.
+ Z1(2)=(COST-SINT)/2.
+ Z2(1)=CONJG((ZI*ZAUIZ-FT*ZMIXSS(1,IZ)*ZI**THIZ)*SINT+
+ $ (ZI*ZBUIZ-FT*ZMIXSS(1,IZ)*(-ZI)**THIZ)*COST)/2.
+ Z2(2)=-CONJG((-ZI*ZAUIZ-FT*ZMIXSS(1,IZ)*ZI**THIZ)*SINT+
+ $ (ZI*ZBUIZ+FT*ZMIXSS(1,IZ)*(-ZI)**THIZ)*COST)/2.
+ CALL SSME3(2,AMT2SS,Z1,Z2)
+ END IF
+100 CONTINUE
+C
+C gluino --> quark + squark mode
+C
+ IF (AMGLSS.GT.(AMULSS+AMUP)) THEN
+ GMQK=ALF3*AMGLSS*(1.+AMUP**2/AMGLSS**2-AMULSS**2/AMGLSS**2)*
+ $ SQRT(SSXLAM(1.,AMUP**2/AMGLSS**2,AMULSS**2/AMGLSS**2))/8.
+ CALL SSSAVE(ISGL,GMQK,-ISUPL,+IDUP,0,0,0)
+ CALL SSSAVE(ISGL,GMQK,+ISUPL,-IDUP,0,0,0)
+ END IF
+ IF (AMGLSS.GT.(AMDLSS+AMDN)) THEN
+ GMQK=ALF3*AMGLSS*(1.+AMDN**2/AMGLSS**2-AMDLSS**2/AMGLSS**2)*
+ $ SQRT(SSXLAM(1.,AMDN**2/AMGLSS**2,AMDLSS**2/AMGLSS**2))/8.
+ CALL SSSAVE(ISGL,GMQK,-ISDNL,+IDDN,0,0,0)
+ CALL SSSAVE(ISGL,GMQK,+ISDNL,-IDDN,0,0,0)
+ END IF
+ IF (AMGLSS.GT.(AMURSS+AMUP)) THEN
+ GMQK=ALF3*AMGLSS*(1.+AMUP**2/AMGLSS**2-AMURSS**2/AMGLSS**2)*
+ $ SQRT(SSXLAM(1.,AMUP**2/AMGLSS**2,AMURSS**2/AMGLSS**2))/8.
+ CALL SSSAVE(ISGL,GMQK,-ISUPR,+IDUP,0,0,0)
+ CALL SSSAVE(ISGL,GMQK,+ISUPR,-IDUP,0,0,0)
+ END IF
+ IF (AMGLSS.GT.(AMDRSS+AMDN)) THEN
+ GMQK=ALF3*AMGLSS*(1.+AMDN**2/AMGLSS**2-AMDRSS**2/AMGLSS**2)*
+ $ SQRT(SSXLAM(1.,AMDN**2/AMGLSS**2,AMDRSS**2/AMGLSS**2))/8.
+ CALL SSSAVE(ISGL,GMQK,-ISDNR,+IDDN,0,0,0)
+ CALL SSSAVE(ISGL,GMQK,+ISDNR,-IDDN,0,0,0)
+ END IF
+C
+ IF (AMGLSS.GT.(AMSLSS+AMST)) THEN
+ GMQK=ALF3*AMGLSS*(1.+AMST**2/AMGLSS**2-AMSLSS**2/AMGLSS**2)*
+ $ SQRT(SSXLAM(1.,AMST**2/AMGLSS**2,AMSLSS**2/AMGLSS**2))/8.
+ CALL SSSAVE(ISGL,GMQK,-ISSTL,+IDST,0,0,0)
+ CALL SSSAVE(ISGL,GMQK,+ISSTL,-IDST,0,0,0)
+ END IF
+ IF (AMGLSS.GT.(AMSRSS+AMST)) THEN
+ GMQK=ALF3*AMGLSS*(1.+AMST**2/AMGLSS**2-AMSRSS**2/AMGLSS**2)*
+ $ SQRT(SSXLAM(1.,AMST**2/AMGLSS**2,AMSRSS**2/AMGLSS**2))/8.
+ CALL SSSAVE(ISGL,GMQK,-ISSTR,+IDST,0,0,0)
+ CALL SSSAVE(ISGL,GMQK,+ISSTR,-IDST,0,0,0)
+ END IF
+C
+ IF (AMGLSS.GT.(AMCLSS+AMCH)) THEN
+ GMQK=ALF3*AMGLSS*(1.+AMCH**2/AMGLSS**2-AMCLSS**2/AMGLSS**2)*
+ $ SQRT(SSXLAM(1.,AMCH**2/AMGLSS**2,AMCLSS**2/AMGLSS**2))/8.
+ CALL SSSAVE(ISGL,GMQK,-ISCHL,+IDCH,0,0,0)
+ CALL SSSAVE(ISGL,GMQK,+ISCHL,-IDCH,0,0,0)
+ END IF
+ IF (AMGLSS.GT.(AMCRSS+AMCH)) THEN
+ GMQK=ALF3*AMGLSS*(1.+AMCH**2/AMGLSS**2-AMCRSS**2/AMGLSS**2)*
+ $ SQRT(SSXLAM(1.,AMCH**2/AMGLSS**2,AMCRSS**2/AMGLSS**2))/8.
+ CALL SSSAVE(ISGL,GMQK,-ISCHR,+IDCH,0,0,0)
+ CALL SSSAVE(ISGL,GMQK,+ISCHR,-IDCH,0,0,0)
+ END IF
+C
+C !!! NEED MIXING ANGLE PIECE LIKE STOPS
+ IF (AMGLSS.GT.(AMB1SS+AMBT)) THEN
+ GMQK=ALF3*AMGLSS*(1.+AMBT**2/AMGLSS**2-AMB1SS**2/AMGLSS**2)*
+ $ SQRT(SSXLAM(1.,AMBT**2/AMGLSS**2,AMB1SS**2/AMGLSS**2))/8.
+ CALL SSSAVE(ISGL,GMQK,-ISBT1,+IDBT,0,0,0)
+ CALL SSSAVE(ISGL,GMQK,+ISBT1,-IDBT,0,0,0)
+ END IF
+C
+ IF (AMGLSS.GT.(AMB2SS+AMBT)) THEN
+ GMQK=ALF3*AMGLSS*(1.+AMBT**2/AMGLSS**2-AMB2SS**2/AMGLSS**2)*
+ $ SQRT(SSXLAM(1.,AMBT**2/AMGLSS**2,AMB2SS**2/AMGLSS**2))/8.
+ CALL SSSAVE(ISGL,GMQK,-ISBT2,+IDBT,0,0,0)
+ CALL SSSAVE(ISGL,GMQK,+ISBT2,-IDBT,0,0,0)
+ END IF
+C
+ IF (AMGLSS.GT.(AMT1SS+AMTP)) THEN
+ GMQK=ALF3*AMGLSS*(1.+AMTP**2/AMGLSS**2-AMT1SS**2/AMGLSS**2-
+ $ 2*SIN(2*THETAT)*AMTP/AMGLSS)*
+ $ SQRT(SSXLAM(1.,AMTP**2/AMGLSS**2,AMT1SS**2/AMGLSS**2))/8.
+ CALL SSSAVE(ISGL,GMQK,-ISTP1,+IDTP,0,0,0)
+ CALL SSSAVE(ISGL,GMQK,+ISTP1,-IDTP,0,0,0)
+ END IF
+C
+ IF (AMGLSS.GT.(AMT2SS+AMTP)) THEN
+ GMQK=ALF3*AMGLSS*(1.+AMTP**2/AMGLSS**2-AMT2SS**2/AMGLSS**2+
+ $ 2*SIN(2*THETAT)*AMTP/AMGLSS)*
+ $ SQRT(SSXLAM(1.,AMTP**2/AMGLSS**2,AMT2SS**2/AMGLSS**2))/8.
+ CALL SSSAVE(ISGL,GMQK,-ISTP2,+IDTP,0,0,0)
+ CALL SSSAVE(ISGL,GMQK,+ISTP2,-IDTP,0,0,0)
+ END IF
+C
+C Decay to gravitino
+C
+ IF (AMGLSS.GT.AMGVSS) THEN
+ WID=AMGLSS**5/48./PI/(AMGVSS*AMPL)**2
+ CALL SSSAVE(ISGL,WID,91,IDGL,0,0,0)
+ END IF
+C
+C Normalize branching ratios
+C
+ CALL SSNORM(ISGL)
+C
+ RETURN
+ END
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff