### 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