ISAJET/isasusy/ssgwt1.F

   1 #include "isajet/pilot.h"
   2         REAL FUNCTION SSGWT1(E)
   3 C-----------------------------------------------------------------------
   4 C          SSGLBF: glss -> wiss + tp + bb
   5 C          Baer's FTBW1
   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,PT,ET,MST1,MST2,TOP,BOT
  15 C
  16         ET=E
  17         MWI=TMP(1)
  18         MG=TMP(2)
  19         MT=TMP(3)
  20         MST1=TMP(6)
  21         MST2=TMP(7)
  22 C
  23 C            Rewrite PT=DSQRT(ET**2-MT**2)
  24         PT=DSQRT((ET-MT)*(ET+MT))
  25         TOP=(MG**2+MT**2-2*MG*ET-MWI**2)**2*ET*PT
  26         BOT=(MG**2+MT**2-2*MG*ET-MST1**2)*(MG**2+MT**2-2*MG*ET-
  27      $      MST2**2)*(MG**2+MT**2-2*ET*MG)
  28         DFTBW=MG*TOP/BOT
  29         SSGWT1=DFTBW
  30         RETURN
  31         END