Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, where all slepton flavors are mass degenerate (co-NLSP).

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the selectron signal model. Selectron ($\tilde{e}_{L, R}$) refers to the scalar superpartners of left- and right-handed electrons.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the superpartners of the left-handed electrons, $\tilde{e}_L$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the superpartners of the right-handed electrons, $\tilde{e}_R$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the smuon signal model. Smuon ($\tilde{\mu}_{L, R}$) refers to the scalar superpartners of left- and right-handed muons.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the superpartners of the left-handed muons, $\tilde{\mu}_L$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the superpartners of the right-handed muons, $\tilde{\mu}_R$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the stau signal model. Staus, $\tilde{\tau}_{1,2}$ are the mixed states of the superpartners of the left- and right-handed $\tau$ leptons, with mixing angle $\sin\theta_{\tilde\tau}=0.95$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for $\tilde{\tau}_L$ production, where $\tilde{\tau}_L$ is pure-state superpartner of the left-handed $\tau$.

The expected and observed yields in the signal regions. Combined statistical and systematic uncertainties are presented. Estimates are truncated at 0 if the size of measured systematic uncertainties would yield a negative result.

Reconstruction efficiency as a function of $|d_{0}|$ and $p_\text{T}$ for signal electrons. Monte Carlo samples with $\tilde{e}$ or $\tilde{\mu}$ with mass 400 GeV and 1 ns lifetime were used. Efficiency is defined as the number of leptons passing all signal requirements and matched to generator-level muons divided by the number of generator level leptons with $p_\text{T} > 65$ GeV, $|d_{0}| >$ 3 mm, and |$\eta$| $<$ 2.47 for electrons. No event-level selections are made. Reconstructed leptons must pass all quality criteria, including the cosmic veto. Electron selection scale factors are included on the reconstructed leptons. The overflow is not shown in these plots.

Reconstruction efficiency as a function of $|d_{0}|$ and $p_\text{T}$ for signal muons. Monte Carlo samples with $\tilde{e}$ or $\tilde{\mu}$ with mass 400 GeV and 1 ns lifetime were used. Efficiency is defined as the number of leptons passing all signal requirements and matched to generator-level muons divided by the number of generator level leptons with $p_\text{T} > 65$ GeV, $|d_{0}| >$ 3 mm, and |$\eta$| $<$ 2.5 for muons. No event-level selections are made. Reconstructed leptons must pass all quality criteria, including the cosmic veto. Muon selection scale factors are included on the reconstructed leptons. The overflow is not shown in these plots.

Acceptance for $\tilde{e}$ of various masses and lifetimes in SR-$ee$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth electrons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Acceptance for $\tilde{\mu}$ of various masses and lifetimes in SR-$\mu\mu$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth muons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Efficiency for $\tilde{e}$ of various masses and lifetimes in SR-$ee$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth electrons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Efficiency for $\tilde{\mu}$ of various masses and lifetimes in SR-$\mu\mu$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth muons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Acceptance as a function of the generator-level $p_\text{T}$ of the leading and subleading lepton in SR-$ee$ from $\tilde{e}$ decays. The plot is made from signal Monte Carlo events with $\tilde{\ell}$ with mass of 400 GeV and lifetime of 1 ns. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47 for electrons, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. The overflow is not shown in these plots.

Acceptance as a function of the generator-level $p_\text{T}$ of the leading and subleading lepton in SR-$\mu\mu$ from $\tilde{\mu}$ decays. The plot is made from signal Monte Carlo events with $\tilde{\ell}$ with mass of 400 GeV and lifetime of 1 ns. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5 for muons, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. The overflow is not shown in these plots.

Acceptance for $\tilde{\tau}$ of various masses and lifetimes in SR-$ee$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth electrons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Acceptance for $\tilde{\tau}$ of various masses and lifetimes in SR-$e\mu$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5 (2.47) for muons (electrons), $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Acceptance for $\tilde{\tau}$ of various masses and lifetimes in SR-$\mu\mu$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth muons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Efficiency for $\tilde{\tau}$ of various masses and lifetimes in SR-$ee$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth electrons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Efficiency for $\tilde{\tau}$ of various masses and lifetimes in SR-$e\mu$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5 (2.47) for muons (electrons), $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Efficiency for $\tilde{\tau}$ of various masses and lifetimes in SR-$\mu\mu$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth muons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Acceptance as a function of the generator-level $p_\text{T}$ of the leading and subleading lepton in SR-$e\mu$ from $\tilde{\tau}$ decays. The plot is made from signal Monte Carlo events with $\tilde{\ell}$ with mass of 400 GeV and lifetime of 1 ns. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5 (2.47) for muons (electrons), $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. The overflow is not shown in these plots.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane where all slepton flavors and chiralities are mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane where all slepton flavors and chiralities are mass degenerate.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{L, R}$) refer to the scalar superpartners of left- and right-handed electrons, which are assumed to be mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{L, R}$) refer to the scalar superpartners of left- and right-handed electrons, which are assumed to be mass degenerate.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{L}$) refer to the scalar superpartners of left-handed electrons.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{L}$) refer to the scalar superpartners of left-handed electrons.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{R}$) refer to the scalar superpartners of right-handed electrons. Purple denotes the region excluded by LEP.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{R}$) refer to the scalar superpartners of right-handed electrons. Purple denotes the region excluded by LEP.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{L, R}$) refer to the scalar superpartners of left- and right-handed muons, which are assumed to be mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{L, R}$) refer to the scalar superpartners of left- and right-handed muons, which are assumed to be mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{L}$) refer to the scalar superpartners of left-handed muons.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{L}$) refer to the scalar superpartners of left-handed muons.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{R}$) refer to the scalar superpartners of right-handed muons. Purple denotes the region excluded by LEP.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{R}$) refer to the scalar superpartners of right-handed muons. Purple denotes the region excluded by LEP.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\tau})$ plane. Staus, $\tilde{\tau}_{1,2}$ are the mixed states of the superpartners of the left- and right-handed $\tau$ leptons, with mixing angle $\sin(\theta_{\tilde\tau})=0.95$, which are assumed to be mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\tau})$ plane. Staus, $\tilde{\tau}_{1,2}$ are the mixed states of the superpartners of the left- and right-handed $\tau$ leptons, with mixing angle $\sin(\theta_{\tilde\tau})=0.95$, which are assumed to be mass degenerate.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\tau}_L)$ plane, where $\tilde{\tau}_L$ is the pure-state super-partner of the left-handed $\tau$. Purple denotes the region excluded by LEP. This result does not present signficant sensitivity to the pure-state superpartner of the right-handed $\tau$.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\tau}_L)$ plane, where $\tilde{\tau}_L$ is the pure-state super-partner of the left-handed $\tau$. Purple denotes the region excluded by LEP. This result does not present signficant sensitivity to the pure-state superpartner of the right-handed $\tau$.

Statistical covariance matrix of the pion form factor $|F_\pi|^2$.

The purity and momentum-resolution corrected correlation functions $C_{true}(k^{*})$ for $p-\Lambda$, $\bar{p}-\bar{\Lambda}$ (a), $\bar{p}-\Lambda$, $p-\bar{\Lambda}$ (b). Curves correspond to fits done using the Lednicky and Lyuboshitz analytical model [12].

$(p-\Lambda) \oplus (\bar{p}-\bar{\Lambda})$ and $(\bar{p}-\Lambda) \oplus (p-\bar{\Lambda})$ combined correlation functions. Correlation functions are corrected for purity and momentum resolution. Curves correspond to fits done using the Lednicky and Lyuboshitz analytical model [12].

The combined $(\bar{p}-\Lambda) \oplus (p-\bar{\Lambda})$ spin-averaged $s$-wave scattering length compared with the previous measurements for the $p-\bar{p}$ system [24-27]. The curves show the one standard deviation contours. Note that for $(\bar{p}-\Lambda) \oplus (p-\bar{\Lambda})$ only, one should read $0.1973 \times \mbox{Im}f_{0}$ instead of $\mbox{Im}f_{0}$ on the x-axis.

The combined $(\bar{p}-\Lambda) \oplus (p-\bar{\Lambda})$ spin-averaged $s$-wave scattering length compared with the previous measurements for the $p-\bar{p}$ system [24-27]. The curves show the one standard deviation contours. Note that for $(\bar{p}-\Lambda) \oplus (p-\bar{\Lambda})$ only, one should read $0.1973 \times \mbox{Im}f_{0}$ instead of $\mbox{Im}f_{0}$ on the x-axis.

Pion source size [5] compared with proton and lambda source sizes as a function of the mean transverse mass $(\langle m_{t}\rangle)$ The curve shows the $\langle m_{t}\rangle^{-1/2}$ dependence with an arbitrary normalization.

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Distribution of $\chi^{2}$/ndf for the best jet combination for simulated 1200\GeV VLQ events (red histogram) and data (black points), for 4-jet events (left), 5-jet events (center), and 6-jet events (right). The simulated signal events and data events are normalized to the same value within the displayed $\chi^{2}$/ndf range.

Distribution of $\chi^{2}$/ndf for the best jet combination for simulated 1200\GeV VLQ events (red histogram) and data (black points), for 4-jet events (left), 5-jet events (center), and 6-jet events (right). The simulated signal events and data events are normalized to the same value within the displayed $\chi^{2}$/ndf range.

Distribution of $\chi^{2}$/ndf for the best jet combination for simulated 1200\GeV VLQ events (red histogram) and data (black points), for 4-jet events (left), 5-jet events (center), and 6-jet events (right). The simulated signal events and data events are normalized to the same value within the displayed $\chi^{2}$/ndf range.

Distributions of the average reconstructed mass of VLQ candidates for the jet combination with the least $\chi^{2}$ in 4-jet (left), 5-jet (center), and 6-jet (right) multiplicity events. The red lines show the exponential fit in the range 1000-2000 GeV. The lower panels show the fractional difference, (data-fit)/fit

Distributions of the average reconstructed mass of VLQ candidates for the jet combination with the least $\chi^{2}$ in 4-jet (left), 5-jet (center), and 6-jet (right) multiplicity events. The red lines show the exponential fit in the range 1000-2000 GeV. The lower panels show the fractional difference, (data-fit)/fit

Distributions of the average reconstructed mass of VLQ candidates for the jet combination with the least $\chi^{2}$ in 4-jet (left), 5-jet (center), and 6-jet (right) multiplicity events. The red lines show the exponential fit in the range 1000-2000 GeV. The lower panels show the fractional difference, (data-fit)/fit

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Expected limits on the VLQ mass at 95% CL as a function of the branching fractions $\mathcal{B}$( B -> bH ) and $\mathcal{B}$( B -> bZ)

Observed limits on the VLQ mass at 95% CL as a function of the branching fractions $\mathcal{B}$( B -> bH ) and $\mathcal{B}$( B -> bZ)

The 95% confidence limit on the cross section for VLQ pair production as a function of VLQ mass for three branching fraction hypotheses: B( B -> bH ) = 100% (upper left), B( B -> bZ ) = 100% (upper right), and and B( B -> bH ) = B( B -> bZ ) = 50% (lower). The solid black line indicates the observed limit and the dashed line indicates the expected limit with 1 sigma (green band) and 2 sigma (yellow band) uncertainties. The theoretical cross section and its uncertainty are shown as the red line and pale red band.

The 95% confidence limit on the cross section for VLQ pair production as a function of VLQ mass for three branching fraction hypotheses: B( B -> bH ) = 100% (upper left), B( B -> bZ ) = 100% (upper right), and and B( B -> bH ) = B( B -> bZ ) = 50% (lower). The solid black line indicates the observed limit and the dashed line indicates the expected limit with 1 sigma (green band) and 2 sigma (yellow band) uncertainties. The theoretical cross section and its uncertainty are shown as the red line and pale red band.

The 95% confidence limit on the cross section for VLQ pair production as a function of VLQ mass for three branching fraction hypotheses: B( B -> bH ) = 100% (upper left), B( B -> bZ ) = 100% (upper right), and and B( B -> bH ) = B( B -> bZ ) = 50% (lower). The solid black line indicates the observed limit and the dashed line indicates the expected limit with 1 sigma (green band) and 2 sigma (yellow band) uncertainties. The theoretical cross section and its uncertainty are shown as the red line and pale red band.

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

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.

Anti-particle to particle ratios, as a function of transverse momentum for pions (a) and protons (b). Data for the four centrality classes show little centrality dependence. Errors are statistical only.

(Color online) (a) Nuclear modification factor, RAA, for charged pions ($\pi^{+}+\pi^{−}$) in Cu+Cu (filled symbols) and Au+Au (open symbols) collisions at $\sqrt{s_{NN}}$=200 GeV. Error bands are shown for most peripheral and most central Cu+Cu data to represent evolution of the systematic uncertainties for this dataset. Error boxes at $R_{AA}$=1 represent Cu+Cu scale uncertainties due to the number of collisions and from $pp$ spectra normalization.

(Color online) (b) Integrated pion $R_{AA}$ over the range $5<p_{T}<8$ GeV/c versus $N_{part}$. The bands represent the systematic uncertainty on ratios. An additional scale error due to $pp$ normalization (∼14%) is not shown.

(Color online) (a) Nuclear modification factor, RAA, for protons and anti-protons ($p + \overline{p}$) in Cu+Cu (filled symbols) and Au+Au (open symbols) collisions at $\sqrt{s_{NN}}$=200 GeV. Error band is shown for most central Cu+Cu data to represent characteristic systematic uncertainties for Cu+Cu data. Error boxes at $R_{AA}$=1 represent Cu+Cu scale uncertainties due to the number of collisions and from pp spectra normalization.

(Color online) (a) Ratio of protons and anti-protons to charged pions, versus transverse momentum, for Cu+Cu collisions at $\sqrt{s_{NN}}$=200 GeV. For clarity, only one systematic error band is shown for Cu+Cu data (most central events), uncertainties on data in other centrality bins are similar in magnitude. Systematic errors for Au+Au data are not shown, for details see [9].

(Color online) (b) Average $(p+\overline{p})/(\pi^{+}+\pi^{-})$ ratio for $3<p_{T}<4$ GeV/c is shown as function of centrality ($N_{part}$) for Cu+Cu and Au+Au data. Error band represents the $pp$ measurement with uncertainty.

The balance function in terms of $\Delta \eta$ for all charged particles with $0.2 < p_{T} < 2.0$ GeV/$c$ from central Au+Au collisions (0-5%) for $\sqrt{s_{NN}}=27$ GeV. The data are the measured balance functions corrected by subtracting balance functions calculated using mixed events. Also shown are balance functions calculated using shuffled events.

The balance function in terms of $\Delta \eta$ for all charged particles with $0.2 < p_{T} < 2.0$ GeV/$c$ from central Au+Au collisions (0-5%) for $\sqrt{s_{NN}}=39$ GeV. The data are the measured balance functions corrected by subtracting balance functions calculated using mixed events. Also shown are balance functions calculated using shuffled events.

The balance function in terms of $\Delta \eta$ for all charged particles with $0.2 < p_{T} < 2.0$ GeV/$c$ from central Au+Au collisions (0-5%) for $\sqrt{s_{NN}}=62.4$ GeV. The data are the measured balance functions corrected by subtracting balance functions calculated using mixed events. Also shown are balance functions calculated using shuffled events.

The balance function in terms of $\Delta \eta$ for all charged particles with $0.2 < p_{T} < 2.0$ GeV/$c$ from central Au+Au collisions (0-5%) for $\sqrt{s_{NN}}=200$ GeV. The data are the measured balance functions corrected by subtracting balance functions calculated using mixed events. Also shown are balance functions calculated using shuffled events.

Energy dependence of the balance function widths compared with the widths of the balance functions calculated using shuffled events. Also shown are the balance function widths calculated using UrQMD. The dashed line represents the width of the balance function calculated using shuffled events for a constant $dN/d\eta$ distribution. Error bars represent the statistical error and the shaded bands represent the systematic error.

Energy dependence of the balance function widths compared with the widths of the balance functions calculated using shuffled events. Also shown are the balance function widths calculated using UrQMD. The dashed line represents the width of the balance function calculated using shuffled events for a constant $dN/d\eta$ distribution. Error bars represent the statistical error and the shaded bands represent the systematic error.

Energy dependence of the balance function widths compared with the widths of the balance functions calculated using shuffled events. Also shown are the balance function widths calculated using UrQMD. The dashed line represents the width of the balance function calculated using shuffled events for a constant $dN/d\eta$ distribution. Error bars represent the statistical error and the shaded bands represent the systematic error.

Energy dependence of the balance function widths compared with the widths of the balance functions calculated using shuffled events. Also shown are the balance function widths calculated using UrQMD. The dashed line represents the width of the balance function calculated using shuffled events for a constant $dN/d\eta$ distribution. Error bars represent the statistical error and the shaded bands represent the systematic error.

Energy dependence of the balance function widths compared with the widths of the balance functions calculated using shuffled events. Also shown are the balance function widths calculated using UrQMD. The dashed line represents the width of the balance function calculated using shuffled events for a constant $dN/d\eta$ distribution. Error bars represent the statistical error and the shaded bands represent the systematic error.

Energy dependence of the balance function widths compared with the widths of the balance functions calculated using shuffled events. Also shown are the balance function widths calculated using UrQMD. The dashed line represents the width of the balance function calculated using shuffled events for a constant $dN/d\eta$ distribution. Error bars represent the statistical error and the shaded bands represent the systematic error.

Energy dependence of the balance function widths compared with the widths of the balance functions calculated using shuffled events. Also shown are the balance function widths calculated using UrQMD. The dashed line represents the width of the balance function calculated using shuffled events for a constant $dN/d\eta$ distribution. Error bars represent the statistical error and the shaded bands represent the systematic error.

Balance function widths for the most central events ($0-5\%$) compared with balance function widths calculated using shuffled events. Also shown are balance function widths calculated using UrQMD and shuffled UrQMD events. The dashed line represents the width of the balance function calculated using shuffled events for a constant $dN/d\eta$ distribution.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ normalized to the most peripheral centrality bin compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ normalized to the most peripheral centrality bin compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ normalized to the most peripheral centrality bin compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ normalized to the most peripheral centrality bin compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ normalized to the most peripheral centrality bin compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ normalized to the most peripheral centrality bin compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ normalized to the most peripheral centrality bin compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ normalized to the most peripheral centrality bin compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

The $m_{W}^{T}$ distribution from data (points) and simulation (shaded histograms) in the 2j1t (left) and 3j1t (right) categories for the muon (upper) and electron (lower) channels. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 2j1t category: the $|\eta|$ of the non-b-tagged jet $|\eta_{j^{'}}|$ (left) and the invariant mass of lepton and b jet momenta system (right), shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distributionfrom the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 2j1t category: the $|\eta|$ of the non-b-tagged jet $|\eta_{j^{'}}|$ (left) and the invariant mass of lepton and b jet momenta system (right), shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distributionfrom the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 2j1t category: the $|\eta|$ of the non-b-tagged jet $|\eta_{j^{'}}|$ (left) and the invariant mass of lepton and b jet momenta system (right), shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distributionfrom the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 2j1t category: the $|\eta|$ of the non-b-tagged jet $|\eta_{j^{'}}|$ (left) and the invariant mass of lepton and b jet momenta system (right), shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distributionfrom the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 3j1t category: the $p_{T}^{miss}$ in the transverse plane (left) and the value of the MVA b tagger discriminator when applied to the extra jet (right) are shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points andthe hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 3j1t category: the $p_{T}^{miss}$ in the transverse plane (left) and the value of the MVA b tagger discriminator when applied to the extra jet (right) are shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points andthe hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 3j1t category: the $p_{T}^{miss}$ in the transverse plane (left) and the value of the MVA b tagger discriminator when applied to the extra jet (right) are shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points andthe hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 3j1t category: the $p_{T}^{miss}$ in the transverse plane (left) and the value of the MVA b tagger discriminator when applied to the extra jet (right) are shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points andthe hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 3j2t category: the $|\eta|$ of the non-b-tagged jet $|\eta_{j^{'}}|$ (left) and the invariant mass of lepton and non-b-tagged jet system (right) are shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 3j2t category: the $|\eta|$ of the non-b-tagged jet $|\eta_{j^{'}}|$ (left) and the invariant mass of lepton and non-b-tagged jet system (right) are shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 3j2t category: the $|\eta|$ of the non-b-tagged jet $|\eta_{j^{'}}|$ (left) and the invariant mass of lepton and non-b-tagged jet system (right) are shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 3j2t category: the $|\eta|$ of the non-b-tagged jet $|\eta_{j^{'}}|$ (left) and the invariant mass of lepton and non-b-tagged jet system (right) are shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distribution of the multivariate discriminators, comparing data to simulation normalised after the fit procedure, for the muon channel on the left and for the electron channel on the right, for 2j1t (upper), 3j1t (middle), and 3j2t (lower). The vertical lines on the points and the hatched bands show the experimental and fit uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the fit.

Distribution of the multivariate discriminators, comparing data to simulation normalised after the fit procedure, for the muon channel on the left and for the electron channel on the right, for 2j1t (upper), 3j1t (middle), and 3j2t (lower). The vertical lines on the points and the hatched bands show the experimental and fit uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the fit.

Distribution of the multivariate discriminators, comparing data to simulation normalised after the fit procedure, for the muon channel on the left and for the electron channel on the right, for 2j1t (upper), 3j1t (middle), and 3j2t (lower). The vertical lines on the points and the hatched bands show the experimental and fit uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the fit.

Distribution of the multivariate discriminators, comparing data to simulation normalised after the fit procedure, for the muon channel on the left and for the electron channel on the right, for 2j1t (upper), 3j1t (middle), and 3j2t (lower). The vertical lines on the points and the hatched bands show the experimental and fit uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the fit.

Distribution of the multivariate discriminators, comparing data to simulation normalised after the fit procedure, for the muon channel on the left and for the electron channel on the right, for 2j1t (upper), 3j1t (middle), and 3j2t (lower). The vertical lines on the points and the hatched bands show the experimental and fit uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the fit.

Distribution of the multivariate discriminators, comparing data to simulation normalised after the fit procedure, for the muon channel on the left and for the electron channel on the right, for 2j1t (upper), 3j1t (middle), and 3j2t (lower). The vertical lines on the points and the hatched bands show the experimental and fit uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the fit.

Values of the third-row elements of the CKM matrix inferred from low-energy measurements with the respective values of the top quark decay branching fractions. The q in bs{V_{tq}} and \mathcal{B}(t o Wq) in the first column refers to b, s, and d quarks, according to the quark label shown in the header row.

For each of the production and decay vertices, the cross section times branching fraction for the corresponding signal process from simulation. The uncertainties shown includethose from the factorisation and renormalisation scales, the PDFs, and any experimental uncer-tainties, where appropriate.

The sources and relative values in percent of the systematic uncertainty in the measurement of the STb,b cross section. The uncertainties are broken up into profiled and nonprofiled sources.

The signal strengths obtained under SM assumptions.

Results obtained under SM assumptions.

Results obtained under first BSM scenario.

Results obtained under second BSM scenario.

The correlation matrix for the njet measurement. The last bin is inclusive.

The fiducial integrated signal strength and cross section The uncertainty breakdown is given in terms of statistical (stat), experimental (exp), theoretical uncertainties on the background (bkg) and on the signal (sig), and the luminosity uncertainty (lumi). The fiducial cross section and its full uncertainty in each bin are also given.

Invariant $p_{T}$-differential spectra of $\Omega- + \bar{\Omega+}$ in p+p and p+Pb at $\sqrt{s}$=5.02 TeV in |$y_{CM}$| < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in |$y_{CM}$| < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in -1.8 < $y_{CM}$ < 0

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0 < $y_{CM}$ < 1.8

Invariant $p_{T}$-differential spectra of ${K_{0}}^{S}$ in p+Pb at $\sqrt{s}$=5.02 TeV in various $y_{CM}$ ranges

Invariant $p_{T}$-differential spectra of $\Lambda + \bar{\Lambda}$ in p+Pb at $\sqrt{s}$=5.02 TeV in various $y_{CM}$ ranges

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.3 < $y_{CM}$ < 0.8

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.8 < $y_{CM}$ < 1.3

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 1.3 < $y_{CM}$ < 1.8

$h^{−}$ pseudorapidity distribution from $5\%$ most central Au+Au collisions for $p_{T} > 100$ MeV/$c$ (filled circles) and all $p_{T}$ (open squares). The correlated systematic error on the overall normalization is estimated to be $6\%$ for $p_{T} > 0.1$ GeV/$c$ and $7\%$ for $p_{T} > 0$.

$(b\rightarrow e$ $R_{AA}) / (c\rightarrow e$ $R_{AA})$ as a function of $p_{T}$.

Charged hadron azimuthal anisotropy v2 as a function of pT in 50-60% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 40-50% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 30-40% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 20-30% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 10-20% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 0-10% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 integrated over pT and eta vs centrality in Cu+Cu collisions at 200 GeV for the various methods.

Charged hadron azimuthal anisotropy v2 integrated over pT and eta vs centrality in Cu+Cu collisions at 62.4 GeV for the various methods.

Charged hadron azimuthal anisotropy v2 as a function of pT in 50-60% Cu+Cu collisions at 62.4 GeV. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 40-50% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 30-40% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 20-30% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 10-20% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 0-10% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 0-60% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.3.

K0S azimuthal anisotropy v2 vs pT in 0-60% central Cu+Cu collisions at 200 GeV.

K0S azimuthal anisotropy v2 vs pT in 0-20% central Cu+Cu collisions at 200 GeV.

K0S azimuthal anisotropy v2 vs pT in 20-60% central Cu+Cu collisions at 200 GeV.

Lambda azimuthal anisotropy v2 vs pT in 0-60% central Cu+Cu collisions at 200 GeV.

Lambda azimuthal anisotropy v2 vs pT in 0-20% central Cu+Cu collisions at 200 GeV.

Lambda azimuthal anisotropy v2 vs pT in 20-60% central Cu+Cu collisions at 200 GeV.

Xi azimuthal anisotropy v2 vs pT in 0-60% central Cu+Cu collisions at 200 GeV.

phi azimuthal anisotropy v2 vs pT in 0-60% central Cu+Cu collisions at 200 GeV.

$\Lambda$ and $\bar\Lambda$ rapidity density as a function of negative hadron multiplicity at mid-rapidity. The open symbols in the figure are data points from an independent analysis of minimum bias data set using event-mixing techniques. Errors shown are statistical only. Systematic errors on hyperon yields $10\%$, on negative hadron yields $6\%$. $\Lambda$ and $\bar\Lambda$ measured for $|y|<0.5$, negative hadrons for $|\eta|<0.1$. Dashed lines correspond to $\Lambda = 0.054 h^{-}$, $\bar\Lambda = 0.040 h^{-}$.

Fit parameters from Boltzmann and exponential fits of the $m_{T}$ spectra for $\Lambda$ and $\bar\Lambda$ at mid-rapidity ($|y|<0.5$). Only statistical errors are presented. The systematic errors on $dN/dy$ and $T$ are approximately $10\%$. The Boltzmann form is adopted in Figure 2.

The $p_{T}$ distributions at $|y|<0.5$ for minimum bias $p$+$p$ and peripheral Au+Au collisions. See text for an explanation of the functions used to fit the data. The errors shown are statistical only and smaller than the symbols.

$\rho^{0}/\pi$ ratios as a function of c.m. system energy. The ratios are from measurements in $e^{+} e^{−}$ collisions at $10.45$ GeV, $29$ GeV and $91$ GeV c.m. system energy, $p$+$p$ at $6.8$ GeV, $19.7$ GeV, $27.5$ GeV, and $52.5$ GeV, Kp at $7.82$ GeV and $\pi^{-}$p at $19.6$ GeV. The errors on the ratios at $\sqrt{s_{NN}}= 200$ GeV are the quadratic sum of the statistical and systematic errors. The ratios at $\sqrt{s_{NN}}= 200$ GeV are offset from one another for clarity.

Nuclear modification factor $R_{dA}$ for $\pi^0$, d + Au, compared to the STAR $\pi^{+-}$ [23,24] and PHENIX $\pi^0$ measurements [27,28]. Nucleon-nucleon inelastic cross section $\sigma^{NN}$ = 42 mb and expected number of collisions $<N_{coll}>$ = 7.5 +- 0.4 from Glauber model. Shaded bands in plot are $p_T$-correlated systematic uncertainties and the error bars are statistical uncertainties. Normalization uncertainty of 14.0% indicated by a shaded band around unity.

Nuclear modification factor $R_{dA}$ for eta, d + Au, compared to PHENIX $\eta$ measurements [27,28]. Nucleon-nucleon inelastic cross section $\sigma^{NN}$ = 42 mb and expected number of collisions $<N_{coll}>$ = 7.5 +- 0.4 from Glauber model. Shaded bands in plot are $p_T$-correlated systematic uncertainties and the error bars are statistical uncertainties. Normalization uncertainty of 14.0% indicated by a shaded band around unity.

Nuclear modification factor $R_{CP}$ for $\pi^0$ measured in d + Au collisions, compared to STAR $\pi^{+-}$ measurements [24]. Peripheral d + Au collisions were used where $R_{CP}$ = ratio of particle production in 40-100% centrality (P) events and 0-20% centrality (C) events. Shaded bands in plot are $p_T$-correlated systematic uncertainties and the error bars are statistical uncertainties. Normalization uncertainty of 11.0% indicated by a shaded band around unity.

Direct photon yield in p + p collisions at $\sqrt{s_{NN}}$ = 200 GeV in terms of the double ration $R_{\gamma}$. Shaded bands in plot are $p_T$-correlated systematic uncertainties and the error bars are statistical uncertainties. The curves correspond to NLO pQCD calculations of the differential cross sections for direct photon [68] and $\pi^0$ [52] production in p + p collisions for different factorization scales $\mu$ (the upper pQCD curve corresponds to $\mu$ = $p_T$/2). The upper limit of the fractional neutral hadron contamination $C_0$ is shown as the shaded band at $R_{\gamma}$ = 0.

Direct photon yield in d + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV in terms of the double ration $R_{\gamma}$. Shaded bands in plot are $p_T$-correlated systematic uncertainties and the error bars are statistical uncertainties. The curves correspond to NLO pQCD calculations of the differential cross sections for direct photon [68] and $\pi^0$ [52] production in p + p collisions for different factorization scales $\mu$ (the upper pQCD curve corresponds to $\mu$ = $p_T$/2). The upper limit of the fractional neutral hadron contamination $C_0$ is shown as the shaded band at $R_{\gamma}$ = 0.

Cross sections for direct photon production at midrapidity in p + p and d + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, compared to the PHENIX measurement [29] and to the NLO pQCD calculation [68], which was scaled with $<T_{dA}>$ [Eq. (23)] in case of d + Au collisions. The error bars are statistical the shaded bands are $p_T$-correlated systematic uncertainties. The arrows correspond to the 95% confidence limits, as defined in the text. Normalization uncertainties of 11.7% for p+p and 5.3% for d+Au are not shown. Separate columns in the table show upper 95% confidence limits. The 95% confidence limits for the cross section are used when $R_{\gamma}$ did not correspond to a significant direct photon signal, assuming that the statistical and systematic errors both followed a Gaussian distribution, and using the fact that $R_{\gamma}\geq 1$ by definition.

Upper panels: same-event (full points) and mixed-event (solid line) $K^{+}K^{-}$ invariant mass distributions at 0.6 < $p_{T}$ < 1.4 GeV/c in p + p 200 GeV collisions (a), 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 62.4 GeV collisions (60–80%) (c), and 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 200 GeV collisions (0–10%) (e). Lower panels: the corresponding $\phi$ meson mass peaks after subtracting the background. Dashed curves show a Breit-Wigner + linear background function fit in (b), (d). In (f), both linear and quadratic backgrounds are shown as dashed and dot-dashed lines, respectively.

Upper panels: same-event (full points) and mixed-event (solid line) $K^{+}K^{-}$ invariant mass distributions at 0.6 < $p_{T}$ < 1.4 GeV/c in p + p 200 GeV collisions (a), 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 62.4 GeV collisions (60–80%) (c), and 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 200 GeV collisions (0–10%) (e). Lower panels: the corresponding $\phi$ meson mass peaks after subtracting the background. Dashed curves show a Breit-Wigner + linear background function fit in (b), (d). In (f), both linear and quadratic backgrounds are shown as dashed and dot-dashed lines, respectively.

Upper panels: same-event (full points) and mixed-event (solid line) $K^{+}K^{-}$ invariant mass distributions at 0.6 < $p_{T}$ < 1.4 GeV/c in p + p 200 GeV collisions (a), 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 62.4 GeV collisions (60–80%) (c), and 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 200 GeV collisions (0–10%) (e). Lower panels: the corresponding $\phi$ meson mass peaks after subtracting the background. Dashed curves show a Breit-Wigner + linear background function fit in (b), (d). In (f), both linear and quadratic backgrounds are shown as dashed and dot-dashed lines, respectively.

Reconstruction efficiency including acceptance of $\phi$ meson as a function of $p_{T}$ in several centrality bins of Au + Au, d + Au, and p + p 200 GeV collisions.

Reconstruction efficiency including acceptance of $\phi$ meson as a function of $p_{T}$ in several centrality bins of Au + Au, d + Au, and p + p 200 GeV collisions.

Reconstruction efficiency including acceptance of $\phi$ meson as a function of $p_{T}$ in several centrality bins of Au + Au, d + Au, and p + p 200 GeV collisions.

Event plane $\Phi_{2}$ resolution as a function of centrality in Au + Au 200 GeV collisions, where the vertical axis starts from 0.3 for clarity.

$\phi−\Phi_{2}$ distribution for $\phi$ meson at 1.5 < $p_{T}$ < 2.0 GeV/c in Au + Au collisions (0–80$\%$) at 200 GeV. The line is the fitting result. Error bars are statistical only.

<cos2($\phi_{K^{+}K^{-}}$-$\Phi_{2}$)> (full red points) and <sin2($\phi_{K^{+}K^{-}}$-$\Phi_{2}$)> (open blue points) as a function of $m_{inv}$ of $K^{+}K^{-}$ pairs at 0.5 < $p_{T}$ < 1.0 GeV/c in Au + Au 200 GeV collisions (0–80\%), where the solid curve is the result of fitting by Eq. (10). The arrow shows the position of the $\phi$ invariant mass peak. The dashed line shows the zero horizontal line.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

The $\Omega/\phi$ ratio vs $p_{T}$ for three centrality bins in $\sqrt{s_{NN}} = 200 GeV$ Au + Au collisions, where the data points for 40–60$\%$ are shifted slightly for clarity. As shown in the legend, the lines represent results from Hwa and Yang [96], Ko et al. [97], and, for SCF, Refs. [98,99].

The $\Omega/\phi$ ratio vs $p_{T}$ for three centrality bins in $\sqrt{s_{NN}} = 200 GeV$ Au + Au collisions, where the data points for 40–60$\%$ are shifted slightly for clarity. As shown in the legend, the lines represent results from Hwa and Yang [96], Ko et al. [97], and, for SCF, Refs. [98,99].

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 200 GeV collisions. The top and bottom panels present $R_{cp}$ from midperipheral and most-peripheral collisions, respectively. See legend for symbol and line designations. The rectangular bands showthe uncertainties of binary and participant scalings. Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 200 GeV collisions. The top and bottom panels present $R_{cp}$ from midperipheral and most-peripheral collisions, respectively. See legend for symbol and line designations. The rectangular bands showthe uncertainties of binary and participant scalings. Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 62.4 and 200 GeV and d + Au 200 GeV collisions, where rectangular bands represent the uncertainties of binary and participant scalings (see legend). Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 62.4 and 200 GeV and d + Au 200 GeV collisions, where rectangular bands represent the uncertainties of binary and participant scalings (see legend). Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 62.4 and 200 GeV and d + Au 200 GeV collisions, where rectangular bands represent the uncertainties of binary and participant scalings (see legend). Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{AB}$ for $\phi$ in Au + Au 200 GeV and d + Au 200 GeV collisions. For comparison, data points for $\pi^{+}+\pi^{-}$ in d + Au 200 GeV and $p+\bar{p}$ in d + Au 200 GeV are also shown (see legend). Rectangular bands show the uncertainties of binary (solid line) and participant (dot-dash line) scalings. Systematic errors are included for $\phi$, $\pi^{+}+\pi^{-}$, and $p+\bar{p}$.

$p_{T}$ dependence of the nuclear modification factor $R_{AB}$ for $\phi$ in Au + Au 200 GeV and d + Au 200 GeV collisions. For comparison, data points for $\pi^{+}+\pi^{-}$ in d + Au 200 GeV and $p+\bar{p}$ in d + Au 200 GeV are also shown (see legend). Rectangular bands show the uncertainties of binary (solid line) and participant (dot-dash line) scalings. Systematic errors are included for $\phi$, $\pi^{+}+\pi^{-}$, and $p+\bar{p}$.

$p_{T}$ dependence of the elliptic flow $v_{2}$ of $\phi$, $\Lambda$, and $K_{s}^{0}$ in Au + Au collisions (0–80$\%$) at 200 GeV. Data points for $\phi$ are from the reaction plane method (full up-triangles) and invariant mass method (full circles), where data points from the reaction plane method are shifted slightly along the $x$ axis for clarity. Vertical error bars represent statistical errors, while the square bands represent systematic uncertainties. The magenta curved band represents the $v_{2}$ of $\phi$ meson from the AMPT model with a string melting mechanism [116]. The dash and dot curves represent parametrizations inspired by number-of-quark scaling ideas from Ref. [117] for NQ = 2 and NQ = 3, respectively.

Elliptic flow $v_{2}$ as a function of $p_{T}$ , $v_{2}$($p_{T}$), for the $\phi$ meson from different centralities. The vertical error bars represent the statistical errors, while the square bands represent the systematic uncertainties. For clarity, data points of 10–40$\%$ are shifted in the $p_{T}$ direction slightly.

Elliptic flow $v_{2}$ as a function of $p_{T}$ , $v_{2}$($p_{T}$), for the $\phi$ meson from different centralities. The vertical error bars represent the statistical errors, while the square bands represent the systematic uncertainties. For clarity, data points of 10–40$\%$ are shifted in the $p_{T}$ direction slightly.

The invariant yield versus transverse momentum for |y| < 1 in 0-60% centrality in Au+Au collisions (solid circles). The results are compared to high-$p_T$ (3 < $p_T$ < 10 GeV/c) results from STAR [9] (solid squares) and PHENIX data [8] (open squares). Also shown is the yield in 0-60% centrality Cu+Cu collisions for low $p_T$ (solid triangles) and high $p_T$ [48] (open triangles). The models are described in the text [49, 50]. The $J/\psi$ cross section in p+p collisions is also shown at STAR (stars) and PHENIX [51] (diamonds), and the scale is indicated on the right axis.

The $J/\psi$ nuclear modification factor versus transverse momentum for |y| < 1 and $p_T$ < 5 GeV/c for 0-60% centrality Au+Au collisions (solid circles). The data are compared with STAR high-$p_T$ (5 < $p_T$ < 10 GeV/c) results [9] and PHENIX results [8] in |y| < 0.35 (open squares). The statistical and systematic uncertainties on the baseline $J/\psi$ cross section in p+p collisions are indicated by the hatched and solid bands, respectively. Boxes on the vertical axes represent the uncertainty on $N_{coll}$ combined with the uncertainty on the inelastic cross section in p+p collisions

The $J/\psi$ nuclear modification factor versus transverse momentum for |y| < 1 and $p_T$ < 5 GeV/c for 0-20% centrality Au+Au collisions (solid circles). The data are compared with STAR high-$p_T$ (5 < $p_T$ < 10 GeV/c) results [9] and PHENIX results [8] in |y| < 0.35 (open squares). The statistical and systematic uncertainties on the baseline $J/\psi$ cross section in p+p collisions are indicated by the hatched and solid bands, respectively. Boxes on the vertical axes represent the uncertainty on $N_{coll}$ combined with the uncertainty on the inelastic cross section in p+p collisions

The $J/\psi$ nuclear modification factor versus transverse momentum for |y| < 1 and $p_T$ < 5 GeV/c for 20-40% centrality Au+Au collisions (solid circles). The data are compared with STAR high-$p_T$ (5 < $p_T$ < 10 GeV/c) results [9] and PHENIX results [8] in |y| < 0.35 (open squares). The statistical and systematic uncertainties on the baseline $J/\psi$ cross section in p+p collisions are indicated by the hatched and solid bands, respectively. Boxes on the vertical axes represent the uncertainty on $N_{coll}$ combined with the uncertainty on the inelastic cross section in p+p collisions

The $J/\psi$ nuclear modification factor versus transverse momentum for |y| < 1 and $p_T$ < 5 GeV/c for 40-60% centrality Au+Au collisions (solid circles). The data are compared with STAR high-$p_T$ (5 < $p_T$ < 10 GeV/c) results [9] and PHENIX results [8] in |y| < 0.35 (open squares). The statistical and systematic uncertainties on the baseline $J/\psi$ cross section in p+p collisions are indicated by the hatched and solid bands, respectively. Boxes on the vertical axes represent the uncertainty on $N_{coll}$ combined with the uncertainty on the inelastic cross section in p+p collisions

The $J/\psi$ invariant yield and nuclear modification factor as a function of centrality for |y| < 1 in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $J/\psi$ invariant yield $Bd^2N/(2\pi p_Tdydp_T) [(GeV/c)^{-2}]$ and nuclear modification factor as a function of transverse momentum for |y| < 1 in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $J/\psi$ invariant yield and nuclear modification factor for |y| < 1 in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV.