Properties of jet fragmentation using charged particles measured with the ATLAS detector in pp collisions at s =13 TeV

This paper presents a measurement of quantities related to the formation of jets from high-energy quarks and gluons (fragmentation). Jets with transverse momentum 100 GeV < p T < 2 . 5 TeV and pseudorapidity j η j < 2 . 1 from an integrated luminosity of 33 fb − 1 of ﬃﬃﬃ s p ¼ 13 TeV proton-proton collisions are reconstructed with the ATLAS detector at the Large Hadron Collider. Charged-particle tracks with p T > 500 MeV and j η j < 2 . 5 are used to probe the detailed structure of the jet. The fragmentation properties of the more forward and the more central of the two leading jets from each event are studied. The data are unfolded to correct for detector resolution and acceptance effects. Comparisons with parton shower Monte Carlo generators indicate that existing models provide a reasonable description of the data across a wide range of phase space, but there are also significant differences. Furthermore, the data are interpreted in the context of quark- and gluon-initiated jets by exploiting the rapidity dependence of the jet flavor fraction. A first measurement of the charged-particle multiplicity using model-independent jet labels (topic modeling) provides a promising alternative to traditional quark and gluon extractions using input from simulation. The simulations provide a reasonable description of the quark-like data across the jet p T range presented in -this measurement, but the gluon-like data have systematically fewer charged particles than the simulation.


I. INTRODUCTION
Jets are collimated sprays of particles resulting from high-energy quark and gluon production. The details of the process that underlies the fragmentation of quarks and gluons with net quantum chromodynamic (QCD) charge into net neutral hadrons are not fully understood. Jet formation is a complex multiscale problem, including important contributions from QCD effects that cannot be described by perturbation theory. Measuring basic quantities related to fragmentation is therefore essential to furthering our understanding of the emergent properties of QCD.
Perturbative and nonperturbative physically inspired models have free parameters that are tuned to data in order to best describe the radiation pattern inside jets [1]. This is in turn an important input to all analyses at the Large Hadron Collider (LHC) due to the ubiquity of jets. Measurements of jet substructure in proton-proton (pp) collisions at a center-of-mass energy of ffiffi ffi s p ¼ 7 TeV [2][3][4][5] have already been used by the ATLAS collaboration for parameter optimizations (tunes) of the PYTHIA 8 Monte Carlo (MC) generator [6]. A measurement of the average number of charged particles inside jets at ffiffi ffi s p ¼ 8 TeV [7] was also used as input to recent developments in the HERWIG 7 MC program [8]. Further measurements of jet constituent multiplicity and energy sharing will provide powerful constraints for future generator optimizations. Quark-and gluon-initiated jets (henceforth quark and gluon jets) have different radiation patterns (see e.g., Ref. [9]). As many analyses at the LHC target either quark-enriched or gluon-enriched processes, these radiation-pattern differences can be useful for jet tagging [10,11]. Measurements of jet structure can be used to calibrate quark-versus-gluon jet taggers. By exploiting the rapidity dependence of the relative quark and gluon jet rates, ATLAS [7] extracted the average chargedparticle multiplicity for quark and gluon jets separately. This was then combined with detector-level systematic uncertainties to provide quark/gluon tagger uncertainties at ffiffi ffi s p ¼ 13 TeV [12]. A more complex tagger based on several jet shapes could be calibrated in a similar manner using extended results. The benefit of a particle-level measurement is that a portion of the calibration can be independent of ATLAS and LHC operating conditions. Uncertainties in detector effects can be updated with the changing detector environment. Adding more observables and measuring their differential distributions will improve this calibration.
Although the full radiation pattern inside jets is not calculable from first principles, the energy dependence of many observables can be calculated in perturbation theory. There have been significant theoretical advances in softcollinear effective theory (SCET) [13][14][15][16] to derive factorization theorems that describe the evolution of universal nonperturbative functions [17][18][19][20]. This was applied to the measurement of jet charge at ffiffi ffi s p ¼ 8 TeV [21]. There have also been predictions and comparisons with the jet transverse momentum (p T ) dependence of the average number of charged particles inside jets (see Ref. [7] and references therein). This quantity does not have a perturbative expansion in the usual sense (as a series in α S ); instead there is a series expansion in ffiffiffiffiffi α S p [22,23]. This behavior is predicted for a wide class of Sudakov safe observables [24]. At least for the case of charged-particle multiplicity, this nonstandard expansion seems to be an excellent model of the data [7]. The goal of this paper is to measure properties of jet fragmentation using charged-particle tracks. Such properties have been measured at many colliders at various centerof-mass energies, including the SPS [25][26][27], PETRA [28,29], PEP [30][31][32][33], TRISTAN [34], CESR [35], LEP [36][37][38][39][40][41][42][43][44][45][46][47], HERA [48,49], and the Tevatron [50][51][52][53] 21,61] in pp collisions and are also compared with jet fragmentation measured in Pb þ Pb collisions [54,56-58, 62,63] and p þ Pb collisions [56]. The measurement presented here represents a significant extension of previous work. In particular, the accessible jet energy range is increased due to the larger ffiffi ffi s p ¼ 13 TeV. There are enough events in the 2016 dataset to probe the substructure of jets with p T up to 2.5 TeV. Next, the precision of the measurement has improved due to advances in track reconstruction inside jets during the long shutdown between LHC runs 1 and 2, including the additional insertable B-layer (IBL) detector [64,65] and new algorithms for tracking inside dense environments [66][67][68]. Furthermore, detailed experimental studies to derive uncertainties in all aspects of tracking inside jets extend the capabilities of previous measurements to a wider region of phase space and also allow differential analyses [67,69]. These new data therefore probe broader and deeper aspects of the radiation pattern inside jets across an extended phase space. The paper is organized as follows. Section II introduces the observables to be measured. Then, following a brief description of the ATLAS detector in Sec. III, the data and simulation samples are documented in Sec. IV. Chargedparticle track, jet, and event reconstruction are detailed in Sec. V. Corrections for detector effects (unfolding) are documented in Sec. VI. A description of the corresponding systematic uncertainties can be found in Sec. VII and the results are presented in Sec. VIII. Section IX provides conclusions and future outlook.

II. OBSERVABLES
This analysis builds upon the previous ATLAS jet structure measurements presented in Refs. [7,21,59]. The fundamental quantity is the fragmentation function D h p ðz; EÞ, which describes the probability of finding a hadron h with energy fraction z of the parton p that has energy E. At a hadron collider, the jet transverse momentum, p T , 1 is a better proxy for the starting scale (μ) of jet evolution. To avoid confusion with previous measurements of similar observables, the transverse momentum fraction is denoted in this paper by the symbol ζ ¼ p particle T =p jet T . The fragmentation function itself, like parton distribution functions (PDFs), cannot be calculated from first principles in perturbation theory. However, it has a Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution [70][71][72] and so the p T dependence of many observables can be calculated. In particular, where P p 0 ←p ðζ; μÞ are the Altarelli-Parisi splitting functions [70] and depend on the scale μ through α S . Charged particles are studied because they provide a way to measure single hadrons inside the jet (as opposed to calorimeter energy deposits, which can result from multiple particles) which gives access to P h D h p . A basic quantity related to the fragmentation function is the charged-particle multiplicity. The average chargedparticle multiplicity is an integral over ζ and a sum over h and p of D h p ðζÞ. An extension of the multiplicity is the set of ζ moments of D. The zeroth moment is the average multiplicity. The full distribution of multiplicity depends on (multihadron) fragmentation functions in a complicated way; a more direct probe of D is to measure hadron production as a function of ζ, which is a sum of D over p and h (but no integral over ζ). Additional observables are also studied in order to probe the angular spread of jet fragmentation beyond the collinear limit. All of the observables are described below. 1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the center of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the center of the LHC ring, and the y-axis points upwards. Cylindrical coordinates ðr; ϕÞ are used in the transverse plane, ϕ being the azimuthal angle around the z-axis. The pseudorapidity is defined in terms of the polar angle θ as η ¼ − ln tanðθ=2Þ. Angular distance is measured in units of ΔR ≡ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi ðΔηÞ 2 þ ðΔϕÞ 2 p .
where f p is the fraction of parton type p at a given jet p T . The multiplicity is not calculable in perturbation theory, but to lowest order in ffiffiffiffiffi α S p , the ratio of the multiplicity for quark-initiated jets to that for gluon-initiated jets is the ratio of color factors C A =C F ¼ 9=4. The fraction of quark jets increases with p T , which decreases the inclusive multiplicity. However, this is compensated by an inherent increase in the multiplicity with p T for both quark and gluon jets [73]. In addition to the mean, the full ð1=N jet ÞdN jet =dn ch distribution is measured.

B. Summed fragmentation function
The distribution of the momentum fraction ζ is studied inside jets summed over charged-hadron types. The quantity that is measured is Fðζ; p jet T Þ ¼ ð1=N jet Þdn ch =dζ. In terms of the fragmentation function, By definition, R dζFðζÞ ¼ hn ch i. In addition to measuring the distribution FðζÞ in bins of jet p T , summary statistics of the FðζÞ distribution are extracted to show how the distribution evolves with jet p T . The following properties of the ζ distribution are extracted: (i) Partial fractions of FðζÞ: R X 0 FðζÞdζ= R FðζÞdζ ¼ n ch ðζ < XÞ=n ch to show how much of the jet energy is carried by particles of a given p T fraction. For illustration, the values considered are X ∈ f0.1; 0.01; 0.001g. As X → 1, these partial fractions become a constant value of 1.0, independent of the jet p T . (ii) Moments of FðζÞ: The distribution of FðζÞ is nearly normally distributed in logðζÞ, which means that it is defined by its first two moments [73]. For this reason, κ ¼ 2 is measured as a function of the jet p T . For illustration, the case κ ¼ 1=2 is also considered. (iii) Weighted sums over the jet: The values considered are κ ∈ f1=2; 2g. The observable P i∈jet ζ 2 i is often called p D T and can be used for quark/gluon jet tagging [74]. For a given jet type, these observables increase monotonically with increasing jet p T for κ ≲ 1 and decrease monotonically for κ ≳ 1 (see Sec. VIII B); the κ values chosen are representative of these trends.
Each of these derived quantities is extracted from the measured FðζÞ distribution. More details about the procedure for unfolding these derived quantities are presented in Sec. VI.
C. Transverse momentum p rel T ≡ p charged particle T sin Δϕ, where Δϕ is the angle between the momentum of the constituent charged particle and the jet axis in the transverse plane. The quantity that is measured is fðp rel T ; p jet T Þ ¼ ð1=N jet Þdn ch =dp rel T . The average value is defined by hp rel

D. Radial profile
The number of charged particles in various annuli around the jet axis. The quantity that is measured is ρ ch ðr; p jet T Þ ¼ ð1=N jet Þdn ch =2πrdr, where r ¼ ΔRðcharged particle; jetÞ. The average value is defined by hri ¼ R rρ ch ðrÞ= R ρ ch ðrÞ. The last two quantities are not simple derivatives of the fragmentation function as they additionally depend on finite opening angles encoded in the dθ=θ emission phase space. Since quantities are measured as a function of jet p T which is defined using charged and neutral particles, the observables are sensitive to the charged-to-neutral fraction inside jets. However, this fraction is robust to mismodeling as isospin is an approximate symmetry of the strong force.

III. ATLAS DETECTOR
The ATLAS detector [75] at the LHC covers nearly the entire solid angle around the collision point. It consists of an inner tracking detector surrounded by a thin superconducting solenoid, electromagnetic and hadronic calorimeters, and a muon spectrometer incorporating three large superconducting toroidal magnets. The inner-detector system (ID) is immersed in a 2 T axial magnetic field and provides charged-particle tracking in the range jηj < 2.5.
The high-granularity silicon pixel detector covers the vertex region and typically provides four measurements per track, the first hit being normally in the IBL. It is followed by the silicon microstrip tracker which usually provides eight measurements per track. These silicon detectors are complemented by the transition radiation tracker, which enables radially extended track reconstruction up to jηj ¼ 2.0.
The calorimeter system covers the pseudorapidity range jηj < 4.9. Within the region jηj < 3.2, electromagnetic calorimetry is provided by barrel and end cap highgranularity lead/liquid-argon (LAr) calorimeters, with an additional thin LAr presampler covering jηj < 1.8, to correct for energy loss in material upstream of the calorimeters. Hadronic calorimetry is provided by the steel/ scintillating-tile calorimeter, segmented into three barrel structures within jηj < 1.7, and two copper/LAr hadronic end cap calorimeters. The solid angle coverage is completed with forward copper/LAr and tungsten/LAr calorimeter modules optimized for electromagnetic and hadronic measurements respectively.
Interesting events are selected to be recorded by the firstlevel trigger system implemented in custom hardware, followed by selections made by algorithms implemented in software in the high-level trigger [76]. The first-level trigger reduces the 40 MHz bunch crossing rate to below 100 kHz, which the high-level trigger further reduces in order to record events to disk at about 1 kHz.

IV. DATASETS AND SIMULATED SAMPLES
These measurements use the dataset of pp collisions recorded by the ATLAS detector in 2016, corresponding to an integrated luminosity of 33 fb −1 at a center-of-mass energy of ffiffi ffi s p ¼ 13 TeV. Events are only considered if they are collected during stable beam conditions and satisfy all data quality requirements. Due to the high instantaneous luminosity and the large total inelastic proton-proton cross section, on average there are about 25 simultaneous (pileup) collisions in each bunch crossing. The measurements presented in this paper use a variety of MC samples for estimating correction factors as well as for comparison with the corrected data. Dijet events were generated at leading order (LO) with PYTHIA 8.186 [77], with the 2 → 2 matrix element convolved with the NNPDF2.3LO PDF set [78] and using the A14 tune of multiple-parton-interaction and shower parameters [6]. PYTHIA uses a p T -ordered parton shower model. Additional dijet events were simulated using different generators, in order to study the impact of modeling uncertainties. SHERPA 2.1 [79] events were generated using multileg 2 → 2 and 2 → 3 matrix elements, which were matched to parton showers following the Catani-Krauss-Kuhn-Webber prescription [80]. These SHERPA events were simulated using the CT10 PDF set [81] and the default SHERPA event tune. HERWIG++ 2.7 [82,83] was used to provide a sample of events with an angle-ordered parton shower model. These events were generated with the 2 → 2 matrix element, convolved with the CTEQ6L1 PDF set [84] and configured with the UE-EE-5 tune [85].
All simulated events were passed through a full simulation of the ATLAS detector [86] implemented in GEANT 4 [87], which describes the interactions of particles with the detector and the subsequent digitization of analog signals. The effects of multiple simultaneous pp collisions were simulated with inelastic pp collisions using the PYTHIA 8.186 generator with the A2 [88] set of tuned parameters and the MSTW2008LO [89] PDF set; these events were overlaid on the nominal dijet events.

V. OBJECT AND EVENT SELECTION
Since the data are unfolded to particle level, it is necessary to define both the particle-level and detector-level objects used in the measurement. The former are chosen to be as close as possible to the latter in order to minimize the model dependence caused by an extrapolation from the measured phase space at detector level to the phase space at particle level. Section VA describes the definition of chargedparticle tracks and jets. Following the discussion of objects, Sec. V B describes the particle-level and detector-level event selection criteria.

A. Object reconstruction
While it is not possible to separate the underlying event from the hard scatter at particle level, it is possible to remove the contribution from pileup. Therefore, the unfolding target is particle-level distributions produced in single proton-proton interactions. However, at detector level, there is ambiguity about which pp collision vertex corresponds to the hard-scatter event. Collision vertices are reconstructed from tracks in the inner detector. Each vertex is required to be associated with at least two tracks with p T > 0.4 GeV. The primary hard-scattering vertex of the event is chosen to be the one with the highest P p T 2 calculated using all tracks associated with the vertex. Particle-level jets are built from MC-simulated stable particles (cτ > 10 mm) excluding muons and neutrinos. By definition, particles from pileup and from interactions with the detectors are not included. These jets are clustered using the anti-k t [90] algorithm with radius parameter R ¼ 0.4 as implemented in FASTJET [91]. Detector-level jets are built from topological calorimeter-cell energy clusters [92] using the same algorithm as is used at particle level. A series of simulation-and data-based correction and calibration factors are applied to ensure that the resulting jet p T is the same as the particle-level value on average [93]. Jets are required to have p T > 60 GeV so that the rate of jets originating from pileup is negligible. The detector-level phase space includes one bin at low jet p T (60-100 GeV) which is not in the fiducial phase space of the measurement due to the large impact of migrations into and out of the acceptance.
Charged particles are used to compute the particle-level definitions of all observables if they are clustered within a particle-level jet and have p T > 500 MeV and jηj < 2.5. The detector-level analog to charged particles is tracks. Tracks are reconstructed from hits in the inner detector (see e.g., Ref. [67]) and a series of quality criteria are applied to the selected tracks to reject those originating from hits due to multiple charged particles (fake tracks) and from pileup. The transverse momentum resolution is approximately σðp T Þ=p T ≈ 0.05% × p T =GeV ⊕ 1%, with a significant degradation in the core of high-p T jets due to challenges associated with pattern recognition. 2 Tracks are required to pass the tight primary selection as well as the 2 For example, the p T resolution is approximately 30% at 100 GeV when five or more particles are within ΔR < 0.015. loose track-to-vertex association [94]. In particular, they must have p T > 500 MeV and jηj < 2.5, and the number of pixel and strip clusters associated with the track is required to be at least 9 (11) for jηj < 1.65 (≥1.65). In addition, the transverse impact parameter d 0 relative to the beam line must be less than 2 mm and the longitudinal impact parameter, z 0 , is required to satisfy jz 0 sin θj < 3 mm. Tracks are matched to jets via ghost association [95]. This matching procedure creates ghost versions of the tracks with the same direction but infinitesimal p T . Jet clustering is repeated and tracks are assigned to the jet that contains their ghosted version. For the isolated, high p T jets used in this measurement, ghost association is nearly identical to a geometric matching based on ΔR < 0.4.

B. Event selection
Particle-level events are required to have at least two jets with jηj < 2.1 (within the tracking detector acceptance) and the leading two such jets must satisfy p lead T =p sublead T < 1.5. This jet-p T balance requirement simplifies the interpretation of the final state in terms of a 2 → 2 scattering process.
Detector-level events are selected using single-jet triggers. Due to the large cross section for jet production, most of the jet triggers are prescaled: events that pass the trigger are randomly discarded with a fixed probability. The trigger used for a particular jet p T is chosen to ensure that the trigger is 100% efficient (for the measurement phase space and prior to prescaling) and has the lowest prescale factor. Events in data are weighted by the prescale. The lowest-threshold unprescaled jet trigger is used for jets with p T > 600 GeV. Detector-level events are required to pass the same selection requirements as particle-level events: there must be at least two off-line calibrated jets with jηj < 2.1 and the leading two of these jets must satisfy p lead T =p sublead T < 1.5. Figure 1 shows the basic kinematic properties of the two leading jets passing this event selection compared with various MC predictions at detector level.
The substructure of the two leading jets is used in the analysis. Figure 2 shows detector-level distributions for a selection of the observables that were introduced in Sec. II. For the jets with p T ∼ 1 TeV shown in Fig. 2, the most probable number of tracks is about 15 and the most probable momentum fraction is about 1%. The radiation pattern is peaked at the center of the jet, so both the p rel T and r distributions are peaked at zero. The PYTHIA, HERWIG++, and SHERPA distributions generally bracket the data and are accurate to within about 20%.
In order to expose differences between quark and gluon jets, the more forward and more central of the two jets are distinguished and measured separately. Figure 3 shows the gluon-jet fraction as a function of jet p T and jet η (more details about quark/gluon definitions are given in Sec. VIII B). For a fixed jet p T , higher-jηj jets are more often quark initiated due to valence quarks scattering off gluons. For a fixed η, the quark fraction increases with jet p T due to the relative increase in valence-quark scattering off a quark or gluon compared with gluon-gluon scattering. Table I summarizes the object and event selections from Secs. VA and V B.

VI. UNFOLDING
The data are corrected for resolution and acceptance effects, and the fiducial phase space of the measurement is described by the particle-level object and event selection in Sec. V. Equation (2) symbolically summarizes the unfolding procedure for a binned distribution x:  The transverse momentum, p T , spectrum for the selected jets; the simulation is normalized to the data. The normalization is dominated by the first bin, which accounts for the overall offset in the other bins for some of the predictions. Right: The pseudorapidity, η, distribution for the selected jets, split into the more forward and the more central of the two jets. Error bars only include the statistical uncertainty.  The gluon-jet fraction as a function of jet transverse momentum p T and jet pseudorapidity η. Right: the fraction of the more forward and the more central jets that are gluon initiated. The error bars in the right plot represent the uncertainty computed from the 100 NNPDF2.3LO replicas [78]. See Sec. VIII B for more details about quark/gluon definitions and uncertainties.
where n jets;unfolded is the unfolded number of forward or central jets (depending on the bin), determined by the number of entries in the n ch unfolding (as there is one entry per jet). The symbols θ and ϵ represent the unfolding matrix and correction factors, described in more detail below. The jet substructure observables are simultaneously unfolded with the jet p T and for the more forward and the more central jets at the same time. For an observable with n bins bins in a given p T bin, this results in a total of N total ¼ 2 × ðn bins Þ × ðp T binsÞ bins. All of these bins are concatenated to form a one-dimensional input. To begin the unfolding, the data are corrected for the fraction of events that pass the detector-level selection but not the particlelevel selection, ϵ reco not true . This also corrects for nondijet events, but their rate is negligible. Then, an iterative Bayesian (IB) unfolding technique [96] is used as a regularized matrix inversion to correct for the detector resolution in events that pass both the detector-level and particle-level selections. The IB method is implemented in the ROOUNFOLD framework [97] with the unfolding matrix θ and one iteration, is chosen to minimize the total uncertainty. After the application of the response matrix, a final correction is applied to account for the fraction of events that pass the particle-level but not detector-level selection, ϵ true not reco . The resulting unfolded measurement is reorganized into individual distributions with n bins per p T bin for each of the more forward and more central jets. The jet p T is also unfolded in parallel and each p T bin of the jet substructure observable is normalized by the number of measured jets in that bin. For n ch , this renders the distributions normalized to unity per jet p T bin; for the other observables, the normalization in each p T bin is (up to acceptance effects) hn ch i, as discussed in Sec. II.
To illustrate the jet-p T dependence of the measured observables, the evolution with the jet p T of various moments (κ) is computed using Eq. (3): where the sum is over all i that correspond to p T bin j. Since the bin center is used to calculate the average, a correction c binning is applied to account for the difference between the bin center and the mean of the distribution within the bin. This correction is calculated using PYTHIA, and is computed by reweighting PYTHIA so that it agrees with the unfolded distribution. For ζ, p rel T , and r, Eq. (3) represents the κ moment for individual particles. For ζ, the jet-based moments are also computed: h P i∈jet ζ κ i i. For these jet-based moments, Eq. (3) is modified by removing the denominator P n bins i¼1 x unfolded;i . By construction, the κ ¼ 0 jet-based moment of ζ is the κ ¼ 1 moment of n ch . The binning correction factor is mostly near unity, deviating by less than 1% for n ch and up to about 10% for the other observables. Figure 4 shows the response matrix normalized per particle-level bin. As stated above, the observable bins are concatenated with the jet p T and for both the more forward and more central jets to form a one-dimensional distribution that is unfolded. A diagonal stripe represents events where the detector-level jet p T is the same as the particle-level value; off-diagonal components represent jet p T migrations. Within a jet p T bin, there is a small dependence on ζ, with a worse resolution at high ζ due to the deteriorating momentum resolution at high track p T . The diagonal strips in the upper left and lower right quadrants correspond to events where the more forward particle-level jet is the more central detector-level jet and vice versa. This migration happens in about 1% of events. Within a given jet p T bin, the migrations to neighboring ζ bins are small. Except at high ζ and high jet p T where the migrations can reach 50%, the off-diagonal components of the response matrix are about 10%.

VII. UNCERTAINTIES
Systematic and statistical uncertainties are assessed for each step of the analysis, including the acceptance Ghost-associated tracks, p T > 500 MeV and jηj < 2.5 Event selection At least two jets, with the leading two satisfying p lead T =p sublead T < 1.5 Jet selection Leading two, separated by η (more forward/central) correction factors, response matrix, and unfolding method. For each uncertainty, some component of the analysis chain is varied and then the entire unfolding procedure is repeated. Data and simulation statistical uncertainties are determined from pseudo-experiments using the bootstrap method [98]. The details of the experimental systematic uncertainties related to track and jet reconstruction are given in Sec. VII A and the uncertainties in the unfolding method and fragmentation modeling are described in Sec. VII B. An additional source of uncertainty arising from binning effects is evaluated when computing the average value of an observable as a function of jet p T . The average values are determined using the bin centers, so the correction described in Sec. VI relies on the simulation for the distribution within a given bin. An uncertainty in the binning correction is estimated by comparing the correction factors derived from PYTHIA with those from SHERPA, where both simulations are reweighted to match the unfolded data distribution. Figure 5 provides an overview of the systematic uncertainties for a selection of observables, using the average value vs p T for illustration. The uncertainty in the rate of fake and secondary tracks is the leading experimental reconstruction uncertainty for hn ch i and hri except at low jet p T where the uncertainties from the inclusive tracking efficiency and the unfolding procedure are larger. The jet energy uncertainties are the most important for ζ, with the tracking uncertainties matching in size in the highest jet-p T bins. The tracking and jet energy uncertainties are about the same size for hp rel T i. Fragmentation modeling uncertainties are large for hn ch i at low jet p T and for hζi at high jet p T . While the size of the binning correction uncertainty is less than 2% for hp rel T i and hri, it is still the dominant uncertainty for these observables. Further details about each uncertainty source are given below and the full covariance matrices, including all correlation information, are made available in Ref. [99].

A. Track and jet reconstruction
Except for ζ, the jet energy is only used to determine the p T bin. Since the fragmentation properties vary slowly with jet p T , the resulting impact of jet energy scale and resolution uncertainties on the analysis is often less important than other sources of uncertainty. Nonetheless, the impact of a 19-parameter decomposition of the jet energy scale uncertainty was evaluated [93]. Six of these 19 components are due to in situ constraints on the jet energy scale from various multiobject balance studies, such as Z þ jets. Additional sources of uncertainty are related to pileup, jet flavor, and extrapolations to high p T . The total uncertainty in the jet energy scale is about 1% for jets with p T between 100 and 1000 GeV and the impact on this measurement is much less than 1% except at high ζ, where it can reach as high as 2%. The impact of the jet energy resolution is determined from an ensemble of event samples with jet energies smeared within the uncertainty.
The most important experimental uncertainties are related to track reconstruction and cover the track reconstruction efficiency, the rate of fake and secondary tracks, the momentum scale, and density effects from pixel and strip cluster merging. In the PYTHIA simulation, approximately 60% of the charged particles/tracks inside jets are charged pions that are well matched, 3 10% are wellmatched kaons, 5% are well-matched protons, 15% are charged particles that are not matched to reconstructed tracks (inefficiency), 5% are secondaries (split equally between photon conversions and nuclear interactions), 1% are not well-matched tracks (fake tracks), and about Oð0.1%Þ are pileup tracks wrongly matched to the primary hard-scattering vertex. The pileup contribution decreases with jet p T and momentum fraction, but increases with jet cone size (reaching 1% at Δ ¼ 0.4). In contrast, the faketrack rate increases slightly with jet p T and has a contribution at high-momentum fraction of a few percent from The probability distribution for the detector-level combined transverse momentum fraction ζ and jet transverse momentum p T distribution normalized in bins of the particlelevel variable using the PYTHIA simulation. The ζ distribution is concatenated with the jet p T so that every 21 bins is a different jet-p T bin. The first 315 bins represent the more central of the two jets and the second 315 bins correspond to the more forward jet. The z-axis is truncated at 10 −3 for visualization only, to aid readability. 3 Reconstructed tracks are matched to charged particles by examining the pattern of sensors where energy was deposited. If over 50% of the weighted number of measurements on a track are due to one charged particle, it is declared matched to the track. The weights are chosen to reflect the amount of information present in each detector and are ten for the pixel detector, five for the strip detector, and one for the straw tube tracker. kinked tracks reconstructed with a very high p T . The reconstruction inefficiency grows with jet p T , and is peaked at both low and high radial distance from the center of the jet and is reduced at high momentum fraction. This is because tracks with a larger radial distance from the jet axis tend to have lower p T (larger material effects and thus lower efficiency), while tracks in the core of the jet suffer from an inefficiency in the pattern recognition in the dense environment.
The uncertainty in the inclusive track reconstruction efficiency is dominated by the uncertainty in the amount of material in the inner detector. Variations in the amount of material that are consistent with detector construction knowledge and measurements from secondary vertices [100] result in an uncertainty of 0.5% for jηj < 0.1, which grows to 2.7% for 2.3 < jηj < 2.5. This uncertainty is applied in the simulation by randomly removing tracks with a p T -and jηj-dependent probability. This uncertainty dominates the n ch measurement for jet p T ≲ 1 TeV.
Since the ATLAS pixel detector measures the charge collected from ionization, it is possible to constrain the inefficiency from density effects by looking for single tracks with pixel charge consistent with two minimumionizing particles [67]. The resulting uncertainty is about 0.4% for tracks with ΔR < 0.1 and is validated with additional studies related to the charged-to-neutral ratio in the jet as well as the geometric orientation of pixel clusters [68]. This uncertainty is most important for the radial energy measurement at small radii from the jet axis and for the n ch measurement in the highest jet-p T bins.
The rate of fake tracks is studied inside jets by inverting some of the track quality criteria such as the fit χ 2 =NDF and  is found to agree between data and simulation at the 30% level [68]. A related source of uncertainty is due to the rate of secondary tracks. These tracks originate from real charged particles, but are the result of interactions in detector material and not direct fragmentation processes. The rate of secondaries is estimated by fitting the track d 0 distribution and is found to agree with simulation within about 30%. These rates are then varied to determine an uncertainty in the measurement. The fake-track rate is the leading source of uncertainty for n ch when p T ≳ 1 TeV and when ζ ∼ 1 or r ≲ 0.05 for all jet-p T bins. Uncertainties related to the modeling of pileup have a negligible impact.
The leading source of uncertainty in the track parameters is in the q=p T (q is the electric charge) from a potential sagitta distortion due to detector misalignment weak modes [94]. This bias is corrected and the uncertainty in the correction is about 0.1=TeV except at ϕ ≈ 0 and jηj ∼ 2.5 where the correction can reach 1=TeV.
The impact on the measurement is smaller than the other tracking uncertainties.

B. Unfolding method and fragmentation modeling
An uncertainty resulting from the unfolding method described in Sec. VI is determined by unfolding the prediction from a reweighted simulation with the nominal procedure. The reweighted simulation is constructed by modifying the nominal PYTHIA 8 particle-level spectrum so that the simulated detector-level spectrum, from propagating the reweighted particle-level spectrum through the response matrix, has significantly improved agreement with the data. The modified detector-level distribution is unfolded with the nominal response matrix and the difference between this and the reweighted particlelevel spectrum is an indication of the bias due to the unfolding method (in particular, the choice of prior) [101]. The weights are chosen by comparing the PYTHIA 8 particle-level spectrum with the unfolded data. After applying the reweighting, the χ 2 =NDF calculated using only the statistical uncertainties improves significantly in each jet p T bin. The resulting systematic uncertainties are generally much smaller than the detector-level differences between the data and simulation, as desired.
The unfolded result depends on the modeling of jet fragmentation through the prior, the response matrix, and the correction factors. Variations in the prior are already accounted for in the data-driven nonclosure uncertainty described above. The other contributions are evaluated by comparing the result using PYTHIA 8 with the result using the alternative HERWIG++ sample described in Sec. IV. A similar uncertainty is obtained when using HERWIG++ or SHERPA as the alternative model. This comparison is decomposed into components corresponding to varying only the response matrix or only the initial/final correction factors, ϵ reco not true and ϵ true not reco in Eq. (2). All three components are added in quadrature to determine the total uncertainty due to fragmentation modeling. Even though these sources of uncertainty are correlated, they were treated as independent because the level of correlation is unknown given that there are only two alternative models. The resulting uncertainty is much smaller than the difference between PYTHIA 8 and HERWIG++ at particle level. For n ch , the response matrix is the dominant contribution to this uncertainty, except in the first jet-p T bin where the correction factors and their uncertainty are also important. For the per-particle observables (ζ, r, p rel T ), the correction factors dominate the uncertainty because acceptance effects are much more important.

VIII. RESULTS
The unfolded data are presented in two ways. Section VIII A focuses on the inclusive spectra for both   jets together, while Sec. VIII B uses the differences between forward and central jets to determine the unique features of quark-initiated and gluon-initiated jets, some of which can be compared with perturbative QCD calculations. These sections show a selection of jet p T bins; a complete set of results can be found in Ref. [99].

A. Inclusive distributions
The unfolded averages of the measured observables are presented as a function of the jet p T in Fig. 6 for the more forward and more central jets separately and then combined in Fig. 7. All other figures in this section combine measurements of both jets. The more central jets show properties that are more gluon-like than the more forward jets: they have a larger charged-particle multiplicity and a softer momentum-fraction spectrum. The modeling of the all-jet spectra is very similar to that of the more forward/backward jets and is described in detail for the alljet spectra only.
As the jet p T increases, the average charged-particle multiplicity increases, the average momentum fraction decreases, the average p rel T increases, and the average multiplicity-weighted radius decreases. Charged-particle multiplicity increases from about 10 at jet p T of 100 GeV to just over 20 at 2.5 TeV. In most cases, PYTHIA 8 and SHERPA bracket the data, and are accurate to better than 10%; HERWIG++ is often between these two and closer to the data. As the distribution of n ch is almost Poissonian, nearly all of the information about the distribution is encoded in the mean value. In contrast, the distribution of ζ is more complicated. 4 The average   The distribution is nearly Gaussian in log ζ, so it is well specified by two parameters instead of one [73]. momentum fraction is about 5% at jet p T of 100 GeV and decreases to about 2.5% at 2.5 TeV (the most probable value, shown below, is lower). The distributions of p rel T and the radial profiles fall steeply (nearly exponentially) away from zero and the average values in Fig. 7 give a sense of how fast they fall (exponential distributions are uniquely specified by their mean). The average p rel T at p jet T ¼ 100 GeV is about 0.35 GeV and increases to about 0.55 GeV at p jet T ¼ 2.5 TeV. If the angular distribution about the jet axis is independent of p T , the average value of p rel T should be proportional to hζiðp jet T Þ × p jet T . This would suggest an increase by a factor of ð2.5%=5%Þ × ð2500=100Þ ∼ 12.5 across the measured range; instead it only increases by a factor of about 1.5. This means that the angular distribution is not independent of p T and in particular, the jets become more collimated. This is also consistent with direct measurement of the radial profile, where the average value drops from about 0.06 at p jet T ¼ 100 GeV to about 0.03 at p jet T ¼ 2.5 TeV. While PYTHIA 8, SHERPA, and HERWIG++ agree well with the data for p rel T , SHERPA provides a poorer model of the average radial profile as a function of the jet p T .
As noted above, the distribution of ζ cannot be described simply by its average value, in contrast to n ch , p rel T and r, which are nearly Poisson or exponentially distributed. Therefore, it is useful to summarize the p jet T dependence of other aspects of the ζ distribution. Figure 8 shows partial integrals of the ζ distribution and Fig. 9 shows the average values of ζ 1=2 , ζ 2 , P i∈jet ζ 1=2 and P i∈jet ζ 2 . Figure 8 illustrates how the average fraction of charged particles  P i∈jet ζ 2 as a function of the jet p T . Section II presents the definitions of both classes of observables. The lower panels show the ratio of various MC predictions to the data, with the total uncertainty band centered on the data at unity. The values are computed from the moments of the unfolded distributions in each jet-p T bin, with a small correction added to account for binning effects. with a given momentum fraction evolves with p jet T . There is no correction for binning effects, as the measured ζ distribution has bin edges which nearly align with 0.1%, 1%, and 10%. In particular, the ζ bins are 1=1.5 n , for n ¼ 0; …; 21, and the fractions in Fig. 8 are estimated as 0.1% ≈ 1=1.5 17 , 1% ≈ 1=1.5 11 , and 10% ≈ 1=1.5 5 . The fraction of particles carrying 10% or less of the momentum changes very little across the entire p jet T range and is also near unity (>90% for all p jet T ). A strong p jet T dependence is introduced when the ζ threshold is lowered to 1% and to 0.1%. Since charged particles are required to have p T > 500 MeV, only jets with p jet T > 500 GeV can have particles with ζ < 0.1%. The fraction of particles with ζ < 1% has a logarithmic increase while the fraction of particles with ζ < 0.1% appears to increase faster than linearly with p jet T . Both of these general trends are reproduced by PYTHIA 8, SHERPA, and HERWIG++, although for example, PYTHIA 8 disagrees with the exact value at low p jet T for ζ < 1% and all p jet T for ζ < 0.1%. For the ζ < 0.1% case, SHERPA and PYTHIA 8 bracket the data, with SHERPA predicting more particles with a lower ζ fraction, while HERWIG++ is much closer to the data. The average values of ffiffi ffi ζ p and ζ 2 for individual particles as a function of jet p jet T in the top panel of Fig. 9 show a decreasing trend that is qualitatively similar to the trend for the average ζ in Fig. 7. For ffiffi ffi ζ p , PYTHIA 8/HERWIG++ and SHERPA bracket the data, although PYTHIA 8 agrees with the data within the uncertainty. SHERPA predicts a significantly higher average ζ 2 than is present in the data. As with hn ch i, the average value of P i∈jet ζ 1=2 increases with the jet p T , while the p T dependence of h P i∈jet ζ 2 i is more complicated as it first decreases and then slowly increases with jet p T . This trend is well reproduced by PYTHIA and HERWIG++, but not by SHERPA.      To present more differential information, the full unfolded distributions for n ch , ζ, p rel T , and r are shown in Figs. 10-13, respectively, for representative p jet T bins. Many of the relevant trends are captured in the above discussion about the p jet T dependence of the moments. However, finer information that may be useful for generator tuning is provided by the differential distributions.

B. Quark and gluon distributions
As discussed in Sec. V B, the more forward and the more central of the two selected jets can be separated to study differences between the radiation patterns within quark and gluon jets. Using the fraction of quark jets f q in the two jet samples (forward f and central c), one can extract the quark (h q i ) and gluon (h g i ) jet fragmentation properties separately by solving a system of equations per bin i of an observable: where f x q is the fraction of quark jets in sample x (see Fig. 3 for the gluon fraction) and the nominal fractions are taken from the default PYTHIA simulation described in Sec. IV. The flavor of a jet is defined as the type of the highestenergy parton from the event record (all partons prior to hadronization) matched to the jet via ghost association. This definition is not unique because quark and gluon labels are not universal due to color connections with other partons in the event. 5       However, for isolated jets, the topology dependence is predicted to be much smaller than the difference between quark and gluon jets [102]. and h c i from the unfolding method, uncertainties in the extracted h g i and h q i distributions arise from the PDF choice, from the matrix elements, from the fragmentation model (due to flavor changing), and from the method nonclosure. The determination of the uncertainty from the choice of PDF uses the NNPDF uncertainty set [NNPDF 2.3 at LO in QCD and QED with α S ðm Z Þ ¼ 0.119] and the matrix-element uncertainty is estimated by comparing the nominal fractions from PYTHIA with those from HERWIG. 6 The nonclosure uncertainty is due to the small (subpercent level) differences between forward and central quark jets, as well as forward and central gluon jets, resulting from an η dependence in the jet fragmentation at a fixed jet p T [102]. When presenting the average properties in bins of jet p T , the binning correction described in Sec. VI is also applied and the corresponding uncertainty contributes to the total uncertainty (though it is smaller than other sources of uncertainty).

ATLAS
The matrix-element uncertainty dominates the total uncertainty in the extraction procedure, resulting in an uncertainty that is about 1% at high jet p T and about 5% at low to moderate jet p T for quark jets, with the inverse trend for gluon jets (low uncertainty at low jet p T and large uncertainty at high jet p T ). The extractions presented here use leading-order matrix elements and leading-logarithm parton shower programs; higher-order effects that modify the fractions f are not included in this leading-order extraction. Figure 14 shows the extracted quark and gluon distributions for jets with 1000 GeV < p jet T < 1200 GeV.   To reinforce the simulation dependence of these extractions, the data distributions are referred to as "extracted quark-like" and "extracted gluon-like."

ATLAS
A key challenge with the extraction based on Eqs. (4) and (5) is that it strongly depends on simulation for the fractions f q and f g . A new approach that does not require the input of any fractions is topic modeling [103,104], which holds great promise for learning about quark-and gluon-like jets with less input from theory. In this approach, one can extract distributions of "topics" T 1 and T 2 : in the limit that min j fh g j =h q j g ¼ min j fh q j =h g j g ¼ 0, h T 1 ¼ h q and h T 2 ¼ h g . When this is not exactly the case, the topics are universal but not pure combinations of quarks and gluons. The extracted topics using n ch in two jet p T bins are shown in Fig. 15. The very low n ch region is dominated by quarks and the very high n ch region is dominated by gluons and therefore n ch nearly has the property that min j fh g j =h q j g ≈ min j fh q j =h g j g ≈ 0. Therefore, the first topic is well aligned with quarks and the second topic is more gluon-like. This alignment is better for quarks than for gluons, but the second topic does converge to the gluon distribution at high jet p T . Other observables aside from n ch are not considered for topic modeling because there are no bins where h g j =h q j ¼ 0 or h q j =h g j ¼ 0 is approximately true and therefore the topics do not align with quark-and gluon-like quantities.       While the full quark and gluon distributions presented in Fig. 14 cannot be predicted from perturbative QCD, it is possible to model the p jet T dependence of the moments of the ζ distribution. Positive moments of the fragmentation function have a perturbative evolution with a proper α S power series via DGLAP-like equations. In general, there are two terms that contribute to the right-hand side of Eq. (1) that prevent an analytic solution: one term proportional to D h g and one term proportional to D h q , where the coefficients for the κ sums are the Mellin transformsP p 0 ←p ðκÞ ¼ R 1 0 dζζ κ P p 0 ←p ðζÞ for p 0 ¼ g and p 0 ¼ q, respectively. For gluon jets (p ¼ g), the g → qq 0 splitting function is finite, 7 so jP g←g j ≫ jP q←g j for κ ≠ 0.8 whereP g←g switches sign. Therefore, away from κ ≈ 0.8 and in the modified leading-logarithm approximation 8 [73,108,109],

ATLAS
where β 0 is the first term in the QCD β-function and Λ is a nonperturbative parameter (of order Λ QCD ). The predictions FIG. 14. The extracted quark-and gluon-like distributions of (top left) charged-particle multiplicity n ch , (top right) transverse momentum fraction ζ, (bottom right) transverse momentum p rel T , and (bottom right) the radial profile in bins of the radial distance r from the jet axis for jets with transverse momentum 1000 GeV < p jet T < 1200 GeV. The quark-and gluon-jet distributions from PYTHIA are also shown for comparison. The uncertainty bands on the data include both statistical and systematic uncertainties while the error bars are due to statistical uncertainties only. 7 The splitting functionP q←q is also finite, but is not numerically small compared withP g←q except when κ is very small so this case is not considered further. 8 This means resummation that includes the leading-order splitting functions and the first-order running of the strong coupling. A more refined calculation [18,19] using SCET [13][14][15][16] and fragmenting jet functions [105][106][107] is possible. However, the deviations from this simple approach are higherorder corrections and do not qualitatively change the comparisons in this section.   FIG. 15. The jet topics extracted using charged-particle multiplicity n ch for jets in two transverse momentum, p T , bins together with the topics and quark and gluon distributions from PYTHIA. The uncertainty bands on the data include both statistical and systematic uncertainties while the error bars are due to statistical uncertainties only. FIG. 16. The dependence on transverse momentum p T of the average extracted gluon-like transverse-momentum-fraction weighted sum P i∈jet ζ κ i for (top left) κ ¼ 0.5, (top right) κ ¼ 1.0, and κ ¼ 2.0 (bottom). For comparison, the results from PYTHIA and a simple leading-logarithm (LL) calculation are also presented. The prediction is normalized to the data in the sixth jet p T bin, called the anchor bin and indicated by an arrow. The uncertainty band on the calculation is from varying Λ in Eq. (6) up and down by a factor of 2 from its nominal value of 400 MeV (in most regions, this band is not much wider than the linewidth and thus not visible). The uncertainty bands on the data include both statistical and systematic uncertainties while the error bars are due to statistical uncertainties only. are scaled to match the data in the sixth jet p T bin (referred to as the "anchor bin"). There is no a priori reason to select any particular bin as the anchor bin so one of the first bins after the lowest-threshold unprescaled jet trigger is selected. Figure 16 shows the distributions of the average P i∈jet ζ κ i for κ ¼ 0.5, 1.0, and 2.0 for gluon jets. As mentioned above,P g←g is predicted to change sign at κ ¼ 0.8, a trend which is supported by the data: for low κ, the average value increases with p T and when κ is large, the average decreases with p T . For κ ¼ 1, momentum conservation and isospin symmetry predict that the average value of P i∈jet ζ i should be constant and approximately 2=3, the ratio of charged pions to all pions. 9 The leadinglogarithm (LL) calculation predictsP g←g ð1Þ ≈ 0 so the p T dependence is already negligible compared with the κ ¼ 0.5 and κ ¼ 2 cases.
When κ → 0, both the quark and gluon fragmentationfunction Mellin transforms diverge and so the analysis with Eq. (6) is not accurate. The κ → 0 limit is hn ch i and there is no known series in α S to describe its p jet T dependence. Despite this, the anomalous dimension for the p jet T dependence of hn ch i has been calculated to "N 3 LO" where the series is in ffiffiffiffiffi α S p instead of α S [22,23]. Figure 17 shows hn ch i as a function of p jet T for both extracted quark-like and gluon-like jets as well as the topic extraction along with the prediction for the pure quark/gluon case. Gluon jets from data deviate significantly from simulation and from the calculation at high jet p T ; this is also true to a lesser extent for quark jets, which seem to have a different slope that is most prominent at low jet p T . A similar trend was first observed in Ref.
[7], albeit with lower precision in the highest p T bins. There are several possibilities for this discrepancy, such as an unaccounted for potential source of bias in the quark/gluon jet fractions. The data in the right panel of Fig. 17 do not yet conclusively support or reject this hypothesis; with more data, it may be possible to determine if the data match topic 2 in PYTHIA or deviate as is the case for gluons in the left panel of Fig. 17.
The p T dependence of the average ζ, p rel T , and r are shown in Fig. 18. Gluon jets have more constituents than quark jets on average so their average ζ is lower. For both quark and gluon jets, hζi decreases with the jet p T in part because constituent multiplicity increases with p T . Gluon jets are wider than quark jets on average, but both quark and gluon jets become denser with increasing jet p T . The data show nearly the same trends as PYTHIA in all cases. The dependence on jet transverse momentum p T of the mean charged-particle multiplicity hn ch i for quark and gluon jets (topic 1 and topic 2) in data and in PYTHIA as well as from a calculation using perturbative QCD. The calculation cannot predict the overall normalization and therefore the prediction is normalized to the data in the sixth p T bin, called the anchor bin and indicated by an arrow. The binning correction is not applied to the average topics, as this correction is very sensitive to fluctuations due to the limited number of simulated events. The uncertainty bands on the data include both statistical and systematic uncertainties while the error bars are due to statistical uncertainties only.

IX. CONCLUSION
This paper documents a measurement of track-based jet fragmentation functions in pp collisions at ffiffi ffi s p ¼ 13 TeV.
The analysis uses a dataset corresponding to an integrated luminosity of 33 fb −1 recorded by the ATLAS detector at the LHC. Multiple jet properties, including the chargedparticle multiplicity, the momentum fraction carried by charged particles, and angular properties of the radiation pattern inside jets are studied. There are key areas where there are significant disagreements between the ATLAS default MC simulation (PYTHIA 8.2 with the A14 tune, HERWIG++, and SHERPA) and the data, especially for the radial profiles and momentum distributions in SHERPA. The radial profile is systematically broader in data than in simulation, but the momentum transverse to the jet axis and the momentum fraction are well modeled within the precision of this measurement. Near 1 TeV in jet p T , these measurements have achieved percent-level uncertainties for a variety of observables. In addition to measuring the forward, central, and combined jet distributions, the forward and central jet spectra are considered separately to study quark-and gluon-like distributions. A first measurement of topic modeling for the charged-particle multiplicity provides a promising alternative to traditional methods of extracting quark-and gluon-jet distributions that use input from simulation. The simulations provide a reasonable description of the quark-like data across the jet p T range presented in this measurement, but the gluon-like data have systematically fewer charged particles than the simulations by about 10%. The unfolded data are made public through HEPDATA to provide input to help improve both perturbative and nonperturbative aspects of fragmentation modeling in the future.

ACKNOWLEDGMENTS
We thank CERN for the very successful operation of the LHC, as well as the support staff from our institutions without whom ATLAS could not be operated efficiently. We acknowledge the support of ANPCyT, Argentina; YerPhI, Armenia; ARC, Australia; BMWFW and FWF, Austria; [7] ATLAS Collaboration, Measurement of the chargedparticle multiplicity inside jets from ffiffi ffi s p ¼ 8 TeV pp tation in ϒ → ggγ and e þ e − → qqγ events at ffiffi ffi s p ¼ 10 GeV, Phys. Rev. D 56, 17 (1997 B 388, 659 (1996).
[42] DELPHI Collaboration, Energy dependence of the differences between the quark and gluon jet fragmentation, Z. Phys. C 70, 179 (1996 [59] ATLAS Collaboration, Measurement of the jet fragmentation function and transverse profile in proton-proton collisions at a center-of-mass energy of 7 TeV with the ATLAS detector, Eur. Phys. J. C 71, 1795 (2011).