gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS

/neutrinoAnalysis [ Modules ]

[ Top ] [ Modules ]

NAME

module neutrinoAnalysis

PURPOSE

  • This module does the analysis of the output of neutrino induced processes.

neutrinoAnalysis/detailed_diff_output [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: detailed_diff_output = .false.

PURPOSE

If .true. then also the detailed output of differential cross sections is produced

neutrinoAnalysis/forPion [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: forPion =.true.

PURPOSE

If .true. then also the detailed output of differential cross sections is produced

neutrinoAnalysis/forEta [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: forEta =.false.

PURPOSE

If .true. then also the detailed output of differential cross sections is produced

neutrinoAnalysis/forKaon [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: forKaon =.false.

PURPOSE

If .true. then also the detailed output of differential cross sections is produced

neutrinoAnalysis/forKaonBar [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: forKaonBar =.false.

PURPOSE

If .true. then also the detailed output of differential cross sections is produced

neutrinoAnalysis/forDmeson [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: forDmeson =.false.

PURPOSE

If .true. then also the detailed output of differential cross sections is produced

neutrinoAnalysis/forDbar [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: forDbar =.false.

PURPOSE

If .true. then also the detailed output of differential cross sections is produced

neutrinoAnalysis/forDs_plus [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: forDs_plus =.false.

PURPOSE

If .true. then also the detailed output of differential cross sections is produced

neutrinoAnalysis/forDs_minus [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: forDs_minus =.false.

PURPOSE

If .true. then also the detailed output of differential cross sections is produced

neutrinoAnalysis/forNucleon [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: forNucleon =.true.

PURPOSE

If .true. then also the detailed output of differential cross sections is produced

neutrinoAnalysis/forLambda [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: forLambda =.false.

PURPOSE

If .true. then also the detailed output of differential cross sections is produced

neutrinoAnalysis/forSigmaResonance [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: forSigmaResonance =.false.

PURPOSE

If .true. then also the detailed output of differential cross sections is produced

neutrinoAnalysis/forXi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: forXi =.false.

PURPOSE

If .true. then also the detailed output of differential cross sections is produced

neutrinoAnalysis/forOmegaResonance [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: forOmegaResonance =.false.

PURPOSE

If .true. then also the detailed output of differential cross sections is produced

neutrinoAnalysis/include_W_dist [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: include_W_dist=.false.

PURPOSE

If .true. then the invariant mass distributions for events with 1 pion and 1 nucleon in the final state are produced

neutrinoAnalysis/dW_Npi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save  :: dW_Npi=0.02

PURPOSE

for dsigma/d(InvariantMass); only work if include_W_dist is .true. set the min, max and steps for various W-distributions

neutrinoAnalysis/Wmin_Npi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save  :: Wmin_Npi=1.08

PURPOSE

for dsigma/d(InvariantMass); only work if include_W_dist is .true. set the min, max and steps for various W-distributions

neutrinoAnalysis/Wmax_Npi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save  :: Wmax_Npi=1.6

PURPOSE

for dsigma/d(InvariantMass); only work if include_W_dist is .true. set the min, max and steps for various W-distributions

neutrinoAnalysis/dW_mupi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save  :: dW_mupi=0.04

PURPOSE

only work if include_W_dist is .true. set the min, max and steps for various W-distributions

neutrinoAnalysis/Wmin_mupi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save  :: Wmin_mupi=0.24

PURPOSE

only work if include_W_dist is .true. set the min, max and steps for various W-distributions

neutrinoAnalysis/Wmax_mupi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save  :: Wmax_mupi=1.2

PURPOSE

only work if include_W_dist is .true. set the min, max and steps for various W-distributions

neutrinoAnalysis/dW_muN [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save  :: dW_muN=0.04

PURPOSE

only work if include_W_dist is .true. set the min, max and steps for various W-distributions

neutrinoAnalysis/Wmin_muN [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save  :: Wmin_muN=1.04

PURPOSE

only work if include_W_dist is .true. set the min, max and steps for various W-distributions

neutrinoAnalysis/Wmax_muN [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save  :: Wmax_muN=2.12

PURPOSE

only work if include_W_dist is .true. set the min, max and steps for various W-distributions

neutrinoAnalysis/numin [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save  :: numin=0.

PURPOSE

for calorimetric analysis: values for transferred energy; only work if calorimetric_analysis is .true. set the min, max and bins for nu distributions

neutrinoAnalysis/numax [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save  :: numax=10.0

PURPOSE

for calorimetric analysis: values for transferred energy; only work if calorimetric_analysis is .true. set the min, max and bins for nu distributions

neutrinoAnalysis/nubin [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save  :: nubin=0.1

PURPOSE

for calorimetric analysis: values for transferred energy; only work if calorimetric_analysis is .true. set the min, max and bins for nu distributions

neutrinoAnalysis/Enumin [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save  :: Enumin=0.

PURPOSE

for calorimetric analysis: values for neutrino energy; only work if calorimetric_analysis is .true. set the min, max and bins for nu distributions

neutrinoAnalysis/Enumax [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save  :: Enumax=10.0

PURPOSE

for calorimetric analysis: values for neutrino energy; only work if calorimetric_analysis is .true. set the min, max and bins for nu distributions

neutrinoAnalysis/Enubin [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save  :: Enubin=0.1

PURPOSE

for calorimetric analysis: values for neutrino energy; only work if calorimetric_analysis is .true. set the min, max and bins for nu distributions

neutrinoAnalysis/EkinMin [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: EkinMin=0.

PURPOSE

if detailed_diff_output is TRUE: minimal kinetic energy for dsigma/dEkin for hadrons

neutrinoAnalysis/EkinMax [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: EkinMax=2.

PURPOSE

if detailed_diff_output is TRUE: Maximal kinetic energy for dsigma/dEkin for hadrons

neutrinoAnalysis/dEkin [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: dEkin=0.01

PURPOSE

if detailed_diff_output is TRUE: Delta(eKin) for dsigma/dEKin for hadrons

neutrinoAnalysis/kineticEnergyDetectionThreshold_nucleon [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: kineticEnergyDetectionThreshold_nucleon=0.0

PURPOSE

kineticEnergyDetectionThreshold lower detection threshold for nucleon kinetic energies

neutrinoAnalysis/kineticEnergyDetectionThreshold_chargedpion [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: kineticEnergyDetectionThreshold_chargedpion=0.0

PURPOSE

kineticEnergyDetectionThreshold

neutrinoAnalysis/kineticEnergyDetectionThreshold_neutralpion [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: kineticEnergyDetectionThreshold_neutralpion=0.0

PURPOSE

kineticEnergyDetectionThreshold

neutrinoAnalysis/kineticEnergyDetectionThreshold_lepton [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: kineticEnergyDetectionThreshold_lepton=0.0

PURPOSE

kineticEnergyDetectionThreshold only lepton kinetic energies above this threshold can be detected This cut affects *all* events, not just the outgoing lepton!

neutrinoAnalysis/AngleUpperDetectionThresholdDegrees_nucleon [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: AngleUpperDetectionThresholdDegrees_nucleon=180.0

PURPOSE

nucleon angles up to this value can be detected

neutrinoAnalysis/AngleUpperDetectionThresholdDegrees_chargedpion [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: AngleUpperDetectionThresholdDegrees_chargedpion=180.0

PURPOSE

charged pion angles up to this value can be detected

neutrinoAnalysis/AngleUpperDetectionThresholdDegrees_neutralpion [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: AngleUpperDetectionThresholdDegrees_neutralpion=180.0

PURPOSE

neutral pion angles angles up to this value can be detected

neutrinoAnalysis/AngleUpperDetectionThresholdDegrees_lepton [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: AngleUpperDetectionThresholdDegrees_lepton=180.0

PURPOSE

lepton angles up to this value can be detected This cut affects *all* events, not just the outgoing lepton!

neutrinoAnalysis/radialScale [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save :: radialScale=1.5

PURPOSE

If radial position of nucleon < radialScale*target radius, then the nucleon is assumed to be bound

neutrinoAnalysis/applyCuts [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  integer, save :: applyCuts = 0

PURPOSE

This parameter encodes 'binary', which cuts should be applied:

Instead of having three indipendent flags (with values=0 or 1) as e.g. labelled 'doLepton', 'doIsBound', 'doBelowThr', applyCuts combines them formally into one number as 

   applyCuts = 1*doLepton + 2*doIsBound + 4*doBelowThr

So by setting any number between 0 and 7, one can individually switch on and off each of these cuts.

'lepton_acceptance' uses the input parameters:

'isBound' tests, whether kinetic energy plus potential is <0

'isBelowThreshold' uses the input parameters:

Some examples:

  • To generate full inclusive output, set the value applyCuts=0
  • To generate output where bound nucleons are dropped, set applyCuts=2
  • To generate output with specific experimental cuts for the outgoing hadrons, set applyCuts=4 or applyCuts=6 and set the corresponding threshold parameters accordingly.
  • If in the experiment also cuts on the outgoing lepton are used, set applyCuts=7 and set the corresponding threshold parameters accordingly.

NOTES

These cuts affects the output into the file "FinalEvents.dat".

The cut 'lepton_acceptance' (de-)selects the full event, while the other two cuts only decide whether a specific particle is accepted or not.

The kinetic energy of a bound nucleon is < 0. Therefore using the default value kineticEnergyDetectionThreshold_nucleon=0.0 also tests, whether the particle is bound or not. Set the parameter to a large negative value to become ineffective.

neutrinoAnalysis/Fissum_analysis [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  Fissum_analysis=.false.

PURPOSE

do analysis with cuts as needed for Fig 25 in Fissum et al, PRC 70, 034606 (2004)

neutrinoAnalysis/ZeroPion_analysis [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  ZeroPion_analysis=.false.

PURPOSE

produce output of xsec for various final states with 0 pions and 2 pions see file see sigma_0pions.dat for the list of the final states

see files neutrino_0pions.dat, neutrino_0pions_QE.dat, neutrino_0pions_Delta.dat, ... for output

neutrinoAnalysis/calorimetric_analysis [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  calorimetric_analysis=.false.

PURPOSE

do calorimetric energy-transfer and neutrino-energy reconstruction (for each QE, Delta, ...) as in the MINOS experiment

neutrinoAnalysis/reconstruct_neutrino_energy [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  reconstruct_neutrino_energy=.false.

PURPOSE

reconstruct neutrino energy for final state in "specificEvent_analysis"

NOTES

.true. must be combined with specificEvent_analysis=.true. and at least one specific event .true.

neutrinoAnalysis/specificEvent_analysis [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  specificEvent_analysis=.false.

PURPOSE

do analysis for specific final states

neutrinoAnalysis/no_pi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  no_pi=.false.

PURPOSE

do analysis for specific final states: specificEvent=1, no_pi (for example, for QE-like MiniBooNE)

neutrinoAnalysis/p_Xn_no_pi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  p_Xn_no_pi=.false.

PURPOSE

do analysis for specific final states: specificEvent=2

neutrinoAnalysis/piplus [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  piplus=.false.

PURPOSE

do analysis for specific final states: specificEvent=3, 1 pi+ X nucleons mesons of other flavor

neutrinoAnalysis/piplus_MULTI [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  piplus_MULTI=.false.

PURPOSE

do analysis for specific final states: specificEvent=4 >=1 pi+ X other pions (incl pi+) X nucleons

neutrinoAnalysis/pi0 [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  pi0=.false.

PURPOSE

do analysis for specific final states: specificEvent=5, 1 pi0 X nucleons, plus mesons of other flavor

neutrinoAnalysis/pi0_MULTI [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  pi0_MULTI=.false.

PURPOSE

do analysis for specific final states: specificEvent=6, >=1 pi0 X other pions X nucleons, (pi0 K2K)

neutrinoAnalysis/piminus [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  piminus=.false.

PURPOSE

do analysis for specific final states: specificEvent=7 1 pi- X other pions X nucleons

neutrinoAnalysis/piminus_MULTI [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  piminus_MULTI=.false.

PURPOSE

do analysis for specific final states: specificEvent=8 >=1 pi- X other pions X nucleons

neutrinoAnalysis/pp_no_pi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  pp_no_pi=.false.

PURPOSE

do analysis for specific final states: specificEvent=9 2 protons, X neutrons, 0 pions

neutrinoAnalysis/pn_no_pi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  pn_no_pi=.false.

PURPOSE

do analysis for specific final states: specificEvent=10 1 neutron, 1 proton, 0 pions

neutrinoAnalysis/nn_no_pi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  nn_no_pi=.false.

PURPOSE

do analysis for specific final states: specificEvent=11 2 neutrons, X protons, 0 pions

neutrinoAnalysis/pp_Xn_no_pi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  pp_Xn_no_pi=.false.

PURPOSE

do analysis for specific final states: specificEvent=12 2 protons, X neutrons, 0 pions

neutrinoAnalysis/nn_Xp_no_pi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  nn_Xp_no_pi=.false.

PURPOSE

do analysis for specific final states: specificEvent=13 2 neutrons, X protons, 0 pions

neutrinoAnalysis/ppp_Xn_no_pi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  ppp_Xn_no_pi=.false.

PURPOSE

do analysis for specific final states: specificEvent=14 3 protons, X neutrons, 0 pions

neutrinoAnalysis/pppp_Xn_no_pi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  pppp_Xn_no_pi=.false.

PURPOSE

do analysis for specific final states: specificEvent=15 4 protons, X neutrons, 0 pions

neutrinoAnalysis/p_no_pi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  p_no_pi=.false.

PURPOSE

do analysis for specific final states: specificEvent=16 1 proton, 0 neutron, 0 pion

neutrinoAnalysis/n_no_pi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  n_no_pi=.false.

PURPOSE

do analysis for specific final states: specificEvent=17 1 neutron, 0 proton, 0 pion

neutrinoAnalysis/Xn_no_pi [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save ::  Xn_no_pi=.false.

PURPOSE

do analysis for specific final states: specificEvent=18, 0 proton, X neutrons, 0 pions

neutrinoAnalysis/excl_hadron [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: excl_hadron=.false.

PURPOSE

do analysis for specific final states: specificEvent=19,20,21 exclusive 1 pion, no other pions or other mesons of different flavor There could be still other mesons which are heavier than the D, Such events (very rare at DUNE energies) could be counted as exclusive single-meson cross section. This could be cured by extending the list of stable mesons

neutrinoAnalysis/full_incl [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: full_incl=.true.

PURPOSE

do analysis for specific final states: specificEvent=22 fully inclusive event, all hadrons in final state

neutrinoAnalysis/QEp [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical ::  QEp=.false.

PURPOSE

if .true, do analysis for specific analysis for QE-like event with 1 mu, 0 pi, X p

neutrinoAnalysis/pcut [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save ::  pcut = 0.0

PURPOSE

determines lower acceptance cut for outgoing protons

values can be changed in the namelist neutrinoAnalysis via the variable kineticEnergyDetectionThreshold_nucleon

neutrinoAnalysis/binsizeQ2 [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save    ::  binsizeQ2=0.01

PURPOSE

do analysis for specific final states: binning for reconstruction of Q2 and Enu

neutrinoAnalysis/binsizeEnu [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save    ::  binsizeEnu=0.02

PURPOSE

do analysis for specific final states: binning for reconstruction of Q2 and Enu

neutrinoAnalysis/maxQ2 [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save    ::  maxQ2=5.0

PURPOSE

do analysis for specific final states: max values for reconstruction of Q2 and Enu

neutrinoAnalysis/maxEnu [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  real, save    ::  maxEnu=5.0

PURPOSE

do analysis for specific final states: max values for reconstruction of Q2 and Enu

neutrinoAnalysis/outputEvents [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save  :: outputEvents = .false.

PURPOSE

If .true. then all events are written to the file 'FinalEvents.dat'.

neutrinoAnalysis/Xsection_analysis [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save  :: Xsection_analysis = .false.

PURPOSE

If .true. then files "..._total_Xsection_..." and "..._dSigmadEkin_..." are printed.

neutrinoAnalysis/doPipe [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save  :: doPipe = .false.

PURPOSE

If .true. then events are not written to the file 'FinalEvents.dat', but insted written into a named pipe (fifo) with the name fileNamePipe.

neutrinoAnalysis/fileNamePipe [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  character(len=1000), save :: fileNamePipe = ""

PURPOSE

name of the pipe to be used

neutrinoAnalysis/DoOutChannels [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

  logical, save :: DoOutChannels = .false.

PURPOSE

switch on/off: reporting of all final state channels

neutrinoAnalysis/readinput [ Subroutines ]

[ Top ] [ neutrinoAnalysis ] [ Subroutines ]

NAME

subroutine readinput

INPUTS

NONE

OUTPUT

NONE

PURPOSE

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

neutrinoAnalysis/NeutrinoAnalysis [ Namelists ]

[ Top ] [ neutrinoAnalysis ] [ Namelists ]

NAME

NAMELIST /NeutrinoAnalysis/

PURPOSE

This namelist includes:

neutrinoAnalysis/W_distributions [ Namelists ]

[ Top ] [ neutrinoAnalysis ] [ Namelists ]

NAME

NAMELIST /W_distributions/

PURPOSE

This Namelist includes:

neutrinoAnalysis/nl_calorimetric_analysis [ Namelists ]

[ Top ] [ neutrinoAnalysis ] [ Namelists ]

NAME

NAMELIST /nl_calorimetric_analysis/

PURPOSE

This Namelist includes:

neutrinoAnalysis/nl_specificEvent [ Namelists ]

[ Top ] [ neutrinoAnalysis ] [ Namelists ]

NAME

NAMELIST /nl_specificEvent/

PURPOSE

This namelist includes:

neutrinoAnalysis/detailed_diff [ Namelists ]

[ Top ] [ neutrinoAnalysis ] [ Namelists ]

NAME

NAMELIST /detailed_diff/

PURPOSE

This namelist includes:

neutrinoAnalysis/neutrino_Analyze [ Subroutines ]

[ Top ] [ neutrinoAnalysis ] [ Subroutines ]

NAME

subroutine neutrino_Analyze(Particles,finalFlag,num_runs_sameEnergy)

INPUTS

  • type(particle), intent(in),dimension(:,:) :: Particles -- Particles which shall be analyzed
  • logical, intent(in) :: finalFlag -- if .true. than the final output for a series of calls will be done
  • integer :: num_runs_sameEnergy

NOTES

  • This subroutine produces output for neutrino-nucleus scattering.

neutrinoAnalysis/Neutrino_total_Xsection_multiplicities.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file Neutrino_total_Xsection_multiplicities.dat

PURPOSE

The file is produced in the runs with eventtype=5=neutrino .

The file shows the cross sections for preselected final states

Units:

  • For process_ID=CC and NC the units 10^{-38} cm^2 for integrated xsec (10^{-38)cm^2/GeV for dsigma/dElepton, 10^{-38)cm^2/GeV^2 for dsigma/dQ^2, and so on)
  • For process_ID=EM the units are nanobarns=10^{-33}cm^2

Columns:

Some columns in detail:

  • #2: 1 pi-, no other pions of any charge and anything else
  • #3: 1 pi0, no other pions of any charge and anything else
  • #4: 1 pi+, no other pions of any charge and anything else
  • #5: 1 eta and anything else
  • #14: 1 neutron, no other nucleons and anything else
  • #15: 1 proton, no other nucleons and anything else
  • #68: 1 nucleon and 1 pion, no other pions or nucleons, anything else
  • #68: 1 proton (no other nucleons) and 0 pions and anything else (QE-like in Cherenkov detector)
  • #69: 0 pions and anything else
  • #70: at least 1 pi-, any number of pi0 and/or pi+ , anything else, each event is counted once
  • #71: at least 1 pi0, any number of pi- and/or pi+ , anything else, each event is counted once
  • #72: at least 1 pi+, any number of pi- and/or pi0 , anything else, each event is counted once
  • #91: at least 1 pi0, any number of pi- and/or pi+ , anything else, each pi0 is counted
  • #96: no nucleons, anything else
  • #97: 5 or more nucleons, anything else
  • #120-#239: errors to Columns #2-#119

neutrinoAnalysis/Neutrino_total_Xsection_multiplicities_QE.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file Neutrino_total_Xsection_multiplicities_QE.dat

PURPOSE

The same as Neutrino_total_Xsection_multiplicities.dat but for QE events (=the first interaction was quasielastic or elastic scattering)

neutrinoAnalysis/neutrino_total_Xsection_multiplicities_Delta.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file Neutrino_total_Xsection_multiplicities_Delta.dat

PURPOSE

The same as Neutrino_total_Xsection_multiplicities.dat but for Delta events (=the first interaction was production of the Delta resonance)

neutrinoAnalysis/Neutrino_total_Xsection_multiplicities_highRES.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file Neutrino_total_Xsection_multiplicities_highRES.dat

PURPOSE

The same as Neutrino_total_Xsection_multiplicities.dat but for highRES events (=the first interaction was production any resonance beyond Delta)

neutrinoAnalysis/Neutrino_total_Xsection_multiplicities_1piBG.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file Neutrino_total_Xsection_multiplicities_1piBG.dat

PURPOSE

The same as Neutrino_total_Xsection_multiplicities.dat but for 1piBG events (=the first interaction was background production of 1-pion final state)

neutrinoAnalysis/Neutrino_total_Xsection_multiplicities_DIS.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file Neutrino_total_Xsection_multiplicities_DIS.dat

PURPOSE

The same as Neutrino_total_Xsection_multiplicities.dat but for DIS events (=the first interaction was DIS)

neutrinoAnalysis/Neutrino_total_Xsection_multiplicities_2p2h.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file Neutrino_total_Xsection_multiplicities_2p2h.dat

PURPOSE

The same as Neutrino_total_Xsection_multiplicities.dat but for 2particle-2hole events

neutrinoAnalysis/Neutrino_total_Xsection_multiplicities_gen0.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file Neutrino_total_Xsection_multiplicities_gen0.dat

PURPOSE

The same as Neutrino_total_Xsection_multiplicities.dat for particles of the 0th generation The definition of generation is given in history.f90, description of the module "history"

neutrinoAnalysis/Neutrino_total_Xsection_multiplicities_gen1.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file Neutrino_total_Xsection_multiplicities_gen1.dat

PURPOSE

The same as Neutrino_total_Xsection_multiplicities.dat for particles of the 1st generation The definition of generation is given in history.f90, description of the module "history"

neutrinoAnalysis/Neutrino_total_Xsection_multiplicities_gen2.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file Neutrino_total_Xsection_multiplicities_gen2.dat

PURPOSE

The same as Neutrino_total_Xsection_multiplicities.dat for particles of the 2nd generation The definition of generation is given in history.f90, description of the module "history"

neutrinoAnalysis/Neutrino_total_Xsection_multiplicities_gen3ormore.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file Neutrino_total_Xsection_multiplicities_gen3ormore.dat

PURPOSE

The same as Neutrino_total_Xsection_multiplicities.dat for particles of the 3rd-or-more generation The definition of generation is given in history.f90, description of the module "history"

neutrinoAnalysis/Neutrino_total_Xsection.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file Neutrino_total_Xsection.dat

PURPOSE

The file is produced in the runs with eventtype=5=neutrino .

The file shows the total cross sections after final state interactions for pions, protons and neutrons The cross sections are semi-inclusive, i.e. for pi+, they may contain events with only one pi+, but also events with pi+pi0.

Units:

  • For process_ID=CC and NC: 10^{-38} cm^2/GeV for integrated xsec and so on ..
  • For process_ID=EM: nanobarns=10^{-33}cm^2/ GeV and so on.
  • All cross sections are given per nucleon (1/A)
  • The various columns in this file are labeld by the particle species

neutrinoAnalysis/Neutrino_kinetic_energy_Xsection.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file Neutrino_kinetic_energy_Xsection.dat

PURPOSE

The file is produced in the runs with eventtype=5=neutrino .

The file shows the kinetic energy differential cross sections after final state interactions for pions, protons and neutrons

Units:

  • For process_ID=CC and NC: 10^{-38} cm^2/GeV for integrated xsec and so on ..
  • For process_ID=EM: nanobarns=10^{-33}cm^2/ GeV and so on.
  • All cross sections are given per nucleon (1/A)

Columns:

  • #1: variable which was raised (e.g. Q^2 for nuXsectionMode=3=dSigmadQs mode, Elepton for nuXsectionMode=2=dSigmadQsdElepton and so on)
  • #2: kinetic energy in GeV
  • #3: xsec for events with at least one pi- in the final state (coincides with column 2 in diff_XXX_dSigma_dEkin_pi_charge_-1_MULTI.dat)
  • #4: xsec for events with at least one pi0 (coincides with column 2 in diff_XXX_dSigma_dEkin_pi_charge_+0_MULTI.dat)
  • #5: xsec for events with at least one pi+ (coincides with column 2 in diff_XXX_dSigma_dEkin_pi_charge_+1_MULTI.dat)
  • #6: xsec for events with at least one neutron (coincides with column 2 in diff_XXX_dSigma_dEkin_N_charge_+0_MULTI.dat)
  • #7: xsec for events with at least one proton (coincides with column 2 in diff_XXX_dSigma_dEkin_N_charge_+1_MULTI.dat)

HERE XXX=000,001, ... is the count of the raise variable

neutrinoAnalysis/Neutrino_Xsection_numbers.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file Neutrino_Xsection_numbers.dat

PURPOSE

The file is produced in the runs with eventtype=5=neutrino .

The file shows number of nucleons produced in one run before and after final state interactions (should be devided by numEnsembles to obtain the average multiplicity per target nucleus and divided further by target_A to obtain the average multiplicity per target nucleon)

Columns:

  • #1: variable which was raised (e.g. Q^2 for nuXsectionMode=3=dSigmadQs mode, Elepton for nuXsectionMode=2=dSigmadQsdElepton and so on)
  • #2: all the protons produced before final state interactions (FSI)
  • #3: protons that made it out of the nucleus after FSI
  • #4: protons produced in generation 0 and made it out of the nucleus after FSI
  • #5: protons produced in generation 1 and made it out of the nucleus after FSI
  • #6: protons produced in generation 2 and made it out of the nucleus after FSI
  • #7: protons produced in generation 3-or-more and made it out of the nucleus after FSI
  • #8: all the neutrons produced before FSI
  • #9: neutrons that made it out of the nucleus after FSI
  • #10: neutrons produced in generation 0 and made it out of the nucleus after FSI
  • #11: neutrons produced in generation 1 and made it out of the nucleus after FSI
  • #12: neutrons produced in generation 2 and made it out of the nucleus after FSI
  • #13: neutrons produced in generation 3-or-more and made it out of the nucleus after FSI
  • #14: all the antiprotons produced before FSI
  • #15: antiprotons that made it out of the nucleus after FSI
  • #16: antiprotons produced in generation 0 and made it out of the nucleus after FSI
  • #17: antiprotons produced in generation 1 and made it out of the nucleus after FSI
  • #18: antiprotons produced in generation 2 and made it out of the nucleus after FSI
  • #19: antiprotons produced in generation 3-or-more and made it out of the nucleus after FSI

neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file diff_ZZZ_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat

with:

  • ZZZ denotes the origin (the first interaction vertex) of the event (see description in neutrino.EprimeCostplaneSX.ZZZ.dat). 

   If ZZZ is missing the total cross section
   (sum over all primary events) is given, otherwise ZZZ=Delta, DIS, ....

* XXX=000, 001, 002 --- the first, second, third and so on values of 

   the "raised variable"
   (e.g. Q^2 for nuXsectionMode=3=dSigmadQs mode,
   Elepton for nuXsectionMode=2=dSigmadQsdElepton  and so on)

* HADRON = pion, N, K, K~, Lambda, SigmaResonance, eta 

   only those are shown, for which the switches in the namelist
   "detailed_diff" are set to .true.
   possible hadrons are contained in the field particleIDs defined in
   AnaEvent.f90

* CHARGE = charge of the outgoing hadron

PURPOSE

The file is produced in runs with:

The file shows the cross sections for ___1-HADRON___ final state versus kinetic energy of the outgoing hadron (one HADRON of a given CHARGE and no HADRONs with same flavor, but different charges; 

  there could be additional hadrons with different flavor)

Units:

Columns:

  • #1: kinetic energy of the outgoing HADRON of given CHARGE [GeV]
  • #2: dsi/dEkin xsec
  • #3: number of events contributed (only for internal use, you can safely neglect it)
  • #4: xsec-statistical-error

neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dEkin_HADRON_charge_CHARGE_1X.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE_1X.dat

PURPOSE

The same as diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat but for ___1-HADRON-X___ final state (one HADRON of a given CHARGE and any number of HADRONs of different charges)

neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dEkin_HADRON_charge_CHARGE_2X.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE_2X.dat

PURPOSE

The same as diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat but for ___2-HADRON-X___ final state (two HADRONs of a given CHARGE and any number of HADRONs of different charges)

neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dEkin_HADRON_charge_CHARGE_MULTI.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE_MULTI.dat

PURPOSE

The same as diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat but for ___MULTI-HADRON___ final state (at least one HADRON of a given CHARGE and any number of HADRONs of different charges)

neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dEcostheta_HADRON_charge_CHARGE.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file diff_XXX_dSigma_dEcostheta_HADRON_charge_CHARGE.dat

file diff_XXX_dSigma_dEcostheta_HADRON_charge_CHARGE_1X.dat

file diff_XXX_dSigma_dEcostheta_HADRON_charge_CHARGE_2X.dat

file diff_XXX_dSigma_dEcostheta_HADRON_charge_CHARGE_MULTI.dat)

PURPOSE

The same as diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat but dsigma/d(E(1-cosTheta))

The file is produced in runs with:

The file shows the cross sections dsigma/d(E(1-cosTheta)) versus E*(1-costheta), where E is the energy (full energy, not kinetic) of the outgoing hadron, costheta its polar scattering angle (recall here that in neutrino runs neutrinos are moving along z-direction)

Columns:

  • #1: E*(1-cosTheta) [GeV]
  • #2: dsi/d(E(1-costheta)) xsec
  • #3: number of events contributed (only for internal use, you can safely neglect it)
  • #4: xsec-statistical-error

neutrinoAnalysis/diff_ZZZ_XXX_dSigma_dTheta_HADRON_charge_CHARGE.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file diff_XXX_dSigma_dTheta_HADRON_charge_CHARGE.dat

file diff_XXX_dSigma_dTheta_HADRON_charge_CHARGE_1X.dat

file diff_XXX_dSigma_dTheta_HADRON_charge_CHARGE_2X.dat

file diff_XXX_dSigma_dTheta_HADRON_charge_CHARGE_MULTI.dat

PURPOSE

The same as diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat but dsigma/dTheta

The file is produced in runs with:

The file shows the cross sections dsigma/dTheta versus Theta, where Theta is the polar scattering angle (in radians) of the outgoing hadron (recall here that in neutrino runs neutrinos are moving along z-direction)

Columns:

  • #1: Theta [radians]
  • #2: dsi/dTheta xsec
  • #3: number of events contributed (only for internal use, you can safely neglect it)
  • #4: xsec-statistical-error

neutrinoAnalysis/diff_XXX_dSigma_dEkin_lepton_PPP.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file diff_XXX_dSigma_dEkin_lepton_PPP.dat

with:

PURPOSE

The file is produced in runs with:

The file shows the cross sections for a specific final state versus kinetic energy of the outgoing lepton

Units:

Columns:

  • #1: kinetic energy of the outgoing lepton [GeV]
  • #2: dsi/dEkin xsec per nucleon
  • #3: number of events contributed (only for internal use, you can safely neglect it)
  • #4: xsec-statistical-error

neutrinoAnalysis/diff_XXX_dSigma_dQ2_lepton_PPP.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file diff_XXX_dSigma_dQ2_lepton_PPP.dat

PURPOSE

The same as diff_XXX_dSigma_dEkin_lepton_PPP.dat but for dsigma/dQ2

The file shows the cross sections for a specific final state versus Q2

Units:

Columns:

  • #1: Q2 [GeV^2] squared transfer momentum (Q2 = -q_mu \cdot q^\mu)
  • #2: dsi/dQ2 xsec
  • #3: number of events contributed (only for internal use, you can safely neglect it)
  • #4: xsec-statistical-error

neutrinoAnalysis/diff_XXX_dSigma_dQ2p_lepton_PPP.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file diff_XXX_dSigma_dQ2p_lepton_PPP.dat

PURPOSE

The same as diff_XXX_dSigma_dEkin_lepton_PPP.dat but for dsigma/dQ2

The file shows the cross sections for a final state with 1 mu, 0 pi , and (at least) 1 p versus Q2 here Q2 is calculated from the kinematics of the outgoing leading proton: Q2prot = (mProt - epsB)**2 - mProt**2 & 

            & + 2*(mNeut - epsB)*(Tp + mProt - mNeut + epsB)

This distribution can be converted into a kinetic energy distribution of protons in 0pion events

Units:

Columns:

  • #1: Q2 [GeV^2] squared transfer momentum, from proton kinematics
  • #2: dsi/dQ2 xsec
  • #3: number of events contributed (only for internal use, you can safely neglect it)
  • #4: xsec-statistical-error

neutrinoAnalysis/diff_XXX_dSigma_dcos_lepton_PPP.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file diff_XXX_dSigma_dcos_lepton_PPP.dat

PURPOSE

The same as diff_XXX_dSigma_dEkin_lepton_PPP.dat but for dsigma/dcos(theta_l)

The file shows the cross sections for a specific final state versus cos of the scattering angle of the outgoin lepton (with respect to neutrino direction)

Units:

Columns:

  • #1: cos(theta_l) cos of the scattering angle of the outgoing lepton
  • #2: dsi/dcos(theta_l) xsec
  • #3: number of events contributed (only for internal use, you can safely neglect it)
  • #4: xsec-statistical-error

neutrinoAnalysis/diff_XXX_d2Sigma_dEdcost_lepton_no_pi.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file diff_XXX_d2Sigma_dEdcost_lepton_no_pi.dat

PURPOSE

Double differential cross section d2sigma/(dE dcostheta) for outgoing lepton for 0-pi events

The file contains the cross section for the outgoing lepton versus cos of the scattering angle of the outgoing lepton (with respect to neutrino direction) and its total energy for 0-pion events

Units:

Columns:

  • #1: total energy of the outgoing lepton (in GeV)
  • #2: cos(theta_l) cos of the scattering angle of the outgoing lepton
  • #3: d2si/(dcos(theta_l) dE_l) dd-xsec

the following columns not implemented yet * #4: number of events contributed (only for internal use, you can 

   safely neglect it)

* #5: xsec-statistical-error

neutrinoAnalysis/diff_XXX_d2Sigma_dpLdpT_lepton_no_pi.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file diff_XXX_d2Sigma_dpLdpT_lepton_no_pi.dat

PURPOSE

Double differential cross section d2sigma/(dpL dpT) for outgoing lepton for 0-pi events

The file contains the cross section for the outgoing lepton vs longitudinal and transverse momentum of the outgoing lepton for 0-pion events

Units:

Columns:

  • #1: p_T longitudinal momentum of the outgoing lepton (in GeV)
  • #2: p_L longitudinal momentum of the outgoing lepton (in GeV)
  • #3: d2si//dpL dpT) dd-xsec

the following columns not implemented yet * #4: number of events contributed (only for internal use, you can safely neglect it) * #5: xsec-statistical-error

neutrinoAnalysis/diff_XXX_dSigma_dW_nucleon_pion_charge_CHARGE.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file diff_XXX_dSigma_dW_nucleon_pion_charge_CHARGE.dat

PURPOSE

Similar to diff_XXX_dSigma_dEkin_HADRON_charge_CHARGE.dat (see notations XXX and CHARGE there) but gives pion-nucleon invariant mass (W) distribution (true invariant mass with the sum of the 4-momenta of the outgoing particles, W^2=(p_1out+p_2out) as opposed to widely used W2=mN2+2*mN*nu-Q2 from lepton kinematics ) for events with ___1pion___ AND ___1nucleon___ in the final state

neutrinoAnalysis/diff_XXX_dSigma_dW_muon_nucleon_charge_CHARGE.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file diff_XXX_dSigma_dW_muon_nucleon_charge_CHARGE.dat

PURPOSE

The same as diff_XXX_dSigma_dW_nucleon_pion_charge_CHARGE.dat but for the muon-nucleon invariant mass distribution

NOTES

Note that CHARGE is still a pion charge

neutrinoAnalysis/diff_XXX_dSigma_dW_muon_pion_charge_CHARGE.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file diff_XXX_dSigma_dW_muon_pion_charge_CHARGE.dat

PURPOSE

The same as diff_XXX_dSigma_dW_nucleon_pion_charge_CHARGE.dat but for the muon-pion invariant mass distribution

NOTES

Note that CHARGE is still a pion charge

neutrinoAnalysis/FinalEvents.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file FinalEvents.dat

PURPOSE

The file contains positions and four-vectors of all outgoing particles from an event.

For a cross-section construction all of the events have to be weighted with the 'perweight' values in col. 5. The final differential cross section dsigma/dE, for example, then is obtained by first binning all the events with respect to the energy. Then the perweights have to be summed within each bit. To obtain the differential cross section the sum of weights in each bin has to be divided by the bin width. In addition, the number of runs at same energy (>1 for better statistics) has to be divided out. The resulting cross section is per nucleon (1/A). Its units are 10^{-38} cm^2 for neutrinos and 10^{-33} cm^2 for electrons.

The output of this file is switched on by the switch 'outputEvents'.

Columns:

NOTES

There is always in each event a particle with weight 0. This is the nucleon on which the initial interaction happened.

In the case of an initial 2p2h process the second initial-state nucleon is not written out. It is chosen to be at the same place as the first initial nucleon (the one with weight=0), with a randomly chosen momentum in the Fermi-sea.

For large target_A, numEnsembles and num_runs_SameEnergy, this file may become very large, e.g.:

  • target_A=12, 1000ens x 20runs, QE and Delta results in 41 Mb
  • target_A=56, 2000ens x 5runs, QE,Delta,highRES,1piBG,DIS results in 500 Mb

neutrinoAnalysis/production_ID [ Global module-variables ]

[ Top ] [ neutrinoAnalysis ] [ Global module-variables ]

SOURCE

PURPOSE

contains info on the very first neutrino-interaction with the nucleus:

  • 1: nucleon (QE)
  • 2-31: non-strange baryon resonance (as in IdTable)
  • 32: pi neutron-background (e.g. nu + n -> mu + pi+ + n)
  • 33: pi proton-background (e.g. nu + n -> mu + pi0 + p)
  • 34: DIS
  • 35: 2p2h QE
  • 36: 2p2h Delta
  • 37: two pion background

neutrinoAnalysis/reconstruction_Enureal_PPP.ZZZ.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file reconstruction_Enureal_PPP.ZZZ.dat

PURPOSE

Similar to reconstruction_Enurec_PPP.ZZZ.dat, but for true neutrino energy

The file shows the event distribution (flux times xsec) versus true neutrino energy for a specific final state The file is produced only for the runs with a neutrino flux (in the namelist "neutrino_induced" nuXsectionMode > 10, nuExp>0)

The file is not affected by cuts on outgoing lepton's angle and energy

Columns:

  • #1: true neutrino energy [GeV]
  • #2: event distribution: normalized flux times xsec
  • #3: number of events contributed (only for internal use, you can safely neglect it)
  • #4: statistical-error of #2

neutrinoAnalysis/reconstruction_Enurec_PPP.ZZZ.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file reconstruction_Enurec_PPP.ZZZ.dat

with:

PURPOSE

The file shows the event distribution (normalized flux times xsec) versus reconstructed neutrino energy (see arXiv:1208.3678 [nucl-th]) for a specific final state

Reconstruction method depends on a specific process; * generally for no_pion events it is based on QE-like kinematics * and for pion events it is based on on-shell-Delta-creation assumption

The file is produced in runs with: 

   The file is not affected by cuts on outgoing lepton's angle and energy

Units:

Columns:

  • #1: reconstructed neutrino energy [GeV]
  • #2: event distribution: flux times xsec)
  • #3: number of events contributed (only for internal use, you can safely neglect it)
  • #4: statistical-error of #2

neutrinoAnalysis/reconstruction_Enu_rec_versus_real_PPP.ZZZ.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file reconstruction_Enu_rec_versus_real_PPP.ZZZ.dat

PURPOSE

Similar to reconstruction_Enurec_PPP.ZZZ.dat, but double-differential

The file shows the 2-D density of the flux times xsec (see arXiv:1208.3678 [nucl-th]) versus true and reconstructed neutrino energies for a specific final state The file is produced only for the runs with a neutrino flux (in the namelist "neutrino_induced" nuXsectionMode > 10, nuExp>0)

Units:

Columns:

  • #1: true neutrino energy [GeV]
  • #2: reconstructed neutrino energy [GeV]
  • #3: flux-folded xsec
  • #4: number of events contributed (only for internal use, you can safely neglect it)
  • #5: xsec-statistical-error

NOTES

up to normalization this is "migration matrix" between true and reconstructed energies The file is not affected by cuts on outgoing lepton's angle and energy

neutrinoAnalysis/reconstruction_Q2real_PPP.ZZZ.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file reconstruction_Q2real_PPP.ZZZ.dat

PURPOSE

Similar to reconstruction_Enurec_PPP.ZZZ.dat, but for the true Q2

The file contains the flux averaged cross section dsigma/dQ2, i.e. integral \int \Phi(E) dsigma/dQ2 (E) dE, versus true Q2 for a specific final state The file is not affected by cuts on outgoing lepton's angle and energy

Units:

Columns:

  • #1: true Q2 [GeV^2] (squared momentum transfer Q2 = - q_mu \cdot q^mu)
  • #2: flux-averaged dsigma/dQ2true
  • #3: number of events contributed (only for internal use, you can safely neglect it)
  • #4: xsec-statistical-error

neutrinoAnalysis/reconstruction_Q2rec_PPP.ZZZ.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file reconstruction_Q2rec_PPP.ZZZ.dat

PURPOSE

Similar to reconstruction_Enurec_PPP.ZZZ.dat, but for the reconstructed Q2 (see arXiv:1208.3678 [nucl-th])

The file contains the flux averaged cross section dsigma/dQ2, i.e. integral \int \Phi(E) dsigma/dQ2 (E) dE, versus reconstructed Q2 for a specific final state The file is not affected by cuts on outgoing lepton's angle and energy

Units:

Columns:

  • #1: reconstructed Q2 [GeV^2] = - q_mu \cdot q^mu)
  • #2: flux-averaged dsigma/dQ2rec
  • #3: number of events contributed (only for internal use, you can safely neglect it)
  • #4: xsec-statistical-error

neutrinoAnalysis/reconstruction_Q2_rec_versus_real_PPP.ZZZ.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file reconstruction_Q2_rec_versus_real_PPP.ZZZ.dat

PURPOSE

Similar to reconstruction_Enurec_PPP.ZZZ.dat, but double-differential in Q2

The file shows the dsigma/dQ2true dQ2rec distribution versus true and reconstructed Q^2 for a specific final state The file is not affected by cuts on outgoing lepton's angle and energy

Units:

Columns:

  • #1: true Q2 [GeV^2]
  • #2: reconstructed Q2 [GeV^2]
  • #3: dsigma/dQ2true dQ2rec
  • #4: number of events contributed (only for internal use, you can safely neglect it)
  • #5: xsec-statistical-error

neutrinoAnalysis/oscillations_UUU_real_PPP.ZZZ.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file oscillations_UUU_real_PPP.ZZZ.dat

with:

  • UUU = "mumu", "mue", "mue_max", "mue_antimax" for oscillation of a given experimental flux for muon_neutrino survival ("mumu"), mu_e appearance with delta_CP=0 ("mue"), delta_CP=pi ("mue_max"), delta_CP=-pi ("mue_antimax")
  • PPP = no_pi, p_Xn_no_pi, piplus, pi0, .... (see namelist nl_specificEvent for the full list) standing for the specific final state under consideration
  • ZZZ=000 - 008 is the origin (the first interaction vertex) of the event (see description in neutrino.EprimeCostplaneXS.ZZZ.dat) 

   The file is not affected by cuts on outgoing lepton's angle and energy

PURPOSE

Similar to reconstruction_Enureal_PPP.ZZZ.dat, but for oscillated neutrino flux

Columns:

  • #1: true neutrino energy [GeV]
  • #2: flux-folded xsec
  • #3: number of events contributed (only for internal use, you can safely neglect it)
  • #4: xsec-statistical-error

neutrinoAnalysis/oscillations_UUU_rec_PPP.ZZZ.dat [ Output files ]

[ Top ] [ neutrinoAnalysis ] [ Output files ]

NAME

file oscillations_UUU_rec_PPP.ZZZ.dat

with:

  • UUU = "mumu", "mue", "mue_max", "mue_antimax" for oscillation of a given experimental flux for muon_neutrino survival ("mumu"), mu_e appearance with delta_CP=0 ("mue"), delta_CP=pi ("mue_max"), delta_CP=-pi ("mue_antimax")
  • PPP = no_pi, p_Xn_no_pi, piplus, pi0, .... (see namelist nl_specificEvent for the full list) standing for the specific final state under consideration
  • ZZZ=000 - 008 is the origin (the first interaction vertex) of the event (see description in neutrino.EprimeCostplaneXS.ZZZ.dat) The file is not affected by cuts on outgoing lepton's angle and energy

PURPOSE

Similar to reconstruction_Enurec_PPP.ZZZ.dat, but for oscillated neutrino flux

Columns:

  • #1: reconstructed neutrino energy [GeV]
  • #2: event rate: flux times xsec
  • #3: number of events contributed (only for internal use, you can safely neglect it)
  • #4: xsec-statistical-error

neutrinoAnalysis/IsBound [ Functions ]

[ Top ] [ neutrinoAnalysis ] [ Functions ]

NAME

logical function IsBound(part)

PURPOSE

return whether particle is bound (.true.) or not (.false.) This cut affects FinalEvents.dat

neutrinoAnalysis/IsBelowThreshold [ Functions ]

[ Top ] [ neutrinoAnalysis ] [ Functions ]

NAME

logical function IsBelowThreshold(part)

PURPOSE

Returns .true. when a particle is below a given detection threshold This routine can be used to remove particles with kinetic energies and angles below specified detection thresholds. This cut affects FinalEvents.dat

neutrinoAnalysis/lepton_acceptance [ Functions ]

[ Top ] [ neutrinoAnalysis ] [ Functions ]

NAME

logical function lepton_acceptance

PURPOSE

Returns .true. when the lepton is above a given kinetic energy and below a given angle, i.e. it is accepted

This routine can be used to remove all events with lepton kinetic energies below and angles below specified detection thresholds. If an event is not accepted it is not added to the full event in subroutine neutrino_Analyze. This cut affects FinalEvents.dat

neutrinoAnalysis/oscillationProbability [ Subroutines ]

[ Top ] [ neutrinoAnalysis ] [ Subroutines ]

NAME

subroutine oscillationProbability(Enu,L,deltaCP, Posc_mumu,Posc_mue,Posc_mue_max,Posc_mue_antimax)

PURPOSE

Prepare Calculation of Oscillation Probability

neutrinoAnalysis/PrintVals10 [ Subroutines ]

[ Top ] [ neutrinoAnalysis ] [ Subroutines ]

NAME

PURPOSE

abbreviation used for many outputfiles: printing the vals stored in sigma, including its error

neutrinoAnalysis/PrintHeader10 [ Subroutines ]

[ Top ] [ neutrinoAnalysis ] [ Subroutines ]

NAME

PURPOSE

abbreviation used for many outputfiles: printing a header line