1 / 55

P article correlations at STAR

P article correlations at STAR. Jan Pluta Heavy Ion Reactions Group (HIRG), Faculty of Physics, Warsaw University of Technology. Some results from the STAR HBT group, presented recently by: Z.Chajecki, A.Kisiel, M.Lisa, M.Lopez-Noriega, S.Panitkin, F.Retiere, P.Szarwas.

jabari
Download Presentation

P article correlations at STAR

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. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Particle correlationsat STAR Jan Pluta Heavy Ion Reactions Group (HIRG), Faculty of Physics, Warsaw University of Technology Some results from the STAR HBT group, presented recently by: Z.Chajecki, A.Kisiel, M.Lisa, M.Lopez-Noriega,S.Panitkin, F.Retiere, P.Szarwas. 3-rd Budapest Winter School on Heavy Ion Physics, 10 XII 2003

  2. Outline: • The STAR experiment • RHIC HBT Puzzle • General analysis • asHBT • Two K0-short correlations • p-Au, d-Au data • Nonidentical particles - emission asymmetry • Plans for future

  3. 12:00 o’clock BRAHMS PHOBOS 2:00 o’clock 10:00 o’clock RHIC PHENIX 8:00 o’clock 4:00 o’clock STAR 6:00 o’clock 9 GeV/u Q = +79 U-line BAF (NASA) m g-2 LINAC BOOSTER HEP/NP AGS 1 MeV/u Q = +32 TANDEMS Relativistic Heavy Ion Collider (RHIC) =3.8 km1740 superconducting magnets • Beam energy up to 100 GeV/A (250 GeV for p) • Two independent rings (asymmetric beam collisions are possible) • Beam species: from p to Au • Six interaction points • Four experiments:STAR, PHENIX, PHOBOS and BRAHMS

  4. Solenoidal Tracker At RHIC

  5. STAR Detector – side view STAR Detector – side view and STAR Collaboration – face view

  6. STAR Collaboration • 500 Collaborators including • ~65 graduate students • ~60 postdocs • 12 countries • 49 institutions • Spokesperson: John Harris 1991 - 2002 Tim Hallman 2002 - now USA, Brazil, China, Croatia Czech Republich, England, France, Germany, India, Netherlands, Poland, Rusia

  7. HBT+FSI Particle correlations The idea: Quantum statistics and Final-State Interaction Space-time sizes and dynamics Correlation function Momenta and momentum difference

  8. Central Event: AuAu 200GeV/A STAR Event 2 4m Real-time track reconstruction Pictures from Level 3 Trigger, online display. Typically 1000 to 2000 tracks per event into the TPC

  9. Minbias trigger STAR PRELIMINARY N ch Event and Particle Selection Au+Au Collisions at Sqrt(SNN)=200GeV • Centrality selection based on number of charged hadrons. three different centralities • Midrapidity -0.5 < y < 0.5 • Particle identification via specific ionization (dE/dx) electron band removed by cuts • Optimum performance for p HBT: 0.150 < pT (GeV/c) < 0.550 for K0s: 0.100 < pT (GeV/c) < 3.500 STAR PRELIMINARY

  10. Some base definitions - to be used for results presentation Two-particle kinematics Main (approximative) relations: Qout <--> Pt Qside <-->  Qlong <-->  KT <--> Pt LCMS: (P1+P2)z=0

  11. HBT Excitation Function Comparison with lower energies for ~ 10% most central events at midrapidity kT ~ 0.17 GeV/c No significant increase in radii with energy RO/RS ~ 1 Gap in energy that needs to be closed

  12. STAR 130 GeV PHENIX 130 GeV p + p - RHIC HBT Puzzle Most “reasonable” models still do not reproduce RHIC sqrt(SNN) = 130GeV HBT radii Hydro + RQMD PHENIX PRL 88 192302 (2002) √SNN = 130GeV • “Blast wave” parameterization (Sollfrank model) can approximately describe data • …but emission duration must be small •  = 0.6 (radial flow) • T = 110 MeV • R = 13.5  1fm (hard-sphere) • emission= 1.5  1fm/c (Gaussian) from spectra, v2

  13. raw Coulomb corrected q* (GeV/c) Sqrt(SNN) = 200 GeV 3Dimensional Pion HBT p –p– Pratt-Bertsch Parameterization LCMS frame: (p1+p2)z=0 Central Events pT = 0.15-0.25 GeV/c Coulomb correction → spherical Gaussian source of 5fm momentum resolution corrected(~1% effect at 200GeV, due to higher B-field) STAR PRELIMINARY Statistical errors only!!

  14. 200GeV Central Midcentral Peripheral Centrality and mT dependence at 200 GeV STAR PRELIMINARY RL varies similar to RO, RS with centrality HBT radii decrease with mT (flow) Roughly parallel mT dependence for different centralities RO/RS ~ 1 (short emission time)

  15. 200GeV - 130 GeV Central Midcentral Peripheral PHENIX Central Comparison: 200 to 130 GeV. Longitudinal radius (fit to STAR Y2 data only) STAR PRELIMINARY Longitudinal radius: at 200GeV identical to 130 GeV

  16. 200GeV - 130 GeV Central Midcentral Peripheral PHENIX Central Evolution timescale from RL Simple Mahklin/Sinyukov fit (assuming boost-invariant longitudinal flow) Makhlin and Sinyukov, Z. Phys. C 39 (1988) 69 STAR PRELIMINARY Assuming TK=110 MeV(from spectra at 130 GeV) (fit to STAR Y2 data only)

  17. 200GeV - 130 GeV Central Midcentral Peripheral PHENIX Central Comparison: 200 to 130 GeV. Transverse radii STAR PRELIMINARY * • Higher B-field  higher pT • Transverse radii: • similar but not identical • low-pT RO, RS larger at 200 GeV • steeper falloff in mT • (PHENIX 130GeV) • Ro falls steeper with mT

  18. P. Kolb and U. Heinz, hep-ph/0204061 Azimuthally sensitive HBT (asHBT) • sensitive to interplay b/t anisotropic geometry & dynamics/evolution • “broken symmetry” for b0 => more detailed, important physics information • another handle on dynamical timescales – likely impt in HBT puzzle

  19. side y K out Lines: projections of 3D Gaussian fit fp x b HBT respect Reaction Plane 1D projections, f=45° √SNN = 130GeV

  20. Data corrected for both event plane resolution and merging systematic HBT(φ) Results – 130 GeV T=100 MeV r0=0.6 ra=0.037, R=11.7 fm, t=2.2 fm/c Minbias events @130GeV Bolstered statistics by summing results of p- and p+ analyses Blast-wave calculation (lines) indicates out-of-plane extended source Star preliminary

  21. RY 0 RX  0 • BW: hydro-inspired parameterization • of freezeout • longitudinal direction • infinite extent geometrically • boost-invariant longitudinal flow • Momentum space • temperature T • transverse rapidity boost ~ r A model of the freezeout - BlastWave • coordinate space • transverse extents RX, RY • freezeout in proper time  • evolution duration 0 • emission duration 

  22. RY RX • BW: hydro-inspired parameterization • of freezeout • Longitudinal direction • infinite extent geometrically • boost-invariant longitudinal flow • Momentum space • temperature T • transverse rapidity boost ~ r A model of the freezeout- BlastWave • Coordinate space • transverse extents RX, RY • freezeout in proper time  • evolution duration 0 • emission duration  7 parameters describing freezeout

  23. Central Midcentral Peripheral T (MeV) 108  3 106  3 95  4 0 0.88  0.01 0.87  0.02 0.81  0.02 a 0.06  0.01 0.05  0.01 0.04  0.01 RX (fm) 12.9  0.3 10.2  0.5 8.0  0.4 RY (fm) 12.8  0.3 11.8  0.6 10.1  0.4 0 (fm/c) 8.9  0.3 7.4  1.2 6.5  0.8  (fm/c) 0.0  1.4 0.8  3.2 0.8  1.9 2 / ndf 80.5 / 101 153.7 / 92 74.3 / 68 central midcentral peripheral BlastWave fits to published RHIC data • reasonable centrality evolution • OOP extended source in non-central collisions

  24. FO = INIT Estimate of initial vs F.O. source shape • estimate INIT from Glauber • from asHBT: • FO < INIT→ dynamic expansion • FO > 1 → source always OOP-extended • constraint on evolution time

  25. asHBT at 200 GeV in STAR – R() vs centrality • 12 (!) -bins b/t 0-180 (kT-integrated) • 72 independent CF’s • clear oscillations observed in transverse radii of symmetry-allowed* type • Ro2, Rs2, Rl2 ~ cos(2) • Ros2 ~ sin(2) • centrality dependence reasonable • oscillation amps higher than 2nd-order ~ 0→ (*) Heinz, Hummel, MAL, Wiedemann, Phys. Rev. C66 044903 (2002)

  26. Pion Correlations d-Au and p-Au • Pion correlation in d – Au : data selection • p-Au selection • 1D correlation function • 3D correlation function • d-Au vs p-Au • KT dependence • Centrality dependence

  27. ZDC-d ZDC-Au p-Au selection: FTPC E -Au Au d All trigger events Using information from ZDC-d STAR can separate events with neutron spectator from deuteron

  28. theoretical CF: Rinv=6 fm, l = 0.5 d*-Au : d-Au without p-Au 1D Correlation Function: Gaussian fit: STAR preliminary • CF is very wide (rel Au-Au) • Coulomb/merging less important • CF looks reasonable • 1D Gaussian fit is not good • needed more deeply study of fit method only statistical error included !

  29. 3D Correlation Function: 3D Gaussian fit: Gaussian parametrization is not perfect but HBT radii characterize the width of CF STAR preliminary Fit results: Rout, Rside sensitive to the number of participants cut on the others Q's components < 30 MeV/c [GeV/c]

  30. KT dependence: STAR preliminary • clear KT dependence • Rout and Rside - sensitive to the number of participants • Rlong – the same KT dependence for dAu and pAu p – Au d – Au

  31. KT dependence: d-Au & Au-Au divided by p-p STAR preliminary • the same trend of KT dependence for d-Au and Au-Au as for p-p • HBT radii are scaled by constant factors

  32. MT dependence of Rlong: STAR preliminary Sinyukov fit: Rlong= const (mT)-a mT = kT2 + massp2 a for different collisions STAR preliminary p-p d-Au Au-Au Au-Au peripheral midcentral

  33. ZDC-d ZDC-Au Centrality definition in d-Au: FTPC-Au: charged primary particle multiplicity in -3.8<<-2.8 most peripheral most central 1 2 3 FTPC E -Au Au d

  34. centrality Centrality dependence: • clear centrality dependence • similar to AuAu • connection to geometry STAR preliminary p – Au d – Au minbias [*] - Glauber calculations (Mike Miller)

  35. K0sK0s Correlations

  36. mt scaling violation? inv Next RHIC HBT puzzle ?

  37. Non-identical particle correlations: The asymmetry analysis Catching up • Interaction time larger • Stronger correlation Moving away • Interaction time smaller • Weaker correlation “Double” ratio • Sensitive to the space-time asymmetry in the emission Kinematics selection on any variable e.g. kOut, kSide, cos(v,k) R.Lednicky, V. L.Lyuboshitz, B.Erazmus, D.Nouais, Phys.Lett. B373 (1996) 30.

  38. Double ratio definitions 2k* = p1 – p2 P = p1 + p2 simulation kside < 0 ksidesign selection arbitrary Correlation functions p1 kout > 0 p2 koutsign selection determined by the direction of the pair momentum P Double ratios kside > 0 kout > 0 p2 2k* [GeV/c] klong is the z component of the momentum of first particle in LCMS p1

  39. We are directly sensitive to time shift, the space shift arises from radial flow – possibility of a new radial flow measurement Out direction Flow velocity observed transverse velocity Side direction thermal velocity What to expect from double ratios • Initial separation in Pair Rest Frame (measured) can come from time shift and/or space shift in Source Frame (what we want to obtain)

  40. Separation between source 1 and 2 in pair rest frame What do we probe? Source of particle 1 Source of particle 2 sr* Boost to pair rest frame Dr Dt <Dr*> Dr* (fm) r (fm) Dr* = gpair (Dr– bpair Dt) Separation due to space and/or time shift • Mean shift (<Dr*>) seen in double ratio • Sigma (sr*) seen in height of CF

  41. Correlation functions and ratios Good agreement for like- sign and unlike-sign pairs points to similar emission process for K+ and K- CF Out Clear sign of emission asymmetry Side Two other ratios done as a double check – expected to be flat Long Preliminary

  42. Results for Pion-Proton 130 AGeV Λ peaks • Similar preliminary analysis done for pion-proton • We observe Lambda peaks at k*~minv of Λ • Good agreement for identical and non-identical charge combinations STAR preliminary

  43. K+ p K- anti-p Best Fit STAR preliminary Preliminary results for Kaon-Proton • Using data from Year2 (200 AGeV) – sufficient statistics • No corrections for momentum resolution done • No error estimation yet – fit indicates theoretical expectations

  44. What do we measure and how to compare it to the models? Is our fitting method working? And if yes, what does it tell us? Need to disentangle flow and time shift Modeling the emission asymmetry • Need models producing strong transverse radial flow: • Blast-wave as a baseline • RQMD • UrQMD • T. Humanic's rescattering model

  45. t Kt = pair Pt R Rside Rout Understanding modelsBlast wave = Flow baseline • Blast wave • Parameterizes source size (source radius) radial flow (average flow rapidity) and momentum distribution (temperature): • No time shift • Only spatial shift due to flow • Infinitely long cyllinder (neglects long contribution) Parameterization of the final state

  46. Blast wave: how does the flow work Average emission points Pion pt = 0.15 GeV/c t = 0.73 Kaon pt = 0.5 GeV/c t = 0.71 Proton pt = 1. GeV/c t = 0.73 Particle momentum Spatial shifts (Dr)

  47. Fitting and quantitative comparisons • Fits assume gaussian source in PRF • r*out distributions have non-gaussian tails • Use the same fitting procedure for models and data - correlation functions constructed with “Lednicky's weights” Example of r*out distribution from RQMD

  48. More points in βt needed to map and discriminate the flow profile – needs STAR upgrades in PID capability (TOF barrel) Comparing models to data • Rescattering models and blast-wave are consistent with data • Blast wave parameters constrained by STAR measurements • In models flow is required to reproduce the data

  49. “traditional” HBT axis STAR HBT Matrix (circa 2003) Analysis in progress published submitted Not shown: 3p Correlations (accepted PRL) asHBT Phase space density Correlations with Cascades dAu, pp Cascades

More Related