Normalized differential cross sections, compared with the SM predictions, as a function of the absolute value of pseudorapidity of the subleading jet in transverse momentum. The SM predictions are obtained using Madgraph5_aMC@NLO version 2.6.5 at NLO in QCD with PYTHIA version 8.240

Normalized differential cross sections, compared with the SM predictions, as a function of the invariant mass of the two jets. The SM predictions are obtained using Madgraph5_aMC@NLO version 2.6.5 at NLO in QCD with PYTHIA version 8.240

Normalized differential cross sections, compared with the SM predictions, as a function of the photon's transverse momentum. The SM predictions are obtained using Madgraph5_aMC@NLO version 2.6.5 at NLO in QCD with PYTHIA version 8.240

Normalized differential cross sections, compared with the SM predictions, as a function of the event centrality as defined in Ref.[13] of the paper. The SM predictions are obtained using Madgraph5_aMC@NLO version 2.6.5 at NLO in QCD with PYTHIA version 8.240

Normalized differential cross sections, compared with the SM predictions, as a function of the Zeppenfeld variable that is defined in Eq.(1) of the paper. The SM predictions are obtained using Madgraph5_aMC@NLO version 2.6.5 at NLO in QCD with PYTHIA version 8.240

Negative of twice in the difference in the log-likelihood as a function of $c_{\mathrm{W}}$, based on 138 $\mathrm{fb}^{-1}$ of CMS data at 13 TeV.

Negative of twice in the difference in the log-likelihood as a function of $c_{\mathrm{HWB}}$, based on 138 $\mathrm{fb}^{-1}$ of CMS data at 13 TeV.

Negative of twice in the difference in the log-likelihood as a function of $c_{\mathrm{HWB}}$ and $c_{\mathrm{W}}$, based on 138 $\mathrm{fb}^{-1}$ of CMS data at 13 TeV.

Predicted and measured EWK $\gamma\mathrm{jj}$ fiducial cross sections.

Normalized differential distributions of the $\Delta R (\gamma,\mathrm{tt})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Absolute differential distributions of the min. $\Delta R (\gamma,\mathrm{t})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Normalized differential distributions of the min. $\Delta R (\gamma,\mathrm{t})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Absolute differential distributions of the $m(\mathrm{tt})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Normalized differential distributions of the $m(\mathrm{tt})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Absolute differential distributions of the leading lepton $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Normalized differential distributions of the leading lepton $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Absolute differential distributions of the photon $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Normalized differential distributions of the photon $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Absolute differential distributions of the $\Delta \phi (\ell,\ell)$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Normalized differential distributions of the $\Delta \phi (\ell,\ell)$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Absolute differential distributions of the ratio between the cross sections of $\mathrm{tt\gamma}$ and $\mathrm{tt}$, as a function of the leading top quark $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Absolute differential distributions of the ratio between the cross sections of $\mathrm{tt\gamma}$ and $\mathrm{tt}$, as a function of the leading lepton $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the leading top quark $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the leading top quark $p_{\mathrm{T}}$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the $\Delta R (\gamma,\mathrm{tt})$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the $\Delta R (\gamma,\mathrm{tt})$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the min. $\Delta R (\gamma,\mathrm{t})$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the min. $\Delta R (\gamma,\mathrm{t})$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the $m(\mathrm{tt})$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the $m(\mathrm{tt})$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the leading lepton $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the leading lepton $p_{\mathrm{T}}$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the photon $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the photon $p_{\mathrm{T}}$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the $\Delta \phi (\ell,\ell)$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the $\Delta \phi (\ell,\ell)$.

Correlation matrix of the total uncertainty in the absolute differential $\mathrm{tt\gamma}$ / $\mathrm{tt}$ cross section ratio as a function of the leading top $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential $\mathrm{tt\gamma}$ / $\mathrm{tt}$ cross section ratio as a function of the leading top $p_{\mathrm{T}}$.

Correlation matrix of the total uncertainty in the absolute differential $\mathrm{tt\gamma}$ / $\mathrm{tt}$ cross section ratio as a function of the leading lepton $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential $\mathrm{tt\gamma}$ / $\mathrm{tt}$ cross section ratio as a function of the leading lepton $p_{\mathrm{T}}$.

The 95% CL upper limits on the inclusive signal strength $ r = \sigma_{HH} /\sigma_{HH} ^\text{SM} $ for each channel and their combination. The inner (green) band and the outer (yellow) band indicate the 68 and 95% CL intervals, respectively, under the background-only hypothesis. The $ b\overline{b}b\overline{b} $ and $ b\overline{b}WW $ contributions have been combined in order to simplify the presentation of results.

The 95$\%$ CL upper limits on the VBF signal strength $r_{VBF HH} = \sigma_{VBF HH}/\sigma_{VBF HH}^{SM}$ for each channel and their combination. The inner (green) band and the outer (yellow) band indicate the 68 and 95$\%$ CL intervals, respectively, under the background-only hypothesis. Not all of the eight channels are used to extract the combined VBF limit. The contributing channels are indicated in the figure. The b\overline{b}b\overline{b} $ and $ b\overline{b}WW $ contributions have been combined in order to simplify the presentation of results.

The 95$\%$ CL upper limits on the inclusive HH cross section as a function of $\kappa_\lambda$. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH GGF and VBF signal cross sections are not considered because here we directly constrain the measured cross section.The inner (green) band and the outer (yellow) band indicate the 68 and 95$\%$ CL intervals, respectively, under the background-only hypothesis. The star shows the limit at the SM value for $\kappa_\lambda$.

The 95$\%$ CL upper limits on the inclusive HH cross section as a function of $\kappa_{2V}$. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH VBF signal cross sections are not considered because here we directly constrain the measured cross section.The inner (green) band and the outer (yellow) band indicate the 68 and 95$\%$ CL intervals, respectively, under the background-only hypothesis. The star shows the limit at the SM value for $\kappa_{2V}$.

The profile likelihood ratio test statistic $-2\Delta\log(L)$ as a function of coupling modifiers $\kappa_\lambda$ for the combination of all channels, when all the other parameters are fixed to their SM values.

The profile likelihood ratio test statistic $-2\Delta\log(L)$ as a function of coupling modifiers $\kappa_{2V}$ for the combination of all channels, when all the other parameters are fixed to their SM values.

The 68 $\%$ CL expected contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ CL expected contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 68 $\%$ CL observed contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ observed CL contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 5 standard deviation expected contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 5 standard deviation observed contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The observed best-fit values in the ($\kappa_\lambda$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 68 $\%$ expected CL contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 68$\%$ expected CL contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ CL expected contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ CL observed contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 5 standard deviation expected contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 5 standard deviation observed contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The observed best-fit values in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 68 $\%$ CL expected contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_t$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 68 $\%$ CL observed contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_t$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ CL expected contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_t$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ CL observed contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_t$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The observed best-fit values in the ($\kappa_\lambda$, $\kappa_t$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

Upper limits on the HH production cross section at 95$\%$ CL for the two sets of HEFT benchmarks. The theoretical uncertainties in the HH GGF signal cross section are not considered because we directly constrain the measured cross section.

Upper limits on the HH cross section as a function of the $C_2$ coupling modifier. The theoretical uncertainties in the HH GGF signal cross section are not considered because we directly constrain the measured cross section.

The profile likelihood ratio test statistic $-2\Delta\log(L)$ as a function of the $c_2$ coupling modifier.

The 68$\%$CL expected contours of $-2\Delta\log(L)$ in the ($c_2$,$\kappa_t$) for the combination of all channels when all the other parameters are fixed to their SM value.

The 68$\%$CL observed contours of $-2\Delta\log(L)$ in the ($c_2$,$\kappa_t$) for the combination of all channels when all the other parameters are fixed to their SM value.

The 95$\%$CL expected contours of $-2\Delta\log(L)$ in the ($c_2$,$\kappa_t$) for the combination of all channels when all the other parameters are fixed to their SM value.

The 95$\%$CL observed contours of $-2\Delta\log(L)$ in the ($c_2$,$\kappa_t$) for the combination of all channels when all the other parameters are fixed to their SM value.

The 68 $\%$ CL expected contours of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) for the combination of all channels when all the other parameters are fixed to their SM value.

The 68 $\%$ CL observed contours of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) for the combination of all channels when all the other parameters are fixed to their SM value.

The 95 $\%$ CL expected contours of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) for the combination of all channels when all the other parameters are fixed to their SM value.

The 95 $\%$ CL observed contours of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) for the combination of all channels when all the other parameters are fixed to their SM value.

The $-2\Delta\log(L)$ scan for the combination of all channels as a function of $|\cos{\alpha}|$ for singlet extensions of the SM with spontaneous $Z_2$ symmetry breaking. The $|\cos{\alpha}|$ is constrained at 95\% CL between 0.79 and 1. The ranges of $|\cos{\alpha}|$ and $\lambda_\text{eff}$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown.

The 95$\%$ CL expected contour of $-2\Delta\log(L)$ in the ($|\cos{\alpha}|$, $\lambda_{eff}$) plane, where $\lambda_{eff}$=$\lambda_\alpha$-$\tan{(\alpha)\frac{m_{2}}{\nu}}$. The ranges of $|\cos{\alpha}|$ and $\lambda_{eff}$ are chosen to guarantee the validity of the models. The excluded regions are to the left of the curves shown.

The 95$\%$ CL observed contour of $-2\Delta\log(L)$ in the ($|\cos{\alpha}|$, $\lambda_{eff}$) plane, where $\lambda_{eff}$=$\lambda_\alpha$-$\tan{(\alpha)\frac{m_{2}}{\nu}}$. The ranges of $|\cos{\alpha}|$ and $\lambda_{eff}$ are chosen to guarantee the validity of the models. The excluded regions are to the left of the curves shown.

The 95% CL expected contour of $-2\Delta\log(L)$ for the combination of all channels in the ($m_{H}$, $\tan\beta$) plane for a fixed value of $|\cos(\beta-\alpha)|=0.2$ in the 2HDM-I model. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $m_{H}$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $m_{H}$ > 800 GeV for all models considered.

The 95$\%$ CL observed contour of $-2\Delta\log(L)$ for the combination of all channels in the ($m_{H}$, $\tan\beta$) plane for a fixed value of $|\cos(\beta-\alpha)|=0.2$ in the 2HDM-I model. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $m_{H}$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $m_{H}$ > 800 GeV for all models considered.

The 95$\%$ CL contour of $-2\Delta\log(L)$ for the combination of all channels in the ($m_{H}$, $\tan\beta$) plane for a fixed value of $|\cos(\beta-\alpha)|=0.2$ for the 2HDM I models. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $m_{H}$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $m_{H}$ > 800 GeV for all models considered.

The 95$\%$ CL contour of $-2\Delta\log(L)$ for the combination of all channels in the ($m_{H}$, $\tan\beta$) plane for a fixed value of $|\cos(\beta-\alpha)|=0.2$ for the 2HDM II models. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $m_{H}$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $m_{H}$ > 800 GeV for all models considered.

The 95$\%$ CL contour of $-2\Delta\log(L)$ for the combination of all channels in the ($m_{H}$, $\tan\beta$) plane for a fixed value of $|\cos(\beta-\alpha)|=0.2$ for the 2HDM X models. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $m_{H}$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $m_{H}$ > 800 GeV for all models considered.

The 95$\%$ CL contour of $-2\Delta\log(L)$ for the combination of all channels in the ($m_{H}$, $\tan\beta$) plane for a fixed value of $|\cos(\beta-\alpha)|=0.2$ for the 2HDM Y models. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $m_{H}$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $m_{H}$ > 800 GeV for all models considered.

The 95$\%$ CL expected contour of $-2\Delta\log(L)$ for the combination of all channels in the ($|\cos(\beta-\alpha)|$, $\tan\beta$) plane for a fixed value of $m_{H}=1000GeV$ in the 2HDM-I model. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $|\cos(\beta-\alpha)|$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $|\cos(\beta-\alpha)|$ > 0.16 and $|\cos(\beta-\alpha)|$ < -0.13 for all models considered.

The 95$\%$ CL observed contour of $-2\Delta\log(L)$ for the combination of all channels in the ($|\cos(\beta-\alpha)|$, $\tan\beta$) plane for a fixed value of $m_{H}=1000GeV$ in the 2HDM-I model. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $|\cos(\beta-\alpha)|$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $|\cos(\beta-\alpha)|$ > 0.16 and $|\cos(\beta-\alpha)|$ < -0.13 for all models considered.

The 95% CL contour of $-2\Delta\log(L)$ for the combination of all channels in the ($|\cos(\beta-\alpha)|$, $\tan\beta$) plane for a fixed value of $m_{H}=1000GeV$ for the 2HDM I models. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $|\cos(\beta-\alpha)|$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $|\cos(\beta-\alpha)|$ > 0.16 and $|\cos(\beta-\alpha)|$ < -0.13 for all models considered.

The 95% CL contour of $-2\Delta\log(L)$ for the combination of all channels in the ($|\cos(\beta-\alpha)|$, $\tan\beta$) plane for a fixed value of $m_{H}=1000GeV$ for the 2HDM II models. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $|\cos(\beta-\alpha)|$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $|\cos(\beta-\alpha)|$ > 0.16 and $|\cos(\beta-\alpha)|$ < -0.13 for all models considered.

The 95% CL contour of $-2\Delta\log(L)$ for the combination of all channels in the ($|\cos(\beta-\alpha)|$, $\tan\beta$) plane for a fixed value of $m_{H}=1000GeV$ for the 2HDM X models. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $|\cos(\beta-\alpha)|$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $|\cos(\beta-\alpha)|$ > 0.16 and $|\cos(\beta-\alpha)|$ < -0.13 for all models considered.

The 95% CL contour of $-2\Delta\log(L)$ for the combination of all channels in the ($|\cos(\beta-\alpha)|$, $\tan\beta$) plane for a fixed value of $m_{H}=1000GeV$ for the 2HDM Y models. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $|\cos(\beta-\alpha)|$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $|\cos(\beta-\alpha)|$ > 0.16 and $|\cos(\beta-\alpha)|$ < -0.13 for all models considered.

The $-2\Delta\log(L)$ scan for the combination of all channels as a function of $\xi$ for MCHM$_4$. The range of $\xi$ is chosen to guarantee the validity of the model. At 95% CL, $\xi$ is constrained to be between 0 and 0.45 in MCHM$_4$.

The $-2\Delta\log(L)$ scan for the combination of all channels as a function of $\xi$ for MCHM$_5$. The range of $\xi$ is chosen to guarantee the validity of the model. At 95% CL, $\xi$ is constrained to be between 0 and 0.26 in MCHM$_5$.

Treatment of most important common systematic uncertainties in the S2 scenario.

The expected upper limits at 95% CL on the HH signal strength from the combination of all the considered channels at different integrated luminosities for the S2 scenario.

The expected upper limits at 95% CL on the HH signal strength from the combination of all the considered channels under different assumptions on the systematic uncertainties for an integrated luminosity of 3000fb$^{-1}$ (right).

The expected $-2\Delta\log(L)$ scan as a function of coupling modifier $\kappa_\lambda$ for the combination of all contributing channels, projected to different integrated luminosities. All other parameters are fixed to their SM values.

The expected $-2\Delta\log(L)$ scan as a function of coupling modifier $\kappa_\lambda$ for the combination of all contributing channels, projected under different assumptions on the systematic uncertainties for an integrated luminosity of 3000fb$^{-1}$. All other parameters are fixed to their SM values.

Expected significance for the HH signal projected to 2000 or 3000 fb$^{-1}$ under different assumptions of systematic uncertainties.

The expected signal significance as a function of integrated luminosity for the nominal scenario of systematic uncertainties S2 and for the scenario with statistical uncertainties only.

The expected signal significance as a function of $\kappa_\lambda$ under different assumptions on the systematic uncertainties for an integrated luminosity of 3000fb$^{-1}$.

The 95$\%$ CL upper limits on the inclusive HH signal strenngth as a function of $\kappa_{\lambda}$ coupling modifier. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH VBF signal cross sections are not considered because here we directly constrain the measured cross section.The inner (green) band and the outer (yellow) band indicate the 68 and 95$\%$ CL intervals, respectively, under the background-only hypothesis.

The 95$\%$ CL upper limits on the inclusive HH signal strenngth as a function of $\kappa_{2V}$ coupling modifier. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH VBF signal cross sections are not considered because here we directly constrain the measured cross section.The inner (green) band and the outer (yellow) band indicate the 68 and 95$\%$ CL intervals, respectively, under the background-only hypothesis.

The $-2\Delta\log(L)$ scan as a function of $\kappa_\lambda$ for all channels, when all the other parameters are fixed to their SM value.

The $-2\Delta\log(L)$ scan as a function of $\kappa_{2V}$ for all channels, when all the other parameters are fixed to their SM value.

The $-2\Delta\log(L)$ scan as a function of $\kappa_\lambda$ for all channels, when all the other parameters are fixed to their SM value.

The 95$\%$ CL upper limits on the inclusive HH cross section as a function of $\kappa_\lambda$ coupling modifier for all channels. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH GGF and VBF signal cross sections are not considered because here we directly constrain the measured cross section.

The 95$\%$ CL upper limits on the inclusive VBF HH cross section as a function of $\kappa_{2V}$ coupling modifier for all channels. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH GGF and VBF signal cross sections are not considered because here we directly constrain the measured cross section.

The 95$\%$ CL upper limits on the inclusive HH signal strength as a function of $\kappa_\lambda$ coupling modifier for all channels. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH GGF and VBF signal cross sections are not considered because here we directly constrain the measured cross section.

The 95$\%$ CL upper limits on the inclusive VBF HH signal strength as a function of $\kappa_{2V}$ coupling modifier for all channels. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH GGF and VBF signal cross sections are not considered because here we directly constrain the measured cross section.

The best fit value for $\kappa_\lambda$ compared to the SM expectation for all channels and their combination, when all the other parameters are fixed to their SM values.

The best fit value for $\kappa_2V$ compared to the SM expectation for all channels and their combination, when all the other parameters are fixed to their SM values.

The 68$\%$ CL expected contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 68$\%$ CL observed contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ CL expected contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ CL observed contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 5 $\sigma$ expected contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 5 $\sigma$ observed contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The best-fit values of ($\kappa_{V}$, $\kappa_{2V}$) for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ CL upper limits on the inclusive VBF HH signal strength as a function of $\kappa_{2V}$ coupling modifier for all channels. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH GGF and VBF signal cross sections are not considered because here we directly constrain the measured cross section. The inner (green) band and the outer (yellow) band indicate the 68 and 95% CL intervals, respectively, under the background-only hypothesis.

The $-2\Delta\log(L)$ scan as a function of the $c_2$ coupling modifier for all channels, when all the other parameters are fixed to their SM values.

The 95% CL upper limits on the inclusive HH cross section as a function of the $c_2$ coupling modifier for all channels. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH ggF signal cross sections are not considered because we directly constrain the measured cross section.

The 95% CL upper limits on the inclusive HH signal strength as a function of the $c_2$ coupling modifier for all channels. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH ggF signal cross sections are considered in this case

The best fit value for the $c_2$ coupling modifier compared to the SM expectation for all channels and their combination, when all the other parameters are fixed to their SM values.

The observed 68 $\%$ CL observed contours ($\kappa_\lambda$ =1.0) of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_t$) plane for the combination of all channels when $\kappa_t$ is allowed to vary. The range of $\kappa_t$ is set between 0.5 and 1.5 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The observed 68 $\%$ CL observed contours (floating $\kappa_\lambda$) of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_t$) plane for the combination of all channels when $\kappa_t$ is allowed to vary. The range of $\kappa_t$ is set between 0.5 and 1.5 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The observed 95 $\%$ CL observed contours ($\kappa_\lambda$ =1.0) of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_t$) plane for the combination of all channels when $\kappa_t$ is allowed to vary. The range of $\kappa_t$ is set between 0.5 and 1.5 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The observed 95 $\%$ CL observed contours (floating $\kappa_\lambda$) of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_t$) plane for the combination of all channels when $\kappa_t$ is allowed to vary. The range of $\kappa_t$ is set between 0.5 and 1.5 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The best fit values ($\kappa_\lambda$ =1.0) of ($c_2$, $\kappa_t$) for the combination of all channels when $\kappa_t$ is allowed to vary. The range of $\kappa_t$ is set between 0.5 and 1.5 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The best fit values (floating $\kappa_\lambda$) of ($c_2$, $\kappa_t$) for the combination of all channels when $\kappa_t$ is allowed to vary. The range of $\kappa_t$ is set between 0.5 and 1.5 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The observed 68$\%$ CL cobsevred contours ($\kappa_{t}$ =1.0 (SM)) of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) plane for the combination of all channels when $\kappa_\lambda$ is allowed to vary. The range of $\kappa_\lambda$ is set between -15 and 15 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The observed 68$\%$ CL cobsevred contours (floating $\kappa_{t}$) of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) plane for the combination of all channels when $\kappa_\lambda$ is allowed to vary. The range of $\kappa_\lambda$ is set between -15 and 15 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The observed 95$\%$ CL cobsevred contours ($\kappa_{t}$ =1.0 (SM)) of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) plane for the combination of all channels when $\kappa_\lambda$ is allowed to vary. The range of $\kappa_\lambda$ is set between -15 and 15 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The observed 95$\%$ CL cobsevred contours (floating $\kappa_{t}$) of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) plane for the combination of all channels when $\kappa_\lambda$ is allowed to vary. The range of $\kappa_\lambda$ is set between -15 and 15 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The best fit values ($\kappa_{t}$ =1.0 (SM)) of ($c_2$, $\kappa_\lambda$) for the combination of all channels when $\kappa_\lambda$ is allowed to vary. The range of $\kappa_\lambda$ is set between -15 and 15 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The best fit values (floating $\kappa_{t}$) for ($c_2$, $\kappa_\lambda$) for the combination of all channels when $\kappa_\lambda$ is allowed to vary. The range of $\kappa_\lambda$ is set between -15 and 15 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The observed 68% CL contours of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) plane for the combination of all channels when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed best fit values in the ($c_2$, $\kappa_\lambda$) plane for the combination of all channels when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed 68% CL contours of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) plane for the $b\bar{b}b\bar{b}$ resolved channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed best fit values in the ($c_2$, $\kappa_\lambda$) plane for the $b\bar{b}b\bar{b}$ resolved channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed 68% CL contours of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) plane for the $b\bar{b}b\bar{b}$ merged channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed best fit values in the ($c_2$, $\kappa_\lambda$) plane for the $b\bar{b}b\bar{b}$ merged channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed 68% CL contours of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) plane for the $b\bar{b}\gamma\gamma$ channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed best fit values in the ($c_2$, $\kappa_\lambda$) plane for the $b\bar{b}\gamma\gamma$ channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed 68% CL contours of $-2\Delta\log(L)$ in in the ($c_2$, $\kappa_\lambda$) plane for the $b\bar{b}W^{+}W^{-}$ channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The best fit values in in the ($c_2$, $\kappa_\lambda$) plane for the $b\bar{b}W^{+}W^{-}$ channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed 68% CL contours of $-2\Delta\log(L)$ in in the ($c_2$, $\kappa_\lambda$) plane for the multilepton channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The best fit values in in the ($c_2$, $\kappa_\lambda$) plane for the multilepton channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed 68% CL contours of $-2\Delta\log(L)$ in in the ($c_2$, $\kappa_\lambda$) plane for the $W^{+}W^{-}\gamma\gamma$ channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The best fit values in in the ($c_2$, $\kappa_\lambda$) plane for the $W^{+}W^{-}\gamma\gamma$ channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

Signal strength and significance for $\text{t}\overline{\text{t}}\text{Z}(\text{Z}\to\text{b}\overline{\text{b}})$ decays with respect to the standard model expectation of unity.

Signal strength and significance for $\text{t}\overline{\text{t}}\text{Z}(\text{Z}\to\text{c}\overline{\text{c}})$ decays with respect to the standard model expectation of unity.

2D likelihood scans for $\kappa_{b}$ and $\kappa_{c}$ for $\text{t}\overline{\text{t}}\text{H}(\text{H}\to\text{c}\overline{\text{c}})$ and $\text{t}\overline{\text{t}}\text{H}(\text{H}\to\text{b}\overline{\text{b}})$.

The nuclear suppression factor $S^{\mathrm{J}/\psi}$ of incoherent $\mathrm{J}/\psi$ meson photoproduction as a function of Bjorken $x$ in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $x$ values are evaluated at the centers of their corresponding rapidity ranges. The theoretical uncertainties are due to the uncertainties in the photon flux.

The ratio between incoherent and coherent $\mathrm{J}/\psi$ meson photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$, measured in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$ values used correspond to the center of each rapidity range. The theoretical uncertainties is due to the uncertainties in the photon flux.

The total covariance matrix of the total incoherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The experimental covariance matrix of the total incoherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The theoretical (photon flux) covariance matrix of the total incoherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

Normalized $\sigma_{\psi(2S)}/\sigma_{J/\psi}$ in $6.5<p_T<30.0\,GeV$ and $ 1 < y_{CM} < 1.935$ as functions of normalized $\text{N}^{{\text{corr.}}}_{\text{track}}$

Normalized $\sigma_{\psi(2S)}/\sigma_{J/\psi}$ in $6.5<p_T<30.0\,GeV$ and $ -2.865 < y_{CM} < 1.935$ as functions of normalized $\text{N}^{{\text{corr.}}}_{\text{track}}$

Normalized $\sigma_{\psi(2S)}/\sigma_{J/\psi}$ in $3.0<p_T<6.5\,GeV$ and $ -2.865 < y_{CM} < -2 $ as functions of normalized $\text{N}^{{\text{corr.}}}_{\text{track}}$

Normalized $\sigma_{\psi(2S)}/\sigma_{J/\psi}$ in $3.0<p_T<6.5\,GeV$ and $ 1 < y_{CM} < 1.935 $ as functions of normalized $\text{N}^{{\text{corr.}}}_{\text{track}}$

Distributions for data and expected signal and background contributions of the MVA score for the µ + jets channel in the 3j1b category, before the maximum likelihood fit. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty in the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-prediction ratio. The first and last bins in each distribution include underflow and overflow events, respectively.

Observed and predicted number of events in each of the eight categories of the signal region, before the maximum likelihood fit. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty in the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-prediction ratio.

Distributions for the e + jets final state before the maximum likelihood fit, MVA score bins for the 3j1b category and ∆R med ( j, j 0 ) bins for the other categories. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty of the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-simulation ratio, where the simulations are adjusted according to their normalization before the fit. The first and last bins in each distribution include underflow and overflow events, respectively.

Distributions for the e + jets final state after the maximum likelihood fit, MVA score bins for the 3j1b category and ∆R med ( j, j 0 ) bins for the other categories. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty of the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-simulation ratio, where the simulations are adjusted according to their normalization after the fit. The first and last bins in each distribution include underflow and overflow events, respectively.

Distributions for the µ + jets final state before the maximum likelihood fit, MVA score bins for the 3j1b category and ∆R med ( j, j 0 ) bins for the other categories. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty of the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-simulation ratio, where the simulations are adjusted according to their normalization before the fit. The first and last bins in each distribution include underflow and overflow events, respectively.

Distributions for the µ + jets final state after the maximum likelihood fit, MVA score bins for the 3j1b category and ∆R med ( j, j 0 ) bins for the other categories. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty of the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-simulation ratio, where the simulations are adjusted according to their normalization after the fit. The first and last bins in each distribution include underflow and overflow events, respectively.

Covariance matrix for the lepton+jets ML fit.

Covariance matrix for the combination ML fit.

Distributions in $S_T$ in the SR for the electron channel, after a background-only fit to the data. The signal distributions are scaled to the cross section predicted by the theory. The hatched bands show the post-fit uncertainty band, combining all sources of uncertainty. The ratio of data to the background predictions is shown in the panels below the distributions.

Correlations of the ttH signal-strength modifiers in the different Higgs pT bins of the STXS measurement.

Observed and expected upper 95% CL limit on the tH signal strength modifier $\mu_{tH}$ with $\mu_{ttH}=1$ for different channels and years, where the uncertainties are uncorrelated between the channels and years, and in their combination. The one and two standard deviation confidence intervals on the expected limit are also listed (sigma).

Likelihood-ratio test statistic as a function of the ttH and tH signal strength modifiers $\mu_{ttH}$ and $\mu_{tH}$. The observed best fit point is $(\mu_{tH}, \mu_{ttH}) = (-3.83, +0.35)$.

Observed and expected likelihood ratio test statistic as a function of $\kappa_{t}$ and $\kappa_{V}$. The observed best fit point is $(\kappa_{t}, \kappa_{V}) = (+0.59, +1.40)$.

Observed and expected likelihood ratio test statistic as a function of $\kappa_{t}$ for $\kappa_{V}=1$. The observed best fit point is $\kappa_{t} = +0.54$.

Observed and expected likelihood ratio test statistic as a function of $\kappa_{t}$ and $\tilde{\kappa}_{t}$ with $\kappa_{V}=1$. The observed best fit point is $(\kappa_{t}, \tilde{\kappa}_{t}) = (+0.53, 0.00)$.

Observed and expected likelihood ratio test statistic as a function of $f_{CP}$. The overall modifier of the top-Higgs coupling strength is profiled such that the likelihood attains its minimum at each point in the plane, and $\kappa_{V}$ is set to 1. The observed best fit point is $f_{CP} = 0.00$.

Observed and expected likelihood ratio test statistic as a function of $\cos(\alpha)$. The overall modifier of the top-Higgs coupling strength is profiled such that the likelihood attains its minimum at each point in the plane, and $\kappa_{V}$ is set to 1. The observed best fit point is $\cos(\alpha) = +1.00$.

Observed and expected likelihood ratio test statistic as a function of $\kappa_{t}$ and $\tilde{\kappa}_{t}$ with $\kappa_{V}=1$, in the $H \rightarrow bb$, $H\rightarrow\gamma\gamma$, $H\rightarrow ZZ$, $H\rightarrow WW$, and $H\rightarrow\tau\tau$ decay channels and in their combination.

Observed and expected likelihood ratio test statistic as a function of $f_{CP}$, in the $H \rightarrow bb$, $H\rightarrow\gamma\gamma$, $H\rightarrow ZZ$, $H\rightarrow WW$, and $H\rightarrow\tau\tau$ decay channel and in their combination. The overall modifier of the top-Higgs coupling strength is profiled such that the likelihood attains its minimum at each point in the plane, and $\kappa_{V}$ is set to 1.

Observed and expected likelihood ratio test statistic as a function of $\cos(\alpha)$, in the $H \rightarrow bb$, $H\rightarrow\gamma\gamma$, $H\rightarrow ZZ$, $H\rightarrow WW$, and $H\rightarrow\tau\tau$ decay channels and in their combination. The overall modifier of the top-Higgs coupling strength is profiled such that the likelihood attains its minimum at each point in the plane, and $\kappa_{V}$ is set to 1.