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_sQuantity

Scale radius, \(r_s\)

alphaQuantity

\(\alpha\)

rho_sQuantity

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__(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)[source]
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, rmaxQuantity

Lower and upper bound of integration range.

**kwargsdict

Keyword arguments passed to integrate_spectrum()

scale_to_local_density()

Scale to local density.