Synopsis
Plot a radial or elliptical profile of the residuals (data - model) to imaging data.
Syntax
prof_resid( [id], [model=None, rstep=None, rmin=None, rmax=None, rlo=None, rhi=None, xpos=None, ypos=None, ellip=None, theta=None, group_counts=None, group_snr=None, label=True, recalc=True, overplot=False, clearwindow=True] )
Description
The prof_resid command calculates the radial - or elliptical - profile of the residuals to imaging data (data - model) and plots it. The profile is defined by the existing model compenents, although it is possible to over-ride these values.
The residual values are calculated as
data - model
The prof_delchi() command can be used if you want to plot
(data - model) / error
Loading the routine
The routine can be loaded into Sherpa by saying:
from sherpa_contrib.profiles import *
Argument options
The argument options are the same as for the prof_data() command, and are described in its ahelp page.
Changing the plot defaults
The get_resid_prof_prefs() returns the current plot preferences used by prof_resid(). Changing these settings will therefore change the appearance of any new plots created by prof_resid(). For example
sherpa> get_resid_prof_prefs()["xlog"] = True
will cause any new residual plots to use logarithmic scaling for the X axis. A full list of the preferences can be found by saying
unix% ahelp get_resid_prof_prefs
Examples
Example 1
sherpa> prof_resid() ... sherpa> prefs = get_resid_prof_prefs() sherpa> prefs["xlog"] = True sherpa> prof_resid()
The preferences are set so that the x axis should be drawn using log scaling. Setting the get_resid_prof_prefs values only affects new plots made after the setting was changed.
Example 2
sherpa> prof_resid(group_snr=15)
The residuals are plotted after the bins have been grouped so that each bin has a signal to noise ratio of 15 or more, where the calculation is done using the data and not the residual values.
Example 3
sherpa> prof_resid(rlo=[0,10,20,40,60,100,200])
Calculate the residual profile in radial bins with low/high limits of 0-10, 10-20, 20-40, 40-60, 60-100, 100-200.
Bugs
See the bugs pages on the Sherpa website for an up-to-date listing of known bugs.
See Also
- contrib
- get_data_prof, get_data_prof_prefs, get_delchi_prof, get_delchi_prof_prefs, get_fit_prof, get_model_prof, get_model_prof_prefs, get_resid_prof, get_resid_prof_prefs, get_source_prof, get_source_prof_prefs, plot_chart_spectrum, plot_marx_spectrum, prof_data, prof_delchi, prof_fit, prof_fit_delchi, prof_fit_resid, prof_model, prof_source, sherpa_profiles
- data
- get_arf_plot, get_bkg_chisqr_plot, get_bkg_delchi_plot, get_bkg_fit_plot, get_bkg_model_plot, get_bkg_plot, get_bkg_ratio_plot, get_bkg_resid_plot, get_bkg_source_plot
- modeling
- normal_sample, t_sample, uniform_sample
- plotting
- get_cdf_plot, get_energy_flux_hist, get_pdf_plot, get_photon_flux_hist, get_pvalue_plot, get_pvalue_results, get_split_plot, plot, plot_arf, plot_bkg, plot_bkg_chisqr, plot_bkg_delchi, plot_bkg_fit, plot_bkg_fit_delchi, plot_bkg_fit_resid, plot_bkg_model, plot_bkg_ratio, plot_bkg_resid, plot_bkg_source, plot_cdf, plot_chisqr, plot_data, plot_delchi, plot_energy_flux, plot_fit, plot_fit_delchi, plot_fit_resid, plot_model, plot_model_component, plot_order, plot_pdf, plot_photon_flux, plot_pvalue, plot_ratio, plot_resid, plot_scatter, plot_source, plot_source_component, plot_trace, set_xlinear, set_xlog, set_ylinear, set_ylog
- psfs
- plot_kernel
- statistics
- get_chisqr_plot, get_delchi_plot
- visualization
- contour_resid