gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS


/pionDelta_resonance [ Modules ]

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