Parton dijet $M_{inv}$ vs $A_{LL}$ values with associated uncertainties, for topology C.

Parton dijet $M_{inv}$ vs $A_{LL}$ values with associated uncertainties, for topology D.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements.The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements coupling with the forward-forward dijet measurements (topology A). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements coupling with the forward-forward dijet measurements (topology B). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements coupling with the forward-forward dijet measurements (topology C). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements coupling with the forward-forward dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) for forward-forward dijet measurements (topology A). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology A) with forward-central dijet measurements (topology B). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology A) with forward-central dijet measurements (topology C). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology A) with forward-central dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) for forward-forward dijet measurements (topology B). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology B) with forward-central dijet measurements (topology C). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology B) with forward-central dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) for forward-forward dijet measurements (topology C). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology C) with forward-central dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) for forward-forward dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} =$ 5 TeV and 13 TeV, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} =$ 5 TeV and 13 TeV, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} =$ 5 TeV and 13 TeV, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the near-side associated peak yields measured in pp collisions at $\sqrt{s} = $5, 7, and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak yields measured in pp collisions at $\sqrt{s}=$ 5 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width measured in pp collisions at $\sqrt{s} = $5, 7, and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width measured in pp collisions at $\sqrt{s}=$ 5 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak yields measured in pp collisions at $\sqrt{s}=$ 5 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak yields measured in pp collisions at $\sqrt{s}=$ 5 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak width measured in pp collisions at $\sqrt{s}=$ 5 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak width measured in pp collisions at $\sqrt{s}=$ 5 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the azimuthal-correlation distributions of of D$^{0}$ mesons with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $p_{\rm T} > 0.3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $p_{\rm T} > 0.3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $p_{\rm T} > 0.3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $p_{\rm T} > 0.3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $2 < p_{\rm T} < 3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $2 < p_{\rm T} < 3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $2 < p_{\rm T} < 3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $2 < p_{\rm T} < 3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the near-side associated peak yields among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $p_{\rm T} > 0.3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak yields among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak yields among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak yields among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak yields among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak width among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

$p_\mathrm{T}$-differential production yield of prompt $\mathrm{D_s^+}$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$. Branching ratio of $\mathrm{D_s^+\to \phi\pi^+\to K^-K^+\pi^+}$: 0.0224

$p_\mathrm{T}$-differential nuclear modification factor $R_\mathrm{AA}$ of prompt $\mathrm{D_s^+}$ mesons in the 0-10% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

$p_\mathrm{T}$-differential nuclear modification factor $R_\mathrm{AA}$ of prompt $\mathrm{D_s^+}$ mesons in the 0-10% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

$p_\mathrm{T}$-differential nuclear modification factor $R_\mathrm{AA}$ of prompt $\mathrm{D_s^+}$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

$p_\mathrm{T}$-differential nuclear modification factor $R_\mathrm{AA}$ of prompt $\mathrm{D_s^+}$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

$p_\mathrm{{T}}$-differential production ratios of prompt $\mathrm{D}_s^+/\mathrm{D}^0$ mesons in the 0-10% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$. Branching ratio of $\mathrm{D_s^+\to \phi\pi^+\to K^-K^+\pi^+}$: 0.0224 Branching ratio of $\mathrm{D^0 \to K^-\pi^+}$: 0.0395

$p_\mathrm{{T}}$-differential production ratios of prompt $\mathrm{D}_s^+/\mathrm{D}^0$ mesons in the 0-10% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$. Branching ratio of $\mathrm{D_s^+\to \phi\pi^+\to K^-K^+\pi^+}$: 0.0224 Branching ratio of $\mathrm{D^0 \to K^-\pi^+}$: 0.0395

$p_\mathrm{{T}}$-differential production ratios of prompt $\mathrm{D}_s^+/\mathrm{D}^0$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$. Branching ratio of $\mathrm{D_s^+\to \phi\pi^+\to K^-K^+\pi^+}$: 0.0224 Branching ratio of $\mathrm{D^0 \to K^-\pi^+}$: 0.0395

$p_\mathrm{{T}}$-differential production ratios of prompt $\mathrm{D}_s^+/\mathrm{D}^0$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$. Branching ratio of $\mathrm{D_s^+\to \phi\pi^+\to K^-K^+\pi^+}$: 0.0224 Branching ratio of $\mathrm{D^0 \to K^-\pi^+}$: 0.0395

$p_\mathrm{{T}}$-differential prompt $\mathrm{D}_s^+/\mathrm{D}^0$ ratios in Pb-Pb collisions in the 0-10% centrality class at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$ divided by those in pp collisions at the same energy.

$p_\mathrm{{T}}$-differential prompt $\mathrm{D}_s^+/\mathrm{D}^0$ ratios in Pb-Pb collisions in the 0-10% centrality class at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$ divided by those in pp collisions at the same energy.

$p_\mathrm{{T}}$-differential prompt $\mathrm{D}_s^+/\mathrm{D}^0$ ratios in Pb-Pb collisions in the 30-50% centrality class at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$ divided by those in pp collisions at the same energy.

$p_\mathrm{{T}}$-differential prompt $\mathrm{D}_s^+/\mathrm{D}^0$ ratios in Pb-Pb collisions in the 30-50% centrality class at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$ divided by those in pp collisions at the same energy.

$p_\mathrm{T}$-differential elliptic flow coefficient $v_2$ of prompt $\mathrm{D_s^+}$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$.

$p_\mathrm{T}$-differential elliptic flow coefficient $v_2$ of prompt $\mathrm{D_s^+}$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$.

Transverse momentum distributions d$N$/(d$y$d$p_{\mathrm{T}}$) of prompt D$^{+}$ (middle) meson in the 30-50% (diamonds) centrality classes in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in |y|<0.5. The uncertainties on the BRs are quoted separately and the horizontal bars represent bin widths. The points in the figure of the paper are scaled for visibility by a factor 10$^{-1}$. Branching ratio of D$^+$ --> K$\pi\pi$ = 0.0938

Transverse momentum distributions d$N$/(d$y$d$p_{\mathrm{T}}$) of prompt D*$^{+}$ (right) meson in the 0-10% (crosses) centrality classes in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in |y|<0.5. The uncertainties on the BRs are quoted separately and the horizontal bars represent bin widths. Branching ratio of D*$^+$ --> D$^0\pi$ = 0.677

Transverse momentum distributions d$N$/(d$y$d$p_{\mathrm{T}}$) of prompt D*$^{+}$ (right) meson in the 30-50% (diamonds) centrality classes in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in |y|<0.5. The uncertainties on the BRs are quoted separately and the horizontal bars represent bin widths. The points in the figure of the paper are scaled for visibility by a factor 10$^{-1}$. Branching ratio of D*$^+$ --> D$^0\pi$ = 0.677

Ratio D$^{+}$/D$^{0}$ in central Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in |y|<0.5 and as a function of $p_\mathrm{T}$. Statistical (bars) and systematic (boxes) uncertainties are shown. Branching ratio of D$^0$ --> K$\pi$ = 0.0395 Branching ratio of D$^+$ --> K$\pi\pi$ = 0.0938

Ratio D$^{+}$/D$^{0}$ in semicentral Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in |y|<0.5 and as a function of $p_\mathrm{T}$. Statistical (bars) and systematic (boxes) uncertainties are shown. Branching ratio of D$^0$ --> K$\pi$ = 0.0395 Branching ratio of D$^+$ --> K$\pi\pi$ = 0.0938

Ratio D*$^{+}$/D$^{0}$ in central Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in |y|<0.5 and as a function of $p_\mathrm{T}$. Statistical (bars) and systematic (boxes) uncertainties are shown. Branching ratio of D$^0$ --> K$\pi$ = 0.0395 Branching ratio of D*$^+$ --> D$^0\pi$ = 0.677

Ratio D*$^{+}$/D$^{0}$ in semicentral Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in |y|<0.5 and as a function of $p_\mathrm{T}$. Statistical (bars) and systematic (boxes) uncertainties are shown. Branching ratio of D$^0$ --> K$\pi$ = 0.0395 Branching ratio of D*$^+$ --> D$^0\pi$ = 0.677

Nuclear modification factor of D$^{0}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 0-10% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Nuclear modification factor of D$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 0-10% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Nuclear modification factor of D*$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 0-10% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Nuclear modification factor of D$^{0}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 30-50% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Nuclear modification factor of D$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 30-50% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Nuclear modification factor of D*$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 30-50% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Average nuclear modification factor of D$^{0}$, D$^{+}$ and D*$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 0-10% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Average nuclear modification factor of D$^{0}$, D$^{+}$ and D*$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 30-50% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

$p_\mathrm{T}$-integrated nuclear modification factor of D$^{0}$ meson in central (0-10%) and semicentral (30-50%) Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in |y|<0.5.

Forward neutron single spin asymmetries in p+p, p+Al and p+Au collisions as a function of pT in bins of xF correlated with hard collision activity

Forward neutron single spin asymmetries in p+p, p+Al and p+Au collisions as a function of pT in bins of xF without hard collision activity

Forward neutron single spin asymmetries in p+p collisions as a function of xF in bins of pT (anti-)correlated with hard collision activity

Forward neutron single spin asymmetries in p+Al collisions as a function of xF in bins of pT (anti-)correlated with hard collision activity

Forward neutron single spin asymmetries in p+Au collisions as a function of xF in bins of pT (anti-)correlated with hard collision activity

Fig. 11 right: The $p_\mathrm{T}$ differential production cross section of charged-particle anti-$k_{T}$ $R=0.4$ b jets in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}~=~5.02$ TeV using the IP method.

Fig. 12 left: The combined differential production cross section of charged-particle anti-$k_{T}$ $R=0.4$ b jets in pp collisions at $\sqrt{s}~=~5.02$ TeV.

Fig. 12 right: The combined differential production cross section of charged-particle anti-$k_{T}$ $R=0.4$ b jets in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}~=~5.02$ TeV.

Fig. 13 left: The b-jet fraction in pp collisions at $\sqrt{s}~=~5.02$ TeV.

Fig. 13 right: The b-jet fraction in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}~=~5.02$ TeV.

Fig. 14 left: The nuclear modification factor $R_{\mathrm{pPb}}$ of b jets as a function of $P_T$ from the SV method.

Fig. 14 left: The nuclear modification factor $R_{\mathrm{pPb}}$ of b jets as a function of $P_T$ from the IP method.

Fig. 14 right: The combined nuclear modification factor $R_{\mathrm{pPb}}$ of b jets as a function of $P_T$.

Covariances $c_{ij}$ for the highest ranked systematic uncertainties in the energy asymmetry measurement differential in $\theta_j$.

Definition of the fiducial phase space with lepton candidate ($\ell$), hadronic top candidate ($j_h$), leptonic top candidate ($j_l$) and associated jet candidate ($j_a$).

Bounds on individual Wilson coefficients $C($TeV$/\Lambda)^2$ from the energy asymmetry, obtained from a combined fit to its values in all three $\theta_j$ bins.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1400-1500 MeV.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1500-1600 MeV.

$M_{p\eta}$ invariant mass distribution under the conditions: $E_\gamma$ = 1400-1500 MeV; 1120 MeV $\leq M_{p\pi^0} \leq$ 1220 MeV; $\theta_p^{lab} \leq$ 25°; 25° $\leq \theta_\gamma^{lab} \leq$ 155°.

$M_{p\eta}$ invariant mass distribution under the conditions: $E_\gamma$ = 1400-1500 MeV; 1120 MeV $\leq M_{p\pi^0} \leq$ 1220 MeV; $\theta_p^{lab} \leq$ 25°; 1° $\leq \theta_\gamma^{lab} \leq$ 155°.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1400-1500 MeV and $M_{p\pi^0} \leq$ 1190 MeV, allowing for all proton and photon laboratory angles.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1170 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1180 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1200 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1210 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1220 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1400-1420 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1460 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1460-1500 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1500-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1540-1600 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{p\eta}$.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{\pi^0\eta}$.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{p\pi^0}$, with a cut on the structure in the $M^2_{p\eta}$ invariant mass distribution.

Excess yield relative to the PWA distributions in the $p-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $p-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $\pi^0-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $\pi^0-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Differential cross section $d\sigma/dM_{\pi^0\eta}$ for $E_\gamma$ = 1540-1600 MeV.

$M_{p\eta}$ invariant mass distribution, taken as difference to the PWA calculations (s. Fig. 11).

$M_{\pi^0\eta}$ invariant mass distribution, taken as difference to the PWA calculations (s. Fig. 11).

$\theta_{p\eta}$ opening angle distribution, taken as difference to the PWA calculations (s. Fig. 12).

$\theta_{\pi^0\eta}$ opening angle distribution, taken as difference to the PWA calculations (s. Fig. 12).

Observed likelihood scan in the ($\kappa_\tau$, $\tilde{\kappa}_\tau$) plane.

Observed likelihood scan in the ($\alpha^{\mathrm{H}\tau\tau}$, $\mu_{\text{ggH}}$) plane. The $\mu_{\text{qqH}}$ signal strength modifier is profiled.

Observed likelihood scan in the ($\alpha^{\mathrm{H}\tau\tau}$, $\mu_{\text{qqH}}$) plane. The $\mu_{\text{ggH}}$ signal strength modifier is profiled.

Observed likelihood scan in the ($\mu_{\text{qqH}}$, $\mu_{\text{ggH}}$) plane. The value of $\alpha^{\mathrm{H}\tau\tau}$ is profiled.

Two-dimensional distributions of the leading electron vs the leading muon $|d_{0}|$, for the events that pass the e$\mu$ preselection. Data events and $\tilde{t} \to b\ell$ signal events with a $\tilde{t}$ mass of 700 GeV and a proper decay length of 10 mm that correspond to 2018 data-taking conditions are shown. In each $|d_{0}|$-$|d_{0}|$ bin, the number of events divided by the bin area is plotted. The inclusive signal region covers the region between 100 $\mu$m and 10 cm in each $|d_{0}|$ variable shown.

Two-dimensional distributions of the leading vs the subleading electron $|d_0|$, for the events that pass the ee preselection. Data events and $\tilde{t} \to b\ell$ signal events with a $\tilde{t}$ mass of 700 GeV and a proper decay length of 10 mm that correspond to 2018 data-taking conditions are shown. In each $|d_{0}|$-$|d_{0}|$ bin, the number of events divided by the bin area is plotted. The inclusive signal region covers the region between 100 $\mu$m and 10 cm in each $|d_{0}|$ variable shown.

Two-dimensional distributions of the leading vs the subleading muon $|d_{0}|$, for the events that pass the $\mu\mu$ preselection. Data events and $\tilde{t} \to b\ell$ signal events with a $\tilde{t}$ mass of 700 GeV and a proper decay length of 10 mm that correspond to 2018 data-taking conditions are shown. In each $|d_{0}|$-$|d_{0}|$ bin, the number of events divided by the bin area is plotted. The inclusive signal region covers the region between 100 $\mu$m and 10 cm in each $|d_{0}|$ variable shown.

The number of observed and estimated background events in each channel and SR, with a representative signal overlaid. For each background estimate and signal yield, the total uncertainty (statistical plus systematic) is given. The distributions shown correspond to the events before the final maximum likelihood fit to the data.

The number of observed and estimated background events in each channel and SR in 2016, with a representative signal overlaid. For each background estimate and signal yield, the total uncertainty is given. The distributions shown correspond to the events before the final maximum likelihood fit to the data.

The number of observed and estimated background events in each channel and SR in 2017 and 2018, with a representative signal overlaid. For each background estimate and signal yield, the total uncertainty is given. The distributions shown correspond to the events before the final maximum likelihood fit to the data.

The observed 95% CL upper limits on the long-lived top squark production cross section, in the $c\tau_0$ -- mass plane, for the three channels combined. The $\tilde{t} \to b\ell $ process is shown. These limits assume $\mathcal{B}(\tilde{t} \to b\ell)$ is 100%, and each lepton has an equal probability of being an electron, a muon, or a tau lepton. The area to the left of the black curve represents the observed exclusion region, and the dashed red lines indicate the expected limits and their 68% confidence intervals.

The observed 95% CL upper limits on the long-lived top squark production cross section, in the $c\tau_0$ -- mass plane, for the three channels combined. The $\tilde{t} \to d\ell$ process is shown. These limits assume $\mathcal{B}(\tilde{t} \to d\ell)$ is 100%, and each lepton has an equal probability of being an electron, a muon, or a tau lepton. The area to the left of the black curve represents the observed exclusion region, and the dashed red lines indicate the expected limits and their 68% confidence intervals.

The observed 95% CL upper limits on the long-lived top squark production cross section, in the $c\tau_0$ -- mass plane, for the three channels combined. The $\tilde{t} \to b\ell $ process is shown. These limits assume $\mathcal{B}(\tilde{t} \to b\ell)$ is 100%, and each lepton has an equal probability of being an electron, a muon, or a tau lepton. The area to the left of the black curve represents the observed exclusion region, and the dashed red lines indicate the expected limits and their 68% confidence intervals.

The observed 95% CL upper limits on the long-lived top squark production cross section, in the $c\tau_0$ -- mass plane, for the three channels combined. The $\tilde{t} \to d\ell$ process is shown. These limits assume $\mathcal{B}(\tilde{t} \to d\ell)$ is 100%, and each lepton has an equal probability of being an electron, a muon, or a tau lepton. The area to the left of the black curve represents the observed exclusion region, and the dashed red lines indicate the expected limits and their 68% confidence intervals.

The 95% CL upper limits on the long-lived top squark proper decay length ($c\tau_0$) as a function of its mass, for the e$\mu$, ee, and $\mu\mu$ channels, and their combination. The $\tilde{t} \to b\ell $ process is shown. These limits assume $\mathcal{B}(\tilde{t} \to b\ell)$ is 100%, and each lepton has an equal probability of being an electron, a muon, or a tau lepton.

The 95% CL upper limits on the long-lived top squark proper decay length ($c\tau_0$) as a function of its mass, for the e$\mu$, ee, and $\mu\mu$ channels, and their combination. The $\tilde{t} \to d\ell$ process is shown. These limits assume $\mathcal{B}(\tilde{t} \to d\ell)$ is 100%, and each lepton has an equal probability of being an electron, a muon, or a tau lepton.

The 95% CL constraints on the long-lived slepton $c\tau_{0}$ and mass. The $\tilde{\tau}$ and co-NLSP limits are shown for the three channels combined, while the $\tilde{e}$ and $\tilde{\mu}$ NLSP limits are shown for the ee and $\mu\mu$ channels, respectively. These limits assume that the superpartners of the left- and right-handed leptons are degenerate in mass and $\mathcal{B}(\tilde{\ell} \to \ell\tilde{G})$ is 100%. The area to the left of the solid curves represents the observed exclusion region, and the dashed lines indicate the expected limits.

The observed 95% CL upper limits on the long-lived slepton production cross section, in the $c\tau_0$ -- mass plane. The co-NLSP limits are shown for the three channels combined. These limits assume that the superpartners of the left- and right-handed leptons are degenerate in mass and $\mathcal{B}(\tilde{\ell} \to \ell\tilde{G})$ is 100%. The area to the left of the black curve represents the observed exclusion region, and the dashed red lines indicate the expected limits and their 68% confidence intervals.

The observed 95% CL upper limits on the long-lived slepton production cross section, in the $c\tau_0$ -- mass plane. The co-NLSP limits are shown for the three channels combined. These limits assume that the superpartners of the left- and right-handed leptons are degenerate in mass and $\mathcal{B}(\tilde{\ell} \to \ell\tilde{G})$ is 100%. The area to the left of the black curve represents the observed exclusion region, and the dashed red lines indicate the expected limits and their 68% confidence intervals.

The 95% CL upper limits on the $\mathrm{H} \to \mathrm{S}\mathrm{S}, \mathrm{S} \to \ell^{+}\ell^{-}$ branching fraction as a function of $c\tau_0$, for a Higgs boson with a mass of 125 GeV and a long-lived scalar with a mass of 30 GeV or 50 GeV, for the three channels combined. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons or two muons. The area above the solid (dashed) curve represents the observed (expected) exclusion region.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S}, \mathrm{S} \to \ell^{+}\ell^{-}$ as a function of $c\tau_0$, for a 300 GeV Higgs boson and several long-lived scalar masses, for the three channels combined. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons or two muons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S}, \mathrm{S} \to \ell^{+}\ell^{-}$ as a function of $c\tau_0$, for a 400 GeV Higgs boson and several long-lived scalar masses, for the three channels combined. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons or two muons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S}, \mathrm{S} \to \ell^{+}\ell^{-}$ as a function of $c\tau_0$, for a 600 GeV Higgs boson and several long-lived scalar masses, for the three channels combined. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons or two muons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S}, \mathrm{S} \to \ell^{+}\ell^{-}$ as a function of $c\tau_0$, for a 800 GeV Higgs boson and several long-lived scalar masses, for the three channels combined. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons or two muons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S}, \mathrm{S} \to \ell^{+}\ell^{-}$ as a function of $c\tau_0$, for a 1000 GeV Higgs boson and several long-lived scalar masses, for the three channels combined. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons or two muons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

The 95% CL upper limits on the $\mathrm{H} \to \mathrm{S}\mathrm{S} \to 4\ell$ branching fraction as a function of $c\tau_0$, for a Higgs boson with a mass of 125 GeV and a long-lived scalar with a mass of 30 GeV or 50 GeV, for the ee channel only. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons. The area above the solid (dashed) curve represents the observed (expected) exclusion region. The curves are not smooth because very few ee channel events pass the preselection for the Higgs boson and scalar masses shown in this figure.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S} \to 4\ell$ as a function of $c\tau_0$, for a 600 GeV Higgs boson and several long-lived scalar masses, for the ee channel only. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S} \to 4\ell$ as a function of $c\tau_0$, for a 1000 GeV Higgs boson and several long-lived scalar masses, for the ee channel only. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two electrons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

The 95% CL upper limits on the $\mathrm{H} \to \mathrm{S}\mathrm{S} \to 4\ell$ branching fraction as a function of $c\tau_0$, for a Higgs boson with a mass of 125 GeV and a long-lived scalar with a mass of 30 GeV or 50 GeV, for the $\mu\mu$ channel only. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two muons. The area above the solid (dashed) curve represents the observed (expected) exclusion region.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S} \to 4\ell$ as a function of $c\tau_0$, for a 600 GeV Higgs boson and several long-lived scalar masses, for the $\mu\mu$ channel only. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two muons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.

The 95% CL upper limits on the product of the cross section and branching fraction $\mathrm{H} \to \mathrm{S}\mathrm{S} \to 4\ell$ as a function of $c\tau_0$, for a 1000 GeV Higgs boson and several long-lived scalar masses, for the $\mu\mu$ channel only. These limits assume that $\mathcal{B}(H \to SS)=100\%$ and each S has a 50% probability of decaying to two muons. The area above the solid (dashed) curves represents the observed (expected) exclusion region.