The measured differences of mass differences between vector and ground states.

The measured ratios of mass differences between vector and ground states.

Normalized $\sigma_{\psi(2S)}/\sigma_{J/\psi}$ in $6.5<p_T<30.0\,GeV$ and $ 1 < y_{CM} < 1.935$ as functions of normalized $\text{N}^{{\text{corr.}}}_{\text{track}}$

Normalized $\sigma_{\psi(2S)}/\sigma_{J/\psi}$ in $6.5<p_T<30.0\,GeV$ and $ -2.865 < y_{CM} < 1.935$ as functions of normalized $\text{N}^{{\text{corr.}}}_{\text{track}}$

Normalized $\sigma_{\psi(2S)}/\sigma_{J/\psi}$ in $3.0<p_T<6.5\,GeV$ and $ -2.865 < y_{CM} < -2 $ as functions of normalized $\text{N}^{{\text{corr.}}}_{\text{track}}$

Normalized $\sigma_{\psi(2S)}/\sigma_{J/\psi}$ in $3.0<p_T<6.5\,GeV$ and $ 1 < y_{CM} < 1.935 $ as functions of normalized $\text{N}^{{\text{corr.}}}_{\text{track}}$

The model-independent 95\% \CL expected and observed upper limits set on ${\sigma(\PP\to 2\Pa+\PX)\mathcal{B}^2(\Pa\to 2\PGm)\alphaGen}$ over the range $0.21 < \MPa < 9\GeV$ for the combined 2016, 2017, and 2018 analyses. Mass ranges that overlap with \JPsi and \PgU resonances are excluded from the search

The 95\% \CL observed upper limits on the effective coupling $\CahLambda$ of the ALP to the SM Higgs assuming the branching fraction of the ALP to muons is 1 (blue) and 0.1 (orange), for both the expected (dashed) and observed (solid) limits

The 95\% \CL observed upper limits on the effective coupling $\cll/\Lambda$ of the ALP to the SM leptons. The orange shaded region represents the parameter space excluded by this search under three choices of $\CahLambda$ $1\TeV^{-2}$ (solid), $0.1\TeV^{-2}$ (dashed), and $0.01\TeV^{-2}$ (dotted)

The 95\% \CL observed upper limits on $\varepsilon^{2} \mathcal{B}(\ZDark \to \sDark \overline{\sDark}) \mathcal{B}^{2}(\sDark \rightarrow 2\PGm)$ for the vector portal model as a function of the dark scalar mass $\MsDark$ and dark vector boson mass $\MZDark$. The yellow and green bands indicate one and two standard deviation values of the limit, respectively. Because the model-independent limits are calculated only up to a dimuon mass of 60\GeV, the \sDark mass considered for these limits is below 60\GeV

The 95\% \CL observed upper limits for $\sigma(\PP\to\PhOneTwo\to 2\PaOne) \mathcal{B}^2(\PaOne\to 2\PGm)$ for the NMSSM as a function of $\MPaOne$ for $\MPhOneTwo$ = 90~GEV. The data presented here reflect the results of the 2017 and 2018 datasets combined with the previously published results of the 2016 dataset in \cite{CMS:2018jid}

The 95\% \CL observed upper limits for $\sigma(\PP\to\PhOneTwo\to 2\PaOne) \mathcal{B}^2(\PaOne\to 2\PGm)$ for the NMSSM as a function of $\MPaOne$ for $\MPhOneTwo$ = 125~GEV. The data presented here reflect the results of the 2017 and 2018 datasets combined with the previously published results of the 2016 dataset in \cite{CMS:2018jid}

The 95\% \CL observed upper limits for $\sigma(\PP\to\PhOneTwo\to 2\PaOne) \mathcal{B}^2(\PaOne\to 2\PGm)$ for the NMSSM as a function of $\MPaOne$ for $\MPhOneTwo$ = 150~GEV. The data presented here reflect the results of the 2017 and 2018 datasets combined with the previously published results of the 2016 dataset in \cite{CMS:2018jid}

The 95\% \CL observed upper limits (black solid curves) from this search as interpreted in the dark SUSY scenario for the process $\PP \to \Ph \to 2\nOne \to 2\gammaDark + 2\nDark \to 4\PGm + \PX$, with \mbox{$\MnOne = 60\GeV$} and $\MnDark = 1\GeV$. The limits are presented in the plane of $\varepsilon$ and $\MgammaDark$. The color gradient represents different branching fraction assumptions for \mbox{$\mathcal{B}(\Ph \to 2\nOne \to 2\gammaDark + 2\nDark)$}, ranging from dark orange (0.05\%) to light orange (10\%). The degradation of the limit around 1\GeV is attributed to the drop in the dimuon branching fraction $\mathcal{B}(\gammaDark \to 2\mu)$ due to the dimuon resonance of the $\phi$ meson~\cite{Batell:2009yf}

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.

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

Corrected mid-rapidity (|y| < 0.5) pT spectra for $K^{+}$, $K^{−}$, $K^{0}_{S}$, Λ, Ξ, and Ω. Λ spectra that have been corrected for feed-down are shown as open symbols in the Λ panel. The dashed lines are fits using Equation 11 except for the $\Omega+\overline{\Omega}$ where the fit uses Equation 9. The error bars displayed include systematic errors while the fits were done using statistical errors only for all species except the charged kaons.

Mid-rapidity (|y| < 0.5) ratios of $\overline{\Lambda}$ to $\Lambda$ and $\overline{\Xi}^{+}$ to $\Xi^{-}$ vs $p_{T}$. The dashed lines in 10(b) and 10(c) are the errorweighted means over the measured $p_{T}$ range, 0.882±0.017 for $\overline{\Lambda}/\Lambda$, 0.921±0.062 for $\overline{\Xi}/\Xi$. Figure 10(a) shows the $p_{T}$ averaged ratio for our measured species compared with measurements from Au+Au. The error bars are statistical only. The dashed line in 10(a) is at unity for reference.

Mid-rapidity (|y| < 0.5) ratios of $\overline{\Lambda}$ to $\Lambda$ and $\overline{\Xi}^{+}$ to $\Xi^{-}$ vs $p_{T}$. The dashed lines in 10(b) and 10(c) are the errorweighted means over the measured $p_{T}$ range, 0.882±0.017 for $\overline{\Lambda}/\Lambda$, 0.921±0.062 for $\overline{\Xi}/\Xi$. Figure 10(a) shows the $p_{T}$ averaged ratio for our measured species compared with measurements from Au+Au. The error bars are statistical only. The dashed line in 10(a) is at unity for reference.

Mid-rapidity (|y| < 0.5) ratios of $\overline{\Lambda}$ to $\Lambda$ and $\overline{\Xi}^{+}$ to $\Xi^{-}$ vs $p_{T}$. The dashed lines in 10(b) and 10(c) are the errorweighted means over the measured $p_{T}$ range, 0.882±0.017 for $\overline{\Lambda}/\Lambda$, 0.921±0.062 for $\overline{\Xi}/\Xi$. Figure 10(a) shows the $p_{T}$ averaged ratio for our measured species compared with measurements from Au+Au. The error bars are statistical only. The dashed line in 10(a) is at unity for reference.

⟨$p_{T}$⟩ vs particle mass for different particles measured by STAR. Error bars include systematic errors. The ISR parameterization is given in reference [32]

⟨$p_{T}$⟩ vs charged multiplicity for $K^{+}$, $K^{−}$, $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi+\overline{\Xi}$. The points for $\Xi+\overline{\Xi}$ have been determined using only the measured region. The error bars are statistical only. See text for more details.

Ratios of multiplicity binned spectra to minimum bias spectra ($R_{pp}$) for $K^{0}_{S}$ and Λ and the ratio of the Λ spectrum to the $K^{0}_{S}$ spectrum in each multiplicity bin. See text for further details.

Ratios of multiplicity binned spectra to minimum bias spectra ($R_{pp}$) for $K^{0}_{S}$ and Λ and the ratio of the Λ spectrum to the $K^{0}_{S}$ spectrum in each multiplicity bin. See text for further details.

Ratios of multiplicity binned spectra to minimum bias spectra ($R_{pp}$) for $K^{0}_{S}$ and Λ and the ratio of the Λ spectrum to the $K^{0}_{S}$ spectrum in each multiplicity bin. See text for further details.

Parameters of ratio data to statistical model fit using THERMUS. Filled circles are ratios from $\sqrt{s}$=200 GeV collisions in STAR. Solid lines are the results from the statis- tical model fit. All ratios to the left of the vertical line were used in the fit. The $(\Omega^{-}+\overline{\Omega^{+}})/\Xi^{-}$ ratio was then predicted from the fit results. The dashed lines in the lower panel are guides for the eye at 1 σ.

Parameters of ratio data to statistical model fit using THERMUS. Filled circles are ratios from $\sqrt{s}$=200 GeV collisions in STAR. Solid lines are the results from the statis- tical model fit. All ratios to the left of the vertical line were used in the fit. The $(\Omega^{-}+\overline{\Omega^{+}})/\Xi^{-}$ ratio was then predicted from the fit results. The dashed lines in the lower panel are guides for the eye at 1 σ.

Transverse momentum distribution for $\pi^+$ production in d+Au collisions with centrality 20-40% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^+$ production in d+Au collisions with centrality 40-100% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^-$ production in d+Au minbias events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^-$ production in p+p NSD events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^-$ production in d+Au collisions with centrality 0-20% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^-$ production in d+Au collisions with centrality 20-40% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^-$ production in d+Au collisions with centrality 40-100% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for proton production in d+Au minbias events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for proton production in p+p NSD events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for proton production in d+Au collisions with centrality 0-20% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for proton production in d+Au collisions with centrality 20-40% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for proton production in d+Au collisions with centrality 40-100% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for pbar production in d+Au minbias events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for pbar production in p+p NSD events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for pbar production in d+Au collisions with centrality 0-20% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for pbar production in d+Au collisions with centrality 20-40% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for pbar production in d+Au collisions with centrality 40-100% in the mid rapidity region, $|y|<0.5$.

Nuclear modification factors $R_{dAu}$ for $\pi^+$ + $\pi^-$ production in the mid rapidity region, $|y|<0.5$. There is an overall additional normalization uncertainty of the order of 17%. For the measurement of $R_{dAu}$ there is an additional uncertainty of 5% due to uncertainty on $N_{bin}$ for minbias collisions.

Nuclear modification factors $R_{dAu}$ for $p + \bar{p}$ production in the mid rapidity region, $|y|<0.5$. There is an overall additional normalization uncertainty of the order of 17%. For the measurement of $R_{dAu}$ there is an additional uncertainty of 5% due to uncertainty on $N_{bin}$ for minbias collisions.

Ratios of $ \pi^-/\pi^+ $ in DEUT AU and P P minbias collisions in the mid rapidity region, $|y|<0.5$. In addition there are approximately 7% correlated errors for the ratios in P P collisions and about 10% in the DEUT AU collisions.

Ratios of $ \bar{p}/p $ in DEUT AU and P P minbias collisions in the mid rapidity region, $|y|<0.5$. In addition there are approximately 7% correlated errors for the ratios in P P collisions and about 10% in the DEUT AU collisions.

Ratios of $ p/\pi^+ $ in DEUT AU and P P minbias collisions in the mid rapidity region, $|y|<0.5$. In addition there are approximately 7% correlated errors for the ratios in P P collisions and about 10% in the DEUT AU collisions.

Ratios of $ \bar{p}/\pi^- $ in DEUT AU and P P minbias collisions in the mid rapidity region, $|y|<0.5$. In addition there are approximately 7% correlated errors for the ratios in P P collisions and about 10% in the DEUT AU collisions.

Upper plot: points are measured $p/E_{tower}$ electron peak position as a function of the distance to the center of the tower. The solid line is from a calculation based on a full GEANT simulation of the detector response to electrons. Lower plot: points show measured energy deposited by electrons in the tower as a function of the momentum for distances to the center of the tower smaller than 2.0 cm. The first point is the electron equivalent energy of the minimum ionizing particles. The solid line is a second order polynomial fit of the data.

Mean values from GEANT simulations of the energy deposited in the EMC by various hadronic species as a function of momentum.

Spatial profiles of energy deposition in the EMC as a function of distance ($d$) from the hit point to the center of the tower for $\pi^{+}$ and $\pi^{-}$ from simulations and for positive and negative hadrons from data. The arrow indicates the distance corresponding to the border of a tower in 0 < $\eta$ < 0.2. An overall agreement between the shapes of the profiles is observed, with a small normalization difference.

Event-by-event ratio of the reconstructed electromagnetic energy for the most central events. The solid line is a gaussian fit.

Total Transverse Energy for 0 < $\eta$ < 1. The minimum bias distribution is presented as well as the distributions for the different centrality bins (see table III). The shaded area corresponds to the 5% most central bin. The main axis scale corresponds to the $E_{T}$ measured in the detector acceptance and the bottom axis is corrected to represent the extrapolation to full azimuthal acceptance.

⟨$dE_{T}/d\eta|_{\eta=0.5}$⟩ per $N_{part}$ pair vs. $N_{part}$. Upper panel: $N_{part}$ is obtained from Monte Carlo Glauber calculations. The lines show calculations using the HIJING model [27](solid), the EKRT saturation model [7] (dotted, Eq. 8) and the two component fit (dashed, see text). Results from WA98 [14] and PHENIX [15] are also shown. The grey bands correspond to overall systematic uncertainties, independent of $N_{part}$. Error bars are the quadrature sum of the errors on the measurements and the uncertainties on $N_{part}$ calculation. Lower panel: the same data are shown as in the upper panel but using and Optical Glauber model calculation for $N_{part}$. The line shows the same result from EKRT model calculation.

⟨$dE_{T}/d\eta|_{\eta=0.5}$⟩ per $N_{part}$ pair vs. $N_{part}$. Upper panel: $N_{part}$ is obtained from Monte Carlo Glauber calculations. The lines show calculations using the HIJING model [27](solid), the EKRT saturation model [7] (dotted, Eq. 8) and the two component fit (dashed, see text). Results from WA98 [14] and PHENIX [15] are also shown. The grey bands correspond to overall systematic uncertainties, independent of $N_{part}$. Error bars are the quadrature sum of the errors on the measurements and the uncertainties on $N_{part}$ calculation. Lower panel: the same data are shown as in the upper panel but using and Optical Glauber model calculation for $N_{part}$. The line shows the same result from EKRT model calculation.

$dE_{T}/dy$ (see text for details) per $N_{part}$ pair vs $sqrt{s_{NN}}$ for central events. In this figure $dE_{T}/dy/(0.5N_{part})$ is seen to grow logarithmicaly with $sqrt{s_{NN}}$. The error bar in the STAR point represents the total systematic uncertainty. The solid line is a EKRT model prediction [7], corrected for $d\eta/dy$, for central Au+Au collisions.

Upper panel - ⟨$dE_{T}/d\eta⟩/⟨dN_{ch}/d\eta$⟩ vs $N_{part}$. Predictions from HIJING simulations for Au+Au at 200 GeV are presented. Results from WA98 [14] and PHENIX [15] are also shown. The grey band corresponds to an overall normalization uncertainty for the STAR measurement. Bottom panel - Charged hadrons mean transverse momentum as a function of $N_{part}$ [37].

Upper panel - ⟨$dE_{T}/d\eta⟩/⟨dN_{ch}/d\eta$⟩ vs $N_{part}$. Predictions from HIJING simulations for Au+Au at 200 GeV are presented. Results from WA98 [14] and PHENIX [15] are also shown. The grey band corresponds to an overall normalization uncertainty for the STAR measurement. Bottom panel - Charged hadrons mean transverse momentum as a function of $N_{part}$ [37].

⟨$dE_{T}/d\eta⟩/⟨dN_{ch}/d\eta$⟩ vs $sqrt{s_{NN}}$ for central events. The error bar in the STAR point corresponds to the systematic uncertaintiy. A constant value of ∼ 800 MeV per charged particle, within errors, characterizes transverse energy production over this full energy range.

Energy dependence of the electromagnetic fraction of the transverse energy for a number of systems spanning SPS to RHIC energy for central events.

Participant number dependence of the electromagnetic fraction of the total transverse energy. The results are consistent with HIJING within errors over the full centrality range.