next up previous contents
Next: Fits and Confidence Intervals Up: Spectral Fitting and XSPEC Previous: The Basics of Spectral

The XSPEC implementation

To summarize the preceding section, the main components of spectral fitting are as follows:

These components are used in the following manner:

This section describes how XSPEC performs these tasks.

C(I): The Observed Spectrum

To obtain each observed spectrum, $C(I)$, XSPEC uses two files: the data (spectrum) file, containing $D(I)$, and the background file, containing $B(I)$. The data file tells XSPEC how many total photon counts were detected by the instrument in a given channel. XSPEC then uses the background file to derive the set of background-subtracted spectra C(I) in units of counts per second. The background-subtracted count rate is given by, for each spectrum:

\begin{displaymath}
C(I) = {D(I)\over{a_{D(I)}t_D}} - {b_{D(I)}\over{b_{B(I)}}} {B(I)\over{a_{b(I)}t_B}}
\end{displaymath} (2.3)

where $D(I)$ and $B(I)$ are the counts in the data and background files; $t_D$ and $t_B$ are the exposure times in the data and background files; $b_{D(I)}$ and $b_{B(I)}$, $a_{D(I)}$ and $a_{B(I)}$ are the background and area scaling values from the spectrum and background respectively, which together refer the background flux to the same area as the observation as necessary. When this is done, XSPEC has an observed spectrum to which the model spectrum can be fit.

R(I,E): The Instrumental Response

Before XSPEC can take a set of parameter values and predict the spectrum that would be detected by a given instrument, XSPEC must know the specific characteristics of the instrument. This information is known as the detector response. Recall that for each spectrum the response $R(I,E)$ is proportional to the probability that an incoming photon of energy $E$ will be detected in channel $I$. As such, the response is a continuous function of $E$. This continuous function is converted to a discrete function by the creator of a response matrix who defines the energy ranges $E_j$ such that:


\begin{displaymath}
R_D(I,J) = {\int_{E_{J-1}}^{E_J} R(I,E) dE\over{E_J - E_{J-1}}}
\end{displaymath} (2.4)

XSPEC reads both the energy ranges, $E_J$, and the response matrix $R_D(I,J)$ from a response file in a compressed format that only stores non-zero elements. XSPEC also includes an option to use an auxiliary response file, which contains an array $A_D(J)$ that is multiplied into $R_D(I,J)$ as follows:

\begin{displaymath}
R_D(I,J) \rightarrow R_D(I,J) \times A_D(J)
\end{displaymath} (2.5)

This array is designed to represent the efficiency of the detector with the response file representing a normalized Redistribution Matrix Function, or RMF. Conventionally, the response is in units of cm2.

M(E): The Model Spectrum

The model spectrum, $M(E)$, is calculated within XSPEC using the energy ranges defined by the response file :


\begin{displaymath}
M_D(J) = \int_{E_{J-1}}^{E_J} M(E) dE
\end{displaymath} (2.6)

and is in units of photons/cm$^2$/s. XSPEC allows the construction of composite models consisting of additive components representing X-ray sources (e.g., power-laws, blackbodies, and so forth), multiplicative components, which modify additive components by an energy-dependent factor (e.g., photoelectric absorption, edges, ...). Convolution and mixing models can then perform sophisticated operations on the result. Models are defined in algebraic notation.

For example, the following expression:

  phabs (power + phabs (bbody))
defines an absorbed blackbody, phabs(bbody), added to a power-law, power. The result then is modified by another absorption component, phabs. For a list of available models, see Chapter 6.


next up previous contents
Next: Fits and Confidence Intervals Up: Spectral Fitting and XSPEC Previous: The Basics of Spectral