Values of the predicted signal yields in each of the 15 signal regions (for $ x=0.75 $).

Pre-fit background-only covariance among different search bins.

Expected $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.25 $).

Top squark - neutralino mass points lying on the expected exclusion contour (for $ x=0.25 $).

Top squark - neutralino mass points lying on the expected-$1\sigma_\text{expected}$ exclusion contour (for $ x=0.25 $).

Top squark - neutralino mass points lying on the expected+$1\sigma_\text{expected}$ exclusion contour (for $ x=0.25 $).

Observed $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.25 $).

Top squark - neutralino mass points lying on the observed exclusion contour (for $ x=0.25 $).

Top squark - neutralino mass points lying on the observed-$1\sigma_\text{observed}$ exclusion contour (for $ x=0.25 $).

Top squark - neutralino mass points lying on the observed+$1\sigma_\text{observed}$ exclusion contour (for $ x=0.25 $).

Expected $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.5 $).

Top squark - neutralino mass points lying on the expected exclusion contour (for $ x=0.5 $).

Top squark - neutralino mass points lying on the expected-$1\sigma_\text{expected}$ exclusion contour (for $ x=0.5 $).

Top squark - neutralino mass points lying on the expected+$1\sigma_\text{expected}$ exclusion contour (for $ x=0.5 $).

Observed $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.5 $).

Top squark - neutralino mass points lying on the observed exclusion contour (for $ x=0.5 $).

Top squark - neutralino mass points lying on the observed-$1\sigma_\text{observed}$ exclusion contour (for $ x=0.5 $).

Top squark - neutralino mass points lying on the observed+$1\sigma_\text{observed}$ exclusion contour (for $ x=0.5 $).

Expected $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.75 $).

Top squark - neutralino mass points lying on the expected exclusion contour (for $ x=0.75 $).

Top squark - neutralino mass points lying on the expected-$1\sigma_\text{expected}$ exclusion contour (for $ x=0.75 $).

Top squark - neutralino mass points lying on the expected+$1\sigma_\text{expected}$ exclusion contour (for $ x=0.75 $).

Observed $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.75 $).

Top squark - neutralino mass points lying on the observed exclusion contour (for $ x=0.75 $).

Top squark - neutralino mass points lying on the observed-$1\sigma_\text{observed}$ exclusion contour (for $ x=0.75 $).

Top squark - neutralino mass points lying on the observed+$1\sigma_\text{observed}$ exclusion contour (for $ x=0.75 $).

The intercept parameter as a function of kT using the HCS method in pp collisions at 13 TeV.

1/Rinv^2 as a function of mT using the HCS method in pp collisions at 13 TeV.

The depth of the anticorrelation as function of multiplicity for kT integrated values.

Correlation function from the double ratio technique, integrated in the range 0 < kT < 1 GeV, for 0=<Ntrkoffline<=4. The error bars represent statistical uncertainties.

Correlation function from the double ratio technique, integrated in the range 0 < kT < 1 GeV, for 10=<Ntrkoffline<=12. The error bars represent statistical uncertainties.

Correlation function from the double ratio technique, integrated in the range 0 < kT < 1 GeV, for 31=<Ntrkoffline<=33. The error bars represent statistical uncertainties.

Correlation function from the double ratio technique, integrated in the range 0 < kT < 1 GeV, for 80=<Ntrkoffline<=84. The error bars represent statistical uncertainties.

Correlation function from the double ratio technique, integrated in the range 0 < kT < 1 GeV, for 105=<Ntrkoffline<=109. The error bars represent statistical uncertainties.

Correlation function from the double ratio technique, integrated in the range 0 < kT < 1 GeV, for 130=<Ntrkoffline<=250. The error bars represent statistical uncertainties.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of Z pt on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of Z pt on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of |y(Z)| on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of |y(Z)| on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of $\phi^{\scriptscriptstyle *}_\eta$ on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of $\phi^{\scriptscriptstyle *}_\eta$ on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of Z pt on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of Z pt on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of |y(Z)| on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of |y(Z)| on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of $\phi^{\scriptscriptstyle *}_\eta$ on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of $\phi^{\scriptscriptstyle *}_\eta$ on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

The measured absolute cross sections in bins of |y(Z)|, using dressed level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of Z pt, using dressed level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of $\phi^{\scriptscriptstyle *}_\eta$, using dressed level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of Z pt different |y(Z)| bins, using dressed level leptons. The cross sections are normalized by the bin width.

The measured normalized cross sections in bins of Z pt, using dressed level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of $\phi^{\scriptscriptstyle *}_\eta$, using dressed level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of |y(Z)|, using dressed level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of Z pt in |y(Z)| bins, using dressed level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of |y(Z)| in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of |y(Z)| in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of |y(Z)| in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of |y(Z)| in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of |y(Z)| in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of |y(Z)| in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of |y(Z)| in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of |y(Z)| in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

The measured absolute cross sections in bins of |y(Z)|, using born level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of Z pt, using born level leptons. The cross sections are normalized by the bin width. The first bin (0-1 GeV) is not shown as large differences were observed in aMC@NLO and POWHEG predictions at Born level.

The measured absolute cross sections in bins of $\phi^{\scriptscriptstyle *}_\eta$, using born level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of Z pt different |y(Z)| bins, using born level leptons. The cross sections are normalized by the bin width. The first bin (0-1 GeV) is not shown as large differences were observed in aMC@NLO and POWHEG predictions at Born level.

The measured normalized cross sections in bins of Z pt, using born level leptons. The first bin (0-1 GeV) is not shown as large differences were observed in aMC@NLO and POWHEG predictions at Born level. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of $\phi^{\scriptscriptstyle *}_\eta$, using born level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of |y(Z)|, using born level leptons.

The measured normalized cross sections (left) in bins of Z pt in |y(Z)| bins, using born level leptons. The first bin (0-1 GeV) is not shown as large differences were observed in aMC@NLO and POWHEG predictions at Born level. The cross sections are divided by the bin width and then the distributions are normalized to 1.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of |y(Z)| in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of |y(Z)| in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of |y(Z)| in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of |y(Z)| in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dilepton final state.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and a $\tilde{\chi}_1^\pm$ that decays to $W^\pm\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction to $q_i\bar{q}_j W\pm\tilde{\chi}_1^0$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to bottom quarks ($\tilde{g}\to b\bar{b}\tilde{\chi}_1^0$). Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to top quarks ($\tilde{g}\to t\bar{t}\tilde{\chi}_1^0$). Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for light-flavor squark pair production, where the squarks decay to $q\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay and 1-fold degeneracy in the light-flavor squarks (corresponding to the inner set of curves in the limit plot). To get the theory cross section for other N-fold degeneracy assumptions (e.g. 8-fold for the outer curves in the limit plot), just multiply by N.

Exclusion limits at 95% CL for bottom squark pair production, where the squarks decay to $b\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $t\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $b\tilde{\chi}_1^\pm$ and the $\tilde{\chi}_1^0$ decay to $W^\pm\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay either to $b\tilde{\chi}_1^\pm\to bW^\pm\tilde{\chi}_1^0$ or to $t\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $c\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for the mono-$\phi$ model, in which a resonantly-produced colored scalar decays to a massive Dirac fermion and a quark. Signal cross sections are calculated at leading order in $\alpha_S$, assuming unity branching fraction for the given decay.

Cross section limits for $\mathrm{LQ}\to\mathrm{q}\nu$, where $q=u,\,d,\,s,\,\mathrm{or}\,c$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratio is assumed to be 100% to $\mathrm{q}\nu$.

Cross section limits for $\mathrm{LQ}\to\mathrm{b}\nu$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratio is assumed to be 100% to $\mathrm{b}\nu$.

Cross section limits for $\mathrm{LQ}\to\mathrm{t}\nu$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratios are assumed to be $\mathcal{B}(\mathrm{LQ}\to\mathrm{t}\nu)=1-\beta$, and $\mathcal{B}(\mathrm{LQ}\to\mathrm{b}\tau)=\beta$.

Predictions and observations for monojet signal regions

Predictions and observations for signal regions with $250 \leq H_\mathrm{T} < 450$ GeV

Predictions and observations for signal regions with $450 \leq H_\mathrm{T} < 575$ GeV and $N_\mathrm{j}<7$

Predictions and observations for signal regions with $450 \leq H_\mathrm{T} < 575$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $N_\mathrm{j}^\mathrm{hi}<4$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $4\leq N_\mathrm{j}^\mathrm{hi}<7$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $N_\mathrm{j}^\mathrm{hi}<4$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $4\leq N_\mathrm{j}^\mathrm{hi}<7$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $H_\mathrm{T} \geq 1500$ GeV and $N_\mathrm{j}<7$

Predictions and observations for signal regions with $H_\mathrm{T} \geq 1500$ GeV and $N_\mathrm{j}\geq7$

Covariance matrix for the 282 signal regions of the inclusive $M_\mathrm{T2}$ search

Correlation matrix for the 282 signal regions of the inclusive $M_\mathrm{T2}$ search

Bin number definitions for the $M_\mathrm{T2}$ covariance and correlation matrices

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

The maximum chargino mass excluded at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the gluino mass is fixed to 1900 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

The maximum chargino mass excluded at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the squark mass is fixed to 900 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

The maximum chargino mass excluded at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the squark mass is fixed to 1500 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, and the eight light squarks' masses are assumed to be degenerate. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the stop mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the gluino mass is fixed to 1600 GeV and the neutralino's mass is fixed to 1575 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the squark mass is fixed to 2000 GeV and the neutralino's mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, and the eight light squarks' masses are assumed to be degenerate.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the stop mass is fixed to 1100 GeV and the neutralino's mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Predictions and observations for 2016 disappearing track signal regions

Predictions and observations for 2017-2018 pixel track signal regions

Predictions and observations for 2017-2018 medium (M) and long (L) length track signal regions

Covariance matrix for the 68 signal regions of the disappearing tracks $M_\mathrm{T2}$ search

Correlation matrix for the 68 signal regions of the disappearing tracks $M_\mathrm{T2}$ search

Observed yields and pre-fit background predictions for Njets 8-9.

Observed yields and pre-fit background predictions for Njets >=10.

Observed yields and pre-fit background predictions for aggregated bins.

Pre-fit background covariance matrix $\sigma_{xy}$ for the 174 analysis bins.

Pre-fit background correlation matrix $\rho_{xy}$ for the 174 analysis bins.

Observed yields and post-fit background predictions for Njets 2-3.

Observed yields and post-fit background predictions for Njets 4-5.

Observed yields and post-fit background predictions for Njets 6-7.

Observed yields and post-fit background predictions for Njets 8-9.

Observed yields and post-fit background predictions for Njets >=10.

Simulated signal yields for representative gluino signal model light-LSP points in the aggregate search regions. The uncertainties shown are statistical.

Simulated signal yields for representative gluino signal model compressed points in the aggregate search regions. The uncertainties shown are statistical.

Simulated signal yields for representative squark signal model light-LSP points in the aggregate search regions. The uncertainties shown are statistical.

Simulated signal yields for representative squark signal model compressed points in the aggregate search regions. The uncertainties shown are statistical.

The 95% CL observed and expected upper limits on the production cross sections of the T1tttt signal model as a function of the gluino and LSP masses.

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

The 95% CL observed and expected upper limits on the production cross sections of the T1bbbb signal model as a function of the gluino and LSP masses.

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

The 95% CL observed and expected upper limits on the production cross sections of the T1qqqq signal model as a function of the gluino and LSP masses.

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

The 95% CL observed and expected upper limits on the production cross sections of the T5qqqqVV signal model as a function of the gluino and LSP masses.

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

The 95% CL observed and expected upper limits on the production cross sections of the T2tt signal model as a function of the top squark and LSP masses.

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

The 95% CL observed and expected upper limits on the production cross sections of the T2bb signal model as a function of the bottom squark and LSP masses.

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

The 95% CL observed and expected upper limits on the production cross sections of the T2qq signal model as a function of the LSP and light squark masses.

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (u d s c)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (u d s c)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (u d s c)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (u d s c)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (u d s c)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (u d s c)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (u d s c)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (single flavor)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (single flavor)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (single flavor)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (single flavor)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (single flavor)

The observed significance of the SMS models for T1tttt as function of the gluino and LSP masses.

The observed significance of the SMS models for T1bbbb as function of the gluino and LSP masses.

The observed significance of the SMS models for T1qqqq as function of the gluino and LSP masses.

The observed significance of the SMS models for T5qqqqVV as function of the gluino and LSP masses.

The observed significance of the SMS models for T2tt as function of the top squark and LSP masses.

The observed significance of the SMS models for T2bb as function of the bottom squark and LSP masses.

The observed significance of the SMS models for T2qq as function of the LSP and light squark masses.

The efficiency times acceptance of the SMS models for T1tttt as function of the gluino and LSP masses. The efficiency times acceptance is given for the total signal across the set of 174 search regions.

The efficiency times acceptance of the SMS models for T1bbbb as function of the gluino and LSP masses. The efficiency times acceptance is given for the total signal across the set of 174 search regions.

The efficiency times acceptance of the SMS models for T1qqqq as function of the gluino and LSP masses. The efficiency times acceptance is given for the total signal across the set of 174 search regions.

The efficiency times acceptance of the SMS models for T5qqqqVV as function of the gluino and LSP masses. The efficiency times acceptance is given for the total signal across the set of 174 search regions.

The efficiency times acceptance of the SMS models for T2tt as function of the top squark and LSP masses. The efficiency times acceptance is given for the total signal across the set of 174 search regions.

The efficiency times acceptance of the SMS models for T2bb as function of the bottom squark and LSP masses. The efficiency times acceptance is given for the total signal across the set of 174 search regions.

The efficiency times acceptance of the SMS models for T2qq as function of the LSP and light squark masses. The efficiency times acceptance is given for the total signal across the set of 174 search regions.

Absolute cumulative efficiencies in percent for each step of the event selection process, listed for four representative gluino pair production signal models with m(gluino) >> m(LSP). Only statistical uncertainties are shown.

Absolute cumulative efficiencies in percent for each step of the event selection process, listed for four representative gluino pair production signal models with m(gluino) ~ m(LSP). Only statistical uncertainties are shown.

Absolute cumulative efficiencies in percent for each step of the event selection process, listed for three representative squark pair production signal models with m(squark) >> m(LSP). Only statistical uncertainties are shown.

Absolute cumulative efficiencies in percent for each step of the event selection process, listed for three representative squark pair production signal models with m(squark) ~ m(LSP). Only statistical uncertainties are shown.

$\mathrm{M_{\mathrm{T}}}$ distribution for the $\mathrm{e}\mu$ control region.

$\mathrm{M_{\mathrm{T}}}$ distribution for the WZ control region.

$\mathrm{M_{\mathrm{T}}}$ distribution for the ZZ control region.

$\mathrm{M_{\mathrm{T}}}$ distribution for the signal region, events with $|\eta_{\gamma}|<1$.

$\mathrm{M_{\mathrm{T}}}$ distribution for the signal region, events with $|\eta_{\gamma}|>1$.

Muon reconstruction, identification, and isolation efficiency as a function of $|\eta|$ and $\mathrm{p_\mathrm{T}}$.

Electron reconstruction, identification, and isolation efficiency as a function of $|\eta|$ and $\mathrm{p_\mathrm{T}}$.

Photon efficiency as a function of $|\eta|$ and $\mathrm{p_\mathrm{T}}$.

Covariance matrix for all bins used in the analysis. There are 45 bins in total, 15 for every data-taking year. For every year, the first bin corresponds to events in the $\mathrm{e}\mu$ control region, the following five bins correspond to events with $|\eta^\gamma|< 1$ in the signal region, the next five bins correspond to events with $|\eta^\gamma|> 1$ in the signal region, the next two bins correspond to events in the WZ control region, and finally the last two bins correspond to events in the ZZ control region.