synthpop.modules.population_density.besancon2003_thindisk
Einasto density profiles, used e.g. for the thin disk in the Besancon Model (Robin et al., 2003).
Classes
Einasto profile |
Module Contents
- class synthpop.modules.population_density.besancon2003_thindisk.Besancon2003Thindisk(e, p0, hrp, hrm, offset=0, power=1, disk_cutoff=14, flare_flag=True, **kwargs)
Bases:
synthpop.modules.population_density._population_density.PopulationDensityEinasto profile
Attributes
- efloat
disc axis ratio
- p0float [Msun/kpc**3]
density at the position of sun
- hrpfloat [kpc]
scale radius of the positive Gaussian
- hrmfloat [kpc]
hrm is the scale radius of the disk hole
- offsetfloat
modification of the exponent
- powerfloat
modification of the exponent for the Besancon thin disk: offset = 0; power = 2 for age<=0.15 Gyr offset = 0.5; power = 1 for age> 0.15 Gyr
- disk_cutofffloat [kpc]
outer edge of the disc
- flare_flagbool
include flare
- population_density_name = 'Besancon2003Thindisk'
- density_unit = 'mass'
- disk_cutoff = 14
- p0
- e
- hrp
- hrm
- offset = 0
- power = 1
- flare_flag = True
- 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.