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

Cross section D(SIG)/ET(P=4) as a function of ET(P=4) for X(C=GAMMA,OBS) <= 0.75 .

Cross section D(SIG)/(ETARAP(P=4)+ETARAP(P=5))/2 as a function of (ETARAP(P=4)+ETARAP(P=5))/2 for X(C=GAMMA,OBS) > 0.75 .

Cross section D(SIG)/(ETARAP(P=4)+ETARAP(P=5))/2 as a function of (ETARAP(P=4)+ETARAP(P=5))/2 for X(C=GAMMA,OBS) <= 0.75 .

Cross section D(SIG)/X(C=P,OBS) as a function of X(C=P,OBS) for X(C=GAMMA,OBS) > 0.75 .

Cross section D(SIG)/X(C=P,OBS) as a function of X(C=P,OBS) for X(C=GAMMA,OBS) <= 0.75 .

Cross section D(SIG)/(PHI(P=4)-PHI(P=5))/2 as a function of (PHI(P=4)-PHI(P=5))/2 for X(C=GAMMA,OBS) > 0.75 .

Cross section D(SIG)/(PHI(P=4)-PHI(P=5))/2 as a function of (PHI(P=4)-PHI(P=5))/2 for X(C=GAMMA,OBS) <= 0.75 .

Cross section D(SIG)/X(C=P,OBS) as a function of X(C=P,OBS) for X(C=GAMMA,OBS) > 0.75 for the High X(C=GAMMA,OBS) 1 optimized set.

Cross section D(SIG)/X(C=P,OBS) as a function of X(C=P,OBS) for X(C=GAMMA,OBS) > 0.75 for the High X(C=GAMMA,OBS) 2 optimized set.

Cross section D(SIG)/X(C=P,OBS) as a function of X(C=P,OBS) for X(C=GAMMA,OBS) > 0.75 for the High X(C=GAMMA,OBS) 3 optimized set.

Cross section D(SIG)/X(C=P,OBS) as a function of X(C=P,OBS) for X(C=GAMMA,OBS) > 0.75 for the High X(C=GAMMA,OBS) 4 optimized set.

Cross section D(SIG)/X(C=P,OBS) as a function of X(C=P,OBS) for X(C=GAMMA,OBS) <= 0.75 for the Low X(C=GAMMA,OBS) 1 optimized set.

Cross section D(SIG)/X(C=P,OBS) as a function of X(C=P,OBS) for X(C=GAMMA,OBS) <= 0.75 for the Low X(C=GAMMA,OBS) 2 optimized set.

Cross section D(SIG)/X(C=P,OBS) as a function of X(C=P,OBS) for X(C=GAMMA,OBS) <= 0.75 for the Low X(C=GAMMA,OBS) 3 optimized set.

Cross section D(SIG)/X(C=P,OBS) as a function of X(C=P,OBS) for X(C=GAMMA,OBS) <= 0.75 for the Low X(C=GAMMA,OBS) 4 optimized set.

Cross section D(SIG)/X(C=GAMMA,OBS) as a function of X(C=GAMMA,OBS) .

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

Differential cross section DSIG/DETARAP(P=4) in bins of ETARAP(P=4) .

Differential cross section DSIG/D(ETARAP(P=5)-ETARAP(P=6)) in bins of ETARAP(P=5)-ETARAP(P=6) .

Differential cross section DSIG/D(ETARAP(P=5)-ETARAP(P=6)) in bins of ETARAP(P=4)-ETARAP(P=5) .

Differential cross section DSIG/D(ETARAP(P=5)-ETARAP(P=6)) in bins of ETARAP(P=4)-ETARAP(P=5) .

Differential cross section DSIG/D(ETARAP(P=5)-ETARAP(P=6)) in bins of ETARAP(P=4)-ETARAP(P=5) .

Differential cross section D3SIG/DQ**2/DET(P=4)**2/DETARAP(P=4) in bins of ETARAP(P=4) for Q**2 from 20 TO 40 GeV*2 and ET(P=4)**2 range from 25 to 35 GeV**2 .

Differential cross section D3SIG/DQ**2/DET(P=4)**2/DETARAP(P=4) in bins of ETARAP(P=4) for Q**2 from 20 TO 40 GeV*2 and ET(P=4)**2 range from 36 to 100 GeV**2 .

Differential cross section D3SIG/DQ**2/DET(P=4)**2/DETARAP(P=4) in bins of ETARAP(P=4) for Q**2 from 20 TO 40 GeV*2 and ET(P=4)**2 range from 100 to 400 GeV**2 .

Differential cross section D3SIG/DQ**2/DET(P=4)**2/DETARAP(P=4) in bins of ETARAP(P=4) for Q**2 from 40 to 100 GeV*2 and ET(P=4)**2 range from 25 to 35 GeV**2 .

Differential cross section D3SIG/DQ**2/DET(P=4)**2/DETARAP(P=4) in bins of ETARAP(P=4) for Q**2 from 40 to 100 GeV*2 and ET(P=4)**2 range from 36 to 100 GeV**2 .

Differential cross section D3SIG/DQ**2/DET(P=4)**2/DETARAP(P=4) in bins of ETARAP(P=4) for Q**2 from 40 to 100 GeV*2 and ET(P=4)**2 range from 100 to 400 GeV**2 .

Cross section D(SIG)/X(C=GAMMA,OBS) as a function of X(C=GAMMA,OBS) in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/Y as a function of Y in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/Y as a function of Y in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/ET(P=4) as a function of ET(P=4) in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/ET(P=5) as a function of ET(P=5) in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/ET(P=6) as a function of ET(P=6) in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/ET(P=4) as a function of ET(P=4) in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/ET(P=5) as a function of ET(P=5) in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/ET(P=6) as a function of ET(P=6) in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/ET(P=7) as a function of ET(P=7) in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/ETARAP(P=4) as a function of ETARAP(P=4) in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/ETARAP(P=5) as a function of ETARAP(P=5) in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/ETARAP(P=6) as a function of ETARAP(P=6) in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/ETARAP(P=4) as a function of ETARAP(P=4) in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/ETARAP(P=5) as a function of ETARAP(P=5) in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/ETARAP(P=6) as a function of ETARAP(P=6) in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/ETARAP(P=7) as a function of ETARAP(P=7) in two jet invariant mass regions, 25 to 50 and > 50 GeV .

Cross section D(SIG)/COS(PSI(3)) as a function of COS(PSI(3)) in two jet invariant mass regions, 25 to 50 and > 50 GeV . PSI(3) is the angle in the 3-jet CM frame between the plane containing the highest energy jet (P=4) and the beam, and the plane containing the three jets .