Measurement of inclusive jet production and nuclear modifications in pPb collisions at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_{_\mathrm {NN}}} =5.02\,\mathrm{TeV} $$\end{document}sNN=5.02TeV

Inclusive jet production in pPb collisions at a nucleon–nucleon (NN) center-of-mass energy of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_{_\mathrm {NN}}} =5.02\,\mathrm{TeV} $$\end{document}sNN=5.02TeV is studied with the CMS detector at the LHC. A data sample corresponding to an integrated luminosity of 30.1 nb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{-1}$$\end{document}-1 is analyzed. The jet transverse momentum spectra are studied in seven pseudorapidity intervals covering the range \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ -2.0<\eta _\mathrm {CM}< 1.5$$\end{document}-2.0<ηCM<1.5 in the NN center-of-mass frame. The jet production yields at forward and backward pseudorapidity are compared and no significant asymmetry about \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _\mathrm {CM} = 0$$\end{document}ηCM=0 is observed in the measured kinematic range. The measurements in the pPb system are compared to reference jet spectra obtained by extrapolation from previous measurements in pp collisions at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}=7\,\mathrm{TeV} $$\end{document}s=7TeV. In all pseudorapidity ranges, nuclear modifications in inclusive jet production are found to be small, as predicted by next-to-leading order perturbative QCD calculations that incorporate nuclear effects in the parton distribution functions.


Introduction
Jet measurements play an important role in the study of the quark gluon plasma (QGP) produced in relativistic heavy ion collisions. A key observable in these studies is the phenomenon of jet quenching [1][2][3][4][5][6], in which the partons produced in hard scattering lose energy through gluon radiation and elastic scattering in the hot and dense partonic medium. Jet quenching was first observed at RHIC through measurements of high transverse momentum ( p T ) hadrons [7] and dihadron correlations [8]. At the LHC, this phenomenon was observed more directly as dijet momentum imbalance [9,10] and photon-jet energy imbalance [11] in PbPb collisions. An important ingredient in understanding how the presence of a hot QCD medium affects the jets is the comparison to reference measurements from collision systems that are not expected to produce the QGP. Most often, pp collisions at the same center-of-mass energy are used as a reference. Modifie-mail: cms-publication-committee-chair@cern.ch cations in jet yields [12,13], shapes [14], and fragmentation patterns [15,16] in PbPb collisions have been found in comparison to expectations based on pp measurements. These modifications are found to depend on the overlap between the colliding nuclei, and are largest in the most central (i.e., largest overlap) PbPb collisions.
The interpretation of the jet modification results in nucleus-nucleus collisions and the understanding of their relation to the properties of the QGP requires detailed knowledge of all nuclear effects that could influence the comparisons with the pp system. Nuclear modifications may already be present at the initial state of the collisions, independently of QGP formation. Such modifications are collectively referred to as cold nuclear matter (CNM) effects and include parton energy loss and multiple scattering before the hard scattering, and modifications of the parton distribution functions in the nucleus (nPDFs) with respect to those of a free nucleon (PDFs). Some nPDF modifications have been previously deduced from measurements of lepton-nucleus deep inelastic scattering and Drell-Yan production of lepton pairs from qq annihilation in proton-nucleus collisions [17]. In addition, measurements of π 0 production in deuteron-gold collisions at RHIC [18] are also included in recent nPDF fits to better constrain the nuclear gluon distributions [19]. There are several ranges in the parton fractional momenta x in which the data show suppression or enhancement in the nPDFs relative to the proton PDFs. At small x ( 0.01), the nPDFs are found to be suppressed, a phenomenon commonly referred to as "shadowing" [20]. In the range 0.02 x 0.2, the nPDFs are enhanced ("antishadowing" [17]), and for x 0.2 a suppression has been seen ("EMC effect" [21]).
Proton-lead (pPb) collisions at the LHC provide an opportunity to evaluate the CNM effects and establish an additional reference for the interpretation of measurements performed in PbPb collisions. The results of several pPb studies involving jets or dijets [22][23][24], electroweak bosons [25,26], and high p T charged particles [27,28] are already available. No significant indication of jet quenching was found so far in the pPb studies of inclusive jet production [22,29], dijet momen-tum balance [23], dijet acoplanarity [23,24], or chargedhadron measurements [27,28]. The shapes of the dijet [23] and Z boson [25] pseudorapidity distributions are found to be in better agreement with EPS09 nPDF predictions [19] than with the free-proton PDFs for measurements inclusive in the impact parameter. Hints of modifications larger than those presently included in the EPS09 nPDFs have also been seen [25][26][27]. In particular, the charged hadron spectra [27] are found to be enhanced at high p T beyond the anti-shadowing included in EPS09. Significant modifications with respect to those included in EPS09 have also been found for impact-parameter-dependent measurements [22,23]. The interpretation of the latter results is more difficult because of the kinematic biases introduced through the event selections [22,23,[30][31][32].
In this paper we present the CMS measurements of inclusive jet production in pPb collisions at a nucleon-nucleon (NN) center-of-mass energy of √ s NN = 5.02 TeV as a function of p T in several pseudorapidity regions in the range −2.0 < η CM < 1.5 in the NN center-of-mass system. No additional event activity selections have been made to avoid the associated kinematic biases. The measurements extend in p T up to 500 GeV/c and are sensitive to nPDF modifications in the anti-shadowing and EMC effect regions. Since presently there are no experimental results available from pp collisions at √ s = 5.02 TeV, pp reference jet spectra in pseudorapidity ranges corresponding to the present measurements are obtained by extrapolating jet measurements at √ s = 7 TeV [33]. The paper is organized as follows: Sect. 2 provides the experimental details, Sect. 3 gives an account of the systematic uncertainties in the measurements, Sect. 4 presents the results, and Sect. 5 summarizes our findings.

Data analysis
This measurement is based on a data sample of pPb collisions corresponding to an integrated luminosity of 30.1 nb −1 collected by the CMS experiment in 2013. The beam energies were 4 TeV for protons and 1.58 TeV per nucleon for lead nuclei, resulting in a center-of-mass energy per nucleon pair of 5.02 TeV. The direction of the higher-energy proton beam was initially set up to be clockwise within CMS conventions, and was reversed after a data set corresponding to an integrated luminosity of 21 nb −1 was recorded. As a result of the energy difference of the colliding beams, the nucleonnucleon center-of-mass in the pPb collisions is shifted with respect to zero rapidity in the laboratory frame. Both portions of the data set are analyzed independently and the results are found to be compatible within their uncertainties. In order to reduce the statistical uncertainties, the two data sets are then combined. Results from the first data taking period are reflected along the z-axis so that in the combined analysis the proton travels in the positive z and pseudorapidity η direction. In the laboratory frame η = − ln[tan(θ/2)], where θ is the polar angle defined with respect to the proton beam direction. The results are presented in this convention, after transformation to the NN center-of-mass frame, which for massless particles is equivalent to a shift in pseudorapidity: η CM = η − 0.465.

Experimental setup
A detailed description of the CMS detector and of its coordinate system can be found in Ref. [34]. It features nearly hermetic calorimetric coverage and high-resolution tracking for the reconstruction of energetic jets and charged particles. The calorimeters consist of a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter (HCAL) with coverage up to |η| = 3. The quartz/steel hadron forward (HF) calorimeters extend the calorimetry coverage in the region 3.0 < |η| < 5.2, and are used in offline event selection. The calorimeter cells are grouped in projective towers of granularity η × φ = 0.087 × 0.087 (where φ is the azimuthal angle in radians) for the central pseudorapidity region used in the present jet measurement, and have coarser segmentation (about twice as large) at forward pseudorapidity. The central calorimeters are enclosed in a superconducting solenoid with 3.8 T magnetic field. Charged particles are reconstructed by the tracking system, located inside the calorimeters and the superconducting coil. It consists of silicon pixel and strip layers covering the range |η| < 2.5, and provides track reconstruction with momentum resolution of about 1.5 % for highp T particles. Muons are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid.

Event selection
The CMS online event selection employs a hardware-based level-1 (L1) trigger and a software-based high-level trigger (HLT). A minimum bias sample is selected by the L1 requirement of a pPb bunch crossing at the interaction point and an HLT requirement of at least one reconstructed track with p T > 0.4 GeV/c in the pixel tracker. This minimum bias trigger was prescaled by a large factor for most of the 5.02 TeV data collection, because of the high instantaneous luminosity of the LHC. In order to increase the p T range of the measurement, additional HLT triggers were used to select events based on the presence of a jet with p T > 20, 40, 60, 80, or 100 GeV/c reconstructed in the calorimeters.
For the offline analysis, an additional selection of hadronic collisions is applied by requiring a coincidence of at least one HF calorimeter tower with more than 3 GeV of total energy on the positive and negative sides of the interaction point. Events are further required to have at least one reconstructed primary vertex with at least two associated tracks [35]. A maximum distance of 15 cm between the primary vertex and the nominal interaction point along the beam line is required to ensure maximum tracking acceptance. Additionally, trackbased selection cuts are applied to suppress of beam-related background events [36]. The instantaneous luminosity of the pPb run in 2013 resulted in a 3 % probability of at least one additional interaction occurring in the same bunch crossing. Events with more than one interaction ("pileup" events) are removed using a rejection algorithm developed in Ref. [27]. The pileup-rejection efficiency of this filter is found to be 90 ±2 % in minimum bias events and it removes a very small fraction (0.01 %) of the events without pileup. In order to combine the spectra measured from the various jet-triggered data samples, the events included in the analysis are weighted according to the individual HLT prescale factors corresponding to the trigger object with maximum p T in the event. The top panel of Fig. 1 shows the prescale-weighted jet spectra that are reconstructed with the anti-k T [37] algorithm from each HLT trigger path and the combined inclusive jet spectrum. The ratios of each HLT-triggered spectrum to the combined jet spectrum are shown in the bottom panel of Fig. 1. In the range of p T where the triggers are fully efficient, this ratio is unity and independent of jet p T .

Jet reconstruction and corrections
The CMS particle-flow (PF) algorithm [38,39] identifies stable particles in an event by combining information from all sub-detector systems, classifying them as electrons, muons, photons, and charged and neutral hadrons. The PF candidates are then clustered into jets using the anti-k T sequential recombination algorithm [37] provided in the FastJet framework [40]. The results in this analysis are obtained using a distance parameter R = 0.3. The underlying event (UE) contribution to the jet energy is subtracted using an iterative procedure described in Refs. [10,41]. The jet energies are then corrected to contain the energy of all final-state jet constituents as described in Ref. [42]. The jet energy corrections are derived using simulated pythia (6.462, Z2 tune) [43,44] events and measurements of the energy balance of dijet and photon + jet pPb collision events are used to correct differences between data and Monte Carlo (MC) distributions [23,42]. In the jet reconstruction process, there is a possibility that the jet energy is estimated incorrectly, or a jet is found in a region where the UE has an upward fluctuation, but no hard scattering has occurred (a "fake" jet). In MC the "real" and "fake" jets can be distinguished by requiring that the reconstructed jet is matched to a generator-level jet. In data, this cannot be done directly, but the contribution of fake jets could be estimated from MC, provided that it is tuned to describe  Top the weighted jet spectra using prescale factors from each HLT-triggered event sample and the combined jet spectrum. A subset of the data is plotted to illustrate the procedure. Bottom the ratios of each individual HLT-triggered jet spectrum to the combined jet spectrum. Statistical uncertainties are shown as vertical bars, and p T bin widths as horizontal bars the data, and specific jet selections are developed to identify and remove the misreconstructed jets. We estimate that about 10 % of the jets reconstructed at p T = 50 GeV/c in pPb collisions are fake, and this fraction quickly drops to a level of 10 −4 at p T ≈ 100 GeV/c. After the jet-identification cuts are applied, we estimate that less than 1 % fake jets remain in the sample.
Because of the finite detector resolution and the steeply falling p T distributions, the measured jet p T spectra are smeared with respect to the true distributions, although the mean value of the reconstructed jet energy is corrected as described above. The jet energy resolution is estimated to be 13 % (8 %) for jet p T = 60 (300) GeV/c. A Bayesian unfold-ing technique [45] is employed to account for such resolution effects, as implemented in the RooUnfold package [46]. The migration of jets in pseudorapidity is not explicitly corrected for; it is instead included as an uncertainty, as discussed in Sect. 3.1. In the unfolding method, a response matrix is built based on MC simulations and is used to obtain the "true" jet p T distribution from the measured one. Jets are first generated with the pythia event generator and then embedded into pPb collisions simulated with the hijing event generator (version 1.383) [47], which have particle multiplicity distributions comparable to the pPb data and can account for additional resolution effects associated with the higher detector occupancy. These embedded MC samples are denoted hereafter by pythia + hijing. The unfolding technique is tested by building the response matrix with detector jets (Reco) and generated jets (Gen) from half of the MC sample and applying it to unfold the other half of the sample. The top panel of Fig. 2 shows the response matrix obtained using the pythia + hijing simulation, while the bottom panel shows the ratio of the jet spectrum reconstructed from the simulation after unfolding to the generator-level jet spectrum. The unfolded MC jet spectrum is compatible with the generatorlevel jet spectrum within the statistical uncertainties. The results reported in this paper are based on the Bayesian unfolding technique that uses four iteration steps. Up to eight iteration steps are used in evaluating the systematic uncertainties as discussed in Sect. 3.1. The generator level pythia jet spectrum is used as a prior in the unfolding. The data points are reported in the center of each p T bin without corrections for binning effects.
The pPb jet cross sections are obtained in several pseudorapidity intervals. To study the evolution of the jet cross section with pseudorapidity, ratios of jet spectra are computed either using symmetric positive and negative pseudorapidity intervals around mid-rapidity, or normalizing the distributions by the mid-rapidity jet spectrum. These ratios are taken in the same p T bin and the values are reported at the center of the bin. To study nuclear effects on jet production, the jet spectra in pPb collisions are compared to pp reference spectra obtained by extrapolation from previous jet cross section measurements in pp collisions at higher center-of-mass energy. The nuclear modification factor, R pPb , evaluated in several pseudorapidity intervals, is defined as where L = 30.1 nb −1 is the effective integrated luminosity in the pPb analysis, corrected for event-selection efficiency and trigger prescales, and A is the mass number of the lead nucleus. Since presently there are no available experimental results from pp collisions at √ s = 5.02 TeV, for this paper Bottom the ratios of the Bayesian unfolded jet p T spectrum reconstructed in the simulation and the generator-level spectrum we use extrapolated, rather than measured, pp reference spectra. Hence we denote the nuclear modification factors as R * pPb .

Proton-proton reference jet spectra
The reference pp spectra are constructed extrapolating previously published inclusive jet spectra measured in pp collisions at √ s = 7 TeV. Measurements performed with the anti-k T jet algorithm with two distance parameters, R = 0.5 and 0.7 [33], are used in the extrapolation. The extrapolation is based on the pythia generator (6.462, Z2 tune) and is performed in two steps. First, the √ s = 7 TeV jet cross section measurements are extrapolated to √ s = 5.02 TeV and then scaled to R = 0.3, since a smaller distance parameter is used in the pPb analysis to minimize the UE background fluctuations. The pythia generator is used to estimate p T -dependent scaling factors. While this scaling is model dependent, the There are several sources of systematic uncertainty in the measurements of the jet spectra, the jet yield asymmetry, and the nuclear modification factors R * pPb . The dominant uncertainties in the spectra measured in pPb collisions come from the unfolding of the spectra and from the jet energy scale (JES) corrections, which are partially correlated since they both aim to correct for the difference between the reconstructed and the true jet energy. The stability of the unfolding procedure and its ability to recover the generator-level jet spectrum have been verified with simulation studies, which included the use of different numbers of iterations (n = 2, 3, 4, 5, 6, 7, 8). In the data, the unfolded spectra for n = 4 are compared to the spectra obtained with different values of n and the difference is included in the systematic uncer-tainty. In addition, since the true jet spectrum may differ in shape from the spectrum in the MC generator, the slope of the prior guess distribution is varied such that the yield at low p T increases or decreases by a factor of 3, while at high p T the yield is changed by about 10-20 %. After this variation the spectra are unfolded and then are compared to the nominal unfolded spectra to estimate the uncertainty due to the nominal input distribution. The uncertainties from unfolding are largest (up to 5 %) in the low p T region and at large absolute pseudorapidity. Uncertainties that arise from the different jet energy resolution in the data and MC simulation are evaluated by smearing the unfolding matrix to account for these differences and then redoing the unfolding. The resulting differences in the final jet spectra are found to be less than 1 %. The JES uncertainty is about 1 % and induces up to 7 % changes at high p T because of the steeply falling jet spectra.
Additional cross checks are performed comparing the spectra obtained with different jet reconstruction algorithms (such as subtracting the UE in the jet algorithm or correcting for it in the transfer matrix), and comparing the unfolded results when the unfolding matrix uses the reconstructed jet p T with or without jet energy corrections. The total uncertainty in the jet spectra due to the JES and unfolding varies from about 5 % at low jet p T at mid-rapidity to about 10 % for high p T and forward rapidity.
The fake jet contribution is estimated on the basis of a MC study of various jet quality variables that are used to identify genuine and misreconstructed jet contributions. In the pythia + hijing embedded samples these variables are optimized to remove misreconstructed jets, while preserving the largest fraction of genuine jets. The uncertainty in the misreconstructed jet contribution in the jet spectra is estimated by varying the jet quality requirements and comparing the resulting spectra in data and in simulation. It is about 1 % for all pseudorapidity ranges.
The unfolding procedure does not correct for possible misreconstruction of the jet axis, and therefore jets may migrate from one pseudorapidity interval to another thus altering the jet spectra measured in different η ranges. The uncertainty associated with the jet pointing resolution is estimated by building the unfolding matrix using either the generated or the reconstructed jet axis, and comparing the resulting unfolded jet spectra. This uncertainty is found to be of the order of 1 % in the central pseudorapidity region and 2 % at large absolute pseudorapidity.
The jet spectra in pPb collisions are also subject to an overall scale uncertainty, due to the uncertainties in the integrated luminosity measurement. The scale uncertainty is estimated to be 3.5 %, as described in Ref. [48].
The systematic uncertainty in the inclusive jet production asymmetry only includes those factors that depend on the jet pseudorapidity, such as the JES, unfolding, and misreconstructed jet contribution uncertainties. The overall scale uncertainty due to the luminosity normalization cancels out. As a cross check, the jet yield asymmetry uncertainties are evaluated using a combination of the two data sets with different beam directions. In that case, the jet yield asymmetry can be measured using detector elements that are only in the positive η or in the negative η ranges in the laboratory frame. Since the detector is symmetric, these regions have similar acceptance and performance and we expect that systematic effects are also similar. Alternatively, the jet yield asymmetry is measured from each portion of the data independently, and the results of this comparison confirm the systematic uncertainty estimate obtained by evaluating each source of uncertainty separately.

Systematic uncertainties in the pp reference
The uncertainties in the extrapolated pp reference spectra take into account the uncertainties in the distance parameter dependence of the cross sections at √ s = 7 TeV and the scaling to the smaller R = 0.3 value, the uncertainty in the √ s dependent scaling, as well as the uncertainties of the input spectra used in the extrapolation. The uncertainties in the inclusive jet measurements from pp collisions at √ s = 7 TeV reported in Ref.
[33] are taken as the upper and lower limits of the cross sections used in the extrapolation, and are reflected in the uncertainties of the resulting reference spectra. The following alternative approaches are used to derive scaling factors and evaluate their uncertainties.
1. pythia 8, CUETP8M1 tune [49,50]: In the kinematic range studied, this tune has a different quark-to-gluon jet ratio and different jet shapes than the pythia 6, Z2 tune used for the nominal result. 2. powheg + pythia event generator [51,52]: The powheg generator is used to compute the cross section at nextto-leading order (NLO) accuracy, and pythia (6.462, Z2 tune) is used to describe the parton showering and hadronization. 3. NLO calculations [53,54] with several different parametrizations of the parton distribution functions [55] and non-perturbative corrections based on pythia ( The jet cross sections for R = 0.3 and R = 0.5 at √ s = 5.02 TeV are evaluated using (1), (2), and (3). Then the ratios between the cross sections obtained with these two distance parameters, in the default pythia calculation (6.462, Z2 tune) and in the alternative methods, are compared to each other, leading to an uncertainty in the distance parameter scaling of around 5 %. The √ s scaling factors are evaluated with (2) and (3) for R = 0.5, and with (2), (3), and (4) for R = 0.7. These scaling factors are compared to the results from pythia (6.462, Z2 tune). The uncertainties in the √ s scaling factors range from 4 % at low jet p T in the mid-rapidity region to 7 % at high p T and at forward rapidity. The total uncertainty in the pp reference extrapolation is found to range between 9 % at mid-rapidity and 11 % at forward rapidity. These uncertainties include a 2.4 % scale uncertainty from the integrated luminosity measurement [33].

Summary of systematic uncertainties
A summary of the systematic uncertainties in the jet spectra in pPb collisions, the jet yield asymmetry measurements in pPb collisions, the reference pp spectra, and the nuclear modification factors R * pPb are listed in Table 1. The uncertainties depend on the jet p T and pseudorapidity, and the table shows representative values in two jet p T and η CM ranges. The uncertainties vary smoothly between these ranges. The total systematic uncertainties listed for the nuclear modification factors R * pPb do not include the scale uncertainty of 4.3 % from the integrated luminosity measurements in pPb (3.5 %) and pp (2.4 %) collisions. The luminosity uncertainties cancel in the measurements of the jet yield asymmetry. The remaining uncertainties are partially correlated in jet p T , with the unfolding uncertainty dominating at low jet p T and the JES uncertainty dominating at high jet p T .

Results and discussion
The inclusive jet differential cross sections in pPb collisions at √ s NN = 5.02 TeV are shown in Fig. 4 for six consecutive η intervals in the range −2.0 < η CM < 1.5 and the range |η CM | < 1.0 for reference purposes. The distributions are scaled by arbitrary factors described in the legend to enhance visibility. These spectra are used to study the pseudorapidity dependence of inclusive jet production in pPb collisions and possible nuclear effects. In symmetric collisions, such as in the pp system, the kinematic range in the fractional momentum x probed with the jets in forward and backward pseudorapidity is the same and the production is symmetric about η CM = 0. In the pPb system, the jets produced at forward pseudorapidity (proton beam direction) correlate with smaller x values from the Pb nucleus than those produced at backward pseudorapidity. Based on a generatorlevel study made with pythia, the average x values from the Pb nucleus (Fig. 5) that are probed in the kinematic range covered by the present measurement are estimated to be in the range 0.03 x Pb 0.5. Values of p T that correspond to x Pb 0.2 are associated with anti-shadowing in the nPDFs. The region x Pb 0.2 is associated with a Table 1 Systematic uncertainties in the measurement of the jet spectra in pPb collisions are shown in the first four lines. The sources and corresponding systematic uncertainties in the extrapolated pp reference are presented in the next four lines. The total uncertainties in the jet spectra in pPb collisions, the reference pp spectra, the jet yield asymmetry in pPb collisions, and R * pPb are shown in the bottom four lines. The uncer-tainties depend on the jet p T and pseudorapidity, and the table shows representative values in two jet p T and η CM ranges. The uncertainties vary smoothly between these two ranges. Total systematic uncertainties listed for the nuclear modification factors R * pPb do not include the scale uncertainty of 4.3 % due to the uncertainty in the integrated luminosity measurements in pPb (3.5 %) and pp (2.4 %) collisions

Source
Jet p T < 80 GeV/c Jet p T > 150 GeV/c  The forward-backward asymmetry of the jet production is evaluated by taking the ratio between the jet yields in the Pbgoing and the proton-going directions for two pseudorapidity  The evolution of the jet spectra with pseudorapidity can also be studied by normalizing each spectrum to the one obtained in the mid-rapidity range (|η CM | < 1). The normalized jet cross section distributions are shown in the top panel of Fig. 7. In the bottom panel of Fig. 7 we examine the pseudorapidity dependence in the normalized jet cross sections in three fixed p T bins. The data points are offset for visibility. No significant pseudorapidity asymmetry is observed as can also be seen by comparing the open and closed stars or open and closed crosses in the top panel. The jet spectra become softer away from the mid-rapidity region, and the pseudorapidity distributions become narrower with increasing jet p T as a result of the softening of the distributions at forward and backward pseudorapidity.
The inclusive jet nuclear modification factors R * pPb as a function of jet p T are shown in Fig. 8 for six center-of-mass pseudorapidity bins, along with an NLO perturbative QCD (pQCD) calculation [57] using the EPS09 nPDFs [19]. For most of the measured p T and η CM ranges, the experimental R * pPb values are systematically above the theoretical prediction. However, this difference is not significant, given the size of the systematic uncertainties and the fact that they are strongly correlated in p T . The R * pPb values are approximately independent of p T . In the theoretical prediction there is a decrease in R pPb with p T in the backward pseudorapidity region, which is associated with the onset of the EMC effect at high values of x in the Pb nucleus. In the range of p T where the measurements probe the anti-shadowing region, the R * pPb values show a hint of an enhancement with respect to the pp reference, e.g. for |η CM | < 0.5 and 56 < p T < 300 GeV/c, R * pPb = 1.17 ± 0.01 (stat) ± 0.12 (syst). This enhancement  The results of the jet R * pPb measurements presented here are consistent with those reported by the ATLAS collaboration [22]. In Fig. 9 we compare our results to the ATLAS measurement at mid-rapidity, |y CM | < 0.3, for the 0-90 %

Summary
The inclusive jet spectra and nuclear modification factors in pPb collisions at √ s NN = 5.02 TeV have been measured. The data, corresponding to an integrated luminosity of 30.1 nb −1 , were collected by the CMS experiment in 2013. The jet transverse momentum spectra were measured for p T > 56 GeV/c in six pseudorapidity intervals covering the range −2 < η CM < 1.5 in the NN center-of-mass system. The jet spectra were found to be softer away from mid-rapidity. The jet production at forward and backward pseudorapidity were compared, and no significant asymmetry about η CM = 0 was observed in the measured kinematic range.
The differential jet cross section results were compared with extrapolated pp reference spectra based on jet measurements in pp collisions at √ s = 7 TeV. The inclusive jet nuclear modification factors R * pPb were observed to have small enhancements compared to the reference pp jet spectra at low jet p T in all η CM ranges. In the anti-shadowing region, for |η CM | < 0.5 and 56 < p T < 300 GeV/c, the value R * pPb = 1.17 ± 0.01 (stat) ± 0.12 (syst) was found. The R * pPb appears to be approximately independent of p T , except in the most backward pseudorapidity range. The R * pPb measurements were found to be compatible with theoretical predictions from NLO pQCD calculations that use EPS09 nPDFs.
Acknowledgments We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centers and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: