Abstract
This paper presents measurements of distributions of charged particles which are produced in proton–proton collisions at a centre-of-mass energy of \(\sqrt{s} = 8 \mathrm {\, TeV}\) and recorded by the ATLAS detector at the LHC. A special dataset recorded in 2012 with a small number of interactions per beam crossing (below 0.004) and corresponding to an integrated luminosity of 160 \(\upmu \mathrm{b}^{-1}\) was used. A minimum-bias trigger was utilised to select a data sample of more than 9 million collision events. The multiplicity, pseudorapidity, and transverse momentum distributions of charged particles are shown in different regions of kinematics and charged-particle multiplicity, including measurements of final states at high multiplicity. The results are corrected for detector effects and are compared to the predictions of various Monte Carlo event generator models which simulate the full hadronic final state.
Similar content being viewed by others
1 Introduction
Measurements of charged-particle spectra probe strong interactions at low momentum transfers. Such measurements have been made in lower-energy \(e^+e^-\), ep and hadron collisions [1–11] and at the CERN Large Hadron Collider (LHC) [12–23]. This paper presents measurements of multiplicity distributions, as well as transverse momentum and pseudorapidity spectra, for primary charged particles produced in pp collisions recorded by the ATLAS experiment [24] at the LHC at \(8 \mathrm {\, TeV}\) centre-of-mass energy.
Although a description of low-energy processes within a perturbative framework is not possible, predictions can be made with phenomenological models inspired by quantum chromodynamics (QCD). Data are used to constrain such models and gain further insight into the particle dynamics of the low transverse momentum regime. Furthermore, low-\(p_{\text {T}}\) processes, arising from pile-up in which there is more than one interaction per beam crossing, may also affect the topologies of events involving an interaction with a high \(p_{\text {T}}\) scale. An understanding of soft QCD processes is therefore important both in its own right and as a means of reducing systematic uncertainties in measurements of high-\(p_{\text {T}}\) phenomena.
The measurements presented in this paper use a methodology similar to that used at lower centre-of-mass energies at ATLAS [18]. Events were selected from data taken in 2012 with a trigger overlapping with the acceptance of the tracking volume. This corresponds to a minimum-bias dataset based on inelastic pp interactions. The term minimum bias is taken to refer to trigger and event selections which are as unrestrictive as possible for the pp-induced final states. The integrated luminosity of the data sample under study is 160 \(\upmu \mathrm{b}^{-1}\). Owing to improvements in understanding the material inside and around the ATLAS inner detector (ID), the uncertainties in the measured spectra are reduced by as much as 30–50 % compared to the analogous measurements at \(7 \mathrm {\, TeV}\) centre-of-mass energy [18].
The following distributions are measured:
Here, \(\eta \) is the particle’s pseudorapidity,Footnote 1 \(p_{\text {T}}\) is the component of the charged-particle momentum which is transverse to the beam direction,Footnote 2 \(n_{\mathrm {ch}}\) is the number of primary charged particles in an event, \(N_{\mathrm {ev}}\) is the event yield for a given event selection, and \(N_{\mathrm {ch}}\) is the total number of primary charged particles in all selected events in the data sample. A primary charged particle is defined as a charged particle with a mean lifetime \(\tau > 300\,\mathrm{ps}\), which is either directly produced in pp interactions or from decays of directly produced particles with \(\tau < 30\, \mathrm{ps}\); particles produced from decays of particles with \(\tau > 30\, \mathrm{ps}\) are considered as secondary particles and are thus excluded. Primary charged particles are furthermore required to satisfy the kinematic selection criteria of \(|\eta | < 2.5\) and either \(p_{\text {T}} > 100\) or \(500 \mathrm {\, MeV}\).
In order to make a more complete study of particle properties in minimum-bias events, results are given for different multiplicity and kinematic selections (termed phase spaces). The most inclusive phase spaces correspond to events with a minimum multiplicity \(n_{\mathrm {ch}} \ge 2\) or 1 and contain primary charged particles possessing a minimum transverse momentum \(p_{\text {T}} > 100\) or \(500 \mathrm {\, MeV}\), respectively. Primary-charged-particle spectra are also shown for higher-multiplicity events (\(p_{\text {T}} > 500 \mathrm {\, MeV}\), \(n_{\mathrm {ch}} \ge 6,20\) and 50) of which the latter two event types have hitherto not been measured by ATLAS. Finally, the average primary-charged-particle densities at central pseudorapidity are compared to existing measurements at different centre-of-mass energies.
2 ATLAS detector
The ATLAS detector covers almost the whole solid angle around the collision point with layers of tracking detectors, calorimeters and muon chambers. The tracking modules and the trigger system are of most relevance for the presented measurements.
The inner detector has full coverage in \(\phi \) and covers the pseudorapidity range \(|\eta | < 2.5\). It comprises a silicon pixel detector (Pixel), a silicon microstrip detector (SCT) and a transition radiation tracker (TRT). These detectors cover a sensitive radial distance from the interaction point of 50.5–150 mm, 299–560 mm and 563–1066 mm, respectively, and are immersed in a 2 T axial magnetic field provided by a solenoid. The inner-detector barrel (end-cap) parts consist of 3 (\(2\times 3\)) Pixel layers, 4 (\(2 \times 9\)) double-layers of single-sided silicon microstrips with a 40 mrad stereo angle, and 73 (\(2 \times 160\)) layers of TRT straws. Typical position resolutions are 10, 17 and 130 \(\upmu \)m for the r–\(\phi \) co-ordinate and, in the case of the Pixel and SCT, 115 and 580 \(\upmu \)m for the second measured co-ordinate. A track from a primary charged particle traversing the barrel detector would typically have 11 silicon hitsFootnote 3 (3 pixel clusters and 8 strip clusters) and more than 30 TRT straw hits.
The ATLAS detector has a three-level trigger system: Level 1 (L1), Level 2 (L2) and Event Filter (EF). For the presented measurements, the trigger relies on the L1 signals from the minimum-bias trigger scintillators (MBTS). The MBTS are positioned at each end of the detector in front of the liquid-argon end-cap calorimeter cryostats at \(z = \pm 3.56\,\mathrm{m}\). They are segmented into eight sectors in azimuth and two rings in pseudorapidity and cover the range \(2.08< |\eta | < 3.75\). The MBTS triggers are configured to require at least one or two hits above threshold from either side of the detector.Footnote 4
3 Monte Carlo simulation
The following Monte Carlo (MC) models of inclusive hadron–hadron interactions were used to generate event samples. These models employ different settings of model parameters (referred to as tunes) which were optimised to reproduce existing experimental data.
-
Pythia 8 [25] and Pythia 6 [26]. In these models, the total inelastic cross-section is separated into non-diffractive (ND) processes, dominated by t-channel gluon exchange, and diffractive processes where a colour-singlet object is exchanged. Multiple parton–parton interactions (MPI) contribute to multiplicity fluctuations and are simulated as part of the ND processes. The diffractive processes consist of single-diffractive dissociation (SD) and double-diffractive dissociation (DD). Pythia 8 is used with the A2 [27] and Monash [28] tunes. The A2 tune was performed on minimum-bias and underlying-event data, utilising the MSTW2008 LO [29] parton distribution function (PDF). The Monash tune was made using a re-analysis of fragmentation-sensitive measurements with \(e^+e^-\) collisions, combined with minimum-bias and underlying-event tuning for hadron–hadron data, utilising the NNPDF23LO PDF. Pythia 6 employs the AMBT2B [30] tune with the CTEQ6L1 [31] PDF. The AMBT2B tune was evaluated using jet and minimum-bias data.
-
Epos [32]. This model implements a parton-based Gribov–Regge theory [33], which is an effective QCD-inspired field theory describing hard and soft scattering simultaneously. Epos has been primarily designed for Pb+Pb interactions and cosmic air showers. The LHC tune [34] is used here, which modifies the modelling of radial flow to be more applicable for small-volume, high-density regions, as are found in pp interactions.
-
Qgsjet-II [35] using the default tune. This model provides a phenomenological treatment of hadronic and nuclear interactions within a Reggeon field theory framework, and includes soft and semi-hard parton processes within the “semi-hard pomeron” approach. Qgsjet-II was also developed for the simulation of cosmic rays. Qgsjet-II and Epos calculations do not rely on the standard PDFs as used in the Pythia generators.
The Pythia 8 A2, Pythia 6 AMBT2B and Epos LHC models were used to generate event samples which were processed by the Geant 4-based [36] ATLAS simulation framework [37]. The simulation also takes into account inactive and inefficient regions of the ATLAS detector. The resulting datasets were used to derive corrections for detector effects and to evaluate systematic uncertainties.
Comparisons to the data corrected to particle level are made with generated events using the Pythia 8 A2 and Monash tunes, the Epos LHC tune, and the default Qgsjet-II tune. These comparisons are shown in Sect. 10.
4 Event selection
A dedicated LHC pp run was used for which the average number of pp interactions per bunch crossing, \(\langle \mu \rangle \), was low (\(0.0028< \mu < 0.004\)). The maximum instantaneous luminosity was approximately \(1.8\times 10^{28}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\). Events were selected for which all subcomponents of the ID were operational and the solenoid magnet was on. Only events from colliding proton bunches in which the MBTS trigger recorded one or more modules above threshold on either side were considered. The MBTS trigger efficiency is described in detail in Sect. 7.1.
The following event selection criteria were applied:
-
A primary vertex with at least two associated tracks constrained to the luminous z-region of the measured beam position (termed beam spot) was required. The tracks were required to possess \(p_{\text {T}} > 100 \mathrm {\, MeV}\) and their transverse distance of closest approach to the beam spot (\(d_\mathrm {0}^{\mathrm {BS}} \)) was restricted such that \(|d_\mathrm {0}^{\mathrm {BS}} | < 4\,\mathrm{mm}\).
-
Events were rejected if they had at least one additional vertex with four or more associated tracks. Following this selection, the estimated fraction of remaining pile-up events with more than one pp interaction, based on \(\langle \mu \rangle \), was about 0.002 %. Events containing additional vertices with less than four tracks are dominated by split vertices, where the vertex reconstruction algorithm wrongly reconstructs two vertices from tracks which actually originate from a single vertex [38], and by secondary interactions being reconstructed as another primary vertex. The fraction of events with split vertices or secondary interactions which are rejected by this criterion was estimated from simulation to be 0.01 %, which is negligible and therefore ignored.
-
Depending on the phase space under study, additional selections were made on track multiplicity given the required minimum transverse momentum possessed by a track. A minimum number of selected tracks \(n_{\mathrm {sel}} \ge 2\) with transverse momentum \(p_{\text {T}} > 100 \mathrm {\, MeV}\), or \(n_{\mathrm {sel}} \ge 1\) with \(p_{\text {T}} > 500 \mathrm {\, MeV}\), which satisfy the constraints given in Sect. 5, was required.
Following the application of the above selections, the event yield is \(9.2 \times 10^6\) for the most inclusive phase space at \(n_{\mathrm {sel}} \ge 2\) and \(p_{\text {T}} > 100 \mathrm {\, MeV}\). The phase space with the lowest number of events (\(\sim \!6.4 \times 10^4\)) corresponds to \(n_{\mathrm {sel}} \ge 50\) and \(p_{\text {T}} > 500 \mathrm {\, MeV}\).
5 Track reconstruction and selection
Tracks were reconstructed using two approaches as in previous studies at \(\sqrt{s} = 7 \mathrm {\, TeV}\) [18]. Firstly, an inside-out algorithm, starting the pattern recognition from clusters in the Pixel detector, was employed. An additional algorithm with relaxed requirements on the number of silicon hits was employed to reconstruct low-momentum tracks from hits which were unused in the first approach. This latter method increases the overall efficiency of finding low-\(p_{\text {T}}\) tracks (mostly \( 100< p_{\text {T}} < 400 \mathrm {\, MeV}\)) by up to a factor of two.
To ensure that well-reconstructed tracks were used at this step, the pseudorapidity and transverse momentum must satisfy \(|\eta | < 2.5\) and \(p_{\text {T}} > 100 \mathrm {\, MeV}\). A number of further quality criteria were also applied. The track must have at least one hit in the pixel detector. A hit in the innermost layer of the pixel detector was required should the extrapolated track have passed through an active region in that layer. At least two, four or six SCT hits are required to be associated with a track for \(100 < p_{\text {T}} \le 200 \mathrm {\, MeV}\), \(200 < p_{\text {T}} \le 300 \mathrm {\, MeV}\), or \(p_{\text {T}} > 300 \mathrm {\, MeV}\), respectively. The SCT hit requirements are relaxed in the event of a track trajectory passing through inactive SCT modules. The probability of the track hypothesis being correct, estimated using the track fit \({\chi }^2\) and \(n_{\mathrm {dof}}\), was required to be greater than 0.01 for \(p_{\text {T}} > 10 \mathrm {\,GeV}\) in order to remove tracks with a mis-measured high \(p_{\text {T}}\) due to interactions with the material or combinatorial fake high-\(p_{\text {T}}\) tracks. The distance of closest approach in the transverse (\(|d_\mathrm {0}^{\mathrm {PV}} |\)) and the longitudinal plane (\(|z_\mathrm {0}^{\mathrm {PV}} \cdot \sin {\theta } |\)) was also required to be less than 1.5 mm with respect to the primary vertex. These constraints reduce the total fraction of non-primary tracks in the data from around 6 to 2 % (see Sect. 6). The average efficiency to reconstruct a track above \(p_{\text {T}} >100 \mathrm {\, MeV}\) is approximately 70 %. The efficiency of the two impact parameter requirements is around 94 %, i.e. applied together they remove approximately 6 % of all tracks that would pass the other track selection criteria.
Figure 1a shows the normalised distribution of all selected tracks as a function of pseudorapidity in the most inclusive phase space. The models reproduce the data well with discrepancies in \(\eta \) at a level of up to 3 %, which stem from the imperfect description of the \(p_{\text {T}}\) spectra (Fig. 1b) by the models, where discrepancies of up to \(\sim \!30~\%\) are visible. Figure 2 shows the normalised distribution of the fraction of all selected events as a function of track multiplicity per event.
The distributions of the average number of hits per reconstructed track in data and MC simulation as a function of pseudorapidity are shown in Fig. 10 in Appendix A, using events selected for the most inclusive phase space.
6 Backgrounds
Background events and tracks can arise from a number of sources, which are described by order of importance.
Corrections were made to the charged-particle spectra to remove the contribution of charged non-primary particles, i.e. those not originating from the pp collision. Non-primary particles are mainly due to hadronic interactions, photon conversions and weak decays. MC simulations of the shape of the \(d_\mathrm {0}^{\mathrm {PV}}\) distributions were used to quantify the fractions of primaries, non-primaries from electrons and other non-primaries which satisfy the track quality criteria. Using the same method as in Ref. [18], fits were made to the data using \(d_\mathrm {0}^{\mathrm {PV}}\) distribution templates, taken from the MC simulation, to assess the fractions of the different classes of charged particles. The fitted impact parameter distributions are shown in Fig. 11 in Appendix B. The total non-primary fraction was about 3 % in the \(100< p_{\text {T}} < 150 \mathrm {\, MeV}\) range and about 2 % at higher \(p_{\text {T}}\) values. The relative contribution of electrons (including those from rare Dalitz decays) to this fraction was about 35 % at \(p_{\text {T}} < 150 \mathrm {\, MeV}\) and dropped below 15 % with rising \(p_{\text {T}}\). Systematic uncertainties due to aspects of the template fitting method as well as the choice of MC models were added in quadrature.
In contrast to previous measurements at lower energies [14, 18], and in line with the \(13 \mathrm {\, TeV}\) measurement [23], charged particles with a mean lifetime \(30< \tau < 300\, \mathrm{ps}\) (mostly charged strange baryons) are considered to be non-stable. The reconstruction efficiency of these short-lived particles and their decay products is strongly momentum-dependent and close to zero for most particles within the measured kinematic range. However, the predicted fraction of the total generated particles associated with charged strange baryons varies with \(p_{\text {T}}\) as well as between MC models. For example, the fractions predicted by Pythia 8 A2 and Epos LHC are 5 and 13 %, respectively for \(p_{\text {T}} \sim \!5 \mathrm {\,GeV}\). To lower the model dependence on the overall track reconstruction efficiency, the contribution of such particles to the distributions under study was excluded from the measurement definition. The residual small contamination of reconstructed tracks, which is less than 0.01 % in \(\eta \) and up to 5 % at high transverse momentum (\(30< p_{\text {T}} < 50 \mathrm {\,GeV}\)), was estimated from simulation using Epos LHC and subtracted, and a systematic uncertainty was assigned following comparisons of the predictions of different MC models.
Fake tracks are reconstructed either due to detector noise or shared hits from more than one charged particle. These were estimated in simulation to be less than 0.1 % of all tracks.
Beam-induced background, i.e. beam–gas interactions and scattering from up-stream collimators, was estimated using unpaired bunches. Beam-induced backgrounds as well as pile-up contamination were reduced to a negligible level by the track-level and event-level criteria described in Sects. 4 and 5. The cosmic-ray background was found to be negligible using the techniques in Ref. [18].
7 Selection efficiencies
In order to obtain inclusive spectra for primary charged particles, the data are corrected from detector level to particle level, using corrections which account for inefficiencies due to trigger selection, vertex and track reconstruction. The methods used to obtain these efficiencies and their systematic uncertainties are described in the following sections.
7.1 Trigger efficiency
The trigger efficiency, \(\varepsilon _\mathrm {trig}\), was measured from a data sample selected using a random control trigger in coincidence with colliding bunches with a minimum requirement of two Pixel and three SCT measurements. For this efficiency, the requirement of a reconstructed primary vertex was removed from the selection of events to account for possible correlations between the trigger and vertex reconstruction efficiencies. The trigger efficiency was therefore parameterised as a function of \(n_{\mathrm {sel}}^{\mathrm {BS}}\), which is defined as the number of tracks in an event that satisfy all track quality criteria; however, instead of the nominal requirements for the impact parameters \(d_\mathrm {0}^{\mathrm {PV}} \) and \(z_\mathrm {0}^{\mathrm {PV}} \), only a constraint on the transverse impact parameter with respect to the beam spot, \(|d_\mathrm {0}^{\mathrm {BS}} | < 1.8\,\mathrm{mm}\), was applied in order to minimise correlations between the trigger and vertex efficiency corrections.
The trigger efficiency was calculated as the ratio of events triggered by the control trigger, in which the MBTS trigger also accepted the event, to the total number of triggered events in the control sample. It was determined separately for the trigger requirement in which the signal in at least one or two of the MBTS modules was above threshold. In order to maximise the recorded data yield, both triggers were combined such that either of the two triggers was required to trigger the event, in which case the corresponding trigger efficiency was applied; this was done due to pre-scale factors that evolved differently for each trigger during the run. The result for the trigger requirement in which the signal in at least one of the MBTS modules was above threshold is presented in Fig. 3a as a function of \(n_{\mathrm {sel}}^{\mathrm {BS}}\) for the most inclusive phase space. In the \(p_{\text {T}} > 500 \mathrm {\, MeV}\) phase space, the efficiency was measured to be above 98 % for \(n_{\mathrm {sel}}^{\mathrm {BS}} = 1\) and it rises more rapidly to 100 % at higher track multiplicities than in the most inclusive phase space. The efficiency for the trigger requirement in which a signal above threshold was required in at least two of the MBTS modules is lower by about 4 % for the (\(n_{\mathrm {sel}}^{\mathrm {BS}} = 2, \ p_{\text {T}} > 100 \mathrm {\, MeV}\)) event category, and about 2 % lower for events with \(n_{\mathrm {sel}}^{\mathrm {BS}} = 1\) and \(p_{\text {T}} > 500 \mathrm {\, MeV}\). It rises more slowly to 100 % as a function of \(n_{\mathrm {sel}}^{\mathrm {BS}}\) for both \(p_{\text {T}}\) requirements. These additional results are shown in Fig. 12 in Appendix C.
The trigger requirement was found to introduce no observable bias in the \(p_{\text {T}}\) and \(\eta \) distributions of selected tracks beyond the statistical uncertainties of the data recorded with the control trigger. The systematic uncertainties shown in Fig. 3a due to beam-induced background and tracks from secondary particles were estimated from differences in the trigger efficiency by varying the impact parameter requirements in the track selection. The total systematic uncertainty on the trigger efficiency in the \(n_{\mathrm {sel}}^{\mathrm {BS}} \ge 2, \ p_{\text {T}} \ge 100 \mathrm {\, MeV}\) phase space was 0.7 % for \(n_{\mathrm {sel}}^{\mathrm {BS}} = 2\), decreasing rapidly at higher track multiplicities. This uncertainty is very small compared to those from other sources.
7.2 Vertex reconstruction efficiency
The vertex reconstruction efficiency, \(\varepsilon _\mathrm {vtx}\), was determined for data and MC simulation from the ratio of selected events which satisfy the trigger requirement and contain a reconstructed vertex to the total number of triggered events. The expected contribution from beam-induced background events is estimated using the same method as described in Ref. [18] and subtracted before measuring the efficiency. The vertex reconstruction efficiency was parameterised as a function of \(n_{\mathrm {sel}}^{\mathrm {BS}}\), using the same track quality criteria with modified impact parameter constraints as for the trigger efficiency.
The result is shown in Fig. 3b as a function of \(n_{\mathrm {sel}}^{\mathrm {BS}}\) for events in the most inclusive phase space with \(p_{\text {T}} > 100 \mathrm {\, MeV}\). The efficiency was measured to be approximately 89 % for \(n_{\mathrm {sel}}^{\mathrm {BS}} = 2\), rapidly rising to 100 % at higher track multiplicities. For the \(p_{\text {T}} > 500 \mathrm {\, MeV}\) phase space, the result is given in Fig. 13a in Appendix D. For events with \(n_{\mathrm {sel}}^{\mathrm {BS}} = 2\) in the \(p_{\text {T}} > 100 \mathrm {\, MeV}\) phase space, the efficiency was parameterised as a function of the minimum difference in longitudinal impact parameter (\(\Delta z^\mathrm {min}_0\)) of track pairs, as well as the minimum transverse momentum (\(p_{\text {T}}^{\mathrm {min}}\)) of selected tracks in the event. For events with \(n_{\mathrm {sel}}^{\mathrm {BS}} = 1\) in the \(p_{\text {T}} > 500 \mathrm {\, MeV}\) phase space, the efficiency was parameterised as a function of \(\eta \) of the single track.
The systematic uncertainty was estimated from adding in quadrature the difference between the nominal vertex reconstruction efficiency, measured with beam-induced background removal, and either (1) the vertex reconstruction efficiency measured without beam-induced background removal, or (2) the vertex reconstruction efficiency with a modified impact parameter constraint. The total uncertainty is below 3 % for \(n_{\mathrm {sel}}^{\mathrm {BS}} = 2\) in the most inclusive phase space, rapidly decreasing at higher track multiplicities. This uncertainty is small compared to those from other sources, except at very low track multiplicities.
7.3 Track reconstruction efficiency
The primary track reconstruction efficiency, \(\varepsilon _\mathrm {trk}\), was determined from MC simulation and parameterised in two-dimensional bins of \(p_{\text {T}}\) and \(\eta \). This efficiency includes the efficiency of the track selection requirements (see Sect. 5). It is defined as the ratio:
where \(p_{\text {T}}^{\mathrm {gen}}\) and \(\eta ^{\mathrm {gen}}\) are properties of the generated particle, \(N_{\mathrm {rec}}^{\mathrm {matched}} (p_{\text {T}}^{\mathrm {gen}},\eta ^{\mathrm {gen}})\) is the number of reconstructed tracks matched to a generated primary charged particle in a \((p_{\text {T}}^{\mathrm {gen}},\eta ^{\mathrm {gen}})\) bin, and \(N_{\mathrm {gen}} (p_{\text {T}}^{\mathrm {gen}},\eta ^{\mathrm {gen}})\) is the number of generated primary charged particles in that bin. A track is considered matched to a generated particle if that particle has the smallest angular distance \(\Delta R\) to the track, if \(\Delta R < 0.15\), and if the particle trajectory is compatible with the position of at least one pixel hit of the associated track.
The resulting reconstruction efficiency as a function of \(\eta \) integrated over \(p_{\text {T}}\) is shown in Fig. 3c for the most inclusive phase space, and in Fig. 13b in Appendix D for the phase space given by \(n_{\mathrm {sel}} \ge 1\) and \(p_{\text {T}} > 500 \mathrm {\, MeV}\). The shape of the \(\eta \) distribution is strongly affected by the amount of material traversed by charged particles, in particular the passive material in supporting structures between the Pixel and SCT detector. A larger amount of passive material is located at high \(|\eta |\) and increases the probability of particles to undergo particle–matter interactions such as hadronic interactions, which reduces the track reconstruction efficiency. The approximately constant efficiency at \(|\eta | \sim \!2.1\) is due to the particles passing through an increasing number of layers in the ID end-cap. Figure 3d shows the efficiency as a function of \(p_{\text {T}}\) integrated over \(\eta \). The \(p_{\text {T}}\) dependence is largely due to the requirement on the minimum number of silicon hits in the track reconstruction algorithms, which is less likely to be fulfilled by lower-\(p_{\text {T}}\) tracks.
As the track reconstruction efficiency is determined from MC simulation, its systematic uncertainties result from model dependencies and from the uncertainty of the detector material description used in the simulation. Since the generated particle composition and the reconstructed track composition differs between MC tunes, a small model-dependence of the track reconstruction efficiency can be observed, leading to an additional systematic uncertainty due to the particle composition. The impact of the choice of physics models for hadronic interactions in Geant 4 simulation is also taken into account.
The amount of material within the ID was constrained to within \(\pm 5~\%\), based on extensive studies of material interactions [39]. The systematic uncertainties on the track reconstruction efficiencies were obtained by comparing the predictions of simulations which assume the nominal ID material distribution with two special simulations in which the assumed material was varied. For one simulation, the amount of non-sensitive ID material was increased by 5 % in terms of radiation length \(X_0\). In the other, the Pixel service material was increased by 10 % in \(X_0\). These studies give rise to an average systematic uncertainty on the track reconstruction efficiency of 1.6–1.7 % in the central region and up to 3.5 % in the forward region, with larger uncertainties up to 8 % for particles with very low transverse momenta of \(p_{\text {T}} < 150 \mathrm {\, MeV}\). This is the dominant contribution to the total systematic uncertainty in most regions of the measured distributions. The reduction of this uncertainty with respect to measurements at \(\sqrt{s} = 7 \mathrm {\, TeV}\), due to our improved knowledge of the ID material distribution, is about 50 % in the central region and rises to as much as 65 % in the forward region.
Systematic uncertainties due to simulation of the efficiency of the requirements on the number of hits associated with a track, the impact parameter requirements, and the efficiency of the track-fit \(\chi ^2\) probability requirement were found by comparing each selection efficiency in data and MC simulation. The sum in quadrature of these uncertainties varies between 0.5 and 1.6 % for all \(\eta \) values and \(p_{\text {T}} < 10 \mathrm {\,GeV}\), and increases to as much as 8 % for high-momentum tracks above \(p_{\text {T}} > 30 \mathrm {\,GeV}\) in the most forward regions.
The systematic uncertainty due to different fractions of positively and negatively charged tracks in data and MC simulation was found to be negligible.
The total uncertainty of the track reconstruction efficiency, shown in Fig. 3c, d, was obtained by adding all effects in quadrature and is dominated by the uncertainty from the material description.
8 Correction procedure
In order to obtain inclusive particle-level distributions, all measured detector-level distributions were corrected by an event-by-event weight, and track distributions were additionally corrected by a track-by-track weight, to compensate for the inefficiencies of the data selection and the reconstruction algorithms, as well as for contaminations due to various sources of background. Furthermore, a Bayesian unfolding procedure [40] was applied to compensate for migration and resolution effects in the observed multiplicity and transverse momentum distributions.
8.1 Event and track weights
All selected events were corrected with an event-by-event weight to compensate for the inefficiencies of the MBTS trigger selection and the vertex reconstruction algorithm. The total event weight \(w_\mathrm {ev}\) is parameterised as:
The parameter x represents a combination of \(p_{\text {T}}^{\mathrm {min}}\) for all selected tracks, the minimum difference in longitudinal impact parameter (\(\Delta z^\mathrm {min}_0\)) for track pairs, and \(\eta \) of a single track (for events with only one selected track), as described in Sect. 7.2. In addition, the MC simulation events were weighted such that the vertex z-distribution agrees with that observed in data.
Furthermore, a track-by-track weight, \(w_\mathrm {trk}(p_{\text {T}},\eta )\), was estimated for each selected track as a function of the transverse momentum and pseudorapidity assigned to the track, based on the track reconstruction efficiency, \(\varepsilon _\mathrm {trk} (p_{\text {T}},\eta )\), the fraction of non-primary tracks, \(f_\mathrm {nonp} (p_{\text {T}},\eta )\), the fraction of tracks associated with a strange baryon, \(f_\mathrm {sb} (p_{\text {T}},\eta )\), and the fraction of additional tracks corresponding to particles outside the kinematic range but migrating into the kinematic region due to resolution effects, \(f_\mathrm {okr} (p_{\text {T}},\eta )\):
The quantities \(\varepsilon _\mathrm {trk} (p_{\text {T}},\eta )\), \(f_\mathrm {sb} (p_{\text {T}},\eta )\) and \(f_\mathrm {okr} (p_{\text {T}},\eta )\) were evaluated using MC models. The quantification of the contamination due to non-primary tracks and strange baryons is described in Sect. 6.
8.2 Correction to \({\mathrm {d} N_{\mathrm {ev}}}/{\mathrm {d} n_{\mathrm {ch}}}\)
Only the event-level corrections for the trigger and vertex inefficiencies were applied to the charged-particle multiplicity distribution. Thereafter, a Bayesian unfolding process was applied to correct the observed multiplicity \(n_{\mathrm {sel}}\) to the true number of primary charged particles \(n_{\mathrm {ch}}\). This is the same procedure as was applied in Ref. [18], using five iterations in the \(p_{\text {T}} > 100 \mathrm {\, MeV}\) phase space, and four iterations for \(p_{\text {T}} > 500 \mathrm {\, MeV}\). After the unfolding, a correction was made to the resulting primary-charged-particle multiplicity distribution to account for events migrating out of the multiplicity range required by the phase space.
The corrected distribution \({\mathrm {d} N_{\mathrm {ev}}}/{\mathrm {d} n_{\mathrm {ch}}}\) was integrated over \(n_{\mathrm {ch}}\) to give the total number of events \(N_{\mathrm {ev}} \). The quantity \(N_{\mathrm {ev}} \) was then used to normalise the distributions \({1}/({2 \pi p_{\text {T}} {N_{\mathrm {ev}}}}) \cdot {\mathrm {d}^2 N_{\mathrm {ch}}}/({\mathrm {d} \eta \, \mathrm {d} p_{\text {T}}})\) and \({1}/{N_{\mathrm {ev}}} \cdot {\mathrm {d} N_{\mathrm {ch}}}/{\mathrm {d} \eta }\), as well as the multiplicity distribution itself, \({1}/{N_{\mathrm {ev}}} \cdot {{\mathrm {d} N_{\mathrm {ev}}}/{\mathrm {d} n_{\mathrm {ch}}}}\).
8.3 Correction to \({1}/({2 \pi p_{\text {T}} {N_{\mathrm {ev}}}}) \cdot {\mathrm {d}^2 N_{\mathrm {ch}}}/({\mathrm {d} \eta \, \mathrm {d} p_{\text {T}}})\) and \({1}/{N_{\mathrm {ev}}} \cdot {\mathrm {d} N_{\mathrm {ch}}}/{\mathrm {d} \eta }\)
Corrections were made for trigger requirements, vertex and track reconstruction inefficiencies, migration effects due to the resolution of reconstructed track parameters, and the influence of non-primary tracks. A Bayesian unfolding method, similar to that used to correct the \(n_{\mathrm {ch}}\) spectra, was then employed to give the \({1}/({2 \pi p_{\text {T}} {N_{\mathrm {ev}}}}) \cdot {\mathrm {d}^2 N_{\mathrm {ch}}}/({\mathrm {d} \eta \, \mathrm {d} p_{\text {T}}})\) distribution, using four iterations in the \(p_{\text {T}} > 100 \mathrm {\, MeV}\) phase space, and up to five iterations for \(p_{\text {T}} > 500 \mathrm {\, MeV}\). Fake high-\(p_{\text {T}}\) tracks are already suppressed by the \(\chi ^2\) probability requirement in the track selection, and remaining fake tracks are also unfolded for by this procedure.
8.4 Correction to \(\langle p_{\text {T}} \rangle \) versus \(n_{\mathrm {ch}}\)
The \(\langle p_{\text {T}} \rangle \) versus \(n_{\mathrm {ch}}\) distribution was evaluated in the following way. Corrections were made to two separate spectra: the distribution of the \(\Sigma _{i} p_{\text {T}} (i)\) (where the summation is made over the transverse momentum of all selected tracks in all events within a certain range of track multiplicity) versus the number of selected tracks per event, and the distribution of the sum of all selected tracks in all events within a certain range of track multiplicity versus the number of selected tracks per event. The distributions were first corrected with the appropriate track weights, which was followed by Bayesian unfolding. Finally, the ratio of the two spectra was taken to obtain the corrected \(\langle p_{\text {T}} \rangle \) versus \(n_{\mathrm {ch}}\) distribution.
9 Systematic uncertainties
In the analysis procedure, most of the individual sources of systematic uncertainty given below were applied separately as variations of the event or track weights, producing new distributions which were used to obtain alternative versions of the final corrected and unfolded results. Other sources were assessed by varying the input distributions (e.g. in \(n_{\mathrm {ch}}\) distributions, the multiplicity of each event was randomly varied with probabilities corresponding to the uncertainties on the track reconstruction efficiencies) or unfolding matrices (using statistical variations, or matrices obtained from different MC generators) which were used for the Bayesian unfolding procedure, thus producing the alternative results. In all these cases, the differences from the nominal distributions were taken as systematic uncertainties.
The following sources of systematic uncertainty in the corrected distributions were considered.
-
Incomplete knowledge of the material distribution in the ID affects the measured spectra by between 1 and 8 %. This source of systematic uncertainty is described in detail in Sect. 7.3. The total uncertainty due to the material distribution is typically less than 5 % over all distributions other than at \(p_{\text {T}} < 150 \mathrm {\, MeV}\) in the transverse momentum spectrum, at \(n_{\mathrm {ch}} \ge 120\) in the multiplicity spectrum of the \(p_{\text {T}} > 100 \mathrm {\, MeV}\) phase space, and at \(n_{\mathrm {ch}} \ge 70\) in the multiplicity spectrum of the \(p_{\text {T}} > 500 \mathrm {\, MeV}\) phase space. This is the dominant uncertainty on \({1}/{N_{\mathrm {ev}}} \cdot {\mathrm {d} N_{\mathrm {ch}}}/{\mathrm {d} \eta }\), and the leading or next-to-leading uncertainty in all other distributions.
-
Different \(p_{\text {T}}\) spectra in the MC models and data lead to differences of up to 2 % in the average track reconstruction efficiency per \(n_{\mathrm {ch}}\) interval. For the final \({\mathrm {d} N_{\mathrm {ev}}}/{\mathrm {d} n_{\mathrm {ch}}}\) spectra in the most inclusive phase space, this effect becomes as large as 12 % at the highest multiplicities.
-
The relative uncertainty on the fraction of non-primary tracks is 15 %, while the relative uncertainty on the fraction of reconstructed strange baryons is 50 %, as described in Sect. 6. The total uncertainty of both sources in the corrected distributions is 3.5 % or smaller and is not a dominant uncertainty in any region.
-
Different particle types have different reconstruction efficiencies. For example, at \(p_{\text {T}} \sim \!1 \mathrm {\,GeV}\) the reconstruction efficiency of charged pions is \(\sim \!82~\%\), whereas for kaons and protons it is \(\sim \)80 and \(\sim \)75 %, respectively. Although the MC generators give consistent efficiencies, the relative fractions of these generated particles vary between the models. For example, Pythia 8 A2 (Epos LHC) gives fractions of 77 % (72 %), 14 % (18 %) and 9 % (10 %) for generated pions, kaons and protons, respectively, at \(p_{\text {T}} \sim \!1 \mathrm {\,GeV}\). Differences in particle composition therefore lead to an uncertainty on the overall track reconstruction efficiency, which varies between 0.2 and 1 % for the corrected distributions. This is not a dominant uncertainty in any region.
-
Systematic uncertainties on the overall track reconstruction efficiency that are associated with the choice of track–particle matching algorithms (0.4 %) and the choice of physics models for MC simulation (0.3 %) are also accounted for, and are not a dominant uncertainty in any region of the corrected distributions.
-
To account for momentum resolution differences between data and MC simulation, which can arise, for example, via imperfect knowledge of the detector alignment, an uncertainty of 5 % was assigned to tracks with \(p_{\text {T}} < 150 \mathrm {\, MeV}\). At higher values of \(p_{\text {T}}\) a one-sided uncertainty of −7 % for \(10< p_{\text {T}} < 30 \mathrm {\,GeV}\) and −9 % for \(p_{\text {T}} > 30 \mathrm {\,GeV}\) tracks was assigned, as in the previous work at \(\sqrt{s} = 7 \mathrm {\, TeV}\) [18], due to the steeply falling \(p_{\text {T}}\) spectrum in combination with the lower momentum resolution in data. This is combined with a one-sided uncertainty due to the estimated fraction of mis-measured high-\(p_{\text {T}}\) tracks, which increases with transverse momentum to as much as −16 % for \(p_{\text {T}} > 30 \mathrm {\,GeV}\) tracks. The effect on the corrected distributions is typically negligible, except in the corrected \(p_{\text {T}}\) spectra.
-
Differences in the efficiencies of track quality criteria between data and MC simulation give rise to systematic uncertainties in the final spectra which are typically below 1 %, except at transverse momenta above \(10 \mathrm {\,GeV}\) and at high multiplicities, reaching as much as 6 and 5 %, respectively. However, this remains a small uncertainty compared to those from other sources in the same regions.
-
Event-level uncertainties on the trigger efficiency and vertex reconstruction efficiency give rise to systematic uncertainties of up to 3 % in the lowest multiplicity intervals of the \({\mathrm {d} N_{\mathrm {ev}}}/{\mathrm {d} n_{\mathrm {ch}}}\) spectra. However, even in these regions this uncertainty is dominated by other sources.
-
For each presented distribution, closure tests were performed. A closure test applies the full nominal correction procedure to reconstructed MC simulation events and quantifies the degree to which the generated particle-level distribution is reproduced. The degree of non-closure is typically less than 1 % and/or below the level of statistical uncertainties. Larger non-closures were found for the lower end of the \(p_{\text {T}}\) spectrum, \(100< p_{\text {T}} < 150 \mathrm {\, MeV}\), where the non-closure is found to be 6 % due to momentum resolution effects, and in the low-multiplicity region of the average transverse momentum \(\langle p_{\text {T}} \rangle \) as a function of \(n_{\mathrm {ch}}\), with up to 4 % non-closure in the \(p_{\text {T}} > 100 \mathrm {\, MeV}\) phase space due to assumptions made in the unfolding procedure. All of these non-closures were taken into account as an additional source of systematic uncertainty.
-
Uncertainties associated with the unfolding technique are estimated as the degree of non-closure following a modified correction procedure, i.e. obtained in corrected multiplicities after varying the input spectra and unfolding matrix. This is the dominant uncertainty on \({1}/({2 \pi p_{\text {T}} {N_{\mathrm {ev}}}}) \cdot {\mathrm {d}^2 N_{\mathrm {ch}}}/({\mathrm {d} \eta \, \mathrm {d} p_{\text {T}}})\) for transverse momentum values of \(p_{\text {T}} > 10 \mathrm {\,GeV}\), for which the uncertainty varies from 6 to 20 %, as well as over the entire range of \(\langle p_{\text {T}} \rangle \) versus \(n_{\mathrm {ch}}\). It is also the largest uncertainty in the low and high multiplicity regions of \({\mathrm {d} N_{\mathrm {ev}}}/{\mathrm {d} n_{\mathrm {ch}}}\), for which it has values between 1 and 12 %.
All sources of systematic uncertainty are added in quadrature, thus yielding the total systematic uncertainties which are shown as shaded areas in the figures in the next section. The total systematic uncertainties in the two most inclusive phase spaces, at \(p_{\text {T}} >100 \mathrm {\, MeV}\) (\(p_{\text {T}} > 500 \mathrm {\, MeV})\), range from 1.8 to 3.6 % (1.3 to 2.1 %) in the final \({1}/{N_{\mathrm {ev}}} \cdot {\mathrm {d} N_{\mathrm {ch}}}/{\mathrm {d} \eta }\) distributions, from 1.6 to 30 % in the final \({1}/({2 \pi p_{\text {T}} {N_{\mathrm {ev}}}}) \cdot {\mathrm {d}^2 N_{\mathrm {ch}}}/({\mathrm {d} \eta \, \mathrm {d} p_{\text {T}}})\) distributions, from 3 to 21 % (2 to 16 %) in the final \({\mathrm {d} N_{\mathrm {ev}}}/{\mathrm {d} n_{\mathrm {ch}}}\) spectra, and from 1.3 to 4 % (0.5 to 2.2 %) in the final \(\langle p_{\text {T}} \rangle \) versus \(n_{\mathrm {ch}}\) distributions. The lowest uncertainties within these ranges are found at central pseudorapidity (\(\eta =0\)), around medium transverse momentum values (\(p_{\text {T}} \sim \!1 \mathrm {\,GeV}\)), and around average multiplicity values of \(n_{\mathrm {ch}} \sim \!20\).
10 Results
Distributions of primary-charged-particle pseudorapidity, \({1}/{N_{\mathrm {ev}}} \cdot {\mathrm {d} N_{\mathrm {ch}}}/{\mathrm {d} \eta }\), are given in Figs. 4a, 5a, 6a, 7a and 8a for all measured phase spaces. The distribution corresponding to the phase space \(n_{\mathrm {ch}} \ge 2\) and \(p_{\text {T}} > 100 \mathrm {\, MeV}\) rises as \(|\eta |\) increases, peaking at \(|\eta | \sim \!2\) before falling. For the phase space \(n_{\mathrm {ch}} \ge 1\) and \(p_{\text {T}} > 500 \mathrm {\, MeV}\), the distribution is approximately constant for \(|\eta |<2\) and falls at higher \(|\eta |\). A similar shape is seen for the phase spaces requiring a higher multiplicity (\(n_{\mathrm {ch}} \ge 6,20,50\)) with the extent of the plateau becoming shorter as the multiplicity threshold is raised. A small apparent structure in the distributions of the central values of the data points occurs at values of \(|\eta | \sim \!1.7\). This is due to systematic effects in the track reconstruction efficiency which arises due to assumptions on the ID material composition, and is thus covered by the total systematic uncertainty (see Sect. 9).
The distribution corresponding to the phase space \(n_{\mathrm {ch}} \ge 2\) and \(p_{\text {T}} > 100 \mathrm {\, MeV}\) is well described by Epos LHC and Pythia 8 Monash but is underestimated by Pythia 8 A2 and Qgsjet-II.Footnote 5 For the phase space \(n_{\mathrm {ch}} \ge 1\) and \(p_{\text {T}} > 500 \mathrm {\, MeV}\), Epos LHC overestimates the distribution at values of \(|\eta | > 1.7\) and describes the data well for the rest of the pseudorapidity range. The data are overestimated by the Qgsjet-II and Pythia 8 Monash calculations and underestimated by the Pythia 8 A2 prediction. All models overestimate the overall yield for the phase spaces \(n_{\mathrm {ch}} \ge 6,20\) although Pythia 8 A2 describes the plateau in the central region well. For the largest multiplicity threshold (\(n_{\mathrm {ch}} \ge 50\)) all of the models overestimate the data at \(|\eta | > 1.7\) but provide a better description in the central region.
Figures 4b, 5b, 6b, 7b and 8b show distributions of primary-charged-particle transverse momentum, \({1}/({2 \pi p_{\text {T}} {N_{\mathrm {ev}}}}) \cdot {\mathrm {d}^2 N_{\mathrm {ch}}}/({\mathrm {d} \eta \, \mathrm {d} p_{\text {T}}})\), for various phase spaces. No model is fully consistent with the distributions, although above \(1 \mathrm {\,GeV}\) the Pythia 8 Monash predictions agree well with the data. This is also the only model which gives a fair description of the data corresponding to the highest multiplicity threshold with \(n_{\mathrm {ch}} \ge 50\) and \(p_{\text {T}} > 500 \mathrm {\, MeV}\), where all other models show large deviations as \(p_{\text {T}}\) increases. The Epos LHC predictions give the best description of the data corresponding to the phase space \(n_{\mathrm {ch}} \ge 2\) and \(p_{\text {T}} > 100 \mathrm {\, MeV}\), particularly at transverse momenta below \(1 \mathrm {\,GeV}\), while the other models underestimate the data at the lowest \(p_{\text {T}}\) values. Epos LHC provides fair predictions for the phase spaces \(n_{\mathrm {ch}} \ge 1, 6\) and \(p_{\text {T}} > 500 \mathrm {\, MeV}\), but for the higher multiplicity thresholds (\(n_{\mathrm {ch}} \ge 20\) and 50) deviations from the data are seen at high transverse momenta. Pythia 8 A2 gives fair descriptions of the data below \(6 \mathrm {\,GeV}\), yet shows deviations of up to 30 % around \(p_{\text {T}} \sim \!10 \mathrm {\,GeV}\). In all measured phase spaces, the Qgsjet-II approach shows large disagreements with the data as \(p_{\text {T}}\) increases.
In Figs. 4c and 5c distributions of primary-charged-particle multiplicity, \({1}/{N_{\mathrm {ev}}} \cdot {{\mathrm {d} N_{\mathrm {ev}}}/{\mathrm {d} n_{\mathrm {ch}}}}\), are shown for minimum transverse momentum thresholds of \(100 \mathrm {\, MeV}\) and \(500 \mathrm {\, MeV}\), respectively. For the lower threshold, the distribution rises until values of \(n_{\mathrm {ch}} \sim \!9\) before falling steeply. For the higher threshold the distribution peaks at \(n_{\mathrm {ch}} \sim \!2\). None of the models are consistent with the data although the Epos LHC model provides a fair description. The two Pythia 8 calculations predict distribution peaks which are at higher \(n_{\mathrm {ch}}\) than those observed and underestimate the event yield at low and high multiplicity. The Qgsjet-II tune overestimates the data at low and high \(n_{\mathrm {ch}}\) values and underestimates the data for intermediate \(n_{\mathrm {ch}}\) values.
The distribution of the average transverse momentum of primary charged particles, \(\langle p_{\text {T}} \rangle \), versus the primary-charged-particle multiplicity, \(n_{\mathrm {ch}}\), is given in Figs. 4d and 5d for transverse momentum thresholds of \(100 \mathrm {\, MeV}\) and \(500 \mathrm {\, MeV}\), respectively. The average \(p_{\text {T}}\) rises with multiplicity although the rise becomes progressively less steep as the multiplicity increases. This is expected due to colour coherence effects in dense parton environments, which are modelled by a colour reconnection mechanism in Pythia 8 or by the hydrodynamical evolution model used in Epos. It is assumed that numerous MPI dominate the high-multiplicity events, and that colour coherence effects thereby lead to fewer additional charged particles produced with every additional MPI, which share a higher average \(p_{\text {T}}\). The Epos LHC and Pythia 8 models provide a fair description of the data. The Qgsjet-II model fails to predict the mean transverse momentum over the entire multiplicity range, as it does not simulate colour coherence effects and therefore shows very little dependence on the multiplicity.
The evolution of the primary-charged-particle multiplicity per unit pseudorapidity at \(\eta =0\) is shown in Fig. 9. It is computed by averaging over \(|\eta |<0.2\) in the \({1}/{N_{\mathrm {ev}}} \cdot {\mathrm {d} N_{\mathrm {ch}}}/{\mathrm {d} \eta }\) distribution. In order to make consistent comparisons with previous measurements, these figures are corrected to the earlier \(\tau > 30\,\mathrm{ps}\) definition of stable particles (to include the fraction of short-lived particles which have been excluded from this study), using a factor \(1.012 \pm 0.004\) in the \(p_{\text {T}} > 100 \mathrm {\, MeV}\) phase space and \(1.025 \pm 0.008\) in the \(p_{\text {T}} > 500 \mathrm {\, MeV}\) phase spaces, derived from predictions of the Epos LHC tune with uncertainties following comparisons of the predictions of different MC models. Results are shown for the phase spaces (\(p_{\text {T}} > 500 \mathrm {\, MeV}\), \(n_{\mathrm {ch}} \ge 1\)), (\(p_{\text {T}} > 500 \mathrm {\, MeV}\), \(n_{\mathrm {ch}} \ge 6\)), and (\(p_{\text {T}} > 100 \mathrm {\, MeV}\), \(n_{\mathrm {ch}} \ge 2\)) along with available results from other ATLAS measurements at \(\sqrt{s} = 0.9, 2.36, 7\) and \(13 \mathrm {\, TeV}\) [14, 18, 23]. It can be seen that the total uncertainty in the measurement at \(\sqrt{s} = 8 \mathrm {\, TeV}\) is about 30–40 % less than for the study with \(7 \mathrm {\, TeV}\) data [18]. This was achieved due to our improved knowledge of the ID material distribution [39], which reduced the dominant source of systematic uncertainty by more than 50 % with respect to the previous \(7 \mathrm {\, TeV}\) measurement. Predictions of various QCD-based models are also shown. The best description of the data is given by Epos LHC. The predictions of the Pythia 8 tunes provide a fair description of the shape of the multiplicity dependence with centre-of-mass energy. As in the case of the other presented distributions, calculations of Qgsjet-II give the worst description.
A full summary of central primary-charged-particle densities at \(\eta =0\) in all measured phase spaces is given in Table 1, showing results obtained with the new as well as the previous fiducial definition.
11 Conclusion
Measurements were made of distributions of primary charged particles produced in minimum-bias pp collisions at \(\sqrt{s} = 8 \mathrm {\, TeV}\) with the ATLAS detector at the LHC. The results are based on a dataset corresponding to an integrated luminosity of \(160~\upmu \mathrm{b}^{-1}\). Distributions of primary-charged-particle multiplicities as well as pseudorapidity and transverse momentum spectra are shown. With the fiducial definition of primary charged particles that was used in this study (\(\tau > 300\,\mathrm{ps}\)), the central primary-charged-particle multiplicity at \(\eta = 0\) per event and unit of pseudorapidity was measured to be \(5.64 \pm 0.10\) in events containing \(n_{\mathrm {ch}} \ge 2\) primary charged particles with transverse momentum \(p_{\text {T}} > 100 \mathrm {\, MeV}\), and \(2.477 \pm 0.031\) in events with \(n_{\mathrm {ch}} \ge 1\) and \(p_{\text {T}} > 500 \mathrm {\, MeV}\). Using an extrapolation factor for short-lived charged particles with a lifetime between \(30< \tau < 300\,\mathrm{ps}\), the central primary-charged-particle multiplicity was measured to be \(5.71 \pm 0.11\) and \(2.54 \pm 0.04\), respectively. The precision of these results is 30–40 % better than for the previous highest precision ATLAS measurements at 0.9 and \(7 \mathrm {\, TeV}\). Compared with earlier studies, this paper also presents ATLAS measurements of final states at high multiplicities of \(n_{\mathrm {ch}} \ge 20\) and \(n_{\mathrm {ch}} \ge 50\). Predictions of various Monte Carlo models were compared with the data, and it was found that the best description is given by the Epos LHC tune, followed by the Pythia 8 A2 and Monash tunes. The measurements presented here are expected to provide valuable constraints for the tuning and further understanding of soft QCD physics models.
Notes
ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upwards. Cylindrical coordinates \((r,\phi )\) are used in the transverse plane, \(\phi \) being the azimuthal angle around the z-axis. The pseudorapidity is defined in terms of the polar angle \(\theta \) as \(\eta = -\ln \tan (\theta /2)\). Angular distance is measured in units of \(\Delta R \equiv \sqrt{(\Delta \eta )^{2} + (\Delta \phi )^{2}}\).
The factor \(2 \pi p_{\text {T}} \) in the \(p_{\text {T}}\) spectrum comes from the Lorentz-invariant definition of the cross-section in terms of \({\mathrm {d}^3}p\). The results could thus be interpreted as the massless approximation to \({\mathrm {d}^3}p\).
A hit is a measurement point in a sensitive detector element which can be assigned to a reconstructed track.
In order to keep the readout rates from the MBTS trigger streams approximately constant during the run in which the data were collected, pre-scale factors that evolved during the run were applied for each trigger. Therefore, a combination of both L1 MBTS triggers was used to maximise the data yield and reduce statistical uncertainties.
The MC models used here were mostly tuned to data in the \(p_{\text {T}} > 500 \mathrm {\, MeV}\) phase space (up to \(\sqrt{s} = 7 \mathrm {\, TeV}\)) and can therefore not necessarily be expected to describe the distributions at \(p_{\text {T}} > 100 \mathrm {\, MeV}\) as well as at \(p_{\text {T}} > 500 \mathrm {\, MeV}\).
References
UA1 Collaboration, G. Arnison et al., Transverse momentum spectra for charged particles at the CERN proton anti-proton collider. Phys. Lett. B 118, 167 (1982)
ABCDHW Collaboration, A. Breakstone et al., Charged multiplicity distribution in pp interactions at ISR energies. Phys. Rev. D 30, 528–535 (1984)
UA5 Collaboration, R. E. Ansorge et al., Diffraction dissociation at the CERN pulsed collider at CM energies of 900 GeV and 200 GeV. Z. Phys. C 33, 175 (1986)
UA5 Collaboration, G. J. Alner et al., UA5: A general study of proton–antiproton physics at \(\sqrt{s}\) = 546 GeV. Phys. Rep. 154, 247–383 (1987)
UA5 Collaboration, R. E. Ansorge et al., Charged particle correlations in \(\bar{p}p\) collisions at c.m. energies of 200 GeV, 546 GeV and 900 GeV. Z. Phys. C 37, 191–213 (1988)
CDF Collaboration, F. Abe et al., Transverse momentum distributions of charged particles produced in \(\bar{p}p\) interactions at \(\sqrt{s}\) = 630 GeV and 1800 GeV. Phys. Rev. Lett. 61, 1819–1822 (1988)
UA5 Collaboration, R. E. Ansorge et al., Charged Particle multiplicity distributions at 200 GeV and 900 GeV center-of-mass energy. Z. Phys. C 43, 357–374 (1989)
UA1 Collaboration, C. Albajar et al., A study of the general characteristics of \(p\bar{p}\) collisions at \(\sqrt{s} = 0.2 TeV\) to 0.9 TeV. Nucl. Phys. B 335, 261–287 (1990)
CDF Collaboration, F. Abe et al., Pseudorapidity distributions of charged particles produced in \(\bar{p}p\) interactions at \(\sqrt{s} = 630\) GeV and 1800 GeV. Phys. Rev. D 41, 2330–2333 (1990)
E735 Collaboration, T. Alexopoulos et al., Multiplicity dependence of transverse momentum spectra of centrally produced hadrons in \(\bar{p}p\) collisions at 0.3 TeV, 0.54 TeV, 0.9 TeV, and 1.8 TeV center-of-mass energy. Phys. Lett. B 336, 599–604 (1994)
CDF Collaboration, T. Aaltonen et al., Measurement of particle production and inclusive differential cross sections in \(p{\bar{p}}\) collisions at \(\sqrt{s} = 1.96\) TeV. Phys. Rev. D 79, 112005 (2009). arXiv:0904.1098 [hep-ex]. [Erratum: Phys. Rev. D 82, 119903 (2010)]
ALICE Collaboration, K Aamodt et al., First proton-proton collisions at the LHC as observed with the ALICE detector: Measurement of the charged particle pseudorapidity density at \(\sqrt{s} = 900\) GeV. Eur. Phys. J. C 65, 111–125 (2010). arXiv:0911.5430 [hep-ex]
CMS Collaboration, Transverse momentum and pseudorapidity distributions of charged hadrons in pp collisions at \(\sqrt{s} = 0.9\) and 2.36 TeV. JHEP 02, 041 (2010). arXiv:1002.0621 [hep-ex]
ATLAS Collaboration, Charged-particle multiplicities in pp interactions at \(\sqrt{s} = 900\) GeV measured with the ATLAS detector at the LHC. Phys. Lett. B 688, 21–42 (2010). arXiv:1003.3124 [hep-ex]
ALICE Collaboration, K. Aamodt et al., Charged-particle multiplicity measurement in proton–proton collisions at \(\sqrt{s} = 7\) TeV with ALICE at LHC. Eur. Phys. J. C 68, 345–354 (2010). arXiv:1004.3514 [hep-ex]
C.M.S. Collaboration, Transverse-momentum and pseudorapidity distributions of charged hadrons in pp collisions at \(\sqrt{s} = 7\) TeV. Phys. Rev. Lett. 105, 022002 (2010). arXiv:1005.3299 [hep-ex]
CMS Collaboration, Charged particle multiplicities in pp interactions at \(\sqrt{s} = 0.9,\) 2.36, and 7 TeV. JHEP 01, 079 (2011). arXiv:1011.5531 [hep-ex]
ATLAS Collaboration, Charged-particle multiplicities in pp interactions measured with the ATLAS detector at the LHC. New J. Phys. 13, 053033 (2011). arXiv:1012.5104 [hep-ex]
LHCb Collaboration, R. Aaij et al., Measurement of charged particle multiplicities and densities in pp collisions at \(\sqrt{s} = 7\) TeV in the forward region. Eur. Phys. J. C 74, 2888 (2014). arXiv:1402.4430 [hep-ex]
CMS and TOTEM Collaborations, S. Chatrchyan et al., Measurement of pseudorapidity distributions of charged particles in proton-proton collisions at \(\sqrt{s} = 8\) TeV by the CMS and TOTEM experiments. Eur. Phys. J. C 74, 3053 (2014). arXiv:1405.0722 [hep-ex]
CMS Collaboration, Pseudorapidity distribution of charged hadrons in proton-proton collisions at \(\sqrt{s} = 13\) TeV. Phys. Lett. B 751, 143–163 (2015). arXiv:1507.05915 [hep-ex]
ALICE Collaboration, J. Adam et al., Pseudorapidity and transverse-momentum distributions of charged particles in proton–proton collisions at \(\sqrt{s} = 13\) TeV. Phys. Lett. B 753, 319–329 (2016). arXiv:1509.08734 [nucl-ex]
ATLAS Collaboration, Charged-particle distributions in \(\sqrt{s} = 13\) TeV pp interactions measured with the ATLAS detector at the LHC. Phys. Lett. B 758, 67–88 (2016). arXiv:1602.01633 [hep-ex]
ATLAS Collaboration, The ATLAS experiment at the CERN Large Hadron Collider. JINST 3, S08003 (2008)
T. Sjostrand, S. Mrenna, P.Z. Skands, A brief introduction to PYTHIA 8.1. Comp. Phys. Comm. 178, 852–867 (2008). arXiv:0710.3820 [hep-ph]
T. Sjostrand, S. Mrenna, P.Z. Skands, PYTHIA 6.4 physics and manual. JHEP 05, 026 (2006). arXiv:hep-ph/0603175 [hep-ph]
ATLAS Collaboration, Further ATLAS tunes of PYTHIA6 and Pythia 8. ATL-PHYS-PUB-2011-014, 2011. https://cdsweb.cern.ch/record/1400677
P. Skands, S. Carrazza, J. Rojo, Tuning PYTHIA 8.1: the monash 2013 tune. Eur. Phys. J. C 74, 3024 (2014). arXiv:1404.5630 [hep-ph]
A.D. Martin et al., Parton distributions for the LHC. Eur. Phys. J. C 63, 189–285 (2009). arXiv:0901.0002 [hep-ph]
ATLAS Collaboration, ATLAS tunes of PYTHIA 6 and Pythia 8 for MC11, ATL-PHYS-PUB-2011-009, 2011. https://cdsweb.cern.ch/record/1363300
J. Pumplin et al., New generation of parton distributions with uncertainties from global QCD analysis. JHEP 07, 012 (2002). arXiv:hep-ph/0201195 [hep-ph]
S. Porteboeuf, T. Pierog, K. Werner, Producing hard processes regarding the complete event: the EPOS event generator. Proceedings, 45th Rencontres de Moriond on Electroweak Interactions and Unified Theories, 2010. arXiv:1006.2967 [hep-ph]
H.J. Drescher et al., Parton based Gribov–Regge theory. Phys. Rep. 350, 93–289 (2001). arXiv:hep-ph/0007198 [hep-ph]
T. Pierog et al., EPOS LHC: test of collective hadronization with data measured at the CERN Large Hadron Collider. Phys. Rev. C 92, 034906 (2015). arXiv:1306.0121 [hep-ph]
S. Ostapchenko, Monte Carlo treatment of hadronic interactions in enhanced Pomeron scheme: I. QGSJET-II model. Phys. Rev. D 83, 014018 (2011). arXiv:1010.1869 [hep-ph]
GEANT4 Collaboration, S. Agostinelli et al., GEANT4: a simulation toolkit. Nucl. Instrum. Methods A 506, 250–303 (2003)
ATLAS Collaboration, The ATLAS simulation infrastructure. Eur. Phys. J. C 70, 823–874 (2010). arXiv:1005.4568 [physics.ins-det]
ATLAS Collaboration, Performance of primary vertex reconstruction in proton–proton collisions at \(\sqrt{s} = 7\) TeV in the ATLAS experiment. ATLAS-CONF-2010-069, 2010. https://cdsweb.cern.ch/record/1281344
W. Lukas, ATLAS inner tracking detectors: Run 1 performance and developments for Run 2. Nucl. Part. Phys. Proc. 273–275, 1134–1140 (2016). doi:10.1016/j.nuclphysbps.2015.09.178. https://cdsweb.cern.ch/record/1956718
G. D’Agostini, A multidimensional unfolding method based on Bayes’ theorem. Nucl. Instrum. Methods A 362, 487–498 (1995)
Acknowledgments
We thank CERN for the very successful operation of the LHC, as well as the support staff from our institutions without whom ATLAS could not be operated efficiently. We acknowledge the support of ANPCyT, Argentina; YerPhI, Armenia; ARC, Australia; BMWFW and FWF, Austria; ANAS, Azerbaijan; SSTC, Belarus; CNPq and FAPESP, Brazil; NSERC, NRC and CFI, Canada; CERN; CONICYT, Chile; CAS, MOST and NSFC, China; COLCIENCIAS, Colombia; MSMT CR, MPO CR and VSC CR, Czech Republic; DNRF and DNSRC, Denmark; IN2P3-CNRS, CEA-DSM/IRFU, France; GNSF, Georgia; BMBF, HGF, and MPG, Germany; GSRT, Greece; RGC, Hong Kong SAR, China; ISF, I-CORE and Benoziyo Center, Israel; INFN, Italy; MEXT and JSPS, Japan; CNRST, Morocco; FOM and NWO, Netherlands; RCN, Norway; MNiSW and NCN, Poland; FCT, Portugal; MNE/IFA, Romania; MES of Russia and NRC KI, Russian Federation; JINR; MESTD, Serbia; MSSR, Slovakia; ARRS and MIZŠ, Slovenia; DST/NRF, South Africa; MINECO, Spain; SRC and Wallenberg Foundation, Sweden; SERI, SNSF and Cantons of Bern and Geneva, Switzerland; MOST, Taiwan; TAEK, Turkey; STFC, United Kingdom; DOE and NSF, United States of America. In addition, individual groups and members have received support from BCKDF, the Canada Council, CANARIE, CRC, Compute Canada, FQRNT, and the Ontario Innovation Trust, Canada; EPLANET, ERC, FP7, Horizon 2020 and Marie Skłodowska-Curie Actions, European Union; Investissements d’Avenir Labex and Idex, ANR, Région Auvergne and Fondation Partager le Savoir, France; DFG and AvH Foundation, Germany; Herakleitos, Thales and Aristeia programmes co-financed by EU-ESF and the Greek NSRF; BSF, GIF and Minerva, Israel; BRF, Norway; the Royal Society and Leverhulme Trust, United Kingdom. The crucial computing support from all WLCG partners is acknowledged gratefully, in particular from CERN and the ATLAS Tier-1 facilities at TRIUMF (Canada), NDGF (Denmark, Norway, Sweden), CC-IN2P3 (France), KIT/GridKA (Germany), INFN-CNAF (Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Taiwan), RAL (UK) and BNL (USA) and in the Tier-2 facilities worldwide.
Author information
Authors and Affiliations
Consortia
Appendices
Appendix
Average number of measurements per track
The distributions of the average number of hits per reconstructed track in data and MC simulation are shown in Fig. 10 for several detector components, using events selected by the (\(n_{\mathrm {sel}} \ge 2\), \(p_{\text {T}} > 100 \mathrm {\, MeV}\)) requirement. The distributions are shown as a function of the pseudorapidity of the reconstructed tracks. The MC simulation distributions, made with Pythia 8 using the A2 tune, have been reweighted to match the reconstructed \(p_{\text {T}}\) spectrum in data.
The level of agreement between data and MC simulation is found to be within ±1 % for the average number of measurements per track in the innermost layer of the Pixel detector (Fig. 10a), and remains within ±0.6 % for the average number of measurements per track in the whole Pixel detector (Fig. 10b) as well as the SCT (Fig. 10c). For the SCT, the sum of average measurements per track and inactive modules per track is shown. This is done in order to account for a mis-modelling of inactive modules in the MC simulation, which was shown to have negligible impact on other results presented in this paper.
Distributions of impact parameters
The normalised distributions of transverse and longitudinal impact parameters of reconstructed tracks in data and MC simulation with respect to the reconstructed primary vertex, using events selected by the (\(n_{\mathrm {sel}} \ge 2\), \(p_{\text {T}} > 100 \mathrm {\, MeV}\)) requirement, are shown in Fig. 11. The fractions of tracks originating from primary and secondary particles in the MC simulation (made with Pythia 8 using the A2 tune), which have been reweighted to match the reconstructed \(p_{\text {T}}\) spectrum as well as the fractions of reconstructed non-primary tracks in data, are also shown.
The level of agreement between data and MC simulation is found to be within ±1.5 % in the signal region of selected tracks used in the analysis, where the impact parameters are within \(\pm 1.5\,\mathrm{mm}\) of the reconstructed primary vertex. In the tail regions, the level of agreement remains within ±4 % for the transverse impact parameter and within ±9 % for the longitudinal impact parameter.
Additional trigger efficiency plots
The efficiency of the trigger requirement in which the signal in at least one of the MBTS modules was above threshold is presented in Fig. 12b for the phase space with \(n_{\mathrm {ch}} \ge 1\) and \(p_{\text {T}} > 500 \mathrm {\, MeV}\). The efficiency of the trigger requirement in which the signal in at least two of the MBTS modules was above threshold is presented in Fig. 12a for the phase space with \(n_{\mathrm {ch}} \ge 2\) and \(p_{\text {T}} > 100 \mathrm {\, MeV}\), and in Fig. 12c for the phase space with \(n_{\mathrm {ch}} \ge 1\) and \(p_{\text {T}} > 500 \mathrm {\, MeV}\). All results are shown as a function of the number of selected tracks per event, \(n_{\mathrm {sel}}^{\mathrm {BS}}\). In the phase space with the higher transverse momentum threshold, the trigger efficiency is higher and rises more quickly to 100 % for both trigger requirements.
Additional vertex and track reconstruction efficiency plots
The vertex reconstruction efficiency for events selected in the phase space with the higher transverse momentum threshold, \(p_{\text {T}} > 500 \mathrm {\, MeV}\), is presented in Fig. 13a as a function of the number of selected tracks per event, \(n_{\mathrm {sel}}^{\mathrm {BS}} \). The systematic uncertainties are found to be small in comparison with those in the most inclusive phase space.
The track reconstruction efficiency for tracks from events selected by the \(n_{\mathrm {sel}} \ge 1\) and \(p_{\text {T}} > 500 \mathrm {\, MeV}\) requirement, which was determined from MC simulation using the Pythia 8 A2 tune, is presented in Fig. 13b as a function of the pseudorapidity.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Funded by SCOAP3
About this article
Cite this article
Aad, G., Abbott, B., Abdallah, J. et al. Charged-particle distributions in pp interactions at \(\sqrt{s}=8 \mathrm {\, TeV}\) measured with the ATLAS detector. Eur. Phys. J. C 76, 403 (2016). https://doi.org/10.1140/epjc/s10052-016-4203-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-016-4203-9