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$.

Ratio of the $\eta/\pi^{0}$ ratio to the $m_{\rm T}$-scaling prediction as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV.

$x_{\rm T}$ scaled $\pi^{0}$ spectrum for pp collisions at $\sqrt{s}$ = 0.9 TeV, with $n = 5.01$.

$x_{\rm T}$ scaled $\pi^{0}$ spectrum for pp collisions at $\sqrt{s}$ = 2.76 TeV, with $n = 5.01$.

$x_{\rm T}$ scaled $\pi^{0}$ spectrum for pp collisions at $\sqrt{s}$ = 7 TeV, with $n = 5.01$.

$x_{\rm T}$ scaled $\pi^{0}$ spectrum for pp collisions at $\sqrt{s}$ = 8 TeV, with $n = 5.01$.

$x_{\rm T}$ scaled $\pi^{0}$ spectrum for pp collisions at $\sqrt{s}$ = 13 TeV, with $n = 5.01$.

$x_{\rm T}$ scaled $\eta$ spectrum for pp collisions at $\sqrt{s}$ = 2.76 TeV, with $n = 4.91$.

$x_{\rm T}$ scaled $\eta$ spectrum for pp collisions at $\sqrt{s}$ = 7 TeV, with $n = 4.91$.

$x_{\rm T}$ scaled $\eta$ spectrum for pp collisions at $\sqrt{s}$ = 8 TeV, with $n = 4.91$.

$x_{\rm T}$ scaled $\eta$ spectrum for pp collisions at $\sqrt{s}$ = 13 TeV, with $n = 4.91$.

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0-100%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 70-100%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 50-70%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 30-50%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 20-30%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 10-20%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 5-10%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 1-5%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0-1%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0.05-0.1%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0-0.1%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0.01-0.05%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0-0.01%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0-100%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 70-100%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 50-70%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 30-50%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 20-30%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 10-20%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 5-10%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 1-5%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0-1%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0.05-0.1%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0-0.1%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0.01-0.05%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0-0.01%. (Scaled for visibility in the figure)

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 70-100% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 50-70% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 30-50% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 20-30% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 10-20% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 5-10% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 1-5% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 0-1% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 0.05-0.1% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 0-0.1% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 0.01-0.05% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 0-0.01% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 70-100% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 50-70% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 30-50% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 20-30% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 10-20% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 5-10% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 1-5% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 0-1% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 0.05-0.1% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 0-0.1% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 0.01-0.05% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 0-0.01% to the inclusive spectrum in the INEL>0 event class.

Slope parameter (n) of a power-law fit ($f(p_{\rm T}) = \alpha p_{\rm T}^{-n}$) for $5 < p_{\rm T} < 15$ GeV/$c$ as a function the charged-particle multiplicity density in units of the average normalized multiplicity density for the INEL>0 event class for the neutral pion.

Slope parameter (k) of an exponential fit ($f(p_{\rm T}) = \alpha e^{-p_{\rm T}/k}$) for $0.4 < p_{\rm T} < 2$ GeV/$c$ as a function the charged-particle multiplicity density in units of the average normalized multiplicity density for the INEL>0 event class for the neutral pion.

Slope parameter (n) of a power-law fit ($f(p_{\rm T}) = \alpha p_{\rm T}^{-n}$) for $5 < p_{\rm T} < 15$ GeV/$c$ as a function the charged-particle multiplicity density in units of the average normalized multiplicity density for the INEL>0 event class for the eta meson.

Slope parameter (k) of an exponential fit ($f(p_{\rm T}) = \alpha e^{-p_{\rm T}/k}$) for $1.5 < p_{\rm T} < 4$ GeV/$c$ as a function the charged-particle multiplicity density in units of the average normalized multiplicity density for the INEL>0 event class for the eta meson.

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 0-100%.

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 70-100% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 50-70%.

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 30-50% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 20-30% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 10-20% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 5-10% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 1-5% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 0-1% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 0.05-0.1% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 0-0.1%.

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 0.01-0.05% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 0-0.01% (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 70-100% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 50-70% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV.

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 30-50% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 20-30% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 10-20% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 5-10% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 1-5% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 0-1% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 0.05-0.1% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 0-0.1% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV.

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 0.01-0.05% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 0-0.01% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Mean values of the ratio of the $\eta/\pi^{0}$ ratio in a certain multiplicity interval to the $\eta/\pi^{0}$ ratio in the INEL>0 event class for $1 < p_{\rm{ T}} < 4 \rm{~GeV}/c$ as a function of the normalized mean charged-particle multiplicity.

Mean values of the ratio of the $\eta/\pi^{0}$ ratio in a certain multiplicity interval to the $\eta/\pi^{0}$ ratio in the INEL>0 event class for $5 < p_{\rm{ T}} < 15 \rm{~GeV}/c$ as a function of the normalized mean charged-particle multiplicity.

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

$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.

Correlation functions $P_2^{\rm LS}$ of charged hadrons measured in minimum bias pp collisions at $\sqrt{s}=13\;\text{TeV}$. Charged hadrons are selected in the range $0.2 < p_{\rm T} < 2.0$ GeV/$c$ and with pseudorapidity $|\eta| < 0.8$.

Correlation functions $R_2^{\rm CI}$ of charged hadrons measured in minimum bias pp collisions at $\sqrt{s}=13\;\text{TeV}$. Charged hadrons are selected in the range $0.2 < p_{\rm T} < 2.0$ GeV/$c$ and with pseudorapidity $|\eta| < 0.8$.

Correlation functions $P_2^{\rm CI}$ of charged hadrons measured in minimum bias pp collisions at $\sqrt{s}=13\;\text{TeV}$. Charged hadrons are selected in the range $0.2 < p_{\rm T} < 2.0$ GeV/$c$ and with pseudorapidity $|\eta| < 0.8$.

Projections onto $\Delta\eta$ of $R_2^{\rm CI}$ correlation functions of charged hadrons calculated in the selected $p_{\rm T}$ range in pp collisions at $\sqrt{s}=13\;\text{TeV}$. The $\Delta\eta$ projection is calculated as averages of the two-dimensional correlations in the range $|\Delta\varphi| \leq \pi$ and $|\Delta\eta| \leq 1.6$, respectively. Simulations using PYTHIA8 with the Monash 2013 tune and color reconnection (CR) enabled are added.

Projections onto $\Delta\varphi$ of $R_2^{\rm CI}$ correlation functions of charged hadrons calculated in the selected $p_{\rm T}$ range in pp collisions at $\sqrt{s}=13\;\text{TeV}$. The $\Delta\varphi$ projection is calculated as averages of the two-dimensional correlations in the range $|\Delta\varphi| \leq \pi$ and $|\Delta\eta| \leq 1.6$, respectively. Simulations using PYTHIA8 with the Monash 2013 tune and color reconnection (CR) enabled are added.

Projections onto $\Delta\eta$ of $P_2^{\rm CI}$ correlation functions of charged hadrons calculated in the selected $p_{\rm T}$ range in pp collisions at $\sqrt{s}=13\;\text{TeV}$. The $\Delta\eta$ projection is calculated as averages of the two-dimensional correlations in the range $|\Delta\varphi| \leq \pi$ and $|\Delta\eta| \leq 1.6$, respectively. Simulations using PYTHIA8 with the Monash 2013 tune and color reconnection (CR) enabled are added.

Projections onto $\Delta\varphi$ of $P_2^{\rm CI}$ correlation functions of charged hadrons calculated in the selected $p_{\rm T}$ range in pp collisions at $\sqrt{s}=13\;\text{TeV}$. The $\Delta\varphi$ projection is calculated as averages of the two-dimensional correlations in the range $|\Delta\varphi| \leq \pi$ and $|\Delta\eta| \leq 1.6$, respectively. Simulations using PYTHIA8 with the Monash 2013 tune and color reconnection (CR) enabled are added.

Correlation functions $R_2^{\rm CD}$ of charged hadrons measured in minimum bias pp collisions at $\sqrt{s}=13\;\text{TeV}$. Charged hadrons are selected in the range $0.2 < p_{\rm T} < 2.0$ GeV/$c$ and with pseudorapidity $|\eta| < 0.8$.

Correlation functions $P_2^{\rm CD}$ of charged hadrons measured in minimum bias pp collisions at $\sqrt{s}=13\;\text{TeV}$. Charged hadrons are selected in the range $0.2 < p_{\rm T} < 2.0$ GeV/$c$ and with pseudorapidity $|\eta| < 0.8$.

Projections onto $\Delta\eta$ of $R_2^{\rm CD}$ correlation functions of charged hadrons calculated in the selected $p_{\rm T}$ range in pp collisions at $\sqrt{s}=13\;\text{TeV}$. The $\Delta\eta$ projection is calculated as averages of the two-dimensional correlations in the range $|\Delta\varphi| \leq \pi$ and $|\Delta\eta| \leq 1.6$, respectively. Simulations using PYTHIA8 with the Monash 2013 tune and color reconnection (CR) enabled are added.

Projections onto $\Delta\varphi$ of $R_2^{\rm CD}$ correlation functions of charged hadrons calculated in the selected $p_{\rm T}$ range in pp collisions at $\sqrt{s}=13\;\text{TeV}$. The $\Delta\varphi$ projection is calculated as averages of the two-dimensional correlations in the range $|\Delta\varphi| \leq \pi$ and $|\Delta\eta| \leq 1.6$, respectively. Simulations using PYTHIA8 with the Monash 2013 tune and color reconnection (CR) enabled are added.

Projections onto $\Delta\eta$ of $P_2^{\rm CD}$ correlation functions of charged hadrons calculated in the selected $p_{\rm T}$ range in pp collisions at $\sqrt{s}=13\;\text{TeV}$. The $\Delta\eta$ projection is calculated as averages of the two-dimensional correlations in the range $|\Delta\varphi| \leq \pi$ and $|\Delta\eta| \leq 1.6$, respectively. Simulations using PYTHIA8 with the Monash 2013 tune and color reconnection (CR) enabled are added.

Projections onto $\Delta\varphi$ of $P_2^{\rm CD}$ correlation functions of charged hadrons calculated in the selected $p_{\rm T}$ range in pp collisions at $\sqrt{s}=13\;\text{TeV}$. The $\Delta\varphi$ projection is calculated as averages of the two-dimensional correlations in the range $|\Delta\varphi| \leq \pi$ and $|\Delta\eta| \leq 1.6$, respectively. Simulations using PYTHIA8 with the Monash 2013 tune and color reconnection (CR) enabled are added.

Width of $R_2^{\rm CI}(\Delta\eta)$ and $P_2^{\rm CI}(\Delta\eta)$ correlation functions along $\Delta\eta$ measured within $|\Delta\varphi| < \pi$ in pp collisions and compared with results from $p-Pb$ and $Pb-Pb$ collisions as a function of $\langle \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta \rangle_{|\eta|<0.5}$.

Width of $R_2^{\rm CD}(\Delta\eta)$ and $P_2^{\rm CD}(\Delta\eta)$ correlation functions along $\Delta\eta$ measured within $|\Delta\varphi| < \pi$ in pp collisions and compared with results from p$-$Pb and Pb$-$Pb collisions as a function of $\langle \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta \rangle_{|\eta|<0.5}$.

Width of $R_2^{\rm CD}(\Delta\varphi)$ and $P_2^{\rm CD}(\Delta\varphi)$ correlation functions along $\Delta\varphi$ measured within $|\Delta\eta| < 1.6$ in pp collisions and compared with results from p$-$Pb and Pb$-$Pb collisions as a function of $\langle \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta \rangle_{|\eta|<0.5}$.

Balance function of charged hadrons, in the selected $p_{\rm T}$ range, obtained using ALICE data in pp collisions at $\sqrt{s}=13\;\text{TeV}$.

Integral of $B$ results for charged hadrons in pp collisions at $\sqrt{s}=13\;\text{TeV}$ compared with the integral of $B$ for identified hadron ($\pi \pi$, $\pi$K, and $\pi$p) pairs in Pb$-$Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$.

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the resolved qqbb signal region.

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the resolved lvbb signal region.

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the merged lvbb signal region

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the merged qqbb signal region

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 5jex3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 5jex4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 6jin3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 6jin4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 5jex3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 5jex4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 6jin3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 6jin4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 5jex3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 5jex4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 6jin3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 6jin4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

The mWh distributions in the Low-NN-score SRs of the merged qqbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Low-NN-score SRs of the merged qqbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged qqbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged qqbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Low-NN-score SRs of the merged lvbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Low-NN-score SRs of the merged lvbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Medium-NN-score SRs of the merged lvbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Medium-NN-score SRs of the merged lvbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged lvbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged lvbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

Cutflow table for a few selected signal samples after applying the selection requirements of the resolved analysis. At the initial cut stage, all events from the tbH ± → W ± h(→ b b) production process (including those from the fully hadronic decay channel) are considered.

Cutflow table for a few selected signal samples after applying the selection requirements of the merged analysis. At the initial cut stage, all events from the tbH ± → W ± h(→ bb) production process (including those from the fully hadronic decay channel) are considered.

Jet angularity $\lambda_{\alpha}$ for $\alpha = 3$. $40<p_{\mathrm{T}}^{\mathrm{ch jet}}<60$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet angularity $\lambda_{\alpha,g}$ for $\alpha = 1$. $40<p_{\mathrm{T}}^{\mathrm{ch jet}}<60$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet angularity $\lambda_{\alpha,g}$ for $\alpha = 1.5$. $40<p_{\mathrm{T}}^{\mathrm{ch jet}}<60$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet angularity $\lambda_{\alpha,g}$ for $\alpha = 2$. $40<p_{\mathrm{T}}^{\mathrm{ch jet}}<60$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet angularity $\lambda_{\alpha,g}$ for $\alpha = 3$. $40<p_{\mathrm{T}}^{\mathrm{ch jet}}<60$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet angularity $\lambda_{\alpha}$ for $\alpha = 1$. $60<p_{\mathrm{T}}^{\mathrm{ch jet}}<80$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet angularity $\lambda_{\alpha}$ for $\alpha = 1.5$. $60<p_{\mathrm{T}}^{\mathrm{ch jet}}<80$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet angularity $\lambda_{\alpha}$ for $\alpha = 2$. $60<p_{\mathrm{T}}^{\mathrm{ch jet}}<80$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet angularity $\lambda_{\alpha}$ for $\alpha = 3$. $60<p_{\mathrm{T}}^{\mathrm{ch jet}}<80$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet angularity $\lambda_{\alpha,g}$ for $\alpha = 1$. $60<p_{\mathrm{T}}^{\mathrm{ch jet}}<80$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet angularity $\lambda_{\alpha,g}$ for $\alpha = 1.5$. $60<p_{\mathrm{T}}^{\mathrm{ch jet}}<80$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet angularity $\lambda_{\alpha,g}$ for $\alpha = 2$. $60<p_{\mathrm{T}}^{\mathrm{ch jet}}<80$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet angularity $\lambda_{\alpha,g}$ for $\alpha = 3$. $60<p_{\mathrm{T}}^{\mathrm{ch jet}}<80$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet angularity $\lambda_{\alpha}$ for $\alpha = 1$. $80<p_{\mathrm{T}}^{\mathrm{ch jet}}<100$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet angularity $\lambda_{\alpha}$ for $\alpha = 1.5$. $80<p_{\mathrm{T}}^{\mathrm{ch jet}}<100$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet angularity $\lambda_{\alpha}$ for $\alpha = 2$. $80<p_{\mathrm{T}}^{\mathrm{ch jet}}<100$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet angularity $\lambda_{\alpha}$ for $\alpha = 3$. $80<p_{\mathrm{T}}^{\mathrm{ch jet}}<100$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet angularity $\lambda_{\alpha,g}$ for $\alpha = 1$. $80<p_{\mathrm{T}}^{\mathrm{ch jet}}<100$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet angularity $\lambda_{\alpha,g}$ for $\alpha = 1.5$. $80<p_{\mathrm{T}}^{\mathrm{ch jet}}<100$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet angularity $\lambda_{\alpha,g}$ for $\alpha = 2$. $80<p_{\mathrm{T}}^{\mathrm{ch jet}}<100$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet angularity $\lambda_{\alpha,g}$ for $\alpha = 3$. $80<p_{\mathrm{T}}^{\mathrm{ch jet}}<100$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet angularity $\lambda_{\alpha}$ for $\alpha = 1$. $100<p_{\mathrm{T}}^{\mathrm{ch jet}}<150$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet angularity $\lambda_{\alpha}$ for $\alpha = 1.5$. $100<p_{\mathrm{T}}^{\mathrm{ch jet}}<150$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet angularity $\lambda_{\alpha}$ for $\alpha = 2$. $100<p_{\mathrm{T}}^{\mathrm{ch jet}}<150$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet angularity $\lambda_{\alpha}$ for $\alpha = 3$. $100<p_{\mathrm{T}}^{\mathrm{ch jet}}<150$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet angularity $\lambda_{\alpha,g}$ for $\alpha = 1$. $100<p_{\mathrm{T}}^{\mathrm{ch jet}}<150$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet angularity $\lambda_{\alpha,g}$ for $\alpha = 1.5$. $100<p_{\mathrm{T}}^{\mathrm{ch jet}}<150$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet angularity $\lambda_{\alpha,g}$ for $\alpha = 2$. $100<p_{\mathrm{T}}^{\mathrm{ch jet}}<150$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet angularity $\lambda_{\alpha,g}$ for $\alpha = 3$. $100<p_{\mathrm{T}}^{\mathrm{ch jet}}<150$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet invariant mass $m_\mathrm{jet}$ in pp data. $40<p_{\mathrm{T}}^{\mathrm{ch jet}}<60$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet invariant mass $m_{\mathrm{jet},g}$ in pp data. $40<p_{\mathrm{T}}^{\mathrm{ch jet}}<60$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet invariant mass $m_\mathrm{jet}$ in pp data. $60<p_{\mathrm{T}}^{\mathrm{ch jet}}<80$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet invariant mass $m_{\mathrm{jet},g}$ in pp data. $60<p_{\mathrm{T}}^{\mathrm{ch jet}}<80$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet invariant mass $m_\mathrm{jet}$ in pp data. $80<p_{\mathrm{T}}^{\mathrm{ch jet}}<100$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet invariant mass $m_{\mathrm{jet},g}$ in pp data. $80<p_{\mathrm{T}}^{\mathrm{ch jet}}<100$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet invariant mass $m_\mathrm{jet}$ in Pb--Pb data. $40<p_{\mathrm{T}}^{\mathrm{ch jet}}<60$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet invariant mass $m_{\mathrm{jet},g}$ in Pb--Pb data. $40<p_{\mathrm{T}}^{\mathrm{ch jet}}<60$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet invariant mass $m_\mathrm{jet}$ in Pb--Pb data. $60<p_{\mathrm{T}}^{\mathrm{ch jet}}<80$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet invariant mass $m_{\mathrm{jet},g}$ in Pb--Pb data. $60<p_{\mathrm{T}}^{\mathrm{ch jet}}<80$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet invariant mass $m_\mathrm{jet}$ in Pb--Pb data. $80<p_{\mathrm{T}}^{\mathrm{ch jet}}<100$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet invariant mass $m_{\mathrm{jet},g}$ in Pb--Pb data. $80<p_{\mathrm{T}}^{\mathrm{ch jet}}<100$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet invariant mass $m_\mathrm{jet}$ in Pb--Pb data. $100<p_{\mathrm{T}}^{\mathrm{ch jet}}<150$ GeV/$c$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Groomed jet invariant mass $m_{\mathrm{jet},g}$ in Pb--Pb data. $100<p_{\mathrm{T}}^{\mathrm{ch jet}}<150$ GeV/$c$, Soft Drop $z_{\mathrm{cut}}=0.2, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding", "random_mass") no correlation information is specified ($\pm$ is always used).

Jet pairing efficiencies over the parameter space for the wide-width heavy resonance signals. The pairing efficiency is evaluated in the 6$b$ region when a correct pairing is possible — that is, the six leading jets are geometrically matched to truth-level b-quarks.

Distributions of scores for nonresDNN in 6b data and the background prediction after background-only fits to the observed data. The Low- and High-Score regions are included. Benchmark signal models (normalized to the background) are overlaid: SM HHH signal and non-resonant TRSM signal with $(m_X, m_S)=(400,200)$ GeV.

Distributions of scores for resDNN in 6b data and the background prediction after background-only fits to the observed data. The Low- and High-Score regions are included. A benchmark signal model (normalized to the background) is overlaid: resonant TRSM signal with $(m_X, m_S)=(500,350)$ GeV.

Distributions of scores for heavyresDNN in 6b data and the background prediction after background-only fits to the observed data. The Low- and High-Score regions are included. Benchmark signal models (normalized to the background) are overlaid: heavy-resonant signals with (mX, mS)=(900,325) GeV with a 1% or 20% width.

Expected 95% CL $HHH$ cross section upper limits for the phase space within the perturbative unitarity bounds of the TRSM signals. While the points were generated under Benchmark point 3 of the TRSM, they can also be interpreted in the DM-CPV model.

Observed 95% CL $HHH$ cross section upper limits for the phase space within the perturbative unitarity bounds of the TRSM signals. While the points were generated under Benchmark point 3 of the TRSM, they can also be interpreted in the DM-CPV model.

Expected 95% CL $HHH$ cross section upper limits for the narrow-width (1%) heavy resonance signals.

Expected 95% CL $HHH$ cross section upper limits for the wide-width (20%) heavy resonance signals.

Observed 95% CL $HHH$ cross section upper limits for the narrow-width (1%) heavy resonance signals.

Observed 95% CL $HHH$ cross section upper limits for the wide-width (20%) heavy resonance signals.

Expected 95% CL constraints on $\kappa_3$ and $\kappa_4$, the ratios of the Higgs tri-linear and quartic self-couplings to their predicted SM values.

Expected 68% CL constraints on $\kappa_3$ and $\kappa_4$, the ratios of the Higgs tri-linear and quartic self-couplings to their predicted SM values.

Observed 95% CL constraints on $\kappa_3$ and $\kappa_4$, the ratios of the Higgs tri-linear and quartic self-couplings to their predicted SM values.

Observed 68% CL constraints on $\kappa_3$ and $\kappa_4$, the ratios of the Higgs tri-linear and quartic self-couplings to their predicted SM values.

Unitarity limits, calculated in Eur. Phys. J. C 84 (2024) 366, are overlaid in the region bounded by the gray dashed line.

Overall acceptance x efficiency for events entering the 6$b$ region over the parameter space for the TRSM signals. Generator weights, and weights associated to pileup modeling corrections are applied. Jet trigger scale factors, NNJVT scale factors and scale factors associated to the $b$-tagging calibration are applied.

Overall acceptance x efficiency for events entering the 6$b$ region over the parameter space for the narrow-width heavy resonance signals. Generator weights, and weights associated to pileup modeling corrections are applied. Jet trigger scale factors, NNJVT scale factors and scale factors associated to the $b$-tagging calibration are applied.

Overall acceptance x efficiency for events entering the 6$b$ region over the parameter space for the wide-width heavy resonance signals. Generator weights, and weights associated to pileup modeling corrections are applied. Jet trigger scale factors, NNJVT scale factors and scale factors associated to the $b$-tagging calibration are applied.

Expected one-dimensional likelihood scan over the Higgs self-coupling modifier $\kappa_3$. Both of the modifiers are profiled in the profile-likelihood fit, and the likelihood scan keeps $\kappa_4$ fixed at one.

Observed one-dimensional likelihood scan over the Higgs self-coupling modifier $\kappa_3$. Both of the modifiers are profiled in the profile-likelihood fit, and the likelihood scan keeps $\kappa_4$ fixed at one.

Expected one-dimensional likelihood scan over the Higgs self-coupling modifier $\kappa_4$. Both of the modifiers are profiled in the profile-likelihood fit, and the likelihood scan keeps $\kappa_3$ fixed at one.

Observed one-dimensional likelihood scan over the Higgs self-coupling modifier $\kappa_4$. Both of the modifiers are profiled in the profile-likelihood fit, and the likelihood scan keeps $\kappa_3$ fixed at one.

Summary of the impact on the analysis of different sources of systematic uncertainty. The impact is estimated as the absolute value of the difference between the expected upper limit with all sources of uncertainty included and the upper limit with a subset of systematic uncertainties excluded. In the left-hand column, more heavily indented rows are a sub-category of the less indented row above.

Cutflow for selected signal points. Generator weights, and weights associated to pileup modeling corrections are applied from the second row onwards. These weight the SM events assuming the SM cross section, whilst the TRSM and Generic resonance events are weighted to 50 fb. Jet trigger scale factors are applied from the third row onwards. NNJVT scale factors are applied from the fourth row onwards. Scale factors associated to the b-tagging calibration are applied in the last three rows.