skymodel.simulate_radio_ps

skymodel.simulate_radio_ps(freq, nside_out=4096, lmax=None, beam_FWHM=None, template='WebSky', unit='cmb')

Computes an all-sky radio point source map at a given frequency and nside based on the simulations provided by Sehgal et al. (2010), which have been recalibrated by the SO collaboration. The original simulations by Sehgal et al. are modelled by a curved power law and can be extrapolated to frequencies beyoned 350 GHz, while the sources in the SO simulations show more complex spectra due to ringing artifacts caused by the processing of the maps with alter_alm(). Due to these complex SEDs, bilinear interpolation is applied for input frequencies between 27 and 353 GHz. For higher frequencies, a null map is returned.

Parameters
freq: float or float array

Frequency of the output map in Hz.

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
radio_ps: float array

Healpix all-sky map of the radio point source emission at the specified frequency.