gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS


neutrinoParms/twopiBG_on [ Global module-variables ]

[ Top ] [ Global module-variables ]

SOURCE

  real, save, public :: twopiBG_on = 1.5 ! GeV

PURPOSE

W for transition from 1piBG to 2piBG for BC parametrization of nr X-section


neutrinoParms/Wtrans [ Global module-variables ]

[ Top ] [ Global module-variables ]

SOURCE

  real, save, public :: Wtrans = 2.7 ! GeV

PURPOSE

W for transition from Bosted-Christy Parametrization to PYTHIA DIS


neutrinoParms/NormpiBG [ Global module-variables ]

[ Top ] [ Global module-variables ]

SOURCE

  real, save, public :: NormpiBG = 1.0

PURPOSE

overall normalization factor for pi BG and Bloom-Gilman X-section, only relevant for neutrinos


neutrinoParms/normRes [ Global module-variables ]

[ Top ] [ Global module-variables ]

PURPOSE

overall normalization factor for neutrino-induced resonance contributions beyond the Delta only relevant for neutrinos


neutrinoParms/normBC [ Global module-variables ]

[ Top ] [ Global module-variables ]

SOURCE

  real, save, public :: normBC = 1.0

PURPOSE

overall normalization factor for neutrino-induced Christy-Bosted contributions between 2 GeV and DIS onset only relevant for neutrinos


neutrinoParms/Q2cut [ Global module-variables ]

[ Top ] [ Global module-variables ]

SOURCE

  real, save, public ::  Q2cut = 4.5 ! GeV^2

PURPOSE

threshold value


neutrinoParms/A00 [ Global module-variables ]

[ Top ] [ Global module-variables ]

SOURCE

  real, save, public ::  A00 = -0.43
  real, save, public ::  A10 = -0.97
  real, save, public ::  A01 = 0.87
  real, save, public ::  A20 = 8.01
  real, save, public ::  A11 = 4.59
  real, save, public ::  A02 = 3.18
  real, save, public ::  A30 = -10.42
  real, save, public ::  A21 = 2.04
  real, save, public ::  A12 = -9.14
  real, save, public ::  A03 = -15.15

  real, save, public ::  W00 = 1.55
  real, save, public ::  x00 = 0.35

PURPOSE

coefficients in Taylor expansion of attenuation function Att(W,x)


neutrinoParms/ME_Version [ Global module-variables ]

[ Top ] [ Global module-variables ]

PURPOSE

indicate the type of matrix element parametrisation

SOURCE

  integer, save, public :: ME_Version = 6

NOTES

possible values:

  • 1: const ME_Norm_XX ! const for CC fitted to MiniBooNE is 1.8e-6
  • 2: constant transverse and decreasing with Enu
  • 3: "Dipole transverse" transverse, fall with Q2 as 4-th power
  • 4: MEC from E. Christy (8/2015), with parametrization for longitudinal
  • 5: MEC from Bosted arXiV:1203.2262, with parametrization for longitudinal
  • 6: MEC additional parametrization, with parametrization for longitudinal not yet implemented

remarks:

  • case 1 is model-I in Lalakulich,Gallmeister,Mosel PRC86(2012)014614
  • case 2 is model-II from Lalakulich,Gallmeister,Mosel PRC86(2012)014614
  • case 3 gives a good description of MiniBooNE data with MA ~ 1.5 GeV


neutrinoParms/ME_Norm_QE [ Global module-variables ]

[ Top ] [ Global module-variables ]

PURPOSE

Overall strength of 2p2h matrix element with 2N out

for (EM,CC,NC)

SOURCE

  real,dimension(1:3), save, public :: ME_Norm_QE    = (/1.0, 1.0, 1.0/)

NOTES

The value == 1 gives the coded strength


neutrinoParms/ME_Norm_Delta [ Global module-variables ]

[ Top ] [ Global module-variables ]

PURPOSE

Overall strength of 2p2h matrix element with NDelta out

for (EM,CC,NC)

SOURCE

  real,dimension(1:3), save, public :: ME_Norm_Delta = (/1.0, 1.0, 1.0/)

NOTES

The value == 1 is a dummy value


neutrinoParms/ME_Mass_QE [ Global module-variables ]

[ Top ] [ Global module-variables ]

PURPOSE

Cutoff-mass in some parametrizations of 2p2h matrix element for NN out

for (EM,CC,NC)

SOURCE

  real,dimension(1:3), save, public :: ME_Mass_QE    = (/1.0, 1.0, 1.0/)

NOTES

The value == 1 is a dummy value


neutrinoParms/ME_Mass_Delta [ Global module-variables ]

[ Top ] [ Global module-variables ]

PURPOSE

Cutoff-mass in some parametrizations of matrix element for NDelta out

for (EM,CC,NC)

SOURCE

  real,dimension(1:3), save, public :: ME_Mass_Delta = (/1.0, 1.0, 1.0/)

NOTES

The value == 1 is a dummy value


neutrinoParms/ME_Transversity [ Global module-variables ]

[ Top ] [ Global module-variables ]

PURPOSE

Parametrisation of structure functions

for (EM,CC,NC)

SOURCE

  real,dimension(1:3), save, public :: ME_Transversity = (/1.0, 1.0, 1.0/)

NOTES

The value = 1 chooses structure function W2 so that 2p2h is pure transverse


neutrinoParms/ME_LONG [ Global module-variables ]

[ Top ] [ Global module-variables ]

PURPOSE

Parametrization of structure functions

for (EM,CC,NC)

SOURCE

  real,dimension(1:3), save, public :: ME_LONG = (/0.0, 0.0, 0.0/)

NOTES

The value = 0 turns any additional longitudinal contribution to structure funct. W2 off


neutrinoParms/ME_W3 [ Global module-variables ]

[ Top ] [ Global module-variables ]

PURPOSE

Overall strength factor for structure function W3

only for (CC,NC)

SOURCE

  real,dimension(1:3), save, public :: ME_W3 = (/0.0, 1.0, 1.0/)

NOTES

overall strength parameter for structure function W3


neutrinoParms/ME_ODW [ Global module-variables ]

[ Top ] [ Global module-variables ]

PURPOSE

switch for choosing the connection between structure functions W1(electron) and W1(neutrino) and W3(neutrino):

  • 1: for expressions from Martini et al.
  • 2: for expressions from O'Connell et al.

only for (CC,NC)

SOURCE

  integer, save, public :: ME_ODW = 2

NOTES

  • O'Connell et al: PR C6 (1972) 719
  • Martini et al: PR C80 (2009) 065501


neutrinoParms/inmedW [ Global module-variables ]

[ Top ] [ Global module-variables ]

PURPOSE

Controls which inv mass W is used in parametrization of 2p2h W1

SOURCE

  integer, save, public :: inmedW = 1

NOTES


neutrinoParms/T [ Global module-variables ]

[ Top ] [ Global module-variables ]

PURPOSE

target isospin, affects only neutrino 2p2h structure function

SOURCE

  real, save, public ::  T = 1

NOTES

  • T = 0, 1 , ...
  • T = 99 gives T = abs((N-Z)/2)


neutrinoParms/Adep [ Global module-variables ]

[ Top ] [ Global module-variables ]

PURPOSE

Switch for A-dependence of 2p2h structure function

SOURCE

  integer, save, public ::  Adep = 2

NOTES

  • 1: A-dependence for zero-range force (Mosel, Gallmeister, 2016)
  • 2: linear A-dependence, normalized to C12


neutrinoParms/new_eN [ Global module-variables ]

[ Top ] [ Global module-variables ]

SOURCE

  logical, save, public :: new_eN = .true.

PURPOSE

switch to choose new electron-Nucleon X-sections with Bosted-Christy parametrization


neutrinoParms/new_eNres [ Global module-variables ]

[ Top ] [ Global module-variables ]

SOURCE

  logical, save, public :: new_eNres = .true.

PURPOSE

switch to choose new electron-Nucleon X-sections for resonances with Bosted-Christy parametrization

NOTES

  • only used if new_eN = .T.
  • currently, .T. should be used for electrons and .F. for neutrinos
  • .T. and .F. can be set as the default value using 'setDefault_Res' giving the process_ID (=EM,CC,NC and anti) as argument
  • any value given in the jobcard will override this default value


neutrinoParms/readInput_neutrino [ Subroutines ]

[ Top ] [ Subroutines ]

NAME

subroutine readInput_neutrino


neutrinoParms/nl_Neutrino2piBack [ Namelists ]

[ Top ] [ Namelists ]

NAME

NAMELIST /nl_Neutrino2piBack/

PURPOSE

Parameters for transition from RES to DIS:

Includes parameters of neutrino 2 pion background:


neutrinoParms/nl_neweN [ Namelists ]

[ Top ] [ Namelists ]

NAME

NAMELIST /nl_neweN/

PURPOSE

Parameters for new eN implementation * new_eN * new_eNres * ME_ODW * T


neutrinoParms/readInput_2p2h [ Subroutines ]

[ Top ] [ Subroutines ]

NAME

subroutine readInput_2p2h


neutrinoParms/Lepton2p2h [ Namelists ]

[ Top ] [ Namelists ]

NAME

NAMELIST /Lepton2p2h/

PURPOSE

Includes parameters for 2p2h events:


neutrinoParms/setDefault_Res [ Subroutines ]

[ Top ] [ Subroutines ]

NAME

subroutine setDefault_Res(process_ID)

PURPOSE

select the default resonance treatment for 'new_eN' according the used process.