gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS


/esample [ Modules ]

[ Top ] [ 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=10000.

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+, 6000 for FASER


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/nl_fluxcuts [ Namelists ]

[ Top ] [ esample ] [ Namelists ]

NAME

NAMELIST /nl_fluxcuts/

PURPOSE

This Namelist includes:


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

RESULT

  • eneut: one sampled neutrino energy