Synopsis
Return the data used to create the radial or elliptical profile of the source model (after any PSF convolution) for imaging data.
Syntax
get_model_prof( [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, recalc=True] )
Description
The get_model_prof() command calculates the radial - or elliptical - profile of the source model of an imaging dataset and returns an object containing the values. The profile is defined by the existing model compenents, although it is possible to over-ride these values.
The model values used to calculate the profile are those after convolution by any PSF model. Please use the get_source_prof() command if you want to plot the un-convolved model values.
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 - except that the overplot and label arguments are not supported - and are described in its ahelp page.
Fields returned by get_model_prof
The following table lists the fields in the object returned by the get_model_prof() command. These values can also be seen by using the print command on the object:
sherpa> d = get_model_prof(rstep=10) sherpa> print (d)
| Field name | Description |
|---|---|
| xlo | The lower-limits of the radial bins |
| xhi | The upper-limits of the radial bins |
| y | The calculated profile value (sum of model / area) |
| xlabel | The label for the X axis |
| ylabel | The label for the Y axis |
| title | The plot title |
| histo_prefs | The preferences for the plot |
Examples
Example 1
sherpa> d10 = get_model_prof(group_snr=10) sherpa> d20 = get_model_prof(group_snr=20)
Two profiles of the same model are calculated; d10 uses a signal-to-noise ratio of 10 and d20 a value of 20.
Example 2
sherpa> prof_model("src1", group_counts=100)
sherpa> da = get_model_prof(recalc=False)
sherpa> prof_model("src2", rlo=da.xlo, rhi=da.xhi, overplot=True)A profile is calculated and displayed for the dataset "src1". The bin edges from this are then used to create a profile for the dataset "src2", which is overlain on the first plot. Note that recalc is set to False in the call to get_model_prof() to avoid unnescessary re-calculations.
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_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_resid, 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