ISAJET/isasusy/ssgx1.F

   1 #include "isajet/pilot.h"
   2         REAL FUNCTION SSGX1(ET)
   3 C-----------------------------------------------------------------------
   4 C          SSGLBF: glss -> ziss + tp + tb
   5 C          Baer's XT1 - PSI- eq. a.6.a of prd45,142 (1992)
   6 C          Modified for t_1 and t_2 eigenstates
   7 C-----------------------------------------------------------------------
   8 #if defined(CERNLIB_IMPNONE)
   9       IMPLICIT NONE
  10 #endif
  11 #include "isajet/sssm.inc"
  12 #include "isajet/sspar.inc"
  13 #include "isajet/sstmp.inc"
  14         REAL ET
  15         DOUBLE PRECISION DET,DMG,DMT,DMZ,DMT1,DMT2,TOP
  16         DOUBLE PRECISION BOT,PT,DXT1,SSDLAM,PI,TOPS
  17         DATA PI/3.14159265D0/
  18         DET=ET
  19         DMG=TMP(1)
  20         DMT=TMP(2)
  21         DMZ=TMP(3)
  22         DMT1=TMP(4)
  23         DMT2=TMP(5)
  24         PT=DSQRT(DET**2-DMT**2)
  25         TOPS=SSDLAM(DMG**2+DMT**2-2*DMG*DET,DMT**2,DMZ**2)
  26         TOP=DSQRT(DMAX1(0.D0,TOPS))
  27         BOT=(DMG**2+DMT**2-2*DMG*DET-DMT1**2)*
  28      $   (DMG**2+DMT**2-2*DMG*DET-DMT2**2)
  29         DXT1=PI**2*DMG*DET*PT*(DMG**2-DMZ**2-2*DMG*DET)*TOP/BOT/
  30      $       (DMG**2-2*DMG*DET+DMT**2)
  31         SSGX1=DXT1
  32         RETURN
  33         END