Example of combinatorial background subtracted invariant mass distributions and the extracted yields as a function of $\cos \theta^*$ for $\phi$ and $K^{*0}$ mesons. \textbf{a)} example of $\phi \rightarrow K^+ + K^-$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{b)} example of $K^{*0} (\overline{K^{*0}}) \rightarrow K^{-} \pi^{+} (K^{+} \pi^{-})$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{c)} extracted yields of $\phi$ as a function of $\cos \theta^*$; \textbf{d)} extracted yields of $K^{*0}$ as a function of $\cos \theta^*$.

Example of combinatorial background subtracted invariant mass distributions and the extracted yields as a function of $\cos \theta^*$ for $\phi$ and $K^{*0}$ mesons. \textbf{a)} example of $\phi \rightarrow K^+ + K^-$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{b)} example of $K^{*0} (\overline{K^{*0}}) \rightarrow K^{-} \pi^{+} (K^{+} \pi^{-})$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{c)} extracted yields of $\phi$ as a function of $\cos \theta^*$; \textbf{d)} extracted yields of $K^{*0}$ as a function of $\cos \theta^*$.

Example of combinatorial background subtracted invariant mass distributions and the extracted yields as a function of $\cos \theta^*$ for $\phi$ and $K^{*0}$ mesons. \textbf{a)} example of $\phi \rightarrow K^+ + K^-$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{b)} example of $K^{*0} (\overline{K^{*0}}) \rightarrow K^{-} \pi^{+} (K^{+} \pi^{-})$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{c)} extracted yields of $\phi$ as a function of $\cos \theta^*$; \textbf{d)} extracted yields of $K^{*0}$ as a function of $\cos \theta^*$.

Efficiency corrected $\phi$-meson yields as a function of cos$\theta$* and corresponding fits with Eq.1 in the method section.

Efficiency and acceptance corrected $K^{*0}$-meson yields as a function of cos$\theta$* and corresponding fits with Eq.4 in the method section.

$\phi$-meson $\rho_{00}$ obtained from 1st- and 2nd-order event planes. The red stars (gray squares) show the $\phi$-meson $\rho_{00}$ as a function of beam energy, obtained with the 2nd-order (1st-order) EP.

$\phi$-meson $\rho_{00}$ with respect to different quantization axes. $\phi$-meson $\rho_{00}$ as a function of beam energy, for the out-of-plane direction (stars) and the in-plane direction (diamonds). Curves are fits based on theoretical calculations with a $\phi$-meson field. The corresponding $G_s^{(y)}$ values obtained from the fits are shown in the legend.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

Global spin alignment measurement of $\phi$ and $K^{*0}$ vector mesons in Au+Au collisions at 0-20\% centrality. The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for $K^{*0}$-meson, obtained with the 2nd-order EP.

Global spin alignment measurement of $\phi$ and $K^{*0}$ vector mesons in Au+Au collisions at 0-20\% centrality. The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for $K^{*0}$-meson, obtained with the 2nd-order EP.

Direct photon spectra(Nucl. Part. Phys. 23, A1 (1997)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in p+p at $\sqrt{s_{NN}}$= 63 GeV.

Direct photon spectra(Sov. J. Nucl. Phys. 51, 836 (1990)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in p+p at $\sqrt{s_{NN}}$= 63 GeV.

Direct photon spectra normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 62.4 GeV.

Direct photon spectra normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 39.0 GeV.

Direct photon spectra(Physical Review C91, 064904 (2015)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 200 GeV.

Direct photon spectra(Physical Review Letters 104, 132301 (2010)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 200 GeV.

Direct photon spectra(Physical Review Letters 109, 152302 (2012)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 200 GeV.

Direct photon spectra(Physical Review C98, 054902 (2018)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Cu+Cu at $\sqrt{s_{NN}}$= 200 GeV.

Direct photon spectra(Physics Letters B754 235 (2016)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Pb+Pb at $\sqrt{s_{NN}}$= 2760 GeV.

Number of binary collisions N_{coll} vs. Charged particle multiplicity in Au+Au at various cm energies.

Charged particle multiplicity vs. centrality in Au+Au at various cm energies.

Number of binary collisions N_{coll} vs. centrality in Au+Au at various cm energies.

Integrated direct photon yield vs. Charged particle multiplicity in various systems at various cm energies.

Charged particle multiplicity vs. centrality in various systems at various cm energies.

Integrated direct photon yield vs. centrality in various systems at various cm energies.

Integrated direct photon yield vs. Charged particle multiplicity in various systems at various cm energies.

Charged particle multiplicity vs. centrality in various systems at various cm energies.

Integrated direct photon yield vs. centrality in various systems at various cm energies.

Differential production cross sections of $\psi(2S)$ as a function of rapidity.

$\psi(2S)$ over $J/\psi$ ratio as a function of $p_{\rm T}$.

$\psi(2S)$ over $J/\psi$ ratio as a function of y.

Differential production cross sections of $J/\psi$ as a function of $p_{\rm T}$.

Differential production cross sections of $J/\psi$ as a function of rapidity.

$J/\psi$ mean transversee momentum vs collision energy. $p_{\rm T}$ integration ranges are 0<$p_{\rm T}$<8 GeV/$c$ at $\sqrt{s}$ =2700 GeV, 0<$p_{\rm T}$<12 GeV/$c$ at $\sqrt{s}$ =5020, 0<$p_{\rm T}$<20 GeV/$c$ at $\sqrt{s}$ =7000, 0<$p_{\rm T}$<20 GeV/$c$ at $\sqrt{s}$ =8000 and 0<$p_{\rm T}$<20 GeV/$c$ at $\sqrt{s}$ =13000.

$J/\psi$ mean transversee momentum square vs collision energy. $p_{\rm T}$ integration ranges are 0<$p_{\rm T}$<8 GeV/$c$ at $\sqrt{s}$ =2700 GeV, 0<$p_{\rm T}$<12 GeV/$c$ at $\sqrt{s}$ =5020, 0<$p_{\rm T}$<20 GeV/$c$ at $\sqrt{s}$ =7000, 0<$p_{\rm T}$<20 GeV/$c$ at $\sqrt{s}$ =8000 and 0<$p_{\rm T}$<20 GeV/$c$ at $\sqrt{s}$ =13000.

$\psi(2S)$ mean transversee momentum vs collision energy. $p_{\rm T}$ integration ranges are 0<$p_{\rm T}$<12 GeV/$c$ at $\sqrt{s}$ =7000, 0<$p_{\rm T}$<12 GeV/$c$ at $\sqrt{s}$ =8000 and 0<$p_{\rm T}$<16 GeV/$c$ at $\sqrt{s}$ =13000.

$\psi(2S)$ mean transversee momentum square vs collision energy. $p_{\rm T}$ integration ranges are 0<$p_{\rm T}$<12 GeV/$c$ at $\sqrt{s}$ =7000, 0<$p_{\rm T}$<12 GeV/$c$ at $\sqrt{s}$ =8000 and 0<$p_{\rm T}$<16 GeV/$c$ at $\sqrt{s}$ =13000.

Differential production cross sections of $J/\psi$ vs collision energy.

Differential production cross sections of $\psi(2S)$ vs collision energy.

Light quark ($uds$) Collins asymmetries obtained by fitting the U/L and U/C double ratios as a function of ($z_1$,$z_2$) for $K\pi$ hadron pairs. In the first column, the $z$ bins and their respective mean values for the hadron ($K$ or $\pi$) in one hemisphere are reported; in the following column, the same variables for the second hadron ($K$ or $\pi$) are shown; in the third column the mean value of $\sin^2\theta_{2}/(1+\cos^2\theta_{2})$ is summarized, calculated in the RF0 frame; in the last two columns the asymmetry results are summarized. The mean values of the quantities reported in the table are calculated by summing the corresponding values for each $K\pi$ pair and dividing by the number of $K\pi$ pairs that fall into each ($z_1$,$z_2$) interval. Note that the $A^{UL}$ and $A^{UC}$ results are strongly correlated since they are obtained by using the same data set.

Light quark ($uds$) Collins asymmetries obtained by fitting the U/L and U/C double ratios as a function of ($z_1$,$z_2$) for pion pairs. In the first column, the $z$ bins and their respective mean values for the pion in one hemisphere are reported; in the following column, the same variables for the second pion are shown; in the third column the mean value of $\sin^2\theta_{th}/(1+\cos^2\theta_{th})$ is summarized, calculated in the RF12 frame; in the last two columns the asymmetry results are summarized. The mean values of the quantities reported in the table are calculated by summing the corresponding values for each $\pi\pi$ pair and dividing by the number of $\pi\pi$ pairs that fall into each ($z_1$,$z_2$) interval. Note that the $A^{UL}$ and $A^{UC}$ results are strongly correlated since they are obtained by using the same data set.

Light quark ($uds$) Collins asymmetries obtained by fitting the U/L and U/C double ratios as a function of ($z_1$,$z_2$) for pion pairs. In the first column, the $z$ bins and their respective mean values for the pion in one hemisphere are reported; in the following column, the same variables for the second pion are shown; in the third column the mean value of $\sin^2\theta_{2}/(1+\cos^2\theta_{2})$ is summarized, calculated in the RF0 frame; in the last two columns the asymmetry results are summarized. The mean values of the quantities reported in the table are calculated by summing the corresponding values for each $\pi\pi$ pair and dividing by the number of $\pi\pi$ pairs that fall into each ($z_1$,$z_2$) interval. Note that the $A^{UL}$ and $A^{UC}$ results are strongly correlated since they are obtained by using the same data set.

The cross section of the process E+ E- --> DEUTBAR X.

The ratio of the cross sections of the processes E+ E- --> DEUTBAR X and E+ E- --> HADRONS.

Integrated cross sections for conventional PI+-, K+- and PBAR/P production. The second (sys) error is the uncertainty due to the model dependence of the extrapolation.

The K$\pi$ pair invariant mass distribution integrated over the $K^{*0}$ $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ =62.4 GeV after mixed-event background subtraction.

The Kπ pair invariant mass distribution for various pT bins (top left) pT = 0.4–0.6 GeV/c in Au+Au collisions at √sNN = 200 GeV after the mixed-event background subtraction.

The Kπ pair invariant mass distribution for various pT bins (top right) pT = 0.6–0.8 GeV/c in Au+Au collisions at √sNN = 62.4 GeV after the mixed-event background subtraction.

The Kπ pair invariant mass distribution for various pT bins (bottom left) pT = 0.8–1.0 GeV/c in Au+Au collisions at √sNN = 200 GeV after the mixed-event background subtraction.

The Kπ pair invariant mass distribution for various pT bins (bottom right) pT = 1.0–1.2 GeV/c in Au+Au collisions at √sNN = 62.4 GeV after the mixed-event background subtraction.

The signal-to-background ratio for $K^{*0}$ measurements as a function of $p_T$ for different collision centrality bins (0-10%, 10-40%, 40-60%, 60-80%) in Au+Au collisions at 200 GeV.

$K^{*0}$ mass as a function of $p_T$ for minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

$K^{*0}$ mass as a function of $p_T$ for minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

$K^{*0}$ mass as a function of $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

$K^{*0}$ mass as a function of $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

$K^{*0}$ width as a function of $p_T$ for minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

$K^{*0}$ width as a function of $p_T$ for minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

$K^{*0}$ width as a function of $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

$K^{*0}$ width as a function of $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

The $K^{*0}$ reconstruction efficiency multiplied by the detector acceptance as a function of $p_T$ in Au+Au (|$\eta$| < 0.8) collisions at 200 GeV for different collision centrality bins (0-20% ,20-40% , 40-60%)

The $K^{*0}$ reconstruction efficiency multiplied by the detector acceptance as a function of $p_T$ in Cu+Cu (|$\eta$| < 1.0) collisions at 200 GeV for different collision centrality bins (0-20% ,20-40% , 40-60%)

Mid-rapidity $K^{*0}$ $p_T$ spectra for various collision centrality bins (0-20%, 20-40%, 40-60%, 60-80%) in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

Mid-rapidity $K^{*0}$ $p_T$ spectra for various collision centrality bins (0-20%, 20-40%, 40-60%) in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

Mid-rapidity $K^{*0}$ $p_T$ spectra for various collision centrality bins (0-20%, 20-40%, 40-60%, 60-80%) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

Mid-rapidity $K^{*0}$ $p_T$ spectra for various collision centrality bins (0-20%, 20-40%, 40-60%) in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity yields dN/dy of $K^{*0}$ as a function of the average number of participating nucleons, $⟨N_{part}⟩$, for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

The mid-rapidity yields dN/dy of $K^{*0}$ as a function of the average number of participating nucleons, $⟨N_{part}⟩$, for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

The mid-rapidity yields dN/dy of $K^{*0}$ as a function of the average number of participating nucleons, $⟨N_{part}⟩$, for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity yields dN/dy of $K^{*0}$ as a function of the average number of participating nucleons, $⟨N_{part}⟩$, for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity $K^{*0}$ $⟨p_T⟩$ as a function $⟨N_{part}⟩$ for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

The mid-rapidity $K^{*0}$ $⟨p_T⟩$ as a function $⟨N_{part}⟩$ for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

The mid-rapidity $K^{*0}$ $⟨p_T⟩$ as a function $⟨N_{part}⟩$ for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity $K^{*0}$ $⟨p_T⟩$ as a function $⟨N_{part}⟩$ for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity $⟨p_T⟩$ of $\pi$, K, p and $K^{*0}$ as a function of $⟨N_{part}⟩$ for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})N(K^-)$ in Au+Au collisions divided by $N(K^{*0})N(K^-)$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$.

Mid-rapidity $N(K^{*0})N(K^-)$ in Cu+Cu collisions divided by $N(K^{*0})N(K^-)$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})N(K^-)$ in d+Au collisions divided by $N(K^{*0})N(K^-)$ ratio in d+Au collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias Au+Au collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias Cu+Cu collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias p+p collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias Au+Au collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias Cu+Cu collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias p+p collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $[N(\phi)/N(K^{*0})]$ in Au+Au collisions divided by $[N(\phi)/N(K^{*0})]$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $[N(\phi)/N(K^{*0})]$ in Cu+Cu collisions divided by $[N(\phi)/N(K^{*0})]$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $[N(\phi)/N(K^{*0})]$ in d+Au collisions divided by $[N(\phi)/N(K^{*0})]$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias Au+Au collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias Cu+Cu collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias p+p collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias Au+Au collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias Cu+Cu collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias p+p collisions as a function of $\sqrt{s_{NN}}$.

The $K^{*0}$ $v_2$ (Run IV) as a function of $p_T$ in minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $K^{*0}$ $v_2$ (Run II) as a function of $p_T$ in minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $K^{*0}$ $R_{CP}$ as a function of $p_T$ in Au+Au collisions at 62.4 and 200 GeV compared to the $R_{CP}$ of $K^0_S$ and $\Lambda$ at 200 GeV.

The $K^{*0}$ $R_{CP}$ as a function of $p_T$ in Au+Au collisions at 62.4 and 200 GeV compared to the $R_{CP}$ of $K^0_S$ and $\Lambda$ at 200 GeV.

The $K^{*0}$ $R_{CP}$ as a function of $p_T$ in Au+Au collisions at 62.4 and 200 GeV compared to the $R_{CP}$ of $K^0_S$ and $\Lambda$ at 200 GeV.

The $K^{*0}$ ~$R_{CP}$~ as a function of $p_T$ in Au+Au collisions at 62.4 and 200 GeV compared to the $R_{CP}$ of $K^0_S$ and $\Lambda$ at 200 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions from Lattice QCD Calculations.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV from Transport Model Calculations.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV from Transport Model Calculations.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations for net-baryon.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations for net-proton (No Decay).

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are STAR data.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are from UrQMD net-proton.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are from AMPT default net-proton.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are from AMPT String Melting net-proton.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are from Therminator net-proton.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are from Hijing net-proton.