synthpop.modules.population_density.besancon2003_halo

Subclass for a spheroidal density distribution, used e.g. for the Halo in the Besancon Model (Robin et al., 2003).

Classes

Besancon2003Halo

Spheroidal density profile:

Module Contents

class synthpop.modules.population_density.besancon2003_halo.Besancon2003Halo(p0=9320.0, e=0.76, ac=0.5, power=-2.44, **kwargs)

Bases: synthpop.modules.population_density._population_density.PopulationDensity

Spheroidal density profile: the density is constant within a radius of ac, and then drops using a power law

Attributes

p0:

density at the position of the sun

efloat

Flattening ratio between x (or y) and z axis

ac: float

Radius until the spheroid has a constant density

power: float

Slope of the decay. Should be less than 0

population_density_name = 'Spheroid'

density_unit = 'mass'

p0 = 9320.0

e = 0.76

ac = 0.5

power = -2.44

density(r: numpy.ndarray, phi_rad: numpy.ndarray, z: numpy.ndarray) → numpy.ndarray

Estimates the density at the given position

Parameters

rndarray [‘kpc’]

Distance to z axis

phi_radndarray [‘rad’]

azimuth angle of the stars. phi_rad = 0 is pointing towards sun.

z : height above the galactic plane (corrected for warp of the galaxy)

Returns

rhondarray [M_sun/kpc^3 or #/kpc^3]

density at the given location, either in number density evolved mass density or initial mass density should be specified in density_unit.