Skip to content

A Pythonic SHT-based quadratic estimator reconstruction code that can work efficiently on cut skies. Primarily useful for CMB lensing, but also for a wide variety of other applications involving QEs including kSZ tomography.

License

Notifications You must be signed in to change notification settings

simonsobs/falafel

Repository files navigation

falafel

This library allows you to reconstruct mode-couplings (e.g. CMB lensing, kSZ) on cut skies. It is SHT-based and uses a cylindrical projection in intermediate steps which allows it to be fast and memory efficient. The implementation of the lensing quadratic estimators closely follows that in the Planck 2018 lensing analysis by Julien Carron and others. SHT-based estimators for shear, point source and mask hardening have also been implemented by Frank Qu.

Only the unnormalized quadratic estimate is provided. For the full-sky normalization for various estimators, we recommend tempura (https://github.com/simonsobs/tempura).

Installation

Use

python setup.py develop

(for in-place development)

or

python setup.py install

(Note that pip install may not work)

A good place to start with the code is tests/simple.py.

Copy input/config_template.yaml to input/config.yaml and edit the latter to point to paths in your filesystem.

About

A Pythonic SHT-based quadratic estimator reconstruction code that can work efficiently on cut skies. Primarily useful for CMB lensing, but also for a wide variety of other applications involving QEs including kSZ tomography.

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published