gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS

/initNeutrino [ Modules ]

[ Top ] [ Modules ]

NAME

module initNeutrino

PURPOSE

This module is the main module for neutrino-induced reactions. It provides the calculation of various inclusive neutrino cross sections, depending only on the exchanged gauge boson momentum. It also sets the produced outgoing particles into the pertubativeVector.

initNeutrino/process_ID [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  integer, save :: process_ID = 2

PURPOSE

Determine the process (cf. module leptonicID):

  • 1 = EM
  • 2 = CC
  • 3 = NC
  • -1 = antiEM
  • -2 = antiCC
  • -3 = antiNC

initNeutrino/flavor_ID [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  integer, save :: flavor_ID = 2

PURPOSE

Determine the lepton flavor:

  • 1 = electron
  • 2 = muon
  • 3 = tau

initNeutrino/nuXsectionMode [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  integer, save :: nuXsectionMode = 0

PURPOSE

To choose which kind of Xsection is calculated. All values set in module neutrino_IDTable.f90

possible values:

calculation for specific experiments taking into account the flux (choose your favorite experiment with flag nuExp):

  • 10 = EXP_dSigmadEnu
  • 11 = EXP_dSigmadCosThetadElepton
  • 12 = EXP_dSigmadQ2dElepton
  • 14 = EXP_dSigmadCosTheta
  • 15 = EXP_dSigmadElepton
  • 16 = EXP_dSigmaMC
  • 17 = EXP_dSigmaMC_dW
  • 13 = EXP_dSigmaMC_dQ2

initNeutrino/nuExp [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  integer, save :: nuExp = 0

PURPOSE

  • 0 = no specific experiment
  • 1 = MiniBooNE neutrino flux (in neutrino mode = positive polarity)
  • 2 = ANL
  • 3 = K2K
  • 4 = BNL
  • 5 = MiniBooNE antienutrino flux (in antineutrino mode = negative polarity)
  • 6 = MINOS muon-neutrino in neutrino mode
  • 7 = MINOS muon-antineutrino in neutrino mode
  • 8 = NOVA neutrino (medium energy NuMI, 14 mrad off-axis), FD
  • 9 = T2K neutrino off-axix 2.5 degrees ( at ND280 detector )
  • 10 = uniform distribution from Eflux_min to Eflux_max (see namelist nl_neutrino_energyFlux in the module expNeutrinoFluxes)
  • 11 = MINOS muon-neutrino in antineutrino mode
  • 12 = MINOS muon-antineutrino in antineutrino mode
  • 13 = MINERvA muon neutrino, old flux
  • 14 = MINERvA muon antineutrino, old flux
  • 15 = LBNF/DUNE in neutrino mode
  • 16 = LBNF/DUNE in antineutrino mode
  • 17 = LBNO neutrino in neutrino mode
  • 18 = NOMAD
  • 19 = BNB nue BNB= Booster Neutrino Beam
  • 20 = BNB nuebar
  • 21 = BNB numu
  • 22 = BNB numubar
  • 23 = NOvA ND
  • 24 = T2K on axis
  • 25 = MINERvA, 2016 flux
  • 26 = FASERnu
  • 99 = user provided input file

initNeutrino/debugFlag [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  logical, parameter :: debugFlag = .false.

PURPOSE

To switch on debugging information

initNeutrino/includeQE [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  logical, save :: includeQE = .true.

PURPOSE

include QE scattering

initNeutrino/includeDelta [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  logical, save :: includeDelta = .true.

PURPOSE

include Delta excitation

initNeutrino/includeRes [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  logical, save :: includeRes = .true.

PURPOSE

include excitation of higher resonances

initNeutrino/include1pi [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  logical, save :: include1pi = .false.

PURPOSE

include one-pion cross section see neutrinoXsection.f90 for details: there one might choose between different models and also whether it is taken as background or as total cross section

initNeutrino/include2pi [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  logical, save :: include2pi = .false.

PURPOSE

include 2 pion background channel

initNeutrino/includeDIS [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  logical, save :: includeDIS = .false.

PURPOSE

include DIS contribution

initNeutrino/include2p2hQE [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  logical, save :: include2p2hQE = .false.

PURPOSE

include 2p2h QE contribution

initNeutrino/include2p2hDelta [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  logical, save :: include2p2hDelta = .false.

PURPOSE

include 2p2h Delta contribution

initNeutrino/realRun [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  logical, save :: realRun = .false.

PURPOSE

Do not initialize the final state particles as perturbative particles but as real ones.

initNeutrino/printAbsorptionXS [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  logical, save :: printAbsorptionXS = .false.

PURPOSE

flag to produce output about inclusive (absorption) cross sections

initNeutrino/printInclHist [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  logical, save :: printInclHist = .true.

PURPOSE

flag to produce additional output about inclusive cross sections

only checked, if printAbsorptionXS = T

initNeutrino/storeNucleon [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  integer, save :: storeNucleon = 2

PURPOSE

indicate which kind of struck nucleon to save:

  • 1: free Nucleon (i.e. potential removed)
  • 2: bound nucleon

NOTES

real check of energy and momentum conservation only possible with '2'

initNeutrino/max_finalstate_ID [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  integer,parameter :: max_finalstate_ID=37
  ! * Parameter determines the reaction mechanism and kind of final states
  ! * Final states are numbered (often IP) by
  ! * 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

initNeutrino/sigmacut [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  real, save :: sigmacut=10e-4
  ! in usual units (10^(-33) cm^2 for electrons, 10^(-38) cm^2 for neutrinos)

PURPOSE

events with a cross section smaller than this value are skipped.

neutrinoAnalysis/cost_min [ Global module-variables ]

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

SOURCE

  real, save :: cost_min= -1.0

PURPOSE

if detailed_diff_output is TRUE: Minimal cos(theta) of outgoing leptons, used in 2D dsigma/dEdcos(theta) This cut affects *only* the outgoing lepton

neutrinoAnalysis/cost_max [ Global module-variables ]

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

SOURCE

  real, save :: cost_max= +1.0

PURPOSE

if detailed_diff_output is TRUE: Maximal cos(theta) of outgoing leptons, used in 2D dsigma/dEdcos(theta) This cut affects *only* the outgoing lepton

neutrinoAnalysis/delta_cost [ Global module-variables ]

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

SOURCE

  real, save :: delta_cost = 0.1      ! for MB,T2K, for higher energies smaller

PURPOSE

if detailed_diff_output is TRUE: stepsize of cos(theta) of outgoing leptons, used in 2D dsigma/dEdcos(theta)

neutrinoAnalysis/Elept_min [ Global module-variables ]

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

SOURCE

  real, save :: Elept_min = 0.0

PURPOSE

if detailed_diff_output is TRUE: minimal energy of outgoing leptons, used in 2D dsigma/dEdcos(theta)

neutrinoAnalysis/Elept_max [ Global module-variables ]

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

SOURCE

  real, save :: Elept_max = 2.0     ! for MB,T2K, for higher energies larger

PURPOSE

if detailed_diff_output or printAbsorption are TRUE: maximal energy of outgoing leptons, used in 2D dsigma/dEdcos(theta)

neutrinoAnalysis/delta_Elept [ Global module-variables ]

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

SOURCE

  real, save :: delta_Elept = 0.01   ! for MB,T2K

PURPOSE

if detailed_diff_output or printAbsorption are TRUE: stepsize of energy of outgoing leptons, used in 2D dsigma/dEdcos(theta)

neutrinoAnalysis/pL_min [ Global module-variables ]

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

SOURCE

  real, save :: pL_min = 0.0

PURPOSE

if detailed_diff or printAbsorption are TRUE: minimal long. momentum of outgoing leptons, used in 2D dsigma/dpLdpT

neutrinoAnalysis/pL_max [ Global module-variables ]

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

SOURCE

  real, save :: pL_max = 20.0

PURPOSE

if detailed_diff_output or printAbsorption are TRUE: maximal long. momentum of outgoing leptons, used in 2D dsigma/dpLdpT

neutrinoAnalysis/delta_pL [ Global module-variables ]

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

SOURCE

  real, save :: delta_pL = 0.25

PURPOSE

if detailed_diff_output or printAbsorption are TRUE: stepsize of long. momentum of outgoing leptons, used in 2D dsigma/dpLdpT

neutrinoAnalysis/pT_min [ Global module-variables ]

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

SOURCE

  real, save :: pT_min = 0.0

PURPOSE

if detailed_diff or printAbsorption are TRUE: minimal transv. momentum of outgoing leptons, used in 2D dsigma/dpLdpT

neutrinoAnalysis/pT_max [ Global module-variables ]

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

SOURCE

  real, save :: pT_max = 2.5

PURPOSE

if detailed_diff_output or printAbsorption are TRUE: maximal transv. momentum of outgoing leptons, used in 2D dsigma/dpLdpT

neutrinoAnalysis/delta_pT [ Global module-variables ]

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

SOURCE

  real, save :: delta_pT = 0.1

PURPOSE

if detailed_diff_output or printAbsorption are TRUE: binwidth of transv. momentum of outgoing leptons, used in 2D dsigma/dpLdpT

neutrinoAnalysis/Q2_Max [ Global module-variables ]

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

SOURCE

  real, save :: Q2_Max = 100.

PURPOSE

maximal value of Q2 in Q2-distribution

initNeutrino/FileNameFlux [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

PURPOSE

The absolute filename of the file containing flux info, if user supplied

possible values:

  • if given, but does not contain '/': default is '[path_To_Input]/[FileNameFlux]'
  • otherwise: filename is absolute, including path ('~' is okay)

NOTE if you want to use the file 'XXX.dat' in the actual directory, give it as './XXX.dat'

SOURCE

  character(1000), save :: FileNameFlux = ''

initNeutrino/Enumax [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

PURPOSE

maximum of neutrino energy in flux distribution, in GeV

SOURCE

  real , save :: Enumax = 100.

initNeutrino/delta_Enumax [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

PURPOSE

bin width of neutrino energy in flux distribution, in GeV

SOURCE

  real , save :: delta_Enumax = 0.02

initNeutrino/equalWeights_Mode [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  integer, save :: equalWeights_Mode = 0

PURPOSE

possible values are:

  • 0: default perweight mode is used (default)
  • 1: default perweight mode is used, but max is printed
  • 2: MC rejection method is used.

In the default mode, the perweights of the final particles are given by cross section/(A * numEnsembles)

If equalWeightsMode==2, then the perweights are given by equalWeights_Max/(A * numEnsembles)

Please check in the output the line "numberOfSuccess = ..." for the number of events actually generated.

initNeutrino/equalWeights_Max [ Global module-variables ]

[ Top ] [ initNeutrino ] [ Global module-variables ]

SOURCE

  real, save :: equalWeights_Max = -1e99

PURPOSE

The maximum value the MC-rejection method is done against.

initNeutrino/readInput [ Subroutines ]

[ Top ] [ initNeutrino ] [ Subroutines ]

NAME

subroutine readInput

PURPOSE

This subroutine reads input out of jobcard from namelist 'neutrino_induced'.

initNeutrino/neutrino_induced [ Namelists ]

[ Top ] [ initNeutrino ] [ Namelists ]

NAME

NAMELIST /neutrino_induced/

PURPOSE

This Namelist includes:

initNeutrino/lepton_bin [ Namelists ]

[ Top ] [ initNeutrino ] [ Namelists ]

NAME

NAMELIST /lepton_bin/

PURPOSE

This Namelist includes:

initNeutrino/init_neutrino [ Subroutines ]

[ Top ] [ initNeutrino ] [ Subroutines ]

NAME

subroutine init_neutrino(realParticles,pertParticles,raiseFlagIn, num_runs_sameEnergy,targetNuc)

PURPOSE

This subroutine initializes a neutrino event on each nucleon in the realparticles vector (given to the routine). The resulting particles are set into the pertParticles vector. The reaction process is determined by the values read in by 'readInput':

The user might also choose which contributions should be included (QE, DELTA and/or higher resonances) and whether the calculation should be done for a specific experiment.

INPUTS

OUTPUT

initNeutrino/neutrino_initialized_energyFlux.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino_initialized_energyFlux.dat

PURPOSE

This file provides the incoming neutrino flux Only energies at which actual sampling has taken place are in this file The density of sampled energies reflects the flux distribution which in turn is reflected in this histogram of sampled energies.

Columns:

  • #1: Enu (in GeV)
  • #2: flux (in 1/GeV), for the calculations of X-sections relevant is only the normalized flux
  • #3: number of events
  • #4: should be ignored

init_neutrino/DoInit [ Subroutines ]

[ Top ] [ init_neutrino ] [ Subroutines ]

NAME

subroutine DoInit

PURPOSE

Doing the actual initialzation steps

initNeutrino/neutrino_absorption_cross_section.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino_absorption_cross_section.dat

PURPOSE

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

The file shows the absorption cross section for lepton ( neutrino or charged lepton) scattering for the sum of all channels which were set to TRUE in the namelist "neutrino_induced" ( QE+Delta+highRES+1pi+DIS if includeQE, includeDelta, includeRes, include1pi, includeDIS were TRUE)

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 nanobars=10^{-33}cm^2 All cross sections are given per nucleon (1/A)

Columns:

initNeutrino/neutrino_absorption_cross_section_ALL.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino_absorption_cross_section_ALL.dat

PURPOSE

More detailed information than in neutrino_absorption_cross_section.dat:

All cross sections are given per nucleon (1/A) Columns:

initNeutrino/neutrino_absorption_cross_section_QE.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino_absorption_cross_section_QE.dat

PURPOSE

The same structure as neutrino_absorption_cross_section.dat, but only for QE events (=the first interaction act was quasielastic or elastis scattering)

initNeutrino/neutrino_absorption_cross_section_2p2h.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino_absorption_cross_section_2p2h.dat

PURPOSE

The same structure as neutrino_absorption_cross_section.dat, but only for 2p2h events (=the first interaction was 2p2h)

initNeutrino/neutrino_absorption_cross_section_Delta.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino_absorption_cross_section_Delta.dat

PURPOSE

The Delta production events (=the first interaction was production of the Delta resonance)

initNeutrino/neutrino_absorption_cross_section_highRES.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino_absorption_cross_section_highRES.dat

PURPOSE

For highRES production events (=the first interaction was production of any resonance beyond Delta)

Columns:

initNeutrino/neutrino_absorption_cross_section_1pi.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino_absorption_cross_section_1pi.dat

PURPOSE

Nearly the same structure as neutrino_absorption_cross_section.dat, but only for nonresonant 1-pion production events (=the first interaction was production of 1-pion final state)

Columns:

  • #2: cross section for the channel with neutron in the final state ("final" here is "after the first interaction act", that is before final state interactions, e.g. nu n \to mu- n pi^0)
  • #3: cross section for the channel with proton in the final state, e.g. nu n \to mu- p pi^-)

initNeutrino/neutrino_absorption_cross_section_2pi.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino_absorption_cross_section_2pi.dat

PURPOSE

The same structure as neutrino_absorption_cross_section.dat, but only for 2piBG production events (=the first interaction was 2piBG)

initNeutrino/neutrino_absorption_cross_section_DIS.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino_absorption_cross_section_DIS.dat

PURPOSE

The same structure as neutrino_absorption_cross_section.dat, but only for DIS production events (=the first interaction was DIS)

initNeutrino/neutrino_absorption_cross_section_numbers.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino_absorption_cross_section_numbers.dat

PURPOSE

Shows the number of test particles which interacted (here we mean the first interaction) via various channels

Columns:

  • #1: variable which was raised (the same as neutrino_absorption_cross_section.dat) (e.g. Q^2 for nuXsectionMode=3=dSigmadQ2 mode, Elepton for nuXsectionMode=2=dSigmadQ2dElepton )
  • #2: QE channel
  • #3: Delta production
  • #4: P_{11}(1440) production
  • #5: S_{11}(1535) production
  • #6: D_{13}(1520) production
  • #7: neutron and pi+ in the final state
  • #8: proton and pi0 in the final state
  • #9: number of test particles that underwent any kind of interaction
  • #10: number of test particles (= number of tries)

init_neutrino/DoWrite [ Subroutines ]

[ Top ] [ init_neutrino ] [ Subroutines ]

NAME

subroutine DoWrite

PURPOSE

Doing the write out of most data

init_neutrino/doStoreNeutrinoInfo [ Subroutines ]

[ Top ] [ init_neutrino ] [ Subroutines ]

PURPOSE

store the neutrino and the nucleon info

init_neutrino/fillHist [ Subroutines ]

[ Top ] [ init_neutrino ] [ Subroutines ]

initNeutrino/neutrino.NuQ2planeXS.ZZZ.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino.NuQ2planeXS.ZZZ.dat, ZZZ=000 - 008

PURPOSE

This files provides absorption cross section in the nu(transferred energy)-Q^2 plane. Format is suitable as input to gnuplots "splot" command

The units are 10^{-38} cm^2/GeV^3 for process_ID=CC, NC and 10^{-33} cm^2/GeV^3 for EM All cross sections are given per nucleon (1/A) The file is not affected by cuts on outgoing lepton's angle and energy

ZZZ = ...:

  • 000: sum over all channels
  • 001: QE cross section
  • 002: Delta
  • 003: highRES
  • 004: 1pi
  • 005: DIS
  • 006: 2p2h QE
  • 007: 2p2h Delta
  • 008: 2pi background

Columns:

  • #1: nu (transfered energy), GeV
  • #2: Q^2, GeV^2
  • #3: cross section
  • #4: number of entries that lead to cross section
  • #5: should be ignored

initNeutrino/neutrino.Q2NuplaneXS.ZZZ.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino.Q2NuplaneXS.ZZZ.dat

PURPOSE

Nearly the same as neutrino.NuQ2planeXS.ZZZ.dat, but:

  • #1: Q^2, GeV^2
  • #2: nu, GeV

initNeutrino/neutrino.EprimeCostplaneXS.ZZZ.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino.EprimeCostplaneXS.ZZZ.dat, ZZZ=000 - 008

PURPOSE

This files provides absorption cross section in the Eprime(energy of the outgoing lepton)-costheta(cos of the angle between the incoming and outgoing leptons) plane. Format is suitable as input to gnuplots "splot" command

The units are 10^{-38} cm^2/GeV for process_ID=CC, NC and 10^{-33} cm^2/GeV for EM All cross sections are given per nucleon (1/A)

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

ZZZ is the origin (the first interaction vertex) of the event:

  • 000: sum over all origins
  • 001: QE cross section
  • 002: Delta
  • 003: highRES
  • 004: 1pi
  • 005: DIS
  • 006: 2p2h QE
  • 007: 2p2h Delta
  • 008: 2pi background

Columns:

  • #1: Eprime (energy of the outgoing lepton), GeV
  • #2: costheta(cos of the angle between the incoming and outgoing leptons)
  • #3: cross section
  • #4: number of entries that lead to cross section
  • #5: should be ignored

initNeutrino/neutrino.pLpTXS.ZZZ.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino.EprimepLpTXS.ZZZ.dat, ZZZ=000 - 008

PURPOSE

This files provides absorption cross section in the longitudinal momentum of the outgoing lepton - transverse momentum of the outgoing lepton plane. Format is suitable as input to gnuplots "splot" command

The units are 10^{-38} cm^2/GeV for process_ID=CC, NC and 10^{-33} cm^2/GeV for EM All cross sections are given per nucleon (1/A)

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

ZZZ is the origin (the first interaction vertex) of the event:

  • 000: sum over all origins
  • 001: QE cross section
  • 002: Delta
  • 003: highRES
  • 004: 1pi
  • 005: DIS
  • 006: 2p2h QE
  • 007: 2p2h Delta
  • 008: 2pi background

Columns:

  • #1: longitudinal momentum of the outgoing lepton, in GeV
  • #2: transverse momentum of the outgoing lepton, in GeV
  • #3: cross section
  • #4: number of entries that lead to cross section
  • #5: should be ignored

initNeutrino/neutrino.XYplaneXS.ZZZ.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino.XYplaneXS.ZZZ.dat ZZZ=000 - 008 is the origin (the first interaction vertex) of the event: 

                (see description in  neutrino.EprimeCostplaneXS.ZZZ.dat)

PURPOSE

This files provides absorption cross section in the x_{Bjorken}-y(=nu/Enu) plane. Format is suitable as input to gnuplots "splot" command

The units are 10^{-38} cm^2 for process_ID=CC, NC and 10^{-33} cm^2 for EM All cross sections are given per nucleon (1/A) The file is not affected by cuts on outgoing lepton's angle and energy

Columns:

  • #1: x_Bjorken
  • #2: y (=nu/Enu)
  • #3: cross section
  • #4: number of entries that lead to cross section
  • #5: should be ignored

initNeutrino/neutrino.NuXS.ZZZ.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino.NuXS.ZZZ.dat ZZZ=000 - 008 is the origin (the first interaction vertex) of the event: 

                (see description in  neutrino.EprimeCostplaneXS.ZZZ.dat)

PURPOSE

This files provides absorption cross section versus nu(=transfered energy) Format is suitable as input to gnuplots "splot" command The file is not affected by cuts on outgoing lepton's angle and energy

The units are 10^{-38} cm^2/GeV for process_ID=CC, NC and 10^{-33} cm^2/GeV for EM All cross sections are given per nucleon (1/A)

Columns:

  • #1: nu (transfered energy), GeV
  • #2: cross section
  • #3: number of entries that lead to cross section
  • #4: should be ignored

initNeutrino/neutrino.Q2XS.ZZZ.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino.NuQ2.ZZZ.dat

PURPOSE

Similar to neutrino.NuXS.ZZZ.dat, but versus Q2

Coulmns:

  • #1: Q2 , GeV^2

initNeutrino/neutrino.XXS.ZZZ.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino.XXS.ZZZ.dat

PURPOSE

Similar to neutrino.NuXS.ZZZ.dat, but versus x_Bjorken

Columns:

  • #1: x_Bjorken

initNeutrino/neutrino.XrecXS.ZZZ.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino.XrecXS.ZZZ.dat

PURPOSE

Similar to neutrino.NuXS.ZZZ.dat, but versus reconstructed x_Bjorken, defined for nucleon at rest: Xrec = Q^2/(2*mN*nu)

Columns:

  • #1: free x_Bjorken

initNeutrino/neutrino.WrecXS.ZZZ.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino.WrecXS.ZZZ.dat

PURPOSE

Similar to neutrino.NuXS.ZZZ.dat, but versus Wrec=qsrt(mN^2+2*mN*nu-Q2) , where mN=0.938 GeV = nucleon mass (see constants.f90), i.e. W reconstructed for nucleon at rest

Columns:

  • #1: Wrec , GeV

initNeutrino/neutrino.WfreeXS.ZZZ.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino.WfreeXS.ZZZ.dat

PURPOSE

Similar to neutrino.NuXS.ZZZ.dat, but versus Wfree, i.e. W value for boson with four-momentum q and Fermi-moving nucleon, without potential

Columns:

  • #1: Wfree , GeV

initNeutrino/neutrino.WXS.ZZZ.dat [ Output files ]

[ Top ] [ initNeutrino ] [ Output files ]

NAME

file neutrino.WXS.ZZZ.dat

PURPOSE

Similar to neutrino.NuXS.ZZZ.dat, but versus W, i.e. W value for boson with four-momentum q and Fermi-moving nucleon in the potential

Columns:

  • #1: Wfree , GeV

init_neutrino/writeInclHist [ Subroutines ]

[ Top ] [ init_neutrino ] [ Subroutines ]

initNeutrino/get_init_namelist [ Subroutines ]

[ Top ] [ initNeutrino ] [ Subroutines ]

NAME

subroutine get_init_namelist(process_ID, flavor_ID, nuXsectionMode, nuExp, debugflag, includeQE, includeDelta, includeRes, include1pi, realRun)

PURPOSE

This subroutine returns any entry of the neutrino init namelist.

OUTPUT

initNeutrino/get_runtime_vars [ Subroutines ]

[ Top ] [ initNeutrino ] [ Subroutines ]

NAME

subroutine get_runtime_vars(sigabsArrFinal,sigabsArr)

PURPOSE

This subroutine returns variables that are changed with every init.

OUTPUT

initNeutrino/isOsc [ Functions ]

[ Top ] [ initNeutrino ] [ Functions ]

NAME

logical function isOsc

PURPOSE

return true, if nuExp points to a neutrino oscllation experiment

initNeutrino/getNuMax [ Functions ]

[ Top ] [ initNeutrino ] [ Functions ]

NAME

real function getNuMax

PURPOSE

return some estimate of a "reasonable" numax