Feature request: Analytical composition which is energy dependent
Implement a function in auger.py that sample charges for the cosmic rays as a function of the energy. I am thinking about 4 analytical fraction functions (H, He, CNO, Fe) as a function of energy, which is fitted to data. A similar approach was done in GAP2016_047 (Figure 5).