usc-cosmolab / hackspace

Materials from CosmoLab group meetings.
MIT License
4 stars 16 forks source link

lya/theory.py #17

Closed maamari closed 4 years ago

maamari commented 4 years ago

PR for theory.py using Pk = Norm k^ns Tk^2 where Tk is that of Eisentein/Hu (https://arxiv.org/pdf/astro-ph/9709112.pdf) and normalization is that of https://ned.ipac.caltech.edu/level5/March02/White/White3.html

timothydmorton commented 4 years ago

This seems more complicated than what we did here: https://github.com/usc-cosmolab/hackspace/blob/master/notebooks/lyman_alpha.ipynb.

Regarding normalization, I don't think there's an objectively right normalization to have; the theory should just match the normalization of the data. Have you plotted the get_lcdm_pk against the data?