$R_{AA}$ as a function of $p_{T}$ from $0-5\%$ and $60-80\%$ central Au+Au events for $\Lambda$. Stat. and syst. errors added in quadrature.The band at unity shows the systematic uncertainty on $⟨N_{bin}⟩$. The dashed line below unity shows the expected value of RAA should the yields scale with $⟨N_{part}⟩$, and the band around it shows the systematic uncertainty on ⟨Npart ⟩. $N_{bin}$ Scaling 1 $\pm$ 0.16 and $N_{part}$ Scaling 0.168 $\pm$ 0.02 .

$R_{AA}$ as a function of $p_{T}$ from $0-5\%$ and $60-80\%$ central Au+Au events for $\Xi^{-} + \bar{\Xi}^{+}$. Stat. and syst. errors added in quadrature. The band at unity shows the systematic uncertainty on $⟨N_{bin}⟩$. The dashed line below unity shows the expected value of RAA should the yields scale with $⟨N_{part}⟩$, and the band around it shows the systematic uncertainty on ⟨Npart ⟩. $N_{bin}$ Scaling 1 $\pm$ 0.16 and $N_{part}$ Scaling 0.168 $\pm$ 0.02 .

$R_{AA}$ as a function of $p_{T}$ from $0-5\%$ and $60-80\%$ central Au+Au events for inclusive $p+\bar{p}$. Stat. and syst. errors added in quadrature. The band at unity shows the systematic uncertainty on $⟨N_{bin}⟩$. The dashed line below unity shows the expected value of RAA should the yields scale with $⟨N_{part}⟩$, and the band around it shows the systematic uncertainty on ⟨Npart ⟩. $N_{bin}$ Scaling 1 $\pm$ 0.16 and $N_{part}$ Scaling 0.168 $\pm$ 0.02 .

$R_{AA}$ from HIJING as a function of $p_{T}$ from $0-5\%$ central Au+Au events for inclusive $p+\bar{p}$.

$R_{AA}$ from HIJING as a function of $p_{T}$ from $0-5\%$ central Au+Au events for $\Lambda$.

$R_{AA}$ from HIJING as a function of $p_{T}$ from $0-5\%$ central Au+Au events for $\Xi^{-}+\bar{\Xi}^{+}$.

$R_{AA}$ from EPOS as a function of $p_{T}$ from $0-5\%$ central Au+Au events for inclusive $p+\bar{p}$.

$R_{AA}$ from EPOS as a function of $p_{T}$ from $0-5\%$ central Au+Au events for $\Lambda$.

$R_{AA}$ from EPOS as a function of $p_{T}$ from $0-5\%$ central Au+Au events for $\Xi^{-}+\bar{\Xi}^{+}$.

$R_{AA}$ from HIJING as a function of $p_{T}$ from $60-80\%$ central Au+Au events for inclusive $p+\bar{p}$.

$R_{AA}$ from HIJING as a function of $p_{T}$ from $60-80\%$ central Au+Au events for $\Lambda$.

$R_{AA}$ from HIJING as a function of $p_{T}$ from $60-80\%$ central Au+Au events for $\Xi^{-}+\bar{\Xi}^{+}$.

$R_{AA}$ from EPOS as a function of $p_{T}$ from $60-80\%$ central Au+Au events for inclusive $p+\bar{p}$.

$R_{AA}$ from EPOS as a function of $p_{T}$ from $60-80\%$ central Au+Au events for $\Lambda$.

$R_{AA}$ from EPOS as a function of $p_{T}$ from $60-80\%$ central Au+Au events for $\Xi^{-}+\bar{\Xi}^{+}$.

Measured and calculated asymmetries vs. di-jet pseudorapidity sum for $+\hat{z}$ (left) and $−\hat{z}$ (right) beams. (a,b): Fraction of the calculated di-jet cross section with a quark (gluon) from the $+\hat{z}$ $(−\hat{z})$ beam. (c,d): Unweighted asymmetries compared with pQCD calculations [20] (histograms) for two models of quark Sivers functions fitted to SIDIS results [8]. (e,f): Asymmetries for $|\sin\zeta|$-weighted yields, compared with calculations [20, 21] based on twist-3 quark-gluon correlations. Vertical (horizontal) bars on the data indicate statistical uncertainties (bin widths). The systematic error bands exclude a $\pm12\%$ beam polarization normalization uncertainty.

(Color online) Global polarization of $\Lambda$–hyperons as a function of centrality given as fraction of the total inelastic hadronic cross section. Filled circles show the results for Au+Au collisions at $\sqrt{s_{NN}}$=200 GeV (centrality region 20-70%) and open squares indicate the results for Au+Au collisions at $\sqrt{s_{NN}}$=62.4 GeV (centrality region 0-80%). Only statistical uncertainties are shown.

(Color online) Global polarization of $\overline{\Lambda}$–hyperons as a function of $\overline{\Lambda}$ transverse momentum $p^{\overline{\Lambda}}_{t}$. Filled circles show the results for Au+Au collisions at $\sqrt{s_{NN}}$=200 GeV (centrality region 20-70%) and open squares indicate the results for Au+Au collisions at $\sqrt{s_{NN}}$=62.4 GeV (centrality region 0-80%). Only statistical uncertainties are shown.

(Color online) Global polarization of $\overline{\Lambda}$–hyperons as a function of $\overline{\Lambda}$ pseudorapidity $\eta^{\overline{\Lambda}}$. Filled circles show the results for Au+Au collisions at $\sqrt{s_{NN}}$=200 GeV (centrality region 20-70%). A constant line fit to these data points yields $P_{\Lambda}=(1.8\pm 10.8)\times 10^{-3}$ with $\chi^{2}/ndf=5.5/10$. Open squares show the results for Au+Au collisions as $\sqrt{s_{NN}}$=62.4 GeV (centrality region 0-80%). A constant line fit gives $P_{\Lambda}=(-17.6\pm 11.1)\times 10^{-3}$ with $\chi^{2}/ndf=8.0/10$. Only statistical uncertainties are shown.

(Color online) Global polarization of $\overline{\Lambda}$–hyperons as a function of as a function of centrality. Filled circles show the results for Au+Au collisions at $\sqrt{s_{NN}}$=200 GeV (centrality region 20-70%) and open squares indicate the result for Au+Au collisions at $\sqrt{s_{NN}}$=62.4 GeV (centrality region 0-80%). Only statistical uncertainties are shown.

(Color online) Acceptance correction function $A_{0}(p^{H}_{t}, \eta^{H})$ defined in (10) as a function of $\Lambda$ (filled circles) and $\overline{\Lambda}$ (open squares) transverse momentum (top) and pseudorapidity (bottom). The deviation of this function from unity affects the overall scale of the measured global polarization according to Eq. 9. See the text for details and discussions on $A_{0}$ $p^{H}_{t}$ and $\eta^{H}$ dependence.

(Color online) Acceptance correction function $A_{0}(p^{H}_{t}, \eta^{H})$ defined in (10) as a function of $\Lambda$ (filled circles) and $\overline{\Lambda}$ (open squares) transverse momentum (top) and pseudorapidity (bottom). The deviation of this function from unity affects the overall scale of the measured global polarization according to Eq. 9. See the text for details and discussions on $A_{0}$ $p^{H}_{t}$ and $\eta^{H}$ dependence.

(Color online) Function $A_{2}(p^{H}_{t}, \eta^{H})$ defined in (11) as a function of $\Lambda$ (filled circles) and $\overline{\Lambda}$ (open squares) transverse momentum (top) and pseudorapidity (bottom). The deviation of this function from zero defines the contribution to the observable (3) from $P^{(2)}_{H} (p^{H}_{t}, \eta^{H})$ in the expansion (6).

(Color online) Function $A_{2}(p^{H}_{t}, \eta^{H})$ defined in (11) as a function of $\Lambda$ (filled circles) and $\overline{\Lambda}$ (open squares) transverse momentum (top) and pseudorapidity (bottom). The deviation of this function from zero defines the contribution to the observable (3) from $P^{(2)}_{H} (p^{H}_{t}, \eta^{H})$ in the expansion (6).

(Color online) $\Lambda$ and $\overline{\Lambda}$ spectra on the deuteron and on the gold side in $d$ + Au minimum bias collisions. The data points on the gold side are multiplied by 2 for better visibility. The statistical errors are smaller than the points marking the measurements. The curves show a fit with a Boltzmann function in transverse mass to the data points.

(Color online) (a) Comparison of the measured $\overline{\Lambda}$ yield with model calculations. Statistical errors are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic er- rors including the overall normalization uncertainty. In both panels the target and projectile beam rapidities are indicated by arrows.

(Color online) (b) Comparison of the net $\Lambda$ yield with model calculations. Statistical errors are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic er- rors including the overall normalization uncertainty. In both panels the target and projectile beam rapidities are indicated by arrows.

(Color online) Comparison of $\overline{\Lambda}$ and net $\Lambda$ yields to model calculations for all three centrality classes. Statistical errors are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic errors. Beam rapidity is indicated by arrows.

(Color online) Comparison of $\overline{\Lambda}$ and net $\Lambda$ yields to model calculations for all three centrality classes. Statistical errors are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic errors. Beam rapidity is indicated by arrows.

(Color online) Comparison of $\overline{\Lambda}$ and net $\Lambda$ yields to model calculations for all three centrality classes. Statistical errors are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic errors. Beam rapidity is indicated by arrows.

(Color online) Comparison of $\overline{\Lambda}$ and net $\Lambda$ yields to model calculations for all three centrality classes. Statistical errors are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic errors. Beam rapidity is indicated by arrows.

(Color online) Comparison of $\overline{\Lambda}$ and net $\Lambda$ yields to model calculations for all three centrality classes. Statistical errors are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic errors. Beam rapidity is indicated by arrows.

(Color online) Comparison of $\overline{\Lambda}$ and net $\Lambda$ yields to model calculations for all three centrality classes. Statistical errors are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic errors. Beam rapidity is indicated by arrows.

(Color online) Minimum bias $\overline{\Lambda}/\Lambda$ ratio compared to model calculations. On the deuteron side HIJING/BB shows the best agreement with the results, while on the Au side only AMPT and EPOS give a satisfactory description of the data.

$\overline{\Lambda}/\Lambda$ ratio and net $\Lambda$ and $\overline{\Lambda}$ yields as a function of collision centrality on both the deuteron (left) and the gold side (right). On the deuteron side, centrality is expressed by the number of collisions per deuteron participant, while on the gold side the number of Au participants is chosen. Only statistical errors are shown. The increase in baryon number transport with centrality, shown by the net $\Lambda$ yield, is matched by the increase of $\overline{\Lambda}$-$\Lambda$ pair production, thus keeping the $\overline{\Lambda}/\Lambda$ ratio constant over a wide centrality range.

$\overline{\Lambda}/\Lambda$ ratio and net $\Lambda$ and $\overline{\Lambda}$ yields as a function of collision centrality on both the deuteron (left) and the gold side (right). On the deuteron side, centrality is expressed by the number of collisions per deuteron participant, while on the gold side the number of Au participants is chosen. Only statistical errors are shown. The increase in baryon number transport with centrality, shown by the net $\Lambda$ yield, is matched by the increase of $\overline{\Lambda}$-$\Lambda$ pair production, thus keeping the $\overline{\Lambda}/\Lambda$ ratio constant over a wide centrality range.

$\overline{\Lambda}/\Lambda$ ratio and net $\Lambda$ and $\overline{\Lambda}$ yields as a function of collision centrality on both the deuteron (left) and the gold side (right). On the deuteron side, centrality is expressed by the number of collisions per deuteron participant, while on the gold side the number of Au participants is chosen. Only statistical errors are shown. The increase in baryon number transport with centrality, shown by the net $\Lambda$ yield, is matched by the increase of $\overline{\Lambda}$-$\Lambda$ pair production, thus keeping the $\overline{\Lambda}/\Lambda$ ratio constant over a wide centrality range.

$\overline{\Lambda}/\Lambda$ ratio and net $\Lambda$ and $\overline{\Lambda}$ yields as a function of collision centrality on both the deuteron (left) and the gold side (right). On the deuteron side, centrality is expressed by the number of collisions per deuteron participant, while on the gold side the number of Au participants is chosen. Only statistical errors are shown. The increase in baryon number transport with centrality, shown by the net $\Lambda$ yield, is matched by the increase of $\overline{\Lambda}$-$\Lambda$ pair production, thus keeping the $\overline{\Lambda}/\Lambda$ ratio constant over a wide centrality range.

(Color online) Net $\Lambda dN/dy$ for central, mid-central and peripheral events on both the deuteron and the Au side of the collision. The data are compared to calculations of the distribution of net baryons obtained with the Multichain model [16] with $\alpha$ = 2.9, scaled by 0.1 to account for the conversion from net baryons to net $\Lambda$. An overall scale uncertainty of 10% on the model curves from this conversion is not shown. See text for details.

(a) The solid curve represents the fraction of $\Delta G$ for GRSV-std that has $x > x_{min}$ for scale $Q^{2}$ = 100 $GeV^{2}/c^{2}$. The histograms show the $x_{gluon}$ sampled in the lowest and highest jet $p_{T}$ bins. (b) Confidence levels for several gluon polarization distributions, characterized by their $\Delta G$ at an input scale of 0.4 $GeV^{2}/c^{2}$ [4, 14, 23].

(a) The solid curve represents the fraction of $\Delta G$ for GRSV-std that has $x > x_{min}$ for scale $Q^{2}$ = 100 $GeV^{2}/c^{2}$. The histograms show the $x_{gluon}$ sampled in the lowest and highest jet $p_{T}$ bins. (b) Confidence levels for several gluon polarization distributions, characterized by their $\Delta G$ at an input scale of 0.4 $GeV^{2}/c^{2}$ [4, 14, 23].

The invariant mass distribution for the coherently produced $\rho^{0}$ candidates obtained from the topology sample with the cut on the $\rho^{0}$ transverse momentum $p_{T}$ < 150 MeV/c. The hatched area is the contribution from the combinatorial background. The solid line corresponds to Eq. 3 which encompasses the Breit-Wigner (dashed), the mass independent contribution from the direct $\pi^{+}\pi^{-}$ production (dash-dotted), and the interference term(dotted).

The ratio |B/A| as a function of $y_{\rho^{0}}$ for the minimum bias data, obtained by fitting Eq.3 to the invariant mass distributions in bins of $y_{\rho^{0}}$.

Coherent $\rho^{0}$ production cross-section for the minimum bias data set as a function of $y_{\rho^{0}}$ (black dots) overlaid by the $d\sigma/dy$ distribution predicted by the KN model [22] (solid line).

$\rho^{0}$ production cross-section as a function of the momentum transfer squared $t$, together with the fit of Eq. 5. The fit parameters are shown in Table I.

Comparison of theoretical predictions to the measured differential cross-section for coherent $\rho^{0}$ production. The statistical errors are shown by the solid vertical line at each data point. The sum of the statistical and systematic error bars is shown by the grey band.

Projections of the two dimensional efficiency corrected $\Phi_{h}$ vs $cos(\Theta_{h})$ distributions obtained with the minimum bias data set. The solid line shows the result of the two-dimensional fit to the data with Eq. 6 and the coefficients given in Tab. III.

Projections of the two dimensional efficiency corrected $\Phi_{h}$ vs $cos(\Theta_{h})$ distributions obtained with the minimum bias data set. The solid line shows the result of the two-dimensional fit to the data with Eq. 6 and the coefficients given in Tab. III.

Analyzing powers A_N(-x_F) in x_F bins at < eta > =3.3 and x_F < 0.

Analyzing powers versus pi^0 transverse momentum (p_T) for events with scaled pi^0 longitudinal momentum x_F > 0.4.

Analyzing powers versus pi^0 transverse momentum (p_T) for events with scaled pi^0 longitudinal momentum x_F < -0.4.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.25 < x_F < 0.30 bin.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.30 < x_F < 0.35 bin.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.35 < x_F < 0.40 bin.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.40 < x_F < 0.47 bin.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.47 < x_F < 0.56 bin.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.56 < x_F < 0.80 bin.

$v_{2}$ as a function of $p_{T}$ for charged hadrons with $|\eta| < 1.0$ in 10–40$%$ $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV, from the Lee-Yang Zero method (solid circles) with Sum Generating Function.

$v_{2}$ as a function of $p_{T}$ for charged hadrons with $|\eta| < 1.0$ in 0–80$%$ $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV, from the Event Plane method.

$p_{T}$ integrated charged hadron $v_{2}$ in the TPC as a function of geometrical cross section. Shown are the Event Plane method ($v_{2}${EP}) (open circles). All data are from $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV. For the TPC, $|\eta|< 1.0$ was used.

$p_{T}$ integrated charged hadron $v_{2}$ in the TPC as a function of geometrical cross section. Shown are the 4 Cumulant method ($v_{2}${4}) (solid squares). All data are from $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV. For the TPC, $|\eta|< 1.0$ was used.

$p_{T}$ integrated charged hadron $v_{2}$ in the TPC as a function of geometrical cross section. Shown are the Lee-Yang Zero method with Sum Generating Function ($v_{2}$) (solid circles). All data are from $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV. For the TPC, $|\eta|< 1.3$ was used.

$p_{T}$ integrated charged hadron $v_{2}$ in the TPC as a function of geometrical cross section. Shown are the Lee-Yang Zero method with Product Generating Function ($v_{2}$) (open stars). All data are from $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV. For the TPC, $|\eta|< 1.3$ was used.

$v_{2}$ as a function of $p_{T}$ for $10–40\%$ centrality using $\eta$-subevent method (open crosses), are shown in (a) for $\Lambda$. All data are from $\sqrt{s_{NN}} = 200$ GeV $Au+Au$ collisions.

$v_{2}$ as a function of $p_{T}$ for $10–40\%$ centrality using Event Plane method (open circles), are shown in (a) for $\Lambda$. All data are from $\sqrt{s_{NN}} = 200$ GeV $Au+Au$ collisions.

$v_{2}$ as a function of $p_{T}$ for $10–40\%$ centrality using $\eta$-subevent method (open crosses), are shown in (b) for $K_{S}^{0}$. All data are from $\sqrt{s_{NN}} = 200$ GeV $Au+Au$ collisions.

$v_{2}$ as a function of $p_{T}$ for $10–40\%$ centrality using Event Plane method (open circles), are shown in (b) for $K_{S}^{0}$. All data are from $\sqrt{s_{NN}} = 200$ GeV $Au+Au$ collisions.

$v_{2}$ as a function of $p_{T}$ for $10–40\%$ centrality using Lee-Yang Zero method with Sum Generating Function (solid circles), are shown in (a) and (b) for $\Lambda$ and $K_{S}^{0}$ respectively. All data are from $\sqrt{s_{NN}} = 200$ GeV $Au+Au$ collisions.

$v_{2}$ of $K_{S}^{0}$ (open circles) as a function of $p_{T}$ for (a) $0–80\%$, (b) $40–80\%$, (c) $10–40\%$ and (d) $0–10\%$ centralities in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Lambda+\bar{\Lambda}$ (open squares) as a function of $p_{T}$ for (a) $0–80\%$, (b) $40–80\%$, (c) $10–40\%$ and (d) $0–10\%$ centralities in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Xi^{-}+\bar{\Xi}^{+}$ (filled triangles) as a function of $p_{T}$ for (a) $0–80\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Xi^{-}+\bar{\Xi}^{+}$ (filled triangles) as a function of $p_{T}$ for (b) $40–80\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Xi^{-}+\bar{\Xi}^{+}$ (filled triangles) as a function of $p_{T}$ for (c) $10–40\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Xi^{-}+\bar{\Xi}^{+}$ (filled triangles) as a function of $p_{T}$ for (d) $0–10\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Omega^{-}+\bar{\Omega}^{+}$ (filled circles) as a function of $p_{T}$ for (a) $0–80\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Omega^{-}+\bar{\Omega}^{+}$ (filled circles) as a function of $p_{T}$ for (d) $0–10\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $p+\bar{p}$ (filled squares) and $\pi^{+}+\pi^{-}$ (filled stars) as a function of $p_{T}$ for (a) $0–80\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

Fig12d inset $\Lambda$ (open squares) expansion at low $p_{T}$

Fig12d inset for $K_{S}^{0}$ (filled circles) expansion at low $p_{T}$

Fig13 Centrality dependence of $v_{2}$ versus number of participants $N_{part}$ for charged hadrons (crosses) in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV. For identified particle's value, $v_{2}$ in Tab. II scaled by participant eccentricity in Tab. I versus number of participant in Tab. I.

Fig14 $p_{T}$ dependence of $v_{2}$ of $K_{S}^{0}$ at 200 GeV.

Fig14 $p_{T}$ dependence of $v_{2}$ of $\Lambda$ at 200 GeV.

Delta++ reconstruction efficiency times detector acceptance as a function of the invariant mass for minimum bias d+Au collisions. The error shown is due to the available statistics in the simulation.

Delta++ reconstruction efficiency times detector acceptance as a function of the invariant mass for p+p collisions. The error shown is due to the available statistics in the simulation.

Lambda* reconstruction efficiency times detector acceptance as a function of pT for minimum bias d+Au collisions. The error shown is due to the available statistics in the simulation.

Sigma*+- reconstruction efficiency times detector acceptance as a function of pT for minimum bias d+Au collisions. The error shown is due to the available statistics in the simulation.

Sigma*+- reconstruction efficiency times detector acceptance as a function of pT for d+Au at 40-100% centrality. The error shown is due to the available statistics in the simulation.

Sigma*+- reconstruction efficiency times detector acceptance as a function of pT for d+Au at 20-40% centrality . The error shown is due to the available statistics in the simulation.

Sigma*+- reconstruction efficiency times detector acceptance as a function of pT for d+Au at 0-20% centrality . The error shown is due to the available statistics in the simulation.

Raw pi+pi- invariant mass distribution at midrapidity after subtraction of the like-sign reference distribution for minimum bias d+Au collisions with the hadronic cocktail fit.

rho0 mass as a function of pT at |y| < 0.5 for d+A at 0-20 % centrality.

rho0 mass as a function of pT at |y| < 0.5 for d+A at 20-40 % centrality.

rho0 mass as a function of pT at |y| < 0.5 for d+A at 40-100 % centrality.

rho0 mass as a function of pT at |y| < 0.5 for minimum bias d+A collisions.

rho0 mass as a function of pT at |y| < 0.5 for p+p collisions.

rho0 mass as a function of pT at |y| < 0.5 for high multiplicity p+p collisions.

rho0 mass as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Invariant pi+pi- mass distribution after background subtraction for minimum bias p+p collisions measured by NA27.

Mixed-event background subtracted K pi raw invariant mass distribution for a particular K*0 pT bin at |y| < 0.5 for minimum bias d+Au interactions.

Mixed-event background subtracted KS*0 pi raw invariant mass distribution integrated over nvariant mass distribution for a particular K*+- pT bin at |y| < 0.5 for minimum bias d+Au interactions.

K*0 mass as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

K*+- mass as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

K*0 width as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

K*+- width as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Mixed-event background subtracted p pi raw invariant mass distribution for a particular pT bin at |y| < 0.5 for minimum bias d+Au interactions.

Mixed-event background subtracted p pi raw invariant mass distribution for a particular pT bin at |y| < 0.5 for minimum bias p+p collisions.

Delta++ mass as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Delta++ width as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Delta++ mass as a function of pT at |y| < 0.5 for minimum bias p+p collisions.

Delta++ width as a function of pT at |y| < 0.5 for minimum bias p+p collisions.

Mixed-event background subtracted Lambda pi raw invariant mass distribution integrated over the Sigma* pT at |y| < 0.5 for minimum bias d+Au interactions.

Mixed-event background subtracted p K raw invariant mass distribution integrated over the Lambda* pT at |y| < 0.5 for minimum bias d+Au interactions.

Sigma* mass as a function of pT at |y| < 0.75 for minimum bias d+Au collisions.

Lambda* mass as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Lambda* width as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

rho0 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 40-100% centrality.

rho0 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 20-40% centrality.

rho0 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 0-20% centrality.

rho0 invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

(K*0+K*0bar)/2 invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

(K*0+K*0bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 0-20% centrality.

(K*0 + K*0bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 20-40% centrality.

(K*0 + K*0bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 40-100% centrality.

(K*+ + K*-)/2 invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

(K*+ +K*-)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 0-20% centrality.

(K*+ +K*-)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 20-40% centrality.

(K*+ +K*-)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 40-100% centrality.

(Delta++ + Delta--bar)/2 invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

(Delta++ + Delta--bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 0-20% centrality.

(Delta++ + Delta--bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 20-40% centrality.

(Delta++ + Delta--bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 40-100% centrality.

(Delta++ + Delta--bar)/2 invariant yields as a function of pT at |y| < 0.5 for minimum bias p+p collisions. The errors shown also include the systematic uncertainties.

Sigma* invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Sigma* invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 0-20% centrality.

Sigma* invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 20-40% centrality.

Sigma* invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 40-100% centrality.

Lambda* invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Rescaled spectra of pion minimum bias d+Au.

Rescaled spectra of p minimum bias d+Au.

Rescaled spectra of rho0 minimum bias d+Au.

Rescaled spectra of K* minimum bias d+Au.

Rescaled spectra of Delta minimum bias d+Au.

Rescaled spectra of Sigma* minimum bias d+Au.

Rescaled spectra of Lambda* minimum bias d+Au.

pion as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

proton as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

K as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

rho0 as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

K* as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

pion as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

proton as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

K as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

Delta++ as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

Sigma* as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

Lambda* as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

rho0/pi ratios in p+p, various centralities in d+Au, and in peripheral Au+Au collisions as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

K*/k- ratios in p+p, various centralities in d+Au and Au+Au collisions as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Delta++/p ratios in p+p, various centralities in d+Au as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Sigma*/Lambda ratios in p+p, various centralities in d+Au, and in central Au+Au collisions as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Lambda*/Lambda ratios in p+p, minimum bias d+Au, and two different centralities in Au+Au collisions as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Delta++/p ratios normalized by their ratio measured in p+p collisions at the same beam energy as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

K*/K- ratios normalized by their ratio measured in p+p collisions at the same beam energy as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

rho0/pi- ratios normalized by their ratio measured in p+p collisions at the same beam energy as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Lambda*/Lambda ratios normalized by their ratio measured in p+p collisions at the same beam energy as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Sigma*/Lambda ratios normalized by their ratio measured in p+p collisions at the same beam energy as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

RdAu for pion.

RdAu for charged hadrons.

RdAu for charged rho0.

RdAu for pion.

RdAu for charged hadrons.

RdAu for charged K*0.

RdAu for charged K*+-.