./code/collisions/twoBodyReactions/baryonBaryon/barBar_barBar.f90

1. /barBar_BarBar

1.1. barBar_BarBar/mat_NR
1.2. barBar_BarBar/mat_DR
1.3. barBar_BarBar/icugnon
1.4. barBar_BarBar/use_ND_ND_model
1.5. barBar_BarBar/new_NR_NR
1.6. barBar_BarBar/NR_NR_massSHIFT
1.7. barBar_BarBar/oldOset_treatment
1.8. barBar_BarBar/etafac
1.9. barBar_BarBar/rhofac
1.10. barBar_BarBar/neufac
1.11. barBar_BarBar/neufac_roper
1.12. barBar_BarBar/deltaN_densityDependence
1.13. barBar_BarBar/alpha
1.14. barBar_BarBar/readInput
1.15. barBar_BarBar/barBar_barBar
1.16. barBar_BarBar/sigmaBB
1.17. barBar_BarBar/m2d16pi
1.18. barBar_BarBar/NR_NR
1.19. barBar_BarBar/nukNuk_nukNuk
1.20. barBar_BarBar/NN_NRes
1.21. barBar_BarBar/NRes_NN
1.22. barBar_BarBar/NN_NDelta
1.23. barBar_BarBar/NDelta_NN
1.24. barBar_BarBar/NN_ResRes
1.25. barBar_BarBar/DRes_NN
1.26. barBar_BarBar/iifac0
1.27. barBar_BarBar/iifac1
1.28. barBar_BarBar/chooseCharge
1.29. barBar_BarBar/delta2Body_inMedium_treatment

/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:

0=old parametrization
1=new parametrization (Alexej Larionov, Cugnon)

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:

false=old parametrization
true =one pion exchange model (Effenberger, Buss)

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/deltaN_densityDependence [ Global module-variables ]

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

SOURCE

  logical, save :: deltaN_densityDependence = .false.

PURPOSE

Switch to turn on the density dependence of the NN <-> N Delta cross section. We use the density dependence from: Song/Ko, arXiv:1403.7363. The strength of the in-medium modification is controlled by the parameter alpha.

barBar_BarBar/alpha [ Global module-variables ]

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

SOURCE

  real, save :: alpha = 1.2

PURPOSE

Parameter which controls the density dependence of the NN <-> N Delta cross section. We use the density dependence from: Song/Ko, arXiv:1403.7363. See also deltaN_densityDependence.

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:

mat_NR
mat_DR
icugnon
use_ND_ND_model
new_NR_NR
NR_NR_massSHIFT
oldOset_treatment
etafac
rhofac
neufac
neufac_roper
deltaN_densityDependence
alpha

barBar_BarBar/sigmaBB [ Functions ]

[ Top ] [ barBar_BarBar ] [ Functions ]

NAME

real function sigmaBB(teilchenIN,idOut,mediumAt,srtfree)

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), intent(in):: teilchenIN -- incoming particles
integer, dimension (1:2), intent(in) :: idOut -- Id's of outgoing particles
real, intent(in) :: srtfree -- SQRT(s)
type(medium),intent(in) :: mediumAt -- mediuma at the collision point

OUTPUT

logical, intent(out) :: 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

real,intent(in) :: pInitial, srtfree ! Initial momentum and sqrt(s)
integer , intent(in) :: resIn, resOut ! Incoming and outgoing resonance ID
integer, intent(in), dimension(1:2) :: idOut,isoIn, chargeIn ! outgoing particles, incoming isospin, incoming charges

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, intent(in) :: srtfree
real, dimension(1:2), intent(in) :: massIn ! masses of incoming particles
integer, dimension (1:2), intent(in):: totalcharge ! total charge of incoming particles

NOTES

The 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, intent(in) :: srtFree
real, intent(in) :: pInitial ! CM impulse of incoming particles
integer, intent(in), dimension(1:2) :: idOut ! IDs of outgoing particles
integer, intent(in), dimension(1:2) :: isoIN ! isospin of incoming particles (times 2!!)
integer, intent(in), dimension(1:2) :: chargeIN ! charges of incoming particles
integer, intent(in), 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, intent(in) :: srtFree
real, intent(in) :: pInitial ! CM impulse of incoming particles
integer, intent(in), dimension(1:2) :: idOut ! IDs of outgoing particles
integer, intent(in), dimension(1:2) :: isoIN ! isospin of incoming particles (times 2!!)
integer, intent(in), dimension(1:2) :: chargeIN ! charges of incoming particles
integer, intent(in), 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, mediumAt)

PURPOSE

Calculates NN -> N Delta scattering cross section

INPUTS

real, intent(in) :: srtFree
real, intent(in) :: pInitial ! CM impulse of incoming particles
integer, intent(in), dimension(1:2) :: idOut
integer, intent(in), dimension(1:2) :: isoIN ! isospin of incoming particles (times 2!!)
integer, intent(in), dimension(1:2) :: chargeIN ! charges of incoming particles
type(medium), intent(in) :: mediumAt

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, mediumAt, massDelta)

PURPOSE

Calculates N Delta -> N N scattering cross section

INPUTS

real, intent(in) :: srtFree
real, intent(in) :: pInitial ! CM impulse of incoming particles
integer, intent(in), dimension(1:2) :: idOut
integer, intent(in), dimension(1:2) :: isoIN ! isospin of incoming particles (times 2!!)
integer, intent(in), dimension(1:2) :: chargeIN ! charges of incoming particles
type(medium), intent(in) :: mediumAt
real, intent(in) :: 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, intent(in) :: srtFree
real, intent(in) :: pInitial ! CM impulse of incoming particles
integer, intent(in), dimension(1:2) :: idOut ! IDs of outgoing particles
integer, intent(in), dimension(1:2) :: isoIN ! isospin of incoming particles (times 2!!)
integer, intent(in), 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, intent(in) :: srtFree
real, intent(in) :: pInitial ! CM impulse of incoming particles
integer, intent(in), dimension(1:2) :: idOut ! IDs of outgoing particles
integer, intent(in), dimension(1:2) :: isoIN ! isospin of incoming particles (times 2!!)
integer, intent(in), dimension(1:2) :: chargeIN ! charges of incoming particles
integer, intent(in), 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 my diploma thesis (M.E.) + factor for identical particles in the final state. Sum over final state charges!

INPUTS

IDout : IDs of final state
is : isospins of initial state times 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

IDout : IDs of final state
is : isospins of initial state times 2
iz : charges of initial state particles
iz3v,iz4v : charges of final state particles

barBar_BarBar/chooseCharge [ Functions ]

[ Top ] [ barBar_BarBar ] [ Functions ]

NAME

function chooseCharge (teilchenIN, idOut) result(chargeOut)

PURPOSE

Two baryons with defined charges and ID's are defined as incoming particles (teilchenIN), 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.

TABLE OF CONTENTS

/barBar_BarBar [ Modules ]

barBar_BarBar/mat_NR [ Global module-variables ]

barBar_BarBar/mat_DR [ Global module-variables ]

barBar_BarBar/icugnon [ Global module-variables ]

barBar_BarBar/use_ND_ND_model [ Global module-variables ]

barBar_BarBar/new_NR_NR [ Global module-variables ]

barBar_BarBar/NR_NR_massSHIFT [ Global module-variables ]

barBar_BarBar/oldOset_treatment [ Global module-variables ]

barBar_BarBar/etafac [ Global module-variables ]

barBar_BarBar/rhofac [ Global module-variables ]

barBar_BarBar/neufac [ Global module-variables ]

barBar_BarBar/neufac_roper [ Global module-variables ]

barBar_BarBar/deltaN_densityDependence [ Global module-variables ]

barBar_BarBar/alpha [ Global module-variables ]

barBar_BarBar/readInput [ Subroutines ]

barBar_BarBar/barBar_barBar [ Namelists ]

barBar_BarBar/sigmaBB [ Functions ]

barBar_BarBar/m2d16pi [ Functions ]

barBar_BarBar/NR_NR [ Functions ]

barBar_BarBar/nukNuk_nukNuk [ Functions ]

barBar_BarBar/NN_NRes [ Functions ]

barBar_BarBar/NRes_NN [ Functions ]

barBar_BarBar/NN_NDelta [ Functions ]

barBar_BarBar/NDelta_NN [ Functions ]

barBar_BarBar/NN_ResRes [ Functions ]

barBar_BarBar/DRes_NN [ Functions ]

barBar_BarBar/iifac0 [ Functions ]

barBar_BarBar/iifac1 [ Functions ]

barBar_BarBar/chooseCharge [ Functions ]

barBar_BarBar/delta2Body_inMedium_treatment [ Subroutines ]