fix neutrino mass parameter in camb #118
Merged
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
Originally, omega_cdm (
omch2
) assumes all matter in cold dark matter and baryons, and ignores massive neutrinos:omch2=(self.cosmo.Om0 - self.cosmo.Ob0) * self.cosmo.h ** 2
I changed the definition to include massive neutrinos:
omch2=(self.cosmo.Om0 - self.cosmo.Ob0 **- self.cosmo.Onu0**) * self.cosmo.h ** 2
where I used astropy's cosmo.Onu0 calculation.
In addition, camb requires the neutrino mass and hierarchy.. I did not fix the code for different hierarchies (normal, inverted, degenerate) since it should be a sub-percent effect that's not observable by near future LSS surveys. So I assumed degenerate hierarchy (i.e. the 3 neutrinos are assumed to have equal mass = 1/3 *total mass):