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

When nbr_split > 1 and Russian Roulette is on, split fluorescent photons resulting from electron impact ionization during Moller. #1230

Open
2 tasks
blakewalters opened this issue Dec 7, 2024 · 0 comments
Assignees

Comments

@blakewalters
Copy link
Contributor

Describe the bug
Currently, when uniform bremsstrahlung splitting (UBS) is turned on with Russian Roulette, we do not split fluorescent photons resulting from atomic relaxations induced by electron impact ionization (EII) during Moller scattering of phat electrons. This results in the emission of phat photons.

To Reproduce

  1. Use BEAMnrc to simulate a low energy (e.g., 225 keV) electron beam incident on a tungsten target using uniform brem splitting (UBS) with a splitting number (nbr_split) of, say, 100, and with Russian Roulette on. Ensure that atomic relaxations and EII are both ON.
  2. Score phase space below the target
  3. Run many histories (>= 10e6)
  4. Observe photons with weight=1 in the phase space data

Expected behavior
Scored photons should all have weight 1/nbr_split. Scored charged particles should all have weight 1.

Operating system

  • [x ] Linux
  • Windows
  • [ x] macOS
  • Other:

EGSnrc version
develop branch

Other
This showed up when testing the egs++ DBS algorithm against BEAMnrc. Using BEAMnrc with UBS and nbr_split=100, I found that when I used a photon weight cutoff of 0.01 (i.e. eliminating phat photons), the energy spectrum with Russian Roulette on was up to ~40% lower than that with Russian Roulette off for energies < 40 keV.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

1 participant