![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0397v2-figure-8f2ca21f-d131-494e-8884-9397ea2e8970.jpg)

Nuclear Physics A

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0397v2-figure-89870cbd-c2b5-4043-978e-f3833626aff8.jpg)

Nuclear Physics A 00 (2014) 1-4

## Transverse energy (E T ) distributions at mid-rapidity in p + p , d + Au and Au + Au collisions at p s NN = 200 GeV and implications for particle production models

M. J. Tannenbaum (for the PHENIX Collaboration) 1

Physics Department, Brookhaven National Laboratory, Upton, NY 11973-5000, USA

## Abstract

Measurements of the mid-rapidity transverse energy distribution d E T = d GLYPH&lt;17&gt; are presented for p + p , d + Au, and Au + Au collisions at p s NN = 62.4-200 GeV. The E T distributions are compared with the number of participants, N part , the number of constituent-quark participants, Nqp , and the number of color-strings (Additive Quark Model - AQM) calculated from a Glauber model. For Au + Au, h d E T = d GLYPH&lt;17&gt; i = (0 5 : N part ) increases with N part , while h d E T = d GLYPH&lt;17&gt; i = Nqp is approximately constant vs. centrality for p s NN GLYPH&lt;21&gt; 62 4 GeV. This : indicates that the two component ansatz, d E T AA = d GLYPH&lt;17&gt; = ( d E T pp = d GLYPH&lt;17&gt; ) [(1 GLYPH&lt;0&gt; x ) N part = 2 + x N coll ], which has been used to represent E T distributions, is simply a proxy for Nqp , and that the N coll term does not represent a hard-scattering component in E T distributions. The d E T = d GLYPH&lt;17&gt; distributions of d + Au, and Au + Au are calculated from the measured p + p E T distribution using two models (AQM and Nqp ) that both reproduce the Au + Au data. For the asymmetric d + Au system, the Nqp model reproduces the data while the AQMdoes not.

Keywords: Transverse energy, Constituent-quarks, RHIC

## 1. Introduction

Recent PHENIX measurements of mid-rapidity transverse energy distributions d E T = d GLYPH&lt;17&gt; (more properly d E T = d GLYPH&lt;17&gt; j GLYPH&lt;17&gt; = 0 , measured in an electromagnetic calorimeter and corrected to total hadronic E T within a reference acceptance of GLYPH&lt;1&gt; GLYPH&lt;17&gt; = 1 0 : ; GLYPH&lt;1&gt; GLYPH&lt;30&gt; = 2 GLYPH&lt;25&gt; [1]) are presented for p + p , d + Au and Au + Au collisions at p s NN = 200 GeV and Au + Au collisions at p s NN = 62.4 and 130 GeV. The transverse energy E T is a multiparticle variable defined as the sum

$$E _ { T } = \sum _ { i } E _ { i } \, \sin \theta _ { i } \quad d E _ { T } ( \eta ) / d \eta = \sin \theta ( \eta ) \, d E ( \eta ) / d \eta,$$

where GLYPH&lt;18&gt; is the polar angle, GLYPH&lt;17&gt; = GLYPH&lt;0&gt; ln tan GLYPH&lt;18&gt;= 2 is the pseudorapidity, Ei is by convention taken as the kinetic energy for baryons, the kinetic energy + 2 mN for antibaryons, and the total energy for all other particles, and the sum is taken over all particles emitted into a fixed solid angle for each event.

The transverse energy, E T , was introduced by high energy physicists [2, 3] as an improved method to detect and study the jets from hard-scattering compared to high pT single particle spectra by which hard-scattering was discovered in p + p collisions and used as a hard-probe in Au + Au collisions at RHIC. However, it didn't work as expected: E T distributions are dominated by soft particles near h pT i [4].

1 A list of members of the PHENIX Collaboration and acknowledgements can be found at the end of this issue.

The significance of systematic measurements of mid-rapidity d E T = d GLYPH&lt;17&gt; and the closely related charged particle multiplicity distributions, dN ch = d GLYPH&lt;17&gt; , in A + B collisions is that they provide excellent characterization of the nuclear geometry (hence centrality) of the reaction on an event-by-event basis, and are also sensitive to the underlying reaction dynamics. For instance, measurements of dN ch = d GLYPH&lt;17&gt; in Au + Au collisions at the Relativistic Heavy Ion Collider (RHIC), as a function of centrality expressed as the number of participating nucleons, N part , do not depend linearly on N part but have a nonlinear increase of h dN ch = d GLYPH&lt;17&gt; i with increasing N part . The nonlinearity has been explained by a two component model [5, 6] proportional to a linear combination of N coll and N part , with the implication that the N coll term represents a contribution from hard-scattering. Alternatively, it has been proposed that dN ch = d GLYPH&lt;17&gt; is linearly proportional to the number of constituent-quark participants (NQP) model [7], without need to introduce a hard-scattering component. For symmetric systems such as Au + Au, the NQP model is identical to another model from the 1980s, the Additive Quark Model (AQM) [8], which is actually a model of particle production by color-strings in which only one colorstring can be attached to a constituent-quark participant. The models can be distinguished for asymmetric systems such as d + Au, since in the AQM the maximum number of color-strings is limited to the number of constituent-quarks in the lighter nucleus, or six for d + Au, while the NQP allows all the quark participants in both nuclei to emit particles. The two-component ansatz, the NQP model and the AQM will be tested with the present data.

## 2. Extreme independent models

The models mentioned above are examples of Extreme Independent Models in which the e GLYPH&lt;11&gt; ect of the nuclear geometry of the interaction can be calculated independently of the dynamics of particle production which can be taken directly from experimental measurements. The nuclear geometry is represented by the relative probability, wn per B + A interaction for a given number n of fundamental elements, in the present case, number of collisions ( N coll ), number of nucleon participants (wounded nucleon model - WNM [10]), number of constituent-quark participants ( Nqp ), number of color-strings (AQM). The dynamics of particle production, the N ch or E T distribution of the fundamental element, is taken from the measured p + p data in the same detector: e.g. the measured N ch distribution for a p + p collision represents: 1 collision; 2 participants (WNM); a predictable convolution of constituent-quark participants (NQP), or projectile-quark participants (AQM). Glauber calculations of the nuclear geometry ( wn ) and the p + p measurement provide a prediction for the B + A measurement in the same detector as the result of particle production by multiple independent fundamental elements. www.annualreviews.org/aronline Annual Reviews

I became acquainted with these models in my first Quark Matter talk (QM1984), in which I presented measurements of E T distributions from p + p and GLYPH&lt;11&gt; + GLYPH&lt;11&gt; interactions at p s NN = 31 GeV at the CERN-ISR (Fig. 1a [11, 12]). It was claimed at the meeting that the deviation from the WNM was due to jets [13], but in both proceedings [12, 13] it was demonstrated that 'there is no . . . sign of jets. This indicates that soft processes are still dominant, and that we are still legimately testing the WNM at these high values of E T .' [13]. As shown in Fig. 1a, the WNM did not follow the data but the AQM did [9]. Jets do appear in E T distributions as a break &lt; GLYPH&lt;24&gt; 10 GLYPH&lt;0&gt; 5 down in cross section (Figs. 1b,c). JETS IN HADRONIC COLLISIONS 111 clear peak  at  A~o12 =  180   is  observed,  which  ndicates that the two  clusters ¡ i are coplanar with the beam irection. d The  emergence  f  two-cluster structures in events with  large ~ ET  is  even o more  dramatically  illustrated by inspecting the transverse energy distribution  over  the calorimeter ells. c Figure  4 shows uch  a distribution for four s typical  events having  ~  ET  &gt; 100 GeV.  The  transverse energy  appears  to  be concentrated within two (or, more  rarely, three)  small angular regions. These  energy  clusters  appear  to be associated with  collimated multiparticle systems (jets), as  found  by reconstructing  the charged-particle tracks these  events  (see  Figure 5).  Furthermore,  longitudinal  shower  develop- ments, as  measured  n  the calorimeter,  are found  to  be inconsistent,

i

in general, with those of single  particles,

but consistent with those of jets containing  both charged  hadrons  and photons  (from  o  d ecay).

Annu. Rev. Nucl. Part. Sci. 1985.35:107-134. Downloaded from arjournals.annualreviews.org the  UA2

central calorimeter.

Figure 1. (a) E T distributions in p + p , GLYPH&lt;11&gt; + GLYPH&lt;11&gt; [11] at p s NN = 31 GeV, with AQM and WNM calculations [9]. (b),(c) E T distributions with breaks indicating jets: (b) p + p p s = 62.3 GeV [14]; (c) d GLYPH&lt;27&gt;= d E T (nb GeV) vs. E / T for ¯ p + p p s = 540 GeV [15].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0397v2-figure-4178cced-efea-4245-b96d-afb928525596.jpg)

## 3. E T and N ch distributions cut on centrality

At RHIC, following the style of the CERN SpS rather than the BNL-AGS fixed target heavy ion program, E T and N ch distributions were not generally shown. The measurements were presented cut in centrality in the form h dN ch AA = d GLYPH&lt;17&gt; i = ( h N part i = 2) vs. h N part i (Fig. 2), which would be a constant equal to h dN ch pp = d GLYPH&lt;17&gt; i if the WNM worked. The measurements clearly deviate from the WNM (Fig. 2a) [16]; so the PHENIX collaboration, inspired by the

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0397v2-figure-61fb7781-a360-4844-b87a-2125af681637.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0397v2-figure-230da0d0-1711-4163-9b9d-1b20bdaa5d2e.jpg)

c)

Figure 2. (a) PHENIX, Au + Au, p s NN = 130 GeV [16]; (b) ALICE, Pb + Pb, p s NN = 2.76 TeV [17]; (c) PHENIX preliminary Au + Au, p s NN = 7.7 GeV compared to the data at larger p s NN [18].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0397v2-figure-fb99f606-bce3-4f6c-b0eb-3255f5ba69b9.jpg)

preceding article in the journal [5], fit their data to the two-component ansatz:

$$d N _ { \text{ch} } ^ { A A } / d \eta = \left ( d N _ { \text{ch} } ^ { p p } / d \eta \right ) \left [ \left ( 1 - x \right ) \langle N _ { \text{part} } \rangle / 2 + x \langle N _ { \text{coll} } \rangle \right ]$$

where the N coll term implied a hard-scattering component for E T and N ch , known to be absent in p-p collisions 2 (recall Fig. 1). A decade later, the first measurement from Pb + Pb collisions with p s NN = 2.76 TeV at the LHC Fig. 2b [17], showed exactly the same shape vs. N part as the RHIC Au + Au data at p s NN = 200 GeV, although h N coll i is a factor of 1.6 larger and the hard-scattering cross section is more than a factor of 20 larger. This strongly argued against a hard-scattering component and for a nuclear geometrical e GLYPH&lt;11&gt; ect, which was reinforced by a PHENIX preliminary measurement in Au + Au at p s NN = 7.7 GeV (Fig. 2c) [18] which also showed the same shape for the evolution of h dN ch AA = d GLYPH&lt;17&gt; i = ( h N part i = 2) with N part as the p s NN = 200 and 2760 GeV measurements. It had previously been proposed that the number of constituent-quark participants provided the nuclear geometry that could explain the RHIC Au + Au data without the need to introduce a hard-scattering component [7]. However an asymmetric system is necessary in order to distinguish the NQP model from the AQM so the two models were applied to the RHIC d + Au data.

## 4. The number of constituent-quark participant model (NQP)

The massive constituent-quarks [19], which form mesons and nucleons (e.g. a proton = uud ), are relevant for static properties and soft physics with pT &lt; 2 GeV c. / They are complex objects or quasiparticles [21] made of the massless partons (valence quarks, gluons and sea quarks) of DIS [20] such that the valence quarks acquire masses GLYPH&lt;25&gt; 1 3 the nucleon mass with radii = GLYPH&lt;25&gt; 0 3 fm when bound in the nucleon. : With smaller resolution one can see inside the bag to resolve the massless partons which can scatter at large angles according to QCD. At RHIC, hardscattering is distinguishable from soft (exponential) particle production only for pT GLYPH&lt;21&gt; 2 GeV c at mid-rapidity, where / Q 2 = 2 p 2 T = 8 (GeV c) 2 / which corresponds to a distance scale (resolution) &lt; 0 07 fm. :

Astandard Monte Carlo Glauber calculation is used to assemble the initial positions of all the nucleons. Then three quarks are distributed around the center of each nucleon according to the proton charge distribution GLYPH&lt;26&gt; ( ~ r ) / e GLYPH&lt;0&gt; p 12 r = r rms , where r rms = 0 81 fm is the rms charge radius of the proton [22]. : The q -q inelastic scattering cross section is adjusted to 9.36 mb at p s = 200 GeV to give the correct p + p inelastic cross section (42 mb) and then used in the B A calculations. + Fig. 3a shows the deconvolution of the p + p E T distribution to the sum of 2-6 constituent-quark participants from which the E T distribution of a constituent-quark is determined and applied to d + Au (Fig. 3b) and

2 It was noted in Ref. [16] that hard-scattering was not a unique interpretation. The shape of the centrality dependences of d E T AA = d GLYPH&lt;17&gt; parameterized as N part GLYPH&lt;11&gt; were the same for Pb + Pb at p s NN = 17.6 GeV at the CERN SpS and Au + Au at p s NN = 130 GeV, with GLYPH&lt;11&gt; = 1 1 and : GLYPH&lt;11&gt; = 1 16 : GLYPH&lt;6&gt; 0 04, : respectively. The LHC data 10 years later [17] gave GLYPH&lt;11&gt; = 1 19 : GLYPH&lt;6&gt; 0 02 for Pb : + Pb at p s NN = 2760 GeV, again the same shape.

]

-1

1

[GeV

(a)

10

-1

200 GeV p+p

T

Au + Au (Fig. 3c) reactions in the same detector. The NQP calculations closely follow the measured d + Au and Au + Au E T distributions in shape and magnitude over a range of more than 1000 in cross section. A complete calculation was also done for the AQM which fails to describe the d + Au data (Fig. 3b). The conclusion is that it is the number of constituent-quark participants that determine the N ch and E T distributions and that the AQM calculation which describes the GLYPH&lt;11&gt; -GLYPH&lt;11&gt; data at p s NN = 31 GeV (Fig. 1a) is equivalent to the NQP in the symmetric system. 0 20 40 60 dY/dE -6 10 -5 10 -4 10 -3 10 -2 10 1 color-string 2 color-strings 3 color-strings

=0.538

AQM

∈

AQM fit sum,

E

[GeV]

T

Figure 3. PHENIX NQP calculations [1] for (a) p + p , (b) d + Au (also AQM), (c) Au + Au E T distributions at p s NN = 200 GeV.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0397v2-figure-c1ab204b-5570-4c9e-9766-0aae90d3abd7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0397v2-figure-aeb5c774-368a-4942-955f-0ffb63f4d07b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0397v2-figure-bd4a9d94-d70d-4070-aaa4-dbf814c2f43d.jpg)

The NQP model was also applied to the centrality cut PHENIX data with the result that d E T = d GLYPH&lt;17&gt; is strictly proportional to Nqp (Fig 4a) so that d E T = d GLYPH&lt;17&gt;= ( Nqp = 2) vs. Nqp is a constant for p s NN &gt; GLYPH&lt;24&gt; 27 GeV (Fig 4b). As a final touch, the ratio of the two-component ansatz with x = 0 08 to the : Nqp as a function of centrality was found to be constant to better than 1% at p s NN = 200 GeV (Fig 4c), which indicates that the ansatz works because the particular linear combination of N part and N coll turns out to be an empirical proxy for Nqp and not because the N coll term implies hard-scattering.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0397v2-figure-47f6ccc4-fe3a-4a75-9a55-33b34cee902c.jpg)

b)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0397v2-figure-565e6793-751d-4a54-aa44-7d1b9631883b.jpg)

c)

Figure 4. PHENIX [1]: (a) d E T = d GLYPH&lt;17&gt; vs. Nqp ; (b) d E T = d GLYPH&lt;17&gt;= ( Nqp = 2) vs. Nqp ; (c) ( Nqp / ansatz) / h Nqp = ansatz i vs. N part .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0397v2-figure-54bf8174-9910-49db-baeb-879d855a1586.jpg)

## References

- [1] PHENIX Collaboration, S.S.Adler et al., Phys. Rev. C 89 (2014) 044905 .
- [2] W. J. Willis, Report CRISP-72-15, BNL-16841, Upton, NY, 1972.
- [3] J.D.Bjorken, Phys. Rev. D 8 (1973) 4098 .
- [4] E.g. see Jan Rak and Michael J. Tannenbaum, High pT Physics in the Heavy Ion Era, Cambridge University Press, Cambridge, 2013.
- [5] X.-N. Wang and M. Gyulassy, Phys. Rev. Lett. 86 (2001) 3496 .
- [6] D. Kharzeev and M. Nardi, Phys. Lett. B 507 (2001) 121 .
- [7] S. Eremin and S. Voloshin, Phys. Rev. C 67 (2003) 064905 .
- [8] A. Bialas, W. Czyz, and L. Lesniak, Phys. Rev. D 25 (1982) 2328 .
- [9] T. Ochiai, Z. Phys. C 35 (1987) 209 .
- [10] A. Białas, M. Bleszy´ nski, and W. Czy˙ z, Nucl. Phys. B 111 (1976) 461 .
- [11] BCMOR Collaboration, A.L.S. Angelis et al., Phys. Lett. B 141 (1984) 140 .
- [12] BCMOR Collaboration, Michael J. Tannenbaum et al., Lect. Notes Phys. 221 (1985) 174 .
- [13] AFS Collaboration, Bruce Callen et al., Lect. Notes Phys. 221 (1985) 133 .
- [14] CMOR Collaboration, A.L.S. Angelis et al., Nucl. Phys. B 244 (1984) 1 .
- [15] L. DiLella, Annu. Rev. Nucl. Part. Sci. 35 (1985) 107 .
- [16] PHENIX Collaboration, K. Adcox et al., Phys. Rev. Lett. 86 (2001) 3500 .
- [17] ALICE Collaboration, K. Aamodt et al., Phys. Rev. Lett. 106 (2011) 032301 .
- [18] PHENIX Collaboration J. T. Mitchell et al., PoS CPOD2013 (2013) 003.
- [19] M. Gell-Mann, Phys. Lett. 8 (1964) 214 .
- [20] M. Breidenbach et al., Phys. Rev. Lett. 23 (1969) 935 .
- [21] E. V. Shuryak, Nucl. Phys. B 203 (1982) 116 .
- [22] R. Hofstadter, F. Bumiller, and M. R. Yerian, Rev. Mod. Phys. 30 (1958) 482 .
M. J. Tannenbaum
Transverse Energy ($E_T$) distributions at mid-rapidity in $p$$+$$p$, $d$$+$Au and Au$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV and implications for particle production models
Measurements of the midrapidity transverse energy distribution $d{\rm
E}_T/d\eta$ are presented for $p$$+$$p$, $d$$+$Au, and Au$+$Au Collisions at
$\sqrt{s_{_{NN}}}$=62.4--200 GeV. The ${\rm E}_T$ distributions are compared
with the number of participants, $N_{\rm part}$, the number of
constituent-quark participants, $N_{qp}$, and the number of color-strings
(Additive Quark Model-AQM) calculated from a Glauber model. For Au$+$Au,
$\langle d{\rm E}_T/d\eta/(0.5 N_{\rm part})\rangle$ increases with $N_{\rm
part}$, while $\langle d{\rm E}_T/d\eta/N_{qp}\rangle$ is approximately
constant vs. centrality for $\sqrt{s_{_{NN}}} \geq 62.4$ GeV. This indicates
that the two component ansatz, $d{\rm E}_T^{\rm AA}/d\eta=(d{\rm E}_T^{\rm
pp}/d\eta)\ [(1-x)\, N_{\rm part}/2 + x\, N_{\rm coll}]$, which has been used
to represent ${\rm E}_T$ distributions, is simply a proxy for $N_{qp}$, and
that the $N_{\rm coll}$ term does not represent a hard-scattering component in
${\rm E}_T$ distributions. The $d{\rm E}_T/d\eta$ distributions of $d$$+$Au,
and Au$+$Au are calculated from the measured $p$$+$$p$ ${\rm E}_T$ distribution
using two models (AQM and $N_{qp}$) that both reproduce the Au$+$Au data. For
the asymmetric $d$$+$Au system, the $N_{qp}$ model reproduces the data while
the AQM does not.