J H E P 0 8 ( 2 0 1 1 ) 0 8 6 Published for SISSA by Springer Received: April 18, 2011 Accepted: August 4, 2011 Published: August 19, 2011 Charged particle transverse momentum spectra in pp collisions at √ s = 0.9 and 7 TeV The CMS collaboration Abstract: The charged particle transverse momentum (pT) spectra are presented for pp collisions at √ s = 0.9 and 7 TeV. The data samples were collected with the CMS detector at the LHC and correspond to integrated luminosities of 231µb−1 and 2.96 pb−1, respectively. Calorimeter-based high-transverse-energy triggers are employed to enhance the statistical reach of the high-pT measurements. The results are compared with leading and next-to- leading order QCD and with an empirical scaling of measurements at different collision energies using the scaling variable xT ≡ 2pT/ √ s over the pT range up to 200 GeV/c. Using a combination of xT scaling and direct interpolation at fixed pT, a reference transverse momentum spectrum at √ s = 2.76 TeV is constructed, which can be used for studying high-pT particle suppression in the dense QCD medium produced in heavy-ion collisions at that centre-of-mass energy. Keywords: Hadron-Hadron Scattering Open Access, Copyright CERN, for the benefit of the CMS collaboration doi:10.1007/JHEP08(2011)086 http://dx.doi.org/10.1007/JHEP08(2011)086 J H E P 0 8 ( 2 0 1 1 ) 0 8 6 Contents 1 Introduction 1 2 The CMS detector 2 3 Event selection 3 4 Primary vertex 5 5 Track selection 7 6 Event classification by leading-jet energy 7 7 Corrections and systematic uncertainties 10 8 Results 12 9 Interpolation to 2.76 TeV 16 10 Summary 17 The CMS collaboration 22 1 Introduction The charged particle transverse momentum (pT) spectrum is an important observable for understanding the fundamental quantum chromodynamic (QCD) interactions involved in proton-proton collisions. While the energy dependence of the bulk of particle production with pT below a few GeV/c is typically described either empirically or with phenomeno- logical models, the rest of the spectrum can be well described by a convolution of parton distribution functions, the hard-scattering cross section from perturbative calculations, and fragmentation functions. Such a prescription has been generally successful over a large range of lower energy pp and pp̄ collisions [1–7]. Along with measurements of the jet production cross section and fragmentation functions, measurements of high-pT spec- tra provide a test of factorised perturbative QCD (pQCD) [8] at the highest collision energy to date. In addition to its relevance to the understanding of pQCD, the charged particle spec- trum in pp collisions will be an important reference for measurements of high-pT particle suppression in the dense QCD medium produced in heavy-ion collisions. At the Rela- tivistic Heavy Ion Collider (RHIC), the sizable suppression of high-pT particle production, compared to the spectrum expected from a superposition of a corresponding number of – 1 – J H E P 0 8 ( 2 0 1 1 ) 0 8 6 pp collisions, was one of the first indications of strong final-state medium effects [9–12]. A similar measurement of nuclear modification to charged particle pT spectra has been one of the first heavy-ion results at the Large Hadron Collider (LHC) [13]. The reference spectrum for the PbPb collisions at √sNN = 2.76 TeV per nucleon can be constrained by interpolating between the pp spectra measured at √ s = 0.9 and 7 TeV. In this paper, the phase-space-invariant differential yield E d3Nch/dp 3 is presented for primary charged particles with energy (E) and momentum (p), averaged over the pseu- dorapidity acceptance of the Compact Muon Solenoid (CMS) tracking system (|η| < 2.4). The pseudorapidity is defined as –ln[tan(θ/2)], with θ being the polar angle of the charged particle with respect to the counterclockwise beam direction. The number of primary charged particles (Nch) is defined to include decay products of particles with proper life- times less than 1 cm. Using the integrated luminosities calculated in refs. [14, 15] with an estimated uncertainty of 11% and 4% at √ s = 0.9 and 7 TeV, respectively, the differential cross sections are constructed and compared to a scaling with the variable xT ≡ 2pT/ √ s. Such a scaling has already been observed for pp̄ measurements at lower collision ener- gies [4, 5, 16, 17]. For consistency with the CDF measurements at √ s = 0.63, 1.8, and 1.96 TeV, the pseudorapidity range of the xT distributions has been restricted to |η| < 1.0. Finally, using the new measurements presented in this paper, as well as previously measured pp and pp̄ cross sections, an estimate of the differential transverse momentum cross section is constructed at the interpolated energy of √ s = 2.76 TeV, corresponding to the nucleon-nucleon centre-of-mass energy of PbPb collisions recorded at the LHC. The paper is organised as follows: section 2 contains a description of the CMS detector; section 3 describes the trigger and event selection; sections 4 and 5 detail the reconstruction and selection of primary vertices and tracks; section 6 explains the characterisation of events based on the leading-jet transverse energy; section 7 describes the various applied corrections and systematic uncertainties; section 8 presents the final invariant differential yields and comparisons to data and simulation; and section 9 discusses the interpolation procedures used to construct a reference spectrum at √ s = 2.76 TeV. 2 The CMS detector A detailed description of the CMS experiment can be found in ref. [18]. The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing an axial magnetic field of 3.8 T. Immersed in the magnetic field are the pixel tracker, the silicon strip tracker, the lead tungstate crystal electromagnetic calorimeter (ECAL), and the brass/scintillator hadron calorimeter (HCAL). Muons are measured in gas ionisation detectors embedded in the steel return yoke. The CMS experiment uses a right-handed coordinate system, with the origin at the nominal interaction point, the x axis pointing to the centre of the LHC ring, the y axis pointing up perpendicular to the plane of the LHC, and the z axis along the counterclock- wise beam direction. The azimuthal angle, φ, is measured in the (x, y) plane. The tracker consists of 1440 silicon pixel and 15 148 silicon strip detector modules and measures charged particle trajectories within the nominal pseudorapidity range |η| < 2.4. – 2 – J H E P 0 8 ( 2 0 1 1 ) 0 8 6 The pixel tracker consists of three 53.3 cm-long barrel layers and two endcap disks on each side of the barrel section. The innermost barrel layer has a radius of 4.4 cm, while for the second and third layers the radii are 7.3 cm and 10.2 cm, respectively. The tracker is designed to provide an impact parameter resolution of about 100µm and a transverse momentum resolution of about 0.7 % for 1 GeV/c charged particles at normal incidence (η = 0) [19]. The tracker was aligned as described in ref. [20] using cosmic ray data prior to the LHC commissioning. The precision achieved for the positions of the detector modules with respect to particle trajectories is 3–4µm in the barrel for the coordinate in the bending plane (φ). Two elements of the CMS detector monitoring system, the beam scintillator counters (BSC) [18, 21] and the beam pick-up timing for the experiments devices (BPTX) [18, 22], were used to trigger the detector readout. The BSCs are located at a distance of 10.86 m from the nominal interaction point (IP), one on each side, and are sensitive in the |η| range from 3.23 to 4.65. Each BSC is a set of 16 scintillator tiles. The BSC elements have a time resolution of 3 ns, an average minimum ionising particle detection efficiency of 95.7%, and are designed to provide hit and coincidence rates. The two BPTX devices, located around the beam pipe at a position of z = ±175 m from the IP, are designed to provide precise information on the bunch structure and timing of the incoming beam, with better than 0.2 ns time resolution. The two steel/quartz-fibre forward calorimeters (HF), which extend the calorimetric coverage beyond the barrel and endcap detectors to the |η| region between 2.9 and 5.2, were used for further offline selection of collision events. The detailed Monte Carlo (MC) simulation of the CMS detector response is based on geant4 [23]. Simulated events were processed and reconstructed in the same manner as collision data. 3 Event selection This analysis uses data samples collected from 0.9 and 7 TeV pp collisions in the first months of the 2010 LHC running, corresponding to integrated luminosities of (231 ± 25)µb−1 and (2.96 ± 0.12) pb−1, respectively [14, 15]. This section gives a brief description of the requirements imposed to select good events for this analysis. A more detailed description of the CMS trigger selections can be found in ref. [24]. First, a minimum bias trigger was used to select events with a signal in any of the BSC tiles, coincident with a signal from either of the two BPTX detectors, indicating the presence of at least one proton bunch crossing the interaction point. From this sample, collision events were selected offline by requiring a coincidence of BPTX signals, indicating the presence of both beams. To select preferentially non-single-diffractive (NSD) events, at least one forward calorimeter (HF) tower with energy deposition E > 3 GeV in each of the forward and backward hemispheres was required. Events with beam-halo muons crossing the detector were identified and rejected based on the time difference between BSC hits on either side – 3 – J H E P 0 8 ( 2 0 1 1 ) 0 8 6 of the interaction point. Beam-induced background events, producing anomalous numbers of low-quality tracks, were rejected by requiring that at least 25% of the charged particles reconstructed in the pixel-silicon tracking system satisfied the highPurity criterion. This criterion, described in ref. [25], consists of numerous selections on the properties of the tracks, including the normalised χ2, the compatibility with the beamline and primary ver- tices, the number of hit layers, the number of ‘3D’ layers, and the number of lost layers. The selection on the fraction of highPurity tracks was only applied to events with more than 10 tracks, providing a clean separation between real pp collisions and beam backgrounds. The remaining non-collision event fraction, determined by applying the same selections to events where only a single beam was crossing the interaction point, is estimated to be less than 2 x 10−5. Events were required to have at least one primary vertex, reconstructed according to the description in the following section from triplets of pixel hits. A further requirement, namely at least one vertex found from fully reconstructed tracks (see next section for details) with number of degrees of freedom (Ndof) greater than four, was im- posed to improve the robustness against triggered events containing multiple pp collisions, i.e., “event pileup”. The loss in event selection efficiency from the fully-reconstructed-track vertex compared to the pixel vertex alone was determined entirely from data, based on a subset of early runs with negligible event pileup. The percentage of events remaining after each selection step is presented in table 1. For a large part of the 7 TeV data collection, the minimum bias trigger paths had to be prescaled by large factors because of the increasing instantaneous luminosity of the LHC. In order to maximise the pT reach of the charged particle transverse momentum measurement at this centre-of-mass energy, two high-level trigger (HLT) paths were used that selected events with minimum uncorrected transverse jet energies (ET) of 15 and 50 GeV, based only on information from the calorimeters. While the higher threshold path was not prescaled during the 7 TeV data-taking period corresponding to the 2.96 pb−1 used in this analysis, the lower threshold path had to be prescaled for a significant fraction of this sample. The 0.9 TeV data sample consists of 6.8 million minimum bias triggered events, while the 7 TeV sample is composed of 18.7 million minimum bias events, and 1.4 (5.6) million events selected with the HLT minimum-ET values of 15 (50) GeV. The selection efficiency for NSD events was determined based on simulated events from the pythia [26] event generator (version 6.420, tune D6T [27]) that were subsequently passed through a Monte Carlo simulation of the CMS detector response. The resulting event selection efficiency as a function of the multiplicity of reconstructed charged particles is shown for 7 TeV collisions in figure 1(a). The corresponding event selection efficiency is calculated by the same technique for the 0.9 TeV data (not shown). Based on events simulated with phojet [28, 29] and pythia, the remaining fraction of single-diffractive (SD) events in the selected sample was estimated to be (5± 1)% and (6± 1)% for the 0.9 and 7 TeV data, respectively. – 4 – J H E P 0 8 ( 2 0 1 1 ) 0 8 6 Collision energy 0.9 TeV 7 TeV Selection Percentage passing each selection cut One BSC + one BPTX 100.0 100.0 BPTX coincidence 94.49 90.05 Beam halo rejection 94.08 89.83 HF coincidence 73.27 83.32 Beam background rejection 73.26 83.32 Valid pixel-track vertex 70.14 82.48 Quality full-track vertex 64.04 77.35 Table 1. Summary of event selection steps applied to the 0.9 and 7 TeV collision data sets and the percentage of events from the original minimum bias samples that remain after each step. Charged particle multiplicity 0 10 20 30 40 50 se le ct ed S D o r f se le ct ed N S D ε 0.0 0.2 0.4 0.6 0.8 1.0 1.2 CMS Simulation PYTHIA 7 TeV NSD selection efficiency pixel vertex (NSD) track vertex (NSD) Selected event SD fraction pixel vertex (SD) track vertex (SD) (a) [cm]0 PVz -15 -10 -5 0 5 10 15 [c m ] 1 P V z -15 -10 -5 0 5 10 15 1 10 210 -1 Ldt = 10.2 nb∫ = 7 TeVsCMS (b) Figure 1. (a) The efficiency (εselectedNSD in eq. (7.2)) for selecting non-single-diffractive (NSD) events as a function of the multiplicity of reconstructed charged particles in the tracker acceptance (|η| < 2.4) after applying the full event selection described in the text, including a single pixel-track vertex (filled circles) and additionally requiring a fully-reconstructed-track vertex with Ndof > 4 (open circles) as described in section 4. Also, the remaining single-diffractive (SD) fraction (f selected SD in eq. (7.2)) as a function of charged particle multiplicity for the same selections (solid and dashed lines). (b) Correlation between the z positions, z0 PV and z1 PV, of the two vertices with the most associated tracks for measured events with more than one fully-reconstructed-track vertex satisfying the quality selections. 4 Primary vertex In this analysis, two separate algorithms are employed to determine the primary vertex position. The first is a highly efficient algorithm based on pixel triplet tracks that requires a minimum of just a single track consistent with the beam-spot position. The position of the beam-spot, taken as the centre of the region where the LHC beams collide, is calculated for each LHC fill based on the average over many events of the three-dimensional fitted – 5 – J H E P 0 8 ( 2 0 1 1 ) 0 8 6 vertex positions [25]. The second vertex-finding algorithm, based on fully reconstructed tracks with hits also in the silicon strip tracker, is less efficient in selecting low-multiplicity events, but more robust in discriminating against event pileup. Since pileup is significant over the majority of the analysed data sample, only the fully-reconstructed-track vertex is used to construct the raw charged particle momentum spectra. The raw spectra are subsequently corrected for the fraction of events with fewer than four tracks (and the fraction of tracks in such low-multiplicity events), based on a subset of the event sample selected with the more efficient pixel-track vertex requirement during collision runs with negligible event pileup. To determine the z position of the pixel vertex in each event, tracks consisting of three pixel hits are constructed with a minimum pT of 75 MeV/c from a region within a transverse distance of 0.2 cm from the beam axis. The x and y positions of the pixel vertex are taken from the transverse position of the beam axis. Fitted tracks are selected based on the requirement that the transverse impact parameter is less than three times the quadratic sum of the transverse errors on the track impact parameter and the beam axis position. The selected tracks are then passed to an agglomerative algorithm [30], which iteratively clusters the tracks into vertex-candidates. The procedure is halted when the distance between nearest clusters, normalised by their respective position uncertainties, reaches 12. Only vertices consisting of at least two tracks are kept, except when the event contains a single reconstructed track, which occurs in 1.67% (0.99%) of the events at √ s = 0.9 (7) TeV. In the case of multiple vertex-candidates, only the vertex with the most associated tracks is kept. While this occurs in as many as 20% of events, the rejected vertex typically has very few associated tracks and is highly correlated in z position to the vertex with the most associated tracks. These characteristics imply that the rejected vertices are not from event pileup, but rather from tracks in the tails of the impact parameter distribution that are not agglomerated into the primary vertex. The fully-reconstructed-track vertex algorithm begins from a set of tracks selected according to their transverse impact parameter to the beam-spot (< 2 cm), number of hits (> 6), and normalised χ2 (< 20). These tracks are passed to an adaptive vertex fitter, in which tracks are assigned a weight between 0 and 1 according to their compatibility with the common vertex [25]. Quality vertices are further required to have more than four degrees of freedom (Ndof), corresponding to at least four tracks with weights of approximately one. For events with multiple reconstructed vertices passing the quality selection, the correlation between the z positions of the two vertices with the most associated tracks is shown in figure 1(b). Other than the diagonal region without multiple vertices, expected from the algorithmic parameter of at least a 1 cm separation, the uncorrelated positions of the two vertices are indicative of random event pileup. The event pileup rate is estimated from the fraction of events with multiple recon- structed vertices, after correcting for vertices that are not found because of their proximity. The beam conditions varied over the analysed minimum bias data samples, such that the corrected fraction of pileup events is in the range (0.4–7.5)%. The uncertainty on the event pileup fraction, determined from the largest correction to the multiple-vertex fraction, is a constant factor of 0.2% and 1.2% for the 0.9 and 7 TeV data, respectively. – 6 – J H E P 0 8 ( 2 0 1 1 ) 0 8 6 5 Track selection This analysis uses tracks from the standard CMS reconstruction algorithm, which consists of multiple iterations of a combinatorial track finder based on various seeding layer pat- terns [31]. After each iteration, hits belonging unambiguously to tracks in the previous step are removed from consideration for subsequent steps. In order to minimise the contribution from misidentified tracks and tracks with poor momentum resolution, a number of quality selections are applied. These include the high- Purity selection mentioned in section 3, the requirement of at least five hits on the track, the normalized χ2 per degree of freedom divided by the number of tracker layers used in the fit less than a maximum value which varies from 0.48 and 0.07 depending on η and pT, and a relative momentum uncertainty of less than 20%. Furthermore, to reject non-primary tracks (i.e., the products of weak decays and secondary interactions with detector mate- rial), only the pixel-seeded tracking iterations are used, and selections are placed on the impact parameter of the tracks with respect to the primary vertex position. Specifically, the transverse and longitudinal impact parameters are required to be less than 0.2 cm and also less than 3 times the sum in quadrature of the uncertainties on the impact parameter and the corresponding vertex position. In the case of multiple quality reconstructed ver- tices in the minimum bias event samples, tracks that pass the impact parameter selections with respect to any vertex are used in the analysis. The number of events, by which the track pT distribution is normalised, is then scaled by a factor to account for the event pileup fraction. In contrast, for the jet-triggered samples, tracks are selected based on the impact parameter with respect to the single vertex responsible for the trigger. The primary vertex of the hard-scattering process is identified as the vertex with the largest value of∑ p2 T for the associated fitted tracks. With the above-mentioned selections applied to the reconstructed tracks, the algo- rithmic efficiency determined from simulated pythia events is greater than 85% (80%) for tracks with transverse momentum above 2.0 (0.4) GeV/c averaged over |η| < 2.4 (fig- ure 2(a)). In the same kinematic region, misidentified and non-primary tracks are each below 1%, while multiple reconstruction occurs for less than 0.01% of tracks. 6 Event classification by leading-jet energy All events in this analysis are classified according to the transverse energy of the most ener- getic reconstructed jet, defined as the leading jet. Jets are reconstructed from calorimeter deposits alone using the anti-kT algorithm [32] with cone radius R = √ (∆φ)2 + (∆η)2 = 0.5. The measured energy of the jet is adjusted according to corrections based on a MC description of the CMS calorimeter response with a 3–6% uncertainty on the jet energy scale [33]. The motivation for classifying events according to the leading-jet transverse energy is twofold. First, the degrading effect of the local-track density on the high-pT tracking performance (e.g., inside a jet) can be parametrised according to this variable. Based on events simulated with pythia in minimum bias and QCD samples with various thresholds – 7 – J H E P 0 8 ( 2 0 1 1 ) 0 8 6 η -2 -1 0 1 2 A lg or ith m ic e ffi ci en cy 0.5 0.6 0.7 0.8 0.9 1 PYTHIA 7 TeV > 0.4 GeV/c T p > 2.0 GeV/c T p CMS Simulation(a) [GeV/c] T p 1 10 210 tr ε × A 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 PYTHIA 7 TeV <20 GeVT 0