Edinburgh Research Explorer Measurement of inclusive jet charged-particle fragmentation functions in Pb+Pb collisions at $\sqrt{s_{NN}}=2.76 TeV with the ATLAS detector

Measurements of charged-particle fragmentation functions of jets produced in ultra-relativistic nuclear collisions can provide insight into the modiﬁcation of parton showers in the hot, dense medium created in the collisions. ATLAS has measured jets in √ s NN = 2 . 76 TeV Pb + Pb collisions at the LHC using a data set recorded in 2011 with an integrated luminosity of 0.14 nb − 1 . Jets were reconstructed using the anti- k t algorithm with distance parameter values R = 0 . 2 , 0 . 3 , and 0 . 4. Distributions of charged-particle transverse momentum and longitudinal momentum fraction are reported for seven bins in collision centrality for R = 0 . 4 jets with p jetT > 100 GeV. Commensurate minimum p T values are used for the other radii. Ratios of fragment distributions in each centrality bin to those measured in the most peripheral bin are presented. These ratios show a reduction of fragment yield in central collisions relative to peripheral collisions at intermediate z values, 0 . 04 (cid:2) z (cid:2) 0 . 2, and an enhancement in fragment yield for z (cid:2) 0 . 04. A smaller, less signiﬁcant enhancement is observed at large z and large p T in central collisions.


Introduction
Collisions between lead nuclei at the LHC are thought to produce a quark-gluon plasma (QGP), a form of strongly interacting matter in which quarks and gluons become locally deconfined. One predicted consequence of QGP formation is the "quenching" of jets generated in hard-scattering processes during the initial stages of the nuclear collisions [1]. Jet quenching refers, collectively, to a set of possible modifications of parton showers by the QGP through interactions of the constituents of the shower with the colour charges in the plasma [2,3]. In particular, quarks and gluons in the shower may be elastically or inelastically scattered resulting in both deflection and energy loss of the constituents of the shower. The deflection and the extra radiation associated with inelastic processes may broaden the parton shower and eject partons out of an experimental jet cone [4][5][6][7][8][9]. As a result, jet quenching can potentially both soften the spectrum of the momentum of hadrons inside the jet and reduce the total energy of the reconstructed jet. A complete characterization of the effects of jet quenching therefore requires measurements of both the single-jet suppression and the jet fragment distributions.
Observations of modified dijet asymmetry distributions [10][11][12], modified balance-jet transverse momentum (p T ) distributions in γ + jet events [13], and suppressed inclusive jet yield in Pb+Pb collisions at the LHC [14,15] are consistent with theoretical calcu-E-mail address: atlas.publications@cern.ch. lations of jet quenching. However, it has been argued that those measurements do not sufficiently discriminate between calculations that make different assumptions regarding the relative importance of the contributions described above [16]. Based on the above arguments, theoretical analyses are incomplete without experimental constraints on the theoretical description of jet fragment distributions.
This Letter presents measurements of charged-particle jet fragmentation functions in √ s NN = 2.76 TeV Pb+Pb collisions using 0.14 nb −1 of data recorded in 2011. The jets used in the measurements were reconstructed with the anti-k t [17] algorithm using distance parameter values R = 0.2, 0.3, and 0.4. Results are presented for the charged-particle transverse momentum ( p ch T ) and longitudinal momentum fraction (z ≡ p ch particle and jet directions. 1 The p jet T thresholds for the three R values were chosen to match the R-dependence of the measured transverse momentum of a typical jet. For simplicity, the terms "fragmentation functions" are used to describe the distributions defined in Eq. (2) with the understanding that D(z) is different from a theoretical fragmentation function, D(z, Q 2 ), calculated using unquenched jet energies and with no restriction on the angles of particles with respect to the jet axis. Earlier measurements by CMS of jet fragmentation functions [18] in Pb+Pb collisions at the LHC show no significant modification, but the uncertainties on that measurement were not sufficient to exclude modifications at the level of ∼ 10%. CMS recently released a new result [19] using higher statistics data from 2011 that show fragmentation function modifications which are consistent with the results presented in this Letter.

Experimental setup
The measurements presented in this Letter were performed using the ATLAS calorimeter, inner detector, muon spectrometer, trigger, and data acquisition systems [20]. The ATLAS calorimeter system consists of a liquid argon (LAr) electromagnetic (EM) calorimeter covering |η| < 3.2, a steel-scintillator sampling hadronic calorimeter covering |η| < 1.7, a LAr hadronic calorimeter covering 1.5 < |η| < 3.2, and two LAr forward calorimeters (FCal) covering 3.2 < |η| < 4.9. The hadronic calorimeter has three sampling layers longitudinal in shower depth and has a η × φ granularity of 0.1 × 0.1 for |η| < 2.5 and 0.2 × 0.2 for 2.5 < |η| < 4.9. 2 The EM calorimeters are segmented longitudinally in shower depth into three compartments with an additional pre-sampler layer. The EM calorimeter has a granularity that varies with layer and pseudorapidity, but which is generally much finer than that of the hadronic calorimeter. The middle sampling layer, which typically has the largest energy deposit in EM showers, has a granularity of 0.025 × 0.025 over |η| < 2.5.
The inner detector [21] measures charged particles within the pseudorapidity interval |η| < 2.5 using a combination of silicon pixel detectors, silicon microstrip detectors (SCT), and a strawtube transition radiation tracker (TRT), all immersed in a 2 T axial magnetic field. All three detectors are composed of a barrel and two symmetrically placed end-cap sections. The pixel detector is composed of 3 layers of sensors with nominal feature size 50 μm × 400 μm. The SCT barrel section contains 4 layers of modules with 80 μm pitch sensors on both sides, while each end-cap consists of nine layers of double-sided modules with radial strips having a mean pitch of 80 μm. The two sides of each SCT layer in both the barrel and the end-caps have a relative stereo angle of 40 mrad. The TRT contains up to 73 (160) layers of staggered straws interleaved with fibres in the barrel (end-cap). Charged particles with p ch T 0.5 GeV typically traverse three layers of pixel sensors, four layers of double-sided SCT sensors, and, in the case of |η| < 2.0, 36 TRT straws.
Minimum bias Pb+Pb collisions were identified using measurements from the zero degree calorimeters (ZDCs) and the minimum-bias trigger scintillator (MBTS) counters [20]. The ZDCs are located symmetrically at z = ±140 m and cover |η| > 8.3. 1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y axis points upward. Cylindrical coordinates (r, φ) are used in the transverse plane, φ being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2). 2 An exception is the third sampling layer that has a segmentation of 0.2 × 0.1 up to |η| = 1.4.
In Pb+Pb collisions the ZDCs measure primarily "spectator" neutrons, which originate from the incident nuclei and do not interact hadronically. The MBTS detects charged particles over 2.1 < |η| < 3.9 using two counters placed at z = ±3.6 m. MBTS counters are divided into 16 modules with 8 different positions in azimuth and covering 2 different |η| intervals. Each counter provides measurement of both the pulse heights and arrival times of ionization energy deposits.
Events used in this analysis were selected for recording by a combination of Level-1 minimum-bias and High Level Trigger (HLT) jet triggers. The Level-1 trigger required a total transverse energy measured in the calorimeter of greater than 10 GeV. The HLT jet trigger ran the offline Pb+Pb jet reconstruction algorithm, described below, for R = 0.2 jets except for the application of the final hadronic energy scale correction. The HLT trigger selected events containing an R = 0.2 jet with transverse energy E T > 20 GeV.

Event selection and data sets
This analysis uses a total integrated luminosity of 0.14 nb −1 of Pb+Pb collisions recorded by ATLAS in 2011. Events selected by the HLT jet trigger were required to have a reconstructed primary vertex and a time difference between hits in the two sides of the MBTS detector of less than 3 ns. The primary vertices were reconstructed from charged-particle tracks with p ch T > 0.5 GeV. The tracks were reconstructed from hits in the inner detector using the ATLAS track reconstruction algorithm described in Ref. [22] with settings optimized for the high hit density in heavy-ion collisions [23]. A total of 14.2 million events passed the described selections.
The centrality of Pb+Pb collisions was characterized by E FCal T , the total transverse energy measured in the forward calorimeters [23]. Jet fragmentation functions were measured in seven centrality bins defined according to successive percentiles of the E FCal T distribution ordered from the most central to the most peripheral collisions: 0-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, and 60-80%. The percentiles were defined after correcting the E FCal T distribution for a 2% minimum-bias trigger inefficiency that affects the most peripheral events which are not included in this analysis. The performance of the ATLAS detector and offline analysis in measuring jets and charged particles in the environment of Pb+Pb collisions was evaluated using a large Monte Carlo (MC) event sample obtained by overlaying simulated [24] PYTHIA [25] pp hard-scattering events at The detector response to the PYTHIA events was simulated using Geant4 [26], and the simulated hits were combined with the data from the minimum-bias Pb+Pb events to produce 1.2 million overlaid events for each p T interval.

Jet and charged-particle analysis
Charged particles included in the fragmentation measurements were required to have at least two hits in the pixel detector, including a hit in the first pixel layer if the track trajectory makes such a hit expected, and seven hits in the silicon microstrip detector. In addition, the transverse (d 0 ) and longitudinal (z 0 sin θ ) impact parameters of the tracks measured with respect to the primary vertex were required to satisfy |d 0 /σ d 0 | < 3 and |z 0 sin θ/σ z | < 3, where σ d 0 and σ z are uncertainties on d 0 and z 0 sin θ , respectively, obtained from the track-fit covariance matrix. Jets were reconstructed using the techniques described in Ref. [14], which are briefly summarized here. The anti-k t algorithm was first run in four-momentum recombination mode, on η × φ = 0.1 × 0.1 logical towers and for three values of the anti-k t distance parameter, R = 0.2, 0.3, and 0.4.
The tower kinematics were obtained by summing electromagneticscale energies of calorimeter cells within the tower boundaries. Then, an iterative procedure was used to estimate a layer-and η-dependent underlying event (UE) energy density while excluding actual jets from that estimate. The UE energy was subtracted from each calorimeter cell within the towers included in the reconstructed jet. The correction takes into account a cos 2φ modulation of the calorimeter response due to elliptic flow of the medium [23] which is estimated by measurement of the amplitude of that modulation in the calorimeter. The final jet kinematics were calculated via a four-momentum sum of all (assumed massless) cells contained within the jets using subtracted E T values. A correction was applied to the reconstructed jet to account for jets not excluded or only partially excluded from the UE estimate. Then, a final jet ηand E T -dependent hadronic energy scale calibration factor was applied.
After the reconstruction, additional selections were applied for the purposes of this analysis. "UE jets" generated by fluctuations in the underlying event, were removed using techniques described in Ref. [14].
To prevent neighbouring jets from distorting the measurement of the fragmentation functions, jets were required to be isolated. The isolation cut required that there be no other jet within R = 1 having p T > p iso T where p iso T , the isolation threshold, is set to half of the analysis threshold for each R value, p iso T = 42.5, 46, and 50 GeV for R = 0.2, 0.3, and 0.4, respectively. To prevent muons from semileptonic heavy-flavour decays from influencing the measured fragmentation functions, all jets with reconstructed muons having p T > 4 GeV within a cone of size R = 0.4 were excluded from the analysis. To prevent inactive regions in the calorimeters from producing artificial high z fragments, jets were required to have more than 90% of their energy contained within fully functional regions of the calorimeter. Finally, all jets included in the analysis were required to match HLT jets reconstructed with transverse momenta greater than the trigger threshold of 20 GeV. The HLT jets were found to be fully efficient for the jet kinematic selection used in this analysis. Table 1 shows the impact of the cuts on the number of measured jets in central (0-10%) and peripheral (60-80%) collisions. All these cuts together retain more than 96% of all jets.

Jet and track reconstruction performance
The performance of the ATLAS detector and analysis procedures in measuring jets was evaluated from the MC sample using the procedures described in Ref. [14]. Reconstructed MC jets were matched to "truth" jets obtained by separately running the anti-k t Table 2 The relationship between the mean truth-jet transverse momenta, p T jet true , and corresponding reconstructed jet transverse momenta, p T   Table 2. The extracted relationships between p T jet rec and p T jet true will be used in the fragmentation analysis to correct for the average shift in the measured jet energy.
MC studies indicate that the efficiency for PYTHIA jets to be reconstructed and to pass UE jet rejection exceeds 98% for p jet T > 60 GeV in the 0-10% centrality bin. For kinematic selection of jets used in this study, the jet reconstruction was fully efficient. The efficiency for reconstructing charged particles within jets in Pb+Pb collisions was evaluated using the MC sample. Fig. 1 shows comparisons of distributions of four important track-quality variables between data and MC simulation for reconstructed tracks over a narrow p ch T interval, 5 < p ch T < 7 GeV, to minimize the impact of differences in MC and data charged-particle p ch T distributions. The ratios of the data to MC distributions also shown in the figure indicate better than 1% agreement in the η dependence of the average number of pixel and SCT hits associated with the tracks. The distributions of d 0 and z 0 sin θ agree to 10% except in the tails of the distributions, which contribute a negligible fraction of the distribution. For the purpose of evaluating the track reconstruction performance and for the evaluation of response matrices that are used in the unfolding (described below), the reference "truth" particles were taken from the set of final-state PYTHIA charged particles. These were matched to re- constructed charged particles using associations between detector hits and truth tracks recorded by the ATLAS Geant4 simulations. Truth particles for which no matching reconstructed particle was found were considered lost due to inefficiency.
The charged-particle reconstruction efficiency, ε(p T , η), was evaluated separately in each of the seven centrality bins used in this analysis for truth particles within Fig. 2 shows the efficiency as a function of truth-particle p T averaged over |η| < 1 (top) and 1 < |η| < 2.5 (bottom) for the 0-10% and 60-80% centrality bins.
For p T < 8 GeV, ε(p T , η) was directly evaluated using fine bins in p T and η. For p T > 8 GeV the p T dependence of the efficiencies were parameterized separately in the two pseudorapidity intervals shown in Fig. 2 using a functional form that describes trends at low p T as well as at high p T . An example of the resulting parameterizations is shown by the solid curves in Fig. 2. A centralitydependent systematic uncertainty in the parameterized efficiencies, shown by the shaded bands in Fig. 2, was evaluated based on both the uncertainties in the parameterization and on observed variations of the efficiency with p T , which largely result from loss of hits in the SCT at higher detector occupancy. Thus, the systematic uncertainty in the 60-80% centrality bin is small because no significant variation of the efficiency is observed at low detector occupancy, while the uncertainties are largest for the 0-10% centrality bin with the largest detector occupancies.
The efficiencies shown in Fig. 2 decrease by about 12% between the |η| < 1 interval covered by the SCT barrel and the 1 < |η| < 2.5 interval covered primarily by the SCT end-cap. More significant localized drops in efficiency of about 20% are observed over 1 < |η| < 1.2 and 2.3 < |η| < 2.5 corresponding to the transition between the SCT barrel and end-cap and the detector edge respectively. To account for this and other localized variations of the high p T reconstruction efficiency with pseudorapidity, the parameterizations in Fig. 2 for p T > 8 GeV are multiplied by an η-dependent factor evaluated in intervals of 0.1 units to produce ε(p T , η).

Fragmentation functions and unfolding
Jets used for the fragmentation measurements presented here were required to have p jet T > 85, 92 and 100 GeV for R = 0.2, 0.3, and 0.4 jets, respectively. The jet thresholds for R = 0.3 and R = 0.2 jets represent the typical energy measured with the smaller jet radii for an R = 0.4 jet with p T = 100 GeV. Jets were also required to have either 0 < |η| < 1 or 1.2 < |η| < 1.9. The restriction of the measurement to |η| < 1.9 avoids the region at the detector edge with reduced efficiency (|η| > 2.3). The exclusion of the range 1 < |η| < 1.2 removes from the measurement jets whose large-z fragments, which are typically collinear with the jet axis, would be detected in the lower-efficiency η region spanning the gap between SCT barrel and end-cap. While this exclusion does when such a matching jet was found. For each charged particle, the longitudinal jet momentum fraction, z, was calculated according to where R here represents the angle between the charged particle and jet directions. 4 Charged particles from the UE contribute a p ch T -and centralitydependent background to the measurement that must be subtracted to obtain the true fragmentation functions. The contribution of the UE background was separately evaluated for R = 0.2, 0.3, and 0.4 jets in events having at least one such jet above the jet p T thresholds using a grid of R = 0.4 cones that spanned the full coverage of the inner detector. Any such cone having a charged particle with p ch T > 6 GeV was assumed to be associated 4 The R is a boost-invariant replacement for the polar angle θ . with a real jet in the event and was excluded from the UE background determination. The threshold of 6 GeV was chosen to be high enough to avoid bias of the UE p ch T distribution. The resulting per-jet UE charged-particle yields, dn UE ch /dp ch T were evaluated over 2 < p ch T < 6 GeV as a function of p ch T , p jet T , and η jet , averaged over all cones in all events within a given centrality bin according to: dn UE ch dp ch Here N cone represents the number of background cones having a jet of a given radius above the corresponding p jet T threshold, and N cone ch represents the number of charged particles in a given p ch T bin in all such cones evaluated for jets with a given p jet T and η jet .
Not shown in Eq. (5) is a correction factor that was applied to each background cone to correct for the difference in the average UE-particle yield at a given p ch T between the η position of the cone and η jet , and a separate correction factor to account for the difference in the elliptic flow modulation at the φ position of the UE cone and φ jet . That correction was based on a parameterization of the p ch T and centrality dependence of previously measured elliptic flow coefficients, v 2 [23].
By evaluating the UE contribution only from events containing jets included in the analysis, the background automatically has the correct distribution of centralities within a given centrality bin.
The dn UE ch /dp ch T is observed to be independent of p jet T both in the data and MC simulation. That observation excludes the possibility that the upfeeding of jets in p jet T due to the finite JER could induce a dependence of the UE on jet p T . However, such upfeeding was observed to induce in the MC events a p jet T -independent, but centrality-dependent mismatch between the extracted dn UE ch /dp ch T and the actual UE contribution to reconstructed jets. That mismatch was found to result from intrinsic correlations between the charged-particle density in the UE and the MC p jet T error, p jet T = p T jet rec − p T jet true . In particular, jets with positive (negative) p jet T are found to have an UE contribution larger (smaller) than jets with p jet T ∼ 0. Due to the net upfeeding on the falling jet spectrum, the selection of jets above a given p jet T threshold causes the UE contribution to be larger than that estimated from the above-described procedure. The average fractional mismatch in the estimated UE background was found to be independent of p ch T and to vary with centrality by factors between 1.04-1.08, 1.07-1.10, and 1.12-1.15 for R = 0.2, 0.3, and 0.4, respectively. The measured dn UE ch /dp ch T values in the data were corrected by these same factors before being subtracted.
Two different sets of charged-particle fragmentation distributions were measured for each centrality bin and R value: and where N jet represents the total number of jets passing the abovedescribed selection cuts in a given centrality bin, and N ch represents the number of measured charged particles within R = 0.4 of the jets in given bins of p ch T and z, respectively. The efficiency correction, 1/ε, was applied on a per-particle basis using the parameterized MC efficiency, ε(p T , η), assuming p T ch true = p T ch rec . While  that assumption is not strictly valid, the efficiency varies sufficiently slowly with p T ch true that the error introduced by this assumption is 1% everywhere. The D meas (p T ) and D meas (z) distributions were unfolded using a one-dimensional Singular Value Decomposition (SVD) method [27] implemented in RooUnfold [28] to remove the effects of charged particle and jet p T resolution. The SVD method implements a regularized matrix-based unfolding that attempts to "invert" the equation b = Ax, where x, is a true spectrum, b is an observed spectrum, and A is the "response matrix" that describes the transformation of x to b. For D(p T ), the unfolding accounts only for the charged-particle p T resolution and uses a response matrix derived from the MC data set that describes the distri-bution of reconstructed p ch T as a function of MC truth p ch T . The response matrix A(p T ch rec , p T ch true ) is filled using the procedures described in Section 5. The D(z) unfolding simultaneously accounts for both charged particle and jet resolution using a response matrix A(z rec , z true ) with z true (z rec ) calculated using purely truth (fully reconstructed) quantities. A cross-check was performed for the D(z) unfolding that included only the jet energy resolution to ensure that the combination of the two sources of resolution in the one-dimensional unfolding did not distort the result. Because the D meas (z) and D meas (p T ) distributions were already corrected for the charged-particle reconstruction efficiency, the response matrices were only populated with truth particles for which a reconstructed particle was obtained and each entry was corrected for reconstruction efficiency so as to not distort the shape of the true distributions.
To ensure that statistical fluctuations in the MC p T jet true or z true distributions do not distort the unfolding, those distributions were smoothed by fitting them to appropriate functional forms. The truth D(p T ) distributions were fit to polynomials in ln(p T ). The truth D(z) distributions were parameterized using an extension of a standard functional form [29], where a, b, c, d i were free parameters of the fit. The non-standard additional parameter "c" was added to improve the description of the truth distribution at large z. When filling the truth spectra and response matrices, the entries were weighted to match the truth spectra to the fit functions.
The SVD unfolding was performed using a regularization parameter obtained from the ninth singular value (k = 9) of the unfolding matrix. Systematic uncertainties in the unfolding due to regularization were evaluated by varying k over the range 5-12 for which the unfolding was observed to be neither significantly biased by regularization nor unstable. The statistical uncertainties in the unfolded spectra were obtained using the pseudo-experiment method [27]. The largest absolute uncertainty obtained over 5 ≤ k ≤ 12 was taken to be the statistical uncertainty in the unfolded The middle and top right panels in Fig. 3 (Fig. 4) show for R = 0.4 and R = 0.2 jets, respectively the ratios of unfolded  ), compared to the ratios before unfolding.
The unfolding reduces the D(z) ratio slightly at low z but otherwise leaves the shapes unchanged. To evaluate the impact of the unfolding on the difference between central and peripheral fragmentation functions, the middle and bottom right panels in Fig. 3  (Fig. 4) show the ratio of R meas D(z) (R meas D(p T ) ) to R D(z) (R D(p T ) ). Except for the lowest z point, the ratio is consistent with one over the entire z range. Thus, the features observed in R meas D(z) (R meas D(p T ) ), namely the enhancement at low z (p T ) in central collisions relative to peripheral collisions, the suppression at intermediate z (p T ), and the rise above one at large z (p T ) are robust with respect to the effects of the charged particle and jet p T resolution.

Systematic uncertainties
Systematic uncertainties in the unfolded D(z) and D(p T ) distributions can arise due to uncertainties in the jet energy scale and jet energy resolution, from systematic uncertainties in the unfolding procedure including uncertainties in the shape of the truth distributions, uncertainties in the charged particle reconstruction, and from the UE subtraction procedure. The systematic uncertainty due to the jet energy scale (JES) has two contributions, an absolute JES uncertainty and an uncertainty in the variation of the JES from peripheral to more central collisions. The absolute JES uncertainty was determined by shifting the transverse momentum of the reconstructed jets according to the evaluation of the jet energy scale uncertainty in Ref. [30]. The typical size of the JES uncertainty for jets used in this study is 2%. The shift in the JES has negligible impact on the ratios between central and peripheral events of D(p T ) and D(z) distributions whereas it has a clear impact on the D(p T ) and D(z) distributions. At high p T or z the resulting uncertainty reaches 15%. The evaluation of centrality-dependent uncertainty on JES uses the estimates from Ref. [14]. The centrality-dependent JES uncertainty is largest for the most central collisions where it reaches 1.5%. The evaluation of the jet energy resolution (JER) uncertainty follows the procedure applied in proton-proton jet measurements [31]. The typical size of JER uncertainty for jets used in the study is less than 2%. This uncertainty is centrality independent since the dijets in MC are overlayed to real data. The resulting combined systematic uncertainty from JER and centrality-dependent JES on the ratios reaches 6% at high p T and 10% at high z and it has a similar size in the case of D(p T ) or D(z) distributions as in the case of their ratios.
The systematic uncertainty associated with the unfolding is connected with the sensitivity of the unfolding procedure to the choice of regularization parameter and to the parameterization of the truth distribution. The uncertainty due to the choice of regularization parameter was evaluated by varying k over the range 5-12. The typical systematic uncertainty is found to be smaller than 3% or 2% for the D(z) or D(p T ), respectively. The systematic uncertainty due to the parameterization of the truth distribution was determined from the statistical uncertainties of the fits to these distributions. This systematic uncertainty is below 1% or 2% for the D(z) or D(p T ), respectively.
The estimate of systematic uncertainty due to the tracking efficiency follows methods of the inclusive charged particle measurement [23]. The uncertainty is quantified using the error of the fit of tracking efficiency and by varying the tracking selection criteria. In the intermediate-p T region the systematic uncertainty is less than 2%. In the low and high p T region the systematic uncertainty is larger, but less than 8%.
An independent evaluation of potential systematic uncertainties in the central-to-peripheral ratios of D(z) and D(p T ), due to all aspects of the analysis, was obtained by evaluating the deviation from unity of the MC central (0-10%) to peripheral (60-80%) ratios of the fragmentation functions. Since there is no jet quenching employed in MC simulation, the ratios are expected not to show any deviation from unity. No deviation from unity is indeed observed, the largest localized deviation is 4%. To quantify the deviations from unity, the MC R D(z) and R D(p T ) ratios were fit by piecewise continuous functions composed of linear functions defined over the z (p T ) ranges z = 0.02-0.06 (p T = 2-6 GeV), z = 0.06-0.3 (p T = 6-30 GeV), and z > 0.3 (p T > 30 GeV) with parameters constrained such that the linear functions match at the boundaries. The resulting fits are used as estimates of the systematic uncertainties on all measured R D(z) and R D(p T ) ratios reported in Section 8. This systematic uncertainty is certainly correlated with and may overlap with other systematic uncertainties described above.

Discussion
To quantify the effects of the modifications observed in Fig Tables 3 and 4 for R = 0.3 and R = 0.2 jets, respectively. Similar results were obtained for R = 0.4 jets but with larger uncertainties. The results presented in the tables indicate an increase in the number of particles with 0.02 < z < 0.04 of less than one particle per jet in the 0-10% centrality bin relative to the 60-80% centrality bin. A decrease of about 1.5 particles per jet is observed for 0.04 < z < 0.2. The differences between the integrals of the fragmentation functions over 0.4 < z < 1 are not significant relative to the uncertainties. The results for D(z)dz shown in the two tables indicate that in the most central collisions a small fraction, < 2%, of the jet transverse momentum is carried by the excess particles in 0.02 < z < 0.04 for central collisions, but that the depletion in fragment yield in 0.04 < z < 0.2 accounts on average for about 14% of p jet T . To better evaluate the significance of the increase in R D (z) and R D(p T ) above one at large z or p T , average R D(z) and R D(p T ) ratios were calculated by summing the central and peripheral D(z) or D(p T ) distributions over different regions corresponding to the last n points in the measured distributions, n = 2-6. For each resulting average ratio, R D(z) or R D(p T ) , the significance of the deviation from one was evaluated as (  Table 3 Differences of D(z) distributions in different centralities with respect to peripheral events for R = 0.3 jets. The errors represent combined statistical and systematic uncertainties.  statistical and systematic uncertainty. Because there is significant cancellation of systematic uncertainties in the ratios, this analysis provides a more sensitive evaluation of the significance of the large-z excess. For R = 0.4 jets the combined R D(z) (R D(p T ) ), differs from one by approximately 1σ (1.5σ ) for any of the n values. For R = 0.2 jets, R D(z) differs from 1 by approximately 1.5σ for all n values, while R D(p T ) differs from one by 2σ for n = 3-6 corresponding to p T > 47.5 GeV through p T > 20 GeV. The greater significance of the deviations of the R = 0.2 R D(p T ) relative to the R = 0.2 R D(z) and the R = 0.4 R D(z) and R D(p T ) can be attributed to the reduced role of the jet energy resolution in influencing the measurement of the central-to-peripheral ratios for large hadron momenta.
Theoretical predictions for medium modifications of fragmentation functions based on radiative energy loss [32][33][34][35] have generally predicted substantial reduction in the yield of high p T , or large-z fragments and an enhancement at low p T or low z. The predicted reduction at large z generically results from the radiative energy loss of the leading partons in the shower and the resulting redistribution of the jet energy to lower z hadrons. Instead of a reduction, an enhanced yield of high z fragments is seen in the data. However, the difference between observed behaviour at large z and expectations from theoretical calculations may be at least partially attributed to the fact that the fragmentation functions presented in this paper were evaluated with respect to the energies of quenched jets. In contrast, theoretical analyses of the fragmentation functions of quenched jets are typically evaluated in terms of the initial, unquenched jet energies. However, some recent theoretical analyses [36,37] of jet fragmentation functions using quenched jet energies have shown that jet quenching calculations can reproduce the general features observed in the results presented in this Letter. In addition to direct modifications of the fragmentation function due to quenching, the quenching may indirectly alter the fragmentation function of inclusive jets by altering the relative fraction of quarks and gluons.
The simultaneous effects of quenching on the hadron constituents of jets and the measured jet energies may explain a relative increase of experimental fragmentation functions in central collisions at large z as suggested by the data. Jets that fragment to large-z hadrons may lose less energy than typical jets due to reduced formation or colour-neutralization time [38]. Thus, the fragmentation function measured for inclusive jets may have a higher proportion of jets with large-z hadrons. The results in Ref. [36] indicate such an effect that is qualitatively similar to the data.

Conclusions
This Letter has presented measurements by ATLAS of chargedparticle fragmentation functions in jets produced in