ISAJET/isasusy/ssgwt3.F

   1 #include "isajet/pilot.h"
   2         REAL FUNCTION SSGWT3(E)
   3 C-----------------------------------------------------------------------
   4 C          SSGLBF: glss -> wiss + tp + bb
   5 C          Drees' G_3
   6 C-----------------------------------------------------------------------
   7 #if defined(CERNLIB_IMPNONE)
   8       IMPLICIT NONE
   9 #endif
  10 #include "isajet/sssm.inc"
  11 #include "isajet/sspar.inc"
  12 #include "isajet/sstmp.inc"
  13         REAL E
  14         DOUBLE PRECISION MWI,MG,MT,EB,BOT,MSB
  15         DOUBLE PRECISION SSDLAM,TERM,SN,DFTBW
  16 C
  17         EB=E
  18         MWI=TMP(1)
  19         MG=TMP(2)
  20         MT=TMP(3)
  21         MSB=TMP(4)
  22         SN=TMP(8)
  23 C
  24
  25         BOT=(MG**2+AMBT**2-2*MG*EB-MSB**2)**2*
  26      $       (MG**2+AMBT**2-2*EB*MG)
  27         TERM=SSDLAM((MG**2+AMBT**2-2*MG*EB),MWI**2,MT**2)
  28         IF(TERM.GT.0.D0) THEN
  29           DFTBW=4*MG*MWI*MT*SN*EB*EB/BOT*DSQRT(TERM)
  30         ELSE
  31           DFTBW=0.D0
  32         ENDIF
  33         SSGWT3=DFTBW
  34         RETURN
  35         END