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

Seltzer-Berger positron correction fix #401

Merged
merged 5 commits into from
Apr 9, 2022

Conversation

amandalund
Copy link
Contributor

@amandalund amandalund commented Apr 5, 2022

I believe we were missing a factor of two pi in the SB positron cross section correction factor -- added here. I've also:

  • Capped the number of iterations in the SB rejection loop for sampling the gamma energy at 100 (same as Geant4)
  • Used the correct maximum differential cross section for positrons in the rejection loop for sampling the gamma energy
  • Updated the geant exporter to use the relativistic pair production model instead of Bethe-Heitler and to use both the SB and relativistic brems models

- Add max number of iterations in SB rejection loop for sampling gamma energy
- Update geant exporter to use both SB and relativistic brems and relativistic pair production
@amandalund amandalund requested a review from sethrj April 5, 2022 22:04
@amandalund amandalund added bug Something isn't working physics Particles, processes, and stepping algorithms labels Apr 5, 2022
Copy link
Member

@sethrj sethrj left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

One name change suggestion. Excellent work!

src/physics/em/detail/SBEnergyDistHelper.hh Outdated Show resolved Hide resolved
@amandalund
Copy link
Contributor Author

Thanks to you for catching on to the other bug!

@sethrj sethrj merged commit 468c930 into celeritas-project:master Apr 9, 2022
@amandalund amandalund deleted the sb-fix branch April 9, 2022 20:50
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working physics Particles, processes, and stepping algorithms
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants