]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | REAL FUNCTION SSGZT(E) | |
3 | C----------------------------------------------------------------------- | |
4 | C SSGLBF: glss -> ziss + tp + tb | |
5 | C Baer's TOPINT | |
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 | C | |
14 | REAL E | |
15 | DOUBLE PRECISION MG,MT,MS,MZ,SSDLAM,P,PSI,XLOG,PHI,SN,C1,C2 | |
16 | DOUBLE PRECISION XT,MUT,MUZ,XMIN,XMAX,EMIN,EMAX | |
17 | C Convert to double precision | |
18 | MG=AMGLSS | |
19 | C1=TMP(1) | |
20 | C2=TMP(2) | |
21 | MS=TMP(3) | |
22 | MT=TMP(4) | |
23 | MZ=TMP(5) | |
24 | SN=TMP(6) | |
25 | C | |
26 | XT=2*E/MG | |
27 | MUT=(MT/MG)**2 | |
28 | MUZ=(MZ/MG)**2 | |
29 | XMIN=((2.D0-XT)*(1.D0+2*MUT-MUZ-XT)-DSQRT((XT**2-4*MUT)* | |
30 | $SSDLAM((1.D0+MUT-XT),MUT,MUZ)))/2.D0/(1.D0-XT+MUT) | |
31 | XMAX=((2.D0-XT)*(1.D0+2*MUT-MUZ-XT)+DSQRT((XT**2-4*MUT)* | |
32 | $SSDLAM((1.D0+MUT-XT),MUT,MUZ)))/2.D0/(1.D0-XT+MUT) | |
33 | EMIN=XMIN*MG/2.D0 | |
34 | EMAX=XMAX*MG/2.D0 | |
35 | P=SQRT(E**2-MT**2) | |
36 | PSI=P*E*(MG**2-MZ**2-2*MG*E)* | |
37 | $DSQRT(SSDLAM((MG**2+MT**2-2*MG*E),MZ**2,MT**2))/MG/ | |
38 | $(MG**2+MT**2-2*MG*E)/(MG**2+MT**2-2*MG*E-MS**2)**2 | |
39 | XLOG=DLOG((MG**2+MT**2-2*MG*EMAX-MS**2)/ | |
40 | $(MG**2+MT**2-2*MG*EMIN-MS**2)) | |
41 | PHI=MZ*(-(EMAX-EMIN)-(2*E*MG+MZ**2-MT**2-MS**2)* | |
42 | $XLOG/2.D0/MG)/ | |
43 | $2.D0/MG/(MG**2+MT**2-MS**2-2*MG*E) | |
44 | SSGZT=C1*PSI+SN*C2*PHI | |
45 | RETURN | |
46 | END |