# Jet energy measurement and its systematic uncertainty in proton-proton collisions at $\sqrt{s} = 7$ TeV with the ATLAS detector

G. Aad, S. Albrand, J. Brown, Q. Buat, B. Clement, J. Collot, S.

Crépé-Renaudin, B. Dechenaux, T. Delemontex, P.A. Delsart, et al.

► To cite this version:

G. Aad, S. Albrand, J. Brown, Q. Buat, B. Clement, et al., Jet energy measurement and its systematic uncertainty in proton-proton collisions at $\sqrt{s} = 7$ TeV with the ATLAS detector. European Physical Journal C: Particles and Fields, 2015, 75, pp.17. 10.1140/epjc/s10052-014-3190-y. in2p3-01003959

HAL Id: in2p3-01003959

https://hal.in2p3.fr/in2p3-01003959

Submitted on 4 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Jet energy measurement and its systematic uncertainty in proton–proton collisions at $\sqrt{s} = 7$ TeV with the ATLAS detector

The ATLAS Collaboration

## Abstract

The jet energy scale (JES) and its systematic uncertainty are determined for jets measured with the ATLAS detector using proton–proton collision data with a centre-of-mass energy of $\sqrt{s} = 7$ TeV corresponding to an integrated luminosity of 4.7 fb$^{-1}$. Jets are reconstructed from energy deposits forming topological clusters of calorimeter cells using the anti-$k_t$ algorithm with distance parameters $R = 0.4$ or $R = 0.6$, and are calibrated using MC simulations. A residual JES correction is applied to account for differences between data and MC simulations. This correction and its systematic uncertainty are estimated using a combination of in situ techniques exploiting the transverse momentum balance between a jet and a reference object such as a photon or a Z boson, for $20 \le p_T^{\text{jet}} < 1000$ GeV and pseudorapidities $|\eta| < 4.5$. The effect of multiple proton–proton interactions is corrected for, and an uncertainty is evaluated using in situ techniques. The smallest JES uncertainty of less than 1% is found in the central calorimeter region ($|\eta| < 1.2$) for jets with $55 \le p_T^{\text{jet}} < 500$ GeV. For central jets at lower $p_T$, the uncertainty is about 3%. A consistent JES estimate is found using measurements of the calorimeter response of single hadrons in proton–proton collisions and test-beam data, which also provide the estimate for $p_T^{\text{jet}} > 1$ TeV. The calibration of forward jets is derived from dijet $p_T$ balance measurements. The resulting uncertainty reaches its largest value of 6% for low-$p_T$ jets at $|\eta| = 4.5$. Additional JES uncertainties due to specific event topologies, such as close-by jets or selections of event samples with an enhanced content of jets originating from light quarks or gluons, are also discussed. The magnitude of these uncertainties depends on the event sample used in a given physics analysis, but typically amounts to 0.5% to 3%.
# Contents

1 Introduction 2

2 The ATLAS detector 2

2.1 Detector description 2

2.2 Calorimeter pile-up sensitivity 3

3 Monte Carlo simulation of jets in the ATLAS detector 3

3.1 Inclusive jet Monte Carlo simulation samples 3

3.2 Z-jet and γ-jet Monte Carlo simulation samples 3

3.3 Top-quark pair Monte Carlo simulation samples 4

3.4 Minimum bias samples 5

3.5 Detector simulation 5

4 Dataset 5

5 Jet reconstruction and calibration with the ATLAS detector 5

5.1 Topological clusters in the calorimeter 5

5.2 Jet reconstruction and calibration 7

5.3 Jet quality selection 7

5.4 Track jets 9

5.5 Truth jets 9

5.6 Jet kinematics and directions 10

6 Jet energy correction for pile-up interactions 10

6.1 Pile-up correction method 10

6.2 Principal pile-up correction strategy 10

6.3 Derivation of pile-up correction parameters 12

6.4 Pile-up validation with in situ techniques and effect of out-of-time pile-up in different calorimeter regions 14

7 In situ transverse momentum balance techniques 14

7.1 Relative in situ calibration between the central and forward rapidity regions 14

7.2 In situ calibration methods for the central rapidity region 15

8 Relative forward-jet calibration using dijet events 15

8.1 Intercalibration using events with dijet topologies 15

8.2 Event selection for dijet analysis 16

8.3 Dijet balance results 17

8.4 Systematic uncertainty 19

8.5 Summary of the η-intercalibration and its uncertainties 21

9 Jet energy calibration using Z-jet events 21

9.1 Description of the $p_T$ balance method 21

9.2 Selection of Z-jet events 21

9.3 Measurement of the $p_T$ balance 23

9.4 Measuring out-of-cone radiation and underlying event contributions 23

9.5 Systematic uncertainties 24

9.6 Summary of the Z-jet analysis 27

10 Jet energy calibration using γ-jet events 27

10.1 In situ jet calibration techniques 27

10.2 Event selection of γ-jet events 27

10.3 Jet response measurement 31

10.4 Systematic uncertainties of photon-jet balance 31

10.5 Summary of the γ-jet analysis 34

11 High-$p_T$ jet energy calibration using multijet events 34

11.1 Multijet balance technique and uncertainty propagation 34

11.2 Selection of multijet events 37

11.3 Multijet balance measurement 37

11.4 Systematic uncertainties on the multijet balance 37

11.5 Summary of multijet analysis 42

12 Forward-jet energy measurement validation using Z-jet and γ-jet data 42

13 Jet energy calibration and uncertainty combination 44

13.1 Overview of the combined JES calibration procedure 44

13.2 Combination technique 44

13.3 Uncertainty sources of the in situ calibration techniques 46

13.4 Combination results 48

13.5 Comparison of the γ-jet calibration methods 49

13.6 Simplified description of the correlations 49

13.7 Jet energy scale correlation scenarios 50

13.8 Alternative reduced configurations 51

14 Comparison to jet energy scale uncertainty from single-hadron response measurements 52

15 Jet energy scale uncertainty from the W boson mass constraint 53

15.1 Event samples 53

15.2 Reconstruction of the W boson 54

15.3 Extraction of the relative light jet scale 54

15.4 Systematic uncertainties 55

15.5 Results 55

16 Systematic uncertainties on corrections for pile-up interactions 55

16.1 Event and object selection 55

16.2 Derivation of the systematic uncertainty 55

16.3 Summary on pile-up interaction corrections 61

17 Close-by jet effects on jet energy scale 62

17.1 Samples and event selection 62

17.2 Non-isolated jet energy scale uncertainty 62

18 Jet response difference for quark and gluon induced jets and associated uncertainty 63

18.1 Event selection 63

18.2 Calorimeter response to quark and gluon induced jets 64

18.3 Discrimination of light-quark and gluon induced jets 66

18.4 Summary of the jet flavour dependence analysis 66

19 Jets with heavy-flavour content 66

19.1 Jet selection and response definition 67

19.2 Track selection 68

19.3 Event selection 68

19.4 MC-based systematic uncertainties on the calorimeter b-jet energy scale 69

19.5 Calorimeter jet energy measurement validation using tracks 69

19.6 Systematic uncertainties 69

19.7 Results 73

19.8 Semileptonic correction and associated uncertainties 73

19.9 Semileptonic neutrino energy validation using dijet balance 74

19.10 Conclusions on heavy-flavour jets 74

20 Jet response in problematic calorimeter regions 75

20.1 Correction algorithms for non-operating calorimeter modules 75

20.2 Performance of the bad calorimeter region corrections 76

21 Summary of the total jet energy scale systematic uncertainty 76

22 Conclusions 84

23 Acknowledgements 84
# 1 Introduction

Jets are the dominant feature of high-energy, hard proton-proton interactions at the Large Hadron Collider (LHC) at CERN. They are key ingredients of many physics measurements and for searches for new phenomena. In this paper, jets are observed as groups of topologically related energy deposits in the ATLAS calorimeters, associated with tracks of charged particles as measured in the inner tacking detector. They are reconstructed with the anti-$k_t$ jet algorithm [1] and are calibrated using Monte Carlo (MC) simulation.

A first estimate of the jet energy scale (JES) uncertainty of about 5% to 9% depending on the jet transverse momentum ($p_T$), described in Ref. [2], is based on information available before the first proton-proton collisions at the LHC, and initial proton-proton collision data taken in 2010. A reduced uncertainty of about 2.5% in the central calorimeter region over a wide $p_T$ range of $60 \lesssim p_T < 800$ GeV was achieved after applying the increased knowledge of the detector performance obtained during the analysis of this first year of ATLAS data taking [3]. This estimation used single-hadron calorimeter response measurements, systematic variations of MC simulation configurations, and in situ techniques, where the jet transverse momentum is compared to the $p_T$ of a reference object. These measurements were performed using the 2010 dataset, corresponding to an integrated luminosity of $38 \text{ pb}^{-1}$ [4].

During the year 2011 the ATLAS detector [5] collected proton-proton collision data at a centre-of-mass energy of $\sqrt{s} = 7$ TeV, corresponding to an integrated luminosity of about $4.7 \text{ fb}^{-1}$.
The larger dataset makes it possible to further improve the precision of the jet energy measurement, and also to apply a correction derived from detailed comparisons of data and MC simulation using in situ techniques. This document presents the results of such an improved calibration of the jet energy measurement and the determination of the uncertainties using the 2011 dataset.

The energy measurement of jets produced in proton-proton and electron-proton collisions is also discussed by other experiments [6–17].

The outline of the paper is as follows. Section 2 describes the ATLAS detector. The Monte Carlo simulation framework is presented in Sect. 3, and the used dataset is described in Sect. 4. Section 5 summarises the jet reconstruction and calibration strategy. The correction method for the effect of additional proton-proton interactions is discussed in Sect. 6. Section 7 provides an overview of the techniques based on $p_T$ balance that are described in detail in Sects. 8 to 11. First the intercalibration between the central and the forward detector using events with two high-$p_T$ jets is presented in Sect. 8. Then, in situ techniques to assess differences of the jet energy measurement between data and Monte Carlo simulation exploiting the $p_T$ balance between a jet and a well-measured reference object are detailed. The reference objects are Z bosons in Sect. 9, photons in Sect. 10, and a system of low-$p_T$ jets in Sect. 11. The validation of the forward-jet energy measurements with $p_T$ balance methods using Z-jet and $\gamma$-jet events follows in Sect. 12. The strategy on how to extract a final jet calibration out of the combination of in situ techniques, and the evaluation strategies for determining the corresponding systematic uncertainties, are discussed in Sect. 13. The same section also shows the final result

of the jet calibration, including its systematic uncertainty, from the combination of the in situ techniques.

Section 14 compares the JES uncertainty as derived from the single-hadron calorimeter response measurements to that obtained from the in situ method based on $p_T$ balance discussed in the preceding sections. Comparisons to JES uncertainties using the W boson mass constraint in final states with hadronically decaying W bosons are presented in Sect. 15.

Additional contributions to the systematic uncertainties of the jet measurement in ATLAS are presented in Sects. 16 to 18, where the correction for the effect of additional proton-proton interactions in the event, the presence of other close-by jets, and the response dependence on the jet fragmentation (jet flavour) are discussed. The uncertainties for explicitly tagged jets with heavy-flavour content are outlined in Sect. 19. A brief discussion of the correction of the calorimeter energy in regions with hardware failures and the associated uncertainty on the jet energy measurement is presented in Sect. 20.

A summary of the total jet energy scale uncertainty is given in Sect. 21. Conclusions follow in Sect. 22. A comparison of the systematic uncertainties of the JES in ATLAS with previous calibrations is presented in Appendix A.

# 2 The ATLAS detector

## 2.1 Detector description

The ATLAS detector consists of a tracking system (Inner Detector, or ID in the following), sampling electromagnetic and hadronic calorimeters and muon chambers. A detailed description of the ATLAS experiment can be found in Ref. [5].

The Inner Detector has complete azimuthal coverage and spans the pseudorapidity<fnref n="1" /> region $|\eta| < 2.5$. It consists of layers of silicon pixel detectors, silicon microstrip detectors and transition radiation tracking detectors, all of which are immersed in a solenoid magnet that provides a uniform magnetic field of 2 T.

Jets are reconstructed using the ATLAS calorimeters, whose granularity and material varies as a function of $\eta$. The electromagnetic calorimetry (EM) is provided by high-granularity liquid-argon sampling calorimeters (LAr), using lead as an absorber. It is divided into one barrel ($|\eta| < 1.475$) and two end-cap ($1.375 < |\eta| < 3.2$) regions. The hadronic calorimetry is divided into three distinct sections. The most central contains the central barrel region ($|\eta| < 0.8$) and two extended barrel regions ($0.8 < |\eta| < 1.7$). These regions are instrumented with scintillator-tile/steel hadronic calorimeters (Tile). Each barrel region consists of 64 modules with individual $\phi$ coverages of $\sim 0.1$ rad. The two hadronic end-cap calorimeters (HEC; $1.5 < |\eta| < 3.2$) feature liquid-argon/copper calorimeter modules. The two forward calorimeters (FCal; $3.1 < |\eta| < 4.9$) are

<footnote n="1">¹ ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upward. Cylindrical coordinates $(r, \phi)$ are used in the transverse plane, $\phi$ being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of the polar angle $\theta$ as $\eta = -\ln \tan(\theta/2)$.</footnote>
instrumented with liquid-argon/copper and liquid-argon/tungsten modules to provide electromagnetic and hadronic energy measurements, respectively.

The muon spectrometer surrounds the ATLAS calorimeter. A system of three large air-core toroids, a barrel and two end-caps, generates a magnetic field in the pseudorapidity range of $|η| < 2.7$. The muon spectrometer measures muon tracks with three layers of precision tracking chambers and is instrumented with separate trigger chambers.

The trigger system for the ATLAS detector consists of a hardware-based Level 1 (L1) and a software-based High Level Trigger (HLT) [18]. At L1, jets are first built from coarse-granularity calorimeter towers using a sliding window algorithm, and then subjected to early trigger decisions. This is refined using jets reconstructed from calorimeter cells in the HLT, with algorithms similar to the ones applied offline.

## 2.2 Calorimeter pile-up sensitivity

One important feature for the understanding of the contribution from additional proton-proton interactions (pile-up) to the signal in the 2011 dataset is the sensitivity of the ATLAS liquid argon calorimeters to the bunch crossing history. In any LAr calorimeter cell, the reconstructed energy is sensitive to the proton-proton interactions occurring in approximately 12 (2011 data, 24 at LHC design conditions) preceding and one immediately following bunch crossings (*out-of-time pile-up*), in addition to pile-up interactions in the current bunch crossing (*in-time pile-up*). This is due to the relatively long charge collection time in these calorimeters (typically 400–600 ns), as compared to the bunch crossing intervals at the LHC (design 25 ns and actually 50 ns in 2011 data). To reduce this sensitivity, a fast, bipolar shaped signal<fnref n="2" /> is used with net zero integral over time.

The signal shapes in the liquid argon calorimeters are optimised for this purpose, leading to cancellation on average of in-time and out-of-time pile-up in any given calorimeter cell. By design of the shaping amplifier, the most efficient suppression is achieved for 25 ns bunch spacing in the LHC beams. It is fully effective in the limit where, for each bunch crossing, about the same amount of energy is deposited in each calorimeter cell.

The 2011 beam conditions, with 50 ns bunch spacing and a relatively low cell occupancy from the achieved instantaneous luminosities, do not allow for full pile-up suppression by signal shaping, in particular in the central calorimeter region. Pile-up suppression is further limited by large fluctuations in the number of additional interactions from bunch crossing to bunch crossing, and in the energy flow patterns of the individual collisions in the time window of sensitivity of approximately 600 ns. Consequently, the shaped signal extracted by digital filtering shows a principal sensitivity to in-time and out-of-time pile-up, in particular in terms of a residual non-zero cell-signal baseline. This baseline can lead to relevant signal offsets once the noise suppression, an important part of the calorimeter signal extraction strategy presented in Sect. 5, is applied.

Corrections mitigating the effect of these signal offsets on the reconstructed jet energy are discussed in the context of the pile-up suppression strategy in Sect. 6.1. All details of the ATLAS liquid argon calorimeter readout and signal processing can be found in Ref. [19].

The Tile calorimeter shows very little sensitivity to pile-up since most of the associated (soft particle) energy flow is absorbed in the LAr calorimeters in front of it. Moreover, out-of-time pile-up is suppressed by a short shaping time with sensitivity to only about 3 bunch crossings [20].

# 3 Monte Carlo simulation of jets in the ATLAS detector

The energy and direction of particles produced in proton-proton collisions are simulated using various MC event generators. An overview of these generators for LHC physics can be found in Ref. [21]. The samples using different event generators and theoretical models are described below. All samples are produced at $\sqrt{s} = 7$ TeV.

## 3.1 Inclusive jet Monte Carlo simulation samples

1. PYTHIA (version 6.425) [22] is used for the generation of the baseline simulation event samples. It models the hard sub-process in the final states of the generated proton-proton collisions using a $2 \to 2$ matrix element at leading order in the strong coupling $\alpha_s$. Additional radiation is modelled in the leading logarithmic (LL) approximation by $p_T$-ordered parton showers [23].

Multiple parton interactions (MPI) [24], as well as fragmentation and hadronisation based on the Lund string model [25], are also generated. Relevant parameters for the modelling of the parton shower and multiple parton interactions in the underlying event (UE) are tuned to LHC data (ATLAS PYTHIA tune AUET2B [26] with the MRST LO** parton density function (PDF) [27]). Data from the LEP collider are included in this tune.

2. HERWIG++ [28] is used to generate samples for evaluating systematic uncertainties. This generator uses a $2 \to 2$ matrix element and angular-ordered parton showers in the LL approximation [29–31]. The cluster model [32] is employed for the hadronisation. The underlying event and soft inclusive interactions are described using a hard and soft MPI model [33]. The parton densities are provided by the MRST LO** PDF set.

3. MADGRAPH [34] with the CTEQ6L1 PDF set [35] is used to generate proton-proton collision samples with up to three outgoing partons from the matrix element and with MLM matching [36] applied in the parton shower, which is performed with PYTHIA using the AUET2B tune.

## 3.2 Z-jet and $\gamma$-jet Monte Carlo simulation samples

1. PYTHIA (version 6.425) is used to produce Z-jet events with the modified leading-order PDF set MRST LO**. The

<footnote n="2">² The shaped pulse has a duration exceeding the charge collection time.</footnote>
<figure><img src="image_2.png" /><figcaption>Fig. 1: The energy-equivalent cell noise in the ATLAS calorimeters on the electromagnetic (EM) scale as a function of the direction $|η|$ in the detector, for the 2010 configuration with (a) $\mu = 0$ and the 2011 configuration with (b) $\mu = 8$. The various colours indicate the noise in the pre-sampler (PS) and the up to three layers of the LAr EM calorimeter, the up to three layers of the Tile calorimeter, the four layers for the hadronic end-cap (HEC) calorimeter, and the three modules of the forward (FCal) calorimeter.</figcaption></figure>

simulation uses a $2 \to 1$ matrix element to model the hard sub-process, and, as for the inclusive jet simulation, $p_T$-ordered parton showers to model additional parton radiation in the LL approximation. In addition, weights are applied to the first branching of the shower, so as to bring agreement with the matrix-element rate in the hard emission region. The same tune and PDF is used as for the inclusive jet sample.

2. The ALPGEN generator (version 2.13) [37] is used to produce Z-jet events, interfaced to HERWIG (version 6.510) [31] for parton shower and fragmentation into particles, and to JIMMY (version 4.31) [38] to model UE contributions using the ATLAS AUET2 tune [39], here with the CTEQ6L1 [35] leading-order PDF set. ALPGEN is a leading-order matrix-element generator for hard multi-parton processes ($2 \to n$) in hadronic collisions. Parton showers are matched to the matrix element with the MLM matching scheme. The CTEQ6L1 PDF set is employed.

3. The baseline $\gamma$-jet sample is produced with PYTHIA (version 6.425). It generates non-diffractive events using a $2 \to 2$ matrix element at leading order in $\alpha_s$ to model the hard sub-process. Again, additional parton radiation is modelled by $p_T$-ordered parton showers in the LL approximation. The modelling of non-perturbative physics effects arising in MPI, fragmentation, and hadronisation is based on the ATLAS AUET2B MRST LO** tune.

4. An alternative $\gamma$-jet event sample is generated with HERWIG (version 6.510) and JIMMY using the ATLAS AUET2 tune and the MRST LO** PDF. It is used to evaluate the systematic uncertainty due to physics modelling.

5. The systematic uncertainty from jets which are misidentified as photons (fake photons) is studied with a dedicated

MC event sample. An inclusive jet sample is generated with PYTHIA (version 6.425) with the same parameter tuning and PDF set as the $\gamma$-jet sample. An additional filter is applied to the jets built from the stable generated particles to select events containing a narrow particle jet, which is more likely to pass photon identification criteria. The surviving events are passed through the same detector simulation software as the MC $\gamma$-jet sample.

### 3.3 Top-quark pair Monte Carlo simulation samples

Top pair ($t\bar{t}$) production samples are relevant for jet reconstruction performance studies, as they are a significant source of hadronically decaying W bosons and therefore important for light-quark jet response evaluations in a radiation environment very different from the inclusive jet and Z-jet/$\gamma$-jet samples discussed above. In addition, they provide jets from a heavy-flavour ($b$-quark) decay, the response to which can be studied in this final state as well.

The nominal $t\bar{t}$ event sample is generated using MC@NLO (version 4.01) [40], which implements a next-to-leading-order (NLO) matrix element for top-pair production. Correspondingly, the CT10 [41] NLO PDF set is used. This matrix-element generator is interfaced to parton showers from HERWIG (version 6.520) [42] and the underlying event modelled by JIMMY (version 4.31), with the CT10 PDF and the ATLAS AUET2 tune.

A number of systematic variation samples use alternative MC generators or different generator parameter sets. Additional $t\bar{t}$ samples are simulated using the POWHEG [43] generator interfaced with PYTHIA, as well as HERWIG and JIMMY.
POWHEG provides alternative implementations of the NLO matrix-element calculation and the interface to parton show- ers. These samples allow comparison of two different parton shower, hadronisation and fragmentation models. In addition, the particular implementations of the NLO matrix-element cal- culations in POWHEG and MC@NLO can be compared. Dif- ferences in the b -hadron decay tables between PYTHIA and HERWIG are also significant enough to provide a conservative uncertainty envelope on the effects of the decay model.

In addition, samples with more or less parton shower activ- ity are generated with the leading-order generator ACERMC [44] interfaced to PYTHIA with the MRST LO** PDF set. The- se are used to estimate the model dependence of the event se- lection. In these samples the initial state radiation (ISR) and the final state radiation (FSR) parameters are varied in value ranges not excluded by the current experimental data, as de- tailed in Refs. [45,46].

3.4 Minimum bias samples

Minimum bias events are generated using PYTHIA8 [47] with the 4C tune [48] and MRST LO** PDF set. These minimum bias events are used to form pile-up events, which are overlaid onto the hard-scatter events following a Poisson distribution around the average number 〈 µ 〉 of additional proton–proton collisions per bunch crossing measured in the experiment. The LHC bunch train structure with 36 proton bunches per train and 50 ns spacing between the bunches, is also modelled by organising the simulated collisions into four such trains. This allows the inclusion of out-of-time pile-up effects driven by the distance of the hard-scatter events from the beginning of the bunch train. The first ten bunch crossings in each LHC bunch train, approximately, are characterised by varying out- of-time pile-up contributions from the collision history, which is getting filled with an increasing number of bunch crossings with proton–proton interactions. For the remaining ≈ 26 bunch crossings in a train, the effect of the out-of-time pile-up con- tribution is stable, i.e. it does not vary with the bunch posi- tion within the bunch train, if the bunch-to-bunch intensity is constant. Bunch-to-bunch fluctuations in proton intensity at the LHC are not included in the simulation.

3.5 Detector simulation

The GEANT4 software toolkit [49] within the ATLAS simula- tion framework [50] propagates the stable particles<fnref n="3" /> produced by the event generators through the ATLAS detector and sim- ulates their interactions with the detector material. Hadronic showers are simulated with the QGSP_BERT model [51–59]. Compared to the simulation used in the context of the 2010 data analysis, a newer version of GEANT4 (version 9.4) is used and a more detailed description of the geometry of the LAr calori- meter absorber structure is available. These geometry changes introduce an increase in the calorimeter response to pions be- low 10 GeV of about 2%.

For the estimation of the systematic uncertainties arising from detector simulation, several samples are also produced with the ATLAS fast (parameterised) detector simulation ATL- FAST2 [60,61].

4 Dataset

The data used in this study were recorded by ATLAS between May and October 2011, with all ATLAS subdetectors opera- tional. The corresponding total integrated luminosity is about 4.7 fb − 1 of proton–proton collisions at a centre-of-mass energy of √$_{s}$ = 7 TeV.

As already indicated in Sect. 3.4, the LHC operated with bunch crossing intervals of 50 ns, and bunches organised in bunch trains. The average number of interactions per bunch crossing ( µ ) as estimated from the luminosity measurement is 3 ≤ µ ≤ 8 until Summer 2011, with an average for this period of 〈 µ 〉 ≈ 6. Between August 2011 and the end of the proton run, µ increased to about 5 ≤ µ ≤ 17, with an average 〈 µ 〉 ≈ 12. The average number of interactions for the whole 2011 dataset is 〈 µ 〉 = 8.

The specific trigger requirements and precision signal ob- ject selections applied to the data are analysis dependent. They are therefore discussed in the context of each analysis presented in this paper.

5 Jet reconstruction and calibration with the ATLAS detector

5.1 Topological clusters in the calorimeter

Clusters of energy deposits in the calorimeter (topo-clusters) are built from topologically connected calorimeter cells that contain a significant signal above noise, see Refs. [3, 62, 63] for details. The topo-cluster formation follows cell signal sig- nificance patterns in the ATLAS calorimeters. The signal sig- nificance is measured by the absolute ratio of the cell signal to the energy-equivalent noise in the cell. The signal-to-noise thresholds for the cluster formation are not changed with re- spect to the settings given in Ref. [3]. However, the noise in the calorimeter increased due to the presence of multiple proton- proton interactions, as discussed in Sect. 2.2, and required the adjustments explained below.

While in ATLAS operations prior to 2011 the cell noise was dominated by electronic noise, the short bunch crossing interval in 2011 LHC running added a noise component from bunch-to-bunch variations in the instantaneous luminosity and in the energy deposited in a given cell from previous colli- sions inside the window of sensitivity of the calorimeters. The cell noise thresholds steering the topo-cluster formation thus needed to be increased from those used in 2010 to accommo- date the corresponding fluctuations, which is done by raising the nominal noise according to

$$
\sigma_{\text{noise}} = \begin{cases} \frac{\sigma_{\text{noise}}^{\text{electronic}}}{\sqrt{(\sigma_{\text{noise}}^{\text{electronic}})^2 + (\sigma_{\text{noise}}^{\text{pile-up}})^2}} & (\text{2010 operations}) \\ \frac{\sigma_{\text{noise}}^{\text{electronic}}}{\sqrt{(\sigma_{\text{noise}}^{\text{electronic}})^2 + (\sigma_{\text{noise}}^{\text{pile-up}})^2}} & (\text{2011 operations}) \end{cases} .
$$

<footnote n="3">³ See the discussion of “truth jets” in Sect. 5.5 for the definition of stable particles.</footnote>
Jet reconstruction

<figure><img src="image_3.png" /><figcaption>Fig. 2: Overview of the ATLAS jet reconstruction. After the jet finding, the jet four momentum is defined as the four momentum sum of its constituents.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>Fig. 3: Overview of the ATLAS jet calibration scheme used for the 2011 dataset. The pile-up, absolute JES and the residual in situ corrections calibrate the scale of the jet, while the origin and the $\eta$ corrections affect the direction of the jet.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>Fig. 4: Average response of simulated jets formed from topo-clusters, calculated as defined in Eq. (1) and shown in (a) for the EM scale ($\mathcal{R}^{\text{EM}}$) and in (b) for the LCW scale ($\mathcal{R}^{\text{LCW}}$). The response is shown separately for various truth-jet energies as function of the uncorrected (detector) jet pseudorapidity $\eta_{\text{det}}$. Also indicated are the different calorimeter regions. The inverse of $\mathcal{R}^{\text{EM}}$ ($\mathcal{R}^{\text{LCW}}$) corresponds to the average jet energy scale correction for EM (LCW) in each $\eta_{\text{det}}$ bin. The results shown are based on the baseline PYTHIA inclusive jet sample.</figcaption></figure>
Here, $\sigma_{\text{noise}}^{\text{electronic}}$ is the electronic noise, and $\sigma_{\text{noise}}^{\text{pile-up}}$ the noise from pile-up, determined with MC simulations and corresponding to an average of eight additional proton-proton interactions per bunch crossing ($\mu = 8$) in 2011. The change of the total nominal noise $\sigma_{\text{noise}}$ and its dependence on the calorimeter region in ATLAS can be seen by comparing Figs. 1(a) and 1(b). In most calorimeter regions, the noise induced by pile-up is smaller than or of the same magnitude as the electronic noise, with the exception of the forward calorimeters, where

$$ \sigma_{\text{noise}}^{\text{pile-up}} \gg \sigma_{\text{noise}}^{\text{electronic}}. $$

The implicit noise suppression implemented by the topological cluster algorithm discussed above leads to significant improvements in the calorimeter performance for e.g. the energy and spatial resolutions in the presence of pile-up. On the other hand, contributions from larger negative and positive signal fluctuations introduced by pile-up can survive in a given event. They thus contribute to the sensitivity to pile-up observed in the jet response, in addition to the cell-level effects mentioned in Sect. 2.2.

## 5.2 Jet reconstruction and calibration

Jets are reconstructed using the anti-$k_t$ algorithm [1] with distance parameters $R = 0.4$ or $R = 0.6$, utilising the FASTJET software package [64,65]. The four-momentum scheme is used at each recombination step in the jet clustering. The total jet four-momentum is therefore defined as the sum of the four-momenta sum of all its constituents. The inputs to the jet algorithm are stable simulated particles (truth jets, see Sect. 5.5 for details), reconstructed tracks in the inner detector (track jets, see Ref. [3] and Sect. 5.4 for details) or energy deposits in the calorimeter (calorimeter jets, see below for details). A schematic overview of the ATLAS jet reconstruction is presented in Fig. 2.

The calorimeter jets are built from the topo-clusters entering as massless particles in the jet algorithm as discussed in the previous section. Only clusters with positive energy are considered. The topo-clusters are initially reconstructed at the EM scale [62,66–73], which correctly measures the energy deposited in the calorimeter by particles produced in electromagnetic showers. A second topo-cluster collection is built by calibrating the calorimeter cell such that the response of the calorimeter to hadrons is correctly reconstructed. This calibration uses the local cell signal weighting (LCW) method that aims at an improved resolution compared to the EM scale by correcting the signals from hadronic deposits, and thus reduces fluctuations due to the non-compensating nature of the ATLAS calorimeter. The LCW method first classifies topo-clusters as either electromagnetic or hadronic, primarily based on the measured energy density and the longitudinal shower depth. Energy corrections are derived according to this classification from single charged and neutral pion MC simulations. Dedicated corrections address effects of calorimeter non-compensation, signal losses due to noise threshold effects, and energy lost in non-instrumented regions close to the cluster [3].

Figure 3 shows an overview of the ATLAS calibration scheme for calorimeter jets used for the 2011 dataset, which restores the jet energy scale to that of jets reconstructed from sta-

ble simulated particles (truth particle level, see Sect. 5.5). This procedure consists of four steps as described below.

### 1 Pile-up correction

Jets formed from topo-clusters at the EM or LCW scale are first calibrated by applying a correction to account for the energy offset caused by pile-up interactions. The effects of pile-up on the jet energy scale are caused by both additional proton collisions in a recorded event (in-time pile-up) and by past and future collisions influencing the energy deposited in the current bunch-crossing (out-of-time pile-up), and are outlined in Sect. 6. This correction is derived from MC simulations as a function of the number of reconstructed primary vertices ($N_{\text{PV}}$, measuring the actual collisions in a given event) and the expected average number of interactions ($\mu$, sensitive to out-of-time pile-up) in bins of jet pseudorapidity and transverse momentum (see Section 6).

### 2 Origin correction

A correction to the calorimeter jet direction is applied that makes the jet pointing back to the primary event vertex instead of the nominal centre of the ATLAS detector.

### 3 Jet calibration based on MC simulations

Following the strategy presented in Ref. [3], the calibration of the energy and pseudorapidity of a reconstructed jet is a simple correction derived from the relation of these quantities to the corresponding ones of the matching truth jet (see Sect. 5.5) in MC simulations. It can be applied to jets formed from topo-clusters at EM or at LCW scale with the resulting jets being referred to as calibrated with the EM+JES or with the LCW+JES scheme. This first JES correction uses isolated jets from an inclusive jet MC sample including pile-up events (the baseline sample described in Sect. 3). Figure 4 shows the average energy response

$$ \Re^{\text{EM(LCW)}} = E_{\text{jet}}^{\text{EM(LCW)}} / E_{\text{jet}}^{\text{truth}}, \quad (1) $$

which is the inverse of the jet energy calibration function, for various jet energies as a function of the jet pseudorapidity $\eta_{\text{det}}$ measured in the detector frame of reference (see Sect. 5.6).

### 4 Residual in situ corrections

A residual correction derived in situ is applied as a last step to jets reconstructed in data. The derivation of this correction is described in in Sect. 7.

## 5.3 Jet quality selection

Jets with high transverse momenta produced in proton-proton collisions must be distinguished from background jet candidates not originating from hard-scattering events. A first strategy to select jets from collisions and to suppress background is presented in Ref. [3].

The main sources of potential background are:

1. Beam–gas events, where one proton of the beam collides with the residual gas within the beam pipe.

2. Beam-halo events, for example caused by interactions in the tertiary collimators in the beam-line far away from the ATLAS detector.
<figure><img src="image_6.png" /><figcaption>Fig. 5: Jet quality selection efficiency for anti-$k_t$ jets with $R = 0.4$ measured with a tag-and-probe technique as a function of $p_T^{jet}$ in various $\eta$ ranges, for the four sets of selection criteria. Only statistical uncertainties are shown. Differences between data and MC simulations are also shown.</figcaption></figure>
3. Cosmic-ray muons overlapping in-time with collision events.

4. Calorimeter noise.

The jet quality selection criteria should efficiently reject jets from these background processes while maintaining high efficiency for selecting jets produced in proton-proton collisions. Since the level and composition of background depend on the event topology and the jet kinematics, four sets of criteria called LOOSER, LOOSE, MEDIUM and TIGHT are introduced in Ref. [74]. They correspond to different levels of fake-jet rejection and jet selection efficiency, with the LOOSER criterion being the one with the highest jet selection efficiency while the TIGHT criterion is the one with the best rejection. The discrimination between jets coming from the collisions and background jet candidates is based on several pieces of experimental information, including the quality of the energy reconstruction at the cell level, jet energy deposits in the direction of the shower development, and reconstructed tracks matched to the jets.

The efficiencies of the jet selection criteria are measured using the tag-and-probe method described in Ref. [3]. The resulting efficiencies for anti-$k_t$ jets with $R = 0.4$ for all selection criteria are shown in Fig. 5. The jet selection efficiency of the LOOSER selection is greater than 99.8% over all calibrated transverse jet momenta $p_T^{jet}$ and $\eta$ bins. A slightly lower efficiency of about 1-2% is measured for the LOOSE selection, in particular at low $p_T^{jet}$ and for $2.5 < |\eta| < 3.6$. The MEDIUM and TIGHT selections have lower jet selection efficiencies mainly due to cuts on the jet charged fraction, which is the ratio of the scalar sum of the $p_T$ of all reconstructed tracks matching the jet, and the jet $p_T$ itself, see Ref. [74] for more details. For jets with $p_T^{jet} \approx 25$ GeV, the MEDIUM and TIGHT selections have inefficiencies of 4% and 15%, respectively. For $p_T^{jet} > 50$ GeV, the MEDIUM and TIGHT selections have efficiencies greater than 99% and 98%, respectively.

The event selection is based on the azimuthal distance between the probe and tag jet $\Delta\phi$(tag, probe) and the significance of the missing transverse momentum $E_T^{miss}$ [75] reconstructed for the event, which is measured by the ratio $E_T^{miss}/\sqrt{\Sigma E_T}$. Here $\Sigma E_T$ is the scalar transverse momentum sum of all particles, jets, and soft signals in the event. The angle $\Delta\phi$(tag, probe), $E_T^{miss}/\sqrt{\Sigma E_T}$, and the TIGHT selection of the reference (tag) jet are varied to study the systematic uncertainties. For the LOOSE and LOOSER selections, the jet selection efficiency is almost unchanged by varying the selection cuts, with variations of less than 0.05%. Slightly larger changes are observed for the two other selections, but they are not larger than 0.1% for the MEDIUM and 0.5% for the TIGHT selection.

The jet selection efficiency is also measured using a MC simulation sample. A very good agreement between data and simulation is observed for the LOOSER and LOOSE selections. Differences not larger than 0.2% and 1% are observed for the MEDIUM and TIGHT selections, respectively, for $p_T^{jet} > 40$ GeV. Larger differences are observed at lower $p_T^{jet}$, but they do not exceed 1%(2%) for the MEDIUM(TIGHT) selection.

## 5.4 Track jets

In addition to the previously described calorimeter jets reconstructed from topo-clusters, track jets in ATLAS are built from reconstructed charged particle tracks associated with the reconstructed primary collision vertex, which is defined by

$$ \sum(p_T^{\text{track}})^2 = \max. $$

Here $p_T^{\text{track}}$ is the transverse momentum of tracks pointing to a given vertex. The tracks associated with the primary vertex are required to have $p_T^{\text{track}} > 500$ MeV and to be within $|\eta| < 2.5$. Additional reconstruction quality criteria are applied, including the number of hits in the pixel detector (at least one) and in the silicon microstrip detector (at least six) of the ATLAS ID system. Further track selections are based on the transverse ($d_0$, perpendicular to the beam axis) and longitudinal ($z_0$, along the beam axis) impact parameters of the tracks measured with respect to the primary vertex ($|d_0| < 1.5$ mm, $|z_0 \sin\theta| < 1.5$ mm). Here $\theta$ is the polar angle of the track.

Generally, track jets used in the studies presented in this paper are reconstructed with the same configurations as calorimeter jets, i.e. using the anti-$k_t$ algorithm with $R = 0.4$ and $R = 0.6$. As only tracks originating from the hardest primary vertex in the collision event are used in the jet finding, the transverse momentum of any of these track jets provides a rather stable kinematic reference for matching calorimeter jets, as it is independent of the pile-up activity. Track jets can of course only be formed within the tracking detector coverage ($|\eta| < 2.5$), yielding an effective acceptance for track jets of $|\eta_{\text{trackjet}}| < 2.5 - R$.

Certain studies may require slight modifications of the track selection and the track-jet formation criteria and algorithms. Those are indicated in the respective descriptions of the applied methods. In particular, track jets may be further selected by requirements concerning the number of clustered tracks, the track-jet $p_T$, and the track-jet direction.

## 5.5 Truth jets

Truth jets can be formed from stable particles generated in MC simulations. In general those are particles with a lifetime $\tau$ defined by $c\tau > 10$ mm [76]. The jet definitions applied are the same as the ones used for calorimeter and track jets (anti-$k_t$ with distance parameters $R = 0.4$ and $R = 0.6$, respectively). If truth jets are employed as a reference for calibrations purposes in MC simulations, neither final-state muons nor neutrinos are included in the stable particles considered for its formation. The simulated calorimeter jets are calibrated with respect to truth jets consisting of stable particles leaving an observable signal (visible energy) in the detector.<fnref n="4" /> This is a particular useful strategy for inclusive jet measurements and the universal jet calibration discussed in this paper, but special truth-jet references including muons and/or neutrinos may be utilised as

<footnote n="4">⁴ Muons can generate an observable signal in some of the ATLAS calorimeters, but it is generally small and usually not proportional to the actual muon energy loss. Their contribution to the truth-jet energy, which can be large, is excluded to avoid biases and tails in the response function due to occasionally occurring high-$p_T$ muons in the MC-simulated calibration samples.</footnote>
well, in particular to understand the heavy-flavour jet response, as discussed in detail in Sect. 19.

## 5.6 Jet kinematics and directions

Kinematic properties of jets relevant for their use in final-state selections and final-state reconstruction are the transverse momentum $p_T$ and the rapidity $y$. The full reconstruction of the jet kinematics including these variables takes into account the physics frame of reference, which in ATLAS is defined event-by-event by the primary collision vertex discussed in Sect. 5.4.

On the other hand, many effects corrected by the various JES calibrations discussed in this paper are highly localised, i.e. they are due to specific detector features and inefficiencies at certain directions or ranges. The relevant directional variable to use as a basis for these corrections is then the detector pseudorapidity $\eta_{det}$, which is reconstructed in the nominal detector frame of reference in ATLAS, and is centred at the nominal collision vertex ($x = 0, y = 0, z = 0$).

Directional relations to jets, and e.g. between the constituents of jet and its principal axis, can then be measured either in the physics or the detector reference frame, with the choice depending on the analysis. In the physics reference frame ((y, $\phi$) space) the distance between any two objects is given by

$$ \Delta R = \sqrt{(\Delta y)^2 + (\Delta \phi)^2}, \qquad (2) $$

where $\Delta y$ is the rapidity distance and $\Delta \phi$ is the azimuthal distance between them. The same distance measured in the detector frame of reference (($\eta$, $\phi$) space) is calculated as

$$ \Delta R = \sqrt{(\Delta \eta)^2 + (\Delta \phi)^2}, \qquad (3) $$

where $\Delta \eta$ is the distance in pseudorapidity between any two objects. In case of jets and their constituents (topo-clusters or tracks), $\eta = \eta_{det}$ is used. All jet clustering algorithms used in ATLAS apply the physics frame distance in Eq. (2) in their distance evaluations, as jets are considered to be massive physical objects, and the jet clustering is intended to follow energy flow patterns introduced by the physics of parton showers, fragmentation, and hadronisation from a common (particle) source. In this context topo-clusters and reconstructed tracks are considered pseudo-particles representing the true particle flow within the limitations introduced by the respective detector acceptances and resolutions.

# 6 Jet energy correction for pile-up interactions

## 6.1 Pile-up correction method

The pile-up correction method applied to reconstructed jets in ATLAS is derived from MC simulations and validated with in situ and simulation based techniques. The approach is to calculate the amount of transverse momentum generated by pile-up in a jet in MC simulation, and subtract this offset $\mathcal{O}$ from the reconstructed jet $p_T^{jet}$ at any given signal scale (EM or LCW).

At least to first order, pile-up contributions to the jet signal can be considered stochastic and diffuse with respect to the true jet signal. Therefore, both in-time and out-of-time pile-up are expected to depend only on the past and present pile-up activity, with linear relations between the amount of activity and the pile-up signal.

## 6.2 Principal pile-up correction strategy

To characterise the in-time pile-up activity, the number of reconstructed primary vertices ($N_{PV}$) is used. The ATLAS tracking detector timing resolution allows the reconstruction of only in-time tracks and vertices, so that $N_{PV}$ provides a good measure of the actual number of proton-proton collisions in a recorded event.

For the out-of-time pile-up activity, the average number of interactions per bunch crossing ($\mu$) at the time of the recorded events provides a good estimator. It is derived by averaging the actual number of interactions per bunch crossing over a rather large window $\Delta t$ in time, which safely encompasses the time interval during which the ATLAS calorimeter signal is sensitive to the activity in the collision history ($\Delta t \gg 600$ ns for the liquid-argon calorimeters). The observable $\mu$ can be reconstructed from the average luminosity $L$ over this period $\Delta t$, the total inelastic proton-proton cross section ($\sigma_{inel} = 71.5$ mb [77]), the number of colliding bunches in LHC ($N_{bunch}$) and the LHC revolution frequency ($f_{LHC}$) (see Ref. [78] for details):

$$ \mu = \frac{L \times \sigma_{inel}}{N_{bunch} \times f_{LHC}}. $$

The MC-based jet calibration is derived for a given (reference) pile-up condition<fnref n="5" /> ($N_{PV}^{ref}$, $\mu^{ref}$) such that $\mathcal{O}(N_{PV} = N_{PV}^{ref}$, $\mu = \mu^{ref}) = 0$. As the amount of energy scattered into a jet by pile-up and the signal modification imposed by the pile-up history determine $\mathcal{O}$, a general dependence on the distances from the reference point is expected. From the nature of pile-up discussed earlier, the linear scaling of $\mathcal{O}$ in both $N_{PV}$ and $\mu$ provides the ansatz for a correction,

$$ \begin{aligned} \mathcal{O}(N_{PV}, \mu, \eta_{det}) &= p_T^{jet}(N_{PV}, \mu, \eta_{det}) - p_T^{truth} \\ &= \frac{\partial p_T}{\partial N_{PV}}(\eta_{det}) (N_{PV} - N_{PV}^{ref}) + \frac{\partial p_T}{\partial \mu}(\eta_{det}) (\mu - \mu^{ref}) \\ &= \alpha(\eta_{det}) \cdot (N_{PV} - N_{PV}^{ref}) + \beta(\eta_{det}) \cdot (\mu - \mu^{ref}) \end{aligned} \qquad (4) $$

Here, $p_T^{jet}(N_{PV}, \mu, \eta_{det})$ is the reconstructed transverse momentum of the jet (without the JES correction described in Section 5.2 applied) in a given pile-up condition ($N_{PV}, \mu$) and at a given direction $\eta_{det}$ in the detector. The true transverse momentum of the jet ($p_T^{truth}$) is available from the generated particle jet matching a reconstructed jet in MC simulations. The coefficients $\alpha(\eta_{det})$ and $\beta(\eta_{det})$ depend on $\eta_{det}$, as both in-time and out-of-time pile-up signal contributions manifest themselves differently in different calorimeter regions, according to the following influences:

<footnote n="5">⁵ The particular choice for a working point, here ($N_{PV}^{ref} = 4.9$, $\mu^{ref} = 5.4$), is arbitrary and bears no consequence for the correction method and its uncertainty.</footnote>
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>Fig. 6: The average reconstructed transverse momentum $p_{\text{T,EM}}^{\text{jet}}$ on EM scale for jets in MC simulations, as function of the number of reconstructed primary vertices $N_{\text{PV}}$ and $7.5 \le \mu < 8.5$, in various bins of truth-jet transverse momentum $p_{\text{T}}^{\text{truth}}$, for jets with (a) $R = 0.4$ and (b) $R = 0.6$. The dependence of $p_{\text{T,EM}}^{\text{jet}}$ on $N_{\text{PV}}$ in data, in bins of track-jet transverse momentum $p_{\text{T}}^{\text{track}}$, is shown in (c) for $R = 0.4$ jets, and in (d) for $R = 0.6$ jets.</figcaption></figure>
<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /><figcaption>Fig. 7: The average reconstructed jet transverse momentum $p_{T,EM}^{jet}$ on EM scale as function of the average number of collisions $\mu$ at a fixed number of primary vertices $N_{PV} = 6$, for truth jets in MC simulation (a) in the lowest bin of $p_{T}^{\text{truth}}$ and (b) in the lowest bin of track jet transverse momentum $p_{T}^{\text{trackjet}}$ considered in data.</figcaption></figure>

1. The energy flow from collisions into that region.

2. The calorimeter granularity and occupancy after topo-cluster reconstruction, leading to different acceptances at cluster level and different probabilities for multiple particle showers to overlap in a single cluster.

3. The effective sensitivity to out-of-time pile-up introduced by different calorimeter signal shapes.

The offset $\mathcal{O}$ can be determined in MC simulation for jets on the EM or the LCW scale by using the corresponding reconstructed transverse momentum on one of those scales, i.e. $p_{T}^{\text{jet}} = p_{T,EM}^{\text{jet}}$ or $p_{T}^{\text{jet}} = p_{T,LCW}^{\text{jet}}$ in Eq. (4), and $p_{T}^{\text{truth}}$. The particular choice of scale affects the magnitude of the coefficients and, therefore, the transverse momentum offset itself,

$$ \begin{aligned} \mathcal{O}^{\text{EM}} &\mapsto \{\alpha^{\text{EM}}(\eta_{\text{det}}), \beta^{\text{EM}}(\eta_{\text{det}})\} \\ \mathcal{O}^{\text{LCW}} &\mapsto \{\alpha^{\text{LCW}}(\eta_{\text{det}}), \beta^{\text{LCW}}(\eta_{\text{det}})\}. \end{aligned} $$

The corrected transverse momentum of the jet at either of the two scales ($p_{T,EM}^{\text{corr}}$ or $p_{T,LCW}^{\text{corr}}$) is then given by

$$ p_{T,EM}^{\text{corr}} = p_{T,EM}^{\text{jet}} - \mathcal{O}^{\text{EM}}(N_{\text{PV}}, \mu, \eta_{\text{det}}) \quad (5) $$

$$ p_{T,LCW}^{\text{corr}} = p_{T,LCW}^{\text{jet}} - \mathcal{O}^{\text{LCW}}(N_{\text{PV}}, \mu, \eta_{\text{det}}). \quad (6) $$

After applying the correction, the original $p_{T,EM}^{\text{jet}}$ and $p_{T,LCW}^{\text{jet}}$ dependence on $N_{PV}$ and $\mu$ is expected to vanish in the corresponding corrected $p_{T,EM}^{\text{corr}}$ and $p_{T,LCW}^{\text{corr}}$.

## 6.3 Derivation of pile-up correction parameters

Figures 6(a) and 6(b) show the dependence of $p_{T,EM}^{\text{jet}}$, and thus $\mathcal{O}^{\text{EM}}$, on $N_{PV}$. In this example, narrow ($R = 0.4, |\eta_{\text{det}}| < 2.1$)

and wide ($R = 0.6, |\eta_{\text{det}}| < 1.9$) central jets reconstructed in MC simulation are shown for events within a given range $7.5 \le \mu < 8.5$. The jet $p_{T}$ varies by $0.277 \pm 0.005$ GeV(in data) and $0.288 \pm 0.003$ GeV(in MC simulations) per primary vertex for jets with $R = 0.4$ and by $0.578 \pm 0.005$ GeV(in data) and $0.601 \pm 0.003$ GeV(in MC simulations) per primary vertex for jets with $R = 0.6$. The slopes $\alpha^{\text{EM}}$ are found to be independent of the true jet transverse momentum $p_{T}^{\text{truth}}$, as expected from the diffuse character of in-time pile-up signal contributions.

A qualitatively similar behaviour can be observed in collision data for calorimeter jets individually matched with track jets, the latter reconstructed as discussed in Sect. 5.4. The $N_{PV}$ dependence of $p_{T,EM}^{\text{jet}}$ can be measured in bins of the track-jet transverse momentum $p_{T}^{\text{trackjet}}$. Jets formed from tracks are much less sensitive to pile-up and can be used as a stable reference to investigate pile-up effects. Figures 6(c) and 6(d) show the results for the same calorimeter regions and out-of-time pile-up condition as for the MC-simulated jets in Figs. 6(a) and 6(b). The results shown in Fig. 6 also confirm the expectation that the contributions from in-time pile-up to the jet signal are larger for wider jets ($\alpha^{\text{EM}}(R = 0.6) > \alpha^{\text{EM}}(R = 0.4)$), but scale only approximately with the size of the jet catchment area [79] determined by the choice of distance parameter $R$ in the anti-$k_t$ algorithm.

The dependence of $p_{T,EM}^{\text{jet}}$ on $\mu$, for a fixed $N_{PV} = 6$, is shown in Fig. 7(a) for MC simulations using truth jets, and in Fig. 7(b) for collision data using track jets. The kinematic bins shown are the lowest bins considered, with $20 < p_{T}^{\text{truth}} < 25$ GeV and $20 < p_{T}^{\text{trackjet}} < 25$ GeV for MC simulations and data, respectively. The jet $p_{T}$ varies by $0.047 \pm 0.003$ GeV (in MC simulations) $0.105 \pm 0.003$ GeV (in data) per primary vertex for jets with $R = 0.4$
<figure><img src="image_13.png" /><figcaption>Fig. 8: The pile-up contribution per additional vertex, measured as $\alpha^{\text{EM}} = \partial p_{\text{T,EM}}^{\text{jet}}/\partial N_{\text{PV}}$, as function of $|\eta_{\text{det}}|$, for the various methods discussed in the text, for (a) $R=0.4$ and (b) $R=0.6$ jets. The contribution from $\mu$, calculated as $\beta^{\text{EM}} = \partial p_{\text{T,EM}}^{\text{jet}}/\partial \mu$ and displayed for the various methods as function of $|\eta_{\text{det}}|$, is shown for the two jet sizes in (c) and (d), respectively. The points for the determination of $\alpha^{\text{EM}}$ and $\beta^{\text{EM}}$ from MC simulations use the offset calculated from the reconstructed $p_{\text{T,EM}}^{\text{jet}}$ and the true (particle level) $p_{\text{T}}^{\text{truth}}$, as indicated in Eq. 4.</figcaption></figure>
The result confirms the expectations that the dependence of $p_{\text{T,EM}}^{\text{jet}}$ on the out-of-time pile-up is linear and significantly less than its dependence on the in-time pile-up contribution scaling with $N_{\text{PV}}$. Its magnitude is still different for jets with $R=0.6$, as the size of the jet catchment area again determines the absolute contribution to $p_{\text{T,EM}}^{\text{jet}}$.

The correction coefficients for jets calibrated with the EM+JES scheme, $\alpha^{\text{EM}}$ and $\beta^{\text{EM}}$, are both determined from MC simulations as functions of the jet direction $\eta_{\text{det}}$. For this, the $N_{\text{PV}}$ dependence of $p_{\text{T,EM}}^{\text{jet}}(\eta_{\text{det}})$ reconstructed in various bins of $\mu$ in the simulation is fitted and then averaged, yielding $\alpha^{\text{EM}}(\eta_{\text{det}})$. Accordingly and independently, the dependence of $p_{\text{T,EM}}^{\text{jet}}$ on $\mu$ is fitted in bins of $N_{\text{PV}}$, yielding the average $\beta^{\text{EM}}(\eta_{\text{det}})$, again using MC simulations. An identical procedure is used to find the correction functions $\alpha^{\text{LCW}}(\eta_{\text{det}})$ and $\beta^{\text{LCW}}(\eta_{\text{det}})$ for jets calibrated with the LCW+JES scheme.

The parameters $\alpha^{\text{EM}}(\alpha^{\text{LCW}})$ and $\beta^{\text{EM}}(\beta^{\text{LCW}})$ can be also measured with in situ techniques. This is discussed in Sect. 6.4.

## 6.4 Pile-up validation with in situ techniques and effect of out-of-time pile-up in different calorimeter regions

The parameters $\alpha^{\text{EM}}(\alpha^{\text{LCW}})$ and $\beta^{\text{EM}}(\beta^{\text{LCW}})$ can be measured in data with respect to a reference that is stable under pile-up using track jets or photons in $\gamma$-jet events as kinematic reference that does not depend on pile-up.

The variation of the $p_T$ balance $p_{\text{T,EM}}^{\text{jet}} - p_T^{\gamma} (p_{\text{T,LCW}}^{\text{jet}} - p_T^{\gamma})$ in $\gamma$-jet events can be used in data and MC simulation (similarly to the strategy discussed in Sect. 10), as a function of $N_{\text{PV}}$ and $\mu$. Figure 8 summarises $\alpha^{\text{EM}}(\eta_{\text{det}})$ and $\beta^{\text{EM}}(\eta_{\text{det}})$ determined with track jets and $\gamma$-jet events, and their dependence on $\eta_{\text{det}}$. Both methods suffer from lack of statistics or large systematic uncertainties in the 2011 data, but are used in data-to-MC comparisons to determine systematic uncertainties of the MC-based method (see the corresponding discussion in Sect. 16.2).

The decrease of $\beta^{\text{EM}}(\eta_{\text{det}})$ towards higher $\eta_{\text{det}}$, as shown in Figs. 8(c) and 8(d), indicates a decreasing signal contribution to $p_{\text{T,EM}}^{\text{jet}}$ per out-of-time pile-up interaction. For jets with $|\eta_{\text{det}}| > 1.5$, the offset is increasingly suppressed in the signal with increasing $\mu$ ($\beta^{\text{EM}}(\eta_{\text{det}}) < 0$). This constitutes a qualitative departure from the behaviour of the pile-up history contribution in the central region of ATLAS, where this out-of-time pile-up leads to systematically increasing signal contributions with increasing $\mu$.

This is a consequence of two effects. First, for $|\eta_{\text{det}}|$ larger than about 1.7 the hadronic calorimetry in ATLAS changes from the Tile calorimeter to the LAr end-cap (HEC) calorimeter. The Tile calorimeter has a unipolar and fast signal shape [20]. It has little sensitivity to out-of-time pile-up, with an approximate shape signal baseline of 150 ns. The out-of-time history manifests itself in this calorimeter as a small positive increase of its contribution to the jet signal with increasing $\mu$.

The HEC, on the other hand, has the typical ATLAS LAr calorimeter bipolar pulse shape with approximately 600 ns baseline. This leads to an increasing suppression of the contribution

from this calorimeter to the jet signal with increasing $\mu$, as more activity from the pile-up history increases the contribution weighted by the negative pulse shape.

Second, for $|\eta_{\text{det}}|$ larger than approximately 3.2, coverage is provided by the ATLAS forward calorimeter (FCal). While still a liquid-argon calorimeter, the FCal features a considerably faster signal due to very thin argon gaps. The shaping function for this signal is bipolar with a net zero integral and a similar positive shape as in other ATLAS liquid-argon calorimeters, but with a shorter overall pulse baseline (approximately 400 ns). Thus, the FCal shaping function has larger negative weights for out-of-time pile-up of up to 70% of the (positive) pulse peak height, as compared to typically 10% to 20% in the other LAr calorimeters [19]. These larger negative weights lead to larger signal suppression with increasing activity in the pile-up history and thus with increasing $\mu$.

## 7 In situ transverse momentum balance techniques

In this section an overview is given on how the data-to-MC differences are assessed using in situ techniques exploiting the transverse momentum balance between the jet and a well-measured reference object.

The calibration of jets in the forward region of the detector relative to jets in the central regions is discussed in more detail in Sect. 8. Jets in the central region are calibrated using photons or Z bosons as reference objects up to a transverse momentum of 800 GeV (see Sect. 9 and Sect. 10). Jets with higher $p_T$ are calibrated using a system of low-$p_T$ jets recoiling against a high-$p_T$ jet (see Sect. 11).

### 7.1 Relative in situ calibration between the central and forward rapidity regions

Transverse momentum balance in dijet events is exploited to study the pseudorapidity dependence of the jet response. A relative $\eta$-intercalibration is derived using the matrix method described in Ref. [3] to correct the jets in data for residual effects not captured by the initial calibration derived from MC simulations and based on truth jets. This method is applied for jets with $20 \le p_T^{\text{jet}} < 1500$ GeV and $|\eta_{\text{det}}| \le 4.5$. Jets up to $|\eta_{\text{det}}| = 2.8$ are calibrated using $|\eta_{\text{det}}| < 0.8$ as a reference region. For jets with $\eta_{\text{det}} > 2.8$ ($\eta_{\text{det}} < -2.8$), for which the uncertainty on the derived calibration becomes large, the calibration determined at $\eta_{\text{det}} = 2.8$ ($\eta_{\text{det}} = -2.8$) is used.<fnref n="6" /> Jets that fall in the reference region receive no additional correction on average. The $\eta$-intercalibration is applied to all jets prior to deriving the absolute calibration of the central region.

The largest uncertainty of the dijet balance technique is due to the modelling of the additional parton radiation altering the $p_T$ balance. This uncertainty is estimated using MC simulations employing the PYTHIA and HERWIG++ generators, respectively.

<footnote n="6">⁶ The relative jet response is measured independently for each $\eta_{\text{det}}$ hemisphere of the detector to accommodate asymmetries introduced by the actual collision vertex position during data taking.</footnote>
## 7.2 In situ calibration methods for the central rapidity region

The energy scale of jets is tested in situ using a well-calibrated object as reference. The following techniques are used for the central rapidity region $\eta_{det} < 1.2$:

### 1 Direct transverse momentum balance between a photon or a Z boson and a jet

Events with a photon or a Z boson and a recoiling jet are used to directly compare the transverse momentum of the jet to that of the photon or the Z boson (direct balance, DB). The data are compared to MC simulations in the jet pseudorapidity range $|\eta_{det}| < 1.2$. The $\gamma$-jet analysis covers a range in photon transverse momentum from 25 to 800 GeV, while the Z-jet analysis covers a range in Z transverse momentum from 15 to 200 GeV. However, only the direct transverse momentum balance between the Z and the jet is used in the derivation of the residual JES correction, as the method employing $p_T$ balance between a photon and the full hadronic recoil, rather than the jet (see item 2 below), is used in place of the direct $\gamma$-jet balance, see Sect. 13.5 for more details.

### 2 Transverse momentum balance between a photon and the hadronic recoil

The photon transverse momentum is balanced against the full hadronic recoil using the projection of the missing transverse momentum onto the photon direction. With this missing transverse momentum projection fraction (MPF) technique, the calorimeter response for the hadronic recoil is measured, which is independent of any jet definition. The comparison is done in the same kinematic region as the direct photon balance method.

### 3 Balance between a low-$p_T$ jet system and a high-$p_T$ jet

Jets at high $p_T$ can be balanced against a recoil system of low $p_T$ jets within $\eta_{det} < 2.8$ if the low $p_T$ jets are well calibrated using $\gamma$-jet or Z-jet in situ techniques. The multijet balance can be iterated several times to increase the non-leading (in terms of $p_T$) jets $p_T$ range beyond the values covered by $\gamma$-jet or Z-jet balance, and reaching higher $p_T$ of the leading jet, until statistical limitations preclude a precise measurement. This method can probe the jet energy scale up to the TeV regime.

In addition to the methods mentioned above, the mean transverse momentum sum of tracks within a cone around the jet direction provides an independent test of the calorimeter energy scale over the entire measured $p_T^{jet}$ range within the tracking acceptance. This method, described in Ref. [3], is used for the 2010 dataset and is also studied for the inclusive jet data sample in 2011. It is also used for b-jets (see Sect. 19). However, because of the relatively large associated systematic uncertainties, it is not included in the JES calibration derived from the combination of in situ methods for inclusive jets in 2011. This calibration can be constrained to much higher quality by applying the three methods described above.

<figure><img src="image_14.png" /></figure>

Fig. 9: Overview of the ($p_T^{avg}$, $\eta_{det}$) bins of the dijet balance measurements for jets reconstructed with distance parameter $R = 0.4$ calibrated using the EM+JES scheme. The solid lines indicate the ($p_T^{avg}$, $\eta_{det}^{probe}$) bin edges, and the points show the average transverse momentum and pseudorapidity of the probe jet within each bin. The measurements within the $\eta_{det}$ range spanned by the two thick, dashed lines are used to derive the residual calibration.

## 8 Relative forward-jet calibration using dijet events

The calibration of the forward detector can be performed by exploiting the transverse momentum balance in events with two jets at high transverse momentum. A well calibrated jet in the central part of the detector is balanced against a jet in the forward region.

Thus the whole detector acceptance in $\eta$ can be equalised as a function of $p_T^{jet}$. In addition to this simple approach, a matrix method is used where jets in all regions (and not only the central one) are used for the $\eta$-intercalibration.

In the following the results for the EM+JES scheme are discussed as an example. While the measured relative response can deviate by a few percent between the EM+JES and the LCW+JES calibration schemes, the ratio between data and Monte Carlo simulation agrees within a few permille.

### 8.1 Intercalibration using events with dijet topologies

#### 8.1.1 Intercalibration using a central reference region

The standard approach for $\eta$-intercalibration with dijet events is to use the central region of the calorimeters as the reference region, as described in Ref. [80]. The relative calorimeter response of jets in other calorimeter regions is measured by the $p_T$ balance between the reference jet (with $p_T^{ref}$) and the probe jet (with $p_T^{probe}$), exploiting the fact that these jets are expected to have equal $p_T$ due to transverse momentum conservation. The $p_T$ balance is expressed in terms of the asymmetry $A$,

$$A = \frac{p_T^{probe} - p_T^{ref}}{p_T^{avg}}, \qquad (7)$$
with $p_T^{\text{avg}} = (p_T^{\text{probe}} + p_T^{\text{ref}})/2$. The reference region is chosen as the central region of the barrel calorimeter, given by $|\eta_{\text{det}}| < 0.8$. If both jets fall into the reference region, each jet is used, in turn, to probe the other. As a consequence, the average asymmetry in the reference region will be zero by construction.

The asymmetry is then used to measure an $\eta$-intercalibration factor $c$ of the probe jet, or its response relative to the reference jet $1/c$, using the relation

$$ \frac{p_T^{\text{probe}}}{p_T^{\text{ref}}} = \frac{2+A}{2-A} = 1/c. \qquad (8) $$

The measurement of $c$ is performed in bins of jet $\eta_{\text{det}}$ and $p_T^{\text{avg}}$, where $\eta_{\text{det}}$ is defined as discussed in Sect. 5.6. Using the standard method outlined above, there is an asymmetry distribution $A_{ik}$ for each probe jet $\eta_{\text{det}}$ bin $i$ and each $p_T^{\text{avg}}$ bin $k$. An overview of the binning is given in Fig. 9 for jets with $R = 0.4$ calibrated with the EM+JES scheme. The same bins are used for jets calibrated with the EM+JES or LCW+JES scheme. However, the binning is changed for jets with $R = 0.6$ to take the different trigger thresholds into account. Intercalibration factors are calculated for each bin according to Eq. (8), resulting in

$$ c_{ik} = \frac{2 - \langle A_{ik} \rangle}{2 + \langle A_{ik} \rangle}, $$

where the $\langle A_{ik} \rangle$ is the mean value of the asymmetry distribution in each bin. The uncertainty on $\langle A_{ik} \rangle$ is taken to be the RMS/√N of each distribution. For the data, $N$ is the number of events in the bin, while for the MC sample the number of effective events $N_{\text{eff}}$ is used ($N = N_{\text{eff}}$) to incorporate MC event weights $w_k$,

$$ N_{\text{eff}} = (\sum w_k)^2 / \sum w_k^2. $$

Here the sums are running over all events of the MC sample. The above procedure is referred to as the *central reference method*.

### 8.1.2 Intercalibration using the matrix method

A disadvantage with the central reference method outlined above is that all events are required to have a jet in the central reference region. This results in a significant loss of event statistics, especially in the forward region, where the dijet cross section drops steeply as the rapidity interval between the jets increases. In order to use the full statistics, one can extend the central reference method by replacing the probe and reference jets by “left” and “right” jets, defined by $\eta_{\text{det}}^{\text{left}} < \eta_{\text{det}}^{\text{right}}$. Equations (7) and (8) then become

$$ A = \frac{p_T^{\text{left}} - p_T^{\text{right}}}{p_T^{\text{avg}}}, \quad \text{and} \quad \mathcal{R} = \frac{p_T^{\text{left}}}{p_T^{\text{right}}} = \frac{c^{\text{right}}}{c^{\text{left}}} = \frac{2+A}{2-A}, $$

where the term $\mathcal{R}$ denotes the ratio of the responses, and $c^{\text{left}}$ and $c^{\text{right}}$ are the $\eta$-intercalibration factors for the left and right jet, respectively.

This approach yields response ratio ($\mathcal{R}_{ijk}$) distributions with an average value $\langle R_{ijk} \rangle$, evaluated for each $\eta_{\text{det}}^{\text{left}}$ bin $i$, $\eta_{\text{det}}^{\text{right}}$ bin $k$.

The relative correction factor $c_{ik}$ for a given jet in $\eta_{\text{det}}$ bin $i$, with $i = 1...N$, and for a fixed $p_T^{\text{avg}}$ bin $k$ is then obtained by a minimisation procedure using a set of $N$ equations,

$$ S(c_{1k}, \dots, c_{Nk}) = \sum_{j=1}^{N} \sum_{i=1}^{j-1} \left( \frac{1}{\Delta \langle \mathcal{R}_{ijk} \rangle} (c_{ik} \langle \mathcal{R}_{ijk} \rangle - c_{jk}) \right)^2 + X(c_{1k}, \dots, c_{Nk}). \qquad (9) $$

Here $\Delta \langle \mathcal{R} \rangle$ is the statistical uncertainty of $\langle \mathcal{R} \rangle$ and the function $X(c_{ik})$ is used to quadratically suppress deviations from unity of the average corrections,<fnref n="7" />

$$ X(c_{1k}, \dots, c_{Nk}) = K \left( N^{-1} \sum_{i=1}^{N} c_{ik} - 1 \right)^2. $$

The value of the constant $K$ does not influence the solution as long as it is sufficiently large, e.g. $K \approx N_{\text{bins}}$, where $N_{\text{bins}}$ is the number of $\eta_{\text{det}}$ bins. The minimisation according to Eq. (9) is performed separately for each $p_T$ bin $k$, and the resulting calibration factors $c_i$ obtained in each $\eta_{\text{det}}$ bin $i$ are scaled such that the average calibration factor in the reference region $|\eta_{\text{det}}| < 0.8$ equals unity. This method is referred to as the *matrix method*.

## 8.2 Event selection for dijet analysis

### 8.2.1 Trigger selection

Events are retained from the calorimeter trigger stream using a combination of central ($|\eta_{\text{det}}| < 3.2$) and forward ($|\eta_{\text{det}}| > 3.1$) jet triggers [18].

The selection is designed such that the trigger efficiency for a specific region of $p_T^{\text{avg}}$ is greater than 99%, and approximately flat as a function of the pseudorapidity of the probe jet. Due to the different prescales for the central and forward jet triggers, the data collected by each trigger correspond to different integrated luminosities. To correctly normalise the data, events are assigned weights depending on the luminosity and the trigger decisions, according to the exclusion method described in Ref. [81].

### 8.2.2 Dataset and jet quality selection

All ATLAS sub-detectors are required to be operational and events are rejected if any data-quality issues are present. The leading two jets are required to fulfil the default set of jet quality criteria (see Sect. 5.3). A dead calorimeter region was present for a subset of the data. To remove any bias from this region, events are removed if any jets are reconstructed close to this region.

<footnote n="7">⁷ This term prevents the minimisation from choosing the trivial solution, which is all $c_{ik} = 0$.</footnote>
### 8.2.3 Dijet topology selection

In order to use the momentum balance of dijet events to measure the jet response, it is important that the events used have a relatively clean $2 \to 2$ topology. If a third jet is produced in the same hard-scatter proton-proton interaction, the balance between the leading (in $p_T$) two jets is affected. To enhance the number of events in the sample that have this $2 \to 2$ topology, selection criteria on the azimuthal angle $\Delta\phi(\text{jet1}, \text{jet2})$ between the two leading jets, and $p_T$ requirements on additional jets are applied. Table 1 summarises these topology selection criteria.

In addition, all jets used for balancing and topology selection have to originate from the hard-scattering vertex, and not from a vertex reconstructed from a pile-up interaction. For this, each jet considered is evaluated with respect to its jet vertex fraction (JVF), a likelihood measure estimating the vertex contribution to a jet [3]. To calculate JVF, reconstructed tracks originating from reconstructed primary vertices $i = 1, \dots, N_{\text{PV}}$ are matched to jets using an angular matching criterion in $(\eta, \phi)$ space of $\Delta R < 0.4$ with respect to the jet axis. The track parameters calculated at the distance of closest approach to the selected hard-scattering vertex are used for this matching. For each jet, the scalar sum of the $p_T$ of these matched tracks, $\Sigma_i$, is calculated for each vertex $i$ contributing to the jet. The JVF variable is then defined as the $p_T$ sum for the hard-scattering vertex, $\Sigma_0$, divided by the sum of $\Sigma_i$ over all primary vertices. Any jet that has $|\eta_{\text{det}}| < 2.5$ and JVF $> 0.6$ is classified as “vertex confirmed” since it is likely to originate from the hard-scattering vertex.

<table><caption>Table 1: Summary of the event topology selection criteria applied in this analysis. The symbols “jet1” and “jet2” refer to the leading two jets (two highest-$p_T$ jets), while “jet3” indicates the highest-$p_T$ sub-leading (third) jet in the event.</caption><thead><tr><td>Variable</td><td>Selection</td></tr></thead><tbody><tr><td>Δφ(jet1, jet2)</td><td>> 2.5 rad</td></tr><tr><td>p<sub>T</sub><sup>jet3</sup>, |η<sub>det</sub>(jet3)| < 2.5</td><td>&lt; max(0.25 p<sub>T</sub><sup>avg</sup>, 12 GeV)</td></tr><tr><td>p<sub>T</sub><sup>jet3</sup>, |η<sub>det</sub>(jet3)| > 2.5</td><td>&lt; max(0.20 p<sub>T</sub><sup>avg</sup>, 10 GeV)</td></tr><tr><td>JVF(jet3), |η<sub>det</sub>(jet3)| < 2.5</td><td>> 0.6</td></tr></tbody></table>

This selection differs from that used in previous studies [3] due to the much higher instantaneous luminosities experienced during data taking and the consequentially increasing pile-up. In the forward region $|\eta_{\text{det}}| > 2.5$, no tracking is available, and events containing any additional forward jet with significant $p_T$ are removed (see the third criteria in Table 1).

## 8.3 Dijet balance results

### 8.3.1 Binning of the balance measurements

An overview of the ($p_T^{\text{avg}}, \eta_{\text{det}}$) bins used in the analysis is presented in Fig. 9. All events falling in a given $p_T^{\text{avg}}$ bin are collected using a dedicated central and forward trigger combination. The statistics in each $p_T^{\text{avg}}$ bin are similar, except for the

highest and lowest bins which contain fewer events. The loss of statistical precision of the measurements for the lower $p_T^{\text{avg}}$ bins is introduced by a larger sensitivity to the inefficiency of the pile-up suppression strategy, which rejects relatively more events due to the kinematic overlap of the hard-scatter jets with jets from pile-up. In addition, the asymmetry distribution broadens due to worsening relative jet $p_T$ resolution, leading to larger fluctuations in this observable.

Each $p_T^{\text{avg}}$ bin is further divided into several $\eta_{\text{det}}$ bins. The $\eta_{\text{det}}$ binning is motivated by detector geometry and statistics.

### 8.3.2 Comparison of intercalibration methods

The relative jet response obtained with the matrix method is compared to the relative jet response obtained using the central reference method. Figures 10(a) and 10(b) show the jet response relative to central jets (1/c) for two $p_T^{\text{avg}}$ bins, $40 \le p_T^{\text{avg}} < 55$ GeV and $220 \le p_T^{\text{avg}} < 300$ GeV. In the most forward region at low $p_T$, the matrix method tends to give a slightly higher relative response compared to the central reference method (see Fig. 10(a)). However, the same relative shift is observed both for data and MC simulations, and consequently the data-to-MC ratios are consistent. The matrix method is therefore used to measure the relative response as it has better statistical precision.

### 8.3.3 Comparison of data with Monte Carlo simulation

Figure 11 shows the relative response obtained using the matrix method as a function of the jet pseudorapidity for data and MC simulations. Four different $p_T^{\text{avg}}$ regions are shown, $22 \le p_T^{\text{avg}} < 30$ GeV, $55 \le p_T^{\text{avg}} < 75$ GeV, $170 \le p_T^{\text{avg}} < 220$ GeV, and $600 \le p_T^{\text{avg}} < 800$ GeV. Figure 12 shows the relative response as a function of $p_T^{\text{avg}}$ for two representative $\eta_{\text{det}}$ bins, namely $-1.2 \le \eta_{\text{det}} < -0.8$ and $2.1 \le \eta_{\text{det}} < 2.8$. The general features of the response in data are reasonably well reproduced by the MC simulations. However, as observed in previous studies [3], the HERWIG++ MC generator predicts a higher relative response than PYTHIA for jets outside the central reference region ($|\eta_{\text{det}}| > 0.8$). Data tend to fall in-between the two predictions. This discrepancy was investigated and is observed both for truth jets built from stable particles (before any detector modelling), and also jets built from partons (before hadronisation). The differences therefore reflect a difference in physics modelling between the event generators, most likely due to the parton showering. The PYTHIA predictions are based upon a $p_T$-ordered parton shower whereas the HERWIG++ predictions are based on an angular-ordered parton shower.

For $p_T > 40$ GeV and $|\eta_{\text{det}}| < 2$, PYTHIA tends to agree better with data than HERWIG++ does. In the more forward region, the spread between the PYTHIA and HERWIG++ response predictions increases and reaches approximately 5% at $|\eta_{\text{det}}| = 4$. In the more forward region ($|\eta_{\text{det}}| > 3$) the relative response prediction of HERWIG++ generally agrees better with data than PYTHIA.
<figure><img src="image_15.png" /><figcaption>Fig. 10: Relative jet response (1/c) for anti-$k_t$ jets with $R = 0.4$ calibrated with the EM+JES scheme as a function of the probe jet pseudorapidity measured using the matrix and the central reference methods. Results are presented in (a) for $40 \le p_T^{\text{avg}} < 55$ GeV and in (b) for $220 \le p_T^{\text{avg}} < 300$ GeV. The lower parts of the figures show the ratios between relative response in data and MC.</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>Fig. 11: Relative jet response (1/c) as a function of the jet pseudorapidity for anti-$k_t$ jets with $R = 0.4$ calibrated with the EM+JES scheme, separately for (a) $22 \le p_T^{\text{avg}} < 30$ GeV, (b) $55 \le p_T^{\text{avg}} < 75$ GeV, (c) $170 \le p_T^{\text{avg}} < 220$ GeV and (d) $600 \le p_T^{\text{avg}} < 800$ GeV. The lower parts of the figures show the ratios between the data and MC relative response. These measurements are performed using the matrix method. The applied correction is shown as a thick line. The line is solid over the range where the measurements is used to constrain the calibration, and dashed in the range where extrapolation is applied.</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>Fig. 12: Relative jet response (1/c) as a function of the average jet $p_T$ of the dijet system for anti-$k_t$ jets with $R = 0.4$ calibrated with the EM+JES scheme, separately for (a) $-1.2 \le \eta_{\text{det}} < -0.8$ and (b) $2.1 \le \eta_{\text{det}} < 2.8$. The lower parts of the figures show the ratios between the data and MC relative response. The applied correction is shown as a thick line.</figcaption></figure>

### 8.3.4 Derivation of the residual correction

The residual calibration is derived from the data/PYTHIA ratio $C_i = c_i^{\text{data}}/c_i^{\text{PYTHIA}}$ of the measured $\eta$-intercalibration factors. PYTHIA is used as the reference as it is also used to obtain the initial (main) calibration, see Sect. 5. The correction is a function of jet $p_T$ and $\eta_{\text{det}}$ ($F_{\text{rel}}(p_T, \eta_{\text{det}})$) and is constructed by combining the $N_{\text{bins}}$ measurements of the ($p_T^{\text{avg}}, \eta_{\text{det}}$) bins using a two-dimensional Gaussian kernel, like

$$F_{\text{rel}}(p_T, \eta_{\text{det}}) = \frac{\sum_{i=1}^{N_{\text{bins}}} C_i w_i}{\sum_{i=1}^{N_{\text{bins}}} w_i},$$

with

$$w_i = \frac{1}{\Delta C_i^2} \times \text{Gaus} \left( \frac{\log p_T - \log \langle p_T^{\text{probe}} \rangle_i}{\sigma_{\log p_T}} \oplus \frac{\eta_{\text{det}} - \langle \eta_{\text{det}} \rangle_i}{\sigma_{\eta}} \right).$$

Here $i$ denotes the index of a ($p_T^{\text{avg}}, \eta_{\text{det}}$)-bin, $\Delta C_i$ is the statistical uncertainty of $C_i$, $\langle p_T^{\text{probe}} \rangle_i$ and $\langle \eta_{\text{det}} \rangle_i$ are the average $p_T$ and $\eta_{\text{det}}$ of the probe jets in the bin (see Fig. 9). The Gaussian function has a central value of zero and a width controlled by $\sigma_{\log p_T}$ and $\sigma_{\eta}$.

Only the measurements with $|\eta_{\text{det}}| < 2.8$ are included in the derivation of the correction function because of the large discrepancy between the modelled response of the MC simulation samples in the more forward region. This $\eta_{\text{det}}$ boundary is indicated by a thick, dashed line in Fig. 9. The residual correction is held fixed for pseudorapidities larger than those of the most forward measurements included ($|\eta_{\text{det}}| \approx 2.4$). All jets with a given $p_T$ and $|\eta_{\text{det}}| > 2.4$ will hence receive the same $\eta$-intercalibration correction. The kernel-width parameters used<fnref n="8" /> are found to capture the shape of the data-to-MC ratio, but at the

same time provide stability against statistical fluctuations. This choice introduces a stronger constraint across $p_T$. The resulting residual correction is shown as a thick line in the lower sections of Figs. 11 and 12. The line is solid over the range where the measurements is used to constrain the calibration, and dashed in the range where extrapolation is applied.

## 8.4 Systematic uncertainty

The observed difference in the relative response between data and MC simulations could be due to mis-modelling of physics or detector effects used in the simulation. Suppression and selection criteria used in the analysis (e.g. topology selection and radiation suppression) can also affect the response through their influence on the mean asymmetry. The systematic uncertainty is evaluated by considering the following effects:

1. Response modelling uncertainty.

2. Additional soft radiation.

3. Response dependence on the $\Delta\phi$ selection between the two leading jets.

4. Uncertainty due to trigger inefficiencies.

5. Influence of pile-up on the relative response.

6. Influence of the jet energy resolution (JER) on the response measurements.

All systematic uncertainties are derived as a function of $p_T$ and $|\eta_{\text{det}}|$. No statistically significant difference is observed for positive and negative $\eta_{\text{det}}$ for any of the uncertainties.

### 8.4.1 Modelling uncertainty

The two generators used for the MC simulation deviate in their predictions of the response for forward jets as discussed in Sect. 8.3.3. Since there is no a priori reason to trust one generator over the other, the full difference between the two predictions is used as the modelling uncertainty. This uncertainty is

<footnote n="8">⁸ A width of $\sigma_{\log p_T} = 0.25$ is used for the $p_T$ interpolation and $\sigma_{\eta} = 0.18$ for the $\eta_{\text{det}}$ interpolation.</footnote>
<figure><img src="image_18.png" /><figcaption>Fig. 13: Summary of uncertainties on the intercalibration as a function of the jet $\eta_{det}$ for anti-$k_t$ jets with $R = 0.4$ calibrated with the EM+JES scheme, separately for (a) $p_T = 35$ GeV and (b) $p_T = 350$ GeV. The individual components are added in quadrature to obtain the total uncertainty. The MC modelling uncertainty is the dominant component.</figcaption></figure>

the largest component of the intercalibration uncertainty. In the reference region ($|\eta_{det}| < 0.8$), no uncertainty is assigned. For $0.8 \le |\eta_{det}| < 2.4$, where data are corrected to the PYTHIA MC predictions, the full difference between PYTHIA and HERWIG is taken as the uncertainty. For $|\eta_{det}| > 2.4$, where the calibration is extrapolated, the uncertainty is taken as the difference between the calibrated data and either PYTHIA or HERWIG, whichever is larger.

### 8.4.2 Sub-leading jet radiation suppression

Additional radiation from sub-leading jets can affect the dijet balance. In order to mitigate these effects, selection criteria are imposed on the $p_T$ of any additional jets in an event as discussed in Sect. 8.2. To assess the uncertainties due to the radiation suppression, the selection criteria are varied for both data and MC simulations, and the calibration is re-evaluated. The uncertainty is taken as the fractional difference between the varied and nominal calibrations. Each of the three selection criteria are varied independently. The JVF requirement is changed by $\pm 0.2$ from its nominal value (0.6) for central jets, and the fractional amount of $p_T$ carried by the third jet relative to $p_T^{avg}$ is varied by $\pm 10\%$. Finally, the minimum $p_T$ cutoff is changed by $\pm 2$ GeV.

### 8.4.3 $\Delta\phi$(jet1, jet2) event selection

The event topology selection requires that the two leading jets have a $\Delta\phi$ separation greater than 2.5 rad. In order to assess the influence of this selection on the $p_T$ balance, the residual calibration is re-derived twice after shifting the selection criterion by $\pm 0.4$ rad ($\Delta\phi$(jet1, jet2) < ($2.5 \pm 0.4$) rad), separately in either direction. The largest difference between the shifted and nominal calibrations is taken as the uncertainty.

### 8.4.4 Trigger efficiencies

Trigger biases can be introduced if the trigger selection, which is applied only to data, is not fully efficient. To assess the uncertainty associated with the small inefficiency in the trigger, the measured efficiencies are applied to the MC samples. The effect on the MC response is found to be negligible in comparison to the other sources, even when exaggerating the effect by shifting the measured efficiency curves to reach the plateau 10% earlier in $p_T$. This uncertainty is hence neglected.

### 8.4.5 Impact of pile-up interactions

The influence of pile-up on the relative response is evaluated. To assess the magnitude of the effect, the differences between low and high pile-up subsets are investigated. Two different selections are used, high and low $\mu$ subsets ($\mu < 7$ and $\mu \ge 7$), and high and low $N_{PV}$ subsets ($N_{PV} < 5$ and $N_{PV} \ge 5$). The discrepancies observed are well within the systematic uncertainty for the pile-up correction itself (see Sect. 16). Therefore, no further contribution from pile-up is included in the evaluation of the full systematic uncertainty of the $\eta$-intercalibration.

### 8.4.6 Jet resolution uncertainty

The jet energy resolution (JER) [82] in the MC simulation is comparable to the resolution observed in data. To assess the impact of the JER on the $p_T$ balance, a smearing factor is applied as a scale factor to the MC jets, which results in an increased jet resolution consistent with the JER measured in data plus its error. It is randomly sampled from a Gaussian with width

$$ \sigma = \sqrt{(\sigma_{\text{data}} + \Delta\sigma_{\text{data}})^2 - \sigma_{\text{data}}^2} \quad (10) $$

where $\sigma_{\text{data}}$ is the measured jet resolution in data and $\Delta\sigma_{\text{data}}$ is the corresponding uncertainty. The difference between the nominal and smeared MC results is taken as the JER systematic uncertainty.
## 8.5 Summary of the η-intercalibration and its uncertainties

The pseudorapidity dependence of the jet response is analysed using dijet pseudorapidity η-intercalibration techniques. A residual $p_T$ and $\eta_{det}$ dependent response correction is derived with a matrix method for jets with $|\eta_{det}| < 2.4$. The correction is applied to data to correct for effects not captured by the default MC-derived calibration. The correction to the jet response is measured to be approximately +1% at $|\eta_{det}| = 1.0$ and falling to -3% and to -1% for $|\eta_{det}| = 2.4$ and beyond. The total systematic uncertainty is obtained as the quadratic sum of the various components mentioned. Figure 13 presents a summary of the uncertainties as a function of $\eta_{det}$ for two representative values of jet transverse momentum, namely $p_T = 35$ GeV and $p_T = 350$ GeV. The uncertainty is parameterised in the same way as the correction as described in Sect. 8.3.4. There is no strong variation of the uncertainties as a function of jet $p_T$. For a $p_T = 25$ GeV jet, the uncertainty is about 1% at $|\eta_{det}| = 1.0$, 3% at $|\eta_{det}| = 2.0$ and about 5% for $|\eta_{det}| > 3.0$. The uncertainty is below 1% for $p_T = 500$ GeV jets with $|\eta_{det}| < 2$.

<figure><img src="image_19.png" /><figcaption>Fig. 14: The Z boson $p_T$ distribution in selected Z events. Data and prediction from the $Z \to ee$ PYTHIA simulation, normalised to the observed number of events, are compared.</figcaption></figure>

## 9 Jet energy calibration using Z-jet events

This section presents results based on events where a Z boson decaying to an $e^+e^-$ pair is produced together with a jet, which balance each other in the transverse plane. The $p_T$ balance is compared in data and in MC simulations, and a study of systematic uncertainties on the data-to-MC ratio is carried out. The results from a similar study with $\gamma$-jet events are discussed in Sect. 10.

The advantage of Z-jet events is the possibility of probing low-$p_T$ jets, which are difficult to reach with $\gamma$-jet events due to trigger thresholds and background contamination in that region. On the other hand, $\gamma$-jet events benefit from larger statistics for $p_T$ above 150 GeV. In the Z-jet and $\gamma$-jet analyses, jets with a pseudorapidity $|\eta_{det}| < 1.2$ are probed.

### 9.1 Description of the $p_T$ balance method

In events where one Z boson and only one jet are produced, the jet recoils against the Z boson, ensuring approximate momentum balance between them in the transverse plane. The direct $p_T$ balance technique exploits this relationship in order to improve the jet energy calibration.

If the Z boson decays into electrons, its four-momentum is reconstructed using the electrons, which are accurately measured in the electromagnetic calorimeter and the inner detector [73]. Ideally, if the jet includes all the particles that recoil against the Z boson, and if the electron energies are perfectly measured, the response of the jet in the calorimeters can be determined by using $p_T^Z$ as the reference truth-jet $p_T$. However, this measurement is affected by the following:

1. Uncertainty on the electron energy measurements.

2. Particles contributing to the $p_T$ balance that are not included in the jet cone (out-of-cone radiation).

3. Additional parton radiation contributing to the recoil against the Z boson.

4. Contribution from the underlying event.

5. In-time and out-of-time pile-up.

Therefore, the direct $p_T$ balance between a Z boson and a jet ($p_T^{\text{jet}}/p_T^{\text{ref}}$) is not used to estimate the jet response, but only to assess how well the MC simulation can reproduce the data.

To at least partly reduce the effect of additional parton radiation perpendicular to the jet axis in the transverse plane, a reference $p_T^{\text{ref}} = p_T^Z \times |\cos(\Delta\phi(\text{jet}, Z))|$ is constructed from the azimuthal angle $\Delta\phi(\text{jet}, Z)$ between the Z boson and the jet, and the Z boson transverse momentum $p_T^Z$.

The jet calibration in the data is then adjusted using the data-to-MC comparison of the $p_T^{\text{jet}}/p_T^{\text{ref}}$ ratio for the two jet calibration schemes EM+JES and LCW+JES described in Sect. 5. The effects altering this ratio are evaluated by changing kinematic and topological selections and MC event generators and other modelling parameters. In particular the extrapolation of the $\Delta\phi(\text{jet}, Z)$ dependence of $p_T^{\text{jet}}/p_T^{\text{ref}}$ to the least radiation-biased regime ($\Delta\phi(\text{jet}, Z) = \pi$) is sensitive to the MC-modeling quality and is investigated with data-to-MC comparisons.

### 9.2 Selection of Z-jet events

Events are selected online using a trigger logic that requires the presence of at least one well-identified electron with transverse energy ($E_T^e$) above 20 GeV (or 22 GeV, depending on the data-taking period) or two well-identified electrons with $E_T^e > 12$ GeV, in the region $|\eta| < 2.5$ [83]. Events are also required to have a primary hard-scattering vertex, as defined in Sect. 5.4, with at least three tracks associated to it. This renders the contribution from fake vertices due to beam backgrounds negligible.
<figure><img src="image_20.png" /><figcaption>Fig. 15: The $p_T^{\text{jet}}/p_T^{\text{ref}}$ distribution in the data for $20 \le p_T^{\text{ref}} < 25$ GeV and $\pi - \Delta\phi(\text{jet},Z) < 0.05$ is shown in (a). The black solid histogram represents the fit function, a Poisson distribution extended to non-integer values, multiplied by a turn-on curve. The value used in each bin is the mean value of that function in the bin. The gray dashed line shows the underlying Poisson distribution, from which the mean is taken as the measurement of the $p_T$ balance. In (b) the measured widths of the $p_T^{\text{jet}}/p_T^{\text{ref}}$ distributions as a function of $p_T^{\text{ref}}$ is shown for data and MC simulations, for events with $\pi - \Delta\phi(\text{jet},Z) < 0.05$. The fitted functional form defined by Eq. (11) is superimposed. In both figures, anti-$k_t$ jets with distance parameter $R = 0.4$ calibrated with the EM+JES scheme are used. Only statistical uncertainties are shown.</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>Fig. 16: Mean $p_T^{\text{jet}}/p_T^{\text{ref}}$ balance vs $\Delta\phi(\text{jet},Z)$ for $20 \le p_T^{\text{ref}} < 25$ GeV in the data and in the simulation. A linear function used to extrapolate the balance to $\Delta\phi = \pi$ is superimposed. Anti-$k_t$ jets with distance parameter $R = 0.4$ calibrated with the EM+JES scheme are used. Only statistical uncertainties are shown.</figcaption></figure>

Details of electron reconstruction and identification can be found in Ref. [73]. Three levels of electron identification qual-

ity are defined, based on different requirements on shower sha-pes, track quality, and track–cluster matching. The intermediate one (“medium”) is used in this analysis.

Events are required to contain exactly two such electron candidates with $E_T^e > 20$ GeV and pseudorapidity in the range $|\eta^e| < 2.47$, where the transition region between calorimeter sections $1.37 < |\eta^e| < 1.52$ is excluded, as well as small re-gions where an accurate energy measurement is not possible due to temporary hardware failures. If these electrons have op-posite-sign charge, and yield a combined invariant mass in the range $66 < M_{e^+e^-} < 116$ GeV, the event is kept and the four-momentum of the Z boson candidate is reconstructed from the four-momenta of the two electrons. The transverse momentum distribution of these Z boson candidates is shown in Fig. 14.

All jets within the full calorimeter acceptance and with a JES-corrected transverse momentum $p_T^{\text{jet}} > 12$ GeV are consid-ered. For each jet the JVF (see Sect. 8.2.3) is used to estimate the degree of pile-up contamination of a jet based on the vertex information. The highest-$p_T$ (leading) jet must pass the quality criteria described in Sect. 5.3, have a JVF > 0.5, and be in the fiducial region $|\eta| < 1.2$.

Furthermore, the leading jet is required to be isolated from the two electrons stemming from the Z boson. The distance $\Delta R$ between the jet and each of the two electrons in $(\eta, \phi)$ space, measured according to Eq. (3) in Sect. 5.6, is required to be $\Delta R > 0.35$ (0.5) for anti-$k_t$ jets with $R = 0.4$ (0.6).

The presence of additional high-$p_T$ parton radiation alter-ing the balance between the Z boson and the leading jet is sup-pressed by requiring that the next-highest-$p_T$ (sub-leading) jet
has a calibrated $p_T$ less than 20% of the $p_T$ of the Z boson, with a minimal $p_T$ of 12 GeV. For sub-leading jets within the tracking acceptance, this cut is only applied if the jet has a JVF > 0.75. A summary of the event selection is presented in Table 2.

<table><caption>Table 2: Summary of the event selection criteria applied in the Z-jet analysis.</caption><thead><tr><th>Variable</th><th>Selection</th></tr></thead><tbody><tr><td>$E_T^e$</td><td>&gt; 20 GeV</td></tr><tr><td>$|\eta^e|$</td><td>&lt; 2.47<br>(excluding calorimeter transition regions)</td></tr><tr><td>$p_T^{jet}$</td><td>&gt; 12 GeV</td></tr><tr><td>$|\eta^{jet}|$</td><td>&lt; 1.2</td></tr><tr><td>$M_{e^+e^-}$</td><td>66 &lt; $M_{e^+e^-}$ &lt; 116 GeV</td></tr><tr><td>$\Delta R(\text{jet, electrons})$</td><td>&gt; 0.35 (0.5), anti-$k_t$ jets with $R = 0.4$ (0.6)</td></tr><tr><td>$p_T^{jet2}/p_T^Z$</td><td>&lt; 0.2</td></tr></tbody></table>

## 9.3 Measurement of the $p_T$ balance

The mean value of the $p_T^{jet}/p_T^{ref}$ ratio distribution is computed in bins of $p_T^{ref}$ and $\Delta\phi(\text{jet}, Z)$. This mean value is obtained with two methods, depending on $p_T^{ref}$.

1. In the low-$p_T^{ref}$ region ($17 \le p_T^{ref} < 35$ GeV), it is obtained by a maximum likelihood fit applied to the distribution of the $p_T^{jet}/p_T^{ref}$ ratio. The function used, hereafter denoted as the “fit function”, is a Poisson distribution extended to non-integer values<fnref n="9" /> and multiplied by a turn-on curve to model the effect of the $p_T^{jet}$ threshold, as depicted in Fig. 15(a). For a given [$p_T^{ref,min}, p_T^{ref,max}$] bin, the turn-on curve is equal to 1 above 12 GeV/$p_T^{ref,min}$ and equal to 0 below 12 GeV/$p_T^{ref,max}$. A linear function is used to interpolate the turn-on between these two values. The mean value of the underlying Poisson distribution is taken as the mean $p_T$ balance. A fit is preferred to an arithmetic mean calculation because of the jet $p_T$ cut, which biases the mean value of the balance distribution at low $p_T^{ref}$ due to the jet energy resolution [84].

2. For larger $p_T^{ref}$ ($p_T^{ref} \ge 35$ GeV), the arithmetic mean calculation is not sensitive to the jet threshold, and it gives results equivalent to those obtained with a fit. In this $p_T^{ref}$ region, an arithmetic mean is therefore used as it leads to smaller uncertainties.

In the region where the fit is used, $17 \le p_T^{ref} < 35$ GeV, the fit is actually performed twice, in order to reduce the impact of statistical fluctuations:

1. In each bin of $p_T^{ref}$ and $\Delta\phi$, the mean and the width of the Poisson distribution are fitted simultaneously.

2. The distribution of the widths is parameterised as a function of $p_T^{ref}$ in each $\Delta\phi$ bin according to:

$$w(p_T^{\text{ref}}) = \frac{a}{p_T^{\text{ref}}} \oplus \frac{b}{\sqrt{p_T^{\text{ref}}}} \oplus c. \quad (11)$$

The parameters $a$, $b$, and $c$ are obtained from a fit to the widths of the fitted Poisson distributions for $p_T < 35$ GeV and to the arithmetic RMS for larger $p_T$ (see Fig. 15(b)). It is emphasised that this measured width can not directly be compared to the resolution determined in Ref. [84], since no extrapolation to a topology without radiation is performed here.

3. The fits to the $p_T^{jet}/p_T^{ref}$ distributions are repeated, but now with the widths fixed to the values resulting from the parameterisations.

In order to estimate the mean balance for a topology where the jet and the Z boson are back-to-back, the mean balances in $\Delta\phi$ bins are extrapolated to $\Delta\phi = \pi$ for each $p_T^{ref}$ bin, using a linear function (see Fig. 16). This extrapolation reduces the sensitivity of the mean balance to additional parton radiation transverse to the leading jet axis, as discussed earlier in Sect. 9.1. The extrapolated mean balances for the data and MC-simulated samples generated by PYTHIA are shown in Fig. 17 for anti-$k_t$ jets with distance parameters of $R = 0.4$ and $R = 0.6$, calibrated with the EM+JES scheme. The mean balance obtained for jets with $R = 0.6$ is larger compared to jets with $R = 0.4$, which is a direct consequence of the larger jet size, and has smaller variations with the transverse momentum.

## 9.4 Measuring out-of-cone radiation and underlying event contributions

The transverse momentum of the Z boson is only approximately equal to the transverse momentum of the truth jet, because of out-of-cone radiation and contributions from the underlying event:

1. The Z boson balances against all particles inside and outside the jet cone, whereas the truth jet clusters particles inside the jet cone only.

2. The truth jet's $p_T$ includes any UE particles that are clustered in the jet, whereas the Z boson's $p_T$ receives almost no such contribution.

These two contributions are estimated by measuring the transverse momentum profile of tracks around the leading jet axis (see Fig. 18). Tracks associated to the hard-scattering vertex are used instead of clusters of calorimeter cells in order to reduce the sensitivity to pile-up interactions. Tracks associated to the two electrons stemming from the Z boson are removed when computing the transverse momentum profiles.

<footnote n="9">⁹ This continuous Poisson function is obtained by extending the discrete Poisson distribution to real values by replacing the factorials in the discrete Poisson function with Euler's Gamma function. This function has only one free parameter ($\lambda$). A linear transformation of the x-scale ($\chi' = a * x$) is introduced and the mean and width of this function are expressed in terms of $\lambda$ and $a$.</footnote>
<figure><img src="image_22.png" /></figure>

<figure><img src="image_23.png" /><figcaption>Fig. 17: Mean $p_T$ balance obtained in the data and with the PYTHIA simulation. Results for anti-$k_t$ jets with distance parameter (a) $R = 0.4$ and (b) $R = 0.6$ calibrated with the EM+JES scheme are shown. Only statistical uncertainties are shown.</figcaption></figure>

<figure><img src="image_24.png" /><figcaption>Fig. 18: Transverse momentum profile of tracks around the lea-ding jet axis for events with $20 \le p_T^{\text{ref}} < 25$ GeV in Z-jet events. The radii $R$ and $R_0$ are those used in Eq. (12) to define the IC and IC+OC regions. The hatched area indicates the contribu-tion from the underlying event (UE). Anti-$k_t$ jets with $R = 0.4$ calibrated with the EM+JES scheme are considered.</figcaption></figure>

A factor is calculated from the out-of-cone and underlying event contributions:

$$k_{\text{OOC}} = \frac{p_{\text{T}}^{\text{IC,ALL}}}{p_{\text{T}}^{\text{IC+OC,ALL}} - p_{\text{T}}^{\text{IC+OC,UE}}}, \quad (12)$$

where $p_T^{\text{IC,ALL}}$ is the average scalar $p_T$ sum of all the tracks inside the jet cone with radius $R$, $p_T^{\text{IC+OC,ALL}}$ is the average

scalar $p_T$ sum of all the tracks inside and outside the jet cone, and $p_T^{\text{IC+OC,UE}}$ is the average contribution of the underlying event to $p_T^{\text{IC+OC,ALL}}$. The transverse momenta $p_T^{\text{IC+OC,ALL}}$ and $p_T^{\text{IC+OC,UE}}$ are estimated in a cone of radius $R_0$, above which only the UE contributes to $p_T^{\text{IC+OC,ALL}}$, and from where the transverse momentum density is constant (see Fig. 18). In practice, $R_0$ is the value where the logarithmic derivative of $k_{\text{OOC}}$ with respect to $R_0$ is equal to 0.05.

## 9.5 Systematic uncertainties

The differences between the balances observed in data and those observed in MC simulations may be due to physics or detector effects directly influencing the calorimeter response to jets (e.g. fragmentation or material in front of the calorimeter), which may not be correctly modelled in the simulation. They can also be due to effects that have an influence on the direct $p_T$ balance method itself, e.g. the estimation of the mean balance or higher-order parton emissions. For a more detailed evaluation of the systematic uncertainties, the following steps are performed:

1. The uncertainty on the width parameterisation is propagated to the mean estimation.

2. The fit range used for the $\Delta\phi$ extrapolation is varied.

3. The cut on sub-leading jets is varied to assess the effect of additional high-$p_T$ parton radiation altering the balance.

4. The effect of soft particles produced outside the jet cone and the underlying event contribution to the jet energy is compared in data and simulation.

5. The impact of additional (pile-up) interactions is studied.

6. The uncertainty in the electron energy measurement is propagated to the $p_T$ balance.
<figure><img src="image_25.png" /><figcaption>Fig. 19: Summary of the Z-jet uncertainties on the data-to-MC ratio of the mean $p_T$ balance, for anti-$k_t$ jets with distance parameter (a) $R = 0.4$ and (b) $R = 0.6$ calibrated with the EM+JES scheme.</figcaption></figure>

7. The results obtained with PYTHIA and ALPGEN+HERWIG are compared.

### 9.5.1 Fitting procedure

For $p_T^{\text{ref}} < 35$ GeV, the mean balance in a given bin of $p_T^{\text{ref}}$ and $\Delta\phi$ is first obtained using the nominal parameterised width given in Eq. (11). The fit is then performed again with a larger and a smaller width according to the uncertainty on the parameterisation. The four differences obtained in the resulting mean balances for the up and down variations in data and Monte Carlo simulation are propagated independently, after $\Delta\phi$ extrapolation, to the data-to-MC ratio. The two positive and two negative deviations are both summed in quadrature and the final uncertainty is taken as the average of the absolute values of the two deviations.

### 9.5.2 Extrapolation procedure

The nominal extrapolated balance is determined with a linear fit from $\Delta\phi = \pi - 0.3$ to $\Delta\phi = \pi$. The lower limit is decreased to $\pi - 0.4$ and increased to $\pi - 0.2$, and the average of the absolute values of the two deviations is taken as a systematic uncertainty on the data-to-MC ratio.

### 9.5.3 Additional radiation suppression

While the extrapolation of the $p_T$ balance in $\Delta\phi$ attempts to reduce the effect of radiation perpendicular to the jet axis at angular scales within the range from $[\pi - 0.3, \pi]$, additional radiation not reflected by the $\Delta\phi$ measurement and extrapolation can still occur and modify the $p_T$ balance between the Z boson

and the leading jet with respect to expectations for true back-to-back topologies. Therefore, events with energetic sub-leading jets are vetoed. Systematic uncertainties associated with this second jet veto are studied, and the mean $p_T$ balances in the data and the simulation are compared when applying different second jet vetoes. The nominal

$$p_T^{\text{jet2,nom}} = \max\{12 \text{ GeV}, 0.2 \times p_T^Z\}$$

is varied up and down to

$$p_T^{\text{jet2,nom}\uparrow} = p_T^{\text{jet2,nom}} + 0.1 \times p_T^Z \quad (\text{down})$$

$$p_T^{\text{jet2,nom}\downarrow} = \max\{10 \text{ GeV}, 0.1 \times p_T^Z\} \quad (\text{up}).$$

The average of the absolute values of the two deviations is taken as a systematic uncertainty on the data-to-MC ratio.

### 9.5.4 Out-of-cone radiation and underlying event

This $k_{\text{OOC}}$ factor defined in Eq. (12) and measured as described in Sect. 9.4 indicates how the Z boson's $p_T$ differs from the truth jet's $p_T$. In order to evaluate the systematic uncertainties coming from out-of-cone radiation and the underlying event, this factor is applied to the Z boson's $p_T$. It is measured in the data and in the simulation in bins of $p_T^{\text{ref}}$. Its value depends on the $p_T$ as well as on the jet size. For jets with $R = 0.4$, $k_{\text{OOC}}$ increases from about 0.93 at low $p_T$ to about 0.99 at high $p_T$. For jets with $R = 0.6$, it varies between 1 and 1.02 without any systematic $p_T$ dependence. A modified data-to-MC ratio of the balance is calculated using the $k_{\text{OOC}}$ factors and the difference with respect to the nominal ratio is taken as a systematic uncertainty.
<figure><img src="image_26.png" /><figcaption>Fig. 20: Data-to-MC ratio of the mean $p_T$ balance for Z-jet events as a function of $p_T^{\text{ref}}$ for anti-$k_t$ jets with distance parameter ((a), (c)) $R = 0.4$ and ((b), (d)) $R = 0.6$ calibrated with the ((a), (b)) EM+JES and the ((c), (d)) LCW+JES schemes. The total uncertainty on this ratio is depicted by grey bands. Dashed lines show the -1%, -2%, and -5% shifts.</figcaption></figure>

### 9.5.5 Impact of additional pile-up interactions

The impact of in-time and out-of-time pile-up is studied by comparing the $p_T$ balance in two samples with different numbers of primary vertices ($N_{PV} \le 5$ and $N_{PV} > 5$), and two samples with different average number of interactions per bunch crossing ($\mu < 8$ and $\mu > 8$). The differences observed between the samples are small compared to the uncertainty on the pile-up offset correction (see Sect. 6.4). Therefore, they are not taken into account in this analysis in order to avoid double-counting between the different steps of the jet calibration procedure.

The direct impact of additional interactions on the leading jet is also studied by relaxing the JVF cut, introduced in Sect.

9.2, for the leading jet. The difference with respect to the nominal result is taken as an additional uncertainty.

### 9.5.6 Electron energy scale

The $p_T$ of the Z boson, measured from the energy of the electrons, is used as a reference to probe the jet energy scale. The electron energy is shifted upwards and downwards according to the uncertainty on its measurement [73], updated using data recorded in 2011.
### 9.5.7 Impact of the Monte Carlo generator

The mean balances are obtained from PYTHIA and ALPGEN samples, using the procedure described in Sect. 9.3. The difference between the data-to-PYTHIA and the data-to-ALPGEN ratios is taken as a systematic uncertainty. The ALPGEN MC generator uses different theoretical models for all steps of the event generation and therefore gives a reasonable estimate of the systematic variations. However, the possible compensation of modelling effects that shift the jet response in opposite directions cannot be excluded. To reduce the impact of statistical fluctuation the first three bins are merged, since they give the same result within their statistical uncertainties.

### 9.5.8 Summary of systematic uncertainties

Additional sources of uncertainties are considered:

1. The main background to Z-jet events is from multijet events, and its fraction in the selected events is only of the order of 3% [85]. Furthermore, jets passing the electron identification cuts contain an important electromagnetic component and the detector response should therefore be similar to the response for true electrons. No additional systematic uncertainty is considered for the contamination of Z-jet events with background events.

2. As already mentioned, the uncertainty on the pile-up offset correction is treated as extra uncertainty (see Sect. 6.4)

3. The uncertainty induced by quark and gluon response differences as well as different quark and gluon compositions in data and in the simulation is addressed in Sect. 18.

In the final evaluation of systematic uncertainties, only effects that are significant with respect to their statistical uncertainties are taken into account [86]. The systematic effects and their statistical uncertainties are first evaluated using the initial binning. Then the results in neighbouring bins are iteratively combined until the observed effects become significant. The quadratic sum of all the components previously described is taken as the overall systematic uncertainty. Figure 19 summarises the different contributions to the total uncertainty, for EM+JES jets, in the whole $p_T$ range. For $R = 0.4$ jets and 25 GeV < $p_T^{\text{ref}}$ < 260 GeV, uncertainties are typically between 1% and 2%, and increase up to 10% for low transverse momenta.

## 9.6 Summary of the Z-jet analysis

The two ATLAS jet energy calibration schemes EM+JES and LCW+JES are probed using the direct $p_T$ balance between a central jet and a Z boson. The responses measured in the data and in the simulation are compared for jets defined by the anti-$k_t$ clustering algorithm with distance parameters of $R = 0.4$ and $R = 0.6$.

Figure 20 shows the data-to-MC ratio of the mean $p_T$ balance for jets calibrated with the EM+JES and the LCW+JES schemes, with statistical and systematic uncertainties. For $R = 0.4$ jets and $p_T^{\text{ref}} > 25$ GeV, this ratio is shifted by at most -4% from unity, and typically by -2% over most of the Z boson

$p_T$ range. Uncertainties are typically between 1% and 2% for $25 < p_T^{\text{ref}} < 260$ GeV, and increase up to 10% for low transverse momenta.

# 10 Jet energy calibration using $\gamma$-jet events

## 10.1 In situ jet calibration techniques

Two in situ techniques probing the calorimeter response to the jet balancing the photon are employed in this analysis:

### 1 Direct $p_T$ balance (DB)

The transverse momentum of the jet with the highest $p_T$ is compared to the transverse momentum of the reference photon ($p_T^\gamma$). The response is then computed as the ratio $p_T^{\text{jet}}/p_T^\gamma$.

### 2 Missing transverse momentum projection fraction (MPF)

The total hadronic recoil is used to estimate the calorimeter response to jets. The hadronic recoil is reconstructed from the vectorial sum of the transverse projections of the energy deposits in the calorimeter projected onto the photon direction. As in the direct $p_T$ balance, the photon $p_T$ serves as reference. The MPF response is defined as

$$ \mathcal{R}_{\text{MPF}} = 1 + \frac{\vec{p}_T^\gamma \cdot \vec{E}_T^{\text{miss}}}{|p_T^\gamma|^2}, $$

where the $\vec{E}_T^{\text{miss}}$ is computed with topo-clusters at the EM or LCW scales. A more detailed description of these two techniques can be found in Ref. [3].

Each technique has different sensitivities to additional soft-parton radiation, as well as to pile-up. The MPF is in general less sensitive to additional particle activity that is symmetric in the transverse plane, like for example pile-up and the underlying event.

The explicit use of jets in the jet response measurement from DB makes this technique clearly dependent on the jet reconstruction algorithm. Conversely, the dependence of the MPF technique on the jet algorithm is relegated to a second-order effect.<fnref n="10" /> Thus, in the following, when presenting the results from the MPF technique, no jet algorithm is in general explicitly mentioned.

## 10.2 Event selection of $\gamma$-jet events

The event selection used in this analysis is basically the same as that described in Ref. [3] for the 2010 analysis, except for changes that are either to adapt to the higher instantaneous luminosity of the 2011 dataset or to the different detector conditions. The event selection proceeds as follows:

1. Events are required to have a primary vertex, as defined in Sect. 5.4, with at least five associated tracks ($N_{\text{vertex}}^\text{tracks} \ge 5$).

<footnote n="10">¹⁰ Any dependence of the MPF response on the jet reconstruction algorithm is introduced solely by the event selection.</footnote>
<figure><img src="image_27.png" /><figcaption>Fig. 21: MPF response distributions in the γ-jet data for (a) $25 \le p_T^\gamma < 45$ GeV and (b) $160 \le p_T^\gamma < 210$ GeV when using topo-clusters at the EM scale. The dashed lines represent the fits with a Gaussian function. The mean value from the fit in each $p_T^\gamma$ bin is the value used as the measured average MPF response.</figcaption></figure>

<figure><img src="image_28.png" /><figcaption>Fig. 22: Jet response distributions in the γ-jet data for (a) $25 \le p_T^\gamma < 45$ GeV and (b) $160 \le p_T^\gamma < 210$ GeV as measured by the DB technique for anti-$k_t$ jets with $R = 0.6$ at the EM+JES scale. The dashed lines represent fits of Gaussian functions, except in the lowest bin ($25 \le p_T^\gamma < 45$ GeV), where the fit function is a Poisson distribution. The mean value from the fit in each $p_T^\gamma$ bin is the value used as the measured average jet response in DB.</figcaption></figure>
<figure><img src="image_29.png" /><figcaption>Fig. 23: Average jet response as determined by the MPF technique in γ-jet events using topo-clusters at the (a) EM and (b) LCW energy scales, for both data and MC simulations, as a function of the photon transverse momentum. The data-to-MC response ratio is shown in the bottom inset of each figure. Only the statistical uncertainties are shown.</figcaption></figure>

<table><caption>Table 3: Summary table of the criteria to select γ-jet events.</caption><thead><tr><th>Variable</th><th>Selection</th></tr></thead><tbody><tr><td>N<sup>vertex</sup><sub>tracks</sub></td><td>&gt; 4</td></tr><tr><td>p<sup>γ</sup><sub>T</sub></td><td>&gt; 25 GeV</td></tr><tr><td>|η<sup>γ</sup>|</td><td>&lt; 1.37</td></tr><tr><td>p<sup>jet</sup><sub>T</sub></td><td>&gt; 12 GeV</td></tr><tr><td>|η<sup>jet</sup>|</td><td>&lt; 1.2</td></tr><tr><td>E<sup>γ</sup><sub>T</sub> <sup>Iso</sup></td><td>&lt; 3 GeV</td></tr><tr><td>Δφ<sub>jet-γ</sub></td><td>&gt; 2.9 rad</td></tr><tr><td>p<sup>jet2</sup><sub>T</sub> / p<sup>γ</sup><sub>T</sub></td><td>&lt; 0.2 for DB (&lt; 0.3 for MPF)</td></tr></tbody></table>

<table><caption>Table 4: Table with the approximate number of selected events in each $p_T^γ$ bin.</caption><thead><tr><th>p<sup>γ</sup><sub>T</sub> [GeV]</th><th>Events</th><th>p<sup>γ</sup><sub>T</sub> [GeV]</th><th>Events</th></tr></thead><tbody><tr><td>25–45</td><td>20480</td><td>210–260</td><td>10210</td></tr><tr><td>45–65</td><td>61220</td><td>260–310</td><td>4650</td></tr><tr><td>65–85</td><td>125040</td><td>310–400</td><td>2770</td></tr><tr><td>85–110</td><td>262220</td><td>400–500</td><td>800</td></tr><tr><td>110–160</td><td>143180</td><td>500–600</td><td>240</td></tr><tr><td>160–210</td><td>32300</td><td>600–800</td><td>100</td></tr></tbody></table>

2. There must be at least one reconstructed photon; the highest-$p_T$ (leading) photon is taken as the hard-process photon and must have $p_T^γ > 25$ GeV.

3. The event is required to pass a single-photon trigger, with trigger $p_T$ threshold depending on the $p_T$ of the leading photon.

4. The leading photon must pass strict identification criteria [87], meaning that the pattern of energy deposition in the calorimeter is consistent with the expected photon showering behaviour.

5. The leading photon must lie in the pseudorapidity range $|η^γ| < 1.37$, meaning it is fully contained within the electromagnetic barrel calorimeter.

6. Jets with high electromagnetic content (e.g., jets fluctuating to a leading $π^0$, with $π^0 \to γγ$) may be misidentified as photons. In order to reduce this background, the leading photon is required to be isolated from other activity in the calorimeter. The isolation variable ($E_T^γ$) [87] is computed in a cone of size $R = 0.4$ around the photon, and corrected for pile-up energy inside the isolation cone. Only photons with $E_T^γ$ < 3 GeV are selected.

7. The photon reconstruction algorithm attempts to retain photons that have converted into an electron-positron pair. While clusters without matching tracks are directly classified as “unconverted” photon candidates, clusters matched to pairs of tracks originating from reconstructed conversion vertices are considered as “converted” photon candidates (double-track conversions). To increase the reconstruction efficiency of converted photons, conversion candidates where only one of the two tracks is reconstructed (single-track conversions) are also retained. Jets that are misidentified as
<figure><img src="image_30.png" /><figcaption>Fig. 24: Average jet response as determined by the DB technique in $\gamma$-jet events for anti-$k_t$ jets with ((a), (b)) $R = 0.4$ and ((c), (d)) $R = 0.6$, calibrated with the ((a), (c)) EM+JES scheme and with the ((b), (d)) LCW+JES scheme, for both data and MC simulations, as a function of the photon transverse momentum.</figcaption></figure>
photons fall more often in the category of converted pho-
tons, because fake photons produce wider showers and have
tracks associated to them. To suppress this background fur-
ther, the ratio of the transverse energy of the photon can-
didate cluster to the scalar sum of the $p_T$ of the matching
tracks ($E_T^\gamma \text{ cluster} / (\sum p_T^{\text{tracks}})$) is required to be in the range
from 0 to 2 for single-track conversions, and from 0.5 to
1.5 for double-track conversions. The fraction of converted
photons is ~30% throughout the $p_T^\gamma$ range under consider-
ation.

8. Only jets with $p_T^{\text{jet}} > 12$ GeV are considered. From those, only jets that pass quality criteria designed to reject fake jets originating from noise bursts in the calorimeters or from non-collision background or cosmics (see Sect. 5.3), are used. After these jet selections, each event is required to have at least one jet.

9. The highest-$p_T$ (leading) jet must be in the region $|\eta^{\text{jet}}| < 1.2$. This choice is motivated by the small $\eta$-intercalibration correction below 1.5% in this region.

10. To suppress soft radiation that would affect the $p_T$ balance between the jet and the photon, the following two conditions are required:

(a) The leading jet must be back-to-back to the photon in the transverse plane ($\Delta\phi(\text{jet}, \gamma) > 2.9$ rad).

(b) The $p_T$ of the sub-leading jet from the hard process ($p_T^{\text{jet2}}$) must be less than 20% (30%) of the $p_T$ of the photon for DB (MPF<fnref n="11" />). In order to distinguish jets from the hard process against jets from pile-up, the sub-leading jet is defined as the highest-$p_T$ jet from the subset of non-leading jets that either have JVF > 0.75 or for which JVF could not be computed because they are outside the region covered by the tracking system. See Sect. 8.2.3 for the explanation of JVF.

11. In the case of DB, the event is rejected if either the lea-
ding jet or the sub-leading jet falls in a region where, for
a certain period, the read-out of the EM calorimeter was
not functioning. For MPF, the condition is extended to all
jets with $p_T^{\text{jet}} > 20$ GeV in the event. A similar condition is
imposed on the reference photon.

A summary of the event selection criteria is given in Table 3.
Table 4 shows the approximate number of selected events per
$p_T^\gamma$ bin.

10.3 Jet response measurement

The calorimeter response to jets is measured in bins of the photon transverse momentum. Distributions of the MPF and the jet responses in the data are shown in Figs. 21 and 22, respectively, for $25 \le p_T^\gamma < 45$ GeV and for $160 \le p_T^\gamma < 210$ GeV. The distributions are fitted with a Gaussian function, except in the lowest $p_T^\gamma$ bin for DB where a Poisson distribution is used to address the issues introduced by the jet reconstruction $p_T$ threshold, as discussed in Sect. 9.3. The mean values from the fits define the average MPF and DB jet responses for each $p_T^\gamma$ bin. Figure 23 presents the results obtained in data and MC simulations for

MPF when the $E_T^{\text{miss}}$ is calculated from topo-clusters at the (a) EM and (b) LCW scales. Figure 24 shows the results for DB for anti-$k_t$ jets with radius parameter $R = 0.4$ and $R = 0.6$ for the EM+JES and LCW+JES calibration schemes.

With increasing jet energies, the particles inside the jet get more energetic as well. Higher incident energies for hadrons in non-compensating calorimeters, like the ones in ATLAS, increase the amount of energy invested in intrinsically induced electromagnetic showers, thus leading to an increase of the calorimeter response [88]. This increase is clearly observed for MPF, especially when topo-clusters at the EM scale are used as for the observations shown in Fig. 23(a). For DB, the effect is masked, because the jets used are already calibrated. DB is, in this case, measuring calibration residuals only.

Furthermore, a comparison of the MPF responses at EM scale in Fig. 23(a) and LCW scale in Fig. 23(b) shows the effect of having applied the LCW calibration to the topo-clusters. The response for jets built from LCW topo-clusters is much closer to unity, because the response differences between electromagnetic and hadronic particles in the jet are largely corrected by LCW at the level of the topo-clusters.

The lower part in Figs. 23 and 24 shows the ratio of the response in data to that in MC simulations. The MC simulation features a response that is 1% to 2% higher than that in data for $p_T^\gamma > 110$ GeV. For lower values of $p_T^\gamma$, the data-to-MC ratio tends to increase. Systematic studies have shown that the increase at low $p_T$ is due to the presence of contamination from multijet background events in the data, the different out-of-cone energy observed in data and in MC simulations, and the different effect of the 12 GeV jet $p_T$ reconstruction threshold (due to differences in the jet $p_T$ spectrum) on the response in data and in MC simulations.

10.4 Systematic uncertainties of photon–jet balance

The following sections briefly describe the procedure to es-
timate the systematic uncertainties of the $\gamma$-jet in situ tech-
niques. The dominant sources of systematic uncertainties, for
$p_T^\gamma \lesssim 75$ GeV, are the purity of the $\gamma$-jet data sample and for DB
also the out-of-cone correction (see Sect. 10.4.7) in the case of
$R = 0.4$ jets. For $p_T^\gamma \gtrsim 75$ GeV, the uncertainty on the energy
scale of the photon dominates.

10.4.1 Influence of pile-up interactions

The influence of in-time pile-up is evaluated by comparing the
response in events with six or more reconstructed primary ver-
tices ($N_{\text{PV}} \ge 6$) to the response in events with one or two recon-
structed primary vertices, inclusively in $\mu$. Similarly, the effect
of out-of-time pile-up is estimated comparing the response in
events with $\mu > 7$ to the response in events with $3.5 < \mu < 5.5$,
inclusively in $N_{\text{PV}}$. Since these two comparisons are highly
correlated, the pile-up uncertainty is estimated in each $p_T^\gamma$
bin as the maximum difference between the two high pile-up re-
sponses and the two low pile-up responses. For MPF, the un-
certainty due to pile-up is typically about $\sim 0.5\%$ or smaller.

In the case of DB however, the jet $p_T$ is already corrected
for the additional energy from pile-up interactions, as detailed

<footnote n="11">¹¹ For MPF, a less strict criterion can be used, since this technique is less sensitive to soft radiation.</footnote>
in Sect. 5. The variations in the data-to-MC response ratio obtained with the procedure explained above are found to be much smaller than other uncertainties on the measurement. They are also well contained within the variations obtained by propagating the uncertainty on the pile-up offset correction (see Sect. 16.2).

### 10.4.2 Soft-radiation suppression

The stability of the data-to-MC response ratio under soft radiation is evaluated in two steps. First, the cut on the $p_T$ of the sub-leading jet is varied, while keeping $\Delta\phi(\text{jet}, \gamma)$ fixed to its nominal cut value, and second, the cut on $\Delta\phi(\text{jet}, \gamma)$ is varied, with the cut on the sub-leading jet fixed to its nominal value. The cut on the sub-leading jet is varied to looser or tighter values as follows:

1. Tight:

$$p_T^{\text{jet2}} < \max\{10 \text{ GeV}, 0.2 \times p_T^{\gamma}\} \text{ for MPF, and} \\ p_T^{\text{jet2}} < \max\{10 \text{ GeV}, 0.1 \times p_T^{\gamma}\} \text{ for DB.}$$

2. Loose:

$$p_T^{\text{jet2}} < \max\{12 \text{ GeV}, 0.3 \times p_T^{\gamma}\} + 0.1 \times p_T^{\gamma} \text{ for MPF, and} \\ p_T^{\text{jet2}} < \max\{12 \text{ GeV}, 0.2 \times p_T^{\gamma}\} + 0.1 \times p_T^{\gamma} \text{ for DB.}$$

The typical variation on the data-to-MC response ratio is of the order of 0.5% for DB and smaller for MPF. Similar variations are observed when the $\Delta\phi(\text{jet}, \gamma)$ cut is relaxed to be $\Delta\phi(\text{jet}, \gamma) > 2.8$ or tightened to be $\Delta\phi(\text{jet}, \gamma) > 3.0$. Other tests of the stability of the data-to-MC response ratio under soft radiation are explored, such as relaxing and tightening the $\Delta\phi(\text{jet}, \gamma)$ and $p_T^{\text{jet2}}$ selection criteria at the same time, and lead to similar results.

### 10.4.3 Background from jet events

The uncertainty on the response due to the presence of jets that are identified as photons (fakes) in the data can be estimated, to first order, as $(1-P) \times (R_{\text{dijet}} - R_{\gamma-\text{jet}})/R_{\gamma-\text{jet}}$, where P is the purity of the $\gamma$-jet sample, and $R_{\gamma-\text{jet}}$ and $R_{\text{dijet}}$ are the responses in signal and background events, respectively.

The difference in response is estimated from MC simulations as in the 2010 analysis [3], using the nominal signal PYTHIA sample, and an inclusive jet PYTHIA sample (see Sect. 3) enriched in events with narrow jets, which are more likely to be misidentified as photons. The comparisons indicate that the relative response differences are below 5% for both techniques, which is taken as a conservative estimate. This is also confirmed by studying the response variation after relaxing the photon identification criterion.

The determination of the purity of the $\gamma$-jet data sample is done in the data using a sideband technique which is described in detail in Refs. [3, 87]. The purity is about 60% at $p_T^{\gamma} = 40$ GeV, rises with $p_T^{\gamma}$, and becomes greater than 95% for $p_T^{\gamma} \gtrsim 200$ GeV. This purity is lower than that measured in the 2010 analysis [3], due to the larger number of pile-up events

in the 2011 data. The effect of pile-up is tested by measuring the purity under the same high and low pile-up conditions used to estimate the uncertainty on the response due to pile-up (see Sect. 6). Variations in the purity of the order of 5% to 10% are found. The systematic uncertainty on the purity measurement is not taken into account in the estimation of the uncertainty due to background events, because it becomes negligible when multiplied by the relative response difference between the signal and background events.

The same purity estimate is used for MPF and DB, since both techniques have the same photon selection. The uncertainty due to background from jet events is ~2.5% at low $p_T^{\gamma}$, and decreases to about 0.1% towards high $p_T^{\gamma}$.

### 10.4.4 Photon energy scale

The electron energy is calibrated in situ using the measurements of the Z mass in $e^-e^+$ decays [73]. The main sources of the electron energy scale uncertainty are the energy loss in the interactions with the material in front of the calorimeter and the leakage of energy transversely to the topo-clusters axis. The calibration factors obtained from the $Z \to e^-e^+$ measurements are also applied to photons, with a corresponding increase in the systematic uncertainty (the difference between the electron and the photon energy scales is caused mainly by the different interaction of electrons and photons with the material in front of the calorimeter). The photon calibration and its uncertainty are propagated to the jet response measurement, leading in both techniques to an uncertainty of approximately +0.8% and -0.5%, independent of $p_T^{\gamma}$.

### 10.4.5 Jet energy resolution

The energy resolution for jets [82] in the MC simulation is very close to the resolution observed in data. The uncertainty on the jet energy resolution measurement in data is propagated as an uncertainty in the response in MC simulations. This is done as described in Sect. 8.4.6 and Eq. (10) therein. The observed difference in response between the varied and the nominal results is defined as the systematic uncertainty due to jet energy resolution.

### 10.4.6 Monte Carlo generator

Uncertainties due to different modelling of the parton shower, jet fragmentation and multiple parton interactions affecting the $p_T$ balance between the photon and the jet, can be estimated using different MC generators which implement different models. The jet response derived with PYTHIA is compared to the response derived using HERWIG. The results are shown in Fig. 25. The central value for the jet response in MC simulations is taken from PYTHIA, since this is the generator used to derive the JES corrections, and the observed full difference between PYTHIA and HERWIG is taken as a (symmetric) systematic uncertainty. The difference in the responses between HERWIG and PYTHIA is maximally about 1%.
<figure><img src="image_31.png" /><figcaption>Fig. 25: Average jet response as determined by the (a) MPF and (b) DB techniques, using anti-$k_t$ jets with $R = 0.4$ at the EM and EM+JES energy scales respectively, for PYTHIA (circles) and HERWIG (squares) MC simulations, as a function of the photon transverse momentum. The HERWIG-to-PYTHIA response ratio is shown in the bottom inset of each figure. Only the statistical uncertainties are shown.</figcaption></figure>

### 10.4.7 Out-of-cone radiation and underlying event

Even in a $2 \to 2$ $\gamma$-jet event, where the outgoing photon and parton (quark or gluon) perfectly balance each other in transverse momentum, the transverse momentum of the photon is only approximately equal to the transverse momentum of the truth jet, formed as described in Sect. 5.5, originating from the parton. The two main reasons for this are the same already described for the Z-jet events in Sect. 9.4, namely the fact that the jet does not capture all particles recoiling from the photon, and the contribution to the jet from the underlying event. The amount of momentum carried by particles outside the jet and by particles coming from soft interactions not contributing to the $p_T$ balance needs to be compared in data and MC simulation.

When averaging over many events, particles not associated to the hard scattering are distributed isotropically, and therefore they do not contribute to the hadronic recoil vector constructed in the MPF method. Thus, their contribution to the MPF response is zero. This is also supported by studies in the MC simulation using the particles produced by the underlying event model. Moreover, in the MPF technique the photon is balanced against the full hadronic recoil, not only against the leading jet. For the DB method the out-of-cone radiation is computed as explained in Sect. 9.4.

The measured $k_{OOC}$ factor (Eq. (12)) is shown as a function of $p_T^\gamma$ in Fig. 26 for anti-$k_t$ jets with $R = 0.4$ (Fig. 26(a)) and with $R = 0.6$ (Fig. 26(b)), for both data and MC simulations.

Systematic uncertainties obtained by varying the parameters in the $k_{OOC}$ factor definition are added in quadrature to the statistical uncertainties. The $k_{OOC}$ varies from 0.92 (0.97) at low $p_T$ to 0.99 (1.01) at high $p_T$ for $R = 0.4(0.6)$, respectively. The data are described by the MC simulation within 1% to 2% at low $p_T$. This deviation is taken as a systematic uncertainty in the DB technique.

### 10.4.8 Summary of systematic uncertainties

A summary of the systematic uncertainties for the MPF and the DB techniques as a function of the photon $p_T$ are presented in Figs. 27 and 28, respectively. The systematic uncertainties are shown for jets calibrated with the EM and LCW schemes for MPF, and with the EM+JES and LCW+JES schemes for DB where also jets with $R = 0.4$ and $R = 0.6$ are considered. The figures also show the statistical uncertainty, and the total uncertainty, which corresponds to the quadratic sum of all individual components (statistical and systematic). Table 5 shows the components of the systematic uncertainty for both methods in two representative $p_T^\gamma$ bins.

For the DB technique, the total uncertainty is as large as 2% to 3% at very low and very high $p_T$ values, and it is around 0.9% in the $p_T$ range from 100 GeV to 500 GeV. The uncertainties are smaller for MPF; the total uncertainty is ~0.7% in the range 100 GeV to 500 GeV and it is dominated by the photon energy scale uncertainty.
<figure><img src="image_32.png" /><figcaption>Fig. 26: Out-of-cone radiation factor $k_{OOC}$ relating the $p_T$ of the photon with the $p_T$ of the truth jet as a function of the photon transverse momentum, measured using charged particles, for anti-$k_t$ jets with (a) $R = 0.4$ and (b) $R = 0.6$, in data and in MC simulations. The data-to-MC response ratio is shown in the bottom inset of each plot. Statistical and systematic uncertainties are added in quadrature.</figcaption></figure>

<table><caption>Table 5: Systematic uncertainties on the data-to-MC ratio of the jet response on the EM scale for both DB and MPF in two representative $p_T^\gamma$ bins.</caption><thead><tr><th></th><th colspan="2">DB, R = 0.6 [%]</th><th colspan="2">MPF [%]</th></tr></thead><tbody><tr><td>p<sub>T</sub><sup>γ</sup> range [GeV]</td><td>45–65</td><td>310–400</td><td>45–65</td><td>310–400</td></tr><tr><td>Event</td><td></td><td></td><td></td><td></td></tr><tr><td>Pile-up</td><td>—</td><td>—</td><td>±0.21</td><td>±0.16</td></tr><tr><td>Radiation</td><td></td><td></td><td></td><td></td></tr><tr><td>p<sub>T</sub><sup>jet2</sup></td><td>±0.43</td><td>±0.28</td><td>±0.09</td><td>±0.10</td></tr><tr><td>Δφ (jet, γ)</td><td>±0.35</td><td>±0.20</td><td>±0.03</td><td>±0.03</td></tr><tr><td>Photon</td><td></td><td></td><td></td><td></td></tr><tr><td>Purity</td><td>±1.18</td><td>±0.15</td><td>±1.18</td><td>±0.15</td></tr><tr><td>Energy</td><td>±0.46</td><td>±0.71</td><td>±0.57</td><td>±0.61</td></tr><tr><td>Jet</td><td></td><td></td><td></td><td></td></tr><tr><td>JER</td><td>±0.01</td><td>±0.11</td><td>±0.04</td><td>±0.01</td></tr><tr><td>Out-of-cone</td><td>±0.60</td><td>±0.00</td><td>—</td><td>—</td></tr><tr><td>Modelling</td><td></td><td></td><td></td><td></td></tr><tr><td>MC generator</td><td>±0.48</td><td>±0.44</td><td>±0.38</td><td>±0.00</td></tr></tbody></table>

## 10.5 Summary of the γ-jet analysis

The average jet response in events with an isolated photon and a jet at high transverse momentum is computed using the 2011 dataset, and compared to the average jet response obtained using MC simulations. Two different techniques are used, the direct $p_T$ balance and the missing-$p_T$ projection fraction methods. Both techniques are highly correlated and show consistent

results within systematic uncertainties. The data-to-MC response ratio is close to 98% for $p_T^\gamma > 85$ GeV. Systematic uncertainties are evaluated for both methods to be of the order of 1% or smaller in most of the $p_T^\gamma$ range under consideration.

# 11 High-$p_T$ jet energy calibration using multijet events

## 11.1 Multijet balance technique and uncertainty propagation

The multijet balance (MJB) technique described in Ref. [3] can be used to verify the energy scale of jets and obtain correction factors that can correct for any non-linearity at very high $p_T$. The method exploits the $p_T$ balance in events where the highest-$p_T$ jet (leading jet) is produced back-to-back to a system composed of non-leading jets, referred to as a “recoil system”. The leading jet is required to have significantly larger $p_T$ than the jets in the recoil system in order to ensure that MJB is testing the absolute high-$p_T$ jet energy scale.

The vectorial sum of the $p_T$ of all non-leading jets defines the transverse momentum of the recoil system ($p_T^{recoil}$) that is expected to approximately balance the $p_T$ of the leading jet. The ratio

$$ \text{MJB} = \frac{|\vec{p}_T^{\text{leading}}|}{|\vec{p}_T^{\text{recoil}}|} $$

thus allows the verification of the JES of the leading jet using the properly calibrated non-leading jets at a lower $p_T$ scale. The asymmetry in the $p_T$ scale between the leading jet and non-leading jets is established by introducing a maximum limit on the ratio between the $p_T$ of the sub-leading (second-highest $p_T$) jet ($p_T^{jet2}$) and $p_T^{recoil}$. The calibration for the non-leading jets in
<figure><img src="image_33.png" /><figcaption>Fig. 27: Systematic uncertainties on the data-to-MC ratio of the jet response, as determined by the MPF technique for $\gamma$-jet events using topo-clusters at the (a) EM and (b) LCW energy scales, as a function of the photon transverse momentum.</figcaption></figure>

<figure><img src="image_34.png" /><figcaption>Fig. 28: Systematic uncertainties on the data-to-MC ratio of the jet response, as determined by the DB technique in $\gamma$-jet events, for anti-$k_t$ jets with ((a),(b)) $R = 0.4$ and ((c), (d)) $R = 0.6$, calibrated with the ((a), (c)) EM+JES scheme and with the ((b), (d)) LCW+JES scheme, as a function of the photon transverse momentum.</figcaption></figure>
<figure><img src="image_35.png" /><figcaption>Fig. 29: Multijet balance as a function of the recoil system $p_T^{\text{recoil}}$ for anti-$k_t$ jets with ((a), (b)) $R = 0.4$ and ((c), (d)) $R = 0.6$, calibrated with the ((a), (c)) EM+JES scheme and with the ((b), (d)) LCW+JES scheme, for both data and MC simulations. The non-leading jets in the data with $p_T < 750$ GeV are corrected by the combination of $\gamma$-jet and Z-jet in situ calibrations as described in Sect. 11.1. The open points in the bottom panel show the ratio of the MJB values between data and MC simulations. The curve in the same panel shows the data-to-MC ratio of the jet $p_T$ relative to the $p_T$ of a photon ($p_T^{\gamma}$) or a Z boson ($p_T^Z$) as a function of the $p_T^{\gamma}$ or $p_T^Z$ in $\gamma$-jet or Z-jet events, obtained in the combination mentioned above. Only the statistical uncertainties are shown.</figcaption></figure>

the recoil system is provided by the combination of the JES corrections derived from the $p_T$ balance in events with a jet and a Z boson (see Sect. 9) or a photon (see Sect. 10) for the absolute jet energy calibration, in addition to the $p_T$ balance in dijet events (see Sect. 8) for the relative ($\eta_{\text{det}}$ dependent) jet energy correction. See later Sect. 13.1 for detailed descriptions of the combination strategies in various $p_T$ ranges.

The MJB measured in data with the corrected non-leading jets (MJB<sup>Data</sup>) is compared with that in the simulation (MJB<sup>MC</sup>) to evaluate the JES calibration for the leading jet and assess the systematic uncertainty for high-$p_T$ jets. The statistical and systematic uncertainties of the $\gamma$-jet and Z-jet measurements are propagated through the combination. They are taken into account, together with the systematic uncertainty of the $\eta$-intercalibration, by fluctuating each sub-leading jet four momentum

within its uncertainties individually, and propagating those to higher $p_T$ as a variation in the MJB measurement. This whole procedure is repeated by increasing the sub-leading jet $p_T$ in steps, and applying the JES calibration derived in the previous step to the new event sample with harder non-leading jets. The MJB-based calibration is then calculated for the specific $p_T$ range and applied in the following increase of the sub-leading jet $p_T$. The procedure terminates once the number of events available for the next step becomes too low for a precise evaluation of MJB with the corresponding sample.

A cut on the ratio between $p_T^{\text{jet2}}$ and $p_T^{\text{recoil}}$, which defines the hard scale for the sub-leading jets, is also relaxed in the repetition sequences to effectively increase the statistics available in the calibration. The convolution of the propagated uncertainties from the JES calibrations applied to the non-leading jets
with systematic uncertainties associated with the MJB method itself, as described in Sect. 11.4, gives rise to a JES systematic uncertainty across the whole jet $p_T$ range accessible in 2011 data.

## 11.2 Selection of multijet events

In order to cover a wide $p_T$ range with enough event statistics, the analysis uses four single-jet triggers, each with a different jet-$p_T$ threshold. The highest $p_T$-threshold trigger that is active for the full dataset requires at least one jet with $p_T > 240$ GeV at the EM scale. The other three triggers are pre-scaled, i.e. only a defined fraction of them are recorded, and they require respective jet-$p_T$ thresholds of 55, 100, and 135 GeV. As shown below, the analysis is not limited by the statistical accuracy even with these pre-scaled jet triggers. In the offline analysis the data collected by a given trigger are used in non-overlapping $p_T^{\text{recoil}}$ ranges where the trigger is > 99% efficient.

Only events containing at least one primary vertex, defined as described in Sect. 5.4 and associated with at least five tracks, are considered. Events are rejected if they contain either an identified lepton (electron or muon) or a photon. Events are also rejected if they contain at least one jet which has $p_T > 20$ GeV that does not pass the jet cleaning criteria discussed in Sect. 5.3 to suppress noise or detector problems and mismeasured jets. For a certain period of time the read-out of a part of the EM calorimeter was not functioning, and events containing jets pointing to the affected region are also rejected. At the last stage of the event pre-selection, events are required to have at least three good-quality jets that have $p_T > 25$ GeV and $|\eta| < 2.8$. The leading jet is required to be within $|\eta| < 1.2$.

In order to select events having one jet produced against a well-defined recoil system, a selection is applied using two angular variables,

1. $\alpha = |\Delta\phi - \pi| < 0.3$ rad, where $\Delta\phi$ is the azimuthal opening angle between the highest-$p_T$ jet and the recoil system, and

2. The azimuthal opening angle between the leading jet and the non-leading jet that is closest in $\phi$ ($\beta$) is required to be $\beta > 1$ rad.

Two more selection criteria ensure that the sub-leading jets have a $p_T$ in the range where the in situ $\gamma$-jet and Z-jet calibrations are available and the leading jet is well above this range. The former is achieved by requiring the sub-leading jet $p_T^{\text{jet2}}$ to be less than 750 GeV and the latter by requiring that the ratio A between $p_T^{\text{jet2}}$ and $p_T^{\text{recoil}}$ satisfies $p_T^{\text{jet2}}/p_T^{\text{recoil}} < 0.6$. These two initial selections are modified when the analysis procedure is repeated as described above.

A summary of all cuts used in the analysis is given in Table 6.

## 11.3 Multijet balance measurement

The multijet balance obtained from the selected events for the EM+JES and LCW+JES calibrated jets with the anti-$k_t$ jet algorithm with $R = 0.4$ or $R = 0.6$ is shown in Fig. 29 for data and the MC simulations with PYTHIA.

<table><caption>Table 6: Summary of the event selection cuts used in the analysis. The first (second, third) values for $p_T^{\text{jet2}}$ and $p_T^{\text{jet2}}/p_T^{\text{recoil}}$ cuts are used in the first (second, subsequent) repetition of the MJB calibration procedure as described in Sec. 11.1.</caption><thead><tr><th>Variable</th><th>Cut value</th></tr></thead><tbody><tr><td>Jet p<sub>T</sub></td><td>&gt; 25 GeV</td></tr><tr><td>Jet rapidity</td><td>|η| &lt; 2.8</td></tr><tr><td>Leading jet rapidity</td><td>|η| &lt; 1.2</td></tr><tr><td>Number of good jets</td><td>≥ 3</td></tr><tr><td>p<sub>T</sub><sup>recoil</sup></td><td>&gt; 210 GeV</td></tr><tr><td>α</td><td>&lt; 0.3 rad</td></tr><tr><td>β</td><td>&gt; 1 rad</td></tr><tr><td>p<sub>T</sub><sup>jet2</sup></td><td>&lt; 750 (1200, 1450) GeV</td></tr><tr><td>p<sub>T</sub><sup>jet2</sup>/p<sub>T</sub><sup>recoil</sup></td><td>&lt; 0.6 (0.8, 0.8)</td></tr></tbody></table>

<table><caption>Table 7: Default values and the range of variation used to evaluate the systematic uncertainty on the analysis cuts.</caption><thead><tr><th>Variable</th><th>Default</th><th>Range</th></tr></thead><tbody><tr><td>Jet p<sub>T</sub></td><td>25 GeV</td><td>20–30 GeV</td></tr><tr><td>α</td><td>0.3 rad</td><td>0.1–0.4 rad</td></tr><tr><td>β</td><td>1.0 rad</td><td>0.50–1.50 rad</td></tr><tr><td>p<sub>T</sub><sup>jet2</sup>/p<sub>T</sub><sup>recoil</sup></td><td>0.6</td><td>0.4–0.7</td></tr></tbody></table>

The MJB decreases slightly at $p_T^{\text{recoil}}$ below 400 GeV, which is a consequence of the broadening of the $p_T^{\text{recoil}}$ distribution that can already be observed for jets formed from truth particles. The ratio between the distributions obtained from the data to the corresponding ones from MC simulations is shown in the lower part of each figure. It is compared with the data-to-MC ratio observed in the $\gamma$-jet and Z-jet in situ measurements. The agreement between data and MC simulations in the $p_T$ range covered by the $\gamma$-jet and Z-jet calibration, evaluated as the average value of the data-to-MC ratio, is within 2% (3%) for jets with $R = 0.4$ (0.6).

## 11.4 Systematic uncertainties on the multijet balance

Two main categories of systematic uncertainties are considered. The first category contains those which affect the reference $p_T$ of the recoil system. The second category includes those that affect the MJB variables used to probe the leading jet $p_T$, introduced mostly by effects from analysis cuts and imperfect MC modelling of the event.

The systematic uncertainty on the recoil system includes the following contributions:

### 1 Absolute JES uncertainty

The standard absolute JES uncertainties obtained from the combination of $\gamma$-jet and Z-jet techniques (see Sect. 13.1)
<figure><img src="image_36.png" /><figcaption>(a) $R = 0.4$, EM+JES</figcaption></figure>

<figure><img src="image_37.png" /><figcaption>(b) $R = 0.4$, LCW+JES</figcaption></figure>

<figure><img src="image_38.png" /><figcaption>(c) $R = 0.6$, EM+JES</figcaption></figure>

<figure><img src="image_39.png" /><figcaption>(d) $R = 0.6$, LCW+JES</figcaption></figure>

Fig. 30: Multijet balance with the nominal and varied Z-jet in situ calibrations as a function of the recoil system $p_T^{\text{recoil}}$ for anti-$k_t$ jets with ((a), (b)) $R = 0.4$ and ((c), (d)) $R = 0.6$, calibrated with the ((a), (c)) EM+JES scheme and with the ((b), (d)) LCW+JES scheme. The varied distributions are obtained by fluctuating the jet energy scale for the non-leading jets by $\pm 1\sigma$ for each of the systematic uncertainties for the Z-jet calibration and repeating the analysis over the data sample. The bottom panel shows the relative variations of the MJB with respect to the nominal case. The uppermost (lowermost) thick line in the bottom panel shows the total variation obtained by adding all the positive (negative) variations in quadrature. The colour coding used in the lower part of the figure is the same as that used in the upper one.

are included for each jet composing the recoil system. Figures 30 and 31 show the MJB variations obtained by scaling the non-leading jet energy and momentum scale by $\pm 1\sigma$ for each of the individual systematic uncertainties in the $\gamma$-jet and Z-jet calibrations, for the four jet calibration schemes. Each source of systematic uncertainty is described in Sect. 9.5 and Sect. 10.4, respectively. In case there are fewer than 10 events in a bin, the uncertainty is taken to be the RMS of the last bin with more than 10 events divided by the square root of the number of events in that bin. The central value of the ratio is unchanged.

This uncertainty ranges from 0.2% to 0.4% for Z-jet and 0.6% to 1.0% for $\gamma$-jet in the jet $p_T$ range of 0.5–1.2 TeV for the two jet sizes of $R = 0.4$ and 0.6.

## 2 Relative JES uncertainty

Relative jet response uncertainties evaluated in the dijet $\eta$-intercalibration (Sect. 8.4) are included in a similar manner for each jet with $|\eta| < 2.8$ in the recoil system.

## 3 Close-by jet uncertainty

The jet response is known to depend on the angular distance to the closest jet in $(\eta, \phi)$ space [3], and the response variation is expected to be more significant for jets belonging to the recoil system. Any discrepancy between MC simulations and data in describing the jet response with close-
<figure><img src="image_40.png" /><figcaption>(a) $R = 0.4$, EM+JES</figcaption></figure>

<figure><img src="image_41.png" /><figcaption>(b) $R = 0.4$, LCW+JES</figcaption></figure>

<figure><img src="image_42.png" /><figcaption>(c) $R = 0.6$, EM+JES</figcaption></figure>

<figure><img src="image_43.png" /><figcaption>(d) $R = 0.6$, LCW+JES</figcaption></figure>

Fig. 31: Multijet balance with the nominal and varied γ-jet in situ calibrations as a function of the recoil system $p_T^{\text{recoil}}$ for anti-$k_t$ jets with ((a), (b)) $R = 0.4$ and ((c), (d)) $R = 0.6$, calibrated with the ((a), (c)) EM+JES scheme and with the ((b), (d)) LCW+JES scheme. The varied distributions are obtained by fluctuating the jet energy scale for the non-leading jets by $\pm 1\sigma$ for each of the systematic uncertainties and repeating the analysis over the data sample. The bottom panel shows the relative variations of the MJB with respect to the nominal case. The uppermost (lowermost) thick line in the bottom panel shows the total variation obtained by adding all the positive (negative) variations in quadrature. The colour coding used in the lower part of the figure is the same as that used in the upper one.

by jets therefore results in an additional systematic uncertainty. The measurement performed to evaluate the effect and the resulting systematic uncertainty are described in Sect. 17. The close-by jet effect on MJB, shown in Fig. 32, is obtained by scaling the jet energy and momentum for each recoil jet using the results in Sect. 17.

The flavour composition of the jets could affect the agreement between MC simulations and data, and in principle cause an additional contribution to the JES uncertainty. Previous studies with 2010 data [3], however, show that the resulting uncertainty on MJB is less than 1%, and is therefore ignored in this evaluation of systematic uncertainties.

The jet response is corrected for energy deposited by additional proton-proton collisions in the same bunch crossings using the pile-up offset correction described in Sect. 6. The residual pile-up effect on MJB is checked by comparing the MJB values using sub-samples of data and MC simulations with different $N_{\text{PV}}$ and $\mu$ values. The result shows that the agreement between MC simulations and data is stable within its statistical uncertainty, and therefore an uncertainty due to pile-up is not considered.

The second systematic uncertainty category includes sources that affect the MJB variable which is used to probe the high-$p_T$ jet energy scale. As said earlier, those are mainly due to
<figure><img src="image_44.png" /><figcaption>Fig. 32: Relative uncertainties on the MJB due to the systematic uncertainty sources considered in the analysis as a function of the recoil system $p_T$ for anti-$k_t$ jets with ((a), (b)) $R = 0.4$ and ((c), (d)) $R = 0.6$, calibrated with the ((a), (c)) EM+JES scheme and with the ((b), (d)) LCW+JES scheme. The black line shows the total uncertainty obtained as a sum of all uncertainties in quadrature.</figcaption></figure>

<table><caption>Table 8: Representative values of systematic uncertainties in the $p_T^{\text{recoil}}$ range $500 \text{ GeV} < p_T^{\text{recoil}} < 1.2 \text{ TeV}$ for all effects considered in the analysis.</caption><thead><tr><th rowspan="2">Source</th><th colspan="2">EM+JES</th><th colspan="2">LCW+JES</th></tr><tr><th>R = 0.4</th><th>R = 0.6</th><th>R = 0.4</th><th>R = 0.6</th></tr></thead><tbody><tr><td>Jet size</td><td></td><td></td><td></td><td></td></tr><tr><td>Absolute JES</td><td>0.8%</td><td>0.7%</td><td>0.7%</td><td>0.7%</td></tr><tr><td>Relative JES</td><td>0.3%</td><td>0.4%</td><td>0.5%</td><td>0.4%</td></tr><tr><td>Close-by jet</td><td>0.6%</td><td>0.3%</td><td>0.6%</td><td>0.4%</td></tr><tr><td>Jet p<sub>T</sub> threshold</td><td></td><td></td><td>&lt; 0.4%</td><td></td></tr><tr><td>α cut</td><td></td><td></td><td>&lt; 0.1%</td><td></td></tr><tr><td>β cut</td><td></td><td></td><td>&lt; 0.2%</td><td></td></tr><tr><td>p<sub>T</sub><sup>jet2</sup>/p<sub>T</sub><sup>recoil</sup> cut</td><td>&lt; 0.1%</td><td>1.5%</td><td>&lt; 0.1%</td><td>1.2%</td></tr><tr><td>UE/radiation model</td><td></td><td></td><td>&lt; 0.5%</td><td></td></tr><tr><td>Fragmentation model</td><td>1.0%</td><td>0.3%</td><td>1.0%</td><td>0.5%</td></tr></tbody></table>

effects from analysis cuts or imperfect MC modelling with the following considerations:

## 1 Analysis cuts

A systematic uncertainty might be induced by event selection cuts on physical quantities that are not perfectly described by the MC simulation. In order to evaluate this systematic uncertainty, all relevant analysis cuts are varied in a range where the corresponding kinematic variables are not strongly biased and can be examined with small statistical fluctuations (see Table 7 for the range of variation). For each value of the cuts, the ratio of the value of MJB in data and simulation is evaluated. The maximum relative deviation of this ratio from the default value is taken as the systematic uncertainty from the source under consideration.

## 2 Jet rapidity acceptance

The analysis uses only jets with $|y| < 2.8$ in order to reduce the impact of the large JES uncertainties in the forward region. This selection, however, can cause additional systematic uncertainty because the fraction of jets produced
<figure><img src="image_45.png" /></figure>

Fig. 33: Multijet balance and systematic uncertainties related to the multijet balance technique and to the propagated uncertainties from the γ-jet and Z-jet balance as a function of the recoil system $p_T^{\text{recoil}}$ for anti-$k_t$ jets with ((a), (b)) $R = 0.4$ and ((c), (d)) $R = 0.6$, calibrated with the ((a), (c)) EM+JES scheme and with the ((b), (d)) LCW+JES scheme. The subleading jets in the data are corrected by the combination of γ-jet and Z-jet in situ calibrations at $p_T < 750$ GeV and MJB calibration at higher $p_T$ as described in Sect. 11.1. The three systematic uncertainty bands are obtained by adding individual systematic uncertainties for each calibration technique in quadrature. Also shown are predictions of the MJB from MC simulations for the three highest $p_T^{\text{recoil}}$-values, together with their systematic uncertainties propagated by using distribution from MC simulations. The bottom panel shows the relative variations of the MJB with respect to the nominal case.

outside the rapidity range can be different in the data and MC simulations, and hence affect the MJB values. This effect is checked, as is done in Ref. [3], by looking at the MJB for events with $p_T^{\text{recoil}} > 210$ GeV, as a function of the total transverse energy ($\Sigma E_T$) summed over all jets with $|y| > 2.8$. The majority of events have a very small $\Sigma E_T$ and the effect turns out to be negligible.

### 3 Underlying event, fragmentation and ISR/FSR modeling

Imperfect modelling of the UE, fragmentation and ISR/FSR may influence the multijet balance by affecting variables used to select events and kinematic properties of the leading jet and the recoil system. The systematic uncertainty for each of the mentioned sources is estimated by evaluating the data-to-MC ratio of the MJB, measured using

the default simulation based on PYTHIA and simulations using alternative MC generators. For the systematic uncertainty contribution from fragmentation, the HERWIG++ samples are used as an alternative. For the underlying event and radiation modelling systematics, the PYTHIA PERUGIA 2011 [89] samples are used. The systematic uncertainty introduced by these effects is 2 % or smaller in all cases except the lowest $p_T^{\text{recoil}}$ bins below 300 GeV.

All systematic uncertainties due to the analysis cuts and event modelling, and the total uncertainty obtained by summing them in quadrature, are shown as a function of $p_T^{\text{recoil}}$ in Fig. 32 for jets with $R = 0.4$ and $R = 0.6$, calibrated with the EM+JES and LCW+JES schemes. The uncertainties due to di-jet $\eta$-intercalibration and close-by jet effects are also included in the figure as well as the total uncertainty. Representative val-
<figure><img src="image_46.png" /><figcaption>Fig. 34: Relative jet response, $1/c$, as a function of the jet $\eta_{\text{det}}$ for anti-$k_t$ jets with $R = 0.4$ calibrated with the EM+JES scheme and in addition the derived $\eta$-intercalibration. Results are shown separately for (a) $55 \le p_T^{\text{avg}} < 75$ GeV and (b) $300 \le p_T^{\text{avg}} < 400$ GeV. For all points included in the original calibration ($|\eta_{\text{det}}| < 2.8$), the data are corrected to be consistent with the response in MC simulations using PYTHIA, as intended. The resulting calibration derived from the already calibrated data is shown as a thick line and is consistent with unity. The lower parts of the figures show the ratios between the relative jet response in data and MC.</figcaption></figure>

ues of the uncertainties in the $p_T^{\text{recoil}}$ range between 0.5 and 1.2 TeV are summarised in Table 8.

The summary of all systematic uncertainties associated with the multijet balance technique and the propagated uncertainties from the $\gamma$-jet and Z-jet in situ techniques overlaid on the data-to-MC ratio of the multijet balance, is shown in Fig. 33, for anti-$k_t$ jets with the distance parameters $R = 0.4$ and 0.6. The JES uncertainty is determined more precisely at jet $p_T$ below ~0.6 TeV by the $\gamma$-jet and Z-jet calibrations than the MJB calibration.

## 11.5 Summary of multijet analysis

The multijet balance technique is used to probe the jet energy scale in the TeV region for anti-$k_t$ jets with distance parameters $R = 0.4$ and $R = 0.6$. Exploiting the $p_T$ balance between the highest-$p_T$ jet and the recoil system composed of jets corrected by the $\gamma$-jet and Z-jet calibrations allows the extension of the in situ JES determination to higher $p_T$, beyond the range covered by the $\gamma$-jet calibration. Propagating systematic uncertainties associated with the $\gamma$-jet, Z-jet and dijet calibrations as well as the systematic uncertainty due to the knowledge of the recoil system transverse momentum in the MJB method (including the close-by jet uncertainty), the total systematic uncertainties for the $\gamma$-jet, Z-jet and MJB calibration methods are obtained to be about 0.6%, 0.3% and 1.5% respectively, for jets with $p_T = 1$ TeV. At high transverse momentum, the main contribution to the systematic uncertainty is due to the uncertainty on the MJB calibration. Considering the statistical uncertainty of the MJB calibration based on the 2011 data, the high-$p_T$ jet energy scale is validated at $p_T > 500$ GeV within 2.4% (2.0%) and 2.2% (3.0%) up to 1.2 TeV for anti-$k_t$ jets with $R = 0.4$ and $R = 0.6$, both calibrated with the EM+JES (LCW+JES) scheme.

# 12 Forward-jet energy measurement validation using Z-jet and $\gamma$-jet data

To test the performance of the forward-jet calibration derived in Sect. 8, this calibration is applied to all jets in the original dataset and the full analysis is repeated. The resulting intercalibration results are within 0.3% of unity across the full ($p_T^{\text{avg}}, \eta_{\text{det}}$) phase space in which the calibration is derived, both for jets with $R = 0.4$ and $R = 0.6$, and for the EM+JES and LCW+JES calibrations. The measured relative response for two representative bins of $p_T^{\text{avg}}$ is shown in Fig. 34.

Similar to the analyses described in Sect. 9 and Sect. 10, the balance between a Z boson decaying to an electron-positron pair and a recoiling forward jet, and the balance between a photon and a forward jet, are used to study the jet response in the forward direction. The results for Z-jet and $\gamma$-jet, as shown in Fig. 35, agree with the calibrations and uncertainty derived from the dijet analysis.

The Z-jet study also includes predictions from the ALPGEN generator, which uses HERWIG for parton shower and fragmentation into particles (see Sect. 3 for generator configuration details). The ALPGEN+HERWIG response predictions generally agree with the expectations within the modelling uncertainty of this analysis (see Sect. 8.4.1). The $\gamma$-jet results include comparisons with PYTHIA events, generated with the same tune and version as the PYTHIA dijet samples used in this analysis, and a sample produced with HERWIG, using the already mentioned ATLAS AUET2B MRST LO** tune and the MRST LO** PDF set (see Sect. 3).
<figure><img src="image_47.png" /><figcaption>(a) $25 \le p_T^{\text{ref}} < 35$ GeV</figcaption></figure>

<figure><img src="image_48.png" /><figcaption>(b) $50 \le p_T^{\text{ref}} < 80$ GeV</figcaption></figure>

<figure><img src="image_49.png" /><figcaption>(c) $85 \le p_T^{\gamma} < 110$ GeV</figcaption></figure>

<figure><img src="image_50.png" /><figcaption>(d) $210 \le p_T^{\gamma} < 260$ GeV</figcaption></figure>

Fig. 35: The ((a), (b)) Z-jet and ((c), (d)) $\gamma$-jet balance for anti-$k_t$ jets with $R = 0.4$, calibrated with the EM+JES scheme. In (a) and (b) the results for events with $25 \le p_T^{\text{ref}} < 35$ GeV and $50 \le p_T^{\text{ref}} < 80$ GeV are shown, respectively. Here $p_T^{\text{ref}}$ is the $p_T$ of the reconstructed Z boson projected onto the axis of the balancing jet. The $p_T$ balance for $\gamma$-jet events with photons with transverse momenta $p_T^{\gamma}$ within $85 \le p_T^{\gamma} < 100$ GeV is shown in (c), while (d) shows the $p_T$ balance for higher photon transverse momenta ($210 \le p_T^{\gamma} < 260$ GeV). As no in situ calibration is applied to these measurements, it is expected that data and MC simulations using PYTHIA are shifted relative to each other by the absolute correction multiplied by the relative ($\eta_{\text{det}}$ dependent) correction presented herein. The resulting JES calibration is shown as a solid line in the lower part of the figures. The dijet modelling uncertainty is shown as a filled band around the in situ correction.

<table><caption>Table 9: Summary of the number of events available for various in situ techniques after all selection cuts. The numbers are given for illustration in specific $p_T^{\text{jet}}$ ranges for anti-$k_t$ jets with $R = 0.4$ reconstructed with the EM+JES scheme. The $\gamma$-jet results are based on the MPF method.</caption><tbody><tr><td>Z-jet method</td><td></td><td></td><td></td></tr><tr><td>$p_T^{\text{jet}}$</td><td>20-25 GeV</td><td>35-45 GeV</td><td>210-260 GeV</td></tr><tr><td>Number of events</td><td>8530</td><td>8640</td><td>309</td></tr><tr><td>γ-jet method</td><td></td><td></td><td></td></tr><tr><td>$p_T^{\text{jet}}$</td><td>25-45 GeV</td><td>45-65 GeV</td><td>210-260 GeV</td><td>600-800 GeV</td></tr><tr><td>Number of events</td><td>20480</td><td>61220</td><td>10210</td><td>100</td></tr><tr><td>multijet method</td><td></td><td></td><td></td></tr><tr><td>$p_T^{\text{jet}}$</td><td></td><td>210-260 GeV</td><td>750-950 GeV</td><td>1.45-1.8 TeV</td></tr><tr><td>Number of events</td><td></td><td>2638</td><td>3965</td><td>48</td></tr></tbody></table>
<figure><img src="image_51.png" /><figcaption>Fig. 36: Weight carried by each in situ technique in the combination to derive the residual jet energy scale calibration as a function of the jet transverse momentum $p_T^{jet}$ for anti-$k_t$ jets with $R = 0.4$ calibrated with the (a) EM+JES and the (b) LCW+JES scheme. The $p_T^{jet}$ dependence of the weights is discussed in Sect. 13.4.</figcaption></figure>

# 13 Jet energy calibration and uncertainty combination

## 13.1 Overview of the combined JES calibration procedure

After the first JES calibration step described in Sect. 5, the jet transverse momenta $p_T^{jet}$ in data and MC simulation are compared using in situ techniques that exploit the balance<fnref n="12" /> between $p_T^{jet}$ and the $p_T$ of a reference object ($p_T^{ref}$):

$$ \mathcal{R}(p_T^{jet}, \eta) = \frac{\langle p_T^{jet} / p_T^{ref} \rangle_{\text{data}}}{\langle p_T^{jet} / p_T^{ref} \rangle_{\text{MC}}} \quad (13) $$

The inverse of this quantity is the residual JES correction factor for jets measured in data, and thus reflects the final JES calibration in ATLAS. It is derived from corrections individually described in Sect. 7. The sequence of these corrections is briefly summarised again below, with references to the corresponding more detailed descriptions:

1. Apply $\eta$-intercalibration to remove the $\eta_{det}$ dependence of the detector response to jets within $0.8 \le |\eta| < 4.5$ by equalising it with the one for jets within $|\eta_{det}| < 0.8$ (see Sect. 7.1).

2. Apply the absolute correction, as derived using a combination of the Z-jet (Sect. 9) and the $\gamma$-jet (Sect. 10) methods, to the the central jet response ($|\eta_{det}| < 1.2$). The slightly

larger $\eta_{det}$ range used here, compared to the one used in $\eta$-intercalibration, provides more statistics while keeping systematic uncertainties small. The corresponding combined JES uncertainty is determined from the uncertainties of each of these techniques, as presented in detail in Sect. 13.3. The absolute scale correction, together with its systematic uncertainties, is also evaluated for jets in the end-cap and forward detector region ($|\eta_{det}| \ge 1.2$), and accordingly applied to those as well.

3. Jets with energies in the TeV regime are calibrated using the multijet transverse momentum balance technique (MJB in Sect. 11). The lower-$p_T$ jets are within $|\eta_{det}| < 2.8$, while the leading jet is required to be within $|\eta_{det}| < 1.2$. The uncertainties derived from $\gamma$-jet, Z-jet and dijet $p_T$ balance for the lower-$p_T$ jets are propagated to the higher-$p_T$ jets (Sect. 11.4).

The in situ JES calibration and the corresponding JES uncertainty for central jets ($|\eta_{det}| < 1.2$) are hence derived by a combination of the data-to-MC ratios $\mathcal{R}$, individually determined as given in Eq. (13), obtained from the $\gamma$-jet, Z-jet and MJB correction methods. The JES uncertainties for forward jets $1.2 < |\eta_{det}| < 4.5$ are then derived from those for central jets using the dijet $\eta$-intercalibration technique.

Table 9 summarises the number of events available for each correction method in various kinematic bins. Details of the combination method, including the full evaluation of the systematic uncertainties and its underlying components (nuisance parameters), are further explained in the remainder of this section.

## 13.2 Combination technique

The data-to-MC response ratios (see Eq. 13) of the various in situ methods are combined using the procedure described in Ref. [3]. The in situ jet response measurements are made in

<footnote n="12">¹² As for all $p_T$ balance evaluations between a reference and a probe object, the expectation value of this balance is not unity, due to physics effects (e.g., ISR) and jet reconstruction inefficiencies (e.g., out-of-cone energy losses). The ability of the MC simulation to reproduce all of these effects is further discussed in the context of the evaluation of the systematic uncertainties in Sect. 13.3.</footnote>
<figure><img src="image_52.png" /><figcaption>(a) EM+JES Z-jet DB</figcaption></figure>

<figure><img src="image_53.png" /><figcaption>(b) LCW+JES Z-jet DB</figcaption></figure>

<figure><img src="image_54.png" /><figcaption>(c) EM+JES γ-jet MPF</figcaption></figure>

<figure><img src="image_55.png" /><figcaption>(d) LCW+JES γ-jet MPF</figcaption></figure>

<figure><img src="image_56.png" /><figcaption>(e) EM+JES Multijet</figcaption></figure>

<figure><img src="image_57.png" /><figcaption>(f) LCW+JES Multijet</figcaption></figure>

Fig. 37: Individual uncertainty sources applicable to the combined response ratio as a function of the jet $p_T$ for the three in situ techniques: (a), (b)) Z-jet direct balance, (c), (d)) γ-jet MPF and (e), (f)) multijet balance for anti-$k_t$ jets with $R = 0.4$ calibrated with the (a), (c), (e)) EM+JES and the (b), (d), (f)) LCW+JES scheme. The systematic uncertainties displayed here correspond to the components listed in Table 10.
<figure><img src="image_58.png" /></figure>

Fig. 38: Ratio of the average jet response $\langle p_T^{\text{jet}}/p_T^{\text{ref}} \rangle$ measured in data to that measured in MC simulations for jets within $|\eta| < 1.2$ as a function of the transverse jet momentum $p_T^{\text{jet}}$. The data-to-MC jet response ratios are shown separately for the three in situ techniques used in the combined calibration: direct balance in Z-jet events, MPF in $\gamma$-jet events, and multijet $p_T$ balance in inclusive jet events. The error bars indicate the statistical and the total uncertainties (adding in quadrature statistical and systematic uncertainties). Results are shown for anti-$k_t$ jets with $R = 0.4$ calibrated with the (a) EM+JES and the (b) LCW+JES scheme. The light band indicates the total uncertainty from the combination of the in situ techniques. The inner dark band indicates the statistical component only.

bins of $p_T^{\text{ref}}$ and within $|\eta_{\text{det}}| < 1.2$, and are evaluated at the barycentre $\langle p_T^{\text{ref}} \rangle$ of each $p_T^{\text{ref}}$ bin, for each $\eta_{\text{det}}$ range.<fnref n="13" />

First, a common, fine $p_T$ binning is introduced for the combination of methods. In each of these $p_T$ bins, and for each in situ method that contributes to that bin, the data-to-MC response ratio is determined using interpolating splines based on second-order polynomials. The combined data-to-MC ratio $R_{\text{extrap}}(\langle p_T^{\text{jet}} \rangle, \eta_{\text{det}})$ is then determined by the weighted average of the interpolated contributions from the various methods. The weights are obtained by a $\chi^2$ minimisation of the response ratios in each $p_T$ bin, and are therefore proportional to the inverse of the square of the uncertainties of the input measurements. The local $\chi^2$ is also used to test the level of agreement between the in situ methods.

Each uncertainty source of the in situ methods is treated as fully correlated across $p_T$ and $\eta_{\text{det}}$, while the individual uncertainty sources inside a method and between the methods are assumed to be independent of each other. The full set of uncertainties is propagated from the in situ methods to the combined result in each $p_T$ bin using pseudo-experiments [3]. For some applications like the combination and comparison of several experimental measurements using jets, it is necessary to understand the contribution of each uncertainty component to the final total uncertainty. For this purpose, each uncertainty component is propagated separately from each in situ method to the combined result. This is achieved by coherently shifting all the correction factors obtained by the in situ methods by one

standard deviation of a given uncertainty component, and re-doing the combination using the same set of averaging weights as in the nominal combination. The comparison of the shifted average correction factors with the nominal ones provides the propagated systematic uncertainty.

To account for potential disagreement between in situ measurements constraining the same term (referred to as measurements which are *in tension*), each uncertainty source is rescaled by the factor $\sqrt{\chi^2/\text{dof}}$, if this factor is larger than 1. This is conservative, as values of $\sqrt{\chi^2/\text{dof}}$ larger than 1 can also be reached due to statistical fluctuations.

$R_{\text{extrap}}(\langle p_T^{\text{jet}} \rangle, \eta_{\text{det}}) = 1/c$ is used as the in situ correction calibration factor and its inverse $c$ is applied to data. The correction factor still contains part of the statistical fluctuations of the in situ measurements. The influence of the statistical fluctuations is reduced by applying a minimal amount of smoothing using a sliding Gaussian kernel to the combined correction factors [3].

Each uncertainty component from the in situ methods is also propagated through the smoothing procedure. Propagating information between close-by $p_T$ regions, the smoothing procedure changes the amplitude of the uncertainties (e.g. reducing them at low $p_T$).

## 13.3 Uncertainty sources of the in situ calibration techniques

The in situ techniques usually rely on assumptions that are only approximately fulfilled. One example is the assumption that the calibrated jet and the reference object are balanced in

<footnote n="13">¹³ Since $\langle p_T^{\text{jet}}/p_T^{\text{ref}} \rangle$ is close to unity for all $p_T^{\text{ref}}$ bins, the bin barycentre $\langle p_T^{\text{ref}} \rangle$ is a good approximation of $\langle p_T^{\text{jet}} \rangle$. In the following $p_T^{\text{jet}}$ is used.</footnote>
<table><caption>Table 10: Summary table of the uncertainty components for each in situ technique (Z-jet (see Section 9), γ-jet (see Section 10), and multijet $p_T$ balance (see Section 11)) used to derive the jet energy scale uncertainty. Shown are the 21 systematic uncertainty components together with the 11, 12 and 10 statistical uncertainty components for each in situ technique. Each uncertainty component is categorised depending on its source as either detector (DETECTOR), physics modelling (MODEL), mixed detector and modelling (MIXED), or as statistics and method (STAT/METH).</caption><thead><tr><th>Name</th><th>Description</th><th>Category</th></tr></thead><tbody><tr><td><strong>Common sources</strong></td><td></td><td></td></tr><tr><td>Electron/photon E scale</td><td>electron or photon energy scale</td><td>DETECTOR</td></tr><tr><td><strong>DB Z-jet p<sub>T</sub> balance</strong></td><td></td><td></td></tr><tr><td>MC generator</td><td>MC generator difference between ALPGEN/HERWIG and PYTHIA</td><td>MODEL</td></tr><tr><td>Radiation suppression</td><td>radiation suppression due to second jet cut</td><td>MODEL</td></tr><tr><td>Extrapolation</td><td>extrapolation in Δφ<sub>jet-Z</sub> between jet and Z boson</td><td>MODEL</td></tr><tr><td>Pile-up jet rejection</td><td>jet selection using jet vertex fraction</td><td>MIXED</td></tr><tr><td>Out-of-cone</td><td>contribution of particles outside the jet cone</td><td>MODEL</td></tr><tr><td>Width</td><td>width variation in Poisson fits to determine jet response</td><td>STAT/METH</td></tr><tr><td>Statistical components</td><td>statistical uncertainty for each of the 11 bins</td><td>STAT/METH</td></tr><tr><td><strong>MPF γ-jet p<sub>T</sub> balance (MPF)</strong></td><td></td><td></td></tr><tr><td>MC generator</td><td>MC generator difference HERWIG and PYTHIA</td><td>MODEL</td></tr><tr><td>Radiation suppression</td><td>radiation suppression due to second jet cut</td><td>MODEL</td></tr><tr><td>Jet resolution</td><td>variation of jet resolution within uncertainty</td><td>DETECTOR</td></tr><tr><td>Photon Purity</td><td>background response uncertainty and photon purity estimation</td><td>DETECTOR</td></tr><tr><td>Pile-up</td><td>sensitivity to pile-up interactions</td><td>MIXED</td></tr><tr><td>Out-of-cone</td><td>contribution of particles outside the jet cone</td><td>MODEL</td></tr><tr><td>Statistical components</td><td>statistical uncertainty for each of the 12 bins</td><td>STAT/METH</td></tr><tr><td><strong>MJB Multijet p<sub>T</sub> balance</strong></td><td></td><td></td></tr><tr><td>α selection</td><td>angle between leading jet and recoil system</td><td>MODEL</td></tr><tr><td>β selection</td><td>angle between leading jet and closest sub-leading jet</td><td>MODEL</td></tr><tr><td>Dijet balance</td><td>dijet balance correction applied for |η| < 2.8</td><td>MIXED</td></tr><tr><td>Close-by, recoil</td><td>JES uncertainty due to close-by jets in the recoil system</td><td>MIXED</td></tr><tr><td>Fragmentation</td><td>jet fragmentation modelling uncertainty</td><td>MIXED</td></tr><tr><td>Jet p<sub>T</sub> threshold</td><td>jet p<sub>T</sub> threshold</td><td>MIXED</td></tr><tr><td>p<sub>T</sub> asymmetry selection</td><td>p<sub>T</sub> asymmetry selection between leading jet and sub-leading jet</td><td>MODEL</td></tr><tr><td>UE,ISR/FSR</td><td>soft physics effects modelling: underlying event and soft radiation</td><td>MIXED</td></tr><tr><td>Statistical components</td><td>statistical uncertainty for each of the 10 bins</td><td>STAT/METH</td></tr></tbody></table>

transverse momentum, while this balance can be altered by the presence of additional high-$p_T$ particles. In order to determine the JES uncertainties, the modelling of physics effects has to be disentangled from detector effects. These effects can be studied by looking at the changes of the data-to-MC response ratios introduced by systematic variation of the event selection criteria. The ability of the MC simulation to describe these changes under large variations of the selection criteria determines the systematic uncertainty in the in situ methods, since physics effects can be suppressed or amplified by these variations. In addition, systematic uncertainties related to the selection, calibration and modelling of the reference object need to be considered.

When performing the variations of the selection criteria, only statistically significant variations of the response ratios are propagated to the systematic uncertainties. This is achieved by evaluating the systematic uncertainties in intervals which

can be larger than the bins used for the measurement of the response ratios, meaning that several bins are iteratively combined until the observed deviations are significant. By doing so, one avoids multiple counting of the statistical uncertainties in the systematics that are evaluated. Using this approach, it is found that the radiation suppression uncertainty for the $\Delta\phi(\text{jet}, \gamma)$ cut on the MPF method (see Sect. 10.4.2) can be dropped.<fnref n="14" />

For the relative $\eta$-intercalibration described in Sect. 7.1 the dominant uncertainty source is due to MC modelling of jets at forward rapidities, where properties differ significantly for the generators under consideration (PYTHIA and HERWIG). Other

<footnote n="14">¹⁴ This uncertainty is very small, and the corresponding variations are not significant, even when the evaluation is performed on the full $p_T$ range.</footnote>
systematic uncertainty sources arise due to the modelling of the jet resolution, the trigger, and dijet topology selection. However, these components are negligible when compared to the MC modelling uncertainty.

The data-to-MC response ratio given in Eq. (13) for the direct balance in Z-jet events, the MPF technique in γ-jet events, and the multijet balance method are combined as described in the previous Sect. 13.2. In this combination, the ability of the MC simulation to describe the data, the individual uncertainties of the in situ techniques and their compatibility, are considered. The uncertainties of the three central in situ methods combined here are described by a set of 54 systematic uncertainty sources listed in Table 10. The photon and electron energy scale uncertainties are treated as being fully correlated at this level. Components directly related to the dijet balance technique are η dependent quantities, and are thus treated differently. Such parameters are not included in the list of the 54 components, although uncertainties related to their propagation through other methods are included. In Table 10, each uncertainty component is assigned to one of four categories, based on its source and correlations:

1. Detector description (DETECTOR)

2. Physics modelling (MODEL)

3. Statistics and method (STAT/METH)

4. Mixed detector and modelling (MIXED).

The motivation for these categories, and to some extend the guidance for assigning the 54 individual components to them, are given by considerations concerning the comparability of jet measurements and their uncertainties in different experiments. For example, the DETECTOR and STAT/METH categories can be considered largely uncorrelated between experiments, while the MODEL category is likely correlated.

<figure><img src="image_59.png" /><figcaption>Fig. 39: Difference between the data-to-MC response ratio $R$ measured using the direct balance (DB) and the missing momentum fraction (MPF) methods for jets reconstructed with the anti-$k_t$ algorithm with $R = 0.4$ calibrated with the EM+JES and LCW+JES schemes. The error bars shown only contain the uncorrelated uncertainties.</figcaption></figure>

## 13.4 Combination results

Figure 36 shows the contribution of each in situ technique to the JES residual calibration, defined to be the fractional weight carried in the combination. In the region $p_T^{\text{jet}} \lesssim 100$ GeV, the Z-jet method has the highest contribution to the overall JES average. The contribution is 100% for $p_T^{\text{jet}}$ below 25 GeV, the region covered only by Z-jet, about 90% at $p_T^{\text{jet}} = 40$ GeV, and decreases to about 50% at $p_T^{\text{jet}} = 100$ GeV. In order to prevent the uncertainties specific to the low-$p_T^{\text{jet}}$ region from propagating to higher $p_T^{\text{jet}}$ in the combination, the Z-jet measurements below and above $p_T^{\text{jet}} = 25$ GeV are treated separately, meaning no interpolation is performed across $p_T^{\text{jet}} = 25$ GeV, although the magnitude of the original systematic uncertainty sources is used, separately, in both regions.

The weaker correlations between the uncertainties of the Z-jet measurements, compared to ones from γ-jet, lead to a faster increase of the extrapolated uncertainties, hence to the reduction of the Z-jet weight in the region between 25 and 40 GeV. In the region $100 \lesssim p_T^{\text{jet}} \lesssim 600$ GeV, the γ-jet method dominates with a weight increasing from 50% at $p_T^{\text{jet}} = 100$ GeV to about 80% at $p_T^{\text{jet}} = 500$ GeV. For $p_T^{\text{jet}} \gtrsim 600$ GeV the measurement based on multijet balance becomes increasingly important and for $p_T^{\text{jet}} \gtrsim 800$ GeV it is the only method contributing to the JES residual calibration. The combination results and the relative uncertainties are considered in the $p_T$ range from 17.5 GeV to 1 TeV, where sufficient statistics are available.

The individual uncertainty components for the final combination results,<fnref n="15" /> are shown in Fig. 37 for anti-$k_t$ jets with $R = 0.4$ for the EM+JES and the LCW+JES calibration scheme and for each in situ technique.

The agreement between the in situ methods is good, with $\chi^2/\text{dof} < 1$ for most $p_T$ bins, and values up to $\chi^2/\text{dof} = 1.5$ in only a few bins. The largest $\chi^2/\text{dof} = 2$ is found for anti-$k_t$ jets with $R = 0.6$ calibrated with the LCW+JES scheme for $p_T^{\text{jet}} = 25$ GeV.

The final JES residual calibration obtained from the combination of the in situ techniques is shown in Fig. 38, together with statistical and systematic uncertainties. A general offset of about -2% is observed in the data-to-MC response ratios for jet transverse momenta below 100 GeV. The offset decreases to about -1% at higher $p_T$ ($p_T^{\text{jet}} \gtrsim 200$). The JES uncertainty from the combination of the in situ techniques is about 2.5% at $p_T^{\text{jet}} = 25$ GeV, and decreases to below 1% for $55 \le p_T^{\text{jet}} < 500$ GeV. The multijet balance method is used up to 1 TeV, as at higher $p_T$ values it has large statistical uncertainties. At 1 TeV the total uncertainty is about 1.5%.

The results for the EM+JES and the LCW+JES calibration schemes for jets with $R = 0.6$ are similar to those for $R = 0.4$.

<footnote n="15">¹⁵ The uncertainties apply to the overall result of the combination of the in situ techniques and differ from the original uncertainties of the in situ methods, as they are convoluted with the corresponding weights.</footnote>