Search for supersymmetry in events with a photon, jets, b\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {b}$$\end{document}-jets, and missing transverse momentum in proton–proton collisions at 13Te\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\text {Te}\text {V}$$\end{document}

A search for supersymmetry is presented based on events with at least one photon, jets, and large missing transverse momentum produced in proton–proton collisions at a center-of-mass energy of 13Te\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\text {Te}\text {V}$$\end{document}. The data correspond to an integrated luminosity of 35.9fb-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\text {fb}^{-1}$$\end{document} and were recorded at the LHC with the CMS detector in 2016. The analysis characterizes signal-like events by categorizing the data into various signal regions based on the number of jets, the number of b\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {b}$$\end{document}-tagged jets, and the missing transverse momentum. No significant excess of events is observed with respect to the expectations from standard model processes. Limits are placed on the gluino and top squark pair production cross sections using several simplified models of supersymmetric particle production with gauge-mediated supersymmetry breaking. Depending on the model and the mass of the next-to-lightest supersymmetric particle, the production of gluinos with masses as large as 2120Ge\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\text {Ge}\text {V}$$\end{document} and the production of top squarks with masses as large as 1230Ge\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\text {Ge}\text {V}$$\end{document} are excluded at 95% confidence level.


Introduction
The standard model (SM) of particle physics successfully describes many phenomena, but lacks several necessary elements to provide a complete description of nature, including a source for the relic abundance of dark matter (DM) [1,2] in the universe. In addition, the SM must resort to fine tuning [3][4][5][6] to explain the hierarchy between the Planck mass scale and the electroweak scale set by the vacuum expectation value of the Higgs field, the existence of which was recently confirmed by the observation of the Higgs boson (H) [7,8]. Supersymmetry (SUSY) [9][10][11][12][13][14][15][16] is an extension of the SM that can provide both a viable DM candidate and additional particles that inherently cancel large quantum corrections to the Higgs boson mass-squared term from the SM fields. e-mail: cms-publication-committee-chair@cern.ch Supersymmetric models predict a bosonic superpartner for each SM fermion and a fermionic superpartner for each SM boson; each new particle's spin differs from that of its SM partner by half a unit. SUSY also includes a second Higgs doublet. New colored states, such as gluinos ( g) and top squarks ( t), the superpartners of the gluon and the top quark, respectively, are expected to have masses on the order of 1 TeV to avoid fine tuning in the SM Higgs boson mass-squared term. In models that conserve R-parity [17], each superpartner carries a conserved quantum number that requires superpartners to be produced in pairs and causes the lightest SUSY particle (LSP) to be stable. The stable LSP can serve as a DM candidate.
The signatures targeted in this paper are motivated by models in which gauge-mediated SUSY breaking (GMSB) is responsible for separating the masses of the SUSY particles from those of their SM counterparts. In GMSB models, the gaugino masses are expected to be proportional to the size of their fundamental couplings. This includes the superpartner of the graviton, the gravitino ( G), whose mass is proportional to M SB /M Pl , where M SB represents the scale of the SUSY breaking interactions and M Pl is the Planck scale where gravity is expected to become strong. GMSB permits a significantly lower symmetry-breaking scale than, e.g., gravity mediation, and therefore generically predicts that the G is the LSP [18][19][20], with a mass often much less than 1 GeV. Correspondingly, the next-to-LSP (NLSP) is typically a neutralino, a superposition of the superpartners of the neutral bosons. The details of the quantum numbers of the NLSP play a large part in determining the phenomenology of GMSB models, including the relative frequencies of the Higgs bosons, Z bosons, and photons produced in the NLSP decay.
The scenario of a natural SUSY spectrum with GMSB and R-parity conservation typically manifests as events with multiple jets, at least one photon, and large p miss T , the magnitude of the missing transverse momentum. Depending on the topology, these jets can arise from either light-flavored quarks (u, d, s, c) or b quarks. We study four simplified models [21][22][23][24][25]; example diagrams depicting these models are shown in Fig. 1. Three models involve gluino pair production (prefixed with T5), and one model involves top squark pair production (prefixed with T6). In the T5qqqqHG model, each gluino decays to a pair of light-flavored quarks (qq) and a neutralino ( χ 0 1 ). The T5bbbbZG and T5ttttZG models are similar to T5qqqqHG, except that the each pair of lightflavored quarks is replaced by a pair of bottom quarks (bb) or a pair of top quarks (tt), respectively. In the T5qqqqHG model, the χ 0 1 decays either to an SM Higgs boson and a G or to a photon and a G. The χ 0 1 → H G branching fraction is assumed to be 50%, and the smallest χ 0 1 mass considered is 127 GeV. In the T5bbbbZG and T5ttttZG models, the neutralinos decay to Z G and γ G with equal probability. The T6ttZG model considers top squark pair production, with each top squark decaying into a top quark and a neutralino. The neutralino can then decay with equal probability to a photon and a G or to a Z boson and a G. For the models involving the decay χ 0 1 → Z G, we probe χ 0 1 masses down to 10 GeV. All decays of SUSY particles are assumed to be prompt. In all models, the mass m G is fixed to be 1 GeV, to be consistent with other published results. For the parameter space explored here, the kinematic properties do not depend strongly on the exact value of m G .
The proton-proton (pp) collision data used in this search correspond to an integrated luminosity of 35.9 fb −1 and were collected with the CMS detector during the 2016 run of the CERN LHC [26]. Signal-like events with at least one photon are classified into signal regions depending on the number of jets N jets , the number of tagged bottom quark jets N b-jets , and the p miss T . The expected yields from SM backgrounds are estimated using a combination of simulation and data control regions. We search for gluino or top squark pair production as an excess of observed data events compared to the expected background yields.
Previous searches for R-parity conserving SUSY with photons in the final state performed by the CMS Collaboration are documented in Refs. [27,28]. Similar searches have also been performed by the ATLAS Collaboration [29][30][31]. This work improves on the previous results by identifying jets from b quarks, which can be produced by all of the signal models shown in Fig. 1. We also include additional signal regions that exploit high jet multiplicities for sensitivity to high-mass gluino models, and we rely more on observed data for the background estimations. These improvements enable us to explore targeted signal models that produce b quarks in the final state and are expected to improve sensitivity to the models explored in Refs. [27][28][29][30][31].
In this paper, a description of the CMS detector and simulation used are presented in Sect. 2. The event reconstruction and signal region selections are presented in Sect. 3. The methods used for predicting the SM backgrounds are presented in Sect. 4. Results are given in Sect. 5. The analysis is summarized in Sect. 6.

Detector and simulation
A detailed description of the CMS detector, along with a definition of the coordinate system and pertinent kinematic variables, is given in Ref. [32]. Briefly, a cylindrical superconducting solenoid with an inner diameter of 6 m provides a 3.8 T axial magnetic field. Within the cylindrical volume are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter (HCAL). The tracking detectors cover the pseudorapidity range |η| < 2.5. The ECAL and HCAL, each composed of a barrel and two endcap sections, cover |η| < 3.0. Forward calorimeters extend the coverage to 3.0 < |η| < 5.0. Muons are detected within |η| < 2.4 by gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid. The detector is nearly hermetic, permitting accurate measurements of p miss T . The CMS trigger is described in Ref. [33].
Monte Carlo (MC) simulation is used to design the analysis, to provide input for background estimation methods that use data control regions, and to predict event rates from simplified models. Simulated SM background processes include jets produced through the strong interaction, referred to as quantum chromodynamics (QCD) multijets, tt+jets, W+jets, Z+jets, γ +jets, ttγ , tγ , and Vγ +jets (V = Z, W). The SM background events are generated using the MadGraph5_amc@nlo v2.2.2 or v2.3.3 generator [34][35][36] at leading order (LO) in perturbative QCD, except ttγ and tγ , which are generated at next-to-leading order (NLO). The cross sections used for normalization are computed at NLO or next-to-NLO [34,[37][38][39]. The QCD multijets, diboson (Vγ ), top quark, and vector boson plus jets events are generated with up to two, two, three, and four additional partons in the matrix element calculations, respectively. Any duplication of events between pairs of related processes -QCD multijets and γ +jets; tt+jets and ttγ ; W+jets and Wγ +jets-is removed using generator information.
The NNPDF3.0 [40] LO (NLO) parton distribution functions (PDFs) are used for samples simulated at LO (NLO). Parton showering and hadronization are described using the pythia 8.212 generator [41] with the CUETP8M1 underlying event tune [42]. Partons generated with Mad-Graph5_amc@nlo and pythia that would otherwise be counted twice are removed using the MLM [43] and FxFx [44] matching schemes in LO and NLO samples, respectively.
The detector response to particles produced in the simulated collisions is modeled with the Geant4 [50] detector simulation package for SM processes. Because of the large number of SUSY signals considered, with various gluino, squark, and neutralino masses, the detector response for these processes is simulated with the CMS fast simulation [51,52]. The results from the fast simulation generally agree with the results from the full simulation. Where there is disagreement, corrections are applied, most notably a correction of up to 10% to adjust for differences in the modeling of p miss T .

Event reconstruction and selection
The CMS particle-flow (PF) algorithm [53] aims to reconstruct every particle in each event, using an optimal combination of information from all detector systems. Particle candidates are identified as charged hadrons, neutral hadrons, electrons, photons, or muons. For electron and photon PF candidates, further requirements are applied to the ECAL shower shape and the ratio of associated energies in the ECAL and HCAL [54,55]. Similarly, for muon PF candidates, further requirements are applied to the matching between track segments in the silicon tracker and the muon detectors [56]. These further requirements improve the quality of the reconstruction. Electron and muon candidates are restricted to |η| < 2.5 and < 2.4, respectively. The p miss T is calculated as the magnitude of the negative vector p T sum of all PF candidates.
After all interaction vertices are reconstructed, the primary pp interaction vertex is selected as the vertex with the largest p 2 T sum of all physics objects. The physics objects used in this calculation are produced by a jet-finding algorithm [57,58] applied to all charged-particle tracks associated to the vertex, plus the corresponding p miss T computed from those jets. To mitigate the effect of secondary pp interactions (pileup), charged-particle tracks associated with vertices other than the primary vertex are not considered for jet clustering or calculating object isolation sums.
Jets are reconstructed by clustering PF candidates using the anti-k T jet algorithm [57,58] with a size parameter of 0.4. To eliminate spurious jets, for example those induced by electronics noise, further jet quality criteria [59] are applied. The jet energy response is corrected for the nonlinear response of the detector [60]. There is also a correction to account for the expected contributions of neutral particles from pileup, which cannot be removed based on association with secondary vertices [61]. Jets are required to have p T > 30 GeV and are restricted to be within |η| < 2.4. The combined secondary vertex algorithm (CSVv2) at the medium working point [62] is applied to each jet to determine if it should be identified as a bottom quark jet. The CSVv2 algorithm at the specified working point has a 55% efficiency to correctly identify b jets with p T ≈ 30 GeV. The corresponding misidentification probabilities are 1.6% for gluon and lightflavor quark jets, and 12% for charm quark jets.
Photons with p T > 100 GeV and |η| < 2.4 are used in this analysis, excluding the ECAL transition region with 1.44 < |η| < 1.56. To suppress jets erroneously identified as photons from neutral hadron decays, photon candidates are required to be isolated. An isolation cone of radius ΔR = √ (Δφ) 2 + (Δη) 2 < 0.2 is used, with no dependence on the p T of the photon candidate. Here, φ is the azimuthal angle in radians. The energy measured in the isolation cone is corrected for contributions from pileup [61].
The shower shape and the fractions of hadronic and electromagnetic energy associated with the photon candidate are required to be consistent with expectations from prompt photons. The candidates matched to a track measured by the pixel detector (pixel seed) are rejected because they are likely to result from electrons that produced electromagnetic showers.
Similarly, to suppress jets erroneously identified as leptons and genuine leptons from hadron decays, electron and muon candidates are also subjected to isolation requirements. The isolation variable I is computed from the scalar p T sum of selected charged hadron, neutral hadron, and photon PF candidates, divided by the lepton p T . PF candidates enter the isolation sum if they satisfy R < R I ( p T ). The cone radius R I decreases with lepton p T because the collimation of the decay products of the parent particle of the lepton increases with the Lorentz boost of the parent [63]. The values used are R I = 0.2 for p T < 50 GeV, R I = 10 GeV/ p T for 50 ≤ p T ≤ 200 GeV, and R I = 0.05 for p T > 200 GeV, where = e, μ. As with photons, the expected contributions from pileup are subtracted from the isolation variable. The isolation requirement is We additionally veto events if they contain PF candidates which are identified as an electron, a muon, or a charged hadron, and satisfy an isolation requirement computed using tracks. Isolated hadronic tracks are common in background events with a tau lepton that decays hadronically. The track isolation variable I track is computed for each candidate from the scalar p T sum of selected other charged-particle tracks, divided by the candidate p T . Other charged-particle tracks are selected if they lie within a cone of radius 0.3 around the candidate direction and come from the primary vertex. The isolation variable must satisfy I track < 0.2 for electrons and muons, and I track < 0.1 for charged hadrons. Isolated tracks are required to satisfy |η| < 2.4, and the transverse mass of each isolated track with p miss where Δφ is the difference in φ between p track T and p miss T , is required to be less than 100 GeV. Signal event candidates were recorded by requiring a photon at the trigger level with a requirement p GeV otherwise. These quantities are computed at the trigger level. The efficiency of this trigger, as measured in data, is (98 ± 2)% after applying the selection criteria described below. Additional triggers, requiring the presence of charged leptons, photons, or minimum H T = Σ p jet T , are used to select control samples employed in the evaluation of backgrounds.
Signal-like candidate events must fulfill one of two requirements, based on the trigger criteria described above: p γ T > 100 GeV and H γ T > 800 GeV, or p γ T > 190 GeV and H γ T > 500 GeV. In addition to these requirements, the events should have at least 2 jets and p miss T > 100 GeV. To reduce backgrounds from the SM processes that pro-duce a leptonically decaying W boson, resulting in p miss T from the undetected neutrino, events are rejected if they have any charged light leptons (e, μ) with p T > 10 GeV or any isolated electron, muon, or charged hadron tracks with p T > 5, 5, 10 GeV, respectively. Events from the γ +jets process typically satisfy the above criteria when the energy of a jet is mismeasured, inducing artificial p miss T . To reject these events, the two highest p T jets are both required to have an angular separation from the p miss T direction in the transverse plane, Δφ 1,2 > 0.3. Events with reconstruction failures, detector noise, or beam halo interactions are rejected using dedicated identification requirements [64].
The selected events are divided into 25 exclusive signal regions, also called signal bins, based on p miss T , the number of jets N jets , and the number of b-tagged jets N b-jets . The signal regions can be grouped into 6 categories based on N jets and N b-jets , whose intervals are defined to be N jets : 2-4, 5-6, ≥7; and N b-jets : 0, ≥1. Within each of the 6 categories, events are further distinguished based on 4 exclusive regions, defined as: 200< p miss T <270, 270< p miss T <350, 350< p miss T < 450, and p miss T >450 GeV. In the lowest N jets , N b-jets category, the highest p miss T bin is further subdivided into two intervals: 450 < p miss T < 750 and p miss T > 750 GeV. Events with 100 < p miss T < 200 GeV are used as a control region for estimating SM backgrounds. These categories in N jets , N b-jets , and p miss T were found to provide good sensitivity to the various signal models described above, while minimizing uncertainties in the background predictions.

Background estimation
There are four main mechanisms by which SM processes can produce events with the target signature of a photon, multiple jets, and p miss T . These mechanisms are: (1) the production of a highp T photon along with a W or Z boson that decays leptonically, and either any resulting electron or muon is "lost" (lost-lepton) or any resulting τ lepton decays hadronically (τ h ); (2) the production of a W boson that decays to eν and the electron is misidentified as a photon; (3) the production of a highp T photon in association with a Z boson that decays to neutrinos; and (4) the production of a photon along with a jet that is mismeasured, inducing high p miss T . QCD multijet events with a jet misidentified as a photon and a mismeasured jet do not contribute significantly to the SM background.
The total event yield from each source of background is estimated separately for each of the 25 signal regions. The methods and uncertainties associated with the background predictions are detailed in the following sections.

Lost-lepton and τ h backgrounds
The lost-lepton background arises from events in which the charged lepton from a leptonically decaying W boson, pro-duced directly or from the decay of a top quark, cannot be identified. This can occur because the lepton is out of acceptance, fails the identification requirements, or fails the isolation requirements. For example, in events with highp T top quarks, the top quark decay products will be collimated, forcing the b jet to be closer to the charged lepton. In this case, the lepton is more likely to fail the isolation requirements. This background is estimated by studying control regions in both data and simulation, obtained by requiring both a wellidentified photon and a light lepton (e, μ). For every signal region, there are two lost lepton control regions that have the exact same definition as the signal region except either exactly one electron or exactly one muon is required.
The τ h background arises from events in which a W boson decays to a τ lepton, which subsequently decays to mesons and a neutrino. These hadronic decays of τ leptons occur approximately 65% of the time. Because of lepton universality, the fraction of events with τ h candidates can be estimated from the yield of events containing a single muon, after correcting for the reconstruction differences and for the τ h branching fraction.
The lost-lepton and τ h background predictions rely on an extrapolation between eγ or μγ event yields and single photon event yields. In all control regions where a single light lepton is required, the dominant SM processes that contribute are Wγ and ttγ . Lost-muon and hadronic tau events are estimated using μγ control regions, while lost-electron events are estimated using eγ control regions. In each control region, exactly one electron or muon is required and the isolated track veto for the selected lepton flavor is removed. In order to reduce the effect of signal contamination and to increase the fraction of SM events in the control sample, events are only selected if the m T of the leptonp miss T system is less than 100 GeV. In SM background events with a single lepton and p miss T , the m T of the system is constrained by the mass of the W boson; this is not the case for signal events, because of the presence of gravitinos. All other kinematic variable requirements for each signal region are applied to the corresponding control regions. Transfer factors are derived using simulated Wγ +jets and ttγ processes, which determine the average number of events expected in the signal region for each eγ or μγ event observed in the control region. The Zγ events in which the Z boson decays leptonically have a negligible contribution to the transfer factors. The transfer factors applied to the μγ control regions account for both lost-μ events and τ h events. They are denoted by the symbol T μ,τ and are typically in the range 0.7 < T μ,τ < 1.0. The transfer factors applied to eγ events account for only the lost-e events. They are denoted by the symbol T e and are typically in the range 0.3 < T e < 0.6. The transfer factors are parameterized versus N jets , N b-jets , and p miss T ; however, for p miss T > 150 GeV, T is found to be independent of p miss T . The parameterization of the transfer   Table 1 factors is validated using simulation by treating eγ or μγ events like data and comparing the predicted lost-lepton and τ h event yields to the true simulated event yields in the signal regions. This comparison is shown in Fig. 2. The prediction in each signal region is N pred = Σ i N i T ,i , where = e, μ and i ranges from 1 to n, where n is the number of transfer factors that contribute in a given signal region.
The dominant uncertainty in the lost-lepton background predictions arises from the limited numbers of events in the eγ and μγ control regions. These uncertainties are modeled in the final statistical interpretations as a gamma distribution whose shape parameter is set by the observed number of events and whose scale parameter is the average transfer factor for that bin. Other systematic uncertainties in the determination of the transfer factors include the statistical uncertainty from the limited number of simulated events, which is typically 5-10% but can be as large as 20%, as well as uncertainties in the jet energy corrections, PDFs, renormalization (μ R ) and factorization (μ F ) scales, and simulation correction factors. The uncertainties in μ R and μ F are obtained by varying each value independently by factors of 0.5 and 2.0 [65,66]. Simulation correction factors are used to account for differences between the observed data and modeling of b-tagging efficiencies, b jet misidentification, and lepton reconstruction efficiencies in simulation. One of the largest uncertainties, apart from the statistical uncertainty in the data control regions and the simulation, comes from mismodeling of photons which are collinear with electrons, which has a 12% effect on the lost-lepton prediction.

Misidentified photon background
Events containing the decay W → eν are the primary source of electrons that are erroneously identified as photons. Photon misidentification can occur when a pixel detector seed fails to be associated with the energy deposit in the ECAL. Given a misidentification rate, which relates events with an erroneously identified photon to events with a well-identified electron, the photon background can be estimated from a single-electron (zero-photon) control region. The misidentification rate is estimated in simulation and corrections are derived from observed data to account for any mismodeling in simulation.
The single-electron control regions are defined by the same kinematic requirements as the single-photon signal regions, except that we require no photons and exactly one electron, and we use the momentum of the electron in place of the momentum of the photon for photon-based variables. As explained in the previous section, in addition to all of the signal region selections, events are required to satisfy m T (e, p miss T ) < 100 GeV. To extrapolate from the event yields in the single-electron control regions to the event yields for the misidentified photon background in the signal regions, we derive a misidentification rate f = N γ /N e using a combination of simulation and data. The misidentification rate is determined as a function of the electron p T and the multiplicity Q mult of charged-particle tracks from the primary vertex in a region around the electron candidate. The charged-track multiplicity is computed by counting the number of charged PF candidates (electrons, muons, hadrons) in the jet closest to the electron candidate. If there is no jet within ΔR < 0.3 of the electron candidate, Q mult is set to zero. A typical event in the single-electron control region has a Q mult of 3-4. The electron p T and Q mult dependence of the misidentification rate is derived using simulated W+jets and tt+jets events. The misidentification rate is on average 1-2%, but can be as low as 0.5% for events with high Q mult .
To account for systematic differences between the misidentification rates in data and simulation, we correct the misidentification rate by measuring it in both simulated and observed Drell-Yan (DY) events. Separate corrections are derived for low Q mult (≤1) and high Q mult (≥2). The DY control region is defined by requiring one electron with p T > 40 GeV and another reconstructed particle, either a photon or an oppositely charged electron, with p T > 100 GeV. A further requirement 50 < (m e + e − or m eγ ) < 130 GeV is applied to ensure the particles are consistent with the decay products of a Z boson, and therefore the photon is likely to be a misidentified electron. The misidentifi-cation rate is computed as the ratio N eγ /N e + e − , where N eγ (N e + e − ) is the number of events in the eγ (e + e − ) control region. It is found to be 15-20% higher in data than in simulation.
The prediction of the misidentified-photon background in the signal region is then given by the weighted sum of the observed events in the control region, where the weight is given by the data-corrected misidentification rate for photons. The dominant uncertainty in the prediction is a 14% uncertainty in the data-to-simulation correction factors, followed by the uncertainty in the limited number of events in the simulation at large values of p miss T . The misidentifiedphoton background prediction also includes uncertainties in the modeling of initial-state radiation (ISR) in the simulation, statistical uncertainties from the limited number of events in the data control regions, uncertainties in the pileup modeling, and uncertainties in the trigger efficiency measurement.

Background from Z(νν)γ events
Decays of the Z boson to invisible particles constitute a major background for events with low N jets , low N b-jets , and high p miss T . The Z(νν)γ background is estimated using Z( + − )γ events. The shape of the distribution of p miss T vs. N jets in Z(νν)γ events is modeled in simulation, while the normalization and the purity of the control region are measured in data.
Events in the + − γ control region are required to have exactly two oppositely charged, same-flavor leptons ( = e or μ) and one photon with p T > 100 GeV. The dilepton invariant mass m is required to be consistent with the Z boson mass, 80 < m < 100 GeV. The charged leptons serve as a proxy for neutrinos, so the event-level kinematic variables, such as p miss T , are calculated after removing charged leptons from the event.
The + − γ control region may contain a small fraction of events from processes other than Z( + − )γ , primarily ttγ . We define the purity of the control region as the percentage of events originating from the Z( + − )γ process. The purity is computed in data by measuring the number of events in the corresponding oppositely charged, different-flavor control region, which has a higher proportion of ttγ events. The purity is found to be (97 ± 3)%. A statistically compatible purity is also measured in the oppositely charged, same-flavor control region. In this region, the m distribution is used to extrapolate from the number of events with m far from the Z boson mass to the number of events with m close to it.
The Z(νν)γ predictions from simulation are scaled to the total Z( + − )γ yield observed according to N Z(νν)γ = β R νν/ N Z( + − )γ , where β is the purity of the Z( + − )γ control region and R νν/ is the ratio between the expected number of Z(νν)γ and Z( + − )γ events. The ratio R νν/ , which accounts for lepton reconstruction effects and the rel-ative branching fractions for Z → νν and Z → + − , is determined from simulation.
The primary uncertainty in the Z(νν)γ prediction arises from uncertainties in the p miss T distribution from the simulation. Other uncertainties include statistical uncertainties from the limited number of events in the simulation and uncertainties in the estimation of the control region purity. The p γ Tdependent NLO electroweak corrections [67] are assigned as additional uncertainties to account for any mismodeling of the photon p T in simulation. This uncertainty has a magnitude of 8% for the lowest p miss T bin and rises to 40% for p miss T > 750 GeV.

Background from γ +jets events
The γ +jets background is dominated by events in which a genuine photon is accompanied by an energetic jet with mismeasured p T , resulting in high p miss T . The QCD multijet events with a jet misidentified as a photon and a mismeasured jet contribute to this background at a much smaller rate; these events are measured together with events from the γ +jets process. Most of these events are removed by requiring that the azimuthal angles between the p miss T and each of the two highest p T jets satisfy Δφ 1,2 > 0.3. Inverting this requirement provides a large control region of low-Δφ events that is used to predict the γ +jets background in the signal regions. The ratio of high-Δφ events to low-Δφ events, R h/l , is derived from the lowp miss T sideband (100 < p miss T < 200 GeV). While most of the events in both the low-Δφ and the lowp miss T control regions are γ +jets events, electroweak backgrounds in which p miss T arises from W or Z bosons decaying to one or more neutrinos, like those discussed previously, will contaminate these control regions. The contamination can be significant for high N jets and N b-jets , where tt events are more prevalent. The rates of these events in the control regions are predicted using the same techniques, as discussed in the previous sections.
is derived from simulated γ +jets events in order to account for the dependence of R h/l on p miss T . To test how well the simulation models κ, we use a zero-photon validation region in which the contribution from events containing a mismeasured jet dominates. To be consistent with the trigger used to select the data in this region, these events are also required to have H T > 1000 GeV. Electroweak contamination in the zero-photon validation region is estimated using simulated Vγ +jets (V = Z, W), ttγ , tt+jets, W+jets, and Z(νν)+jets events. The comparison of κ in data and simulation is shown in Fig. 3. The level of disagreement is found to be less than 20%.
Event yields for the γ +jets background are computed from the highp miss T , low-Δφ control regions according to  N γ +jets = κ N low-Δφ R h/l . N low-Δφ is the event yield in the highp miss T , low-Δφ control region after removing contributions from electroweak backgrounds.
Uncertainties in the γ +jets prediction are dominated by the statistical uncertainties either from the limited number of events in the low-Δφ control regions or from the predictions of the electroweak contamination. The <20% disagreement between the κ values in data and simulation in the zerophoton validation region is included as an additional uncertainty. Uncertainties in the b-tagging correction factors are a minor contribution to the uncertainty in the γ +jets prediction.

Results and interpretations
The predicted background and observed yields are shown in Table 1 and Fig. 4. The largest deviation is found in bin 2 (2 ≤ N jets ≤ 4, N b-jets = 0, and 270 < p miss T < 350 GeV), where the background is predicted to be 91 events with 51 events observed. The local significance of this single bin was computed to be around 2 standard deviations below the SM expectation. This calculation does not account for the lookelsewhere effect associated with the use of 25 exclusive signal regions, which is expected to reduce this significance. In general, a large deviation in a single bin is inconsistent with the expected distributions of events from the signal models considered here. The observations in all other bins are consistent with the SM expectations within one standard deviation.
Limits are evaluated for the production cross sections of the signal scenarios discussed in Sect. 1 using a maximum likelihood fit for the SUSY signal strength, the yields of the five classes of background events shown in Fig. 4, and various  nuisance parameters. The SUSY signal strength μ is defined to be the ratio of the observed signal cross section to the predicted cross section. A nuisance parameter refers to a variable not of interest in this search, such as the effect of parton distribution function uncertainties in a background prediction. The nuisance parameters are constrained by observed data in the fit. The uncertainties in the predicted signal yield arise from the uncertainties in renormalization and factorization scales, ISR modeling, jet energy scale, b-tagging efficiency and misidentification rate, corrections to simulation, limited numbers of simulated events, and the integrated luminosity measurement [26]. The largest uncertainty comes from the ISR modeling; it ranges from 4 to 30% depending on the signal region and the signal parameters, taking higher values for regions with large N jets or for signals with Δm ≈ 0. Here, Δm is the difference in mass between the gluino or squark and its decay products, e.g. Δm = m g − (m χ 0 1 + 2m t ) for the T5ttttZG model when on-shell top quarks are produced. The second-largest uncertainty comes from the correction for differences between Geant4 and the fast simulation in p miss T modeling, with a maximum value of 10%. The procedures used to evaluate the systematic uncertainties in the signal predictions in the context of this search are described in Ref. [68].
For the models of gluino pair production considered here, the limits are derived as a function of m g and m χ 0 1 , while for the model of top squark pair production, the limits are a function of m t and m χ 0 1 . The likelihood used for the statistical interpretation models the yield in each of the signal regions as a Poisson distribution, multiplied by constraints which account for the uncertainties in the background predictions and signal yields. For the predictions in which an observed event yield in a control region is scaled, a gamma distribution is used to model the Poisson uncertainty of the observed control region yield. All other uncertainties are modeled as log-normal distributions. The test statistic is q μ = −2 ln L μ /L max , where L max is the maximum likelihood determined by leaving all parameters as free, including the signal strength, and L μ is the maximum likelihood for a fixed value of μ. Limits are determined   Table 1. The lower panel shows the ratio of the observed events to the predicted SM background events. The error bars in the lower panel are the quadrature sum of the statistical uncertainty in the observed data and the systematic uncertainty in the predicted backgrounds before the adjustments based on a maximum likelihood fit to data assuming no signal strength using an approximation of the asymptotic form of the test statistic distribution [69] in conjunction with the CL s criterion [70,71]. Expected upper limits are derived by varying observed yields according to the expectations from the background-only hypothesis.
Using the statistical procedure described above, 95% confidence level (CL) upper limits are computed on the signal cross section for each simplified model and each mass hypothesis. Exclusion limits are defined by comparing observed upper limits to the predicted NLO+NLL signal cross section. The signal cross sections are also varied according to theoretical uncertainties to give a ±1 standard deviation variation on the observed exclusion contour. The 95% CL cross section limits and exclusion contours for the four models considered, T5qqqqHG, T5bbbbZG, T5ttttZG, and T6ttZG, are shown in Fig. 5.
Generally, the limits degrade at both high and low m χ 0 1 . For m χ 0 1 ≈ m g (m t ), the quarks from the decay of gluinos (top squarks) have low p T . Correspondingly, the H γ T , N jets , and N b-jets distributions tend toward lower values, reducing the signal efficiency and causing signal events to populate regions with higher background yields. For small m χ 0 1 , the quarks produced in the decay of gluinos or top squarks have high p T but lower p miss T on average. For all models except T5qqqqHG, when the NLSP mass drops below the mass of the Z boson, the kinematics of the NLSP decay require the Z boson to be far off-shell. As the Z boson mass is forced to be lower, the LSP will carry a larger fraction of the momentum of the NLSP, producing larger p miss T . This causes a slight increase in the sensitivity when the NLSP mass is near the Z boson mass. While a similar effect would happen for the T5qqqqHG model, the simulation used here does not probe the region of parameter space where the Higgs boson would be forced to have a mass far off-shell. Similarly, the limits for top squark production improve slightly at very high m χ 0 1 , when the top quarks become off-shell. In this case, the χ 0 1 carries a larger fraction of the top squark momentum, increasing the p miss T . For moderate m χ 0 1 , gluino masses as large as 2090, 2120, and 1970 GeV are excluded for the T5qqqqHG, T5bbbbZG, and T5ttttZG models, respectively. Top squark masses as large as 1230 GeV are excluded for the T6ttZG model. For small m χ 0 1 , gluino masses as large as 1920, 1950, and 1800 GeV are excluded for the T5qqqqHG, T5bbbbZG, and T5ttttZG models, respectively. Top squark masses as large as 1110 GeV are excluded for the T6ttZG model. There is close agreement between the observed and expected limits.

Summary
A search for gluino and top squark pair production is presented, based on a proton-proton collision dataset at a centerof-mass energy of 13 TeV recorded with the CMS detector in 2016. The data correspond to an integrated luminosity of 35.9 fb −1 . Events are required to have at least one isolated photon with transverse momentum p T > 100 GeV, two jets with p T > 30 GeV and pseudorapidity |η| < 2.4, and missing transverse momentum p miss T > 200 GeV. The data are categorized into 25 exclusive signal regions based on the number of jets, the number of b-tagged jets, and p miss T . Background yields from the standard model processes are predicted using simulation and data control regions. The observed event yields are found to be consistent with expectations from the standard model processes within the uncertainties.
Results are interpreted in the context of simplified models. Four such models are studied, three of which involve gluino pair production and one of which involves top squark pair production. All models assume a gauge-mediated supersymmetry (SUSY) breaking scenario, in which the lightest SUSY particle is a gravitino ( G). We consider scenarios in which the gluino decays to a neutralino χ 0 1 and a pair of light-flavor quarks (T5qqqqHG), bottom quarks (T5bbbbZG), or top quarks (T5ttttZG). In the T5qqqqHG model, the χ 0 1 decays with equal probability either to a photon and a G or to a Higgs boson and a G. In the T5bbbbZG and T5ttttZG models, the χ 0 1 decays with equal probability either to a photon and a G or to a Z boson and a G. In the top squark pair production model (T6ttZG), top squarks decay to a top quark and χ 0 1 , and the χ 0 1 decays with equal probability either to a photon and a G or to a Z boson and a G.   Using the cross sections for SUSY pair production calculated at next-to-leading order plus next-to-leading logarithmic accuracy, we place 95% confidence level lower limits on the gluino mass as large as 2120 GeV, depending on the model and the m χ 0 1 value, and limits on the top squark mass as large as 1230 GeV, depending on the m χ 0 1 value. These results significantly improve upon those from previous searches for SUSY with photons.
Acknowledgements 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. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agen-

Data Availability Statement
This manuscript has no associated data or the data will not be deposited [Authors' comment: Release and preservation of data used by the CMS Collaboration as the basis for publications is guided by the CMS policy as written in its document "CMS data preservation, re-use and open access policy" (https://cms-docb.cern.ch/ cgi-bin/PublicDocDB/RetrieveFile?docid=6032&filename=CMSDat aPolicyV1.2.pdf&version=2).] Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3