Last modified: December 2020

URL: https://cxc.cfa.harvard.edu/sherpa/ahelp/load_template_model.html
Jump to: Description · Examples · PARAMETERS · Notes · Bugs · See Also


AHELP for CIAO 4.13 Sherpa v1

load_template_model

Context: modeling

Synopsis

Load a set of templates and use it as a model component.

Syntax

load_template_model(modelname, templatefile, dstype=Data1D, sep=' ',
comment='#', method=linear_interp, template_interpolator_name='default')

modelname - str
templatefile - str
dstype - data class to use, optional
sep - str, optional
comment - str, optional
method - func
template_interpolator_name - str

Description

A template model can be considered to be an extension of the table model supported by `load_table_model` . In the template case, a set of models (the "templates") are read in and then compared to the data, with the best-fit being used to return a set of parameters.


Examples

Example 1

Load in the templates from the file "index.tmpl" as the model component "kerr", and set it as the source model for the default data set. The optimisation method is switched to use a grid search for the parameters of this model.

>>> load_template_model("kerr", "index.tmpl")
>>> set_source(kerr)
>>> set_method('gridsearch')
>>> set_method_opt('sequence', kerr.parvals)
>>> fit()

Example 2

Fit a constant plus the templates, using the neville scheme for integrating the template onto the data grid. The Monte-Carlo based optimiser is used.

>>> load_template_model('tbl', 'table.lis',
...                     sherpa.utils.neville)
>>> set_source(tbl + const1d.bgnd)
>>> set_method('moncar')

PARAMETERS

The parameters for this function are:

Parameter Definition
modelname The identifier for this table model.
templatefile The name of the file to read in. This file lists the template data files.
dstype What type of data is to be used. Supported values include `Data1D` (the default) and `Data1DInt` .
sep The separator character. The default is ' ' .
comment The comment character. The default is '#' .
method The interpolation method to use to map the input data onto the coordinate grid of the data set. Linear, nearest-neighbor, and polynomial schemes are provided in the sherpa.utils module.
template_interpolator_name The method used to interpolate within the set of templates. The default is default . A value of none turns off the interpolation; in this case the grid-search optimiser must be used to fit the data.

Notes

Examples of interpolation schemes provided by `sherpa.utils` are: `linear_interp` , `nearest_interp` , and `neville` .

The template index file is the argument to `load_template_model` , and is used to list the data files. It is an ASCII file with one line per template, and each line containing the model parameters (numeric values), followed by the MODELFLAG column and then the file name for the data file (its name must begin with FILE). The MODELFLAG column is used to indicate whether a file should be used or not; a value of 1 means that the file should be used, and a value of 0 that the line should be ignored. The parameter names are set by the column names.

The data file - the last column of the template index file - is read in and the first two columns used to set up the x and y values ( `Data1D` ) or xlo, xhi, and y values ( `Data1DInt` ). These files must be in ASCII format.

The method parameter determines how the template data values are interpolated onto the source data grid.

The template_interpolator_name parameter determines how the dependent axis (Y) values are interpolated when the parameter values are varied. This interpolation can be turned off by using a value of none , in which case the grid-search optimiser must be used. See `load_template_interpolator` for how to create a valid interpolator. The "default" interpolator uses `sherpa.models.KNNInterpolator` with k=2 and order=2.


Bugs

See the bugs pages on the Sherpa website for an up-to-date listing of known bugs.

See Also

data
dataspace1d, dataspace2d, datastack, fake, 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, load_arf, load_arrays, load_ascii, load_bkg, load_bkg_arf, load_bkg_rmf, load_data, load_grouping, load_image, load_multi_arfs, load_multi_rmfs, load_pha, load_quality, load_rmf, load_staterror, load_syserror, load_table, pack_image, pack_pha, pack_table, unpack_arf, unpack_arrays, unpack_ascii, unpack_bkg, unpack_data, unpack_image, unpack_pha, unpack_rmf, unpack_table
filtering
load_filter
info
get_default_id, list_bkg_ids, list_data_ids, list_model_ids, show_bkg_model, show_bkg_source
modeling
add_model, add_user_pars, clean, create_model_component, delete_bkg_model, delete_model, delete_model_component, get_model, get_model_autoassign_func, get_model_component, get_model_component_image, get_model_component_plot, get_model_plot, get_num_par, get_num_par_frozen, get_num_par_thawed, get_order_plot, get_par, get_pileup_model, get_source, get_source_component_image, get_source_component_plot, get_source_contour, get_source_image, get_source_plot, get_xsabund, get_xscosmo, get_xsxsect, get_xsxset, image_model, image_model_component, image_source, image_source_component, integrate, link, load_table_model, load_template_interpolator, load_user_model, normal_sample, reset, save_model, save_source, set_bkg_model, set_bkg_source, set_full_model, set_model, set_model_autoassign_func, set_pileup_model, set_source, set_xsabund, set_xscosmo, set_xsxsect, set_xsxset, t_sample, uniform_sample
plotting
get_cdf_plot, get_pdf_plot, get_pvalue_plot, get_pvalue_results, plot_cdf, plot_model, plot_model_component, plot_pdf, plot_pvalue, plot_scatter, plot_source, plot_source_component, plot_trace
psfs
delete_psf, load_conv
saving
save_arrays, save_data, save_delchi, save_error, save_filter, save_grouping, save_image, save_pha, save_quality, save_resid, save_staterror, save_syserror, save_table
statistics
load_user_stat
utilities
calc_chisqr, calc_energy_flux, calc_model_sum, calc_photon_flux, calc_source_sum, calc_stat, eqwidth
visualization
contour_model, contour_ratio, contour_resid