SpDustPol

class pysm.SpDustPol(map_I, freq_ref_I, emissivity, freq_peak, freq_ref_peak, pol_frac, angle_Q, angle_U, nside, unit_I=None, map_dist=None)

Bases: pysm.SpDust

SpDust2 model with Polarized emission

Methods Summary

get_emission(self, freqs[, weights])

This function evaluates the component model at a either a single frequency, an array of frequencies, or over a bandpass.

Methods Documentation

get_emission(self, freqs: Unit("GHz"), weights=None)

This function evaluates the component model at a either a single frequency, an array of frequencies, or over a bandpass.

Parameters:

freqs: scalar or array astropy.units.Quantity

Frequency at which the model should be evaluated, in a frequency which can be converted to GHz using astropy.units. If an array of frequencies is provided, integrate using trapz with a equal weighting, i.e. simulate a top-hat bandpass.

weights: np.array, optional

Array of weights describing the frequency response of the instrument, i.e. the bandpass. Weights are normalized and applied in Jy/sr.

Returns:

output : astropy.units.Quantity

Simulated map at the given frequency or integrated over the given bandpass. The shape of the output is (3,npix) for polarized components, (1,npix) for temperature-only components. Output is in uK_RJ.