srcmodel.polarization

ximpol.srcmodel.polarization.constant(C)[source]

Simple wrapper returning a constant, independently of the input arguments.

class ximpol.srcmodel.polarization.xPolarizationMap(xmap_file_path, ymap_file_path)[source]

Read-mode interface for the polarization maps.

build_grid_sample(ra0, dec0, num_points=25, radius=5.0)[source]

Calculate the polarization components on a rectangular grid.

Warning

This could be probably vectorized.

build_random_sample(ra0, dec0, num_points=1000, radius=5.0)[source]

Calculate the polarization components on a random set of points.

Warning

This could be probably vectorized.

overlay_arrows(fig, markers=True)[source]

Overlay the polarization map arrows over an existing aplpy figure.

plot_polarization_angle(degrees=True, show=True)[source]
plot_polarization_degree(show=True)[source]
plot_xmap(overlay=True, show=False)[source]

Plot the x polarization map.

plot_ymap(overlay=True, show=False)[source]

Plot the y polarization map.

polarization_angle(ra, dec)[source]

Return the polarization angle for a given direction in the sky.

polarization_degree(ra, dec)[source]

Return the polarization degree for a given direction in the sky.

Warning

Note that we’re calling the polarization_vector() function here and in the polarization_angle() method, while we could in principle get away with just one function call. The issue is that downstream we need the polarization degree and angle separately, and if we want to optimize things here, we would have to implement a generic polarization() interface in the model components that is called consistently in rvs_event_list().

polarization_vector(ra, dec)[source]

Return the polarization vector for a given (array of) RA and Dec values.