![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0466v1-figure-3dc071a9-16ac-4ab4-94dc-2b06dfec9fdd.jpg)

Nuclear Physics A

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0466v1-figure-8a3f0460-8821-4d7d-ad05-3afafe2b5940.jpg)

Nuclear Physics A 00 (2014) 1-4

## Systematic studies of the centrality dependence of soft photon production in Au + Au collision with PHENIX

Benjamin Bannier (for the PHENIX Collaboration) 1

Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794, USA

## Abstract

Since the earliest days of Heavy Ion Physics thermal soft photon radiation emitted during the reaction had been theorized as a smoking gun signal for formation of a quark-gluon plasma and as a tool to characterize its properties. In recent years the existence of excess photon radiation in heavy ion collisions over the expectation from initial hard interactions has been confirmed at both RHIC and LHC energies by PHENIX and ALICE respectively. There the radiation has been found to exhibit elliptic flow v 2 well above what can currently be reconciled with a picture of early emission from a plasma phase. During the 2007 and 2010 Au + Au runs PHENIX has measured a high purity sample of soft photons down to pT &gt; 0 4 GeV : = c using an external conversion method. We present recent systematic studies by PHENIX from that sample on the centrality dependence of the soft photon yield, and elliptic and triangular flow v 2 and v 3 in Au + Au collisions which fill in the experimental picture and enable discrimination of competing soft photon production scenarios.

Keywords:

## 1. Introduction

Collisions of heavy ions create dense and hot states of hadronic medium like the Quark Gluon Plasma (QGP), for which direct photons, i.e. photons not produced in decays of hadrons, provide an excellent probe: they can be produced during all stages of the interaction, and leave the medium virtually undisturbed due their vanishing interaction cross section with the hadronic medium. They can thus probe the full space-time evolution of the system. Here photons with transverse momenta pT = q p 2 x + p 2 y glyph[greaterorsimilar] 2 5 GeV : = c are produced predominately in early, hard interactions and called hard photons, while photons with smaller transverse momentum are called soft and thought to be produced from the medium. The photon yield gives experimental access to the rates of their di GLYPH&lt;11&gt; erent production processes, and their correlation with the event geometry, i.e. their elliptical and triangular flow coe GLYPH&lt;14&gt; cients v 2 and v 3 , are sensitive to the dynamics of the medium. While direct photon yield and flow integrate over all production channels, both quantities measured together strongly constrain possible production scenarios.

Measurements of soft photons are notoriously di GLYPH&lt;14&gt; cult in electromagnetic calorimeters due to large contamination from misidentified hadrons and a deteriorating energy resolution. Instead we here reconstruct photons from external conversions to electrons and positrons. For the 2007 and 2010 RHIC runs the HBD detector was installed in the PHENIX detector providing spatially well-defined conversion locations on a cylindrical shell R = 60 cm from the beam pipe with X X = 0 = 2 to 3% [1]. We identify conversion pairs by their characteristic apparent pair masses

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

(a) Correlation between pair mass assuming production of electron and positron at the nominal interaction point M cgl , and at the HBD shell M atm . The concentration of photons with 10 MeV = c 2 &lt; M cgl &lt; 15 MeV = c 2 and M atm &lt; 5 MeV = c 2 corresponds to conversions in the HBD shell; the opposite concentration at 10 MeV = c 2 &lt; M atm &lt; 15 MeV = c 2 and M cgl &lt; 5 MeV = c 2 is due to decay photons from GLYPH&lt;25&gt; 0 Dalitz decays, GLYPH&lt;25&gt; 0 ! GLYPH&lt;13&gt; ( ee ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0466v1-figure-57daa3d8-7fce-4c8e-940e-d515c78d1429.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0466v1-figure-e1a07800-e307-4664-a331-1ad0302acf72.jpg)

(b) R GLYPH&lt;13&gt; for di GLYPH&lt;11&gt; erent centrality classes from the 2007 and 2010 runs, compared to the virtual photon result from PHENIX [2]. Here and later statistical uncertainties are shown as bars, systematic uncertainties as boxes. Both results agree within uncertainties.

Figure 1.

(opening angles) at the nominal interaction vertex and at the HBD detector shell. Momenta of electrons and positrons can be calculated assuming they came from either the vertex or the HBD shell, and their invariant pair mass (opening angle) at the vertex and the HBD shell can be compared, see Fig. 1a. Applying a simultaneous selection on both mass variables allows a clean separation of electron-positron pairs from Dalitz decays and external photon conversions with the level of background &lt; 1% in the conversion sample while maintaining a good photon momentum resolution.

## 2. Direct photon yield

While by selecting on the two apparent pair masses we have already arrived at a very pure photon sample its relation to the actual photon production rate Y incl GLYPH&lt;13&gt; depends on detector-specific quantities, namely the photon conversion probability p conv , the geometrical acceptance of the detector for electron-positron pairs aee , and the e GLYPH&lt;14&gt; ciency of the used experimental cuts " ee , with the exact values only known inside relatively large systematic uncertainties. Instead of attempting to accurately determine these correction factors we follow a di GLYPH&lt;11&gt; erent approach: in addition to the raw yield of photons N incl GLYPH&lt;13&gt; , we measure a raw yield of photons from the decay GLYPH&lt;25&gt; 0 ! GLYPH&lt;13&gt;GLYPH&lt;13&gt; , N GLYPH&lt;25&gt; 0 GLYPH&lt;13&gt; which we extract by pairing one photon reconstructed in a conversion pair with another photon reconstructed in the PHENIX calorimeters with very loose cuts and estimation of the combinatorial background with mixed-event photon-photon pairs. The raw GLYPH&lt;25&gt; 0 -tagged yield N GLYPH&lt;25&gt; 0 GLYPH&lt;13&gt; is related to the yield Y GLYPH&lt;25&gt; 0 GLYPH&lt;13&gt; by the same detector-dependent factors as the inclusive photon yield, and additionally a conditional acceptance factor h " f i quantifying the probability to reconstruct both photons from a GLYPH&lt;25&gt; 0 decay, given that one photon was already reconstructed in a conversion pair. Since we use only very loose cuts to select calorimeter photons we can trade less dependence on systematic uncertainties for reduced statistical significance via the signal-to-background ratio in the GLYPH&lt;25&gt; 0 -tagged sample. The shared factors then drop out in the ratio of both quantities so that we can formulate a quantity R GLYPH&lt;13&gt; ,

$$R _ { \gamma } = \frac { Y _ { \gamma } ^ { \text{incl} } } { Y _ { \gamma } ^ { \text{decay} } } = \frac { Y _ { \gamma } ^ { \text{incl} } / Y _ { \gamma } ^ { \text{n} ^ { 0 } } } { Y _ { \gamma } ^ { \text{decay} } / Y _ { \gamma } ^ { \text{n} ^ { 0 } } } = \frac { \frac { N _ { \gamma } ^ { \text{incl} } / P _ { \text{com} ^ { a _ { e _ { e _ { e _ { e _ { e _ { e _ { f } } } } } } } } } { Y _ { \gamma } ^ { \text{decay} } } } = \frac { \langle \varepsilon f \rangle \frac { N _ { \gamma } ^ { \text{incl} } } { N _ { \gamma } ^ { \text{n} ^ { 0 } } } } { Y _ { \gamma } ^ { \text{decay} } }$$

(a) Direct photon pT spectra for di GLYPH&lt;11&gt; erent centrality classes. The shaded bands indicate N coll -scaled fits to PHENIX pp data. While the direct photon yield varies over two orders of magnitude between centralities the shape shows little change, also see the text.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0466v1-figure-80449eb7-d010-468c-b7ea-ec9e892815bf.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0466v1-figure-4acea833-a965-488e-aa58-43a283f7c405.jpg)

(b) The pT -integrated direct photon yield for di GLYPH&lt;11&gt; erent lower integration limits, note logarithmic axes. The dashed lines indicate independent fits to each set up measurements. The centrality-dependence of the integrated direct photon yield shows no dependence on the lower integration limit outside of uncertainties.

Figure 2. Direct photon yield

Here we have used the yield of photons from the decay of any hadron Y decay GLYPH&lt;13&gt; which can be calculated from the known hadron yields and their branching ratios to photons. The numerator of the RHS of Eq. (1) depends only on measured raw yields and the conditional acceptance h " f i which has to be determined in a Monte Carlo simulation of the detector; the denominator depends on known yields and branching ratios and can be calculate in e.g. a simple phase space simulation. With these definition any measurement R GLYPH&lt;13&gt; &gt; 1 corresponds to a direct photon signal. Our results for R GLYPH&lt;13&gt; are show in Fig. 1b. We observe a substantial direct photon signal.

From R GLYPH&lt;13&gt; we can calculate the direct photon yield shown in Fig. 2a,

$$Y _ { \gamma } ^ { \text{direct} } = ( R _ { \gamma } - 1 ) Y _ { \gamma } ^ { \text{decay} }$$

and analyze its centrality-dependence. We find that the direct photon excess over the N coll -scaled pp yield has inverse slopes roughly independent of centrality, (239 GLYPH&lt;6&gt; 25 GLYPH&lt;6&gt; 7) MeV = c (0-20%), (260 GLYPH&lt;6&gt; 33 GLYPH&lt;6&gt; 8) MeV = c (20-40%), (225 GLYPH&lt;6&gt; 28 GLYPH&lt;6&gt; 6) MeV = c (40-60%), and (238 GLYPH&lt;6&gt; 50 GLYPH&lt;6&gt; 6) MeV = c (60-92%), and the pT -integrated yield of the excess over the N coll -scaled pp yield has a power-law dependence on the number of participants N part , N GLYPH&lt;13&gt; / N GLYPH&lt;11&gt; part with a power larger than that of hadrons, GLYPH&lt;11&gt; = 1 48 : GLYPH&lt;6&gt; 0 08(stat) : GLYPH&lt;6&gt; 0 04(syst), see Fig. 2b. :

## 3. Direct photon v 2 and v 3

The elliptical and triangular flow coe GLYPH&lt;14&gt; cients v 2 and v 3 of direct photons can be calculated from the raw inclusive photon flow coe GLYPH&lt;14&gt; cients v incl n , the expectation for photons from hadron decays v decay n , and the known composition of

v

Figure 3. The direct photon v 2 ( top ) and v 3 ( bottom ) as a function of the photon pT in di GLYPH&lt;11&gt; erent centrality classes. The results from this analysis are shown with circle markers, from a preliminary calorimeter analysis as squares. Both results are consistent within systematic uncertainties.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0466v1-figure-7e7216ae-3641-4935-8141-876c32b904b4.jpg)

the inclusive photon sample quantified by R GLYPH&lt;13&gt; ,

$$v _ { n } ^ { \text{direct} } = \frac { R _ { \gamma } v _ { n } ^ { \text{incl} } - v _ { n } ^ { \text{decay} } } { R _ { \gamma } - 1 }$$

Here the vn can be calculated from the angles between the n -th order event plane measured at forward rapidities 1 0 : &lt; GLYPH&lt;17&gt; j j &lt; 2 8, : n , and the photon direction GLYPH&lt;30&gt; n with a Fourier decomposition, v 0 n = h cos 2( GLYPH&lt;30&gt; n GLYPH&lt;0&gt; n ) i . To obtain the actual vn we perform resolution corrections of the raw v 0 n with the 3-subevent method [3], taking the di GLYPH&lt;11&gt; erence to results from a 2-subevent resolution correction into account the systematic uncertainties. The v decay n can be calculated in a phase space simulation from the known vn and yields of the parent hadrons and their branching ratios to photons by measuring the decay photon GLYPH&lt;30&gt; n against the known event planes, i.e. with perfect resolution. Our result for the direct photon v 2 and v 3 are shown in Fig. 3. We observe markedly positive, non-zero coe GLYPH&lt;14&gt; cients which remain large down to low pT across all centralities.

## 4. Conclusions

We have extracted a high-purity sample of soft photons and simultaneously measured the direct photon yield and the direct photon elliptical and triangular flow coe GLYPH&lt;14&gt; cients v 2 and v 3 , and extended the measurements in the soft regime. We find a substantial direct photon signal consistent with an earlier measurement using virtual photons [2], and with pT -integrated yields growing with the number of participants N part faster than the yield of soft hadrons. The shape of the direct photon spectra shows no changes outside of uncertainties across centralities. The coe GLYPH&lt;14&gt; cients of the elliptical flow measured in the same direct photon sample show markedly positive values, consistent with results from a virtual photon analysis [4] and preliminary results from an analysis of photons measured in the PHENIX calorimeters. While as one might expected from the more eccentric collision geometry we find increasing v 2 values when going to more peripheral collision, the v 2 of direct photons appears to show less pT dependence than that of soft hadrons, even with an indication of flattening towards the smallest pT as already indicated by earlier measurements [5].

## References

- [1] W. Anderson, B. Azmoun, A. Cherlin, C. Chi, Z. Citron, et al., Nucl.Instrum.Meth. A646 , 35 (2011), 1103.4277 .
- [3] A. M. Poskanzer and S. Voloshin, Phys.Rev. C58 , 1671 (1998), nucl-ex/9805001
- [2] A. Adare et al. (PHENIX Collaboration), Phys.Rev.Lett. 104 , 132301 (2010), 0804.4168 .
- [4] A. Adare et al. (PHENIX Collaboration), Phys.Rev.Lett. 109 , 122302 (2012), 1105.4126 .

.

- [5] A. Drees (PHENIX collaboration), Nucl.Phys. A910-911 , 179 (2013).
Benjamin Bannier
Systematic studies of the centrality dependence of soft photon production in Au+Au collision with PHENIX
Since the earliest days of Heavy Ion Physics thermal soft photon radiation
emitted during the reaction had been theorized as a smoking gun signal for
formation of a quark-gluon plasma and as a tool to characterize its properties.
In recent years the existence of excess photon radiation in heavy ion
collisions over the expectation from initial hard interactions has been
confirmed at both RHIC and LHC energies by PHENIX and ALICE respectively. There
the radiation has been found to exhibit elliptic flow $v_2$ well above what can
currently be reconciled with a picture of early emission from a plasma phase.
During the 2007 and 2010 Au+Au runs PHENIX has measured a high purity sample of
soft photons down to $p_T>0.4\,\text{GeV}/c$ using an external conversion
method. We present recent systematic studies by PHENIX from that sample on the
centrality dependence of the soft photon yield, and elliptic and triangular
flow $v_2$ and $v_3$ in Au+Au collisions which fill in the experimental picture
and enable discrimination of competing soft photon production scenarios.