]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | REAL FUNCTION SSGWT1(E) | |
3 | C----------------------------------------------------------------------- | |
4 | C SSGLBF: glss -> wiss + tp + bb | |
5 | C Baer's FTBW1 | |
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 MWI,MG,MT,DFTBW,PT,ET,MST1,MST2,TOP,BOT | |
15 | C | |
16 | ET=E | |
17 | MWI=TMP(1) | |
18 | MG=TMP(2) | |
19 | MT=TMP(3) | |
20 | MST1=TMP(6) | |
21 | MST2=TMP(7) | |
22 | C | |
23 | C Rewrite PT=DSQRT(ET**2-MT**2) | |
24 | PT=DSQRT((ET-MT)*(ET+MT)) | |
25 | TOP=(MG**2+MT**2-2*MG*ET-MWI**2)**2*ET*PT | |
26 | BOT=(MG**2+MT**2-2*MG*ET-MST1**2)*(MG**2+MT**2-2*MG*ET- | |
27 | $ MST2**2)*(MG**2+MT**2-2*ET*MG) | |
28 | DFTBW=MG*TOP/BOT | |
29 | SSGWT1=DFTBW | |
30 | RETURN | |
31 | END |