gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS

/parametrizationsBarMes [ Modules ]

[ Top ] [ Modules ]

NAME

module parametrizationsBarMes

PURPOSE

Includes routines which are parametrizations of "baryon meson -> X" data.

parametrizationsBarMes/piN_elastic [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

real function piN_elastic(charge,srts,success)

PURPOSE

Returns the cross section for elastic pion proton scattering.

These are own brute force fits to data for plab=0..20 GeV, i.e. for sqrt(s)<6.2 GeV !!!

INPUTS

OUTPUT

  • logical :: success -- .true. if result is useful
  • funcion value :: sigma(pi N -> pi N) in mb

NOTES

  • Fitted to data in the range of plab=0..20 GeV (only pi- and pi+ !!).
  • Brute force fits to data

parametrizationsBarMes/piN_chargeExchange [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

real function piN_chargeExchange(srts)

PURPOSE

Returns the cross section for pion nucleon charge exchange, i.e.:

  • p pi- <-> n pi0
  • n pi+ <-> p pi0

Fitted to data in the range of plab=0..1.4 GeV (i.e. sqrt(s)<1.9 GeV)

INPUTS

  • integer :: srts -- sqrt(s) in GeV

OUTPUT

  • function value -- charge exchange Xsection in mb

NOTES

See buuinput/pionNucleon/chargeExchange/fit.nb

parametrizationsBarMes/pionquasibg [ Subroutines ]

[ Top ] [ parametrizationsBarMes ] [ Subroutines ]

NAME

subroutine pionquasibg(nukcharge,pioncharge,sqrts,elastic,chargeEx)

PURPOSE

Subroutine evaluates quasielastic background for low energy pi N->pi N

INPUTS

  • integer,intent(in) :: nukcharge,pioncharge --- charges
  • real,intent(in) ::sqrts --- SQRT(s)

OUTPUT

  • real,intent(out) :: elastic, chargeEx --background crosssection in mB, elastic and charge exchange

parametrizationsBarMes/sigelkp [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

real function sigelkp(p,n)

PURPOSE

Calculate:

  • K+ p -> K+ p for n=1
  • K+ n -> K+ n for n=2

INPUTS

  • real :: p -- kaon momentum in the rest frame of the proton
  • integer :: n -- tells which cross section is calculated (see above)

parametrizationsBarMes/sigCEXkp [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

real function sigCEXkp(srts)

PURPOSE

Calculate * K+ n -> K0 p

INPUTS

  • real :: srts -- invariant energy (GeV)

parametrizationsBarMes/sigCEXk0 [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

real function sigCEXk0(srts)

PURPOSE

Calculate * K0 p -> K+ n

INPUTS

  • real :: srts -- invariant energy (GeV)

parametrizationsBarMes/kppbg [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

real function kppbg(plab2)

PURPOSE

Calculate * K+ p -> p K pion

INPUTS

  • real :: plab2

parametrizationsBarMes/kpnbg [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

real function kpnbg(plab2)

PURPOSE

Calculate:

  • K0 p -> n K pion

INPUTS

  • real :: plab2

parametrizationsBarMes/kaonbg [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

real function kaonbg(srts,ichin,ichout,isospinc)

PURPOSE

Calculate K- p and K- n cross section im mb

INPUTS

  • real, intent(in) :: srts -- invariant energy (GeV)
  • integer, intent(in) :: ichin -- initial channel
  • integer, intent(in) :: ichout -- final channel
  • integer, intent(in) :: isospinc -- switch for isospin of incoming nucleon: 0 -> K- p, 1 -> K- n

NOTES

Enumeration of incoming and outgoing channels is as follows:

  • 1 --- K- p or K- n
  • 2 --- Kbar0 n
  • 3 --- Lambda pi0
  • 4 --- Sigma+ pi- or Sigma- pi0 (Sigma0 pi-)
  • 5 --- Sigma- pi+
  • 6 --- Sigma0 pi0

parametrizationsBarMes/kaonbgt [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

real function kaonbgt(plab2)

PURPOSE

Calculate K- p total cross section im mb

parametrizationsBarMes/huang [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

function huang(srts) result(sig)

PURPOSE

cross sections for pi p -> K Sigma in mb:

  • sig(1) = pi^{+} p -> K^{+} Sigma+
  • sig(2) = pi^{0} p -> K^{+} Sigma0 or pi^{-} p -> K^{0} Sigma0
  • sig(3) = pi^{-} p -> K^{0} Sigma0
  • sig(4) = pi^{-} p -> K^{+} Sigma-

NOTES

See Effenberger PHD, page 236

parametrizationsBarMes/huanglam [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

real function huanglam(srts)

PURPOSE

cross section for pi- p -> Lambda K0 in mb

parametrizationsBarMes/huanglamd [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

real function huanglamd(srts)

PURPOSE

cross section for pi- Delta++ -> Lambda K+ in mb

Clebsch for (pi^- D++ -> I=1/2 I_z=1/2) is 1/2

Therefore : Xsection( pi- D++ -> Lambda k+)=1/2 Xsection( pion Delta -> Lambda kaon)

=> Xsection( pion Delta -> R(I=1/2)-> Sigma kaon) = 2 * Xsection( pi- D++ -> S0 k+)

parametrizationsBarMes/huangd [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

real function huangd(srts)

PURPOSE

cross section for pi- Delta++ -> Sigma0 K+ in mb

Assumption: all processes happen via I=1/2 resonance

Clebsch for (pi^- D++ -> I=1/2 I_z=1/2) is 1/2

Clebsch for (k^+ S0 -> I=1/2 I_z=1/2) is 1/3

Therefore : Xsection( pi- D++ -> S0 k+)=1/2*1/3* Xsection( pion Delta -> R(I=1/2)-> Sigma kaon)

=> Xsection( pion Delta -> R(I=1/2)-> Sigma kaon) = 6 * Xsection( pi- D++ -> S0 k+)

parametrizationsBarMes/sibirtpi [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

real function sibirtpi(srts)

INPUTS

real :: srts --- sqrt(s)

PURPOSE

cross section for pi N -> N K Kbar

parametrizationsBarMes/golub_omega [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

function golub_omega (srts) result(sigma)

INPUTS

  • real :: srts = sqrt(s) in GeV
  • real :: m_V = mass of vector meson V in GeV, only used for sigma(3) and sigma(4) (unused!)

NOTES

Calculates the cross sections (in mb) for interactions of omegas with nucleons. * sigma(1): pi- p -> V n * sigma(2): pi- p -> V pi N * sigma(3): V N -> V N (unused!) * sigma(4): V N -> inelastic (unused!)

Original parameterizations are taken from:

  • J. Cugnon et al., PRC 41, 1701 (1990) for sigma(pi^- p -> omega n);
  • A. Sibirtsev et al., Z. Phys. A 358, 357 (1997) for sigma(pi^- p -> phi n), sigma(pi^+ p -> omega X) and sigma(pi^+ p -> phi X).
  • Some of these parameterizations are also used in Golubeva NPA625 (1997) 832.
  • See also Effenberger Diss. (A.20)...(A.25)

This routine needs a careful check yet (see subtractions between various cross sections below)!!!!

parametrizationsBarMes/golub_phi [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

function golub_phi (srts, m_V) result(sigma)

INPUTS

  • real :: srts = sqrt(s) in GeV
  • real :: m_V = mass of vector meson V in GeV, only used for sigma(3) and sigma(4)

NOTES

Calculates the cross sections (in mb) for interactions of phis with nucleons. * sigma(1): pi- p -> V n * sigma(2): pi- p -> V pi N * sigma(3): V N -> V N * sigma(4): V N -> inelastic

Original parameterizations are taken from:

  • A. Sibirtsev et al., Z. Phys. A 358, 357 (1997) for sigma(pi^- p -> phi n), sigma(pi^+ p -> omega X) and sigma(pi^+ p -> phi X).
  • Some of these parameterizations are also used in Golubeva NPA625 (1997) 832.
  • See also Effenberger Diss. (A.26)...(A.31)

This routine needs a careful check yet (see subtractions between various cross sections below)!!!!

parametrizationsBarMes/omegaN_lykasov [ Functions ]

[ Top ] [ parametrizationsBarMes ] [ Functions ]

NAME

real function omegaN_lykasov(srts,m,ichannel)

INPUTS

  • real, intent(in) :: srts --- sqrt(s) in GeV
  • real, intent(in) :: m --- mass of the omega meson in GeV
  • integer, intent(in) :: ichannel --- 1=elastic, 2=inelastic

NOTES

calculates omega-nucleon cross section in mb:

  • ichannel=1 : elastic
  • ichannel=2 : inelastic

References:

  • Lykasov,Cassing,Sibirtsev,Rzyanin, Eur. Phys. J. A6 (1999) 71, nucl-th/9811019.
  • see also Muehlich Diss. (9.1),(9.2)

parametrizationsBarMes/JPsiN [ Subroutines ]

[ Top ] [ parametrizationsBarMes ] [ Subroutines ]

NAME

subroutine JPsiN(srts,sigma)

INPUTS

  • real :: srts = sqrt(s) in GeV

NOTES

subroutine for calculation of cross sections for interactions of a J/Psi with a nucleon

OUTPUT

cross sections in mb:

  • sigma(1)=J N -> J N
  • sigma(2)=J N -> Lambda_c Dbar
  • sigma(3)=J N -> Lambda_c D*bar
  • sigma(4)=J N -> N D Dbar

Parameterizations of the boson-exchange model calculations from A. Sibirtsev et al., PRC 63, 044906 (2001)

parametrizationsBarMes/weideta [ Subroutines ]

[ Top ] [ parametrizationsBarMes ] [ Subroutines ]

NAME

subroutine weideta(srts,sig)

INPUTS

  • real :: srts =sqrt(s)

OUTPUT

  • real sig : parameterization of pi- p -> eta n from Weidmann diploma thesis

parametrizationsBarMes/piN_to_strangeBaryon_kaon_pion [ Subroutines ]

[ Top ] [ parametrizationsBarMes ] [ Subroutines ]

NAME

subroutine piN_to_strangeBaryon_kaon_pion(sqrtS,crossSections_piMinus_p,crossSections_piPlus_p)

INPUTS

  • real :: srts =sqrt(s)

RESULT

  • real, dimension(1:4) :: crossSections_piPlus_p
  • real, dimension(1:7) :: crossSections_piMinus_p

crossSections_piPlus_p gives the cross section for pi^+ p -> ...:

  • 1 : Lambda K^+ pi^+
  • 2 : Sigma^0 K^+ pi^+
  • 3 : Sigma^+ K^+ pi^0
  • 4 : Sigma^+ K^0 pi^+

crossSections_piMinus_p gives the cross section for pi^- p -> ...:

  • 1 : Lambda K^0 pi^0
  • 2 : Lambda K^+ pi^-
  • 3 : Sigma^0 K^0 pi^0
  • 4 : Sigma^0 K^+ pi^-
  • 5 : Sigma^- K^+ pi^0
  • 6 : Sigma^- K^0 pi^+
  • 7 : Sigma^+ K^0 pi^-

All cross sections in units of mb!

NOTES

This parametrization is based on a fit by K. Gallmeister (12/2008) * Region in which the fit was performed: p_lab<3.5 GeV. * Data error bars have not been considered in the fit.

parametrizationsBarMes/piN_to_strangeBaryon_kaon_pion_matrix [ Subroutines ]

[ Top ] [ parametrizationsBarMes ] [ Subroutines ]

NAME

subroutine piN_to_strangeBaryon_kaon_pion_matrix(sqrtS,pionCharge,nukCharge,matrix_lambda,matrix_sigma)

INPUTS

  • real :: srts =sqrt(s)
  • integer, intent(in) :: pionCharge -- charge of incoming pion
  • integer, intent(in) :: nucCharge -- charge of incoming nucleon

RESULT

Cross section for kaon pion Lambda production :

  • real, dimension(0:1,-1:1) , intent(out) :: matrix_lambda ! 1st Index: kaon charge, 2nd: pion charge

Cross section for kaon pion Sigma production (Charge of Sigma is fixed by total charge):

  • real, dimension(0:1,-1:1) , intent(out) :: matrix_sigma ! 1st Index: kaon charge, 2nd: pion charge

All cross sections in units of mb!

NOTES

parametrizationsBarMes/sigma_KbarToXi [ Subroutines ]

[ Top ] [ parametrizationsBarMes ] [ Subroutines ]

NAME

real function sigma_KbarToXi(srts,ichannel)

INPUTS

  • real :: srts =sqrt(s)
  • integer :: ichannel = outgoing channel

NOTES

Considered channels:

  • (ichannel = 1) K^- p --> Xi^- K^+
  • (ichannel = 2) K^- p --> Xi^0 K^0
  • (ichannel = 3) K^- p --> Xi^- K^0 \pi^+
  • (ichannel = 4) K^- p --> Xi^- K^+ \pi^0
  • (ichannel = 5) K^- p --> Xi^0 K^+ \pi^-