Regular Article -Experimental Physics
Performance ofmissing transverse momentum reconstruction with theATLASdetector usingproton–proton collisions
√
at s = 13TeV
ATLAS Collaboration
CERN, 1211 Geneva 23, Switzerland Received: 23 February 2018 / Accepted: 27 September 2018 / Published online: 8 November 2018
© CERN for the benefit of the ATLAS collaboration 2018
Abstract The performance of the missing transverse momentum (Emiss) reconstruction with the ATLAS detector
T is evaluated using data collected in proton–proton collisions at the LHC at a centre-of-mass energy of 13 TeV in 2015. To reconstruct Emiss, fully calibrated electrons, muons, pho-
T tons, hadronically decaying τ -leptons, and jets reconstructed from calorimeter energy deposits and charged-particle tracks are used. These are combined with the soft hadronic activity measured by reconstructed charged-particle tracks not associated with the hard objects. Possible double counting of contributions from reconstructed charged-particle tracks from the inner detector, energy deposits in the calorimeter, and reconstructed muons from the muon spectrometer is avoided by applying a signal ambiguity resolution procedure which rejects already used signals when combining the various ETmiss contributions. The individual terms as well as the overall reconstructed ETmiss are evaluated with various performance metrics for scale (linearity), resolution, and sensitivity to the data-taking conditions. The method developed to determine the systematic uncertainties of the Emiss
T scale and resolution is discussed. Results are shown based on the full 2015 data sample corresponding to an integrated luminosity of 3 2fb−1.
Contents
1
Introduction ..................... 2
2
ATLASdetector ................... 2 Emiss
3 reconstruction ................. 2
T Emiss
3.1 basics ................... 3
T
Emiss3.2 terms ................... 4
T 3.3Objectselection ................. 4
3.3.1Electronselection ............. 4
3.3.2Photonselection ............. 6
3.3.3 τ
-Leptonselection ............ 6
3.3.4Muonselection .............. 6
* e-mail: atlas.publications@cern.ch
3.3.5Jetselection ................ 6
3.3.6 Muonoverlapwithjets .......... 6 Emiss
3.4 softterm ................. 7
T
3.4.1 Track and vertex selection ........ 8
3.4.2Tracksoftterm .............. 8
4 Dataandsimulationsamples ............. 8
4.1Datasamples .................. 8
4.2MonteCarlosamples .............. 8
4.3Pile-up ...................... 9
5 Eventselection .................... 9
5.1
Z →
ͣ
eventselection ............ 9
5.2
W →
eυ
eventselection ............ 9
5.3 tt¯eventselection ................ 10 6 Performance of ETmiss reconstruction in data and MonteCarlosimulation ............... 10 Emiss
6.1 T modelling in Monte Carlo simulations .. 11 Emiss
6.2 T responseandresolution .......... 11
6.2.1 ETmiss scale determination ......... 15
6.2.2 Measuring the Emissresponse ...... 15
T
6.2.3 Determination of the ETmiss resolution .. 17
6.2.4 ETmiss resolution measurements ...... 17
6.2.5 ETmiss resolution in final states with neutrinos .18 Emiss6.3 tails .................... 18
T 7 Systematicuncertainties ............... 20 7.1Methodology .................. 21
7.1.1Observables ................ 22
7.1.2Procedures ................ 22
7.2 Systematic uncertainties in Emiss response and
T resolution .................... 23
8
Missing transverse momentum reconstruction variants 23
8.1 Calorimeter-based ETmiss ............ 23 Emiss
8.2 fromtracks ................ 25
T
8.3 Performance evaluations for ETmiss variants ... 25
8.3.1 Comparisons of Emiss resolution ..... 25
T
8.3.2 Comparisons of Emissscale ....... 26
T
8.3.3 Summary of performance ......... 28
9 Conclusion ...................... 28
AppendixA:Glossaryofterms ............. 29
123
Appendix B: Alternative Emiss composition ...... 29
T AppendixC:Jetselection ................ 30
References........................ 32
1 Introduction
The missing transverse momentum (Emiss) is an important
T observable serving as an experimental proxy for the trans-verse 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 Emiss is challenging because
T 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 Emiss reconstruction focussing
T on the minimisation of effects introduced by pile-up for the data recorded between 2010 and 2012 (LHC Run1) [2,3]. These approaches are the basis for the Emiss reconstruction
T developed for the data collected in 2015 (LHC 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 Emiss
T 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 selec
tion is outlined in Sect. 5, followed by results for Emiss per-
T formance in Sect. 6. Section 7 comprises a discussion of methods used to determine systematic uncertainties associated with the Emiss measurement, and the presentation of
T the corresponding results. Section 8 describes variations of the Emiss
T 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 Emiss-related variables and descriptors can be
T found in Appendix A, while the composition of Emiss recon-
T struction variants is presented in Appendix B. An evaluation
of the effect of alternative jet selections on the Emiss recon-
T struction 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) coverage in solid angle.1 It consists of an inner detector (ID) tracking system in a 2 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 ||Ϫ 2 5, and consists of a silicon pixel detector, a silicon microstrip detector and a transition radiation tracker for ||Ϫ 2 0. During the LHC shutdown between Run 1 and 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 || Ϫ 3 2. The regions || Ϫ 137 and 15 Ϫ || Ϫ 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 ||Ϫ 1 7. LAr technology is also used for the hadronic calorimeters in the endcap region 1 5 Ϫ ||Ϫ 3 2 and for electromagnetic and hadronic energy measurements in the forward calorimeters covering 3 2 Ϫ ||Ϫ 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 ||= 2 7, and fast detectors for triggering in the region ||Ϫ 2 4.
A two-level trigger system is used to select events [5]. A low-level hardware trigger reduces the data rate, and a highlevel software trigger selects events with interesting final states. More details of the ATLAS detector can be found in Ref. [1].
3 Emiss
reconstruction
T
The reconstructed Emiss in ATLAS is characterised by two
T 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, -leptons, and muons. While muons are recon
1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upwards. Cylindrical coordinates r ψ are used in the transverse plane, ψ being the azimuthal angle around the z-axis. The pseudorapidity is defined in terms of the polar angle ω as χ =−ln tan2 . Angular distance is measured in units of
R ↓χ 2 +ψ 2.
123
structed from ID and MS tracks, electrons and τ -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 ETmiss 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 Emiss. This issue is addressed by the explicit
T signal ambiguity resolution in the object-based Emiss recon-
T struction 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 Emiss
are available and discussed in detail in Sect. 8.One
T of these alternative options is the calorimeter-based Emiss
T 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 Emiss
T only in the use of tracks in place of jets. It is described in
more detail in Sect. 8.2.
3.1 Emiss
basics
T
The missing transverse momentum reconstruction provides a set of observables constructed from the components pxy of the transverse momentum vectors (pT) of the various contributions. The missing transverse momentum components
Emiss
xy serve as the basic input for most of these observables. They are given by
Emiss
xy =−pxyi −pxyj (1) i↔{hard objects} j↔{soft signals}
The set of observables constructed from Emiss is
xy
Emiss = Emiss Emiss
(2)
T xy
J
Emiss =|Emiss|= Emiss 2 + Emiss 2 (3)
TT xy
miss EmissEmiss
= tan−1(4)
yx
The vector Emiss provides the amount of the missing trans-
T verse momentum via its magnitude Emiss, and its direction
T in the transverse plane in terms of the azimuthal angle miss.
Consequently, Emiss is non-negative by definition. However,
T in an experimental environment where not all relevant pT
from the hard-scatter interaction can be reconstructed and used in Eq. (1), and the reconstructed pT from each contribution is affected by the limited resolution of the detector, an observationbias towards non-vanishing values for Emiss
T is introduced even for final states without genuine missing transverse momentum generated by undetectable particles. The scalar sum of all transverse momenta (pT =|pT|) from the objects contributing to Emiss reconstruction is given
T by
L ET =pT i +pT j (5) i↔{hard objects} j↔{soft signals}
In the context of Emiss reconstruction, L ET is calculated in
T addition to the sum given in Eq. (1), and the derived quanti
ties defining Emiss given in Eqs. (2)–(4). It provides a useful
T 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 Emiss reconstruc-
T tion performance studies.
In the calculation of Emiss and L ET the contributing
xy
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 Emiss
T reconstruction sequence for the hard-object contribution starts with electrons (e), followed by photons (γ ), then hadronically decaying τ -leptons (had), and finally jets. Muons () 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 ETmiss reconstruction first. The lowerpriority reconstructed particles (γ , had) are fully rejected if they share their calorimeter signal with a higher-priority object that has already entered the Emiss reconstruction.
T Muons experience energy loss in the calorimeters, but only non-isolated muons overlap with other hard objects, most likely jets or τ -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 Emiss
T reconstruction in the case of partial or marginal overlap, and
123
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 Emiss, ID tracks from the hard-
T scatter collision vertex are used to construct the soft-event
signal for the results presented in this paper.
3.2 Emiss
terms
T
Particle and jet selections in a given analysis should be reflected in Emiss and L ET for a consistent interpretation of
T 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 Emiss
T and L ET, 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 Emiss. The applied selec-
T tions are presented in Sect. 3.3, and summarised in Table 1.
In ATLAS the flexibility needed to recalculate Emiss and
T L ET under changing analysis requirements for the same
event is implemented using dedicated variables correspond
ing to specific object contributions. In this approach the full Emiss
T is the vectorial sum of missing transverse momentum terms EmissT p, with p ↔{e had ͣ jet}reconstructed from the pT =( px py of accepted particles and jets, and the corresponding soft term Emiss soft from the soft-event sig-
T nals introduced in Sect. 3.1 and further specified in Sect. 3.4.
This yields2
Emiss e had
=− p− p − p − p −
T TTT T
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, L ET is calculated from the scalar sums of the transverse momenta of hard objects entering the Emiss
T reconstruction and the soft term,
L ET e had ͣ jet
= pT + pT + p + pT + p
TT selected accepted accepted selected acceptedelectrons photons τ -leptons muons jets
hardterm
track
+ pT (7)
unused
tracks
softterm
The hard term in both ETmiss and L ET 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 hardscatter 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
jet track
p− p(6)
TT
selected accepted accepted selected accepted unused electrons photons τ -leptons muons jets tracks
Emiss e Emiss ͣ Emiss soft
Emiss γ miss had Emiss jet
T TT
E
TT T
softtermhardterm
The Emiss
T and miss observables can be constructed accordingtoEqs.(3) and (4), respectively, for the overall missing transverse momentum (from Emiss) as well as for each
T
individual term indicated in Eq. (6). In the priority-ordered reconstruction sequence for Emiss, contributions are defined
T
by a combination of analysis-dependent selections and a possible rejection due to the applied signal ambiguity resolution.
2 In this formula the notion of selected, which is only applicable to electrons and muons, means that the choice of reconstructed particles is purely given by a set of criteria similar to those given in Sects. 3.3.1 and 3.3.4, respectively, with possible modifications imposed by a given analysis. The notion of accepted indicates a modification of the set of initially selected objects imposed by the signal ambiguity resolution.
in Sects. 6–8. Generally, these selections require refinements to achieve optimal Emiss reconstruction performance in the
T
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
123
Table1 Overview of the contributions to Emiss and L ET from hard
T
objects such as electrons (e), photons (γ ), hadronically decaying τ -leptons (had), muons (), 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 Pof consideration for Emiss reconstruction, with
T
(1) being the first and (5) being the last calculated hard-object contri
bution. The soft-event contribution is constructed at the lowest priority (6), after all hard objects are considered. The transverse (longitudinal) impact parameter d0(z0 sinω ) used to select the ID tracks contributing
to Emiss soft and L Esoft
T Tin P = 6 is measured relative to the hardscatter vertex. All variables are explained in Sect. 3.2. The angular distance Rbetween objects is defined as R = χ 2 + ψ 2
P Objects contributing to Emiss and L ET
T
Type Selections Variables Comments
(1) e || Ϫ 1 37 or 1 52 Ϫ || Ϫ 2 47 pT Ϭ 10 GeV Emiss e T L Ee T All e
(2) γ || Ϫ 1 37 or 1 52 Ϫ || Ϫ 2 47 pT Ϭ 25 GeV Emiss γ T L Eγ T
(3) had || Ϫ 1 37 or 1 52 Ϫ || Ϫ 2 47 pT Ϭ 20 GeV Emiss had T L Ehad T
(4) ͣ || Ϫ 2 7 pT Ϭ 10 GeV Emiss ͣ T L Eͣ T
(5) Jet ||
Ϫ
4 5 pT Ϭ
60 GeV Emiss jet T L
Ejet T Sect. 3.3.6
or
2 4 Ϫ || Ϫ 4 5 20 GeV Ϫ pT Ϫ 60 GeV
or
|| Ϫ 2 4 20 GeV Ϫ pT Ϫ 60 GeV JVT Ϭ 0 59
(6) ID track pT Ϭ 400 MeV Emiss soft T
|d0| Ϫ 1 5mm |z0 sinω | Ϫ 1 5mm Rtrack e − γ cluster Ϭ 0 05 Rtrack had Ϭ 0 2 L Esoft T
± passing medium reconstruction quality and kinematic selections
All γ passing tight quality and kinematic selections in reconstruction, and without signal overlap with (1)
All had passing medium reconstruction quality and kinematic selections, and without signal overlap with (1) and (2)
All ͣ passing medium quality and kinematic selections in reconstruction; for the discussion of the –jet overlap removal see
All jets passing reconstruction quality (jet cleaning) and kinematic selections, and without signal overlapa with (1)–(3); for the dedicated overlap removal strategy with ͣ from (4) see Sect. 3.3.6
All ID tracks from the hard-scatter vertex passing reconstruction quality and kinematic selections, and not associated with any particle from (1), (3) or (4), or ghost-associated with a jet from (5)
aWhile for single reconstructed particles no overlap is accepted at all, jets with a signal overlap fraction E Ϫ 50% can still contribute their associated tracks to Emiss soft if those pass the selections for P = 6 , as discussed in Sect. 3.3.5. The definition of E is giveninEq. (8)
T
Ref. [7]. To be considered for Emiss reconstruction, elec-deposit by electrons within 1 37 Ϫ || Ϫ 1 52 is likely recon-
T
trons passing the reconstruction quality requirements are in structed as a jet and enters the Emiss reconstruction as such, if
T
addition required to have pT Ϭ 10 GeV and || Ϫ 1 37 or this jet meets the corresponding selection criteria discussed 1 52 Ϫ || Ϫ 2 47, to avoid the transition region between the in Sect. 3.3.5. central and endcap electromagnetic calorimeters. Any energy
123
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 pT Ϭ 25 GeV and || Ϫ 137 or 152 Ϫ || Ϫ 2 37 to be included in the Emiss reconstruction. Similarly to electrons,
T photons emitted within 1 37 Ϫ || Ϫ 1 52 may contribute to Emiss
T as a jet.
3.3.3 τ -Lepton selection
Hadronically decaying τ -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 pT Ϭ 20 GeV and || Ϫ 1 37 or 1 52 Ϫ || Ϫ 2 47. Any τ -lepton not satisfying these τ identification criteria may contribute to Emiss when passing
T the jet selection.
3.3.4 Muon selection
Muons are reconstructed within || Ϫ 2 5 employing a combined MS and ID track fit. Outside of the ID coverage, muons are reconstructed within 2 5 Ϫ || Ϫ 2 7 from a track fit to MS track segments alone. Muons are further selected for
Emiss
T reconstruction by requiring the medium reconstruction quality defined in Ref. [12], pT Ϭ 10 GeV, and an association with the hard-scatter vertex for those within || Ϫ 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.3 The anti-kt algorithm [13], as provided by the FastJet toolkit [14], is employed with a radius param-eter 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 pT Ϭ 20 GeV after the full calibration. The jet contribution to
Emiss
T and L ET 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
3 On this scale the energy deposited in the calorimeter by electrons and photons is represented well. The hadron signal at the EM scale is not corrected for the non-compensating signal features of the ATLAS calorimeters.
variable JVT, ranging from 0 (pile-up-like) to 1 (hard-scatterlike), for each jet with matched tracks.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 E of the electron(photon) energy EEM to the jet energy EEM,
eγ jet
EEM E = eγ (8)
EEM
jet
with both energies calibrated at the EM scale. In the case of E ↕ 50%, the jet is included in Emiss reconstruction,
T
with its pT scaled by 1 − E.For 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 pT Ϭ 20 GeV and || Ϭ 2 4, or with pT ↖ 60 GeV and || Ϫ 4 5, are always accepted for Emiss reconstruction.
T
Jets reconstructed with 20 GeV Ϫ pT Ϫ 60 GeV and || Ϫ 2 4 are only accepted if they are tagged by JVT Ϭ 0 59. In both cases, the jet pT thresholds are applied to the jet pT before applying the E correction. Additional configurations for selecting jets used in Emiss reconstruction are discussed in
T
Appendix A, together with the effect of the variation of these selection criteria on the Emiss reconstruction performance.
T
3.3.6 Muon overlap with jets
Jets overlapping with a reconstructed muon affect the Emiss
T
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 hardscatter interaction due to the muon ID track, they need to be rejected for Emiss reconstruction. On the other hand, jets
T
reconstructed from final state radiation (FSR) off the muon need to be included into ETmiss reconstruction.
In all cases, the muon–jet overlap is determined by ghostassociating 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
4 In the calculation of JVT the total amount of pT carried by tracks from the hard-scatter vertex matched to the given jet is related to the total amount of pT carried by all matched tracks, among other inputs, to tag jets from the hard-scatter interaction.
123
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 Emiss reconstruction and rejects jets
T 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 pT from the hard-scatter vertex and thus increases JVT. As a consequence of this fake tag, the pile-up jet pT contributes to Emiss, and thus degrades
T both the Emiss response and resolution due to the stochastic
T 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 Emiss reconstruction leads to double-counting of the
T transverse momentum associated with the muon energy loss, as the fully reconstructed muon pT 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:
ͣ jet
• p T track Ϭ 0 8 – the transverse momentum of the
T trackϝ p
ID track associated with the muon (pT track) represents a jet
significant fraction of the transverse momentum p
T track, the sum of the transverse momenta of all ID tracks asso
ciated with the jet;
jet ͣ jet
• pT ϝ pT track Ϫ 2 – the overall transverse momentum p
T
of the jet is not too large compared to p
T track;
• NPV
track Ϫ 5 – the total number of tracks NPV
track 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 Emiss reconstruction. The
T muons are retained for the ETmiss reconstruction. Another consideration for muon contributions to ETmiss 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 Emiss
T reconstruction. Such jets are characterised by the following
selections, which are highly indicative of a photon in the ATLAS calorimeter:
• NPV
Ϫ 3 – the jet has low charged-particle content,
track indicated by a very small number of tracks from the hard
scatter vertex;
• fEMC Ϭ 0 9 – the jet energy Ejet is largely deposited in the electromagnetic calorimeter (EMC), as expected for photons and measured by the corresponding energy
Ejet
fraction fEMC = EMCEjet;
jet
• pT PS Ϭ 2 5 GeV – the transverse momentum contribujet jet
tion pT PS from presampler signals to pindicates an
T
early starting point for the shower;
• jet Ϫ 0 1 – the jet is narrow, with a width jet comparable to a dense electromagnetic shower; jet is reconstructed according to
i RipT i
jet =
ipT i
where Ri = i 2 + i 2 is the angular distance of topo-cluster i from the jet axis, and pT i is the trans-verse momentum of this cluster;
jet
• pϬ 0 8 – the transverse momentum
T trackϝ pT track
jet
pT track carried by all tracks associated with the jet is
close to p
T track.
Jets are accepted for Emiss reconstruction when consistent
T 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 Emiss
soft term
T
The soft term introduced in Sect. 3.2 is exclusively reconstructed from ID tracks from the hard-scatter vertex, thus only using the pT-flow from soft charged particles. It is an important contribution to Emiss for the improvement of both
T
the Emiss
T 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 ETmiss soft term including signals from soft neutral particles uses calorimeter topo-clusters. The reconstruction performance using the calorimeter-based Emiss soft calo is inferior to the track-only-
T
based Emiss soft
T , mostly due to a larger residual depen
dence on pile-up. More details of the topo-cluster-based Emiss soft calo
reconstruction are discussed in Sect. 8.1.
T
123
3.4.1 Trackand 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 NPV Ϭ 1 of reconstructed primary vertices.
Tracks are required to have pT Ϭ 400 MeV and || Ϫ 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 |d0| Ϫ 1 5mm (|z0 sinω | Ϫ 1 5 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 Tracksoft term
The track sample contributing to Emiss soft for each recon-
T structed event is collected from high-quality tracks emerging from the hard-scatter vertex but not associated with any electron, τ -lepton, muon, or jet contributing to Emiss recon-
T struction. The applied signal-overlap resolution removes
•
ID tracks with Rtrack,electron/photon cluster 0 05;
•
ID tracks with Rtrack τ -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 Emiss soft .
T Since only reconstructed tracks associated with the hard
scatter vertex are used, the track-based Emiss soft is largely
T 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 Emiss reconstruction performance
T uses selected final states without (Emiss true = 0) and
T with genuine missing transverse momentum from neutrinos (Emiss true
= pT). Samples with Emiss true = 0 are composed
TT of leptonic Zboson decays (Z → ee and Z → ) 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 → eυ and W → υ decays. In addition, samples with neutrinos and higher hard-object multiplicity were collected from top-quark pair (tt¯) production with at least either the t or the t¯ decaying semi-leptonically.
4.1 Data samples
The data sample used corresponds to a total integrated luminosity of 3 2fb−1, collected with a proton bunchcrossing interval of 25 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-oftime calorimeter deposits from cosmic-ray and beam backgrounds.
4.2 Monte Carlo samples
The Z → λ and W → υ 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 (partonlevel) matrix element output was interfaced to Pythia8 [21],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 tt¯-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 → λ and W → υ final state samples are the production of dibosons, single top quarks, and multijets. Dibosons were generated using Sherpa [28–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 (Wt) production, all interfaced to the PS and UE from the same Pythia6 configuration used for tt¯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
5 Version 8.186 was used for all final states generated with Pythia8.
123
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
Emiss
in T reconstruction, an alternative inclusive sample of Z → ͣ events was generated using the Mad-Graph_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 particles6 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 hardscatter 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 Emiss.
T The in-time pile-up activity is measured by the number of reconstructed primary collision vertices NPV. The out-oftime pile-up is proportional to the number of collisions per bunch crossing , 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 ͣ 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 ͣ distribution observed in data.
6 In ATLAS stable particles are those with an expected laboratory life-time τ corresponding to cϬ 10 mm.
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 → ͣ event selection
The Z → ͣ final state is ideal for the evaluation of ETmiss 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
Emiss
T are indicative of the reconstruction quality and reflect limitations introduced by both the detector and the ambiguity resolution procedure. The well-defined expectation value
Emiss true
T = 0 allows the reconstruction quality to be deter
mined in both data and MC simulations. The reconstructed Emiss
in this final state is also sensitive to the effectiveness
T 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 Emiss.
T Events must pass one of three high-level muon triggers
with different pT thresholds and isolation requirements. The isolation is determined by the ratio of the scalar sum of pT of reconstructed tracks other than the muon track itself, in a cone of size R = 0 2 around the muon track ( pTcone), to
ͣ
pT . The individual triggers require (1) p Ϭ 20 GeV and
T
ͣ
cone cone
pϝ p Ϫ 0 12, or (2) p Ϭ 24 GeV and pϝ p
TTT TT
0 06, or (3) pT Ϭ 50 GeV without isolation requirement.
Theofflineselectionof Z → ͣ 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) pT Ϭ 25 GeV, and (3) the reconstructed invariant mass mͣ of the dimuon system is consistent with the mass mZ of the Z boson, |mͣ − mZ| Ϫ 25 GeV.
5.2 W → eυ event selection
Events with W → eυ or W → υ in the final state pro-vide a well-defined topology with neutrinos produced in the hard-scatter interaction. In combination with Z → , the effectiveness of signal ambiguity resolution and lepton energy reconstruction for both the electrons and muons can be observed. The W → eυ events in particular provide a good
metric with Emiss true = pϬ 0 to evaluate and validate
TT the scale, resolution and direction (azimuth) of the recon
structed Emiss,asthe Emiss reconstruction is sensitive to the
TT
123
electron–jet overlap resolution performance. This metric is
only available in MC simulations where pT is known. Candidate W → eυ events are required to pass the high-level electron trigger with pT Ϭ 17 GeV. Electron candidates are selected according to criteria described in Sect. 3.3.1.Only events containing exactly one electron are considered.
Further selections using Emiss and the reconstructed trans-
T verse mass mT, given by
J
T eEmiss
mT = 2 pT 1 − cos ψ
are applied to reduce the multijet background with one jet
emulating an isolated electron from the W boson. Here Emiss
T is calculated as presented in Sect. 3. The transverse
e
momentum of the electron is denoted by pT, and ψ is the distance between miss and the azimuth of the electron. Selected events are required to have Emiss Ϭ 25 GeV and
T mT Ϭ 50 GeV.
5.3 tt¯event selection
Events with tt¯in the final state allow the evaluation of the Emiss
T 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, respec
tively, and are required to have pT Ϭ 25 GeV.
The final tt¯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-kt 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 R Ϭ 0 4 from the electron.
6 Performance of ETmissreconstructionin data and MonteCarlo 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, Emiss reconstruc-
T tion yields a non-linear response, especially in regions of phase space where the observation bias discussed in Sect. 3.1 dominates the reconstructed ETmiss. In addition, the ETmiss resolution functions are characterised by a high level of
complexity, due to the composite character of the observable. Objects with different pT-resolutions contribute, and the Emiss
T composition can fluctuate significantly for events from the same final state. Due to the dependence of the Emiss
T response on the resolution, both performance characteris
tics 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 ETmiss.
Emiss
The T reconstruction performance is therefore assessed by comparing a set of reconstructed Emiss-related
T observables in data and MC simulations for the same finalstate selection, with the same object and event selections applied. Systematic uncertainties in the Emiss response and
T 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 Emiss
T 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 Emiss perfor-
T mance evaluation are obtained by applying selections accord
ing to Sect. 5.1 for a final state without genuine Emiss
T (Z → ), and according to Sect. 5.2 for a final state with
genuine Emiss (W → e). From these, specific exclusive
T samples are extracted by applying conditions on the number of jets reconstructed. In particular, zero jet (Njet = 0) samples without any jet with pT Ϭ 20 GeV (fully calibrated) and || Ϫ 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 pT Ϭ 20 GeV are useful for evaluating the contribution of jets to Emiss.
T While the pT response of jets is fully calibrated and provides a better measurement of the overall event pT-flow, the pT resolution for jets is affected by pile-up and can introduce a detrimental effect on Emiss reconstruction performance.
T 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 Emiss recon-
T struction sequence discussed in Sect. 3.1. It employs the hard-object selections defined in Sects. 3.3.1–3.3.4, with jets selected according to the prescriptions given in Sect. 3.3.5.
123
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 Emiss
T modelling in Monte Carlo simulations
, Emiss, Emiss
The quality of the MC modelling of Emiss and
xy T L ET, reconstructed as given in Eqs. (1), (3) and (5), is eval
uated for an inclusive sample of Z → ͣ events by comparing the distributions of these observables to data. The results are presented in Fig. 1. The data and MC simula
tions agree within 20% for the bulk of the Emiss distribu-
T tion shown in Fig. 1a, with larger differences not accommo
dated by the total (systematic and statistical) uncertainties of the distributions for high Emiss. These differences sug-
T gest a mismodelling in tt¯ events, the dominant background in the tail regime [43]. The L ET distributions compared between data and MC simulations in Fig. 1b show discrep
ancies significantly larger than the overall uncertainties for 200 GeV L ET Ϫ 1 2 TeV. These reflect the level of mismodelling of the final state mostly in terms of hard-object and Emiss
composition in MC simulations. The Emiss spec
xy
tra shown in Fig. 1c, d, respectively, show good agreement between data and MC simulations for the bulk of the distributions within |Emiss| Ϫ 100 GeV, with larger differences
xy
observed outside of this range still mostly within the uncertainties. The distributions of individual contributions to ETmiss from jets (Emiss jet
), muons (Emiss ), and the soft term (Emiss soft ),
TT T as defined in Eq. (6), are compared between data and MC sim
ulations for the same inclusive Z → ͣ sample in Fig. 2. Agreement between data and MC simulations for Emiss jet in
T Fig. 2a is of the order of ±20% and within the total uncer
tainties for Emiss jet � 120 GeV, but beyond those for higher
T
Emiss jet
. A similar observation holds for Emiss ͣ in Fig. 2b,
TT where data and MC simulations agree within the uncertain
ties for low Emiss ͣ but significantly beyond them for larger
T
Emiss ͣ
T . Agreement between data and MC simulations is bet-ter for the soft term Emiss soft , with differences up to 10% for
T
Emiss soft
T � 30 GeV, as seen in Fig. 2c. Larger differences for larger Emiss soft
are still found to be within the uncertainties.
T The peak around Emiss jet = 20 GeV indicates the onset
T of single-jet events at the threshold pT = 20 GeV for jets
contributing to Emiss jet . Larger values of Emiss jet arise from
TT events with one or more high-pT jets balancing the pT of the
Z boson.
For the W → eυ sample with genuine missing transverse
momentum given by pT, both the total reconstructed Emiss
T and the soft term are compared between data and MC simu
lations in Fig. 3. The level of agreement between the Emiss
T distributions for data and MC simulations shown in Fig. 3a
for the inclusive event sample is at ±20%, similar to that observed for the Z → ͣ 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 ETmiss distributions observed with the exclusive Njet = 0 sample shown in Fig. 3b are well below 20%, but show a trend to larger discrepancies for decreasing Emiss � 40 GeV. This
T 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 Emiss required in the W → eυ selection described
T in Sect. 5.2. The comparison of the Emiss soft distributions
T from data and MC simulations shown in Fig. 3c yields agree
ment well within the uncertainties, for Emiss soft 2 10 GeV.
T The rising deficiencies observed in the MC distribution for
decreasing Emiss soft � 10 GeV are expected to be related to
T the missing multijet contribution.
6.2 Emiss
T response and resolution
The response in the context of Emiss reconstruction is deter-
T mined by the deviation of the observed Emiss from the expec-
T tation value for a given final state. This deviation sets the scale
for the observed Emiss. If this deviation is independent of the
T genuine missing transverse momentum, or any other hard
pT indicative of the overall hard-scatter activity, the Emiss
T response is linear. In this case, a constant bias in the recon
structed Emiss is still possible due to detector inefficiencies
T and coverage (acceptance) limitations. Final states balanced in transverse momentum are expected to show a non-linear Emiss response at low event activity, as
T the response in this case suffers from the observation bias in Emiss
reconstruction discussed in Sect. 3.1. With increas-
T ing momentum transfers in the hard-scatter interaction, the Emiss
T response becomes increasingly dominated by a wellmeasured hadronic recoil and thus more linear. In the case of final states with genuine missing transverse momentum,
the Emiss
response is only linear once Emiss true exceeds the
TT observation bias. These features are discussed in Sect. 6.2.1 and explored in Sect. 6.2.2. Contributions to the fluctuations in the ETmiss 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
123
L ET, c Emiss and d Emiss
Fig. 1 Distributions of a ETmiss, b xy for an inclusive sample of Z → ͣ 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
coverage of || Ϫ 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 Emiss to
T
the reconstructed and accepted e, γ , had and , and those
pT 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 cand 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
being part of a reconstructed and accepted jet. In addition, the pT-flow of not explicitly reconstructed charged particles emerging from the hard-scatter vertex is represented by ID tracks contributing to Emiss soft given in Eqs. (6) and (7), but
T
123
Fig. 2 Distributions of athe jet term Emiss jet , bthe muon term Emiss ͣ ,
TT and c the soft term Emiss soft for the inclusive samples of Z →
T
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
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 ETmiss reconstruction.
respective systematic uncertainties from athe jet, bthe muon, and cthe 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
Like for the Emiss response, resolution-related aspects
T
of Emiss
reconstruction are understood from data-to-MC
T
simulations comparisons. The scales used for the correspond
123
Fig. 3 Distributions of the total Emiss in a theinclusive caseand b
T the Njet = 0 case, as well as c the soft term Emiss soft reconstructed
T
in Njet = 0eventswith W → eυ in data. The expectation from MC simulation is superimposed and includes all relevant background final states passing the event selection. The inclusive Emiss distribution from
T
MC simulations contains a small contribution from multijet final states at low Emiss, which is absent for the Njet = 0 selection. The shaded
T
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
123
ing evaluations are the overall event activity represented by L ET, and the pile-up activity measured by NPV. The measurement of the ETmiss resolution is discussed in Sect. 6.2.3 and results are presented in Sect. 6.2.4.
Emiss
6.2.1 scale determination
T
In events with Z →ͣ decays, the transverse momentum of the Z boson ( pT Z) is an indicator of the hardness of the interaction. It provides a useful scale for the evaluation of the
Emiss
T response for this final state without genuine missing
transverse momentum. The direction of the corresponding Z Z
boson transverse momentum vector pT defines an axis AZ in the transverse plane of the collision, which is reconstructed from the pT of the decay products by
+ −
p +ppZ
AZ = TT = T (9)
+ − Z
p
p +p T
TT
The magnitude of the component of Emiss parallel to AZ is
T
PZ Emiss
=·AZ (10)
f T
This projection is sensitive to any limitation in Emiss recon-
T
struction, in particular with respect to the contribution from
Z
the hadronic recoil against pT , 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 Emiss response and the associated systematic
T
uncertainties.
The expectation value for a balanced interaction producing a Z boson against a hadronic recoil is E[PfZ]=0. Any observed deviation from this value represents a bias in the Emiss reconstruction. For PfZϪ 0, the reconstructed
T hadronic activity recoiling against pT Z is too small, while for PfZϬ 0 too much hadronic recoil is reconstructed. The evo
lution of PfZas a function of the hardness of the Z boson production can be measured by evaluating the mean (PfZ)in
hard Z
bins of the hard-scatter scale pT =pT. In addition to measuring the Emiss response in data and
T
MC simulation without genuine Emiss, its linearity can be
T
determined using samples of final states with genuine Emiss
T
in MC simulations. This is done by evaluating the relative deviation lin of the reconstructed Emiss from the expected
TT
Emiss true
Ϭ 0 as a function of Emiss true ,
TT
Emiss −Emiss true
lin TT
T Emiss true
= (11)
T Emiss true
T
6.2.2 MeasuringtheE miss response
T
Z
Figure 4 shows (PfZ)as a function of pfor the Njet =0
T and the inclusive Z →ͣ sample, respectively. MC simulations compare well with the data for Njet =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.
Z
The steep decrease of (PfZ)with increasing pin the
T Njet =0 sample seen in Fig. 4a reflects the inherent under
estimation of the soft term, as in this case the hadronic recoil is exclusively represented by ID tracks with pT Ϭ 400 MeV within ||Ϫ 2 5. It thus does not contain any signal from (1) neutral particles, (2) charged particles produced with ||Ϭ 2 5, and (3) charged particles produced within ||Ϫ 2 5but with pT 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 Emiss Z
response is recovered better as pT increases, since an
T increasing number of events enter the sample with a reconstructed recoil containing fully calibrated jets. These pro-vide a more complete representation of the hadronic trans-verse momentum flow. The residual offsets in (PfZ)of about 8 GeV in data and 6 GeV in MC simulations observed for
Z
pT 2 40 GeV in Fig. 4b agree within the uncertainties of this measurement.
The persistent bias in (PfZ)is further explored in Fig. 5, which compares variations of (PfZ)respectively using the full Emiss
, the soft-term contribution Emiss soft only, the hard-
TT
−Emiss soft
term contribution EmissT T , and the true soft term Emiss true soft Z
only, as a function of pT ,for the Z →
T sample from MC simulations. In particular the difference and Emiss soft
between the projections using Emiss true soft
TT indicates the lack of reconstructed hadronic response, when Emiss soft Emiss true soft
T =T is expected for a fully measured recoil. The parallel projection using only the soft terms is larger than zero for all pT Z due to the missing Z-boson con-
Z
tribution to Emiss given by −p
T T. The deviation from linearity in Emiss reconstruction, mea-
T sured by lin given in Eq. (11), is shown as a function
T
of Emiss true
T for MC simulations of W →e, W →υ and tt¯production in Fig. 6. The observed lin Ϭ 0at
T
low Emiss true
indicates an overestimation of Emiss true by
TT the reconstructed ETmiss due to the observation biases arising from the finite Emiss resolution, as discussed in Sect. 3.1.
T This bias overcompensates the lack of reconstructed pT-flow from the incompletely measured hadronic recoil in W →eυ and W →υ events for ETmiss true 40 GeV with an increasing non-linearity observed with decreasing Emiss true .
T
For Emiss true
2 70 GeV the Emiss response is directly pro-
TT
123
Fig. 4 The average projection of Emiss onto the direction AZ of the Z
T
boson’s transverse momentum vector pT Z,asgiven in Eq.(10), is shown
ZZ
as a function of p=|pT |in Z →ͣ events from athe Njet =0sam
T
ple and from bthe inclusive sample. In both cases data are compared to MC simulations. The ratio of the averages from data and MC simula-
ZZ
as a function of p=|pT |in Z →ͣ events from the inclusive MC
T 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
portional to Emiss true , with the reconstructed recoil being
T approximately 2% too small. The W →eυ and W →υ T Emiss true
final states show very similar lin , thus indicating
T the universality of the recoil reconstruction and the indepen
dence on the lepton flavour of the reconstructed Emiss in a
T low-multiplicity final state with Emiss true Ϭ 0.
T
tions 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
Fig. 6 The deviation of the Emiss response from linearity, measured
T by lin
as a function of the expected Emiss true in Eq. (11), in W →e,
TT W →,and tt¯final states in MC simulations. The lower plot shows
a zoomed-in view on the lin dependence on Emiss true with a highly
TT suppressed ordinate
In tt¯final-state reconstruction, resolution effects tend at Emiss true
to dominate lin 120 GeV. Compared to
TT the W → eυ and W → υ final states, a signifi
123
cantly poorer Emiss resolution is observed in this kinematic
T region, due to the presence of at least four jets with relatively low pT and high sensitivity to pile-up-induced fluctuations in each event of the tt¯sample. For Emiss true
T
T Emiss true
120 GeV, lin ≈ 2% indicates a proportional
T
Emiss
T response with a systematic shift similar to the one observed in inclusive W-boson production.
6.2.3 Determination oftheE Tmiss resolution
The Emiss
T resolution is determined by the width of the com
bined distribution of the differences between the measured Emiss
xy and the components of the true missing transverse Emiss true Emiss true
momentum vector Emiss true = xy .The
T width is measured in terms of the RMS, with
⎧
Emiss − Emiss true t sample Emiss true
⎨ RMS W → eυ or t¯Ϭ 0
xyxy T
RMSmiss
=
xy ��
⎩ Emiss true Z → ͣ sample Emiss true
RMS =0
xy T
(12)
This metric does not capture all of the effects driving the fluctuations in ETmiss reconstruction, such as biases between individual Emiss terms or the behaviour of outliers, but it is
T an appropriate general measure of how well Emiss represents
T
Emiss true
T.
Using the Z → ͣ sample allows direct comparisons of RMSmiss
between data and MC simulations, as Emiss true = 0
xy T in this case. The resolution in final states with genuine Emiss
T is determined with MC simulations alone. For W → eυ and
tt¯final states, Emiss true = pis used.
xyxy
Emiss
6.2.4 resolution measurements
T
The Emiss
resolution measured by RMSmiss is evaluated as
T xy a function of the event activity measured by L ET given in
Eq. (7). For the inclusive Z → ͣ sample, Fig. 7ashows RMSmiss
xy quickly rising from less than 5 GeV to about 10 GeV with increasing L ET within 50 GeV ↕ L ET Ϫ 70 GeV.7 This is due to the fact that in this range the two muons are the dominant hard objects contributing, with a pT resolution
proportional to ( p 2. A convolution of the muon resolu-
T tion with a small contribution from Emiss soft is possible for
T L ET Ϭ 50 GeV. This component is on average about 60% Z
of pT , and subject to the stochastic fluctuations further discussed below. The increase of Z → +1 jet topologies in the Z → ͣ sample leads to an additional source of fluctuations affecting
7 This lower boundary of this range is given by the muon selection with
ͣ
pT Ϭ 25 GeV, as described in Sect. 5.1, assuming no other hard-scatter vertex tracks, i.e. Emiss soft = 0. The upper boundary indicates the lower
T jet
limit of L ET to accommodate at least one jet with pT Ϭ 20 GeV in addition the two muons (for the jet selection see Sect. 3.3.5).
RMSmiss
xy L ET for 70 GeV L ET 180 GeV. In gen
eral the Z → ͣ sample collected for this study covers Z
pT 140 GeV with relevant statistics. At this limit it is expected that the hadronic recoil contains two reconstructed jets, with the onset of this contribution at L ET of about 180 GeV. The corresponding change of the dominant final state composition for L ET Ϭ 180 GeV leads to a change of shape of RMSmissL ET , as the transverse momentum of
xy
the individual jets rises and the number of contributing jets slowly increases. The expected RMSmissL ET ∝√ L ET
xy scaling driven by the jet-pT resolution [44] therefore dom
inates RMSmiss at these higher L ET. The MC predictions
xy for RMSmissL ET agree with the data within a few percent
xy
and well within the total uncertainties of this measurement. A tendency for slightly poorer resolution in MC simulations is observed, in particular for L ET Ϭ 200 GeV.
Any contribution from pile-up to RMSmiss is expected
xy
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 jetpT measurement. The dependence of RMSmiss on the pile-up
xy activity measured by NPV is shown in Fig. 7b. Data show a
less steep slope of RMSmiss than MC simulations, but
xy NPV with about 10% worse resolution in the low pile-up region of NPV 5. The resolution in data is better than in MC simulations by about 10% for the region of higher pile-up activity at NPV ≈ 20. The differences between data and MC simulations seen in RMSmissL ET for the inclusive Z → ͣ sample can
xy
be further analysed by splitting the sample according to the value of Njet. Figure 8a shows the dependence of RMSmiss
xy on L ET for Z → ͣ events with Njet = 0. The dominant source of fluctuations other than the muon-pT resolution is in this case introduced by the incomplete reconstruction of the hadronic recoil. These fluctuations increase with increasing pT Z, which in turn means higher overall event activity measured by L ET. For this sample RMSmiss in data compares
xy
well to MC simulations, at a level of a few percent, without any observed dependence on L ET.
The exclusive Njet ↖ 1 samples extracted from Z → ͣ data and MC simulations show the expected RMSmiss ∝
xy
√
L ET scaling in Fig. 8b. The resolution in data is well rep
resented by MC simulations, at the level of a few percent. The slightly better resolution observed in data with increasing L ET 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 RMSmiss on NPV shown in Fig. 8c
xy indicates that the Emiss resolution is basically independent
T of pile-up, for the Njet = 0 sample. This is expected from the
123
Fig. 7 The RMS width of the Emiss distributions a in bins of L ET
xy
and bin bins of the number of primary vertices in an inclusive sample of Z → ͣ events. Predictions from MC simulations are overlaid on
exclusive Emiss composition comprising the (track-based)
T
Emiss ͣ and Emiss soft
T T terms only. Data and MC simulations compare well within a few percent, and without any observable dependence on NPV. Figure 8dshows the NPV dependence of RMSmiss for the Njet ↖ 1 sample. Comparing this
xy
result to Fig. 7b confirms that all pile-up dependence of the Emiss
T resolution is arising from the jet term. Both trend and magnitude of the data-to-MC comparison follow the observation from the inclusive analysis.
Emiss
6.2.5 resolutionin final states with neutrinos
T
The Emiss
T resolution for final states with ETmiss true Ϭ 0is measured by RMSmiss according to Eq. (12) and evaluated
xy using dedicated inclusive W → eυ and W → υ samples from MC simulations, and the inclusive tt¯MC sample defined in Sect. 5.3. For these samples, RMSmiss can be
xy
determined as a function of Emiss true = pT. The dedicated
T W → eυ and W → υ samples are obtained with an event selection based on the description in Sect. 5.2, but omitting both the Emiss-based and the mT-based selections.
T Figure 9 shows RMSmiss evaluated as a function of
xy Emiss true
T for these samples. The universality of the response to the hadronic recoil observed in Fig. 6, together with the dif
ferent but subdominant contributions from the pT resolutions of the electrons and muons, yield a very similar Emiss resolu-
T tion for W → eυ and W → υ final states. Generally, poorer
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
resolution is observed in tt¯ final states. The deviation from
J
Emiss true
the expected RMSmiss ∝ scaling behaviour for
xy T W → υ at lower Emiss true
reflects the kinematic features of
T the W boson and its decay. Events with low pT W, and therefore small hadronic recoil, lie predominantly in the region
25 GeV pT 50 GeV. Since the hadronic recoil is generally the poorly measured component of an event and the reconstructed Emiss is dominated by the lepton pT in this
T region, the Emiss resolution tends to be better here than for
T
events with larger hadronic recoil populating pT 25 GeV
and pT 2 50 GeV.
6.3 Emiss
tails
T
Large reconstructed Emiss is an indicator for the production
T 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 Emiss components for final states with
T well-known expectation values for Emiss are indicative of
T such inefficiencies. Non-Gaussian shapes in the distribution arise from a combination of object selection inefficiencies and potentially non-Gaussian resolutions of the ETmiss constituents. Even for a well-defined final state, event-by-event fluctuations in terms of which particles, jets, and soft tracks enter the Emiss recon-
T struction, and with which pT, lead to deviations from a nor,Emiss
mally distributed (Emiss ) response.
xy
123
The Emiss
Fig. 8 resolution RMSmiss determined for a an exclusive
T xy
Z → ͣ sample without jets with pT Ϭ 20GeV (Njet = 0) and for b an exclusive sample with at least one jet above this threshold (Njet ↖ 1), as a function of L ET in data and MC simulations. The dependence of
,Emiss
Figure 10 shows the combined (Emiss ) distribution
xy
for the inclusive Z → ͣ 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, ±1 × RMS and ±2 × RMS.
RMSmiss
xy on the pile-up activity, as measured by NPV,for thesetwo samples is shown in c and d, respectively. The shaded bands indicate the combined statistical and systematic uncertainties associated with the measurement
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 Emiss with near exponen
xy
tial slopes. The result of this comparison supports the choice of RMSmiss defined in Eq. (12) in Sect. 6.2.3 for the deter
xy mination of the Emiss resolution, rather than using any of
T
123
and Emiss
xy
Z → ͣ from simulation. Gaussian fits limited to the ±1 × RMS and ±2 × RMS ranges around the centre of the distribution are shown, together with the respective differences between the fitted functions and the actual distribution
the widths measured by fitting Gauss functions in selected ranges of the distribution.
The tails in this shape are reflected in the distribution of Emiss
T itself and can be estimated by measuring the fraction Ϭ Emiss threshold
of events with ETmiss T ,
1 ∞
1 hEmiss dEmiss
ftail = TT
HEmiss threshold
T
1
∞
hEmiss dEmiss
with H = (13)
TT
0
Here hEmiss is the Emiss distribution for a given event sam-
TT
ple, and Emiss threshold
T is a threshold set to estimate tails. Any decrease of ftail at a fixed integral H indicates an improvement of the Emiss resolution, and is more sensitive to particu-
T lar improvements than e.g. RMSmiss. For example, improving
xy the Emiss soft
T reconstruction by rejecting ID tracks from the hard-scatter vertex with poor reconstruction quality yields a significantly smaller ftail for the same event sample.
The tails in the ETmiss distributions for the final states considered for this study are quantified by the fraction of events above a certain Emiss threshold using MC simulations. Fig-
T ure 11a shows that the Z → λ events (λ = e or λ = ) with Emiss true
T = 0 have significantly reduced tails when compared to W → υ and tt¯ with this metric, and that the tails do not depend on the lepton flavour. A modification of this metric, taking into account Emiss true such that the fraction
T
− Emiss true
of events with |Emiss | above a given threshold is
TT determined, shows the universality of the hadronic recoil in Z → λ and W → , as can be seen in Fig. 11b. Another finding of this study is that the tail in the |Emiss −
T ETmiss true | distribution for the higher L ETjet tt¯sample is con
siderably larger than for the low-L Ejet samples with Z →
T or W → υ final states. As can be seen in Fig. 11c, the tails are much more consistent between Z → ͣ and tt¯samples when the distribution for the Z → ͣ sample is reweighted such that it follows the same L Ejet distribution as the tt¯
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 Emiss are provided for the response (Emiss scale) as
TT 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 Emiss terms need to be extracted
T from the scale and resolution uncertainties for the individual contributing objects comprising electrons, photons, muons, τ -leptons, and jets. In the corresponding propagations, cor
123
T tions of Z → , W → ,and tt¯final states. The tail fraction in terms
− Emiss true
of a threshold applied to |Emiss |, the distance between the
TT reconstructed (Emiss) and the expected (Emiss true ) vectors, is shown in
TT
relations 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 Emiss reconstruction are uncorrelated. The
T determination of the ETmiss scale and resolution uncertain
ties arising from the soft term Emiss soft is described in this
T section.
bfor all considered final states. The same fraction is shown in c for the Emiss
T distributions for Z → ͣ before and after a reweighting following the L ET distribution for tt¯is applied, together with ftail from the tt¯final state
7.1 Methodology
The extraction of the systematic uncertainties for the reconstructed Emiss is based on data-to-MC comparisons of spectra
T of observables measuring the contribution of Emiss soft to the
T overall ETmiss.
123
7.1.1 Observables
The vector sum of the transverse momentum vectors of all particles and jets emerging from a hard-scatter interaction (pHS)isgiven by
T
HS e ͣ jet
p= ppppp
TT + T + T + T + T
pT (observable)
+ pυ T
inv
obs
p(not observable)
T
Here pυ T generally represents the transverse momenta of noninv
observable particles, which are summed up to form pT .All
other transverse momenta are carried by particles that are obs
observable in principle, and sum up to pT . Momentum con-HS HS
servation dictates p=|p|=0.
TT Due to detector acceptance limitations and inefficien
cies in hard-object reconstruction and calibration, and all hard
other effects discussed in Sect. 3, only a proxy (pT )for the observable-particle contribution pobs can be measured.
T The reconstructed hard final-state objects entering Emiss as
T described in Sect. 3.2 are used to measure phard as
T
hard e had
p= pT + pT + p
TT contributing contributing contributing electrons photons τ -leptons
ͣ jet
+ pT + p
T contributing contributing
muons jets
hard hard hard obs
The expectation is that p=|p|Ϭ 0 and p =p.
TT TT
=−Emiss soft , with Emiss soft
Adding psoft defined in Eq. (6),
TT T
hard
to pT yields an improved estimate of the net transverse momentum carried by undetectable particles, as some of the experimental inefficiencies are mitigated.8
In the Z →ͣ final state without genuine missing trans-hard
verse momentum the expectation is that EmissT =−p+
T
soft
pT =0. While this expectation does not hold due to the
experimental inefficiencies, it nevertheless raises the expecsoft
tation that for events without jets pT points into the direction of the hadronic recoil, i.e. opposite to phard in the transverse-
T momentum plane. The deviation from this expectation is
measured in terms of the parallel (Pf) and perpendicular soft hard
(P⊥) projections of pT onto pT . Figure 12 schematically shows these projections for Z+0-jet and Z+1-jet topologies. The average (Pf)in a given bin kof phase space defined by
hard hard
pT measures the Emiss soft response, with (Pf)=(pT )k
T indicating a perfect response in this bin. The Emiss resolution
T
8 As discussed in Sect. 3.4, the soft term represents only charged parti
cles with pT Ϭ 400 MeV not associated with fully identified and reconsoft
structed particles or jets. Therefore, including pT can only recover a part of the actual soft pT-flow of the interaction.
contribution from Emiss soft reconstruction is measured by
T two components, the fluctuations in response (RMSf2) and
the fluctuations of the (transverse) angular deflection around hard
the paxis, measured by RMS⊥2 . These fluctuations are
T expressed in terms of variances, with
��2
22
RMSf 2 =Pf− Pf and RMS⊥ 2 =P⊥
7.1.2 Procedures
The extraction of the systematic uncertainties introduced into the Emiss
measurement by the Emiss soft term is based on
TT hard
data-to-MC-simulations comparisons of (Pf)( pfor the
T
hard hard
response, and of RMSf2( pand RMS⊥2 ( pfor the res-
TT olution. 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 Emiss soft on Emiss
T T , the exclusive Z → ͣ selection with Njet = 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
hard Z
the reconstructed Z boson, i.e. pT = pT as shown in Fig. 12a.
The uncertainties are determined by comparing Pf and P⊥ spectra from data and MC simulations, in bins of phard.
T For Pf, the smearing of the response and the width both yield scale and (longitudinal) resolution offsets. In the case of P⊥, only smearing of the width is applied to provide trans-verse 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, Gaus-sian smearing is applied in simulation to the longitudinal and transverse components of Emiss soft relative to the direction
T
hard
of pT . To account for differences in response between data and simulation, the longitudinal component of Emiss soft is
T 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.
123
(b)
soft hard
projections of pon pfor Z →ͣ events without genuine Emiss,
TT T
for a a final state without any jets and ba final state with one jet. The
7.2 Systematic uncertainties in Emiss response and
T
resolution
Emiss
The result for the systematic uncertainty of the T scale, determined as discussed in the previous section, is summarised in Fig. 13. The average longitudinal projec
soft hard hard
tion of ponto p, (Pf), as a function of pis
TT T
shown in Fig. 13a which compares data to both the stan
dard 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 (Pf)MC within the uncertainties of the data. The lower panel of Fig. 13a confirms that the ratio (Pf)MC(Pf)data lies within the systematic uncertainty band over the full phard
T
range. The systematic uncertainty for the Emiss resolution is
T extracted from the variances of the parallel (RMSf2) and per-
hard
pendicular (RMS⊥2 ) projections of Emiss onto pdefined in
TT
hard
Sect. 7.1.2. Figure 13bshows the pdependence of RMS2
T f
measured for the exclusive Z → ͣ sample (Njet = 0) in data and two MC simulations. The variances RMSf 2 MC calculated for both sets of simulations agree within the systematic uncertainties of RMSf 2 data with the data, as illustrated in the lower panel of the figure, where the ratio
hard
RMSf 2 MCRMSf 2 data is shown as a function of p.The
T
Z
expectation values for a perfect Emiss reconstruction are E[Pf]=p
TT
hard
for Njet =0and E[Pf]=pfor Njet ↖1, with E[P⊥]=0in all
T cases. a Z +0 jet topology. b Z +1 jet topology
results of the evaluation of the variances RMS⊥2 of the perpendicular projections as a function of pThard are shown in Fig. 13c, together with the resulting pThard dependence of the ratio RMS⊥ 2 MCRMS⊥ 2 data. The systematic uncertainties of the data cover all differences to MC simulations.
8Missing transversemomentumreconstruction variants
8.1 Calorimeter-based ETmiss
The Emiss
T soft term from the calorimeter ETmiss soft 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 Emiss soft calo
T is calculated from topo-clusters calibrated at the LCW scale. Only topo-clusters with a calibrated energy ELCW Ϭ 0, not contributing to the recon
clus struction of the hard objects used to calculate the hard term given in Eq. (6), are considered for Emiss soft calo . In addi-
T tion, topo-clusters that are formed at the same location as the
hard object signals are not considered for Emiss soft calo even
T if their signals are not directly contributing to the reconstruc
tion of the hard objects. The fully reconstructed Emiss using
T Emiss soft calo is Emiss calo
T T.
123
24
140
22 20
120
18 16
100
MC/Data
[GeV]
10 8
40 2
60
6
4
RMS2 [GeV2]
14 12
80
0.8
0.8 102030405060708090100 102030405060708090100
phard [GeV] phard [GeV]
(a) T(b) T
MC/Data
0
20
1.2 1.2
1.1 1.1
1
1
0.9 0.9
102030 40
the b variance RMS2 of the longitudinal projection Pf of ponto
f T phard for Z → ͣ event with Njet = 0, for data and two different
T hard
MC simulations, shown as a function of pT . The variance RMS⊥ 2 of
and L ET, Emiss calo
Compared to the reference Emiss and
TT
L Ecalo
T have an enhanced dependence on pile-up, mostly introduced by the soft term. To partly compensate for the irreducible contribution of pT-flow reconstructed from topoclusters generated by pile-up to Emiss calo , a modified jet
T 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 pT Ϭ 20 GeV. The contribution of these jets to Emiss calo and L Ecalo
T T , defined in terms of momentum compo
50 60 70 80 90100
phard [GeV]
(c) T
the perpendicular projection P⊥ is shown in c for the same event samples. The shaded band indicates the systematic uncertainties derived as described in the text
nents ( px py , depends on the overlap with already accepted reconstructed particles,
00 E ↖ 50% (large overlap)
( px py = jet jet
1 − E × ( px py E Ϫ 50% (small or no overlap)
(14)
The overlap fraction E is given in Eq. (8). Jets with E ↖ 50% are not used at all. The JVT-based tagging of non-pileup jets is omitted. It is found that this strategy reduces the fluctuations in the ETmiss calo reconstruction. The transverse
123
momentum contribution of groups of clusters representing a jet-like pT-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-pT threshold applied in Emiss reconstruction, is reduced
T if all jets and jet fragments, including those from pile-up, are included.
8.2 Emiss
from tracks
T
The reference track-based soft term Emiss soft is largely insen-
T sitive to pile-up, as indicated by the dependence of the Emiss
T resolution RMSmiss on NPV in the exclusive Z → ͣ sam
xy ple (Njet = 0) shown in Fig. 7c. As discussed in Sect. 6.2.4 and from the comparison of Fig. 7c, d, the pile-up depen
dence of RMSmiss in the inclusive Z → ͣ sample is largely
xy
introduced by the jet contribution. This contribution suffers from (1) the lack of pile-up suppression for forward jets with || Ϭ 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 Emiss
T 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 ( pmiss) all jets and reconstructed particles are ignored,
T
i.e. the pTmiss reconstruction does not include any calorimeter or MS signals. The pmiss resolution is then inherently
T immune to pile-up, while the pTmiss response is low as all neutral pT-flow in || Ϫ 2 5aswellasall pT-flow outside of this region is excluded.
8.3 Performance evaluations for ETmiss variants
The main motivation to study Emiss-reconstruction variants
T is to improve some combination of the Emiss resolution,
T scale, and stability against pile-up. As with the composition
of objects entering Emiss reconstruction in general, the par-
T ticular 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 Emiss calo miss
T and pT to the corresponding measurements using the reference Emiss illustrates their principal features for the
T Z → ͣ and tt¯production final state.
8.3.1 Comparisons ofE miss resolution
T miss
Figure 14 compares the Emiss calo and presolutions with
TT the one obtained from the reference Emiss, for the inclusive
T Z → ͣ sample in data. Each is shown as a function of L ET corresponding to the reference Emiss, giving an esti-
T mate of the total hard-scatter activity. The low-L ET region
is dominated by events with Njet = 0, where the contribution of Emiss soft calo in Emiss calo
T T yields a poorer resolution than
for Emiss miss
, and where Emiss and phave identical perfor-
T TT mance. The high-L ET region is dominated by events with
higher jet multiplicity, where pmiss resolution is degraded rel-
T ative to the reference Emiss by the incomplete measurement
T of jets. miss
Figure 15a compares the Emiss calo and presolution as
TT functions of the pile-up activity measured by NPV, with the one obtained from the reference ETmiss for the exclusive Z → ͣ samples with Njet = 0 in data. The Emiss calo resolution
T is dominated by pile-up and shows significantly degraded miss
performance relative to pand the reference Emiss.The
TT exclusive use of only tracks from the hard-scatter vertex for miss and Emiss
both pT T yields the same stability against pile-up.
In events with jet activity, the degraded pTmiss resolution is observable, especially outside the region of high-est pile-up activity, as seen in Fig. 15bfor the Emiss res-
T olution obtained with the inclusive Z → ͣ sample in data for NPV 15. This is even more obvious in final states with relatively high jet multiplicity and genuine missing transverse momentum, like for the tt¯-production sam
ple from MC simulations. As shown in Fig. 15c for this final state, both the reference ETmiss and the calorimeterbased Emiss calo
T have a significantly better resolution than pmiss, at the price of some sensitivity to pile-up, which is
T absent for pmiss.The NPV dependence of the resolution is
T enhanced in Emiss calo , due to the increased contribution from
T soft calorimeter signals without pile-up suppression at higher
NPV.
123
T
olutions of the track-only-based variant pmiss described in Sect. 8.2,
T and the reconstruction variant Emiss calo employing a calorimeter-based
T
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
8.3.2 Comparisons ofE miss scale
T
Following the description in Sect. 6.2.1,the Emiss response is
T
, Emiss calo miss
evaluated for the reference Emiss , and pusing
TT T
, Emiss calo miss
the respective projections of Emiss , and ponto
TT T
the direction of pT Z, according to Eqs. (9) and (10). Fig-
Z
ure 16a shows the average projection as a function of pT for the exclusive Z → ͣ sample with Njet = 0 in data.
miss
Both Emiss and pshow the same increasingly incomplete
TT
reconstruction of the hadronic recoil in this sample for rising
measured in terms of the number of reconstructed vertices NPV for a an exclusive Z → ͣ sample without jets with pT Ϭ 20 GeV and ban inclusive Z → ͣ sample, both selected from data. In c, the resolution of the Emiss reconstruction-variants in a final state with significant jet
T
activity and pT Ϭ 0 is compared using MC simulations of tt¯production
Z
pT . This reconstruction is slightly improved for Emiss calo ,
T
but still insufficient at higher pT Z . In the inclusive Z → ͣ sample, shown in Fig. 16b, the
T is that Emiss calo
indication at lower pZ has a higher response
T
and thus a better representation of the hadronic recoil, due to the more complete Emiss soft reconstruction and the lack
T
of a JVT-tagging requirement. This effect is partly due to the observation bias in the response introduced by the relatively poor Emiss calo resolution, as discussed in Sect. 3.1.
T
and Emiss calo
Both ETmiss T show comparable response for
123
T
Emiss calo miss
T and track-only-based pT response in an a exclusive and
an b inclusive Z → ͣ sample from data. The projections of the , Emiss calo miss Z
respective Emiss ,and ponto the direction of pT , calcu-
TT T
Z
pT 2 60 GeV, owing to the JVT cut-off at 60 GeV. The slightly larger Emiss calo response of about 1 GeV reflects
T the contribution from neutral signals to the soft term. The
degraded response associated with pTmiss related to the exclu
sion of hard objects is clearly visible in this figure. Emiss
Figure 16c shows the linearity of the various T reconstruction approaches as a function of Emiss true for
T the tt¯-production sample from MC simulations. Beyond Emiss true and Emiss calo
≈ 60 GeV both the reference Emiss
T TT show the same good linearity, while the lack of a jet con
lated according to Eqs. (9)and (10), are shown as a function of pT Z.In c, , Emiss calo miss
the linearity of the reference ETmiss T ,and pT scales, calculated according to Eq. (11), is shown as a function of the true Emiss true
T for the tt¯-production MC simulation sample
tribution to pmiss shows a loss of response up to about 50%
T at higher Emiss true . The overestimation of Emiss true by all
TT three reconstruction variants at lower ETmiss true reflects the observation bias in the response introduced by the resolution. The poorer resolution associated with pTmiss observed in Fig. 15c for this sample leads to a faster rise of the response with decreasing Emiss true than for the reference
T
Emiss and Emiss calo
T T , which show a very similar dependence on Emiss true
T.
123
8.3.3 Summary of performance
Both Emiss calo miss
T and pT offer alternative measures for
Emiss
. The calorimeter-based Emiss calo uses topo-clusters
TT
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 Emiss.This
T introduces a pile-up dependence into Emiss calo , due to the
T
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 Emiss reconstruction, to allow the cancel-
T lation of jet-like pT-flow from pile-up in Emiss soft calo by
T pile-up jets in its hard term. The Emiss calo response in the
T
inclusive Z → ͣ sample is better than the reference Emiss
T
response, in particular in the region of small hadronic recoil
Z
( pT 20 GeV). It is comparable to the reference in tt¯ final states. The observed RMSmiss, in particular in Z →
xy
without jets, is significantly more affected by pile-up than
miss
is the reference Emiss or the track-only-based p. In final
TT
t, Emiss calo
states with a considerable number of jets, like t¯T performs nearly as well as the reference Emiss, with a slight
T
degradation of the Emiss resolution at highest pile-up activi-
T
ties. This variant is useful for physics analyses least sensitive to the soft-term contribution to Emiss resolution but requiring
T
a linear Emiss response.
T
The track-only-based pmiss displays a degraded response
T
for the inclusive Z → ͣ 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 refer
miss
ence Emiss algorithm. Nevertheless, pprovides a stable
TT
observable for event and phase-space selections in analyses sensitive to ETmiss 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 s = 13 TeV and corresponding to about 3 2fb−1, are evaluated for selected event samples with (W → e, W → , tt¯) and without (Z → ) genuine Emiss. The comparison
T
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 ETmiss is determined with Z → ͣ 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 EmissT onto the vector sum of the transverse momenta of the hard objects phard. The parallel projections yield the
T uncertainty of the Emiss scale, evaluated as a function of the
T total transverse momentum of the hard objects (pThard). The
widths of the distributions of the parallel and perpendicular
projections yield the respective systematic uncertainties of the Emiss
resolution. Simulation tends to underestimate the
T perpendicular resolution and overestimate the scale and par-allel resolution, in each case differing from data by at most 10%. The performance evaluation of Emiss response and res-
T olution for the inclusive Z → ͣ sample shows that data and MC simulations agree within the systematic uncertainties. The Emiss response shows an underestimation of the soft
T contributions to Emiss. A degradation of the Emiss resolution
TT is observed for increasing L ET and NPV, due to pile-up and
detector resolution effects. Additional performance measures
considered in these studies include the estimate of tails in the Emiss
T distribution. As expected from the universality of the hadronic recoil, the integral tail fraction of the Emiss distri-
T bution is identical for inclusive Z and W boson production, independent of the leptonic decay mode. The tt¯ 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 Emiss, in terms of the inclusion of pile-up
T 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 Emiss reconstruction in ATLAS, which was
T developed for LHC Run 1 and used in a large number of physics analyses, can be used with the discussed refinements and adjustments for Run 2 as well.
Acknowledgements We thank CERN for the very successful operation of the LHC, as well as the support staff from our institutions without whom ATLAS could not be operated efficiently. We acknowledge the support of ANPCyT, Argentina; YerPhI, Armenia; ARC, Australia; BMWFW and FWF, Austria; ANAS, Azerbaijan; SSTC, Belarus; CNPq and FAPESP, Brazil; NSERC, NRC and CFI, Canada; CERN; CONICYT, Chile; CAS, MOST and NSFC, China; COLCIENCIAS, Colombia; MSMT CR, MPO CR and VSC CR, Czech Republic; DNRF and DNSRC, Denmark; IN2P3-CNRS, CEA-DRF/IRFU, France; SRNSFG, Georgia; BMBF, HGF, and MPG, Germany; GSRT, Greece; RGC, Hong Kong SAR, China; ISF, I-CORE and Benoziyo Center, Israel; INFN, Italy; MEXT and JSPS, Japan; CNRST, Morocco; NWO, Netherlands; RCN, Norway; MNiSW and NCN, Poland; FCT, Portugal; MNE/IFA, Romania; MES of Russia and NRC KI, Russian Federation; JINR; MESTD, Serbia; MSSR, Slovakia; ARRS and MIZŠ, Slovenia; DST/NRF, South Africa; MINECO, Spain; SRC and Wallenberg Foundation, Sweden; SERI, SNSF and Cantons of Bern and Geneva, Switzerland; MOST, Taiwan; TAEK, Turkey; STFC, United Kingdom; DOE and NSF, United States of America. In addition, individual groups and members have received support from BCKDF, the Canada Council, CANARIE, CRC, Compute Canada, FQRNT, and the Ontario Innovation Trust, Canada; EPLANET, ERC, ERDF, FP7, Horizon 2020 and Marie Skłodowska-Curie Actions, European Union; Investissements d’Avenir Labex and Idex, ANR, Région Auvergne and Fondation Partager le Savoir, France; DFG and AvH Foundation, Germany; Herakleitos, Thales and Aristeia programmes co-financed by
123
EU-ESF and the Greek NSRF; BSF, GIF and Minerva, Israel; BRF, Norway; CERCA Programme Generalitat de Catalunya, Generalitat Valenciana, Spain; the Royal Society and Leverhulme Trust, United Kingdom. The crucial computing support from all WLCG partners is acknowledged gratefully, in particular from CERN, the ATLAS Tier-1 facilities at TRIUMF (Canada), NDGF (Denmark, Norway, Sweden), CC-IN2P3 (France), KIT/GridKA (Germany), INFN-CNAF (Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Taiwan), RAL (UK) and BNL (USA), the Tier-2 facilities worldwide and large non-WLCG resource providers. Major contributors of computing resources are listed in Ref. [46].
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP3.
AppendixA:Glossaryof terms
In this paper several acronyms and qualifiers are used to describe the reconstruction of ETmiss and related observables. This brief glossary of terms is intended to help with the nomenclature. All terms should be interpreted in the context
of Emiss
T reconstruction and may have other interpretations in other contexts.
Emiss
T reconstruction: Using this nomenclature usually encompasses the reconstruction of a set of observables comprising the missing transverse momentum vector
Emiss
T , its components Exmiss y , and its absolute value ETmiss. In addition, the scalar sum L ET of the pT of all kinematic objects contributing to Emiss is calculated. The calcula-
T tion of these variables is described in detail in Eqs. (1)–(4) in Sect. 3.1. Hard scatter and primary vertices: The hardest pp interaction in a given event is referred to as the hard scatter. It is normally associated with a reconstructed hard-scatter vertex, which is considered the hardest vertex among all reconstructed primary vertices in this event. The hard-scatter vertex is defined as the one with the largest sum of pT2 of tracks associated with it. The other primary vertices are assumed to be produced by in-time pile-up interactions. The variable NPV denotes the number of reconstructed primary collision vertices in the event.
Hard event and hard term: The reconstruction of hard objects includes individual particles such as electrons, photons, muons and τ -leptons, and jets. In all cases the final objects are characterised by a kinematic threshold and reconstruction quality requirements. Both the reconstructed charged-particle tracks from the ID and topo-clusters from the calorimeter are used as the input signals for these objects. In the context
Emiss
of T reconstruction, the use of the same detector signals by different hard objects is excluded, see details in Sect. 3. The finally accepted hard event objects give rise to the hard term in Emiss reconstruc-
T tion. Soft event and soft term: All detector signals recorded for one triggered event and not used by the hard objects discussed above can be considered as soft signals contributing to Emiss. They include signals or signal traces
T from scattered soft particles arising from the underlying event accompanying the hard-scatter interaction, or from statistically completely independent pile-up interactions producing diffuse particle emissions in the same bunch crossing. In addition, signals from particles and jets which do not satisfy the hard-object quality criteria, or are below the kinematic threshold, can be included in the soft event. The reference ETmiss reconstruction configuration for the results presented in Sects. 6 and 7 uses reconstructed ID tracks from the soft event to form the soft term in ETmiss reconstruction, with the track selection details outlined under priority (6) in Table 1 in Sect. 3. Alternative configurations employ topo-clusters from the calorimeter, see Table 2 in Appendix A.
AppendixB:Alternative Emiss composition
T
Table 2 summarises the Emiss reconstruction configurations
T employing only ID tracks, or using topo-clusters for the soft
term.
123
Table2 Representations of Emiss and L ET calculated from (1) reconstructed charged-particle tracks from the ID or (2) using a soft term from
T
topo-clusters in the calorimeter only
# Objects contributing to Emiss and L ET
T
Type Selections Variables Comments
miss
(1) ID track || Ϫ 25 ppT Ϭ 400 MeV T pT
|d0| Ϫ 1 5mm |z0 sinω | Ϫ 1 5mm
Emiss soft calo
ELCW
(2) Topo-cluster >0(softterm) T
clus L Esoft calo T
Charged-particle-based estimators for Emiss,
T L ET using all reconstructed tracks from the hard-scatter vertex passing requirements for high-quality reconstruction in addition to kinematic selections
(L Ecalo
Variant reconstructing Emiss calo )
TT (L Esoft calo
using a soft term Emiss soft calo )
TT reconstructed from topo-clusters not used by, or not overlapping with, the hard objects used for the hard term composed of items (1)–(5) in Table 1, with the jet selection described in Sect. 8.1 applied
AppendixC:Jet selection
As discussed in Sect. 3.3.5, jets that are not rejected by the sig
nal ambiguity resolution and have pT Ϭ 60 GeV contribute to Emiss
reconstruction. Jets with less transverse momentum
T that fall within || Ϫ 2 4 are subjected to further selection based on JVT calculated by a track-based jet vertex tagger [16]. Three values for JVT, each representing a different efficiency for the reconstruction of non-pile-up jets, were considered in the course of the optimisation of the JVT-based selection,
JVTtight Ϭ 0 11 tight selection with high pile-up
rejection power at lower signal efficiency;
JVTmedium Ϭ 0 59 medium selection with good signal efficiency and pile-up rejection power;
JVTloose Ϭ 0 91 loose selection with lower pile-up
rejection power and higher signal efficiency.
Emiss
The effects of these selections on the resolution
T
RMSmiss
xy and response are shown in Fig. 17,for Z → ͣ events in MC simulation. Figure 17a shows that the pile-up dependence of RMSmiss is not significantly affected by
xy the choice for JVT. However, the Emiss response mea-
T sured by the projection given in Eq. (10) in Sect. 6.2.1
Z
and evaluated as a function of pT in the same sample, shows significant sensitivity to the choice of JVT, as seen in Fig. 17b.
In addition to the signal ambiguity resolution and the choice for JVT, the contribution from jets in Emiss recon-
T struction is controlled by a kinematic threshold requiring the
transverse momentum of the jet to be pT Ϭ 20 GeV. The effects of variations of this threshold on RMSmiss and the
xy
Emiss
T response are shown in Fig. 18. Increasing the thresh
old to 30 GeV for all jets satisfying the JVTmedium condition reduces the pile-up dependence of the resolution shown in Fig. 17a, but leads to significant loss of Emiss response, as
T
seen in Fig. 17b. Depending on the sensitivities observed in a given physics analysis, the pT-threshold choice for the jet contribution to Emiss reconstruction needs to be adjusted to
T
meet the required performance.
Extending the pT threshold studies with the option of regional thresholds yields the performance results presented in Fig. 19. In this case jets within || Ϫ 2 4 are sub-
central jet
jected to the pT Ϭ 20 GeV selection, while jets
forward jet
outside of this χ range are filtered using pT Ϭ {20 2530} GeV. This leads to the improvements in the pile-up dependence of RMSmiss shown in Fig. 19a, which are
xy
very similar to the ones observed in Fig. 18a for a global jet-pT threshold variation. The comparison indicates that the main pile-up contribution to RMSmiss is introduced by for
xy
ward jets, for which no JVT-based pile-up-mitigation is available.
Increasing the pT threshold only for forward jets reduces the average loss of response observed in case of the global pT threshold increase. This can be seen by comparing the results shown in Fig. 19b for regional pT thresholds with the ones shown in Fig. 18b for global thresholds. Like for the JVT threshold selection, the choice of the appropriate global or regional pT-threshold depends on the Emiss reconstruc-
T
tion performance required in the context of a given analysis.
123
Emiss Z
T discussed in the text, as measured for Z →ͣ events in MC simulation. global selection of pT Ϭ 20 GeV is applied to the transverse momentum The a Emiss central jet
T resolution is shown as function of the pile-up activity mea-of jets within ||Ϫ 24(p) and for forward jets with ||↖24
T sured by the number of primary vertices NPV,and the bEmiss response forward jet
T(pT)
Emiss Z
T central jet forward jet
shown as function of the pile-up activity measured by the number of ||Ϫ 24(p) and for forward jets with ||↖24(p),
TT primary vertices NPV,and the b ETmiss response is shown as function with pcentral jet pforward jet Ϭ {20 2530}GeV
TT
123
Emiss
T
shown as function of the pile-up activity measured by the number of primary vertices NPV,and the b Emiss response is shown as function
T
References
1.
ATLAS Collaboration, The ATLAS experiment at the CERN large hadron collider. JINST 3, S08003 (2008)
2.
ATLAS Collaboration, Performance of missing transverse momentum reconstruction in proton–proton collisions at 7 TeV with ATLAS. Eur. Phys. J. C 72, 1844 (2012). arXiv:1108.5602 [hep-ex]
3.
ATLAS Collaboration, Performance of algorithms that reconstruct
√
missing transverse momentum in s =8 TeV proton–proton collisions in the ATLAS detector. Eur. Phys. J. C 77, 241 (2017). arXiv:1609.09324 [hep-ex]
4.
ATLAS Collaboration, ATLAS insertable B-layer technical design report. ATLAS-TDR-19 (2010). https://cds.cern.ch/record/ 1291633 [ATLAS insertable B-layer technical design report addendum. ATLAS-TDR-19-ADD-1 (2012). https://cds.cern.ch/record/ 1451888]
5.
ATLAS Collaboration, Performance of the ATLAS trigger system in 2015. Eur. Phys. J. C 77, 317 (2017). arXiv:1611.09661 [hep-ex]
6.
ATLAS Collaboration, Topological cell clustering in the ATLAS calorimeters and its performance in LHC Run 1. Eur. Phys. J. C 77, 490 (2017). arXiv:1603.02934 [hep-ex]
7.
ATLAS Collaboration, Electron and photon energy calibration with the ATLAS detector using LHC Run 1 data. Eur. Phys. J. C 74, 3071 (2014) arXiv:1407.5063 [hep-ex]
8.
ATLAS Collaboration, Electron efficiency measurements with the ATLAS detector using the 2012 LHC proton–proton collision data. ATLAS-CONF-2014-032 (2014). https://cds.cern.ch/record/ 1706245
9.
ATLAS Collaboration, Measurement of the photon identification efficiencies with the ATLAS detector using LHC Run-1 data. Eur. Phys. J. C 76, 666 (2016). arXiv:1606.01813 [hep-ex]
10.
ATLAS Collaboration, Reconstruction of hadronic decay products of tau leptons with the ATLAS experiment. Eur. Phys. J. C 76, 295 (2016). arXiv:1512.05955 [hep-ex]
11.
ATLAS Collaboration, Identification and energy calibration of hadronically decaying tau leptons with the ATLAS experiment in
of the transverse momentum pT Z of the Z boson. The same threshold is applied to the transverse momentum of jets within || Ϫ 24
central jet
(pϬ 20 GeV), while for forward jets with ||↖2 4 threshold
T forward jet
variations (pϬ {20 25 30}GeV) are studied
T
√
pp collisions at s = 8 TeV. Eur. Phys. J. C 75, 303 (2015). arXiv:1412.7086 [hep-ex]
12. ATLAS Collaboration, Muon reconstruction performance of the
√
ATLAS detector in proton–proton collision data at s =13 TeV. Eur. Phys. J. C 76, 292 (2016). arXiv:1603.05598 [hep-ex]
13.
M. Cacciari, G.P. Salam, G. Soyez, The anti −kt jet clustering algorithm. JHEP 04, 063 (2008). arXiv:0802.1189 [hep-ph]
14.
M. Cacciari, G.P. Salam, G. Soyez, FastJet user manual. Eur. Phys.
J. C 72, 1896 (2012). arXiv:1111.6097
15.
ATLAS Collaboration, Jet energy scale measurements and their
√
systematic uncertainties in proton–proton collisions at s = 13 TeV with the ATLAS detector. Phys. Rev. D 96, 072002 (2017). arXiv:1703.09665 [hep-ex]
16. ATLAS Collaboration, Performance of pile-up mitigation tech
√
niques for jets in pp collisions at s =8 TeV using the ATLAS detector. Eur. Phys. J. C 76, 581 (2016). arXiv:1510.03823 [hep-ex]
17.
M. Cacciari, G.P. Salam, G. Soyez, The catchment area of jets. JHEP 04, 005 (2008). arXiv:0802.1188 [hep-ph]
18.
ATLAS Collaboration, Performance of the ATLAS track reconstruction algorithms in dense environments in LHC Run 2. Eur. Phys.J.C 77, 673 (2017). arXiv:1704.07983 [hep-ex]
19.
ATLAS Collaboration, Data-driven determination of the energy scale and resolution of jets reconstructed in the ATLAS calorime
√
ters using dijet and multijet events at s =8 TeV. ATLAS-CONF2015-017 (2015). https://cds.cern.ch/record/2008678
20.
S. Frixione, P. Nason, C. Oleari, Matching NLO QCD computations with Parton Shower simulations: the POWHEG method. JHEP 11, 070 (2007). arXiv:0709.2092 [hep-ph]
21.
T. Sjöstrand, S. Mrenna, P.Z. Skands, A brief introduction to PYTHIA 8.1. Comput. Phys. Commun. 178, 852 (2008). arXiv:0710.3820 [hep-ph]