$\pi^0$ spectra from figure 2a from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ spectra from figure 2a from 60-80% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 2b from minimum bias U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 2b from 0-20% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 2b from 20-40% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 2b from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 2b from 60-80% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

Ratio of $\pi^0$ yield to the fit from figure 2c from minimum bias U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\pi^0$ yield to the fit from figure 2c from 0-20% U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\pi^0$ yield to the fit from figure 2c from 60-80% U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\eta$ yield to the fit from figure 2d from minimum bias U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\eta$ yield to the fit from figure 2d from 0-20% U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\eta$ yield to the fit from figure 2d from 40-60% U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

$\eta/\pi^0$ ratio from figure 3 from minimum bias U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from 0-20% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from 20-40% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from 60-80% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from minimum bias Au+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\pi^0$ $R_{AA}$ from figure 4a from 0-20% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\eta$ $R_{AA}$ from figure 4a from 0-20% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\pi^0$ $R_{AA}$ from figure 4b from 20-40% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\eta$ $R_{AA}$ from figure 4b from 20-40% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\pi^0$ $R_{AA}$ from figure 4c from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\eta$ $R_{AA}$ from figure 4c from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\pi^0$ $R_{AA}$ from figure 4d from minimum bias U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\eta$ $R_{AA}$ from figure 4d from minimum bias U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\pi^0$ $R_{AA}$ from figure 5a from 0-20% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 2 set.

$\pi^0$ $R_{AA}$ from figure 5a from 0-5% Au+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AA}$ from figure 5a from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 2 set.

$\pi^0$ $R_{AA}$ from figure 5a from 40-50% Au+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ integrated $R_{AA}$ for $p_T>5$ GeV/$c$ in U+U collisions as a function of centrality from figure 6a. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\eta$ integrated $R_{AA}$ for $p_T>5$ GeV/$c$ in U+U collisions as a function of centrality from figure 6a. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\pi^0$ integrated $R_{AA}$ for $p_T>5$ GeV/$c$ in Au+Au collisions as a function of centrality from figure 6a. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ integrated $R_{AA}$ for $p_T>5$ GeV/$c$ in U+U collisions as a function of centrality from figure 6b. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 2 set.

Number of participant nucleons of in dep U$+$U collisions

Number of participant nucleons of in dep U$+$U collisions

Number of participant nucleons of in dep Au$+$Au collisions

$\lambda_{\phi}$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with NRQCD predictions in the Helicity and Collins-Soper frames.

$\lambda_{\theta \phi}$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with NRQCD predictions in the Helicity and Collins-Soper frames.

$\lambda_{\theta \phi}$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with NRQCD predictions in the Helicity and Collins-Soper frames.

$\bar{\lambda}$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with CSM and NRQCD predictions in the Helicity frame.

$\bar{\lambda}$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with CSM and NRQCD predictions in the Helicity frame.

Invariant $F$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with CSM and NRQCD predictions in the Helicity frame.

Invariant $F$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with CSM and NRQCD predictions in the Helicity frame.

$b\bar{b}$ invariant yield measured via B meson decay to like-sign dimuons as a function of pair $p_{T}$. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

Integrated away-side $\gamma_{dir}$-h per-trigger $I_{AA}$ and $I_{dA}$ of Au+Au, d+Au, and p+p, as a function of $\xi$.

$I_{AA}$ vs $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7 GeV/$c$, 7-9 GeV/$c$, and 9-12 GeV/$c$.

$I_{AA}$ vs $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7 GeV/$c$, 7-9 GeV/$c$, and 9-12 GeV/$c$.

$I_{AA}$ vs $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7 GeV/$c$, 7-9 GeV/$c$, and 9-12 GeV/$c$.

$I_{AA}$ as a function of $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7, 7-9, and 9-12 GeV/$c$. Three away-side integration ranges are chosen to calculate the per-trigger yield and the corresponding $I_{AA}$, $|\Delta\phi-\pi|<\pi/2$, $|\Delta\phi-\pi|<\pi/3$, and $|\Delta\phi-\pi|<\pi/6$.

$I_{AA}$ as a function of $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7, 7-9, and 9-12 GeV/$c$. Three away-side integration ranges are chosen to calculate the per-trigger yield and the corresponding $I_{AA}$, $|\Delta\phi-\pi|<\pi/2$, $|\Delta\phi-\pi|<\pi/3$, and $|\Delta\phi-\pi|<\pi/6$.

$I_{AA}$ as a function of $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7, 7-9, and 9-12 GeV/$c$. Three away-side integration ranges are chosen to calculate the per-trigger yield and the corresponding $I_{AA}$, $|\Delta\phi-\pi|<\pi/2$, $|\Delta\phi-\pi|<\pi/3$, and $|\Delta\phi-\pi|<\pi/6$.

Ratios of $I_{AA}$ as a function of direct photon $p_T$ for three different away-side integration ranges.

Measured $I_{AA}$ for direct photon $p_T$ of 5-7, 7-9, and 9-12 GeV/$c$, as a function of $\xi$, are compared with theoretical model calculations.

Measured $I_{AA}$ for direct photon $p_T$ of 5-7, 7-9, and 9-12 GeV/$c$, as a function of $\xi$, are compared with theoretical model calculations.

Measured $I_{AA}$ for direct photon $p_T$ of 5-7, 7-9, and 9-12 GeV/$c$, as a function of $\xi$, are compared with theoretical model calculations.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\nu\nu jj$ SR.

The expected exclusion limit contours at 95% CL for the $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ production in the $m(\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2})$-$m(\tilde{\chi}^{0}_{1})$ plane.

The +1 $\sigma$ expected exclusion limit contours at 95% CL for the $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ production in the $m(\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2})$-$m(\tilde{\chi}^{0}_{1})$ plane.

The -1 $\sigma$ expected exclusion limit contours at 95% CL for the $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ production in the $m(\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2})$-$m(\tilde{\chi}^{0}_{1})$ plane.

The +1 $\sigma$ observed exclusion limit contours at 95% CL for the $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ production in the $m(\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2})$-$m(\tilde{\chi}^{0}_{1})$ plane.

The -1 $\sigma$ observed exclusion limit contours at 95% CL for the $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ production in the $m(\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2})$-$m(\tilde{\chi}^{0}_{1})$ plane.

The observed and expected exclusion limit at 95% CL for the $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ production in the $m(\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2})$-$m(\tilde{\chi}^{0}_{1})$ plane.

Expected and observed 95% CL exclusion upper limits on the higgsino production ($\tilde{\chi}\tilde{\chi} \equiv \tilde{\chi}^{0}_{1}\tilde{\chi}^{0}_{2}, \tilde{\chi}^{0}_{1}\tilde{\chi}^{\pm}_{1}, \tilde{\chi}^{0}_{2}\tilde{\chi}^{\pm}_{1}, \tilde{\chi}^{\pm}_{1}\tilde{\chi}^{\mp}_{1}$) cross-section in the signal of $h \tilde{G} h \tilde{G} $ as a function of the higgsino mass.

The distribution of $S_{E_\mathrm{T}^{\mathrm{miss}}}$ after the selection of diphoton candidates with $120~GeV < m_{\gamma\gamma} < 130~GeV$. Expected distributions are shown for the $\tilde{\chi}^{\pm}_{1} \tilde{\chi}^{0}_{2} \to W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1} $ signal with $m(\tilde{\chi}^{\pm}_{1} \tilde{\chi}^{0}_{2})=200~GeV$ and $m(\tilde{\chi}^{0}_{1})=0.5~GeV$, and the $h\tilde{G} h\tilde{G}$ signal with $m(\tilde{\chi}^{0}_{1})=150~GeV$ and $m(\tilde{G})=1~MeV$. These overlaid signal points are representative of the model kinematics. The sum in quadrature of the MC statistical and experimental systematic uncertainties in the total background is shown as the hatched bands, while the theoretical uncertainties in the background normalisation are not included. The $t\bar{t}\gamma$ and $t\bar{t}\gamma\gamma$ processes have a negligible contribution and are not represented. Overflow events are included in the rightmost bin. The lower panel shows the ratio of the data to the background.

Acceptance and efficiency for the $W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1}$ signal in each category.

Acceptance and efficiency for the $h \tilde{G} h \tilde{G}$ signal in each category.

Cut flow for the $W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1}$ signal sample at $m(\tilde{\chi}^{\pm}_{1} \tilde{\chi}^{0}_{2})=200~GeV$, $m(\tilde{\chi}^{0}_{1})=0.5~GeV$, with 260000 entries.

Cut flow for the $h \tilde{G} h \tilde{G}$ signal sample at $m(\tilde{\chi}^{0}_{1})=150~GeV$, $m(\tilde{G})=1~MeV$, with 50000 entries.

Event yields in the range 120 <$m\gamma\gamma$< 130 GeV for data, signal models, the SM Higgs boson background and non-resonant background in each analysis category after the fit.

Event yields in the range 120 <$m\gamma\gamma$< 130 GeV for data, signal models, the SM Higgs boson background and non-resonant background in each analysis category before the fit.

Acceptance for the $W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1}$ signals in the leptonic channel.

Efficiency for the $W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1}$ signals in the leptonic channel.

Acceptance for the $W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1}$ signals in the hadronic channel.

Efficiency for the $W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1}$ signals in the hadronic channel.

Acceptance for the $W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1}$ signals in the rest channel.

Efficiency for the $W^{\pm} \tilde{\chi}^{0}_{1} h \tilde{\chi}^{0}_{1}$ signals in the rest channel.

Acceptance for the $h \tilde{G} h \tilde{G}$ signals in the each channel.

Efficiency for the $h \tilde{G} h \tilde{G}$ signals in the each channel.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 0.33 < ln(R/#DeltaR) < 0.67.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 0.67 < ln(R/#DeltaR) < 1.00.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 1.00 < ln(R/#DeltaR) < 1.33.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 1.33 < ln(R/#DeltaR) < 1.67.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 1.67 < ln(R/#DeltaR) < 2.00.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 2.00 < ln(R/#DeltaR) < 2.33.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 2.33 < ln(R/#DeltaR) < 2.67.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 2.67 < ln(R/#DeltaR) < 3.00.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 3.00 < ln(R/#DeltaR) < 3.33.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 3.33 < ln(R/#DeltaR) < 3.67.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 3.67 < ln(R/#DeltaR) < 4.00.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 4.00 < ln(R/#DeltaR) < 4.33.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 0.69 < ln(1/z) < 0.97.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 0.97 < ln(1/z) < 1.25.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 1.25 < ln(1/z) < 1.52.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 1.52 < ln(1/z) < 1.80.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 1.80 < ln(1/z) < 2.08.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 2.08 < ln(1/z) < 2.36.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 2.36 < ln(1/z) < 2.63.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 2.63 < ln(1/z) < 2.91.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 2.91 < ln(1/z) < 3.19.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 3.19 < ln(1/z) < 3.47.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 3.47 < ln(1/z) < 3.74.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 3.74 < ln(1/z) < 4.02.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 4.02 < ln(1/z) < 4.30.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 4.30 < ln(1/z) < 4.57.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 4.57 < ln(1/z) < 4.85.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 4.85 < ln(1/z) < 5.13.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 5.13 < ln(1/z) < 5.41.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 5.41 < ln(1/z) < 5.68.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 5.68 < ln(1/z) < 5.96.

The summed covariance matrix of all systematic and statistical uncertainties associated with the measurement in bins of $\ln{(1/z)} \times \ln{(R/\Delta R)}$.

The summed covariance matrix of all statistical uncertainties associated with the measurement in bins of $\ln{(1/z)} \times \ln{(R/\Delta R)}$.

Track multiplicity $n_{Tracks}$ for preselected DVs in muon-triggered events with at least one muon passing the full selection. Along with the data shown with black markers, the stacked filled histograms represent the background estimates, and predictions for signal scenarios are overlaid with dashed lines. The errors include statistical and systematic uncertainties and are indicated by hatched bands. The DV full selection requirements, $n_{Tracks} \geq 3$ and $m_{DV} > 20$ GeV are visualized with a black arrow.

Invariant mass $m_{DV}$ for the highest-mass preselected DV with at least three associated tracks in MET-triggered events with at least one muon passing the full selection. Along with the data shown with black markers, the stacked filled histograms represent the background estimates, and predictions for signal scenarios are overlaid with dashed lines. The errors include statistical and systematic uncertainties and are indicated by hatched bands. The DV full selection requirements, $n_{Tracks} \geq 3$ and $m_{DV} > 20$ GeV are visualized with a black arrow.

Invariant mass $m_{DV}$ for the highest-mass preselected DV with at least three associated tracks in muon-triggered events with at least one muon passing the full selection. Along with the data shown with black markers, the stacked filled histograms represent the background estimates, and predictions for signal scenarios are overlaid with dashed lines. The errors include statistical and systematic uncertainties and are indicated by hatched bands. The DV full selection requirements, $n_{Tracks} \geq 3$ and $m_{DV} > 20$ GeV are visualized with a black arrow.

The observed event yields in the control, validation and signal regions are shown for the MET Trigger selections, along with the predicted background yields. The bottom panel shows the ratio of observed events to the total background yields. The errors represent the total uncertainty of the backgrounds prediction, including the statistical and systematic uncertainties added in quadrature.

The observed event yields in the control, validation and signal regions are shown for the Muon Trigger selections, along with the predicted background yields. The bottom panel shows the ratio of observed events to the total background yields. The errors represent the total uncertainty of the backgrounds prediction, including the statistical and systematic uncertainties added in quadrature.

Expected exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Expected (1 sigma band) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Expected (2 sigma band) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed (+1 sigma) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed (-1 sigma) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Exclusion limits on the production cross section as a function of m($\tilde{t}$) are shown for several values of $\tau(\tilde{t})$ along with the nominal signal production cross section and its theoretical uncertainty.

Parameterized event selection efficiencies for the $E_{T}^{miss}$ Trigger SR. The event-level efficiencies for each SR are extracted for all events passing the acceptance of the corresponding SR.

Parameterized event selection efficiencies for the Muon Trigger SR. The event-level efficiencies for each SR are extracted for all events passing the acceptance of the corresponding SR.

Parameterized muon-level reconstruction efficiencies as a function of the muon $p_{T}$ and $d_{0}$. The muon-level efficiencies are extracted using muons passing the muon acceptance criteria.

Parameterized vertex-level reconstruction efficiencies as a function of the radial position of the truth vertex. The efficiency is calculated independent of the muons originating from this truth vertex.

Parameterized vertex-level reconstruction efficiencies as a function of the radial position of the truth vertex. The efficiency is calculated only for truth vertices which have a muon originating from them which is matched to a reconstructed muon.

The $p_{T}$ distribution of all muons originating from LLP decays in the samples used to calculate and validate the efficiencies.

The invariant mass and multiplicity of selected decay products of all truth vertices used in the calculation and validation of the reconstructed efficiencies.