13 data ( px(i), i= 1 , np )/
14 1 0.499, 0.596, 0.769, 0.919, 0.967,
15 2 1.071, 1.233, 1.492, 2.001, 2.423,
16 3 3.087, 3.883, 5.140, 7.929, 13.895,
17 4 20.630, 32.647, 51.663, 80.719, 124.520,
18 5 207.359, 392.263, 811.319, 2005.980, 4959.750,
22 data ( mb(i), i= 1 , np )/
23 1 10.185, 11.888, 14.742, 17.207, 18.529,
24 2 19.290, 19.418, 18.897, 18.389, 17.896,
25 3 17.652, 17.647, 17.642, 17.752, 17.860,
26 4 17.852, 18.084, 18.319, 18.683, 19.055,
27 5 19.828, 20.910, 21.901, 23.088, 24.670,
34 elseif( p .gt. 0.6)
then 55 real(8),
parameter:: Pnorm=50.0
57 real(8):: xspel, xspt, xspit, PPP
61 logical,
save:: first = .
true.
65 data ( px(i), i= 1 , np )/
66 1 0.499, 0.619, 0.759, 0.908, 0.992,
67 2 1.142, 1.348, 1.695, 2.053, 2.717,
68 3 4.086, 5.990, 8.451, 11.331, 15.387,
69 4 21.710, 33.066, 49.724, 75.734, 112.445,
70 5 146.967, 197.049, 271.020, 445.603, 661.604,
71 6 982.309, 1678.060, 2903.360, 5023.390, 7749.220
73 data ( mb(i), i= 1 , np )/
74 1 5.160,5.942,6.661,7.318,7.567,
75 2 7.316,6.884,6.263,5.510,4.718,
76 3 4.150,3.505,3.146,2.901,2.711,
77 4 2.550,2.432,2.366,2.365,2.461,
78 5 2.477,2.561,2.648,2.793,2.926,
79 6 3.025,3.169,3.342,3.408,3.668
85 if( p .gt. pnorm)
then 101 xs = xspit*xspel/xspt
103 elseif(p .gt. px(1))
then 113 real(8),
intent(in)::p
114 real(8),
intent(out)::xs
119 xs =max( txs - exs, 0.
d0)
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
real(8), parameter, public m
subroutine cpptotxs(p, xs)
real(8) function, public csigmat(a, b, p)
block data cblkElemag data *AnihiE ! Eposi< 1 TeV, anihilation considered *X0/365.667/, ! radiation length of air in kg/m2 *Ecrit/81.e-3/, ! critical energy of air in GeV *MaxComptonE/1./, ! compton is considered below 1 GeV *MaxPhotoE/1.e-3/, ! above this, PhotoElectric effect neg. *MinPhotoProdE/153.e-3/, ! below 153 MeV, no gp --> hadrons ! scattering const not MeV *Knockon true
subroutine kpolintplogxyfe(xa, xstep, ya, ystep, nt, m, logxy, x, y, error)
subroutine ckpntotxs(p, xs)
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
subroutine cppelaxs(p, xs)
subroutine ckpnelaxs(p, xs)
subroutine ckpninelaxs(p, xs)
subroutine kpolintpfe(xa, xstep, ya, ystep, nt, m, x, y, error)
block data cblkElemag data *AnihiE ! Eposi< 1 TeV, anihilation considered *X0/365.667/, ! radiation length of air in kg/m2 *Ecrit/81.e-3/, ! critical energy of air in GeV *MaxComptonE/1./, ! compton is considered below 1 GeV *MaxPhotoE/1.e-3/, ! above this, PhotoElectric effect neg. *MinPhotoProdE/153.e-3/, ! below 153 MeV, no gp --> hadrons ! scattering const not MeV *Knockon ! knockon is considered Obsolete *PhotoProd false