1 / 47

Jaiby Joseph Ajish 11/9/2011

Charmed Meson measurements using a Silicon Tracker in Au+Au collisions at √ s NN = 200 GeV in STAR experiment at RHIC. Jaiby Joseph Ajish 11/9/2011. Outline. Introduction Why collide nuclei at high energies? RHIC and STAR Physics at RHIC Important observations Heavy quark sector

morrie
Télécharger la présentation

Jaiby Joseph Ajish 11/9/2011

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. Charmed Meson measurements using a Silicon Tracker in Au+Au collisions at √sNN = 200 GeV in STAR experiment at RHIC Jaiby Joseph Ajish 11/9/2011

  2. Outline • Introduction • Why collide nuclei at high energies? • RHIC and STAR • Physics at RHIC • Important observations • Heavy quark sector • Charm measurement using Silicon Tracker • Secondary vertexing • Proof of principle with Ks0 • Results • New results using TOF detector • Future

  3. Why collide nuclei at high energies? Phase Diagram Study the Strong Interaction at high temperatures/densities Understand how matter behaved at the dawn of the Universe (~10-6 s) Create and study the properties of the Quark-Gluon Plasma (QGP), a new phase of nuclear matter. (Net Baryon) Density Quantum Chromo Dynamics (QCD) is the theory of strong interactions. A phase transition is predicted at high temperatures and/or densities. Collisions of heavy ions at relativistic speeds creates extreme temperatures/densities: Nuclear MatterQuark-Gluon Plasma (deconfinedpartonic matter) Lattice QCD predicts the phase transition at: Tc ~ 150 -170 MeV and ρ ~ 1GeV/fm3

  4. PHOBOS BRAHMS RHIC PHENIX STAR AGS TANDEMS Relativistic Heavy Ion Collider (RHIC) 1 km exploringnuclear matteratextreme conditions over thelastdecade2000-2010

  5. RHIC Collisions Initial Conditions STAR Detector view of the event Initial high Q2 interactions Hadronization Freeze-Out Partonic matter QGP Experimental approach to induce the QCD phase transition: collide nuclei like Au+Au How to vary the T ? the Volume ? vary energy, Nr of participant Nucleons of Colliding Nuclei Collision systems used at RHIC are: Au+Au, Cu+Cu, d+Au and p+p at energies (7.7 GeV to 500 GeV for p+p)

  6. STAR Detector (in 2007) • The tracking system consisted of : • TPC : provides momentum, particle identification • Silicondetectors : • 1 layer of silicon strip detectors (SSD) and 3 layers of silicon drift detectors (SVT). • higher spatial resolution : pointing resolution of 250µm in transverse direction (at 1GeV) was achieved.

  7. MTD EMC barrel MRPC ToF barrel EMC End Cap FMS BBC FPD TPC FHC PMD Completed DAQ1000 HLT Ongoing HFT FGT STAR Detector (in 2010) Time of Flight –Full Barrel (Excellent Particle ID) Previous generation Silicon Detectors are removed and Heavy Flavor Tracker is being built for exclusive charm measurement

  8. schematic view of jet production hadrons Leading particle q q hadrons leading particle Physics @ RHICNew with Heavy Ions • Hard Parton Scattering • Jets and mini-jets (from hard-scattering of partons) 30 - 50 % of particle production • high pt leading particles • Extends into perturbative regime • Calculations reliable In p-p collisions hard parton scattering will lead to two jets emerging back-to-back with about equal energy

  9. schematic view of jet production hadrons leading particle • Scattered partonsthat propagate through hot and dense nuclear matter will radiate (lose) energy in colored medium • interaction of parton with partonicmatter will lead to: • suppression of angular correlations • suppression of high pTparticles • aka “parton energy loss” or “jet quenching” q q hadrons leading particle vacuum QGP Physics @ RHICNew with Heavy Ions

  10. Physics @ RHICImportant observations (light quarks [u,d,s]) • Suppression of angular correlations • In central Au+Au collisions the light hadrons in away-side jets are suppressed. • Not the casein p+p and d+Au • partons lose energy via gluon radiation Energy loss depends on properties of medium (gluon densities, size) and properties of “probe” (color charge, mass) Medium created at RHIC has very high opacity

  11. Physics @ RHICImportant observations (light quarks [u,d,s]) • suppression of high pT particles • In addition, a measurement of energy loss of high pTpartons using RAB shows significant suppression • Not the case in d+Au. NOT a cold nuclear matter effect. • partons lose energy via gluon radiation Nuclear modification factor RAB → energy loss in partonic mater RAB= (A-B pTspectra)/(p-p pTspectra * “volume”) Medium created at RHIC has very high opacity

  12. Physics @ RHIC –Heavy Quark Sector • Heavy flavor is mostly produced at the earlier stages of the collision via gluon fusion : • not affected by chiral symmetry restoration(i.e. mass is the same in-medium and vacuum) • production cross-section found to binary scale • ideal (calibrated) probe the medium createdin heavy ion collision • Look at the energy loss (RAA) [and elliptic flow (v2)] of heavy quarks. • Theoretical models predicted gluon radiative energy loss for heavy quarks to besmaller than of light quarks, which is not experimentally observed. 1) Non-photonic electrons (NPE) Method - decayed from charm and beauty hadrons 2) At pT ≥ 6 GeV/c, RAA(NPE) ~ RAA(h±) !!! 3) Surprising Results: contradicts pQCD predictions challenges our understanding of the energy loss mechanism Needs Direct measurement of D and B mesons

  13. Measurement via Semi-leptonic (indirect) channels • Indirect measurement through Semi-leptonic decay channels: • D0 e+ + X (BR : 6.9 %) • D+/-e+/- + X (BR : 17.2%) • ✔Large pT range. • Use of specific triggers • Relative contribution of electrons from B and D mesons are unknown. Measurement using azimuthal correlation of D mesons with e- Azimuthal correlations of e-h and e-D0 can be utilized to disentangle the charm and bottom contributions 13 ✔ Triggers on high pT electrons

  14. Measurement via hadronic (direct) channels • Direct measurement using a combinatorial method (combining K and π tracks) • Measurement of hadronic decay modes via invariant mass analysis. • D0 (D0bar)K-+(K+-) BR : 3.8 % • D+/-K BR : 9.2% • ✔ C and B contributions separated. • Limited to low momentum range. • No triggers, no decay vertex reconstruction • Challenging due to small decay length Cu+Cu 200 GeV TPC Only (Low pT)

  15. Measurement using Silicon Vertex Detector and decay vertex fit ✔ SVT/SSD not designed (thickness, geometry) for charm measurement ✔ Full reconstruction/fit of the decay vertex by combining K and π tracks – Some particle ID capabilities obtained from TPC dE/dx dE/dx bands of Kaon, Pion Challenges: Very short lived particles (average decay-length 70 μm) coupled to marginal SVT resolution (the pointing resolution is about 250 μm/GeV) Poor PID: Lack of TOF+SVT data sets, dE/dx has limited resolving power. Poor PID

  16. D0 decay length (Simulation) For <pT> ~ 1GeV/c, βγ ~ 0.54 average decay length ~ 65μm (in the transverse plane) un-boost in the collider!

  17. Distance ofClosest Approach resolution • run 7 Au+Au@200GeV (MinBias trigger). • DCA resolution as a function of inverse momentum. • Reflect the resolution and Multiple Coulomb Scattering. STAR preliminary • Including the silicon detectors in the tracking improves the pointing resolution. • with 4 silicon hits, the pointing resolution to the interaction point ~ 250 μmat P = 1GeV/c.

  18. D0 Decay Topology • We introduced the Full reconstruction/fit of the decay vertex • We used the full track error matrix (inside the beam pipe) for best error estimates • Cut Optimization is based on Monte Carlo studies • Better resolution in secondary vertex position is achieved with the fit method compared to usual helix swimming methods

  19. Secondary vertex fit(MC Data - pure D0 Events) Mean of the difference reconstructed -MC Rms of the difference reconstructed -MC Correlation between the reconstructed decay and MC • There is no systematic shift in reconstructed quantities. • The standard deviation of the distribution is flat at ~ 250 m , which is • of the order of the resolution of (SSD+SVT).

  20. Reconstructed Quantities (example) (MC Data - pure D0 Events) pT Resolution Invariant Mass Reconstructed Decay Vertex Resolution • Resolution : • Inv Mass sigma ~ 13 MeV (0.7%, after a gauss fit) • Decay Vertex Coordinates ~ 220 μm (transverse) • ~ 200 μm (z-direction) • The reconstructed parameters behave as • expected with the current detector resolution.

  21. Proof of principle with K0s • Test with K0s decay reconstruction (about 200x the decay distance of D0s): • c = 2.68 cm • Signed decay length : • an excess can be observed on the positive side of the decay length distribution, indicating the presence of long-lived decays. • use the decay length significanceL/L to improve the signal. • more appropriate because of the momentum dependence of the decay length. Before cut After cut Signal+background background After using a cutSL > 10, a clear peak at the K0S mass is observed.

  22. D0+D0bar Signal (in 2007 Data) • 24 Million Au+Au @ 200 GeV/c events are used for this analysis. • 3rd degree polynomial fit is used for background estimation. • Fit yielded an apparent signal significance of 10-σ (combined D0+D0bar signal) • Signal remained stable as cuts are varied. Pol3 + gaus Pol3 gaus fit Gaussian Mean = 1864.19 ± 10 MeV

  23. D0 and D0bar separately D0bar/D0 Ratio ~ 1.18 ± 0.24 Attempt to extract physics with polynomial fit method revealed some problems: A robust background estimation method needed to see if the peak observed was an artifact A Multi Variate analysis is in progress to measure signal using same sign background estimation method. The fact that about a third of the SVT/SSD system was dead during Run-7, combined with the marginal resolution of the previous generation silicon detector and combinatorial background is challenging to our efforts. ( compatible with unity indicating a vanishing μB)

  24. Ongoing Analysis with Multivariate Analysis (TMVA) • TMVA is a tool for simultaneous optimization of many correlated cuts. • Training samples for signal (pure D0) and background (`HIJING Au+Au’) are provided. It then produces a classifier output with weight files for signal and background. • After training, testing can be done with Data sample (MC Embedding/Real) MC D0 Embedding 2007 Au+Au Data (1-2% of available data) • Preliminary results looks promising but… work in progress.

  25. Recent charm measurements with Time Of Flight (TOF) Detector • STAR Time Of Flight (TOF) detector provides better particle ID (measures particle velocity β[β=v/c]) • dE/dx + TOF offers excellent K, π separation up to p ~ 1.5 - 2 GeV/c • New results use ~ 250 Million Au+Au Events from year 2010 and ~ 210 Million p+p events from 2009

  26. Recent charm measurements with Time Of Flight (TOF) Detector Corrected pT Spectra RAA RAA shows charm suppression at ~ 4 GeV/c Cross-section is found to binary scale, indicating its production via initial hard scattering at RHIC Charm Cross Section

  27. Future Heavy Flavor Tracker (HFT) • STAR is undergoing a detector upgrade for the unambiguous measurement of charm – The Heavy Flavor Tracker (HFT) • Low mass detector designed to identify mid-rapidity Charm and Beauty mesons and baryons through direct reconstruction, with unprecedented pointing resolution. • CMOS sensors will provide single track resolution ~ 20-30 µm at <pT> ~ 1GeV/c.

  28. Key Measurements of the HFT (2) Elliptic flow, v2 (1) Energy loss of direct D0 - Rcp (3) Charmed Baryon to Meson Enhancement • The methods we have developed • here are directly applicable in HFT

  29. Summary • For the 1st time in STAR, a secondary vertex reconstruction with full fit is developed using silicon vertex detectors. • Results obtained looks interesting, work in progress to wrap up the analysis with the most advanced tools in High Energy Physics – Multi Variate Analysis • Pioneering work that is directly applicable to the upcoming upgrade to STAR – Heavy Flavor Tracker (HFT)

  30. Thank you

  31. Back-Up

  32. Physics @ RHICImportant observations (Light flavors) Partonic Collectivity • Substantial elliptic flow (v2) signal observed for a variety of particle species. • Rapid Thermalization • v2 scaled by the number of valance quarks shows an apparent scaling • Development of anisotropy in the partonic stage of collision Behaves like an ideal fluid

  33. Charm Cross-Section Comparison at 200 GeV STAR and PHENIX do not agree about total charm production x-section NLO Ref: R. Vogt, arXiv:0709.2531v1 [hep-ph] Need precise, exclusive measurements

  34. Secondary Vertex Resolution Plots (x,y,z) Fit Method (central region) σXY ~ 55μm σZ ~ 25μm Z Y X Helix Swimming Method – using global parameters (central region) σXY ~ 150μm σZ ~ 135μm Helix Swimming Method – using DCAGeometry (central region) σXY ~ 140μm σZ ~ 125μm Simulation results shows that a factor of two was gained in secondary vertex Resolution

  35. Strategy of Reconstruction Cuts are applied in the analysis code to reduce background and to increase the candidate pool Select Event – Apply Event Level Cuts Select Trigger Cuts on Z-Vertex Position and its error Loop over Tracks – Apply Track Level Cuts Number of Silicon Hits Transverse DCA (DCAXY) Track Momentum etc. Pair Association - D0 Candidate Level Cuts rapidity, Cosine of Kaon decay angle etc. Decay Vertex Fit – Decay fit Level Cuts probability of fit, decay length error of decay length etc. Particle Identification – Apply PID Cuts |nσK|, |nσπ| Output Saved for offline Analysis

  36. Attempt to extract physics • Uncorrected pT spectra: • A normalized pT Spectra corrected for • acceptance and efficiency would be used to: • - extract total charm cross-section, freeze out parameters etc. • - calculate energy loss RAA The results with a polynomial background estimate seem to be inconsistent. A robust background estimation method needed to see if the peak observed was an artifact – a “same sign” background subtraction method was performed

  37. 2007 Production MinBias Cuts in 1st Production • cut changed • new cut Cuts in 2nd Production • EVENT level • triggerId : 200001, 200003, 200013 • Primary vertex position along the beam axis : • |zvertex| < 10 cm • Resolution of the primary vertex position along the beam axis: • |zvertex|< 200µm • TRACKS level • Number of hits in the vertex detectors: SiliconHits>1 • Transverse Momentum of tracks: • pT >.5GeV/c • Momentum of tracks • p >.8GeV/c • Ratio TPC hits Fitted/Possible > 0.51 • Pseudo-rapidity :||<1.2 • dEdxTrackLength>40 cm • DCA to Primary vertex (transverse), • DCAxy< .2 cm • Radius of first hit on track : • < 9 cm if number of silicon hits =2 • < 13 cm else • EVENT level • triggerId : 200001, 200003, 200013 • Primary vertex position along the beam axis : • |zvertex| < 10 cm • Resolution of the primary vertex position along the beam axis: • |zvertex|< 200µm • TRACKSlevel • Number of hits in the vertex detectors : • SiliconHits>2 (tracks with sufficient DCA resolution) • Transverse Momentum of tracks: • pT >.5GeV/c • Momentum of tracks: • p >.5GeV/c • Number of fitted: • TPC hits > 20 • Pseudo-rapidity :||<1(SSD acceptance) • dEdxTrackLength>40 cm • DCA to Primary vertex (transverse), • DCAxy< .1 cm

  38. Continued.. Cuts from Previous production Cuts in New Production • D0 candidate • |y(D0)|<1 • |cos(*)|<0.8 • DECAY FIT level • Probability of fit >0.01 && |sLength|<.1cm • Particle ID : ndEdx :|nK|<2.5, |nπ|<2.5 DECAY FIT level Probability of fit >0.1 && |sLength|<.1cm Particle ID : ndEdx :|nK|<2, |nπ|<2 In both productions we made a pico file for further analysis. Cuts Used for making a pico file Previous Production New Production • ndEdx :|nK|<2, |nπ|<2 • |cos(*)|<0.6 • DCA daughters < 300 µm |D0Eta|<1.85 |Cos(θ*)<0.6

  39. Uncorrected pT Spectra

  40. Physics @ RHIC Important observations (Light flavors) Partonic Energy Loss Partonic Collectivity • Substantial elliptic flow (v2) signal observed for a variety of particle species. • Rapid Thermalization • v2 scaled by the number of valance quarks shows an apparent scaling • Development of anisotropy in the partonic stage of collision • In central Au+Au collisions the light hadrons in away-side jets are suppressed. • Different for p+p and d+Au • In addition, a measurement of energy loss of high pTpartons using RAB shows significant suppression • partons lose energy via gluon radiation Medium created at RHIC has very high opacity Behaves like an ideal fluid

  41. ‘Same sign’ background subtraction

  42. Heavy Quark Energy Loss Puzzle – NPE Method Still the main method at RHIC STAR: Phys. Rew. Lett, 98, 192301(2007) and nucl-ex/0607012v3 1) Non-photonic electrons (NPE) decayed from - charm and beauty hadrons 2) At pT ≥ 6 GeV/c, RAA(NPE) ~ RAA(h±) !!! Contradicts naïve pQCD predictions Surprising results - - challenge our understanding of the energy loss mechanism - force us to RE-think about the elastic-collisions energy loss - Requiresdirect measurements of c- and b-hadrons.

  43. Measurement via Semi leptonic (indirect) channels • Indirect measurement through Semi-leptonic decay channels: • D0 e+ + X (BR : 6.9 %) • D+/-e+/- + X (BR : 17.2%) • ✔Large pT range. • Relative contribution of electrons from B and D mesons are unknown. • Use of specific triggers Measurement using azimuthal correlation of D mesons with e- Azimuthal correlation of open charm mesons with non-photonic Electroncan be utilized to disentangle the charm and bottom contributions[3] ✔ Triggers on high pT electrons Any information from direct reconstruction of D and B-mesons would help

  44. My contributions Charm Analysis • QA, Problem fixing, Resolution studies • Detailed studies for ‘online’/‘offline’ Cut optimization • Data Productions (Micro/Pico-DSTs) • First observation charmed meson signal in real data (from 2007 Au+Au dataset) • Signal extraction, optimization, fitting, pT binning • Embedding QA, Study of Systematics and Physics Analysis Service Work • acceptance of D-mesons with a prototype designfor the HFT upgrade

  45. An explanation for the non-consistent physics results Same cuts are used to produce this picture that were used in the polynomial fit case Red = Signal Blue = [(++) + (--)] • This was a big surprise!!! • The fact that about a third of the SVT/SSD system was dead during Run-7, combined with the marginal resolution of the previous generation silicon detector and combinatorial background limits our efforts. • A final effort to measure the signal using a multivariate analysis is in progress.

  46. Charm cross section vs √sNN

  47. Charm and beauty from e-D0 azimuthal correlations CHARM PYTHIA BEAUTY sign(e) ≠ sign(K) Like Sign Δϕ ~ π Like Sign Δϕ ~ 0 sign(e) = sign(K) Unlike Sign Δϕ ~ π A Mischke, Phys. Lett. B671, 361 (2009) Kent State University, 9 November 2010, S. Kabana

More Related