gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS

/barBar_BarHypKaon [ Modules ]

[ Top ] [ Modules ]

NAME

module barBar_BarHypKaon

PURPOSE

Administrates the calculation of cross sections for strangeness production from BB->BYK collisions (with B=N,Delta and Y=Lambda,Sigma,Lambda*,Sigma*). The properties of the final-state channels are calculated as well.

NOTES

  • All isospin channels are included.
  • Ref.: K. Tsushima et al., PRC59 (1999) 369.
  • The original Tsushima model was extended and updated in arXiv:1404.7011.
  • Useful for reactions near and above threshold (up to sqrt(s)=3.6 GeV), not appropriate for sqrt(s) far below threshold (cross section does not fit well existing data & perturbative treatment of kaon production necessary).

barBar_BarHypKaon/enable [ Global module-variables ]

[ Top ] [ barBar_BarHypKaon ] [ Global module-variables ]

SOURCE

  logical, save :: enable = .true.

PURPOSE

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

barBar_BarHypKaon/parameter_set [ Global module-variables ]

[ Top ] [ barBar_BarHypKaon ] [ Global module-variables ]

SOURCE

  integer, save :: parameter_set = 2

PURPOSE

Select a particular parameter set for BB->BYK collisions. Possible values:

  • 1 = original Tsushima model: Tsushima et al., PRC59 (1999) 369
  • 2 = extended/adjusted model, fitted to HADES data: Agakishiev et al., arXiv:1404.7011
  • 3 = custom parameters based on Tsushima values (as given by the array 'a' in the jobcard; those values not given in the jobcard are adopted from Tsushima, i.e. parameter set 1)
  • 4 = custom parameters based on HADES values (as given by the array 'a' in the jobcard; those values not given in the jobcard are adopted from HADES, i.e. parameter set 2)

barBar_BarHypKaon/a [ Global module-variables ]

[ Top ] [ barBar_BarHypKaon ] [ Global module-variables ]

SOURCE

  real, dimension(1:Nch), save :: a = -1.

PURPOSE

This array contains the "a" parameters (in microbarn) for the 30 primary channels, see:

  • Tsushima et al., PRC59 (1999) 369, table III
  • Agakishiev et al., arXiv:1404.7011, chapter IV

Note: The values given in the jobcard are only used for parameter_set = 3 and 4.

barBar_BarHypKaon/BB_BYK [ Namelists ]

[ Top ] [ barBar_BarHypKaon ] [ Namelists ]

NAME

NAMELIST BB_BYK

PURPOSE

Includes the switches:

barBar_BarHypKaon/barBar_barBarMeson_strange [ Functions ]

[ Top ] [ barBar_BarHypKaon ] [ Functions ]

NAME

function barBar_barBarMeson_strange(srts, partIN) result (sigma_BYK)

PURPOSE

XSections for strangeness production from BB->BYK collisions.

INPUTS

  • real :: srts -- sqrt(s) of the incoming channel
  • type(particle),dimension(1:2) :: partIN -- the incoming particles

OUTPUT

  • real, dimension(1:18) :: sigma_BYK -- cross section for each channel (mb)

NOTES

  • This routine calculates also the properties of the finalstates, needed for the construction of the Type(PreEvent) variable.
  • B stands for nucleon or a \Delta or a N* resonance. Other resonances not included. N* resonances are trated as nucleons.
  • In RMF mode sqrt(s*) already contains the self energy of the incoming channel. See its definition in generateFinalState in master_2Body.f90-module variable "srtS_XS".

barBar_BarHypKaon/xs_param [ Functions ]

[ Top ] [ barBar_BarHypKaon ] [ Functions ]

NAME

real function xs_param (s, j)

PURPOSE

Calculate the parametrized cross section for the primary channels (without isospin factor).

INPUTS

  • real, intent(in) :: s --- Mandelstam s in GeV**2
  • integer, intent(in) :: j --- channel number

OUTPUT

  • Cross section in mb.

barBar_BarHypKaon/get_Channels_BYK [ Functions ]

[ Top ] [ barBar_BarHypKaon ] [ Functions ]

NAME

function get_Channels_BYK (InChan) result (parts)

PURPOSE

Defines the elements of the type(PreEvent) function

INPUTS

integer, intent(in) :: InChan -- Number of final channel

OUTPUT

type(preEvent), dimension(1:3) :: parts

barBar_BarHypKaon/Get_Charge [ Functions ]

[ Top ] [ barBar_BarHypKaon ] [ Functions ]

NAME

function Get_Charge

PURPOSE

Convert the variables (between GiBUU & Munich-RBUU) related to the isospin states.

INPUTS

integer :: direction : 1--GiBUU-->RBUU / 2--RBUU-->GiBUU. integer :: species : particle ID. integer :: iso_in : charge of incoming particle.

NOTES

Only for nucleons, Delta's, N*'s, hyperons (Lambda,Sigma) and kaons (K^{+,0}).

barBar_BarHypKaon/Get_ID [ Functions ]

[ Top ] [ barBar_BarHypKaon ] [ Functions ]

NAME

Integer Function Get_ID (species, iso)

PURPOSE

Convert identification number (ID) between GiBUU and Munich-RBUU.

INPUTS

integer :: species : particle ID. integer :: iso_in : charge of incoming particle.

NOTES

Only for nucleons, Delta's, N*'s, hyperons (Lambda,Sigma) and kaons (K^{+,0}).

barBar_BarHypKaon/get_ChannelParameters [ Subroutines ]

[ Top ] [ barBar_BarHypKaon ] [ Subroutines ]

NAME

subroutine get_ChannelParameters

PURPOSE

Returns the number of final states & their properties.

INPUTS

OUTPUT

  • integer :: nexit -- Number of final states for BB-->BYK^{+}
  • integer :: iexit -- Relation to table (see above)
  • integer :: fac -- isospin factor for each final state in BYK^{+}
  • integer :: isohf -- isospin-state of hyperon Y in BYK^{+}
  • integer :: isobf -- isospin-state of the baryon B in BYK^{+}
  • integer :: nexit0 -- Number of final states for BB-->BYK^{0}
  • integer :: iexit0 -- Relation to table (see above)
  • integer :: fac0 -- isospin factor for each final state in BYK^{0}
  • integer :: isohf0 -- isospin-state of hyperon Y in BYK^{0}
  • integer :: isobf0 -- isospin-state of the baryon B in BYK^{0}

NOTES

  • Notation: B : p,n,\DELTA^{-},\DELTA^{0},\DELTA^{+},\DELTA^{++} Y : \Lambda^{0},\Sigma^{-},\Sigma^{0},\Sigma^{+} K : Kaon^{+},Kaon^{0}
  • This is a partial copy of the Munich-RBUU code. For this reason we use here - and only here - local variables for the different isospin-states of the considered particle species, which may differ from those used in the GiBUU code. Please, DO NOT USE THE FOLLOWING PRESCRIPTIONS OUTSIDE OF THIS ROUTINE!!!
  • Relation between local variables and baryon/meson isospin states, used in this routine:
  • PROTONS = 1
  • NEUTRONS = 0
  • \DELTA^{-} = 10
  • \DELTA^{0} = 11
  • \DELTA^{+} = 12
  • \DELTA^{++} = 13
  • N*^{0} = 20
  • N*^{+} = 21
  • \Lambda^{0} = 0
  • \Sigma^{-} = 1
  • \Sigma^{0} = 2
  • \Sigma^{+} = 3
  • Kaon^{+} = 1
  • Kaon^{0} = 0