Abstract
Abstract We present the XFaster analysis package, a fast, iterative angular power spectrum estimator based on a diagonal approximation to the quadratic Fisher matrix estimator. It uses Monte Carlo simulations to compute noise biases and filter transfer functions and is thus a hybrid of both Monte Carlo and quadratic estimator methods. In contrast to conventional pseudo-C ℓ –based methods, the algorithm described here requires a minimal number of simulations and does not require them to be precisely representative of the data to estimate accurate covariance matrices for the bandpowers. The formalism works with polarization-sensitive observations and also data sets with identical, partially overlapping, or independent survey regions. The method was first implemented for the analysis of BOOMERanG data and also used as part of the Planck analysis. Here we describe the full, publicly available analysis package, written in Python, as developed for the analysis of data from the 2015 flight of the Spider instrument. The package includes extensions for self-consistently estimating null spectra and estimating fits for Galactic foreground contributions. We show results from the extensive validation of XFaster using simulations and its application to the Spider data set.
Published Version (Free)
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.