Charged particle mean multiplicities measured in 4x4 kinematic bins of $Q^2$ and $y$.

Charged particle mean multiplicities measured with the additional restriction to select only particles from the current region of the Breit frame $0<\eta^{*}<4$, in 4x4 kinematic bins of $Q^2$ and $y$.

Variance of the charged particle multiplicity measured in 4x4 kinematic bins of $Q^2$ and $y$.

Variance of the charged particle multiplicity measured with the additional restriction to select only particles from the current region of the Breit frame $0<\eta^{*}<4$, in 4x4 kinematic bins of $Q^2$ and $y$.

KNO function for charged particles neasured with the additional restriction to select only particles from the current region of the Breit frame $0<\eta^{*}<4$, in 4x4 kinematic bins of $Q^2$ and $y$.

Hadron entropy for charged particles selected in a sliding pseudorapidity window of 1.4 units size with transverse momenta $P_{T lab}>150$ MeV, in the accessible kinematic bins of momentum transfer $Q^2$ and $y$. The data are presented as a function of $Q^2$ and the corresponding average $x$-Bjorken.

Hadron entropy for charged particles neasured with the additional restriction to select only particles from the current region of the Breit frame $0<\eta^{*}<4$, in 4x4 kinematic bins of $Q^2$ and $y$. The data are presented as a function of $Q^2$ and the corresponding average $x$-Bjorken.

The TPC mid-rapidity multiplicity distributions ($|\eta|$ < 0.5) for the corresponding E-FTPC selected centrality bins.

Uncorrected charged particle multiplicity distribution measured in the TPC within $|\eta|$ < 0.5 in pp collisions at 200 GeV.

The ratio of the transverse overlap area $S_{⊥}$ to $(N_{part}/2)^{2/3}$ versus $(N_{part}/2)^{2/3}$.

The ratio of the transverse overlap area $S_{⊥}$ to $(N_{part}/2)^{2/3}$ versus $(N_{part}/2)^{2/3}$.

The ratio of the transverse overlap area $S_{⊥}$ to $(N_{part}/2)^{2/3}$ versus $(N_{part}/2)^{2/3}$.

The ratio of the charged pion multiplicity to the transverse overlap area $dN_{\pi}/dy/S_{⊥}$ versus $N_{coll}/N_{part}$. Errors shown are total errors, dominated by systematic uncertainties. The systematic uncertainties are correlated between Npart, Ncoll, and S⊥, and are largely canceled in the plotted ratio quantities.

The ratio of the charged pion multiplicity to the transverse overlap area $dN_{\pi}/dy/S_{⊥}$ versus $N_{coll}/N_{part}$. Errors shown are total errors, dominated by systematic uncertainties. The systematic uncertainties are correlated between Npart, Ncoll, and S⊥, and are largely canceled in the plotted ratio quantities.

The ratio of the charged pion multiplicity to the transverse overlap area $dN_{\pi}/dy/S_{⊥}$ versus $N_{coll}/N_{part}$. Errors shown are total errors, dominated by systematic uncertainties. The systematic uncertainties are correlated between Npart, Ncoll, and S⊥, and are largely canceled in the plotted ratio quantities.

Energy loss effect for $\pi^{+-}$ (a), $K^{+-}$ (b), and p and pbar (c) at mid-rapidity (|y| < 0.1) as a function of particle momentum magnitude in 200 GeV pp and 62.4 GeV central 0-5% Au+Au collisions. Only negative particles are shown; energy loss for particles and antiparticles are the same. Errors shown are statistical only. The pion energy loss is already corrected by the track reconstruction algorithm.

Energy loss effect for $\pi^{+-}$ (a), $K^{+-}$ (b), and p and pbar (c) at mid-rapidity (|y| < 0.1) as a function of particle momentum magnitude in 200 GeV pp and 62.4 GeV central 0-5% Au+Au collisions. Only negative particles are shown; energy loss for particles and antiparticles are the same. Errors shown are statistical only. The pion energy loss is already corrected by the track reconstruction algorithm.

Energy loss effect for $\pi^{+-}$ (a), $K^{+-}$ (b), and p and pbar (c) at mid-rapidity (|y| < 0.1) as a function of particle momentum magnitude in 200 GeV pp and 62.4 GeV central 0-5% Au+Au collisions. Only negative particles are shown; energy loss for particles and antiparticles are the same. Errors shown are statistical only. The pion energy loss is already corrected by the track reconstruction algorithm.

Vertex-finding efficiency (ǫvtx) and fake vertex rate (δfake) as a function of the number of good global tracks (a) and the number of good primary tracks (b). Errors shown or smaller than the point size are statistical only.

Vertex-finding efficiency (ǫvtx) and fake vertex rate (δfake) as a function of the number of good global tracks (a) and the number of good primary tracks (b). Errors shown or smaller than the point size are statistical only.

The $p_{T}$ dependent correction to particle spectra due to fake vertex events, $\epsilon_{fake}(p_{T})$, in 200 GeV minimum bias pp and d+Au collisions. Errors shown are statistical only.

The $p_{T}$ dependent correction to particle spectra due to fake vertex events, $\epsilon_{fake}(p_{T})$, in 200 GeV minimum bias pp and d+Au collisions. Errors shown are statistical only.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in pp (a) and d+Au collisions (b) at 200 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in pp (a) and d+Au collisions (b) at 200 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in 70-80% peripheral Au+Au (a) and 0-5% central Au+Au collisions (b) at 62.4 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in 70-80% peripheral Au+Au (a) and 0-5% central Au+Au collisions (b) at 62.4 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in 70-80% peripheral Au+Au (a) and 0-5% central Au+Au collisions (b) at 62.4 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in 70-80% peripheral Au+Au (a) and 0-5% central Au+Au collisions (b) at 62.4 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

The dca distributions of protons and antiprotons for 0.40 < p⊥ < 0.45 GeV/c in 200 GeV minimum bias d+Au. Errors shown are statistical only.

The dca distributions of protons and antiprotons for 0.70 < p⊥ < 0.75 GeV/c in 200 GeV minimum bias d+Au. Errors shown are statistical only.

The dca distributions of protons and antiprotons for 0.40 < p⊥ < 0.45 GeV/c in 62.4 GeV 0-5% central Au+Au collisions. Errors shown are statistical only.

The dca distributions of protons and antiprotons for 0.70 < p⊥ < 0.75 GeV/c in 62.4 GeV 0-5% central Au+Au collisions. Errors shown are statistical only.

Pion background fraction from weak decays ($\Lambda, K_S^0$) and $\mu^{+-}$ contamination as a function of p⊥ in minimum bias d+Au collisions at 200 GeV. Errors shown are statistical only.

Pion background fraction from weak decays ($\Lambda, K_S^0$) and $\mu^{+-}$ contamination as a function of p⊥ in minimum bias d+Au collisions at 200 GeV. Errors shown are statistical only.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity (|y| < 0.1) identified particle spectra in d+Au collisions at 200 GeV. The p and pbar spectra are inclusive, including weak decay products. Spectra are plotted for three centrality bins and for minimum bias events. Spectra from top to bottom are for 0-20% scaled by 4, 20-40% scaled by 2, minimum bias not scaled, and 40-100% scaled by 1/2. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the minimum bias data; the normalizations of the curves are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in d+Au collisions at 200 GeV. The p and pbar spectra are inclusive, including weak decay products. Spectra are plotted for three centrality bins and for minimum bias events. Spectra from top to bottom are for 0-20% scaled by 4, 20-40% scaled by 2, minimum bias not scaled, and 40-100% scaled by 1/2. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the minimum bias data; the normalizations of the curves are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in d+Au collisions at 200 GeV. The p and pbar spectra are inclusive, including weak decay products. Spectra are plotted for three centrality bins and for minimum bias events. Spectra from top to bottom are for 0-20% scaled by 4, 20-40% scaled by 2, minimum bias not scaled, and 40-100% scaled by 1/2. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the minimum bias data; the normalizations of the curves are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in Au+Au collisions at 62.4 GeV. The p and p spectra are inclusive, including weak decay products. Spectra are plotted for nine centrality bins, from top to bottom, 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70%, and 70-80%. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the spectra; the normalizations of the curves in (a,b) are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in Au+Au collisions at 62.4 GeV. The p and p spectra are inclusive, including weak decay products. Spectra are plotted for nine centrality bins, from top to bottom, 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70%, and 70-80%. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the spectra; the normalizations of the curves in (a,b) are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in Au+Au collisions at 62.4 GeV. The p and p spectra are inclusive, including weak decay products. Spectra are plotted for nine centrality bins, from top to bottom, 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70%, and 70-80%. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the spectra; the normalizations of the curves in (a,b) are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in Au+Au collisions at 62.4 GeV. The p and p spectra are inclusive, including weak decay products. Spectra are plotted for nine centrality bins, from top to bottom, 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70%, and 70-80%. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the spectra; the normalizations of the curves in (a,b) are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified pion spectra in Au+Au collisions at 130 GeV. Spectra are plotted for eight centrality bins, from top to bottom, 0-6%, 6-11%, 11-18%, 18- 26%, 26-34%, 34-45%, 45-58%, and 58-85%. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but they are smaller than the data point size. The curves are the Bose-Einstein fits to the spectra; the normalizations of the curves are fixed by the corresponding negative particle spectra.

Estimate of the product of the Bjorken energy density and the formation time ($\epsilon_{B_j}\tau$) as a function of centrality Npart. Errors shown are the quadratic sum of statistical and systematic uncertainties.

Estimate of the product of the Bjorken energy density and the formation time ($\epsilon_{B_j}\tau$) as a function of centrality Npart. Errors shown are the quadratic sum of statistical and systematic uncertainties.

Estimate of the product of the Bjorken energy density and the formation time ($\epsilon_{B_j}\tau$) as a function of centrality Npart. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K−/π− ratio as a function of the number of participants Npart in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K−/π− ratio as a function of the number of participants Npart in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K−/π− ratio as a function of the number of participants Npart in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K−/π− ratio as a function of dNπ/dy/S⊥ in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K−/π− ratio as a function of dNπ/dy/S⊥ in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K−/π− ratio as a function of dNπ/dy/S⊥ in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

Differential cross-sections obtained from the optical and MC Glauber calculations for Au+Au collisions at 200 GeV. Statistical errors are smaller than the point size.

Differential cross-sections obtained from the optical and MC Glauber calculations for Au+Au collisions at 200 GeV. Statistical errors are smaller than the point size.

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 4. (Color online.) Top panel: The number of photons divided by $\left\langle N_{\text{part}}\right\rangle / 2$ as a function of average number of participating nucleons for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s} _{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$ for $-3.7<\eta<-2.3 .$ Errors shown are systematic only and include those for $\left\langle N_{\text {part }}\right\rangle .$ Results from HIJING are shown as lines (solid for $\mathrm{Au}+\mathrm{Au}$ and dashed for $\mathrm{Cu}+\mathrm{Cu}$ ). Bottom panel: Same as above, for both photons and charged particles from PHOBOS [10]. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 4. (Color online.) Top panel: The number of photons divided by $\left\langle N_{\text{part}}\right\rangle / 2$ as a function of average number of participating nucleons for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s} _{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$ for $-3.7<\eta<-2.3 .$ Errors shown are systematic only and include those for $\left\langle N_{\text {part }}\right\rangle .$ Results from HIJING are shown as lines (solid for $\mathrm{Au}+\mathrm{Au}$ and dashed for $\mathrm{Cu}+\mathrm{Cu}$ ). Bottom panel: Same as above, for both photons and charged particles from PHOBOS [10]. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 5. (Color online.) Photon pseudorapidity distributions normalized by the average number of participating nucleon pairs for different collision centralities are plotted as a function of pseudorapidity shifted by the beam rapidity (-5.36 for $200 \mathrm{GeV}$ and -4.19 for $62.4 \mathrm{GeV}$ ) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ collisions at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$. Errors are systematic only, statistical errors are negligible in comparison. For clarity of presentation, results for only four centralities are shown. The $\mathrm{Cu}+$ Cu data are shifted by 0.1 unit in $\eta-y_{\text {beam }} .$ The solid line is a second order polynomial fit to the data (see text for details).

Fig. 5. (Color online.) Photon pseudorapidity distributions normalized by the average number of participating nucleon pairs for different collision centralities are plotted as a function of pseudorapidity shifted by the beam rapidity (-5.36 for $200 \mathrm{GeV}$ and -4.19 for $62.4 \mathrm{GeV}$ ) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ collisions at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$. Errors are systematic only, statistical errors are negligible in comparison. For clarity of presentation, results for only four centralities are shown. The $\mathrm{Cu}+$ Cu data are shifted by 0.1 unit in $\eta-y_{\text {beam }} .$ The solid line is a second order polynomial fit to the data (see text for details).

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for NSD collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for NSD collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -2.0 to 2.0 for NSD collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -2.4 to 2.4 for NSD collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.0 to 3.0 for NSD collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for NSD collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for NSD collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -2.0 to 2.0 for NSD collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -2.4 to 2.4 for NSD collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.0 to 3.0 for NSD collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for NSD collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for NSD collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -2.0 to 2.0 for INEL collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -2.4 to 2.4 for INEL collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -3.0 to 3.0 for INEL collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for INEL collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for INEL collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -2.0 to 2.0 for INEL collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -2.4 to 2.4 for INEL collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.0 to 3.0 for INEL collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for INEL collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for INEL collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -2.0 to 2.0 for INEL collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -2.4 to 2.4 for INEL collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.0 to 3.0 for INEL collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for INEL collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for INEL collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -2.0 to 2.0 for INEL>0 collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -2.4 to 2.4 for INEL>0 collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -3.0 to 3.0 for INEL>0 collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for INEL>0 collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for INEL>0 collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -2.0 to 2.0 for INEL>0 collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -2.4 to 2.4 for INEL>0 collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.0 to 3.0 for INEL>0 collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for INEL>0 collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for INEL>0 collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -2.0 to 2.0 for INEL>0 collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -2.4 to 2.4 for INEL>0 collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.0 to 3.0 for INEL>0 collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for INEL>0 collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for INEL>0 collisions at a centre-of-mass energy of 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in NSD collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in NSD collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in NSD collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in NSD collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in NSD collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in NSD collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in NSD collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in NSD collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in NSD collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in NSD collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in NSD collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in NSD collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in NSD collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in NSD collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in NSD collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in INEL collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in INEL collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in INEL collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in INEL collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in INEL collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in INEL collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in INEL collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in INEL collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in INEL collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in INEL collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in INEL collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in INEL collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in INEL collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in INEL collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in INEL collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in INEL>0 collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in INEL>0 collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in INEL>0 collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in INEL>0 collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in INEL>0 collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in INEL>0 collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in INEL>0 collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in INEL>0 collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in INEL>0 collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in INEL>0 collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in INEL>0 collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in INEL>0 collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in INEL>0 collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in INEL>0 collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in INEL>0 collisions at center-of-mass energy 8000 GeV.

The mean, $\langle n_{\rm gluon}^{\rm ch.} \rangle$, and first two non-trivial normalized factorial moments, $F_{\rm 2,\,gluon}$ and $F_{\rm 3,\,gluon}$, of the charged particle multiplicity distribution of gluon jets. The data have been corrected for detector acceptance and resolution, for event selection, and for gluon jet impurity.

The charged particle fragmentation function of gluon jets, ${1/N \,({\rm d}n_{\rm gluon}^{\rm ch.}/{\rm d}}x_E^*)$, for $E_{\rm g}^*$$\,=\,$14.24 and 17.72 GeV. The data have been corrected for detector acceptance and resolution, for event selection, and for gluon jet impurity.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV from Transport Model Calculations.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV from Transport Model Calculations.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations for net-baryon.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations for net-proton (No Decay).

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

$\langle p_\perp \rangle$ vs. multiplicity at $\sqrt{s} = 1800~\text{GeV}$, $|\eta| < 1$, $p_T > 0.4~\text{GeV}$.

Et EMC distributions for sqrt(sNN) = 62.4 GeV Au+Au collisions shown in 5% wide centrality bins.

Et EMC distributions for sqrt(sNN) = 62.4 GeV Au+Au collisions shown in 5% wide centrality bins.

Et EMC distributions for sqrt(sNN) = 62.4 GeV Au+Au collisions shown in 5% wide centrality bins.

Et EMC distributions for sqrt(sNN) = 62.4 GeV Au+Au collisions shown in 5% wide centrality bins.

Et EMC distributions for sqrt(sNN) = 62.4 GeV Au+Au collisions shown in 5% wide centrality bins.

Et EMC distributions for sqrt(sNN) = 62.4 GeV Au+Au collisions shown in 5% wide centrality bins.

Et EMC distributions for sqrt(sNN) = 62.4 GeV Au+Au collisions shown in 5% wide centrality bins.

Et EMC distributions for sqrt(sNN) = 62.4 GeV Au+Au collisions shown in 5% wide centrality bins.

Corrected dE_T/deta distributions at sqrt(sNN) = 200 GeV for five sectors of PbSc in Au+Au collisions.

Corrected dE_T/deta distributions at sqrt(sNN) = 200 GeV for five sectors of PbSc in p+p collisions.

Corrected dE_T/deta distributions at sqrt(sNN) = 200 GeV for five sectors of PbSc in d+Au collisions.

The number of quark participants as a function of the number of nucleon participants.

The number of quark participants as a function of the number of nucleon participants.

The number of quark participants as a function of the number of nucleon participants.

The ratio of the number of quark participants to the number of nucleon participants as a function of the number of nucleon participants.

The ratio of the number of quark participants to the number of nucleon participants as a function of the number of nucleon participants.

The ratio of the number of quark participants to the number of nucleon participants as a function of the number of nucleon participants.

dE_T/deta normalized by the number of participant pairs as a function of the number of participants.

dE_T/deta normalized by the number of participant pairs as a function of the number of participants.

dE_T/deta normalized by the number of participant pairs as a function of the number of participants.

dE_T/deta normalized by the number of participant quark pairs as a function of the number of nucleon participants.

dE_T/deta normalized by the number of participant quark pairs as a function of the number of nucleon participants.

dE_T/deta normalized by the number of participant quark pairs as a function of the number of nucleon participants.

dE_T/deta as a function of the number of quark participants.

dE_T/deta as a function of the number of quark participants.

dE_T/deta as a function of the number of quark participants.

Corrected dE_T/deta distributions at sqrt(sNN) = 200 GeV for five sectors of PbSc in p+p collisions.

Distribution of the number of quark participants in 200 GeV Au+Au collsions.

Corrected dE_T/deta distributions at sqrt(sNN) = 200 GeV for five sectors of PbSc in p+p collisions.

Corrected dE_T/deta distributions at sqrt(sNN) = 200 GeV for five sectors of PbSc in Au+Au collisions.

Corrected dE_T/deta distributions at sqrt(sNN) = 200 GeV for five sectors of PbSc in d+Au collisions.

Corrected dE_T/deta distributions at sqrt(sNN) = 200 GeV for five sectors of PbSc in d+Au collisions.

Corrected dE_T/deta distributions at sqrt(sNN) = 200 GeV for five sectors of PbSc in Au+Au collisions.

Corrected dE_T/deta distributions at sqrt(sNN) = 200 GeV for five sectors of PbSc in d+Au collisions.

Corrected dE_T/deta distributions at sqrt(sNN) = 200 GeV for five sectors of PbSc in Au+Au collisions.

Corrected dE_T/deta distributions at sqrt(sNN) = 200 GeV for five sectors of PbSc in Au+Au collisions.

Corrected dE_T/deta distributions at sqrt(sNN) = 200 GeV for five sectors of PbSc in d+Au collisions.

Corrected dE_T/deta distributions at sqrt(sNN) = 200 GeV for five sectors of PbSc in Au+Au collisions.

Corrected dE_T/deta distributions at sqrt(sNN) = 200 GeV for five sectors of PbSc in p+p collisions.

Charged particle multiplicity distribution.

Average charged particle transverse momentum as a function of the number of charged particles in the event.