alpha_prof

clustertools.analysis.profiles.alpha_prof(cluster, mmin=None, mmax=None, nmass=10, rmin=None, rmax=None, nrad=20, vmin=None, vmax=None, emin=None, emax=None, kwmin=0, kwmax=1, npop=None, indx=None, bins=None, projected=False, normalize=False, r_lower=None, r_upper=None, aerror=False, mcorr=None, plot=False, **kwargs)

Measure the radial variation in the mass function

• measure the delta_alpha parameter from Webb & Vesperini 2016

• Webb, J.J. & Vesperini, E. 2016, MNRAS, 463, 2383

Parameters:

clusterclass

StarCluster

mmin/mmaxfloat

minimum and maximum stellar mass

nmassint

number of mass bins (default: 10)

rmin/rmaxfloat

minimum and maximum stellar radii

nradint

number of radial bins

vmin/vmaxfloat

minimum and maximum stellar velocity

emin/emaxfloat

minimum and maximum stellar energy

kwmin/kwmaxfloat

minimum and maximum stellar type (kw)

npopint

population number

indxfloat

user defined boolean array from which to extract the subset

binsfloat

User defined bins in the form of (rlower,rmean,rupper) (default: None)

projectedbool

use projected values and constraints (default:False)

normalizebool

normalize radial bins by cluster’s half-mass radius (default: True)

r_lowerfloat

lower limits to radial bins

r_upperfloat

upper limits to radial bins

aerrorbool

return error in alpha calculations (default:True)

mcorrfloat

completeness correction for masses (default: None)

plotbool

plot the alpha profile (default: False)

Returns:

rprofnfloat

natural log of each radius bin (normalized by half-mass radius)

aproffloat

slope of the mass function in each bin

dalphafloat

radial variation in alpha calculated as delta_alpha = d(alpha)/d(ln(r/rm)

edalphafloat

error in dalpha

ydalphafloat

y-intercept in fit to alpha vs ln(r/rm)

eydalphafloat

error in ydalpha

rbinerrorfloat

if mcorr is not None, output lowest corrected mass fraction at each radius

Other Parameters:

kwragsstr

key word arguments for plotting