]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | REAL FUNCTION SSGX6(ET) | |
3 | C----------------------------------------------------------------------- | |
4 | C SSGLBF: glss -> ziss + tp + tb | |
5 | C Baer's XT6 - ZETA- eq. a.6.f 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,T1,T2,DXT6 | |
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 | T1=DMG**2+DMT**2-2*DMG*DET-DMT1**2 | |
37 | T2=DMG**2+DMT**2-2*DMG*DET-DMT2**2 | |
38 | DXT6=PI**2*(EMAX-EMIN)/T1/T2 | |
39 | SSGX6=DXT6 | |
40 | RETURN | |
41 | END |