Removing the hard-wired particle masses (B. Hippolyte)

[u/mrichter/AliRoot.git] / STEER / AliQADataMakerSteer.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id$ */
///////////////////////////////////////////////////////////////////////////////
//                                                                           //
// class for running the QA makers                                           //
//                                                                           //
//   AliQADataMakerSteer qas;                                                //
//   qas.Run(AliQA::kRAWS, rawROOTFileName);                                 //
//   qas.Run(AliQA::kHITS);                                                  //
//   qas.Run(AliQA::kSDIGITS);                                               //
//   qas.Run(AliQA::kDIGITS);                                                //
//   qas.Run(AliQA::kRECPOINTS);                                             //
//   qas.Run(AliQA::kESDS);                                                  //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

#include <TKey.h>
#include <TFile.h>
#include <TFileMerger.h>
#include <TPluginManager.h>
#include <TROOT.h>
#include <TString.h>
#include <TSystem.h>

#include "AliCDBManager.h"
#include "AliCDBEntry.h"
#include "AliCDBId.h"
#include "AliCDBMetaData.h"
#include "AliCodeTimer.h"
#include "AliCorrQADataMakerRec.h"
#include "AliDetectorRecoParam.h"
#include "AliESDEvent.h"
#include "AliGeomManager.h"
#include "AliGlobalQADataMaker.h"
#include "AliHeader.h"
#include "AliLog.h"
#include "AliModule.h"
#include "AliQA.h"
#include "AliQADataMakerRec.h"
#include "AliQADataMakerSim.h"
#include "AliQADataMakerSteer.h" 
#include "AliRawReaderDate.h"
#include "AliRawReaderFile.h"
#include "AliRawReaderRoot.h"
#include "AliRun.h"
#include "AliRunLoader.h"

ClassImp(AliQADataMakerSteer) 

//_____________________________________________________________________________
AliQADataMakerSteer::AliQADataMakerSteer(const Char_t * mode, const Char_t* gAliceFilename, const Char_t * name, const Char_t * title) :
	TNamed(name, title), 
	fCurrentEvent(0),  
	fCycleSame(kFALSE),
	fDetectors("ALL"), 
	fDetectorsW("ALL"), 
	fESD(NULL), 
	fESDTree(NULL),
	fGAliceFileName(gAliceFilename), 
	fFirstEvent(0),        
	fMaxEvents(0),   
  fMode(mode), 
	fNumberOfEvents(999999), 
  fRecoParam(),
	fRunNumber(0), 
	fRawReader(NULL), 
	fRawReaderDelete(kTRUE), 
	fRunLoader(NULL), 
  fTasks(""), 
  fEventSpecie(AliRecoParam::kDefault)
{
	// default ctor
	fMaxEvents = fNumberOfEvents ; 
	for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
		if (IsSelected(AliQA::GetDetName(iDet))) {
			fLoader[iDet]      = NULL ;
			fQADataMaker[iDet] = NULL ;
			fQACycles[iDet]    = 999999 ;
      fQAWriteExpert[iDet] = kTRUE ;
		}
	}	
}

//_____________________________________________________________________________
AliQADataMakerSteer::AliQADataMakerSteer(const AliQADataMakerSteer & qas) : 
	TNamed(qas), 
	fCurrentEvent(qas.fCurrentEvent),  
	fCycleSame(kFALSE),
	fDetectors(qas.fDetectors), 
	fDetectorsW(qas.fDetectorsW), 
	fESD(NULL), 
	fESDTree(NULL), 
	fGAliceFileName(qas.fGAliceFileName), 
	fFirstEvent(qas.fFirstEvent),        
	fMaxEvents(qas.fMaxEvents),    
	fMode(qas.fMode), 
	fNumberOfEvents(qas.fNumberOfEvents), 
	fRecoParam(),		
	fRunNumber(qas.fRunNumber), 
	fRawReader(NULL), 
	fRawReaderDelete(kTRUE), 
	fRunLoader(NULL), 
  fTasks(qas.fTasks),
  fEventSpecie(qas.fEventSpecie)
{
	// cpy ctor
	for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
		fLoader[iDet]         = qas.fLoader[iDet] ;
		fQADataMaker[iDet]    = qas.fQADataMaker[iDet] ;
		fQACycles[iDet]       = qas.fQACycles[iDet] ;	
    fQAWriteExpert[iDet] = qas.fQAWriteExpert[iDet] ;
	}
}

//_____________________________________________________________________________
AliQADataMakerSteer & AliQADataMakerSteer::operator = (const AliQADataMakerSteer & qas) 
{
	// assignment operator
  this->~AliQADataMakerSteer() ;
  new(this) AliQADataMakerSteer(qas) ;
  return *this ;
}

//_____________________________________________________________________________
AliQADataMakerSteer::~AliQADataMakerSteer() 
{
	// dtor
	for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
		if (IsSelected(AliQA::GetDetName(iDet))) {
		  fLoader[iDet] = NULL;
		  if (fQADataMaker[iDet]) {
			  (fQADataMaker[iDet])->Finish() ; 
				delete fQADataMaker[iDet] ;
		  }
		}
	}
	if (fRawReaderDelete) { 
		fRunLoader = NULL ;
		delete fRawReader ;
		fRawReader = NULL ;
	}
}

//_____________________________________________________________________________
Bool_t AliQADataMakerSteer::DoIt(const AliQA::TASKINDEX_t taskIndex)
{
	// Runs all the QA data Maker for every detector
		
	Bool_t rv = kFALSE ;
    // Fill QA data in event loop 
	for (UInt_t iEvent = fFirstEvent ; iEvent < (UInt_t)fMaxEvents ; iEvent++) {
		fCurrentEvent++ ; 
		// Get the event
		if ( iEvent%10 == 0  ) 
			AliInfo(Form("processing event %d", iEvent));
		if ( taskIndex == AliQA::kRAWS ) {
			if ( !fRawReader->NextEvent() )
				break ;
		} else if ( taskIndex == AliQA::kESDS ) {
			if ( fESDTree->GetEntry(iEvent) == 0 )
				break ;
		} else {
			if ( fRunLoader->GetEvent(iEvent) != 0 )
				break ;
		}
		// loop  over active loaders
		for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
			if (IsSelected(AliQA::GetDetName(iDet))) {
				AliQADataMaker * qadm = GetQADataMaker(iDet) ;
				if (!qadm) continue; // This detector doesn't have any QA (for example, HLT)
				if ( qadm->IsCycleDone() ) {
          qadm->EndOfCycle(taskIndex) ;
				}
				TTree * data = NULL ; 
				AliLoader* loader = GetLoader(qadm->GetUniqueID());
				switch (taskIndex) {
					case AliQA::kNULLTASKINDEX : 
						break ; 
					case AliQA::kRAWS :
						qadm->Exec(taskIndex, fRawReader) ; 
						break ; 
					case AliQA::kHITS :
            if( loader ) {
              loader->LoadHits() ; 
              data = loader->TreeH() ; 
              if ( ! data ) {
                AliWarning(Form(" Hit Tree not found for  %s", AliQA::GetDetName(iDet))) ; 
                break ; 
              } 
            } 
            qadm->Exec(taskIndex, data) ;
						break ;
						case AliQA::kSDIGITS :
            if( loader ) {      
              loader->LoadSDigits() ; 
              data = loader->TreeS() ; 
              if ( ! data ) {
                AliWarning(Form(" SDigit Tree not found for  %s", AliQA::GetDetName(iDet))) ; 
                break ; 
              } 
            }
            qadm->Exec(taskIndex, data) ; 
						break; 
						case AliQA::kDIGITS :
            if( loader ) {      
              loader->LoadDigits() ; 
              data = loader->TreeD() ; 
              if ( ! data ) {
                AliWarning(Form(" Digit Tree not found for  %s", AliQA::GetDetName(iDet))) ; 
                break ; 
              } 
            }
            qadm->Exec(taskIndex, data) ;
						break; 
						case AliQA::kRECPOINTS :
            if( loader ) {      
              loader->LoadRecPoints() ; 
              data = loader->TreeR() ; 
              if (!data) {
                AliWarning(Form("RecPoints not found for %s", AliQA::GetDetName(iDet))) ; 
                break ; 
              } 
            }
            qadm->Exec(taskIndex, data) ; 
            break; 
						case AliQA::kTRACKSEGMENTS :
						break; 
						case AliQA::kRECPARTICLES :
						break; 
						case AliQA::kESDS :
						qadm->Exec(taskIndex, fESD) ;
						break; 
						case AliQA::kNTASKINDEX :
						break; 
				} //task switch
			}
		} // detector loop
    Increment() ; 
	} // event loop	
	// Save QA data for all detectors
	rv = Finish(taskIndex) ;
	
	if ( taskIndex == AliQA::kRAWS ) 
		fRawReader->RewindEvents() ;

	return rv ; 
}

//_____________________________________________________________________________
Bool_t AliQADataMakerSteer::Finish(const AliQA::TASKINDEX_t taskIndex) 
{
	// write output to file for all detectors
	for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
		if (IsSelected(AliQA::GetDetName(iDet))) {
			AliQADataMaker * qadm = GetQADataMaker(iDet) ;
			if (qadm) 
        qadm->EndOfCycle(taskIndex) ;
		}
	}
	return kTRUE ; 
}

//_____________________________________________________________________________
TObjArray * AliQADataMakerSteer::GetFromOCDB(AliQA::DETECTORINDEX_t det, AliQA::TASKINDEX_t task, const char * year) const 
{
	// Retrieve the list of QA data for a given detector and a given task 
	TObjArray * rv = NULL ;
	if ( !strlen(AliQA::GetQARefStorage()) ) { 
		AliError("No storage defined, use AliQA::SetQARefStorage") ; 
		return NULL ; 
	}	
	AliCDBManager* man = AliCDBManager::Instance() ; 
	if ( ! man->IsDefaultStorageSet() ) {
		TString tmp(AliQA::GetQARefDefaultStorage()) ; 
		tmp.Append(year) ; 
		tmp.Append("/") ; 
		man->SetDefaultStorage(tmp.Data()) ; 		
		man->SetSpecificStorage(Form("%s/*", AliQA::GetQAName()), AliQA::GetQARefStorage()) ;
	}
	TString detOCDBDir(Form("%s/%s/%s", AliQA::GetQAName(), AliQA::GetDetName((Int_t)det), AliQA::GetRefOCDBDirName())) ; 
	AliInfo(Form("Retrieving reference data from %s/%s for %s", AliQA::GetQARefStorage(), detOCDBDir.Data(), AliQA::GetTaskName(task).Data())) ; 
	AliCDBEntry* entry = man->Get(detOCDBDir.Data(), 0) ; //FIXME 0 --> Run Number
	TList * listDetQAD = dynamic_cast<TList *>(entry->GetObject()) ;
	if ( listDetQAD ) 
		rv = dynamic_cast<TObjArray *>(listDetQAD->FindObject(AliQA::GetTaskName(task))) ; 
	return rv ; 
}

//_____________________________________________________________________________
AliLoader * AliQADataMakerSteer::GetLoader(Int_t iDet)
{
	// get the loader for a detector

	if ( !fRunLoader || iDet == AliQA::kCORR) 
		return NULL ; 
	
	TString detName = AliQA::GetDetName(iDet) ;
    fLoader[iDet] = fRunLoader->GetLoader(detName + "Loader");
	if (fLoader[iDet]) 
		return fLoader[iDet] ;
	
	// load the QA data maker object
	TPluginManager* pluginManager = gROOT->GetPluginManager() ;
	TString loaderName = "Ali" + detName + "Loader" ;

	AliLoader * loader = NULL ;
	// first check if a plugin is defined for the quality assurance data maker
	TPluginHandler* pluginHandler = pluginManager->FindHandler("AliLoader", detName) ;
	// if not, add a plugin for it
	if (!pluginHandler) {
		AliDebug(1, Form("defining plugin for %s", loaderName.Data())) ;
		TString libs = gSystem->GetLibraries() ;
		if (libs.Contains("lib" + detName + "base.so") || (gSystem->Load("lib" + detName + "base.so") >= 0)) {
			pluginManager->AddHandler("AliQADataMaker", detName, loaderName, detName + "loader", loaderName + "()") ;
		} else {
			pluginManager->AddHandler("AliLoader", detName, loaderName, detName, loaderName + "()") ;
		}
		pluginHandler = pluginManager->FindHandler("AliLoader", detName) ;
	}
	if (pluginHandler && (pluginHandler->LoadPlugin() == 0)) {
		loader = (AliLoader *) pluginHandler->ExecPlugin(0) ;
	}
	if (loader) 
		fLoader[iDet] = loader ;
	return loader ;
}

//_____________________________________________________________________________
AliQADataMaker * AliQADataMakerSteer::GetQADataMaker(const Int_t iDet) 
{
	// get the quality assurance data maker for a detector
	
	if (fQADataMaker[iDet]) {
    fQADataMaker[iDet]->SetEventSpecie(fEventSpecie) ; 
		return fQADataMaker[iDet] ;
  }
	
	AliQADataMaker * qadm = NULL ;
	
	if (iDet == AliQA::kGLOBAL) { //Global QA
		qadm = new AliGlobalQADataMaker();
		qadm->SetName(AliQA::GetDetName(iDet));
		qadm->SetUniqueID(iDet);
		fQADataMaker[iDet] = qadm;
    qadm->SetEventSpecie(fEventSpecie) ; 
		return qadm;
	}

	if (iDet == AliQA::kCORR) { //the data maker for correlations among detectors
    qadm = new AliCorrQADataMakerRec(fQADataMaker) ; 
		qadm->SetName(AliQA::GetDetName(iDet));
		qadm->SetUniqueID(iDet);
		fQADataMaker[iDet] = qadm;
    qadm->SetEventSpecie(fEventSpecie) ; 
		return qadm;
  }

	// load the QA data maker object
	TPluginManager* pluginManager = gROOT->GetPluginManager() ;
	TString detName = AliQA::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(1, 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(AliQA::GetDetName(iDet));
		qadm->SetUniqueID(iDet);
		fQADataMaker[iDet] = qadm ;
    qadm->SetEventSpecie(fEventSpecie) ; 
	}

  return qadm ;
}

//_____________________________________________________________________________
void  AliQADataMakerSteer::EndOfCycle(TObjArray * detArray) 
{
	// End of cycle QADataMakers 
	
	for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
		if (IsSelected(AliQA::GetDetName(iDet))) {
			AliQADataMaker * qadm = GetQADataMaker(iDet) ;
			if (!qadm) 
				continue ;	
			// skip non active detectors
			if (detArray) {
				AliModule* det = static_cast<AliModule*>(detArray->FindObject(AliQA::GetDetName(iDet))) ;
				if (!det || !det->IsActive())  
					continue ;
			}
			for (UInt_t taskIndex = 0; taskIndex < AliQA::kNTASKINDEX; taskIndex++) {
				if ( fTasks.Contains(Form("%d", taskIndex)) ) 
					qadm->EndOfCycle(AliQA::GetTaskIndex(AliQA::GetTaskName(taskIndex))) ;
			}
			qadm->Finish();
		}
	}
}

//_____________________________________________________________________________
void  AliQADataMakerSteer::EndOfCycle(TString detectors) 
{
	// End of cycle QADataMakers 
	
	for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
		if (IsSelected(AliQA::GetDetName(iDet))) {
			AliQADataMaker * qadm = GetQADataMaker(iDet) ;
			if (!qadm) 
				continue ;	
			// skip non active detectors
      if (!detectors.Contains(AliQA::GetDetName(iDet))) 
        continue ;
   		for (UInt_t taskIndex = 0; taskIndex < AliQA::kNTASKINDEX; taskIndex++) {
				if ( fTasks.Contains(Form("%d", taskIndex)) ) 
					qadm->EndOfCycle(AliQA::GetTaskIndex(AliQA::GetTaskName(taskIndex))) ;
			}
			qadm->Finish();
		}
	}
}

//_____________________________________________________________________________
void AliQADataMakerSteer::Increment()
{
  // Increments the cycle counter for all QA Data Makers
 	for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
		if (IsSelected(AliQA::GetDetName(iDet))) {
			AliQADataMaker * qadm = GetQADataMaker(iDet) ;
			if (qadm) 
        qadm->Increment() ;
    }
  }
}
  
//_____________________________________________________________________________
Bool_t AliQADataMakerSteer::Init(const AliQA::TASKINDEX_t taskIndex, const  char * input )
{
	// Initialize the event source and QA data makers
	
	//fTasks = AliQA::GetTaskName(taskIndex) ; 

	if (taskIndex == AliQA::kRAWS) { 
		if (!fRawReader) {
		        fRawReader = AliRawReader::Create(input);
		}
		if ( ! fRawReader ) 
			return kFALSE ; 
		fRawReaderDelete = kTRUE ; 
		fRawReader->NextEvent() ; 
		fRunNumber = fRawReader->GetRunNumber() ; 
		AliCDBManager::Instance()->SetRun(fRunNumber) ; 
		fRawReader->RewindEvents();
		fNumberOfEvents = 999999 ;
		if ( fMaxEvents < 0 ) 
			fMaxEvents = fNumberOfEvents ; 
		} else if (taskIndex == AliQA::kESDS) {
			fTasks = AliQA::GetTaskName(AliQA::kESDS) ; 
      if (!gSystem->AccessPathName("AliESDs.root")) { // AliESDs.root exists
        TFile * esdFile = TFile::Open("AliESDs.root") ;
        fESDTree = dynamic_cast<TTree *> (esdFile->Get("esdTree")) ; 
        if ( !fESDTree ) {
          AliError("esdTree not found") ; 
          return kFALSE ; 
        } else {
          fESD     = new AliESDEvent() ;
          fESD->ReadFromTree(fESDTree) ;
          fESDTree->GetEntry(0) ; 
          fRunNumber = fESD->GetRunNumber() ; 
          fNumberOfEvents = fESDTree->GetEntries() ;
          if ( fMaxEvents < 0 ) 
            fMaxEvents = fNumberOfEvents ; 
        }
      } else {
        AliError("AliESDs.root not found") ; 
        return kFALSE ; 
      }			
    } else {
      if ( !InitRunLoader() ) { 
        AliWarning("No Run Loader not found") ; 
      } else {
        fNumberOfEvents = fRunLoader->GetNumberOfEvents() ;
        if ( fMaxEvents < 0 ) 
          fMaxEvents = fNumberOfEvents ; 
      }
    }

  // Get Detectors 
  TObjArray* detArray = NULL ; 
	if (fRunLoader) // check if RunLoader exists 
		if ( fRunLoader->GetAliRun() ) { // check if AliRun exists in gAlice.root
			detArray = fRunLoader->GetAliRun()->Detectors() ;
			fRunNumber = fRunLoader->GetHeader()->GetRun() ; 
		}

	// Initialize all QA data makers for all detectors
	fRunNumber = AliCDBManager::Instance()->GetRun() ; 
	if ( !  AliGeomManager::GetGeometry() ) 
		AliGeomManager::LoadGeometry() ; 
	
	InitQADataMaker(fRunNumber, detArray) ; //, fCycleSame, kTRUE, detArray) ; 
	return kTRUE ; 
}

//_____________________________________________________________________________
void  AliQADataMakerSteer::InitQADataMaker(UInt_t run, TObjArray * detArray) 
{
	// Initializes The QADataMaker for all active detectors and for all active tasks 
	for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
		if (IsSelected(AliQA::GetDetName(iDet))) {
			AliQADataMaker * qadm = GetQADataMaker(iDet) ;
			if (!qadm) {
				AliError(Form("AliQADataMaker not found for %s", AliQA::GetDetName(iDet))) ; 
				fDetectorsW.ReplaceAll(AliQA::GetDetName(iDet), "") ; 
			} else {
        if (fQAWriteExpert[iDet])
          qadm->SetWriteExpert() ; 
				AliDebug(1, Form("Data Maker found for %s", qadm->GetName())) ; 
				// skip non active detectors
				if (detArray) {
					AliModule* det = static_cast<AliModule*>(detArray->FindObject(AliQA::GetDetName(iDet))) ;
					if (!det || !det->IsActive())  
						continue ;
				}
				if (fQAWriteExpert[iDet]) qadm->SetWriteExpert() ; 
	      // Set default reco params
        Bool_t sameCycle = kFALSE ; 
				for (UInt_t taskIndex = 0; taskIndex < AliQA::kNTASKINDEX; taskIndex++) {
					if ( fTasks.Contains(Form("%d", taskIndex)) ) {
						qadm->Init(AliQA::GetTaskIndex(AliQA::GetTaskName(taskIndex)), GetQACycles(qadm->GetUniqueID())) ;
            qadm->StartOfCycle(AliQA::GetTaskIndex(AliQA::GetTaskName(taskIndex)), run,  sameCycle) ;
            sameCycle = kTRUE ;
					}
				}
			}
		}
	}
}


//_____________________________________________________________________________
Bool_t AliQADataMakerSteer::InitRunLoader()
{
	// get or create the run loader
	if (fRunLoader) {
		fCycleSame = kTRUE ; 
	} else {
		if (!gSystem->AccessPathName(fGAliceFileName.Data())) { // galice.root exists
			// load all base libraries to get the loader classes
			TString libs = gSystem->GetLibraries() ;
			for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
				if (!IsSelected(AliQA::GetDetName(iDet))) 
					continue ; 
				TString detName = AliQA::GetDetName(iDet) ;
				if (detName == "HLT") 
					continue;
				if (libs.Contains("lib" + detName + "base.so")) 
					continue;
				gSystem->Load("lib" + detName + "base.so");
			}
			fRunLoader = AliRunLoader::Open(fGAliceFileName.Data());
			if (!fRunLoader) {
				AliError(Form("no run loader found in file %s", fGAliceFileName.Data()));
				return kFALSE;
			}
			fRunLoader->CdGAFile();
			if (fRunLoader->LoadgAlice() == 0) {
				gAlice = fRunLoader->GetAliRun();
			}

			if (!gAlice) {
				AliError(Form("no gAlice object found in file %s", fGAliceFileName.Data()));
				return kFALSE;
			}

		} else {               // galice.root does not exist
			AliError(Form("the file %s does not exist", fGAliceFileName.Data()));
			return kFALSE;
		}
	}

	if (!fRunNumber) { 
		fRunLoader->LoadHeader();
		fRunNumber = fRunLoader->GetHeader()->GetRun() ; 
	}
	return kTRUE;
}

//_____________________________________________________________________________
Bool_t AliQADataMakerSteer::IsSelected(const char * det) 
{
	// check whether detName is contained in detectors
	// if yes, it is removed from detectors
	
	Bool_t rv = kFALSE;
	const TString detName(det) ;
  // always activates Correlation
  if ( detName.Contains(AliQA::GetDetName(AliQA::kCORR))) {
    rv = kTRUE ; 
  } else {
    // check if all detectors are selected
    if (fDetectors.Contains("ALL")) {
      fDetectors = "ALL";
      rv = kTRUE;
    } else if ((fDetectors.CompareTo(detName) == 0) ||
               fDetectors.BeginsWith(detName+" ") ||
               fDetectors.EndsWith(" "+detName) ||
               fDetectors.Contains(" "+detName+" ")) {
      rv = kTRUE;
    }
  }
	return rv ;
}

//_____________________________________________________________________________
Bool_t AliQADataMakerSteer::Merge(const Int_t runNumber) const
{
	// Merge data from all the cycles from all detectors in one single file per run
	// Merge the QA results from all the data chunks in one run 
 Bool_t rv = MergeData(runNumber) ; 
 rv *= MergeResults(runNumber) ;
 return rv ; 
}
	
	
//_____________________________________________________________________________
Bool_t AliQADataMakerSteer::MergeData(const Int_t runNumber) const
{
	// Merge all the cycles from all detectors in one single file per run
	TString cmd ;
	if (runNumber == -1) 
		cmd = Form(".! ls *%s*.%d.root > tempo.txt", AliQA::GetQADataFileName(), runNumber) ; 
	else 
		cmd = Form(".! ls *%s*.*.root > tempo.txt", AliQA::GetQADataFileName()) ; 
	gROOT->ProcessLine(cmd.Data()) ;
	ifstream in("tempo.txt") ; 
	const Int_t runMax = 10 ;  
	TString file[AliQA::kNDET*runMax] ;
	
	Int_t index = 0 ; 
	while ( 1 ) {
		in >> file[index] ; 
		if ( !in.good() ) 
			break ; 
		AliInfo(Form("index = %d file = %s", index, (file[index]).Data())) ; 
		index++ ;
	}
	
	if ( index == 0 ) { 
		AliError(Form("run number %d not found", runNumber)) ; 
		return kFALSE ; 
	}

  TFileMerger merger ; 
  TString outFileName(Form("Merged.%s.Data.%d.root",AliQA::GetQADataFileName(),runNumber)); 
  merger.OutputFile(outFileName.Data()) ; 
  for (Int_t ifile = 0 ; ifile < index-1 ; ifile++) {
    TString pattern(Form("%s.%d.", AliQA::GetQADataFileName(), runNumber)); 
    TString tmp(file[ifile]) ; 
    if (tmp.Contains(pattern)) {
      merger.AddFile(tmp) ; 
    }
	}
  merger.Merge() ; 
	return kTRUE ; 
}

//_____________________________________________________________________________
Bool_t AliQADataMakerSteer::MergeResults(const Int_t runNumber) const
{
	// Merge the QA result from all the data chunks in a run 
	TString cmd ;
	cmd = Form(".! ls %s*.root > tempo.txt", AliQA::GetQADataFileName()) ; 
	gROOT->ProcessLine(cmd.Data()) ;
	ifstream in("tempo.txt") ; 
	const Int_t chunkMax = 100 ;  
	TString fileList[chunkMax] ;
	
	Int_t index = 0 ; 
	while ( 1 ) {
		TString file ; 
		in >> fileList[index] ; 
		if ( !in.good() ) 
			break ; 
		AliInfo(Form("index = %d file = %s", index, (fileList[index].Data()))) ; 
		index++ ;
	}
	
	if ( index == 0 ) { 
		AliError("No QA Result File found") ; 
		return kFALSE ; 
	}
	
	TFileMerger merger ; 
	TString outFileName(Form("Merged.%s.Result.%d.root", AliQA::GetQADataFileName(), runNumber));		
	merger.OutputFile(outFileName.Data()) ; 
	for (Int_t ifile = 0 ; ifile < index ; ifile++) {
		TString file = fileList[ifile] ; 
		merger.AddFile(file) ; 
	}
	merger.Merge() ; 
	
	return kTRUE ; 
}

//_____________________________________________________________________________
void AliQADataMakerSteer::Reset(const Bool_t sameCycle)
{
	// Reset the default data members

	for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
		if (IsSelected(AliQA::GetDetName(iDet))) {
			AliQADataMaker * qadm = GetQADataMaker(iDet);
			qadm->Reset();
		}
	} 
	if (fRawReaderDelete) { 
		delete fRawReader ;
		fRawReader      = NULL ;
	}

	fCycleSame      = sameCycle ; 
	fESD            = NULL ; 
	fESDTree        = NULL ; 
	//fFirst          = kTRUE ;   
	fNumberOfEvents = 999999 ;  
}

//_____________________________________________________________________________
TString AliQADataMakerSteer::Run(const char * detectors, AliRawReader * rawReader, const Bool_t sameCycle) 
{
	//Runs all the QA data Maker for Raws only
	
	fCycleSame       = sameCycle ;
	fRawReader       = rawReader ;
	fDetectors       = detectors ; 
	fDetectorsW      = detectors ; 	
	
	AliCDBManager* man = AliCDBManager::Instance() ; 

	if ( man->GetRun() == -1 ) {// check if run number not set previously and set it from raw data
		rawReader->NextEvent() ; 
		man->SetRun(fRawReader->GetRunNumber()) ;
		rawReader->RewindEvents() ;
	}	
	
	if (!fCycleSame) 
    if ( !Init(AliQA::kRAWS) ) 
      return kFALSE ; 
  fRawReaderDelete = kFALSE ; 

	DoIt(AliQA::kRAWS) ; 
	return 	fDetectorsW ;
}

//_____________________________________________________________________________
TString AliQADataMakerSteer::Run(const char * detectors, const char * fileName, const Bool_t sameCycle) 
{
	//Runs all the QA data Maker for Raws only

	fCycleSame       = sameCycle ;
	fDetectors       = detectors ; 
	fDetectorsW      = detectors ; 	
	
	AliCDBManager* man = AliCDBManager::Instance() ; 
	if ( man->GetRun() == -1 ) { // check if run number not set previously and set it from AliRun
		AliRunLoader * rl = AliRunLoader::Open("galice.root") ;
		if ( ! rl ) {
			AliFatal("galice.root file not found in current directory") ; 
		} else {
			rl->CdGAFile() ; 
			rl->LoadgAlice() ;
			if ( ! rl->GetAliRun() ) {
				AliFatal("AliRun not found in galice.root") ;
			} else {
				rl->LoadHeader() ;
				man->SetRun(rl->GetHeader()->GetRun());
			}
		}
	}
	
	if (!fCycleSame) 
    if ( !Init(AliQA::kRAWS, fileName) ) 
      return kFALSE ; 
	
	DoIt(AliQA::kRAWS) ; 
	return 	fDetectorsW ;
}

//_____________________________________________________________________________
TString AliQADataMakerSteer::Run(const char * detectors, const AliQA::TASKINDEX_t taskIndex, Bool_t const sameCycle, const  char * fileName ) 
{
	// Runs all the QA data Maker for every detector
	
	fCycleSame       = sameCycle ;
	fDetectors       = detectors ; 
	fDetectorsW      = detectors ; 		
	
	AliCDBManager* man = AliCDBManager::Instance() ; 	
	if ( man->GetRun() == -1 ) { // check if run number not set previously and set it from AliRun
		AliRunLoader * rl = AliRunLoader::Open("galice.root") ;
		if ( ! rl ) {
			AliFatal("galice.root file not found in current directory") ; 
		} else {
			rl->CdGAFile() ; 
			rl->LoadgAlice() ;
			if ( ! rl->GetAliRun() ) {
				AliInfo("AliRun not found in galice.root") ;
			} else {
				rl->LoadHeader() ;
				man->SetRun(rl->GetHeader()->GetRun()) ;
			}
		}
	}
	

	if ( taskIndex == AliQA::kNULLTASKINDEX) { 
		for (UInt_t task = 0; task < AliQA::kNTASKINDEX; task++) {
			if ( fTasks.Contains(Form("%d", task)) ) {
        if (!fCycleSame)
          if ( !Init(AliQA::GetTaskIndex(AliQA::GetTaskName(task)), fileName) ) 
            return kFALSE ;
        DoIt(AliQA::GetTaskIndex(AliQA::GetTaskName(task))) ;
			}
		}
	} else {
    if (! fCycleSame )
      if ( !Init(taskIndex, fileName) ) 
        return kFALSE ; 
      DoIt(taskIndex) ; 
  } 		
	
	return fDetectorsW ;

}

//_____________________________________________________________________________
void AliQADataMakerSteer::RunOneEvent(AliRawReader * rawReader) 
{
	//Runs all the QA data Maker for Raws only and on one event only (event loop done by calling method)
  if ( ! rawReader ) 
    return ; 
	AliCodeTimerAuto("") ;
  if (fTasks.Contains(Form("%d", AliQA::kRAWS))){
    for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
      if (!IsSelected(AliQA::GetDetName(iDet))) 
        continue;
      AliQADataMaker *qadm = GetQADataMaker(iDet);  
      if (!qadm) 
        continue;
      if ( qadm->IsCycleDone() ) {
        qadm->EndOfCycle() ;
      }
      AliCodeTimerStart(Form("running RAW quality assurance data maker for %s", AliQA::GetDetName(iDet))); 
      qadm->SetEventSpecie(fEventSpecie) ; 
			qadm->Exec(AliQA::kRAWS, rawReader) ;
      AliCodeTimerStop(Form("running RAW quality assurance data maker for %s", AliQA::GetDetName(iDet)));
		}
  }
}

//_____________________________________________________________________________
void AliQADataMakerSteer::RunOneEvent(AliESDEvent *& esd) 
{
	//Runs all the QA data Maker for ESDs only and on one event only (event loop done by calling method)
	
  AliCodeTimerAuto("") ;
  if (fTasks.Contains(Form("%d", AliQA::kESDS))) {
    for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
      if (!IsSelected(AliQA::GetDetName(iDet))) 
        continue;
      AliQADataMaker *qadm = GetQADataMaker(iDet);  
      if (!qadm) 
        continue;
      if ( qadm->IsCycleDone() ) {
        qadm->EndOfCycle() ;
      }
      AliCodeTimerStart(Form("running ESD quality assurance data maker for %s", AliQA::GetDetName(iDet)));
			qadm->Exec(AliQA::kESDS, esd) ;
      AliCodeTimerStop(Form("running ESD quality assurance data maker for %s", AliQA::GetDetName(iDet)));
		}
	}
}

//_____________________________________________________________________________
void AliQADataMakerSteer::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)
	AliCodeTimerAuto("") ;
  if (fTasks.Contains(Form("%d", AliQA::kRECPOINTS))) {
    if (IsSelected(AliQA::GetDetName(det))) {
      AliQADataMaker *qadm = GetQADataMaker(det);  
      if (qadm) { 
        if ( qadm->IsCycleDone() ) {
          qadm->EndOfCycle() ;
        }
        AliCodeTimerStart(Form("running RecPoints quality assurance data maker for %s", AliQA::GetDetName(det)));
        qadm->Exec(AliQA::kRECPOINTS, tree) ;
        AliCodeTimerStop(Form("running RecPoints quality assurance data maker for %s", AliQA::GetDetName(det)));
      }
    }
  }
}

//_____________________________________________________________________________
Bool_t AliQADataMakerSteer::Save2OCDB(const Int_t runNumber, AliRecoParam::EventSpecie_t es, const char * year, const char * detectors) const
{
	// take the locasl QA data merge into a single file and save in OCDB 
	Bool_t rv = kTRUE ; 
	TString tmp(AliQA::GetQARefStorage()) ; 
	if ( tmp.IsNull() ) { 
		AliError("No storage defined, use AliQA::SetQARefStorage") ; 
		return kFALSE ; 
	}
	if ( !(tmp.Contains(AliQA::GetLabLocalOCDB()) || tmp.Contains(AliQA::GetLabAliEnOCDB())) ) {
		AliError(Form("%s is a wrong storage, use %s or %s", AliQA::GetQARefStorage(), AliQA::GetLabLocalOCDB().Data(), AliQA::GetLabAliEnOCDB().Data())) ; 
		return kFALSE ; 
	}
	TString sdet(detectors) ; 
	sdet.ToUpper() ;
	TFile * inputFile ; 
	if ( sdet.Contains("ALL") ) {
		rv = Merge(runNumber) ; 
		if ( ! rv )
			return kFALSE ; 
		TString inputFileName(Form("Merged.%s.Data.%d.root", AliQA::GetQADataFileName(), runNumber)) ; 
		inputFile = TFile::Open(inputFileName.Data()) ; 
		rv = SaveIt2OCDB(runNumber, inputFile, year, es) ; 
	} else {
		for (Int_t index = 0; index < AliQA::kNDET; index++) {
			if (sdet.Contains(AliQA::GetDetName(index))) {
				TString inputFileName(Form("%s.%s.%d.root", AliQA::GetDetName(index), AliQA::GetQADataFileName(), runNumber)) ; 
				inputFile = TFile::Open(inputFileName.Data()) ; 			
				rv *= SaveIt2OCDB(runNumber, inputFile, year, es) ; 
			}
		}
	}
	return rv ; 
}

//_____________________________________________________________________________
Bool_t AliQADataMakerSteer::SaveIt2OCDB(const Int_t runNumber, TFile * inputFile, const char * 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 ;
	AliInfo(Form("Saving TH1s in %s to %s", inputFile->GetName(), AliQA::GetQARefStorage())) ; 
	AliCDBManager* man = AliCDBManager::Instance() ; 
	if ( ! man->IsDefaultStorageSet() ) {
		TString tmp( AliQA::GetQARefStorage() ) ; 
		if ( tmp.Contains(AliQA::GetLabLocalOCDB()) ) 
			man->SetDefaultStorage(AliQA::GetQARefStorage()) ;
		else {
			TString tmp1(AliQA::GetQARefDefaultStorage()) ; 
			tmp1.Append(year) ; 
			tmp1.Append("?user=alidaq") ; 
			man->SetDefaultStorage(tmp1.Data()) ; 
		}
	}
	man->SetSpecificStorage("*", AliQA::GetQARefStorage()) ; 
	if(man->GetRun() < 0) 
		man->SetRun(runNumber);

	AliCDBMetaData mdr ;
	mdr.SetResponsible("yves schutz");

	for ( Int_t detIndex = 0 ; detIndex < AliQA::kNDET ; detIndex++) {
		TDirectory * detDir = inputFile->GetDirectory(AliQA::GetDetName(detIndex)) ; 
		if ( detDir ) {
			AliInfo(Form("Entering %s", detDir->GetName())) ;
      AliQA::SetQARefDataDirName(es) ;
			TString detOCDBDir(Form("%s/%s/%s", AliQA::GetDetName(detIndex), AliQA::GetRefOCDBDirName(), AliQA::GetRefDataDirName())) ; 
			AliCDBId idr(detOCDBDir.Data(), runNumber, AliCDBRunRange::Infinity())  ;
			TList * listDetQAD = new TList() ;
			TString listName(Form("%s QA data Reference", AliQA::GetDetName(detIndex))) ; 
			mdr.SetComment(Form("%s QA stuff", AliQA::GetDetName(detIndex)));
			listDetQAD->SetName(listName) ; 
			TList * taskList = detDir->GetListOfKeys() ; 
			TIter nextTask(taskList) ; 
			TKey * taskKey ; 
			while ( (taskKey = dynamic_cast<TKey*>(nextTask())) ) {
				TDirectory * taskDir = detDir->GetDirectory(taskKey->GetName()) ; 
        TDirectory * esDir   = taskDir->GetDirectory(AliRecoParam::GetEventSpecieName(es)) ; 
				AliInfo(Form("Saving %s", esDir->GetName())) ; 
				TObjArray * listTaskQAD = new TObjArray(100) ; 
				listTaskQAD->SetName(Form("%s/%s", taskKey->GetName(), AliRecoParam::GetEventSpecieName(es))) ;
				listDetQAD->Add(listTaskQAD) ; 
				TList * histList = esDir->GetListOfKeys() ; 
				TIter nextHist(histList) ; 
				TKey * histKey ; 
				while ( (histKey = dynamic_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())) ;
					} else {
            if ( AliQA::GetExpert() == histKey->GetName() ) {
              TDirectory * expertDir   = esDir->GetDirectory(histKey->GetName()) ; 
              TList * expertHistList = expertDir->GetListOfKeys() ; 
              TIter nextExpertHist(expertHistList) ; 
              TKey * expertHistKey ; 
              while ( (expertHistKey = dynamic_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())) ;
                } else {
                  AliInfo(Form("Adding %s", expertHistKey->GetName())) ;
                  if ( expertOdata->IsA()->InheritsFrom("TH1") ) {
                    AliInfo(Form("Adding %s", expertHistKey->GetName())) ;
                    TH1 * hExpertdata = static_cast<TH1*>(expertOdata) ; 
                    listTaskQAD->Add(hExpertdata) ; 
                  }                  
                }                
              }
            }
						AliInfo(Form("Adding %s", histKey->GetName())) ;
						if ( odata->IsA()->InheritsFrom("TH1") ) {
							AliInfo(Form("Adding %s", histKey->GetName())) ;
							TH1 * hdata = static_cast<TH1*>(odata) ; 
							listTaskQAD->Add(hdata) ; 
						}
					}
				}
			}
			man->Put(listDetQAD, idr, &mdr) ;
		}
	}
	return rv ; 
}	

//_____________________________________________________________________________
void AliQADataMakerSteer::SetEventSpecie(AliRecoParam::EventSpecie_t es) 
{
  // set the current event specie and inform AliQA that this event specie has been encountered
  fEventSpecie = es ;
  AliQA::Instance()->SetEventSpecie(es) ; 
}

//_____________________________________________________________________________
void AliQADataMakerSteer::SetRecoParam(const Int_t det, const AliDetectorRecoParam *par) 
{
  // Set custom reconstruction parameters for a given detector
  // Single set of parameters for all the events
  GetQADataMaker(det)->SetRecoParam(par) ; 
}