ISAJET/isasusy/ssglbf.F

   1 #include "isajet/pilot.h"
   2         SUBROUTINE SSGLBF
   3 C-----------------------------------------------------------------------
   4 C
   5 C       This subroutine gives gluino branching fractions to gauginos
   6 C       according to Baer,Barger,Karatas,Tata (Phys.Rev.D36,96(1987)
   7 C       (Now includes gluino->gluon+zino1,2,3,4 loop decays. Jan 1990)
   8 C       Also includes contribution due to non-degenerate t1-t2 stops
   9 C       Also includes contribution due to non-degenerate b_L and b_R
  10 C       Updated to include mixed sbottom states b1 and b2: 10/9/96
  11 C
  12 C       Auxiliary functions are called SSGxyi, where normally x
  13 C       indicates the SUSY particle, y the SM particle(s), and i is
  14 C       a counter.
  15 C
  16 C       Baer's GLUBF
  17 C
  18 C-----------------------------------------------------------------------
  19 #if defined(CERNLIB_IMPNONE)
  20       IMPLICIT NONE
  21 #endif
  22 #include "isajet/sslun.inc"
  23 #include "isajet/ssmode.inc"
  24 #include "isajet/sssm.inc"
  25 #include "isajet/sspar.inc"
  26 #include "isajet/sstmp.inc"
  27 #include "isajet/sstype.inc"
  28 C
  29       EXTERNAL SSGWQ1,SSGWQ2,SSGZG1,SSGZG2,SSGZG3,SSGZT
  30       EXTERNAL SSGWT1,SSGWT2,SSGWT3,SSGWT4,SSGWT5,SSGWT6,
  31      $SSGWT7,SSGWT8
  32       EXTERNAL SSGX1,SSGX2,SSGX3,SSGX4,SSGX5,SSGX6,SSGX7,SSGX8
  33      $,SSGX9
  34       REAL WIDU,WIDD,WIDI,WIDL,WIDR
  35       REAL SSXINT,SSXLAM,XUPPER,XTCON,SUM,GMQK,C,G,FT,GP
  36      $,UPPER,GS,FB,TANB,PI,SR2,UPPR,ALF3
  37       REAL MW1,MW2,SNW1,SNW2,XM,YM,THX,THY,COST,SINT,COSB,SINB
  38       REAL MZIZ,AUIZ,ADIZ,BUIZ,BDIZ
  39       REAL WID,SNIZ,FACT
  40       REAL XT1,XT2,XT3,XT4,XT5,XT6,XT7,XT8,XT9,AL1,AL2,BE1,BE2
  41       REAL XLL,XRR,XL1R1,XL2R2,XL1R2,XL2R1,XL1L2,XR1R2,XLR1
  42       REAL XMST1,XMST2,XMST12,XLR2,XMSB1,XMSB2,XMSB12
  43       REAL BW1,BW2,GT1,GT2,GT1T2,GB1,GB2,GT1B1,GT1B2,GT2B1,GT2B2
  44       REAL KUL,KUR,KDL,KDR,KCL,KCR,KSL,KSR,KBL,KBR,KT1,KT2
  45       REAL XKUL,XKUR,XKDL,XKDR,XKSL,XKSR,XKCL,XKCR,XKBL,XKBR,
  46      $XKT1,XKT2
  47       REAL XI1UL,XI1UR,XI1DL,XI1DR,XI1SL,XI1SR,XI1CL,XI1CR,
  48      $XI1BL,XI1BR,XI1T1,XI1T2,XIT1,XIT2
  49       REAL ALT1,ART1,ALT2,ART2,ALB1,ARB1,ALB2,ARB2
  50       REAL BETA,ASMB,MBMB,MBQ,ASMT,MTMT,MTQ,SUALFS,AMPL
  51       INTEGER IZ,ISZI(4),THIZ
  52       COMPLEX ZONE,ZI,ZAT1(2),ZAT2(2),ZADW1,ZADW2,ZAUW1,ZAUW2
  53       COMPLEX ZAUIZ,ZADIZ,ZBUIZ,ZBDIZ,Z1(2),Z2(2)
  54       DOUBLE PRECISION SSALFS,SSMQCD
  55       SAVE ZONE,ZI
  56       DATA ZONE,ZI/(1.,0.),(0.,1.)/
  57
  58 C
  59 C          Partly duplicated from SSMASS.
  60 C
  61       AMPL=2.4E18
  62       PI=4.*ATAN(1.)
  63       SR2=SQRT(2.)
  64       G=SQRT(4*PI*ALFAEM/SN2THW)
  65       GP=G*SQRT(SN2THW/(1.-SN2THW))
  66       ALF3=SSALFS(DBLE(AMGLSS**2))
  67       GS=SQRT(4.*PI*ALF3)
  68       TANB=1./RV2V1
  69       BETA=ATAN(TANB)
  70       ASMB=SUALFS(AMBT**2,.36,AMTP,3)
  71       MBMB=AMBT*(1.-4*ASMB/3./PI)
  72       MBQ=SSMQCD(DBLE(MBMB),DBLE(AMGLSS))
  73       ASMT=SUALFS(AMTP**2,.36,AMTP,3)
  74       MTMT=AMTP/(1.+4*ASMT/3./PI+(16.11-1.04*(5.-6.63/AMTP))*
  75      $(ASMT/PI)**2)
  76       MTQ=SSMQCD(DBLE(MTMT),DBLE(AMGLSS))
  77       FB=G*MBQ/SR2/AMW/COS(BETA)
  78       FT=G*MTQ/SR2/AMW/SIN(BETA)
  79       MW1=ABS(AMW1SS)
  80       MW2=ABS(AMW2SS)
  81       SNW1=SIGN(1.,AMW1SS)
  82       SNW2=SIGN(1.,AMW2SS)
  83       XM=1./TAN(GAMMAL)
  84       YM=1./TAN(GAMMAR)
  85       THX=SIGN(1.,XM)
  86       THY=SIGN(1.,YM)
  87       COST=COS(THETAT)
  88       SINT=SIN(THETAT)
  89       COSB=COS(THETAB)
  90       SINB=SIN(THETAB)
  91 C
  92       ZADW1=ZI*G*SNW1*SIN(GAMMAR)
  93       ZAUW1=ZI*G*SIN(GAMMAL)
  94       ZADW2=ZI*G*SNW2*COS(GAMMAR)*THY
  95       ZAUW2=ZI*G*COS(GAMMAL)*THX
  96       BW1=-FT*SNW1*COS(GAMMAR)
  97       BW2=FT*SNW2*SIN(GAMMAR)*THY
  98       ZAT1(1)=ZADW1*COST+ZI*BW1*SINT
  99       ZAT2(1)=ZADW1*SINT-ZI*BW1*COST
 100       ZAT1(2)=ZADW2*COST+ZI*BW2*SINT
 101       ZAT2(2)=ZADW2*SINT-ZI*BW2*COST
 102       C=SSALFS(DBLE(AMGLSS**2))*AMGLSS/8./PI**2
 103 C
 104 C          gluino --> w1 + qk + qb
 105 C
 106 C     Now includes sbottom as well as stop mixing/Yukawa effects
 107 C     as of 3/31/97, thanks to M. Drees
 108       TMP(1)=MW1
 109       UPPR=(AMGLSS**2-MW1**2)/2./AMGLSS
 110       IF (AMGLSS.GT.(MW1+AMUP+AMDN)) THEN
 111         IF (AMGLSS.LT.AMULSS.AND.AMGLSS.LT.AMDLSS) THEN
 112           TMP(2)=AMULSS
 113           WIDU=ZADW1*CONJG(ZADW1)*SSXINT(0.,SSGWQ1,UPPR)
 114           TMP(2)=AMDLSS
 115           WIDD=ZAUW1*CONJG(ZAUW1)*SSXINT(0.,SSGWQ1,UPPR)
 116           WIDI=2*REAL(ZAUW1*ZADW1)*SSXINT(0.,SSGWQ2,UPPR)
 117           WID=GS**2/2./AMGLSS/(2*PI)**5*(WIDU+WIDD+WIDI)
 118         ELSE IF (AMGLSS.LT.AMULSS.AND.AMGLSS.GE.AMDLSS) THEN
 119           TMP(2)=AMULSS
 120           WIDU=ZADW1*CONJG(ZADW1)*SSXINT(0.,SSGWQ1,UPPR)
 121           WID=GS**2/2./AMGLSS/(2*PI)**5*WIDU
 122         ELSE IF (AMGLSS.LT.AMDLSS.AND.AMGLSS.GE.AMULSS) THEN
 123           TMP(2)=AMDLSS
 124           WIDD=ZAUW1*CONJG(ZAUW1)*SSXINT(0.,SSGWQ1,UPPR)
 125           WID=GS**2/2./AMGLSS/(2*PI)**5*WIDD
 126         ELSE
 127           WID=0.
 128         END IF
 129         CALL SSSAVE(ISGL,WID,+ISW1,+IDDN,-IDUP,0,0)
 130         CALL SSSAVE(ISGL,WID,-ISW1,+IDUP,-IDDN,0,0)
 131         Z1(1)=1.
 132         Z1(2)=-Z1(1)
 133         Z2(1)=G*SIN(GAMMAL)
 134         Z2(2)=Z2(1)
 135         CALL SSME3(2,AMDLSS,Z1,Z2)
 136         Z1(1)=G*SIN(GAMMAR)
 137         Z1(2)=-Z1(1)
 138         Z2(1)=1.
 139         Z2(2)=Z2(1)
 140         CALL SSME3(3,AMULSS,Z1,Z2)
 141       END IF
 142 C
 143       IF (AMGLSS.GT.(MW1+AMST+AMCH)) THEN
 144         IF (AMGLSS.LT.AMCLSS.AND.AMGLSS.LT.AMSLSS) THEN
 145           TMP(2)=AMCLSS
 146           WIDU=ZADW1*CONJG(ZADW1)*SSXINT(0.,SSGWQ1,UPPR)
 147           TMP(2)=AMSLSS
 148           WIDD=ZAUW1*CONJG(ZAUW1)*SSXINT(0.,SSGWQ1,UPPR)
 149           WIDI=2*REAL(ZAUW1*ZADW1)*SSXINT(0.,SSGWQ2,UPPR)
 150           WID=GS**2/2./AMGLSS/(2*PI)**5*(WIDU+WIDD+WIDI)
 151         ELSE IF (AMGLSS.LT.AMCLSS.AND.AMGLSS.GE.AMSLSS) THEN
 152           TMP(2)=AMCLSS
 153           WIDU=ZADW1*CONJG(ZADW1)*SSXINT(0.,SSGWQ1,UPPR)
 154           WID=GS**2/2./AMGLSS/(2*PI)**5*WIDU
 155         ELSE IF (AMGLSS.LT.AMSLSS.AND.AMGLSS.GE.AMCLSS) THEN
 156           TMP(2)=AMSLSS
 157           WIDD=ZAUW1*CONJG(ZAUW1)*SSXINT(0.,SSGWQ1,UPPR)
 158           WID=GS**2/2./AMGLSS/(2*PI)**5*WIDD
 159         ELSE
 160           WID=0.
 161         END IF
 162         CALL SSSAVE(ISGL,WID,+ISW1,+IDST,-IDCH,0,0)
 163         CALL SSSAVE(ISGL,WID,-ISW1,+IDCH,-IDST,0,0)
 164         Z1(1)=1.
 165         Z1(2)=-Z1(1)
 166         Z2(1)=G*SIN(GAMMAL)
 167         Z2(2)=Z2(1)
 168         CALL SSME3(2,AMSLSS,Z1,Z2)
 169         Z1(1)=G*SIN(GAMMAR)
 170         Z1(2)=-Z1(1)
 171         Z2(1)=1.
 172         Z2(2)=Z2(1)
 173         CALL SSME3(3,AMCLSS,Z1,Z2)
 174       END IF
 175 C
 176       IF (AMGLSS.GT.(MW1+AMTP+AMBT)) THEN
 177         ALT1=-G*SIN(GAMMAR)*COST+FT*COS(GAMMAR)*SINT
 178         ART1=-FB*COS(GAMMAL)*COST
 179         ALT2=-G*SIN(GAMMAR)*SINT-FT*COS(GAMMAR)*COST
 180         ART2=-FB*COS(GAMMAL)*SINT
 181         ALB1=-G*SIN(GAMMAL)*COSB+FB*COS(GAMMAL)*SINB
 182         ARB1=-FT*COS(GAMMAR)*COSB
 183         ALB2=-G*SIN(GAMMAL)*SINB-FB*COS(GAMMAL)*COSB
 184         ARB2=-FT*COS(GAMMAR)*SINB
 185         UPPER=(AMGLSS**2+AMTP**2-(MW1+AMBT)**2)/2./AMGLSS
 186         FACT=GS**2*PI**2/(2*PI)**5/2./AMGLSS
 187         TMP(1)=MW1
 188         TMP(2)=AMGLSS
 189         TMP(3)=AMTP
 190         IF (AMGLSS.LT.(AMTP+AMT1SS)) THEN
 191           TMP(6)=AMT1SS
 192           TMP(7)=AMT1SS
 193           GT1=(ALT1**2+ART1**2)*(SSXINT(AMTP,SSGWT1,UPPER)+
 194      ,        SIN(2*THETAT)*SSXINT(AMTP,SSGWT8,UPPER))
 195           TMP(7)=AMT2SS
 196           GT1T2=2*(ALT1*ALT2+ART1*ART2)*(SINT**2-COST**2)*
 197      ,          SSXINT(AMTP,SSGWT8,UPPER)
 198         ELSE
 199           GT1=0.
 200           GT1T2=0.
 201         END IF
 202         IF (AMGLSS.LT.(AMTP+AMT2SS)) THEN
 203           TMP(6)=AMT2SS
 204           TMP(7)=AMT2SS
 205           GT2=(ALT2**2+ART2**2)*(SSXINT(AMTP,SSGWT1,UPPER)-
 206      ,        SIN(2*THETAT)*SSXINT(AMTP,SSGWT8,UPPER))
 207         ELSE
 208           GT2=0.
 209         END IF
 210         IF (AMGLSS.LT.(AMB1SS+AMBT)) THEN
 211           TMP(4)=AMB1SS
 212           TMP(8)=SNW1
 213 C          Rewrite UPPR=(AMGLSS**2+AMBT**2-(AMTP+MW1)**2)/2./AMGLSS
 214           UPPR=((AMGLSS-AMTP-MW1)*(AMGLSS+AMTP+MW1)+AMBT**2)/(2*AMGLSS)
 215           GB1=(ALB1**2+ARB1**2)*SSXINT(AMBT,SSGWT2,UPPR)-
 216      ,        ALB1*ARB1*SSXINT(AMBT,SSGWT3,UPPR)
 217         ELSE
 218           GB1=0.
 219         END IF
 220         IF (AMGLSS.LT.(AMB2SS+AMBT)) THEN
 221           TMP(4)=AMB2SS
 222           TMP(8)=SNW1
 223 C          Rewrite UPPR=(AMGLSS**2+AMBT**2-(AMTP+MW1)**2)/2./AMGLSS
 224           UPPR=((AMGLSS-AMTP-MW1)*(AMGLSS+AMTP+MW1)+AMBT**2)/(2*AMGLSS)
 225           GB2=(ALB2**2+ARB2**2)*SSXINT(AMBT,SSGWT2,UPPR)-
 226      ,        ALB2*ARB2*SSXINT(AMBT,SSGWT3,UPPR)
 227         ELSE
 228           GB2=0.
 229         END IF
 230         IF (AMGLSS.LT.(AMTP+AMT1SS).AND.AMGLSS.LT.
 231      $    (AMB1SS+AMBT)) THEN
 232           TMP(4)=AMB1SS
 233           TMP(6)=AMT1SS
 234           TMP(8)=SNW1
 235           GT1B1=(COST*SINB*ART1*ALB1+SINT*COSB*ALT1*ARB1)*
 236      ,           SSXINT(AMTP,SSGWT6,UPPER)-
 237      ,          (COST*COSB*ALT1*ALB1+SINT*SINB*ART1*ARB1)*
 238      ,           SSXINT(AMTP,SSGWT4,UPPER)+
 239      ,          (COST*COSB*ALT1*ARB1+SINT*SINB*ART1*ALB1)*
 240      ,           SSXINT(AMTP,SSGWT5,UPPER)-
 241      ,          (COST*SINB*ART1*ARB1+SINT*COSB*ALT1*ALB1)*
 242      ,           SSXINT(AMTP,SSGWT7,UPPER)
 243         ELSE
 244           GT1B1=0.
 245         END IF
 246         IF (AMGLSS.LT.(AMTP+AMT1SS).AND.AMGLSS.LT.
 247      $    (AMB2SS+AMBT)) THEN
 248           TMP(4)=AMB2SS
 249           TMP(6)=AMT1SS
 250           TMP(8)=SNW1
 251           GT1B2=(-COST*COSB*ART1*ALB2+SINT*SINB*ALT1*ARB2)*
 252      ,           SSXINT(AMTP,SSGWT6,UPPER)-
 253      ,          (COST*COSB*ALT1*ALB2-SINT*COSB*ART1*ARB2)*
 254      ,           SSXINT(AMTP,SSGWT4,UPPER)+
 255      ,          (COST*SINB*ALT1*ARB2-SINT*COSB*ART1*ALB2)*
 256      ,           SSXINT(AMTP,SSGWT5,UPPER)-
 257      ,          (-COST*COSB*ART1*ARB2+SINT*SINB*ALT1*ALB2)*
 258      ,           SSXINT(AMTP,SSGWT7,UPPER)
 259         ELSE
 260           GT1B2=0.
 261         END IF
 262         IF (AMGLSS.LT.(AMTP+AMT2SS).AND.AMGLSS.LT.
 263      $    (AMB1SS+AMBT)) THEN
 264           TMP(4)=AMB1SS
 265           TMP(6)=AMT2SS
 266           TMP(8)=SNW1
 267           GT2B1=(SINT*SINB*ART2*ALB1-COST*COSB*ALT2*ARB1)*
 268      ,           SSXINT(AMTP,SSGWT6,UPPER)-
 269      ,          (SINT*COSB*ALT2*ALB1-COST*SINB*ART2*ARB1)*
 270      ,           SSXINT(AMTP,SSGWT4,UPPER)+
 271      ,          (SINT*COSB*ALT2*ARB1-COST*SINB*ART2*ALB1)*
 272      ,           SSXINT(AMTP,SSGWT5,UPPER)-
 273      ,          (SINT*SINB*ART2*ARB1-COST*COSB*ALT2*ALB1)*
 274      ,           SSXINT(AMTP,SSGWT7,UPPER)
 275         ELSE
 276           GT2B1=0.
 277         END IF
 278         IF (AMGLSS.LT.(AMTP+AMT2SS).AND.AMGLSS.LT.
 279      $    (AMB2SS+AMBT)) THEN
 280           TMP(4)=AMB2SS
 281           TMP(6)=AMT2SS
 282           TMP(8)=SNW1
 283           GT2B2=(-SINT*COSB*ART2*ALB2-COST*SINB*ALT2*ARB2)*
 284      ,           SSXINT(AMTP,SSGWT6,UPPER)-
 285      ,          (SINT*SINB*ALT2*ALB2+COST*COSB*ART2*ARB2)*
 286      ,           SSXINT(AMTP,SSGWT4,UPPER)+
 287      ,          (SINT*SINB*ALT2*ARB2+COST*COSB*ART2*ALB2)*
 288      ,           SSXINT(AMTP,SSGWT5,UPPER)-
 289      ,          (-SINT*COSB*ART2*ARB2-COST*SINB*ALT2*ALB2)*
 290      ,           SSXINT(AMTP,SSGWT7,UPPER)
 291         ELSE
 292           GT2B2=0.
 293         END IF
 294         WID=GT1+GT2+GT1T2+GB1+GB2+GT1B1+GT1B2+GT2B1+GT2B2
 295         WID=FACT*WID
 296         IF (WID.GT.0.) THEN
 297         CALL SSSAVE(ISGL,WID,+ISW1,+IDBT,-IDTP,0,0)
 298         Z1(1)=SINT-COST
 299         Z1(2)=SINT+COST
 300         Z2(1)=ALT1+SNW1*ART1
 301         Z2(2)=ALT1-SNW1*ART1
 302         CALL SSME3(2,AMT1SS,Z1,Z2)
 303         Z1(1)=-COST-SINT
 304         Z1(2)=-COST+SINT
 305         Z2(1)=ALT2+SNW1*ART2
 306         Z2(2)=ALT2-SNW1*ART2
 307         CALL SSME3(2,AMT2SS,Z1,Z2)
 308         Z1(1)=SNW1*ALB1+ARB1
 309         Z1(2)=-SNW1*ALB1+ARB1
 310         Z2(1)=-COSB+SINB
 311         Z2(2)=-COSB-SINB
 312         CALL SSME3(3,AMB1SS,Z1,Z2)
 313         Z1(1)=SNW1*ALB2+ARB2
 314         Z1(2)=-SNW1*ALB2+ARB2
 315         Z2(1)=-SINB-COSB
 316         Z2(2)=-SINB+COSB
 317         CALL SSME3(3,AMB2SS,Z1,Z2)
 318         CALL SSSAVE(ISGL,WID,-ISW1,+IDTP,-IDBT,0,0)
 319         Z1(1)=SINB-COSB
 320         Z1(2)=SINB+COSB
 321         Z2(1)=ALB1+SNW1*ARB1
 322         Z2(2)=ALB1-SNW1*ARB1
 323         CALL SSME3(2,AMB1SS,Z1,Z2)
 324         Z1(1)=-COSB-SINB
 325         Z1(2)=-COSB+SINB
 326         Z2(1)=ALB2+SNW1*ARB2
 327         Z2(2)=ALB2-SNW1*ARB2
 328         CALL SSME3(2,AMB2SS,Z1,Z2)
 329         Z1(1)=SNW1*ALT1+ART1
 330         Z1(2)=-SNW1*ALT1+ART1
 331         Z2(1)=-COST+SINT
 332         Z2(2)=-COST-SINT
 333         CALL SSME3(3,AMT1SS,Z1,Z2)
 334         Z1(1)=SNW1*ALT2+ART2
 335         Z1(2)=-SNW1*ALT2+ART2
 336         Z2(1)=-SINT-COST
 337         Z2(2)=-SINT+COST
 338         CALL SSME3(3,AMT2SS,Z1,Z2)
 339         END IF
 340       END IF
 341 C
 342 C          gluino --> w2 + qk + qb
 343 C
 344       TMP(1)=MW2
 345       UPPR=(AMGLSS**2-MW2**2)/2./AMGLSS
 346       IF (AMGLSS.GT.(MW2+AMUP+AMDN)) THEN
 347         IF (AMGLSS.LT.AMULSS.AND.AMGLSS.LT.AMDLSS) THEN
 348           TMP(2)=AMULSS
 349           WIDU=ZADW2*CONJG(ZADW2)*SSXINT(0.,SSGWQ1,UPPR)
 350           TMP(2)=AMDLSS
 351           WIDD=ZAUW2*CONJG(ZAUW2)*SSXINT(0.,SSGWQ1,UPPR)
 352           WIDI=2*REAL(ZAUW2*ZADW2)*SSXINT(0.,SSGWQ2,UPPR)
 353           WID=GS**2/2./AMGLSS/(2*PI)**5*(WIDU+WIDD+WIDI)
 354         ELSE IF (AMGLSS.LT.AMULSS.AND.AMGLSS.GE.AMDLSS) THEN
 355           TMP(2)=AMULSS
 356           WIDU=ZADW2*CONJG(ZADW2)*SSXINT(0.,SSGWQ1,UPPR)
 357           WID=GS**2/2./AMGLSS/(2*PI)**5*WIDU
 358         ELSE IF (AMGLSS.LT.AMDLSS.AND.AMGLSS.GE.AMULSS) THEN
 359           TMP(2)=AMDLSS
 360           WIDD=ZAUW2*CONJG(ZAUW2)*SSXINT(0.,SSGWQ1,UPPR)
 361           WID=GS**2/2./AMGLSS/(2*PI)**5*WIDD
 362         ELSE
 363           WID=0.
 364         END IF
 365         CALL SSSAVE(ISGL,WID,+ISW2,+IDDN,-IDUP,0,0)
 366         CALL SSSAVE(ISGL,WID,-ISW2,+IDUP,-IDDN,0,0)
 367         Z1(1)=1.
 368         Z1(2)=-Z1(1)
 369         Z2(1)=G*THX*COS(GAMMAL)
 370         Z2(2)=Z2(1)
 371         CALL SSME3(2,AMDLSS,Z1,Z2)
 372         Z1(1)=G*THY*COS(GAMMAR)
 373         Z1(2)=-Z1(1)
 374         Z2(1)=1.
 375         Z2(2)=Z2(1)
 376         CALL SSME3(3,AMULSS,Z1,Z2)
 377       END IF
 378 C
 379       IF (AMGLSS.GT.(MW2+AMST+AMCH)) THEN
 380         IF (AMGLSS.LT.AMCLSS.AND.AMGLSS.LT.AMSLSS) THEN
 381           TMP(2)=AMCLSS
 382           WIDU=ZADW2*CONJG(ZADW2)*SSXINT(0.,SSGWQ1,UPPR)
 383           TMP(2)=AMSLSS
 384           WIDD=ZAUW2*CONJG(ZAUW2)*SSXINT(0.,SSGWQ1,UPPR)
 385           WIDI=2*REAL(ZAUW2*ZADW2)*SSXINT(0.,SSGWQ2,UPPR)
 386           WID=GS**2/2./AMGLSS/(2*PI)**5*(WIDU+WIDD+WIDI)
 387         ELSE IF (AMGLSS.LT.AMCLSS.AND.AMGLSS.GE.AMSLSS) THEN
 388           TMP(2)=AMCLSS
 389           WIDU=ZADW2*CONJG(ZADW2)*SSXINT(0.,SSGWQ1,UPPR)
 390           WID=GS**2/2./AMGLSS/(2*PI)**5*WIDU
 391         ELSE IF (AMGLSS.LT.AMSLSS.AND.AMGLSS.GE.AMCLSS) THEN
 392           TMP(2)=AMSLSS
 393           WIDD=ZAUW2*CONJG(ZAUW2)*SSXINT(0.,SSGWQ1,UPPR)
 394           WID=GS**2/2./AMGLSS/(2*PI)**5*WIDD
 395         ELSE
 396           WID=0.
 397         END IF
 398         CALL SSSAVE(ISGL,WID,+ISW2,+IDST,-IDCH,0,0)
 399         CALL SSSAVE(ISGL,WID,-ISW2,+IDCH,-IDST,0,0)
 400         Z1(1)=1.
 401         Z1(2)=-Z1(1)
 402         Z2(1)=G*THX*COS(GAMMAL)
 403         Z2(2)=Z2(1)
 404         CALL SSME3(2,AMSLSS,Z1,Z2)
 405         Z1(1)=G*THY*COS(GAMMAR)
 406         Z1(2)=-Z1(1)
 407         Z2(1)=1.
 408         Z2(2)=Z2(1)
 409         CALL SSME3(3,AMCLSS,Z1,Z2)
 410       END IF
 411 C
 412       IF (AMGLSS.GT.(MW2+AMTP+AMBT)) THEN
 413         ALT1=-G*THY*COS(GAMMAR)*COST-FT*THY*SIN(GAMMAR)*SINT
 414         ART1=FB*THX*SIN(GAMMAL)*COST
 415         ALT2=-G*THY*COS(GAMMAR)*SINT+FT*THY*SIN(GAMMAR)*COST
 416         ART2=FB*THX*SIN(GAMMAL)*SINT
 417         ALB1=-G*THX*COS(GAMMAL)*COSB-FB*THX*SIN(GAMMAL)*SINB
 418         ARB1=FT*THY*SIN(GAMMAR)*COSB
 419         ALB2=-G*THX*COS(GAMMAL)*SINB+FB*THX*SIN(GAMMAL)*COSB
 420         ARB2=FT*THY*SIN(GAMMAR)*SINB
 421         UPPER=(AMGLSS**2+AMTP**2-(MW2+AMBT)**2)/2./AMGLSS
 422         FACT=GS**2*PI**2/(2*PI)**5/2./AMGLSS
 423         TMP(1)=MW2
 424         TMP(2)=AMGLSS
 425         TMP(3)=AMTP
 426         IF (AMGLSS.LT.(AMTP+AMT1SS)) THEN
 427           TMP(6)=AMT1SS
 428           TMP(7)=AMT1SS
 429           GT1=(ALT1**2+ART1**2)*(SSXINT(AMTP,SSGWT1,UPPER)+
 430      ,        SIN(2*THETAT)*SSXINT(AMTP,SSGWT8,UPPER))
 431           TMP(7)=AMT2SS
 432           GT1T2=2*(ALT1*ALT2+ART1*ART2)*(SINT**2-COST**2)*
 433      ,          SSXINT(AMTP,SSGWT8,UPPER)
 434         ELSE
 435           GT1=0.
 436           GT1T2=0.
 437         END IF
 438         IF (AMGLSS.LT.(AMTP+AMT2SS)) THEN
 439           TMP(6)=AMT2SS
 440           TMP(7)=AMT2SS
 441           GT2=(ALT2**2+ART2**2)*(SSXINT(AMTP,SSGWT1,UPPER)-
 442      ,        SIN(2*THETAT)*SSXINT(AMTP,SSGWT8,UPPER))
 443         ELSE
 444           GT2=0.
 445         END IF
 446         IF (AMGLSS.LT.(AMB1SS+AMBT)) THEN
 447           TMP(4)=AMB1SS
 448           TMP(8)=SNW2
 449 C          Rewrite UPPR=(AMGLSS**2+AMBT**2-(AMTP+MW2)**2)/2./AMGLSS
 450           UPPR=((AMGLSS-AMTP-MW2)*(AMGLSS+AMTP+MW2)+AMBT**2)/(2*AMGLSS)
 451           GB1=(ALB1**2+ARB1**2)*SSXINT(AMBT,SSGWT2,UPPR)-
 452      ,        ALB1*ARB1*SSXINT(AMBT,SSGWT3,UPPR)
 453         ELSE
 454           GB1=0.
 455         END IF
 456         IF (AMGLSS.LT.(AMB2SS+AMBT)) THEN
 457           TMP(4)=AMB2SS
 458           TMP(8)=SNW2
 459           UPPR=((AMGLSS-AMTP-MW2)*(AMGLSS+AMTP+MW2)+AMBT**2)/(2*AMGLSS)
 460           GB2=(ALB2**2+ARB2**2)*SSXINT(AMBT,SSGWT2,UPPR)-
 461      ,        ALB2*ARB2*SSXINT(AMBT,SSGWT3,UPPR)
 462         ELSE
 463           GB2=0.
 464         END IF
 465         IF (AMGLSS.LT.(AMTP+AMT1SS).AND.AMGLSS.LT.
 466      $    (AMB1SS+AMBT)) THEN
 467           TMP(4)=AMB1SS
 468           TMP(6)=AMT1SS
 469           TMP(8)=SNW2
 470           GT1B1=(COST*SINB*ART1*ALB1+SINT*COSB*ALT1*ARB1)*
 471      ,           SSXINT(AMTP,SSGWT6,UPPER)-
 472      ,          (COST*COSB*ALT1*ALB1+SINT*SINB*ART1*ARB1)*
 473      ,           SSXINT(AMTP,SSGWT4,UPPER)+
 474      ,          (COST*COSB*ALT1*ARB1+SINT*SINB*ART1*ALB1)*
 475      ,           SSXINT(AMTP,SSGWT5,UPPER)-
 476      ,          (COST*SINB*ART1*ARB1+SINT*COSB*ALT1*ALB1)*
 477      ,           SSXINT(AMTP,SSGWT7,UPPER)
 478         ELSE
 479           GT1B1=0.
 480         END IF
 481         IF (AMGLSS.LT.(AMTP+AMT1SS).AND.AMGLSS.LT.
 482      $    (AMB2SS+AMBT)) THEN
 483           TMP(4)=AMB2SS
 484           TMP(6)=AMT1SS
 485           TMP(8)=SNW2
 486           GT1B2=(-COST*COSB*ART1*ALB2+SINT*SINB*ALT1*ARB2)*
 487      ,           SSXINT(AMTP,SSGWT6,UPPER)-
 488      ,          (COST*COSB*ALT1*ALB2-SINT*COSB*ART1*ARB2)*
 489      ,           SSXINT(AMTP,SSGWT4,UPPER)+
 490      ,          (COST*SINB*ALT1*ARB2-SINT*COSB*ART1*ALB2)*
 491      ,           SSXINT(AMTP,SSGWT5,UPPER)-
 492      ,          (-COST*COSB*ART1*ARB2+SINT*SINB*ALT1*ALB2)*
 493      ,           SSXINT(AMTP,SSGWT7,UPPER)
 494         ELSE
 495           GT1B2=0.
 496         END IF
 497         IF (AMGLSS.LT.(AMTP+AMT2SS).AND.AMGLSS.LT.
 498      $    (AMB1SS+AMBT)) THEN
 499           TMP(4)=AMB1SS
 500           TMP(6)=AMT2SS
 501           TMP(8)=SNW2
 502           GT2B1=(SINT*SINB*ART2*ALB1-COST*COSB*ALT2*ARB1)*
 503      ,           SSXINT(AMTP,SSGWT6,UPPER)-
 504      ,          (SINT*COSB*ALT2*ALB1-COST*SINB*ART2*ARB1)*
 505      ,           SSXINT(AMTP,SSGWT4,UPPER)+
 506      ,          (SINT*COSB*ALT2*ARB1-COST*SINB*ART2*ALB1)*
 507      ,           SSXINT(AMTP,SSGWT5,UPPER)-
 508      ,          (SINT*SINB*ART2*ARB1-COST*COSB*ALT2*ALB1)*
 509      ,           SSXINT(AMTP,SSGWT7,UPPER)
 510         ELSE
 511           GT2B1=0.
 512         END IF
 513         IF (AMGLSS.LT.(AMTP+AMT2SS).AND.AMGLSS.LT.
 514      $    (AMB2SS+AMBT)) THEN
 515           TMP(4)=AMB2SS
 516           TMP(6)=AMT2SS
 517           TMP(8)=SNW2
 518           GT2B2=(-SINT*COSB*ART2*ALB2-COST*SINB*ALT2*ARB2)*
 519      ,           SSXINT(AMTP,SSGWT6,UPPER)-
 520      ,          (SINT*SINB*ALT2*ALB2+COST*COSB*ART2*ARB2)*
 521      ,           SSXINT(AMTP,SSGWT4,UPPER)+
 522      ,          (SINT*SINB*ALT2*ARB2+COST*COSB*ART2*ALB2)*
 523      ,           SSXINT(AMTP,SSGWT5,UPPER)-
 524      ,          (-SINT*COSB*ART2*ARB2-COST*SINB*ALT2*ALB2)*
 525      ,           SSXINT(AMTP,SSGWT7,UPPER)
 526         ELSE
 527           GT2B2=0.
 528         END IF
 529         WID=GT1+GT2+GT1T2+GB1+GB2+GT1B1+GT1B2+GT2B1+GT2B2
 530         WID=FACT*WID
 531         IF (WID.GT.0.) THEN
 532         CALL SSSAVE(ISGL,WID,+ISW2,+IDBT,-IDTP,0,0)
 533         Z1(1)=SINT-COST
 534         Z1(2)=SINT+COST
 535         Z2(1)=ALT1+SNW2*ART1
 536         Z2(2)=ALT1-SNW2*ART1
 537         CALL SSME3(2,AMT1SS,Z1,Z2)
 538         Z1(1)=-COST-SINT
 539         Z1(2)=-COST+SINT
 540         Z2(1)=ALT2+SNW2*ART2
 541         Z2(2)=ALT2-SNW2*ART2
 542         CALL SSME3(2,AMT2SS,Z1,Z2)
 543         Z1(1)=SNW2*ALB1+ARB1
 544         Z1(2)=-SNW2*ALB1+ARB1
 545         Z2(1)=-COSB+SINB
 546         Z2(2)=-COSB-SINB
 547         CALL SSME3(3,AMB1SS,Z1,Z2)
 548         Z1(1)=SNW2*ALB2+ARB2
 549         Z1(2)=-SNW2*ALB2+ARB2
 550         Z2(1)=-SINB-COSB
 551         Z2(2)=-SINB+COSB
 552         CALL SSME3(3,AMB2SS,Z1,Z2)
 553         CALL SSSAVE(ISGL,WID,-ISW2,+IDTP,-IDBT,0,0)
 554         Z1(1)=SINB-COSB
 555         Z1(2)=SINB+COSB
 556         Z2(1)=ALB1+SNW2*ARB1
 557         Z2(2)=ALB1-SNW2*ARB1
 558         CALL SSME3(2,AMB1SS,Z1,Z2)
 559         Z1(1)=-COSB-SINB
 560         Z1(2)=-COSB+SINB
 561         Z2(1)=ALB2+SNW2*ARB2
 562         Z2(2)=ALB2-SNW2*ARB2
 563         CALL SSME3(2,AMB2SS,Z1,Z2)
 564         Z1(1)=SNW2*ALT1+ART1
 565         Z1(2)=-SNW2*ALT1+ART1
 566         Z2(1)=-COST+SINT
 567         Z2(2)=-COST-SINT
 568         CALL SSME3(3,AMT1SS,Z1,Z2)
 569         Z1(1)=SNW2*ALT2+ART2
 570         Z1(2)=-SNW2*ALT2+ART2
 571         Z2(1)=-SINT-COST
 572         Z2(2)=-SINT+COST
 573         CALL SSME3(3,AMT2SS,Z1,Z2)
 574         END IF
 575       END IF
 576 C
 577 C       gluino --> zi decays, zi = z1, z2, z3, z4
 578 C       the auiz etc, below are Atilde's etc. of PRD 42,1568 (1990)
 579 C
 580       ISZI(1)=ISZ1
 581       ISZI(2)=ISZ2
 582       ISZI(3)=ISZ3
 583       ISZI(4)=ISZ4
 584       DO 100 IZ=1,4
 585         MZIZ=ABS(AMZISS(IZ))
 586         AUIZ=G/SR2*ZMIXSS(3,IZ)+GP/3./SR2*ZMIXSS(4,IZ)
 587         ADIZ=-G/SR2*ZMIXSS(3,IZ)+GP/3./SR2*ZMIXSS(4,IZ)
 588         BUIZ=4*GP*ZMIXSS(4,IZ)/3./SR2
 589         BDIZ=-2*GP/3./SR2*ZMIXSS(4,IZ)
 590         SNIZ=SIGN(1.,AMZISS(IZ))
 591         THIZ=0
 592         IF (AMZISS(IZ).LT.0.) THIZ=1
 593         ZAUIZ=ZI**(THIZ-1)*SNIZ
 594      $  *(-G/SR2*ZMIXSS(3,IZ)-GP/3./SR2*ZMIXSS(4,IZ))
 595         ZBUIZ=ZI**(THIZ-1)*4*GP*ZMIXSS(4,IZ)/3./SR2
 596         ZADIZ=ZI**(THIZ-1)*SNIZ
 597      $  *(G/SR2*ZMIXSS(3,IZ)-GP/3./SR2*ZMIXSS(4,IZ))
 598         ZBDIZ=-2*ZI**(THIZ-1)*GP*ZMIXSS(4,IZ)/3./SR2
 599 C          Radiative gluino --> gluon + zi loop decay
 600         IF (AMGLSS.GT.MZIZ) THEN
 601           IF (AMGLSS.LT.(AMT1SS+AMTP)) THEN
 602             TMP(1)=AMTP
 603             TMP(2)=MZIZ
 604             TMP(3)=AMT1SS
 605             XIT1=SSXINT(0.,SSGZG1,1.)
 606             XI1T1=SSXINT(0.,SSGZG2,1.)
 607             XKT1=SSXINT(0.,SSGZG3,1.)
 608           ELSE
 609             XIT1=0.
 610             XI1T1=0.
 611             XKT1=0.
 612           END IF
 613           IF (AMGLSS.LT.(AMT2SS+AMTP)) THEN
 614             TMP(1)=AMTP
 615             TMP(2)=MZIZ
 616             TMP(3)=AMT2SS
 617             XIT2=SSXINT(0.,SSGZG1,1.)
 618             XI1T2=SSXINT(0.,SSGZG2,1.)
 619             XKT2=SSXINT(0.,SSGZG3,1.)
 620           ELSE
 621             XIT2=0.
 622             XI1T2=0.
 623             XKT2=0.
 624           END IF
 625 C         !!! NEEDS UPDATE FOR SBOTTOM MIXING !!!
 626           IF (AMGLSS.LT.(AMB1SS+AMBT)) THEN
 627             TMP(1)=AMBT
 628             TMP(2)=MZIZ
 629             TMP(3)=AMB1SS
 630             XI1BL=SSXINT(0.,SSGZG2,1.)
 631             XKBL=SSXINT(0.,SSGZG3,1.)
 632           ELSE
 633             XI1BL=0.
 634             XKBL=0.
 635           END IF
 636           IF (AMGLSS.LT.(AMB2SS+AMBT)) THEN
 637             TMP(1)=AMBT
 638             TMP(2)=MZIZ
 639             TMP(3)=AMB2SS
 640             XI1BR=SSXINT(0.,SSGZG2,1.)
 641             XKBR=SSXINT(0.,SSGZG3,1.)
 642           ELSE
 643             XI1BR=0.
 644             XKBR=0.
 645           END IF
 646           IF (AMGLSS.LT.(AMULSS+AMUP)) THEN
 647             TMP(1)=AMUP
 648             TMP(2)=MZIZ
 649             TMP(3)=AMULSS
 650             XI1UL=SSXINT(0.,SSGZG2,1.)
 651             XKUL=SSXINT(0.,SSGZG3,1.)
 652           ELSE
 653             XI1UL=0.
 654             XKUL=0.
 655           END IF
 656           IF (AMGLSS.LT.(AMURSS+AMUP)) THEN
 657             TMP(1)=AMUP
 658             TMP(2)=MZIZ
 659             TMP(3)=AMURSS
 660             XI1UR=SSXINT(0.,SSGZG2,1.)
 661             XKUR=SSXINT(0.,SSGZG3,1.)
 662           ELSE
 663             XI1UR=0.
 664             XKUR=0.
 665           END IF
 666           IF (AMGLSS.LT.(AMDLSS+AMDN)) THEN
 667             TMP(1)=AMDN
 668             TMP(2)=MZIZ
 669             TMP(3)=AMDLSS
 670             XI1DL=SSXINT(0.,SSGZG2,1.)
 671             XKDL=SSXINT(0.,SSGZG3,1.)
 672           ELSE
 673             XI1DL=0.
 674             XKDL=0.
 675           END IF
 676           IF (AMGLSS.LT.(AMDRSS+AMDN)) THEN
 677             TMP(1)=AMDN
 678             TMP(2)=MZIZ
 679             TMP(3)=AMDRSS
 680             XI1DR=SSXINT(0.,SSGZG2,1.)
 681             XKDR=SSXINT(0.,SSGZG3,1.)
 682           ELSE
 683             XI1DR=0.
 684             XKDR=0.
 685           END IF
 686           IF (AMGLSS.LT.(AMSLSS+AMST)) THEN
 687             TMP(1)=AMST
 688             TMP(2)=MZIZ
 689             TMP(3)=AMSLSS
 690             XI1SL=SSXINT(0.,SSGZG2,1.)
 691             XKSL=SSXINT(0.,SSGZG3,1.)
 692           ELSE
 693             XI1SL=0.
 694             XKSL=0.
 695           END IF
 696           IF (AMGLSS.LT.(AMSRSS+AMST)) THEN
 697             TMP(1)=AMST
 698             TMP(2)=MZIZ
 699             TMP(3)=AMSRSS
 700             XI1SR=SSXINT(0.,SSGZG2,1.)
 701             XKSR=SSXINT(0.,SSGZG3,1.)
 702           ELSE
 703             XI1SR=0.
 704             XKSR=0.
 705           END IF
 706           IF (AMGLSS.LT.(AMCLSS+AMCH)) THEN
 707             TMP(1)=AMCH
 708             TMP(2)=MZIZ
 709             TMP(3)=AMCLSS
 710             XI1CL=SSXINT(0.,SSGZG2,1.)
 711             XKCL=SSXINT(0.,SSGZG3,1.)
 712           ELSE
 713             XI1CL=0.
 714             XKCL=0.
 715           END IF
 716           IF (AMGLSS.LT.(AMCRSS+AMCH)) THEN
 717             TMP(1)=AMCH
 718             TMP(2)=MZIZ
 719             TMP(3)=AMCRSS
 720             XI1CR=SSXINT(0.,SSGZG2,1.)
 721             XKCR=SSXINT(0.,SSGZG3,1.)
 722           ELSE
 723             XI1CR=0.
 724             XKCR=0.
 725           END IF
 726           KUL=AUIZ*(XKUL*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1UL)
 727           KUR=-BUIZ*(XKUR*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1UR)
 728           KDL=ADIZ*(XKDL*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1DL)
 729           KDR=-BDIZ*(XKDR*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1DR)
 730           KCL=AUIZ*(XKCL*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1CL)
 731           KCR=-BUIZ*(XKCR*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1CR)
 732           KSL=ADIZ*(XKSL*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1SL)
 733           KSR=-BDIZ*(XKSR*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1SR)
 734           KBL=ADIZ*(XKBL*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1BL)
 735           KBR=-BDIZ*(XKBR*(MZIZ-SNIZ*AMGLSS)+MZIZ*XI1BR)
 736           KT1=(MZIZ*(XKT1+XI1T1)*(AUIZ*COST-FT*ZMIXSS(1,IZ)*SINT)
 737      $      -SNIZ*AMGLSS*XKT1*(AUIZ*COST-FT*ZMIXSS(1,IZ)*SINT)+
 738      $      SNIZ*AMTP*XIT1*(BUIZ*SINT+FT*ZMIXSS(1,IZ)*COST))*COST+
 739      $      (MZIZ*(XKT1+XI1T1)*(-BUIZ*SINT-FT*ZMIXSS(1,IZ)*COST)
 740      $      -SNIZ*AMGLSS*XKT1*(-BUIZ*SINT-FT*ZMIXSS(1,IZ)*COST)-
 741      $      SNIZ*AMTP*XIT1*(AUIZ*COST-FT*ZMIXSS(1,IZ)*SINT))*SINT
 742           KT2=(MZIZ*(XKT2+XI1T2)*(AUIZ*SINT+FT*ZMIXSS(1,IZ)*COST)
 743      $      -SNIZ*AMGLSS*XKT2*(AUIZ*SINT+FT*ZMIXSS(1,IZ)*COST)+
 744      $      SNIZ*AMTP*XIT2*(-BUIZ*COST+FT*ZMIXSS(1,IZ)*SINT))*SINT+
 745      $      (-MZIZ*(XKT2+XI1T2)*(BUIZ*COST-FT*ZMIXSS(1,IZ)*SINT)
 746      $      +SNIZ*AMGLSS*XKT2*(BUIZ*COST-FT*ZMIXSS(1,IZ)*SINT)+
 747      $      SNIZ*AMTP*XIT2*(AUIZ*SINT+FT*ZMIXSS(1,IZ)*COST))*COST
 748           SUM=(KUL+KUR+KDL+KDR+KSL+KSR+KCL+KCR+KBL+KBR+KT1+
 749      $         KT2)**2/AMGLSS**2
 750           WID=ALF3**2*AMGLSS*(1.-MZIZ**2/AMGLSS**2)/256./PI**3*SUM
 751           CALL SSSAVE(ISGL,WID,ISZI(IZ),IDGL,0,0,0)
 752         END IF
 753 C          3 body gluino --> q + qb + zi decay, q=u,d
 754         UPPR=(AMGLSS**2-MZIZ**2)/2./AMGLSS
 755         IF (AMGLSS.GT.(MZIZ+2*AMUP)) THEN
 756           IF (AMGLSS.LT.AMULSS.AND.AMGLSS.LT.AMURSS) THEN
 757             TMP(1)=MZIZ
 758             TMP(2)=AMULSS
 759             WIDL=2*AUIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
 760      $          SSXINT(0.,SSGWQ2,UPPR))
 761             TMP(2)=AMURSS
 762             WIDR=2*BUIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
 763      $          SSXINT(0.,SSGWQ2,UPPR))
 764             WID=WIDL+WIDR
 765           ELSE IF (AMGLSS.LT.AMULSS.AND.AMGLSS.GE.AMURSS) THEN
 766             TMP(2)=AMULSS
 767             WID=2*AUIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
 768      $          SSXINT(0.,SSGWQ2,UPPR))
 769           ELSE IF (AMGLSS.LT.AMURSS.AND.AMGLSS.GE.AMULSS) THEN
 770             TMP(2)=AMURSS
 771             WID=2*BUIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
 772      $          SSXINT(0.,SSGWQ2,UPPR))
 773           ELSE
 774             WID=0.
 775           END IF
 776           WID=GS**2/AMGLSS/2./(2*PI)**5*WID
 777           IF (WID.GT.0.) THEN
 778             CALL SSSAVE(ISGL,WID,ISZI(IZ),IDUP,-IDUP,0,0)
 779 C           Enter decay matrix element info
 780             Z1(1)=1.
 781             Z1(2)=-Z1(1)
 782             Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)*(THIZ+1)*AUIZ)
 783             Z2(2)=Z2(1)
 784             CALL SSME3(2,AMULSS,Z1,Z2)
 785             Z1(1)=1.
 786             Z1(2)=Z1(1)
 787             Z2(1)=-CONJG(ZI**(THIZ-1)*BUIZ)
 788             Z2(2)=-Z2(1)
 789             CALL SSME3(2,AMURSS,Z1,Z2)
 790             Z1(1)=ZI**(THIZ-1)*(-1.)*(THIZ+1)*AUIZ
 791             Z1(2)=-Z1(1)
 792             Z2(1)=1.
 793             Z2(2)=Z2(1)
 794             CALL SSME3(3,AMULSS,Z1,Z2)
 795             Z1(1)=ZI**(THIZ-1)*BUIZ
 796             Z1(2)=Z1(1)
 797             Z2(1)=1.
 798             Z2(2)=-Z2(1)
 799             CALL SSME3(3,AMURSS,Z1,Z2)
 800           END IF
 801         END IF
 802 C
 803         IF (AMGLSS.GT.(MZIZ+2*AMDN)) THEN
 804           IF (AMGLSS.LT.AMDLSS.AND.AMGLSS.LT.AMDRSS) THEN
 805             TMP(1)=MZIZ
 806             TMP(2)=AMDLSS
 807             WIDL=2*ADIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
 808      $          SSXINT(0.,SSGWQ2,UPPR))
 809             TMP(2)=AMDRSS
 810             WIDR=2*BDIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
 811      $          SSXINT(0.,SSGWQ2,UPPR))
 812             WID=WIDL+WIDR
 813           ELSE IF (AMGLSS.LT.AMDLSS.AND.AMGLSS.GE.AMDRSS) THEN
 814             TMP(2)=AMDLSS
 815             WID=2*ADIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
 816      $          SSXINT(0.,SSGWQ2,UPPR))
 817           ELSE IF (AMGLSS.LT.AMDRSS.AND.AMGLSS.GE.AMDLSS) THEN
 818             TMP(2)=AMDRSS
 819             WID=2*BDIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
 820      $          SSXINT(0.,SSGWQ2,UPPR))
 821           ELSE
 822             WID=0.
 823           END IF
 824           WID=GS**2/AMGLSS/2./(2*PI)**5*WID
 825           IF (WID.GT.0.) THEN
 826           CALL SSSAVE(ISGL,WID,ISZI(IZ),IDDN,-IDDN,0,0)
 827 C           Enter decay matrix element info
 828             Z1(1)=1.
 829             Z1(2)=-Z1(1)
 830             Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)*(THIZ+1)*ADIZ)
 831             Z2(2)=Z2(1)
 832             CALL SSME3(2,AMDLSS,Z1,Z2)
 833             Z1(1)=1.
 834             Z1(2)=Z1(1)
 835             Z2(1)=-CONJG(ZI**(THIZ-1)*BDIZ)
 836             Z2(2)=-Z2(1)
 837             CALL SSME3(2,AMDRSS,Z1,Z2)
 838             Z1(1)=ZI**(THIZ-1)*(-1.)*(THIZ+1)*ADIZ
 839             Z1(2)=-Z1(1)
 840             Z2(1)=1.
 841             Z2(2)=Z2(1)
 842             CALL SSME3(3,AMDLSS,Z1,Z2)
 843             Z1(1)=ZI**(THIZ-1)*BDIZ
 844             Z1(2)=Z1(1)
 845             Z2(1)=1.
 846             Z2(2)=-Z2(1)
 847             CALL SSME3(3,AMDRSS,Z1,Z2)
 848           END IF
 849         END IF
 850 C          3 body gluino --> q + qb + zi decay, q=s
 851         IF (AMGLSS.GT.(MZIZ+2*AMST)) THEN
 852           IF (AMGLSS.LT.AMSLSS.AND.AMGLSS.LT.AMSRSS) THEN
 853             TMP(1)=MZIZ
 854             TMP(2)=AMSLSS
 855             WIDL=2*ADIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
 856      $          SSXINT(0.,SSGWQ2,UPPR))
 857             TMP(2)=AMSRSS
 858             WIDR=2*BDIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
 859      $          SSXINT(0.,SSGWQ2,UPPR))
 860             WID=WIDL+WIDR
 861           ELSE IF (AMGLSS.LT.AMSLSS.AND.AMGLSS.GE.AMSRSS) THEN
 862             TMP(2)=AMSLSS
 863             WID=2*ADIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
 864      $          SSXINT(0.,SSGWQ2,UPPR))
 865           ELSE IF (AMGLSS.LT.AMSRSS.AND.AMGLSS.GE.AMSLSS) THEN
 866             TMP(2)=AMSRSS
 867             WID=2*BDIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
 868      $          SSXINT(0.,SSGWQ2,UPPR))
 869           ELSE
 870             WID=0.
 871           END IF
 872           WID=GS**2/AMGLSS/2./(2*PI)**5*WID
 873           IF (WID.GT.0.) THEN
 874           CALL SSSAVE(ISGL,WID,ISZI(IZ),IDST,-IDST,0,0)
 875 C           Enter decay matrix element info
 876             Z1(1)=1.
 877             Z1(2)=-Z1(1)
 878             Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)*(THIZ+1)*ADIZ)
 879             Z2(2)=Z2(1)
 880             CALL SSME3(2,AMDLSS,Z1,Z2)
 881             Z1(1)=1.
 882             Z1(2)=Z1(1)
 883             Z2(1)=-CONJG(ZI**(THIZ-1)*BDIZ)
 884             Z2(2)=-Z2(1)
 885             CALL SSME3(2,AMDRSS,Z1,Z2)
 886             Z1(1)=ZI**(THIZ-1)*(-1.)*(THIZ+1)*ADIZ
 887             Z1(2)=-Z1(1)
 888             Z2(1)=1.
 889             Z2(2)=Z2(1)
 890             CALL SSME3(3,AMDLSS,Z1,Z2)
 891             Z1(1)=ZI**(THIZ-1)*BDIZ
 892             Z1(2)=Z1(1)
 893             Z2(1)=1.
 894             Z2(2)=-Z2(1)
 895             CALL SSME3(3,AMDRSS,Z1,Z2)
 896           END IF
 897         END IF
 898 C          3 body gluino --> q + qb + zi decay, q=c
 899         IF (AMGLSS.GT.(MZIZ+2*AMCH)) THEN
 900           IF (AMGLSS.LT.AMCLSS.AND.AMGLSS.LT.AMCRSS) THEN
 901             TMP(1)=MZIZ
 902             TMP(2)=AMCLSS
 903             WIDL=2*AUIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
 904      $          SSXINT(0.,SSGWQ2,UPPR))
 905             TMP(2)=AMCRSS
 906             WIDR=2*BUIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
 907      $          SSXINT(0.,SSGWQ2,UPPR))
 908             WID=WIDL+WIDR
 909           ELSE IF (AMGLSS.LT.AMCLSS.AND.AMGLSS.GE.AMCRSS) THEN
 910             TMP(2)=AMCLSS
 911             WID=2*AUIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
 912      $          SSXINT(0.,SSGWQ2,UPPR))
 913           ELSE IF (AMGLSS.LT.AMCRSS.AND.AMGLSS.GE.AMCLSS) THEN
 914             TMP(2)=AMCRSS
 915             WID=2*BUIZ**2*(SSXINT(0.,SSGWQ1,UPPR)-SNIZ*
 916      $          SSXINT(0.,SSGWQ2,UPPR))
 917           ELSE
 918             WID=0.
 919           END IF
 920           WID=GS**2/AMGLSS/2./(2*PI)**5*WID
 921           IF (WID.GT.0.) THEN
 922           CALL SSSAVE(ISGL,WID,ISZI(IZ),IDCH,-IDCH,0,0)
 923 C           Enter decay matrix element info
 924             Z1(1)=1.
 925             Z1(2)=-Z1(1)
 926             Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)*(THIZ+1)*AUIZ)
 927             Z2(2)=Z2(1)
 928             CALL SSME3(2,AMULSS,Z1,Z2)
 929             Z1(1)=1.
 930             Z1(2)=Z1(1)
 931             Z2(1)=-CONJG(ZI**(THIZ-1)*BUIZ)
 932             Z2(2)=-Z2(1)
 933             CALL SSME3(2,AMURSS,Z1,Z2)
 934             Z1(1)=ZI**(THIZ-1)*(-1.)*(THIZ+1)*AUIZ
 935             Z1(2)=-Z1(1)
 936             Z2(1)=1.
 937             Z2(2)=Z2(1)
 938             CALL SSME3(3,AMULSS,Z1,Z2)
 939             Z1(1)=ZI**(THIZ-1)*BUIZ
 940             Z1(2)=Z1(1)
 941             Z2(1)=1.
 942             Z2(2)=-Z2(1)
 943             CALL SSME3(3,AMURSS,Z1,Z2)
 944           END IF
 945         END IF
 946 C          3 body gluino --> q + qb + zi decay, q=b
 947         XTCON=ALF3/8./PI**4/AMGLSS
 948         IF (AMGLSS.GT.(MZIZ+2*AMBT).AND.AMGLSS.LT.
 949      $     (AMB1SS+AMBT)) THEN
 950           TMP(1)=AMGLSS
 951           TMP(2)=AMBT
 952           TMP(3)=MZIZ
 953           TMP(4)=AMB1SS
 954           TMP(5)=AMB1SS
 955           XUPPER=(AMGLSS**2+AMBT**2-(AMBT+MZIZ)**2)/2./AMGLSS
 956           XT1=SSXINT(AMBT,SSGX1,XUPPER)
 957           XT2=SSXINT(AMBT,SSGX2,XUPPER)
 958           XT3=SSXINT(AMBT,SSGX3,XUPPER)
 959           XT4=SSXINT(AMBT,SSGX4,XUPPER)
 960           XT5=SSXINT(AMBT,SSGX5,XUPPER)
 961           XT6=SSXINT(AMBT,SSGX6,XUPPER)
 962           XT7=SSXINT(AMBT,SSGX7,XUPPER)
 963           XT8=SSXINT(AMBT,SSGX8,XUPPER)
 964           XT9=SSXINT(AMBT,SSGX9,XUPPER)
 965           AL1=ADIZ*COSB-FB*ZMIXSS(2,IZ)*SINB
 966           BE1=FB*ZMIXSS(2,IZ)*COSB+BDIZ*SINB
 967 C         ---- here, al2 is (-) al2 of tata notes-----
 968           AL2=BDIZ*SINB+FB*ZMIXSS(2,IZ)*COSB
 969           BE2=-FB*ZMIXSS(2,IZ)*SINB+ADIZ*COSB
 970           XLL=(AL1**2+BE1**2)*XT1-4*AMBT*MZIZ*SNIZ*AL1*
 971      $     BE1*XT3-AMGLSS*(SNIZ*MZIZ*(AL1**2*XT2/AMGLSS/
 972      $     MZIZ+BE1**2*AMBT**2*XT5)-AL1*BE1*AMBT*(XT4-
 973      $     MZIZ**2*XT5))
 974           XRR=(AL2**2+BE2**2)*XT1-4*AMBT*MZIZ*SNIZ*AL2*
 975      $     BE2*XT3-AMGLSS*(SNIZ*MZIZ*(AL2**2*XT2/AMGLSS/
 976      $     MZIZ+BE2**2*AMBT**2*XT5)-AL2*BE2*AMBT*(XT4-
 977      $     MZIZ**2*XT5))
 978           XL1R1=-2*AMGLSS*AMBT*((AL1*AL2+BE1*BE2)*SNIZ*AMBT*
 979      $     MZIZ*XT6-(AL2*BE1+AL1*BE2)*XT7)
 980           XL2R2=XL1R1
 981           XL1R2=BE1*BE2*XT8+AL1*AL2*AMBT**2*XT4-AMBT*MZIZ*
 982      $     SNIZ*(AL1*BE2+AL2*BE1)*XT9
 983           XL2R1=XL1R2
 984           XMSB1=COSB**2*XLL+SINB**2*XRR-SINB*COSB*(XL1R1+XL1R2+
 985      $     XL2R1+XL2R2)
 986         ELSE
 987           XMSB1=0.
 988         END IF
 989         IF (AMGLSS.GT.(MZIZ+2*AMBT).AND.AMGLSS.LT.
 990      $     (AMB2SS+AMBT)) THEN
 991           TMP(1)=AMGLSS
 992           TMP(2)=AMBT
 993           TMP(3)=MZIZ
 994           TMP(4)=AMB2SS
 995           TMP(5)=AMB2SS
 996           XUPPER=(AMGLSS**2+AMBT**2-(AMBT+MZIZ)**2)/2./AMGLSS
 997           XT1=SSXINT(AMBT,SSGX1,XUPPER)
 998           XT2=SSXINT(AMBT,SSGX2,XUPPER)
 999           XT3=SSXINT(AMBT,SSGX3,XUPPER)
1000           XT4=SSXINT(AMBT,SSGX4,XUPPER)
1001           XT5=SSXINT(AMBT,SSGX5,XUPPER)
1002           XT6=SSXINT(AMBT,SSGX6,XUPPER)
1003           XT7=SSXINT(AMBT,SSGX7,XUPPER)
1004           XT8=SSXINT(AMBT,SSGX8,XUPPER)
1005           XT9=SSXINT(AMBT,SSGX9,XUPPER)
1006           AL1=ADIZ*SINB+FB*ZMIXSS(2,IZ)*COSB
1007           BE1=FB*ZMIXSS(2,IZ)*SINB-BDIZ*COSB
1008 C         ---- here, al2 is (-) al2 of tata notes-----
1009           AL2=-BDIZ*COSB+FB*ZMIXSS(2,IZ)*SINB
1010           BE2=FB*ZMIXSS(2,IZ)*COSB+ADIZ*SINB
1011           XLL=(AL1**2+BE1**2)*XT1-4*AMBT*MZIZ*SNIZ*AL1*
1012      $     BE1*XT3-AMGLSS*(SNIZ*MZIZ*(AL1**2*XT2/AMGLSS/
1013      $     MZIZ+BE1**2*AMBT**2*XT5)-AL1*BE1*AMBT*(XT4-
1014      $     MZIZ**2*XT5))
1015           XRR=(AL2**2+BE2**2)*XT1-4*AMBT*MZIZ*SNIZ*AL2*
1016      $     BE2*XT3-AMGLSS*(SNIZ*MZIZ*(AL2**2*XT2/AMGLSS/
1017      $     MZIZ+BE2**2*AMBT**2*XT5)-AL2*BE2*AMBT*(XT4-
1018      $     MZIZ**2*XT5))
1019           XL1R1=-2*AMGLSS*AMBT*((AL1*AL2+BE1*BE2)*SNIZ*AMBT*
1020      $     MZIZ*XT6-(AL2*BE1+AL1*BE2)*XT7)
1021           XL2R2=XL1R1
1022           XL1R2=BE1*BE2*XT8+AL1*AL2*AMBT**2*XT4-AMBT*MZIZ*
1023      $     SNIZ*(AL1*BE2+AL2*BE1)*XT9
1024           XL2R1=XL1R2
1025           XMSB2=SINB**2*XLL+COSB**2*XRR+SINB*COSB*(XL1R1+XL1R2+
1026      $     XL2R1+XL2R2)
1027         ELSE
1028           XMSB2=0.
1029         END IF
1030 C       ----cross term between b_1 and b_2 graphs -----------
1031         IF (AMGLSS.GT.(MZIZ+2*AMBT).AND.AMGLSS.LT.
1032      $     (AMB1SS+AMBT)) THEN
1033           TMP(1)=AMGLSS
1034           TMP(2)=AMBT
1035           TMP(3)=MZIZ
1036           TMP(4)=AMB1SS
1037           TMP(5)=AMB2SS
1038           XUPPER=(AMGLSS**2+AMBT**2-(AMBT+MZIZ)**2)/2./AMGLSS
1039           XT1=SSXINT(AMBT,SSGX1,XUPPER)
1040           XT2=SSXINT(AMBT,SSGX2,XUPPER)
1041           XT3=SSXINT(AMBT,SSGX3,XUPPER)
1042           XT4=SSXINT(AMBT,SSGX4,XUPPER)
1043           XT5=SSXINT(AMBT,SSGX5,XUPPER)
1044           XT6=SSXINT(AMBT,SSGX6,XUPPER)
1045           XT7=SSXINT(AMBT,SSGX7,XUPPER)
1046           XT8=SSXINT(AMBT,SSGX8,XUPPER)
1047           XT9=SSXINT(AMBT,SSGX9,XUPPER)
1048           AL1=ADIZ*COSB-FB*ZMIXSS(2,IZ)*SINB
1049           AL2=ADIZ*SINB+FB*ZMIXSS(2,IZ)*COSB
1050           BE1=FB*ZMIXSS(2,IZ)*COSB+BDIZ*SINB
1051           BE2=FB*ZMIXSS(2,IZ)*SINB-BDIZ*COSB
1052           XL1L2=COSB*SINB*(2*(AL1*AL2+BE1*BE2)*XT1-4*SNIZ*AMBT*
1053      $     MZIZ*(AL1*BE2+AL2*BE1)*XT3-AMGLSS*(2*MZIZ*
1054      $     SNIZ*(AL1*AL2*XT2/AMGLSS/MZIZ+BE1*BE2*AMBT**2*
1055      $     XT5)-(AL1*BE2+AL2*BE1)*AMBT*(XT4-MZIZ**2*XT5)))
1056           AL1=-BDIZ*SINB-FB*ZMIXSS(2,IZ)*COSB
1057           AL2=BDIZ*COSB-FB*ZMIXSS(2,IZ)*SINB
1058           BE1=-FB*ZMIXSS(2,IZ)*SINB+ADIZ*COSB
1059           BE2=FB*ZMIXSS(2,IZ)*COSB+ADIZ*SINB
1060           XR1R2=-COSB*SINB*(2*(AL1*AL2+BE1*BE2)*XT1+4*SNIZ*AMBT*
1061      $     MZIZ*(AL1*BE2+AL2*BE1)*XT3-AMGLSS*(2*MZIZ*
1062      $     SNIZ*(AL1*AL2*XT2/AMGLSS/MZIZ+BE1*BE2*AMBT**2*
1063      $     XT5)+(AL1*BE2+AL2*BE1)*AMBT*(XT4-MZIZ**2*XT5)))
1064           AL1=ADIZ*COSB-FB*ZMIXSS(2,IZ)*SINB
1065           AL2=BDIZ*COSB-FB*ZMIXSS(2,IZ)*SINB
1066           BE1=FB*ZMIXSS(2,IZ)*COSB+BDIZ*SINB
1067           BE2=FB*ZMIXSS(2,IZ)*COSB+ADIZ*SINB
1068           XL1R1=2*AMGLSS*AMBT*COSB**2*(SNIZ*(AL1*AL2-BE1*BE2)*
1069      $     AMBT*MZIZ*XT6-(AL2*BE1-AL1*BE2)*XT7)
1070           XL1R2=COSB**2*(BE1*BE2*XT8-AL1*AL2*AMBT**2*XT4+AMBT*
1071      $     MZIZ*SNIZ*XT9*(-AL1*BE2+BE1*AL2))
1072           XLR1=2*(XL1R1+XL1R2)
1073           AL1=ADIZ*SINB+FB*ZMIXSS(2,IZ)*COSB
1074           AL2=-BDIZ*SINB-FB*ZMIXSS(2,IZ)*COSB
1075           BE1=FB*ZMIXSS(2,IZ)*SINB-BDIZ*COSB
1076           BE2=-FB*ZMIXSS(2,IZ)*SINB+ADIZ*COSB
1077           TMP(4)=AMB2SS
1078           TMP(5)=AMB1SS
1079           XT8=SSXINT(AMBT,SSGX8,XUPPER)
1080           XT9=SSXINT(AMBT,SSGX9,XUPPER)
1081           XL1R1=2*AMGLSS*AMBT*SINB**2*(SNIZ*(-AL1*AL2+BE1*BE2)*
1082      $     AMBT*MZIZ*XT6+(AL2*BE1-AL1*BE2)*XT7)
1083           XL1R2=-SINB**2*(BE1*BE2*XT8-AL1*AL2*AMBT**2*XT4+AMBT*
1084      $     MZIZ*SNIZ*XT9*(-AL1*BE2+BE1*AL2))
1085           XLR2=2*(XL1R1+XL1R2)
1086           XMSB12=XL1L2+XR1R2+XLR1+XLR2
1087         ELSE
1088           XMSB12=0.
1089         END IF
1090         WID=XTCON*(XMSB1+XMSB2+XMSB12)
1091         IF (WID.GT.0.) THEN
1092           CALL SSSAVE(ISGL,WID,ISZI(IZ),IDBT,-IDBT,0,0)
1093           Z1(1)=((ZI*ZADIZ-FB*ZMIXSS(2,IZ)*ZI**THIZ)*COSB-
1094      $           (ZI*ZBDIZ-FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINB)/2.
1095           Z1(2)=((-ZI*ZADIZ-FB*ZMIXSS(2,IZ)*ZI**THIZ)*COSB-
1096      $           (ZI*ZBDIZ+FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINB)/2.
1097           Z2(1)=(COSB-SINB)/2.
1098           Z2(2)=-(COSB+SINB)/2.
1099           CALL SSME3(3,AMB1SS,Z1,Z2)
1100           Z1(1)=(COSB-SINB)/2.
1101           Z1(2)=-(COSB+SINB)/2.
1102           Z2(1)=CONJG((ZI*ZADIZ-FB*ZMIXSS(2,IZ)*ZI**THIZ)*COSB-
1103      $           (ZI*ZBDIZ-FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINB)/2.
1104           Z2(2)=-CONJG((-ZI*ZADIZ-FB*ZMIXSS(2,IZ)*ZI**THIZ)*COSB-
1105      $           (ZI*ZBDIZ+FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINB)/2.
1106           CALL SSME3(2,AMB1SS,Z1,Z2)
1107           Z1(1)=((ZI*ZADIZ-FB*ZMIXSS(2,IZ)*ZI**THIZ)*SINB+
1108      $           (ZI*ZBDIZ-FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSB)/2.
1109           Z1(2)=((-ZI*ZADIZ-FB*ZMIXSS(2,IZ)*ZI**THIZ)*SINB+
1110      $           (ZI*ZBDIZ+FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSB)/2.
1111           Z2(1)=(COSB+SINB)/2.
1112           Z2(2)=(COSB-SINB)/2.
1113           CALL SSME3(3,AMB2SS,Z1,Z2)
1114           Z1(1)=(COSB+SINB)/2.
1115           Z1(2)=(COSB-SINB)/2.
1116           Z2(1)=CONJG((ZI*ZADIZ-FB*ZMIXSS(2,IZ)*ZI**THIZ)*SINB+
1117      $           (ZI*ZBDIZ-FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSB)/2.
1118           Z2(2)=-CONJG((-ZI*ZADIZ-FB*ZMIXSS(2,IZ)*ZI**THIZ)*SINB+
1119      $           (ZI*ZBDIZ+FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSB)/2.
1120           CALL SSME3(2,AMB2SS,Z1,Z2)
1121         END IF
1122 C          3 body gluino --> q + qb + zi decay, q=t
1123         IF (AMGLSS.GT.(MZIZ+2*AMTP).AND.AMGLSS.LT.
1124      $     (AMT1SS+AMTP)) THEN
1125           TMP(1)=AMGLSS
1126           TMP(2)=AMTP
1127           TMP(3)=MZIZ
1128           TMP(4)=AMT1SS
1129           TMP(5)=AMT1SS
1130           XUPPER=(AMGLSS**2+AMTP**2-(AMTP+MZIZ)**2)/2./AMGLSS
1131           XT1=SSXINT(AMTP,SSGX1,XUPPER)
1132           XT2=SSXINT(AMTP,SSGX2,XUPPER)
1133           XT3=SSXINT(AMTP,SSGX3,XUPPER)
1134           XT4=SSXINT(AMTP,SSGX4,XUPPER)
1135           XT5=SSXINT(AMTP,SSGX5,XUPPER)
1136           XT6=SSXINT(AMTP,SSGX6,XUPPER)
1137           XT7=SSXINT(AMTP,SSGX7,XUPPER)
1138           XT8=SSXINT(AMTP,SSGX8,XUPPER)
1139           XT9=SSXINT(AMTP,SSGX9,XUPPER)
1140           AL1=AUIZ*COST-FT*ZMIXSS(1,IZ)*SINT
1141           BE1=FT*ZMIXSS(1,IZ)*COST+BUIZ*SINT
1142 C         ---- here, al2 is (-) al2 of tata notes-----
1143           AL2=BUIZ*SINT+FT*ZMIXSS(1,IZ)*COST
1144           BE2=-FT*ZMIXSS(1,IZ)*SINT+AUIZ*COST
1145           XLL=(AL1**2+BE1**2)*XT1-4*AMTP*MZIZ*SNIZ*AL1*
1146      $     BE1*XT3-AMGLSS*(SNIZ*MZIZ*(AL1**2*XT2/AMGLSS/
1147      $     MZIZ+BE1**2*AMTP**2*XT5)-AL1*BE1*AMTP*(XT4-
1148      $     MZIZ**2*XT5))
1149           XRR=(AL2**2+BE2**2)*XT1-4*AMTP*MZIZ*SNIZ*AL2*
1150      $     BE2*XT3-AMGLSS*(SNIZ*MZIZ*(AL2**2*XT2/AMGLSS/
1151      $     MZIZ+BE2**2*AMTP**2*XT5)-AL2*BE2*AMTP*(XT4-
1152      $     MZIZ**2*XT5))
1153           XL1R1=-2*AMGLSS*AMTP*((AL1*AL2+BE1*BE2)*SNIZ*AMTP*
1154      $     MZIZ*XT6-(AL2*BE1+AL1*BE2)*XT7)
1155           XL2R2=XL1R1
1156           XL1R2=BE1*BE2*XT8+AL1*AL2*AMTP**2*XT4-AMTP*MZIZ*
1157      $     SNIZ*(AL1*BE2+AL2*BE1)*XT9
1158           XL2R1=XL1R2
1159           XMST1=COST**2*XLL+SINT**2*XRR-SINT*COST*(XL1R1+XL1R2+
1160      $     XL2R1+XL2R2)
1161         ELSE
1162           XMST1=0.
1163         END IF
1164         IF (AMGLSS.GT.(MZIZ+2*AMTP).AND.AMGLSS.LT.
1165      $     (AMT2SS+AMTP)) THEN
1166           TMP(1)=AMGLSS
1167           TMP(2)=AMTP
1168           TMP(3)=MZIZ
1169           TMP(4)=AMT2SS
1170           TMP(5)=AMT2SS
1171           XUPPER=(AMGLSS**2+AMTP**2-(AMTP+MZIZ)**2)/2./AMGLSS
1172           XT1=SSXINT(AMTP,SSGX1,XUPPER)
1173           XT2=SSXINT(AMTP,SSGX2,XUPPER)
1174           XT3=SSXINT(AMTP,SSGX3,XUPPER)
1175           XT4=SSXINT(AMTP,SSGX4,XUPPER)
1176           XT5=SSXINT(AMTP,SSGX5,XUPPER)
1177           XT6=SSXINT(AMTP,SSGX6,XUPPER)
1178           XT7=SSXINT(AMTP,SSGX7,XUPPER)
1179           XT8=SSXINT(AMTP,SSGX8,XUPPER)
1180           XT9=SSXINT(AMTP,SSGX9,XUPPER)
1181           AL1=AUIZ*SINT+FT*ZMIXSS(1,IZ)*COST
1182           BE1=FT*ZMIXSS(1,IZ)*SINT-BUIZ*COST
1183 C         ---- here, al2 is (-) al2 of tata notes-----
1184           AL2=-BUIZ*COST+FT*ZMIXSS(1,IZ)*SINT
1185           BE2=FT*ZMIXSS(1,IZ)*COST+AUIZ*SINT
1186           XLL=(AL1**2+BE1**2)*XT1-4*AMTP*MZIZ*SNIZ*AL1*
1187      $     BE1*XT3-AMGLSS*(SNIZ*MZIZ*(AL1**2*XT2/AMGLSS/
1188      $     MZIZ+BE1**2*AMTP**2*XT5)-AL1*BE1*AMTP*(XT4-
1189      $     MZIZ**2*XT5))
1190           XRR=(AL2**2+BE2**2)*XT1-4*AMTP*MZIZ*SNIZ*AL2*
1191      $     BE2*XT3-AMGLSS*(SNIZ*MZIZ*(AL2**2*XT2/AMGLSS/
1192      $     MZIZ+BE2**2*AMTP**2*XT5)-AL2*BE2*AMTP*(XT4-
1193      $     MZIZ**2*XT5))
1194           XL1R1=-2*AMGLSS*AMTP*((AL1*AL2+BE1*BE2)*SNIZ*AMTP*
1195      $     MZIZ*XT6-(AL2*BE1+AL1*BE2)*XT7)
1196           XL2R2=XL1R1
1197           XL1R2=BE1*BE2*XT8+AL1*AL2*AMTP**2*XT4-AMTP*MZIZ*
1198      $     SNIZ*(AL1*BE2+AL2*BE1)*XT9
1199           XL2R1=XL1R2
1200           XMST2=SINT**2*XLL+COST**2*XRR+SINT*COST*(XL1R1+XL1R2+
1201      $     XL2R1+XL2R2)
1202         ELSE
1203           XMST2=0.
1204         END IF
1205 C       ----cross term between t_1 and t_2 graphs -----------
1206         IF (AMGLSS.GT.(MZIZ+2*AMTP).AND.AMGLSS.LT.
1207      $     (AMT1SS+AMTP)) THEN
1208           TMP(1)=AMGLSS
1209           TMP(2)=AMTP
1210           TMP(3)=MZIZ
1211           TMP(4)=AMT1SS
1212           TMP(5)=AMT2SS
1213           XUPPER=(AMGLSS**2+AMTP**2-(AMTP+MZIZ)**2)/2./AMGLSS
1214           XT1=SSXINT(AMTP,SSGX1,XUPPER)
1215           XT2=SSXINT(AMTP,SSGX2,XUPPER)
1216           XT3=SSXINT(AMTP,SSGX3,XUPPER)
1217           XT4=SSXINT(AMTP,SSGX4,XUPPER)
1218           XT5=SSXINT(AMTP,SSGX5,XUPPER)
1219           XT6=SSXINT(AMTP,SSGX6,XUPPER)
1220           XT7=SSXINT(AMTP,SSGX7,XUPPER)
1221           XT8=SSXINT(AMTP,SSGX8,XUPPER)
1222           XT9=SSXINT(AMTP,SSGX9,XUPPER)
1223           AL1=AUIZ*COST-FT*ZMIXSS(1,IZ)*SINT
1224           AL2=AUIZ*SINT+FT*ZMIXSS(1,IZ)*COST
1225           BE1=FT*ZMIXSS(1,IZ)*COST+BUIZ*SINT
1226           BE2=FT*ZMIXSS(1,IZ)*SINT-BUIZ*COST
1227           XL1L2=COST*SINT*(2*(AL1*AL2+BE1*BE2)*XT1-4*SNIZ*AMTP*
1228      $     MZIZ*(AL1*BE2+AL2*BE1)*XT3-AMGLSS*(2*MZIZ*
1229      $     SNIZ*(AL1*AL2*XT2/AMGLSS/MZIZ+BE1*BE2*AMTP**2*
1230      $     XT5)-(AL1*BE2+AL2*BE1)*AMTP*(XT4-MZIZ**2*XT5)))
1231           AL1=-BUIZ*SINT-FT*ZMIXSS(1,IZ)*COST
1232           AL2=BUIZ*COST-FT*ZMIXSS(1,IZ)*SINT
1233           BE1=-FT*ZMIXSS(1,IZ)*SINT+AUIZ*COST
1234           BE2=FT*ZMIXSS(1,IZ)*COST+AUIZ*SINT
1235           XR1R2=-COST*SINT*(2*(AL1*AL2+BE1*BE2)*XT1+4*SNIZ*AMTP*
1236      $     MZIZ*(AL1*BE2+AL2*BE1)*XT3-AMGLSS*(2*MZIZ*
1237      $     SNIZ*(AL1*AL2*XT2/AMGLSS/MZIZ+BE1*BE2*AMTP**2*
1238      $     XT5)+(AL1*BE2+AL2*BE1)*AMTP*(XT4-MZIZ**2*XT5)))
1239           AL1=AUIZ*COST-FT*ZMIXSS(1,IZ)*SINT
1240           AL2=BUIZ*COST-FT*ZMIXSS(1,IZ)*SINT
1241           BE1=FT*ZMIXSS(1,IZ)*COST+BUIZ*SINT
1242           BE2=FT*ZMIXSS(1,IZ)*COST+AUIZ*SINT
1243           XL1R1=2*AMGLSS*AMTP*COST**2*(SNIZ*(AL1*AL2-BE1*BE2)*
1244      $     AMTP*MZIZ*XT6-(AL2*BE1-AL1*BE2)*XT7)
1245           XL1R2=COST**2*(BE1*BE2*XT8-AL1*AL2*AMTP**2*XT4+AMTP*
1246      $     MZIZ*SNIZ*XT9*(-AL1*BE2+BE1*AL2))
1247           XLR1=2*(XL1R1+XL1R2)
1248           AL1=AUIZ*SINT+FT*ZMIXSS(1,IZ)*COST
1249           AL2=-BUIZ*SINT-FT*ZMIXSS(1,IZ)*COST
1250           BE1=FT*ZMIXSS(1,IZ)*SINT-BUIZ*COST
1251           BE2=-FT*ZMIXSS(1,IZ)*SINT+AUIZ*COST
1252           TMP(4)=AMT2SS
1253           TMP(5)=AMT1SS
1254           XT8=SSXINT(AMTP,SSGX8,XUPPER)
1255           XT9=SSXINT(AMTP,SSGX9,XUPPER)
1256           XL1R1=2*AMGLSS*AMTP*SINT**2*(SNIZ*(-AL1*AL2+BE1*BE2)*
1257      $     AMTP*MZIZ*XT6+(AL2*BE1-AL1*BE2)*XT7)
1258           XL1R2=-SINT**2*(BE1*BE2*XT8-AL1*AL2*AMTP**2*XT4+AMTP*
1259      $     MZIZ*SNIZ*XT9*(-AL1*BE2+BE1*AL2))
1260           XLR2=2*(XL1R1+XL1R2)
1261           XMST12=XL1L2+XR1R2+XLR1+XLR2
1262         ELSE
1263           XMST12=0.
1264         END IF
1265           WID=XTCON*(XMST1+XMST2+XMST12)
1266         IF (WID.GT.0.) THEN
1267           CALL SSSAVE(ISGL,WID,ISZI(IZ),IDTP,-IDTP,0,0)
1268           Z1(1)=((ZI*ZAUIZ-FT*ZMIXSS(1,IZ)*ZI**THIZ)*COST-
1269      $           (ZI*ZBUIZ-FT*ZMIXSS(1,IZ)*(-ZI)**THIZ)*SINT)/2.
1270           Z1(2)=((-ZI*ZAUIZ-FT*ZMIXSS(1,IZ)*ZI**THIZ)*COST-
1271      $           (ZI*ZBUIZ+FT*ZMIXSS(1,IZ)*(-ZI)**THIZ)*SINT)/2.
1272           Z2(1)=(COST-SINT)/2.
1273           Z2(2)=-(COST+SINT)/2.
1274           CALL SSME3(3,AMT1SS,Z1,Z2)
1275           Z1(1)=(COST-SINT)/2.
1276           Z1(2)=-(COST+SINT)/2.
1277           Z2(1)=CONJG((ZI*ZAUIZ-FT*ZMIXSS(1,IZ)*ZI**THIZ)*COST-
1278      $           (ZI*ZBUIZ-FT*ZMIXSS(1,IZ)*(-ZI)**THIZ)*SINT)/2.
1279           Z2(2)=-CONJG((-ZI*ZAUIZ-FT*ZMIXSS(1,IZ)*ZI**THIZ)*COST-
1280      $           (ZI*ZBUIZ+FT*ZMIXSS(1,IZ)*(-ZI)**THIZ)*SINT)/2.
1281           CALL SSME3(2,AMT1SS,Z1,Z2)
1282           Z1(1)=((ZI*ZAUIZ-FT*ZMIXSS(1,IZ)*ZI**THIZ)*SINT+
1283      $           (ZI*ZBUIZ-FT*ZMIXSS(1,IZ)*(-ZI)**THIZ)*COST)/2.
1284           Z1(2)=((-ZI*ZAUIZ-FT*ZMIXSS(1,IZ)*ZI**THIZ)*SINT+
1285      $           (ZI*ZBUIZ+FT*ZMIXSS(1,IZ)*(-ZI)**THIZ)*COST)/2.
1286           Z2(1)=(COST+SINT)/2.
1287           Z2(2)=(COST-SINT)/2.
1288           CALL SSME3(3,AMT2SS,Z1,Z2)
1289           Z1(1)=(COST+SINT)/2.
1290           Z1(2)=(COST-SINT)/2.
1291           Z2(1)=CONJG((ZI*ZAUIZ-FT*ZMIXSS(1,IZ)*ZI**THIZ)*SINT+
1292      $           (ZI*ZBUIZ-FT*ZMIXSS(1,IZ)*(-ZI)**THIZ)*COST)/2.
1293           Z2(2)=-CONJG((-ZI*ZAUIZ-FT*ZMIXSS(1,IZ)*ZI**THIZ)*SINT+
1294      $           (ZI*ZBUIZ+FT*ZMIXSS(1,IZ)*(-ZI)**THIZ)*COST)/2.
1295           CALL SSME3(2,AMT2SS,Z1,Z2)
1296         END IF
1297 100   CONTINUE
1298 C
1299 C          gluino --> quark + squark mode
1300 C
1301       IF (AMGLSS.GT.(AMULSS+AMUP)) THEN
1302         GMQK=ALF3*AMGLSS*(1.+AMUP**2/AMGLSS**2-AMULSS**2/AMGLSS**2)*
1303      $  SQRT(SSXLAM(1.,AMUP**2/AMGLSS**2,AMULSS**2/AMGLSS**2))/8.
1304         CALL SSSAVE(ISGL,GMQK,-ISUPL,+IDUP,0,0,0)
1305         CALL SSSAVE(ISGL,GMQK,+ISUPL,-IDUP,0,0,0)
1306       END IF
1307       IF (AMGLSS.GT.(AMDLSS+AMDN)) THEN
1308         GMQK=ALF3*AMGLSS*(1.+AMDN**2/AMGLSS**2-AMDLSS**2/AMGLSS**2)*
1309      $  SQRT(SSXLAM(1.,AMDN**2/AMGLSS**2,AMDLSS**2/AMGLSS**2))/8.
1310         CALL SSSAVE(ISGL,GMQK,-ISDNL,+IDDN,0,0,0)
1311         CALL SSSAVE(ISGL,GMQK,+ISDNL,-IDDN,0,0,0)
1312       END IF
1313       IF (AMGLSS.GT.(AMURSS+AMUP)) THEN
1314         GMQK=ALF3*AMGLSS*(1.+AMUP**2/AMGLSS**2-AMURSS**2/AMGLSS**2)*
1315      $  SQRT(SSXLAM(1.,AMUP**2/AMGLSS**2,AMURSS**2/AMGLSS**2))/8.
1316         CALL SSSAVE(ISGL,GMQK,-ISUPR,+IDUP,0,0,0)
1317         CALL SSSAVE(ISGL,GMQK,+ISUPR,-IDUP,0,0,0)
1318       END IF
1319       IF (AMGLSS.GT.(AMDRSS+AMDN)) THEN
1320         GMQK=ALF3*AMGLSS*(1.+AMDN**2/AMGLSS**2-AMDRSS**2/AMGLSS**2)*
1321      $  SQRT(SSXLAM(1.,AMDN**2/AMGLSS**2,AMDRSS**2/AMGLSS**2))/8.
1322         CALL SSSAVE(ISGL,GMQK,-ISDNR,+IDDN,0,0,0)
1323         CALL SSSAVE(ISGL,GMQK,+ISDNR,-IDDN,0,0,0)
1324       END IF
1325 C
1326       IF (AMGLSS.GT.(AMSLSS+AMST)) THEN
1327         GMQK=ALF3*AMGLSS*(1.+AMST**2/AMGLSS**2-AMSLSS**2/AMGLSS**2)*
1328      $  SQRT(SSXLAM(1.,AMST**2/AMGLSS**2,AMSLSS**2/AMGLSS**2))/8.
1329         CALL SSSAVE(ISGL,GMQK,-ISSTL,+IDST,0,0,0)
1330         CALL SSSAVE(ISGL,GMQK,+ISSTL,-IDST,0,0,0)
1331       END IF
1332       IF (AMGLSS.GT.(AMSRSS+AMST)) THEN
1333         GMQK=ALF3*AMGLSS*(1.+AMST**2/AMGLSS**2-AMSRSS**2/AMGLSS**2)*
1334      $  SQRT(SSXLAM(1.,AMST**2/AMGLSS**2,AMSRSS**2/AMGLSS**2))/8.
1335         CALL SSSAVE(ISGL,GMQK,-ISSTR,+IDST,0,0,0)
1336         CALL SSSAVE(ISGL,GMQK,+ISSTR,-IDST,0,0,0)
1337       END IF
1338 C
1339       IF (AMGLSS.GT.(AMCLSS+AMCH)) THEN
1340         GMQK=ALF3*AMGLSS*(1.+AMCH**2/AMGLSS**2-AMCLSS**2/AMGLSS**2)*
1341      $  SQRT(SSXLAM(1.,AMCH**2/AMGLSS**2,AMCLSS**2/AMGLSS**2))/8.
1342         CALL SSSAVE(ISGL,GMQK,-ISCHL,+IDCH,0,0,0)
1343         CALL SSSAVE(ISGL,GMQK,+ISCHL,-IDCH,0,0,0)
1344       END IF
1345       IF (AMGLSS.GT.(AMCRSS+AMCH)) THEN
1346         GMQK=ALF3*AMGLSS*(1.+AMCH**2/AMGLSS**2-AMCRSS**2/AMGLSS**2)*
1347      $  SQRT(SSXLAM(1.,AMCH**2/AMGLSS**2,AMCRSS**2/AMGLSS**2))/8.
1348         CALL SSSAVE(ISGL,GMQK,-ISCHR,+IDCH,0,0,0)
1349         CALL SSSAVE(ISGL,GMQK,+ISCHR,-IDCH,0,0,0)
1350       END IF
1351 C
1352 C     !!! NEED MIXING ANGLE PIECE LIKE STOPS
1353       IF (AMGLSS.GT.(AMB1SS+AMBT)) THEN
1354         GMQK=ALF3*AMGLSS*(1.+AMBT**2/AMGLSS**2-AMB1SS**2/AMGLSS**2)*
1355      $  SQRT(SSXLAM(1.,AMBT**2/AMGLSS**2,AMB1SS**2/AMGLSS**2))/8.
1356         CALL SSSAVE(ISGL,GMQK,-ISBT1,+IDBT,0,0,0)
1357         CALL SSSAVE(ISGL,GMQK,+ISBT1,-IDBT,0,0,0)
1358       END IF
1359 C
1360       IF (AMGLSS.GT.(AMB2SS+AMBT)) THEN
1361         GMQK=ALF3*AMGLSS*(1.+AMBT**2/AMGLSS**2-AMB2SS**2/AMGLSS**2)*
1362      $  SQRT(SSXLAM(1.,AMBT**2/AMGLSS**2,AMB2SS**2/AMGLSS**2))/8.
1363         CALL SSSAVE(ISGL,GMQK,-ISBT2,+IDBT,0,0,0)
1364         CALL SSSAVE(ISGL,GMQK,+ISBT2,-IDBT,0,0,0)
1365       END IF
1366 C
1367       IF (AMGLSS.GT.(AMT1SS+AMTP)) THEN
1368         GMQK=ALF3*AMGLSS*(1.+AMTP**2/AMGLSS**2-AMT1SS**2/AMGLSS**2-
1369      $   2*SIN(2*THETAT)*AMTP/AMGLSS)*
1370      $  SQRT(SSXLAM(1.,AMTP**2/AMGLSS**2,AMT1SS**2/AMGLSS**2))/8.
1371         CALL SSSAVE(ISGL,GMQK,-ISTP1,+IDTP,0,0,0)
1372         CALL SSSAVE(ISGL,GMQK,+ISTP1,-IDTP,0,0,0)
1373       END IF
1374 C
1375       IF (AMGLSS.GT.(AMT2SS+AMTP)) THEN
1376         GMQK=ALF3*AMGLSS*(1.+AMTP**2/AMGLSS**2-AMT2SS**2/AMGLSS**2+
1377      $   2*SIN(2*THETAT)*AMTP/AMGLSS)*
1378      $  SQRT(SSXLAM(1.,AMTP**2/AMGLSS**2,AMT2SS**2/AMGLSS**2))/8.
1379         CALL SSSAVE(ISGL,GMQK,-ISTP2,+IDTP,0,0,0)
1380         CALL SSSAVE(ISGL,GMQK,+ISTP2,-IDTP,0,0,0)
1381       END IF
1382 C
1383 C     Decay to gravitino
1384 C
1385       IF (AMGLSS.GT.AMGVSS) THEN
1386         WID=AMGLSS**5/48./PI/(AMGVSS*AMPL)**2
1387         CALL SSSAVE(ISGL,WID,91,IDGL,0,0,0)
1388       END IF
1389 C
1390 C          Normalize branching ratios
1391 C
1392       CALL SSNORM(ISGL)
1393 C
1394       RETURN
1395       END