`powspec_analysis` - Full analysis¶

class powspechi.powspec_analysis(nside, angs, etacut=0.9, detcorr=False, multcorr_norm=None, multcorr_obs=None)¶

Perform the full analysis and extract the final angular power spectrum \(\langle S_{\ell}^{m \neq 0} \rangle\).

Parameters

nside (int, scalar) – Chosen map resolution.
angs (float, array_like) – List or array of angular coordinates of all events. Each of said events should be a 2-D array with shape (length, 2), where columns 0 and 1 stand for the azimuthal \(\phi\) and polar \(\theta\) coordinates, respectively.
etacut (float, scalar, optional) – Limit imposed on pseudorapidity, i.e., \(|\eta|\) < etacut. If there is no limit, just set it to None. Default: 0.9.
detcorr (bool, optional) – Determines whether to divide the observed maps by their ensemble sum. Usually used when there is a need to correct for detector’s non-uniform efficiency or to compare a heavy-ion model to the published angular power spectrum 1. Default: False.
multcorr_norm (filenamestr, optional) – File containing the averaged spectrum which corrects for event multiplicity the normalized spectrum. Format: four columns respectively representing the full averaged spectrum, its mean error, the averaged spectrum for \(m\neq0\) and its mean error. First line will be skipped. Also, it can only be activated if detcorr = True. Default: None.
multcorr_obs (filenamestr, optional) – File containing the averaged spectrum which corrects for event multiplicity the observed spectrum. Format is the same as in multcorr_norm. Default: None.

Variables

tmap (float, ndarray) – Sum of all maps within the event ensemble.
avcl_obs (dict) – Observed angular power spectrum averaged over all events. Its keys indicate the full spectrum and its mean error as well as the \(m\neq0\) spectrum and its mean error.
avcl_norm (dict, optional) – Normalized angular power spectrum averaged over all events. Its format is the same as avcl_obs. Activated when detcorr = True.
avcl_normcorr (dict, optional) – Normalized angular power spectrum corrected by event multiplicity. Its format is the same as avcl_obs. Activated when both detcorr = True and a file is given as multcorr_norm.
avcl_obscorr (dict, optional) – Observed angular power spectrum corrected by event multiplicity. Its format is the same as avcl_obs. Activated when a file is given as multcorr_obs.

Raises

SupmapError – When one wants to divide a single map by itself.
IsomapError – When the desired supmap_iso*.fits file does not exist. See maps_manip.getsupmapiso for more details.
PowSpecError – When one tries to correct a normalized spectrum when it does not exist.

powspec_analysis - Full analysis¶