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 |
File attributes
Modifcation date: | Fri Jun 6 08:19:26 2008 |
Lines: | 275 |
Docformat: | rst rst |