Accounting for Calibration Uncertainties in X-ray Analysis: Effective Areas in Spectral Fitting

Abstract

While considerable advance has been made to account for statistical uncertainties in astronomical analyses, systematic instrumental uncertainties have been generally ignored. This can be crucial to a proper interpretation of analysis results because instrumental calibration uncertainty is a form of systematic uncertainty. Ignoring it can underestimate error bars and introduce bias into the fitted values of model parameters. Accounting for such uncertainties currently requires extensive case-specific simulations if using existing analysis packages. Here we present general statistical methods that incorporate calibration uncertainties into spectral analysis of high-energy data. We first present a method based on multiple imputation that can be applied with any fitting method, but is necessarily approximate. We then describe a more exact Bayesian approach that works in conjunction with a Markov chain Monte Carlo based fitting. We explore methods for improving computational efficiency, and in particular detail a method of summarizing calibration uncertainties with a principal component analysis of samples of plausible calibration files. This method is implemented using recently codified Chandra effective area uncertainties for low-resolution spectral analysis and is verified using both simulated and actual Chandra data. Our procedure for incorporating effective area uncertainty is easily generalized to other types of calibration uncertainties.

61 pages double spaced, 8 figures, accepted for publication in ApJ

References in corpus (7)

• Supernova Constraints and Systematic Uncertainties from the First 3 Years of the Supernova Legacy Survey
• The Bolocam Galactic Plane Survey: Survey Description and Data Reduction
• Calibration of the LIGO Gravitational Wave Detectors in the Fifth Science Run
• Type Ia Supernova Light Curve Inference: Hierarchical Bayesian Analysis in the Near Infrared
• X-ray Properties of the GigaHertz-Peaked and Compact Steep Spectrum Sources
• Searching for Narrow Emission Lines in X-ray Spectra: Computation and Methods
• An X-ray photometry system I: Chandra ACIS

Cited by in corpus (26)

• A NICER View of PSR J0030+0451: Millisecond Pulsar Parameter Estimation
• Stellar Population Inference with Prospector
• The Fermi Large Area Telescope On Orbit: Event Classification, Instrument Response Functions, and Calibration
• The stellar rotation--activity relationship in fully convective M dwarfs
• Constraining the Mass and Radius of Neutron Stars in Globular Clusters
• Lensing and X-ray mass estimates of clusters (SIMULATION)
• Strong Lens Time Delay Challenge: II. Results of TDC1
• Neutron Star Mass-Radius Constraints of the Quiescent Low-mass X-ray Binaries X7 and X5 in the Globular Cluster 47 Tuc
• Deep Chandra observations of Pictor A
• Metropolis-Hastings within Partially Collapsed Gibbs Samplers
• Estimating the properties of hard X-ray solar flares by constraining model parameters
• Conservative cosmology: combining data with allowance for unknown systematics
• dart_board: Binary Population Synthesis with Markov Chain Monte Carlo
• Incorporating Uncertainties in Atomic Data Into the Analysis of Solar and Stellar Observations: A Case Study in Fe XIII
• Spatially Resolved Chandra Spectroscopy of the Large Magellanic Cloud Supernova Remnant N132D
• Disentangling Overlapping Astronomical Sources using Spatial and Spectral Information
• Time Variability of Nonthermal X-ray Stripes in Tycho's Supernova Remnant with Chandra
• Chandra Imaging of the Western Hotspot in the Radio Galaxy Pictor A: Image Deconvolution and Variability Analysis
• The Proper Motion of the Pulsar J1124-5916 in the Galactic Supernova Remnant G292.0+1.8
• Concordance: In-flight Calibration of X-ray Telescopes without Absolute References
• Long-term X-ray Variability of the Symbiotic System RT Cru based on Chandra Spectroscopy
• Systematic errors in the maximum likelihood regression of Poisson count data: introducing the overdispersed chi-square distribution
• X-ray Spectral Analysis of the Jet Termination Shock in Pictor A on Sub-Arcsecond Scales with Chandra
• Estimating Global Ejecta Mass Ratios in Tycho's SNR
• Estimating Ejecta Mass Ratios in Kepler's SNR: Global X-Ray Spectral Analysis Including Suzaku Systematics and Emitting Volume Uncertainties
• Six Maxims of Statistical Acumen for Astronomical Data Analysis