$\xi$ distributions for different jet $p_T$ bins.

$\xi$ distributions for different jet $p_T$ bins.

The $j_T/p_{T}^{\rm jet}$ distributions for different jet $p_T$ bins.

The $j_T/p_{T}^{\rm jet}$ distributions for different jet $p_T$ bins.

$r$ distributions for different jet $p_T$ bins.

Normalized $p_{\mathrm{T}}$-spectrum ratio as a function as a function of centrality in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}} \in$ $0.5 \leq |\eta|\leq 0.8$.

Normalized $p_{\mathrm{T}}$-spectrum ratio as a function as a function of centrality in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{tracklets}} \in$ $0.5 \leq |\eta|\leq 0.8$.

Normalized $p_{\mathrm{T}}$-spectrum ratio as a function as a function of centrality in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}} \in$ $-3.7<\eta<-1.7$ and $2.8 < \eta <5.1$.

Event fraction distribution as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}} \in -3.7<\eta<-1.7$ and $2.8 < \eta <5.1$.

Event fraction distribution as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}} \in |\eta|\leq 0.8$.

Event fraction distribution as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}} \in |\eta|\leq 0.8$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $|\eta| \leq 0.8$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $0.5 \leq |\eta| \leq 0.8$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}}$ $\in$ $|\eta| \leq 0.8$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}}$ $\in$ $0.5 \leq |\eta| \leq 0.8$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{tracklets}}$ $\in$ $|\eta| \leq 0.8$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{tracklets}}$ $\in$ $0.5 \leq |\eta| \leq 0.8$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{tracklets}}$ $\in$ $0.3 \leq |\eta| \leq 0.6$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{tracklets}}$ $\in$ $0.7 \leq |\eta| \leq 1$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $-3.7<\eta<-1.7$ and $2.8 < \eta <5.1$.

Extracted $c_{\mathrm{s}}^{2}$ as a function of the minimum pseudorapidity gap between the centrality estimation and transverse-momentum spectra measurement in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$.

$k_{2}/k_{2}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $|\eta| \leq 0.8$.

$k_{2}/k_{2}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}}$ $\in$ $|\eta| \leq 0.8$.

$k_{2}/k_{2}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $-3.7<\eta<-1.7$ and $2.8 < \eta <5.1$.

$k_{3}/k_{3}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $|\eta| \leq 0.8$.

$k_{3}/k_{3}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}}$ $\in$ $|\eta| \leq 0.8$.

$k_{3}/k_{3}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $-3.7<\eta<-1.7$ and $2.8 < \eta <5.1$.

$\gamma_{\langle [p_{\mathrm{T}}] \rangle}/\gamma_{\langle [p_{\mathrm{T}}] \rangle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $|\eta| \leq 0.8$.

$\gamma_{\langle [p_{\mathrm{T}}] \rangle}/\gamma_{\langle [p_{\mathrm{T}}] \rangle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}}$ $\in$ $|\eta| \leq 0.8$.

$\gamma_{\langle [p_{\mathrm{T}}] \rangle}/\gamma_{\langle [p_{\mathrm{T}}] \rangle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $-3.7<\eta<-1.7$ and $2.8 < \eta <5.1$.

$\Gamma_{\langle [p_{\mathrm{T}}] \rangle}/\Gamma_{\langle [p_{\mathrm{T}}] \rangle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $|\eta| \leq 0.8$.

$\Gamma_{\langle [p_{\mathrm{T}}] \rangle}/\Gamma_{\langle [p_{\mathrm{T}}] angle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}}$ $\in$ $|\eta| \leq 0.8$.

$\Gamma_{\langle [p_{\mathrm{T}}] \rangle}/\Gamma_{\langle [p_{\mathrm{T}}] \rangle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $-3.7<\eta<-1.7$ and $2.8 < \eta <5.1$.

$\kappa_{\langle [p_{\mathrm{T}}] \rangle}/\kappa_{\langle [p_{\mathrm{T}}] \rangle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $|\eta| \leq 0.8$.

$\kappa_{\langle [p_{\mathrm{T}}] \rangle}/\kappa_{\langle [p_{\mathrm{T}}] \rangle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}}$ $\in$ $|\eta| \leq 0.8$.

$\kappa_{\langle [p_{\mathrm{T}}] \rangle}/\kappa_{\langle [p_{\mathrm{T}}] \rangle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $-3.7<\eta<-1.7$ and $2.8 < \eta <5.1$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $0.5 \leq |\eta| \leq 0.8$ with $0.45 \leq p_{\mathrm{T}} \leq 50~\mathrm{GeV}/c.$.

$R_{\rm p}$ vs. $m_{\rm T}$ for $0-10\%$ centrality for Pb-Pb collisions at 5020 GeV.

$R_{\rm p}$ vs. $m_{\rm T}$ for $10-30\%$ centrality for Pb-Pb collisions at 5020 GeV.

$R_{\rm p}$ vs. $m_{\rm T}$ for $30-50\%$ centrality for Pb-Pb collisions at 5020 GeV.

$R_{\rm N}$ vs. $m_{\rm T}$ for $0-10\%$ centrality for Pb-Pb collisions at 5020 GeV.

$R_{\rm N}$ vs. $m_{\rm T}$ for $10-30\%$ centrality for Pb-Pb collisions at 5020 GeV.

$R_{\rm N}$ vs. $m_{\rm T}$ for $30-50\%$ centrality for Pb-Pb collisions at 5020 GeV.

$R_{\rm d}$ vs. $m_{\rm T}$ for $0-10\%$ centrality for Pb-Pb collisions at 5020 GeV.

$R_{\rm d}$ vs. $m_{\rm T}$ for $10-30\%$ centrality for Pb-Pb collisions at 5020 GeV.

$R_{\rm d}$ vs. $m_{\rm T}$ for $30-50\%$ centrality for Pb-Pb collisions at 5020 GeV.

$R_{\rm p}$ vs. '<dN_{ch}/d \eta>^{1/3}' for mT bin 0 for Pb-Pb collisions at 5020 GeV.

$R_{\rm p}$ vs. '<dN_{ch}/d \eta>^{1/3}' for mT bin 1 for Pb-Pb collisions at 5020 GeV.

$R_{\rm p}$ vs. '<dN_{ch}/d \eta>^{1/3}' for mT bin 2 for Pb-Pb collisions at 5020 GeV.

$R_{\rm p}$ vs. '<dN_{ch}/d \eta>^{1/3}' for mT bin 3 for Pb-Pb collisions at 5020 GeV.

$R_{\rm p}$ vs. '<dN_{ch}/d \eta>^{1/3}' for mT bin 4 for Pb-Pb collisions at 5020 GeV.

$R_{\rm p}$ vs. '<dN_{ch}/d \eta>^{1/3}' for mT bin 5 for Pb-Pb collisions at 5020 GeV.

$R_{\rm p}$ vs. '<dN_{ch}/d \eta>^{1/3}' for mT bin 6 for Pb-Pb collisions at 5020 GeV.

$R_{\rm p}$ vs.'<dN_{ch}/d \eta>^{1/3}' for mT bin 7 for Pb-Pb collisions at 5020 GeV.

proton-deuteron (same charge) correlation function for centrality 0-10% from Pb-Pb collisions at 5020 GeV

proton-deuteron (same charge) correlation function for centrality 10-30% from Pb-Pb collisions at 5020 GeV

proton-deuteron (same charge) correlation function for centrality 30-50% from Pb-Pb collisions at 5020 GeV

The extracted kinetic decoupling temperature is derived from $\pi^{-}$–d correlation functions.

This is the mixed event distribution for $\pi^{+}$–d, required for the fitting

This is the mixed event distribution for $\pi^{-}$–d, required for the fitting

The momentum smearing matrix, required for the fitting, used both for $\pi^{-}$–d and $\pi^{+}$–d

$\Delta R_{\rm axis}$ distribution for ${\rm STD-WTA / STD-D^0}$ for $\rm D^0$-tagged jets of $R=0.4$, in the intervals $10<p_{\rm T}^{\rm ch \ jet}<20 \ {\rm GeV}/c$ and $5<p_{\rm T}^{\rm D^0}<20 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for $\rm STD-D^0$ for $\rm D^0$-tagged jets of $R=0.4$, in the intervals $10<p_{\rm T}^{\rm ch \ jet}<20 \ {\rm GeV}/c$ and $5<p_{\rm T}^{\rm D^0}<20 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for $\rm SD-D^0$ with grooming setting ($z_{\rm cut} = 0.1, \beta = 0$) for $\rm D^0$-tagged jets of $R=0.4$, in the intervals $10<p_{\rm T}^{\rm ch \ jet}<20 \ {\rm GeV}/c$ and $5<p_{\rm T}^{\rm D^0}<20 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for $\rm SD-D^0$ with grooming setting ($z_{\rm cut} = 0.2, \beta = 0$) for $\rm D^0$-tagged jets of $R=0.4$, in the intervals $10<p_{\rm T}^{\rm ch \ jet}<20 \ {\rm GeV}/c$ and $5<p_{\rm T}^{\rm D^0}<20 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for $\rm SD-D^0$ with grooming setting ($z_{\rm cut} = 0.1, \beta = 0$) for $\rm D^0$-tagged jets of $R=0.4$, in the intervals $10<p_{\rm T}^{\rm ch \ jet}<20 \ {\rm GeV}/c$ and $5<p_{\rm T}^{\rm D^0}<20 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for $\rm WTA-SD$ with grooming setting ($z_{\rm cut} = 0.1, \beta = 0$) for $\rm D^0$-tagged jets of $R=0.4$, in the intervals $10<p_{\rm T}^{\rm ch \ jet}<20 \ {\rm GeV}/c$ and $5<p_{\rm T}^{\rm D^0}<20 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for ${\rm WTA-SD / SD-D^0}$ with grooming setting ($z_{\rm cut} = 0.1, \beta = 0$) for $\rm D^0$-tagged jets of $R=0.4$, in the intervals $10<p_{\rm T}^{\rm ch \ jet}<20 \ {\rm GeV}/c$ and $5<p_{\rm T}^{\rm D^0}<20 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for $\rm WTA-SD$ with grooming setting ($z_{\rm cut} = 0.2, \beta = 0$) for $\rm D^0$-tagged jets of $R=0.4$, in the intervals $10<p_{\rm T}^{\rm ch \ jet}<20 \ {\rm GeV}/c$ and $5<p_{\rm T}^{\rm D^0}<20 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for ${\rm WTA-SD / SD-D^0}$ with grooming setting ($z_{\rm cut} = 0.2, \beta = 0$) for $\rm D^0$-tagged jets of $R=0.4$, in the intervals $10<p_{\rm T}^{\rm ch \ jet}<20 \ {\rm GeV}/c$ and $5<p_{\rm T}^{\rm D^0}<20 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for $\rm STD-SD$ with grooming setting ($z_{\rm cut} = 0.1, \beta = 0$) for $\rm D^0$-tagged jets of $R=0.4$, in the intervals $10<p_{\rm T}^{\rm ch \ jet}<20 \ {\rm GeV}/c$ and $5<p_{\rm T}^{\rm D^0}<20 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for $\rm STD-SD$ with grooming setting ($z_{\rm cut} = 0.2, \beta = 0$) for $\rm D^0$-tagged jets of $R=0.4$, in the intervals $10<p_{\rm T}^{\rm ch \ jet}<20 \ {\rm GeV}/c$ and $5<p_{\rm T}^{\rm D^0}<20 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for the ratio of $\rm STD-SD$ ($z_{\rm cut} = 0.2, \beta=0$) over $\rm STD-SD$ ($z_{\rm cut} = 0.1, \beta=0$) for $\rm D^0$-tagged jets of $R=0.4$, in the intervals $10<p_{\rm T}^{\rm ch \ jet}<20 \ {\rm GeV}/c$ and $5<p_{\rm T}^{\rm D^0}<20 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for the ratio of $\rm SD-D^0$ ($z_{\rm cut} = 0.2, \beta=0$) over $\rm SD-D^0$ ($z_{\rm cut} = 0.1, \beta=0$) for $\rm D^0$-tagged jets of $R=0.4$, in the intervals $10<p_{\rm T}^{\rm ch \ jet}<20 \ {\rm GeV}/c$ and $5<p_{\rm T}^{\rm D^0}<20 \ {\rm GeV}/c$.

$p_{T}$ dependence of $v_{2}$ for $\pi^{+}$ in Au+Au collisions at 3.9 GeV

$p_{T}$ dependence of $v_{2}$ for $\pi^{+}$ in Au+Au collisions at 4.5 GeV

$p_{T}$ dependence of $v_{2}$ for $\pi^{-}$ in Au+Au collisions at 3 GeV

$p_{T}$ dependence of $v_{2}$ for $\pi^{-}$ in Au+Au collisions at 3.2 GeV

$p_{T}$ dependence of $v_{2}$ for $\pi^{-}$ in Au+Au collisions at 3.5 GeV

$p_{T}$ dependence of $v_{2}$ for $\pi^{-}$ in Au+Au collisions at 3.9 GeV

$p_{T}$ dependence of $v_{2}$ for $\pi^{-}$ in Au+Au collisions at 4.5 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{+}$ in Au+Au collisions at 3 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{+}$ in Au+Au collisions at 3.2 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{+}$ in Au+Au collisions at 3.5 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{+}$ in Au+Au collisions at 3.9 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{+}$ in Au+Au collisions at 4.5 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{-}$ in Au+Au collisions at 3 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{-}$ in Au+Au collisions at 3.2 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{-}$ in Au+Au collisions at 3.5 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{-}$ in Au+Au collisions at 3.9 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{-}$ in Au+Au collisions at 4.5 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{0}_{S}$ in Au+Au collisions at 3 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{0}_{S}$ in Au+Au collisions at 3.2 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{0}_{S}$ in Au+Au collisions at 3.5 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{0}_{S}$ in Au+Au collisions at 3.9 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{0}_{S}$ in Au+Au collisions at 4.5 GeV

$p_{T}$ dependence of $v_{2}$ for $p$ in Au+Au collisions at 3 GeV

$p_{T}$ dependence of $v_{2}$ for $p$ in Au+Au collisions at 3.2 GeV

$p_{T}$ dependence of $v_{2}$ for $p$ in Au+Au collisions at 3.5 GeV

$p_{T}$ dependence of $v_{2}$ for $p$ in Au+Au collisions at 3.9 GeV

$p_{T}$ dependence of $v_{2}$ for $p$ in Au+Au collisions at 4.5 GeV

$p_{T}$ dependence of $v_{2}$ for $\Lambda$ in Au+Au collisions at 3 GeV

$p_{T}$ dependence of $v_{2}$ for $\Lambda$ in Au+Au collisions at 3.2 GeV

$p_{T}$ dependence of $v_{2}$ for $\Lambda$ in Au+Au collisions at 3.5 GeV

$p_{T}$ dependence of $v_{2}$ for $\Lambda$ in Au+Au collisions at 3.9 GeV

$p_{T}$ dependence of $v_{2}$ for $\Lambda$ in Au+Au collisions at 4.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{+}$ at 3 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{+}$ at 3.2 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{+}$ at 3.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{+}$ at 3.9 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{+}$ at 4.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{-}$ at 3 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{-}$ at 3.2 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{-}$ at 3.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{-}$ at 3.9 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{-}$ at 4.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{+}$ at 3 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{+}$ at 3.2 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{+}$ at 3.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{+}$ at 3.9 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{+}$ at 4.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{-}$ at 3 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{-}$ at 3.2 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{-}$ at 3.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{-}$ at 3.9 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{-}$ at 4.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{0}_{S}$ at 3 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{0}_{S}$ at 3.2 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{0}_{S}$ at 3.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{0}_{S}$ at 3.9 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{0}_{S}$ at 4.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $p$ at 3 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $p$ at 3.2 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $p$ at 3.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $p$ at 3.9 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $p$ at 4.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\Lambda$ at 3 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\Lambda$ at 3.2 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\Lambda$ at 3.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\Lambda$ at 3.9 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\Lambda$ at 4.5 GeV

The energy dependence of $p_T$ integrated $v_2$ for $\pi^{\pm}$, $K^{\pm}$, $K^{0}_{S}$, $p$, and $\Lambda$ in 10-40\% centrality from Au+Au collisions at $\sqrt{s_{\text{NN}}}$ = 3.0, 3.2, 3.5, 3.9, and 4.5 GeV.

The energy dependence of $n_{q}$ scaled $v_{2}$ ratios of $v_{2}^{q}(\pi^{+}) / v_{2}^{q}(K^{+})$ and $v_{2}^{q}(p) / v_{2}^{q}(K^{+})$ at ($m_T - m_0$)/$n_q$ = 0.4 GeV/$c^2$

UrQMD model calculations for net-proton cumulant ratios and proton factorial cumulant ratios in Au+Au collisions from $\sqrt{s_{NN}}$ = 7.7 - 200 GeV for 0-5$\%$ centrality. Same choices of centrality definition are used as used for experimental data.

Significance of deviations of central 0-5$\%$ data from various choices of baselines or references.

Cumulant ratios C2/C1 and C4/C2 of net-proton distributions in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV as a function of multiplicity (RefMult3X).

Cumulant ratios C2/C1 and C4/C2 of net-proton distributions in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV as a function of average multiplicity (RefMult3X) for 8 centrality bins from 0-10$\%$ to 70-80$\%$. Data with and without CBWC are given.

Cumulant ratios C2/C1 and C4/C2 of net-proton distributions in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV as a function of average multiplicity (RefMult3X) for 16 centrality bins from 0-5$\%$ to 75-80$\%$. Data with and without CBWC are given.

Reference multiplicity distributions in Au+Au collisions at $\sqrt{s_{NN}}$ = 11.5 GeV from BES-II. Three choices of multiplicity distributions (RefMult3 scaled by efficiency factor of $\epsilon$ = 0.83 or 83$\%$ to mimic RefMult3 of BES-I, RefMult3, and RefMult3X) are presented.

Cumulant ratios C2/C1 and C4/C2 of net-proton distributions in Au+Au collisions at $\sqrt{s_{NN}}$ = 11.5 GeV from BES-II as a function of centrality obtained with various centrality definitions.

Reference multiplicity distributions (RefMult3, RefMult3X, and RefMult3A) in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV from UrQMD model. The RefMult3 and RefMult3X distributions have been scaled by an efficiency factor of $\epsilon$ = 77.9$\%$ to mimic the distribution in experimenal data.

Cumulant ratios C2/C1 and C4/C2 of net-proton distributions in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV from UrQMD model as a function of centrality obtained with various centrality definitions (RefMult3, RefMult3X, and RefMult3A). The RefMult3 and RefMult3X distributions have been scaled by an efficiency factor of $\epsilon$ = 77.9$\%$ to mimic the distribution in experimenal data.

Collision energy dependence of net-proton C4/C2 in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7 - 27 GeV from BES-II obtained with RefMult3 centrality definition.

Difference between net-proton C4/C2 (0-5$\%$) data from various choices of baselines or references as a function of collision energy from $\sqrt{s_{NN}}$ = 7.7 - 200 GeV.

The distributions of the variables used in the simultaneous fit for the SR. The black points with error bars represent the data and their statistical uncertainties, whereas the shaded band represents the predicted uncertainties. The bottom panel in each figure shows the ratio of the number of events observed in data to that of the total SM prediction. The last bin of each plot has been extended to include the overflow contribution.

Expected and observed 95% upper limits on the product of the cross section and branching fraction$\sigma(\mathrm{pp} \ ightarrow W\mathrm{a})\mathcal{B}(W ightarrow \ell^{+} v_{\ell})\mathcal{B}(\mathrm{a} \rightarrow Z\gamma)\mathcal{B}(Z \rightarrow \ell^{+} \ell^{-})$ s a function of the ALP mass from 110 to 400 GeV

Expected and observed 95% upper limits on the photophobic ALP model parameter $1/f_{\mathrm{a}}$ as a function of ALP mass reinterpreted from $1/f_{\mathrm{a}}=2~\mathrm{TeV}^{-1}$