ISAJET/isasusy/ssgx10.F

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