Abstract
Two searches for new phenomena in final states containing a same-flavour opposite-sign lepton (electron or muon) pair, jets, and large missing transverse momentum are presented. These searches make use of proton–proton collision data, collected during 2015 and 2016 at a centre-of-mass energy \(\sqrt{s}=13\) \(\text {TeV}\) by the ATLAS detector at the large hadron collider, which correspond to an integrated luminosity of \(14.7~\mathrm {fb}^{-1}\). Both searches target the pair production of supersymmetric particles, squarks or gluinos, which decay to final states containing a same-flavour opposite-sign lepton pair via one of two mechanisms: a leptonically decaying Z boson in the final state, leading to a peak in the dilepton invariant-mass distribution around the Z boson mass; and decays of neutralinos (e.g. \(\tilde{\chi }_2^0 \rightarrow \ell ^+\ell ^- \tilde{\chi }_1^0\)), yielding a kinematic endpoint in the dilepton invariant-mass spectrum. The data are found to be consistent with the Standard Model expectation. Results are interpreted in simplified models of gluino-pair (squark-pair) production, and provide sensitivity to gluinos (squarks) with masses as large as 1.70 \(\text {TeV}\) (980 \(\text {GeV}\)).
Similar content being viewed by others
1 Introduction
Supersymmetry (SUSY) [1,2,3,4,5,6,7] is an extension of the Standard Model (SM) that introduces partner particles (called sparticles) that differ by half a unit of spin from their SM counterparts. The squarks (\(\tilde{q}\)) and sleptons (\(\tilde{\ell }\)) are the scalar partners of the quarks and leptons, respectively, and the gluinos (\(\tilde{g}\)) are the fermionic partners of the gluons. The charginos (\({\tilde{\chi }}_{i}^{\pm }\)) and neutralinos (\({\tilde{\chi }}_{i}^{0}\)) are the mass eigenstates (where the index i is ordered from the lightest to the heaviest) formed from the linear superpositions of the SUSY partners of the Higgs bosons (higgsinos) and electroweak gauge bosons.
If the masses of the gluino, higgsinos, and top squarks are close to the \(\text {TeV}\) scale, SUSY may offer a solution to the SM hierarchy problem [8,9,10,11]. In this case, strongly interacting sparticles should be produced at a high enough rate to be detected by the experiments at the large hadron collider (LHC). For models with R-parity conservation [12], such sparticles would be pair-produced and are expected to decay into jets, perhaps leptons, and the lightest stable SUSY particle (LSP). The LSP is assumed to be only weakly interacting and therefore escapes the detector, resulting in events with potentially large missing transverse momentum (\({\varvec{p}}_{\mathrm {T}}^\mathrm {miss}\), with magnitude \(E_{\text {T}}^{\text {miss}}\)). In such a scenario the LSP could be a dark-matter candidate [13, 14].
Final states containing pairs of leptons may arise from the cascade decays of squarks and gluinos via several mechanisms. In this paper, two search channels are considered that target scenarios with same-flavour (SF) opposite-sign (OS) lepton (electron or muon) pairs. The first channel requires a lepton pair with an invariant mass \(m_{\ell \ell }\) that is consistent with the Z boson mass \(m_Z\) (“on-shell Z” channel), while the second channel considers all SFOS lepton pairs (“edge” channel). The presence of two leptons in the final state suppresses large SM backgrounds from, e.g., QCD multijet and \(W+\mathrm {jets}\) production, providing a clean environment in which to search for new physics. As discussed further below, in such events the distribution of dilepton mass \(m_{\ell \ell }\) may be used to characterise the nature of the SUSY particle decay and constrain mass differences between SUSY particles.
The SFOS lepton pairs may be produced in the decay \(\tilde{\chi }_{2}^{0} \rightarrow \ell ^{+}\ell ^{-} \tilde{\chi }_{1}^{0}\) (or, in models of generalised gauge mediation with a gravitino LSP [15,16,17], via \(\tilde{\chi }_{1}^{0} \rightarrow \ell ^{+}\ell ^{-} \tilde{G}\)). The properties of the \(\tilde{\chi }_2^0\) decay depend on the mass difference \(\Delta m_\chi \equiv m_{\tilde{\chi }_{2}^{0}} - m_{\tilde{\chi }_{1}^{0}}\), the mixing of the charginos and neutralinos, and on whether there are additional sparticles with masses less than \(m_{\tilde{\chi }_{2}^{0}}\) that may be produced in the decay of the \(\tilde{\chi }_2^0\) particle. For \(\Delta m_\chi >m_Z\), SFOS lepton pairs may be produced in the decay \(\tilde{\chi }_{2}^{0} \rightarrow Z \tilde{\chi }_{1}^{0} \rightarrow \ell ^{+}\ell ^{-} \tilde{\chi }_{1}^{0}\), leading to a peak in the invariant-mass distribution near \(m_{\ell \ell } \approx m_Z\). Such models are the target of the on-shell Z search. For \(\Delta m_\chi < m_Z\), the decay \(\tilde{\chi }_{2}^{0} \rightarrow Z^* \tilde{\chi }_{1}^{0} \rightarrow \ell ^{+}\ell ^{-} \tilde{\chi }_{1}^{0}\) leads to a rising \(m_{\ell \ell }\) distribution that is truncated at a kinematic endpoint, whose position is given by \(m_{\ell \ell }^{\text {max}}=\Delta m_\chi < m_Z\), below the Z boson mass peak. If there are sleptons with masses less than \(m_{\tilde{\chi }_{2}^{0}}\), the \(\tilde{\chi }_2^0\) particle may decay as \(\tilde{\chi }_{2}^{0} \rightarrow \tilde{\ell }^{\pm }\ell ^{\mp } \rightarrow \ell ^{+}\ell ^{-} \tilde{\chi }_{1}^{0}\), also leading to a kinematic endpoint but with a different shape and \(m_{\ell \ell }\) endpoint position, given by \(m_{\ell \ell }^{\mathrm {max}} = \sqrt{ (m^2_{\tilde{\chi }_2^0}-m^2_{\tilde{\ell }})(m^2_{\tilde{\ell }}-m^2_{\tilde{\chi }_1^0}) / m^2_{\tilde{\ell }}}\), which may occur below, on, or above the Z boson mass peak. The latter two scenarios are targeted by the “edge” search channel, which considers the full \(m_{\ell \ell }\) range.
This paper reports on a search for SUSY in the same-flavour dilepton final state with \(14.7~\mathrm {fb}^{-1}\) of pp collision data at \(\sqrt{s}=13\) \(\text {TeV}\) recorded in 2015 and 2016 by the ATLAS detector at the LHC. Searches for SUSY in the \(Z+\mathrm {jets}+E_{\text {T}}^{\text {miss}} \) final state have previously been performed at \(\sqrt{s}=8\) \(\text {TeV}\) by the CMS [18, 19] and ATLAS [20] collaborations using Run-1 LHC data. In the ATLAS analysis performed with 20.3 fb\(^{-1}\) of \(\sqrt{s}=8\) \(\text {TeV}\) data reported in Ref. [20], an excess of events above the SM background with a significance of 3.0 standard deviations was observed. The event selection criteria for the on-shell Z search in this paper are almost identical, differing only in the details of the analysis object definitions and missing transverse momentum. CMS performed a search with \(\sqrt{s}=13\) \(\text {TeV}\) data in a similar kinematic region but did not observe evidence to corroborate this excess [21].
Searches for an edge in the \(m_{\ell \ell }\) distribution in events with \(2\ell +\mathrm {jets}+E_{\text {T}}^{\text {miss}} \) have been performed by the CMS [19, 22] and ATLAS [20] collaborations. In Ref. [19], CMS reported an excess above the SM prediction with a significance of 2.6 standard deviations. In a similar search region, however, the Run-1 ATLAS analysis [20] and Run-2 CMS analysis [21] observed results consistent with the SM prediction.
2 ATLAS detector
The ATLAS detector [23] is a general-purpose detector with almost \(4\pi \) coverage in solid angle.Footnote 1 The detector comprises an inner tracking detector, a system of calorimeters, and a muon spectrometer.
The inner tracking detector (ID) is immersed in a 2 T magnetic field provided by a superconducting solenoid and allows charged-particle tracking out to \(|\eta |=2.5\). It includes silicon-pixel and silicon-strip tracking detectors inside a straw-tube tracking detector. In 2015 the detector received a new innermost layer of silicon pixels, which improves the track impact parameter resolution by almost a factor of two in both the transverse and longitudinal directions [24].
High-granularity electromagnetic and hadronic calorimeters cover the region \(|\eta |<4.9\). All the electromagnetic calorimeters, as well as the endcap and forward hadronic calorimeters, are sampling calorimeters with liquid argon as the active medium and lead, copper, or tungsten as the absorber. The central hadronic calorimeter is a sampling calorimeter with scintillator tiles as the active medium and steel as the absorber.
The muon spectrometer uses several detector technologies to provide precision tracking out to \(|\eta |=2.7\) and triggering in \(|\eta |<2.4\), making use of a system of three toroidal magnets.
The ATLAS detector incorporates a two-level trigger system, with the first level implemented in custom hardware and the second level implemented in software. This trigger system selects events of interest at an output rate of about 1 kHz.
3 SUSY signal models
SUSY-inspired simplified models are considered as signal scenarios for these analyses. In all of these models, squarks or gluinos are directly pair-produced, decaying via an intermediate neutralino, \(\tilde{\chi }_2^0\), into the LSP (\(\tilde{\chi }_1^0\)). All sparticles not directly involved in the decay chains considered are effectively decoupled. Two example decay topologies are shown in Fig. 1. For all models with gluino-pair production, a three-body decay for \(\tilde{g}\rightarrow q \bar{q} \tilde{\chi }_2^0\) is used. Signal models are generated in a grid over a two-dimensional space, varying the gluino or squark mass and the mass of either the \(\tilde{\chi }_2^0\) or the \(\tilde{\chi }_1^0\).
Three models, one with squark-pair production and two with gluino-pair production, which result exclusively in events with two on-shell Z bosons in the final state are considered for the on-shell search. For two of these models, signal mass points are generated across the \(\tilde{g}\)–\(\tilde{\chi }_2^0\) (or \(\tilde{q}\)–\(\tilde{\chi }_2^0\)) plane. These models are produced following the decays \(\tilde{g}\rightarrow q \bar{q} \tilde{\chi }_2^0\) or \(\tilde{q}\rightarrow q \tilde{\chi }_2^0\), with the \(\tilde{\chi }_1^0\) (LSP) mass set to 1 \(\text {GeV}\), inspired by SUSY scenarios with a low-mass LSP (e.g. generalised gauge mediation). These two models are referred to here as the \(\tilde{g}\)–\(\tilde{\chi }_2^0\) on-shell and \(\tilde{q}\)–\(\tilde{\chi }_2^0\) on-shell grids, respectively, and are summarised in Table 1. The third model is based on MSSM-inspired topologies [25,26,27] with potentially higher mass LSPs. Signal points are generated across the \(\tilde{g}\)–\(\tilde{\chi }_1^0\) plane, and this model is thus referred to as the \(\tilde{g}\)–\(\tilde{\chi }_1^0\) on-shell grid. In this case the \(\tilde{\chi }_2^0\) mass is set to be 100 \(\text {GeV}\) above the \(\tilde{\chi }_1^0\) mass, which in many models maximises the branching fraction of the \(\tilde{\chi }_2^0\) decay to Z bosons. For the two models with gluino-pair production, since the gluino coupling to \(q\tilde{q}\) is flavour independent and the corresponding flavours of squarks are assumed to be mass degenerate, the branching fractions for \(q=u,d,c,s\) are each 25%. Other ATLAS searches are dedicated to final states with two leptons and heavy flavour jets [28, 29]. For the model involving squark-pair production, the superpartners of the u, d, c and s quarks have the same mass, with the superpartners of the b and t quarks being decoupled.
The edge search considers two scenarios, both of which involve the direct pair production of gluinos and differ by the decay mode of the \(\tilde{\chi }_{2}^{0}\). These signal models are also summarised in Table 1. In the \(Z^{(*)}\) model the \(\tilde{\chi }_{2}^{0}\) decays as \(\tilde{\chi }_{2}^{0} \rightarrow Z^{(*)}\tilde{\chi }_{1}^{0}\). For \(\Delta m_\chi = m(\tilde{\chi }_2^0) - m(\tilde{\chi }_1^0) > m_Z\), the Z boson is on-shell, leading to a peak in the \(m_{\ell \ell }\) distribution at \(m_Z\), while for \(\Delta m_\chi < m_Z\), the Z boson is off-shell, leading to an edge in the dilepton mass distribution with a position below \(m_Z\). The slepton model assumes that the sleptons are lighter than the \(\tilde{\chi }_{2}^{0}\), which decays as \(\tilde{\chi }_{2}^{0} \rightarrow \tilde{\ell }^{\mp }\ell ^{\pm }\) with \(\tilde{\ell } \rightarrow \ell \tilde{\chi }_{1}^{0}\) or as \(\tilde{\chi }_{2}^{0} \rightarrow \tilde{\nu }\nu \) with \(\tilde{\nu } \rightarrow \nu \tilde{\chi }_{1}^{0}\), each with a branching fraction of 50%, where \(\tilde{\ell }=\tilde{e},\tilde{\mu },\tilde{\tau }\) and \(\tilde{\nu }=\tilde{\nu _e},\tilde{\nu _\mu },\tilde{\nu _\tau }\). The endpoint position can occur at any mass, highlighting the need to search over the full dilepton mass distribution. The gluino decays as \(\tilde{g}\rightarrow q \bar{q} \tilde{\chi }_2^0\), and both models have equal branching fractions for \(q=u,d,c,s,b\). The \(\tilde{\chi }_2^0\) mass is set to the average of the gluino and \(\tilde{\chi }_1^0\) masses. For the slepton model, the masses of the superpartners of the left-handed leptons are set as the average of the \(\tilde{\chi }_2^0\) and \(\tilde{\chi }_1^0\) masses, while the superpartners of the right-handed leptons are decoupled. The three slepton flavours are mass-degenerate. In both these models the \(\tilde{g}\) and \(\tilde{\chi }_1^0\) masses are free parameters that are varied to produce the two-dimensional signal grid. The mass splittings are chosen to maximise the differences between these simplified models and other models with only one intermediate particle between the gluino and the LSP [30].
4 Data and Monte Carlo samples
The data used in this analysis were collected by ATLAS during 2015 and 2016, with a mean number of additional pp interactions per bunch crossing (pile-up) of approximately 14 in 2015 and 21 in 2016, and a centre-of-mass collision energy of 13 \(\text {TeV}\). Following requirements based on beam and detector conditions and data quality, the data set corresponds to an integrated luminosity of \(14.7~\mathrm {fb}^{-1}\). The uncertainty in the combined 2015 and 2016 integrated luminosity is \(\pm 2.9\)%. It is derived, following a methodology similar to that detailed in Refs. [31] and [32], from a preliminary calibration of the luminosity scale using x–y beam-separation scans performed in August 2015 and May 2016.
Data events are collected using a combination of single-lepton and dilepton triggers [33], in order to maximise the signal acceptance. The dielectron, dimuon, and electron–muon triggers have leading-lepton \(p_{\text {T}} \) thresholds in the range 12–24 \(\text {GeV}\). Additional single-electron (single-muon) triggers are also used, with trigger \(p_{\text {T}} \) thresholds of 60 (50) \(\text {GeV}\), to increase the trigger efficiency for models with high-\(p_{\text {T}}\) leptons. Events are required to contain at least two selected leptons with \(p_{\text {T}} >25\) \(\text {GeV}\), making the selection fully efficient with respect to the trigger \(p_{\text {T}}\) thresholds.
An additional control sample of events containing photons is collected using a set of single-photon triggers with \(p_{\text {T}}\) thresholds in the range 20–140 \(\text {GeV}\). All triggers except for the one with threshold \(p_{\text {T}} =120\) \(\text {GeV}\) in 2015, or the one with \(p_{\text {T}} =140\) \(\text {GeV}\) in 2016, are prescaled. Events are required to contain a selected photon with \(p_{\text {T}} >37\) \(\text {GeV}\), so that they are selected efficiently by the lowest available \(p_{\text {T}}\) trigger in 2015, which had a threshold of \(p_{\text {T}} ^{\gamma }=35\) \(\text {GeV}\).
Simulated event samples are used to aid in the estimation of SM backgrounds, validate the analysis techniques, optimise the event selection, and provide predictions for SUSY signal processes. All SM background samples used are listed in Table 2, along with the parton distribution function (PDF) set, the configuration of underlying-event and hadronisation parameters (underlying-event tune) and the cross-section calculation order in \(\alpha _{\text {S}}\) used to normalise the event yields for these samples.
Samples simulated using MG5_aMC@NLO v2.2.2 [34], interfaced with Pythia 8.186 [35] with the A14 underlying-event tune [36] to simulate the parton shower and hadronisation, are generated at leading order in \(\alpha _{\text {S}}\) (LO) with the NNPDF23LO PDF set [37]. For samples generated using Powheg Box V2 [38,39,40], Pythia 6.428 [41] is used to simulate the parton shower, hadronisation, and the underlying event. The CTEQ6L1 PDF set is used with the corresponding Perugia2012 [42] tune. In the case of both the MG5_aMC@NLO and Powheg samples, the EvtGen v1.2.0 program [43] is used for properties of the bottom and charm hadron decays. Sherpa 2.1.1 [44] simulated samples use the CT10 PDF set with Sherpa’s own internal parton shower [45] and hadronisation methods, as well as the Sherpa default underlying-event tune. Diboson processes with four charged leptons, three charged leptons and a neutrino or two charged leptons and two neutrinos are simulated using the Sherpa 2.1.1 generator. Matrix elements contain all diagrams with four electroweak vertices. They are calculated for up to one (\(4\ell \), \(2\ell +2\nu \)) or zero (\(3\ell +1\nu \)) partons at next-to-leading order in \(\alpha _{\text {S}}\) (NLO) and up to three partons at LO using the Comix [46] and OpenLoops [47] matrix element generators and merged with the Sherpa parton shower using the ME+PS@NLO prescription [48]. For the \(Z/\gamma ^{*}+\text {jets}\) background, Sherpa 2.1.1 is used to generate a sample with up to two additional partons at NLO and up to four at LO. For Monte Carlo (MC) closure studies, \(\gamma +\text {jets}\) events are generated at LO with up to four additional partons using Sherpa 2.1.1. Additional MC simulation samples of events with a leptonically decaying vector boson and photon (\(V\gamma \), where \(V=W,Z\)) are generated at LO using Sherpa 2.1.1. Matrix elements including all diagrams with three electroweak couplings are calculated with up to three partons. These samples are used to estimate backgrounds with real \(E_{\text {T}}^{\text {miss}}\) in \(\gamma +\text {jets} \) event samples.
The SUSY signal samples are produced at LO using MG5_aMC@NLO with the NNPDF2.3LO PDF set, interfaced with Pythia 8.186. The scale parameter for CKKW-L matching [49, 50] is set at a quarter of the mass of the gluino. Up to one additional parton is included in the matrix element calculation. The underlying event is modelled using the A14 tune for all signal samples, and EvtGen is adopted to describe the properties of bottom and charm hadron decays. Signal cross sections are calculated at NLO in \(\alpha _{\text {S}}\). This includes the resummation of soft gluon emission at next-to-leading-logarithm accuracy (NLO+NLL) [51,52,53,54,55].
All of the SM background MC samples are subject to a full ATLAS detector simulation [56] using GEANT4 [57]. A fast simulation [56], which uses a combination of a parameterisation of the response of the ATLAS electromagnetic and hadronic calorimeters and GEANT4, is used in the case of signal MC samples. This fast simulation is validated by comparing a few chosen signal samples to some fully simulated points. Minimum-bias interactions are generated and overlaid on the hard-scattering process to simulate the effect of multiple pp interactions occurring during the same (in-time) or a nearby (out-of-time) bunch-crossing (pile-up). These are produced using Pythia 8.186 with the A2 tune [58] and MSTW 2008 PDF set [59]. The pile-up distribution in MC samples is simulated to match that in data during 2015 and 2016 pp data-taking.
5 Analysis object identification and selection
All analysis objects are categorised as either “baseline” or “signal” based on various quality and kinematic requirements. Baseline objects are used in the calculation of missing transverse momentum and to disambiguate between the analysis objects in the event, while the jets and leptons entering the final analysis selection must pass more stringent signal requirements. The selection criteria for both the baseline and signal objects differ from the requirements used in the Run-1 ATLAS \(Z+\mathrm {jets}+E_{\text {T}}^{\text {miss}} \) search reported in Ref. [20], owing to the new silicon-pixel tracking layer and significant changes to the reconstruction software since 2012 data-taking. In particular, improvements in the lepton identification criteria have reduced the background due to hadrons misidentified as electrons. The primary vertex in each event is defined as the reconstructed vertex [74] with the highest \(\sum p_{\text {T}}^2\), where the summation includes all particle tracks with \(p_{\text {T}} >400\) \(\text {MeV}\) associated with the vertex.
Electron candidates are reconstructed from energy clusters in the electromagnetic calorimeter matched to ID tracks. Baseline electrons are required to have transverse energy \(E_{\text {T}}>10\) \(\text {GeV}\), satisfy the “loose likelihood” criteria described in Ref. [75] and reside within the region \(|\eta |<2.47\). Signal electrons are further required to have \(p_{\text {T}} >25\) \(\text {GeV}\), satisfy the “medium likelihood” criteria of Ref. [75], and be consistent with originating from the primary vertex. The signal electrons must originate from within \(|z_0\sin \theta | = 0.5\) mm of the primary vertex along the direction of the beamline.Footnote 2 The transverse-plane distance of closest approach of the electron to the beamline, divided by the corresponding uncertainty, must be \(|d_0/\sigma _{d_0}|<5\). These electrons must also be isolated with respect to other objects in the event, according to a \(p_{\text {T}}\)-dependent isolation requirement. The isolation uses calorimeter- and track-based information to obtain 95% efficiency at \(p_{\text {T}} =25\) \(\text {GeV}\), rising to 99% efficiency at \(p_{\text {T}} =60\) \(\text {GeV}\).
Baseline muons are reconstructed from either ID tracks matched to muon segments (collections of hits in a single muon spectrometer layer) or combined tracks formed from the ID and muon spectrometer [76]. They must satisfy the “medium” selection criteria described in Ref. [76], and to satisfy \(p_{\text {T}}>10\) \(\text {GeV}\) and \(|\eta |<2.5\). Signal muon candidates are further required to have \(p_{\text {T}} >25\) \(\text {GeV}\), be isolated, and have \(|z_0\sin \theta | < 0.5\) mm and \(|d_0/\sigma _{d_0}|<3\). Calorimeter- and track-based isolation criteria are used to obtain 95% efficiency at \(p_{\text {T}} =25\) \(\text {GeV}\), rising to 99% efficiency at \(p_{\text {T}} =80\) \(\text {GeV}\) [76]. Further, the relative uncertainties in the q / p of each of the ID track alone and muon spectrometer track alone are required to be less than 80% of the uncertainty in the q / p of the combined track. This reduces the already low rate of grossly mismeasured muons. The combined isolation and identification efficiency for single leptons, after the trigger requirements, is about 70% (80%) for electrons (muons) with \(p_{\text {T}} \sim 25\) \(\text {GeV}\), rising to about 90% for \(p_{\text {T}} >200\) \(\text {GeV}\).
Jets are reconstructed from topological clusters of energy [77] in the calorimeter using the anti-\(k_{t}\) algorithm [78, 79] with a radius parameter of 0.4. Calibration corrections are applied to the jets based on a comparison to jets made of stable particles (those with lifetimes \(\tau > 0.3 \times 10^{-10}\) s) in the MC simulation. A residual correction is applied to jets in data, based on studies of \(p_{\text {T}}\) balance between jets and well-calibrated objects in the MC simulation and data [80, 81]. Baseline jet candidates are required to have \(p_{\text {T}} >20\) \(\text {GeV}\) and reside within the region \(|\eta |<4.5\). Signal jets are further required to satisfy \(p_{\text {T}} >30\) \(\text {GeV}\) and reside within the region \(|\eta |<2.5\). Jets with \(p_{\text {T}} <60\) \(\text {GeV}\) and \(|\eta |<2.4\) must meet additional criteria designed to select jets from the hard-scatter interaction and reject those originating from pile-up. This is enforced by using the jet vertex tagger described in Ref. [82]. Finally, events containing a jet that does not pass specific jet quality requirements are vetoed from the analysis selection in order to remove events impacted by detector noise and non-collision backgrounds [83, 84]. The MV2c10 boosted decision tree algorithm [85, 86] identifies jets with \(|\eta |<2.5\) containing b-hadrons (b-jets) based on quantities such as the impact parameters of associated tracks and any reconstructed secondary vertices. A selection that provides 77% efficiency for tagging b-jets in simulated \(t\bar{t}\) events is used. These tagged jets are called b-tagged jets.
Photon candidates must satisfy “tight” selection criteria described in Ref. [87], have \(p_{\text {T}} >25\) \(\text {GeV}\) and reside within the region \(|\eta |<2.37\), excluding the transition region \(1.37<|\eta |<1.6\) where there is a discontinuity in the calorimeter. Signal photons are further required to have \(p_{\text {T}} >37\) \(\text {GeV}\) and to be isolated from other objects in the event, using \(p_{\text {T}}\)-dependent requirements on both track- and calorimeter-based isolation.
To avoid the duplication of analysis objects in more than one baseline selection, an overlap removal procedure is applied. Any baseline jet within \(\Delta R=0.2\) of a baseline electron is removed, unless the jet is b-tagged, in which case the electron is identified as originating from a heavy-flavour decay and is removed instead. Remaining electrons residing within \(\Delta R=0.4\) of a baseline jet are then removed from the event. Subsequently, any baseline muon residing within \(\Delta R=0.2\) of a remaining baseline b-tagged jet is discarded. If such a jet is not b-tagged then the jet is removed instead. Any remaining muon found within \(\mathrm { min } (0.04+(10~\text {GeV})/p_{\text {T}},0.4)\) of a jet is also discarded. This stage of the overlap removal procedure differs from that used in Ref. [20]. It was improved to retain muons near jet candidates mostly containing calorimeter energy from final-state radiation from muons, while still rejecting muons from heavy-flavour decays. Finally, to remove electron candidates originating from muon bremsstrahlung, any baseline electron within \(\Delta R=0.01\) of any remaining baseline muon is removed from the event. Photons are removed if they reside within \(\Delta R=0.4\) of a baseline electron, and any jet within \(\Delta R=0.4\) of any remaining photon is discarded.
The \(E_{\text {T}}^{\text {miss}}\) is defined as the magnitude of the negative vector sum, \({\varvec{p}}_{\mathrm {T}}^\mathrm {miss}\), of the transverse momenta of all baseline electrons, muons, jets, and photons [88, 89]. Low-momentum contributions from particle tracks from the primary vertex that are not associated with reconstructed analysis objects are included in the calculation of \(E_{\text {T}}^{\text {miss}}\). This contribution to the \(E_{\text {T}}^{\text {miss}}\) is referred to as the “soft term”.
Models with large hadronic activity are targeted by placing additional requirements on the quantity \(H_\text {T}\), defined as the scalar sum of the \(p_{\text {T}} \) values of all signal jets, or on \(H_{\mathrm {T}}^{\mathrm {incl}}\), the scalar sum of the \(p_{\text {T}} \) values of all signal jets and the two leptons with largest \(p_{\text {T}} \).
All MC samples have correction factors applied to take into account small differences between data and MC simulation in identification, reconstruction and trigger efficiencies for leptons. The \(p_{\text {T}} \) values of leptons in MC samples are additionally smeared to match the momentum resolution in data.
6 Event selection
For each search channel, signal regions (SRs) are designed to target events from specific SUSY signal models. Control regions (CRs) are defined to be depleted in SUSY signal events and enriched in specific SM backgrounds, and they are used to assist in estimating these backgrounds in the SRs. To validate the background estimation procedures, various validation regions (VRs) are defined to be analogous to the CRs and SRs, but with less stringent requirements than the SRs on \(E_{\text {T}}^{\text {miss}}\), \(H_{\mathrm {T}}^{\mathrm {incl}}\) or \(H_{\text {T}}\). Other VRs with additional requirements on the number of leptons are used to validate the modelling of backgrounds in which more than two leptons are expected.
Events in SRs are required to contain at least two signal leptons (electrons or muons). If more than two signal leptons are present in a given event, the selection process continues based on the two leptons with the highest \(p_{\text {T}}\) values in the event.
The selected events must pass at least one of the leptonic triggers. If an event is selected by a dilepton trigger, the two leading, highest \(p_{\text {T}}\), leptons must be matched to one of the objects that triggered the event. These leptons must also have \(p_{\text {T}}\) higher than the threshold of the trigger in question. For events selected by a single-lepton trigger, at least one of the two leading leptons must be matched to the trigger object in the same way. The leading two leptons in the event must have \(p_{\text {T}} >25\) \(\text {GeV}\), and form an SFOS pair.
As at least two jets are expected in all signal models studied, selected events are further required to contain at least two signal jets. Furthermore, events in which the azimuthal opening angle between either of the leading two jets and the \(E_{\text {T}}^{\text {miss}}\) satisfies \(\Delta \phi (\text {jet}_{12},{\varvec{p}}_{\mathrm {T}}^\mathrm {miss})<0.4\) are rejected so as to remove events with \(E_{\text {T}}^{\text {miss}}\) from jet mismeasurements. This requirement also suppresses \(t\bar{t}\) events in which the top quark, the anti-top quark, or the entire \(t\bar{t}\) system has a large Lorentz boost.
The various methods used predict the background in the SRs are discussed in Sect. 7. The selection criteria for the CRs, VRs, and SRs are summarised in Tables 3 and 4. The most important of these regions are shown graphically in Fig. 2.
For the on-shell Z search, the leading-lepton \(p_{\text {T}}\) threshold is raised to 50 \(\text {GeV}\) to increase the sensitivity to signal models with final-state Z bosons. This is an increased leading-lepton \(p_{\text {T}}\) threshold relative to Ref. [20] and is found to better reject fake-lepton candidates from misidentified jets, photon conversions and b-hadron decays, while retaining high efficiency for signal events, which tend to produce boosted Z bosons. To select events containing a leptonically decaying Z boson, the invariant mass of the dilepton system is required to be \(81<m_{\ell \ell }<101\) \(\text {GeV}\). In the CRs and VRs that use the Z mass sidebands, only events with \(m_{\ell \ell }>45\) \(\text {GeV}\) are used to reject the lower \(m_{\ell \ell }\) region dominated by Drell–Yan (DY) production. In Ref. [20] an “on-Z” SR, denoted SRZ, is defined requiring \(E_{\text {T}}^{\text {miss}} >225\) \(\text {GeV}\) and \(H_{\mathrm {T}}^{\mathrm {incl}}> 600\) \(\text {GeV}\). The region is motivated by SUSY signals with high gluino or squark mass and high jet activity. Since b-jets are not always expected in the simplified models considered here, no requirement is placed on b-tagged jet multiplicity (\(n_{b\mathrm {-jets}}\)) so as to be as inclusive as possible and to be consistent with Ref. [20]. Dedicated CRs are defined, with selection criteria similar to those of SRZ, to estimate the contribution from the dominant SM backgrounds in SRZ. These CRs are discussed in more detail in Sect. 7.
The edge selection requires at least two leptons with \(p_{\text {T}} >25\) \(\text {GeV}\). The search is performed across the full \(m_{\ell \ell }\) spectrum, with the exception of the region with \(m_{\ell \ell }<12\) GeV, which is vetoed to reject low-mass DY events and the \(J/\psi \) and \(\Upsilon \) resonances. Three regions are defined to target signal models with low, medium and high values of \(\Delta m_{\tilde{g}} = m(\tilde{g}) - m(\tilde{\chi }^0_1)\), denoted SR-low, SR-medium, and SR-high, respectively. All these regions require \(E_{\text {T}}^{\text {miss}} >200\) \(\text {GeV}\). SR-medium and SR-high also include the requirements \(H_{\text {T}} >400\) \(\text {GeV}\) and \(H_{\text {T}} >700\) \(\text {GeV}\), respectively, to further isolate high-\(\Delta m_{\tilde{g}}\) events. Here the leptons are not included in the \(H_{\text {T}}\) definition to avoid introducing any bias in the \(m_{\ell \ell }\) distribution. Events selected in SR-low, SR-medium and SR-high are further grouped into non-orthogonal \(m_{\ell \ell }\) windows, which represent the search regions used in the edge analysis. The dilepton mass ranges of these are chosen to maximise sensitivity to the targeted signal models, with the window boundaries being motivated by the dilepton mass endpoints of generated signal points. In total, 24 \(m_{\ell \ell }\) windows are defined by selecting ranges with the best expected sensitivity to signal models. Of these windows, nine are in SR-low, eight are in SR-medium and seven are in SR-high. Details of the \(m_{\ell \ell }\) definitions in these regions are given along with the results in Sect. 9. Models without light sleptons are targeted by windows with \(m_{\ell \ell }<60\) \(\text {GeV}\) or \(m_{\ell \ell }<80\) \(\text {GeV}\) for \(\Delta m_\chi < m_Z\) leading to off-shell Z bosons, and by the window with \(81<m_{\ell \ell }<101\) \(\text {GeV}\) for \(\Delta m_\chi > m_Z\) leading to on-shell Z bosons. Models with light sleptons are targeted by the remaining \(m_{\ell \ell }\) windows, which cover the full \(m_{\ell \ell }\) range. The edge selection and on-shell Z selection are not orthogonal. In particular, SR-medium in the range \(81<m_{\ell \ell }<101\) \(\text {GeV}\) overlaps significantly with SRZ.
For the combined \(ee+\mu \mu \) channels, the typical signal acceptance times efficiency values for the signal models considered in SRZ are 2–8%. They are 8–40%, 3–35%, and 1–35%, inclusively in \(m_{\ell \ell }\), for SR-low, SR-medium and SR-high, respectively. The on-shell Z and edge analyses are each optimised for different signal models. There are models in which signal contamination in CRs or VRs can become significant. For example, CRT in Table 3 is used to normalise the \(t\bar{t}\) MC sample to data as a cross-check in the on-shell Z search, but it is a region where the signal contamination from signal models targeted by the edge search can be up to 80% relative to the expected background. In addition, the contamination from on-shell Z signals in the region used to validate the \(Z/\gamma ^{*}+\text {jets}\) and flavour-symmetric estimates, VR-S, is up to 60% for models with \(m(\tilde{g})<1\) \(\text {TeV}\). The signal contamination from the slepton models in the DF regions used to estimate the flavour-symmetric backgrounds in the edge search, CR-FS-low/medium/high in Table 4, is less than 20% for models with \(m(\tilde{g})>600\) \(\text {GeV}\). It is only the contamination in these \(e\mu \) CRs that is relevant in terms of the model-dependent interpretation of the results, and its impact is further discussed in Sect. 10. In general, for models giving substantial contamination in the CRs, the signal-to-background ratio in the SRs is found to be large enough for this contamination to have negligible impact on the sensitivity of the search.
7 Background estimation
The dominant background processes in the SRs are “flavour-symmetric” (FS) backgrounds, where the ratio of ee, \(\mu \mu \) and \(e\mu \) dileptonic branching fractions is 1:1:2 because the two leptons originate from independent \(W\rightarrow \ell \nu \) decays. This background is dominated by \(t\bar{t}\) (50–70%) and also includes WW, Wt, and \(Z\rightarrow \tau \tau \) processes. The FS background constitutes 60–90% of the expected background in the SRs, and is estimated using control samples of \(e\mu \) events.
As all the SRs have a high-\(E_{\text {T}}^{\text {miss}}\) requirement, \(Z/\gamma ^{*}+\text {jets}\) events only enter the SRs when there is large \(E_{\text {T}}^{\text {miss}}\) originating from instrumental effects or from neutrinos in jet fragments. This background is generally small, but it is difficult to model with MC simulation and can mimic signal, particularly for the on-shell Z search. This background is estimated using a control sample of \(\gamma +\text {jets}\) events in data, which are kinematically similar to \(Z/\gamma ^{*}+\text {jets}\) and have similar sources of \(E_{\text {T}}^{\text {miss}}\).
The production of WZ / ZZ dibosons contributes approximately 30% of the SM background in SRZ and up to 20% of the background in the edge SR \(m_{\ell \ell }\) windows. These backgrounds are estimated from MC simulation, after validation in dedicated \(3\ell \) (WZ) and \(4\ell \) (ZZ) VRs. Rare top backgrounds, which include \(t\bar{t}W\), \(t\bar{t}Z\) and \(t\bar{t}WW\) processes, constitute <5% of the expected SM background in all SRs, and are estimated from MC simulation. The contribution from events with fake or misidentified leptons is at most 15% (in one of the edge \(m_{\ell \ell }\) ranges in SR-low), but is generally <5% of the expected SM background in most SRs.
7.1 Flavour-symmetric backgrounds
The flavour-symmetric background is dominant in all SRs. To estimate the contribution of this background to each SR, the so-called “flavour-symmetry” method, detailed in Ref. [20], is used. In this method, data events from a DF control sample, which is defined with the same kinematic requirements as the SR, are used to determine the expected event yields in the SF channels. In the on-shell Z analysis, the method is used to predict the background yield in the Z mass window, defined as \(81<m_{\ell \ell }<101\) \(\text {GeV}\). In the edge analysis, the method is extended to predict the full dilepton mass shape, such that a prediction can be extracted in any of the predefined \(m_{\ell \ell }\) windows.
For the edge search, the flavour-symmetric contribution to each \(m_{\ell \ell }\) bin of the signal regions is predicted using data from the corresponding bin in an \(e\mu \) control region. All edge CR-FS regions (definitions can be seen in Table 4) are 88–97% pure in flavour-symmetric processes (this purity is calculated from MC simulation).
For the on-shell search, this method is complicated slightly by a widening of the \(m_{\ell \ell }\) window used in CR-FS, the \(e\mu \) control region (defined in Table 3). The window is enlarged to \(61<m_{\ell \ell }<121\) \(\text {GeV}\) to approximately triple the amount of data in the control region and thus increase the statistical precision of the method. This results in a region that is \(\sim \)95% pure in flavour-symmetric processes (the expected composition of this 95% is \(\sim \)80% \(t\bar{t}\), \(\sim \)10% Wt, \(\sim \)10% WW and <1% \(Z \rightarrow \tau \tau \)).
Apart from the \(m_{\ell \ell }\) widening in CR-FS, the method used is identical for the on-shell and edge regions. Events in the control regions are subject to lepton \(p_{\text {T}} \)- and \(\eta \)-dependent correction factors measured in data and MC simulation. Because the triggers used are not identical in 2015 and 2016, these factors are measured separately for each year and account for the different identification and reconstruction efficiencies for electrons and muons, as well as the different trigger efficiencies for the dielectron, dimuon and electron–muon selections. The estimated numbers of events in the SF channels, \(N^\text {est}_{ee/\mu \mu }\), are given by:
where \(N_{e\mu }^\text {data}\) is the number of data events observed in a given control region, \(\alpha (p_{\text {T}} ^i, \eta ^i)\) accounts for the different trigger efficiencies for SF and DF events, and \(k_{e}(p_{\text {T}} ^{i,\mu }, \eta ^{i,\mu })\) and \(k_{\mu }(p_{\text {T}} ^{i,e}, \eta ^{i,e})\) are electron and muon selection efficiency factors for the kinematics of the lepton being replaced, in event i. The trigger and selection efficiency correction factors are derived from the events in an inclusive on-Z selection (\(81<m_{\ell \ell }<101\) \(\text {GeV}\), \(\ge 2\) jets), according to:
where \(\epsilon ^\text {trig}_{ee/\mu \mu }\) is the trigger efficiency and \(N_{ee/\mu \mu }^{\text {meas}}\) is the number of \(ee/\mu \mu \) events in the inclusive on-Z region outlined above. Here \(k_{e}(p_{\text {T}}, \eta )\) and \(k_{\mu }(p_{\text {T}}, \eta )\) are calculated separately for leading and sub-leading leptons, while \(\alpha \) is calculated for the leading lepton, \(\ell _1\). The correction factors are typically within 10% of unity, except in the region \(|\eta |<0.1\) where, because of the lack of coverage by the muon spectrometer, they are up to 50% from unity. For all background estimates based on the flavour-symmetry method, results are computed separately for ee and \(\mu \mu \) and then summed to obtain the combined predictions. The resulting estimates from the DF channels are scaled according to the fraction of flavour-symmetric backgrounds in each \(e\mu \) control sample, \(f_{\mathrm {FS}}\) (\(95\%\) in CR-FS), which is determined by subtracting non-flavour-symmetric backgrounds taken from MC simulation from the data observed in the corresponding \(e\mu \) region. In the on-shell case, the result is also scaled by the fraction of events in CR-FS expected to be contained within \(81<m_{\ell \ell }<101\) \(\text {GeV}\), \(f_{Z \text {-mass}}\) (\(38\%\)), which is otherwise set to 100% for the edge regions. The validity of extrapolating in \(m_{\ell \ell }\) between CR-FS and SRZ was checked by comparing the \(m_{\ell \ell }\) shapes in data and MC simulation in a region similar to VR-S, but with the \(m_{\ell \ell }\) requirement relaxed and \(H_{\mathrm {T}}^{\mathrm {incl}}\) \( > 300\) \(\text {GeV}\) to obtain a sample with a large number of events. The resulting on-Z fractions in MC simulation were found to agree with data within statistical uncertainties, which are summed in quadrature to assign a systematic uncertainty. In the case of the edge search the full \(m_{\ell \ell }\) distribution is validated by applying a flavour-symmetry method to \(t\bar{t}\) MC evnets in VR-low, VR-medium and VR-high. This procedure results in good closure, which is further discussed in Sect. 7.5. The difference between the prediction and the observed distribution is used to assign an MC non-closure uncertainty to the estimate.
The flavour-symmetry method in SRZ is further cross-checked by performing a profile likelihood fit [90] of MC yields to data in the Z-mass sidebands (\(m_{\ell \ell }\notin [81, 101]\) \(\text {GeV}\)), the region denoted CRT in Table 3, which is dominated by \(t\bar{t}\) (with a purity of >75%) and contains 273 events in data. The other flavour-symmetric processes in this region contribute \(\sim \)12% (Wt), 10% (WW) and <1% (\(Z\rightarrow \tau \tau \)). All SM background processes are taken directly from MC simulation in this cross-check, including backgrounds also estimated using the flavour-symmetry method. The normalisation of the dominant \(t\bar{t}\) background is a free parameter and is the only parameter affected by the fit. For this cross-check, the contamination from Beyond Standard Model processes in the Z-mass sidebands is assumed to be negligible. The fit results in a scale factor of 0.64 for the \(t\bar{t}\) yield predicted by simulation. This result is extrapolated from the Z-mass sidebands to SRZ and gives a prediction of \( 29 \pm 7 \) events, which is consistent with the nominal flavour-symmetry background estimate of \(33 \pm 4\) in this region.
The sideband fit is repeated at lower \(E_{\text {T}}^{\text {miss}}\) in VRT, with the results being propagated to VR-S, so as to test the \(m_{\ell \ell }\) extrapolation used in the sideband fit method. The normalisation to data in this region, which is at lower \(E_{\text {T}}^{\text {miss}}\) relative to CRT, results in a scale factor of 0.80 for the \(t\bar{t}\) yield predicted by simulation. The number of FS events predicted in VR-S using the sideband fit in VRT is compatible with the number estimated by applying the FS method to data in VR-FS. The results of the background estimate in both VR-S and SRZ obtained from the flavour-symmetry method are compared with the values obtained by the sideband fit cross-check in Table 5. The methods result in consistent estimates in both regions. Further results in the edge VRs are discussed in Sect. 7.5.
A potential cause of the low scale factors obtained from the sideband fit at large \(H_{\text {T}}\) and \(E_{\text {T}}^{\text {miss}}\) is mismodelling of the top-quark \(p_{\text {T}}\) distribution, where measurements of \(t\bar{t}\) differential cross sections by the ATLAS and CMS experiments indicate that the top-quark \(p_{\text {T}}\) distribution predicted by most generators is harder than that observed in data [91, 92]. Corrections to the MC predictions according to NNLO calculations provided in Ref. [93] indicate an improvement in the top-quark pair modelling at high \(H_{\text {T}}\), which should lead to scale factors closer to unity. With the data-driven method used to estimate \(t\bar{t}\) contributions in this analysis, the results do not depend on these corrections. They are therefore not applied to the \(t\bar{t}\) MC sample for the sideband-fit cross-check.
7.2 \(Z/\gamma ^{*}+\text {jets}\) background
The \(Z/\gamma ^{*}+\text {jets}\) background estimate is based on a data-driven method that uses \(\gamma +\text {jets}\) events in data to model the \(E_{\text {T}}^{\text {miss}}\) distribution of \(Z/\gamma ^{*}+\text {jets}\). The \(\gamma +\text {jets}\) and \(Z/\gamma ^{*}+\text {jets}\) processes have similar event topologies, with a well-measured object recoiling against a hadronic system, and both tend to have \(E_{\text {T}}^{\text {miss}}\) that stems from jet mismeasurements and neutrinos in hadronic decays. In this method, which has been used by CMS in a search in this final state [18], a sample of data events containing at least one photon and no leptons is constructed using the same kinematic selection as each of the SRs, without the \(E_{\text {T}}^{\text {miss}}\) and \(\Delta \phi (\text {jet}_{12},{\varvec{p}}_{\mathrm {T}}^\mathrm {miss})\) requirements (the CR\(\gamma \) regions defined in Tables 3, 4).
The requirement \(\Delta \phi (\text {jet}_{12},{\varvec{p}}_{\mathrm {T}}^\mathrm {miss}) > 0.4\) applied in the SRs suppresses \(E_{\text {T}}^{\text {miss}}\) from jet mismeasurements and increases the relative contributions to \(E_{\text {T}}^{\text {miss}}\) from the photon, electrons, and muons. The difference in resolution between photons, electrons, and muons can be significant at high \(p_{\text {T}} \). Therefore, before the \(\Delta \phi (\text {jet}_{12},{\varvec{p}}_{\mathrm {T}}^\mathrm {miss}) > 0.4\) requirement is applied, the photon \(p_{\text {T}}\) is smeared according to a \(Z\rightarrow ee\) or \(Z\rightarrow \mu \mu \) resolution function. The smearing function is derived by comparing the \(E_{\text {T}}^{\text {miss}}\)-projection along the boson momentum in \(Z/\gamma ^{*}+\text {jets}\) and \(\gamma +\text {jets}\) MC events in a 1-jet control region with no other event-level kinematic requirements. A deconvolution is applied to avoid including the photon resolution in the Z resolution. For each event, a photon \(p_{\text {T}}\) smearing \(\Delta p_{\text {T}} \) is obtained by sampling the smearing function. The photon \(p_{\text {T}}\) is shifted by \(\Delta p_{\text {T}} \), with the parallel component of the \(E_{\text {T}}^{\text {miss}}\) being correspondingly adjusted by \(-\Delta p_{\text {T}} \).
The smeared \(\gamma +\text {jets}\) events are then reweighted to match the boson \(p_{\text {T}} \) distribution of the \(Z/\gamma ^{*}+\text {jets}\) events. This reweighting is applied separately in each region and accounts for small differences between the \(\gamma +\text {jets}\) events and \(Z/\gamma ^{*}+\text {jets}\) events, which arise mainly from the mass of the Z boson. The reweighting is done using \(Z/\gamma ^{*}+\text {jets}\) events in data, and is checked using \(Z/\gamma ^{*}+\text {jets}\) MC simulation in an MC closure test, as described further below. Following this smearing and reweighting procedure, the \(E_{\text {T}}^{\text {miss}}\) of each \(\gamma +\text {jets}\) event is recalculated, and the final \(E_{\text {T}}^{\text {miss}}\) distribution is obtained after applying the \(\Delta \phi (\text {jet}_{12},{\varvec{p}}_{\mathrm {T}}^\mathrm {miss} ) > 0.4\) requirement. For each SR, the resulting \(E_{\text {T}}^{\text {miss}}\) distribution is normalised to data in a CRZ with the same requirements except that the SR \(E_{\text {T}}^{\text {miss}}\) requirement is replaced by \(E_{\text {T}}^{\text {miss}} <60\) \(\text {GeV}\).
The shape of the \(Z/\gamma ^{*}+\text {jets}\) \(m_{\ell \ell }\) distribution is extracted from MC simulation and validated by comparing to data in events with lower \(E_{\text {T}}^{\text {miss}}\) requirements and a veto on b-tagged jets, to suppress the background from \(t\bar{t}\). The \(m_{\ell \ell }\) distribution is modelled by parameterising the \(m_{\ell \ell }\) in \(Z/\gamma ^{*}+\text {jets}\) events as a function of the difference between reconstructed and true Z boson \(p_{\text {T}}\) in MC simulation. This parameterization ensures that the correlation between lepton momentum mismeasurement and observed \(m_{\ell \ell }\) values far from the Z boson mass is preserved. Each photon event is assigned an \(m_{\ell \ell }\) via a random sampling of the corresponding distribution, equating photon \(\Delta p_{\text {T}} \) and the difference between true and reconstructed Z boson \(p_{\text {T}}\). The resulting \(m_{\ell \ell }\) distribution in \(\gamma +\text {jets}\) MC simulation is compared to that extracted from \(Z/\gamma ^{*}+\text {jets}\) MC simulation and the difference is assessed as a systematic uncertainty in the background prediction for each \(m_{\ell \ell }\) bin.
The full smearing, reweighting, and \(m_{\ell \ell }\) assignment procedure is applied to the \(V\gamma \) MC sample in parallel with the \(\gamma +\text {jets}\) data sample. After applying all corrections to both samples, the \(V\gamma \) contribution to the \(\gamma +\text {jets}\) data sample is subtracted to remove contamination from backgrounds with real \(E_{\text {T}}^{\text {miss}}\). Contamination by events with fake photons in these \(\gamma +\text {jets}\) data samples is small, and this contribution is therefore neglected.
In the \(H_{\text {T}}\)-inclusive region corresponding to VR-low, there is a non-negligible contribution expected from \(Z/\gamma ^{*}+\text {jets}\) events with \(p_{\text {T}} ^{Z}<37\) \(\text {GeV}\). Given the photon trigger strategy discussed in Sect. 4, no photons with \(p_{\text {T}} <37\) \(\text {GeV}\) are included in the event selection. To account for this photon \(p_{\text {T}}\) threshold, a boson-\(p_{\text {T}}\) correction of up to 50% is applied as a function of \(E_{\text {T}}^{\text {miss}}\) in VR-low. This correction uses the fraction of \(Z/\gamma ^{*}+\text {jets}\) events in a given \(E_{\text {T}}^{\text {miss}}\) bin expected to have \(p_{\text {T}} ^{Z}<37\) \(\text {GeV}\), according to MC simulation. The \(\gamma +\text {jets}\) data are then scaled according to this fraction, as a function of \(E_{\text {T}}^{\text {miss}}\), to correct for the missing \(p_{\text {T}} ^{Z}<37\) \(\text {GeV}\) contribution. The correction is found to be negligible in all signal regions.
The distribution of \(E_{\text {T}}^{\text {miss}}\) obtained in Sherpa \(Z/\gamma ^{*}+\text {jets}\) MC simulation is compared to that obtained by applying this background estimation technique to Sherpa \(\gamma +\text {jets}\) MC samples. In this check the \(\gamma +\text {jets}\) MC simulation is reweighted according to the \(p_{\text {T}}\) distribution given by the \(Z/\gamma ^{*}+\text {jets}\) MC simulation. The result of this MC closure check is shown in Fig. 3a for events in VRZ (without an upper \(E_{\text {T}}^{\text {miss}}\) cut), where good agreement between \(Z/\gamma ^{*}+\text {jets}\) and corrected \(\gamma +\text {jets}\) MC simulation can be seen across the entire \(E_{\text {T}}^{\text {miss}}\) spectrum. A comparison between the full \(E_{\text {T}}^{\text {miss}}\) spectrum in data and the \(Z/\gamma ^{*}+\text {jets}\) background estimated via the \(\gamma +\text {jets}\) method is also shown in Fig. 3b for events in VRZ. The systematic uncertainties associated with this method are described in Sect. 8.
7.3 Fake-lepton background
Semileptonic \(t\bar{t}\), \(W\rightarrow \ell \nu \) and single top (s- and t-channel) events enter the dilepton channels via “fake” leptons. These can include misidentified hadrons, converted photons or non-prompt leptons from b-hadron decays. The extent of this background is estimated using the matrix method, detailed in Ref. [94]. Its contribution in regions with high lepton \(p_{\text {T}}\) and dilepton invariant mass is negligible, but in the edge search, where lower-\(p_{\text {T}}\) leptons are selected and events can have low \(m_{\ell \ell }\), the fake-lepton background can make up to 15% of the total background. In this method a control sample is constructed using baseline leptons, thereby enhancing the probability of selecting a fake lepton due to the looser lepton selection and identification criteria relative to the signal lepton selection. For each relevant CR, VR or SR, the region-specific kinematic requirements are placed upon this sample of baseline leptons. The number of events in this sample in which the selected leptons subsequently pass (\(N_{\text {pass}}\)) or fail (\(N_{\text {fail}}\)) the signal lepton requirements in Sect. 5 are then counted. In the case of a one-lepton selection, the number of fake-lepton events in a given region is then estimated according to:
Here \(\epsilon ^{\text {real}}\) is the relative identification efficiency (from baseline to signal) for genuine, prompt (“real”) leptons and \(\epsilon ^{\text {fake}}\) is the relative identification efficiency (again from baseline to signal) with which non-prompt leptons or jets might be misidentified as prompt leptons. This principle is then expanded to a dilepton selection by using a four-by-four matrix to account for the various possible real–fake combinations for the two leading leptons in an event.
The real-lepton efficiency, \(\epsilon ^{\text {real}}\), is measured in \(Z\rightarrow \ell \ell \) data events using a tag-and-probe method in CR-real, defined in Table 4. In this region the \(p_{\text {T}}\) of the leading lepton is required to be \(>\)40 \(\text {GeV}\), and only events with exactly two SFOS leptons are selected. The fake-lepton efficiency, \(\epsilon ^{\text {fake}}\), is measured in CR-fake, a region enriched with fake leptons by requiring same-sign lepton pairs. The lepton \(p_{\text {T}}\) requirements are the same as those in CR-real, with the leading lepton being tagged as the “real” lepton and the fake efficiency being evaluated based on the sub-leading lepton in the event. An \(E_{\text {T}}^{\text {miss}}\) requirement of <125 \(\text {GeV}\) is used to reduce possible contamination from Beyond Standard Model processes. In this region the background due to prompt-lepton production, estimated from MC simulation, is subtracted from the total data contribution. Prompt-lepton production makes up 7% (11%) of the baseline electron (muon) sample and 10% (61%) of the signal electron (muon) sample in CR-fake. From the resulting data sample the fraction of events in which the baseline leptons pass a signal-like selection yields the fake efficiency. Both the real- and fake-lepton efficiencies are binned as a function of lepton \(p_{\text {T}}\) and calculated separately for the 2015 and 2016 data sets.
This method is validated by checking the closure in MC simulation and data–background agreement in VR-fake.
7.4 Diboson and rare top processes
The remaining SM background contribution in the SRs is due to WZ / ZZ diboson production and rare top processes (\(t\bar{t} Z\), \(t\bar{t} W\) and \(t\bar{t} WW\)). The rare top processes compose <5% of the expected SM background in the SRs and are taken directly from MC simulation.
Production of WZ / ZZ dibosons constitutes about 30% of the expected background in SRZ and up to 20% in some edge SR \(m_{\ell \ell }\) windows. In SRZ, this background is composed of roughly 70% WZ, about 40% of which is \(WZ\rightarrow \ell \ell \tau \nu \). This is the largest background contribution that is estimated from MC simulation, and must be carefully validated, especially because these backgrounds contain Z bosons and can thus mimic a signal by producing a peak at \(m_{\ell \ell }\approx m_Z\). To validate the MC modelling of these backgrounds, VRs with three leptons (VR-WZ) and four leptons (VR-ZZ) are defined (selection shown in Table 3). In VR-WZ, from the three selected leptons in an event, the SFOS pair with \(m_{\ell \ell }\) most consistent with the Z mass is indentified as the Z candidate. The transverse mass of the remaining lepton and the \(E_{\text {T}}^{\text {miss}}\), \(m_{\text {T}}(\ell _{3},E_{\text {T}}^{\text {miss}})\), is then required to be <100 \(\text {GeV}\), forming the W candidate. In VR-ZZ an \(E_{\text {T}}^{\text {miss}} <100\) \(\text {GeV}\) requirement is used to suppress WZ and top processes. The yields and kinematic distributions observed in these regions are well-modelled by MC simulation. In particular, the \(E_{\text {T}}^{\text {miss}}\), \(H_\text {T}\), jet multiplicity, and boson \(p_{\text {T}} \) distributions show good agreement. An additional three-lepton VR (VR-3L) is defined to provide validation of the diboson background in a region of phase space closer to the SR; good agreement is observed in this region as well.
7.5 Results in validation regions
The expected background yields in VR-S are shown in Table 6 and compared with the observed data yield. Agreement between the data and the expected Standard Model background is observed. The expected background yields in the three diboson VRs are also shown in Table 6. The data are consistent with the expected background. Similar information for the edge VRs is provided in Table 7. Data and background estimates are in agreement within uncertainties.
Figure 4 shows the observed and expected \(m_{\ell \ell }\) distributions in the same edge VRs. The same background estimation methods are applied to both MC simulation and data. In the MC studies, the flavour-symmetry method of Sect. 7.1 is applied to \(t\bar{t}\) MC simulation, and the observed SF \(m_{\ell \ell }\) distribution is compared to the prediction based on DF events. In the data studies, the observed SF \(m_{\ell \ell }\) distribution is compared to the sum of FS backgrounds from the extended flavour-symmetry method, the \(Z/\gamma ^{*}+\text {jets}\) background from the \(\gamma +\text {jets}\) method, and the WZ / ZZ diboson, rare top, and fake-lepton backgrounds.
The observed MC closure is good in all validation regions. The data agree with the expected background in the validation regions as well. No significant discrepancies or trends are apparent.
8 Systematic uncertainties
The data-driven background estimates are subject to uncertainties associated with the methods employed and the limited number of events used in their estimation. The dominant uncertainty (10%) for the flavour-symmetry-based background estimate in SRZ is due to the limited number of events in CR-FS. Other systematic uncertainties assigned to this background estimate include those due to MC closure (3%), the measurement of the efficiency correction factors (3%) and the extrapolation in \(m_{\ell \ell }\) (1%). In the case of the edge SRs the statistical uncertainty is also the dominant uncertainty in the flavour-symmetric background estimate in the case of SR-high, but for both SR-medium and SR-low the uncertainties from the MC non-closure and efficiency correction factors are comparable in size, or in some cases larger. These uncertainties can contribute up to 5% in SR-low and SR-medium and 10% in SR-high.
Several sources of systematic uncertainty are assessed for the \(Z/\gamma ^{*}+\text {jets}\) background. The boson \(p_{\text {T}} \) reweighting procedure is assigned an uncertainty based on a comparison of the nominal results with those obtained by reweighting using three other kinematic variables, namely \(H_{\text {T}}\), Z-boson \(E_{\text {T}}\) and jet multiplicity. For the smearing function, which is measured using MC events in a 1-jet control region, an uncertainty is derived by comparing the results obtained using the nominal smearing function with those obtained using a smearing function from a 2-jet sample of MC events, and also using a smearing function measured in a 1-jet data sample. An uncertainty of between 40–100% is assigned to account for different reweighting procedures and between 20–100% for the smearing procedure applied to \(\gamma +\text {jets}\) events. The smearing uncertainty dominates in SR-high, while the reweighting uncertainty dominates in SR-low and SR-medium, with both being around 60% in SRZ. The full reweighting and smearing procedure is carried out using \(\gamma +\text {jets}\) MC events such that an MC non-closure uncertainty can be derived by comparing the resulting \(\gamma +\text {jets}\) MC \(E_{\text {T}}^{\text {miss}}\) distribution to that in \(Z/\gamma ^{*}+\text {jets}\) MC events. The resulting uncertainty of up to 35% is calculated in the VRs, so as to maximise the number of events that contribute. An uncertainty of 16% is assessed for the \(V\gamma \) backgrounds, based on data-to-MC agreement in a \(V\gamma \)-enriched control region. This uncertainty is propagated to the final \(Z/\gamma ^{*}+\text {jets}\) estimate following the subtraction of the \(V\gamma \) background. In VR-low, a correction is applied to the \(E_{\text {T}}^{\text {miss}}\) distribution in \(\gamma +\text {jets}\) events to account for the fraction of \(Z/\gamma ^{*}+\text {jets}\) events in this \(H_{\text {T}}\)-inclusive region expected to have boson \(p_{\text {T}}\) less than 37 \(\text {GeV}\). The full size of this correction (up to 50% for \(E_{\text {T}}^{\text {miss}} =150\) \(\text {GeV}\)) is applied as a systematic uncertainty. The \(m_{\ell \ell }\) distribution assigned to \(\gamma +\text {jets}\) MC events is compared to that of \(Z/\gamma ^{*}+\text {jets}\) MC events, and the relative difference in a given \(m_{\ell \ell }\) bin is assigned as an uncertainty. Finally, the statistical precision of the estimate also enters as a systematic uncertainty of \(\sim \)10% in the final background estimate. After applying the correction procedure, differences in the number of b-tagged jets between \(Z/\gamma ^{*}+\text {jets}\) and \(\gamma +\text {jets}\) are found to be negligible, indicating good agreement in heavy-flavour content.
The uncertainties in the fake-lepton background stem from the number of events in the regions used to measure the real- and fake-lepton efficiencies, the limited size of the inclusive loose-lepton sample, and from varying the region used to measure the fake-lepton efficiency. The nominal fake-lepton efficiency is compared with those measured in a region with b-tagged jets and a region with a b-jet, as well as a region with the prompt-lepton subtraction varied by \(20\%\). Varying the sample composition via b-jet tagging gives the largest uncertainty. The uncertainty for the edge SRs from the statistical component of the lepton efficiencies is 30–45%, and from varying the region for the fake-lepton efficiency it is 50–75%. The uncertainties in SRZ are generally larger due to the small number of events contributing to the estimate in this region.
Theoretical and experimental uncertainties are taken into account for the signal models, as well as background processes that rely on MC simulation. The estimated uncertainty in the luminosity measurement is \(2.9\%\) [31, 32]. The jet energy scale is subject to uncertainties associated with the jet flavour composition, the pile-up and the jet and event kinematics [81]. Uncertainties in the jet energy resolution are included to account for differences between data and MC simulation [81]. An uncertainty in the \(E_{\text {T}}^{\text {miss}}\) soft-term resolution and scale is taken into account [88], and uncertainties due to the lepton energy scales and resolutions, as well as trigger, reconstruction, and identification efficiencies, are also considered.
The WZ / ZZ processes are assigned a cross-section uncertainty of \(6\%\) and an additional uncertainty based on comparisons between Sherpa and Powheg MC samples, which is up to \(50\%\) in the SRs. Uncertainties due to the choice of factorisation and renormalisation scales are calculated by varying the nominal values up and down by a factor of two and can be up to \(23\%\). For rare top processes, a \(13\%\) PDF and scale variation uncertainty is applied [34] in addition to a \(22\%\) cross-section uncertainty [61,62,63].
For signal models, the nominal cross section and the uncertainty are taken from an envelope of cross-section predictions using different PDF sets and factorisation and renormalisation scales, as described in Refs. [95, 96]. These are calculated at next-to-leading-logarithm accuracy (NLO + NLL) [51,52,53,54,55], and the resulting uncertainties range from 16 to 30%.
A breakdown of the dominant uncertainties in the background prediction in the SRs is provided in Table 8 for the on-shell Z and edge searches. Here these uncertainties are quoted relative to the total background. In the case of the edge regions a range is quoted, taking into account the relative contribution of the given uncertainty in each of the \(m_{\ell \ell }\) ranges in SR-low, SR-medium and SR-high. The largest uncertainties in the signal regions are due to the size of the \(e\mu \) data sample in CR-FS, used to provide the flavour-symmetric background estimate, the combined systematic uncertainty in the same background, the systematic uncertainty in \(\gamma +\text {jets}\), or, in the case of SRZ, the WZ / ZZ generator uncertainty. The statistical component of the uncertainty from the flavour-symmetry estimate is largest for the edge analysis in SR-medium and SR-high in the highest \(m_{\ell \ell }\) regions. In the edge SRs the uncertainty in the WZ / ZZ background tends to be highest in the \(m_{\ell \ell }\) ranges that include the Z window. The uncertainty in the fake-lepton background is largest in SR-high, where fake leptons can compose a larger fraction of the background. Experimental uncertainties have a far lower impact on the systematic uncertainty of the total background (<2%).
9 Results
9.1 Results in SRZ
For the on-shell Z search, the expected background and observed yields in the SR are shown in Table 9. A total of 60 events are observed in data with a predicted background of \(53.5 \pm 9.3\) events. There are 35 events observed in data in the ee channel, and 25 events observed in the \(\mu \mu \) channel. The probability for the background to produce a fluctuation greater than or equal to that observed in the data, called the significance when expressed in terms of the number of standard deviations, corresponds to \(0.47\sigma \) (details of the significance calculation are presented in Sect. 10). The level of agreement between the observed event yields in data and the background predictions in the VRs, shown previously in Table 6, is also displayed in Fig. 5, along with the results in SRZ.
The dilepton invariant-mass distribution for the \(ee+\mu \mu \) and \(e\mu \) channels with the kinematic requirements of SRZ, but over the full \(m_{\ell \ell }\) range, is shown in Fig. 6. Here the data are consistent with the expected background over the full \(m_{\ell \ell }\) range. The dilepton invariant-mass, jet and b-tagged jet multiplicity, \(E_{\text {T}}^{\text {miss}}\), \(H_{\mathrm {T}}^{\mathrm {incl}}\) and \(p_{\text {T}} ^{\ell \ell }\) distributions in SRZ are shown in Fig. 7. The shapes of the background distributions in these figures are obtained from MC simulation, where the MC simulation is normalised according to the data-driven estimates in the SR. Here two representative examples of \(\tilde{g}\)–\(\tilde{\chi }_2^0 \) on-shell signal models, with \((m(\tilde{g}),m(\tilde{\chi }^{0}_{2}))=(1095, 205)\) \(\text {GeV}\) and \((m(\tilde{g}),m(\tilde{\chi }^{0}_{2}))=(1240, 960)\) \(\text {GeV}\), are overlaid. To demonstrate the modelling of the \(Z/\gamma ^{*}+\text {jets}\) background in VR-S and SRZ, Fig. 8 shows the minimum \(\Delta \phi (\text {jet}_{12},{\varvec{p}}_{\mathrm {T}}^\mathrm {miss})\) distribution over the full range, where \(\Delta \phi (\text {jet}_{12},{\varvec{p}}_{\mathrm {T}}^\mathrm {miss})>0.4\) is required in VR-S and SRZ. Here the \(Z/\gamma ^{*}+\text {jets}\) distribution is modelled using the full data-driven prediction from \(\gamma +\text {jets}\). Two of the events in the SR contain a third signal lepton.
9.2 Results in the edge SRs
The integrated yields in the edge signal regions are compared to the expected background in Table 10. To allow for the visualisation of a potential edge, the full \(m_{\ell \ell }\) distributions in the three search regions are compared to the expected background in Fig. 9. In addition, the observed \(m_{\ell \ell }\) distributions are compared to the predictions from MC simulation in Fig. 10, in which the \(t\bar{t}\) background is scaled such that the total MC expected yield matches the data in the \(e\mu \) CR. The \(t\bar{t}\) normalisation factors are \(\mu _{t\bar{t}}=0.85\pm 0.03\), \(0.75\pm 0.04\), and \(0.57\pm 0.07\) in SR-low, SR-medium, and SR-high, respectively, where the uncertainty is the data statistical uncertainty. The data-driven flavour-symmetry prediction is used for the quantitative results of the analysis. This prediction does not rely on the \(t\bar{t}\) normalisation scale factors discussed above. The MC-based cross-check method is used to examine the \(m_{\ell \ell }\) distribution in finer bins than can be achieved with the flavour-symmetry method, due to the limited statistical precision of the \(e\mu \) CR.
As signal models may produce kinematic endpoints at any value of \(m_{\ell \ell }\), any excess must be searched for across the \(m_{\ell \ell }\) distribution. To do this a “sliding window” approach is used. The binning in the SRs, shown in Fig. 9, defines many possible dilepton mass windows. The 24 \(m_{\ell \ell }\) ranges (9 for SR-low, 8 for SR-medium, and 7 for SR-high) are chosen because they are the most sensitive for at least one grid point in the signal model parameter space. Some of the ranges overlap. The results in these regions are summarised in Fig. 11, and the expected and observed yields in the combined \(ee+\mu \mu \) channel for all 24 \(m_{\ell \ell }\) ranges are presented in Table 11. In SR-low and SR-medium, the data are consistent with the expected background across the full \(m_{\ell \ell }\) range. In SR-high the data show a slight excess above the background at low \(m_{\ell \ell }\). Of these 24 \(m_{\ell \ell }\) ranges, the largest excess is observed in SR-high with \(12<m_{\ell \ell }<101\) \(\text {GeV}\). Here a total of 90 events are observed in data, compared to an expectation of \(65\pm 10\) events, corresponding to a local significance of \(1.7\sigma \).
10 Interpretation
In this section, exclusion limits are shown for the SUSY models detailed in Sect. 3. The asymptotic \(CL_{\text {S}}\) prescription [90, 98], implemented in the HistFitter program [97], is used to determine cross-section upper limits at \(95\%\) confidence level (CL) for the on-Z search. For the edge search, pseudo-experiments are used to evaluate the cross-section upper limits. A Gaussian model for nuisance parameters is used for all signal and background uncertainties. Exceptions are the statistical uncertainties of the flavour-symmetry method, \(\gamma +\text {jets}\) method and MC-based backgrounds, all of which are treated as Poissonian nuisance parameters. The different experimental uncertainties are treated as correlated between signal and background events. The theoretical uncertainty of the signal cross section is not accounted for in the limit-setting procedure. Instead, following the initial limit determination, the impact of varying the signal cross section within its uncertainty is evaluated separately and indicated in the exclusion results. Limits are based on the combined \(ee+\mu \mu \) results. Possible signal contamination in the CRs is neglected in the limit-setting procedure; the contamination is found to be negligible for signal points near the exclusion boundaries. Far from the exclusion boundary, although the signal contamination can be significant, the number of events appearing in the signal region is large enough that the points are still excluded, due to the relative branching fractions for the signal in the CR and SR. For example, for models with signal contamination of 50% in CR-FS the signal-to-background ratio in SRZ is \(\sim \)10.
The results of the on-shell Z search are interpreted in a simplified model with gluino-pair production, where each gluino decays as \(\tilde{g} \rightarrow q\bar{q} \tilde{\chi }^{0}_{2}, \tilde{\chi }^{0}_{2} \rightarrow Z \tilde{\chi }^{0}_{1}\) and the \(\tilde{\chi }^{0}_{1}\) mass is set to 1 \(\text {GeV}\). The expected and observed exclusion contours for this \(\tilde{g}\)–\(\tilde{\chi }_2^0 \) on-shell grid are shown in the \(m(\tilde{g})\)–\(m(\tilde{\chi }^{0}_{2})\) plane in Fig. 12. The expected (observed) lower limit on the gluino mass is about 1.35 \(\text {TeV}\) (1.30 \(\text {TeV}\)) for a \(\tilde{\chi }_2^0\) with a mass of 1.1 \(\text {TeV}\) in this model. The impact of the systematic uncertainties in the background and the experimental uncertainties in the signal, shown with a coloured band, is about 100 \(\text {GeV}\) on the gluino mass limit. The systematic uncertainty of the signal cross section, shown as dotted lines around the observed contour, has an impact of about 40 \(\text {GeV}\). Figure 12 also shows the expected and observed exclusion limits for the \(\tilde{q}\)–\(\tilde{\chi }_2^0 \) on-shell model. This is a simplified model with squark-pair production, where each squark decays to a quark and a neutralino, with the neutralino subsequently decaying to a Z boson and an LSP with a mass of 1 \(\text {GeV}\). In this model, exclusion is expected (observed) for squarks with masses below 1040 \(\text {GeV}\) (980 \(\text {GeV}\)) for a \(\tilde{\chi }^{0}_{2}\) mass of 600 \(\text {GeV}\).
Figure 13 shows the expected and observed exclusion contours for the \(\tilde{g}\)–\(\tilde{\chi }_1^0 \) on-shell model, in which the produced gluinos follow the same decay chain as in the model above. In this case the mass difference \(\Delta m = m(\tilde{\chi }^{0}_{2})-m(\tilde{\chi }^{0}_{1})\) is set to 100 \(\text {GeV}\).
The results of the edge search are interpreted in two simplified models with gluino-pair production, in which each gluino decays as \(\tilde{g} \rightarrow q\bar{q} \tilde{\chi }^{0}_{2}\). For each point in the signal-model parameter space, limits on the signal strength are calculated using the \(m_{\ell \ell }\) window with the best expected sensitivity. Details of the windows are described in Sect. 9.
The excluded regions in the \(m(\tilde{g})\)–\(m(\tilde{\chi }^{0}_{1})\) plane are presented in Fig. 14 for the slepton model. In this model, pair-produced gluinos each decay as \(\tilde{g}\rightarrow q \bar{q} \tilde{\chi }^{0}_2, \tilde{\chi }^{0}_2 \rightarrow \ell ^\pm \tilde{\ell }^\mp , \tilde{\ell }^\mp \rightarrow \ell ^\mp \tilde{\chi }^{0}_{1}\). Here, the results exclude gluinos with masses as large as 1.7 \(\text {TeV}\), with an expected limit of 1.75 \(\text {TeV}\) for small \(m(\tilde{\chi }^{0}_{1})\). The results probe kinematic endpoints as small as \(m_{\ell \ell }^{\text {max}} = m(\tilde{\chi }^{0}_{2})-m(\tilde{\chi }^{0}_{1}) = 1/2(m(\tilde{g})-m(\tilde{\chi }^{0}_{1}) ) = 50\) \(\text {GeV}\).
The \(Z^{(*)}\) exclusion limits from the results in the edge SRs are compared with the same limits derived using the results in SRZ in Fig. 15. In this model, pair-produced gluinos each decay as \(\tilde{g} \rightarrow q\bar{q} \tilde{\chi }^{0}_{2}, \tilde{\chi }^{0}_{2} \rightarrow Z^{(*)} \tilde{\chi }^{0}_{1}\), and the mass splitting between the \(\tilde{\chi }^{0}_{2}\) and the \(\tilde{\chi }^{0}_{1}\) determines whether the Z boson is produced on-shell. Here the edge limits extend into the more compressed region, whereas the expected SRZ exclusion probes higher \(\tilde{\chi }^{0}_{1}\) masses in the on-shell regime. At high gluino masses, the edge SRs provide stronger limits. For the \(Z^{(*)}\) model, the expected and observed gluino mass limits are 1.4 \(\text {TeV}\) and 1.34 \(\text {TeV}\) (1.35 and 1.3 \(\text {TeV}\) for the on-Z signal region), respectively, for \(\tilde{\chi }_1^0\) masses below 400 \(\text {GeV}\). The sensitivity in the \(Z^{(*)}\) model is smaller than that of the slepton model because the leptonic branching fraction of the Z boson suppresses the signal production rate.
Model-independent upper limits at 95% CL on the number of events that could be attributed to non-SM sources (\(S^{95}\)) for SRZ are derived using the \(CL_{\text {S}}\) prescription and neglecting possible signal contamination in the CRs. For these upper limits, pseudo-experiments are used rather than the asymptotic approximation. The expected and observed upper limits are given in Table 9. The same information is given for the 24 \(m_{\ell \ell }\) ranges of the edge search in Table 11.
11 Conclusion
This paper presents two searches for new phenomena in final states containing a same-flavour opposite-sign lepton (electron or muon) pair, jets, and large missing transverse momentum using \(14.7~\mathrm {fb}^{-1}\) of ATLAS data collected during 2015 and 2016 at the LHC at \(\sqrt{s}=13\) \(\text {TeV}\). The first search (on-shell Z search) targets lepton pairs consistent with Z boson decay, while the second search (edge search) targets a kinematic endpoint feature in the dilepton mass distribution. For the edge search, a set of 24 mass ranges are considered, with different requirements on \(E_{\text {T}}^{\text {miss}}\) and \(H_{\text {T}} \), and different kinematic endpoint values in the dilepton invariant-mass distribution. The data in both searches are found to be consistent with the Standard Model prediction. The results are interpreted in simplified models of gluino-pair production and squark-pair production, and exclude gluinos (squarks) with masses as large as 1.7 \(\text {TeV}\) (980 \(\text {GeV}\)).
Notes
ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upward. Cylindrical coordinates \((r,\phi )\) are used in the transverse plane, \(\phi \) being the azimuthal angle around the z-axis. The pseudorapidity is defined in terms of the polar angle \(\theta \) as \(\eta =-\ln \tan (\theta /2)\) and the rapidity is defined as \(y=1/2 \cdot \ln [(E+p_{z})/(E-p_{z})])\), where E is the energy and \(p_{z}\) the longitudinal momentum of the object of interest. The opening angle between two analysis objects in the detector is defined as \(\Delta R=\sqrt{(\Delta y)^2+(\Delta \phi )^2}\).
The distance of closest approach between a particle object and the primary vertex (beamline) in the longitudinal (transverse) plane is denoted by \(z_0\) (\(d_0\)).
References
Y.A. Golfand, E.P. Likhtman, Extension of the algebra of Poincare group generators and violation of p invariance. JETP Lett. 13, 323–326 (1971)
Y.A. Golfand, E.P. Likhtman, Pisma. Zh. Eksp. Teor. Fiz. 13, 452 (1971)
D.V. Volkov, V.P. Akulov, Is the neutrino a goldstone particle? Phys. Lett. B 46, 109–110 (1973). doi:10.1016/0370-2693(73)90490-5
J. Wess, B. Zumino, Supergauge transformations in four-dimensions. Nucl. Phys. B 70, 39–50 (1974). doi:10.1016/0550-3213(74)90355-1
J. Wess, B. Zumino, Supergauge invariant extension of quantum electrodynamics. Nucl. Phys. B 78, 1 (1974). doi:10.1016/0550-3213(74)90112-6
S. Ferrara, B. Zumino, Supergauge invariant Yang–Mills theories. Nucl. Phys. B 79, 413 (1974). doi:10.1016/0550-3213(74)90559-8
A. Salam, J.A. Strathdee, Supersymmetry and nonabelian gauges. Phys. Lett. B 51, 353–355 (1974). doi:10.1016/0370-2693(74)90226-3
N. Sakai, Naturalness in supersymmetric GUTS. Z. Phys. C 11, 153 (1981). doi:10.1007/BF01573998
S. Dimopoulos, S. Raby, F. Wilczek, Supersymmetry and the scale of unification. Phys. Rev. D 24, 1681–1683 (1981). doi:10.1103/PhysRevD.24.1681
L.E. Ibanez, G.G. Ross, Low-energy predictions in supersymmetric grand unified theories. Phys. Lett. B 105, 439 (1981). doi:10.1016/0370-2693(81)91200-4
S. Dimopoulos, H. Georgi, Softly broken supersymmetry and SU(5). Nucl. Phys. B 193, 150 (1981). doi:10.1016/0550-3213(81)90522-8
G.R. Farrar, P. Fayet, Phenomenology of the production, decay, and detection of new hadronic states associated with supersymmetry. Phys. Lett. B 76, 575–579 (1978). doi:10.1016/0370-2693(78)90858-4
H. Goldberg, Constraint on the photino mass from cosmology. Phys. Rev. Lett. 50, 1419 (1983). doi:10.1103/PhysRevLett.50.1419. [Erratum: Phys. Rev. Lett. 103, 099905 (2009)]
J.R. Ellis et al., Supersymmetric relics from the big bang. Nucl. Phys. B 238, 453–476 (1984). doi:10.1016/0550-3213(84)90461-9
M. Dine, W. Fischler, A Phenomenological model of particle physics based on supersymmetry. Phys. Lett. B 110, 227 (1982). doi:10.1016/0370-2693(82)91241-2
L. Alvarez-Gaume, M. Claudson, M.B. Wise, Low-energy supersymmetry. Nucl. Phys. B 207, 96 (1982). doi:10.1016/0550-3213(82)90138-9
C.R. Nappi, B.A. Ovrut, Supersymmetric extension of the SU(3) x SU(2) x U(1) model. Phys. Lett. B 113, 175 (1982). doi:10.1016/0370-2693(82)90418-X
CMS Collaboration, Search for physics beyond the standard model in events with a \(Z\) boson, jets, and missing transverse energy in \(pp\) collisions at \(\sqrt{s}=7\) TeV. Phys. Lett. B 716, 260–284 (2012). doi:10.1016/j.physletb.2012.08.026. arXiv:1204.3774 [hep-ex]
CMS Collaboration, Search for physics beyond the standard model in events with two leptons, jets, and missing transverse momentum in pp collisions at \(\sqrt{s}\) = 8 TeV. JHEP 04, 124 (2015). doi:10.1007/JHEP04(2015)124. arXiv:1502.06031 [hep-ex]
ATLAS Collaboration, Search for supersymmetry in events containing a same-flavour opposite-sign dilepton pair, jets, and large missing transverse momentum in \(\sqrt{s} = 8\;TeV\) \(pp\) collisions with the ATLAS detector. Eur. Phys. J. C 75, 318 (2015). doi:10.1140/epjc/s10052-015-3518-2. arXiv:1503.03290 [hep-ex]
CMS Collaboration, Search for new physics in final states with two opposite-sign same-flavor leptons, jets and missing transverse momentum in pp collisions at \(\sqrt{s}=13\) TeV. JHEP 12, 013 (2016). doi:10.1007/JHEP12(2016)013. arXiv:1607.00915 [hep-ex]
ATLAS Collaboration, Search for new physics in events with opposite-sign leptons, jets, and missing transverse energy in \(pp\) collisions at \(\sqrt{s}=7\) TeV. Phys. Lett. B 718, 815–840 (2013). doi:10.1016/j.physletb.2012.11.036. arXiv:1206.3949 [hep-ex]
ATLAS Collaboration, The ATLAS experiment at the CERN large hadron collider. JINST. 3, S08003 (2008). doi:10.1088/1748-0221/3/08/S08003
ATLAS Collaboration, Early inner detector tracking performance in the 2015 data at \(\sqrt{s}\) = 13 TeV. ATL-PHYS-PUB-2015-051, 2015. http://cds.cern.ch/record/2110140
P. Fayet, Supersymmetry and weak, electromagnetic and strong interactions. Phys. Lett. B 64, 159 (1976). doi:10.1016/0370-2693(76)90319-1
P. Fayet, Spontaneously broken supersymmetric theories of weak, electromagnetic and strong interactions. Phys. Lett. B 69, 489 (1977). doi:10.1016/0370-2693(77)90852-8
M. Cahill-Rowley et al., ATLAS Z+ missing transverse energy excess in the MSSM. Phys. Rev. D 92, 075029 (2015). doi:10.1103/PhysRevD.92.075029. arXiv:1506.05799 [hep-ph]
ATLAS Collaboration, Search for direct top squark pair production in events with a \(Z\) boson, \(b\)-jets and missing transverse momentum in \(\sqrt{s} = 8\;TeV\) \(pp\) collisions with the ATLAS detector. Eur. Phys. J. C 74, 2883 (2014). doi:10.1140/epjc/s10052-014-2883-6. arXiv:1403.5222 [hep-ex]
ATLAS Collaboration, Search for direct top-squark pair production in final states with two leptons in \(pp\) collisions at \(\sqrt{s} = 8\;TeV\) with the ATLAS detector. JHEP 06, 124 (2014). doi:10.1007/JHEP06(2014)124. arXiv:1403.4853 [hep-ex]
D. Alves, E. Izaguirre, J. Wacker, Where the sidewalk ends: jets and missing energy search strategies for the 7 TeV LHC. JHEP 10, 012 (2011). doi:10.1007/JHEP10(2011)012. arXiv:1102.5338 [hep-ph]
ATLAS Collaboration, Improved luminosity determination in pp collisions at \(\sqrt{s}\) = 7 TeV using the ATLAS detector at the LHC. Eur. Phys. J. C 73, 2518 (2013). doi:10.1140/epjc/s10052-013-2518-3. arXiv:1302.4393 [hep-ex]
ATLAS Collaboration, Luminosity determination in pp collisions at \(\sqrt{s}=8\) TeV using the ATLAS detector at the LHC. Eur. Phys. J. C 76, 653 (2016). doi:10.1140/epjc/s10052-016-4466-1. arXiv:1608.03953 [hep-ex]
ATLAS Collaboration, 2015 start-up trigger menu and initial performance assessment of the ATLAS trigger using Run-2 data. ATL-DAQ-PUB-2016-001, 2016. http://cds.cern.ch/record/2136007
J. Alwall et al., The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations. JHEP 07, 079 (2014). doi:10.1007/JHEP07(2014)079. arXiv:1405.0301 [hep-ph]
T. Sjöstrand, S. Mrenna, P. Skands, A brief introduction to PYTHIA 8.1. Comput. Phys. Commun. 178, 852 (2008). doi:10.1016/j.cpc.2008.01.036. arXiv:0710.3820 [hep-ph]
ATLAS Collaboration, ATLAS Pythia 8 tunes to 7 TeV data. ATL-PHYS-PUB-2014-021, 2014. http://cdsweb.cern.ch/record/1966419
R.D. Ball et al., Parton distributions with LHC data. Nucl. Phys. B 867, 244–289 (2013). doi:10.1016/j.nuclphysb.2012.10.003. arXiv:1207.1303 [hep-ph]
P. Nason, A new method for combining NLO QCD with shower Monte Carlo algorithms. JHEP 11, 040 (2004). arXiv:hep-ph/0409146 [hep-ph]
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]
S. Alioli et al., A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG BOX. JHEP 06, 043 (2010). arXiv:1002.2581 [hep-ph]
T. Sjöstrand, S. Mrenna, P. Skands, PYTHIA 6.4 physics and manual. JHEP 05, 026 (2006). doi:10.1088/1126-6708/2006/05/026. arXiv:hep-ph/0603175 [hep-ph]
B. Cooper et al., Monte Carlo tuning in the presence of matching. Eur. Phys. J. C 72, 2078 (2011). doi:10.1140/epjc/s10052-012-2078-y. arXiv:1109.5295 [hep-ph]
D.J. Lange, The EvtGen particle decay simulation package. Nucl. Instrum. Methods A 462, 152 (2001). doi:10.1016/S0168-9002(01)00089-4
T. Gleisberg et al., Event generation with Sherpa 1.1. JHEP 02, 007 (2009). arXiv:0811.4622 [hep-ph]
S. Schumann, F. Krauss, A Parton shower algorithm based on Catani–Seymour dipole factorisation. JHEP 03, 038 (2008). doi:10.1088/1126-6708/2008/03/038. arXiv:0709.1027 [hep-ph]
T. Gleisberg, S. Höche, Comix, a new matrix element generator. JHEP 12, 039 (2008). doi:10.1088/1126-6708/2008/12/039. arXiv:0808.3674 [hep-ph]
F. Cascioli, P. Maierhofer, S. Pozzorini, Scattering amplitudes with open loops. Phys. Rev. Lett. 108, 111601 (2012). doi:10.1103/PhysRevLett. 108.111601. arXiv:1111.5206 [hep-ph]
S. Höche et al., QCD matrix elements + parton showers: the NLO case. JHEP 04, 027 (2013). doi:10.1007/JHEP04(2013) 027. arXiv:1207.5030 [hep-ph]
S. Catani et al., QCD Matrix elements + Parton showers. JHEP 11, 063 (2001). doi:10.1088/1126-6708/2001/11/063. arXiv:hep-ph/0109231
L. Lönnblad, Correcting the colour-dipole cascade model with fixed order matrix elements. JHEP 05, 046 (2002). doi:10.1088/1126-6708/2002/05/046. arXiv:hep-ph/0112284
W. Beenakker et al., Squark and gluino production at hadron colliders. Nucl. Phys. B 492, 51–103 (1997). doi:10.1016/S0550-3213(97)00084-9. arXiv:hep-ph/9610490 [hep-ph]
A. Kulesza, L. Motyka, Threshold resummation for squark–antisquark and gluino-pair production at the LHC. Phys. Rev. Lett. 102, 111802 (2009). doi:10.1103/PhysRevLett.102.111802. arXiv:0807.2405 [hep-ph]
A. Kulesza, L. Motyka, Soft gluon resummation for the production of gluino–gluino and squark–antisquark pairs at the LHC. Phys. Rev. D 80, 095004 (2009). doi:10.1103/PhysRevD.80.095004. arXiv:0905.4749 [hep-ph]
W. Beenakker et al., Soft-gluon resummation for squark and gluino hadroproduction. JHEP 12, 041 (2009). doi:10.1088/1126-6708/2009/12/041. arXiv:0909.4418 [hep-ph]
W. Beenakker et al., Squark and gluino hadroproduction. Int. J. Mod. Phys. A 26, 2637–2664 (2011). doi:10.1142/S0217751X11053560. arXiv:1105.1110 [hep-ph]
ATLAS Collaboration, The ATLAS simulation infrastructure. Eur. Phys. J. C 70, 823–874 (2010). doi:10.1140/epjc/s10052-010-1429-9. arXiv:1005.4568 [physics.ins-det]
S. Agostinelli et al., GEANT4: a simulation toolkit. Nucl. Instrum. Methods A 506, 250–303 (2003). doi:10.1016/S0168-9002(03)01368-8
ATLAS Collaboration, Summary of ATLAS Pythia 8 tunes. ATL-PHYS-PUB-2012-003, 2012. http://cds.cern.ch/record/1474107
G. Watt, R.S. Thorne, Study of Monte Carlo approach to experimental uncertainty propagation with MSTW 2008 PDFs. JHEP 08, 052 (2012). doi:10.1007/JHEP08(2012)052. arXiv:1205.4024 [hep-ph]
ATLAS Collaboration, Modelling of the \(t\bar{t}H\) and \(t\bar{t}V\) \((V=W,Z)\) processes for \(\sqrt{s}=13\) TeV ATLAS analyses. ATL-PHYS-PUB-2016-005, 2016. http://cds.cern.ch/record/2120826
M.V. Garzelli et al., \(t\bar{t}\) \(W^{+-}\) and \(t\bar{t}Z\) Hadroproduction at NLO accuracy in QCD with Parton Shower and Hadronization effects. JHEP 11, 056 (2012). doi:10.1007/JHEP11(2012)056. arXiv:1208.2665 [hep-ph]
J.M. Campbell, R.K. Ellis, \(t\bar{t}W\) production and decay at NLO. JHEP 07, 052 (2012). arXiv:1204.5678 [hep-ph]
A. Lazopoulos et al., Next-to-leading order QCD corrections to \(t\bar{t}Z\) production at the LHC. Phys Lett. B 666, 62 (2008). arXiv:0804.2220 [hep-ph]
ATLAS Collaboration, Simulation of top quark production for the ATLAS experiment at \(\sqrt{s}\) = 13 TeV. ATL-PHYS-PUB-2016-004, 2016. http://cds.cern.ch/record/2120417
M. Czakon, P. Fiedler, A. Mitov, Total top-quark pair-production cross section at hadron colliders through \(O(\alpha _s^4)\). Phys. Rev. Lett. 110, 252004 (2013). arXiv:1303.6254 [hep-ph]
M. Czakon, A. Mitov, Top++: a program for the calculation of the top-pair cross-section at hadron colliders. Comput. Phys. Commun. 185, 2930 (2014). doi:10.1016/j.cpc.2014.06.021. arXiv:1112.5675 [hep-ph]
N. Kidonakis, Two-loop soft anomalous dimensions for single top quark associated production with a \(W^-\) or \(H^-\). Phys. Rev. D 82, 054018 (2010). doi:10.1103/PhysRevD.82.054018. arXiv:1005.4451 [hep-ph]
ATLAS Collaboration, Multi-Boson simulation for 13 TeV ATLAS analyses. ATL-PHYS-PUB-2016-002, 2016. http://cds.cern.ch/record/2119986
J.M. Campbell, R.K. Ellis, An update on vector boson pair production at hadron colliders. Phys. Rev. D 60, 113006 (1999). arXiv:hep-ph/9905386 [hep-ph]
J.M. Campbell, R.K. Ellis, C. Williams, Vector boson pair production at the LHC. JHEP 07, 018 (2011). arXiv:1105.0020 [hep-ph]
ATLAS Collaboration, Monte Carlo generators for the production of a \(W\) or \(Z/\gamma ^*\) boson in association with jets at ATLAS in Run 2. ATL-PHYS-PUB-2016-003, 2016. http://cds.cern.ch/record/2120133
S. Catani et al., Vector boson production at hadron colliders: a fully exclusive QCD calculation at NNLO. Phys. Rev. Lett. 103, 082001 (2009). arXiv:0903.2120 [hep-ph]
S. Catani, M. Grazzini, An NNLO subtraction formalism in hadron collisions and its application to Higgs boson production at the LHC. Phys. Rev. Lett. 98, 222002 (2007). arXiv:hep-ph/0703012 [hep-ph]
ATLAS Collaboration, Vertex reconstruction performance of the ATLAS detector at \(\sqrt{s}\) = 13 TeV. ATL-PHYS-PUB-2015-026, 2015. http://cds.cern.ch/record/2037717
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
ATLAS Collaboration, Muon reconstruction performance of the ATLAS detector in proton–proton collision data at \(\sqrt{s} = 13 TeV\). Eur. Phys. J. C 76, 292 (2016). doi:10.1140/epjc/s10052-016-4120-y. arXiv:1603.05598 [hep-ex]
ATLAS Collaboration, Topological cell clustering in the ATLAS calorimeters and its performance in LHC Run 1. (2016). arXiv:1603.02934 [hep-ex]
M. Cacciari, G.P. Salam, G. Soyez, The anti-k\(_t\) jet clustering algorithm. JHEP 04, 063 (2008). doi:10.1088/1126-6708/2008/04/063. arXiv:0802.1189 [hep-ph]
M. Cacciari, G.P. Salam, Dispelling the N\(^{3}\) myth for the Kt jet-finder. Phys. Lett. B 641, 57–61 (2006). doi:10.1016/j.physletb.2006.08.037. arXiv:hep-ph/0512210
ATLAS Collaboration, Jet energy measurement and its systematic uncertainty in proton–proton collisions at \(\sqrt{s} = 7 TeV\) with the ATLAS detector. Eur. Phys. J. C 75, 17 (2015). doi:10.1140/epjc/s10052-014-3190-y. arXiv:1406.0076 [hep-ex]
ATLAS Collaboration, Jet calibration and systematic uncertainties for jets reconstructed in the ATLAS detector at \(\sqrt{s}=13\) TeV. ATL-PHYS-PUB-2015-015, 2015. http://cds.cern.ch/record/2037613
ATLAS Collaboration, Tagging and suppression of pileup jets with the ATLAS detector. ATLAS-CONF-2014-018, 2014. http://cds.cern.ch/record/1700870
ATLAS Collaboration, Characterisation and mitigation of beam-induced backgrounds observed in the ATLAS detector during the 2011 proton–proton run. JINST 8, P07004 (2013). doi:10.1088/1748-0221/8/07/P07004. arXiv:1303.0223 [hep-ex]
ATLAS Collaboration, Selection of jets produced in 13 TeV proton–proton collisions with the ATLAS detector. ATLAS-CONF-2015-029, 2015. http://cds.cern.ch/record/2037702
ATLAS Collaboration, Performance of \(b\)-Jet Identification in the ATLAS experiment. JINST 11, P04008 (2016). doi:10.1088/1748-0221/11/04/P04008. arXiv:1512.01094 [hep-ex]
ATLAS Collaboration, Optimisation of the ATLAS b-tagging performance for the 2016 LHC Run. ATL-PHYS-PUB-2016-012, 2016. http://cds.cern.ch/record/2160731
ATLAS Collaboration, Electron and photon energy calibration with the ATLAS detector using LHC Run 1 data. Eur. Phys. J. C 74, 3071 (2014). doi:10.1140/epjc/s10052-014-3071-4. arXiv:1407.5063 [hep-ex]
ATLAS Collaboration, Expected performance of missing transverse momentum reconstruction for the ATLAS detector at \(\sqrt{s}=13\) TeV. ATL-PHYS-PUB-2015-023, 2015. http://cds.cern.ch/record/2037700
ATLAS Collaboration, Performance of missing transverse momentum reconstruction for the ATLAS detector in the first proton–proton collisions at at \(\sqrt{s}\) = 13 TeV. ATL-PHYS-PUB-2015-027, 2015. http://cds.cern.ch/record/2037904
G. Cowan et al., Asymptotic formulae for likelihood-based tests of new physics. Eur. Phys. J. C 71, 1554 (2011). doi:10.1140/epjc/s10052-011-1554-0. arXiv:1007.1727 [physics.data-an]
ATLAS Collaboration, Measurement of the differential cross-section of highly boosted top quarks as a function of their transverse momentum in \(\sqrt{s}\) = 8 TeV proton–proton collisions using the ATLAS detector. Phys. Rev. D 93, 032009 (2016). doi:10.1103/PhysRevD.93.032009. arXiv:1510.03818 [hep-ex]
C.M.S. Collaboration, Measurement of the integrated and differential t-tbar production cross sections for high-\(p_{\text{T}}\) top quarks in pp collisions at \(\sqrt{s}\) = 8 TeV. Phys. Rev. D 94, 072002 (2016). doi:10.1103/PhysRevD.94.072002. arXiv:1605.00116 [hep-ex]
M. Czakon, D. Heymes, A. Mitov. Dynamical scales for multi-TeV top-pair production at the LHC. (2016). arXiv:1606.03350 [hep-ph]
ATLAS Collaboration, Search for squarks and gluinos in events with isolated leptons, jets and missing transverse momentum at \(\sqrt{s} = 8 TeV\) with the ATLAS detector. JHEP 04, 116 (2015). doi:10.1007/JHEP04(2015)116. arXiv:1501.03555 [hep-ex]
M. Krämer et al., Supersymmetry production cross sections in \(pp\) collisions at \(\sqrt{s} = 7~TeV\). (2012). arXiv:1206.2892 [hep-ph]
C. Borschensky et al., Squark and gluino production cross sections in pp collisions at \(\sqrt{s}\) = 13, 14, 33 and 100 TeV Eur. Phys. J. C 74, 3174 (2014). doi:10.1140/epjc/s10052-014-3174-y. arXiv:1407.5066 [hep-ph]
M. Baak et al., HistFitter software framework for statistical data analysis Eur. Phys. J. C 75, 153 (2014). doi:10.1140/epjc/s10052-015-3327-7. arXiv:1410.1280 [hep-ex]
A. Read, Presentation of search results: the CLs technique. J. Phys. G Nucl. Part. Phys. 28, 2693–2704 (2002). doi:10.1088/0954-3899/28/10/313
ATLAS Collaboration, ATLAS computing acknowledgements 2016–2017. ATL-GEN-PUB-2016-002, 2016. https://cds.cern.ch/record/2202407
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-DSM/IRFU, France; GNSF, 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; FOM and 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 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. [99].
Author information
Authors and Affiliations
Consortia
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.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
About this article
Cite this article
Aaboud, M., Aad, G., Abbott, B. et al. Search for new phenomena in events containing a same-flavour opposite-sign dilepton pair, jets, and large missing transverse momentum in \(\varvec{\sqrt{s}=13}\) \(\text {TeV}\) \(\varvec{pp}\) collisions with the ATLAS detector. Eur. Phys. J. C 77, 144 (2017). https://doi.org/10.1140/epjc/s10052-017-4700-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-017-4700-5