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

Fix check for no angle distribution in elastic scattering

Browse source
This commit is contained in:
Paul Romano 2020-01-14 13:34:17 -06:00
parent d0d1b1db12
commit f417f35809
1 changed files with 3 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

6
src/physics.cpp
View file

@ -696,12 +696,12 @@ void elastic_scatter(int i_nuclide, const Reaction& rx, double kT,
// Find speed of neutron in CM
vel = v_n.norm();
// Sample scattering angle, checking if it is an ncorrelated angle-energy
// distribution
// Sample scattering angle, checking if angle distribution is present (assume
// isotropic otherwise)
double mu_cm;
auto& d = rx.products_[0].distribution_[0];
auto d_ = dynamic_cast<UncorrelatedAngleEnergy*>(d.get());
if (d_) {
if (!d_->angle().empty()) {
mu_cm = d_->angle().sample(p->E_, p->current_seed());
} else {
mu_cm = 2.0*prn(p->current_seed()) - 1.0;