TABLE OF CONTENTS
- 1. /sigmaNucleon
- 1.1. sigmaNucleon/sigmaNuc
- 1.1.1. sigmaNuc/makeOutput
/sigmaNucleon [ Modules ]
NAME
module sigmaNucleon
PURPOSE
Includes the cross sections for sigma-nucleon scattering in the resonance regime. Public routines:
sigmaNucleon/sigmaNuc [ Subroutines ]
[ Top ] [ sigmaNucleon ] [ Subroutines ]
NAME
subroutine sigmaNuc(srts,teilchenIN,mediumATcollision,momentumLRF,teilchenOUT,sigmaTot,sigmaElast,useHiEnergy,HiEnergySchwelle,plotFlag)
PURPOSE
Evaluates sigma Nucleon -> anything cross sections and returns also a "preevent"
INPUTS
- real, intent(in) :: srts ! sqrt(s) in the process
- type(particle),dimension(1:2), intent(in) :: teilchenIn ! colliding particles
- type(medium), intent(in) :: mediumATcollision ! Medium informations at the position of the collision
- real, intent(in) ,dimension(0:3) :: momentumLRF ! Total Momentum in LRF
High energy matching:
- logical,intent(in) :: useHiEnergy
- .true. if High-Energy cross sections are given by paramBarMesHE
- real,intent(in) :: HiEnergySchwelle
- threshold sqrt(s) for paramBarMesHE, i.e. at which energy the cross sections of paramBarMesHE are used
Debugging:
- logical, intent(in),optional :: plotFlag ! Switch on plotting of the Xsections
RESULT
- real, intent(out) :: sigmaTot ! total Xsection
- real, intent(out) :: 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), intent(out) :: teilchenOut ! colliding particles
The cross sections are based upon a parametrization by Golubeva. See routine golub in parametrizationBarMes.
NOTES
Possible final states are :
sigmaNuc/makeOutput [ Subroutines ]
[ Top ] [ sigmaNuc ] [ Subroutines ]
NAME
subroutine makeOutput
PURPOSE
Writes all cross sections to file as function of srts and plab [GeV]. Filenames:
- 'sigmaN_sigTotElast.dat' : sigmaTot, sigmaElast
- 'sigmaN_resProd.dat' : Baryon resonance production
- 'sigmaN_nonStrange_nuk.dat' : non-strange meson with nucleon in final state