]>
Commit | Line | Data |
---|---|---|
0795afa3 | 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 |