Resonance production in Pomeron-Pomeron collisions at the LHC

A Regge pole model for Pomeron-Pomeron total cross section in the resonance region $\sqrt{M^2}\le$ 5 GeV is presented. The cross section is saturated by direct-channel contributions from the Pomeron as well as from two different $f$ trajectories, accompanied by the isolated f$_0(500)$ resonance which dominates the $\sqrt{M^{2}}\lesssim 1$ GeV region. A slowly varying background is taken into account. The calculated Pomeron-Pomeron total cross section cannot be measured directly, but is an essential part of central diffractive processes. In preparation of future calculations of central resonance production at the hadron level, and corresponding measurements at the LHC, we normalize the Pomeron-Pomeron cross section at large masses $\sigma_{t}^{PP} (\sqrt{M^2}\rightarrow \infty) \approx$ 1 mb as suggested by QCD-motivated estimates.


I. INTRODUCTION
Central production in proton-proton collisions has been studied from the low energy range √ s = 12.7-63 GeV at the ISR at CERN up to the presently highest energy of √ s = 13 TeV available in Run II at the LHC.Ongoing data analysis of central production events include data taken by the COMPASS Collaboration at the SPS [1], the CDF Collaboration at the TEVATRON [2], the STAR Collaboration at RHIC [3], and the ALICE and LHCb Collaborations at the LHC [4,5].A comprehensive survey of central exclusive production is given in a recent review article [6].
The analysis of central production necessitates the simulation of such events to study the acceptance and efficiency of the complex large detector systems.With the existing detector upgrade programmes for central production measurements at RHIC and at the LHC, much larger data samples are expected in the next few years which will allow the analysis of differential distributions.The purpose of the study presented here is the development of a Regge pole model for simulating such differential distributions.
The study of central production, in particular at the soft scale, is interesting for a variety of reasons.Here, we refer to central production as arising from the fusion of two stronlgy interacting colour singlet objects, and do not discuss any contributions due to photon exchange.The absence of a hard scale precludes a perturbative QCD description.The traditional framework for studying soft hadronic processes has been the Regge formalism.In this formalism, bound states are associated to Regge trajectories.The classification of mesons by means of nonlinear Regge trajectories has spectroscopic value by its own.At high energies, the hadronic interaction is dominated by the exchange of a leading trajectory, the pomeron.Within QCD, it is conjectured that this trajectory represents the exchange of purely gluonic objects.The study of central production at high energies allows to identify the contribution from the Pomeron trajectory.The dynamics of the corresponding multi-gluon colour-singlet exchange is presently only poorly understood within QCD, and such studies will hence contribute to an improved QCD-based understanding of Regge phenomenology.The fusion of multi-gluon objects is characterized by a gluon-dominated environment with highly suppressed quark degrees of freedom, and the evolution of this initial state is expected to populate with increased probability gluon-rich hadronic states, glueballs, and hybrids.The analysis of these centrally produced resonances by a Partial Wave Analysis reveals the quantum numbers J P C of these resonances.Of particular interest is the search for states with exotic quantum numbers which cannot be q q-mesons, and hence must be exotic such as of tetra-quark nature (q q + qq), or gluonic hybrid (q q + gluon).Moreover, the decomposition into states of known quantum numbers will shed new light also on the existence of numerous states in the scalar sector, a topic of fundamental interest in hadron spectroscopy [7].This article is organized as follows.In the introduction in Sect. 1, the study of central production at hadron colliders is motivated.In Sect.2, central production is reviewed.In Sect.3, the dual resonance model of Pomeron-Pomeron scattering is analysed.Non-linear complex meson trajectories are introduced in Sect. 4. Two f trajectories, relevant for the calculation of the Pomeron-Pomeron cross section, are discussed in Sect.5, while in Sect.6 the Pomeron trajectory is presented.In Sect.7, the f 0 (500) resonance is examined.The Pomeron-Pomeron total cross section is investigated in Sect.8.A summary and an outlook for more detailed studies of the topic presented here arXiv:1512.04977v2[hep-ph] 31 Jan 2016 In Fig. 1, the differential cross section dσ / dP t dM is shown for exclusive pion pair production in the CDF Run II at the TEVATRON.Clearly seen in this picture is the resonance structure associated to the f 2 (1270).The complete kinematical determination of the final state of centrally produced pion or kaon pairs requires the measurement of the 3-momentum (p x ,p y ,p z ) of both the positive and the negative partner of the pair.The experimental single track acceptance is, however, limited by finite detector coverage in pseudorapidity, as well as by a cut-off in minimum transverse momentum p T .This single track acceptance translates into missing acceptance for pairs of low mass and low transverse momentum.This missing pair acceptance is visible in Fig. 1 in the data analyzed by the CDF Collaboration.For pairs of low masses, only the high end tail of the transverse momentum distribution can be measured.The extrapolation of the cross section to low transverse momenta is, however, possible based on models which are able to reproduce the resonance structures in the part of phase space covered by the detector acceptance.[9], CDF Collaboration in the middle [10], and ALICE Collaboration on the right [11].
The distributions of pion-pair invariant masses measured in proton-proton collisions by the COMPASS, the CDF and the ALICE Collaboration are shown in Fig. 2 on the left, in the middle and on the right, respectively.In all these measurements, clear resonance structures are seen.At the energy of the COMPASS measurement √ s = 18.9 GeV, Reggeons still contribute significantly to central production as evidenced in the prominent ρ-peak.In addition, peaks associated to the f 0 (980) and the f 2 (1270) are seen, with a broad continuum extending to the two-pion threshold.At the higher TEVATRON energy of √ s = 1.96TeV, the distribution at pair masses M < 900 MeV/c 2 is significantly affected by the p T dependence of the acceptance as shown in Fig. 1.Full acceptance down to p T = 0 is reached for masses M > 900 MeV/c 2 .A clear resonance structure consistent with the f 2 (1270) is seen at around 1270 MeV/c 2 .A similar distribution is measured by the ALICE Collaboration at the LHC energy √ s = 7 TeV as shown in Fig. 2 on the right.In addition, the ALICE double gap measurement is compared to the pion-pair invariant mass distribution from no-gap events, i.e. from inclusive production.In inclusive production, the ρ as well as the K 0 s -signal are seen.These two signals are absent in the double-gap events, corroborating Pomeron-Pomeron dominance at the LHC energies discussed below.The double-gap topology of central production and the relevant kinematics are shown in Fig. 3.This figure shows central production when the incoming protons remain in the ground state on the left, when one of the protons gets diffractively excited in the middle, and when both protons get excited on the right.All these reactions proceed by the exchange of Regge trajectories α(t 1 ) and α(t 2 ) which collide in the central region to produce a system of mass M x .The total energy s of the reaction is shared by the subenergies s 1 and s 2 associated to the trajectories α(t 1 ) and α(t 2 ), respectively.The LHC energies of √ s = 7, 8 and 13 TeV are sufficient to provide Pomeron dominance and allow for the neglect of Reggeon exchange which was not the case at the energies of previous accelerators.The scope of the study presented here is the central part of the diagrams shown in Fig. 3, i.e.Pomeron-Pomeron scattering producing mesonic states of mass M x .We isolate the Pomeron-Pomeron-meson vertex shown in Fig. 4, and calculate the Pomeron-Pomeron total cross section as a function of the centrally produced system of mass M x .The emphasis in this study is the behaviour in the low mass resonance region where perturbative QCD approaches are not applicable.Instead, similar to [12,13], we use the pole decomposition of a dual amplitude with relevant direct-channel trajectories α(M 2 ) for fixed values of Pomeron virtualities, t 1 = t 2 = const.Due to Regge factorization, the calculated Pomeron-Pomeron cross section will enter the measurable proton-proton cross section [14].
The nature of the Pomeron exchange is of fundamental interest for QCD-based studies of exchange amplitudes.An effective vectorial-exchange is very successful in reproducing the energy dependence of hadron-hadron cross sections [15].Such an approach results in opposite signs for proton-proton and proton-antiproton amplitudes.Pomeron exchange, however, must yield the same sign for these two reaction channels.Recent studies on soft high-energy scattering solve this problem in terms of effective propagators and vertices for the Pomeron exchange [16].Within this model, the Pomeron exchange is decribed as an effective rank-two tensor exchange [17,18].

III. DUAL RESONANCE MODEL OF POMERON-POMERON SCATTERING
The study of Pomeron-Pomeron (P P ) scattering is related to photon-photon scattering, the main difference being the positive and negative C-parity of the Pomeron and photon, respectively.High-virtuality γ * γ * scattering is a favourite process in the framework of perturbative QCD, where the total cross section was calculated in Ref. [19].In the leading-order BFKL where is the photon virtuality and Y = ln s s0 .The quantity σ 0 is a free parameter and the exponent λ ≡ α BF KL P is the familiar BFKL eigenvalue N c α s 4 ln 2/π.We recall that the transition from photon-photon to central Pomeron-Pomeron scattering is accompanied by the change of variables Q 2 → t and s → M 2 .The above result was improved in Refs.[20,21].
Most of the studies on diffraction dissociation, single, double and central, use the triple Reggeon formalism.This approach is useful in the smooth Regge region, beyond the resonance region, but is not applicable for the production of low masses which is dominated by resonances.We solve this problem by using a dual model.
The one-by-one account of single resonances is possible, but not economic for the calculation of cross section, to which a sequence of many resonances contributes at low masses.These resonances overlap and gradually disappear in the continuum at higher masses.An approach to account for many resonances, based on the idea of duality with a limited number of resonances lying on nonlinear Regge trajectories, was suggested in Ref. [22].Later on, this approach was used in Refs.[12,13] to calculate low mass single-and double-diffractive dissociation at the LHC.The main idea behind this approach is illustrated in Fig. 5, realized by dual amplitudes with Mandelstam analyticity (DAMA) [23].For s → ∞ and fixed t it is Regge-behaved.Contrary to the Veneziano model, DAMA not only allows for, but rather requires the use of nonlinear complex trajectories providing the resonance widths via the imaginary part of the trajectory.In the case of limited real part of the trajectory, a finite number of resonances is produced.More specifically, the asymptotic rise of the trajectories in DAMA is limited by the condition, in accordance with an important upper bound, For our purpose of central production, the direct-channel pole decomposition of the dual amplitude A(M 2 X , t) is relevant.Different trajectories α i (M 2 X ) contribute to this amplitude, with α i (M 2 X ) a nonlinear, complex Regge trajectory in the Pomeron-Pomeron system, The pole decomposition of the dual amplitude A(M 2 X , t) is shown in Eq. ( 3), with t the squared momentum transfer in the P P → P P reaction.The index i sums over the trajectories which contribute to the amplitude.Within each trajectory, the second sum extends over the bound states of spin J.The prefactor a in Eq. (3) has the numerical value a = 1 GeV −2 = 0.389 mb.
The pole residue f (t) appearing in the P P → P P system is fixed by the dual model, in particular by the compatibility of its Regge asymptotics with Bjorken scaling and reads where t 0 is a parameter to be fitted to the data.However, due to the absence of data so far, we set t 0 = 0.71 GeV 2 for the moment as in the proton elastic form factor.Note that the residue enters with a power (J +2) in Eq. (3), thereby strongly damping higher spin resonance contributions.
The imaginary part of the amplitude A(M 2 X , t) given in Eq. ( 3) is defined by ( For the P P total cross section we use the norm and recall that the amplitude A and the cross section σ t carry dimensions of mb due to the dimensional parameter a discussed above.The Pomeron-Pomeron channel, P P → M 2 X , couples to the Pomeron and f channels dictated by conservation of quantum numbers.In order to calculate the P P cross section, we therefore take into account the trajectories associated to the f 0 (980) and to the f 2 (1270) resonance, and the Pomeron trajectory.

IV. NONLINEAR, COMPLEX MESON REGGE TRAJECTORIES
A non-trivial task for analytic models of Regge trajectories consists in deriving the imaginary part of the trajectory from the seemingly linearly increasing real part [24][25][26].The importance of the nonlinearity of the real part was studied in Refs.[27][28][29].A dispersion relation connects the real and imaginary part of the trajectory.We follow Ref. [25] to relate the nearly linear real part of the meson trajectory to its imaginary part, In Eq. ( 7), PV denotes the Cauchy Principal Value of the integral.The imaginary part is related to the decay width by The quantity α in Eq. ( 8) denotes the derivative of the real part, α = d e α(s) ds .The relation between Γ(M ) and m α(s) requires m α(s) > 0. In a simple analytical model, the imaginary part is chosen as a sum of single threshold terms [25] The imaginary part of the trajectory shown in Eq. ( 9) has the correct threshold and asymptotic behaviour.Since m α(s) > 0, all the expansion coefficients c n must be positive.The values of s n represent kinematical thresholds of decay channels.The highest threshold, higher than all the resonance masses lying on the trajectory, is chosen as an effective threshold.This highest threshold ensures that e α(s) tends to a constant value for s → ∞.
The parameterisation of the real and imaginary part of a meson trajectory, and the extraction of the expansion coefficients c n shown in Eq. ( 9), are derived in Appendix A for the case of the ρ-a trajectory.
V. TWO f TRAJECTORIES Apart from the Pomeron trajectory discussed below, the direct-channel f trajectory is essential in the PP system.Guided by conservation of quantum numbers, we include two f trajectories, labelled f 1 and f 2 , with mesons lying on these trajectories as specified in Table I.

VI. THE POMERON TRAJECTORY
While ordinary meson trajectories can be fitted both in the resonance and scattering region corresponding to positive and negative values of the argument, the parameters of the Pomeron trajectory can only be determined in the scattering region M 2 < 0. The poles of this trajectory at M 2 > 0 are identified with glueball candidates.An extensive literature on such candidates exists, including theoretical predictions and experimental identification.The status of glueballs is, however, controversial and a topic of ongoing discussions and debate; see Refs.[7,30] and references therein.In this study, we associate the bound states of the Pomeron trajectory to glueball candidates, as has been done previously in Refs.[29,[31][32][33].
A comprehensive fit to high-energy pp and pp of the nonlinear Pomeron trajectory is discussed in Ref. [14] α with ε = 0.08, α = 0.25 GeV −2 and s 0 the two pion threshold s 0 = 4m 2 π .The value of c is taken as c = α /10 = 0.025.
In order to be consistent with the mesonic trajectories shown above, the linear term in Eq. ( 10) is replaced by a heavy threshold mimicking linear behaviour in the mass region of interest (M < 5 GeV), where M H is an effective heavy threshold set at M = 3.5 GeV.The coefficients α 0 , α 1 and α 2 are chosen such that the Pomeron trajectory of Eq. ( 11) has a low energy behaviour as defined by Eq. (10).The experimental data on central exclusive pion-pair production measured at the energies of the ISR, RHIC, TEVATRON and the LHC collider all show a broad continuum for pair masses m π + π − < 1 GeV/c 2 .This mass region is experimentally difficult to access due to the missing acceptance for pion-pairs of low mass and low transverse momentum p T as discussed above.The population of this mass region is attributed to the f 0 (500), a resonance which has been controversial for many decades.In the 2010 edition of the Review of Particle Physics (RPP), this resonance is listed as f 0 (600) with a mass M 0 in the range 400 < M 0 < 1200 MeV, and a width Γ in the range 600 < Γ < 1000 MeV.Since the RPP edition of 2012, this resonance is listed as f 0 (500) with mass in the range 400 < M 0 < 550 MeV, and width in the range 400 < Γ < 700 MeV [34].
The f 0 (500) resonance is of prime importance for the understanding of the attractive part of the nucleon-nucleon interaction, as well as for the mechanism of spontaneous breaking of chiral symmetry.The nature of the f 0 (500) is a topic of ongoing studies and discussions, it is however generally agreed that it cannot be interpreted as a predominant q q-state.The non-ordinary nature of the f 0 (500) resonance is corroborated by the fact that it does not fit into the Regge description of classifying q q-states into trajectories [35].A possible interpretation of the f 0 (500) is a tetraquark configuration consisting of two valence and two antiquarks in colour-neutral state.It was shown that such a configuration can give rise to a nonet of light scalar-isoscalar mesons [36].Different approaches interpret the f 0 (500) as arising from an inner tetra-quark structure and changing to an outer structure of a pion-pion state [37].There is strong evidence that this f 0 (500) state belongs to a SU(3) nonet composed of the f 0 (500), f 0 (980), a 0 (980) and K * 0 (800).
In spite of the complexity of the f 0 (500) resonance, and the controversy on its interpretation and description, we take here the practical but simple-minded approach of a Breit-Wigner resonance [38] In Eq. ( 12), the parameterisation of the relativistic Breit-Wigner amplitude is shown with M 0 and Γ the mass and width, respectively.Here, the prefactor a is added for consistency with the definition of the amplitude shown in Eq. ( 3).The Breit-Wigner amplitude of Eq. ( 12) is used below for calculating the contribution of the f 0 (500) resonance to the Pomeron-Pomeron cross section.

VIII. POMERON-POMERON TOTAL CROSS SECTION
The Pomeron-Pomeron cross section is calculated from the imaginary part of the amplitude by use of the optical theorem In Eq. ( 13), the index i sums over the trajectories which contribute to the cross section, in our case the f 1 , f 2 and the Pomeron trajectory discussed above.Within each trajectory, the summation extends over the bound states of spin J as expressed by the second summation sign.The value f i (0) = f i (t) t=0 is not known a priori.The analysis of relative strengths of the states of trajectory i will, however, allow to extract a numerical value for f i (0) from the experimental data .
In addition to the contributions discussed above, a background term is added to the PP cross section.This background is of form [39] σ P P backgr.(M 2 ) = c * (0.1 + log(M 2 )) mb, (15) with the numerical value of the parameter c fitted to data.In Fig. 9, the different contributions to the PP total cross section are shown.The contribution of the f 0 (500) resonance according to Eq. ( 14) is displayed by the dashed cyan line.The contribution of the f 1 trajectory indicated by the solid green line clearly shows the f 0 (980) and the f 1 (1420) resonance.The higher mass states, the f 2 (1810) and the f 4 (2300), are barely visible due to their reduced cross section and much larger width.Similarly, the contribution of the f 2 trajectory indicated by the dashed blue line shows peaks for the f 2 (1270) and the f 4 (2050) resonances, with the f 6 (2510) barely visible.The contributions of both f 1 and f 2 trajectory show a kink at about M = 5.5 GeV due to the heavy threshold s 2 = 30 GeV 2 .The contribution from the Pomeron trajectory is displayed in Fig. 9 by the dashed magenta line.Visible here is the resonance structure due to the J = 2, 4 and 6 states on the trajectory labelled by gb(J=2), gb(J=4) and gb(J=6), respectively.Beyond the heavy threshold, M = 3.5 GeV, the transition to the continuum is seen, reflecting the behaviour of the real and imaginary part of the trajectory as shown in Fig. 8.The background contribution to the PP cross section is shown in Fig. 9 by the dashed black line, and is normalized here to represent approximately 10% of the signal at M = 7 GeV.
The Pomeron-Pomeron total cross section is calculated by summing over the contributions discussed above, and is shown in Fig. 9 by the solid black line.The prominent structures seen in the total cross section are labelled by the resonances generating the peaks.The model presented here does not specify the relative strength of the different contributions shown in Fig. 9.A Partial Wave Analysis of experimental data on central production events will be able to extract the quantum numbers of these resonances, and will hence allow to associate each resonance to its trajectory.The relative strengths of the contributing trajectories need to be determined from the experimental data.

IX. SUMMARY AND OUTLOOK
A Regge pole model is presented for calculating the Pomeron-Pomeron total cross section in the resonance region √ M 2 ≤ 5 GeV.The direct-channel contributions of the Pomeron and two f trajectories, including a background, are presented.The resonance region √ M 2 ≤ 1 GeV is characterised by a Breit-Wigner parameterisation of the f 0 (500) resonance.The relative strength of these contributions cannot be specified within the model, and must hence be determined from the analysis of experimental data.The model presented allows an extension to central production of strangeonia and charmonia states by taking into account the direct-channel contribution of the respective trajectories.Moreover, this model can be extended to lower beam energies where not only Pomeron-Pomeron, but also Pomeron-Reggeon and Reggeon-Reggeon diagrams need to be considered.The result of the presented work is only the starting point for a comprehensive study of central exclusive production.To make measurable predictions for the LHC, all the diagrams shown in Fig. 3 must be calculated.The results presented here are necessary and essential input for such calculations.Anticipating further studies, we remind of the possible reference points that can be used as a guide.For the absolute value we use the asymptotic value σ t ≈ 1 mb, compatible with both QCD-inspired and phenomenological estimates [40,41].The Pomeron-Pomeron total cross section depends also on Pomerons' virtualities, t 1 , t 2 .We ignored this dependence for two reasons: First, this dependence is known at best at their high values, where perturbative QCD results, such as that of Eq. ( 1), may be valid.Second, for simplicity, we fix this dependence, including it as part of the normalisation factor.Varying the t dependence and the partition between t 1 and t 2 may be attempted to account for by, following Eq.( 1), simply dividing Eq. ( 3) by √ t 1 t 2 , this may be true only for high values of t i , beyond diffraction.

X. ACKNOWLEDGEMENTS Appendix A
Although the ρ-a trajectory does not couple to the P P channel, and hence does not contribute to the P P total cross section, we list it here to present the procedure for extracting the coefficients c n in the expansion of the trajectory.For convenience, we repeat here the Ansatz for the imaginary part of the trajectory expressed in Eq. ( 9) The real part e α(s) of the trajectory expressed by Eq. ( 7) can be calculated as The equations ( 16) -( 18) are used to calculate the parameters c n of e α, e α and m α by a χ 2 -fit procedure.A linear fit to the real part provides start values for α(0) and e α , from which start values for m α(M 2 R ) are calculated.
The fit of the parameters c n is done in three steps.First, the c n are fitted to the expression of the imaginary part m α(s) as given in Eq. (16).Second, the coefficients α(0) and c 2 are extracted by using the parameterisation of the real part e α(s) as defined by Eq. (17).Third, new values are calculated for e α(s n ), and the imaginary part m α(s) is updated according to Eq. ( 16).These three steps can be repeated if necessary until convergence of the values α(0) and c n is reached.

FIG. 2 :
FIG.2: Invariant pion pair masses from COMPASS Collaboration on the left[9], CDF Collaboration in the middle[10], and ALICE Collaboration on the right[11].

FIG. 5 :
FIG. 5: Connection, through unitarity (generalized optical theorem) and Veneziano-duality, between the Pomeron-Pomeron cross section and the sum of direct-channel resonances.

FIG. 6 :
FIG. 6: Real part of f1 trajectory on the left, width function Γ(M 2 ) on the right.

FIG. 7 :
FIG. 7: Real part of f2 trajectory on the left, width function Γ(M 2 ) on the right.

FIG. 8 :
FIG. 8: Real part of pomeron trajectory on the left, imaginary part on the right.

FIG. 9 :
FIG. 9: Contributions of the f0(500) resonance, the f1, f2 and the Pomeron trajectory, and of the background to PP total cross section.

TABLE I :
Parameters of resonances belonging to the f1 and f2 trajectories.

TABLE II :
Parameters of ρ-and a-resonances.