skymodel.simulate_cmb

skymodel.simulate_cmb(freq, cl_file=None, lensed=True, nside_out=4096, lmax=None, beam_FWHM=None, template='WebSky', unit='cmb')

Function that computes a CMB map from a power spectrum.

Parameters

freq: float or float array

Frequency of the output map in Hz.

cl_file: str or array, optional

Name of the .dat contaning the values of the power spectrum given by CAMB. Or array containing the power spectrum to generate random maps in Kelvin. Default: None

lensed: bool, optional

If True select the lensed CMB map when possible. This is only possible for ‘WebSky’ and ‘Sehgal.’ Default: True

nside_out: float, optional

Healpix nside parameter of the output map. Must be a valid value for nside. Default: 4096

lmax: float, optional

Maximum value of the multipolemoment at which the atmospheric power spectrum wil be computed. Default: 3*nside_out-1

beam_FWHM: bool, optional

If set, the output will be convolved with a gaussian. The FWHM of the Gaussian in units of arcmin is given by the provided value. Default: None

template: bool, optional

Determines the all-sky foregrounds templates to be used to build the sky model. If ‘Sehgal’ is chosen, simulations by Sehgal et al. (2010) are used. If ‘SO’ is chosen, the Simons Observatory sky model provided by Colin Hill and based on the simulations by Sehgal et al. (2010) is used. If ‘WebSky’ is chosen, the used templates will be based on the WebSky Extragalactic CMB Mocks provided by CITA. If ‘SO_reproduced’ is chosen, the SO sky model is reproduced directly from the Sehgal et al. (2010) data. Default: ‘WebSky’

unit: bool, optional

Determines the units of the output map. The available units are ‘mjy’ –> MJy/sr (specific intensity), ‘cmb’ –> K_CMB (thermodynamic temperature), and ‘rj’ –> K_RJ (brightness temperature). Default: ‘cmb’.

Returns

float array

Healpix allsky map contaning the CMB map at a given frequency.