plancklens.bandpowers
Band-powers construction module.
This is used to construct anistropy band-powers from a parameter file.
- class plancklens.bandpowers.ffp10_binner(k1, k2, parfile, btype, ksource='p')[source]
Band-power construction library using the FFP10 fiducial cosmology.
This combines lensing (or anisotropy) estimates spectra to build band-powers in the exact same way than the Planck 2018 lensing analysis.
- This uses the various QE and QE spectra libraries defined in a parameter file, in particular:
qcls_dd (for data band-powers, to build covariance matrix, Monte-Carlo and point-source correction)
qcls_ds (for RDN0 and point-source correction)
qcls_ss (for MCN0, RDN0, Monte-Carlo and point-source correction)
qresp_dd (for the estimator normalization)
n1_dd (for the N1 bias subtraction)
nhl_dd (to build the semi-analytical covariance matrix)
ivfs (for the N1 bias and point-source correction)
In each of the methods defined here (e.g. MCN0, RDN0…), if the relevant QE, QE spectra, etc cannot be found precomputed, this will be performed on the fly. Hence in a realistic configuration it is always advisable to build them all previously (for example with examples/run_qlms.py)
This library can be used to build the cross power spectra of two anisotropy estimators, calculates biases, obtain MC and point-source corrections and the covariance matrix.
- Parameters:
k1 – First quadratic estimator key. See the qest.py module for the key definitions.
k2 – Second quadratic estimator key. See the qest.py module for the key definitions.
parfile – parameter file where the relevant QE libraries are defined
btype – bin type descriptor (‘consext8’ or ‘arg2’ were the Planck 2018 lensing analysis defaults)
ksource – anisotropy source (defaults to ‘p’, lensing)
- get_ampl_x_input(mc_sims=None)[source]
Returns cross-correlation of phi-maps to input lensing maps.
Uses qlms_x_i library of parfile
- get_bamc(wn1=True)[source]
Binned additive MC correction, with crude error bars.
This compares the reconstruction on the simulations to the FFP10 input lensing spectrum.
Note
the approximate error corrections to the additive MC correction variance follows Appendix C of https://arxiv.org/abs/1807.06210, check this for more details on its validity.
- get_bmmc(mc_sims_dd=None, mc_sims_ss=None, wN1=True)[source]
Binned multiplicative MC correction.
This compares the reconstruction on the simulations to the FFP10 input lensing spectrum.
- get_dat_bandpowers()[source]
Returns data raw band-powers, prior to any biases subtraction or correction.
- get_dat_nhl()[source]
Returns N0 lensing bias, semi-analytical version.
This is not highly accurate on the cut-sky
- get_mcn0_cov(mc_sims_dd=None)[source]
Covariance matrix obtained from the realization-independent debiaser.
- get_n1(k1=None, k2=None, unnormed=False)[source]
Returns analytical N1 lensing bias.
This uses the analyical approximation to the QE pair filtering as input.