Transverse momentum dependence of charmonium production in heavy ion collisions.
1 / 41

- PowerPoint PPT Presentation

  • Updated On :

Transverse momentum dependence of charmonium production in heavy ion collisions. N. Topilskaya, A.Kurepin – INR, Moscow . 3rd INT. WORKSHOP ON HIGH-PT PHYSICS AT LHC TOKAJ, HUNGARY March, 16-19, 2008 Tokaj , Hungary. Charmonium.

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
Download Presentation

PowerPoint Slideshow about '' - RexAlvis

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
Slide1 l.jpg

Transverse momentum dependence of charmonium production in heavy ion collisions.

N. Topilskaya, A.Kurepin – INR, Moscow


Slide2 l.jpg

Charmonium heavy ion collisions.

  • 33 years ago:discovery of J/ψ, 21 years ago: Matsui & Satz

  • colour screening in deconfined matter →J/ψ suppression

  • →possible signature of QGP formation

  • Experimental and theoretical progress since then → situation is much more complicated

    • cold nuclear matter / initial state effects

      • “normal” absorption in cold matter

      • (anti)shadowing

      • saturation, color glass condensate

    • suppression via comovers

    • feed down from cc, y’

    • sequential screening (first: cc, y’, J/y only well above Tc)

    • regeneration via statistical hadronization or charm coalescence

      • important for “large” charm yield, i.e. RHIC and LHC

Tokaj, N.Topilskaya, March 16-19, 2008

Slide3 l.jpg

NA50 experimental setup heavy ion collisions.

The J/ is detected via its decay into muon pairs

Dimuon spectrometer: Centrality detectors: EM calorimeter (1.1< lab<2.3)

2.92 < ylab< 3.92 ZDC calorimeter (lab> 6.3)

cos CS < 0.5 Multiplicity detector (1.9<lab<4.2)

Pb-Pb 158 GeV/c p – A 400 GeV/c 2000 year

Data period Subtargets Number of J/ Target Number of J/

1995 7 50000 Be 38000

1996 7 190000 Al 48000

1998 1 49000 Cu 45000

2000 1 in vacuum 129000 Ag 41000

W 49000

Pb 69000

J/y suppression is generally considered as one of the most direct signatures of QGP formation (Matsui-Satz 1986)

Slide4 l.jpg

Fit to the mass spectrum heavy ion collisions.

Slide5 l.jpg

  • Light systems and peripheral Pb-Pb collisions: heavy ion collisions.J/ψ is absorpted by nuclear matter . The scaling variable -L (length of nuclear matter crossed by the J/ψ)

  •  (J/ψ) ~ exp( -abs L)

  • Central Pb-Pb collisions:the L scaling is broken - anomalous suppression

J/ψ suppression from p-A to Pb-Pb collisions

J/ψ production has been extensively studied inp-A, S-UandPb-Pbcollisions by the NA38 and NA50 experiments at the CERN SPS




J/y normal nuclear absorption curve

NA60 : is anomalous suppression present also in lighter In-In

nuclear systems ?Scaling variable- L, Npart, ε ?

Slide6 l.jpg

MWPC’s heavy ion collisions.


~ 1m

Muon Spectrometer



Hadron absorber

Toroidal Magnet

Target area



Trigger Hodoscopes

Dipole field2.5 T





Matching in coordinate and in momentum space







not to scale

  • Origin of muons can be accurately determined

  • Improved dimuon mass resolution

allows studies vs.

collision centrality


NA60 experimental setup

High granularity and radiation-hard silicon tracking telescope in the vertex region before the absorber

Slide7 l.jpg

Comparison of J/ heavy ion collisions./DY results

An “anomalous suppression” is presented already in In-In

The normal absorption curve is based on NA50 results. Its uncertainty (~ 8%) at 158 GeV is dominated by the (model dependent) extrapolation from the 400 and 450 GeV

p-A data.

need p-A measurements at 158 GeV

Slide8 l.jpg

Direct J/ heavy ion collisions. in In-In

Data are compared with a theoretical J/ distribution, obtained within the Glauber model, taking into account the nuclear absorption.

The ratio Measured / Expected is normalized to thestandard analysis



Anomalous suppression begins in the range 80 < NPart < 100

Large systematic errors


Slide9 l.jpg

С heavy ion collisions.omparison J/resultsvesus Npart

NA50: Npart ftom Et (left) and from Ezdc (right, as in NA60)

J/ysuppression inIn-Inis in agreement withPb-Pb

S-Uhas different behaviour

Slide10 l.jpg

(J/)/DY = heavy ion collisions.29.2  2.3

L = 3.4 fm

Сomparison of J//DY

Preliminary NA60 results on p-A at 158 GeV show that rescaling from 400 and 450 GeV to 158 GeV is correct.

Results on abs will appear soon HP08- crucial to confirm (or modify) the anomalous suppression pattern

Slide11 l.jpg

Preliminary! heavy ion collisions.

’ suppression(NA38, NA50, NA60)

abs=8±1 mb

abs~20 mb

Small statistics in NA60 In-In for’ (~300)

The most peripheral point (Npart~60)– normal nuclearabsorption

Slide12 l.jpg

Suppression by produced hadrons (“comovers”) heavy ion collisions.

The model takes into account nuclear absorption and comovers interaction with σco = 0.65 mb (Capella-Ferreiro) EPJ C42(2005) 419

In-In 158 GeV

J/y / NColl

nuclear absorption

comover + nuclear absorption

(E. Ferreiro, private communication)

Pb-Pb 158 GeV

NA60 In-In 158 GeV

Slide13 l.jpg

QGP + hadrons + regeneration + in-medium effects heavy ion collisions.

The model simultaneously takes into account dissociation and regeneration processes in both QGP and hadron gas (Grandchamp, Rapp, Brown EPJ C43 (2005) 91)

In-In 158 GeV

fixed thermalization time

centrality dependent thermalization time


Nuclear Absorption

Suppression + Regeneration

QGP+hadronic suppression


Number of participants

Pb-Pb 158 GeV

centrality dependent thermalization time

fixed thermalization time

NA60 In-In 158 GeV

Slide14 l.jpg

Suppression due to a percolation phase transition heavy ion collisions.

Model based on percolation (Digal-Fortunato-Satz)

Eur.Phys.J.C32 (2004) 547.

Prediction: sharp onset (due to the disappearance of the cc meson) at Npart ~ 125 for Pb-Pb and

~ 140 for In-In

Pb-Pb 158 GeV

NA60 In-In 158 GeV

The dashed line includes the smearing due to the resolution

Slide15 l.jpg

Maximal hadronic heavy ion collisions.absorption

  • Comparison J/ production

  • with calculations

    • nuclear absorption---

    • maximal possible __

    • absorption in a hadron

    • gas(T = 180 MeV)

L. Maiani et al.,

Nucl.Phys. A748(2005) 209

F. Becattini et al.,Phys. Lett. B632(2006) 233

l –transverse size of fire-ball

  • Pb-Pb and In-In (in lower order)showextra suppression

Slide16 l.jpg

Comparison of experimental SPS data. heavy ion collisions.


J/ and - normal nuclear absorption


J/ - normal nuclear absorption

 - anomalous suppression


J/ - onset of anomalous suppression

- anomalous suppression ~ S-U


J/ - onset of anomalous suppression

- anomalous suppression < S-U

Open question: S-U vs In-In ? Theoretical description?

Slide17 l.jpg

J/ heavy ion collisions.transverse momentum distribution

Study <pT2> and T dependence on centrality

NA60 In-In

Slide18 l.jpg

J/ heavy ion collisions.transverse momentum distribution

NA50 and NA38

Fitting: <pT2>(L) = <pT2>pp + αgN L

Simultanious fit with an energy dependent pT2pp and

a common slope:

gN= 0.081±0.002 (GeV/c)2/fm-1

Then model dependent

extrapolation of all

data to 158 GeV

Slide19 l.jpg

J/ heavy ion collisions.transverse momentum distribution

<pT2> versus L

Fitting: <pT2>(L) = <pT2>pp + αgN L

<pT2>pp= 1.08 ± 0.02 GeV2/c2

χ2= 0.85

 αgN = 0.083 ± 0.002 GeV2/c2fm-1

The observed dependence could simply

result from parton initial state multiple


Slide20 l.jpg

J/ heavy ion collisions.transverse momentum distribution in p-A

<pT2> versus L

NA60 p-A at 158 GeV/c-

the same energy and kinematical domain as

Pb-Pb and In-In

New 158 GeV/c data show that at SPS gN

depends on

theenergy of the collision

Slide21 l.jpg

J/ heavy ion collisions.transverse momentum distribution in p-A and A-A

<pT2> versus L

NA60p-A and In-In and

NA50 Pb-Pb - at 158 GeV

and in the same kinematical


  • pT2 increases linearly with L in both p-A, In-In and Pb-Pb

  • However, the scaling of pT2 with L is broken moving from p-A to A-A

  • On one hand comparing p-A and peripheral In-In the suppression scales with L

  • On the other hand the J/ pT distributions do not scale with L

Slide22 l.jpg

NA50 and NA38 heavy ion collisions.Teff rescalculated to

158 GeV vs energy density

T(=0) =( 182)2 MeV

Tslope = ( 20.16  1.04)  10-3 fm3

Tslope(cent Pb-Pb)=(8.87  2.07) 10-3 fm3

R(slopes)=2.27 +/- 0.54

InNA38 and NA50 TJ/ ψ grows linearly with the energy density and with L.

Model dependent recalculation 400 and 200 GeV data to 158 GeV- scaling.

For the most central Pb-Pb collisions more flat behaviour could be seen.

Slide23 l.jpg

С heavy ion collisions.omparisonT(J/ψ) at 158 GeV

Fitting functions:

dN/dMT ~ MT2K1(MT/T) – NA50

dN/dMT ~ MT exp(-MT/T) – NA60 –

gives slightly lower temperature ~ 7 MeV

Fitting functions

No scaling with L for p-A and A-A

Slide24 l.jpg

J/ψ suppression versus p heavy ion collisions.T.

F=(J//DY>4.2 )acc vs pT in 5 ET bins


NA50 Pb-Pb 2000


Et bins in GeV

1. 5 - 20

2. 20 - 40

3. 40 - 70

4. 70 - 100

5. >100


Slide25 l.jpg

J/ψ suppression versus E heavy ion collisions.T.

F=(J//DY>4.2 )acc vs ET in 11 pT bins

5 Et bins

NA50 Pb-Pb 2000

log scale

Clear centrality dependence for low pt.

Much weaker dependence for high pt.

Slide26 l.jpg

R heavy ion collisions.cp = (J/ψi(pT)/DYi>4.2)/(J/ψ1(pT)/DY1>4.2)

Pb-Pb NA50

5 Et bins

The ratios to the most peripheral E 1 bin.

The suppression vs the most peripheral events is significant mainly at low pT where it strongly increases with centrality. For central events the suppression exists over the whole pT range.

Slide27 l.jpg

Suppression vs p heavy ion collisions.T for p-A, S-U and Pb-Pb





Cronin effect- enhancement at pT>2 GeV/c

Pb-Pb 2000


Et bins GeV

5 - 40

40 - 80

80 – 125

Slide28 l.jpg

R heavy ion collisions.AA


1.5-5 %







pT (GeV/c)

NA60 In-In

Nuclear modification factor RAA=NAA/(Npp*<Ncoll>)

J/ pT distribution for pp was calculated in the form 1/pT dN/dpT ~ MTK1(MT/T) – systematic error 11%

Enhancement (Cronin effect) at pT > 2 GeV/c

Slide29 l.jpg

1.5-5% heavy ion collisions.








pT (GeV/c)

Rcp vs pT.

NA60 In-In

Rcp = (J/ψi(pT)/Ncolli)/(J/ψ1(pT)/Ncoll1)

The ratios to the peripheral i=1 (47-57%)bin.

Large suppression at low pT, growing with centrality- as in RAA NA60

and in Rcp NA50.

Slide30 l.jpg

Summary for SPS data heavy ion collisions.

  • The J/y shows an anomalous suppression discovered in Pb-Pb

  • and existing already in In-In

  • None of the available models properly describes the observed suppression pattern simultaneously in Pb-Pb and In-In

  • The transverse momentum dependence for p-A and A-A at 158 GeV

  • shows no L scaling in <pT2>

    The suppression in Pb-Pb and In-In is significant mainly at low

    pT where it strongly increases with centrality.

    For central events the Rcp suppression exists over the whole

    pT range in Pb-Pb and In-In.

    In p-A, S-U, peripheral Pb-Pb events and in RAA In-In the

    enhancement for pT> 2 GeV (Cronin effect) is seen.

Slide31 l.jpg

J/ heavy ion collisions. in PHENIX

J/  e+e–

identified in RICH and EMCal

  • |y| < 0.35

  • Pe > 0.2 GeV/c

  •  = 


    identified in 2 fwd spectrometers

    South :

    • -2.2 < y < -1.2

      North :

    • 1.2 < y < 2.4

  • P > 2 GeV/c

  •  = 2 

    Event centrality and vertex given by

    BBC in 3<||<3.9 (+ZDC)

Centrality is calculated to Npart (Ncoll) using Glauber model

Slide32 l.jpg

Yan, Zhuang, Xu heavy ion collisions.


All models for y=0







Suppression RAA vs Npart at RHIC.

PHENIX Au-Au data

Models for mid-rapidity Au-Au data

Without regeneration

With regeneration

Slide33 l.jpg

( heavy ion collisions.dN/dy)AuAu

(dN/dy)pp x<Ncol>


Suppression RAA vs Npart at RHIC.

Au+Au: A. Adare et al. (PHENIX) PRL 98 232301 (2007)

Cu+Cu: A. Adare et al. (PHENIX) arXiv:0801.0220

  • Cold Nuclear Matter (CNM) effects

    • Nuclear absorption

    • Gluons shadowing

  • Evaluated from J/ψ production

  • in d+Au collisions.

  • A.Adare et al. (PHENIX) arXiv:0711.3917





  • J/y suppression at mid-rapidity

  • at RHIC is compatible to

  • CNM effects except most central Au+Au collisions.

  • Stronger suppression at forward

  • rapidity than CNM effects.





Slide34 l.jpg

J heavy ion collisions./ψ suppression (SPS and RHIC)

J/ψ yield vs Npart,

normalized on Ncoll.

Unexpected good scaling.

Coherent interpretation-

problem for theory.

Work start - : Karsch, Kharzeev and Satz., PRL637(2006)75

Slide35 l.jpg

arXiv:0801.0220 [nucl-ex] heavy ion collisions.

PHENIX invariant cross sections of J/y

J/y was measured from pT=0GeV/c to beyond pT =5GeV/c.

Slide36 l.jpg

arXiv:0801.0220 [nucl-ex] heavy ion collisions.

J/ψ suppression RAA vs pT at PHENIX.




For low pT suppression grows with centrality.

Slide37 l.jpg

Comparison SPS (NA60) and RHIC (PHENIX) data heavy ion collisions.

The same suppression at

low pT.

Larger values of <pT2> at


Slide38 l.jpg

Suppression R heavy ion collisions.AA in Au-Au (PHENIX) vs pT.


J/ψ up to only 5 GeV

Central events

The same RAA for

0,  at all pT

and J/ (up to 4 GeV/c).

RAA for  is higher.

RAA for direct  <1 for

high pT.

Slide39 l.jpg

J/ψ suppression R heavy ion collisions.AA at RHIC.

PHENIX and STAR Cu-Cu data

  • Data consistent with no suppression at

    high pT: RAA(pT > 5 GeV/c) = 0.9 ± 0.2

  • At low-pT RAA: 0.5—0.6 (PHENIX)

  • RAA increase from low pT to high pT

  • Most models expect a decrease RAA at

    high pT:

    X. Zhao and R. Rapp, hep-ph/07122407

    H. Liu, K. Rajagopal and U.A. Wiedemann,

    PRL 98, 182301(2007) and hep-ph/0607062

  •  But some models predict an increase RAA

  • at high pT:

  • K.Karch and R.Petronzio, 193(1987105;

  • J.P.Blaizot and J.Y.Ollitrault, PRL (1987)499

Slide40 l.jpg

Conclusions heavy ion collisions.

  • At SPS energiesthe J/y shows an anomalous suppression discovered in Pb-Pb and existing already in In-In

  • None of the available models properly describes the observed suppression pattern simultaneously in Pb-Pb and In-In

  • The  shows an anomalous suppression for S-U, In-In

  • and Pb-Pb

  • At RHIC energies the J/suppression is of the same order as at SPS

  • None of the theoretical model could describe all the data

  • The transverse momentum dependence of J/ψ suppression shows

  • suppression mainly ay low pT, growing with centrality

  • Need information at high pT.

Slide41 l.jpg

Hope heavy ion collisions.- measurement at LHC with high values of energy density

and transverse momentum pT.

Need- high statistic pp, p-A and A-A data at the same conditions.

Work for theory.