mirror of
https://github.com/openmc-dev/openmc.git
synced 2026-07-28 14:15:42 -04:00
61 lines
2.1 KiB
C++
61 lines
2.1 KiB
C++
//! \file secondary_correlated.h
|
|
//! Correlated angle-energy distribution
|
|
|
|
#ifndef OPENMC_SECONDARY_CORRELATED_H
|
|
#define OPENMC_SECONDARY_CORRELATED_H
|
|
|
|
#include <vector>
|
|
|
|
#include "hdf5.h"
|
|
#include "xtensor/xtensor.hpp"
|
|
|
|
#include "openmc/angle_energy.h"
|
|
#include "openmc/endf.h"
|
|
#include "openmc/distribution.h"
|
|
|
|
namespace openmc {
|
|
|
|
//==============================================================================
|
|
//! Correlated angle-energy distribution corresponding to ACE law 61 and ENDF
|
|
//! File 6, LAW=1, LANG!=2.
|
|
//==============================================================================
|
|
|
|
class CorrelatedAngleEnergy : public AngleEnergy {
|
|
public:
|
|
//! Outgoing energy/angle at a single incoming energy
|
|
struct CorrTable {
|
|
int n_discrete; //!< Number of discrete lines
|
|
Interpolation interpolation; //!< Interpolation law
|
|
xt::xtensor<double, 1> e_out; //!< Outgoing energies [eV]
|
|
xt::xtensor<double, 1> p; //!< Probability density
|
|
xt::xtensor<double, 1> c; //!< Cumulative distribution
|
|
std::vector<UPtrDist> angle; //!< Angle distribution
|
|
};
|
|
|
|
explicit CorrelatedAngleEnergy(hid_t group);
|
|
|
|
//! Sample distribution for an angle and energy
|
|
//! \param[in] E_in Incoming energy in [eV]
|
|
//! \param[out] E_out Outgoing energy in [eV]
|
|
//! \param[out] mu Outgoing cosine with respect to current direction
|
|
void sample(double E_in, double& E_out, double& mu) const override;
|
|
|
|
// energy property
|
|
std::vector<double>& energy() { return energy_; }
|
|
const std::vector<double>& energy() const { return energy_; }
|
|
|
|
// distribution property
|
|
std::vector<CorrTable>& distribution() { return distribution_; }
|
|
const std::vector<CorrTable>& distribution() const { return distribution_; }
|
|
private:
|
|
int n_region_; //!< Number of interpolation regions
|
|
std::vector<int> breakpoints_; //!< Breakpoints between regions
|
|
std::vector<Interpolation> interpolation_; //!< Interpolation laws
|
|
std::vector<double> energy_; //!< Energies [eV] at which distributions
|
|
//!< are tabulated
|
|
std::vector<CorrTable> distribution_; //!< Distribution at each energy
|
|
};
|
|
|
|
} // namespace openmc
|
|
|
|
#endif // OPENMC_SECONDARY_CORRELATED_H
|