CORSIKA implementation of heavy quark production and propagation in extensive air showers ?
Heavy quarks are commonly produced in current accelerator experiments. Hence it is natural to think that they should be likewise created in collisions with larger center of mass energies like the ones involving ultra-high energy cosmic rays and atmospheric nuclei. Despite this fact, a detailed treatment of heavy hadrons is missing in Monte Carlo generators of Extensive Air Showers (EAS). It is a must to improve the description of how heavy flavors appear and evolve in atmospheric showers. With this goal in mind, we study two different models for heavy quark production in protonair collisions. We also analyze a dedicated treatment of heavy hadrons interactions with atmospheric nuclei. This paper shows how those models have been implemented as new options available in CORSIKA, one of the most used EAS simulators. This new computational tool allows us to analyze the effects that the propagation of heavy hadrons has in the EAS development. Program title: corsika-6990-Heavy Catalogue identifier: AERC_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AERC_v1_0.html Program obtainable from: CPC Program Library, Queens University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 7681 No. of bytes in distributed program, including test data, etc.: 91906 Distribution format: tar.gz Programming language: FORTRAN77. Computer: Intel-Pentium based Personal Computers. Operating system: Linux. RAM: 373 Mb Classification: 1.1. External routines: Corsika (http://www-ik.fzk.de/~corsika/) Nature of problem: Charmed and bottom hadron production and propagation inside Extensive Air Showers. Solution method: Heavy quarks are produced according to two different production models. Propagation in the atmosphere is handled using already present CORSIKA subroutines. New subroutines are written to simulate their interactions with air nuclei. Restrictions: Heavy quark production has only been implemented in the first interaction. Running time: From 4.2h (output file size ?120Mb) at 1019 eV to 4.7h (output file size ?170Mb) at 1019.75 eV, with thinning 10?6 E(GeV).
corsika heavy quarks eas simulation