Single hadron response measurement and calorimeter jet energy scale uncertainty with the ATLAS detector at the LHC

The uncertainty on the calorimeter energy response to jets of particles is derived for the ATLAS experiment at the Large Hadron Collider (LHC). First, the calorimeter response to single isolated charged hadrons is measured and compared to the Monte Carlo simulation using proton-proton collisions at centre-of-mass energies of sqrt(s) = 900 GeV and 7 TeV collected during 2009 and 2010. Then, using the decay of K_s and Lambda particles, the calorimeter response to specific types of particles (positively and negatively charged pions, protons, and anti-protons) is measured and compared to the Monte Carlo predictions. Finally, the jet energy scale uncertainty is determined by propagating the response uncertainty for single charged and neutral particles to jets. The response uncertainty is 2-5% for central isolated hadrons and 1-3% for the final calorimeter jet energy scale.


Introduction
Partons scattered in proton-proton interactions are measured with the ATLAS detector as collimated jets of hadrons.The uncertainty on the jet energy scale is the largest source of detector-related systematic uncertainty for many physics analyses carried out by the ATLAS Collaboration, from the di-jet cross-section and top mass measurements, to searches for new physics with jets in the final state.It is thus the subject of an extensive and detailed study [1].
The jet energy measured by the calorimeter is corrected for calorimeter non-compensation and energy loss in dead material.The corresponding jet energy scale CERN, Geneva, Switzerland correction factor is referred to as the JES.The JES is derived from Monte Carlo (MC) simulations by comparing the calorimeter energy of an isolated reconstructed jet to that of the particle jet1 that points to it [1].
The uncertainty on the calorimeter energy response is a significant component of the total uncertainty on the JES.It is derived in this paper by convolving the measured uncertainty on the single charged hadron energy response and the estimated uncertainty on the neutral particle energy response with the expected particle spectrum within a jet.
The calorimeter response to single isolated charged hadrons, and the accuracy of its Monte Carlo simulation description, can be evaluated from the ratio of the calorimeter energy E to the associated isolated track momentum p.The aim of the measurement is to estimate the systematic uncertainty on jet calorimeter response and therefore the focus is on data-to-MC comparison.The ratio E/p is measured using proton-proton collisions at centre-of-mass energies of √ s = 900 GeV and 7 TeV over a wide range of track momenta in the central region of the calorimeter.Possible additional uncertainties introduced by certain particle species are addressed by measuring the response to hadrons identified through the reconstruction of known short-lived particles.
This paper is organised as follows.The relevant features of the ATLAS detector are summarised in Section 2. Then the measurement of the single charged hadron response, at centre-of-mass energies of √ s = 900 GeV and 7 TeV, is discussed in Section 3. Section 4 covering 1.5 < |η| < 3.2) are LAr sampling calorimeters with copper absorbers.The copper/tungsten-LAr forward calorimeters (FCal) provide both electromagnetic and hadronic energy measurements and extend the coverage to |η| < 4.9.
The data used for this analysis were triggered using the minimum bias trigger scintillators (MBTS) [3].The MBTS are mounted at each end of the detector in front of the LAr end-cap calorimeter cryostats at z = ±3.56m and are segmented into eight sectors in azimuth and two rings in pseudorapidity (2.09 < |η| < 2.82 and 2.82 < |η| < 3.84).

Single particle response for charged hadrons
In this section, the response of the calorimeters to isolated charged hadrons is compared to the predictions from MC simulation.The ratio of the energy, E, deposited by an isolated charged particle in the calorimeter to the track momentum, p, is used to evaluate the uncertainty on the calorimeter response modelling in the MC simulation.

Event selection
Events are required to have at least one reconstructed vertex with at least four associated tracks.The total number of events satisfying this selection is about 25 million collected in 2009 (approximately one million events at √ s = 900 GeV) and 2010 (approximately 24 million events at √ s = 7 TeV).For every selected event, each track candidate is extrapolated to the second longitudinal layer of the EM calorimeter.A track is defined as isolated if its impact point has a distance ∆R = (∆η)2 + (∆φ) 2 > 0.4 from all other track candidate impact points.The isolation conditions are studied in detail in Section 4. 5.
The isolated tracks must also have: a transverse momentum of p T > 500 MeV3 , -a minimum of one hit in the pixel detector and six hits in the SCT, and small transverse and longitudinal impact parameters (defined in Ref. [4]) computed with respect to the primary vertex 4 , |d 0 | < 1.5 mm and |z 0 | sin θ < 1.5 mm.
The above requirements ensure a good quality of the track and reduce contributions from fake tracks to a negligible level.

Monte Carlo simulation
A sample of about 10 (20) million non-diffractive protonproton collision events at √ s = 900 GeV ( √ s = 7 TeV) are generated using Pythia 6.421 [5] with the ATLAS minimum bias tune 1 (AMBT1) [4].For isolated, highmomentum (p T > 15 GeV) tracks this corresponds to ∼ 60% of the available events in data.All the events are run through a full detector simulation [6] based on Geant4 [7].The set of Geant4 physics models used is QGSP BERT [8].The reconstruction and analysis software used for the MC simulation is the same as for the data.

Defining the E/p observable
The sum of the energy deposits in layers of calorimeter cells associated to a selected track is computed using topological clusters [9] at the electromagnetic scale, i.e., without applying any correction for the calorimeter non-compensation or for energy loss in dead material.
The purpose of the topological clustering algorithm is to identify areas of connected energy deposits in the calorimeter, based on the significance of the energy deposits in cells with respect to the expected noise level.Topological clusters are formed around cells with energy |E cell | > 4σ noise ("seeds"), where σ noise is the RMS of that cell noise.Then, iteratively, the cluster is expanded by adding all neighbouring cells with |E cell | > 2σ noise .Finally, the cells surrounding the resulting cluster are added, regardless of their energy.The η − φ position of a cluster i in a given calorimeter layer j, (η ij cl , φ ij cl ) is computed as the energy-weighted position of the cells in layer j belonging to the cluster.
The position of the track k extrapolated to the layer j is (η kj tr , φ kj tr ).The energy of a cluster in the layer j (E j ) is associated to the track if: The parameter R coll is set to 0.2 based on a tradeoff between maximising the particle shower containment and minimising the background contribution coming from neutral particles produced close to the track.Roughly 90% of the shower energy is collected in a cone of such size.
The energy E associated to a track is computed as the sum of the associated energy deposits in all layers, E = j E j , and the ratio E/p is formed with the reconstructed track momentum.Note that because of calorimeter noise fluctuations, E (and therefore E/p) can assume negative values.

E/p distributions
The E/p distributions in two representative regions of η and track momentum are shown in Fig. 1.The large number of entries with E/p = 0 corresponds to isolated tracks that have no associated cluster in the calorimeter.Several effects may be responsible for this: -Particles can interact hadronically before reaching the calorimeter (in the ID, cryostat or solenoid magnet).Such particles can change their direction, or produce a large number of low momentum secondary particles.-A cluster is created only if a seed is found.Hadrons with low momentum and an extended shower topology sometimes do not have a single cell energy deposit large enough to seed a topological cluster.
The cases where the calorimeter response is compatible with zero have been further studied.The probability that the calorimeter response is suppressed by noise threshold requirements, P (E = 0), is shown in Fig. 2a as a function of the amount of material (in nuclear interaction lengths) in front of the active volume in the central (|η| < 1.0) calorimeter region.When the hadron passes through more material, the probability that no energy is associated to the track increases.
Figure 2b shows P (E = 0) as a function of the track momentum in the central region of the calorimeter (|η| < 0.6).In this region, the dead material in front of the calorimeter is approximately constant.The probability decreases with increasing track momentum.In general, P (E = 0) is well predicted by the MC simulation.

Background subtraction
The energy measured inside the cone of ∆R < R coll = 0.2 centered around the track impact point may be contaminated by energy deposits from the showers of close-by particles produced in the proton-proton collision.The track isolation requirement suppresses possible shower contamination from charged particles.There is no obvious way to suppress shower contamination from photons, mostly produced in π 0 → γγ decays, and neutral hadrons.The neutral particle background contribution to the E/p measurement in the MC simulation depends on the event generator settings of the parameters governing non-perturbative QCD processes and on the modelling of the calorimeter response to low momentum neutral particles, and it is therefore difficult to model correctly.The neutral background is thus subtracted from the measured response using an in situ background estimate.In the following, the expression "EM (HAD) energy" will refer to the energy deposited in the electromagnetic (hadronic) calorimeter.
The background subtraction relies on the assumption that the EM energy from photons and neutral hadrons is independent of the energy deposited by the selected track.Charged hadrons are selected that behave like minimum ionising particles in the EM calorimeter and start their shower in the hadronic calorimeter (late-showering hadrons).Excluding a narrow region around the late-showering hadron track, the remaining EM energy is mainly due to showers from neutral particles.The strategy is sketched in Fig. 3.
Late-showering hadrons are selected by requiring a small amount of EM energy in a cone of ∆R < 0.1, E 0.1 EM < 1.1 GeV, and a large HAD energy fraction, E 0.1 HAD /p > 0.4.The background is measured in the EM calorimeter in an annulus around the late-showering charged hadrons.The mean of the background distribution over many events in a given momentum and pseudorapidity bin estimates the energy deposition of pho- tons and neutral hadrons showering in the EM calorimeter: where E 0.2 EM is the EM energy in a cone of ∆R < 0.2.The background contribution from neutral hadrons depositing their energy in the hadronic calorimeter was estimated with a similar technique applied to different hadronic calorimeter layers and found to be negligibly small.
Figures 4 and 5 show E/p BG as a function of p in two bins of pseudorapidity at √ s = 900 GeV and √ s = 7 TeV.The background E/p BG at √ s = 7 TeV in both bins is about 0.04−0.08 for p < 10 GeV (0.02−0.08 at √ s = 900 GeV) and decreases to ∼ 0.04 for track momenta of ∼20 GeV.The general trend is confirmed by the MC simulation, although the two show some significant differences.The discrepancy is attributed to an imperfect modelling of non-perturbative QCD processes by the Pythia settings used for the MC event simulation.
The quantity that is studied in detail in the next section as a function of the track momentum and calorimeter impact point pseudorapidity is: where E/p raw is the mean value of the E/p distribution before background subtraction (cf.Fig. 1).The background (Eq.2) is rescaled by the factor of 4/3, which is the ratio of the area of the full ∆R < 0.2 cone to that of an annulus with 0.1 ≤ ∆R < 0.2, to take into account the background contribution in E 0.1 EM , since the background is only measured in this annulus.A uniform energy density of the background in the ∆R < 0.2 cone is assumed.The assumption has been validated by comparing its prediction with that of a more complex procedure, which is described in Appendix A.
The background subtraction procedure is applied in the following to all E/p measurements, both at √ s = 900 GeV and √ s = 7 TeV.Because the background contribution increases with √ s, the JES uncertainty estimation relies on the 900 GeV data for the measurement of calorimeter response to tracks with momenta below 2.2 GeV.

Results
Several systematic uncertainties on the measurement are estimated.Each is taken to be completely correlated between all pseudorapidity and momentum bins in the E/p measurement: -Track selection: The dependence of E/p on the track selection requirements in Section 3.1 has been estimated by varying the number of silicon tracker hits required and the impact parameter selection with respect to the primary vertex within reasonable ranges.The MC-to-data ratio of E/p was found to be almost unaffected by variations in the track selection.The maximum variation found in the ratio (0.5%) is taken as a systematic uncertainty.-Track momentum scale: The uncertainty on the momentum scale p as measured by the inner detector is negligibly small for p < 5 GeV [10].For p > 5 GeV, a conservative 1% uncertainty has been assumed on the momentum scale.-Background subtraction: The difference between the background estimate obtained with the method described in Section 3.5 and with the validation method described in Appendix A is taken as a systematic uncertainty.This results in a 1% uncertainty on the E/p measurements.
The mean E/p value after background subtraction is evaluated in bins of momentum and pseudorapidity.Figures 6 and 7 show E/p as a function of the track momentum, in two different |η| bins up to |η| = 1.1, at √ s = 900 GeV and √ s = 7 TeV.The lower parts of the figures present the ratio of MC simulation to data.The data with √ s = 900 GeV and √ s = 7 TeV agree within the statistical uncertainty.The maximum momentum that can be probed with the data considered is approximately 30 GeV.The agreement between data and MC simulation is within ∼2% for particles with momenta in the 1-10 GeV range, and it is around 5% for momenta in the 10-30 GeV range.Below 1 GeV, where tracks are just at the kinematic threshold of entering the calorimeter volume, large differences of ∼10% or more between  The black markers represent the background estimated from collision data, while the green rectangles represent the MC prediction, with the vertical width showing its statistical uncertainty.The lower panes show the ratio of the MC prediction to collision data.The MC prediction of the background energy deposit is obtained using the same procedure as applied to the data.data and MC simulation are visible.However, due to the low absolute calorimeter response to very low momentum particles, these differences are not critical for the JES determination.

Calorimeter response to identified hadrons
The extrapolation of the previous single particle response studies into the environment of a jet requires understanding of two additional effects.A jet includes a variety of hadrons that may differ from the inclusive sample of isolated hadrons.Therefore, measuring the average response to different species of particles is valuable to ensure that the Monte Carlo correctly models all aspects of the jet shower.Additionally, the hadrons in a jet are not isolated.Threshold effects and hadronic shower widths affect the calorimeter response to the multi-hadron system.In order to address these two points, this study adds two additional features to help complete the understanding of calorimeter response.To minimise the impact of the background and to simplify its estimation, the measurements are presented as a difference between the ratio E/p for two different particle types or as a ratio to the inclusive measurement.
First, single hadrons are identified using decays of K S (for positive and negative pions), Λ (for protons), and Λ (for anti-protons) particles.These single hadrons are required to be isolated from all other charged particles in the event.Single pions will have manifestly lower energy response distributions from those of antiprotons, because of the eventual annihilation of the anti-proton.By identifying and isolating single pions, single protons, and single anti-protons, the effects of hadronic interactions and annihilation can be separated at low to moderate energies, where they are most important.
Second, when the mother particle is highly boosted, the decay products are more collimated.For track momenta of ∼2-6 GeV from K S decays, the range of opening angles of the decay products allows a measurement of calorimeter response as more tightly collimated pairs of pions are selected.There are two related effects probed by such a measurement.Energy deposited in the calorimeter may fall below the thresholds for reconstruction and be neglected as consistent with noise.As two showers overlap, the addition of energies that might have individually been below threshold can produce a signal above these noise thresholds.The increase of the signal is related to both the thresholds themselves and the width of the hadronic shower.For this study, the π + π − system is required to be isolated from other charged particles in the calorimeter.

Event selection and observable definition
The event selection for this measurement follows closely the event selection of the inclusive measurement already presented, except that a secondary vertex is required in the event.The events are collected using random triggers at √ s = 7 TeV and correspond to an integrated luminosity of about 800 µb −1 .The same primary vertex requirements are applied, and the same requirements are placed on tracks entering the measurement of E/p, except the requirements on the tracks' impact parameters.The Monte Carlo used is described in Section 3.2.
The definition of the ratio E/p and the isolation requirement for tracks is described in Section 3.3.
The ratio E/p is presented as a function of available energy, E a .For pions, this is simply the particle's total energy: E a = p 2 + m 2 .For protons, only the kinetic energy is included: For anti-protons, the kinetic energy plus double the restmass is included, in order to take into account annihilation, E a = p 2 + m 2 + m.The available energy is therefore calculated using the information from the tracker.

Reconstruction of short-lived particle candidates
The reconstruction and selection of long-lived particles is based on previous ATLAS results [10].The decay K S → π + π − , which dominates the K S decays to other charged particles by several orders of magnitude, is used to identify pions.The decay Λ → π − p (Λ → π + p), which dominates the Λ (Λ) decays to other charged particles by several orders of magnitude, is used to identify protons (anti-protons).Both decay product tracks are required to have p T > 100 MeV, and the tracks used for the E/p measurement must have p T > 500 MeV.In the case of Λ (Λ) candidate decays, the charge of the higher momentum track is used to label the candidate as a Λ (positively-charged higher momentum track) or Λ (negatively-charged higher momentum track), as is kinematically favored.The tracks entering the E/p distributions are additionally required to have a pseudorapiditydependent number of hits in the TRT, the outermost tracking system.This requirement suppresses the effect of nuclear interactions in the material of the inner detector, particularly in the outer layers of the silicon tracker.
The individual tracks entering the response distributions are divided into two bins of pseudorapidity, |η| < 0.6 and 0.6 ≤ |η| < 1.1.Several higher pseudorapidity bins show consistent results, albeit with significantly lower statistics.Example distributions of reconstructed mass for candidates with at least one cen-  The signal purity for a given pseudorapidity interval, where the binning is always in terms of the kinematic properties of the track of interest, is estimated using a fit to the signal over a cubic polynomial fit to the background, following Ref.[10].For the K S , a double Gaussian signal function is used, constrained to have identical central values for both Gaussians, and for the Λ and Λ, a modified Gaussian is used of the form where , m is the reconstructed mass, and a, b and c are free parameters of the fit.The peak mean, b, is stable to within a few hundred keV over all pseudorapidity bins.The fitted width of the mass peak, c, increases with track pseudorapidity, due to the track resolution.Example fit results for the Λ are shown in Fig. 8 for positively charged tracks (proton candidates) and negatively charged tracks (pion candidates) in the central pseudorapidity bin for the E/p measurement.The negatively charged tracks are softer, which is reflected in the smaller statistics for the negatively charged track mass distribution.
For each mass peak and each bin of track pseudorapidity, an acceptance window is constructed in order to optimise both signal purity and statistics.The width of the window is set at three times the width of the narrower Gaussian for the K S and three times the width of the modified Gaussian for the Λ (Λ).The purities are over 97% for the K S candidates and above 92% for the central Λ candidates.The purities in the MC simulation are within 2% of those in the data.The number of tracks from signal candidates, extracted from the fits, in the data and 20 million MC simulation events are shown in Table 1.Roughly twice as many K S candidates and three to four times more Λ and Λ candidates are found in the data than in the MC simulation, because of the limited MC simulation statistics.The difference in K S and Λ yields between data and MC simulation are discussed further in Ref. [10], and the response measurements here are insensitive to these differences.Only the positively and negatively charged tracks that form short-lived particle candidates within the defined acceptance window are considered in the remainder of this section.
Table 1: The number of signal candidate tracks, extracted from the fits described in the text, in data and MC simulation.
The ATLAS Collaboration

Background subtraction
There are three sources of charged backgrounds entering the E/p distributions.Nuclear interactions in the material of the inner detector that fake short-lived particle candidates are suppressed by the narrow mass acceptance window.The smoothness of the distribution of secondary vertex positions confirms that these are at the level of a few percent.The tracks may also undergo nuclear interactions prior to entering the calorimeter.These interactions are suppressed by requiring that the track have hits in the TRT.An additional charged background comes from combinatorics, particularly in high track-multiplicity events.Because the purity of K S candidates is typically 5-10% higher than that for Λ or Λ candidates, and because the majority of the background is from charged pions, charged background corrections are applied only for calculating proton and anti-proton response.The response to pions is taken from the K S candidates and is used to correct the response to protons and anti-protons.When charged pions enter the proton energy response distribution, they are given an incorrect track mass hypothesis.The track extrapolation to the calorimeter does not calculate energy loss differently for pions and (anti-)protons, so that the measured associated energy in the calorimeter remains the same.However, the available energy changes significantly when the different mass hypothesis is used.Thus, in order to subtract the charged pion background from the (anti-)proton distribution, the pions from K S candidates are given the (anti-)proton mass hypothesis and their response is re-calculated.The (anti-)proton response in a bin i of available energy and η, E/p i , is then given by where ǫ i is the purity of the sample in that pseudorapidity bin, E/p raw i is the measured response to the (anti-)protons, and E/p π i is the measured response to the charged pions with a proton mass hypothesis.
The corrections are small in the region |η| < 1.1, falling from ∼5% at low available energy to ∼1% at high available energy for both protons and anti-protons.The corrections are derived independently for data and MC simulation to ensure that the differences in purity and any differences in response are taken into account.
As the momenta of the daughter tracks increase, the position resolution of the secondary vertex broadens.This introduces a larger combinatorial background in the high-momentum bins.The purities in MC simulation follow, to ∼5%, those of the data with increasing momentum.The systematic uncertainty on the corrected value of E/p introduced by the p T -dependence of the purity is determined to be well below 1% and is therefore neglected.
There is an additional contribution to the (anti-) proton response measurement for Λ (Λ) decays faking Λ (Λ) decays.Particularly when the Λ or Λ has low energy, the pion may be the higher-momentum track.These "fakes" are suppressed by the kinematic cuts applied in the candidate selection, but are still present at some level, particularly in the lowest bin of available energy.This background is dominated by well-understood two-body decay kinematics.It should be well-described by the MC simulation and is therefore taken into account using MC predictions.
As discussed in Section 3.5, there is an additional contribution to E/p from neutral particles in the event, since isolation is only defined relative to charged particles.Only differences in response between particle species are reported here, because the neutral background should cancel in the difference (except insofar as the K S or Λ production occurs near additional activity).The charged particle isolation criterion should be sufficient to ensure that the neutral background is uncorrelated with any jet-like activity in the event.The cancellation of the neutral background is tested using a simulation sample of single particles and is found to be valid up to the available statistical accuracy.Therefore, no additional correction or uncertainty is added for this background.Thus, background systematic uncertainties are found to be negligible with respect to the statistical errors of the MC simulation sample.

Isolated identified single particle response
The uncorrected distribution of E/p for π − , p, π + and p in a single bin of available energy and pseudorapidity, 2.2 ≤ E a < 2.8 GeV and |η| < 0.6, is shown in Fig. 9 (cf.Fig. 1 for inclusive hadrons in a different p range).All the distributions have a small negative tail from noise in the calorimeter and a long positive tail from the neutral background.There is a much more prominent positive tail in the p response distribution, where annihilation plays a significant role.A significant fraction of tracks have E = 0, for which no cluster of energy was found near the track in the calorimeter, as discussed in Section 3.4.The average response, E/p , is defined as the arithmetic mean of the full distribution.The response of protons and positively charged pions is well modeled in the MC simulation.The response to pions is significantly lower than that to anti-protons in this available energy range, although fewer pions have E = 0.The agreement in the fraction of pions with E = 0 builds confidence in the modelling of the material in front of the calorimeter.In order to reduce the effects of neutral background subtraction and focus on the differences in response amongst particle species, the difference in response between several pairs of particles is measured.Figure 10 shows the difference in response between π + and π − in the central (|η| < 0.6) and forward (0.6 ≤ |η| < 1.1) pseudorapidity bins.The response to π + is higher than that to π − at low available energy, in agreement with Ref. [11], wherein this difference is attributed to a material dependent charge-exchange effect.The difference may also be related to the "Barkas correction" described in Ref. [12].The MC simulation is, however, reasonably consistent with the data.
Figure 11 shows the difference in response between π + and p in the central (|η| < 0.6) and forward (0.6 ≤ |η| < 1.1) pseudorapidity bins.Again, the difference between π + and proton response in MC simulation is consistent with the data over the entire range of available energies.At low available energy, a larger average fraction of the initial hadron energy is converted into an electromagnetic shower for pions than for protons [12].This leads to a lower response for protons, particularly at low available energy, due to the non-compensation of the calorimeter.
Figure 12 shows the difference in response between π − and p in the central (|η| < 0.6) and forward (0.6 ≤ |η| < 1.1) pseudorapidity bins.The difference shows a ∼30%× E/p disagreement between data and MC simulation for E a < 3 GeV, though they are consistent for E a > 4 GeV.At these low available energies, the anti-proton response should be dominated by the annihilation and the subsequent shower.The difference indicates large contributions from processes not wellmodelled by Geant4.To test the contribution from the calorimeter acceptance to this effect, the response is constructed excluding the tracks with E = 0.The pions and protons show the same level of agreement between data and MC simulation, and the disagreement in the response difference between pions and anti-protons remains near 10%.

Response to nearby particles
To study the effects of calorimeter thresholds and noise suppression as more tightly collimated pairs of particles are selected, isolation is not required with respect to the other daughter of the particle candidate.That is, when constructing the energy response for the π + from a K S decay, the π + is required to be isolated from all tracks in the event other than the π − from the K S decay.The response is then examined as a function of the distance between the two pions in η − φ after extrapolation to the second layer of the electromagnetic calorimeter.No isolation criteria are applied to the π − .Because the pions are required to come from a particle with a narrow width, the kinematics of the second pion are constrained, thus removing one possible source of discrepancy between data and MC simulation.
The ratio E/p as calculated here for a single pion should increase as another particle approaches it and more of the second particle's energy is included in the calorimeter energy, E. As shown in Fig. 13 and 14, the response does not vary with separation for large extrapolated distances, where the single pion is isolated.The distributions are normalised to the average E/p in the same η and p bin in order to remove any differences from the average response.In Fig. 14, pions of both charges are considered in order to increase the available statistics.The response rises below ∆R ≈ 0.3, which confirms that the isolation criterion used in Sections 3.1 and 4.4 does not lead to a bias.It increases significantly for ∆R < 0.2, when the second pion is inside the cone in which the energy is counted when constructing E/p.The height of the peak at low ∆R is related to the kinematics of the K S decay.In particu-lar, some of the low ∆R and low available energy bins contain mostly asymmetric decays.
The height of the response distribution for closeby pions and small opening angles is connected to the threshold effects, since it is a measurement of the out-ofcluster energy, to first order.The slope of the response curve is closely related to the single-hadron shower width.The agreement between data and MC is good over most of the range of opening angles, but the peak at low opening angles is wider in data, indicating a somewhat broader shower.This discrepancy is qualitatively Fig. 13: E/p for pions as a function of the extrapolated distance between that pion and a pion of the opposite sign.In all cases the pions are daughters of a reconstructed K S candidate and are required to be isolated from all tracks in the event except the other daughter of the decay.The response is shown for low energy (2.2 ≤ E a < 2.8 GeV) (a) π + and (b) π − in the central pseudorapidity bin (|η| < 0.6).The MC simulation has no tracks passing the selection in the smallest bin (smallest two bins) of opening distance for π + (π − ). in agreement with the discrepancies observed in lateral hadronic shower shape, seen in test-beam studies and elsewhere [8,13,14,15].

Calorimeter jet energy scale uncertainty
The jet energy scale calibration (JES) corrects the measured jet energy for several effects, including calorime-ter non-compensation and energy loss in dead material.The calibration itself is derived from MC simulation.The calorimeter uncertainty on the JES is calculated from the uncertainty on the energy response of all particles contributing to a jet.Within the MC simulation, the energy contribution of each individual particle to a given jet can be separated.The convolution of the uncertainty on the single particle energy response with the Fig. 14: E/p for pions as a function of the extrapolated distance between that pion and a pion of the opposite sign.In all cases the pions are daughters of a reconstructed K S candidate and are required to be isolated from all tracks in the event except the other daughter of the decay.The response is shown for (a) more forward (0.6 ≤ |η| < 1.1) pions with 2.2 ≤ E a < 2.8 GeV and for central (|η| < 0.6) pions in three bins of available energy: (b) 2.8 ≤ E a < 3.6 GeV, (c) 3.6 ≤ E a < 4.6 GeV and (d) 4.6 ≤ E a < 6.0 GeV.The MC simulation has no tracks passing the selection in the smallest opening distance bin.
MC jet particle composition is then used to calculate the calorimeter uncertainty on the jet energy scale.
The calorimeter JES uncertainty is derived for the well-understood central region of the calorimeter.The main reasons for the restriction to the central calorimeter region are the smaller amount of material in front of the calorimeter and the existence of combined test beam measurements with a setup very similar to the final ATLAS configuration.Therefore, the calorimeter JES uncertainty is only evaluated using the single particle response for |η| < 0.8.The total JES uncertainty, including all calorimeter regions and all effects not re-lated to the calorimeter energy response, is discussed in Ref. [1].
The JES uncertainty is determined for jets reconstructed from topological clusters [9] with the anti-k t jet algorithm [16] implemented in the FASTJET package [17] for R = 0.4 and R = 0.6 jet sizes.
The analysis is performed with inclusive di-jet Monte Carlo events simulated with Pythia.Jets are selected requiring a separation of ∆R > 2.0 to any other jet with EM scale 5 p T (EM) > 7 GeV.The jet transverse mo-mentum is calibrated from the EM scale to the hadronic scale using a MC-based JES [1].
The numerical evaluation of the uncertainty on the jet energy scale is performed with Monte Carlo pseudoexperiments.In each pseudo-experiment, the jet energy scale is calculated after randomly changing the Monte Carlo single particle energy response within the appropriate uncertainty range given by the measured data/MC ratio.The final uncertainty on the jet energy scale is then given by the spread of the distribution of the jet energy scale over all pseudo-experiments.Within each pseudo-experiment, all randomly changed parameters are kept fixed, such that the energy response correlations are properly taken into account.
The uncertainties on the particle energy response functions entering the calculation are taken from the E/p measurements described in Sections 3 and 4, ATLAS combined test beam (CTB) measurements [13] and Geant4 Monte Carlo predictions.The details are described in the following subsections.

Additional E/p uncertainties contributing to the jet energy scale
When the single charged hadron response from E/p is used to assess the uncertainty on the jet energy scale, further systematic uncertainties that might affect the propagation of the response to the jet have to be taken into account: -E/p acceptance: it was shown in Section 3.3 that the probability to find E/p = 0 is strongly correlated to the amount of upstream material and hence a good measure of the E/p acceptance.A fully correlated (in p and η) 28% uncertainty is derived from the maximal observed difference between data and MC simulation in this probability.-Topological clustering effect: given an amount of energy released in the calorimeter, the energy collected in the topological clusters may differ from the energy released in the calorimeter, depending on how isolated the energy deposit is.This can introduce a bias in the particle response in a jet with respect to the response measured for isolated hadrons.A conservative systematic uncertainty has been computed by comparing the results of the E/p measurement obtained using clusters as described in Section 3 to those obtained by repeating the measurement using all calorimeter cells in a cone of size ∆R < 0.2.The double ratio of data/MC cluster response to data/MC cell response is used to estimate the uncertainty.The relevant result for the central calorimeter region is presented in Fig. 15, showing discrep-ancies of ∼ 5% at low p that disappear within the statistical uncertainties for p 10 GeV.-Out-of-cone ∆R > 0.2 energy deposits: statistical consistent results were found for the E/p data to Monte Carlo ratio for calorimeter energy measured in ∆R < 0.3.Hence no additional uncertainty is assumed.

Additional uncertainty in the jet energy scale
The E/p measurements only cover the response of charged hadronic particles with momenta less than ∼ 20 GeV.However, depending on the jet momentum, on average between 35% and 90% of the energy in jets is carried by particles that are not measured in situ using the isolated track analysis (mostly photons from π 0 decays, neutral hadrons and high momentum charged hadrons).Hence, the uncertainty on the energy response to these particles is needed in order to obtain the total calorimeter uncertainty on the jet energy scale.

High momentum charged particles
In 2004, an ATLAS Combined Test Beam (CTB) program was carried out at CERN.A "slice" of the ATLAS detector composed of the final versions of all sub-detectors in the barrel region was exposed to Super Proton Synchrotron (SPS) test beams.The layout of the sub-detectors was designed to be as close to that of ATLAS as possible.The setup was used to measure the combined calorimeter response to single charged pions of energies between 20 and 350 GeV for pseudorapidity values of 0.20, 0.25, 0.35, 0.45, 0.55 and 0.65 [13,18,19].From the measurements in Ref. [13], the ratio of data to MC simulation predictions is used to supplement the E/p measurements in Section 3 with a larger energy range.However, since these measurements are not made with the same detector, additional systematic uncertainties from the test beam have to be taken into account [13]: a fully correlated (in p and η) 0.7% total energy scale uncertainty for the energies in the LAr calorimeters; -a fully correlated 0.5% total energy scale uncertainty for the energies in the TileCal; -an additional 0.4% uncertainty on the LAr uniformity for all measured energies at the same η points; -an additional 1.5% uncertainty on the TileCal uniformity for all measured energies at the same η points; and a 1% uncertainty from the material in front of the calorimeter.
The ATLAS Collaboration For single particle momenta above 400 GeV no direct measurements in a test beam or in situ exist.Therefore an additional uncertainty of 10% is added in quadrature to that of the 350 GeV measurement uncertainty in order to cover possible effects from calorimeter nonlinearities at high energy densities and longitudinal leakage [20].

Absolute calorimeter energy scale
The absolute electromagnetic energy scale in ATLAS has been established using Z → ee decays for the electromagnetic LAr calorimeters and using the energy loss of minimum ionising muons in the TileCal.For the bulk of the electromagnetic LAr barrel calorimeter, the uncertainty on the cell energy measurement is 1.5%, and for the LAr presampler the uncertainty is 5% [21].For the TileCal, the scale uncertainty is 3% [22].This uncertainty does not affect charged particles with E/p measured in situ, but needs to be considered for all other particles contributing to jets.

Baryons and neutral hadronic particles
Test beam measurements of protons [13,18,19,20] and E/p measurements for identified pions and protons (Section 4) have shown that the agreement between data and MC simulation for protons is similar to the data to MC agreement for charged pions.The most significant data to MC difference in the energy response is 25% for low energy anti-protons.However, on average, these low energy anti-protons contribute no more than 0.5% to the jet energy.Hence no additional uncertainty for charged baryons is assumed.
In addition, the charged particle composition of the tracks (mainly π ± , p ± , K ± ) used for the inclusive E/p measurements in minimum bias events and the charged particle composition in jets is found to be sufficiently similar to cause no additional systematic uncertainty.Thus, even if some identified particle measurements show differences (as for low energy anti-protons), the inclusive E/p measures, on average, the charged particle response that is needed for an application to jets.The related uncertainty is proportional to the MC simulation modelling of the difference in the charged track composition between minimum bias events and jets and found to be negligible.
No test beam measurements for neutral hadronic particles have been carried out.Moreover, the Geant4 models have large uncertainties.On average 10-12% of the jet energy is carried by neutral hadrons, mostly K S , K L and neutrons.Most of the K S decay to pions before they reach the calorimeter.Hence the E/p and CTB measurements can be used for K S .
For neutrons and anti-neutrons Geant4 studies (see Fig. 16a and 16b) comparing alternative Geant4 hadronic physics models to the ATLAS-default hadronic physics model QGSP BERT [8] show that the (anti-)neutron to (anti-)proton response ratio is determined at the 10% level for particle momenta below 3 GeV and at the 5% level for higher momenta.Hence the (anti-)neutron response can be related to the sufficiently well simulated (anti-)proton response with these additional uncertainties.
Few measurements are available for kaon interaction cross sections in materials.While the response of charged kaons is covered by the inclusive E/p measurements, the K L response has to rely on Monte Carlo simulation predictions.Geant4 studies comparing alternative Geant4 hadronic physics models show uncertainties of ∼10% with respect to the ATLAS-default hadronic physics model QGSP BERT (see Fig. 16c).However, because of the limited availability of measurements, a more conservative uncertainty of 20% on the K L calorimeter response is added in quadrature to the uncertainty for charged particles [23].

Jet energy scale uncertainty estimation
This section combines the single particle uncertainties discussed in Sections 3, 5.1 and 5.2 into an expected shift and uncertainty on the calorimeter jet energy response with respect to the MC simulation prediction.Because of the CTB measurements used for high momentum particles, the estimation is limited to the pseudorapidity region |η| < 0.8.
Figure 17 shows the individual contributions to the uncertainty on the jet energy scale for |η| < 0. The E/p response and CTB single particle measurements show a shift of less than 1% between data and Monte Carlo simulation, which propagates into an expected shift of ∼ 0.5-1% in the JES (Fig. 17a and 17c).Due to the high precision of the measurements, the uncertainty is small.The visible drop in the response due to CTB measurements is caused by a systematic difference between the CTB data and Monte Carlo simulation for high momenta.
The effects of the E/p acceptance and of the clustering are small (Fig. 17b and Fig. 17e).The influence of the absolute electromagnetic energy scale uncertainty is small for low p T jets, where most of the energy is carried by particles with E/p measured in situ.However, high p T jets have only a small fraction of energy in the particle momentum range of the E/p analysis and are hence fully affected by this uncertainty (Fig. 17d).
The uncertainty on the calorimeter response to neutral hadrons was estimated in Section 5.2.3, resulting in a ∼ 1% contribution to the total JES uncertainty (Fig. 17f).
These individual contributions to the uncertainty are summarised in Fig. 18 together with the total expected shift and uncertainty.No single component is dominant and depending on the jet momentum several components contribute at approximately the same level to the total uncertainty.
Finally, the total calorimeter uncertainty on the jet energy scale is shown in Fig. 19 for anti-k t jets with R = 0.4 and R = 0.6 in the pseudorapidity range |η| < 0.8.For both jet sizes the maximum expected shift in the jet energy scale is ∼ 1% with an uncertainty of 1-3%.The envelope of the shift and uncertainty on the calorimeter JES is taken as the contribution to the total JES uncertainty discussed in Ref. [1].

Jet energy scale uncertainty correlations
The use of pseudo-experiments for the determination of the JES uncertainty allows a direct extraction of the correlation of uncertainties between different jet momenta, pseudorapidities or algorithms by correlating fluctuations of different quantities within each pseudoexperiment.Figure 20 shows the correlation of the JES uncertainty between different η and p bins.As expected, the correlation between neighboring bins in |η| and p T is almost 100%, while widely separated bins show only a ∼ 30% correlation.This remaining ∼ 30% correlation is mostly caused by the calorimeter energy scale and neutral hadron uncertainty, which contribute identically to all jets.
One specific use for the correlation of the uncertainty between jets is to compare the relative calorimeter response between jets selected from two different categories (different reconstruction algorithm, different originating parton, etc.).For each jet category, the procedure of Section 5.3 is repeated and the expected shift and uncertainty on the JES with respect to the MC simulation is derived 6 .For the determination of the uncertainty on the relative calorimeter response, the ratio of two jet categories is taken, cancelling out all common shifts and uncertainties.Fig. 16: Ratio of the average calorimeter response for alternative Geant4 hadronic physics models FTF BIC, FTFP BERT and QGSP BERT CHIPS [8] to the ATLAS default hadronic physics model QGSP BERT as function of the particle kinetic energy E kin in the range between 100 MeV and 50 GeV for (a) the ratio of neutrons n to protons p, (b) the ratio of anti-neutrons n to anti-protons p and (c) neutral kaons K L .
In Fig. 21a the relative calorimeter uncertainty on the JES between anti-k t jets reconstructed with R = 0.4 and R = 0.6 is shown.
Figure 21b shows the relative calorimeter uncertainty for b-tagged jets in t t events with respect to jets in inclusive di-jet events.Identified b-jets are associated with a displaced secondary vertex reconstructed by the SV0 algorithm [24], with a weight greater than 5.72.This weight gives an average b-tag efficiency of 50% for b-jets from t t events.Finally, Fig. 21c shows the relative uncertainty for jets initiated mostly by quarks with respect to jets in inclusive di-jet events.Jets initiated mostly by quarks are selected in γ+jet event simulations using the se-lection criteria applied in the method of the direct p T balance between γ and jets [1].
Because of the small differences in particle composition between these different jet samples, the relative uncertainty on the pure calorimeter response is found to be below 0.5% for low p T and negligible for high p T .

Conclusions
The average calorimeter response to isolated hadrons with respect to the track momentum E/p has been measured in minimum bias events at √ s = 900 GeV and √ s = 7 TeV.A background from neutral hadrons of 4-8% is subtracted from data and Monte Carlo with   Fig. 19: Expected shift and total uncertainty on the relative calorimeter jet response with respect to the MC simulation for jets reconstructed with the anti-k t jet algorithm (R = 0.4 and R = 0.6) in the range |η| < 0.3 and 0.3 ≤ |η| < 0.8 as function of the jet transverse momentum.The x-axis is the jet transverse momentum calibrated from the EM scale to the hadronic scale using an MC-based JES calibration factor [1].
a systematic uncertainty of below 1%.After the background has been removed, the agreement between data and Monte Carlo simulation is within 2% for particles with momenta up to 10 GeV and is around 5% for momenta in the 10-30 GeV range, where the statistical uncertainty dominates the total uncertainty.
The calorimeter response of identified single charged hadrons has been measured using short-lived particle decays.A good agreement between data measured at √ s = 7 TeV and the Monte Carlo simulation is found for charged pions and protons.However, a disagreement of up to 10% is found between the Monte Carlo simulation and data for the difference of responses to low momentum anti-pions and anti-protons ( E/p π − − E/p p).This difference is attributed to the poor modelling of the anti-proton response in the Monte Carlo simulation.
The ATLAS calorimeter jet energy scale uncertainty has been determined for the well understood central detector region by propagating the energy response uncertainty of all particles contributing to a jet.For charged hadron momenta below 20 GeV, the single charged hadron response has been used, while for higher momenta, the response measured in the ATLAS combined test beam has been included.For the response to neutral pions (π 0 → γγ), the uncertainty on the electromagnetic calorimeter energy scale is dominant, while for all other neutral hadrons an additional uncertainty due to the limited  An uncertainty of 1-3% on the response to jets in the calorimeter is determined for jets in the central region of the detector with |η| < 0.8 and a transverse momentum between 15 GeV and 2.5 TeV.An analysis of the correlation in the jet energy scale uncertainty yields an almost complete correlation for jets close in |η| or p T , and a correlation of ∼ 30% for jets with a large transverse momentum separation.The relative jet energy scale uncertainty between anti-k t jets of size R = 0.4 and R = 0.6, between b-tagged jets and inclusive jets in di-jet events and between mostly quark-initiated jets and inclusive jets in the same event sample was found to be small.

A Validation of the background subtraction
The background subtraction procedure relies on the assumption that the energy density deposited by other particles in the cone around a late showering track is constant as a function of the distance with respect to the track impact point on the calorimeter.Although this assumption is valid for low track momenta, the validity of the assumption for high track momenta is questionable, as jet-like structures are expected to emerge.
In order to address this issue, a different background subtraction procedure has been developed.The background is still measured in events where the energy in the EM calorimeter associated to the track is compatible with that of a late showering hadron (MIP-tagged track, in the following).In this case, however, the energy density of the background in small annuli between a cone of radius 0.1 and 0.2 is used to linearly extrapolate the background inside the cone of size 0.1.The differential energy density ρ i = d 2 E/dηdφ for an annulus of inner radius R i = i × 0.025 and outer radius R i+1 = (i + 1) × 0.025 is defined as The ATLAS Collaboration (a) between jets reconstructed with the antik t jet algorithm with R = 0.4 and jets reconstructed with R = 0.6, (b) between b-tagged jets in t t events and jets in inclusive di-jet events, and (c) between mostly quark-initiated jets in γ+jet events and jets in inclusive di-jet events.The x-axis is the jet transverse momentum calibrated from the EM scale to the hadronic scale using an MC-based JES calibration factor [1]. where The dependence of ρ i on the distance from the track impact point on the calorimeter has been studied by making use of the true energy deposits as predicted by the Geant4 simulation.The true energy deposited by the track and by the neutral background for late showering tracks, and by the neutral background in all selected events is shown in Fig. 22 for two different track momentum bins.The energy density of the MIP-tagged tracks is very narrow in both momentum bins, as expected.The background energy density is constant and independent of the distance to the track impact point at low track momenta.This is consistent with the assumption of constant energy density made to evaluate the background with the baseline method.As the track momentum increases, the energy density close to the track impact point is higher, perhaps indicating a jet-like structure around the track.In both track momentum bins, the background energy density in the annulus 0.1 ≤ R < 0.2 is the same for MIP-tagged and all tracks.
The same features can be observed in Fig. 23.This time, the reconstructed energy density is shown for data and MC simulation (the points corresponding to the true energy density associated to the MIP-tagged track are kept for reference).A comparison of Fig. 22 and Fig. 23 reveals that the main features observed for the true energy densities still hold for reconstructed energy densities: the density in the halo of the cone is a good measurement of the true background density, while the core of the cone mainly contains energy from the MIP track.
The alternative estimate of the background is obtained by fitting the measured energy density with a linear function between R i = 0.1 and R i = 0.25 and integrating the obtained function between R i = 0 and R i = 0.25: The obtained background estimation gives a total of 10% more background at high track momenta with respect to the baseline background subtraction procedure, consistently for data and MC simulation.The reason for the small difference is that the two background estimates differ only in the core of the cone, which is a region of small area.
Since the background itself is a 10% correction to the E/p measurement, the total effect introduced by the linear parameterisation of the energy density is of the order of 1% on E/p .Since data and MC points are affected in a similar way, the effect on the agreement between MC and data is negligible.

Fig. 1 :Fig. 2 :
Fig.1:(a) E/p distribution for isolated tracks with an impact point in the region |η| < 0.6 and with a momentum in the range 1.2 ≤ p < 1.8 GeV.(b) E/p distribution for tracks with impact points in the region 1.9 ≤ |η| < 2.3 and with momenta in the range 2.8 ≤ p < 3.6 GeV.

Fig. 3 :
Fig.3: Sketch of the background subtraction method, selecting isolated charged hadrons (ICH) that pass the EM calorimeter as minimum ionising particles (MIP) and shower in the hadronic calorimeter.

Fig. 4 :Fig. 5 :
Fig. 4: E/p BG as a function of the track momentum at √ s = 900 GeV for (a) |η| < 0.6 and (b) 0.6 ≤ |η| < 1.1.The black markers represent the background estimated from collision data, while the green rectangles represent the MC prediction, with the vertical width showing its statistical uncertainty.The lower panes show the ratio of the MC prediction to collision data.The MC prediction of the background energy deposit is obtained using the same procedure as applied to the data.

Fig. 6 :
Fig. 6: E/p at √ s = 900 GeV as a function of the track momentum for (a) |η| < 0.6 and (b) 0.6 ≤ |η| < 1.1.The black markers represent the collision data, while the green rectangles represent the MC prediction, with the vertical width showing its statistical uncertainty.The lower panes show the ratio of the MC simulation prediction to collision data.The grey band indicates the size of the systematic uncertainty on the measurement.The dotted lines are placed at ±5% of unity and at unity.

Fig. 7 :Fig. 8 :
Fig. 7: E/p at √ s = 7 TeV as a function of the track momentum for (a) |η| < 0.6 and (b) 0.6 ≤ |η| < 1.1.The black markers represent the collision data, while the green rectangles represent the MC prediction, with the vertical width showing its statistical uncertainty.The lower panes show the ratio of the MC simulation prediction to collision data.The grey band indicates the size of the systematic uncertainty on the measurement.The dotted lines are placed at ±5% of unity and at unity.

Fig. 9 :
Fig. 9: The uncorrected E/p distribution for (a) single π − and p and (b) π + and p in a single bin of available energy and pseudorapidity, 2.2 ≤ E a < 2.8 GeV and |η| < 0.6.

Fig. 15 :
Fig. 15: Ratio of the E/p measurement obtained with topological clusters to that obtained using all calorimeter cells in (a) the central (|η| < 0.6) and (b) forward (0.6 ≤ |η| < 1.1) regions.The inset shows the ratio of the MC prediction to data.

Fig. 17 :Fig. 18 :
Fig.17: Expected shift (black markers) and uncertainty (error bars) on the relative calorimeter jet response with respect to the MC simulation for jets reconstructed with the anti-k t jet algorithm (R = 0.6) in the range |η| < 0.3 as function of the jet transverse momentum from the uncertainty on (a) the E/p response, (b) the uncertainty on the E/p acceptance, (c) the uncertainty from the CTB response, (d) the uncertainty from the absolute energy scale, (e) the uncertainty from the clustering effect, and (f) the uncertainty from neutral hadrons.The x-axis is the jet transverse momentum calibrated from the EM scale to the hadronic scale using an MC-based JES calibration factor[1].

1 ATLASFig. 20 :
Fig. 20: Correlation coefficient of the total uncertainty on the calorimeter jet response for jets reconstructed with the anti-k t jet algorithm (R = 0.6) in the range |η| < 0.8 and jet transverse momenta between 15 GeV≤ p T <2.5 TeV.The bins on the x-axis are identical to the bins on the y-axis.A dark area indicates highly correlated uncertainties, a light area almost uncorrelated uncertainties.

Fig. 21 :
Fig.21: Expected shift (black dots) and uncertainty (error bars) on the relative calorimeter jet response as function of the jet transverse momentum for jets in the range 0.3 ≤ |η| < 0.8: (a) between jets reconstructed with the antik t jet algorithm with R = 0.4 and jets reconstructed with R = 0.6, (b) between b-tagged jets in t t events and jets in inclusive di-jet events, and (c) between mostly quark-initiated jets in γ+jet events and jets in inclusive di-jet events.The x-axis is the jet transverse momentum calibrated from the EM scale to the hadronic scale using an MC-based JES calibration factor[1].

Fig. 22 :Fig. 23 :
Fig. 22: Energy density, ρ, as a function of the radial distance R in η × φ space from the track impact point for tracks with (a) 1.5 ≤ p < 1.8 GeV and (b) 3.6 ≤ p < 4.6 GeV.The blue rectangles with diagonal hatching show the true energy density associated with MIP tracks.The yellow rectangles with horizontal hatching (red rectangles with vertical hatching) show the true energy density in the EM calorimeter associated to background particles for MIP-tagged (all) tracks.
Fig. 12: The difference in E/p between π − from K S candidates and p from Λ candidates in tracks with (a) |η| < 0.6 and (b) 0.6 ≤ |η| < 1.1.The responses to the anti-protons are corrected for charged particle backgrounds (see text).