$\Delta\phi$ projection of balance function in low $p_{T}$ in 0-10% centrality

$\Delta\phi$ projection of balance function in low $p_{T}$ in 30-40% centrality

1D $\Delta\eta$ projection of balance function in low $p_{T}$ in 70-80% centrality

$\Delta\eta$ projection of balance function in low $p_{T}$ in 0-40 $N_{trk}$

$\Delta\eta$ projection of balance function in low $p_{T}$ in 120-150 $N_{trk}$

$\Delta\eta$ projection of balance function in low $p_{T}$ in 270-300 $N_{trk}$

$\Delta\phi$ projection of balance function in low $p_{T}$ in 0-40 $N_{trk}$

$\Delta\phi$ projection of balance function in low $p_{T}$ in 120-150 $N_{trk}$

$\Delta\phi$ projection of balance function in low $p_{T}$ in 270-300 $N_{trk}$

The width of balance function in $<|\Delta\eta|>$ in low $p_{T}$ for PbPb collisions

Balance function width in low $p_{T}$ PbPb relative $<|\Delta\eta|>/<|\Delta\eta|>_{N_{ch}<65}$

The width of balance function in $<|\Delta\eta|>$ in low $p_{T}$ for pPb collisions

Balance function width in low $p_{T}$ in pPb relative $<|\Delta\eta|>/<|\Delta\eta|>_{N_{ch}<24}$

The width of the balance function in $<|\Delta\phi|>$ in low $p_{T}$ for PbPb collisions

Balance function width in low $p_{T}$ in PbPb for relative $<|\Delta\phi|>/<|\Delta\phi|>_{N_{ch}<65}$

The width of the balance function in $<|\Delta\phi|>$ in low $p_{T}$ for pPb collisions

Balance function width in low $p_{T}$ in pPb for relative $|\Delta\phi|/<|\Delta\phi|>_{N_{ch}<24}$

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in pPb for 0-40 $N_{trk}$

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in pPb for 120-150 $N_{trk}$

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in pPb for 270-300 $N_{trk}$

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in pPb for $N_{trk}$ 0-40

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in pPb for $N_{trk}$ 120-150

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in pPb for $N_{trk}$ 270-300

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in pPb for $N_{trk}$ 0-40

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in pPb for $N_{trk}$ 120-150

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in pPb for $N_{trk}$ 270-300

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for $N_{trk}$ 0-40

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for $N_{trk}$ 120-150

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for $N_{trk}$ 270-300

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for low $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for high $p_{T}$

Average width of the balance function in $\Delta \phi$ with $p_{T}$ for low $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for high $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for low $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for high $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for low $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for high $p_{T}$

The nuclear modification factor $\mathrm{R_{AA}}$ versus $p_{\mathrm{T}}$ for prompt $\Lambda^+_c$ production in centrality regions of 0-90, 0-10, 10-30, 30-50 and 50-90% in PbPb collisions. The global uncertainty includes the uncertainties for the luminosity of pp collisions, number of MB events in PbPb collisions, and tracking efficiency. The global uncertainty for $\mathrm{R_{AA}}$ is 16.5%.

The ratio of the production cross sections of prompt $\Lambda^+_c$ to prompt $\mathrm{D_0}$ versus $p_{\mathrm{T}}$ from pp collisions. The global normalization uncertainty is 6.6%.

The ratio of the $\mathrm{T_{AA}}$-scaled prompt $\Lambda^+_c$ yields to prompt $\mathrm{D_0}$ versus $p_{\mathrm{T}}$ from PbPb collisions for 0-90 and 0-10% centrality classes. The global normalization uncertainty is 7.3%.

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.

Reconstructed mass of the Z' candidate in CR1, which consists of flavor sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Pre-fit (post-fit) results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR1, which consists of flavor sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Pre-fit (post-fit) results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR1, which consists of flavor sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Pre-fit (post-fit) results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

The observed and expected 95\% CL upper limits on the product of \Zp boson cross section times branching fraction are shown for the case of ${m}_{N}$ = ${m}_{Z'}/4$ (upper row) and ${m}_{N}$ = 100 GeV (lower row) for dielectron (left column) and dimuon (right column) channels. The green and yellow bands indicate the 68\% and 95\% CL regions around the expected limit.

The observed and expected 95\% CL upper limits on the product of \Zp boson cross section times branching fraction are shown for the case of ${m}_{N}$ = ${m}_{Z'}/4$ (upper row) and ${m}_{N}$ = 100 GeV (lower row) for dielectron (left column) and dimuon (right column) channels. The green and yellow bands indicate the 68\% and 95\% CL regions around the expected limit.

The observed and expected 95\% CL upper limits on the product of \Zp boson cross section times branching fraction are shown for the case of ${m}_{N}$ = ${m}_{Z'}/4$ (upper row) and ${m}_{N}$ = 100 GeV (lower row) for dielectron (left column) and dimuon (right column) channels. The green and yellow bands indicate the 68\% and 95\% CL regions around the expected limit.

The observed and expected 95\% CL upper limits on the product of \Zp boson cross section times branching fraction are shown for the case of ${m}_{N}$ = ${m}_{Z'}/4$ (upper row) and ${m}_{N}$ = 100 GeV (lower row) for dielectron (left column) and dimuon (right column) channels. The green and yellow bands indicate the 68\% and 95\% CL regions around the expected limit.

Observed and expected exclusion regions at 95% CL in the 2D phase space of ${m}_{Z'}$ vs. ${m}_{N}$ for dielectron (left) and dimuon (right) channels. One standard deviation (s.d.) limits are also shown. Both opposite- and same-sign (OS and SS) dileptons from decays of heavy neutrino pairs are considered.

Observed and expected exclusion regions at 95% CL in the 2D phase space of ${m}_{Z'}$ vs. ${m}_{N}$ for dielectron (left) and dimuon (right) channels. One standard deviation (s.d.) limits are also shown. Both opposite- and same-sign (OS and SS) dileptons from decays of heavy neutrino pairs are considered.

The Levy scale parameter $R$ versus the transverse mass $m_{\mathrm{T}}$ in different centrality classes for negatively and positively charged hadron pairs.

The $1/R^2$ distribution vs. transverse mass $m_{\mathrm{T}}$ in different centrality classes for negatively and positively charged hadron pairs.

The linear-fit parameter $A$ versus $\langle N_{\text{part}}\rangle$ for positively and negatively charged hadron pairs. The bin range on $\langle N_{\text{part}}\rangle$ represents the 68% confidence interval associated to its uncertainty.

The linear-fit parameter $B$ versus $\langle N_{\text{part}}\rangle$ for positively and negatively charged hadron pairs. The bin range on $\langle N_{\text{part}}\rangle$ represents the 68% confidence interval associated to its uncertainty.

The Levy scale parameter $R$ versus $\langle N_{\text{part}}\rangle^{1/3}$ in different $m_{\mathrm{T}}$ classes for negatively and positively charged hadron pairs. The bin range on $\langle N_{\text{part}}\rangle^{1/3}$ represents the 68% confidence interval associated to its uncertainty.

The Levy stability index $\alpha$ versus the transverse mass $m_{\mathrm{T}}$ in different centrality classes for negatively and positively charged hadron pairs.

The average Levy stability index $\langle \alpha \rangle$ versus $\langle N_{\text{part}}\rangle$ for both positively and negatively charged hadron pairs. The bin range on $\langle N_{\text{part}}\rangle$ represents the 68% confidence interval associated to its uncertainty.

The correlation strength $\lambda$ versus the transverse mass $m_{\mathrm{T}}$ in different centrality classes for negatively and positively charged hadron pairs.

The correlation strength (re-scaled with the square of the pion fraction) $\lambda ^*$ versus the transverse mass $m_{\mathrm{T}}$ in different centrality classes for positively and negatively charged hadron pairs.

The $m_{e\mu}$ distribution of the observed data is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in red solid line, and the background component of the fit in red dotted line. Events and fit in each category are weighted by $S/(S+B)$. The one and two standard deviation uncertainty bands of the background component are shown in green and yellow. The lower panel shows the residuals after subtracting the background component of the fit from data.

Observed (expected) 95% confidence level upper limits on $\sigma(p p \to X(146) \to e \mu)$ for each individual analysis category (as shown in the left axis label) and for the combination of all analysis categories assuming the relative SM-like production cross sections of the ggH and VBF production modes.

The $m_{e\mu}$ distribution of the observed data in the category ggH cat 0 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category ggH cat 1 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category ggH cat 2 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category ggH cat 3 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category VBF cat 0 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category VBF cat 1 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

Best fit of $\sigma(p p \to X(146) \to e \mu)$ for each analysis category (black point) compared to the best fit for the combination of all analysis categories (blue line). The one standard deviation uncertainty of the per-category best fit is shown in red line and the one standard deviation uncertainty of the combined best fit is shown in green band.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.0.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.1.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.2.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.3.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.4.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.5.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.6.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.7.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.8.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.9.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 1.0.

The $v_{3}$ values as a function of centrality for prompt and nonprompt J/$\psi$ mesons.The kinematic range is 6.5 $< p_{\text{T}} <$ 50 GeV/c and $|y| < 2.4$.

The $v_{2}$ values of prompt J/$\psi$ mesons as a function of $p_{\text{T}}$ in the 10–60% centrality range. The results for 3 $< p_{\text{T}} <$ 6.5 and 6.5 $< p_{\text{T}} <$ 50 GeV/c are studied in the rapidity range of 1.6 $< |y| <$ 2.4 and $|y| <$ 2.4, respectively.

The $v_{2}$ values as functions of $p_{\text{T}}$ for prompt J/$\psi$ and prompt $\psi$(2S) mesons in the 10–60% centrality range. The results for $p_{\text{T}}$ < 6.5 and $p_{\text{T}}$ > 6.5 GeV/c are studied in the rapidity range of 1.6 $< |y| <$ 2.4 and $|y| <$ 2.4, respectively.

The $v_{2}$ values as functions of centrality for prompt J/$\psi$ and prompt $\psi$(2S) mesons. The kinematic range is 6.5 $< p_{\text{T}} <$ 50 GeV/c and $|y| <$ 2.4.

The $v_{3}$ values as functions of $p_{\text{T}}$ for prompt J/$\psi$ and prompt $\psi$(2S) mesons in the 10–60% centrality range. The results for $p_{\text{T}}$ < 6.5 and $p_{\text{T}}$ > 6.5 GeV/c are studied in the rapidity range of 1.6 $< |y| <$ 2.4 and $|y| <$ 2.4, respectively.

The $v_{3}$ values as functions of centrality for prompt J/$\psi$ and prompt $\psi$(2S) mesons. The kinematic range is 6.5 $< p_{\text{T}} <$ 50 GeV/c and $|y| <$ 2.4.

Measured min-$d_{xy}$ distribution in the 2-match category, after requiring the $d_{xy}$ muon to pass the isolation requirement $I_{\mathrm{rel}}^{\mathrm{PF}} <0.25$ (i. e., the B and D bins of the ABCD plane). Overlaid with a red histogram is the background predicted from the region of the ABCD plane failing the same requirement (the A and C bins), as well as three signal benchmark hypotheses (as defined in the legends), assuming $\alpha_D = \alpha_{\mathrm{EM}}$ (the fine-structure constant). The red hatched bands correspond to the background prediction uncertainty. The last bin includes the overflow.

Two-dimensional exclusion surfaces for $\Delta = 0.1 \, m_1$ as a function of the DM mass $m_1$ and the signal strength $y$, with $m_{A'} = 3 \, m_1$. Filled histograms denote observed limits on $\sigma(\mathrm{pp} \to A' \to \chi_2 \chi_1) \, \mathcal{B}(\chi_2 \to \chi_1 \mu^+ \mu^-)$. Solid (dashed) curves denote the observed (expected) exclusion limits at 95% CL, with 68% CL uncertainty bands around the expectation. Regions above the curves are excluded, depending on the $\alpha_D$ hypothesis: $\alpha_{\mathrm{D}} = \alpha_{\mathrm{EM}}$ (dark blue) or 0.1 (light magenta).

Two-dimensional exclusion surfaces for $\Delta = 0.4 \, m_1$ as a function of the DM mass $m_1$ and the signal strength $y$, with $m_{A'} = 3 \, m_1$. Filled histograms denote observed limits on $\sigma(\mathrm{pp} \to A' \to \chi_2 \chi_1) \, \mathcal{B}(\chi_2 \to \chi_1 \mu^+ \mu^-)$. Solid (dashed) curves denote the observed (expected) exclusion limits at 95% CL, with 68% CL uncertainty bands around the expectation. Regions above the curves are excluded, depending on the $\alpha_D$ hypothesis: $\alpha_{\mathrm{D}} = \alpha_{\mathrm{EM}}$ (dark blue) or 0.1 (light magenta).