Observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance. The acceptance $\mathcal{A}$ is defined as $\mathcal{A} = N$(events with $m_{X}^{GEN} > 85\% m_{X}^{input}$) / $N$(events generated in the full phase space defined by the CMS default generator settings). Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Mass exclusion ranges of the benchmark signal scenarios are depicted with hatched areas inside the black contours.

Observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance. The acceptance $\mathcal{A}$ is defined as $\mathcal{A} = N$(events with $m_{X}^{GEN} > 85\% m_{X}^{input}$) / $N$(events generated in the full phase space defined by the CMS default generator settings). Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Mass exclusion ranges of the benchmark signal scenarios are depicted with hatched areas inside the black contours.

Efficiencies of the selection requirements on the benchmark signal processes: $Z_{R} \to ggg$ with nominal width ($\Gamma_{X}/m_{X}\sim 3\%$). A value of -1 means that the corresponding efficiency is not calculated for this year.

Efficiencies of the selection requirements on the benchmark signal processes: $Z_{R} \to ggg$ with narrow width ($\Gamma_{X}/m_{X}\sim 0.01\%$). A value of -1 means that the corresponding efficiency is not calculated for this year.

Efficiencies of the selection requirements on the benchmark signal processes: $G_{KK} \to {\varphi}(gg)g$, where $\varphi$ is the radion. A value of -1 means that the corresponding efficiency is not calculated for this year.

Efficiencies of the selection requirements on the benchmark signal processes: $q^{*} \to V(qq)q$, where $V$ is a beyond-the-SM vector boson. A value of -1 means that the corresponding efficiency is not calculated for this year.

Acceptance of the signal selection requirement $m_{X}^{\text{GEN}}/m_{X}^{\text{input}} > 85\%$ on the benchmark signal process $Z_{R} \to ggg$. The acceptance is defined as $\mathcal{A} = N$(events with $m_{X}^{\text{GEN}}/m_{X}^{\text{input}} > 85\%$)/$N$(events generated in the full phase space defined by the CMS default generator settings).

Acceptance of the signal selection requirement $m_{X}^{\text{GEN}}/m_{X}^{\text{input}} > 85\%$ on the benchmark signal process $G_{KK} \to {\varphi}(gg)g$.

Acceptance of the signal selection requirement $m_{X}^{\text{GEN}}/m_{X}^{\text{input}} > 85\%$ on the benchmark signal process $q^{*} \to V(qq)q$.

Observed local significance for a 3-body decay $ggg$ resonance, shown for resonances with nominal width (blue solid line) and narrow width (red dashed line).

Observed local significance for a cascade decay $ggg$ resonance.

Observed local significance for a cascade decay $qqq$ resonance.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.2$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.3$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.4$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.5$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.6$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.7$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.8$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance $m_{Y} / m_{X} = 0.2$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.3$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.4$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.5$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.6$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.7$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.8$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal

Cut flow table of the selection requirments for the $Z_{R}$ model scenarios. Values shown are absolute efficiencies, e.g., values shown in the column of 'Efficiency ($\Delta_{R}^{max} < 3.0$)' represent the cumulative efficiencies achieved through all selection requirements applied up to and including the current selection criterion.

Cut flow table of the selection requirments for the $G_{KK}$ model scenarios. Values shown are absolute efficiencies, e.g., values shown in the column of 'Efficiency ($\Delta_{R}^{max} < 3.0$)' represent the cumulative efficiencies achieved through all selection requirements applied up to and including the current selection criterion.

Cut flow table of the selection requirments for the $q^{*}$ model scenarios. Values shown are absolute efficiencies, e.g., values shown in the column of 'Efficiency ($\Delta_{R}^{max} < 3.0$)' represent the cumulative efficiencies achieved through all selection requirements applied up to and including the current selection criterion.

Predicted production cross section of the benchmark signal process $pp \to Z_{R} \to ggg$.

Predicted production cross section of the benchmark signal process $pp \to G_{KK} \to \varphi(gg)g$.

Predicted production cross section of the benchmark signal process $pp \to q^{*} \to V(qq)q$.

Correlation matrix for the differential branching fraction extraction between different $q^2$ bins in the simultaneous fit.

The product of acceptance and efficiency ($A\epsilon$) of the $\mathcal{B}(\mathrm{B}^+\to\mathrm{K}^+\mu^+\mu^-)$ channel, as a function of the muon pair $q^2$, as measured in simulated signal events, after all the corrections applied. In the figure, regions corresponding to resonances are displayed with red markers.

Same as Table 5, but for the bins containing the J$\psi$ and $\psi$(2S) resonances.

Integrated branching fraction $\mathcal{B}(\mathrm{B}^\pm \to \mathrm{K}^\pm\mu^+\mu^-)$ in the low-$q^2$ region.

Ratio of branching fractions $\mathcal{B}(\text{B}^\pm \to \text{K}^\pm\mu^+\mu^-)$ to $\mathcal{B}(\text{B}^\pm \to \text{K}^\pm\text{e}^+\text{e}^-)$, determined from the measured double ratio $R$(K) of these decays to the respective branching fractions of the decay chains $\text{B}^\pm\to\text{J/}\psi\text{K}^\pm$ with $\text{J/}\psi\to\mu^+\mu^-$ and $\text{e}^+\text{e}^-$.

Negative log likelihood function from the fit profiled as a function of the inverse ratio $R(\mathrm{K})^{-1}$.

Postfit score of the $\text{t}\bar{\text{t}}\text{Z}$ output node from the multiclass classifier in $\text{SR}_\text{3l,3j}$ for events with at least 2 b jets

Postfit score of the tWZ output node from the binary classifier in $\text{SR}_\text{3l,2j}$

$\text{p}_\text{T}$ of the leading jet in the $\text{SR}_\text{3l,2j}$

$\text{p}_\text{T}$ of the leading lepton in the $\text{SR}_\text{3l,2j}$

Max $\text{p}_\text{T}$ of the lepton-jet in the $\text{SR}_\text{3l,3j}$

Invariant mass of the system of leptons, jets, and $\vec{\text{p}}^\text{miss}_\text{T}$ in the $\text{SR}_\text{3l,3j}$

Two-dimensional likelihood scan of the signal strenghts of tWZ ($\mu_\text{tWZ}$) vs. $\text{t}\bar{\text{t}}\text{Z}$ ($\mu_{\text{t}\bar{\text{t}}\text{Z}}$). The blue lines show the 68%, 95%, and 99% confidence level contours. The black cross represents the best-fit value, while the red diamond the SM expectation

Two-dimensional likelihood scan of the signal strenghts of tWZ ($\mu_\text{tWZ}$) vs. $\text{t}\bar{\text{t}}\text{Z}$ ($\mu_{\text{t}\bar{\text{t}}\text{Z}}$). The blue lines show the 68%, 95%, and 99% confidence level contours. The black cross represents the best-fit value, while the red diamond the SM expectation

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet no-btag category with the BDT requirements selecting the most signal-like events. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet no-btag category with the BDT requirements selecting the most signal-like events. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet no-btag category with the BDT requirements selecting the most signal-like events. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to b quarks and $\tau$ leptons, using $\lambda = 1.5$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to light-flavor quarks and $\tau$ leptons, using $\lambda = 1.5$.

Expected upper limit at 95% CL on the coupling strength $\lambda$ of a scalar LQ to b quarks and $\tau$ leptons.

Expected upper limit at 95% CL on the coupling strength $\lambda$ of a scalar LQ to light-flavor quarks and $\tau$ leptons.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet btag category with the intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet no-btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet no-btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet no-btag category with low BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet no-btag category with low BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet no-btag category with low BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet no-btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet no-btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to b quarks and $\tau$ leptons, using $\lambda = 1.0$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to b quarks and $\tau$ leptons, using $\lambda = 2.0$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to b quarks and $\tau$ leptons, using $\lambda = 3.0$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to light-flavor quarks and $\tau$ leptons, using $\lambda = 0.5$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to light-flavor quarks and $\tau$ leptons, using $\lambda = 1.0$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to light-flavor quarks and $\tau$ leptons, using $\lambda = 2.0$.

Observed and expected distributions of the jet multiplicity in the e+jet btag control region.

Observed and expected distributions of the jet multiplicity in the $\mu$+jet btag control region.

Example description

Example description

Example description

Comparison between the data and SM background predictions for the number of muons in the long track final state.

Comparison between the data and SM background predictions for the m_{DTk; stopping power} in the long track final state.

Comparison between the data and SM background predictions for the number of electrons in the long track final state.

Comparison between the data and SM background predictions for the number of b-tagged jets in the short track final state.

Comparison between the data and SM background predictions for the number of jets in the short track final state.

Comparison between the data and SM background predictions for the hard missing transverse momentum in the short track final state.

Comparison between the data and SM background predictions for the number of muons in the short track final state.

Comparison between the data and SM background predictions for the m_{DTk; stopping power} in the short track final state.

Comparison between the data and SM background predictions for the number of electrons in the short track final state.

Comparison between the data and post-fit SM background predictions for the 49 search regions in the global search region.

95% CL upper limit on the cross section in pb for the T6bt 10 cm model.

Limit in the mass plane of the T6bt 10 cm model

95% CL upper limit on the cross section in pb for the T6bt 200 cm model.

Limit in the mass plane of the T6bt 200 cm model

95% CL upper limit on the cross section in pb for the T6tb 10 cm model.

Limit in the mass plane of the T6tb 10 cm model

95% CL upper limit on the cross section in pb for the T6tb 200 cm model.

Limit in the mass plane of the T6tb 200 cm model

95% CL upper limit on the cross section in pb for the T5btbt 10 cm model.

Limit in the mass plane of the T5btbt 10 cm model

95% CL upper limit on the cross section in pb for the T5btbt 200 cm model.

Limit in the mass plane of the T5btbt 200 cm model

95% CL upper limit on the cross section in pb for the pure Wino model.

Limit in the mass plane of the pure Wino model

95% CL upper limit on the cross section in pb for the pure Higgsino model.

Limit in the mass plane of the pure Higgsino model

Cut flow table showing object-level efficiency at each stage of the track selection. The denominator of the efficiency is the number of chargino-matched tracks that have been reconstructed and satisfy the pT requirements, for which the chargino decays within the pixel detector volume; the numerator is the number of such tracks which also satisfy all requirements up to and including the associated criterion.

Cut flow table showing object-level efficiency at each stage of the track selection. The denominator of the efficiency is the number of chargino-matched tracks that have been reconstructed and satisfy the pT requirements, for which the chargino decays within the strip detector volume; the numerator is the number of such tracks which also satisfy all requirements up to and including the associated criterion.

Percent efficiencies for sum of all SR bins, each SR bin, SR bin groups:

Percent efficiencies for sum of all SR bins, each SR bin, SR bin groups:

Percent efficiencies for sum of all SR bins, each SR bin, SR bin groups:

Percent efficiencies for sum of all SR bins, each SR bin, SR bin groups:

Percent efficiencies for sum of all SR bins, each SR bin, SR bin groups:

Percent efficiencies for sum of all SR bins, each SR bin, SR bin groups:

Feynman diagrams of $\mathrm{Z'}\to\mu^{-}\mu^{+}$ with a $\mathrm{Z'}$ boson produced via $\mathrm{s}\overline{\mathrm{b}}\to\mathrm{Z'}$, with one $\mathrm{b}$ quark and one $\overline{\mathrm{s}}$ quark in the final state.

Distribution of $\min(m_{\mu\mathrm{b}})$ as obtained from simulation in events with $N_{\mathrm{b}} \geq 1$ passing all the other selection requirements. In this search, we require $\min(m_{\mu\mathrm{b}})>175~\mathrm{GeV}$. The stacked histogram displays the expected distribution from the simulation of the SM backgrounds, while the overlaid open histograms illustrate the size and shape of the $\mathrm{Z'}$ contribution from the LFU model described in Eq. (1), for several $\mathrm{Z'}$ mass hypotheses. For illustrative purposes, we choose couplings $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$. The contribution of background processes other than DY and $\mathrm{t}\overline{\mathrm{t}}$ is so small that it is only barely visible at the bottom of the stacked histogram. The hatched region indicates the statistical uncertainty arising from the limited size of the SM simulated samples. Histograms are normalized to unit area.

Distributions of $m_{\mu\mu}$ in the $N_{\mathrm{b}}=1$ event category. The stacked histogram displays the expected distribution from the SM background simulation. The overlaid open distributions illustrate the $\mathrm{Z'}$ contribution from the LFU model at $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$ for a variety of $\mathrm{Z'}$ mass hypotheses. The observed data are shown as black points with statistical error bars. The hatched region indicates the statistical uncertainty arising from the limited size of the SM simulated samples. The size of the bins increases as a function of $m_{\mu\mu}$. In extracting the results of the search, the background is estimated directly from data, so the SM background simulation is only illustrative in these distributions.

Distributions of $m_{\mu\mu}$ in the $N_{\mathrm{b}}\geq 2$ event category. The stacked histogram displays the expected distribution from the SM background simulation. The overlaid open distributions illustrate the $\mathrm{Z'}$ contribution from the LFU model at $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$ for a variety of $\mathrm{Z'}$ mass hypotheses. The observed data are shown as black points with statistical error bars. The hatched region indicates the statistical uncertainty arising from the limited size of the SM simulated samples. The size of the bins increases as a function of $m_{\mu\mu}$. In extracting the results of the search, the background is estimated directly from data, so the SM background simulation is only illustrative in these distributions.

Invariant mass $m_{\mu\mu}$ distributions in the $N_{\mathrm{b}} = 1$ category, shown together with the corresponding selected background functional forms used as input to the $discrete~profiling$ method when probing the $m_{\mathrm{Z'}}=500~\mathrm{GeV}$ hypothesis. The expected signal distribution for the LFU model described in Eq. (1), with couplings $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$, is overlaid. The displayed mass range corresponds to the fit window used for this $m_{\mathrm{Z'}}$ hypothesis, which is $\pm 10\, \sigma_{\mathrm{mass}}$ around the probed $m_{\mathrm{Z'}}$ value. While the likelihood fits are performed on unbinned data, here we present the data in binned histograms with binning chosen to reflect the size of $\sigma_{\mathrm{mass}}$.

Invariant mass $m_{\mu\mu}$ distributions in the $N_{\mathrm{b}} \geq 2$ category, shown together with the corresponding selected background functional forms used as input to the $discrete~profiling$ method when probing the $m_{\mathrm{Z'}}=500~\mathrm{GeV}$ hypothesis. The expected signal distribution for the LFU model described in Eq. (1), with couplings $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$, is overlaid. The displayed mass range corresponds to the fit window used for this $m_{\mathrm{Z'}}$ hypothesis, which is $\pm 10\, \sigma_{\mathrm{mass}}$ around the probed $m_{\mathrm{Z'}}$ value. While the likelihood fits are performed on unbinned data, here we present the data in binned histograms with binning chosen to reflect the size of $\sigma_{\mathrm{mass}}$.

Exclusion limits at $95\%$ CL on the number of selected BSM events with $N_{\mathrm{b}} \geq 1$ as functions of $m_{\mathrm{Z'}}$ for $f_{2\mathrm{b}}=0$. The quantity $f_{2\mathrm{b}}$ is the fraction of BSM events passing the analysis selection that have at least two $\mathrm{b}$ quark jets. The solid black (dashed red) curve represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing $68~(95)\%$ of the distribution of limits expected under the background-only hypothesis.

Exclusion limits at $95\%$ CL on the number of selected BSM events with $N_{\mathrm{b}} \geq 1$ as functions of $m_{\mathrm{Z'}}$ for $f_{2\mathrm{b}}=0.25$. The quantity $f_{2\mathrm{b}}$ is the fraction of BSM events passing the analysis selection that have at least two $\mathrm{b}$ quark jets. The solid black (dashed red) curve represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing $68~(95)\%$ of the distribution of limits expected under the background-only hypothesis.

Exclusion limits at $95\%$ CL on the number of selected BSM events with $N_{\mathrm{b}} \geq 1$ as functions of $m_{\mathrm{Z'}}$ for $f_{2\mathrm{b}}=0.75$. The quantity $f_{2\mathrm{b}}$ is the fraction of BSM events passing the analysis selection that have at least two $\mathrm{b}$ quark jets. The solid black (dashed red) curve represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing $68~(95)\%$ of the distribution of limits expected under the background-only hypothesis.

Exclusion limits at $95\%$ CL on the number of selected BSM events with $N_{\mathrm{b}} \geq 1$ as functions of $m_{\mathrm{Z'}}$ for $f_{2\mathrm{b}}=1$. The quantity $f_{2\mathrm{b}}$ is the fraction of BSM events passing the analysis selection that have at least two $\mathrm{b}$ quark jets. The solid black (dashed red) curve represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing $68~(95)\%$ of the distribution of limits expected under the background-only hypothesis.

Observed (solid) and median expected (dashed) exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution, marked by the dotted curves. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed (solid) and median expected (dashed) exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution, marked by the dotted curves. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The solid black (dashed red) curves represent the observed (median expected) exclusions. The dotted curves denote the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and dotted curves is (expected to be) excluded. The dotted curve for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$ lies beyond the displayed $|g_{\mathrm{Z'}}|$ range and is, therefore, not shown. The shaded blue area represents the region preferred from the global fit in JHEP 04 (2023) 033 at $95\%$ CL. The region above the green dash-dotted curve is incompatible at $95\%$ CL with the measurement of the mass difference between the mass eigenstates of the neutral $\mathrm{B}_\mathrm{s}$ mesons.

Observed exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$.

Expected exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The solid black (dashed red) curves represent the observed (median expected) exclusions. The dotted curves denote the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and dotted curves is (expected to be) excluded. The shaded blue area represents the region preferred from the global fit in JHEP 04 (2023) 033 at $95\%$ CL. The region above the green dash-dotted curve is incompatible at $95\%$ CL with the measurement of the mass difference between the mass eigenstates of the neutral $\mathrm{B}_\mathrm{s}$ mesons.

Observed exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$.

Expected exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The solid black (dashed red) curves represent the observed (median expected) exclusions. The dotted curves denote the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and dotted curves is (expected to be) excluded. The shaded blue area represents the region preferred from the global fit in JHEP 04 (2023) 033 at $95\%$ CL. The region above the green dash-dotted curve is incompatible at $95\%$ CL with the measurement of the mass difference between the mass eigenstates of the neutral $\mathrm{B}_\mathrm{s}$ mesons.

Observed exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$.

Expected exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The solid black (dashed red) curves represent the observed (median expected) exclusions. The dotted curves denote the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and dotted curves is (expected to be) excluded. The shaded blue area represents the region preferred from the global fit in JHEP 04 (2023) 033 at $95\%$ CL. The region above the green dash-dotted curve is incompatible at $95\%$ CL with the measurement of the mass difference between the mass eigenstates of the neutral $\mathrm{B}_\mathrm{s}$ mesons.

Observed exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$.

Expected exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{Z'}}|$--$m_{\mathrm{Z'}}$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for a fixed value of $\theta_{23}=0$. The solid black (dashed red) curve represents the observed (median expected) exclusion. The dotted curve denotes the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and the dotted curves is (expected to be) excluded. The shaded blue area represents the $|g_{\mathrm{Z'}}|$ range preferred from a global fit in JHEP 04 (2023) 033 at $95\%$ CL.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{Z'}}|$--$m_{\mathrm{Z'}}$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for a fixed value of $\theta_{23}=0$.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{Z'}}|$--$m_{\mathrm{Z'}}$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for a fixed value of $\theta_{23}=0$.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{Z'}}|$--$m_{\mathrm{Z'}}$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for a fixed value of $\theta_{23}=0$. The coupling values are shown at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid.

Post-fit distribution of the $M(ll)$ variable for the ultrahigh-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '2l soft' signal region of the '2/3l soft' analysis.

Post-fit distribution of the $M(ll)$ variable for the low-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '3l soft' signal region of the the '2/3l soft' analysis.

Post-fit distribution of the $M(ll)$ variable for the medium-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '3l soft' signal region of the the '2/3l soft' analysis.

Post-fit distribution of the $m_{\mathrm{T2}}(ll)$ variable for low-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '2l soft' signal region of the '2/3l soft' analysis.

Post-fit distribution of the $m_{\mathrm{T2}}(ll)$ variable for medium-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '2l soft' signal region of the '2/3l soft' analysis.

Post-fit distribution of the $m_{\mathrm{T2}}(ll)$ variable for high-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '2l soft' signal region of the '2/3l soft' analysis.

Post-fit distribution of the $m_{\mathrm{T2}}(ll)$ variable for ultrahigh-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '2l soft' signal region of the '2/3l soft' analysis.

2SS $\ell/{\geq}\,3\ell$ search: observed and expected yields across the SRs in category A, events with three light leptons of which at least two form an OSSF pair, after the requirement that the leading-lepton $p_{\mathrm{T}}$ be greater than 30 GeV is applied.

2SS $\ell/{\geq}\,3\ell$ search: observed and expected yields across the SRs of the '${\geq}\ 3\ell$' search in category B, events with three light leptons and no OSSF pair, after the requirement that the leading-lepton $p_{\mathrm{T}}$ be greater than 30 GeV is applied.

Wino-bino model: cross section limits in the model parameter space, for wino-like chargino-neutralino production in the WZ topology for the full parameter space.

Wino-bino model: cross section limits in the model parameter space, for wino-like chargino-neutralino production in the WZ topology for the compressed space.

Wino-bino model: cross section limits in the model parameter space, for wino-like chargino-neutralino production in the WH topology for the full parameter space.

Wino-bino model: cross section limits in the model parameter space, for wino-like chargino-neutralino production with mixed topology with equal branching fraction to WZ and WH.

Wino-bino model: exclusion contours from the individual and combined analyses targeting WZ topology for the full parameter space. For visualization of the exclusion contours, linear interpolation is employed to account for the limited granularity of the available signal samples.

Wino-bino model: exclusion contours from the individual and combined analyses targeting the corresponding compressed region. For visualization of the exclusion contours, linear interpolation is employed to account for the limited granularity of the available signal samples.

Wino-bino model: exclusion contours from the individual and combined analyses targeting the WH topology for the full parameter space. For visualization of the exclusion contours, linear interpolation is employed to account for the limited granularity of the available signal samples.

Wino-bino model: exclusion contours from the individual and combined analyses targeting combined contours for these two topologies. For visualization of the exclusion contours, linear interpolation is employed to account for the limited granularity of the available signal samples.

GMSB model: expected and observed cross section limits for the neutralino-neutralino production for the ZZ topology.

GMSB model: expected and observed cross section limits for the neutralino-neutralino production for the HH topology.

GMSB model: expected and observed cross section limits for the neutralino-neutralino production for the mixed topology with equal branching fraction to H and Z.

GMSB model: cross section limits for neutralino-neutralino production as a function of the NSLP mass and the branching fraction to the H boson for the combination of the searches.

GMSB model: exclusion limit for neutralino-neutralino production as a function of the NSLP mass and the branching fraction to the H boson for the combination of the searches along with the input searches. For visualization of the exclusion contours, linear interpolation is employed to account for the limited granularity of the available signal samples.

Cross section upper limit(s) in the mass plane of NLSP and LSP masses for the higgsino-bino model.

Mass plane cross section upper limit for direct slepton pair production, with observed and expected exclusion limits in the full mass plane from the combination.

Mass plane cross section upper limit for direct slepton pair production, with observed and expected exclusion limits in the compressed region from '2/3l' soft search.