1 Introduction

The missing transverse momentum (\(E_{\mathrm{T}}^{\mathrm{miss}}\)) is an important observable serving as an experimental proxy for the transverse momentum carried by undetected particles produced in proton–proton (\( pp \)) collisions measured with the ATLAS detector [1] at the Large Hadron Collider (LHC). It is reconstructed from the signals of detected particles in the final state. A value incompatible with zero may indicate not only the production of Standard Model (SM) neutrinos but also the production of new particles suggested in models for physics beyond the SM that escape the ATLAS detector without being detected. The reconstruction of \(E_{\mathrm{T}}^{\mathrm{miss}}\) is challenging because it involves all detector subsystems and requires the most complete and unambiguous representation of the hard interaction of interest by calorimeter and tracking signals. This representation is obscured by limitations introduced by the detector acceptance and by signals and signal remnants from additional \( pp \) interactions occurring in the same, previous and subsequent LHC bunch crossings (pile-up) relative to the triggered hard-scattering. ATLAS has developed successful strategies for a high-quality \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction focussing on the minimisation of effects introduced by pile-up for the data recorded between 2010 and 2012 (LHC \(\text {Run}\,1\)) [2, 3]. These approaches are the basis for the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction developed for the data collected in 2015 (LHC \(\text {Run}\,2\)) that is described in this paper, together with results from performance evaluations and the determination of systematic uncertainties.

This paper is organised as follows. The subsystems forming the ATLAS detector are described in Sect. 2. The \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction is discussed in Sect. 3. The extraction of the data samples and the generation of the Monte Carlo (MC) simulation samples are presented in Sect. 4. The event selection is outlined in Sect. 5, followed by results for \(E_{\mathrm{T}}^{\mathrm{miss}}\) performance in Sect. 6. Section 7 comprises a discussion of methods used to determine systematic uncertainties associated with the \(E_{\mathrm{T}}^{\mathrm{miss}}\) measurement, and the presentation of the corresponding results. Section 8 describes variations of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction using calorimeter signals for the soft hadronic event activity, or reconstructed charged-particle tracks only. The paper concludes with a summary and outlook in Sect. 9. The nomenclature and conventions used by ATLAS for \(E_{\mathrm{T}}^{\mathrm{miss}}\)-related variables and descriptors can be found in Appendix A, while the composition of \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction variants is presented in Appendix B. An evaluation of the effect of alternative jet selections on the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction performance is given in Appendix C.

2 ATLAS detector

The ATLAS experiment at the LHC features a multi-purpose particle detector with a forward–backward symmetric cylindrical geometry and a nearly full (\(4\pi \)) coverage in solid angle.Footnote 1 It consists of an inner detector (ID) tracking system in a \(2\,\text {T}\) axial magnetic field provided by a superconducting solenoid. The solenoid is surrounded by electromagnetic and hadronic calorimeters, and a muon spectrometer (MS). The ID covers the pseudorapidity range \(|\eta | < 2.5\), and consists of a silicon pixel detector, a silicon microstrip detector and a transition radiation tracker for \(|\eta | < 2.0\). During the LHC shutdown between \(\text {Run}\,1\) and \(\text {Run}\,2\), a new tracking layer, known as the insertable B-layer [4], was added between the previous innermost pixel layer and a new, narrower beam pipe.

The high-granularity lead/liquid-argon (LAr) sampling electromagnetic calorimeter covers the region \(|\eta | < 3.2\). The regions \(|\eta | < 1.37\) and \(1.5< |\eta | < 1.8\) are instrumented with presamplers in front of the LAr calorimeter in the same cryostat. A steel/scintillator-tile calorimeter (Tile) provides hadronic coverage in the central pseudorapidity range \(|\eta | < 1.7\). LAr technology is also used for the hadronic calorimeters in the endcap region \(1.5< |\eta | < 3.2\) and for electromagnetic and hadronic energy measurements in the forward calorimeters covering \(3.2< |\eta | < 4.9\).

The MS surrounds the calorimeters. It consists of three large superconducting air-core toroidal magnets, precision tracking chambers providing precise muon tracking out to \(|\eta | = 2.7\), and fast detectors for triggering in the region \(|\eta | < 2.4\).

A two-level trigger system is used to select events [5]. A low-level hardware trigger reduces the data rate, and a high-level software trigger selects events with interesting final states. More details of the ATLAS detector can be found in Ref. [1].

3 \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction

The reconstructed \(E_{\mathrm{T}}^{\mathrm{miss}}\) in ATLAS is characterised by two contributions. The first one is from the hard-event signals comprising fully reconstructed and calibrated particles and jets (hard objects). The reconstructed particles are electrons, photons, \(\tau \text {-leptons}\), and muons. While muons are reconstructed from ID and MS tracks, electrons and \(\tau \text {-leptons}\) are identified combining calorimeter signals with tracking information. Photons and jets are principally reconstructed from calorimeter signals, with possible signal refinements from reconstructed tracks. The second contribution to \(E_{\mathrm{T}}^{\mathrm{miss}}\) is from the soft-event signals consisting of reconstructed charged-particle tracks (soft signals) associated with the hard-scatter vertex defined in Appendix A but not with the hard objects.

ATLAS carries out a dedicated reconstruction procedure for each kind of particle as well as for jets, casting a particle or jet hypothesis on the origin of (a group of) detector signals. These procedures are independent of one another. This means that e.g. the same calorimeter signal used to reconstruct an electron is likely also used to reconstruct a jet, thus potentially introducing double counting of the same signal when reconstructing \(E_{\mathrm{T}}^{\mathrm{miss}}\). This issue is addressed by the explicit signal ambiguity resolution in the object-based \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction originally introduced in Refs. [2, 3], and by its 2015 implementation described in Sects. 3.1 and 3.2.

Additional options for the set of signals used to reconstruct \(E_{\mathrm{T}}^{\mathrm{miss}}\) are available and discussed in detail in Sect. 8. One of these alternative options is the calorimeter-based \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction discussed in Sect. 8.1, which uses a soft event built from clusters of topologically connected calorimeter cells (topo-clusters) [6]. Another option is the track-based missing transverse momentum, which differs from \(E_{\mathrm{T}}^{\mathrm{miss}}\) only in the use of tracks in place of jets. It is described in more detail in Sect. 8.2.

3.1 \(E_{\mathrm{T}}^{\mathrm{miss}}\) basics

The missing transverse momentum reconstruction provides a set of observables constructed from the components \(p_{x(y)}\) of the transverse momentum vectors (\(\mathbf p _{\mathrm{T}}\)) of the various contributions. The missing transverse momentum components \(E_{x(y)}^{{\text {miss}}}\) serve as the basic input for most of these observables. They are given by

$$\begin{aligned} E_{x(y)}^{{\text {miss}}} = - \sum _{i\in \{\text {hard objects}\}} p_{x(y),i} - \sum _{j\in \{\text {soft signals}\}} p_{x(y),j} . \end{aligned}$$
(1)

The set of observables constructed from \(E_{x(y)}^{{\text {miss}}}\) is

$$\begin{aligned} \mathbf E _{\mathrm{T}}^{{\text {miss}}}= & {} (E_{x}^{{\text {miss}}},E_{y}^{{\text {miss}}}), \end{aligned}$$
(2)
$$\begin{aligned} E_{\mathrm{T}}^{\mathrm{miss}}= & {} |\mathbf E _{\mathrm{T}}^{{\text {miss}}} | = \sqrt{(E_{x}^{{\text {miss}}})^{2}+(E_{y}^{{\text {miss}}})^{2}}, \end{aligned}$$
(3)
$$\begin{aligned} \phi ^{{\text {miss}}}= & {} \tan ^{-1}(E_{y}^{{\text {miss}}}/E_{x}^{{\text {miss}}}). \end{aligned}$$
(4)

The vector \(\mathbf E _{\mathrm{T}}^{{\text {miss}}}\) provides the amount of the missing transverse momentum via its magnitude \(E_{\mathrm{T}}^{\mathrm{miss}}\), and its direction in the transverse plane in terms of the azimuthal angle \(\phi ^{{\text {miss}}}\). Consequently, \(E_{\mathrm{T}}^{\mathrm{miss}}\) is non-negative by definition. However, in an experimental environment where not all relevant \(p_{\mathrm{T}}\) from the hard-scatter interaction can be reconstructed and used in Eq. (1), and the reconstructed \(p_{\mathrm{T}}\) from each contribution is affected by the limited resolution of the detector, an observation bias towards non-vanishing values for \(E_{\mathrm{T}}^{\mathrm{miss}}\) is introduced even for final states without genuine missing transverse momentum generated by undetectable particles.

The scalar sum of all transverse momenta (\(p_{\mathrm{T}} = |\mathbf p _{\mathrm{T}} |\)) from the objects contributing to \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction is given by

$$\begin{aligned} \Sigma E_{\mathrm{T}} = \sum _{i\in \{\text {hard objects}\}} p_{\mathrm{T},i} + \sum _{j\in \{\text {soft signals}\}} p_{\mathrm{T},j}. \end{aligned}$$
(5)

In the context of \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction, \(\Sigma E_{\mathrm{T}}\) is calculated in addition to the sum given in Eq. (1), and the derived quantities defining \(E_{\mathrm{T}}^{\mathrm{miss}}\) given in Eqs. (2)–(4). It provides a useful overall scale for evaluating the hardness of the hard-scatter event in the transverse plane, and thus provides a measure for the event activity in physics analyses and \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction performance studies.

In the calculation of \(E_{x(y)}^{{\text {miss}}}\) and \(\Sigma E_{\mathrm{T}}\) the contributing objects need to be reconstructed from mutually exclusive detector signals. This rule avoids multiple inclusions of the same signal in all constructed observables. The implementation of this rule in terms of the signal ambiguity resolution requires the definition of a sequence for selected contributions, in addition to a rejection mechanism based on common signal usage between different objects. Similarly to the analysis presented in Ref. [3], the most commonly used order for the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction sequence for the hard-object contribution starts with electrons (\(e\)), followed by photons (\(\gamma \)), then hadronically decaying \(\tau \text {-leptons}\) (\(\tau _{\mathrm{had}}\)), and finally jets. Muons (\(\mu \)) are principally reconstructed from ID and MS tracks alone, with corrections based on their energy loss in the calorimeter, leading to little or no signal overlap with the other reconstructed particles in the calorimeter.

In the sequence discussed here, all electrons passing the selection enter the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction first. The lower-priority reconstructed particles (\(\gamma \), \(\tau _{\mathrm{had}}\)) are fully rejected if they share their calorimeter signal with a higher-priority object that has already entered the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction. Muons experience energy loss in the calorimeters, but only non-isolated muons overlap with other hard objects, most likely jets or \(\tau \text {-leptons}\). In this case the muon’s energy deposit in the calorimeter cannot be separated from the overlapping jet-like objects with the required precision, and the calorimeter-signal-overlap resolution based on the shared use of topo-clusters cannot be applied. A discussion of the treatment of isolated and non-isolated muons is given in Sect. 3.3.4.

Generally, jets are rejected if they overlap with accepted higher-priority particles. To avoid signal losses for \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction in the case of partial or marginal overlap, and to suppress the accidental inclusion of jets reconstructed from calorimeter signals from large muon energy losses or pile-up, the more refined overlap resolution strategies described in Sects. 3.3.5 and 3.3.6 are applied.

Excluding ID tracks associated with any of the accepted hard objects contributing to \(E_{\mathrm{T}}^{\mathrm{miss}}\), ID tracks from the hard-scatter collision vertex are used to construct the soft-event signal for the results presented in this paper.

Table 1 Overview of the contributions to \(E_{\mathrm{T}}^{\mathrm{miss}}\) and \(\Sigma E_{\mathrm{T}}\) from hard objects such as electrons (\(e\)), photons (\(\gamma \)), hadronically decaying \(\tau \text {-leptons}\) (\(\tau _{\mathrm{had}}\)), muons (\(\mu \)), and jets, together with the signals for the soft term. The configuration shown is the one used as reference for the performance evaluations presented in this paper. The table is ordered descending in priority \(\mathbb {P}\) of consideration for \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction, with (1) being the first and (5) being the last calculated hard-object contribution. The soft-event contribution is constructed at the lowest priority (6), after all hard objects are considered. The transverse (longitudinal) impact parameter \(d_{0}\) (\(z_{0} \sin (\theta )\)) used to select the ID tracks contributing to \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) and \(\Sigma E_{\mathrm{T}}^{{\text {soft}}}\) in \(\mathbb {P} = (6)\) is measured relative to the hard-scatter vertex. All variables are explained in Sect. 3.2. The angular distance \(\Delta R\) between objects is defined as \(\Delta R = \sqrt{(\Delta \eta )^{2}+(\Delta \phi )^{2}}\)

3.2 \(E_{\mathrm{T}}^{\mathrm{miss}}\) terms

Particle and jet selections in a given analysis should be reflected in \(E_{\mathrm{T}}^{\mathrm{miss}}\) and \(\Sigma E_{\mathrm{T}}\) for a consistent interpretation of a given event. Each reconstructed particle and jet has its own dedicated calibration translating the detector signals into a fully corrected four-momentum. This means that e.g. rejecting certain electrons in a given analysis can change both \(E_{\mathrm{T}}^{\mathrm{miss}}\) and \(\Sigma E_{\mathrm{T}}\), if the corresponding calorimeter signal is included and calibrated as a jet or a significant part of a jet. This also means that systematic uncertainties for the different particles can be consistently propagated to \(E_{\mathrm{T}}^{\mathrm{miss}}\). The applied selections are presented in Sect. 3.3, and summarised in Table 1.

In ATLAS the flexibility needed to recalculate \(E_{\mathrm{T}}^{\mathrm{miss}}\) and \(\Sigma E_{\mathrm{T}}\) under changing analysis requirements for the same event is implemented using dedicated variables corresponding to specific object contributions. In this approach the full \(\mathbf E _{\mathrm{T}}^{{\text {miss}}}\) is the vectorial sum of missing transverse momentum terms \(\mathbf E _{\mathrm{T}}^{{\text {miss}},p}\), with \(p \in \{e,\gamma ,\tau _{\mathrm{had}},\mu ,\text {jet}\}\) reconstructed from the \(\mathbf p _{\mathrm{T}} = (p_{x},p_{y})\) of accepted particles and jets, and the corresponding soft term \(\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {soft}}\) from the soft-event signals introduced in Sect. 3.1 and further specified in Sect. 3.4. This yieldsFootnote 2

(6)

The \(E_{\mathrm{T}}^{\mathrm{miss}}\) and \(\phi ^{{\text {miss}}}\) observables can be constructed according to Eqs. (3) and (4), respectively, for the overall missing transverse momentum (from \(\mathbf E _{\mathrm{T}}^{{\text {miss}}}\)) as well as for each individual term indicated in Eq. (6). In the priority-ordered reconstruction sequence for \(E_{\mathrm{T}}^{\mathrm{miss}}\), contributions are defined by a combination of analysis-dependent selections and a possible rejection due to the applied signal ambiguity resolution. The muon and electron contributions are typically not subjected to the signal overlap resolution and are thus exclusively defined by the selection requirements. Unused tracks in Eq. (6) refers to those tracks associated with the hard-scatter vertex but not with any hard object. Neutral particle signals from the calorimeter suffer from significant contributions from pile-up and are not included in the soft term.

Correspondingly, \(\Sigma E_{\mathrm{T}}\) is calculated from the scalar sums of the transverse momenta of hard objects entering the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction and the soft term,

(7)

The hard term in both \(E_{\mathrm{T}}^{\mathrm{miss}}\) and \(\Sigma E_{\mathrm{T}}\) is characterised by little dependence on pile-up, as it includes only fully calibrated objects, where the calibration includes a pile-up correction and objects tagged as originating from pile-up are removed. The particular choice of using only tracks from the hard-scatter vertex for the soft term strongly suppresses pile-up contributions to this term as well. The observed residual pile-up dependencies are discussed with the performance results shown in Sect. 6.

3.3 Object selection

The following selections are applied to reconstructed particles and jets used for the performance evaluations presented in Sects. 68. Generally, these selections require refinements to achieve optimal \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction performance in the context of a given physics analysis, and the selections performed in this study are an example set of criteria.

3.3.1 Electron selection

Reconstructed electrons are selected on the basis of their shower shapes in the calorimeter and how well their calorimeter cell clusters are matched to ID tracks [7]. Both are evaluated in a combined likelihood-based approach [8]. Electrons with at least medium reconstruction quality are selected. They are calibrated using the default calibration given in Ref. [7]. To be considered for \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction, electrons passing the reconstruction quality requirements are in addition required to have \(p_{\mathrm{T}} > 10\,{\text {Ge}\text {V}} \) and \(|\eta | < 1.37\) or \(1.52< |\eta | < 2.47\), to avoid the transition region between the central and endcap electromagnetic calorimeters. Any energy deposit by electrons within \(1.37< |\eta | < 1.52\) is likely reconstructed as a jet and enters the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction as such, if this jet meets the corresponding selection criteria discussed in Sect. 3.3.5.

3.3.2 Photon selection

The identification and reconstruction of photons exploits the distinctive evolution of their electromagnetic showers in the calorimeters [9]. Photons are selected and calibrated using the tight selection criteria given in Ref. [7]. In addition to the reconstruction quality requirements, photons must have \(p_{\mathrm{T}} > 25\,{\text {Ge}\text {V}} \) and \(|\eta | < 1.37\) or \(1.52< |\eta | < 2.37\) to be included in the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction. Similarly to electrons, photons emitted within \(1.37< |\eta | < 1.52\) may contribute to \(E_{\mathrm{T}}^{\mathrm{miss}}\) as a jet.

3.3.3 \(\tau \)-Lepton selection

Hadronically decaying \(\tau \text {-leptons}\) are reconstructed from narrow jets with low associated track multiplicities [10]. Candidates must pass the medium quality selection given in Ref. [11], and in addition have \(p_{\mathrm{T}} > 20\,{\text {Ge}\text {V}} \) and \(|\eta | < 1.37\) or \(1.52< |\eta | < 2.47\). Any \(\tau \text {-lepton}\) not satisfying these \(\tau \)-identification criteria may contribute to \(E_{\mathrm{T}}^{\mathrm{miss}}\) when passing the jet selection.

3.3.4 Muon selection

Muons are reconstructed within \(|\eta | < 2.5\) employing a combined MS and ID track fit. Outside of the ID coverage, muons are reconstructed within \(2.5< |\eta | < 2.7\) from a track fit to MS track segments alone. Muons are further selected for \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction by requiring the medium reconstruction quality defined in Ref. [12], \(p_{\mathrm{T}} > 10\,{\text {Ge}\text {V}} \), and an association with the hard-scatter vertex for those within \(|\eta | < 2.5\).

3.3.5 Jet selection

Jets are reconstructed from clusters of topologically connected calorimeter cells (topo-clusters), described in Ref. [6]. The topo-clusters are calibrated at the electromagnetic (EM) energy scale.Footnote 3 The anti-\(k_{t}\) algorithm [13], as provided by the FastJet toolkit [14], is employed with a radius parameter \(R = 0.4\) to form jets from these topo-clusters. The jets are fully calibrated using the EM+JES scheme [15] including a correction for pile-up [16]. They are required to have \(p_{\mathrm{T}} > 20\,{\text {Ge}\text {V}} \) after the full calibration. The jet contribution to \(E_{\mathrm{T}}^{\mathrm{miss}}\) and \(\Sigma E_{\mathrm{T}}\) is primarily defined by the signal ambiguity resolution.

Jets not rejected at that stage are further filtered using a tagging algorithm to select hard-scatter jets (“jet vertex tagging”) [16]. This algorithm provides the jet vertex tagger variable JVT, ranging from 0 (pile-up-like) to 1 (hard-scatter-like), for each jet with matched tracks.Footnote 4 The matching of tracks with jets is done by ghost association, where tracks are clustered as ghost particles into the jet, as described in Ref. [3] and based on the approach outlined in Ref. [17].

The overlap resolution can result in a partial overlap of the jet with an electron or photon, in terms of the fraction of common signals contributing to the respective reconstructed energy. This is measured by the ratio \(\kappa _{E}\) of the electron(photon) energy \(E_{e(\gamma )}^{\text {EM}}\) to the jet energy \(E_{\text {jet}}^{\text {EM}}\),

$$\begin{aligned} \kappa _{E} = \dfrac{E_{e(\gamma )}^{\text {EM}}}{E_{\text {jet}}^{\text {EM}}}, \end{aligned}$$
(8)

with both energies calibrated at the EM scale. In the case of \(\kappa _{E} \le 50\)%, the jet is included in \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction, with its \(p_{\mathrm{T}}\) scaled by \(1-\kappa _{E} \). For \(\kappa _{E} > 50\)%, only the tracks associated with the jet, excluding the track(s) associated with the overlapping particle if any, contribute to the soft term as discussed in Sect. 3.4.

Jets not rejected by the signal ambiguity resolution and with \(p_{\mathrm{T}} > 20\,{\text {Ge}\text {V}} \) and \(|\eta | > 2.4\), or with \(p_{\mathrm{T}} \ge 60\,{\text {Ge}\text {V}} \) and \(|\eta | < 4.5\), are always accepted for \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction. Jets reconstructed with \(20\,{\text {Ge}\text {V}}< p_{\mathrm{T}} < 60\,{\text {Ge}\text {V}} \) and \(|\eta | < 2.4\) are only accepted if they are tagged by \(\text {JVT} > 0.59\). In both cases, the jet \(p_{\mathrm{T}}\) thresholds are applied to the jet \(p_{\mathrm{T}}\) before applying the \(\kappa _{E}\) correction. Additional configurations for selecting jets used in \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction are discussed in Appendix A, together with the effect of the variation of these selection criteria on the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction performance.

3.3.6 Muon overlap with jets

Jets overlapping with a reconstructed muon affect the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction in a manner that depends on their origin. If these jets represent a significant (catastrophic) energy loss along the path of the muon through the calorimeter, or if they are pile-up jets tagged by JVT as originating from the hard-scatter interaction due to the muon ID track, they need to be rejected for \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction. On the other hand, jets reconstructed from final state radiation (FSR) off the muon need to be included into \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction.

In all cases, the muon–jet overlap is determined by ghost-associating the muon with the jet. For this, each muon enters the jet clustering as ghost particle with infinitesimal small momentum, together with the EM-scale topo-clusters from the calorimeter. If a given ghost particle becomes part of a jet, the corresponding muon is considered overlapping with this jet. This procedure is very similar to the track associations with jets mentioned in Sect. 3.3.5.

Tagging jets using JVT efficiently retains those from the hard-scatter vertex for \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction and rejects jets generated by pile-up. A muon overlapping with a pile-up jet can lead to a mis-tag, because the ID track from the muon represents a significant amount of \(p_{\mathrm{T}}\) from the hard-scatter vertex and thus increases JVT. As a consequence of this fake tag, the pile-up jet \(p_{\mathrm{T}}\) contributes to \(E_{\mathrm{T}}^{\mathrm{miss}}\), and thus degrades both the \(E_{\mathrm{T}}^{\mathrm{miss}}\) response and resolution due to the stochastic nature of its contribution.

A jet that is reconstructed from a catastrophic energy loss of a muon tends to be tagged as a hard-scatter jet as well. This jet is reconstructed from topo-clusters in close proximity to the extrapolated trajectory of the ID track associated with the muon bend in the axial magnetic field. Inclusion of such a jet into \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction leads to double-counting of the transverse momentum associated with the muon energy loss, as the fully reconstructed muon \(p_{\mathrm{T}}\) is already corrected for this effect.

To reject contributions from pile-up jets and jets reconstructed from muon energy loss, the following selection criteria are applied:

  • \(p_{\mathrm{T},\text {track}}^{\mu }/p_{\mathrm{T},\text {track}}^{\text {jet}} > 0.8\) – the transverse momentum of the ID track associated with the muon (\(p_{\mathrm{T},\text {track}}^{\mu }\)) represents a significant fraction of the transverse momentum \(p_{\mathrm{T},\text {track}}^{\text {jet}}\), the sum of the transverse momenta of all ID tracks associated with the jet;

  • \(p_{\mathrm{T}}^{\text {jet}}/p_{\mathrm{T},\text {track}}^{\mu } < 2\) – the overall transverse momentum \(p_{\mathrm{T}}^{\text {jet}}\) of the jet is not too large compared to \(p_{\mathrm{T},\text {track}}^{\mu }\);

  • \(N_{\text {track}}^{\mathrm{PV}} < 5\) – the total number of tracks \(N_{\text {track}}^{\mathrm{PV}}\) associated with the jet and emerging from the hard-scatter vertex is small.

All jets with overlapping muons meeting these criteria are understood to be either from pile-up or a catastrophic muon energy loss and are rejected for \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction. The muons are retained for the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction.

Another consideration for muon contributions to \(E_{\mathrm{T}}^{\mathrm{miss}}\) is FSR. Muons can radiate hard photons at small angles, which are typically not reconstructed as such because of the nearby muon ID track violating photon isolation requirements. They are also not reconstructed as electrons, due to the mismatch between the ID track momentum and the energy measured by the calorimeter. Most likely the calorimeter signal generated by the FSR photon is reconstructed as a jet, with the muon ID track associated. As the transverse momentum carried by the FSR photon is not recovered in muon reconstruction, jets representing this photon need to be included in the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction. Such jets are characterised by the following selections, which are highly indicative of a photon in the ATLAS calorimeter:

  • \(N_{\text {track}}^{\mathrm{PV}} < 3\) – the jet has low charged-particle content, indicated by a very small number of tracks from the hard-scatter vertex;

  • \(f_{\mathrm{EMC}} > 0.9\) – the jet energy \(E^{\text {jet}}\) is largely deposited in the electromagnetic calorimeter (EMC), as expected for photons and measured by the corresponding energy fraction \(f_{\mathrm{EMC}} = E^{\text {jet}}_{{\text {EMC}}}/E^{\text {jet}} \);

  • \(p_{\mathrm{T},\text {PS}}^{\text {jet}} > 2.5\,{\text {Ge}\text {V}} \) – the transverse momentum contribution \(p_{\mathrm{T},\text {PS}}^{\text {jet}}\) from presampler signals to \(p_{\mathrm{T}}^{\text {jet}}\) indicates an early starting point for the shower;

  • \(w_{\text {jet}} < 0.1\) – the jet is narrow, with a width \(w_{\text {jet}}\) comparable to a dense electromagnetic shower; \(w_{\text {jet}}\) is reconstructed according to

    $$\begin{aligned} w_{\text {jet}} = \dfrac{\sum _{i} \Delta R_{i} p_{\mathrm{T},i}}{\sum _{i}p_{\mathrm{T},i}}, \end{aligned}$$

    where \(\Delta R_{i}=\sqrt{(\Delta \eta _{i})^{2}+(\Delta \phi _{i})^{2}}\) is the angular distance of topo-cluster i from the jet axis, and \(p_{\mathrm{T},i} \) is the transverse momentum of this cluster;

  • \(p_{\mathrm{T},\text {track}}^{\text {jet}}/p_{\mathrm{T},\text {track}}^{\mu } > 0.8\) – the transverse momentum \(p_{\mathrm{T},\text {track}}^{\text {jet}}\) carried by all tracks associated with the jet is close to \(p_{\mathrm{T},\text {track}}^{\mu }\).

Jets are accepted for \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction when consistent with an FSR photon defined by the ensemble of these selection criteria, with their energy scale set to the EM scale, to improve the calibration.

3.4 \(E_{\mathrm{T}}^{\mathrm{miss}}\) soft term

The soft term introduced in Sect. 3.2 is exclusively reconstructed from ID tracks from the hard-scatter vertex, thus only using the \(p_{\mathrm{T}}\)-flow from soft charged particles. It is an important contribution to \(E_{\mathrm{T}}^{\mathrm{miss}}\) for the improvement of both the \(E_{\mathrm{T}}^{\mathrm{miss}}\) scale and resolution, in particular in final states with a low hard-object multiplicity. In this case it is indicative of (hadronic) recoil, comprising the event components not otherwise represented by reconstructed and calibrated particles or jets.

The more inclusive reconstruction of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) soft term including signals from soft neutral particles uses calorimeter topo-clusters. The reconstruction performance using the calorimeter-based \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}},\text {calo}}\) is inferior to the track-only-based \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\), mostly due to a larger residual dependence on pile-up. More details of the topo-cluster-based \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}},\text {calo}}\) reconstruction are discussed in Sect. 8.1.

3.4.1 Track and vertex selection

Hits in the ID are used to reconstruct tracks pointing to a particular collision vertex [18]. Both the tracks and vertices need to pass basic quality requirements to be accepted. Each event typically has a number \(N_{\mathrm{PV}} > 1\) of reconstructed primary vertices.

Tracks are required to have \(p_{\mathrm{T}} > 400\,{\text {Me}\text {V}} \) and \(|\eta | < 2.5\), in addition to the reconstruction quality requirements given in Ref. [19]. Vertices are constructed from at least two tracks passing selections on the transverse (longitudinal) impact parameter \(|d_{0}| < 1.5\,{\text {mm}} \) (\(|z_{0}\sin (\theta )| < 1.5\,{\text {mm}} \)) relative to the vertex candidate. These tracks must also pass requirements on the number of hits in the ID. The hard-scatter vertex is identified as described in Appendix A.

3.4.2 Track soft term

The track sample contributing to \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) for each reconstructed event is collected from high-quality tracks emerging from the hard-scatter vertex but not associated with any electron, \(\tau \text {-lepton}\), muon, or jet contributing to \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction. The applied signal-overlap resolution removes

  • ID tracks with \(\Delta R(\text {track,electron/photon cluster}) < 0.05\);

  • ID tracks with \(\Delta R(\text {track},\tau \text {-lepton}) < 0.2\);

  • ID tracks associated with muons;

  • ID tracks ghost-associated with fully or partially contributing jets.

ID tracks from the hard-scatter vertex that are associated with jets rejected by the overlap removal or are associated with jets that are likely from pile-up, as tagged by the JVT procedure discussed in Sect. 3.3.5, contribute to \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\).

Since only reconstructed tracks associated with the hard-scatter vertex are used, the track-based \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) is largely insensitive to pile-up effects. It does not include contributions from any soft neutral particles, including those produced by the hard-scatter interaction.

4 Data and simulation samples

The determination of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction performance uses selected final states without (\(E_{\mathrm{T}}^{{\text {miss}},\text {true}} = 0\)) and with genuine missing transverse momentum from neutrinos (\(E_{\mathrm{T}}^{{\text {miss}},\text {true}} = p_{\mathrm{T}}^{\nu } \)). Samples with \(E_{\mathrm{T}}^{{\text {miss}},\text {true}} = 0\) are composed of leptonic \(Z\) boson decays (\(Z \rightarrow ee\) and \(Z \rightarrow \mu \mu \)) collected by a trigger and event selection that do not depend on the particular pile-up conditions, since both the electron and muon triggers as well as the corresponding reconstructed kinematic variables are only negligibly affected by pile-up. Also using lepton triggers, samples with neutrinos were collected from \(W \rightarrow e\nu \) and \(W \rightarrow \mu \nu \) decays. In addition, samples with neutrinos and higher hard-object multiplicity were collected from top-quark pair (\(t\bar{t} \,\)) production with at least either the t or the \(\bar{t}\) decaying semi-leptonically.

4.1 Data samples

The data sample used corresponds to a total integrated luminosity of \(3.2~\hbox {fb}^{-1}\), collected with a proton bunch-crossing interval of \(25\,{\text {ns}}\). Only high-quality data with a well-functioning calorimeter, inner detector and muon spectrometer are analysed. The data-quality criteria are applied, which reduce the impact of instrumental noise and out-of-time calorimeter deposits from cosmic-ray and beam backgrounds.

4.2 Monte Carlo samples

The \(Z \rightarrow \ell \ell \) and \(W \rightarrow \ell \nu \) samples were generated using Powheg-Box  [20] (version v1r2856) employing a matrix element calculation at next-to-leading order (NLO) in perturbative QCD. To generate the particle final state, the (parton-level) matrix element output was interfaced to Pythia8  [21],Footnote 5 which generated the parton shower (PS) and the underlying event (UE) using the AZNLO tuned parameter set [22]. Parton distribution functions (PDFs) were taken from the CTEQ6L1 PDF set [23].

The \(t\bar{t}\)-production sample was generated with a Powheg NLO kernel (version v2r3026) interfaced to Pythia6  [24] (version 6.428) with the Perugia2012 set of tuned parameters [25] for the PS and UE generation. The CT10 NLO PDF set [26] was employed. The resummation of soft-gluon terms in the next-to-next-to-leading-logarithmic (NNLL) approximation with top++  2.0  [27] was included.

Additional processes contributing to the \(Z \rightarrow \ell \ell \) and \(W \rightarrow \ell \nu \) final state samples are the production of dibosons, single top quarks, and multijets. Dibosons were generated using Sherpa [28,29,30,31] version v2.1.1 employing the CT10 PDF set. Single top quarks were generated using Powheg version v1r2556 with the CT10 PDF set for the t-channel production and Powheg version v1r2819 for the s-channel and the associated top quark (\(W t\)) production, all interfaced to the PS and UE from the same Pythia6 configuration used for \(t\bar{t}\) production. Multijet events were generated using Pythia8 with the NNPDF23LO PDF set [32] and the A14 set of tuned PS and UE parameters described in Ref. [33].

Minimum bias (MB) events were generated using Pythia8 with the MSTW2008LO PDF set [34] and the A2 tuned parameter set [35] for PS and UE. These MB events were used to model pile-up, as discussed in Sect. 4.3.

For the determination of the systematic uncertainties in \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction, an alternative inclusive sample of \(Z \rightarrow \mu \mu \) events was generated using the MadGraph_aMC@NLO (version v2.2.2) matrix element generator [36] employing the CTEQ6L1 PDF set. Both PS and UE were generated using Pythia8 with the NNPDF23LO PDF set and the A14 set of tuned parameters.

The MC-generated events were processed with the Geant4 software toolkit [37], which simulates the propagation of the generated stable particlesFootnote 6 through the ATLAS detector and their interactions with the detector material [38].

4.3 Pile-up

The calorimeter signals are affected by pile-up and the short bunch-crossing period at the LHC. In 2015, an average of about 13 pile-up collisions per bunch crossing was observed. The dominant contribution of the additional \( pp \) collisions to the detector signals of the recorded event arises from a diffuse emission of soft particles superimposed to the hard-scatter interaction final state (in-time pile-up). In addition, the LAr calorimeter signals are sensitive to signal remnants from up to 24 previous bunch crossings and one following bunch crossing (out-of-time pile-up), as discussed in Refs. [6, 39]. Both types of pile-up affect signals contributing to \(E_{\mathrm{T}}^{\mathrm{miss}}\).

The in-time pile-up activity is measured by the number of reconstructed primary collision vertices \(N_{\mathrm{PV}}\). The out-of-time pile-up is proportional to the number of collisions per bunch crossing \(\mu \), measured as an average over time periods of up to two minutes by integrated signals from the luminosity detectors in ATLAS [40].

To model in-time pile-up in MC simulations, a number of generated pile-up collisions was drawn from a Poisson distribution around the value of \(\mu \) recorded in data. The collisions were randomly collected from the MB sample discussed in Sect. 4.2. The particles emerging from them were overlaid onto the particle-level final state of the generated hard-scatter interaction and converted into detector signals before event reconstruction. The event reconstruction then proceeds as for data.

Similar to the LHC proton-beam structure, events in MC simulations are organised in bunch trains, where the structure in terms of bunch-crossing interval and gaps between trains is taken into account to model the effects of out-of-time pile-up. The fully reconstructed events in MC simulation samples are finally weighted such that the distribution of the number of overlaid collisions over the whole sample corresponds to the \(\mu \) distribution observed in data.

The effect of pile-up on the signal in the Tile calorimeter is reduced due to its location behind the electromagnetic calorimeter and its fast time response [41]. Reconstructed ID and MS tracks are largely unaffected by pile-up.

5 Event selection

5.1 \(Z \rightarrow \mu \mu \) event selection

The \(Z \rightarrow \mu \mu \) final state is ideal for the evaluation of \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction performance, since it can be selected with a high signal-to-background ratio and the \(Z\) kinematics can be measured with high precision, even in the presence of pile-up. Neutrinos are produced only through very rare heavy-flavour decays in the hadronic recoil. This channel can therefore be considered to have no genuine missing transverse momentum. Thus, the scale and resolution for the reconstructed \(E_{\mathrm{T}}^{\mathrm{miss}}\) are indicative of the reconstruction quality and reflect limitations introduced by both the detector and the ambiguity resolution procedure. The well-defined expectation value \(E_{\mathrm{T}}^{{\text {miss}},\text {true}} = 0\) allows the reconstruction quality to be determined in both data and MC simulations. The reconstructed \(E_{\mathrm{T}}^{\mathrm{miss}}\) in this final state is also sensitive to the effectiveness of the muon–jet overlap resolution, which can be explored in this low-multiplicity environment in both data and MC simulations, with a well-defined \(E_{\mathrm{T}}^{\mathrm{miss}}\).

Events must pass one of three high-level muon triggers with different \(p_{\mathrm{T}}^{\mu }\) thresholds and isolation requirements. The isolation is determined by the ratio of the scalar sum of \(p_{\mathrm{T}}\) of reconstructed tracks other than the muon track itself, in a cone of size \(\Delta R = 0.2\) around the muon track (\(p_{\mathrm{T}}^{\text {cone}}\)), to \(p_{\mathrm{T}}^{\mu }\). The individual triggers require (1) \(p_{\mathrm{T}}^{\mu } > 20\,{\text {Ge}\text {V}} \) and \(p_{\mathrm{T}}^{\text {cone}}/p_{\mathrm{T}}^{\mu } < 0.12\), or (2) \(p_{\mathrm{T}}^{\mu } > 24\,{\text {Ge}\text {V}} \) and \(p_{\mathrm{T}}^{\text {cone}}/p_{\mathrm{T}}^{\mu } < 0.06\), or (3) \(p_{\mathrm{T}}^{\mu } > 50\,{\text {Ge}\text {V}} \) without isolation requirement.

The offline selection of \(Z \rightarrow \mu \mu \) events requires exactly two muons, each selected as defined in Sect. 3.3.4, with the additional criteria that (1) the muons must have opposite charge, (2) \(p_{\mathrm{T}}^{\mu } > 25\,{\text {Ge}\text {V}} \), and (3) the reconstructed invariant mass \(m_{\mu \mu }\) of the dimuon system is consistent with the mass \(m_{Z}\) of the \(Z\) boson, \(|m_{\mu \mu }- m_{Z} | < 25\,{\text {Ge}\text {V}} \).

5.2 \(W \rightarrow e\nu \) event selection

Events with \(W \rightarrow e\nu \) or \(W \rightarrow \mu \nu \) in the final state provide a well-defined topology with neutrinos produced in the hard-scatter interaction. In combination with \(Z \rightarrow \mu \mu \), the effectiveness of signal ambiguity resolution and lepton energy reconstruction for both the electrons and muons can be observed. The \(W \rightarrow e\nu \) events in particular provide a good metric with \(E_{\mathrm{T}}^{{\text {miss}},\text {true}} = p_{\mathrm{T}}^{\nu } > 0\) to evaluate and validate the scale, resolution and direction (azimuth) of the reconstructed \(E_{\mathrm{T}}^{\mathrm{miss}}\), as the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction is sensitive to the electron–jet overlap resolution performance. This metric is only available in MC simulations where \(p_{\mathrm{T}}^{\nu }\) is known. Candidate \(W \rightarrow e\nu \) events are required to pass the high-level electron trigger with \(p_{\mathrm{T}} > 17\,{\text {Ge}\text {V}} \). Electron candidates are selected according to criteria described in Sect. 3.3.1. Only events containing exactly one electron are considered.

Further selections using \(E_{\mathrm{T}}^{\mathrm{miss}}\) and the reconstructed transverse mass \(m_{\mathrm{T}}\), given by

$$\begin{aligned} m_{\mathrm{T}} = \sqrt{ 2p_{\mathrm{T}}^{e} E_{\mathrm{T}}^{\mathrm{miss}} (1-\cos {\Delta \phi })}, \end{aligned}$$

are applied to reduce the multijet background with one jet emulating an isolated electron from the \(W\) boson. Here \(E_{\mathrm{T}}^{\mathrm{miss}}\) is calculated as presented in Sect. 3. The transverse momentum of the electron is denoted by \(p_{\mathrm{T}}^{e}\), and \(\Delta \phi \) is the distance between \(\phi ^{{\text {miss}}}\) and the azimuth of the electron. Selected events are required to have \(E_{\mathrm{T}}^{\mathrm{miss}} > 25\,{\text {Ge}\text {V}} \) and \(m_{\mathrm{T}} > 50\,{\text {Ge}\text {V}} \).

5.3 \(t\bar{t}\) event selection

Events with \(t\bar{t}\) in the final state allow the evaluation of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) performance in interactions with a large jet multiplicity. Electrons and muons used to define these samples are reconstructed as discussed in Sects. 3.3.1 and 3.3.4, respectively, and are required to have \(p_{\mathrm{T}} > 25\,{\text {Ge}\text {V}} \).

The final \(t\bar{t}\) sample is selected by imposing additional requirements. Each event must have exactly one electron and no muons passing the selections described above. In addition, at least four jets reconstructed by the anti-\(k_{t}\) algorithm with \(R = 0.4\) and selected following the description in Sect. 3.3.5 are required. At least one of the jets needs to be b-tagged using the tagger configuration for a 77% efficiency working point described in Ref. [42]. All jets are required to be at an angular distance of \(\Delta R > 0.4\) from the electron.

6 Performance of \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction in data and Monte Carlo simulation

Unlike for fully reconstructed and calibrated particles and jets, and in the case of the precise reconstruction of charged particle kinematics provided by ID tracks, \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction yields a non-linear response, especially in regions of phase space where the observation bias discussed in Sect. 3.1 dominates the reconstructed \(E_{\mathrm{T}}^{\mathrm{miss}}\). In addition, the \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution functions are characterised by a high level of complexity, due to the composite character of the observable. Objects with different \(p_{\mathrm{T}}\)-resolutions contribute, and the \(E_{\mathrm{T}}^{\mathrm{miss}}\) composition can fluctuate significantly for events from the same final state. Due to the dependence of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) response on the resolution, both performance characteristics change as a function of the total event activity and are affected by pile-up. There is no universal way of mitigating these effects, due to the inability to validate in data a stable and universal calibration reference for \(E_{\mathrm{T}}^{\mathrm{miss}}\).

The \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction performance is therefore assessed by comparing a set of reconstructed \(E_{\mathrm{T}}^{\mathrm{miss}}\)-related observables in data and MC simulations for the same final-state selection, with the same object and event selections applied. Systematic uncertainties in the \(E_{\mathrm{T}}^{\mathrm{miss}}\) response and resolution are derived from these comparisons and are used to quantify the level of understanding of the data from the physics models. The quality of the detector simulation is independently determined for all reconstructed jets, particles and ID tracks, and can thus be propagated to the overall \(E_{\mathrm{T}}^{\mathrm{miss}}\) uncertainty for any given event. Both the distributions of observables as well as their average behaviour with respect to relevant scales measuring the overall kinematic activity of the hard-scatter event or the pile-up activity are compared. To focus on distribution shapes rather than statistical differences in these comparisons, the overall distribution of a given observable obtained from MC simulations is normalised to the integral of the corresponding distribution in data.

As the reconstructed final state can be produced by different physics processes, the individual process contributions in MC simulations are scaled according to the cross section of the process. This approach is taken to both show the contribution of a given process to the overall distribution, and to identify possible inadequate modelling arising from any individual process, or a subset of processes, by its effect on the overall shape of the MC distribution.

Inclusive event samples considered for the \(E_{\mathrm{T}}^{\mathrm{miss}}\) performance evaluation are obtained by applying selections according to Sect. 5.1 for a final state without genuine \(E_{\mathrm{T}}^{\mathrm{miss}}\) (\(Z \rightarrow \mu \mu \)), and according to Sect. 5.2 for a final state with genuine \(E_{\mathrm{T}}^{\mathrm{miss}}\) (\(W \rightarrow e\nu \)). From these, specific exclusive samples are extracted by applying conditions on the number of jets reconstructed. In particular, zero jet (\(N_{\mathrm{jet}} = 0\)) samples without any jet with \(p_{\mathrm{T}} > 20\,{\text {Ge}\text {V}} \) (fully calibrated) and \(|\eta |<4.9\) are useful for exclusively studying the performance of the soft term. Samples with events selected on the basis of a non-zero number of reconstructed jets with \(p_{\mathrm{T}} > 20\,{\text {Ge}\text {V}} \) are useful for evaluating the contribution of jets to \(E_{\mathrm{T}}^{\mathrm{miss}}\). While the \(p_{\mathrm{T}}\) response of jets is fully calibrated and provides a better measurement of the overall event \(p_{\mathrm{T}}\)-flow, the \(p_{\mathrm{T}}\) resolution for jets is affected by pile-up and can introduce a detrimental effect on \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction performance.

Missing transverse momentum and its related observables presented in Sect. 3.1 are reconstructed for the performance evaluations shown in the following sections using a standard reconstruction configuration. This configuration implements the signal ambiguity resolution in the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction sequence discussed in Sect. 3.1. It employs the hard-object selections defined in Sects. 3.3.13.3.4, with jets selected according to the prescriptions given in Sect. 3.3.5. The overlap resolution strategy for jets and muons described in Sect. 3.3.6 is applied. The soft term is formed from ID tracks according to Sect. 3.4.

6.1 \(E_{\mathrm{T}}^{\mathrm{miss}}\) modelling in Monte Carlo simulations

Fig. 1
figure 1

Distributions of a \(E_{\mathrm{T}}^{\mathrm{miss}}\), b \(\Sigma E_{\mathrm{T}}\), c \(E_{x}^{{\text {miss}}}\) and d \(E_{y}^{{\text {miss}}}\) for an inclusive sample of \(Z \rightarrow \mu \mu \) events extracted from data and compared to MC simulations including all relevant backgrounds. The shaded areas indicate the total uncertainty for MC simulations, including the overall statistical uncertainty combined with systematic uncertainties from the \(p_{\mathrm{T}}\) scale and resolution which are contributed by muons, jets, and the soft term. The last bin of each distribution includes the overflow, and the first bin contains the underflow in c and d. The respective ratios between data and MC simulations are shown below the distributions, with the shaded areas showing the total uncertainties for MC simulations

Fig. 2
figure 2

Distributions of a the jet term \(E_{\mathrm{T}}^{{\text {miss}},\text {jet}}\), b the muon term \(E_{\mathrm{T}}^{{\text {miss}},\mu }\), and c the soft term \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) for the inclusive samples of \(Z \rightarrow \mu \mu \) events in data, compared to MC simulations including all relevant backgrounds. The shaded areas indicate the total uncertainty from MC simulations, including the overall statistical uncertainty combined with the respective systematic uncertainties from a the jet, b the muon, and c the soft term. The last bin of each distribution includes the overflow entries. The respective ratios between data and MC simulations are shown below the distributions, with the shaded areas showing the corresponding total uncertainties from MC simulations

The quality of the MC modelling of \(E_{x}^{{\text {miss}}}\), \(E_{y}^{{\text {miss}}}\), \(E_{\mathrm{T}}^{\mathrm{miss}}\) and \(\Sigma E_{\mathrm{T}}\), reconstructed as given in Eqs. (1), (3) and (5), is evaluated for an inclusive sample of \(Z \rightarrow \mu \mu \) events by comparing the distributions of these observables to data. The results are presented in Fig. 1. The data and MC simulations agree within 20% for the bulk of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) distribution shown in Fig. 1a, with larger differences not accommodated by the total (systematic and statistical) uncertainties of the distributions for high \(E_{\mathrm{T}}^{\mathrm{miss}}\). These differences suggest a mismodelling in \(t\bar{t}\) events, the dominant background in the tail regime [43]. The \(\Sigma E_{\mathrm{T}}\) distributions compared between data and MC simulations in Fig. 1b show discrepancies significantly larger than the overall uncertainties for \(200\,{\text {Ge}\text {V}}< \Sigma E_{\mathrm{T}} < 1.2\,{\text {Te}\text {V}} \). These reflect the level of mismodelling of the final state mostly in terms of hard-object composition in MC simulations. The \(E_{x}^{{\text {miss}}}\) and \(E_{y}^{{\text {miss}}}\) spectra shown in Fig. 1c, d, respectively, show good agreement between data and MC simulations for the bulk of the distributions within \(|E_{x(y)}^{{\text {miss}}} |<100\,{\text {Ge}\text {V}} \), with larger differences observed outside of this range still mostly within the uncertainties.

Fig. 3
figure 3

Distributions of the total \(E_{\mathrm{T}}^{\mathrm{miss}}\) in a the inclusive case and b the \(N_{\mathrm{jet}} = 0\) case, as well as c the soft term \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) reconstructed in \(N_{\mathrm{jet}} = 0\) events with \(W \rightarrow e\nu \) in data. The expectation from MC simulation is superimposed and includes all relevant background final states passing the event selection. The inclusive \(E_{\mathrm{T}}^{\mathrm{miss}}\) distribution from MC simulations contains a small contribution from multijet final states at low \(E_{\mathrm{T}}^{\mathrm{miss}}\), which is absent for the \(N_{\mathrm{jet}} = 0\) selection. The shaded areas indicate the total uncertainty for MC simulations, including the overall statistical uncertainty combined with systematic uncertainties comprising contributions from the electron, jet, and the soft term. The last bins contain the respective overflows. The respective ratios between data and MC simulations are shown below the distributions, with the shaded areas indicating the total uncertainties for MC simulations

The distributions of individual contributions to \(E_{\mathrm{T}}^{\mathrm{miss}}\) from jets (\(E_{\mathrm{T}}^{{\text {miss}},\text {jet}}\)), muons (\(E_{\mathrm{T}}^{{\text {miss}},\mu }\)), and the soft term (\(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\)), as defined in Eq. (6), are compared between data and MC simulations for the same inclusive \(Z \rightarrow \mu \mu \) sample in Fig. 2. Agreement between data and MC simulations for \(E_{\mathrm{T}}^{{\text {miss}},\text {jet}}\) in Fig. 2a is of the order of \(\pm 20\)% and within the total uncertainties for \(E_{\mathrm{T}}^{{\text {miss}},\text {jet}} \lesssim 120\,{\text {Ge}\text {V}} \), but beyond those for higher \(E_{\mathrm{T}}^{{\text {miss}},\text {jet}}\). A similar observation holds for \(E_{\mathrm{T}}^{{\text {miss}},\mu }\) in Fig. 2b, where data and MC simulations agree within the uncertainties for low \(E_{\mathrm{T}}^{{\text {miss}},\mu }\) but significantly beyond them for larger \(E_{\mathrm{T}}^{{\text {miss}},\mu }\). Agreement between data and MC simulations is better for the soft term \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\), with differences up to 10% for \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}} \lesssim 30\,{\text {Ge}\text {V}} \), as seen in Fig. 2c. Larger differences for larger \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) are still found to be within the uncertainties.

The peak around \(E_{\mathrm{T}}^{{\text {miss}},\text {jet}} = 20\,{\text {Ge}\text {V}} \) indicates the onset of single-jet events at the threshold \(p_{\mathrm{T}} = 20\,{\text {Ge}\text {V}} \) for jets contributing to \(E_{\mathrm{T}}^{{\text {miss}},\text {jet}}\). Larger values of \(E_{\mathrm{T}}^{{\text {miss}},\text {jet}}\) arise from events with one or more high-\(p_{\mathrm{T}}\) jets balancing the \(p_{\mathrm{T}}\) of the \(Z\) boson.

For the \(W \rightarrow e\nu \) sample with genuine missing transverse momentum given by \(p_{\mathrm{T}}^{\nu }\), both the total reconstructed \(E_{\mathrm{T}}^{\mathrm{miss}}\) and the soft term are compared between data and MC simulations in Fig. 3. The level of agreement between the \(E_{\mathrm{T}}^{\mathrm{miss}}\) distributions for data and MC simulations shown in Fig. 3a for the inclusive event sample is at \(\pm 20\%\), similar to that observed for the \(Z \rightarrow \mu \mu \) sample in Fig. 1a, except that for this final state it is found to be within the total uncertainties of the measurement. The differences between the \(E_{\mathrm{T}}^{\mathrm{miss}}\) distributions observed with the exclusive \(N_{\mathrm{jet}} = 0\) sample shown in Fig. 3b are well below 20%, but show a trend to larger discrepancies for decreasing \(E_{\mathrm{T}}^{\mathrm{miss}} \lesssim 40\,{\text {Ge}\text {V}} \). This trend is due to the missing background contribution in MC simulations from multijet final states. The extraction of this contribution is very inefficient and only possible with large statistical uncertainties. Even very large MC samples of multijet final states provide very few events with only one jet that is accidentally reconstructed as an electron, and with the amount of \(E_{\mathrm{T}}^{\mathrm{miss}}\) required in the \(W \rightarrow e\nu \) selection described in Sect. 5.2. The comparison of the \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) distributions from data and MC simulations shown in Fig. 3c yields agreement well within the uncertainties, for \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}} \gtrsim 10\,{\text {Ge}\text {V}} \). The rising deficiencies observed in the MC distribution for decreasing \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}} \lesssim 10\,{\text {Ge}\text {V}} \) are expected to be related to the missing multijet contribution.

6.2 \(E_{\mathrm{T}}^{\mathrm{miss}}\) response and resolution

The response in the context of \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction is determined by the deviation of the observed \(E_{\mathrm{T}}^{\mathrm{miss}}\) from the expectation value for a given final state. This deviation sets the scale for the observed \(E_{\mathrm{T}}^{\mathrm{miss}}\). If this deviation is independent of the genuine missing transverse momentum, or any other hard \(p_{\mathrm{T}}\) indicative of the overall hard-scatter activity, the \(E_{\mathrm{T}}^{\mathrm{miss}}\) response is linear. In this case, a constant bias in the reconstructed \(E_{\mathrm{T}}^{\mathrm{miss}}\) is still possible due to detector inefficiencies and coverage (acceptance) limitations.

Final states balanced in transverse momentum are expected to show a non-linear \(E_{\mathrm{T}}^{\mathrm{miss}}\) response at low event activity, as the response in this case suffers from the observation bias in \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction discussed in Sect. 3.1. With increasing momentum transfers in the hard-scatter interaction, the \(E_{\mathrm{T}}^{\mathrm{miss}}\) response becomes increasingly dominated by a well-measured hadronic recoil and thus more linear. In the case of final states with genuine missing transverse momentum, the \(E_{\mathrm{T}}^{\mathrm{miss}}\) response is only linear once \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\) exceeds the observation bias. These features are discussed in Sect. 6.2.1 and explored in Sect. 6.2.2.

Contributions to the fluctuations in the \(E_{\mathrm{T}}^{\mathrm{miss}}\) measurement arise from (1) the limitations in the detector acceptance not allowing the reconstruction of the complete transverse momentum flow from the hard interaction, (2) the irreducible intrinsic signal fluctuations in the detector response, and from (3) the additional response fluctuations due to pile-up. In particular (1) introduces fluctuations driven by the large variations of the particle composition of the final state with respect to their types, momenta and directions. The limited detector coverage of \(|\eta | < 4.9\) for all particles, together with the need to suppress the pile-up-induced signal fluctuations as much as possible, restricts the contribution of particles to \(E_{\mathrm{T}}^{\mathrm{miss}}\) to the reconstructed and accepted \(e\), \(\gamma \), \(\tau _{\mathrm{had}}\) and \(\mu \), and those being part of a reconstructed and accepted jet. In addition, the \(p_{\mathrm{T}}\)-flow of not explicitly reconstructed charged particles emerging from the hard-scatter vertex is represented by ID tracks contributing to \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) given in Eqs. (6) and (7), but only in the phase space defined by the selections given in Sect. 3.4.1. All other charged and neutral particles do not contribute to \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction.

Like for the \(E_{\mathrm{T}}^{\mathrm{miss}}\) response, resolution-related aspects of \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction are understood from data-to-MC-simulations comparisons. The scales used for the corresponding evaluations are the overall event activity represented by \(\Sigma E_{\mathrm{T}}\), and the pile-up activity measured by \(N_{\mathrm{PV}}\). The measurement of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution is discussed in Sect. 6.2.3 and results are presented in Sect. 6.2.4.

6.2.1 \(E_{\mathrm{T}}^{\mathrm{miss}}\) scale determination

In events with \(Z \rightarrow \mu \mu \) decays, the transverse momentum of the \(Z\) boson (\(p_{\mathrm{T}}^{Z}\)) is an indicator of the hardness of the interaction. It provides a useful scale for the evaluation of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) response for this final state without genuine missing transverse momentum. The direction of the corresponding \(Z\) boson transverse momentum vector \(\mathbf p _{\mathrm{T}}^{Z}\) defines an axis \(\mathbf A _{Z}\) in the transverse plane of the collision, which is reconstructed from the \(\mathbf p _{\mathrm{T}}\) of the decay products by

(9)

The magnitude of the component of \(\mathbf E _{\mathrm{T}}^{{\text {miss}}}\) parallel to \(\mathbf A _{Z}\) is

$$\begin{aligned} \mathcal {P} _{\parallel }^{\,\!Z} = \mathbf E _{\mathrm{T}}^{{\text {miss}}} \cdot \mathbf A _{Z} . \end{aligned}$$
(10)

This projection is sensitive to any limitation in \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction, in particular with respect to the contribution from the hadronic recoil against \(\mathbf p _{\mathrm{T}}^{Z}\), both in terms of response and resolution. Because it can be determined both for data and MC simulations, it provides an important tool for the validation of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) response and the associated systematic uncertainties.

The expectation value for a balanced interaction producing a \(Z\) boson against a hadronic recoil is \(\text {E}[\mathcal {P} _{\parallel }^{\,\!Z} ] = 0\). Any observed deviation from this value represents a bias in the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction. For \(\mathcal {P} _{\parallel }^{\,\!Z} < 0\), the reconstructed hadronic activity recoiling against \(\mathbf p _{\mathrm{T}}^{Z}\) is too small, while for \(\mathcal {P} _{\parallel }^{\,\!Z} > 0\) too much hadronic recoil is reconstructed. The evolution of \(\mathcal {P} _{\parallel }^{\,\!Z}\) as a function of the hardness of the \(Z\) boson production can be measured by evaluating the mean \(\langle \mathcal {P} _{\parallel }^{\,\!Z} \rangle \) in bins of the hard-scatter scale \(p_{\mathrm{T}}^{{\text {hard}}} = p_{\mathrm{T}}^{Z} \).

In addition to measuring the \(E_{\mathrm{T}}^{\mathrm{miss}}\) response in data and MC simulation without genuine \(E_{\mathrm{T}}^{\mathrm{miss}}\), its linearity can be determined using samples of final states with genuine \(E_{\mathrm{T}}^{\mathrm{miss}}\) in MC simulations. This is done by evaluating the relative deviation \(\Delta _{\mathrm{T}}^{\mathrm{lin}}\) of the reconstructed \(E_{\mathrm{T}}^{\mathrm{miss}}\) from the expected \(E_{\mathrm{T}}^{{\text {miss}},\text {true}} > 0\) as a function of \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\),

$$\begin{aligned} \Delta _{\mathrm{T}}^{\mathrm{lin}} (E_{\mathrm{T}}^{{\text {miss}},\text {true}}) = \dfrac{E_{\mathrm{T}}^{\mathrm{miss}}- E_{\mathrm{T}}^{{\text {miss}},\text {true}}}{E_{\mathrm{T}}^{{\text {miss}},\text {true}}}. \end{aligned}$$
(11)

6.2.2 Measuring the \(E_{\mathrm{T}}^{\mathrm{miss}}\) response

Fig. 4
figure 4

The average projection of \(\mathbf E _{\mathrm{T}}^{{\text {miss}}}\) onto the direction \(\mathbf A _{Z}\) of the \(Z\) boson’s transverse momentum vector \(\mathbf p _{\mathrm{T}}^{Z}\), as given in Eq. (10), is shown as a function of \(p_{\mathrm{T}}^{Z} = |\mathbf p _{\mathrm{T}}^{Z} |\) in \(Z \rightarrow \mu \mu \) events from a the \(N_{\mathrm{jet}} = 0\) sample and from b the inclusive sample. In both cases data are compared to MC simulations. The ratio of the averages from data and MC simulations are shown below the plots. The shaded areas indicate the overall statistical uncertainty combined with systematic uncertainties comprising contributions from the muon and soft-term systematic uncertainties in a, and including the additional jet systematic uncertainties in b, for MC simulations

Figure 4 shows \(\langle \mathcal {P} _{\parallel }^{\,\!Z} \rangle \) as a function of \(p_{\mathrm{T}}^{Z}\) for the \(N_{\mathrm{jet}} = 0\) and the inclusive \(Z \rightarrow \mu \mu \) sample, respectively. MC simulations compare well with the data for \(N_{\mathrm{jet}} = 0\), but show larger deviations up to 30% for the inclusive selection. Nevertheless, these differences are still found to be within the total uncertainty of the measurement.

The steep decrease of \(\langle \mathcal {P} _{\parallel }^{\,\!Z} \rangle \) with increasing \(p_{\mathrm{T}}^{Z}\) in the \(N_{\mathrm{jet}} = 0\) sample seen in Fig. 4a reflects the inherent underestimation of the soft term, as in this case the hadronic recoil is exclusively represented by ID tracks with \(p_{\mathrm{T}} > 400\,{\text {Me}\text {V}} \) within \(|\eta | < 2.5\). It thus does not contain any signal from (1) neutral particles, (2) charged particles produced with \(|\eta | > 2.5\), and (3) charged particles produced within \(|\eta | < 2.5\) but with \(p_{\mathrm{T}}\) below threshold, rejected by the track quality requirements, or not represented by a track at all due to insufficient signals in the ID (e.g., lack of hits for track fitting).

In the case of the inclusive sample shown in Fig. 4b, the \(E_{\mathrm{T}}^{\mathrm{miss}}\) response is recovered better as \(p_{\mathrm{T}}^{Z}\) increases, since an increasing number of events enter the sample with a reconstructed recoil containing fully calibrated jets. These provide a more complete representation of the hadronic transverse momentum flow. The residual offsets in \(\langle \mathcal {P} _{\parallel }^{\,\!Z} \rangle \) of about \(8\,{\text {Ge}\text {V}}\) in data and \(6\,{\text {Ge}\text {V}}\) in MC simulations observed for \(p_{\mathrm{T}}^{Z} \gtrsim 40\,{\text {Ge}\text {V}} \) in Fig. 4b agree within the uncertainties of this measurement.

The persistent bias in \(\langle \mathcal {P} _{\parallel }^{\,\!Z} \rangle \) is further explored in Fig. 5, which compares variations of \(\langle \mathcal {P} _{\parallel }^{\,\!Z} \rangle \) respectively using the full \(\mathbf E _{\mathrm{T}}^{{\text {miss}}}\), the soft-term contribution \(\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {soft}}\) only, the hard-term contribution \(\mathbf E _{\mathrm{T}}^{{\text {miss}}}-\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {soft}} \), and the true soft term \(\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {true soft}}\) only, as a function of \(p_{\mathrm{T}}^{Z}\), for the \(Z \rightarrow \mu \mu \) sample from MC simulations. In particular the difference between the projections using \(\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {true soft}}\) and \(\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {soft}}\) indicates the lack of reconstructed hadronic response, when \(\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {soft}} = \mathbf E _{\mathrm{T}}^{{\text {miss}},\text {true soft}} \) is expected for a fully measured recoil. The parallel projection using only the soft terms is larger than zero for all \(p_{\mathrm{T}}^{Z}\) due to the missing \(Z\)-boson contribution to \(\mathbf E _{\mathrm{T}}^{{\text {miss}}}\) given by \(-\mathbf p _{\mathrm{T}}^{Z} \).

Fig. 5
figure 5

The average projection of \(\mathbf E _{\mathrm{T}}^{{\text {miss}}}\) onto the direction \(\mathbf A _{Z}\) of the \(Z\) boson’s transverse momentum vector \(\mathbf p _{\mathrm{T}}^{Z}\), as given in Eq. (10), is shown as a function of \(p_{\mathrm{T}}^{Z} = |\mathbf p _{\mathrm{T}}^{Z} |\) in \(Z \rightarrow \mu \mu \) events from the inclusive MC sample. The average projection of the soft term and the true soft term are also shown, to demonstrate the source of the deviation from zero

Fig. 6
figure 6

The deviation of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) response from linearity, measured as a function of the expected \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\) by \(\Delta _{\mathrm{T}}^{\mathrm{lin}}\) in Eq. (11), in \(W \rightarrow e\nu \), \(W \rightarrow \mu \nu \), and \(t\bar{t}\) final states in MC simulations. The lower plot shows a zoomed-in view on the \(\Delta _{\mathrm{T}}^{\mathrm{lin}}\) dependence on \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\) with a highly suppressed ordinate

The deviation from linearity in \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction, measured by \(\Delta _{\mathrm{T}}^{\mathrm{lin}}\) given in Eq. (11), is shown as a function of \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\) for MC simulations of \(W \rightarrow e\nu \), \(W \rightarrow \mu \nu \) and \(t\bar{t}\) production in Fig. 6. The observed \(\Delta _{\mathrm{T}}^{\mathrm{lin}} > 0\) at low \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\) indicates an overestimation of \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\) by the reconstructed \(E_{\mathrm{T}}^{\mathrm{miss}}\) due to the observation biases arising from the finite \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution, as discussed in Sect. 3.1. This bias overcompensates the lack of reconstructed \(p_{\mathrm{T}}\)-flow from the incompletely measured hadronic recoil in \(W \rightarrow e\nu \) and \(W \rightarrow \mu \nu \) events for \(E_{\mathrm{T}}^{{\text {miss}},\text {true}} \lesssim 40\,{\text {Ge}\text {V}} \) with an increasing non-linearity observed with decreasing \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\). For \(E_{\mathrm{T}}^{{\text {miss}},\text {true}} \gtrsim 70\,{\text {Ge}\text {V}} \) the \(E_{\mathrm{T}}^{\mathrm{miss}}\) response is directly proportional to \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\), with the reconstructed recoil being approximately 2% too small. The \(W \rightarrow e\nu \) and \(W \rightarrow \mu \nu \) final states show very similar \(\Delta _{\mathrm{T}}^{\mathrm{lin}} (E_{\mathrm{T}}^{{\text {miss}},\text {true}})\), thus indicating the universality of the recoil reconstruction and the independence on the lepton flavour of the reconstructed \(E_{\mathrm{T}}^{\mathrm{miss}}\) in a low-multiplicity final state with \(E_{\mathrm{T}}^{{\text {miss}},\text {true}} > 0\).

In \(t\bar{t}\) final-state reconstruction, resolution effects tend to dominate \(\Delta _{\mathrm{T}}^{\mathrm{lin}}\) at \(E_{\mathrm{T}}^{{\text {miss}},\text {true}} \lesssim 120\,{\text {Ge}\text {V}} \). Compared to the \(W \rightarrow e\nu \) and \(W \rightarrow \mu \nu \) final states, a significantly poorer \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution is observed in this kinematic region, due to the presence of at least four jets with relatively low \(p_{\mathrm{T}}\) and high sensitivity to pile-up-induced fluctuations in each event of the \(t\bar{t}\) sample. For \(E_{\mathrm{T}}^{{\text {miss}},\text {true}} > 120\,{\text {Ge}\text {V}} \), \(\Delta _{\mathrm{T}}^{\mathrm{lin}} (E_{\mathrm{T}}^{{\text {miss}},\text {true}}) \approx 2\)% indicates a proportional \(E_{\mathrm{T}}^{\mathrm{miss}}\) response with a systematic shift similar to the one observed in inclusive \(W\)-boson production.

6.2.3 Determination of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution

The \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution is determined by the width of the combined distribution of the differences between the measured \(E_{x(y)}^{{\text {miss}}}\) and the components of the true missing transverse momentum vector \(\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {true}} = (E_{x}^{{\text {miss}},{\text {true}}},E_{y}^{{\text {miss}},{\text {true}}})\). The width is measured in terms of the \({\text {RMS}}\), with

$$\begin{aligned} \mathrm {RMS}^{\mathrm{miss}}_{x(y)} = \left\{ \begin{array}{ll} {\text {RMS}}(E_{x(y)}^{{\text {miss}}} - E_{x(y)}^{{\text {miss}},\text {true}}) &{} W \rightarrow e\nu \text { or } t\bar{t} \text { sample } (E_{\mathrm{T}}^{{\text {miss}},\text {true}} > 0) \\ {\text {RMS}}(E_{x(y)}^{{\text {miss}}}) &{} Z \rightarrow \mu \mu \text { sample } (E_{\mathrm{T}}^{{\text {miss}},\text {true}} = 0) \end{array}\right. . \end{aligned}$$
(12)

This metric does not capture all of the effects driving the fluctuations in \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction, such as biases between individual \(E_{\mathrm{T}}^{\mathrm{miss}}\) terms or the behaviour of outliers, but it is an appropriate general measure of how well \(E_{\mathrm{T}}^{\mathrm{miss}}\) represents \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\).

Using the \(Z \rightarrow \mu \mu \) sample allows direct comparisons of \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) between data and MC simulations, as \(E_{\mathrm{T}}^{{\text {miss}},\text {true}} = 0\) in this case. The resolution in final states with genuine \(E_{\mathrm{T}}^{\mathrm{miss}}\) is determined with MC simulations alone. For \(W \rightarrow e\nu \) and \(t\bar{t}\) final states, \(E_{x(y)}^{{\text {miss}},\text {true}} = p_{x(y)}^{\nu } \) is used.

6.2.4 \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution measurements

Fig. 7
figure 7

The RMS width of the \(E_{x(y)}^{{\text {miss}}}\) distributions a in bins of \(\Sigma E_{\mathrm{T}}\) and b in bins of the number of primary vertices in an inclusive sample of \(Z \rightarrow \mu \mu \) events. Predictions from MC simulations are overlaid on the data points, and the ratios are shown below the respective plot. The shaded bands indicate the combined statistical and systematic uncertainties of the resolution measurements

Fig. 8
figure 8

The \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) determined for a an exclusive \(Z \rightarrow \mu \mu \) sample without jets with \(p_{\mathrm{T}} > 20\,{\text {Ge}\text {V}} \) (\(N_{\mathrm{jet}} = 0\)) and for b an exclusive sample with at least one jet above this threshold (\(N_{\mathrm{jet}} \ge 1\)), as a function of \(\Sigma E_{\mathrm{T}}\) in data and MC simulations. The dependence of \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) on the pile-up activity, as measured by \(N_{\mathrm{PV}}\), for these two samples is shown in c and d, respectively. The shaded bands indicate the combined statistical and systematic uncertainties associated with the measurement

The \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution measured by \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) is evaluated as a function of the event activity measured by \(\Sigma E_{\mathrm{T}}\) given in Eq. (7). For the inclusive \(Z \rightarrow \mu \mu \) sample, Fig. 7a shows \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) quickly rising from less than \(5\,{\text {Ge}\text {V}}\) to about \(10\,{\text {Ge}\text {V}}\) with increasing \(\Sigma E_{\mathrm{T}}\) within \(50\,{\text {Ge}\text {V}} \le \Sigma E_{\mathrm{T}} < 70\,{\text {Ge}\text {V}} \).Footnote 7 This is due to the fact that in this range the two muons are the dominant hard objects contributing, with a \(p_{\mathrm{T}}\) resolution proportional to \((p_{\mathrm{T}}^{\mu })^{2}\). A convolution of the muon resolution with a small contribution from \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) is possible for \(\Sigma E_{\mathrm{T}} > 50\,{\text {Ge}\text {V}} \). This component is on average about 60% of \(p_{\mathrm{T}}^{Z}\), and subject to the stochastic fluctuations further discussed below.

The increase of \(Z \rightarrow \mu \mu + 1\) jet topologies in the \(Z \rightarrow \mu \mu \) sample leads to an additional source of fluctuations affecting \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)} (\Sigma E_{\mathrm{T}})\) for \(70\,{\text {Ge}\text {V}} < \Sigma E_{\mathrm{T}} \lesssim 180\,{\text {Ge}\text {V}} \). In general the \(Z \rightarrow \mu \mu \) sample collected for this study covers \(p_{\mathrm{T}}^{Z} \lesssim 140\,{\text {Ge}\text {V}} \) with relevant statistics. At this limit it is expected that the hadronic recoil contains two reconstructed jets, with the onset of this contribution at \(\Sigma E_{\mathrm{T}}\) of about \(180\,{\text {Ge}\text {V}}\). The corresponding change of the dominant final state composition for \(\Sigma E_{\mathrm{T}} > 180\,{\text {Ge}\text {V}} \) leads to a change of shape of \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)} (\Sigma E_{\mathrm{T}})\), as the transverse momentum of the individual jets rises and the number of contributing jets slowly increases. The expected \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)} (\Sigma E_{\mathrm{T}}) \propto \sqrt{\Sigma E_{\mathrm{T}}}\) scaling driven by the jet-\(p_{\mathrm{T}}\) resolution [44] therefore dominates \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) at these higher \(\Sigma E_{\mathrm{T}}\). The MC predictions for \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)} (\Sigma E_{\mathrm{T}})\) agree with the data within a few percent and well within the total uncertainties of this measurement. A tendency for slightly poorer resolution in MC simulations is observed, in particular for \(\Sigma E_{\mathrm{T}} > 200\,{\text {Ge}\text {V}} \).

Any contribution from pile-up to \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) is expected to be associated with the jets. While dedicated corrections applied to the jets largely suppress pile-up contributions in the jet response, residual irreducible fluctuations introduced into the calorimeter signals by pile-up lead to a degradation of the jet energy resolution and thus poorer resolution in the jet-\(p_{\mathrm{T}}\) measurement. The dependence of \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) on the pile-up activity measured by \(N_{\mathrm{PV}}\) is shown in Fig. 7b. Data show a less steep slope of \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)} (N_{\mathrm{PV}})\) than MC simulations, but with about 10% worse resolution in the low pile-up region of \(N_{\mathrm{PV}} \lesssim 5\). The resolution in data is better than in MC simulations by about 10% for the region of higher pile-up activity at \(N_{\mathrm{PV}} \approx 20\).

The differences between data and MC simulations seen in \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)} (\Sigma E_{\mathrm{T}})\) for the inclusive \(Z \rightarrow \mu \mu \) sample can be further analysed by splitting the sample according to the value of \(N_{\mathrm{jet}}\). Figure 8a shows the dependence of \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) on \(\Sigma E_{\mathrm{T}}\) for \(Z \rightarrow \mu \mu \) events with \(N_{\mathrm{jet}} = 0\). The dominant source of fluctuations other than the muon-\(p_{\mathrm{T}}\) resolution is in this case introduced by the incomplete reconstruction of the hadronic recoil. These fluctuations increase with increasing \(p_{\mathrm{T}}^{Z}\), which in turn means higher overall event activity measured by \(\Sigma E_{\mathrm{T}}\). For this sample \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) in data compares well to MC simulations, at a level of a few percent, without any observed dependence on \(\Sigma E_{\mathrm{T}}\).

The exclusive \(N_{\mathrm{jet}} \ge 1\) samples extracted from \(Z \rightarrow \mu \mu \) data and MC simulations show the expected \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)} \propto \sqrt{\Sigma E_{\mathrm{T}}}\) scaling in Fig. 8b. The resolution in data is well represented by MC simulations, at the level of a few percent. The slightly better resolution observed in data with increasing \(\Sigma E_{\mathrm{T}}\) follows the trend observed in Fig. 7a. The similar trends are expected as this kinematic region is largely affected by the jet contribution.

The dependence of \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) on \(N_{\mathrm{PV}}\) shown in Fig. 8c indicates that the \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution is basically independent of pile-up, for the \(N_{\mathrm{jet}} = 0\) sample. This is expected from the exclusive \(E_{\mathrm{T}}^{\mathrm{miss}}\) composition comprising the (track-based) \(E_{\mathrm{T}}^{{\text {miss}},\mu }\) and \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) terms only. Data and MC simulations compare well within a few percent, and without any observable dependence on \(N_{\mathrm{PV}}\). Figure 8d shows the \(N_{\mathrm{PV}}\) dependence of \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) for the \(N_{\mathrm{jet}} \ge 1\) sample. Comparing this result to Fig. 7b confirms that all pile-up dependence of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution is arising from the jet term. Both trend and magnitude of the data-to-MC comparison follow the observation from the inclusive analysis.

6.2.5 \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution in final states with neutrinos

Fig. 9
figure 9

The \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution measured by \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) as a function of the true missing transverse momentum \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\) for the \(W \rightarrow e\nu \), \(W \rightarrow \mu \nu \), and \(t\bar{t}\) samples from MC simulations

The \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution for final states with \(E_{\mathrm{T}}^{{\text {miss}},\text {true}} > 0\) is measured by \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) according to Eq. (12) and evaluated using dedicated inclusive \(W \rightarrow e\nu \) and \(W \rightarrow \mu \nu \) samples from MC simulations, and the inclusive \(t\bar{t}\) MC sample defined in Sect. 5.3. For these samples, \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) can be determined as a function of \(E_{\mathrm{T}}^{{\text {miss}},\text {true}} = p_{\mathrm{T}}^{\nu } \). The dedicated \(W \rightarrow e\nu \) and \(W \rightarrow \mu \nu \) samples are obtained with an event selection based on the description in Sect. 5.2, but omitting both the \(E_{\mathrm{T}}^{\mathrm{miss}}\)-based and the \(m_{\mathrm{T}}\)-based selections.

Figure 9 shows \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) evaluated as a function of \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\) for these samples. The universality of the response to the hadronic recoil observed in Fig. 6, together with the different but subdominant contributions from the \(p_{\mathrm{T}}\) resolutions of the electrons and muons, yield a very similar \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution for \(W \rightarrow e\nu \) and \(W \rightarrow \mu \nu \) final states. Generally, poorer resolution is observed in \(t\bar{t}\) final states. The deviation from the expected \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)} \propto \sqrt{(E_{\mathrm{T}}^{{\text {miss}},\text {true}})}\) scaling behaviour for \(W \rightarrow \ell \nu \) at lower \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\) reflects the kinematic features of the \(W\) boson and its decay. Events with low \(p_{\mathrm{T}}^{W}\), and therefore small hadronic recoil, lie predominantly in the region \(25\,{\text {Ge}\text {V}} \lesssim p_{\mathrm{T}}^{\nu } \lesssim 50\,{\text {Ge}\text {V}} \). Since the hadronic recoil is generally the poorly measured component of an event and the reconstructed \(E_{\mathrm{T}}^{\mathrm{miss}}\) is dominated by the lepton \(p_{\mathrm{T}}\) in this region, the \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution tends to be better here than for events with larger hadronic recoil populating \(p_{\mathrm{T}}^{\nu } \lesssim 25\,{\text {Ge}\text {V}} \) and \(p_{\mathrm{T}}^{\nu } \gtrsim 50\,{\text {Ge}\text {V}} \).

6.3 \(E_{\mathrm{T}}^{\mathrm{miss}}\) tails

Large reconstructed \(E_{\mathrm{T}}^{\mathrm{miss}}\) is an indicator for the production of (potentially new) undetectable particles, but can also be generated by detector problems and/or poor reconstruction of the objects used for its reconstruction. Enhanced tails in the distribution of the \(\mathbf E _{\mathrm{T}}^{{\text {miss}}}\) components for final states with well-known expectation values for \(E_{\mathrm{T}}^{\mathrm{miss}}\) are indicative of such inefficiencies.

Non-Gaussian shapes in the distribution arise from a combination of object selection inefficiencies and potentially non-Gaussian resolutions of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) constituents. Even for a well-defined final state, event-by-event fluctuations in terms of which particles, jets, and soft tracks enter the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction, and with which \(p_{\mathrm{T}}\), lead to deviations from a normally distributed (\(E_{x}^{{\text {miss}}}\),\(E_{y}^{{\text {miss}}}\)) response.

Figure 10 shows the combined (\(E_{x}^{{\text {miss}}}\),\(E_{y}^{{\text {miss}}}\)) distribution for the inclusive \(Z \rightarrow \mu \mu \) sample from MC simulations. To illustrate its symmetric nature and its deviation from a normal distribution in particular with respect to the tails, Gaussian functions are fitted to two limited ranges around the centre of the distribution, \(\pm 1\times \text {RMS}\) and \(\pm 2\times \text {RMS}\). The differences between these functions and the data distribution (lower panel of Fig. 10) indicate a more peaked shape around the most probable value for \(E_{x(y)}^{{\text {miss}}}\) with near exponential slopes. The result of this comparison supports the choice of \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) defined in Eq. (12) in Sect. 6.2.3 for the determination of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution, rather than using any of the widths measured by fitting Gauss functions in selected ranges of the distribution.

The tails in this shape are reflected in the distribution of \(E_{\mathrm{T}}^{\mathrm{miss}}\) itself and can be estimated by measuring the fraction of events with \(E_{\mathrm{T}}^{\mathrm{miss}} > E_{\mathrm{T}}^{{\text {miss}},\text {threshold}} \),

$$\begin{aligned}&f_{\mathrm{tail}} = \dfrac{1}{\mathcal {H}}{\int _{E_{\mathrm{T}}^{{\text {miss}},\text {threshold}}}^{\infty } h(E_{\mathrm{T}}^{\mathrm{miss}}) \text {d}E_{\mathrm{T}}^{\mathrm{miss}}}, \nonumber \\&\quad \text {with }\mathcal {H} = \int _{0}^{\infty } h(E_{\mathrm{T}}^{\mathrm{miss}}) \text {d}E_{\mathrm{T}}^{\mathrm{miss}}. \end{aligned}$$
(13)

Here \(h(E_{\mathrm{T}}^{\mathrm{miss}})\) is the \(E_{\mathrm{T}}^{\mathrm{miss}}\) distribution for a given event sample, and \(E_{\mathrm{T}}^{{\text {miss}},\text {threshold}}\) is a threshold set to estimate tails. Any decrease of \(f_{\mathrm{tail}}\) at a fixed integral \(\mathcal {H}\) indicates an improvement of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution, and is more sensitive to particular improvements than e.g. \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\). For example, improving the \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) reconstruction by rejecting ID tracks from the hard-scatter vertex with poor reconstruction quality yields a significantly smaller \(f_{\mathrm{tail}}\) for the same event sample.

Fig. 10
figure 10

The combined distribution of \(E_{x}^{{\text {miss}}}\) and \(E_{y}^{{\text {miss}}}\) for an inclusive \(Z \rightarrow \mu \mu \) from simulation. Gaussian fits limited to the \(\pm 1\times \text {RMS}\) and \(\pm 2\times \text {RMS}\) ranges around the centre of the distribution are shown, together with the respective differences between the fitted functions and the actual distribution

Fig. 11
figure 11

In a the integral tail fraction \(f_{\mathrm{tail}}\) given in Eq. (13) is shown as a function of the integration threshold \(E_{\mathrm{T}}^{{\text {miss}},\text {threshold}}\), for MC simulations of \(Z \rightarrow \ell \ell \), \(W \rightarrow \ell \nu \), and \(t\bar{t}\) final states. The tail fraction in terms of a threshold applied to \(|\mathbf E _{\mathrm{T}}^{{\text {miss}}}-\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {true}} |\), the distance between the reconstructed (\(\mathbf E _{\mathrm{T}}^{{\text {miss}}}\)) and the expected (\(\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {true}}\) ) vectors, is shown in b for all considered final states. The same fraction is shown in c for the \(E_{\mathrm{T}}^{\mathrm{miss}}\) distributions for \(Z \rightarrow \mu \mu \) before and after a reweighting following the \(\Sigma E_{\mathrm{T}}\) distribution for \(t\bar{t}\) is applied, together with \(f_{\mathrm{tail}}\) from the \(t\bar{t}\) final state

The tails in the \(E_{\mathrm{T}}^{\mathrm{miss}}\) distributions for the final states considered for this study are quantified by the fraction of events above a certain \(E_{\mathrm{T}}^{\mathrm{miss}}\) threshold using MC simulations. Figure 11a shows that the \(Z \rightarrow \ell \ell \) events (\(\ell = e \) or \(\ell = \mu \)) with \(E_{\mathrm{T}}^{{\text {miss}},\text {true}} = 0\) have significantly reduced tails when compared to \(W \rightarrow \ell \nu \) and \(t\bar{t}\) with this metric, and that the tails do not depend on the lepton flavour. A modification of this metric, taking into account \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\) such that the fraction of events with \(|\mathbf E _{\mathrm{T}}^{{\text {miss}}}- \mathbf E _{\mathrm{T}}^{{\text {miss}},\text {true}} |\) above a given threshold is determined, shows the universality of the hadronic recoil in \(Z \rightarrow \ell \ell \) and \(W \rightarrow \ell \nu \), as can be seen in Fig. 11b.

Another finding of this study is that the tail in the \(|\mathbf E _{\mathrm{T}}^{{\text {miss}}}-\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {true}} |\) distribution for the higher \(\Sigma E_{\mathrm{T}}^{\text {jet}}\)\(t\bar{t}\) sample is considerably larger than for the low-\(\Sigma E_{\mathrm{T}}^{\text {jet}}\) samples with \(Z \rightarrow \ell \ell \) or \(W \rightarrow \ell \nu \) final states. As can be seen in Fig. 11c, the tails are much more consistent between \(Z \rightarrow \mu \mu \) and \(t\bar{t}\) samples when the distribution for the \(Z \rightarrow \mu \mu \) sample is reweighted such that it follows the same \(\Sigma E_{\mathrm{T}}^{\text {jet}}\) distribution as the \(t\bar{t}\) sample. The enhanced tails are thus likely introduced by the jet response and multiplicity, which has a residual sensitivity to pile-up.

7 Systematic uncertainties

The systematic uncertainties associated with the measurement of \(E_{\mathrm{T}}^{\mathrm{miss}}\) are provided for the response (\(E_{\mathrm{T}}^{\mathrm{miss}}\) scale) as well as for the resolution. They depend on the composition of the hard terms and on the magnitude of the corresponding soft term. As the hard-term composition is generally defined by optimisations implemented in the context of a given analysis, the contributions of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) terms need to be extracted from the scale and resolution uncertainties for the individual contributing objects comprising electrons, photons, muons, \(\tau \text {-leptons}\), and jets. In the corresponding propagations, correlations between systematic uncertainties for the same type of object are typically taken into account. However, it is assumed that systematic uncertainties of the different object types entering \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction are uncorrelated. The determination of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) scale and resolution uncertainties arising from the soft term \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) is described in this section.

7.1 Methodology

The extraction of the systematic uncertainties for the reconstructed \(E_{\mathrm{T}}^{\mathrm{miss}}\) is based on data-to-MC comparisons of spectra of observables measuring the contribution of \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) to the overall \(E_{\mathrm{T}}^{\mathrm{miss}}\).

7.1.1 Observables

The vector sum of the transverse momentum vectors of all particles and jets emerging from a hard-scatter interaction (\(\mathbf p _{\mathrm{T}}^{\text {HS}}\)) is given by

Here \(\mathbf p _{\mathrm{T}}^{\nu }\) generally represents the transverse momenta of non-observable particles, which are summed up to form \(\mathbf p _{\mathrm{T}}^{\text {inv}}\). All other transverse momenta are carried by particles that are observable in principle, and sum up to \(\mathbf p _{\mathrm{T}}^{\text {obs}}\). Momentum conservation dictates \(p_{\mathrm{T}}^{\text {HS}} = |\mathbf p _{\mathrm{T}}^{\text {HS}} | = 0\).

Due to detector acceptance limitations and inefficiencies in hard-object reconstruction and calibration, and all other effects discussed in Sect. 3, only a proxy (\(\mathbf p _{\mathrm{T}}^{{\text {hard}}}\)) for the observable-particle contribution \(\mathbf p _{\mathrm{T}}^{\text {obs}}\) can be measured. The reconstructed hard final-state objects entering \(E_{\mathrm{T}}^{\mathrm{miss}}\) as described in Sect. 3.2 are used to measure \(\mathbf p _{\mathrm{T}}^{{\text {hard}}}\) as

$$\begin{aligned} \mathbf p _{\mathrm{T}}^{{\text {hard}}}= & {} \sum _{\begin{array}{c} \text {contributing} \\ \text {electrons} \end{array}} \mathbf p _{\mathrm{T}}^{e} + \sum _{\begin{array}{c} \text {contributing} \\ \text {photons} \end{array}} \mathbf p _{\mathrm{T}}^{\gamma } + \sum _{\begin{array}{c} \text {contributing} \\ \tau \text {-leptons} \end{array}} \mathbf p _{\mathrm{T}}^{\tau _{\mathrm{had}}} \\&+\, \sum _{\begin{array}{c} \text {contributing} \\ \text {muons} \end{array}} \mathbf p _{\mathrm{T}}^{\mu } + \sum _{\begin{array}{c} \text {contributing} \\ \text {jets} \end{array}} \mathbf p _{\mathrm{T}}^{\text {jet}}. \end{aligned}$$

The expectation is that \(p_{\mathrm{T}}^{{\text {hard}}} = |\mathbf p _{\mathrm{T}}^{{\text {hard}}} | > 0\) and \(\mathbf p _{\mathrm{T}}^{{\text {hard}}} \ne \mathbf p _{\mathrm{T}}^{\text {obs}} \). Adding \(\mathbf p _{\mathrm{T}}^{{\text {soft}}} = -\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {soft}} \), with \(\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {soft}}\) defined in Eq. (6), to \(\mathbf p _{\mathrm{T}}^{{\text {hard}}}\) yields an improved estimate of the net transverse momentum carried by undetectable particles, as some of the experimental inefficiencies are mitigated.Footnote 8

In the \(Z \rightarrow \mu \mu \) final state without genuine missing transverse momentum the expectation is that \(\mathbf E _{\mathrm{T}}^{{\text {miss}}} = -(\mathbf p _{\mathrm{T}}^{{\text {hard}}} + \mathbf p _{\mathrm{T}}^{{\text {soft}}}) = \mathbf {0}\). While this expectation does not hold due to the experimental inefficiencies, it nevertheless raises the expectation that for events without jets \(\mathbf p _{\mathrm{T}}^{{\text {soft}}}\) points into the direction of the hadronic recoil, i.e. opposite to \(\mathbf p _{\mathrm{T}}^{{\text {hard}}}\) in the transverse-momentum plane. The deviation from this expectation is measured in terms of the parallel (\(\mathcal {P} _{\parallel }\)) and perpendicular (\(\mathcal {P} _{\!\perp }\)) projections of \(\mathbf p _{\mathrm{T}}^{{\text {soft}}}\) onto \(\mathbf p _{\mathrm{T}}^{{\text {hard}}}\). Figure 12 schematically shows these projections for \(Z + 0\)-jet and \(Z + 1\)-jet topologies.

The average \(\langle \mathcal {P} _{\parallel } \rangle \) in a given bin k of phase space defined by \(p_{\mathrm{T}}^{{\text {hard}}}\) measures the \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) response, with \(\langle \mathcal {P} _{\parallel } \rangle = \langle p_{\mathrm{T}}^{{\text {hard}}} \rangle _{k}\) indicating a perfect response in this bin. The \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution contribution from \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) reconstruction is measured by two components, the fluctuations in response (\(\mathrm {RMS}^{2}_{\!\parallel }\)) and the fluctuations of the (transverse) angular deflection around the \(\mathbf p _{\mathrm{T}}^{{\text {hard}}}\) axis, measured by \(\mathrm {RMS}^{2}_{\!\perp }\). These fluctuations are expressed in terms of variances, with

$$\begin{aligned} \mathrm {RMS}^{2}_{\!\parallel } = \left\langle (\mathcal {P} _{\parallel })^{2}\right\rangle - \left\langle \mathcal {P} _{\parallel } \right\rangle ^{2} \quad \text {and}\quad \mathrm {RMS}^{2}_{\!\perp } = \left\langle (\mathcal {P} _{\!\perp })^{2}\right\rangle . \end{aligned}$$
Fig. 12
figure 12

Schematic view of the parallel (\(\mathcal {P} _{\parallel }\)) and perpendicular (\(\mathcal {P} _{\!\perp }\)) projections of \(\mathbf p _{\mathrm{T}}^{{\text {soft}}}\) on \(\mathbf p _{\mathrm{T}}^{{\text {hard}}}\) for \(Z \rightarrow \mu \mu \) events without genuine \(E_{\mathrm{T}}^{\mathrm{miss}}\), for a a final state without any jets and b a final state with one jet. The expectation values for a perfect \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction are \(\text {E}[\mathcal {P} _{\parallel } ] = p_{\mathrm{T}}^{Z} \) for \(N_{\mathrm{jet}} = 0\) and \(\text {E}[\mathcal {P} _{\parallel } ] = p_{\mathrm{T}}^{{\text {hard}}} \) for \(N_{\mathrm{jet}} \ge 1\), with \(\text {E}[\mathcal {P} _{\!\perp } ] = 0\) in all cases. a \(Z + 0\,\text {jet} \) topology. b \(Z + 1\,\text {jet} \) topology

7.1.2 Procedures

The extraction of the systematic uncertainties introduced into the \(E_{\mathrm{T}}^{\mathrm{miss}}\) measurement by the \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) term is based on data-to-MC-simulations comparisons of \(\langle \mathcal {P} _{\parallel } \rangle (p_{\mathrm{T}}^{{\text {hard}}})\) for the response, and of \(\mathrm {RMS}^{2}_{\!\parallel } (p_{\mathrm{T}}^{{\text {hard}}})\) and \(\mathrm {RMS}^{2}_{\!\perp } (p_{\mathrm{T}}^{{\text {hard}}})\) for the resolution. Alternative MC samples are considered, with variations of either the event generator or the detector simulation (description and shower models). For the highest impact of \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) on \(E_{\mathrm{T}}^{\mathrm{miss}}\), the exclusive \(Z \rightarrow \mu \mu \) selection with \(N_{\mathrm{jet}} = 0\) is the basis for the determination of the systematic uncertainty components for both data and all MC simulations. In this case, the only hard contribution is from the reconstructed \(Z\) boson, i.e. \(\mathbf p _{\mathrm{T}}^{{\text {hard}}} = \mathbf p _{\mathrm{T}}^{Z} \) as shown in Fig. 12a.

Fig. 13
figure 13

The a average value of the longitudinal projection \(\langle \mathcal {P} _{\parallel } \rangle \) and the b variance \(\mathrm {RMS}^{2}_{\!\parallel }\) of the longitudinal projection \(\mathcal {P} _{\parallel }\) of \(\mathbf p _{\mathrm{T}}^{{\text {soft}}}\) onto \(\mathbf p _{\mathrm{T}}^{{\text {hard}}}\) for \(Z \rightarrow \mu \mu \) event with \(N_{\mathrm{jet}} = 0\), for data and two different MC simulations, shown as a function of \(p_{\mathrm{T}}^{{\text {hard}}}\). The variance \(\mathrm {RMS}^{2}_{\!\perp }\) of the perpendicular projection \(\mathcal {P} _{\!\perp }\) is shown in c for the same event samples. The shaded band indicates the systematic uncertainties derived as described in the text

The uncertainties are determined by comparing \(\mathcal {P} _{\parallel }\) and \(\mathcal {P} _{\!\perp }\) spectra from data and MC simulations, in bins of \(p_{\mathrm{T}}^{{\text {hard}}}\). For \(\mathcal {P} _{\parallel }\), the smearing of the response and the width both yield scale and (longitudinal) resolution offsets. In the case of \(\mathcal {P} _{\!\perp }\), only smearing of the width is applied to provide transverse resolution offsets. These fitted offsets, determined for the various MC configurations, provide the systematic uncertainties with respect to a specific MC modelling configuration. In practice, to account for the resolution offsets, Gaussian smearing is applied in simulation to the longitudinal and transverse components of \(\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {soft}}\) relative to the direction of \(\mathbf p _{\mathrm{T}}^{{\text {hard}}}\). To account for differences in response between data and simulation, the longitudinal component of \(\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {soft}}\) is scaled up and down to give an uncertainty band.

In order to generate the required number of simulated events, some analyses in ATLAS may have to use the fast detector simulation ATLFAST2 [38, 45] for the calorimeter response. It employs parameterisations for electromagnetic and hadronic showers, instead of the explicit simulation of the particle tracking through matter and the energy-loss mechanisms in a detailed detector geometry. An additional uncertainty is assigned to effects introduced by ATLFAST2. This uncertainty contribution only needs to be considered in analyses using this fast simulation, and does not apply for the results presented in this paper. In analyses where it is applicable, it is added in quadrature to the standard uncertainties.

7.2 Systematic uncertainties in \(E_{\mathrm{T}}^{\mathrm{miss}}\) response and resolution

The result for the systematic uncertainty of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) scale, determined as discussed in the previous section, is summarised in Fig. 13. The average longitudinal projection of \(\mathbf p _{\mathrm{T}}^{{\text {soft}}}\) onto \(\mathbf p _{\mathrm{T}}^{{\text {hard}}}\), \(\langle \mathcal {P} _{\parallel } \rangle \), as a function of \(p_{\mathrm{T}}^{{\text {hard}}}\) is shown in Fig. 13a which compares data to both the standard Powheg +Pythia8-based simulations and the alternative MC simulation employing MadGraph, as described in Sect. 4.2. All MC simulation results are expected to have \(\langle \mathcal {P} _{\parallel } \rangle _{\mathrm{MC}}\) within the uncertainties of the data. The lower panel of Fig. 13a confirms that the ratio \(\langle \mathcal {P} _{\parallel } \rangle _{\mathrm{MC}}/\langle \mathcal {P} _{\parallel } \rangle _{\mathrm{data}}\) lies within the systematic uncertainty band over the full \(p_{\mathrm{T}}^{{\text {hard}}}\) range.

The systematic uncertainty for the \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution is extracted from the variances of the parallel (\(\mathrm {RMS}^{2}_{\!\parallel }\)) and perpendicular (\(\mathrm {RMS}^{2}_{\!\perp }\)) projections of \(\mathbf E _{\mathrm{T}}^{{\text {miss}}}\) onto \(\mathbf p _{\mathrm{T}}^{{\text {hard}}}\) defined in Sect. 7.1.2. Figure 13b shows the \(p_{\mathrm{T}}^{{\text {hard}}}\) dependence of \(\mathrm {RMS}^{2}_{\!\parallel }\) measured for the exclusive \(Z \rightarrow \mu \mu \) sample (\(N_{\mathrm{jet}} = 0\)) in data and two MC simulations. The variances \((\mathrm {RMS}^{2}_{\!\parallel })_{\mathrm{MC}}\) calculated for both sets of simulations agree within the systematic uncertainties of \((\mathrm {RMS}^{2}_{\!\parallel })_{\mathrm{data}}\) with the data, as illustrated in the lower panel of the figure, where the ratio \((\mathrm {RMS}^{2}_{\!\parallel })_{\mathrm{MC}}/(\mathrm {RMS}^{2}_{\!\parallel })_{\mathrm{data}}\) is shown as a function of \(p_{\mathrm{T}}^{{\text {hard}}}\). The results of the evaluation of the variances \(\mathrm {RMS}^{2}_{\!\perp }\) of the perpendicular projections as a function of \(p_{\mathrm{T}}^{{\text {hard}}}\) are shown in Fig. 13c, together with the resulting \(p_{\mathrm{T}}^{{\text {hard}}}\) dependence of the ratio \((\mathrm {RMS}^{2}_{\!\perp })_{\mathrm{MC}}/(\mathrm {RMS}^{2}_{\!\perp })_{\mathrm{data}}\). The systematic uncertainties of the data cover all differences to MC simulations.

8 Missing transverse momentum reconstruction variants

8.1 Calorimeter-based \(E_{\mathrm{T}}^{\mathrm{miss}}\)

The \(E_{\mathrm{T}}^{\mathrm{miss}}\) soft term from the calorimeter \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}},\text {calo}}\) is reconstructed from topo-clusters. As discussed in Ref. [6], each topo-cluster provides a basic EM scale signal as well as a calibrated signal reconstructed using local cell weighting (LCW), and \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}},\text {calo}}\) is calculated from topo-clusters calibrated at the LCW scale. Only topo-clusters with a calibrated energy \(E_{{\text {clus}}}^{\text {LCW}} > 0\), not contributing to the reconstruction of the hard objects used to calculate the hard term given in Eq. (6), are considered for \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}},\text {calo}}\). In addition, topo-clusters that are formed at the same location as the hard object signals are not considered for \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}},\text {calo}}\) even if their signals are not directly contributing to the reconstruction of the hard objects. The fully reconstructed \(E_{\mathrm{T}}^{\mathrm{miss}}\) using \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}},\text {calo}}\) is \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\).

Compared to the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\) and \(\Sigma E_{\mathrm{T}}\), \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) and \(\Sigma E_{\mathrm{T}}^{\text {calo}}\) have an enhanced dependence on pile-up, mostly introduced by the soft term. To partly compensate for the irreducible contribution of \(p_{\mathrm{T}}\)-flow reconstructed from topo-clusters generated by pile-up to \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\), a modified jet selection and ambiguity resolution is applied in their reconstruction. The considered jets are reconstructed following the prescription in Sect. 3.3.5, and required to have a fully calibrated \(p_{\mathrm{T}} > 20\,{\text {Ge}\text {V}} \). The contribution of these jets to \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) and \(\Sigma E_{\mathrm{T}}^{\text {calo}}\), defined in terms of momentum components \((p_{x},p_{y})\), depends on the overlap with already accepted reconstructed particles,

$$\begin{aligned} (p_{x},p_{y}) = \left\{ \begin{array}{lll} (0,0) &{} \kappa _{E} \ge 50\% &{} \text {(large overlap)} \\ (1-\kappa _{E})\times (p_{x}^{\text {jet}},p_{y}^{\text {jet}}) &{} \kappa _{E} < 50\% &{} \text {(small or no overlap)} \\ \end{array} \right. . \end{aligned}$$
(14)

The overlap fraction \(\kappa _{E}\) is given in Eq. (8). Jets with \(\kappa _{E} \ge 50\%\) are not used at all. The JVT-based tagging of non-pile-up jets is omitted. It is found that this strategy reduces the fluctuations in the \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) reconstruction. The transverse momentum contribution of groups of clusters representing a jet-like \(p_{\mathrm{T}}\)-flow e.g. from pile-up in a given direction that are not reconstructed and calibrated as a jet, or do not pass the jet-\(p_{\mathrm{T}}\) threshold applied in \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction, is reduced if all jets and jet fragments, including those from pile-up, are included.

8.2 \(E_{\mathrm{T}}^{\mathrm{miss}}\) from tracks

The reference track-based soft term \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) is largely insensitive to pile-up, as indicated by the dependence of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) on \(N_{\mathrm{PV}}\) in the exclusive \(Z \rightarrow \mu \mu \) sample (\(N_{\mathrm{jet}} = 0\)) shown in Fig. 7c. As discussed in Sect. 6.2.4 and from the comparison of Fig. 7c, d, the pile-up dependence of \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\) in the inclusive \(Z \rightarrow \mu \mu \) sample is largely introduced by the jet contribution. This contribution suffers from (1) the lack of pile-up suppression for forward jets with \(|\eta | > 2.4\), (2) any inefficiency connected with the JVT-based tagging, and (3) irreducible pile-up-induced fluctuations in the calorimeter jet signals. Using a representation of \(E_{\mathrm{T}}^{\mathrm{miss}}\) employing only reconstructed ID tracks from the primary vertex increases stability against pile-up as long as the tracking and vertex resolution is not affected by it. In this representation (\(p_{\mathrm{T}}^{\mathrm{miss}}\)) all jets and reconstructed particles are ignored, i.e. the \(p_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction does not include any calorimeter or MS signals. The \(p_{\mathrm{T}}^{\mathrm{miss}}\) resolution is then inherently immune to pile-up, while the \(p_{\mathrm{T}}^{\mathrm{miss}}\) response is low as all neutral \(p_{\mathrm{T}}\)-flow in \(|\eta | < 2.5\) as well as all \(p_{\mathrm{T}}\)-flow outside of this region is excluded.

8.3 Performance evaluations for \(E_{\mathrm{T}}^{\mathrm{miss}}\) variants

The main motivation to study \(E_{\mathrm{T}}^{\mathrm{miss}}\)-reconstruction variants is to improve some combination of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution, scale, and stability against pile-up. As with the composition of objects entering \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction in general, the particular choice of variant used for a given analysis strongly depends on the performance requirements for this analysis. The comparison of both the resolution and response of \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) and \(p_{\mathrm{T}}^{\mathrm{miss}}\) to the corresponding measurements using the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\) illustrates their principal features for the \(Z \rightarrow \mu \mu \) and \(t\bar{t}\) production final state.

Fig. 14
figure 14

Comparison of the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution with the resolutions of the track-only-based variant \(p_{\mathrm{T}}^{\mathrm{miss}}\) described in Sect. 8.2, and the reconstruction variant \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) employing a calorimeter-based soft term, as discussed in Sect. 8.1. The resolutions are determined as described in Sect. 6.2.3 and shown as a function of the \(\Sigma E_{\mathrm{T}}\). For consistency, for all three variants, the \(\Sigma E_{\mathrm{T}}\) value is taken from \(E_{\mathrm{T}}^{\mathrm{miss}}\)

Fig. 15
figure 15

Comparison of the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution with the resolutions of the track-only-based variant \(p_{\mathrm{T}}^{\mathrm{miss}}\) described in Sect. 8.2, and the reconstruction variant \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) employing a calorimeter-based soft term, as discussed in Sect. 8.1. The resolutions are determined as described in Sect. 6.2.3 and shown as a function of the pile-up activity measured in terms of the number of reconstructed vertices \(N_{\mathrm{PV}}\) for a an exclusive \(Z \rightarrow \mu \mu \) sample without jets with \(p_{\mathrm{T}} > 20\,{\text {Ge}\text {V}} \) and b an inclusive \(Z \rightarrow \mu \mu \) sample, both selected from data. In c, the resolution of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction-variants in a final state with significant jet activity and \(p_{\mathrm{T}}^{\nu } > 0\) is compared using MC simulations of \(t\bar{t}\) production

8.3.1 Comparisons of \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution

Figure 14 compares the \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) and \(p_{\mathrm{T}}^{\mathrm{miss}}\) resolutions with the one obtained from the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\), for the inclusive \(Z \rightarrow \mu \mu \) sample in data. Each is shown as a function of \(\Sigma E_{\mathrm{T}}\) corresponding to the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\), giving an estimate of the total hard-scatter activity. The low-\(\Sigma E_{\mathrm{T}}\) region is dominated by events with \(N_{\mathrm{jet}} = 0\), where the contribution of \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}},\text {calo}}\) in \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) yields a poorer resolution than for \(E_{\mathrm{T}}^{\mathrm{miss}}\), and where \(E_{\mathrm{T}}^{\mathrm{miss}}\) and \(p_{\mathrm{T}}^{\mathrm{miss}}\) have identical performance. The high-\(\Sigma E_{\mathrm{T}}\) region is dominated by events with higher jet multiplicity, where \(p_{\mathrm{T}}^{\mathrm{miss}}\) resolution is degraded relative to the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\) by the incomplete measurement of jets.

Figure 15a compares the \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) and \(p_{\mathrm{T}}^{\mathrm{miss}}\) resolution as functions of the pile-up activity measured by \(N_{\mathrm{PV}}\), with the one obtained from the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\) for the exclusive \(Z \rightarrow \mu \mu \) samples with \(N_{\mathrm{jet}} = 0\) in data. The \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) resolution is dominated by pile-up and shows significantly degraded performance relative to \(p_{\mathrm{T}}^{\mathrm{miss}}\) and the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\). The exclusive use of only tracks from the hard-scatter vertex for both \(p_{\mathrm{T}}^{\mathrm{miss}}\) and \(E_{\mathrm{T}}^{\mathrm{miss}}\) yields the same stability against pile-up.

In events with jet activity, the degraded \(p_{\mathrm{T}}^{\mathrm{miss}}\) resolution is observable, especially outside the region of highest pile-up activity, as seen in Fig. 15b for the \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution obtained with the inclusive \(Z \rightarrow \mu \mu \) sample in data for \(N_{\mathrm{PV}} \lesssim 15\). This is even more obvious in final states with relatively high jet multiplicity and genuine missing transverse momentum, like for the \(t\bar{t}\)-production sample from MC simulations. As shown in Fig. 15c for this final state, both the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\) and the calorimeter-based \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) have a significantly better resolution than \(p_{\mathrm{T}}^{\mathrm{miss}}\), at the price of some sensitivity to pile-up, which is absent for \(p_{\mathrm{T}}^{\mathrm{miss}}\). The \(N_{\mathrm{PV}}\) dependence of the resolution is enhanced in \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\), due to the increased contribution from soft calorimeter signals without pile-up suppression at higher \(N_{\mathrm{PV}}\).

8.3.2 Comparisons of \(E_{\mathrm{T}}^{\mathrm{miss}}\) scale

Fig. 16
figure 16

Comparison of the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\), the calorimeter-based \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) and track-only-based \(p_{\mathrm{T}}^{\mathrm{miss}}\) response in an a exclusive and an b inclusive \(Z \rightarrow \mu \mu \) sample from data. The projections of the respective \(\mathbf E _{\mathrm{T}}^{{\text {miss}}}\), \(\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {calo}}\), and \(\mathbf p _{\mathrm{T}}^{\mathrm{miss}}\) onto the direction of \(\mathbf p _{\mathrm{T}}^{Z}\), calculated according to Eqs. (9) and (10), are shown as a function of \(p_{\mathrm{T}}^{Z}\). In c, the linearity of the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\), \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\), and \(p_{\mathrm{T}}^{\mathrm{miss}}\) scales, calculated according to Eq. (11), is shown as a function of the true \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\) for the \(t\bar{t}\)-production MC simulation sample

Following the description in Sect. 6.2.1, the \(E_{\mathrm{T}}^{\mathrm{miss}}\) response is evaluated for the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\), \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\), and \(p_{\mathrm{T}}^{\mathrm{miss}}\) using the respective projections of \(\mathbf E _{\mathrm{T}}^{{\text {miss}}}\), \(\mathbf E _{\mathrm{T}}^{{\text {miss}},\text {calo}}\), and \(\mathbf p _{\mathrm{T}}^{\mathrm{miss}}\) onto the direction of \(\mathbf p _{\mathrm{T}}^{Z}\), according to Eqs. (9) and (10). Figure 16a shows the average projection as a function of \(p_{\mathrm{T}}^{Z}\) for the exclusive \(Z \rightarrow \mu \mu \) sample with \(N_{\mathrm{jet}} = 0\) in data. Both \(E_{\mathrm{T}}^{\mathrm{miss}}\) and \(p_{\mathrm{T}}^{\mathrm{miss}}\) show the same increasingly incomplete reconstruction of the hadronic recoil in this sample for rising \(p_{\mathrm{T}}^{Z}\). This reconstruction is slightly improved for \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\), but still insufficient at higher \(p_{\mathrm{T}}^{Z}\).

In the inclusive \(Z \rightarrow \mu \mu \) sample, shown in Fig. 16b, the indication at lower \(p_{\mathrm{T}}^{Z}\) is that \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) has a higher response and thus a better representation of the hadronic recoil, due to the more complete \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}}}\) reconstruction and the lack of a JVT-tagging requirement. This effect is partly due to the observation bias in the response introduced by the relatively poor \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) resolution, as discussed in Sect. 3.1. Both \(E_{\mathrm{T}}^{\mathrm{miss}}\) and \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) show comparable response for \(p_{\mathrm{T}}^{Z} \gtrsim 60\,{\text {Ge}\text {V}} \), owing to the JVT cut-off at \(60\,{\text {Ge}\text {V}} \). The slightly larger \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) response of about \(1\,{\text {Ge}\text {V}}\) reflects the contribution from neutral signals to the soft term. The degraded response associated with \(p_{\mathrm{T}}^{\mathrm{miss}}\) related to the exclusion of hard objects is clearly visible in this figure.

Figure 16c shows the linearity of the various \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction approaches as a function of \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\) for the \(t\bar{t}\)-production sample from MC simulations. Beyond \(E_{\mathrm{T}}^{{\text {miss}},\text {true}} \approx 60\,{\text {Ge}\text {V}} \) both the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\) and \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) show the same good linearity, while the lack of a jet contribution to \(p_{\mathrm{T}}^{\mathrm{miss}}\) shows a loss of response up to about 50% at higher \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\). The overestimation of \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\) by all three reconstruction variants at lower \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\) reflects the observation bias in the response introduced by the resolution. The poorer resolution associated with \(p_{\mathrm{T}}^{\mathrm{miss}}\) observed in Fig. 15c for this sample leads to a faster rise of the response with decreasing \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\) than for the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\) and \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\), which show a very similar dependence on \(E_{\mathrm{T}}^{{\text {miss}},\text {true}}\).

8.3.3 Summary of performance

Both \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) and \(p_{\mathrm{T}}^{\mathrm{miss}}\) offer alternative measures for \(E_{\mathrm{T}}^{\mathrm{miss}}\). The calorimeter-based \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) uses topo-clusters calibrated at the LCW scale for the soft term, which are neither part of the signal nor otherwise overlapping with the signals of other hard objects contributing to \(E_{\mathrm{T}}^{\mathrm{miss}}\). This introduces a pile-up dependence into \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\), due to the lack of pile-up suppression of calorimeter signals outside of reconstructed hard objects. It features a slightly modified jet contribution without the JVT-based selection used in case of the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction, to allow the cancellation of jet-like \(p_{\mathrm{T}}\)-flow from pile-up in \(E_{\mathrm{T}}^{{\text {miss}},{\text {soft}},\text {calo}}\) by pile-up jets in its hard term. The \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) response in the inclusive \(Z \rightarrow \mu \mu \) sample is better than the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\) response, in particular in the region of small hadronic recoil (\(p_{\mathrm{T}}^{Z} \lesssim 20\,{\text {Ge}\text {V}} \)). It is comparable to the reference in \(t\bar{t}\) final states. The observed \(\mathrm {RMS}^{\mathrm{miss}}_{x(y)}\), in particular in \(Z \rightarrow \mu \mu \) without jets, is significantly more affected by pile-up than is the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\) or the track-only-based \(p_{\mathrm{T}}^{\mathrm{miss}}\). In final states with a considerable number of jets, like \(t\bar{t}\), \(E_{\mathrm{T}}^{{\text {miss}},\text {calo}}\) performs nearly as well as the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\), with a slight degradation of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution at highest pile-up activities. This variant is useful for physics analyses least sensitive to the soft-term contribution to \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution but requiring a linear \(E_{\mathrm{T}}^{\mathrm{miss}}\) response.

The track-only-based \(p_{\mathrm{T}}^{\mathrm{miss}}\) displays a degraded response for the inclusive \(Z \rightarrow \mu \mu \) sample, which is expected from the exclusive use of hard-scatter-vertex tracks. As expected, resolution is not affected by pile-up in the considered final states, but is poorer than, or at most comparable to, the reference \(E_{\mathrm{T}}^{\mathrm{miss}}\) algorithm. Nevertheless, \(p_{\mathrm{T}}^{\mathrm{miss}}\) provides a stable observable for event and phase-space selections in analyses sensitive to \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution.

9 Conclusion

The performance and features of the missing transverse momentum reconstruction in \( pp \) collision data at the LHC, acquired in 2015 with the ATLAS detector at \(\sqrt{s} = 13\,{\text {Te}\text {V}} \) and corresponding to about \(3.2~\hbox {fb}^{-1}\), are evaluated for selected event samples with (\(W \rightarrow e\nu \), \(W \rightarrow \mu \nu \), \(t\bar{t}\)) and without (\(Z \rightarrow \mu \mu \)) genuine \(E_{\mathrm{T}}^{\mathrm{miss}}\). The comparison of the data from the detector with the corresponding MC simulations generally yields good agreement in the covered phase space. The systematic uncertainty contribution from the soft event to the reconstructed \(E_{\mathrm{T}}^{\mathrm{miss}}\) is determined with \(Z \rightarrow \mu \mu \) final states without jets. It is calculated from the data-to-MC-simulations comparison of the parallel and perpendicular projections of the missing transverse momentum vector \(\mathbf E _{\mathrm{T}}^{{\text {miss}}}\) onto the vector sum of the transverse momenta of the hard objects \(\mathbf p _{\mathrm{T}}^{{\text {hard}}}\). The parallel projections yield the uncertainty of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) scale, evaluated as a function of the total transverse momentum of the hard objects (\( p_{\mathrm {T}}^{\mathrm {hard}}\)). The widths of the distributions of the parallel and perpendicular projections yield the respective systematic uncertainties of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution. Simulation tends to underestimate the perpendicular resolution and overestimate the scale and parallel resolution, in each case differing from data by at most 10%.

The performance evaluation of \(E_{\mathrm{T}}^{\mathrm{miss}}\) response and resolution for the inclusive \(Z \rightarrow \mu \mu \) sample shows that data and MC simulations agree within the systematic uncertainties. The \(E_{\mathrm{T}}^{\mathrm{miss}}\) response shows an underestimation of the soft contributions to \(E_{\mathrm{T}}^{\mathrm{miss}}\). A degradation of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) resolution is observed for increasing \(\Sigma E_{\mathrm{T}}\) and \(N_{\mathrm{PV}}\), due to pile-up and detector resolution effects. Additional performance measures considered in these studies include the estimate of tails in the \(E_{\mathrm{T}}^{\mathrm{miss}}\) distribution. As expected from the universality of the hadronic recoil, the integral tail fraction of the \(E_{\mathrm{T}}^{\mathrm{miss}}\) distribution is identical for inclusive \(Z\) and \(W\) boson production, independent of the leptonic decay mode. The \(t\bar{t}\) final states feature a higher jet multiplicity and show larger tails reflecting a higher sensitivity to residual pile-up surviving in the jet contribution to \(E_{\mathrm{T}}^{\mathrm{miss}}\), in terms of the inclusion of pile-up jets as well as the increased fluctuations of the jet response introduced by pile-up.

From the performance studies presented in this paper, the object-based \(E_{\mathrm{T}}^{\mathrm{miss}}\) reconstruction in ATLAS, which was developed for LHC \(\text {Run}\,1\) and used in a large number of physics analyses, can be used with the discussed refinements and adjustments for \(\text {Run}\,2\) as well.