diff --git a/src/secondary_correlated.cpp b/src/secondary_correlated.cpp index 8e0f23d520..476d783609 100644 --- a/src/secondary_correlated.cpp +++ b/src/secondary_correlated.cpp @@ -259,7 +259,7 @@ void CorrelatedAngleEnergy::sample(double E_in, double& E_out, double& mu, } // Find correlated angular distribution for closest outgoing energy bin - if (r1 - c_k < c_k1 - r1) { + if (r1 - c_k < c_k1 - r1 || distribution_[l].interpolation == Interpolation::histogram) { mu = distribution_[l].angle[k]->sample(seed); } else { mu = distribution_[l].angle[k + 1]->sample(seed);