Jet properties in PbPb and pp collisions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV

Modifications of the properties of jets in PbPb collisions, relative to those in pp collisions, are studied at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV via correlations of charged particles with the jet axis in relative pseudorapidity ($\Delta \eta$), relative azimuth ($\Delta \phi$), and relative angular distance from the jet axis $\Delta \mathrm{r} = \sqrt{{(\Delta\eta)^{2}+(\Delta\phi)^{2}}}$. This analysis uses data collected with the CMS detector at the LHC, corresponding to integrated luminosities of 404 $\mu$b$^{-1}$ and 27.4 pb$^{-1}$ for PbPb and pp collisions, respectively. Charged particle number densities, jet fragmentation functions, and jet shapes are presented as a function of PbPb collision centrality and charged-particle track transverse momentum, providing a differential description of jet modifications due to interactions with the quark-gluon plasma.


Introduction
Jets can be used as proxies for partons produced in the initial hard scatterings in heavy ion collisions to probe the properties of the quark-gluon plasma (QGP), a new state of matter characterized by an increase in the color degrees of freedom. One of the well-established properties of the QGP is its high opacity to such penetrating probes, resulting in significant energy loss of partons traversing the hot nuclear matter. Parton energy loss manifests itself in a variety of experimental observables, including suppression of high transverse momentum (p T ) hadrons and jets, as well as modifications of the properties of parton showers. These phenomena, collectively referred to as jet quenching [1], were first observed at the BNL RHIC [2, 3], and subsequently at the CERN LHC [4][5][6][7]. The LHC experiments have previously demonstrated that the medium also affects the structure of a jet, as observed from measurements of the jet fragmentation pattern [8,9] and the distribution of charged-particle transverse momenta (p trk T ) as a function of the relative angular distance ∆r from the jet axis [10], where lowercase r is used explicitly to avoid conflict with the jet clustering distance parameter R [11]. The distance ∆r is given by ∆r = √ (∆η) 2 + (∆φ) 2 , where ∆η and ∆φ denote the relative pseudorapidity and azimuthal angle (in radians) with respect to the jet axis, respectively. These modifications extend to large values of ∆η and ∆φ [12][13][14]. Various theoretical models have since attempted to account for these modifications [15][16][17][18][19], and while most models reproduce the modification effects close to the jet axis, the large modifications far from the jet axis (∆r > 0.5) are not yet understood.
This paper describes modifications to jet structure in PbPb collisions at a nucleon-nucleon center of mass energy √ s NN = 5.02 TeV relative to pp collisions at the same energy, extending previous results based on 2.76 TeV data [10,13]. At the higher collision energy, an increase in the magnitude of jet quenching is expected because of the greater medium density and temperature [20], potentially increasing the size of the modification effects. The data were collected by the CMS detector at the LHC and correspond to integrated luminosities of 404 µb −1 and 27.4 pb −1 for PbPb and pp collisions, respectively. The distributions of charged-particle tracks with respect to the jet axis are studied as a function of ∆η, ∆φ, and ∆r. The jet shapes ρ(∆r), defined as the distribution of particle yields in ∆r weighted by p trk T , are also examined. The results are presented differentially in p trk T and as a function of the overlap of the colliding Pb nuclei (centrality), with head-on collisions defined as most central. Compared to the size of the data samples in Ref. [13], the present study uses a much larger data set and, hence, has a greater statistical precision. This also allows the measurements to be extended to larger distances with respect to the jet axis.

The CMS detector
The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. 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 barrel and endcap sections. Two hadronic forward (HF) steel and quartz-fiber calorimeters complement the barrel and endcap detectors, extending the calorimeter from the range |η| < 3.0 provided by the barrel and endcap out to |η| < 5.2. The scalar p T sum of calorimeter towers in the HF region (4.0 < |η| < 5.2) is used to define the event centrality in PbPb events and to divide the event sample into centrality classes, each representing a percentage of the total nucleus-nucleus hadronic interaction cross section. A detailed description of the centrality determination can be found in Ref. [6].

Event selection and simulated event samples
The pp and PbPb data are selected with a calorimeter-based trigger that uses the anti-k T jet clustering algorithm with distance parameter of R = 0.4 [24]. The trigger requires events to contain at least one jet with p T > 80 GeV. This trigger is fully efficient for events containing jets with reconstructed p T > 90 GeV. For both PbPb and pp collisions, the data selected by this trigger are referred to as "jet-triggered", and are used to study the jet-related particle yields. While jet-triggered samples are used in pp collisions to correct for the limited jet and track acceptance via an event mixing technique described in Section 5, an additional data sample collected with a minimum-bias trigger [25] is used for this correction in PbPb collisions in order to properly capture the long-range correlated particle yields in PbPb data. To reduce contamination from noncollision events, including calorimeter noise and beam-gas collisions, vertex and noise filters are applied to both the pp and PbPb data as described in previous analyses [5,6]. The filters include the requirements that events contain at least 3 GeV of energy in each of the two HF calorimeters and that a primary vertex with at least two tracks be present within 15 cm of the center of the nominal interaction region along the beam axis (|v z | < 15 cm).
Monte Carlo (MC) simulated event samples are used to evaluate the performance of the event reconstruction, particularly the track reconstruction efficiency and the jet energy response and resolution. The MC samples use the PYTHIA (version 6.424, tune Z2 [26,27]) event generator to describe the hard scattering, parton showering, and hadronization of the partons. The GEANT4 [28] toolkit is used to simulate the CMS detector response. To account for the soft underlying PbPb event component, the hard PYTHIA interactions are embedded into simulated minimum-bias PbPb events, produced with the HYDJET 1.383 [29] event generator. We refer to this latter sample as PYTHIA +HYDJET.
In PbPb collisions, jets are produced more frequently in central events than in noncentral events because of the large number of binary collisions per nuclear interaction. Since the HYDJET event generator simulates minimum-bias PbPb collisions only, a centrality-based reweighting is applied to the PYTHIA +HYDJET sample in order to match the centrality distribution of the jettriggered PbPb data. An additional reweighting procedure is performed to match the simulated v z distributions to data for both the pp and PbPb samples.

Jet and track reconstruction
The jet reconstruction in PbPb and pp events is performed with the anti-k T jet algorithm with a distance parameter of R = 0.4, as implemented in the FASTJET framework [11], with individually calibrated calorimeter towers as input. Only calorimeter information is used in order to avoid biases due to an interplay between track reconstruction efficiency and the jet energy scale. In PbPb collisions, the contributions of the underlying event are subtracted using a variant of the iterative "noise/pedestal subtraction" technique described in Ref. [30]. Following the subtraction, jets are calibrated such that the calorimeter response is uniform as a function of jet p T and η. A reweighting procedure based on the number of combined calorimeter towers and associated charged-particle tracks with p T > 2 GeV is then applied to correct for the variation in detector response with the total number of jet constituents. This latter calibration corrects for a difference in the simulated calorimetric jet energy response between quark and gluon jets and reduces the difference in this response between the two jet flavors, as a function of jet p T , from 10% to around 3%. After reconstruction and offline jet energy calibration, jets are required to have p T > 120 GeV and |η| < 1.6. Within this selection, it is possible for multiple jets to be selected from the same event. Roughly 25% of pp events contain multiple jets that satisfy all kinematic selection criteria and no additional selections are made to distinguish between leading and subleading jets in this measurement.
For pp data and simulations, charged-particle tracks are reconstructed using an iterative tracking method [22] that allows the reconstruction of charged-particle tracks within |η| < 2.4 down to p T = 0.1 GeV. For the PbPb case, an alternative iterative charged-particle reconstruction procedure is employed because of the large track multiplicity in such collisions, as discussed in earlier heavy ion analyses [10,31]. This reconstruction algorithm is based on hit information from both the pixel and silicon strip subdetectors and is capable of reconstructing charged-particle tracks down to p T = 0.4 GeV. The tracking efficiency in pp collisions ranges from approximately 80% at p T = 0.5 GeV to 90% or higher for p T > 10 GeV. Track reconstruction is more difficult in the heavy ion environment because of the large track multiplicity and so the tracking efficiency ranges from approximately 30% at p T = 0.5 GeV to about 70% at p T = 10 GeV [32]. Corrections for the tracking efficiency and related effects are derived as a function of p trk T , η, and φ using PYTHIA simulations, and, additionally for PbPb events, as a function of centrality using PYTHIA +HYDJET.

Jet-track angular correlations
Correlations between reconstructed jets and charged-particle tracks are studied by forming a two-dimensional array of the ∆η and ∆φ values of the tracks relative to the jet axis. Events in PbPb collisions are divided into four centrality intervals based on the fraction of the total recorded energy that is collected within the HF calorimeter, given by 0-10% (most central, corresponding to the largest overlap of the colliding nuclei), 10-30%, 30-50%, and 50-100% (most peripheral), and into eight bins of p trk T bounded by the values 0.7, 1, 2, 3, 4, 8, 16, 20, and 300 GeV. After classification, the jet-track correlation yields are normalized by the number of jets in the sample and corrected for tracking efficiency on a per-track basis. The normalization creates a per-jet averaged ∆η-∆φ distribution of charged particles about the jet axis for each p trk T and centrality interval. For the p trk T measurements, the two-dimensional correlations are weighted by p trk T on a per-track basis, producing a per-jet averaged ∆η-∆φ distribution of p trk T about the jet axis.
Following construction of the two-dimensional correlations described above, the remaining analysis procedure consists of the following three steps: first, an acceptance correction is derived using a mixed-event method, which accounts for the effects of the limited detector acceptance. Second, the backgrounds from tracks unrelated to jets are subtracted using a region in ∆η far from the jet axis. Third, simulation-derived corrections are applied to account for jet reconstruction biases. Each of these steps is described in detail below.
As a consequence of the shape of the jet and track η distributions and the different acceptance requirements for jets (|η| < 1.6) and tracks (|η| < 2.4), the jet-track correlation yield decreases with increasing ∆η. This decrease is due to a combination of the η-dependence of jet and particle production, as well as the limited detector acceptance. Jet-track pairs, therefore, have a higher probability to be reconstructed at central values of jet and/or track η than in the forward or backward regions. To correct for this pair-acceptance effect, a mixed-event distribution is constructed by creating jet-track pairs using jets from the jet-triggered event sample and tracks from a sample of minimum-bias events, matched in the vertex position along the beam axis (within 1 cm) and collision centrality (within 2.5%). This procedure is analogous to that employed for two-particle correlations in Refs. [33][34][35] and for jet-track correlations in Refs. [13,14,36]. In the following, N jets denotes the total number of jets satisfying the selection criteria, either from pp or PbPb collisions. The per-jet associated yield is corrected for the acceptance effects via the following relation: where the signal pair distribution, S(∆η, ∆φ), represents the yield of jet-track pairs from the same event, normalized by N jets , and the mixed-event pair distribution ME(∆η, ∆φ) is The ratio ME(0, 0)/ME(∆η, ∆φ) is the normalized correction factor; ME(0, 0) is the mixedevent yield for jet-track pairs that are approximately collinear. These collinear pairs have the maximum pair acceptance.
The acceptance-corrected distributions resulting from Eq. (1) exhibit a Gaussian-like peak confined to a fairly narrow ∆η, ∆φ range on top of a large background from unrelated jet-track pairs and pairs connected through long-range correlations, e.g., azimuthal anisotropies. To model this background, the ∆φ distribution averaged over 1.5 < |∆η| < 2.5 is used to estimate the ∆φ-dependence of the combinatoric contribution to the correlations over the entire |∆η| < 4.0 region and is subtracted from the acceptance-corrected yields of Eq.
(1). The use of a narrow range in ∆η to model the background automatically accounts for the magnitude of the azimuthal anisotropy contribution, without the need for an explicit measurement.
Finally, simulation-based corrections are applied to account for jet position resolution and two biases in the jet reconstruction: a bias toward selecting jets with a harder constituent p T spectrum (affecting PbPb and pp events similarly), and a bias toward selecting jets that are affected by upward fluctuations in the p T values of the objects in the soft underlying event (relevant for PbPb events only). Jets with a harder constituent p T spectrum are more likely to be successfully reconstructed than jets with a softer constituent p T spectrum because the calorimeter response does not scale linearly with incident particle energy, resulting in a bias toward the selection of jets with fewer associated tracks. This bias is reduced, but not eliminated, by applying the jet energy corrections mentioned above based on the number of jet constituents. A residual correction for this bias is derived following the method described in Refs. [12][13][14], by comparing per-jet yields of generated particles correlated to reconstructed jets relative to those correlated to generated jets. For pp events, this correction is derived using the PYTHIA simulation. For PbPb events, we consider only generated particles from the embedded PYTHIA hard process, excluding particles from the underlying event.
For PbPb events, there is an additional jet reconstruction bias toward the selection of jets that are produced in the vicinity of upward fluctuations in the underlying event. Since the jet p T spectrum falls steeply, more jets on upward fluctuations are included in the sample than jets on downward fluctuations are excluded. To estimate and then account for this bias, we follow a similar procedure to that outlined in Refs. [13,14,37]. We consider correlations in the PYTHIA +HYDJET sample between reconstructed jets and generated particles from only the HYDJET underlying event, excluding particles from the embedded hard process. To avoid propagating low-p T HYDJET fluctuations to the data, the per-jet excess yields from the HYDJET underlying event are fit with Gaussian functions in both ∆η and ∆φ and applied as a correction to the PbPb data.

Systematic uncertainties
A number of sources of systematic uncertainty are considered, including effects from the tracking efficiency, acceptance corrections, background subtraction, and jet reconstruction. Where relevant, systematic uncertainties related to tracking and jet reconstruction are determined by comparing properties of reconstructed PYTHIA and PYTHIA +HYDJET events to their generated counterparts. Jet reconstruction related sources of systematic uncertainty include the two biases in jet reconstruction discussed in Section 5 and the jet energy scale (JES).
To evaluate the uncertainty related to the jet reconstruction biases and the jet position resolution, three variations of the analysis are performed to gauge differences of the jet energy calorimetric response in data relative to the MC simulation. First, the collision reaction plane dependence of the jet reconstruction performance is tested by determining the corrections independently for in-plane and out-of-plane jets, where the reaction plane is the plane containing the beam axis and the projection onto the transverse plane of the line connecting the centers of the colliding nuclei. The difference in the corrections and hence the resulting uncertainty is found to be negligible. Second, the uncertainties associated with the evaluation of the magnitude of both jet reconstruction biases are obtained by varying the corresponding acceptance correction and background subtraction within their statistical uncertainties. For the jet constituent hardness bias, the resulting uncertainty in the correction factor is found to be negligible. Conversely, for the correction accounting for upward underlying event fluctuations, the uncertainty is found to vary between 1 and 18% of the total correction, strongly depending on p trk T and the centrality. Third, the dependence on the relative numbers of quark and gluon jets is studied. By fitting distributions of the quark and gluon jet constituent multiplicities in simulation, we obtain templates that are used in a fit to estimate the fraction of quark and gluon jets in data. We observe that the quark jet fraction in data is centrality dependent, ranging between 49 and 56%, while the quark jet fraction in simulation is constant at approximately 42%. The jet reconstruction bias corrections are reweighted in simulation to correspond to the measured ratio of quark to gluon jets. The reweighting has a negligible impact on the correction account-ing for upward fluctuations in the underlying event, but affects the corrections that account for the bias toward jets with a harder constituent p T spectrum on the order of 10%, resulting in an uncertainty of 2% in the correlated yields. In addition to the uncertainties evaluated from these three variations of the analysis, an uncertainty associated with the Gaussian fitting procedure described at the end of Section 5 is given by the uncertainty in this fit, though this uncertainty is negligible compared to the overall reconstruction bias uncertainty.
The JES is subject to three sources of uncertainty, all of roughly equal magnitude. All are related to potential differences between simulation and data. First, differences in the relative fraction of quark and gluon jets can affect the overall energy scale of the jet sample. Second, jet quenching effects in PbPb events are not simulated in HYDJET. Third, residual differences in the calorimeter response exist between data and simulation. All three sources are conservatively accounted for by a 5% variation in the jet energy scale. Template fits are used to estimate the quark and gluon jet fractions in pp data, absent of quenching effects, and are found to be consistent to within 5%. Studies of jet quenching result in an estimate of a 7-10 GeV shift in jet energy, roughly half of which remains in the jet cone. Finally, residual differences in calorimeter response are estimated using gamma-jet correlations [12]. These studies yield observed discrepancies of up to 5% in the calorimeter response between photons and jets. To evaluate the effect of all these uncertainties on the results, we vary the jet p T threshold of 120 GeV up and down by 6 GeV, corresponding to a shift of 5%. The resulting uncertainty in the correlated yields is found to be around 2%, because the jet multiplicities vary only slowly as a function of jet p T .
The uncertainty related to the tracking efficiency is estimated from simulation by taking the ratio between the corrected reconstructed yields and the corresponding generated yields. This ratio is referred to as the "closure" and its deviation from unity defines the systematic uncertainty. The systematic uncertainty in the tracking correction is found to be 1% in both PbPb and pp events, noting that the PbPb closure is derived using simple two-dimensional efficiency tables in η and φ, while the pp closure uses a multidimensional parameterization of the tracking efficiency. An uncertainty is also evaluated to account for possible differences in track reconstruction between data and simulation, including the erroneous reconstruction of tracks, and is found to be 5% (4%) in PbPb (pp) events.
The uncertainty arising from long-range ∆η-dependent asymmetries of the mixed-event acceptance correction is estimated by considering the sideband asymmetry of the correlated yield after dividing by the mixed-event background. Ideally, the mixed-event acceptance corrected signal is composed of a correlated jet-track yield sitting on top of an uncorrelated background. Far from the jet axis, the uncorrelated background should dominate the mixed event corrected signal and, thus, the signal is expected to be uniform for |∆η| > 1.5. To probe for discrepancies from this ideal case, each sideband region of the final ∆η distribution (−2.5 < ∆η < −1.5 and 1.5 < ∆η < 2.5) is independently fit with a constant to estimate the nonuniformity of the mixed event correction, with respect to any residual ∆η asymmetry. The difference between these two fits is assigned as a systematic uncertainty, and is found to be about 5% for the lowest p trk T bin, where such effects are largest.
Uncertainties associated with the background subtraction are evaluated by considering the average point-to-point difference between two sideband regions (1.5 < |∆η| < 2.0 and 2.0 < |∆η| < 2.5) following the background subtraction. In central events (0-10%), the background subtraction uncertainty is found to be roughly 4% for the lowest p trk T bin, where the ratio of signal to background events is lowest, and decreases as a function of centrality and p trk T . The systematic uncertainties from all sources are added in quadrature. Table 1 lists the ranges of the estimated contributions from the individual sources.  T < 3 GeV) in the PbPb data, relative to the pp data, that becomes more pronounced with increasing centrality. This low-p trk T excess is relatively symmetric in ∆η and ∆φ and remains significant even for large ∆η and ∆φ. The particle yield excess in the PbPb data decreases with increasing p trk T , such that no significant enhancement is observed for p trk T > 3 GeV. Figure 3 shows the jet-track correlations as a function of ∆r. As for the ∆η and ∆φ distributions, an enhancement for p trk T > 3 GeV is seen in the PbPb data relative to the pp data, which increases with increasing centrality. The change of slope in the figures at ∆r = 0.4 is due to the nature of the anti-k T algorithm, where low-p T particles just inside the jet radius parameter are generally clustered within the jet and particles just outside the radius are not.
To further illustrate the p trk T dependence of the results, Fig. 4 (top row) shows the chargedparticle track yields in PbPb and pp events as a function of p trk T , integrated over ∆η and ∆φ. The excess observed in PbPb events is seen to decrease smoothly with increasing p trk T , in each centrality interval. For p trk T > 3 GeV, the PbPb results are consistent with those from the pp collisions. Figure 4 also shows a comparison of the results at 5.02 TeV with those previously obtained at 2.76 TeV [13] in order to show the dependence of the excess on the center-of-mass energy. The results from the two center-of-mass energies are consistent within one standard deviation.
Measurements of the jet shapes ρ(∆r) are obtained by examining the distribution of chargedparticle tracks in annular rings of width ∆r = 0.05 around the jet axis, with each particle weighted by its p trk T value. In contrast to Figs. 1-4, tracks with p trk T > 20 GeV are included in these distributions as they make a significant contribution to the jet momentum, even though their rate is small. The resulting transverse momentum profile P(∆r) of the jet is defined as: where ∆r a and ∆r b define the annular edges of ∆r, and δr = ∆r b − ∆r a . This profile is normalized to unity within ∆r = 1 to produce ρ(∆r): ( The plots in the top left and middle row of Fig. 5 show the P(∆r) distribution for pp and PbPb events, respectively. The bottom row shows the ratio of the PbPb to the pp data for three different ranges in p trk T , namely p trk T > 0.7, 2.0, and 4.0 GeV. These results demonstrate a large excess of soft particles in PbPb events relative to pp events at intermediate to large angles from the jet axis, compensated for by a relative depletion at all angles of tracks at high-p trk T . For charged-particle tracks with p trk T > 0.7 GeV, the difference between the PbPb and pp data  reaches nearly a factor of two for large ∆r in central collisions, as seen from the bottom right plot in Fig. 5. This behavior is possibly due to a combination of jet quenching in the medium and the so called "backreaction" response of the medium to the jet, where the jet induces spallationtype effects in the soft underlying event, increasing low-p trk T contributions over a wide range in ∆r. Out of a number of theoretical models at 2.76 TeV [19,[38][39][40] only those that include such effects in combination with other processes, like jet quenching and jet broadening, are able to reproduce the large low-p trk T excesses at very large ∆r observed in the momentum profiles. The analogous results for ρ(∆r) are shown in Fig. 6. The ratios of the PbPb to the pp data, presented in the bottom row of Fig. 6, are shown for the inclusive range 0.7 < p trk T < 300 GeV. A redistribution of energy from small to large angles relative to the jet axis is evident from these data, as seen (for example) from the dip and then rise in the PbPb/pp ratio distributions as ∆r increases. The effect is more pronounced in central events. From Fig. 6, we conclude that this large enhancement at large ∆r is dominated by the contribution from low-p trk T tracks, has pronounced centrality dependence, and, as shown in Fig. 4, is slightly larger at higher centerof-mass energy.   Note that for the 2.76 TeV data, shown in Fig. 2 of Ref.
[14], the PbPb/pp jet shape ratios for leading jets at large ∆r are larger than those shown for inclusive jets from 5.02 GeV in Fig. 6 of this manuscript. This is mainly because the jet shape in the 2.76 TeV pp reference data falls more steeply as a function of ∆r than at 5.02 TeV, resulting in a larger value for the 2.76 TeV ratio even though the magnitude of the modification effects is similar at the two energies. This difference in the shape of the two pp reference samples is also present in the PYTHIA simulation, where the effect is due to differences in the relative fraction of quark and gluon jets at the two energies.

Summary
In this paper, measurements are presented of the modifications to charged-particle track yields and jet shapes in PbPb collisions with respect to pp collisions using data collected with the CMS detector at the LHC at a nucleon-nucleon center-of-mass energy of √ s NN = 5.02 TeV. The correlations of charged particles having transverse momentum p trk T > 0.7 GeV and pseudorapidity |η| < 2.4 with the jet axis of jets having p T > 120 GeV and |η| < 1.6 are studied. Charged particle yields associated with jets are shown as a function of relative angular distance ∆r = √ (∆η) 2 + (∆φ) 2 from the jet axis, as well as individually in ∆η and ∆φ. In these studies, a strong enhancement of tracks with p trk T < 3 GeV, extending to large angles, is found in PbPb collisions with respect to pp collisions. This low p T excess remains correlated with the jet axis, but the distribution is broader in PbPb than that observed in the corresponding p T bin of pp collisions, which could indicate additional in-medium gluon radiation and/or a medium backreaction, i.e., a wake-like response of the QGP to the propagating parton. For PbPb events with centrality 0-10%, the correlated yield is increased relative to pp collisions at low p T by up to a factor of two.
In addition to these charged-particle yields, we examine the jet transverse momentum profile P(∆r) and the jet shape ρ(∆r) variables, defined using the distribution of charged-particle tracks in annular rings around the jet axis, with each particle weighted by its p trk T value. A redistribution of energy from small to large angles from the jet axis is observed for PbPb relative to pp events, with the most pronounced effect seen for central collisions. The energy flow within the jet is modified by shifting the momentum away from the jet axis out to large relative angular distances such that in central PbPb collisions, the ratio of the jet shapes in PbPb to pp collisions approaches a value of two. These measurements provide a comprehensive picture of the modifications of the parton shower evolution and the quark-gluon plasma response to the propagating jets in 5 TeV PbPb collisions.

Acknowledgments
We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centres and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: BMWFW and FWF (Aus