GiBUU documentation

Table of Contents Sourcefiles Index Functions Global module-variables Modules Namelists Output files Program Subroutines Types

TABLE OF CONTENTS

• 1. neutrinoParms/VAfact
• 2. neutrinoParms/setDefault_Res
• 3. neutrinoParms/Lepton2p2h
• 4. neutrinoParms/readInput_2p2h
• 5. neutrinoParms/nl_neweN
• 6. neutrinoParms/nl_Neutrino2piBack
• 7. neutrinoParms/readInput_neutrino
• 8. neutrinoParms/new_eNres
• 9. neutrinoParms/new_eN
• 10. neutrinoParms/Adep
• 11. neutrinoParms/T
• 12. neutrinoParms/inmedW
• 13. neutrinoParms/ME_ODW
• 14. neutrinoParms/ME_W3
• 15. neutrinoParms/ME_LONG
• 16. neutrinoParms/ME_Transversity
• 17. neutrinoParms/ME_Mass_Delta
• 18. neutrinoParms/ME_Mass_QE
• 19. neutrinoParms/ME_Norm_Delta
• 20. neutrinoParms/ME_Norm_QE
• 21. neutrinoParms/ME_Version
• 22. neutrinoParms/A00
• 23. neutrinoParms/Q2cut
• 24. neutrinoParms/normBC
• 25. neutrinoParms/normRes
• 26. neutrinoParms/NormpiBG
• 27. neutrinoParms/twopiBG_on

---

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

SOURCE

  real, save, public :: NormRes = 1.0

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.
• 3: ???

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

• 1: W = static inv. mass in 2p2h parametrization of W1
• 2: W = inv mass for Fermi moving nucleons in potential
• 3: W = inv mass for Fermi moving nucleons without potential

---

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:

• NormpiBG
• normRES
• normBC

Includes parameters of neutrino 2 pion background:

• Q2cut
• A00
• A10
• A01
• A02
• A11
• A20
• A03
• A12
• A21
• A30
• W00
• x00
• twopiBG_on

---

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:

• ME_Version
• ME_Norm_QE
• ME_Norm_Delta
• ME_Mass_QE
• ME_Mass_Delta
• ME_Transversity
• ME_LONG
• ME_W3
• inmedW
• Adep

---

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.

---

neutrinoParms/VAfact [ Subroutines ]

[ Top ] [ Subroutines ]

NAME

subroutine VAfact(ME_ODW,qvecsq,nu,kinfact)

PURPOSE

???

NOTES

• should be converted into a function (even pure?)
• ME_ODW does not need to be a parameter, since the module variable could be used

Generated from ./code/init/neutrino/neutrinoParms.f90 with ROBODoc V4.99.42 on Thu Nov 21 2024 11:32:54