Anthony Smith's IDL routines

Galaxies etc.

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

gals/

ajs_lf_plot.pro

This procedure produces a plot of the luminosity function.

Routines

result = ajs_lf_plot_xrange(bincentres=bincentres, loglum=loglum, schechter=schechter, double_schechter=double_schechter)

Choose and set xrange

result = ajs_lf_plot_yrange(phi=phi, input_logphi=input_logphi, plot_logphi=plot_logphi, schechter=schechter, double_schechter=double_schechter)

Choose and set yrange

ajs_lf_plot_legend, legend_text=legend_text, legend_pos=legend_pos, arg_legend_pos=arg_legend_pos, schechter=schechter, double_schechter=double_schechter, color=color, psym=psym, _REF_EXTRA=_REF_EXTRA

Add legend text to the plot

ajs_lf_plot [, bincentres=bincentres] [, phi=phi] [, err_phi=err_phi] [, neg_err_phi=neg_err_phi] [, pos_err_phi=pos_err_phi] [, lower_limit_where=lower_limit_where] [, schechter=schechter] [, sch_range=sch_range] [, double_schechter=double_schechter] [, double_cov_mat=double_cov_mat] [, abs_mag=abs_mag] [, weights=weights] [, jackknife=jackknife] [, nbins=nbins] [, datfile=datfile] [, overplot=overplot] [, plotfile=plotfile] [, show_plot=show_plot] [, open=open] [, close=close] [, xrange=xrange] [, yrange=yrange] [, input_logphi=input_logphi] [, plot_logphi=plot_logphi] [, loglum=loglum] [, ylog=ylog] [, psym=psym] [, symcat_thick=symcat_thick] [, band_name=band_name] [, color=color] [, xstyle=xstyle] [, ystyle=ystyle] [, title=title] [, xtitle=xtitle] [, ytitle=ytitle] [, legend_pos=fltarr(2)] [, legend_text=legend_text] [, cov_mat=cov_mat] [, _REF_EXTRA=_REF_EXTRA]

This procedure produces a plot of the luminosity function.

Routine details

topajs_lf_plot_xrange

result = ajs_lf_plot_xrange(bincentres=bincentres, loglum=loglum, schechter=schechter, double_schechter=double_schechter)

Choose and set xrange

Keywords

bincentres
loglum
schechter
double_schechter

topajs_lf_plot_yrange

result = ajs_lf_plot_yrange(phi=phi, input_logphi=input_logphi, plot_logphi=plot_logphi, schechter=schechter, double_schechter=double_schechter)

Choose and set yrange

Keywords

phi
input_logphi
plot_logphi
schechter
double_schechter

topajs_lf_plot_legend

ajs_lf_plot_legend, legend_text=legend_text, legend_pos=legend_pos, arg_legend_pos=arg_legend_pos, schechter=schechter, double_schechter=double_schechter, color=color, psym=psym, _REF_EXTRA=_REF_EXTRA

Add legend text to the plot

Keywords

legend_text
legend_pos
arg_legend_pos
schechter
double_schechter
color
psym
_REF_EXTRA

topajs_lf_plot

ajs_lf_plot [, bincentres=bincentres] [, phi=phi] [, err_phi=err_phi] [, neg_err_phi=neg_err_phi] [, pos_err_phi=pos_err_phi] [, lower_limit_where=lower_limit_where] [, schechter=schechter] [, sch_range=sch_range] [, double_schechter=double_schechter] [, double_cov_mat=double_cov_mat] [, abs_mag=abs_mag] [, weights=weights] [, jackknife=jackknife] [, nbins=nbins] [, datfile=datfile] [, overplot=overplot] [, plotfile=plotfile] [, show_plot=show_plot] [, open=open] [, close=close] [, xrange=xrange] [, yrange=yrange] [, input_logphi=input_logphi] [, plot_logphi=plot_logphi] [, loglum=loglum] [, ylog=ylog] [, psym=psym] [, symcat_thick=symcat_thick] [, band_name=band_name] [, color=color] [, xstyle=xstyle] [, ystyle=ystyle] [, title=title] [, xtitle=xtitle] [, ytitle=ytitle] [, legend_pos=fltarr(2)] [, legend_text=legend_text] [, cov_mat=cov_mat] [, _REF_EXTRA=_REF_EXTRA]

This procedure produces a plot of the luminosity function.

Data may be entered in four forms:

(1) Values at particular absolute magnitudes (bincentres, phi[, err_phi])

(2) (double) Schechter function parameters (always plotted if supplied)

(3) Array of galaxy absolute magnitudes and weights (e.g., 1/Vmax) - ignored if (1) set

(4) From a file, to be read using ajs_lf_read - ignored if (1) or (3) set

In addition, a (double) Schechter function fit can be performed.

Keywords

bincentres in optional
phi in optional
err_phi in optional
neg_err_phi in optional

Negative error, for asymmetric

pos_err_phi in optional

Positive error, for asymmetric

lower_limit_where in optional

Indices of values of phi where the upper error bar is a lower limit (i.e., plot arrowhead and empty circle, if psym=16 set)

schechter in optional

Plot Schechter function. Set to undefined variable to calculate best-fitting Schechter function and return the parameters. Or set /schechter to calculate and print the parameters

sch_range in optional

Range of absolute magnitudes for (double) Schechter fit (ajs_schechter_fit)

double_schechter in out optional

Plot double Schechter function. Set to undefined variable to calculate best-fitting double Schechter function and return the parameters. Or set /double_schechter to calculate and print the parameters

double_cov_mat out optional

Covariance matrix of double Schechter function fit

abs_mag in optional

Array of absolute magnitudes (use with weights)

weights in optional

Array of weights (e.g., 1/Vmax) for estimating luminosity function

jackknife in optional

Used with absmag, this is an integer for each galaxy giving the jackknife region in which the galaxy lies

nbins in optional

Number of bins for abs_mag & weights (use with xrange or bincentres)

datfile in optional

Input data: binary format file created by ajs_lf_write

overplot in optional

Set /overplot to overplot

plotfile in optional

EPS file name

show_plot in optional

Set /show_plot to open the EPS file for viewing

open in optional

Set /open to open the EPS file and leave it for overplotting

close in optional

Set /close to close the EPS plotfile

xrange in out optional
yrange in out optional
input_logphi in optional

Set /input_logphi to indicate data values are log(phi) rather than phi

plot_logphi in optional

Set /plot_logphi to plot log(phi) rather than phi (redundant if overplotting)

loglum in optional

Set /loglum for log luminosities rather than magnitudes

ylog in optional default=1-plot_logphi
psym in optional

Uses D. Fanning's symcat (e.g., 16 = filled circle) for non-standard psym values

symcat_thick in optional

Thick keyword for symcat

band_name in optional

For the x-axis: M_ - 5 log h

color in optional
xstyle in optional

Default /xstyle

ystyle in optional

Default /ystyle

title in optional

Title for plot

xtitle in optional

Title for x-axis

ytitle in optional

Title for y-axis

legend_pos in out optional type=fltarr(2)

[xpos, ypos] for legend text. Output is position for next legend label, for overplotting. varies from [0, 0] = [left, bottom] to [1, 1] = [right, top].

legend_text in optional

Legend text

cov_mat out optional

Covariance matrix of Schechter function fit

_REF_EXTRA in optional

Extra keywords for plots

Examples

bincentres = [-20, -22, -24]

phi = [1e-2, 1e-2, 1e-3]

err_phi = [1e-2, 1e-3, 1e-4]

ajs_lf_plot, bincentres=bincentres, phi=phi

ajs_lf_plot, bincentres=bincentres, phi=phi, err_phi=err_phi

ajs_lf_plot, bincentres=bincentres, phi=phi, err_phi=err_phi, psym=16

ajs_lf_plot, schechter=[-24,-1,1e-2]

ajs_lf_plot, schechter=[-24,-1,1e-2], color=2, bincentres=bincentres, phi=phi, err_phi=err_phi, xrange=[-16,-26], psym=16, plotfile='~/lf.eps', /show_plot

Author information

History:

13 Feb 2008 Created, Anthony Smith

5 Mar 2008 err_phi no longer required

6 Mar 2008 Added abs_mag, weights and nbins

18 Mar 2008 Added jackknife

25 Mar 2008 Added legend_text and legend_pos

4 Apr 2008 Added sch_range (ajs_schechter_fit range)

7 Apr 2008 Schechter function info included in call to ajs_vmax_lf

21 Apr 2008 double Schechter function input

23 Apr 2008 psym now uses D. Fanning's symcat routine; added lower_limit_where keyword

1 May 2008 Uses ajs_plot_pos for legend position

File attributes

Modifcation date: Wed Oct 29 11:00:51 2008
Lines: 715
Docformat: rst rst