** Next:** Data Analysis Concepts for the Next Generation Space Telescope ** Up:** Data Analysis Applications** Previous:** News on the ISOPHOT Interactive Analysis PIA

Table of Contents -- Index -- PS reprint -- PDF reprint

ASP Conference Series, Vol. 145, 1998

Editors: R. Albrecht, R. N. Hook and H. A. Bushouse

**A. Alexov, W. McLaughlin and D. Huenemoerder**

AXAF Science Center; Smithsonian Astrophysical Observatories,
TRW, and MIT, Cambridge, MA 02138

Initial data processing applies transformations from detector to sky coordinates for each photon in the data set. However, with grating data, additional event processing must be performed before data analysis can commence.

In order to find where the dispersed photons lie, the center of the zero order source position on the sky must be determined, since this is the origin of the spectral coordinates. Source detection must not be fooled by emission lines located away from the center of the field; several methods exist to discriminate bright emission lines from zero order sources. Primarily, emission lines have a small variance in the PHA (energy) spectrum, while zero order has a large variance since it encompasses all of the source energies. Alternatively, the PSF (Point Spread Function) can be used instead of PHA, to weed out emission lines. Once zero order source positions are found, the target source is identified by matching the sky positions with an observation target list.

Grating data events may be categorized by their part of the spectrum. This is done by creating mask regions for each grating part and by checking every photon for inclusion within the mask regions. For the HETG grating data, the relevant parts for each source are: zero order, MEG photons, HEG photons, and background photons. Regions are defined as rectangles in diffraction coordinates, in which the zero order is the origin, and one axis is parallel to the spectrum (dispersion direction), and the other perpendicular to it (cross-dispersion direction). The width of the rectangle is calculated using the effective PSF (mirror psf, Rowland geometry astigmatism, instrument de-focus, and aspect). These regions are translated into sky coordinates using each zero order source position, and by rotating the region by the grating angle (known) plus the mean observed telescope roll angle. Any event inside the rectangle is then assigned to that part of the spectrum (HEG or MEG, or LEG). Zero order photons are assigned by being within some radius of the zero order centroid. Each photon is tagged with one or more source ID's as well as grating part(s), depending on the number of regions into which it falls (overlaps).

To calculate the diffraction angle (*r*) of a photon, work needs to be
done in the
geometric system in which the photon was diffracted. The sky coordinates,
which are referenced to the telescope mirror node and therefore independent
of the grating reference node, are only useful for the imaging of zero
order and the filtering of grating arm photons.
Grating diffraction coordinates are
referenced to the grating assembly node. Reverting back to the original
chip coordinates (chip number, chip position, grating node and zero order
position) allows grating diffraction coordinates to be calculated for each
time interval.

Now that diffraction angles have been determined parallel to the
dispersion angle, the basic grating equation, ,can be applied to determine for each photon.

Here,

The ACIS detector provides moderate spectral resolution via a
pulse-height (PHA, or PI) for every photon. This resolution is enough
to determine, with high confidence, what the spectral order is at any r
or , since only integral multiples are allowed at any given position
(grating equation, diagram below). Order sorting is also useful for
photons in overlapping mask regions. These events can be resolved as
belonging to one arm/source versus another through order sorting. Since the
HRC detector has no energy resolution, a method has yet to be
determined for resolving orders or any overlapping regions in data
taken by this instrument.
The grating equation is used to calculate an estimate of the grating
order (), using physical constants and previously calculated values:
*m*_{est} = *Pr*(*PI*)/*hc*.
Here, *hc* is a physical constant; *PI* is the ``Pulse Invariant"
energy from the CCD;
*r* is the calculated diffraction angle; and, *P* is the known grating
period.

The ACIS instrument response determines the distribution of PI.
This information permits the calculation of an
allowed range of *m* for each photon. If is within this
range, then the
grating order is equal to the rounded value of . Otherwise,
the order is unresolved.
For ACIS overlap regions, order is calculated for each region. If the
order can be uniquely identified as belonging to a single region then it
is resolved. Otherwise, the order is left as unresolved.

Now that the data are complete with identifiers for the part of the spectrum, diffraction angles, orders, wavelengths, and counts spectra vs can be created. These spectra are used for further data analysis (i.e., emission line identification and flux analysis).

Processing grating data comprised of multiple sources is a non-trivial task. Many design tradeoffs were considered in attempting to create efficient and effective analysis software. This section identifies some of the considerations taken into account when designing the gratings software.

The possibility that parts of different sources may overlap imposes a great challenge to the data analysis software. All overlapping regions need to be identified and the software must be able to decompose these regions into the component source parts. While standard IRAF region software can easily support overlapping regions of a few sources via region algebra, a mechanism is still necessary for keeping track of the source parts in a given overlap region. As the total number of regions (including overlaps) has an exponential growth (5 to the n for HETG), this method is impractical when dealing with several sources. For instance, the worst case scenario for 10 HETG sources is over 9.5 million source parts and overlap combinations to track. Realistic scenarios may contain approximately 10 sources, but these sources will typically be more spread out (i.e., Orion Trapezium).

To circumvent this problem, the software has been designed to maintain tables for all of the parts of each source (resulting in 5n tables for HETG). Photons are checked for inclusion in each table and source bits are set to indicate the sources and parts to which the photon may belong. Since the geometry of the regions being utilized is limited to circles and rectangles, the tables simply contain sky coordinate range boundaries. To save memory, the axis spanning the minimum distance for any given source part is offset to zero and utilized as the index. For cases where the instrument roll angle equals 45 degrees, the index axis is arbitrarily chosen as neither of the axes provides an advantage.

AXAF data fields may contain multiple sources. In order to be able to detect, mask out, coordinate transform, and order sort all these sources correctly is quite a challenge. We have made the software flexible for the users to be able to specify source specific characteristics/mask widths and response matrices. We hope that these extra features will allow severely crowded fields to be scientifically analyzed in the same context as more common single source fields.

For more details on the AXAF grating flight processing software, see:

http://space.mit.edu/ASC/analysis/L1.5_overview/L15.html

We are grateful to Maureen Conroy for initial design work on the XRCF grating software, which has lead to the flight version. This project is supported by NASA contract NAS8-39073 (ASC).

© Copyright 1998 Astronomical Society of the Pacific, 390 Ashton Avenue, San Francisco, California 94112, USA

** Next:** Data Analysis Concepts for the Next Generation Space Telescope ** Up:** Data Analysis Applications** Previous:** News on the ISOPHOT Interactive Analysis PIA

Table of Contents -- Index -- PS reprint -- PDF reprint