ISAJET/isasusy/ssgx8.F

   1 #include "isajet/pilot.h"
   2         REAL FUNCTION SSGX8(ET)
   3 C-----------------------------------------------------------------------
   4 C          SSGLBF: glss -> ziss + tp + tb
   5 C          Baer's XT8 - Y- eq. a.6.h 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,DXT8
  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-DMT2**2
  37         BOT=DMG**2-2*DMG*EMIN+DMT**2-DMT2**2
  38         DXT8=((EMAX-EMIN)*(DMG**2-2*DMG*DET+DMT**2)+
  39      $        (DMZ**2*DMG**2-DMT2**2*DMG**2+2*DMT2**2*DMG*DET+DMT**4
  40      $         -DMT2**2*DMT**2)/2.D0/DMG*DLOG(TOP/BOT))*PI**2/2.D0/
  41      $         (DMG**2+DMT**2-2*DMG*DET-DMT1**2)
  42         SSGX8=DXT8
  43         RETURN
  44         END