Using peak distribution of the cosmic microwave background for WMAP and Planck data analysis: Formalism and simulations

¹ Max-Planck Institute für Astrophysik, Postfach 1317, 85741 Garching bei München, Germany e-mail: chm@mpa-garching.mpg.de
² SSAI, Code 685, Goddard Space Flight Center, Greenbelt, MD 20771, USA e-mail: kash@haiti.gsfc.nasa.gov
³ Facultad de Ciencias, Universidad de Salamanca. Pza de la Merced s/n, 37008, Spain e-mail: atrio@usal.es

Corresponding author: C. Hernández–Monteagudo, chm@mpa-garching.mpg.de

Received: 25 June 2003
Accepted: 7 October 2003

Abstract

We implement and further refine the recently proposed method at Kashlinsky et al. ([CITE]) (KHA) for a time efficient extraction of the power spectrum from future cosmic microwave background (CMB) maps. The method is based on the clustering properties of peaks and troughs of the Gaussian CMB sky. The procedure takes only steps where is the fraction of pixels with standard deviations in the map of N pixels. We use the new statistic introduced in KHA, , which characterizes spatial clustering of the CMB sky peaks of progressively increasing thresholds. The tiny fraction of the remaining pixels (peaks and troughs) contains the required information on the CMB power spectrum of the entire map. The threshold ν is the only parameter that determines the accuracy of the final spectrum. We performed detailed numerical simulations for parameters of the two-year WMAP and Planck CMB sky data including cosmological signal, inhomogeneous noise and foreground residuals. In all cases we find that the method can recover the power spectrum out to the Nyquist scale of the experiment channel. We discuss how the error bars scale with ν allowing to decide between accuracy and speed. The method can determine with significant accuracy the CMB power spectrum from the upcoming CMB maps in only ~ operations.