Expected limit contour for GGM scenario

Observed limit contour for GGM scenario

Expected limit contour for GGM scenario

Observed limit contour for GGM scenario

Expected limit contour for Neutralino BF scenario

Observed limit contour for Neutralino BF scenario

Expected limit contour for Chargino BF scenario

Observed limit contour for Chargino BF scenario

Expected limit contour for T5Wg

Observed limit contour for T5Wg

Expected limit contour for Gluino BF scenario

Observed limit contour for Gluino BF scenario

Covariance matrix for all bins (additional material)

Correlation matrix for all bins (additional material)

The observed and expected 95% CL upper limits on model-independent visible cross-sections, along with the observed $p0$ values, for the stable signal region, as a function of different mass windows, for which the lower bound is shown. The upper boundary on the mass window is 5 TeV for all windows.

The observed and expected 95% CL upper limits on model-independent visible cross-sections, along with the observed $p0$ values, for the metastable signal region, as a function of different mass windows, for which the lower bound is shown. The upper boundary on the mass window is 5 TeV for all windows.

For each gluino lifetime and mass in the signal samples, the lower boundary of the mass window in which at least $70\%$ of the reconstructed signal appears. The upper boundary for all mass windows is 5 TeV.

Acceptance and efficiency for a representative set of pair-produced gluino signal samples. The mass of the gluino ($m(\tilde{g})$), its lifetime ($\tau(\tilde{g})$) and the mass of the neutralino ($m(\tilde{\chi}^{0})$) are given in the first three columns. The Pythia 6.4.27 signal samples shown in this table are not reweighted to match the transverse momentum of the gluino-gluino system as simulated by MadGraph5_aMC@NLO. The acceptance is defined as the fraction of events passing a loose set of fiducial requirements. The full simulation efficiency (Full sim. $\epsilon$) is defined as the ratio of the number of reconstructed events, as expected by the full ATLAS simulation, and the number of events passing the fiducial requirements. The parameterised simulation efficiency (Param. sim. $\epsilon$) is defined as the ratio of the number of events estimated using a set of parametrised efficiencies (see auxiliary Figures 9,10,11,12) and the number of events passing the fiducial requirements alone.

The reconstructed candidate track mass distributions for observed data, predicted background, and the expected contribution from two signal models in the metastable R-hadron signal region. The yellow band around the background estimation includes both the statistical and systematic uncertainties.

The reconstructed candidate track mass distributions for observed data, predicted background, and the expected contribution from two signal models in the stable R-hadron signal region. The yellow band around the background estimation includes both the statistical and systematic uncertainties.

The 95% CL upper limit on the cross-section as a function of mass for gluinos with lifetime $\tau = 10$ ns decaying into $q\bar{q}$ and a 100 GeV neutralino, with the observed limit shown as a solid black line. The predicted production cross-section values are shown in purple along with their uncertainty. The expected upper limit in the case of only background is shown by the dashed black line, with a green $\pm 1\sigma$ and a yellow $\pm 2\sigma$ band.

The 95% CL upper limit on the cross-section as a function of mass for stable gluino $R$-hadrons, with the observed limit shown as a solid black line. The predicted production cross-section values are shown in purple along with their uncertainty. The expected upper limit in the case of only background is shown by the dashed black line, with a green $\pm 1\sigma$ and a yellow $\pm 2\sigma$ band.

Observed 95% lower limits on the gluino mass in the gluino lifetime--mass plane. The excluded area is to the left of the curves.

Expected 95% lower limits on the gluino mass in the gluino lifetime--mass plane. The excluded area is to the left of the curves.

The 95% CL upper limit on the cross-section as a function of mass for gluinos with lifetime $\tau = 1$ ns decaying into $q\bar{q}$ and a 100 GeV neutralino, with the observed limit shown as a solid black line. The predicted production cross-section values are shown in purple along with their uncertainty. The expected upper limit in the case of only background is shown by the dashed black line, with a green $\pm 1\sigma$ and a yellow $\pm 2\sigma$ band.

The 95% CL upper limit on the cross-section as a function of mass for gluinos with lifetime $\tau = 3$ ns decaying into $q\bar{q}$ and a 100 GeV neutralino, with the observed limit shown as a solid black line. The predicted production cross-section values are shown in purple along with their uncertainty. The expected upper limit in the case of only background is shown by the dashed black line, with a green $\pm 1\sigma$ and a yellow $\pm 2\sigma$ band.

The 95% CL upper limit on the cross-section as a function of mass for gluinos with lifetime $\tau = 30$ ns decaying into $q\bar{q}$ and a 100 GeV neutralino, with the observed limit shown as a solid black line. The predicted production cross-section values are shown in purple along with their uncertainty. The expected upper limit in the case of only background is shown by the dashed black line, with a green $\pm 1\sigma$ and a yellow $\pm 2\sigma$ band.

The 95% CL upper limit on the cross-section as a function of mass for gluinos with lifetime $\tau = 50$ ns decaying into $q\bar{q}$ and a 100 GeV neutralino, with the observed limit shown as a solid black line. The predicted production cross-section values are shown in purple along with their uncertainty. The expected upper limit in the case of only background is shown by the dashed black line, with a green $\pm 1\sigma$ and a yellow $\pm 2\sigma$ band.

The relationship between generated and reconstructed mass for gluino $R$-hadrons. Above 1500 GeV, the reconstructed mass falls below the generated mass due to bias in the reconstructed momentum. The uncertainty on the reconstructed mass is dominated by momentum uncertainty. The black dots represent the reconstructed mass computed as the most probable value of a Gaussian fit function, with the error bars showing its statistical uncertainty, while the orange band is the full-width at half maximum of the reconstructed mass distribution.

The parameterised efficiency for events to pass metastable event selections (including trigger, E$_{T}^{miss}$, and event cleaning requirements) as a function of the true E$_{T}^{miss}$ in the system, which is calculated at generator level. Event-level efficiencies are evaluated for events which have at least true E$_{T}^{miss} > 50$ GeV. The metastable event efficiencies are evaluated for different radial regions depending on the smallest radial distance, R, at which an R-hadron decays in the detector.

The parameterised efficiency for events to pass metastable event selections (including trigger, E$_{T}^{miss}$, and event cleaning requirements) as a function of the true E$_{T}^{miss}$ in the system, which is calculated at generator level. Event-level efficiencies are evaluated for events which have at least true E$_{T}^{miss} > 50$ GeV. The stable event efficiencies are evaluated for samples in which no R-hadron decays within the detector.

The parameterised efficiency for particles to pass full track selections in the metastable signal region, as function of the particle’s $\beta$, in different bins of transverse momentum, $p_{T}$, and for different radial decay positions of the particle. The efficiency is evaluated for particles which pass a loose set of fiducial requirements at generator level.

The parameterised efficiency for particles to pass full track selections in the metastable signal region, as function of the particle’s $\beta$, in different bins of transverse momentum, $p_{T}$, and for different radial decay positions of the particle. The efficiency is evaluated for particles which pass a loose set of fiducial requirements at generator level.

The parameterised efficiency for particles to pass full track selections in the metastable signal region, as function of the particle’s $\beta$, in different bins of transverse momentum, $p_{T}$, and for different radial decay positions of the particle. The efficiency is evaluated for particles which pass a loose set of fiducial requirements at generator level.

The parameterised efficiency for particles to pass full track selections in the metastable signal region, as function of the particle’s $\beta$, in different bins of transverse momentum, $p_{T}$, and for different radial decay positions of the particle. The efficiency is evaluated for particles which pass a loose set of fiducial requirements at generator level.

The parameterised efficiency for particles to pass full track selections in the metastable signal region, as function of the particle’s $\beta$, in different bins of transverse momentum, $p_{T}$, and for different radial decay positions of the particle. The efficiency is evaluated for particles which pass a loose set of fiducial requirements at generator level.

The parameterised efficiency for particles to pass full track selections in the metastable signal region, as function of the particle’s $\beta$, in different bins of transverse momentum, $p_{T}$. The stable efficiency is evaluated for samples which do not decay within the detector. The efficiency is evaluated for particles which pass a loose set of fiducial requirements at generator level.

The $m_{\mathrm{eff}}$ distribution for events passing the signal region requirements except the $m_{\mathrm{eff}}$ requirement in SR2. Distributions for data, the estimated SM backgrounds, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. Both the statistical and systematic uncertainties in the SM background are included in the shaded band. The red arrows indicate the $m_{\mathrm{eff}}$ selections in the signal region.

Expected 95% CL exclusion limits on wino $W/Z$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{12k} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on wino $W/Z$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{12k} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Expected 95% CL exclusion limits on wino $W/Z$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{i33} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on wino $W/Z$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{i33} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Expected 95% CL exclusion limits on wino $W/h$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{12k} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on wino $W/h$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{12k} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Expected 95% CL exclusion limits on wino $W/h$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{i33} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on wino $W/h$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{i33} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Expected 95% CL exclusion limits on $\tilde{\ell}/\tilde{\nu}$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{12k} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on $\tilde{\ell}/\tilde{\nu}$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{12k} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Expected 95% CL exclusion limits on $\tilde{\ell}/\tilde{\nu}$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{i33} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on $\tilde{\ell}/\tilde{\nu}$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{i33} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Expected 95% CL exclusion limits on gluino NSLP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{12k} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on gluino NSLP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{12k} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Expected 95% CL exclusion limits on gluino NSLP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{i33} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on gluino NSLP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{i33} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Expected 95% CL exclusion limits on the higgsino GGM models. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on the higgsino GGM models. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0$ masses in the context of the wino $W/Z$ RPV scenario with $\lambda_{12k} \neq 0$ couplings. The signal region used to obtain the limits is SR0A. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0$ masses in the context of the wino $W/Z$ RPV scenario with $\lambda_{i33} \neq 0$ couplings. The signal region used to obtain the limits is the combination of SR0A, SR1 and SR2. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0$ masses in the context of the wino $W/Z$ RPV scenario with $\lambda_{12k} \neq 0$ couplings. The signal region used to obtain the limits is SR0B. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0$ masses in the context of the wino $W/Z$ RPV scenario with $\lambda_{i33} \neq 0$ couplings. The signal region used to obtain the limits is the combination of SR0B, SR1 and SR2. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0$ masses in the context of the wino $W/h$ RPV scenario with $\lambda_{12k} \neq 0$ couplings. The signal region used to obtain the limits is SR0A. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0$ masses in the context of the wino $W/h$ RPV scenario with $\lambda_{i33} \neq 0$ couplings. The signal region used to obtain the limits is the combination of SR0A, SR1 and SR2. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0$ masses in the context of the wino $W/h$ RPV scenario with $\lambda_{12k} \neq 0$ couplings. The signal region used to obtain the limits is SR0B. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0$ masses in the context of the wino $W/h$ RPV scenario with $\lambda_{i33} \neq 0$ couplings. The signal region used to obtain the limits is the combination of SR0B, SR1 and SR2. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the slepton/sneutrino masses in the context of the slepton/sneutrino RPV scenario with $\lambda_{12k} \neq 0$ couplings. The signal region used to obtain the limits is SR0A. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the slepton/sneutrino masses in the context of the slepton/sneutrino RPV scenario with $\lambda_{i33} \neq 0$ couplings. The signal region used to obtain the limits is the combination of SR0A, SR1 and SR2. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the slepton/sneutrino masses in the context of the slepton/sneutrino RPV scenario with $\lambda_{12k} \neq 0$ couplings. The signal region used to obtain the limits is SR0B. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the slepton/sneutrino masses in the context of the slepton/sneutrino RPV scenario with $\lambda_{i33} \neq 0$ couplings. The signal region used to obtain the limits is the combination of SR0B, SR1 and SR2. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the gluino masses in the context of the gluino RPV scenario with $\lambda_{12k} \neq 0$ couplings. The signal region used to obtain the limits is SR0A. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the gluino masses in the context of the gluino RPV scenario with $\lambda_{i33} \neq 0$ couplings. The signal region used to obtain the limits is the combination of SR0A, SR1 and SR2. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the gluino masses in the context of the gluino RPV scenario with $\lambda_{12k} \neq 0$ couplings. The signal region used to obtain the limits is SR0B. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the gluino masses in the context of the gluino RPV scenario with $\lambda_{i33} \neq 0$ couplings. The signal region used to obtain the limits is the combination of SR0B, SR1 and SR2. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}\tilde{\chi}_1^{0}$ masses in the context of the higgsino GGM scenario in SR0C. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}\tilde{\chi}_1^{0}$ masses in the context of the higgsino GGM scenario in SR0D. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

The best expected exclusion power between signal regions at each signal point as is adopted in the limit calculation of the wino $W/Z$ model for $\lambda_{12k} \neq 0$ RPV couplings. For the $\lambda_{12k} \neq 0$ case, the results from SR0B are adopted everywhere in the final exclusion limit contours, as is found to be the most powerful signal region among SR0A and SR0B in the majority of the signal grid points of this model.

The best expected exclusion power between signal regions at each signal point as is adopted in the limit calculation of the wino $W/Z$ model for $\lambda_{i33} \neq 0$ RPV couplings. A combination of SR0A + SR1 + SR2 or SR0B + SR1 + SR2 is tested to detect the best expected limit; shown on the plot is the region (SR0A or SR0B) which is chosen to provide the best expected limit in combination with SR1 and SR2. The results from the combination SR0B + SR1 + SR2 are finally adopted everywhere in the final exclusion limit contour since they provide the best expected limit.

The best expected exclusion power between signal regions at each signal point as is adopted in the limit calculation of the wino $W/h$ model for $\lambda_{12k} \neq 0$ RPV couplings. For the $\lambda_{12k} \neq 0$ case, the results from SR0B are adopted everywhere in the final exclusion limit contours, as is found to be the most powerful signal region among SR0A and SR0B in the majority of the signal grid points of this model.

The best expected exclusion power between signal regions at each signal point as is adopted in the limit calculation of the wino $W/h$ model for $\lambda_{i33} \neq 0$ RPV couplings. A combination of SR0A + SR1 + SR2 or SR0B + SR1 + SR2 is tested to detect the best expected limit; shown on the plot is the region (SR0A or SR0B) which is chosen to provide the best expected limit in combination with SR1 and SR2. The results from the combination SR0B + SR1 + SR2 are finally adopted everywhere in the final exclusion limit contour since they provide the best expected limit.

The best expected exclusion power between signal regions at each signal point as is adopted in the limit calculation of the slepton/sneutrino model for $\lambda_{12k} \neq 0$ RPV couplings.

The best expected exclusion power between signal regions at each signal point as is adopted in the limit calculation of the slepton/sneutrino model for $\lambda_{i33} \neq 0$ RPV couplings. A combination of SR0A + SR1 + SR2 or SR0B + SR1 + SR2 is tested to detect the best expected limit; shown on the plot is the region (SR0A or SR0B) which is chosen to provide the best expected limit in combination with SR1 and SR2.

The best expected exclusion power between signal regions at each signal point as is adopted in the limit calculation of the gluino model for $\lambda_{12k} \neq 0$ RPV couplings.

The best expected exclusion power between signal regions at each signal point as is adopted in the limit calculation of the gluino model for $\lambda_{i33} \neq 0$ RPV couplings. A combination of SR0A + SR1 + SR2 or SR0B + SR1 + SR2 is tested to detect the best expected limit; shown on the plot is the region (SR0A or SR0B) which is chosen to provide the best expected limit in combination with SR1 and SR2.

The best expected exclusion power between the overlapping SR0C and SR0D at each signal point as is adopted in the limit combination of the GGM higgsino model.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0A. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ $W/Z$ $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0A. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ $W/h$ $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0A. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the slepton/sneutrino $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0A. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the gluino $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0A. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0B. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ $W/Z$ $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0B. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ $W/h$ $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0B. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the slepton/sneutrino $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0B. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the gluino $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0B. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the GGM Higgsino models with BR($\tilde{\chi}_1^0 \rightarrow Z \tilde{G}$)=100% and BR($\tilde{\chi}_1^0 \rightarrow Z \tilde{G}$)=50% fulfilling the selection criteria of SR0C. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the GGM Higgsino models with BR($\tilde{\chi}_1^0 \rightarrow Z \tilde{G}$)=100% and BR($\tilde{\chi}_1^0 \rightarrow Z \tilde{G}$)=50% fulfilling the selection criteria of SR0D. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR1. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ $W/Z$ $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR1. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ $W/h$ $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR1. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the slepton/sneutrino $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR1. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the gluino $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR1. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR2. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ $W/Z$ $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR2. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ $W/h$ $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR2. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the slepton/sneutrino $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR2. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the gluino $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR2. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Cutflow event yields in regions SR0A and SR0B for RPV models with the $\lambda_{12k} \neq 0$ coupling. All yields correspond to weighted events, so that effects from lepton reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalized to the integrated luminosity of the data sample, $\int L dt = 36.1$ fb$^{-1}$. The "Initial" number indicates the weighted number of events before any selection cut is applied. The last entries show the normalized number of events surviving the selection requirements of SR0A and SR0B.

Cutflow event yields in regions SR0C and SR0D for GGM Higgsino models with BR($\tilde{\chi}_1^0 \rightarrow Z \tilde{G}$)=100% or BR($\tilde{\chi}_1^0 \rightarrow Z \tilde{G}$)=50%, and $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0\tilde{\chi}_1^0$ mass of 400 GeV. All yields correspond to weighted events, so that effects from lepton reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalized to the integrated luminosity of the data sample, $\int L dt = 36.1$ fb$^{-1}$. The "Initial" number indicates the weighted number of events before any selection cut is applied. The "Generator Filter" step is applied during the MC generation of the simulated events; the BR=100% sample has a generator filter of $\geq 4e/\mu$ leptons with $p_{\mathrm{T}}>4$ GeV and $|\eta|<2.8$, and the BR=50% sample has a generator filter of $\geq 4 e/\mu/\tau_{\mathrm{had-vis}}$ leptons with $p_{\mathrm{T}}(e,\mu)>4$ GeV, $p_{\mathrm{T}}(\tau_{\mathrm{had-vis}}^{\mathrm{visible}})>15$ GeV and $|\eta|<2.8$. The $ZZ$ selection cutflow step refers to the mass window cut for the leading and subleading $Z$ boson candidate between $81.2-101.2$ GeV and $61.2-101.2$ GeV, respectively. The last entries show the efficiency of events surviving the selection requirements defined in SR0C and SR0D.

Cutflow event yields in region SR1 for RPV models with the $\lambda_{i33} \neq 0$ coupling. All yields correspond to weighted events, so that effects from lepton reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalized to the integrated luminosity of the data sample, $\int L dt = 36.1$ fb$^{-1}$. The "Initial" number indicates the weighted number of events before any selection cut is applied. The last entries show the normalized number of events surviving the selection requirements of SR1.

Cutflow event yields in region SR2 for RPV models with the $\lambda_{i33} \neq 0$ coupling. All yields correspond to weighted events, so that effects from lepton reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalized to the integrated luminosity of the data sample, $\int L dt = 36.1$ fb$^{-1}$. The "Initial" number indicates the weighted number of events before any selection cut is applied. The last entries show the normalized number of events surviving the selection requirements of SR2.

Cross sections for the slepton/snuetrino model for different NLSP masses.

Distribution of the missing transverse momentum $E_{\mathrm{T}}^{\mathrm{miss}}$ for the sample satisfying all requirements of the $\mathrm{SR}^{\gamma\gamma}_{W-L}$ selection except the $E_{\mathrm{T}}^{\mathrm{miss}}$ requirement itself. Overlaid are the expected SM backgrounds, separated into the various contributing sources. Also shown are the signal expectations for the ($m_{\tilde{W}}$, $m_{\tilde{\chi}^{0}_{1}}$) = (1000,100) GeV and ($m_{\tilde{W}}$, $m_{\tilde{\chi}^{0}_{1}}$) = (1000,800) GeV models. The vertical dashed lines and right-pointing arrows show the region of the $E_{\mathrm{T}}^{\mathrm{miss}}$ observable selected for inclusion in $\mathrm{SR}^{\gamma\gamma}_{W-L}$ and $\mathrm{SR}^{\gamma\gamma}_{W-H}$. The lower panels show the ratio of observed data to the combined SM expectation. For these plots, the band represents the range of combined statistical and systematic uncertainty in the SM expectation. Events outside the range of the displayed region are included in the highest-value bin.

Distribution of the missing transverse momentum $E_{\mathrm{T}}^{\mathrm{miss}}$ for the sample satisfying all requirements of the $\mathrm{SR}^{\gamma\gamma}_{W-H}$ selection except the $E_{\mathrm{T}}^{\mathrm{miss}}$ requirement itself. Overlaid are the expected SM backgrounds, separated into the various contributing sources. Also shown are the signal expectations for the ($m_{\tilde{W}}$, $m_{\tilde{\chi}^{0}_{1}}$) = (1000,100) GeV and ($m_{\tilde{W}}$, $m_{\tilde{\chi}^{0}_{1}}$) = (1000,800) GeV models. The vertical dashed lines and right-pointing arrows show the region of the $E_{\mathrm{T}}^{\mathrm{miss}}$ observable selected for inclusion in $\mathrm{SR}^{\gamma\gamma}_{W-L}$ and $\mathrm{SR}^{\gamma\gamma}_{W-H}$. The lower panels show the ratio of observed data to the combined SM expectation. For these plots, the band represents the range of combined statistical and systematic uncertainty in the SM expectation. Events outside the range of the displayed region are included in the highest-value bin.

Comparisons between expected and observed content of the validation and signal regions for the photon+jets analysis. The uncertainties in the expected numbers of events are the combined statistical and systematic uncertainties. The lower panel shows the pull (difference between observed and expected event counts normalized by the uncertainty) for each region.

Distribution of the missing transverse momentum $E_{\mathrm{T}}^{\mathrm{miss}}$ for the sample satisfying all requirements of the $\mathrm{SR}^{\gamma j}_{H}$ selection except the $E_{\mathrm{T}}^{\mathrm{miss}}$ requirement itself. Overlaid are the expected SM backgrounds, separated into the various contributing sources. Also shown are the signal expectations for points in the $m_{\tilde{g}}-m_{\tilde{\chi}^{0}_{1}}$ parameter space of the GGM model relevant to the photon+jets analysis (mass values in GeV). The value of the gluino mass arises from the choice $M_{3}$ = 1900 GeV. The $\tilde{\chi}^{0}_{1}$ mass values of 1868, 1920, 442 and 652 GeV arise from the choices $\mu$ = 1810, 1868, 400 and 600 GeV, respectively, combined with the constraint that the branching fraction of $\tilde{\chi}^{0}_{1}$ $\to \gamma \tilde{G}$ be 50%. The vertical dashed lines and right-pointing arrows show the region of the $E_{\mathrm{T}}^{\mathrm{miss}}$ observable selected for inclusion in $\mathrm{SR}^{\gamma j}_{H}$ and $\mathrm{SR}^{\gamma j}_{L}$ for $\mathrm{SR}^{\gamma j}_{L200}$, the $E_{\mathrm{T}}^{\mathrm{miss}}$ requirement is 200 GeV rather than 300 GeV. The lower panels show the ratio of observed data to the combined SM expectation. For these plots, the band represents the range of statistical uncertainty in the SM expectation. Events outside the range of the displayed region are included in the highest-value bin.

Distribution of the missing transverse momentum $E_{\mathrm{T}}^{\mathrm{miss}}$ for the sample satisfying all requirements of the $\mathrm{SR}^{\gamma j}_{L}$ or $\mathrm{SR}^{\gamma j}_{L200}$ selection except the $E_{\mathrm{T}}^{\mathrm{miss}}$ requirement itself. Overlaid are the expected SM backgrounds, separated into the various contributing sources. Also shown are the signal expectations for points in the $m_{\tilde{g}}-m_{\tilde{\chi}^{0}_{1}}$ parameter space of the GGM model relevant to the photon+jets analysis (mass values in GeV). The value of the gluino mass arises from the choice $M_{3}$ = 1900 GeV. The $\tilde{\chi}^{0}_{1}$ mass values of 1868, 1920, 442 and 652 GeV arise from the choices $\mu$ = 1810, 1868, 400 and 600 GeV, respectively, combined with the constraint that the branching fraction of $\tilde{\chi}^{0}_{1} \to \gamma \tilde{G}$ be 50%. The vertical dashed lines and right-pointing arrows show the region of the $E_{\mathrm{T}}^{\mathrm{miss}}$ observable selected for inclusion in $\mathrm{SR}^{\gamma j}_{H}$ and $\mathrm{SR}^{\gamma j}_{L}$ for $\mathrm{SR}^{\gamma j}_{L200}$, the $E_{\mathrm{T}}^{\mathrm{miss}}$ requirement is 200 GeV rather than 300 GeV. The lower panels show the ratio of observed data to the combined SM expectation. For these plots, the band represents the range of statistical uncertainty in the SM expectation. Events outside the range of the displayed region are included in the highest-value bin.

Expected exclusion limits in the gluino-bino mass plane, using the $\mathrm{SR}^{\gamma\gamma}_{S-H}$ analysis for $m_{\tilde{\chi}^{0}_{1}} > 1600$ GeV and the $\mathrm{SR}^{\gamma\gamma}_{S-L}$ analysis for $m_{\tilde{\chi}^{0}_{1}} < 1600$ GeV.

Observed exclusion limits in the gluino--bino mass plane, using the $\mathrm{SR}^{\gamma\gamma}_{S-H}$ analysis for $m_{\tilde{\chi}^{0}_{1}} > 1600$ GeV and the $\mathrm{SR}^{\gamma\gamma}_{S-L}$ analysis for $m_{\tilde{\chi}^{0}_{1}} < 1600$ GeV.

Expected exclusion limit in the squark-bino mass plane, using the $\mathrm{SR}^{\gamma\gamma}_{S-H}$ analysis for $m_{\tilde{\chi}^{0}_{1}} > 900$ GeV and the $\mathrm{SR}^{\gamma\gamma}_{S-L}$ analysis for $m_{\tilde{\chi}^{0}_{1}} < 900$ GeV.

Observed exclusion limit in the squark--bino mass plane, using the $\mathrm{SR}^{\gamma\gamma}_{S-H}$ analysis for $m_{\tilde{\chi}^{0}_{1}} > 900$ GeV and the $\mathrm{SR}^{\gamma\gamma}_{S-L}$ analysis for $m_{\tilde{\chi}^{0}_{1}} < 900$ GeV.

Expected exclusion limit in the wino-bino mass plane, using the $\mathrm{SR}^{\gamma\gamma}_{W-H}$ analysis for $m_{\tilde{\chi}^{0}_{1}} > 400$ GeV and the $\mathrm{SR}^{\gamma\gamma}_{W-L}$ analysis for $m_{\tilde{\chi}^{0}_{1}}$ < 400$ GeV.

Observed exclusion limit in the wino-bino mass plane, using the $\mathrm{SR}^{\gamma\gamma}_{W-H}$ analysis for $m_{\tilde{\chi}^{0}_{1}} > 400$ GeV and the $\mathrm{SR}^{\gamma\gamma}_{W-L}$ analysis for $m_{\tilde{\chi}^{0}_{1}}$ < 400$ GeV.

Expected exclusion limits for the $\mu > 0$ higgsino-bino GGM model explored by the photon+jets analysis.

Observed exclusion limits for the $\mu > 0$ higgsino-bino GGM model explored by the photon+jets analysis.

Distribution of the transverse momentum $p_{\mathrm{T}} (\ell\gamma\gamma)$ of events in the $\ell\gamma\gamma$ control region (except without a cut on $p_{\mathrm{T}} (\ell\gamma\gamma)$). Also shown is the expected contribution from various SM sources, including $W(\to\ell\nu) + \gamma\gamma$ production itself. The displayed uncertainties are a combination of those from all SM sources except $W(\to\ell\nu) + \gamma\gamma$ production, and include statistical and systematic uncertainties.

Distribution of $E_{\mathrm{T}}^{\mathrm{miss}}$ for diphoton events in a validation region defined by a requirement of $H_{\mathrm{T}} > 1750$ GeV. Also shown is the expected contribution from various SM sources, as well as their combined statistical and systematic uncertainties.

Distribution of $H_{\mathrm{T}}$ for diphoton events in a validation region defined by requirement of $E_{\mathrm{T}}^{\mathrm{miss}} > 100$ GeV. Also shown is the expected contribution from various SM sources, as well as their combined statistical and systematic uncertainties.

Distribution of $m_{\mathrm{eff}}$ for events satisfying all requirements $\mathrm{SR}^{\gamma j}_{H}$ save the $m_{\mathrm{eff}}$ requirement itself. Also shown is the expected contribution from various SM sources, and their combined statistical uncertainties.

Distribution of $m_{\mathrm{eff}}$ for events satisfying all requirements $\mathrm{SR}^{\gamma j}_{L}$ save the $m_{\mathrm{eff}}$ requirement itself. Also shown is the expected contribution from various SM sources, and their combined statistical uncertainties.

Derived exclusion limits for the gluino-bino GGM model explored by the diphoton analysis. For each point in the gluino-bino parameter space, the SR ($\mathrm{SR}^{\gamma\gamma}_{S-L}$ or $\mathrm{SR}^{\gamma\gamma}_{S-H}$) that provides the best expected sensitivity is used to estimate the exclusion likelihood. The model dependent upper limits on cross-section (fb) are shown by grey numbers for each signal point.

Derived exclusion limits for the squark-bino GGM model explored by the diphoton analysis. For each point in the squark-bino parameter space, the SR ($\mathrm{SR}^{\gamma\gamma}_{S-L}$ or $\mathrm{SR}^{\gamma\gamma}_{S-H}$) that provides the best expected sensitivity is used to estimate the exclusion likelihood. The model dependent upper limits on cross-section (fb) are shown by grey numbers for each signal point.

Derived exclusion limits for the wino-bino GGM model explored by the diphoton analysis. For each point in the wino-bino parameter space, the SR ($\mathrm{SR}^{\gamma\gamma}_{W-L}$ or $\mathrm{SR}^{\gamma\gamma}_{W-H}$) that provides the best expected sensitivity is used to estimate the exclusion likelihood. The model dependent upper limits on cross-section (fb) are shown by grey numbers for each signal point.

Derived exclusion limits for the $\mu > 0$ higgsino-bino GGM model explored by the photon+jets analysis. For each point in the higgsino-bino parameter space, the SR ($\mathrm{SR}^{\gamma j}_{L}$ or $\mathrm{SR}^{\gamma j}_{H}$) that provides the best expected sensitivity is used to estimate the exclusion likelihood. The model dependent upper limits on cross-section (fb) are shown by grey numbers for each signal point.

Derived exclusion limits for the gluino-bino GGM model explored by the diphoton analysis. For each point in the gluino-bino parameter space, the SR ($\mathrm{SR}^{\gamma\gamma}_{S-L}$ or $\mathrm{SR}^{\gamma\gamma}_{S-H}$) that provides the best expected sensitivity is used to estimate the exclusion likelihood. The labels indicate the best-expected signal region for each point, where SL and SH mean $\mathrm{SR}^{\gamma\gamma}_{S-L}$ and $\mathrm{SR}^{\gamma\gamma}_{S-H}$, respectively.

Derived exclusion limits for the squark-bino GGM model explored by the diphoton analysis. For each point in the squark-bino parameter space, the SR ($\mathrm{SR}^{\gamma\gamma}_{S-L}$ or $\mathrm{SR}^{\gamma\gamma}_{S-H}$) that provides the best expected sensitivity is used to estimate the exclusion likelihood. The labels indicate the best-expected signal region for each point, where SL and SH mean $\mathrm{SR}^{\gamma\gamma}_{S-L}$ and $\mathrm{SR}^{\gamma\gamma}_{S-H}$, respectively.

Derived exclusion limits for the wino--bino GGM model explored by the diphoton analysis. For each point in the wino-bino parameter space, the SR ($\mathrm{SR}^{\gamma\gamma}_{W-L}$ or $\mathrm{SR}^{\gamma\gamma}_{W-H}$) that provides the best expected sensitivity is used to estimate the exclusion likelihood. The labels indicate the best-expected signal region for each point, where WL and WH mean $\mathrm{SR}^{\gamma\gamma}_{W-L}$ and $\mathrm{SR}^{\gamma\gamma}_{W-H}$, respectively.

Derived exclusion limits for the $\mu > 0$ higgsino-bino GGM model explored by the photon+jets analysis. For each point in the higgsino-bino parameter space, the SR ($\mathrm{SR}^{\gamma j}_{L}$ or $\mathrm{SR}^{\gamma j}_{H}$) that provides the best expected sensitivity is used to estimate the exclusion likelihood. The labels indicate the best-expected signal region for each point, where L and H mean $\mathrm{SR}^{\gamma j}_{L}$ and $\mathrm{SR}^{\gamma j}_{H}$, respectively.

Acceptance and efficiency for $\mathrm{SR}^{\gamma\gamma}_{S-L}$ for the signal models of the gluino-bino GGM grid.

Acceptance and efficiency for $\mathrm{SR}^{\gamma\gamma}_{S-H}$ for the signal models of the gluino-bino GGM grid.

Acceptance and efficiency for $\mathrm{SR}^{\gamma\gamma}_{S-L}$ for the signal models of the squark-bino GGM grid.

Acceptance and efficiency for $\mathrm{SR}^{\gamma\gamma}_{S-H}$ for the signal models of the squark-bino GGM grid.

Acceptance and efficiency for $\mathrm{SR}^{\gamma\gamma}_{W-L}$ for the signal models of the wino-bino GGM grid.

Acceptance and efficiency for $\mathrm{SR}^{\gamma\gamma}_{W-H}$ for the signal models of the wino-bino GGM grid.

Acceptance and efficiency for $\mathrm{SR}^{\gamma j}_{L}$ for the signal models of the photon+jets GGM grid.

Acceptance and efficiency for $\mathrm{SR}^{\gamma j}_{H}$ for the signal models of the photon+jets GGM grid.

Cutflow for the $\mathrm{SR}^{\gamma\gamma}_{S-L}$ selection for one relevant signal point in the gluino-bino model, where the gluinos have mass of 1900 GeV and the $\tilde{\chi}^{0}_{1}$ has a mass of 300 GeV (10000 generated events). The numbers are normalized to a luminosity of 36.1 fb$^{-1}$.

Cutflow for the $\mathrm{SR}^{\gamma\gamma}_{S-H}$ selection for one relevant signal point in the gluino-bino model, where the gluinos have mass of 1900 GeV and the $\tilde{\chi}^{0}_{1}$ has a mass of 1700 GeV (10000 generated events). The numbers are normalized to a luminosity of 36.1 fb$^{-1}$.

Cutflow for the $\mathrm{SR}^{\gamma\gamma}_{S-L}$ selection for one relevant signal point in the squark-bino model, where the squarks have mass of 1700 GeV and the $\tilde{\chi}^{0}_{1}$ has a mass of 200 GeV (10000 generated events). The numbers are normalized to a luminosity of 36.1 fb$^{-1}$.

Cutflow for the $\mathrm{SR}^{\gamma\gamma}_{S-H}$ selection for one relevant signal point in the squark-bino model, where the squarks have mass of 1700 GeV and the $\tilde{\chi}^{0}_{1}$ has a mass of 1600 GeV (10000 generated events). The numbers are normalized to a luminosity of 36.1 fb$^{-1}$.

Cutflow for the $\mathrm{SR}^{\gamma\gamma}_{W-L}$ selection for one relevant signal point in the wino-bino model, where the winos have mass of 1000 GeV and the $\tilde{\chi}^{0}_{1}$ has a mass of 200 GeV (10000 generated events). The numbers are normalized to a luminosity of 36.1 fb$^{-1}$.

Cutflow for the $\mathrm{SR}^{\gamma\gamma}_{W-H}$ selection for one relevant signal point in the wino-bino model, where the winos have mass of 1000 GeV and the $\tilde{\chi}^{0}_{1}$ has a mass of 800 GeV (10000 generated events). The numbers are normalized to a luminosity of 36.1 fb$^{-1}$.

Cutflow for the $\mathrm{SR}^{\gamma j}_{L}$ selection, for two relevant signal points in the higgsino-bino model, where the gluinos have mass of 1974 GeV and the $\tilde{\chi}^{0}_{1}$ has a mass of 442 GeV (10000 generated events), and 652 GeV (10000 generated events). The numbers are normalized to a luminosity of 36.1 fb$^{-1}$.

Cutflow for the $\mathrm{SR}^{\gamma j}_{H}$ selection, for two relevant signal points in the higgsino-bino model, where the gluinos have mass of 1974 GeV and the $\tilde{\chi}^{0}_{1}$ has a mass of 1868 GeV (10000 generated events), and 1920 GeV (10000 generated events). The numbers are normalized to a luminosity of 36.1 fb$^{-1}$.

Figure 8. Gluino pair production cross section.

Table 1. Post-fit yields for the background-only fit, observed data, and expected yields for mgluino = 1600 GeV in each search bin.

Cutflow table for a pair produced coloron of 1500 GeV decaying into two quarks.

The observed number of data, background and top squark signal events in each of the signal regions of the inclusive selection

The observed number of data, background and coloron signal events in each of the signal regions of the inclusive selection

The observed number of data, background and top squark signal events in each of the signal regions of the b-tagged selection

Signal acceptance and efficiency (in %) as a function of M(STOP), before mass windows

Signal acceptance (in %) and efficiency as a function of M(STOP), after mass windows

Signal acceptance and efficiency (in %) as a function of M(RHO), before mass windows

Signal acceptance and efficiency (in %) as a function of M(RHO), after mass windows

Signal acceptance (in %) and efficiency as a function of M(STOP), before mass windows

Signal acceptance (in %) and efficiency as a function of M(STOP), after mass windows

Cross section excluded at 95% CL as a function of the top squark mass, for a pair produced top squark with decays into a pair of light-quarks.

Cross section excluded at 95% CL as a function of the cooron mass, for a pair produced coloron with decays into a pair of light-quarks.

Cross section excluded at 95% CL as a function of the top squark mass, for a pair produced top squark with decays into a b- and an s-quark.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 700 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 800 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 800 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 900 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 900 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1000 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1000 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1050 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1050 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1100 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1100 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1150 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1150 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1200 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1200 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1250 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1250 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1300 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1300 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1350 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1350 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1400 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1400 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1450 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1450 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1500 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1500 GeV stop. All limits are computed at 95% CL.

Signal acceptance (in %) in the (BRe,BRtau) plane for a 800 GeV stop, for the SR800 signal region.

Signal acceptance (in %) in the (BRe,BRtau) plane for a 800 GeV stop, for the SR1100 signal region.

Signal efficiency (in %) in the (BRe,BRtau) plane for a 800 GeV stop, for the SR800 signal region.

Signal efficiency (in %) in the (BRe,BRtau) plane for a 800 GeV stop, for the SR1100 signal region.

Signal acceptance (in %) in the (BRe,BRtau) plane for a 1200 GeV stop, for the SR800 signal region.

Signal acceptance (in %) in the (BRe,BRtau) plane for a 1200 GeV stop, for the SR1100 signal region.

Signal efficiency (in %) in the (BRe,BRtau) plane for a 1200 GeV stop, for the SR800 signal region.

Signal efficiency (in %) in the (BRe,BRtau) plane for a 1200 GeV stop, for the SR1100 signal region.

Signal acceptance (in %) in the (BRe,BRtau) plane for a 1500 GeV stop, for the SR800 signal region.

Signal acceptance (in %) in the (BRe,BRtau) plane for a 1500 GeV stop, for the SR1100 signal region.

Signal efficiency (in %) in the (BRe,BRtau) plane for a 1500 GeV stop, for the SR800 signal region.

Signal efficiency (in %) in the (BRe,BRtau) plane for a 1500 GeV stop, for the SR1100 signal region.

$m_{bl}^{0}$ distribution in SR800. All selection criteria are applied, except the selection on the variable that is displayed in the plot. The SM backgrounds are normalized to the values determined in the background-only fit. The last bin includes overflows.

$m_{bl}^\mathrm{asym}$ distribution in SR800. All selection criteria are applied, except the selection on the variable that is displayed in the plot. The SM backgrounds are normalized to the values determined in the background-only fit. The last bin includes overflows.

$H_\mathrm{T}$ distribution in SR800. All selection criteria are applied, except the selection on the variable that is displayed in the plot. The SM backgrounds are normalized to the values determined in the background-only fit. The last bin includes overflows.

$m_{ll}$ distribution in SR800. All selection criteria are applied, except the selection on the variable that is displayed in the plot. The SM backgrounds are normalized to the values determined in the background-only fit. The last bin includes overflows.

$m_{bl}^{1}$(rej) distribution in SR800. All selection criteria are applied, except the selection on the variable that is displayed in the plot. The SM backgrounds are normalized to the values determined in the background-only fit. The last bin includes overflows.

Full list of event selections and MC generator-weighted yields and efficiencies in the inclusive SR800 and SR1100 signal regions for several signal samples of varying stop mass with decay into b-electron, b-muon or b-tau at 1/3 branching ratio.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 600 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 700 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 800 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 900 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1000 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1050 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1100 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1150 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1200 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1250 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1300 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1350 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1400 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1450 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1500 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1550 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1600 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1350 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1400 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1450 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1500 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1550 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1600 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 600 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 700 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 800 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 900 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1000 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1050 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1100 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1150 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1200 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1250 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1300 GeV stop. All limits are computed at 95% CL.

A summary of benchmark simplified models, the most sensitive $n_{\mathrm{jet}}$ categories, and representative values for the corresponding experimental acceptance times efficiency ($\mathcal{A}\varepsilon$), the dominant systematic uncertainties, the theoretical production cross section ($\sigma_{\mathrm{theory}}$), and the expected and observed upper limits on the production cross section, expressed in terms of the signal strength parameter ($\mu$).

Summary of the mass limits obtained for the fourteen classes of simplified models. The limits indicate the strongest observed and expected (in parentheses) mass exclusions in $\tilde{g}$, $\tilde{q}$, $\tilde{b}$, $\tilde{t}$, and $\tilde{\chi}_1^0$. The quoted values have uncertainties of $\pm 25$ and $\pm 10$ GeV for models involving the pair production of, respectively, gluinos and squarks.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1qqqq model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1qqqq model. This line is the observed exclusion.

Covariance matrix for the SM background estimates obtained using the simplified binning scheme, determined from a simultaneous fit to data in the control regions only (CR-only fit). The uncertainties in the background estimates are correlated in such a way that the covariance is typically positive. Small positive values, as well as the few negative values, are not shown.

Correlation matrix for the SM background estimates obtained using the simplified binning scheme, determined from a simultaneous fit to data in the control regions only (CR-only fit).

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{q} } }, m_{\tilde{\chi}^0_1 })$ plane for the T2qq model. This line is the expected exclusion.

Observed data counts and "post-fit" background expectations based on the result of a combined fit to the signal region and multiple control regions under the SM-only hypothesis for the "monojet" event category. The rows labelled SM "pre-fit' show the background expectations when excluding the signal region from the fit. The uncertainties include statistical as well as systematic contributions.

Observed data counts and "post-fit" background expectations based on the result of a combined fit to the signal region and multiple control regions under the SM-only hypothesis for the "asymmetric" event categories. The rows labelled SM "pre-fit" show the background expectations when excluding the signal region from the fit. The uncertainties include statistical as well as systematic contributions.

Observed data counts and "post-fit" background expectations based on the result of a combined fit to the signal region and multiple control regions under the SM-only hypothesis for the "symmetric" event categories. The rows labelled SM "pre-fit" show the background expectations when excluding the signal region from the fit. The uncertainties include statistical as well as systematic contributions.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{q} } }, m_{\tilde{\chi}^0_1 })$ plane for the T2qq model. This line is the observed exclusion.

Summary of a simplified binning scheme comprising $H_{\mathrm{T}}^{\mathrm{miss}}$ shapes for eight representative event topologies, defined in terms of requirements on $n_{\mathrm{jet}}$ both symmetric and asymmetric categories) and $n_{\mathrm{b}}$. For each topology, event yields are integrated over the full $H_{\mathrm{T}}$ range, $H_{\mathrm{T}} > 200$ GeV, and are categorised according to eight bins in $H_{\mathrm{T}}^{\mathrm{miss}}$ -- $130 < H_{\mathrm{T}}^{\mathrm{miss}} < 200$ GeV, six bins 100 GeV wide in the region $200 < H_{\mathrm{T}}^{\mathrm{miss}} < 800$ GeV, and an open final bin, $H_{\mathrm{T}}^{\mathrm{miss}} > 800$ GeV. This categorisation scheme leads to 64 bins that are exclusive, contiguous, and provide complete coverage of the signal region. The corresponding SM background estimates for these 64 bins are obtained using the nominal likelihood model. The $H_{\mathrm{T}}^{\mathrm{miss}}$ template for each topology is determined from simulation and the normalisation for each template is given by an aggregate of estimates from several ($n_{\mathrm{jet}}$, $n_{\mathrm{b}}$, $H_{\mathrm{T}}$) bins, defined according to the nominal binning scheme.

Observed data yields and "CR-only fit" SM background expectations using the simplified binning scheme, as a function of topology and $H_{\mathrm{T}}^{\mathrm{miss}}$. The SM backgrounds are determined from a simultaneous fit to data in the control regions. The data yields in the signal region are masked and excluded from the fit. The uncertainties include statistical as well as systematic contributions.

Expected ($\mu_{\mathrm{exp}}$) and observed ($\mu_{\mathrm{obs}}$) upper limits on the production cross section, expressed in terms of the signal strength parameter, obtained using both the nominal and simplified binning schema. Comparable values are obtained from the two schema for these benchmark SUSY models because the signal events typically populate bins at high values of $H_{\mathrm{T}}$ and $H_{\mathrm{T}}^{\mathrm{miss}}$ and so are at or near the limit of the search sensitivity. Larger differences in $\mu$ can be observed for models that populate the core of these distributions.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1bbbb model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1bbbb model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1tttt model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1tttt model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1ttbb model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1qqqq model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1ttbb model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1qqqq model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1qqqq model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5tttt_DM175 model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1qqqq model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{q} } }, m_{\tilde{\chi}^0_1 })$ plane for the T2qq model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{q} } }, m_{\tilde{\chi}^0_1 })$ plane for the T2qq model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5tttt_DM175 model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{q} } }, m_{\tilde{\chi}^0_1 })$ plane for the T2qq model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{q} } }, m_{\tilde{\chi}^0_1 })$ plane for the T2qq model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1bbbb model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5ttcc model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1bbbb model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1bbbb model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1bbbb model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5ttcc model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1tttt model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1tttt model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1tttt model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ b } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bb model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1tttt model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1ttbb model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1ttbb model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ b } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bb model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1ttbb model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T1ttbb model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5tttt_DM175 model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tb model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5tttt_DM175 model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5tttt_DM175 model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5tttt_DM175 model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tb model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5ttcc model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5ttcc model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5ttcc model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{g} } }, m_{\tilde{\chi}^0_1 })$ plane for the T5ttcc model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ b } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bb model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ b } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bb model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ b } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bb model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ b } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bb model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tb model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2cc model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tb model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tb model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tb model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2cc model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt_degen model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2cc model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2cc model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt_degen model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2cc model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2cc model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt_degen model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2_mixed model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt_degen model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt_degen model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2tt_degen model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2_mixed model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2_mixed model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2_mixed model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2_mixed model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bW model. This line is the expected exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2_mixed model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bW model. This line is the $+1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bW model. This line is the $-1\sigma $ expected exclusion due to experimental uncertainties.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bW model. This line is the observed exclusion.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bW model. This line is the $+1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Upper limit on the cross section in the $(m_{ \tilde{ \mathrm{ t } } }, m_{\tilde{\chi}^0_1 })$ plane for the T2bW model. This line is the $-1\sigma $ observed exclusion due to theoretical uncertainties on the signal cross section.

Kinematic distributions for "extended SR selections" of the 2-tau channel, for mTsum in the Compressed SR selection without the mTsum>1400 GeV requirement. The last bin includes overflow events. Cited uncertainties are statistical uncertainties only. Signal predictions are overlaid for several benchmark models, normalised to their predicted cross sections. For the simplified model, "MM" refers to a medium mass splitting, with m(gluino)=1145 GeV and m(neutralino)=265 GeV; "HM" denotes a high mass splitting scenario, with m(gluino)=1305 GeV and m(neutralino)=105 GeV. The GMSB benchmark model corresponds to Lambda = 90 TeV and tanbeta = 40.

Kinematic distributions for "extended SR selections" of the 2-tau channel, for mTtau1+mTtau2 in the High-Mass SR selection without the mTtau1+mTtau2>350GeV requirement. The last bin includes overflow events. Cited uncertainties are statistical uncertainties only. Signal predictions are overlaid for several benchmark models, normalised to their predicted cross sections. For the simplified model, "MM" refers to a medium mass splitting, with m(gluino)=1145 GeV and m(neutralino)=265 GeV; "HM" denotes a high mass splitting scenario, with m(gluino)=1305 GeV and m(neutralino)=105 GeV. The GMSB benchmark model corresponds to Lambda = 90 TeV and tanbeta = 40.

Kinematic distributions for "extended SR selections" of the 2-tau channel, for HT in the GMSB SR selection without the HT > 1700 GeV requirement. The last bin includes overflow events. Cited uncertainties are statistical uncertainties only. Signal predictions are overlaid for several benchmark models, normalised to their predicted cross sections. For the simplified model, "MM" refers to a medium mass splitting, with m(gluino)=1145 GeV and m(neutralino)=265 GeV; "HM" denotes a high mass splitting scenario, with m(gluino)=1305 GeV and m(neutralino)=105 GeV. The GMSB benchmark model corresponds to Lambda = 90 TeV and tanbeta = 40.

Expected exclusion contour at the 95% confidence level for the simplified model of gluino pair production, based on the combined results from the 1tau and 2tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data.

Observed exclusion contour at the 95% confidence level for the simplified model of gluino pair production, based on the combined results from the 1tau and 2tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data.

Expected exclusion contour at the 95% confidence level for the simplified model of gluino pair production, based on results from the 2tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data.

Expected exclusion contour at the 95% confidence level for the simplified model of gluino pair production, based on results from the 1tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data.

Observed exclusion contours at the 95% confidence level for the gauge-mediated supersymmetry-breaking model, based on results from the 2 tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data. Additional model parameters are M(mess) = 250 TeV, N5 = 3, mu>0 and Cgrav =1.

Expected exclusion contours at the 95% confidence level for the gauge-mediated supersymmetry-breaking model, based on results from the 2 tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data. Additional model parameters are M(mess) = 250 TeV, N5 = 3, mu>0 and Cgrav =1.

Observed upper cross section limits in pb for the simplified model of gluino pair production for the combination of all SRs.

Best expected signal region for the simplified model of gluino pair production. The respective SR has been used in the combination of the results.

Acceptance for the gluino production simplified model grid in the Compressed 1tau signal region.

Efficiency for the gluino production simplified model grid in the Compressed 1tau signal region.

Acceptance times Efficiency for the gluino production simplified model grid in the Compressed 1tau signal region.

Acceptance for the gluino production simplified model grid in the medium mass 1tau signal region.

Efficiency for the gluino production simplified model grid in the medium mass 1tau signal region.

Acceptance times Efficiency for the gluino production simplified model grid in the medium mass 1tau signal region.

Acceptance for the gluino production simplified model grid in the high mass 1tau signal region.

Efficiency for the gluino production simplified model grid in the high mass 1tau signal region.

Acceptance times Efficiency for the gluino production simplified model grid in the high mass 1tau signal region.

Acceptance for the gluino production simplified model grid in the compressed 2tau signal region.

Efficiency for the gluino production simplified model grid in the compressed 2tau signal region.

Acceptance times Efficiency for the gluino production simplified model grid in the compressed 2tau signal region.

Acceptance for the gluino production simplified model grid in the high mass 2tau signal region.

Efficiency for the gluino production simplified model grid in the high mass 2tau signal region.

Acceptance times Efficiency for the gluino production simplified model grid in the high mass 2tau signal region.

Acceptance for the GMSB model grid in the 2tau signal region.

Efficiency for the GMSB model grid in the 2tau signal region.

Acceptance times Efficiency for the GMSB model grid in the 2tau signal region.