Skip to content

Commit

Permalink
Add a brief description for the operator()
Browse files Browse the repository at this point in the history
  • Loading branch information
whokion committed Mar 6, 2024
1 parent e3093fb commit c0cd4e9
Showing 1 changed file with 13 additions and 0 deletions.
13 changes: 13 additions & 0 deletions src/celeritas/neutron/interactor/ChipsNeutronElasticInteractor.hh
Original file line number Diff line number Diff line change
Expand Up @@ -103,6 +103,19 @@ CELER_FUNCTION ChipsNeutronElasticInteractor::ChipsNeutronElasticInteractor(
//---------------------------------------------------------------------------//
/*!
* Sample the final state of the neutron-nucleus elastic scattering.
*
* The scattering angle (\f$ cos \theta \f$) in the elastic neutron-nucleus
* scattering is expressed in terms of the momentum transfer (\f$ Q^{2} \f$),
* \f[
* cos \theta = 1 - \frac{Q^{2}}{2 |\vec{k}_i|^{2}}
* \f]
* where \f$ \vec{k}_i \f$ is the momentum of the incident neutron in the
* center of mass frame and the momentum transfer (\f$ Q^{2} \f$) is calculated
* according to the CHIPS (Chiral Invariant Phase Space) model (see references
* in model/ChipsNeutronElasticModel.cc and detail/MomentumTransferSampler.hh).
* The final direction of the scattered neutron in the laboratory frame is
* then transformed by the Lorentz boost of the initial four vector of the
* neutron-nucleus system.
*/
template<class Engine>
CELER_FUNCTION Interaction ChipsNeutronElasticInteractor::operator()(Engine& rng)
Expand Down

0 comments on commit c0cd4e9

Please sign in to comment.