Figure 1b, lower panel, subevent

Figure 2a, full-event

Figure 2b, subevent

Figure 3b

Figure 3c

Figure 3e

Figure 3f

Figure 4a

Figure 4b

Figure 5a, full-event

Figure 5b, subevent

Figure 5c, full-event

Figure 5d, subevent

Figure 6a

Figure 6b

Figure 7a, full-event

Figure 7b, subevent

Figure 7c, full-event

Figure 7d, subevent

Figure 8a, full-event

Figure 8b, subevent

Figure 9a, full-event

Figure 9b, subevent

Figure 10a, full-event

Figure 10b, subevent

Figure 11a, full-event

Figure 11b, full-event

Figure 11c, full-event

Figure 11d, subevent

Figure 11e, subevent

Figure 11f, subevent

Figure 12a, full-event

Figure 12b, subevent

Figure 13

Figure 14a

Figure 14b

Invariant mass distributions of $^3\hbox{He}+\pi^-$ (A), $^3\overline{\hbox{He}}+\pi^+$ (B), $^4\hbox{He}+\pi^-$ (C) and $^4\overline{\hbox{He}}+\pi^+$ (D). The solid bands mark the signal invariant mass regions. The obtained signal count ($N_{\rm Sig}$), background count ($N_{\rm Bg}$), and signal significance are listed in each panel.

(A) (anti)hypernuclei yields versus $L/\beta\gamma$. The vertical error bars represent the statistical error only.

(A) (anti)hypernuclei yields versus $L/\beta\gamma$. The vertical error bars represent the statistical error only.

Our measured (anti)hypernuclei lifetimes compared with world data and theoretical predictions (solid triangles). Error bars and boxes represent statistical and systematic uncertainties, respectively. Vertical solid lines with shaded regions show the average lifetime of (anti)hypernuclei and their corresponding uncertainties. The vertical gray line shows the lifetime of the free $\Lambda$.

Production yield ratios of various particles with the same baryon number. Results combining all collision systems in this work are shown by filled stars. Open stars show results with only U+U and Au+Au collisions, while quadrangular stars show results with only Zr+Zr and Ru+Ru collisions. Statistical uncertainties and systematic uncertainties are shown by vertical bars and boxes, respectively. Previous measurement results and thermal model predictions are also shown for comparison. Ratio Typle Index, 0:$\frac{{}^{3}\bar{He}}{{}^{3}He}$ 1:$\frac{{}^{3}\bar{He}}{{}^{3}He}\times\frac{\bar{p}}{p}$ 2:$\frac{{}^{4}\bar{He}}{{}^{4}He}$ 3:$\frac{{}^{3}_{\bar{\Lambda}}\bar{H}}{{}^{3}_{\Lambda}H}$ 4:$\frac{{}^{3}_{\bar{\Lambda}}\bar{H}}{{}^{3}_{\Lambda}H}\times\frac{\bar{p}}{p}$ 5:$\frac{{}^{4}_{\bar{\Lambda}}\bar{H}}{{}^{4}_{\Lambda}H}$ 6:$\frac{{}^{3}_{\Lambda}H}{{}^{3}He}$ 7:$\frac{{}^{4}_{\Lambda}H}{{}^{4}He}$ 8:$\frac{{}^{3}_{\bar{\Lambda}}\bar{H}}{^{3}\bar{He}}$ 9:$\frac{{}^{4}_{\bar{\Lambda}}\bar{H}}{^{4}\bar{He}}$

Production yield ratios of various particles with the same baryon number. Results combining all collision systems in this work are shown by filled stars. Open stars show results with only U+U and Au+Au collisions, while quadrangular stars show results with only Zr+Zr and Ru+Ru collisions. Statistical uncertainties and systematic uncertainties are shown by vertical bars and boxes, respectively. Previous measurement results and thermal model predictions are also shown for comparison. Ratio Typle Index, 0:$\frac{{}^{3}\bar{He}}{{}^{3}He}$ 1:$\frac{{}^{3}\bar{He}}{{}^{3}He}\times\frac{\bar{p}}{p}$ 2:$\frac{{}^{4}\bar{He}}{{}^{4}He}$ 3:$\frac{{}^{3}_{\bar{\Lambda}}\bar{H}}{{}^{3}_{\Lambda}H}$ 4:$\frac{{}^{3}_{\bar{\Lambda}}\bar{H}}{{}^{3}_{\Lambda}H}\times\frac{\bar{p}}{p}$ 5:$\frac{{}^{4}_{\bar{\Lambda}}\bar{H}}{{}^{4}_{\Lambda}H}$ 6:$\frac{{}^{3}_{\Lambda}H}{{}^{3}He}$ 7:$\frac{{}^{4}_{\Lambda}H}{{}^{4}He}$ 8:$\frac{{}^{3}_{\bar{\Lambda}}\bar{H}}{^{3}\bar{He}}$ 9:$\frac{{}^{4}_{\bar{\Lambda}}\bar{H}}{^{4}\bar{He}}$

Production yield ratios of various particles with the same baryon number. Results combining all collision systems in this work are shown by filled stars. Open stars show results with only U+U and Au+Au collisions, while quadrangular stars show results with only Zr+Zr and Ru+Ru collisions. Statistical uncertainties and systematic uncertainties are shown by vertical bars and boxes, respectively. Previous measurement results and thermal model predictions are also shown for comparison. Ratio Typle Index, 0:$\frac{{}^{3}\bar{He}}{{}^{3}He}$ 1:$\frac{{}^{3}\bar{He}}{{}^{3}He}\times\frac{\bar{p}}{p}$ 2:$\frac{{}^{4}\bar{He}}{{}^{4}He}$ 3:$\frac{{}^{3}_{\bar{\Lambda}}\bar{H}}{{}^{3}_{\Lambda}H}$ 4:$\frac{{}^{3}_{\bar{\Lambda}}\bar{H}}{{}^{3}_{\Lambda}H}\times\frac{\bar{p}}{p}$ 5:$\frac{{}^{4}_{\bar{\Lambda}}\bar{H}}{{}^{4}_{\Lambda}H}$ 6:$\frac{{}^{3}_{\Lambda}H}{{}^{3}He}$ 7:$\frac{{}^{4}_{\Lambda}H}{{}^{4}He}$ 8:$\frac{{}^{3}_{\bar{\Lambda}}\bar{H}}{^{3}\bar{He}}$ 9:$\frac{{}^{4}_{\bar{\Lambda}}\bar{H}}{^{4}\bar{He}}$

Production yield ratios of various particles with the same baryon number. Results combining all collision systems in this work are shown by filled stars. Open stars show results with only U+U and Au+Au collisions, while quadrangular stars show results with only Zr+Zr and Ru+Ru collisions. Statistical uncertainties and systematic uncertainties are shown by vertical bars and boxes, respectively. Previous measurement results and thermal model predictions are also shown for comparison. Ratio Typle Index, 0:$\frac{{}^{3}\bar{He}}{{}^{3}He}$ 1:$\frac{{}^{3}\bar{He}}{{}^{3}He}\times\frac{\bar{p}}{p}$ 2:$\frac{{}^{4}\bar{He}}{{}^{4}He}$ 3:$\frac{{}^{3}_{\bar{\Lambda}}\bar{H}}{{}^{3}_{\Lambda}H}$ 4:$\frac{{}^{3}_{\bar{\Lambda}}\bar{H}}{{}^{3}_{\Lambda}H}\times\frac{\bar{p}}{p}$ 5:$\frac{{}^{4}_{\bar{\Lambda}}\bar{H}}{{}^{4}_{\Lambda}H}$ 6:$\frac{{}^{3}_{\Lambda}H}{{}^{3}He}$ 7:$\frac{{}^{4}_{\Lambda}H}{{}^{4}He}$ 8:$\frac{{}^{3}_{\bar{\Lambda}}\bar{H}}{^{3}\bar{He}}$ 9:$\frac{{}^{4}_{\bar{\Lambda}}\bar{H}}{^{4}\bar{He}}$

Reconstruction efficiency as a function of $L/(\beta\gamma)$ obtained from the embedding Monte Carlo technique. Hypernuclei have stricter topological cuts than antihypernuclei to suppress knock-out $^3\hbox{He}$ and $^4\hbox{He}$, resulting in lower efficiencies.

Reconstruction efficiency as a function of $L/(\beta\gamma)$ obtained from the embedding Monte Carlo technique. Hypernuclei have stricter topological cuts than antihypernuclei to suppress knock-out $^3\hbox{He}$ and $^4\hbox{He}$, resulting in lower efficiencies.

H3L, antiH3L, H4L and antiH4L candidate invariant-mass distributions in different $L/\beta\gamma$ intervals.

H3L, antiH3L, H4L and antiH4L candidate invariant-mass distributions in different $L/\beta\gamma$ intervals.

H3L, antiH3L, H4L and antiH4L candidate invariant-mass distributions in different $L/\beta\gamma$ intervals.

H3L, antiH3L, H4L and antiH4L candidate invariant-mass distributions in different $L/\beta\gamma$ intervals.

$dN/d(L/\beta\gamma)$ as a function of $L/\beta\gamma$ for $\Lambda$ and $\bar{\Lambda}$, and exponential fits to obtain their lifetimes.

Efficiency corrected $p_{T}$ spectra for $^3\hbox{He}$, $^3\overline{\hbox{He}}$, \htl, and \ahtl. The spectra are in the phase space of $|y|<0.7$ with only MB-triggered events. The spectra are not normalized by the number of events. The lines represent the BW-function fits.

Efficiency corrected $p_{T}$ spectra for $^3\hbox{He}$, $^3\overline{\hbox{He}}$, \htl, and \ahtl. The spectra are in the phase space of $|y|<0.7$ with only MB-triggered events. The spectra are not normalized by the number of events. The lines represent the BW-function fits.

Transverse dilepton-$p_{T}^{miss}$ mass pre-fit distribution for selected events in SR2 of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton mass pre-fit distribution for selected events in SR3 of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton $p_T$ pre-fit distribution for selected events in SR3 of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dijet $p_T$ pre-fit distribution for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Lepton-dijet $p_T$ pre-fit distribution for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Maximum $p_T$ of the visible particles divided by the missing $p_T$, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between jets, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and missing $p_T$, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and the dijet system, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta R$ distance between the leptons, for selected events in the non-prompt validation region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the di-lepton system and the missing $p_T$, for selected events in the non-prompt validation region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and the di-jet system, for selected events in the non-prompt validation region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton-dijet system and the missing $p_T$, for selected events in the non-prompt validation region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Leading lepton $p_T$ pre-fit distribution for selected events in the top control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Trailing lepton $p_T$ pre-fit distribution for selected events in the top control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton mass pre-fit distribution for selected events in the Drell-Yan control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Missing $p_T$ pre-fit distribution for selected events in the Drell-Yan control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Missing $p_T$ pre-fit distribution for selected events in the $WW$ control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton $p_T$ pre-fit distribution for selected events in the $WW$ control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Minimum $p_T$ of the visible particles divided by the missing $p_T$, for selected events in the top control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between jets, for selected events in the top control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and missing $p_T$, for selected events in the $W$+jets control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and the dijet system, for selected events in the $W$+jets control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Unrolled $m_{ll}-m_{T}^{l min, p_{T}^{miss}}$ post-fit distributions in the di-leptonic channel in SR1 for the full data set. The histogram bins are spaced uniformly. Each group of five bins (from left to right) corresponds to $m_{T}^{l min, p_{T}^{miss}}$ distribution in a $m_{ll}$ region, placed in ascending order. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$.

Unrolled $m_{ll}-m_{T}^{l min, p_{T}^{miss}}$ post-fit distributions in the di-leptonic channel in SR2 for the full data set. The histogram bins are spaced uniformly. Each group of five bins (from left to right) corresponds to $m_{T}^{l min, p_{T}^{miss}}$ distribution in a $m_{ll}$ region, placed in ascending order. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$.

Unrolled $m_{ll}-m_{T}^{l min, p_{T}^{miss}}$ post-fit distributions in the di-leptonic channel in SR3 for the full data set. The histogram bins are spaced uniformly. Each group of five bins (from left to right) corresponds to $m_{T}^{l min, p_{T}^{miss}}$ distribution in a $m_{ll}$ region, placed in ascending order. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$.

Post-fit BDT distributions in the semi-leptonic channel for the full data set in the top control region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Post-fit BDT distributions in the semi-leptonic channel for the full data set in the $W$+jets control region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Post-fit BDT distributions in the semi-leptonic channel in the 2016 signal region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Post-fit BDT distributions in the semi-leptonic channel in the 2017+2018 signal region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

The distributions of the invariant mass of the dilepton-photon system after the signal region selection before the fit to data.The shaded band represents the Monte Carlo statistical unc

The distributions of the invariant mass of the transverse mass of the WW system after the signal region selection before the fit to data.The shaded band represents the Monte Carlo statistical unc

The unrolled two-dimensional mTWW-mll$\gamma$ distributions with category 0 jet after the fit to data. The data are compared to the sum of the signal and expected background. The black points with error bars represent the data and their statistical uncertainties, whereas the hatched bands represent the total uncertainties in the predictions.

The unrolled two-dimensional mTWW-mll$\gamma$ distributions with category 1more jets after the fit to data. The data are compared to the sum of the signal and expected background. The black points with error bars represent the data and their statistical uncertainties, whereas the hatched bands represent the total uncertainties in the predictions.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 350 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 400 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 450 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 500 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 550 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 600 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 650 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 700 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 750 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 800 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 850 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 900 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 950 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 1000 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Full experimental (statistical and systematic) correlations (in \%) of the normalized partial decay rates over the $\overline{B}^0\to D^{*+} e^- \bar\nu_e$ and $\overline{B}^0\to D^{*+} \mu^- \bar\nu_\mu$ decays.

Measured partial decay rates $\Delta\Gamma$ (in units of $10^{-15}$ GeV) using the matrix inversion unfolding method

Average of normalized decay rates that are determined using the matrix inversion unfolding method

Full experimental (statistical and systematic) correlations (in \%) of the partial decay rates determined with the matrix inversion unfolding method for the $\overline{B}^0\to D^{*+} e^- \bar\nu_e$ and $\overline{B}^0\to D^{*+} \mu^- \bar\nu_\mu$ decays.

Full experimental (statistical and systematic) correlations (in \%) of the normalized partial decay rates determined using the matrix inversion unfolding method

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.

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

Observed cross 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

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

Observed cross 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

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

Observed cross 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

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

Observed cross 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

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

Observed cross 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

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

Observed cross 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

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

Observed cross 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

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

Observed cross 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:pp -> chi1p-chi1m / chi10-chi1pmHiggsino model: ctau calculated as eq. 12 in arXiv:1703.09675.Individual events weighted according to the resulting decay lengths.

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:pp -> chi1p-chi1m / chi10-chi1pmPure Wino model: ctau calculated as eq. 61 in arXiv:1410.4549.Individual events weighted according to the resulting decay lengths.

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:pp -> chi20-chi1pm / chi1p-chi1m / chi10-chi1pm; dm(chi20 - chi10) = 2dm(chi1pm - chi10)Higgsino model: ctau calculated as eq. 12 in arXiv:1703.09675.Individual events weighted according to the resulting decay lengths.

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:pp -> chi20-chi1pm / chi1p-chi1m / chi10-chi1pm; dm(chi20 - chi10) = 2dm(chi1pm - chi10)Pure Higgsino model: ctau calculated as eq. 12 in arXiv:1703.09675.Individual events weighted according to the resulting decay lengths.

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

Efficiency as a function of transverse (x-y) decay length of the chargino for short track selection in Run 2.

Efficiency as a function of transverse (x-y) decay length of the chargino for long track selection in Run 2.

Efficiency as a function of transverse (x-y) decay length of the chargino for combined track selection in Run 2.

Cut flow for model T2tb (m(~b)=1 TeV, m(~chi10)=900 GeV, c#tau=10 cm), showing events passing each cut stage.

Cut flow for model T1btbt (m(~g)=1.5 TeV, m(~chi10)=1.1 TeV, c#tau=200 cm), showing events passing each cut stage.

Cut flow for model T2tb (m(~b)=1 TeV, m(~chi10)=900 GeV, c#tau=200 cm), showing events passing each cut stage.

Cut flow for model T2tb (m(~b)=1 TeV, m(~chi10)=900 GeV, c#tau=200 cm), showing events passing each cut stage.

Cut flow for model T1btbt (m(~g)=1.5 TeV, m(~chi10)=1.1 TeV, c#tau=200 cm), showing events passing each cut stage.

Cut flow for model T1btbt (m(~g)=1.5 TeV, m(~chi10)=1.1 TeV, c#tau=10 cm), showing events passing each cut stage.

Cut flow for model T2tb (m(~b)=1 TeV, m(~chi10)=900 GeV, c#tau=10 cm), showing events passing each cut stage.

Cut flow for model T1btbt (m(~g)=1.5 TeV, m(~chi10)=1.1 TeV, c#tau=10 cm), showing events passing each cut stage.