The Virtual Monte Carlo

The concept of Virtual Monte Carlo (VMC) has been developed by the ALICE Software Project to allow different Monte Carlo simulation programs to run without changing the user code, such as the geometry definition, the detector response simulation or input and output formats. Recently, the VMC classes have been integrated into the ROOT framework, and the other relevant packages have been separated from the AliRoot framework and can be used individually by any other HEP project. The general concept of the VMC and its set of base classes provided in ROOT will be presented. Existing implementations for Geant3, Geant4 and FLUKA and simple examples of usage will be described.Comment: Talk from the 2003 Computing in High Energy and Nuclear Physics (CHEP03), La Jolla, Ca, USA, March 2003, 8 pages, LaTeX, 6 eps figures. PSN THJT006. See http://root.cern.ch/root/vmc/VirtualMC.htm

Search for a simultaneous signal from small transient events in the Pierre Auger Observatory and the Tupi muon telescopes

We present results of a search for a possible signal from small scale solar transient events (such as flares and interplanetary shocks) as well as possible counterparts to Gamma-Ray Burst (GRB) observed simultaneously by the Tupi muon telescope Niteroi-Brazil, 22.90S, 43.20W, 3 m above sea level) and the Pierre Auger Observatory surface detectors (Malargue-Argentina, 69.30S, 35.30W, altitude 1400 m). Both cosmic ray experiments are located inside the South Atlantic Anomaly (SAA) region. Our analysis of several examples shows similarities in the behavior of the counting rate of low energy (above 100 MeV) particles in association with the solar activity (solar flares and interplanetary shocks). We also report an observation by the Tupi experiment of the enhancement of muons at ground level with a significance higher than 8 sigma in the 1-sec binning counting rate (raw data) in close time coincidence (T-184 sec) with the Swift-BAT GRB110928B (trigger=504307). The GRB 110928B coordinates are in the field of view of the vertical Tupi telescope, and the burst was close to the MAXI source J1836-194. The 5-min muon counting rate in the vertical Tupi telescope as well as publicly available data from Auger (15 minutes averages of the scaler rates) show small peaks above the background fluctuations at the time following the Swift-BAT GRB 110928B trigger. In accordance with the long duration trigger, this signal can possibly suggest a long GRB, with a precursor narrow peak at T-184 sec.Comment: 9 pages, 13 figure

The physics models of FLUKA: status and recent development

A description of the intermediate and high energy hadronic interaction models used in the FLUKA code is given. Benchmarking against experimental data is also reported in order to validate the model performances. Finally the most recent developments and perspectives for nucleus-nucleus interactions are described together with some comparisons with experimental data.Comment: talk from the 2003 Computing in High Energy and Nuclear Physics (CHEP03), La Jolla, Ca, USA, March 2003, 10 pages, p

Calculation of atmospheric neutrino flux using the interaction model calibrated with atmospheric muon data

Using the ``modified DPMJET-III'' model explained in the previous paper, we calculate the atmospheric neutrino flux. The calculation scheme is almost the same as HKKM04 \cite{HKKM2004}, but the usage of the ``virtual detector'' is improved to reduce the error due to it. Then we study the uncertainty of the calculated atmospheric neutrino flux summarizing the uncertainties of individual components of the simulation. The uncertainty of $K$-production in the interaction model is estimated by modifying FLUKA'97 and Fritiof 7.02 so that they also reproduce the atmospheric muon flux data correctly, and the calculation of the atmospheric neutrino flux with those modified interaction models. The uncertainties of the flux ratio and zenith angle dependence of the atmospheric neutrino flux are also studied

First observations of beam losses due to bound-free pair production in a heavy-ion collider

We report the first observations of beam losses due to bound-free pair production at the interaction point of a heavy-ion collider. This process is expected to be a major luminosity limit for the Large Hadron Collider (LHC) when it operates with 208Pb82+ ions because the localized energy deposition by the lost ions may quench superconducting magnet coils. Measurements were performed at the Relativistic Heavy Ion Collider (RHIC) during operation with 100 GeV/nucleon 63Cu29+ ions. At RHIC, the rate, energy and magnetic field are low enough so that magnet quenching is not an issue. The hadronic showers produced when the single-electron ions struck the RHIC beampipe were observed using an array of photodiodes. The measurement confirms the order of magnitude of the theoretical cross section previously calculated by others.Comment: 4 pages, 5 figures. Added journal ref. Corrected typos. Fixed fig 1. Minor improvements to fig. 1,3,4. Rephrased a small number of sentences (p1,3,4). Added numerical values of the aperture and the displacement for Au (p 2). Changed reference 5, added name in acknowledgments (p 4

FLUKA and PENELOPE simulations of 10 keV to 10 MeV photons in LYSO and soft tissue

n/

Heavy-ion collisions described by a new QMD code interfaced to FLUKA: model validation by comparisons with experimental data concerning neutron and charged fragment production

A new code, based on the Quantum Molecular Dynamics theoretical approach, has been developed and interfaced to the FLUKA evaporation/fission/Fermi break-up module. At present, this code is undergoing a series of validation tests. In this paper its predictions are compared to measured charged fragment yields and double differential neutron spectra in thin target heavy-ion reactions, at bombarding energies of about 100 MeV/A. The comparisons with the predictions of a modified version of RQMD 2.4 originally developed in Frankfurt, already available in the FLUKA code, are presented and potential improvements are briefly sketched.A new code, based on the Quantum Molecular Dynamics theoretical approach, has been developed and interfaced to the FLUKA evaporation/fission/Fermi break-up module. At present, this code is undergoing a series of validation tests. In this paper its predictions are compared to measured charged fragment yields and double differential neutron spectra in thin target heavy-ion reactions, at bombarding energies of about 100 MeV/A. The comparisons with the predictions of a modified version of RQMD 2.4 originally developed in Frankfurt, already available in the FLUKA code, are presented and potential improvements are briefly sketched

Monte Carlo calculations for the ATLAS cavern background

A new application for simulating the ATLAS cavern background was developed. This was done using FLUGG, software that allows Geant4 geometry to be used within the FLUKA simulation framework. A Geant4 description of the ATLAS detector including its cavern was built from scratch for this application. In order to gain computing performance, our geometry is less detailed than that of GeoModel which is used in the full detector simulation, but good enough for the investigation of cavern background. Our geometry can also be used in a standalone Geant4 simulation. Thus it is possible to perform unbiased comparisons between Geant4 and FLUKA using the same complex geometry. We compared neutron and photon fluxes using the FLUKA-FLUGG application with the result of Geant4 simulations based on the QGSP_BERT and QGSP_BERT_HP physics lists. In all cases the same set of initial collision 4-vectors produced by the PHOJET event generator was used. The result from the QGSP_BERT_HP physics list, which uses the High Precision (HP) neutron model, is similar to the result of FLUKA-FLUGG and the differences in the fluxes between them are within 40% in most regions of the ATLAS cavern. The result from the QGSP_BERT physics list, which does not include the HP model, does not agree with either of the previous two results

Charm Production in DPMJET

In this work, charm production in the {\sc dpmjet} hadronic jet simulation is compared to experimental data. Since the major application of {\sc dpmjet} is the simulation of cosmic ray-induced air showers, the version of the code integrated in the CORSIKA simulation package has been used for the comparison. Wherever necessary, adjustments have been made to improve agreement between simulation and data. With the availability of new muon/neutrino detectors that combine a large fiducial volume with large amounts of shielding, investigation of prompt muons and neutrinos from cosmic ray interactions will be feasible for the first time. Furthermore, above $\gtrsim 100$ TeV charmed particle decay becomes the dominant background for diffuse extraterrestrial neutrino flux searches. A reliable method to simulate charm production in high-energy proton-nucleon interactions is therefore required.Comment: 10 pages, to be published in JCA

The influence of the geomagnetic field and of the uncertainties in the primary spectrum on the development of the muon flux in the atmosphere

In this paper we study the sensitivity of the flux of atmospheric muons to uncertainties in the primary cosmic ray spectrum and to the treatment of the geomagnetic field in a calculation. We use the air shower simulation program AIRES to make the calculation for two different primary spectra and under several approximations to the propagation of charged particles in the geomagnetic field. The results illustrate the importance of accurate modelling of the geomagnetic field effects. We propose a high and a low fit of the proton and helium fluxes, and calculate the muon fluxes with these different inputs. Comparison with measurements of the muon flux by the CAPRICE experiment shows a slight preference for the higher primary cosmic ray flux parametrization.Comment: 24 pages, 13 figures, submitted to Phys.Rev.