Dependence of $\overline{p}$-$\overline{d}$ correlation on pseudorapidity acceptance of $\overline{p}$ and $\overline{d}$ selection in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. Results are for 20.0--30.0$\%$ collision centrality.

Dependence of $\overline{p}$-$\overline{d}$ correlation on pseudorapidity acceptance of $\overline{p}$ and $\overline{d}$ selection in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. Results are for 50.0--60.0$\%$ collision centrality.

Fig. 6b: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Away and Toward regions after the subtraction of Number density $N_{\rm ch}$ and $\Sigma p_{\rm T}$ distributions in the transverse region for p-Pb collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. 7a: $<p_{\rm T}>$ as a function of $p_{\rm T}^{\rm trig}$ in Toward (left) and Away (right) regions after the subtraction of transverse region for pp collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. 7b: $<p_{\rm T}>$ as a function of $p_{\rm T}^{\rm trig}$ in Toward (left) and Away (right) regions after the subtraction of transverse region for p-Pb collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A1: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 0.15 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A2: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 1.0 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A3: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 0.15 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A4: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 1.0 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Differential excess $e^+e^-$ yield in 70--90\% Pb--Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV as a function of $m_{\rm ee}$ for $p_{\rm T,ee} < 0.1$ GeV/$c$. Electrons are measured within $|\eta_{\rm e}| < 0.8$ and $p_{\rm T,e} > 0.2$ GeV/$c$.

Differential $e^+e^-$ yield in 70--90\% Pb--Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV as a function of $p_{\rm T,ee}$ for 0.4 $< m_{\rm ee} < 0.7$ GeV/$c^{\rm 2}$. Electrons are measured within $|\eta_{\rm e}| < 0.8$ and $p_{\rm T,e} > 0.2$ GeV/$c$.

Differential $e^+e^-$ yield in 70--90\% Pb--Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV as a function of $p_{\rm T,ee}$ for 0.7 $< m_{\rm ee} < 1.1$ GeV/$c^{\rm 2}$. Electrons are measured within $|\eta_{\rm e}| < 0.8$ and $p_{\rm T,e} > 0.2$ GeV/$c$. The quoted upper limits correspond to a 90% confidence level.

Differential $e^+e^-$ yield in 70--90\% Pb--Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV as a function of $p_{\rm T,ee}$ for 1.1 $< m_{\rm ee} < 2.7$ GeV/$c^{\rm 2}$. Electrons are measured within $|\eta_{\rm e}| < 0.8$ and $p_{\rm T,e} > 0.2$ GeV/$c$. The quoted upper limits correspond to a 90% confidence level.

Differential excess $e^+e^-$ yield in 70--90\% Pb--Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV as a function of $p^{2}_{\rm T,ee}$ for 0.4 $< m_{\rm ee} < 0.7$ GeV/$c^{\rm 2}$. Electrons are measured within $|\eta_{\rm e}| < 0.8$ and $p_{\rm T,e} > 0.2$ GeV/$c$. The quoted upper limits correspond to a 90% confidence level.

Differential excess $e^+e^-$ yield in 70--90\% Pb--Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV as a function of $p^{2}_{\rm T,ee}$ for 0.7 $< m_{\rm ee} < 1.1$ GeV/$c^{\rm 2}$. Electrons are measured within $|\eta_{\rm e}| < 0.8$ and $p_{\rm T,e} > 0.2$ GeV/$c$. The quoted upper limits correspond to a 90% confidence level.

Differential excess $e^+e^-$ yield in 70--90\% Pb--Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV as a function of $p^{2}_{\rm T,ee}$ for 1.1 $< m_{\rm ee} < 2.7$ GeV/$c^{\rm 2}$. Electrons are measured within $|\eta_{\rm e}| < 0.8$ and $p_{\rm T,e} > 0.2$ GeV/$c$.

$p_{\rm T}$-differential yields of $\Lambda$(1520) (sum of particle and anti-particle states) in p--Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ $\mathrm{=}$ 5.02 TeV in multiplicity interval 20--40\%. The uncertainty 'sys,$p_{\rm T}$-correlated' indicates the systematic uncertainty after removing the contributions of $p_{\rm T}$-uncorrelated uncertainty.

$p_{\rm T}$-differential yields of $\Lambda$(1520) (sum of particle and anti-particle states) in p--Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ $\mathrm{=}$ 5.02 TeV in multiplicity interval 40--60\%. The uncertainty 'sys,$p_{\rm T}$-correlated' indicates the systematic uncertainty after removing the contributions of $p_{\rm T}$-uncorrelated uncertainty.

$p_{\rm T}$-differential yields of $\Lambda$(1520) (sum of particle and anti-particle states) in p--Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ $\mathrm{=}$ 5.02 TeV in multiplicity interval 60--100\%. The uncertainty 'sys,$p_{\rm T}$-correlated' indicates the systematic uncertainty after removing the contributions of $p_{\rm T}$-uncorrelated uncertainty.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to charged $\pi$ ($\pi^{+} + \pi^{-}$) at midrapidity as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in inelastic pp collisions at $\sqrt{s}$ $\mathrm{=}$ 7 TeV.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to charged $\pi$ ($\pi^{+} + \pi^{-}$) as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in p--Pb collisions at $\sqrt{s}$ $\mathrm{=}$ 5.02 TeV. The uncertainty 'sys,uncorrelated' indicates the multiplicity uncorrelated systematic uncertainty.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to K$^{-}$ at midrapidity as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in inelastic pp collisions at $\sqrt{s}$ $\mathrm{=}$ 7 TeV.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to K$^{-}$ as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in p--Pb collisions at $\sqrt{s}$ $\mathrm{=}$ 5.02 TeV. The uncertainty 'sys,uncorrelated' indicates the multiplicity uncorrelated systematic uncertainty.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to its ground state particle $\Lambda$ ($\Lambda + \bar{\Lambda}$) at midrapidity as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in inelastic pp collisions at $\sqrt{s}$ $\mathrm{=}$ 7 TeV.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to its ground state particle $\Lambda$ ($\Lambda + \bar{\Lambda}$) as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in p--Pb collisions at $\sqrt{s}$ $\mathrm{=}$ 5.02 TeV. The uncertainty 'sys,uncorrelated' indicates the multiplicity uncorrelated systematic uncertainty.

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.

Correlation matrix for the individual sources of systematic uncertainty considered in the analysis, for the combined opposite sign (OS) and same sign (SS) binned-template profile likelihood fit to data. The OS or SS refers to the charge signs of the primary lepton and the soft muon. The gamma parameters are NPs used to describe the effect of the limited statistics of the samples.

Best fit values of $\sigma(gg\to H\to ZZ)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{ZZ}$ relative to their SM prediction from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The results are shown for the combination of ATLAS and CMS, and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The expected uncertainties in the measurements are also shown.

Best fit values of $\sigma(gg\to H\to WW)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{WW}$ from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The values involving cross sections are given for $\sqrt{s}$=8 TeV, assuming the SM values for $\sigma_i(7 \mathrm{TeV})/\sigma_i(8 \mathrm{TeV})$. The results are shown for the combination of ATLAS and CMS, and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The expected uncertainties in the measurements are also shown.

Best fit values of $\sigma(gg\to H\to WW)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{WW}$ relative to their SM prediction from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The results are shown for the combination of ATLAS and CMS, and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The expected uncertainties in the measurements are also shown.

Best fit values of $\kappa_{gZ}= \kappa_{g}\cdot\kappa_{Z} / \kappa_{H}$ and of the ratios of coupling modifiers, as defined in the most generic parameterisation described in the context of the $\kappa$ framework, from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The results are shown for the combination of ATLAS and CMS and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The uncertainties in $\lambda_{tg}$ and $\lambda_{WZ}$, for which a negative solution is allowed, are calculated around the overall best fit value. The expected total uncertainties in the measurements are also shown. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Measured global signal strength $\mu$ and its total uncertainty, together with the breakdown of the uncertainty into its four components as defined in the text. The results are shown for the combination of ATLAS and CMS, and separately for each experiment. The expected uncertainty, with its breakdown, is also shown.

Measured signal strengths $\mu$ and their total uncertainties for different Higgs boson production processes. The results are shown for the combination of ATLAS and CMS, and separately for each experiment, for the combined $\sqrt{s}$=7 and 8 TeV data. The expected uncertainties in the measurements are also displayed. These results are obtained assuming that the Higgs boson branching fractions are the same as in the SM.

Measured signal strengths $\mu$ and their total uncertainties for different Higgs boson decay channels. The results are shown for the combination of ATLAS and CMS, and separately for each experiment, for the combined $\sqrt{s}$=7 and 8 TeV data. The expected uncertainties in the measurements are also displayed. These results are obtained assuming that the Higgs boson production process cross sections at $\sqrt{s}$ = 7 and 8 TeV are the same as in the SM.

Results of the ten-parameter fit of $\mu_F^f = \mu_{gg\mathrm{F}+ttH}^f$ and $\mu_V^f = \mu_{\mathrm{VBF}+VH}^f$ for each of the five decay channels. The results are shown for the combination of ATLAS and CMS, together with their measured and expected uncertainties. The measured results are also shown separately for each experiment.

Results of the six-parameter fit of the global ratio $\mu_V/\mu_F = \mu_{\mathrm{VBF}+VH}/\mu_{gg\mathrm{F}+ttH}$ together with $\mu_F^f$ for each of the five decay channels. The results are shown for the combination of ATLAS and CMS, together with their measured and expected uncertainties. The measured results are also shown separately for each experiment.

Fit results allowing BSM loop couplings, assuming that $|\kappa_{V}| \le$ 1, where $\kappa_{V}$ denotes $\kappa_{Z}$ or $\kappa_{W}$, and that $\mathrm{B_{BSM}} \ge$ 0. The results for the combination of ATLAS and CMS are reported with their measured and expected uncertainties. Also shown are the results from each experiment. For the parameters with both signs allowed, the other 1$\sigma$ interval is shown on a second line. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Fit results allowing BSM loop couplings, assuming that there are no additional BSM contributions to the Higgs boson width, i.e. $\mathrm{B_{BSM}}=0$. The results for the combination of ATLAS and CMS are reported with their measured and expected uncertainties. Also shown are the results from each experiment. For the parameters with both signs allowed, the other 1$\sigma$ interval is shown on a second line. When a parameter is constrained and reaches a boundary, namely $\mathrm{B_{BSM}}$ = 0, the uncertainty is not indicated. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Fit results for the parameterisation assuming the absence of BSM particles in the loops ($\mathrm{B_{BSM}}$=0). The results with their measured and expected uncertainties are reported for the combination of ATLAS and CMS, together with the individual results from each experiment. For the parameters with both signs allowed, the other 1$\sigma$ CL interval is shown on a second line. When a parameter is constrained and reaches a boundary, namely $|\kappa_\mu|$ = 0, the uncertainty is not indicated. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Summary of fit results for a parameterisation probing the ratios of coupling modifiers for up-type versus down-type fermions. The results for the combination of ATLAS and CMS are reported together with their measured and expected uncertainties. Also shown are the results from each experiment. The parameter $\kappa_{uu}$ is positive definite since $\kappa_H$ is always assumed to be positive. For the parameter $\lambda_{du}$, for which both signs are allowed, the other 1$\sigma$ CL interval is shown on a second line. Negative values for the parameter $\lambda_{Vu}$ are excluded by more than 4$\sigma$.

Summary of fit results for a parameterisation probing the ratios of coupling modifiers for leptons versus quarks. The results for the combination of ATLAS and CMS are reported together with their measured and expected uncertainties. Also shown are the results from each experiment. The parameter $\kappa_{qq}$ is positive definite since $\kappa_H$ is always assumed to be positive. For the parameter $\lambda_{lq}$, for which there is no sensitivity to the sign, only the absolute values are shown. Negative values for the parameter $\lambda_{Vq}$ are excluded by more than 4$\sigma$.

Correlation matrix obtained from the fit combining the ATLAS and CMS data using the generic parameterisation with 23 parameters, $\sigma_i\cdot\mathrm{B}^f$ for each specific channel $i\to H\to f$. Only 20 parameters are shown because the other three, corresponding to the $H\to ZZ$ decay channel for the $WH$, $ZH$, and $ttH$ production processes, are not measured with a meaningful precision.

Correlation matrix obtained from the fit combining the ATLAS and CMS data using the generic parameterisation with nine parameters, $\sigma(gg\to H\to ZZ)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{ZZ}$.

Correlation matrix obtained from the fit combining the ATLAS and CMS data using the generic parameterisation with seven parameters, $\kappa_{gZ}=\kappa_{g}\cdot\kappa_{Z} / \kappa_{H}$ and the ratios of coupling modifiers.

Expected correlation matrix obtained from the fit combining ATLAS and CMS pre-fit Asimov data sets using the generic parameterisation with 23 parameters, $\sigma_i\cdot\mathrm{B}^f$ for each specific channel $i\to H\to f$. Only 20 parameters are shown because the other three, corresponding to the $H\to ZZ$ decay channel for the $WH$, $ZH$, and $ttH$ production processes, are not measured with a meaningful precision.

Expected correlation matrix obtained from the fit combining ATLAS and CMS pre-fit Asimov data sets using the generic parameterisation with nine parameters, $\sigma(gg\to H\to ZZ)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{ZZ}$.

Expected correlation matrix obtained from the fit combining ATLAS and CMS pre-fit Asimov data sets using the generic parameterisation with seven parameters, $\kappa_{gZ}=\kappa_{g}\cdot\kappa_{Z} / \kappa_{H}$ and the ratios of coupling modifiers.

ATLAS+CMS combined best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{g}$,$\kappa_{\gamma}$) plane.

ATLAS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{g}$,$\kappa_{\gamma}$) plane.

CMS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{g}$,$\kappa_{\gamma}$) plane.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{F}$,$\kappa_{V}$) plane.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\gamma\gamma$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow ZZ$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow WW$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\tau\tau$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow bb$ channel.

ATLAS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{F}$,$\kappa_{V}$) plane.

CMS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{F}$,$\kappa_{V}$) plane.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\gamma\gamma$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow ZZ$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow WW$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\tau\tau$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow bb$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow\gamma\gamma$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow ZZ$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow WW$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow\tau\tau$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow bb$ channel.

The +1 $\sigma$ band of the observed exclusion contor for the $Z^{\prime}_{B}$ model in the $m_{\chi}$-$m_{Z^{\prime}_{B}}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the $Z^{\prime}_{B}$ model in the $m_{\chi}$-$m_{Z^{\prime}_{B}}$ plane.

A comparison of the inferred limits to the constraints from direct detection experiments on the spin-independent DM--nucleon cross section in the context of the $Z'_B$ simplified model with vector couplings. Limits are shown at 90% CL.

The observed exclusion contor for the $Z^{\prime}$-2HDM model in the $m_{A}$-$m_{Z^{\prime}}$ plane.

The expected exclusion contor for the $Z^{\prime}$-2HDM model in the $m_{A}$-$m_{Z^{\prime}}$ plane.

The +1 $\sigma$ band of the observed exclusion contor for the $m_{A}$-$Z^{\prime}$-2HDM model in the $m_{Z^{\prime}}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the $m_{A}$-$Z^{\prime}$-2HDM model in the $m_{Z^{\prime}}$ plane.

The observed exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The expected exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The +1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The observed exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The expected exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The +1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The exclusion limits at 95% CL for the 2HDM+a model as a function of $\sin \theta$ for $m_{A,H^{\pm},H}$= 600GeV, $m_a$ = 200GeV, $\tan \beta$ = 1.0$.

The exclusion limits at 95% CL for the 2HDM+a model as a function of $\sin \theta$ for $m_{A,H^{\pm},H}$= 1000GeV, $m_a$ = 350GeV, $\tan \beta$ = 1.0$.

Breakdown of the dominant systematic uncertainties.

Event yields in the range of 120 $<m_{\gamma\gamma}<$ 130 GeV for data, signal models, the SM Higgs boson background and non-resonant background in each analysis category, for an integrated luminosity of $139$fb$^{-1}$.

Detailed background contributions from the SM Higgs boson and continuum background for each cut

Detailed contributions from the signals for each cut.

Acceptance times efficiency for several signals in each category.

Fitted event-level b-tagging selection efficiency as a function of the dijet invariant mass. The uncertainties are evaluated by propagating the uncertainties on the fit parameters to the fitted efficiency.