]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - STEER/AliQA.cxx
Fixing ompilation warnings
[u/mrichter/AliRoot.git] / STEER / AliQA.cxx
index 69825003cdea6537fd7e1950f6fa1bc755bcf1ba..47365f1c0dfddcae6713d7bfd0af4dc07413b84a 100644 (file)
 // Quality Assurance Object//_________________________________________________________________________
 // Quality Assurance object. The QA status is held in one word per detector,
 // each bit corresponds to a different status.
-// bit 0-3  : QA raised during simulation      (SIM)
-// bit 4-7  : QA raised during reconstruction  (REC)
-// bit 8-11 : QA raised during ESD checking    (ESD)
-// bit 12-15: QA raised during analysis        (ANA)
+// bit 0-3   : QA raised during simulation      (RAW)
+// bit 4-7   : QA raised during simulation      (SIM)
+// bit 8-11  : QA raised during reconstruction  (REC)
+// bit 12-15 : QA raised during ESD checking    (ESD)
+// bit 16-19 : QA raised during analysis        (ANA)
 // Each of the 4 bits corresponds to a severity level of increasing importance
 // from lower to higher bit (INFO, WARNING, ERROR, FATAL)
 //
@@ -31,6 +32,7 @@
 //////////////////////////////////////////////////////////////////////////////
 
 
+#include <cstdlib>
 // --- ROOT system ---
 #include <TFile.h>
 #include <TSystem.h>
 ClassImp(AliQA)
 AliQA    * AliQA::fgQA                   = 0x0 ;
 TFile    * AliQA::fgQADataFile           = 0x0 ;   
-TString    AliQA::fgQADataFileName       = "QA" ;  // will transform into Det.QA.run.cycle.root  
+TString    AliQA::fgQADataFileName       = "QA" ;  // will transform into Det.QA.run.root  
 TFile    * AliQA::fgQARefFile            = 0x0 ;   
-TString    AliQA::fgQARefDirName            = "" ; 
+TString    AliQA::fgQARefDirName              = "" ; 
 TString    AliQA::fgQARefFileName        = "QA.root" ;
 TFile    * AliQA::fgQAResultFile         = 0x0 ;  
 TString    AliQA::fgQAResultDirName      = "" ;  
 TString    AliQA::fgQAResultFileName     = "QA.root" ; 
 TString    AliQA::fgDetNames[]           = {"ITS", "TPC", "TRD", "TOF", "PHOS", "HMPID", "EMCAL", "MUON", "FMD",
-                                                                                   "ZDC", "PMD", "T0", "VZERO", "ACORDE", "HLT", "Global"} ;   
+                                            "ZDC", "PMD", "T0", "VZERO", "ACORDE", "HLT", "Global", "CORR"} ;   
+TString    AliQA::fgGRPPath              = "GRP/GRP/Data" ; 
 TString       AliQA::fgTaskNames[]       = {"Raws", "Hits", "SDigits", "Digits", "RecPoints", "TrackSegments", "RecParticles", "ESDs"} ;   
-const TString AliQA::fkgLabLocalFile     = "file://"  ; 
-const TString AliQA::fkgLabLocalOCDB     = "local://" ;  
-const TString AliQA::fkgLabAliEnOCDB     = "alien://" ;  
-const TString AliQA::fkgRefFileName      = "QA.root" ; 
-const TString AliQA::fkgQAName           = "QA"  ; 
-const TString AliQA::fkgRefOCDBDirName   = "Ref"  ; 
-TString AliQA::fkgRefDataDirName   = "Data"  ; 
-const TString AliQA::fkgQARefOCDBDefault = "alien://folder=/alice/QA/20"  ; 
+const TString AliQA::fgkLabLocalFile     = "file://"  ; 
+const TString AliQA::fgkLabLocalOCDB     = "local://" ;  
+const TString AliQA::fgkLabAliEnOCDB     = "alien://" ;  
+const TString AliQA::fgkRefFileName      = "QA.root" ; 
+const TString AliQA::fgkQAName           = "QA"  ; 
+const TString AliQA::fgkQACorrNtName     = "CorrQA" ;  
+const TString AliQA::fgkRefOCDBDirName   = "QA"  ; 
+TString AliQA::fgRefDataDirName                 = ""  ; 
+const TString AliQA::fgkQARefOCDBDefault = "alien://folder=/alice/QA/20"  ; 
+const TString AliQA::fgkExpert           = "Expert" ; 
+const UInt_t  AliQA::fgkExpertBit        = 16 ; 
+const UInt_t  AliQA::fgkQABit            = 17 ; 
+
 //____________________________________________________________________________
 AliQA::AliQA() : 
   TNamed("", ""), 
   fNdet(kNDET), 
-  fQA(new ULong_t[fNdet]), 
+  fNEventSpecies(AliRecoParam::kNSpecies), 
+  fLengthQA(fNdet*fNEventSpecies),
+  fQA(new ULong_t[fLengthQA]), 
   fDet(kNULLDET),
-  fTask(kNULLTASK)
-       
+  fTask(kNULLTASK), 
+  fEventSpecie(AliRecoParam::kDefault), 
+  fEventSpecies(new Bool_t[fNEventSpecies])
 {
   // default constructor
-  // beware singleton: not to be used
-  
-  for (Int_t index = 0 ; index < fNdet ; index++) 
-       fQA[index] = 0 ; 
+  memset(fQA,0,fLengthQA*sizeof(ULong_t));
+  memset(fEventSpecies,kFALSE,fNEventSpecies*sizeof(Bool_t));
 }
 
 //____________________________________________________________________________
 AliQA::AliQA(const AliQA& qa) :
   TNamed(qa),
   fNdet(qa.fNdet), 
-  fQA(qa.fQA), 
+  fNEventSpecies(qa.fNEventSpecies), 
+  fLengthQA(qa.fLengthQA),
+  fQA(new ULong_t[fLengthQA]), 
   fDet(qa.fDet),
-  fTask(qa.fTask)
+  fTask(qa.fTask), 
+  fEventSpecie(qa.fEventSpecie), 
+  fEventSpecies(new Bool_t[fNEventSpecies])
 { 
   // cpy ctor
+  memcpy(fQA,qa.fQA,fLengthQA*sizeof(ULong_t));
+  memcpy(fEventSpecies,qa.fEventSpecies,fNEventSpecies*sizeof(Bool_t));
 }
 
 //_____________________________________________________________________________
 AliQA& AliQA::operator = (const AliQA& qa)
 {
-// assignment operator
-
-  this->~AliQA();
-  new(this) AliQA(qa);
+  // assignment operator
+  if(&qa != this) {
+    TNamed::operator=(qa);
+    fNdet          = qa.fNdet;
+    fNEventSpecies = qa.fNEventSpecies; 
+    fLengthQA      = qa.fLengthQA;
+
+    if(fQA) delete [] fQA;
+    fQA = new ULong_t[fLengthQA];
+    memcpy(fQA,qa.fQA,fLengthQA*sizeof(ULong_t));
+
+    fDet = qa.fDet;
+    fTask = qa.fTask;
+    fEventSpecie = qa.fEventSpecie; 
+    if(fEventSpecies) delete [] fEventSpecies;
+    fEventSpecies = new Bool_t[fNEventSpecies];
+    memcpy(fEventSpecies,qa.fEventSpecies,fNEventSpecies*sizeof(Bool_t));
+  }  
   return *this;
 }
 
+//_______________________________________________________________
+AliQA::AliQA(const Int_t qalength, ULong_t * qa, const Int_t eslength, Bool_t * es) :
+TNamed("QA", "Quality Assurance status"),
+fNdet(kNDET), 
+fNEventSpecies(eslength), 
+fLengthQA(qalength),
+fQA(new ULong_t[fLengthQA]), 
+fDet(kNULLDET),
+fTask(kNULLTASK), 
+fEventSpecie(AliRecoParam::kDefault), 
+fEventSpecies(new Bool_t[fNEventSpecies])
+{
+  // constructor to be used
+  memcpy(fQA, qa, fLengthQA*sizeof(ULong_t));
+  memcpy(fEventSpecies, es, fNEventSpecies*sizeof(Bool_t));
+}
+
 //_______________________________________________________________
 AliQA::AliQA(const DETECTORINDEX_t det) :
   TNamed("QA", "Quality Assurance status"),
-  fNdet(kNDET),  
-  fQA(new ULong_t[fNdet]), 
+  fNdet(kNDET), 
+  fNEventSpecies(AliRecoParam::kNSpecies), 
+  fLengthQA(fNdet*fNEventSpecies),
+  fQA(new ULong_t[fLengthQA]), 
   fDet(det),
-  fTask(kNULLTASK) 
+  fTask(kNULLTASK), 
+  fEventSpecie(AliRecoParam::kDefault), 
+  fEventSpecies(new Bool_t[fNEventSpecies])
 {
   // constructor to be used
-  if (! CheckRange(det) ) {
-    fDet = kNULLDET ; 
-    return ;
-  } 
-  Int_t index ; 
-  for (index = 0; index < fNdet; index++) 
-    fQA[index] = 0 ; 
+  if (! CheckRange(det) ) fDet = kNULLDET ; 
+  memset(fQA,0,fLengthQA*sizeof(ULong_t));
+  memset(fEventSpecies,kFALSE,fNEventSpecies*sizeof(Bool_t));
 }
   
 //_______________________________________________________________
 AliQA::AliQA(const ALITASK_t tsk) :
-  TNamed("QA", "Quality Assurance status"), 
-  fNdet(kNDET),
-  fQA(new ULong_t[fNdet]), 
+  TNamed("QA", "Quality Assurance status"),
+  fNdet(kNDET), 
+  fNEventSpecies(AliRecoParam::kNSpecies), 
+  fLengthQA(fNdet*fNEventSpecies),
+  fQA(new ULong_t[fLengthQA]), 
   fDet(kNULLDET),
-  fTask(tsk)
+  fTask(tsk), 
+  fEventSpecie(AliRecoParam::kDefault), 
+  fEventSpecies(new Bool_t[fNEventSpecies])
 {
   // constructor to be used in the AliRoot module (SIM, REC, ESD or ANA)
-  if (! CheckRange(tsk) ) {
-    fTask = kNULLTASK ; 
-    return ;
-  } 
-  Int_t index ; 
-  for (index = 0; index < fNdet; index++) 
-    fQA[index] = 0 ; 
+  if (! CheckRange(tsk) ) fTask = kNULLTASK ; 
+  memset(fQA,0,fLengthQA*sizeof(ULong_t));
+  memset(fEventSpecies,kFALSE,fNEventSpecies*sizeof(Bool_t));
 }
 
 //____________________________________________________________________________
 AliQA::~AliQA() 
 {
   // dtor  
-  delete[] fQA ;
+  delete [] fQA;
+  delete [] fEventSpecies;
 }
 
 //_______________________________________________________________
@@ -160,18 +207,18 @@ void AliQA::Close()
 } 
 
 //_______________________________________________________________
-const Bool_t AliQA::CheckFatal() const
+Bool_t AliQA::CheckFatal() const
 {
   // check if any FATAL status is set
   Bool_t rv = kFALSE ;
   Int_t index ;
   for (index = 0; index < kNDET ; index++)
-    rv = rv || IsSet(DETECTORINDEX_t(index), fTask, kFATAL) ;
+    rv = rv || IsSet(DETECTORINDEX_t(index), fTask, fEventSpecie, kFATAL) ;
   return rv ;
 }
 
 //_______________________________________________________________
-const Bool_t AliQA::CheckRange(DETECTORINDEX_t det) const
+Bool_t AliQA::CheckRange(DETECTORINDEX_t det) const
 { 
   // check if detector is in given detector range: 0-kNDET
 
@@ -182,7 +229,7 @@ const Bool_t AliQA::CheckRange(DETECTORINDEX_t det) const
 }
 
 //_______________________________________________________________
-const Bool_t AliQA::CheckRange(ALITASK_t task) const
+Bool_t AliQA::CheckRange(ALITASK_t task) const
 { 
   // check if task is given taskk range: 0:kNTASK
   Bool_t rv = ( task < kRAW || task > kNTASK )  ? kFALSE : kTRUE ;
@@ -192,7 +239,7 @@ const Bool_t AliQA::CheckRange(ALITASK_t task) const
 }
 
 //_______________________________________________________________
-const Bool_t AliQA::CheckRange(QABIT_t bit) const
+Bool_t AliQA::CheckRange(QABIT_t bit) const
 { 
   // check if bit is in given bit range: 0-kNBit
 
@@ -202,40 +249,92 @@ const Bool_t AliQA::CheckRange(QABIT_t bit) const
   return rv ;
 }
 
-
+//_______________________________________________________________
+Bool_t AliQA::CheckRange(AliRecoParam::EventSpecie_t es) const
+{ 
+  // check if bit is in given bit range: 0-kNBit
+  Bool_t rv = kFALSE ; 
+  switch (es) {
+    case AliRecoParam::kDefault: 
+      rv = kTRUE ; 
+      break ; 
+    case AliRecoParam::kLowMult: 
+      rv = kTRUE ; 
+      break ; 
+    case AliRecoParam::kHighMult: 
+      rv = kTRUE ; 
+      break ; 
+    case AliRecoParam::kCosmic: 
+      rv = kTRUE ; 
+      break ; 
+    case AliRecoParam::kCalib: 
+      rv = kTRUE ; 
+      break ; 
+  }
+  if (!rv)
+    AliFatal(Form("Event Specie %d is not valid", es)) ;
+  return rv ;
+}
 
 //_______________________________________________________________
 const char * AliQA::GetAliTaskName(ALITASK_t tsk)
 {
-  // returns the char name corresponding to module index
-  TString tskName ;
-  switch (tsk) {
-  case kNULLTASK:
-    break ; 
-  case kRAW:
-    tskName = "RAW" ;
-    break ;  
-  case kSIM:
-    tskName = "SIM" ;
-    break ;
-  case kREC:
-    tskName = "REC" ;
-    break ;
-  case kESD:
-    tskName = "ESD" ;
-    break ;
-  case kANA:
-    tskName = "ANA" ;
-    break ;
-  default:
-    tsk = kNULLTASK ; 
-    break ;
-  }
-  return tskName.Data() ;
+       // returns the char name corresponding to module index
+       TString tskName ;
+       switch (tsk) {
+               case kNULLTASK:
+                       break ; 
+               case kRAW:
+                       tskName = "RAW" ;
+                       break ;  
+               case kSIM:
+                       tskName = "SIM" ;
+                       break ;
+               case kREC:
+                       tskName = "REC" ;
+                       break ;
+               case kESD:
+                       tskName = "ESD" ;
+                       break ;
+               case kANA:
+                       tskName = "ANA" ;
+                       break ;
+               default:
+                       tsk = kNULLTASK ; 
+                       break ;
+       }
+       return tskName.Data() ;
+}
+
+//_______________________________________________________________
+const char * AliQA::GetBitName(QABIT_t bit) const
+{
+       // returns the char name corresponding to bit 
+       TString bitName ;
+       switch (bit) {
+               case kNULLBit:
+                       break ; 
+               case kINFO:
+                       bitName = "INFO" ;
+                       break ;  
+               case kWARNING:
+                       bitName = "WARNING" ;
+                       break ;
+               case kERROR:
+                       bitName = "ERROR" ;
+                       break ;
+               case kFATAL:
+                       bitName = "FATAL" ;
+                       break ;
+               default:
+                       bit = kNULLBit ; 
+                       break ;
+       }
+       return bitName.Data() ;
 }
 
 //_______________________________________________________________
-const AliQA::DETECTORINDEX_t AliQA::GetDetIndex(const char * name) 
+AliQA::DETECTORINDEX_t AliQA::GetDetIndex(const char * name) 
 {
        // returns the detector index corresponding to a given name
        TString sname(name) ; 
@@ -261,10 +360,10 @@ const char * AliQA::GetDetName(Int_t det)
 }
 
 //_______________________________________________________________
-TFile * AliQA::GetQADataFile(const char * name, const Int_t run, const Int_t cycle
+TFile * AliQA::GetQADataFile(const char * name, Int_t run
 {
   // opens the file to store the detectors Quality Assurance Data Maker results
-       const char * temp = Form("%s.%s.%d.%d.root", name, fgQADataFileName.Data(), run, cycle) ; 
+       const char * temp = Form("%s.%s.%d.root", name, fgQADataFileName.Data(), run) ; 
        TString opt ; 
        if (! fgQADataFile ) {     
                if  (gSystem->AccessPathName(temp))
@@ -313,8 +412,8 @@ TFile * AliQA::GetQAResultFile()
        fgQAResultFile = 0x0 ; 
 //     if (!fgQAResultFile) { 
                TString dirName(fgQAResultDirName) ; 
-               if ( dirName.Contains(fkgLabLocalFile)) 
-                       dirName.ReplaceAll(fkgLabLocalFile, "") ;
+               if ( dirName.Contains(fgkLabLocalFile)) 
+                       dirName.ReplaceAll(fgkLabLocalFile, "") ;
                TString fileName(dirName + fgQAResultFileName) ; 
                TString opt("") ; 
                if ( !gSystem->AccessPathName(fileName) )
@@ -331,27 +430,112 @@ TFile * AliQA::GetQAResultFile()
 }
 
 //_______________________________________________________________
-const Bool_t AliQA::IsSet(DETECTORINDEX_t det, ALITASK_t tsk, QABIT_t bit) const
+AliQA::TASKINDEX_t AliQA::GetTaskIndex(const char * name) 
+{
+       // returns the detector index corresponding to a given name
+       TString sname(name) ; 
+       TASKINDEX_t rv = kNULLTASKINDEX ; 
+       for (Int_t tsk = 0; tsk < kNTASKINDEX ; tsk++) {
+               if ( GetTaskName(tsk) == sname ) {
+                       rv = TASKINDEX_t(tsk) ; 
+                       break ; 
+               }
+       }
+       return rv ;             
+}
+
+//_______________________________________________________________
+Bool_t AliQA::IsSet(DETECTORINDEX_t det, ALITASK_t tsk, Int_t ies, QABIT_t bit) const 
 {
   // Checks is the requested bit is set
+   
+  const AliRecoParam::EventSpecie_t es = AliRecoParam::Convert(ies) ; 
+  return IsSet(det, tsk, es, bit) ; 
+  
+}  
 
+//_______________________________________________________________
+Bool_t AliQA::IsSet(DETECTORINDEX_t det, ALITASK_t tsk, AliRecoParam::EventSpecie_t es, QABIT_t bit) const
+{
+  // Checks is the requested bit is set
+       
   CheckRange(det) ; 
   CheckRange(tsk) ;
   CheckRange(bit) ;
-
+  CheckRange(es) ;
+       
   ULong_t offset = Offset(tsk) ;
-  ULong_t status = GetStatus(det) ;
+  ULong_t status = GetStatus(det, es) ;
   offset+= bit ;
   status = (status & 1 << offset) != 0 ;
   return status ;
 }
 
+//_______________________________________________________________
+Bool_t AliQA::IsSetAny(DETECTORINDEX_t det, ALITASK_t tsk, AliRecoParam::EventSpecie_t es) const
+{
+  // Checks is the requested bit is set
+       
+  CheckRange(det) ; 
+  CheckRange(tsk) ;
+  CheckRange(es) ;
+       
+  ULong_t offset = Offset(tsk) ;
+  ULong_t status = GetStatus(det, es) ;
+       ULong_t st = 0 ; 
+       for ( Int_t bit = 0 ; bit < kNBIT ; bit++) {
+               offset+= bit ;
+               st += (status & 1 << offset) != 0 ;             
+       }
+       if ( st == 0 ) 
+               return kFALSE ; 
+       else 
+               return kTRUE ;
+}
+//_______________________________________________________________
+Bool_t AliQA::IsSetAny(DETECTORINDEX_t det, AliRecoParam::EventSpecie_t es) const
+{
+  // Checks is the requested bit is set
+       
+  CheckRange(det) ; 
+  CheckRange(es) ; 
+       
+       ULong_t status = GetStatus(det, es) ;
+       ULong_t st = 0 ; 
+       for ( Int_t tsk = 0 ; tsk < kNTASK ; tsk++) {
+               ULong_t offset = Offset(ALITASK_t(tsk)) ;
+               for ( Int_t bit = 0 ; bit < kNBIT ; bit++) {
+                       offset+= bit ;
+                       st += (status & 1 << offset) != 0 ;             
+               }
+       }
+       if ( st == 0 ) 
+               return kFALSE ; 
+       else 
+               return kTRUE ;
+}
+
 //_______________________________________________________________
 AliQA * AliQA::Instance()
 {
-  // Get an instance of the singleton.
-  // Object must have been instantiated with Instance(ALITASK) first
+  // Get an instance of the singleton. The only authorized way to call the ctor
+
+  if ( ! fgQA) {
+    TFile * f = GetQAResultFile() ; 
+    fgQA = dynamic_cast<AliQA *>(f->Get("QA")) ; 
+    if ( ! fgQA ) 
+      fgQA = new AliQA() ;
+  }    
+  return fgQA ;
+}
 
+//_______________________________________________________________
+AliQA * AliQA::Instance(const Int_t qalength, ULong_t * qa, const Int_t eslength, Bool_t * es)
+{
+  // Get an instance of the singleton. The only authorized way to call the ctor
+  
+  if ( ! fgQA) 
+    fgQA = new AliQA(qalength, qa, eslength, es) ;
   return fgQA ;
 }
 
@@ -373,7 +557,7 @@ AliQA * AliQA::Instance(const DETECTORINDEX_t det)
 //_______________________________________________________________
 AliQA * AliQA::Instance(const ALITASK_t tsk)
 {
-  // get an instance of the singleton.
+  // Get an instance of the singleton. The only authorized way to call the ctor
 
   if ( ! fgQA)
     switch (tsk) {
@@ -422,7 +606,29 @@ AliQA *  AliQA::Instance(const TASKINDEX_t tsk)
 }
 
 //_______________________________________________________________
-const ULong_t AliQA::Offset(ALITASK_t tsk) const
+void AliQA::Merge(TCollection * list) {
+       // Merge the QA resuls in the list into this single AliQA object
+       
+       for (Int_t det = 0 ; det < kNDET ; det++) {
+               Set(DETECTORINDEX_t(det)) ; 
+               for (Int_t task = 0 ; task < kNTASK ; task++) {
+                       Set(ALITASK_t(task)) ; 
+                       for (Int_t bit = 0 ; bit < kNBIT ; bit++) {
+                               TIter next(list) ;
+                               AliQA * qa ; 
+                               while ( (qa = (AliQA*)next() ) ) {
+          for (Int_t es = 0 ; es < fNEventSpecies ; es++) {
+            if (qa->IsSet(DETECTORINDEX_t(det), ALITASK_t(task), es, QABIT_t(bit)))
+              Set(QABIT_t(bit), es) ; 
+          }
+                               } // qa list
+                       } // bit
+               } // task
+       } // detector
+}
+
+//_______________________________________________________________
+ULong_t AliQA::Offset(ALITASK_t tsk) const
 {
   // Calculates the bit offset for a given module (SIM, REC, ESD, ANA)
 
@@ -455,11 +661,26 @@ const ULong_t AliQA::Offset(ALITASK_t tsk) const
 }
 
 //_______________________________________________________________
-void AliQA::Set(QABIT_t bit)
+void AliQA::ResetStatus(DETECTORINDEX_t det) 
+{ 
+  // reset the status of det for all event specie
+  for (Int_t es = 0 ; es < fNEventSpecies ; es++)
+    fQA[det*fNdet+es] = 0 ; 
+}
+
+//_______________________________________________________________
+void AliQA::Set(QABIT_t bit, Int_t ies)
+{
+  // Set the status bit of the current detector in the current module and for the current event specie 
+   Set(bit, AliRecoParam::Convert(ies)) ;
+}
+
+//_______________________________________________________________
+void AliQA::Set(QABIT_t bit, AliRecoParam::EventSpecie_t es)
 {
-  // Set the status bit of the current detector in the current module
+  // Set the status bit of the current detector in the current module and for the current event specie 
   
-  SetStatusBit(fDet, fTask, bit) ;
+  SetStatusBit(fDet, fTask, es, bit) ;
 }
 
 //_____________________________________________________________________________
@@ -468,12 +689,12 @@ void AliQA::SetQARefStorage(const char * name)
        // Set the root directory where the QA reference data are stored
 
        fgQARefDirName = name ; 
-       if ( fgQARefDirName.Contains(fkgLabLocalFile) )
-               fgQARefFileName =  fkgRefFileName ; 
-       else if ( fgQARefDirName.Contains(fkgLabLocalOCDB) )
-               fgQARefFileName =  fkgQAName ; 
-       else if ( fgQARefDirName.Contains(fkgLabAliEnOCDB) )
-               fgQARefFileName =  fkgQAName ; 
+       if ( fgQARefDirName.Contains(fgkLabLocalFile) )
+               fgQARefFileName =  fgkRefFileName ; 
+       else if ( fgQARefDirName.Contains(fgkLabLocalOCDB) )
+               fgQARefFileName =  fgkQAName ; 
+       else if ( fgQARefDirName.Contains(fgkLabAliEnOCDB) )
+               fgQARefFileName =  fgkQAName ; 
 
   else {
          printf("ERROR: %s is an invalid storage definition\n", name) ; 
@@ -491,71 +712,125 @@ void AliQA::SetQAResultDirName(const char * name)
 
   fgQAResultDirName.Prepend(name) ; 
   printf("AliQA::SetQAResultDirName: QA results are in  %s\n", fgQAResultDirName.Data()) ;
-  if ( fgQAResultDirName.Contains(fkgLabLocalFile)) 
-    fgQAResultDirName.ReplaceAll(fkgLabLocalFile, "") ;
+  if ( fgQAResultDirName.Contains(fgkLabLocalFile)) 
+    fgQAResultDirName.ReplaceAll(fgkLabLocalFile, "") ;
   fgQAResultFileName.Prepend(fgQAResultDirName) ;
 }
 
 //_______________________________________________________________
-void AliQA::SetStatusBit(DETECTORINDEX_t det, ALITASK_t tsk, QABIT_t bit)
+void AliQA::SetStatusBit(DETECTORINDEX_t det, ALITASK_t tsk, AliRecoParam::EventSpecie_t es, QABIT_t bit)
 {
  // Set the status bit for a given detector and a given task
 
   CheckRange(det) ;
   CheckRange(tsk) ;
   CheckRange(bit) ;
+  CheckRange(es) ;
 
   ULong_t offset = Offset(tsk) ;
-  ULong_t status = GetStatus(det) ;
+  ULong_t status = GetStatus(det, es) ;
   offset+= bit ;
   status = status | 1 << offset ;
-  SetStatus(det, status) ;
+  SetStatus(det, es, status) ;
+}
+
+//_______________________________________________________________
+void AliQA::Show(DETECTORINDEX_t det) const 
+{ 
+  // dispplay the QA status word
+  if ( det == kNULLDET) 
+    det = fDet ;  
+  for (Int_t ies = 0 ; ies < fNEventSpecies ; ies++) {
+    const Bool_t what = IsEventSpecieSet(ies) ;
+    if ( what )
+      ShowStatus(det, kNULLTASK, AliRecoParam::Convert(ies)) ; 
+  }
 }
 
 //_______________________________________________________________
-void AliQA::ShowAll() const
+void AliQA::ShowAll() const 
 {
   // dispplay the QA status word
   Int_t index ;
-  for (index = 0 ; index < kNDET ; index++)
-    ShowStatus(DETECTORINDEX_t(index)) ;
+  for (index = 0 ; index < kNDET ; index++) {
+               for (Int_t tsk = kRAW ; tsk < kNTASK ; tsk++) {
+      for (Int_t ies = 0 ; ies < fNEventSpecies ; ies++) {
+        const Bool_t what = IsEventSpecieSet(ies) ;
+        if ( what )
+          ShowStatus(DETECTORINDEX_t(index), ALITASK_t(tsk), AliRecoParam::Convert(ies)) ;
+      }
+    }
+       }
 }
 
 //_______________________________________________________________
-void AliQA::ShowStatus(DETECTORINDEX_t det) const
+void AliQA::ShowStatus(DETECTORINDEX_t det, ALITASK_t tsk, AliRecoParam::EventSpecie_t es) const
 {
-  // Prints the full QA status of a given detector
-  CheckRange(det) ;
-  ULong_t status = GetStatus(det) ;
-  ULong_t rawStatus = status & 0x0000f ;
-  ULong_t simStatus = status & 0x000f0 ;
-  ULong_t recStatus = status & 0x00f00 ;
-  ULong_t esdStatus = status & 0x0f000 ;
-  ULong_t anaStatus = status & 0xf0000 ;
+       // Prints the full QA status of a given detector
+       CheckRange(det) ;
+       CheckRange(es) ;
+       ULong_t status = GetStatus(det, es) ;
+       ULong_t tskStatus[kNTASK] ; 
+       tskStatus[kRAW] = status & 0x0000f ;
+       tskStatus[kSIM] = status & 0x000f0 ;
+       tskStatus[kREC] = status & 0x00f00 ;
+       tskStatus[kESD] = status & 0x0f000 ;
+       tskStatus[kANA] = status & 0xf0000 ;
+
+       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), 
+                                tskStatus[kRAW], tskStatus[kSIM], tskStatus[kREC], tskStatus[kESD], tskStatus[kANA] )) ;
+       if (tsk == kNULLTASK) {
+               for (Int_t itsk = kRAW ; itsk < kNTASK ; itsk++) {
+                       ShowASCIIStatus(es, det, ALITASK_t(itsk), tskStatus[itsk]) ; 
+               } 
+       } else {
+                       ShowASCIIStatus(es, det, tsk, tskStatus[tsk]) ; 
+       }
+}
 
-  AliInfo(Form("QA Status for %8s raw =0x%x, sim=0x%x, rec=0x%x, esd=0x%x, ana=0x%x\n", GetDetName(det).Data(), rawStatus, simStatus, recStatus, esdStatus, anaStatus )) ;
+//_______________________________________________________________
+void AliQA::ShowASCIIStatus(AliRecoParam::EventSpecie_t es, DETECTORINDEX_t det, ALITASK_t tsk, const ULong_t status) const 
+{
+       // print the QA status in human readable format
+       TString text; 
+       for (Int_t bit = kINFO ; bit < kNBIT ; bit++) {
+               if (IsSet(det, tsk, es, QABIT_t(bit))) {
+                       text = GetBitName(QABIT_t(bit)) ; 
+                       text += " " ; 
+               }
+       }
+       if (! text.IsNull())
+               printf("           %8s %8s %4s 0x%4lx, Problem signalled: %8s \n", AliRecoParam::GetEventSpecieName(es), GetDetName(det).Data(), GetAliTaskName(tsk), status, text.Data()) ; 
+}
+
+//_______________________________________________________________
+void AliQA::UnSet(QABIT_t bit, Int_t ies)
+{
+       // UnSet the status bit of the current detector in the current module
+               UnSet(bit, AliRecoParam::Convert(ies)) ;
 }
 
 //_______________________________________________________________
-void AliQA::UnSet(QABIT_t bit)
+void AliQA::UnSet(QABIT_t bit, AliRecoParam::EventSpecie_t es)
 {
        // UnSet the status bit of the current detector in the current module
        
-       UnSetStatusBit(fDet, fTask, bit) ;
+       UnSetStatusBit(fDet, fTask, es, bit) ;
 }
 
 //_______________________________________________________________
-void AliQA::UnSetStatusBit(DETECTORINDEX_t det, ALITASK_t tsk, QABIT_t bit)
+void AliQA::UnSetStatusBit(DETECTORINDEX_t det, ALITASK_t tsk, AliRecoParam::EventSpecie_t es, QABIT_t bit)
 {
        // UnSet the status bit for a given detector and a given task
        
        CheckRange(det) ;
        CheckRange(tsk) ;
        CheckRange(bit) ;
+       CheckRange(es) ;
        
        ULong_t offset = Offset(tsk) ;
-       ULong_t status = GetStatus(det) ;
+       ULong_t status = GetStatus(det, es) ;
        offset+= bit ;
        status = status & 0 << offset ;
-       SetStatus(det, status) ;
+       SetStatus(det, es, status) ;
 }