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

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.

Balance function $B^{\pi\pi}$ projected onto the $\Delta y$ axis for particle pairs within the full range $|\Delta\varphi|\leq \pi$ in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{{\rm K}\pi}$ projected onto the $\Delta y$ axis for particle pairs within the full range $|\Delta\varphi|\leq \pi$ in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{{\rm p}\pi}$ projected onto the $\Delta y$ axis for particle pairs within the full range $|\Delta\varphi|\leq \pi$ in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{\pi{\rm K}}$ projected onto the $\Delta y$ axis for particle pairs within the full range $|\Delta\varphi|\leq \pi$ in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{{\rm KK}}$ projected onto the $\Delta y$ axis for particle pairs within the full range $|\Delta\varphi|\leq \pi$ in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{{\rm pK}}$ projected onto the $\Delta y$ axis for particle pairs within the full range $|\Delta\varphi|\leq \pi$ in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{\pi{\rm p}}$ projected onto the $\Delta y$ axis for particle pairs within the full range $|\Delta\varphi|\leq \pi$ in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{{\rm Kp}}$ projected onto the $\Delta y$ axis for particle pairs within the full range $|\Delta\varphi|\leq \pi$ in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{{\rm pp}}$ projected onto the $\Delta y$ axis for particle pairs within the full range $|\Delta\varphi|\leq \pi$ in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{\pi\pi}$ projected onto the $\Delta \varphi$ axis in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{{\rm K}\pi}$ projected onto the $\Delta \varphi$ axis in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{{\rm p}\pi}$ projected onto the $\Delta \varphi$ axis in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{\pi{\rm K}}$ projected onto the $\Delta \varphi$ axis in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{{\rm KK}}$ projected onto the $\Delta \varphi$ axis in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{{\rm pK}}$ projected onto the $\Delta \varphi$ axis in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{\pi{\rm p}}$ projected onto the $\Delta \varphi$ axis in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{{\rm Kp}}$ projected onto the $\Delta \varphi$ axis in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Balance function $B^{{\rm pp}}$ projected onto the $\Delta \varphi$ axis in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$ for selected centrality ranges ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Longitudinal ($\Delta y$) $\sigma$ widths of balance functions of the full species matrix of $\pi^{\pm}$, ${\rm K}^{\pm}$, and ${\rm p/\bar{p}}$ with centrality. Pairs ${\rm K}\pi$, ${\rm p}\pi$, and ${\rm pK}$ have the same values with $\pi{\rm K}$, $\pi{\rm p}$, and ${\rm Kp}$, respectively. The relative azimuthal angle range for all the species pairs is the full azimuth range $|\Delta\varphi|\leq\pi$ ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Azimuthal ($\Delta\varphi$) $\sigma$ widths of balance functions of the full species matrix of $\pi^{\pm}$, ${\rm K}^{\pm}$, and ${\rm p/\bar{p}}$ with centrality. Pairs ${\rm K}\pi$, ${\rm p}\pi$, and ${\rm pK}$ have the same values with $\pi{\rm K}$, $\pi{\rm p}$, and ${\rm Kp}$, respectively. The relative rapidity range used for all species pairs is $|\Delta y|\leq 1.2$, with the exception of $|\Delta y|\leq 1.4$ for $\pi\pi$ and $|\Delta y|\leq 1.0$ for $\rm pp$ ($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

Integrals of balance functions of the full species matrix of $\pi^{\pm}$, ${\rm K}^{\pm}$, and ${\rm p/\bar{p}}$ with centrality($\pi,{\rm K}: 0.2 \leq p_{\rm T} \leq 2.0\;{\rm GeV}/c$; ${\rm p}: 0.5 \leq p_{\rm T} \leq 2.5\;{\rm GeV}/c$).

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

Invariant mass distribution of the diphoton pairs for the elastic selection region with events satisfying a < 0.005 (signal contribution is magnified by a factor 5000).

Predicted number of events having an elastic diphoton pair in association with a pair of protons observed within the range where the proton detectors have a radiation inefficiency less than 10%. The yields where the two-photon and two-proton systems mass and rapidity are matching at 2 and 3$\sigma$ are also quoted. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

Predicted number of events having an elastic diphoton pair in association with a pair of protons observed within the range where the proton detectors have a radiation inefficiency less than 10%. The yields where the two-photon and two-proton systems mass and rapidity are matching at 2 and 3$\sigma$ are also quoted. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

95% CL expected and observed frequentist upper limit for the LbL cross section within the $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$ search region.

95% CL expected and observed frequentist upper limit for the LbL cross section within the $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$ search region.

95% expected and observed one-dimensional limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters, when the other parameter is set to zero. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

95% expected and observed one-dimensional limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters, when the other parameter is set to zero. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

95% expected and observed limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters. The parametric elliptic form is assumed: $a_0\zeta_1^2+a_1\zeta_1\zeta_2+a_2\zeta_2^2=1$. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

95% expected and observed limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters. The parametric elliptic form is assumed: $a_0\zeta_1^2+a_1\zeta_1\zeta_2+a_2\zeta_2^2=1$. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

Correlation coefficients $\rho_{i,j}$ for the statistical and systematic uncertainties between the $i$-th and $j$-th bin of the differential $A_E$ measurement as a function of the jet scattering angle $\theta_j$

The effects on the energy asymmetry of $1\sigma$ variations in its influencing nuisance parameters for the three $\theta_j$ bins. These are extracted from the samples of the posterior distribution with $\sigma_i^{(j)} = c_{ij}/\sqrt{c_{jj}}$ being the estimated shift of bin $i$ in conjunction with a shift $\Delta\theta_j$ of nuisance parameter $j$. The data statistical (Data stat.) uncertainty is obtained from running the unfolding with all nuisance parameters being fixed to their post-marginalised values, the MC statistical uncertainty on the response matrix ($t\bar{t}$ response MC stat.) is evaluated using a bootstrapping method from the covariance matrix of the ensemble of repeated unfolding results with varied response matrices. The $\gamma$ variations denote the statistical uncertainties of the background predictions in the corresponding bin of the $\Delta E$ vs $\theta_{j}$ distribution. The numbers appended to the $W$+jets PDF variations denote the corresponding NNPDF3.0 PDF sets.

The effects on the energy asymmetry of $1\sigma$ variations in its influencing nuisance parameters for the three $\theta_j$ bins. These are extracted from the samples of the posterior distribution with $\sigma_i^{(j)} = c_{ij}/\sqrt{c_{jj}}$ being the estimated shift of bin $i$ in conjunction with a shift $\Delta\theta_j$ of nuisance parameter $j$. The data statistical (Data stat.) uncertainty is obtained from running the unfolding with all nuisance parameters being fixed to their post-marginalised values, the MC statistical uncertainty on the response matrix ($t\bar{t}$ response MC stat.) is evaluated using a bootstrapping method from the covariance matrix of the ensemble of repeated unfolding results with varied response matrices. The $\gamma$ variations denote the statistical uncertainties of the background predictions in the corresponding bin of the $\Delta E$ vs $\theta_{j}$ distribution. The numbers appended to the $W$+jets PDF variations denote the corresponding NNPDF3.0 PDF sets.

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

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

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.

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.

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.

Upper limits on $\mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ at 95% CL including the BDT cut as a function of the signal mass.

Upper limits on $\mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ at 95% CL without the BDT cut as a function of the signal mass.

The number of events when applying each cut in the analysis for $m_a$ = 16 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 20 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 30 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 40 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 50 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 60 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

Post-fit event yields in the signal and control regions obtained from a scan over $\epsilon_{VV}$ exploiting both shape and rate information. The quoted uncertainties include the theoretical and experimental systematic sources and those due to sample statistics. The fit constrains the total expected yield to the observed yield. The diboson background is split into $W W$ and non-$W W$ contributions.

Best-fit values and their uncertainties as obtained from the shape-only and shape-plus-rate likelihood fits to the Asimov dataset and to ATLAS data. Results of both shape-only and shape+rate fits for $a_L$ and $a_T$ are shown. Results of fits to one parameter with the other one fixed or profiled are presented.

Best-fit values and their uncertainties as obtained from the shape-only and shape-plus-rate likelihood fits to the Asimov dataset and to ATLAS data. Results of both shape-only and shape+rate fits for $\epsilon_{VV}$ and $\kappa_{VV}$ are shown. Results of fits to one parameter with the other one fixed or profiled are presented.

The contributions of the leading individual systematic uncertainties together with the data statistical uncertainties, in the one dimensional fit for the pseudo-observables $\kappa_{VV}$ (a) and $\epsilon_{VV}$ (b) for electroweak-boson polarisation in the VBF $H\to WW$ channel. Both shape and rate informations are exploited in the fit. The theoretical and experimental uncertainties are subdivided further into categories.

The contributions of the leading individual systematic uncertainties together with the data statistical uncertainties, in the one dimensional fit for the pseudo-observables $\kappa_{VV}$ (a) and $\epsilon_{VV}$ (b) for electroweak-boson polarisation in the VBF $H\to WW$ channel. Both shape and rate informations are exploited in the fit. The theoretical and experimental uncertainties are subdivided further into categories.

Post-fit distribution of the BDT response observable presented in the four $|\Delta \eta jj|$ categories of the ggF +2 jets signal region, with signal and background yields fixed from the fit for $\mu^{\text{ggF+2jets}}$. The distributions of the ggF + 2 jets and VBF processes are overlaid with their respective contributions multiplied by 50.

The weighted $\Delta \Phi_{jj}$ post-fit distribution in the ggF +2 jets signal region, with signal and background yields fixed from the fit to $\tan \alpha$ using shape and rate information.

Expected and observed likelihood curves for scans over $\tan \alpha$ where only the shape is taken into account in the fit, $\mu_{VBF}$ is fixed.

Expected and observed likelihood curves for scans over $\tan \alpha$ where both shape and normalisation are taken into account in the fit, $\mu_{VBF}$ is fixed.

68% and 95% CL two-dimensional likelihood contours of the CP-even and CP-odd coupling parameters $K_{gg} \cos(\alpha)$ and $K_{gg} \sin(\alpha)$. The minima are represented by black stars, while the SM value is shown as a red star.

The weighted $\Delta \Phi jj$ distribution in the VBF signal region, with signal and background yields fixed from the fit for $\epsilon_{VV}$ using shape and rate information.

Likelihood scans over the transversally polarised couplings. The fit is using shape-only information. All relevant experimental and modelling systematic uncertainties are considered in the fit.

Likelihood scans over the transversally polarised couplings. The fit is using shape + rate information. All relevant experimental and modelling systematic uncertainties are considered in the fit.

Likelihood scans over the longitudinally polarised couplings. The fit is using shape + rate information. All relevant experimental and modelling systematic uncertainties are considered in the fit.

Likelihood scans over $\kappa_{VV}$ with the $\epsilon_{VV}$ profiled. The fit is performed using both shape and rate information. All relevant experimental and theoretical systematic uncertainties are considered in the fit.

Likelihood scans over $\epsilon_{VV}$ with the $\kappa_{VV}$ profiled. The fit is performed using both shape and rate information. All relevant experimental and theoretical systematic uncertainties are considered in the fit.

The contributions of the leading individual systematic uncertainties together with the data statistical uncertainties, in the one dimensional fit for electroweak-boson polarisation in the VBF $H\to WW$ channel, using (aL, aT) parametrisation. Both shape and rate informations are exploited in the fit. The theoretical and experimental uncertainties are subdivided further into categories.

The contributions of the leading individual systematic uncertainties together with the data statistical uncertainties, in the one dimensional fit for electroweak-boson polarisation in the VBF $H\to WW$ channel, using (aL, aT) parametrisation.. Both shape and rate informations are exploited in the fit. The theoretical and experimental uncertainties are subdivided further into categories.