so_window - a module for window function generation

routines for window function generation

pspy.so_window.apod_C1(binary, radius, rmax=None)

Create a C1 apodisation as defined in https://arxiv.org/pdf/0903.2350.pdf

Parameters:

• binary (so_map) – a so_map with binary data (1 is observed, 0 is masked)

• apo_radius (float) – the apodisation radius in degrees

pspy.so_window.apod_C2(binary, radius, rmax=None)

Create a C2 apodisation as defined in https://arxiv.org/pdf/0903.2350.pdf

Parameters:

• binary (so_map) – a so_map with binary data (1 is observed, 0 is masked)

• apo_radius (float) – the apodisation radius in degrees

pspy.so_window.create_apodization(binary, apo_type, apo_radius_degree, use_rmax=False)

Create a apodized window function from a binary mask.

Parameters:

• binary (so_map) – a so_map with binary data (1 is observed, 0 is masked)

• apo_type (string) – the type of apodisation you want to use (“C1”,”C2” or “Rectangle”)

• apo_radius (float) – the apodisation radius in degrees

pspy.so_window.get_distance(binary, rmax=None)

Get the distance to the closest masked pixels for CAR and healpix so_map binary.

Parameters:

• binary (so_map) –

• observed (a so_map with binary data (1 is) –

• masked) (0 is) –

pspy.so_window.get_spinned_windows(w, lmax, niter)

Compute the spinned window functions (for pure B modes method)

Parameters:

• w (so_map) – map of the window function

• lmax (integer) – maximum value of the multipole for the harmonic transform

• niter (integer) – number of iteration for the harmonic transform

pspy.so_window.get_survey_solid_angle(window, naive=False)

return Omega, the effective solid angle convered by the window function we use the window correction as displayed in equation 17 of https://arxiv.org/pdf/astro-ph/0105302.pdf to get f_sky just do:

fsky = Omega / (4 * np.pi)

Parameters:

• window (so_map) – map of the window function

• naive (boolean) – if True just compute the area covered by the non-zero pixels