Arxiv: Production of Strange Particles in Hadronic Interactions

Home , Antiproton

arXiv:1810.09973v1 [hep-ph] 23 Oct 2018 aawt rdcin fvrosMneCromdl PS1. EPOS – models Carlo Monte various of predictions with [ data GeV/c 350 and GeV/c 158 momenta had various initial two of production at the interactions about NA61/ paper the a Recently, published appeared. ration [4] GeV/c 31 momentum initial A1SIEclaoain[] n21,teN6/HN coll on NA61/SHINE data experimental the detailed 2017, sented In we data [1]. experimental collaboration useful NA61/SHINE and interesting very Recently, Introduction 1 in t18GVc[]wr ulse yteclaoain A collaboration. the of by spectra published the were 2016, [3] GeV/c on 158 data at experimental tions fr 2016, range In energy wide [2]. a GeV/c in interactions proton-proton in tion rdcino tag atce nhadronic interactions in particles strange of Production Soﬁa Press, Heron Minkov, N. Gaidarov, M. eds. 37 Vol. THEORY, NUCLEAR uhr fteppr efre oprsno h measured the of comparison a performed papers the of Authors eta AR(S,Germany). (GSI, FAIR at ment Λ lsv rs etosof experime sections reproduces cross well clusive also model FTF Geant4 improved The i rcse ihtopril rdcin – productions t particle one two allows with model processes FTF tic modiﬁed The energies. various at tions tigmse.Det hs hne,w ecietedt ft and of pC, data pp, the in describe production particle we on changes, collaboration these to Due baryon-antibar masses. and string intr mesons strange and impro of have masses probabilities we small of model, cies with FTF strings the quark-gluon In of t mentation problems. (Fritio of FTF these the All of of most version solves others. latest The and problems. GiBUU, vario various with EPOS, have data UrQMD, these calculations: compared model has collaboration The GeV/c. eosi pitrcin t2,3,4,8 n 5 e/,adi and GeV/c, 158 and 80 40, 31, 20, at interactions pp in perons t3 e/.Teclaoainas rsne experimental – presented particles also of collaboration The GeV/c. 31 at Abstract. 2 1 .Galoyan A. indt npouto of production on data sion EN eea Switzerland Geneva, CERN, Russia Dubna, JINR, Σ ¯ and , pp ¯ h A1SIEclaoainhsrcnl ulse ihp high published recently has collaboration NA61/SHINE The → π π 1 ± , ± Σ , .Ribon A. , K Σ ¯ K hc ilb tde nteftr ytePNAexperi- PANDA the by future the in studied be will which , ± ± , , K p Λ ± s 0 , π , , Λ ¯ 2 ± Λ ρ , 0 and and .Uzhinsky V. n rtn npoo-abnitrcin at interactions proton-carbon in protons and , π ω ± K K and , 0 ± K Λ rdcini nirtnpoo interac- antiproton-proton in production K eos rtn,atpoosand antiprotons protons, mesons, ± yeo rdcini - interac- p-p in production hyperon rtn n nirtn produc- antiprotons and protons , ∗ 0 pp ¯ in 1 , → π 2 − Λ C π Λ − ¯ olsosa 5 n 350 and 158 at collisions , interactions. C m2 e/ o158 to GeV/c 20 om pp ¯ epbihdb the by published re o arsceto on creation pair’s yon osi pion-carbon in rons h einn of beginning the t → )mdlo Geant4 of model f) aao production on data dcddependen- oduced HN collabo- SHINE 5]. 9[] UrQMD [6], 99 eNA61/SHINE he iuaerealis- simulate o Cinteractions pC n K sMneCarlo Monte us brto pre- aboration tldt nin- on data ntal experimental + e h frag- the ved K emodels he − , pp Λ ¯ reci- hy- → 1 A. Galoyan, A. Ribon, V. Uzhinsky

3.4 [7] [8], GiBUU 1.6 [9], [10] and Geant4 [11] hadronic models: FTFB and QGSM. Main problems were observed in model descriptions of strange particles production. The aim of our work was to improve the Geant4 FTF model to reproduce inclusive spectra of strange particles. The Geant4 FTF model is based on the Fritiof model [12], [13] of the LUND university. We have modiﬁed the FTF model using the experimentaldata on strange particle and antiproton productions in proton-proton interactions at initial momenta larger than 20 GeV/c [2]. Then, we used the improved FTF model for describing the experimental data on pC and π−C interactions given in the papers [3], [4], [5]. We also used the new FTF model to reproduce differential cross sections of Λ-Λ¯ and K+-K− productions in antiproton-proton interactions [14], [15], [16], and kinematical properties of ¯ 0 Λ, Λ, KS produced in antiptoton-proton interactions in a wide initial energy range – from 360 MeV/c up to 100 GeV/c [17], [18], [19], [20], [21].

2 Inclusive spectra of particles produced in proton-proton interactions at momenta of 20, 31, 40, 80 and 158 GeV/c

It is assumed in the FTF model that hadrons turn into excited states in hadron- hadron inelastic interactions. If only one hadron is excited, the process is called diffraction dissociation. The excited states of hadrons are considered as quark- gluon strings. The Lund fragmentation model is used for quark-gluon string fragmentation. Strange particles are produced by the fragmentation of quark- gluon strings containing a strange quark or antiquark coming from the creation of ss¯ pairs from the vacuum. We have introduced some improvements in the FTF model for a better description of the experimental data [2] on π±, proton, antiproton, and K± productions in pp-interactions at Plab ≥ 20 GeV/c. In Figure 1, the rapidity distributions of K+ and K− mesons produced in pp-interactions at momenta 20, 30, 40, 80 and 158 GeV/c are given. The experimental data are presented by points, the red lines are calculations with the FTF model in the Geant4 10.03.09 development release. As one can see, the calculations overestimate the K meson production at low energy, and describe the data at high energies. The overestimation at low energies is connected with too large probability of strange quark production at string decays. We suppressed the probability of string qq¯ pair production according to the formula:

4 Pss¯ =0.108 · (1 − (mth/Mstr) ), (1)

where 0.108 is an asymptotical value of the probability, mth is a threshold mass equals to the sum of two strange meson masses, Mstr is a string mass. At this probability, the improved FTF model of Geant4 describes well the experimental data at all energies. Only at 40 GeV/c the new FTF model cannot describe K

2 Production of strange particles in hadronic interactions

0.1 0.04 PP-> K+ +X 20 GeV/c PP->K- +X 20 GeV/c 0.08 New FTF New FTF 0.03 10-03-ref-09 FTF 10-03-ref-09 FTF 0.06 exp. data exp. data 0.02 dn/dY

0.04 dn/dY 0.01 0.02

00 0.5 1 1.5 2 2.5 00 0.5 1 1.5 2 2.5 Y Y

0.1 0.04

0.08 PP->K+ +X 30 GeV/c PP->K- +X 30 GeV/c 0.03 0.06 0.02 dn/dY

0.04 dn/dY 0.01 0.02

0 0 0 0.5 1 1.5 2 2.5 0 0.5 1 1.5 2 2.5 Y Y

0.1 0.04

0.08 PP->K+ +X 40 GeV/c PP->K- +X 40 GeV/c 0.03 0.06 0.02 dn/dY

0.04 dn/dY 0.01 0.02

0 0 0 0.5 1 1.5 2 2.5 0 0.5 1 1.5 2 2.5 Y Y

0.1 0.06

0.08 PP->K+ +X 80 GeV/c PP->K- +X 80 GeV/c 0.04 0.06 dn/dY

0.04 dn/dY 0.02 0.02

00 0.5 1 1.5 2 2.5 00 0.5 1 1.5 2 2.5 Y Y

0.1 0.06

0.08 PP->K+ +X 158 GeV/c PP->K- +X 158 GeV/c 0.04 0.06 dn/dY

0.04 dn/dY 0.02 0.02

00 0.5 1 1.5 2 2.5 00 0.5 1 1.5 2 2.5 Y Y

Figure 1. Rapidity distributions of K+ (left side) and K− (right side) mesons in pp interactions at various momenta. meson spectra. We note that these experimental data do not seem to follow the common trend of all the rest of the data. In Figure 2, rapidity distributions of protons and antiprotons produced in pp- interaction at momenta 20, 30, 40, 80 and 158 GeV/c [2] are shown in a comparison with FTF model calculations. Looking at the antiproton spectra calculated by the FTF model of Geant4 10.03.09 release, it is seen that the calculations overestimate antiproton production at low energies. To suppress this over production, we introduced a probability of diquark-antidiquark pair creation at a string decay given by the formula:

2 Pqq−q¯q¯ =0.15 · (1 − (1400 · Nb/Mstr) ), (2) where 0.15 is the asymptotical value of the probability, Mstr is the string mass, and Nb is a maximumnumber of baryonsand antibaryonswhat can be produced. The FTF model calculations with this improvement are shown in the ﬁgure by the blue lines. As seen, the calculations reproduce well the antiproton distributions. It is very important that these changes of the FTF model don’t affect the

3 A. Galoyan, A. Ribon, V. Uzhinsky

proton inclusive spectra. The FTF model describes the proton rapidity distributions quite well. The same situation is observed for π meson spectra.

PP->P + X 20 GeV/c 0.8 New FTF 10-03-ref-09 FTF 0.6 exp. data

0.4 dn/dY

0.2

00 0.5 1 1.5 2 2.5 Y

0.02 0.8 PP->P + X 30 GeV/c PP->Pbar + X 30 GeV/c 0.015 0.6 0.01 0.4 dn/dY dn/dY

0.2 0.005

0 0 0 0.5 1 Y 1.5 2 2.5 0 0.5Y 1 1.5 2

0.02 0.8 PP->P +X 40 GeV/c PP->Pbar +X 40 GeV/c 0.015 0.6

0.01 0.4 dn/dY dn/dY

0.2 0.005

0 0 0 0.5 1 Y 1.5 2 2.5 0 0.5Y 1 1.5 2

0.02 0.8 PP->P + X 80 GeV/c PP->Pbar + X 80 GeV/c 0.015 0.6

0.01 0.4 dn/dY dn/dY

0.2 0.005

00 0.5 1 1.5 2 2.5 00 0.5 1 1.5 2 Y Y

0.02 0.8 PP->P + X 158 GeV/c PP->Pbar + X 158 GeV/c 0.015 0.6

0.01 0.4 dn/dY dn/dY

0.2 0.005

00 0.5 1 1.5 2 2.5 00 0.5 1 1.5 2 Y Y

Figure 2. Rapidity distributions of protons (left side) and antiprotons (right side) in pp interactions at various momenta.

The NA61/SHINE collaboration also published various specrtra of Λ hyperons, especially, rapidity and Feyman X distributions [3] in proton-proton interactions at Plab=158 GeV/c. The authors of the paper compared these data with various model calculations – UrQMD, Fritiof 7.0 and EPOS 1.99. Only EPOS model described the experimental data. Now, the FTF model also describes well these data, as seen in Figure 3. We conclude that the improved FTF model works as well as the EPOS 1.99 model.

3 Production of strange particles in hadron-nucleus interactions

The FTF model uses the Glauber approach for estimating the number of participating hadrons in hadron-nucleus collisions. In the model, the participating hadrons transform into quark-gluon strings which fragment into secondary particles, including the strange particles. The NA61/SHINE collaboration presented measurements of spectra of charged 0 pions, kaons, and protons as well as Ks and Λ hyperons in pC interactions at 31

4 Production of strange particles in hadronic interactions

exp. dat exp. dat New FTF New FTF 150 10-03-ref-09 FTF 10-03-ref-09 FTF 40

100

F dn/dY 20 dn/dX 50

0 0 −2 −1 0 1 2 −1 −0.5 0 0.5 1 Y XF

Figure 3. Rapidity (left side) and Feyman X (right side) distributions of Λ-hyperons in pp interactions at 158 GeV/c.

GeV/c in the paper[4]. In that paper,data are comparedwith calculations of FTF + Binary Cascade model, Venus, EPOS1.99, and GiBUU models. The conclusion of the paper was that none of the models provided a satisfactory description of all the spectra.

20 20 PC->K++X 0 <θ <20 mrad 18 18 PC->K++X 20 <θ< 40 mrad 16 New FTF 16 14 14 12 10-03-ref-09 FTFP 12 10 exp.data 31 GeV/c 10 8 8 6 6 /dp (mb/GeV/c) /dp (mb/GeV/c)

σ 4 σ 4 d 2 d 2 00 5 10 15 20 00 5 10 15 20 p (GeV/c) P (GeV/c)

40 40 + + 35 PC->K +X 40 <θ< 60 mrad 35 PC->K +X 60 <θ< 100 mrad 30 30

25 25

20 20

15 15

/dp (mb/GeV/c) 10 /dp (mb/GeV/c) 10 σ σ

d 5 d 5

0 0 0 5 10 15 20 0 5 10 15 20 p (GeV/c) P (GeV/c)

40 40 35 PC->K++X 100 <θ< 140 mrad 35 PC->K++X 140 <θ< 180 mrad 30 30 25 25 20 20 15 15

/dp (mb/GeV/c) 10 /dp (mb/GeV/c) 10 σ σ

d 5 d 5

00 5 10 15 20 00 5 10 15 p (GeV/c) p (GeV/c)

40 40 35 PC->K++X 180 <θ< 240 mrad 35 PC->K++X 240 <θ< 300 mrad 30 30 25 25 20 20 15 15

/dp (mb/GeV/c) 10 /dp (mb/GeV/c) 10 σ σ

d 5 d 5

00 5 10 15 00 2 4 6 8 10 P (GeV/c) P (GeV/c)

Figure 4. Momentum distributions of K+ mesons in pC interactions at 31 GeV/c at various angles.

In Figure 4, we present the differential cross sections of K+ mesons produced in pC interactions as functions of laboratory momentum of the mesons. The experimental data are given at various emission angles of K+ – 0-20 mrad, 20-40 mrad, 40-60 mrad and so on up to 240-300 mrad. Before the improve-

5 A. Galoyan, A. Ribon, V. Uzhinsky

ment, the FTF model calculations (red lines) overestimated the maxima of the distributions by about 50%. After the improvement, the FTF model reproduces the inclusive spectra of K+ mesons quite well. The same situation takes place − 0 for the momentum distributions of K and KS mesons produced in pC interactions at 31 GeV/c. Now, we have reached a good agreement between the FTF calculations and the experimental inclusive spectra of mesons and Λ hyperons. Very interesting experimental data on π−C interactions have been published by the NA61/SHINE collaboration in Ref. [5].

1.4 π- + C -> π+ + X 1.4 π π 1.2 New FTF 1.2 - + C -> - + X 10-03-ref-09 FTFP 1 exp.dat 158 GeV/c 1 0.8 0.8

0.6 0.6 p dn/dp p dn/dp 0.4 0.4

0.2 0.2

0 0 1 P (GeV/c) 10 1 P (GeV/c) 10

0.15 0.15 π- + C -> K+ + X π- + C -> K- + X

0.1 0.1

p dn/dp 0.05 p dn/dp 0.05

0 0 1 P (GeV/c)10 1 P (GeV/c)10

0.2 0.06 π- + C -> P + X π- + C -> Pbar + X 0.15 0.04

0.1

p dn/dp p dn/dp 0.02 0.05

0 0 1 P (GeV/c) 10 1 P (GeV/c)10

Figure 5. Momentum distributions of π+, π−, K+, K−, p, p¯ in π−C interactions at 158 GeV/c.

Figure 5 shows the experimental data on π± , K± , proton and antiproton momentum distributions in π−C interactions at 158 GeV/c. The authors of the paper [5] showed also calculations of different models: QGSJet, Sybill, DpmJet and EPOS. The main difference between the models was observed in the calculations of K+, K− mesons and antiprotons. In Figure 5, we present the FTF model calculations. In general, the improved FTF model calculations (blue lines) describe better the experimental data than other models presented in Ref. [5]. The FTF model results are in agreement with the data on π±, protons, and K± mesons. At the same time, the FTF calculations are slightly lower than the data for antiprotons with momenta larger than 10 GeV/c. We have a similar situation for the description of the experimental data on pion, kaon, proton and antiproton production in π−C interactions at 350 GeV/c [5]. Summing up, we conclude that the FTF model with the new probabilities describes quite well the inclusive spectra of protons, antiprotons, π±, K±, K0 mesons and Λ-hyperons produced in pp, pC, and π−C interactions.

6 Production of strange particles in hadronic interactions

4 Production of strange particles in antiproton-proton interactions

A detailed experimental study of antiproton-protonand antiproton-nucleus interactions is foreseen at the future FAIR facilities (GSI, Darmstadt, Germany) by the PANDA Collaboration [22]. Simulations of antiproton interactions with protons and nuclei are needed for cosmic space experiments like PAMELA, BESS, AMS, CAPRISE, which are going to search for anti-nuclei in cosmic rays. For the simulation of antiproton-proton interactions, the Geant4 FTF model uses the main assumptions of the Quark-Gluon-String Model [23]. The FTF model assumes production and fragmentation of quark-anti-quark and diquark-antidiquark strings in the mentioned interactions. Main ingredients of the FTF model are cross sections of string creation processes [24] and the use of the LUND string fragmentation algorithm. They allow to achieve a satisfactory description of a large set of experimental data, in particular, strange particle production, such as Λ hyperons and K mesons [25]. To develop the hypernuclei program of the Panda experiment, it is needed to estimate multiplicity of hyperons and correctly reproduce their kinematical characteristics. Thus, we calculate 0 with the FTF model the average multiplicities of Λ hyperons and KS mesons in anti-proton-proton interactions presented in Figure 6.

102 102 Pbar+P -> Λ + X Pbar+P -> Ks0 + X New FTF New FTF 10-03-ref-09 FTF 10-03-ref-09 FTF exp. data exp.data 10 10

1 1 (mb) (mb) σ σ

− − 10 1 10 1

− − 10 2 10 2 2 1 10 10 1 10 102 Plab (GeV/c) Plab (GeV/c)

0 Figure 6. Inclusive cross sections of Λ-hyperon and Ks -meson production in antiproton- proton interactions as function of projectile momentum.

The Figure shows experimental data [17] on inclusive cross section of Λ hy- 0 peron and KS meson in antiproton-protoninteractions as a function of projectile 0 momentum. As seen, the FTF model describes well KS meson cross sections at 0 low energies, and overestimates, a little bit, the yield of KS mesons at momenta more than 10 GeV/c. Calculations by the old version of FTF underestimated the production cross section of Λ hyperons at low momenta, less than 2 GeV/c, and were in agreement with the experimental cross section at higher energies up to 30 GeV/c. The developed FTF model calculations are in agreement with Λ production cross section at low energies up to 3 GeV.

7 A. Galoyan, A. Ribon, V. Uzhinsky

Distributions of the strange particles on kinematical variables are very im- 0 portant for experiments. Distributions of Λ hyperons and KS mesons on rapidi- ties, Feyman X, and squared transverse mometum in pp¯ interactions at momenta of 12, 24 and 100 GeV/c calculated with the FTF model are presented in Figure 7.

800 104 Pbar+P->Ks0 12 GeV/c

New FTF 4 Pbar+P->Λ 12 GeV/c 10 Pbar+P -> Ks0 12 GeV/c Pbar+P->Λ 12 GeV/c 400

10-03-ref-09 FTF 600 exp. data 2 2 3 10

3 10 /dY(b)

400 σ /dY(b) mb/(GeV/c) mb/(GeV/c) σ 200 2 2 d 2

d 10 10/dPt 2 /dPt 200 σ σ d d

0 1 2 3 4 0 1 2 3 0 0.5 1 0 0.5 1 2 2 2 2 Y Y Pt (GeV/c) Pt (GeV/c)

0.8 0.5

10 Pbar+P->Ks0 22.4 GeV/c Pbar+P->Λ 22.4 GeV/c 10 Pbar+P->Ks0 22.4 GeV/c Pbar+P->Λ 22.4 GeV/c 0.4

0.6 2 2

0.3 1 1 ) ) f f 0.4 F(X F(X mb/(GeV/c) mb/(GeV/c) 0.2 2 2

− 1

/dPt − /dPt 10 1 10 0.2 σ σ 0.1 d d

− 2 − 2 − − − − 10 10 1 0.5 0 0.5 1 1 0.5 0 0.5 1 0 0.52 2 1 0 0.52 2 1 Xf Xf Pt (GeV/c) Pt (GeV/c)

3 2

Pbar+P->Ks0 100 GeV/c Pbar+P->Λ 100 GeV/c Pbar+P->Ks0 100 GeV/c Pbar+P->Λ 100 GeV/c 2.5 10

1.5 2 10 2

1.5 1 mb/(GeV/c) mb/(GeV/c) /dY(mb) /dY(mb) 1 2 2

σ σ 1

d 1 d /dPt /dPt

0.5 σ σ d d 0.5 − 10 1

− 10 1 −2 0 2 −2 0 2 0 0.5 1 0 0.5 1 2 2 2 2 Y Y Pt (GeV/c) Pt (GeV/c)

2 0 Figure 7. Rapidity, Feyman X, Pt distributions of Ks -mesons and Λ-hyperon in pp¯ interactions at 12, 22.4, 100 GeV/c.

As seen, there is a good agreement between the improved FTF model calculations (blue lines) and the experimental data on the inclusive spectra at 12 GeV/c [19] and 22.4 GeV/c [20]. When we consider antiproton-protoninteractions at very high momenta, 100 GeV/c, annihilation processes are absent. Here, only diffraction processes of 0 FTF are working. At this energy, FTF gives a quite good description of KS and 2 Λ rapidity and PT spectra. Let us consider now binary reactions with strange particle production in antiproton-proton annihilations. There are a lot of experimental data on differential cross sections of Λ-Λ¯ production in antiproton-proton interactions. We have gathered most of them, performed corresponding simulations by the FTF model and compared with the experimental data. For example, in Figure 8 we present differential cross sections of Λ-Λ¯ at various initial momenta: 1.56, 1.76, 1.88 and 2.06 GeV/c.

8 Production of strange particles in hadronic interactions

Pbar+P -> Λ + ΛBar Pbar+P -> Λ + ΛBar New FTF exp.dat 1.76 GeV/c 10-03-ref-09 FTF

2 102 exp.dat 1.56 GeV/c 2 102

10 10 /dT mb/(GeV/c) /dT mb/(GeV/c) σ σ d d

1 1 0 0.2 0.4 0.6 0.8 0 0.2 0.4 0.62 0.8 1 T (GeV/c)2 T (GeV/c)

3 10 Pbar+P -> Λ + ΛBar Pbar+P -> Λ + ΛBar exp.dat 1.88 GeV/c exp.dat 2.06 GeV/c

2 102 2 102

10 10 /dT mb/(GeV/c) /dT mb/(GeV/c) σ σ d d

1 1 0 0.5 1 1.5 0 0.5 1 1.5 T (GeV/c)2 T (GeV/c)2

Figure 8. Differential cross section of Λ-hyperon production as function of transferred momentum in pp¯ interactions at various energies.

To describe the data [15], we introduced in the FTF model the spatial rota- tion of quark-gluon strings. The FTF model without rotating strings does not describe the experimental data [15]. The FTF model results with the rotating strings are in agreement with the experimental distributions.

5 Conclusion

A new formula for the probability of strange qq¯ pair production at string decay is proposed and implemented in the FTF model of Geant4. The new probability improves considerably the description of the strange particle production in the FTF model. A new formula for the probability of diquark-antidiquark pair production at string decay is proposed and implemented in the FTF model. The new probability helps to improve the description of antiproton and Λ hyperon production in the FTF model. Good agreementbetween the improvedFTF model calculations and the NA61/ ± 0 SHINE experimental data is reached on K , KS , and Λ production in proton- proton, proton-Carbon and π−-meson-Carbon collisions in a wide energy range. 0 Kinematical properties of Λ hyperons and KS mesons produced in pp¯ reactions are calculated in the improved FTF model with the rotating strings, and compared with the experimental data at various initial momenta. Good description of the experimental data is obtained in the FTF model with the new probabilities.

9 A. Galoyan, A. Ribon, V. Uzhinsky

The differential cross sections of the processes pp¯ → ΛΛ¯ are calculated in various versions of the FTF model. Reasonable description of the experimental data is reached in the FTF model with the rotating strings and the new probabilities.

6 Acknowledgements

A. Galoyan is thankful to heterogeneous computing team of LIT JINR (Hy- briLIT) for support of calculations and to Professor Cheuk-Yin Wong for helpful discussions.

References

[1] URL:https://home.cern/about/experiments/na61shine. [2] NA61/SHINE Collaboration, Eur.Phys.J. C77 (2017) no.10, 671. [3] NA61/SHINE Collaboration Eur. Phys. J. C76 (2016) no.4, 198. [4] NA61/SHINE Collaboration Eur. Phys. J. C76 (2016) no.2, 84 [5] NA61/SHINE Collaboration (Raul R. Prado (Sao Paulo U., Sao Carlos) Jul 25, 2017. Conference: C17-07-12 Proceedings. arXiv:1707.07902 [hep-ex] [6] K. Werner, Nucl. Phys. Proc. Suppl. 175-176 (2008) 81-87. [7] S. Bass et al., Prog.Part.Nucl.Phys. 41 (1998) 255, arXiv:9803035 [nucl-th]. [8] M. Bleicher et al., J.Phys. G25 (1999) 1859, arXiv:9909407 [hep-ph]. [9] O. Buss, T. Gaitanos, K. Gallmeister, H. van Hees, M. Kaskulov, et al., Phys. Rept. 512 (2012) 1. [10] K. Gallmeister, U. Mosel, Nucl. Phys. A826 (2009) 151. [11] J.Allison et al., Nucl.Instrum.Meth. A835 (2016) 186. [12] B. Andersson et al., Nucl. Phys. B281 (1987) 289. [13] B. Nilsson-Almquist and E. Stenlund, Comp. Phys. Commun. 43 (1987) 387. [14] P.D. Barnes et al., Nuclear Physics A558 (1993) 277. [15] B. Jayet et al., Nuov. Cim. A45 (1978) 371. [16] H. Becker et al., Nucl. Phys. B141 (1978) 48. [17] S. Banerjee et al., TIFR-BC-78-8 [18] A.M. Cooper et al., Nucl.Phys. B136 (1978) 365. [19] D. Bertrand et al., Nucl. Phys. B128 (1977) 365. [20] B.V. Batyunya et al., Z. Phys. C25 (1984) 213. [21] D.R. Ward et al., Phys. Lett. 62B (1976) 237. [22] PANDA Collaboration (M.F.M. Lutz et al.) arXiv:0903.3905 [hep-ex] (2009). [23] A. Capella, U. Sukhatme, C-I Tan, J. Tran Thanh Van, Phys. Rept. 236 (1994) 225; A.B. Kaidalov and K.A. Ter-Martirosian, Phys. Lett. B 117 (1982) 247. [24] A. Galoyan, A. Ribon and V. Uzhinsky PoS BaldinISHEPPXXII (2015) 049. [25] A. Galoyan, A. Ribon, V. Uzhinsky arXiv:1610.08341 [nucl-th] (2016).