xsclumin

Context: models

Synopsis

The XSPEC clumin convolution model: calculate luminosity

Syntax

xsclumin

The xsclumin model is a convolution model component.

Description

The model is described at [1] .

Examples

Example 1

>>> create_model_component("xsclumin", "mdl")
>>> print(mdl)

Create a component of the xsclumin model and display its default parameters. The output is:

mdl
   Param        Type          Value          Min          Max      Units
   -----        ----          -----          ---          ---      -----
   mdl.Emin     frozen          0.5            0        1e+06        keV
   mdl.Emax     frozen           10            0        1e+06        keV
   mdl.Redshift frozen            0       -0.999           10           
   mdl.lg10Lum  thawed           40         -100          100        cgs

Example 2

With the following definitions:

>>> xsclumin.clumin
>>> xsphabs.absmdl
>>> xspowerlaw.plmdl
>>> xsgaussian.gmdl
>>> srcmdl = plmdl + gmdl

then the model can be applied in a number of ways, such as:

>>> mdl1 = clumin(absmdl * srcmdl)
>>> mdl2 = absmdl * clumin(srcmdl)
>>> mdl3 = absmdl * (plmdl + clumin(gmdl))

ATTRIBUTES

The attributes for this object are:

Attribute	Definition
Emin	Minimum energy over which the luminosity is calculated.
Emax	Maximum energy over which the luminosity is calculated.
Redshift	redshift of the source
lg10Lum	log (base 10) of the luminosity in erg/s

Notes

Unlike XSPEC, the convolution model is applied directly to the model, or models, rather than using the multiplication symbol.

See [1] for the meaning and restrictions, in particular the necessity of freezing the amplitude, or normalization, of the emission component (or components) at 1.

References

[1]

Changes in CIAO

Changed in CIAO 4.14

The default lg10Lum parameter value has changed from -40 to 40 to match XSPEC.

Added in CIAO 4.13

XSPEC version

CIAO 4.18 comes with support for version 12.14.0k of the XSPEC models. This can be checked with the following:

% python -c 'from sherpa.astro import xspec;
print(xspec.get_xsversion())'
12.14.0k

Bugs

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