1 / 24

G. Battistoni, A. Margiotta, S. Muraro , M. Sioli (University and INFN of Bologna and Milano)

The FLUKA high energy cosmic ray generator: predictions for the charge ratio of muons detected underground. G. Battistoni, A. Margiotta, S. Muraro , M. Sioli (University and INFN of Bologna and Milano) for the FLUKA collaboration. 44 th Rencontres de Moriond

brandy
Download Presentation

G. Battistoni, A. Margiotta, S. Muraro , M. Sioli (University and INFN of Bologna and Milano)

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. The FLUKA high energy cosmic ray generator:predictions for the charge ratio of muons detected underground G. Battistoni, A. Margiotta, S. Muraro, M. Sioli (University and INFN of Bologna and Milano) for the FLUKA collaboration 44th Rencontres de Moriond Very High Energy Phenomena in the Universe

  2. In the environment of the FLUKA Monte Carlo code application to cosmic rays physics, a new generator for high energy cosmic rays is under development, with the aim of extend the existing FLUKA cosmic rays library to include the TeV region. • The application of FLUKA in cosmic ray physics arises from the interest inapplied physics topics (radioprotection in space or in atmosphere) and in basic research (calculation of atmospheric neutrino fluxes). • Generator dedicated to: • physics of high energy underground muons • exploiting the full integration in the calculation of both air shower development and muons transport in the rock. • Aim: predict multiple muon rates for different primary masses and energy within the framework of a unique simulation model.

  3. Work under way within the ICARUS and OPERA collaborations at Gran Sasso. First application: Analyze the predictions for the charge ratio of underground muons. Compare preliminary results with data from an ongoing experiment (MINOS).

  4. Authors: A. Fasso`1, A. Ferrari2, J. Ranft3, P.R. Sala4 1 SLAC Stanford, 2 CERN, 3 Siegen University, 4 INFN Milan Interaction and Transport Monte Carlo code FLUKA http://www.fluka.org • FLUKA is a general purpose tool for calculations of particle transport and interactions with matter, covering an extended range of applications (Shielding, Radiobiology, High energy physics, Cosmic Ray physics, Nuclear and reactor physics). • Built and maintained with the aim of including the best possible physical models in terms of completeness and precision. • Continuously benchmarked with a wide set of experimental data from well controlled accelerator experiments. Hadronic interaction models: based on a theoretical microscopic approach (no parametrizations). Free parameters are set (thin target experiments at accelerators) and kept fixed for all projectile-target combinations and energies. => High predictivity also in regions where experimental data are not available.

  5. Relevant for HE C.R. physics The FLUKA hadronic Modelsfor a detailed study of the validity for CR studies see hep-ph/0612075 and 0711.2044 > 5 GeV Elab DPM:soft physics based on (multi)Pomeron exchange DPMJET:soft physics of DPM plus 2+2 processes from pQCD

  6. The underground muons generator: main features • Geometry setup • Earth: • sphere of radius R = 6378.14 km • Atmospheric geometry & profile: • 100 concentric spherical shells whose density and composition is varied according to the U.S. Standard Atmospheric Model. • Gran Sasso mountain: • spherical body whose radius is dynamically changed, according to the primary direction and to the Gran Sasso mountain map. • LNGS laboratory: • experimental underground halls • ICARUS and OPERA detectors volumes • rock box where muon–induced secondary are activated • (e.m. & hadronic showers from photo-nuclear interaction).

  7. Primary spectrum Sampled from a primary mass composition model (a description of the relative abundances of cosmic rays and their energy spectra), at present, derived from the analysis of the MACRO experiment at Gran Sasso. For each primary nucleus and for each amount of rock to be crossed, we compute the minimum energy required to produce at least one muon underground (probability < 10-5 to survive). Hadronic interaction model High energy cosmic ray interactions relevant for this work are treated in FLUKA by means of the interface to DPMJET.II.5. (hadron-hadron, hadron-nucleus and also nucleus-nucleus collision by means of the Glauber-Gribov mechanism)

  8. Muon bundle from primary iron nuclei (E ≈ 105 TeV) in the ICARUS T600 detector

  9. First application: prediction for thecharge ratio of underground muons • Muons that reach the Earth • come from mesonswith enough energy: • to reflect the forward fragmentation region • of the primary initiated interaction and • to “remember” the nature of the projectile • (there are more protons than neutrons • in the primary spectrum) • The muon charge ratio reflects the excess of π+ over π- and K+ over K−. • NOTE: • π and K hadronic production are affected by • uncertainties up to 20%

  10. Validation of the DPMJET-III hadronic models: Comparison with the NA49 experiment Data from the NA49 experiment at CERN SPS particle production by p beams on p, C targets: 158 GeV/c beam momentum First published results: Eur. Phys. J. 45 (2006), 343 hep-ex/0606028 hep-ex/0606029 + , - production p + p p + C + , - production as a function of Feynman-x

  11. SPY experiment (CERN North Area) Nucl. Instr. Meth. A449, 609 (2000) p + Be => π+ + X p + Be => π- + X Ecm = 450 GeV

  12. Benchmark for the CNGS beam construction. Limited phase space for cosmic rays physics. Nucl. Instr. Meth. A449, 609 (2000) p + Be => K+ + X p + Be => K- + X Ecm = 450 GeV

  13. FLUKA for Cosmic Rays validationat low energy (Eμ< 100 GeV) • FLUKA simulations comparison with the experimental data of atmospheric muons taken at the top of Mt. Norikura, Japan, with the BESS detector. (Phys. Lett. B 564 (2003), 8 – 20) • 2770 m above sea level • Geomagnetic Cut-off: 11.2 GV • cone of ~11o • The energy range for muons extends up to 100 GeV. Results within 20%. S.Muraro PhD Thesis Black points: exp. Data Open symb: FLUKA Simulation m- m+

  14. The muons result from pions and kaons that decay before they interact in the atmosphere. Instead of π production, because of the strangeness of K, inclusive cross section for K+ production is bigger then inclusive cross section for K- → the K+/K− ratio is larger than the π+/π− ratio. As energy increases, the fraction of muons from kaon decays also increases: the longer-lived pions (π± : cτ0= 780 cm, ε = 115 GeV) start to interact more before decaying than the shorter-lived kaons (K± : cτ0 = 371 cm, ε = 850 GeV). critical energy ε: beyond this energy interaction process dominates on decay.

  15. FLUKA in the TeV region:μ/All fromπand K μ(from π) / All μ (from K) / All

  16. As energy increases, kaon decays became a more important contribution to the muon charge ratio. Since Nμ+ Nμ+ (fromK) > (fromπ) Nμ- Nμ− the total muons charge ratio is expected to increase with energy

  17. FLUKA for Cosmic Rays validation (Eμ < 1 TeV) FLUKA simulations comparison with the experimental data of atmospheric muons charge ratio from L3+COSMIC experiment (hep-ex/0408114). Vertical 0.975 < cosθ < 1. RFLUKA = 1.295  0.0482 Rexp = 1.285  0.484 Black points: exp. Data Open symb: FLUKA At large angle 0.525 < cosθ < 0.6 (S.Muraro PhD Thesis)

  18. FLUKA in the TeV region:muon charge ratio from K fromπ μ+/ μ-FLUKA 1.312  0.0155 μ+/ μ-FLUKAfromπ 1.22017 0.0163948 μ+/ μ-FLUKAfrom K 1.78605 0.0458639 PRELIMINARY

  19. MINOS Charge Ratio at the Surface = 1.371± 0.003 hep-ex 0705.3815 RFLUKA = 1.295  0.0482 RFLUKA = 1.312  0.0155 FLUKA prediction ~5% lower then MINOS exp. data => Possible lack of K production in FLUKA

  20. Conclusions FLUKA models have been benchmarked with experimental data from accelerator experiments and from atmospheric muons experiments (BESS: Eμ< 100 GeV; L3+C: Eμ < 1 TeV) The FLUKA charge ratio prediction in the TeV region is ~5% lower then MINOS experimental results Uncertenties on π and K production in hadronic models can reach 20% because of lack of data from K production experiments at high energy => Possible lack of K production in FLUKA We are waiting for OPERA data

  21. Thank you

  22. Muon charge ratio VS muon bundle multiplicity In the primary heavy elements the ratio of primary protons to neutrons decreases with respect to primary protons → the muon charge ratio is expected to decrease. Muon bundle high multiplicity ↕ High primary energy and High primary mass number Muon charge ratio is expected decreases with growing multiplicity

  23. xlab = Ej/Ei ratio of the total energies of the secondary particle j over the primary particle i dNij/dxlab differential multiplicity distributions of secondary j as produced by primary i in collisions with air nuclei as a function of xlab ”spectrum weighted moments” Zij : the multiplicity of secondary particles j as produced by primary particles i in interaction, weighted for the primary spectrum. Strictly bound to inclusive cross sections. γ= 2.7 approximate spectral index of the differential cosmic ray spectrum.

  24. For isospin symmetry: On the other hand: where N is a nucleon. K+ and 0 (S = +1), can be produced in association with a leading Λ or Σ barion, whereas production of K requires production of a strange-antistrange pair from the sea in addition to the leading nucleon So the K+/K− ratio is larger than the π+/π− ratio. Spectrum weighted moments (γ = 2.7) for secondary particles produced in p-air collisions as a function of the projectile kinetic energy in the FLUKA code.

More Related