Jet energy resolution in proton-proton 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{\mathrm{s}}=7\mbox{ TeV}$\end{document} recorded in 2010 with the ATLAS detector

The measurement of the jet energy resolution is presented using data recorded with the ATLAS detector in proton-proton 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\mbox{ TeV}$\end{document}. The sample corresponds to an integrated luminosity of 35 pb−1. Jets are reconstructed from energy deposits measured by the calorimeters and calibrated using different jet calibration schemes. The jet energy resolution is measured with two different in situ methods which are found to be in agreement within uncertainties. The total uncertainties on these measurements range from 20 % to 10 % for jets within |y|<2.8 and with transverse momenta increasing from 30 GeV to 500 GeV. Overall, the Monte Carlo simulation of the jet energy resolution agrees with the data within 10 %.


Introduction
Precise knowledge of the jet energy resolution is of key importance for the measurement of the cross-sections of inclusive jets, dijets, multijets or vector bosons accompanied by jets [1][2][3][4], top-quark cross-sections and mass measurements [5], and searches involving resonances decaying to jets [6,7]. The jet energy resolution also has a direct impact on the determination of the missing transverse energy, which plays an important role in many searches for new physics with jets in the final state [8,9]. This article presents the determination with the ATLAS detector [10, 11] of the jet energy resolution in proton-proton collisions at a centre-ofmass energy of √ s = 7 TeV. The data sample was collected during 2010 and corresponds to 35 pb −1 of integrated luminosity delivered by the Large Hadron Collider (LHC) [12] at CERN. The jet energy resolution is determined by exploiting the transverse momentum balance in events containing jets with large transverse momenta (p T ). This article is structured as follows: Sect. 2 describes the ATLAS detector. Sections 3,4 and 5 respectively introduce the Monte Carlo simulation, the event and jet selection criteria, and the jet calibration methods. The two techniques to estimate the jet energy resolution from calorimeter observables, the dijet balance method [13] and the bisector method [14], are discussed respectively in Sects. 6 and 7. These methods rely on somewhat different assumptions, which can be validated in data and are sensitive to different sources of systematic uncertainty. As such, the use of these two independent in situ measurements of the jet energy resolution is important to validate the Monte Carlo simulation. Section 8 presents the results obtained for data and simulation for the default jet energy calibration scheme implemented in ATLAS. Section 9 compares the resolutions obtained by applying the two in situ methods to the Monte Carlo simulation and the resolutions determined by comparing the jet energy at calorimeter and particle level. This comparison will be referred to as a closure test. Sources of systematic uncertainty on the jet energy resolution estimated using the available Monte Carlo simulations and collision data are discussed in Sect. 10. The results for other jet energy calibration schemes are discussed in Sects. 11 and 12, and the conclusions can be found in Sect. 13.

The ATLAS detector
The ATLAS detector is a multi-purpose detector designed to observe particles produced in high energy proton-proton collisions. A detailed description can be found in Refs. [10,11]. The Inner (tracking) Detector has complete azimuthal coverage and spans the pseudorapidity region |η| < 2.5. 1 The Inner Detector consists of layers of silicon pixel, silicon microstrip and transition radiation tracking detectors. These sub-detectors are surrounded by a superconducting solenoid that produces a uniform 2 T axial magnetic field.
The calorimeter system is composed of several subdetectors. A high-granularity liquid-argon (LAr) electromagnetic sampling calorimeter covers the |η| < 3.2 range, and it is split into a barrel (|η| < 1.475) and two end-caps (1.375 < |η| < 3.2). Lead absorber plates are used over its 1 The ATLAS reference system is a Cartesian right-handed coordinate system, with the nominal collision point at the origin. The anticlockwise beam direction defines the positive z-axis, with the x-axis pointing to the centre of the LHC ring. The angle φ defines the direction in the plane transverse to the beam (x, y). The pseudorapidity is given by η = − ln tan θ 2 , where the polar angle θ is taken with respect to the positive z direction. The rapidity is defined as y = 0.5 × ln[(E + p z )/(E − p z )], where E denotes the energy and p z is the component of the momentum along the z-axis. full coverage. The hadronic calorimetry in the barrel is provided by a sampling calorimeter using steel as the absorber material and scintillating tiles as active material in the range |η| < 1.7. This tile hadronic calorimeter (Tilecal) is separated into a large barrel (|η| < 0.8) and two smaller extended barrel cylinders, one on either side of the central barrel. In the end-caps, copper/LAr technology is used for the hadronic end-cap calorimeters (HEC), covering the range 1.5 < |η| < 3.2. The copper-tungsten/LAr forward calorimeters (FCal) provide both electromagnetic and hadronic energy measurements, extending the coverage to |η| = 4.9.
The trigger system consists of a hardware-based Level 1 (L1) and a two-tier, software-based High Level Trigger (HLT). The L1 jet trigger uses a sliding window algorithm with coarse-granularity calorimeter towers. This is then refined using jets reconstructed from calorimeter cells in the HLT.

Event generators
Data are compared to Monte Carlo (MC) simulations of jets with large transverse momentum produced via strong interactions described by Quantum Chromodynamics (QCD) in proton-proton collisions at a centre-of-mass energy of √ s = 7 TeV. The jet energy resolution is derived for several simulation models in order to study its dependence on the event generator, on the parton showering and hadronisation models, and on tunes of other soft model parameters, such as those of the underlying event. The event generators used for this analysis are described below.
1. PYTHIA 6.4 MC10 tune: The event generator PYTHIA [15] simulates non-diffractive proton-proton collisions using a 2 → 2 matrix element at the leading order (LO) of the strong coupling constant to model the hard subprocess, and uses p T -ordered parton showers to model additional radiation in the leading-logarithm approximation [

Simulated pile-up samples
The nominal MC simulation does not include additional proton-proton interactions (pile-up). In order to study its effect on the jet energy resolution, two additional MC samples are used. The first one simulates additional protonproton interactions in the same bunch crossing (in-time pileup) while the second sample in addition simulates effects on calorimeter cell energies from close-by bunches (out-oftime pile-up). The average number of interactions per event is 1.7 (1.9) for the in-time (in-time plus out-of-time) pile-up samples, which is a good representation of the 2010 data.

Event and jet selection
The status of each sub-detector and trigger, as well as reconstructed physics objects in ATLAS is continuously assessed by inspection of a standard set of distributions, and dataquality flags are recorded in a database for each luminosity block (of about two minutes of data-taking). This analysis selects events satisfying data-quality criteria for the Inner Detector and the calorimeters, and for track, jet, and missing transverse energy reconstruction [39]. For each event, the reconstructed primary vertex position is required to be consistent with the beamspot, both transversely and longitudinally, and to be reconstructed from at least five tracks with transverse momentum p track T > 150 MeV associated with it. The primary vertex is defined as the one with the highest associated sum of squared track transverse momenta Σ(p track T ) 2 , where the sum runs over all tracks used in the vertex fit. Events are selected by requiring a specific OR combination of inclusive single-jet and dijet calorimeter-based triggers [40,41]. The combinations are chosen such that the trigger efficiency for each p T bin is greater than 99 %. For the lowest p T bin (30-40 GeV), this requirement is relaxed, allowing the lowest-threshold calorimeter inclusive single-jet trigger to be used with an efficiency above 95 %.
Jets are reconstructed with the anti-k t jet algorithm [42] using the FastJet software [43] with radius parameters R = 0.4 or R = 0.6, a four-momentum recombination scheme, and three-dimensional calorimeter topological clusters [44] as inputs. Topological clusters are built from calorimeter cells with a signal at least four times higher than the rootmean-square (RMS) of the noise distribution (seed cells). Cells neighbouring the seed which have a signal to RMSnoise ratio ≥ 2 are then iteratively added. Finally, all nearest neighbour cells are added to the cluster without any threshold.
Jets from non-collision backgrounds (e.g. beam-gas events) and instrumental noise are removed using the selection criteria outlined in Ref. [39].
Jets are categorized according to their reconstructed rapidity in four different regions to account for the differently instrumented parts of the calorimeter: -Central region (|y| < 0.8).
Events are selected only if the transverse momenta of the two leading jets are above a jet reconstruction threshold of 7 GeV at the electromagnetic scale (see Sect. 5) and within |y| ≤ 2.8, at least one of them being in the central region. The analysis is restricted to |y| ≤ 2.8 because of the limited number of jets at higher rapidities. Monte Carlo simulated "particle jets" are defined as those built using the same jet algorithm as described above, but using instead as inputs the stable particles from the event generator (with a lifetime longer than 10 ps), excluding muons and neutrinos.

Jet energy calibration
Calorimeter jets are reconstructed from calorimeter energy deposits measured at the electromagnetic scale (EM-scale), the baseline signal scale for the energy deposited by electromagnetic showers in the calorimeter. Their transverse momentum is referred to as p EM-scale T . For hadrons this leads to a jet energy measurement that is typically 15-55 % lower than the true energy, due mainly to the non-compensating nature of the ATLAS calorimeter [45]. Fluctuations of the hadronic shower, in particular of its electromagnetic content, as well as energy losses in the dead material lead to a degraded resolution and jet energy response compared to particles interacting only electromagnetically. The jet response is defined as the ratio of calorimeter jet p T and particle jet p T (see Sect. 4), reconstructed with the same algorithm, and matched in η − φ space (see Sect. 9). Several complementary jet calibration schemes with different levels of complexity and different sensitivity to systematic effects have been developed to understand the jet energy measurements. The jet calibration is performed by applying corrections derived from Monte Carlo simulations to restore the jet response to unity. This is referred to as determining the jet energy scale (JES).
The analysis presented in this article aims to determine the jet energy resolution for jets reconstructed using various JES strategies. A simple calibration, referred to as the EM + JES calibration scheme, has been chosen for the first physics analysis of the 2010 data [39]. It allows a direct evaluation of the systematic uncertainties from single-hadron response measurements and is therefore suitable for first physics analyses. More sophisticated calibration techniques to improve the jet resolution and reduce partonic flavour response differences have also been developed. They are the Local Cluster Weighting (LCW), the Global Cell Weighting (GCW) and the Global Sequential (GS) methods [39]. In addition to these calorimeter calibration schemes, a Track-Based Jet Correction (TBJC) has been derived to adjust the response and reduce fluctuations on a jet-by-jet basis without changing the average jet energy scale. These calibration techniques are briefly described below.

The EM + JES calibration
For the analysis of the first proton-proton collisions, a simple Monte Carlo simulation-based correction is applied as the default to restore the hadronic energy scale on average. The EM + JES calibration scheme applies corrections as a function of the jet transverse momentum and pseudorapidity to jets reconstructed at the electromagnetic scale. The main advantage of this approach is that it allows the most direct evaluation of the systematic uncertainties. The uncertainty on the absolute jet energy scale was determined to be less than ±2.5 % in the central calorimeter region (|y| < 0.8) and ±14 % in the most forward region (3.2 ≤ |y| < 4.5) for jets with p T > 30 GeV [39]. These uncertainties were evaluated using test-beam results, single hadron response in situ measurements, comparison with jets built from tracks, p T balance in dijet and γ + jet events, estimations of pile-up energy deposits, and detailed Monte Carlo comparisons.

The Local Cluster Weighting (LCW) calibration
The LCW calibration scheme uses properties of clusters to calibrate them individually prior to jet finding and reconstruction. The calibration weights are determined from Monte Carlo simulations of charged and neutral pions according to the cluster topology measured in the calorimeter. The cluster properties used are the energy density in the cells forming them, the fraction of their energy deposited in the different calorimeter layers, the cluster isolation and its depth in the calorimeter. Corrections are applied to the cluster energy to account for the energy deposited in the calorimeter but outside of clusters and energy deposited in material before and in between the calorimeters. Jets are formed from calibrated clusters. A final jet-level energy correction based on the same procedure as for the EM + JES case is applied to attain unity response, but with corrections that are numerically smaller. The resulting jet energy calibration is denoted as LCW + JES.

The Global Cell Weighting (GCW) calibration
The GCW calibration scheme attempts to compensate for the different calorimeter response to hadronic and electromagnetic energy deposits at cell level. The hadronic signal is characterized by low cell energy densities and, thus, a positive weight is applied. The weights, which depend on the cell energy density and the calorimeter layer only, are determined by minimizing the jet resolution evaluated by comparing reconstructed and particle jets in Monte Carlo simulation. They correct for several effects at once (calorimeter non-compensation, dead material, etc.). A jet-level correction is applied to jets reconstructed from weighted cells to account for global effects. The resulting jet energy calibration is denoted as GCW + JES.

The Global Sequential (GS) calibration
The GS calibration scheme uses the longitudinal and transverse structure of the jet calorimeter shower to compensate for fluctuations in the jet energy measurement. In this scheme the jet energy response is first calibrated with the EM + JES calibration. Subsequently, the jet properties are used to exploit the topology of the energy deposits in the calorimeter to characterize fluctuations in the hadronic shower development. These corrections are applied such that the mean jet energy is left unchanged, and each correction is applied sequentially. This calibration is designed to improve the jet energy resolution without changing the average jet energy scale.

Track-based correction to the jet calibration
Regardless of the inputs, algorithms and calibration methods chosen for calorimeter jets, more information on the jet topology can be obtained from reconstructed tracks associated to the jet. Calibrated jets have an average energy response close to unity. However, the energy of an individual jet can be over-or underestimated depending on several factors, for example: the ratio of the electromagnetic and hadronic components of the jet; the fraction of energy lost in dead material, in either the inner detector, the solenoid, the cryostat before the LAr, or the cryostat between the LAr and the TileCal. The reconstructed tracks associated to the jet are sensitive to some of these effects and therefore can be used to correct the calibration on a jet-by-jet basis.
In the method referred to as Track-Based Jet Correction (TBJC) [45], the response is adjusted depending on the number of tracks associated with the jet. The jet energy response is observed to decrease with increasing track multiplicity of the jets, mainly because the ratio of the electromagnetic to the hadronic component decreases on average as the number of tracks increases. In effect, a low charged-track multiplicity typically indicates a predominance of neutral hadrons, in particular π 0 s which yield electromagnetic deposits in the calorimeter with R 1. A large number of charged particles, on the contrary, signals a more dominant hadronic component, with a lower response due to the non-compensating nature of the calorimeter (h/e < 1). The TBJC method is designed to be applied as an option in addition to any JES calibration scheme, since it does not change the average response, to reduce the jet-to-jet energy fluctuations and improve the resolution.

In situ jet resolution measurement using the dijet balance method
Two methods are used in dijet events to measure in situ the fractional jet p T resolution, σ (p T )/p T , which at fixed rapidity is equivalent to the fractional jet energy resolution, σ (E)/E. The first method, presented in this section, relies on the approximate scalar balance between the transverse momenta of the two leading jets and measures the sensitivity of this balance to the presence of extra jets directly from data. The second one, presented in the next section, uses the projection of the vector sum of the leading jets' transverse momenta on the coordinate system bisector of the azimuthal angle between the transverse momentum vectors of the two jets. It takes advantage of the very different sensitivities of each of these projections to the underlying physics of the dijet system and to the jet energy resolution.

Measurement of resolution from asymmetry
The dijet balance method for the determination of the jet p T resolution is based on momentum conservation in the transverse plane. The asymmetry between the transverse momenta of the two leading jets A(p T,1 , p T,2 ) is defined as where p T,1 and p T,2 refer to the randomly ordered transverse momenta of the two leading jets. The width σ (A) of a Gauss distribution fitted to A(p T,1 , p T,2 ) is used to characterize the asymmetry distribution and determine the jet p T resolutions. For events with exactly two particle jets that satisfy the hypothesis of momentum balance in the transverse plane, and requiring both jets to be in the same rapidity region, the relation between σ (A) and the fractional jet resolution is given by where σ (p T,1 ) = σ (p T,2 ) = σ (p T ), since both jets are in the same y region. If one of the two leading jets (j ) is in the rapidity bin being probed and the other one (i) in a reference y region where the resolution may be different, the fractional jet p T resolution is given by where A (i,j ) is measured in a topology with the two jets in different rapidity regions and where (i) ≡ (i, i) denotes both jets in the same y region.
The back-to-back requirement is approximated by an azimuthal angle cut between the leading jets, φ(j 1 , j 2 ) ≥ 2.8, and a veto on the third jet momentum, p EM-scale T,3 < 10 GeV, with no rapidity restriction. The resulting asymmetry distribution is shown in Fig. 1 for ap T ≡ (p T,1 + p T,2 )/2 bin of 60 GeV ≤p T < 80 GeV, in the central region (|y| < 0.8). Reasonable agreement in the bulk is observed between data and Monte Carlo simulation.

Soft radiation correction
Although requirements on the azimuthal angle between the leading jets and on the third jet transverse momentum are designed to enrich the purity of the back-to-back jet sample, it is important to account for the presence of additional soft particle jets not detected in the calorimeter.
In order to estimate the value of the asymmetry for a pure particle dijet event, → 0, in order to estimate the expected resolution for an ideal dijet topology The dependence of the jet p T resolution on the presence of a third jet is illustrated in Fig. 2. The linear fits and their extrapolations for ap T bin of 60 ≤p T < 80 GeV are shown. Note that the resolutions become systematically broader as the p EM-scale T,3 cut increases. This is a clear indication that the jet resolution determined from two-jet topologies depends on the presence of additional radiation and on the underlying event.
A soft radiation (SR) correction factor, K soft (p T ), is obtained from the ratio of the values of the linear fit at 0 GeV This multiplicative correction is applied to the resolutions extracted from the dijet asymmetry for p EM-scale T,3 < 10 GeV events. The correction varies from 25 % for events withp T of 50 GeV down to 5 % forp T of 400 GeV. In order to limit the statistical fluctuations, K soft (p T ) is fit with a parameterization of the form K soft (p T ) = a + b/(logp T ) 2 , which was found to describe the distribution well, within uncertainties. The differences in the resolution due to other parameterizations of K were studied and treated as a systematic uncertainty, resulting in a relative uncertainty of about 6 % (see Sect. 10).

Particle balance correction
The p T difference between the two calorimeter jets is not solely due to resolution effects, but also to the balance between the respective particle jets, The measured difference (left side) is decomposed into resolution fluctuations (the first two terms on the right side) plus a particle-level balance (PB) term that originates from out-of-jet showering in the particle jets. In order to correct for this contribution, the particle-level balance is estimated Fig. 3 Fractional jet resolution obtained in simulation using the dijet balance method, shown as a function ofp T , both before (circles) and after the particle-balance (PB) correction (triangles). Also shown is the dijet PB correction itself (squares) and, in the lower panel, its relative size with respect to the fractional jet resolution. The curves correspond to fits with the functional form in Eq. (9). The errors shown are only statistical using the same technique (asymmetry plus soft radiation correction) as for calorimeter jets. The contribution of the dijet PB after the SR correction is subtracted in quadrature from the in situ resolution for both data and Monte Carlo simulation. The result of this procedure is shown for simulated events in the central region in Fig. 3. The relative size of the particle-level balance correction with respect to the measured resolutions is of the order of 5 %.
7 In situ jet resolution measurement using the bisector method

Bisector rationale
The bisector method [14] is based on a transverse balance vector, P T , defined as the sum of the momenta of the two leading jets in dijet events, p T,1 and p T,2 . This vector is projected along an orthogonal coordinate system in the transverse plane, (ψ, η), where η is chosen in the direction that bisects φ 12 , the angle formed by p T,1 and p T,2 . This is illustrated in Fig. 4. For a perfectly balanced dijet event, P T = 0. There are of course a number of sources that give rise to significant fluctuations around this value, and thus to a non-zero variance of its ψ and η components, denoted σ 2 ψ and σ 2 η , respectively. At particle level, P part T receives contributions mostly from initial-state radiation. This effect is expected to be isotropic in (ψ, η), leading to similar fluctuations in both components, At calorimeter level, P calo T will further differ from zero due to detector effects. Its ψ component, P calo T,2ψ , can be decomposed into three contributions, where the first two terms correspond to fluctuations due to the detector p T resolution, and the last one to the particle jet imbalance. Taking the variance of the sum of these three independent terms yields where the following relations have been used Here σ (p T ) corresponds to σ (p T,1 ) σ (p T,2 ), as both jets have the same p T resolution since they belong to the same y region. A relation similar to Eq. (5) holds for the η component: Subtracting Eq. (6) from Eq. (5), and using σ part ψ = σ part η , yields where the fractional jet p T resolution, σ (p T )/p T , is expressed in terms of calorimeter observables only. The contribution from soft radiation and the underlying event is minimised by subtracting in quadrature σ η from σ ψ . If one of the leading jets (j ) belongs to the rapidity region being probed, and the other one (i) to a previously measured reference y region, then The dispersions σ ψ and σ η are extracted from Gaussian fits to the P Tψ and P Tη distributions in bins ofp T . There is no φ cut imposed between the leading jets, but it is implicitly limited by a p EM-scale T,3 < 10 GeV requirement on the third jet, as discussed in the next section. Figure  cut, for anti-k t jets with R = 0.6. The two leading jets are required to be in the same rapidity region, |y| < 0.8, while there is no rapidity restriction for the third jet. As expected, both components increase due to the contribution from soft radiation as the p T,3 cut is increased. Also shown as a function of the p EM-scale T,3 cut is the square-root of the difference between their variances, which yields the fractional momentum resolution when divided by 2 p 2 T cos φ . It is observed in Fig. 7  20 GeV for 160 ≤p T < 260 GeV. The same behaviour is observed for otherp T ranges. This cancellation demonstrates that the isotropy assumption used for the bisector method is consistent with the data over a wide range of choices of p EM-scale T,3 without the need for requiring an explicit φ cut between the leading jets. The precision with which it can be ascertained that the data is consistent with

Performance for the EM + JES calibration
The performances of the dijet balance and bisector methods are compared for both data and Monte Carlo simulation as a function of jet p T for jets reconstructed in the central region with the anti-k t algorithm with R = 0.6 and using the EM + JES calibration scheme. The results are shown in Fig. 8. The resolutions obtained from the two independent in situ methods are in good agreement with each other within the statistical uncertainties. The agreement between data and Monte Carlo simulation is also good within the statistical precision.
The resolutions for the three jet rapidity bins with |y| > 0.8, the Extended Tile Barrel, the Transition and the End-Cap regions, are measured using Eqs. (3) and (8), taking the central region as the reference. The results for the bisector method are shown in Fig. 9. Within statistical errors the resolutions obtained for data and Monte Carlo simulation are in agreement within ±10% over most of the p T -range in the various regions. Figure 9 shows that dependences are well described by fits to the standard functional form expected for calorimeterbased resolutions, with three independent contributions, the effective noise (N ), stochastic (S) and constant (C) terms. The N term is due to external noise contributions that are not (or only weakly) dependent on the jet p T , and include the electronics and detector noise, and contributions from pile-up. It is expected to be significant in the low-p T region, below ∼30 GeV. The C term encompasses the fluctuations that are a constant fraction of the jet p T , assumed at this early stage of data-taking to be due to real signal lost in passive material (e.g. cryostats and solenoid coil), to nonuniformities of response across the calorimeter, etc. It is expected to dominate the high-p T region, above 400 GeV. For intermediate values of the jet p T , the statistical fluctuations, represented by the S term, become the limiting factor of the resolution. With the present data sample that covers a restricted p T range, 30 GeV ≤ p T < 500 GeV, there is a high degree of correlation between the fitted parameters and it is not possible to unequivocally disentangle their contributions.

Closure test using Monte Carlo simulation
The Monte Carlo simulation expected resolution is derived considering matched particle and calorimeter jets in the event, with no back-to-back geometry requirements. Matching is done in η-φ space, and jets are associated if R = ( η) 2 + ( φ) 2  correspond to the transverse momentum of the reconstructed jet and its matched particle jet, respectively. The jet response distribution is modelled by a fitted Gauss distribution, and its standard deviation is defined as the truth jet p T resolution.
The Monte Carlo simulation truth jet p T resolution is compared to the results obtained from the dijet balance and  Fig. 10. This comparison will be referred to as the closure test. The in situ and truth resolutions agree within 10 %, with the truth results typically 10 % lower. This result confirms the validity of the physical assumptions discussed in Sects. 6 and 7 and the inference that the observables derived for the in situ MC dijet balance and bisector methods provide reliable estimates of the jet energy resolution. The systematic uncertainties on these estimates are of the order

Jet energy resolution uncertainties
There are three kind of systematic uncertainties to be considered. Section 10.1 discusses the experimental uncertainties that affect the in situ measurements. Section 10.2 addresses the method uncertainties, that is the precision with which the in situ methods in data describe the truth resolution. Finally, Sect. 10.3 studies the truth resolution uncertainty due to event modelling in the Monte Carlo simulation.

Experimental in situ uncertainties
The squares (circles) in Fig. 11 show the experimental relative systematic uncertainty in the dijet balance (bisector) method as a function ofp T . The different contributions are discussed below. The shaded area corresponds to the larger of the two systematic uncertainties for eachp T bin.
For the dijet balance method, systematic uncertainties take into account the variation in resolution when applying different φ cuts (varied from 2.6 to 3.0), resulting in a 2-3 % effect for 30 ≤ p T < 60 GeV, and when varying the parameterization of K soft (p T ) (see Sect. 6.2), which contributes up to 6 % at p T ≈ 30 GeV. For the bisector method, the relative systematic uncertainty is about 4-5 %, and is derived from the precision with which it can be verified that The contribution from the JES uncertainties [39] is common to both methods. It is 1-2 %, determined by recalculating the jet resolutions after varying the JES within Fig. 11 The experimental systematic uncertainty on the dijet balance (squares) and bisector (circles) methods as a function ofp T , for jets with |y| < 0.8. Also shown is the absolute value of the relative difference between the two methods in each p T bin for data (dot-dashed lines) and for Monte Carlo simulation (dashed lines) its uncertainty in a fully correlated way. The resolution has also been studied in simulated events with added pile-up events (i.e. additional interactions as explained in Sect. 3.3), as compared to events with one hard interaction only. The sensitivity of the resolution to pile-up is found to be less than 1 % for an average number of vertices per event of 1.9.
In summary, the overall relative uncertainty from the in situ methods decreases from about 7 % at p T = 30 GeV down to 4 % at p T = 500 GeV. Figure 11 also shows the absolute value of the relative difference between the two in situ methods, for both data and Monte Carlo simulation. They are found to be in agreement within 4 % up to 500 GeV, and consistent with these systematic uncertainties.

Uncertainties on the measured resolutions
The uncertainties in the measured resolutions are dominated by the systematic uncertainties, which are shown in Table 1 as a percentage of the resolution for the four rapidity regions and the two jet sizes considered, and for characteristic ranges, low (∼ 50 GeV), medium (∼ 150 GeV) and high (∼ 400 GeV) p T . The results are similar for the four calibration schemes.
The dominant sources of systematic uncertainty are the closure and the data/MC agreement. The experimental systematic uncertainties, discussed in Sect. 10.1, are significantly smaller. The closure uncertainty (see Sect. 9), defined as the precision with which in simulation the resolution determined using the in situ method reproduces the truth jet resolution, is larger for R = 0.4 than for R = 0.6, smaller at high p T than at low p T , and basically independent of the rapidity. The data/MC agreement uncertainty, the precision with which the MC simulation describes the data, is independent of R, larger at low and high p T than at medium p T , and it grows with rapidity because of the increasingly limited statistical accuracy with which checks can be performed to assess it.  The systematic uncertainties in Table 1 for jets with R = 0.4 are dominated by the contribution from the closure test. They decrease with increasing p T and are constant for the highest three rapidity bins. The systematic uncertainties for jets with R = 0.6 are consistently smaller than for the R = 0.4 case, and receive comparable contributions from closure and data/MC agreement. They tend to increase with rapidity and are slightly lower in the medium p T range. The uncertainty increases at high p T for the end-cap, 2.1 ≤ |y| < 2.8, because of the limited number of events in this region.

Uncertainties due to the event modelling in the Monte Carlo generators
Although not relevant for the in situ measurements of the jet energy resolution, physics analyses sensitive to the expected resolution have to consider its systematic uncertainty arising from the simulation of the event. The expected jet p T resolution is calculated for several Monte Carlo simulations in order to assess its dependence on different generator models (ALPGEN and HERWIG++), PYTHIA tunes (PERUGIA2010), and other systematic variations (PARP90; see Sect. 3.1). Differences between the nominal Monte Carlo simulation and PYTHIA8 [24] have also been considered. These effects, displayed in Fig. 12, never exceed 4 %. The total modelling uncertainty is estimated from the sum in quadrature of the different cases considered here. This is shown by the shaded area in Fig. 12 and found to be at most 5 %.

Jet energy resolution for other calibration schemes
The resolution performance for anti-k t jets with R = 0.6 reconstructed from calorimeter topological clusters for the Local Cluster Weighting (LCW + JES), the Global Cell Weighting (GCW + JES) and the Global Sequential (GS) calibration strategies (using the bisector method) is presented in Fig. 13 for the Central, Extended Tile Barrel, Transition and End-Cap regions. The top panel shows the resolutions determined from data, whereas the bottom part compares data and Monte Carlo simulation results. The three more sophisticated calibration techniques improve the resolution σ (p T )/p T with respect to the EM + JES calibrated jets by approximately 0.02 over the whole p T range. The relative improvement ranges from 10 % at low p T up to 40 % at high p T for all four rapidity regions. Figure 14 displays the resolutions for the two in situ methods applied to data and Monte Carlo simulation for |y| < 0.8 (left plots). It can be observed that the results from the two methods agree, within uncertainties. The Monte Carlo simulation reproduces the data within 10 %. The figures on the right show the results of a study of the closure for each case, where the truth resolution is compared to that obtained from the in situ methods applied to Monte Carlo simulation data. The agreement is within 10 %. Overall, comparable agreement in resolution is observed in data

Improvement in jet energy resolution using tracks
The addition of tracking information to the calorimeterbased energy measurement is expected to compensate for the jet-by-jet fluctuations and improve the jet energy resolution (see Sect. 5.5).
The performance of the Track-Based Jet Correction method (TBJC) is studied by applying it to both the EM + JES and LCW + JES calibration schemes, in the central region. The measured resolution for anti-k t jets with R = 0.6 (R = 0.4) is presented as a function of the average jet transverse momentum in the top (bottom) plot of Fig. 15. The relative improvement in resolution due to the addition of tracking information is larger at low p T and more important for the EM + JES calibration scheme. It ranges from 22 % (10 %) at low p T to 15 % (5 %) at high p T for the EM + JES (LCW + JES) calibration. For p T < 70 GeV, jets calibrated with the EM + JES + TBJC scheme show a similar performance to those calibrated with the LCW + JES + TBJC scheme. Overall, jets with LCW + JES + TBJC show the best fractional energy resolution over the full p T range.

Summary
The jet energy resolution for various JES calibration schemes has been measured using two in situ methods with a data sample corresponding to an integrated luminosity of 35 pb −1 collected in 2010 by the ATLAS experiment at √ s = 7 TeV. The Monte Carlo simulation describes the jet energy resolution measured in data within 10 % for jets with p T values between 30 GeV and 500 GeV in the rapidity range |y| < 2.8.
The resolutions obtained applying the in situ techniques to Monte Carlo simulation are in agreement within 10 % with the resolutions determined by comparing jets at calorimeter and particle level. Overall, the results measured with the two in situ methods have been found to be consistent within systematic uncertainties.