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

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

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/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

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

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

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

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) =the For process_ID=EM the units are nanobars=10^{-33}cm^2

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:

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