Anthony Smith's IDL routines

Galaxies etc.

single page | use frames     summary     class     fields     routine details     file attributes

gals/

ajs_vmax_lf.pro

This procedure calculates the 1/Vmax luminosity function from input arrays of absolute magnitude and weights (e.g., 1/Vmax)

Routines

result = ajs_vmax_lf(absmag [, weights] [, bincentres=bincentres] [, err_phi=err_phi] [, ninbin=ninbin] [, nbins=nbins] [, xrange=xrange] [, jackknife=jackknife] [, lum_dens=lum_dens] [, err_lum_dens=err_lum_dens] [, loglum=loglum] [, schechter=schechter] [, cov_mat=cov_mat] [, sch_range=sch_range] [, sch_lum_dens=sch_lum_dens] [, sch_err_lum_dens=sch_err_lum_dens] [, double_schechter=double_schechter] [, double_cov_mat=double_cov_mat] [, double_sch_lum_dens=double_sch_lum_dens] [, double_sch_err_lum_dens=double_sch_err_lum_dens])

This procedure calculates the 1/Vmax luminosity function from input arrays of absolute magnitude and weights (e.g., 1/Vmax)

ajs_vmax_lf_test

Test ajs_vmax_lf

Routine details

topajs_vmax_lf

result = ajs_vmax_lf(absmag [, weights] [, bincentres=bincentres] [, err_phi=err_phi] [, ninbin=ninbin] [, nbins=nbins] [, xrange=xrange] [, jackknife=jackknife] [, lum_dens=lum_dens] [, err_lum_dens=err_lum_dens] [, loglum=loglum] [, schechter=schechter] [, cov_mat=cov_mat] [, sch_range=sch_range] [, sch_lum_dens=sch_lum_dens] [, sch_err_lum_dens=sch_err_lum_dens] [, double_schechter=double_schechter] [, double_cov_mat=double_cov_mat] [, double_sch_lum_dens=double_sch_lum_dens] [, double_sch_err_lum_dens=double_sch_err_lum_dens])

This procedure calculates the 1/Vmax luminosity function from input arrays of absolute magnitude and weights (e.g., 1/Vmax)

Bins must have equal width

Return value

fltarr/dblarr Returns phi, value of luminosity function

Parameters

absmag in required
weights in optional

E.g., 1./Vmax. Set to a single value to give all galaxies that weight, or leave blank to give all galaxies a weight of 1.

Keywords

bincentres in out optional

Centre of each bin in absolute magnitude

err_phi out optional

Poisson errors = phi / sqrt(n)

ninbin out optional

Number of galaxies in each bin

nbins in optional

Number of bins (ignored if bincentres set)

xrange in optional

Minimum and maximum absmag (ignored if bincentres set)

jackknife in optional

Integer for each galaxy giving the jackknife region in which the galaxy lies

lum_dens out optional

Luminosity density. Return value is sum(lum * weight), j = ... x 10^(0.4 M_sun) h L_sun Mpc^-3

err_lum_dens out optional

Error on luminosity density, calculated from errors on phi, or from jackknife resampling

loglum in optional

For Schechter fits: /loglum indicates log(luminosity) inputs

schechter out optional

Returns best-fitting Schechter function (using jackknife if present)

cov_mat out optional

Returns covariance matrix of Schechter function fit (using jackknife if present)

sch_range in optional

Range of absolute magnitudes for (double) Schechter function fit

sch_lum_dens out optional

As lum_dens, but for the Schechter function fit

sch_err_lum_dens out optional

Error on Schechter luminosity density

double_schechter out optional

Returns best-fitting double Schechter function (using jackknife if present)

double_cov_mat out optional

Returns covariance matrix of double Schechter function fit (using jackknife if present)

double_sch_lum_dens out optional

As lum_dens, but for the double Schechter function fit

double_sch_err_lum_dens out optional

Error on double Schechter luminosity density

Author information

History:

5 Sep 2007 Created (Anthony Smith)

6 Mar 2008 Re-written

18 Mar 2008 Added jackknife

7 Apr 2008 Jackknife estimation of Schechter function fit

15 Apr 2008 Added luminosity density

21 Apr 2008 Added double Schechter function

topajs_vmax_lf_test

ajs_vmax_lf_test

Test ajs_vmax_lf

File attributes

Modifcation date: Fri Jun 6 08:19:26 2008
Lines: 275
Docformat: rst rst