ISAJET/isasusy/sszibf.F

   1 #include "isajet/pilot.h"
   2         SUBROUTINE SSZIBF
   3 C-----------------------------------------------------------------------
   4 C       This subroutine calculates the neutralino (zi) subset of
   5 C       SSWZBF, which was too long.
   6 C       Valid for all scalar masses (functions in double precision)
   7 C       Includes Higgs sector radiative corrections (Aug. 31)
   8 C
   9 C       Auxiliary functions are called SSWxyi, SSZxyi, where normally
  10 C       x indicates the SUSY particle, y the SM particle(s), and i is
  11 C       a counter.
  12 C
  13 C       Part of Baer's GAUGBF
  14 C
  15 C-----------------------------------------------------------------------
  16 #if defined(CERNLIB_IMPNONE)
  17       IMPLICIT NONE
  18 #endif
  19 #include "isajet/sslun.inc"
  20 #include "isajet/ssmode.inc"
  21 #include "isajet/sssm.inc"
  22 #include "isajet/sspar.inc"
  23 #include "isajet/sstype.inc"
  24 #include "isajet/sstmp.inc"
  25 #include "isajet/sspols.inc"
  26 C
  27       EXTERNAL SSZWF1,SSZZF1,SSZZF2,SSZZF3,SSWZF1,SSWZF2,SSWZF3
  28      $,SSWZF4,SSWZF5,SSWZF6,SSWZF7,SSWWF1,SSZZF4,SSZZF5,SSGX1
  29      $,SSGX2,SSGX8
  30 C
  31       REAL AUI(4),BUI(4),ADI(4),BDI(4),ALI(4),BLI(4),ANI(4),BNI(4)
  32      $,WIJ(4,4),AUWI(4),ADWI(4),ANWI(4),ALWI(4),XIM(4),YIM(4)
  33      $,XIP(4),YIP(4),SNIJ(4,4),XLIJ(4,4),HIJ(4,4)
  34      $,V1I(4),V2I(4),V3I(4),V4I(4),XHIJ(4,4),XPIJ(4,4),AMWISS(2)
  35 C
  36       INTEGER ISZ(4),THIZ,THJZ
  37 C
  38       REAL MWIW,SL,PP,SP,PL,MZIZ,PH,SH,COSB,SINB,COSA,SINA
  39      $,MZJZ,FAC,COSBE,SINBE
  40       REAL STHW,CTHW
  41       REAL EF,A,B,TANB,FB,FT,SR2,G,PI,GP,FL
  42      $,MW2,SNW1,MW1,YM,BE,SNW2,XM,THX,THY
  43      $,BTN,APD,APL,APU,BTD,BTL,APN,BTU
  44       REAL TANW,COTW,XWINO,YWINO,SN,SNIW,SNJZ
  45       REAL SSXINT,SSXLAM
  46       REAL WID,WIDEE,E,
  47      $TERMLL,TERMRR,TERMZZ,TERMLZ,TERMRZ,TERMLH
  48       REAL FACT,ALIZ1,ALJZ1,ALIZ2,ALJZ2,BEIZ1,BEJZ1,BEIZ2,BEJZ2,
  49      $SGNIJ,XUPPER,GLLF1,GRRF1,GLRF1,GF1,GLLF2,GRRF2,GLRF2,GF2,
  50      $GLGL,GRGR,GLGR,GRGL,GF12,GF,TERMHL,TERMHH,TERMHA,
  51      $TERM1Z,TERM2Z,TERM1L,TERM2L,TERM1H,TERM2H,TERM1A,TERM2A
  52       REAL COSL,SINL,BPWI(2),BPLWI(2)
  53       REAL BWI(2),AS,BS,COST,SINT,SNIZ
  54       REAL GLLF1L,GRRF1R,GLGLL,GRGRR,GLLF2L,GRRF2R,TMZZRL,TMZZLR,
  55      $TM1ZRL,TM1ZLR,TM2ZRL,TM2ZLR
  56       REAL ASMB,MBMB,MBQ,ASMT,MTMT,MTQ,SUALFS
  57       REAL SUALFE,MTAMTA,MTAMB,MTAMZ,AMPL
  58       REAL FUDGE
  59       DOUBLE PRECISION SSMQCD
  60       COMPLEX ZI,ZONE,ZA,ZB,ZPP,ZPM,ZAUIZ,ZBUIZ,ZADIZ,ZBDIZ,
  61      $ZALIZ,ZBLIZ,Z1(2),Z2(2)
  62       INTEGER IW,JZ,IZ,ISZIZ,ISWIW
  63       DATA FUDGE/1.0/
  64       DATA ZONE/(1.,0.)/,ZI/(0.,1.)/
  65 C
  66 C          Constants from neutralino mass matrix
  67 C
  68       AMPL=2.4E18
  69       PI=4.*ATAN(1.)
  70       SR2=SQRT(2.)
  71       G=SQRT(4*PI*ALFAEM/SN2THW)
  72       GP=G*SQRT(SN2THW/(1.-SN2THW))
  73       E=SQRT(4*PI/128.)
  74 C
  75       TANW=SQRT(SN2THW/(1.-SN2THW))
  76       COTW=1./TANW
  77       STHW=SQRT(SN2THW)
  78       CTHW=SQRT(1.-SN2THW)
  79       APL=.25*(3*TANW-COTW)
  80       BTL=.25*(COTW+TANW)
  81       APN=.25*(TANW+COTW)
  82       BTN=-.25*(COTW+TANW)
  83       APU=-5*TANW/12.+COTW/4.
  84       BTU=-.25*(COTW+TANW)
  85       APD=-COTW/4.+TANW/12.
  86       BTD=.25*(COTW+TANW)
  87 C
  88       TANB=1./RV2V1
  89       BE=ATAN(1./RV2V1)
  90       SINBE=SIN(BE)
  91       COSBE=COS(BE)
  92       XM=1./TAN(GAMMAL)
  93       YM=1./TAN(GAMMAR)
  94       THX=SIGN(1.,XM)
  95       THY=SIGN(1.,YM)
  96       ASMB=SUALFS(AMBT**2,.36,AMTP,3)
  97       MBMB=AMBT*(1.-4*ASMB/3./PI)
  98       MBQ=SSMQCD(DBLE(MBMB),DBLE(AMTP))
  99       ASMT=SUALFS(AMTP**2,.36,AMTP,3)
 100       MTMT=AMTP/(1.+4*ASMT/3./PI+(16.11-1.04*(5.-6.63/AMTP))*
 101      $(ASMT/PI)**2)
 102       MTQ=SSMQCD(DBLE(MTMT),DBLE(AMTP))
 103       FB=G*MBQ/SR2/AMW/COS(BE)
 104       FT=G*MTQ/SR2/AMW/SIN(BE)
 105       MTAMTA=AMTAU*(1.-SUALFE(AMTAU**2)/PI)
 106       MTAMB=MTAMTA*(SUALFE(AMBT**2)/SUALFE(AMTAU**2))**(-27./76.)
 107       MTAMZ=MTAMB*(SUALFE(AMZ**2)/SUALFE(AMBT**2))**(-27./80.)
 108       FL=G*MTAMZ/SR2/AMW/COS(BE)
 109       SNW1=SIGN(1.,AMW1SS)
 110       SNW2=SIGN(1.,AMW2SS)
 111       AMWISS(1)=AMW1SS
 112       AMWISS(2)=AMW2SS
 113       BWI(1)=-FT*SNW1*COS(GAMMAR)
 114       BWI(2)=FT*SNW2*THY*SIN(GAMMAR)
 115       BPWI(1)=-FB*COS(GAMMAL)
 116       BPWI(2)=FB*THX*SIN(GAMMAL)
 117       BPLWI(1)=-FL*COS(GAMMAL)
 118       BPLWI(2)=FL*THX*SIN(GAMMAL)
 119       MW1=ABS(AMW1SS)
 120       MW2=ABS(AMW2SS)
 121       XWINO=.5*(THX*SIN(GAMMAL)*COS(GAMMAL)
 122      $-THY*SIN(GAMMAR)*COS(GAMMAR))
 123       YWINO=.5*(THX*SIN(GAMMAL)*COS(GAMMAL)
 124      $+THY*SIN(GAMMAR)*COS(GAMMAR))
 125       COST=COS(THETAT)
 126       SINT=SIN(THETAT)
 127       COSB=COS(THETAB)
 128       SINB=SIN(THETAB)
 129       COSL=COS(THETAL)
 130       SINL=SIN(THETAL)
 131 C
 132 C          Constants from Higgs mass matrix
 133 C
 134       SINA=SIN(ALFAH)
 135       COSA=COS(ALFAH)
 136 C
 137 C          Gaugino couplings
 138 C
 139       DO 100 IZ=1,4
 140         AUI(IZ)=G/SR2*ZMIXSS(3,IZ)+GP/3./SR2*ZMIXSS(4,IZ)
 141         BUI(IZ)=4.*GP/3./SR2*ZMIXSS(4,IZ)
 142         ADI(IZ)=-G/SR2*ZMIXSS(3,IZ)+GP/3./SR2*ZMIXSS(4,IZ)
 143         BDI(IZ)=-2.*GP/3./SR2*ZMIXSS(4,IZ)
 144         ALI(IZ)=G/SR2*ZMIXSS(3,IZ)+GP/SR2*ZMIXSS(4,IZ)
 145         BLI(IZ)=-SR2*GP*ZMIXSS(4,IZ)
 146         ANI(IZ)=G/SR2*ZMIXSS(3,IZ)-GP/SR2*ZMIXSS(4,IZ)
 147         BNI(IZ)=0.0
 148 100   CONTINUE
 149 C
 150       DO 110 IZ=1,4
 151       DO 110 JZ=1,4
 152         IF(IZ.LT.JZ) THEN
 153           WIJ(IZ,JZ)=SQRT(G**2+GP**2)
 154      $    *(ZMIXSS(1,IZ)*ZMIXSS(1,JZ)-ZMIXSS(2,IZ)*ZMIXSS(2,JZ))/4.
 155         ELSEIF(IZ.GT.JZ) THEN
 156           WIJ(IZ,JZ)=-SQRT(G**2+GP**2)
 157      $    *(ZMIXSS(1,IZ)*ZMIXSS(1,JZ)-ZMIXSS(2,IZ)*ZMIXSS(2,JZ))/4.
 158         ENDIF
 159 110   CONTINUE
 160 C
 161       AUWI(2)=G*THX*COS(GAMMAL)
 162       ADWI(2)=SNW2*G*THY*COS(GAMMAR)
 163       ALWI(2)=ADWI(2)
 164       ANWI(2)=AUWI(2)
 165       AUWI(1)=G*SIN(GAMMAL)
 166       ADWI(1)=SNW1*G*SIN(GAMMAR)
 167       ALWI(1)=ADWI(1)
 168       ANWI(1)=AUWI(1)
 169 C
 170       DO 120 IZ=1,4
 171         XIM(IZ)=.5*(SNW1*SIGN(1.,AMZISS(IZ))*(COS(GAMMAR)
 172      $  *ZMIXSS(1,IZ)/SR2+SIN(GAMMAR)*ZMIXSS(3,IZ))-COS(GAMMAL)
 173      $  *ZMIXSS(2,IZ)/SR2+SIN(GAMMAL)*ZMIXSS(3,IZ))
 174         YIM(IZ)=.5*(-SNW1*SIGN(1.,AMZISS(IZ))*(COS(GAMMAR)
 175      $  *ZMIXSS(1,IZ)/SR2+SIN(GAMMAR)*ZMIXSS(3,IZ))-COS(GAMMAL)
 176      $  *ZMIXSS(2,IZ)/SR2+SIN(GAMMAL)*ZMIXSS(3,IZ))
 177         XIP(IZ)=.5*(SNW2*SIGN(1.,AMZISS(IZ))*THY*(-SIN(GAMMAR)
 178      $  *ZMIXSS(1,IZ)/SR2+COS(GAMMAR)*ZMIXSS(3,IZ))+THX*(SIN(GAMMAL)
 179      $  *ZMIXSS(2,IZ)/SR2+COS(GAMMAL)*ZMIXSS(3,IZ)))
 180         YIP(IZ)=.5*(-SNW2*SIGN(1.,AMZISS(IZ))*THY*(-SIN(GAMMAR)
 181      $  *ZMIXSS(1,IZ)/SR2+COS(GAMMAR)*ZMIXSS(3,IZ))+THX*(SIN(GAMMAL)
 182      $  *ZMIXSS(2,IZ)/SR2+COS(GAMMAL)*ZMIXSS(3,IZ)))
 183 120   CONTINUE
 184 C
 185       DO 130 IZ=1,4
 186       DO 130 JZ=1,4
 187         IF(IZ.NE.JZ) THEN
 188           SNIJ(IZ,JZ)=-1.*SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ))
 189           XLIJ(IZ,JZ)=-SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ))
 190      $    *(ZMIXSS(2,IZ)*SINA-ZMIXSS(1,IZ)*COSA)
 191      $    *(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/2.
 192           XHIJ(IZ,JZ)=-SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ))
 193      $    *(ZMIXSS(2,IZ)*COSA+ZMIXSS(1,IZ)*SINA)
 194      $    *(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/2.
 195           XPIJ(IZ,JZ)=SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ))
 196      $    *(ZMIXSS(2,IZ)*SINBE-ZMIXSS(1,IZ)*COSBE)
 197      $    *(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/2.
 198           HIJ(IZ,JZ)=-SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ))
 199      $    *(ZMIXSS(2,IZ)*COSA+ZMIXSS(1,IZ)*SINA)
 200      $    *(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/2.
 201         ENDIF
 202 130   CONTINUE
 203 C
 204       SP=-.5*(-THY*SNW2*COSBE*SIN(GAMMAL)*SIN(GAMMAR)+
 205      $THY*SNW2*SINBE*COS(GAMMAL)*COS(GAMMAR)-
 206      $THX*SNW1*COSBE*COS(GAMMAL)*COS(GAMMAR)+
 207      $THX*SNW1*SINBE*SIN(GAMMAL)*SIN(GAMMAR))
 208       PP=-.5*(-THY*SNW2*COSBE*SIN(GAMMAL)*SIN(GAMMAR)+
 209      $THY*SNW2*SINBE*COS(GAMMAL)*COS(GAMMAR)+
 210      $THX*SNW1*COSBE*COS(GAMMAL)*COS(GAMMAR)-
 211      $THX*SNW1*SINBE*SIN(GAMMAL)*SIN(GAMMAR))
 212 C
 213       SL=.5*(THY*SNW2*SINA*COS(GAMMAL)*COS(GAMMAR)-
 214      $THY*SNW2*COSA*SIN(GAMMAL)*SIN(GAMMAR)+
 215      $THX*SNW1*COSA*COS(GAMMAL)*COS(GAMMAR)-
 216      $THX*SNW1*SINA*SIN(GAMMAL)*SIN(GAMMAR))
 217       PL=.5*(THY*SNW2*SINA*COS(GAMMAL)*COS(GAMMAR)-
 218      $THY*SNW2*COSA*SIN(GAMMAL)*SIN(GAMMAR)-
 219      $THX*SNW1*COSA*COS(GAMMAL)*COS(GAMMAR)+
 220      $THX*SNW1*SINA*SIN(GAMMAL)*SIN(GAMMAR))
 221 C
 222       SH=.5*(THY*SNW2*COSA*COS(GAMMAL)*COS(GAMMAR)+
 223      $THY*SNW2*SINA*SIN(GAMMAL)*SIN(GAMMAR)-
 224      $THX*SNW1*SINA*COS(GAMMAL)*COS(GAMMAR)-
 225      $THX*SNW1*COSA*SIN(GAMMAL)*SIN(GAMMAR))
 226       PH=.5*(THY*SNW2*COSA*COS(GAMMAL)*COS(GAMMAR)+
 227      $THY*SNW2*SINA*SIN(GAMMAL)*SIN(GAMMAR)+
 228      $THX*SNW1*SINA*COS(GAMMAL)*COS(GAMMAR)+
 229      $THX*SNW1*COSA*SIN(GAMMAL)*SIN(GAMMAR))
 230 C
 231       DO 140 IZ=1,4
 232         V1I(IZ)=-SIN(GAMMAR)/SR2*(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))
 233      $  -G*COS(GAMMAR)*ZMIXSS(1,IZ)
 234         V2I(IZ)=COS(GAMMAR)/SR2*(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))
 235      $  -G*SIN(GAMMAR)*ZMIXSS(1,IZ)
 236         V3I(IZ)=-SIN(GAMMAL)/SR2*(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))
 237      $  +G*COS(GAMMAL)*ZMIXSS(2,IZ)
 238         V4I(IZ)=COS(GAMMAL)/SR2*(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))
 239      $  +G*SIN(GAMMAL)*ZMIXSS(2,IZ)
 240 140   CONTINUE
 241 C-----------------------------------------------------------------------
 242 C          Generate Neutralino zi Branching Fractions
 243 C-----------------------------------------------------------------------
 244       ISZ(1)=ISZ1
 245       ISZ(2)=ISZ2
 246       ISZ(3)=ISZ3
 247       ISZ(4)=ISZ4
 248 C
 249       DO 200 IZ=2,4
 250         MZIZ=ABS(AMZISS(IZ))
 251         SNIZ=SIGN(1.,AMZISS(IZ))
 252         IF (SNIZ.EQ.1.) THEN
 253            THIZ=0
 254         ELSE
 255            THIZ=1
 256         END IF
 257         ISZIZ=ISZ(IZ)
 258 C          z2 --> z1 + photon
 259 C          !!! NEEDS UPDATING
 260         IF (IZ.EQ.2.AND.(ABS(ZMIXSS(3,1)).LE.1.E-4).AND.
 261      $  (ABS(ZMIXSS(4,1)).LE.1.E-4)) THEN
 262           CALL SSSAVE(ISZ2,1.,ISZ1,IDGM,0,0,0)
 263           GOTO 200
 264         END IF
 265 C
 266 C          zi --> wi + x partial widths
 267 C
 268         DO 205 IW=1,2
 269 C          Loop over w1, w2
 270           IF(IW.EQ.1) THEN
 271             MWIW=MW1
 272             SNIW=SNW1
 273             ISWIW=ISW1
 274           ELSE
 275             MWIW=MW2
 276             SNIW=SNW2
 277             ISWIW=ISW2
 278           ENDIF
 279 C
 280 C          zi -> wj + w
 281 C
 282           IF(MZIZ.GT.MWIW+AMW) THEN
 283             EF=MZIZ**2+MWIW**2-AMW**2+((MZIZ**2-MWIW**2)**2
 284      $         -AMW**4)/AMW/AMW
 285             WID=G*G*SQRT(SSXLAM(MZIZ**2,AMW**2,MWIW**2))/32./PI/
 286      $          MZIZ**3*(2.*EF*(XIM(IZ)**2+YIM(IZ)**2)-12.*
 287      $          MZIZ*MWIW*(XIM(IZ)**2-YIM(IZ)**2))
 288             CALL SSSAVE(ISZIZ,WID,+ISWIW,-IDW,0,0,0)
 289             CALL SSSAVE(ISZIZ,WID,-ISWIW,+IDW,0,0,0)
 290 C
 291 C          zi -> wj + f + fbar (w forbidden) ONLY W CONTRIBUTION INCLUDED!
 292 C
 293           ELSEIF(MZIZ.GT.FUDGE*MWIW) THEN
 294             TMP(1)=MWIW
 295             TMP(2)=MZIZ
 296             IF(IW.EQ.1) THEN
 297               TMP(3)=XIM(IZ)
 298               TMP(4)=YIM(IZ)
 299             ELSE
 300               TMP(3)=XIP(IZ)
 301               TMP(4)=YIP(IZ)
 302             ENDIF
 303             WID=G**4/96./PI**3/MZIZ
 304      $      *SSXINT(MWIW,SSZWF1,(MWIW**2+MZIZ**2)/2./MZIZ)
 305             CALL SSSAVE(ISZIZ,3.*WID,-ISWIW,IDUP,-IDDN,0,0)
 306             Z1(1)=(-ZI)**THIZ*G*TMP(3)
 307             Z1(2)=(-ZI)**THIZ*G*TMP(4)
 308             Z2(1)=G/2./SR2
 309             Z2(2)=-Z2(1)
 310             CALL SSME3(1,AMW,Z1,Z2)
 311             CALL SSSAVE(ISZIZ,WID,-ISWIW,IDNE,-IDE,0,0)
 312             CALL SSME3(1,AMW,Z1,Z2)
 313             CALL SSSAVE(ISZIZ,WID,-ISWIW,IDNM,-IDMU,0,0)
 314             CALL SSME3(1,AMW,Z1,Z2)
 315             CALL SSSAVE(ISZIZ,3.*WID,ISWIW,IDDN,-IDUP,0,0)
 316             CALL SSME3(1,AMW,Z1,Z2)
 317             CALL SSSAVE(ISZIZ,WID,ISWIW,IDE,-IDNE,0,0)
 318             CALL SSME3(1,AMW,Z1,Z2)
 319             CALL SSSAVE(ISZIZ,WID,ISWIW,IDMU,-IDNM,0,0)
 320             CALL SSME3(1,AMW,Z1,Z2)
 321             IF (MZIZ.GT.(MWIW+AMCH+AMST)) THEN
 322               CALL SSSAVE(ISZIZ,3.*WID,-ISWIW,IDCH,-IDST,0,0)
 323               CALL SSME3(1,AMW,Z1,Z2)
 324               CALL SSSAVE(ISZIZ,WID,-ISWIW,IDNT,-IDTAU,0,0)
 325               CALL SSME3(1,AMW,Z1,Z2)
 326               CALL SSSAVE(ISZIZ,3.*WID,ISWIW,IDST,-IDCH,0,0)
 327               CALL SSME3(1,AMW,Z1,Z2)
 328               CALL SSSAVE(ISZIZ,WID,ISWIW,IDTAU,-IDNT,0,0)
 329               CALL SSME3(1,AMW,Z1,Z2)
 330             ENDIF
 331           ENDIF
 332 C
 333 C          zi --> wj + hc
 334 C
 335           IF (MZIZ.GT.(MWIW+AMHC)) THEN
 336           IF (IW.EQ.1) THEN
 337             A=(SNW1*COSBE*V2I(IZ)-SIGN(1.,AMZISS(IZ))
 338      $      *SINBE*V4I(IZ))/2.
 339             B=(SNW1*COSBE*V2I(IZ)+SIGN(1.,AMZISS(IZ))
 340      $      *SINBE*V4I(IZ))/2.
 341           ELSE
 342             A=(THY*SNW2*COSBE*V1I(IZ)-SIGN(1.,AMZISS(IZ))
 343      $      *THX*SINBE*V3I(IZ))/2.
 344             B=(THY*SNW2*COSBE*V1I(IZ)+SIGN(1.,AMZISS(IZ))
 345      $      *THX*SINBE*V3I(IZ))/2.
 346           END IF
 347             WID=SQRT(MWIW**4+MZIZ**4+AMHC**4-2.*(MWIW*MZIZ)**2
 348      $      -2*(MWIW*AMHC)**2-2*(MZIZ*AMHC)**2)/8./PI/MZIZ**3
 349      $      *((A*A+B*B)*(MWIW*MWIW+MZIZ*MZIZ-AMHC*AMHC)/2.
 350      $      +(A*A-B*B)*MWIW*MZIZ)
 351             CALL SSSAVE(ISZIZ,WID,+ISWIW,-ISHC,0,0,0)
 352             CALL SSSAVE(ISZIZ,WID,-ISWIW,+ISHC,0,0,0)
 353           ENDIF
 354 205     CONTINUE
 355 C
 356 C          zi --> zj + z
 357 C
 358 C          Note that if m(zi) > m(zj) + m(z), then the z terms are
 359 C          omitted from the zi -> zj + f + fbar calculation, so there
 360 C          is no double counting with zi -> zj + z.
 361 C
 362         DO 210 JZ=1,IZ-1
 363           MZJZ=ABS(AMZISS(JZ))
 364           IF(MZIZ.GT.(MZJZ+AMZ)) THEN
 365             SN=-1.*SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ))
 366             WID=WIJ(JZ,IZ)**2/(2.*PI)/(MZIZ)**3*SQRT(MZIZ**4+MZJZ**4
 367      $      +AMZ**4-2.*(MZIZ*MZJZ)**2-2.*(MZIZ*AMZ)**2-2.*(MZJZ*AMZ)**2)
 368      $      *((MZIZ**2+MZJZ**2-AMZ**2)/2.+((MZIZ**2-MZJZ**2)**2-AMZ**4)/
 369      $      2./AMZ**2-3.*SN*MZIZ*MZJZ)
 370             CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDZ,0,0,0)
 371           END IF
 372 210     CONTINUE
 373 C
 374 C          zi --> zj + f + fbar
 375 C
 376         DO 220 JZ=1,IZ-1
 377           MZJZ=ABS(AMZISS(JZ))
 378           SNJZ=SIGN(1.,AMZISS(JZ))
 379           IF (SNJZ.EQ.1.) THEN
 380            THJZ=0
 381           ELSE
 382            THJZ=1
 383           END IF
 384           IF(MZIZ.LT.FUDGE*MZJZ) GO TO 220
 385           FAC=1./2./MZIZ/(2.*PI)**5*PI**2*MZIZ**2
 386 C          Leptons -- Z decay allowed, so omit Z
 387             TMP(1)=MZIZ
 388             TMP(2)=MZJZ
 389             TMP(4)=-SNIJ(JZ,IZ)
 390             TMP(6)=0.
 391 C          zi -> zj + u + ubar
 392           IF (MZIZ.GT.(MZJZ+2*AMUP)) THEN
 393             IF (MZIZ.LT.AMULSS) THEN
 394             TMP(3)=AMULSS
 395             TERMLL=3*2*AUI(JZ)**2*AUI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.)
 396             ELSE
 397             TERMLL=0.
 398             END IF
 399             IF (MZIZ.LT.AMURSS) THEN
 400             TMP(3)=AMURSS
 401             TERMRR=3*2*BUI(JZ)**2*BUI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.)
 402             ELSE
 403             TERMRR=0.
 404             END IF
 405             IF (MZIZ.LT.(MZJZ+AMZ)) THEN
 406             TMP(4)=+SNIJ(JZ,IZ)
 407             TERMZZ=3*(APU**2+BTU**2)*E**2*WIJ(JZ,IZ)**2/MZIZ
 408      $             *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ)
 409             ELSE
 410             TERMZZ=0.
 411             END IF
 412             IF (MZIZ.LT.AMULSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
 413             TMP(4)=+SNIJ(JZ,IZ)
 414             TMP(3)=AMULSS
 415             TERMLZ=3*8*E*(APU-BTU)*WIJ(JZ,IZ)*AUI(JZ)*AUI(IZ)/MZIZ
 416      $      /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2)
 417             ELSE
 418             TERMLZ=0.
 419             END IF
 420             IF (MZIZ.LT.AMURSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
 421             TMP(4)=+SNIJ(JZ,IZ)
 422             TMP(3)=AMURSS
 423             TERMRZ=-3*8*E*(APU+BTU)*WIJ(JZ,IZ)*BUI(JZ)*BUI(IZ)/MZIZ
 424      $      /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2)
 425             ELSE
 426             TERMRZ=0.
 427             END IF
 428             WID=TERMLL+TERMRR+TERMZZ+TERMLZ+TERMRZ
 429             CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDUP,-IDUP,0,0)
 430 C     Enter information for decay matrix element
 431             IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN
 432              Z1(1)=0.
 433              Z1(2)=2*E*WIJ(JZ,IZ)
 434             ELSE
 435              Z1(1)=2*E*WIJ(JZ,IZ)
 436              Z1(2)=0.
 437             END IF
 438             Z2(1)=APU
 439             Z2(2)=BTU
 440             CALL SSME3(1,AMZ,Z1,Z2)
 441             Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*AUI(IZ)/2.
 442             Z1(2)=-Z1(1)
 443             Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*AUI(JZ)/2.)
 444             Z2(2)=Z2(1)
 445             CALL SSME3(2,AMULSS,Z1,Z2)
 446             Z1(1)=ZI**(THIZ-1)*BUI(IZ)/2.
 447             Z1(2)=Z1(1)
 448             Z2(1)=CONJG(ZI**(THJZ-1)*BUI(JZ)/2.)
 449             Z2(2)=-Z2(1)
 450             CALL SSME3(2,AMURSS,Z1,Z2)
 451             Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*AUI(JZ)/2.
 452             Z1(2)=-Z1(1)
 453             Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*AUI(IZ)/2.)
 454             Z2(2)=Z2(1)
 455             CALL SSME3(3,AMULSS,Z1,Z2)
 456             Z1(1)=ZI**(THJZ-1)*BUI(JZ)/2.
 457             Z1(2)=Z1(1)
 458             Z2(1)=-CONJG(ZI**(THIZ-1)*BUI(IZ)/2.)
 459             Z2(2)=-Z2(1)
 460             CALL SSME3(3,AMURSS,Z1,Z2)
 461           END IF
 462 C          zi -> zj + d + dbar
 463           TMP(4)=-SNIJ(JZ,IZ)
 464           IF (MZIZ.GT.(MZJZ+2*AMDN)) THEN
 465             IF (MZIZ.LT.AMDLSS) THEN
 466             TMP(3)=AMDLSS
 467             TERMLL=3*2*ADI(JZ)**2*ADI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.)
 468             ELSE
 469             TERMLL=0.
 470             END IF
 471             IF (MZIZ.LT.AMDRSS) THEN
 472             TMP(3)=AMDRSS
 473             TERMRR=3*2*BDI(JZ)**2*BDI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.)
 474             ELSE
 475             TERMRR=0.
 476             END IF
 477             IF (MZIZ.LT.(MZJZ+AMZ)) THEN
 478             TMP(4)=+SNIJ(JZ,IZ)
 479             TERMZZ=3*(APD**2+BTD**2)*E**2*WIJ(JZ,IZ)**2/MZIZ
 480      $             *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ)
 481             ELSE
 482             TERMZZ=0.
 483             END IF
 484             IF (MZIZ.LT.AMDLSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
 485             TMP(4)=+SNIJ(JZ,IZ)
 486             TMP(3)=AMDLSS
 487             TERMLZ=3*8*E*(APD-BTD)*WIJ(JZ,IZ)*ADI(JZ)*ADI(IZ)/MZIZ
 488      $      /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2)
 489             ELSE
 490             TERMLZ=0.
 491             END IF
 492             IF (MZIZ.LT.AMDRSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
 493             TMP(4)=+SNIJ(JZ,IZ)
 494             TMP(3)=AMDRSS
 495             TERMRZ=-3*8*E*(APD+BTD)*WIJ(JZ,IZ)*BDI(JZ)*BDI(IZ)/MZIZ
 496      $      /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2)
 497             ELSE
 498             TERMRZ=0.
 499             END IF
 500             WID=TERMLL+TERMRR+TERMZZ+TERMLZ+TERMRZ
 501             CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDDN,-IDDN,0,0)
 502 C     Enter information for decay matrix element
 503             IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN
 504              Z1(1)=0.
 505              Z1(2)=2*E*WIJ(JZ,IZ)
 506             ELSE
 507              Z1(1)=2*E*WIJ(JZ,IZ)
 508              Z1(2)=0.
 509             END IF
 510             Z2(1)=APD
 511             Z2(2)=BTD
 512             CALL SSME3(1,AMZ,Z1,Z2)
 513             Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*ADI(IZ)/2.
 514             Z1(2)=-Z1(1)
 515             Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ADI(JZ)/2.)
 516             Z2(2)=Z2(1)
 517             CALL SSME3(2,AMDLSS,Z1,Z2)
 518             Z1(1)=ZI**(THIZ-1)*BDI(IZ)/2.
 519             Z1(2)=Z1(1)
 520             Z2(1)=CONJG(ZI**(THJZ-1)*BDI(JZ)/2.)
 521             Z2(2)=-Z2(1)
 522             CALL SSME3(2,AMDRSS,Z1,Z2)
 523             Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ADI(JZ)/2.
 524             Z1(2)=-Z1(1)
 525             Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*ADI(IZ)/2.)
 526             Z2(2)=Z2(1)
 527             CALL SSME3(3,AMDLSS,Z1,Z2)
 528             Z1(1)=ZI**(THJZ-1)*BDI(JZ)/2.
 529             Z1(2)=Z1(1)
 530             Z2(1)=-CONJG(ZI**(THIZ-1)*BDI(IZ)/2.)
 531             Z2(2)=-Z2(1)
 532             CALL SSME3(3,AMDRSS,Z1,Z2)
 533           END IF
 534 C          zi -> zj + s + sbar
 535           TMP(4)=-SNIJ(JZ,IZ)
 536           IF (MZIZ.GT.(MZJZ+2*AMST)) THEN
 537             IF (MZIZ.LT.AMSLSS) THEN
 538             TMP(3)=AMSLSS
 539             TERMLL=3*2*ADI(JZ)**2*ADI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.)
 540             ELSE
 541             TERMLL=0.
 542             END IF
 543             IF (MZIZ.LT.AMSRSS) THEN
 544             TMP(3)=AMSRSS
 545             TERMRR=3*2*BDI(JZ)**2*BDI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.)
 546             ELSE
 547             TERMRR=0.
 548             END IF
 549             IF (MZIZ.LT.(MZJZ+AMZ)) THEN
 550             TMP(4)=+SNIJ(JZ,IZ)
 551             TERMZZ=3*(APD**2+BTD**2)*E**2*WIJ(JZ,IZ)**2/MZIZ
 552      $             *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ)
 553             ELSE
 554             TERMZZ=0.
 555             END IF
 556             IF (MZIZ.LT.AMSLSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
 557             TMP(4)=+SNIJ(JZ,IZ)
 558             TMP(3)=AMSLSS
 559             TERMLZ=3*8*E*(APD-BTD)*WIJ(JZ,IZ)*ADI(JZ)*ADI(IZ)/MZIZ
 560      $      /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2)
 561             ELSE
 562             TERMLZ=0.
 563             END IF
 564             IF (MZIZ.LT.AMSRSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
 565             TMP(4)=+SNIJ(JZ,IZ)
 566             TMP(3)=AMSRSS
 567             TERMRZ=-3*8*E*(APD+BTD)*WIJ(JZ,IZ)*BDI(JZ)*BDI(IZ)/MZIZ
 568      $      /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2)
 569             ELSE
 570             TERMRZ=0.
 571             END IF
 572             WID=TERMLL+TERMRR+TERMZZ+TERMLZ+TERMRZ
 573             CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDST,-IDST,0,0)
 574 C     Enter information for decay matrix element
 575             IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN
 576              Z1(1)=0.
 577              Z1(2)=2*E*WIJ(JZ,IZ)
 578             ELSE
 579              Z1(1)=2*E*WIJ(JZ,IZ)
 580              Z1(2)=0.
 581             END IF
 582             Z2(1)=APD
 583             Z2(2)=BTD
 584             CALL SSME3(1,AMZ,Z1,Z2)
 585             Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*ADI(IZ)/2.
 586             Z1(2)=-Z1(1)
 587             Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ADI(JZ)/2.)
 588             Z2(2)=Z2(1)
 589             CALL SSME3(2,AMSLSS,Z1,Z2)
 590             Z1(1)=ZI**(THIZ-1)*BDI(IZ)/2.
 591             Z1(2)=Z1(1)
 592             Z2(1)=CONJG(ZI**(THJZ-1)*BDI(JZ)/2.)
 593             Z2(2)=-Z2(1)
 594             CALL SSME3(2,AMSRSS,Z1,Z2)
 595             Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ADI(JZ)/2.
 596             Z1(2)=-Z1(1)
 597             Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*ADI(IZ)/2.)
 598             Z2(2)=Z2(1)
 599             CALL SSME3(3,AMSLSS,Z1,Z2)
 600             Z1(1)=ZI**(THJZ-1)*BDI(JZ)/2.
 601             Z1(2)=Z1(1)
 602             Z2(1)=-CONJG(ZI**(THIZ-1)*BDI(IZ)/2.)
 603             Z2(2)=-Z2(1)
 604             CALL SSME3(3,AMSRSS,Z1,Z2)
 605           END IF
 606 C          zi -> zj + c + cbar
 607           TMP(4)=-SNIJ(JZ,IZ)
 608           IF (MZIZ.GT.(MZJZ+2*AMCH)) THEN
 609             IF (MZIZ.LT.AMCLSS) THEN
 610             TMP(3)=AMCLSS
 611             TERMLL=3*2*AUI(JZ)**2*AUI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.)
 612             ELSE
 613             TERMLL=0.
 614             END IF
 615             IF (MZIZ.LT.AMCRSS) THEN
 616             TMP(3)=AMCRSS
 617             TERMRR=3*2*BUI(JZ)**2*BUI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.)
 618             ELSE
 619             TERMRR=0.
 620             END IF
 621             IF (MZIZ.LT.(MZJZ+AMZ)) THEN
 622             TMP(4)=+SNIJ(JZ,IZ)
 623             TERMZZ=3*(APU**2+BTU**2)*E**2*WIJ(JZ,IZ)**2/MZIZ
 624      $             *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ)
 625             ELSE
 626             TERMZZ=0.
 627             END IF
 628             IF (MZIZ.LT.AMCLSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
 629             TMP(4)=+SNIJ(JZ,IZ)
 630             TMP(3)=AMCLSS
 631             TERMLZ=3*8*E*(APU-BTU)*WIJ(JZ,IZ)*AUI(JZ)*AUI(IZ)/MZIZ
 632      $      /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2)
 633             ELSE
 634             TERMLZ=0.
 635             END IF
 636             IF (MZIZ.LT.AMCRSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
 637             TMP(4)=+SNIJ(JZ,IZ)
 638             TMP(3)=AMCRSS
 639             TERMRZ=-3*8*E*(APU+BTU)*WIJ(JZ,IZ)*BUI(JZ)*BUI(IZ)/MZIZ
 640      $      /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2)
 641             ELSE
 642             TERMRZ=0.
 643             END IF
 644             WID=TERMLL+TERMRR+TERMZZ+TERMLZ+TERMRZ
 645             CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDCH,-IDCH,0,0)
 646 C     Enter information for decay matrix element
 647             IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN
 648              Z1(1)=0.
 649              Z1(2)=2*E*WIJ(JZ,IZ)
 650             ELSE
 651              Z1(1)=2*E*WIJ(JZ,IZ)
 652              Z1(2)=0.
 653             END IF
 654             Z2(1)=APU
 655             Z2(2)=BTU
 656             CALL SSME3(1,AMZ,Z1,Z2)
 657             Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*AUI(IZ)/2.
 658             Z1(2)=-Z1(1)
 659             Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*AUI(JZ)/2.)
 660             Z2(2)=Z2(1)
 661             CALL SSME3(2,AMCLSS,Z1,Z2)
 662             Z1(1)=ZI**(THIZ-1)*BUI(IZ)/2.
 663             Z1(2)=Z1(1)
 664             Z2(1)=CONJG(ZI**(THJZ-1)*BUI(JZ)/2.)
 665             Z2(2)=-Z2(1)
 666             CALL SSME3(2,AMCRSS,Z1,Z2)
 667             Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*AUI(JZ)/2.
 668             Z1(2)=-Z1(1)
 669             Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*AUI(IZ)/2.)
 670             Z2(2)=Z2(1)
 671             CALL SSME3(3,AMCLSS,Z1,Z2)
 672             Z1(1)=ZI**(THJZ-1)*BUI(JZ)/2.
 673             Z1(2)=Z1(1)
 674             Z2(1)=-CONJG(ZI**(THIZ-1)*BUI(IZ)/2.)
 675             Z2(2)=-Z2(1)
 676             CALL SSME3(3,AMCRSS,Z1,Z2)
 677           END IF
 678 C          zi -> zj + b+ bbar ; mixing/yukawa effects now included!
 679 C                              thanks to M. Drees
 680           FACT=1./2./(2*PI)**5/2./MZIZ
 681           ALIZ1=ADI(IZ)*COSB-FB*ZMIXSS(2,IZ)*SINB
 682           ALJZ1=ADI(JZ)*COSB-FB*ZMIXSS(2,JZ)*SINB
 683           ALIZ2=ADI(IZ)*SINB+FB*ZMIXSS(2,IZ)*COSB
 684           ALJZ2=ADI(JZ)*SINB+FB*ZMIXSS(2,JZ)*COSB
 685           BEIZ1=BDI(IZ)*SINB+FB*ZMIXSS(2,IZ)*COSB
 686           BEJZ1=BDI(JZ)*SINB+FB*ZMIXSS(2,JZ)*COSB
 687           BEIZ2=-BDI(IZ)*COSB+FB*ZMIXSS(2,IZ)*SINB
 688           BEJZ2=-BDI(JZ)*COSB+FB*ZMIXSS(2,JZ)*SINB
 689           SGNIJ=-SNIJ(JZ,IZ)
 690           XUPPER=(MZIZ**2+AMBT**2-(AMBT+MZJZ)**2)/2./MZIZ
 691           IF (MZIZ.GT.(MZJZ+2*AMBT)) THEN
 692             TMP(1)=MZIZ
 693             TMP(2)=AMBT
 694             TMP(3)=MZJZ
 695             TMP(4)=AMB1SS
 696             TMP(5)=AMB1SS
 697             TMP(6)=AMBT
 698             IF (MZIZ.LT.AMB1SS) THEN
 699               GLLF1=4*ALIZ1**2*((ALJZ1**2+BEJZ1**2)*
 700      ,             SSXINT(AMBT,SSGX1,XUPPER)+SGNIJ*ALJZ1**2*
 701      ,             SSXINT(AMBT,SSGX2,XUPPER))
 702               GRRF1=4*BEIZ1**2*((ALJZ1**2+BEJZ1**2)*
 703      ,             SSXINT(AMBT,SSGX1,XUPPER)+SGNIJ*BEJZ1**2*
 704      ,             SSXINT(AMBT,SSGX2,XUPPER))
 705               GLRF1=-8*ALIZ1*BEIZ1*ALJZ1*BEJZ1*
 706      ,               SSXINT(AMBT,SSGX8,XUPPER)
 707               GF1=GLLF1+GRRF1+GLRF1
 708             ELSE
 709               GF1=0.
 710             END IF
 711             IF (MZIZ.LT.AMB1SS) THEN
 712             TMP(4)=AMB1SS
 713             TMP(5)=AMB2SS
 714               GLGL=8*ALIZ1*ALIZ2*(ALJZ1*ALJZ2+BEJZ1*BEJZ2)*
 715      ,             SSXINT(AMBT,SSGX1,XUPPER)+SGNIJ*8*ALIZ1*ALIZ2*
 716      ,             ALJZ1*ALJZ2*SSXINT(AMBT,SSGX2,XUPPER)
 717               GRGR=8*BEIZ1*BEIZ2*(ALJZ1*ALJZ2+BEJZ1*BEJZ2)*
 718      ,             SSXINT(AMBT,SSGX1,XUPPER)+SGNIJ*8*BEIZ1*BEIZ2*
 719      ,             BEJZ1*BEJZ2*SSXINT(AMBT,SSGX2,XUPPER)
 720               GLGR=-8*ALIZ1*BEIZ2*ALJZ2*BEJZ1*
 721      ,              SSXINT(AMBT,SSGX8,XUPPER)
 722               GRGL=-8*ALIZ2*BEIZ1*ALJZ1*BEJZ2*
 723      ,              SSXINT(AMBT,SSGX8,XUPPER)
 724               GF12=GLGL+GRGR+GLGR+GRGL
 725             ELSE
 726               GF12=0.
 727             END IF
 728             IF (MZIZ.LT.AMB2SS) THEN
 729             TMP(4)=AMB2SS
 730             TMP(5)=AMB2SS
 731               GLLF2=4*ALIZ2**2*((ALJZ2**2+BEJZ2**2)*
 732      ,              SSXINT(AMBT,SSGX1,XUPPER)+SGNIJ*ALJZ2**2*
 733      ,              SSXINT(AMBT,SSGX2,XUPPER))
 734               GRRF2=4*BEIZ2**2*((ALJZ2**2+BEJZ2**2)*
 735      ,              SSXINT(AMBT,SSGX1,XUPPER)+SGNIJ*BEJZ2**2*
 736      ,              SSXINT(AMBT,SSGX2,XUPPER))
 737               GLRF2=-8*ALIZ2*BEIZ2*ALJZ2*BEJZ2*
 738      ,              SSXINT(AMBT,SSGX8,XUPPER)
 739               GF2=GLLF2+GRRF2+GLRF2
 740             ELSE
 741               GF2=0.
 742             END IF
 743             GF=FACT*(GF1+GF2+GF12)
 744             TMP(2)=MZJZ
 745             IF (MZIZ.LT.(MZJZ+AMZ)) THEN
 746             TMP(4)=+SNIJ(JZ,IZ)
 747             TERMZZ=(APD**2+BTD**2)*E**2*WIJ(JZ,IZ)**2/MZIZ
 748      $       *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2-4*AMBT**2)/2./MZIZ)
 749             ELSE
 750             TERMZZ=0.
 751             END IF
 752             IF (MZIZ.LE.(MZJZ+AMHL)) THEN
 753             TMP(3)=AMHL
 754             TMP(5)=AMHL
 755             TMP(4)=-SNIJ(JZ,IZ)
 756             TERMHL=G**2/64./PI**3/MZIZ*(MBQ*SINA*(XLIJ(JZ,IZ)+
 757      $        XLIJ(IZ,JZ))/AMW/COSBE)**2*
 758      $        SSXINT(MZJZ,SSZZF4,(MZIZ**2+MZJZ**2-4*AMBT**2)/2./MZIZ)
 759             ELSE
 760             TERMHL=0.
 761             END IF
 762             IF (MZIZ.LE.(MZJZ+AMHH)) THEN
 763             TMP(3)=AMHH
 764             TMP(5)=AMHH
 765             TMP(4)=-SNIJ(JZ,IZ)
 766             TERMHH=G**2/64./PI**3/MZIZ*(MBQ*COSA*(XHIJ(JZ,IZ)+
 767      $       XHIJ(IZ,JZ))/AMW/COSBE)**2*
 768      $       SSXINT(MZJZ,SSZZF4,(MZIZ**2+MZJZ**2-4*AMBT**2)/2./MZIZ)
 769             ELSE
 770             TERMHH=0.
 771             END IF
 772             IF (MZIZ.LE.(MZJZ+AMHH).AND.MZIZ.LE.(MZJZ+AMHL)) THEN
 773             TMP(3)=AMHL
 774             TMP(5)=AMHH
 775             TMP(4)=-SNIJ(JZ,IZ)
 776             TERMLH=2*G**2/64./PI**3/MZIZ*(MBQ/AMW/COSBE)**2*
 777      $            (COSA*(XHIJ(JZ,IZ)+XHIJ(IZ,JZ))*SINA*
 778      $       (XLIJ(JZ,IZ)+XLIJ(IZ,JZ)))*
 779      $       SSXINT(MZJZ,SSZZF4,(MZIZ**2+MZJZ**2-4*AMBT**2)/2./MZIZ)
 780             ELSE
 781             TERMLH=0.
 782             END IF
 783             IF (MZIZ.LE.(MZJZ+AMHA)) THEN
 784             TMP(3)=AMHA
 785             TMP(5)=AMHA
 786             TMP(4)=SNIJ(JZ,IZ)
 787             TERMHA=G**2*TANB**2/64./PI**3/MZIZ*(MBQ*(XPIJ(JZ,IZ)+
 788      $            XPIJ(IZ,JZ))/AMW)**2*
 789      $       SSXINT(MZJZ,SSZZF4,(MZIZ**2+MZJZ**2-4*AMBT**2)/2./MZIZ)
 790             ELSE
 791             TERMHA=0.
 792             END IF
 793             IF (MZIZ.LT.AMB1SS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
 794             TMP(4)=-SGNIJ
 795             TMP(3)=AMB1SS
 796             TERM1Z=8*E*WIJ(JZ,IZ)*(ALJZ1*ALIZ1*(APD-BTD)-
 797      $       BEJZ1*BEIZ1*(APD+BTD))/MZIZ
 798      $      /(2*PI)**5*SSXINT(4*AMBT**2,SSZZF3,(MZIZ-MZJZ)**2)
 799             ELSE
 800             TERM1Z=0.
 801             END IF
 802             IF (MZIZ.LT.AMB2SS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
 803             TMP(4)=-SGNIJ
 804             TMP(3)=AMB2SS
 805             TERM2Z=8*E*WIJ(JZ,IZ)*(ALJZ2*ALIZ2*(APD-BTD)-
 806      $       BEJZ2*BEIZ2*(APD+BTD))/MZIZ
 807      $      /(2*PI)**5*SSXINT(4*AMBT**2,SSZZF3,(MZIZ-MZJZ)**2)
 808             ELSE
 809             TERM2Z=0.
 810             END IF
 811             IF (MZIZ.LT.AMB1SS.AND.MZIZ.LT.(MZJZ+AMHL)) THEN
 812             TMP(3)=AMB1SS
 813             TMP(4)=SGNIJ
 814             TMP(5)=AMHL
 815             TERM1L=2*PI*PI*G*MBQ*SINA*SGNIJ*(XLIJ(IZ,JZ)+
 816      $        XLIJ(JZ,IZ))*(ALIZ1*BEJZ1+ALJZ1*BEIZ1)
 817      $       /MZIZ/AMW/COSBE*SSXINT(4*AMBT**2,SSZZF5,(MZIZ-MZJZ)**2)
 818             TERM1L=FACT*TERM1L
 819             ELSE
 820             TERM1L=0.
 821             END IF
 822             IF (MZIZ.LT.AMB2SS.AND.MZIZ.LT.(MZJZ+AMHL)) THEN
 823             TMP(3)=AMB2SS
 824             TMP(4)=SGNIJ
 825             TMP(5)=AMHL
 826             TERM2L=2*PI*PI*G*MBQ*SINA*SGNIJ*(XLIJ(IZ,JZ)+
 827      $        XLIJ(JZ,IZ))*(ALIZ2*BEJZ2+ALJZ2*BEIZ2)
 828      $       /MZIZ/AMW/COSBE*SSXINT(4*AMBT**2,SSZZF5,(MZIZ-MZJZ)**2)
 829             TERM2L=FACT*TERM2L
 830             ELSE
 831             TERM2L=0.
 832             END IF
 833             IF (MZIZ.LT.AMB1SS.AND.MZIZ.LT.(MZJZ+AMHH)) THEN
 834             TMP(3)=AMB1SS
 835             TMP(4)=SGNIJ
 836             TMP(5)=AMHH
 837             TERM1H=2*PI*PI*G*MBQ*COSA*SGNIJ*(XHIJ(IZ,JZ)+
 838      $        XHIJ(JZ,IZ))*(ALIZ1*BEJZ1+ALJZ1*BEIZ1)
 839      $       /MZIZ/AMW/COSBE*SSXINT(4*AMBT**2,SSZZF5,(MZIZ-MZJZ)**2)
 840             TERM1H=FACT*TERM1H
 841             ELSE
 842             TERM1H=0.
 843             END IF
 844             IF (MZIZ.LT.AMB2SS.AND.MZIZ.LT.(MZJZ+AMHH)) THEN
 845             TMP(3)=AMB2SS
 846             TMP(4)=SGNIJ
 847             TMP(5)=AMHH
 848             TERM2H=2*PI*PI*G*MBQ*COSA*SGNIJ*(XHIJ(IZ,JZ)+
 849      $        XHIJ(JZ,IZ))*(ALIZ2*BEJZ2+ALJZ2*BEIZ2)
 850      $       /MZIZ/AMW/COSBE*SSXINT(4*AMBT**2,SSZZF5,(MZIZ-MZJZ)**2)
 851             TERM2H=FACT*TERM2H
 852             ELSE
 853             TERM2H=0.
 854             END IF
 855             IF (MZIZ.LT.AMB1SS.AND.MZIZ.LT.(MZJZ+AMHA)) THEN
 856             TMP(3)=AMB1SS
 857             TMP(4)=-SGNIJ
 858             TMP(5)=AMHA
 859             TERM1A=-2*PI*PI*G*MBQ*TANB*SGNIJ*(XPIJ(IZ,JZ)+
 860      $        XPIJ(JZ,IZ))*(ALIZ1*BEJZ1+ALJZ1*BEIZ1)
 861      $        /MZIZ/AMW*SSXINT(4*AMBT**2,SSZZF5,(MZIZ-MZJZ)**2)
 862             TERM1A=FACT*TERM1A
 863             ELSE
 864             TERM1A=0.
 865             END IF
 866             IF (MZIZ.LT.AMB2SS.AND.MZIZ.LT.(MZJZ+AMHA)) THEN
 867             TMP(3)=AMB2SS
 868             TMP(4)=-SGNIJ
 869             TMP(5)=AMHA
 870             TERM2A=-2*PI*PI*G*MBQ*TANB*SGNIJ*(XPIJ(IZ,JZ)+
 871      $        XPIJ(JZ,IZ))*(ALIZ2*BEJZ2+ALJZ2*BEIZ2)
 872      $        /MZIZ/AMW*SSXINT(4*AMBT**2,SSZZF5,(MZIZ-MZJZ)**2)
 873             TERM2A=FACT*TERM2A
 874             ELSE
 875             TERM2A=0.
 876             END IF
 877             WID=3*(GF+TERMZZ+TERMHL+TERMHH+TERMLH+TERMHA+TERM1Z+
 878      $        TERM2Z+TERM1L+TERM2L+TERM1H+TERM2H+TERM1A+TERM2A)
 879             CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDBT,-IDBT,0,0)
 880 C     Enter information for decay matrix element
 881             IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN
 882              Z1(1)=0.
 883              Z1(2)=2*E*WIJ(JZ,IZ)
 884             ELSE
 885              Z1(1)=2*E*WIJ(JZ,IZ)
 886              Z1(2)=0.
 887             END IF
 888             Z2(1)=APD
 889             Z2(2)=BTD
 890             CALL SSME3(1,AMZ,Z1,Z2)
 891             Z1(1)=((ZI**(THIZ)*(-1.)**(THIZ+1)*ADI(IZ)-
 892      $            FB*ZMIXSS(2,IZ)*ZI**THIZ)*COSB-(ZI**THIZ*BDI(IZ)-
 893      $            FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINB)/2.
 894             Z1(2)=((-ZI**(THIZ)*(-1.)**(THIZ+1)*ADI(IZ)-
 895      $            FB*ZMIXSS(2,IZ)*ZI**THIZ)*COSB-(ZI**THIZ*BDI(IZ)+
 896      $            FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINB)/2.
 897             Z2(1)=CONJG(((ZI**(THJZ)*(-1.)**(THJZ+1)*ADI(JZ)-
 898      $            FB*ZMIXSS(2,JZ)*ZI**THJZ)*COSB-(ZI**THJZ*BDI(JZ)-
 899      $            FB*ZMIXSS(2,JZ)*(-ZI)**THJZ)*SINB)/2.)
 900             Z2(2)=-CONJG(((-ZI**(THJZ)*(-1.)**(THJZ+1)*ADI(JZ)-
 901      $            FB*ZMIXSS(2,JZ)*ZI**THJZ)*COSB-(ZI**THJZ*BDI(JZ)+
 902      $            FB*ZMIXSS(2,JZ)*(-ZI)**THJZ)*SINB)/2.)
 903             CALL SSME3(2,AMB1SS,Z1,Z2)
 904             Z1(1)=((ZI**(THIZ)*(-1.)**(THIZ+1)*ADI(IZ)-
 905      $            FB*ZMIXSS(2,IZ)*ZI**THIZ)*SINB+(ZI**THIZ*BDI(IZ)-
 906      $            FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSB)/2.
 907             Z1(2)=((-ZI**(THIZ)*(-1.)**(THIZ+1)*ADI(IZ)-
 908      $            FB*ZMIXSS(2,IZ)*ZI**THIZ)*SINB+(ZI**THIZ*BDI(IZ)+
 909      $            FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSB)/2.
 910             Z2(1)=CONJG(((ZI**(THJZ)*(-1.)**(THJZ+1)*ADI(JZ)-
 911      $            FB*ZMIXSS(2,JZ)*ZI**THJZ)*SINB+(ZI**THJZ*BDI(JZ)-
 912      $            FB*ZMIXSS(2,JZ)*(-ZI)**THJZ)*COSB)/2.)
 913             Z2(2)=-CONJG(((-ZI**(THJZ)*(-1.)**(THJZ+1)*ADI(JZ)-
 914      $            FB*ZMIXSS(2,JZ)*ZI**THJZ)*SINB+(ZI**THJZ*BDI(JZ)+
 915      $            FB*ZMIXSS(2,JZ)*(-ZI)**THJZ)*COSB)/2.)
 916             CALL SSME3(2,AMB2SS,Z1,Z2)
 917             Z1(1)=((ZI**(THJZ)*(-1.)**(THJZ+1)*ADI(JZ)-
 918      $            FB*ZMIXSS(2,JZ)*ZI**THJZ)*COSB-(ZI**THJZ*BDI(JZ)-
 919      $            FB*ZMIXSS(2,JZ)*(-ZI)**THJZ)*SINB)/2.
 920             Z1(2)=((-ZI**(THJZ)*(-1.)**(THJZ+1)*ADI(JZ)-
 921      $            FB*ZMIXSS(2,JZ)*ZI**THJZ)*COSB-(ZI**THJZ*BDI(JZ)+
 922      $            FB*ZMIXSS(2,JZ)*(-ZI)**THJZ)*SINB)/2.
 923             Z2(1)=-CONJG(((ZI**(THIZ)*(-1.)**(THIZ+1)*ADI(IZ)-
 924      $            FB*ZMIXSS(2,IZ)*ZI**THIZ)*COSB-(ZI**THIZ*BDI(IZ)-
 925      $            FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINB)/2.)
 926             Z2(2)=CONJG(((-ZI**(THIZ)*(-1.)**(THIZ+1)*ADI(IZ)-
 927      $            FB*ZMIXSS(2,IZ)*ZI**THIZ)*COSB-(ZI**THIZ*BDI(IZ)+
 928      $            FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINB)/2.)
 929             CALL SSME3(3,AMB1SS,Z1,Z2)
 930             Z1(1)=((ZI**(THJZ)*(-1.)**(THJZ+1)*ADI(JZ)-
 931      $            FB*ZMIXSS(2,JZ)*ZI**THJZ)*SINB+(ZI**THJZ*BDI(JZ)-
 932      $            FB*ZMIXSS(2,JZ)*(-ZI)**THJZ)*COSB)/2.
 933             Z1(2)=((-ZI**(THJZ)*(-1.)**(THJZ+1)*ADI(JZ)-
 934      $            FB*ZMIXSS(2,JZ)*ZI**THJZ)*SINB+(ZI**THJZ*BDI(JZ)+
 935      $            FB*ZMIXSS(2,JZ)*(-ZI)**THJZ)*COSB)/2.
 936             Z2(1)=-CONJG(((ZI**(THIZ)*(-1.)**(THIZ+1)*ADI(IZ)-
 937      $            FB*ZMIXSS(2,IZ)*ZI**THIZ)*SINB+(ZI**THIZ*BDI(IZ)-
 938      $            FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSB)/2.)
 939             Z2(2)=CONJG(((-ZI**(THIZ)*(-1.)**(THIZ+1)*ADI(IZ)-
 940      $            FB*ZMIXSS(2,IZ)*ZI**THIZ)*SINB+(ZI**THIZ*BDI(IZ)+
 941      $            FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSB)/2.)
 942             CALL SSME3(3,AMB2SS,Z1,Z2)
 943             IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN
 944              Z1(1)=(-ZI)**(THIZ+THJZ)*(XLIJ(IZ,JZ)+XLIJ(JZ,IZ))
 945              Z1(2)=0.
 946             ELSE
 947              Z1(1)=0.
 948              Z1(2)=(-ZI)**(THIZ+THJZ)*(XLIJ(IZ,JZ)+XLIJ(JZ,IZ))
 949             END IF
 950             Z2(1)=-G*MBQ*SINA/2./AMW/COSBE
 951             Z2(2)=0.
 952             CALL SSME3(4,AMHL,Z1,Z2)
 953             IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN
 954              Z1(1)=(-ZI)**(THIZ+THJZ)*(XHIJ(IZ,JZ)+XHIJ(JZ,IZ))
 955              Z1(2)=0.
 956             ELSE
 957              Z1(1)=0.
 958              Z1(2)=(-ZI)**(THIZ+THJZ)*(XHIJ(IZ,JZ)+XHIJ(JZ,IZ))
 959             END IF
 960             Z2(1)=-G*MBQ*COSA/2./AMW/COSBE
 961             Z2(2)=0.
 962             CALL SSME3(4,AMHH,Z1,Z2)
 963             IF (-SNIJ(JZ,IZ).GT.0.) THEN
 964              Z1(1)=0.
 965              Z1(2)=(-ZI)**(THIZ+THJZ+1)*(XPIJ(IZ,JZ)+XPIJ(JZ,IZ))
 966             ELSE
 967              Z1(1)=(-ZI)**(THIZ+THJZ+1)*(XPIJ(IZ,JZ)+XPIJ(JZ,IZ))
 968              Z1(2)=0.
 969             END IF
 970             Z2(1)=0.
 971             Z2(2)=ZI*G*MBQ*TANB/2./AMW
 972             CALL SSME3(4,AMHA,Z1,Z2)
 973           END IF
 974 C          zi -> zj + e + ebar
 975           TMP(2)=MZJZ
 976           TMP(4)=-SNIJ(JZ,IZ)
 977           TMP(6)=0.
 978           IF (MZIZ.GT.(MZJZ+2*AME)) THEN
 979             IF (MZIZ.LT.AMELSS) THEN
 980             TMP(3)=AMELSS
 981             TERMLL=2*ALI(JZ)**2*ALI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.)
 982             ELSE
 983             TERMLL=0.
 984             END IF
 985             IF (MZIZ.LT.AMERSS) THEN
 986             TMP(3)=AMERSS
 987             TERMRR=2*BLI(JZ)**2*BLI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.)
 988             ELSE
 989             TERMRR=0.
 990             END IF
 991             IF (MZIZ.LT.(MZJZ+AMZ)) THEN
 992             TMP(4)=+SNIJ(JZ,IZ)
 993             TERMZZ=(APL**2+BTL**2)*E**2*WIJ(JZ,IZ)**2/MZIZ
 994      $             *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ)
 995             ELSE
 996             TERMZZ=0.
 997             END IF
 998             IF (MZIZ.LT.AMELSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
 999             TMP(4)=+SNIJ(JZ,IZ)
1000             TMP(3)=AMELSS
1001             TERMLZ=8*E*(APL-BTL)*WIJ(JZ,IZ)*ALI(JZ)*ALI(IZ)/MZIZ
1002      $      /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2)
1003             ELSE
1004             TERMLZ=0.
1005             END IF
1006             IF (MZIZ.LT.AMERSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
1007             TMP(4)=+SNIJ(JZ,IZ)
1008             TMP(3)=AMERSS
1009             TERMRZ=-8*E*(APL+BTL)*WIJ(JZ,IZ)*BLI(JZ)*BLI(IZ)/MZIZ
1010      $      /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2)
1011             ELSE
1012             TERMRZ=0.
1013             END IF
1014             WID=TERMLL+TERMRR+TERMZZ+TERMLZ+TERMRZ
1015             CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDE,-IDE,0,0)
1016 C     Enter information for decay matrix element
1017             IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN
1018              Z1(1)=0.
1019              Z1(2)=2*E*WIJ(JZ,IZ)
1020             ELSE
1021              Z1(1)=2*E*WIJ(JZ,IZ)
1022              Z1(2)=0.
1023             END IF
1024             Z2(1)=APL
1025             Z2(2)=BTL
1026             CALL SSME3(1,AMZ,Z1,Z2)
1027             Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*ALI(IZ)/2.
1028             Z1(2)=-Z1(1)
1029             Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ALI(JZ)/2.)
1030             Z2(2)=Z2(1)
1031             CALL SSME3(2,AMELSS,Z1,Z2)
1032             Z1(1)=ZI**(THIZ-1)*BLI(IZ)/2.
1033             Z1(2)=Z1(1)
1034             Z2(1)=CONJG(ZI**(THJZ-1)*BLI(JZ)/2.)
1035             Z2(2)=-Z2(1)
1036             CALL SSME3(2,AMERSS,Z1,Z2)
1037             Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ALI(JZ)/2.
1038             Z1(2)=-Z1(1)
1039             Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*ALI(IZ)/2.)
1040             Z2(2)=Z2(1)
1041             CALL SSME3(3,AMELSS,Z1,Z2)
1042             Z1(1)=ZI**(THJZ-1)*BLI(JZ)/2.
1043             Z1(2)=Z1(1)
1044             Z2(1)=-CONJG(ZI**(THIZ-1)*BLI(IZ)/2.)
1045             Z2(2)=-Z2(1)
1046             CALL SSME3(3,AMERSS,Z1,Z2)
1047           END IF
1048 C          zi -> zj + mu + mubar
1049           TMP(4)=-SNIJ(JZ,IZ)
1050           IF (MZIZ.GT.(MZJZ+2*AMMU)) THEN
1051             IF (MZIZ.LT.AMMLSS) THEN
1052             TMP(3)=AMMLSS
1053             TERMLL=2*ALI(JZ)**2*ALI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.)
1054             ELSE
1055             TERMLL=0.
1056             END IF
1057             IF (MZIZ.LT.AMMRSS) THEN
1058             TMP(3)=AMMRSS
1059             TERMRR=2*BLI(JZ)**2*BLI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.)
1060             ELSE
1061             TERMRR=0.
1062             END IF
1063             IF (MZIZ.LT.(MZJZ+AMZ)) THEN
1064             TMP(4)=+SNIJ(JZ,IZ)
1065             TERMZZ=(APL**2+BTL**2)*E**2*WIJ(JZ,IZ)**2/MZIZ
1066      $             *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ)
1067             ELSE
1068             TERMZZ=0.
1069             END IF
1070             IF (MZIZ.LT.AMMLSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
1071             TMP(4)=+SNIJ(JZ,IZ)
1072             TMP(3)=AMMLSS
1073             TERMLZ=8*E*(APL-BTL)*WIJ(JZ,IZ)*ALI(JZ)*ALI(IZ)/MZIZ
1074      $      /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2)
1075             ELSE
1076             TERMLZ=0.
1077             END IF
1078             IF (MZIZ.LT.AMMRSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
1079             TMP(4)=+SNIJ(JZ,IZ)
1080             TMP(3)=AMMRSS
1081             TERMRZ=-8*E*(APL+BTL)*WIJ(JZ,IZ)*BLI(JZ)*BLI(IZ)/MZIZ
1082      $      /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2)
1083             ELSE
1084             TERMRZ=0.
1085             END IF
1086             WID=TERMLL+TERMRR+TERMZZ+TERMLZ+TERMRZ
1087             CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDMU,-IDMU,0,0)
1088 C     Enter information for decay matrix element
1089             IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN
1090              Z1(1)=0.
1091              Z1(2)=2*E*WIJ(JZ,IZ)
1092             ELSE
1093              Z1(1)=2*E*WIJ(JZ,IZ)
1094              Z1(2)=0.
1095             END IF
1096             Z2(1)=APL
1097             Z2(2)=BTL
1098             CALL SSME3(1,AMZ,Z1,Z2)
1099             Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*ALI(IZ)/2.
1100             Z1(2)=-Z1(1)
1101             Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ALI(JZ)/2.)
1102             Z2(2)=Z2(1)
1103             CALL SSME3(2,AMMLSS,Z1,Z2)
1104             Z1(1)=ZI**(THIZ-1)*BLI(IZ)/2.
1105             Z1(2)=Z1(1)
1106             Z2(1)=CONJG(ZI**(THJZ-1)*BLI(JZ)/2.)
1107             Z2(2)=-Z2(1)
1108             CALL SSME3(2,AMMRSS,Z1,Z2)
1109             Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ALI(JZ)/2.
1110             Z1(2)=-Z1(1)
1111             Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*ALI(IZ)/2.)
1112             Z2(2)=Z2(1)
1113             CALL SSME3(3,AMMLSS,Z1,Z2)
1114             Z1(1)=ZI**(THJZ-1)*BLI(JZ)/2.
1115             Z1(2)=Z1(1)
1116             Z2(1)=-CONJG(ZI**(THIZ-1)*BLI(IZ)/2.)
1117             Z2(2)=-Z2(1)
1118             CALL SSME3(3,AMMRSS,Z1,Z2)
1119           END IF
1120 C          zi -> zj + tau + taubar.
1121 C          Mixing/yukawa effects now included thanks to M. Drees
1122           ALIZ1=-ALI(IZ)*COSL-FL*ZMIXSS(2,IZ)*SINL
1123           ALJZ1=-ALI(JZ)*COSL-FL*ZMIXSS(2,JZ)*SINL
1124           ALIZ2=-ALI(IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
1125           ALJZ2=-ALI(JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
1126           BEIZ1=BLI(IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
1127           BEJZ1=BLI(JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
1128           BEIZ2=-BLI(IZ)*COSL+FL*ZMIXSS(2,IZ)*SINL
1129           BEJZ2=-BLI(JZ)*COSL+FL*ZMIXSS(2,JZ)*SINL
1130           SGNIJ=-SNIJ(JZ,IZ)
1131           XUPPER=(MZIZ**2+AMTAU**2-(AMTAU+MZJZ)**2)/2./MZIZ
1132 C          Polarization for stau_i -> z2ss+tau, z3ss+tau, z4ss+tau.
1133 C          See below for z1ss+tau.
1134           IF(JZ.EQ.1) THEN
1135             PTAU1(IZ)=(BEIZ1**2-ALIZ1**2)/(BEIZ1**2+ALIZ1**2)
1136             PTAU2(IZ)=(BEIZ2**2-ALIZ2**2)/(BEIZ2**2+ALIZ2**2)
1137           ENDIF
1138           IF (MZIZ.GT.(MZJZ+2*AMTAU)) THEN
1139             TMP(1)=MZIZ
1140             TMP(2)=AMTAU
1141             TMP(3)=MZJZ
1142             TMP(4)=AML1SS
1143             TMP(5)=AML1SS
1144             TMP(6)=AMTAU
1145             IF (MZIZ.LT.AML1SS) THEN
1146               GLLF1=4*ALIZ1**2*((ALJZ1**2+BEJZ1**2)*
1147      ,             SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*ALJZ1**2*
1148      ,             SSXINT(AMTAU,SSGX2,XUPPER))
1149               GLLF1L=FACT*4*ALIZ1**2*(ALJZ1**2*
1150      ,             SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*ALJZ1**2*
1151      ,             SSXINT(AMTAU,SSGX2,XUPPER))
1152               GRRF1=4*BEIZ1**2*((ALJZ1**2+BEJZ1**2)*
1153      ,             SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*BEJZ1**2*
1154      ,             SSXINT(AMTAU,SSGX2,XUPPER))
1155               GRRF1R=FACT*4*BEIZ1**2*(BEJZ1**2*
1156      ,             SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*BEJZ1**2*
1157      ,             SSXINT(AMTAU,SSGX2,XUPPER))
1158               GLRF1=-8*ALIZ1*BEIZ1*ALJZ1*BEJZ1*
1159      ,               SSXINT(AMTAU,SSGX8,XUPPER)
1160               GF1=GLLF1+GRRF1+GLRF1
1161             ELSE
1162               GLLF1L=0.
1163               GRRF1R=0.
1164               GF1=0.
1165             END IF
1166             IF (MZIZ.LT.AML1SS) THEN
1167             TMP(4)=AML1SS
1168             TMP(5)=AML2SS
1169               GLGL=8*ALIZ1*ALIZ2*(ALJZ1*ALJZ2+BEJZ1*BEJZ2)*
1170      ,             SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*8*ALIZ1*ALIZ2*
1171      ,             ALJZ1*ALJZ2*SSXINT(AMTAU,SSGX2,XUPPER)
1172               GLGLL=FACT*(8*ALIZ1*ALIZ2*ALJZ1*ALJZ2*
1173      ,             SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*8*ALIZ1*ALIZ2*
1174      ,             ALJZ1*ALJZ2*SSXINT(AMTAU,SSGX2,XUPPER))
1175               GRGR=8*BEIZ1*BEIZ2*(ALJZ1*ALJZ2+BEJZ1*BEJZ2)*
1176      ,             SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*8*BEIZ1*BEIZ2*
1177      ,             BEJZ1*BEJZ2*SSXINT(AMTAU,SSGX2,XUPPER)
1178               GRGRR=FACT*(8*BEIZ1*BEIZ2*BEJZ1*BEJZ2*
1179      ,             SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*8*BEIZ1*BEIZ2*
1180      ,             BEJZ1*BEJZ2*SSXINT(AMTAU,SSGX2,XUPPER))
1181               GLGR=-8*ALIZ1*BEIZ2*ALJZ2*BEJZ1*
1182      ,              SSXINT(AMTAU,SSGX8,XUPPER)
1183               GRGL=-8*ALIZ2*BEIZ1*ALJZ1*BEJZ2*
1184      ,              SSXINT(AMTAU,SSGX8,XUPPER)
1185               GF12=GLGL+GRGR+GLGR+GRGL
1186             ELSE
1187               GLGLL=0.
1188               GRGRR=0.
1189               GF12=0.
1190             END IF
1191             IF (MZIZ.LT.AML2SS) THEN
1192             TMP(4)=AML2SS
1193             TMP(5)=AML2SS
1194               GLLF2=4*ALIZ2**2*((ALJZ2**2+BEJZ2**2)*
1195      ,              SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*ALJZ2**2*
1196      ,              SSXINT(AMTAU,SSGX2,XUPPER))
1197               GLLF2L=FACT*4*ALIZ2**2*(ALJZ2**2*
1198      ,              SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*ALJZ2**2*
1199      ,              SSXINT(AMTAU,SSGX2,XUPPER))
1200               GRRF2=4*BEIZ2**2*((ALJZ2**2+BEJZ2**2)*
1201      ,              SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*BEJZ2**2*
1202      ,              SSXINT(AMTAU,SSGX2,XUPPER))
1203               GRRF2R=FACT*4*BEIZ2**2*(BEJZ2**2*
1204      ,              SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*BEJZ2**2*
1205      ,              SSXINT(AMTAU,SSGX2,XUPPER))
1206               GLRF2=-8*ALIZ2*BEIZ2*ALJZ2*BEJZ2*
1207      ,              SSXINT(AMTAU,SSGX8,XUPPER)
1208               GF2=GLLF2+GRRF2+GLRF2
1209             ELSE
1210               GLLF2L=0.
1211               GRRF2R=0.
1212               GF2=0.
1213             END IF
1214             GF=FACT*(GF1+GF2+GF12)
1215             TMP(2)=MZJZ
1216             IF (MZIZ.LT.(MZJZ+AMZ)) THEN
1217             TMP(4)=+SNIJ(JZ,IZ)
1218             TERMZZ=(APL**2+BTL**2)*E**2*WIJ(JZ,IZ)**2/MZIZ
1219      $      *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2-4*AMTAU**2)/2./MZIZ)
1220             TMZZRL=TERMZZ*(APL+BTL)**2/2./(APL**2+BTL**2)
1221             TMZZLR=TERMZZ*(APL-BTL)**2/2./(APL**2+BTL**2)
1222             ELSE
1223             TERMZZ=0.
1224             TMZZRL=0.
1225             TMZZLR=0.
1226             END IF
1227             IF (MZIZ.LE.(MZJZ+AMHL)) THEN
1228             TMP(3)=AMHL
1229             TMP(5)=AMHL
1230             TMP(4)=-SNIJ(JZ,IZ)
1231             TERMHL=G**2/64./PI**3/MZIZ*(MTAMZ*SINA*(XLIJ(JZ,IZ)+
1232      $         XLIJ(IZ,JZ))/AMW/COSBE)**2*
1233      $      SSXINT(MZJZ,SSZZF4,(MZIZ**2+MZJZ**2-4*AMTAU**2)/2./MZIZ)
1234             ELSE
1235             TERMHL=0.
1236             END IF
1237             IF (MZIZ.LE.(MZJZ+AMHH)) THEN
1238             TMP(3)=AMHH
1239             TMP(5)=AMHH
1240             TMP(4)=-SNIJ(JZ,IZ)
1241             TERMHH=G**2/64./PI**3/MZIZ*(MTAMZ*COSA*(XHIJ(JZ,IZ)+
1242      $            XHIJ(IZ,JZ))/AMW/COSBE)**2*
1243      $      SSXINT(MZJZ,SSZZF4,(MZIZ**2+MZJZ**2-4*AMTAU**2)/2./MZIZ)
1244             ELSE
1245             TERMHH=0.
1246             END IF
1247             IF (MZIZ.LE.(MZJZ+AMHH).AND.MZIZ.LE.(MZJZ+AMHL)) THEN
1248             TMP(3)=AMHL
1249             TMP(5)=AMHH
1250             TMP(4)=-SNIJ(JZ,IZ)
1251             TERMLH=2*G**2/64./PI**3/MZIZ*(MTAMZ/AMW/COSBE)**2*
1252      $            (COSA*(XHIJ(JZ,IZ)+XHIJ(IZ,JZ))*SINA*
1253      $             (XLIJ(JZ,IZ)+XLIJ(IZ,JZ)))*
1254      $      SSXINT(MZJZ,SSZZF4,(MZIZ**2+MZJZ**2-4*AMTAU**2)/2./MZIZ)
1255             ELSE
1256             TERMLH=0.
1257             END IF
1258             IF (MZIZ.LE.(MZJZ+AMHA)) THEN
1259             TMP(3)=AMHA
1260             TMP(5)=AMHA
1261             TMP(4)=SNIJ(JZ,IZ)
1262             TERMHA=G**2*TANB**2/64./PI**3/MZIZ*(MTAMZ*(XPIJ(JZ,IZ)+
1263      $            XPIJ(IZ,JZ))/AMW)**2*
1264      $      SSXINT(MZJZ,SSZZF4,(MZIZ**2+MZJZ**2-4*AMTAU**2)/2./MZIZ)
1265             ELSE
1266             TERMHA=0.
1267             END IF
1268             IF (MZIZ.LT.AML1SS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
1269             TMP(4)=-SGNIJ
1270             TMP(3)=AML1SS
1271             TERM1Z=8*E*WIJ(JZ,IZ)*(ALJZ1*ALIZ1*(APL-BTL)-
1272      $       BEJZ1*BEIZ1*(APL+BTL))/MZIZ
1273      $      /(2*PI)**5*SSXINT(4*AMTAU**2,SSZZF3,(MZIZ-MZJZ)**2)
1274             TM1ZRL=-8*E*WIJ(JZ,IZ)*BEJZ1*BEIZ1*(APL+BTL)/MZIZ
1275      $      /(2*PI)**5*SSXINT(4*AMTAU**2,SSZZF3,(MZIZ-MZJZ)**2)
1276             TM1ZLR=8*E*WIJ(JZ,IZ)*ALJZ1*ALIZ1*(APL-BTL)/MZIZ
1277      $      /(2*PI)**5*SSXINT(4*AMTAU**2,SSZZF3,(MZIZ-MZJZ)**2)
1278             ELSE
1279             TERM1Z=0.
1280             TM1ZRL=0.
1281             TM1ZLR=0.
1282             END IF
1283             IF (MZIZ.LT.AML2SS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
1284             TMP(4)=-SGNIJ
1285             TMP(3)=AML2SS
1286             TERM2Z=8*E*WIJ(JZ,IZ)*(ALJZ2*ALIZ2*(APL-BTL)-
1287      $       BEJZ2*BEIZ2*(APL+BTL))/MZIZ
1288      $      /(2*PI)**5*SSXINT(4*AMTAU**2,SSZZF3,(MZIZ-MZJZ)**2)
1289             TM2ZRL=-8*E*WIJ(JZ,IZ)*BEJZ2*BEIZ2*(APL+BTL)/MZIZ
1290      $      /(2*PI)**5*SSXINT(4*AMTAU**2,SSZZF3,(MZIZ-MZJZ)**2)
1291             TM2ZLR=8*E*WIJ(JZ,IZ)*ALJZ2*ALIZ2*(APL-BTL)/MZIZ
1292      $      /(2*PI)**5*SSXINT(4*AMTAU**2,SSZZF3,(MZIZ-MZJZ)**2)
1293             ELSE
1294             TERM2Z=0.
1295             TM2ZRL=0.
1296             TM2ZLR=0.
1297             END IF
1298             IF (MZIZ.LT.AML1SS.AND.MZIZ.LT.(MZJZ+AMHL)) THEN
1299             TMP(3)=AML1SS
1300             TMP(4)=SGNIJ
1301             TMP(5)=AMHL
1302             TERM1L=2*PI*PI*G*MTAMZ*SINA*SGNIJ*(XLIJ(IZ,JZ)+
1303      $        XLIJ(JZ,IZ))*(ALIZ1*BEJZ1+ALJZ1*BEIZ1)
1304      $      /MZIZ/AMW/COSBE*SSXINT(4*AMTAU**2,SSZZF5,(MZIZ-MZJZ)**2)
1305             TERM1L=FACT*TERM1L
1306             ELSE
1307             TERM1L=0.
1308             END IF
1309             IF (MZIZ.LT.AML2SS.AND.MZIZ.LT.(MZJZ+AMHL)) THEN
1310             TMP(3)=AML2SS
1311             TMP(4)=SGNIJ
1312             TMP(5)=AMHL
1313             TERM2L=2*PI*PI*G*MTAMZ*SINA*SGNIJ*(XLIJ(IZ,JZ)+
1314      $        XLIJ(JZ,IZ))*(ALIZ2*BEJZ2+ALJZ2*BEIZ2)
1315      $      /MZIZ/AMW/COSBE*SSXINT(4*AMTAU**2,SSZZF5,(MZIZ-MZJZ)**2)
1316             TERM2L=FACT*TERM2L
1317             ELSE
1318             TERM2L=0.
1319             END IF
1320             IF (MZIZ.LT.AML1SS.AND.MZIZ.LT.(MZJZ+AMHH)) THEN
1321             TMP(3)=AML1SS
1322             TMP(4)=SGNIJ
1323             TMP(5)=AMHH
1324             TERM1H=2*PI*PI*G*MTAMZ*COSA*SGNIJ*(XHIJ(IZ,JZ)+
1325      $        XHIJ(JZ,IZ))*(ALIZ1*BEJZ1+ALJZ1*BEIZ1)
1326      $      /MZIZ/AMW/COSBE*SSXINT(4*AMTAU**2,SSZZF5,(MZIZ-MZJZ)**2)
1327             TERM1H=FACT*TERM1H
1328             ELSE
1329             TERM1H=0.
1330             END IF
1331             IF (MZIZ.LT.AML2SS.AND.MZIZ.LT.(MZJZ+AMHH)) THEN
1332             TMP(3)=AML2SS
1333             TMP(4)=SGNIJ
1334             TMP(5)=AMHH
1335             TERM2H=2*PI*PI*G*MTAMZ*COSA*SGNIJ*(XHIJ(IZ,JZ)+
1336      $        XHIJ(JZ,IZ))*(ALIZ2*BEJZ2+ALJZ2*BEIZ2)
1337      $    /MZIZ/AMW/COSBE*SSXINT(4*AMTAU**2,SSZZF5,(MZIZ-MZJZ)**2)
1338             TERM2H=FACT*TERM2H
1339             ELSE
1340             TERM2H=0.
1341             END IF
1342             IF (MZIZ.LT.AML1SS.AND.MZIZ.LT.(MZJZ+AMHA)) THEN
1343             TMP(3)=AML1SS
1344             TMP(4)=-SGNIJ
1345             TMP(5)=AMHA
1346             TERM1A=-2*PI*PI*G*MTAMZ*TANB*SGNIJ*(XPIJ(IZ,JZ)+
1347      $        XPIJ(JZ,IZ))*(ALIZ1*BEJZ1+ALJZ1*BEIZ1)
1348      $       /MZIZ/AMW*SSXINT(4*AMTAU**2,SSZZF5,(MZIZ-MZJZ)**2)
1349             TERM1A=FACT*TERM1A
1350             ELSE
1351             TERM1A=0.
1352             END IF
1353             IF (MZIZ.LT.AML2SS.AND.MZIZ.LT.(MZJZ+AMHA)) THEN
1354             TMP(3)=AML2SS
1355             TMP(4)=-SGNIJ
1356             TMP(5)=AMHA
1357             TERM2A=-2*PI*PI*G*MTAMZ*TANB*SGNIJ*(XPIJ(IZ,JZ)+
1358      $        XPIJ(JZ,IZ))*(ALIZ2*BEJZ2+ALJZ2*BEIZ2)
1359      $       /MZIZ/AMW*SSXINT(4*AMTAU**2,SSZZF5,(MZIZ-MZJZ)**2)
1360             TERM2A=FACT*TERM2A
1361             ELSE
1362             TERM2A=0.
1363             END IF
1364             WID=GF+TERMZZ+TERMHL+TERMHH+TERMLH+TERMHA+TERM1Z+
1365      $        TERM2Z+TERM1L+TERM2L+TERM1H+TERM2H+TERM1A+TERM2A
1366 C              tau polarization for 3-body z2 -> z1 tau tau
1367             IF (IZ.EQ.2.AND.JZ.EQ.1.AND.WID.GT.0.) THEN
1368               PTAUZZ=(GRRF1R+GRGRR+GRRF2R+TMZZRL+TM1ZRL+TM2ZRL-
1369      $               (GLLF1L+GLGLL+GLLF2L+TMZZLR+TM1ZLR+TM2ZLR))
1370      $               /WID
1371             END IF
1372             CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDTAU,-IDTAU,0,0)
1373 C     Enter information for decay matrix element
1374             IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN
1375              Z1(1)=0.
1376              Z1(2)=2*E*WIJ(JZ,IZ)
1377             ELSE
1378              Z1(1)=2*E*WIJ(JZ,IZ)
1379              Z1(2)=0.
1380             END IF
1381             Z2(1)=APL
1382             Z2(2)=BTL
1383             CALL SSME3(1,AMZ,Z1,Z2)
1384             Z1(1)=((ZI**(THIZ)*(-1.)**(THIZ+1)*ALI(IZ)-
1385      $            FL*ZMIXSS(2,IZ)*ZI**THIZ)*COSL-(ZI**THIZ*BLI(IZ)-
1386      $            FL*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINL)/2.
1387             Z1(2)=((-ZI**(THIZ)*(-1.)**(THIZ+1)*ALI(IZ)-
1388      $            FL*ZMIXSS(2,IZ)*ZI**THIZ)*COSL-(ZI**THIZ*BLI(IZ)+
1389      $            FL*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINL)/2.
1390             Z2(1)=CONJG(((ZI**(THJZ)*(-1.)**(THJZ+1)*ALI(JZ)-
1391      $            FL*ZMIXSS(2,JZ)*ZI**THJZ)*COSL-(ZI**THJZ*BLI(JZ)-
1392      $            FL*ZMIXSS(2,JZ)*(-ZI)**THJZ)*SINL)/2.)
1393             Z2(2)=-CONJG(((-ZI**(THJZ)*(-1.)**(THJZ+1)*ALI(JZ)-
1394      $            FL*ZMIXSS(2,JZ)*ZI**THJZ)*COSL-(ZI**THJZ*BLI(JZ)+
1395      $            FL*ZMIXSS(2,JZ)*(-ZI)**THJZ)*SINL)/2.)
1396             CALL SSME3(2,AML1SS,Z1,Z2)
1397             Z1(1)=((ZI**(THIZ)*(-1.)**(THIZ+1)*ALI(IZ)-
1398      $            FL*ZMIXSS(2,IZ)*ZI**THIZ)*SINL+(ZI**THIZ*BLI(IZ)-
1399      $            FL*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSL)/2.
1400             Z1(2)=((-ZI**(THIZ)*(-1.)**(THIZ+1)*ALI(IZ)-
1401      $            FL*ZMIXSS(2,IZ)*ZI**THIZ)*SINL+(ZI**THIZ*BLI(IZ)+
1402      $            FL*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSL)/2.
1403             Z2(1)=CONJG(((ZI**(THJZ)*(-1.)**(THJZ+1)*ALI(JZ)-
1404      $            FL*ZMIXSS(2,JZ)*ZI**THJZ)*SINL+(ZI**THJZ*BLI(JZ)-
1405      $            FL*ZMIXSS(2,JZ)*(-ZI)**THJZ)*COSL)/2.)
1406             Z2(2)=-CONJG(((-ZI**(THJZ)*(-1.)**(THJZ+1)*ALI(JZ)-
1407      $            FL*ZMIXSS(2,JZ)*ZI**THJZ)*SINL+(ZI**THJZ*BLI(JZ)+
1408      $            FL*ZMIXSS(2,JZ)*(-ZI)**THJZ)*COSL)/2.)
1409             CALL SSME3(2,AML2SS,Z1,Z2)
1410             Z1(1)=((ZI**(THJZ)*(-1.)**(THJZ+1)*ALI(JZ)-
1411      $            FL*ZMIXSS(2,JZ)*ZI**THJZ)*COSL-(ZI**THJZ*BLI(JZ)-
1412      $            FL*ZMIXSS(2,JZ)*(-ZI)**THJZ)*SINL)/2.
1413             Z1(2)=((-ZI**(THJZ)*(-1.)**(THJZ+1)*ALI(JZ)-
1414      $            FL*ZMIXSS(2,JZ)*ZI**THJZ)*COSL-(ZI**THJZ*BLI(JZ)+
1415      $            FL*ZMIXSS(2,JZ)*(-ZI)**THJZ)*SINL)/2.
1416             Z2(1)=-CONJG(((ZI**(THIZ)*(-1.)**(THIZ+1)*ALI(IZ)-
1417      $            FL*ZMIXSS(2,IZ)*ZI**THIZ)*COSL-(ZI**THIZ*BLI(IZ)-
1418      $            FL*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINL)/2.)
1419             Z2(2)=CONJG(((-ZI**(THIZ)*(-1.)**(THIZ+1)*ALI(IZ)-
1420      $            FL*ZMIXSS(2,IZ)*ZI**THIZ)*COSL-(ZI**THIZ*BLI(IZ)+
1421      $            FL*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINL)/2.)
1422             CALL SSME3(3,AML1SS,Z1,Z2)
1423             Z1(1)=((ZI**(THJZ)*(-1.)**(THJZ+1)*ALI(JZ)-
1424      $            FL*ZMIXSS(2,JZ)*ZI**THJZ)*SINL+(ZI**THJZ*BLI(JZ)-
1425      $            FL*ZMIXSS(2,JZ)*(-ZI)**THJZ)*COSL)/2.
1426             Z1(2)=((-ZI**(THJZ)*(-1.)**(THJZ+1)*ALI(JZ)-
1427      $            FL*ZMIXSS(2,JZ)*ZI**THJZ)*SINL+(ZI**THJZ*BLI(JZ)+
1428      $            FL*ZMIXSS(2,JZ)*(-ZI)**THJZ)*COSL)/2.
1429             Z2(1)=-CONJG(((ZI**(THIZ)*(-1.)**(THIZ+1)*ALI(IZ)-
1430      $            FL*ZMIXSS(2,IZ)*ZI**THIZ)*SINL+(ZI**THIZ*BLI(IZ)-
1431      $            FL*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSL)/2.)
1432             Z2(2)=CONJG(((-ZI**(THIZ)*(-1.)**(THIZ+1)*ALI(IZ)-
1433      $            FL*ZMIXSS(2,IZ)*ZI**THIZ)*SINL+(ZI**THIZ*BLI(IZ)+
1434      $            FL*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSL)/2.)
1435             CALL SSME3(3,AML2SS,Z1,Z2)
1436             IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN
1437              Z1(1)=(-ZI)**(THIZ+THJZ)*(XLIJ(IZ,JZ)+XLIJ(JZ,IZ))
1438              Z1(2)=0.
1439             ELSE
1440              Z1(1)=0.
1441              Z1(2)=(-ZI)**(THIZ+THJZ)*(XLIJ(IZ,JZ)+XLIJ(JZ,IZ))
1442             END IF
1443             Z2(1)=-G*MTAMZ*SINA/2./AMW/COSBE
1444             Z2(2)=0.
1445             CALL SSME3(4,AMHL,Z1,Z2)
1446             IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN
1447              Z1(1)=(-ZI)**(THIZ+THJZ)*(XHIJ(IZ,JZ)+XHIJ(JZ,IZ))
1448              Z1(2)=0.
1449             ELSE
1450              Z1(1)=0.
1451              Z1(2)=(-ZI)**(THIZ+THJZ)*(XHIJ(IZ,JZ)+XHIJ(JZ,IZ))
1452             END IF
1453             Z2(1)=-G*MTAMZ*COSA/2./AMW/COSBE
1454             Z2(2)=0.
1455             CALL SSME3(4,AMHH,Z1,Z2)
1456             IF (-SNIJ(JZ,IZ).GT.0.) THEN
1457              Z1(1)=0.
1458              Z1(2)=(-ZI)**(THIZ+THJZ+1)*(XPIJ(IZ,JZ)+XPIJ(JZ,IZ))
1459             ELSE
1460              Z1(1)=(-ZI)**(THIZ+THJZ+1)*(XPIJ(IZ,JZ)+XPIJ(JZ,IZ))
1461              Z1(2)=0.
1462             END IF
1463             Z2(1)=0.
1464             Z2(2)=ZI*G*MTAMZ*TANB/2./AMW
1465             CALL SSME3(4,AMHA,Z1,Z2)
1466           END IF
1467 C          zi -> zj + nu_e + nu_e bar
1468           TMP(6)=0.
1469           IF (MZIZ.GT.MZJZ) THEN
1470             IF (MZIZ.LT.AMN1SS) THEN
1471             TMP(4)=-SNIJ(JZ,IZ)
1472             TMP(3)=AMN1SS
1473             TERMLL=2*ANI(JZ)**2*ANI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.)
1474             ELSE
1475             TERMLL=0.
1476             END IF
1477             IF (MZIZ.LT.(MZJZ+AMZ)) THEN
1478             TMP(4)=+SNIJ(JZ,IZ)
1479             TERMZZ=(APN**2+BTN**2)*E**2*WIJ(JZ,IZ)**2/MZIZ
1480      $             *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ)
1481             ELSE
1482             TERMZZ=0.
1483             END IF
1484             IF (MZIZ.LT.AMN1SS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
1485             TMP(4)=+SNIJ(JZ,IZ)
1486             TMP(3)=AMN1SS
1487             TERMLZ=8*E*(APN-BTN)*WIJ(JZ,IZ)*ANI(JZ)*ANI(IZ)/MZIZ
1488      $      /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2)
1489             ELSE
1490             TERMLZ=0.
1491             END IF
1492             WID=TERMLL+TERMZZ+TERMLZ
1493             CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDNE,-IDNE,0,0)
1494 C     Enter information for decay matrix element
1495             IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN
1496              Z1(1)=0.
1497              Z1(2)=2*E*WIJ(JZ,IZ)
1498             ELSE
1499              Z1(1)=2*E*WIJ(JZ,IZ)
1500              Z1(2)=0.
1501             END IF
1502             Z2(1)=APN
1503             Z2(2)=BTN
1504             CALL SSME3(1,AMZ,Z1,Z2)
1505             Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*ANI(IZ)/2.
1506             Z1(2)=-Z1(1)
1507             Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2.)
1508             Z2(2)=Z2(1)
1509             CALL SSME3(2,AMN1SS,Z1,Z2)
1510             Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2.
1511             Z1(2)=-Z1(1)
1512             Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*ANI(IZ)/2.)
1513             Z2(2)=Z2(1)
1514             CALL SSME3(3,AMN1SS,Z1,Z2)
1515           END IF
1516 C          zi -> zj + nu_mu + nu_mu bar
1517           IF (MZIZ.GT.MZJZ) THEN
1518             IF (MZIZ.LT.AMN2SS) THEN
1519             TMP(4)=-SNIJ(JZ,IZ)
1520             TMP(3)=AMN2SS
1521             TERMLL=2*ANI(JZ)**2*ANI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.)
1522             ELSE
1523             TERMLL=0.
1524             END IF
1525             IF (MZIZ.LT.(MZJZ+AMZ)) THEN
1526             TMP(4)=+SNIJ(JZ,IZ)
1527             TERMZZ=(APN**2+BTN**2)*E**2*WIJ(JZ,IZ)**2/MZIZ
1528      $             *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ)
1529             ELSE
1530             TERMZZ=0.
1531             END IF
1532             IF (MZIZ.LT.AMN2SS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
1533             TMP(4)=+SNIJ(JZ,IZ)
1534             TMP(3)=AMN2SS
1535             TERMLZ=8*E*(APN-BTN)*WIJ(JZ,IZ)*ANI(JZ)*ANI(IZ)/MZIZ
1536      $      /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2)
1537             ELSE
1538             TERMLZ=0.
1539             END IF
1540             WID=TERMLL+TERMZZ+TERMLZ
1541             CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDNM,-IDNM,0,0)
1542 C     Enter information for decay matrix element
1543             IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN
1544              Z1(1)=0.
1545              Z1(2)=2*E*WIJ(JZ,IZ)
1546             ELSE
1547              Z1(1)=2*E*WIJ(JZ,IZ)
1548              Z1(2)=0.
1549             END IF
1550             Z2(1)=APN
1551             Z2(2)=BTN
1552             CALL SSME3(1,AMZ,Z1,Z2)
1553             Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*ANI(IZ)/2.
1554             Z1(2)=-Z1(1)
1555             Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2.)
1556             Z2(2)=Z2(1)
1557             CALL SSME3(2,AMN2SS,Z1,Z2)
1558             Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2.
1559             Z1(2)=-Z1(1)
1560             Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*ANI(IZ)/2.)
1561             Z2(2)=Z2(1)
1562             CALL SSME3(3,AMN2SS,Z1,Z2)
1563           END IF
1564 C          zi -> zj + nu_tau + nu_tau bar
1565           IF (MZIZ.GT.MZJZ) THEN
1566             IF (MZIZ.LT.AMN3SS) THEN
1567             TMP(4)=-SNIJ(JZ,IZ)
1568             TMP(3)=AMN3SS
1569             TERMLL=2*ANI(JZ)**2*ANI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.)
1570             ELSE
1571             TERMLL=0.
1572             END IF
1573             IF (MZIZ.LT.(MZJZ+AMZ)) THEN
1574             TMP(4)=+SNIJ(JZ,IZ)
1575             TERMZZ=(APN**2+BTN**2)*E**2*WIJ(JZ,IZ)**2/MZIZ
1576      $             *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ)
1577             ELSE
1578             TERMZZ=0.
1579             END IF
1580             IF (MZIZ.LT.AMN3SS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN
1581             TMP(4)=+SNIJ(JZ,IZ)
1582             TMP(3)=AMN3SS
1583             TERMLZ=8*E*(APN-BTN)*WIJ(JZ,IZ)*ANI(JZ)*ANI(IZ)/MZIZ
1584      $      /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2)
1585             ELSE
1586             TERMLZ=0.
1587             END IF
1588             WID=TERMLL+TERMZZ+TERMLZ
1589             CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDNT,-IDNT,0,0)
1590 C     Enter information for decay matrix element
1591             IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN
1592              Z1(1)=0.
1593              Z1(2)=2*E*WIJ(JZ,IZ)
1594             ELSE
1595              Z1(1)=2*E*WIJ(JZ,IZ)
1596              Z1(2)=0.
1597             END IF
1598             Z2(1)=APN
1599             Z2(2)=BTN
1600             CALL SSME3(1,AMZ,Z1,Z2)
1601             Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*ANI(IZ)/2.
1602             Z1(2)=-Z1(1)
1603             Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2.)
1604             Z2(2)=Z2(1)
1605             CALL SSME3(2,AMN3SS,Z1,Z2)
1606             Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2.
1607             Z1(2)=-Z1(1)
1608             Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*ANI(IZ)/2.)
1609             Z2(2)=Z2(1)
1610             CALL SSME3(3,AMN3SS,Z1,Z2)
1611           END IF
1612 220     CONTINUE
1613 C
1614 C          zi --> zj + higgs
1615 C
1616         DO 230 JZ=1,IZ-1
1617 C          zi --> zj + hl
1618           MZJZ=ABS(AMZISS(JZ))
1619           IF (MZIZ.GT.(MZJZ+AMHL)) THEN
1620             SN=SIGN(1.,AMZISS(JZ))*SIGN(1.,AMZISS(IZ))
1621             WID=(XLIJ(JZ,IZ)+XLIJ(IZ,JZ))**2/8./PI/(MZIZ)**3
1622      $      *SQRT(MZIZ**4+MZJZ**4+AMHL**4-2.*(MZIZ*MZJZ)**2
1623      $      -2.*(MZIZ*AMHL)**2-2.*(MZJZ*AMHL)**2)*((MZIZ**2+MZJZ**2
1624      $      -AMHL**2)/2.+SN*MZIZ*MZJZ)
1625             CALL SSSAVE(ISZIZ,WID,ISZ(JZ),ISHL,0,0,0)
1626           END IF
1627 C          zi --> zj + hh
1628           IF (MZIZ.GT.(MZJZ+AMHH)) THEN
1629             SN=SIGN(1.,AMZISS(JZ))*SIGN(1.,AMZISS(IZ))
1630             WID=(HIJ(JZ,IZ)+HIJ(IZ,JZ))**2/8./PI/(MZIZ)**3
1631      $      *SQRT(MZIZ**4+MZJZ**4+AMHH**4-2.*(MZIZ*MZJZ)**2
1632      $      -2.*(MZIZ*AMHH)**2-2.*(MZJZ*AMHH)**2)
1633      $      *((MZIZ**2+MZJZ**2-AMHH**2)/2.+SN*MZIZ*MZJZ)
1634             CALL SSSAVE(ISZIZ,WID,ISZ(JZ),ISHH,0,0,0)
1635           END IF
1636 C          zi --> zj + ha
1637           IF (MZIZ.GT.(MZJZ+AMHA)) THEN
1638             SN=-SIGN(1.,AMZISS(JZ))*SIGN(1.,AMZISS(IZ))
1639             WID=(XPIJ(IZ,JZ)+XPIJ(JZ,IZ))**2/8./PI/(MZIZ)**3
1640      $      *SQRT(MZIZ**4+MZJZ**4+AMHA**4-2.*(MZIZ*MZJZ)**2
1641      $      -2.*(MZIZ*AMHA)**2-2.*(MZJZ*AMHA)**2)*((MZIZ**2+MZJZ**2
1642      $      -AMHA**2)/2.+SN*MZIZ*MZJZ)
1643             CALL SSSAVE(ISZIZ,WID,ISZ(JZ),ISHA,0,0,0)
1644           END IF
1645 230     CONTINUE
1646 200   CONTINUE
1647 C
1648 C          zi --> squark + qbar; enlarge to include Z1 decays
1649 C                 in case of models with light gravitino
1650 C
1651       DO 245 IZ=1,4
1652         MZIZ=ABS(AMZISS(IZ))
1653         SNIZ=SIGN(1.,AMZISS(IZ))
1654         IF (SNIZ.EQ.1.) THEN
1655            THIZ=0
1656         ELSE
1657            THIZ=1
1658         END IF
1659         ISZIZ=ISZ(IZ)
1660 C          left squarks
1661         IF (MZIZ.GT.(AMULSS+AMUP)) THEN
1662           WID=3*AUI(IZ)**2*(MZIZ**2+AMUP**2-AMULSS**2)/MZIZ**3/
1663      $         32./PI*SQRT(SSXLAM(MZIZ**2,AMUP**2,AMULSS**2))
1664           CALL SSSAVE(ISZIZ,WID,ISUPL,-IDUP,0,0,0)
1665           CALL SSSAVE(ISZIZ,WID,-ISUPL,IDUP,0,0,0)
1666         ENDIF
1667         IF (MZIZ.GT.(AMDLSS+AMDN)) THEN
1668           WID=3*ADI(IZ)**2*(MZIZ**2+AMDN**2-AMDLSS**2)/MZIZ**3/
1669      $         32./PI*SQRT(SSXLAM(MZIZ**2,AMDN**2,AMDLSS**2))
1670           CALL SSSAVE(ISZIZ,WID,ISDNL,-IDDN,0,0,0)
1671           CALL SSSAVE(ISZIZ,WID,-ISDNL,IDDN,0,0,0)
1672         END IF
1673         IF (MZIZ.GT.(AMSLSS+AMST)) THEN
1674           WID=3*ADI(IZ)**2*(MZIZ**2+AMST**2-AMSLSS**2)/MZIZ**3/
1675      $         32./PI*SQRT(SSXLAM(MZIZ**2,AMST**2,AMSLSS**2))
1676           CALL SSSAVE(ISZIZ,WID,ISSTL,-IDST,0,0,0)
1677           CALL SSSAVE(ISZIZ,WID,-ISSTL,IDST,0,0,0)
1678         END IF
1679         IF (MZIZ.GT.(AMCLSS+AMCH)) THEN
1680           WID=3*AUI(IZ)**2*(MZIZ**2+AMCH**2-AMCLSS**2)/MZIZ**3/
1681      $         32./PI*SQRT(SSXLAM(MZIZ**2,AMCH**2,AMCLSS**2))
1682           CALL SSSAVE(ISZIZ,WID,ISCHL,-IDCH,0,0,0)
1683           CALL SSSAVE(ISZIZ,WID,-ISCHL,IDCH,0,0,0)
1684         ENDIF
1685 C          right squarks
1686         IF (MZIZ.GT.(AMURSS+AMUP)) THEN
1687           WID=3*BUI(IZ)**2*(MZIZ**2+AMUP**2-AMURSS**2)/MZIZ**3/
1688      $         32./PI*SQRT(SSXLAM(MZIZ**2,AMUP**2,AMURSS**2))
1689           CALL SSSAVE(ISZIZ,WID,ISUPR,-IDUP,0,0,0)
1690           CALL SSSAVE(ISZIZ,WID,-ISUPR,IDUP,0,0,0)
1691         END IF
1692         IF (MZIZ.GT.(AMDRSS+AMDN)) THEN
1693           WID=3*BDI(IZ)**2*(MZIZ**2+AMDN**2-AMDRSS**2)/MZIZ**3/
1694      $         32./PI*SQRT(SSXLAM(MZIZ**2,AMDN**2,AMDRSS**2))
1695           CALL SSSAVE(ISZIZ,WID,ISDNR,-IDDN,0,0,0)
1696           CALL SSSAVE(ISZIZ,WID,-ISDNR,IDDN,0,0,0)
1697         END IF
1698         IF (MZIZ.GT.(AMSRSS+AMST)) THEN
1699           WID=3*BDI(IZ)**2*(MZIZ**2+AMST**2-AMSRSS**2)/MZIZ**3/
1700      $         32./PI*SQRT(SSXLAM(MZIZ**2,AMST**2,AMSRSS**2))
1701           CALL SSSAVE(ISZIZ,WID,ISSTR,-IDST,0,0,0)
1702           CALL SSSAVE(ISZIZ,WID,-ISSTR,IDST,0,0,0)
1703         END IF
1704         IF(MZIZ.GT.(AMCRSS+AMCH)) THEN
1705           WID=3*BUI(IZ)**2*(MZIZ**2+AMCH**2-AMCRSS**2)/MZIZ**3/
1706      $         32./PI*SQRT(SSXLAM(MZIZ**2,AMCH**2,AMCRSS**2))
1707             CALL SSSAVE(ISZIZ,WID,ISCHR,-IDCH,0,0,0)
1708             CALL SSSAVE(ISZIZ,WID,-ISCHR,IDCH,0,0,0)
1709         ENDIF
1710 C          z_i --> sbottom + bottom
1711 C
1712         ZAUIZ=ZI**(THIZ-1)*SNIZ
1713      $  *(-G/SR2*ZMIXSS(3,IZ)-GP/3./SR2*ZMIXSS(4,IZ))
1714         ZBUIZ=ZI**(THIZ-1)*4*GP*ZMIXSS(4,IZ)/3./SR2
1715         ZADIZ=ZI**(THIZ-1)*SNIZ
1716      $  *(G/SR2*ZMIXSS(3,IZ)-GP/3./SR2*ZMIXSS(4,IZ))
1717         ZBDIZ=-2*ZI**(THIZ-1)*GP*ZMIXSS(4,IZ)/3./SR2
1718         ZALIZ=ZI**(THIZ-1)*SNIZ
1719      $  *(G/SR2*ZMIXSS(3,IZ)+GP/SR2*ZMIXSS(4,IZ))
1720         ZBLIZ=-1*ZI**(THIZ-1)*SR2*GP*ZMIXSS(4,IZ)
1721         ZPP=ZI**THIZ
1722         ZPM=(-ZI)**THIZ
1723         IF(MZIZ.GT.(AMB1SS+AMBT)) THEN
1724           ZA=((ZI*ZADIZ-ZPP*FB*ZMIXSS(2,IZ))*COSB
1725      $     -(ZI*ZBDIZ-ZPM*FB*ZMIXSS(2,IZ))*SINB)/2.
1726           ZB=((-ZI*ZADIZ-ZPP*FB*ZMIXSS(2,IZ))*COSB
1727      $     -(ZI*ZBDIZ+ZPM*FB*ZMIXSS(2,IZ))*SINB)/2.
1728           AS=ZA*CONJG(ZA)
1729           BS=ZB*CONJG(ZB)
1730           WID=3*(AS*((AMBT+MZIZ)**2-AMB1SS**2)+BS*((MZIZ-AMBT)**2-
1731      $     AMB1SS**2))/16./PI/MZIZ**3*
1732      $     SQRT(SSXLAM(MZIZ**2,AMB1SS**2,AMBT**2))
1733             CALL SSSAVE(ISZIZ,WID,ISBT1,-IDBT,0,0,0)
1734             CALL SSSAVE(ISZIZ,WID,-ISBT1,IDBT,0,0,0)
1735         ENDIF
1736         IF(MZIZ.GT.(AMB2SS+AMBT)) THEN
1737           ZA=((ZI*ZADIZ-ZPP*FB*ZMIXSS(2,IZ))*SINB
1738      $     +(ZI*ZBDIZ-ZPM*FB*ZMIXSS(2,IZ))*COSB)/2.
1739           ZB=((-ZI*ZADIZ-ZPP*FB*ZMIXSS(2,IZ))*SINB
1740      $     +(ZI*ZBDIZ+ZPM*FB*ZMIXSS(2,IZ))*COSB)/2.
1741           AS=ZA*CONJG(ZA)
1742           BS=ZB*CONJG(ZB)
1743           WID=3*(AS*((AMBT+MZIZ)**2-AMB2SS**2)+BS*((MZIZ-AMBT)**2-
1744      $     AMB2SS**2))/16./PI/MZIZ**3*
1745      $     SQRT(SSXLAM(MZIZ**2,AMB2SS**2,AMBT**2))
1746             CALL SSSAVE(ISZIZ,WID,ISBT2,-IDBT,0,0,0)
1747             CALL SSSAVE(ISZIZ,WID,-ISBT2,IDBT,0,0,0)
1748         ENDIF
1749 C          z_i --> stop + top
1750 C
1751         IF(MZIZ.GT.AMT1SS+AMTP) THEN
1752           ZA=((ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*COST
1753      $     -(ZI*ZBUIZ-ZPM*FT*ZMIXSS(1,IZ))*SINT)/2.
1754           ZB=((-ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*COST
1755      $     -(ZI*ZBUIZ+ZPM*FT*ZMIXSS(1,IZ))*SINT)/2.
1756           AS=ZA*CONJG(ZA)
1757           BS=ZB*CONJG(ZB)
1758           WID=3*(AS*((AMTP+MZIZ)**2-AMT1SS**2)+BS*((MZIZ-AMTP)**2-
1759      $     AMT1SS**2))/16./PI/MZIZ**3*
1760      $     SQRT(SSXLAM(MZIZ**2,AMT1SS**2,AMTP**2))
1761           CALL SSSAVE(ISZIZ,WID,ISTP1,-IDTP,0,0,0)
1762           CALL SSSAVE(ISZIZ,WID,-ISTP1,IDTP,0,0,0)
1763         ENDIF
1764         IF(MZIZ.GT.AMT2SS+AMTP) THEN
1765           ZA=((ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*SINT
1766      $     +(ZI*ZBUIZ-ZPM*FT*ZMIXSS(1,IZ))*COST)/2.
1767           ZB=((-ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*SINT
1768      $     +(ZI*ZBUIZ+ZPM*FT*ZMIXSS(1,IZ))*COST)/2.
1769           AS=ZA*CONJG(ZA)
1770           BS=ZB*CONJG(ZB)
1771           WID=3*(AS*((AMTP+MZIZ)**2-AMT2SS**2)+BS*((MZIZ-AMTP)**2-
1772      $     AMT2SS**2))/16./PI/MZIZ**3*
1773      $     SQRT(SSXLAM(MZIZ**2,AMT2SS**2,AMTP**2))
1774           CALL SSSAVE(ISZIZ,WID,ISTP2,-IDTP,0,0,0)
1775           CALL SSSAVE(ISZIZ,WID,-ISTP2,IDTP,0,0,0)
1776         ENDIF
1777 C
1778 C          zi --> slepton + lepton
1779 C
1780         IF(MZIZ.GT.(AMELSS+AME)) THEN
1781           WID=ALI(IZ)**2*(MZIZ**2+AME**2-AMELSS**2)/MZIZ**3/
1782      $         32./PI*SQRT(SSXLAM(MZIZ**2,AME**2,AMELSS**2))
1783           CALL SSSAVE(ISZIZ,WID,ISEL,-IDE,0,0,0)
1784           CALL SSSAVE(ISZIZ,WID,-ISEL,IDE,0,0,0)
1785         END IF
1786         IF(MZIZ.GT.(AMMLSS+AMMU)) THEN
1787           WID=ALI(IZ)**2*(MZIZ**2+AMMU**2-AMMLSS**2)/MZIZ**3/
1788      $         32./PI*SQRT(SSXLAM(MZIZ**2,AMMU**2,AMMLSS**2))
1789           CALL SSSAVE(ISZIZ,WID,ISMUL,-IDMU,0,0,0)
1790           CALL SSSAVE(ISZIZ,WID,-ISMUL,IDMU,0,0,0)
1791         END IF
1792         IF(MZIZ.GT.(AMERSS+AME)) THEN
1793           WID=BLI(IZ)**2*(MZIZ**2+AME**2-AMERSS**2)/MZIZ**3/
1794      $         32./PI*SQRT(SSXLAM(MZIZ**2,AME**2,AMERSS**2))
1795           CALL SSSAVE(ISZIZ,WID,ISER,-IDE,0,0,0)
1796           CALL SSSAVE(ISZIZ,WID,-ISER,IDE,0,0,0)
1797         END IF
1798         IF(MZIZ.GT.(AMMRSS+AMMU)) THEN
1799           WID=BLI(IZ)**2*(MZIZ**2+AMMU**2-AMMRSS**2)/MZIZ**3/
1800      $         32./PI*SQRT(SSXLAM(MZIZ**2,AMMU**2,AMMRSS**2))
1801           CALL SSSAVE(ISZIZ,WID,ISMUR,-IDMU,0,0,0)
1802           CALL SSSAVE(ISZIZ,WID,-ISMUR,IDMU,0,0,0)
1803         END IF
1804         IF(MZIZ.GT.(AML1SS+AMTAU)) THEN
1805           ZA=((ZI*ZALIZ-ZPP*FL*ZMIXSS(2,IZ))*COSL
1806      $     -(ZI*ZBLIZ-ZPM*FL*ZMIXSS(2,IZ))*SINL)/2.
1807           ZB=((-ZI*ZALIZ-ZPP*FL*ZMIXSS(2,IZ))*COSL
1808      $     -(ZI*ZBLIZ+ZPM*FL*ZMIXSS(2,IZ))*SINL)/2.
1809           AS=ZA*CONJG(ZA)
1810           BS=ZB*CONJG(ZB)
1811           WID=(AS*((AMTAU+MZIZ)**2-AML1SS**2)+BS*((MZIZ-AMTAU)**2-
1812      $     AML1SS**2))/16./PI/MZIZ**3*
1813      $     SQRT(SSXLAM(MZIZ**2,AML1SS**2,AMTAU**2))
1814             CALL SSSAVE(ISZIZ,WID,ISTAU1,-IDTAU,0,0,0)
1815             CALL SSSAVE(ISZIZ,WID,-ISTAU1,IDTAU,0,0,0)
1816         END IF
1817         IF(MZIZ.GT.(AML2SS+AMTAU)) THEN
1818           ZA=((ZI*ZALIZ-ZPP*FL*ZMIXSS(2,IZ))*SINL
1819      $     +(ZI*ZBLIZ-ZPM*FL*ZMIXSS(2,IZ))*COSL)/2.
1820           ZB=((-ZI*ZALIZ-ZPP*FL*ZMIXSS(2,IZ))*SINL
1821      $     +(ZI*ZBLIZ+ZPM*FL*ZMIXSS(2,IZ))*COSL)/2.
1822           AS=ZA*CONJG(ZA)
1823           BS=ZB*CONJG(ZB)
1824           WID=(AS*((AMTAU+MZIZ)**2-AML2SS**2)+BS*((MZIZ-AMTAU)**2-
1825      $     AML2SS**2))/16./PI/MZIZ**3*
1826      $     SQRT(SSXLAM(MZIZ**2,AML2SS**2,AMTAU**2))
1827             CALL SSSAVE(ISZIZ,WID,ISTAU2,-IDTAU,0,0,0)
1828             CALL SSSAVE(ISZIZ,WID,-ISTAU2,IDTAU,0,0,0)
1829         ENDIF
1830         IF(MZIZ.GT.AMN1SS) THEN
1831           WID=(MZIZ**2-AMN1SS**2)**2/32./PI/MZIZ**3
1832           CALL SSSAVE(ISZIZ,ANI(IZ)**2*WID,ISNEL,-IDNE,0,0,0)
1833           CALL SSSAVE(ISZIZ,ANI(IZ)**2*WID,-ISNEL,IDNE,0,0,0)
1834         END IF
1835         IF(MZIZ.GT.AMN2SS) THEN
1836           WID=(MZIZ**2-AMN2SS**2)**2/32./PI/MZIZ**3
1837           CALL SSSAVE(ISZIZ,ANI(IZ)**2*WID,ISNML,-IDNM,0,0,0)
1838           CALL SSSAVE(ISZIZ,ANI(IZ)**2*WID,-ISNML,IDNM,0,0,0)
1839         END IF
1840         IF(MZIZ.GT.AMN3SS) THEN
1841           WID=(MZIZ**2-AMN3SS**2)**2/32./PI/MZIZ**3
1842           CALL SSSAVE(ISZIZ,ANI(IZ)**2*WID,ISNTL,-IDNT,0,0,0)
1843           CALL SSSAVE(ISZIZ,ANI(IZ)**2*WID,-ISNTL,IDNT,0,0,0)
1844         END IF
1845 245   CONTINUE
1846 C
1847 C          Compute decays to gravitino for GMSB models
1848 C
1849       DO 250 IZ=1,4
1850         MZIZ=ABS(AMZISS(IZ))
1851         ISZIZ=ISZ(IZ)
1852         IF (MZIZ.GT.AMGVSS) THEN
1853         WID=(ZMIXSS(4,IZ)*CTHW+ZMIXSS(3,IZ)*STHW)**2*MZIZ**5/
1854      ,       48./PI/(AMGVSS*AMPL)**2
1855         CALL SSSAVE(ISZIZ,WID,ISGRAV,IDGM,0,0,0)
1856 C          Dalitz decay
1857           IF (MZIZ.GT.(AMGVSS+2*AME)) THEN
1858             WIDEE=WID*2*ALFAEM/(3*PI)*LOG(MZIZ/AME)
1859             CALL SSSAVE(ISZIZ,WIDEE,ISGRAV,IDE,-IDE,0,0)
1860           END IF
1861         END IF
1862         IF (MZIZ.GT.(AMZ+AMGVSS)) THEN
1863           WID=(2*(ZMIXSS(4,IZ)*STHW-ZMIXSS(3,IZ)*CTHW)**2+
1864      ,        (ZMIXSS(1,IZ)*SINBE-ZMIXSS(2,IZ)*COSBE)**2)*
1865      ,        (MZIZ**2-AMZ**2)**4/96./PI/MZIZ**3/(AMGVSS*AMPL)**2
1866           CALL SSSAVE(ISZIZ,WID,ISGRAV,IDZ,0,0,0)
1867         END IF
1868         IF (MZIZ.GT.(AMHL+AMGVSS)) THEN
1869           WID=(ZMIXSS(1,IZ)*COSA+ZMIXSS(2,IZ)*SINA)**2/6./
1870      ,        (AMGVSS*AMPL)**2*(MZIZ**2-AMHL**2)**4/16./PI/MZIZ**3
1871           CALL SSSAVE(ISZIZ,WID,ISGRAV,ISHL,0,0,0)
1872         END IF
1873         IF (MZIZ.GT.(AMHH+AMGVSS)) THEN
1874           WID=(-ZMIXSS(1,IZ)*SINA+ZMIXSS(2,IZ)*COSA)**2/6./
1875      ,        (AMGVSS*AMPL)**2*(MZIZ**2-AMHH**2)**4/16./PI/MZIZ**3
1876           CALL SSSAVE(ISZIZ,WID,ISGRAV,ISHH,0,0,0)
1877         END IF
1878         IF (MZIZ.GT.(AMHA+AMGVSS)) THEN
1879           WID=(ZMIXSS(1,IZ)*COSBE+ZMIXSS(2,IZ)*SINBE)**2/6./
1880      ,        (AMGVSS*AMPL)**2*(MZIZ**2-AMHA**2)**4/16./PI/MZIZ**3
1881           CALL SSSAVE(ISZIZ,WID,ISGRAV,ISHA,0,0,0)
1882         END IF
1883 250   CONTINUE
1884 C
1885 C          Normalize zi branching ratios
1886 C
1887       CALL SSNORM(ISZ1)
1888       CALL SSNORM(ISZ2)
1889       CALL SSNORM(ISZ3)
1890       CALL SSNORM(ISZ4)
1891 C
1892       RETURN
1893       END