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

cross correlation > 1 because of forced shot noise subtraction #16

Open
adrianbayer opened this issue Sep 26, 2024 · 1 comment
Open

Comments

@adrianbayer
Copy link
Contributor

There are cases when the cross correlation coefficient computed using MeshFFTCorrelator.__call__ is larger than 1 which is mathematically questionable. It happens because shot noise is removed from the autospectra. It would be nice to have an option to not remove the shot noise, or more generally, use the kwargs when computing the ratio in __call__. I believe all that would need to be changed is this line
https://github.com/cosmodesi/pyrecon/blob/26f5e828fe0619afe303fa4e7d137095c63f7f15/pyrecon/metrics.py#L148C9-L148C14
from
power = self.ratio
to
power = self.get_ratio(complex=complex, **kwargs)
Let me know what you think, thanks!

@adematti
Copy link
Member

Hi Adrian, looks good --- please file a PR, or merge directly, as you prefer!
Thanks!

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