crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-27 05:35:49 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

Break up photon heating by nuclide

Browse source
This commit is contained in:
Paul Romano 2022-06-02 14:48:16 -05:00
parent ce035884ac
commit daf7bc4955
1 changed files with 15 additions and 10 deletions
Download patch file Download diff file Expand all files Collapse all files

25
src/tallies/tally_scoring.cpp
View file

@ -963,18 +963,23 @@ void score_general_ce_nonanalog(Particle& p, int i_tally, int start_index,
score = score_neutron_heating(
p, tally, flux, HEATING, i_nuclide, atom_density);
} else {
// The energy deposited is the difference between the pre-collision and
// post-collision energy...
score = E - p.E();
if (i_nuclide < 0 || i_nuclide == p.event_nuclide()) {
// The energy deposited is the difference between the pre-collision
// and post-collision energy...
score = E - p.E();
// ...less the energy of any secondary particles since they will be
// transported individually later
const auto& bank = p.secondary_bank();
for (auto it = bank.end() - p.n_bank_second(); it < bank.end(); ++it) {
score -= it->E;
// ...less the energy of any secondary particles since they will be
// transported individually later
const auto& bank = p.secondary_bank();
for (auto it = bank.end() - p.n_bank_second(); it < bank.end();
++it) {
score -= it->E;
}
score *= p.wgt_last();
} else {
score = 0.0;
}
score *= p.wgt_last();
}
break;