galpy.df.quasiisothermaldf.density

quasiisothermaldf.density(R, z, nsigma=None, mc=False, nmc=10000, gl=True, ngl=10, **kwargs)[source]

Calculate the density at R,z by marginalizing over velocity.

Parameters:

  • R (float or Quantity) – Radius at which to calculate the density.

  • z (float or Quantity) – Height at which to calculate the density.

  • nsigma (float, optional) – Number of sigma to integrate the velocities over.

  • mc (bool, optional) – If True, calculate using Monte Carlo integration.

  • nmc (int, optional) – If mc, use nmc samples.

  • gl (bool, optional) – If True, calculate using Gauss-Legendre integration.

  • ngl (int, optional) – If gl, use ngl-th order Gauss-Legendre integration for each dimension.

  • **kwargs (dict, optional) – scipy.integrate.tplquad kwargs epsabs and epsrel.

Returns:

Density at (R,z).

Return type:

float

Notes

  • 2012-07-26 - Written - Bovy (IAS@MPIA)