Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Improve lss_theory #4

Open
1 of 9 tasks
damonge opened this issue Mar 23, 2017 · 3 comments
Open
1 of 9 tasks

Improve lss_theory #4

damonge opened this issue Mar 23, 2017 · 3 comments

Comments

@damonge
Copy link
Collaborator

damonge commented Mar 23, 2017

A possibly not complete list of things that should be implemented in lss_theory.py is:

  • generalize to correlation functions. This is not currently in CCL, but will soon be
  • convolve multipoles with window function (currently C_ell only computed at band centres)
  • add flags for RSDs and magnification (magnification bias should then be passed as parameters too)
  • implement a cleverer/more general way of passing parameters (?)
  • make bias interpolation order a parameter (?)
  • make sure we're happy with the way bias parameters are currently being passed
  • do we want to start implementing other tracers (e.g. cmb, lensing...)?
  • benchmark likelihood evaluation
  • include simple photo-z systematics (?)
@damonge
Copy link
Collaborator Author

damonge commented Mar 23, 2017

I imagine the answer to the penultimate one is no.

Regarding the last one: currently the bottleneck is the calls to CCL (which is good news). Within CCL, if you use a fast power spectrum (e.g. Eisenstein & Hu or BBKS), the bottleneck is the LOS integrals. If you use CLASS, the bottleneck is actually the call to CLASS.

@jeremyneveu
Copy link
Contributor

What does the second point mean exactly ? Sky window or redshift window ?

@damonge
Copy link
Collaborator Author

damonge commented Apr 18, 2017

Basically, after having deconvolved any sky mask/weights/contaminants you'll be left with an estimate of the power spectrum at certain ells which is a convolution of the true C_ell with some window function (in ell). SACC has some capability to store these, but we're not using them (we just assume it's enough to evaluate the theoretical C_ell at the centre of each bandpower).
This is probably not a super bad approximation, so I would focus on the other tasks for the time being.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants