ISAJET/isasusy/sswibf.F

   1 #include "isajet/pilot.h"
   2         SUBROUTINE SSWIBF
   3 C-----------------------------------------------------------------------
   4 C       This subroutine calculates the chargino (wi) 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),THJZ
  37 C
  38       REAL MWIW,SL,PP,SP,PL,PH,SH,COSB,SINB,COSA,SINA
  39      $,FACTOR,MZJZ,ULIM,XINTGL,COSBE,SINBE
  40      $,UPPER,CONST
  41       REAL W21ZU,W21ZN,W21ZL,W21ZD,W21U,W21D,W21S,W21C,W21N1
  42      $,W21N2,W21N3,W21E,W21M,W21L,STHW,CTHW
  43       REAL T3,XI2,CC,PSIINT,T2,T1,PHIINT,XI1,EF,A,Z,B
  44      $,TANB,FB,FT,SR2,G,PI,GP,FL
  45      $,MW2,SNW1,MW1,YM,BE,SNW2,XM,THX,THY
  46      $,BTN,APD,APL,APU,BTD,BTL,APN,BTU,Y,MZ1,FPI
  47       REAL TANW,COTW,XWINO,YWINO,SNIW,SNJZ
  48       REAL SSXINT,SSXLAM
  49       REAL WID,TERM1,TERM2,TERM3,TERM4,E,TERMH
  50       REAL FACT,ALJZ1,ALJZ2,BEJZ1,BEJZ2
  51       REAL ALIW1,ALIW2,AHCJZ,BHCJZ,TERMW,TERMN,TERM12,TERMN1,
  52      $TERMN2,TERMWN,TERMW1,TERMW2,TERMH1,TERMH2,TERMHN
  53       REAL XIPM,YIPM,COSL,SINL,BPWI(2),BPLWI(2)
  54       REAL BWI(2),AS,COST,SINT
  55       REAL POLNL,POLNR,POL1L,POL1R,POL2L,POL2R,POL12L,POL12R,
  56      $POLN1L,POLN1R,POLN2L,POLN2R
  57       REAL ASMB,MBMB,MBQ,ASMT,MTMT,MTQ,SUALFS
  58       REAL SUALFE,MTAMTA,MTAMB,MTAMZ,AMPL,AMPI
  59       REAL FUDGE
  60       DOUBLE PRECISION SSMQCD
  61       COMPLEX ZI,ZONE,Z1(2),Z2(2)
  62       INTEGER IW,JZ,IZ,ISZJZ,ISWIW
  63       DATA FUDGE/1.0/,AMPI/.140/
  64       DATA ZONE/(1.,0.)/,ZI/(0.,1.)/,FPI/.1315/
  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       MZ1=ABS(AMZ1SS)
 122       XWINO=.5*(THX*SIN(GAMMAL)*COS(GAMMAL)
 123      $-THY*SIN(GAMMAR)*COS(GAMMAR))
 124       YWINO=.5*(THX*SIN(GAMMAL)*COS(GAMMAL)
 125      $+THY*SIN(GAMMAR)*COS(GAMMAR))
 126       COST=COS(THETAT)
 127       SINT=SIN(THETAT)
 128       COSB=COS(THETAB)
 129       SINB=SIN(THETAB)
 130       COSL=COS(THETAL)
 131       SINL=SIN(THETAL)
 132 C
 133 C          Constants from Higgs mass matrix
 134 C
 135       SINA=SIN(ALFAH)
 136       COSA=COS(ALFAH)
 137 C
 138 C          Gaugino couplings
 139 C
 140       DO 100 IZ=1,4
 141         AUI(IZ)=G/SR2*ZMIXSS(3,IZ)+GP/3./SR2*ZMIXSS(4,IZ)
 142         BUI(IZ)=4.*GP/3./SR2*ZMIXSS(4,IZ)
 143         ADI(IZ)=-G/SR2*ZMIXSS(3,IZ)+GP/3./SR2*ZMIXSS(4,IZ)
 144         BDI(IZ)=-2.*GP/3./SR2*ZMIXSS(4,IZ)
 145         ALI(IZ)=G/SR2*ZMIXSS(3,IZ)+GP/SR2*ZMIXSS(4,IZ)
 146         BLI(IZ)=-SR2*GP*ZMIXSS(4,IZ)
 147         ANI(IZ)=G/SR2*ZMIXSS(3,IZ)-GP/SR2*ZMIXSS(4,IZ)
 148         BNI(IZ)=0.0
 149 100   CONTINUE
 150 C
 151       DO 110 IZ=1,4
 152       DO 110 JZ=1,4
 153         IF(IZ.LT.JZ) THEN
 154           WIJ(IZ,JZ)=SQRT(G**2+GP**2)
 155      $    *(ZMIXSS(1,IZ)*ZMIXSS(1,JZ)-ZMIXSS(2,IZ)*ZMIXSS(2,JZ))/4.
 156         ELSEIF(IZ.GT.JZ) THEN
 157           WIJ(IZ,JZ)=-SQRT(G**2+GP**2)
 158      $    *(ZMIXSS(1,IZ)*ZMIXSS(1,JZ)-ZMIXSS(2,IZ)*ZMIXSS(2,JZ))/4.
 159         ENDIF
 160 110   CONTINUE
 161 C
 162       AUWI(2)=G*THX*COS(GAMMAL)
 163       ADWI(2)=SNW2*G*THY*COS(GAMMAR)
 164       ALWI(2)=ADWI(2)
 165       ANWI(2)=AUWI(2)
 166       AUWI(1)=G*SIN(GAMMAL)
 167       ADWI(1)=SNW1*G*SIN(GAMMAR)
 168       ALWI(1)=ADWI(1)
 169       ANWI(1)=AUWI(1)
 170 C
 171       DO 120 IZ=1,4
 172         XIM(IZ)=.5*(SNW1*SIGN(1.,AMZISS(IZ))*(COS(GAMMAR)
 173      $  *ZMIXSS(1,IZ)/SR2+SIN(GAMMAR)*ZMIXSS(3,IZ))-COS(GAMMAL)
 174      $  *ZMIXSS(2,IZ)/SR2+SIN(GAMMAL)*ZMIXSS(3,IZ))
 175         YIM(IZ)=.5*(-SNW1*SIGN(1.,AMZISS(IZ))*(COS(GAMMAR)
 176      $  *ZMIXSS(1,IZ)/SR2+SIN(GAMMAR)*ZMIXSS(3,IZ))-COS(GAMMAL)
 177      $  *ZMIXSS(2,IZ)/SR2+SIN(GAMMAL)*ZMIXSS(3,IZ))
 178         XIP(IZ)=.5*(SNW2*SIGN(1.,AMZISS(IZ))*THY*(-SIN(GAMMAR)
 179      $  *ZMIXSS(1,IZ)/SR2+COS(GAMMAR)*ZMIXSS(3,IZ))+THX*(SIN(GAMMAL)
 180      $  *ZMIXSS(2,IZ)/SR2+COS(GAMMAL)*ZMIXSS(3,IZ)))
 181         YIP(IZ)=.5*(-SNW2*SIGN(1.,AMZISS(IZ))*THY*(-SIN(GAMMAR)
 182      $  *ZMIXSS(1,IZ)/SR2+COS(GAMMAR)*ZMIXSS(3,IZ))+THX*(SIN(GAMMAL)
 183      $  *ZMIXSS(2,IZ)/SR2+COS(GAMMAL)*ZMIXSS(3,IZ)))
 184 120   CONTINUE
 185 C
 186       DO 130 IZ=1,4
 187       DO 130 JZ=1,4
 188         IF(IZ.NE.JZ) THEN
 189           SNIJ(IZ,JZ)=-1.*SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ))
 190           XLIJ(IZ,JZ)=-SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ))
 191      $    *(ZMIXSS(2,IZ)*SINA-ZMIXSS(1,IZ)*COSA)
 192      $    *(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/2.
 193           XHIJ(IZ,JZ)=-SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ))
 194      $    *(ZMIXSS(2,IZ)*COSA+ZMIXSS(1,IZ)*SINA)
 195      $    *(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/2.
 196           XPIJ(IZ,JZ)=SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ))
 197      $    *(ZMIXSS(2,IZ)*SINBE-ZMIXSS(1,IZ)*COSBE)
 198      $    *(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/2.
 199           HIJ(IZ,JZ)=-SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ))
 200      $    *(ZMIXSS(2,IZ)*COSA+ZMIXSS(1,IZ)*SINA)
 201      $    *(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/2.
 202         ENDIF
 203 130   CONTINUE
 204 C
 205       SP=-.5*(-THY*SNW2*COSBE*SIN(GAMMAL)*SIN(GAMMAR)+
 206      $THY*SNW2*SINBE*COS(GAMMAL)*COS(GAMMAR)-
 207      $THX*SNW1*COSBE*COS(GAMMAL)*COS(GAMMAR)+
 208      $THX*SNW1*SINBE*SIN(GAMMAL)*SIN(GAMMAR))
 209       PP=-.5*(-THY*SNW2*COSBE*SIN(GAMMAL)*SIN(GAMMAR)+
 210      $THY*SNW2*SINBE*COS(GAMMAL)*COS(GAMMAR)+
 211      $THX*SNW1*COSBE*COS(GAMMAL)*COS(GAMMAR)-
 212      $THX*SNW1*SINBE*SIN(GAMMAL)*SIN(GAMMAR))
 213 C
 214       SL=.5*(THY*SNW2*SINA*COS(GAMMAL)*COS(GAMMAR)-
 215      $THY*SNW2*COSA*SIN(GAMMAL)*SIN(GAMMAR)+
 216      $THX*SNW1*COSA*COS(GAMMAL)*COS(GAMMAR)-
 217      $THX*SNW1*SINA*SIN(GAMMAL)*SIN(GAMMAR))
 218       PL=.5*(THY*SNW2*SINA*COS(GAMMAL)*COS(GAMMAR)-
 219      $THY*SNW2*COSA*SIN(GAMMAL)*SIN(GAMMAR)-
 220      $THX*SNW1*COSA*COS(GAMMAL)*COS(GAMMAR)+
 221      $THX*SNW1*SINA*SIN(GAMMAL)*SIN(GAMMAR))
 222 C
 223       SH=.5*(THY*SNW2*COSA*COS(GAMMAL)*COS(GAMMAR)+
 224      $THY*SNW2*SINA*SIN(GAMMAL)*SIN(GAMMAR)-
 225      $THX*SNW1*SINA*COS(GAMMAL)*COS(GAMMAR)-
 226      $THX*SNW1*COSA*SIN(GAMMAL)*SIN(GAMMAR))
 227       PH=.5*(THY*SNW2*COSA*COS(GAMMAL)*COS(GAMMAR)+
 228      $THY*SNW2*SINA*SIN(GAMMAL)*SIN(GAMMAR)+
 229      $THX*SNW1*SINA*COS(GAMMAL)*COS(GAMMAR)+
 230      $THX*SNW1*COSA*SIN(GAMMAL)*SIN(GAMMAR))
 231 C
 232       DO 140 IZ=1,4
 233         V1I(IZ)=-SIN(GAMMAR)/SR2*(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))
 234      $  -G*COS(GAMMAR)*ZMIXSS(1,IZ)
 235         V2I(IZ)=COS(GAMMAR)/SR2*(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))
 236      $  -G*SIN(GAMMAR)*ZMIXSS(1,IZ)
 237         V3I(IZ)=-SIN(GAMMAL)/SR2*(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))
 238      $  +G*COS(GAMMAL)*ZMIXSS(2,IZ)
 239         V4I(IZ)=COS(GAMMAL)/SR2*(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))
 240      $  +G*SIN(GAMMAL)*ZMIXSS(2,IZ)
 241 140   CONTINUE
 242 C
 243 C-----------------------------------------------------------------------
 244 C          Generate Chargino Branching Fractions
 245 C-----------------------------------------------------------------------
 246       ISZ(1)=ISZ1
 247       ISZ(2)=ISZ2
 248       ISZ(3)=ISZ3
 249       ISZ(4)=ISZ4
 250 C     FIRST TRY EXCLUSIVE DECAY TO SINGLE PION
 251       IF (MW1.GT.(MZ1+AMPI).AND.MW1.LT.(MZ1+1.)) THEN
 252         WID=G**4*FPI**2*SQRT(SSXLAM(MW1**2,MZ1**2,AMPI**2))/
 253      $      128./MW1**3/PI/AMW**4*((XIM(1)**2+YIM(1)**2)*
 254      $      (MW1-MZ1)**2*(MW1+MZ1)**2-AMPI**2*(XIM(1)**2*(MW1-MZ1)**2+
 255      $      YIM(1)**2*(MW1+MZ1)**2))
 256         CALL SSSAVE(ISW1,WID,ISZ(1),120,0,0,0)
 257       END IF
 258       DO 300 IW=1,2
 259 C          Loop over w1, w2
 260         IF(IW.EQ.1) THEN
 261             MWIW=MW1
 262             SNIW=SNW1
 263             ISWIW=ISW1
 264         ELSE
 265             MWIW=MW2
 266             SNIW=SNW2
 267             ISWIW=ISW2
 268         ENDIF
 269 C
 270 C          Decays to zj
 271 C
 272         DO 310 JZ=1,4
 273           MZJZ=ABS(AMZISS(JZ))
 274           ISZJZ=ISZ(JZ)
 275           SNJZ=SIGN(1.,AMZISS(JZ))
 276           THJZ=0.
 277           IF (AMZISS(JZ).LT.0.) THJZ=1.
 278           IF(MWIW.LE.FUDGE*MZJZ) GOTO 310
 279 C          Couplings
 280           IF(IW.EQ.1) THEN
 281             XIPM=XIM(JZ)
 282             YIPM=YIM(JZ)
 283           ELSE
 284             XIPM=XIP(JZ)
 285             YIPM=YIP(JZ)
 286           ENDIF
 287 C
 288 C          wi --> f + fbar + zj
 289 C
 290           IF (MWIW.GT.(MZJZ+AMUP+AMDN)) THEN
 291           IF (IW.EQ.1.AND.JZ.EQ.1.AND.(MW1.LT.(MZ1+1.))) THEN
 292 C           EXIT HADRONIC MODE IF EXCLUSIVE PION DECAY ACTIVE
 293             GO TO 200
 294           END IF
 295             IF (MWIW.LT.(AMW+MZJZ)) THEN
 296             TMP(1)=XIPM**2+YIPM**2
 297             TMP(2)=XIPM**2-YIPM**2
 298             TMP(3)=MWIW
 299             TMP(4)=MZJZ
 300             TERM1=SSXINT(MZJZ,SSWZF1,(MWIW**2+MZJZ**2)/2./MWIW)/
 301      $            2./MWIW/(2*PI)**5*2*G**4*PI**2/3.
 302             ELSE
 303             TERM1=0.
 304             END IF
 305             ULIM=MWIW/2.*(1.-MZJZ**2/MWIW**2)
 306             TMP(1)=MWIW
 307             TMP(3)=MZJZ
 308             IF (MWIW.LT.AMULSS) THEN
 309             TMP(2)=AMULSS
 310             PSIINT=SSXINT(0.,SSWZF2,ULIM)
 311             T1=AUI(JZ)**2*ADWI(IW)**2*PSIINT
 312             ELSE
 313             T1=0.
 314             END IF
 315             IF (MWIW.LT.AMDLSS) THEN
 316             TMP(2)=AMDLSS
 317             PSIINT=SSXINT(0.,SSWZF2,ULIM)
 318             T2=ADI(JZ)**2*AUWI(IW)**2*PSIINT
 319             ELSE
 320             T2=0.
 321             END IF
 322             IF (MWIW.LT.AMDLSS.AND.MWIW.LT.AMULSS) THEN
 323             TMP(2)=0.
 324             TMP(4)=AMDLSS
 325             TMP(5)=AMULSS
 326             PHIINT=SSXINT(0.,SSGX2,ULIM)
 327             CC=2*SIGN(1.,AMZISS(JZ))*AUWI(IW)*ADWI(IW)*AUI(JZ)*ADI(JZ)
 328             T3=CC*PHIINT
 329             ELSE
 330             T3=0.
 331             END IF
 332             TERM2=(T1+T2+T3)/2./MWIW/(2*PI)**5
 333             FACTOR=1./2./MWIW/(2*PI)**5*2*SR2*G**2
 334             TMP(1)=MWIW
 335             TMP(3)=MZJZ
 336             IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMULSS) THEN
 337             TMP(2)=AMULSS
 338             XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
 339             XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
 340             TERM3=FACTOR*ADWI(IW)*AUI(JZ)*((XIPM-YIPM)*XI1
 341      $      -(XIPM+YIPM)*XI2)*SIGN(1.,AMZISS(JZ))
 342             ELSE
 343             TERM3=0.
 344             END IF
 345             IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMDLSS) THEN
 346             TMP(2)=AMDLSS
 347             XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
 348             XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
 349             TERM4=-FACTOR*AUWI(IW)*ADI(JZ)*((XIPM+YIPM)*XI1
 350      $      -(XIPM-YIPM)*XI2)
 351             ELSE
 352             TERM4=0.
 353             END IF
 354             WID=3*(TERM1+TERM2+TERM3+TERM4)
 355             CALL SSSAVE(ISWIW,WID,ISZJZ,IDUP,-IDDN,0,0)
 356 C         Enter information for decay matrix element
 357             Z1(1)=ZI**THJZ*G*XIPM
 358             Z1(2)=ZI**THJZ*G*YIPM
 359             Z2(1)=G/2./SR2
 360             Z2(2)=-G/2./SR2
 361             CALL SSME3(1,AMW,Z1,Z2)
 362             Z1(1)=ZI*AUWI(IW)/2.
 363             Z1(2)=Z1(1)
 364             Z2(1)=-CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ADI(JZ))/2.
 365             Z2(2)=-Z2(1)
 366             CALL SSME3(2,AMDLSS,Z1,Z2)
 367             Z1(1)=CONJG(ZI*ADWI(IW))/2.
 368             Z1(2)=-Z1(1)
 369             Z2(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*AUI(JZ)/2.
 370             Z2(2)=Z2(1)
 371             CALL SSME3(3,AMULSS,Z1,Z2)
 372           END IF
 373 C          wi --> c + sbar + zj
 374 200       IF (MWIW.GT.(MZJZ+AMCH+AMST)) THEN
 375             IF (MWIW.LT.(AMW+MZJZ)) THEN
 376             TMP(1)=XIPM**2+YIPM**2
 377             TMP(2)=XIPM**2-YIPM**2
 378             TMP(3)=MWIW
 379             TMP(4)=MZJZ
 380             TERM1=SSXINT(MZJZ,SSWZF1,(MWIW**2+MZJZ**2)/2./MWIW)/
 381      $            2./MWIW/(2*PI)**5*2*G**4*PI**2/3.
 382             ELSE
 383             TERM1=0.
 384             END IF
 385             ULIM=MWIW/2.*(1.-MZJZ**2/MWIW**2)
 386             TMP(1)=MWIW
 387             TMP(3)=MZJZ
 388             IF (MWIW.LT.AMCLSS) THEN
 389             TMP(2)=AMCLSS
 390             PSIINT=SSXINT(0.,SSWZF2,ULIM)
 391             T1=AUI(JZ)**2*ADWI(IW)**2*PSIINT
 392             ELSE
 393             T1=0.
 394             END IF
 395             IF (MWIW.LT.AMSLSS) THEN
 396             TMP(2)=AMSLSS
 397             PSIINT=SSXINT(0.,SSWZF2,ULIM)
 398             T2=ADI(JZ)**2*AUWI(IW)**2*PSIINT
 399             ELSE
 400             T2=0.
 401             END IF
 402             IF (MWIW.LT.AMSLSS.AND.MWIW.LT.AMCLSS) THEN
 403             TMP(2)=0.
 404             TMP(4)=AMSLSS
 405             TMP(5)=AMCLSS
 406             PHIINT=SSXINT(0.,SSGX2,ULIM)
 407             CC=2*SIGN(1.,AMZISS(JZ))*AUWI(IW)*ADWI(IW)*AUI(JZ)*ADI(JZ)
 408             T3=CC*PHIINT
 409             ELSE
 410             T3=0.
 411             END IF
 412             TERM2=(T1+T2+T3)/2./MWIW/(2*PI)**5
 413             FACTOR=1./2./MWIW/(2*PI)**5*2*SR2*G**2
 414             TMP(1)=MWIW
 415             TMP(3)=MZJZ
 416             IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMCLSS) THEN
 417             TMP(2)=AMCLSS
 418             XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
 419             XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
 420             TERM3=FACTOR*ADWI(IW)*AUI(JZ)*((XIPM-YIPM)*XI1
 421      $      -(XIPM+YIPM)*XI2)*SIGN(1.,AMZISS(JZ))
 422             ELSE
 423             TERM3=0.
 424             END IF
 425             IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMSLSS) THEN
 426             TMP(2)=AMSLSS
 427             XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
 428             XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
 429             TERM4=-FACTOR*AUWI(IW)*ADI(JZ)*((XIPM+YIPM)*XI1
 430      $      -(XIPM-YIPM)*XI2)
 431             ELSE
 432             TERM4=0.
 433             END IF
 434             WID=3*(TERM1+TERM2+TERM3+TERM4)
 435             CALL SSSAVE(ISWIW,WID,ISZJZ,IDCH,-IDST,0,0)
 436 C         Enter information for decay matrix element
 437             Z1(1)=ZI**THJZ*G*XIPM
 438             Z1(2)=ZI**THJZ*G*YIPM
 439             Z2(1)=G/2./SR2
 440             Z2(2)=-G/2./SR2
 441             CALL SSME3(1,AMW,Z1,Z2)
 442             Z1(1)=ZI*AUWI(IW)/2.
 443             Z1(2)=Z1(1)
 444             Z2(1)=-CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ADI(JZ))/2.
 445             Z2(2)=-Z2(1)
 446             CALL SSME3(2,AMSLSS,Z1,Z2)
 447             Z1(1)=CONJG(ZI*ADWI(IW))/2.
 448             Z1(2)=-Z1(1)
 449             Z2(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*AUI(JZ)/2.
 450             Z2(2)=Z2(1)
 451             CALL SSME3(3,AMCLSS,Z1,Z2)
 452           END IF
 453 C          wi -> t + bbar + zj neglected since 2-body modes should dominate
 454 C          wi --> nu_e + e + zj
 455           IF (MWIW.GT.(MZJZ+AME)) THEN
 456             IF (MWIW.LT.(AMW+MZJZ)) THEN
 457             TMP(1)=XIPM**2+YIPM**2
 458             TMP(2)=XIPM**2-YIPM**2
 459             TMP(3)=MWIW
 460             TMP(4)=MZJZ
 461             TERM1=SSXINT(MZJZ,SSWZF1,(MWIW**2+MZJZ**2)/2./MWIW)/
 462      $            2./MWIW/(2*PI)**5*2*G**4*PI**2/3.
 463             ELSE
 464             TERM1=0.
 465             END IF
 466             ULIM=MWIW/2.*(1.-MZJZ**2/MWIW**2)
 467             TMP(1)=MWIW
 468             TMP(3)=MZJZ
 469             IF (MWIW.LT.AMN1SS) THEN
 470             TMP(2)=AMN1SS
 471             T1=ANI(JZ)**2*ALWI(IW)**2*SSXINT(0.,SSWZF2,ULIM)
 472             ELSE
 473             T1=0.
 474             END IF
 475             IF (MWIW.LT.AMELSS) THEN
 476             TMP(2)=AMELSS
 477             T2=ALI(JZ)**2*ANWI(IW)**2*SSXINT(0.,SSWZF2,ULIM)
 478             ELSE
 479             T2=0.
 480             END IF
 481             IF (MWIW.LT.AMELSS.AND.MWIW.LT.AMN1SS) THEN
 482             TMP(2)=0.
 483             TMP(4)=AMELSS
 484             TMP(5)=AMN1SS
 485             T3=-2*SIGN(1.,AMZISS(JZ))*ANWI(IW)*ALWI(IW)*ANI(JZ)*
 486      $         ALI(JZ)*SSXINT(0.,SSGX2,ULIM)
 487             ELSE
 488             T3=0.
 489             END IF
 490             TERM2=(T1+T2+T3)/2./MWIW/(2*PI)**5
 491             FACTOR=1./2./MWIW/(2*PI)**5*2*SR2*G**2
 492             TMP(1)=MWIW
 493             TMP(3)=MZJZ
 494             IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMN1SS) THEN
 495             TMP(2)=AMN1SS
 496             XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
 497             XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
 498             TERM3=FACTOR*ALWI(IW)*ANI(JZ)*((XIPM-YIPM)*XI1
 499      $      -(XIPM+YIPM)*XI2)*SIGN(1.,AMZISS(JZ))
 500             ELSE
 501             TERM3=0.
 502             END IF
 503             IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMELSS) THEN
 504             TMP(2)=AMELSS
 505             XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
 506             XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
 507             TERM4=FACTOR*ANWI(IW)*ALI(JZ)*((XIPM+YIPM)*XI1
 508      $      -(XIPM-YIPM)*XI2)
 509             ELSE
 510             TERM4=0.
 511             END IF
 512             WID=TERM1+TERM2+TERM3+TERM4
 513             CALL SSSAVE(ISWIW,WID,ISZJZ,-IDE,IDNE,0,0)
 514 C         Enter information for decay matrix element
 515             Z1(1)=ZI**THJZ*G*XIPM
 516             Z1(2)=ZI**THJZ*G*YIPM
 517             Z2(1)=G/2./SR2
 518             Z2(2)=-G/2./SR2
 519             CALL SSME3(1,AMW,Z1,Z2)
 520             Z1(1)=ZI*ANWI(IW)/2.
 521             Z1(2)=Z1(1)
 522             Z2(1)=-CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ALI(JZ))/2.
 523             Z2(2)=-Z2(1)
 524             CALL SSME3(2,AMELSS,Z1,Z2)
 525             Z1(1)=CONJG(ZI*ALWI(IW))/2.
 526             Z1(2)=-Z1(1)
 527             Z2(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2.
 528             Z2(2)=Z2(1)
 529             CALL SSME3(3,AMN1SS,Z1,Z2)
 530           END IF
 531 C          wi --> nu_mu + mu + zj
 532           IF (MWIW.GT.(MZJZ+AMMU)) THEN
 533             IF (MWIW.LT.(AMW+MZJZ)) THEN
 534             TMP(1)=XIPM**2+YIPM**2
 535             TMP(2)=XIPM**2-YIPM**2
 536             TMP(3)=MWIW
 537             TMP(4)=MZJZ
 538             TERM1=SSXINT(MZJZ,SSWZF1,(MWIW**2+MZJZ**2)/2./MWIW)/
 539      $            2./MWIW/(2*PI)**5*2*G**4*PI**2/3.
 540             ELSE
 541             TERM1=0.
 542             END IF
 543             ULIM=MWIW/2.*(1.-MZJZ**2/MWIW**2)
 544             TMP(1)=MWIW
 545             TMP(3)=MZJZ
 546             IF (MWIW.LT.AMN2SS) THEN
 547             TMP(2)=AMN2SS
 548             T1=ANI(JZ)**2*ALWI(IW)**2*SSXINT(0.,SSWZF2,ULIM)
 549             ELSE
 550             T1=0.
 551             END IF
 552             IF (MWIW.LT.AMMLSS) THEN
 553             TMP(2)=AMMLSS
 554             T2=ALI(JZ)**2*ANWI(IW)**2*SSXINT(0.,SSWZF2,ULIM)
 555             ELSE
 556             T2=0.
 557             END IF
 558             IF (MWIW.LT.AMMLSS.AND.MWIW.LT.AMN2SS) THEN
 559             TMP(2)=0.
 560             TMP(4)=AMMLSS
 561             TMP(5)=AMN2SS
 562             T3=-2*SIGN(1.,AMZISS(JZ))*ANWI(IW)*ALWI(IW)*ANI(JZ)*
 563      $         ALI(JZ)*SSXINT(0.,SSGX2,ULIM)
 564             ELSE
 565             T3=0.
 566             END IF
 567             TERM2=(T1+T2+T3)/2./MWIW/(2*PI)**5
 568             FACTOR=1./2./MWIW/(2*PI)**5*2*SR2*G**2
 569             TMP(1)=MWIW
 570             TMP(3)=MZJZ
 571             IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMN2SS) THEN
 572             TMP(2)=AMN2SS
 573             XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
 574             XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
 575             TERM3=FACTOR*ALWI(IW)*ANI(JZ)*((XIPM-YIPM)*XI1
 576      $      -(XIPM+YIPM)*XI2)*SIGN(1.,AMZISS(JZ))
 577             ELSE
 578             TERM3=0.
 579             END IF
 580             IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMMLSS) THEN
 581             TMP(2)=AMMLSS
 582             XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
 583             XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
 584             TERM4=FACTOR*ANWI(IW)*ALI(JZ)*((XIPM+YIPM)*XI1
 585      $      -(XIPM-YIPM)*XI2)
 586             ELSE
 587             TERM4=0.
 588             END IF
 589             WID=TERM1+TERM2+TERM3+TERM4
 590             CALL SSSAVE(ISWIW,WID,ISZJZ,-IDMU,IDNM,0,0)
 591 C         Enter information for decay matrix element
 592             Z1(1)=ZI**THJZ*G*XIPM
 593             Z1(2)=ZI**THJZ*G*YIPM
 594             Z2(1)=G/2./SR2
 595             Z2(2)=-G/2./SR2
 596             CALL SSME3(1,AMW,Z1,Z2)
 597             Z1(1)=ZI*ANWI(IW)/2.
 598             Z1(2)=Z1(1)
 599             Z2(1)=-CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ALI(JZ))/2.
 600             Z2(2)=-Z2(1)
 601             CALL SSME3(2,AMMLSS,Z1,Z2)
 602             Z1(1)=CONJG(ZI*ALWI(IW))/2.
 603             Z1(2)=-Z1(1)
 604             Z2(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2.
 605             Z2(2)=Z2(1)
 606             CALL SSME3(3,AMN2SS,Z1,Z2)
 607           END IF
 608 C          wi --> nu_tau + tau + zj ; includes mixing and Yukawas
 609           FACT=1./2./MWIW/(2*PI)**5
 610           ALJZ1=-ALI(JZ)*COSL-FL*ZMIXSS(2,JZ)*SINL
 611           ALJZ2=-ALI(JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
 612           BEJZ1=BLI(JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
 613           BEJZ2=-BLI(JZ)*COSL+FL*ZMIXSS(2,JZ)*SINL
 614           SNJZ=SIGN(1.,AMZISS(JZ))
 615 C         Change ALWI def'ns in accord with Drees note
 616           ALWI(1)=-G*SIN(GAMMAR)
 617           ALWI(2)=-G*THY*COS(GAMMAR)
 618 C          Polarization for stau_1 -> z1ss+tau.
 619 C          See above for other cases.
 620           IF(IW.EQ.1.AND.JZ.EQ.1) THEN
 621             PTAU1(JZ)=(BEJZ1**2-ALJZ1**2)/(BEJZ1**2+ALJZ1**2)
 622             PTAU2(JZ)=(BEJZ2**2-ALJZ2**2)/(BEJZ2**2+ALJZ2**2)
 623           ENDIF
 624           IF (IW.EQ.1) THEN
 625             ALIW1=-G*SIN(GAMMAL)*COSL+FL*COS(GAMMAL)*SINL
 626             ALIW2=-G*SIN(GAMMAL)*SINL-FL*COS(GAMMAL)*COSL
 627             AHCJZ=COSBE*V2I(JZ)
 628             BHCJZ=-SINBE*V4I(JZ)
 629           ELSE IF (IW.EQ.2) THEN
 630             ALIW1=(-G*COS(GAMMAL)*COSL-FL*SIN(GAMMAL)*SINL)*THX
 631             ALIW2=(-G*COS(GAMMAL)*SINL+FL*SIN(GAMMAL)*COSL)*THX
 632             AHCJZ=COSBE*V1I(JZ)*THY
 633             BHCJZ=-SINBE*V3I(JZ)*THX
 634           END IF
 635           IF (MWIW.GT.(MZJZ+AMTAU)) THEN
 636             IF (MWIW.LT.(AMW+MZJZ)) THEN
 637             TMP(1)=XIPM**2+YIPM**2
 638             TMP(2)=XIPM**2-YIPM**2
 639             TMP(3)=MWIW
 640             TMP(4)=MZJZ
 641             TERMW=SSXINT(MZJZ,SSWZF1,(MWIW**2+MZJZ**2)/2./MWIW)*
 642      $            FACT*2*G**4*PI**2/3.
 643             ELSE
 644             TERMW=0.
 645             END IF
 646             ULIM=(MWIW**2-MZJZ**2)/2./MWIW
 647             TMP(1)=MWIW
 648             TMP(3)=MZJZ
 649             IF (MWIW.LT.AMN3SS) THEN
 650             TMP(2)=AMN3SS
 651             POLNL=FACT*ANI(JZ)**2*ALWI(IW)**2*
 652      $            SSXINT(0.,SSWZF2,ULIM)
 653             POLNR=POLNL*BPLWI(IW)**2/ALWI(IW)**2
 654             TERMN=POLNL+POLNR
 655             ELSE
 656             POLNL=0.
 657             POLNR=0.
 658             TERMN=0.
 659             END IF
 660             IF (MWIW.LT.AML1SS) THEN
 661             TMP(2)=AML1SS
 662             POL1L=FACT*ALJZ1**2*ALIW1**2*SSXINT(0.,SSWZF2,ULIM)
 663             POL1R=POL1L*BEJZ1**2/ALJZ1**2
 664             TERM1=POL1L+POL1R
 665             ELSE
 666             POL1L=0.
 667             POL1R=0.
 668             TERM1=0.
 669             END IF
 670             IF (MWIW.LT.AML2SS) THEN
 671             TMP(2)=AML2SS
 672             POL2L=FACT*ALJZ2**2*ALIW2**2*SSXINT(0.,SSWZF2,ULIM)
 673             POL2R=POL2L*BEJZ2**2/ALJZ2**2
 674             TERM2=POL2L+POL2R
 675             ELSE
 676             POL2L=0.
 677             POL2R=0.
 678             TERM2=0.
 679             END IF
 680             IF (MWIW.LT.AML1SS) THEN
 681             TMP(2)=0.
 682             TMP(3)=MZJZ
 683             TMP(4)=AML1SS
 684             TMP(5)=AML2SS
 685             POL12L=FACT*2*ALIW1*ALIW2*ALJZ1*ALJZ2*
 686      $             SSXINT(0.,SSGX1,ULIM)
 687             POL12R=POL12L*BEJZ1*BEJZ2/ALJZ1/ALJZ2
 688             TERM12=POL12L+POL12R
 689             ELSE
 690             POL12L=0.
 691             POL12R=0.
 692             TERM12=0.
 693             END IF
 694             IF (MWIW.LT.(AMHC+MZJZ)) THEN
 695             TMP(2)=AMHC
 696             TMP(3)=MZJZ
 697             TMP(4)=AHCJZ
 698             TMP(5)=BHCJZ
 699             TMP(6)=SIGN(1.,AMZISS(JZ))*SIGN(1.,AMWISS(IW))
 700             TERMH=FACT*PI**2*MWIW*(G*MTAMZ*TANB/AMW)**2/2.*
 701      $         SSXINT(MZJZ,SSWZF6,(MWIW**2+MZJZ**2)/2./MWIW)
 702             ELSE
 703             TERMH=0.
 704             END IF
 705             IF (MWIW.LT.AML1SS.AND.MWIW.LT.AMN3SS) THEN
 706             TMP(2)=0.
 707             TMP(3)=MZJZ
 708             TMP(4)=AMN3SS
 709             TMP(5)=AML1SS
 710             POLN1L=+2*FACT*ANI(JZ)*ALIW1*SNJZ*SNIW*ALWI(IW)*
 711      $             ALJZ1*SSXINT(0.,SSGX2,ULIM)
 712             POLN1R=-2*FACT*ANI(JZ)*ALIW1*BPLWI(IW)*BEJZ1*
 713      $             SSXINT(0.,SSGX8,ULIM)
 714             TERMN1=POLN1L+POLN1R
 715             ELSE
 716             POLN1L=0.
 717             POLN1R=0.
 718             TERMN1=0.
 719             END IF
 720             IF (MWIW.LT.AML2SS.AND.MWIW.LT.AMN3SS) THEN
 721             TMP(2)=0.
 722             TMP(3)=MZJZ
 723             TMP(4)=AMN3SS
 724             TMP(5)=AML2SS
 725             POLN2L=+2*FACT*ANI(JZ)*ALIW2*SNJZ*SNIW*ALWI(IW)*ALJZ2*
 726      $             SSXINT(0.,SSGX2,ULIM)
 727             POLN2R=-2*FACT*ANI(JZ)*ALIW2*BPLWI(IW)*BEJZ2*
 728      $             SSXINT(0.,SSGX8,ULIM)
 729             TERMN2=POLN2L+POLN2R
 730             ELSE
 731             POLN2L=0.
 732             POLN2R=0.
 733             TERMN2=0.
 734             END IF
 735             TMP(1)=MWIW
 736             TMP(3)=MZJZ
 737             IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMN3SS) THEN
 738             TMP(2)=AMN3SS
 739             XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
 740             XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
 741             TERMWN=2*SR2*G**2*FACT*ALWI(IW)*ANI(JZ)*((XIPM-
 742      $         YIPM)*XI1-(XIPM+YIPM)*XI2)*SIGN(1.,AMZISS(JZ))
 743             ELSE
 744             TERMWN=0.
 745             END IF
 746             IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AML1SS) THEN
 747             TMP(2)=AML1SS
 748             XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
 749             XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
 750             TERMW1=2*SR2*G**2*FACT*ALIW1*ALJZ1*
 751      $             ((XIPM+YIPM)*XI1-(XIPM-YIPM)*XI2)
 752             ELSE
 753             TERMW1=0.
 754             END IF
 755             IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AML2SS) THEN
 756             TMP(2)=AML2SS
 757             XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
 758             XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
 759             TERMW2=2*SR2*G**2*FACT*ALIW2*ALJZ2*
 760      $             ((XIPM+YIPM)*XI1-(XIPM-YIPM)*XI2)
 761             ELSE
 762             TERMW2=0.
 763             END IF
 764             TMP(2)=MZJZ
 765             TMP(3)=AMHC
 766             TMP(5)=AHCJZ
 767             TMP(6)=BHCJZ
 768             TMP(7)=SNJZ*SNIW
 769             IF (MWIW.LT.(AMHC+MZJZ).AND.MWIW.LT.AML1SS) THEN
 770             TMP(4)=AML1SS
 771             TERMH1=PI**2/2./MWIW*FACT*SR2*ALIW1*BEJZ1*G*MTAMZ*
 772      $             TANB/AMW*SSXINT(0.,SSWZF7,(MWIW-MZJZ)**2)
 773             ELSE
 774             TERMH1=0.
 775             END IF
 776             IF (MWIW.LT.(AMHC+MZJZ).AND.MWIW.LT.AML2SS) THEN
 777             TMP(4)=AML2SS
 778             TERMH2=PI**2/2./MWIW*FACT*SR2*ALIW2*BEJZ2*G*MTAMZ*
 779      $             TANB/AMW*SSXINT(0.,SSWZF7,(MWIW-MZJZ)**2)
 780             ELSE
 781             TERMH2=0.
 782             END IF
 783             IF (MWIW.LT.(AMHC+MZJZ).AND.MWIW.LT.AMN3SS) THEN
 784             TMP(4)=AMN3SS
 785             TERMHN=PI**2/2./MWIW*FACT*SR2*ANI(JZ)*BPLWI(IW)*G*
 786      $         MTAMZ*TANB/AMW*SSXINT(0.,SSWZF7,(MWIW-MZJZ)**2)
 787             ELSE
 788             TERMHN=0.
 789             END IF
 790             WID=TERMW+TERMN+TERM1+TERM2+TERMH+TERMWN+TERMW1+
 791      $          TERMW2+TERM12+TERMN1+TERMN2+TERMH1+TERMH2+
 792      $          TERMHN
 793 C          tau polarization for 3-body w1 -> z1 tau nu
 794             IF (IW.EQ.1.AND.JZ.EQ.1.AND.WID.GT.0.) THEN
 795               PTAUWZ=(POLNR+POL1R+POL2R+POL12R+TERMH+POLN1R+
 796      $               POLN2R+TERMHN+TERMH1+TERMH2-(TERMW+POLNL+
 797      $               POL1L+POL2L+POL12L+POLN1L+POLN2L+TERMWN+
 798      $               TERMW1+TERMW2))/WID
 799             END IF
 800             CALL SSSAVE(ISWIW,WID,ISZJZ,-IDTAU,IDNT,0,0)
 801             Z1(1)=ZI**THJZ*G*XIPM
 802             Z1(2)=ZI**THJZ*G*YIPM
 803             Z2(1)=-G/2./SR2
 804             Z2(2)=-Z2(1)
 805             CALL SSME3(1,AMW,Z1,Z2)
 806             Z1(1)=-(ZI)**THJZ*ALIW1/2.
 807             Z1(2)=-Z1(1)
 808             Z2(1)=(ALJZ1+(-1.)**THJZ*BEJZ1)/2.
 809             Z2(2)=(ALJZ1-(-1.)**THJZ*BEJZ1)/2.
 810             CALL SSME3(2,AML1SS,Z1,Z2)
 811             Z1(1)=-(ZI)**THJZ*ALIW2/2.
 812             Z1(2)=-Z1(1)
 813             Z2(1)=(ALJZ2+(-1.)**THJZ*BEJZ2)/2.
 814             Z2(2)=(ALJZ2-(-1.)**THJZ*BEJZ2)/2.
 815             CALL SSME3(2,AML2SS,Z1,Z2)
 816             Z1(1)=(-1.)**THJZ*ANI(JZ)/2.
 817             Z1(2)=-Z1(1)
 818             Z2(1)=(SNIW*ALWI(IW)+BPLWI(IW))/2.
 819             Z2(2)=(SNIW*ALWI(IW)-BPLWI(IW))/2.
 820             CALL SSME3(3,AMN3SS,Z1,Z2)
 821             Z1(1)=ZI**THJZ*G*MTAMZ*TANB/SR2/AMW/2.
 822             Z1(2)=-Z1(1)
 823             Z2(1)=(SNIW*AHCJZ-SNJZ*BHCJZ)/2.
 824             Z2(2)=(SNIW*AHCJZ+SNJZ*BHCJZ)/2.
 825             CALL SSME3(4,AMHC,Z1,Z2)
 826           END IF
 827 C
 828 C          wi --> w + zj
 829 C
 830           IF (MWIW.GT.(MZJZ+AMW)) THEN
 831             EF=MWIW**2+MZJZ**2-AMW**2+((MWIW**2-MZJZ**2)**2
 832      $         -AMW**4)/AMW/AMW
 833             WID=G*G*SQRT(SSXLAM(MWIW**2,MZJZ**2,AMW**2))/32./PI/
 834      $          MWIW**3*(2.*EF*(XIPM**2+YIPM**2)-12*MWIW*MZJZ*
 835      $          (XIPM**2-YIPM**2))
 836             CALL SSSAVE(ISWIW,WID,ISZJZ,IDW,0,0,0)
 837           END IF
 838 C
 839 C          wi --> h+ + zj
 840 C
 841           IF (MWIW.GT.(MZJZ+AMHC)) THEN
 842             IF(IW.EQ.1) THEN
 843               A=(SNW1*COSBE*V2I(JZ)
 844      $        -SIGN(1.,AMZISS(JZ))*SINBE*V4I(JZ))/2.
 845               B=(SNW1*COSBE*V2I(JZ)
 846      $        +SIGN(1.,AMZISS(JZ))*SINBE*V4I(JZ))/2.
 847             ELSE
 848               A=(THY*SNW2*COSBE*V1I(JZ)
 849      $        -THX*SIGN(1.,AMZISS(JZ))*SINBE*V3I(JZ))/2.
 850               B=(THY*SNW2*COSBE*V1I(JZ)
 851      $        +THX*SIGN(1.,AMZISS(JZ))*SINBE*V3I(JZ))/2.
 852             ENDIF
 853             WID=SQRT(MWIW**4+MZJZ**4+AMHC**4-2.*(MWIW*MZJZ)**2-
 854      $      2.*(MWIW*AMHC)**2-2.*(MZJZ*AMHC)**2)/8./PI/MWIW**3*
 855      $      ((A*A+B*B)*(MWIW*MWIW+MZJZ*MZJZ-AMHC*AMHC)/2.
 856      $      +(A*A-B*B)*MWIW*MZJZ)
 857             CALL SSSAVE(ISWIW,WID,ISZJZ,ISHC,0,0,0)
 858           ENDIF
 859 310     CONTINUE
 860 C
 861 C          wi --> quark + squark
 862 C
 863         IF(MWIW.GT.(AMULSS+AMDN)) THEN
 864           WID=3.*ADWI(IW)**2/32./PI/MWIW*(1.+AMDN**2/MWIW**2-
 865      $     AMULSS**2/MWIW**2)*SQRT(SSXLAM(MWIW**2,AMDN**2,AMULSS**2))
 866           CALL SSSAVE(ISWIW,WID,+ISUPL,-IDDN,0,0,0)
 867         END IF
 868         IF(MWIW.GT.(AMDLSS+AMUP)) THEN
 869           WID=3.*AUWI(IW)**2/32./PI/MWIW*(1.+AMUP**2/MWIW**2-
 870      $     AMDLSS**2/MWIW**2)*SQRT(SSXLAM(MWIW**2,AMUP**2,AMDLSS**2))
 871           CALL SSSAVE(ISWIW,WID,-ISDNL,+IDUP,0,0,0)
 872         END IF
 873         IF(MWIW.GT.(AMCLSS+AMST)) THEN
 874           WID=3.*ADWI(IW)**2/32./PI/MWIW*(1.+AMST**2/MWIW**2-
 875      $     AMCLSS**2/MWIW**2)*SQRT(SSXLAM(MWIW**2,AMST**2,AMCLSS**2))
 876           CALL SSSAVE(ISWIW,WID,+ISCHL,-IDST,0,0,0)
 877         END IF
 878         IF(MWIW.GT.(AMSLSS+AMCH)) THEN
 879           WID=3.*AUWI(IW)**2/32./PI/MWIW*(1.+AMCH**2/MWIW**2-
 880      $     AMCLSS**2/MWIW**2)*SQRT(SSXLAM(MWIW**2,AMCH**2,AMCLSS**2))
 881           CALL SSSAVE(ISWIW,WID,-ISSTL,+IDCH,0,0,0)
 882         ENDIF
 883         IF(MWIW.GT.(AMT1SS+AMBT)) THEN
 884           AS=(-ADWI(IW)*COST-BWI(IW)*SINT)**2
 885           WID=3*((AS+BPWI(IW)**2*COST**2)*(MWIW**2+AMBT**2-AMT1SS**2)
 886      $        +4*SQRT(AS)*BPWI(IW)*COST*AMBT*MWIW)/32./PI/MWIW**3*
 887      $      SQRT(SSXLAM(MWIW**2,AMBT**2,AMT1SS**2))
 888           CALL SSSAVE(ISWIW,WID,ISTP1,-IDBT,0,0,0)
 889         ENDIF
 890         IF(MWIW.GT.(AMT2SS+AMBT)) THEN
 891           AS=(-ADWI(IW)*SINT+BWI(IW)*COST)**2
 892           WID=3*((AS+BPWI(IW)**2*SINT**2)*(MWIW**2+AMBT**2-AMT2SS**2)
 893      $        +4*SQRT(AS)*BPWI(IW)*SINT*AMBT*MWIW)/32./PI/MWIW**3*
 894      $      SQRT(SSXLAM(MWIW**2,AMBT**2,AMT2SS**2))
 895           CALL SSSAVE(ISWIW,WID,ISTP2,-IDBT,0,0,0)
 896         ENDIF
 897         IF(MWIW.GT.(AMB1SS+AMTP)) THEN
 898           AS=(-AUWI(IW)*COSB-BPWI(IW)*SINB)**2
 899           WID=3*((AS+BWI(IW)**2*COSB**2)*(MWIW**2+AMTP**2-AMB1SS**2)
 900      $        +4*SQRT(AS)*BWI(IW)*COSB*AMTP*MWIW)/32./PI/MWIW**3*
 901      $      SQRT(SSXLAM(MWIW**2,AMTP**2,AMB1SS**2))
 902           CALL SSSAVE(ISWIW,WID,-ISBT1,IDTP,0,0,0)
 903         ENDIF
 904         IF(MWIW.GT.(AMB2SS+AMTP)) THEN
 905           AS=(-AUWI(IW)*SINB+BPWI(IW)*COSB)**2
 906           WID=3*((AS+BWI(IW)**2*SINB**2)*(MWIW**2+AMTP**2-AMB2SS**2)
 907      $        +4*SQRT(AS)*BWI(IW)*SINB*AMTP*MWIW)/32./PI/MWIW**3*
 908      $      SQRT(SSXLAM(MWIW**2,AMTP**2,AMB2SS**2))
 909           CALL SSSAVE(ISWIW,WID,-ISBT2,IDTP,0,0,0)
 910         ENDIF
 911 C
 912 C          wi --> lepton + slepton
 913 C
 914         IF(MWIW.GT.(AMN1SS+AME)) THEN
 915           AS=(-ALWI(IW))**2
 916           WID=AS*(MWIW**2+AME**2-AMN1SS**2)/32./PI/MWIW**3*
 917      $        SQRT(SSXLAM(MWIW**2,AME**2,AMN1SS**2))
 918           CALL SSSAVE(ISWIW,WID,ISNEL,-IDE,0,0,0)
 919         END IF
 920         IF(MWIW.GT.(AMN2SS+AMMU)) THEN
 921           AS=(-ALWI(IW))**2
 922           WID=AS*(MWIW**2+AMMU**2-AMN2SS**2)/32./PI/MWIW**3*
 923      $        SQRT(SSXLAM(MWIW**2,AMMU**2,AMN2SS**2))
 924           CALL SSSAVE(ISWIW,WID,ISNML,-IDMU,0,0,0)
 925         END IF
 926         IF(MWIW.GT.(AMN3SS+AMTAU)) THEN
 927           AS=(-ALWI(IW))**2
 928           WID=((AS+BPLWI(IW)**2)*(MWIW**2+AMTAU**2-AMN3SS**2)+
 929      $        4*SQRT(AS)*BPLWI(IW)*AMTAU*MWIW)/32./PI/MWIW**3*
 930      $        SQRT(SSXLAM(MWIW**2,AMTAU**2,AMN3SS**2))
 931           CALL SSSAVE(ISWIW,WID,ISNTL,-IDTAU,0,0,0)
 932         ENDIF
 933         IF(MWIW.GT.AMELSS) THEN
 934           WID=ANWI(IW)**2*(MWIW**2-AMELSS**2)**2/32./PI/MWIW**3
 935           CALL SSSAVE(ISWIW,WID,-ISEL,IDNE,0,0,0)
 936         ENDIF
 937         IF(MWIW.GT.AMMLSS) THEN
 938           WID=ANWI(IW)**2*(MWIW**2-AMMLSS**2)**2/32./PI/MWIW**3
 939           CALL SSSAVE(ISWIW,WID,-ISMUL,IDNM,0,0,0)
 940         END IF
 941         IF(MWIW.GT.AML1SS) THEN
 942           AS=(-ANWI(IW)*COSL-BPLWI(IW)*SINL)**2
 943           WID=AS*(MWIW**2-AML1SS**2)**2/32./PI/MWIW**3
 944           CALL SSSAVE(ISWIW,WID,-ISTAU1,IDNT,0,0,0)
 945         END IF
 946         IF(MWIW.GT.AML2SS) THEN
 947           AS=(-ANWI(IW)*SINL+BPLWI(IW)*COSL)**2
 948           WID=AS*(MWIW**2-AML2SS**2)**2/32./PI/MWIW**3
 949           CALL SSSAVE(ISWIW,WID,-ISTAU2,IDNT,0,0,0)
 950         END IF
 951 300   CONTINUE
 952 C
 953 C          w2 --> w1 + z
 954 C          w2 --> w1 + f + fbar
 955 C
 956       IF (MW2.GT.(MW1+AMZ)) THEN
 957         EF=MW2**2+MW1**2-AMZ**2+((MW2**2-MW1**2)**2-
 958      $     AMZ**4)/AMZ/AMZ
 959         Y=(THX*SIN(GAMMAL)*COS(GAMMAL)-THY*SIN(GAMMAR)*COS(GAMMAR))/2.
 960         Z=(THX*SIN(GAMMAL)*COS(GAMMAL)+THY*SIN(GAMMAR)*COS(GAMMAR))/2.
 961         WID=(COTW+TANW)**2*SQRT(SSXLAM(MW2**2,MW1**2,AMZ**2))/
 962      $      128./32./MW2**3*(2*EF*(Y*Y+Z*Z)+
 963      $      12*MW2*MW1*(Y*Y-Z*Z)*SNW2*SNW1)
 964         CALL SSSAVE(ISW2,WID,ISW1,IDZ,0,0,0)
 965         W21ZL=0.
 966         W21ZN=0.
 967         W21ZU=0.
 968         W21ZD=0.
 969 C          ...w1 + f + fbar
 970       ELSE
 971         CONST=E**4*(COTW+TANW)**2/96./PI**3/MW2
 972         UPPER=(MW2**2+MW1**2)/2./MW2
 973         TMP(1)=MW2
 974         TMP(2)=MW1
 975         TMP(3)=AMZ
 976         TMP(4)=SNW1*SNW2
 977         TMP(5)=XWINO
 978         TMP(6)=YWINO
 979         XINTGL=SSXINT(MW1,SSWWF1,UPPER)
 980         W21ZL=(APL**2+BTL**2)*CONST*XINTGL
 981         W21ZN=(APN**2+BTN**2)*CONST*XINTGL
 982         W21ZU=3*(APU**2+BTU**2)*CONST*XINTGL
 983         W21ZD=3*(APD**2+BTD**2)*CONST*XINTGL
 984       END IF
 985 C     do w2 ->w1+q+qbar via sq'
 986       TMP(1)=MW2
 987       TMP(3)=MW1
 988       UPPER=(MW2**2-MW1**2)/2./MW2
 989       IF (MW2.LT.AMULSS) THEN
 990         TMP(2)=AMULSS
 991         XINTGL=SSXINT(0.,SSWZF2,UPPER)
 992         W21D=3*(ADWI(2)*ADWI(1))**2*XINTGL/2./MW2/(2*PI)**5
 993       ELSE
 994         W21D=0.
 995       END IF
 996       IF (MW2.LT.AMDLSS) THEN
 997         TMP(2)=AMDLSS
 998         XINTGL=SSXINT(0.,SSWZF2,UPPER)
 999         W21U=3*(AUWI(2)*AUWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1000       ELSE
1001         W21U=0.
1002       END IF
1003       IF (MW2.LT.AMCLSS) THEN
1004         TMP(2)=AMCLSS
1005         XINTGL=SSXINT(0.,SSWZF2,UPPER)
1006         W21S=3*(ADWI(2)*ADWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1007       ELSE
1008         W21S=0.
1009       END IF
1010       IF (MW2.LT.AMSLSS) THEN
1011         TMP(2)=AMSLSS
1012         XINTGL=SSXINT(0.,SSWZF2,UPPER)
1013         W21C=3*(AUWI(2)*AUWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1014       ELSE
1015         W21C=0.
1016       END IF
1017       IF (MW2.LT.AMN1SS) THEN
1018         TMP(2)=AMN1SS
1019         XINTGL=SSXINT(0.,SSWZF2,UPPER)
1020         W21E=(ALWI(2)*ALWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1021       ELSE
1022         W21E=0.
1023       END IF
1024       IF (MW2.LT.AMN2SS) THEN
1025         TMP(2)=AMN2SS
1026         XINTGL=SSXINT(0.,SSWZF2,UPPER)
1027         W21M=(ALWI(2)*ALWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1028       ELSE
1029         W21M=0.
1030       END IF
1031       IF (MW2.LT.AMN3SS) THEN
1032         TMP(2)=AMN3SS
1033         XINTGL=SSXINT(0.,SSWZF2,UPPER)
1034         W21L=(ALWI(2)*ALWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1035       ELSE
1036         W21L=0.
1037       END IF
1038       IF (MW2.LT.AMELSS) THEN
1039         TMP(2)=AMELSS
1040         XINTGL=SSXINT(0.,SSWZF2,UPPER)
1041         W21N1=(ANWI(2)*ANWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1042       ELSE
1043         W21N1=0.
1044       END IF
1045       IF (MW2.LT.AMMLSS) THEN
1046         TMP(2)=AMMLSS
1047         XINTGL=SSXINT(0.,SSWZF2,UPPER)
1048         W21N2=(ANWI(2)*ANWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1049       ELSE
1050         W21N2=0.
1051       END IF
1052 C     !!! W2->W1+NU_TAU+NU_TAUBAR NEEDS UPDATING FOR STAU MIXING
1053       IF (MW2.LT.AML1SS) THEN
1054         TMP(2)=AML1SS
1055         XINTGL=SSXINT(0.,SSWZF2,UPPER)
1056         W21N3=(ANWI(2)*ANWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1057       ELSE
1058         W21N3=0.
1059       END IF
1060 C-----WINO-2 ->WINO-1 +BBBAR NEEDS UPDATE -------------------
1061 C-----WINO-2 ->WINO-1 +TTBAR NEEDS UPDATE -------------------
1062 C-----THESE ALL LACK INTERFERENCE TERMS AS WELL
1063       W21D=W21D+W21ZD
1064       W21U=W21U+W21ZU
1065       W21S=W21S+W21ZD
1066       W21C=W21C+W21ZU
1067       W21N1=W21N1+W21ZN
1068       W21N2=W21N2+W21ZN
1069       W21N3=W21N3+W21ZN
1070       W21E=W21E+W21ZL
1071       W21M=W21M+W21ZL
1072       W21L=W21L+W21ZL
1073       IF(W21D.GT.0.) THEN
1074         CALL SSSAVE(ISW2,W21D,ISW1,IDDN,-IDDN,0,0)
1075       END IF
1076       IF(MW2.GT.(MW1+2*AMST+1.)) THEN
1077         CALL SSSAVE(ISW2,W21S,ISW1,IDST,-IDST,0,0)
1078       END IF
1079 C      IF (MW2.GT.(MW1+2*AMBT+2.)) THEN
1080 C        CALL SSSAVE(ISW2,W21D,ISW1,IDBT,-IDBT,0,0)
1081 C      END IF
1082       IF(W21U.GT.0.) THEN
1083         CALL SSSAVE(ISW2,W21U,ISW1,IDUP,-IDUP,0,0)
1084       END IF
1085       IF (MW2.GT.(MW1+2*AMCH+1.)) THEN
1086         CALL SSSAVE(ISW2,W21C,ISW1,IDCH,-IDCH,0,0)
1087       ENDIF
1088 C        IF (MW2.GT.(MW1+2*AMTP+2.)) THEN
1089 C          CALL SSSAVE(ISW2,W21U,ISW1,IDTP,-IDTP,0,0)
1090 C        END IF
1091       IF(W21N1.GT.0.) THEN
1092         CALL SSSAVE(ISW2,W21N1,ISW1,IDNE,-IDNE,0,0)
1093       ENDIF
1094       IF(W21N2.GT.0.) THEN
1095         CALL SSSAVE(ISW2,W21N2,ISW1,IDNM,-IDNM,0,0)
1096       ENDIF
1097       IF(W21N3.GT.0.) THEN
1098         CALL SSSAVE(ISW2,W21N3,ISW1,IDNT,-IDNT,0,0)
1099       ENDIF
1100       IF(W21E.GT.0.) THEN
1101         CALL SSSAVE(ISW2,W21E,ISW1,IDE,-IDE,0,0)
1102       ENDIF
1103       IF(MW2.GT.(MW1+2*AMMU+1.)) THEN
1104         CALL SSSAVE(ISW2,W21M,ISW1,IDMU,-IDMU,0,0)
1105       ENDIF
1106       IF(MW2.GT.(MW1+2*AMTAU+1.)) THEN
1107         CALL SSSAVE(ISW2,W21L,ISW1,IDTAU,-IDTAU,0,0)
1108       ENDIF
1109 C
1110 C          w2 --> w1 + higgs
1111 C
1112 C          w2 --> w1 + hl
1113       IF(MW2.GT.(MW1+AMHL)) THEN
1114         WID=G*G/16./PI/MW2**3*SQRT(SSXLAM(MW2**2,MW1**2,AMHL**2))*
1115      $  ((SL*SL+PL*PL)*(MW2*MW2+MW1*MW1-AMHL*AMHL)/2.+
1116      $  (SL*SL-PL*PL)*MW2*MW1)
1117         CALL SSSAVE(ISW2,WID,ISW1,ISHL,0,0,0)
1118       ENDIF
1119 C          w2 --> w1 + hh
1120       IF(MW2.GT.(MW1+AMHH)) THEN
1121         WID=G*G/16./PI/MW2**3*SQRT(SSXLAM(MW2**2,MW1**2,AMHH**2))*
1122      $  ((SH*SH+PH*PH)*(MW2*MW2+MW1*MW1-AMHH*AMHH)/2.+
1123      $  (SH*SH-PH*PH)*MW2*MW1)
1124         CALL SSSAVE(ISW2,WID,ISW1,ISHH,0,0,0)
1125       ENDIF
1126 C          w2 --> w1 + ha
1127       IF(MW2.GT.(MW1+AMHA)) THEN
1128         WID=G*G/16./PI/MW2**3*SQRT(SSXLAM(MW2**2,MW1**2,AMHA**2))*
1129      $  ((SP*SP+PP*PP)*(MW2*MW2+MW1*MW1-AMHA*AMHA)/2.+
1130      $  (SP*SP-PP*PP)*MW2*MW1)
1131         CALL SSSAVE(ISW2,WID,ISW1,ISHA,0,0,0)
1132       END IF
1133 C
1134 C          Normalize wi branching ratios
1135 C
1136       CALL SSNORM(ISW1)
1137       CALL SSNORM(ISW2)
1138 C
1139       RETURN
1140       END