]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | REAL FUNCTION SSGX10(ET) | |
3 | C----------------------------------------------------------------------- | |
4 | C SSGLBF: glss -> ziss + tp + tb | |
5 | C Baer's XT10 | |
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 MG,MT,MZ,MST1,MST2,ET | |
14 | DOUBLE PRECISION DET,DMG,DMT,DMZ,DMT1,DMT2,TOP,BOT,DXT10 | |
15 | DOUBLE PRECISION XT,MUT,MUZ,XMIN,XMAX,EMIN,EMAX,SSDLAM | |
16 | DOUBLE PRECISION PI | |
17 | DATA PI/3.14159265D0/ | |
18 | MG=TMP(1) | |
19 | MT=TMP(2) | |
20 | MZ=TMP(3) | |
21 | MST1=TMP(4) | |
22 | MST2=TMP(5) | |
23 | DET=ET | |
24 | DMG=TMP(1) | |
25 | DMT=TMP(2) | |
26 | DMZ=TMP(3) | |
27 | DMT1=TMP(4) | |
28 | DMT2=TMP(5) | |
29 | XT=2*ET/MG | |
30 | MUT=(MT/MG)**2 | |
31 | MUZ=(MZ/MG)**2 | |
32 | XMIN=((2.D0-XT)*(1.D0+2*MUT-MUZ-XT)-DSQRT(DMAX1(0.D0, | |
33 | $ (XT**2-4*MUT)*SSDLAM((1.D0+MUT-XT),MUT,MUZ)))) | |
34 | $ /2.D0/(1.D0-XT+MUT) | |
35 | XMAX=((2.D0-XT)*(1.D0+2*MUT-MUZ-XT)+DSQRT(DMAX1(0.D0, | |
36 | $ (XT**2-4*MUT)*SSDLAM((1.D0+MUT-XT),MUT,MUZ)))) | |
37 | $ /2.D0/(1.D0-XT+MUT) | |
38 | EMIN=XMIN*MG/2.D0 | |
39 | EMAX=XMAX*MG/2.D0 | |
40 | TOP=DMG**2-2*DMG*EMAX+DMT**2-DMT2**2 | |
41 | BOT=DMG**2-2*DMG*EMIN+DMT**2-DMT2**2 | |
42 | DXT10=(EMAX-EMIN-(DMG**2-2*DMG*DET+DMT2**2-DMT**2)/2.D0/DMG* | |
43 | $ DLOG(TOP/BOT))*PI**2/2.D0/(DMG**2+DMT**2-2*DMG*DET-DMT1**2) | |
44 | SSGX10=DXT10 | |
45 | RETURN | |
46 | END |