galpy.df.quasiisothermaldf.sigmaT2

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

Calculate sigma_T^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 is True, use nmc samples.

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

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

  • **kwargs – scipy.integrate.tplquad kwargs epsabs and epsrel.

Returns:

sigma_T^2.

Return type:

float

Notes

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