CMake: removing qpythia from the depedencies
[u/mrichter/AliRoot.git] / TOF / AliTOFQADataMaker.cxx
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 //  Produces the data needed to calculate the TOF quality assurance. //
19 //  QA objects are 1 & 2 Dimensional histograms.                     //
20 //  author: S.Arcelli                                                //
21 //                                                                   //
22 ///////////////////////////////////////////////////////////////////////
23
24 #include <TClonesArray.h>
25 #include <TFile.h> 
26 #include <TH1I.h> 
27 #include <TH1F.h> 
28 #include <TH2F.h> 
29 #include "AliESDEvent.h"
30 #include "AliESDtrack.h"
31 #include "AliTOFcluster.h"
32 #include "AliTOFdigit.h"
33 #include "AliTOFSDigit.h"
34 #include "AliTOFhitT0.h"
35 #include "AliTOFQADataMaker.h"
36 #include "AliQAChecker.h"
37 #include "AliRawReader.h"
38 #include "AliTOFRawStream.h"
39 #include "AliTOFrawData.h"
40 #include "AliLog.h"
41
42 ClassImp(AliTOFQADataMaker)
43            
44 //____________________________________________________________________________ 
45   AliTOFQADataMaker::AliTOFQADataMaker() : 
46   AliQADataMaker(AliQAv1::GetDetName(AliQAv1::kTOF), "TOF Quality Assurance Data Maker"),
47     fTOFRawStream(AliTOFRawStream())
48 {
49   //
50   // ctor
51   //
52 }
53
54 //____________________________________________________________________________ 
55 AliTOFQADataMaker::AliTOFQADataMaker(const AliTOFQADataMaker& qadm) :
56   AliQADataMaker(),
57   fTOFRawStream(AliTOFRawStream(qadm.fTOFRawStream))
58 {
59   //
60   //copy ctor 
61   //
62   SetName((const char*)qadm.GetName()) ; 
63   SetTitle((const char*)qadm.GetTitle()); 
64 }
65
66 //__________________________________________________________________
67 AliTOFQADataMaker& AliTOFQADataMaker::operator = (const AliTOFQADataMaker& qadm )
68 {
69   //
70   // assignment operator.
71   //
72   this->~AliTOFQADataMaker();
73   new(this) AliTOFQADataMaker(qadm);
74   return *this;
75 }
76  
77 //____________________________________________________________________________ 
78 void AliTOFQADataMaker::InitHits()
79 {
80   //
81   // create Hits histograms in Hits subdir
82   //
83
84   Bool_t expert = kFALSE;
85
86   TH1F * h0 = new TH1F("hTOFHits",    "Number of TOF Hits ",301, -1.02, 5.) ;
87   h0->Sumw2() ;
88   h0->GetXaxis()->SetTitle("TOF hit number [10 power]");
89   Add2HitsList(h0, 0, expert) ;
90
91   TH1F * h1  = new TH1F("hTOFHitsTime", "Hits Time Spectrum in TOF (ns)", 2000, 0., 200) ;
92   h1->Sumw2() ;
93   h1->GetXaxis()->SetTitle("Simulated TOF time [ns]");
94   Add2HitsList(h1, 1, expert) ;
95
96   TH1F * h2  = new TH1F("hTOFHitsLength", "Length Spectrum in TOF (cm)", 500, 0., 500) ;
97   h2->Sumw2() ;
98   h2->GetXaxis()->SetTitle("Track length from primary vertex till hit TOF pad [cm]");
99   Add2HitsList(h2, 2, expert) ;
100
101   TH2F * h3  = new TH2F("hTOFHitsClusMap","Hits vs TOF eta-phi",183, -0.5, 182.5,865,-0.5,864.5) ;
102   h3->Sumw2() ;
103   h3->GetXaxis()->SetTitle("2*strip+padz (eta)");
104   h3->GetYaxis()->SetTitle("48*sector+padx (phi)");
105   h3->GetYaxis()->SetTitleOffset(1.15);
106   Add2HitsList(h3, 3, expert) ;
107   //
108   ClonePerTrigClass(AliQAv1::kHITS); // this should be the last line
109 }
110
111 //____________________________________________________________________________ 
112 void AliTOFQADataMaker::InitDigits()
113 {
114   //
115   // create Digits histograms in Digits subdir
116   //
117
118   Bool_t expert = kFALSE;
119
120   TH1F * h0 = new TH1F("hTOFDigits",    "Number of TOF Digits ",301, -1.02, 5.) ;
121   h0->Sumw2() ;
122   h0->GetXaxis()->SetTitle("TOF digit number [10 power]");
123   Add2DigitsList(h0, 0, expert) ;
124
125   TH1F * h1  = new TH1F("hTOFDigitsTime", "Digits Time Spectrum in TOF (ns)", 2000, 0., 200) ;
126   h1->Sumw2() ;
127   h1->GetXaxis()->SetTitle("Digitized TOF time [ns]");
128   Add2DigitsList(h1, 1, expert) ;
129
130   TH1F * h2  = new TH1F("hTOFDigitsToT", "Digits ToT Spectrum in TOF (ns)", 500, 0., 50) ;
131   h2->Sumw2() ;
132   h2->GetXaxis()->SetTitle("Digitized TOT [ns]");
133   h2->GetYaxis()->SetTitle("Digitized TOF time [ns]");
134   Add2DigitsList(h2, 2, expert) ;
135
136   TH2F * h3  = new TH2F("hTOFDigitsClusMap","Digits vs TOF eta-phi",183, -0.5, 182.5,865,-0.5,864.5) ;
137   h3->Sumw2() ;
138   h3->GetXaxis()->SetTitle("2*strip+padz (eta)");
139   h3->GetYaxis()->SetTitle("48*sector+padx (phi)");
140   h3->GetYaxis()->SetTitleOffset(1.15);
141   Add2DigitsList(h3, 3, expert) ;
142   //
143   ClonePerTrigClass(AliQAv1::kDIGITS); // this should be the last line
144 }
145
146 //____________________________________________________________________________ 
147 void AliTOFQADataMaker::InitSDigits()
148 {
149   //
150   // create SDigits histograms in SDigits subdir
151   //
152
153   Bool_t expert = kFALSE;
154
155   TH1F * h0 = new TH1F("hTOFSDigits",    "Number of TOF SDigits ",301, -1.02, 5.) ;
156   h0->Sumw2() ;
157   h0->GetXaxis()->SetTitle("TOF sdigit number [10 power]");
158   Add2SDigitsList(h0, 0, expert) ;
159
160   TH1F * h1  = new TH1F("hTOFSDigitsTime", "SDigits Time Spectrum in TOF (ns)", 2000, 0., 200) ;
161   h1->Sumw2() ;
162   h1->GetXaxis()->SetTitle("SDigitized TOF time [ns]");
163   Add2SDigitsList(h1, 1, expert) ;
164
165   TH2F * h2  = new TH2F("hTOFSDigitsClusMap","SDigits vs TOF eta-phi",183, -0.5, 182.5,865,-0.5,864.5) ;
166   h2->Sumw2() ;
167   h2->GetXaxis()->SetTitle("2*strip+padz (eta)");
168   h2->GetYaxis()->SetTitle("48*sector+padx (phi)");
169   h2->GetYaxis()->SetTitleOffset(1.15);
170   Add2SDigitsList(h2, 2, expert) ;
171   //
172   ClonePerTrigClass(AliQAv1::kSDIGITS); // this should be the last line
173 }
174
175 //____________________________________________________________________________ 
176 void AliTOFQADataMaker::InitRaws()
177 {
178   //
179   // create Raws histograms in Raws subdir
180   //
181
182   Bool_t expert = kFALSE;
183
184   TH1F * h0 = new TH1F("hTOFRaws",    "Number of TOF Raws ",301, -1.02, 5.) ;
185   h0->Sumw2() ;
186   h0->GetXaxis()->SetTitle("TOF raw number [10 power]");
187   Add2RawsList(h0, 0, expert) ;
188
189   TH1F * h1  = new TH1F("hTOFRawsTime", "Raws Time Spectrum in TOF (ns)", 20000, 0., 2000) ;
190   h1->Sumw2() ;
191   h1->GetXaxis()->SetTitle("Measured TOF time [ns]");
192   Add2RawsList(h1, 1, expert) ;
193
194   TH1F * h2  = new TH1F("hTOFRawsToT", "Raws ToT Spectrum in TOF (ns)", 500, 0., 50) ;
195   h2->Sumw2() ;
196   h2->GetXaxis()->SetTitle("Measured TOT [ns]");
197   h2->GetYaxis()->SetTitle("Measured TOF time [ns]");
198   Add2RawsList(h2, 2, expert) ;
199
200   TH2F * h3  = new TH2F("hTOFRawsClusMap","Raws vs TOF eta-phi",183, -0.5, 182.5,865,-0.5,864.5) ;
201   h3->Sumw2() ;
202   h3->GetXaxis()->SetTitle("2*strip+padz (eta)");
203   h3->GetYaxis()->SetTitle("48*sector+padx (phi)");
204   h3->GetYaxis()->SetTitleOffset(1.15);
205   Add2RawsList(h3, 3, expert) ;
206   //
207   ClonePerTrigClass(AliQAv1::kRAWS); // this should be the last line
208 }
209
210 //____________________________________________________________________________ 
211 void AliTOFQADataMaker::InitRecPoints()
212 {
213   //
214   // create RecPoints histograms in RecPoints subdir
215   //
216
217   Bool_t expert = kFALSE;
218
219   TH1F * h0 = new TH1F("hTOFRecPoints",    "Number of TOF RecPoints ",301, -1.02, 5.) ;
220   h0->Sumw2() ;
221   h0->GetXaxis()->SetTitle("TOF recPoint number [10 power]");
222   Add2RecPointsList(h0, 0, expert) ;
223
224   TH1F * h1  = new TH1F("hTOFRecPointsTime", "RecPoints Time Spectrum in TOF (ns)", 20000, 0., 2000) ;
225   h1->Sumw2() ;
226   h1->GetXaxis()->SetTitle("Calibrated TOF time [ns]");
227   Add2RecPointsList(h1, 1, expert) ;
228
229   TH1F * h2  = new TH1F("hTOFRecPointsRawTime", "RecPoints raw Time Spectrum in TOF (ns)", 20000, 0., 2000) ;
230   h2->Sumw2() ;
231   h2->GetXaxis()->SetTitle("Measured TOT [ns]");
232   h2->GetYaxis()->SetTitle("Measured TOF time [ns]");
233   Add2RecPointsList(h2, 2, expert) ;
234
235   TH1F * h3  = new TH1F("hTOFRecPointsToT", "RecPoints ToT Spectrum in TOF (ns)", 500, 0., 50) ;
236   h3->Sumw2() ;
237   h3->GetXaxis()->SetTitle("Measured TOT [ns]");
238   Add2RecPointsList(h3, 3, expert) ;
239
240   TH2F * h4  = new TH2F("hTOFRecPointsClusMap","RecPoints vs TOF phi-eta",183, -0.5, 182.5,865,-0.5,864.5) ;
241   h4->Sumw2() ;
242   h4->GetXaxis()->SetTitle("2*strip+padz (eta)");
243   h4->GetYaxis()->SetTitle("48*sector+padx (phi)");
244   h4->GetYaxis()->SetTitleOffset(1.15);
245   Add2RecPointsList(h4, 4, expert) ;
246   //
247   ClonePerTrigClass(AliQAv1::kRECPOINTS); // this should be the last line
248 }
249
250 //____________________________________________________________________________ 
251 void AliTOFQADataMaker::InitESDs()
252 {
253   //
254   //create ESDs histograms in ESDs subdir
255   //
256
257   Bool_t expert = kFALSE;
258
259   TH1F * h0 = new TH1F("hTOFESDs",    "Number of matched TOF tracks over ESDs",       250, -1., 4.) ;
260   h0->Sumw2() ; 
261   h0->GetXaxis()->SetTitle("Number of TOF matched ESD tracks [10 power]");
262   Add2ESDsList(h0, 0, expert) ;
263
264   TH1F * h1  = new TH1F("hTOFESDsTime", "Time Spectrum in TOF (ns)", 20000, 0., 2000) ;
265   h1->Sumw2() ;
266   h1->GetXaxis()->SetTitle("Calibrated TOF time [ns]");
267   Add2ESDsList(h1, 1, expert) ;
268
269   TH1F * h2  = new TH1F("hTOFESDsRawTime", "raw Time Spectrum in TOF (ns)", 20000, 0., 2000) ;
270   h2->Sumw2() ;
271   h2->GetXaxis()->SetTitle("Measured TOF time [ns]");
272   Add2ESDsList(h2, 2, expert) ;
273
274   TH1F * h3  = new TH1F("hTOFESDsToT", "ToT Spectrum in TOF (ns)", 500, 0., 50) ;
275   h3->Sumw2() ;
276   h3->GetXaxis()->SetTitle("Measured TOT [ns]");
277   Add2ESDsList(h3, 3, expert) ;
278
279   TH1F * h4 = new TH1F("hTOFESDsPID",    "Fraction of matched TOF tracks with good PID flag", 101, 0., 101.) ;
280   h4->Sumw2() ;
281   h4->GetXaxis()->SetTitle("Fraction of TOF matched ESD tracks with good flag [%]");
282   Add2ESDsList(h4, 4, expert) ;
283   //
284   ClonePerTrigClass(AliQAv1::kESDS); // this should be the last line
285 }
286
287
288 //____________________________________________________________________________
289 void AliTOFQADataMaker::MakeHits(TClonesArray * hits)
290 {
291   //
292   //make QA data from Hits
293   //
294
295   Int_t in[5];
296   Int_t out[5];
297
298   Int_t nentries=hits->GetEntriesFast();
299   if(nentries<=0) {
300     FillHitsData(0,-1.) ; 
301   } else{
302     FillHitsData(0,TMath::Log10(nentries)) ; 
303   }
304   TIter next(hits) ; 
305   AliTOFhitT0 * hit ; 
306   while ( (hit = dynamic_cast<AliTOFhitT0 *>(next())) ) {
307
308     FillHitsData(1, hit->GetTof()*1.E9) ;//in ns
309     FillHitsData(2, hit->GetLen()) ;//in cm
310   
311     in[0] = hit->GetSector();
312     in[1] = hit->GetPlate();
313     in[2]= hit->GetStrip();
314     in[3]= hit->GetPadx();
315     in[4]= hit->GetPadz();
316     GetMapIndeces(in,out);
317     FillHitsData(3, out[0],out[1]) ;//hit map
318   }
319
320 }
321
322
323 //____________________________________________________________________________
324 void AliTOFQADataMaker::MakeHits(TTree * hitTree)
325 {
326   //
327   // make QA data from Hit Tree
328   //
329   if(!hitTree){
330     AliError("can't get the tree with TOF hits !");
331     return;
332   }     
333
334   TBranch * branch = hitTree->GetBranch("TOF") ;
335
336   if (!branch ) {
337     AliError("TOF branch in Hit Tree not found") ; 
338     return;
339   }
340
341   TClonesArray * hits = new TClonesArray("AliTOFhitT0", 1000);
342   TClonesArray * dummy = new TClonesArray("AliTOFhitT0", 1000);
343   branch->SetAddress(&dummy);
344   Int_t index = 0 ;  
345   for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) {
346     branch->GetEntry(ientry) ; 
347     for (Int_t ihit = 0 ; ihit < dummy->GetEntries() ; ihit++) {
348       AliTOFhitT0 * hit = dynamic_cast<AliTOFhitT0 *> (dummy->At(ihit)) ; 
349       new((*hits)[index]) AliTOFhitT0(*hit) ; 
350       index++ ;
351     } 
352   }     
353
354   dummy->Delete();
355   delete dummy;
356   MakeHits(hits) ; 
357   //
358   IncEvCountCycleHits();
359   IncEvCountTotalHits();
360   //
361 }
362
363 //____________________________________________________________________________
364 void AliTOFQADataMaker::MakeDigits(TClonesArray * digits)
365 {
366   //
367   // makes data from Digits
368   //
369   Double_t tdc2ns=AliTOFGeometry::TdcBinWidth()*1E-3;
370   Double_t tot2ns=AliTOFGeometry::ToTBinWidth()*1E-3;
371   Int_t in[5];
372   Int_t out[5];
373
374   Int_t nentries=digits->GetEntriesFast();
375   if(nentries<=0){
376     FillDigitsData(0,-1.) ; 
377   }else{
378     FillDigitsData(0, TMath::Log10(nentries)) ; 
379   } 
380
381   TIter next(digits) ; 
382   AliTOFdigit * digit ; 
383   while ( (digit = dynamic_cast<AliTOFdigit *>(next())) ) {
384     
385     FillDigitsData(1, digit->GetTdc()*tdc2ns) ;//in ns
386     FillDigitsData(2, digit->GetToT()*tot2ns) ;//in ns
387
388     in[0] = digit->GetSector();
389     in[1] = digit->GetPlate();
390     in[2] = digit->GetStrip();
391     in[3] = digit->GetPadx();
392     in[4]= digit->GetPadz();
393     GetMapIndeces(in,out);
394     FillDigitsData(3, out[0],out[1]) ;//digit map
395   }
396
397 }
398
399
400 //____________________________________________________________________________
401 void AliTOFQADataMaker::MakeDigits(TTree * digitTree)
402 {
403   //
404   // makes data from Digit Tree
405   //
406   TClonesArray * digits = new TClonesArray("AliTOFdigit", 1000) ; 
407   
408   TBranch * branch = digitTree->GetBranch("TOF") ;
409   if ( ! branch ) {
410     AliError("TOF branch in Digit Tree not found") ; 
411     return;
412   }
413   branch->SetAddress(&digits) ;
414   branch->GetEntry(0) ; 
415   MakeDigits(digits) ; 
416   //
417   IncEvCountCycleDigits();
418   IncEvCountTotalDigits();
419   //
420 }
421
422 //____________________________________________________________________________
423 void AliTOFQADataMaker::MakeSDigits(TClonesArray * sdigits)
424 {
425   //
426   // makes data from SDigits
427   //
428
429   Double_t tdc2ns=AliTOFGeometry::TdcBinWidth()*1E-3;
430   Int_t in[5];
431   Int_t out[5];
432
433   Int_t nentries=sdigits->GetEntriesFast();
434   if(nentries<=0){
435     FillSDigitsData(0,-1.) ; 
436   }else{
437     FillSDigitsData(0,TMath::Log10(nentries)) ; 
438   } 
439
440   TIter next(sdigits) ; 
441   AliTOFSDigit * sdigit ; 
442   while ( (sdigit = dynamic_cast<AliTOFSDigit *>(next())) ) {
443     
444     for(Int_t i=0;i<sdigit->GetNDigits();i++){
445       FillSDigitsData(1, sdigit->GetTdc(i)*tdc2ns) ;//in ns
446     }
447
448     in[0] = sdigit->GetSector();
449     in[1] = sdigit->GetPlate();
450     in[2] = sdigit->GetStrip();
451     in[3] = sdigit->GetPadx();
452     in[4]= sdigit->GetPadz();
453     GetMapIndeces(in,out);
454     FillSDigitsData(2, out[0],out[1]) ;//sdigit map
455   }
456 }
457
458 //____________________________________________________________________________
459 void AliTOFQADataMaker::MakeSDigits(TTree * sdigitTree)
460 {
461   //
462   // makes data from SDigit Tree
463   //
464   TClonesArray * sdigits = new TClonesArray("AliTOFSDigit", 1000) ; 
465   
466   TBranch * branch = sdigitTree->GetBranch("TOF") ;
467   if ( ! branch ) {
468     AliError("TOF branch in SDigit Tree not found") ; 
469     return;
470   }
471   branch->SetAddress(&sdigits) ;
472   branch->GetEntry(0) ; 
473   MakeSDigits(sdigits) ; 
474   //
475   IncEvCountCycleSDigits();
476   IncEvCountTotalSDigits();
477   //
478 }
479
480 //____________________________________________________________________________
481 void AliTOFQADataMaker::MakeRaws(AliRawReader* rawReader)
482 {
483   //
484   // makes data from Raws
485   //
486
487   Double_t tdc2ns=AliTOFGeometry::TdcBinWidth()*1E-3;
488   Double_t tot2ns=AliTOFGeometry::ToTBinWidth()*1E-3;
489
490
491   Int_t ntof = 0 ; 
492   Int_t in[5];
493   Int_t out[5];
494
495   TClonesArray * clonesRawData;
496   //AliTOFRawStream tofInput(rawReader);
497   fTOFRawStream.SetRawReader(rawReader);
498   for (Int_t iDDL = 0; iDDL < AliTOFGeometry::NDDL()*AliTOFGeometry::NSectors(); iDDL++){
499     rawReader->Reset();
500     fTOFRawStream.LoadRawData(iDDL);
501     clonesRawData = (TClonesArray*)fTOFRawStream.GetRawData();
502     for (Int_t iRawData = 0; iRawData<clonesRawData->GetEntriesFast(); iRawData++) {
503       AliTOFrawData *tofRawDatum = (AliTOFrawData*)clonesRawData->UncheckedAt(iRawData);
504       //if (!tofRawDatum->GetTOT() || !tofRawDatum->GetTOF()) continue;
505       if (!tofRawDatum->GetTOF()) continue;
506       ntof++;
507       FillRawsData(1, tofRawDatum->GetTOF()*tdc2ns) ;//in ns
508       FillRawsData(2, tofRawDatum->GetTOT()*tot2ns) ;//in ns
509
510       fTOFRawStream.EquipmentId2VolumeId(iDDL, 
511                                          tofRawDatum->GetTRM(), 
512                                          tofRawDatum->GetTRMchain(),
513                                          tofRawDatum->GetTDC(), 
514                                          tofRawDatum->GetTDCchannel(), 
515                                          in);
516     
517       GetMapIndeces(in,out);
518       FillRawsData(3, out[0],out[1]) ;//raw map
519       
520     } // while loop
521     
522     clonesRawData->Clear();
523     
524   } // DDL Loop
525   
526   fTOFRawStream.Clear();
527
528   Int_t nentries=ntof;
529   if(nentries<=0){
530     FillRawsData(0,-1.) ; 
531   }else{
532     FillRawsData(0,TMath::Log10(nentries)) ; 
533   }
534   //
535   IncEvCountCycleRaws();
536   IncEvCountTotalRaws();
537   //
538 }
539
540 //____________________________________________________________________________
541 void AliTOFQADataMaker::MakeRecPoints(TTree * clustersTree)
542 {
543   //
544   // Make data from Clusters
545   //
546
547   Double_t tdc2ns=AliTOFGeometry::TdcBinWidth()*1E-3;
548   Double_t tot2ns=AliTOFGeometry::ToTBinWidth()*1E-3;
549
550   Int_t in[5];
551   Int_t out[5];
552
553   TBranch *branch=clustersTree->GetBranch("TOF");
554   if (!branch) { 
555     AliError("can't get the branch with the TOF clusters !");
556     return;
557   }
558
559   TClonesArray dummy("AliTOFcluster",10000), *clusters=&dummy;
560   branch->SetAddress(&clusters);
561
562   // Import the tree
563   clustersTree->GetEvent(0);  
564   
565   Int_t nentries=clusters->GetEntriesFast();
566   if(nentries<=0){
567     FillRecPointsData(0,-1.) ; 
568   }else{
569     FillRecPointsData(0,TMath::Log10(nentries)) ; 
570   } 
571  
572   TIter next(clusters) ; 
573   AliTOFcluster * c ; 
574   while ( (c = dynamic_cast<AliTOFcluster *>(next())) ) {
575     FillRecPointsData(1,c->GetTDC()*tdc2ns);
576     FillRecPointsData(2,c->GetTDCRAW()*tdc2ns);
577     FillRecPointsData(3,c->GetToT()*tot2ns);
578     
579     in[0] = c->GetDetInd(0);
580     in[1] = c->GetDetInd(1);
581     in[2] = c->GetDetInd(2);
582     in[3] = c->GetDetInd(4); //X and Z indeces inverted in RecPoints
583     in[4] = c->GetDetInd(3); //X and Z indeces inverted in RecPoints
584     
585     GetMapIndeces(in,out);
586     FillRecPointsData(4,out[0],out[1]);
587   }
588   //
589   IncEvCountCycleRecPoints();
590   IncEvCountTotalRecPoints();
591   //
592 }
593
594 //____________________________________________________________________________
595 void AliTOFQADataMaker::MakeESDs(AliESDEvent * esd)
596 {
597   //
598   // make QA data from ESDs
599   //  
600   Int_t ntrk = esd->GetNumberOfTracks() ; 
601   Int_t ntof=0;
602   Int_t ntofpid=0;
603   while (ntrk--) {
604     AliESDtrack *track=esd->GetTrack(ntrk);
605     Double_t tofTime=track->GetTOFsignal()*1E-3;//in ns
606     Double_t tofTimeRaw=track->GetTOFsignalRaw()*1E-3;//in ns
607     Double_t tofToT=track->GetTOFsignalToT(); //in ns
608     if(!(tofTime>0))continue;
609     ntof++;
610     FillESDsData(1,tofTime);
611     FillESDsData(2,tofTimeRaw); 
612     FillESDsData(3,tofToT);
613     //check how many tracks where ESD PID is ok 
614     UInt_t status=track->GetStatus();
615     if (((status&AliESDtrack::kESDpid)==0) || 
616         ((status&AliESDtrack::kTOFpid)==0)) continue;
617     ntofpid++;
618   }
619   
620   Int_t nentries=ntof;
621   if(nentries<=0){
622     FillESDsData(0,-1.) ;
623   }else{
624     FillESDsData(0,TMath::Log10(nentries)) ;
625   }
626
627   Float_t ratio = 0.;
628   if(ntof>0) ratio=(Float_t)ntofpid/(Float_t)ntof*100.;
629   FillESDsData(4,ratio) ;
630   //
631   IncEvCountCycleESDs();
632   IncEvCountTotalESDs();
633   //
634 }
635
636 //____________________________________________________________________________ 
637 void AliTOFQADataMaker::StartOfDetectorCycle()
638 {
639   //
640   //Detector specific actions at start of cycle
641   //to be implemented  
642 }
643
644 //____________________________________________________________________________ 
645 void AliTOFQADataMaker::EndOfDetectorCycle(AliQAv1::TASKINDEX task, TObjArray * list)
646 {
647   //Detector specific actions at end of cycle
648   // do the QA checking
649
650   AliQAChecker::Instance()->Run(AliQAv1::kTOF, task, list) ;  
651 }
652 //____________________________________________________________________________
653 void AliTOFQADataMaker::GetMapIndeces(Int_t* in , Int_t* out)
654 {
655   //
656   //return appropriate indeces for the theta-phi map
657   //
658
659   Int_t npadX = AliTOFGeometry::NpadX();
660   Int_t npadZ = AliTOFGeometry::NpadZ();
661   Int_t nStripA = AliTOFGeometry::NStripA();
662   Int_t nStripB = AliTOFGeometry::NStripB();
663   Int_t nStripC = AliTOFGeometry::NStripC();
664
665   Int_t isector = in[0];
666   Int_t iplate = in[1];
667   Int_t istrip = in[2];
668   Int_t ipadX = in[3]; 
669   Int_t ipadZ = in[4]; 
670   
671   Int_t stripOffset = 0;
672   switch (iplate) {
673   case 0:
674     stripOffset = 0;
675       break;
676   case 1:
677     stripOffset = nStripC;
678     break;
679   case 2:
680     stripOffset = nStripC+nStripB;
681     break;
682   case 3:
683     stripOffset = nStripC+nStripB+nStripA;
684     break;
685   case 4:
686     stripOffset = nStripC+nStripB+nStripA+nStripB;
687     break;
688   default:
689     AliError(Form("Wrong plate number in TOF (%d) !",iplate));
690     break;
691   };
692   Int_t zindex=npadZ*(istrip+stripOffset)+(ipadZ+1);
693   Int_t phiindex=npadX*isector+ipadX+1;
694   out[0]=zindex;  
695   out[1]=phiindex;  
696   
697 }