Tools module

tools.get_blsqinv(beamval, el, make_2d=0, mapparams=None)

Get the inverse of the beam window function sqaured.

Parameters:

• beamval (float) – Beam FWHM in arcmins.

• el (array) – Multipoles over which the window function must be defined.

• make_2d (bool) – Convert to 2D if desired. Default is False.

• mapparams (list) – [nx, ny, dx, dy] where ny, nx = flatskymap.shape; and dy, dx are the pixel resolution in arcminutes. for example: [100, 100, 0.5, 0.5] is a 50’ x 50’ flatskymap that has dimensions 100 x 100 with dx = dy = 0.5 arcminutes.

Returns:

blsqinv – 1/Bl^2 either in 1d or 2D.

Return type:

array.

tools.get_nl(noiseval_in_ukarcmin, el, beamval=None, elknee_t=-1, alpha_knee=0)

Get the noise power spectra: White and 1/f spectrum. Can return beam deconvolved nl if desired.

\[P(f) = A^2 \left[ 1+ \left( \frac{\ell_{\rm knee}}{\ell}\right)^{\alpha_{\rm knee}} \right].\]

Parameters:

• noiseval_in_ukarcmin (float) – White noise level in uK-arcmin.

• el (array) – Multipoles over which the window function must be defined.

• beamval (float) – Beam FWHM in arcmins. Default is None. If supplied, bl^2 will be divided from nl.

• elknee_t (float) – Knee frequency for 1/f (el_knee in the above equation).

• alpha_knee (float) – Slope for 1/f (alpha_knee in the above equation).

Returns:

nl – (Beam deconvoled) noise power spectrum.

Return type:

array.