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

Proton/Neutron/Electron ratio calculation incorrectly accounts for material composition #11

Open
austinschneider opened this issue Feb 28, 2021 · 0 comments
Assignees
Labels
bug Something isn't working

Comments

@austinschneider
Copy link
Member

In the SetMaterialModel function, the PNE ratios are computed by examining the PDG particle codes for each isotope in the material which lists the number of baryons and the number of protons. This calculation should account for the relative abundance of each isotope to compute the PNE ratios. However, the current calculation effectively assumes equal abundance between the components in the material.

In practice this does not change the nucleon density used for injection because these numbers are not referenced for nucleon density. However, this bug does result in approximately a 5% overestimation in the electron density, which will in turn affect the injection of Glashow resonance events.

The loop over as material's isotopes is here:

for (it2 = mats.begin(); it2 != mats.end(); ++it2) {

@austinschneider austinschneider self-assigned this Feb 28, 2021
@austinschneider austinschneider added the bug Something isn't working label Feb 28, 2021
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

1 participant