[u/mrichter/AliRoot.git] / STEER / AliQAChecker.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 Quality Assurance Checker
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

#include "AliCDBEntry.h"
#include "AliQAManager.h"
#include "AliCDBStorage.h"
#include "AliRunInfo.h" 
#include "AliLog.h"
#include "AliModule.h" 
#include "AliQAv1.h"
#include "AliQAChecker.h"
#include "AliQACheckerBase.h"
#include "AliCorrQAChecker.h"
#include "AliGlobalQAChecker.h"
#include "AliGRPObject.h"

#include <TKey.h>
#include <TObjArray.h>
#include <TObjString.h>
#include <TPluginManager.h> 
#include <TROOT.h>
#include <TStopwatch.h> 
#include <TString.h> 
#include <TSystem.h> 
#include <TList.h>
#include <TNtupleD.h>

ClassImp(AliQAChecker)
  AliQAChecker * AliQAChecker::fgQAChecker = 0x0 ;

//_____________________________________________________________________________
AliQAChecker::AliQAChecker(const char* name, const char* title) :
  TNamed(name, title),
  fDataFile(0x0), 
  fRunInfo(0x0), 
  fRunInfoOwner(kFALSE), 
  fRefFile(0x0), 
  fFoundDetectors("."), 
  fEventSpecie(AliRecoParam::kDefault), 
  fRun(0)
{
  // ctor: initialise checkers and open the data file   
  for (Int_t det = 0 ; det < AliQAv1::kNDET ; det++) 
    fCheckers[det] = NULL ; 
}

//_____________________________________________________________________________
AliQAChecker::AliQAChecker(const AliQAChecker& qac) :
  TNamed(qac),
  fDataFile(qac.fDataFile), 
  fRunInfo(qac.fRunInfo), 
  fRunInfoOwner(kFALSE),   
  fRefFile(qac.fRefFile), 
  fFoundDetectors(qac.fFoundDetectors),
  fEventSpecie(qac.fEventSpecie),
  fRun(qac.fRun)
{
  // copy constructor
  
  for (Int_t det = 0 ; det < AliQAv1::kNDET ; det++) 
    fCheckers[det] = NULL ; 
}

//_____________________________________________________________________________
AliQAChecker& AliQAChecker::operator = (const AliQAChecker& qac)
{
// assignment operator

  this->~AliQAChecker();
  new(this) AliQAChecker(qac);
  return *this;
}

//_____________________________________________________________________________
AliQAChecker::~AliQAChecker()
{
// clean up
  if (fRunInfo)
    delete fRunInfo ; 
  delete [] fCheckers ; 
  AliQAv1::Close() ; 
}

//_____________________________________________________________________________
  AliQACheckerBase * AliQAChecker::GetDetQAChecker(Int_t det)
{
	// Gets the Quality Assurance checker for the detector specified by its name
	
	if (fCheckers[det]) 
    return fCheckers[det];

	AliQACheckerBase * qac = NULL ;

	TString detName(AliQAv1::GetDetName(det)) ; 
	
	if (det == AliQAv1::kGLOBAL) {
		qac = new AliGlobalQAChecker() ; 
	} else if (det == AliQAv1::kCORR) {
		qac = new AliCorrQAChecker() ; 
	} else {
		AliDebugClass(AliQAv1::GetQADebugLevel(), Form("Retrieving QA checker for %s", detName.Data())) ; 
		TPluginManager* pluginManager = gROOT->GetPluginManager() ;
		TString qacName = "Ali" + detName + "QAChecker" ;

		// first check if a plugin is defined for the quality assurance checker
		TPluginHandler* pluginHandler = pluginManager->FindHandler("AliQAChecker", detName.Data());
		// if not, add a plugin for it
		if (!pluginHandler) {
			//AliDebug(AliQAv1::GetQADebugLevel(), Form("defining plugin for %s", qacName.Data()));
			TString libs = gSystem->GetLibraries();
		
			if (libs.Contains("lib" + detName + "base.so") || (gSystem->Load("lib" + detName + "base.so") >= 0))
				pluginManager->AddHandler("AliQAChecker", detName, qacName, detName + "qac", qacName + "()");
			else 
				pluginManager->AddHandler("AliQAChecker", detName, qacName, detName, qacName + "()");

			pluginHandler = pluginManager->FindHandler("AliQAChecker", detName);	

			if (pluginHandler && (pluginHandler->LoadPlugin() == 0)) 
				qac = (AliQACheckerBase *) pluginHandler->ExecPlugin(0);
  
		}
	}
	if (qac) 
		fCheckers[det] = qac ;
	
	return qac ; 
}
 
//_____________________________________________________________________________
void AliQAChecker::GetRefSubDir(const char * det, const char * task, TDirectory *& dirFile, TObjArray **& dirOCDB)     
{ 
  // Opens and returns the file with the reference data 
  dirFile = NULL ; 
  TString refStorage(AliQAv1::GetQARefStorage()) ;
  if (!refStorage.Contains(AliQAv1::GetLabLocalOCDB()) && !refStorage.Contains(AliQAv1::GetLabAliEnOCDB())) {
    AliError(Form("%s is not a valid location for reference data", refStorage.Data())) ; 
    return ; 
  } else {
    AliQAManager* manQA = AliQAManager::QAManager(AliQAv1::GetTaskIndex(task)) ;
    dirOCDB = new TObjArray*[AliRecoParam::kNSpecies] ;	
    for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
      dirOCDB[specie] = NULL ; 
      if ( !AliQAv1::Instance()->IsEventSpecieSet(specie) ) 
        continue ; 
      AliQAv1::SetQARefDataDirName(specie) ;
      if ( ! manQA->GetLock() ) { 
        manQA->SetDefaultStorage(AliQAv1::GetQARefStorage()) ; 
        manQA->SetSpecificStorage("*", AliQAv1::GetQARefStorage()) ;
        manQA->SetRun(AliCDBManager::Instance()->GetRun()) ; 
        manQA->SetLock() ; 
      }
      char * detOCDBDir = Form("%s/%s/%s", det, AliQAv1::GetRefOCDBDirName(), AliQAv1::GetRefDataDirName()) ; 
      AliCDBEntry * entry = manQA->Get(detOCDBDir, manQA->GetRun()) ;
      if (entry) {
        TList * listDetQAD =static_cast<TList *>(entry->GetObject()) ;
        if ( strcmp(listDetQAD->ClassName(), "TList") != 0 ) {
          AliError(Form("Expected a Tlist and found a %s for detector %s", listDetQAD->ClassName(), det)) ; 
          listDetQAD = NULL ; 
          continue ; 
        } 
        if ( listDetQAD ) {
          TIter next(listDetQAD) ;
          TObjArray * ar ; 
          while ( (ar = (TObjArray*)next()) ) 
            dirOCDB[specie] = static_cast<TObjArray *>(listDetQAD->FindObject(Form("%s/%s", task, AliRecoParam::GetEventSpecieName(specie)))) ;             
        }
      }
    }
  }
}

//_____________________________________________________________________________
AliQAChecker * AliQAChecker::Instance()
{
	// returns unique instance of the checker
  if ( ! fgQAChecker ) 
   fgQAChecker = new AliQAChecker() ; 
  return fgQAChecker ;  
}

//_____________________________________________________________________________
void AliQAChecker::LoadRunInfoFromGRP() 
{
  AliCDBManager* man = AliCDBManager::Instance() ;
  AliCDBEntry* entry = man->Get(AliQAv1::GetGRPPath().Data());
  AliGRPObject* grpObject = 0x0;
  if (entry) {

	  TMap* m = static_cast<TMap*>(entry->GetObject());  // old GRP entry

	  if (m) {
	    AliDebug(AliQAv1::GetQADebugLevel(), "It is a map");
	    //m->Print();
	    grpObject = new AliGRPObject();
	         grpObject->ReadValuesFromMap(m);
    }

    else {
	    AliDebug(AliQAv1::GetQADebugLevel(), "It is a new GRP object");
        grpObject = static_cast<AliGRPObject*>(entry->GetObject());  // new GRP entry
    }

    entry->SetOwner(0);
    AliCDBManager::Instance()->UnloadFromCache("GRP/GRP/Data");
  }

  if (!grpObject) {
     AliFatal("No GRP entry found in OCDB!");
  }

  TString lhcState = grpObject->GetLHCState();
  if (lhcState==AliGRPObject::GetInvalidString()) {
    AliError("GRP/GRP/Data entry:  missing value for the LHC state ! Using UNKNOWN");
    lhcState = "UNKNOWN";
  }

  TString beamType = grpObject->GetBeamType();
  if (beamType==AliGRPObject::GetInvalidString()) {
    AliError("GRP/GRP/Data entry:  missing value for the beam type ! Using UNKNOWN");
    beamType = "UNKNOWN";
  }

  Float_t beamEnergy = grpObject->GetBeamEnergy();
  if (beamEnergy==AliGRPObject::GetInvalidFloat()) {
    AliError("GRP/GRP/Data entry:  missing value for the beam energy ! Using 0");
    beamEnergy = 0;
  }

  TString runType = grpObject->GetRunType();
  if (runType==AliGRPObject::GetInvalidString()) {
    AliError("GRP/GRP/Data entry:  missing value for the run type ! Using UNKNOWN");
    runType = "UNKNOWN";
  }

  Int_t activeDetectors = grpObject->GetDetectorMask();
  if (activeDetectors==AliGRPObject::GetInvalidInt()) {
    AliError("GRP/GRP/Data entry:  missing value for the detector mask ! Using 1074790399");
    activeDetectors = 1074790399;
  }

  fRunInfo = new AliRunInfo(lhcState, beamType, beamEnergy, runType, activeDetectors);

  fRunInfoOwner = kTRUE ; 

  // set the event specie
  fEventSpecie = AliRecoParam::kDefault ;
  if (strcmp(runType,"PHYSICS")) {
    // Not a physics run, the event specie is set to kCalib
    fEventSpecie = AliRecoParam::kCalib ;
    return;
  }
  if (strcmp(lhcState,"STABLE_BEAMS") == 0) {
    // Heavy ion run (any beam tha is not pp, the event specie is set to kHighMult
    fEventSpecie = AliRecoParam::kHighMult ;
    if ((strcmp(beamType,"p-p") == 0) ||
        (strcmp(beamType,"p-")  == 0) ||
        (strcmp(beamType,"-p")  == 0) ||
        (strcmp(beamType,"P-P") == 0) ||
        (strcmp(beamType,"P-")  == 0) ||
        (strcmp(beamType,"-P")  == 0)) {
      // Proton run, the event specie is set to kLowMult
      fEventSpecie = AliRecoParam::kLowMult ;
    }
    else if (strcmp(beamType,"-") == 0) {
      // No beams, we assume cosmic data
      fEventSpecie = AliRecoParam::kCosmic ;
    }
    else if (strcmp(beamType,"UNKNOWN") == 0) {
      // No LHC beam information is available, we use the default event specie
      fEventSpecie = AliRecoParam::kDefault ;
    }
  }
}

//_____________________________________________________________________________
Bool_t AliQAChecker::Run(const char * fileName, AliDetectorRecoParam * recoParam)
{
  // run the Quality Assurance Checker for all tasks Hits, SDigits, Digits, DigitsR, RecPoints, TrackSegments, RecParticles and ESDs
  // starting from data in file  
  TStopwatch stopwatch;
  stopwatch.Start();
  
  //search for all detectors QA directories
  TList * detKeyList = AliQAv1::GetQADataFile(fileName)->GetListOfKeys() ; 
  TIter nextd(detKeyList) ; 
  TKey * detKey ; 
  while ( (detKey = static_cast<TKey *>(nextd()) ) ) {
    AliDebug(AliQAv1::GetQADebugLevel(), Form("Found %s", detKey->GetName())) ;
    //Check which detector
    TString detName ; 
    TString detNameQA(detKey->GetName()) ; 
    Int_t det ; 
    for ( det = 0; det < AliQAv1::kNDET ; det++) {
      detName = AliQAv1::GetDetName(det) ; 
      if (detNameQA.Contains(detName)) {
        fFoundDetectors+=detName ; 
        fFoundDetectors+="." ;		
        break ; 
      }
    } 
    TDirectory * detDir = AliQAv1::GetQADataFile(fileName)->GetDirectory(detKey->GetName()) ; 
    TList * taskKeyList = detDir->GetListOfKeys() ;
    TIter nextt(taskKeyList) ; 
    TKey * taskKey ; 
    // now search for the tasks dir
    while ( (taskKey = static_cast<TKey *>(nextt()) ) ) {
      TString taskName( taskKey->GetName() ) ; 
      AliDebug(AliQAv1::GetQADebugLevel(), Form("Found %s", taskName.Data())) ;
      TDirectory * taskDir = detDir->GetDirectory(taskName.Data()) ; 
      taskDir->cd() ; 
      AliQACheckerBase * qac = GetDetQAChecker(det) ; 
      if (qac)
        AliDebug(AliQAv1::GetQADebugLevel(), Form("QA checker found for %s", detName.Data())) ; 
      if (!qac)
        AliFatal(Form("QA checker not found for %s", detName.Data())) ; 
      AliQAv1::ALITASK_t index = AliQAv1::kNULLTASK ; 
      if ( taskName == AliQAv1::GetTaskName(AliQAv1::kHITS) ) 
        index = AliQAv1::kSIM ; 
      if ( taskName == AliQAv1::GetTaskName(AliQAv1::kSDIGITS) ) 
        index = AliQAv1::kSIM ; 
      if ( taskName == AliQAv1::GetTaskName(AliQAv1::kDIGITS) ) 
        index = AliQAv1::kSIM ; 
      if ( taskName == AliQAv1::GetTaskName(AliQAv1::kRAWS) ) 
        index = AliQAv1::kRAW ;       
      if ( taskName == AliQAv1::GetTaskName(AliQAv1::kDIGITSR) ) 
        index = AliQAv1::kREC ; 
      if ( taskName == AliQAv1::GetTaskName(AliQAv1::kRECPOINTS) ) 
        index = AliQAv1::kREC ; 
      if ( taskName == AliQAv1::GetTaskName(AliQAv1::kTRACKSEGMENTS) ) 
        index = AliQAv1::kREC ; 
      if ( taskName == AliQAv1::GetTaskName(AliQAv1::kRECPARTICLES) ) 
        index = AliQAv1::kREC ; 
      if ( taskName == AliQAv1::GetTaskName(AliQAv1::kESDS) ) 
        index = AliQAv1::kESD ; 
      qac->Init(AliQAv1::DETECTORINDEX_t(det)) ; 
      
      TDirectory * refDir     = NULL ; 
      TObjArray ** refOCDBDir = NULL ;	
      GetRefSubDir(detNameQA.Data(), taskName.Data(), refDir, refOCDBDir) ;
		  qac->SetRefandData(refDir, refOCDBDir, taskDir) ;
		  qac->Run(index, recoParam) ; 
    }
  }
  TString detList ; 
  for ( Int_t det = 0; det < AliQAv1::kNDET; det++) {
    if (fFoundDetectors.Contains(AliQAv1::GetDetName(det))) {
      detList += AliQAv1::GetDetName(det) ; 
      detList += " " ; 
      fFoundDetectors.ReplaceAll(AliQAv1::GetDetName(det), "") ; 
      AliQAv1::Instance()->Show(AliQAv1::GetDetIndex(AliQAv1::GetDetName(det))) ; 
    }	
  }
  AliInfo(Form("QA performed for following detectors: %s", detList.Data())) ; 
  return kTRUE ; 
}

//_____________________________________________________________________________
Bool_t AliQAChecker::Run(AliQAv1::DETECTORINDEX_t det, AliQAv1::TASKINDEX_t task, TObjArray ** list, AliDetectorRecoParam * recoParam)
{
	// run the Quality Assurance Checker for detector det, for task task starting from data in list

	AliQACheckerBase * qac = GetDetQAChecker(det) ; 
	if (qac)
		AliDebug(AliQAv1::GetQADebugLevel(), Form("QA checker found for %s", AliQAv1::GetDetName(det).Data())) ;
	if (!qac)
		AliError(Form("QA checker not found for %s", AliQAv1::GetDetName(det).Data())) ; 
  
	AliQAv1::ALITASK_t index = AliQAv1::kNULLTASK ; 
	if ( task == AliQAv1::kRAWS ) 
		index = AliQAv1::kRAW ; 
	else if ( task == AliQAv1::kHITS ) 
		index = AliQAv1::kSIM ; 
	else if ( task == AliQAv1::kSDIGITS ) 
		index = AliQAv1::kSIM ; 
	else if ( task == AliQAv1::kDIGITS ) 
		index = AliQAv1::kSIM ; 
	else if ( task == AliQAv1::kDIGITSR ) 
		index = AliQAv1::kREC ; 
	else if ( task == AliQAv1::kRECPOINTS ) 
		index = AliQAv1::kREC ; 
	else if ( task == AliQAv1::kTRACKSEGMENTS ) 
		index = AliQAv1::kREC ; 
	else if ( task == AliQAv1::kRECPARTICLES ) 
		index = AliQAv1::kREC ; 
	else if ( task == AliQAv1::kESDS ) 
		index = AliQAv1::kESD ; 

	TDirectory * refDir     = NULL ; 
	TObjArray ** refOCDBDir = NULL  ;	
  qac->Init(det) ; 
  GetRefSubDir(AliQAv1::GetDetName(det), AliQAv1::GetTaskName(task), refDir, refOCDBDir) ;
  qac->SetRefandData(refDir, refOCDBDir) ; 
  qac->Run(index, list, recoParam) ; 
  
  // make the image 
  qac->MakeImage(list, task, AliQAv1::Mode(task)) ; 
  
	return kTRUE ; 
}

//_____________________________________________________________________________
Bool_t AliQAChecker::Run(AliQAv1::DETECTORINDEX_t det, AliQAv1::TASKINDEX_t task, TNtupleD ** list, AliDetectorRecoParam * recoParam)
{
	// run the Quality Assurance Checker for detector det, for task task starting from data in list
  
	AliQACheckerBase * qac = GetDetQAChecker(det) ; 
	if (qac)
		AliDebug(AliQAv1::GetQADebugLevel(), Form("QA checker found for %s", AliQAv1::GetDetName(det).Data())) ;
	if (!qac)
		AliError(Form("QA checker not found for %s", AliQAv1::GetDetName(det).Data())) ; 
  
	AliQAv1::ALITASK_t index = AliQAv1::kNULLTASK ; 
	if ( task == AliQAv1::kRAWS ) 
		index = AliQAv1::kRAW ; 
	else if ( task == AliQAv1::kHITS ) 
		index = AliQAv1::kSIM ; 
	else if ( task == AliQAv1::kSDIGITS ) 
		index = AliQAv1::kSIM ; 
	else if ( task == AliQAv1::kDIGITS ) 
		index = AliQAv1::kSIM ; 
	else if ( task == AliQAv1::kDIGITSR ) 
		index = AliQAv1::kREC ; 
	else if ( task == AliQAv1::kRECPOINTS ) 
		index = AliQAv1::kREC ; 
	else if ( task == AliQAv1::kTRACKSEGMENTS ) 
		index = AliQAv1::kREC ; 
	else if ( task == AliQAv1::kRECPARTICLES ) 
		index = AliQAv1::kREC ; 
	else if ( task == AliQAv1::kESDS ) 
		index = AliQAv1::kESD ; 
  
	TDirectory * refDir     = NULL ; 
	TObjArray ** refOCDBDir = NULL ;	
  qac->Init(det) ; 
  GetRefSubDir(AliQAv1::GetDetName(det), AliQAv1::GetTaskName(task), refDir, refOCDBDir) ;
  qac->SetRefandData(refDir, refOCDBDir) ; 
  qac->Run(index, list, recoParam) ; 

  return kTRUE ; 
}
Commit	Line	Data
421ab0fb	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	/* $Id: */
	17
	18	///////////////////////////////////////////////////////////////////////////////
	19	// //
	20	// class for running the Quality Assurance Checker
	21	// //
	22	///////////////////////////////////////////////////////////////////////////////
	23
4edbc5bc	24	#include "AliCDBEntry.h"
b03591ab	25	#include "AliQAManager.h"
96d67a8d	26	#include "AliCDBStorage.h"
1aaf4118	27	#include "AliRunInfo.h"
421ab0fb	28	#include "AliLog.h"
421ab0fb	29	#include "AliModule.h"
4e25ac79	30	#include "AliQAv1.h"
2e42b4d4	31	#include "AliQAChecker.h"
2e42b4d4	32	#include "AliQACheckerBase.h"
a2b64fbd	33	#include "AliCorrQAChecker.h"
401b84d6	34	#include "AliGlobalQAChecker.h"
44e45fac	35	#include "AliGRPObject.h"
421ab0fb	36
a5fa6165	37	#include <TKey.h>
421ab0fb	38	#include <TObjArray.h>
cbae8032	39	#include <TObjString.h>
a5fa6165	40	#include <TPluginManager.h>
a5fa6165	41	#include <TROOT.h>
421ab0fb	42	#include <TStopwatch.h>
421ab0fb	43	#include <TString.h>
a5fa6165	44	#include <TSystem.h>
a4976ef3	45	#include <TList.h>
a2b64fbd	46	#include <TNtupleD.h>
421ab0fb	47
2e42b4d4	48	ClassImp(AliQAChecker)
2e42b4d4	49	AliQAChecker * AliQAChecker::fgQAChecker = 0x0 ;
421ab0fb	50
421ab0fb	51	//_____________________________________________________________________________
2e42b4d4	52	AliQAChecker::AliQAChecker(const char* name, const char* title) :
421ab0fb	53	TNamed(name, title),
a5fa6165	54	fDataFile(0x0),
7e88424f	55	fRunInfo(0x0),
7e88424f	56	fRunInfoOwner(kFALSE),
4edbc5bc	57	fRefFile(0x0),
57acd2d2	58	fFoundDetectors("."),
634696f5	59	fEventSpecie(AliRecoParam::kDefault),
634696f5	60	fRun(0)
421ab0fb	61	{
a5fa6165	62	// ctor: initialise checkers and open the data file
4e25ac79	63	for (Int_t det = 0 ; det < AliQAv1::kNDET ; det++)
a5fa6165	64	fCheckers[det] = NULL ;
421ab0fb	65	}
	66
	67	//_____________________________________________________________________________
2e42b4d4	68	AliQAChecker::AliQAChecker(const AliQAChecker& qac) :
421ab0fb	69	TNamed(qac),
a5fa6165	70	fDataFile(qac.fDataFile),
7e88424f	71	fRunInfo(qac.fRunInfo),
7e88424f	72	fRunInfoOwner(kFALSE),
4edbc5bc	73	fRefFile(qac.fRefFile),
57acd2d2	74	fFoundDetectors(qac.fFoundDetectors),
634696f5	75	fEventSpecie(qac.fEventSpecie),
634696f5	76	fRun(qac.fRun)
421ab0fb	77	{
a5fa6165	78	// copy constructor
a5fa6165	79
4e25ac79	80	for (Int_t det = 0 ; det < AliQAv1::kNDET ; det++)
a5fa6165	81	fCheckers[det] = NULL ;
421ab0fb	82	}
	83
	84	//_____________________________________________________________________________
2e42b4d4	85	AliQAChecker& AliQAChecker::operator = (const AliQAChecker& qac)
421ab0fb	86	{
	87	// assignment operator
	88
2e42b4d4	89	this->~AliQAChecker();
2e42b4d4	90	new(this) AliQAChecker(qac);
421ab0fb	91	return *this;
	92	}
	93
	94	//_____________________________________________________________________________
2e42b4d4	95	AliQAChecker::~AliQAChecker()
421ab0fb	96	{
421ab0fb	97	// clean up
7e88424f	98	if (fRunInfo)
7e88424f	99	delete fRunInfo ;
a4976ef3	100	delete [] fCheckers ;
4e25ac79	101	AliQAv1::Close() ;
421ab0fb	102	}
	103
	104	//_____________________________________________________________________________
2e42b4d4	105	AliQACheckerBase * AliQAChecker::GetDetQAChecker(Int_t det)
421ab0fb	106	{
7ce961eb	107	// Gets the Quality Assurance checker for the detector specified by its name
	108
	109	if (fCheckers[det])
a5fa6165	110	return fCheckers[det];
a5fa6165	111
401b84d6	112	AliQACheckerBase * qac = NULL ;
a5fa6165	113
4e25ac79	114	TString detName(AliQAv1::GetDetName(det)) ;
401b84d6	115
4e25ac79	116	if (det == AliQAv1::kGLOBAL) {
401b84d6	117	qac = new AliGlobalQAChecker() ;
4e25ac79	118	} else if (det == AliQAv1::kCORR) {
a2b64fbd	119	qac = new AliCorrQAChecker() ;
401b84d6	120	} else {
634696f5	121	AliDebugClass(AliQAv1::GetQADebugLevel(), Form("Retrieving QA checker for %s", detName.Data())) ;
401b84d6	122	TPluginManager* pluginManager = gROOT->GetPluginManager() ;
	123	TString qacName = "Ali" + detName + "QAChecker" ;
	124
	125	// first check if a plugin is defined for the quality assurance checker
	126	TPluginHandler* pluginHandler = pluginManager->FindHandler("AliQAChecker", detName.Data());
	127	// if not, add a plugin for it
	128	if (!pluginHandler) {
5379c4a3	129	//AliDebug(AliQAv1::GetQADebugLevel(), Form("defining plugin for %s", qacName.Data()));
401b84d6	130	TString libs = gSystem->GetLibraries();
7ce961eb	131
401b84d6	132	if (libs.Contains("lib" + detName + "base.so") \|\| (gSystem->Load("lib" + detName + "base.so") >= 0))
	133	pluginManager->AddHandler("AliQAChecker", detName, qacName, detName + "qac", qacName + "()");
	134	else
	135	pluginManager->AddHandler("AliQAChecker", detName, qacName, detName, qacName + "()");
a5fa6165	136
401b84d6	137	pluginHandler = pluginManager->FindHandler("AliQAChecker", detName);
a5fa6165	138
401b84d6	139	if (pluginHandler && (pluginHandler->LoadPlugin() == 0))
401b84d6	140	qac = (AliQACheckerBase *) pluginHandler->ExecPlugin(0);
a5fa6165	141
401b84d6	142	}
7ce961eb	143	}
401b84d6	144	if (qac)
	145	fCheckers[det] = qac ;
	146
	147	return qac ;
a5fa6165	148	}
cbae8032	149
a5fa6165	150	//_____________________________________________________________________________
57acd2d2	151	void AliQAChecker::GetRefSubDir(const char * det, const char * task, TDirectory & dirFile, TObjArray *& dirOCDB)
a5fa6165	152	{
a5fa6165	153	// Opens and returns the file with the reference data
cbae8032	154	dirFile = NULL ;
4e25ac79	155	TString refStorage(AliQAv1::GetQARefStorage()) ;
4e25ac79	156	if (!refStorage.Contains(AliQAv1::GetLabLocalOCDB()) && !refStorage.Contains(AliQAv1::GetLabAliEnOCDB())) {
cf13d879	157	AliError(Form("%s is not a valid location for reference data", refStorage.Data())) ;
	158	return ;
	159	} else {
634696f5	160	AliQAManager* manQA = AliQAManager::QAManager(AliQAv1::GetTaskIndex(task)) ;
958ed851	161	dirOCDB = new TObjArray*[AliRecoParam::kNSpecies] ;
57acd2d2	162	for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
958ed851	163	dirOCDB[specie] = NULL ;
4e25ac79	164	if ( !AliQAv1::Instance()->IsEventSpecieSet(specie) )
57acd2d2	165	continue ;
4e25ac79	166	AliQAv1::SetQARefDataDirName(specie) ;
b03591ab	167	if ( ! manQA->GetLock() ) {
4e25ac79	168	manQA->SetDefaultStorage(AliQAv1::GetQARefStorage()) ;
4e25ac79	169	manQA->SetSpecificStorage("*", AliQAv1::GetQARefStorage()) ;
b03591ab	170	manQA->SetRun(AliCDBManager::Instance()->GetRun()) ;
b03591ab	171	manQA->SetLock() ;
57acd2d2	172	}
4e25ac79	173	char * detOCDBDir = Form("%s/%s/%s", det, AliQAv1::GetRefOCDBDirName(), AliQAv1::GetRefDataDirName()) ;
b03591ab	174	AliCDBEntry * entry = manQA->Get(detOCDBDir, manQA->GetRun()) ;
57acd2d2	175	if (entry) {
eca4fa66	176	TList * listDetQAD =static_cast<TList *>(entry->GetObject()) ;
	177	if ( strcmp(listDetQAD->ClassName(), "TList") != 0 ) {
	178	AliError(Form("Expected a Tlist and found a %s for detector %s", listDetQAD->ClassName(), det)) ;
0e5b0b02	179	listDetQAD = NULL ;
eca4fa66	180	continue ;
0e5b0b02	181	}
	182	if ( listDetQAD ) {
	183	TIter next(listDetQAD) ;
	184	TObjArray * ar ;
	185	while ( (ar = (TObjArray*)next()) )
958ed851	186	dirOCDB[specie] = static_cast<TObjArray *>(listDetQAD->FindObject(Form("%s/%s", task, AliRecoParam::GetEventSpecieName(specie)))) ;
0e5b0b02	187	}
57acd2d2	188	}
cbae8032	189	}
cbae8032	190	}
421ab0fb	191	}
421ab0fb	192
a4976ef3	193	//_____________________________________________________________________________
2e42b4d4	194	AliQAChecker * AliQAChecker::Instance()
a4976ef3	195	{
	196	// returns unique instance of the checker
	197	if ( ! fgQAChecker )
2e42b4d4	198	fgQAChecker = new AliQAChecker() ;
a4976ef3	199	return fgQAChecker ;
	200	}
	201
cbae8032	202	//_____________________________________________________________________________
7e88424f	203	void AliQAChecker::LoadRunInfoFromGRP()
cbae8032	204	{
cbae8032	205	AliCDBManager* man = AliCDBManager::Instance() ;
4e25ac79	206	AliCDBEntry* entry = man->Get(AliQAv1::GetGRPPath().Data());
44e45fac	207	AliGRPObject* grpObject = 0x0;
	208	if (entry) {
	209
eca4fa66	210	TMap* m = static_cast<TMap*>(entry->GetObject()); // old GRP entry
44e45fac	211
44e45fac	212	if (m) {
5379c4a3	213	AliDebug(AliQAv1::GetQADebugLevel(), "It is a map");
44e45fac	214	//m->Print();
	215	grpObject = new AliGRPObject();
	216	grpObject->ReadValuesFromMap(m);
	217	}
	218
	219	else {
5379c4a3	220	AliDebug(AliQAv1::GetQADebugLevel(), "It is a new GRP object");
eca4fa66	221	grpObject = static_cast<AliGRPObject*>(entry->GetObject()); // new GRP entry
44e45fac	222	}
	223
	224	entry->SetOwner(0);
	225	AliCDBManager::Instance()->UnloadFromCache("GRP/GRP/Data");
	226	}
	227
	228	if (!grpObject) {
	229	AliFatal("No GRP entry found in OCDB!");
cbae8032	230	}
44e45fac	231
	232	TString lhcState = grpObject->GetLHCState();
	233	if (lhcState==AliGRPObject::GetInvalidString()) {
	234	AliError("GRP/GRP/Data entry: missing value for the LHC state ! Using UNKNOWN");
	235	lhcState = "UNKNOWN";
	236	}
	237
	238	TString beamType = grpObject->GetBeamType();
	239	if (beamType==AliGRPObject::GetInvalidString()) {
	240	AliError("GRP/GRP/Data entry: missing value for the beam type ! Using UNKNOWN");
	241	beamType = "UNKNOWN";
cbae8032	242	}
44e45fac	243
	244	Float_t beamEnergy = grpObject->GetBeamEnergy();
	245	if (beamEnergy==AliGRPObject::GetInvalidFloat()) {
	246	AliError("GRP/GRP/Data entry: missing value for the beam energy ! Using 0");
	247	beamEnergy = 0;
cbae8032	248	}
44e45fac	249
	250	TString runType = grpObject->GetRunType();
	251	if (runType==AliGRPObject::GetInvalidString()) {
	252	AliError("GRP/GRP/Data entry: missing value for the run type ! Using UNKNOWN");
	253	runType = "UNKNOWN";
7e88424f	254	}
44e45fac	255
	256	Int_t activeDetectors = grpObject->GetDetectorMask();
	257	if (activeDetectors==AliGRPObject::GetInvalidInt()) {
	258	AliError("GRP/GRP/Data entry: missing value for the detector mask ! Using 1074790399");
	259	activeDetectors = 1074790399;
cbae8032	260	}
44e45fac	261
	262	fRunInfo = new AliRunInfo(lhcState, beamType, beamEnergy, runType, activeDetectors);
	263
7e88424f	264	fRunInfoOwner = kTRUE ;
57acd2d2	265
	266	// set the event specie
	267	fEventSpecie = AliRecoParam::kDefault ;
	268	if (strcmp(runType,"PHYSICS")) {
	269	// Not a physics run, the event specie is set to kCalib
	270	fEventSpecie = AliRecoParam::kCalib ;
	271	return;
	272	}
	273	if (strcmp(lhcState,"STABLE_BEAMS") == 0) {
	274	// Heavy ion run (any beam tha is not pp, the event specie is set to kHighMult
	275	fEventSpecie = AliRecoParam::kHighMult ;
	276	if ((strcmp(beamType,"p-p") == 0) \|\|
	277	(strcmp(beamType,"p-") == 0) \|\|
	278	(strcmp(beamType,"-p") == 0) \|\|
	279	(strcmp(beamType,"P-P") == 0) \|\|
	280	(strcmp(beamType,"P-") == 0) \|\|
	281	(strcmp(beamType,"-P") == 0)) {
	282	// Proton run, the event specie is set to kLowMult
	283	fEventSpecie = AliRecoParam::kLowMult ;
	284	}
	285	else if (strcmp(beamType,"-") == 0) {
	286	// No beams, we assume cosmic data
	287	fEventSpecie = AliRecoParam::kCosmic ;
	288	}
	289	else if (strcmp(beamType,"UNKNOWN") == 0) {
	290	// No LHC beam information is available, we use the default event specie
	291	fEventSpecie = AliRecoParam::kDefault ;
	292	}
	293	}
cbae8032	294	}
cbae8032	295
421ab0fb	296	//_____________________________________________________________________________
9c980fe5	297	Bool_t AliQAChecker::Run(const char * fileName, AliDetectorRecoParam * recoParam)
421ab0fb	298	{
44ed7a66	299	// run the Quality Assurance Checker for all tasks Hits, SDigits, Digits, DigitsR, RecPoints, TrackSegments, RecParticles and ESDs
57acd2d2	300	// starting from data in file
421ab0fb	301	TStopwatch stopwatch;
421ab0fb	302	stopwatch.Start();
c39ee44c	303
a5fa6165	304	//search for all detectors QA directories
4e25ac79	305	TList * detKeyList = AliQAv1::GetQADataFile(fileName)->GetListOfKeys() ;
a5fa6165	306	TIter nextd(detKeyList) ;
a5fa6165	307	TKey * detKey ;
eca4fa66	308	while ( (detKey = static_cast<TKey *>(nextd()) ) ) {
5379c4a3	309	AliDebug(AliQAv1::GetQADebugLevel(), Form("Found %s", detKey->GetName())) ;
a5fa6165	310	//Check which detector
	311	TString detName ;
	312	TString detNameQA(detKey->GetName()) ;
	313	Int_t det ;
4e25ac79	314	for ( det = 0; det < AliQAv1::kNDET ; det++) {
4e25ac79	315	detName = AliQAv1::GetDetName(det) ;
a5fa6165	316	if (detNameQA.Contains(detName)) {
57acd2d2	317	fFoundDetectors+=detName ;
	318	fFoundDetectors+="." ;
	319	break ;
a5fa6165	320	}
a5fa6165	321	}
4e25ac79	322	TDirectory * detDir = AliQAv1::GetQADataFile(fileName)->GetDirectory(detKey->GetName()) ;
a5fa6165	323	TList * taskKeyList = detDir->GetListOfKeys() ;
	324	TIter nextt(taskKeyList) ;
	325	TKey * taskKey ;
	326	// now search for the tasks dir
	327	while ( (taskKey = static_cast<TKey *>(nextt()) ) ) {
	328	TString taskName( taskKey->GetName() ) ;
5379c4a3	329	AliDebug(AliQAv1::GetQADebugLevel(), Form("Found %s", taskName.Data())) ;
a5fa6165	330	TDirectory * taskDir = detDir->GetDirectory(taskName.Data()) ;
a5fa6165	331	taskDir->cd() ;
2e42b4d4	332	AliQACheckerBase * qac = GetDetQAChecker(det) ;
a5fa6165	333	if (qac)
5379c4a3	334	AliDebug(AliQAv1::GetQADebugLevel(), Form("QA checker found for %s", detName.Data())) ;
a5fa6165	335	if (!qac)
57acd2d2	336	AliFatal(Form("QA checker not found for %s", detName.Data())) ;
4e25ac79	337	AliQAv1::ALITASK_t index = AliQAv1::kNULLTASK ;
	338	if ( taskName == AliQAv1::GetTaskName(AliQAv1::kHITS) )
	339	index = AliQAv1::kSIM ;
	340	if ( taskName == AliQAv1::GetTaskName(AliQAv1::kSDIGITS) )
	341	index = AliQAv1::kSIM ;
	342	if ( taskName == AliQAv1::GetTaskName(AliQAv1::kDIGITS) )
	343	index = AliQAv1::kSIM ;
d62b07c8	344	if ( taskName == AliQAv1::GetTaskName(AliQAv1::kRAWS) )
d62b07c8	345	index = AliQAv1::kRAW ;
44ed7a66	346	if ( taskName == AliQAv1::GetTaskName(AliQAv1::kDIGITSR) )
44ed7a66	347	index = AliQAv1::kREC ;
4e25ac79	348	if ( taskName == AliQAv1::GetTaskName(AliQAv1::kRECPOINTS) )
	349	index = AliQAv1::kREC ;
	350	if ( taskName == AliQAv1::GetTaskName(AliQAv1::kTRACKSEGMENTS) )
	351	index = AliQAv1::kREC ;
	352	if ( taskName == AliQAv1::GetTaskName(AliQAv1::kRECPARTICLES) )
	353	index = AliQAv1::kREC ;
	354	if ( taskName == AliQAv1::GetTaskName(AliQAv1::kESDS) )
	355	index = AliQAv1::kESD ;
	356	qac->Init(AliQAv1::DETECTORINDEX_t(det)) ;
57acd2d2	357
	358	TDirectory * refDir = NULL ;
	359	TObjArray ** refOCDBDir = NULL ;
	360	GetRefSubDir(detNameQA.Data(), taskName.Data(), refDir, refOCDBDir) ;
4edbc5bc	361	qac->SetRefandData(refDir, refOCDBDir, taskDir) ;
9c980fe5	362	qac->Run(index, recoParam) ;
57acd2d2	363	}
4edbc5bc	364	}
c39ee44c	365	TString detList ;
4e25ac79	366	for ( Int_t det = 0; det < AliQAv1::kNDET; det++) {
2a022cb2	367	if (fFoundDetectors.Contains(AliQAv1::GetDetName(det))) {
2a022cb2	368	detList += AliQAv1::GetDetName(det) ;
c39ee44c	369	detList += " " ;
2a022cb2	370	fFoundDetectors.ReplaceAll(AliQAv1::GetDetName(det), "") ;
2a022cb2	371	AliQAv1::Instance()->Show(AliQAv1::GetDetIndex(AliQAv1::GetDetName(det))) ;
a5fa6165	372	}
421ab0fb	373	}
c39ee44c	374	AliInfo(Form("QA performed for following detectors: %s", detList.Data())) ;
57acd2d2	375	return kTRUE ;
421ab0fb	376	}
421ab0fb	377
a4976ef3	378	//_____________________________________________________________________________
9c980fe5	379	Bool_t AliQAChecker::Run(AliQAv1::DETECTORINDEX_t det, AliQAv1::TASKINDEX_t task, TObjArray ** list, AliDetectorRecoParam * recoParam)
a4976ef3	380	{
d6372ce8	381	// run the Quality Assurance Checker for detector det, for task task starting from data in list
a4976ef3	382
d6372ce8	383	AliQACheckerBase * qac = GetDetQAChecker(det) ;
d6372ce8	384	if (qac)
5379c4a3	385	AliDebug(AliQAv1::GetQADebugLevel(), Form("QA checker found for %s", AliQAv1::GetDetName(det).Data())) ;
d6372ce8	386	if (!qac)
4e25ac79	387	AliError(Form("QA checker not found for %s", AliQAv1::GetDetName(det).Data())) ;
a4976ef3	388
4e25ac79	389	AliQAv1::ALITASK_t index = AliQAv1::kNULLTASK ;
	390	if ( task == AliQAv1::kRAWS )
	391	index = AliQAv1::kRAW ;
	392	else if ( task == AliQAv1::kHITS )
	393	index = AliQAv1::kSIM ;
	394	else if ( task == AliQAv1::kSDIGITS )
	395	index = AliQAv1::kSIM ;
	396	else if ( task == AliQAv1::kDIGITS )
	397	index = AliQAv1::kSIM ;
44ed7a66	398	else if ( task == AliQAv1::kDIGITSR )
44ed7a66	399	index = AliQAv1::kREC ;
4e25ac79	400	else if ( task == AliQAv1::kRECPOINTS )
	401	index = AliQAv1::kREC ;
	402	else if ( task == AliQAv1::kTRACKSEGMENTS )
	403	index = AliQAv1::kREC ;
	404	else if ( task == AliQAv1::kRECPARTICLES )
	405	index = AliQAv1::kREC ;
	406	else if ( task == AliQAv1::kESDS )
	407	index = AliQAv1::kESD ;
4edbc5bc	408
57acd2d2	409	TDirectory * refDir = NULL ;
	410	TObjArray ** refOCDBDir = NULL ;
	411	qac->Init(det) ;
4e25ac79	412	GetRefSubDir(AliQAv1::GetDetName(det), AliQAv1::GetTaskName(task), refDir, refOCDBDir) ;
57acd2d2	413	qac->SetRefandData(refDir, refOCDBDir) ;
9c980fe5	414	qac->Run(index, list, recoParam) ;
634696f5	415
	416	// make the image
	417	qac->MakeImage(list, task, AliQAv1::Mode(task)) ;
	418
d6372ce8	419	return kTRUE ;
31316597	420	}
57acd2d2	421
57acd2d2	422	//_____________________________________________________________________________
9c980fe5	423	Bool_t AliQAChecker::Run(AliQAv1::DETECTORINDEX_t det, AliQAv1::TASKINDEX_t task, TNtupleD ** list, AliDetectorRecoParam * recoParam)
57acd2d2	424	{
	425	// run the Quality Assurance Checker for detector det, for task task starting from data in list
	426
	427	AliQACheckerBase * qac = GetDetQAChecker(det) ;
	428	if (qac)
5379c4a3	429	AliDebug(AliQAv1::GetQADebugLevel(), Form("QA checker found for %s", AliQAv1::GetDetName(det).Data())) ;
57acd2d2	430	if (!qac)
4e25ac79	431	AliError(Form("QA checker not found for %s", AliQAv1::GetDetName(det).Data())) ;
57acd2d2	432
4e25ac79	433	AliQAv1::ALITASK_t index = AliQAv1::kNULLTASK ;
	434	if ( task == AliQAv1::kRAWS )
	435	index = AliQAv1::kRAW ;
	436	else if ( task == AliQAv1::kHITS )
	437	index = AliQAv1::kSIM ;
	438	else if ( task == AliQAv1::kSDIGITS )
	439	index = AliQAv1::kSIM ;
	440	else if ( task == AliQAv1::kDIGITS )
	441	index = AliQAv1::kSIM ;
44ed7a66	442	else if ( task == AliQAv1::kDIGITSR )
44ed7a66	443	index = AliQAv1::kREC ;
4e25ac79	444	else if ( task == AliQAv1::kRECPOINTS )
	445	index = AliQAv1::kREC ;
	446	else if ( task == AliQAv1::kTRACKSEGMENTS )
	447	index = AliQAv1::kREC ;
	448	else if ( task == AliQAv1::kRECPARTICLES )
	449	index = AliQAv1::kREC ;
	450	else if ( task == AliQAv1::kESDS )
	451	index = AliQAv1::kESD ;
57acd2d2	452
	453	TDirectory * refDir = NULL ;
	454	TObjArray ** refOCDBDir = NULL ;
	455	qac->Init(det) ;
4e25ac79	456	GetRefSubDir(AliQAv1::GetDetName(det), AliQAv1::GetTaskName(task), refDir, refOCDBDir) ;
57acd2d2	457	qac->SetRefandData(refDir, refOCDBDir) ;
9c980fe5	458	qac->Run(index, list, recoParam) ;
634696f5	459
57acd2d2	460	return kTRUE ;
57acd2d2	461	}