Search for supersymmetry in proton-proton collisions at $\sqrt{s} =$ 13 TeV in events with high-momentum Z bosons and missing transverse momentum

A search for new physics in events with two highly Lorentz-boosted Z bosons and large missing transverse momentum is presented. The analyzed proton-proton collision data, corresponding to an integrated luminosity of 137 fb$^{-1}$, were recorded at $\sqrt{s}= $ 13 TeV by the CMS experiment at the CERN LHC. The search utilizes the substructure of jets with large radius to identify quark pairs from Z boson decays. Backgrounds from standard model processes are suppressed by requirements on the jet mass and the missing transverse momentum. No significant excess in the event yield is observed beyond the number of background events expected from the standard model. For a simplified supersymmetric model in which the Z bosons arise from the decay of gluinos, an exclusion limit of 1920 GeV on the gluino mass is set at 95% confidence level. This is the first search for beyond-standard-model production of pairs of boosted Z bosons plus large missing transverse momentum.


Introduction
The discovery of a Higgs boson in 2012 by the ATLAS and CMS experiments [1][2][3] at the CERN LHC fulfilled the predicted particle content of the standard model (SM). However, within the SM as a quantum field theory, the measured Higgs boson mass of around 125 GeV presents a special challenge as the calculated mass is unstable against corrections from loop processes when the theory is extended to higher mass scales. In the absence of extreme fine tuning [4][5][6][7] that would precisely cancel the divergent terms, the mass value can run up to the ultraviolet cutoff of the model at the Planck scale. This instability of the Higgs boson mass and the entire electroweak scale is known as the gauge hierarchy problem.
One widely studied extension of the SM is supersymmetry (SUSY) [8-10], which posits a partner for each SM particle differing in spin by one-half unit. For example, squarks q and gluinos g are the SUSY partners of quarks and gluons, respectively. Depending on the mass hierarchy of these new particles, they could resolve the gauge hierarchy problem by providing necessary radiative corrections to partly cancel the SM contributions. Furthermore, in R-parity conserving models [11,12], SUSY particles are produced in pairs, while the lightest of them is neutral, stable, and weakly interacting. This lightest SUSY particle (LSP) provides a suitable candidate for dark matter [12], which is not described in the SM. The typical experimental signatures of pair-produced SUSY particles with R-parity conserving decay chains are jets, leptons, and large missing transverse momentum (p miss T ). As gluinos and squarks carry color charge, like their SM partners, they can be produced via the strong interaction; therefore among SUSY particles they have the highest production cross sections at hadron colliders for a given mass. Searches for direct decays of gluinos to quarks and the LSP have excluded m( g ) 2 TeV [13][14][15][16], depending on the model. The search described in this paper focuses on gluino decay cascades to Z bosons and the LSP via the next-to-lightest SUSY particle (NLSP). We consider a picture in which the NLSP and LSP are respectively the neutralinos χ 0 2 and χ 0 1 , mixed states of SUSY partners of the neutral Higgs and gauge bosons. Such a situation arises in SUSY scenarios like those described in Ref.
[17] that seek to preserve "naturalness," that is, minimal fine tuning of the SM to solve the gauge hierarchy problem, by admitting large mass splittings among the neutralinos (and charginos), leading to experimental signatures with vector bosons and p miss T in the final state. Figure 1 shows our signal process, expressed within the framework of simplified models [18][19][20][21], and referred to as T5ZZ. We further assume a heavy χ 0 2 , (with mass below that of the g), and a light χ 0 1 . This gives rise to energetic Z bosons along with large p miss T and additional soft quarks in the final state. In our model calculations we set the branching fraction for χ 0 2 → Z χ 0 1 to 100%, the χ 0 1 mass to 1 GeV, and the difference in mass between the g and χ 0 2 to 50 GeV, though any set of mass parameters with a large [O(TeV)] mass difference between the χ 0 2 and χ 0 1 will result in highly energetic Z bosons. For the dominant Z → qq decay at large momentum, the decay products can be contained in a single reconstructed jet with a large angular radius (wide-cone jet).
In this paper, we present a search in proton-proton (pp) collisions at √ s = 13 TeV for events with two highly Lorentz-boosted, hadronically decaying Z bosons and large p miss T . The analysis is based on the LHC Run 2 data set with an integrated luminosity of 137 fb −1 , recorded by the CMS experiment during 2016-2018. The signature for a signal is a pair of wide-cone jets, each having a reconstructed mass consistent with the Z boson mass. This selection, in combination with large p miss T , greatly suppresses backgrounds from SM processes. The assumed small mass splitting between the g and χ 0 2 implies a massive χ 0 2 . We further assume a 100% branching fraction for the χ 0 2 decay to the Z boson and χ 0 1 , leading to an energetic Z boson and large p miss T .

The CMS detector and trigger
A detailed description of the CMS detector and the associated coordinate system and kinematic variables is given in Ref. [22]. The main components of the apparatus are briefly discussed here. The core of CMS is a cylindrical superconducting solenoid with an inner diameter of 6 m that provides a 3.8 T axial magnetic field. A silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter, and a brass and scintillator hadron calorimeter are placed within the volume enclosed by the solenoid. Gas-ionization detectors are embedded in the steel fluxreturn yoke outside the solenoid to identify muons. The detector is nearly hermetic, permitting accurate measurements of p miss T . The CMS trigger system is described in Ref. [23]. For this analysis, signal candidate events were recorded by requiring p miss T at the trigger level to exceed a threshold that varied between 100 and 120 GeV, depending on the LHC instantaneous luminosity. The efficiency of this trigger is measured in data to be greater than 97% for events satisfying the selection criteria described in Section 5. Additional triggers based on an isolated lepton or photon are used to select control samples for the background predictions.

Simulated event samples
The estimation of yields for the most prominent backgrounds is based on data in orthogonal signal-depleted control regions and is described in Section 6. Samples of Monte Carlo (MC) simulated events are used to test the background estimation, as well as to optimize the selection criteria. These samples include events with top quark pair production (tt), and photon, W boson, or Z boson production accompanied by jets, denoted γ+jets, W+jets, or Z+jets, respectively.
The SM production of tt, γ+jets, W+jets, Z+jets, and quantum chromodynamics (QCD) multijet events is simulated using the MADGRAPH5 aMC@NLO 2.2.2 [24,25] generator for 2016 samples and MADGRAPH5 aMC@NLO 2.4.2 for 2017 and 2018 samples, all with leading order (LO) precision. The tt events are generated with up to three additional partons in the matrix element calculations, while the γ+jets, W+jets, and Z+jets events are generated with up to four additional partons. Single top quark events produced via the s channel, diboson events originating from WW, ZZ, or ZH production, and events from ttW, ttZ, and WWZ production, are generated with MADGRAPH5 aMC@NLO 2.2.2 at next-to-leading order (NLO) [26], except that WW events in which both W bosons decay leptonically are generated using POWHEG 2.0 [27][28][29][30][31] at NLO. The POWHEG generator is also used to describe t-channel production of single top quarks as well as tW events. Normalization of the simulated background samples is derived from the most accurate cross section calculations available [24,[30][31][32][33][34][35][36][37][38][39][40], which generally correspond to NLO or next-to-NLO (NNLO) precision.
All simulated samples make use of PYTHIA 8.205 (2016) or 8.230 (2017 and 2018) [53] to describe parton showering and hadronization. The CUETP8M1 [54] tune was used to simulate both the SM background and signal samples for the 2016 simulation. To generate the 2017 and 2018 samples, PYTHIA was used, with the CP5 tune [55] for the backgrounds and the CP2 tune [55] for signals. Simulated samples generated at LO (NLO) with the CUETP8M1 tune use the NNPDF3.0LO (NNPDF3.0NLO) [56] PDF set, while those generated with the CP2 or CP5 tune use the NNPDF3.1LO (NNPDF3.1NNLO) [57] PDF set. Here PDF refers to the parton distribution function. The detector response is modeled with GEANT4 [58]. The simulated events are generated with a distribution of pp interactions per bunch crossing ("pileup") that is adjusted to match the corresponding distribution measured in data.
To improve the description of initial-state radiation (ISR), the MADGRAPH5 aMC@NLO prediction of the jet multiplicity distribution is compared with data in a control sample enriched in tt events [13]. A correction factor derived therefrom is subsequently applied to the simulated tt and signal events. The correction is found to be unnecessary for tt samples that are generated with the CP5 tune, so it is not applied to those samples.

Event reconstruction
Individual particles are reconstructed with the CMS particle-flow (PF) algorithm [59], which identifies them as photons, charged or neutral hadrons, electrons, or muons. These objects are characterized kinematically by their transverse momentum p T , pseudorapidity η, and azimuthal angle φ. Photon and electron candidates are required to satisfy |η| < 2.5, and muon candidates |η| < 2.4, within the fiducial coverage of the tracking and muon system, respectively.
The missing transverse momentum vector p miss T is computed as the negative vector sum of the p T of all of the PF candidates in an event, and its magnitude is denoted as p miss T [60]. The p miss T is modified to account for corrections to the energy scale of the reconstructed jets in the event.
The reconstructed vertex with the largest value of summed physics-object p 2 T is taken to be the primary pp interaction vertex, where the physics objects are the jets, clustered using the antik T algorithm [61,62] with the charged particle tracks assigned to the vertex as inputs, and the associated missing transverse momentum, taken as the negative vector p T sum of those jets. Charged particle tracks associated with vertices other than the primary vertex are removed from further consideration.
Jets are defined as clusters of PF candidates formed by the anti-k T algorithm with a distance parameter of 0.4 or 0.8. Quality criteria [63,64] are imposed to suppress jets from spurious sources such as electronics noise in the calorimeters. The jet energies are corrected for the nonlinear response of the detector [65]. Jets with p T > 30 GeV, |η| < 2.4, and a distance parameter of 0.4 (AK4) are used as specified in Section 5 to calculate some of the selection variables. For these jets, charged particles that emerge from vertices other than the primary one are removed from the list of PF candidates used for the jet clustering. The expected contribution from neutral particles from pileup is removed using the effective area technique [64,66].
The hadronically decaying Z boson candidates are reconstructed as wide-cone jets with a distance parameter of 0.8 (AK8). These AK8 jets are reclustered from their original constituents using the "soft drop" method [67] to remove soft, wide-angle radiation that can adversely impact the mass measurement of the jet. Contributions from pileup in these jets are removed with the PUPPI technique [68]. The soft drop mass m jet is then used to identify jets from Z → qq decays. No requirements on their flavor content are imposed.
The identification of b jets (b jet tagging) is performed by applying, to the AK4 jets, a version of the combined secondary vertex algorithm based on deep neural networks [69] (DeepCSV). A working point ("medium") of this algorithm is used that has a tagging efficiency for b jets of 68%, and a misidentification probability of approximately 1% for gluon and light-flavor quark jets and 12% for charm quark jets.
As described in Section 5, events with leptons or photons are vetoed in the search sample selection. Electron and muon candidates are identified as described in Refs.
[70] and [71], respectively. To suppress jets erroneously identified as leptons or genuine leptons from hadron decays, electron and muon candidates are subjected to an isolation requirement. The isolation criterion is based on a variable I, which is the scalar p T sum of charged hadron, neutral hadron, and photon PF candidates within a cone of radius ∆R = √ (∆η) 2 + (∆φ) 2 around the lepton direction, divided by the lepton p T . The expected contributions of neutral particles from pileup are subtracted [64,66]. The radius of the cone, in radians, is 0.2 for lepton p T < 50 GeV, 10 GeV/p T for 50 ≤ p T ≤ 200 GeV, and 0.05 for p T > 200 GeV. The decrease in cone size with increasing lepton p T accounts for the increased collimation of the decay products from the lepton's parent particle as the Lorentz boost of the latter increases [72]. The isolation requirement is I < 0.1 (0.2) for electrons (muons).
To further suppress events with leptons from hadron decays and single-prong hadronic τ lepton decays, the event selection veto is extended to include isolated charged-particle tracks not identified as electrons or muons by the criteria of the previous paragraph. For these candidates the scalar p T sum of all other charged-particle tracks within ∆R = 0.3 around the track direction, divided by the track p T , is required to be less than 0.2 if the track is identified as a PF electron or muon, and less than 0.1 otherwise. Isolated tracks are required to satisfy |η| < 2.4.
Photon candidates are identified as described in Ref.
[73], using the "loose" working point, and with an isolation requirement based on the individual sums of energy from charged and neutral hadrons and electromagnetically interacting particles, excluding the photon candidate itself, within ∆R = 0.3 around the direction of the photon candidate. Each of the three individual sums, corrected for pileup, is required not to exceed a threshold that depends on the calorimeter geometry.

Event selection
We select events with large jet activity and p miss T , no leptons or photons, and wide-cone jets from Lorentz-boosted, hadronically decaying Z bosons. Control regions for the determination of backgrounds are also defined.
The observables used to characterize candidate events are: • N jet , the number of AK4 jets in the event; • ∆φ j, H miss T , the azimuthal angle between the p T of the j th AK4 jet and H miss • ∆R Z,b , the angular separation between a wide-cone jet and a b-tagged jet.
The following requirements define the event selection: 3) for the first two (up to next two, if N jet > 2) AK4 jets ranked in descending order of p T ; 5. no identified isolated photon, electron, or muon candidate with p T > 10 GeV; 6. no isolated track with m T < 100 GeV and p T > 5 GeV if the track is identified as a PF electron or muon, 10 GeV otherwise.
7. at least two AK8 jets with p T > 200 GeV; 8. m jet of the two highest p T AK8 jets between 40 and 140 GeV; 9. ∆R Z,b > 0.8, for the second-highest p T AK8 jet and any b-tagged jet.
The ∆φ j, H miss T requirements suppress background from QCD multijet events, as well as those from hadronic Z and W boson decay, for which H miss T is usually aligned along a jet direction. The m T requirement restricts the isolated track veto to situations consistent with a W boson decay.
The first six requirements define an inclusive "hadronic baseline" selection, and the last three specify the further selection of events with jet pairs that include pairs of hadronically decaying Z boson candidates. The accepted range in m jet is chosen to reject the bulk of nonresonant SM processes on the low side, and the peak from boosted top quark jets on the high side, while including sidebands around the Z boson peak to facilitate the determination of the background. The ∆R Z,b requirement suppresses backgrounds from tt and single top quark events in which a top quark is reconstructed as a b-tagged jet together with a W boson reconstructed as an AK8 jet. Figure 2 shows the simulated SM background components and two example signal mass points for events selected without and with the three Z boson requirements. The main sources of SM background are Z+jets, W+jets, and tt, which can yield large p miss T accompanied by AK8 jets formed from random combinations of hadrons. In the case of Z+jets, large p miss T comes from the Z → νν decay. For W+jets and tt, p miss T arises from a leptonically decaying W boson where the charged lepton is undetected. Smaller background contributions arise from the QCD multijet

Background estimation
This section focuses on the estimation of SM backgrounds in each p miss T bin. We first describe the method based on control samples in data, then follow with a description of the performance of the method in simulation (MC closure), and lastly deal with the uncertainty in the p miss T dependence (shape uncertainty) based on the data observed in the validation samples.

Background estimation method
Control regions (CRs) are formed from the events in which one or both of the highest p T (leading) and second-highest p T (subleading) jets lie in the m jet sideband [40,70] ∨ [100, 140] GeV. Figure 3 shows the definition of the SR and CRs in the plane of jet masses of the leading and subleading jets. In addition, validation samples are selected by inverting the lepton or photon veto requirement.
The first step of the method is to determine the background normalization B norm integrated over all p miss T bins above 300 GeV. We fit the m jet distribution for the leading jet in the leadingjet mass sideband, defined as the sample having the subleading jet m jet within, and the leading jet m jet outside, the Z signal window. The bulk of the background is from nonresonant SM contributions, which can be modeled with a smoothly falling shape. The nominal fit is performed with a linear function, as shown in Fig. 4.   Figure 4: Leading AK8 jet m jet shape fit in the mass sidebands. The Z candidate selection is applied and the subleading AK8 jet m jet value is required to lie in the Z signal window. The blue hatched region represents the ±1 standard deviation uncertainty in the fit to the mass sideband performed with a linear function, which is indicated by the blue line. The stacked histogram shows the background from simulation scaled to the data. Expected signal contributions for two example mass points are also shown.

CMS
The uncertainties in B norm include a statistical component from the fit, and a systematic one due to the choice of the fitting function. To obtain the statistical uncertainty due to the interpolation of the fit into the SR, pseudo-experiments generated from the background model are fitted using a linear function with free slope and normalization. The Gaussian width of the resulting distribution of the yields in the Z signal window, 10.7 events, is taken as the statistical uncertainty in the total background prediction.
To test if the linear function is adequate to represent the m jet distribution, we consider higherorder polynomials as alternative functions. We check Chebyshev polynomials of up to the fourth order. The largest variation in the fitted yield with respect to the nominal one, 10.9 events, comes from a fit with a third-order Chebyshev polynomial, and is taken as an additional uncertainty attributable to the fit shape. Considering the statistical uncertainty described above, this results in B norm = 325 ± 15.
To determine the distribution of background events in the p miss T bins, we rely on an underlying assumption that p miss T and m jet have minimal correlation. To derive the p miss T shape in the SR, a nonoverlapping CR is used in which both leading and subleading AK8 jets have m jet in the mass sideband. This is referred to as the p miss T CR (Fig. 3). In each of the six p miss T bins, we calculate the background prediction as where N CR i is the yield in p miss T bin i in the p miss T CR, and the transfer factor, scales the p miss T CR yield to that of the SR. The uncertainty in T includes both statistical and systematic uncertainties in B norm .

Background closure in simulation
The background estimation method based on control samples in data is tested by applying the procedure to MC simulation. We perform this closure test in two steps.
The main assumption to verify is the lack of correlation between the AK8 jet mass and p miss T shape. Figure 5 shows the results of a test of this assumption, where the simulated sample size permits a distribution in relatively fine steps. The plots compare the p miss T shape in the search and control regions, for the two main background processes. In both cases we see that the p miss T shapes are consistent between the two regions.
For the closure test of the background estimation method we calculate the background prediction in each p miss T bin [Eq. (1)] and compare these predictions with the background yields taken directly from simulation. The results of this test, shown in Fig. 6, demonstrate good agreement within the statistical precision of the test. To account for the uncertainties in the comparison, we assign the relative difference between the prediction and direct observation as a nonclosure systematic uncertainty in the p miss T shape. This difference ranges from 1 to 20%, where the variations in the four lower p miss T bins are treated as being anti-correlated with those in the higher p miss T bins to give a systematic uncertainty in the p miss T shape that does not affect the overall normalization of the background estimation.

T shape uncertainty
While the background estimation method is shown to close well in simulation, we additionally verify in data how well the p miss T CR models the p miss T shape in the Z signal window. In   Figure 6: Results of the closure test in which the background estimation method based on control samples in data is applied to simulation and compared with the direct yield, in the analysis search bins. Expected signal contribution for one example mass point is also shown. The lower panel shows the ratio of the prediction to the direct yield. The gray band shows the statistical uncertainty in the direct yield, and the error bars on the points represent the total uncertainty in the prediction. particular, two validation samples are used to compare the p miss T shape obtained from the p miss T CR with the one obtained in the Z signal window, used to define our SR, for the main background components. A photon validation sample is used as a proxy for the Z+jets background component, while a single-lepton sample is used to validate the modeling of tt and W+jets combined.
We select the photon validation sample from events recorded with a single-photon trigger, replacing the photon veto with the requirement of exactly one photon, defined as in Section 4. The photon p T is used to emulate the p miss T from the Z boson when the latter decays to neutrinos. The lower-p T trigger threshold for the photon compared with the p miss T threshold in the signal trigger allows us to consider the photon validation sample down to 200 GeV in photon p T as a proxy for p miss T . To enhance the event count in this sample, we do not require a threshold on ∆R Z,b since there is a low risk of heavy flavor contamination. All other event selection requirements are the same as for the SR of the analysis.
For the single-lepton sample, the same p miss T trigger is used as for the SR. The same offline criteria are also applied, with the exception that the p miss T requirement is relaxed to 200 GeV to gain a longer lever arm for the p miss T shape comparison, and the lepton vetoes are applied only after selecting exactly one electron or muon. Figure 7 shows the p miss T shape comparison for the photon and single-lepton data. Both ratios are consistent with being independent of p miss T , as expected from the MC closure test, albeit within the limited statistical precision of the data. To account for possible shape differences between the search and control regions, we apply a systematic uncertainty in the p miss T shape calculated using the photon and single-lepton samples. The uncertainty is the difference with respect to a uniform distribution of a fit to the SR/CR distribution with a linear function having a free slope parameter. This results in uncertainties ranging from 0-33% in the Z+jets background based on the photon validation sample, and 1-14% in the combined tt and W+jets background based on the single-lepton validation sample. Weighting these by the proportions of those components in the total background yields uncertainties of 2-30%, depending on the p miss T bin.

Systematic uncertainties
The uncertainties in the SM background prediction are described in Section 6, along with the description of the background estimation method. The uncertainties in the background normalization include the statistical uncertainty from the mass sideband fit interpolation as well as the systematic one derived from alternative fit functions. The uncertainties in the p miss T shape include the statistical uncertainties of the p miss T CR. The systematic uncertainties only affect the p miss T shape without changing the background normalization. These are derived from the MC closure test and data validation samples. All of these systematic uncertainties are summarized in the upper section of Table 1.
The sources of uncertainty in the signal efficiency affect the signal normalization, the signal p miss T shape, or both, as indicated in Table 1 , and 2018, respectively. The trigger, lepton veto, and isolated-track veto efficiencies are measured in data validation samples and their statistical uncertainties propagated to the signal yields. The ISR modeling in the simulation is adjusted to match the efficiencies measured in data events enriched in dileptonic tt production and decay, and the uncertainty in this correction is propagated to the signal yields. To evaluate the uncertainty associated with the renormalization (µ R ) and factorization (µ F ) scales, each  scale is varied independently by a factor of 2.0 and 0.5 [78,79]. Uncertainties in the simulation of pileup are found to be of the order of 0.02%; thus no associated uncertainty is applied.
The jet momenta in MC samples are smeared to match the jet energy resolution (JER) in data. The jet energy corrections (JECs) are varied using p T -and η-dependent uncertainties. Both effects are propagated to the jet-dependent variables, including p miss T , H T , and ∆φ j, H miss T , and are varied within the uncertainty of the corrections to derive a systematic uncertainty in the signal yields. The efficiency of the jet quality requirements used to suppress events with misreconstructed jets is found to differ by 1% between data and simulation, and this is applied as a systematic uncertainty. The difference in the resolution of m jet between data and simulation is applied as a smearing factor to the MC events, and the statistical uncertainty in the size of the correction is included as a systematic uncertainty in the corresponding selection efficiency. Lastly, the statistical precision due to the limited event count in the simulated samples is accounted for as an uncertainty.
The systematic uncertainties associated with the signal yields are evaluated assuming that the contributions from the three years of data taking are fully correlated. The total systematic uncertainties in the signal yields range from 0.2 to 6%.

Results
The background predictions and observed yields for each p miss T bin are shown in Fig. 8 and Table 2. The table also gives the inputs to the prediction calculation, Eq. (1). The observations are found to be consistent with the SM predictions within uncertainties, and no evidence for SUSY is observed. We calculate upper limits on the gluino pair-production cross section using a maximum-likelihood fit in which the free parameters are the signal strength µ and the nuisance parameters associated with the systematic uncertainties in the background and signal model. The uncertainty in the normalization of the background is represented with a lognormal function correlated across all p miss T bins, while the p miss T CR statistical uncertainties are assigned as uncorrelated. The MC closure and data-MC agreement uncertainties are assigned as correlated across p miss T bins.
We evaluate 95% confidence level (CL) upper limits based on the asymptotic form of a likelihood ratio test statistic [80], in conjunction with the CL s criterion described in Refs. [81][82][83]. The test statistic is q(µ) = −2 ln(L µ /L max ), where L µ is the maximum likelihood for fixed µ, and L max is the same determined by allowing all parameters, including µ, to vary.
Expected and observed 95% CL upper limits, and the predicted gluino pair-production cross sections, are shown in Fig. 9, taking m( χ 0 1 ) = 1 GeV and m( g ) − m( χ 0 2 ) = 50 GeV. The observed (expected) gluino mass limits reach as high as 1920 (2060) GeV. The observed limit is 1.4 standard deviations weaker than the expected one due to the mild excesses observed in the two highest p miss T bins. The sensitivity of the search is independent of m( χ 0 1 ) values that are small compared with m( χ 0 2 ), and of m( χ 0 2 ) values large enough to ensure Lorentz-boosted Z boson daughters. A gradual loss of signal efficiency occurs with increasing ∆m( g, χ 0 2 ) as quarks from the gluino decay that form AK8 jets with p T above the 200 GeV threshold displace Z jets as leading or subleading in p T .

Summary
Results are presented of a search for events with two hadronically decaying, highly energetic Z bosons and large transverse momentum imbalance, in proton-proton collisions at √ s = 13 TeV. The sample corresponds to an integrated luminosity of 137 fb −1 . The signature for a Z boson candidate is a wide-cone jet having a measured mass compatible with the Z boson mass. Yields   2 s.d. ± Figure 9: The 95% CL upper limit on the production cross section for the T5ZZ signal model as a function of the gluino mass. The solid black curve shows the observed exclusion limit. The dashed black curve presents the expected limit while the green and yellow bands represent the ±1 and ±2 standard deviation uncertainty ranges. The approximate-NNLO+NNLL cross sections [41][42][43][44][45] are shown in the solid blue curve while the dashed blue curves show their theoretical uncertainties [84]. The T5ZZ model assumes a 100% branching fraction for the χ 0 2 to decay to the Z boson and χ 0 1 .
from standard model background processes, which are small for events with the largest transverse momentum imbalance, are estimated from the data in jet mass sidebands. No evidence for physics beyond the standard model is observed. The reach of the search is interpreted in a simplified supersymmetric model of gluino pair production in which each gluino decays to a low-momentum quark pair and the next-to-lightest supersymmetric particle (NLSP), and the latter decays to a Z boson and the lightest supersymmetric particle (LSP). With the further assumption of a large mass splitting between the NLSP and LSP, the data exclude gluino masses below 1920 GeV at 95% confidence level. This is the first search for beyond-standard-model production of pairs of boosted Z bosons plus large missing transverse momentum.

Acknowledgments
We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centers and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses.  [9] H. P. Nilles, "Supersymmetry, supergravity and particle physics", Phys. Rept. 110 (1984) 1, doi:10.1016/0370-1573(84)90008-5.