cosmic_ray_flux

gammapy.spectrum.cosmic_ray_flux(energy, particle='proton')[source]

Cosmic ray flux at Earth.

These are the spectra assumed in this CTA study: Table 3 in http://adsabs.harvard.edu/abs/2013APh….43..171B

The hadronic spectra are simple power-laws, the electron spectrum is the sum of a power law and a log-normal component to model the “Fermi shoulder”.

Parameters:

energy : Quantity

Particle energy

particle : {‘electron’, ‘proton’, ‘He’, ‘N’, ‘Si’, ‘Fe’}

Particle type

Returns:

flux : Quantity

Cosmic ray flux in unit m^-2 s^-1 TeV^-1 sr^-1