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.

Upper limits on $\sigma\times B$ as a function of $m_{e^*}$ in the $e\nu J$ channel. The $\pm 1(2)\sigma$ uncertainty bands around the expected limit represent all sources of systematic and statistical uncertainties.

Lower limits on $\Lambda$ as a function of $m_{e^*}$ for the $eejj$ channel. The $\pm\ 1(2)\sigma$ uncertainty bands around the expected limit represent all sources of systematic and statistical uncertainties. The limits for $m_{e^*}\gt4$ TeV are the result of extrapolation.

Lower limits on $\Lambda$ as a function of $m_{e^*}$ for the $e\nu J$ channel. The $\pm\ 1(2)\sigma$ uncertainty bands around the expected limit represent all sources of systematic and statistical uncertainties.

Lower limits on $\Lambda$ as a function of $m_{e^*}$ combined for the $eejj$ and the $e\nu J$ channels. The $\pm\ 1(2)\sigma$ uncertainty bands around the expected limit represent all sources of systematic and statistical uncertainties. The limits for $m_{e^*}\gt4$ TeV are the result of extrapolation.

Ratio of unfolded width of azimuthal angular correlation distributions (P PB/ P P). Different colors correspond to different combinations of p_{T,1} and p_{T,2}

Ratio of unfolded Dijet conditional yields (P PB/ P P). Different colors correspond to different combinations of p_{T,1} and p_{T,2}

Unfolded width of azimuthal angular correlation distributions (Delta p_{T} > 3). Full markers represent p+Pb, open markers p+p

Unfolded Dijet conditional yields (Delta p_{T} > 3). Full markers represent p+Pb, open markers p+p

Ratio of unfolded width of azimuthal angular correlation distributions (P PB/ P P) (Delta p_{T} > 3). Different colors correspond to different combinations of p_{T,1} and p_{T,2}

Ratio of unfolded Dijet conditional yields (P PB/ P P) (Delta p_{T} > 3). Different colors correspond to different combinations of p_{T,1} and p_{T,2}

Unfolded azimuthal angular correlation distributions. Black markers represent p+Pb, red markers p+p

Unfolded azimuthal angular correlation distributions (Delta p_{T} > 3). Black markers represent p+Pb, red markers p+p

The symmetric cumulant $sc_{2\,3}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The symmetric cumulant $sc_{2\,3}\{4\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The symmetric cumulant $sc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The symmetric cumulant $sc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The symmetric cumulant $sc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The symmetric cumulant $sc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The symmetric cumulant $sc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The symmetric cumulant $sc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The symmetric cumulant $sc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The asymmetric cumulant $ac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The asymmetric cumulant $ac_{2}\{3\}$results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The asymmetric cumulant $ac_{2}\{3\}$results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The asymmetric cumulant $ac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The asymmetric cumulant $ac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The asymmetric cumulant $ac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The normalized symmetric cumulant $nsc_{2\,3}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The normalized symmetric cumulant $nsc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The normalized asymmetric cumulant $nac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The normalized symmetric cumulant $nsc_{2\,3}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The normalized symmetric cumulant $nsc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The normalized asymmetric cumulant $nac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 13 TeV

The normalized symmetric cumulant $nsc_{2\,3}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The normalized symmetric cumulant $nsc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The normalized asymmetric cumulant $nac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The normalized symmetric cumulant $nsc_{2\,3}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The normalized symmetric cumulant $nsc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The normalized asymmetric cumulant $nac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The normalized symmetric cumulant $nsc_{2\,3}\{4\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The normalized symmetric cumulant $nsc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The normalized asymmetric cumulant $nac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The normalized symmetric cumulant $nsc_{2\,3}\{4\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The normalized symmetric cumulant $nsc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The normalized asymmetric cumulant $nac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The $v_{2}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The $v_{3}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The $v_{4}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The $v_{2}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The $v_{3}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The $v_{4}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The $v_{2}\{2\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The $v_{3}\{2\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The $v_{4}\{2\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The $v_{2}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The $v_{3}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The $v_{4}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The $v_{2}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The $v_{3}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The $v_{4}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The $v_{2}\{2\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The $v_{3}\{2\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The $v_{4}\{2\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The symmetric cumulant $ac_{2}\{3\}$ in Pb+Pb from different methods

The symmetric cumulant $ac_{2}\{3\}$ in Pb+Pb from different methods

The symmetric cumulant $ac_{2}\{3\}$ in p+Pb from different methods

The symmetric cumulant $ac_{2}\{3\}$ in p+Pb from different methods

The symmetric cumulant $ac_{2}\{3\}$ in pp from different methods

The symmetric cumulant $ac_{2}\{3\}$ in pp from different methods

Differential cross-sections for EW $Zjj$ production as a function of $\Delta\phi_{jj}$ with breakdown of associated uncertainties. The statistical uncertainty is correlated across bins according to the statistical cross correlation matrix presented in Table 21.

Differential cross-sections for inclusive $Zjj$ production in the EW $Zjj$ signal region ($N^\mathrm{gap}_\mathrm{jets}=0$, $\xi_Z < 0.5$) as a function of $m_{jj}$ with breakdown of associated uncertainties.

Differential cross-sections for inclusive $Zjj$ production in the EW $Zjj$ signal region ($N^\mathrm{gap}_\mathrm{jets}=0$, $\xi_Z < 0.5$) as a function of $\Delta y_{jj}$ with breakdown of associated uncertainties.

Differential cross-sections for inclusive $Zjj$ production in the EW $Zjj$ signal region ($N^\mathrm{gap}_\mathrm{jets}=0$, $\xi_Z < 0.5$) as a function of $p_{\mathrm{T},\ell\ell}$ with breakdown of associated uncertainties.

Differential cross-sections for inclusive $Zjj$ production in the EW $Zjj$ signal region ($N^\mathrm{gap}_\mathrm{jets}=0$, $\xi_Z < 0.5$) as a function of $\Delta\phi_{jj}$ with breakdown of associated uncertainties.

Differential cross-sections for inclusive $Zjj$ production in control region CRa ($N^\mathrm{gap}_\mathrm{jets} \geq 1$, $\xi_Z < 0.5$) as a function of $m_{jj}$ with breakdown of associated uncertainties.

Differential cross-sections for inclusive $Zjj$ production in control region CRa ($N^\mathrm{gap}_\mathrm{jets} \geq 1$, $\xi_Z < 0.5$) as a function of $\Delta y_{jj}$ with breakdown of associated uncertainties.

Differential cross-sections for inclusive $Zjj$ production in control region CRa ($N^\mathrm{gap}_\mathrm{jets} \geq 1$, $\xi_Z < 0.5$) as a function of $p_{\mathrm{T},\ell\ell}$ with breakdown of associated uncertainties.

Differential cross-sections for inclusive $Zjj$ production in control region CRa ($N^\mathrm{gap}_\mathrm{jets} \geq 1$, $\xi_Z < 0.5$) as a function of $\Delta\phi_{jj}$ with breakdown of associated uncertainties.

Differential cross-sections for inclusive $Zjj$ production in control region CRb ($N^\mathrm{gap}_\mathrm{jets} \geq 1$, $\xi_Z > 0.5$) as a function of $m_{jj}$ with breakdown of associated uncertainties.

Differential cross-sections for inclusive $Zjj$ production in control region CRb ($N^\mathrm{gap}_\mathrm{jets} \geq 1$, $\xi_Z > 0.5$) as a function of $\Delta y_{jj}$ with breakdown of associated uncertainties.

Differential cross-sections for inclusive $Zjj$ production in control region CRb ($N^\mathrm{gap}_\mathrm{jets} \geq 1$, $\xi_Z > 0.5$) as a function of $p_{\mathrm{T},\ell\ell}$ with breakdown of associated uncertainties.

Differential cross-sections for inclusive $Zjj$ production in control region CRb ($N^\mathrm{gap}_\mathrm{jets} \geq 1$, $\xi_Z > 0.5$) as a function of $\Delta\phi_{jj}$ with breakdown of associated uncertainties.

Differential cross-sections for inclusive $Zjj$ production in control region CRc ($N^\mathrm{gap}_\mathrm{jets} = 0$, $\xi_Z > 0.5$) as a function of $m_{jj}$ with breakdown of associated uncertainties.

Differential cross-sections for inclusive $Zjj$ production in control region CRc ($N^\mathrm{gap}_\mathrm{jets} = 0$, $\xi_Z > 0.5$) as a function of $\Delta y_{jj}$ with breakdown of associated uncertainties.

Differential cross-sections for inclusive $Zjj$ production in control region CRc ($N^\mathrm{gap}_\mathrm{jets} = 0$, $\xi_Z > 0.5$) as a function of $p_{\mathrm{T},\ell\ell}$ with breakdown of associated uncertainties.

Differential cross-sections for inclusive $Zjj$ production in control region CRc ($N^\mathrm{gap}_\mathrm{jets} = 0$, $\xi_Z > 0.5$) as a function of $\Delta\phi_{jj}$ with breakdown of associated uncertainties.

Statistical correlation between bins of the differential cross-section measurements of EW $Zjj$ production in the signal region ($N^\mathrm{gap}_\mathrm{jets} = 0$, $\xi_Z < 0.5$). The associated measured central values, statistical uncertainty magnitude and bin ranges are presented in Tables 1-4. The bins are presented in order such that for example mjj_bin_1 spans $m_{jj} \in (1000,1500)$ GeV (see Table 1), while pT_ll_bin7 spans $p_\mathrm{T}^{ll} \in (200,275)$ GeV (see Table 3).

Differential fiducial cross section of the Z boson $|y|$ in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the leading b-jet $|y|$ in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta \phi$ between Z boson and leading $b$-jet in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta y$ between Z boson and leading $b$-jet in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta R$ between Z boson and leading $b$-jet in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta \phi$ between the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta y$ between the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta R$ between the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the invariant mass of the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the Z boson $p_{\text{T}}$ in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $p_{\text{T}}$ of the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the ratio between the $p_{\text{T}}$ and the invariant mass of the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Transverse profile for 70 GeV < pT < 100 GeV.

Longitudinal profile for pT > 100 GeV.

Transverse profile for pT > 100 GeV.

Average longitudinal profile

Transverse profile for pT > 100 GeV.

Expected and observed 95% CL upper limits on $\sigma(pp\to H)B(H\to Z(\eta_c~\text{or}~J/\psi))/$pb. The smaller (larger) quoted ranges around the expected limits represent $\pm 1\sigma$ ($\pm 2\sigma$) fluctuations.

Post-fit BDT distributions after the VBF event selection for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. The VBF signal is shown for $\mu = 0.73$ and $\tilde d = -0.01$. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given.

The observed $\Delta\mathrm{NLL}$ curve as a function of $\tilde d$ values. For comparison, expected $\Delta\mathrm{NLL}$ curves are also shown. The constraints on the nuisance parameters and normalization factors are first determined in a CR-only fit, and then a fit to pseudo-data corresponding to these nuisance parameters, normalization factors, and to $\tilde d=0, \mu = 1$ or $\tilde d =0, \mu = 0.73$ is performed to obtain these $\Delta\mathrm{NLL}$ curves. Moreover, a pre-fit expected $\Delta\mathrm{NLL}$ is shown, using pseudo-data corresponding to $\tilde d =0$ and $\mu = 1$ in the signal and control regions.

The expected $\Delta\mathrm{NLL}$ curves comparing the sensitivity of the fit with and without systematic uncertainties. For comparison, other curves are shown which remove the effect of jet-based systematic uncertainties, $\tau$-based systematic uncertainties, and MC statistical uncertainties.

The observed $\Delta\mathrm{NLL}$ curves for each analysis channel as a function of $\tilde d$, compared to the combined result. For the individual analysis channel $\Delta\mathrm{NLL}$ curves, only event yield information in the other SRs is used, ensuring that the Optimal Observable distributions in the other SRs do not influence the preferred value of $\tilde d$. The signal strength is constrained to be positive in these individual channel $\Delta\mathrm{NLL}$ curves. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit BDT distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit BDT distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit BDT distributions in the $Z\to \ell\ell$ CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit Optimal Observable distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit Optimal Observable distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit Optimal Observable distributions in the $Z\to \ell\ell$ CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit distribution of weighted event yields as a function of the Optimal Observable for all four SRs combined. The contributions of the different SRs are weighted by a factor of ln(1 + S/B), where S and B are the post-fit expected numbers of signal and background events in that region, respectively. The size of the combined statistical, experimental, and theoretical uncertainties is given.