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

Neutrino telescopes are opening new opportunities in observational high energy astrophysics. In these detectors, atmospheric muons from primary cosmic ray interactions in the atmosphere play an important role. They provide the most abundant source of events for calibration and for testing the reconstruction algorithms. On the other hand, they represent the major background source. The simulation of a statistically significant number of muons in large volume neutrino telescopes requires a big effort in terms of computing time. Some parameterizations are currently available, but they do not explicitly take into account the arrival of muons in bundles. A fast Monte Carlo generator (MUPAGE) was developed to generate single and multiple atmospheric muon events in underwater/ice neutrino telescopes. The code reduces the computing time for the simulation of atmospheric muons significantly. The event kinematics is produced on the surface of a user-defined cylinder, virtually surrounding the detector volume. The flux of muon bundles at different depths and zenith angles, the lateral spread and the energy spectrum of the muons in the bundles are based on parametric formulas. These formulas were obtained according to a specific primary cosmic ray flux model and constrained by the measurements of the muon flux in the MACRO underground experiment. Some MUPAGE applications are presented.

More than 30 million of high-energy muons collected with the MACRO detector at the underground Gran Sasso Laboratory have been used to search for flux variations of different natures. Two kinds of studies were carried out: search for periodic variations and for the occurrence of clusters of events. Different analysis methods, including Lomb-Scargle spectral analysis and Scan Test statistics have been applied to the data.

We present the result of a search for correlations in the arrival times of high energy muons collected from 1995 till 2000 with the streamer tube system of the complete MACRO detector at the underground Gran Sasso Lab. Large samples of single muons (8.6 million), double muons (0.46 million) and multiple muons with multiplicities from 3 to 6 (0.08 million) were selected. These samples were used to search for time correlations of cosmic ray particles coming from the whole upper hemisphere or from selected space cones. The results of our analyses confirm with high statistics a random arrival time distribution of high energy cosmic rays.