gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS

/barBar_BarBar [ Modules ]

[ Top ] [ Modules ]

NOTES

This module includes the calculation of all baryon baryon -> baryon baryon cross sections.

barBar_BarBar/mat_NR [ Global module-variables ]

[ Top ] [ barBar_BarBar ] [ Global module-variables ]

PURPOSE

Squared matrix elements M**2/16pi for N N -> N R. See http://arxiv.org/abs/1203.3557.

SOURCE

  real, dimension(Delta:F37_1950), save :: mat_NR =  &
       (/14., & ! Delta   =02
         70., & ! P11_1440=03
          8., & ! S11_1535=04, eta/rho
          4., & ! S11_1650=05, rho
          0., & ! S11_2090=06
          4., & ! D13_1520=07, rho
          0., & ! D13_1700=08
          0., & ! D13_2080=09
         17., & ! D15_1675=10
          0., & ! G17_2190=11, omega/rho
          0., & ! P11_1710=12
          0., & ! P11_2100=13
          4., & ! P13_1720=14, rho
          0., & ! P13_1900=15, omega/rho
          4., & ! F15_1680=16, rho
          0., & ! F15_2000=17
          0., & ! F17_1990=18
          7., & ! S31_1620=19, rho
          0., & ! S31_1900=20
          7., & ! D33_1700=21, rho
          0., & ! D33_1940=22
          0., & ! D35_1930=23
          0., & ! D35_2350=24
          0., & ! P31_1750=25
         14., & ! P31_1910=26
         14., & ! P33_1600=27
          0., & ! P33_1920=28
          0., & ! F35_1750=29
          7., & ! F35_1905=30, rho
         14. /) ! F37_1950=31

barBar_BarBar/mat_DR [ Global module-variables ]

[ Top ] [ barBar_BarBar ] [ Global module-variables ]

PURPOSE

Squared matrix elements M**2/16pi for N N -> Delta R. See http://arxiv.org/abs/1203.3557.

SOURCE

  real, dimension(Delta:F37_1950), save :: mat_DR = &
       (/210., & ! Delta   =02
           0., & ! P11_1440=03
          60., & ! S11_1535=04, eta/rho
          12., & ! S11_1650=05, rho
           0., & ! S11_2090=06
          12., & ! D13_1520=07, rho
           0., & ! D13_1700=08
           0., & ! D13_2080=09
           0., & ! D15_1675=10
           0., & ! G17_2190=11, omega/rho
           0., & ! P11_1710=12
           0., & ! P11_2100=13
          12., & ! P13_1720=14, rho
           0., & ! P13_1900=15, omega/rho
          12., & ! F15_1680=16, rho
           0., & ! F15_2000=17
           0., & ! F17_1990=18
          21., & ! S31_1620=19, rho
           0., & ! S31_1900=20
          21., & ! D33_1700=21, rho
           0., & ! D33_1940=22
           0., & ! D35_1930=23
           0., & ! D35_2350=24
           0., & ! P31_1750=25
           0., & ! P31_1910=26
           0., & ! P33_1600=27
           0., & ! P33_1920=28
           0., & ! F35_1750=29
          21., & ! F35_1905=30, rho
           0. /) ! F37_1950=31

barBar_BarBar/icugnon [ Global module-variables ]

[ Top ] [ barBar_BarBar ] [ Global module-variables ]

SOURCE

  integer, save :: icugnon=1

PURPOSE

Switch for nucleon nucleon -> nucleon nucleon cross sections:

barBar_BarBar/use_ND_ND_model [ Global module-variables ]

[ Top ] [ barBar_BarBar ] [ Global module-variables ]

SOURCE

  logical, save :: use_ND_ND_model=.false.

PURPOSE

Switch for delta nucleon -> delta nucleon cross sections:

barBar_BarBar/new_NR_NR [ Global module-variables ]

[ Top ] [ barBar_BarBar ] [ Global module-variables ]

SOURCE

  logical, save :: new_NR_NR=.true.

PURPOSE

  • .false.= Switch off the NR-> NR improvement (improvement= better NN<->NN fit is being used)
  • only for debugging or comparing

barBar_BarBar/NR_NR_massSHIFT [ Global module-variables ]

[ Top ] [ barBar_BarBar ] [ Global module-variables ]

SOURCE

  logical, save :: NR_NR_massSHIFT=.false.

PURPOSE

  • .true.= Shift the srts in NR-> NR elastic collisions.

barBar_BarBar/oldOset_treatment [ Global module-variables ]

[ Top ] [ barBar_BarBar ] [ Global module-variables ]

SOURCE

  logical, save :: oldOset_treatment=.false.

PURPOSE

  • .true.= Use the old treatment for the Oset Delta width: Put everything into 3-body.
  • only for debugging or comparing

barBar_BarBar/etafac [ Global module-variables ]

[ Top ] [ barBar_BarBar ] [ Global module-variables ]

SOURCE

  real, save :: etafac = 6.5

PURPOSE

Parameter for enhancement of p n -> N*(1535) N, relative to p p -> N*(1535) N, in order to enhance eta production in pn collisions. See Calen et al., PRC 58 (1998) 2667.

barBar_BarBar/rhofac [ Global module-variables ]

[ Top ] [ barBar_BarBar ] [ Global module-variables ]

SOURCE

  real, save :: rhofac = 1.

PURPOSE

Parameter for enhancement of p n -> N*(1520) N, relative to p p -> N*(1520) N, in order to enhance rho production in p n collisions.

barBar_BarBar/neufac [ Global module-variables ]

[ Top ] [ barBar_BarBar ] [ Global module-variables ]

SOURCE

  real, save :: neufac = 1.

PURPOSE

Parameter for enhancement of p n -> N R, relative to p p -> N R, affecting all resonances.

barBar_BarBar/neufac_roper [ Global module-variables ]

[ Top ] [ barBar_BarBar ] [ Global module-variables ]

SOURCE

  real, save :: neufac_roper = 2.

PURPOSE

Parameter for enhancement of p n -> N N*(1440), relative to p p -> N N*(1440). See http://arxiv.org/abs/1203.3557.

barBar_BarBar/readInput [ Subroutines ]

[ Top ] [ barBar_BarBar ] [ Subroutines ]

NAME

subroutine readInput

PURPOSE

Reads input in jobcard out of namelist "barBar_barBar".

barBar_BarBar/barBar_barBar [ Namelists ]

[ Top ] [ barBar_BarBar ] [ Namelists ]

NAME

NAMELIST barBar_barBar

PURPOSE

Includes the switches:

barBar_BarBar/sigmaBB [ Functions ]

[ Top ] [ barBar_BarBar ] [ Functions ]

NAME

real function sigmaBB(partIN,idOut,mediumAt,srtfree_in, pauliIncluded)

NOTES

  • This function calculates all baryon-baryon -> baryon-baryon cross sections.
  • The cross sections are given in mB. Note that we sum over the channels with different charges of the final state particles.

INPUTS

  • type(particle), dimension (1:2) :: partIN -- incoming particles
  • integer, dimension (1:2) :: idOut -- Id's of outgoing particles
  • real :: srtfree_in -- SQRT(s)
  • type(medium) :: mediumAt -- medium at collision point

OUTPUT

  • logical :: pauliIncluded -- true = cross section includes Pauli blocking

NOTES

Included Xsections:

  • NN -> NN
  • NN <-> Delta Delta
  • NR <-> NR'
  • NN <-> NR

where N is a nucleon, R and R' are non-strange non-charmed baryons.

barBar_BarBar/m2d16pi [ Functions ]

[ Top ] [ barBar_BarBar ] [ Functions ]

NAME

real function m2d16pi(resID, totalCharge)

PURPOSE

Matrix Elements for NN -> NR, see: 

  * Teis et al., Pion-production in heavy-ion collisions at SIS energies,
    Z. Phys. A 356 (1997) 421-435
  * Weil et al., http://arxiv.org/abs/1203.3557.

INPUTS

  • integer, intent(in) :: resID --- resonance ID
  • integer, intent(in) :: totalCharge --- total charge of the reaction

OUTPUT

  • matrixelement**2 / 16 pi

barBar_BarBar/NR_NR [ Functions ]

[ Top ] [ barBar_BarBar ] [ Functions ]

NAME

real function NR_NR(resIN,resOut,srtFree,pinitial,idOut, isoIn, chargeIn)

PURPOSE

Evaluates NR -> NR scattering, where the resonances do not need to be identical. Assuming that it goes through a NN intermediate state.

INPUTS

  • integer :: resIn -- Id of incoming resonance
  • integer :: resOut -- Id of outgoing resonance
  • real :: pInitial -- initial momentum
  • real :: srtfree -- initial sqrt(s)
  • integer, dimension(1:2) :: idOut -- Ids of outgoing particles
  • integer, dimension(1:2) :: isoIn -- isospin of incoming particles
  • integer, dimension(1:2) :: chargeIn -- charges of incoming particles
  • real, OPTIONAL :: srtFree_withoutMass --

barBar_BarBar/nukNuk_nukNuk [ Functions ]

[ Top ] [ barBar_BarBar ] [ Functions ]

NAME

real function nukNuk_nukNuk(srtfree, massIn, totalcharge)

PURPOSE

Calculates NN -> N N scattering cross section

INPUTS

  • real :: srtfree -- sqrt(s)
  • real, dimension(1:2) :: massIn -- masses of incoming particles
  • integer :: totalcharge -- total charge of incoming particles

NOTES

The global parameter icugnon can be used to swith between two different modes of the cross section. References:

  • Cugnon old: J. Cugnon, T. Mizutani and J. Vandermeulen, Nucl. Phys. A 352, 505 (1981).
  • Cugnon new: J. Cugnon, J. Vandermeulen, D. L’hote, Nuclear Instrum. Methods B 111 (1996) 215–220.
  • PDG: L. Montanet et al., Phys. Rev. D 50 (1994) 1173–1814

barBar_BarBar/NN_NRes [ Functions ]

[ Top ] [ barBar_BarBar ] [ Functions ]

NAME

real function NN_NRes(srtFree,pInitial,idOut,isoIn,chargeIn,integerInput)

PURPOSE

Calculates NN -> N Res scattering cross section

INPUTS

  • real :: srtFree
  • real :: pInitial -- CM impulse of incoming particles
  • integer, dimension(1:2) :: idOut -- IDs of outgoing particles
  • integer, dimension(1:2) :: isoIn -- isospin of incoming particles (x2!!)
  • integer, dimension(1:2) :: chargeIn -- charges of incoming particles
  • integer, optional :: integerInput

NOTES

integerInput defines the choice of the regarded final state of the resonance:

  • integerInput=1 -- All final states allowed
  • integerInput=2 -- final state of Res is pion N
  • integerInput=3 -- final state of Res is eta N
  • integerInput=4 -- final state of Res is rho N
  • integerInput=5 -- final state of Res is omega N

Xsection is summed over all possible final state charge configurations.

barBar_BarBar/NRes_NN [ Functions ]

[ Top ] [ barBar_BarBar ] [ Functions ]

NAME

real function NRes_NN (srtFree,pInitial,idOut,isoIn,chargeIn,resonanceIn)

PURPOSE

Calculates NRes -> NN scattering cross section

INPUTS

  • real :: srtFree
  • real :: pInitial -- CM impulse of incoming particles
  • integer, dimension(1:2) :: idOut -- IDs of outgoing particles
  • integer, dimension(1:2) :: isoIn -- isospin of incoming particles (x2!!)
  • integer, dimension(1:2) :: chargeIn -- charges of incoming particles
  • integer, optional :: resonanceIn -- ID of incoming resonance

NOTES

Xsection is summed over all possible final state charge configurations.

barBar_BarBar/NN_NDelta [ Functions ]

[ Top ] [ barBar_BarBar ] [ Functions ]

NAME

real function NN_NDelta(srtFree, pInitial, idOut, isoIn, chargeIn)

PURPOSE

Calculates NN -> N Delta scattering cross section

INPUTS

  • real :: srtFree
  • real :: pInitial -- CM impulse of incoming particles
  • integer, dimension(1:2) :: idOut
  • integer, dimension(1:2) :: isoIn -- isospin of incoming particles (x2!!)
  • integer, dimension(1:2) :: chargeIn ! charges of incoming particles

NOTES

Xsection is summed over all possible final state charge configurations.

barBar_BarBar/NDelta_NN [ Functions ]

[ Top ] [ barBar_BarBar ] [ Functions ]

NAME

real function NDelta_NN(srtFree, pInitial, idOut, isoIn, chargeIn, massDelta)

PURPOSE

Calculates N Delta -> N N scattering cross section

INPUTS

  • real :: srtFree
  • real :: pInitial -- CM impulse of incoming particles
  • integer, dimension(1:2) :: idOut
  • integer, dimension(1:2) :: isoIn -- isospin of incoming particles (x2!!)
  • integer, dimension(1:2) :: chargeIn -- charges of incoming particles
  • real :: massDelta -- mass of the Delta resonance

NOTES

Xsection is summed over all possible final state charge configurations.

barBar_BarBar/NN_ResRes [ Functions ]

[ Top ] [ barBar_BarBar ] [ Functions ]

NAME

real function NN_ResRes(srtFree,pInitial,idOut,isoIn,chargeIn)

PURPOSE

Calculates NN -> Res Res scattering cross section, assuming constant matrix elements.

INPUTS

  • real :: srtFree
  • real :: pInitial -- CM impulse of incoming particles
  • integer, dimension(1:2) :: idOut -- IDs of outgoing particles
  • integer, dimension(1:2) :: isoIn -- isospin of incoming particles (x2!!)
  • integer, dimension(1:2) :: chargeIn -- charges of incoming particles

NOTES

The cross section is summed over all possible final state charge configurations. Currently two channels are implemented: Delta-Delta and Delta-N*(1535). For Delta-Delta see Effenberger, PhD, appendix A.1.2. The Delta-N*(1535) channel was fitted to roughly match Pythia's pi-eta production XS around sqrts ~ 3 GeV, in order to describe HADES dilepton data for p+p@3.5GeV.

barBar_BarBar/DRes_NN [ Functions ]

[ Top ] [ barBar_BarBar ] [ Functions ]

NAME

real function DRes_NN(srtFree, pInitial, idOut, isoIn, chargeIn, resonanceIN)

PURPOSE

Calculates Delta Res -> N N scattering cross section.

INPUTS

  • real :: srtFree
  • real :: pInitial -- CM impulse of incoming particles
  • integer, dimension(1:2) :: idOut -- IDs of outgoing particles
  • integer, dimension(1:2) :: isoIN -- isospin of incoming particles (x2!!)
  • integer, dimension(1:2) :: chargeIN -- charges of incoming particles
  • integer, optional :: resonanceIN -- ID of incoming resonance

NOTES

Xsection is summed over all possible final state charge configurations.

barBar_BarBar/iifac0 [ Functions ]

[ Top ] [ barBar_BarBar ] [ Functions ]

NAME

real function iifac0(IDout,is,iz)

NOTES

This function calculates the isospin-factors given in table 3.4 of M.E.'s diploma thesis + factor for identical particles in the final state. Sum over final state charges!

INPUTS

  • integer, dimension(1:2) :: IDout -- IDs of final state
  • integer, dimension(1:2) :: is -- isospins of initial state times 2
  • integer, dimension(1:2) :: iz -- charges of initial state particles

barBar_BarBar/iifac1 [ Functions ]

[ Top ] [ barBar_BarBar ] [ Functions ]

NAME

real function iifac1 (IDout,is,iz,iz3v,iz4v)

NOTES

This function calculates the isospin-factors given by eq. (A.5) in Effenberger's PhD thesis + factor for identical particles in the final state. No sum over final state charges!

INPUTS

  • integer, dimension(1:2) :: IDout -- IDs of final state
  • integer, dimension(1:2) :: is -- isospins of initial state times 2
  • integer, dimension(1:2) :: iz -- charges of initial state particles
  • integer :: iz3v,iz4v -- charges of final state particles

barBar_BarBar/chooseCharge [ Functions ]

[ Top ] [ barBar_BarBar ] [ Functions ]

NAME

function chooseCharge (partIn, idOut) result(chargeOut)

PURPOSE

Two baryons with defined charges and ID's are defined as incoming particles (partIN), and we also pass the outgoing particles to this routine (idOUT). This subroutine then makes a Monte-Carlo decision for the charges of the outgoing particles, using the isospin factors provided by iifac or the NDelta -> NDelta model.

barBar_BarBar/delta2Body_inMedium_treatment [ Subroutines ]

[ Top ] [ barBar_BarBar ] [ Subroutines ]

NAME

logical function delta2Body_inMedium_treatment()

PURPOSE

This function decides whether the condition is fulfilled such that the in-medium width of the delta can be partially implemented via two body processes.