FORSCHUNGSZENTRUM KARLSRUHE
Technik und Umwelt
Wissenschaftliche Berichte
FZKA 6097B
Upgrade of the Monte Carlo Co de
CORSIKA to Simulate Extensive
20
Air Showers with Energies > 10 eV
1y
D. Heck and J. Knapp
Institut fur Kernphysik
Contribution to: Fruhjahrstagung DPG, Sektion Teilchenphysik, March 23 - 27 1998 Freiburg
1
Universitat Karlsruhe, Institut fur Exp erimentelle Kernphysik, Postfach 3640, D-76021 Karlsruhe
y
now at: Department of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, United
Kingdom
Forschungszentrum Karlsruhe GmbH, Karlsruhe 1998
Abstract
Upgrade of the Monte Carlo Co de CORSIKA to Simulate Extensive Air
20
Showers with Energies > 10 eV.
To simulate extensive air showers at the highest observed energies the Monte Carlo
program CORSIKA has b een upgraded. Some of the uncertainties of the hadronic
interaction mo dels are discussed. The e ects of the Landau-Pomeranchuk-Migdal
e ect onto the electromagnetic interactions are demonstrated. Sp ecial attention is
paid to the reduction of computing time by the so-called thin sampling pro cedure.
Various steps of realization are describ ed. The in uence of the thin sampling level
on the computing time, the particle numb ers to b e stored, and the statistical uncer-
tainty of observable quantities is shown.
Zusammenfassung
Erweiterung des Monte Carlo Programms CORSIKA zur Simulation von
20
ausgedehnten Luftschauern bei den hochsten Energien > 10 eV.
Um ausgedehnte Luftschauer b ei den hochsten b eobachteten Energien simulieren zu
konnen, wurde das Monte Carlo Programm CORSIKA erweitert. Einige der Un-
sicherheiten, die von den hadronischen Wechselwirkungsmo dellen herruhren, werden
diskutiert. Der Ein u des Landau-Pomeranchuk-Migdal-E ekts auf die elektromag-
netischen Wechselwirkungen wird gezeigt. Besondere Aufmerksamkeit wird der Ver-
minderung der Rechenzeit durch das sogenannte thin sampling Ausdunnen gewid-
met. Die verschiedenen Ausbaustufen der Verwirklichung werden b eschrieb en. Der
Ein u der Ausdunn-Energiesc hwelle auf die Rechenzeit, auf die Zahl der zu sp ei-
chernden Teilchen und auf die statistische Unsicherheit b eobachtbarer Groen wird
demonstriert. i ii
Contents
1 Intro duction 1
2 Features of Hadronic Interaction Mo dels 1
3 Landau-Pomeranchuk-Migdal E ect 2
4 Computing Times without Thin Sampling 6
5 Stages of Thin Sampling 7
6 In uence of Thin Sampling Level 14
7 Conclusions 16
Acknowledgement 17
Bibliography 18 iii iv
1 Intro duction
15
The greatest mysteries of cosmic rays arise at the highest energies > 10 eV where
the ux b ecomes so small that very extended detectors on the Earth surface are
needed to collect statistically sucient events. Such exp eriments cannot observe
the primary cosmic ray directly, rather they lo ok for the reaction pro ducts which
are formed within the atmosphere - so-called extensive air showers EAS - and
reach the ground. The reconstruction of parameters like primary energy, direction,
and mass is not straight forward and is only p ossible by comparison with detailed
simulations of the EAS development. Sp ecial interest is given to cosmic rays at
19
energies > 10 eV as their existence cannot to be explained consistently with
known acceleration mechanisms and large distance traveling through space without
o
interaction with the 3 K background radiation. The exp erimental determination of
the precise form of the sp ectrum at those energies would b e a great step towards the
understanding of the highest energy particles in the universe.
Some exp eriments to explore this energy range are now op erating or are to be real-
ized in the near future. Such exp eriments are AGASA [1], HiRes [2], Yakutsk [3],
and AUGER [4]. Understanding the exp erimental data and designing new exp eri-
ments needs Monte Carlo programs to simulate EAS at the highest energies. These
requests pushed the upgrade of the EAS simulation program CORSIKA COsmic
20
Ray SImulation for KASCADE for primary energies > 10 eV .
2 Features of Hadronic Interaction Mo dels
The most essential ingredient to simulate EAS at the highest energies is a hadronic
interaction mo del which is able to extrap olate from lower energies covered by man
made accelerators to predict the results of hadronic collisions at higher energies.
This extrap olation is sub ject to considerable uncertainties. Within CORSIKA [5] we
have ve high energy hadronic interaction mo dels available: DPMJET [6], HDPM
[5, 7], QGSJET [8], SIBYLL [9], and VENUS [10]. A comparison of these mo dels
is given in Ref. [11, 12, 13] which includes a discussion of their in uences onto the
development of EAS. Three of these mo dels DPMJET, QGSJET, and SIBYLL
handle hadronic interactions at the highest energies. Besides other features they
predict also the hadronic interaction cross-section which determines the free path
between twointeractions. The predicted inelastic proton-air cross-sections are shown
in Fig. 1. Despite the extrap olation, the cross-section values of the three mo dels do
not di er by more than 15 at the highest energies. Similarly the average
transverse momentum which is carried away by the charged secondary particles is
19 5
predicted in accordance up to p =10 eV corresp onding with E 10 GeV ,
lab cm
see Fig. 8 of Ref. [11] and the di erences between the considered mo dels are
mo derate. Other extrap olations like the average multiplicity of secondary charged
particles emerging from a hadronic collision are a ected byamuch larger uncertainty 1 650 (mb)
p-air inel 600 Mielke et al. σ Yodh et al. 550 Gaisser et al. Baltrusaitis et al. 500 Honda et al.
450
400 VENUS 350 QGSJET SIBYLL 300 HDPM 250 DPMJET
2 3 4 5 6 7 8 9 10 11 10 10 10 10 10 10 10 10 10 10
plab (GeV/c)
Figure 1: Inelastic p-air cross-sections for hadronic interaction mo dels available in
CORSIKA together with exp erimental values derived from air shower exp eriments
5
[14, 15, 16, 17, 18]. The shaded band represents a t to the data b elow p < 10 GeV
lab
p air
2
according to the function = a log p+b log p+c .
inel
of more than a factor 2.5, as maybe deduced from Fig. 5 of Ref. [11].
3 Landau-Pomeranchuk-Migdal E ect
When an ultrarelativistic electron emits a low-energy photon via bremsstrahlung, the
longitudinal momentum transfer b etween electron and the involved target nucleus can
be very small. Because of the uncertainty principle this implies a long distance within
which the momentum transfer takes place, known as the formation length. If anything
happ ens to one of the two electromagnetic particles while traveling this formation
length, the emission can b e disrupted. This is known as Landau-Pomeranchuk-Migdal
LPM e ect [19, 20] and is caused predominantly by electron multiple scattering.
This e ect b ecomes noticeable for interactions in air at normal temp erature and
16
pressure at energies ab ove 10 eV . Similarly in pair pro duction it suppresses
+
predominantly those events in which the energy is shared symmetrically by the e e -
pair.
The magnitude of this suppression dep ends on the pro duct E of the energy
em air
of inducing em-particle and of the air density at the interaction p oint. Consequently
for pair pro duction and bremsstrahlung the LPM-e ect is most pronounced at high
energies and low altitudes. In Figs. 2 and 3 the in uences of the LPM-e ect are shown 2 1 v•dN/dv
-1 10
E = 1016eV e 17 Ee= 10 eV E = 3.1*1017eV e 18 Ee= 10 eV -2 E = 3.1*1018eV 10 LPM bremsstrahlung e 19 Ee= 10 eV altitude = 20 km E = 3.1*1019eV e 20 Ee= 10 eV
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
fractional photon energy v = Eγ/E0
1 v•dN/dv
-1 10
E = 1016eV e 17 Ee= 10 eV E = 3.1*1017eV e 18 Ee= 10 eV -2 E = 3.1*1018eV 10 LPM bremsstrahlung e 19 Ee= 10 eV altitude = 10 km E = 3.1*1019eV e 20 Ee= 10 eV
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
fractional photon energy v = Eγ/E0
1 v•dN/dv
-1 10
E = 1016eV e 17 Ee= 10 eV E = 3.1*1017eV e 18 Ee= 10 eV -2 E = 3.1*1018eV 10 LPM bremsstrahlung e 19 Ee= 10 eV altitude = 0 km E = 3.1*1019eV e 20 Ee= 10 eV
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
fractional photon energy v = Eγ/E0
Figure 2: LPM-e ect for bremsstrahlung at various altitudes and electron energies.
On the ordinate the pro duct of fractional photon energy v times the rate p er fractional
photon energy dN =dv is given in arbitrary units. For the lowest energies it coincides
with the Bethe-Heitler case.
16 20
for 3 altitudes and various energies b etween E =10 eV and 10 eV . These curves
have b een established with a routine taken from the FORTRAN program AIRES 3 1 dN/du
-1 10 E = 1017eV γ 18 Eγ= 10 eV LPM pair formation E = 3.1*1018eV γ 19 Eγ= 10 eV E = 3.1*1019eV altitude = 20 km γ 20 Eγ= 10 eV
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
fractional electron energy u=Ee/E0
1 dN/du
-1 10 E = 1017eV γ 18 Eγ= 10 eV LPM pair formation E = 3.1*1018eV γ 19 Eγ= 10 eV E = 3.1*1019eV altitude = 10 km γ 20 Eγ= 10 eV
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
fractional electron energy u=Ee/E0
1 dN/du
-1 10 E = 1017eV γ 18 Eγ= 10 eV LPM pair formation E = 3.1*1018eV γ 19 Eγ= 10 eV E = 3.1*1019eV altitude = 0 km γ 20 Eγ= 10 eV
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
fractional electron energy u=Ee/E0
Figure 3: LPM-e ect for pair pro duction at various altitudes and gamma-ray en-
ergies. On the ordinate the rate per fractional electron energy dN =du is given in
arbitrary units. For the lowest energies it coincides with the Bethe-Heitler case.
[21] which corresp onds with the PASCAL program MOCCA [22] and adapted to
CORSIKA.
According to Migdal [20] we use the fractional form to express the energies in the 4
cross-section formulas following Ref. [23]. The probabilityper radiation length that
an electron with energy E radiates a photon with energy k between v = k=E and
v + dv is given by
n h i o
2 2
W E; vdv = s v Gs+2 1+1 v s dv =3v : 1
b
Similarly the probability p er radiation length that a photon with energy K pro duces
+
an e e pair with electron energy p between u = p=K and u + du is expressed as
n h i o
2 2
s Gs+2 u +1 u s du=3 : 2 W K; udu =
p
The fractional energy of the accompanying p ositron is 1 u. The functions Gs,
s, and s are given by
"
1
2
X
1
2 3 2
Gs= 12s 48s 1= s + l + + s 3
2
l =0
1
i h
X
2 3 2 2
4 s=6s 6s +24s 1= s + l + s
l =0
8
>
2 : s s
1
>
>
<
1 + ln s= ln s : s
s= 5
1 1
>
>
>
:
1 : s 1
2
1=3
where s = Z =190 and Z is the atomic numb er of the medium. The parameters
1
s in Eqs. 1 and 2 are given by s and
"
v 1
2 2
t = mc =E 6 s =
0 e
16 1 v s
"
1 1
2 2
t = mc =K 7 s =
0 e
16 u 1 u s
where is the ne structure constant, t is the density dep endent radiation length,
0
2
is the Compton wavelength of the electron, and mc is the electron rest mass.
e
1
Both t and are measured in cm. The value of the factor t = amounts to
0 e 0 e
16
3 2
1:371 10 t . For s>1 rsp. s>1 which corresp onds with a small energy of the
0
primary particle, Eqs. 1 and 2 transform into the familiar Bethe-Heitler cross-sections
BH BH
W and W as for this case G =1,=1,and =1.
b p
In CORSIKA the simulation p erforms a bremsstrahlung or pair pro duction event
always in the ordinary Bethe-Heitler case. Before the secondary particles are stored,
the LPM-routine is applied to check, whether the interaction should be discarded.
2 2 2
First s rsp. s is calculated according to Eq. 6 rsp. 7, assuming = 1. If s 1
2
rsp. s 1 then we are in the ordinary low energy case and the parameters of the 5
secondary particles are completed and stored, the LPM-e ect is not e ective.
2 2
If s < 1 rsp. s < 1 then is determined according Eq. 5 and s rsp. s are corrected
for the contribution of . In the next step Gs Eq. 3 and s Eq. 4 are
calculated using the approximations
0 0
Gs= g =1 + g
with
8
2
<
s [14:1+2:36=s +0:1] : s<0:1
0
g =
2
:
s [24:0+0:0394=s 0:08] : s 0:1
and
8
2
<
6s 16s : s<0:1
s=
2
:
6s +24s [=4 arctan 0:944 + 0:59=s] : s 0:1 :
Finally, the fraction of the LPM cross-section according to Eq. 1 or 2 relative to
LP M BH LP M BH
the Bethe-Heitler cross-section W =W rsp. W =W is compared with
b b p p
a random numb er, distributed uniformly between 0 and 1. If the random number
exceeds the fraction, the interaction is rejected and the primary particle is propagated
further on without interaction.
By the LPM-suppression electromagnetic interactions with extreme asymmetric en-
ergy distribution of the secondary particles are favoured. Together with the de-
creasing electromagnetic cross-sections the net e ect onto EAS results in a retarded
shower development accompanied by enhanced uctuations, as demonstrated in Ref.
[23]. Our tests with and without LPM-e ect revealed, that for -induced showers of
20 2
10 eV the shower maximum is pushed deep er into the atmosphere by 100 g =cm
slant depth when switching on the LPM-e ect.
4 Computing Times without Thin Sampling
The most serious problem in simulations of EAS at the highest energies comes from
the nearly linear growth of CPU computing time with primary energy. This growth
is shown in Fig. 4 for CORSIKA using EGS full Monte Carlo simulation for electro-
magnetic particles. Even with to day's fast work stations, CPU-times of more than
20
one year per shower for E > 10 eV exceed the limits of feasibility.
0
A similar problem comes from the huge number of secondaries to be stored when
arriving at a detector. The number of 1 million particles for a vertical shower of
15
10 eV primary energy will explo de to ab out 200 billion particles per shower at
20
10 eV , which need a storage of 6 TB. 6 1y
10 7 CPU time (sec)
1m
10 6
5 10 1d
10 4
1h
10 3 DEC ALPHAstation 500/500 MHz
10 2 14 15 16 17 18 19 20 10 10 10 10 10 10 10
primary energy (eV)
Figure 4: CPU-time as function of the proton primary energy at vertical incidence
for full simulation of the electromagnetic subshowers.
5 Stages of Thin Sampling
An elegantway out of the enormous CPU-time consumption has b een demonstrated
by M. Hillas [22, 24] using the so-called thin sampling algorithm. It resembles the
variance reduction metho d of Ref. [25].
If an EAS has develop ed far enough so that very many particles are available, out
of the particles emerging from an interaction only one is selected at random to be
representative and considered for the further development of the EAS, while the
other secondary particles are discarded. The advantage of this pro cedure is a drastic
reduction of CPU-time, as now only few particles have to be followed down to the
ground, while the huge bulk of low energy particles not essentially contributing to
the gross development of an EAS maybe dropp ed.
The decision at which p oint of the EAS development the thin sampling should start
is given by the thinning level energy E , usually expressed by the fraction " of
thin thin 7 a) b) c) d)
< > > > E0 Ethin E0 Ethin E0 Ethin E0 Ethin
p ~E /ΣE for ΣE
pi~Ei/Ethin for Ei Ethin
Figure 5: Various development stages of thin sampling in CORSIKA. The horizontal
bar b elow the vertex symb olizes the thin sampling algorithm which passes only the
selected particles and discards the others. The summation in case a and b runs
always over all i secondary particles with E < E . Particles with E > E are
i thin thin
represented by thick solid lines, particles selected by the thin sampling pro cedure are
given by dashed-dotted and dotted lines.
the primary energy E
0
" = E =E :
thin thin 0
The representative particle is selected at random resp ecting a survival probability
according to its energy fraction of the energy summed over all secondary particles
emerging from the interaction under consideration. This is drawn schematically in
Fig. 5a. After the selection pro cess the weight of this particle is increased by
the weight factor w = 1=p . This rst stage a of development has b een used in
i i
CORSIKA version 5.20.
Results showed, that still many low energy particles had to be followed. They orig-
inate from interactions induced by particles with energy E >E . Typically such
thin
low energy particles with E < E emerge in the target rapidity range together
i thin
with high energy particles with energies E >E in the pro jectile rapidity range.
i thin
Dep ending on the energy sum of those i secondaries with E < E the survival
i thin
probability p is calculated to
i
P P
p = E = E for E < E and
i i i i thin
i i
P
p = E =E for E E :
i i thin i thin
i 8 E E 10 10 10 5 10 5
4 4 + - 10 γ 10 e e d(N)/dlog d(N)/dlog 10 3 10 3
10 2 10 2
10 10
1 1
-1 -1 10 10
-2 -2 10 10 -2 0 2 4 -2024 log10(energy) (GeV) log10(energy) (GeV) E E 10 10 10 4 10 4 + - 3 µ µ 3 h d(N)/dlog 10 d(N)/dlog 10
10 2 10 2
10 10
1 1
-1 -1 10 10
-2 -2 10 10 024 024
log10(energy) (GeV) log10(energy) (GeV)
Figure 6: Unweighted energy sp ectra resulting from various stages of thin sampling
+
for gamma rays top left, e and e top right, muons b ottom left, and hadrons
b ottom right. The gradual reduction of the numb er of particles reaching the detec-
tor level is demonstrated for the cumulative application of the development stages a
solid line, b dashed line, c dashed-dotted line, and d thin solid line. Proton
6
induced EAS with E =10 GeV at vertical incidence have b een thin sampled b elow
0
3
E =10 GeV for this gure using CORSIKA default parameters. The dotted line
thin
on top of each plot is without thin sampling.
P
In the second case the summation over all p results in p > 1 indicating that
i i
eventually more than one secondary is retained for further treatment. This case is
represented by Fig. 5b. 9 E E 10 10 10 5 10 5
4 4 + - 10 γ 10 e e
d(N•w)/dlog 10 3 d(N•w)/dlog 10 3
10 2 10 2
10 10
1 1
-1 -1 10 10
-2 -2 10 10 -2 0 2 4 -2024 log10(energy) (GeV) log10(energy) (GeV) E E 10 10 10 4 10 4 µ+µ- 10 3 10 3 h d(N•w)/dlog d(N•w)/dlog 10 2 10 2
10 10
1 1
-1 -1 10 10
-2 -2 10 10 024 024
log10(energy) (GeV) log10(energy) (GeV)
Figure 7: Weighted energy sp ectra resulting from various stages of thin sampling.
The good coincidence of thin sampled sp ectra in all stages of development with
sp ectra without thin sampling is obvious. Parameters as in Fig. 6.
A sp ecial form of case b o ccurs mainly in the electromagnetic cascade, where only
two particles emerge from an interaction which is shown in Fig. 5c. Here the
incoming particle carries an energy E >E while b oth secondaries fall b elow the
thin
threshold E . Again one particle is kept, while the chance p for retaining also the
thin j
other secondary j is given by
P
X X
E
i
i
p = p 1= 1 > 0 as E >E :
j i i thin
E
thin
i i
An additional variant of case b is shown in Fig. 5d frequently o ccurring in brems- 10 w w 10 γ 10 + - 10 3 10 3 e e dN/dlog dN/dlog 10 2 10 2
10 10
1 1
-1 -1 10 10
-2 -2 10 10 02460246 log10(weight) log10(weight) w w 10 µ+µ- 10 10 3 10 3 h dN/dlog dN/dlog 10 2 10 2
10 10
1 1
-1 -1 10 10
-2 -2 10 10 02460246
log10(weight) log10(weight)
Figure 8: Distribution of particle weights resulting from various stages of thin sam-
pling. Parameters as in Fig. 6.
strahlung pro cesses. The electron or p ositron keeps the dominant p ortion of the
incoming energy E >E , while a low energy photon is radiated o . The chance p
thin
to retain the photon is given by its energy E relativeto the thinning level energy
ph
p = E =E :
ph thin
All these cases a to d have b een realized successively within CORSIKA, each con-
tributing to the improvement of the simulation sp eed. In the following the results
gained with the higher development stages always comprise the use of the earlier
development stages, i.e. stage d comprises stages a to c. 11 6 6 γ + - 5 5 e e log10(weight) log10(weight) 4 4
3 3
2 2
1 1
0 0 -2024 -2 0 2 4 log10(energy) (GeV) log10(energy) (GeV)
6 6 µ+µ- 5 5 h log10(weight) log10(weight) 4 4
3 3
2 2
1 1
0 0 024 024
log10(energy) (GeV) log10(energy) (GeV)
Figure 9: Two dimensional histogram of energy and weight of di erent particle sp ecies
after thin sampling with stage a. Parameters as in Fig. 6.
In Fig. 6 the energy sp ectra of particles arriving at the detector level 110 m a.s.l.
are plotted for di erent particle sp ecies separately, each arriving particle counted
once. The sp ectra without thin sampling are the dotted ones on top of each plot. The
reduction in all particle sp ecies pro ceeds with the cumulative activation of stages a to
d. The e ectiveness of stage d on the electromagnetic particles without signi cant
in uence onto muonic and hadronic sp ectra is obvious. If the energy sp ectra are
formed resp ecting the particle weights the sp ectra coincide well with sp ectra without
thin sampling, as is demonstrated in Fig. 7.
The distribution of particle weights is shown in Fig. 8. With increasing improvement 12 6 6 γ + - 5 5 e e log10(weight) log10(weight) 4 4
3 3
2 2
1 1
0 0 -2024 -2 0 2 4 log10(energy) (GeV) log10(energy) (GeV)
6 6 µ+µ- 5 5 h log10(weight) log10(weight) 4 4
3 3
2 2
1 1
0 0 024 024
log10(energy) (GeV) log10(energy) (GeV)
Figure 10: Two dimensional histogram of energy and weight of di erent particle
sp ecies after thin sampling including stages a to d. Parameters as in Fig. 6.
of thin sampling the numb er of particles with lowweight is reduced. The o ccurrence
of weight w = 1, as is carried by particles which did not undergo thin sampling, is
caused by high energy particles with E > E p enetrating down to the detector
thin
level. This is observed esp ecially for muons and hadrons.
The scatter plots Figs. 9 and 10 showhowby the stages of improvement the number
of low energy particles with lowweight lower left corner of each plot is reduced and
only the particles along the line w E = const survive the optimized thin sampling
pro cedure. In Fig. 10 still some few electromagnetic and muonic particles in the
low energy range fall b elow this line. These particles originate from weak interaction 13
3
without thin thin sampled b elow 10 GeV
sampling a b c d
1
CPU-time 8 h 5'44" 1'20" 31" 18"
particles 550960 8136 2187 357 323
1
for DEC 3000/600 AXP 175 MHz
Table 1: E ects of various stages a to d of thin sampling. The numb ers are deter-
6
mined for proton primaries of 10 GeV at vertical incidence, QGSJET interaction
mo del, and CORSIKA default parameters.
decays with emission of an unobservable neutrino. In these decays the energy sum of
outgoing observable particles di ers from the incoming energy, which results in a
lackof energy for a xed weight. Therefore this e ect do es not happ en for hadrons.
The deviation from the straight line in the low energy tail of the hadronic particles
visible in Fig. 10 results from the rest mass included into the energy used in the
formulas to calculate the weight of surviving particles.
The gain in computing time and reduced storage requirements for the di erent stages
a to d of thin sampling development is summarized in Table 1. Besides the di erent
stages also the thin sampling level in uences the required storage for particles and
the required CPU-time as is discussed b elow.
6 In uence of Thin Sampling Level
Of great imp ortance is the energy level b elow which thin sampling b ecomes active.
Generally with decreasing thin sampling level the CPU-time is increasing and also the
numb er of particles reaching the detector. This is demonstrated with proton induced
19
EAS of 10 eV primary energy in Fig. 11 for the two development stages a and
d. But the gain in CPU-time and the reduced storage requirement with increasing
thin sampling level is paied for by an increasing uncertaintyininteresting observable
quantities. This uncertainty is caused by the additional statistical uctuations which
stem from the thin sampling pro cess. The dep endence of the uncertainty on the thin
sampling level is shown in Fig. 12 for the number of muons arriving at the detector
level.
As an example we compare in Fig. 12 the thin sampling levels of di erent develop-
ment stages which lead to the same uncertaintyofmuon numb er: To reach=