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