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 handles the calculation of various neutrino cross sections, depending on the choice of the user for each entry in the particle vector and finally, it sets the produced output in 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_dSigmadQsdElepton
  • 13 = EXP_dSigmadQs
  • 14 = EXP_dSigmadCosTheta
  • 15 = EXP_dSigmadElepton
  • 16 = EXP_dSigmaMC
  • 17 = EXP_dSigmadW


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 neutrino in neutrino mode
  • 16 = LBNF/DUNE antineutrino 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
  • 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/storeNucleon [ Global module-variables ]

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

SOURCE

  integer, save :: storeNucleon = 1

PURPOSE

indicate which kind of struck nucleon to save:

  • 1: free Nuecleon (i.e. potential and Fermi momentum 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

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

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

stepsize of min. 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

stepsize of 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

stepsize of energy of outgoing leptons, used in 2D dsigma/dEdcos(theta)


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

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

SOURCE

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

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:

parameters that control the binning in the 2D cross sections d2sigma/(dE dcostheta) E limits for lepton energies: Elept_min = minimal value of lepton energy Elept_max = maximal value of lepton energy delta_Elept = bin width in Elept cos(theta) = cost limits for lepton angles: cost_min = minimal value of cost cost_max = maximal value of cost delta_cost = bin width in cost


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

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


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)

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)

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.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)

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 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/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_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_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=dSigmadQs mode, Elepton for nuXsectionMode=2=dSigmadQsdElepton )
  • #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


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