Correlations between jets and charged particles in PbPb and pp collisions at sqrt(s[NN]) = 2.76 TeV

The quark-gluon plasma is studied via medium-induced changes to correlations between jets and charged particles in PbPb collisions compared to pp reference data. This analysis uses data sets from PbPb and pp collisions with integrated luminosities of 166 inverse microbarns and 5.3 inverse picobarns, respectively, collected at sqrt(s[NN]) = 2.76 TeV. The angular distributions of charged particles are studied as a function of relative pseudorapidity (Delta eta) and relative azimuthal angle (Delta phi) with respect to reconstructed jet directions. Charged particles are correlated with all jets with transverse momentum (pt) above 120 GeV, and with the leading and subleading jets (the highest and second-highest in pt, respectively) in a selection of back-to-back dijet events. Modifications in PbPb data relative to pp reference data are characterized as a function of PbPb collision centrality and charged particle pt. A centrality-dependent excess of low-pt particles is present for all jets studied, and is most pronounced in the most central events. This excess of low-pt particles follows a Gaussian-like distribution around the jet axis, and extends to large relative angles of Delta eta approximately 1 and Delta phi approximately 1.


Introduction
The quark-gluon plasma (QGP) produced in ultra-relativistic heavy ion collisions may be probed in situ via partons produced in the initial hard-scatterings, which carry high transverse momentum (p T ) compared to most of the particles in the event.In the QGP, partons are expected to suffer energy loss in the medium, a phenomenon known as "jet quenching" [1].This effect was discovered at RHIC via observables including suppression of high-p T particle production [2] and charged particle correlations [3].Jet quenching has been studied at the CERN LHC via high-p T particle suppression [4][5][6], and via the momentum balance of reconstructed back-to-back dijets.In these latter studies, dijet transverse momentum balance was investigated in PbPb, pPb and pp collisions [7][8][9][10].Significant imbalance was found in central PbPb collisions, consistent with a pathlength-dependent energy loss in the QGP medium.In peripheral PbPb collisions and in pPb collisions, the dijet momentum balance is comparable to the one measured in pp collisions, thus confirming that the energy loss in central PbPb collisions is not due to initial state cold nuclear matter effects.In the QGP, the interaction of the hard-scattered partons (and their fragmentation products) with the medium leads to a redistribution of energy carried by the produced particles.Comparing the charged-particle distributions from heavy ion data to the pp reference can help to differentiate between energy loss models and ultimately constrain the properties of the QGP [11][12][13][14].
Early RHIC studies of two-particle correlations involving a leading high-p T (8)(9)(10)(11)(12)(13)(14)(15) GeV) particle did not find a significant modification in the distribution of the associated particles at small angles from the leading particle.A quenching effect was found in the distribution of particles opposite in azimuth to the leading particle, observed as a reduction in the associated yield [15][16][17][18][19]. These results could be interpreted as an in-medium energy loss (in which associated particles fully thermalize and do not retain correlation to the jet direction) followed by a vacuum-like fragmentation of the remaining jet [20].Compared to the capabilities of the LHC, these studies at RHIC are significantly limited by the lower production rates for hard probes.Subsequent high-precision measurements at the LHC [21][22][23] have shown that the detailed jet structures within a jet cone radius of 0.3 are modified by the medium in terms of both p T and angular distributions.However, these observed in-cone changes only explain a small fraction of the dijet momentum imbalance, indicating that a large amount of energy is radiated outside of the jet cone or transferred to particles with very low momentum.Direct measurements of energy redistribution between event hemispheres containing subleading and leading jets were made by CMS via the "missing-p T " observable [10,24].These studies found the overall energy flow to be modified in PbPb collisions out to large radial distances from the dijet axis, and various theoretical models have since attempted to describe the result [25][26][27].Extending measurements of jet structure modifications to similarly large angles is crucial to properly constrain the energy loss mechanism.
In this paper, we use 166 µb −1 of PbPb collisions taken during the 2011 LHC heavy ion run at a nucleon-nucleon center-of-mass energy of √ s NN = 2.76 TeV.For the reference measurement, we use pp data taken in 2013 at the same center-of-mass energy corresponding to an integrated luminosity of 5.3 pb −1 .Two-dimensional angular correlations were previously studied in CMS for pairs of charged particles [28].In the present analysis, this technique is applied to correlate jets with charged particles.For each charged particle and reconstructed jet, p T and pseudorapidity (η) are measured with respect to the beam axis, and azimuthal angle (φ) is measured in the transverse plane.From these measurements, the relative pseudorapidity (∆η = η track − η jet ) and relative azimuth (∆φ = φ track − φ jet ) between jets and charged-particle tracks is determined.These relative angles are used to construct two-dimensional ∆η-∆φ charged particle density distributions, which we will refer to as "jet-track correlations".The jet-track correla-tions are then studied as a function of centrality (defined as a percentile of the total inelastic cross section, with 0% indicating collisions with impact parameter b = 0) and charged-particle transverse momentum (p trk T ).In order to extend these measurements to low p trk T , where soft particles resulting from energy loss mechanisms such as gluon radiation are expected to appear, an analysis must carefully handle both large combinatorial backgrounds typical for the heavy ion environment, and longrange correlations arising from hydrodynamic expansion [29].Taking advantage of the kinematic reach of hard probes and the availability of detailed characterization of the event bulk properties, the CMS detector permits the statistical separation of the medium-related modifications of jet-track correlations from the long-range hydrodynamic background.Using this technique, this study captures jet-related energy flow both inside and outside of the jet cone, extending measurements of intrinsic jet properties to large relative angles in ∆η and ∆φ.

The CMS detector
The central feature of the CMS apparatus is a superconducting solenoid with an internal diameter of 6 m, providing an axial uniform magnetic field of 3.8 T. Muons are measured in gasionization detectors embedded in the steel flux-return yoke outside the solenoid.Within the solenoid volume are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter (HCAL), each composed of a barrel and two endcap sections.Extensive hadronic forward (HF) steel and quartz fiber calorimetry complements the barrel and endcap detectors, providing coverage to |η| < 5.In this analysis, the collision centrality is determined using the total sum of transverse energy (E T ) from calorimeter towers in the HF region (covering 2.9 < |η| < 5.2).For the forward region 2.9 < |η| < 5.0 relevant for collision centrality determination, HF fibers are bundled in towers with widths of 0.175 × 0.175 (∆η × ∆φ) [30].The E T distribution is used to divide the event sample into bins, each representing 0.5% of the total nucleus-nucleus hadronic interaction cross section.A detailed description of centrality determination can be found in Ref. [8].
Jet reconstruction for this analysis relies on calorimeter information from the ECAL and HCAL.For the central region |η| < 1.6 in which jets are selected for this analysis, the HCAL cells have widths of 0.087 in both η and φ.In the η-φ plane, and for |η| < 1.48, the HCAL cells map on to 5 × 5 ECAL crystals arrays to form calorimeter towers projecting radially outwards from close to the nominal interaction point.The barrel section of the ECAL has an energy resolution of 1%-2.5%, while the endcaps have an energy resolution of 2.5-4% [31].Within each tower, the energy deposits in ECAL and HCAL cells are summed to define the calorimeter tower energies, subsequently used to provide the energies and directions of hadronic jets.When combining information from the entire detector, the jet energy resolution amounts typically to 15% at 10 GeV, 8% at 100 GeV, and 4% at 1 TeV, to be compared to about 40%, 12%, and 5% obtained when the ECAL and HCAL calorimeters alone are used [30].
Accurate particle tracking is critical for measurements of charged-hadron yields.The CMS silicon tracker measures charged particles within the region |η| < 2.5.For particles of 1 < p T < 10 GeV and |η| < 1.4, the track resolutions are typically 1.5% in p T and 25-90 (45-150) µm in the transverse (longitudinal) impact parameter with respect to the collision vertex [32].A detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [33].

Jet and track reconstruction and corrections
For both pp and PbPb collisions, jet reconstruction is performed with the anti-k T algorithm, encoded in the FASTJET framework [34,35].Following previous CMS studies [10,[21][22][23], a narrow jet reconstruction distance parameter, R = 0.3, is chosen due to the large underlying event in heavy ion collisions.Jet p T and direction in η and φ are determined based on iterative clustering of energy deposits in the CMS calorimeters.For PbPb collisions, the CMS algorithm "HF/Voronoi" is used to subtract the heavy ion underlying event [36].This algorithm estimates the underlying event contribution to the E T in each calorimeter tower by performing a singular value decomposition of the particle distributions.The average E T , as a function of η and φ, is subtracted from each calorimeter tower, and then the energy is redistributed between neighboring calorimeter towers so that no tower contains non-physical negative E T .For pp collisions no underlying event subtraction is employed, as the effect of the underlying event on the jet energy is small relative to the jet energy scale (JES) uncertainty.
Monte Carlo (MC) event generators have been used for evaluation of the jet and track reconstruction performance, in particular for determining the tracking efficiency as well as the jet energy response and resolution.Jet events are generated with PYTHIA [37] (version 6.423, tune Z2 [38]).These generated PYTHIA events are propagated through the CMS detector using the GEANT4 package [39] to simulate the detector response.In order to account for the influence of the underlying PbPb events, the PYTHIA events are embedded into fully simulated PbPb events that are generated by HYDJET [40] (version 1.8), which is tuned to reproduce the total particle multiplicities, charged-hadron spectra, and elliptic flow at all centralities.The embedding is performed by mixing the simulated digital signal information from PYTHIA and HYDJET, hereafter referred to as PYTHIA+HYDJET.No simulation of jet quenching is applied in this PYTHIA+HYDJET simulation.These events are then propagated through the same reconstruction and analysis procedures used for data events.The JES is established using PYTHIA and PYTHIA+HYDJET events in bins of event centrality as a function of p T , η, and number of charged particles inside the jet cone.For studies of pp data and PYTHIA simulation, charged particles are reconstructed using the same iterative method [32] as in the previous CMS analyses of pp collisions.For PbPb data and PYTHIA+HYDJET MC, an iterative charged-particle reconstruction similar to that of earlier heavy ion analyses [5,22] is employed, as described in detail in Ref. [24].

Data samples and triggers
The first level (L1) of the CMS trigger system, composed of custom hardware processors, uses information from the calorimeters and muon detectors to select the most interesting events in a fixed time interval of less than 4 µs.The high-level trigger (HLT) processor farm further decreases the event rate from around 100 kHz to less than 1 kHz before data storage.The events for this analysis were selected using an HLT that selects all events containing at least one calorimeter jet with p T > 80 GeV.The HLT is fully efficient for events containing offline reconstructed jets with p T > 120 GeV.In order to suppress noise due to noncollision sources such as cosmic rays and beam backgrounds, the events used in this analysis were further required to satisfy offline selection criteria as documented in Refs.[8,41].These criteria include selecting only events with a reconstructed vertex including at least two tracks and a z position within 15 cm of the detector center, and requiring energy deposits in at least 3 forward calorimeter towers on either side of the interaction point.
The offline selection of events begins with jets reconstructed as described in Section 3. To study the jet-track correlations, the events are then categorized into two samples: an inclusive selection of high-p T jets and a selection of back-to-back dijet events.For the inclusive sample, jets are required to have p T > 120 GeV and to fall within |η| < 1.6, with multiple jets from the same event permitted in this inclusive jet sample.These inclusive selection criteria match previous CMS studies [22,23] that measured jet fragmentation functions and jet shapes within the jet cone (∆R < 0.3), and allow this analysis to extend comparable measurements to large angles from the jet axis.We then separately select a dijet sample using criteria matched to those of a previous analysis that explores dijet energy balance [10].In this dijet selection, events are first required to contain a highest p T (leading) calorimeter jet in the range of |η| < 2, with a corrected jet p T > 120 GeV and a next-highest p T (subleading) jet of p T > 50 GeV, also in |η| < 2. The azimuthal angle between the leading and subleading jets is required to be at least 5π/6.No explicit requirement is made either on the presence or absence of a third jet in the event.To ensure stable jet reconstruction performance, only events in which both leading and subleading jets fall within |η| < 1.6 are included in the final data sample.

Jet-track angular correlations
Charged tracks in the event with p trk T above 1 GeV are used to construct two-dimensional ∆η-∆φ correlations with respect to the individually measured jet axis for inclusive jets and for leading and subleading jets in dijet events.The jet-track correlations are constructed according to the procedure established in Ref. [42], for the following bins in p trk T : 1-2, 2-3, 3-4, and 4-8 GeV.This work does not attempt to construct correlations below 1 GeV, where the jet-related signal is very small compared to the combinatorial and long-range correlated background, or for p trk T > 8 GeV where the statistical power becomes limited.The correlations are corrected for tracking efficiency and misreconstruction on a per-track basis, using an efficiency parametrization defined as a function of centrality, p trk T , η, φ, and radial distance from the nearest jet with p T > 50 GeV [24].
Correlations are formed by measuring angular distances to the inclusive, leading and subleading jet axes for each p trk T range.The signal pair distribution, S(∆η, ∆φ), represents the per-track efficiency-corrected yield of jet-track pairs N same from the same event normalized by the total number of jets: To correct for pair acceptance effects, we use the mixed event technique [28,43,44] to determine the geometrical ∆η-∆φ shape that arises from selecting jets and tracks from within our respective acceptances of |η jet | < 1.6 and |η track | < 2.4.In this technique, a mixed event distribution, ME(∆η, ∆φ), is constructed by correlating the reconstructed jet axis direction from a selected signal event to tracks from events in a minimum bias data sample.For each signal event, 40 minimum bias events are selected to have a similar vertex position (within 1 cm) and event centrality (within 2.5%) to the jet event.Mixed event correlations are corrected for tracking efficiency and misreconstruction on a per-track basis applying the same efficiency parametrization used to correct signal correlations.The distribution of such mixed event jet-track pairs N mix is denoted: This distribution ME(∆η, ∆φ) is normalized to unity at (∆η = 0, ∆φ = 0), where jets and tracks are close together and therefore have full pair acceptance.Correlations are corrected for pair acceptance effects by dividing them by the normalized mixed event distribution ME(∆η, ∆φ)/ME(0, 0).The resulting yield of associated tracks per jet is defined as: This process is illustrated in Fig. 1: the raw correlations to the leading and subleading jets in dijet events are shown on the left for the lowest p trk T bin (1-2 GeV) and 0-10% centrality range.The jet-like peak at (∆η, ∆φ) = (0, 0) is visible about both the leading and the subleading jets despite the high background levels in these most central events.An away-side peak at (∆η, ∆φ) = (0, π) is also visible in both leading and subleading jet correlations, corresponding in the leading jet correlation to the ∆η-smeared subleading jet peak, and likewise corresponding to the ∆η-smeared leading jet peak in the subleading jet correlation.The middle panel shows the shape of the pair acceptance correction determined using the mixed event technique.Finally, on the right, we present the acceptance-corrected correlations for the same p trk T bin before the subtraction of long-range correlation terms.T < 2 GeV in 0-10% central events (left), and corresponding mixed event shape ME(∆η, ∆φ) for the same centrality and p trk T bin (center).Their ratio gives the acceptance-corrected yield (right).The top row shows the correlation between leading jets (with p T,jet1 > 120 GeV) and all tracks, while the bottom row shows the correlation between subleading jets (with p T,jet2 > 50 GeV) and all tracks.
To subtract the random combinatorial backgrounds and long-range correlations (dominated by hydrodynamic flow in PbPb and momentum conservation constraints in pp events), we employ a sideband subtraction technique in which these backgrounds are approximated by the measured two-dimensional correlations in the range 1.5 < |∆η| < 3.0.Based on a CMS study that shows no appreciable variation of the elliptic flow for charged particles with p trk T > 1 GeV in the ∆η interval of ±3.0 relevant for the present analysis [45], the Fourier harmonics are assumed to be constant in ∆η.This background distribution in relative azimuthal angle (integrated over 1.5 < |∆η| < 3.0) is then fitted with a function modeling harmonic flow plus a term to capture the (Gaussian or sharper) peak at ∆φ = π associated with the (smeared) jet opposite to the jet under study: where B 0 is the overall background level; V 1 , V 2 , and V 3 are Fourier coefficients modeling harmonic flow; and A AS , α, and β are respectively the magnitude, width, and shape parameters of the away-side peak.We find that at low p trk T the long-range azimuthal sideband distributions are exhausted by the first three Fourier coefficients (V 1 , V 2 , V 3 ), while at high p trk T , V 1 and V 2 are sufficient to describe the background level within statistical uncertainties.Figure 2 illustrates the background subtraction process.The long-range contributions of the full 2D correlation (left) are estimated by the ∆φ projection (shown in the middle panel) of this correlation over the range 1.5 < |∆η| < 3.0.The fit to this ∆φ distribution is propagated uniformly in ∆η, and subtracted from the acceptance-corrected yield.The short-range correlations remaining after this background subtraction are shown on the right panel, again for p trk T = 1-2 GeV.Jet-track correlations obtained from PbPb data are compared with those obtained from the pp reference data.To ensure that the kinematic range of the jets included in this comparison is the same, correlations are reweighted on a jet-by-jet basis so that the resulting jet p T spectrum matches that of PbPb data for a given centrality class.Weighting factors are derived from the ratio of the normalized PbPb to pp jet spectra in bins of 10 GeV.The reference pp jet p T spectrum is also smeared to account for jet energy resolution differences between the PbPb and pp samples.Reference correlations in ∆η and ∆φ are then constructed and analyzed following the procedure described above for PbPb data.

Corrections and systematic uncertainties
An analysis of PYTHIA and PYTHIA+HYDJET MC simulated events is performed to evaluate and correct for the effects of two jet reconstruction biases on the measured correlated yield: a bias toward the selection of jets with harder fragmentation, and a bias toward the selection of jets that coincide with upward fluctuations in the background.The first correction addresses a jet fragmentation function (JFF) bias in which the jet energy is over-estimated for jets with hard fragmentation and under-estimated for those with soft fragmentation, resulting in a preferred selection of jets with harder fragmentation.This bias affects pp and PbPb data similarly, and results in a reduction in the charged-particle correlated yield.To correct for this effect in pp data, we compare correlations between reconstructed versus generated jets and generated particles in PYTHIA simulations, and subtract the difference (reconstructed minus generated) from data correlations.For the corresponding PbPb correction, we consider PYTHIA jets embedded into and reconstructed within a HYDJET-simulated environment, comparing correlations between generated versus reconstructed jets and the generated particles corresponding to the embedded PYTHIA hard-scattering.We note that this procedure also corrects for jet axis smearing in  reconstruction, which is found to have no significant effect on the total integral of the correlation, and to affect the correlation shapes only within ∆η < 0.2 and ∆φ < 0.2.The magnitude of this correction (relative to the total correlated yield) ranges from 3 to 6% in pp data, and from 3 to 7% in PbPb data.In PYTHIA+HYDJET, this JFF bias correction is found to be centralityindependent (and very similar to that for PYTHIA), and is applied as a single correction for all centrality bins.Maximum variations between centrality bins are used to evaluate the systematic uncertainty in this correction, which is found to be within 2% of the correlated yield at low p trk T and decreasing to zero at high p trk T .The second correction evaluates and subtracts the measured charged-particle yield resulting from the preferential selection of jets that coincide with upward fluctuations in the background as detailed in Ref. [23].The selections of inclusive and leading jets with p T > 120 GeV and subleading jets with p T > 50 GeV are sensitive to fluctuations in the background.Lower-energy jets that coincide with upward fluctuations in the background are included in the sample, while higher-energy jets that coincide with downward background fluctuations are excluded.Because the inclusive and leading jet p T spectra are both steeply falling, the inclusion in the sam-ple of a jet coinciding with an upward fluctuation in the background is much more common than the exclusion of a jet coinciding with a downward fluctuation, resulting in an excess of background tracks near the jet axis.To quantify this effect, we performed the full analysis using a sample of PYTHIA jets embedded into a HYDJET heavy ion environment, and then extracted the correlated yield (with respect to reconstructed jets) comprised of particles originating from the HYDJET background.This correction was also checked with a data-driven technique using minimum bias PbPb events to confirm that HYDJET appropriately reproduces fluctuations in the PbPb background, and the resulting upward bias in charged-particle correlated yields when these fluctuations contribute to the reconstructed jet energy.This background fluctuation bias strongly depends on event centrality and p trk T , with a magnitude of up to 24% of the corrected signal for the lowest p T tracks in the most central collisions, decreasing to within 3% for high-p T tracks and to a negligible contribution in peripheral collisions.This correlated yield due to the background fluctuation bias is subtracted to correct PbPb data, and half its magnitude is applied as a systematic uncertainty.
In addition to the systematic uncertainty associated with these two jet-reconstruction-related corrections, other sources of systematic uncertainty in this analysis include the JES determination, track reconstruction, and the procedures applied to correct for pair acceptance effects and subtract the uncorrelated and long-range backgrounds.The correlated yield uncertainty associated with the JES is assessed by varying the inclusive and leading jet p T selection threshold up and down by 3% (according to the JES uncertainty and also including differences in quark versus gluon JES [24]).The resulting maximum variations in total correlated particle yield are found to be within 3% in all cases, and we assign a 3% systematic uncertainty to account for this effect.The uncertainties of the p T -dependent tracking efficiency and misidentified track corrections are found to be within 3-4% in PbPb and pp collisions, and are independent of the centrality of the collisions.To account for the possible track reconstruction differences in data and simulation, a residual 5% uncertainty is applied based on observed variations in corrected to initial track p T and η spectra for different track quality selections [24].
We evaluate pair acceptance uncertainties by considering differences in the background levels measured separately in each of the two sideband regions of our acceptance-corrected correlations (−3.0 < ∆η < −1.5 and 1.5 < ∆η < 3.0).This results in an uncertainty within the range of 5-9%.The overall systematic uncertainty due to background subtraction is calculated by varying all fit parameters up and down by their respective uncertainties and calculating the maximum resulting differences in background level, and by considering the deviation from the "0" level after background subtraction in the sideband region 1.5 < |∆η| < 3.0.In more central events (0-10%), the background subtraction uncertainty is found to be within 2-5% for the lowest p trk T bin where the background is most significant compared to the signal level.All systematic uncertainties, evaluated as a function of p trk T and event centrality are summarized in Table 1 as fractions of the total measured yield.The range of uncertainties listed presents the variation with track transverse momentum, with larger uncertainty values corresponding to the lowest p trk T bin (1-2 GeV) for all sources.The systematic uncertainties from all seven sources are added in quadrature to obtain the total systematic uncertainty, which is quoted as a fraction of the total charged-particle yield associated with the jet under study.

Results
In this analysis, jet-track correlations are studied differentially in centrality and p trk T .Correlations are projected in ∆η and ∆φ to probe possible differences between azimuthal and pseudorapidity distributions.Figures 3 and 4 show inclusive jet correlations projected on the ∆η (over |∆φ| < 1.0) and ∆φ (over |∆η| < 1.0) axes respectively for the lowest p trk T selection.The upper panels of each figure present the centrality evolution of the correlations for inclusive jets with p T > 120 GeV, together with a reference measurement from pp data at the same collision energy shown with open symbols.To better visualize the PbPb to pp comparisons, the difference of the PbPb and pp correlation distributions is presented in the bottom panel for all centralities.Correlations are symmetrized in ∆η and ∆φ for clarity.
For the most peripheral events studied (centrality 50-100%), the PbPb correlations at low transverse momentum, 1 < p trk T < 2 GeV, show a very small excess (at most slightly larger than the uncertainties) relative to the pp reference data.This excess grows with collision centrality, with the most significant excess present in the most central collisions.The shape of this excess in the low-p trk T per-jet particle yields is found to be similar in the ∆η and ∆φ distributions, and in both dimensions exhibits a Gaussian-like shape that extends to large relative angles ∆η ≈ 1 and ∆φ ≈ 1.We note that these results are consistent with previous CMS studies of jet-shape modifications [22] and fragmentation functions [23] within the previously studied small ∆R < 0.3 region, while extending measurements to individually study ∆η and ∆φ distributions over the full range ∆η and ∆φ < 1.5.
The next two figures present the results of the jet-track correlation measurements for dijets with leading jet p T > 120 GeV and subleading jet p T > 50 GeV, obeying the back-to-back angular selection criteria previously described.Figure 5 presents the projection of jet-track correlations measured for charged tracks with p trk T between 1 and 2 GeV on the ∆η axis for the leading (upper panel) and subleading (middle panel) jets, while Fig. 6 shows the corresponding projections on the ∆φ axis.Again pp data are included for comparison, and for the most peripheral (50-100% central) PbPb events the correlations are similar to the pp reference for the leading jets, and differ only slightly for the subleading jets.As in the case of inclusive jets, differences of correlations between pp and PbPb collisions gradually increase from peripheral to central collisions, and are most pronounced in the 0-10% central events for both leading and subleading jets.We note that there is little difference between the leading and inclusive jet correlated-yield distributions, indicating that the requirement that leading jets have the highest p T in the event does not significantly bias the selection of jets with p T > 120 GeV.
For this lowest p trk T bin shown, we observe that (as for the inclusive jet selection) the excess of correlated yield extends significantly beyond the typical jet reconstruction radius for both leading and subleading jets.The soft excess is more pronounced on the more "quenched" subleading side, but is also present on the leading side.This indicates that leading jets, although surface-biased toward shorter path-lengths through the medium, also experience quenching in   central PbPb collisions.To better illustrate both subleading and leading modifications, the last row of Figs. 5 and 6 shows the differences (PbPb minus pp) of the correlations in the two upper panels.
To quantify the total per-jet excess yield observed in the PbPb data with respect to the pp reference, we plot the integrals of the excess yields (PbPb minus pp) as a function of p trk T and collision centrality in Fig. 7.As the figure shows, in both leading and subleading jets, the excess yield diminishes for higher momentum tracks until the yield becomes similar to the pp reference for the highest p trk T bin of 4-8 GeV.As seen in previous figures, central collisions exhibit the largest low-p trk T excesses.This demonstrates the expected trend corresponding to quenching of both the leading and the subleading jets, as energy from particles with higher p trk T is redistributed into particles with lower p trk T via interactions with the medium.In order to characterize the angular widths of the charged-particle distributions in ∆η and ∆φ, we fit the measured correlations with a double Gaussian function (which was found to best describe the overall correlation shapes).The width is defined as the region around zero in |∆η| or |∆φ| that contains 67% of the total correlated yield.Width uncertainties are calculated by repeating the measurement for the ∆η and ∆φ distributions varied by their respective systematic uncertainties, which are treated as fully correlated for the purposes of this determination.Widths for leading and subleading jet correlations in ∆η and ∆φ are presented as a function of p trk T in Figs.8-11.Distributions of low-p T tracks correlated with either of the two jets are found to be significantly broader in central PbPb events compared to those in pp data in both ∆η and ∆φ dimensions.This broadening is greatest for the low-p T tracks and in the most central events, and diminishes quickly with increasing track momenta.Above 4 GeV, the widths measured in PbPb and pp events are the same within the systematic uncertainties.We note that the width of the PbPb minus pp excess yield is similar for leading and subleading jets.In pp data, however, the peak associated with the subleading jet is softer and broader than the peak associated with the leading jet.There is therefore a larger difference in peak width when comparing PbPb leading jet peaks to the narrow pp leading jet peaks (Figs.8-9), and a smaller difference when comparing PbPb subleading jet peaks to the broader pp subleading jet peaks (Figs.[10][11].

Summary
In this analysis, jet-track correlations have been studied as a function of ∆η and ∆φ with respect to the jet axis in PbPb and pp collisions at √ s NN = 2.76 TeV.Two-dimensional angular correlations have been considered for charged particles with p trk T > 1 GeV as a function of p trk T and collision centrality for two jet selections.A sample of inclusive jets above the jet momentum threshold of 120 GeV was studied, as well as a sample of dijet events selected to include a leading jet with p T > 120 GeV and a subleading jet with p T > 50 GeV.In all cases, an excess of soft particle yields was observed in central PbPb collisions with respect to pp reference data, similar for inclusive and leading jet samples and larger for the (more-quenched) subleading jet sample.The low-p trk T (1-3 GeV) excess-yield distributions were studied individually and, in both ∆η and ∆φ, they exhibit similar Gaussian-like distributions out to large relative angles (∆η ≈ 1 and ∆φ ≈ 1) from the jet axis.The excess was found to be largest at the lowest p trk T (1-2 GeV) in the most central (0-10%) PbPb data, and to decrease gradually with centrality.For peripheral (50-100%) PbPb collisions, correlated low-p trk T particle yields are only slightly larger than those for the pp reference.The excess also gradually decreases with increasing p trk T until yields of particles with p trk T > 4 GeV are similar to pp reference data, consistent with the results of a previous CMS jet quenching study.This new correlation analysis provides a comprehensive evaluation of medium effects on jet properties, extending information about jet shapes to

Figure 1 :
Figure1: Jet-track correlation signal shape S(∆η, ∆φ) for tracks with 1 < p trk T < 2 GeV in 0-10% central events (left), and corresponding mixed event shape ME(∆η, ∆φ) for the same centrality and p trk T bin (center).Their ratio gives the acceptance-corrected yield (right).The top row shows the correlation between leading jets (with p T,jet1 > 120 GeV) and all tracks, while the bottom row shows the correlation between subleading jets (with p T,jet2 > 50 GeV) and all tracks.

Figure 2 :
Figure 2: Acceptance-corrected 2D jet-track correlation yield (left) is projected over the range 1.5 < |∆η| < 3.0, producing a 1D background distribution (center).The fit to this distribution (indicated with a red dark line) is subtracted from the total yield to obtain the 2D backgroundsubtracted yield shown on the right (for tracks with 1< p trk T < 2 GeV).The black dashed line shows the background level and Fourier flow harmonic components of the fit only, excluding the away-side peak.Yellow lines in the B(∆φ) plot (middle panel) indicate the systematic uncertainty assigned to the background subtraction.

Figure 3 :
Figure 3: Symmetrized ∆η distributions (projected over |∆φ| < 1) of background-subtracted particle yields correlated to PbPb and pp inclusive jets with p T > 120 GeV are shown in the top panels for tracks with 1 < p trk T < 2 GeV.The difference in PbPb and pp per-jet yields is shown in the bottom panels.The total systematic uncertainties are shown as shaded boxes, and statistical uncertainties are shown as vertical bars (often smaller than the symbol size).

Figure 4 :
Figure 4: Symmetrized ∆φ distributions (projected over |∆η| < 1) of background-subtracted particle yields correlated to PbPb and pp inclusive jets with p T > 120 GeV are shown in the top panels for tracks with 1 < p trk T < 2 GeV.The difference in PbPb and pp per-jet yields is shown in the bottom panels.The total systematic uncertainties are shown as shaded boxes, and statistical uncertainties are shown as vertical bars (often smaller than the symbol size).

Figure 5 :
Figure 5: The top panels show the ∆η distributions (projected over |∆φ| < 1) of chargedparticle background-subtracted yields correlated to PbPb and pp leading jets with p T,jet1 > 120 GeV.The middle panels show the same distributions for subleading jets with p T,jet2 > 50 GeV, and the bottom panels show the difference PbPb minus pp for both leading and subleading jets.The total systematic uncertainties are shown as shaded boxes, and statistical uncertainties are shown as vertical bars (often smaller than the symbol size).

Figure 6 :
Figure 6: The top panels show the ∆φ distributions (projected over |∆η < 1) of charged-particle background-subtracted yields correlated to PbPb and pp leading jets with p T,jet1 > 120 GeV.The middle panels show the same distributions for subleading jets with p T,jet2 > 50 GeV, and the bottom panels show the difference PbPb minus pp for both leading and subleading jets.The total systematic uncertainties are shown as shaded boxes, and statistical uncertainties are shown as vertical bars (often smaller than the symbol size).

Figure 7 :
Figure 7: Total excess correlated yield observed in the PbPb data with respect to the reference measured in pp collisions, shown as a function of track p T in four different centrality intervals (0-10%, 10-30%, 30-50%, 50-100%) for both leading jets with p T,jet1 >120 GeV and subleading jets with p T,jet2 > 50 GeV.The total systematic uncertainties are shown as shaded boxes, and statistical uncertainties are shown as vertical bars (often smaller than the symbol size).

Figure 8 :
Figure 8: Comparison of the widths in PbPb and pp of the ∆η charged-particle distributions correlated to leading jets with p T,jet1 > 120 GeV, as a function of p trk T .The bottom row shows the difference of the widths in PbPb and pp data.The shaded band corresponds to systematic uncertainty, and statistical uncertainties are smaller than symbol size.

Figure 9 :
Figure 9: Comparison of the widths in PbPb and pp of the ∆φ charged-particle distributions correlated to leading jets with p T,jet1 > 120 GeV, as a function of p trk T .The bottom row shows the difference of the widths in PbPb and pp data.The shaded band corresponds to systematic uncertainty, and statistical uncertainties are smaller than symbol size.

Figure 10 :
Figure 10: Comparison of the widths in PbPb and pp of the ∆η charged-particle distributions correlated to leading jets with p T,jet2 > 50 GeV, as a function of p trk T .The bottom row shows the difference of the widths in PbPb and pp data.The shaded band corresponds to systematic uncertainty, and statistical uncertainties are smaller than symbol size.

Figure 11 :
Figure 11: Comparison of the widths in PbPb and pp of the ∆φ charged-particle distributions correlated to leading jets with p T,jet2 > 50 GeV, as a function of p trk T .The bottom row shows the difference of the widths in PbPb and pp data.The shaded band corresponds to systematic uncertainty, and statistical uncertainties are smaller than symbol size.

Table 1 :
Systematic uncertainties in the measurement of the jet-track correlations in PbPb and pp collisions, as percentage of the total measured correlated yield.The numbers presented in this table summarize the range of values of systematic uncertainty (as a function of p trk