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