Upper limits at 95% CL on $\sigma$B(pp→X→Y(bb)H) for combination as a function of m$_Y$.

Expected upper limits at 95% CL on the X → HH cross section as a function of mX from the combination of the three channels bbγγ, bbττ, and bbbb for an integrated luminosity of 3000 fb−1 with different assumptions on the systematic uncertainties.

Expected upper limits on the X → YH cross section times B(Y → bb) for mX ≤ 1000 GeV at 95 % CL as a function of m$_X$ and m$_Y$ from the combination of all considered channels projected to an integrated luminosity of 3000 fb$^{-1}$. Systematic uncertainties according to the ”S2” scenario are assumed.

Example description

Example description

Example description

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.33, 0.33, 0.33) and in the Majorana scenario.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.0, 1.0, 0.0) and in the Dirac-like scenario.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.0, 0.5, 0.5) and in the Dirac-like scenario.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.5, 0.5, 0.0) and in the Dirac-like scenario.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.33, 0.33, 0.33) and in the Dirac-like scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 1.0 GeV in the Majorana scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 1.5 GeV in the Majorana scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 2.0 GeV in the Majorana scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 1.0 GeV in the Dirac-like scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 1.5 GeV in the Dirac-like scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 2.0 GeV in the Dirac-like scenario.

Measured natural width

The measured inclusive ratio of production cross sections

Comparison of the observed (black points) and expected (histograms) numbers of events in nonoverlapping $m^{corr}_{\mu\mu}$ intervals in the STA-STA dimuon category, in the signal region optimized for the RPV SUSY model. Yellow and green stacked filled histograms represent mean expected background contributions from QCD and DY, respectively, while statistical uncertainties in the total expected background are shown as hatched histograms. Signal contributions expected from simulated signals indicated in the legends are shown in red and blue. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Distribution of min($d_0 / \sigma_{d_0}$) for TMS-TMS dimuons with $|\Delta\Phi| < \pi/30$, for events in all mass intervals combined, for both the validation (min($d_0 / \sigma_{d_0}$) < 6) and signal (min($d_0 / \sigma_{d_0}$) > 6) regions. The number of observed events (black circles) is overlaid with the stacked histograms showing the expected numbers of QCD (yellow) and DY (green) background events. Statistical uncertainties in the total expected background are shown as hatched histograms. Signal contributions expected from simulated signals indicated in the legends are shown in red and blue. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. Events are required to satisfy all nominal selection criteria with the exception of the $d_0 / \sigma_{d_0}$ requirement. The last bin includes events in the histogram overflow.

Distribution of min($d_0 / \sigma_{d_0}$) for TMS-TMS dimuons with $|\Delta\Phi| < \pi/4$, for events in all mass intervals combined, for both the validation (min($d_0 / \sigma_{d_0}$) < 6) and signal (min($d_0 / \sigma_{d_0}$) > 6) regions. The number of observed events (black circles) is overlaid with the stacked histograms showing the expected numbers of QCD (yellow) and DY (green) background events. Statistical uncertainties in the total expected background are shown as hatched histograms. Signal contributions expected from simulated signals indicated in the legends are shown in red and blue. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. Events are required to satisfy all nominal selection criteria with the exception of the $d_0 / \sigma_{d_0}$ requirement. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in bins of $m^{corr}_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the RPV SUSY model. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m^{corr}_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: 6-10. Hatched histograms show statistical uncertainties in the total expected background. Contributions expected from signal events predicted by the RPV SUSY model with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in bins of $m^{corr}_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the RPV SUSY model. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m^{corr}_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: 10-20. Hatched histograms show statistical uncertainties in the total expected background. Contributions expected from signal events predicted by the RPV SUSY model with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in bins of $m^{corr}_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the RPV SUSY model. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m^{corr}_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: >20. Hatched histograms show statistical uncertainties in the total expected background. Contributions expected from signal events predicted by the RPV SUSY model with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in bins of $m_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the HAHM. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: 6-10. Hatched histograms show statistical uncertainties in the total expected background. Signal contributions expected from simulated $H \rightarrow Z_DZ_D$ events with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in bins of $m_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the HAHM. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: 10-20. Hatched histograms show statistical uncertainties in the total expected background. Signal contributions expected from simulated $H \rightarrow Z_DZ_D$ events with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in bins of $m_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the HAHM. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: >20. Hatched histograms show statistical uncertainties in the total expected background. Signal contributions expected from simulated $H \rightarrow Z_DZ_D$ events with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 10\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 20\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 30\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 40\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 50\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 60\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 10\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 20\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 30\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 40\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 50\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 60\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 125\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The predicted cross section for $m(\tilde{q}) = 125\ GeV$ is 7200 pb, and falls outside the y-axis range.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 200\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The predicted cross section for $m(\tilde{q}) = 200 GeV$ is 840 pb, and falls outside the y-axis range.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 350\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The predicted cross section for $m(\tilde{q}) = 350\ GeV$ is 50 pb, and falls outside the y-axis range.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 700\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The gray horizontal line indicates the theoretical value of the squark-antisquark production cross section with the uncertainties shown as the gray shaded band.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 1150\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The gray horizontal line indicates the theoretical value of the squark-antisquark production cross section with the uncertainties shown as the gray shaded band.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 1600\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The gray horizontal line indicates the theoretical value of the squark-antisquark production cross section with the uncertainties shown as the gray shaded band.

Fractions of signal events with zero (green), one (blue), and two (red) STA muons matched to TMS muons by the STA to TMS association procedure, as a function of generated $L_{xy}$, in all HAHM signal samples combined.

Efficiencies of the Run 2 and Run 3 displaced dimuon triggers as a function of $c\tau$ for the HAHM signal events with $m(Z_D) = 50\ GeV$. The efficiency is defined as the fraction of simulated events that satisfy the requirements of the following sets of trigger paths: the Run 2 (2018) triggers (dashed black); the Run 3 (2022, L3) triggers (blue); the Run 3 (2022, L2) triggers (red); and the OR of all these triggers (Run 3 (2022), black). The lower panel shows the ratio of the overall Run 3 (2022) efficiency to the Run 2 (2018) efficiency.

Efficiencies of the Run 2 (2018) (red) and Run 3 (2022) (black) sets of displaced dimuon triggers as a function of $m(Z_D)$ for the HAHM signal events with $c\tau = 1\ cm$. The efficiency is defined as the fraction of simulated events that satisfy the detector acceptance and the requirements of the indicated set of trigger paths. The lower panel shows the ratio of the Run 3 (2022) efficiency to the Run 2 (2018) efficiency.

Efficiencies of the Run 2 (2018) (red) and Run 3 (2022) (black) sets of displaced dimuon triggers as a function of $m(Z_D)$ for the HAHM signal events with $c\tau = 10\ m$. The efficiency is defined as the fraction of simulated events that satisfy the detector acceptance and the requirements of the indicated set of trigger paths. The lower panel shows the ratio of the Run 3 (2022) efficiency to the Run 2 (2018) efficiency.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM signal with $m(Z_D) = 20\ GeV$ in different years of data taking. Efficiencies are computed as the ratios of the number of simulated signal events in which at least one dimuon candidate passes all 2016 (dashed green), 2018 (dashed red), and 2022 (solid black) trigger and offline selection criteria to the total number of simulated signal events. The lower panel shows the ratio of the 2022 efficiency to the 2018 efficiency (dashed red) and to the 2016 efficiency (dashed green).

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM signal with $m(Z_D) = 50\ GeV$ in different years of data taking. Efficiencies are computed as the ratios of the number of simulated signal events in which at least one dimuon candidate passes all 2016 (dashed green), 2018 (dashed red), and 2022 (solid black) trigger and offline selection criteria to the total number of simulated signal events. The lower panel shows the ratio of the 2022 efficiency to the 2018 efficiency (dashed red) and to the 2016 efficiency (dashed green).

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 10\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 20\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 30\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 40\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 60\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 125\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 200\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 350\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 700\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 1150\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 1600\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 700\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 500\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 1150\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 500\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 1600\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 500\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the TMS-TMS dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ smaller than 20 cm in the HAHM signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ smaller than 20 cm in the HAHM signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the TMS-TMS dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 20-70 cm in the HAHM signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 20-70 cm in the HAHM signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 70-500 cm in the HAHM signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the TMS-TMS dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ smaller than 20 cm in the RPV SUSY signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ smaller than 20 cm in the RPV SUSY signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the TMS-TMS dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 20-70 cm in the RPV SUSY signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 20-70 cm in the RPV SUSY signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 70-500 cm in the RPV SUSY signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Truth $P_{T}^{Z}$ distribution after reweghting to data.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Unfolded $M(l^{+}\gamma\gamma)$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons.

Unfolded $M(l^{+}\gamma\gamma)$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for leading photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with bare leptons for leading photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for second photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with bare leptons for second photon.

Unfolded dRgg distribution for $Z \to ll\gamma\gamma$ process with dressed leptons.

Unfolded dRgg distribution for $Z \to ll\gamma\gamma$ process with bare leptons.

Unfolded $P_{T}^{\gamma 1}$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for leading photon.

Unfolded $P_{T}^{\gamma 1}$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons for leading photon.

Unfolded $P_{T}^{\gamma 2}$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for second photon.

Unfolded $P_{T}^{\gamma 2}$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons for second photon.

Acceptance times efficiency for signal grid in merged two-prong signal region.

Acceptance times efficiency for signal grid in resolved two-prong signal region.

Acceptance times efficiency for signal grid in resolved two-prong signal region.

The obtained p-value from comparing data to background estimation across all bins, defined by windows in the X and Y particle masses, in the anomaly signal region.

The obtained p-value from comparing data to background estimation across all bins, defined by windows in the X and Y particle masses, in the anomaly signal region.

The expected 95% CL limits on the cross-section $\sigma(pp \rightarrow Y \rightarrow XH \rightarrow q\bar{q}b\bar{b}$) in pb in the two-dimensional space of $m_Y$ versus $m_X$, obtained from a simultaneous fit of both merged and resolved two-prong signal regions with all statistical and systematic uncertainties.

The expected 95% CL limits on the cross-section $\sigma(pp \rightarrow Y \rightarrow XH \rightarrow q\bar{q}b\bar{b}$) in pb in the two-dimensional space of $m_Y$ versus $m_X$, obtained from a simultaneous fit of both merged and resolved two-prong signal regions with all statistical and systematic uncertainties.

The observed 95% CL limits on the cross-section $\sigma(pp \rightarrow Y \rightarrow XH \rightarrow q\bar{q}b\bar{b}$) in pb in the two-dimensional space of $m_Y$ versus $m_X$, obtained from a simultaneous fit of both merged and resolved two-prong signal regions with all statistical and systematic uncertainties.

The observed 95% CL limits on the cross-section $\sigma(pp \rightarrow Y \rightarrow XH \rightarrow q\bar{q}b\bar{b}$) in pb in the two-dimensional space of $m_Y$ versus $m_X$, obtained from a simultaneous fit of both merged and resolved two-prong signal regions with all statistical and systematic uncertainties.