Measurement of inclusive jet and dijet production in pp collisions at sqrt(s) = 7 TeV using the ATLAS detector

Inclusive jet and dijet cross sections have been measured in proton-proton collisions at a centre-of-mass energy of 7 TeV using the ATLAS detector at the Large Hadron Collider. The cross sections were measured using jets clustered with the anti-kT algorithm with parameters R=0.4 and R=0.6. These measurements are based on the 2010 data sample, consisting of a total integrated luminosity of 37 inverse picobarns. Inclusive jet double-differential cross sections are presented as a function of jet transverse momentum, in bins of jet rapidity. Dijet double-differential cross sections are studied as a function of the dijet invariant mass, in bins of half the rapidity separation of the two leading jets. The measurements are performed in the jet rapidity range |y|<4.4, covering jet transverse momenta from 20 GeV to 1.5 TeV and dijet invariant masses from 70 GeV to 5 TeV. The data are compared to expectations based on next-to-leading order QCD calculations corrected for non-perturbative effects, as well as to next-to-leading order Monte Carlo predictions. In addition to a test of the theory in a new kinematic regime, the data also provide sensitivity to parton distribution functions in a region where they are currently not well-constrained.

The ATLAS Collaboration has published a first measurement of inclusive jet and dijet production at √ s = 7 TeV, using an integrated luminosity of 17 nb −1 [21]. This measurement considered only jets with transverse momentum larger than 60 GeV and in a rapidity interval |y| < 2.8 1 .
The analysis presented here extends the previous measurement using the 2010 data sample of (37.3±1.2) pb −1 , an integrated luminosity more than 2000 times larger than that of the previous study. This more than doubles the kinematic reach at high jet transverse momentum and large dijet invariant mass. There are strong physics reasons to extend the measurement to jets of 1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y axis points upward. Cylindrical coordinates (r, φ) are used in the transverse plane, φ being the azimuthal angle around the beam pipe, referred to the x-axis. The pseudorapidity is defined in terms of the polar angle θ with respect to the beamline as η = − ln tan(θ/2). When dealing with massive jets and particles, the rapidity y = 1 2 ln E+pz E−pz is used, where E is the jet energy and pz is the z-component of the jet momentum.
lower transverse momentum and larger rapidity as well. Jets at lower p T are more sensitive to non-perturbative effects from hadronisation and the underlying event, and forward jets may be sensitive to different dynamics in QCD than central jets. Moreover, LHC experiments have much wider rapidity coverage than those at the Tevatron, so forward jet measurements at the LHC cover a phase space region that has not been explored before. The kinematic reach of this analysis is compared to that of the previous ATLAS study in Fig. 1. This data sample extends the existing inclusive jet p T measurement from 700 GeV to 1.5 TeV and the existing dijet mass measurement from 1.8 TeV to 5 TeV. Thus this analysis probes next-to-leading order (NLO) perturbative QCD (pQCD) and parton distribution functions (PDFs) in a new kinematic regime. The results span approximately 7 × 10 −5 < x < 0.9 in x, the fraction of the proton momentum carried by each of the partons involved in the hard interaction.

II. THE ATLAS DETECTOR
The ATLAS detector is described in detail in Ref. [22]. In this analysis, the tracking detectors are used to define candidate collision events by constructing vertices from tracks, and the calorimeters are used to reconstruct jets.
The inner detector used for tracking and particle identification has complete azimuthal coverage and spans the region |η| < 2.5. It consists of layers of silicon pixel detectors, silicon microstrip detectors, and transition radiation tracking detectors, surrounded by a solenoid magnet that provides a uniform field of 2 T.
The electromagnetic calorimetry is provided by the liquid argon (LAr) calorimeters that are split into three arXiv:1112.6297v2 [hep-ex] 16 Apr 2012 Kinematic reach of the inclusive jet cross section measured in this analysis compared to that of the previous study [21] for jets identified using the anti-kt algorithm with R = 0.6. The kinematic limit for the center-of-mass energy of 7 TeV is also shown.

III. CROSS SECTION DEFINITION
Jet cross sections are defined using the anti-k t jet algorithm [23] implemented in the FASTJET [24] package. Two different values are used for the clustering parameter R (0.4 and 0.6), which can be seen intuitively as the radius of a circular jet in the plane (φ, y) of azimuthal angle and rapidity. The jet cross section measurements are corrected for all experimental effects, and thus are defined for the "particle-level" final state of a proton-proton collision [25]. Particle-level jets in the Monte Carlo simulation are identified using the anti-k t algorithm and are built from stable particles, which are defined as those with a proper lifetime longer than 10 ps. This definition includes muons and neutrinos from decaying hadrons.
Inclusive jet double-differential cross sections are measured as a function of jet p T in bins of y, in the region p T > 20 GeV, |y| < 4.4. The term "inclusive jets" is used in this paper to indicate that all jets in each event are considered in the cross section measurement. Dijet double-differential cross sections are measured as a function of the invariant mass of the two leading (highest p T ) jets, which is given as m 12 = (E 1 + E 2 ) 2 − ( p 1 + p 2 ) 2 , where E 1,2 and p 1,2 are the energies and momenta of the two leading jets. The cross sections are binned in the variable y * , defined as half the rapidity difference of the two leading jets, y * = |y 1 − y 2 |/2. The quantity y * is the rapidity in the two-parton centre-of-mass frame (in the massless particle limit), where it is determined by the polar scattering angle with respect to the beamline, θ * : For the dijet measurement, the two leading jets are selected to lie in the |y| < 4.4 region, where the leading jet is required to have p T > 30 GeV and the sub-leading jet p T > 20 GeV. Restricting the leading jet to higher p T improves the stability of the NLO calculation [26]. Theory calculations are used in the same kinematic range as the measurement.

IV. MONTE CARLO SAMPLES
The PYTHIA 6.423 generator [27] with the MRST LO* PDF set [28] was used to simulate jet events in proton-proton collisions at a centre-of-mass energy of √ s = 7 TeV and to correct for detector effects. This generator utilizes leading-order perturbative QCD matrix elements (ME) for 2 → 2 processes, along with a leading-logarithmic, p T -ordered parton shower (PS), an underlying event simulation with multiple parton interactions, and the Lund string model for hadronisation. Samples were generated using the ATLAS Minimum Bias Tune 1 (AMBT1) set of parameters [29], in which the model of non-diffractive scattering has been tuned to ATLAS measurements of charged particle production at √ s = 900 GeV and √ s = 7 TeV. The particle four-vectors from these generators were passed through a full simulation [30] of the ATLAS detector and trigger that is based on GEANT4 [31]. Finally, the simulated events were reconstructed and jets were calibrated using the same reconstruction chain as the data.

V. THEORETICAL PREDICTIONS
A. Fixed-Order calculations

NLO Predictions
The measured jet cross sections are compared to fixed-order NLO pQCD predictions, with corrections for non-perturbative effects applied. For the hard scattering, both the NLOJET++ 4.1. 2 [32] package and the POWHEG generator [33,34] were used, the latter in a specific configuration where the parton shower was switched off and calculations were performed using NLO matrix elements. The two programs have been used with the CT10 [35] NLO parton distribution functions, and the same value of normalisation and factorisation scale, corresponding to the transverse momentum of the leading jet, p max T : For POWHEG, p max T is evaluated at leading order and is denoted p Born T . Using this scale choice, the cross section results of the two NLO codes are compatible at the few percent level for inclusive jets over the whole rapidity region. They are also consistent for dijet events where both jets are in the central region, while they differ substantially when the two leading jets are widely separated in rapidity (y * 3). In these regions, NLO-JET++ gives an unstable and much smaller cross section than POWHEG that is even negative for some rapidity separations. POWHEG remains positive over the whole region of phase space. It should be noted that the forward dijet cross section predicted by NLOJET++ in this region has a very strong scale dependence, which however is much reduced for larger values of scale than that of Eq. 2.
The forward dijet cross section for NLOJET++ is much more stable if instead of a scale fixed entirely by p T , a scale that depends on the rapidity separation between the two jets is used. The values chosen for each y * -bin follow the formula: and are indicated by the histogram in Fig. 2. These values are motivated by the formula (shown by the dotdashed curve): that is suggested in Ref. [36], and are in a region where the cross section predictions are more stable as a function of scale (they reach a "plateau"). At small y * , the scale in Eq. 3 reduces to the leading jet p T (dotted line), which is used for the inclusive jet predictions. With this scale choice, NLOJET++ is again in reasonable agreement with POWHEG, which uses the scale from Eq. 2. The NLOJET++ predictions are used as a baseline for both inclusive jet and dijet calculations, with the scale choice from Eq. 2 for the former and that from Eq. 3 for the latter. The POWHEG scale used for both inclusive jets and dijets, p Born T , is given by Eq. 2 but evaluated at leading order. Despite using different scale choices, the dijet theory predictions from NLOJET++ and POWHEG are stable with respect to relatively small scale variations and give consistent results.
The main uncertainties on the NLO prediction come from the uncertainties on the PDFs, the choice of factorisation and renormalisation scales, and the uncertainty on used for the dijet predictions obtained using NLOJET++, as a function of y * , half the rapidity separation between the two leading jets. This is motivated by the scale choice suggested in Ref. [36] (dot-dashed line), and is also compared to the scale choice used for the inclusive jet predictions (dotted line).
the value of the strong coupling constant α s . To allow for fast and flexible evaluation of PDF and scale uncertainties, the APPLGRID [41] software was interfaced with NLOJET++ in order to calculate the perturbative coefficients once and store them in a look-up table. The PDF uncertainties are defined at 68% CL and evaluated following the prescriptions given for each PDF set. They account for the data uncertainties, tension between input data sets, parametrisation uncertainties, and various theoretical uncertainties related to PDF determination.
To estimate the uncertainty on the NLO prediction due to neglected higher-order terms, each observable was recalculated while varying the renormalisation scale by a factor of two with respect to the default choice. Similarly, to estimate the sensitivity to the choice of scale where the PDF evolution is separated from the matrix element, the factorisation scale was separately varied by a factor of two. Cases where the two scales are simultaneously varied by a factor 2 in opposite directions were not considered due to the presence of logarithmic factors in the theory calculation that become large in these configurations. The envelope of the variation of the observables was taken as a systematic uncertainty. The effect of the uncertainty on the value of the strong coupling constant, α s , is evaluated following the recommendation of the CTEQ group [42], in particular by using different PDF sets that were derived using the positive and negative variations of the coupling from its best estimate.
Electro-weak corrections were not included in the theory predictions and may be non-negligible [43]. Non-perturbative correction factors for inclusive jets identified using the anti-kt algorithm with distance parameters R = 0.4 and R = 0.6 in the rapidity region |y| < 0.3, derived using various Monte Carlo generators. The correction derived using PYTHIA 6.425 with the AUET2B CTEQ6L1 tune is used for the fixed-order NLO calculations presented in this analysis.

Non-Perturbative Corrections
The fixed-order NLO calculations predict partonlevel cross sections, which must be corrected for nonperturbative effects to be compared with data. This is done by using leading-logarithmic parton shower generators. The corrections are derived by using PYTHIA 6.425 with the AUET2B CTEQ6L1 tune [44] to evaluate the bin-wise ratio of cross sections with and without hadronisation and the underlying event. Each bin of the parton-level cross section is then multiplied by the corresponding correction. The uncertainty is estimated as the maximum spread of the correction factors obtained from PYTHIA 6.425 using the AUET2B LO * * , AUET2 LO * * , AMBT2B CTEQ6L1, AMBT1, Perugia 2010, and Perugia 2011 tunes (PYTUNE 350), and the PYTHIA 8.150 tune 4C [44][45][46][47], as well as those obtained from the HERWIG++ 2.5.1 [48] tune UE7000-2 [44]. The AMBT2B CTEQ6L1 and AMBT1 tunes, which are based on observables sensitive to the modeling of minimum bias interactions, are included to provide a systematically different estimate of the underlying event activity.
The corrections depend strongly on the jet size; therefore separate sets of corrections and uncertainties were derived for jets with R = 0.4 and R = 0.6. The correction factors and their uncertainties depend on the interplay of the hadronisation and the underlying event for the dif-ferent jet sizes, and they have a significant influence at low p T and low dijet mass. For R = 0.4, the correction factors are dominated by the effect of hadronisation and are approximately 0.95 at jet p T = 20 GeV, increasing closer to unity at higher p T . For R = 0.6, the correction factors are dominated by the underlying event and are approximately 1.6 at jet p T = 20 GeV, decreasing to between 1.0-1.1 for jets above p T = 100 GeV. Fig. 3 shows the non-perturbative corrections for inclusive jets with rapidity in the interval |y| < 0.3, for jet clustering parameters R = 0.4 and R = 0.6. The correction factors for the other rapidity bins become closer to unity as the jet rapidity increases, as can be seen in Fig. 20 in Appendix A.
Non-perturbative corrections have been evaluated for the dijet measurement as well, as a function of the dijet mass and the rapidity interval y * , for each of the two jet sizes. These follow a similar behaviour to those for inclusive jets, with the corrections becoming smaller for large invariant masses and rapidity differences.

B. NLO Matrix Element + Parton Shower
The measured jet cross sections are also compared to POWHEG [49], an NLO parton shower Monte Carlo generator that has only recently become available for inclusive jet and dijet production. POWHEG uses the POWHEG BOX package [50][51][52] and allows one to use either PYTHIA or HERWIG [53] + JIMMY [54] to shower the partons, hadronise them, and model the underlying event. The ATLAS underlying event tunes, AUET2B for PYTHIA and AUET2 [55] for HERWIG, are derived from the standalone versions of these event generators, with no optimisation for the POWHEG predictions. The showering portion of POWHEG uses the PDFs from PYTHIA or HERWIG as part of the specific tune chosen.
In the POWHEG algorithm, each event is built by first producing a QCD 2 → 2 partonic scattering. The renormalisation and factorisation scales are set to be equal to the transverse momentum of the outgoing partons, p Born T , before proceeding to generate the hardest partonic emission in the event. 2 The CT10 NLO PDF set is used in this step of the simulation. Then the event is evolved to the hadron level using a parton shower event generator, where the radiative emissions in the parton showers are required to be softer than the hardest partonic emission generated by POWHEG.
The coherent simulation of the parton showering, hadronisation, and the underlying event with the NLO matrix element is expected to produce a more accurate theoretical prediction. In particular, the nonperturbative effects are modeled in the NLO parton shower simulation itself, rather than being derived separately using a LO parton shower Monte Carlo generator as described in Sec. V A 2.

A. Dataset
The inclusive jet and dijet cross section measurements use the full ATLAS 2010 data sample from proton-proton collisions at √ s = 7 TeV. For low-p T jets, only the first 17 nb −1 of data taken are considered since the instantaneous luminosity of the accelerator was low enough that a large data sample triggered with a minimum bias trigger (see Sec. VI B) could be recorded. This provides an unbiased sample for reconstructing jets with p T between 20-60 GeV, below the lowest jet trigger threshold. In addition, during this period there were negligible contributions from "pile-up" events, in which there are multiple proton-proton interactions during the same or neighbouring bunch crossings. Thus this period provides a well-measured sample of lowp T jets. The first data taking period was not used for forward jets with |y| > 2.8 and p T > 60 GeV because the forward jet trigger was not yet commissioned.
For all events considered in this analysis, good operation status was required for the first-level trigger, the solenoid magnet, the inner detector, the calorimeters and the luminosity detectors, as well as for tracking and jet reconstruction. In addition, stable operation was required for the high-level trigger during the periods when this system was used for event rejection.

B. Trigger
Three different triggers have been used in this measurement: the minimum bias trigger scintillators (MBTS); the central jet trigger, covering |η| < 3.2; and the forward jet trigger, spanning 3.1 < |η| < 4.9. The MBTS trigger requires at least one hit in the minimum bias scintillators located in front of the endcap cryostats, covering 2.09 < |η| < 3.84, and is the primary trigger used to select minimum-bias events in ATLAS. It has been demonstrated to have negligible inefficiency for the events of interest for this analysis [56] and is used to select events with jets having transverse momenta in the range 20-60 GeV. The central and forward jet triggers are com-posed of three consecutive levels: Level 1 (L1), Level 2 (L2) and Event Filter (EF). In 2010, only L1 information was used to select events in the first 3 pb −1 of data taken, while both the L1 and L2 stages were used for the rest of the data sample. The jet trigger did not reject events at the EF stage in 2010.
The central and forward jet triggers independently select data using several thresholds for the jet transverse energy (E T ≡ E sin θ), each of which requires the presence of a jet with sufficient E T at the electromagnetic (EM) scale. 3 For each L1 threshold, there is a corresponding L2 threshold that is generally 15 GeV above the L1 value. Each such L1+L2 combination is referred to as an L2 trigger chain. Fig. 4 shows the efficiency for L2 jet trigger chains with various thresholds as a function of the reconstructed jet p T for anti-k t jets with R = 0.6 for both the central and forward jet triggers. Similar efficiencies are found for jets with R = 0.4, such that the same correspondence between transverse momentum regions and trigger chains can be used for the two jet sizes. The highest trigger chain does not apply a threshold at L2, so its L1 threshold is listed.
As the instantaneous luminosity increased throughout 2010, it was necessary to prescale triggers with lower E T thresholds, while the central jet trigger with the highest E T threshold remained unprescaled. As a result, the vast majority of the events where the leading jet has transverse momentum smaller than about 100 GeV have been taken in the first period of data-taking, under conditions with a low amount of pile-up, while the majority of the high-p T events have been taken during the second datataking period, with an average of 2-3 interactions per bunch crossing. For each p T -bin considered in this analysis, a dedicated trigger chain is chosen that is fully efficient (> 99%) while having as small a prescale factor as possible. For inclusive jets fully contained in the central or in the forward trigger region, only events taken by this fully efficient trigger are considered. For inclusive jets in the HEC-FCal transition region 2.8 ≤ |y| < 3.6, neither the central nor the forward trigger is fully efficient. Instead, the logical OR of the triggers is used, which is fully efficient at sufficiently high jet p T (see Fig. 5

).
A specific strategy is used to account for the various prescale combinations for inclusive jets in the HEC-FCal transition region, which can be accepted either by the central jet trigger only, by the forward jet trigger only, or by both. A similar strategy is used for dijet events in a given (m 12 , y * )-bin, which can be accepted by several jet triggers depending on the transverse momenta and pseudorapidities of the two leading jets. Events that can 3 The electromagnetic scale is the basic calorimeter signal scale for the ATLAS calorimeters. It has been established using testbeam measurements for electrons and muons to give the correct response for the energy deposited in electromagnetic showers, while it does not correct for the lower response of the calorimeter to hadrons.
[GeV] FIG. 4. Combined L1+L2 jet trigger efficiency as a function of reconstructed jet pT for anti-kt jets with R = 0.6 in the central region |y| < 0.3 (a), the barrel-endcap transition region 1.2 ≤ |y| < 2.1 (b) and the FCal region 3.6 ≤ |y| < 4.4 (c) for the different L2 trigger thresholds used in the analysis. The trigger thresholds are at the electromagnetic scale, while the jet pT is at the calibrated scale (see Sec. VI C). The highest trigger chain used for |y| < 2.8 does not apply a threshold at L2, so its L1 threshold is listed. The efficiency in the |y| > 3.2 rapidity range is not expected to reach 100% due to the presence of a dead FCal trigger tower that spans 0.9% of the (η, φ)-acceptance. This inefficiency is assigned as a systematic uncertainty on the trigger efficiency in the measurement. Efficiencies for the central and forward jet triggers with a L1 ET threshold of 10 GeV, and for their logical OR, as a function of the rapidity y of the reconstructed jet in the transition region between the two trigger systems. The logical OR is used for the inclusive jet measurement to collect data in the 2.8 ≤ |y| < 3.6 rapidity slice. be accepted by more than one trigger chain have been divided into several categories according to the trigger combination that could have accepted the events. For inclusive jets in the transition region, these correspond to central and forward triggers with a similar threshold; for dijets the trigger combination depends on the position and transverse momenta of the two leading jets, each of which is "matched" to a trigger object using angular criteria. Corrections are applied for any trigger ineffi-ciencies, which are generally below 1%. The equivalent luminosity of each of the categories of events is computed based on the prescale values of these triggers throughout the data-taking periods, and all results from the various trigger combinations are combined together according to the prescription given in Ref. [57].

C. Jet Reconstruction and Calibration
Jets are reconstructed at the electromagnetic scale using the anti-k t algorithm. The input objects to the jet algorithm are three-dimensional topological clusters [58] built from calorimeter cells. The four-momentum of the uncalibrated, EM-scale jet is defined as the sum of the four-momenta of its constituent calorimeter energy clusters. Additional energy due to multiple proton-proton interactions within the same bunch crossing ("pile-up") is subtracted by applying a correction derived as a function of the number of reconstructed vertices in the event using minimum bias data. The energy and the position of the jet are next corrected for instrumental effects such as dead material and non-compensation. This jet energy scale (JES) correction is calculated using isolated jets 4 in the Monte Carlo simulation as a function of energy and pseudorapidity of the reconstructed jet. The JES correction factor ranges from about 2.1 for low-energy jets with p T = 20 GeV in the central region |y| < 0.3 to less than 1.2 for high-energy jets in the most forward region 3.6 ≤ |y| < 4.4. The corrections are cross-checked using in-situ techniques in collision data (see below) [59].

D. Uncertainties in Jet Calibration
The uncertainty on the jet energy scale is the dominant uncertainty for the inclusive jet and dijet cross section measurements. Compared to the previous analysis [21], this uncertainty has been reduced by up to a factor of two, primarily due to the improved calibration of the calorimeter electromagnetic energy scale obtained from Z → ee events [60], as well as an improved determination of the single particle energy measurement uncertainties from in-situ and test-beam measurements [61]. This improvement is confirmed by independent measurements, including studies of the momenta of tracks associated to jets, as well as the momentum balance observed in γ+jet, dijet, and multijet events [59].
In the central barrel region (|η| < 0.8), the dominant source of the JES uncertainty is the knowledge of the calorimeter response to hadrons. This uncertainty is obtained by measuring the response to single hadrons using proton-proton and test-beam data, and propagating the uncertainties to the response for jets. Additional uncertainties are evaluated by studying the impact on the calorimeter response from varying settings for the parton shower, hadronization, and underlying event in the Monte Carlo simulation. The estimate of the uncertainty is extended from the central calorimeter region to the endcap and forward regions, the latter of which lies outside the tracking acceptance, by exploiting the transverse momentum balance between a central and a forward jet in events where only two jets are produced.
In the central region (|η| < 0.8), the uncertainty is lower than 4.6% for all jets with p T > 20 GeV, which decreases to less than 2.5% uncertainty for jet transverse momenta between 60 and 800 GeV. The JES uncertainty is the largest for low-p T (∼20 GeV) jets in the most forward region |η| > 3.6, where it is about 11-12%. Details of the JES determination and its uncertainty are given in Ref. [59].

Event Selection
To reject events due to cosmic-ray muons and other non-collision backgrounds, events are required to have at least one primary vertex that is consistent with the beamspot position and that has at least five tracks associated to it. The efficiency for collision events to pass these vertex requirements, as measured in a sample of events passing all selections of this analysis, is well over 99%.

Jet Selection
For the inclusive jet measurements, jets are required to have p T > 20 GeV and to be within |y| < 4.4. They must also pass the specific fully-efficient trigger for each p T -and |y|-bin, as detailed in Sec. VI B. For the dijet measurements, events are selected if they have at least one jet with p T > 30 GeV and another jet with p T > 20 GeV, both within |y| < 4.4. Corrections are applied for inefficiencies in jet reconstruction, which are generally less than a few percent.
Jet quality criteria first established with early collision data are applied to reject jets reconstructed from calorimeter signals that do not originate from a protonproton collision, such as those due to noisy calorimeter cells [59]. For this analysis, various improvements to the jet quality selection have been made due to increased experience with a larger data set and evolving beam conditions, including the introduction of new criteria for the forward region.
The main sources of fake jets were found to be: noise bursts in the hadronic endcap calorimeter electronics; coherent noise from the electromagnetic calorimeter; cosmic rays; and beam-related backgrounds.
Quality selection criteria were developed for each of these categories by studying jet samples classified as real or fake energy depositions. This classification was performed by applying criteria on the magnitude and direction of the missing transverse momentum, E miss T . Following this, about a dozen events with | E miss T | > 500 GeV were found that pass the standard analysis selection. These events were visually scanned and were generally found to be collision events with mostly low p T jets and a muon escaping at low scattering angle.
The efficiency for identifying real jets was measured using a tag-and-probe method. A "probe jet" sample was selected by requiring the presence of a "tag jet" that is within |η| < 2.0, fulfills the jet quality criteria, and is back-to-back (∆φ > 2.6) and well-balanced with a probe jet (|p T1 −p T2 |/p avg T < 0.4, with p avg T = (p T1 +p T2 )/2 and where p T1,2 are the transverse momenta of the tag and probe jets). The jet quality criteria were then applied to the probe jet, measuring as a function of its |η| and p T the fraction of jets that are not rejected.
The efficiency to select a jet is shown in Fig. 6 for an example rapidity region, along with the systematic uncertainty on this efficiency.
The jet quality selection efficiency is greater than 96% for jets with p T = 20 GeV and quickly increases with jet p T . The efficiency is above 99% for jet p T > 60 GeV in all rapidity regions. The inclusive jet and dijet cross sections are corrected for these inefficiencies in regions where the efficiency is less than 99%. The systematic uncertainty on the efficiency is taken as a systematic uncertainty on the cross section. Efficiency for jet quality selection as a function of pT for anti-kt jets with R = 0.6 in the rapidity region |y| < 0.3. The black circles indicate the efficiency measured in-situ using a tag-probe method. The blue squares indicate the fit to the parameterisation (pT) = A − e −(Bp T −C) used in this analysis, where A, B, and C are fitted constants, and the shaded band indicates the systematic uncertainty on the efficiency obtained by varying the tag jet selection. The turnon is due to more stringent jet quality selection at low jet pT.

F. Background, Vertex Position, and Pile-Up
Background contributions from sources other than proton-proton collisions were evaluated using events from cosmic-ray runs, as well as unpaired proton bunches in the accelerator, in which no real collision candidates are expected. Based on the duration of the cosmic-ray runs and the fact that only one event satisfied the selection criteria, the non-collision background rates across the entire data period are considered to be negligible.
The primary vertices span the luminous region around the nominal beamspot. To determine the systematic uncertainty due to possibly incorrect modeling of the event vertex position, the jet p T spectrum was studied as a function of the |z| position of the primary vertex with the largest p 2 T of associated tracks. The fraction of events with |z| > 200 mm is 0.06%, and the difference in the p T spectrum compared to events with |z| < 100 mm is small. Consequently, the uncertainty from mis-modeling of the vertex position was taken to be negligible.
The p T of each jet is corrected for additional energy from soft pile-up interactions in the event (see Sec. VI C). An uncertainty associated to this pile-up offset correction is assigned that is dependent on the number of reconstructed primary vertices, as described in Sec. VIII A. The jet measurements are then compared to the Monte Carlo simulation without pile-up.

G. Luminosity
The integrated luminosity is calculated by measuring interaction rates using several ATLAS devices, where the absolute calibration is derived using van der Meer scans [62].
The uncertainty on the luminosity is 3.4% [63]. The calculation of the effective luminosity for each bin of the observable for inclusive jets follows the trigger scheme described in Sec. VI B. The integrated luminosity for each individual trigger is derived using separate prescale factors for each luminosity block (an interval of luminosity with homogeneous data-taking conditions, which is typically two minutes). For dijets, each bin receives contributions from several trigger combinations, for which the luminosity is calculated independently. The luminosity that would be obtained without correction for trigger prescale is (37.3 ± 1.2) pb −1 . Since the central jet trigger with the largest transverse momentum threshold was always unprescaled, this is the effective luminosity taken for jets with transverse momentum above about 220 GeV.

A. Technique used
Aside from the jet energy scale correction, all other corrections for detector inefficiencies and resolutions are performed using an iterative unfolding, based on a transfer matrix that relates the particle-level and reconstructionlevel observable, with the same binning as the final distribution. The unfolding is performed separately for each bin in rapidity since the migrations across rapidity bins are negligible compared to those across jet p T (dijet mass) bins. A similar procedure is applied for inclusive jets and dijets, with the following description applying specifically to the inclusive jet case.
The Monte Carlo simulation described in Sec. IV is used to derive the unfolding matrices. Particle-level and reconstructed jets are matched together based on geometrical criteria and used to derive a transfer matrix. This matrix contains the expected number of jets within each bin of particle-level and reconstructed jet p T . A folding matrix is constructed from the transfer matrix by normalising row-by-row so that the sum of the elements corresponding to a given particle-level jet p T is unity. Similarly, an unfolding matrix is constructed by normalising column-by-column so that the sum of the elements corresponding to a specific reconstructed jet p T is unity. Thus each element of the unfolding matrix reflects the probability for a reconstructed jet in a particular p T bin to originate from a specific particle-level p T bin, given the assumed input particle-level jet p T spectrum. The spectra of unmatched particle-level and reconstructed jets are also derived from the simulated sample. The ratio between the number of matched jets and the total number of jets provides the matching efficiency both for particlelevel jets, ptcl,i , and for reconstructed jets, reco,j .
The data are unfolded to particle level using a threestep procedure, with the final results being given by the equation: where i and j are the particle-level and reconstructed bin indices, respectively, and A ptcl,i reco,j is an unfolding matrix refined through iteration, as discussed below.
The first step is to multiply the reconstructed jet spectrum in data by the matching efficiency reco,j , such that it can be compared to the matched reconstructed spectrum from the Monte Carlo simulation. In the second step, the iterated unfolding matrix A ptcl,i reco,j is determined using the Iterative, Dynamically Stabilised (IDS) method [64]. This procedure improves the transfer matrix through a series of iterations, where the particlelevel distribution is reweighted to the shape of the corrected data spectrum, while leaving the folding matrix unchanged. The main difference with respect to previous iterative unfolding techniques [65] is that, when performing the corrections, regularisation is provided by the use of the significance of the data-MC differences in each bin. The third step is to divide the spectrum obtained after the iterative unfolding by the matching efficiency at particle level, thus correcting for the jet reconstruction inefficiency.
The statistical uncertainties on the spectrum are propagated through the unfolding by performing pseudoexperiments. An ensemble of pseudo-experiments is created in which each bin of the transfer matrix is varied according to its statistical uncertainty. A separate set of pseudo-experiments is performed where the data spectrum is varied while respecting correlations between jets produced in the same event. The unfolding is then applied to each pseudo-experiment, and the resulting ensembles are used to calculate the covariance matrix of the corrected spectrum.
As a cross-check, the results obtained from the iterative unfolding have been compared to those using a simpler bin-by-bin correction procedure, as well as the "singular value decomposition" (SVD) method implemented in TSVDUnfold [66, 67]. These methods use different regularisation procedures and rely to different degrees on the Monte Carlo simulation modelling of the shape of the spectrum. The unfolding techniques have been tested using a data-driven closure test [64]. In this test the particle-level spectrum in the Monte Carlo simulation is reweighted and convolved through the folding matrix such that a significantly improved agreement between the data and the reconstructed spectrum from the Monte Carlo simulation is attained. The reweighted, reconstructed spectrum in the Monte Carlo simulation is then unfolded using the same procedure as for the data. The comparison of the result with the reweighted particle- The jet shape ρ(r) measured using calorimeter energy clusters for anti-kt jets with R = 0.6 in the rapidity interval 2.8 ≤ |y| < 3.6, compared to PYTHIA with tune AMBT1 (used for unfolding), and for jets with transverse momenta in the range 30 < pT < 45 GeV. The statistical error bars are smaller than the size of the markers, while systematic errors are not shown.
level spectrum from the Monte Carlo simulation provides the estimation of the bias.
The bin-by-bin method gives results consistent with those obtained using the IDS technique, but requires the application of an explicit correction for the NLO k-factor to obtain good agreement. A somewhat larger bias is observed for the SVD method.

B. Cross-check with jet shapes
The use of Monte Carlo simulation to derive the transfer matrix in the unfolding procedure requires that the simulation models the jet properties well. The modelling of the energy flow around the jet core provides a useful test of this. The energy and momentum flow within a jet can be expressed in terms of the differential jet shape, defined for a jet with radius parameter R, as the fraction ρ(r) = 1 ∆r T is the transverse momentum within a radius R of the jet centre, and p r T is the transverse momentum contained within a ring of thickness ∆r = 0.1 at a radius r = (∆y) 2 + (∆φ) 2 from the jet centre.
Jet shape measurements using calorimeter energy clusters and tracks were performed with 3 pb −1 of data [68], and show good agreement with the PYTHIA and HER-WIG + JIMMY Monte Carlo simulations in the kinematic region 30 GeV < p T < 600 GeV and rapidity |y| < 2.8. Using the same technique, the uncorrected jet shapes in the forward rapidity region 2.8 ≤ |y| < 4.4 have been studied in the context of the present analysis. As an example, the results for the HEC-FCal transition region 2.8 ≤ |y| < 3.6, the most difficult detector region to model, are shown in Fig. 7. The maximum disagreement in shape between data and the Monte Carlo simulation is approximately 20%, demonstrating that the distribution of energy within the jets is reasonably well-modeled even in this worst case. Any bias from mis-modeling of the jet shape is included in the unfolding uncertainties described below, so this jet shape study serves only as a cross-check.

VIII. SYSTEMATIC UNCERTAINTIES AND CORRELATIONS
A. Uncertainty sources from jet reconstruction and calibration The uncertainty on the jet reconstruction efficiency for |y| < 2.1 (within the tracking acceptance) is evaluated using track jets, which are used to play the role of "truth jets". In this paper, truth jets are defined to be jets at the particle level, but excluding muons and neutrinos. The efficiency to reconstruct a calorimeter jet given a track jet nearby is studied in both data and the MC simulation. The data versus MC comparison of this efficiency is used to infer the degree to which the calorimeter jet reconstruction efficiency may be mis-modeled in the Monte Carlo simulation. The disagreement was found to be 2% for calorimeter jets with p T of 20 GeV and less than 1% for those with p T > 30 GeV. The disagreement for jets with |y| < 2.1 is taken as a systematic uncertainty for all jets in the rapidity range |y| < 4.4. This is expected to be a conservative estimate in the forward region where the jets have higher energy for a given p T .
The JES uncertainty was evaluated as described in Sec. VI D and in Ref. [59]. The jet energy and angular resolutions are estimated from the Monte Carlo simulation using truth jets that have each been matched to a reconstructed calorimeter jet. The jet energy resolution (JER) in the Monte Carlo simulation is compared to that obtained in data using two in-situ techniques, one based on dijet balance and the other using a bisector method [69]. In general the two resolutions agree within 14%, and the full difference is taken as a contribution to the uncertainty on the unfolding corrections, which propagates to a systematic uncertainty on the measured cross section as described in Sec. VIII B. The angular resolution is estimated from the angle between each calorimeter jet and its matched truth-level jet. The associated systematic uncertainty is assessed by varying the requirement that the jet is isolated.
The JES uncertainty due to pile-up is proportional to (N PV −1)/p T , where N PV is the number of reconstructed vertices. The total pile-up uncertainty for a given (p T , y)bin is calculated as the average of the uncertainties for each value of N PV weighted by the relative frequency of that number of reconstructed vertices in the bin.

B. Uncertainty propagation
The uncertainty of the measured cross section due to jet energy scale and jet energy and angular resolutions has been estimated using the Monte Carlo simulation by repeating the analysis after systematically varying these effects. The jet energy scale applied to the reconstructed jets in MC is varied separately for each JES uncertainty source both up and down by one standard deviation. The resulting p T spectra are unfolded using the nominal unfolding matrix, and the relative shifts with respect to the nominal unfolded spectra are taken as uncertainties on the cross section. The effects of the jet energy and angular resolutions are studied by smearing the reconstructed jets such that these resolutions are increased by one standard deviation of their respective uncertainties (see Sec. VIII A). For each such variation, a new transfer matrix is constructed, which is used to unfold the reconstructed jet spectrum of the nominal MC sample. The relative shift of this spectrum with respect to the nominal unfolded spectra is taken as the uncertainty on the cross section.
The impact of possible mis-modeling of the cross section shape in the Monte Carlo simulation is assessed by shape variations of the particle-level jet spectra introduced to produce reconstructed-level spectra in agreement with data as discussed in Sec. VII.
The total uncertainty on the unfolding corrections is defined as the sum in quadrature of the uncertainties on the jet energy resolution, jet angular resolution, and the simulated shape. It is approximately 4-5% at low and high p T (except for the lowest p T -bin at 20 GeV, where it reaches 20%), and is smaller at intermediate p T values. This uncertainty is dominated by the component from the jet energy resolution.

C. Summary of the magnitude of the systematic uncertainties
The largest systematic uncertainty for this measurement arises from the jet energy scale. Even with the higher precision achieved recently as described in Sec. VI D, the very steeply falling jet p T spectrum, especially for large rapidities, translates even relatively modest uncertainties on the transverse momentum into large changes for the measured cross section.
As described in Sec. VI G, the luminosity uncertainty is 3.4%. The detector unfolding uncertainties have been discussed in the previous subsection. Various other sources of systematic uncertainties were considered and were found to have a small impact on the results. The jet energy and angular resolutions, as well as the jet reconstruction efficiency, also contribute to the total uncertainty through the unfolding corrections.
The dominant systematic uncertainties for the measurement of the inclusive jet p T spectrum in representative p T and y regions for anti-k t jets with R = 0.6 are shown in Table I. Similarly, the largest systematic uncertainties for the dijet mass measurement are given for a few representative m 12 and y * regions in Table II. An example of the breakdown of the systematic uncertainties as a function of the jet transverse momentum for the various rapidity bins used in the inclusive jet measurement is shown in Fig. 8.

D. Correlations
The behaviour of various sources of systematic uncertainty in different parts of the detector has been studied in detail in order to understand their correlations across various p T , m 12 and rapidity bins. As shown in Tables III and IV, 22 independent sources of systematic uncertainty have been identified, including luminosity, jet energy scale and resolution, and theory effects such as the uncertainty of the modeling of the underlying event and the QCD showering. For example, the sources labeled "JES 7-13" in these tables correspond to the calorimeter response to hadrons, which dominates the JES uncertainty in the central region. After examining the rapidity dependence of all 22 sources, it was found that 87 independent nuisance parameters are necessary to describe the correlations over the whole phase space. The systematic effect on the cross section measurement associated with each nuisance parameter in its range of use is completely correlated in p T and y (dijet mass and y * ). These parameters represent correlations between the uncertainties of the various bins. Since many of the systematic effects are not symmetric, it is not possible to provide a covariance matrix containing the full information. For symmetric uncertainties corresponding to independent sources, the total covariance matrix is given by: where λ is an index running over the nuisance parameters, and Γ λi is the one-standard-deviation amplitude of the systematic effect due to source λ in bin i. The full list of relative uncertainties, γ λ , where each uncertainty may be asymmetric, is given for all sources λ and bins of this analysis in Tables V-XVIII and XIX-XXXVI. Fig. 8 shows the magnitude and approximate bin-to-bin correlations of the total systematic uncertainty of the inclusive jet cross section measurement. The correlation matrix is here converted from the covariance matrix, which is obtained using Eq. 6, after symmetrising the uncertainties: The inclusive jet and dijet data should not be used simultaneously for PDF fits due to significant correlations between the two measurements.

A. Inclusive Jet Cross Sections
The inclusive jet double-differential cross section is shown in Figs. 9 and 10 and Tables V-XVIII in Appendix B for jets reconstructed with the anti-k t algorithm with R = 0.4 and R = 0.6. The measurement extends from jet transverse momentum of 20 GeV to almost 1.5 TeV, spanning two orders of magnitude in p T and ten orders of magnitude in the value of the cross section. The measured cross sections have been corrected for all detector effects using the unfolding procedure described in Sec. VII. The results are compared to NLOJET++ predictions (using the CT10 PDF set) corrected for nonperturbative effects, where the theoretical uncertainties from scale variations, parton distribution functions, and non-perturbative corrections have been accounted for.
In Figs. 11-13, the inclusive jet results are presented in terms of the ratio with respect to the NLOJET++ predictions using the CT10 PDF set. Fig. 11 compares the current results to the previous measurements published by ATLAS [21], for jets reconstructed with the anti-k t algorithm with parameter R = 0.6. This figure is limited to the central region, but similar conclusions can be drawn in all rapidity bins. In particular the two measurements are in good agreement, although the new results cover a much larger kinematic range with much reduced statistical and systematic uncertainties. Fig. 12 shows the ratio of the measured cross sections to the NLOJET++ theoretical predictions for various PDF sets. Predictions obtained using CT10, MTSW 2008, NNPDF 2.1, and HERAPDF 1.5, including uncertainty bands, are compared to the measured cross sections, where data and theoretical predictions are normalised to the prediction from the CT10 PDF set. The data show a marginally smaller cross section than the predictions from each of the PDF sets. This trend is more Description of bin-to-bin uncertainty correlation for the inclusive jet measurement. Each number corresponds to a nuisance parameter for which the corresponding uncertainty is fully correlated versus pT. Bins with the same nuisance parameter are treated as fully correlated, while bins with different nuisance parameters are uncorrelated. The sources indicated by the letter "u" are uncorrelated both between pT-and |y|-bins. The one-standard-deviation amplitude of the systematic effect associated with each nuisance parameter is detailed in Tables V-XVIII in Appendix B. The JES uncertainties for jets with |y| ≥ 0.8 are determined relative to the JES of jets with |y| < 0.8. As a consequence, several of the uncertainties that are determined using jets with |y| < 0.8 are also propagated to the more forward rapidities (such as the E/p uncertainties).  Noise threshold  1  1  2  3  4  4  5  6  6  JES 2: Theory UE  7  7  8  9  10  10  11  12  12  JES 3: Theory Showering  13  13  14  15  16  16  17  18  18  JES 4: Non-closure  19  19  20  21  22  22  23  24  24  JES 5: Dead material  25  25  26  27  28  28  29  30   pronounced for the measurements corresponding to the anti-k t algorithm with parameter R = 0.4, compared to R = 0.6. The description becomes worse for large jet transverse momenta and rapidities, where the MSTW 2008 PDF set follows the measured trend better. However, the differences between the measured cross section and the prediction of each PDF set are of the same order as the total systematic uncertainty on the measurement, including both experimental and theoretical uncertainty sources. A χ 2 test of the compatibility between data and the PDF curves, accounting for correlations between bins, provides reasonable probabilities for all sets, with non-significant differences between them. 5 The comparison of the data with the POWHEG prediction, using the CT10 NLO PDF set, is shown for anti-k t jets with R = 0.4 and R = 0.6 in different rapidity regions in Fig. 13. The data are compared with four theory curves, all of which are normalised to the same common denominator of the NLOJET++ prediction corrected for non-perturbative effects: POWHEG showered with PYTHIA with the default AUET2B tune; the same with the Perugia 2011 tune; POWHEG showered with HERWIG; and POWHEG run in "pure NLO" mode (fixed-order calculation), without matching to parton shower, after application of soft corrections calculated using PYTHIA and the AUET2B tune. Scale uncertainties are not shown for the POWHEG curves, but they have been found to be similar to those obtained with NLOJET++.
Good agreement at the level of a few percent is observed between NLO fixed-order calculations based on NLOJET++ and POWHEG, as described in Sec. V A 1. However, significant differences reaching O(30%) are observed if POWHEG is interfaced to different showering and soft physics models, particularly at low p T and forward rapidity, but also at high p T . These differences exceed the uncertainties on the non-perturbative corrections, which are not larger than 10% for the inclusive jet measurements with R = 0.4, thus indicating a significant impact of the parton shower. The Perugia 2011 tune tends to produce a consistently larger cross section than the standard AUET2B tune over the full rapidity range. The technique of correcting fixed-order calcula-  . Inclusive jet double-differential cross section as a function of jet pT in different regions of |y| for jets identified using the anti-kt algorithm with R = 0.4. For convenience, the cross sections are multiplied by the factors indicated in the legend. The data are compared to NLO pQCD calculations using NLOJET++ to which non-perturbative corrections have been applied. The error bars, which are usually smaller than the symbols, indicate the statistical uncertainty on the measurement. The dark-shaded band indicates the quadratic sum of the experimental systematic uncertainties, dominated by the jet energy scale uncertainty. There is an additional overall uncertainty of 3.4% due to the luminosity measurement that is not shown. The theory uncertainty, shown as the light, hatched band, is the quadratic sum of uncertainties from the choice of the renormalisation and factorisation scales, parton distribution functions, αs(MZ ), and the modeling of non-perturbative effects, as described in the text.
tions for non-perturbative effects remains the convention to define the baseline theory prediction until NLO parton shower generators become sufficiently mature to describe data well. The corrected NLO result predicts a consistently larger cross section than that seen in the data. Good agreement in normalisation is found between the data and the prediction from POWHEG showered with the default tune of PYTHIA. These results are confirmed by a χ 2 test of the compatibility of the POWHEG results with the data, where the curve obtained using the HER-WIG shower results in a much worse χ 2 after all error correlations have been accounted for.
[GeV]  Inclusive jet double-differential cross section as a function of jet pT in different regions of |y| for jets identified using the anti-kt algorithm with R = 0.6. For convenience, the cross sections are multiplied by the factors indicated in the legend. The data are compared to NLO pQCD calculations using NLOJET++ to which non-perturbative corrections have been applied. The theoretical and experimental uncertainties indicated are calculated as described in Fig. 9.
Ratio of inclusive jet cross section to the theoretical prediction obtained using NLOJET++ with the CT10 PDF set. The ratio is shown as a function of jet pT in the rapidity region |y| < 0.3, for jets identified using the anti-kt algorithm with R = 0.6. The current result is compared to that published in Ref. [21].       13. Ratios of inclusive jet double-differential cross section to the theoretical prediction obtained using NLOJET++ with the CT10 PDF set. The ratios are shown as a function of jet pT in different regions of |y| for jets identified using the anti-kt algorithm with R = 0.4 (upper plots) and R = 0.6 (lower plots). The ratios of POWHEG predictions showered using either PYTHIA or HERWIG to the NLOJET++ predictions corrected for non-perturbative effects are shown and can be compared to the corresponding ratios for data. Only the statistical uncertainty on the POWHEG predictions is shown. The total systematic uncertainties on the theory and the measurement are indicated. The NLOJET++ prediction and the POWHEG ME calculations use the CT10 PDF set. Statistically insignificant data points at large pT are omitted in the ratio.

B. Dijet Cross Sections
The dijet double-differential cross section has been measured as a function of the dijet invariant mass for various bins of the variable y * , which is the rapidity in the two-parton centre-of-mass frame. The quantity y * is calculated as half the absolute value of the rapidity difference of the two leading jets, ranging from 0 to 4.4. The results are shown in Figs. 14 and 15 and Tables XIX-XXXVI in Appendix C for anti-k t jets with R = 0.4 and R = 0.6. The cross section measurements extend from dijet masses of 70 GeV to almost 5 TeV, covering two orders of magnitude in invariant mass and nine orders of magnitude in the cross section. The dijet measurements are fully corrected for detector effects and are compared to NLOJET++ predictions calculated using the scale defined in Eq. 3 (see Sec. V A 1) and the CT10 PDF set, with non-perturbative corrections applied to the theory prediction. The theoretical uncertainties have been assessed as described for the inclusive jet measurements in Sec. IX A.
The dijet data are also compared with NLOJET++ predictions obtained using the MSTW 2008, NNPDF 2.1, and HERAPDF 1.5 PDF sets. Figs. 16 and 17 show the dijet mass spectra for anti-k t jets with R = 0.4 and R = 0.6 respectively, where both the data and the predictions from the above-mentioned PDF sets have been normalised to the CT10 prediction. The data for R = 0.6 exhibit a slight falling slope with respect to the CT10 prediction and appear to be described better by other PDF sets, a similar behaviour to that observed in the inclusive jet data. However, in all cases, the differences between the data and each PDF set lie well within the systematic and theory uncertainties, indicating a reasonable agreement with the dijet data, particularly in the kinematic region at low y * .
The data are also compared with POWHEG predictions produced using the CT10 PDF set and showered with different tunes of the PYTHIA or HERWIG generator. These comparisons are shown for R = 0.4 and R = 0.6 respectively in Figs. 18 and 19, where the data and all theory predictions have been normalised to the NLOJET++ prediction with CT10. The NLOJET++ prediction has been corrected for non-perturbative effects calculated using the PYTHIA MC with the AUET2B tune. The POWHEG predictions shown are interfaced to the PYTHIA parton shower with the AUET2B or Pe-rugia2011 tune, and to the HERWIG parton shower using the AUET2 tune. The data are also compared to the POWHEG fixed-order NLO prediction (corrected for non-perturbative effects), where the POWHEG prediction has been calculated using a scale choice of The data are in best agreement with the POWHEG prediction showered with PYTHIA using the AUET2B tune. The other POWHEG showered predictions exhibit discrepancies at low dijet mass in all y * slices, where they predict larger cross sections than are observed in the data.

X. CONCLUSIONS
Cross section measurements have been presented for inclusive jets and dijets reconstructed with the anti-k t algorithm using two values of the clustering parameter (R = 0.4 and R = 0.6). Inclusive jet production has been measured as a function of jet transverse momentum, in bins of jet rapidity. Dijet production has been measured as a function of the invariant mass of the two leading jets, in bins of half their rapidity difference. These results are based on the data sample collected with the ATLAS detector during 2010, which corresponds to (37.3 ± 1.2) pb −1 of integrated luminosity.
Two different sizes of the jet clustering parameter have been used in order to probe the relative effects of the parton shower, hadronisation, and the underlying event. The measurements have been corrected for all detector effects to the particle level so that they can be compared to any theoretical calculation. In this paper, they have been compared to fixed-order NLO pQCD calculations corrected for non-perturbative effects, as well as to parton shower Monte Carlo simulations with NLO matrix elements. The latter predictions have only recently become available for inclusive jet and dijet production.
The current results reflect a number of significant experimental accomplishments: • The cross section measurements extend to 1.5 TeV in jet transverse momentum and 5 TeV in dijet invariant mass, the highest ever measured. These results probe NLO pQCD in a large, new kinematic regime.
• Using data taken with minimum bias and forward jet triggers, these measurements extend to both the low-p T region (down to jet transverse momentum of 20 GeV and dijet invariant mass of 70 GeV) and to the forward region (out to rapidities of |y| = 4.4). The forward region, in particular, has never been explored before with such precision at a hadronhadron collider.
• High-precision measurements of the data collected during LHC beam position scans have determined the uncertainty on the collected luminosity to 3.4%.
• Detailed understanding of the detector performance has precisely determined systematic uncertainties, in particular those arising from the jet energy scale. In the central region (|η| < 0.8) the JES uncertainty is lower than 4.6% for all jets with p T > 20 GeV, while for jet transverse momenta between 60 and 800 GeV the JES uncertainty is below 2.5%.
• The correlations of the cross section measurement across various p T , m 12 , and rapidity bins have been studied for 22 independent sources of systematic uncertainty. These have been provided in the form of 87 nuisance parameters, each of which is fully correlated in p T and y (dijet mass and y * ), for use in PDF fits.
The experimental uncertainties achieved are similar in size to the theoretical uncertainties in some regions of phase space, thereby providing some sensitivity to different theoretical predictions. The measurements are compared to fixed-order NLO pQCD calculations, as well as to new calculations in which NLO pQCD matrix elements are matched to leading-logarithmic parton showers. Overall, both sets of calculations agree with the data over many orders of magnitude, although the cross sections predicted by the theory tend to be larger than the measured values at large jet transverse momentum and dijet invariant mass. The matched NLO parton shower calculations predict significant effects of the parton shower in some regions of phase space, in some cases improving and in others degrading the agreement with data with respect to the fixed-order calculations.
These measurements probe and may constrain the largely unexplored area of parton distribution functions at large x and high momentum transfer. The results reported here constitute a comprehensive test of QCD across a large kinematic regime.
[TeV]  16 10 × y* < 4.4 ( ≤ 4.0 ) 14 10 × y* < 4.0 ( ≤ 3.5 ) 12 10 × y* < 3.5 ( ≤ 3.0 ) 10 10 × y* < 3.0 ( ≤ 2.5 ) 8 10 14. Dijet double-differential cross section as a function of dijet mass, binned in half the rapidity separation between the two leading jets, y * = |y1 − y2|/2. The results are shown for jets identified using the anti-kt algorithm with R = 0.4. For convenience, the cross sections are multiplied by the factors indicated in the legend. The data are compared to NLO pQCD calculations using NLOJET++ to which non-perturbative corrections have been applied. The error bars, which are usually smaller than the symbols, indicate the statistical uncertainty on the measurement. The dark-shaded band indicates the quadratic sum of the experimental systematic uncertainties, dominated by the jet energy scale uncertainty. There is an additional overall uncertainty of 3.4% due to the luminosity measurement that is not shown. The theory uncertainty, shown as the light, hatched band, is the quadratic sum of uncertainties from the choice of the renormalisation and factorisation scales, parton distribution functions, αs(MZ ), and the modeling of non-perturbative effects, as described in the text.
[TeV]  16 10 × y* < 4.4 ( ≤ 4.0 ) 14 10 × y* < 4.0 ( ≤ 3.5 ) 12 10 × y* < 3.5 ( ≤ 3.0 ) 10 10 × y* < 3.0 ( ≤ 2.5 ) 8 10 15. Dijet double-differential cross section as a function of dijet mass, binned in half the rapidity separation between the two leading jets, y * = |y1 − y2|/2. The results are shown for jets identified using the anti-kt algorithm with R = 0.6. For convenience, the cross sections are multiplied by the factors indicated in the legend. The data are compared to NLO pQCD calculations using NLOJET++ to which non-perturbative corrections have been applied. The theoretical and experimental uncertainties indicated are calculated as described in Fig. 14 18. Ratios of dijet double-differential cross section to the theoretical prediction obtained using NLOJET++ with the CT10 PDF set. The ratios are shown as a function of dijet mass, binned in half the rapidity separation between the two leading jets, y * = |y1 −y2|/2. The results are shown for jets identified using the anti-kt algorithm with R = 0.4. The ratios of POWHEG predictions showered using either PYTHIA or HERWIG to the NLOJET++ predictions corrected for non-perturbative effects are shown and can be compared to the corresponding ratios for data. Only the statistical uncertainty on the POWHEG predictions is shown. The total systematic uncertainties on the theory and the measurement are indicated. The NLOJET++ prediction and the POWHEG ME calculations use the CT10 PDF set.   19. Ratios of dijet double-differential cross section to the theoretical prediction obtained using NLOJET++ with the CT10 PDF set. The ratios are shown as a function of dijet mass, binned in half the rapidity separation between the two leading jets, y * = |y1 −y2|/2. The results are shown for jets identified using the anti-kt algorithm with R = 0.6. The ratios of POWHEG predictions showered using either PYTHIA or HERWIG to the NLOJET++ predictions corrected for non-perturbative effects are shown and can be compared to the corresponding ratios for data. Only the statistical uncertainty on the POWHEG predictions is shown. The total systematic uncertainties on the theory and the measurement are indicated. The NLOJET++ prediction and the POWHEG ME calculations use the CT10 PDF set. FIG. 20. Non-perturbative correction factors for inclusive jets identified using the anti-kt algorithm with distance parameters R = 0.4 and R = 0.6 in various rapidity regions, derived using various Monte Carlo generators. The correction derived using PYTHIA 6.425 with the AUET2B CTEQ6L1 tune is used for the fixed-order NLO calculations presented in this analysis.