]> git.uio.no Git - u/mrichter/AliRoot.git/blame - TOF/AliTOFtrackerV1.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
A new method DrawPMDModule is added
[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// //
18// AliTOFtrackerV1 Class //
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"
41#include "AliLog.h"
42#include "AliTrackPointArray.h"
5c7c93fa 43#include "AliGeomManager.h"
44#include "AliCDBManager.h"
d686d062 45
46#include "AliTOFcalib.h"
47#include "AliTOFRecoParam.h"
48#include "AliTOFcluster.h"
ca270b6d 49#include "AliTOFGeometry.h"
d686d062 50#include "AliTOFtrackerV1.h"
51#include "AliTOFtrack.h"
5c7c93fa 52#include "AliTOFpidESD.h"
d0c2872c 53
d686d062 54extern TROOT *gROOT;
55
56ClassImp(AliTOFtrackerV1)
57
58//_____________________________________________________________________________
59AliTOFtrackerV1::AliTOFtrackerV1():
60 fRecoParam(0x0),
d686d062 61 fPid(0x0),
62 fN(0),
63 fNseeds(0),
64 fNseedsTOF(0),
65 fngoodmatch(0),
66 fnbadmatch(0),
67 fnunmatch(0),
68 fnmatch(0),
b42a5ea1 69 fTracks(new TClonesArray("AliTOFtrack")),
70 fSeeds(new TClonesArray("AliESDtrack")),
d686d062 71 fHDigClusMap(0x0),
72 fHDigNClus(0x0),
73 fHDigClusTime(0x0),
74 fHDigClusToT(0x0),
75 fHRecNClus(0x0),
76 fHRecChi2(0x0),
77 fHRecDistZ(0x0),
78 fHRecSigYVsP(0x0),
79 fHRecSigZVsP(0x0),
80 fHRecSigYVsPWin(0x0),
81 fHRecSigZVsPWin(0x0)
82 {
83 //AliTOFtrackerV1 main Ctor
84
d686d062 85 // Read the reconstruction parameters from the OCDB
ca270b6d 86 AliTOFcalib *calib = new AliTOFcalib();
d686d062 87 fRecoParam = (AliTOFRecoParam*)calib->ReadRecParFromCDB("TOF/Calib",-1);
88 if(fRecoParam->GetApplyPbPbCuts())fRecoParam=fRecoParam->GetPbPbparam();
d686d062 89 Double_t parPID[2];
90 parPID[0]=fRecoParam->GetTimeResolution();
91 parPID[1]=fRecoParam->GetTimeNSigma();
92 fPid=new AliTOFpidESD(parPID);
93 InitCheckHists();
ca270b6d 94 delete calib;
d686d062 95
96}
97//_____________________________________________________________________________
98AliTOFtrackerV1::AliTOFtrackerV1(const AliTOFtrackerV1 &t):
99 AliTracker(),
100 fRecoParam(0x0),
d686d062 101 fPid(0x0),
102 fN(0),
103 fNseeds(0),
104 fNseedsTOF(0),
105 fngoodmatch(0),
106 fnbadmatch(0),
107 fnunmatch(0),
108 fnmatch(0),
b42a5ea1 109 fTracks(new TClonesArray("AliTOFtrack")),
110 fSeeds(new TClonesArray("AliESDtrack")),
d686d062 111 fHDigClusMap(0x0),
112 fHDigNClus(0x0),
113 fHDigClusTime(0x0),
114 fHDigClusToT(0x0),
115 fHRecNClus(0x0),
116 fHRecChi2(0x0),
117 fHRecDistZ(0x0),
118 fHRecSigYVsP(0x0),
119 fHRecSigZVsP(0x0),
120 fHRecSigYVsPWin(0x0),
121 fHRecSigZVsPWin(0x0)
122 {
123 //AliTOFtrackerV1 copy Ctor
124
125 fNseeds=t.fNseeds;
126 fNseeds=t.fNseeds;
127 fNseedsTOF=t.fNseedsTOF;
128 fngoodmatch=t.fngoodmatch;
129 fnbadmatch=t.fnbadmatch;
130 fnunmatch=t.fnunmatch;
131 fnmatch=t.fnmatch;
132 fRecoParam=t.fRecoParam;
d686d062 133 fPid=t.fPid;
134 fSeeds=t.fSeeds;
135 fTracks=t.fTracks;
136 fN=t.fN;
137}
138
139//_____________________________________________________________________________
140AliTOFtrackerV1& AliTOFtrackerV1::operator=(const AliTOFtrackerV1 &t)
141{
142 //AliTOFtrackerV1 assignment operator
143
144 this->fNseeds=t.fNseeds;
145 this->fNseedsTOF=t.fNseedsTOF;
146 this->fngoodmatch=t.fngoodmatch;
147 this->fnbadmatch=t.fnbadmatch;
148 this->fnunmatch=t.fnunmatch;
149 this->fnmatch=t.fnmatch;
150 this->fRecoParam = t.fRecoParam;
d686d062 151 this->fPid = t.fPid;
152 this->fSeeds=t.fSeeds;
153 this->fTracks=t.fTracks;
154 this->fN=t.fN;
155 return *this;
156
157}
158//_____________________________________________________________________________
159AliTOFtrackerV1::~AliTOFtrackerV1() {
160 //
161 // Dtor
162 //
163
164 SaveCheckHists();
165
5d456bcd 166 if(!(AliCDBManager::Instance()->GetCacheFlag())){
167 delete fRecoParam;
168 }
d686d062 169 delete fPid;
170 delete fHDigClusMap;
171 delete fHDigNClus;
172 delete fHDigClusTime;
173 delete fHDigClusToT;
174 delete fHRecNClus;
175 delete fHRecChi2;
176 delete fHRecDistZ;
177 delete fHRecSigYVsP;
178 delete fHRecSigZVsP;
179 delete fHRecSigYVsPWin;
180 delete fHRecSigZVsPWin;
b42a5ea1 181 if (fTracks){
182 fTracks->Delete();
183 delete fTracks;
184 fTracks=0x0;
185 }
186 if (fSeeds){
187 fSeeds->Delete();
188 delete fSeeds;
189 fSeeds=0x0;
190 }
d686d062 191}
192//_____________________________________________________________________________
193Int_t AliTOFtrackerV1::PropagateBack(AliESDEvent* event) {
194 //
195 // Gets seeds from ESD event and Match with TOF Clusters
196 //
197
198
199 //Initialise some counters
200
201 fNseeds=0;
202 fNseedsTOF=0;
203 fngoodmatch=0;
204 fnbadmatch=0;
205 fnunmatch=0;
206 fnmatch=0;
207
208 Int_t ntrk=event->GetNumberOfTracks();
209 fNseeds = ntrk;
d686d062 210 TClonesArray &aESDTrack = *fSeeds;
211
212
213 //Load ESD tracks into a local Array of ESD Seeds
214
215 for (Int_t i=0; i<fNseeds; i++) {
216 AliESDtrack *t=event->GetTrack(i);
217 new(aESDTrack[i]) AliESDtrack(*t);
218 }
219
220 //Prepare ESD tracks candidates for TOF Matching
221 CollectESD();
222
223 //Matching Step
224 MatchTracks();
225
226 AliInfo(Form("Number of matched tracks: %d",fnmatch));
227 AliInfo(Form("Number of good matched tracks: %d",fngoodmatch));
228 AliInfo(Form("Number of bad matched tracks: %d",fnbadmatch));
229
230 //Update the matched ESD tracks
231
232 for (Int_t i=0; i<ntrk; i++) {
233 AliESDtrack *t=event->GetTrack(i);
234 AliESDtrack *seed =(AliESDtrack*)fSeeds->UncheckedAt(i);
235 if(seed->GetTOFsignal()>0){
236 t->SetTOFsignal(seed->GetTOFsignal());
237 t->SetTOFcluster(seed->GetTOFcluster());
238 t->SetTOFsignalToT(seed->GetTOFsignalToT());
239 t->SetTOFsignalRaw(seed->GetTOFsignalRaw());
240 t->SetTOFsignalDz(seed->GetTOFsignalDz());
241 t->SetTOFCalChannel(seed->GetTOFCalChannel());
242 Int_t tlab[3]; seed->GetTOFLabel(tlab);
243 t->SetTOFLabel(tlab);
244 AliTOFtrack *track = new AliTOFtrack(*seed);
245 t->UpdateTrackParams(track,AliESDtrack::kTOFout);
246 delete track;
247 }
248 }
249
250 //Handle Time Zero information
251
252 Double_t timeZero=0.;
253 Double_t timeZeroMax=99999.;
254 Bool_t usetimeZero = fRecoParam->UseTimeZero();
255 Bool_t timeZeroFromT0 = fRecoParam->GetTimeZerofromT0();
256 Bool_t timeZeroFromTOF = fRecoParam->GetTimeZerofromTOF();
257
258 AliDebug(1,Form("Use Time Zero?: %d",usetimeZero));
259 AliDebug(1,Form("Time Zero from T0? : %d",timeZeroFromT0));
260 AliDebug(1,Form("Time Zero From TOF? : %d",timeZeroFromTOF));
261
262 if(usetimeZero){
263 if(timeZeroFromT0){
264 timeZero=GetTimeZerofromT0(event);
265 }
266 if(timeZeroFromTOF && (timeZero>timeZeroMax || !timeZeroFromT0)){
267 timeZero=GetTimeZerofromTOF(event);
268 }
269 }
270 AliDebug(2,Form("time Zero used in PID: %f",timeZero));
271 //Make TOF PID
272 fPid->MakePID(event,timeZero);
273
b42a5ea1 274 fSeeds->Clear();
275 fTracks->Clear();
d686d062 276 return 0;
277
278}
279//_________________________________________________________________________
280void AliTOFtrackerV1::CollectESD() {
281 //prepare the set of ESD tracks to be matched to clusters in TOF
282
283 Int_t seedsTOF1=0;
284 Int_t seedsTOF2=0;
285
d686d062 286 TClonesArray &aTOFTrack = *fTracks;
287 for (Int_t i=0; i<fNseeds; i++) {
288
289 AliESDtrack *t =(AliESDtrack*)fSeeds->UncheckedAt(i);
290 if ((t->GetStatus()&AliESDtrack::kTPCout)==0)continue;
291
8dacd1bb 292 // TRD 'good' tracks, already propagated at 372 cm
d686d062 293
294 AliTOFtrack *track = new AliTOFtrack(*t); // New
295 Double_t x = track->GetX(); //New
296
297 if (((t->GetStatus()&AliESDtrack::kTRDout)!=0 ) &&
ca270b6d 298 ( x >= AliTOFGeometry::RinTOF()) ){
d686d062 299 track->SetSeedIndex(i);
300 t->UpdateTrackParams(track,AliESDtrack::kTOFout);
301 new(aTOFTrack[fNseedsTOF]) AliTOFtrack(*track);
302 fNseedsTOF++;
303 seedsTOF1++;
304 delete track;
305 }
306
307 // Propagate the rest of TPCbp
308
309 else {
310 if(track->PropagateToInnerTOF()){
311 track->SetSeedIndex(i);
312 t->UpdateTrackParams(track,AliESDtrack::kTOFout);
313 new(aTOFTrack[fNseedsTOF]) AliTOFtrack(*track);
314 fNseedsTOF++;
315 seedsTOF2++;
316 }
317 delete track;
318 }
319 }
320
321 AliInfo(Form("Number of TOF seeds %i",fNseedsTOF));
322 AliInfo(Form("Number of TOF seeds Type 1 %i",seedsTOF1));
323 AliInfo(Form("Number of TOF seeds Type 2 %i",seedsTOF2));
324
325 // Sort according uncertainties on track position
326 fTracks->Sort();
327
328}
329//_________________________________________________________________________
330void AliTOFtrackerV1::MatchTracks( ){
331 //
332 //Match ESD tracks to clusters in TOF
333 //
334
335
336 // Parameters regulating the reconstruction
ca270b6d 337 Float_t dY=AliTOFGeometry::XPad();
338 Float_t dZ=AliTOFGeometry::ZPad();
d686d062 339
3da9d64f 340 const Int_t kncmax = 100;
d686d062 341 Float_t sensRadius = fRecoParam->GetSensRadius();
342 Float_t scaleFact = fRecoParam->GetWindowScaleFact();
343 Float_t dyMax=fRecoParam->GetWindowSizeMaxY();
344 Float_t dzMax=fRecoParam->GetWindowSizeMaxZ();
345 Double_t maxChi2=fRecoParam->GetMaxChi2();
346 Bool_t timeWalkCorr = fRecoParam->GetTimeWalkCorr();
347 AliDebug(1,"++++++++++++++TOF Reconstruction Parameters:++++++++++++ \n");
348 AliDebug(1,Form("TOF sens radius: %f",sensRadius));
349 AliDebug(1,Form("TOF Window scale factor: %f",scaleFact));
350 AliDebug(1,Form("TOF Window max dy: %f",dyMax));
351 AliDebug(1,Form("TOF Window max dz: %f",dzMax));
352 AliDebug(1,Form("TOF Max Chi2: %f",maxChi2));
353 AliDebug(1,Form("Time Walk Correction? : %d",timeWalkCorr));
354
355
356 //The matching loop
357
d686d062 358 for (Int_t iseed=0; iseed<fNseedsTOF; iseed++) {
359
360 AliTOFtrack *track =(AliTOFtrack*)fTracks->UncheckedAt(iseed);
361 AliESDtrack *t =(AliESDtrack*)fSeeds->UncheckedAt(track->GetSeedIndex());
362 if(t->GetTOFsignal()>0. ) continue;
363 AliTOFtrack *trackTOFin =new AliTOFtrack(*track);
364
365 // Determine a window around the track
366 Double_t x,par[5]; trackTOFin->GetExternalParameters(x,par);
367 Double_t cov[15]; trackTOFin->GetExternalCovariance(cov);
368
369 Double_t z = par[1];
370 Double_t dz = scaleFact*3.*TMath::Sqrt(cov[2]+dZ*dZ/12);
371 Double_t dphi = scaleFact*3.*TMath::Sqrt(cov[0]+dY*dY/12.)/sensRadius;
372
373 Double_t phi=TMath::ATan2(par[0],x) + trackTOFin->GetAlpha();
374 if (phi<-TMath::Pi())phi+=2*TMath::Pi();
375 if (phi>=TMath::Pi())phi-=2*TMath::Pi();
376
377 //upper limit on window's size.
378 if(dz> dzMax) dz=dzMax;
379 if(dphi*sensRadius> dyMax) dphi=dyMax/sensRadius;
380
381
382 // find the clusters inside the selected window
383 Int_t nc=0;
3da9d64f 384 AliTOFcluster *clusters[kncmax]; // pointers to the clusters in the window
385 Int_t index[kncmax];//to keep track of the cluster index
d686d062 386 for (Int_t k=FindClusterIndex(z-dz); k<fN; k++) {
387 AliTOFcluster *c=fClusters[k];
3da9d64f 388 if(nc>kncmax)break;
d686d062 389 if(c->GetZ() > z+dz) break;
390 if(c->IsUsed()) continue;
391 if(!c->GetStatus()) continue; // skip bad channels as declared in OCDB
392 Float_t xyz[3]; c->GetGlobalXYZ(xyz);
393 Double_t clPhi=TMath::ATan2(xyz[1],xyz[0]);
394 Double_t dph=TMath::Abs(clPhi-phi);
395 if (dph>TMath::Pi()) dph-=2.*TMath::Pi();
396 if (TMath::Abs(dph)>dphi) continue;
397 clusters[nc]=c;
398 index[nc] = k;
399 nc++;
400 }
401
8dacd1bb 402 //start propagation: go to the average TOF pad middle plane at ~379.5 cm
d686d062 403
404 Float_t xTOF = sensRadius;
ca270b6d 405 Double_t ymax = xTOF*TMath::Tan(0.5*AliTOFGeometry::GetAlpha());
d686d062 406 Bool_t skip = kFALSE;
407 Double_t ysect = trackTOFin->GetYat(xTOF,skip);
408 if (skip) break;
409 if (ysect > ymax) {
ca270b6d 410 if (!trackTOFin->Rotate(AliTOFGeometry::GetAlpha())) {
d686d062 411 break;
412 }
413 } else if (ysect <-ymax) {
ca270b6d 414 if (!trackTOFin->Rotate(AliTOFGeometry::GetAlpha())) {
d686d062 415 break;
416 }
417 }
418 if(!trackTOFin->PropagateTo(xTOF)) {
419 break;
420 }
421
422
423 AliTOFcluster *bestCluster=0;
424 Double_t bestChi2=maxChi2;
425 Int_t idclus=-1;
426 for (Int_t i=0; i<nc; i++){
427 AliTOFcluster *c=clusters[i]; // one of the preselected clusters
428 Double_t chi2=trackTOFin->GetPredictedChi2((AliCluster3D*)c);
429 if (chi2 >= bestChi2) continue;
430 bestChi2=chi2;
431 bestCluster=c;
432 idclus=index[i];
433 }
434
435 if (!bestCluster) { // no matching , go to the next track
436 fnunmatch++;
437 delete trackTOFin;
438 continue;
439 }
440
441 fnmatch++;
442 bestCluster->Use();
443 if (
444 (bestCluster->GetLabel(0)==TMath::Abs(trackTOFin->GetLabel()))
445 ||
446 (bestCluster->GetLabel(1)==TMath::Abs(trackTOFin->GetLabel()))
447 ||
448 (bestCluster->GetLabel(2)==TMath::Abs(trackTOFin->GetLabel()))
449 ) {
450 fngoodmatch++;
451 AliDebug(2,Form(" track label good %5i",trackTOFin->GetLabel()));
452
453 }
454 else{
455 fnbadmatch++;
456 AliDebug(2,Form(" track label bad %5i",trackTOFin->GetLabel()));
457 }
458
459 //Propagate the track to the best matched cluster
460 trackTOFin->PropagateTo(bestCluster);
461
462
463 //now take the local distance in Z from the pad center for time walk correction
ca270b6d 464 Float_t tiltangle = AliTOFGeometry::GetAngles(bestCluster->GetDetInd(1),bestCluster->GetDetInd(2))*TMath::DegToRad();
d686d062 465 Double_t dzTW=trackTOFin->GetZ()-bestCluster->GetZ(); // in cm
466 dzTW/=TMath::Cos(tiltangle);
467
468 //update the ESD track and delete the TOFtrack
469 t->UpdateTrackParams(trackTOFin,AliESDtrack::kTOFout);
470 delete trackTOFin;
471
472 // Store quantities to be used in the TOF Calibration
ca270b6d 473 Float_t tToT=AliTOFGeometry::ToTBinWidth()*bestCluster->GetToT()*1E-3; // in ns
d686d062 474 t->SetTOFsignalToT(tToT);
ca270b6d 475 Float_t rawTime=AliTOFGeometry::TdcBinWidth()*bestCluster->GetTDCRAW()+32; // RAW time,in ps
d686d062 476 t->SetTOFsignalRaw(rawTime);
477 t->SetTOFsignalDz(dzTW);
478 AliDebug(2,Form(" Setting TOF raw time: %f z distance: %f time: %f = ",rawTime,dzTW));
479 Int_t ind[5];
480 ind[0]=bestCluster->GetDetInd(0);
481 ind[1]=bestCluster->GetDetInd(1);
482 ind[2]=bestCluster->GetDetInd(2);
483 ind[3]=bestCluster->GetDetInd(3);
484 ind[4]=bestCluster->GetDetInd(4);
ca270b6d 485 Int_t calindex = AliTOFGeometry::GetIndex(ind);
d686d062 486 t->SetTOFCalChannel(calindex);
487
488 // keep track of the track labels in the matched cluster
489 Int_t tlab[3];
490 tlab[0]=bestCluster->GetLabel(0);
491 tlab[1]=bestCluster->GetLabel(1);
492 tlab[2]=bestCluster->GetLabel(2);
493 AliDebug(2,Form(" tdc time of the matched track %i = ",bestCluster->GetTDC()));
ca270b6d 494 Double_t tof=AliTOFGeometry::TdcBinWidth()*bestCluster->GetTDC()+32; // in ps
d686d062 495 AliDebug(2,Form(" tof time of the matched track: %f = ",tof));
496 Double_t tofcorr=tof;
497 if(timeWalkCorr)tofcorr=CorrectTimeWalk(dzTW,tof);
498 AliDebug(2,Form(" tof time of the matched track, after TW corr: %f = ",tofcorr));
499 //Set TOF time signal and pointer to the matched cluster
500 t->SetTOFsignal(tofcorr);
501 t->SetTOFcluster(idclus); // pointing to the recPoints tree
502 t->SetTOFLabel(tlab);
503
504 Double_t mom=t->GetP();
505 // Fill Reco-QA histos for Reconstruction
506 fHRecNClus->Fill(nc);
507 fHRecChi2->Fill(bestChi2);
508 fHRecDistZ->Fill(dzTW);
509 fHRecSigYVsP->Fill(mom,TMath::Sqrt(cov[0]));
510 fHRecSigZVsP->Fill(mom,TMath::Sqrt(cov[2]));
511 fHRecSigYVsPWin->Fill(mom,dphi*sensRadius);
512 fHRecSigZVsPWin->Fill(mom,dz);
513
514 // Fill Tree for on-the-fly offline Calibration
515 // no longer there - all info is in the ESDs now
516
517 }
d686d062 518
519}
520//_________________________________________________________________________
521Int_t AliTOFtrackerV1::LoadClusters(TTree *cTree) {
522 //--------------------------------------------------------------------
523 //This function loads the TOF clusters
524 //--------------------------------------------------------------------
525
ca270b6d 526 Int_t npadX = AliTOFGeometry::NpadX();
527 Int_t npadZ = AliTOFGeometry::NpadZ();
528 Int_t nStripA = AliTOFGeometry::NStripA();
529 Int_t nStripB = AliTOFGeometry::NStripB();
530 Int_t nStripC = AliTOFGeometry::NStripC();
d686d062 531
532 TBranch *branch=cTree->GetBranch("TOF");
533 if (!branch) {
534 AliError("can't get the branch with the TOF clusters !");
535 return 1;
536 }
537
b42a5ea1 538 static TClonesArray dummy("AliTOFcluster",10000);
539 dummy.Clear();
540 TClonesArray *clusters=&dummy;
d686d062 541 branch->SetAddress(&clusters);
542
543 cTree->GetEvent(0);
544 Int_t nc=clusters->GetEntriesFast();
545 fHDigNClus->Fill(nc);
546
b42a5ea1 547 AliInfo(Form("Number of clusters: %d",nc));
548
d686d062 549 for (Int_t i=0; i<nc; i++) {
550 AliTOFcluster *c=(AliTOFcluster*)clusters->UncheckedAt(i);
551 fClusters[i]=new AliTOFcluster(*c); fN++;
552
553 // Fill Digits QA histos
554
555 Int_t isector = c->GetDetInd(0);
556 Int_t iplate = c->GetDetInd(1);
557 Int_t istrip = c->GetDetInd(2);
558 Int_t ipadX = c->GetDetInd(4);
559 Int_t ipadZ = c->GetDetInd(3);
560
ca270b6d 561 Float_t time =(AliTOFGeometry::TdcBinWidth()*c->GetTDC())*1E-3; // in ns
562 Float_t tot = (AliTOFGeometry::TdcBinWidth()*c->GetToT())*1E-3;//in ns
d686d062 563
564 Int_t stripOffset = 0;
565 switch (iplate) {
566 case 0:
567 stripOffset = 0;
568 break;
569 case 1:
570 stripOffset = nStripC;
571 break;
572 case 2:
573 stripOffset = nStripC+nStripB;
574 break;
575 case 3:
576 stripOffset = nStripC+nStripB+nStripA;
577 break;
578 case 4:
579 stripOffset = nStripC+nStripB+nStripA+nStripB;
580 break;
581 default:
582 AliError(Form("Wrong plate number in TOF (%d) !",iplate));
583 break;
584 };
585 Int_t zindex=npadZ*(istrip+stripOffset)+(ipadZ+1);
586 Int_t phiindex=npadX*isector+ipadX+1;
587 fHDigClusMap->Fill(zindex,phiindex);
588 fHDigClusTime->Fill(time);
589 fHDigClusToT->Fill(tot);
590 }
591
592
593 return 0;
594}
595//_________________________________________________________________________
596void AliTOFtrackerV1::UnloadClusters() {
597 //--------------------------------------------------------------------
598 //This function unloads TOF clusters
599 //--------------------------------------------------------------------
600 for (Int_t i=0; i<fN; i++) {
601 delete fClusters[i];
602 fClusters[i] = 0x0;
603 }
604 fN=0;
605}
606
607//_________________________________________________________________________
608Int_t AliTOFtrackerV1::FindClusterIndex(Double_t z) const {
609 //--------------------------------------------------------------------
610 // This function returns the index of the nearest cluster
611 //--------------------------------------------------------------------
612 //MOD
613 //Here we need to get the Z in the tracking system
614
615 if (fN==0) return 0;
616 if (z <= fClusters[0]->GetZ()) return 0;
617 if (z > fClusters[fN-1]->GetZ()) return fN;
618 Int_t b=0, e=fN-1, m=(b+e)/2;
619 for (; b<e; m=(b+e)/2) {
620 if (z > fClusters[m]->GetZ()) b=m+1;
621 else e=m;
622 }
623 return m;
624}
625
626//_________________________________________________________________________
627Bool_t AliTOFtrackerV1::GetTrackPoint(Int_t index, AliTrackPoint& p) const
628{
629 // Get track space point with index i
630 // Coordinates are in the global system
631 AliTOFcluster *cl = fClusters[index];
632 Float_t xyz[3];
633 cl->GetGlobalXYZ(xyz);
634 Float_t phi=TMath::ATan2(xyz[1],xyz[0]);
635 Float_t phiangle = (Int_t(phi*TMath::RadToDeg()/20.)+0.5)*20.*TMath::DegToRad();
636 Float_t sinphi = TMath::Sin(phiangle), cosphi = TMath::Cos(phiangle);
ca270b6d 637 Float_t tiltangle = AliTOFGeometry::GetAngles(cl->GetDetInd(1),cl->GetDetInd(2))*TMath::DegToRad();
d686d062 638 Float_t sinth = TMath::Sin(tiltangle), costh = TMath::Cos(tiltangle);
ca270b6d 639 Float_t sigmay2 = AliTOFGeometry::XPad()*AliTOFGeometry::XPad()/12.;
640 Float_t sigmaz2 = AliTOFGeometry::ZPad()*AliTOFGeometry::ZPad()/12.;
d686d062 641 Float_t cov[6];
642 cov[0] = sinphi*sinphi*sigmay2 + cosphi*cosphi*sinth*sinth*sigmaz2;
643 cov[1] = -sinphi*cosphi*sigmay2 + sinphi*cosphi*sinth*sinth*sigmaz2;
644 cov[2] = -cosphi*sinth*costh*sigmaz2;
645 cov[3] = cosphi*cosphi*sigmay2 + sinphi*sinphi*sinth*sinth*sigmaz2;
646 cov[4] = -sinphi*sinth*costh*sigmaz2;
647 cov[5] = costh*costh*sigmaz2;
648 p.SetXYZ(xyz[0],xyz[1],xyz[2],cov);
649
650 // Detector numbering scheme
ca270b6d 651 Int_t nSector = AliTOFGeometry::NSectors();
652 Int_t nPlate = AliTOFGeometry::NPlates();
653 Int_t nStripA = AliTOFGeometry::NStripA();
654 Int_t nStripB = AliTOFGeometry::NStripB();
655 Int_t nStripC = AliTOFGeometry::NStripC();
d686d062 656
657 Int_t isector = cl->GetDetInd(0);
658 if (isector >= nSector)
659 AliError(Form("Wrong sector number in TOF (%d) !",isector));
660 Int_t iplate = cl->GetDetInd(1);
661 if (iplate >= nPlate)
662 AliError(Form("Wrong plate number in TOF (%d) !",iplate));
663 Int_t istrip = cl->GetDetInd(2);
664
665 Int_t stripOffset = 0;
666 switch (iplate) {
667 case 0:
668 stripOffset = 0;
669 break;
670 case 1:
671 stripOffset = nStripC;
672 break;
673 case 2:
674 stripOffset = nStripC+nStripB;
675 break;
676 case 3:
677 stripOffset = nStripC+nStripB+nStripA;
678 break;
679 case 4:
680 stripOffset = nStripC+nStripB+nStripA+nStripB;
681 break;
682 default:
683 AliError(Form("Wrong plate number in TOF (%d) !",iplate));
684 break;
685 };
686
687 Int_t idet = (2*(nStripC+nStripB)+nStripA)*isector +
688 stripOffset +
689 istrip;
690 UShort_t volid = AliGeomManager::LayerToVolUID(AliGeomManager::kTOF,idet);
691 p.SetVolumeID((UShort_t)volid);
692 return kTRUE;
693}
694//_________________________________________________________________________
695void AliTOFtrackerV1::InitCheckHists() {
696
697 //Init histos for Digits/Reco QA and Calibration
698
699 //Digits "QA"
700 fHDigClusMap = new TH2F("TOFDig_ClusMap", "",182,0.5,182.5,864, 0.5,864.5);
701 fHDigNClus = new TH1F("TOFDig_NClus", "",200,0.5,200.5);
702 fHDigClusTime = new TH1F("TOFDig_ClusTime", "",2000,0.,200.);
703 fHDigClusToT = new TH1F("TOFDig_ClusToT", "",500,0.,100);
704
705 //Reco "QA"
706 fHRecNClus =new TH1F("TOFRec_NClusW", "",50,0.5,50.5);
707 fHRecDistZ=new TH1F("TOFRec_DistZ", "",50,0.5,10.5);
708 fHRecChi2=new TH1F("TOFRec_Chi2", "",100,0.,10.);
709 fHRecSigYVsP=new TH2F("TOFDig_SigYVsP", "",40,0.,4.,100, 0.,5.);
710 fHRecSigZVsP=new TH2F("TOFDig_SigZVsP", "",40,0.,4.,100, 0.,5.);
711 fHRecSigYVsPWin=new TH2F("TOFDig_SigYVsPWin", "",40,0.,4.,100, 0.,50.);
712 fHRecSigZVsPWin=new TH2F("TOFDig_SigZVsPWin", "",40,0.,4.,100, 0.,50.);
713}
714
715//_________________________________________________________________________
716void AliTOFtrackerV1::SaveCheckHists() {
717
718 //write histos for Digits/Reco QA and Calibration
719
720 TDirectory *dir = gDirectory;
721 TFile *logFile = 0;
722 TFile *logFileTOF = 0;
723
724 TSeqCollection *list = gROOT->GetListOfFiles();
725 int n = list->GetEntries();
726 for(int i=0; i<n; i++) {
727 logFile = (TFile*)list->At(i);
728 if (strstr(logFile->GetName(), "AliESDs.root")) break;
729 }
730
731 Bool_t isThere=kFALSE;
732 for(int i=0; i<n; i++) {
733 logFileTOF = (TFile*)list->At(i);
734 if (strstr(logFileTOF->GetName(), "TOFQA.root")){
735 isThere=kTRUE;
736 break;
737 }
738 }
739
740 logFile->cd();
741 fHDigClusMap->Write(fHDigClusMap->GetName(), TObject::kOverwrite);
742 fHDigNClus->Write(fHDigNClus->GetName(), TObject::kOverwrite);
743 fHDigClusTime->Write(fHDigClusTime->GetName(), TObject::kOverwrite);
744 fHDigClusToT->Write(fHDigClusToT->GetName(), TObject::kOverwrite);
745 fHRecNClus->Write(fHRecNClus->GetName(), TObject::kOverwrite);
746 fHRecChi2->Write(fHRecChi2->GetName(), TObject::kOverwrite);
747 fHRecDistZ->Write(fHRecDistZ->GetName(), TObject::kOverwrite);
748 fHRecSigYVsP->Write(fHRecSigYVsP->GetName(), TObject::kOverwrite);
749 fHRecSigZVsP->Write(fHRecSigZVsP->GetName(), TObject::kOverwrite);
750 fHRecSigYVsPWin->Write(fHRecSigYVsPWin->GetName(), TObject::kOverwrite);
751 fHRecSigZVsPWin->Write(fHRecSigZVsPWin->GetName(), TObject::kOverwrite);
752 logFile->Flush();
753
754 if(!isThere)logFileTOF = new TFile( "TOFQA.root","RECREATE");
755 logFileTOF->cd();
756 fHDigClusMap->Write(fHDigClusMap->GetName(), TObject::kOverwrite);
757 fHDigNClus->Write(fHDigNClus->GetName(), TObject::kOverwrite);
758 fHDigClusTime->Write(fHDigClusTime->GetName(), TObject::kOverwrite);
759 fHDigClusToT->Write(fHDigClusToT->GetName(), TObject::kOverwrite);
760 fHRecNClus->Write(fHRecNClus->GetName(), TObject::kOverwrite);
761 fHRecChi2->Write(fHRecChi2->GetName(), TObject::kOverwrite);
762 fHRecDistZ->Write(fHRecDistZ->GetName(), TObject::kOverwrite);
763 fHRecSigYVsP->Write(fHRecSigYVsP->GetName(), TObject::kOverwrite);
764 fHRecSigZVsP->Write(fHRecSigZVsP->GetName(), TObject::kOverwrite);
765 fHRecSigYVsPWin->Write(fHRecSigYVsPWin->GetName(), TObject::kOverwrite);
766 fHRecSigZVsPWin->Write(fHRecSigZVsPWin->GetName(), TObject::kOverwrite);
767 logFileTOF->Flush();
768
769 dir->cd();
770 }
771//_________________________________________________________________________
772Float_t AliTOFtrackerV1::CorrectTimeWalk( Float_t dist, Float_t tof) {
773
774 //dummy, for the moment
775 Float_t tofcorr=0.;
ca270b6d 776 if(dist<AliTOFGeometry::ZPad()*0.5){
d686d062 777 tofcorr=tof;
778 //place here the actual correction
779 }else{
780 tofcorr=tof;
781 }
782 return tofcorr;
783}
784//_________________________________________________________________________
785Float_t AliTOFtrackerV1::GetTimeZerofromT0(AliESDEvent *event) const {
786
787 //Returns TimeZero as measured by T0 detector
788
789 return event->GetT0();
790}
791//_________________________________________________________________________
792Float_t AliTOFtrackerV1::GetTimeZerofromTOF(AliESDEvent * /*event*/) const {
793
794 //dummy, for the moment. T0 algorithm using tracks on TOF
795 {
796 //place T0 algo here...
797 }
798 return 0.;
799}
128563f6 800//_________________________________________________________________________
801
802void AliTOFtrackerV1::FillClusterArray(TObjArray* arr) const
803{
804 //
805 // Returns the TOF cluster array
806 //
807
808 if (fN==0)
809 arr = 0x0;
810 else
811 for (Int_t i=0; i<fN; ++i) arr->Add(fClusters[i]);
812
813}