crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 14:15:42 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Make sure photons don't contribute to many tally scores

Browse source
This commit is contained in:
Paul Romano 2020-07-07 08:57:56 -05:00
parent 7594e7570c
commit e0e2891f2a
1 changed files with 13 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

15
src/tallies/tally_scoring.cpp
View file

@ -588,6 +588,8 @@ score_general_ce(Particle& p, int i_tally, int start_index, int filter_index,
case SCORE_INVERSE_VELOCITY:
if (p.type_ != Particle::Type::neutron) continue;
if (tally.estimator_ == TallyEstimator::ANALOG) {
// All events score to an inverse velocity bin. We actually use a
// collision estimator in place of an analog one since there is no way
@ -609,6 +611,8 @@ score_general_ce(Particle& p, int i_tally, int start_index, int filter_index,
case SCORE_SCATTER:
if (p.type_ != Particle::Type::neutron) continue;
if (tally.estimator_ == TallyEstimator::ANALOG) {
// Skip any event where the particle didn't scatter
if (p.event_ != TallyEvent::SCATTER) continue;
@ -620,14 +624,15 @@ score_general_ce(Particle& p, int i_tally, int start_index, int filter_index,
score = (p.neutron_xs_[i_nuclide].total
- p.neutron_xs_[i_nuclide].absorption) * atom_density * flux;
} else {
score = (p.macro_xs_.total
- p.macro_xs_.absorption) * flux;
score = (p.macro_xs_.total - p.macro_xs_.absorption) * flux;
}
}
break;
case SCORE_NU_SCATTER:
if (p.type_ != Particle::Type::neutron) continue;
// Only analog estimators are available.
// Skip any event where the particle didn't scatter
if (p.event_ != TallyEvent::SCATTER) continue;
@ -649,6 +654,8 @@ score_general_ce(Particle& p, int i_tally, int start_index, int filter_index,
case SCORE_ABSORPTION:
if (p.type_ != Particle::Type::neutron) continue;
if (tally.estimator_ == TallyEstimator::ANALOG) {
if (settings::survival_biasing) {
// No absorption events actually occur if survival biasing is on --
@ -1205,6 +1212,8 @@ score_general_ce(Particle& p, int i_tally, int start_index, int filter_index,
case ELASTIC:
if (p.type_ != Particle::Type::neutron) continue;
if (tally.estimator_ == TallyEstimator::ANALOG) {
// Check if event MT matches
if (p.event_mt_ != ELASTIC) continue;
@ -1244,6 +1253,8 @@ score_general_ce(Particle& p, int i_tally, int start_index, int filter_index,
case N_GAMMA:
case N_P:
case N_A:
if (p.type_ != Particle::Type::neutron) continue;
if (tally.estimator_ == TallyEstimator::ANALOG) {
// Check if the event MT matches
if (p.event_mt_ != score_bin) continue;