Systematic covariance matrix for the 6 bins of the normalized differential cross section as a function of $\Delta|y_\mathrm{t}|$.

Values of the 12 bins of the normalized differential cross section as a function of $\Delta|\eta_{\ell}|$.

Statistical covariance matrix for the 12 bins of the normalized differential cross section as a function of $\Delta|\eta_{\ell}|$.

Systematic covariance matrix for the 12 bins of the normalized differential cross section as a function of $\Delta|\eta_{\ell}|$.

Values of $A_{\mathrm{C}}$ in the 3 bins of $M_\mathrm{t\bar{t}}$, and inclusively (bottom row). The value 9999 is used as a placeholder for infinity.

Statistical correlation matrix for $A_{\mathrm{C}}$ in the 3 bins of $M_\mathrm{t\bar{t}}$.

Systematic correlation matrix for $A_{\mathrm{C}}$ in the 3 bins of $M_\mathrm{t\bar{t}}$.

Fraction of events in each of the $6\times3$ bins of the normalized differential cross section as a function of $\Delta|y_\mathrm{t}|$ and $M_\mathrm{t\bar{t}}$. The value 9999 is used as a placeholder for infinity.

Statistical covariance matrix for the $6\times3$ bins of the normalized differential cross section as a function of $\Delta|y_\mathrm{t}|$ and $M_\mathrm{t\bar{t}}$.

Systematic covariance matrix for the $6\times3$ bins of the normalized differential cross section as a function of $\Delta|y_\mathrm{t}|$ and $M_\mathrm{t\bar{t}}$.

Values of $A^{\mathrm{lep}}_{\mathrm{C}}$ in the 3 bins of $M_\mathrm{t\bar{t}}$, and inclusively (bottom row). The value 9999 is used as a placeholder for infinity.

Statistical correlation matrix for $A^{\mathrm{lep}}_{\mathrm{C}}$ in the 3 bins of $M_\mathrm{t\bar{t}}$.

Systematic correlation matrix for $A^{\mathrm{lep}}_{\mathrm{C}}$ in the 3 bins of $M_\mathrm{t\bar{t}}$.

Fraction of events in each of the $12\times3$ bins of the normalized differential cross section as a function of $\Delta|\eta_{\ell}|$ and $M_\mathrm{t\bar{t}}$. The value 9999 is used as a placeholder for infinity.

Statistical covariance matrix for the $12\times3$ bins of the normalized differential cross section as a function of $\Delta|\eta_{\ell}|$ and $M_\mathrm{t\bar{t}}$.

Systematic covariance matrix for the $12\times3$ bins of the normalized differential cross section as a function of $\Delta|\eta_{\ell}|$ and $M_\mathrm{t\bar{t}}$.

Values of $A_{\mathrm{C}}$ in the 3 bins of $p_\mathrm{T}^\mathrm{t\bar{t}}$, and inclusively (bottom row). The value 9999 is used as a placeholder for infinity.

Statistical correlation matrix for $A_{\mathrm{C}}$ in the 3 bins of $p_\mathrm{T}^\mathrm{t\bar{t}}$.

Systematic correlation matrix for $A_{\mathrm{C}}$ in the 3 bins of $p_\mathrm{T}^\mathrm{t\bar{t}}$.

Fraction of events in each of the $6\times3$ bins of the normalized differential cross section as a function of $\Delta|y_\mathrm{t}|$ and $p_\mathrm{T}^\mathrm{t\bar{t}}$. The value 9999 is used as a placeholder for infinity.

Statistical covariance matrix for the $6\times3$ bins of the normalized differential cross section as a function of $\Delta|y_\mathrm{t}|$ and $p_\mathrm{T}^\mathrm{t\bar{t}}$.

Systematic covariance matrix for the $6\times3$ bins of the normalized differential cross section as a function of $\Delta|y_\mathrm{t}|$ and $p_\mathrm{T}^\mathrm{t\bar{t}}$.

Values of $A^{\mathrm{lep}}_{\mathrm{C}}$ in the 3 bins of $p_\mathrm{T}^\mathrm{t\bar{t}}$, and inclusively (bottom row). The value 9999 is used as a placeholder for infinity.

Statistical correlation matrix for $A^{\mathrm{lep}}_{\mathrm{C}}$ in the 3 bins of $p_\mathrm{T}^\mathrm{t\bar{t}}$.

Systematic correlation matrix for $A^{\mathrm{lep}}_{\mathrm{C}}$ in the 3 bins of $p_\mathrm{T}^\mathrm{t\bar{t}}$.

Fraction of events in each of the $12\times3$ bins of the normalized differential cross section as a function of $\Delta|\eta_{\ell}|$ and $p_\mathrm{T}^\mathrm{t\bar{t}}$. The value 9999 is used as a placeholder for infinity.

Statistical covariance matrix for the $12\times3$ bins of the normalized differential cross section as a function of $\Delta|\eta_{\ell}|$ and $p_\mathrm{T}^\mathrm{t\bar{t}}$.

Systematic covariance matrix for the $12\times3$ bins of the normalized differential cross section as a function of $\Delta|\eta_{\ell}|$ and $p_\mathrm{T}^\mathrm{t\bar{t}}$.

Values of $A_{\mathrm{C}}$ in the 3 bins of $|y_\mathrm{t\bar{t}}|$, and inclusively (bottom row). The value 9999 is used as a placeholder for infinity.

Statistical correlation matrix for $A_{\mathrm{C}}$ in the 3 bins of $|y_\mathrm{t\bar{t}}|$.

Systematic correlation matrix for $A_{\mathrm{C}}$ in the 3 bins of $|y_\mathrm{t\bar{t}}|$.

Fraction of events in each of the $6\times3$ bins of the normalized differential cross section as a function of $\Delta|y_\mathrm{t}|$ and $|y_\mathrm{t\bar{t}}|$. The value 9999 is used as a placeholder for infinity.

Statistical covariance matrix for the $6\times3$ bins of the normalized differential cross section as a function of $\Delta|y_\mathrm{t}|$ and $|y_\mathrm{t\bar{t}}|$.

Systematic covariance matrix for the $6\times3$ bins of the normalized differential cross section as a function of $\Delta|y_\mathrm{t}|$ and $|y_\mathrm{t\bar{t}}|$.

Values of $A^{\mathrm{lep}}_{\mathrm{C}}$ in the 3 bins of $|y_\mathrm{t\bar{t}}|$, and inclusively (bottom row). The value 9999 is used as a placeholder for infinity.

Statistical correlation matrix for $A^{\mathrm{lep}}_{\mathrm{C}}$ in the 3 bins of $|y_\mathrm{t\bar{t}}|$.

Systematic correlation matrix for $A^{\mathrm{lep}}_{\mathrm{C}}$ in the 3 bins of $|y_\mathrm{t\bar{t}}|$.

Fraction of events in each of the $12\times3$ bins of the normalized differential cross section as a function of $\Delta|\eta_{\ell}|$ and $|y_\mathrm{t\bar{t}}|$. The value 9999 is used as a placeholder for infinity.

Statistical covariance matrix for the $12\times3$ bins of the normalized differential cross section as a function of $\Delta|\eta_{\ell}|$ and $|y_\mathrm{t\bar{t}}|$.

Systematic covariance matrix for the $12\times3$ bins of the normalized differential cross section as a function of $\Delta|\eta_{\ell}|$ and $|y_\mathrm{t\bar{t}}|$.

$\Upsilon$(nS) polarization parameter $\tilde\lambda$ in the HX frame for $p_T$ of 10 to 15 GeV. The global uncertainties, independent of state and $N_{ch}$ bin, are also indicated.

$\Upsilon$(nS) polarization parameter $\lambda_\vartheta$ in the HX frame for $p_T$ of 15 to 35 GeV. The global uncertainties, independent of state and $N_{ch}$ bin, are also indicated.

$\Upsilon$(nS) polarization parameter $\lambda_\varphi$ in the HX frame for $p_T$ of 15 to 35 GeV. The global uncertainties, independent of state and $N_{ch}$ bin, are also indicated.

$\Upsilon$(nS) polarization parameter $\lambda_{\vartheta\varphi}$ in the HX frame for $p_T$ of 15 to 35 GeV. The global uncertainties, independent of state and $N_{ch}$ bin, are also indicated.

$\Upsilon$(nS) polarization parameter $\tilde\lambda$ in the HX frame for $p_T$ of 15 to 35 GeV. The global uncertainties, independent of state and $N_{ch}$ bin, are also indicated.

Observed event yields and estimated backgrounds in the low- and high-mass control regions. The uncertainty in the background yield is the sum in quadrature of the statistical and systematic uncertainties.

Summary of the relative systematic uncertainties in heavy Majorana neutrino signal yields and the background from prompt same-sign leptons, both estimated from simulation. The relative systematic uncertainties assigned to the data-driven backgrounds for the misidentified lepton background and mismeasured charge background are also shown. The uncertainties are given for the low-mass selections.

Summary of the relative systematic uncertainties in heavy Majorana neutrino signal yields and the background from prompt same-sign leptons, both estimated from simulation. The relative systematic uncertainties assigned to the data-driven backgrounds for the misidentified lepton background and mismeasured charge background are also shown. The uncertainties are given for the high-mass selections.

Summary of contributions to the systematic uncertainty related to the prompt same-sign leptons background, misidentified lepton background, and mismeasured charge background on the total background uncertainty for the case of $m_{\rm N}$ = 100 and 500 GeV.

Observed event yields and estimated backgrounds after the application of all selection, except for the final optimization. The background predictions from prompt same-sign leptons (Prompt bkgd.), misidentified leptons (Misid. bkgd.), mismeasured charge (Charge mismeas. bkgd.) and the total background (Total bkgd.) are shown together withthe number of events observed in data. The uncertainties shown are the statistical and systematic components, respectively.

Dielectron channel results after final optimization. The background predictions from prompt same-sign leptons, misidentified leptons, and mismeasured charge are shown along with the total background estimate and the number of events observed in data. The uncertainties shown are the statistical and systematic components, respectively.

Electron-muon channel results after final optimization. The background predictions from prompt same-sign leptons and misidentified leptons are shown along with the total background estimate and the number of events observed in data. The uncertainties shown are the statistical and systematic components, respectively.

Correlation matrix of systematic uncertainties for $\sigma^{\pm}$. The statistical and integrated luminosity uncertainties are not included. The first and the last eleven indices correspond to the systematic uncertainites of $W^+$ ($\sigma^+_{1-11}$) and $W^-$ ($\sigma^-_{1-11}$) cross sections, respectively, in 11 muon pseudorapidity bins. The values are expressed as percentages.

Correlation matrix of statistical uncertainties for $\sigma^{\pm}$. The first and the last eleven indices correspond to the statistical uncertainites of $W^+$ ($\sigma^+_{1-11}$) and $W^-$ ($\sigma^-_{1-11}$) cross sections, respectively, in 11 muon pseudorapidity bins. The values are expressed as percentages.

Correlation matrix of systematic uncertainties for $\mathcal{A}$. The values are expressed as percentages.

Observed and expected 95% CL upper limits on the product of cross section of a narrow resonance and the branching fraction $\sigma(gg \to X) B(X \to HH \to bbbb)$ for HPHP category. Theory curves corresponding to WED models with radion are superimposed.

Observed and expected 95% CL upper limits on the product of cross section of a narrow resonance and the branching fraction $\sigma(gg \to X) B(X \to HH \to bbbb)$ for HPLP category. Theory curves corresponding to WED models with radion are superimposed.

Observed and expected 95% CL upper limits on the product of cross section of a narrow resonance and the branching fraction $\sigma(gg \to X) B(X \to HH \to bbbb)$ for LPHP category. Theory curves corresponding to WED models with radion are superimposed.

Observed and expected 95% CL upper limits on the product of cross section and the branching fraction $\sigma(gg \to X) B(X \to HH \to bbbb)$ for HPHP, HPLP and LPHP categories combined. The one standard deviation on the 95% CL upper limit is also provided.

Cut flow of the baseline selection for the benchmark models $\tilde{g}\tilde{g}$, $\tilde{g}\rightarrow\text{b}\bar{\text{b}}\tilde{\chi}_{1}^{0}$ ($m_{\tilde{g}} = 1500$ GeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV), $\tilde{g}\tilde{g}$, $\tilde{g}\rightarrow\text{t}\bar{\text{t}}\tilde{\chi}_{1}^{0}$ ($m_{\tilde{g}} = 1500$ GeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV), and $\tilde{g}\tilde{g}$, $\tilde{g}\rightarrow\text{q}\bar{\text{q}}\tilde{\chi}_{1}^{0}$ ($m_{\tilde{g}} = 1000$ GeV, $m_{\tilde{\chi}_{1}^{0}} = 800$ GeV).

Expected numbers of signal events in search bins with $4\geq N_{\text{jets}}\geq6$ region for the benchmark models $\tilde{g}\tilde{g}$, $\tilde{g}\rightarrow\text{b}\bar{\text{b}}\tilde{\chi}_{1}^{0}$ ($m_{\tilde{g}} = 1500$ GeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV), $\tilde{g}\tilde{g}$, $\tilde{g}\rightarrow\text{t}\bar{\text{t}}\tilde{\chi}_{1}^{0}$ ($m_{\tilde{g}} = 1500$ GeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV), and $\tilde{g}\tilde{g}$, $\tilde{g}\rightarrow\text{q}\bar{\text{q}}\tilde{\chi}_{1}^{0}$ ($m_{\tilde{g}} = 1000$ GeV, $m_{\tilde{\chi}_{1}^{0}} = 800$ GeV).

Expected numbers of signal events in search bins with $7\geq N_{\text{jets}}\geq8$ region for the benchmark models $\tilde{g}\tilde{g}$, $\tilde{g}\rightarrow\text{b}\bar{\text{b}}\tilde{\chi}_{1}^{0}$ ($m_{\tilde{g}} = 1500$ GeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV), $\tilde{g}\tilde{g}$, $\tilde{g}\rightarrow\text{t}\bar{\text{t}}\tilde{\chi}_{1}^{0}$ ($m_{\tilde{g}} = 1500$ GeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV), and $\tilde{g}\tilde{g}$, $\tilde{g}\rightarrow\text{q}\bar{\text{q}}\tilde{\chi}_{1}^{0}$ ($m_{\tilde{g}} = 1000$ GeV, $m_{\tilde{\chi}_{1}^{0}} = 800$ GeV).

Expected numbers of signal events in search bins with $N_{\text{jets}}\geq9$ region for the benchmark models $\tilde{g}\tilde{g}$, $\tilde{g}\rightarrow\text{b}\bar{\text{b}}\tilde{\chi}_{1}^{0}$ ($m_{\tilde{g}} = 1500$ GeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV), $\tilde{g}\tilde{g}$, $\tilde{g}\rightarrow\text{t}\bar{\text{t}}\tilde{\chi}_{1}^{0}$ ($m_{\tilde{g}} = 1500$ GeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV), and $\tilde{g}\tilde{g}$, $\tilde{g}\rightarrow\text{q}\bar{\text{q}}\tilde{\chi}_{1}^{0}$ ($m_{\tilde{g}} = 1000$ GeV, $m_{\tilde{\chi}_{1}^{0}} = 800$ GeV).