FLUKA is a multipurpose Monte Carlo code, which can transport particles over a wide range of energies in user-defined geometries. Here we present a new FLUKA library, which allows the interaction and propagation of high energy cosmic rays in the Earth atmosphere and the transport of high energy muons in underground/underwater environments

Primary GCR interact with the Earth's atmosphere originating atmospheric showers, thus giving rise to fluxes of secondary particles in the atmosphere. Electromagnetic and hadronic interactions interplay in the production of these particles, whose detection is performed by means of complementary techniques in different energy ranges and at different depths in the atmosphere, down to the Earth's surface. Monte Carlo codes are essential calculation tools which can describe the complexity of the physics of these phenomena, thus allowing the analysis of experimental data. However, these codes are affected by important uncertainties, concerning, in particular, hadronic physics at high energy. In this paper we shall report some results concerning inclusive particle fluxes and atmospheric shower properties as obtained using the FLUKA transport and interaction code. Some emphasis will also be given to the validation of the physics models of FLUKA involved in these calculations.

The 3-dimensional (3-D) calculation of the atmospheric neutrino flux by means of the FLUKA Monte Carlo model is here described in all details, starting from the latest data on primary cosmic ray spectra. The importance of a 3-D calculation and of its consequences have been already debated in a previous paper. Here instead the focus is on the absolute flux. We stress the relevant aspects of the hadronic interaction model of FLUKA in the atmospheric neutrino flux calculation. This model is constructed and maintained so to provide a high degree of accuracy in the description of particle production. The accuracy achieved in the comparison with data from accelerators and cross checked with data on particle production in atmosphere certifies the reliability of shower calculation in atmosphere. The results presented here can be already used for analysis by current experiments on atmospheric neutrinos. However they represent an intermediate step towards a final release, since this calculation does not yet include the bending of charged particles in atmosphere. On the other hand this last aspect, while requiring a considerable effort in a fully 3-D description of the Earth, if a high level of accuracy has to be maintained, does not affect in a significant way the analysis of atmospheric neutrino events.

Substantial fluxes of protons and leptons with energies below the geomagnetic cutoff have been measured by the AMS experiment at altitudes of 370-390 Km, in the latitude interval +/- 51.7 degrees. The production mechanisms of the observed trapped fluxes are investigated in detail by means of the FLUKA Monte Carlo simulation code. All known processes involved in the interaction of the cosmic rays with the atmosphere (detailed descriptions of the magnetic field and atmospheric density, as well as the electromagnetic and nuclear interaction processes) are included in the simulation. The results are presented and compared with the experimental data, indicating good agreement with the observed fluxes. The impact of secondary proton flux on particle production in atmosphere is briefly discussed.

In view of the interpretation of data collected by large, deep underground experiments in terms of primary Cosmic Ray physics, this work is focused on the study of the production of of TeV muons in Extensive Air Showers. The review tries to point out those features of hadronic interactions that mostly affect the production of the high energy muons. A few different Monte Carlo codes are compared, with a particular attention to those based on the Regge-Gribov framework. The possibility of performing experimental tests of the proposed models is also briefly discussed.

We present first results about the measurement of the characteristics of charged hadrons production by atmospheric muons in the rock above MACRO. Selection criteria which allow to discriminate hadron cascades from e.m. showers generated by muons are described. A comparison between the measured rate with that expected from a Monte Carlo simulation which treats the process as dominated by photo-nuclear interaction is presented. These data can be used to validate such models aiming to the evaluation of hadron background from cosmic muons in different experimental environments.

A detailed simulation of vertical showers in atmosphere produced by primary gammas and protons, in the energy range 1-100 TeV, has been performed by means of the FLUKA Monte Carlo code, with the aim of studying the time structure of the shower front at different detector heights. It turns out that the time delay distribution can be fitted using few parameters coincident with the distribution central moments. Such parameters exhibit a smooth behaviour as a function of energy. These results can be used both for detector design and for the interpretation of the existing measurements. Differences in the time structure between gamma and proton induced showers are found and explained in terms of the non-relativistic component of extensive air showers.