Measurement of the energy density as a function of pseudorapidity in proton-proton collisions at $\sqrt{s}=$ 13 TeV

A measurement of the energy density in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=$ 13 TeV is presented. The data have been recorded with the CMS experiment at the LHC during low luminosity operations in 2015. The energy density is studied as a function of pseudorapidity in the ranges $-$6.6$<\eta<-$5.2 and 3.15$<|\eta|<$5.20. The results are compared with the predictions of several models. All the models considered suggest a different shape of the pseudorapidity dependence compared to that observed in the data. A comparison with LHC proton-proton collision data at $\sqrt{s}=$ 0.9 and 7 TeV confirms the compatibility of the data with the hypothesis of limiting fragmentation.


Introduction
In the framework of quantum chromodynamics (QCD), inelastic proton-proton collisions are described by a combination of hard and soft exchanges between the constituents of the protons. Hard collisions between one or multiple pairs of partons are complemented by soft parton scattering from Multiple Parton Interactions (MPI) [1][2][3][4], parton shower effects including initialand final-state radiation, which, along with projectile fragmentation, constitute the underlying event (cf. Ref. [5]). At the CERN LHC these effects can be studied at the highest possible centreof-mass energies covering a very large angular phase space. The measurement of the average energy per event in different pseudorapidity (η) regions probes our general understanding of QCD multiparticle production. Moreover, because of the extended calorimetric instrumentation of the CMS experiment beyond |η| > 3, covering the full range from −6.6 to +5.2 in pseudorapidity, smaller scattering angles may be accessed compared to other measurements.
In this paper, a measurement of the energy density in proton-proton collisions at the centreof-mass energy √ s = 13 TeV within the pseudorapidity ranges −6.6 < η < −5.2 and 3.15 < |η| < 5.20 is presented. This measurement extends the √ s and pseudorapidity range covered by previous results from the CMS [6], ATLAS [7], and LHCb [8] Collaborations. The average energy density per event is defined as where ∑ i E i is the summed energy measurements of all calorimeter towers i within a bin of pseudorapidity having a width ∆η, c(η) is the η-dependent conversion factor from the calorimeter measurements to a stable-particle level energy, and N evt is the number of selected events corrected for the contributions from noise and simultaneous pp collisions occurring in the same event (pileup). To investigate various aspects of MPIs in high-energy proton-proton collisions the measurement is performed for several different categories of events, each category defined by a specific event selection.
Moreover, the data collected at √ s = 13 TeV are analysed together with data collected at 0.9 and 7 TeV [6]. This is interesting since projectile fragmentation can then be studied in the regions close to the beam rapidity, y beam = acosh( √ s/2m p ), where m p is the mass of the projectile particle, i.e. a proton in the present case. At √ s = 13 TeV, y beam ≈ 9.5, while at √ s = 0.9 TeV it is just ≈6.8. Thus, the detectors of CMS, although located at fixed η, cover a very wide range in η = η − y beam when data recorded at different centre-of-mass energies are combined. The hypothesis of limiting fragmentation [9] suggests that particle production reveals longitudinal scaling, i.e. the dependence of very forward particle production on the centre-of-mass energy vanishes in the region η ≈ 0 [10]. In this paper, the hypothesis of limiting fragmentation is tested in collisions at √ s from 0.9 to 13 TeV.
Measurements of the energy density at collider energies are an important reference necessary for extrapolating to even higher centre-of-mass energies. The results reported here provide valuable input for the tuning of Monte Carlo models used to describe the highest energy hadronic interactions needed for the interpretation of cosmic ray measurements [11,12].

The CMS detector
At the heart of the CMS detector is a superconducting solenoid of 6 m internal diameter, providing a strong magnetic field of 3.8 T. The data used for this paper were taken in June 2015 during a period without magnetic field. Within the CMS magnet volume are an inner silicon pixel and strip tracker that measure charged particles in the range |η| < 2.5, a homogeneous lead tungstate crystal electromagnetic calorimeter, and a brass and scintillator hadron calorimeter. The corresponding endcap detectors instrument the pseudorapidity range up to |η| 3 with tracking and calorimetry. Forward Cherenkov calorimeters extend the coverage beyond |η| 3. Muons are measured in gas-ionization detectors embedded in the steel return yoke.
The hadron forward (HF) calorimeters cover the region 2.9 < |η| < 5.2 and consist of 2×432 readout towers, each containing a long and a short quartz fiber embedded within a steel absorber running parallel to the beam. The long fibers run the entire depth of the HF calorimeter (165 cm, or approximately 10 interaction length), while the short fibers start at a depth of 22 cm from the front of the detector. The response of each tower is determined from the sum of signal in the corresponding long and short fiber. There are 13 rings of towers in |η|, each with a size of ∆η 0.175, except for the lowest and highest |η| rings, which have a size ∆η 0.11 and ∆η 0.30, respectively. The azimuthal segmentation of all towers is 10 • , except for the one at highest |η|, which has ∆ϕ = 20 • .
The very forward angles on one side of CMS (−6.6 < η < −5.2) are covered by the CASTOR calorimeter. It has 16 azimuthal towers, each built from 14 longitudinal modules. The 2 front modules form the electromagnetic section, and the 12 rear modules form the hadronic section. The calorimeter is made of stacks of tungsten and quartz plates, read out by PMTs, in two halfcylindrical mechanical structures, and is placed around the beam pipe at a distance of −14.4 m away from the nominal interaction point. The overall longitudinal depth of both CASTOR and HF corresponds to 10 hadronic interaction lengths. The CASTOR calorimeter is only operated during periods of low LHC luminosity (L inst < 10 30 cm −2 s −1 ) since it cannot distinguish the secondaries from simultaneous pileup collisions.
The present analysis is restricted to the range of pseudorapidity covered by the HF and CAS-TOR calorimeters, excluding the two lowest |η| segments of the HF calorimeters because they are partially located in the shadow of the endcap calorimeters. This corresponds to a combined pseudorapidity range of 3.15 < |η| < 5.2 and −6.6 < η < −5.2. The analysis is performed using a data sample corresponding to an integrated luminosity of 0.06 nb −1 recorded with an average number of proton-proton interactions per bunch crossing of about 0.05.
A more detailed description of the CMS detector can be found in Ref. [13].

Monte Carlo models
In this paper, various Monte Carlo event generators are used to correct the data from detectorto stable-particle level and to compare with the experimental results.
The PYTHIA8 [14] generator is a general purpose Monte Carlo package that builds most of its predictive power upon hard-scattering matrix elements calculated in perturbative QCD and parton showering according to the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) [15][16][17][18][19] equations. The string fragmentation model [20] is used for hadronization. The free parameters of the simulations can be adjusted to describe measurements at different centre-of-mass energies, resulting in the production of different so-called tunes of the model [21].
In this analysis, PYTHIA8 (version 8.212) is used together with the CUETP8M1 [21], CUETP8S1 [21], and MONASH 2013 [22] tunes, as well as with the MBR model [23] combined with the 4C [24] and CUETP8M1 tunes. In the CUETP8M1 and CUETP8S1 tunes, which are based on the MONASH 2013 and 4C tunes, the parameters are adjusted to describe underlying event measurements from the Fermilab Tevatron and the LHC. The tunes are constructed using different parton distribution function sets (NNPDF2.3LO [25]) and CTEQ6L1 [26], respectively).
The EPOS-LHC [27] and QGSJETII.04 [28] generators are commonly used to describe extensive air showers in the atmosphere initiated by cosmic ray particles, where soft physics is of primary importance. A combination of Gribov-Regge multiple scattering [29], perturbative QCD, and string fragmentation are the cornerstones of both models. While QGSJETII.04 includes a small number of fundamental parameters, the phenomenology implemented in EPOS-LHC offers more opportunities for tuning. In EPOS-LHC a hydrodynamic, or collective, component is included in a parametrised form [27].
The events simulated with the MONASH and MBR tunes of PYTHIA8, as well as with the EPOS-LHC and QGSJETII.04 event generators, have been processed with a detailed simulation of the full CMS detector based on GEANT4 [30] and reconstructed using the same software sequence that is used for collision events. These four models are used to correct for detector effects.

Event selection, and definition of event classes
Events are selected online in an unbiased way by triggering the data acquisition system with the Beam Pick-up-Timing for the eXperiments (BPTX) devices [31]. Three different categories of inelastic collision events are defined offline: an inclusive inelastic (INEL) selection to be as inclusive as possible, a non-single-diffractive-enhanced (NSD-enhanced) selection, where single diffractive dissociation contributions are suppressed, and a single-diffractive-enhanced (SDenhanced) selection enriched in single diffractive dissociation events. These selections are achieved by requiring an energy deposit in the HF calorimeters above noise level either on at least one side (for the INEL events) or on both sides (for the NSD-enhanced events), with respect to the nominal interaction point of CMS. The SD-enhanced selection is defined by requiring activity in one of the calorimeters on exactly one side, with a veto condition being applied to the other side.
Energy deposition in the HF calorimeters is characterised by the calorimeter tower with the highest energy in the negative (positive) pseudorapidity region, E HF− (E HF+ ), considering all towers, except those belonging to the two rings closest to the endcap (i.e. at smallest |η|). The Table 1: Summary of the event selections used for the different event categories in data at the detector level and in simulations at the stable-particle level.
Class Detector level NSD-enhanced E HF+ > 5 GeV and E HF− > 5 GeV at least one stable particle with E > 5 GeV in −5.20 < η < −3.15 and 3.15 < η < 5.20 SD-enhanced E HF+ > 5 GeV and E HF− < 5 GeV or E HF+ < 5 GeV and E HF− > 5 GeV at least one stable particle with E > 5 GeV in 3.15 < |η| < 5.20 on one side, vetoing particles with E > 5 GeV on the other side Limiting fragmentation study E HF+ > 4 GeV and E HF− > 4 GeV one stable particle in −4.4 < η < −3.9 and 3.9 < η < 4.4 Table 2: Selection factors and purities for various event categories. The last two rows are derived using the formulae in the text. The event selection probability is determined from simulations, and the value quoted here is the average value from all event generators, with a maximal model dependence of 2%. The rightmost column quantifies the combined correction due to noise and pileup. All statistical uncertainties are negligible.
energy thresholds for event selection are determined from a study of events without beam and are optimised to effectively reduce the contribution from detector noise, while still allowing a high selection efficiency. In Fig. 1, the measured distributions for E HF− and E HF+ from collision data are shown together with the noise distributions obtained from data without the presence of LHC beams. This is achieved at the trigger level by requiring prescaled triggers where the two BPTX detectors are silent. In Fig. 1 simulated events are also shown. Events are selected for the INEL class if max(E HF− , E HF+ ) > E threshold , and for the NSD-enhanced class if min(E HF− , E HF+ ) > E threshold . An energy threshold of E threshold = 5 GeV is found to be optimal to suppress the noise contribution in both event classes for simulated and measured events. For the NSDenhanced events, the threshold could in principle be lowered down to about 3 GeV without increasing the noise contribution, but for consistency a unified threshold of 5 GeV is used for all event classes. The noise peaks in collision data contain significantly more events than in simulations, since the interaction probability in the data is about 5% while the simulated samples are produced with no pileup (each simulated event corresponds to one collision). The detector noise distribution as measured from empty-beam data are also overlaid as shaded areas.
In simulated events particle four-momenta are used to build sums of energies. At the stableparticle level (i.e. for particles with proper decay length >1 cm), simulated events are selected and M X and M Y are the invariant masses of the particle systems on the negative and positive side of the largest rapidity gap in the event, respectively. This particular criterion for stableparticle level is identical within a few percent with the INEL detector level selection [32].
The NSD-enhanced events are selected at the stable-particle level with a requirement of at least one stable particle (either charged or neutral) within the pseudorapidity acceptance of the HF calorimeters 3.15 < |η| < 5.2 on both sides of the interaction point.
The SD-enhanced events at the stable-particle level are defined by the presence of at least one stable particle with energy E > 5 GeV within the pseudorapidity range 3.15 < |η| < 5.2 on one side, whereas the other side must be devoid of particles with energy E > 5 GeV.
The phase space definitions for INEL, NSD-enhanced and SD-enhanced events at the detector and stable-particle level are summarised in Table 1. The last row of the table indicates the event selection needed for the limiting fragmentation study. This is chosen to be identical to that used in previously published data [6] to allow a direct comparison of the results.
The energy density is measured with the HF and CASTOR calorimeters by summing up all the energy deposits in the calorimeter towers above noise threshold. The value of the threshold was determined by measuring the detector noise and beam backgrounds using empty-beam triggers (see Fig. 1 for HF results) and is chosen to be 5 GeV in HF and 2.5 GeV in CASTOR. The energy density measurement is performed as a function of |η|. In the range 3.15 < |η| < 5.2 the corresponding measurements at positive and negative pseudorapidities in HF are averaged, while for −6.6 < η < −5.2 the energy in CASTOR is used. For the SD-enhanced measurement only the side on which the HF calorimeter is above noise level (thus, opposite to the forward rapidity gap) is used for the measurement.

Data analysis
The measurement of the energy density according to Eq. (1) requires the determination of the number of selected events N evt and the energy sum, ∑ i E i .

Event counting, noise, and pileup
The number of selected events in the analysis, N sel , is corrected to eliminate the residual contribution from detector noise to yield the corrected number of events, N corr , containing only signal and no noise events. With N ZB and N EB being the number of events collected with the unbiased and empty-beam triggers, respectively, and f ZB and f EB the corresponding fractions of offline-selected events, we can define the number of selected collision events N sel = N ZB f ZB , and the number of noise events in the same data sample N noise = N ZB f EB . The latter contains N sig+noise = N corr f EB events that are selected because towers in the same event are above threshold due to signal and noise fluctuations. Thus, the corrected number of events containing collisions is where we define the purity as p = (1 − f EB / f ZB )/(1 − f EB ). The purity of the data used in this analysis is found to be above 99%. The noise contribution depends weakly on the event selection criteria.
The reconstructed number of events is also corrected for the effect of pileup. The number of proton-proton interactions per bunch crossing n follows a Poisson distribution with a mean value λ. The probability to have no interaction is given by e −λ = 1 − N corr /N ZB , which allows λ to be determined from inelastic collisions in data. Here we find λ = − ln(1 − f ZB p) = 0.05, using the value of f ZB determined from the INEL event selection.
The number of visible collisions in N tot bunch crossings is N coll = N tot ∑ ∞ n=0 n Pois (n; λ) = N tot λ, where is the probability for each collision to be observed and selected. In the presence of pileup another important quantity is the probability for the observation of events with exactly n simultaneous collisions, n = 1 − (1 − ) n . Thus the average number of actually observed events is N obs = N tot ∑ ∞ n=0 n Pois (n; λ). Using this result we can correct for pileup using the factor For the data analysis we use the corrected number of collisions N evt = N coll = N ZB f ZB p f PU for counting inelastic events in Eq. (1). In general, the impact of pileup depends on the event selection procedure. In particular, an exclusivity criterion as used in the SD-enhanced event selection leads to fewer selected events  Table 2. The combined corrections for each event selections are at the level of 1%. The value quoted for is the average obtained from the different event generators with a maximum discrepancy between the model predictions of about 2%, which induces a maximum effect on p f PU of less than < 10 −3 .

Energy measurement
The measured response from the calorimeters is corrected to the stable-particle level to provide a well-defined event classification and energy quantification for comparisons to the model predictions. The corrections are applied explicitly for each range in pseudorapidity. There is no relevant migration or detector smearing in pseudorapidity; it is basically the characteristic response of the calorimeters, as well as the event selection acceptance and inefficiency, that is corrected. These corrections are determined with the PYTHIA8 tune MONASH 2013, PYTHIA8 tune 4C with MBR model, EPOS-LHC, and QGSJETII.04 simulated event samples. The corrections are evaluated from the ratio of the predictions at the stable-particle level to the predictions at the detector level for every |η| bin. The final correction is the average of the four different simulated samples. The magnitude of the correction varies from 1.5 to around 2.5 depending on the value of |η| and the selection criteria applied at the stable-particle level. The main contribution to the correction is related to the extrapolation of observed detector-level energy above the calorimeter noise threshold to the energy with no threshold applied at the stable-particle level.

Uncertainties
The energy scales for the HF and CASTOR calorimeters are known to within an accuracy of 10% [6] and 17% [33], respectively. These are the dominant sources of experimental uncertainty in this analysis.
The impact of the energy scale uncertainty on the measurement of the energy density is estimated by scaling the tower energies up and down by the energy scale uncertainties in the data while keeping the simulated correction factors constant. The resulting impact is 10% for HF and 17% for CASTOR as expected.
To assess the residual impact of noise on the event selection, the thresholds in the event selection at detector level are increased from 5 to 5.5 GeV for all INEL, NSD-enhanced, and SDenhanced events. This corresponds to an improved noise rejection at the expense of larger correction factors. The resulting uncertainties are about 0.7, 0.01, and 5% for INEL, NSD-enhanced, and SD-enhanced events, respectively.
Furthermore, to study the impact of the energy threshold on the energy measurement, the threshold for the tower energy sum is increased by the energy scale uncertainty, which leads to uncertainties of 1% for all three event selections categories.
The systematic uncertainty due to model dependence is estimated from the maximum variation of the correction factor values obtained using the event generators PYTHIA8 with MONASH and 4C+MBR tunes, EPOS-LHC, and QGSJETII.04. The resulting uncertainty is below 3.5% for INEL and NSD-enhanced events, while for SD-enhanced events it varies from 16 to 37%, depending on η.
The statistical uncertainty is < 1%, which is significantly smaller than the systematic uncertainties.
The individual contributions for each |η| bin are assumed to contribute quadratically to the total systematic uncertainty since the contributions are not correlated within a bin; the systematic uncertainties are, however, highly correlated between different |η| bins. All uncertainties are summarised in Table 3.

Results
The measured energy density, dE/dη, in the range −6.6 < η < −5.2 and 3.15 < |η| < 5.20, corrected to the stable-particle level, is presented in Figs. Figure 2: Energy density at the stable-particle level for the INEL (upper row), NSD-enhanced (middle row), and SD-enhanced (lower row) event selections categories compared to predictions from PYTHIA8 MONASH, EPOS-LHC, and QGSJETII.04. The gray band shows the total systematic uncertainty. The right panels show the ratio of model predictions to measured data. 3 Figure 3: Energy density at the stable-particle level for the INEL (upper row), NSD-enhanced (middle row), and SD-enhanced (lower row) event selections categories compared to predictions from PYTHIA8 with the tunes CUETP8M1, CUETP8M1+MBR, and CUETP8S1. The gray band shows the total systematic uncertainty. The band around PYTHIA8 CUETP8S1 corresponds to the uncertainties of the tune parameters. The right panels show the ratio of model predictions to measured data. beam y − η  Figure 4: A comparison of the measurements of the transverse energy density, dE T /dη , at √ s = 13 TeV, as a function of shifted pseudorapidity, η = η − y beam , to the predictions and to earlier proton-proton data [6] for NSD-enhanced events at several different centre-of-mass energies. The error bars indicate the total systematic uncertainties. The beam rapidities y beam are about 9.5, 8.9, and 6.8 at √ s of 13, 7 and 0.9 TeV, respectively.

CMS NSD-enhanced
shown in Figs. 2 (upper) and 3 (upper). The gray band represents the total systematic uncertainty correlated across |η| bins. The statistical uncertainties are <1% and are not shown. In the left panel the comparison of the distribution in data and simulation is shown, while in the right panel the ratio quantifies the agreement between them. While the cosmic ray models (EPOS-LHC and QGSJETII.04) and the PYTHIA8 MONASH tune describe the data well at |η| < 4 and in the CASTOR region, they overshoot the data around |η| ≈ 4.5. This effect is most pronounced in QGSJETII.04. The PYTHIA8 CUET tunes describe the data slightly better, but have a tendency to undershoot the data towards |η| < 3.5. The band around PYTHIA8 CUETP8S1 in Fig. 3 (upper) indicates the typical uncertainties due to the tune parameters. The best description of the data is provided by the PYTHIA8 tune CUETP8S1. When MPIs are switched off in PYTHIA8 more than half of the measured energy is missing, with a slight dependence on η.
In Figs. 2 (middle) and 3 (middle) the energy density measurements are compared with predictions for the NSD-enhanced event selection. The differences between the model predictions are smaller compared with the INEL category. The EPOS-LHC and QGSJETII.04 hadronic event generators overshoot the measurement only at |η| ≈ 4.5 and otherwise show a good description of the data. The PYTHIA8 tune CUETP8S1 at the upper limit of its uncertainties provides the best overall description of the data. Figure 2 (lower) shows a comparison of the energy density measurements as a function of η for SD-enhanced events to predictions from PYTHIA8 MONASH, EPOS-LHC, and QGSJETII.04. The comparison of the same data to the different PYTHIA8 tunes is shown in Fig. 3 (lower). For the SD-enhanced events the model spread becomes significantly larger. It is interesting that the EPOS-LHC and QGSJETII.04 models are both compatible with the data only at the very lower limit of the systematic uncertainties, while all PYTHIA8 tunes are consistent with the data within the uncertainties. Furthermore, the shape of all the model predictions is very similar and, in contrast to the INEL and NSD-enhanced data, consistent with the data. Finally, we observe that for the SD-enhanced events switching off MPIs in simulations has almost no impact on the model predictions. This is an indication that the influence of MPIs within the diffractive system is small, whereas MPIs between the colliding protons will quickly destroy the singlediffractive-enhanced signature. Thus, the SD-enhanced event selection is an effective way to minimise MPI effects.
For a detailed comparison to previously published energy density results at lower centre-ofmass energies [6], the event selection is adapted to match the one previously used at detector and stable-particle levels. The whole measurement is repeated for the NSD-enhanced event selection with the requirement of at least one charged particle on both sides of the interaction point in the pseudorapidity range 3.9 < |η| < 4.4. This is combined with a reduced energy threshold of 4 GeV to ensure consistency. Finally, for all calculations the transverse energy E T = E cosh(η) per tower is used instead of just the tower energy E. In Fig. 4 the resulting corrected transverse energy density, dE T /dη , is compared to earlier published CMS data at lower √ s and to model predictions, as a function of the shifted pseudorapidity variable η = η − y beam for NSD-enhanced events. The analysis presented here uses the latest CMS detector description in the simulations, which includes an improved knowledge of the HF nonuniformity due to nonsensitive areas [34], that was not present in the original publication [6]. In order to facilitate the direct comparison of the current analysis with earlier results [6], corrections are applied to the published data that cause the results in the HF to be shifted in an η-dependent way; from about −2% at |η| = 3 to about −15% at |η| = 5, which is within the experimental uncertainties of these data.
A comparison of the model predictions and data at different √ s is shown in Fig. 4. Both the data and the model predictions are shifted by the beam rapidity to η = η − y beam . The observed behaviour is consistent with the measurements of earlier experiments in proton-proton and heavy ion collisions (e.g. [34]). At η ≈ 0 the transverse energy density does not depend on √ s, which is in agreement with the hypothesis of limiting fragmentation.

Summary
The energy density, dE/dη, is measured in the pseudorapidity range −6.6 < η < −5.2 and 3.15 < |η| < 5.20. Special low-luminosity data recorded by the CMS experiment during proton-proton collisions at the centre-of-mass energy √ s = 13 TeV are analysed for this purpose. The data are presented at the stable-particle level to allow a straightforward comparison to any theory prediction or model simulation. The measurements are compared to models tuned to describe high-energy hadronic interactions (PYTHIA8) and to the predictions of models used in cosmic ray physics (EPOS-LHC, QGSJETII.04) for inclusive inelastic (INEL), non-singlediffractive-enhanced (NSD-enhanced) , and single-diffractive-enhanced (SD-enhanced) event selections.
It is shown that the INEL and NSD-enhanced data are extremely sensitive to multi-parton interactions, while the SD-enhanced are essentially unaffected. The shape of the measured η dependencies suggest a difference in the models compared to the data. However, the predictions of PYTHIA8 tune CUETP8S1 are in satisfactory agreement with all measurements when the experimental and tune uncertainties are combined. The EPOS-LHC and QGSJETII.04 models exhibit the largest differences when compared to the single-diffractive-enhanced results.
At high energies, the hypothesis of limiting fragmentation [9, 10] assumes a longitudinal scaling behaviour in terms of shifted pseudorapidity η = η − y beam (where y beam is the beam rapidity) and thus soft-particle production in the projectile fragmentation region, η ≈ 0, is predicted to be independent of the centre-of-mass energy. This is studied by measuring the transverse energy density dE T /dη, with E T = E cosh(η), and comparing it to measurements performed in proton-proton collisions at different centre-of-mass energies. The predictions of the EPOS-LHC and QGSJETII.04 models nicely describe the combined data in the forward pseudorapidity range close to the projectile fragmentation region. The result supports the mechanism of limiting fragmentation. Since this predicts the independence of very forward particle production on the energy of the projectile particle, these data are very important for the modelling of ultra-high energy interactions that typically occur in cosmic ray collisions.

Acknowledgments
We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centres and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses.  [5] CMS Collaboration, "First measurement of the underlying event activity at the LHC with √ s = 0.9 TeV", Eur. Phys. J. C 70 (2010) 555, doi:10.1140/epjc/s10052-010-1453-9, arXiv:1006.2083.