1 #include "isajet/pilot.h"
2 REAL FUNCTION SSGWT4(E)
3 C-----------------------------------------------------------------------
4 C SSGLBF: glss -> wiss + tp + bb
5 C Baer's FTBW12; EQ. A.3D OF BTW, MODIFIED FOR MB=/0
7 C-----------------------------------------------------------------------
8 #if defined(CERNLIB_IMPNONE)
11 #include "isajet/sssm.inc"
12 #include "isajet/sspar.inc"
13 #include "isajet/sstmp.inc"
15 DOUBLE PRECISION ET,MWI,MG,MT,MSB,MST,PT,EBMX
16 DOUBLE PRECISION BOT,TOP,DFTBW,EBMN,XX
17 DOUBLE PRECISION MB,Z,RDL,DEN,R1,R2,R3,SSDLAM,SN
28 C Rewrite PT=DSQRT(ET**2-MT**2)
29 PT=DSQRT((ET-MT)*(ET+MT))
30 Z=(MG**2+MT**2-2*MG*ET+MB**2-MWI**2)/2.
31 R1=1.D0+MT**2/MG**2-2.D0*ET/MG
34 RDL=DSQRT(DMAX1(0.D0,SSDLAM(R1,R2,R3)))
35 DEN=MG**2+MT**2-2*ET*MG
36 EBMX=(2*Z*(1.D0-ET/MG)+PT*MG*RDL)*MG/2.D0/DEN
37 EBMN=(2*Z*(1.D0-ET/MG)-PT*MG*RDL)*MG/2.D0/DEN
38 XX=(MSB**2+2*MG*EBMX-MG**2)/(MSB**2+2*MG*EBMN-MG**2)
39 TOP=EBMX-EBMN-(MSB**2-2*MG*ET+MT**2-MWI**2)*DLOG(XX)/2.D0/MG
40 BOT=(MG**2+MT**2-2*MG*ET-MST**2)
41 DFTBW=-SN*MG*MWI*TOP/BOT