MARLEY (Model of Argon Reaction Low Energy Yields)  v1.1.0 A Monte Carlo event generator for tens-of-MeV neutrino-nucleus interactions in liquid argon
marley::DecayAtRestNeutrinoSource Class Reference

Muon decay-at-rest neutrino source. More...

#include <NeutrinoSource.hh>

Inheritance diagram for marley::DecayAtRestNeutrinoSource:

## 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) override
Probability density function describing the incident neutrino energy distribution. More...

Public Member Functions inherited from marley::NeutrinoSource
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)
Samples an incident neutrino energy and loads pdg with the PDG code of the appropriate neutrino type. More...

Static Public Member Functions inherited from marley::NeutrinoSource
static bool pdg_is_allowed (const int pdg)

Protected Attributes inherited from marley::NeutrinoSource
int pid_
PDG particle ID for the neutrinos produced by this source.

## Detailed Description

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*}

Other neutrino species are not allowed.

## Constructor & Destructor Documentation

 marley::DecayAtRestNeutrinoSource::DecayAtRestNeutrinoSource ( int particle_id = marley_utils::ELECTRON_NEUTRINO )
Parameters
 particle_id PDG particle ID for the neutrinos produced by this source

## Member Function Documentation

 double marley::DecayAtRestNeutrinoSource::pdf ( double E )
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

Parameters
 E neutrino energy (MeV)
Returns
Probability density (MeV -1)
Todo:
Refine the approximate decay-at-rest Michel spectra to use more exact expressions.

Implements marley::NeutrinoSource.

The documentation for this class was generated from the following files: