ISAJET/isasusy/ssgx9.F

   1 #include "isajet/pilot.h"
   2         REAL FUNCTION SSGX9(ET)
   3 C-----------------------------------------------------------------------
   4 C          SSGLBF: glss -> ziss + tp + tb
   5 C          Baer's XT9 - XI'- eq. a.6.i 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,BOT,DXT9
  16         DOUBLE PRECISION XT,MUT,MUZ,XMIN,XMAX,EMIN,EMAX,SSDLAM
  17         DOUBLE PRECISION PI
  18         DATA PI/3.14159265D0/
  19         DET=ET
  20         DMG=TMP(1)
  21         DMT=TMP(2)
  22         DMZ=TMP(3)
  23         DMT1=TMP(4)
  24         DMT2=TMP(5)
  25         XT=2*DET/DMG
  26         MUT=(DMT/DMG)**2
  27         MUZ=(DMZ/DMG)**2
  28         XMIN=((2.D0-XT)*(1.D0+2*MUT-MUZ-XT)-DSQRT(DMAX1(0.D0,
  29      $   (XT**2-4*MUT)*SSDLAM((1.D0+MUT-XT),MUT,MUZ))))
  30      $   /2.D0/(1.D0-XT+MUT)
  31         XMAX=((2.D0-XT)*(1.D0+2*MUT-MUZ-XT)+DSQRT(DMAX1(0.D0,
  32      $   (XT**2-4*MUT)*SSDLAM((1.D0+MUT-XT),MUT,MUZ))))
  33      $   /2.D0/(1.D0-XT+MUT)
  34         EMIN=XMIN*DMG/2.D0
  35         EMAX=XMAX*DMG/2.D0
  36         TOP=DMG**2-2*DMG*EMAX+DMT**2-DMT1**2
  37         BOT=DMG**2-2*DMG*EMIN+DMT**2-DMT1**2
  38         DXT9=-PI**2*DET*DLOG(TOP/BOT)/2.D0/(DMG**2-2*DMG*DET+DMT**2
  39      $         -DMT2**2)
  40         SSGX9=DXT9
  41         RETURN
  42         END