gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS

/resProd_lepton [ Modules ]

[ Top ] [ 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

real function dSigmadOmega_fdE_f_resProd_eN(eN,resID,pout,bareMass,processID)

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

  • function value -- dsigma/dOmega/dE in units of mb/(GeV sr)
  • real :: bareMass -- bare resonance mass
  • real, dimension(0:3) :: pout -- resonance momentum

resProd_lepton/sigma_pipi_res_vac [ Functions ]

[ Top ] [ resProd_lepton ] [ Functions ]

NAME

function sigma_pipi_res_vac(targetN,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:

INPUTS

  • type(particle) :: targetN -- 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(targetN,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:

INPUTS

  • type(particle) :: targetN -- 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(targetN,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:

INPUTS

  • type(particle) :: targetN -- 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(targetN,q,IDbar,IDmes) result(sigma)

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:

INPUTS

  • type(particle) :: targetN -- 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 (cf. IDmes)

resProd_lepton/dSdOmega_k_med_res [ Functions ]

[ Top ] [ resProd_lepton ] [ Functions ]

NAME

function dSdOmega_k_med_res(targetN,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) :: targetN -- 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: qPion

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: qPion

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(targetN,resID,q,baremass) result(xSection)

PURPOSE

Evaluates the cross section for gamma N -> R scattering. The return value is sigma in [mb].

Assumptions:

INPUTS

  • type(particle) :: targetN -- Target nucleon
  • real, dimension (0:3) :: q -- Virtual photon 4-momentum
  • integer :: resId -- ID of resonance

OUTPUT

  • real :: baremass -- bare mass of resonance (mass without scalar potential)
  • real :: xSection -- sigma