The task sxirmf creates an SXI RMF file that may be used for spectral analysis with an input spectrum file. The output RMF contains a matrix that, for a given incoming X-ray photon energy, describes the normalized probability of detecting that photon in any PI (pulse invariant) channel. This photon can be redistributed to PI channels away from the expected energy due to interactions with different materials in the detector. For SXI, the response is composed of: (1) a primary Gaussian component, centered at the expected energy, for photons that interact and are detected completely in the sensitive layer of the detector; (2) a secondary Gaussian sub-peak at slightly lower energy, for photons that lose some charge to pixels below the split threshold; (3) a Gaussian silicon escape peak centered at 1.74 keV below the primary peak, for photons that have some energy absorbed by a CCD Si atom; (4) a Gaussian silicon fluorescence peak at constant energy 1.74 keV, if the Si atom excited in the interaction with the X-ray photon emits a Si K alpha photon that is then detected in a different region; and (5) a constant continuum component extending from the peak energy to 0 eV, produced when a photon interacts near the insensitive detector layer and some charge is lost.
The SXI response has been measured on the ground using sources of monochromatic X-rays at a number of energies. However, increasing charge transfer inefficiency (CTI) due to radiation damage changes the response, especially the width of the Gaussian features. Therefore, the RMF depends on both observation time and position on the detector. This information is extracted from the header keywords and the weighted map (WMAP) image, respectively, in the input spectrum file. The WMAP is a small counts image of the extracted events, which is used to weight the response parameters in each region before combining to produce a weighted response. Response parameters are defined in the calibration database (CalDB) at a number of energies, and interpolation is performed to place them on the energy grid in the MATRIX extension before the response is calculated.
The output RMF contains a MATRIX extension, which contains the response matrix, and an EBOUNDS extension listing the output energy boundaries of each PI channel.
The energy bin widths of the input energy grid of the RMF may be set with the 'dein' parameter, with 'eminin' and 'nchanin' determining the energy range for those bins (see below for more detail). The recommended energy bin widths are 2 eV and 24 eV for the energies below and above 12 keV, respectively. The default parameter values of 'eminin=200', 'dein=2,24', and 'nchanin=5900,500' are the recommended values. Furthermore, the output energy grid, EBOUNDS, must match the PI energy grid of the spectrum file. Therefore the default parameter values of 'deout=6', 'eminout=0', and 'nchanout=4096' should be used, unless the user wants to change them for their specific purpose. Changing any of the hidden parameter values for this task is not recommended and must be done with care.
sxirmf infile=spectrum.pha outfile=response.rmf
sxirmf infile=spectrum.pha outfile=response.rmf dein=2 nchanin=11900
June 2, 2023