Measurement of Charged Hadron Production in Z-Tagged Jets in Proton-Proton Collisions at s=8 TeV

The production of charged hadrons within jets recoiling against a Z boson is measured in proton-proton collision data at ffiffi s p 1⁄4 8 TeV recorded by the LHCb experiment. The charged-hadron structure of the jet is studied longitudinally and transverse to the jet axis for jets with transverse momentum pT > 20 GeV and in the pseudorapidity range 2.5 < η < 4. These are the first measurements of jet hadronization at these forward rapidities and also the first where the jet is produced in association with a Z boson. In contrast to previous hadronization measurements at the Large Hadron Collider, which are dominated by gluon jets, these measurements probe predominantly light-quark jets which are found to be more longitudinally and transversely collimated with respect to the jet axis when compared to the previous gluon dominated measurements. Therefore, these results provide valuable information on differences between quarks and gluons regarding nonperturbative hadronization dynamics.

The production of charged hadrons within jets recoiling against a Z boson is measured in proton-proton collision data at ffiffi ffi s p ¼ 8 TeV recorded by the LHCb experiment. The charged-hadron structure of the jet is studied longitudinally and transverse to the jet axis for jets with transverse momentum p T > 20 GeV and in the pseudorapidity range 2.5 < η < 4. These are the first measurements of jet hadronization at these forward rapidities and also the first where the jet is produced in association with a Z boson. In contrast to previous hadronization measurements at the Large Hadron Collider, which are dominated by gluon jets, these measurements probe predominantly light-quark jets which are found to be more longitudinally and transversely collimated with respect to the jet axis when compared to the previous gluon dominated measurements. Therefore, these results provide valuable information on differences between quarks and gluons regarding nonperturbative hadronization dynamics. DOI: 10.1103/PhysRevLett.123.232001 Quantum chromodynamics (QCD), the theory of the strong interaction, is unique amongst the fundamental forces due to the nonperturbative processes that confine quarks and gluons, collectively referred to as partons, within boundstate hadrons. The parton structure of protons has been the focus of intense research efforts; however, the understanding of how hadrons arise from scattered partons is limited in comparison. Perturbative QCD calculations utilize fragmentation functions to determine cross sections of hadron production from scattered partons. Fragmentation functions describe the probability for a particular parton to transform into a particular hadron [1][2][3]. Several global fits to experimental data have provided parametrized fragmentation functions (see, e.g., Ref. [4] and references therein). However, there is a significant lack of understanding in the mechanisms through which hadrons are formed in the nonperturbative hadronization process and therefore additional data are required.
Fragmentation function studies have been performed using inclusive hadron production at e þ e − colliders, which benefit from a simpler environment free of initial-state gluon radiation [5][6][7][8][9][10][11][12][13]. Semi-inclusive deep-inelastic-scattering measurements have also been used to constrain fragmentation functions at smaller values of Q 2 , the momentum transfer [14,15]. Additionally, inclusive hadron production measurements have been used to study fragmentation functions in the more complex environment, relative to interactions involving leptons, of proton-proton (pp) collisions [16][17][18]. However, such measurements are limited by the lack of an explicit way to relate the scattered parton to the final-state hadron. Measuring fragmentation functions with respect to high transverse momentum (p T ) jets offers a unique opportunity to study hadron production relative to an object that is correlated to the scattered parton. For example, the transverse profile, in addition to the longitudinal dynamics of hadrons within jets, can be used to study fragmentation functions in the longitudinal and transverse directions with respect to the jet axis. Such multidimensional measurements that go beyond inclusive hadrons, or those that consider correlations between particles, have the potential to answer unique questions within QCD related to universality, factorization, and the importance of color-charge flow [19,20].
This Letter reports a study of charged hadrons produced in jets recoiling against a Z boson, also referred to as Z-tagged jets, in the forward region of pp collisions (throughout this Letter the notation Z includes both the Z 0 and virtual γ Ã contributions). The longitudinal momentum fraction, z, the momentum transverse to the jet axis, j T , and the radial distribution, r, of charged hadrons are measured with respect to the jet axis in the laboratory frame, defined as and r ≡ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi Here, p is the three-momentum vector, ϕ is the azimuthal angle, and y is the rapidity. The data sample is selected from an integrated luminosity of approximately 2 fb −1 collected at a center-of-mass energy ffiffi ffi s p ¼ 8 TeV with the LHCb detector in 2012. Events with only one reconstructed primary vertex are analyzed to better identify signatures of a hard twoto-two partonic scattering. Jets are clustered with the anti-k T algorithm [21] using a distance parameter R ¼ 0.5 and are measured differentially in p T for p T > 20 GeV, and in the pseudorapidity range 2.5 < η < 4. (In this Letter, natural units (c ¼ ℏ ¼ 1) are used.) Charged hadrons within the jet are required to have p T > 0.25 GeV, momentum p > 4 GeV, and to lie within the jet cone such that ΔR < 0.5, where tions are unfolded to account for the detector response and to facilitate comparisons with theoretical and numerical predictions. This is the first measurement of charged hadrons within jets produced in association with a Z boson, as well as the first measurement of charged hadrons in jets at these forward pseudorapidities. The Z þ jet process is primarily sensitive to light quark jets, as demonstrated by PYTHIA in this kinematic range [22,23]. Thus, these data provide new and complementary information to previous jet substructure measurements in the inclusive jet channel at midrapidity in hadronic collisions, which are sensitive to primarily gluon jets [24][25][26][27][28][29]. Recent results at midrapidity in the isolated photon-jet channel can also probe fragmentation differences when a photon, rather than a massive vector boson, is measured opposite the jet [30]. The LHCb detector is a single-arm forward spectrometer covering the pseudorapidity range 2 < η < 5, described in detail in Refs. [31,32]. Simulations are used to evaluate the detector performance with regard to the jet reconstruction, track-in-jet reconstruction, and to validate the analysis methods. The simulated pp → Z þ jet þ X events are generated using PYTHIA8 [23] with a specific LHCb configuration [33]. Decays of hadronic particles are described by EVTGEN [34], and final-state radiation in the simulation is generated using PHOTOS [35]. Finally, the GEANT4 toolkit [36] is used to simulate the interactions of the particles with the detector, as described in Ref. [37].
This analysis uses the same data set as that used for the Z þ jet cross section measurement, where events are selected and Z bosons are measured via their dimuon decay as described in Ref. [38]. Candidate events are required to pass a trigger [39], which selects muons with p T > 10 GeV. Only events that contain two high-p T muons are retained. The muons are required to satisfy track-reconstruction and muon-identification criteria, as in Refs. [38,40], and are also required to fall within the fiducial region of 2 < η < 4.5, where the detector performance is well understood. Finally, the dimuon system must have an invariant mass, M μμ , within the range 60 < M μμ < 120 GeV.
Jet reconstruction is performed using a particle flow algorithm [41], where the charged and neutral particles are clustered using the anti-k T algorithm as implemented in Ref. [42]. Reconstructed jets with p T > 15 GeV that lie within 2.5 < η < 4 are analyzed. The 15 < p T < 20 GeV region is included to avoid inefficiencies at the lower p T limit of the measurement in the unfolding procedure. The pseudorapidity requirement ensures that the full jet cone lies within the fiducial area of the LHCb detector. Selection requirements are placed on the jets to reduce the rate of jets not associated with the partonic process producing the Z boson, which is already suppressed by the requirement of a single reconstructed primary vertex in the event. Additionally, the decay muons from the Z boson must not be contained within the jet cone. Only jets that are on the azimuthal away-side of the Z boson, defined by Δϕ Z−jet ≡ jϕ Z − ϕ jet j > 7π=8, are analyzed. The jet energy calibrations are the same as those used in Ref. [38]. Charged hadrons within the jet are identified by the particle flow algorithm utilizing the particle-identification systems and several track-quality criteria [32]. The charged hadrons must also satisfy ΔR < 0.5, which ensures that the corresponding tracks fall within the tracking acceptance.
The methods used to determine the chargedhadron fragmentation distributions are as described in Refs. [24,27,43]. The fragmentation distributions are corrected for tracking inefficiencies and two-dimensionally unfolded for resolution effects, which primarily occur due to the jet energy resolution. The unfolded fragmentation distributions are then normalized by the total number of Z þ jet events in a given jet p T bin, which is determined separately from the hadron-in-jet unfolding procedure. A test of the method described here and below, performed with the reconstructed simulation samples, confirmed that the generated distributions were reproduced for all observables studied, within the statistical uncertainties of the simulated sample. In this analysis, the Z -boson p T is integrated to provide the statistical precision to measure the fragmentation as a function of jet p T . The integral of the fragmentation distributions then corresponds to the mean multiplicity of charged hadrons within the jet.
The number of Z þ jet pairs in each jet p T bin is corrected to account for reconstruction and selection inefficiencies, and is determined independently from and normalizes the fragmentation distributions. The muon detection efficiencies are determined in data using the technique employed in the inclusive weak boson cross section measurements of LHCb [40,44]. The jet reconstruction efficiency is evaluated from simulation, and is greater than 90% for jets with p T > 20 GeV. A correction is also applied to account for differences PHYSICAL REVIEW LETTERS 123, 232001 (2019) 232001-2 between the number of events produced and measured in a given p T bin due to the jet p T resolution. This correction is determined from simulation, and is less than 10%. The method described above is cross checked by comparing the results to a full Bayesian unfolding [45] as implemented in Ref. [46]. The two methods agree to within 1%.
Simulation is used to determine the tracking efficiency and to account for effects from misreconstructed tracks that are incorrectly measured inside or outside of the jet cone. The efficiency is evaluated as a function of momentum and pseudorapidity and applied on a per-track basis. The efficiency decreases for p > 150 GeV due to a requirement on the uncertainty of the track bending radius which is part of the particle flow algorithm and has a larger effect at high momentum. The efficiencies were found to be independent of the observables z, j T and r in the simulated sample. To validate the efficiency corrections procedure, the simulation sample is split in half and the efficiencies are determined with one half and applied to the other. Good recovery of the generated charged hadron distributions in p and η is observed. Within the statistical precision of the sample the tracking efficiency does not depend on the jet p T .
The effects of bin migration in jet p T and in the fragmentation observables on the fragmentation distributions, primarily due to the jet energy and momentum resolutions, are corrected using the two-dimensional Bayesian unfolding method. Response matrices are constructed for each fragmentation observable using simulated samples that study the correlations between the generated and reconstructed yields in bins of ½z; p jet T , ½j T ; p jet T , and ½r; p jet T . Typically the bin migration is less than 5%; however, it can be larger for more extreme values of the fragmentation variables, for example at large z. The number of iterations in the Bayesian unfolding procedure is selected to be the minimum number for which the relative change in the fragmentation functions at z ≈ 0.05 is smaller than 0.2% per additional iteration in all of the jet p T bins. Based on this criterion, the unfolding is iterated seven times for each observable.
Systematic uncertainties that arise from the uncertainties on the various efficiencies are assigned to the number of Z þ jet pairs measured in each jet p T bin. The uncertainty in the muon reconstruction efficiency is negligible. Systematic uncertainties on the jet reconstruction are evaluated as in Ref. [38] by comparing the jet reconstruction quality requirements in simulation and data. Similarly to the muon efficiencies, the precision with which the uncertainty of the jet reconstruction corrections are determined, due to the limited simulation sample size, is also evaluated; however, this is found to be negligible compared to the jet reconstruction quality requirement uncertainty of 1.9%. The normalization is not corrected for Z þ jet background events, and thus a systematic uncertainty of 1.7% is assigned for the impurity of both Z bosons and jets in the measurement, as determined in Refs. [38,40]. Effects from pile up are also studied and found to be negligible. The total normalization uncertainty of 2.7% is determined by adding these components in quadrature.
The jet-energy scale and resolution are also considered as sources of systematic uncertainty. The jet-energy scale and its uncertainty have been studied in previous measurements of the Z þ jet cross section [38,41]. To estimate these effects, the scale is varied by one standard deviation of its uncertainty. New unfolding matrices are constructed with this modification, and the difference in the fragmentation distributions determined with the modified and nominal response matrices is taken as a systematic uncertainty. Similarly, the systematic uncertainty due to the jet-energy resolution is evaluated by smearing each component of the jet momentum by an additional term corresponding to the uncertainty on the jet resolution and constructing new response matrices. The difference between the nominal and smeared unfolded charged hadron-in-jet distributions is taken as the systematic uncertainty on the jet-energy resolution.
The unfolding method is validated with two different tests. The first test is performed by splitting the simulated sample in two and using one half to generate the response matrices with which the other half is unfolded. Recovery of the generator-level fragmentation distributions is observed and average deviations from perfect agreement are 2%, which is assigned as an uncertainty related to the unfolding procedure. A second test is performed by splitting the simulated sample in half by Z -boson p T , and performing a similar test to the previous one to check for any uncertainty associated with the assumed prior. The results again deviate from perfect agreement by about 2%, confirming that a 2% systematic uncertainty on the unfolding procedure is appropriate.
The track selection requirements, tracking efficiency, and charged-hadron identification are also studied as sources of systematic uncertainty. The track selection uncertainty is assigned by requiring a tight fake-track removal criterion and repeating the analysis. The differences in the final fragmentation distributions with and without this requirement are taken as systematic uncertainties. The track selection uncertainty is typically less than 5%; however, it reaches a maximum of approximately 8% at some values of z, j T , and r in the highest jet p T bin studied. The systematic uncertainty on the tracking efficiency is determined by smoothing the two-dimensional efficiency and repeating the analysis, which accounts for the statistical precision with which the efficiency is determined. The resulting distributions are compared to the nominal distributions and the differences are taken as uncertainties on the tracking efficiency; these are generally less than 3% but rise up to 10% in some bins. Uncertainties associated to misidentifying charged hadrons are also considered by comparing the nominal fragmentation PHYSICAL REVIEW LETTERS 123, 232001 (2019) functions to those obtained when hadron-to-lepton (and vice versa) misidentification probabilities are considered. These uncertainties are less than 5%, except at large z where the charged pion-to-electron misidentification probability becomes larger [32]. Figure 1 shows the distributions of z in three jet p T bins. These illustrate that the longitudinal momentum fraction is approximately constant as a function of jet p T at high z. At low z the fragmentation functions differ, which is a kinematic effect due to the requirement that the track momentum be greater than 4 GeV; therefore, higher p T jets can probe smaller z. This also reflects that the charged hadron multiplicity increases with jet p T . Comparing these forward measurements to inclusive jet measurements at central rapidity from ATLAS [24] indicates that the fragmentation functions are not as steeply falling at high z [47]. This may reflect differences between light-quark and gluon fragmentation. Figures 2 and 3 show the j T and r distributions of charged hadrons within jets. The j T profiles show a rounded peak at small j T which transitions to a perturbative tail at larger j T as also seen in Ref. [48]. This is indicative of an observable that can be treated in the so-called transversemomentum-dependent framework [1][2][3]49], where sensitivity to both a large and small transverse momentum scale is necessary. The radial profiles show that the number of charged hadrons at small r is highly dependent on jet p T ; however, the values are relatively constant as a function of jet p T at nearly all other values of r. Interestingly, the j T fragmentation distributions are similar to the central pseudorapidity inclusive jet results; however, these measurements are more collimated in r than the inclusive jet measurements [47]. This behavior in r is correlated to the flatter fragmentation in z and may be a reflection of the different pseudorapidity region or differences in light-quark and gluon fragmentation. We note that the comparisons to the measurements by ATLAS should be qualitative in nature, rather than quantitative, due to the slightly different kinematic criteria placed on each of the measurements. The distributions in j T and r offer the opportunity to study the interplay between perturbative parton shower and nonperturbative hadronization dynamics. For example, the steeply falling tail of the j T distributions results from a combination of perturbative radiation and nonperturbative hadronization processes.
The fragmentation functions are compared to predictions from PYTHIA8 Z þ jet events, where the details of the PYTHIA8 configuration can be found in Ref. [47]. These comparisons are made since the specific LHCb tune contains realistic experimental conditions [33] and also shows that the unfolding procedure is not simply correcting the measured distributions to the predictions from PYTHIA8. An example of the comparison as a function of z is shown in Fig. 4; all of the comparisons described in this text can be found in Ref. [47]. In general, PYTHIA8 underestimates the number of charged hadrons at high z; PYTHIA8 also underestimates the number of charged hadrons at small r.  Comparisons of the data to predictions from PYTHIA8 as a function of j T show a consistent shape, but in general PYTHIA underestimates the number of charged hadrons in each bin by approximately 20%. The integral of the ratio of the PYTHIA8 predictions to the data is always less than unity, which is a reflection of the underestimation of the mean charged hadron multiplicity in PYTHIA8. In summary, the production of charged hadrons in jets recoiling against a Z boson is measured in ffiffi ffi s p ¼ 8 TeV pp collisions by the LHCb experiment. The jets are measured in the fiducial region of 20 < p T < 100 GeV and 2.5 < η < 4, while the hadrons are required to have p T > 0.25 GeV, p > 4 GeV, and to be located within the jet cone of distance parameter R ¼ 0.5. The longitudinal momentum fraction, momentum transverse to the jet axis, and radial profile of the charged hadrons are measured with respect to the jet axis. These results provide insight into hadronization mechanisms as they probe a new kinematic regime. They additionally probe a high fraction of light-quark jets versus gluon jets when compared to midrapidity inclusive jet measurements in the same jet p T range. The results are compared to predictions from the PYTHIA8 event generator with a specific LHCb configuration, and show that the simulation underestimates the number of high momentum hadrons. Additionally, comparisons with inclusive midrapidity gluon-dominated jet measurements indicate that light quark-dominated jets recoiling against a Z boson at forward rapidity are more collimated in both z and r [47]. This work lays the foundation for a broader hadronization research program at LHCb, utilizing the excellent tracking, particle identification, and heavy-flavor jet tagging capabilities already demonstrated by the LHCb detector [32,50].
We express our gratitude to our colleagues in the CERN accelerator departments for the excellent performance of the LHC. We thank the technical and administrative staff at the LHCb institutes. We acknowledge support from CERN and from the national agencies: CAPES, CNPq, FAPERJ, and FINEP (Brazil)