Документ взят из кэша поисковой машины. Адрес оригинального документа : http://hea-www.harvard.edu/AstroStat/etc/CCW2007_109_drkeijmsv.pdf Дата изменения: Tue Oct 30 16:53:43 2007 Дата индексирования: Tue Oct 2 04:30:04 2012 Кодировка: Поисковые слова: ngc 4303

Monte Carlo Methods for Including Calibration Uncertainties in Model Fitting Analyses
Jeremy J. Drakea , Peter Ratzlaffa , Vinay Kashyapa , Richard Edgara , Rima Izemb , Diab Jeriusa , Herman Marshallc , Aneta Siemiginowskaa and Alexey Vikhlinina
a Smithsonian Astrophysical Observatory, Cambridge MA 02138 Department of Statistics, Harvard University, Cambridge, MA 02138 MIT Kavli Institute for Astrophysics and Space Research, Cambridge, MA 02139 b c

SUMMARY · Instrument response uncertainties almost universally ignored in astrophysical X-ray data analyses, yet for good quality observations can be dominant source of error. · Response uncertainties are correlated; both understanding and specifying the uncertainties is technically challenging. Moreover, there is no standard set of procedures for incorp orating complicated correlated systematic uncertainties in non-linear parameter estimation (eg XSPEC fitting): the approaches used for treating indep endent errors simply do not apply. · We have developed Monte Carlo methods to treat calibration uncertainties for the Chandra High Resolution Mirror Assembly (HRMA) and Advanced CCD Imaging Spectrograph (ACIS). The code and ancilliary data will be released to Chandra Users upon acceptance for publication of the article describing this work. CIAO Sherpa methods are also under development to utilise these techniques (see accompanying poster by Kashyap et al).
Main Uncertainties in Instrument Response: Chandra ACIS-S
HRMA
- - - - Geometry Obscuration Reflectivity Scattering

Figure 2. Left: Illustration of the relative change in the HRMA effective area caused by different hydrocarbon contamination layers. The range shown corresponds to the nominal adopted 22 ± 6 ° layer thickness. Middle: Relative A changes in the model ACIS S3 QE caused by ±20% differences in the model CCD depletion depth (note that ±13 % was adopted here, which corresponds to a range of about ±10% in the QE at 10 keV). Right: QE changes caused by the adopted ±20% differences in CCD SiO2 thickness.

hi
Perturbation Factor

1

lo

Energy

edge
E hi

E lo

Figure 1. Left: Truncated normal distribution (product of a Gaussian with
variance 2 and a rectangular Step function with unit density between ± ), P ( ), representing the distribution of calibration uncertainties used in the perturbation function and Monte Carlo draws. Right: Illustration of a perturbation function segment used to apply random deviations from a nominal subassembly response within a given energy range. Within each energy range, Elo ­Ehi, a random low-order polynomial ( 3) is generated that is constrained to lie within the grey shaded region defined by the uncertainties lo and hi, and also to join up with neighbouring segments within the edge constraints edge. The deviation from unity for a large sample of vectors corresponds to P ( ).

ACIS OBF
- Transmittance - Contamination - Uniformity

ACIS CCD
- QE - Gain - Pulse Height Distribution

Figure 3. The nominal "seed" Chandra ACIS-S effective area (black) compared with a sample of 30 effective areas generated using the Monte Carlo modification method described in the text (grey).

METHODS
Construct different realisations of instrument response by a combination of (1) randomly varying input parameters describing subassembly performance and (2) random multiplicative perturbation functions, µ(E ), designed to sample subassembly responses with their assessed uncertainties (Fig. 1). Adopt "curtailed Gaussian" probability distribution P ( ) for Monte Carlo draws (Fig. 1a). The different subassemblies were treated as follows:

HRMA On-Axis:

A combination of perturbation functions, µH (E ), running within prominent Ir edges, and raytrace-derived effective areas sampling the effects of different hydrocarbon contamination layers with the measured range of allowed values (Fig. 2).

HRMA Vignetting Function: For off-axis angles (in
arcmin) the area perturbation function was multiplied by a combination of a fractional uncertainty of the azimuthally-averaged vignetting function, Ї V (), and an expression involving the ratio of Debye-Waller to perfect mirror reflectivities: µv (E , ) = P (v )(1 - V Ї)) + P (s)(1 - RDW /R); v , s = 0.2. (

ACIS OBF and Contamination Layer: OBF transmittance uncertainty employed perturbation functions, µOB F (E ), between C-N and N-O edges, and O-10 keV constrained by different allowed maximum deviations and relative edge transmittance discontinuity errors. The contamination perturbation function is: µC L(E ) = e-P (C )C +P (O )O +P (F )F +P (F l )F l ; µC L(0.7k eV ) < 0.05 where C , O , F and F l are the fractional uncertainties in the optical depths C, O, F and Fluffium at a fiducial date (2003.29).

Figure 4. Example frequency distributions of best-fit parameters obtained
for typical blackbody, thermal plasma and powerlaw models from XSPEC for synthetic data sets containing 104 (upper panels) and 105 (lower panels). Black histograms are distributions resulting from 1000 Monte Carlo samplings of the synthetic data allowing Poisson noise variations alone. Red histograms are the distributions of parameters resulting from fits to a single synthetic data set using 1000 Monte Carlo-generated effective areas and response matrices.

ACIS QE: QE uncertainty base on combination of perturbation
functions, µQE (E ), with boundaries at O and Si K edges and ACIS QE model predictions for uncertainties of 13% in CCD depletion depth and 20% in SiO2 thickness.
Figure 5. Modes and highest posterior density ±95% confidence intervals obtained for the blackbody, thermal plasma and powerlaw models investigated using XSPEC for synthetic data sets containing 105 counts. The y-axes correspond in all cases to the ratio of the input parameter to that retrieved in the model fit. Black error bars correspond to 1000 Monte Carlo samplings of the synthetic data and show the effects of Poisson noise variations alone. Dashed error bars correspond to fits to a single synthetic data set using 1000 Monte Carlo-generated effective areas and response matrices.

ACIS Gain and Pulse Height Distribution: Program calcrmf2 used to generate CCD Gain and pulse height response matrix files (RMF) for P (G) variations in gain and pulse height width; G =1% @0.7 keV, 0.5% @1.5 keV, and 0.2% @ 4 keV.

ESTIMATING EFECTS OF CALIBRATION UNCERTAINTIES
Perturb nominal effective area, A(E ), with combined perturbation functions, µ(E ), to make A(E ) (Fig. 3): A (E ) = µH (E )µV (E )µ
CL

(E )µ

OB F

(E )µ

QE

(E )A(E )

For each effective area and RMF, find the parameters of the best-fit model for a synthetic Chandra ACIS observation computed using the nominal instrument response. Compare with parameters found from fits to 1000 synthetic spectra differing only by Poisson noise and generated using the nominal area and RMF (Figs. 4 and 5). This step utilized the XSPEC fitting engine and investigated Blackbody, optically-thin thermal plasma, and power law continuum models. Rep eat 1000 times.

· Limiting accuracy of Chandra ACIS reached in spectra with 10

4

counts. Beyond this, errors in in best-fit parameters due to calibration uncertainties completely dominate those due to photon noise.