mirror of
https://github.com/openmc-dev/openmc.git
synced 2026-07-21 06:25:30 -04:00
1578698129
Co-authored-by: Your Name <you@example.com> Co-authored-by: GuySten <62616591+GuySten@users.noreply.github.com> Co-authored-by: GuySten <guyste@post.bgu.ac.il> Co-authored-by: Eliezer214 <110336440+Eliezer214@users.noreply.github.com> Co-authored-by: Paul Romano <paul.k.romano@gmail.com>
78 lines
3.1 KiB
C++
78 lines
3.1 KiB
C++
//! \file reaction_product.h
|
|
//! Data for a reaction product
|
|
|
|
#ifndef OPENMC_REACTION_PRODUCT_H
|
|
#define OPENMC_REACTION_PRODUCT_H
|
|
|
|
#include "hdf5.h"
|
|
|
|
#include "openmc/angle_energy.h"
|
|
#include "openmc/chain.h"
|
|
#include "openmc/endf.h"
|
|
#include "openmc/memory.h" // for unique_ptr
|
|
#include "openmc/particle_type.h"
|
|
#include "openmc/vector.h" // for vector
|
|
|
|
namespace openmc {
|
|
|
|
//==============================================================================
|
|
//! Data for a reaction product including its yield and angle-energy
|
|
//! distributions, each of which has a given probability of occurring for a
|
|
//! given incoming energy. In general, most products only have one angle-energy
|
|
//! distribution, but for some cases (e.g., (n,2n) in certain nuclides) multiple
|
|
//! distinct distributions exist.
|
|
//==============================================================================
|
|
|
|
class ReactionProduct {
|
|
public:
|
|
//! Emission mode for product
|
|
enum class EmissionMode {
|
|
prompt, // Prompt emission of secondary particle
|
|
delayed, // Yield represents total emission (prompt + delayed)
|
|
total // Delayed emission of secondary particle
|
|
};
|
|
|
|
using Secondary = unique_ptr<AngleEnergy>;
|
|
|
|
//! Construct reaction product from HDF5 data
|
|
//! \param[in] group HDF5 group containing data
|
|
explicit ReactionProduct(hid_t group);
|
|
|
|
//! Construct reaction product for decay photon from chain nuclide product
|
|
//! \param[in] product Chain nuclide product
|
|
explicit ReactionProduct(const ChainNuclide::Product& product);
|
|
|
|
//! Sample an outgoing 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
|
|
//! \param[inout] seed Pseudorandom seed pointer
|
|
void sample(double E_in, double& E_out, double& mu, uint64_t* seed) const;
|
|
|
|
//! Select which angle-energy distribution to sample
|
|
//! \param[in] E_in Incoming energy in [eV]
|
|
//! \param[inout] seed Pseudorandom seed pointer
|
|
//! \return Reference to the selected angle-energy distribution
|
|
AngleEnergy& sample_dist(double E_in, uint64_t* seed) const;
|
|
|
|
//! Sample an outgoing energy and evaluate the angular PDF
|
|
//! \param[in] E_in Incoming energy in [eV]
|
|
//! \param[in] mu Scattering cosine with respect to current direction
|
|
//! \param[out] E_out Outgoing energy in [eV]
|
|
//! \param[inout] seed Pseudorandom seed pointer
|
|
//! \return Probability density for the scattering cosine
|
|
double sample_energy_and_pdf(
|
|
double E_in, double mu, double& E_out, uint64_t* seed) const;
|
|
|
|
ParticleType particle_; //!< Particle type
|
|
EmissionMode emission_mode_; //!< Emission mode
|
|
double decay_rate_; //!< Decay rate (for delayed neutron precursors) in [1/s]
|
|
unique_ptr<Function1D> yield_; //!< Yield as a function of energy
|
|
vector<Tabulated1D> applicability_; //!< Applicability of distribution
|
|
vector<Secondary> distribution_; //!< Secondary angle-energy distribution
|
|
int parent_nuclide_ = -1; //!< Index of chain nuclide that is parent
|
|
};
|
|
|
|
} // namespace openmc
|
|
|
|
#endif // OPENMC_REACTION_PRODUCT_H
|