Suppl. Fig. 7. The reconstructed shower angle for all single photon events. The individual signal and background event type categories added together form the unconstrained prediction.

Suppl. Fig. 7. The constrained covariance matrix for the reconstructed shower angle for all single-photon events. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Suppl. Fig. 7. The unconstrained covariance matrix for reconstructed shower angle for all single-photon events. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Fig. 3. The reconstructed number of protons, with a proton KE energy threshold of 35 MeV, for all single-photon events. The individual signal and background event type categories added together form the unconstrained prediction.

Fig. 4, Suppl. Fig. 8. The reconstructed shower energy for zero proton events binned from 0 to 600 MeV. The individual signal and background event type categories added together form the unconstrained prediction.

Fig. 4. The constrained covariance matrix for the reconstructed shower energy for zero proton events binned from 0 to 600 MeV. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Fig. 4. The unconstrained covariance matrix for reconstructed shower energy for zero proton events binned from 0 to 600 MeV. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Fig. 6. The reconstructed shower energy for zero proton events binned from 150 to 1250 MeV, with LEE prediction. The individual signal and background event type categories added together form the unconstrained prediction. The "LEE" category is the unfolded MiniBooNE LEE model prediction scaled under the neutrino-target nucleon interaction assumption.

Fig. 6. The constrained covariance matrix for the reconstructed shower energy for zero proton events binned from 150 to 1250 MeV. Both with and without the LEE model added. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Fig. 6. The unconstrained covariance matrix for reconstructed shower energy for zero proton events binned from 150 to 1250 MeV. Both with and without the LEE model added. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Suppl. Fig. 10. The reconstructed shower backwards projected distance for zero proton events. The individual signal and background event type categories added together form the unconstrained prediction.

Suppl. Fig. 9. The reconstructed shower angle for zero proton events. Includes LEE model prediction. The individual signal and background event type categories added together form the unconstrained prediction. The "LEE" category is the unfolded MiniBooNE LEE model prediction scaled under the neutrino-target nucleon interaction assumption.

Suppl. Fig. 9. The constrained covariance matrix for the reconstructed shower angle for zero proton events. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Suppl. Fig. 9. The unconstrained covariance matrix for reconstructed shower angle for zero proton events. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

1$\gamma$Xp shower energy and angle selection efficiencies for true 1$\gamma$Xp events. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. Energy and Angle values represent bin lower edges.

1$\gamma$0p shower energy and angle selection efficiencies for true single-photon events with no other true particles of any energy. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. Energy and Angle values represent bin lower edges.

Covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Pearson data statistical uncertainties have been included, and include small correlations due to events which can be selected by both WC and Pandora. The first four bins are the WC $1\gamma Np$, WC $1\gamma 0p$, Pandora $1\gamma 1p$, and Pandora $1\gamma 0p$ channels, and the next four sets of 16 bins are the NC $\pi^0 Np$, NC $\pi^0 0p$, $\nu_\mu$CC $Np$, and $\nu_\mu$CC $0p$ channels. This corresponds to Fig. 1 of the paper and Fig. 6 of the Supplemental Material.

Constrained WC $1\gamma Np$ energy distribution. There are 6 bins from 0 to 600 MeV of reconstructed shower energy, with an additional overflow bin. Constrained systematic uncertainties on the predictions are illustrated, and a more detailed covariance matrix is included in the Constrained Energy Distributions Covariance Matrix tab. This corresponds to Fig. 2 of the paper.

Constrained WC $1\gamma 0p$ energy distribution. There are 14 bins from 0 to 700 MeV of reconstructed shower energy, with an additional overflow bin. Constrained systematic uncertainties on the predictions are illustrated, and a more detailed covariance matrix is included in the Constrained Energy Distributions Covariance Matrix tab. This corresponds to Fig. 2 of the paper.

Constrained Pandora $1\gamma 1p$ energy distribution. There are 6 bins from 0 to 600 MeV of reconstructed shower energy, with an additional overflow bin. Constrained systematic uncertainties on the predictions are illustrated, and a more detailed covariance matrix is included in the Constrained Energy Distributions Covariance Matrix tab. This corresponds to Fig. 2 of the paper.

Constrained Pandora $1\gamma 0p$ energy distribution. There are 14 bins from 0 to 700 MeV of reconstructed shower energy, with an additional overflow bin. Constrained systematic uncertainties on the predictions are illustrated, and a more detailed covariance matrix is included in the Constrained Energy Distributions Covariance Matrix tab. This corresponds to Fig. 2 of the paper.

Constrained covariance matrix showing constrained uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Pearson data statistical uncertainties have been included, and include small correlations due to events which can be selected by both WC and Pandora. The first seven bins are the WC $1\gamma Np$ reconstructed shower energy from 0-600 MeV with an overflow bin, the next 15 bins are the WC $1\gamma 0p$ reconstructed shower energy from 0-700 MeV with an overflow bin, the next seven bins are the Pandora $1\gamma 1p$ reconstructed shower energy from 0-600 MeV with an overflow bin, and the last 15 bins are the Pandora $1\gamma 0p$ reconstructed shower energy from 0-700 MeV with an overflow bin. This corresponds to Fig. 2 of the paper.

2D efficiency matrix showing the efficiency of the WC $1\gamma Np$ selection with respect to all true NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 15 bins from 0 to 1500 MeV in true photon energy. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. This corresponds to Fig. 9 of the Supplemental Material.

2D efficiency matrix showing the efficiency of the WC $1\gamma 0p$ selection with respect to all true NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 30 bins from 0 to 1500 MeV in true photon energy. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. This corresponds to Fig. 9 of the Supplemental Material.

2D efficiency matrix showing the efficiency of the Pandora $1\gamma 1p$ selection with respect to all true NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 15 bins from 0 to 1500 MeV in true photon energy, with lower bin edges indicated in the table. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. This corresponds to Fig. 9 of the Supplemental Material.

2D efficiency matrix showing the efficiency of the Pandora $1\gamma 0p$ selection with respect to all true NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 30 bins from 0 to 1500 MeV in true photon energy, with lower bin edges indicated in the table. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. This corresponds to Fig. 9 of the Supplemental Material.

2D efficiency matrix showing the efficiency of the WC $1\gamma 0p$ selection with respect to all true zero-primary-final-state-proton NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 30 bins from 0 to 1500 MeV in true photon energy, with lower bin edges indicated in the table. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied.

2D efficiency matrix showing the efficiency of the Pandora $1\gamma 0p$ selection with respect to all true zero-primary-final-state-proton NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 30 bins from 0 to 1500 MeV in true photon energy, with lower bin edges indicated in the table. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied.

Observed and expected 95 % CL exclusion limits on $\tan\beta$ as a function of $m_{H^{\pm}}$, shown in the context of the mh125 scenario, for $m_{H^{\pm}}>150$ GeV and $(1 \leq \tan\beta \leq 60)$. The surrounding shaded bands correspond to the 1$\sigma$ and 2$\sigma$ confidence intervals around the expected limit.

Observed and expected 95 % CL exclusions on $\tan\beta$ as a function of $m_{H^{\pm}}$ derived from exclusion limits on $\mathrm{\cal{B}}(t\to bH^+)\times \mathrm{\cal{B}}(H^+ \to \tau \nu)$, for Type I 2HDM low mass scenario. The surrounding shaded bands correspond to the 1$\sigma$ and 2$\sigma$ confidence intervals around the expected limit.

Observed and expected 95 % CL exclusions on $\tan\beta$ as a function of $m_{H^{\pm}}$ derived from exclusion limits on $\mathrm{\cal{B}}(t\to bH^+)\times \mathrm{\cal{B}}(H^+ \to \tau \nu)$, for Lepton Specific 2HDM low mass scenario. The surrounding shaded bands correspond to the 1$\sigma$ and 2$\sigma$ confidence intervals around the expected limit.

The measured differential cross-sections and the predictions from the GM-VFNS calculation for the $D_s^\pm$ meson in bins of $p_T$ for $|\eta| < 2.5$. The statistical, systematic (excluding branching ratio) and branching ratio uncertainties are shown separately for data, while the total theory uncertainties are shown for GM-VFNS.

Inclusive $D^\pm$ meson production cross-sections in different fiducial volumes defined by $|\eta|<2.5$ and different $p_T$ regions. For the ATLAS measurements, the statistical, systematic (excluding branching ratio) and branching ratio uncertainties are shown separately; for the theoretical predictions, the total theoretical uncertainties are shown.

Inclusive $D_s^\pm$ meson production cross-sections in different fiducial volumes defined by $|\eta|<2.5$ and different $p_T$ regions. For the ATLAS measurements, the statistical, systematic (excluding branching ratio) and branching ratio uncertainties are shown separately; for the theoretical predictions, the total theoretical uncertainties are shown.

The distribution of the $p_\mathrm{T}{}_{\mu\mathrm{-had}}^\mathrm{col}$ in the VR. `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top-quark, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic sources.

The distributions of $m_{\mu\mathrm{-had}}^\mathrm{col}$ corresponding to the SR selection criteria defined in the text. `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top-quark, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic sources.

The distributions of $m_{\mu\mathrm{-had}}^\mathrm{col}$ corresponding to the SR_\text{std} selection criteria defined in the text. `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top-quark, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic sources.

The distributions of $m_{\mu\mathrm{-had}}^\mathrm{col}$ corresponding to the SR_\text{tight} selection criteria defined in the text. `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top-quark, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic sources.

The distributions of $m_{\mu\mathrm{-had}}^\mathrm{col}$ corresponding to the SR_\text{tight}^\text{std} selection criteria defined in the text. `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top-quark, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic sources.

The distribution of $\Delta\phi^\mathrm{signed}_{\mu-MET}$ with all other event selection criteria applied in the SR; the straight lines indicate the positions of the selection requirements. `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic errors.

The distribution of $\Delta\phi^\mathrm{signed}_{\mu-MET}$ with all other event selection criteria applied in the VR; the straight lines indicate the positions of the selection requirements. `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic errors.

The distribution of the $m_{\mu\mathrm{-had}}^\mathrm{col}$ in the SR. `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic errors.

The distribution of the $m_{\mu\mathrm{-had}}^\mathrm{col}$ in the VR. `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic errors.

The distribution of the $p_\mathrm{T}$ of the $\tau_\mathrm{had}^{\mu\mkern-10mu\backslash}$ in the SR. `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic errors.

The distribution of the $p_\mathrm{T}$ of the $\tau_\mathrm{had}^{\mu\mkern-10mu\backslash}$ in the VR. `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic errors.

The distribution of the $p_\mathrm{T}$ of the muon removed from the $\tau_\mathrm{had}^{\mu\mkern-10mu\backslash}$ in the SR `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic errors.

The distribution of the $p_\mathrm{T}$ of the muon removed from the $\tau_\mathrm{had}^{\mu\mkern-10mu\backslash}$ in the VR `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic errors.

The distribution of the $p_\mathrm{T}$ of the leading jet in the SR. `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic errors.

The distribution of the $p_\mathrm{T}$ of the leading jet in the VR. `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic errors.

The distribution of the $\Delta{R}_{\tau\tau}^\mathrm{reco}$ between the muon and the $\tau_\mathrm{had}^{\mu\mkern-10mu\backslash}$ in the SR `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic errors.

The distribution of the $\Delta{R}_{\tau\tau}^\mathrm{reco}$ between the muon and the $\tau_\mathrm{had}^{\mu\mkern-10mu\backslash}$ in the VR `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic errors.

The distribution of the reconstructed $E_\mathrm{T}^\mathrm{miss}$ in the SR `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic errors.

The distribution of the reconstructed $E_\mathrm{T}^\mathrm{miss}$ in the VR `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic errors.

The observed and expected ($\pm2\sigma$) upper limits at 95% CL on the branching ratio for $H\rightarrow aa\rightarrow \gamma\gamma\tau\tau$ as a function of the resonance mass hypothesis $m_{a}$.

The polynomial parameterisation of the total signal selection efficiency as a function of $m_{a}$, defined as the fraction of the reconstructed and selected events out of the total number of events in the signal simulated samples.

Fig. 3 bottom figure - Distributions of MC simulation compared with data for reconstructed shower energy in the 1$e$0$p$0$\pi$ signal channel, along with the LEE Signal Model 2. The signal and background event categories are summed to form the unconstrained prediction (excluding LEE). Signal events correspond to $\nu_e$ CC events. Background events include $\nu$ with $\pi^0$ events, $\nu$ other events, and cosmic ray events. In Fig. 3, the LEE component is plotted on top of the constrained prediction (excluding LEE) for illustrative purposes. In all statistical tests (results summarized in Table I), the prediction under an LEE hypothesis corresponds to a constrained prediction including LEE. The statistical uncertainties of data use a combined Neyman-Pearson (CNP) version (Eq.(19) in https://doi.org/10.1016/j.nima.2020.163677).

Fig. 4 top figure - Distributions of MC simulation compared with data for reconstructed shower cos($\theta$) in the 1$e$N$p$0$\pi$ signal channel, along with the LEE Signal Model 2. The signal and background event categories are summed to form the unconstrained prediction (excluding LEE). Signal events correspond to $\nu_e$ CC events. Background events include $\nu$ with $\pi^0$ events, $\nu$ other events, and cosmic ray events. In Fig. 4, the LEE component is plotted on top of the constrained prediction (excluding LEE) for illustrative purposes. In all statistical tests (results summarized in Table I), the prediction under an LEE hypothesis corresponds to a constrained prediction including LEE. The statistical uncertainties of data use a combined Neyman-Pearson (CNP) version (Eq.(19) in https://doi.org/10.1016/j.nima.2020.163677).

Fig. 4 bottom figure - Distributions of MC simulation compared with data for reconstructed shower cos($\theta$) in the 1$e$0$p$0$\pi$ signal channel, along with the LEE Signal Model 2. The signal and background event categories are summed to form the unconstrained prediction (excluding LEE). Signal events correspond to $\nu_e$ CC events. Background events include $\nu$ with $\pi^0$ events, $\nu$ other events, and cosmic ray events. In Fig. 4, the LEE component is plotted on top of the constrained prediction (excluding LEE) for illustrative purposes. In all statistical tests (results summarized in Table I), the prediction under an LEE hypothesis corresponds to a constrained prediction including LEE. The statistical uncertainties of data use a combined Neyman-Pearson (CNP) version (Eq.(19) in https://doi.org/10.1016/j.nima.2020.163677).

Fig. 1 top figure - Reconstructed neutrino energy distributions of the 1$\mu$N$p$0$\pi$ control sample (sideband). The signal and background event categories are summed to form the unconstrained prediction. Signal events correspond to $\nu_\mu$ CC events. Background events include $\nu_e$ CC events, NC without $\pi^0$ events, NC with $\pi^0$ events, and cosmic ray events. The statistical uncertainties of data use a combined Neyman-Pearson (CNP) version (Eq.(19) in https://doi.org/10.1016/j.nima.2020.163677).

Fig. 1 middle figure - Reconstructed neutrino energy distributions of the 1$\mu$0$p$0$\pi$ control sample (sideband). The signal and background event categories are summed to form the unconstrained prediction. Signal events correspond to $\nu_\mu$ CC events. Background events include $\nu_e$ CC events, NC without $\pi^0$ events, NC with $\pi^0$ events, and cosmic ray events. The statistical uncertainties of data use a combined Neyman-Pearson (CNP) version (Eq.(19) in https://doi.org/10.1016/j.nima.2020.163677).

Fig. 1 bottom figure - Reconstructed neutrino energy distributions of the NC$\pi^{0}$ control sample (sideband). The signal and background event categories are summed to form the unconstrained prediction. Signal events correspond to NC with $\pi^0$ events. Background events include $\nu_\mu$ CC events, $\nu_e$ CC events, NC without $\pi^0$ events, and cosmic ray events. Only bins up to 1 GeV in this sideband are released and are used in the constraint procedure due to low statistics above this energy. The statistical uncertainties of data use a combined Neyman-Pearson (CNP) version (Eq.(19) in https://doi.org/10.1016/j.nima.2020.163677).

The covariance matrix includes both signal channels and sideband channels, with the LEE component either excluded or included in the signal channels. All channels are binned in reconstructed neutrino energy. Fig. 2 in Supplemental Material shows the corresponding correlation matrix. The constraint procedure uses the correlations with sideband channels to update the prediction and its uncertainty in the signal region. The matrix accounts for uncertainties and correlations between bins due to flux, cross-section, hadron reinteraction, detector systematics, MC statistics, and cosmic ray statistics. Data statistical uncertainties are not included.

The covariance matrix includes both signal channels and sideband channels, with the LEE component either excluded or included in the signal channels. Signal channels are binned shower energy and sideband channels are binned in reconstructed neutrino energy. The constraint procedure uses the correlations with sideband channels to update the prediction and its uncertainty in the signal region. The matrix accounts for uncertainties and correlations between bins due to flux, cross-section, hadron reinteraction, detector systematics, MC statistics, and cosmic ray statistics. Data statistical uncertainties are not included.

The covariance matrix includes both signal channels and sideband channels, with the LEE component either excluded or included in the signal channels. Signal channels are binned shower cos($\theta$) and sideband channels are binned in reconstructed neutrino energy. The constraint procedure uses the correlations with sideband channels to update the prediction and its uncertainty in the signal region. The matrix accounts for uncertainties and correlations between bins due to flux, cross-section, hadron reinteraction, detector systematics, MC statistics, and cosmic ray statistics. Data statistical uncertainties are not included.

The covariance matrices of the sideband channels, binned in reconstructed neutrino energy. The systematic covariance matrix accounts for uncertainties and correlations between bins due to flux, cross-section, hadron reinteraction, detector systematics, MC statistics, and cosmic ray statistics. The total covariance matrix additionally includes data statistical uncertainties (CNP version). The $\chi^2$ values shown in Fig. 1 are calculated using this total covariance matrix, together with the data and MC predictions from the released sideband distribtions.

The unconstrained covariance matrices of the signal channels, binned in reconstructed neutrino energy, with the LEE component either excluded or included. The systematic covariance matrix accounts for uncertainties and correlations between bins due to flux, cross-section, hadron reinteraction, detector systematics, MC statistics, and cosmic ray statistics. The total covariance matrix additionally includes data statistical uncertainties (CNP version). The unconstrained systematic uncertainty in neutrino energy, shown in Supplemental Material Fig. 19, is derived from this systematic covariance matrix (excluding LEE).

The constrained covariance matrices of the signal channels, binned in reconstructed neutrino energy, with the LEE component either excluded or included. The constrained systematic uncertainty in neutrino energy, shown in Supplemental Material Fig. 19, is derived from the systematic covariance matrix (excluding LEE). The $\chi^2$ values shown in Fig. 2 and Table I are calculated using these total covariance matrices, together with the data and constrained predictions from the released signal distributions under the $H_0$ hypothesis (LEE excluded) and $H_1$ hypothesis (LEE included).

The unconstrained covariance matrices of the signal channels, binned in reconstructed shower energy, with the LEE component either excluded or included. The systematic covariance matrix accounts for uncertainties and correlations between bins due to flux, cross-section, hadron reinteraction, detector systematics, MC statistics, and cosmic ray statistics. The total covariance matrix additionally includes data statistical uncertainties (CNP version). The unconstrained systematic uncertainty in shower energy, shown in Supplemental Material Fig. 20, is derived from this systematic covariance matrix (excluding LEE).

The constrained covariance matrices of the signal channels, binned in reconstructed shower energy, with the LEE component either excluded or included. The constrained systematic uncertainty in shower energy, shown in Supplemental Material Fig. 20, is derived from the systematic covariance matrix (excluding LEE). The $\chi^2$ values shown in Fig. 3 and Table I are calculated using these total covariance matrices, together with the data and constrained predictions from the released signal distributions under the $H_0$ hypothesis (LEE excluded) and $H_1$ hypothesis (LEE included).

The unconstrained covariance matrices of the signal channels, binned in reconstructed shower cos($\theta$), with the LEE component either excluded or included. The systematic covariance matrix accounts for uncertainties and correlations between bins due to flux, cross-section, hadron reinteraction, detector systematics, MC statistics, and cosmic ray statistics. The total covariance matrix additionally includes data statistical uncertainties (CNP version). The unconstrained systematic uncertainty in shower cos($\theta$), shown in Supplemental Material Fig. 21, is derived from this systematic covariance matrix (excluding LEE).

The constrained covariance matrices of the signal channels, binned in reconstructed shower cos($\theta$), with the LEE component either excluded or included. The constrained systematic uncertainty in shower cos($\theta$), shown in Supplemental Material Fig. 21, is derived from the systematic covariance matrix (excluding LEE). The $\chi^2$ values shown in Fig. 4 and Table I are calculated using these total covariance matrices, together with the data and constrained predictions from the released signal distributions under the $H_0$ hypothesis (LEE excluded) and $H_1$ hypothesis (LEE included).

Breakdown of relative systematic uncertainties in the measured ${t\bar{t}}$ cross-section. The quoted uncertainties are obtained by repeating the fit with a group of nuisance parameters fixed to their fitted values and subtracting in quadrature the resulting total uncertainty from the uncertainty of the complete fit. However, the total uncertainty is not the quadratic sum of the grouped impacts, as this approach neglects the correlation among the different groups.

The figure shows the measurement of the top quark pair production cross-section, $\sigma_{t\bar{t}}$, in lead-lead (Pb+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV, based on 1.9 nb$^{-1}$ of data.

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the resolved qqbb signal region.

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the resolved lvbb signal region.

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the merged lvbb signal region

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the merged qqbb signal region

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 5jex3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 5jex4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 6jin3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 6jin4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 5jex3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 5jex4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 6jin3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 6jin4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 5jex3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 5jex4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 6jin3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 6jin4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

The mWh distributions in the Low-NN-score SRs of the merged qqbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Low-NN-score SRs of the merged qqbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged qqbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged qqbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Low-NN-score SRs of the merged lvbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Low-NN-score SRs of the merged lvbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Medium-NN-score SRs of the merged lvbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Medium-NN-score SRs of the merged lvbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged lvbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged lvbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

Cutflow table for a few selected signal samples after applying the selection requirements of the resolved analysis. At the initial cut stage, all events from the tbH ± → W ± h(→ b b) production process (including those from the fully hadronic decay channel) are considered.

Cutflow table for a few selected signal samples after applying the selection requirements of the merged analysis. At the initial cut stage, all events from the tbH ± → W ± h(→ bb) production process (including those from the fully hadronic decay channel) are considered.