8 !
is not included here.
9 ! unit of field strength
is in
T (1 gauss = 10**-4
T)
14 real*8
x, ! in earth_center coordinate
22 endmap ! same
as x,
y,
z 25 real*8
n, ! north
com.
31 * v, ! vertical comp.(down
is +)
32 *
a ! deflection angle (deg. east
is +)
38 character*4
sys ! which
system.
'xyz',
'ned',
'hva'
dE dx *! Nuc Int sampling table e
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz sys
real(4), dimension(:), allocatable, save h
! Parameters used for hadronic cascade shower is generated newline ! For you may give as as or em as(order/case/separator insensitive) is to generate EM-cascade and AS. \newline ! Generate
! Zobs h header file for observation sites definition ! integer horizontal
max ptcl codes in the system(except krare) integer kphoton
dE dx *! Nuc Int sampling table d
! type magfield sequence ! Note that position vector where the magnetic field is given ! is not included here ! unit of field strength is in T(1 gauss=10 **-4 T) union map real *8 x
! type magfield sequence ! Note that position vector where the magnetic field is given ! is not included here ! unit of field strength is in ! in earth_center coordinate * y
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos Zfirst pos *Zfirst pos *Zfirst Zfirst vec w *Zfirst vec w Zfirst vec *Zfirst Zfirst where
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below this
block data cblkIncident data *Za1ry is
! structure defining a particle at production ! Basic idea of what is to be contained in ! the particle structue is that dynamical ones should be included those derivable from the particle code ! is not included ******************************************************type fmom momentum sequence union map real e endmap map real * x
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 Note