TABLE OF CONTENTS
- 1. /parametrizationBarMes
- 1.1. parametrizationBarMes/piN_elastic
- 1.2. parametrizationBarMes/piN_chargeExchange
- 1.3. parametrizationBarMes/pionquasibg
- 1.4. parametrizationBarMes/sigelkp
- 1.5. parametrizationBarMes/sigCEXkp
- 1.6. parametrizationBarMes/sigCEXk0
- 1.7. parametrizationBarMes/kppbg
- 1.8. parametrizationBarMes/kpnbg
- 1.9. parametrizationBarMes/kaonbg
- 1.10. parametrizationBarMes/kaonbgt
- 1.11. parametrizationBarMes/huang
- 1.12. parametrizationBarMes/huanglam
- 1.13. parametrizationBarMes/huanglamd
- 1.14. parametrizationBarMes/huangd
- 1.15. parametrizationBarMes/sibirtpi
- 1.16. parametrizationBarMes/golub_omega
- 1.17. parametrizationBarMes/golub_phi
- 1.18. parametrizationBarMes/omegaN_lykasov
- 1.19. parametrizationBarMes/JPsiN
- 1.20. parametrizationBarMes/weideta
- 1.21. parametrizationBarMes/piN_to_strangeBaryon_kaon_pion
- 1.22. parametrizationBarMes/piN_to_strangeBaryon_kaon_pion_matrix
- 1.23. parametrizationBarMes/sigma_KbarToXi
/parametrizationBarMes [ Modules ]
NAME
module parametrizationBarMes
PURPOSE
Includes routines which are parametrizations of "baryon meson -> X" data.
parametrizationBarMes/piN_elastic [ Functions ]
[ Top ] [ parametrizationBarMes ] [ 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
parametrizationBarMes/piN_chargeExchange [ Functions ]
[ Top ] [ parametrizationBarMes ] [ Functions ]
NAME
real function piN_chargeExchange(srts)
PURPOSE
Returns the cross section for pion nucleon charge exchange, i.e.:
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
parametrizationBarMes/pionquasibg [ Subroutines ]
[ Top ] [ parametrizationBarMes ] [ 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
parametrizationBarMes/sigelkp [ Functions ]
[ Top ] [ parametrizationBarMes ] [ 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)
parametrizationBarMes/sigCEXkp [ Functions ]
[ Top ] [ parametrizationBarMes ] [ Functions ]
NAME
real function sigCEXkp(srts)
PURPOSE
Calculate * K+ n -> K0 p
INPUTS
- real :: srts -- invariant energy (GeV)
parametrizationBarMes/sigCEXk0 [ Functions ]
[ Top ] [ parametrizationBarMes ] [ Functions ]
NAME
real function sigCEXk0(srts)
PURPOSE
Calculate * K0 p -> K+ n
INPUTS
- real :: srts -- invariant energy (GeV)
parametrizationBarMes/kppbg [ Functions ]
[ Top ] [ parametrizationBarMes ] [ Functions ]
NAME
real function kppbg(plab2)
PURPOSE
Calculate * K+ p -> p K pion
INPUTS
- real :: plab2
parametrizationBarMes/kpnbg [ Functions ]
[ Top ] [ parametrizationBarMes ] [ Functions ]
NAME
real function kpnbg(plab2)
PURPOSE
Calculate:
- K0 p -> n K pion
INPUTS
- real :: plab2
parametrizationBarMes/kaonbg [ Functions ]
[ Top ] [ parametrizationBarMes ] [ 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
parametrizationBarMes/kaonbgt [ Functions ]
[ Top ] [ parametrizationBarMes ] [ Functions ]
NAME
real function kaonbgt(plab2)
PURPOSE
Calculate K- p total cross section im mb
parametrizationBarMes/huang [ Functions ]
[ Top ] [ parametrizationBarMes ] [ 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
parametrizationBarMes/huanglam [ Functions ]
[ Top ] [ parametrizationBarMes ] [ Functions ]
NAME
real function huanglam(srts)
PURPOSE
cross section for pi- p -> Lambda K0 in mb
parametrizationBarMes/huanglamd [ Functions ]
[ Top ] [ parametrizationBarMes ] [ 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+)
parametrizationBarMes/huangd [ Functions ]
[ Top ] [ parametrizationBarMes ] [ 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+)
parametrizationBarMes/sibirtpi [ Functions ]
[ Top ] [ parametrizationBarMes ] [ Functions ]
NAME
real function sibirtpi(srts)
INPUTS
real :: srts --- sqrt(s)
PURPOSE
cross section for pi N -> N K Kbar
parametrizationBarMes/golub_omega [ Functions ]
[ Top ] [ parametrizationBarMes ] [ 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)!!!!
parametrizationBarMes/golub_phi [ Functions ]
[ Top ] [ parametrizationBarMes ] [ 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)!!!!
parametrizationBarMes/omegaN_lykasov [ Functions ]
[ Top ] [ parametrizationBarMes ] [ 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)
parametrizationBarMes/JPsiN [ Subroutines ]
[ Top ] [ parametrizationBarMes ] [ 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)
parametrizationBarMes/weideta [ Subroutines ]
[ Top ] [ parametrizationBarMes ] [ Subroutines ]
NAME
subroutine weideta(srts,sig)
INPUTS
- real :: srts =sqrt(s)
OUTPUT
- real sig : parameterization of pi- p -> eta n from Weidmann diploma thesis
parametrizationBarMes/piN_to_strangeBaryon_kaon_pion [ Subroutines ]
[ Top ] [ parametrizationBarMes ] [ 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.
parametrizationBarMes/piN_to_strangeBaryon_kaon_pion_matrix [ Subroutines ]
[ Top ] [ parametrizationBarMes ] [ 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 :
Cross section for kaon pion Sigma production (Charge of Sigma is fixed by total charge):
- real, dimension(0:1,-1:1) :: matrix_sigma
-- 1st Index: kaon charge, 2nd: pion charge
All cross sections in units of mb!
NOTES
- This parametrization is based on the subroutine piN_to_strangeBaryon_kaon_pion
- No Xsections for pi^0 induced reactions!!
parametrizationBarMes/sigma_KbarToXi [ Subroutines ]
[ Top ] [ parametrizationBarMes ] [ 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^-