ISAJET/isasusy/sslpbf.F

   1 #include "isajet/pilot.h"
   2         SUBROUTINE SSLPBF
   3 C--------------------------------------------------------
   4 C
   5 C        This program gives slepton branching fractions to gauginos
   6 C        according to Baer,Bartl,Karatas,Majerotto,Tata
   7 C        (Int. J. Mod. Phys. A4,4111 (1989); updated 10/21/96
   8 C
   9 C
  10 C-------------------------------------------------------
  11 #if defined(CERNLIB_IMPNONE)
  12       IMPLICIT NONE
  13 #endif
  14 #include "isajet/sslun.inc"
  15 #include "isajet/ssmode.inc"
  16 #include "isajet/sssm.inc"
  17 #include "isajet/sspar.inc"
  18 #include "isajet/sstype.inc"
  19 #include "isajet/sstmp.inc"
  20 C
  21       EXTERNAL SSLRT1,SSL1ST,SSN1ST,SSSNWS
  22       REAL SSXLAM,SSXINT,SSLRT1,SSN1ST,SSL1ST,SSSNWS
  23       REAL WID,SNZI,THIZ,XM,YM,THX,THY,BPLWI(2),CS2THW
  24       REAL VS,PI,SR2,G,GP,VP,V,TANB,ANWI(2),ALWI(2)
  25       REAL ANIZ,AEIZ,BEIZ,MZIZ,SINL,COSL,BETA,FL,AMPL
  26       REAL MTAMTA,MTAMB,MTAMZ,SUALFE
  27       REAL MW1,MW2,TN2THW,SNW1,SNW2,AS,BS,BH,COSA,SINA,A
  28       REAL XLO,SUM1,SUM2,AEJZ,TERM,WID1,WID2,SINB,COSB,COS2B,ANJZ,EMAX
  29       REAL MEME,MEMB,MEMZ,MMMM,MMUMB,MMUMZ,TANTHE,TANTHM,THETAE,THETAM
  30       REAL AAE,FE,AAM,FM,AME1,AML1,APE1,APL1,TM1,TM2,TM3,AMM1,APM1
  31       REAL BME1,BPE1
  32       INTEGER IZ,JZ
  33       INTEGER ISZIZ(4)
  34       COMPLEX ZI,ZONE,ZA,ZB,ZALIZ,ZBLIZ,ZPP,ZPM
  35       DATA ZONE/(1.,0.)/,ZI/(0.,1.)/
  36 C
  37 C          Partly duplicated from SSMASS
  38 C
  39       AMPL=2.4E18
  40       PI=4.*ATAN(1.)
  41       SR2=SQRT(2.)
  42       G=SQRT(4*PI*ALFAEM/SN2THW)
  43       GP=G*SQRT(SN2THW/(1.-SN2THW))
  44       CS2THW=1.-SN2THW
  45       TN2THW=SN2THW/CS2THW
  46       VS=2*AMW**2/G/G/(1.+RV2V1**2)
  47       V=SQRT(VS)
  48       VP=RV2V1*V
  49       TANB=1./RV2V1
  50       BETA=ATAN(TANB)
  51       SINA=SIN(ALFAH)
  52       COSA=COS(ALFAH)
  53       SINB=SIN(BETA)
  54       COSB=COS(BETA)
  55       COS2B=COS(2*BETA)
  56       XM=1./TAN(GAMMAL)
  57       YM=1./TAN(GAMMAR)
  58       THX=SIGN(1.,XM)
  59       THY=SIGN(1.,YM)
  60       MEME=AME*(1.-SUALFE(AME**2)/PI)
  61       MEMB=MEME*(SUALFE(AMBT**2)/SUALFE(AME**2))**(-27./76.)
  62       MEMZ=MEMB*(SUALFE(AMZ**2)/SUALFE(AMBT**2))**(-27./80.)
  63       FE=G*MEMZ/SR2/AMW/COS(BETA)
  64 C     SINCE A_e not defined in ISAJET, use A_tau as approximation
  65       AAE=AAL
  66       TANTHE=(AMERSS**2-MEMZ**2+AMZ**2*COS2B*(.5-SN2THW)-
  67      $AMELSS**2)/MEMZ/(TWOM1*SINB/COSB+AAE)
  68       THETAE=ATAN(TANTHE)
  69       MMMM=AMMU*(1.-SUALFE(AMMU**2)/PI)
  70       MMUMB=MMMM*(SUALFE(AMBT**2)/SUALFE(AMMU**2))**(-27./76.)
  71       MMUMZ=MMUMB*(SUALFE(AMZ**2)/SUALFE(AMBT**2))**(-27./80.)
  72       FM=G*MMUMZ/SR2/AMW/COS(BETA)
  73       AAM=AAL
  74       TANTHM=(AMMRSS**2-MMUMZ**2+AMZ**2*COS2B*(.5-SN2THW)-
  75      $AMMLSS**2)/MMUMZ/(TWOM1*SINB/COSB+AAM)
  76       THETAM=ATAN(TANTHM)
  77       MTAMTA=AMTAU*(1.-SUALFE(AMTAU**2)/PI)
  78       MTAMB=MTAMTA*(SUALFE(AMBT**2)/SUALFE(AMTAU**2))**(-27./76.)
  79       MTAMZ=MTAMB*(SUALFE(AMZ**2)/SUALFE(AMBT**2))**(-27./80.)
  80       FL=G*MTAMZ/SR2/AMW/COS(BETA)
  81       SINL=SIN(THETAL)
  82       COSL=COS(THETAL)
  83       SNW1=SIGN(1.,AMW1SS)
  84       SNW2=SIGN(1.,AMW2SS)
  85       BPLWI(1)=-FL*COS(GAMMAL)
  86       BPLWI(2)=FL*THX*SIN(GAMMAL)
  87       ANWI(1)=G*SIN(GAMMAL)
  88       ALWI(1)=SNW1*G*SIN(GAMMAR)
  89       ANWI(2)=G*THX*COS(GAMMAL)
  90       ALWI(2)=SNW2*G*THY*COS(GAMMAR)
  91 C          Reconstruct masses from SSMASS
  92       MW1=ABS(AMW1SS)
  93       MW2=ABS(AMW2SS)
  94 C
  95 C          Compute slepton branching fractions to zi
  96 C
  97       ISZIZ(1)=ISZ1
  98       ISZIZ(2)=ISZ2
  99       ISZIZ(3)=ISZ3
 100       ISZIZ(4)=ISZ4
 101       DO 100 IZ=1,4
 102         SNZI=SIGN(1.,AMZISS(IZ))
 103         IF (SNZI.EQ.1.) THEN
 104            THIZ=0.
 105         ELSE
 106            THIZ=1.
 107         END IF
 108         MZIZ=ABS(AMZISS(IZ))
 109         ANIZ=G/SR2*ZMIXSS(3,IZ)-GP/SR2*ZMIXSS(4,IZ)
 110         AEIZ=G/SR2*ZMIXSS(3,IZ)+GP/SR2*ZMIXSS(4,IZ)
 111         BEIZ=SR2*GP*ZMIXSS(4,IZ)
 112         ZALIZ=ZI**(THIZ-1.)*SNZI
 113      $  *(G/SR2*ZMIXSS(3,IZ)+GP/SR2*ZMIXSS(4,IZ))
 114         ZBLIZ=-1*ZI**(THIZ-1.)*SR2*GP*ZMIXSS(4,IZ)
 115         ZPP=ZI**THIZ
 116         ZPM=(-ZI)**THIZ
 117 C          sLEPTON --> LEPTON + zi
 118         IF (AMELSS.GT.(MZIZ+AME)) THEN
 119           WID=AEIZ**2*(AMELSS**2-MZIZ**2-AME**2)/AMELSS**3
 120      $        /16./PI*SQRT(SSXLAM(AMELSS**2,MZIZ**2,AME**2))
 121           CALL SSSAVE(ISEL,WID,ISZIZ(IZ),IDE,0,0,0)
 122         ENDIF
 123         IF (AMMLSS.GT.(MZIZ+AMMU)) THEN
 124           WID=AEIZ**2*(AMMLSS**2-MZIZ**2-AMMU**2)/AMMLSS**3
 125      $        /16./PI*SQRT(SSXLAM(AMMLSS**2,MZIZ**2,AMMU**2))
 126           CALL SSSAVE(ISMUL,WID,ISZIZ(IZ),IDMU,0,0,0)
 127         ENDIF
 128         IF (AMN1SS.GT.MZIZ) THEN
 129           WID=ANIZ**2*AMN1SS*(1.-MZIZ**2/AMN1SS**2)**2/16./PI
 130           CALL SSSAVE(ISNEL,WID,ISZIZ(IZ),IDNE,0,0,0)
 131         ENDIF
 132         IF (AMN2SS.GT.MZIZ) THEN
 133           WID=ANIZ**2*AMN2SS*(1.-MZIZ**2/AMN2SS**2)**2/16./PI
 134           CALL SSSAVE(ISNML,WID,ISZIZ(IZ),IDNM,0,0,0)
 135         ENDIF
 136         IF (AMN3SS.GT.MZIZ) THEN
 137           WID=ANIZ**2*AMN3SS*(1.-MZIZ**2/AMN3SS**2)**2/16./PI
 138           CALL SSSAVE(ISNTL,WID,ISZIZ(IZ),IDNT,0,0,0)
 139         ENDIF
 140         IF (AMERSS.GT.(MZIZ+AME)) THEN
 141           WID=BEIZ**2*(AMERSS**2-MZIZ**2-AME**2)/AMERSS**3
 142      $        /16./PI*SQRT(SSXLAM(AMERSS**2,MZIZ**2,AME**2))
 143           CALL SSSAVE(ISER,WID,ISZIZ(IZ),IDE,0,0,0)
 144         ENDIF
 145         IF (AMMRSS.GT.(MZIZ+AMMU)) THEN
 146           WID=BEIZ**2*(AMMRSS**2-MZIZ**2-AMMU**2)/AMMRSS**3
 147      $        /16./PI*SQRT(SSXLAM(AMMRSS**2,MZIZ**2,AMMU**2))
 148           CALL SSSAVE(ISMUR,WID,ISZIZ(IZ),IDMU,0,0,0)
 149         ENDIF
 150 C          sTAU_1 --> TAU + zi
 151         IF (AML1SS.GT.(MZIZ+AMTAU)) THEN
 152           ZA=((ZI*ZALIZ-ZPP*FL*ZMIXSS(2,IZ))*COSL
 153      $     -(ZI*ZBLIZ-ZPM*FL*ZMIXSS(2,IZ))*SINL)/2.
 154           ZB=((-ZI*ZALIZ-ZPP*FL*ZMIXSS(2,IZ))*COSL
 155      $     -(ZI*ZBLIZ+ZPM*FL*ZMIXSS(2,IZ))*SINL)/2.
 156           AS=ZA*CONJG(ZA)
 157           BS=ZB*CONJG(ZB)
 158           WID=(AS*(AML1SS**2-(AMTAU+MZIZ)**2)+BS*(AML1SS**2-
 159      $     (MZIZ-AMTAU)**2))/8./PI/AML1SS**3*
 160      $     SQRT(SSXLAM(AML1SS**2,MZIZ**2,AMTAU**2))
 161           CALL SSSAVE(ISTAU1,WID,ISZIZ(IZ),IDTAU,0,0,0)
 162         END IF
 163         IF (AML2SS.GT.(MZIZ+AMTAU)) THEN
 164           ZA=((ZI*ZALIZ-ZPP*FL*ZMIXSS(2,IZ))*SINL
 165      $     +(ZI*ZBLIZ-ZPM*FL*ZMIXSS(2,IZ))*COSL)/2.
 166           ZB=((-ZI*ZALIZ-ZPP*FL*ZMIXSS(2,IZ))*SINL
 167      $     +(ZI*ZBLIZ+ZPM*FL*ZMIXSS(2,IZ))*COSL)/2.
 168           AS=ZA*CONJG(ZA)
 169           BS=ZB*CONJG(ZB)
 170           WID=(AS*(AML2SS**2-(AMTAU+MZIZ)**2)+BS*(AML2SS**2-
 171      $     (MZIZ-AMTAU)**2))/8./PI/AML2SS**3*
 172      $     SQRT(SSXLAM(AML2SS**2,MZIZ**2,AMTAU**2))
 173           CALL SSSAVE(ISTAU2,WID,ISZIZ(IZ),IDTAU,0,0,0)
 174         END IF
 175 100   CONTINUE
 176 C
 177 C           Compute branching fractions to wi --- theta-C = 0
 178 C
 179       IF (AMELSS.GT.MW1) THEN
 180         WID=ANWI(1)**2*AMELSS*(1.-MW1**2/AMELSS**2)**2/16./PI
 181         CALL SSSAVE(ISEL,WID,-ISW1,IDNE,0,0,0)
 182       END IF
 183       IF (AMMLSS.GT.MW1) THEN
 184         WID=ANWI(1)**2*AMMLSS*(1.-MW1**2/AMMLSS**2)**2/16./PI
 185         CALL SSSAVE(ISMUL,WID,-ISW1,IDNM,0,0,0)
 186       END IF
 187       IF (AML1SS.GT.MW1) THEN
 188         AS=(-ANWI(1)*COSL-BPLWI(1)*SINL)**2
 189         WID=AS*AML1SS*(1.-MW1**2/AML1SS**2)**2/16./PI
 190         CALL SSSAVE(ISTAU1,WID,-ISW1,IDNT,0,0,0)
 191       END IF
 192       IF (AML2SS.GT.MW1) THEN
 193         AS=(-ANWI(1)*SINL+BPLWI(1)*COSL)**2
 194         WID=AS*AML2SS*(1.-MW1**2/AML2SS**2)**2/16./PI
 195         CALL SSSAVE(ISTAU2,WID,-ISW1,IDNT,0,0,0)
 196       END IF
 197 C
 198       IF (AMN1SS.GT.(MW1+AME)) THEN
 199         WID=ALWI(1)**2*(AMN1SS**2-MW1**2-AME**2)*
 200      $      SQRT(SSXLAM(AMN1SS**2,MW1**2,AME**2))/16./PI/AMN1SS**3
 201         CALL SSSAVE(ISNEL,WID,ISW1,IDE,0,0,0)
 202       END IF
 203       IF (AMN2SS.GT.(MW1+AMMU)) THEN
 204         WID=ALWI(1)**2*(AMN2SS**2-MW1**2-AMMU**2)*
 205      $      SQRT(SSXLAM(AMN2SS**2,MW1**2,AMMU**2))/16./PI/AMN2SS**3
 206         CALL SSSAVE(ISNML,WID,ISW1,IDMU,0,0,0)
 207       END IF
 208 C
 209       IF (AMN3SS.GT.(MW1+AMTAU)) THEN
 210         WID=((ALWI(1)**2+BPLWI(1)**2)*(AMN3SS**2-MW1**2-AMTAU**2)+
 211      $     4*MW1*AMTAU*BPLWI(1)*ALWI(1))*
 212      $     SQRT(SSXLAM(AMN3SS**2,MW1**2,AMTAU**2))/16./PI/AMN3SS**3
 213         CALL SSSAVE(ISNTL,WID,ISW1,IDTAU,0,0,0)
 214       END IF
 215 C
 216 C
 217       IF (AMELSS.GT.MW2) THEN
 218         WID=ANWI(2)**2*AMELSS*(1.-MW2**2/AMELSS**2)**2/16./PI
 219         CALL SSSAVE(ISEL,WID,-ISW2,IDNE,0,0,0)
 220       END IF
 221       IF (AMMLSS.GT.MW2) THEN
 222         WID=ANWI(2)**2*AMMLSS*(1.-MW2**2/AMMLSS**2)**2/16./PI
 223         CALL SSSAVE(ISMUL,WID,-ISW2,IDNM,0,0,0)
 224       END IF
 225       IF (AML1SS.GT.MW2) THEN
 226         AS=(-ANWI(2)*COSL-BPLWI(2)*SINL)**2
 227         WID=AS*AML1SS*(1.-MW2**2/AML1SS**2)**2/16./PI
 228         CALL SSSAVE(ISTAU1,WID,-ISW2,IDNT,0,0,0)
 229       END IF
 230       IF (AML2SS.GT.MW2) THEN
 231         AS=(-ANWI(2)*SINL+BPLWI(2)*COSL)**2
 232         WID=AS*AML2SS*(1.-MW2**2/AML2SS**2)**2/16./PI
 233         CALL SSSAVE(ISTAU2,WID,-ISW2,IDNT,0,0,0)
 234       END IF
 235 C
 236       IF (AMN1SS.GT.(MW2+AME)) THEN
 237         WID=ALWI(2)**2*(AMN1SS**2-MW2**2-AME**2)*
 238      $      SQRT(SSXLAM(AMN1SS**2,MW2**2,AME**2))/16./PI/AMN1SS**3
 239         CALL SSSAVE(ISNEL,WID,ISW2,IDE,0,0,0)
 240       END IF
 241       IF (AMN2SS.GT.(MW2+AMMU)) THEN
 242         WID=ALWI(2)**2*(AMN2SS**2-MW2**2-AMMU**2)*
 243      $      SQRT(SSXLAM(AMN2SS**2,MW2**2,AMMU**2))/16./PI/AMN2SS**3
 244         CALL SSSAVE(ISNML,WID,ISW2,IDMU,0,0,0)
 245       END IF
 246 C
 247       IF (AMN3SS.GT.(MW2+AMTAU)) THEN
 248         WID=((ALWI(2)**2+BPLWI(2)**2)*(AMN3SS**2-MW2**2-AMTAU**2)+
 249      $     4*MW2*AMTAU*BPLWI(2)*ALWI(2))*
 250      $     SQRT(SSXLAM(AMN3SS**2,MW2**2,AMTAU**2))/16./PI/AMN3SS**3
 251         CALL SSSAVE(ISNTL,WID,ISW2,IDTAU,0,0,0)
 252       END IF
 253 C
 254 C      DECAYS TO VECTOR BOSONS
 255 C
 256       IF (AMELSS.GT.(AMN1SS+AMW)) THEN
 257         WID=G*G*(SSXLAM(AMELSS**2,AMN1SS**2,AMW**2))**1.5/
 258      $      32./PI/AMELSS**3/AMW**2
 259         CALL SSSAVE(ISEL,WID,-IDW,ISNEL,0,0,0)
 260       END IF
 261 C
 262       IF (AMMLSS.GT.(AMN2SS+AMW)) THEN
 263         WID=G*G*(SSXLAM(AMMLSS**2,AMN2SS**2,AMW**2))**1.5/
 264      $      32./PI/AMMLSS**3/AMW**2
 265         CALL SSSAVE(ISMUL,WID,-IDW,ISNML,0,0,0)
 266       END IF
 267 C
 268       IF (AML1SS.GT.(AMN3SS+AMW)) THEN
 269         WID=G*G*(SSXLAM(AML1SS**2,AMN3SS**2,AMW**2))**1.5/
 270      $      32./PI/AML1SS**3/AMW**2*COSL**2
 271         CALL SSSAVE(ISTAU1,WID,-IDW,ISNTL,0,0,0)
 272       END IF
 273 C
 274       IF (AML2SS.GT.(AMN3SS+AMW)) THEN
 275         WID=G*G*(SSXLAM(AML2SS**2,AMN3SS**2,AMW**2))**1.5/
 276      $      32./PI/AML2SS**3/AMW**2*SINL**2
 277         CALL SSSAVE(ISTAU2,WID,-IDW,ISNTL,0,0,0)
 278       END IF
 279 C
 280       IF (AMN1SS.GT.(AMELSS+AMW)) THEN
 281         WID=G*G*(SSXLAM(AMN1SS**2,AMELSS**2,AMW**2))**1.5/
 282      $      32./PI/AMN1SS**3/AMW**2
 283         CALL SSSAVE(ISNEL,WID,IDW,ISEL,0,0,0)
 284       END IF
 285 C
 286       IF (AMN2SS.GT.(AMMLSS+AMW)) THEN
 287         WID=G*G*(SSXLAM(AMN2SS**2,AMMLSS**2,AMW**2))**1.5/
 288      $      32./PI/AMN2SS**3/AMW**2
 289         CALL SSSAVE(ISNML,WID,IDW,ISMUL,0,0,0)
 290       END IF
 291 C
 292       IF (AMN3SS.GT.(AML1SS+AMW)) THEN
 293         WID=G*G*(SSXLAM(AMN3SS**2,AML1SS**2,AMW**2))**1.5/
 294      $      32./PI/AMN3SS**3/AMW**2*COSL**2
 295         CALL SSSAVE(ISNTL,WID,IDW,ISTAU1,0,0,0)
 296       END IF
 297 C
 298       IF (AMN3SS.GT.(AML2SS+AMW)) THEN
 299         WID=G*G*(SSXLAM(AMN3SS**2,AML2SS**2,AMW**2))**1.5/
 300      $      32./PI/AMN3SS**3/AMW**2*SINL**2
 301         CALL SSSAVE(ISNTL,WID,IDW,ISTAU2,0,0,0)
 302       END IF
 303 C
 304       IF (AML2SS.GT.(AML1SS+AMZ)) THEN
 305         WID=G*G*(SSXLAM(AML2SS**2,AML1SS**2,AMZ**2))**1.5/
 306      $      64./PI/AML2SS**3/CS2THW/AMZ**2*SINL**2*COSL**2
 307         CALL SSSAVE(ISTAU2,WID,IDZ,ISTAU1,0,0,0)
 308       END IF
 309 C
 310 C     3-body decay of l_R -> l+tau+stau_1
 311 C
 312       XLO=(AML1SS+AMTAU)**2
 313       SUM1=0.
 314       SUM2=0.
 315       IF (AMERSS.GT.(AML1SS+AMTAU+AME)) THEN
 316         DO IZ=1,4
 317           DO JZ=IZ,4
 318             TMP(1)=AMERSS
 319             TMP(2)=-SR2*GP*ZMIXSS(4,IZ)
 320             TMP(3)=-SR2*GP*ZMIXSS(4,JZ)
 321             AEIZ=-(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))/SR2
 322             AEJZ=-(G*ZMIXSS(3,JZ)+GP*ZMIXSS(4,JZ))/SR2
 323             TMP(4)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
 324             TMP(5)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
 325             TMP(6)=TMP(2)*SINL+FL*ZMIXSS(2,IZ)*COSL
 326             TMP(7)=TMP(3)*SINL+FL*ZMIXSS(2,JZ)*COSL
 327             TMP(8)=AMZISS(IZ)
 328             TMP(9)=AMZISS(JZ)
 329             IF (AMERSS.LT.ABS(AMZISS(IZ)).AND.
 330      ,          AMERSS.LT.ABS(AMZISS(JZ))) THEN
 331               TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMERSS**2)/AMERSS
 332             ELSE
 333               TERM=0.
 334             END IF
 335             IF (IZ.EQ.JZ) TERM=TERM/2.
 336             SUM1=SUM1+TERM
 337             TMP(4)=TMP(6)
 338             TMP(5)=TMP(7)
 339             TMP(6)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
 340             TMP(7)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
 341             IF (AMERSS.LT.ABS(AMZISS(IZ)).AND.
 342      ,          AMERSS.LT.ABS(AMZISS(JZ))) THEN
 343               TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMERSS**2)/AMERSS
 344             ELSE
 345               TERM=0.
 346             END IF
 347             IF (IZ.EQ.JZ) TERM=TERM/2.
 348             SUM2=SUM2+TERM
 349           END DO
 350         END DO
 351         WID1=SUM1/2./AMERSS/(2*PI)**5
 352         WID2=SUM2/2./AMERSS/(2*PI)**5
 353         CALL SSSAVE(ISER,WID1,ISTAU1,IDE,-IDTAU,0,0)
 354         CALL SSSAVE(ISER,WID2,-ISTAU1,IDE,IDTAU,0,0)
 355       END IF
 356       SUM1=0.
 357       SUM2=0.
 358       IF (AMMRSS.GT.(AML1SS+AMTAU+AMMU)) THEN
 359         DO IZ=1,4
 360           DO JZ=IZ,4
 361             TMP(1)=AMMRSS
 362             TMP(2)=-SR2*GP*ZMIXSS(4,IZ)
 363             TMP(3)=-SR2*GP*ZMIXSS(4,JZ)
 364             AEIZ=-(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))/SR2
 365             AEJZ=-(G*ZMIXSS(3,JZ)+GP*ZMIXSS(4,JZ))/SR2
 366             TMP(4)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
 367             TMP(5)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
 368             TMP(6)=TMP(2)*SINL+FL*ZMIXSS(2,IZ)*COSL
 369             TMP(7)=TMP(3)*SINL+FL*ZMIXSS(2,JZ)*COSL
 370             TMP(8)=AMZISS(IZ)
 371             TMP(9)=AMZISS(JZ)
 372             IF (AMMRSS.LT.ABS(AMZISS(IZ)).AND.
 373      ,          AMMRSS.LT.ABS(AMZISS(JZ))) THEN
 374               TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMMRSS**2)/AMMRSS
 375             ELSE
 376               TERM=0.
 377             END IF
 378             IF (IZ.EQ.JZ) TERM=TERM/2.
 379             SUM1=SUM1+TERM
 380             TMP(4)=TMP(6)
 381             TMP(5)=TMP(7)
 382             TMP(6)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
 383             TMP(7)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
 384             IF (AMMRSS.LT.ABS(AMZISS(IZ)).AND.
 385      ,          AMMRSS.LT.ABS(AMZISS(JZ))) THEN
 386               TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMMRSS**2)/AMMRSS
 387             ELSE
 388               TERM=0.
 389             END IF
 390             IF (IZ.EQ.JZ) TERM=TERM/2.
 391             SUM2=SUM2+TERM
 392           END DO
 393         END DO
 394         WID1=SUM1/2./AMMRSS/(2*PI)**5
 395         WID2=SUM2/2./AMMRSS/(2*PI)**5
 396         CALL SSSAVE(ISMUR,WID1,ISTAU1,IDMU,-IDTAU,0,0)
 397         CALL SSSAVE(ISMUR,WID2,-ISTAU1,IDMU,IDTAU,0,0)
 398       END IF
 399 C
 400 C     3-body decay of l_L -> l+tau+stau_1
 401 C
 402       SUM1=0.
 403       SUM2=0.
 404       IF (AMELSS.GT.(AML1SS+AMTAU+AME)) THEN
 405         DO IZ=1,4
 406           DO JZ=IZ,4
 407             TMP(1)=AMELSS
 408             AEIZ=-(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))/SR2
 409             AEJZ=-(G*ZMIXSS(3,JZ)+GP*ZMIXSS(4,JZ))/SR2
 410             TMP(2)=AEIZ
 411             TMP(3)=AEJZ
 412             TMP(6)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
 413             TMP(7)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
 414             TMP(4)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
 415             TMP(5)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
 416             TMP(8)=AMZISS(IZ)
 417             TMP(9)=AMZISS(JZ)
 418             IF (AMELSS.LT.ABS(AMZISS(IZ)).AND.
 419      ,          AMELSS.LT.ABS(AMZISS(JZ))) THEN
 420               TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMELSS**2)/AMELSS
 421             ELSE
 422               TERM=0.
 423             END IF
 424             IF (IZ.EQ.JZ) TERM=TERM/2.
 425             SUM1=SUM1+TERM
 426             TMP(4)=TMP(6)
 427             TMP(5)=TMP(7)
 428             TMP(6)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
 429             TMP(7)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
 430             IF (AMELSS.LT.ABS(AMZISS(IZ)).AND.
 431      ,          AMELSS.LT.ABS(AMZISS(JZ))) THEN
 432               TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMELSS**2)/AMELSS
 433             ELSE
 434               TERM=0.
 435             END IF
 436             IF (IZ.EQ.JZ) TERM=TERM/2.
 437             SUM2=SUM2+TERM
 438           END DO
 439         END DO
 440         WID1=SUM1/2./AMELSS/(2*PI)**5
 441         WID2=SUM2/2./AMELSS/(2*PI)**5
 442         CALL SSSAVE(ISEL,WID1,ISTAU1,IDE,-IDTAU,0,0)
 443         CALL SSSAVE(ISEL,WID2,-ISTAU1,IDE,IDTAU,0,0)
 444       END IF
 445       SUM1=0.
 446       SUM2=0.
 447       IF (AMMLSS.GT.(AML1SS+AMTAU+AMMU)) THEN
 448         DO IZ=1,4
 449           DO JZ=IZ,4
 450             TMP(1)=AMMLSS
 451             AEIZ=-(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))/SR2
 452             AEJZ=-(G*ZMIXSS(3,JZ)+GP*ZMIXSS(4,JZ))/SR2
 453             TMP(2)=AEIZ
 454             TMP(3)=AEJZ
 455             TMP(6)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
 456             TMP(7)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
 457             TMP(4)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
 458             TMP(5)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
 459             TMP(8)=AMZISS(IZ)
 460             TMP(9)=AMZISS(JZ)
 461             IF (AMMLSS.LT.ABS(AMZISS(IZ)).AND.
 462      ,          AMMLSS.LT.ABS(AMZISS(JZ))) THEN
 463               TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMMLSS**2)/AMMLSS
 464             ELSE
 465               TERM=0.
 466             END IF
 467             IF (IZ.EQ.JZ) TERM=TERM/2.
 468             SUM1=SUM1+TERM
 469             TMP(4)=TMP(6)
 470             TMP(5)=TMP(7)
 471             TMP(6)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
 472             TMP(7)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
 473             IF (AMMLSS.LT.ABS(AMZISS(IZ)).AND.
 474      ,          AMMLSS.LT.ABS(AMZISS(JZ))) THEN
 475               TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMMLSS**2)/AMMLSS
 476             ELSE
 477               TERM=0.
 478             END IF
 479             IF (IZ.EQ.JZ) TERM=TERM/2.
 480             SUM2=SUM2+TERM
 481           END DO
 482         END DO
 483         WID1=SUM1/2./AMMLSS/(2*PI)**5
 484         WID2=SUM2/2./AMMLSS/(2*PI)**5
 485         CALL SSSAVE(ISMUL,WID1,ISTAU1,IDMU,-IDTAU,0,0)
 486         CALL SSSAVE(ISMUL,WID2,-ISTAU1,IDMU,IDTAU,0,0)
 487       END IF
 488 C
 489 C     3-body decay of nu_eL -> nu_e+tau+stau_1
 490 C
 491       SUM1=0.
 492       SUM2=0.
 493       IF (AMN1SS.GT.(AML1SS+AMTAU)) THEN
 494         DO IZ=1,4
 495           DO JZ=IZ,4
 496             TMP(1)=AMN1SS
 497             ANIZ=(G*ZMIXSS(3,IZ)-GP*ZMIXSS(4,IZ))/SR2
 498             ANJZ=(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/SR2
 499             AEIZ=-(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))/SR2
 500             AEJZ=-(G*ZMIXSS(3,JZ)+GP*ZMIXSS(4,JZ))/SR2
 501             TMP(2)=ANIZ
 502             TMP(3)=ANJZ
 503             TMP(6)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
 504             TMP(7)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
 505             TMP(4)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
 506             TMP(5)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
 507             TMP(8)=AMZISS(IZ)
 508             TMP(9)=AMZISS(JZ)
 509             IF (AMN1SS.LT.ABS(AMZISS(IZ)).AND.
 510      ,          AMN1SS.LT.ABS(AMZISS(JZ))) THEN
 511               TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMN1SS**2)/AMN1SS
 512             ELSE
 513               TERM=0.
 514             END IF
 515             IF (IZ.EQ.JZ) TERM=TERM/2.
 516             SUM1=SUM1+TERM
 517             TMP(4)=TMP(6)
 518             TMP(5)=TMP(7)
 519             TMP(6)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
 520             TMP(7)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
 521             IF (AMN1SS.LT.ABS(AMZISS(IZ)).AND.
 522      ,          AMN1SS.LT.ABS(AMZISS(JZ))) THEN
 523               TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMN1SS**2)/AMN1SS
 524             ELSE
 525               TERM=0.
 526             END IF
 527             IF (IZ.EQ.JZ) TERM=TERM/2.
 528             SUM2=SUM2+TERM
 529           END DO
 530         END DO
 531         WID1=SUM1/2./AMN1SS/(2*PI)**5
 532         WID2=SUM2/2./AMN1SS/(2*PI)**5
 533         CALL SSSAVE(ISNEL,WID1,ISTAU1,IDNE,-IDTAU,0,0)
 534         CALL SSSAVE(ISNEL,WID2,-ISTAU1,IDNE,IDTAU,0,0)
 535       END IF
 536 C
 537 C     3-body decay of nu_muL -> nu_mu+tau+stau_1
 538 C
 539       SUM1=0.
 540       SUM2=0.
 541       IF (AMN2SS.GT.(AML1SS+AMTAU)) THEN
 542         DO IZ=1,4
 543           DO JZ=IZ,4
 544             TMP(1)=AMN2SS
 545             ANIZ=(G*ZMIXSS(3,IZ)-GP*ZMIXSS(4,IZ))/SR2
 546             ANJZ=(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/SR2
 547             AEIZ=-(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))/SR2
 548             AEJZ=-(G*ZMIXSS(3,JZ)+GP*ZMIXSS(4,JZ))/SR2
 549             TMP(2)=ANIZ
 550             TMP(3)=ANJZ
 551             TMP(6)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
 552             TMP(7)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
 553             TMP(4)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
 554             TMP(5)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
 555             TMP(8)=AMZISS(IZ)
 556             TMP(9)=AMZISS(JZ)
 557             IF (AMN2SS.LT.ABS(AMZISS(IZ)).AND.
 558      ,          AMN2SS.LT.ABS(AMZISS(JZ))) THEN
 559               TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMN2SS**2)/AMN2SS
 560             ELSE
 561               TERM=0.
 562             END IF
 563             IF (IZ.EQ.JZ) TERM=TERM/2.
 564             SUM1=SUM1+TERM
 565             TMP(4)=TMP(6)
 566             TMP(5)=TMP(7)
 567             TMP(6)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
 568             TMP(7)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
 569             IF (AMN2SS.LT.ABS(AMZISS(IZ)).AND.
 570      ,          AMN2SS.LT.ABS(AMZISS(JZ))) THEN
 571               TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMN2SS**2)/AMN2SS
 572             ELSE
 573               TERM=0.
 574             END IF
 575             IF (IZ.EQ.JZ) TERM=TERM/2.
 576             SUM2=SUM2+TERM
 577           END DO
 578         END DO
 579         WID1=SUM1/2./AMN2SS/(2*PI)**5
 580         WID2=SUM2/2./AMN2SS/(2*PI)**5
 581         CALL SSSAVE(ISNML,WID1,ISTAU1,IDNM,-IDTAU,0,0)
 582         CALL SSSAVE(ISNML,WID2,-ISTAU1,IDNM,IDTAU,0,0)
 583       END IF
 584 C
 585 C     3-body decay of nu_tauL -> nu_tau+tau+stau_1
 586 C     Here, we include decay via Z_i, W_i and W*, but neglect
 587 C     some interference terms (Z_i-W* and W_i-W* and Z_i-W_i)
 588 C     until a future date. HB 8/24/98
 589 C
 590       AME1=-G*SIN(GAMMAL)*COS(THETAE)+FE*COS(GAMMAL)*SIN(THETAE)
 591       AML1=-G*SIN(GAMMAL)*COS(THETAL)+FL*COS(GAMMAL)*SIN(THETAL)
 592       APE1=-G*THX*COS(GAMMAL)*COS(THETAE)-FE*THX*SIN(GAMMAL)*
 593      $       SIN(THETAE)
 594       APL1=-G*THX*COS(GAMMAL)*COS(THETAL)-FL*THX*SIN(GAMMAL)*
 595      $       SIN(THETAL)
 596       SUM1=0.
 597       SUM2=0.
 598       IF (AMN3SS.GT.(AML1SS+AMTAU)) THEN
 599         DO IZ=1,4
 600           DO JZ=IZ,4
 601             TMP(1)=AMN3SS
 602             ANIZ=(G*ZMIXSS(3,IZ)-GP*ZMIXSS(4,IZ))/SR2
 603             ANJZ=(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/SR2
 604             AEIZ=-(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))/SR2
 605             AEJZ=-(G*ZMIXSS(3,JZ)+GP*ZMIXSS(4,JZ))/SR2
 606             TMP(2)=ANIZ
 607             TMP(3)=ANJZ
 608             TMP(6)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
 609             TMP(7)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
 610             TMP(4)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
 611             TMP(5)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
 612             TMP(8)=AMZISS(IZ)
 613             TMP(9)=AMZISS(JZ)
 614             IF (AMN3SS.LT.ABS(AMZISS(IZ)).AND.
 615      ,          AMN3SS.LT.ABS(AMZISS(JZ))) THEN
 616               TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMN3SS**2)/AMN3SS
 617             ELSE
 618               TERM=0.
 619             END IF
 620             IF (IZ.EQ.JZ) TERM=TERM/2.
 621             SUM1=SUM1+TERM
 622             TMP(4)=TMP(6)
 623             TMP(5)=TMP(7)
 624             TMP(6)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
 625             TMP(7)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
 626             IF (AMN3SS.LT.ABS(AMZISS(IZ)).AND.
 627      ,          AMN3SS.LT.ABS(AMZISS(JZ))) THEN
 628               TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMN3SS**2)/AMN3SS
 629             ELSE
 630               TERM=0.
 631             END IF
 632             IF (IZ.EQ.JZ) TERM=TERM/2.
 633             SUM2=SUM2+TERM
 634           END DO
 635         END DO
 636         WID1=SUM1/2./AMN3SS/(2*PI)**5
 637         WID2=SUM2/2./AMN3SS/(2*PI)**5
 638 C-----COMPUTE SNU_TAU --> TAU +STAU_1BAR +NU_TAU via W_i DECAYS --------------
 639         BME1=-FL*COS(GAMMAL)
 640         BPE1=FL*THX*SIN(GAMMAL)
 641         TMP(1)=MW1
 642         TMP(2)=MW1
 643         TMP(3)=AML1SS
 644         TMP(4)=AMN3SS
 645         IF (AMN3SS.LT.MW1) THEN
 646         TM1=AML1**2*(AME1**2*MW1**2*SSXINT(AML1SS**2,SSN1ST,AMN3SS**2)
 647      $+BME1**2*SSXINT(AML1SS**2,SSL1ST,AMN3SS**2))
 648         TMP(2)=MW2
 649         TM2=2*AML1*APL1*(SNW1*SNW2*AME1*APE1*
 650      $ SSXINT(AML1SS**2,SSN1ST,AMN3SS**2)+BME1*BPE1*
 651      $ SSXINT(AML1SS**2,SSL1ST,AMN3SS**2))
 652         ELSE
 653         TM1=0.
 654         TM2=0.
 655         END IF
 656         TMP(1)=MW2
 657         TMP(2)=MW2
 658         IF (AMN3SS.LT.MW2) THEN
 659         TM3=APL1**2*(APE1**2*MW2**2*SSXINT(AML1SS**2,SSN1ST,AMN3SS**2)
 660      $+BPE1**2*SSXINT(AML1SS**2,SSL1ST,AMN3SS**2))
 661         ELSE
 662         TM3=0.
 663         END IF
 664         WID=PI**2*(TM1+TM2+TM3)/8./2./AMN3SS/(2.*PI)**5
 665         CALL SSSAVE(ISNTL,WID1+WID,ISTAU1,IDNT,-IDTAU,0,0)
 666 C-----COMPUTE SNU_TAU --> TAUBAR +STAU_1 +NU_TAU via W* DECAYS --
 667         EMAX=(AML1SS**2+AMN3SS**2)/2./AMN3SS
 668         TMP(1)=AMW
 669         TMP(2)=AMW
 670         TMP(3)=AML1SS
 671         TMP(4)=AMN3SS
 672         IF (AMN3SS.LT.(AMW+AML1SS)) THEN
 673         TM1=G**4*COS(THETAL)**2*SSXINT(AML1SS,SSSNWS,EMAX)
 674         ELSE
 675         TM1=0.
 676         END IF
 677         WID=2*PI**2*TM1/3./2./AMN3SS/(2.*PI)**5
 678         CALL SSSAVE(ISNTL,WID2+WID,-ISTAU1,IDNT,IDTAU,0,0)
 679       END IF
 680 C-----Now impose all sneutrino 3 body decays via W* ------------
 681         TMP(1)=AMW
 682         TMP(2)=AMW
 683         TMP(3)=AML1SS
 684         TMP(4)=AMN3SS
 685         EMAX=(AML1SS**2+AMN3SS**2)/2./AMN3SS
 686       IF (AMN3SS.GT.(AML1SS+AME).AND.AMN3SS.LT.(AMW+AML1SS)) THEN
 687         TM1=G**4*COS(THETAL)**2*SSXINT(AML1SS,SSSNWS,EMAX)
 688         WID=2*PI**2*TM1/3./2./AMN3SS/(2.*PI)**5
 689         CALL SSSAVE(ISNTL,WID,ISTAU1,IDNE,-IDE,0,0)
 690       END IF
 691       IF (AMN3SS.GT.(AML1SS+AMMU).AND.AMN3SS.LT.(AMW+AML1SS)) THEN
 692         TM1=G**4*COS(THETAL)**2*SSXINT(AML1SS,SSSNWS,EMAX)
 693         WID=2*PI**2*TM1/3./2./AMN3SS/(2.*PI)**5
 694         CALL SSSAVE(ISNTL,WID,ISTAU1,IDNM,-IDMU,0,0)
 695       END IF
 696       IF (AMN3SS.GT.(AML1SS+AMUP+AMDN).AND.AMN3SS.LT.(AMW+AML1SS)) THEN
 697         TM1=G**4*COS(THETAL)**2*SSXINT(AML1SS,SSSNWS,EMAX)
 698         WID=2*PI**2*TM1/3./2./AMN3SS/(2.*PI)**5
 699         CALL SSSAVE(ISNTL,3*WID,ISTAU1,IDUP,-IDDN,0,0)
 700       END IF
 701       IF (AMN3SS.GT.(AML1SS+AMCH+AMST).AND.AMN3SS.LT.(AMW+AML1SS)) THEN
 702         TM1=G**4*COS(THETAL)**2*SSXINT(AML1SS,SSSNWS,EMAX)
 703         WID=2*PI**2*TM1/3./2./AMN3SS/(2.*PI)**5
 704         CALL SSSAVE(ISNTL,3*WID,ISTAU1,IDCH,-IDST,0,0)
 705       END IF
 706       TMP(4)=AMN2SS
 707       EMAX=(AMMLSS**2+AMN2SS**2)/2./AMN2SS
 708       IF (AMN2SS.GT.(AMMLSS+AME).AND.AMN2SS.LT.(AMW+AML1SS)) THEN
 709         TM1=G**4*COS(THETAL)**2*SSXINT(AMMLSS,SSSNWS,EMAX)
 710         WID=2*PI**2*TM1/3./2./AMN2SS/(2.*PI)**5
 711         CALL SSSAVE(ISNML,WID,ISMUL,IDNE,-IDE,0,0)
 712       END IF
 713       IF (AMN2SS.GT.(AMMLSS+AMMU).AND.AMN2SS.LT.(AMW+AML1SS)) THEN
 714         TM1=G**4*COS(THETAL)**2*SSXINT(AMMLSS,SSSNWS,EMAX)
 715         WID=2*PI**2*TM1/3./2./AMN2SS/(2.*PI)**5
 716         CALL SSSAVE(ISNML,WID,ISMUL,IDNM,-IDMU,0,0)
 717       END IF
 718       IF (AMN2SS.GT.(AMMLSS+AMUP+AMDN).AND.AMN2SS.LT.(AMW+AML1SS)) THEN
 719         TM1=G**4*COS(THETAL)**2*SSXINT(AMMLSS,SSSNWS,EMAX)
 720         WID=2*PI**2*TM1/3./2./AMN2SS/(2.*PI)**5
 721         CALL SSSAVE(ISNML,3*WID,ISMUL,IDUP,-IDDN,0,0)
 722       END IF
 723       IF (AMN2SS.GT.(AMMLSS+AMCH+AMST).AND.AMN2SS.LT.(AMW+AML1SS)) THEN
 724         TM1=G**4*COS(THETAL)**2*SSXINT(AMMLSS,SSSNWS,EMAX)
 725         WID=2*PI**2*TM1/3./2./AMN2SS/(2.*PI)**5
 726         CALL SSSAVE(ISNML,3*WID,ISMUL,IDCH,-IDST,0,0)
 727       END IF
 728       TMP(4)=AMN1SS
 729       EMAX=(AMELSS**2+AMN1SS**2)/2./AMN1SS
 730       IF (AMN1SS.GT.(AMELSS+AME).AND.AMN1SS.LT.(AMW+AML1SS)) THEN
 731         TM1=G**4*COS(THETAL)**2*SSXINT(AMELSS,SSSNWS,EMAX)
 732         WID=2*PI**2*TM1/3./2./AMN1SS/(2.*PI)**5
 733         CALL SSSAVE(ISNEL,WID,ISEL,IDNE,-IDE,0,0)
 734       END IF
 735       IF (AMN1SS.GT.(AMELSS+AMMU).AND.AMN1SS.LT.(AMW+AML1SS)) THEN
 736         TM1=G**4*COS(THETAL)**2*SSXINT(AMELSS,SSSNWS,EMAX)
 737         WID=2*PI**2*TM1/3./2./AMN1SS/(2.*PI)**5
 738         CALL SSSAVE(ISNEL,WID,ISEL,IDNM,-IDMU,0,0)
 739       END IF
 740       IF (AMN1SS.GT.(AMELSS+AMUP+AMDN).AND.AMN1SS.LT.(AMW+AML1SS)) THEN
 741         TM1=G**4*COS(THETAL)**2*SSXINT(AMELSS,SSSNWS,EMAX)
 742         WID=2*PI**2*TM1/3./2./AMN1SS/(2.*PI)**5
 743         CALL SSSAVE(ISNEL,3*WID,ISEL,IDUP,-IDDN,0,0)
 744       END IF
 745       IF (AMN1SS.GT.(AMELSS+AMCH+AMST).AND.AMN1SS.LT.(AMW+AML1SS)) THEN
 746         TM1=G**4*COS(THETAL)**2*SSXINT(AMELSS,SSSNWS,EMAX)
 747         WID=2*PI**2*TM1/3./2./AMN1SS/(2.*PI)**5
 748         CALL SSSAVE(ISNEL,3*WID,ISEL,IDCH,-IDST,0,0)
 749       END IF
 750 C-----COMPUTE ER AND MUR--> STAU_1+NEUTRINOS DECAYS --------------
 751       IF (AMERSS.LT.AMELSS.AND.AMERSS.GT.AML1SS) THEN
 752         TMP(1)=MW1
 753         TMP(2)=MW1
 754         TMP(3)=AML1SS
 755         TMP(4)=AMERSS
 756         IF (AMERSS.LT.MW1) THEN
 757         TM1=AML1**2*AME1**2*SSXINT(AML1SS**2,SSL1ST,AMERSS**2)
 758         TMP(2)=MW2
 759         TM2=2*AML1*AME1*APL1*APE1*SSXINT(AML1SS**2,SSL1ST,AMERSS**2)
 760         ELSE
 761         TM1=0.
 762         TM2=0.
 763         END IF
 764         IF (AMERSS.LT.MW2) THEN
 765         TMP(1)=MW2
 766         TMP(2)=MW2
 767         TM3=APL1**2*APE1**2*SSXINT(AML1SS**2,SSL1ST,AMERSS**2)
 768         ELSE
 769         TM3=0.
 770         END IF
 771         WID=PI**2*(TM1+TM2+TM3)/8./2./AMERSS/(2.*PI)**5
 772         CALL SSSAVE(ISER,WID,ISTAU1,IDNE,-IDNT,0,0)
 773       END IF
 774       IF (AMMRSS.LT.AMMLSS.AND.AMMRSS.GT.AML1SS) THEN
 775         AMM1=-G*SIN(GAMMAL)*COS(THETAM)+FM*COS(GAMMAL)*SIN(THETAM)
 776         AML1=-G*SIN(GAMMAL)*COS(THETAL)+FL*COS(GAMMAL)*SIN(THETAL)
 777         APM1=-G*THX*COS(GAMMAL)*COS(THETAM)-FM*THX*SIN(GAMMAL)*
 778      $       SIN(THETAM)
 779         APL1=-G*THX*COS(GAMMAL)*COS(THETAL)-FL*THX*SIN(GAMMAL)*
 780      $       SIN(THETAL)
 781         TMP(1)=MW1
 782         TMP(2)=MW1
 783         TMP(3)=AML1SS
 784         TMP(4)=AMMRSS
 785         IF (AMMRSS.LT.MW1) THEN
 786         TM1=AML1**2*AMM1**2*SSXINT(AML1SS**2,SSL1ST,AMMRSS**2)
 787         TMP(2)=MW2
 788         TM2=2*AML1*AMM1*APL1*APM1*SSXINT(AML1SS**2,SSL1ST,AMMRSS**2)
 789         ELSE
 790         TM1=0.
 791         TM2=0.
 792         END IF
 793         IF (AMMRSS.LT.MW2) THEN
 794         TMP(1)=MW2
 795         TMP(2)=MW2
 796         TM3=APL1**2*APM1**2*SSXINT(AML1SS**2,SSL1ST,AMMRSS**2)
 797         ELSE
 798         TM3=0.
 799         END IF
 800         WID=PI**2*(TM1+TM2+TM3)/8./2./AMMRSS/(2.*PI)**5
 801         CALL SSSAVE(ISMUR,WID,ISTAU1,IDNM,-IDNT,0,0)
 802       END IF
 803 C-----COMPUTE SNU_E --> E +STAU_1BAR +NU_TAU DECAYS --------------
 804       AME1=-G*SIN(GAMMAL)
 805       AML1=-G*SIN(GAMMAL)*COS(THETAL)+FL*COS(GAMMAL)*SIN(THETAL)
 806       APE1=-G*THX*COS(GAMMAL)
 807       APL1=-G*THX*COS(GAMMAL)*COS(THETAL)-FL*THX*SIN(GAMMAL)*
 808      $       SIN(THETAL)
 809       IF (AMN1SS.GT.(AME+AML1SS)) THEN
 810         BME1=-FE*COS(GAMMAL)
 811         BPE1=FE*THX*SIN(GAMMAL)
 812         TMP(1)=MW1
 813         TMP(2)=MW1
 814         TMP(3)=AML1SS
 815         TMP(4)=AMN1SS
 816         IF (AMN1SS.LT.MW1) THEN
 817         TM1=AML1**2*(AME1**2*MW1**2*SSXINT(AML1SS**2,SSN1ST,AMN1SS**2)
 818      $+BME1**2*SSXINT(AML1SS**2,SSL1ST,AMN1SS**2))
 819         TMP(2)=MW2
 820         TM2=2*AML1*APL1*(SNW1*SNW2*AME1*APE1*
 821      $ SSXINT(AML1SS**2,SSN1ST,AMN1SS**2)+BME1*BPE1*
 822      $ SSXINT(AML1SS**2,SSL1ST,AMN1SS**2))
 823         ELSE
 824         TM1=0.
 825         TM2=0.
 826         END IF
 827         TMP(1)=MW2
 828         TMP(2)=MW2
 829         IF (AMN1SS.LT.MW2) THEN
 830         TM3=APL1**2*(APE1**2*MW2**2*SSXINT(AML1SS**2,SSN1ST,AMN1SS**2)
 831      $+BPE1**2*SSXINT(AML1SS**2,SSL1ST,AMN1SS**2))
 832         ELSE
 833         TM3=0.
 834         END IF
 835         WID=PI**2*(TM1+TM2+TM3)/8./2./AMN1SS/(2.*PI)**5
 836         CALL SSSAVE(ISNEL,WID,-ISTAU1,IDE,IDNT,0,0)
 837       END IF
 838 C-----COMPUTE SNU_M --> MU +STAU_1BAR +NU_TAU DECAYS --------------
 839       IF (AMN2SS.GT.(AMMU+AML1SS)) THEN
 840         BME1=-FM*COS(GAMMAL)
 841         BPE1=FM*THX*SIN(GAMMAL)
 842         TMP(1)=MW1
 843         TMP(2)=MW1
 844         TMP(3)=AML1SS
 845         TMP(4)=AMN2SS
 846         IF (AMN2SS.LT.MW1) THEN
 847         TM1=AML1**2*(AME1**2*MW1**2*SSXINT(AML1SS**2,SSN1ST,AMN2SS**2)
 848      $+BME1**2*SSXINT(AML1SS**2,SSL1ST,AMN2SS**2))
 849         TMP(2)=MW2
 850         TM2=2*AML1*APL1*(SNW1*SNW2*AME1*APE1*
 851      $ SSXINT(AML1SS**2,SSN1ST,AMN2SS**2)+BME1*BPE1*
 852      $ SSXINT(AML1SS**2,SSL1ST,AMN2SS**2))
 853         ELSE
 854         TM1=0.
 855         TM2=0.
 856         END IF
 857         TMP(1)=MW2
 858         TMP(2)=MW2
 859         IF (AMN2SS.LT.MW2) THEN
 860         TM3=APL1**2*(APE1**2*MW2**2*SSXINT(AML1SS**2,SSN1ST,AMN2SS**2)
 861      $+BPE1**2*SSXINT(AML1SS**2,SSL1ST,AMN2SS**2))
 862         ELSE
 863         TM3=0.
 864         END IF
 865         WID=PI**2*(TM1+TM2+TM3)/8./2./AMN2SS/(2.*PI)**5
 866         CALL SSSAVE(ISNML,WID,-ISTAU1,IDMU,IDNT,0,0)
 867       END IF
 868 C
 869 C
 870 C     decay to neutral higgs bosons
 871 C
 872       IF (AML2SS.GT.(AMHL+AML1SS)) THEN
 873         BH=G*AMW*SIN(BETA-ALFAH)*(-1.+3*TN2THW)*SINL*COSL/2.+G*
 874      $  AMTAU*(TWOM1*COSA+AAL*SINA)*COS(2*THETAL)/2./AMW/COS(BETA)
 875         WID=BH**2*SQRT(SSXLAM(AML2SS**2,AMHL**2,AML1SS**2))/
 876      $      16./PI/AML2SS**3
 877         CALL SSSAVE(ISTAU2,WID,ISHL,ISTAU1,0,0,0)
 878       ENDIF
 879 C
 880       IF (AML2SS.GT.(AMHA+AML1SS)) THEN
 881         BH=G*AMTAU*(TWOM1-AAL*TANB)/2./AMW
 882         WID=BH**2*SQRT(SSXLAM(AML2SS**2,AMHA**2,AML1SS**2))/
 883      $      16./PI/AML2SS**3
 884         CALL SSSAVE(ISTAU2,WID,ISHA,ISTAU1,0,0,0)
 885       ENDIF
 886 C
 887       IF (AML2SS.GT.(AMHH+AML1SS)) THEN
 888         BH=-G*AMW*COS(BETA-ALFAH)*(-1.+3*TN2THW)*SINL*COSL/2.+G*
 889      $  AMTAU*(-TWOM1*SINA+AAL*COSA)*COS(2*THETAL)/2./AMW/COS(BETA)
 890         WID=BH**2*SQRT(SSXLAM(AML2SS**2,AMHH**2,AML1SS**2))/
 891      $      16./PI/AML2SS**3
 892         CALL SSSAVE(ISTAU2,WID,ISHH,ISTAU1,0,0,0)
 893       ENDIF
 894 C
 895 C     stau_i -> H^- stau_j
 896 C
 897       IF (AMN3SS.GT.(AML1SS+AMHC)) THEN
 898         A=G/SR2/AMW*((AMTAU**2*TANB-AMW**2*SIN(2*BETA))*COSL-
 899      $    AMTAU*(TWOM1-AAL*TANB)*SINL)
 900         WID=A*A*SQRT(SSXLAM(AMN3SS**2,AML1SS**2,AMHC**2))/
 901      $      16./PI/AMN3SS**3
 902         CALL SSSAVE(ISNTL,WID,ISHC,ISTAU1,0,0,0)
 903       END IF
 904 C
 905       IF (AMN3SS.GT.(AML2SS+AMHC)) THEN
 906         A=G/SR2/AMW*((AMTAU**2*TANB-AMW**2*SIN(2*BETA))*SINL+
 907      $    AMTAU*(TWOM1-AAL*TANB)*COSL)
 908         WID=A*A*SQRT(SSXLAM(AMN3SS**2,AML2SS**2,AMHC**2))/
 909      $      16./PI/AMN3SS**3
 910         CALL SSSAVE(ISNTL,WID,ISHC,ISTAU2,0,0,0)
 911       END IF
 912 C
 913       IF (AML1SS.GT.(AMN3SS+AMHC)) THEN
 914         A=G/SR2/AMW*((AMTAU**2*TANB-AMW**2*SIN(2*BETA))*COSL-
 915      $    AMTAU*(TWOM1-AAL*TANB)*SINL)
 916         WID=A*A*SQRT(SSXLAM(AML1SS**2,AMN3SS**2,AMHC**2))/
 917      $      16./PI/AML1SS**3
 918         CALL SSSAVE(ISTAU1,WID,-ISHC,ISNTL,0,0,0)
 919       END IF
 920 C
 921       IF (AML2SS.GT.(AMN3SS+AMHC)) THEN
 922         A=G/SR2/AMW*((AMTAU**2*TANB-AMW**2*SIN(2*BETA))*SINL+
 923      $    AMTAU*(TWOM1-AAL*TANB)*COSL)
 924         WID=A*A*SQRT(SSXLAM(AML2SS**2,AMN3SS**2,AMHC**2))/
 925      $      16./PI/AML2SS**3
 926         CALL SSSAVE(ISTAU2,WID,-ISHC,ISNTL,0,0,0)
 927       END IF
 928 C
 929 C     Gauge mediated SUSY breaking model decays to gravitino
 930       IF (AMERSS.GT.(AME+AMGVSS)) THEN
 931         WID=(AMERSS**2-AME**2)**4/48./PI/AMERSS**3/(AMPL*AMGVSS)**2
 932         CALL SSSAVE(ISER,WID,IDE,91,0,0,0)
 933       END IF
 934       IF (AMMRSS.GT.(AMMU+AMGVSS)) THEN
 935         WID=(AMMRSS**2-AMMU**2)**4/48./PI/AMMRSS**3/(AMPL*AMGVSS)**2
 936         CALL SSSAVE(ISMUR,WID,IDMU,91,0,0,0)
 937       END IF
 938       IF (AML1SS.GT.(AMTAU+AMGVSS)) THEN
 939         WID=(AML1SS**2-AMTAU**2)**4/48./PI/AML1SS**3/(AMPL*AMGVSS)**2
 940         CALL SSSAVE(ISTAU1,WID,IDTAU,91,0,0,0)
 941       END IF
 942       IF (AMELSS.GT.(AME+AMGVSS)) THEN
 943         WID=(AMELSS**2-AME**2)**4/48./PI/AMELSS**3/(AMPL*AMGVSS)**2
 944         CALL SSSAVE(ISEL,WID,IDE,91,0,0,0)
 945       END IF
 946       IF (AMMLSS.GT.(AMMU+AMGVSS)) THEN
 947         WID=(AMMLSS**2-AMMU**2)**4/48./PI/AMMLSS**3/(AMPL*AMGVSS)**2
 948         CALL SSSAVE(ISMUL,WID,IDMU,91,0,0,0)
 949       END IF
 950       IF (AMN1SS.GT.AMGVSS) THEN
 951         WID=AMN1SS**5/48./PI/(AMPL*AMGVSS)**2
 952         CALL SSSAVE(ISNEL,WID,IDNE,91,0,0,0)
 953       END IF
 954       IF (AMN2SS.GT.AMGVSS) THEN
 955         WID=AMN2SS**5/48./PI/(AMPL*AMGVSS)**2
 956         CALL SSSAVE(ISNML,WID,IDNM,91,0,0,0)
 957       END IF
 958       IF (AMN3SS.GT.AMGVSS) THEN
 959         WID=AMN3SS**5/48./PI/(AMPL*AMGVSS)**2
 960         CALL SSSAVE(ISNTL,WID,IDNT,91,0,0,0)
 961       END IF
 962 C
 963 C          Normalize branching ratios
 964 C
 965       CALL SSNORM(ISEL)
 966       CALL SSNORM(ISMUL)
 967       CALL SSNORM(ISTAU1)
 968       CALL SSNORM(ISER)
 969       CALL SSNORM(ISMUR)
 970       CALL SSNORM(ISTAU2)
 971       CALL SSNORM(ISNEL)
 972       CALL SSNORM(ISNML)
 973       CALL SSNORM(ISNTL)
 974 C
 975        RETURN
 976        END