Let's consider the spectral case first.
Let be and
the number of photon counts for each
spectral bin in the source+background and the background spectra
created via the Projection tasks. The two spectra are assumed to have
the same bin structure. The background subtracted spectrum
is then defined as follows:
where and
are the areas of the sky portions from which the
source+background and the background spectra were obtained. The exposure
time for the source+background spectra is assumed by default to be the same
as the exposure time for the background one.
The corresponding statistical error is defined as
This formula is modified consistently for data which have been already corrected in the Projection run in correction mode.
When data were not corrected within the Projection run,
it might be necessary to keep into account the following:
in a pointed observation different sky regions affect
systematically different detector portions. Therefore
when the source+background and the background spectra
are taken from sky regions which are relatively far from each other,
it may be necessary to keep into account the differences between the
corresponding spectral correction vectors. Let be
and
the effective areas
corresponding to each amplitude channel (ampl) of the two spectra (they
are obtained from the corresponding correction vectors, which are
effective areas as a function of energy). In this specific case
the background subtraction is performed as
All this is not applicable -- because not necessary -- whenever Projection was run in correction mode, since all spectra are processed in such a way that their correction vectors are all identical to the effective area on-axis.
As a last step the result spectra can be optionally rebinned,
either interactively or automatically. The automatic rebin is
performed in order to obtain for each of the new bins
a signal to noise ratio not less than a specified value.
Once the signal to noise ratio is fixed to a value
,
the new binsize of a spectral bin starting from a given amplitude
channel
is defined as
where is a multiple of the original spectra binsize
(i.e. the binsize declared in the Projection command
BIN AMPL which was used in creating the spectra).
The sequence of binsizes in the rebinned spectrum is defined
by the recursive
The results of the spectral
background subtraction are always converted
into counts s bin
.
Let's now briefly consider the light curve background subtraction case. Let be
and
the two sequences of time bins, and
,
the corresponding (optional) correction vectors
(they are not produced when Projection is run in correction mode).
The two light curves are assumed to
have the same structure, i.e. the same binsizes, the same gaps and the
same extension in time. The corrected and background subtracted light
curve is then given by
The results of the light curve
background subtraction are always converted
into counts s .
All the operations described above are automatically performed either via the Correction commands or (partly) within Projection run in correction mode (applicable only for PSPC data). In the following a detailed description of the correction commands syntax is given: