ceb444ff49a25096d93eef71d69f7be8095146de
[u/mrichter/AliRoot.git] / STEER / STEER / AliQAv1.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 /* $Id: AliQAv1.cxx 31503 2009-03-16 11:01:16Z schutz $ */
16
17 //////////////////////////////////////////////////////////////////////////////
18 //
19 // Quality Assurance Object//_________________________________________________________________________
20 // Quality Assurance object. The QA status is held in one word per detector,
21 // each bit corresponds to a different status.
22 // bit 0-3   : QA raised during simulation      (RAW)
23 // bit 4-7   : QA raised during simulation      (SIM)
24 // bit 8-11  : QA raised during reconstruction  (REC)
25 // bit 12-15 : QA raised during ESD checking    (ESD)
26 // bit 16-19 : QA raised during analysis        (ANA)
27 // Each of the 4 bits corresponds to a severity level of increasing importance
28 // from lower to higher bit (INFO, WARNING, ERROR, FATAL)
29 //
30 //*-- Yves Schutz CERN, July 2007 
31 //////////////////////////////////////////////////////////////////////////////
32
33
34 #include <cstdlib>
35 // --- ROOT system ---
36 #include <TClass.h>
37 #include <TFile.h>
38 #include <TH1.h>
39 #include <TSystem.h>
40 #include <TROOT.h>
41
42 // --- Standard library ---
43
44 // --- AliRoot header files ---
45 #include "AliLog.h"
46 #include "AliQAv1.h"
47
48
49 ClassImp(AliQAv1)
50 AliQAv1  *     AliQAv1::fgQA                   = 0x0 ;
51 TFile    *     AliQAv1::fgQADataFile           = 0x0 ;   
52 TString        AliQAv1::fgQADataFileName       = "QA" ;  // will transform into Det.QA.run.root  
53 TFile    *     AliQAv1::fgQARefFile            = 0x0 ;   
54 TString        AliQAv1::fgQARefDirName         = "" ; 
55 TString        AliQAv1::fgQARefFileName        = "QA.root" ;
56 TFile    *     AliQAv1::fgQAResultFile         = 0x0 ;  
57 TString        AliQAv1::fgQAResultDirName      = "" ;  
58 TString        AliQAv1::fgQAResultFileName     = "QA.root" ; 
59 TString        AliQAv1::fgDetNames[]           = {"ITS", "TPC", "TRD", "TOF", "PHOS", "HMPID", "EMCAL", "MUON", "FMD",
60                                                   "ZDC", "PMD", "T0", "VZERO", "ACORDE", "HLT", "Global", "CORR", "MFT" } ;   
61 TString        AliQAv1::fgGRPPath              = "GRP/GRP/Data" ; 
62 TString        AliQAv1::fgTaskNames[]          = {"Raws", "Hits", "SDigits", "Digits", "DigitsR", "RecPoints", "TrackSegments", "RecParticles", "ESDs"} ;   
63 TString        AliQAv1::fgModeNames[]          = {"", "Sim", "Rec", "QA"} ;   
64 const TString  AliQAv1::fgkLabLocalFile        = "file://"  ; 
65 const TString  AliQAv1::fgkLabLocalOCDB        = "local://" ;  
66 const TString  AliQAv1::fgkLabAliEnOCDB        = "alien://" ;  
67 const TString  AliQAv1::fgkRefFileName         = "QA.root" ; 
68 const TString  AliQAv1::fgkQAName              = "QA"  ; 
69 const TString  AliQAv1::fgkQACorrNtName        = "CorrQA" ;  
70 const TString  AliQAv1::fgkRefOCDBDirName      = "QA"  ; 
71 TString AliQAv1::fgRefDataDirName                    = ""  ; 
72 const TString  AliQAv1::fgkQARefOCDBDefault    = "alien://folder=/alice/QA/20"  ; 
73 const TString  AliQAv1::fgkExpert              = "Expert" ; 
74 //
75 const UInt_t   AliQAv1::fgkQAClonedBit         = BIT(14); // RS: bits set to the histos or arrays of cloned histos
76 const UInt_t   AliQAv1::fgkForbidCloningBit    = BIT(15); // to flag the state (cloned or not) of the 
77 const UInt_t   AliQAv1::fgkOrigHistoKeptBit    = BIT(16); // object at the given slot of fXXXQAList[]
78 //
79 const UInt_t   AliQAv1::fgkExpertBit           = BIT(19);  
80 const UInt_t   AliQAv1::fgkQABit               = BIT(20) ; 
81 const UInt_t   AliQAv1::fgkImageBit            = BIT(21) ; 
82 const Int_t    AliQAv1::fgkQADebugLevel        = 99 ; 
83 const TString  AliQAv1::fgkImageFileName       = "QAImage" ; 
84 const TString  AliQAv1::fImageFileFormat       = "ps" ; 
85 const UShort_t AliQAv1::fgkMaxQAObjects        = 10000 ; 
86
87 //____________________________________________________________________________
88 AliQAv1::AliQAv1() : 
89   TNamed("", ""), 
90   fNdet(kNDET), 
91   fNEventSpecies(AliRecoParam::kNSpecies), 
92   fLengthQA(fNdet*fNEventSpecies),
93   fQA(new ULong_t[fLengthQA]), 
94   fDet(kNULLDET),
95   fTask(kNULLTASK), 
96   fEventSpecie(AliRecoParam::kDefault), 
97   fEventSpecies(new Bool_t[fNEventSpecies])
98 {
99   // default constructor
100   memset(fQA,0,fLengthQA*sizeof(ULong_t));
101   memset(fEventSpecies,kFALSE,fNEventSpecies*sizeof(Bool_t));
102 }
103
104 //____________________________________________________________________________
105 AliQAv1::AliQAv1(const AliQAv1& qa) :
106   TNamed(qa),
107   fNdet(qa.fNdet), 
108   fNEventSpecies(qa.fNEventSpecies), 
109   fLengthQA(qa.fLengthQA),
110   fQA(new ULong_t[fLengthQA]), 
111   fDet(qa.fDet),
112   fTask(qa.fTask), 
113   fEventSpecie(qa.fEventSpecie), 
114   fEventSpecies(new Bool_t[fNEventSpecies])
115
116   // cpy ctor
117   memcpy(fQA,qa.fQA,fLengthQA*sizeof(ULong_t));
118   memcpy(fEventSpecies,qa.fEventSpecies,fNEventSpecies*sizeof(Bool_t));
119 }
120
121 //_____________________________________________________________________________
122 AliQAv1& AliQAv1::operator = (const AliQAv1& qa)
123 {
124   // assignment operator
125   if(&qa != this) {
126     TNamed::operator=(qa);
127     fNdet          = qa.fNdet;
128     fNEventSpecies = qa.fNEventSpecies; 
129     fLengthQA      = qa.fLengthQA;
130
131     if(fQA) delete [] fQA;
132     fQA = new ULong_t[fLengthQA];
133     memcpy(fQA,qa.fQA,fLengthQA*sizeof(ULong_t));
134
135     fDet = qa.fDet;
136     fTask = qa.fTask;
137     fEventSpecie = qa.fEventSpecie; 
138     if(fEventSpecies) delete [] fEventSpecies;
139     fEventSpecies = new Bool_t[fNEventSpecies];
140     memcpy(fEventSpecies,qa.fEventSpecies,fNEventSpecies*sizeof(Bool_t));
141   }  
142   return *this;
143 }
144
145 //_______________________________________________________________
146 AliQAv1::AliQAv1(const Int_t qalength, ULong_t * qa, const Int_t eslength, Bool_t * es) :
147 TNamed("QA", "Quality Assurance status"),
148 fNdet(kNDET), 
149 fNEventSpecies(eslength), 
150 fLengthQA(qalength),
151 fQA(new ULong_t[fLengthQA]), 
152 fDet(kNULLDET),
153 fTask(kNULLTASK), 
154 fEventSpecie(AliRecoParam::kDefault), 
155 fEventSpecies(new Bool_t[fNEventSpecies])
156 {
157   // constructor to be used
158   memcpy(fQA, qa, fLengthQA*sizeof(ULong_t));
159   memcpy(fEventSpecies, es, fNEventSpecies*sizeof(Bool_t));
160 }
161
162 //_______________________________________________________________
163 AliQAv1::AliQAv1(const DETECTORINDEX_t det) :
164   TNamed("QA", "Quality Assurance status"),
165   fNdet(kNDET), 
166   fNEventSpecies(AliRecoParam::kNSpecies), 
167   fLengthQA(fNdet*fNEventSpecies),
168   fQA(new ULong_t[fLengthQA]), 
169   fDet(det),
170   fTask(kNULLTASK), 
171   fEventSpecie(AliRecoParam::kDefault), 
172   fEventSpecies(new Bool_t[fNEventSpecies])
173 {
174   // constructor to be used
175   if (! CheckRange(det) ) fDet = kNULLDET ; 
176   memset(fQA,0,fLengthQA*sizeof(ULong_t));
177   memset(fEventSpecies,kFALSE,fNEventSpecies*sizeof(Bool_t));
178 }
179   
180 //_______________________________________________________________
181 AliQAv1::AliQAv1(const ALITASK_t tsk) :
182   TNamed("QA", "Quality Assurance status"),
183   fNdet(kNDET), 
184   fNEventSpecies(AliRecoParam::kNSpecies), 
185   fLengthQA(fNdet*fNEventSpecies),
186   fQA(new ULong_t[fLengthQA]), 
187   fDet(kNULLDET),
188   fTask(tsk), 
189   fEventSpecie(AliRecoParam::kDefault), 
190   fEventSpecies(new Bool_t[fNEventSpecies])
191 {
192   // constructor to be used in the AliRoot module (SIM, REC, ESD or ANA)
193   if (! CheckRange(tsk) ) fTask = kNULLTASK ; 
194   memset(fQA,0,fLengthQA*sizeof(ULong_t));
195   memset(fEventSpecies,kFALSE,fNEventSpecies*sizeof(Bool_t));
196 }
197
198 //____________________________________________________________________________
199 AliQAv1::~AliQAv1() 
200 {
201   // dtor  
202   delete [] fQA;
203   delete [] fEventSpecies;
204 }
205
206 //_______________________________________________________________
207 void AliQAv1::Close() 
208 {
209         // close the open files
210         if (fgQADataFile) 
211                 if (fgQADataFile->IsOpen())
212                         fgQADataFile->Close() ; 
213         if (fgQAResultFile) 
214                 if (fgQAResultFile->IsOpen()) 
215                         fgQAResultFile->Close() ;
216         if (fgQARefFile)
217                 if (fgQARefFile->IsOpen())
218                         fgQARefFile->Close() ; 
219
220
221 //_______________________________________________________________
222 Bool_t AliQAv1::CheckFatal() const
223 {
224   // check if any FATAL status is set
225   Bool_t rv = kFALSE ;
226   Int_t index ;
227   for (index = 0; index < kNDET ; index++)
228     rv = rv || IsSet(DETECTORINDEX_t(index), fTask, fEventSpecie, kFATAL) ;
229   return rv ;
230 }
231
232 //_______________________________________________________________
233 Bool_t AliQAv1::CheckRange(DETECTORINDEX_t det) const
234
235   // check if detector is in given detector range: 0-kNDET
236
237   Bool_t rv = ( det < 0 || det > kNDET )  ? kFALSE : kTRUE ;
238   if (!rv)
239     AliFatal(Form("Detector index %d is out of range: 0 <= index <= %d", det, kNDET)) ;
240   return rv ;
241 }
242
243 //_______________________________________________________________
244 Bool_t AliQAv1::CheckRange(ALITASK_t task) const
245
246   // check if task is given taskk range: 0:kNTASK
247   Bool_t rv = ( task < kRAW || task > kNTASK )  ? kFALSE : kTRUE ;
248   if (!rv)
249     AliFatal(Form("Module index %d is out of range: 0 <= index <= %d", task, kNTASK)) ;
250   return rv ;
251 }
252
253 //_______________________________________________________________
254 Bool_t AliQAv1::CheckRange(QABIT_t bit) const
255
256   // check if bit is in given bit range: 0-kNBit
257
258   Bool_t rv = ( bit < 0 || bit > kNBIT )  ? kFALSE : kTRUE ;
259   if (!rv)
260     AliFatal(Form("Status bit %d is out of range: 0 <= bit <= %d", bit, kNBIT)) ;
261   return rv ;
262 }
263
264 //_______________________________________________________________
265 Bool_t AliQAv1::CheckRange(AliRecoParam::EventSpecie_t es) const
266
267   // check if bit is in given bit range: 0-kNBit
268   Bool_t rv = kFALSE ; 
269   switch (es) {
270     case AliRecoParam::kDefault: 
271       rv = kTRUE ; 
272       break ; 
273     case AliRecoParam::kLowMult: 
274       rv = kTRUE ; 
275       break ; 
276     case AliRecoParam::kHighMult: 
277       rv = kTRUE ; 
278       break ; 
279     case AliRecoParam::kCosmic: 
280       rv = kTRUE ; 
281       break ; 
282     case AliRecoParam::kCalib: 
283       rv = kTRUE ; 
284       break ; 
285   }
286   if (!rv)
287     AliFatal(Form("Event Specie %d is not valid", es)) ;
288   return rv ;
289 }
290
291 //_______________________________________________________________
292 const char * AliQAv1::GetAliTaskName(ALITASK_t tsk)
293 {
294         // returns the char name corresponding to module index
295         switch (tsk) {
296                 case kNULLTASK:
297                         break ; 
298                 case kRAW:
299                         return "RAW" ;
300                         break ;  
301                 case kSIM:
302                         return "SIM" ;
303                         break ;
304                 case kREC:
305                         return "REC" ;
306                         break ;
307                 case kESD:
308                         return "ESD" ;
309                         break ;
310                 case kANA:
311                         return "ANA" ;
312                         break ;
313                 default:
314       return "" ; 
315                         break ;
316         }
317   return "" ;
318 }
319
320 //_______________________________________________________________
321 const char * AliQAv1::GetBitName(QABIT_t bit) const
322 {
323         // returns the char name corresponding to bit 
324         TString bitName ;
325         switch (bit) {
326                 case kNULLBit:
327                         bitName = "NONE" ;
328                         break ; 
329                 case kINFO:
330                         bitName = "INFO" ;
331                         break ;  
332                 case kWARNING:
333                         bitName = "WARNING" ;
334                         break ;
335                 case kERROR:
336                         bitName = "ERROR" ;
337                         break ;
338                 case kFATAL:
339                         bitName = "FATAL" ;
340                         break ;
341                 default:
342                         bit = kNULLBit ; 
343                         break ;
344         }
345         return bitName.Data() ;
346 }
347
348 //_______________________________________________________________
349 TH1 * AliQAv1::GetData(TObjArray** list, Int_t index, AliRecoParam::EventSpecie_t eventSpecie)
350 {
351     // retrieve QA data from the list at a given index and for a given event specie 
352   TH1 * rv = NULL ; 
353   Int_t esindex = AliRecoParam::AConvert(eventSpecie) ; 
354   TObjArray * arr = list[esindex] ;
355   if (arr) {
356     if ( index > AliQAv1::GetMaxQAObj() ) {
357                         AliErrorClass(Form("Max number of authorized QA objects is %d", AliQAv1::GetMaxQAObj())) ; 
358                 } else {
359       if ( arr->At(index) )  {
360         rv = static_cast<TH1*>(arr->At(index)) ; 
361       }         
362     }
363   }  
364   return rv ; 
365 }
366
367 //_______________________________________________________________
368 AliQAv1::DETECTORINDEX_t AliQAv1::GetDetIndex(const char * name) 
369 {
370         // returns the detector index corresponding to a given name
371         TString sname(name) ; 
372         DETECTORINDEX_t rv = kNULLDET ; 
373         for (Int_t det = 0; det < kNDET ; det++) {
374                 if ( GetDetName(det) == sname ) {
375                         rv = DETECTORINDEX_t(det) ; 
376                         break ; 
377                 }
378         }
379         return rv ;             
380 }
381
382 //_______________________________________________________________
383 const char * AliQAv1::GetDetName(Int_t det) 
384 {
385         // returns the detector name corresponding to a given index (needed in a loop)
386         
387         if ( det >= 0 &&  det < kNDET) 
388                 return (fgDetNames[det]).Data() ; 
389         else 
390                 return NULL ; 
391 }
392
393 //_______________________________________________________________
394 TFile * AliQAv1::GetQADataFile(const char * name, Int_t run) 
395 {
396   // opens the file to store the detectors Quality Assurance Data Maker results
397         const char * temp = Form("%s.%s.%d.root", name, fgQADataFileName.Data(), run) ; 
398         TString opt ; 
399         if (! fgQADataFile ) {     
400                 if  (gSystem->AccessPathName(temp))
401                         opt = "NEW" ;
402                 else 
403                         opt = "UPDATE" ; 
404                 fgQADataFile = TFile::Open(temp, opt.Data()) ;
405         } else {
406                 if ( strcmp(temp, fgQADataFile->GetName()) != 0 ) {
407                         fgQADataFile = static_cast<TFile *>(gROOT->FindObject(temp)) ;
408                         if ( !fgQADataFile ) {
409           if  (gSystem->AccessPathName(temp))
410             opt = "NEW" ;
411           else 
412             opt = "UPDATE" ; 
413                                 fgQADataFile = TFile::Open(temp, opt.Data()) ;
414                         }
415                 }
416   }
417         return fgQADataFile ;
418
419
420 //_____________________________________________________________________________
421 AliQAv1::MODE_t AliQAv1::Mode(TASKINDEX_t task) {
422   // return "rec" or "sim" depending on the task
423   
424   switch (task) {
425     case AliQAv1::kRAWS:
426       return kRECMODE ; 
427       break;
428     case AliQAv1::kHITS:
429       return kSIMMODE ; 
430       break;
431     case AliQAv1::kSDIGITS:
432       return kSIMMODE ; 
433       break;
434     case AliQAv1::kDIGITS:
435       return kSIMMODE ; 
436       break;
437     case AliQAv1::kDIGITSR:
438       return kRECMODE ; 
439       break;
440     case AliQAv1::kRECPOINTS:
441       return kRECMODE ; 
442       break ; 
443     case AliQAv1::kTRACKSEGMENTS:
444       return kRECMODE ; 
445       break;
446     case AliQAv1::kRECPARTICLES:
447       return kRECMODE ; 
448       break;
449     case AliQAv1::kESDS:
450       return kRECMODE ; 
451       break;
452     default:
453       break;
454   }
455   return AliQAv1::kNULLMODE;
456 }
457
458 //_____________________________________________________________________________
459 TFile * AliQAv1::GetQADataFile(const char * fileName)
460 {
461   // Open if necessary the Data file and return its pointer
462
463   if (!fgQADataFile) {
464     if (!fileName) 
465       fileName = AliQAv1::GetQADataFileName() ; 
466     if  (!gSystem->AccessPathName(fileName)) {
467       fgQADataFile =  TFile::Open(fileName) ;
468     } else {
469       AliFatalClass(Form("File %s not found", fileName)) ;
470     }
471   }
472   return fgQADataFile ; 
473 }
474
475 //_______________________________________________________________
476 TFile * AliQAv1::GetQAResultFile() 
477 {
478   // opens the file to store the  Quality Assurance Data Checker results        
479   if (fgQAResultFile) {
480    if (fgQAResultFile->IsOpen()) 
481     fgQAResultFile->Close();
482    delete fgQAResultFile;
483   }
484   TString dirName(fgQAResultDirName) ; 
485   if ( dirName.Contains(fgkLabLocalFile)) 
486     dirName.ReplaceAll(fgkLabLocalFile, "") ;
487   TString fileName(dirName + fgQAResultFileName) ; 
488   TString opt("") ; 
489   if ( !gSystem->AccessPathName(fileName) )
490     opt = "UPDATE" ; 
491   else { 
492     if ( gSystem->AccessPathName(dirName) )
493       gSystem->mkdir(dirName) ; 
494     opt = "NEW" ; 
495   }
496   fgQAResultFile = TFile::Open(fileName, opt) ;   
497   return fgQAResultFile ; 
498 }
499
500 //_______________________________________________________________
501 AliQAv1::QABIT_t AliQAv1::GetQAStatusBit(AliRecoParam::EventSpecie_t es, DETECTORINDEX_t det, ALITASK_t tsk) const
502 {
503     // returns the QA bit set
504   QABIT_t rv = kNULLBit ; 
505   if ( es == AliRecoParam::kDefault) 
506     es = fEventSpecie ; 
507   if ( det == kNULLDET ) 
508     det = fDet ; 
509   if ( tsk == kNULLTASK ) 
510     tsk = fTask ; 
511   for (Int_t bit = kINFO ; bit < kNBIT ; bit++) {
512                 if (IsSet(det, tsk, es, QABIT_t(bit))) 
513       rv = QABIT_t(bit) ;
514         }
515   return rv ; 
516 }
517
518 //_______________________________________________________________
519 AliQAv1::TASKINDEX_t AliQAv1::GetTaskIndex(const char * name) 
520 {
521         // returns the detector index corresponding to a given name
522         TString sname(name) ; 
523         TASKINDEX_t rv = kNULLTASKINDEX ; 
524         for (Int_t tsk = 0; tsk < kNTASKINDEX ; tsk++) {
525                 if ( GetTaskName(tsk) == sname ) {
526                         rv = TASKINDEX_t(tsk) ; 
527                         break ; 
528                 }
529         }
530         return rv ;             
531 }
532
533 //_______________________________________________________________
534 Bool_t AliQAv1::IsSet(DETECTORINDEX_t det, ALITASK_t tsk, Int_t ies, QABIT_t bit) const 
535 {
536   // Checks is the requested bit is set
537    
538   const AliRecoParam::EventSpecie_t es = AliRecoParam::Convert(ies) ; 
539   return IsSet(det, tsk, es, bit) ; 
540   
541 }  
542
543 //_______________________________________________________________
544 Bool_t AliQAv1::IsSet(DETECTORINDEX_t det, ALITASK_t tsk, AliRecoParam::EventSpecie_t es, QABIT_t bit) const
545 {
546   // Checks is the requested bit is set
547         
548   CheckRange(det) ; 
549   CheckRange(tsk) ;
550   CheckRange(bit) ;
551   CheckRange(es) ;
552         
553   ULong_t offset = Offset(tsk) ;
554   ULong_t status = GetStatus(det, es) ;
555   offset+= bit ;
556   status = (status & 1 << offset) != 0 ;
557   return status ;
558 }
559
560 //_______________________________________________________________
561 Bool_t AliQAv1::IsSetAny(DETECTORINDEX_t det, ALITASK_t tsk, AliRecoParam::EventSpecie_t es) const
562 {
563   // Checks is the requested bit is set
564         
565   CheckRange(det) ; 
566   CheckRange(tsk) ;
567   CheckRange(es) ;
568         
569   ULong_t offset = Offset(tsk) ;
570   ULong_t status = GetStatus(det, es) ;
571         ULong_t st = 0 ; 
572         for ( Int_t bit = 0 ; bit < kNBIT ; bit++) {
573                 offset+= bit ;
574                 st += (status & 1 << offset) != 0 ;             
575         }
576         if ( st == 0 ) 
577                 return kFALSE ; 
578         else 
579                 return kTRUE ;
580 }
581 //_______________________________________________________________
582 Bool_t AliQAv1::IsSetAny(DETECTORINDEX_t det, AliRecoParam::EventSpecie_t es) const
583 {
584   // Checks is the requested bit is set
585         
586   CheckRange(det) ; 
587   CheckRange(es) ; 
588         
589         ULong_t status = GetStatus(det, es) ;
590         ULong_t st = 0 ; 
591         for ( Int_t tsk = 0 ; tsk < kNTASK ; tsk++) {
592                 ULong_t offset = Offset(ALITASK_t(tsk)) ;
593                 for ( Int_t bit = 0 ; bit < kNBIT ; bit++) {
594                         offset+= bit ;
595                         st += (status & 1 << offset) != 0 ;             
596                 }
597         }
598         if ( st == 0 ) 
599                 return kFALSE ; 
600         else 
601                 return kTRUE ;
602 }
603
604 //_______________________________________________________________
605 AliQAv1 * AliQAv1::Instance()
606 {
607   // Get an instance of the singleton. The only authorized way to call the ctor
608
609   if ( ! fgQA) {
610     GetQAResultFile() ; 
611     fgQA = static_cast<AliQAv1 *>(fgQAResultFile->Get("QA")) ; 
612     if ( ! fgQA ) 
613       fgQA = new AliQAv1() ;
614   }
615   return fgQA ;
616 }
617
618 //_______________________________________________________________
619 AliQAv1 * AliQAv1::Instance(const Int_t qalength, ULong_t * qa, const Int_t eslength, Bool_t * es)
620 {
621   // Get an instance of the singleton. The only authorized way to call the ctor
622   
623   if ( ! fgQA) 
624     fgQA = new AliQAv1(qalength, qa, eslength, es) ;
625   return fgQA ;
626 }
627
628 //_______________________________________________________________
629 AliQAv1 * AliQAv1::Instance(const DETECTORINDEX_t det)
630 {
631   // Get an instance of the singleton. The only authorized way to call the ctor
632   
633   if ( ! fgQA) {
634     GetQAResultFile() ; 
635     fgQA = static_cast<AliQAv1 *>(fgQAResultFile->Get(GetQAName())) ; 
636     if ( ! fgQA ) 
637       fgQA = new AliQAv1(det) ;
638   }             
639   fgQA->Set(det) ;
640   return fgQA ;
641 }
642
643 //_______________________________________________________________
644 AliQAv1 * AliQAv1::Instance(const ALITASK_t tsk)
645 {
646   // Get an instance of the singleton. The only authorized way to call the ctor
647
648   if ( ! fgQA)
649     switch (tsk) {
650     case kNULLTASK:
651       break ;
652         case kRAW:
653       fgQA = new AliQAv1(tsk) ;
654       break ;
655         case kSIM:
656       fgQA = new AliQAv1(tsk) ;
657       break ;
658     case kREC:
659       AliInfoClass("fgQA = gAlice->GetQA()") ;
660       break ;
661     case kESD:
662       AliInfoClass("fgQA = static_cast<AliQAv1 *> (esdFile->Get(\"QA\")") ;
663       break ;
664     case kANA:
665       AliInfoClass("fgQA = static_cast<AliQAv1 *> (esdFile->Get(\"QA\")") ;
666       break ;
667     case kNTASK:
668       break ;
669     }
670   if (fgQA) 
671     fgQA->Set(tsk) ;
672   return fgQA ;
673 }
674
675 //_______________________________________________________________
676 AliQAv1 *  AliQAv1::Instance(const TASKINDEX_t tsk) 
677 {
678         // get an instance of the singleton.
679         ALITASK_t index = kNULLTASK ; 
680
681         if ( tsk == kRAWS )
682                 index = kRAW ;
683         else if (tsk < kDIGITS)
684                 index = kSIM ;
685         else if (tsk < kRECPARTICLES)
686                 index = kREC ; 
687         else if (tsk == kESDS) 
688                 index = kESD ; 
689
690         return Instance(index) ; 
691 }
692
693 //_______________________________________________________________
694 Long64_t AliQAv1::Merge(const TCollection * list) {
695         // Merge the QA resuls in the list into this single AliQAv1 object
696         
697   Long64_t nmerge = 0;
698   for (Int_t det = 0 ; det < kNDET ; det++) {
699     Set(DETECTORINDEX_t(det)) ; 
700     for (Int_t task = 0 ; task < kNTASK ; task++) {
701       Set(ALITASK_t(task)) ; 
702       for (Int_t bit = 0 ; bit < kNBIT ; bit++) {
703         TIter next(list) ;
704         AliQAv1 * qa ; 
705         while ( (qa = (AliQAv1*)next() ) ) {
706           for (Int_t es = 0 ; es < fNEventSpecies ; es++) {
707             ++nmerge;
708             if (qa->IsSet(DETECTORINDEX_t(det), ALITASK_t(task), es, QABIT_t(bit)))
709               Set(QABIT_t(bit), es) ; 
710           }
711         } // qa list
712       } // bit
713     } // task
714   } // detector
715   return nmerge;
716 }
717
718 //_______________________________________________________________
719 ULong_t AliQAv1::Offset(ALITASK_t tsk) const
720 {
721   // Calculates the bit offset for a given module (SIM, REC, ESD, ANA)
722
723   CheckRange(tsk) ; 
724
725   ULong_t offset = 0 ;
726   switch (tsk) {
727   case kNULLTASK:
728     break ;
729   case kRAW:
730     offset+= 0 ;
731     break ;
732   case kSIM:
733     offset+= 4 ;
734     break ;
735   case kREC:
736     offset+= 8 ;
737     break ;
738   case kESD:
739     offset+= 12 ;
740     break ;
741   case kANA:
742     offset+= 16 ;
743     break ;
744   case kNTASK:
745     break ;
746   }
747
748   return offset ;
749 }
750
751 //_______________________________________________________________
752 void AliQAv1::Reset(DETECTORINDEX_t det, ALITASK_t tsk, AliRecoParam::EventSpecie_t es)
753 {
754     // resets all the bits for a given tsk and event specie
755   CheckRange(det) ; 
756   CheckRange(tsk) ;
757   CheckRange(es) ;
758         
759   ULong_t offset = Offset(tsk) ;
760   ULong_t status = GetStatus(det, es) ;
761         for ( Int_t bit = 0 ; bit < kNBIT ; bit++) {
762                 offset+= bit ;
763                 status = status & ~1 << offset ;                
764         }  
765   SetStatus(det, es, status) ;
766 }
767
768 //_______________________________________________________________
769 void AliQAv1::ResetStatus(DETECTORINDEX_t det) 
770
771   // reset the status of det for all event specie
772   for (Int_t es = 0 ; es < fNEventSpecies ; es++)
773     fQA[det*fNdet+es] = 0 ; 
774 }
775
776 //_______________________________________________________________
777 void AliQAv1::Set(QABIT_t bit, Int_t ies)
778 {
779   // Set the status bit of the current detector in the current module and for the current event specie 
780    Set(bit, AliRecoParam::Convert(ies)) ;
781 }
782
783 //_______________________________________________________________
784 void AliQAv1::Set(QABIT_t bit, AliRecoParam::EventSpecie_t es)
785 {
786   // Set the status bit of the current detector in the current module and for the current event specie 
787   
788   SetStatusBit(fDet, fTask, es, bit) ;
789 }
790
791 //_____________________________________________________________________________
792 void AliQAv1::SetQARefStorage(const char * name)
793 {
794         // Set the root directory where the QA reference data are stored
795
796         fgQARefDirName = name ; 
797         if ( fgQARefDirName.Contains(fgkLabLocalFile) )
798                 fgQARefFileName =  fgkRefFileName ; 
799         else if ( fgQARefDirName.Contains(fgkLabLocalOCDB) )
800                 fgQARefFileName =  fgkQAName ; 
801         else if ( fgQARefDirName.Contains(fgkLabAliEnOCDB) )
802                 fgQARefFileName =  fgkQAName ; 
803
804   else {
805           AliErrorClass(Form("ERROR: %s is an invalid storage definition\n", name)) ; 
806           fgQARefDirName  = "" ; 
807           fgQARefFileName = "" ; 
808   }     
809         TString tmp(fgQARefDirName) ; // + fgQARefFileName) ;
810         AliInfoClass(Form("AliQAv1::SetQARefDir: QA references are in  %s\n", tmp.Data() )) ;
811 }
812
813 //_____________________________________________________________________________
814 void AliQAv1::SetQAResultDirName(const char * name)
815 {
816   // Set the root directory where to store the QA status object
817
818   fgQAResultDirName.Prepend(name) ; 
819   AliInfoClass(Form("AliQAv1::SetQAResultDirName: QA results are in  %s\n", fgQAResultDirName.Data())) ;
820   if ( fgQAResultDirName.Contains(fgkLabLocalFile)) 
821     fgQAResultDirName.ReplaceAll(fgkLabLocalFile, "") ;
822   fgQAResultFileName.Prepend(fgQAResultDirName) ;
823 }
824
825 //_______________________________________________________________
826 void AliQAv1::SetStatusBit(DETECTORINDEX_t det, ALITASK_t tsk, AliRecoParam::EventSpecie_t es, QABIT_t bit)
827 {
828  // Set the status bit for a given detector and a given task
829
830   CheckRange(det) ;
831   CheckRange(tsk) ;
832   CheckRange(bit) ;
833   CheckRange(es) ;
834   
835   Reset(det, tsk, es) ; 
836   
837   ULong_t offset = Offset(tsk) ;
838   ULong_t status = GetStatus(det, es) ;
839   offset+= bit ;
840   status = status | 1 << offset ;
841   SetStatus(det, es, status) ;
842 }
843
844 //_______________________________________________________________
845 void AliQAv1::Show(DETECTORINDEX_t det) const 
846
847   // dispplay the QA status word
848   if ( det == kNULLDET) 
849     det = fDet ;  
850   for (Int_t ies = 0 ; ies < fNEventSpecies ; ies++) {
851     if ( IsEventSpecieSet(ies) )
852       ShowStatus(det, kNULLTASK, AliRecoParam::ConvertIndex(ies)) ; 
853   }
854 }
855
856 //_______________________________________________________________
857 void AliQAv1::ShowAll() const 
858 {
859   // dispplay the QA status word
860   Int_t index ;
861   for (index = 0 ; index < kNDET ; index++) {
862                 for (Int_t tsk = kRAW ; tsk < kNTASK ; tsk++) {
863       for (Int_t ies = 0 ; ies < fNEventSpecies ; ies++) {
864         if ( IsEventSpecieSet(ies) )
865           ShowStatus(DETECTORINDEX_t(index), ALITASK_t(tsk), AliRecoParam::ConvertIndex(ies)) ;
866       }
867     }
868         }
869 }
870
871 //_______________________________________________________________
872 void AliQAv1::ShowStatus(DETECTORINDEX_t det, ALITASK_t tsk, AliRecoParam::EventSpecie_t es) const
873 {
874         // Prints the full QA status of a given detector
875         CheckRange(det) ;
876         CheckRange(es) ;
877         ULong_t status = GetStatus(det, es) ;
878         ULong_t tskStatus[kNTASK] ; 
879         tskStatus[kRAW] = status & 0x0000f ;
880         tskStatus[kSIM] = status & 0x000f0 ;
881         tskStatus[kREC] = status & 0x00f00 ;
882         tskStatus[kESD] = status & 0x0f000 ;
883         tskStatus[kANA] = status & 0xf0000 ;
884
885         AliInfo(Form("====> QA Status for %8s %8s raw =0x%lx, sim=0x%lx, rec=0x%lx, esd=0x%lx, ana=0x%lx", GetDetName(det).Data(), AliRecoParam::GetEventSpecieName(es), 
886                                  tskStatus[kRAW], tskStatus[kSIM], tskStatus[kREC], tskStatus[kESD], tskStatus[kANA] )) ;
887         if (tsk == kNULLTASK) {
888                 for (Int_t itsk = kRAW ; itsk < kNTASK ; itsk++) {
889                         ShowASCIIStatus(es, det, ALITASK_t(itsk), tskStatus[itsk]) ; 
890                 } 
891         } else {
892                         ShowASCIIStatus(es, det, tsk, tskStatus[tsk]) ; 
893         }
894 }
895
896 //_______________________________________________________________
897 void AliQAv1::ShowASCIIStatus(AliRecoParam::EventSpecie_t es, DETECTORINDEX_t det, ALITASK_t tsk, const ULong_t status) const 
898 {
899         // print the QA status in human readable format
900   const QABIT_t bit = GetQAStatusBit(es, det, tsk) ; 
901   if ( bit != kNULLBit ) {
902           AliInfoClass(Form("           %8s %8s %4s  \n", AliRecoParam::GetEventSpecieName(es), GetDetName(det).Data(), GetAliTaskName(tsk))) ; 
903           AliInfoClass(Form("           %8s %8s %4s 0x%4lx, Problem signalled: %8s \n", AliRecoParam::GetEventSpecieName(es), GetDetName(det).Data(), GetAliTaskName(tsk), status, GetBitName(bit))) ; 
904   }
905 }
906
907 //_______________________________________________________________
908 void AliQAv1::UnSet(QABIT_t bit, Int_t ies)
909 {
910         // UnSet the status bit of the current detector in the current module
911                 UnSet(bit, AliRecoParam::Convert(ies)) ;
912 }
913
914 //_______________________________________________________________
915 void AliQAv1::UnSet(QABIT_t bit, AliRecoParam::EventSpecie_t es)
916 {
917         // UnSet the status bit of the current detector in the current module
918         
919         UnSetStatusBit(fDet, fTask, es, bit) ;
920 }
921
922 //_______________________________________________________________
923 void AliQAv1::UnSetStatusBit(DETECTORINDEX_t det, ALITASK_t tsk, AliRecoParam::EventSpecie_t es, QABIT_t bit)
924 {
925         // UnSet the status bit for a given detector and a given task
926         
927         CheckRange(det) ;
928         CheckRange(tsk) ;
929         CheckRange(bit) ;
930         CheckRange(es) ;
931         
932         ULong_t offset = Offset(tsk) ;
933         ULong_t status = GetStatus(det, es) ;
934         offset+= bit ;
935         status = status & 0 << offset ;
936         SetStatus(det, es, status) ;
937 }