Note

You are not reading the most up to date version of Gammapy documentation.
Access the latest stable version v1.3 or the list of Gammapy releases.

EinastoProfile

class gammapy.astro.darkmatter.profiles.EinastoProfile(r_s=None, alpha=None, rho_s=<Quantity 1. GeV / cm3>)

Bases: gammapy.astro.darkmatter.profiles.DMProfile

Einasto Profile.

$$\rho(r) = \rho_s \exp{ \left(-\frac{2}{\alpha}\left[ \left(\frac{r}{r_s}\right)^{\alpha} - 1\right] \right)}$$

Parameters:

r_s : Quantity

Scale radius, $r_s$

alpha : Quantity

$\alpha$

rho_s : Quantity

Characteristic density, $\rho_s$

References

• 1965TrAlm…5…87E
• 2011JCAP…03..051

Attributes Summary

DEFAULT_ALPHA	Default scale radius as given in reference 2
DEFAULT_SCALE_RADIUS	Default scale radius as given in reference 2
DISTANCE_GC
LOCAL_DENSITY

Methods Summary

__call__(radius)	Call evaluate method of derived classes.
evaluate(radius, r_s, alpha, rho_s)
integral(rmin, rmax, **kwargs)	Integrate squared dark matter profile numerically.
scale_to_local_density()	Scale to local density.

Attributes Documentation

DEFAULT_ALPHA = 0.17

Default scale radius as given in reference 2

DEFAULT_SCALE_RADIUS = <Quantity 28.44 kpc>

Default scale radius as given in reference 2

DISTANCE_GC = <Quantity 8.33 kpc>

LOCAL_DENSITY = <Quantity 0.3 GeV / cm3>

Methods Documentation

__call__(radius)

Call evaluate method of derived classes.

static evaluate(radius, r_s, alpha, rho_s)

integral(rmin, rmax, **kwargs)

Integrate squared dark matter profile numerically.

$$F(r_{min}, r_{max}) = \int_{r_{min}}^{r_{max}}\rho(r)^2 dr$$

Parameters:

rmin, rmax : Quantity

Lower and upper bound of integration range.

**kwargs : dict

Keyword arguments passed to integrate_spectrum()

scale_to_local_density()

Scale to local density.