TABLE OF CONTENTS
- 1. neutrinoParms/Lepton2p2h
- 2. neutrinoParms/twoptwoh_readInput
- 3. neutrinoParms/nl_neweN
- 4. neutrinoParms/nl_Neutrino2piBack
- 5. neutrinoParms/neutrino_readInput
- 6. neutrinoParms/new_eNres
- 7. neutrinoParms/new_eN
- 8. neutrinoParms/Adep
- 9. neutrinoParms/T
- 10. neutrinoParms/inmedW
- 11. neutrinoParms/ME_ODW
- 12. neutrinoParms/ME_W3
- 13. neutrinoParms/ME_LONG
- 14. neutrinoParms/ME_Transversity
- 15. neutrinoParms/ME_Mass_Delta
- 16. neutrinoParms/ME_Mass_QE
- 17. neutrinoParms/ME_Norm_Delta
- 18. neutrinoParms/ME_Norm_QE
- 19. neutrinoParms/ME_Version
- 20. neutrinoParms/A00
- 21. neutrinoParms/Q2cut
- 22. neutrinoParms/normCB
- 23. neutrinoParms/normRes
- 24. neutrinoParms/NormpiBG
- 25. neutrinoParms/Wtrans
- 26. neutrinoParms/twopiBG_on
neutrinoParms/twopiBG_on [ Global module-variables ]
[ Top ] [ Global module-variables ]
SOURCE
real, save:: 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:: Wtrans = 2.7 ! GeV
PURPOSE
W for transition from Bosted-Christie Parametrization to PYTHIA DIS
neutrinoParms/NormpiBG [ Global module-variables ]
[ Top ] [ Global module-variables ]
SOURCE
real, save:: 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/normCB [ Global module-variables ]
[ Top ] [ Global module-variables ]
SOURCE
real, save:: normCB = 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:: Q2cut = 4.5 ! GeV^2
PURPOSE
threshold value
neutrinoParms/A00 [ Global module-variables ]
[ Top ] [ Global module-variables ]
SOURCE
real, save:: A00 = -0.43 real, save:: A10 = -0.97 real, save:: A01 = 0.87 real, save:: A20 = 8.01 real, save:: A11 = 4.59 real, save:: A02 = 3.18 real, save:: A30 = -10.42 real, save:: A21 = 2.04 real, save:: A12 = -9.14 real, save:: A03 = -15.15 real, save:: W00 = 1.55 real, save:: 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 :: 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 :: 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 :: 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 :: 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 :: 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 :: 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 :: 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 :: ME_W3 = (/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 :: 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 :: 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 :: T = 0
NOTES
- T = 0, 1 , ...
- T = 99 gives T = (N-Z)/2
neutrinoParms/Adep [ Global module-variables ]
[ Top ] [ Global module-variables ]
PURPOSE
Switch for A-dependence of 2p2h structure function
SOURCE
integer, save :: 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 :: 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 :: new_eNres = .true.
PURPOSE
switch to choose new electron-Nucleon X-sections for resonances with Bosted-Christy parametrization
neutrinoParms/neutrino_readInput [ Subroutines ]
[ Top ] [ Subroutines ]
NAME
subroutine neutrino_readInput
neutrinoParms/nl_Neutrino2piBack [ Namelists ]
NAME
NAMELIST /nl_Neutrino2piBack/
PURPOSE
Parameters for transition from RES to DIS:
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 ]
NAME
NAMELIST /nl_neweN/
PURPOSE
Parameters for new eN implementation * new_eN * new_eNres * ME_ODW * T
neutrinoParms/twoptwoh_readInput [ Subroutines ]
[ Top ] [ Subroutines ]
NAME
subroutine twoptwoh_readInput
neutrinoParms/Lepton2p2h [ Namelists ]
NAME
NAMELIST /Lepton2p2h/
PURPOSE
Includes parameters for 2p2h events: