The study in the primary energy range eV  eV with the Muon Tracking Detector in the KASCADEGrande experiment
Abstract
The KASCADEGrande Muon Tracking Detector enables with high accuracy the measurement of directions of EAS muons with energy above 0.8 GeV and up to 700 m distance from the shower centre. Reconstructed muon tracks are used to investigate muon pseudorapidity () distributions. These distributions are nearly identical to the pseudorapidity distributions of their parent mesons produced in hadronic interactions. Comparison of the distributions from measured and simulated showers can be used to test the quality of the high energy hadronic interaction models. In this context a comparison of the QGSJetII2 and QGSJetII4 model will be shown. The pseudorapidity distributions reflect the longitudinal development of EAS and, as such, are sensitive to the mass of the cosmic rays primary particles. With various parameters of the distribution, obtained from the MTD data, it is possible to calculate the mean logarithmic mass of CRs. The results of the analysis in the primary energy range eV  eV with the quartile (Q1) of distribution will be presented.
from MTD in KASCADEGrande
\authors
P. Łuczak,
W.D. Apel,
J.C. ArteagaVelázquez,
K. Bekk,
M. Bertaina,
J. Blümer,
H. Bozdog,
I.M. Brancus,
E. Cantoni,
A. Chiavassa,
F. Cossavella,
C. Curcio,
K. Daumiller,
V. de Souza,
F. Di Pierro,
P. Doll,
R. Engel,
J. Engler,
B. Fuchs,
D. Fuhrmann,
H.J. Gils,
R. Glasstetter,
C. Grupen,
A. Haungs,
D. Heck,
J.R. Hörandel,
D. Huber,
T. Huege,
K.H. Kampert,
D. Kang,
H.O. Klages,
K. Link,
M. Ludwig,
H.J. Mathes,
H.J. Mayer,
M. Melissas,
J. Milke,
B. Mitrica,
C. Morello,
J. Oehlschläger,
S. Ostapchenko,
N. Palmieri,
M. Petcu,
T. Pierog,
H. Rebel,
M. Roth,
H. Schieler,
S. Schoo,
F.G. Schröder,
O. Sima,
G. Toma,
G.C. Trinchero,
H. Ulrich,
A. Weindl,
J. Wochele,
J. Zabierowski
KASCADEGrande Collaboration
\afiliations
National Centre for Nuclear Research, Department of Cosmic Ray Physics, Łódź, Poland
Institut für Kernphysik, KIT  Karlsruher Institut für Technologie, Germany
Universidad Michoacana, Instituto de Física y Matemáticas, Morelia, Mexico
Dipartimento di Fisica, Università degli Studi di Torino, Italy
Institut für Experimentelle Kernphysik, KIT  Karlsruher Institut für Technologie, Germany
Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest, Romania
Osservatorio Astrofisico di Torino, INAF Torino, Italy
Universidade So Paulo, Instituto de Física de São Carlos, Brasil
Fachbereich Physik, Universität Wuppertal, Germany
Department of Physics, Siegen University, Germany
Dept. of Astrophysics, Radboud University Nijmegen, The Netherlands
Department of Physics, University of Bucharest, Bucharest, Romania
now at: Istituto Nazionale di Ricerca Metrologia, INRIM, Torino;
now at: MaxPlanckInstitut Physik, München, Germany;
now at: University of DuisburgEssen, Duisburg, Germany;
now at: University of Trondheim, Norway.
\emailluczak@zpk.u.lodz.pl
\keywordsKASCADEGrande, extensive air shower, Muon Tracking Detector, muon pseudorapidity, mass composition, model tests
1 Introduction
The Muon Tracking Detector (MTD) [1] is one of the detector components in the KASCADEGrande EAS experiment [2] (see figure 1), operated at the Karlsruhe Institute of Technology (KIT)  Campus North, in Germany, by an international collaboration.
The MTD is a group of 16 muon telescopes made out of streamer tube gas detectors grouped in four modules (3 horizontal, 1 vertical). Each module is 2 meters wide and 4 meters long, the total detection area for vertical particles is 128 m. The telescopes are located in a concrete tunnel, under the ground level, covered by the concretesoiliron absorbent which absorbs a large fraction of accompanying lowenergy particles, thus enhancing the identification of muons with an energy exceeding 0.8 GeV. The MTD measures directions of muon tracks in extensive air showers with excellent angular resolution of .
These directional data enables to investigate the longitudinal development of the muonic component in air showers which is a signature of the development of the hadronic EAS core, being in turn dependent on the mass of the primary cosmic ray particle initiating a shower. Such studies can be done either by the determination of a mean muon production height [3] or by using the mean pseudorapidity () of EAS muons, expressed in terms of their tangential () and radial () angles, quantities reconstructed from muon tracks obtained with the MTD [4], [5].
In this work we compare the distributions for two high energy hadronic interaction models and investigate to what extent one can use the value of the first quartile of the muon pseudorapidity distribution for the determination of the mean logarithm of mass of cosmic rays.
2 Comparison of the distributions for QGSJetII2 and QGSJetII4
In this analysis EAS initiated by the proton and iron primary CR particles were simulated with the CORSIKA[6] code. The showers were simulated with two combinations of hadronic interaction models QGSJetII2
Figures 2 and 3 show the comparison between distribution for showers simulated with the QGSJetII2+FLUKA (solid line) and the QGSJetII4+FLUKA (dashed line) for proton and iron initiated showers. Both distributions are obtained for 0.5 and 0 in the muontoshoweraxis distance range 250370 meters, for showers with zenith angle up to 18 and in the energy range . Due to the low statistics, there are no statistically significant differences between the distributions for both models which was checked by performing a KolmogorovSmirnov test for the two distributions. The shapes of the distributions are similar for both primary particles, but one can notice that the mean of the proton distribution is lower for the QGSJetII4 model. This shift in proton distribution should be confirmed with more statistics of showers. If it remains, it will improve the separation between proton, data and iron distributions, leading to more reliable description of the measured data and higher values of the calculated from the parameters of the distributions.
3 The mass sensitivity of the EAS muon pseudorapidity
In the KASCADEGrande experiment muons can be registered up to 700 meters from the shower core. In case of the presented analysis, the muontoshoweraxis (R) distance ranges (Table 1) were limited to the distances where the mass composition of the detected showers is constant (taking into account the / ratio) in each energy range, not being affected by experimental inefficiencies. In this analysis EAS initiated by the proton and iron primary CR particles were simulated with the CORSIKA code (QGSJetII2 as high energy hadronic (HE) interaction model and FLUKA2002.4 as low energy (LE) hadronic interaction model). The QGSJet model was used to simulate hadronic interactions of particles with energy above 200 GeV while the FLUKA model was used to simulate ones below this energy. In the measured and simulated data only showers with zenith angle up to 18 are analysed.
R [m]  

10 GeV  
7.07.3  1.340.01  250370 
7.0  2.070.01  250370 
7.37.6  2.690.01  280400 
7.67.9  5.340.02  280430 
7.6  7.150.07  280430 
It was shown in [10] that the pseudorapidity of EAS muons is a parameter sensitive to the mass of the primary cosmic ray particles and can be used to calculate the of CRs.
It was observed that the distributions from the measured data are bracketed by the distributions from simulated showers and that and are linearly dependent. The main conclusion of the analysis is that is increasing with the energy, but its values are lower than expected in the energy range eV eV. The reason for this are the distortions of the distributions in simulations caused by the large number of high muons that are shifting the of the distributions towards higher values. As a result of this shift, the distributions from the measurement are close to these from proton initiated showers. This behaviour is caused by the QGSJet model which creates too many high muons that reach the observation level.
This conclusion was confirmed by the investigation of the muon production height (h)[3]. In this analysis, h of muons, from proton and iron initiated showers, were compared with muons from detected showers divided into these from light and heavy primary CR particles. It occurred that the measured values of the h above 4 km, where the interactions are described by the HE interaction model are not well described by the simulations. This is not the case for h4 km, where interactions are described by the LE interaction model, and where the measurements are compatible with the simulations.
In the simulations it is possible to distinguish between muons that stem from hadrons originating from HE or LE hadronic interaction models. Dividing the muon sample into those originating from grandparent hadrons with energy above or below 200 GeV (HE and LE muon sample) we can analyse contributions of these two groups to the combined distribution which can be obtained with the MTD data. In the R distance range that is valid in case of the MTD analysis (250400 meters), there are about 70% muons from the LE sample and 30% from the HE sample. Most of the muons from the latter contribute to the of the distribution, creating a long tail of high muons. In this pseudorapidity range the number of muons from this HE sample is larger than from the LE sample (with ratio about 60% to 40%, respectively).
The differences in the LE and HE muon sample contributions suggest that calculating the with the quartile of the distribution, where LE interaction model dominates, will remove the influence of the HE interaction model. However, the value of the quartile is affected by the tail from the HE muon sample and the analysis is biased. To eliminate this bias the radial and tangential angle cuts were applied. The effect of such angular cuts is depicted in figure 4. Here the distributions from the QGSJet and FLUKA muon sample (solid line and dotted line, respectively) are compared before and after the angular cuts (thin and bold solid lines and dotted lines, respectively). However, in the experiment conditions, it is necessary to take into account the statistics of available simulated and measured data. That is why it was not possible to eliminate completely the tails from the HE muon sample.
The angular cuts which provide as small as possible decrease in muon statistics and suppress tails from HE sample without significant distortion of the shape of the pseudorapidity distribution are: 0.75 and 0.2.
4 Results and conclusions
Comparison between distributions from H and Fe showers simulated with QGSJetII2 and QGSJetII4 model revealed that with the currently available statistics it is too early to draw conclusions about the quality of the new version of the QGSJet model. However, the shift in the of the proton distribution indicates that with larger number of simulated showers we can expect that a change in separation between proton and iron distributions, will lead to the better description of the measured data by the simulations.
Combining angular cuts and the quartile analysis, the calculated values of the rise with the energy at a similar rate as the ones from electron/muon detectors in KASCADE, and are significantly higher than in the previous analysis [10]. These values are similar to the results obtained by KASCADE with the e/m ratio analysis using QGSJet model [11].
The results of the calculations are presented in Table 2, and compared with the results from the previous analysis where the values of the distributions were used. Figure 5 shows the results from the table compared with the results from electron/muon detectors in the KASCADE experiment.
, QGSJetII2  

10 GeV  
7.0  7.3  1.340.01  1.980.55  0.710.24 
7.0  2.070.01  2.680.36  1.150.14 
7.3  7.6  2.690.01  2.790.87  1.870.34 
7.6  7.9  5.340.02  3.171.12  2.130.48 
7.6  7.150.07  3.860.81  2.450.39 
Acknowledgment:The authors would like to thank the members of the engineering and technical staff of the KASCADEGrande collaboration, who contributed to the success of the experiment. The KASCADEGrande experiment is supported in Germany by the BMBF and by the ÂHelmholtz Alliance for Astroparticle Physics  HAPÂ funded by the Initiative and Networking Fund of the Helmholtz Association, by the MIUR and INAF of Italy, the Polish Ministry of Science and Higher Education, and the Romanian Authority for Scientific Research UEFISCDI (grants PNIIIDEI 271/2011 and RUPD 17/2011).
Footnotes
References
 P. Doll et al., Nucl. Instrum. Methods Phys. Res., Sect. A 488 (2002) 517.
 W.D. Apel et al. (KASCADEGrande Collaboration), Nucl. Instrum. Methods Phys. Res., Sect. A 620 (2010) 202.
 W.D. Apel et al., (KASCADEGrande Collaboration), Astroparticle Physics 34 (2011) 476.
 J. Zabierowski, K. Daumiller and P. Doll, Nuclear Physics B (Proc. Suppl.) 122 (2003) 275.
 J. Zabierowski, P. Łuczak, P. Doll et al. (KASCADEGrande Collaboration), in Proceedings of the 32 ICRC, vol. 1, p. 220, Beijing, China, 2011.
 D. Heck et al., Report FZKA 6019, Forschungszentrum Karlsruhe, 1998.
 S.S. Ostapchenko, Nucl. Phys. B (Proc. Suppl.) 151 (2006) 143&147; S. Ostapchenko, Phys. Rev. D 74 (2006) 014026.
 A. Fassò et al., Report CERN200510, INFN/TC05/11, SLACR773 (2005).
 S. Ostapchenko, Physical Review D 81 (2010) 114028; S. Ostapchenko, Physical Review D 83 (2011) 014018.

P. Łuczak, PhD thesis, National Centre for Nuclear Research, Poland, 2012
 T. Antoni et al. , Astroparticle Physics 24 (2005) 1â25.