5d60c8f9 |
1 | // $Id$ |
2 | |
3 | //___________________________________________________________________________ |
4 | ///////////////////////////////////////////////////////////////////////////// |
5 | // // |
6 | // File reader for HLT tracks ESD // |
7 | // // |
8 | // loizides@ikf.uni-frankfurt.de // |
9 | ///////////////////////////////////////////////////////////////////////////// |
10 | |
11 | #include <Riostream.h> |
12 | #include <TMath.h> |
13 | #include <TString.h> |
14 | #include <TObjString.h> |
15 | #include <TTree.h> |
16 | #include <TFile.h> |
17 | #include <AliESD.h> |
18 | #include <AliESDtrack.h> |
19 | #include <AliESDHLTtrack.h> |
75a0c43e |
20 | #include <AliL3Track.h> |
46a140bc |
21 | #include <AliL3Vertex.h> |
04a02430 |
22 | #include <AliKalmanTrack.h> |
5d60c8f9 |
23 | #include <AliJetEventParticles.h> |
24 | #include "AliJetParticlesReaderHLT.h" |
25 | |
26 | ClassImp(AliJetParticlesReaderHLT) |
27 | |
28 | AliJetParticlesReaderHLT::AliJetParticlesReaderHLT(Bool_t bMapper, const Char_t* esdfilename) : |
d9a296b1 |
29 | AliJetParticlesReaderESD(0,esdfilename), |
04a02430 |
30 | fTrackerType(bMapper), |
31 | fMinHits(0), |
32 | fMinWeight(0) |
5d60c8f9 |
33 | { |
34 | //constructor |
35 | } |
36 | |
37 | /********************************************************************/ |
38 | |
39 | AliJetParticlesReaderHLT::AliJetParticlesReaderHLT( |
40 | Bool_t bMapper, |
41 | TObjArray* dirs, |
42 | const Char_t* esdfilename) : |
d9a296b1 |
43 | AliJetParticlesReaderESD(0,dirs,esdfilename), |
04a02430 |
44 | fTrackerType(bMapper), |
45 | fMinHits(0), |
46 | fMinWeight(0) |
47 | |
5d60c8f9 |
48 | { |
49 | //constructor |
50 | } |
51 | |
52 | /********************************************************************/ |
53 | |
54 | AliJetParticlesReaderHLT::~AliJetParticlesReaderHLT() |
55 | { |
56 | //desctructor |
57 | } |
58 | |
59 | Int_t AliJetParticlesReaderHLT::ReadESD(AliESD* esd) |
60 | { |
61 | //Reads one ESD |
62 | |
63 | if (esd == 0) |
64 | { |
65 | Error("ReadESD","ESD is NULL"); |
66 | return kFALSE; |
67 | } |
68 | |
04a02430 |
69 | Float_t mf = esd->GetMagneticField(); |
70 | if (mf <= 0.0) |
71 | { |
72 | Error("ReadESD","Magnetic Field is 0. Skipping to next event."); |
73 | return kFALSE; |
74 | } |
75 | AliKalmanTrack::SetMagneticField(mf/10.); |
76 | |
a86edc4e |
77 | Info("ReadESD","Reading Event %d",fCurrentDir*1000+fCurrentEvent); |
5d60c8f9 |
78 | if((!fOwner) || (fEventParticles==0)) |
79 | fEventParticles = new AliJetEventParticles(); |
80 | |
04a02430 |
81 | Int_t ntr=0; |
82 | if(fTrackerType){ |
83 | ntr =esd->GetNumberOfHLTHoughTracks(); |
84 | Info("ReadESD","Found %d conformal tracks.",ntr); |
85 | } else { |
86 | ntr=esd->GetNumberOfHLTHoughTracks(); |
87 | Info("ReadESD","Found %d hough tracks.",ntr); |
88 | } |
89 | fEventParticles->Reset(ntr); |
90 | |
91 | TString headdesc=""; |
92 | headdesc+="Run "; |
93 | headdesc+=esd->GetRunNumber(); |
94 | headdesc+=": Ev "; |
95 | headdesc+=esd->GetEventNumber(); |
96 | fEventParticles->SetHeader(headdesc); |
97 | |
5d60c8f9 |
98 | Double_t vertexpos[3];//vertex position |
99 | const AliESDVertex* kvertex = esd->GetVertex(); |
100 | if (kvertex == 0) |
101 | { |
102 | Info("ReadESD","ESD returned NULL pointer to vertex - assuming (0.0,0.0,0.0)"); |
103 | vertexpos[0] = 0.0; |
104 | vertexpos[1] = 0.0; |
105 | vertexpos[2] = 0.0; |
106 | } |
107 | else |
108 | { |
109 | kvertex->GetXYZ(vertexpos); |
110 | } |
111 | |
112 | fEventParticles->SetVertex(vertexpos[0],vertexpos[1],vertexpos[2]); |
74e5a625 |
113 | //cout << vertexpos[0] << " " << vertexpos[1] << " " << vertexpos[2] << endl; |
5d60c8f9 |
114 | |
75a0c43e |
115 | AliL3Track l3; |
46a140bc |
116 | AliL3Vertex v; |
117 | v.SetX(vertexpos[0]); |
c5368c62 |
118 | v.SetY(vertexpos[1]); |
119 | v.SetZ(vertexpos[2]); |
46a140bc |
120 | Double_t xc=0.,yc=0.,zc=0.; |
04a02430 |
121 | for (Int_t i = 0;i<ntr; i++) { |
122 | AliESDHLTtrack *kesdtrack; |
123 | if(fTrackerType){ |
124 | kesdtrack=esd->GetHLTConfMapTrack(i); |
125 | } else { |
126 | kesdtrack=esd->GetHLTHoughTrack(i); |
127 | } |
5d60c8f9 |
128 | |
04a02430 |
129 | if (kesdtrack == 0) |
5d60c8f9 |
130 | { |
131 | Error("ReadESD","Can not get track %d", i); |
132 | continue; |
133 | } |
134 | |
46a140bc |
135 | //const Float_t kpid=kesdtrack->GetPID(); |
136 | const Int_t knhits=kesdtrack->GetNHits(); |
137 | const Int_t kweight=kesdtrack->GetWeight(); |
138 | //cout << i << " " << kweight << " " << knhits << endl; |
139 | if((fMinHits>0) && (knhits<fMinHits)) continue; |
140 | if(kweight>1000) continue; //avoid ghosts |
141 | if((fMinWeight>0) && (kweight<fMinWeight)) continue; |
142 | |
143 | Float_t px=kesdtrack->GetPx(); |
144 | Float_t py=kesdtrack->GetPy(); |
145 | Float_t pz=kesdtrack->GetPz(); |
146 | |
75a0c43e |
147 | if(fTrackerType){ |
148 | //if(!kesdtrack->ComesFromMainVertex()) continue; |
74e5a625 |
149 | //cout << kesdtrack->GetPx() << " " << kesdtrack->GetPy() << " " << kesdtrack->GetPz() << endl; |
75a0c43e |
150 | l3.SetFirstPoint(kesdtrack->GetFirstPointX(),kesdtrack->GetFirstPointY(),kesdtrack->GetFirstPointZ()); |
151 | l3.SetLastPoint(kesdtrack->GetLastPointX(),kesdtrack->GetLastPointY(),kesdtrack->GetLastPointZ()); |
152 | l3.SetCharge(kesdtrack->GetCharge()); |
153 | l3.SetPt(kesdtrack->GetPt()); |
154 | l3.SetTgl(kesdtrack->GetTgl()); |
155 | l3.SetPsi(kesdtrack->GetPsi()); |
156 | l3.CalculateHelix(); |
46a140bc |
157 | l3.GetClosestPoint(&v,xc,yc,zc); |
158 | if(TMath::Abs(zc)>10.) continue; |
159 | l3.SetFirstPoint(vertexpos[0],vertexpos[1],vertexpos[2]); |
74e5a625 |
160 | //l3.CalculateHelix(); |
75a0c43e |
161 | l3.UpdateToFirstPoint(); |
46a140bc |
162 | px=l3.GetPx(); |
163 | py=l3.GetPy(); |
164 | pz=l3.GetPz(); |
74e5a625 |
165 | //cout << px << " " << py << " " << pz << endl; |
75a0c43e |
166 | } |
167 | |
46a140bc |
168 | const Float_t kpt=TMath::Sqrt(px*px+py*py); |
169 | const Float_t kp=TMath::Sqrt(pz*pz+kpt*kpt); |
170 | const Float_t keta=0.5*TMath::Log((kp+pz+1e-30)/(kp-pz+1e-30)); |
171 | const Float_t kphi=TMath::Pi()+TMath::ATan2(-py,-px); |
5d60c8f9 |
172 | |
04a02430 |
173 | if(IsAcceptedParticle(kpt,kphi,keta)) |
46a140bc |
174 | fEventParticles->AddParticle(px,py,pz,kp,i,kesdtrack->GetMCid(),knhits,kpt,kphi,keta); |
04a02430 |
175 | } //loop over tracks |
5d60c8f9 |
176 | |
5d60c8f9 |
177 | return kTRUE; |
178 | } |