ISAJET/isasusy/ssgwt2.F

   1 #include "isajet/pilot.h"
   2         REAL FUNCTION SSGWT2(E)
   3 C-----------------------------------------------------------------------
   4 C          SSGLBF: glss -> wiss + tp + bb
   5 C          Drees' G_2
   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,DFTBW,EB,TOP,BOT,MSB
  15         DOUBLE PRECISION SSDLAM,TERM
  16 C
  17         EB=E
  18         MWI=TMP(1)
  19         MG=TMP(2)
  20         MT=TMP(3)
  21         MSB=TMP(4)
  22 C
  23         TOP=EB**2*(MG**2+AMBT**2-2*MG*EB-MWI**2-MT**2)
  24         BOT=(MG**2+AMBT**2-2*MG*EB-MSB**2)**2*
  25      $       (MG**2+AMBT**2-2*EB*MG)
  26         TERM=SSDLAM((MG**2+AMBT**2-2*MG*EB),MWI**2,MT**2)
  27         IF(TERM.GT.0.D0) THEN
  28           DFTBW=MG*TOP/BOT*DSQRT(TERM)
  29         ELSE
  30           DFTBW=0.D0
  31         ENDIF
  32         SSGWT2=DFTBW
  33         RETURN
  34         END