EinastoProfile

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

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

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)[source]
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, rmax : Quantity

Lower and upper bound of integration range.

**kwargs : dict

Keyword arguments passed to integrate_spectrum()

scale_to_local_density(self)

Scale to local density.