# How robust is a thermal photon interpretation of the ALICE low- data?

###### Abstract

We present a rigorous theoretical analysis of the ALICE measurement of low- direct-photon production in central lead-lead collisions at the LHC with a centre-of-mass energy of TeV. Using NLO QCD, we compute the relative contributions to prompt-photon production from different initial and final states and the theoretical uncertainties coming from independent variations of the renormalisation and factorisation scales, the nuclear parton densities and the fragmentation functions. Based on different fits to the unsubtracted and prompt-photon subtracted ALICE data, we consistently find 58 MeV and MeV for the effective temperature of the quark-gluon plasma (or hot medium) at GeV and GeV as well as a power-law () behavior for GeV as predicted by QCD hard scattering.

MS-TP-13-18

## 1 Introduction

One of the goals of the physics program at the CERN Large Hadron Collider (LHC) is the study of deconfined, strongly interacting matter, which existed in the early universe and which can today be re-created and investigated in heavy-ion collisions. An important probe for this so-called quark-gluon plasma (QGP) are photons emitted from the deconfined partons before thermalisation, in the thermal bath, during expansion and cooling of the QGP, and finally from the thermal hadron gas. The transverse momentum distribution of these photons can be used to estimate the temperature of the system, although the exact interpretation is complicated by these different phases, the radial expansion, more generally the temporal evolution of phase space, and co-existing states of matter. Experimentally, an effective temperature is usually extracted from an exponential fit to the excess of direct-photon production at low transverse momentum () above the expectation from vacuum production.

In Ref. Wilde:2012wc (), the first observation of a low- direct-photon
signal at the LHC has been reported by the ALICE collaboration.^{1}^{1}1In nuclear
collision experiments, photons originating from meson decays are usually
distinguished from direct photons. The latter are divided into thermal and
non-thermal prompt (plus medium-induced)
photons, and these again into photons produced directly
in the hard collision and those coming from quark and gluon fragmentation.
The double use of the word direct can sometimes lead to confusion.
There, an inverse slope parameter of MeV has been
extracted from an exponential
fit to the photon spectrum in central (0-40%) lead-lead
collisions at TeV and low transverse momenta of
GeV. The inverse slope parameter of this
measurement is significantly higher than the one obtained previously
by the PHENIX collaboration in 0-20% central gold-gold collisions with
GeV at RHIC, which resulted in MeV
Adare:2008ab (). The latter was higher than
the transition temperature to the QGP of about MeV, but
1.5 to 3 times smaller
than the initial temperature of the dense matter due to the
space-time evolution following its initial formation;
in hydrodynamical models, which describe the data at =
200 GeV, ranges from 600 to 300 MeV depending on the
formation time, assumed to lie between and
0.6 fm d'Enterria:2005vz ().

The two general-purpose LHC experiments ATLAS and CMS have also recently measured prompt-photon production in lead-lead collisions at TeV with and 20 GeV, respectively ATLAS:2012zla (); Chatrchyan:2012vq (). In this transverse momentum range, photons (like electroweak and bosons Steinberg:2013laa (); GranierdeCassagnac:2013maa ()) are expected to be dominantly produced in hard partonic collisions and to be unaffected by the strongly interacting medium, so that they can be described with next-to-leading order (NLO) QCD programs such as JETPHOX Aurenche:2006vj ().

However, photons from electromagnetic decays of neutral mesons (mostly ) must first be removed using meson data Wilde:2012wc (), Monte Carlo simulations Adare:2008ab () and/or by applying an isolation criterion, which e.g. limits the energy of all other particles in a cone of size or 0.4 around the photon to or 5 GeV ATLAS:2012zla (); Chatrchyan:2012vq (). The remaining “prompt” photon sample then still contains both direct photons, produced at tree-level through the processes and , and photons emitted in the collinear fragmentation of quarks and gluons produced in partonic scattering processes. Isolation criteria also reduce the contribution of fragmentation photons, since these are generally produced inside a high- jet of hadrons. The separation of these direct and fragmentation photons becomes ambiguous, i.e. scheme- and scale-dependent, at NLO of perturbative QCD, but their interplay and the factorisation of the associated infrared singularities have been studied intensively theoretically Klasen:2002xb (). As a result, NLO calculations describe very well prompt-photon production in collisions at the LHC with TeV down to GeV Aad:2011tw (); ATLAS:2013ema (); Khachatryan:2010fm (); Chatrchyan:2011ue ().

At smaller transverse momenta, a transition from hard to thermal photon production is expected to occur in central heavy-ion collisions, in contrast to peripheral heavy-ion and collisions. For this reason, the transverse momentum distribution of inclusive photons has been measured by the ALICE collaboration in lead-lead collisions at TeV from 16 GeV down to 0.4 GeV. In addition, the same has been done in collisions at TeV down to 0.3 GeV. In both cases, contributions from decay photons were subtracted using neutral pion data and transverse mass () scaling Khandai:2011cf (). The data for the 0-40% most central heavy-ion collisions showed a clear deviation from the perturbative QCD prediction below -values of 4 GeV which was not observed for and peripheral (40-80%) heavy-ion collisions Wilde:2012wc (). Similar observations had been made in a preceding analysis of the PHENIX collaboration in minimum bias, 0-20% and 20-40% central gold-gold collisions at RHIC with GeV Adare:2008ab ().

In this paper we will scrutinise the thermal photon interpretation of the ALICE low- data with respect to theoretical uncertainties in the perturbative QCD prediction. Not only do these uncertainties arise from variations of the renormalisation and factorisation scales in the incoming and outgoing hadrons and photons, which have so far only been varied simultaneously by a factor of two about the photon , but they are also caused by uncertainties in the factorisation scheme Klasen:1996yk (), the nuclear parton density functions (nPDFs) and the photon fragmentation functions (FFs).

Nuclear PDFs are subject to intrinsic uncertainties of the proton PDFs, in particular of the gluon Stump:2003yu (), and cold nuclear modification effects, such as nuclear shadowing Hammon:1998gq (). Both are expected to be significantly larger at the LHC with its small values of partonic momentum fractions than at RHIC and at low scales of GeV rather than at large scales. Extra -broadening (the Cronin effect) is expected to become small at high (RHIC and LHC) centre-of-mass energies and small Dumitru:2001jx (). The latter has been confirmed experimentally at RHIC by the measurement of direct photons in +Au reactions Adare:2012vn () and at the LHC via the comparison of charged hadron production in proton-proton and proton-nucleus collisions ALICE:2012mj (). In general, when going from proton-proton to nucleus-nucleus collisions, one must also account for isospin suppression of prompt-photon production due the presence of neutrons (126) and not only protons (82) in the lead nucleus. This is, however, mostly relevant in the valence-quark region at large . Saturation effects appear to set in only at very small values of . The influence of these effects has been analysed in several global fits of nPDFs Eskola:2009uj (); Hirai:2007sx (); deFlorian:2011fp (); Schienbein:2009kk ().

The photon fragmentation functions are also subject to considerable uncertainties, in particular the one for gluons. In the absence of sufficiently detailed experimental data from collisions, they are typically modeled on vector meson fragmentation Bourhis:1997yu (). As an alternative, it has been proposed to use slightly virtual photons decaying into low-mass lepton pairs, as these do not have fragmentation contributions and need not be isolated experimentally Berger:1998ev (); Berger:1999es (); Brandt:2013hoa ().

The remainder of this paper is organised as follows: In Sec. 2, we will study the relative importance of the different partonic subprocesses contributing to prompt photon production in lead-lead collisions. Sections 3 and 4 give a brief review of the current uncertainty on the nuclear PDFs and photon FFs. Our main results are presented in Sec. 5, where we compare our NLO QCD predictions to the ALICE data and estimate the total theoretical uncertainty coming from the various sources defined above. Our conclusions are given in Sec. 6. For future use, we present in App. LABEL:app:a our NLO QCD calculations for the invariant prompt-photon yields in tabular form.

## 2 Partonic subprocesses

As stated in the introduction, the goal of this paper is to test how robust a thermal photon interpretation of the low- ALICE data is. We will explore to what extent the observed excess of photons with transverse momenta below 4 GeV can be explained by the production of hard photons in the initial partonic collision alone, taking into account all theoretical uncertainties.

The production of these photons is calculated perturbatively in NLO QCD using JETPHOX 1.2 Aurenche:2006vj (). As is well known, the truncation of the perturbative series at this order leads to an artificial dependence on the unphysical renormalisation scale and the initial- and final-state factorisation scales and . These have so far been varied only simultaneously, and we will later vary them also individually by a factor of two around the central scale . In addition, uncertainties in the determination of the PDFs must also be taken into account, and in the case of heavy-ion collisions not only for the free-proton PDFs, but also for their nuclear modification, which is commonly parametrised as a modification factor of the density for parton . Furthermore, at low transverse momenta and in the absence of a photon isolation criterion, fragmentation processes are expected to dominate, thus introducing an additional uncertainty from the photon fragmentation function.

Before we address these sources of uncertainty in the following chapters, it is illustrative to learn about the relative contributions of the various partonic subprocesses leading to prompt-photon production in lead-lead collisions. These depend in particular on the momentum fractions and of the partons in the PDFs and FFs , but also on the factorisation scales and . Since we are interested in the low- region, we set these scales here to 2 GeV, which is already in the perturbative regime, but only slightly larger than the starting scale GeV Eskola:2009uj (); Hirai:2007sx (); deFlorian:2011fp (); Schienbein:2009kk (). The relevant momentum fractions in the initial state can be computed from the photon’s transverse momentum by (the true at central rapidity), leading to typical values of . Our simulations show that the momentum fractions in the final-state fragmentation range from () to 1 at (16) GeV. In this section, we use the best fit of EPS09 Eskola:2009uj () for the nPDFs and the BFG set II Bourhis:1997yu () for the FFs as our baseline.

In Fig. 1, we plot the relative contributions to the total NLO prompt-photon cross section in lead-lead collisions with TeV. Individually shown are the fragmentation (full) and direct (dashed) subprocesses with initial gluons as well as those with initial quarks and antiquarks only (dotted and dot-dashed). As expected, fragmentation contributions dominate at low transverse momenta (up to 12 GeV), i.e. over almost the entire -range measured by ALICE. Only at larger transverse momenta the direct QCD “Compton” process takes over. In the initial state, the gluon contribution dominates in fact throughout, while pure quark and antiquark initial states make up for at most ten percent. We therefore expect the uncertainty on the gluon in bound nucleons at small to play an important role.

Looking at the final state, we plot in Fig. 2

the subprocesses with gluon (full) and quark (dashed) fragmentation as well as the direct subprocesses (dotted). As can be seen, the very low- region (below 4 GeV) is dominated by gluon fragmentation into photons, which is very poorly constrained experimentally. Quark fragmentation gives an almost constant contribution of about 40% over the entire -range. Direct photons become larger than quark fragmentation at transverse momenta of 8 GeV and than all fragmentation processes above 12 GeV (see above).

## 3 Nuclear parton density uncertainties

In this paper, we use the EPS09 nPDFs as our baseline Eskola:2009uj (). There, the bound-proton PDFs are defined in terms of nuclear modifications applied to CTEQ6.1M Stump:2003yu () free-proton PDFs ,

(1) |

The bound-neutron PDFs are obtained by assuming isospin symmetry. E.g. the total up-quark () distribution per nucleon in a nucleus with protons is

(2) |

where corresponds to the down-quark and the superscripts and refer to valence and sea quark contributions, respectively. The parametrisation of the nuclear modifications is performed at the charm quark mass () threshold imposing the momentum and baryon number sum rules

(3) |

for each nucleus separately. At higher scales, the nPDFs are obtained by solving the DGLAP evolution equations. This approach results in an excellent fit to the different types of nuclear hard-process data Eskola:2009uj (), suggesting that factorisation works well in the energy range studied and that the extracted nPDFs are universal in the region , GeV.

The theoretical uncertainties of the free-proton PDFs are obtained using the 40 error sets of the CTEQ6.1M parametrisation Stump:2003yu (). An additional 30 error sets are assigned pairwise to the uncorrelated eigendirections of the 15 parameters fitted to the nuclear collision data sets. A total uncertainty band at 90% confidence level is then calculated from the 71 sets defined by fixing either to the best fit value and varying the free-proton PDFs or fixing the latter to its best fit value and varying the former. These variations then contribute pairwise to the size of the upper and lower errors via

(4) | |||||

(5) |

As the authors acknowledge, this factorised approach represents a simplification, violating, e.g., in some cases momentum conservation, so that strictly speaking the free- and bound-proton PDF uncertainty analyses should not be separated Eskola:2009uj ().

Note that the minimal value of constrained by experimental data in the EPS09 fit lies above the value of expected to be relevant for the low- ALICE data. The nPDFs therefore rely in this region on an extrapolation from higher -values. Experimental data down to exist from the BRAHMS collaboration Arsene:2004ux (), who have measured charged-hadron production at forward rapidity and 3.2 in and +Au collisions at GeV. However, they have not been included in any of the nPDFs and would in fact lead to a poor description by the EPS09 fits, in particular for transverse momenta below 2 GeV at and below 4 GeV at (see Fig. 14 in Eskola:2009uj ()).

In order to estimate the bias from different underlying free-proton PDFs, parametrisations of the nuclear modification, and fitted nuclear data sets, we also study the best fits of the HKN07 Hirai:2007sx (), DSSZ deFlorian:2011fp () and nCTEQ Schienbein:2009kk () collaborations. In particular, the HKN07 nPDFs are based on the MRST1998 free-proton set Martin:1998sq () and those of DSSZ on the more recent MSTW2008 set Martin:2009iq (). The nCTEQ parametrisation is (so far) the only one that does not rely on a factorisation into a nuclear modification factor and free-proton PDFs. Instead it introduces an explicit -dependence in the coefficients of the -dependent functional form of the PDFs at the starting scale. So only the technical framework of the CTEQ6M analysis is used here Pumplin:2002vw ().

In Fig. 3, we plot the total parton densities defined in Eqs. (1) and (2) for protons bound in lead ions and as obtained in the EPS09 (full) Eskola:2009uj (), HKN07 (dashed) Hirai:2007sx (), DSSZ (dotted) deFlorian:2011fp () and nCTEQ (dot-dashed) Schienbein:2009kk () fits. As it is usually done, we plot times the nPDF, i.e. the momentum distribution of the partons in the bound proton. Error bands are only shown for the EPS09 analysis. In descending order of importance, the gluon, up- and strange-quark densities are shown separately. Those for down quarks differ only very little from the up-quark densities, mostly in the valence (large-) region, and are not shown separately. Whereas the up- (and down-) quark densities show the characteristic valence behavior and dominate at large , the gluon density becomes already dominant for -values below and in particular by more than an order of magnitude in the low- region around relevant for the prompt-photon data at low . The uncertainty estimated by EPS09 amounts there to almost a factor of two. Our calculations show that the uncertainty due to the proton PDF alone accounts for less than half of the uncertainty. If also the other parametrisations are taken into account in order to eliminate the theoretical bias, the uncertainty increases to more than a factor of two, due in particular to the different shapes of the HKN07 and nCTEQ gluon fits. We stress again that we are considering here a region of extrapolation, where no experimental constraints are taken into account. In the sea-quark region at low , the up- and strange-quark densities follow roughly the shape of the gluon as expected, but are smaller by about an order of magnitude. The uncertainties for the up quark (and down quark) are somewhat smaller than those of the gluon. This is probably due to better constraints on the large- (valence) region, which influence the low- region through the sum rules. However, the uncertainty for the strange quark amounts again to at least a factor of two. We therefore expect a considerable impact of the nPDF uncertainties, in particular from the gluon and from the different parametrisations, on the description of the low- prompt photon data.

## 4 Fragmentation function uncertainties

Fragmentation photons typically contribute less than 20% to the total yield in fixed-target experiments, but can become dominant at collider energies, in particular at low transverse momenta and in the absence of isolation criteria, as shown in Sec. 2. Note that the removal of any isolation criterion is imperative in the search for thermal photon sources in heavy-ion collisions. The hadronic input needed for the non-perturbative fragmentation functions would best be determined from inclusive photon production in annihilation, but unfortunately the experimental data are very limited and are furthermore dominated by the pointlike quark-photon fragmentation function. Therefore, all existing parametrisations use vector-meson dominance to model the photon fragmentation at low scales. The two most recent BFG parametrisations have, e.g., been obtained from ALEPH and HRS data on production Bourhis:1997yu (). They leave a large uncertainty in the gluon fragmentation to photons, which is parametrised with the two sets I and II. Heavy quarks are included above their production thresholds with zero boundary conditions.

Fig. 4 shows the FFs as obtained in the two BFG sets. Again, we show the momentum distribution, i.e. times the FF with GeV as the factorisation scale. In descending size, the up-, down- and strange-quark FFs are shown separately, together with the one for the gluon. Due to its larger charge (), which enters quadratically, the up-quark contribution exceeds the one for down and strange quarks with charge by about a factor of four. The difference between down and strange quarks is due to quark mass effects. The gluon is very weakly constrained (essentially only by the momentum sum rule), and its uncertainty amounts to almost an order of magnitude. Only when comparing the ratio of gluon and up-quark FFs into pseudoscalar mesons to the one into vector mesons and assuming that they should be of the same order due to the reduction of non-perturbative effects, BFG come to the conclusion that set II may be preferable. As is customary, we have followed BFG in this choice for our baseline, but stress again that the uncertainty on the gluon FF into photons remains large.

Despite the fact that photons are colour-neutral objects, they may also be affected by final-state effects, similarly, but perhaps not quite as strongly as mesons. In particular, the fragmentation of partons into photons is expected to be sensitive to the presence of a strongly interacting medium, since the final partons produced in the hard collision undergo multiple scattering and suffer energy loss from multiple-gluon radiation during the fragmentation process. This is usually modeled by rescaling the energy in a medium-modified fragmentation function Arleo:2006xb ()

(6) |

where and where denotes the probability that the leading parton with energy has lost an energy in the medium. The soft-gluon emission is usually assumed to follow a Poisson distribution. Alternatively, an improved Double Log Approximation (DLA) Albino:2009hu (), the Modified Leading-Log Approximation (MLLA) Borghini:2005em (), or parton shower and hadronisation models in Monte Carlo generators can be employed.

One would expect the modification of the fragmentation function to be largest for low values of the momentum fraction , where its hadronic component is most important, just as the photon structure function is dominated by its hadronic component at low Albino:2002ck (). As discussed above, the values of probed in the ALICE experiment range from () to 1 at (16) GeV, i.e. we expect the largest modification at low values of . This is also confirmed by a simple argument about the fragmentation time, which can be estimated with fm for GeV and thus becomes comparable to the formation time of the QGP 0.6 fm d'Enterria:2005vz () only at low values of .

Since nothing is known experimentally about light vector-meson (in particular ) fragmentation in the medium (and very little in vacuo, see above), one must resort to pseudoscalar (pion) fragmentation. It is then possible to simulate both initial- (isospin effects, shadowing) and final-state effects on inclusive photon production and compare, e.g., to PHENIX data Adler:2005ig (). The fragmentation modification is then found to be strongest at low values of (down to 4 GeV), but the experimental errors are still too large to identify it unambiguously. It has therefore been proposed to study not only inclusive photon production, but also photon-pion correlations to improve the separation of initial- and final-state effects Arleo:2006xb ().

## 5 Direct-photon production in heavy-ion collisions

In this section, we compute the transverse-momentum distribution of prompt (non-thermal) photons in NLO QCD and compare them to the experimental data for direct photons (not coming from decays) of the ALICE collaboration. Particular emphasis is put on the theoretical errors coming from the different sources described above, i.e. an individual variation of the renormalisation and factorisation scales by a factor of two up and down the central value , the nuclear PDFs from the EPS09 fits as described in Sec. 3, and the FFs from the BFG I and II fits as described in Sec. 4. Our goal is to establish to what extent the excess of the direct-photon data over the perturbative QCD prediction survives, once it is subtracted and all theoretical errors are taken into account, and how this procedure and these errors affect the determination of the effective temperature of the QGP.

In Fig. 5 we compare the transverse-momentum distribution of the ALICE data obtained in central (0-40%) lead-lead collisions at TeV (black data points Wilde:2012wc ()) to our NLO QCD predictions for prompt photons (red) as well as their fragmentation (blue) and direct (green) contributions. The centrality of the collision enters our calculations as an overall normalisation factor d'Enterria:2003qs ()

(7) |

which relates the invariant differential cross section (after integration over the azimuthal angle and for a finite rapidity interval ) to the invariant yield through d'Enterria:2003qs ()

(8) |

If we use for Eq. (7) the integration boundaries and averaging integrals of the ALICE experiment given in Tab. 1 of Ref. Abelev:2013qoq (), we find mb. If we then apply this normalisation to our theoretical predictions for the transverse-momentum spectrum in Fig. 5 (color) and compare to the experimental data (black) in the highest measured bin (8 GeV 11 GeV), we find very good agreement as expected in this perturbative region.

It is well known that the fragmentation contribution (blue) is considerably softer in than the direct contribution (green), i.e. the former falls off much more steeply than the latter. In Fig. 5, the fragmentation contribution dominates in all significant -bins, and in the lowest -bins it is more than an order of magnitude larger than the direct contribution. This leads to a large uncertainty from the fragmentation function. E.g., in the lowest -bins the theoretical prediction with BFG I FFs (magenta) is about a factor of three smaller than the one with BFG II FFs (blue). On the basis of BFG I FFs, one would thus attribute a larger excess to thermal photon production than with BFG II FFs. In the higher -bins, the FF uncertainty is smaller due to the larger contributions of both direct and quark-fragmentation photons (see Sec. 4).

The scale variation uncertainty (yellow) also turns out to be important. When all three scales are varied individually, the cross section changes by more than a factor of two up and down, much more than if all scales would be varied simultaneously. For low scales, where the cross section is larger, one would thus attribute a smaller excess to thermal photon production. Note that for GeV, we vary the factorisation scales only from to to stay above the nPDF and FF starting scales, which means that the true scale uncertainty there is even larger than the one shown. Only at larger values of , the scale uncertainty becomes eventually smaller. It could in principle be reduced by resumming the logarithmic parts of the higher-order corrections to all orders. This has been attempted, but unfortunately the corresponding codes for transverse-momentum Fink:2001bt () and joint Laenen:2000de (); Sterman:2004yk () resummation are not publicly available. The low- region is furthermore known to be sensitive to the Gaussian smearing parameters entering the non-perturbative factor required there Landry:2002ix ().

The uncertainties arising from the nPDFs (red, green and blue bands) are less important than the FF and scale errors, but are still well visible in Fig. 5. Their relative size remains constant over many low- bins, corresponding to low momentum fractions of the partons in the lead nucleus and the shadowing and free gluon density uncertainties present there, but eventually also becomes somewhat smaller at larger or . This was expected from Sec. 3, since the shadowing and gluon uncertainties are smaller and the quark contribution is larger at higher values of . As the nPDF fits by the other collaborations mentioned in Sec. 3 mostly fall wihin the EPS09 uncertainty band, we do not show the corresponding predictions here.

Considering statistical (vertical bars) and systematic (horizontal brackets) experimental as well as all theoretical errors described above individually, we find a photon excess (i.e. non-overlapping error bars) for GeV and the bin 2.6 GeV 3.0 GeV, i.e. at smaller -values than the ALICE collaboration, which did not take into account theoretical uncertainties, but used a larger theoretical normalisation (see Fig. 6 of Ref. Wilde:2012wc ()). Fitting this excess in the range 0.8 GeV 2.2 GeV with an exponential form , we find an inverse slope parameter of 58 MeV (304 51 GeV in Ref. Wilde:2012wc ()), where our (slightly larger) error has been obtained by adding statistical and systematic errors in quadrature and where both turn out to contribute with equal weight. Note, however, that these results have been obtained from data which include still both thermal and prompt photon contributions.

The thermal (plus medium-induced) photon production rate alone can be estimated by first subtracting from the experimental data the theoretically computed prompt photon contribution and then fitting again the exponential slope of the difference. This can, however, only be done for values larger than the starting scales of the nPDFs and FFs. In Fig. 6, the data points (full circles) now show the difference of data versus theory. The error bars have been obtained by propagating statistical and systematic experimental errors as well as scale and nPDF errors in quadrature. This leads to an additional photon excess in the bin 3.0 GeV 3.5 GeV. The (discrete) FF uncertainty is estimated by performing the same procedure with the (less favoured) BFG I fit (open circles) instead of the BFG II fit. Here, only the experimental errors were taken into account. Nevertheless, it falls within the total error of the other uncertainties.

Since the limitation of the previous fit to GeV is somewhat arbitrary, we perform two fits to the differential data points in Fig. 6. The first one for 1.5 GeV 2.2 GeV (blue) includes only data points from the overlap of the experimentally fitted and our theoretically calculated -bins. As the data point in the lowest -bin is somewhat higher than the three data points below and above it (see Fig. 5), this fit leads to a lower effective temperature of MeV ( MeV for BFG I FFs). Within its larger error, it is, however, still in agreement with the unsubtracted fit value of MeV obtained above. If we include all data points that show a non-zero difference from prompt-photon production, i.e. 1.5 GeV GeV, we find instead MeV ( MeV for BFG I FFs), which is again much closer to the unsubtracted fit value. The additional data points thus counterbalance the weight of the first one and bring the result into good agreement with the one based on Fig. 5.

This result may at first sight seem surprising and coincidental. One must, however, take into account the relative size of the prompt-photon contribution, which falls off like as predicted by our calculations, to the measured direct-photon rate in the different regions, shown in Fig. 7. Indeed, we find that prompt photons contribute less than 20% for GeV (less than 10% for BFG I FFs), i.e. their subtraction and the theoretical error do not strongly modify the exponential fall-off of the experimental data in this region. Above 4 GeV, we find instead an almost constant ratio of prompt over direct photons that is consistent with one within the uncertainties, i.e. these observed photons are produced in hard scatterings. In the intermediate region, it was already clear from Fig. 5 that the prompt-photon contribution is substantial, larger than 20% (10% for BFG I FFs), and must indeed first be subtracted from the experimental data. Only then can the excess be described again by an exponential fit, which we found to be in good agreement with the fit to the unsubtracted data at low .

Several attempts have been made to describe the transition from thermal to hard photons in the intermediate region theoretically. E.g., bremsstrahlung from a fast charged parton may be induced by the medium Zakharov:2004bi (), or a hard quark or gluon and a thermal parton may convert to a hard photon through annihilation or QCD Compton scattering Liu:2008zb (). As discussed above, we found that current data can well be described by an exponential fit up to GeV, once the prompt-photon contribution is subtracted. If we include the data points above values of 3.5 GeV, we obtain a very similar effective temperature of MeV, as these data points no longer carry much statistical weight. An almost identical effective temperature of MeV is also obtained when we vary all scales simultaneously instead of individually. This clearly demonstrates that the scale variation error is strongly correlated among different -bins. It affects mostly the normalisation constant rather than the inverse slope parameter .

## 6 Conclusions

In conclusion, we have performed in this work a first rigorous theoretical analysis of the ALICE measurement of direct-photon production in central lead-lead collisions with a centre-of-mass energy of TeV at low values of of 0.8 to 14 GeV. Based on a next-to-leading order QCD calculation, we found that prompt photon production in this region was largely induced by initial gluons and dominated by fragmentation contributions (for GeV those of a final gluon) due to the absence of an isolation criterion. This resulted in large theoretical uncertainties from independent variations of the renormalisation and factorisation scales, nuclear parton densities and fragmentation functions.

Nevertheless, we were able to confirm that the experimental data are well fitted in the region GeV by an exponential form with an effective temperature of 58 MeV, with a slightly larger uncertainty than the 51 MeV quoted by the ALICE collaboration. The reason is that in this region prompt photons contribute less than 10 to 20% to the total direct photon rate, so that their subtraction and theoretical error do not influence the fit result very much. We also verified that already for values of GeV the experimental data fall off with as predicted by perturbative QCD. In this region, the theoretical uncertainty is under reasonable control and strongly correlated among the different -bins, in particular the one coming from the scale variation.

In the intermediate -region from 1.5 to 3.5 GeV, the prompt-photon contribution had to be subtracted from the experimental data before a sensible exponential fit could be performed. We were able to verify an exponential fall-off with a very similar effective temperature of MeV. This result did not change significantly when extended to values of 7 GeV or when all scales were varied simultaneously, confirming the strong correlation of the theoretical error among different bins. Based on current ALICE data, a global thermal description thus seems to hold also in the intermediate region, and the presence of additional mechanisms, such as medium-induced bremsstrahlung and/or thermal parton conversion, do not need to be invoked to reproduce the measured photon spectra.

## Appendix A Prompt-photon yields and their uncertainties in NLO QCD

For future use, we present in Tab. 1 our NLO QCD calculations for the invariant prompt-photon yield in central (0-40%) lead-lead collisions at TeV and their uncertainties as a function of transverse momentum in tabular form. All numbers have been rescaled according to Eq. (8) using mb.

(GeV) | BFG set II | PDF error | Scale error | Fragm. | Direct | BFG set I | ||
---|---|---|---|---|---|---|---|---|

+ | - | + | - | |||||

1.5 - 1.6 | 2.68e-02 | 7.2e-03 | 8.4e-03 | 2.8e-02 | 1.6e-02 | 2.51e-02 | 1.71e-03 | 8.39e-03 |

1.6 - 1.8 | 1.98e-02 | 4.7e-03 | 6.4e-03 | 1.4e-02 | 1.1e-02 | 1.83e-02 | 1.51e-03 | 6.90e-03 |

1.8 - 2.0 | 1.19e-02 | 4.3e-03 | 2.8e-03 | 5.9e-03 | 6.2e-03 | 1.07e-02 | 1.19e-03 | 4.66e-03 |

2.0 - 2.2 | 7.77e-03 | 2.4e-03 | 1.7e-03 | 3.7e-03 | 3.8e-03 | 6.87e-03 | 8.99e-04 | 3.22e-03 |

2.2 - 2.4 | 4.96e-03 | 1.9e-03 | 6.8e-04 | 2.3e-03 | 2.0e-03 | 4.26e-03 | 7.05e-04 | 2.32e-03 |

2.4 - 2.6 | 3.98e-03 | 5.0e-04 | 2.4e-03 | 8.7e-04 | 2.0e-03 | 3.41e-03 | 5.68e-04 | 1.72e-03 |

2.6 - 3.0 | 2.20e-03 | 3.6e-04 | 7.8e-04 | 4.4e-04 | 8.5e-04 | 1.83e-03 | 3.66e-04 | 1.01e-03 |

3.0 - 3.5 | 1.04e-03 | 2.5e-04 | 2.7e-04 | 1.0e-03 | 7.0e-04 | 8.15e-04 | 2.27e-04 | 5.73e-04 |

3.5 - 4.0 | 4.94e-04 | 2.4e-04 | 1.5e-04 | 4.4e-04 | 2.9e-04 | 3.70e-04 | 1.23e-04 | 3.06e-04 |

4.0 - 4.5 | 2.87e-04 | 7.4e-05 | 3.6e-05 | 1.8e-04 | 1.5e-04 | 2.09e-04 | 7.80e-05 | 1.91e-04 |

4.5 - 5.0 | 1.72e-04 | 3.0e-05 | 2.5e-05 | 9.6e-05 | 9.0e-05 | 1.19e-04 | 5.22e-05 | 1.21e-04 |

5.0 - 5.5 | 1.09e-04 | 9.8e-06 | 4.0e-05 | 4.6e-05 | 5.2e-05 | 7.33e-05 | 3.53e-05 | 7.37e-05 |

5.5 - 6.0 | 6.58e-05 | 1.2e-05 | 1.5e-05 | 3.5e-05 | 2.8e-05 | 4.35e-05 | 2.23e-05 | 4.64e-05 |

6.0 - 6.5 | 4.65e-05 | 7.4e-06 | 5.6e-06 | 1.5e-05 | 2.0e-05 | 3.01e-05 | 1.64e-05 | 3.49e-05 |

6.5 - 7.0 | 3.05e-05 | 6.1e-06 | 3.1e-06 | 1.3e-05 | 1.0e-05 | 1.87e-05 | 1.18e-05 | 2.50e-05 |

7.0 - 7.5 | 2.28e-05 | 1.4e-06 | 7.7e-06 | 5.8e-06 | 8.6e-06 | 1.42e-05 | 8.55e-06 | 1.71e-05 |

7.5 - 8.0 | 1.50e-05 | 3.5e-06 | 1.5e-06 | 7.0e-06 | 4.4e-06 | 8.50e-06 | 6.54e-06 | 1.31e-05 |

8.0 - 8.5 | 1.21e-05 | 9.5e-07 | 2.3e-06 | 3.5e-06 | 3.8e-06 | 7.12e-06 | 4.94e-06 | 9.84e-06 |

8.5 - 9.0 | 8.74e-06 | 1.0e-06 | 1.5e-06 | 1.9e-06 | 2.7e-06 | 5.16e-06 | 3.59e-06 | 7.31e-06 |

9.0 - 9.5 | 6.47e-06 | 1.6e-06 | 9.8e-07 | 2.1e-06 | 1.8e-06 | 3.51e-06 | 2.97e-06 | 5.78e-06 |

9.5 - 10.0 | 5.00e-06 | 1.6e-06 | 3.8e-07 | 1.3e-06 | 1.4e-06 | 2.68e-06 | 2.31e-06 | 4.42e-06 |

10.0 - 10.5 | 4.10e-06 | 4.4e-07 | 7.0e-07 | 1.0e-06 | 1.1e-06 | 2.16e-06 | 1.94e-06 | 3.58e-06 |

10.5 - 11.0 | 3.13e-06 | 6.0e-07 | 2.9e-07 | 7.8e-07 | 7.1e-07 | 1.75e-06 | 1.38e-06 | 2.71e-06 |

11.0 - 11.5 | 2.48e-06 | 3.0e-07 | 3.1e-07 | 5.1e-07 | 6.5e-07 | 1.30e-06 | 1.17e-06 | 2.28e-06 |

11.5 - 12.0 | 2.00e-06 | 2.7e-07 | 3.3e-07 | 4.9e-07 | 5.0e-07 | 1.08e-06 | 9.21e-07 | 1.76e-06 |

12.0 - 12.5 | 1.63e-06 | 3.5e-07 | 3.1e-08 | 4.2e-07 | 3.5e-07 | 8.40e-07 | 7.92e-07 | 1.45e-06 |

12.5 - 13.0 | 1.36e-06 | 7.9e-08 | 2.4e-07 | 2.5e-07 | 3.1e-07 | 6.96e-07 | 6.61e-07 | 1.23e-06 |

13.0 - 13.5 | 1.13e-06 | 3.5e-08 | 2.9e-07 | 2.0e-07 | 3.0e-07 | 5.72e-07 | 5.53e-07 | 1.03e-06 |

13.5 - 14.0 | 8.89e-07 | 2.2e-07 | 1.3e-07 | 2.4e-07 | 1.8e-07 | 4.57e-07 | 4.32e-07 | 8.29e-07 |

14.0 - 14.5 | 7.78e-07 | 9.2e-08 | 8.5e-08 | 1.8e-07 | 1.7e-07 | 3.82e-07 | 3.96e-07 | 7.35e-07 |

14.5 - 15.0 | 6.51e-07 | 7.5e-08 | 6.6e-08 | 1.1e-07 | 1.5e-07 | 3.17e-07 | 3.34e-07 | 6.00e-07 |

15.0 - 15.5 | 5.09e-07 | 1.2e-07 | 2.8e-08 | 1.1e-07 | 8.0e-08 | 2.44e-07 | 2.65e-07 | 4.84e-07 |

15.5 - 16.0 | 4.69e-07 | 3.0e-08 | 1.3e-07 | 5.5e-08 | 1.1e-07 | 2.41e-07 | 2.28e-07 | 4.13e-07 |

The work of C. Klein-Bösing was supported by the Helmholtz Alliance Program of the Helmholtz Association, contract HA216/EMMI “Extremes of Density and Temperature: Cosmic Matter in the Laboratory”. We thank the ALICE collaboration for making their preliminary data available to us and J.P. Guillet, J. Owens, W. Vogelsang and M. Wilde for useful discussions.

## References

- (1) M. Wilde [ALICE Collaboration], Nucl. Phys. A904-905 2013 (2013) 573c [arXiv:1210.5958 [hep-ex]].
- (2) A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 104 (2010) 132301 [arXiv:0804.4168 [nucl-ex]].
- (3) D. G. d’Enterria and D. Peressounko, Eur. Phys. J. C 46 (2006) 451 [nucl-th/0503054] and references therein.
- (4) [ATLAS Collaboration], ATLAS-CONF-2012-051.
- (5) S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 710 (2012) 256 [arXiv:1201.3093 [nucl-ex]].
- (6) P. Steinberg [ATLAS Collaboration], Nucl. Phys. A 904-905 (2013) 233c.
- (7) R. Granier de Cassagnac [CMS Collaboration], Nucl. Phys. A 904-905 (2013) 241c.
- (8) P. Aurenche, M. Fontannaz, J. -P. Guillet, E. Pilon and M. Werlen, Phys. Rev. D 73 (2006) 094007 [hep-ph/0602133].
- (9) M. Klasen, Rev. Mod. Phys. 74 (2002) 1221 [hep-ph/0206169].
- (10) G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 706 (2011) 150 [arXiv:1108.0253 [hep-ex]].
- (11) [ATLAS Collaboration], ATLAS-CONF-2013-022.
- (12) V. Khachatryan et al. [CMS Collaboration], Phys. Rev. Lett. 106 (2011) 082001 [arXiv:1012.0799 [hep-ex]].
- (13) S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. D 84 (2011) 052011 [arXiv:1108.2044 [hep-ex]].
- (14) P. K. Khandai, P. Shukla and V. Singh, Phys. Rev. C 84 (2011) 054904 [arXiv:1110.3929 [hep-ph]].
- (15) M. Klasen and G. Kramer, Phys. Lett. B 386 (1996) 384 [hep-ph/9605210].
- (16) D. Stump, J. Huston, J. Pumplin, W. -K. Tung, H. L. Lai, S. Kuhlmann and J. F. Owens, JHEP 0310 (2003) 046 [hep-ph/0303013].
- (17) N. Hammon, A. Dumitru, H. Stoecker and W. Greiner, Phys. Rev. C 57 (1998) 3292 [hep-ph/9801292].
- (18) A. Dumitru, L. Frankfurt, L. Gerland, H. Stoecker and M. Strikman, Phys. Rev. C 64 (2001) 054909 [hep-ph/0103203].
- (19) A. Adare et al. [PHENIX Collaboration], arXiv:1208.1234 [nucl-ex].
- (20) B. Abelev et al. [ALICE Collaboration], Phys. Rev. Lett. 110 (2013) 082302 [arXiv:1210.4520 [nucl-ex]].
- (21) K. J. Eskola, H. Paukkunen and C. A. Salgado, JHEP 0904 (2009) 065 [arXiv:0902.4154 [hep-ph]].
- (22) M. Hirai, S. Kumano and T. -H. Nagai, Phys. Rev. C 76 (2007) 065207 [arXiv:0709.3038 [hep-ph]].
- (23) D. de Florian, R. Sassot, P. Zurita and M. Stratmann, Phys. Rev. D 85 (2012) 074028 [arXiv:1112.6324 [hep-ph]].
- (24) I. Schienbein, J. Y. Yu, K. Kovarik, C. Keppel, J. G. Morfin, F. Olness and J. F. Owens, Phys. Rev. D 80 (2009) 094004 [arXiv:0907.2357 [hep-ph]].
- (25) L. Bourhis, M. Fontannaz and J. P. Guillet, Eur. Phys. J. C 2 (1998) 529 [hep-ph/9704447].
- (26) E. L. Berger, L. E. Gordon and M. Klasen, Phys. Rev. D 58 (1998) 074012 [hep-ph/9803387].
- (27) E. L. Berger, L. E. Gordon and M. Klasen, Phys. Rev. D 62 (2000) 014014 [hep-ph/9909446].
- (28) M. Brandt and M. Klasen, arXiv:1305.5677 [hep-ph].
- (29) I. Arsene et al. [BRAHMS Collaboration], Phys. Rev. Lett. 93 (2004) 242303 [nucl-ex/0403005].
- (30) A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thorne, Eur. Phys. J. C 4 (1998) 463 [hep-ph/9803445].
- (31) A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. C 63 (2009) 189 [arXiv:0901.0002 [hep-ph]].
- (32) J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. M. Nadolsky and W. K. Tung, JHEP 0207 (2002) 012 [hep-ph/0201195].
- (33) F. Arleo, JHEP 0609 (2006) 015 [hep-ph/0601075].
- (34) S. Albino, B. A. Kniehl and R. Perez-Ramos, Nucl. Phys. B 819 (2009) 306 [arXiv:0903.4825 [hep-ph]].
- (35) N. Borghini and U. A. Wiedemann, hep-ph/0506218.
- (36) S. Albino, M. Klasen and S. Soldner-Rembold, Phys. Rev. Lett. 89 (2002) 122004 [hep-ph/0205069].
- (37) S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 94 (2005) 232301 [nucl-ex/0503003].
- (38) D. G. d’Enterria, nucl-ex/0302016.
- (39) B. Abelev et al. [ALICE Collaboration], arXiv:1301.4361 [nucl-ex].
- (40) C. E. Fink, hep-ph/0105276.
- (41) E. Laenen, G. F. Sterman and W. Vogelsang, Phys. Rev. Lett. 84 (2000) 4296 [hep-ph/0002078].
- (42) G. F. Sterman and W. Vogelsang, Phys. Rev. D 71 (2005) 014013 [hep-ph/0409234].
- (43) F. Landry, R. Brock, P. M. Nadolsky and C. P. Yuan, Phys. Rev. D 67 (2003) 073016 [hep-ph/0212159].
- (44) B. G. Zakharov, JETP Lett. 80 (2004) 1 [Pisma Zh. Eksp. Teor. Fiz. 80 (2004) 3] [hep-ph/0405101].
- (45) W. Liu and R. J. Fries, Phys. Rev. C 77 (2008) 054902 [arXiv:0801.0453 [nucl-th]].