- Optimally combines data from different bands using noise + foreground signals in different bands.
- Supports standard / constrained / partial internal linear combinations.
- Computes lensing noise curves using residual noise.
- Combines delensed CMB spectra and lensing spectra to forecast cosmological parameter constraints using Fisher formalism.
- Delensing reference: Green, Meyers, van Engelen 2016, arXiv: 1609.08143.
- Delensing + Fisher code reference: Hotinli, Meyers, Trendafilova, Green, van Engelen 2021, arXiv: 2111.15036.
- Simple Fisher code repo without delensing.
- Estimates biases in cosmological parameters due to residual foregrounds also using Fisher formalism (yet to be integrated into this repo).
- CLASS_delens https://github.com/selimhotinli/class_delens (9944e0ec5f0b734b178436a8652b186046f63fdb)
This code uses a wrapper for the CLASS_delens code to facilitate Fisher forecasting of cosmological parameter constraints from CMB spectra.
- Selim C. Hotinli, Joel Meyers, Cynthia Trendafilova, Daniel Green, Alexander van Engelen
- Clone the repo first:
git clone [email protected]:sriniraghunathan/CMB-S4_DRAFT.git
- list the submodules:
git submodule
. This will show theCLASS_delens
and the specific commit9944e0ec5f0b734b178436a8652b186046f63fdb
- pull the submodule:
git submodule update --init --recursive
- Install CLASS:
cd CLASS_delens; make class
- Ensure CLASS is working properly by running the following document within the
CLASS_delens
folder:./class explanatory.ini
Parent script: perform_forecasts and this consists of two parts which can be run seprately as well.
- ILC code:
python3 get_ilc_weights_and_residuals.py -opfname results//s4_cmb_ilc.npy
or look into this ilc code example notebook - Delensing + Fisher code (and this uses the outout from the previous step):
python3 get_fisher_withdelensing.py -ilc_fname results//s4_cmb_ilc.npy -opfname results//s4_cmb_fisher.npy
- Results are stored within results folder.
- Look into (read_ilc.py)[https://github.com/sriniraghunathan/CMB-S4_DRAFT/tree/main/results/read_ilc.py] to read the ILC residuals.
- Look into perform_forecasts to read the Fisher matrices.
Joel Meyers, Cynthia Trendafilova, and Benjamin Wallisch.
- Selim C. Hotinli, Joel Meyers, Cynthia Trendafilova, Daniel Green, Alexander van Engelen
- Chilean LAT: S4-Wide
- Delensing SouthPole LAT: S4-Ultra deep
- Extragalactic foregrounds: Radio, CIB, tSZ and kSZ power spectra from SPT measurements (George et al. 2015, arXiv: 1408.3161 and Reichardt et al. 2020, arXiv: 2002.06197).
- Assumed polarisation fractions: CIB = 2%; Radio = 3%; tSZ/kSZ = 0. But these are configurable. Look into params.ini.
- Galactic foregrounds: Dust and Synchrotron power spectra obtained from pySM3 simulations.