6 ! individual particle info.
is output.
13 ! These must be
double. because addition
is done 1 by 1
14 ! and reach 10^7 or more.(at least
Ng,
Ne, at E0
22 * nptcls(
nrbin, 4, MaxNoOfASsites),
23 * nrfaiRec(
nrbin,
nfai, 4, MaxNoOfASsites),
24 * nrfaiAll(
nrbin,
nfai, 4, MaxNoOfASsites),
31 common /Zprivatec/ SumEloss, CosRot, SinRot,
32 * rbin, recprob, nptcls, nrfaiRec, nrfaiAll,
34 *
histdep, indivdep, ansites, hnsites,
*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
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
integer maxn parameter(maxn=50000) real *8 wzmin
real ** Ng(maxNoOfASSites)
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
! 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