A CGC/saturation approach for angular correlations in proton-proton scattering

TL;DR

This paper extends a model to calculate azimuthal angle correlations in proton-proton scattering, showing that Bose-Einstein correlations can explain observed angular correlations without invoking quark-gluon plasma formation.

Contribution

The authors generalize their model to describe azimuthal correlations in proton-proton collisions, emphasizing the role of Bose-Einstein correlations over final state interactions.

Findings

Calculated $v_n$ values agree with experimental data

Energy and multiplicity independence of correlations

Bose-Einstein correlations explain angular correlations

Abstract

We generalized our model for the description of hard processes, and calculate the value of the azimuthal angle correlations ( Fourier harmonics $v_n$ ), for proton-proton scattering.The energy and multiplicity independence, as well as the value of $v_n$, turn out the be in accord with the experimental data, or slightly larger. Therefore, before making extreme assumptions on proton-proton collisions, such as the production of quark-gluon plasma in the large multiplicity events, we need to explain how these effect theBose-Einstein correlations which are so large, that have to be taken into account, and which are able to describe the angular correlations in proton-proton collisions, without including final state interactions

Figures

Tables

Table 1. Table 1: Fitted parameters of the model. Fit I: parameters for the soft interaction at high energy are taken from Ref. [ 25 ] . The additional parameters for DIS were found by fitting to the F 2 subscript 𝐹 2 F_{2} structure function (see below). Fit II: joint fit to the soft interaction data at high energy and the DIS data.

model	$\lambda$	${\varphi }_0$ ($Ge{V}^{-2}$)	$g_1$ ($Ge{V}^{-1}$)	$g_2$ ($Ge{V}^{-1}$)	$m\left(GeV\right)$	$m_1\left(GeV\right)$	$m_2\left(GeV\right)$	$\beta$
I(soft int.)	0.38	0.0019	110.2	11.2	5.25	0.92	1.9	0.58
II:(soft + DIS)	0.38	0.0022	96.9	20.96	5.25	0.86	1.76	0.66

Table 2. Table 2: The values of cross sections versus energy. σ sd smd subscript superscript 𝜎 smd sd \sigma^{\rm smd}_{\rm sd} and σ dd smd subscript superscript 𝜎 smd dd \sigma^{\rm smd}_{\rm dd} denote the cross sections for diffraction dissociation in the small mass region, for single and double diffraction, which stem from the Good-Walker mechanism. While σ sd lmd subscript superscript 𝜎 lmd sd \sigma^{\rm lmd}_{\rm sd} and σ dd lmd subscript superscript 𝜎 lmd dd \sigma^{\rm lmd}_{\rm dd} denote high mass diffraction, coming from the dressed Pomeron contributions. The predictions of fit II, are shown in brackets.

W	${\sigma }_{tot}$	${\sigma }_{el}$(mb)	$B_{el}$	single	diffraction	double	diffraction
(TeV)	(mb)	(mb)	$\left(Ge{V}^{-2}\right)$	${\sigma }_{\mathrm{sd}}^{\mathrm{smd}}$ (mb)	${\sigma }_{\mathrm{sd}}^{\mathrm{lmd}}$ (mb)	${\sigma }_{\mathrm{dd}}^{\mathrm{smd}}$ (mb)	${\sigma }_{\mathrm{dd}}^{\mathrm{lmd}}$ (mb)
0.576	62.3(60.7)	12.9(13.1)	15.2(15.17)	5.64(4.12)	1.85(1.79)	0.7(0.39)	0.46 (0.50)
0.9	69.2(68.07)	15(15.05)	16(15.95)	6.254.67)	2.39(2.35)	0.77(0.46)	0.67(0.745)
1.8	79.2(78.76)	18.2(19.1)	17.1(17.12)	7.1(5.44)	3.35(3.28)	0.89(0.56)	1.17 (1.30)
2.74	85.5(85.44)	20.2(21.4)	17.8(17.86)	7.6(5.91)	4.07(4.02)	0.97(0.63)	1.62(1.79)
7	99.8(100.64)	25(26.7)	19.5(19.6)	8.7(6.96)	6.2(6.17)	1.15(0.814)	3.27(3.67)
8	101.8(102.8)	25.7(27.4)	19.7(19.82)	8.82(7.1)	6.55(6.56)	1.17(0.841)	3.63(4.05)
13	109.3(111.07)	28.3(30.2)	20.6(20.74)	9.36(7.64)	8.08(8.11)	1.27(0.942)	5.11(5.74)
14	110.5(111.97)	28.7(30.6)	20.7(20.88)	9.44(7.71)	8.34(8.42)	1.27(0.96)	5.4(6.06)
57	131.7(134.0)	36.2(38.5)	23.1(23.0)	10.85(9.15)	15.02(15.01)	1.56(1.26)	13.7(15.6)

Table 3. Table 3: Fitted parameters for DIS. The description of fit I and fit II is given in section 2.3, and in the caption of Table 1.

model	$m_q(GeV$	$m_c\left(GeV\right)$	${\alpha }_S\left(\mu \right)$	$\mu$ ($GeV$)	$A_{IR}(Ge{V}^2$	${\alpha }_{IR}\left(0\right))$
I	0.3	1.25	0.263	1.2	2.34	0.55
II	0.2	1.2	0.34	1.25	5.44	0.56

Equations

$$Y\le \frac{2}{{\Delta }_{\text{mbox}BFKL}}\ln \left(\frac{1}{{\Delta }_{\text{mbox}BFKL}^2}\right)$$

$$G_{IP}^{\text{mbox}dressed}\left(Y,r,R;b\right)=$$

$$\displaystyle\,\,\int\prod_{i=1}d^{2}r_{i}\,d^{2}b_{i}\,d^{2}r^{\prime}_{i}\,d^{2}b^{\prime}_{i}\,N\left(Y-Y^{\prime},r,\{r_{i},b-b_{i}\}\right)\framebox{$A^{\rm BA}_{\mbox{\tiny dipole-dipole}}\left(r_{i},r^{\prime}_{i},\vec{b}_{i}-\vec{b^{\prime}}_{i}\right)$}\,N\left(Y^{\prime},R,\{r^{\prime}_{i},b^{\prime}_{i}\}\right)$$

$$N\left(Y-Y^{\prime },r,\{r_i,b_i\}\right)=N^{\mathrm{BK}}\left(Y-Y^{\prime },r,\{r_i,b_i\}\right)$$

$$=\sum _{n=1}^{\infty }{\left(-1\right)}^{n+1}{\tilde{C}}_n\left(r\right)\prod _{i=1}^n{G}_{IP}\left(Y-Y^{\prime };r,r_i,b_i\right)=\sum _{n=1}^{\infty }{\left(-1\right)}^{n+1}{\tilde{C}}_n\left(r\right)\prod _{i=1}^n{G}_{IP}\left(z-z_i\right).$$

$$G_{IP}\left(z_i\right)={\varphi }_0{\left(r_i^2{Q}_s^2\left(Y,b_i\right)\right)}^{1-{\gamma }_{cr}}$$

$$N\left(G_{IP}\left({\varphi }_0,z\right)\right)=\sum _{n=1}^{\infty }{\left(-1\right)}^{n+1}{C}_n\left({\varphi }_0\right)G_{IP}^n\left({\varphi }_0,z\right).$$

$${\tilde{C}}_n\left(r\right)=C_n\left({\varphi }_0\right).$$

$$N^{\mathrm{BK}}\left(G_{IP}\left({\varphi }_0,z\right)\right)=a\left(1-\exp \left(-G_{IP}\left({\varphi }_0,z\right)\right)\right)+\left(1-a\right)\frac{G_{IP}\left({\varphi }_0,z\right)}{1+G_{IP}\left({\varphi }_0,z\right)},$$

$$\frac{{\alpha }_S^2}{4\pi }{G}_{IP}\left(Y-0,r,R;b\right)=$$

$$\int d^2{r}^{\prime }{d}^2{b}^{\prime }{d}^2{r}^{\prime \prime }{d}^2{b}^{\prime \prime }{G}_{IP}\left(Y-Y^{\prime },r,r^{\prime },\vec{b}-{\vec{b}}^{{}^{\prime }}\right)G_{IP}\left(Y^{\prime }{r}^{\prime \prime },R,{\vec{b}}^{{}^{\prime \prime }}\right)A_{\text{mbox}dipole-dipole}^{\mathrm{BA}}\left(r^{\prime },r^{\prime \prime },\overrightarrow{b^{\prime \prime }}-\overrightarrow{b^{\prime }}\right)$$

$$G^{\text{mbox}dressed}\left(T\right)$$

$$G\left(T\right)=1-\frac{1}{T}\exp \left(\frac{1}{T}\right)\Gamma \left(0,\frac{1}{T}\right)$$

$$T\left(r_{\perp },Y=\ln \left(s/s_0\right),b\right)={\varphi }_0{\left(r_{\perp }^2{Q}_s^2\left(Y,b\right)\right)}^{\bar{\gamma }}$$

$$Q_s^2\left(b,Y\right)=Q_{0s}^2\left(b,Y_0\right)e^{\lambda (Y-Y_0)}$$

$$Q_{0s}^2\left(b,Y_0\right)={\left(m^2\right)}^{1-1/\bar{\gamma }}{\left(S\left(b,m\right)\right)}^{1/\bar{\gamma }}S\left(b,m\right)=\frac{m^2}{2\pi }{e}^{-mb}\text{mbox}and\bar{\gamma }=0.63$$

$${\psi }_h=\alpha {\Psi }_1+\beta {\Psi }_2;{\psi }_D=-\beta {\Psi }_1+\alpha {\Psi }_2;\text{mbox}where{\alpha }^2+{\beta }^2=1;$$

$$A_{i,k}^{i^{\prime }{k}^{\prime }}=<{\psi }_i{\psi }_k|\mathbf{T}|{\psi }_{i^{\prime }}{\psi }_{k^{\prime }}>=A_{i,k}{\delta }_{i,i^{\prime }}{\delta }_{k,k^{\prime }}.$$

$$2\text{mbox}Im{A}_{i,k}\left(s,b\right)=|A_{i,k}\left(s,b\right){|}^2+G_{i,k}^{in}(s,b),$$

$$A_{i,k}(s,b)=i\left(1-\exp \left(-{\Omega }_{i,k}(s,b)\right)\right),$$

$$G_{i,k}^{in}(s,b)=1-\exp \left(-2{\Omega }_{i,k}(s,b)\right).$$

$${\Omega }_{i,k}(r_{\perp },Y-Y_0,b)=\int d^2{b}^{\prime }{d}^2{b}^{\prime \prime }{g}_i\left({\vec{b}}^{\prime },m_i\right)G^{\text{mbox}dressed}\left(T\left(r_{\perp },Y-Y_0,{\vec{b}}^{\prime \prime }\right)\right)g_k\left(\vec{b}-{\vec{b}}^{\prime }-{\vec{b}}^{\prime \prime },m_k\right)$$

$$G_{3{I\!\!P}}\Big{/}g_{i}(b=0)\,\ll\,\,1;~{}~{}~{}~{}~{}~{}~{}~{}m\,\gg\,m_{1}~{}\mbox{and}~{}m_{2}$$

$${\Omega }_{i,k}(r_{\perp },Y-Y_0,b)$$

$$\Omega \left(r,Y-Y_0;b\right)$$

$$g_i\left(b\right)$$

$$S_p\left(b,m_i\right)=\frac{1}{4\pi }{m}_i^{3}b{K}_1\left(m_{i}b\right)\stackrel{\text{mbox}Fourierimage}{\to }\frac{1}{{\left(1+Q_T^2/m_i^2\right)}^2}$$

$${\tilde{G}}^{\text{mbox}dressed}\left(r,Y-Y_0\right)=\int d^{2}b{G}^{\text{mbox}dressed}\left(T\left(r,Y-Y_0,b\right)\right)$$

$$N^{BK}\left(G_{IP}^i\left(r_{\perp },Y,b\right)\right)=a\left(1-\exp \left(-G_{IP}^i\left(r_{\perp },Y,b\right)\right)\right)+\left(1-a\right)\frac{G_{IP}^i\left(r_{\perp },Y,b\right)}{1+G_{IP}^i\left(r_{\perp },Y,b\right)},$$

$$N_{T,L}\left(Q,Y;b\right)=\int \frac{d^{2}r}{4\pi }{\int }_0^{1}dz|{\Psi }_{T,L}^{{\gamma }^{\ast }}\left(Q,r,z\right){|}^{2}N\left(r,Y;b\right)$$

$$N\left(r,Y;b\right)={\alpha }^2{N}_1^{BK}\left(g_{1}S\left(b,m_1\right){\tilde{G}}_{IP}\left(r;Y\right)\right)+{\beta }^2{N}_2^{BK}\left(g_{2}S\left(b,m_2\right){\tilde{G}}_{IP}\left(r;Y\right)\right)$$

$$({\Psi }^{\ast }\Psi {)}_T^{{\gamma }^{\ast }}$$

$$({\Psi }^{\ast }\Psi {)}_L^{{\gamma }^{\ast }}$$

Full text

Soft Pomeron, BFKL Pomeron, Diffractive Cross Sections, Rapidity Correlations, Hard processes, Azimuthal angle Correlations

††dedicated: PACS: 13.85.-t, 13.85.Hd, 11.55.-m, 11.55.Bq††preprint: TAUP - 3008/16

1 Introduction

The experimental data on azimuthal angle correlations (see Refs. [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11] show suprizing similarities between different processes: nucleus-nucleus, hadron-nucleus and hadron-hadron collisions. The popular explanation is related to elliptic flow, and stems from the interaction in the final state. In the framework of such an approach, we have to assume that the proton interactions are similar to nucleus scattering, at least for events with large multiplicity. However, the ATLAS data [9] show that $v_{2,2}$,$v_{3,3}$ and $v_{4,4}$ do not depend on multiplicity at W=13 TeV and at W= 2.76 TeV.

In this paper we will discuss these correlations from a different point of view. We believe that the general origin of the azimuthal angle correlations in all reactions, stems from the Bose-Einstein correlations (BEC) of the produced gluons, which originate from the gluon wave function in the initial state[12, 13, 15, 14, 16, 17, 21, 18, 19, 20]. The attractive feature of this idea is that, BEC have a general source that characterizes the volume of the interaction[22, 23]. Therefore, the main dimensional parameters of the interaction that manifest themselves in diffraction scattering, and in inclusive production, should determine the BEC. In other words, that in spite of the embryonic stage of our understanding of the confinement of quarks and gluons, we can develop a quantitative approach for the BEC in the framework of a model for soft interactions at high energy. To accomplish this, we need to construct such a model, that will allow us to discuss soft and hard processes on the same footing.

The main goal of this paper is to develop such model. Fortunately, we have built a model which provides a good description of all the soft data[24, 25, 26, 27, 28, 29], including, total, inelastic, elastic and diffractive cross sections, the $t$-dependence of these cross sections as well as the inclusive production and rapidity correlations. In this paper we expand this model to include the hard interactions mostly using the geometric scaling behavior of the scattering amplitude[30] for the hard kinematic region in the Colour Glass Condensate (CGC)/saturation approach[31].

The idea of BEC being the main source of the azimuthal angle correlations, is marred by the observation[15, 16], that the process of the central diffraction production of colourless gluon dijets, gives a contribution which is equal to that of the BEC. In this case $v_{n,n}$ with odd $n$ are equal to zero, while $v_{n,n}$ with even $n$, are twice larger. We will not discuss this problem here. Our main goal is to obtain reliable quantitative estimates for $v_{n,n}$. However, we believe that due to the Sudakov suppression in Double Log Approximation of perturbative QCD, the dijet contribution is negligibly small [32].

The paper is organized as follows. In the next section we give a brief review of our model which is based on CGC/saturation approach. We discuss what we have taken from the theory in our approach, and what we have considered from a pure phenomenological approach. We will attempt to clarify the physical meaning of the introduced phenomenological parameters, and show how we include three dimensional sizes, which have been used to describe the scattering amplitude. In the language of the Constituent Quark Model these three sizes are the hadron radius, the size of the constituent quark, and the saturation momentum, which is a typical scale for the high energy amplitude.

In section 3 we generalize our model including the construction of an amplitude at short distances, which is able to describe the deep inelastic scattering (DIS). We compare our amplitude with HERA experimental data [33]. In section four we calculate the inclusive cross section, and show that we obtain good agreement with the experimental data. This is very important for our calculation, since it demonstrates that we are able to describe the experimental data for inclusive production not only at short distances, but also at long distances.

In section 5 we calculate the value of BEC correlations, its energy and multiplicity dependence. We obtain the value of $v_{n}$ which are a bit larger than the experimental ones, with a mild dependence on energy and multiplicity. We consider these estimates as the first quantitative prediction for $v_{n}$ in proton-proton scattering, which are in agreement with the values of the inclusive cross sections, and the cross sections for the hard processes.

In section 6 we draw our conclusions and outline the problems for future investigations.

2 The model

2.1 Theoretical input from the CGC/saturation approach

In this section we generalize our model for soft interactions at high energy[24, 25, 26, 27, 28, 29] to include a description of hard processes. This model incorporates two ingredients: the achievements of the CGC/saturation approach, which is an effective theory for QCD at high energy; and the pure phenomenological treatment of the long distance non-perturbative physics, due to the lack of the theoretical understanding of confinement of quark and gluons.

We wish to stress that the most of this section does not contain new results, it reviews our approach, and it is included in the paper only for the sake of completeness of presentation. One can find more details in Refs.[24, 25, 26, 27, 28, 29].

The effective theory for QCD at high energies exists in two different formulations: the CGC/saturation approach [34, 35, 36, 37], and the BFKL Pomeron calculus [38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51]. In building our model we rely on the BFKL Pomeron calculus, as the relation to diffractive physics and soft processes in general, is more transparent in this approach. However, we believe the CGC/saturation approach produces a more general pattern [49, 50] for the treatment of high energy QCD. Fortunately, in Ref.[50] it was shown, that these two approaches are equivalent for

[TABLE]

where $\Delta_{\mbox{\tiny BFKL}}$ denotes the intercept of the BFKL Pomeron. As we will see, in our model $\Delta_{\mbox{\tiny BFKL}}\,\approx\,0.2-0.25$ leading to $Y_{max}=20-30$, which covers all accessible energies.

The main ingredient, that we need to find, is the resulting (dressed) BFKL Pomeron Green function, which can be calculated using $t$-channel unitarity constraints:

[TABLE]

where $N\left(Y-Y^{\prime},r,\{r_{i},b-b_{i}\}\right)$ denotes the amplitude for the production in the $t$-channel of the set of dipoles with $Y=Y^{\prime}$ and with the size $r_{i}$, at the impact parameters $b_{i}$. $A^{BA}_{\mbox{\tiny dipole-dipole}}$ denotes the dipole-dipole scattering amplitude in the Born approximation of perturbative QCD, which are indicated by red circles in Fig. 1-a. In addition, in Ref.[50] it is shown that for such $Y$, we can safely use the Mueller-Patel-Salam-Iancu (MPSI) approach[52]. In this approximation we estimate the amplitudes $N$ in Eq. (2.2), using BFKL Pomeron ’fan’ diagrams (see Fig. 1-a for examples of such diagrams). In other words, we can use the parton cascade of the Baslitsky-Kovchegov[36] equation, to find the amplitude for the production of dipoles of size $r_{i}$ at impact parameters $b_{i}$. This amplitude can be written as (see Fig. 1-c)

[TABLE]

$G_{I\!\!P}$ denotes the Green function of the BFKL Pomeron. In the last equation we used the fact that in the saturation region this Green function has geometric scaling behavior, and so it depends on one variable: $z_{i}\,=\,\ln\left(Q^{2}_{s}(Y^{\prime})r^{2}_{i}\right)$, where $Q_{s}\left(Y^{\prime}\right)$, is the saturation scale, in the vicinity of the saturation scale[53]

[TABLE]

where $\gamma_{cr}=0.37$.

In Ref.[51], it was shown that, the solution to the non-linear BK equation has the following general form

[TABLE]

Comparing Eq. (2.3) with Eq. (2.5) we see

[TABLE]

Coefficients $C_{n}$ can be determined from the solution to the Balitsky-Kovchegov equation [36], in the saturation region. The numerical solution has been found in Ref.[51] for the simplified BFKL kernel in which only the leading twist contribution was taken into account:

[TABLE]

with $a$ = 0.65. Eq. (2.7) is a convenient parameterization of the numerical solution, with an accuracy of better than 5%. Having $C_{n}$ we can calculate the Green function of the dressed BFKL Pomeron using Eq. (2.2), and the property of the BFKL Pomeron exchange:

[TABLE]

Carrying out the integrations in Eq. (2.2), we obtain the Green function of the dressed Pomeron in the following form:

[TABLE]

where $\Gamma\left(s,z\right)$ is the upper incomplete gamma function (see Ref.[54] formula 8.35), and $T$ denotes the BFKL Pomeron in the vicinity of the saturation scale ( see Eq. (2.4))

[TABLE]

The Green function of Eq. (2.1) depends on the size of the dipoles, and we will use it for discussing the hard processes. In our analysis of the soft interaction we fixed $r=1/m$, and $m$ was a fitting parameter.

2.2 Phenomenology: assumptions and new small parameters

Unfortunately, due to the embryonic stage of theoretical understanding of the confinement of quarks and gluons, it is necessary to use pure phenomenalogical ideas to fix two major problems in high energy scattering: the structure of hadrons, and the large impact parameter behavior of the scattering amplitude[55]. The main idea to correct the large impact parameter behaviour, is to assume that the saturation momentum has the following dependence on the impact parameter $b$:

[TABLE]

with

[TABLE]

We have introduced a new phenomenological parameter $m$ to describe the large $b$ behaviour. The $Y$ dependence as well as $r^{2}$ dependence, can be found from CGC/saturation approach [31], since $\phi_{0}$ and $\lambda$ can be calculated in the leading order of perturbative QCD. However, since the higher order corrections turn out to be large [56], we treat them as parameters to be fitted. $m$ is a non-perturbative parameter, which determines the typical sizes of dipoles within the hadrons. In Table 1, we show that from the fit, $m$ = 5.25 GeV, supporting our main assumption that we can apply the BFKL Pomeron calculus, based on perturbative QCD, to the soft interaction since $m\,\gg\,\mu_{soft}$, where $\mu_{soft}$ is the scale of soft interaction, which is of the order of the mass of pion or $\Lambda_{\rm QCD}$.

The idea to absorb the non-perturbative $b$ dependence into the saturation scale, stems both from the success of this idea in the description of the hard processes in framework of the saturation model [57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 77, 76], and from the semi-classical solution to the BK equation[78], as well as from the analytical solution deep in the saturation domain[79].

The second unsolved problem for which we need a phenomenological input, is the structure of the scattering hadrons. We use a two channel model, which allows us to calculate the diffractive production in the region of small masses. In this model, we replace the rich structure of the diffractively produced states, by a single state with the wave function $\psi_{D}$, a la Good-Walker [80]. The observed physical hadronic and diffractive states are written in the form

[TABLE]

Functions $\psi_{1}$ and $\psi_{2}$ form a complete set of orthogonal functions $\{\psi_{i}\}$ which diagonalize the interaction matrix $T$

[TABLE]

The unitarity constraints take the form

[TABLE]

where $G^{in}_{i,k}$ denotes the contribution of all non diffractive inelastic processes, i.e. it is the summed probability for these final states to be produced in the scattering of a state $i$ off a state $k$. In Eq. (2.15) $\sqrt{s}=W$ denotes the energy of the colliding hadrons, and $b$ the impact parameter. A simple solution to Eq. (2.15) at high energies, has the eikonal form with an arbitrary opacity $\Omega_{ik}$, where the real part of the amplitude is much smaller than the imaginary part.

[TABLE]

[TABLE]

Eq. (2.17) implies that $P^{S}_{i,k}=\exp\left(-2\,\Omega_{i,k}(s,b)\right)$, is the probability that the initial projectiles $(i,k)$ reach the final state interaction unchanged, regardless of the initial state re-scatterings.

The first approach is to use the eikonal approximation for $\Omega$ in which

[TABLE]

where $m_{i}$ denote the masses, which is introduced phenomenologically to determine the $b$ dependence of $g_{i}$ (see below).

We propose a more general approach, which takes into account the new small parameters, that are determined by fitting to the experimental data (see Table 1 and Fig. 1 for notation):

[TABLE]

The second equation in Eq. (2.19) leads to the fact that $b^{\prime\prime}$ in Eq. (2.18) is much smaller than $b$ and $b^{\prime}$, therefore, Eq. (2.18) can be re-written in a simpler form

[TABLE]

Using the first small parameter of Eq. (2.19), we see that the main contribution stems from the net diagrams shown in Fig. 1-b. The sum of these diagrams[25] leads to the following expression for $\Omega_{i,k}(s,b)$

[TABLE]

where

[TABLE]

[TABLE]

where $T\left(r,Y-Y_{0},b\right)$ is given by Eq. (2.10).

The impact parameter dependence of $S_{p}\left(b,m_{i}\right)$ is purely phenomenological, however, Eq. (2.23) which has a form of the electromagnetic proton form factor, leads to the correct ($\exp\left(-\mu b\right)$) behavior at large $b$[82], and has correct behavior at large $Q_{T}$, which has been calculate in the framework of perturbative QCD [83]. We wish to draw the reader’s attention to the fact that $m_{1}$ and $m_{2}$ are the two dimensional scales in a hadron, which in the framework of the constituent quark model, we assign to the size of the hadron ($R_{h}\propto 1/m_{1}$), and the size of the constituent quark ($R_{Q}\propto 1/m_{2}$).

Note that $\tilde{G}^{\mbox{\tiny dressed}}\left(Y-Y_{0}\right)$ does not depend on $b$. In all previous formulae, the value of the triple BFKL Pomeron vertex is known: $G_{3{I\!\!P}}=1.29\,GeV^{-1}$.

For further discussion, we introduce the notation

[TABLE]

with $a=0.65$ . Eq. (2.25) is an analytical approximation to the numerical solution for the BK equation[51]. $G^{i}_{I\!\!P}\left(r_{\bot},Y;b\right)\,=\,\,g_{i}\left(b\right)\,\tilde{G}^{\mbox{\tiny dressed}}\left(r_{\bot},Y-Y_{0}\right)$. We recall that the BK equation sums the ‘fan’ diagrams.

2.3 Results of the fit

In this paper we make two fits. In the first one (fit I in Table 1 and Table 3) we do not change the parameters that govern the soft interactions in our model, and are shown in Table 1. The additional parameters that we need for the description of the deep inelastic data, and which we will discuss in the next section (see Table 3), were fitted using the HERA data on the deep inelastic structure function $F_{2}$. The second fit, is a joint fit to the soft strong interaction data and the DIS data. In Fig. 2 we show the results of our model compared with the HERA data. The model predictions are in accord with the data for $0.85\leq Q^{2}\leq 27GeV^{2}$, while for higher values of $Q^{2}$ and of $x$, the model values are slightly larger than the data.

In Table 2 we present our predictions for the soft interaction observables, in general the values obtained in the model for the soft interactions agree with the published LHC data, as well as the new preliminary TOTEM values at W =2.7,7,8,13 TeV (see Ref.[81]). We are in very good agreement with the data for $\sigma_{tot}$, $\sigma_{el}$ and $B_{el}$. Regarding $\sigma_{sd}$ and $\sigma_{dd}$, a problem exists when attempting to compare with the experimetal results. This is due to the difficulties of measuring diffractive events at LHC energies, the different experiments have different cuts on the values of the diffractive mass measured, making it problematic when attempting to compare the model predictions with the experimental results.

In Table 2 we show the results of the two fits, the results are close to one another, the main difference shows up only at high energies. Indeed, in fit I the cross section for single diffraction is equal to 14.9 mb, while in fit II this value is smaller (13.1 mb). The smaller value of the diffraction cross sections is closer to TOTEM and CMS data.

3 Deep inelastic scattering

3.1 Generalities

In this section, we compare our amplitude with the experimental data on deep inelastic scattering (DIS). In the framework of our approach, the observables of DIS can be re-written using

[TABLE]

where $Y\,=\,\ln\left(1/x_{Bj}\right)$ and $x_{Bj}$ is the Bjorken $x$. $z$ is the fraction of energy carried by quark. $Q$ is the photon virtuality. $b$ denotes the impact parameter for the scattering of the colorless dipole of size $r$ with the proton. $N\left(r,Y;b\right)$ is the scattering amplitude of this dipole, which in our model can be written in the following form:

[TABLE]

In Eq. (3.26) $|\Psi^{\gamma^{*}}_{T,L}\left(Q,r,z\right)|^{2}$ is the probability to find a dipole of size $r$ in a photon with the virtuality $Q$, and with transverse or longitudinal polarization. The wave functions are known (see Ref.[31] and reference therein) and they are equal to the following expressions:

[TABLE]

where $\epsilon^{2}\,=\,m^{2}_{f}+z(1-z)Q^{2}$.

Finally, the physical observables take the form:

[TABLE]

[TABLE]

3.2 Modification to also include DIS

First, we need to include the mild violation of the geometric scaling behavior of the scattering amplitude. We use the same procedure as has been suggested in Res. [57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 77, 76]: we change $\bar{\gamma}$ in Eq. (2.10)

[TABLE]

where $\kappa\,=\,\chi^{\prime\prime}\left(\gamma_{cr}\right)/\chi^{\prime}\left(\gamma_{cr}\right)=9.8$. $\chi\left(\gamma\right)$ is the BFKL kernel which has the foliowing form

[TABLE]

where $\psi(z)=d\Gamma(z)/dz$ is the Euler $\psi$-function (see Ref.[54] formula 8.360).

Since we take into account the contribution of the heavy $c$-quark we introduce a correction due to large mass of this quark:

[TABLE]

In describing the saturation phenomena and fitting the strong interaction data, we assumed that the QCD coupling is frozen at some value of momentum $\mu_{\rm soft}$. However, for DIS we take into account the running QCD coupling, replacing Eq. (3.31) by the following expression

[TABLE]

where $\mu$ denotes the typical mass of the soft strong interaction $\mu\sim 1\,GeV$ and

[TABLE]

with $\beta=3/4$.

We consider the strong interaction data for energies $W\geq 0.546\,TeV$, while the experimental data from HERA were measured for lower energies. Therefore, we need to include the contribution of the secondary Reggeons which give a substantial contribution[84].

[TABLE]

with $Q_{0}=1\,GeV$.

The final equation for $F_{2}$ takes a form:

[TABLE]

3.3 The description of the HERA data

We introduce in Eq. (3.38), a set of new parameters for DIS: $m_{q}$-mass of the light quark, which we hope will be of the order of the constituent quark mass ($\sim 300\,MeV$), the mass of charm quark ($m_{c}=1.2\div 1.5\,GeV$), $\mu$ which we believe will be of the order of 1 GeV, and we introduce two new parameters $A$ and $\alpha_{I\!\!R}\left(0\right)$ for the secondary Reggeon contribution. For $A$ there is only one restriction that at $x_{bj}=4\,10^{-6}$ $\sigma^{\rm light\,q}_{{I\!\!R}}\leq 0.02\sigma_{tot}$, while $\alpha_{I\!\!R}\left(0\right)=0.4\div 0.6$.

In Table 3 we display the parameters that were determined by fitting to the data, Fig. 2 shows the quality of our fit to the DIS HERA data.

We consider the fit shown in Fig. 2 to be in very good agreement with the experimental data, and to demonstrate that our model is able to describe the hard processes to within an accuracy of 5%.

4 Inclusive production

The cross section of the inclusive production is a very important observable for our estimates, since it indicates how well, we can describe the multi particle generation processes in our model. We have described the experimental data in our soft interaction model [26], we now recalculate using our generalization of the model, that we have discussed above. Ref.[65] showed that the CGC/saturation approach is able to describe the LHC data on inclusive production. In this section we re-visit these calculations, using our model, which we can now apply both to soft and to hard processes.

The expression for the inclusive cross section takes the form[86, 65, 31] (see Fig. 3 for notation)

[TABLE]

where the scattering amplitudes $N^{h_{i}}_{G}$ can be found from the dipole amplitude[86]

[TABLE]

and $r$ denotes the dipole size. $C_{F}=(N^{2}_{c}-1)/2N_{c}$. For further discussion it is convenient to introduce two more observables

[TABLE]

where

[TABLE]

Taking for $N\left(y_{i};r;b\right)$ in Eq. (4.40) the amplitude of Eq. (3.27) we obtain

[TABLE]

The values of $\sigma_{\rm NSD}\,=\,\sigma_{\rm tot}-\sigma_{\rm el}-\sigma_{\rm single\,diffraction}$ we take from the description of total and diffraction cross section in our model[25]. One can see that integral over $p_{T}$ is logarithmical divergent at small $p_{T}$. As shown in Ref.[87] this divergence is regularized by the mass of produced gluon jet at $y=0$. In Fig. 4 we plotted our estimates for \frac{dN}{dy}\Big{|}_{y=0} using the value of this mass as was taken in Ref.[65] $m_{\rm jet}=350\,MeV$. The agreement with the experimental data is good and it gives us confidence that our model is able to discuss the typical process of many particle production. We do not need to discuss the rapidity distribution of the single inclusive cross section, since it has been discussed in Refs.[26, 65], where it is shown that this distribution agrees with the experimental data.

5 Azimuthal angle correlations

5.1 Double inclusive cross section

The Mueller diagram[85] for double inclusive cross section is shown in Fig. 3-b. Using Eq. (4) this cross section can be written in the form:

[TABLE]

We can re-write Eq. (5.44) in a different form if we introduce

[TABLE]

Note, that $\vec{Q}_{T}$ denotes the transverse momentum carried by the BFKL Pomeron, which emits gluons with momentum $\vec{p}_{1,T}$ or $\vec{p}_{2,T}$.

Plugging Eq. (5.46) in Eq. (5.44) the expression for the double inclusive cross section takes the form

[TABLE]

Therefore, using either Eq. (5.44) or Eq. (5.1) and the decomposition of Eq. (5.45), one can calculate the double inclusive cross section.

5.2 Bose-Einstein correlation: energy dependence

The double inclusive cross section of two identical gluons has the following general form:

[TABLE]

where $C\left(L_{c}|\vec{p}_{T2}-\vec{p}_{T1}|\right)$ denotes the correlation function, and $L_{c}$ the correlation length. The first term in Eq. (5.48) is given by Eq. (5.44) or Eq. (5.1), while the second term describes the interference diagram for the identical gluons (see Fig. 3-c and Refs. [17, 32] for details). The expression for the interference term is more transparent in momentum representation, where it has the form

[TABLE]

Eq. (5.2) takes into account that the lower BFKL Pomerons in Fig. 3-c, carry momenta $\vec{Q}_{T}-\vec{p}_{12,T}$, where $\vec{p}_{12,T}\,\equiv\,\vec{p}_{1,T}\,-\,\vec{p}_{2,T}$.

Eq. (5.2) can be re-written in the impact parameter representation using Eq. (4)

[TABLE]

Finally, using the decomposition of Eq. (5.45), we can calculate the correlation function. In Fig. 5 we show the calculated correlation function

[TABLE]

From this figure we note that the correlation function does not depend on energy. This is an expected result. Indeed, the production of two parton showers, which is taken into account in Fig. 3-b and Fig. 3-c, leads to the correlation function, which does not depend on $y_{12}=|y_{1}-y_{2}|$ (long range rapidity correlations(LRC)). This happens in our approach where the structure of one parton shower cannot be reduced to the exchange of the one BFKL Pomeron. Fig. 5 illustrates that the dependence on energy also cancels in the ratio of Eq. (5.51).

5.3 Bose-Einstein correlation: values of $v_{n}$ and its

multiplicity dependence

We first introduce $v_{n}$, that can be defined it terms of the following representation of the double inclusive cross section

[TABLE]

where $\varphi$ is the angle between $\vec{p}_{T1}$ and $\vec{p}_{T2}$. $v_{n}$ is determined from $v_{n,n}\left(p_{T1},p_{T2}\right)$

[TABLE]

Eq. (5.53)-1 and Eq. (5.53)-2 depict two methods of how the values of $v_{n}$ have been extracted from the experimentally measured $v_{n,n}\left(p_{T1},p_{T2}\right)$. Where $p^{\rm Ref}_{T}$ denotes the momentum of the reference trigger. These two definitions are equivalent if $v_{n,n}\left(p_{T1},p_{T2}\right)$ can be factorized as $v_{n,n}\left(p_{T1},p_{T2}\right)\,=\,v_{n}\left(p_{T1}\right)\,v_{n}\left(p_{T2}\right)$. In this paper we use Eq. (5.53)-1 definition.

Taking into account Eq. (5.51) and Eq. (5.52) we obtain

[TABLE]

Eq. (5.54) gives the prescription for the calculation of $v_{n}$ that is measured as a sum of the events with all possible multiplicities of the secondary hadrons. However, in practice, only events with multiplicities larger than $2\bar{n}$, where $\bar{n}$ is the average multiplicity which are measured in single inclusive experiments. Fig. 6 shows our calculations for W= 13 TeV.

The dependence of $v_{n}$ on the multiplicity of the event has been discussed in Ref.[32]. Eq. (5.50) takes a different form:

[TABLE]

where $\Omega\left(r=1/mY,b\right)$ is given by Eq. (2.21). Eq. (5.55) describes the Bose-Einstein correlations in the event whose multiplicity is large than $5\bar{n}$ ($N\geq 5\bar{n}$), where $\bar{n}$ is the average multiplicity.

Comparing Fig. 6 and Fig. 7 one can see that in the framework of our approach, $v_{n}$ do not depend on the multiplicity of the event. This independence is in excellent agreement with the experimental data (see Ref.[9] and Fig. 8). Note, that $v_{n}$ do not depend on $N$ only for proton-proton scattering, while for hadron-nucleus collisions, such dependence is considerable.

Fig. 5 shows that the correlation length $L_{c}\approx 1/m_{1}$ (the typical momentum is about $m_{1}$). From Table 1, the technical reason for this is clear, since the component with such characteristic momentum makes the largest contribution. In more general language the correlation length depends on the non-perturbative hadron structure. In terms of the processes, this typical transverse momentum is responsible for diffractive scattering with the production of hadrons with small masses. Intuitively, we expect that diffractive production of large masses, which depend on the saturation scale, can lead to larger typical momenta (smaller correlation length). We will discuss these processes in the next section.

5.4 Bose-Einstein correlation: contribution of the semi-enhanced

and enhanced diagrams (diffraction production of large masses)

In Fig. 9 we show the diagrams in our model that have not been taken into account. They correspond to single diffraction in the region of large masses(Fig. 9-a and Fig. 9-b), and to double diffraction in two bunches of particles with large masses (Fig. 9-c and Fig. 9-d).

One can see from Fig. 9, that all these diagrams contain the integration over $y^{\prime}$ . This integration is concentrated in the region $Y-y^{\prime}\propto 1/\Delta_{\rm BFKL}$, where $\Delta_{\rm BFKL}$ is the intercept of the BFKL Pomeron. Preforming this integration, we reduce the diagrams of the upper part of Fig. 9 to almost the same expression as it was used in the previous section, but instead of $g_{i}\left(b\right)$ we need to insert the $b$- dependence of the triple Pomeron vertex, which in our model has the following form:

[TABLE]

Bearing this in mind we can re-write $\frac{d\sigma}{dy\,d^{2}p_{T}\,d^{2}B\,d^{2}b}$ of Eq. (4) in the form

[TABLE]

In Eq. (5.58) we restrict ourselves, by accounting only for interaction with the state $|1>$, as $g_{1}(0)\,\gg\,g_{2}(0)$. Since in our model we have $m\,\gg\,m_{1}$, we can put $b=0$. and reduce Eq. (5.58) to

[TABLE]

for diagrams of Fig. 9-a and Fig. 9-b.

For the diagrams of Fig. 9 -c and Fig. 9-d which correspond to double diffraction in large masses, we obtain

[TABLE]

Plugging Eq. (5.59) and Eq. (5.4) into Eq. (5.45), we can calculate the double inclusive cross sections. Using them and plugging in Eq. (5.50) and Eq. (5.51), we obtain the correlation function and $v_{n}$. The result of these calculations is shown in Fig. 10. One can see that contributions of semi-enhanced and enhanced diagrams, which are closely related to the processes of diffractive production of large masses in single diffraction (LMD-SD) and in double diffraction (LSD-DD), increase the typical transverse momentum in $v_{n}$ dependence, on transverse momenta.

Such behaviour is a direct consequence of the fact that typical momenta in the LSD contribution are of the order of $Q_{s}$, which is larger than $m_{1}$ and $m_{2}$, which determine the hadron structure (see Fig. 11)

Comparing Fig. 10-a and Fig. 10-b, shows that the typical momentum for the sum of the diagrams, is larger than for the non-enhanced diagrams. Fig. 11 displays the dependence of $v_{n,n}$ in the semi-enhanced and enhanced diagrams. Comparing this figure with Fig. 10-b, we note that the contribution of these diagrams are larger than the non-enhanced one, leading to the explanation of $p_{T}$ dependence in the experimental data of Fig. 13. Therefore, in our model the typical momentum is close to $Q_{s}$.

5.5 Comparison with the experiment

. In Fig. 13 and Fig. 13 we plot the experimental data [9] and the results of our calculations. One can see that we predict the values and $p_{T}$ dependence of $v_{n}$ which are in agreement with the experimental data. We wish to stress that we used Eq. (5.52)-1 for the estimates of the values of $v_{n}$, but one can see that our prediction for $v_{n,n}$ are also in accord with the data. As we have mentioned the semi-enhanced and enhanced diagrams are closely related to the processes of large mass diffraction. On the other hand, these processes give only about 30% contributions (see Table 2). Indeed, at $W=13\,TeV$ $R^{\rm lmd}_{\rm sd}=\sigma^{\rm lmd}_{\rm sd}/(\sigma_{\rm el}+\sigma^{\rm smd}_{\rm sd}+\sigma^{\rm smd}_{\rm dd})\,=\,0.26$ and $R^{\rm lmd}_{\rm dd}=\sigma^{\rm lmd}_{\rm sd}/(\sigma_{\rm el}+\sigma^{\rm smd}_{\rm sd}+\sigma^{\rm smd}_{\rm dd})\,=\,0.16$.

Such essential difference stems from the fact that the cross sections of diffractive production should be multiplied by the survival probability factor $\exp\left(-2\Omega(r,Y-Y_{0},b\right)$ (see Eq. (2.21) and Ref.[25]). This factor results in substantial suppression of the diffractive production, however, it is absent in the double inclusive cross sections.

6 Conclusions

In this paper we generalized our model to include the hard processes and presented our estimates for $v_{n}$ for proton-proton collisions at high energy. Our main result can be formulated shortly: the model predicts Bose-Einstein correlations which lead to the values of $v_{n}$, that are in accord with the experimental values. Our estimates are obtained from a model which is able to describe the typical soft observables for diffractive production, such as total and elastic cross section and cross section of diffraction production, inclusive cross sections, long range rapidity correlations and the deep inelastic $F_{2}$ structure function. In spite of being a phenomenological model which parameterizes the data rather than gives a theoretical interpretation, we believe that our model leads to reliable predictions for $v_{n}$ at high energies. This belief is based not only on the fact that the model describes both diffractive processes and processes of the multi-particle generation, but also on the fact that it includes all that we know from CGC on the behavior of the scattering amplitude in the saturation region. We showed that the angular correlations do not depend on energy and multiplicity, in accord with the experimental data.

Therefore, before making extreme assumptions on proton-proton collisions, such as the production of quark-gluon plasma in the large multiplicity events, we need to explain what happens to the Bose-Einstein correlations which are so large, that they are able to describe the angular correlations in the proton-proton scattering, without taking into account interactions in the final state.

Acknowledgements We thank our colleagues at Tel Aviv University and UTFSM for encouraging discussions. Our special thanks go to Carlos Cantreras, Alex Kovner and Michel Lublinsky for elucidating discussions on the subject of this paper.

This research was supported by the BSF grant 2012124, by Proyecto Basal FB 0821(Chile) , Fondecyt (Chile) grant 1140842, and by CONICYT grant PIA ACT1406.