Open charm reconstruction in the ALICE experiment Elena Bruna Supervisor: Prof. Massimo Masera

Download Report

Transcript Open charm reconstruction in the ALICE experiment Elena Bruna Supervisor: Prof. Massimo Masera

Open charm reconstruction in
the ALICE experiment
Elena Bruna
Supervisor: Prof. Massimo Masera
Seminar for the end of 2nd year (XIX) – Torino, Dec 2nd 2005
Outline
• Physics motivations of open charm analysis in Heavy Ion
Collisions
• D+ → K-π+π+ : overview of the kinematics
• Measurement of open charm in the ALICE experiment
• Exclusive reconstruction of D+ → K-π+π+:
• Event generation and reconstruction
• Reconstruction of the secondary vertex
• Selection strategy
• Perspectives for the measurement of D+ elliptic flow
• Summary and work plans
Elena Bruna
2
Motivations for the Open Charm
physics in Heavy Ion Collisions
Elena Bruna
3
Heavy quarks as probes of
nuclear medium /1
• charm, bottom produced at early stages of the collision
(timescale ~ 1/mQ < QGP ~ 10 fm at LHC)
Studies of initial state effects: nuclear shadowing
Because of the very low x down to ~10-4 at LHC the so many gluons merge
together, affecting the partons densities at low x w.r.t. protons partons ones.
• thermal production
 The c quark might be produced in the plasma phase:
mc (~ 1.2 GeV) comparable with predicted Tplasma (~
0.6-0.8 GeV)
• open QQ production (not Drell-Yan) natural normalization
for QQ studies
 Quarkonia enhancement at low PT and suppression at
high PT.
Elena Bruna
4
Heavy quarks as probes of
nuclear medium /2
• charm, bottom have long lifetime (> QGP ) and can probe
the bulk, strongly interacting phase
Studies of final state effects:
1) radiative energy loss
Hard partons radiate gluons in the medium, lose energy and become
quenched. Heavy quarks are expected to lose less energy than light
quarks.
Nuclear modification factor
RD ( pT )
Rc.h. ( pT )
yield AA
RAA ( pT ) 
coll
yield pp  Nbin
High E  suppression of the produced particles (at high PT)  RAA≠1
It depends on the properties of the medium (gluon density,
temperature and volume), it provides information on such properties.
Elena Bruna
5
Heavy quarks as probes of
nuclear medium /3
Studies of final state effects:
2) anisotropic flow on the transverse plane
Elliptic Flow = collective motion of particles (due to high pressure arising from
compression and heating of nuclear matter) superimposed on top of the
thermal motion
Correlation between azimuthal angles  of outgoing particles
and the direction of the impact parameter (REACTION PLANE RP)
dX X 0
 1  2v1 cos(   RP )  2v2 cos 2  RP   .... 
d 2




v
c
o
s
2



2
R
P
Elliptic flow coefficient
High opacity of the medium (strongly interacting) 
high anisotropic flow  high v2
v2 provides information
onBruna
the opacity of the medium.
Elena
6
Few experimental results from RHIC /1
radiative energy loss - RAA of the D mesons
( PT spectra of e+e- from D semileptonic decays )
from QM05
from QM05
q = 0 GeV2/fm
dNg / dy = 1000
q = 4 GeV2/fm
q = 14 GeV2/fm
 Charm is suppressed! Suppression is
approximately the same as for hadrons.
 Challenge for energy loss models. Elena Bruna
Also pp and pA data are
needed as reference!7
Few experimental results from RHIC /2
anisotropic flow – v2 of the D mesons
(  spectra of e+e- from D semileptonic decays )
from QM05
from QM05
 Significant flow of charm quark as for light quarks
 Strong coupling of charm quark to the medium
 Indication for reduction of v2 at pT > 2 GeV/c
(PHENIX)
Elena Bruna
Also pp and pA data are
needed as reference!
8
D+ → K-++ : overview of the
kinematics
Elena Bruna
9
Why D+ → K-++ ?
Advantages…
1. D+ has a “long” mean life (~311mm compared to ~123
mm of the D0)
2. D+ → K-++ is a 3-charge body decay  the most
promising from an experimental point of view
3. D+ → K-++ has a relatively large branching ratio
(BR=9.2% compared to 3.8% for D0 → K-+).
…drawbacks
1. Combinatorial background for this 3-body channel is
larger than for D0 → K-+.
2. The average PT of the decay product is softer (~ 0.7
GeV/c compared to ~ 1 GeV/c)
Elena Bruna
10
Hadronic 3-charge-body decays of D+
D±
I(JP) = ½ (0-)
m = 1869.4 MeV/c2
c = 311.8 mm
(PDG ’04)
D+K-++
BR = 9.2 %
D+→K-++
Non Resonant BR = 8.8 %
D+→K*0(892)+→K-++
Resonant
BR = 1.3 %
D+→K*0(1430)+→K-++ Resonant
BR = 2.3 %
D+→K*0(1680)+→K-++ Resonant
BR = 3.8·10-3 %
Elena Bruna
11
Kinematics (1)
PT distributions of the generated
particles (ONLY PYTHIA generation,
NO propagation and reconstruction in
the detector)
K
Mean = 0.87 GeV/c
(nonresonant events)
D
Mean = 1.66 GeV/c

Mean = 0.67 GeV/c
Knowledge of the PT shapes of the
decay products important at the level of
Elena Bruna
the selection strategy
12
Kinematics (2)

Comparing with Pb-Pb central
events (ONLY HIJING generation,
NO propagation and reconstruction
in the detector):
Mean = 0.67 GeV/c
Mean = 0.50 GeV/c
PT distributions:
nonresonant D+ decay
K
HIJING central (normalized)
Mean = 0.87 GeV/c
Mean = 0.65 GeV/c
Elena Bruna
K and  from D+ are
harder than K and 
produced in a Pb-Pb event
13
Dalitz Plots: Kinematics (3)
Non resonant
Resonant
Sharp borders due to PYTHIA cut
off on the tails of distributions
Elena Bruna
14
Measurement of open charm in
the ALICE experiment
Elena Bruna
15
ALICE @ LHC
setup
Time Projection Chamber (TPC)
Tracking, PID (dE/dx)
-0.9<<0.9
HMPID
TRD
MUON SPECTR..
PHOS
Inner Tracking System
(ITS):
6 SILICON layers
(pixel, drift, strip)
Vertices reconstruction,
PID (dE/dx)
-0.9<<0.9
Time Of Flight (TOF)
Tracking, PID (time)
-0.9<<0.9 Elena Bruna
Size:
16 x 26 m
Weight: ~10,000 16
tons
Track Impact Parameter d0
SIGMA (fit)
expected d0 resolution (s)
d0 – d0 sim
MEAN (fit)
Elena Bruna
17
Track Impact Parameter : d0 pull
SIGMA (fit)
Calculate the pull
pull 
d 0  d 0 sim
err (d 0 )
MEAN (fit)
Elena Bruna
18
Exclusive reconstruction of
D + → K-  + +
Elena Bruna
19
Simulation strategy
Our purpose: exclusive reconstruction of D± in the ALICE barrel
(Inner Tracking System employed in the search for secondary vertexes)
Too large statistics
Central Pb-Pb
(108 events) would
~ 9 D+/Devent (b<3.5 fm,
be required to
in |y|<1
dN/dy = 6000,
study the signal!!
√s=5.5 TeV)
Signal and background events separately
generated with the Italian GRID
• 5’000 signal events with only D± decaying in K (using PYTHIA):
 Check the kinematics and the reconstruction
 Optimize the vertexing algorithm
• 20’000 background events (central Pb-Pb events using HIJING):
 cc pairs merged in addition in order to reproduce the charm yield
predicted by NLO pQCD calculations (≈ 118 per event)
 Tune the cuts (impact parameter cut,…) on the tracks to be
analyzed by the vertexing algorithm
 Evaluate the combinatorial background
Elena Bruna
20
Reconstructed signal events:
Dalitz Plots
(
From reconstructed tracks
: the info given by the generation are taken into account)
This is done as an internal cross-check procedure
Non resonant
Resonant
Elena Bruna
21
Reconstructed signal events:
D+ invariant mass
Mean
MEAN = 1.867 GeV/c2
RMS = 0.019 GeV/c2
this is not a complete
reconstruction of the
signal: tracks are grouped by means
of info. stored at generation time.
Knowledge of MINV resolution vs PT
is important when selecting the
signal candidates Elena Bruna
MINV Resolution (SIGMA
of the gaussian fit)
22
Reconstruction of the secondary
vertex for
D+ → K-++
 First idea: adapting and improving the method already written
for the primary vertex finding and fitting in p-p
 Second idea: writing a new secondary vertex finder and
comparing its performace with the previous ones
Elena Bruna
23
Vertex finder
•
Originally developed to find the primary vertex in p-p
 Based on the Straight Line Approximation of a track
(helix)
•
Main steps
1. The method receives N (N=3 in our case) tracks as input
2. Each track is approximated by a straight line in the
vicinity of the primary vertex
3. An estimation of the secondary vertex from each pair of
tracks is obtained evaluating the crossing point between
the 2 straight lines
4. The coordinates of secondary vertex are determined
averaging among all the track pairs:
x found 
1
N pairs
 xij
ij
y found 
1
N pairs
 yij
Elena Bruna
ij
z found 
1
N pairs
z
ij
ij
24
Improving the Straight Line
Vertex Finder
1.
Add a cut on the distance of closest approach (DCA)
between the two straight lines
 A pair of tracks is not used for the vertex estimation if
their distance of closest approach is > fDCAcut
2.
Use a weighted mean of the 2 DCA points
 In order to take into account the errors on the tracks
parameters
3.
Calculate a parameter representing the dispersion of the
vertices given by the track pairs (fSigma)
Elena Bruna
25
DCA cut effect
No DCAcut
X coord
RMS=179 μm
Finder- MC (mm)
Y coord
RMS=183 μm
Finder- MC (mm)
Z coord
RMS=166 μm
Finder- MC (mm)
fDCAcut = 1.5 mm
RMS=179 μm
Finder- MC (mm)
RMS=182 μm
Finder- MC (mm)
RMS=165 μm
Elena Bruna
FinderMC (mm)
fDCAcut = 0.7 mm
RMS=178 μm
Finder- MC (mm)
RMS=181 μm
Finder- MC (mm)
RMS=163 μm
Finder- MC (mm) 26
Weighted mean effect
Arithmetic mean
X coord
Weighted mean
RMS=179 μm
RMS=179 μm
Finder- MC (mm)
Y coord
Finder- MC (mm)
RMS=183 μm
RMS=183 μm
Finder- MC (mm)
Z coord
Finder- MC (mm)
RMS=166 μm
Finder- MC (mm)
Improved
resolution
on Z
RMS=160 μm
Elena Bruna
Finder- MC (mm)
27
Vertices dispersion
Dispersion fSigma = standard deviation of the 3 vertex
estimations obtained from each track pair
The DCA cut (at
0.7 mm) reduces
the dispersion
fSigma (cm)
Elena Bruna
28
All events
Cutting on fSigma
RMS=700 μm
Finder- MC (mm)
fSigma < 0.4 cm
• A cut fSigma < 0.4 cm cuts 0.5% of
the events and ≈30% of the overflows
and underflows (i.e. events for which
the VertexFinder misses the true
vertex by more than 1 mm)
RMS=224 μm
Finder- MC (mm)
• A cut fSigma < 0.07 cm (700 mm) cuts
6.4% of the events and gives a RMS
of 151 mm (for X coordinate)
fSigma < 0.07 cm
RMS=151 μm
Finder- MC (mm)
Elena Bruna
29
Another improvement:
Helix vertex finder
•
Based on the Distance of Closest Approach (DCA) between
helices

•
Does not use a Straight Line Approximation as the old one
Main steps
 The method receives N (N=3 in our case) tracks as input
 For each pair of tracks, the coordinates of the 2 points of closest
approach are calculated
 An estimation of the secondary vertex from each pair of tracks is
obtained averaging the coordinates of the points defining the DCA.
Two different implemetations: arithmetic vs. wieghted mean
 The coordinates of secondary vertex are determined averaging among
all the track pairs:
x found 
1
N pairs
 xij
ij
y found 
1
N pairs
 yij
ij
z found 
1
N pairs
z
ij
ij
 The dispersion of the vertices given by the track pairs is calculated
Elena Bruna
30
Results from the helix finder
Straight Line Finder
X coord
Helix Finder
RMS=169 μm
RMS=179 μm
Finder- MC (mm)
Finder- MC (mm)
Y coord
RMS=171 μm
RMS=183 μm
Finder- MC (mm)
Finder- MC (mm)
Z coord
RMS=162 μm
RMS=166 μm
Finder- MC (mm)
Helix finder
has better
resolution
and also a
lower number
of overflows
and
underflows
(≈400 instead
of ≈650)
Elena Bruna
Finder- MC (mm)
31
DCA cut effect on helix finders
fDCAcut=1 cm
X coord
RMS=169 μm
Finder- MC (mm)
Y coord
RMS=171 μm
Finder- MC (mm)
Z coord
RMS=162 μm
Finder- MC (mm)
fDCAcut=1.5 mm
RMS=168 μm
Finder- MC (mm)
RMS=170 μm
Finder- MC (mm)
RMS=161 μm
Elena MC
Bruna
Finder(mm)
fDCAcut=0.7 mm
RMS=167 μm
Finder- MC (mm)
RMS=169 μm
Finder- MC (mm)
RMS=158 μm
Finder- MC (mm)
32
Weighted mean effect on helix finder
Weighted mean
Arithmetic mean
X coord
RMS=168 μm
RMS=169 μm
Finder- MC (mm)
Y coord
Finder- MC (mm)
RMS=171 μm
RMS=169 μm
Finder- MC (mm)
Z coord
Finder- MC (mm)
RMS=162 μm
Finder- MC (mm)
Improved
resolution
on Z
RMS=154 μm
Elena Bruna
Finder- MC (mm)
33
Vertices dispersion on Helix Finder
The DCA cut
reduces the
dispersion
fSigma (cm)
Same distribution as for Straight Line finder
Elena Bruna
34
All events
Cutting on fSigma
RMS=480 μm
Finder- MC (mm)
fSigma < 0.4 cm
• A cut fSigma < 0.4 cm cuts 0.5% of
the events and ≈35% of the overflows
and underflows (i.e. events for which
the VertexFinder misses the true
vertex by more than 1 mm)
RMS=209 μm
Finder- MC (mm)
fSigma < 0.07 cm
• A cut fSigma < 0.07 cm (700 mm) cuts
5.6% of the events and gives a RMS
of 140 mm (for X coordinate)
RMS=140 μm
Finder- MC (mm)
Elena Bruna
35
New secondary vertex finder
Straight Line Approximation used → analytic method
Vertex coordinates (x0,y0,z0)
from minimization of:
D  d  d 2  d3
2
2
1
2
2
Where: d1,d2,d3 are the distances (weighted with the
errors on the tracks) of the vertex from the 3 tracks:
 x1  x0 

d1  
 sx 
2
2
2
 y1  y0   z1  z0 
 



 s

y

  sz 
σx = σy
2
P1 (x1,y1,z1)
d1
Elena Bruna
Secondary
Vertex
(x0,y0,z0)
36
Resolution of the vertex finder
RMS x
RMS y
At high Pt of D+ (Pt>5-6 GeV/c), the RMS in the bending plane increases,
instead of going down to ~15µm (spatial pixel resolution) as expected.
RMS z
Conclusion
New method improves RMS of
~40μm for PtD+ ~ 2GeV/c for x, y
and z with respect to previous
Helix vertex finder based on DCA
of Bruna
pairs of tracks.
Elena
37
Resolution at high Pt /1
Checks with events only made of pions show that the RMS on the bending
plane:
Decreases down to 50 µm if the 3 tracks have Pt ~ 2 GeV/c
Reaches a value of ~20 µm (in agreement with spatial pixel resolution) if
the 3 tracks have Pt =100 GeV/c
3 pion vertex: RMS in the bending plane vs. Pt
Elena Bruna
38
Resolution at high Pt /2
In the signal events, as the Pt of the D+ increases, the “daughters” become
more and more co-linear, resulting in a worse resolution along the D+
direction.
y
y’
π+
x’
π+
K-
bending plane
rotated
D+
Elena Bruna
x
39
Resolution in the rotated frame /1
Along the Pt of the D+ (x’ coord.)
Orthogonal to the Pt of the D+ (y’ coord.)
→ Along the Pt of the D+: as Pt increases (for Pt>5-6 GeV/c) the angles between
the decay tracks become smaller: in this coordinate the RMS increases
→ Orthogonal to the Pt of the D+: the RMS decreases as expected
Elena Bruna
40
Resolution in the rotated frame /2
Ratios:
RMS along Pt
RMS orthog Pt
RMS along Pt
RMS z
RMS orthog Pt
RMS z
Elena Bruna
41
Vertices dispersions/1
fSigma  d1  d 2  d3
2
2
2
Δx = XVertex FOUND – XVertex MC
Δx < 1000 μm
1000<Δx <3000 μm
3000<Δx <5000 μm
Δx > 5000 μm
fSigma bigger for
bad vertices
fSigma (cm)
Elena Bruna
42
Vertices dispersions/2
Vertices taken / Vertices Tot (“True” vertices)
“Fake” vertices (tracks coming from 3
different D+ vertices)
• fSigma < 0.7 cm cuts ~1% of the events
and gives a RMS of 130 μm
• fSigma < 0.5 cm cuts ~6% of the events
and gives a RMS of 110 μm
RMS x (μm)
(for X coordinate)
Mean x (μm)
Cut on fSigma
Elena Bruna
43
Conclusions on the finders
The Straight Line vertex finder:
♣ DCA cut: negligible effect on the RMS of the residual distributions, slightly
reduced number of overflows and underflows
♣ The use of a weighted mean: improves Z resolution by ≈6 mm
♣ Cutting on the dispersion fSigma: removes the events for which the
VertexFinder misses the true vertex by more than 1 mm and improves the
resolution
The Helix vertex finder:
♦ Has better resolution w.r.t. Straight Line finder (by approximately 10 mm)
♦ Has less overflows and underflows w.r.t. Straight Line finder
♦ DRAWBACK: the DCA between helices is obtained by minimization
♦ DCA cut, weighted mean and fSigma cut: improve the resolution
The Minimum Distance vertex finder:
♥ Has better resolution w.r.t. Helix finder (by approximately 30 mm)
♥ Has less overflows and underflows w.r.t. previous finders
♥ Is an analytic method
♥ Weighted mean and fSigma cut: improve the resolution
♥ Is presently THE candidate for first D+ analysis
A cut on fSigma has to be tuned (it can be done at analysis level)
Elena Bruna
44
D+ selection strategy
Elena Bruna
45
Tuning the cuts
GOAL: tune the cuts on both signal and background events and find the
cuts giving the best S/B. (S/B = 11% was found for the D0K-+)
CUT TIPOLOGIES:
1. On the single tracks used to “feed” the vertexer (Particle Identification,
pT, track impact parameter)
 reduce the number af all the possible combinations of track-triplets
in a central Pb-Pb collision (~ 1010 without any initial cut!!). It MUST be
cut by 4-5 orders of magnitude before using the more time-consuming
vertexer.
In progress.
2. Once the triplets are combined, additional cuts (invariant mass and
eventually pT, impact parameter) are mandatory before using the
vertexer. These cuts are done on the triplets.
To be done.
3. The third kind of cuts is applied on the quality of the secondary vertices
found (vertex dispersion-fSigma, pointing angle,…)
To be done.
Elena Bruna
46
Single track cuts /1
GOAL: find a compromise between the
number of background triplets and the
number of signals we want to take
HOW: for each triplet (both signal and bkg)
a loop on all the possible cuts (d0,Pt ,Pt K)
is done
%
SIGNAL
TAKEN
Cut on the track
impact parameter (d0)
Particle Id.
given by the
generation:
initial approach
Pt cut 
(GeV/c)
Pt cut K
(GeV/c)
MINIMUM
Triplet BKG
taken
d0 cut
(mm)
1-2
1,200
1,175
120
131
3-4
0,875
0,775
95
77.000
4-5
1,400
1,150
0
250.000
15 - 20
1,000
0,800
0
7.600.000
25 - 30
0,750
0,550
0
100.000.000
45 - 50
0,525
0,350
0
1.000.000.000
75 - 80
0,350
0,325
0
6.000.000.000
90 - 95
0,275
0,300
0
10.000.000.000
The number of BKG triplets is reduced by a
factor of ~100 when doing the cut on the
Invariant Mass within 3s (see slide 22)
Elena Bruna
47
Single track cuts /2
Triplet BKG
Bkg=Triplets
No cut on the track
impact parameter (d0)
Particle Id.
given by the
generation:
initial approach
Pt cut 
(GeV/c)
Pt cut K
(GeV/c)
d0 cut
(mm)
% MAX
signal taken
101 - 102
1,325
1,200
105
0,9
104 - 105
0,900
0,800
85
3,1
105 - 106
1,225
1,000
0
6,0
106 - 107
0,975
0,775
0
11,0
107 - 108
0,750
0,600
0
19,4
1010 – 1011
0,000
0,000
0
100,0
The number of BKG triplets is reduced by a factor of
~100 when doing the cut on the Invariant Mass
within 3s (see slide 22)
Cut on d0  lower
cuts on Pt (useful
up to Bkg ~105)
Elena Bruna
48
Tuning the single track cuts /2
When tuning a cut, one has to
keep in mind how the Pt
distribution of the D+ is modified
Pt reconstructed D+
Mean=2.5 GeV/c
Pt reconstructed D+
Pt cut (p) = 0.75 GeV/c
Pt cut (K) = 0.6 GeV/c
Mean=1.8 GeV/c
Ratio:
With cut / Wo cut
Elena Bruna
49
Perspectives for the
measurement of D+ elliptic flow
Elena Bruna
50
Measurement of v2
Elliptic Flow = correlation of particle emission angles with the reaction plane
(i.e. w.r.t the impact parameter of the collision)
 Calculate the 2nd order coefficient of
Fourier expansion of particle azimuthal
distribution relative to the reaction plane
v2  cos(2(  RP ))
 The reaction plane is unknown.
 Estimate the reaction plane from particle
azimuthal anisotropy:
 n = Event plane = estimator of the
unknown reaction plane
 w sin i 
1
1   i

n  tan


n
w
cos


i
i


 Calculate particle distribution relative to
the event plane
v'2  cos(2(  2 ))
 Correct for event plane resolution
v'2
v2 
cos22  RP 
 Resolution contains the unknown RP
 Can be extracted from sub-events
Event plane resolution
Elena Bruna
51
Motivation and Method

•
•

GOAL: Evaluate the statistical error bars for
measurements of v2 for D± mesons decaying in K
v2 vs. centrality (pT integrated)
v2 vs. pT in different centrality bins
TOOL: fast simulation
•
•
•
Assume to have only events with signal
Generate ND±(b, pT) events with 1 D± per event
For each event
1.
2.
3.
4.
5.
Generate a random reaction plane (fixed RP=0)
Get an event plane (with correct event plane resolution)
Generate the D+ azimuthal angle (φD) according to the probability
distribution p(φ)  1 + 2v2 cos [2(φ-RP)]
Smear φD with the experimental resolution on D± azimuthal angle
Calculate v′2(D+), event plane resolution and v2(D+)
Elena Bruna
52
D+ azimuthal angle resolution
MEAN
RMS
Average  resolution =
8 mrad = 0.47 degrees
Elena Bruna
53
D+ statistics
bmin-bmax
(fm)
s inel
Pb-Pb
(%)
Nevents
(106)
Ncc
D±
per
event
yield
per
event
0-3
3.6
0.72
118
45.8
3-6
11
2.2
82
31.8
6-9
18
3.6
42
16.3
9-12
25.4
5.1
12.5
4.85
12-18
42
8.4
1.2
0.47
Nevents for 2·107 Minimum Bias
triggers (without any requirement on
the impact parameter of the collision)
D+ selected after all the cuts is
still missing: for the time being
 e=1.5% (same as D0)
ND±(b, pT) selected = e × D+ reconstructed
 Total number of ND±(b, pT) selected
Normalized to 2·107 Minimum Bias Events
Elena Bruna
54
Results: v2 vs. centrality
2·107 Minimum
Bias events
bmin-bmax N(D±)selected sv2)
0-3
1070
0.024
3-6
2270
0.015
6-9
1900
0.016
9-12
800
0.026
12-18
125
0.09
• Error bars quite large


Would be larger in a scenario with worse event plane resolution
May prevent to draw conclusions in case of small anisotropy of D mesons
Elena Bruna
55
Results: v2 vs. pT
2·107 MB events
pT limits N(D±)sel
sv2)
pT limits N(D±)sel
sv2)
pT limits N(D±)sel
sv2)
0-0.5
140
0.06
0-0.5
120
0.06
0-0.5
50
0.10
0.5-1
280
0.04
0.5-1
230
0.05
0.5-1
100
0.07
1-1.5
390
0.04
1-1.5
330
0.04
1-1.5
140
0.06
1.5-2
360
0.04
1.5-2
300
0.04
1.5-2
125
0.06
2-3
535
0.03
2-3
450
0.03
2-3
190
0.05
3-4
250
0.05
3-4
210
0.05
3-4
90
0.07
4-8
265
0.05
4-8
220
0.05
4-8
95
0.07
8-15
50
0.11
8-15
40
0.11
8-15
20
0.15
Elena Bruna
56
Summary and work plans
• Preparatory checks on the kinematics and on the
reconstructed signal events: completed
• Secondary Vertex: completed
 the method of the Minimum Distance of 3 tracks is presently THE
candidate for first D+ analysis
 cuts on fSigma will be tuned at the analysis level
• D+ analysis cuts
 the work on the cuts on “single tracks” to feed the vertexer is in
progress: Pt, impact parameter, PID.
 The work on the cuts in the “triplets” and on the secondary vertices has
to be done.
• Analysis on D+ elliptic flow: in progress
Elena Bruna
57