"Beta-fit" neutrino source More...

#include <NeutrinoSource.hh>

Inheritance diagram for marley::BetaFitNeutrinoSource:

Public Member Functions
	BetaFitNeutrinoSource (int particle_id=marley_utils::ELECTRON_NEUTRINO, double Emin=0., double Emax=50., double Emean=13., double beta=4.5)

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

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) const
	Samples an incident neutrino energy and loads pdg with the PDG code of the appropriate neutrino type. More...

Protected Attributes
double	beta_
	pinching parameter

double	C_
	dimensionless normalization constant

double	Emax_

double	Emean_
	mean neutrino energy More...

double	Emin_
	minimum neutrino energy (MeV)

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

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

Detailed Description

"Beta-fit" neutrino source

Neutrino energies from this source are sampled from a "beta-fit" spectrum (see, for example, equation 7 in this preprint) with mean neutrino energy \(E_\text{mean}\) and fit parameter \(\beta\). The probability density function is given by

\begin{align*} P(E) &= C\left(E/\,E_\text{mean}\right)^{\beta - 1} \exp\left(-\beta\,E/\,E_\text{mean}\right) & \text{E}_\text{min} \leq E \leq \text{E}_\text{max} \end{align*}

where \(C\) is a normalization constant.

Constructor & Destructor Documentation

◆ BetaFitNeutrinoSource()

marley::BetaFitNeutrinoSource::BetaFitNeutrinoSource	(	int	particle_id = `marley_utils::ELECTRON_NEUTRINO`,
		double	Emin = `0.`,
		double	Emax = `50.`,
		double	Emean = `13.`,
		double	beta = `4.5`
	)

Parameters

particle_id	PDG particle ID for the neutrinos produced by this source
Emin	minimum allowed neutrino energy (MeV)
Emax	maximum allowed neutrino energy (MeV)
Emean	mean neutrino energy (MeV)
beta	dimensionless fit parameter

Member Function Documentation

◆ pdf()

double marley::BetaFitNeutrinoSource::pdf ( double E ) const

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)

Implements marley::NeutrinoSource.

Member Data Documentation

◆ Emax_

double marley::BetaFitNeutrinoSource::Emax_

protected

maximum neutrino energy (MeV)

◆ Emean_

double marley::BetaFitNeutrinoSource::Emean_

protected

mean neutrino energy

Note: This is exactly the mean neutrino energy only if \( E_\text{min} = 0 \) and \( Emax = \infty \). Truncating the distribution may alter the mean value appreciably.

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

/Users/gardiner/marley/include/marley/NeutrinoSource.hh
/Users/gardiner/marley/src/NeutrinoSource.cc

Public Member Functions

Protected Attributes

Additional Inherited Members