| 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 |
git://git.uio.no / u / mrichter / AliRoot.git / blame