Distributions of p_T(Z) × ∆m in the SR after the fit to data with a (mA , mH) signal hypothesis of (1000, 850) GeV. The pre-fit signal is arbitrarily scaled and therefore omitted.

Expected 95% CL upper limits on the production cross section times branching ratio of the A → ZH → Ztt process in the (mA , mH) plane.

Expected 95% CL upper limits on the production cross section times branching ratio of the A → ZH → Ztt process in the (mA , mH) plane.

Observed 95% CL upper limits on the production cross section times branching ratio of the A → ZH → Ztt process in the (mA , mH) plane.

Observed 95% CL upper limits on the production cross section times branching ratio of the A → ZH → Ztt process in the (mA , mH) plane.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Observed local significance of the production cross section times branching ratio of the A → ZH → Ztt process in the (mA , mH) plane.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tan β, mA) parameter space at mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tan β, mA) parameter space at mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tan β, mA) parameter space at mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tanβ, cos(β−α)) parameter space at mA = 600 GeV and mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tan β, mA) parameter space at mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tanβ, cos(β−α)) parameter space at mA = 600 GeV and mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tanβ, cos(β−α)) parameter space at mA = 600 GeV and mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tanβ, cos(β−α)) parameter space at mA = 600 GeV and mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hq}^{(3)}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hu}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hu}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hu}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hu}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $g^{(2)}_{ZZ}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $g^{(2)}_{ZZ}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

The modified energy correlation function, for data, background (pre- and post-reweighting) and three $H\rightarrow Za$ signal hypotheses (for $a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection but not the classification NN requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

$Z$ boson transverse momentum, for data, background (pre- and post-reweighting) and three $H\rightarrow Za$ signal hypotheses (for $a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection but not the classification NN requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

$Z$ boson transverse momentum, for data, background (pre- and post-reweighting) and three $H\rightarrow Za$ signal hypotheses (for $a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection but not the classification NN requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Invariant mass of the lepton pair plus jet system, for data, background (pre- and post-reweighting) and three $H\rightarrow Za$ signal hypotheses (for $a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection but not the classification NN requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Invariant mass of the lepton pair plus jet system, for data, background (pre- and post-reweighting) and three $H\rightarrow Za$ signal hypotheses (for $a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection but not the classification NN requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Regression NN output variable, for data, reweighted background, and three $H\rightarrow Za$ signal hypotheses ($a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection, including a $120<m_{\ell\ell j}<140$ GeV requirement, but not the classification NN output variable requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Regression NN output variable, for data, reweighted background, and three $H\rightarrow Za$ signal hypotheses ($a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection, including a $120<m_{\ell\ell j}<140$ GeV requirement, but not the classification NN output variable requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Classification NN output variable, for data, reweighted background, and three $H\rightarrow Za$ signal hypotheses ($a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection, including a $120<m_{\ell\ell j}<140$ GeV requirement, but not the classification NN output variable requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Classification NN output variable, for data, reweighted background, and three $H\rightarrow Za$ signal hypotheses ($a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection, including a $120<m_{\ell\ell j}<140$ GeV requirement, but not the classification NN output variable requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Invariant mass of the lepton pair plus jet system for data, reweighted background, and various signal hypotheses. Events are required to pass the complete event selection (preselection and classification NN requirement). The pre-fit background distribution is shown along with three $H\rightarrow Za$ signal hypotheses coming from the pure MC sample ($a\rightarrow q\bar{q}/gg$ inclusively). The background normalization is set equal to that of the data for events passing the preselection and being in the region 100-180 GeV. The error bars represent the data sample's statistical uncertainty and the grey bands represent the full background model uncertainty before the fit (assuming all the uncertainties are uncorrelated). The lower panel shows the ratio of the data to the pre-fit background prediction.

Invariant mass of the lepton pair plus jet system for data, reweighted background and its systematic uncertainties (only "up" variations are shown). Events are required to pass the complete event selection (preselection and classification NN requirement). The background normalization is set equal to that of the data for events passing this selection. The experimental uncertainty combines all the relevant jet uncertainties assuming they are completely uncorrelated. The error bars represent the data sample's statistical uncertainty and the grey bands represent the background statistical uncertainty. The lower panel shows the ratio of the data to the pre-fit background prediction.

Invariant mass of the lepton pair plus jet system for data, reweighted background, and 2.0 GeV signal hypotheses. Events are required to pass the complete event selection (preselection and classification NN requirement). The background shape and normalization come from the fit to the data. The signal includes only $gg$ decays, its shape comes from the fit, and it is normalized to $\mathcal{B}(H\rightarrow Za)=100\%$. The error bars represent the data sample's statistical uncertainty and the grey bands represent the full background model uncertainty after the fit. The lower panel shows the per-bin signal significance, $(N_\text{data}-N_\text{MC})\,\text{\Large{/}}\sqrt{\sigma_\text{data}^2+\sigma_\text{MC}^2}$. The post-fit background uncertainty is calculated from a fit under the background-only hypothesis. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Invariant mass of the lepton pair plus jet system for data, reweighted background, and 0.5 GeV signal hypotheses. Events are required to pass the complete event selection (preselection and classification NN requirement). The background shape and normalization come from the fit to the data. The signal includes only $gg$ decays, its shape comes from the fit, and it is normalized to $\mathcal{B}(H\rightarrow Za)=100\%$. The error bars represent the data sample's statistical uncertainty and the grey bands represent the full background model uncertainty after the fit. The lower panel shows the per-bin signal significance, $(N_\text{data}-N_\text{MC})\,\text{\Large{/}}\sqrt{\sigma_\text{data}^2+\sigma_\text{MC}^2}$. The post-fit background uncertainty is calculated from a fit under the background-only hypothesis. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Observed 95$\%$ CL upper limits on $\mathcal{B}(H \to Za)=100\%$ for $\mathcal{B}(a \to gg)=100\%$, and the expectation under the background-only hypothesis together with its $\pm1\sigma$ and $\pm2\sigma$ intervals. The weaker limits from the previous version of the analysis (PRL 125(2020) 221802) are also shown. Linear interpolation is used to set limits between the mass points for which MC signal samples were generated.

Observed 95$\%$ CL upper limits on $\mathcal{B}(H \to Za)=100\%$ for $\mathcal{B}(a \to gg)=100\%$, and the expectation under the background-only hypothesis together with its $\pm1\sigma$ and $\pm2\sigma$ intervals. The weaker limits from the previous version of the analysis (PRL 125(2020) 221802) are also shown. Linear interpolation is used to set limits between the mass points for which MC signal samples were generated.

Observed 95$\%$ confidence level upper limits on $\mathcal{B}(H \to Za)=100\%$ for $\mathcal{B}(a\rightarrow q\bar{q})=100\%$, and the expectation under the background-only hypothesis together with its $\pm1\sigma$ and $\pm2\sigma$ intervals. The weaker limits from the previous version of the analysis (PRL 125(2020) 221802) are also shown. Linear interpolation is used to set limits between the mass points for which MC signal samples were generated.

Observed 95$\%$ confidence level upper limits on $\mathcal{B}(H \to Za)=100\%$ for $\mathcal{B}(a\rightarrow q\bar{q})=100\%$, and the expectation under the background-only hypothesis together with its $\pm1\sigma$ and $\pm2\sigma$ intervals. The weaker limits from the previous version of the analysis (PRL 125(2020) 221802) are also shown. Linear interpolation is used to set limits between the mass points for which MC signal samples were generated.

The direct-photon fraction r in high-multiplicity pp collisions at $\sqrt{s}$ = 13 TeV as a function of transverse momentum for 1 < $p_{\rm T}$ < 6 GeV/$c$. r is the ratio of direct GAMMA to inclusive GAMMA.

The direct-photon inelastic cross section in pp collisions at $\sqrt{s}$ = 13 TeV as a function of transverse momentum for 1 < $p_{\rm T}$ < 6 GeV/$c$

The direct-photon invariant differential yield in high-multiplicity pp collisions at $\sqrt{s}$ = 13 TeV as a function of transverse momentum for 1 < $p_{\rm T}$ < 6 GeV/$c$

The integrated photon yield in inelastic and high-multiplicity pp collisions at $\sqrt{s}$ = 13 TeV as a function of charged particle multiplicity

The dielectron cross section in inelastic pp collisions at $\sqrt{s}$ = 13 TeV as a function of invariant mass for 1 < $p_{\rm T,ee}$ < 2 GeV/$c$.

The dielectron cross section in high-multiplicity pp collisions at $\sqrt{s}$ = 13 TeV as a function of invariant mass for 1 < $p_{\rm T,ee}$ < 2 GeV/$c$.

The dielectron cross section in inelastic pp collisions at $\sqrt{s}$ = 13 TeV as a function of invariant mass for 3 < $p_{\rm T,ee}$ < 6 GeV/$c$.

The dielectron cross section in high-multiplicity pp collisions at $\sqrt{s}$ = 13 TeV as a function of invariant mass for 3 < $p_{\rm T,ee}$ < 6 GeV/$c$.

Results for the $P_3$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the $P_4^\prime$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the $P_5^\prime$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the $P_6^\prime$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the $P_8^\prime$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the CP averaged observables $F_L$ and $P_1$–$P_8^\prime$. The first uncertainties are statistical and the second systematic.

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 1.1 < $q^2$ < 2 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 2 < $q^2$ < 4.3 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 4.3 < $q^2$ < 6 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 6 < $q^2$ < 8.68 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 10.09 < $q^2$ < 12.86 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 14.18 < $q^2$ < 16 GeV$^2$

Comparison between the values of the near-side associated peak yield for D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and $\Lambda_\mathrm{c}^{+}$ correlated with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Comparison between the values of the near-side associated peak yield for D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and $\Lambda_\mathrm{c}^{+}$ correlated with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ 0.3 < p_{\rm T} < 1 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Comparison between the values of the near-side associated peak yield for D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and $\Lambda_\mathrm{c}^{+}$ correlated with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ p_{\rm T} > 1.0 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Ratio between the values of the near-side associated peak yield of the $\Lambda_\mathrm{c}^{+}$ and D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) azimuthal correlation distribution with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Ratio between the values of the near-side associated peak yield of the $\Lambda_\mathrm{c}^{+}$ and D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) azimuthal correlation distribution with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ 0.3 < p_{\rm T} < 1 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Ratio between the values of the near-side associated peak yield of the $\Lambda_\mathrm{c}^{+}$ and D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) azimuthal correlation distribution with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ p_{\rm T} > 1.0 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Comparison between the values of the near-side associated peak width for D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and $\Lambda_\mathrm{c}^{+}$ correlated with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Comparison between the values of the near-side associated peak width for D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and $\Lambda_\mathrm{c}^{+}$ correlated with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ 0.3 < p_{\rm T} < 1.0 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Comparison between the values of the near-side associated peak width for D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and $\Lambda_\mathrm{c}^{+}$ correlated with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ p_{\rm T} > 1.0 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Ratio between the values of the near-side associated peak width of the $\Lambda_\mathrm{c}^{+}$ and D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) azimuthal correlation distribution with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Ratio between the values of the near-side associated peak width of the $\Lambda_\mathrm{c}^{+}$ and D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) azimuthal correlation distribution with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ 0.3 < p_{\rm T} < 1.0 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Ratio between the values of the near-side associated peak width of the $\Lambda_\mathrm{c}^{+}$ and D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) azimuthal correlation distribution with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ p_{\rm T} > 1.0 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Comparison between the values of the away-side associated peak yield for D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and $\Lambda_\mathrm{c}^{+}$ correlated with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Comparison between the values of the away-side associated peak yield for D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and $\Lambda_\mathrm{c}^{+}$ correlated with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ 0.3 < p_{\rm T} < 1 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Comparison between the values of the away-side associated peak yield for D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and $\Lambda_\mathrm{c}^{+}$ correlated with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ p_{\rm T} > 1.0 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Ratio between the values of the away-side associated peak yield of the $\Lambda_\mathrm{c}^{+}$ and D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) azimuthal correlation distribution with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Ratio between the values of the away-side associated peak yield of the $\Lambda_\mathrm{c}^{+}$ and D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) azimuthal correlation distribution with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ 0.3 < p_{\rm T} < 1 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

Ratio between the values of the away-side associated peak yield of the $\Lambda_\mathrm{c}^{+}$ and D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) azimuthal correlation distribution with charged particles, as a function of the charm-hadron $p_{\rm T}$, for associated particle $ p_{\rm T} > 1.0 $ GeV/$c$. The rapidity range for the charm hadrons is $|y^{\rm D,\Lambda_{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$.

$v_2${2, $1.1<|\Delta\eta|< 7.8$} of $\rm K^0_\rm S$ as a function of $p_{\mathrm{T}}$ in high-multiplicity p--Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV.

$v_2${2, $1.1<|\Delta\eta|< 7.8$} of $\Lambda$ + $\overline \Lambda$ as a function of $p_{\mathrm{T}}$ in high-multiplicity p--Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV.

$v_2${2, $1.1<|\Delta\eta|< 6.4$} of $\mathrm{\pi}^{\pm}$ as a function of $p_{\mathrm{T}}$ in high-multiplicity pp collisions at $\sqrt{s_{\mathrm{NN}}} = 13$ TeV.

$v_2${2, $1.1<|\Delta\eta|< 6.4$} of $\rm K^+$ + $\rm K^-$ as a function of $p_{\mathrm{T}}$ in high-multiplicity pp collisions at $\sqrt{s_{\mathrm{NN}}} = 13$ TeV.

$v_2${2, $1.1<|\Delta\eta|< 6.4$} of p + $\rm\overline p$ as a function of $p_{\mathrm{T}}$ in high-multiplicity pp collisions at $\sqrt{s_{\mathrm{NN}}} = 13$ TeV.

$v_2${2, $1.1<|\Delta\eta|< 6.4$} of $\rm K^0_\rm S$ as a function of $p_{\mathrm{T}}$ in high-multiplicity pp collisions at $\sqrt{s_{\mathrm{NN}}} = 13$ TeV.

$v_2${2, $1.1<|\Delta\eta|< 6.4$} of $\Lambda$ + $\overline \Lambda$ as a function of $p_{\mathrm{T}}$ in high-multiplicity pp collisions at $\sqrt{s_{\mathrm{NN}}} = 13$ TeV.

Normalized transverse momentum correlator, $\sqrt{ \langle\langle \Delta p_{{\rm T}1}\Delta p_{{\rm T}2} \rangle\rangle }/\langle\langle p_{\rm T} \rangle\rangle $, as a function of the average number of participant nucleons, $\langle{ N_{\rm part}}\rangle$, in Pb--Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV

Normalized transverse momentum correlator, $\sqrt{ \langle\langle \Delta p_{{\rm T}1}\Delta p_{{\rm T}2} \rangle\rangle }/\langle\langle p_{\rm T} \rangle\rangle $, as a function of the average number of participant nucleons, $\langle{ N_{\rm part}}\rangle$, in Xe--Xe collisions at $\sqrt{s_{\rm NN}}$ = 5.44 TeV

Normalized transverse momentum correlator, $\sqrt{ \langle\langle \Delta p_{{\rm T}1}\Delta p_{{\rm T}2} \rangle\rangle }/\langle\langle p_{\rm T} \rangle\rangle $, as a function of the average number of participant nucleons, $\langle{ N_{\rm part}}\rangle$, in Pb--Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV

Normalized transverse momentum correlator, $\sqrt{ \langle\langle \Delta p_{{\rm T}1}\Delta p_{{\rm T}2} \rangle\rangle }/\langle\langle p_{\rm T} \rangle\rangle $ as a function of the charged-particle multiplicity density for spherocity-integrated events in pp collisions at $\sqrt{s}=5.02$ TeV

Normalized transverse momentum correlator, $\sqrt{ \langle\langle \Delta p_{{\rm T}1}\Delta p_{{\rm T}2} \rangle\rangle }/\langle\langle p_{\rm T} \rangle\rangle $ as a function of the charged-particle multiplicity density for jetty events in pp collisions at $\sqrt{s}=5.02$ TeV

Normalized transverse momentum correlator, $\sqrt{ \langle\langle \Delta p_{{\rm T}1}\Delta p_{{\rm T}2} \rangle\rangle }/\langle\langle p_{\rm T} \rangle\rangle $ as a function of the charged-particle multiplicity density for isotropic events in pp collisions at $\sqrt{s}=5.02$ TeV

Normalized transverse momentum correlator, $\sqrt{ \langle\langle \Delta p_{{\rm T}1}\Delta p_{{\rm T}2} \rangle\rangle }/\langle\langle p_{\rm T} \rangle\rangle $ as a function of the charged-particle multiplicity density for spherocity-integrated events in pp collisions at $\sqrt{s}=13$ TeV

Normalized transverse momentum correlator, $\sqrt{ \langle\langle \Delta p_{{\rm T}1}\Delta p_{{\rm T}2} \rangle\rangle }/\langle\langle p_{\rm T} \rangle\rangle $ as a function of the charged-particle multiplicity density for jetty events in pp collisions at $\sqrt{s}=13$ TeV

Normalized transverse momentum correlator, $\sqrt{ \langle\langle \Delta p_{{\rm T}1}\Delta p_{{\rm T}2} \rangle\rangle }/\langle\langle p_{\rm T} \rangle\rangle $ as a function of the charged-particle multiplicity density for isotropic events in pp collisions at $\sqrt{s}=13$ TeV