TABLE OF CONTENTS
- 1. /lepton2p2h
- 1.1. lepton2p2h/ME_Version
- 1.2. lepton2p2h/ME_Norm_QE
- 1.3. lepton2p2h/ME_Norm_Delta
- 1.4. lepton2p2h/ME_Mass_QE
- 1.5. lepton2p2h/ME_Mass_Delta
- 1.6. lepton2p2h/ME_Transversity
- 1.7. lepton2p2h/ME_LONG
- 1.8. lepton2p2h/ME_W3
- 1.9. lepton2p2h/ME_ODW
- 1.10. lepton2p2h/inmedW
- 1.11. lepton2p2h/T
- 1.12. lepton2p2h/readInput
- 1.13. lepton2p2h/Lepton2p2h
- 1.14. lepton2p2h/lepton2p2h_DoQE
- 1.15. lepton2p2h/lepton2p2h_DoDelta
- 1.16. lepton2p2h/lepton2p2h_SelectN2
- 1.17. lepton2p2h/lepton2p2h_FinalState
- 1.18. lepton2p2h/lepton2p2h_XS
- 1.18.1. lepton2p2h_XS/ME_Dipole_transverse
- 1.18.2. lepton2p2h_XS/ME_W1W2W3
/lepton2p2h [ Modules ]
NAME
module lepton2p2h
PURPOSE
Do all the internals for 2p2h scattering:
EM:
- ell N1 N2 --> ell' N1' N2' == gamma* N1 N2 --> N1' N2'
- ell N1 N2 --> ell' N Delta == gamma* N1 N2 --> N Delta
NC:
CC:
- nu N1 N2 --> ell- N1' N2' (sum of hadronic charges increases by +1)
- nu N1 N2 --> ell- N Delta ( -- " -- )
antiEM, antiNC and antiCC are the same as EM, NC, CC.
cases 1 - 3 give parametrizations for 2p2h part of structure function W1 in terms of Q^2, no distinction for neutrinos and antineutrinos Cases 1 and 2 are those in: Lalakulich Gallmeister Mosel PRC86(2012)014614 Case 3 gives a reasonable description of MiniBooNE dd neutrino data cases 4 - 5 give parametrization for MEC part of W1 from Christy and Bosted In this case also W3 is related to W1 (acc. Martini and Ericsson) Case 4 describes double-differential data from MiniBooNE for neutrino and antineutrino scattering. It also describes the dd inclusive Xsection for neutrinos from T2K.
lepton2p2h/ME_Version [ Global module-variables ]
[ Top ] [ lepton2p2h ] [ Global module-variables ]
PURPOSE
indicate the type of matrix element parametrisation
SOURCE
integer, save :: ME_Version = 4 ! ! 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
lepton2p2h/ME_Norm_QE [ Global module-variables ]
[ Top ] [ lepton2p2h ] [ 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
lepton2p2h/ME_Norm_Delta [ Global module-variables ]
[ Top ] [ lepton2p2h ] [ 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
lepton2p2h/ME_Mass_QE [ Global module-variables ]
[ Top ] [ lepton2p2h ] [ 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
lepton2p2h/ME_Mass_Delta [ Global module-variables ]
[ Top ] [ lepton2p2h ] [ 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
lepton2p2h/ME_Transversity [ Global module-variables ]
[ Top ] [ lepton2p2h ] [ 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 turns structure function W2 on so that 2p2h is pure transverse
lepton2p2h/ME_LONG [ Global module-variables ]
[ Top ] [ lepton2p2h ] [ 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 longitudinal contribution to structure funct. W2 off
lepton2p2h/ME_W3 [ Global module-variables ]
[ Top ] [ lepton2p2h ] [ 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
lepton2p2h/ME_ODW [ Global module-variables ]
[ Top ] [ lepton2p2h ] [ 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 = 1
NOTES
- O'Connell et al: PR C6 (1972) 719
- Martini et al: PR C80 (2009) 065501
lepton2p2h/inmedW [ Global module-variables ]
[ Top ] [ lepton2p2h ] [ Global module-variables ]
PURPOSE
Controls which inv mass W is used in parametrization of 2p2h W1
SOURCE
integer, save :: inmedW = 2
NOTES
- inmedW = 1 : W = free, static inv. mass in 2p2h parametrization of W1
- inmedW = 2 : W = inv mass for Fermi moving nucleons in potential
- inmedW = 3 : W = inv mass for Fermi moving nucleons without potential
lepton2p2h/T [ Global module-variables ]
[ Top ] [ lepton2p2h ] [ Global module-variables ]
PURPOSE
Isospin of target nucleus, factor in Eq. (A10) in O'Connell, Donnelly, Walecka, PR C6 (1972) 719
SOURCE
real, save :: T = 1
NOTES
This is isospin T of NN pair, for 2p2h processes
lepton2p2h/readInput [ Subroutines ]
[ Top ] [ lepton2p2h ] [ Subroutines ]
NAME
subroutine readInput
lepton2p2h/Lepton2p2h [ Namelists ]
[ Top ] [ lepton2p2h ] [ 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
- ME_ODW
- inmedW
- T
lepton2p2h/lepton2p2h_DoQE [ Subroutines ]
[ Top ] [ lepton2p2h ] [ Subroutines ]
NAME
subroutine lepton2p2h_DoQE(eN,outPart,XS)
PURPOSE
Do all the electron induced 2p2h-QE scattering gamma* N1 N2 -> N1' N2'
INPUTS
- type(electronNucleon_event) :: eN -- electron-Nucleon event info
OUTPUT
- type(particle), dimension(:) :: OutPart -- the two produced nucleons
- real :: XS -- the cross section
lepton2p2h/lepton2p2h_DoDelta [ Subroutines ]
[ Top ] [ lepton2p2h ] [ Subroutines ]
NAME
subroutine lepton2p2h_DoDelta(eN,outPart,XS)
PURPOSE
Do all the electron induced 2p2h-QE scattering gamma* N1 N2 -> N Delta
INPUTS
- type(electronNucleon_event) :: eN -- electron-Nucleon event info
OUTPUT
lepton2p2h/lepton2p2h_SelectN2 [ Subroutines ]
[ Top ] [ lepton2p2h ] [ Subroutines ]
NAME
subroutine lepton2p2h_SelectN2(eN)
PURPOSE
Finds the second nucleon for the 2p2h collision
INPUTS
- type(electronNucleon_event) :: eN -- electron-Nucleon event info
OUTPUT
- type(electronNucleon_event) :: eN -- a second nucleon is added
NOTES
- The seond particle is generated analytically, not by selecting a testparticle from the real particle vector.
- This is at a very basic level. You may add more sophisticated features as eq. two-particle correlatione etc.
- A threshold check Wfree>(2*mN+1MeV) is performed
lepton2p2h/lepton2p2h_FinalState [ Subroutines ]
[ Top ] [ lepton2p2h ] [ Subroutines ]
NAME
subroutine lepton2p2h_FinalState(eN,outPart,DoQE,flagOK)
PURPOSE
Generate the final state of the electron 2p2h event
lepton2p2h/lepton2p2h_XS [ Functions ]
[ Top ] [ lepton2p2h ] [ Functions ]
NAME
real function lepton2p2h_XS(eN,outPart,DoQE)
PURPOSE
calculate the electron induced 2p2h-QE cross section
INPUTS
- type(electronNucleon_event) :: eN -- electron-Nucleon event info
- type(particle),dimension(:) :: OutPart -- the outgoing particles
- logical :: DoQE -- .true. for NN final state, .false. for N Delta
OUTPUT
- the function value
NOTES
- One has to give a realistic parametrization of the matrix element
- If one randomly selects the position of the second particle, one has to respect this in the XS calculation (and not multiply it with the density at the position)
lepton2p2h_XS/ME_Dipole_transverse [ Functions ]
[ Top ] [ lepton2p2h_XS ] [ Functions ]
NAME
real function ME_Dipole_transverse(eN)
PURPOSE
calculate the 2p2h matrix element according to W_1(g_munu -q_um q_nu /Q2) * L^munu so that the contribution is only to the transverse part
NOTES
You have full access to all incoming and outgoing particles:
- eN%lepton_in -- incoming lepton
- eN%nucleon -- incoming nucleon 1
- eN%nucleon2 -- incoming nucleon 2
exchanged boson:
- eN%boson -- exchanged boson
even without considering the final state particles, you know the kind of process via 'eN%idProcess', which may take the values EM,NC,CC and also antiEM,antiNC,antiCC
lepton2p2h_XS/ME_W1W2W3 [ Functions ]
[ Top ] [ lepton2p2h_XS ] [ Functions ]
NAME
real function ME_W1W2W3(eN)
PURPOSE
to calculate the 2p2h contribution to the inclusive cross sections for electrons and neutrinos, cross section depends on all 3 structure functs
NOTES
You have full access to all incoming and outgoing particles:
- eN%lepton_in -- incoming lepton
- eN%nucleon -- incoming nucleon 1
- eN%nucleon2 -- incoming nucleon 2
exchanged boson:
- eN%boson -- exchanged boson
even without considering the final state particles, you know the kind of process via 'eN%idProcess', which may take the values EM,NC,CC and also antiEM,antiNC,antiCC