]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | REAL FUNCTION SSGX9(ET) | |
3 | C----------------------------------------------------------------------- | |
4 | C SSGLBF: glss -> ziss + tp + tb | |
5 | C Baer's XT9 - XI'- eq. a.6.i of prd45,142 (1992) | |
6 | C Modified for t_1 and t_2 eigenstates | |
7 | C----------------------------------------------------------------------- | |
8 | #if defined(CERNLIB_IMPNONE) | |
9 | IMPLICIT NONE | |
10 | #endif | |
11 | #include "isajet/sssm.inc" | |
12 | #include "isajet/sspar.inc" | |
13 | #include "isajet/sstmp.inc" | |
14 | REAL ET | |
15 | DOUBLE PRECISION DET,DMG,DMT,DMZ,DMT1,DMT2,TOP,BOT,DXT9 | |
16 | DOUBLE PRECISION XT,MUT,MUZ,XMIN,XMAX,EMIN,EMAX,SSDLAM | |
17 | DOUBLE PRECISION PI | |
18 | DATA PI/3.14159265D0/ | |
19 | DET=ET | |
20 | DMG=TMP(1) | |
21 | DMT=TMP(2) | |
22 | DMZ=TMP(3) | |
23 | DMT1=TMP(4) | |
24 | DMT2=TMP(5) | |
25 | XT=2*DET/DMG | |
26 | MUT=(DMT/DMG)**2 | |
27 | MUZ=(DMZ/DMG)**2 | |
28 | XMIN=((2.D0-XT)*(1.D0+2*MUT-MUZ-XT)-DSQRT(DMAX1(0.D0, | |
29 | $ (XT**2-4*MUT)*SSDLAM((1.D0+MUT-XT),MUT,MUZ)))) | |
30 | $ /2.D0/(1.D0-XT+MUT) | |
31 | XMAX=((2.D0-XT)*(1.D0+2*MUT-MUZ-XT)+DSQRT(DMAX1(0.D0, | |
32 | $ (XT**2-4*MUT)*SSDLAM((1.D0+MUT-XT),MUT,MUZ)))) | |
33 | $ /2.D0/(1.D0-XT+MUT) | |
34 | EMIN=XMIN*DMG/2.D0 | |
35 | EMAX=XMAX*DMG/2.D0 | |
36 | TOP=DMG**2-2*DMG*EMAX+DMT**2-DMT1**2 | |
37 | BOT=DMG**2-2*DMG*EMIN+DMT**2-DMT1**2 | |
38 | DXT9=-PI**2*DET*DLOG(TOP/BOT)/2.D0/(DMG**2-2*DMG*DET+DMT**2 | |
39 | $ -DMT2**2) | |
40 | SSGX9=DXT9 | |
41 | RETURN | |
42 | END |