Measurement of pseudorapidity distributions of charged particles in proton – proton collisions at s = 8 TeV by the CMS and TOTEM experiments

Pseudorapidity (η) distributions of charged particles produced in proton–proton collisions at a centre-of-mass energy of 8 TeV are measured in the ranges |η| < 2.2 and 5.3 < |η| < 6.4 covered by the CMS and TOTEM detectors, respectively. The data correspond to an integrated luminosity of L = 45 μb−1. Measurements are presented for three event categories. The most inclusive category is sensitive to 91–96 % of the total inelastic proton–proton cross section. The other two categories are disjoint subsets of the inclusive sample that are either enhanced or depleted in single diffractive dissociation events. The data are compared to models used to describe high-energy hadronic interactions. None of the models considered provide a consistent description of the measured distributions.


Introduction
Measurements of the yields and kinematic distributions of particles produced in proton-proton (pp) collisions at the Large Hadron Collider (LHC) can provide a better understanding of the mechanisms of hadron production in highenergy hadronic interactions. Two types of processes contribute to the production of most of the final-state particles at LHC energies. Semi-hard (multi)parton scattering, with exchanged momenta of a few GeV, is the dominant contribution. Diffractive scattering in more peripheral pp interactions, where one or both protons survive the interaction and/or are excited into a low-mass state, accounts for 15-40 % [1,2] of the pp inelastic cross section. As the particle multiplicity produced in these processes is modelled phenomenologically in the existing Monte Carlo (MC) event generators of hadronic interactions, experimental results provide an important input for tuning of the models. * e-mails: cms-publication-committee-chair@cern.ch; Totem-EB@cern.ch The results presented here focus on the charged-particle multiplicity density (dN ch /dη, also referred to as the pseudorapidity distribution) corrected down to zero transverse momentum ( p T ), in the pseudorapidity ranges |η| < 2.2 and 5.3 < |η| < 6.4, where η is defined as − ln[tan(θ/2)], with θ being the polar angle of the particle trajectory with respect to the anticlockwise-beam direction. Inclusive measurements of the η and p T distributions of charged particles have previously been performed in pp and pp collisions for different centre-of-mass energies and phase space regions [3][4][5][6][7][8][9][10][11][12][13][14].
In this paper, the data samples were collected with a minimum bias trigger generated by at least one arm of the TOTEM T2 telescopes, which also triggered the readout of the Compact Muon Solenoid (CMS). Three event samples with different final-state topologies are selected offline: a sample of inclusive inelastic pp events, a sample dominated by nonsingle diffractive dissociation (NSD) events and a sample enriched in single diffractive dissociation (SD) events. The measured data are compared to the predictions of MC event generators that model pp collider data and high-energy cosmic ray hadronic interactions.

Experimental setup
The CMS and TOTEM experiments use a right-handed coordinate system, with the origin at the nominal interaction point (IP), the x-axis pointing to the centre of the LHC ring, the y-axis pointing upwards and the z-axis pointing along the anticlockwise-beam direction. The azimuthal angle, φ, is measured in the (x, y) plane, where φ = 0 is the +x and φ = π/2 is the +y direction.
A complete description of the CMS detector can be found in Ref. [15]. The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a uniform magnetic field of 3.8 T parallel to the beam axis. Inside the magnetic field are the pixel tracker, the siliconstrip tracker, the lead tungstate electromagnetic calorime-ter and the brass/scintillator hadron calorimeter. Muons are measured in gas-ionisation detectors embedded in the steel return yoke outside the solenoid. In addition to the barrel and endcap detectors, which extend up to |η| = 3.0, the steel/quartz-fibre hadron forward calorimeters (HF) cover the region 2.9 < |η| < 5.2. The tracking detector consists of 1440 silicon pixel and 15 148 silicon strip detector modules. The barrel is composed of 3 pixel and 10 strip layers around the interaction point at radii from 4.4 to 110 cm. The forward and backward endcaps each consist of 2 pixel disks and 12 strip disks in up to 9 rings. Three of the strip rings and four of the barrel strip layers contain an additional plane, with a stereo angle of 100 mrad, to provide a measurement of the r -coordinate and z-coordinate, respectively. The tracker is designed to provide a longitudinal and transverse impact parameter resolution of about 100 µm and a p T resolution of about 0.7 % for 1 GeV/c charged particles at η = 0 [16].
The TOTEM experiment [17,18] is composed of three subdetectors: the Roman pots and the T1 and T2 telescopes. Minimum bias events are triggered by the two T2 telescopes, which are placed symmetrically on each side of the IP at about |z| =14 m. They detect charged particles produced in the polar angular range of ≈3-10 mrad (5.3 < |η| < 6.5), with full azimuthal acceptance. Each telescope consists of two half-arms, with each half-arm composed of 10 semicircular planes of triple-GEM (gas electron multiplier) chambers, arranged within a 40 cm space along the z-axis. Each chamber provides two-dimensional information on the track position covering 192 • in azimuthal angle with a small vertical overlap region between chambers of two neighbouring half-arms. Every chamber has a double-layered readout board containing two columns of 256 concentric strips (400 µm pitch, 80 µm width) to measure the radial coordinate and a matrix of 1560 pads, each covering Δη × Δφ ≈ 0.06 × 0.018 rad, to measure the azimuthal coordinate, and for triggering. The radial and azimuthal coordinate resolutions are about 110 µm and 1 • , respectively [19]. Local angles of single tracks are reconstructed with an average resolution of 0.5 mrad, and the track pseudorapidity resolution for charged particles is better than 0.05 [20], once the track is identified as coming from the vertex.
The detailed MC simulations of the CMS and TOTEM detectors are based on Geant4 [21]. Simulated events are processed and reconstructed in the same manner as collision data.

Monte Carlo models
Various MC event generators for hadronic collisions are used for data corrections and for comparison with the final, fully corrected results. The pythia6 (version 6.426) [22] generator is used with tune Z2* and pythia8 (version 8.153) [23] with tune 4C [24]. These programs provide different descriptions of the diffractive component and they both use a model [25] in which multiple partonic interactions are interleaved with parton showering. The Z2* tune [26] uses the CTEQ6L1 [27] parton distribution function (PDF) set. Tune 4C of pythia8 is based on early LHC data [24]. Parton showers in pythia are modelled according to the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) [28][29][30] prescription and hadronisation is based on the Lund string fragmentation model [31]. Diffractive cross sections are described by the Schuler-Sjöstrand model [32]. In pythia6, particle production from a low-mass state, X , with M X < 1 GeV/c 2 is treated as an isotropic two-body decay, while for high-mass states it is based on the string model. In pythia8, the same model is used to generate the cross section and the diffractive mass, but particle production is modelled differently. For lowmass states, the string model is used, but for higher masses (M X > 10 GeV/c 2 ) a perturbative description of pomeronproton scattering is introduced, based on diffractive PDFs [33][34][35], which represent probability distributions for partons in the proton under the constraint that the proton emerges intact from the collision. The non-perturbative string model introduces a mass dependence on the relative probability for a pomeron to couple to a quark or a gluon [36]. The perturbative treatment of pomeron-proton scattering results in harder p T spectra and higher multiplicity for diffractive events generated with pythia8 than for those obtained with pythia6.
The herwig++ (version 2.5.0) [37] MC event generator, with a recent tune to LHC data (UE-EE-3 with CTEQ6L1 PDFs [27,38]), is also used for comparison with the data. This generator is based on matrix element calculations similar to those used in pythia. However, herwig++ features DGLAP-based parton showers ordered in angle and uses cluster fragmentation for the hadronisation [39]. The description of hard diffractive processes also makes use of diffractive PDFs; however, soft diffraction is not implemented.
The data are also compared to predictions from two MC event generators used in cosmic ray physics . Both models include contributions from soft-and hard-parton dynamics. The soft component is described in terms of the exchange of virtual quasi-particle states, as in Gribov's Reggeon field theory [44], with multi-pomeron exchanges accounting for underlying-event effects. At higher energies, the interaction is described in terms of the same degrees of freedom, but generalised to include hard processes via hard-pomeron scattering diagrams, which are equivalent to a leading-order perturbative Quantum Chromodynamics (QCD) approach with DGLAP evolution. These models are retuned to LHC data [45], including cross section measurements by TOTEM, and charged-particle multiplicity measurements in the central region, by ALICE and CMS, at √ s = 7 TeV.

Datasets
The data were collected in July 2012 during a dedicated run with low probability (∼4 %) of overlapping pp interactions in the same bunch crossing (pileup) and a non-standard β * = 90 m optics configuration, where β * is the amplitude function of the beam at the interaction point. These data correspond to an integrated luminosity of L = MC samples were used to determine the event selection efficiency and the tracking performance. The efficiency corrections and related uncertainties for the CMS tracker are based on the pythia6, pythia8, and epos samples. The MC-based corrections for the TOTEM T2 detector and the corresponding uncertainties were determined by using pythia8 and epos, which were found to bracket the measured dN ch /dη distributions in the forward region.

Event selection and track reconstruction
The T2 track reconstruction is based on a Kalman filter-like algorithm, simplified thanks to the small amount of material in the GEM planes and the weak magnetic field in the T2 region [20]. In these conditions, the particle trajectory can be successfully reconstructed with a straight-line fit. Single tracks are reconstructed with almost 100 % efficiency for p T > 20 MeV/c, but because of multiple scattering and the magnetic field, tracks can be identified as coming from the IP with an efficiency that increases as a function of p T and is greater than 80 % for p T > 40 MeV/c [18]. The pseudorapidity of a track in T2 is defined as the average pseudorapidity of all T2 track hits, calculated from the angle between the z-axis and the line joining the hit and the nominal IP. This definition is adopted on the basis of MC simulation studies and gives an optimal estimation of the pseudorapidity of the selected primary (i.e. produced at the IP) particle. Because of the small scattering angle of the particles reconstructed in T2, the position of the vertex does not affect significantly the reconstruction of the track pseudorapidity. Due to the limited number of primary particles in T2 per event, no vertex reconstruction based on T2 information is used for the analysis. The pseudorapidity region 5.4 < |η| < 5.6 is not included in the analysis because of the large effect that the beam pipe cone at |η| ≈ 5.5 has on the propagation of the primary particles.
About 80 % of the reconstructed tracks in T2 are due to non-primary particles, hereafter referred to as secondary particles, that are mainly electrons and positrons generated by photon conversions at the edge of the HF calorimeter of CMS and in the conical section of the beam pipe at |η| ≈ 5.5. It is therefore important to discriminate these particles from the primary charged particles. The most effective primary/secondary-particle separation is achieved by using the z impact track parameter (see Fig. 1), which is defined as the z coordinate of the intersection point between the track and the plane "π 2". This is the plane which contains the z-axis and is orthogonal to the plane "π 1" defined by the z-axis and containing the track entry point in T2 [9]. This parameter is found to be stable against residual misalignment biases. Simulation studies demonstrate that the z impact distribution in the primary-particle region can be described by the sum of two Gaussian distributions plus an exponential distribution. The Gaussians, hereafter referred to as a "double-Gaussian" distribution, are mainly due to primary particles, whereas the exponential distribution accounts for most secondary particles. Simulations predict a contamination of the double-Gaussian distribution by secondary particles of about 20 %, primarily given by photons converted in the material between the IP and T2, with a contribution from decay products of strange particles; the z impact distribution of these particles is a Gaussian centred around z impact = 0. Figure 2 shows the z impact distribution at the median η of the inclusive forward pseudorapidity distribution. A com-  Fig. 2 The z impact parameter distribution measured in the data, for tracks reconstructed in one T2 half-arm in the range −5.8 < η < −5.75. A global (double-Gaussian + exponential function) fit, performed in the range from −8 to 15 m is shown by the solid curve. The dotted curve represents the exponential component from secondary particles, while the dashed curve is the double-Gaussian component, mainly due to primary tracks bined fit is performed for each η bin of the dN ch /dη distribution with the sum of a double-Gaussian and an exponential function, yielding standard deviations (amplitudes) of both Gaussian functions that increase (decrease) with η. The mean, required to be the same for both Gaussian distributions, the standard deviations and the amplitudes of the two Gaussian functions, as well as the mean and the amplitude of the exponential, are left free in the fit. The widths of the double-Gaussian distributions are consistent with the observed angular resolution of about 0.5 mrad for the T2 track reconstruction. The relative abundance of secondary particles is found to be smaller for higher |η|. The fit of the z impact distribution is also repeated by using a second degree polynomial for the description of the background. The results are found to be stable with respect to the choice of the background function. The integral of the fitting function approximates the area of the z impact distribution to within 1 %.
The T2 tracks are considered "primary candidates" if they are in the z impact range corresponding to 96 % of the area of the double-Gaussian, taken symmetrically around the mean.
The standard CMS track reconstruction algorithm is based on a combinatorial track finder (CTF) [16]. The collection of reconstructed tracks is produced by multiple iterations of the CTF track reconstruction sequence, in a process called iterative tracking. The reconstruction of the interaction vertices in the event uses the available reconstructed track collection. Prompt tracks are selected based on given quality criteria, and the selected tracks are then clustered in z using a "deterministic annealing" (DA) algorithm [48]. After identifying candidate vertices based on the DA clustering, the candidates containing at least two tracks are fitted by means of an adaptive vertex fit, where tracks in the vertex are assigned a weight between 0 and 1 based on their compatibility with the common vertex. In the central region, covered by the CMS tracker, high-purity primary tracks [49] are selected with p T > 0.1 GeV/c and relative p T uncertainty less than 10 %. To maximise the track-vertex association efficiency, a selection is applied on the track impact parameter, both along the z-axis and in the transverse plane. The impact parameter with respect to the beam spot in the transverse plane, d xy , is required to be |d xy /σ xy | < 3, while for the point of closest approach to the primary vertex along the z-direction, d z , the requirement |d z /σ z | < 3 is imposed. Here σ xy and σ z denote the uncertainties in d xy and d z , respectively. The analysis is restricted to |η| < 2.2, to avoid effects from tracks close to the geometric edge of the tracker. Events with more than one reconstructed vertex are discarded, thus reducing the effect of pileup to a negligible level (<1 %).
The pseudorapidity distributions of the charged particles are measured in the central and forward regions for three different event samples, with topologies corresponding to three different event selection criteria. An inclusive sample of events is selected by requiring at least one primarytrack candidate in T2. Event samples enhanced in non-single diffractive dissociation (NSD) and single diffractive dissociation (SD) events are also selected, the former defined by requiring a least one primary candidate in each of the two T2 telescopes and the latter by selecting events with at least one primary candidate in one T2 telescope and none in the other. Therefore, the intersection of the NSD-enhanced and SDenhanced samples is empty, while the union is the inclusive sample.
The inclusive sample includes ∼99 % of non-diffractive events. The reconstruction efficiency for diffractive events is 50 % for a diffractive mass M ∼ 3.6 GeV/c 2 and increases rapidly to 99 % for M > 10 GeV/c 2 . Most of the nondiffractive and double diffractive events, as well as the single diffractive events with masses larger than 1.3 TeV/c 2 , produce particles in both T2 telescopes and are therefore included in the NSD-enhanced sample. Simulation studies based on pythia8 and epos show that the fraction of NSD events in the SD-enhanced sample amounts to 45-65 %.

Data analysis
The pseudorapidity density measurements presented here refer to "stable" primary charged particles, with an average lifetime longer than 3×10 −11 s, either directly produced in pp collisions or from decays of particles with shorter lifetimes. Such a definition, consistent with that of previous studies [3-7,9], considers as secondary particles the decay products of K 0 S and Λ hadrons and all of the charged particles gener-

Fig. 3
Trigger efficiency as a function of the total track multiplicity in the T2 telescopes for single-sided events with primary candidates in only the z > 0 (+) or z < 0 (−) telescope and for double-sided events with primary candidates in both telescopes ated by interactions with the material in front and around the detectors.

Trigger efficiency correction
The trigger efficiency is determined with a zero bias event sample separately for events with primary-track candidates reconstructed offline in both arms of T2 and in only one arm. All zero bias data taken were used to determine the trigger inefficiency. The inefficiency of the trigger is mainly due to non-operating and noisy channels. For each event category, the trigger efficiency is calculated as a function of the total number of tracks reconstructed offline in T2, n T2 , as trig = N (n T2 ) trig /N (n T2 ) zb , where N (n T2 ) trig is the number of events with the total T2 track multiplicity n T2 passing the trigger selection and N (n T2 ) zb is the number of events with n T2 tracks selected with the zero bias trigger. The measured trigger efficiency is shown in Fig. 3.
The trigger efficiency correction, 1/ trig , is applied separately for the three event categories and is significant for events with n T2 = 1, while it approaches unity for n T2 ≥ 3. The overall trigger inefficiency results in a 0.1-0.2 % relative correction to the total dN ch /dη distributions for the three event categories. The value of the associated systematic uncertainty is conservatively assumed to be equal to this relative correction.

Event selection correction
In order to take into account the differences between the sample defined at the MC-particle level ("gen selected") and the one selected based on the reconstructed tracks ("reco selected"), a correction factor needs to be introduced. This correction is calculated for each η bin and event category from the ratio where the numerator is the pseudorapidity density obtained from the MC simulation for events selected based on the charged particles within the T2 acceptance, while the denominator is the density obtained by selecting the simulated events according to the topology of the primary candidates in T2, as explained in Sect. 5. In general, C sel is different from unity because of migrations between the different event categories. Because of misidentification of secondary particles as primaries and of track reconstruction inefficiencies, an event can be classified (according to the configuration of its reconstructed tracks) in a category that does not reflect its true charged-particle content.
For the inclusive and NSD-enhanced analysis, C sel is evaluated with pythia8 and epos. Moreover, to quantify possible biases related to this correction, the analyses are repeated with the same selection method defined in Sect. 5 but without the primary candidate z impact requirement.
As the SD-enhanced multiplicity is smaller than the NSD-enhanced multiplicity, a larger selection inefficiency is expected for the former class of events. Moreover, as the SD-enhanced sample represents only 26 % of the inclusive sample, the NSD events that are wrongly identified as SD can introduce a large bias in the measured SD-enhanced dN ch /dη distributions. Additional studies of the event selection strategy for the SD-enhanced analysis are therefore performed. The analysis is repeated with different event selection strategies and C sel is reevaluated for each as a function of η and of the track multiplicity in T2. The selection methods differ in the maximum number of tracks and of primary candidates reconstructed in T2 on each side of the IP. Simulation studies show that, depending on the method and the MC generator used, the selection efficiency is in the range of 70-90 %. The purity of the SD-enhanced sample, defined as the fraction of the selected events with primary charged particles in only one arm of T2, varies between 66 and 81 %. The dependence of the SD-enhanced dN ch /dη distributions on the selection methods is used to evaluate the systematic uncertainty related to the event selection. More details on the numerical values of C sel are given in Sects. 6.3 and 6.4.

Measurement of dN ch /dη in the central region
The charged-particle pseudorapidity distributions in the central region are obtained from the raw distributions of charged tracks after applying a number of corrections according to the formula: where S is the sample of tracks that pass the selection criteria for a given η bin, n CMS and n T2 is the total number of reconstructed tracks in the CMS tracker and T2, respectively, N evt is the number of triggered events in the corresponding track multiplicity bins, ω evt corrects for the trigger and the vertex reconstruction efficiencies, ω trk corrects for the tracking efficiency and the effect of non-primary tracks and C sel corrects for the effect of the event and primary-track selection with T2. The bin width in η is Δη = 0.2. The event correction, ω evt , depends on the track multiplicity in T2, n T2 , as well as on the multiplicity in the CMS tracker because of the minimum number of tracks required in the vertex reconstruction. It is given by where trig is the trigger efficiency (Fig. 3) and vtx is the primary vertex reconstruction and selection efficiency, calculated with pythia6 as the ratio of the number of reconstructed events satisfying the primary vertex selection to the total number of generated events. The correction for the tracking efficiency and non-primary tracks, ω trk (n CMS , p T , η), is defined as: . (4) Here, trk corrects for the geometric detector acceptance and the reconstruction tracking efficiency; the correction factor f np accounts for the fraction of non-primary tracks, i.e. secondary and misidentified tracks, while f m corrects for the effect of single charged particles that are reconstructed multiple times. These quantities are obtained from a detector simulation in bins of n CMS , p T and η. The effect of bin migrations is found to be negligible. Reconstructed events are required to pass the event selection and the generated particles are matched to the reconstructed tracks by using spatial and momentum information.
The tracking efficiency, trk , is determined as the ratio of the number of all reconstructed tracks that are matched to generated particles and satisfy the track selection criteria in an (n CMS , p T , η) bin to the number of generated primary charged particles in that bin. As shown in Fig. 4 (top), trk approaches unity for tracks with p T > 0.5 GeV/c and |η| < 1.5.
The correction for non-primary tracks, f np , is estimated as the ratio of the number of reconstructed tracks not matched  Fig. 4 Top Tracking efficiency, trk , as a function of p T and η and averaged over all multiplicity bins (n CMS ), for tracks with p T > 0.1 GeV/c and |η| < 2.2. Bottom Correction factor, f np , for non-primary tracks as a function of p T and η and averaged over all multiplicity bins (n CMS ), for tracks with p T > 0.1 GeV/c and |η| < 2.2 to a generated primary particle in a n CMS , p T , η bin to all reconstructed tracks in that bin. The correction varies as a function of η and p T of the tracks, as shown in Fig. 4 (bottom), and reaches its lowest values of about 2 % for |η| < 1.5 and p T > 0.5 GeV/c. It becomes as large as 20 % at very low transverse momentum ( p T < 0.2 GeV/c) and large pseudorapidity (|η| > 1.8).
The correction factor for multiply reconstructed particles, f m , is estimated as the ratio of the number of generated primary charged particles that are associated to multiply reconstructed tracks in a given n CMS , p T , η bin to the number of generated charged particles in that bin. It is found to be below 1 %.
The model dependence of the corrections is determined by using pythia6, epos and pythia8. The corrected data, based on correction factors derived from each generator indepen- dently, are found to differ by 1-4 % depending on the pseudorapidity bin. This amount is taken as a systematic uncertainty. The average correction factor for the event selection, C sel (η), defined in Sect. 6.2, is found to be 1.01, 1.025, and 0.94 for the inclusive, NSD-enhanced and SD-enhanced samples, respectively, independent of pseudorapidity. The correction factor is obtained from epos and pythia8. The average value of the correction factors from the two generators is applied to the data, while the relative difference between the two generators is taken as a systematic uncertainty. In addition, the event selection bias and the corresponding systematic uncertainty is estimated for each η bin as described in Sect. 6.2 by comparing the pseudorapidity distributions obtained with different methods. For the inclusive and NSDenhanced samples, the systematic uncertainty in the event selection is found to be 3-5 and 4-6 %, respectively, while for the SD-enhanced sample it is 9-16 %.
Corrections of 4-6 % are applied to the final results to extrapolate to p T = 0. The corrections are determined from the dN /d p T spectrum of primary charged particles predicted by pythia6 and pythia8. The corrections obtained from the two MC generators differ by about 3 %, resulting in a systematic uncertainty in the dN ch /dη distributions of about 0.2 %. The same corrections are also estimated from Tsallis fits [50] to the p T distributions for each η bin, giving consistent results.
A summary of the systematic uncertainties is given in Table 1. The uncertainties associated with the tracking efficiency are treated as uncorrelated between the η bins. For the inclusive and the NSD-enhanced samples, the most significant systematic uncertainties are those due to the uncertainty in the tracking efficiency and the event selection. The model dependence of C sel and the uncertainty in the event selection are dominant for the SD-enhanced sample. The total uncertainty in the tracking efficiency is estimated to be 3.9 % from a comparison of two-body and four-body D 0 decays in data and simulated samples [51]. The uncertainties related to the primary vertex selection, the trigger efficiency and pileup events are found to be around 0.1 % and are neglected.

Measurement of dN ch /dη in the forward region
The pseudorapidity density in the forward region is measured for each T2 half-arm independently, thus providing a consistency check, as each half-arm differs in its alignment and track reconstruction efficiency. The number of primary-track candidates passing the z impact parameter selection criteria is calculated for each η bin as a function of the z impact -value with the double-Gaussian and exponential fits described in Sect. 5. The fraction of tracks associated to the double-Gaussian distribution core ranges from about 74 % (lower |η| bins) to about 92 % (higher |η| bins), and is used to weight each track by the probability for the track to be primary.
Each track is also weighted by the primary-track efficiency, which depends on η and on the average pad cluster multiplicity (APM) per plane in the corresponding half-arm. This efficiency, evaluated from MC generators, is defined as the probability to successfully reconstruct, with a z impact parameter within the allowed region, a generated primary particle that falls within the acceptance of the detector. As shown in Fig. 5 for one of the T2 half-arms, the tracking efficiency decreases with increasing APM, since the reconstruction of tracks with sufficient number of hits becomes more difficult with larger occupancy. The average primary track efficiency ranges from about 73 % (lower |η| bins) to about 87 % (higher |η| bins), as shown in Fig. 5 for one of the T2 half-arms. The APM probability is a rapidly decreasing distribution, with average ∼ 24 and rms ∼ 29, for the inclusive selection. The rate of multiple associations of reconstructed tracks to the primary particle is negligible (<0.4 %) once the z impact requirement is imposed.
Conversion of photons from π 0 decays in the material between the IP and T2, as well as decay products of strange particles, also contribute to the double-Gaussian peak. The overall non-primary contribution, to be subtracted from the double-Gaussian peak, was estimated as a function of η with pythia8, epos and sibyll 2.1 [52]. The first two generators bracket the data in the forward region, while sibyll is introduced to also enclose the measurements from the LHCf The correction factors for the event selection bias (C sel (η)) are about 1.05, 1.06 and 1.00 for the inclusive, NSDenhanced and SD-enhanced event samples, respectively.
Bin migration effects in η are corrected for with pythia8, which gives the best description of the slope of the measured dN ch /dη distribution. The effects are typically at the few percent level.
Events characterised by a high T2 hit multiplicity, typically due to showers generated by particles interacting with the material before T2, are not included in the analysis. These events, for which track reconstruction capability is limited, are characterised by an APM value larger than 60 and constitute 13.5, 16.5 and 6.3 % of the inclusive, NSD-enhanced and SD-enhanced samples, respectively. The effect of removing these events is evaluated in a MC study, which yields overall correction factors of about 1.05, 1.04 and 1.06 for the inclusive, NSD-enhanced and SD-enhanced samples, respectively. To verify the stability of this correction, the analysis is also repeated after excluding events with APM values larger than 45 and re-evaluating the corresponding MC corrections. The results of the two analyses agree within 1 %. In addition, this correction is also estimated by extrapolating the measured average multiplicity obtained as a function of the maximum APM included in the sample and without correcting for the missing fraction of the sample. The extrapolation, performed with a second degree polynomial, gives a correction that is within the MC uncertainty.
The fully corrected dN ch /dη distributions in each η bin are determined from: where S j is the sample of tracks with η j − Δη/2 < η j < η j +Δη/2 satisfying the selection criteria above, Δη = 0.05 is the bin width, C sel is the correction factor related to the event selection defined in Sect. 6.2, B i j is the bin migration correction associated with the jth bin in η. In addition, Δφ/2π = 192 • /360 • is the azimuthal acceptance of each T2 half-arm, ω evt (n T2 ) ≡ 1/ trig (n T2 ) is the trigger efficiency correction, N evt (n T2 ) is the number of selected events with track multiplicity n T2 , and ω trk is defined as: where P prim is the probability for a track to be primary. Here, is the primary-track efficiency, S np is the correction factor for the non-primary contribution to the double-Gaussian peak and C mult is the correction factor for the exclusion of events with APM values above 60.
The dN ch /dη distribution thus obtained refers to charged particles with p T > 40 MeV/c, corresponding to the nominal p T acceptance of T2. A MC-based estimation obtained with epos LHC and pythia8 4C is used to correct the measurement down to p T = 0. This correction, taken from the average of the two MC predictions, is about 2 %.
The evaluation of the systematic uncertainties for the dN ch /dη distributions is performed similarly to that discussed in [9]. Details are given in the following only for the most significant contributions.
The systematic uncertainty in the P prim function, which is of order 4-5 %, is evaluated by including four effects: (a) the sensitivity to the misalignment corrections, quantified by varying the corrections within their uncertainties, (b) the sensitivity to the z impact parameter fitting range, which was changed by either one or two meters depending on the η bin, (c) the sensitivity to the background parametrisation, obtained by replacing the exponential function with a second degree polynomial and (d) the difference between the area estimated by the fitting function and the integral of the z impact distribution.
The systematic uncertainty in the primary-track efficiency is evaluated in studies where tracks are reconstructed with a set of five detector planes (out of the total of ten) in a single T2 half-arm. The track reconstruction efficiency is determined with the other set of detector planes in the same half-arm. The 5-6 % difference between the results obtained from simulation and from data is taken as an estimate of the systematic uncertainty. Table 2 Systematic and statistical uncertainties of the dN ch /dη measurement with the T2 detector in the forward region. The first two contributions are half-arm dependent and partly η-uncorrelated, while the remaining, excluding the statistical one, are half-arm independent and correlated across bins in η. The given ranges indicate the η dependence of the uncertainties Source Inclusive (%) NSD-enhanced (%) SD-enhanced (%) Tracking efficiency data-MC discrepancy 5-6 5-6 5-6 Primary selection (including alignment) 4-5 4-5 4-5 Non-primaries in the double-Gaussian peak 5 5 5 Material effects 3-6 3-6 3-6 High-multiplicity events 3 3 3 Event selection 2-3 2-3 13-15 Tracking efficiency dependence on energy spectrum and magnetic field The systematic uncertainty due to non-primary tracks included in the double-Gaussian once the exponential contribution has been removed, S np , is evaluated by considering the range of the predictions of the epos, pythia8 and sibyll MC generators, and is about 5 %.
The uncertainty in the correction for the exclusion of events with high secondary-particle multiplicity (C mult ) is taken as the difference between the epos and pythia8 predictions, which is about 3 %. The uncertainty on the correction for the event selection (C sel ) is evaluated by taking into account both the dependence of the correction from the MCs mentioned above and the dependence of the dN ch /dη results on the different event selection strategies discussed earlier. This uncertainty is 13-15 % for the SD-enhanced sample and 2-3 % for the inclusive and NSD-enhanced samples.
The possible bias due to the material uncertainty and therefore on the production of secondary tracks is evaluated as a function of η from the MC vs. data discrepancy of the ratio between the number of tracks contained in the 96 % double-Gaussian area and all the tracks in the same range. The average discrepancy is in the range of 2-6 %. Simulation studies are also performed by varying the thickness of the material in front of T2 by 40 %. This part of the material is the main source of secondary tracks that contribute to the double-Gaussian. The effect of the change of the material results in a possible bias of less than 3 %. Table 2 shows the uncertainties due to the corrections. To compute the total systematic uncertainty the errors are first separated into half-arm-correlated and uncorrelated parts and a weighted average between the four half-arms is taken. The dN ch /dη measurements obtained for the different T2 halfarms are found to be compatible.
The first two systematic uncertainties in Table 2 vary as a function of η and contribute to the uncorrelated bin-by-bin uncertainties. Conversely, the remaining systematic uncertainties affect all η bins in the same direction. The effect that systematic uncertainties might introduce in the difference of the dN ch /dη values at the beginning and at the end of the T2 acceptance is estimated to be at most 7 %.
The total uncertainty is obtained by adding in quadrature the η-correlated uncertainty and the η-uncorrelated one and the (negligible) statistical uncertainty.

Results
The combined CMS-TOTEM charged-particle pseudorapidity distributions are presented in Fig. 6 for the three event selection samples shown in Table 3. The results are derived in the central region by averaging the data points in the corresponding ±η bins and in the forward region by averaging over the four T2 half-arms. The uncertainty band represents the total uncertainty, while the bars show the statistical and uncorrelated systematics between neighbouring bins.
In the central region, the pseudorapidity density at η = 0 is 5.35 ± 0.36 for the inclusive sample, 6.20 ± 0.46 for the NSD-enhanced sample and 1.94 + 0.26 − 0.23 for the SD-enhanced sample, with negligible statistical uncertainties. The predictions from various MC event generators differ from the data by up to 20 % for the inclusive and NSD-enhanced samples, with even larger discrepancies for the SD-enhanced sample. The data are well described by pythia6 and qgsjetII-04 for the inclusive selection. For the NSD-enhanced sample, the predictions obtained from pythia6 and qgsjetII-04 agree with the data for most η bins. A good description of the mea- show the combined systematic and statistical uncertainties. The measurements are compared to results from pythia6, tune Z2*, pythia8, tune 4C, herwig++, tune UE-EE-3 with CTEQ6L1 PDFs, epos, tune LHC and qgsjetII-04 Table 3 Event selection criteria applied at the stable-particle level in the MC simulation Inclusive sample N charged particles > 0 i n 5 .3 < η < 6.5 o r −6.5 < η < −5.3, p T > 0 NSD-enhanced sample N charged particles > 0 i n 5 .3 < η < 6.5 a n d −6.5 < η < −5.3, p T > 0 SD-enhanced sample N charged particles > 0 in only 5.3 < η < 6.5 or only in −6.5 < η < −5.3, p T > 0 surement for the SD-enhanced sample is provided by both epos and pythia6.
The forward pseudorapidity density decreases with |η|. In the inclusive sample, dN ch /dη is 3.85 ± 0.49 at η = 5.375 and 2.61 ± 0.28 at η = 6.350, with negligible statistical uncertainty. The pseudorapidity density of the NSD-enhanced sample ranges between 4.80 ± 0.62 and 3.17 ± 0.35, while for the SD-enhanced sample it is in the range of 1.49 ± 0.27 to 1.20 ± 0.20. The MC predictions for the three samples differ from the data by up to about ±30 %. For the inclusive and NSD-enhanced samples, the data in the forward region are in agreement with the prediction from qgsjetII-04 and are between the epos and pythia8 results. For the SD-enhanced selection, the TOTEM data points are close to the pythia8 and her-wig++ predictions, while qgsjetII-04 underestimates the data. The change in the slope of the MC curves close to η = 5.3, more visible for the NSD-and SD-enhanced distributions, is due to the event selection requirement of at least one charged particle in the pseudorapidity region of T2.
The centre-of-mass energy dependence of the pseudorapidity distribution at η ≈ 0 is shown in Fig. 7, which includes data from various other experiments for NSD events in pp and pp collisions. Although the different experiments do not use identical event selection criteria, they all include a large fraction of NSD events. Particle production at η ≈ 0 is expected to follow a power-law dependence, dN ch /dη η=0 ∝ s , with in the range 0.14-0. 24 [40]. The result of fitting the high-energy pp and pp central-pseudorapidity particle densities with this function is shown in Fig. 7. A value of = 0.23 ± 0.01 is obtained.

Summary
Measurements of charged-particle densities over a large pseudorapidity range are presented for proton-proton (pp) collisions at a centre-of-mass energy of 8 TeV. The data were collected concurrently with the CMS and TOTEM detectors during a dedicated run with low probability for overlapping pp interactions in the same bunch crossing and correspond to an integrated luminosity of L = 45 µb −1 . Pseudorapidity distributions of charged particles within |η| < 2.2 and 5.3 < |η| < 6.4 have been measured for three event samples with different final state topologies: a sample of inclusive inelastic pp events, a sample dominated by non-single diffractive dissociation (NSD) events and a sample enriched in single diffractive dissociation (SD) events. The data are compared to theoretical predictions obtained from five different MC event generators and tunes (pythia6 Z2*, pythia8 4C, herwig++ UE-EE-3, epos LHC tune, and qgsjetII-04).
In the central region, the inclusive and NSD-enhanced samples are well described by pythia6 and qgsjetII. For the SD-enhanced sample a good description of the data is provided by both pythia6 and epos. In the forward region, the pseudorapidity distributions for the inclusive and NSDenhanced samples are between the pythia8 and epos predictions. The qgsjetII predictions are compatible with the data. The pseudorapidity distribution in the SD-enhanced sample, affected by a larger systematic uncertainty, is best described by pythia8 and herwig++.
The charged-particle densities obtained in this paper span the largest pseudorapidity interval ever measured at the LHC and have the unique potential to probe the correlation between particle production in the central region and that in the forward region. With the tunes used, none of the MC event generators are able to consistently describe the data over the whole η region and for all event samples.
Acknowledgments We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC. We especially thank the beam optics development team for the design and the successful commissioning of the high β-optics and the LHC machine coordinators for scheduling the dedicated fills. We are grateful to the technical and administrative staffs at CERN and at other CMS and TOTEM institutes for their contributions to the success of the CMS and TOTEM efforts. 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. We acknowledge the enduring support for the construction and operation of the LHC and the CMS and TOTEM detectors provided by our affiliated institutions as included in the lists of authors and the following funding agencies:  M. Macrí 6b , T. Mäki 3a , A. Mercadante 5a , N. Minafra 5b,8 , S. Minutoli 6b , F. Nemes 4,d , H. Niewiadomski 8 , E. Oliveri 7b , F. Oljemark 3a,3b , R. Orava 3a,3b , M. Oriunno 13 , K. Österberg 3a,3b , P. Palazzi 7b , Z. Petroutka 11 , J. Procházka 1a , M. Quinto 5a,5b , E. Radermacher 8 , E. Radicioni 5a , F. Ravotti 8 , L. Ropelewski 8 , G. Ruggiero 8 , H. Saarikko 3a,3b , A. Scribano 7b , J. Smajek 8 , W. Snoeys 8 , J. Sziklai 4 , C. Taylor 9 , N. Turini 7b , V. Vacek 1b , J. Welti 3a,3b , J. Whitmore 14 , P. Wyszkowski 10 , K. Zielinski 10