gibuu is hosted by Hepforge, IPPP Durham
GiBUU

## /twoBodyStatistics [ Modules ]

[ Top ] [ Modules ]

NAME

module twoBodyStatistics

PURPOSE

This module contains routines for collecting 2-body statistics, e.g. sqrt(s) distributions and collision rates.

## twoBodyStatistics/flag_sqrts [ Global module-variables ]

PURPOSE

If .true., then the calculation and output of the sqrts distributions from subroutine sqrts_distribution will be done

SOURCE

logical, save :: flag_sqrts=.false.

## twoBodyStatistics/flag_rate [ Global module-variables ]

PURPOSE

If .true., then the calculation and output of the collision rates from subroutine rate will be done

SOURCE

logical, save :: flag_rate=.false.

## twoBodyStatistics/flag_varirate [ Global module-variables ]

PURPOSE

If .true., then the calculation and output of the collision rates from subroutine varirate will be done

SOURCE

logical, save :: flag_varirate=.false.

## twoBodyStatistics/varirate_chargeZero [ Global module-variables ]

PURPOSE

If this flag to .true., then all charge states are combined together

SOURCE

logical, save :: varirate_chargeZero=.true.

## twoBodyStatistics/varirate_size [ Global module-variables ]

PURPOSE

size of array to hold all rates

SOURCE

integer, save :: varirate_size=100

## twoBodyStatistics/sqrts_mode [ Global module-variables ]

PURPOSE

This flag determines the way how sqrt(s) is calculated (if flag_sqrts = .true.). 1 = use vacuum sqrt(s) 2 = use in-medium, i.e. full sqrts

SOURCE

integer, save ::  sqrts_mode = 1

## twoBodyStatistics/init [ Subroutines ]

[ Top ] [ twoBodyStatistics ] [ Subroutines ]

NAME

subroutine init

PURPOSE

INPUTS

• (none)

OUTPUT

• Initializes global module variables

## twoBodyStatistics/ColStat [ Namelists ]

[ Top ] [ twoBodyStatistics ] [ Namelists ]

NAME

NAMELIST ColStat

PURPOSE

Namelist which includes the input switches:

## twoBodyStatistics/sqrts_distribution [ Subroutines ]

[ Top ] [ twoBodyStatistics ] [ Subroutines ]

NAME

subroutine sqrts_distribution (pair, itype, flag)

PURPOSE

Computes the distribution of colliding pairs of particles over sqrt(s).

INPUTS

• type(particle), dimension(1:2), intent(in) :: pair -- incoming pair of particles
• integer, intent(in) :: itype
• logical, optional, intent(in) :: flag -- .true. do output

Possible values for itype:

• 1 -- separate collision,
• 2 -- collision in presence of a particle nearby before sqrt(s) correction,
• 3 -- collision in presence of a particle nearby after sqrt(s) correction,
• 4 -- elastic binary collisions
• 5 -- inelastic binary collisions

OUTPUT

"dNdsqrts_Bm.dat" and "dNdsqrts_mm.dat".

## twoBodyStatistics/dNdsqrts_BB.dat [ Output files ]

[ Top ] [ twoBodyStatistics ] [ Output files ]

NAME

file dNdsqrts_BB.dat

PURPOSE

Contains the sqrt(s) distribution of baryon-baryon collisions. Enabled by flag_sqrts.

## twoBodyStatistics/dNdsqrts_Bm.dat [ Output files ]

[ Top ] [ twoBodyStatistics ] [ Output files ]

NAME

file dNdsqrts_Bm.dat

PURPOSE

Contains the sqrt(s) distribution of baryon-meson collisions. Enabled by flag_sqrts.

## twoBodyStatistics/dNdsqrts_mm.dat [ Output files ]

[ Top ] [ twoBodyStatistics ] [ Output files ]

NAME

file dNdsqrts_mm.dat

PURPOSE

Contains the sqrt(s) distribution of meson-meson collisions. Enabled by flag_sqrts.

## twoBodyStatistics/rate [ Subroutines ]

[ Top ] [ twoBodyStatistics ] [ Subroutines ]

NAME

subroutine rate(teilchenIn,teilchenOut,time,flag)

PURPOSE

Computes the collision rates for the different types of collisions.

INPUTS

• type(particle), intent(in), dimension(:) :: teilchenIn -- incoming particles
• type(particle), intent(in), dimension(:) :: teilchenOut -- outgoing particles
• real, intent(in) :: time -- current time step
• logical, optional, intent(in) :: flag -- .true. do output of the accumulated collision numbers

OUTPUT

• file "rate.dat"

## twoBodyStatistics/varirate [ Subroutines ]

[ Top ] [ twoBodyStatistics ] [ Subroutines ]

NAME

subroutine varirate(partIn,partOut,time)

PURPOSE

Computes the collision rates for alltypes of collisions.

Here the list of input and output channels are generated on the fly, thus this routine calculates the "variate-rate", or short: "varirate".

INPUTS

• type(particle), intent(in), dimension(:) :: partIn -- incoming particles
• type(particle), intent(in), dimension(:) :: partOut -- outgoing particles
• real, intent(in) :: time -- current time step

OUTPUT

• file "VariRate.particles.dat", "VariRate.rates.dat"