calc_model_sum

sherpa.astro.ui.calc_model_sum(lo=None, hi=None, id=None, bkg_id=None)

Sum up the fitted model over a pass band.

Sum up M(E) over a range of bins, where M(E) is the per-bin model value after it has been convolved with any instrumental response (e.g. RMF and ARF or PSF). This is intended for one-dimensional data sets: use calc_model_sum2d for two-dimensional data sets. The calc_source_sum function is used to calculate the sum of the model before any instrumental response is applied.

Parameters
  • lo (number, optional) – If both are None or both are set then sum up over the given band. If only one is set then use the model value in the selected bin. The units for lo and hi are given by the current analysis setting.

  • hi (number, optional) – If both are None or both are set then sum up over the given band. If only one is set then use the model value in the selected bin. The units for lo and hi are given by the current analysis setting.

  • id (int or str, optional) – Use the source expression associated with this data set. If not given then the default identifier is used, as returned by get_default_id.

  • bkg_id (int or str, optional) – If set, use the model associated with the given background component rather than the source model.

Returns

signal – The model value (sum or individual bin).

Return type

number

See also

calc_data_sum

Sum up the observed counts over a pass band.

calc_energy_flux

Integrate the unconvolved source model over a pass band.

calc_photon_flux

Integrate the unconvolved source model over a pass band.

calc_source_sum

Sum up the source model over a pass band.

set_model

Set the source model expression for a data set.

Notes

The units of lo and hi are determined by the analysis setting for the data set (e.g. get_analysis). The summation occurs over those points in the data set that lie within this range, not the range itself.

Any existing filter on the data set - e.g. as created by ignore or notice - is ignored by this function.

The units of the answer depend on the model components used in the source expression and the axis or axes of the data set.

Examples

Calculate the model evaluated over the full data set (all points or pixels of the independent axis) for the default data set, and compare it to the sum for th first background component:

>>> tsrc = calc_model_sum()
>>> tbkg = calc_model_sum(bkg_id=1)

Sum up the model over the data range 0.5 to 2 for the default data set, and compared to the data over the same range:

>>> calc_model_sum(0.5, 2)
404.97796489631639
>>> calc_data_sum(0.5, 2)
745.0

Calculate the model sum, evaluated over the range 20 to 22 Angstroms, for the first background component of the “histate” data set:

>>> set_analysis("histate", "wavelength")
>>> calc_model_sum(20, 22, "histate", bkg_id=1)

In the following example, a small data set is created, covering the axis range of -5 to 5, and an off-center gaussian model created (centered at 1). The model is evaluated over the full data grid and then a subset of pixels. As the summation is done over those points in the data set that lie within the requested range, the sum for lo=-2 to hi=1 is the same as that for lo=-1.5 to hi=1.5:

>>> load_arrays('test', [-5, -2.5, 0, 2.5, 5], [2, 5, 12, 7, 3])
>>> set_source('test', gauss1d.gmdl)
>>> gmdl.pos = 1
>>> gmdl.fwhm = 2.4
>>> gmdl.ampl = 10
>>> calc_model_sum(id='test')
9.597121089731253
>>> calc_model_sum(-2, 1, id='test')
6.179472329646446
>>> calc_model_sum(-1.5, 1.5, id='test')
6.179472329646446