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/2/2011
Feb 24, 2016
Charmed Meson measurements using a Silicon Tracker in Au+Au
collisions at √sNN = 200 GeV in STAR experiment at RHIC
Jaiby Joseph Ajish11/2/2011
2
OutlineQuick summary of my contributionsIntroduction
Why collide nuclei at high energies?RHIC and STAR
Physics at RHICImportant observationsHeavy quark sector
Charm measurement using Silicon TrackerSecondary vertexingProof of principle with Ks
0
ResultsNew results using TOF detectorFuture
3
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 design for the HFT upgrade
4
Why collide nuclei at high energies?
Quantum Chromo Dynamics (QCD) is the theory of strong interactions. A phase transition is predicted at sufficient high temperatures and/or densities.
Collisions of heavy ions at relativistic speeds creates extreme temperatures/densities: Nuclear Matter Quark-Gluon Plasma (deconfined
partonic matter)
Lattice QCD predicts the phase transition at:
Study the Strong Interaction at high temperatures/densities
Understand how matter behaved at the dawn of the Universe
Create and study the properties of the Quark-Gluon Plasma (QGP) phase of nuclear matter.
Phase Diagram
(Net Baryon) Density
Tc ~ 150 -170 MeV and ρ ~ 1GeV/fm3
5
RHIC BRAHMSPHOBOSPHENIX
STAR
AGS
TANDEMS
1 km
Relativistic Heavy Ion Collider (RHIC)
6
RHIC Collisions
Collision systems used at RHIC are: Au+Au, Cu+Cu, d+Au and p+p at different energies (7.7 GeV to 500 GeV for p+p)
Initial Conditions Initial high Q2
interactions Partonic matter QGP
HadronizationFreeze-Out
STAR Detector view of the event
7
STAR Detector (in 2007)
The tracking system consisted of : TPC : provides momentum, particle identification Silicon detectors :
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.
• Scattered partons that propagate through matter radiate energy (~ GeV/fm) in colored medium • interaction of parton with partonic matter• suppression of high pt particles
aka “parton energy loss” or “jet quenching”• suppression of angular correlation
vacu
umQG
P
• 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 hadrons
q
q
hadrons leadingparticle
leading particle
schematic view of jet production
Physics @ RHIC New with Heavy Ions
9
Physics @ RHIC Important observations (Light flavors)
Partonic Energy Loss
Medium created at RHIC has very high opacity
In central Au+Au collisions the light hadrons in away-side jets are suppressed.
Not the case in p+p and d+Au
In addition, a measurement of energy loss of high pT partons using RAB shows significant suppression
partons lose energy via gluon radiation
Nuclear modification factor RAA
→ energy loss in partonic materRAA = (A-A pT spectra)/(p-p pT spectra *
“volume”)
10
Physics @ RHIC – Heavy Quark Sector Heavy flavor is produced at the earlier stages of the collision via gluon fusion :
not affected by chiral symmetry restoration (i.e. mass is the same in/out of medium) production cross-section found to binary scale
ideal (calibrated) probe the medium created in heavy ion collision
Theoretical models predicted gluon radiative energy loss for heavy quarks to be smaller 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 predictionschallenges our understanding of the energy loss mechanism
Needs Direct measurement of D and B mesons
11
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
Poor PIDChallenges:1. Very short lived particles (average decay-length
70 μm) coupled to marginal SVT resolution (the pointing resolution is about 250 μm/GeV)
2. Poor PID: Lack of TOF+SVT data sets, dE/dx has limited resolving power.
dE/dx bands of Kaon, Pion
12
D0 Decay Topology
✿ We introduced the Full reconstruction/fit of the decay vertex
✿ We introduced and Used the full track error matrix for best error estimates
✿ Optimization of cuts based on MC studies
✿ Better resolution in secondary vertex position is achieved with the fit method compared to usual helix swimming methods.
Mean of the difference reconstructed -MC
Rms of the difference reconstructed -MC
Reco - MC [cm]
13
Reconstructed Quantities (example) (MC Data (pure D0) Events)
Invariant Mass Reconstructed pT Resolution
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.
14
After cut
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 significance L/L to improve the signal.– more appropriate because of the momentum dependence of the decay length.
Proof of principle with K0s
Before cut
After using a cut SL > 10, a clear peak at the K0S mass is observed.
backgroundSignal+background
15
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
Gaussian Mean = 1864.19 ± 10 MeV
gaus fit
16
D0 and D0bar separately
Attempt to extract physics with polynomial fit method revealed some problems:• Inconsistent pT spectra and x-section scaling depending on the production,
cuts etc.• A robust background estimation method needed to see if the peak
observed was an artifact• – a “same sign” background subtraction method was performed
17
• 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.
Same cuts are used to produce this picture that were used in the polynomial fit case
An explanation for the non-consistent physics results
Red = SignalBlue = [(++) + (--)]
18
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.
Recent charm measurements with Time Of Flight (TOF) Detector
• STAR Time Of Flight (TOF) detector provides better particle ID (measures particle velocity β)
• 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 p+p events from 2009
Corrected pT Spectrum and RAA in AuAuCharm Cross Section
Future Heavy Flavor Tracker (HFT)• STAR is undergoing a detector upgrade for the unambiguous measurement of
charm – The Heavy Flavor Tracker (HFT)
Key Measurements of HFT include:(1) Rcp (2) Elliptic flow, v2
(3) Charmed Baryon to Meson Enhancement
The methods we have developed are directly applicable in HFT
21
Thank you
22
Back-Up
Attempt to extract physics with 1st production
✿ 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
23
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
24
Distance of Closest Approach resolution
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 μm at P =
1GeV/c.
• run 7 Au+Au@200GeV (MinBias trigger).
• DCA resolution as a function of inverse momentum.
• Reflect the resolution and Multiple Coulomb Scattering.
25
Secondary Vertex fit – Simulation
✿ 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).
Reco vs. MC [cm] Mean of the difference reconstructed -MC
Rms of the difference reconstructed -MC
Reco - MC [cm]
26
Strategy of ReconstructionSelect Event – Apply Event Level Cuts
Select TriggerCuts on Z-Vertex Position and its error
Loop over Tracks – Apply Track Level CutsNumber of Silicon Hits
Transverse DCA (DCAXY)Track Momentum etc.
Pair Association - D0 Candidate Level Cutsrapidity, Cosine of Kaon decay angle etc.
Decay Vertex Fit – Decay fit Level Cutsprobability of fit, decay length
error of decay length etc.
Particle Identification – Apply PID Cuts|nσK|, |nσπ|
Cuts are applied in the analysis code to reduce background and to increase the candidate pool
Output Saved for offline Analysis
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.
27
Measurement using azimuthal correlation of D mesons with e-
Azimuthal correlation of open charm mesons with non-photonic Electron can be utilized to disentangle the charm and bottom contributions
28
Measurement via hadronic (direct) channelsDirect measurement using a combinatorial method– Measurement of hadronic decay modes via invariant mass analysis.
• D0 (D0)K-+(K+-) BR : 3.8 %• D+/-K BR : 9.2%
Results using STAR Time-Of-Flight (TOF) Detector(TOF+TPC offers better PID)TPC Only (Low pT)
✔ C and B contributions separated.✔ Limited to low momentum range.✔ No triggers, no decay vertex reconstruction✔ Challenging for charm mesons due to small decay length
Cuts in 1st Production
2007 Production MinBiasCuts in 2nd Production
EVENT leveltriggerId : 200001, 200003, 200013Primary vertex position along the beam axis : |zvertex| < 10 cmResolution of the primary vertex position along the beam axis: |zvertex|< 200µm
TRACKS levelNumber of hits in the vertex detectors :SiliconHits>2 (tracks with sufficient DCA resolution)Transverse Momentum of tracks: pT >.5GeV/cMomentum of tracks: p >.5GeV/c Number of fitted:TPC hits > 20 Pseudo-rapidity :||<1 (SSD acceptance)dEdxTrackLength>40 cmDCA to Primary vertex (transverse),DCAxy< .1 cm
EVENT leveltriggerId : 200001, 200003, 200013Primary vertex position along the beam axis : |zvertex| < 10 cmResolution of the primary vertex position along the beam axis: |zvertex|< 200µm
TRACKS levelNumber 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.51Pseudo-rapidity :||<1.2dEdxTrackLength>40 cmDCA to Primary vertex (transverse),DCAxy< .2 cmRadius of first hit on track :
< 9 cm if number of silicon hits =2 < 13 cm else
• cut changed• new cut
29
Cuts from Previous production
Continued..Cuts in New Production
DECAY FIT level
Probability of fit >0.1 && |sLength|<.1cm
Particle ID : ndEdx :|nK|<2, |nπ|<2
D0 candidate |y(D0)|<1|cos(*)|<0.8
DECAY FIT level
Probability of fit >0.01 && |sLength|<.1cm
Particle ID : ndEdx :|nK|<2.5, |nπ|<2.5
ndEdx :|nK|<2, |nπ|<2|cos(*)|<0.6DCA daughters < 300 µm
In both productions we made a pico file for further analysis.
Cuts Used for making a pico filePrevious Production New Production
|D0Eta|<1.85|Cos(θ*)<0.6
30
31
Uncorrected pT Spectra
32
Physics @ RHIC Important observations (Light flavors)
Partonic CollectivityPartonic Energy Loss
Behaves like an ideal fluidMedium created at RHIC has very high opacity
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 pT partons using RAB shows significant suppression
partons lose energy via gluon radiation
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
33
‘Same sign’ background subtraction
34
Heavy Quark Energy Loss Puzzle – NPE Method
Surprising results - - challenge our understanding of the energy loss mechanism - force us to RE-think about the elastic-collisions energy loss - Requires direct measurements of c- and b-hadrons.
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
STAR: Phys. Rew. Lett, 98, 192301(2007)
and nucl-ex/0607012v3
Still the main method at RHIC
Charm Cross-Section Comparison at 200 GeV
NLO Ref: R. Vogt, arXiv:0709.2531v1 [hep-ph]
STAR and PHENIX do not agree about total charm production x-section
Need precise, exclusive measurements
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 Electron can be utilized to disentangle the charm and bottom contributions[3]
36
✔ Triggers on high pT electrons
Any information from direct reconstruction of D and B-mesons would help