galpy.df.quasiisothermaldf.sigmaz2

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

Calculate sigma_z^2 by marginalizing over velocity.

Parameters:
  • R (float or Quantity) – Radius at which to calculate this.

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

  • nsigma (int, 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:

sigma_z^2.

Return type:

float

Notes

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