ISAJET/isasusy/sslrt1.F

   1 #include "isajet/pilot.h"
   2         REAL FUNCTION SSLRT1(SS)
   3 C-----------------------------------------------------------------------
   4 C          SSLRT1: l_R -> l+tau+stau_1
   5 C-----------------------------------------------------------------------
   6 #if defined(CERNLIB_IMPNONE)
   7       IMPLICIT NONE
   8 #endif
   9 #include "isajet/sssm.inc"
  10 #include "isajet/sspar.inc"
  11 #include "isajet/sstmp.inc"
  12         REAL SS
  13         DOUBLE PRECISION ETMX,ETMN,S,MT1,MT,MLR,BEZI,BEZJ,TM,
  14      ,AI,AJ,BI,BJ,MZI,MZJ,SNZI,SNZJ,XL,BK1,BK2,BK3,BK,WID,SSDLAM
  15         S=SS
  16         MT1=AML1SS
  17         MT=AMTAU
  18         MLR=TMP(1)
  19         BEZI=TMP(2)
  20         BEZJ=TMP(3)
  21         AI=TMP(4)
  22         AJ=TMP(5)
  23         BI=TMP(6)
  24         BJ=TMP(7)
  25         MZI=ABS(TMP(8))
  26         MZJ=ABS(TMP(9))
  27         SNZI=SIGN(1.0,TMP(8))
  28         SNZJ=SIGN(1.0,TMP(9))
  29         TM=SSDLAM(S,MT**2,MT1**2)
  30         XL=DSQRT(MAX(0.D0,TM))
  31         ETMN=(S+MT**2-MT1**2-XL*(MLR**2-S)/(MLR**2+S))*(MLR**2+S)/
  32      ,       (2*S)/(2*MLR)
  33         ETMX=(S+MT**2-MT1**2+XL*(MLR**2-S)/(MLR**2+S))*(MLR**2+S)/
  34      ,       (2*S)/(2*MLR)
  35         BK1=-(ETMX-ETMN)*((ETMX+ETMN)*MLR*S-
  36      ,       (S+MT**2-MT1**2)*MLR**2)/2.D0
  37         BK2=(ETMX-ETMN)*((ETMX+ETMN)*MLR-S-MT**2+MT1**2)/2.D0
  38         BK3=SNZJ*BI*AJ*MZJ+SNZI*BJ*AI*MZI
  39         BK=BI*BJ*BK1+AI*AJ*MZI*MZJ*SNZI*SNZJ*BK2+BK3*MT*(MLR**2-S)*
  40      ,     (ETMX-ETMN)/2.D0
  41         WID=BEZI*BEZJ*BK/(S-MZI**2)/(S-MZJ**2)
  42         SSLRT1=WID
  43         RETURN
  44         END