]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - STEER/AliQAManager.cxx
1.The QA data created on demand according to the event species at filling time. 2...
[u/mrichter/AliRoot.git] / STEER / AliQAManager.cxx
index d813b525ba9a70e92b083eecdb402f39a07f2699..c1d73b925934fec3f952b2f5ecabeda9f7acd9cf 100644 (file)
@@ -28,6 +28,7 @@
 //                                                                           //
 ///////////////////////////////////////////////////////////////////////////////
 
+#include <TCanvas.h>
 #include <TKey.h>
 #include <TFile.h>
 #include <TFileMerger.h>
@@ -85,7 +86,9 @@ AliQAManager::AliQAManager() :
   fRawReader(NULL), 
   fRawReaderDelete(kTRUE), 
   fRunLoader(NULL), 
-  fTasks("") 
+  fTasks(""),  
+  fEventSpecie(AliRecoParam::kDefault), 
+  fPrintImage(kTRUE) 
 {
        // default ctor
        fMaxEvents = fNumberOfEvents ; 
@@ -118,7 +121,9 @@ AliQAManager::AliQAManager(const Char_t * mode, const Char_t* gAliceFilename) :
        fRawReader(NULL), 
        fRawReaderDelete(kTRUE), 
        fRunLoader(NULL), 
-  fTasks("") 
+  fTasks(""), 
+  fEventSpecie(AliRecoParam::kDefault), 
+  fPrintImage(kTRUE) 
 {
        // default ctor
        fMaxEvents = fNumberOfEvents ; 
@@ -130,6 +135,10 @@ AliQAManager::AliQAManager(const Char_t * mode, const Char_t* gAliceFilename) :
     }
   }
   SetWriteExpert() ; 
+  if (fMode.Contains("sim")) 
+               fMode.ReplaceAll("s", "S") ; 
+  else if (fMode.Contains("rec")) 
+    fMode.ReplaceAll("r", "R") ; 
 }
 
 //_____________________________________________________________________________
@@ -151,7 +160,10 @@ AliQAManager::AliQAManager(const AliQAManager & qas) :
        fRawReader(NULL), 
        fRawReaderDelete(kTRUE), 
        fRunLoader(NULL), 
-  fTasks(qas.fTasks)
+  fTasks(qas.fTasks), 
+  fEventSpecie(qas.fEventSpecie), 
+  fPrintImage(qas.fPrintImage) 
+
 {
        // cpy ctor
        for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
@@ -189,6 +201,9 @@ AliQAManager::~AliQAManager()
                delete fRawReader ;
                fRawReader = NULL ;
        }
+  TCanvas fakeCanvas ; 
+  if (fPrintImage) 
+    fakeCanvas.Print(Form("%s%s%d.%s]", AliQAv1::GetImageFileName(), GetMode(), fRunNumber, AliQAv1::GetImageFileFormat())); 
 }
 
 //_____________________________________________________________________________
@@ -237,8 +252,8 @@ Bool_t AliQAManager::DoIt(const AliQAv1::TASKINDEX_t taskIndex)
                 AliWarning(Form(" Hit Tree not found for  %s", AliQAv1::GetDetName(iDet))) ; 
                 break ; 
               } 
+              qadm->Exec(taskIndex, data) ;
             } 
-            qadm->Exec(taskIndex, data) ;
                                                break ;
                                                case AliQAv1::kSDIGITS :
             if( loader ) {      
@@ -248,19 +263,30 @@ Bool_t AliQAManager::DoIt(const AliQAv1::TASKINDEX_t taskIndex)
                 AliWarning(Form(" SDigit Tree not found for  %s", AliQAv1::GetDetName(iDet))) ; 
                 break ; 
               } 
+              qadm->Exec(taskIndex, data) ; 
             }
-            qadm->Exec(taskIndex, data) ; 
                                                break; 
                                                case AliQAv1::kDIGITS :
-            if( loader ) {      
-              loader->LoadDigits() ; 
-              data = loader->TreeD() ; 
-              if ( ! data ) {
-                AliWarning(Form(" Digit Tree not found for  %s", AliQAv1::GetDetName(iDet))) ; 
-                break ; 
-              } 
-            }
-            qadm->Exec(taskIndex, data) ;
+              if( loader ) {      
+                loader->LoadDigits() ; 
+                data = loader->TreeD() ; 
+                if ( ! data ) {
+                  AliWarning(Form(" Digit Tree not found for  %s", AliQAv1::GetDetName(iDet))) ; 
+                  break ; 
+                } 
+                qadm->Exec(taskIndex, data) ;
+              }
+            break;
+            case AliQAv1::kDIGITSR :
+              if( loader ) {      
+                loader->LoadDigits() ; 
+                data = loader->TreeD() ; 
+                if ( ! data ) {
+                  AliWarning(Form(" Digit Tree not found for  %s", AliQAv1::GetDetName(iDet))) ; 
+                  break ; 
+                } 
+                qadm->Exec(taskIndex, data) ;
+              }
                                                break; 
                                                case AliQAv1::kRECPOINTS :
             if( loader ) {      
@@ -270,8 +296,8 @@ Bool_t AliQAManager::DoIt(const AliQAv1::TASKINDEX_t taskIndex)
                 AliWarning(Form("RecPoints not found for %s", AliQAv1::GetDetName(iDet))) ; 
                 break ; 
               } 
+              qadm->Exec(taskIndex, data) ; 
             }
-            qadm->Exec(taskIndex, data) ; 
             break; 
                                                case AliQAv1::kTRACKSEGMENTS :
                                                break; 
@@ -311,7 +337,7 @@ Bool_t AliQAManager::Finish(const AliQAv1::TASKINDEX_t taskIndex)
 }
 
 //_____________________________________________________________________________
-TObjArray * AliQAManager::GetFromOCDB(AliQAv1::DETECTORINDEX_t det, AliQAv1::TASKINDEX_t task, const char * year) const 
+TObjArray * AliQAManager::GetFromOCDB(AliQAv1::DETECTORINDEX_t det, AliQAv1::TASKINDEX_t task, const Char_t * year) const 
 {
        // Retrieve the list of QA data for a given detector and a given task 
        TObjArray * rv = NULL ;
@@ -329,12 +355,23 @@ TObjArray * AliQAManager::GetFromOCDB(AliQAv1::DETECTORINDEX_t det, AliQAv1::TAS
        TString detOCDBDir(Form("%s/%s/%s", AliQAv1::GetQAName(), AliQAv1::GetDetName((Int_t)det), AliQAv1::GetRefOCDBDirName())) ; 
      AliDebug(AliQAv1::GetQADebugLevel(), Form("Retrieving reference data from %s/%s for %s", AliQAv1::GetQARefStorage(), detOCDBDir.Data(), AliQAv1::GetTaskName(task).Data())) ; 
        AliCDBEntry* entry = QAManager()->Get(detOCDBDir.Data(), 0) ; //FIXME 0 --> Run Number
-       TList * listDetQAD = dynamic_cast<TList *>(entry->GetObject()) ;
+       TList * listDetQAD = static_cast<TList *>(entry->GetObject()) ;
        if ( listDetQAD ) 
-               rv = dynamic_cast<TObjArray *>(listDetQAD->FindObject(AliQAv1::GetTaskName(task))) ; 
+               rv = static_cast<TObjArray *>(listDetQAD->FindObject(AliQAv1::GetTaskName(task))) ; 
        return rv ; 
 }
 
+//_____________________________________________________________________________
+TCanvas ** AliQAManager::GetImage(Char_t * detName)
+{
+  // retrieves QA Image for the given detector
+  TCanvas ** rv = NULL ; 
+  Int_t detIndex = AliQAv1::GetDetIndex(detName) ; 
+  AliQADataMaker * qadm = GetQADataMaker(detIndex) ; 
+  rv = qadm->GetImage() ;
+  return rv ; 
+}
+
 //_____________________________________________________________________________
 AliLoader * AliQAManager::GetLoader(Int_t iDet)
 {
@@ -378,13 +415,13 @@ AliLoader * AliQAManager::GetLoader(Int_t iDet)
 AliQAv1 * AliQAManager::GetQA(UInt_t run, UInt_t evt) 
 {
 // retrieves the QA object stored in a file named "Run{run}.Event{evt}_1.ESD.tag.root"  
-  char * fileName = Form("Run%d.Event%d_1.ESD.tag.root", run, evt) ; 
+  Char_t * fileName = Form("Run%d.Event%d_1.ESD.tag.root", run, evt) ; 
   TFile * tagFile = TFile::Open(fileName) ;
   if ( !tagFile ) {
     AliError(Form("File %s not found", fileName)) ;
     return NULL ; 
   }
-  TTree * tagTree = dynamic_cast<TTree *>(tagFile->Get("T")) ; 
+  TTree * tagTree = static_cast<TTree *>(tagFile->Get("T")) ; 
   if ( !tagTree ) {
     AliError(Form("Tree T not found in %s", fileName)) ; 
     tagFile->Close() ; 
@@ -403,65 +440,69 @@ AliQADataMaker * AliQAManager::GetQADataMaker(const Int_t iDet)
 {
        // get the quality assurance data maker for a detector
        
-       if (fQADataMaker[iDet]) {
-    fQADataMaker[iDet]->SetEventSpecie(fQADataMaker[iDet]->GetRecoParam()->GetEventSpecie()) ; 
-               return fQADataMaker[iDet] ;
-  }
-       
-       AliQADataMaker * qadm = NULL ;
-       
-       if (iDet == AliQAv1::kGLOBAL) { //Global QA
+  AliQADataMaker * qadm =  fQADataMaker[iDet] ; 
+  
+       if (qadm) {
+    qadm->SetEventSpecie(fEventSpecie) ;  
+    if ( qadm->GetRecoParam() ) 
+      if ( AliRecoParam::Convert(qadm->GetRecoParam()->GetEventSpecie()) != AliRecoParam::kDefault)
+        qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ; 
+
+  } else if (iDet == AliQAv1::kGLOBAL) { //Global QA
+
                qadm = new AliGlobalQADataMaker();
                qadm->SetName(AliQAv1::GetDetName(iDet));
                qadm->SetUniqueID(iDet);
                fQADataMaker[iDet] = qadm;
-    qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ; 
-               return qadm;
-       }
+    qadm->SetEventSpecie(fEventSpecie) ;  
+    if ( qadm->GetRecoParam() ) 
+      if ( AliRecoParam::Convert(qadm->GetRecoParam()->GetEventSpecie()) != AliRecoParam::kDefault)  
+        qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ; 
 
-       if (iDet == AliQAv1::kCORR) { //the data maker for correlations among detectors
+       }       else if (iDet == AliQAv1::kCORR && strcmp(GetMode(), "Rec") == 0 ) { //the data maker for correlations among detectors
     qadm = new AliCorrQADataMakerRec(fQADataMaker) ; 
                qadm->SetName(AliQAv1::GetDetName(iDet));
                qadm->SetUniqueID(iDet);
                fQADataMaker[iDet] = qadm;
-    qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ; 
-               return qadm;
-  }
-
-       // load the QA data maker object
-       TPluginManager* pluginManager = gROOT->GetPluginManager() ;
-       TString detName = AliQAv1::GetDetName(iDet) ;
-       TString tmp(fMode) ; 
-       if (tmp.Contains("sim")) 
-               tmp.ReplaceAll("s", "S") ; 
-       else if (tmp.Contains("rec")) 
-               tmp.ReplaceAll("r", "R") ; 
-
-       TString qadmName = "Ali" + detName + "QADataMaker" + tmp ;
-
-       // first check if a plugin is defined for the quality assurance data maker
-       TPluginHandler* pluginHandler = pluginManager->FindHandler("AliQADataMaker", detName) ;
-       // if not, add a plugin for it
-       if (!pluginHandler) {
-               AliDebug(AliQAv1::GetQADebugLevel(), Form("defining plugin for %s", qadmName.Data())) ;
-               TString libs = gSystem->GetLibraries() ;
-               if (libs.Contains("lib" + detName + fMode + ".so") || (gSystem->Load("lib" + detName + fMode + ".so") >= 0)) {
-                       pluginManager->AddHandler("AliQADataMaker", detName, qadmName, detName + "qadm", qadmName + "()") ;
-               } else {
-                       pluginManager->AddHandler("AliQADataMaker", detName, qadmName, detName, qadmName + "()") ;
-               }
-               pluginHandler = pluginManager->FindHandler("AliQADataMaker", detName) ;
-       }
-       if (pluginHandler && (pluginHandler->LoadPlugin() == 0)) {
-               qadm = (AliQADataMaker *) pluginHandler->ExecPlugin(0) ;
-       }
-       if (qadm) {
-               qadm->SetName(AliQAv1::GetDetName(iDet));
-               qadm->SetUniqueID(iDet);
-               fQADataMaker[iDet] = qadm ;
-    qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ; 
-       }
+    qadm->SetEventSpecie(fEventSpecie) ;  
+    if ( qadm->GetRecoParam() ) 
+      if ( AliRecoParam::Convert(qadm->GetRecoParam()->GetEventSpecie()) != AliRecoParam::kDefault)  
+        qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ; 
 
+  } else {
+    
+    // load the QA data maker object
+    TPluginManager* pluginManager = gROOT->GetPluginManager() ;
+    TString detName = AliQAv1::GetDetName(iDet) ;
+    TString qadmName = "Ali" + detName + "QADataMaker" + GetMode() ;
+    
+    // first check if a plugin is defined for the quality assurance data maker
+    TPluginHandler* pluginHandler = pluginManager->FindHandler("AliQADataMaker", detName) ;
+    // if not, add a plugin for it
+    if (!pluginHandler) {
+      AliDebug(AliQAv1::GetQADebugLevel(), Form("defining plugin for %s", qadmName.Data())) ;
+      TString libs = gSystem->GetLibraries() ;
+      if (libs.Contains("lib" + detName + GetMode() + ".so") || (gSystem->Load("lib" + detName + GetMode() + ".so") >= 0)) {
+        pluginManager->AddHandler("AliQADataMaker", detName, qadmName, detName + "qadm", qadmName + "()") ;
+      } else {
+        pluginManager->AddHandler("AliQADataMaker", detName, qadmName, detName, qadmName + "()") ;
+      }
+      pluginHandler = pluginManager->FindHandler("AliQADataMaker", detName) ;
+    }
+    if (pluginHandler && (pluginHandler->LoadPlugin() == 0)) {
+      qadm = (AliQADataMaker *) pluginHandler->ExecPlugin(0) ;
+    }
+    if (qadm) {
+      qadm->SetName(AliQAv1::GetDetName(iDet));
+      qadm->SetUniqueID(iDet);
+      fQADataMaker[iDet] = qadm ;
+      qadm->SetEventSpecie(fEventSpecie) ;  
+      if ( qadm->GetRecoParam() ) 
+        if ( AliRecoParam::Convert(qadm->GetRecoParam()->GetEventSpecie()) != AliRecoParam::kDefault)  
+          qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ; 
+    }
+  }
   return qadm ;
 }
 
@@ -470,6 +511,10 @@ void  AliQAManager::EndOfCycle(TObjArray * detArray)
 {
        // End of cycle QADataMakers 
        
+  if (fPrintImage) {
+    TCanvas fakeCanvas ; 
+    fakeCanvas.Print(Form("%s%s%d.%s[", AliQAv1::GetImageFileName(), GetMode(), fRunNumber, AliQAv1::GetImageFileFormat())) ; 
+  }
        for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
                if (IsSelected(AliQAv1::GetDetName(iDet))) {
                        AliQADataMaker * qadm = GetQADataMaker(iDet) ;
@@ -481,6 +526,8 @@ void  AliQAManager::EndOfCycle(TObjArray * detArray)
                                if (!det || !det->IsActive())  
                                        continue ;
                        }
+      qadm->SetPrintImage(fPrintImage) ;
+      
                        for (UInt_t taskIndex = 0; taskIndex < AliQAv1::kNTASKINDEX; taskIndex++) {
                                if ( fTasks.Contains(Form("%d", taskIndex)) ) 
                                        qadm->EndOfCycle(AliQAv1::GetTaskIndex(AliQAv1::GetTaskName(taskIndex))) ;
@@ -495,7 +542,11 @@ void  AliQAManager::EndOfCycle(TString detectors)
 {
        // End of cycle QADataMakers 
        
-       for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
+  if (fPrintImage) {
+    TCanvas fakeCanvas ; 
+    fakeCanvas.Print(Form("%s%s%d.%s[", AliQAv1::GetImageFileName(), GetMode(), fRunNumber, AliQAv1::GetImageFileFormat())) ; 
+  }
+  for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
                if (IsSelected(AliQAv1::GetDetName(iDet))) {
                        AliQADataMaker * qadm = GetQADataMaker(iDet) ;
                        if (!qadm) 
@@ -503,6 +554,7 @@ void  AliQAManager::EndOfCycle(TString detectors)
                        // skip non active detectors
       if (!detectors.Contains(AliQAv1::GetDetName(iDet))) 
         continue ;
+      qadm->SetPrintImage(fPrintImage) ;
                for (UInt_t taskIndex = 0; taskIndex < AliQAv1::kNTASKINDEX; taskIndex++) {
                                if ( fTasks.Contains(Form("%d", taskIndex)) ) 
                                        qadm->EndOfCycle(AliQAv1::GetTaskIndex(AliQAv1::GetTaskName(taskIndex))) ;
@@ -526,7 +578,7 @@ void AliQAManager::Increment()
 }
   
 //_____________________________________________________________________________
-Bool_t AliQAManager::InitQA(const AliQAv1::TASKINDEX_t taskIndex, const  char * input )
+Bool_t AliQAManager::InitQA(const AliQAv1::TASKINDEX_t taskIndex, const  Char_t * input )
 {
        // Initialize the event source and QA data makers
        
@@ -550,7 +602,7 @@ Bool_t AliQAManager::InitQA(const AliQAv1::TASKINDEX_t taskIndex, const  char *
                        fTasks = AliQAv1::GetTaskName(AliQAv1::kESDS) ; 
       if (!gSystem->AccessPathName("AliESDs.root")) { // AliESDs.root exists
         TFile * esdFile = TFile::Open("AliESDs.root") ;
-        fESDTree = dynamic_cast<TTree *> (esdFile->Get("esdTree")) ; 
+        fESDTree = static_cast<TTree *> (esdFile->Get("esdTree")) ; 
         if ( !fESDTree ) {
           AliError("esdTree not found") ; 
           return kFALSE ; 
@@ -591,6 +643,10 @@ Bool_t AliQAManager::InitQA(const AliQAv1::TASKINDEX_t taskIndex, const  char *
                AliGeomManager::LoadGeometry() ; 
        
        InitQADataMaker(fRunNumber, detArray) ; //, fCycleSame, kTRUE, detArray) ; 
+  if (fPrintImage) {
+    TCanvas fakeCanvas ; 
+    fakeCanvas.Print(Form("%s%s%d.%s[", AliQAv1::GetImageFileName(), GetMode(), fRunNumber, AliQAv1::GetImageFileFormat())) ;    
+  }    
        return kTRUE ; 
 }
 
@@ -598,6 +654,7 @@ Bool_t AliQAManager::InitQA(const AliQAv1::TASKINDEX_t taskIndex, const  char *
 void  AliQAManager::InitQADataMaker(UInt_t run, TObjArray * detArray) 
 {
        // Initializes The QADataMaker for all active detectors and for all active tasks 
+  fRunNumber = run ; 
        for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
                if (IsSelected(AliQAv1::GetDetName(iDet))) {
                        AliQADataMaker * qadm = GetQADataMaker(iDet) ;
@@ -678,7 +735,7 @@ Bool_t AliQAManager::InitRunLoader()
 }
 
 //_____________________________________________________________________________
-Bool_t AliQAManager::IsSelected(const char * det) 
+Bool_t AliQAManager::IsSelected(const Char_t * det) 
 {
   // check whether detName is contained in detectors
        // if yes, it is removed from detectors
@@ -704,19 +761,20 @@ Bool_t AliQAManager::IsSelected(const char * det)
 }
 
 //_____________________________________________________________________________
-Bool_t AliQAManager::Merge(Int_t runNumber) const
+Bool_t AliQAManager::Merge(Int_t runNumber, const char *fileName) const
 {
        // Merge data from all detectors from a given run in one single file 
-       // Merge the QA results from all the data chunks in one run 
+       // Merge the QA results from all the data chunks in one run
+  // The 'fileName' is name of the output file with merged QA data  
  if ( runNumber == -1)
    runNumber = fRunNumber ; 
- Bool_t rv = MergeData(runNumber) ; 
+ Bool_t rv = MergeData(runNumber,fileName) ; 
  //rv *= MergeResults(runNumber) ; // not needed for the time being
  return rv ; 
 }
        
 //______________________________________________________________________
-Bool_t AliQAManager::MergeXML(const char * collectionFile, const char * subFile, const char * outFile) 
+Bool_t AliQAManager::MergeXML(const Char_t * collectionFile, const Char_t * subFile, const Char_t * outFile) 
 {
   // merges files listed in a xml collection 
   // usage Merge(collection, outputFile))
@@ -735,15 +793,13 @@ Bool_t AliQAManager::MergeXML(const char * collectionFile, const char * subFile,
    return kFALSE ; 
  
   // Open the file collection 
-  printf("*** Create Collection       ***\n");
-  printf("***  Wk-Dir = |%s|             \n",gSystem->WorkingDirectory());
-  printf("***  Coll   = |%s|             \n",collectionFile);                  
+  AliInfoClass(Form("*** Create Collection       ***\n***  Wk-Dir = |%s|             \n***  Coll   = |%s|             \n",gSystem->WorkingDirectory(), collectionFile));               
   
   TGridCollection * collection = (TGridCollection*)gROOT->ProcessLine(Form("TAlienCollection::Open(\"%s\")",collectionFile));
   TGridResult* result = collection->GetGridResult("", 0, 0);
   
   Int_t index = 0  ;
-  const char * turl ;
+  const Char_t * turl ;
   TFileMerger merger(kFALSE) ; 
   if (!outFile) {
     TString tempo(collectionFile) ; 
@@ -756,7 +812,7 @@ Bool_t AliQAManager::MergeXML(const char * collectionFile, const char * subFile,
   merger.OutputFile(outFile) ; 
   
   while ( (turl = result->GetKey(index, "turl")) ) {
-    char * file ;
+    Char_t * file ;
     if ( subFile )
       file = Form("%s#%s", turl, subFile) ; 
     else 
@@ -799,7 +855,7 @@ void AliQAManager::MergeCustom() const
       loRun = cuRun ; 
     if (cuRun > hiRun)
       hiRun = cuRun ; 
-    while ( (srun = dynamic_cast<TObjString *> (nextRun())) ) {
+    while ( (srun = static_cast<TObjString *> (nextRun())) ) {
       if ( cuRun == (srun->String()).Atoi() ) {
         runExist = kTRUE ; 
         break ; 
@@ -827,15 +883,15 @@ void AliQAManager::MergeCustom() const
       }
     }
   }
-  while ( (srun = dynamic_cast<TObjString *> (nextRun())) ) {
+  while ( (srun = static_cast<TObjString *> (nextRun())) ) {
     runNumber = (srun->String()).Atoi() ; 
     hisRun->Fill(runNumber) ; 
     AliDebug(AliQAv1::GetQADebugLevel(), Form("Merging run number %d", runNumber)) ; 
     // search all QA files for runNumber in the current directory
-    char * fileList[AliQAv1::kNDET] ;
+    Char_t * fileList[AliQAv1::kNDET] ;
     index = 0 ; 
     for (Int_t iDet = 0; iDet < AliQAv1::kNDET ; iDet++) {
-      char * file = gSystem->Which(gSystem->WorkingDirectory(), Form("%s.%s.%d.root", AliQAv1::GetDetName(iDet), AliQAv1::GetQADataFileName(), runNumber)); 
+      Char_t * file = gSystem->Which(gSystem->WorkingDirectory(), Form("%s.%s.%d.root", AliQAv1::GetDetName(iDet), AliQAv1::GetQADataFileName(), runNumber)); 
       if (file) 
         fileList[index++] = file ;
     }
@@ -878,11 +934,11 @@ void AliQAManager::MergeCustom() const
                   TList * listOfData = directoryEventSpecie->GetListOfKeys() ; 
                   TIter nextData(listOfData) ; 
                   TKey * key ; 
-                  while ( (key = dynamic_cast<TKey *>(nextData())) ) {
+                  while ( (key = static_cast<TKey *>(nextData())) ) {
                     TString className(key->GetClassName()) ; 
                     if (  className.Contains("TH") || className.Contains("TProfile") ) {
-                      TH1 * histIn = dynamic_cast<TH1*> (key->ReadObj()) ; 
-                      TH1 * histOu = dynamic_cast<TH1*> (mergedDirectory->FindObjectAny(histIn->GetName())) ; 
+                      TH1 * histIn = static_cast<TH1*> (key->ReadObj()) ; 
+                      TH1 * histOu = static_cast<TH1*> (mergedDirectory->FindObjectAny(histIn->GetName())) ; 
                       AliDebug(AliQAv1::GetQADebugLevel(), Form("%s %x %x\n", key->GetName(), histIn, histOu)) ; 
                       mergedDirectory->cd() ; 
                       if ( ! histOu ) {
@@ -899,11 +955,11 @@ void AliQAManager::MergeCustom() const
                       TList * listOfExpertData = dirExpert->GetListOfKeys() ; 
                       TIter nextExpertData(listOfExpertData) ; 
                       TKey * keykey ; 
-                      while ( (keykey = dynamic_cast<TKey *>(nextExpertData())) ) {
+                      while ( (keykey = static_cast<TKey *>(nextExpertData())) ) {
                         TString classNameExpert(keykey->GetClassName()) ; 
                         if (classNameExpert.Contains("TH")) {
-                          TH1 * histInExpert = dynamic_cast<TH1*> (keykey->ReadObj()) ; 
-                          TH1 * histOuExpert = dynamic_cast<TH1*> (mergedDirectory->FindObjectAny(histInExpert->GetName())) ; 
+                          TH1 * histInExpert = static_cast<TH1*> (keykey->ReadObj()) ; 
+                          TH1 * histOuExpert = static_cast<TH1*> (mergedDirectory->FindObjectAny(histInExpert->GetName())) ; 
                           mergedDirectoryExpert->cd() ; 
                           if ( ! histOuExpert ) {
                             histInExpert->Write() ; 
@@ -935,15 +991,16 @@ void AliQAManager::MergeCustom() const
 }
 
 //_____________________________________________________________________________
-Bool_t AliQAManager::MergeData(const Int_t runNumber) const
+Bool_t AliQAManager::MergeData(const Int_t runNumber, const char *fileName) const
 {
        // Merge QA data from all detectors for a given run in one single file 
   
-  TFileMerger merger ; 
-  TString outFileName = Form("Merged.%s.Data.root",AliQAv1::GetQADataFileName()) ;
+  TFileMerger merger(kFALSE) ; 
+  TString outFileName = fileName;
+  if (outFileName.IsNull()) outFileName.Form("Merged.%s.Data.root",AliQAv1::GetQADataFileName());
   merger.OutputFile(outFileName.Data()) ; 
   for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
-    char * file = gSystem->Which(gSystem->WorkingDirectory(), Form("%s.%s.%d.root", AliQAv1::GetDetName(iDet), AliQAv1::GetQADataFileName(), runNumber)); 
+    Char_t * file = gSystem->Which(gSystem->WorkingDirectory(), Form("%s.%s.%d.root", AliQAv1::GetDetName(iDet), AliQAv1::GetQADataFileName(), runNumber)); 
     if (file) 
       merger.AddFile(file) ; 
   }
@@ -1023,6 +1080,8 @@ AliQAManager * AliQAManager::QAManager(const Char_t * mode, TMap *entryCache, In
   // returns AliQAManager instance (singleton)
   
        if (!fgQAInstance) {
+    if ( (strcmp(mode, "sim") != 0) && (strcmp(mode, "rec") != 0) )
+      AliFatalClass("You must specify sim or rec") ; 
     fgQAInstance = new AliQAManager(mode) ;  
     if (!entryCache)
                  fgQAInstance->Init();
@@ -1033,7 +1092,7 @@ AliQAManager * AliQAManager::QAManager(const Char_t * mode, TMap *entryCache, In
 }
 
 //_____________________________________________________________________________
-TString AliQAManager::Run(const char * detectors, AliRawReader * rawReader, const Bool_t sameCycle) 
+TString AliQAManager::Run(const Char_t * detectors, AliRawReader * rawReader, const Bool_t sameCycle) 
 {
        //Runs all the QA data Maker for Raws only
        
@@ -1060,7 +1119,7 @@ TString AliQAManager::Run(const char * detectors, AliRawReader * rawReader, cons
 }
 
 //_____________________________________________________________________________
-TString AliQAManager::Run(const char * detectors, const char * fileName, const Bool_t sameCycle) 
+TString AliQAManager::Run(const Char_t * detectors, const Char_t * fileName, const Bool_t sameCycle) 
 {
        //Runs all the QA data Maker for Raws only
 
@@ -1094,7 +1153,7 @@ TString AliQAManager::Run(const char * detectors, const char * fileName, const B
 }
 
 //_____________________________________________________________________________
-TString AliQAManager::Run(const char * detectors, const AliQAv1::TASKINDEX_t taskIndex, Bool_t const sameCycle, const  char * fileName ) 
+TString AliQAManager::Run(const Char_t * detectors, const AliQAv1::TASKINDEX_t taskIndex, Bool_t const sameCycle, const  Char_t * fileName ) 
 {
        // Runs all the QA data Maker for every detector
        
@@ -1118,9 +1177,7 @@ TString AliQAManager::Run(const char * detectors, const AliQAv1::TASKINDEX_t tas
                        }
                }
        }
-       
-
-       if ( taskIndex == AliQAv1::kNULLTASKINDEX) { 
+  if ( taskIndex == AliQAv1::kNULLTASKINDEX) { 
                for (UInt_t task = 0; task < AliQAv1::kNTASKINDEX; task++) {
                        if ( fTasks.Contains(Form("%d", task)) ) {
         if (!fCycleSame)
@@ -1135,9 +1192,7 @@ TString AliQAManager::Run(const char * detectors, const AliQAv1::TASKINDEX_t tas
         return kFALSE ; 
       DoIt(taskIndex) ; 
   }            
-       
        return fDetectorsW ;
-
 }
 
 //_____________________________________________________________________________
@@ -1158,7 +1213,10 @@ void AliQAManager::RunOneEvent(AliRawReader * rawReader)
         qadm->EndOfCycle() ;
       }
       AliCodeTimerStart(Form("running RAW quality assurance data maker for %s", AliQAv1::GetDetName(iDet))); 
-      qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ; 
+      qadm->SetEventSpecie(fEventSpecie) ;  
+      if ( qadm->GetRecoParam() ) 
+        if ( AliRecoParam::Convert(qadm->GetRecoParam()->GetEventSpecie()) != AliRecoParam::kDefault)  
+          qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ; 
                        qadm->Exec(AliQAv1::kRAWS, rawReader) ;
       AliCodeTimerStop(Form("running RAW quality assurance data maker for %s", AliQAv1::GetDetName(iDet)));
                }
@@ -1178,6 +1236,10 @@ void AliQAManager::RunOneEvent(AliESDEvent *& esd)
       AliQADataMaker *qadm = GetQADataMaker(iDet);  
       if (!qadm) 
         continue;
+      qadm->SetEventSpecie(fEventSpecie) ;  
+      if ( qadm->GetRecoParam() ) 
+        if ( AliRecoParam::Convert(qadm->GetRecoParam()->GetEventSpecie()) != AliRecoParam::kDefault)  
+          qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ; 
       if ( qadm->IsCycleDone() ) {
         qadm->EndOfCycle() ;
       }
@@ -1192,24 +1254,34 @@ void AliQAManager::RunOneEvent(AliESDEvent *& esd)
 void AliQAManager::RunOneEventInOneDetector(Int_t det, TTree * tree) 
 {
        // Runs all the QA data Maker for ESDs only and on one event only (event loop done by calling method)
+  
+  TString test(tree->GetName()) ; 
        AliCodeTimerAuto("") ;
   if (fTasks.Contains(Form("%d", AliQAv1::kRECPOINTS))) {
     if (IsSelected(AliQAv1::GetDetName(det))) {
       AliQADataMaker *qadm = GetQADataMaker(det);  
       if (qadm) { 
+        qadm->SetEventSpecie(fEventSpecie) ;  
+        if ( qadm->GetRecoParam() ) 
+          if ( AliRecoParam::Convert(qadm->GetRecoParam()->GetEventSpecie()) != AliRecoParam::kDefault)  
+            qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ;         
         if ( qadm->IsCycleDone() ) {
           qadm->EndOfCycle() ;
         }
         AliCodeTimerStart(Form("running RecPoints quality assurance data maker for %s", AliQAv1::GetDetName(det)));
-        qadm->Exec(AliQAv1::kRECPOINTS, tree) ;
-        AliCodeTimerStop(Form("running RecPoints quality assurance data maker for %s", AliQAv1::GetDetName(det)));
+        if (test.Contains("TreeD")) {
+          qadm->Exec(AliQAv1::kDIGITSR, tree) ;
+        } else  if (test.Contains("TreeR")) {
+          qadm->Exec(AliQAv1::kRECPOINTS, tree) ;
+          AliCodeTimerStop(Form("running RecPoints quality assurance data maker for %s", AliQAv1::GetDetName(det)));
+        }
       }
     }
   }
 }
 
 //_____________________________________________________________________________
-Bool_t AliQAManager::Save2OCDB(const Int_t runNumber, AliRecoParam::EventSpecie_t es, const char * year, const char * detectors) const
+Bool_t AliQAManager::Save2OCDB(const Int_t runNumber, AliRecoParam::EventSpecie_t es, const Char_t * year, const Char_t * detectors) const
 {
        // take the locasl QA data merge into a single file and save in OCDB 
        Bool_t rv = kTRUE ; 
@@ -1245,7 +1317,7 @@ Bool_t AliQAManager::Save2OCDB(const Int_t runNumber, AliRecoParam::EventSpecie_
 }
 
 //_____________________________________________________________________________
-Bool_t AliQAManager::SaveIt2OCDB(const Int_t runNumber, TFile * inputFile, const char * year, AliRecoParam::EventSpecie_t es) const
+Bool_t AliQAManager::SaveIt2OCDB(const Int_t runNumber, TFile * inputFile, const Char_t * year, AliRecoParam::EventSpecie_t es) const
 {
        // reads the TH1 from file and adds it to appropriate list before saving to OCDB
        Bool_t rv = kTRUE ;
@@ -1282,7 +1354,7 @@ Bool_t AliQAManager::SaveIt2OCDB(const Int_t runNumber, TFile * inputFile, const
                        TList * taskList = detDir->GetListOfKeys() ; 
                        TIter nextTask(taskList) ; 
                        TKey * taskKey ; 
-                       while ( (taskKey = dynamic_cast<TKey*>(nextTask())) ) {
+                       while ( (taskKey = static_cast<TKey*>(nextTask())) ) {
                                TDirectory * taskDir = detDir->GetDirectory(taskKey->GetName()) ; 
         TDirectory * esDir   = taskDir->GetDirectory(AliRecoParam::GetEventSpecieName(es)) ; 
                                AliDebug(AliQAv1::GetQADebugLevel(), Form("Saving %s", esDir->GetName())) ; 
@@ -1292,7 +1364,7 @@ Bool_t AliQAManager::SaveIt2OCDB(const Int_t runNumber, TFile * inputFile, const
                                TList * histList = esDir->GetListOfKeys() ; 
                                TIter nextHist(histList) ; 
                                TKey * histKey ; 
-                               while ( (histKey = dynamic_cast<TKey*>(nextHist())) ) {
+                               while ( (histKey = static_cast<TKey*>(nextHist())) ) {
                                        TObject * odata = esDir->Get(histKey->GetName()) ; 
                                        if ( !odata ) {
                                                AliError(Form("%s in %s/%s returns a NULL pointer !!", histKey->GetName(), detDir->GetName(), taskDir->GetName())) ;
@@ -1302,7 +1374,7 @@ Bool_t AliQAManager::SaveIt2OCDB(const Int_t runNumber, TFile * inputFile, const
               TList * expertHistList = expertDir->GetListOfKeys() ; 
               TIter nextExpertHist(expertHistList) ; 
               TKey * expertHistKey ; 
-              while ( (expertHistKey = dynamic_cast<TKey*>(nextExpertHist())) ) {
+              while ( (expertHistKey = static_cast<TKey*>(nextExpertHist())) ) {
                 TObject * expertOdata = expertDir->Get(expertHistKey->GetName()) ; 
                 if ( !expertOdata ) {
                   AliError(Form("%s in %s/%s/Expert returns a NULL pointer !!", expertHistKey->GetName(), detDir->GetName(), taskDir->GetName())) ;
@@ -1335,6 +1407,7 @@ Bool_t AliQAManager::SaveIt2OCDB(const Int_t runNumber, TFile * inputFile, const
 void AliQAManager::SetEventSpecie(AliRecoParam::EventSpecie_t es) 
 {
   // set the current event specie and inform AliQAv1 that this event specie has been encountered
+  fEventSpecie = es ; 
   AliQAv1::Instance()->SetEventSpecie(es) ; 
 }
 
@@ -1355,3 +1428,15 @@ void AliQAManager::SetWriteExpert()
       fQAWriteExpert[iDet] = kTRUE ;
   }
 }  
+
+//_____________________________________________________________________________
+void AliQAManager::Destroy() {
+// delete AliQAManager instance and
+// all associated objects
+
+  if (fgQAInstance) {
+    delete fgQAInstance ;
+    fgQAInstance = NULL ;
+  }
+}
+