Model¶

class
pysm.
Model
(nside, map_dist=None, dataurl=None)[source] [edit on github]¶ Bases:
object
This is the template object for PySM objects.
If a MPI communicator is passed as input and
pixel_indices
is None, the class automatically distributes the maps across processes. You can implement your own pixel distribution passing both a MPI communicator andpixel_indices
, however that won’t support smoothing withlibsharp
. Iflibsharp
is available, the rings are distributed as expected bylibsharp
to perform distributed spherical harmonics transforms, seepysm.mpi.distribute_rings_libsharp()
, thelibsharp
grid object is saved inself.libsharp_grid
. If libsharp is not available, pixels are distributed uniformly across processes, seepysm.mpi.distribute_pixels_uniformly()
Parameters:  mpi_comm: object
MPI communicator object (optional, default=None).
 nside: int
Resolution parameter at which this model is to be calculated.
 smoothing_lmax : int
\(\ell_{max}\) for the smoothing step, by default \(2*N_{side}\)
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. read_map
(self, path[, unit, field])Wrapper of the PySM read_map function that automatically uses nside, pixel_indices and mpi_comm defined in this Model read_txt
(self, path, \*\*kwargs)Methods Documentation

get_emission
(self, freqs: Unit("GHz"), weights=None) > Unit("uK_RJ")[source] [edit on github]¶ 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 tophat 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 temperatureonly components. Output is in
uK_RJ
.

read_map
(self, path, unit=None, field=0)[source] [edit on github]¶ Wrapper of the PySM read_map function that automatically uses nside, pixel_indices and mpi_comm defined in this Model

read_txt
(self, path, **kwargs)[source] [edit on github]¶