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_sQuantity

Scale radius, \(r_s\)

alphaQuantity

\(\alpha\)

rho_sQuantity

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__(self, radius)	Call evaluate method of derived classes.
evaluate(radius, r_s, alpha, rho_s)
integral(self, rmin, rmax, \*\*kwargs)	Integrate squared dark matter profile numerically.
scale_to_local_density(self)	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__(self, radius)

Call evaluate method of derived classes.

static evaluate(radius, r_s, alpha, rho_s)

integral(self, 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, rmaxQuantity

Lower and upper bound of integration range.

**kwargsdict

Keyword arguments passed to integrate_spectrum()

scale_to_local_density(self)

Scale to local density.