CMS-PAS-SUS-17-012

CMS-PAS-SUS-17-012
Search for supersymmetry using events with a photon, a lepton, and missing transverse momentum in pp collisions at $\sqrt{s} =$ 13 TeV
CMS Collaboration
July 2018

Abstract: Results of a search for supersymmetry are presented using events with a photon, an electron or muon, and large missing transverse momentum. The analysis is based on a data sample corresponding to an integrated luminosity of 35.9 fb $^{-1}$ of proton-proton collisions at $\sqrt{s} =$ 13 TeV, produced by the CERN LHC and collected with the CMS detector in 2016. Models of supersymmetry with gauge-mediated supersymmetry breaking yield events with photons in the final state as well as electroweak gauge bosons decaying to leptons. Searches for events with both a photon and a lepton are sensitive probes of these models. No excess of events is observed beyond expectations from standard model processes. The results of the search are interpreted in the context of simplified models inspired by gauge-mediated supersymmetry breaking. These models are used to derive upper limits on the production cross sections of supersymmetric processes and set bounds on masses of supersymmetric particles.
Links: CDS record (PDF) ; CADI line (restricted) ; These preliminary results are superseded in this paper, JHEP 01 (2019) 154. The superseded preliminary plots can be found here.

Figures & Tables	Summary	References	CMS Publications

Figures
png pdf	Figure 1: Feynman diagrams showing the production and decay modes of the signal models T5Wg (left), T6Wg (center), and TChiWg (right) considered in this analysis.
png pdf	Figure 1-a: Feynman diagrams showing the production and decay modes of the signal models T5Wg (left), T6Wg (center), and TChiWg (right) considered in this analysis.
png pdf	Figure 1-b: Feynman diagrams showing the production and decay modes of the signal models T5Wg (left), T6Wg (center), and TChiWg (right) considered in this analysis.
png pdf	Figure 1-c: Feynman diagrams showing the production and decay modes of the signal models T5Wg (left), T6Wg (center), and TChiWg (right) considered in this analysis.
png pdf	Figure 2: Verification of the e-to- ${\gamma}$ misidentification estimation method using simulated data. The predicted ${{p_{\mathrm {T}}} ^\text {miss}}$ distribution for events with misidentified photons in the e $\gamma$ (left) and $\mu {\gamma}$ (right) channel from direct simulation (points) and from using the proxy sample estimation method from data (histograms). The vertical bars on the points show the statistical uncertainty in the simulation, while the horizontal bars give the bin widths. The lower panels show the ratio of the predictions from direct simulation to those from the proxy sample method. The vertical bars on the points show again the statistical uncertainty, and the hatched areas give the quadrature sum of the statistical and systematic uncertainties in the simulated background.
png pdf	Figure 2-a: Verification of the e-to- ${\gamma}$ misidentification estimation method using simulated data. The predicted ${{p_{\mathrm {T}}} ^\text {miss}}$ distribution for events with misidentified photons in the e $\gamma$ (left) and $\mu {\gamma}$ (right) channel from direct simulation (points) and from using the proxy sample estimation method from data (histograms). The vertical bars on the points show the statistical uncertainty in the simulation, while the horizontal bars give the bin widths. The lower panels show the ratio of the predictions from direct simulation to those from the proxy sample method. The vertical bars on the points show again the statistical uncertainty, and the hatched areas give the quadrature sum of the statistical and systematic uncertainties in the simulated background.
png pdf	Figure 2-b: Verification of the e-to- ${\gamma}$ misidentification estimation method using simulated data. The predicted ${{p_{\mathrm {T}}} ^\text {miss}}$ distribution for events with misidentified photons in the e $\gamma$ (left) and $\mu {\gamma}$ (right) channel from direct simulation (points) and from using the proxy sample estimation method from data (histograms). The vertical bars on the points show the statistical uncertainty in the simulation, while the horizontal bars give the bin widths. The lower panels show the ratio of the predictions from direct simulation to those from the proxy sample method. The vertical bars on the points show again the statistical uncertainty, and the hatched areas give the quadrature sum of the statistical and systematic uncertainties in the simulated background.
png pdf	Figure 3: The post-fit $\Delta \phi (\ell, {{p_{\mathrm {T}}} ^\text {miss}})$ distributions for the $V {\gamma}$ (dashed-green) and misidentified-lepton (solid-red) backgrounds. The black points show the data in the 40 $< {{p_{\mathrm {T}}} ^\text {miss}} <$ 70 GeV control region with the fit result overlaid for the e $\gamma$ (left) and $\mu {\gamma}$ (right) channels, while the grey-blue distribution gives the fit result and the hatched area indicates the fit uncertainty. The vertical bars on the points represent the statistical uncertainty. The lower panels show the ratio of the fit result to the data. The vertical bars on the points show the statistical uncertainty in the data, and the hatched areas give the fit uncertainty.
png pdf	Figure 3-a: The post-fit $\Delta \phi (\ell, {{p_{\mathrm {T}}} ^\text {miss}})$ distributions for the $V {\gamma}$ (dashed-green) and misidentified-lepton (solid-red) backgrounds. The black points show the data in the 40 $< {{p_{\mathrm {T}}} ^\text {miss}} <$ 70 GeV control region with the fit result overlaid for the e $\gamma$ (left) and $\mu {\gamma}$ (right) channels, while the grey-blue distribution gives the fit result and the hatched area indicates the fit uncertainty. The vertical bars on the points represent the statistical uncertainty. The lower panels show the ratio of the fit result to the data. The vertical bars on the points show the statistical uncertainty in the data, and the hatched areas give the fit uncertainty.
png pdf	Figure 3-b: The post-fit $\Delta \phi (\ell, {{p_{\mathrm {T}}} ^\text {miss}})$ distributions for the $V {\gamma}$ (dashed-green) and misidentified-lepton (solid-red) backgrounds. The black points show the data in the 40 $< {{p_{\mathrm {T}}} ^\text {miss}} <$ 70 GeV control region with the fit result overlaid for the e $\gamma$ (left) and $\mu {\gamma}$ (right) channels, while the grey-blue distribution gives the fit result and the hatched area indicates the fit uncertainty. The vertical bars on the points represent the statistical uncertainty. The lower panels show the ratio of the fit result to the data. The vertical bars on the points show the statistical uncertainty in the data, and the hatched areas give the fit uncertainty.
png pdf	Figure 4: Distributions of ${{p_{\mathrm {T}}} ^\text {miss}}$ (a, b), ${p_{\mathrm {T}}} ^{{\gamma}}$ (c, d), and ${H_{\mathrm {T}}}$ (e, f) from data (points) and simulated SM predictions (stacked histograms) for the e $\gamma$ (left) and $\mu {\gamma}$ (right) channels. Simulated signal distributions from from the TChiWg model (dotted) with $M_{\tilde{\chi}^0 / \tilde{\chi}^\pm} =$ 800 GeV and the T5Wg model (solid) with $M_{{\tilde{g}}} =$ 1700 GeV are overlaid. The ${{p_{\mathrm {T}}} ^\text {miss}}$ distribution includes all events with ${M_{\mathrm {T}}} >$ 100 GeV, while the ${p_{\mathrm {T}}} ^{{\gamma}}$ and ${H_{\mathrm {T}}}$ distributions only include events with ${M_{\mathrm {T}}} >$ 100 GeV and ${{p_{\mathrm {T}}} ^\text {miss}} >$ 120 GeV. The vertical bars on the points give the statistical uncertainty in the data. The horizontal bars show the bin widths. The lower panels display the ratio of the data to the total background prediction. The vertical bars on the points show the statistical uncertainty in the data, and the hatched areas give the quadrature sum of the statistical and systematic uncertainties in the simulated background.
png pdf	Figure 4-a: Distributions of ${{p_{\mathrm {T}}} ^\text {miss}}$ (a, b), ${p_{\mathrm {T}}} ^{{\gamma}}$ (c, d), and ${H_{\mathrm {T}}}$ (e, f) from data (points) and simulated SM predictions (stacked histograms) for the e $\gamma$ (left) and $\mu {\gamma}$ (right) channels. Simulated signal distributions from from the TChiWg model (dotted) with $M_{\tilde{\chi}^0 / \tilde{\chi}^\pm} =$ 800 GeV and the T5Wg model (solid) with $M_{{\tilde{g}}} =$ 1700 GeV are overlaid. The ${{p_{\mathrm {T}}} ^\text {miss}}$ distribution includes all events with ${M_{\mathrm {T}}} >$ 100 GeV, while the ${p_{\mathrm {T}}} ^{{\gamma}}$ and ${H_{\mathrm {T}}}$ distributions only include events with ${M_{\mathrm {T}}} >$ 100 GeV and ${{p_{\mathrm {T}}} ^\text {miss}} >$ 120 GeV. The vertical bars on the points give the statistical uncertainty in the data. The horizontal bars show the bin widths. The lower panels display the ratio of the data to the total background prediction. The vertical bars on the points show the statistical uncertainty in the data, and the hatched areas give the quadrature sum of the statistical and systematic uncertainties in the simulated background.
png pdf	Figure 4-b: Distributions of ${{p_{\mathrm {T}}} ^\text {miss}}$ (a, b), ${p_{\mathrm {T}}} ^{{\gamma}}$ (c, d), and ${H_{\mathrm {T}}}$ (e, f) from data (points) and simulated SM predictions (stacked histograms) for the e $\gamma$ (left) and $\mu {\gamma}$ (right) channels. Simulated signal distributions from from the TChiWg model (dotted) with $M_{\tilde{\chi}^0 / \tilde{\chi}^\pm} =$ 800 GeV and the T5Wg model (solid) with $M_{{\tilde{g}}} =$ 1700 GeV are overlaid. The ${{p_{\mathrm {T}}} ^\text {miss}}$ distribution includes all events with ${M_{\mathrm {T}}} >$ 100 GeV, while the ${p_{\mathrm {T}}} ^{{\gamma}}$ and ${H_{\mathrm {T}}}$ distributions only include events with ${M_{\mathrm {T}}} >$ 100 GeV and ${{p_{\mathrm {T}}} ^\text {miss}} >$ 120 GeV. The vertical bars on the points give the statistical uncertainty in the data. The horizontal bars show the bin widths. The lower panels display the ratio of the data to the total background prediction. The vertical bars on the points show the statistical uncertainty in the data, and the hatched areas give the quadrature sum of the statistical and systematic uncertainties in the simulated background.
png pdf	Figure 4-c: Distributions of ${{p_{\mathrm {T}}} ^\text {miss}}$ (a, b), ${p_{\mathrm {T}}} ^{{\gamma}}$ (c, d), and ${H_{\mathrm {T}}}$ (e, f) from data (points) and simulated SM predictions (stacked histograms) for the e $\gamma$ (left) and $\mu {\gamma}$ (right) channels. Simulated signal distributions from from the TChiWg model (dotted) with $M_{\tilde{\chi}^0 / \tilde{\chi}^\pm} =$ 800 GeV and the T5Wg model (solid) with $M_{{\tilde{g}}} =$ 1700 GeV are overlaid. The ${{p_{\mathrm {T}}} ^\text {miss}}$ distribution includes all events with ${M_{\mathrm {T}}} >$ 100 GeV, while the ${p_{\mathrm {T}}} ^{{\gamma}}$ and ${H_{\mathrm {T}}}$ distributions only include events with ${M_{\mathrm {T}}} >$ 100 GeV and ${{p_{\mathrm {T}}} ^\text {miss}} >$ 120 GeV. The vertical bars on the points give the statistical uncertainty in the data. The horizontal bars show the bin widths. The lower panels display the ratio of the data to the total background prediction. The vertical bars on the points show the statistical uncertainty in the data, and the hatched areas give the quadrature sum of the statistical and systematic uncertainties in the simulated background.
png pdf	Figure 4-d: Distributions of ${{p_{\mathrm {T}}} ^\text {miss}}$ (a, b), ${p_{\mathrm {T}}} ^{{\gamma}}$ (c, d), and ${H_{\mathrm {T}}}$ (e, f) from data (points) and simulated SM predictions (stacked histograms) for the e $\gamma$ (left) and $\mu {\gamma}$ (right) channels. Simulated signal distributions from from the TChiWg model (dotted) with $M_{\tilde{\chi}^0 / \tilde{\chi}^\pm} =$ 800 GeV and the T5Wg model (solid) with $M_{{\tilde{g}}} =$ 1700 GeV are overlaid. The ${{p_{\mathrm {T}}} ^\text {miss}}$ distribution includes all events with ${M_{\mathrm {T}}} >$ 100 GeV, while the ${p_{\mathrm {T}}} ^{{\gamma}}$ and ${H_{\mathrm {T}}}$ distributions only include events with ${M_{\mathrm {T}}} >$ 100 GeV and ${{p_{\mathrm {T}}} ^\text {miss}} >$ 120 GeV. The vertical bars on the points give the statistical uncertainty in the data. The horizontal bars show the bin widths. The lower panels display the ratio of the data to the total background prediction. The vertical bars on the points show the statistical uncertainty in the data, and the hatched areas give the quadrature sum of the statistical and systematic uncertainties in the simulated background.
png pdf	Figure 4-e: Distributions of ${{p_{\mathrm {T}}} ^\text {miss}}$ (a, b), ${p_{\mathrm {T}}} ^{{\gamma}}$ (c, d), and ${H_{\mathrm {T}}}$ (e, f) from data (points) and simulated SM predictions (stacked histograms) for the e $\gamma$ (left) and $\mu {\gamma}$ (right) channels. Simulated signal distributions from from the TChiWg model (dotted) with $M_{\tilde{\chi}^0 / \tilde{\chi}^\pm} =$ 800 GeV and the T5Wg model (solid) with $M_{{\tilde{g}}} =$ 1700 GeV are overlaid. The ${{p_{\mathrm {T}}} ^\text {miss}}$ distribution includes all events with ${M_{\mathrm {T}}} >$ 100 GeV, while the ${p_{\mathrm {T}}} ^{{\gamma}}$ and ${H_{\mathrm {T}}}$ distributions only include events with ${M_{\mathrm {T}}} >$ 100 GeV and ${{p_{\mathrm {T}}} ^\text {miss}} >$ 120 GeV. The vertical bars on the points give the statistical uncertainty in the data. The horizontal bars show the bin widths. The lower panels display the ratio of the data to the total background prediction. The vertical bars on the points show the statistical uncertainty in the data, and the hatched areas give the quadrature sum of the statistical and systematic uncertainties in the simulated background.
png pdf	Figure 4-f: Distributions of ${{p_{\mathrm {T}}} ^\text {miss}}$ (a, b), ${p_{\mathrm {T}}} ^{{\gamma}}$ (c, d), and ${H_{\mathrm {T}}}$ (e, f) from data (points) and simulated SM predictions (stacked histograms) for the e $\gamma$ (left) and $\mu {\gamma}$ (right) channels. Simulated signal distributions from from the TChiWg model (dotted) with $M_{\tilde{\chi}^0 / \tilde{\chi}^\pm} =$ 800 GeV and the T5Wg model (solid) with $M_{{\tilde{g}}} =$ 1700 GeV are overlaid. The ${{p_{\mathrm {T}}} ^\text {miss}}$ distribution includes all events with ${M_{\mathrm {T}}} >$ 100 GeV, while the ${p_{\mathrm {T}}} ^{{\gamma}}$ and ${H_{\mathrm {T}}}$ distributions only include events with ${M_{\mathrm {T}}} >$ 100 GeV and ${{p_{\mathrm {T}}} ^\text {miss}} >$ 120 GeV. The vertical bars on the points give the statistical uncertainty in the data. The horizontal bars show the bin widths. The lower panels display the ratio of the data to the total background prediction. The vertical bars on the points show the statistical uncertainty in the data, and the hatched areas give the quadrature sum of the statistical and systematic uncertainties in the simulated background.
png pdf	Figure 5: The number of signal events from data (points) and from simulated background (stacked histograms) for the 18 search bins in ${{p_{\mathrm {T}}} ^\text {miss}}$ , ${H_{\mathrm {T}}}$ , and ${p_{\mathrm {T}}} ^{{\gamma}}$ in the $\mu + {\gamma}$ channel (left) and the $e+ {\gamma}$ channel (right). For each ${{p_{\mathrm {T}}} ^\text {miss}}$ range, the first, second, and last bins correspond to the ${H_{\mathrm {T}}}$ ranges [0,100], [100,400] and $>$ 400 GeV, respectively. The lower panel displays the ratio of the data to the background predictions. The vertical bars on the points show the statistical uncertainty in the data, and the hatched areas give the quadrature sum of the statistical and systematic uncertainties in the simulated background.
png pdf	Figure 6: The observed and expected 95% CL upper limits on the production cross sections for the TChiWg simplified model, together with the theoretical cross sections. The inner (darker) band and outer (lighter) band around the expected upper limits indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The band around the theoretical cross section gives the $\pm$ 1 standard deviation uncertainty in the cross section.
png pdf	Figure 7: The observed and expected 95% CL exclusion contours for $M_{{\tilde{g}}}/ {\mathrm {\tilde{q}}}$ versus $M_{\tilde{\chi}}$ (regions to the left of the curves are excluded) and the 95% CL upper limits on the pair production cross sections for the (a) T5Wg and (b) T6Wg simplified models. The upper limits on the cross section assume a 50% branching fraction for $\tilde{g}(\tilde{q})\rightarrow \tilde{\chi}^0/\tilde{\chi}^\pm \mathrm{q\bar{q}}$ . The bands around the observed and expected exclusion contours indicate the $\pm$ 1 standard deviation range when including the experimental and theoretical uncertainties, respectively.
png pdf	Figure 7-a: The observed and expected 95% CL exclusion contours for $M_{{\tilde{g}}}/ {\mathrm {\tilde{q}}}$ versus $M_{\tilde{\chi}}$ (regions to the left of the curves are excluded) and the 95% CL upper limits on the pair production cross sections for the (a) T5Wg and (b) T6Wg simplified models. The upper limits on the cross section assume a 50% branching fraction for $\tilde{g}(\tilde{q})\rightarrow \tilde{\chi}^0/\tilde{\chi}^\pm \mathrm{q\bar{q}}$ . The bands around the observed and expected exclusion contours indicate the $\pm$ 1 standard deviation range when including the experimental and theoretical uncertainties, respectively.
png pdf	Figure 7-b: The observed and expected 95% CL exclusion contours for $M_{{\tilde{g}}}/ {\mathrm {\tilde{q}}}$ versus $M_{\tilde{\chi}}$ (regions to the left of the curves are excluded) and the 95% CL upper limits on the pair production cross sections for the (a) T5Wg and (b) T6Wg simplified models. The upper limits on the cross section assume a 50% branching fraction for $\tilde{g}(\tilde{q})\rightarrow \tilde{\chi}^0/\tilde{\chi}^\pm \mathrm{q\bar{q}}$ . The bands around the observed and expected exclusion contours indicate the $\pm$ 1 standard deviation range when including the experimental and theoretical uncertainties, respectively.

Tables
png pdf	Table 1: Summary of the systematic uncertainties affecting the SUSY signals and SM background estimates given in percent.

Summary

A search for supersymmetry with general gauge mediation in events with a photon, an electron or muon, and large missing transverse momentum is presented. This analysis is based on a sample of proton-proton collisions at

$\sqrt{s} =$ 13 TeV corresponding to an integrated luminosity of 35.9 fb

$^{-1}$ recorded by the CMS experiment in 2016. The data are examined in bins of the photon transverse energy, the magnitude of the missing transverse momentum, and

${H_{\mathrm{T}}}$ , the scalar sum of jet energies. The standard model background is evaluated primarily using control samples in the data, with simulation used to evaluate backgrounds from electroweak processes. The data is found to agree with the standard model expectation, without significant excess in the search region. The results of the search are interpreted as a 95% confidence level upper limits on the production cross section of supersymmetry particles in the context of simplified models motivated by gauge-mediated supersymmetry breaking. The TChiWg simplified model, based on direct electroweak production of a neutralino and chargino, is excluded for next-to-lightest supersymmetric particle masses below 900 GeV. For strong production models, such as the T5Wg simplified model of gluino pair production and the T6Wg model of squark pair production, this search excludes gluinos (squarks) with masses of up to 1700 (1400) GeV in the T5Wg (T6Wg) scenarios.

References
1	M. Dine and W. Fischler	A phenomenological model of particle physics based on supersymmetry	PLB 110 (1982) 227
2	L. Alvarez-Gaume, M. Claudson, and M. B. Wise	Low-energy supersymmetry	NPB 207 (1982) 96
3	C. R. Nappi and B. A. Ovrut	Supersymmetric extension of the SU(3)xSU(2)xU(1) model	PLB 113 (1982) 175
4	M. Dine and A. E. Nelson	Dynamical supersymmetry breaking at low-energies	PRD 48 (1993) 1277	hep-ph/9303230
5	M. Dine, A. E. Nelson, and Y. Shirman	Low energy dynamical supersymmetry breaking simplified	PRD 51 (1995) 1362	hep-ph/9408384
6	M. Dine, A. E. Nelson, Y. Nir, and Y. Shirman	New tools for low-energy dynamical supersymmetry breaking	PRD 53 (1996) 2658	hep-ph/9507378
7	G. R. Farrar and P. Fayet	Phenomenology of the production, decay, and detection of new hadronic states associated with supersymmetry	PLB 76 (1978) 575
8	S. Dimopoulos, G. F. Giudice, and A. Pomarol	Dark matter in theories of gauge mediated supersymmetry breaking	PLB 389 (1996) 37	hep-ph/9607225
9	S. P. Martin	Generalized messengers of supersymmetry breaking and the sparticle mass spectrum	PRD 55 (1997) 3177	hep-ph/9608224
10	E. Poppitz and S. P. Trivedi	Some remarks on gauge mediated supersymmetry breaking	PLB 401 (1997) 38	hep-ph/9703246
11	P. Meade, N. Seiberg, and D. Shih	General gauge mediation	Prog. Theor. Phys. Suppl. 177 (2009) 143	0801.3278
12	M. Buican, P. Meade, N. Seiberg, and D. Shih	Exploring general gauge mediation	JHEP 03 (2009) 016	0812.3668
13	S. Abel, M. J. Dolan, J. Jaeckel, and V. V. Khoze	Phenomenology of pure general gauge mediation	JHEP 12 (2009) 001	0910.2674
14	L. M. Carpenter, M. Dine, G. Festuccia, and J. D. Mason	Implementing general gauge mediation	PRD 79 (2009) 035002	0805.2944
15	T. T. Dumitrescu, Z. Komargodski, N. Seiberg, and D. Shih	General messenger gauge mediation	JHEP 05 (2010) 096	1003.2661
16	ATLAS Collaboration	Search for photonic signatures of gauge-mediated supersymmetry in 8 TeV pp collisions with the ATLAS detector	PRD 92 (2015) 072001	1507.05493
17	CMS Collaboration	Search for supersymmetry in events with a lepton, a photon, and large missing transverse energy in pp collisions at $\sqrt{s} =$ 7 TeV	JHEP 06 (2011) 093	CMS-SUS-11-002 1105.3152
18	CMS Collaboration	Search for supersymmetry with photons in $pp$ collisions at $\sqrt{s} =$ 8 TeV	PRD 92 (2015) 072006	CMS-SUS-14-004 1507.02898
19	CMS Collaboration	Search for supersymmetry in events with a photon, a lepton, and missing transverse momentum in pp collisions at $\sqrt{s} =$ 8 TeV	PLB 757 (2016) 6	CMS-SUS-14-013 1508.01218
20	CMS Collaboration	Interpretation of searches for supersymmetry with simplified models	PRD 88 (2013) 052017	CMS-SUS-11-016 1301.2175
21	CMS Collaboration	Performance of photon reconstruction and identification with the CMS detector in proton-proton collisions at $\sqrt{s} =$ 8 TeV	JINST 10 (2015) P08010	CMS-EGM-14-001 1502.02702
22	CMS Collaboration	Performance of electron reconstruction and selection with the CMS detector in proton-proton collisions at $\sqrt{s} =$ 8 TeV	JINST 10 (2015) P06005	CMS-EGM-13-001 1502.02701
23	CMS Collaboration	Performance of CMS muon reconstruction in pp collision events at $\sqrt{s} =$ 7 TeV	JINST 7 (2012) P10002	CMS-MUO-10-004 1206.4071
24	CMS Collaboration	The CMS experiment at the CERN LHC	JINST 3 (2008) S08004	CMS-00-001
25	CMS Collaboration	Particle-flow reconstruction and global event description with the CMS detector	JINST 12 (2017) P10003	CMS-PRF-14-001 1706.04965
26	CMS Collaboration	Performance of the CMS muon detector and muon reconstruction with proton-proton collisions at $\sqrt{s} =$ 13 TeV		CMS-MUO-16-001 1804.04528
27	M. Cacciari, G. P. Salam, and G. Soyez	The anti- $k_t$ jet clustering algorithm	JHEP 04 (2008) 063	0802.1189
28	M. Cacciari, G. P. Salam, and G. Soyez	FastJet user manual	EPJC 72 (2012) 1896	1111.6097
29	M. Cacciari and G. P. Salam	Pileup subtraction using jet areas	PLB 659 (2008) 119	0707.1378
30	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 (2014) 079	1405.0301
31	T. Sjostrand et al.	An introduction to PYTHIA 8.2	CPC 191 (2015) 159	1410.3012
32	CMS Collaboration	Event generator tunes obtained from underlying event and multiparton scattering measurements	EPJC 76 (2016) 155	CMS-GEN-14-001 1512.00815
33	G. Bozzi et al.	Production of Drell--Yan lepton pairs in hadron collisions: Transverse-momentum resummation at next-to-next-to-leading logarithmic accuracy	PLB 696 (2011) 207	1007.2351
34	GEANT4 Collaboration	GEANT4---a simulation toolkit	NIMA 506 (2003) 250
35	CMS Collaboration	Search for top-squark pair production in the single-lepton final state in pp collisions at $\sqrt{s} =$ 8 TeV	EPJC 73 (2013) 2677	CMS-SUS-13-011 1308.1586
36	CMS Collaboration	The CMS trigger system	JINST 12 (2017) P01020	CMS-TRG-12-001 1609.02366
37	Particle Data Group, C. Patrignani et al.	Review of particle physics	CPC 40 (2016) 100001
38	CMS Collaboration	Measurement of the inclusive W and Z production cross sections in pp collisions at $\sqrt{s} =$ 7 TeV	JHEP 10 (2011) 132	CMS-EWK-10-005 1107.4789
39	CMS Collaboration	Measurement of the semileptonic $\mathrm{t\bar{t}}$ + $\gamma$ production cross section in pp collisions at $\sqrt{s} =$ 8 TeV	JHEP 10 (2017) 006	CMS-TOP-14-008 1706.08128
40	CMS Collaboration	Measurement of the $t\bar{t}$ production cross section using events in the e $\mu$ final state in pp collisions at $\sqrt{s} =$ 13 TeV	EPJC 77 (2017) 172	CMS-TOP-16-005 1611.04040
41	CMS Collaboration	The fast simulation of the CMS detector at LHC	J. Phys. Conf. Ser. 331 (2011) 032049
42	T. Junk	Confidence level computation for combining searches with small statistics	NIMA 434 (1999) 435	hep-ex/9902006
43	A. L. Read	Presentation of search results: The CLs technique	JPG 28 (2002) 2693