Higherorder Responses for HRCS/LETG Spectra
CIAO 4.16 Science Threads
Overview
Synopsis:
Because of the low energy resolution in the HRCS, the PHA2 file contains only two rows (negative and positive) containing all the spectral orders. While it is not possible to separate the overlapping orders, CIAO can make response files (gRMFs and gARFs) for higher orders for modeling and fitting.
Deciding how many orders to include may not be simple. The best advice is to use the empirical approach, and keep adding one or two orders at a time to your analysis until variations in fit results reach the level of other uncertainties. If nothing else, this will provide an enhanced appreciation for the level of systematic uncertainties in your analysis. This is also a good approach to use because the computational load steadily increases as orders are added.
Purpose:
This thread describes how to estimate the importance of higher orders when analyzing LETG/HRCS spectra, so that one can at least determine the minimum number of orders with which to start. Keep in mind that the spectrum of order m has m times the dispersion of 1st order, i.e. continuum is more spread out in terms of mm/Å on the detector. Likewise, resolved lines will become broader (in mm) as the order increases.
Last Update: 30 Jan 2017  Reviewed for CIAO 4.9. Updated the equations to use mathjax.
Contents
Explanation
The number of counts at the single detector position mλ is equal to the sum of counts from 1st order (at emitted wavelength λ), 2nd order (at λ/2), 3rd order (λ/3), etc. This can be written as
\[ I(m\lambda) = f(\lambda)A_1(\lambda) + f(\lambda/2)A_2(\lambda/2) + f(\lambda/3)A_3(\lambda/3) + \cdots = \sum_m f(\lambda/m)A_m(\lambda/m) \]where I is the intensity (counts/s/Å) observed at a given location on the detector, f is the source flux (photons/s/cm^{2}/Å) at a specific wavelength, and A_{m} is the LETG/HRCS effective area (cm^{2}) for order m at that wavelength.
This can be rewritten for integrated line intensities (in cts/s) as
\[ L(m\lambda) = L(\lambda)+L(\lambda/2)\left[\frac{A_2(\lambda/2)}{A_1(\lambda)}\right] + L(\lambda/3)\left[\frac{A_3(\lambda/3)}{A_1(\lambda)}\right] + \cdots = \sum_m L(\lambda/m)\left[ \frac{A_m(\lambda/m)}{A_1(\lambda)} \right] \]and for continuum (cts/s/Å) as
\[ C(m\lambda) = C(\lambda) + \frac{C(\lambda/2)}{2}\left[\frac{A_2(\lambda/2)}{A_1(\lambda)}\right] + \frac{C(\lambda/3)}{3}\left[\frac{A_3(\lambda/3)}{A_1(\lambda)}\right] + \cdots = \sum_m \frac{C(\lambda/m)}{m}\left[ \frac{A_m(\lambda/m)}{A_1(\lambda)} \right] \](Recall that the extra 1/m factor for the continuum equation is because continuum intensity is measured per unit wavelength, and the dispersion of higherorder spectra is proportional to m. Line intensities are integrated, so broadening doesn't matter for them.)
In the above equations, the only directly measured intensity is L(mλ) or C(mλ), but as long as most of the counts at a given mλ come from 1st order, one can estimate the contribution of each order in the mλ spectrum, as shown in the examples. The effective areas, A_{m}, can be read from these plots.
[PostScript file of this plot]
Note that the two plots contain the same information, but present it differently: mλ vs λ for the xaxis. The shaded area in the top plot marks the plate gaps around mλ = 53Å and +64Å.
At any wavelength, the combined effective area of orders beyond m=10 is at most a few percent of the 1st order EA unless nearly all of the source's emission is between 6 and 17 Å, in which case the observer would almost certainly want to use the HETG with ACISS.
Example Analyses
Continuum
Here we use ObsId 460 (LETG/HRCS, 3C 273) to illustrate the steps. This plot of the positiveorder spectrum (Figure 1) has the wavelengths we will focus on marked with a red "X".
Figure 1: Positiveorder spectrum of ObsID 460

3rd order often provides the strongest higherorder contamination, so we begin with that, looking at mλ = 90 Å, which has contributions from 1st order (90 Å), 2nd order (45 Å), 3rd order (30 Å), etc.

The 3rd order contribution is approximately equal to (1/3)C(λ/3)[A_{3}(λ/3)/A_{1}(λ)]. There are about 4 counts/bin at 90 Å (C(λ) and 60 at 30 Å (C(λ/3). From the effective area plots, A_{1}(90 Å) = 8 cm^{2} and A_{3}(30 Å) = 1.1 cm^{2}, so 3rd order contributes (1/3)(60 counts)(1.1 cm^{2}/8 cm^{2}) = 2.8 of the total of 4 counts/bin.

Similarly, the number of counts/bin from other orders is:
2nd: (1/2)(60 counts at 45.0 Å)(0.8 cm^{2}/8 cm^{2}) = 3.0
4th: (1/4)(70 counts at 22.5 Å)(0.5 cm^{2}/8 cm^{2}) = 1.1
5th: (1/5)(100 counts at 18.0 Å)(0.5 cm^{2}/8 cm^{2}) = 1.2
6th: (1/6)(120 counts at 15.0 Å)(0.7 cm^{2}/8 cm^{2}) = 1.8
7th: (1/7)(120 counts at 12.9 Å)(0.3 cm^{2}/8 cm^{2}) = 0.6
8th: (1/8)(120 counts at 11.2 Å)(0.6 cm^{2}/8 cm^{2}) = 1.1
9th: (1/9)(130 counts at 10.0 Å)(0.35 cm^{2}/8 cm^{2}) = 0.6
10th: (1/10)(150 counts at 9.0 Å)(0.35 cm^{2}/8 cm^{2}) = 0.7
In this case, orders 210 contribute an estimated total of more than 12 counts/bin. The fact that this is higher than the observed total of 4 counts/bin simply means that higher orders are also important at shorter wavelengths, thus inflating the apparent spectral intensity at λ/m and leading to overestimates of higherorder contributions at λ. In this extreme example, higher orders contribute virtually all the flux at long wavelengths (beyond roughly 70 Å), and spectral fitting at those wavelengths will be of dubious quality.
Line Spectrum
For this example we refer to the coronal line spectrum of Capella, as shown in Figure 9.25 of the POG (from Brinkman et al. 2000, ApJ, 530, L111).

Say we want to analyze the region around 45 Å, which may contain 2nd order features from 22 Å, 3rd order from 15 Å, etc. There is a strong line at 15.02 Å so we focus on 3rd order.

The 3rd order contribution is approximately equal to L(15.02 Å)[A_{3}(15.02 Å)/A_{1}(45.06 Å)]. From Brinkman et al., the intensity of the 15.02 Å line is 44.2 cts/ks. From the effective area plots, A_{1}(45 Å) = 25 cm^{2} and A_{3}(15 Å) = 2.5 cm^{2}, so the 3rd order line at mλ should have an intensity of (44.2 cts/ks)(2.5 cm^{2}/25 cm^{2}) = 4.4 cts/ks. Indeed, Brinkman et al. measure 4.2 cts/ks, which is nearly as large as the 1st order intensity of the Si XII line at 44.16 Å.

Similarly, 5th order of the 9.169 Å Mg XI line will appear at 45.84 Å with (6.2 cts/ks)(1.2 cm^{2}/25 cm^{2}) = 0.3 cts/ks, probably blending with the 1.9cts/ksec Si XII line at 45.68 Å. the Si XII line at 44.16 Å.

The Si XIII line at 6.648 Å is the shortest wavelength feature of any significance, and it would appear in 7th order at 46.54 Å with (5.1 cts/ks)(0.35 cm^{2}/25 cm^{2}) = 0.07 cts/ks, which is small but perhaps nonnegligible depending on the analysis goals.
Related Threads
The following threads show how to create response files for HRCS/LETG data:
History
10 May 2005  original version, new for CIAO 3.2 
06 Dec 2005  reviewed for CIAO 3.3: no changes 
02 Nov 2006  added note about orders beyond m=10 in Explanation section 
01 Dec 2006  reviewed for CIAO 3.4: no changes 
11 Jan 2008  reviewed for CIAO 4.0: no changes 
06 Jun 2008  edited sentence in Overview to remove "(m up to 10)"; refer to content for thread for information on m values. 
13 Feb 2009  updated for CIAO 4.1: image is inline 
13 Jan 2010  reviewed for CIAO 4.2: no changes 
12 Jan 2011  reviewed for CIAO 4.3: no changes 
15 Dec 2011  reviewed for CIAO 4.4: the figure has been updated to reflect the latest calibration of the HRCS QE and LETG higher order efficiencies. 
03 Dec 2012  Review for CIAO 4.5 
14 Dec 2013  Finished CIAO 4.5 review; updated link in POG. 
03 Dec 2013  Review for CIAO 4.6; no changes. 
18 Dec 2014  Reviewed for CIAO 4.7; no changes. 
30 Jan 2017  Reviewed for CIAO 4.9. Updated the equations to use mathjax. 