MARLEY (Model of Argon Reaction Low Energy Yields)
v1.2.0
A Monte Carlo event generator for tens-of-MeV neutrino interactions
|
Muon decay-at-rest neutrino source. More...
#include <NeutrinoSource.hh>
Public Member Functions | |
DecayAtRestNeutrinoSource (int particle_id=marley_utils::ELECTRON_NEUTRINO) | |
virtual double | get_Emax () const override |
Get the maximum neutrino energy (MeV) that can be sampled by this source. | |
virtual double | get_Emin () const override |
Get the minimum neutrino energy (MeV) that can be sampled by this source. | |
virtual double | pdf (double E) const override |
Probability density function describing the incident neutrino energy distribution. More... | |
![]() | |
NeutrinoSource (int particle_id) | |
virtual int | get_pid () const |
Get the PDG particle ID for the neutrino type produced by this source. | |
virtual double | sample_incident_neutrino (int &pdg, marley::Generator &gen) const |
Samples an incident neutrino energy and loads pdg with the PDG code of the appropriate neutrino type. More... | |
Additional Inherited Members | |
![]() | |
static bool | pdg_is_allowed (const int pdg) |
![]() | |
int | pid_ |
PDG particle ID for the neutrinos produced by this source. | |
Muon decay-at-rest neutrino source.
Neutrino energies from this source are sampled from the appropriate Michel spectrum for muon decay-at-rest. For a \(\nu_e\) source, the probability density function is given by
\begin{align*} P(E) &= 96E^2m_\mu^{-4}(m_\mu - 2E) & 0 < E < m_\mu/2 \end{align*}
where \(m_\mu\) is the muon mass. For a \(\bar{\nu}_\mu\) source, the probability density function is given by
\begin{align*} P(E) &= 16E^2m_\mu^{-4}(3m_\mu - 4E) & 0 < E < m_\mu/2. \end{align*}
marley::DecayAtRestNeutrinoSource::DecayAtRestNeutrinoSource | ( | int | particle_id = marley_utils::ELECTRON_NEUTRINO | ) |
particle_id | PDG particle ID for the neutrinos produced by this source |
|
overridevirtual |
Probability density function describing the incident neutrino energy distribution.
The neutrino spectrum produced by this source will be folded with the relevant cross sections by a Generator object during event creation
E | neutrino energy (MeV) |
Implements marley::NeutrinoSource.