Last modified: December 2023

Jump to: Description · Examples · PARAMETERS · Notes · Bugs · See Also

AHELP for CIAO 4.16 Sherpa


Context: modeling


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


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


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.


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')


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. This is currently unused.
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.


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.


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

See Also

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
get_default_id, list_bkg_ids, list_data_ids, list_model_ids, show_bkg_model, show_bkg_source
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
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
delete_psf, load_conv
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
calc_chisqr, calc_energy_flux, calc_model_sum, calc_photon_flux, calc_source_sum, calc_stat, eqwidth
contour_model, contour_ratio, contour_resid