Pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in pp collisions at $\sqrt{s} = 5020~\mathrm{GeV}$.

Pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in Pb-p collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in pp collisions at various LHC energies.

$\langle N_{\rm \gamma} \rangle$ measured within $2.3<\eta_{\rm lab}<3.9$ as a function of collision energy in pp collisions.

Centrality-dependent pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$ (Centrality estimator - CL1).

Centrality-dependent pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$ (Centrality estimator - ZNA).

$\langle N_{\rm \gamma} \rangle$ measured within $2.3<\eta_{\rm lab}<3.9$ as a function of $\langle N_{\rm part} \rangle$ for CL1 centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

$\langle N_{\rm \gamma} \rangle$ measured within $2.3<\eta_{\rm lab}<3.9$ as a function of $\langle N_{\rm part} \rangle$ for ZNA (high-$p_{\rm T}$) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

$\langle N_{\rm \gamma} \rangle$ measured within $2.3<\eta_{\rm lab}<3.9$ as a function of $\langle N_{\rm part} \rangle$ for ZNA (Pb-side) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm part} \rangle$ for CL1 centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm part} \rangle$ for ZNA (high-$p_{\rm T}$) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm part} \rangle$ for ZNA (Pb-side) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm \gamma} \rangle$ for CL1 centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm \gamma} \rangle$ for ZNA (high-$p_{\rm T}$) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm \gamma} \rangle$ for ZNA (Pb-side) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

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.

Fraction of bottom hadron decayed electrons in 0-80% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV.

Fraction of bottom hadron decayed electrons in 0-20% and 20-40% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV.

Fraction of bottom hadron decayed electrons in 40-80% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV.

Nuclear modification factor $R_{\rm AA}$ of inclusive heavy flavor decayed electrons in 0-80% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV. The 8% global uncertainty from the $p$+$p$ reference and a 8% gloabl uncertainty from $N_{\rm coll.}$ are not included in the table values.

Nuclear modification factor $R_{\rm AA}$ of bottom- and charm-decay electrons in 0-80% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV. The 8% global uncertainty from the $p$+$p$ reference and a 8% gloabl uncertainty from $N_{\rm coll.}$ are not included in the table values.

Nuclear modification factor $R_{\rm AA}$ ratio of bottom- to charm-decay electrons in 0-80% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV.

Null hypothesis for the nuclear modification factor $R_{\rm AA}$ ratio of bottom- to charm-decay electrons in 0-80% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV.

Nuclear modification factor $R_{\rm CP}$ ratio of bottom- and charm-decay electrons in 0-80% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV.

Null hypotheses for the nuclear modification factor $R_{\rm CP}$ ratio of bottom- and charm-decay electrons in 0-80% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV.

Transverse momentum spectrum of $K^{*0} + \overline{K^{*0}}$ measured at midrapidity ($|y|<0.5$) in inelastic pp collisions at $\sqrt{s}$ = 13 TeV. The normalization uncertainty of $\pm2.6\%$ is excluded.

Transverse momentum spectrum of $\phi$ measured at midrapidity ($|y|<0.5$) in inelastic pp collisions at $\sqrt{s}$ = 13 TeV. The normalization uncertainty of $\pm2.6\%$ is excluded.

Transverse momentum spectrum of $p + \overline{p}$ measured at midrapidity ($|y|<0.5$) in inelastic pp collisions at $\sqrt{s}$ = 13 TeV. The normalization uncertainty of $\pm2.6\%$ is excluded.

Transverse momentum spectrum of $\Lambda + \overline{\Lambda}$ measured at midrapidity ($|y|<0.5$) in inelastic pp collisions at $\sqrt{s}$ = 13 TeV. The normalization uncertainty of $\pm2.6\%$ is excluded.

Transverse momentum spectrum of $\Xi^{-} + \overline{\Xi}^{+}$ measured at midrapidity ($|y|<0.5$) in inelastic pp collisions at $\sqrt{s}$ = 13 TeV. The normalization uncertainty of $\pm2.6\%$ is excluded.

Transverse momentum spectrum of $\Omega^{-} + \overline{\Omega}^{+}$ measured at midrapidity ($|y|<0.5$) in inelastic pp collisions at $\sqrt{s}$ = 13 TeV. The normalization uncertainty of $\pm2.6\%$ is excluded.

Transverse momentum spectrum of $K^{0}_{S}$ measured at midrapidity ($|y|<0.5$) in inelastic pp collisions at $\sqrt{s}$ = 7 TeV. The normalization uncertainty of $^{+7.3}_{-3.5}\%$ is excluded.

Transverse momentum spectrum of $\Lambda + \overline{\Lambda}$ measured at midrapidity ($|y|<0.5$) in inelastic pp collisions at $\sqrt{s}$ = 7 TeV. The normalization uncertainty of $^{+7.3}_{-3.5}\%$ is excluded.

Ratios of the transverse-momentum spectra of $\pi^{+} + \pi^{-}$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

Ratios of the transverse-momentum spectra of $K^{+} + K^{-}$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

Ratios of the transverse-momentum spectra of $K^{0}_{S}$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

Ratios of the transverse-momentum spectra of $K^{*0} + \overline{K^{*0}}$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

Ratios of the transverse-momentum spectra of $\phi$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

Ratios of the transverse-momentum spectra of $p + \overline{p}$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

Ratios of the transverse-momentum spectra of $\Lambda + \overline{\Lambda}$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

Ratios of the transverse-momentum spectra of $\Xi^{-} + \overline{\Xi}^{+}$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

Ratios of the transverse-momentum spectra of $\Omega^{-} + \overline{\Omega}^{+}$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

$(K^{+} + K^{-})/(\pi^{+}+\pi^{-})$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$(K^{*0} + \overline{K^{*0}})/(\pi^{+}+\pi^{-})$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$K^{0}_{S}/(\pi^{+}+\pi^{-})$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$\phi/(\pi^{+}+\pi^{-})$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$\phi/(\pi^{+}+\pi^{-})$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 7 TeV.

Average transverse momentum of $\pi^{+}+\pi^{-}$ as a function of the center-of-mass energy.

Average transverse momentum of $K^{+} + K^{-}$ as a function of the center-of-mass energy.

Average transverse momentum of $K^{0}_{S}$ as a function of the center-of-mass energy.

Average transverse momentum of $K^{*0} + \overline{K^{*0}}$ as a function of the center-of-mass energy.

Average transverse momentum of $\phi$ as a function of the center-of-mass energy.

Average transverse momentum of $p + \overline{p}$ as a function of the center-of-mass energy.

Average transverse momentum of $\Lambda+\overline{\Lambda}$ as a function of the center-of-mass energy.

Average transverse momentum of $\Xi^{-} + \overline{\Xi}^{+}$ as a function of the center-of-mass energy.

Average transverse momentum of $\Omega^{-} + \overline{\Omega}^{+}$ as a function of the center-of-mass energy.

Ratio of $p_{\rm T}$-integrated yields of $K^{+} + K^{-}$ to $\pi^{+} + \pi^{-}$ as a function of the center-of-mass energy.

Ratio of $p_{\rm T}$-integrated yields of $p + \overline{p}$ to $\pi^{+} + \pi^{-}$ as a function of the center-of-mass energy.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda+\overline{\Lambda}$ to $\pi^{+} + \pi^{-}$ as a function of the center-of-mass energy.

Ratio of $p_{\rm T}$-integrated yields of $\Xi^{-} + \overline{\Xi}^{+}$ to $\pi^{+} + \pi^{-}$ as a function of the center-of-mass energy.

Ratio of $p_{\rm T}$-integrated yields of $\Omega^{-} + \overline{\Omega}^{+}$ to $\pi^{+} + \pi^{-}$ as a function of the center-of-mass energy.

Ratio of $p_{\rm T}$-integrated yields of $\phi$ to $K^{+} + K^{-}$ as a function of the center-of-mass energy.

Ratio of $p_{\rm T}$-integrated yields of $K^{*0} + \overline{K^{*0}}$ to $K^{+} + K^{-}$ as a function of the center-of-mass energy.

$(p + \overline{p})$/($\pi^{+}+\pi^{-}$) particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$(\Lambda + \overline{\Lambda})/K^{0}_{S}$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$(\Omega^{-} + \overline{\Omega}^{+})/\phi$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$(\Xi^{-} + \overline{\Xi}^{+})/\phi$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$(\Omega^{-} + \overline{\Omega}^{+})/\phi$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 7 TeV.

$(\Xi^{-} + \overline{\Xi}^{+})/\phi$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 7 TeV.

The results of fits of the two-particle double-ratio correlation function $R_{2} (Q)$ for $\sqrt{s} = 0.9\ TeV$ events with the unlike-charge reference sample for various multiplicity intervals $n_{ch}$ for the exponential parametrization $\Omega^{(E)}$. The uncertainties for parameters $\lambda$ and $R$ are square root from the quadratic sum of statistical and systematic errors. Where only one error is shown for parameters $C_{0}$ and $\epsilon$, this represents the statistical uncertainty only.

The results of fits of the two-particle double-ratio correlation function $R_{2} (Q)$ for $\sqrt{s} = 7\ TeV$ events with the unlike-charge reference sample for various multiplicity intervals $n_{ch}$ for the exponential parametrization $\Omega^{(E)}$. The uncertainties for parameters $\lambda$ and $R$ are square root from the quadratic sum of statistical and systematic errors. The statistical uncertainties at $7\ TeV$ were corrected on the $\sqrt{\chi^{2}/ndf}$ and they more smaller than systematic uncertainties. Where only one error is shown for parameters $C_{0}$ and $\epsilon$, this represents the statistical uncertainty only.

The results of fits of the two-particle double-ratio correlation function $R_{2} (Q)$ for $\sqrt{s} = 7\ TeV\ HM$ events with the unlike-charge reference sample for various multiplicity intervals $n_{ch}$ for the exponential parametrization $\Omega^{(E)}$. The uncertainties for parameters $\lambda$ and $R$ are square root from the quadratic sum of statistical and systematic errors. Where only one error is shown for parameters $C_{0}$ and $\epsilon$, this represents the statistical uncertainty only.

The results of fits of the two-particle double-ratio correlation function $R_{2} (Q)$ for $\sqrt{s} = 0.9\ TeV$ ($n_{ch} \ge 2$) events with the unlike-charge reference sample for various $k_{T}$ intervals for the exponential parametrization $\Omega^{(E)}$. The uncertainties for parameters $\lambda$ and $R$ are square root from the quadratic sum of statistical and systematic errors. The statistical uncertainties were corrected on the $\sqrt{\chi^{2}/ndf}$ when $\chi^{2}/ndf > 1$. Where only one error is shown for parameters $C_{0}$ and $\epsilon$, this represents the statistical uncertainty only.

The results of fits of the two-particle double-ratio correlation function $R_{2} (Q)$ for $\sqrt{s} = 7\ TeV$ ($n_{ch} \ge 2$) events with the unlike-charge reference sample for various $k_{T}$ intervals for the exponential parametrization $\Omega^{(E)}$. The uncertainties for parameters $\lambda$ and $R$ are square root from the quadratic sum of statistical and systematic errors. Where only one error is shown for parameters $C_{0}$ and $\epsilon$ , this represents the statistical uncertainty only. The statistical uncertainties were corrected on the $\sqrt{\chi^{2}/ndf}$ when $\chi^{2}/ndf > 1$. The statistical uncertainties at $7\ TeV$ more smaller than systematic uncertainties.

The results of fits of the two-particle double-ratio correlation function $R_{2} (Q)$ for $\sqrt{s} = 7\ TeV\ HM$ ($n_{ch} \ge 150$) events with the unlike-charge reference sample for various $k_{T}$ intervals for the exponential parametrization $\Omega^{(E)}$. The uncertainties for parameters $\lambda$ and $R$ are square root from the quadratic sum of statistical and systematic errors. The statistical uncertainties were corrected on the $\sqrt{\chi^{2}/ndf}$ when $\chi^{2}/ndf > 1$.Where only one error is shown for parameters $C_{0}$ and $\epsilon$, this represents the statistical uncertainty only.

The results of fits of the two-particle double-ratio correlation function $R_{2} (Q)$ for $\sqrt{s} = 7\ TeV$ events with the unlike-charge reference sample for various $k_{T}$ intervals of different multiplicity regions $n_{ch}$ for the exponential parametrization $\Omega^{(E)}$. The statistical uncertainties were corrected on the $\sqrt{\chi^{2}/ndf}$ when $\chi^{2}/ndf > 1$. The statistical uncertainties smaller than systematic uncertainties. The uncertainties for parameters $\lambda$ and $R$ are square root from the quadratic sum of statistical and systematic errors. Where only one error is shown for parameters $C_{0}$ and $\epsilon$, this represents the statistical uncertainty only.

The $Q$ distribution measured at $0.9\ TeV$ for unlike-sign pairs for $p_{T} > 100\ MeV$ and $|\eta|<2.5$, renormalized to the total number of charged-particle pairs. The error bars represents the statistical uncertainties.

The $Q$ distribution measured at $0.9\ TeV$ for like-sign pairs for $p_{T} > 100\ MeV$ and $|\eta|<2.5$, renormalized to the total number of charged-particle pairs. The error bars represents the statistical uncertainties.

The $Q$ distribution measured at $7\ TeV$ for unlike-sign pairs for $p_{T} > 100\ MeV$ and $|\eta|<2.5$, renormalized to the total number of charged-particle pairs. The error bars represents the statistical uncertainties.

The $Q$ distribution measured at $7\ TeV$ for like-sign pairs for $p_{T} > 100\ MeV$ and $|\eta|<2.5$, renormalized to the total number of charged-particle pairs. The error bars represents the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $0.9\ TeV$ using unlike-charge particle reference sample for different $n_{ch}$ intervals $p_{T} > 100\ MeV$ and $|\eta|<2.5$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $0.9\ TeV$ using unlike-charge particle reference sample for different $k_{T}$ intervals with multiplicity $n_{ch} \ge 2$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $7\ TeV$ using unlike-charge particle reference sample for different multiplicity intervals $n_{ch}$ with $p_{T} > 100\ MeV$ and $|\eta| < 2.5$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $7\ TeV$ using unlike-charge particle reference sample for different $k_{T}$ intervals with multiplicity $n_{ch} \ge 2$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $7\ TeV\ HM$ using unlike-charge particle reference sample for different $n_{ch}$ intervals $p_{T} > 100\ MeV$ and $|\eta|<2.5$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured with $7\ TeV\ HM$ data using unlike-charge particle reference sample for different $k_{T}$ intervals with multiplicity $n_{ch} \ge 150$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $7\ TeV$ using unlike-charge particle reference sample for different $k_{T}$ intervals within multiplicity interval $n_{ch} = 2-9$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $7\ TeV$ using unlike-charge particle reference sample for different $k_{T}$ intervals within multiplicity interval $n_{ch} = 10-24$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $7\ TeV$ using unlike-charge particle reference sample for different $k_{T}$ intervals within multiplicity interval $n_{ch} = 25-80$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $7\ TeV$ using unlike-charge particle reference sample for different $k_{T}$ intervals within multiplicity interval $n_{ch} = 81-125$. The error bars represents only the statistical uncertainties.

High transverse-momentum spectra ($p_{T} > 2.5$ GeV/c) of charged pions, protons, and antiprotons for the rapidity regions $|y| < 0.5$ (solid symbols) and $0.5 < |y| < 1.0$ (open symbols) for $d+Au$ collisions and various event centrality classes at $\sqrt{s_{NN}}=200$ GeV.

High transverse-momentum spectra ($p_{T} > 2.5$ GeV/c) of charged pions, protons, and antiprotons for the rapidity regions $|y| < 0.5$ (solid symbols) and $0.5 < |y| < 1.0$ (open symbols) for $d+Au$ collisions and various event centrality classes at $\sqrt{s_{NN}}=200$ GeV.

High transverse-momentum spectra ($p_{T} > 2.5$ GeV/c) of charged pions, protons, and antiprotons for the rapidity regions $|y| < 0.5$ (solid symbols) and $0.5 < |y| < 1.0$ (open symbols) for $d+Au$ collisions and various event centrality classes at $\sqrt{s_{NN}}=200$ GeV.

High transverse-momentum spectra ($p_{T} > 2.5$ GeV/c) of charged pions, protons, and antiprotons for the rapidity regions $|y| < 0.5$ (solid symbols) and $0.5 < |y| < 1.0$ (open symbols) for $d+Au$ collisions and various event centrality classes at $\sqrt{s_{NN}}=200$ GeV.

High transverse-momentum spectra ($p_{T} > 2.5$ GeV/c) of charged pions, protons, and antiprotons for the rapidity regions $|y| < 0.5$ (solid symbols) and $0.5 < |y| < 1.0$ (open symbols) for $d+Au$ collisions and various event centrality classes at $\sqrt{s_{NN}}=200$ GeV.

High transverse-momentum rapidity asymmetry factor ($Y_{\scriptsize{\mbox{Asym}}}$) for $\pi^{+}+\pi^{-}$ and $p+\bar{p}$ for $|y| < 0.5$ and $0.5 < |y| < 1.0$ for minimum bias $d+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV. For comparison the inclusive charged hadron results from STAR [5] are also shown by the curves.

Centrality dependence of high transverse-momentum rapidity asymmetry factor ($Y_{\scriptsize{\mbox{Asym}}}$) for $\pi^{+}+\pi^{-}$ and $p+\bar{p}$ for $|y| < 0.5$ (left panels) and $0.5 < |y| < 1.0$ (right panels) for $d+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV.

Centrality dependence of high transverse-momentum rapidity asymmetry factor ($Y_{\scriptsize{\mbox{Asym}}}$) for $\pi^{+}+\pi^{-}$ and $p+\bar{p}$ for $|y| < 0.5$ (left panels) and $0.5 < |y| < 1.0$ (right panels) for $d+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV.

Variation of nuclear modification factor ($R_{CP}$) for $\pi^{+}+\pi^{-}$ and $p+\bar{p}$ with rapidity for $p_{T} > 2.5$ GeV/c for $d+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV. Also shown for comparison are the $R_{CP}$ values for inclusive charged hadrons as a function of pseudorapidity [5]. The errors shown as boxes are the systematic errors. The error due to number of binary collisions is $\sim 14\%$ and is not shown in the figure.

Variation of $\pi^{-} / \pi^{+}$ and $\bar{p} / p$ with rapidity for $2.5 < p_{T} < 10$ GeV/c for peripheral ($40-100\%$) and n-tag events for $d+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV. Also shown for comparison are the $\bar{p} / p$ for minimum bias $p+p$ collisions. The boxes are the systematic errors. The ratios for n-tag events are shifted by 0.05 units in rapidity and those for $p+p$ collisions by -0.05 units in rapidity for clarity of presentation.

Variation of $\pi^{-} / \pi^{+}$ and $\bar{p} / p$ with rapidity for $2.5 < p_{T} < 10$ GeV/c for peripheral ($40-100\%$) and n-tag events for $d+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV. Also shown for comparison are the $\bar{p} / p$ for minimum bias $p+p$ collisions. The boxes are the systematic errors. The ratios for n-tag events are shifted by 0.05 units in rapidity and those for $p+p$ collisions by -0.05 units in rapidity for clarity of presentation.

Variation of double ratio $(\pi^{-} / \pi^{+})_{\scriptsize{\mbox{n-tag}}}/(\pi^{-} / \pi^{+})_{40-100\%}$ and $(\bar{p} / p)_{\scriptsize{\mbox{n-tag}}}/(\bar{p} / p)_{40-100\%}$ with rapidity for $2.5 < p_{T} < 10$ GeV/c for $d+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV. The boxes are the systematic errors.

K0sK0s correlation function

K0sK0s theoretical correlation function

K0sK0s correlation function with fit functions

K0sK0s radius versus pair purity

K0sK0s lambda versus pair purity

K0sK0s radius versus Mean KT

K0sK0s radius versus MT

FIG. 2: a) left side: The correlation data for the ratio of the histograms of same-event-pairs to mixed-event-pairs for unlike-sign charged pairs, shown in a two-dimensional (2-D) perspective plot $\Delta\phi$ - $\Delta\eta$. The plot was normalized to a mean of 1. b) right side: The similar correlation data for like-sign charge pairs.

FIG. 3: a) left side: Folded correlation data for unlike-sign charge pairs on $\Delta\phi$ and $\Delta\eta$ based on demonstrated symmetries in the data. This increases the statistics in a typical bin by a factor of four. Henceforth we will be dealing with folded data only. b) right side: Similar folded data for like-sign charge pairs.

FIG. 3: a) left side: Folded correlation data for unlike-sign charge pairs on $\Delta\phi$ and $\Delta\eta$ based on demonstrated symmetries in the data. This increases the statistics in a typical bin by a factor of four. Henceforth we will be dealing with folded data only. b) right side: Similar folded data for like-sign charge pairs.

FIG. 4: “(Color online)”a) left side: The 2 dimensional (2-D) approximately gaussian signal shape equation (4) from the fit to the normalized and folded unlike-sign charge pairs signal data. b) right side: The unlike-sign charge pairs signal data corresponding to the adjoining signal function on the left side.

FIG. 4: “(Color online)”a) left side: The 2 dimensional (2-D) approximately gaussian signal shape equation (4) from the fit to the normalized and folded unlike-sign charge pairs signal data. b) right side: The unlike-sign charge pairs signal data corresponding to the adjoining signal function on the left side.

FIG. 5 :“(Color online)”a) left side: A perspective plot of the fit to the 2-D like-sign charge pairs signal shape of equation (5). b) right side: Like-sign charge pairs signal data corresponding to the adjacent signal fit on the left side.

FIG. 5: “(Color online)”a) left side: A perspective plot of the fit to the 2-D like-sign charge pairs signal shape of equation (5). b) right side: Like-sign charge pairs signal data corresponding to the adjacent signal fit on the left side.

FIG. 6: “(Color online)”a) left side: The charge dependent (CD) signal shape is a 2 dimensional (2-D) approximate gaussian which is symmetrical. The average gaussian width is 27.5$^{\circ}$ $\pm$ 3.2$^{\circ}$ (0.48 $\pm$ 0.056 rad) both in $\Delta\phi$ and $\Delta\eta$. Thus we are observing a greater probability for unlike-sign charge pairs of particles than for like-sign charge pairs emitted randomly on 2-D $\eta\phi$ directions. b) right side: The CD signal data corresponding to the adjacent signal fit on the left.

FIG. 6: “(Color online)”a) left side: The charge dependent (CD) signal shape is a 2 dimensional (2-D) approximate gaussian which is symmetrical. The average gaussian width is 27.5$^{\circ}$ $\pm$ 3.2$^{\circ}$ (0.48 $\pm$ 0.056 rad) both in $\Delta\phi$ and $\Delta\eta$. Thus we are observing a greater probability for unlike-sign charge pairs of particles than for like-sign charge pairs emitted randomly on 2-D $\eta\phi$ directions. b) right side: The CD signal data corresponding to the adjacent signal fit on the left.

FIG. 7: “(Color online)”a) left side: The fit to the charge independent (CI) signal shape which is the sum of the fits of like-sign plus unlike-sign charge pairs signals. The 2-D gaussian equivalent RMS signal has a $\Delta\phi$ width of about 32$^{\circ}$ (0.52 rad), and $\Delta\eta$ width of about 66$^{\circ}$ (1.15 rad) which is about double the $\Delta\phi$ width. b) right side: The CI signal data corresponding to the adjacent fit on the left.

FIG. 7: “(Color online)”a) left side: The fit to the charge independent (CI) signal shape which is the sum of the fits of like-sign plus unlike-sign charge pairs signals. The 2-D gaussian equivalent RMS signal has a $\Delta\phi$ width of about 32$^{\circ}$ (0.52 rad), and $\Delta\eta$ width of about 66$^{\circ}$ (1.15 rad) which is about double the $\Delta\phi$ width. b) right side: The CI signal data corresponding to the adjacent fit on the left.

FIG. 8: a) left side: The charge independent (CI) signal data with the lower range of elliptic flow amplitude and the best $\chi^2$ (the same as Figure 7b) plotted as a 2-D perspective plot on $\Delta\phi$ vs. $\Delta\eta$. b) right side: The charge independent (CI) signal data with our maximum elliptic flow amplitude used in our analysis ($\chi^2$ was worse by 1$\sigma$).

FIG. 8: a) left side: The charge independent (CI) signal data with the lower range of elliptic flow amplitude and the best $\chi^2$ (the same as Figure 7b) plotted as a 2-D perspective plot on $\Delta\phi$ vs. $\Delta\eta$. b) right side: The charge independent (CI) signal data with our maximum elliptic flow amplitude used in our analysis ($\chi^2$ was worse by 1$\sigma$).

FIG. 7. Left: Parameters $n$ and $n_{h} / \tilde{n}_{s}$ vs $\hat{n}_{\mathrm{ch}}$ from free $\chi^{2}$ fits (solid symbols) in Table I. The open symbols represent similar free fits but with $\bar{y}_{t}$ held fixed at $2.65$ (see text for discussion). The bands represent the corresponding two-component fixed parametrizations with their stated errors also given in Table $I$. Right: Ratio $n_{h}\left(\hat{n}_{\mathrm{ch}}\right) / n_{s}\left(\hat{n}_{\mathrm{ch}}\right)$ derived from integrals in Fig. 4 (left panel) (open squares) and from Gaussian amplitudes in Fig. 5 (left panel) (solid dots). The dashed and dotted lines have slopes $0.105$ and $0.095$, respectively. The solid curve is described in the text. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

FIG. 8. $\left\langle p_{t}\right\rangle\left(\hat{n}_{\mathrm{ch}}\right)$ derived from the two-component $H_{0}$ Gaussian amplitudes (solid dots), from the running integrals (open squares), and from power-law fits to STAR and UA1 data (open FIG. 8. $\left\langle p_{t}\right\rangle\left(\hat{n}_{\mathrm{ch}}\right)$ derived from the two-component $H_{0}$ Gaussian amplitudes (solid dots), from the running integrals (open squares), and from power-law fits to STAR and UA1 data (open circles, triangles). The solid and dotted lines correspond to Eq. (6) with $\alpha=0.0095$ and $0.015$, respectively.circles, triangles). The solid and dotted lines correspond to Eq. (6) with $\alpha=0.0095$ and $0.015$, respectively. NOTE 1: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty. NOTE 2: The "ext" uncertainty corresponds to the systematic uncertainty of the extrapolation of the particle yield below $p_T < 0.2 GeV/c$ (hatched region).

(a) The corrected integrated yield $dN/dy$ at mid-rapidity for $\overline{\Xi}^{+}$ , $\overline{\Lambda}$ and $\overline{p}$ divided by $N_{part}$, normalized to the most peripheral centrality interval (60 − 80%), plotted as a function of $N_{part}$. The gray, black and dashed bands represent the errors on the normalization to the most peripheral bin for the $\overline{\Xi}^{+}$, $\overline{\Lambda}$ and $\overline{p}$. Other errors shown are statistical only. (b) $\gamma_{s}$ as a function of $N_{part}$ calculated from thermal model fits to the measured particle yields ($\pi$,K,p [24], $\Lambda$, $\Xi$, $\Omega$ and their anti-particles) at 200 GeV. Values for $e^{+}+e^{−}$ and $p+\overline{p}$ collisions at $\sqrt{s_{NN}}$ = 91 and 200 GeV respectively and for Pb+Pb SPS collisions at $\sqrt{s_{NN}}$ = 17.2 GeV are shown for comparison [26, 27, 28].

(a) The corrected integrated yield $dN/dy$ at mid-rapidity for $\overline{\Xi}^{+}$ , $\overline{\Lambda}$ and $\overline{p}$ divided by $N_{part}$, normalized to the most peripheral centrality interval (60 − 80%), plotted as a function of $N_{part}$. The gray, black and dashed bands represent the errors on the normalization to the most peripheral bin for the $\overline{\Xi}^{+}$, $\overline{\Lambda}$ and $\overline{p}$. Other errors shown are statistical only. (b) $\gamma_{s}$ as a function of $N_{part}$ calculated from thermal model fits to the measured particle yields ($\pi$,K,p [24], $\Lambda$, $\Xi$, $\Omega$ and their anti-particles) at 200 GeV. Values for $e^{+}+e^{−}$ and $p+\overline{p}$ collisions at $\sqrt{s_{NN}}$ = 91 and 200 GeV respectively and for Pb+Pb SPS collisions at $\sqrt{s_{NN}}$ = 17.2 GeV are shown for comparison [26, 27, 28].

$R_{CP}$ for $\Xi^{-}+\overline{\Xi}^{+}$ and $\Omega^{-}+\overline{\Omega}^{+}$ at mid-rapidity (centrality interval : 0 − 5% vs. 40 − 60%). A dashed line for charged hadrons and gray band for $\Lambda^{-}+\overline{\Lambda}$ are shown as comparison. The gray rectangles represent participant and binary scalings.