4 !
sys=
"xyz": origin
is center of the Earth.
5 !
x: directed
to longitude 0, latitude 0
6 !
y: directed
to longitude 90 deg, latitude 0
7 !
z: center
to the North pole
8 ! ****************************************************************
10 ! ****************************************************************
29 *
long, ! longitude in deg. +
is to the east.
34 *
phi, ! azimuthal angle
35 * radius ! radial distance
39 character*4
sys ! which
system.
'xyz',
'llh',
'sph' real(8), parameter, public m
latitude latitude this system is used *****************************************************************! type coord sequence union map real * r(3) endmap map real *8 x
latitude latitude this system is used *****************************************************************! type coord sequence union map real z z in m endmap xyz map real * lat
latitude latitude this system is used *****************************************************************! type coord sequence union map real z z in m endmap xyz map real ! latitude in deg is to the north * long
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz sys
real(4), dimension(:), allocatable, save h
integer maxbin nregion c minsave drx2 ! drawing region real maxdep integer maxpos integer to
max ptcl codes in the system(except krare) integer kphoton
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos Zfirst pos height
latitude latitude this system is used *****************************************************************! type coord sequence union map real y
latitude latitude this system is used *****************************************************************! type coord sequence union map real z z in m endmap xyz map real ! latitude in deg is to the north ! longitude in deg is to the east *h ! height in m endmap llh map real * theta
! common parameters for tracking(-> ------------------------------------- real *8 Deltpp !2 p-p xsection increases as $E^{Deltpp}$(E$>$ 100GeV) real *8 Deltpip !2 pi-p xsection increases as $E^{Deltpip}$(E$>$ 100GeV) real *8 Deltkp !2 k-p xsection increases as $E^{Deltkp}$(E$>$ 100GeV) real *8 IncreaseXsec !2 how the xsection increases. 1.0$\rightarrow$ power of E ! above ones are obsolete !<-) ----------------------------------- common/Zxsectionp/*Deltpp
latitude latitude this system is used *****************************************************************! type coord sequence union map real z z in m endmap xyz map real ! latitude in deg is to the north ! longitude in deg is to the east *h ! height in m endmap llh map real ! polar angle * phi
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