Anthony Smith's IDL routines

Galaxies etc.

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

gals/

ajs_choloniewski_fit.pro

This function finds the best-fitting Choloniewski function for input arrays of absolute magnitudes, surface brightness, phi and phi_err.

Routines

ajs_choloniewski_fit_phierr, absmag, sb, phi, phierr, vars=vars, zmin=zmin, zmax=zmax, area=area, q0=q0, q1=q1

Estimate phierr for bins with zero phierr

result = ajs_choloniewski_fit(absmag, sb, phi, phierr [, params_init] [, txtfile=txtfile] [, xrange=fltarr(2)] [, yrange=fltarr(2)] [, vars=vars] [, zmin=zmin] [, zmax=zmax] [, area=area] [, q0=q0] [, q1=q1] [, cov_mat=cov_mat] [, _REF_EXTRA=_REF_EXTRA])

This function finds the best-fitting Choloniewski function for input arrays of absolute magnitudes, surface brightness, phi and phi_err.

ajs_choloniewski_fit_test

Test ajs_choloniewski_fit

Routine details

topajs_choloniewski_fit_phierr

ajs_choloniewski_fit_phierr, absmag, sb, phi, phierr, vars=vars, zmin=zmin, zmax=zmax, area=area, q0=q0, q1=q1

Estimate phierr for bins with zero phierr

Parameters

absmag
sb
phi
phierr

Keywords

vars
zmin
zmax
area
q0
q1

topajs_choloniewski_fit

result = ajs_choloniewski_fit(absmag, sb, phi, phierr [, params_init] [, txtfile=txtfile] [, xrange=fltarr(2)] [, yrange=fltarr(2)] [, vars=vars] [, zmin=zmin] [, zmax=zmax] [, area=area] [, q0=q0] [, q1=q1] [, cov_mat=cov_mat] [, _REF_EXTRA=_REF_EXTRA])

This function finds the best-fitting Choloniewski function for input arrays of absolute magnitudes, surface brightness, phi and phi_err.

Return value

Choloniewski parameters [mstar,alpha,phistar,sbstar,sigmasb,beta]

Parameters

absmag in required

[n1] Bin centres for absolute magnitude

sb in required

[n2] Bin centres for surface brightness

phi in required

[n1, n2] Value of space density at absmag and sb

phierr in out required

[n1, n2] If any values of phierr are zero, either raises an error message or estimates (and changes input) values using vars etc.

params_init in optional

Initial guess for [mstar,alpha,phistar,sbstar,sigmasb,beta] (or default)

Keywords

txtfile in optional

Output text file for Choloniewski parameters

xrange in optional type=fltarr(2)

Range in m1 for fitting

yrange in optional type=fltarr(2)

Range in m2 for fitting

vars in optional

For estimating errors on bins with zero phi or phierr, for Choloniewski fit. Array of structures with fields {abs_name:} giving the column name in table corresponding to the absolute value, {band:} giving the name of the band, {type:} giving 'M'agnitudes, 'S'urface brightness or 'R'adius and {app_min:, app_max: abs_min: abs_max:} giving the limits. The array must be ordered so that vars[0] corresponds to absmag and vars[1] corresponds to sb.

zmin in optional

Minimum redshift. Use with vars when performing Choloniewski fit with phi = 0 and phierr = 0 in some bins

zmax in optional

Maximum redshift. Use with vars when performing Choloniewski fit with phi = 0 and phierr = 0 in some bins

area in optional

Area in square degrees. Use with vars when performing Choloniewski fit with phi = 0 and phierr = 0 in some bins

q0 in optional

Evolution parameters: for ajs_vmax

q1 in optional

Evolution parameters: for ajs_vmax

cov_mat out optional

Covariance matrix

_REF_EXTRA in optional

Extra keywords to pass to mpfit2dfun

Author information

History:

17 Jan 2008 Created, Anthony Smith

4 Apr 2008 Added xrange and yrange

Other attributes

Uses:

mpfit2d (idlutils) ajs_choloniewski

topajs_choloniewski_fit_test

ajs_choloniewski_fit_test

Test ajs_choloniewski_fit

File attributes

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