TABLE OF CONTENTS
- 1. /expNeutrinofluxes
- 1.1. expNeutrinofluxes/Eb
- 1.2. expNeutrinofluxes/Eflux_min
- 1.3. expNeutrinofluxes/Eflux_max
- 1.4. expNeutrinofluxes/readinput
- 1.5. expNeutrinofluxes/nl_neutrino_energyFlux
- 1.6. expNeutrinofluxes/getFluxEnu
- 1.7. expNeutrinofluxes/MiniBooNEenergy
- 1.8. expNeutrinofluxes/MiniBooNEenergyBARNU
- 1.9. expNeutrinofluxes/CCQE_recQs
- 1.10. expNeutrinofluxes/CCQE_recEnergy
- 1.11. expNeutrinofluxes/CCQE_recEnergy_Delta
- 1.12. expNeutrinofluxes/ANLenergy
- 1.13. expNeutrinofluxes/BNLenergy
- 1.14. expNeutrinofluxes/K2Kenergy
- 1.15. expNeutrinofluxes/K2K_recEnergy
- 1.16. expNeutrinofluxes/T2Kenergy_ND
- 1.17. expNeutrinofluxes/MINOSenergyNU_fluxNU
- 1.18. expNeutrinofluxes/MINOSenergyBARNU_fluxNU
- 1.19. expNeutrinofluxes/MINOSenergyNU_fluxBARNU
- 1.20. expNeutrinofluxes/MINOSenergyBARNU_fluxBARNU
- 1.21. expNeutrinofluxes/uniformFlux
- 1.22. expNeutrinofluxes/MINERVAenergyNU
- 1.23. expNeutrinofluxes/MINERVAenergyBARNU
- 1.24. expNeutrinofluxes/DUNEenergyNU
- 1.25. expNeutrinofluxes/DUNEenergyBARNU
- 1.26. expNeutrinofluxes/LBNOenergyNU
- 1.27. expNeutrinofluxes/NOMADenergyNU
- 1.28. expNeutrinofluxes/BNBenergyNUe
- 1.29. expNeutrinofluxes/BNBenergyNUebar
- 1.30. expNeutrinofluxes/BNBenergyNUmu
- 1.31. expNeutrinofluxes/BNBenergyNUmubar
- 1.32. expNeutrinofluxes/NOVAenergy_ND
- 1.33. expNeutrinofluxes/NOVAenergyNU_FD
- 1.34. expNeutrinofluxes/MINERvAenergy
- 1.35. expNeutrinofluxes/T2Konaxisenergy
- 1.36. expNeutrinofluxes/FASERenergy
- 1.37. expNeutrinofluxes/userFlux
- 1.38. expNeutrinofluxes/userFluxEmax
/expNeutrinofluxes [ 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