1 #define MAX_NO_OF_SITES 50 2 #define MAX_NO_OF_AS_SITES 50 3 #define MAX_SEGMENTS 60 4 #define MAX_NO_OF_COMPS 11 7 #define MAX_STACK 50000 9 #define MAX_USERHOOKC 5 10 #define MAX_USERHOOKI 10 11 #define MAX_USERHOOKR 10 12 ! make next
as large
as 1~2 milion
for practical
case 13 ! (
for general MPI application)
17 ! The
user may change the definition of the max
18 ! size of arrays by the above variables. It
is better not
24 !
MAX_SEGMENTS : Max number of segments when expressing primary
25 ! energy spectrum by
a segments of lines
for one
28 ! definable at
a time.
29 !
MAX_AZM 37: Max number of azimuthal angle bins
for a one rigidity
30 !
cut table can contain.
31 !
MAX_ZEN 21: Max number of zenith angle bins
for a one rigidity
32 !
cut table can contain.
34 ! no multi-particle produciton,
this may be
as small
35 !
as ~30, but
if multipicity of
n must be treated,
36 !
this must be
as large
as a few times *
n.
37 !
MAX_PTCL 8000: Max number of partilces at
a particle production.
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
#define MAX_NO_OF_AS_SITES
! 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
integer Charge2heavyG charge of heavy rightarrow heavy group index conversion array integer HeavyG2massN heavy group index rightarrow mass number conversion array integer HeavyG2charge heavy group index rightarrow charge of heavy conversion array integer HeavyG2code heavy group index rightarrow particle code conversion array integer Code2massN particle code rightarrow mass number conversion array integer Code2heavyG particle code rightarrow heavy group index conversion array real *FragmentTbl the number of interacting nucleons among a projectile heavy nucleus is ! determined as the number of first collision of each interacting nucleon inside ! the nucleus If
integer maxbin nregion c minsave drx2 ! drawing region real maxdep integer maxpos integer to
real cut integer ncpi0 common Zuserc cut
! 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
block data cblkIncident data *Za1ry is
1.8.13