An implementation quadratic estimator for lensing extraction of full-sky CMB data in python, following Okamoto & Hu's paper CMB Lensing Reconstruction on the Full Sky.
The module requires an installation of HEALPix C++ (version 3.70). Furthermore, the python modules numpy, cython and healpy should be installed.
The module can be easily installed with pip by executing
pip install . [--user]
in your local lensquest directory.
maps: CMB maps (T or list of T,Q,U) or corresponding harmonic coefficients (T or list of T,E,B) in healpy formatwcl: Power spectra used in the weights of the quadratic estimator (array of TT or TT,EE,BB,TE power spectra)dcl: Power spectra used in the Wiener-filter of the input fields (array of TT or TT,EE,BB,TE power spectra)
import lensquest
questhelper=lensquest.quest(maps, wcl, dcl, lmin=2, lmax=None, lminCMB=2, lmaxCMB=None)
questhelper.grad(XY)
# returns a_lm^Phi XY, where XY=TT,TE,EE,TB,EB or BB
lensquest.quest_norm(wcl, dcl, lmin=2, lmax=None, lminCMB=2, lmaxCMB=None, bias=False)
# returns dictionary of A_L (and N_L if bias=True) of TT or TT,TE,EE,TB,EBDominic Beck: dbeck [at] apc.in2p3.fr
This project is licensed under the MIT License - see the LICENSE.md file for details