Central exclusive diffractive production of two-pion continuum at hadron colliders

Calculations of central exclusive diffractive di-pion continuum production are presented in the Regge-eikonal approach. Data from ISR, STAR, CDF and CMS were analyzed and compared with theoretical description. We also consider theoretical predictions for LHC, possible nuances and problems of calculations and prospects of investigations at present and future hadron colliders.


Introduction
In previous papers [1], [2] the general properties and calculations of the Central Exclusive Diffractive Production (CEDP) were considered. It was shown, especially in [2], that diffractive patterns (differential crosssections) of CEDP play a significant role in model verification.
Here we partially continue the subject of [2] and investigate in detail the process of low mass CEDP (LM CEDP) with production of two pions. This process is one of the "standard candles" for LM CEDP. Why do we need exact calculations and predictions for this process?
-Di-pion LM CEDP is the basic background process for CEDP of resonances (like f 2 or f 0 ), since one of the basic hadronic decay modes for these resonances is the two pion one. -We can use LM CEDP to fix the procedure of calculations of "rescattering" (unitarity) corrections. In a e-mail: Roman.Rioutine@cern.ch the case of di-pion LM CEDP there are two kinds of corrections, in the proton-proton and the pionproton subamplitudes. They will be considered in the present work. -The pion is the most fundamental particle in the strong interactions, and LM CEDP gives us a powerful tool to thoroughly investigate its properties, especially to investigate the form factor and scattering amplitudes for the off-shell pion. -LM CEDP has rather large cross-sections. It is very important for an exclusive process, since in the special low luminosity runs (of the LHC) we need more time to get enough statistics. -As was proposed in [3], it is possible to extract some Reggeon-hadron cross-sections. In the case of single and double dissociation it was the Pomeronproton one. Here, in the LM CEDP of the di-pion we can analyze properties of the Pomeron-Pomeron to pion-pion exclusive cross-section, and also check again the predictions of the covariant reggeization method [3]. -Diffractive patterns of this process are very sensitive to different approaches (subamplitudes, form factors, unitarization, reggeization procedure), especially differential cross-sections in t and φ pp (azimuthal angle between final protons), and also the M ππ dependence. That is why this process is used to verify different models of diffraction. -All the above items are additional advantages provided by the LM CEDP of two pions, which has the usual properties of CEDP: a clear signature with two final protons and two large rapidity gaps (LRG) [4], [5] and the possibility to use the missing mass method [6].
Processes of the LM CEDP were calculated in some other work [7]- [13] devoted to the most popular models for the LM CEDP of di-mesons. All authors have considered nonperturbative approach in reggeon-reggeon collision subprocess. For example, in the Durham model [7]- [9] (see Fig. 1), they take the Born Regge term for the amplitude of the process p+p → p+π + +π − +p with reggeized propagator of the off-shell pion, then they take into account unitarity (rescattering) corrections (in the initial proton-proton state and also in the so-called "enhanced" one). The authors of [7]- [9] point out, however, that "enhanced" corrections are negligible due to the small triple Pomeron vertex. In addition, the possibility to include reggeization of the virtual pion propagator is not obvious, since the effect of this is expected to be small and, moreover, it is not even clear that we are in the relevant kinematic region (|t| ŝ = M 2 ππ ) to include such corrections for central production. This is therefore not included by default into their calculations (see [9], for example). For example, the authors of [7]- [9] use the replacement which gives the correct "reggeized" behavior in the relevant kinematical region, and the usual "bare" pion propagator behavior for a small difference between rapidities of pions. The authors of [10]- [13] use a phenomenological expression for the virtual pion propagator (see (3) for notations, and also (3.25) and (3.26) of [12]) like to take into account possible non-Regge behavior for t ∼ŝ/2, i.e. for small rapidity separation ∆Y between the final pions. In this paper we do not use such a "mixed" method, and consider these two cases ("bare" or "reggeized" virtual pion propagator) separately, especially to see the relevance of the Regge approach in the kinematical regiont ∼ŝ/2 ∼ 1 GeV 2 . The Regge model really does not work in this area or it needs to be modified (as was done, for example, in Refs. [10]- [13], with empirical formulas or additional assumptions), and also taking into account, for example, possible significant contribution of the background integral in the Sommerfeld-Watson transform (see [14]) or taking Legendre polynomials of order α π (t) instead of the classical Regge term (see (3)) in the "reggeized" virtual pion propagator. All of this should be verified in future research.
In Refs. [10]- [13] the authors do not introduce "enhanced" corrections, but they take into account pionproton interactions in the final state (see Fig. 2). There still is an issue in this approach, since they mix a partially Regge approach with continuous complex spin and a model with fixed "Pomeron spin" (exactly vector or tensor Pomerons). It may be convenient for calculations and gives results close to reality, but we have no clear physical explanation. On the one hand we have collision of two particles with fixed spin (1 or 2), and on the other hand we use the Regge expression, where the spin is replaced by the complex Regge trajectory.
In the first part of the present work we introduce the framework for calculations of double pion LM CEDP (kinematics, amplitudes, differential cross-sections) in the Regge-eikonal approach.
In the second part we analyze the experimental data on the process at different energies, find the best approach and make some predictions for LHC experiments.
In the final part we discuss the possibilities to extract Pomeron-Pomeron cross-sections from the data and analyze the present situation. Also we show some nuances of the calculations, which we should take into account (elastic amplitudes for virtual particles, offshell pion form factor, pion-pion elastic amplitude at low energies, and nonlinearity of the pion trajectory). Fig. 1 Amplitude of the process of double pion LM CEDP p+p → p+π + +π − +p in the KMR approach [7]- [9]. The central part of the diagram is the Born amplitude (with Pomeron and two reggeons). A reggeized off-shell pion propagator is shown as a dashed zigzag line. Proton-proton rescattering is depicted as the S pp -blob, and "enhanced" corrections are also shown as the S enh -blob. The off-shell pion form factor is presented as a black circle. Fig. 2 Amplitude of the process of double pion LM CEDP p + p → p + π + + π − + p in the approach [10]- [13]. The central part of the diagram is the Born amplitude (with Pomeron and two reggeons). A mixture of reggeized and bare off-shell pion propagators is shown as a dashed zigzag line plus dashed straight line. Proton-proton rescattering is depicted as the S pp -blob, and pion-proton rescattering corrections are also shown as S πp -blobs. The off-shell pion form factor is presented as a black circle. Fig. 3 Amplitude of the process of double pion LM CEDP p + p → p + π + + π − + p in the Regge-eikonal approach for different cases. PB: central part of the diagram is the Born amplitude (with Pomeron and two reggeons) plus an off-shell bare pion propagator depicted as a dashed straight line. RB: the same as PB, but with the reggeized off-shell pion propagator depicted as a dahshed zigzag line. PF: central part of the diagram contains full eikonalized pion-proton amplitudes plus off-shell bare pion propagator depicted as a dashed straight line. RF: the same as before but with reggeized off-shell pion propagator depicted as a dashed zigzag line. Proton-proton rescatterings in the initial and final states are depicted as V pp -and V pp -blobs, respectively, and pion-proton rescattering corrections are also shown as S πp -blobs. The offshell pion form factor is presented as a black circle.

General framework for calculations of LM CEDP
LM CEDP is the first exclusive two to four process which is driven by the Pomeron-Pomeron fusion subprocess. That is why it serves as a basic background for LM CEDP of resonances like f 0 (980), f 2 (1270). At the moment, for low central pion-pion masses (less than ∼ 3 GeV), it is a huge problem to use a perturbative approach, that is why we apply the Regge-eikonal method for all the calculations. For proton-proton and protonpion elastic amplitudes we use the model of [23], [24], which describes all the available experimental data on elastic scattering.

Components of the framework
LM CEDP process can be calculated in the following scheme (see Fig. 3): 1. We calculate the primary amplitude of the process, which is depicted as the central part of diagrams in Fig. 3. Here we consider four cases to show that only one of them gives the best description of the data on this process. The case PB represents the Born term (the letter P means that we use the usual "bare" virtual pion propagator 1/(t − m 2 π ), and B means the Born term (28) in each shoulder or pionproton elastic subamplitudes of the central primary amplitude). The expressions of pion-proton elastic subamplitudes can be found in Appendix B. The case RB is similar to the previous one, and the letter R (instead of P) means that the bare off-shell pion propagator is replaced by the reggeized one whereŝ is the di-pion mass squared andt is the square of the momentum transfer between a Pomeron and a pion in the Pomeron-Pomeron fusion process (see Appendix A for details). The cases PF (RF) can be obtained from PB (RB) if we replace the Born pion-proton elastic amplitudes to full eikonalized expressions (which is reflected in the replacement of B to F in the notation for these cases), which could be found in Appendix B. 2. After the calculation of the primary LM CEDP amplitude we have to take into account all possible corrections in proton-proton and proton-pion elastic channels due to the unitarization procedure (the so called "soft survival probability" or "rescattering corrections"), which are depicted as V pp , V pp and S πp blobs in Fig. 3. For the proton-proton and proton-pion elastic amplitudes we use the model of [23], [24] (see Appendix B). A possible final pionpion interaction is not shown in Fig. 3, since we neglect it in the present calculations. RB and PB cases of Fig. 3 are similar to the one in Fig. 2.
In this article we do not consider the so called "enhanced" corrections [7]- [9], since they give nonleading contributions in our model due to the smallness of the triple Pomeron vertex. Also we have no possible absorptive corrections in the pion-pion final elastic channel, since the central mass is low, and also there is a lack of data on this process to define parameters of the model. Nevertheless we will consider these corrections in further work, as was done by some authors recently [25], since they could play significant role for masses less than 1 GeV.
Exact kinematics of the two to four process is outlined in Appendix A.
Here we use the model, presented in Appendix B for example. One can use another one, which has been proved to describe well all the available data on protonproton and proton-pion elastic processes. But it is difficult to find now more than a couple of models which have more or less predictable power (see [26] for a detailed discussion). That is why we use the model, which has been proved to be good in data fitting, especially in the kinematical region of our interest.
The final expression for the amplitude with protonproton and pion-proton "rescattering" corrections can be written as where the functions are defined in (30)-(34) of Appendix B, and the sets of vectors are {p} ≡ {p a − q, p b + q; p 1 + q + q 1 , and The off-shell pion form factor is equal to unity on the mass shell,t = m 2 π , and taken as exponential where Λ π is taken from the fits to LM CEDP of two pions at low energies (see next section). In this paper we use only exponential form, but it is possible to use other parametrizations (see [7]- [13]). The xponential form shows more appropriate results in the data fitting.
Other functions are defined in Appendix B. Then we can use Eq. (23) to calculate the differential crosssection of the process.

Nuances of calculations
In the next section one can see that there are some difficulties in the data fitting, which have also been presented in other work [10]- [13]. In this section let us discuss some nuances of calculations, which could change the situation.
We have to pay special attention to the amplitudes, where one or more external particles are off their mass shell. The example of such an amplitude is the pionproton one T π + p (T π − p ), which is the part of the CEDP amplitude (see (4)). For this amplitude in the present paper we use the Regge-eikonal model with the eikonal function in the classical Regge form. The "off-shell" condition for one of the pions is taken into account by an additional phenomenological form factorF π (t). But there are at least two other possibilities.
The first one was considered in [27]. For the amplitude with one particle off-shell the formula was used. In our case δ πp is the eikonal function (see (26)). This is similar to the introduction of the additional form factor, but in a more consistent way, which takes into account the unitarity condition. The second one arises from the covariant reggeization method, which is considered in Appendix C. For the case of conserved hadronic currents we have the definite structure in the Legendre function (50), which is transformed in a natural way to the case of the offshell amplitude. But in this case the off-shell amplitude shows a specific behavior at low t values (see Fig. 4 and [2] for details). As was shown in [2], unitarity corrections can mask this behavior. To check this we need to make all the calculations and fitting of the data for the process p + p → p + π + π + p, but with an amplitude like (50) instead of (33). This will be done in further work on the subject.
In the present calculations we use the linear pion trajectory 0.7(t − m 2 π ). The nonlinear case was also verified, and the difference in the final result is not significant.

Data from hadron colliders versus results of calculations
Our basic task is to extract the fundamental information on the interaction of hadrons from different crosssections ("diffractive patterns"): from t-distributions we can obtain the size and shape of the interaction region; the distribution on the azimuthal angle between the final protons gives the quantum numbers of the produced system (see [2],[28] and the references therein); from M c (here M c = M ππ ) dependence and its influence on the t-dependence we can draw some conclusions about the interaction at different space-time scales and the interrelation between them.
The process p + p → p + π + π + p is the first "standard candle", which we can use to estimate other LM CEDP processes, like resonance production [28], [29]. In this section we consider the experimental data on the process and its description for different model cases.

STAR collaboration data versus model cases
In this section the data of the STAR collaboration [17], [18] and model curves for different cases of Fig. 3 are presented. In our approach we have only one free parameter, Λ π , that is why all the distributions are depicted for its different values. Also in every case we consider two possibilities: fitting data by formulas with all rescattering corrections (two upper pictures) and also in the approach with protonproton rescattering only (i.e. fitting the data by formulas without pion-proton interactions in the final state, two lower pictures). We change Λ π and try to get the best description. As you can see from Figs. 5-8, the best description is given in the RF case for both possibilities (see Fig. 5). Since the final pion-proton interaction can give rather large suppression (about 10-20%, as in Fig. 14), in our further calculations we use the full amplitude as depicted in Fig. 3 for the RF case. The RF case without pion-proton interactions in the final state (with its own values of Λ π for the best data description) we will show just for cheking of this possibility.  Fig. 3). The data on the process p + p → p + π + + π − + p at √ s = 200 GeV, |η π | < 1, |η ππ | < 2, p T π > 0.15 GeV, 0.005 < −t 1,2 < 0.03 GeV 2 , (STAR collaboration [17], [18]). a, b The result with all p p and π p rescattering corrections, c, d show the result, when we try to fit the data by formulas without π p rescattering terms. Curves from up to down correspond to different values of the parameter Λ π in the off-shell pion form factor (11) Fig. 3). The data on the process p + p → p + π + + π − + p at √ s = 200 GeV, |η π | < 1, |η ππ | < 2, p T π > 0.15 GeV, 0.005 < −t 1,2 < 0.03 GeV 2 , (STAR collaboration [17], [18]). a, b The result with all p p and π p rescattering corrections, c, d show the result when we try to fit the data by formulas without π p rescattering terms. Curves from up to down correspond to different values of the parameter Λ π in the off-shell pion form factor (11): a, b Λ π = 1.6, 1.2, 0.8 GeV, c, d Λ π = 1, 0.8, 0.6 GeV.

ISR and CDF data versus RF case of the model
Let us look at the ISR [15], [16] and CDF [19], [20] data with parameter Λ π , which we use to describe the data from STAR collaboration. Different cases are depicted on Figs. 9-12.  Fig. 3). The data on the process p + p → p + π + + π − + p at √ s = 200 GeV, |η π | < 1, |η ππ | < 2, p T π > 0.15 GeV, 0.005 < −t 1,2 < 0.03 GeV 2 , (STAR collaboration [17], [18]). a, b The result with all p p and π p rescattering corrections, c, d show the result when we try to fit the data by formulas without π p rescattering terms. Curves from up to down correspond to different values of the parameter Λ π in the off-shell pion form factor (11) We see an underestimation of the ISR data. For these low energies we have to take into account possible corrections to pion-proton amplitudes, since our approach describes data well only for energies greater than ∼ 3 GeV. In each shoulder (T πp amplitude in Fig. 3 RF) the energy can be less than 3 GeV.
As to the CDF data (Figs. 11, 12), which is overestimated for M ππ < 1.5 GeV, we can say that there are corrections (destructive interference terms) from reso-  Fig. 3). The data on the process p + p → p + π + + π − + p at √ s = 200 GeV, |η π | < 1, |η ππ | < 2, p T π > 0.15 GeV, 0.005 < −t 1,2 < 0.03 GeV 2 , (STAR collaboration [17], [18]). a, b The result with all p p and π p rescattering corrections, c, d show the result when we try to fit the data by formulas without π p rescattering terms. Curves from up to down correspond to different values of the parameter Λ π in the off-shell pion form factor (11) nances to the amplitude, like in Fig. 3 of [10], [11] and other effects for low M ππ , for example, the interference with γγ or γO fusion in the central production process. Also we have to take into account effects related to the irrelevance and possible modifications of the Regge approach (for the virtual pion exchange) in this kinematical region, as was discussed in the introduction.  Fig. 9 RF case of the model (see Fig. 3) with all rescattering corrections (a, b) and also when we try to fit the data by formulas without final pion-proton rescattering (c, d). The data on the process p+p → p+π + +π − +p at √ s = 63 GeV, |y π | < 1, ξ p > 0.9, (ISR and ABCDHW collaborations [15], [16]). Curves from up to down correspond to different values of the parameter Λ π in the off-shell pion form factor (11): a, b Λ π = 4, 3, 1.6, 1.2 GeV, c, d Λ π = 1.2, 1, 0.8 GeV.

CMS data and predictions
In Fig. 13 one can see the recent data from the CMS collaboration and curves of our model. The upper curve, which corresponds to the parameter Λ π , which better fits the STAR data on φ pp distribution (but gives higher values for M ππ > 1 GeV as depicted in Fig. 5b), also describes the data of CMS collaboration well (but overestimates the CDF data, as was shown in Figs. 11, 12). The lower curve underestimates the data from STAR, ISR and CMS, but it is close to the CDF data. Interference with resonances and modifications of the model in the region |t| ∼ŝ/2 ∼ 1 GeV 2 can change the picture, a) b) Fig. 10 RF case of the model (see Fig. 3) with all rescattering corrections (a) and also when we try to fit the data by formulas without final pion-proton rescattering (b). The data on the process p+p → p+π + +π − +p at √ s = 62 GeV, |y π | < 1.5, ξ p > 0.9, (ISR and ABCDHW collaborations [15], [16]). Curves from up to down correspond to different values of the parameter Λ π in the off-shell pion form factor (11): a Λ π = 4, 3, 1.6, 1.2 GeV, b Λ π = 1.2, 1, 0.8 GeV. a) b) Fig. 11 RF case of the model (see Fig. 3) with all rescattering corrections (a) and also when we try to fit the data by formulas without final pion-proton rescattering (b). The data on the process p +p → p + π + + π − +p at √ s = 1.96 TeV, |η π | < 1.3, |y ππ | < 1, p T ,π > 0.4 GeV, (CDF collaboration [19], [20]). Curves from up to down correspond to different values of the parameter Λ π in the off-shell pion form factor (11): a Λ π = 4, 3, 1.6, 1.2 GeV, b Λ π = 1.2, 1, 0.8 GeV.

Summary
After the experimental data analysis we have several facts: -In our approach the best description is given by the case RF (Fig. 3). That is why effects from rescattering (unitarity) corrections are very important. -The result is crucially dependent on the choice of Λ π in the off-shell pion form factor, i.e. ont (virtuality of the pion) dependence. -If we try to fit the data from STAR [17], [18], we find that the best description gives overestimation Fig. 14 RF case of the model (see Fig. 3) for the CMS energy. Curves from up to down correspond to the Born term, the amplitude with proton-proton rescattering corrections only and the one with all the corrections (proton-proton and pion-proton). Λ π = 3 GeV.
of the CDF data [19], [20] (especially in the region M ππ < 1.5 GeV) and underestimation of the ISR data [15], [16]. This is due to effects like the interference with resonance contributions or γγ → ππ and γO → ππ processes, effects related to the irrelevance and possible modifications of the Regge approach (for the virtual pion exchange) in this kinematical region, as was discussed in the introduction, corrections to pion-pion scattering at low M ππ , and corrections to T πp (s, t) for √ s < 3 GeV; -Predictions for CMS are close to the data, if we use the best fit to the STAR data on φ pp distribution (see Fig. 5b). We need also an estimate of the interference with resonance terms to see the full picture and draw final conclusions.

Pomeron-Pomeron to pion-pion cross-section
Another interesting question, which we can discuss here, concerns the Pomeron-Pomeron cross-section. As was shown in [3], it is possible to extract the Pomeronproton cross-section from the data on single (SD) and double (DD) dissociation, and this numerical value appears to be of the order of typical hadron-hadron crosssections. It was done by the use of covariant reggeization method with conserved spin-J meson currents, which helps to solve the old problem of the very small Pomeron-proton cross-section extracted by other authors [30]- [31]. Reggeon-hadron and reggeon-reggeon scattering can be considered as a scattering of all possible real mesons lying on the Regge trajectory of hadrons. Conceptually it is similar to hydrogen-hadron or hydrogen-hydrogen scattering, since hydrogen has the spectrum of the states, and each of them has its own probability to scatter on a hadron or another hydrogen atom. A specific feature is that we deal in this case with "off-shell atoms". There is some misunderstanding concerning the physical nature of reggeons. Actually, the reggeon is quite a general notion to qualitatively describe a generic quantum composite system. Conceptually a hadronic reggeon differs from the familiar hydrogen atom only by constituent content: electron-proton in the latter case and, say, quark-antiquark in the former. Energy levels of both are given by the corresponding Regge trajectories. So the reggeon is a full fledged composite particle in the same sense as an atom. Certainly the specific properties are different due to different binding forces.
As was shown also in [3], absorptive corrections play the crucial role in high energy scattering, and make the extraction procedure rather complicated and model dependent. We should propose some appropriate parametrization for Pomeron-hadron cross-section, then we apply unitarization procedure to obtain real SD or DD cross-sections.
It is possible to perform a similar procedure to extract the Pomeron-Pomeron cross-section. In further work we shall consider the extraction of the total Pomeron-Pomeron cross-section. Here we restrict ourselves by the extraction of the Pomeron-Pomeron to pion-pion one. Let us use the parametrization (4) for CEDP di-pion production to fix one parameter Λ π from the experimental data. After that we can simply estimate the Pomeron-Pomeron to di-pion cross-section (we use the covariant method from Appendix C): Here Π where all the functions are defined in Appendix B. If we use the approach (16), where all terms like in (50) are absorbed into the residue, we have to multiply the result by the additional factor F(t 1 ,t)F(t 2 ,t), where for conserved currents, and is the leading term for the case of non-conserved currents (see (53)).
Results of calculations (15)- (18) are shown in Fig. 15. In the case of conserved meson currents we obtain Pomeron-Pomeron to pion-pion cross-section 10 ÷ 100 times higher than in the case of non-conserved currents and with more specific and strong dependence on the Pomeron virtuality. In the old work [32]- [36] the extracted Pomeron-Pomeron total cross-section was of the order 100 µb at √ŝ < 3 GeV and almost independent on Pomeron's virtuality. σ PP→ππ should be at least less than this number. 1 Our calculations in the same kinematical region give numbers of the order 0.1 − 5 µb for non-conserved currents and 0.3−100 µb for the case of conserved currents. There are strong contributions of other processes (especially production of resonances) in this region, that is why we should have σ PP→ππ σ tot PP , which is obvious in the case of non-conserved currents where we have both extracted numbers to compare. 1 Let us note that authors of [32]- [36] extracted σ tot PP in the classical approach which corresponds to the case of nonconserved currents in the present paper, that is why we should compare their result with the solid curve on Fig. 15b, which gives σ PP→ππ ∼ 0.1÷5 µb for √ s < 3 GeV and |t| = 0.1 GeV 2 as found in [32]- [36]. To compare our results for the case of conserved currents we need to use the same approach also when extract σ tot In Sect. 2 it is shown that RF and, possibly, PF modes (see Fig. 3 for notations) give an appropriate description of the data, i. e. we have to take into account all rescattering corrections (even in the T πp amplitudes). If to use the RB mode, as some other authors do [7]- [13], 2 it is possible to extract the Pomeron-Po-meron cross-section more easily ("almost model independent method", as was done, for example, for the pion-proton cross-section [37]).

Conclusions
In this paper we have considered the process LM CEDP of di-pions and its description in the framework of the Regge-eikonal approach. Here we summarize all the facts and conclusions: -After calculations of several cases (see Fig. 3) we can see that the RF case is better suited to describe the data on the process p + p → p + π + + π − + p. -When we try to fit the data from STAR collaboration [17], [18] with different values of Λ π (different behavior of the virtual pion form factor), we obtain an underestimation of the ISR data [15], [16] and an overestimation of the CDF [19], [20]. Possible explanations are: interference terms with resonances, corrections to T πp (s, t) for √ s < 3 GeV, off-shell pion effects, and some other mechanisms in Pomeron-Pomeron to pion-pion process at low M ππ .
-We have rather good predictions to the CMS data, when we use the fit to the STAR data on φ pp distribution depicted on Fig. 5b, but we have to take into account interference with resonances to see the full picture. These main open problems regarding the model parameters are related to interference terms (we have to know all couplings of pions to resonances), which require full spectrum simulation comparisons with data simultaneously in several differential observables. This should be done in further work. -After estimations of the Pomeron-Pomeron to pionpion cross-section in the framework of covariant reggeization approach we obtain cross-sections which are (at least in the case of non-conserved Pomeron currents) much lower (∼ 0.1 ÷ 5 µb) than the total Pomeron-Pomeron cross-section, estimated by other authors [32]- [36] (∼ 100−300 µb), which shows strong contributions from other processes (especially from resonance production). For the case of conserved currents we need re-evaluate old results [32]- [36] in the same approach.
In further work we will take into account possible modifications of the model (amplitudes for resonances in LM CEDP and their interference with di-pion one, pion-pion cross-section, additional off-shell effects in subamplitudes and so on) for the best description of the data. This model will be implemented to the Monte-Carlo event generator ExDiff [38]. It is possible to calculate LM CEDP for other di-hadron final states (pp for "Odderon" hunting, K + K − , ηη and so on), which are also very informative for our understanding of diffractive mechanisms in strong interactions.
Here the y i (η i ) are the rapidities (pseudorapidities) of the final pions.
The phase space of the process in terms of the above variables is the following p i⊥ dp i⊥ dφ i · dy 3 dy 4 · J ; where the p i⊥ = |p i |,p 1,2z are the appropriate roots of the system Here λ(x, y, z) = x 2 + y 2 + z 2 − 2xy − 2xz − 2yz, and then J = λ 1/2 0 /2. For the differential cross-section we have dσ 2→4 3 i=1 dp i⊥ dφ i · dy 3 dy 4 The pseudorapidity is more convenient experimental variable, and we can use the transform to get the differential cross-section in the pseudorapidities.
T pp (s, q 2 ) = T πp (s, q 2 ) = Here we take S πp (s, t) = S π + p (s, t) = S π − p (s, t) and T el π + p (s, t) = T el π − p (s, t) = 4πsT πp (s, t) The approach (28) describes the data on pion-proton scattering better even at low energies, that is why we use it instead of the one presented in [23].
The functionsT pp andT πp are convenient for numerical calculations, since the oscillations are not so strong. and the same conditions are imposed on each group of indices in the tensors W and Π. Let us note (as was done in [39]), that conditions (42)-(44) are valid only on the mass shell of the spin-J meson. And when we go to the phase space of the scattering region, these conditions could be relevant only for conserved hadronic currents. However, this may well not be the case.