TABLE OF CONTENTS
- 1. /resProd_lepton
- 1.1. resProd_lepton/debug
- 1.2. resProd_lepton/dSigmadOmega_fdE_f_resProd_eN
- 1.3. resProd_lepton/sigma_pipi_res_vac
- 1.4. resProd_lepton/sigma_pipi_res
- 1.5. resProd_lepton/sigma_barMes_res_vac
- 1.6. resProd_lepton/sigma_barMes_res
- 1.7. resProd_lepton/dSdOmega_k_med_res
- 1.8. resProd_lepton/dSdO_fdE_fdO_k_med_res_EN
- 1.9. resProd_lepton/dOmegaCM_dOmega
- 1.10. resProd_lepton/sigma_resProd
/resProd_lepton [ Modules ]
NAME
module resProd_lepton
PURPOSE
- Evaluates cross sections for gamma N -> R
- For details see the notes about this in the work of Oliver Buss
resProd_lepton/debug [ Global module-variables ]
[ Top ] [ resProd_lepton ] [ Global module-variables ]
SOURCE
logical, parameter :: debug=.false.
PURPOSE
resProd_lepton/dSigmadOmega_fdE_f_resProd_eN [ Functions ]
[ Top ] [ resProd_lepton ] [ Functions ]
NAME
function dSigmadOmega_fdE_f_resProd_eN(eN,resID,pout, bareMass,processID) result(xSection)
PURPOSE
- Evaluates cross sections for e N -> R
- For details see the notes about this in the work of Oliver Buss
INPUTS
- type(electronNucleon_event) :: eN -- underlying electron nucleon event
- integer :: resID -- resonance ID
- integer, optional :: processID
OUTPUT
- real :: bareMass -- bare resonance mass
- real, dimension(0:3) :: pout -- resonance momentum
- real :: xSection -- dsigma/dOmega/dE in units of mb/(GeV sr)
NOTES
- Enhances dSigmadOmega_fdE_f_resProd by allowing arbitrary electron momentum directions.
resProd_lepton/sigma_pipi_res_vac [ Functions ]
[ Top ] [ resProd_lepton ] [ Functions ]
NAME
function sigma_pipi_res_vac(targetNuc,q) result(sigma)
PURPOSE
Evaluates the resonance contribution to double-pion production in gamma N -> R -> N 2Pi scattering. The return value is sigma in [mb]. Converts target nucleon first to vacuum nucleon!!! dsigma/dOmega_electron/dE_electron/dOmega_pion
Assumptions:
- No interferences among resonances.
INPUTS
- type(particle) :: targetNuc -- Target nucleon
- real, dimension (0:3) :: q -- Virtual photon 4-momentum
OUTPUT
- sigma(0:3) Cross sections for gamma+nucleon->nucleon+2Pi production
- sigma(1) -> nucleon piMinus piPlus
- sigma(2) -> nucleon piPlus piNull or nucleon piMinus piNull
- sigma(3) -> nucleon piNull piNull
- sigma(0) Total Xsection into nucleon+2 Pions
resProd_lepton/sigma_pipi_res [ Functions ]
[ Top ] [ resProd_lepton ] [ Functions ]
NAME
function sigma_pipi_res(targetNuc,q) result(sigma)
PURPOSE
Evaluates the resonance contribution to double-pion production in gamma N -> R -> N 2Pi scattering. The return value is sigma in [mb].
Assumptions:
- No interferences among resonances.
INPUTS
- type(particle) :: targetNuc -- Target nucleon
- real, dimension (0:3) :: q -- Virtual photon 4-momentum
OUTPUT
- sigma(0:3) Cross sections for gamma+nucleon->nucleon+2Pi production
- sigma(1) -> nucleon piMinus piPlus
- sigma(2) -> nucleon piPlus piNull or nucleon piMinus piNull
- sigma(3) -> nucleon piNull piNull
- sigma(0) Total Xsection into nucleon+2 Pions
resProd_lepton/sigma_barMes_res_vac [ Functions ]
[ Top ] [ resProd_lepton ] [ Functions ]
NAME
function sigma_barMes_res_vac(targetNuc,q,IDbar,IDmes) result (sigma)
PURPOSE
Evaluates the resonance contribution of gamma N -> R -> B m^0 scattering (where X may be a nucleon or Delta, while m^0 is a neutral meson). The return value is sigma in [mb]. Converts target nucleon first to vacuum nucleon!!!
Assumptions:
- No interferences among resonances.
INPUTS
- type(particle) :: targetNuc -- Target nucleon
- real, dimension (0:3) :: q -- Virtual photon 4-momentum
- integer :: IDbar -- ID of produced baryon N (nucleon or Delta)
- integer :: IDmes -- array containing the IDs of produced mesons m
OUTPUT
- sigma -- Cross sections for gamma N -> R -> B m^0 production, for all mesons you asked for
resProd_lepton/sigma_barMes_res [ Functions ]
[ Top ] [ resProd_lepton ] [ Functions ]
NAME
function sigma_barMes_res(targetNuc,q,IDbar,IDmes) result(sigma_VM)
PURPOSE
Evaluates the resonance contribution of gamma N -> R -> B m^0 scattering (where B may be a nucleon or Delta, while m^0 is a neutral meson). The return value is sigma in [mb]. The cross section is calculated separately for all mesons which are passed in IDmes.
Assumptions:
- No interferences among resonances.
INPUTS
- type(particle) :: targetNuc -- Target nucleon
- real, dimension (0:3) :: q -- Virtual photon 4-momentum
- integer :: IDbar -- ID of produced baryon B (nucleon or Delta)
- integer :: IDmes -- array containing the IDs of produced mesons m
OUTPUT
- sigma -- Cross sections for gamma N -> R -> B m^0 production, for all mesons you asked for
resProd_lepton/dSdOmega_k_med_res [ Functions ]
[ Top ] [ resProd_lepton ] [ Functions ]
NAME
function dSdOmega_k_med_res(targetNuc,q,k,pf) result(sigma_dOmega)
PURPOSE
Evaluates the resonance contribution to pion production in gamma R->eNPi scattering. The return value is dsigma/dOmega(pion) in [mub/Sr].
Assumptions:
- No interferences among resonances.
- Isotropic decay of the resonance in its rest-frame.
- In the vacuum.
INPUTS
- type(particle) :: targetNuc -- Target nucleon
- real, dimension (0:3) :: q -- Virtual photon 4-momentum
- real, dimension (0:3) :: k -- pion 4-momentum
- real, dimension (0:3) :: pf -- Outgoing nucleon 4-momentum
OUTPUT
- logical :: success -- flag
- real,dimension(-1:1) :: sigma_dOmega -- dsigma/dOmega_pion; Index: PionCharge
resProd_lepton/dSdO_fdE_fdO_k_med_res_EN [ Functions ]
[ Top ] [ resProd_lepton ] [ Functions ]
NAME
function dSdO_fdE_fdO_k_med_res_EN(eN,k,pf,processID) result(sigma_dOmega)
PURPOSE
Evaluates the resonance contribution to pion production in eN->eR->eNPi scattering. The return value is dsigma/dOmega(electron)/dE(electron)/dOmega(pion) in [mb/GeV/Sr**2].
Assumptions:
- No interferences among resonances.
- Isotropic decay of the resonance in its rest-frame.
- In the vacuum.
INPUTS
- type(electronNucleon_event) :: eN -- electron-nucleon scattering event
- real, dimension (0:3) :: pf -- Outgoing nucleon 4-momentum
- real, dimension (0:3) :: k -- Outgoing pion 4-momentum
- integer, optional :: processID -- See module leptonicID for usage
- integer, optional :: pionNucleonSystem -- If this parameter is set to 1, then we evaluate dOmega_pion in the calculation frame. If it's 2 then it is evaluated in the cm frame of the outgoing pion and nucleon.
OUTPUT
- real,dimension(-1:1):: sigma_dOmega -- dsigma/dOmega_electron/dE_electron/dOmega_pion; Index: PionCharge
NOTES
- Enhances dSdO_fdE_fdO_k_med_res by allowing arbitrary electron momentum directions
resProd_lepton/dOmegaCM_dOmega [ Functions ]
[ Top ] [ resProd_lepton ] [ Functions ]
NAME
real function dOmegaCM_dOmega()
PURPOSE
Evaluates the Jacobian for dOmega_CM(pion)/dOmega_lab(pion)
resProd_lepton/sigma_resProd [ Functions ]
[ Top ] [ resProd_lepton ] [ Functions ]
NAME
function sigma_resProd(targetNucleon,resID,q,baremass) result(xSection)
PURPOSE
Evaluates the cross section for gamma N -> R scattering. The return value is sigma in [mb].
Assumptions:
- No interferences among resonances.
INPUTS
- type(particle) :: targetNuc -- Target nucleon
- real, dimension (0:3) :: q -- Virtual photon 4-momentum
- integer :: resId -- ID of resonance
OUTPUT