gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS

/barBar_Main [ Modules ]

[ Top ] [ Modules ]

NAME

module barBar_Main

barBar_Main/BBtoBYK_flag [ Global module-variables ]

[ Top ] [ barBar_Main ] [ Global module-variables ]

SOURCE

  logical, save :: BBtoBYK_flag = .true.

PURPOSE

Include BB -> B Hyperon Kaon channels. B=Nucleon^{0,1},Delta^{-,0.+,++}; Hyperon=Lambda^{0},Sigma^{0,-,+}; Kaon=K^{+,0}

barBar_Main/NNpi_BG [ Global module-variables ]

[ Top ] [ barBar_Main ] [ Global module-variables ]

SOURCE

  integer, save :: NNpi_BG = 2

PURPOSE

Switch for the N N -> N N pi background:

  • 0 = no BG
  • 1 = BG according to Teis
  • 2 = BG according to Buss
  • 3 = BG according to Weil

barBar_Main/NNomega_BG [ Global module-variables ]

[ Top ] [ barBar_Main ] [ Global module-variables ]

SOURCE

  logical, save :: NNomega_BG = .true.

PURPOSE

Incude a N N -> N N omega background term (in addition to possible resonance contributions).

barBar_Main/isofac_omega [ Global module-variables ]

[ Top ] [ barBar_Main ] [ Global module-variables ]

SOURCE

  real, save :: isofac_omega = 1.

PURPOSE

Isospin enhancement factor for p n -> p n omega, relative to p p -> p p omega. Data indicate that this is around 2, while theory predicts even larger values (up to 5). Reference: Barsov et al., EPJ A21 (2004) 521-527.

barBar_Main/XsectionBarBar [ Subroutines ]

[ Top ] [ barBar_Main ] [ Subroutines ]

NAME

subroutine XsectionBarBar (srts, teilchenIN, mediumATcollision, teilchenOUT, sigmaTot, sigmaElast, pauliIncluded, plotFlag)

PURPOSE

This routine is the main routine for baryon-baryon scattering and its cross sections. Determines total and elastic cross section and makes a Monte-Carlo decision for a special reaction channnel. This leads to a definition of ID and charge of teilchenOut, which is the final state vector.

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 information at the position of the collision
  • character(len=*), intent(in), optional :: plotFlag -- Switch on plotting of the Xsections

OUTPUT

  • type(particle), dimension(1:3), intent(out) :: teilchenOut -- outgoing particles
  • real, intent(out) :: sigmaTot -- total Xsection
  • real, intent(out) :: sigmaElast -- elastic Xsection
  • logical, intent(out) :: pauliIncluded -- true = cross section includes Pauli blocking

NOTES

plotFlag=.true. causes to make output to the files:

  • 'BaryonBaryon_Crosssection_1.dat'
  • 'BaryonBaryon_Crosssection_2.dat'

The content is explained in the files.

barBar_Main/baryonBaryon [ Namelists ]

[ Top ] [ barBar_Main ] [ Namelists ]

NAME

NAMELIST baryonBaryon

PURPOSE

Includes the switches:

barBar_Main/eta_deuteron [ Functions ]

[ Top ] [ barBar_Main ] [ Functions ]

NAME

real function eta_deuteron (srts)

PURPOSE

This routine returns the cross section for "n p -> d eta" in mb, as a function of sqrt(s). We use a tabulated spline fit of the data points, which is being read from an input file ("eta_deuteron_spline.txt"). The experimental data is taken from: H. Calen et al., PRL 79 (1997) 2642.

INPUTS

  • real, intent(in) :: srts -- sqrt(s) in the process

OUTPUT

  • returns the cross section in mb