Astrophysical source population models#
gammapy.astro.population module provides a simple framework for population
synthesis of gamma-ray sources, which is useful in the context of surveys and
The following example illustrates how to simulate a basic catalog including a spiral arm model.
import astropy.units as u from gammapy.astro.population import make_base_catalog_galactic max_age = 1E6 * u.yr SN_rate = 3. / (100. * u.yr) n_sources = int(max_age * SN_rate) table = make_base_catalog_galactic( n_sources=n_sources, rad_dis='L06', vel_dis='F06B', max_age=max_age, spiralarms=True, )
The total number of sources is determined assuming a maximum age and a supernova
rate. The table returned is an instance of
can be used for further processing. The example population with spiral-arms is
illustrated in the following plot.
Galactocentric spatial distributions#
Here is a comparison plot of all available radial distribution functions of the surface density of pulsars and related objects used in literature:
TODO: add illustration of Galactocentric z-distribution model and combined (r, z) distribution for the Besancon model.
Spiral arm models#
Here is a comparison plot of all available velocity distribution functions: