PHYSICAL REVIEW D 98, 032009 (2018)
Search for supersymmetrypffiffi in events with four or more leptons in s =13 TeV pp collisions with ATLAS
M. Aaboud et al.* (ATLAS Collaboration)
(Received 11 April 2018; published 15 August 2018)
Results from a search for supersymmetry in events with four or more charged leptons (electrons, muons 36 1 −1 – and taus) are presented. The analysis uses a data samplepffiffiffi corresponding to . fb of proton proton collisions delivered by the Large Hadron Collider at s ¼ 13 TeV and recorded by the ATLAS detector. Four-lepton signal regions with up to two hadronically decaying taus are designed to target a range of supersymmetric scenarios that can be either enriched in or depleted of events involving the production and decay of a Z boson. Data yields are consistent with Standard Model expectations and results are used to set upper limits on the event yields from processes beyond the Standard Model. Exclusion limits are set at the 95% confidence level in simplified models of general gauge mediated supersymmetry, where Higgsino masses are excluded up to 295 GeV. In R-parity-violating simplified models with decays of the lightest supersymmetric particle to charged leptons, lower limits of 1.46, 1.06, and 2.25 TeV are placed on wino, slepton and gluino masses, respectively.
DOI: 10.1103/PhysRevD.98.032009
I. INTRODUCTION lifetime [11]. This conflict can be avoided by imposing ð−1Þ3ðB−LÞþ2S Supersymmetry (SUSY) [1–6] is a space-time symmetry the conservation of R-parity [12], defined as , that postulates the existence of new particles with spin where S is spin, or by explicitly conserving either B or L in R differing by one half-unit from their Standard Model (SM) the Lagrangian in -parity-violating (RPV) scenarios. In RPV models, the lightest SUSY particle (LSP) is unstable partners. In supersymmetric extensions of the SM, each SM and decays to SM particles, including charged leptons fermion (boson) is associated with a SUSY boson (fermion), and neutrinos when violating L but not B.InR-parity- having the same quantum numbers as its partner except for conserving (RPC) models, the LSP is stable and leptons can spin. The introduction of these new SUSY particles provides originate from unstable weakly interacting sparticles a potential solution to the hierarchy problem [7–10]. decaying into the LSP. Both the RPV and RPC SUSY The scalar superpartners of the SM fermions are called l˜ scenarios can therefore result in signatures with high lepton sfermions (comprising the charged sleptons, , the sneu- multiplicities and substantial missing transverse momen- ν˜ ˜ trinos, , and the squarks, q), while the gluons have tum, selections on which can be used to suppress SM ˜ fermionic superpartners called gluinos (g). The bino, wino background processes effectively. and Higgsino fields are fermionic superpartners of the This paper presents a search for new physics in final ð2Þ ð1Þ SU ×U gauge fields of the SM, and the two complex states with at least four isolated, charged leptons (electrons, scalar doublets of a minimally extended Higgs sector, muons or taus) where up to two hadronically decaying taus respectively. Their mass eigenstates are referred to as are considered. The analysis exploits the full proton–proton χ˜ ¼1 χ˜0 ¼1 charginos i (i , 2) and neutralinos j (j , 2, 3, 4), data set collected by the ATLAS experiment during the numbered in order of increasing mass. 2015 and 2016 data-taking periods, corresponding to an In the absence of a protective symmetry, SUSY proc- integrated luminosity of 36.1 fb−1 at a center-of-mass esses not conserving lepton number (L) and baryon number energy of 13 TeV. The search itself is optimized using (B) could result in proton decay at a rate that is in conflict several signal models but is generally model independent, with the tight experimental constraints on the proton using selections on the presence or absence of Z bosons in the event and loose requirements on effective mass or *Full author list given at end of the article. missing transverse momentum. Results are presented in terms of the number of events from new physics processes Published by the American Physical Society under the terms of with a four charged lepton signature, and also in terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to RPV and RPC SUSY models. the author(s) and the published article’s title, journal citation, Previous searches for SUSY particles using signatures and DOI. Funded by SCOAP3. with three or more leptons were carried out at the Tevatron
2470-0010=2018=98(3)=032009(31) 032009-1 © 2018 CERN, for the ATLAS Collaboration M. AABOUD et al. PHYS. REV. D 98, 032009 (2018)
– 0 collider [13 18], and at the LHC by the ATLAS experiment TABLE I. Decay modes and branching ratios for the χ˜1 LSP in [19–22] and the CMS experiment [23–27]. This analysis the RPV models, where ν denotes neutrinos or antineutrinos of closely follows the 7 TeV [19] and 8 TeV [22] ATLAS any lepton generation. analyses. 0 Scenario χ˜1 branching ratios ¯ 12 þ −ν 1 4 μ∓ν 1 2 μþμ−ν 1 4 II. SUSY SCENARIOS LLE k e e ( = ) e ( = ) ( = ) LLEi¯ 33 e τ∓ν (1=4) τþτ−ν (1=2) μ τ∓ν (1=4) SUSY models are used for signal region optimization and to interpret the results of this analysis. Models of both RPV SUSY and RPC SUSY are considered here, as they every signal event contains a minimum of four charged each require a different approach for signal selection, as leptons and two neutrinos. 1 λ χ˜0 discussed in Sec. V. In principle, the nine ijk RPV couplings allow the 1 to In all scenarios, the light CP-even Higgs boson, h, of the decay to every possible combination of charged lepton pairs, minimal supersymmetric extension of the SM [28,29] where the branching ratio for each combination differs for λ λ ≠ 0 Higgs sector is assumed to be practically identical to the each ijk. For example, for 121 the branching ratios for 0 0 0 SM Higgs boson [30], with the same mass and couplings as χ˜1 → eμν, χ˜1 → eeν and χ˜1 → μμν are 50%, 50% and 0% measured at the LHC [31–33]. In addition, the decoupling respectively, whereas for λ122 ≠ 0 the corresponding limit is used, which is defined by mA ≫ mZ, while the branching ratios are 50%, 0% and 50%. In Ref. [22],it CP-odd (A), the neutral CP-even (H), and the two charged was found that the four-charged-lepton search sensitivity is (H ) Higgs bosons are considered to be very heavy and comparable in the cases of λ121 ≠ 0 or λ122 ≠ 0, and for λ ≠ 0 λ ≠ 0 thus considerably beyond the kinematic reach of the LHC. 133 or 233 . Since the analysis reported here uses similar techniques, the number of L-violating RPV scenar- ios studied is reduced by making no distinction between the A. RPV SUSY scenarios 0 electron and muon decay modes of the χ˜1. Two extremes of In generic SUSY models with minimal particle content, the λijk RPV couplings are considered: the superpotential includes terms that violate conservation ¯ (i) LLE12k (k ∈ 1, 2) scenarios, where λ12k ≠ 0 and of L and B [34,35]: only decays to electrons and muons are included, ¯ 33 ∈ 1 λ ≠ 0 1 1 (ii) LLEi (i , 2) scenarios, where i33 and λ ¯ þ λ0 ¯ þ λ00 ¯ ¯ ¯ þκ only decays to taus and either electrons or muons are ijkLiLjEk ijkLiQjDk ijkUi Dj Dk iLiH2; 2 2 included. In both cases, all other RPV couplings are assumed to be where L and Q indicate the lepton and quark SU(2)- i i zero. The branching ratios for the χ˜0 decay in the LLE¯ 12k doublet superfields, respectively, and E¯ , U¯ and D¯ are the 1 i i i and LLEi¯ 33 are shown in Table I. The sensitivity to λ corresponding singlet superfields. Quark and lepton gen- couplings not considered here (e.g., λ123) is expected to be erations are referred to by the indices i, j and k, while the between that achieved in the LLE¯ 12k and LLEi¯ 33 scenarios. Higgs field that couples to up-type quarks is represented by 0 0 0 0 For the pure-bino χ˜1 considered here, the χ˜1χ˜1 production the Higgs SU(2)-doublet superfield H2. The λ , λ and ijk ijk cross section is found to be vanishingly small, thus models λ00 parameters are three sets of new Yukawa couplings, ijk that include one or more next-to-lightest SUSY particles κ while the i parameters have dimensions of mass. (NLSP) are considered in order to obtain a reasonably large Simplified models of RPV scenarios are considered, cross section. The choice of NLSP in the LLE¯ 12k and χ˜0 where the LSP is a bino-like neutralino ( 1) and decays via LLEi¯ 33 scenarios determines the cross section of the SUSY an RPV interaction. The LSP decay is mediated by the scenario, and can impact the signal acceptance to a lesser following lepton-number-violating superpotential term: extent. In all cases, the NLSP is pair produced in an RPC 1 interaction, and decays to the LSP (which itself undergoes an ¼ λ ¯ WLLE¯ 2 ijkLiLjEk: RPV decay). Three different possibilities are considered for the NLSP in the LLE¯ 12k and LLEi¯ 33 scenarios: This RPV interaction allows the following decay of the (i) Wino NLSP: Mass-degenerate wino-like charginos χ˜þχ˜− neutralino LSP: and neutralinos are produced in association ( 1 1 0 0 or χ˜1 χ˜2). The χ˜1 (χ˜ 2) decays to the LSP while emitting 0 ∓ χ˜1 → l l νj=i; ð1Þ a W (Z or h) boson, as shown in Figs. 1(a) and 1(b). k i=j l˜ ν˜ (ii) L= NLSP: Mass-degenerate left-handed sleptons through a virtual slepton or sneutrino, with the allowed and sneutrinos of all three generations are produced lepton flavors depending on the indices of the associated λ 1 ijk couplings [36]. The complex conjugate of the decay in The 27 λijk RPV couplings are reduced to 9 by the Eq. (1) is also allowed. Thus, in the case of pair production, antisymmetry requirement λijk ¼ −λjik.
032009-2 SEARCH FOR SUPERSYMMETRY IN EVENTS WITH FOUR … PHYS. REV. D 98, 032009 (2018)
(a) (b)
FIG. 2. Diagrams of the processes in the SUSY RPC GGM 0 Higgsino models. The W =Z produced in the χ˜1 =χ˜2 decays are off-shell (m ∼ 1 GeV) and their decay products are usually not reconstructed.
Planck-scale-mediated SUSY breaking scenario the grav- itino G˜ is the fermionic superpartner of the graviton, and its mass is comparable to the masses of the other SUSY particles, m ∼ 100 GeV [42,43]. General gauge mediated (GGM) SUSY models [44] predict the G˜ is nearly massless and offer an opportunity to study light Higgsinos. The FIG. 1. Diagrams of the benchmark SUSY models of RPC NLSP decays of the Higgsinos to the LSP G˜ would lead to on- pair production of (a) and (b) a wino, (c) slepton/sneutrino and 0 shell Z=h, and the decay products can be reconstructed. (d) gluino, followed by the RPV decay of the χ˜1 LSP. The LSP is assumed to decay as χ˜0 → llν with 100% branching ratio. Simplified RPC models inspired by GGM are considered 1 here, where the only SUSY particles within reach of the LHC are an almost mass-degenerate Higgsino triplet χ˜1 , l˜ l˜ ν˜ ν˜ l˜ ν˜ l˜ ν˜ 0 0 ˜ in association ( L L, , L ). The L ( ) decays to χ˜1, χ˜2 and a massless G. To ensure the SUSY decays are 0 0 the LSP while emitting a charged lepton (neutrino) prompt, the χ˜1 and χ˜2 masses are set to 1 GeVabove the χ˜1 as seen in Fig. 1(c). mass, and due to their weak coupling with the gravitino ˜ 0 (iii) g NLSP: Gluino pair production, where the gluino always decay to the χ˜1 via virtual Z=W bosons (which in decays to the LSP while emitting a quark-antiquark turn decay to very soft final states that are not recon- pair (u, d, s, c, b only, with equal branching ratios), 0 structed). The χ˜ 1 decays promptly to a gravitino plus a Z or as seen in Fig. 1(d). 0 h boson, χ˜ → Z=h þ G˜ , where the leptonic decays of the For the RPV models, the LSP mass is restricted to the range 1 Z=h are targeted in this analysis. Four production processes 10 GeV ≤ mðLSPÞ ≤ mðNLSPÞ − 10 GeV to ensure that þ − 0 are included in this Higgsino GGM model: χ˜1 χ˜1 , χ˜1 χ˜ 1, both the RPC cascade decay and the RPV LSP decay are χ˜ χ˜0 χ˜0χ˜0 χ˜0 1 2 and 1 2, as shown in Fig. 2, and the total SUSY cross prompt. Nonprompt decays of the 1 in similar models were 0 0 section is dominated by χ˜1 χ˜1 and χ˜1 χ˜2 production. The previously studied in Ref. [37]. 0 ˜ χ˜1 → ZG branching ratio is a free parameter of the GGM Higgsino scenarios, and so offers an opportunity to study B. RPC SUSY scenarios 4l signatures with one or more Z candidates. 0 0 RPC scenarios with light χ˜1, χ˜2 and χ˜1 Higgsino states are well motivated by naturalness [38,39]. However, they III. THE ATLAS DETECTOR can be experimentally challenging, as members of the 0 The ATLAS detector [45] is a multipurpose particle Higgsino triplet are close in mass and decays of the χ˜ 2=χ˜1 χ˜0 physics detector with forward-backward symmetric cylin- to a 1 LSP result in low-momentum decay products that 2 are difficult to reconstruct efficiently. Searches for drical geometry. The inner tracking detector (ID) covers Higgsino-like χ˜1 in approximately mass-degenerate sce- 2 narios were performed by the LEP experiments, where ATLAS uses a right-handed coordinate system with its origin chargino masses below 103.5 GeV were excluded [40] at the nominal interaction point (IP) in the center of the detector and the z-axis along the beam pipe. The x-axis points from the IP (reduced to 92 GeV for chargino-LSP mass differences to the center of the LHC ring, and the y-axis points upward. between 0.1 and 3 GeV). Recently, the ATLAS experiment Cylindrical coordinates ðr; ϕÞ are used in the transverse plane, ϕ has excluded Higgsino-like χ˜0 up to masses ∼145 GeV and being the azimuthal angle around the z-axis. The pseudorapidity 2 θ η ¼ − ðθ 2Þ down to χ˜0-LSP mass differences of 2.5 GeV [41] for is defined in terms of the polar angle as ln tan = . 2 Rapidity is defined as y ¼ 0.5 ln ½ðE þ p Þ=ðE − p Þ , where E scenarios where the χ˜ mass is assumed to be halfway z z 1 denotes the energy and pz is the component of the momentum between the two lightest neutralino masses. In the along the beam direction.
032009-3 M. AABOUD et al. PHYS. REV. D 98, 032009 (2018) jηj < 2.5 and consists of a silicon pixel detector, a semi- in perturbative QCD used for yield normalization, are conductor microstrip detector, and a transition radiation shown in Table II. Signal cross sections were calculated tracker. The innermostp pixelffiffiffi layer, the insertable B-layer to next-to-leading order in the strong coupling constant, [46], was added for the s ¼ 13 TeV running period of the adding the resummation of soft gluon emission at next-to- LHC. The ID is surrounded by a thin superconducting leading-logarithmic accuracy (NLO þ NLL) [48–55]. The solenoid providing a 2 T axial magnetic field. A high- nominal signal cross section and its uncertainty were taken granularity lead/liquid-argon sampling calorimeter mea- from an envelope of cross section predictions using differ- sures the energy and the position of electromagnetic ent parton distribution function (PDF) sets and factorization showers within jηj < 3.2. Sampling calorimeters with and renormalization scales, as described in Ref. [56]. liquid argon as the active medium are also used to measure The dominant irreducible background processes that can hadronic showers in the end cap (1.5 < jηj < 3.2) and produce four prompt and isolated charged leptons are ZZ, forward (3.1 < jηj < 4.9) regions, while a steel/scintillator ttZ¯ , VVV and Higgs production (where V ¼ W, Z, and tile calorimeter measures hadronic showers in the central includes off-shell contributions). For the simulated ZZ region (jηj < 1.7). The muon spectrometer (MS) surrounds production, the matrix elements contain all diagrams with the calorimeters and consists of three large superconducting four electroweak vertices, and they were calculated for up to air-core toroid magnets, each with eight coils, a system of one extra parton at NLO, and up to three extra partons at LO. precision tracking chambers (jηj < 2.7), and fast trigger The production of top quark pairs with an additional Z boson chambers (jηj < 2.4). A two-level trigger system [47] was simulated with the cross section normalized to NLO. selects events to be recorded for off-line analysis. Simulated triboson (VVV) production includes the processes ZZZ, WZZ and WWZ with four to six charged leptons, and IV. MONTE CARLO SIMULATION was generated at NLO with additional LO matrix elements for up to two extra partons. The simulation of Higgs Monte Carlo (MC) generators were used to simulate SM processes includes Higgs production via gluon-gluon fusion processes and new physics signals. The SM processes (ggH) and vector-boson fusion (VBF), and associated considered are those that can lead to signatures with at least production with a boson (WH, ZH) or a top-antitop pair four reconstructed charged leptons. Details of the signal (ttH¯ ). Other irreducible background processes with small and background MC simulation samples used in this cross sections are grouped into a category labeled “Other,” analysis, as well as the order of cross section calculations which contains the tWZ, ttWW¯ , tttW and ttt¯ t¯ processes.
TABLE II. Summary of the simulated SM background samples used in this analysis, where V ¼ W, Z, and includes off-shell contributions. “Tune” refers to the set of tuned parameters used by the generator. The sample marked with a † is used for a cross-check of yields and for studies of systematic uncertainties.
Cross-section Process Generator(s) Simulation calculation Tune PDF set
WZ, WW SHERPA 2.2.1 [57] Full NLO [58] SHERPA default NNPDF30NNLO [59] ZZ SHERPA 2.2.2 [57] Full NLO [58] SHERPA default NNPDF30NNLO [59]
VVV SHERPA 2.2.1 Full NLO [58] SHERPA default NNPDF30NNLO ggH, VBF, ggZH POWHEG v2 [60] + PYTHIA 8.186 [61] Full NNLO þ NNLL [62] AZNLO [63] CT10 [64] ZH, WH PYTHIA 8.186 Full NNLO þ NNLL A14 [65] NNPDF23LO ttH¯ MADGRAPH5_AMC@NLO 2.3.2 [66] Full NLO [67] A14 NNPDF23LO [68] + PYTHIA 8.186 ttZ¯ , ttW¯ , ttWW¯ MADGRAPH_5_AMC@NLO 2.2.2 [69] Full NLO [67] A14 NNPDF23LO + PYTHIA 8.186 † ttZ¯ SHERPA 2.2.1 Fast NLO [67] SHERPA default NNPDF30NNLO tWZ aMC@NLO 2.3.2 + PYTHIA 8.186 Full NLO [67] A14 NNPDF23LO tttW, ttt¯ t¯ MADGRAPH5_AMC@NLO 2.2.2 Full NLO [69] A14 NNPDF23LO + PYTHIA 8.186 tt¯ POWHEG v2 + PYTHIA 6.428 [70] Full NNLO þ NNLL [71] Perugia2012 [72] CT10
Z þ jets, W þ jets MADGRAPH5_AMC@NLO 2.2.2 Full NNLO [73] A14 NNPDF23LO + PYTHIA 8.186
SUSY signal MADGRAPH5_AMC@NLO 2.2.2 Fast NLO þ NLL [48–55] A14 NNPDF23LO + PYTHIA 8.186
032009-4 SEARCH FOR SUPERSYMMETRY IN EVENTS WITH FOUR … PHYS. REV. D 98, 032009 (2018)
For all MC simulation samples, the propagation of are rejected to suppress events with large calorimeter noise particles through the ATLAS detector was modeled with or noncollision backgrounds. GEANT4 [74] using the full ATLAS detector simulation The visible part of hadronically decaying tau leptons, τ [75], or a fast simulation using a parametrization of the denoted as had-vis and conventionally referred to as taus response of the electromagnetic and hadronic calorimeters throughout this paper, is reconstructed [82] using jets as jηj 2 47 10 [75] and GEANT4 elsewhere. The effect of multiple proton- described above with < . and pT > GeV. The τ proton collisions in the same or nearby bunch crossings, in- had-vis reconstruction algorithmpffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi uses information about the time and out-of-time pileup, is incorporated into the tracks within ΔR ≡ ðΔϕÞ2 þðΔηÞ2 ¼ 0.2 of the jet simulation by overlaying additional minimum-bias events direction, in addition to the electromagnetic and hadronic generated with PYTHIA8 [61] onto hard-scatter events. τ shower shapes in the calorimeters. Preselected had-vis Simulated events are reconstructed in the same manner candidates are required to have one or three associated as data, and are weighted to match the distribution of the tracks (prongs), because taus predominantly decay to either expected mean number of interactions per bunch crossing one or three charged hadrons together with a neutrino and τ in data. The simulated MC samples are corrected to account often additional neutral hadrons. The preselected had-vis are 20 for differences from the data in the triggering efficiencies, required to have pT > GeV and unit total charge of their lepton reconstruction efficiencies, and the energy and constituent tracks. In order to suppress electrons misiden- τ momentum measurements of leptons and jets. tified as preselected had-vis, taus are vetoed using transition radiation and calorimeter information. The preselected τ τ V. EVENT SELECTION had-vis candidates are corrected to the had-vis energy scale η using an - and pT-dependent calibration. A boosted After the application of beam, detector and data-quality decision tree algorithm (BDT) uses discriminating track requirements, the total integrated luminosity considered in and cluster variables to optimize τ identification, −1 had-vis this analysis corresponds to 36.1 1.2 fb . Events where “loose,”“medium” and “tight” working points are recorded during stable data-taking conditions are used in defined [83], but not used to preselect tau leptons. In this the analysis if the reconstructed primary vertex has at least analysis, kinematic variables built with hadronically 400 two tracks with transverse momentum pT > MeV decaying taus use only their visible decay products. associated with it. The primary vertex of an event is identified The missing transverse momentum, Emiss, is the magni- Σ 2 T as the vertex with the highest pT of associated tracks. tude of the negative vector sum of the transverse momenta Preselected electrons are required to have jηj < 2.47 and of all identified physics objects (electrons, photons, muons 7 η pT > GeV, where the pT and are determined from the and jets) and an additional soft term [84]. Taus are included calibrated clustered energy deposits in the electromagnetic miss as jets in the ET . The soft term is constructed from the calorimeter and the matched ID track, respectively. tracks matched to the primary vertex, but not associated Electrons must satisfy “loose” criteria of the likelihood- with identified physics objects, which allows the soft term based identification algorithm [76], with additional track to be nearly independent of pileup. requirements based on the innermost pixel layer. To avoid potential ambiguities among identified physics Preselected muons are reconstructed by combining tracks objects, preselected charged leptons and jets must survive in the ID with tracks in the MS [77], and are required to “overlap removal,” applied in the following order: jηj 2 7 5 have < . and pT > GeV. Muons must satisfy (1) Any tau within ΔR ¼ 0.2 of an electron or muon is “medium” identification requirements based on the number removed. of hits in the different ID and MS subsystems, and the (2) Any electron sharing an ID track with a muon is significance of the charge-to-momentum ratio, defined in removed. Ref. [77]. Events containing one or more muons that have a (3) Jets within ΔR ¼ 0.2 of a preselected electron are transverse impact parameter relative to the primary vertex discarded. jd0j > 0.2 mm or a longitudinal impact parameter relative (4) Electrons within ΔR ¼ 0.4 of a preselected jet are to the primary vertex jz0j > 1 mm are rejected to suppress discarded, to suppress electrons from semileptonic the cosmic-ray muon background. decays of c- and b-hadrons. Jets are reconstructed with the anti-kt algorithm [78] (5) Jets with fewer than three associated tracks are with a radius parameter of R ¼ 0.4. Three-dimensional discarded either if a preselected muon is within calorimeter energy clusters are used as input to the jet ΔR ¼ 0.2 or if the muon can be matched to a track reconstruction, and jets are calibrated following Ref. [79]. associated with the jet. jηj 2 8 20 Jets must have < . and pT > GeV. To reduce (6) Muons within ΔR ¼ 0.4 of a preselected jet are 60 jηj 2 4 pileup effects, jets with pT < GeV and < . must discarded to suppress muons from semileptonic satisfy additional criteria using the jet vertex tagging decays of c- and b-hadrons. algorithm described in Ref. [80]. Events containing jets (7) Jets within ΔR ¼ 0.4 of a preselected tau passing failing to satisfy the quality criteria described in Ref. [81] medium identification requirements are discarded.
032009-5 M. AABOUD et al. PHYS. REV. D 98, 032009 (2018)
TABLE III. The triggers used in the analysis of 2015 and 2016 VI. SIGNAL REGIONS data. The off-line pT thresholds are required only for recon- μ τ structed charged leptons which match to the trigger signatures. Events with four or more signal leptons (e, , had-vis) are Trigger thresholds for data recorded in 2016 are higher than in selected and are classified according to the number of light 2015 due to the increase in beam luminosity, and “or” denotes a signal leptons (L ¼ e, μ) and signal taus (T) required: at move to a higher-threshold trigger during data taking. least four light leptons and exactly zero taus 4L0T, exactly three light leptons and at least one tau 3L1T, or exactly two Off-line pT threshold [GeV] light leptons and at least two taus 2L2T. Trigger 2015 2016 Events are further classified according to whether they Single isolated e 25 27 are consistent with a leptonic Z boson decay or not. The Z Single nonisolated e 61 61 requirement selects events where any SFOS LL pair Single isolated μ 21 25 or 27 combination has an invariant mass close to the Z boson Single nonisolated μ 41 41 or 51 mass, in the range 81.2–101.2 GeV. A second Z candidate LL Double e 13, 13 18, 18 may be identified if a second SFOS pair is present and satisfies 61.2