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

mocks for a non-gaussian field. #660

Open
Jayashree-Behera opened this issue Dec 13, 2021 · 7 comments
Open

mocks for a non-gaussian field. #660

Jayashree-Behera opened this issue Dec 13, 2021 · 7 comments

Comments

@Jayashree-Behera
Copy link

Hi, I am working on a project based on non-gaussian fields. So while generating mocks using nbodykit (let's say lognormal mocks), one requires to give the matter power spectrum(as a functional argument ) and bias (as a number) as input (among other parameters). This is fine as long as we are working are in a gaussian field but when it comes to non-gaussian field, bias is no longer an independent number but a function of k, hence a functional argument. I was wondering if this feature could be added where one could generate mocks in non-gaussian fields. I have some working code that I could share if you are interested. :)
Thanks,
Shree

@Jayashree-Behera Jayashree-Behera changed the title Galaxy power spectrum for a non-gaussian field. mocks for a non-gaussian field. Dec 13, 2021
@rainwoodman
Copy link
Member

Thanks for reaching out. Yes code example will be very helpful!

Did you mean making the b parameter of LogNormalCatalog scale dependent?
The b parameter in LogNormalCatalog is applied on configuration space -- how do you get around that?

Or did you mean introducing a scale dependent bias to the linear field (gaussian field), before Log Normal transformation? In that case, the scale dependent biased linear field just need a different powerspectrum or transfer function -- the interaction with the b parameter of LogNormalCatalog will be even more intriguing in this set up.

@Jayashree-Behera
Copy link
Author

Jayashree-Behera commented Dec 20, 2021

Hi... Thanks for replying.
My code basically lies in the lines of this "..Or did you mean introducing a scale-dependent bias to the linear field (gaussian field), before Log-Normal transformation? In that case, the scale-dependent biased linear field just needs a different powerspectrum or transfer function --..."

I wrote a script for the galaxy power spectrum in a non-gaussian field (which incorporates the scale-dependent bias). Then I generated the mocks by giving the galaxy power spectrum as the input parameter setting bias =1.
Please let me know if this wasn't clear.

Have a good day,
Shree

@rainwoodman
Copy link
Member

rainwoodman commented Dec 22, 2021 via email

@Jayashree-Behera
Copy link
Author

Sure. Should I make a pull request and you can check it out there?

@Jayashree-Behera
Copy link
Author

Hi.. Did you get a chance to look at the pull request? Does it make sense?

@rainwoodman
Copy link
Member

rainwoodman commented Jan 12, 2022 via email

@Jayashree-Behera
Copy link
Author

Hi..Did you get a chance to review the updated version in PR? I have made some modifications and added the unit test file.

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