TABLE OF CONTENTS
- 1. /pionDelta_resonance
- 1.1. pionDelta_resonance/pionDelta
- 1.1.1. pionDelta/evaluateXsections
- 1.1.2. pionDelta/makeDecision
- 1.1.3. pionDelta/makeOutput
/pionDelta_resonance [ Modules ]
NAME
module pionDelta_resonance
PURPOSE
Includes the cross sections for pion-delta scattering in the resonance regime Implemented are the following reactions:
* pion delta -> X
Public routines:
* pionDelta
NOTE pi Delta -> K Kbar N not yet implemented
pionDelta_resonance/pionDelta [ Modules ]
[ Top ] [ pionDelta_resonance ] [ Modules ]
NAME
subroutine pionDelta(srts,teilchenIN,mediumATcollision,momentumLRF,teilchenOUT,sigmaTot,sigmaElast,plotFlag)
PURPOSE
Evaluates pion Delta -> anything cross sections and returns also a "preevent"
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
NOTES
Possible final states are :
NOTE pi Delta -> K Kbar N not yet implemented
pionDelta/evaluateXsections [ Subroutines ]
[ Top ] [ pionDelta ] [ Subroutines ]
NAME
subroutine evaluateXsections
PURPOSE
Evaluates all Xsections for pi N scattering
pionDelta/makeDecision [ Subroutines ]
[ Top ] [ pionDelta ] [ Subroutines ]
NAME
subroutine MakeDecision
PURPOSE
Decides on the final state which is returned via teilchenOut by Monte-Carlo.
* Assigns charges and ID's. * Only for resonance-production also the mass is assigned, since the mass of the resonance needed to be calculated earlier.
The Monte-Carlo routine is adding up channels until the sum is exceeding x*sigma(total). x has a flat distribution in [0,1]. The last added channel is then the one which is chosen for the event. After choosing a channel, the subroutine is returning to the calling routine.
pionDelta/makeOutput [ Subroutines ]
[ Top ] [ pionDelta ] [ Subroutines ]
NAME
subroutine makeOutput
PURPOSE
Writes all cross sections to file as function of srts and plab [GeV]. Filenames:
- 'piDelta_sigTotElast.dat' : sigmaTot, sigmaElast
- 'piDelta_strangeProd.dat' : strangeness production
- 'piDelta_resProd.dat' : Baryon resonance production