TABLE OF CONTENTS
- 1. /expNeutrinofluxes
- 1.1. expNeutrinofluxes/Eb
- 1.2. expNeutrinofluxes/Eflux_min
- 1.3. expNeutrinofluxes/Eflux_max
- 1.4. expNeutrinofluxes/T2K_oscillated
- 1.5. expNeutrinofluxes/readinput
- 1.6. expNeutrinofluxes/nl_neutrino_energyFlux
- 1.7. expNeutrinofluxes/MiniBooNEenergy
- 1.8. expNeutrinofluxes/MiniBooNEenergyBARNU
- 1.9. expNeutrinofluxes/MiniBooNE_recQs
- 1.10. expNeutrinofluxes/MiniBooNE_recEnergy
- 1.11. expNeutrinofluxes/MiniBooNE_recEnergy_Delta
- 1.12. expNeutrinofluxes/ANLenergy
- 1.13. expNeutrinofluxes/BNLenergy
- 1.14. expNeutrinofluxes/K2Kenergy
- 1.15. expNeutrinofluxes/K2K_recEnergy
- 1.16. expNeutrinofluxes/T2Kinput
- 1.17. expNeutrinofluxes/T2K_energyFlux
- 1.18. expNeutrinofluxes/T2K_OA25_energy
- 1.19. expNeutrinofluxes/MINOSenergyNU_fluxNU
- 1.20. expNeutrinofluxes/MINOSenergyBARNU_fluxNU
- 1.21. expNeutrinofluxes/MINOSenergyNU_fluxBARNU
- 1.22. expNeutrinofluxes/MINOSenergyBARNU_fluxBARNU
- 1.23. expNeutrinofluxes/NOVAenergyNU
- 1.24. expNeutrinofluxes/uniformFlux
- 1.25. expNeutrinofluxes/MINERVAenergyNU
- 1.26. expNeutrinofluxes/MINERVAenergyBARNU
- 1.27. expNeutrinofluxes/LBNEenergyNU
- 1.28. expNeutrinofluxes/LBNEenergyBARNU
- 2. /esample
- 2.1. esample/read_fluxfile
- 2.2. esample/eneut
/esample [ Modules ]
NAME
module esample
This module contains 2 files to read in the file with flux data for a given neutrino exeriment and for the sampling of neutrino energies from that flux
esample/read_fluxfile [ Functions ]
[ Top ] [ esample ] [ Functions ]
NAME
subroutine read_fluxfile
PURPOSE
Reads the input file to obtain the energies and corresponding flux. The fluxfile must contain line by line: energy (middle of bin) and flux. The file can contain comments in the first few lines, starting with # the subroutine returns the number of elements in the fluxfile, as well as the cumulative flux distribution (sum of fluxes).
SYNOPSIS
subroutine read_fluxfile(fluxfilename,n_E,enu,flux,sumflux)
INPUTS
fluxfilename: name of input flux file, with energy and flux pairwise in different lines.
RESULT
n_E: number of records in data file (not counting comments). enu: vector of neutrino energies. flux: vector of flux values. sumflux: vector with partial sums of flux (cumulative flux).
esample/eneut [ Functions ]
[ Top ] [ esample ] [ Functions ]
NAME
function eneut
PURPOSE
returns one energy value by sampling the flux distribution using discrete cumulative sampling
SYNOPSIS
eneut(n_E,sumflux,enu)
INPUTS
n_E: number of elements in input flux file and in sumflux and enu. sumflux: vector of cumulative flux distributions. enu: vector of neutrino energies.
RESULT
eneut: one sampled neutrino energy
/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/T2K_oscillated [ Global module-variables ]
[ Top ] [ expNeutrinofluxes ] [ Global module-variables ]
SOURCE
logical, save :: T2K_oscillated=.false.
PURPOSE
variables used for T2K
if true, use oscillated flux
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/MiniBooNEenergy [ Functions ]
[ Top ] [ expNeutrinofluxes ] [ Functions ]
NAME
real function MiniBooNEenergy()
PURPOSE
This function gives back 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)
expNeutrinofluxes/MiniBooNEenergyBARNU [ Functions ]
[ Top ] [ expNeutrinofluxes ] [ Functions ]
NAME
real function MiniBooNEenergyBARNU()
PURPOSE
This function gives back 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)
expNeutrinofluxes/MiniBooNE_recQs [ Functions ]
[ Top ] [ expNeutrinofluxes ] [ Functions ]
NAME
real function MiniBooNE_recQs(k_out)
PURPOSE
This function gives back 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/MiniBooNE_recEnergy [ Functions ]
[ Top ] [ expNeutrinofluxes ] [ Functions ]
NAME
real function MiniBooNE_recEnergy(k_out)
PURPOSE
This function gives back 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/MiniBooNE_recEnergy_Delta [ Functions ]
[ Top ] [ expNeutrinofluxes ] [ Functions ]
NAME
real function MiniBooNE_recEnergy_Delta(k_out)
PURPOSE
This function gives back 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 gives back the neutrino energy for the ANL experiment for QE events Flux is taken from PRD 16, 3103 (1977), Fig. 7.
expNeutrinofluxes/BNLenergy [ Functions ]
[ Top ] [ expNeutrinofluxes ] [ Functions ]
NAME
real function BNLenergy()
PURPOSE
This function gives back 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
expNeutrinofluxes/K2Kenergy [ Functions ]
[ Top ] [ expNeutrinofluxes ] [ Functions ]
NAME
real function K2Kenergy()
PURPOSE
This function gives back the neutrino energy for the K2K experiment. Flux is taken from PLB 619 (2005), Fig. 1
expNeutrinofluxes/K2K_recEnergy [ Functions ]
[ Top ] [ expNeutrinofluxes ] [ Functions ]
NAME
real function K2K_recEnergy(k_out)
PURPOSE
This function gives back 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 ArXiV 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/T2Kinput [ Subroutines ]
[ Top ] [ expNeutrinofluxes ] [ Subroutines ]
NAME
subroutine T2Kinput
PURPOSE
This subroutine reads out the namelist "T2K_energyFlux". Only called once to initialize the module.
expNeutrinofluxes/T2K_energyFlux [ Namelists ]
[ Top ] [ expNeutrinofluxes ] [ Namelists ]
NAME
NAMELIST T2K_energyFlux
PURPOSE
This Namelist includes:
expNeutrinofluxes/T2K_OA25_energy [ Functions ]
[ Top ] [ expNeutrinofluxes ] [ Functions ]
NAME
real function T2K_OA25_energy()
PURPOSE
This function gives back the neutrino energy for the T2K ND280 experiment. Flux is 2.5 degrees off-axis flux for the ND280 detector implemented is ND280_horn_205kA taken from the http://t2k-experiment.org/results/
expNeutrinofluxes/MINOSenergyNU_fluxNU [ Functions ]
[ Top ] [ expNeutrinofluxes ] [ Functions ]
NAME
real function MINOSenergyNU_fluxNU()
PURPOSE
This function gives back 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)
expNeutrinofluxes/MINOSenergyBARNU_fluxNU [ Functions ]
[ Top ] [ expNeutrinofluxes ] [ Functions ]
NAME
real function MINOSenergyBARNU_fluxNU()
PURPOSE
This function gives back 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 gives back 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 gives back the muon-antineutrino energy NUMI low-energy flux in antineutrino mode. Flux sent to us by LAr team (Ornella Palamara)
expNeutrinofluxes/NOVAenergyNU [ Functions ]
[ Top ] [ expNeutrinofluxes ] [ Functions ]
NAME
real function MINOSenergyNU()
PURPOSE
This function gives back the neutrino energy for the NOVA neutrino experiment. (NuMI medium-energy 14mrad off-axis) Flux is taken from http://www-nova.fnal.gov/nova_beam_anu.html
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 samples the antineutrino energy for the MINERvA experiment in neutrino mode. The sampling uses the discrete inversion method. Flux is obtained from B. Tice, June 2013
expNeutrinofluxes/MINERVAenergyBARNU [ Functions ]
[ Top ] [ expNeutrinofluxes ] [ Functions ]
NAME
real function MINERVAenergyBARNU()
PURPOSE
This function samples the antineutrino energy for the MINERvA experiment in antineutrino mode. The sampling uses the discrete inversion method. Flux is obtained from B. Tice, June 2013
expNeutrinofluxes/LBNEenergyNU [ Functions ]
[ Top ] [ expNeutrinofluxes ] [ Functions ]
NAME
real function LBNEenergyNU()
PURPOSE
This function samples the neutrino energy for the LBNE experiment in neutrino mode. The sampling uses the discrete inversion method. Flux is obtained from P. Huber, Virginia Tech, March 2013
expNeutrinofluxes/LBNEenergyBARNU [ Functions ]
[ Top ] [ expNeutrinofluxes ] [ Functions ]
NAME
real function LBNEenergyBARNU()
PURPOSE
This function samples the antineutrino energy for the LBNE experiment in antineutrino mode. The sampling uses the discrete inversion method. Flux is obtained from P. Huber, Virginia Tech, March 2013