ModifiedBlackBodyLayers

class pysm3.ModifiedBlackBodyLayers(map_layers, freq_ref, map_mbb_index, map_mbb_temperature, nside, color_correction=1, max_nside=None, num_layers=1, unit_layers=None, unit_mbb_temperature=None, map_dist=None)

Bases: Model

Modified Black Body model with multiple layers

Used for the MKD 3D dust model by Ginés Martínez-Solaeche, Ata Karakci, Jacques Delabrouille: https://arxiv.org/abs/1706.04162

Parameters:

map_layers: `pathlib.Path`, str or ndarray

Path or string with a templated layer number {layer} to download or access locally a IQU map for each layer (1-based layer number) Alternatively an array of shape (num_layers, 3, num_pix)

num_layers: int

Number of layers

unit_*str or Unit

Unit string or Unit object for all input FITS maps, if None, the input file should have a unit defined in the FITS header.

freq_ref: Quantity or string

Reference frequencies at which the intensity and polarization templates are defined. They should be a astropy Quantity object or a string (e.g. “1500 MHz”) compatible with GHz.

map_mbb_index, map_mbb_temperature: `pathlib.Path`, str or ndarray

Path or string with a templated 1-based layer number {layer} with the spectra index or the blackbody temperature. Alternatively an array of shape (num_layers, num_pix)

nside: int

Resolution parameter at which this model is to be calculated (with ud_grade)

color_correction: float

Scalar correction factor multiplied to the maps, implemented to add a color correction factor to Planck HFI 353 GHz maps

Methods Summary

get_emission(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(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:

outputastropy.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.