]>
Commit | Line | Data |
---|---|---|
dadc7068 | 1 | // $Id$ |
2 | ||
3 | //************************************************************************** | |
4 | //* This file is property of and copyright by the ALICE HLT Project * | |
5 | //* ALICE Experiment at CERN, All rights reserved. * | |
6 | //* * | |
7 | //* Primary Authors: Gaute Ovrebekk <ovrebekk@ift.uib.no> * | |
8 | //* for The ALICE HLT Project. * | |
9 | //* * | |
10 | //* Permission to use, copy, modify and distribute this software and its * | |
11 | //* documentation strictly for non-commercial purposes is hereby granted * | |
12 | //* without fee, provided that the above copyright notice appears in all * | |
13 | //* copies and that both the copyright notice and this permission notice * | |
14 | //* appear in the supporting documentation. The authors make no claims * | |
15 | //* about the suitability of this software for any purpose. It is * | |
16 | //* provided "as is" without express or implied warranty. * | |
17 | //************************************************************************** | |
18 | ||
19 | /** @file AliHLTTPCTrackHistoComponent.cxx | |
20 | @author Gaute Ovrebekk, Matthias Richter | |
21 | @date | |
22 | @brief The TPC conformal mapping tracker component. | |
23 | */ | |
24 | ||
25 | ||
26 | #if __GNUC__>= 3 | |
27 | using namespace std; | |
28 | #endif | |
29 | ||
30 | #include "AliHLTTPCTrackHistoComponent.h" | |
31 | #include "AliHLTTPCTransform.h" | |
32 | #include "AliHLTTPCClusterDataFormat.h" | |
33 | #include "AliHLTTPCTrackletDataFormat.h" | |
34 | #include "AliHLTTPCMemHandler.h" | |
35 | #include "AliHLTTPCDefinitions.h" | |
36 | #include <TFile.h> | |
37 | #include <TString.h> | |
a267c68e | 38 | #include "TNtuple.h" |
dadc7068 | 39 | #include "TObjString.h" |
40 | #include "TObjArray.h" | |
a267c68e | 41 | #include "AliHLTTPCTrackArray.h" |
42 | #include "AliHLTTPCTrack.h" | |
dadc7068 | 43 | |
44 | //#include "AliHLTTPC.h" | |
45 | //#include <stdlib.h> | |
46 | //#include <cerrno> | |
47 | ||
dadc7068 | 48 | AliHLTTPCTrackHistoComponent gAliHLTTPCTrackHistoComponent; |
49 | ||
50 | /** ROOT macro for the implementation of ROOT specific class methods */ | |
51 | ClassImp(AliHLTTPCTrackHistoComponent) | |
52 | ||
53 | AliHLTTPCTrackHistoComponent::AliHLTTPCTrackHistoComponent() | |
54 | : | |
a267c68e | 55 | fClusters(NULL), |
56 | fTracks(NULL), | |
57 | fTracksArray(NULL) | |
dadc7068 | 58 | { |
59 | ||
60 | // see header file for class documentation | |
61 | // or | |
62 | // refer to README to build package | |
63 | // or | |
64 | // visit http://web.ift.uib.no/~kjeks/doc/alice-hlt | |
65 | ||
66 | } | |
67 | ||
68 | AliHLTTPCTrackHistoComponent::~AliHLTTPCTrackHistoComponent() | |
69 | { | |
70 | // see header file for class documentation | |
71 | } | |
72 | ||
73 | // Public functions to implement AliHLTComponent's interface. | |
74 | // These functions are required for the registration process | |
75 | ||
76 | const char* AliHLTTPCTrackHistoComponent::GetComponentID() | |
77 | { | |
78 | // see header file for class documentation | |
79 | ||
80 | return "TPCTrackHisto"; | |
81 | } | |
82 | ||
83 | void AliHLTTPCTrackHistoComponent::GetInputDataTypes(AliHLTComponentDataTypeList& list) | |
84 | { | |
85 | // see header file for class documentation | |
86 | list.clear(); | |
87 | list.push_back( AliHLTTPCDefinitions::fgkClustersDataType ); | |
88 | list.push_back( AliHLTTPCDefinitions::fgkTrackSegmentsDataType ); | |
89 | list.push_back( AliHLTTPCDefinitions::fgkTracksDataType ); | |
90 | } | |
91 | ||
92 | AliHLTComponentDataType AliHLTTPCTrackHistoComponent::GetOutputDataType() | |
93 | { | |
94 | // see header file for class documentation | |
a267c68e | 95 | return kAliHLTDataTypeTNtuple; |
dadc7068 | 96 | |
97 | } | |
98 | ||
99 | void AliHLTTPCTrackHistoComponent::GetOutputDataSize( unsigned long& constBase, double& inputMultiplier ) | |
100 | { | |
101 | // see header file for class documentation | |
102 | // XXX TODO: Find more realistic values. | |
103 | constBase = 0; | |
104 | inputMultiplier = 1; | |
105 | } | |
106 | ||
107 | AliHLTComponent* AliHLTTPCTrackHistoComponent::Spawn() | |
108 | { | |
109 | // see header file for class documentation | |
110 | return new AliHLTTPCTrackHistoComponent; | |
111 | } | |
112 | ||
113 | int AliHLTTPCTrackHistoComponent::DoInit( int argc, const char** argv ) | |
114 | { | |
a267c68e | 115 | |
116 | fClusters = new TNtuple("fCluster", "fCluster", "charge:qmax:residualY:residualZ:used:event"); | |
117 | fTracks = new TNtuple("fTracks", "fTracks", "pt:eta:psi:nclusters:event"); | |
118 | fTracksArray=new AliHLTTPCTrackArray(); | |
dadc7068 | 119 | |
120 | int iResult=0; | |
121 | TString configuration=""; | |
122 | TString argument=""; | |
123 | for (int i=0; i<argc && iResult>=0; i++) { | |
124 | argument=argv[i]; | |
125 | if (!configuration.IsNull()) configuration+=" "; | |
126 | configuration+=argument; | |
127 | } | |
128 | ||
129 | if (!configuration.IsNull()) { | |
130 | iResult=Configure(configuration.Data()); | |
131 | } | |
132 | return iResult; | |
133 | } | |
134 | ||
135 | int AliHLTTPCTrackHistoComponent::DoDeinit() | |
136 | { | |
137 | // see header file for class documentation | |
138 | ||
a267c68e | 139 | delete fClusters; |
140 | delete fTracks; | |
141 | delete fTracksArray; | |
142 | ||
dadc7068 | 143 | return 0; |
144 | } | |
145 | ||
146 | int AliHLTTPCTrackHistoComponent::DoEvent(const AliHLTComponentEventData& /*evtData*/, AliHLTComponentTriggerData& /*trigData*/) | |
147 | { | |
a267c68e | 148 | if(!fTracksArray){fTracksArray=new AliHLTTPCTrackArray();} |
149 | ||
dadc7068 | 150 | const AliHLTComponentBlockData* iter = NULL; |
151 | ||
152 | if ( GetFirstInputBlock( kAliHLTDataTypeSOR ) || GetFirstInputBlock( kAliHLTDataTypeEOR ) ) | |
153 | return 0; | |
154 | ||
155 | Int_t TotalTrack = 0; | |
156 | ||
157 | //Reading Merged Tracks | |
158 | for ( iter = GetFirstInputBlock(AliHLTTPCDefinitions::fgkTracksDataType); iter != NULL; iter = GetNextInputBlock() ) { | |
159 | if(iter->fDataType!=AliHLTTPCDefinitions::fgkTracksDataType){continue;} | |
160 | ReadTracks(iter,TotalTrack); | |
161 | } | |
162 | ||
163 | //Reading Tracks form slice | |
164 | for ( iter = GetFirstInputBlock(AliHLTTPCDefinitions::fgkTrackSegmentsDataType); iter != NULL; iter = GetNextInputBlock() ) { | |
165 | if(iter->fDataType!=AliHLTTPCDefinitions::fgkTrackSegmentsDataType){continue;} | |
166 | ReadTracks(iter,TotalTrack); | |
167 | } | |
168 | ||
169 | int TotalSpacePoint = 0; | |
170 | int nClustersUsed=0; | |
171 | ||
172 | for ( iter = GetFirstInputBlock(AliHLTTPCDefinitions::fgkClustersDataType); iter != NULL; iter = GetNextInputBlock() ) { | |
173 | ||
174 | if(iter->fDataType!=AliHLTTPCDefinitions::fgkClustersDataType){continue;} | |
175 | ||
a267c68e | 176 | AliHLTUInt8_t slice = AliHLTTPCDefinitions::GetMinSliceNr( *iter ); |
177 | AliHLTUInt8_t patch = AliHLTTPCDefinitions::GetMinPatchNr( *iter ); | |
dadc7068 | 178 | |
179 | //HLTDebug ( "Input Data - TPC cluster - Slice/Patch: %d/%d.", slice, patch ); | |
180 | const AliHLTTPCClusterData* clusterData = (const AliHLTTPCClusterData*) iter->fPtr; | |
181 | Int_t nSpacepoint = (Int_t) clusterData->fSpacePointCnt; | |
182 | TotalSpacePoint += nSpacepoint; | |
183 | //HLTInfo("TrackHisto found %d Spacepoints in slice %d patch %d", nSpacepoint, slice, patch); | |
184 | AliHLTTPCSpacePointData *clusters = (AliHLTTPCSpacePointData*) clusterData->fSpacePoints; | |
a267c68e | 185 | |
186 | if (fClustersArray[slice][patch]!=NULL) { | |
187 | delete(fClustersArray[slice][patch]); | |
188 | fClustersArray[slice][patch]=NULL; | |
189 | } | |
190 | Int_t arraysize=nSpacepoint*sizeof(AliHLTTPCSpacePointData); | |
191 | fClustersArray[slice][patch] = (AliHLTTPCSpacePointData*)new Byte_t[arraysize]; | |
192 | if (fClustersArray[slice][patch]) { | |
193 | memcpy(fClustersArray[slice][patch], clusters, arraysize); | |
194 | fNcl[slice][patch]=nSpacepoint; | |
195 | } else { | |
196 | fNcl[slice][patch]=nSpacepoint; | |
197 | HLTError ( "Memory allocation failed!" ); | |
198 | } | |
dadc7068 | 199 | |
200 | for(int i=0;i<nSpacepoint;i++){ | |
201 | UInt_t idCluster = clusters[i].fID; | |
202 | Int_t sliceCl = (idCluster>>25) & 0x7f; | |
203 | Int_t patchCl = (idCluster>>22) & 0x7; | |
204 | UInt_t pos = idCluster&0x3fffff; | |
a267c68e | 205 | Int_t used = 0; |
206 | Float_t resy = 0, resz = 0; | |
dadc7068 | 207 | for(UInt_t id=0;id<fTrackClusterID[sliceCl][patchCl].size();id++){ |
208 | if(fTrackClusterID[sliceCl][patchCl][id]==pos){ | |
209 | clusters[i].fUsed=kTRUE; | |
210 | nClustersUsed++; | |
a267c68e | 211 | used=1; |
212 | FillResidual(pos,sliceCl,patchCl,resy,resz); | |
dadc7068 | 213 | } |
a267c68e | 214 | } |
215 | if(used==1){ | |
216 | fClusters->Fill(clusters[i].fCharge,clusters[i].fQMax,resy,resz,used,GetEventId()); | |
217 | } | |
218 | else{ | |
219 | fClusters->Fill(clusters[i].fCharge,clusters[i].fQMax,-100,-100,used,GetEventId()); | |
220 | } | |
dadc7068 | 221 | } |
222 | } | |
a9f47f83 | 223 | |
dadc7068 | 224 | HLTInfo("TrackHisto found %d Spacepoints",TotalSpacePoint); |
225 | HLTInfo("TrackHisto found %d Tracks",TotalTrack); | |
226 | ||
227 | PushHisto(); | |
228 | ||
a267c68e | 229 | delete fTracksArray; |
230 | fTracksArray=NULL; | |
231 | ||
a9f47f83 | 232 | for(UInt_t i=0;i<36;i++){ |
233 | for(UInt_t j=0;j<6;j++){ | |
234 | fTrackClusterID[i][j].clear(); | |
235 | } | |
236 | } | |
dadc7068 | 237 | |
238 | return 0; | |
239 | } | |
240 | ||
241 | int AliHLTTPCTrackHistoComponent::Configure(const char* arguments) | |
242 | { | |
243 | ||
244 | int iResult=0; | |
245 | if (!arguments) return iResult; | |
246 | ||
247 | TString allArgs=arguments; | |
248 | TString argument; | |
249 | ||
250 | TObjArray* pTokens=allArgs.Tokenize(" "); | |
dadc7068 | 251 | if (pTokens) { |
252 | for (int i=0; i<pTokens->GetEntries() && iResult>=0; i++) { | |
253 | argument=((TObjString*)pTokens->At(i))->GetString(); | |
254 | if (argument.IsNull()) continue; | |
255 | ||
a267c68e | 256 | HLTError("unknown argument %s", argument.Data()); |
257 | iResult=-EINVAL; | |
258 | break; | |
dadc7068 | 259 | } |
260 | delete pTokens; | |
261 | } | |
262 | ||
263 | return iResult; | |
264 | } | |
265 | ||
266 | void AliHLTTPCTrackHistoComponent::ReadTracks(const AliHLTComponentBlockData* iter,Int_t &tt){ | |
267 | ||
a267c68e | 268 | AliHLTUInt8_t slice = AliHLTTPCDefinitions::GetMinSliceNr( *iter ); |
dadc7068 | 269 | //AliHLTUInt8_t patch = AliHLTTPCDefinitions::GetMinPatchNr( *iter ); |
270 | ||
271 | //HLTDebug ( "Input Data - TPC cluster - Slice/Patch: %d/%d.", slice, patch ); | |
a267c68e | 272 | AliHLTTPCTrackletData* trackData = (AliHLTTPCTrackletData*) iter->fPtr; |
dadc7068 | 273 | AliHLTUInt32_t nTracks = trackData->fTrackletCnt; |
a267c68e | 274 | fTracksArray->FillTracksChecked(trackData->fTracklets,trackData->fTrackletCnt,iter->fSize,slice,true); |
dadc7068 | 275 | tt += nTracks; |
276 | //HLTInfo("TrackHisto found %d Tracks in slice %d patch %d", nTracks, slice, patch); | |
277 | AliHLTTPCTrackSegmentData *tracks = (AliHLTTPCTrackSegmentData*) trackData->fTracklets; | |
278 | ||
279 | for(AliHLTUInt32_t i=0;i<nTracks;i++){ | |
dadc7068 | 280 | UInt_t nHits = tracks->fNPoints; |
a267c68e | 281 | |
282 | fTracks->Fill(tracks->fPt,tracks->fPsi,tracks->fTgl,nHits,GetEventId()); | |
283 | ||
dadc7068 | 284 | const UInt_t *hitnum = tracks->fPointIDs; |
285 | for(UInt_t h=0; h<nHits; h++){ | |
286 | UInt_t idTrack = hitnum[h]; | |
287 | Int_t sliceTrack = (idTrack>>25) & 0x7f; | |
288 | Int_t patchTrack = (idTrack>>22) & 0x7; | |
289 | UInt_t pos = idTrack&0x3fffff; | |
290 | fTrackClusterID[sliceTrack][patchTrack].push_back(pos); | |
291 | } | |
292 | UChar_t *tmpP = (UChar_t*)tracks; | |
293 | tmpP += sizeof(AliHLTTPCTrackSegmentData)+tracks->fNPoints*sizeof(UInt_t); | |
294 | tracks = (AliHLTTPCTrackSegmentData*)tmpP; | |
295 | } | |
296 | } | |
297 | ||
298 | void AliHLTTPCTrackHistoComponent::PushHisto(){ | |
299 | ||
dadc7068 | 300 | AliHLTUInt32_t fSpecification = AliHLTTPCDefinitions::EncodeDataSpecification(0,35,0,5); |
a267c68e | 301 | PushBack( (TObject*) fTracks,kAliHLTDataTypeTNtuple, fSpecification); |
302 | PushBack( (TObject*) fClusters,kAliHLTDataTypeTNtuple, fSpecification); | |
303 | ||
304 | } | |
305 | void AliHLTTPCTrackHistoComponent::FillResidual( UInt_t pos,AliHLTUInt8_t slice,AliHLTUInt8_t patch,Float_t& resy,Float_t& resz){ | |
306 | ||
307 | AliHLTTPCSpacePointData *cl = &fClustersArray[slice][patch][pos]; | |
308 | if(!cl){return;} | |
309 | ||
310 | AliHLTTPCTrack *gtrack = NULL; | |
311 | ||
af69d367 | 312 | for(int i=0;i<fTracksArray->GetNTracks();i++){ |
a267c68e | 313 | AliHLTTPCTrack *tt = fTracksArray->GetCheckedTrack(i); |
314 | UInt_t *hitnum =tt->GetHitNumbers(); | |
315 | Int_t nHits = tt->GetNHits(); | |
316 | for(Int_t h=0; h<nHits; h++){ | |
317 | UInt_t id=hitnum[h]; | |
318 | Int_t Tslice = (id>>25) & 0x7f; | |
319 | Int_t Tpatch = (id>>22) & 0x7; | |
320 | UInt_t Tpos = id&0x3fffff; | |
321 | if(Tslice==slice && Tpatch==patch && Tpos==pos) { | |
322 | gtrack = tt; | |
323 | break; | |
324 | } | |
dadc7068 | 325 | } |
dadc7068 | 326 | } |
a267c68e | 327 | |
328 | if(!gtrack){return;} | |
329 | ||
330 | Int_t tslice = gtrack->GetSector(); | |
331 | Double_t radius = gtrack->GetRadius(); // radius | |
332 | Double_t kappa = gtrack->GetKappa(); // curvature = 1/R , signed | |
333 | Double_t lambda = atan( gtrack->GetTgl() ); // dipAngle lambda | |
334 | ||
335 | // ------------------------------------ | |
336 | // ++ Get first/last point of the track | |
337 | ||
338 | Double_t xyzL[3]; // lastpoint of track | |
339 | Double_t xyzF[3]; // firstpoint of track | |
340 | ||
341 | xyzF[0] = gtrack->GetFirstPointX(); | |
342 | xyzF[1] = gtrack->GetFirstPointY(); | |
343 | xyzF[2] = gtrack->GetFirstPointZ(); | |
344 | ||
345 | xyzL[0] = gtrack->GetLastPointX(); | |
346 | xyzL[1] = gtrack->GetLastPointY(); | |
347 | xyzL[2] = gtrack->GetLastPointZ(); | |
348 | ||
349 | // -------------------------- | |
350 | // ++ Calculate length of the track | |
351 | ||
352 | Double_t s = 0.; // length of the track | |
353 | if ( AliHLTTPCTransform::GetBFieldValue() == 0. || kappa == 0 ) | |
354 | s = sqrt ( (xyzL[0] - xyzF[0])*(xyzL[0] - xyzF[0]) + (xyzL[1] - xyzF[1])*(xyzL[1] - xyzF[1]) ); | |
355 | else { | |
356 | // Calculate the length of the track. If it is to flat in in s,z plane use sxy, otherwise use sz | |
357 | if (fabs(lambda) > 0.05){ | |
358 | // length of track calculated out of z | |
359 | s = fabs( (xyzL[2] - xyzF[2]) / sin(lambda) ); // length of track calculated out of z | |
360 | } | |
361 | else { | |
362 | Double_t d = (xyzL[0] - xyzF[0])*(xyzL[0] - xyzF[0]) + (xyzL[1] - xyzF[1])*(xyzL[1] - xyzF[1]); | |
363 | // length of track calculated out of xy | |
364 | s = fabs ( acos( 0.5 * (2 - (d / (radius*radius)))) / ( kappa * cos(lambda) ) ); | |
365 | } | |
a9f47f83 | 366 | } |
dadc7068 | 367 | |
a267c68e | 368 | gtrack->Rotate(tslice,kTRUE); |
369 | ||
af69d367 | 370 | //Double_t padrows = 0; |
a267c68e | 371 | |
372 | Float_t xyzC[3]; // cluster tmp | |
373 | Float_t xyzTtmp[3]; // track tmp | |
374 | ||
375 | xyzC[0] = cl->fX; | |
376 | xyzC[1] = cl->fY; | |
377 | xyzC[2] = cl->fZ; | |
378 | ||
379 | Int_t padrow = AliHLTTPCTransform::GetPadRow(cl->fX); | |
380 | ||
381 | xyzTtmp[0] = gtrack->GetFirstPointX(); | |
382 | ||
383 | if(gtrack->GetCrossingPoint(padrow,xyzTtmp)) { | |
384 | // ---------------------- | |
385 | // ++ Calculate Residuals | |
386 | ||
387 | Float_t deltaY = ( xyzC[1] - xyzTtmp[1] ); | |
388 | Float_t deltaZ = ( xyzC[2] - xyzTtmp[2] ); | |
389 | ||
390 | resy = deltaY; | |
391 | resz = deltaZ; | |
392 | } | |
393 | else{ | |
394 | resy = -1000; | |
395 | resz = -1000; | |
5a90a5a4 | 396 | } |
dadc7068 | 397 | } |