Midrapidity inclusive densities in high energy $pp$ collisions in additive quark model

High energy (CERN SPS and LHC) inelastic $pp$ ($p\bar p$) scattering is treated in the framework of additive quark model together with Pomeron exchange theory. We extract the midrapidity inclusive density of the charged secondaries produced in a single quark-quark collision and investigate its energy dependence. Predictions for the $\pi p$ collisions are presented.


Introduction
Regge theory provides a useful tool for the phenomenological description of high energy hadron collisions [1,2,3,4]. The quantitative predictions of Regge calculus are essentially dependent on the assumed coupling of participating hadrons to the Pomeron. In our previous papers [5,6] we described elastic pp (pp) scattering and diffractive dissociation processes including the recent LHC data in terms of a simple Regge exchange approach in the framework of the additive quark model (AQM) [7,9], or constituent quarks model as it is also referred to. It has been successfully applied to pp scattering processes at LHC energies [8]. In the present paper we consider the inclusive densities of produced secondaries in the midrapidity region in the same approach.
In AQM baryon is treated as a system of three spatially separated compact objects -the constituent quarks. Each constituent quark is colored and has an internal quarkgluon structure and a finite radius that is much less than the radius of the proton, r 2 q ≪ r 2 p . This picture is in good agreement both with SU(3) symmetry of the strong interaction and the quark-gluon structure of the proton [10,11,12]. The constituent quarks play the roles of incident particles in terms of which pp scattering is described in AQM.
In the case of inelastic pp collisions the secondary particles are produced in AQM in one or several qq collisions, so it opens the possibility to investigate inclusive densities of the secondaries produced in the single qq collision at the different initial energies. After that we can calculate the central inclusive densities in πp collisions without any new parameters.

High energy pp interactions in AQM
Elastic amplitudes for large energy s = (p 1 + p 2 ) 2 and small momentum transfer t are dominated by Pomeron exchange. We neglect the small difference in pp and pp scattering coming from the exchange of negative signature Reggeons, Odderon (see e.g. [13] and references therein), ω-Reggeon etc., since their contributions are suppressed by s.
The single t-channel exchange results into the amplitude of constituent quarks scattering, where α P (t) = α P (0) + α ′ P · t is the Pomeron trajectory specified by the intercept and slope values α P (0) and α ′ P , respectively. The Pomeron signature factor, determines the complex structure of the amplitude. The factor γ qq (t) = g 1 (t) · g 2 (t) has the meaning of the Pomeron coupling to the beam and target particles, the functions g 1,2 (t) being the vertices of the constituent quark-Pomeron interaction (filled circles in Fig. 1). It is worth to emphasize that the qq interaction is described here by single effective Pomeron exchange between each qq pair. Generally it may include the contributions of several Gribov bare Pomerons [14] and the parameters of the effective Pomeron could be different from those of the bare Pomerons. At the same time the one quark interaction with the two different target quarks is mediated by the exchange of the two effective Pomerons. In this respect there is a close resemblance between the nucleon-nucleon scattering in AQM and the nucleus-nucleus scattering in Glauber theory. The qq interaction plays the same role in the first case as the NN interaction in the second.
The elastic pp scattering amplitude (or pp scattering amplitude, here we do not distinguish between the two) is basically expressed as In this formula ψ(k i ) ≡ ψ(k 1 , k 2 , k 3 ), is the initial proton wave function in terms of the quarks' transverse momenta stands for the multipomeron exchange, Q k and Q ′ l are the momenta transferred to the target quark k or beam quark l by the Pomerons attached to them, Q is the total momentum transferred in the scattering, The scattering amplitude is presented in AQM as a sum over the terms with a given number of Pomerons, where the amplitudes M (n) pp collect all diagrams comprising various connections of the beam and target quark lines with n Pomerons. Similar to Glauber theory [15,16] one has to rule out the multiple interactions between the same quark pair. AQM permits the Pomeron to connect any two quark lines only once. It crucially decreases the combinatorics, leaving the diagrams with no more than n = 9 effective Pomerons. Several AQM diagrams are shown in Fig. 1. In the following we assume the Pomeron trajectory to have the simplest form, The value r 2 q defines the radius of the quark-quark interaction, while S 0 = (9 GeV) 2 has the meaning of typical energy scale in Regge theory.
In the first order there are nine equal quark-quark contributions due to one Pomeron exchange between qq pairs. The amplitude (2) reduces to a single term with expressed through the overlap function The function F P (Q, 0, 0) plays the role of a proton form factor for the strong interaction in AQM.
The quarks' wave function has been taken in the simple form of Gaussian packets, normalized to unity. The parametrization by the single exponent is unable to reproduce the minimum in dσ/dt distribution evidently seen in the experimental data for √ s = 7 TeV [17,18]. The two exponential fit used in our previous papers [5,6] reproduces this minimum but gives too low values of dσ/dt at |t| ∼ 0.7 − 0.8 Gev 2 . In the present paper the wave function is parameterized by the sum of three Gaussian exponents, which allows for the better description of dσ/dt.

Now the parameters read
One has to remark here that we do not claim the real matter distribution inside the proton to be close to the Gaussian shape, this form is suitable only to perform all the integrals analytically. The value a 1 is quite compatible with the large proton size assumed above, whereas a 2,3 values manifest the presence of the small radius components in the proton wave function. However, their relative weights are small so the total wave function (6) matches the condition r 2 p ≫ r 2 q , which is effectively fulfilled for the mean radii that are important for the calculations validity.
The higher orders elastic terms are expressed through the functions (5) integrated over Pomerons' momenta, × e −r 2 The sum in this formula refers to all distinct ways to connect the beam and target quark lines with n Pomerons in the scattering diagram. The set of momenta Q i and Q ′ l the quarks acquire from the attached Pomerons is particular for each connection pattern. A more detailed description can be found in [5].
With the amplitude (3) the differential cross section in the normalization adopted here is evaluated as The optical theorem, which relates the total cross section and the imaginary part of the amplitude, in this normalization reads The condition for the AQM applicability, r 2 q /r 2 p ≪ 1, holds rather well since r 2 p ≃ 12 Gev −2 whereas r 2 q ≃ 1.5 Gev −2 at √ s ≈ 7 TeV.
The resulting differential cross sections for pp scattering at √ s = 7 TeV are presented in Fig. 2 together with the predictions for πp scattering at the same energy (see below).

AGK cuts and inclusive densities in pp and qq interactions
All amplitudes of the inelastic processes in high energy pp collisions can be treated as the sum of various absorptive parts of elastic pp amplitude; see AGK cutting rules [20].
In the AQM the diagram with a single qq interaction, Fig. 3, has only one absorptive part. The one-Pomeron cut in the left hand side of Fig. 3 corresponds to the multiperipheral ladder of the produced secondaries in the right hand side of Fig. 3. The resulting cross section is σ (1) .
In the case of double interaction in Fig. 1b the imaginary part is given by the sum of three different absorptive parts presented in Fig. 4. The first one, the cut between Pomerons, is shown in Fig. 4a. It describes the elastic and diffractive dissociation processes without production of secondaries in the central (midrapidity) region. The second absorptive part, Fig. 4b, corresponds to the cut of one Pomeron and gives the first rescattering correction to the processes in Fig. 3. The multipheripheral ladders in Figs. 3 and 4b are practically the same and have the same midrapidity inclusive  [17,18]. Right Pseudorapidity distribution of the secondaries dN NSD /dη in pp scattering for the non-single diffractive events. The solid line shows the AQM estimates by Eq. (12). The experimental points have been taken from [19] densities dN qq /dy. The absorptive part in Fig. 4b has the numerical factor -4 due to the combinatorics [20]. The third absorptive part is shown in Fig. 4c, where the cut slices both Pomerons. It means the simultaneous production of two multipheripheral ladders of the secondaries. These ladders are also practically identical to those in Figs. 3 and 4b and result in equal inclusive densities dN qq /dy (we neglect the very small numerical difference coming from the energy conservation). The combinatorial factor here is +2 [20]. The contribution to the inclusive density of the secondaries from Fig. 4c is 4dN qq /dy. It is compensated by the negative contribution −4dN qq /dy from the process in Fig. 4b. Finally, the sum of all absorptive parts collected in Fig. 4 yields zero contribution to the inclusive densities of secondary particles in complete agreement with the AGK cutting rules [20].
Similarly, there is no contribution to the resulting inclusive density due to diagrams with a larger number of quark-quark interaction, therefore it is only the impulse approximation diagrams in Fig. 3 that provide the inclusive density of the secondaries produced in pp collisions in the midrapidity region, where σ (1) qq is the first order contribution to the total cross section. This equation is true as well for the pseudorapidity inclusive densities, which amounts to the replacement dy → dη.
The cross section of pp interaction used in (10) depends on the way the value of the inclusive density is fixed. It is determined as the number of the secondaries divided by the number of events in the small interval dy in the midrapidity region. In the diffractive dissociation the secondaries are produced practically only in the fragmentation regions; therefore, the number of secondary particles in the midrapidity region does not change whether or not we include diffractive dissociation events in our sample. However, the number of events, i.e. the denominator in the definition of dN/dy, differs for these two cases, so the net value dN pp /dy in the left hand side of (10) has to be multiplied by σ inel P P , if we take all inelastic events, or by σ nondif f r P P if we take the events without diffractive dissociation.
Let us try to estimate the energy dependence of dN qq /dy and dN qq /dη using the existing data. There are several available experimental points for dN pp /dη at the different energies measured in all inelastic events. They are shown in Fig. 5 together with the calculations of dN pp /dy and dN pp /dη in the quark-gluon string model (QGSM) [21,22]. Actually QGSM output is employed here only to extrapolate the existing experimental data.
The inelastic cross section entering Eq. (10) can be obtained from the identity where the total cross section is evaluated through the optical theorem (9) by summing up all nine orders of AQM diagrams (7), while the elastic cross section, σ elastic pp , is obtained by integrating differential cross section (8) over t. The cross section σ (1) qq is given by the first order of the AQM amplitude (4).
Equations (10) and (11) allow one to find dN qq /dη values presented in Fig. 5. At the LHC energies √ s > 0.9 TeV dN qq /dη becomes independent on the initial energy within our theoretical accuracy ∼ 10%.  [23,24]. Right The total cross sections of pp (solid line) and πp (dashed line) as the initial energy functions. The experimental pp points are taken from [25,26] Experimental papers often present the data for the energy dependence of the pseudorapidity distribution of the secondaries, dN NSD /dη, measured in the non-single diffractive events. It can be obtained by the formula where σ SD pp is the cross section of the single diffractive pp scattering (from one side). To make a quick estimate we have used σ SD pp values calculated in AQM in our previous paper [6]. The results are shown in Fig. 2 (right panel) together with the existing experimental points. We get a reasonable agreement, the ratio dN NSD dη / dN dη ∼ 1.1 ÷ 1.15 at the LHC energies.

Predictions for πp collisions
To obtain the predictions of midrapidity inclusive densities in πp collisions one needs to know the total πp cross section. It has not been measured experimentally at the very high energies but can be calculated in AQM. In our approach the interaction of quarks and antiquarks constituting the pion are the same as those in the proton (so far as only Pomeron exchange is encountered). The amplitude of the elastic πp collision is evaluated in AQM similarly to the elastic pp one, see (7), πp (s, t) = i n−1 γ qq η P (t n )e ∆·ξ n d 2 q 1 π · · · d 2 q n π π δ (2) (q 1 + . . . + q n − Q) (13) × e −r 2 Here F π and F P are the pion and proton form factors while all other variables are the same as those for pp scattering.
The quark combinatorics is more simple for πp collisions compared to pp case. In particular, there are only six orders of the admissible diagrams. The sum for the first order contribution reduces to a single term, 6F π (Q, 0) F P (Q, 0, 0), Q 2 = −t. The second order sum includes three types of diagrams, = 6 F π (Q, 0)F P (q 1 , q 2 , 0) + 3 F π (q 1 , q 2 )F P (Q, 0, 0) + 6 F π (q 1 , q 2 )F P (q 1 , q 2 , 0), where the first two terms come from the diagram with both Pomerons coupled to the same quark line in the pion (first term, Fig. 6a) and in the proton (second term, Fig. 6b); in the third term they connect different quark lines (Fig. 6c). The numerical coefficients encounter the number of connections resulting in equal expressions after variables changing in the integrals (13). The rest orders have a similar structure derived from the combinatorics to redistribute q 1 , . . . , q n momenta among Q i and Q ′ i groups. In the highest order, containing six effective Pomerons, = F π (q 1 + q 2 + q 3 , q 4 + q 5 + q 6 )F P (q 1 + q 4 , q 2 + q 5 , q 3 + q 6 ), each quark from the proton interacts with both the quarks from the pion.
The differential and the total πp cross sections are evaluated via Eq. (8) and the optical theorem (9) respectively. Our results for the elastic πp scattering show a minimum at √ s = 7 TeV placed at −t ≈ 0.65 GeV 2 (Fig 2). The ratio of the total pp and πp cross section in the optical approximation of AQM is well known to be 3/2 [7]. With the multiple rescattering included this value changes depending on the ratio of proton and pion radii. The experimental data [27] gives for the ratio r 2 π /r 2 p ≈ 0.57, so we take the pion wave function in the same form, (6), rescaling all radius parameters as a π 1,2,3 = 0.57a p 1,2,3 . Actually the dependence σ tot πp on the parameters of a 1,2,3 and C 1,2 is rather weak. We get the ratio σ tot pp /σ tot πp ≈ 1.2 ÷ 1.3. Unfortunately there are no experimental data on the πp scattering at the LHC energies, the AQM results for them are presented in Fig. 5 together with the predictions for the pp case. Note here as well that AQM predictions for dσ pp /dt(t = 0) are in good agreement with the data [5].
The obtained values σ tot πp allow to find the midrapidity inclusive density in πp collisions. The results for dN πp /dη(η = 0) as a function of the initial energy are presented in Fig. 5. The obtained data can be used for the calculation of particle production at the very high energies; in particular, in cosmic ray physics.

Conclusion
In the framework of AQM we have extracted the inclusive density of the secondaries in qq interactions in the midrapidity region. We used these values to get prediction for πp collisions at high energies. These quantities can be useful to estimate the secondary production at the very high energies, say, in cosmic ray physics.
The applicability of AQM requires the contribution from the multipomeron qq interactions to be small compared to the interaction between different quarks responsible in this approach for the pp scattering. This is valid for the soft processes, whose amplitude is practically pure imaginary so that the qq cross section does not exceed the geometrical limit ∼ r 2 q . On the other hand there are additional combinatorial factors increasing the pp cross section, so it can always be assumed to be larger than the qq one. A reasonable description of the elastic pp scattering has been reached without appealing to the enhanced diagrams with interacted Pomerons. It provides evidence that AQM is at work up to LHC energies. However, for the energies essentially above the LHC ones the multipomeron interactions would begin to play an important role, which could modify our results for asymptotically high energies. This work has been supported by RSCF grant No 14 -22 -00281.