ISAJET/isasusy/ssgwt6.F

   1 #include "isajet/pilot.h"
   2         REAL FUNCTION SSGWT6(E)
   3 C-----------------------------------------------------------------------
   4 C          SSGLBF: glss -> wiss + tp + bb
   5 C          Baer's FTBW23/Drees' G_6
   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 ET,MWI,MG,MT,MSB,MST,PT,EBMX
  15         DOUBLE PRECISION DEN,T1,DFTBW,EBMN,XX,XL
  16         DOUBLE PRECISION MB,Z,RDL,DENO,R1,R2,R3,SSDLAM
  17 C
  18         ET=E
  19         MWI=TMP(1)
  20         MG=TMP(2)
  21         MT=TMP(3)
  22         MSB=TMP(4)
  23         MST=TMP(6)
  24         MB=AMBT
  25 C
  26 C            Rewrite PT=DSQRT(ET**2-MT**2)
  27         PT=DSQRT((ET-MT)*(ET+MT))
  28         Z=(MG**2+MT**2-2*MG*ET+MB**2-MWI**2)/2.
  29         R1=1.D0+MT**2/MG**2-2.D0*ET/MG
  30         R2=MB**2/MG**2
  31         R3=MWI**2/MG**2
  32         RDL=DSQRT(DMAX1(0.D0,SSDLAM(R1,R2,R3)))
  33         DENO=MG**2+MT**2-2*ET*MG
  34         EBMX=(2*Z*(1.D0-ET/MG)+PT*MG*RDL)*MG/2.D0/DENO
  35         EBMN=(2*Z*(1.D0-ET/MG)-PT*MG*RDL)*MG/2.D0/DENO
  36         XX=(MSB**2+2*MG*EBMX-MG**2)/(MSB**2+2*MG*EBMN-MG**2)
  37         XL=DLOG(XX)
  38         DEN=MG**2-2*MG*ET+MT**2-MST**2
  39         T1=(MG*(MG**2+MT**2-2*MG*ET-MWI**2)-(MSB**2-MG**2)*
  40      $      (2*ET*MG-MT**2-MG**2)/MG)*XL+2*(2*ET*MG-MT**2-MG**2)*
  41      $      (EBMX-EBMN)
  42         DFTBW=.5D0*T1/DEN
  43         SSGWT6=DFTBW
  44         RETURN
  45         END