TABLE OF CONTENTS
- 1. /omegaNucleon
- 1.1. omegaNucleon/omegaNuc
- 1.2. omegaNucleon/makeOutput
/omegaNucleon [ Modules ]
NAME
module omegaNucleon
PURPOSE
Includes the cross sections for omega-nucleon scattering in the resonance regime.
omegaNucleon/omegaNuc [ Subroutines ]
[ Top ] [ omegaNucleon ] [ Subroutines ]
NAME
subroutine omegaNuc(srts,teilchenIN,mediumATcollision,momentumLRF, teilchenOUT,sigmaTot,sigmaElast,useHiEnergy,HiEnergySchwelle,plotFlag, sigmaArr)
PURPOSE
Evaluates omega Nucleon -> anything cross sections and returns also a "preevent"
INPUTS
- real :: srts --- sqrt(s) in the process
- type(particle),dimension(1:2) :: teilchenIn --- colliding particles
- type(medium) :: mediumATcollision --- Medium informations at the position of the collision
- real, dimension(0:3) :: momentumLRF --- Total Momentum in LRF
High energy matching:
- logical :: useHiEnergy --- .true. if High-Energy cross sections are given by paramBarMesHE
- real :: HiEnergySchwelle --- threshold sqrt(s) for paramBarMesHE, i.e. at which energy the cross sections of paramBarMesHE are used
Debugging:
- logical,optional :: plotFlag --- Switch on plotting of the Xsections
RESULT
- real :: sigmaTot --- total Xsection
- real :: sigmaElast --- elastic Xsection
This routine does a Monte-Carlo-decision according to the partial cross sections to decide on a final state with maximal 3 final state particles. These are returned in the vector teilchenOut. The kinematics of these teilchen is only fixed in the case of a single produced resonance. Otherwise the kinematics still need to be established. The result is:
- type(preEvent),dimension(1:3) :: teilchenOut --- particles
- real, dimension(6), optional :: sigmaArr -- partial cross sections
The cross sections are based upon a parametrization by Golubeva. See routine golub_omega in parametrizationBarMes.
NOTES
Possible final states are :
- 1-particle : baryon resonances
- 2-particle : pi N, omega N, K Lambda, K Sigma
- 3-particle : pi pi N
omegaNucleon/makeOutput [ Subroutines ]
[ Top ] [ omegaNucleon ] [ Subroutines ]
NAME
subroutine makeOutput
PURPOSE
Writes all cross sections to file as function of srts and plab [GeV] . Filenames:
- 'omegaN_sigTotElast.dat' : sigmaTot, sigmaElast
- 'omegaN_resProd.dat' : Baryon production (Resonances with ID's 1:40)
- 'omegaN_nonStrange_nuk.dat' : non-strange meson with nucleon in final state
- 'omegaN_strangeProd.dat' : Kaon and hyperon in final state