TABLE OF CONTENTS
- 1. /esample
- 1.1. esample/Enu_upper_cut
- 1.2. esample/Enu_lower_cut
- 1.3. esample/energylimit_for_Qsrec
- 1.4. esample/read_fluxfile
- 1.5. esample/eneut
/esample [ Modules ]
NAME
module esample
PURPOSE
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/Enu_upper_cut [ Global module-variables ]
[ Top ] [ esample ] [ Global module-variables ]
SOURCE
real, save, public :: Enu_upper_cut=200.
PURPOSE
cut events with neutrino energy above Enu_upper_cut; for ANL experiment, for example, Enu_upper_cut=1.5 for ppi0 and npi+ final state, but 5.98 for ppi+
esample/Enu_lower_cut [ Global module-variables ]
[ Top ] [ esample ] [ Global module-variables ]
SOURCE
real, save, public :: Enu_lower_cut=0.
PURPOSE
cut events with neutrino energy below Enu_lower_cut; for ANL experiment, for example, Enu_lower_cut=0.5 for some analysis of ppi+
esample/energylimit_for_Qsrec [ Global module-variables ]
[ Top ] [ esample ] [ Global module-variables ]
SOURCE
logical, save, public :: energylimit_for_Qsrec = .False.
PURPOSE
switch for using the energylimits Enu_upper_cut and Enu_lower_cut for the Q^2 reconstruction; values: .true. or .false. (default: .false.)
esample/read_fluxfile [ Functions ]
[ Top ] [ esample ] [ Functions ]
NAME
subroutine read_fluxfile(NDIM,fluxfilename,n_E,enu,flux,sumflux)
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).
INPUTS
- fluxfilename: name of input flux file, with energy and flux pairwise in different lines.
- NDIM: maximal dimension of flux file
OUTPUT
- 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(NDIM,n_E,sumflux,enu)
PURPOSE
returns one energy value by sampling the flux distribution using discrete cumulative sampling
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