meq_function¶
- clustertools.analysis.functions.meq_function(cluster, mmin=None, mmax=None, nmass=10, rmin=None, rmax=None, vmin=None, vmax=None, emin=None, emax=None, kwmin=0, kwmax=1, npop=1, indx=None, projected=False, plot=False, **kwargs)[source]¶
NAME: Find meq from velocity dispersion versus mass
mass bins are set up so that there are an equal number of stars in each bin
As per Bianchini, P. et al. 2016, MNRAS, 458, 3644, velocity dispersion versus mass is fit with the following: sigma(m)= sigma e^(-1/2 m/meq) if m<= meq
= sigma0 e^(-1/2) (m/meq)^-1/2 if m > meq
- Parameters:
- clusterclass
StarCluster instance
- mmin/mmaxfloat
specific mass range
- nmass
number of mass bins used to calculate alpha
- rmin/rmax
specific radial range
- vmin/vmaxfloat
specific velocity range
- emin/emaxfloat
specific energy range
- kwmin/kwmaxint
specific stellar evolution type range
- npopint
population number
- indxbool
specific subset of stars
- projectedbool
use projected values (default: False)
- plotbool
plot the mass function
- Returns:
- m_meanfloat
mean mass in each bin
- sigvmfloat
velocity dispersion of stars in each bin
- meqfloat
Bianchini fit to sigvm vs m
- emeqfloat
error in Bianchini fit to sigvm vs m
- sigma0float
Bianchini fit to sigvm vs m
- esigma0float
error in Bianchini fit to sigvm vs m
- Other Parameters:
- kwargsstr
key words for plotting