Coding convention violations: suppression
[u/mrichter/AliRoot.git] / TOF / AliTOFtrackerV1.cxx
CommitLineData
d686d062 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14**************************************************************************/
15
16//--------------------------------------------------------------------//
17// //
3a646035 18// AliTOFtrackerV1 Class //
d686d062 19// Task: Perform association of the ESD tracks to TOF Clusters //
20// and Update ESD track with associated TOF Cluster parameters //
21// //
22// -- Authors : S. Arcelli, C. Zampolli (Bologna University and INFN) //
23// -- Contacts: Annalisa.De.Caro@cern.ch //
24// -- : Chiara.Zampolli@bo.infn.it //
25// -- : Silvia.Arcelli@bo.infn.it //
26// //
27//--------------------------------------------------------------------//
28
29#include <Rtypes.h>
30#include <TROOT.h>
31
32#include <TClonesArray.h>
33#include <TTree.h>
34#include <TFile.h>
35#include <TH1F.h>
36#include <TH2F.h>
5c7c93fa 37#include <TSeqCollection.h>
d686d062 38
d686d062 39#include "AliESDtrack.h"
40#include "AliESDEvent.h"
10d100d4 41#include "AliESDpid.h"
d686d062 42#include "AliLog.h"
43#include "AliTrackPointArray.h"
5c7c93fa 44#include "AliGeomManager.h"
45#include "AliCDBManager.h"
d686d062 46
d686d062 47#include "AliTOFRecoParam.h"
3a646035 48#include "AliTOFReconstructor.h"
d686d062 49#include "AliTOFcluster.h"
ca270b6d 50#include "AliTOFGeometry.h"
d686d062 51#include "AliTOFtrackerV1.h"
52#include "AliTOFtrack.h"
d0c2872c 53
d686d062 54extern TROOT *gROOT;
55
56ClassImp(AliTOFtrackerV1)
57
58//_____________________________________________________________________________
59AliTOFtrackerV1::AliTOFtrackerV1():
60 fRecoParam(0x0),
d686d062 61 fN(0),
62 fNseeds(0),
63 fNseedsTOF(0),
64 fngoodmatch(0),
65 fnbadmatch(0),
66 fnunmatch(0),
67 fnmatch(0),
b42a5ea1 68 fTracks(new TClonesArray("AliTOFtrack")),
69 fSeeds(new TClonesArray("AliESDtrack")),
d686d062 70 fHDigClusMap(0x0),
71 fHDigNClus(0x0),
72 fHDigClusTime(0x0),
73 fHDigClusToT(0x0),
74 fHRecNClus(0x0),
75 fHRecChi2(0x0),
76 fHRecDistZ(0x0),
77 fHRecSigYVsP(0x0),
78 fHRecSigZVsP(0x0),
79 fHRecSigYVsPWin(0x0),
80 fHRecSigZVsPWin(0x0)
81 {
82 //AliTOFtrackerV1 main Ctor
3a646035 83
d686d062 84 InitCheckHists();
85
86}
87//_____________________________________________________________________________
d686d062 88AliTOFtrackerV1::~AliTOFtrackerV1() {
89 //
90 // Dtor
91 //
92
93 SaveCheckHists();
94
5d456bcd 95 if(!(AliCDBManager::Instance()->GetCacheFlag())){
96 delete fRecoParam;
97 }
d686d062 98 delete fHDigClusMap;
99 delete fHDigNClus;
100 delete fHDigClusTime;
101 delete fHDigClusToT;
102 delete fHRecNClus;
103 delete fHRecChi2;
104 delete fHRecDistZ;
105 delete fHRecSigYVsP;
106 delete fHRecSigZVsP;
107 delete fHRecSigYVsPWin;
108 delete fHRecSigZVsPWin;
b42a5ea1 109 if (fTracks){
110 fTracks->Delete();
111 delete fTracks;
112 fTracks=0x0;
113 }
114 if (fSeeds){
115 fSeeds->Delete();
116 delete fSeeds;
117 fSeeds=0x0;
118 }
d686d062 119}
120//_____________________________________________________________________________
10d100d4 121void AliTOFtrackerV1::GetPidSettings(AliESDpid *esdPID) {
122 //
123 // Sets TOF resolution from RecoParams
124 //
125 if (fRecoParam)
126 esdPID->GetTOFResponse().SetTimeResolution(fRecoParam->GetTimeResolution());
127 else
128 AliWarning("fRecoParam not yet set; cannot set PID settings");
129}
130//_____________________________________________________________________________
d686d062 131Int_t AliTOFtrackerV1::PropagateBack(AliESDEvent* event) {
132 //
133 // Gets seeds from ESD event and Match with TOF Clusters
134 //
135
3a646035 136 // initialize RecoParam for current event
bafe290d 137 AliDebug(1,"Initializing params for TOF");
3a646035 138
139 fRecoParam = AliTOFReconstructor::GetRecoParam(); // instantiate reco param from STEER...
90b234fe 140
3a646035 141 if (fRecoParam == 0x0) {
142 AliFatal("No Reco Param found for TOF!!!");
143 }
144 //fRecoParam->Dump();
a825d829 145 //if(fRecoParam->GetApplyPbPbCuts())fRecoParam=fRecoParam->GetPbPbparam();
90b234fe 146 //fRecoParam->PrintParameters();
147
10d100d4 148 //Handle Time Zero information
149
150 Double_t timeZero=0.;
151 Double_t timeZeroMax=99999.;
d686d062 152
153 //Initialise some counters
154
155 fNseeds=0;
156 fNseedsTOF=0;
157 fngoodmatch=0;
158 fnbadmatch=0;
159 fnunmatch=0;
160 fnmatch=0;
161
162 Int_t ntrk=event->GetNumberOfTracks();
163 fNseeds = ntrk;
d686d062 164 TClonesArray &aESDTrack = *fSeeds;
165
166
167 //Load ESD tracks into a local Array of ESD Seeds
168
169 for (Int_t i=0; i<fNseeds; i++) {
170 AliESDtrack *t=event->GetTrack(i);
171 new(aESDTrack[i]) AliESDtrack(*t);
172 }
173
174 //Prepare ESD tracks candidates for TOF Matching
175 CollectESD();
176
177 //Matching Step
178 MatchTracks();
179
6ef25b64 180 AliInfo(Form("Number of matched tracks = %d (good = %d, bad = %d)",fnmatch,fngoodmatch,fnbadmatch));
d686d062 181
182 //Update the matched ESD tracks
183
184 for (Int_t i=0; i<ntrk; i++) {
185 AliESDtrack *t=event->GetTrack(i);
186 AliESDtrack *seed =(AliESDtrack*)fSeeds->UncheckedAt(i);
6819758a 187
af61c656 188 if ( (seed->GetStatus()&AliESDtrack::kTOFin)!=0 ) {
189 t->SetStatus(AliESDtrack::kTOFin);
190 //if(seed->GetTOFsignal()>0){
191 if ( (seed->GetStatus()&AliESDtrack::kTOFout)!=0 ) {
192 t->SetStatus(AliESDtrack::kTOFout);
193 t->SetTOFsignal(seed->GetTOFsignal());
194 t->SetTOFcluster(seed->GetTOFcluster());
195 t->SetTOFsignalToT(seed->GetTOFsignalToT());
196 t->SetTOFsignalRaw(seed->GetTOFsignalRaw());
197 t->SetTOFsignalDz(seed->GetTOFsignalDz());
46d7d82e 198 t->SetTOFsignalDx(seed->GetTOFsignalDx());
af61c656 199 t->SetTOFCalChannel(seed->GetTOFCalChannel());
200 Int_t tlab[3]; seed->GetTOFLabel(tlab);
201 t->SetTOFLabel(tlab);
24322e5c 202
e81ea7b6 203 Double_t alphaA = dynamic_cast<AliExternalTrackParam*>(t)->GetAlpha();
204 Double_t xA = dynamic_cast<AliExternalTrackParam*>(t)->GetX();
205 Double_t yA = dynamic_cast<AliExternalTrackParam*>(t)->GetY();
206 Double_t zA = dynamic_cast<AliExternalTrackParam*>(t)->GetZ();
207 Double_t p1A = dynamic_cast<AliExternalTrackParam*>(t)->GetSnp();
208 Double_t p2A = dynamic_cast<AliExternalTrackParam*>(t)->GetTgl();
209 Double_t p3A = dynamic_cast<AliExternalTrackParam*>(t)->GetSigned1Pt();
210 const Double_t *covA = dynamic_cast<AliExternalTrackParam*>(t)->GetCovariance();
211
212 // Make attention, please:
213 // AliESDtrack::fTOFInfo array does not be stored in the AliESDs.root file
214 // it is there only for a check during the reconstruction step.
24322e5c 215 Float_t info[10]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
216 seed->GetTOFInfo(info);
217 t->SetTOFInfo(info);
218 AliDebug(2,Form(" distance=%f; residual in the pad reference frame: dX=%f, dZ=%f", info[0],info[1],info[2]));
219
e81ea7b6 220 // Check done:
221 // by calling the AliESDtrack::UpdateTrackParams,
222 // the current track parameters are changed
223 // and it could cause refit problems.
224 // We need to update only the following track parameters:
225 // the track length and expected times.
226 // Removed AliESDtrack::UpdateTrackParams call
227 // Called AliESDtrack::SetIntegratedTimes(...) and
228 // AliESDtrack::SetIntegratedLength() routines.
229 /*
230 AliTOFtrack *track = new AliTOFtrack(*seed);
231 t->UpdateTrackParams(track,AliESDtrack::kTOFout); // to be checked - AdC
232 delete track;
233 Double_t time[10]; t->GetIntegratedTimes(time);
234 */
235
236 Double_t time[10]; seed->GetIntegratedTimes(time);
237 t->SetIntegratedTimes(time);
238
239 Double_t length = seed->GetIntegratedLength();
240 t->SetIntegratedLength(length);
241
242 Double_t alphaB = dynamic_cast<AliExternalTrackParam*>(t)->GetAlpha();
243 Double_t xB = dynamic_cast<AliExternalTrackParam*>(t)->GetX();
244 Double_t yB = dynamic_cast<AliExternalTrackParam*>(t)->GetY();
245 Double_t zB = dynamic_cast<AliExternalTrackParam*>(t)->GetZ();
246 Double_t p1B = dynamic_cast<AliExternalTrackParam*>(t)->GetSnp();
247 Double_t p2B = dynamic_cast<AliExternalTrackParam*>(t)->GetTgl();
248 Double_t p3B = dynamic_cast<AliExternalTrackParam*>(t)->GetSigned1Pt();
249 const Double_t *covB = dynamic_cast<AliExternalTrackParam*>(t)->GetCovariance();
250 AliDebug(2,"Track params -now(before)-:");
251 AliDebug(2,Form(" X: %f(%f), Y: %f(%f), Z: %f(%f) --- alpha: %f(%f)",
252 xB,xA,
253 yB,yA,
254 zB,zA,
255 alphaB,alphaA));
256 AliDebug(2,Form(" p1: %f(%f), p2: %f(%f), p3: %f(%f)",
257 p1B,p1A,
258 p2B,p2A,
259 p3B,p3A));
260 AliDebug(2,Form(" cov1: %f(%f), cov2: %f(%f), cov3: %f(%f)"
261 " cov4: %f(%f), cov5: %f(%f), cov6: %f(%f)"
262 " cov7: %f(%f), cov8: %f(%f), cov9: %f(%f)"
263 " cov10: %f(%f), cov11: %f(%f), cov12: %f(%f)"
264 " cov13: %f(%f), cov14: %f(%f), cov15: %f(%f)",
265 covB[0],covA[0],
266 covB[1],covA[1],
267 covB[2],covA[2],
268 covB[3],covA[3],
269 covB[4],covA[4],
270 covB[5],covA[5],
271 covB[6],covA[6],
272 covB[7],covA[7],
273 covB[8],covA[8],
274 covB[9],covA[9],
275 covB[10],covA[10],
276 covB[11],covA[11],
277 covB[12],covA[12],
278 covB[13],covA[13],
279 covB[14],covA[14]
280 ));
281 AliDebug(3,Form(" %6d %f %f %f %f % %6d %3d %f %f %f %f %f %f",
af61c656 282 i,
283 t->GetTOFsignalRaw(),
284 t->GetTOFsignal(),
285 t->GetTOFsignalToT(),
286 t->GetTOFsignalDz(),
46d7d82e 287 t->GetTOFsignalDx(),
af61c656 288 t->GetTOFCalChannel(),
289 t->GetTOFcluster(),
290 t->GetIntegratedLength(),
291 time[0], time[1], time[2], time[3], time[4]
292 )
293 );
294 }
d686d062 295 }
296 }
297
b42a5ea1 298 fSeeds->Clear();
299 fTracks->Clear();
d686d062 300 return 0;
301
302}
303//_________________________________________________________________________
304void AliTOFtrackerV1::CollectESD() {
305 //prepare the set of ESD tracks to be matched to clusters in TOF
306
307 Int_t seedsTOF1=0;
308 Int_t seedsTOF2=0;
309
d686d062 310 TClonesArray &aTOFTrack = *fTracks;
311 for (Int_t i=0; i<fNseeds; i++) {
312
313 AliESDtrack *t =(AliESDtrack*)fSeeds->UncheckedAt(i);
314 if ((t->GetStatus()&AliESDtrack::kTPCout)==0)continue;
315
d686d062 316 AliTOFtrack *track = new AliTOFtrack(*t); // New
af61c656 317 Float_t x = (Float_t)track->GetX(); //New
318
319 // TRD 'good' tracks, already propagated at 371 cm
320 if ( ( (t->GetStatus()&AliESDtrack::kTRDout)!=0 ) &&
321 ( x >= AliTOFGeometry::Rmin() ) ) {
322 if ( track->PropagateToInnerTOF() ) {
323
324 AliDebug(1,Form(" TRD propagated track till rho = %fcm."
325 " And then the track has been propagated till rho = %fcm.",
326 x, (Float_t)track->GetX()));
327
328 track->SetSeedIndex(i);
329 t->UpdateTrackParams(track,AliESDtrack::kTOFin);
330 new(aTOFTrack[fNseedsTOF]) AliTOFtrack(*track);
331 fNseedsTOF++;
332 seedsTOF1++;
333 }
d686d062 334 delete track;
335 }
336
af61c656 337 // Propagate the rest of TPCbp
d686d062 338 else {
af61c656 339 if ( track->PropagateToInnerTOF() ) {
340
341 AliDebug(1,Form(" TRD propagated track till rho = %fcm."
342 " And then the track has been propagated till rho = %fcm.",
343 x, (Float_t)track->GetX()));
344
d686d062 345 track->SetSeedIndex(i);
af61c656 346 t->UpdateTrackParams(track,AliESDtrack::kTOFin);
d686d062 347 new(aTOFTrack[fNseedsTOF]) AliTOFtrack(*track);
348 fNseedsTOF++;
349 seedsTOF2++;
350 }
351 delete track;
352 }
353 }
354
6ef25b64 355 AliInfo(Form("Number of TOF seeds = %d (Type 1 = %d, Type 2 = %d)",fNseedsTOF,seedsTOF1,seedsTOF2));
d686d062 356
357 // Sort according uncertainties on track position
358 fTracks->Sort();
359
360}
361//_________________________________________________________________________
362void AliTOFtrackerV1::MatchTracks( ){
363 //
364 //Match ESD tracks to clusters in TOF
365 //
366
367
368 // Parameters regulating the reconstruction
ca270b6d 369 Float_t dY=AliTOFGeometry::XPad();
370 Float_t dZ=AliTOFGeometry::ZPad();
d686d062 371
6819758a 372 const Float_t kTimeOffset = 32.; // time offset for tracking algorithm [ps]
373
3da9d64f 374 const Int_t kncmax = 100;
d686d062 375 Float_t sensRadius = fRecoParam->GetSensRadius();
376 Float_t scaleFact = fRecoParam->GetWindowScaleFact();
377 Float_t dyMax=fRecoParam->GetWindowSizeMaxY();
378 Float_t dzMax=fRecoParam->GetWindowSizeMaxZ();
379 Double_t maxChi2=fRecoParam->GetMaxChi2();
380 Bool_t timeWalkCorr = fRecoParam->GetTimeWalkCorr();
381 AliDebug(1,"++++++++++++++TOF Reconstruction Parameters:++++++++++++ \n");
382 AliDebug(1,Form("TOF sens radius: %f",sensRadius));
383 AliDebug(1,Form("TOF Window scale factor: %f",scaleFact));
384 AliDebug(1,Form("TOF Window max dy: %f",dyMax));
385 AliDebug(1,Form("TOF Window max dz: %f",dzMax));
386 AliDebug(1,Form("TOF Max Chi2: %f",maxChi2));
387 AliDebug(1,Form("Time Walk Correction? : %d",timeWalkCorr));
388
389
390 //The matching loop
d686d062 391 for (Int_t iseed=0; iseed<fNseedsTOF; iseed++) {
392
393 AliTOFtrack *track =(AliTOFtrack*)fTracks->UncheckedAt(iseed);
394 AliESDtrack *t =(AliESDtrack*)fSeeds->UncheckedAt(track->GetSeedIndex());
af61c656 395 //if ( t->GetTOFsignal()>0. ) continue;
396 if ( (t->GetStatus()&AliESDtrack::kTOFout)!=0 ) continue;
d686d062 397 AliTOFtrack *trackTOFin =new AliTOFtrack(*track);
398
399 // Determine a window around the track
400 Double_t x,par[5]; trackTOFin->GetExternalParameters(x,par);
401 Double_t cov[15]; trackTOFin->GetExternalCovariance(cov);
402
53884c34 403 if (cov[0]<0. || cov[2]<0.) {
404 AliWarning(Form("Very strange track (%d)! At least one of its covariance matrix diagonal elements is negative!",iseed));
ecd795d8 405 //delete trackTOFin;
406 //continue;
53884c34 407 }
408
d686d062 409 Double_t z = par[1];
ecd795d8 410 Double_t dz = scaleFact*3.*TMath::Sqrt(TMath::Abs(cov[2])+dZ*dZ/12.);
411 Double_t dphi = scaleFact*3.*TMath::Sqrt(TMath::Abs(cov[0])+dY*dY/12.)/sensRadius;
d686d062 412
413 Double_t phi=TMath::ATan2(par[0],x) + trackTOFin->GetAlpha();
414 if (phi<-TMath::Pi())phi+=2*TMath::Pi();
415 if (phi>=TMath::Pi())phi-=2*TMath::Pi();
416
417 //upper limit on window's size.
418 if(dz> dzMax) dz=dzMax;
419 if(dphi*sensRadius> dyMax) dphi=dyMax/sensRadius;
420
d686d062 421 // find the clusters inside the selected window
422 Int_t nc=0;
3da9d64f 423 AliTOFcluster *clusters[kncmax]; // pointers to the clusters in the window
424 Int_t index[kncmax];//to keep track of the cluster index
d686d062 425 for (Int_t k=FindClusterIndex(z-dz); k<fN; k++) {
426 AliTOFcluster *c=fClusters[k];
c0b60ae7 427 // if(nc>kncmax)break; /* R+ fix (buffer overflow) */
53884c34 428 if (nc>=kncmax) {
429 AliWarning("No more matchable clusters can be stored! Please, increase the corresponding vectors size.");
430 break; /* R+ fix (buffer overflow protection) */
431 }
d686d062 432 if(c->GetZ() > z+dz) break;
433 if(c->IsUsed()) continue;
17149e6b 434 if(!c->GetStatus()) {
53884c34 435 AliDebug(1,"Cluster in channel declared bad!");
436 continue; // skip bad channels as declared in OCDB
17149e6b 437 }
d686d062 438 Float_t xyz[3]; c->GetGlobalXYZ(xyz);
439 Double_t clPhi=TMath::ATan2(xyz[1],xyz[0]);
440 Double_t dph=TMath::Abs(clPhi-phi);
441 if (dph>TMath::Pi()) dph-=2.*TMath::Pi();
442 if (TMath::Abs(dph)>dphi) continue;
443 clusters[nc]=c;
444 index[nc] = k;
445 nc++;
446 }
447
53884c34 448 AliDebug(1,Form(" Number of matchable TOF clusters for the track number %d: %d",iseed,nc));
449
8dacd1bb 450 //start propagation: go to the average TOF pad middle plane at ~379.5 cm
d686d062 451
452 Float_t xTOF = sensRadius;
ca270b6d 453 Double_t ymax = xTOF*TMath::Tan(0.5*AliTOFGeometry::GetAlpha());
d686d062 454 Bool_t skip = kFALSE;
455 Double_t ysect = trackTOFin->GetYat(xTOF,skip);
456 if (skip) break;
457 if (ysect > ymax) {
ca270b6d 458 if (!trackTOFin->Rotate(AliTOFGeometry::GetAlpha())) {
d686d062 459 break;
460 }
461 } else if (ysect <-ymax) {
13c769b4 462 if (!trackTOFin->Rotate(-AliTOFGeometry::GetAlpha())) {
d686d062 463 break;
464 }
465 }
466 if(!trackTOFin->PropagateTo(xTOF)) {
467 break;
468 }
469
470
471 AliTOFcluster *bestCluster=0;
472 Double_t bestChi2=maxChi2;
473 Int_t idclus=-1;
c0b60ae7 474 // for (Int_t i=0; i<nc; i++){ /* R+ fix (unsafe) */
475 for (Int_t i=0; i<nc && i<kncmax; i++){ /* R+ fix (buffer overflow protection) */
d686d062 476 AliTOFcluster *c=clusters[i]; // one of the preselected clusters
477 Double_t chi2=trackTOFin->GetPredictedChi2((AliCluster3D*)c);
478 if (chi2 >= bestChi2) continue;
479 bestChi2=chi2;
480 bestCluster=c;
481 idclus=index[i];
482 }
483
484 if (!bestCluster) { // no matching , go to the next track
485 fnunmatch++;
486 delete trackTOFin;
487 continue;
488 }
489
490 fnmatch++;
32ead898 491
3a646035 492 AliDebug(2, Form("%7i %7i %10i %10i %10i %10i %7i",
32ead898 493 iseed,
3a646035 494 fnmatch-1,
32ead898 495 TMath::Abs(trackTOFin->GetLabel()),
496 bestCluster->GetLabel(0),
497 bestCluster->GetLabel(1),
3a646035 498 bestCluster->GetLabel(2),
499 idclus)); // AdC
32ead898 500
d686d062 501 bestCluster->Use();
502 if (
503 (bestCluster->GetLabel(0)==TMath::Abs(trackTOFin->GetLabel()))
504 ||
505 (bestCluster->GetLabel(1)==TMath::Abs(trackTOFin->GetLabel()))
506 ||
507 (bestCluster->GetLabel(2)==TMath::Abs(trackTOFin->GetLabel()))
508 ) {
509 fngoodmatch++;
6819758a 510 AliDebug(2,Form(" track label good %5d",trackTOFin->GetLabel()));
d686d062 511
512 }
513 else{
514 fnbadmatch++;
6819758a 515 AliDebug(2,Form(" track label bad %5d",trackTOFin->GetLabel()));
d686d062 516 }
517
518 //Propagate the track to the best matched cluster
519 trackTOFin->PropagateTo(bestCluster);
520
77a9ea9a 521 // Fill the track residual histograms.
522 FillResiduals(trackTOFin,bestCluster,kFALSE);
d686d062 523
524 //now take the local distance in Z from the pad center for time walk correction
ca270b6d 525 Float_t tiltangle = AliTOFGeometry::GetAngles(bestCluster->GetDetInd(1),bestCluster->GetDetInd(2))*TMath::DegToRad();
6819758a 526 Double_t dzTW=trackTOFin->GetZ()-bestCluster->GetZ(); // in cm - in the ALICE RF -
527 dzTW/=TMath::Cos(tiltangle); // from ALICE/tracking RF to pad RF (1)
528 dzTW=-dzTW; // from ALICE/tracking RF to pad RF (2)
529 if (tiltangle!=0.) AliDebug(2,Form(" rho_track = %f --- rho_cluster = %f ",trackTOFin->GetX(),bestCluster->GetX()));
d686d062 530
531 //update the ESD track and delete the TOFtrack
af61c656 532 t->UpdateTrackParams(trackTOFin,AliESDtrack::kTOFout);
d686d062 533
534 // Store quantities to be used in the TOF Calibration
ca270b6d 535 Float_t tToT=AliTOFGeometry::ToTBinWidth()*bestCluster->GetToT()*1E-3; // in ns
d686d062 536 t->SetTOFsignalToT(tToT);
6819758a 537 Float_t rawTime=AliTOFGeometry::TdcBinWidth()*bestCluster->GetTDCRAW()+kTimeOffset; // RAW time,in ps
d686d062 538 t->SetTOFsignalRaw(rawTime);
539 t->SetTOFsignalDz(dzTW);
6819758a 540
24322e5c 541 Float_t deltaY = trackTOFin->GetY()-bestCluster->GetY();
46d7d82e 542 t->SetTOFsignalDx(deltaY);
24322e5c 543
544 Float_t distR = (trackTOFin->GetX()-bestCluster->GetX())*
545 (trackTOFin->GetX()-bestCluster->GetX());
546 distR+=deltaY*deltaY;
547 distR+=dzTW*dzTW;
548 distR = TMath::Sqrt(distR);
549 Float_t info[10] = {distR, deltaY, dzTW,
550 0.,0.,0.,0.,0.,0.,0.};
551 t->SetTOFInfo(info);
552
d686d062 553 Int_t ind[5];
554 ind[0]=bestCluster->GetDetInd(0);
555 ind[1]=bestCluster->GetDetInd(1);
556 ind[2]=bestCluster->GetDetInd(2);
557 ind[3]=bestCluster->GetDetInd(3);
558 ind[4]=bestCluster->GetDetInd(4);
ca270b6d 559 Int_t calindex = AliTOFGeometry::GetIndex(ind);
d686d062 560 t->SetTOFCalChannel(calindex);
561
562 // keep track of the track labels in the matched cluster
563 Int_t tlab[3];
564 tlab[0]=bestCluster->GetLabel(0);
565 tlab[1]=bestCluster->GetLabel(1);
566 tlab[2]=bestCluster->GetLabel(2);
6819758a 567 AliDebug(2,Form(" tdc time of the matched track %6d = ",bestCluster->GetTDC()));
568 Double_t tof=AliTOFGeometry::TdcBinWidth()*bestCluster->GetTDC()+kTimeOffset; // in ps
569 AliDebug(2,Form(" tof time of the matched track: %f = ",tof));
d686d062 570 Double_t tofcorr=tof;
571 if(timeWalkCorr)tofcorr=CorrectTimeWalk(dzTW,tof);
572 AliDebug(2,Form(" tof time of the matched track, after TW corr: %f = ",tofcorr));
573 //Set TOF time signal and pointer to the matched cluster
574 t->SetTOFsignal(tofcorr);
575 t->SetTOFcluster(idclus); // pointing to the recPoints tree
576 t->SetTOFLabel(tlab);
577
6819758a 578 AliDebug(2,Form(" Setting TOF raw time: %f z distance: %f corrected time: %f",rawTime,dzTW,tofcorr));
579
d686d062 580 Double_t mom=t->GetP();
53884c34 581 AliDebug(2,Form(" Momentum for track %d -> %f", iseed,mom));
d686d062 582 // Fill Reco-QA histos for Reconstruction
583 fHRecNClus->Fill(nc);
584 fHRecChi2->Fill(bestChi2);
585 fHRecDistZ->Fill(dzTW);
586 fHRecSigYVsP->Fill(mom,TMath::Sqrt(cov[0]));
587 fHRecSigZVsP->Fill(mom,TMath::Sqrt(cov[2]));
588 fHRecSigYVsPWin->Fill(mom,dphi*sensRadius);
589 fHRecSigZVsPWin->Fill(mom,dz);
590
591 // Fill Tree for on-the-fly offline Calibration
592 // no longer there - all info is in the ESDs now
593
af61c656 594 delete trackTOFin;
d686d062 595 }
d686d062 596
597}
598//_________________________________________________________________________
599Int_t AliTOFtrackerV1::LoadClusters(TTree *cTree) {
600 //--------------------------------------------------------------------
601 //This function loads the TOF clusters
602 //--------------------------------------------------------------------
603
ca270b6d 604 Int_t npadX = AliTOFGeometry::NpadX();
605 Int_t npadZ = AliTOFGeometry::NpadZ();
606 Int_t nStripA = AliTOFGeometry::NStripA();
607 Int_t nStripB = AliTOFGeometry::NStripB();
608 Int_t nStripC = AliTOFGeometry::NStripC();
d686d062 609
610 TBranch *branch=cTree->GetBranch("TOF");
611 if (!branch) {
612 AliError("can't get the branch with the TOF clusters !");
613 return 1;
614 }
615
b42a5ea1 616 static TClonesArray dummy("AliTOFcluster",10000);
617 dummy.Clear();
618 TClonesArray *clusters=&dummy;
d686d062 619 branch->SetAddress(&clusters);
620
621 cTree->GetEvent(0);
622 Int_t nc=clusters->GetEntriesFast();
623 fHDigNClus->Fill(nc);
624
b42a5ea1 625 AliInfo(Form("Number of clusters: %d",nc));
626
d686d062 627 for (Int_t i=0; i<nc; i++) {
628 AliTOFcluster *c=(AliTOFcluster*)clusters->UncheckedAt(i);
16a2e36a 629//PH fClusters[i]=new AliTOFcluster(*c); fN++;
630 fClusters[i]=c; fN++;
d686d062 631
632 // Fill Digits QA histos
633
634 Int_t isector = c->GetDetInd(0);
635 Int_t iplate = c->GetDetInd(1);
636 Int_t istrip = c->GetDetInd(2);
637 Int_t ipadX = c->GetDetInd(4);
638 Int_t ipadZ = c->GetDetInd(3);
639
ca270b6d 640 Float_t time =(AliTOFGeometry::TdcBinWidth()*c->GetTDC())*1E-3; // in ns
641 Float_t tot = (AliTOFGeometry::TdcBinWidth()*c->GetToT())*1E-3;//in ns
d686d062 642
643 Int_t stripOffset = 0;
644 switch (iplate) {
645 case 0:
646 stripOffset = 0;
647 break;
648 case 1:
649 stripOffset = nStripC;
650 break;
651 case 2:
652 stripOffset = nStripC+nStripB;
653 break;
654 case 3:
655 stripOffset = nStripC+nStripB+nStripA;
656 break;
657 case 4:
658 stripOffset = nStripC+nStripB+nStripA+nStripB;
659 break;
660 default:
661 AliError(Form("Wrong plate number in TOF (%d) !",iplate));
662 break;
663 };
664 Int_t zindex=npadZ*(istrip+stripOffset)+(ipadZ+1);
665 Int_t phiindex=npadX*isector+ipadX+1;
666 fHDigClusMap->Fill(zindex,phiindex);
667 fHDigClusTime->Fill(time);
668 fHDigClusToT->Fill(tot);
669 }
670
671
672 return 0;
673}
674//_________________________________________________________________________
675void AliTOFtrackerV1::UnloadClusters() {
676 //--------------------------------------------------------------------
677 //This function unloads TOF clusters
678 //--------------------------------------------------------------------
679 for (Int_t i=0; i<fN; i++) {
16a2e36a 680//PH delete fClusters[i];
d686d062 681 fClusters[i] = 0x0;
682 }
683 fN=0;
684}
685
686//_________________________________________________________________________
687Int_t AliTOFtrackerV1::FindClusterIndex(Double_t z) const {
688 //--------------------------------------------------------------------
689 // This function returns the index of the nearest cluster
690 //--------------------------------------------------------------------
691 //MOD
692 //Here we need to get the Z in the tracking system
693
694 if (fN==0) return 0;
695 if (z <= fClusters[0]->GetZ()) return 0;
696 if (z > fClusters[fN-1]->GetZ()) return fN;
697 Int_t b=0, e=fN-1, m=(b+e)/2;
698 for (; b<e; m=(b+e)/2) {
699 if (z > fClusters[m]->GetZ()) b=m+1;
700 else e=m;
701 }
702 return m;
703}
704
705//_________________________________________________________________________
706Bool_t AliTOFtrackerV1::GetTrackPoint(Int_t index, AliTrackPoint& p) const
707{
708 // Get track space point with index i
709 // Coordinates are in the global system
710 AliTOFcluster *cl = fClusters[index];
711 Float_t xyz[3];
712 cl->GetGlobalXYZ(xyz);
713 Float_t phi=TMath::ATan2(xyz[1],xyz[0]);
714 Float_t phiangle = (Int_t(phi*TMath::RadToDeg()/20.)+0.5)*20.*TMath::DegToRad();
715 Float_t sinphi = TMath::Sin(phiangle), cosphi = TMath::Cos(phiangle);
ca270b6d 716 Float_t tiltangle = AliTOFGeometry::GetAngles(cl->GetDetInd(1),cl->GetDetInd(2))*TMath::DegToRad();
d686d062 717 Float_t sinth = TMath::Sin(tiltangle), costh = TMath::Cos(tiltangle);
ca270b6d 718 Float_t sigmay2 = AliTOFGeometry::XPad()*AliTOFGeometry::XPad()/12.;
719 Float_t sigmaz2 = AliTOFGeometry::ZPad()*AliTOFGeometry::ZPad()/12.;
d686d062 720 Float_t cov[6];
721 cov[0] = sinphi*sinphi*sigmay2 + cosphi*cosphi*sinth*sinth*sigmaz2;
722 cov[1] = -sinphi*cosphi*sigmay2 + sinphi*cosphi*sinth*sinth*sigmaz2;
723 cov[2] = -cosphi*sinth*costh*sigmaz2;
724 cov[3] = cosphi*cosphi*sigmay2 + sinphi*sinphi*sinth*sinth*sigmaz2;
725 cov[4] = -sinphi*sinth*costh*sigmaz2;
726 cov[5] = costh*costh*sigmaz2;
727 p.SetXYZ(xyz[0],xyz[1],xyz[2],cov);
728
729 // Detector numbering scheme
ca270b6d 730 Int_t nSector = AliTOFGeometry::NSectors();
731 Int_t nPlate = AliTOFGeometry::NPlates();
732 Int_t nStripA = AliTOFGeometry::NStripA();
733 Int_t nStripB = AliTOFGeometry::NStripB();
734 Int_t nStripC = AliTOFGeometry::NStripC();
d686d062 735
736 Int_t isector = cl->GetDetInd(0);
737 if (isector >= nSector)
738 AliError(Form("Wrong sector number in TOF (%d) !",isector));
739 Int_t iplate = cl->GetDetInd(1);
740 if (iplate >= nPlate)
741 AliError(Form("Wrong plate number in TOF (%d) !",iplate));
742 Int_t istrip = cl->GetDetInd(2);
743
744 Int_t stripOffset = 0;
745 switch (iplate) {
746 case 0:
747 stripOffset = 0;
748 break;
749 case 1:
750 stripOffset = nStripC;
751 break;
752 case 2:
753 stripOffset = nStripC+nStripB;
754 break;
755 case 3:
756 stripOffset = nStripC+nStripB+nStripA;
757 break;
758 case 4:
759 stripOffset = nStripC+nStripB+nStripA+nStripB;
760 break;
761 default:
762 AliError(Form("Wrong plate number in TOF (%d) !",iplate));
763 break;
764 };
765
766 Int_t idet = (2*(nStripC+nStripB)+nStripA)*isector +
767 stripOffset +
768 istrip;
769 UShort_t volid = AliGeomManager::LayerToVolUID(AliGeomManager::kTOF,idet);
770 p.SetVolumeID((UShort_t)volid);
771 return kTRUE;
772}
773//_________________________________________________________________________
774void AliTOFtrackerV1::InitCheckHists() {
775
776 //Init histos for Digits/Reco QA and Calibration
777
3a646035 778 TDirectory *dir = gDirectory;
779 TFile *logFileTOF = 0;
780
781 TSeqCollection *list = gROOT->GetListOfFiles();
782 int n = list->GetEntries();
783 Bool_t isThere=kFALSE;
784 for(int i=0; i<n; i++) {
785 logFileTOF = (TFile*)list->At(i);
786 if (strstr(logFileTOF->GetName(), "TOFQA.root")){
787 isThere=kTRUE;
788 break;
789 }
790 }
791
792 if(!isThere)logFileTOF = new TFile( "TOFQA.root","RECREATE");
793 logFileTOF->cd();
794
d686d062 795 //Digits "QA"
796 fHDigClusMap = new TH2F("TOFDig_ClusMap", "",182,0.5,182.5,864, 0.5,864.5);
797 fHDigNClus = new TH1F("TOFDig_NClus", "",200,0.5,200.5);
798 fHDigClusTime = new TH1F("TOFDig_ClusTime", "",2000,0.,200.);
799 fHDigClusToT = new TH1F("TOFDig_ClusToT", "",500,0.,100);
800
801 //Reco "QA"
802 fHRecNClus =new TH1F("TOFRec_NClusW", "",50,0.5,50.5);
803 fHRecDistZ=new TH1F("TOFRec_DistZ", "",50,0.5,10.5);
804 fHRecChi2=new TH1F("TOFRec_Chi2", "",100,0.,10.);
805 fHRecSigYVsP=new TH2F("TOFDig_SigYVsP", "",40,0.,4.,100, 0.,5.);
806 fHRecSigZVsP=new TH2F("TOFDig_SigZVsP", "",40,0.,4.,100, 0.,5.);
807 fHRecSigYVsPWin=new TH2F("TOFDig_SigYVsPWin", "",40,0.,4.,100, 0.,50.);
808 fHRecSigZVsPWin=new TH2F("TOFDig_SigZVsPWin", "",40,0.,4.,100, 0.,50.);
3a646035 809
810 dir->cd();
811
d686d062 812}
813
814//_________________________________________________________________________
815void AliTOFtrackerV1::SaveCheckHists() {
816
817 //write histos for Digits/Reco QA and Calibration
818
819 TDirectory *dir = gDirectory;
820 TFile *logFile = 0;
821 TFile *logFileTOF = 0;
822
823 TSeqCollection *list = gROOT->GetListOfFiles();
824 int n = list->GetEntries();
825 for(int i=0; i<n; i++) {
826 logFile = (TFile*)list->At(i);
827 if (strstr(logFile->GetName(), "AliESDs.root")) break;
828 }
829
830 Bool_t isThere=kFALSE;
831 for(int i=0; i<n; i++) {
832 logFileTOF = (TFile*)list->At(i);
833 if (strstr(logFileTOF->GetName(), "TOFQA.root")){
834 isThere=kTRUE;
835 break;
836 }
837 }
838
3a646035 839 if(!isThere) {
840 AliError(Form("File TOFQA.root not found!! not wring histograms...."));
841 return;
842 }
d686d062 843 logFile->cd();
844 fHDigClusMap->Write(fHDigClusMap->GetName(), TObject::kOverwrite);
845 fHDigNClus->Write(fHDigNClus->GetName(), TObject::kOverwrite);
846 fHDigClusTime->Write(fHDigClusTime->GetName(), TObject::kOverwrite);
847 fHDigClusToT->Write(fHDigClusToT->GetName(), TObject::kOverwrite);
848 fHRecNClus->Write(fHRecNClus->GetName(), TObject::kOverwrite);
849 fHRecChi2->Write(fHRecChi2->GetName(), TObject::kOverwrite);
850 fHRecDistZ->Write(fHRecDistZ->GetName(), TObject::kOverwrite);
851 fHRecSigYVsP->Write(fHRecSigYVsP->GetName(), TObject::kOverwrite);
852 fHRecSigZVsP->Write(fHRecSigZVsP->GetName(), TObject::kOverwrite);
853 fHRecSigYVsPWin->Write(fHRecSigYVsPWin->GetName(), TObject::kOverwrite);
854 fHRecSigZVsPWin->Write(fHRecSigZVsPWin->GetName(), TObject::kOverwrite);
855 logFile->Flush();
856
d686d062 857 dir->cd();
3a646035 858
d686d062 859 }
860//_________________________________________________________________________
861Float_t AliTOFtrackerV1::CorrectTimeWalk( Float_t dist, Float_t tof) {
862
863 //dummy, for the moment
864 Float_t tofcorr=0.;
ca270b6d 865 if(dist<AliTOFGeometry::ZPad()*0.5){
d686d062 866 tofcorr=tof;
867 //place here the actual correction
868 }else{
869 tofcorr=tof;
870 }
871 return tofcorr;
872}
873//_________________________________________________________________________
874Float_t AliTOFtrackerV1::GetTimeZerofromT0(AliESDEvent *event) const {
875
876 //Returns TimeZero as measured by T0 detector
877
878 return event->GetT0();
879}
880//_________________________________________________________________________
881Float_t AliTOFtrackerV1::GetTimeZerofromTOF(AliESDEvent * /*event*/) const {
882
883 //dummy, for the moment. T0 algorithm using tracks on TOF
884 {
885 //place T0 algo here...
886 }
887 return 0.;
888}
128563f6 889//_________________________________________________________________________
890
891void AliTOFtrackerV1::FillClusterArray(TObjArray* arr) const
892{
893 //
894 // Returns the TOF cluster array
895 //
896
897 if (fN==0)
898 arr = 0x0;
899 else
900 for (Int_t i=0; i<fN; ++i) arr->Add(fClusters[i]);
901
902}