1 #include "isajet/pilot.h"
2 REAL FUNCTION SSGWT8(E)
3 C-----------------------------------------------------------------------
4 C SSGLBF: glss -> wiss + tp + bb
5 C Baer's FTBW13/ Drees' G_8
6 C-----------------------------------------------------------------------
7 #if defined(CERNLIB_IMPNONE)
10 #include "isajet/sssm.inc"
11 #include "isajet/sspar.inc"
12 #include "isajet/sstmp.inc"
14 DOUBLE PRECISION ET,MWI,MG,MT,MST1,MST2,PT,EBMX,EBMN
15 DOUBLE PRECISION TOP,BOT,DFTBW
16 DOUBLE PRECISION MB,Z,RDL,DEN,R1,R2,R3,SSDLAM
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
32 RDL=DSQRT(DMAX1(0.D0,SSDLAM(R1,R2,R3)))
33 DEN=MG**2+MT**2-2*ET*MG
34 EBMX=(2*Z*(1.D0-ET/MG)+PT*MG*RDL)*MG/2.D0/DEN
35 EBMN=(2*Z*(1.D0-ET/MG)-PT*MG*RDL)*MG/2.D0/DEN
36 TOP=(MG**2+MT**2-2*MG*ET-MWI**2)*(EBMX-EBMN)
37 BOT=(MG**2+MT**2-2*MG*ET-MST1**2)*
38 $ (MG**2+MT**2-2*MG*ET-MST2**2)