6 ! individual particle info.
is output.
12 ! These must be
double. because addition
is done 1 by 1
13 ! and reach 10^7 or more.(at least
Ng,
Ne, at E0
16 * Ng(maxNoOfASSites), Ne(maxNoOfASSites),
Nmu(maxNoOfASSites),
17 * Nhad(maxNoOfASSites)
21 * nptcls(nrbin, 4, MaxNoOfASsites),
22 * nrfaiRec(nrbin,
nfai, 4, MaxNoOfASsites),
23 * nrfaiAll(nrbin, nfai, 4, MaxNoOfASsites),
24 * dErfai(nrbin, nfai, MaxNoOfASsites)
25 logical tklat ! take lateral or not
26 logical tkelosslat ! take energy loss lateral or not
51 common /Zprivatec/ Ng, Ne,
Nmu, Nhad, SumEloss,
52 * rbin, recprob, nptcls, nrfaiRec, nrfaiAll,
54 *
histdep, indivdep, ansites, hnsites,
55 * recxy, tkelosslat, tklat,
tkarspec,
56 * tkrespec, tkrzspec, tkzfspec, tkrfspec, tkefspec,
58 * tkretspec, tkrezspec, tkrzfspec, tkrefspec
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos depth
integer lengid integer lengdir character *dir integer kgetenv2 character *numb character *execid character *msg logical takehist save do nsites histdep(i)=0indivdep(i)=0enddo leng
real Nmu(maxNoOfASSites) common/Zprivatec0/Ng
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz r
logical KeepWeight ! see setupenv sh logical tkarspec ! get lateral histo in a web fai bin logical tkrtspec ! get time histo in each web bin logical SeeLowdE common Zprivatec2 tkrtspec
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
integer maxn parameter(maxn=50000) real *8 wzmin
logical KeepWeight ! see setupenv sh logical tkarspec ! get lateral histo in a web fai bin logical tkrtspec ! get time histo in each web bin logical SeeLowdE common Zprivatec2 tkarspec
real ** Ng(maxNoOfASSites)
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
! 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
! Zobs h header file for observation sites definition ! integer maxNoOfASSites
! Parameters used for hadronic cascade shower is generated newline ! For you may give as as or em quick generation of AS for heavy primaries is tried See chookASbyH f character *Generate2 don t touch this for skeleton flesh use integer MagBrem no magnetic bremsstrahlung is considered newline ! if and Ee energy loss due to magnetic brems is considered newline ! if and Ee real sampling of gamma is performed WaitRatio ! must be made small so that WaitRatio *E0 sim MagBremEmin integer MagPair no magnetic pair creation is considered newline ! if and Eg real sampling is tried(note, actually upsilon is referred further). To see these magnetic effects
block data cblkIncident data *Za1ry is