##################################################
#### Atmospheric E-M showers distributions   #####
##################################################
L. Darme
20-05-2022

## --- General description --- ##

This is the companion dataset to the publication "Atmospheric resonant production for light dark sectors" by the same author.
It includes the track-length fluxes for electrons, positrons and photons generated by cosmic ray showers, as well as light mesons fluxes obtained by the same process.

A simple python example ReadTL.py is additionally provided.

## --- Practical implementation  --- ##

The distributions  are provided in two data files which can be used "per se". Note that we have packaged the "track length flux"
from the main paper into two different files !

- AverageTrackLength_CR.dat: contains the differential track length fluxes (TLF) for photon, electron and positron, for which the cosmic ray flux has been
  integrated on a half hemisphere (2 Pi solid angle). This data is thus in (m^-2.s^-1).GeV^-1.cm where the (m^-2.s^-1) comes from the CR flux,
  the cm from the track lengths and the energy is in GeV.
- AverageRho_CR.dat: contains the average mass density (NOT normalised to the ground level atmospheric density 1.2e-3 g.cm-3). Note that the average is weighted by the track length of each
    track (so an electron propagating two times longer than another will contribute more) to reflect the fact that the density is used to obtain the number of potential targets

The product of the  above two distributions and subsequent normalisation by ground level atmospheric density (1.2e-3 g.cm-3) gives the normalised track length flux as defined in Eq. (1)
of the paper (see the example python loader ReadTL.py).

The mesons distributions are distributed in another file:

- MesonsDist_CR.dat: contains the differential meson flux for pi0 and eta mesons, for which the cosmic ray flux has been integrated on a half hemisphere (2 Pi solid angle)
  The data is thus in (m^-2.s^-1).GeV^-1 where the (m^-2.s^-1) comes from the CR flux and the energy is in GeV.


## --- Additional details  --- ##

The energy is always the full energy of the particle, not the kinetic energy. Thus the distributions for the mesons stop at the meson mass.


## --- Version history  --- ##

- 1.0.0: Public release alongside the ArXiV Publication