This repository contains some tools to perform the last steps of a cosmological analysis of the 1D power spectrum (P1D) of the Lyman alpha forest.
It uses the LaCE emulator (https://github.com/igmhub/LaCE), and some extra tools to run MCMC analyses on cosmological and IGM parameters for a mock P1D measurement.
If you would like to collaborate, please email Andreu Font-Ribera (afont@ifae.es)
Start a fresh environment:
conda create -n cup1d python=3.10
conda activate cup1d
Clone the cup1d repo, set a CUP1D_PATH environment variable pointing to it:
cd $CUP1D_PATH
pip install -e .
Note: This will replace your local LaCE installation if you have any. To use your editable LaCE installation, uninstall LaCE and reinstall your own copy.
pip uninstall lace
cd $LACE_PATH
pip install -e .
You can start by plotting the many P1D measurements stored in the repo, by looking at notebooks/p1d_measurements
You can also redo old neutrino mass constraints by importance sampling WMAP and Planck chains, following notebooks/planck
You can also play with the LaCE emulator with the notebooks in notebooks/emulator
Finally, you can run your own analysis on mock data following the notebooks in notebooks/likelihood
You can run the script under scripts/forecast.py to forecast the constraints on linear power parameters for a given P1D covariance.
It marginalizes over 8 nuisance parameters, so it might take a while to run!