Covariance for spin correlation coefficients in beam basis for $p_{T}(t)$ vs. $|cos(\theta)|$ bins

Results for Eigenstate decomposition for $|\Phi^{+} \rangle$ in Helicity basis in $m(t\bar{t})$ vs. $|cos(\theta)|$ bins

Results for Eigenstate decomposition for $|\Phi^{-} \rangle$ in Helicity basis in $m(t\bar{t})$ vs. $|cos(\theta)|$ bins

Results for Eigenstate decomposition for $|\Psi^{+} \rangle$ in Helicity basis in $m(t\bar{t})$ vs. $|cos(\theta)|$ bins

Results for Eigenstate decomposition for $|\Psi^{-} \rangle$ in Helicity basis in $m(t\bar{t})$ vs. $|cos(\theta)|$ bins

Results for Eigenstate decomposition for $|\uparrow \uparrow \rangle$ in Beam basis in $m(t\bar{t})$ vs. $|cos(\theta)|$ bins

Results for Eigenstate decomposition for $|\downarrow \downarrow \rangle$ in Beam basis in $m(t\bar{t})$ vs. $|cos(\theta)|$ bins

Results for Eigenstate decomposition for $|\Psi^{+} \rangle$ in Beam basis in $m(t\bar{t})$ vs. $|cos(\theta)|$ bins

Results for Eigenstate decomposition for $|\Psi^{-} \rangle$ in Beam basis in $m(t\bar{t})$ vs. $|cos(\theta)|$ bins

Results for Purity in Helicity basis in $m(t\bar{t})$ vs. $|cos(\theta)|$ bins

Results for Purity in Beam basis in $m(t\bar{t})$ vs. $|cos(\theta)|$ bins

Results for Entropy in Helicity basis in $m(t\bar{t})$ vs. $|cos(\theta)|$ bins

Results for Entropy in Beam basis in $m(t\bar{t})$ vs. $|cos(\theta)|$ bins

Results for Entanglement in Helicity basis in $m(t\bar{t})$ vs. $|cos(\theta)|$ bins

Results for Entanglement in Beam basis in $m(t\bar{t})$ vs. $|cos(\theta)|$ bins

Results for Eigenstate decomposition for $|\Phi^{+} \rangle$ in Helicity basis in $p_{T}(t)$ vs. $|cos(\theta)|$ bins

Results for Eigenstate decomposition for $|\Phi^{-} \rangle$ in Helicity basis in $p_{T}(t)$ vs. $|cos(\theta)|$ bins

Results for Eigenstate decomposition for $|\Psi^{+} \rangle$ in Helicity basis in $p_{T}(t)$ vs. $|cos(\theta)|$ bins

Results for Eigenstate decomposition for $|\Psi^{-} \rangle$ in Helicity basis in $p_{T}(t)$ vs. $|cos(\theta)|$ bins

Results for Eigenstate decomposition for $|\uparrow \uparrow \rangle$ in Beam basis in $p_{T}(t)$ vs. $|cos(\theta)|$ bins

Results for Eigenstate decomposition for $|\downarrow \downarrow \rangle$ in Beam basis in $p_{T}(t)$ vs. $|cos(\theta)|$ bins

Results for Eigenstate decomposition for $|\Psi^{+} \rangle$ in Beam basis in $p_{T}(t)$ vs. $|cos(\theta)|$ bins

Results for Eigenstate decomposition for $|\Psi^{-} \rangle$ in Beam basis in $p_{T}(t)$ vs. $|cos(\theta)|$ bins

Results for Purity in Helicity basis in $p_{T}(t)$ vs. $|cos(\theta)|$ bins

Results for Purity in Beam basis in $p_{T}(t)$ vs. $|cos(\theta)|$ bins

Results for Entropy in Helicity basis in $p_{T}(t)$ vs. $|cos(\theta)|$ bins

Results for Entropy in Beam basis in $p_{T}(t)$ vs. $|cos(\theta)|$ bins

Results for Entanglement in Helicity basis in $p_{T}(t)$ vs. $|cos(\theta)|$ bins

Results for Entanglement in Beam basis in $p_{T}(t)$ vs. $|cos(\theta)|$ bins

The postfit pDNN distributions in the SR $\mu$ 3b3j assuming $m_{H^\pm} = 600$ GeV. Postfit signal for $m_{H^\pm} = 600$ GeV is also shown. Beneath plot the ratio of data to predictions is shown.

The postfit pDNN distributions in the SR e 2b2j assuming $m_{H^\pm} = 1$ TeV. The signal for $m_{H^\pm} = 1$ TeV is also shown before the fit, assuming a cross section of 1 pb.Beneath plot the ratio of data to predictions is shown.

The postfit pDNN distributions in the SR $\mu$ 2b2j assuming $m_{H^\pm} = 1$ TeV. The signal for $m_{H^\pm} = 1$ TeV is also shown before the fit, assuming a cross section of 1 pb.Beneath plot the ratio of data to predictions is shown.

The postfit pDNN distributions in the SR e 3b3j assuming $m_{H^\pm} = 1$ TeV. The signal for $m_{H^\pm} = 1$ TeV is also shown before the fit, assuming a cross section of 1 pb.Beneath plot the ratio of data to predictions is shown.

The postfit pDNN distributions in the SR $\mu$ 3b3j assuming $m_{H^\pm} = 1$ TeV. The signal for $m_{H^\pm} = 1$ TeV is also shown before the fit, assuming a cross section of 1 pb.Beneath plot the ratio of data to predictions is shown.

Observed and expected 95% CL limits on the cross section times branching ratio for different $H^{\pm}$ mass hypotheses with $ρ_{tc} = 0.4$, $ρ_{tt} = 0.6$ . The inner (green) band and the outer (yellow) band represent the regions containing 68% and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. Theoretical prediction with different $ρ_{tc}$ and $ρ_{tt}$ couplingsare shown with the grey bands representing the corresponding uncertainties due to the factorization and renormalization scales and the PDFs.

Observed excluded phase-space regions, based on 95% CL limits, as a function of $m_{H^{\pm}}$ and $ρ_{tc}$ for various assumed $ρ_{tt}$ values represented with different colors. The limits are extracted from the pDNN distributions based on g2HDM assuming all the extra Yukawa couplings except $ρ_{tt}$ and $ρ_{tc}$ are zero. The results are obtained from the 95% CL limits on the cross section times branching fraction in Fig 6.

Expected excluded phase-space regions, based on 95% CL limits, as a function of $m_{H^{\pm}}$ and $ρ_{tc}$ for various assumed $ρ_{tt}$ values represented with different colors. The limits are extracted from the pDNN distributions based on g2HDM assuming all the extra Yukawa couplings except $ρ_{tt}$ and $ρ_{tc}$ are zero. The results are obtained from the 95% CL limits on the cross section times branching fraction in Fig 6.

Observed excluded phase-space regions, based on 95% CL limits, as a function of $m_{H^{\pm}}$ and $ρ_{tc}$ for various assumed $ρ_{tt}$ values represented with different colors. The limits are extracted from the pDNN distributions based on g2HDM assuming all the extra Yukawa couplings except $ρ_{tt}$ and $ρ_{tc}$ are zero. The results are obtained from the 95% CL limits on the cross section times branching fraction in Fig 6.

Expected excluded phase-space regions, based on 95% CL limits, as a function of $m_{H^{\pm}}$ and $ρ_{tc}$ for various assumed $ρ_{tt}$ values represented with different colors. The limits are extracted from the pDNN distributions based on g2HDM assuming all the extra Yukawa couplings except $ρ_{tt}$ and $ρ_{tc}$ are zero. The results are obtained from the 95% CL limits on the cross section times branching fraction in Fig 6.

Observed excluded phase-space regions, based on 95% CL limits, as a function of $m_{H^{\pm}}$ and $ρ_{tc}$ for various assumed $ρ_{tt}$ values represented with different colors. The limits are extracted from the pDNN distributions based on g2HDM assuming all the extra Yukawa couplings except $ρ_{tt}$ and $ρ_{tc}$ are zero. The results are obtained from the 95% CL limits on the cross section times branching fraction in Fig 6.

Expected excluded phase-space regions, based on 95% CL limits, as a function of $m_{H^{\pm}}$ and $ρ_{tc}$ for various assumed $ρ_{tt}$ values represented with different colors. The limits are extracted from the pDNN distributions based on g2HDM assuming all the extra Yukawa couplings except $ρ_{tt}$ and $ρ_{tc}$ are zero. The results are obtained from the 95% CL limits on the cross section times branching fraction in Fig 6.

Observed excluded phase-space regions, based on 95% CL limits, as a function of $m_{H^{\pm}}$ and $ρ_{tc}$ for various assumed $ρ_{tt}$ values represented with different colors. The limits are extracted from the pDNN distributions based on g2HDM assuming all the extra Yukawa couplings except $ρ_{tt}$ and $ρ_{tc}$ are zero. The results are obtained from the 95% CL limits on the cross section times branching fraction in Fig 6.

Expected excluded phase-space regions, based on 95% CL limits, as a function of $m_{H^{\pm}}$ and $ρ_{tc}$ for various assumed $ρ_{tt}$ values represented with different colors. The limits are extracted from the pDNN distributions based on g2HDM assuming all the extra Yukawa couplings except $ρ_{tt}$ and $ρ_{tc}$ are zero. The results are obtained from the 95% CL limits on the cross section times branching fraction in Fig 6.

Observed excluded phase-space regions, based on 95% CL limits, as a function of $m_{H^{\pm}}$ and $ρ_{tc}$ for various assumed $ρ_{tt}$ values represented with different colors. The limits are extracted from the pDNN distributions based on g2HDM assuming all the extra Yukawa couplings except $ρ_{tt}$ and $ρ_{tc}$ are zero. The results are obtained from the 95% CL limits on the cross section times branching fraction in Fig 6.

Expected excluded phase-space regions, based on 95% CL limits, as a function of $m_{H^{\pm}}$ and $ρ_{tc}$ for various assumed $ρ_{tt}$ values represented with different colors. The limits are extracted from the pDNN distributions based on g2HDM assuming all the extra Yukawa couplings except $ρ_{tt}$ and $ρ_{tc}$ are zero. The results are obtained from the 95% CL limits on the cross section times branching fraction in Fig 6.

Two times the difference of the negative log-likelihood (NLL) as a function of σ(pp →(q$^{light}$)H$^{\pm}$) B(H$^{\pm}$ → tb, t → blν) and σ(pp → bH$^{\pm}$)B(H$^{\pm}$ → tb, t → blν) with the best fit point extracted from the pDNN distribution for the $H^{\pm}$ mass $m_{H^{\pm}} = 600$ GeV. The solid and dashed contours represent the observed and expected limits, respectively. The points along the contours represent 68% and 95% CL regions, extracted from the 2∆NLL values of 2.30 and 5.99, respectively. The g2HDM predictions are also shown. The points along the g2HDM prediction line represent different $ρ_{tt}$, $ρ_{tc}$ coupling sets, with all other extra Yukawa couplings assumed to be zero.

Two times the difference of the negative log-likelihood (NLL) as a function of σ(pp →(q$^{light}$)H$^{\pm}$) B(H$^{\pm}$ → tb, t → blν) and σ(pp → bH$^{\pm}$)B(H$^{\pm}$ → tb, t → blν) with the best fit point extracted from the pDNN distribution for the $H^{\pm}$ mass $m_{H^{\pm}} = 600$ GeV. The solid and dashed contours represent the observed and expected limits, respectively. The points along the contours represent 68% and 95% CL regions, extracted from the 2∆NLL values of 2.30 and 5.99, respectively. The g2HDM predictions are also shown. The points along the g2HDM prediction line represent different $ρ_{tt}$, $ρ_{tc}$ coupling sets, with all other extra Yukawa couplings assumed to be zero.

Two times the difference of the negative log-likelihood (NLL) as a function of σ(pp →(q$^{light}$)H$^{\pm}$) B(H$^{\pm}$ → tb, t → blν) and σ(pp → bH$^{\pm}$)B(H$^{\pm}$ → tb, t → blν) with the best fit point extracted from the pDNN distribution for the $H^{\pm}$ mass $m_{H^{\pm}} = 1000$ GeV. The solid and dashed contours represent the observed and expected limits, respectively. The points along the contours represent 68% and 95% CL regions, extracted from the 2∆NLL values of 2.30 and 5.99, respectively. The g2HDM predictions are also shown. The points along the g2HDM prediction line represent different $ρ_{tt}$, $ρ_{tc}$ coupling sets, with all other extra Yukawa couplings assumed to be zero.

Two times the difference of the negative log-likelihood (NLL) as a function of σ(pp →(q$^{light}$)H$^{\pm}$) B(H$^{\pm}$ → tb, t → blν) and σ(pp → bH$^{\pm}$)B(H$^{\pm}$ → tb, t → blν) with the best fit point extracted from the pDNN distribution for the $H^{\pm}$ mass $m_{H^{\pm}} = 1000$ GeV. The solid and dashed contours represent the observed and expected limits, respectively. The points along the contours represent 68% and 95% CL regions, extracted from the 2∆NLL values of 2.30 and 5.99, respectively. The g2HDM predictions are also shown. The points along the g2HDM prediction line represent different $ρ_{tt}$, $ρ_{tc}$ coupling sets, with all other extra Yukawa couplings assumed to be zero.

Observed (solid black line) and median expected (dashed black line) upper limits at 95% CL on the signal cross section times branching fraction $\sigma (\mathrm{pp} \to \mathrm{H})\,\mathcal{B}(\mathrm{H} \to \mathcal{A} \mathcal{A} \to 4 \gamma)$ at various $m_A$ mass points. The 68 and 95% confidence intervals around the median expected limit are shown as the green and yellow bands, respectively.

Relative cutflow fractions for representative signal samples, normalized to the number of events before any selection (All events = 1.0). The selections correspond to the trigger requirement, photon multiplicity and identification, kinematic requirements on the three-photon system, a selection on the angular separation (dR) of photons consistent with a merged topology, and a multivariate (MVA) discriminator and a charged-hadron isolation requirement applied to photon candidates used to suppress jet background. The final SR corresponds to the combination of the $m_H$ and $m_A$ signal regions used in the final fit.

Expected and observed 95% CL upper limits on the V' boson production cross section as functions of the resonance mass mV' in the HVT model A.

Expected and observed 95% CL upper limits on the V' boson production cross section as functions of the resonance mass mV' in the HVT model B.

Expected and observed 95% CL upper limits on the V' boson production cross section as functions of the resonance mass mV' in the HVT model C.

Expected and observed 95% CL exclusion contour in gF-gH parameter space for the DIQ channel at 3 TeV. Here, model A corresponds to (0.556, 0.562) and model B to (-2.928, 0.146)

Expected and observed 95% CL exclusion contour in gF-gH parameter space for the DIQ channel at 4 TeV.

Expected and observed 95% CL exclusion contour in gF-gH parameter space for the DIL channel at 3 TeV.

Expected and observed 95% CL exclusion contour in gF-gH parameter space for the DIL channel at 4 TeV.

Expected and observed 95% CL exclusion contour in gH-gF parameter space for the DIB channel at 3 TeV.

Expected and observed 95% CL exclusion contour in gH-gF parameter space for the DIB channel at 4 TeV.

Expected and observed 95% CL exclusion contour in gF-gH parameter space for the Full Combined channel at 3 TeV.

Expected and observed 95% CL exclusion contour in gF-gH parameter space for the Full Combined channel at 4 TeV.

Expected and observed 95% CL exclusion contour in g_q12-g_q3 parameter space at m_V' = 3 TeV. Uncertainty bands shown as separate contour lines. Channels not in combination are shown separately.

Expected and observed 95% CL exclusion contour in g_q12-g_q3 parameter space at m_V' = 4 TeV. Uncertainty bands shown as separate contour lines. Channels not in combination are shown separately.

Expected and observed 95% CL upper limits on the W' production cross section in the HVT Model A benchmark, compared with the theoretical prediction.

Expected and observed 95% CL upper limits on the W' production cross section in the HVT Model B benchmark, compared with the theoretical prediction.

Expected and observed 95% CL upper limits on the W' production cross section in the HVT Model C benchmark, compared with the theoretical prediction.

Expected and observed 95% CL upper limits on the Z' production cross section in the HVT Model A benchmark, compared with the theoretical prediction.

Expected and observed 95% CL upper limits on the Z' production cross section in the HVT Model B benchmark, compared with the theoretical prediction.

Expected and observed 95% CL upper limits on the Z' production cross section in the HVT Model C benchmark, compared with the theoretical prediction.

Measured and predicted cross sections (fb) in bins of photon transverse momentum $p_{T}^{\gamma}$. The fiducial phase space definition follows the analysis selection in the paper.

Expected and observed 95% CL intervals for anomalous coupling parameters. All other anomalous coupling parameters are fixed to zero when extracting each interval.