Pixel window functions

A pixelated signal f(p) is the average within each pixel p (with surface area $\Omega_{\mathrm{pix}}$) of the underlying signal

f(p) $\textstyle \myequal $ $\displaystyle \int d\u w_p(\u )f(\u )$ (29)

where wp is equal to $1/\Omega_{\mathrm{pix}}$ within the pixel, and equal to 0 outside, so that $\int d\u w_p(\u ) = 1$. Eq. (4) then becomes
f(p) $\textstyle \myequal $ $\displaystyle \sum_{\ell =0}^{\ell_{\mathrm{max}}}\sum_{m}a_{\ell m}w_{\ell m}(p),$ (30)

where
$\displaystyle w_{\ell m}(p)$ $\textstyle \myequal $ $\displaystyle \int d\u w_p(\u ) Y_{\ell m}(\u ),$ (31)

is the Spherical Harmonic Transform of the pixel p.

However, complete analysis of a pixelated map with the exact $w_{\ell m}(p)$ defined above would be computationally intractable (because of azimutal variation of pixel shape over the polar caps of the HEALPix grid), and some simplifying asumptions have to be made. If the pixel is small compared to the signal correlation length (determined by the beam size), the exact structure of the pixel can be ignored in the subsequent analysis and we can assume

$\displaystyle w_{\ell m}(p)$ $\textstyle \myequal $ $\displaystyle w_\ell(p) Y_{\ell m}(p)$ (32)

where we introduced the m-averaged window function
$\displaystyle w_{\ell}(p)$ $\textstyle \myequal $ $\displaystyle \left(\frac{4 \pi}{2\ell+1}\sum_{m=-\ell}^{\ell} \left\vert w_{\ell m}(p)\right\vert^2\right)^{1/2},$ (33)

which is independent of the pixel location on the sky.

If we assume all the pixels to be identical, the power spectrum of the pixelated map, $C_{\ell}^{\mathrm{pix}}$, is related to the hypothetical unpixelated one, $C_{\ell}^{\mathrm{unpix}}$, by

$\displaystyle C_{\ell}^\mathrm{pix}$ $\textstyle \myequal $ $\displaystyle w^2_{\ell} C_{\ell}^\mathrm{unpix}$ (34)

where the effective pixel window function $w_{\ell}$ is defined as
$\displaystyle w_{\ell}$ $\textstyle \myequal $ $\displaystyle \left(\frac{1}{N_{\mathrm{pix}}}\sum_{p=0}^{N_{\mathrm{pix}}-1} w^2_{\ell}(p)\right)^{1/2}.$ (35)

This function is provided with the HEALPix package for $\ell\le 4N_{\mathrm{side}}$ for each resolution parameter Nside.

The pixel window functions are now available for both temperature and polarization.

For $N_{\mathrm{side}}\le 128$, those window functions are computed exactly using Eqs. (33) and (35). For Nside> 128 the calculations are too costly to be done exactly at all $\ell$. The temperature windows are extrapolated from the case Nside= 128 assuming a scaling in $\ell$ similar to the one exhibited by the window of a tophat pixel. The polarization windows are assumed to be proportional to those for temperature, with a proportionality factor given by the exact calculation of $w_{\ell}$ at low $\ell$.

Because of a change of the extrapolation scheme used, the temperature window functions provided with HEALPix 1.2 and higher for Nside> 128 are slighty different from those provided with HEALPix 1.1. For a given Nside, the relative difference increases almost linearly with $\ell$, and is of the order of $\Delta w/w < 7\ 10^{-4}$ at $\ell=2N_{\mathrm{side}}$ and $\Delta w/w < 1.7\ 10^{-3}$ at $\ell=4N_{\mathrm{side}}$.

Version 3.50, 2018-12-10