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
 All Classes Functions Variables Enumerations Pages
Public Member Functions | List of all members
marley::DecayAtRestNeutrinoSource Class Reference

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

#include <NeutrinoSource.hh>

Inheritance diagram for marley::DecayAtRestNeutrinoSource:
marley::NeutrinoSource

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...
 

Additional Inherited Members

- 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_idPDG 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
Eneutrino 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: