galpy.df.quasiisothermaldf.tilt

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

NAME:

tilt

PURPOSE:

calculate the tilt of the velocity ellipsoid by marginalizing over velocity

INPUT:

R - radius at which to calculate this (can be Quantity)

z - height at which to calculate this (can be Quantity)

OPTIONAL INPUT:

nsigma - number of sigma to integrate the velocities over

scipy.integrate.tplquad kwargs epsabs and epsrel

mc= if True, calculate using Monte Carlo integration

nmc= if mc, use nmc samples

gl= if True, calculate using Gauss-Legendre integration

ngl= if gl, use ngl-th order Gauss-Legendre integration for each dimension

OUTPUT:

tilt in rad

HISTORY:

2012-12-23 - Written - Bovy (IAS)

2017-10-28 - Changed return unit to rad - Bovy (UofT)