Bayesian upper exclusion limits on the ratio $g_{W'}/g_{W}$ for an SSM-like W$\prime$ boson are shown. The unity coupling ratio (blue dotted curve) corresponds to the SSM common benchmark. The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian lower exclusion limits on the NUGIM G(221) mixing angle $\cot(\theta_{E})$ are shown as a function of the W$\prime$ boson mass. The theoretically excluded region is shaded in grey. The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian upper exclusion limits at 95% CL on the product of the production cross section and branching fraction of a QBH in an associated $ au$ lepton and neutrino final state. Minimum threshold masses $m_{th}$ of up to 6.6 TeV are excluded at 95% CL. The observed limit (solid line) is obtained from the intersection with the LO QBH cross section (blue dotted curve). The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian upper limits at 95% CL on the cross section of the process $pp\rightarrow\tau\nu$ mediated via LQ exchange in the t-channel. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The predicted LQ cross section at LO in the three coupling benchmark scenarios is depicted in different colors for $g_{U}=1$. The uncertainty bandy correspond to the sum in quadrature of PDF and scale variations. The first benchmark scenario considers only couplings to left-handed SM fermions (i.e. $\beta_{\text{R}}^{ij}=0$) and is referred to as "best fit LH". The second benchmark, referred to as "best fit LH+RH", considers $|\beta_{\text{R}}^{\text{b}\tau}|=1$ and all other $\beta_{\text{R}}^{ij}=0$. In the third "democratic" benchmark, equal couplings only to LH fermions are assumed, i.e. $\beta_{\text{L}}^{ij}=1$ and $\beta_{\text{R}}^{ij}=0$ for all $i$ and $j$.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the LH scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the LH+RH scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the democratic scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Bayesian upper exclusion limits at 95% CL on each of the Wilson coefficients described by the EFT model based on 2016-2018 data. The three different coupling types represent a left-handed vector coupling ($\epsilon^{cb}_{L}$), tensor-like coupling ($\epsilon^{cb}_{T}$), and scalar-tensor-like coupling ($\epsilon^{cb}_{S_{L}}$). The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Summary of exclusion limits (expected and observed) calculated at 95% CL for full Run-2 CMS data.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "year_binnumber", following the binning from Figure 4.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "year_binnumber", following the binning used in Figure 4.

Predicted signal yields for the 2017 data-taking period, corresponding to 41.3 fb$^{-1}$, after the application of each search requirement (cumulative) for various signal hypothesis. The requirements listed are presented as total efficiencies w.r.t. the previous selection step.

Upper limits on the production cross section times branching fraction of the b* RH hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the b* VL hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the B+b hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the B+t hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Dijet signal efficiencies as a function of the signal mass and lifetime for events satisfying all event and vertex requirements, with corrections based on systematic differences in the vertex reconstruction efficiency between data and simulation.

The distribution of $d_{\mathrm{BV}}$ for $\geq$5-track one-vertex events in data and three simulated multijet signal samples each with a mass of 1600 GeV. The production cross section for each signal model is assumed to be the lower limit excluded by CMS-EXO-17-018, corresponding to values of 0.8, 0.25, and 0.15 fb for the samples with $c\tau =$ 0.3, 1.0, and 10 mm, respectively. The last bin includes the overflow events. This bin includes one event in data with a vertex with large $d_{\mathrm{BV}}$ that appears to arise from tracks originating from separate pp interaction vertices, consistent with background.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data with two 3-track vertices. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data which have exactly one 4-track vertex and one 3-track vertex. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data with two 4-track vertices. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized using $\geq$5-track one-vertex event information. The two vertical red dashed lines separate the regions used in the fit.

Predicted yields for the background-only normalized template, predicted yields for three simulated multijet signals each with a mass of 1600 GeV, and the observed yield in each $d_{\mathrm{VV}}$ bin. The production cross section for each signal model is assumed to be the lower limit excluded by CMS-EXO-17-018, corresponding to values of 0.8, 0.25, and 0.15 fb for samples with $c\tau =$ 0.3, 1.0, and 10 mm, respectively. The uncertainties in the signal yields and the systematic uncertainties in the background prediction reflect the systematic uncertainties given in the text.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the dijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the top squark with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the dijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the top squark with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 300um in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 300um in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 1 mm in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 1 mm in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 10 mm in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 10 mm in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 800 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 800 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 1600 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 1600 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 2400 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals. For $m$ = 2400 GeV, the expected neutralino cross section is $\approx 8\times 10^{-5}$ fb and is not shown.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 2400 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Data-to-simulation efficiency correction factors, shown for multijet and dijet signal topologies in several ranges of $c\tau$. Note that these correction factors account for the two long-lived particles in the simulated events, and are therefore the total correction factors used to scale event yields rather than the correction factors one would apply to individual vertices.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 3-track one-vertex events in data, and is normalized to the number of 3-track two-vertex events in data.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 4-track and 3-track one-vertex events in data, and is normalized to the number of two-vertex events in data which have exactly one 4-track vertex and one 3-track vertex.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 4-track one-vertex events in data, and is normalized to the number of 4-track two-vertex events in data.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from $\geq$5-track one-vertex events in data, and is normalized using one-vertex event information. No $\geq$5-track two-vertex data events pass the selection.

Expected and observed 95 percent CL upper limits on BR(H to etau) for each individual category and combined from collinear mass fit analysis

Summary of the observed and expected upper limits at the 95percent CL and the best fit branching fractions in percent for the H->mutau and H->etau processes, for the main analysis (BDT fit) and the cross check (Mcol fit) method.

95 percent CL observed upper limit on the Yukawa couplings for the main analysis (BDT fit) and the cross check (Mcol fit) method

Normalized inclusive 2-jet cross section differential in $\Delta\phi_{1,2}$ for $500 < p_{T}^{max} < 600$ GeV

Normalized inclusive 2-jet cross section differential in $\Delta\phi_{1,2}$ for $600 < p_{T}^{max} < 700$ GeV

Normalized inclusive 2-jet cross section differential in $\Delta\phi_{1,2}$ for $700 < p_{T}^{max} < 800$ GeV

Normalized inclusive 2-jet cross section differential in $\Delta\phi_{1,2}$ for $800 < p_{T}^{max} < 1000$ GeV

Normalized inclusive 2-jet cross section differential in $\Delta\phi_{1,2}$ for $1000 < p_{T}^{max} < 1200$ GeV

Normalized inclusive 2-jet cross section differential in $\Delta\phi_{1,2}$ for $p_{T}^{max} > 1200$ GeV

Normalized inclusive 3-jet cross section differential in $\Delta\phi_{1,2}$ for $200 < p_{T}^{max} < 300$ GeV

Normalized inclusive 3-jet cross section differential in $\Delta\phi_{1,2}$ for $300 < p_{T}^{max} < 400$ GeV

Normalized inclusive 3-jet cross section differential in $\Delta\phi_{1,2}$ for $400 < p_{T}^{max} < 500$ GeV

Normalized inclusive 3-jet cross section differential in $\Delta\phi_{1,2}$ for $500 < p_{T}^{max} < 600$ GeV

Normalized inclusive 3-jet cross section differential in $\Delta\phi_{1,2}$ for $600 < p_{T}^{max} < 700$ GeV

Normalized inclusive 3-jet cross section differential in $\Delta\phi_{1,2}$ for $700 < p_{T}^{max} < 800$ GeV

Normalized inclusive 3-jet cross section differential in $\Delta\phi_{1,2}$ for $800 < p_{T}^{max} < 1000$ GeV

Normalized inclusive 3-jet cross section differential in $\Delta\phi_{1,2}$ for $p_{T}^{max} > 1000$ GeV

Normalized inclusive 4-jet cross section differential in $\Delta\phi_{1,2}$ for $200 < p_{T}^{max} < 300$ GeV

Normalized inclusive 4-jet cross section differential in $\Delta\phi_{1,2}$ for $300 < p_{T}^{max} < 400$ GeV

Normalized inclusive 4-jet cross section differential in $\Delta\phi_{1,2}$ for $400 < p_{T}^{max} < 500$ GeV

Normalized inclusive 4-jet cross section differential in $\Delta\phi_{1,2}$ for $500 < p_{T}^{max} < 600$ GeV

Normalized inclusive 4-jet cross section differential in $\Delta\phi_{1,2}$ for $600 < p_{T}^{max} < 700$ GeV

Normalized inclusive 4-jet cross section differential in $\Delta\phi_{1,2}$ for $700 < p_{T}^{max} < 800$ GeV

Normalized inclusive 4-jet cross section differential in $\Delta\phi_{1,2}$ for $800 < p_{T}^{max} < 1000$ GeV

Normalized inclusive 4-jet cross section differential in $\Delta\phi_{1,2}$ for $p_{T}^{max} > 1000$ GeV

Normalized inclusive 3-jet cross section differential in $\Delta\phi_{2j}^{min}$ for $200 < p_{T}^{max} < 300$ GeV

Normalized inclusive 3-jet cross section differential in $\Delta\phi_{2j}^{min}$ for $300 < p_{T}^{max} < 400$ GeV

Normalized inclusive 3-jet cross section differential in $\Delta\phi_{2j}^{min}$ for $400 < p_{T}^{max} < 500$ GeV

Normalized inclusive 3-jet cross section differential in $\Delta\phi_{2j}^{min}$ for $500 < p_{T}^{max} < 600$ GeV

Normalized inclusive 3-jet cross section differential in $\Delta\phi_{2j}^{min}$ for $600 < p_{T}^{max} < 700$ GeV

Normalized inclusive 3-jet cross section differential in $\Delta\phi_{2j}^{min}$ for $700 < p_{T}^{max} < 800$ GeV

Normalized inclusive 3-jet cross section differential in $\Delta\phi_{2j}^{min}$ for $800 < p_{T}^{max} < 1000$ GeV

Normalized inclusive 3-jet cross section differential in $\Delta\phi_{2j}^{min}$ for $p_{T}^{max} > 1000$ GeV

Normalized inclusive 4-jet cross section differential in $\Delta\phi_{2j}^{min}$ for $200 < p_{T}^{max} < 300$ GeV

Normalized inclusive 4-jet cross section differential in $\Delta\phi_{2j}^{min}$ for $300 < p_{T}^{max} < 400$ GeV

Normalized inclusive 4-jet cross section differential in $\Delta\phi_{2j}^{min}$ for $400 < p_{T}^{max} < 500$ GeV

Normalized inclusive 4-jet cross section differential in $\Delta\phi_{2j}^{min}$ for $500 < p_{T}^{max} < 600$ GeV

Normalized inclusive 4-jet cross section differential in $\Delta\phi_{2j}^{min}$ for $600 < p_{T}^{max} < 700$ GeV

Normalized inclusive 4-jet cross section differential in $\Delta\phi_{2j}^{min}$ for $700 < p_{T}^{max} < 800$ GeV

Normalized inclusive 4-jet cross section differential in $\Delta\phi_{2j}^{min}$ for $800 < p_{T}^{max} < 1000$ GeV

Normalized inclusive 4-jet cross section differential in $\Delta\phi_{2j}^{min}$ for $p_{T}^{max} > 1000$ GeV

Unfolded distributions of $\Sigma p_{T}$ density in Z events, as a function of $p_{T}^{\mu\mu}$ in the towards ($\Delta\phi< 60^{\circ}$) region. Error bars represent the statistical and systematic uncertainties added in quadrature.

Unfolded distributions of $\Sigma p_{T}$ density in Z events, as a function of $p_{T}^{\mu\mu}$ in the transverse ($60^{\circ} <\Delta\phi< 120^{\circ}$) region. Error bars represent the statistical and systematic uncertainties added in quadrature.

Unfolded distributions of $\Sigma p_{T}$ density in Z events, as a function of $p_{T}^{\mu\mu}$ in the away ($\Delta\phi> 120^{\circ}$) region. Error bars represent the statistical and systematic uncertainties added in quadrature.

Figure 3 (upper). Distribution of $M_{T2}(ll)$ in a control region enriched in $t\bar{t}$ events defined by requiring at least 2 jets, at least one of them b-tagged, and missing transverse momentum below 80 GeV. For the plot, simulated yields are normalized to data using the yields in the $M_{T2}(ll)$ < 100 GeV region.

Figure 4 (upper). Expected and observed yields in the five $t\bar{t}Z$ control regions, which are defined by different requirements on the number of reconstructed jets and b jets, before the fit.

Figure 4 (lower). Expected and observed yields in the five $t\bar{t}Z$ control regions, which are defined by different requirements on the number of reconstructed jets and b jets, after the fit.

Figure 5 (left). Distribution of $M_{T2}(ll)$ of SF events falling within in the Z boson mass window with at least two jets and no b jet, $p_{T}^{miss} > 80$ GeV and $M_{T2}(ll) > 100$ GeV.

Figure 5 (center). Distribution of $M_{T2}(blbl)$ of SF events falling within in the Z boson mass window with at least two jets and no b jet, $p_{T}^{miss} > 80$ GeV and $M_{T2}(ll) > 100$ GeV.

Figure 5 (right). Distribution of $p_{T}^{miss}$ of SF events falling within in the Z boson mass window with at least two jets and no b jet, $p_{T}^{miss} > 80$ GeV and $M_{T2}(ll) > 100$ GeV.

Figure 6. Event yields in the 13 Drell-Yan and multiboson control regions for events with SF leptons falling within the Z boson mass window and no b jets, after renormalizing with the scale factors obtained from the fit procedure described in the text.

Figure 7 (left). Distributions of $M_{T2}(ll)$ for observed events in the dimuon channel compared to the predicted SM background.

Figure 7 (center). Distributions of $M_{T2}(ll)$ for observed events in the dielectron channel compared to the predicted SM background.

Figure 7 (right). Distributions of $M_{T2}(ll)$ for observed events in the electron-muon channel compared to the predicted SM background.

Figure 8 (left). Distributions of $M_{T2}(ll)$ for all lepton flavors for the defined selection.

Figure 8 (center). Distributions of $M_{T2}(blbl)$ for all lepton flavors for the defined selection, after requiring $M_{T2}(ll) > 100$ GeV.

Figure 8 (right). Distributions of $p_{T}^{miss}$ for all lepton flavors for the defined selection, after requiring $M_{T2}(ll) > 100$ GeV.

Figure 9 (upper). Predicted backgrounds and observed yields in the dielectron and dimuon search regions. The yields are post fit and correspond to the numbers quoted in Tab. 4 of the publication.

Figure 9 (lower). Predicted backgrounds and observed yields in the electron-muon search regions. The yields are post fit and correspond to the numbers quoted in Tab. 4 of the publication.

Figure 10. Predicted backgrounds and observed yields in all search regions combined. The yields are post fit and correspond to the numbers quoted in Tab. 4 of the publication.

Figure 11 (upper). Observed limits for the T2tt model in the top squark - LSP mass plane. The numbers indicate the 95% CL upper limit on the cross section times the square of the branching fraction at each point of the plane.

Figure 11 (upper). Observed exclusion region at 95% CL assuming 100% branching fraction.

Figure 11 (upper). Expected exclusion region at 95% CL assuming 100% branching fraction.

Figure 11 (lower). Observed limits for the T2bW model in the top squark - LSP mass plane. The numbers indicate the 95% CL upper limit on the cross section times the square of the branching fraction at each point of the plane.

Figure 11 (lower). Observed exclusion region at 95% CL assuming 100% branching fraction.

Figure 11 (lower). Expected exclusion region at 95% CL assuming 100% branching fraction.

Figure 12 (upper left). Observed limits for the T8bbllnunu model in the top squark - LSP mass plane. The numbers indicate the 95% CL upper limit on the cross section times the square of the branching fraction at each point of the plane. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.05 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 12 (upper left). Observed exclusion region at 95% CL assuming 100% branching fraction. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.05 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 12 (upper left). Expected exclusion region at 95% CL assuming 100% branching fraction. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.05 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 12 (upper left). Observed limits for the T8bbllnunu model in the top squark - LSP mass plane. The numbers indicate the 95% CL upper limit on the cross section times the square of the branching fraction at each point of the plane. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.5 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 12 (upper left). Observed exclusion region at 95% CL assuming 100% branching fraction. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.5 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 12 (upper left). Expected exclusion region at 95% CL assuming 100% branching fraction. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.5 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 12 (lower). Observed limits for the T8bbllnunu model in the top squark - LSP mass plane. The numbers indicate the 95% CL upper limit on the cross section times the square of the branching fraction at each point of the plane. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.95 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 12 (lower). Observed exclusion region at 95% CL assuming 100% branching fraction. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.95 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 12 (lower). Expected exclusion region at 95% CL assuming 100% branching fraction. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.95 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 13. Observed limits for the T2tt model in the top squark - LSP mass plane combining the dilepton final state with single-lepton and the all-hadronic final states. The numbers indicate the 95% CL upper limit on the cross section times the square of the branching fraction at each point of the plane.

Figure 13. Observed exclusion region at 95% CL assuming 100% branching fraction.

Figure 13. Expected exclusion region at 95% CL assuming 100% branching fraction.

Figure 14. Observed limits for the T2bW model in the top squark - LSP mass plane combining the dilepton final state with single-lepton and the all-hadronic final states. The numbers indicate the 95% CL upper limit on the cross section times the square of the branching fraction at each point of the plane.

Figure 14. Observed exclusion region at 95% CL assuming 100% branching fraction.

Figure 14. Expected exclusion region at 95% CL assuming 100% branching fraction.

Additional Figure 3. Observed limits for the T2tt model in the top squark - LSP mass plane using aggregate signal regions. The numbers indicate the 95% CL upper limit on the cross section times the square of the branching fraction at each point of the plane.

Additional Figure 3. Observed exclusion region at 95% CL assuming 100% branching fraction.

Additional Figure 3. Expected exclusion region at 95% CL assuming 100% branching fraction.

Table 5. Ratios of $\mu=\frac{\sigma$}{\sigma_{theory}}$ of the 95% CL expected and observed limits to simplified model expectations for different DM particle masses and mediator masses.

Table 6. Expected and observed event yields, summed over all lepton flavors, in aggregate signal regions.

Table 7 (upper). Covariance matrix for the background prediction in aggregate signal regions.

Table 7 (lower). Correlation matrix for the background prediction in aggregate signal regions.

Additional Figure 1. Correlation matrix of the 13 signal regions, split into same-flavor (SF) and different-flavor (DF) regions.

Additional Figure 2. Covariance matrix of the 13 signal regions, split into same-flavor (SF) and different-flavor (DF) regions.

Double-differential cross section times the dimuon branching fraction of the Y(2S) meson for different ranges of pT in bins of |y| and for the full |y| < 1.2 range, for the unpolarized decay hypothesis. The global uncertainty in the integrated luminosity of 2.3% is not included in the systematic uncertainties.

Double-differential cross section times the dimuon branching fraction of the Y(3S) meson for different ranges of pT in bins of |y| and for the full |y| < 1.2 range, for the unpolarized decay hypothesis. The global uncertainty in the integrated luminosity of 2.3% is not included in the systematic uncertainties.

Multiplicative scaling factors to obtain the J/psi differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, +0.10, -1) from the unpolarized cross section measurements given in Table 1.

Multiplicative scaling factors to obtain the psi(2S) differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, +0.03, -1) from the unpolarized cross section measurements given in Table 2.

Multiplicative scaling factors to obtain the Y(1S) differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, k, -1) from the unpolarized cross section measurements given in Table 3. The parameter k corresponds to a linear interpolation of the CMS measured value as a function of pT for pT < 50 GeV. For pT > 50 GeV, where no measurements exist, k is taken as the average of all the measured values for pT < 50 GeV.

Multiplicative scaling factors to obtain the Y(2S) differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, k, -1) from the unpolarized cross section measurements given in Table 3. The parameter k corresponds to a linear interpolation of the CMS measured value as a function of pT for pT < 50 GeV. For pT > 50 GeV, where no measurements exist, k is taken as the average of all the measured values for pT < 50 GeV.

Multiplicative scaling factors to obtain the Y(3S) differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, k, -1) from the unpolarized cross section measurements given in Table 3. The parameter k corresponds to a linear interpolation of the CMS measured value as a function of pT for pT < 50 GeV. For pT > 50 GeV, where no measurements exist, k is taken as the average of all the measured values for pT < 50 GeV.

Ratio of differential cross section times branching fractions of the prompt psi(2S) to J/psi mesons assuming unpolarized productions as a function of pT for |y| < 1.2.

Ratios of differential cross section times branching fractions of the prompt Y(2S) to Y(1S) and Y(3S) to Y(1S) mesons assuming unpolarized productions as a function of pT for |y| < 1.2.