gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS


/expNeutrinofluxes [ Modules ]

[ Top ] [ Modules ]

NAME

module expNeutrinofluxes

PURPOSE

This module provides specific experimental neutrino fluxes and it selects the neutrino energy according to the experimental flux. For MiniBooNE and K2K, it also extracts the reconstructed neutrino energy and Qs as it is done in the experiment.


expNeutrinofluxes/Eb [ Global module-variables ]

[ Top ] [ expNeutrinofluxes ] [ Global module-variables ]

SOURCE

  real, save    :: Eb=0.034

PURPOSE

contant binding energy used for energy and Q2 reconstruction based on QE scattering kinematics


expNeutrinofluxes/Eflux_min [ Global module-variables ]

[ Top ] [ expNeutrinofluxes ] [ Global module-variables ]

SOURCE

  real, save    :: Eflux_min=0.2

PURPOSE

minimum energy for uniform flux distribution

minimum and maximum energies for the uniform neutrino flux (nuExp=10 in the namelist neutrino_induced) can be changed in the namelist nl_neutrino_energyFlux


expNeutrinofluxes/Eflux_max [ Global module-variables ]

[ Top ] [ expNeutrinofluxes ] [ Global module-variables ]

SOURCE

  real, save    :: Eflux_max=2.5

PURPOSE

maximum energy for uniform flux distribution

minimum and maximum energies for the uniform neutrino flux (nuExp=10 in the namelist neutrino_induced) can be changed in the namelist nl_neutrino_energyFlux


expNeutrinofluxes/readinput [ Subroutines ]

[ Top ] [ expNeutrinofluxes ] [ Subroutines ]

NAME

subroutine readinput

INPUTS

NONE

OUTPUT

NONE

PURPOSE

This subroutine reads out the namelist "nl_neutrino_energyFlux". Only called once to initialize the module.


expNeutrinofluxes/nl_neutrino_energyFlux [ Namelists ]

[ Top ] [ expNeutrinofluxes ] [ Namelists ]

NAME

NAMELIST nl_neutrino_energyFlux

PURPOSE

This Namelist includes:


expNeutrinofluxes/getFluxEnu [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function getFluxEnu(nuExp,iFlavor,iProcess,FileNameFlux)

PURPOSE

Return the neutrino energy randomly sampled according the given flux.

This is a wrapper around the different routines for the different experiments.

INPUTS

  • integer :: nuExp -- number of experiment
  • integer :: iFlavor -- kind of beam (e,mu,tau)
  • integer :: iProcess -- kind of process (EM,CC,NC)
  • character*(*) :: FileNameFlux -- filename for user defined fluxes

OUTPUT

function value :: neutrino energy (GeV)


expNeutrinofluxes/MiniBooNEenergy [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function MiniBooNEenergy()

PURPOSE

This function returns the neutrino energy for the MiniBooNE experiment. It determines the energy randomly weighted with the flux. Flux is taken from http://www-boone.fnal.gov/for_physicists/data_release/flux/ and normalized to 1

paper for reference:

  • A.A.Aguilar-Arevalo et al., "The Neutrino Flux Prediction at MiniBooNE" Phys. Rev. D. 79, 072002 (2009)

NOTES

  • this routine is outdated


expNeutrinofluxes/MiniBooNEenergyBARNU [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function MiniBooNEenergyBARNU()

PURPOSE

This function returns the antineutrino energy for the MiniBooNE experiment in antineutrino mode (=negartive polarity). It determines the energy randomly weighted with the flux. Flux is taken from http://www-boone.fnal.gov/for_physicists/data_release/flux/

paper for reference:

  • A.A.Aguilar-Arevalo et al., "The Neutrino Flux Prediction at MiniBooNE" Phys. Rev. D. 79, 072002 (2009)

NOTES

  • this routine is outdated


expNeutrinofluxes/CCQE_recQs [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function CCQE_recQs(k_out)

PURPOSE

This function returns the reconstruced Qs. The reconstruction is done as in the experiment, neglecting Fermi motion (see arXiv:0706.0926v1, eq.(3)), where k_out is the "real" outgoing lepton momentum.


expNeutrinofluxes/CCQE_recEnergy [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function CCQE_recEnergy(k_out)

PURPOSE

This function returns the reconstruced neutrino energy. The reconstruction is done as in the experiment, neglecting Fermi motion (see arXiv:0706.0926v1, eq.(4)), where k_out is the "real" outgoing lepton momentum.


expNeutrinofluxes/CCQE_recEnergy_Delta [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function CCQE_recEnergy_Delta(k_out)

PURPOSE

This function returns the reconstruced neutrino energy for pions. The reconstruction is done as in the experiment, neglecting Fermi motion and binding (see PRL 103, 081801 (2009) eq.(1)), where k_out is the "real" outgoing lepton momentum.


expNeutrinofluxes/ANLenergy [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function ANLenergy()

PURPOSE

This function returns the neutrino energy for the ANL experiment for QE events

Flux is taken from PRD 16, 3103 (1977), Fig. 7.

NOTES

  • this routine is outdated


expNeutrinofluxes/BNLenergy [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function BNLenergy()

PURPOSE

This function returns the neutrino energy for the BNL experiment. Flux is taken from K. Furuno, NUINT02 proceedings, available at http://www.ps.uci.edu/~nuint/proceedings/furuno.pdf or see Baker et al Phys Rev D23 (1981) 2499, fig.7

NOTES

  • enumin is for the whole flux and is good for calculating event histograms; for calculating absolute cross sections BNL used enumin=0.5, which should be used here
  • this routine is outdated


expNeutrinofluxes/K2Kenergy [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function K2Kenergy()

PURPOSE

This function returns the neutrino energy for the K2K experiment. Flux is taken from PLB 619 (2005), Fig. 1

NOTES

  • this routine is outdated


expNeutrinofluxes/K2K_recEnergy [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function K2K_recEnergy(k_out)

PURPOSE

This function returns the reconstruced neutrino energy. The reconstruction is done as in the experiment, neglecting Fermi motion (see PRL 90, 041801 (2003), eq.(1)), where k_out is the "real" outgoing lepton momentum.

another formula is used in recent CC-pi0/QE measurements, see https://arxiv.org/abs/1012.1794

it will be used if an optional parameter W is given (default K2K W=1.483 GeV) for W=MN the formular coincides with the old one


expNeutrinofluxes/T2Kenergy_ND [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function T2Kenergy_ND(iFlavor,iProcess)

PURPOSE

This function returns the neutrino energy for the T2K ND280 experiment. Flux is 2.5 degrees off-axis flux for the ND280 detector implemented is ND280_horn_250kA taken from http://t2k-experiment.org/results/


expNeutrinofluxes/MINOSenergyNU_fluxNU [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function MINOSenergyNU_fluxNU()

PURPOSE

This function returns the muon-neutrino energy for the MINOS neutrino experiment (NUMI low-energy flux) in neutrino mode. Flux sent to us by Minerva team (Steve Dytman)

NOTES

  • this routine is outdated


expNeutrinofluxes/MINOSenergyBARNU_fluxNU [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function MINOSenergyBARNU_fluxNU()

PURPOSE

This function returns the muon-antineutrino energy for the MINOS neutrino experiment (NUMI low-energy flux) in neutrino mode. Flux sent to us by LAr team (Ornella Palamara)


expNeutrinofluxes/MINOSenergyNU_fluxBARNU [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function MINOSenergyNU_fluxBARNU()

PURPOSE

This function returns the muon-neutrino energy for NUMI low-energy flux in antineutrino mode. Flux sent to us by LAr team (Ornella Palamara)


expNeutrinofluxes/MINOSenergyBARNU_fluxBARNU [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function MINOSenergyBARNU_fluxBARNU()

PURPOSE

This function returns the muon-antineutrino energy NUMI low-energy flux in antineutrino mode. Flux sent to us by LAr team (Ornella Palamara)


expNeutrinofluxes/uniformFlux [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function uniformFlux()

PURPOSE

generated uniform flux from Eflux_min to Eflux_max (see namelist nl_neturino_energyFlux)


expNeutrinofluxes/MINERVAenergyNU [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function MINERVAenergyNU()

PURPOSE

This function returns the sampled neutrino energy for the MINERvA experiment in neutrino mode. Flux is obtained from B. Tice, June 2013, now outdated


expNeutrinofluxes/MINERVAenergyBARNU [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function MINERVAenergyBARNU()

PURPOSE

This function returns the sampled antineutrino energy for the MINERvA experiment in antineutrino mode. Flux is obtained from B. Tice, June 2013, now outdated


expNeutrinofluxes/DUNEenergyNU [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function DUNEenergyNU(iFlavor,iProcess)

PURPOSE

This function returns the sampled neutrino energies for the DUNE experiment in neutrino mode for both mu- and e-neutrinos Flux is obtained from http://home.fnal.gov/~ljf26/DUNEFluxes/ "Engineered Optimized 3 Horn Design"


expNeutrinofluxes/DUNEenergyBARNU [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function DUNEenergyBARNU(iFlavor,iProcess)

PURPOSE

This function returns the sampled neutrino energies for the DUNE experiment in antineutrino mode for both mu- and e-neutrinos Flux is obtained from http://home.fnal.gov/~ljf26/DUNEFluxes/ "Engineered Optimized 3 Horn Design"


expNeutrinofluxes/LBNOenergyNU [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function LBNOenergyNU()

PURPOSE

This function returns the sampled neutrino energy for the LBNO experiment in antineutrino mode.

Flux is obtained from Dario Autiero, May 2014


expNeutrinofluxes/NOMADenergyNU [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function NOMADenergyNU()

PURPOSE

This function returns the sampled neutrino energy for the NOMAD experiment in antineutrino mode.

Flux is obtained from Roberto Petti, Jan. 2014


expNeutrinofluxes/BNBenergyNUe [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function BNBenergyNUe()

PURPOSE

This function returns the sampled neutrino energy for the BNB experiment in electron-neutrino mode.

Flux is obtained from Anne Schukraft, December 2015


expNeutrinofluxes/BNBenergyNUebar [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function BNBenergyNUebar()

PURPOSE

This function returns the sampled neutrino energy for the BNB experiment in electron-antineutrino mode.

Flux is obtained from Anne Schukraft, December 2015


expNeutrinofluxes/BNBenergyNUmu [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function BNBenergyNUmu()

PURPOSE

This function returns the sampled neutrino energy for the BNB experiment in muon-neutrino mode.

Flux is obtained from Anne Schukraft, December 2015


expNeutrinofluxes/BNBenergyNUmubar [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function BNBenergyNUmubar()

PURPOSE

This function returns the sampled neutrino energy for the BNB experiment in muon-antineutrino mode.

Flux is obtained from Anne Schukraft, December 2015


expNeutrinofluxes/NOVAenergy_ND [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function NOVAenergy_ND(iFlavor,iProcess)

PURPOSE

This function returns the sampled neutrino energy for the NOvA experiment at the Near Detector, using the FHC files for neutrinos and RHC files for antineutrinos

Flux is obtained from Jonathan Paley, March 2016


expNeutrinofluxes/NOVAenergyNU_FD [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function NOVAenergyNU_FD(iFlavor,iProcess)

PURPOSE

This function returns the sampled neutrino energy for the NOvA experiment at the Far Detector, using the FHC files for neutrinos and RHC files for antineutrinos Flux is obtained from Jonathan Paley, March 2016


expNeutrinofluxes/MINERvAenergy [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function MINERvAenergy(iFlavor,iProcess)

PURPOSE

This function returns the sampled neutrino energy for the MINERvA experiment at the Far Detector, flux is obtained from Phys.Rev. D94 (2016) no.9, 092005, Addendum: Phys.Rev. D95 (2017) no.3, 039903


expNeutrinofluxes/T2Konaxisenergy [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function T2Konaxisenergy(iFlavor,iProcess)

PURPOSE

This function returns the sampled neutrino energy for the T2K on axis experiment


expNeutrinofluxes/FASERenergy [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function FASERenergy(iFlavor,iProcess)

PURPOSE

This function returns the sampled neutrino energy for the FASERnu experiment at CERN, flux is obtained from from simulation first presented in 2105.08270 generated with EPOS-LHC and POWHEG+Pythia8, updated to Run 4 and with updated charm production from 2309.12793" , Max Fieg, priv. comm, Jan 2024, these flux files have been rebinned


expNeutrinofluxes/userFlux [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function userFlux(FileNameFlux)

PURPOSE

This function returns the sampled neutrino energy according the user provided file.

This file must be a .txt file with two columns: 1st col: energy, 2nd col: flux, the energies must be equidistant


expNeutrinofluxes/userFluxEmax [ Functions ]

[ Top ] [ expNeutrinofluxes ] [ Functions ]

NAME

real function userFluxEMax

PURPOSE

This function returns the maximal energy value of the user provided flux