Enhanced production of direct photons in Au+Au collisions at =200 GeV

1. Enhanced production of direct photons in Au+Au collisions at =200 GeV Y. Akiba (RIKEN/RBRC) for PHENIX Collaboration 2008.04.25

2. Thermal Photons from the hot matter thermal: Decay photons hard: High energy density matter is formed at RHIC If the matter is thermailzed, it should emit “thermal radiation”, The temperature of the matter can directly be measured from the spectrum of thermal photon. Thermal photons can be the dominant source of direct photon for 1<pT<3 GeV/c at RHIC energies.

3. Thermal photons (theory prediciton) S.Turbide et al PRC 69 014903 • It is predicted that themal photons from QGP can be the dominant source of direct photons for 1<pT<3 GeV/c • Higher pT: pQCD photon • Lower pT: from hadronic phase • Recently, other sources, such as jet-medium interaction are discussed • Measurement is difficut since the expected signal is only 1/10 of photons from hadron decays

4. Photon measurement in PHENIX Extended in RUN5 data p+p Au+Au PHENIX measured direct photons both in p+p and Au+Au Good agreement with NLO pQCD and Ncoll scaling at high pT Measurement is limited to pT > 4-5 GeV/c

5. Alternative method --- meaure virtual photon • Source of real photon should also be able to emit virtual photon • If the Q2 (=m2) of virtual photon is sufficiently small, the source strength should be the same • The ratio of real photon and quasi-real photon can be calculated by QED  Real photon yield can be measured from virtual photon yield, which is observed as low mass e+e- pairs

6. Virtual Photon Measurement • Any source of real g can emit g* with very low mass. • Relation between the g* yield and real photon yield is known. Eq. (1) S : Process dependent factor • Case of Hadrons • Obviously S = 0 at Mee > Mhadron • Case of direct g* • If pT2>>Mee2 • Possible to separate hadron decay components from real signal in the proper mass window.

7. Not a new idea J.H.Cobb, et al, PL 78B, 519 (1978) g/p0 = 10% Dalitz g/p0 = 0.53 ±0.92% (2< pT < 3 GeV/c) The idea of measuring direct photon via low mass lepton pair is not new one. It is as old as the concept of direct photon. This method is first tried at CERN ISR in search for direct photon in p+p at 55GeV. They look for e+e- pairs for 200<m<500 MeV, and they set one of the most stringent limit on direct photon production at low pT Later, UA1 measured low mass muon pairs and deduced the direct photon cross section.

8. Measurement of low mass electron pairs arXiv: 0802.0050 arXiv: 0706.3034 p+p Au+Au • Real signal • di-electron continuum • Background sources • Combinatorial background • Material conversion pairs • Additional correlated background • Visible in p+p collisions • Cross pairs from decays with 4 electrons in the final state • Pairs in same jet or back-to-back jet

9. Enhancement of almost real photon pp Au+Au (MB) • Kinematic region of e+e- pairs • m<300 MeV and 1<pT<5 GeV/c • In this kinematic region, the S/B of the continuum is at least 10% (combinatorial BG is not a problem) • p+p • Good agreement of p+p data and hadronic decay cocktail • Small excess in p+p at large mee and high pT • Au+Au • Clear enhancement visible above for all pT 1 < pT < 2 GeV 2 < pT < 3 GeV 3 < pT < 4 GeV 4 < pT < 5 GeV

10. Possible sources of the excess • Internal conversion of direct photon A source of real photon should also produce quasi-real virtual photon • However, presence of virtual photon does not necessarily mean that it is related to real photon Example: q+q  e+e- , p +p-  e+e- BUT if they contribute for m<300MeV, the effective mass quark should be smaller than 150MeV and pion mass should be strongly modified…

11. Is excess low mass enhancement? 0 < pT < 8 GeV/c 0 < pT < 0.7 GeV/c 0.7 < pT < 1.5 GeV/c 1.5 < pT < 8 GeV/c Normalized by the yield in mee < 100MeV • Au+Au • p+p PHENIX Preliminary • The low mass enhancement decreases with higher pT • We see no significant indication that this low mass enhancement contribute to m<300 MeV and pT>1 GeV/c (see next slide) • We assume that excess is entirely due to internal conversion of direct g

12. Determination of g* fraction, r Direct g*/inclusive g* is determined by fitting the following function for each pT bin. Reminder : fdirect is given by Eq.(1) with S = 1. r : direct g*/inclusive g* • Fit in 80-300MeV gives • Assuming direct g* mass shape • c2/NDF=13.8/10 • Assuming h shape instead of direct g* shape • c2/NDF=21.1/10 •  Assumption of direct g* is favorable. • the mass spectrum follows the expectation for m>300 MeV •  No significant contribution from “low mass enhancement”

13. Fit pp 1.0<pT<1.5 GeV/c Fit range: 0-300, 30-300, 50-300, 80-300, 100-300, 120-300 There is little direct photon component in this pT bin. For the last three fit ranges, chi**2/DOF ~ 1 Variation of the fit results is included in sys. error in r

14. Fit pp 3.0-4.0 GeV/c For higher pT, small direct photon contribution is revealed

15. Fit AuAu MB. 1.0<pT<1.5 Au+Au data has much larger excess

16. Systematic uncertainties • Function fc(m) and fdir(m) are both normalized to the data for m<30 MeV so the sys. error. in the absolute normalization cancels • Sources of the systematic uncertainties are • Particle composition in the hadronic background coktail • The subtraction of background • Combinatorial background • Correlated background (Jet pairs and cross pairs) • Distortion of the mass spectrum due to dead area of the detector • Efficiency correction • We evaluate the distortion of mass spectrum due to these sources, and then evaluated their effect on the direct photon fraction r.

17. Sys. error. Cocktail hadron Ratio • mee spectra with upper/lower value of particle ratio normalized to mee<30MeV • ratio: upper/nominal, lower/nominal • in the next slide we’ll show these ratios for different pT bins • fit is repeated with the distorted cocktail fc(m). The change in the fraction r is taken as the systematic error. upper lower

18. Fraction of direct photons Fraction of direct photons Compared to direct photons from pQCD p+p Consistent with NLO pQCD favors small μ Au+Au Clear excess above pQCD Au+Au (MB) p+p μ = 0.5pT μ = 1.0pT μ = 2.0pT NLO pQCD calculation is provided by Werner Vogelsang

19. Inclusive photon • To convert the direct photon fraction r to direct photon yield, we need the invariant yield of inclusive photon. • We measure inclusive photon yield from the yield of low mass electron pairs at Dalitz peak (m<30 MeV) • We remind: • For small Mee, the process dependent factor S becomes unity. This means that electron pair yield in Dalitz peak is proportional to inclusive photon yield. • C(Mmax) is the same within a few % for any photon source for Mmax = 30 MeV

20. Inclusive photon (continued) The source of the systematic error is the error in the acceptance correction de, which is estimated to be 14% The sys. uncertainty of the cocktail does not contribute to the sys. error in the inclusive photon.

21. Direct photon spectra p+p Au+Au (MB) Direct photon yield is determined as

22. Fit to the p+p spectrum This analysis PHENIX EMCal (PRL98, 012202) Fit function: To characterize the p+p data, a modifed power-law function is fit to the spectrum The fit is repeated for the upper/lower systematic error of the spectrum (mostly common-mode) to determine the systematic uncertainty of the fit. A simple power-law gives worse c2/DOF

23. Fit to Au+Au Au+Au (MB) sys. error of pp fit p+p spectrum scaled by TAA exponential scaled pp To characterize the excess of Au+Au spectrum over the TAA scaled p+p spectrum, exponentail + scaled pp is fit to the Au+Au data

24. Direct  via  * for p+p, Au+Au • New p+p result with * method agrees with NLO pQCD predictions, and with the measurement by the calorimeter • For Au+Au there is a significant low pT excess above p+p expectations • The excess above TAA scaled p+p spectrum is characterized by the exponential fit explained in the previous slides. The inverse slope and the yield of the exponential is determined.

25. Direct  via  * for p+p, Au+Au • New p+p result with * method agrees with NLO pQCD predictions, and with the measurement by the calorimeter • For Au+Au there is a significant low pT excess above p+p expectations • The excess above TAA scaled p+p spectrum is characterized by the exponential fit explained in the previous slides. The inverse slope and the yield of the exponential is determined. exp + TAA scaled pp Fit to pp NLO pQCD (W. Vogelsang)

26. Summary of the fit Significant yield of the exponential component (excess over the scaled p+p) The inverse slope is ~220MeV. (If power-law is used for the pp component, the value of T would increase by ~24MeV)

27. Theory comparison Thery compilation by D. d’Enterria and D. Peressounko EPJC46, 451 (2006) • Hydrodynamical models are compared with the data D.d’Enterria &D.Peressounko T=590MeV, t0=0.15fm/c S. Rasanen et al. T=580MeV, t0=0.17fm/c D. K. Srivastava T=450-600MeV, t0=0.2fm/c S. Turbide et al. T=370MeV, t0=0.33fm/c J. Alam et al. T=300MeV, t0=0.5fm/c • Hydrodynamical models are in qualitative agreement with the data

28. Summary and conclusion • We have measured e+e- pairs for m<300MeV and 1<pT<5 GeV/c • Excess above hadronic background is observed • Excess is much greater in Au+Au than in p+p • Treating the excess as internal conversion of direct photons, the yield of direct photon is dedued. • Direct photon yield in pp is consistent with a NLO pQCD • Direct photon yield in Au+Au is much larger. • Spectrum shape above TAA scaled pp is exponential, with inverse slope T=221 ±23(stat)±18(sys) MeV • Hydrodynamical models with Tinit=300-600MeV at t0=0.6-0.15 fm/c are in qualitative agreement with the data.