1 integer nsites ! # of sites
to be histogramed.
10 ! These must be
double. because addition
is done 1 by 1
11 ! and reach 10^7 or more.(at least
Ng,
Ne, at E0
14 * Ng(maxNoOfSites), Ne(maxNoOfSites),
15 *
Nmu(maxNoOfSites), Nhad(maxNoOfSites)
22 character*192 basefilename, basefilename2, filename
27 * fnoT=42, fnoL=43, fnoN=44, fnoB=45)
29 * webmin(nrbin, nfai, MaxNoOfSites)
30 real dErfai(nrbin, nfai, MaxNoOfSites)
31 real dECent( MaxNoOfSites )
36 common /Zprivatec/ Ng, Ne,
Nmu, Nhad, SumEloss,
40 common/ Zprivatec2/ basefilename
*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
! Zobs h header file for observation sites definition ! integer maxNoOfSites
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)
integer maxbin nregion c minsave drx2 ! drawing region real maxdep integer maxpos integer to
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
! 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