Distribution of the BDT output score in the the 3l channel signal region.

Distribution of the BDT output score in the 2lSC channel signal region.

Distribution of the BDT output score in the 2lSC+$\tau_{had}$ channel signal region.

Distribution of the BDT output score in the the 2l+2$\tau_{had}$ channel signal region.

Distribution of the BDT output score in the l+2$\tau_{had}$ channel signal region.

Invariant mass of the diphoton system in the $\gamma\gamma$+2(l,$\tau_{had}$) channel signal region.

Invariant mass of the diphoton system in the $\gamma\gamma$+$e/\mu$ channel in the loose BDT signal region.

Invariant mass of the diphoton system in the $\gamma\gamma$+$e/\mu$ channel in the medium BDT signal region.

Invariant mass of the diphoton system in the $\gamma\gamma$+$e/\mu$ channel in the tight BDT signal region.

Invariant mass of the diphoton system in the $\gamma\gamma$+$\tau$ channel in the loose BDT signal region.

Invariant mass of the diphoton system in the $\gamma\gamma$+$\tau$ channel in the medium BDT signal region.

Invariant mass of the diphoton system in the $\gamma\gamma$+$\tau$ channel in the tight BDT signal region.

Observed and expected 95% CL upper limits on the signal strength for the ML channels and the $\gamma\gamma$+ML channels.

Expected values of -2ln$\Lambda$ as a function of $\kappa_\lambda$ for the combination. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Observed values of -2ln$\Lambda$ as a function of $\kappa_\lambda$ for the combination. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Expected values of -2ln$\Lambda$ as a function of $\kappa_{2V}$ for the combination. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Observed values of -2ln$\Lambda$ as a function of $\kappa_{2V}$ for the combination. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Expected values of -2ln$\Lambda$ as a function of $\kappa_\lambda$ for the ML channels. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Observed values of -2ln$\Lambda$ as a function of $\kappa_\lambda$ for the ML channels. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Expected values of -2ln$\Lambda$ as a function of $\kappa_\lambda$ for the $\gamma\gamma$+ML channels. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Observed values of -2ln$\Lambda$ as a function of $\kappa_\lambda$ for the $\gamma\gamma$+ML channels. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Expected values of -2ln$\Lambda$ as a function of $\kappa_{2V}$ for the ML channels. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Observed values of -2ln$\Lambda$ as a function of $\kappa_{2V}$ for the ML channels. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Expected values of -2ln$\Lambda$ as a function of $\kappa_{2V}$ for the $\gamma\gamma$+ML channels. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Observed values of -2ln$\Lambda$ as a function of $\kappa_{2V}$ for the $\gamma\gamma$+ML channels. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

The expected upper limits on the production cross-section times the product of branching ratios for the benchmark signal scenario involving a scalar particle $X$ with narrow width decaying via $X\rightarrow aa\rightarrow 4\gamma$, $\sigma_X\times B(X\rightarrow aa)\times B(a\rightarrow\gamma\gamma)^2$. The limits for $m_{a}$ = 5 GeV and 10 GeV do not cover as large a range as the other mass points, since the region of interest is limited to $ m_{a} < 0.01 \times m_{X}$. Additionally, the expected limits are not provided for a small number of points, indicated with a hyphen, because of a technical failure with the computation.

The observed upper limits on the production cross-section times the product of branching ratios for the benchmark signal scenario involving a scalar particle $X$ with narrow width decaying via $X\rightarrow aa\rightarrow 6\pi^0$, $\sigma_X\times B(X\rightarrow aa)\times B(a\rightarrow 3\pi^0)^2$. The limits for $m_{a}$ = 5 GeV and 10 GeV do not cover as large a range as the other mass points, since the region of interest is limited to $ m_{a} < 0.01 \times m_{X}$.

The expected upper limits on the production cross-section times the product of branching ratios for the benchmark signal scenario involving a scalar particle $X$ with narrow width decaying via $X\rightarrow aa\rightarrow 6\pi^0$, $\sigma_X\times B(X\rightarrow aa)\times B(a\rightarrow 3\pi^0)^2$. The limits for $m_{a}$ = 5 GeV and 10 GeV do not cover as large a range as the other mass points, since the region of interest is limited to $ m_{a} < 0.01 \times m_{X}$. Additionally, the expected limits are not provided for a small number of points, indicated with a hyphen, because of a technical failure with the computation.

Observed 95% CL upper limits on the visible cross section as a function of $m_X$ and the fraction of events in the low-$\Delta E$ category.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 0.1 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 0.5 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 0.7 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 1 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 2 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 5 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 10 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 0.5 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 0.7 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 1 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 2 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 5 GeV.

Selection efficiency for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 10 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 0.1 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 0.5 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 0.7 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 1 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 2 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 5 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 2\gamma$ with $m_a$ = 10 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 0.5 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 0.7 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 1 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 2 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 5 GeV.

Fraction of reconstructed photons with a value of shower shape variable $\Delta E$ lower than the threshold, for reconstructed photons originating from the decay $a\rightarrow 3\pi^0\rightarrow 6\gamma$ with $m_a$ = 10 GeV.

Selection efficiency for photons originating from the BSM process $X\rightarrow\gamma\gamma$, where the $X$ particle is a high-mass narrow-width scalar particle originating from the gluon--gluon fusion process.

Fraction of photons with a value of shower shape variable $\Delta E$ lower than the threshold, for photons originating from the BSM process $X\rightarrow\gamma\gamma$, where the $X$ particle is a high-mass narrow-width scalar particle originating from the gluon--gluon fusion process.