/CRFluxModels

Models/parametrizations of the cosmic ray flux at the top of the atmosphere

Primary LanguageJupyter NotebookMIT LicenseMIT

CRFluxModels --- models of the high-energy cosmic ray flux

Historically, this module was part of the research code for a paper A. Fedynitch, J. Becker Tjus, and P. Desiati, Phys. Rev. D 86, 114024 (2012), where we compared the effects of different primary Cosmic Ray Flux models - therefore the name. Later the module became an integral part of my research tools and has been recently cleaned up and documented. It

  • implements numerical models of high energy cosmic ray fluxes,
  • converts from all-particle into all-nucleon flux,
  • includes some convenience functions for semi-analytical atmospheric lepton flux calculations.

Since this module never was a performance bottle-neck, the formulae are written in a human-readable awy without any attempt to optimize for speed or elegance. If you find it slow and improve something, feel free to send a pull request or just branch/fork and let me know.

The module is used as a submodule in the Matrix Cascade Equation app.

Documentation

Please follow this link to the documentation.

Installation

The module was developed using Python2.7 and NumPy. It doesn't use any kind of fancy functionality. It should therefore work with all versions. Some models require SciPy.

Example

Run python CRFluxModels.py from a shell to generate a set of standard plots.