///////////////////////////////////////////////////////////////////////////////
#include "AliCDBEntry.h"
-#include "AliCDBManager.h"
+#include "AliQAManager.h"
#include "AliCDBStorage.h"
-#include "AliRunInfo.h"
+#include "AliRunInfo.h"
#include "AliLog.h"
#include "AliModule.h"
-#include "AliQA.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>
fRunInfo(0x0),
fRunInfoOwner(kFALSE),
fRefFile(0x0),
- fFoundDetectors(".")
+ fFoundDetectors("."),
+ fEventSpecie(AliRecoParam::kDefault),
+ fRun(0)
{
// ctor: initialise checkers and open the data file
- for (Int_t det = 0 ; det < AliQA::kNDET ; det++)
+ for (Int_t det = 0 ; det < AliQAv1::kNDET ; det++)
fCheckers[det] = NULL ;
}
fRunInfo(qac.fRunInfo),
fRunInfoOwner(kFALSE),
fRefFile(qac.fRefFile),
- fFoundDetectors(qac.fFoundDetectors)
+ fFoundDetectors(qac.fFoundDetectors),
+ fEventSpecie(qac.fEventSpecie),
+ fRun(qac.fRun)
{
// copy constructor
- for (Int_t det = 0 ; det < AliQA::kNDET ; det++)
+ for (Int_t det = 0 ; det < AliQAv1::kNDET ; det++)
fCheckers[det] = NULL ;
}
// clean up
if (fRunInfo)
delete fRunInfo ;
- delete [] fCheckers ;
- AliQA::Close() ;
+ for (Int_t det=0; det<AliQAv1::kNDET; det++)
+ delete fCheckers[det] ;
+ AliQAv1::Close() ;
}
//_____________________________________________________________________________
AliQACheckerBase * qac = NULL ;
- TString detName(AliQA::GetDetName(det)) ;
+ TString detName(AliQAv1::GetDetName(det)) ;
- if (det == AliQA::kGLOBAL) {
+ if (det == AliQAv1::kGLOBAL) {
qac = new AliGlobalQAChecker() ;
- } else if (det == AliQA::kCORR) {
+ } else if (det == AliQAv1::kCORR) {
qac = new AliCorrQAChecker() ;
} else {
- AliDebug(1, Form("Retrieving QA checker for %s", detName.Data())) ;
+ AliDebugClass(AliQAv1::GetQADebugLevel(), Form("Retrieving QA checker for %s", detName.Data())) ;
TPluginManager* pluginManager = gROOT->GetPluginManager() ;
TString qacName = "Ali" + detName + "QAChecker" ;
TPluginHandler* pluginHandler = pluginManager->FindHandler("AliQAChecker", detName.Data());
// if not, add a plugin for it
if (!pluginHandler) {
- //AliInfo(Form("defining plugin for %s", qacName.Data()));
+ //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))
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 ;
- dirOCDB = NULL ;
- TString refStorage(AliQA::GetQARefStorage()) ;
- //refStorage += AliQA::GetQARefFileName() ;
- if (refStorage.Contains(AliQA::GetLabLocalFile())) {
- refStorage.ReplaceAll(AliQA::GetLabLocalFile(), "") ;
- if ( fRefFile )
- if ( fRefFile->IsOpen() )
- fRefFile->Close() ;
- fRefFile = TFile::Open(refStorage.Data()) ;
- if (!fRefFile) {
- AliError(Form("Cannot find reference file %s", refStorage.Data())) ;
- dirFile = NULL ;
- }
- dirFile = fRefFile->GetDirectory(det) ;
- if (!dirFile) {
- AliWarning(Form("Directory %s not found in %d", det, refStorage.Data())) ;
- } else {
- dirFile = dirFile->GetDirectory(task) ;
- if (!dirFile)
- AliWarning(Form("Directory %s/%s not found in %s", det, task, refStorage.Data())) ;
- }
- } else if (refStorage.Contains(AliQA::GetLabLocalOCDB()) || refStorage.Contains(AliQA::GetLabAliEnOCDB())) {
- AliCDBManager* man = AliCDBManager::Instance() ;
- if ( strcmp(AliQA::GetRefDataDirName(), "") == 0 ) { // the name of the last level of the directory is not set (RUNTYPE)
- // Get it from RunInfo
- if (!fRunInfo) // not yet set, get the info from GRP
- LoadRunInfoFromGRP() ;
- AliQA::SetQARefDataDirName(fRunInfo->GetRunType()) ;
- }
- if ( ! man->GetLock() ) {
- man->SetDefaultStorage(AliQA::GetQARefStorage()) ;
- man->SetSpecificStorage("*", AliQA::GetQARefStorage()) ;
- }
- char * detOCDBDir = Form("%s/%s/%s", det, AliQA::GetRefOCDBDirName(), AliQA::GetRefDataDirName()) ;
- AliInfo(Form("Reference QA data are taken from %s", detOCDBDir)) ;
- AliCDBEntry * entry = man->Get(detOCDBDir, man->GetRun()) ;
- if (entry) {
- TList * listDetQAD = dynamic_cast<TList *>(entry->GetObject()) ;
- if ( listDetQAD )
- dirOCDB = dynamic_cast<TObjArray *>(listDetQAD->FindObject(task)) ;
- }
- }
-}
//_____________________________________________________________________________
AliQAChecker * AliQAChecker::Instance()
void AliQAChecker::LoadRunInfoFromGRP()
{
AliCDBManager* man = AliCDBManager::Instance() ;
- AliCDBEntry* entry = man->Get(AliQA::GetGRPPath().Data());
- TMap * data = 0x0 ;
- if (entry)
- data = dynamic_cast<TMap*>(entry->GetObject());
- if (!data) {
- AliFatal("No GRP entry found in OCDB!");
+ 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";
}
- TObjString *lhcState=
- dynamic_cast<TObjString*>(data->GetValue("fLHCState"));
- if (!lhcState) {
- AliWarning(Form("%s entry: missing value for the LHC state ! Using UNKNOWN", AliQA::GetGRPPath().Data()));
+
+ TString beamType = grpObject->GetBeamType();
+ if (beamType==AliGRPObject::GetInvalidString()) {
+ AliError("GRP/GRP/Data entry: missing value for the beam type ! Using UNKNOWN");
+ beamType = "UNKNOWN";
}
- TObjString *beamType=
- dynamic_cast<TObjString*>(data->GetValue("fAliceBeamType"));
- if (!beamType) {
- AliWarning(Form("%s entry: missing value for the LHC state ! Using UNKNOWN", AliQA::GetGRPPath().Data()));
+
+ Float_t beamEnergy = grpObject->GetBeamEnergy();
+ if (beamEnergy==AliGRPObject::GetInvalidFloat()) {
+ AliError("GRP/GRP/Data entry: missing value for the beam energy ! Using 0");
+ beamEnergy = 0;
}
- TObjString *beamEnergyStr=
- dynamic_cast<TObjString*>(data->GetValue("fAliceBeamEnergy"));
- if (!beamEnergyStr) {
- AliWarning(Form("%s entry: missing value for the beam energy ! Using 0", AliQA::GetGRPPath().Data()));
+
+ TString runType = grpObject->GetRunType();
+ if (runType==AliGRPObject::GetInvalidString()) {
+ AliError("GRP/GRP/Data entry: missing value for the run type ! Using UNKNOWN");
+ runType = "UNKNOWN";
}
- TObjString *runType=
- dynamic_cast<TObjString*>(data->GetValue("fRunType"));
- if (!runType) {
- AliWarning(Form("%s entry: missing value for the run type ! Using UNKNOWN", AliQA::GetGRPPath().Data())); }
- TObjString *activeDetectors=
- dynamic_cast<TObjString*>(data->GetValue("fDetectorMask"));
- if (!activeDetectors) {
- AliWarning(Form("%s entry: missing value for the detector mask ! Using 1074790399", AliQA::GetGRPPath().Data()));
+
+ 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 ? lhcState->GetString().Data() : "UNKNOWN",
- beamType ? beamType->GetString().Data() : "UNKNOWN",
- beamEnergyStr ? beamEnergyStr->GetString().Atof() : 0,
- runType ? runType->GetString().Data() : "UNKNOWN",
- activeDetectors ? activeDetectors->GetString().Atoi() : 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)
+Bool_t AliQAChecker::Run(const char * fileName, AliDetectorRecoParam * recoParam)
{
- // run the Quality Assurance Checker for all tasks Hits, SDigits, Digits, recpoints, tracksegments, recparticles and ESDs
- // starting from data in file
-
- Bool_t rv = kFALSE ;
-
+ // 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 = AliQA::GetQADataFile(fileName)->GetListOfKeys() ;
+ TList * detKeyList = AliQAv1::GetQADataFile(fileName)->GetListOfKeys() ;
TIter nextd(detKeyList) ;
TKey * detKey ;
- while ( (detKey = dynamic_cast<TKey *>(nextd()) ) ) {
- AliDebug(1, Form("Found %s", detKey->GetName())) ;
+ 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 < AliQA::kNDET ; det++) {
- detName = AliQA::GetDetName(det) ;
+ for ( det = 0; det < AliQAv1::kNDET ; det++) {
+ detName = AliQAv1::GetDetName(det) ;
if (detNameQA.Contains(detName)) {
- fFoundDetectors+=detName ;
- fFoundDetectors+="." ;
- break ;
+ fFoundDetectors+=detName ;
+ fFoundDetectors+="." ;
+ break ;
}
}
- TDirectory * detDir = AliQA::GetQADataFile(fileName)->GetDirectory(detKey->GetName()) ;
+ 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() ) ;
- AliInfo(Form("Found %s", taskName.Data())) ;
+ AliDebug(AliQAv1::GetQADebugLevel(), Form("Found %s", taskName.Data())) ;
TDirectory * taskDir = detDir->GetDirectory(taskName.Data()) ;
taskDir->cd() ;
AliQACheckerBase * qac = GetDetQAChecker(det) ;
if (qac)
- AliInfo(Form("QA checker found for %s", detName.Data())) ;
+ AliDebug(AliQAv1::GetQADebugLevel(), Form("QA checker found for %s", detName.Data())) ;
if (!qac)
- AliFatal(Form("QA checker not found for %s", detName.Data())) ;
- AliQA::ALITASK_t index = AliQA::kNULLTASK ;
- if ( taskName == AliQA::GetTaskName(AliQA::kHITS) )
- index = AliQA::kSIM ;
- if ( taskName == AliQA::GetTaskName(AliQA::kSDIGITS) )
- index = AliQA::kSIM ;
- if ( taskName == AliQA::GetTaskName(AliQA::kDIGITS) )
- index = AliQA::kSIM ;
- if ( taskName == AliQA::GetTaskName(AliQA::kRECPOINTS) )
- index = AliQA::kREC ;
- if ( taskName == AliQA::GetTaskName(AliQA::kTRACKSEGMENTS) )
- index = AliQA::kREC ;
- if ( taskName == AliQA::GetTaskName(AliQA::kRECPARTICLES) )
- index = AliQA::kREC ;
- if ( taskName == AliQA::GetTaskName(AliQA::kESDS) )
- index = AliQA::kESD ;
- qac->Init(AliQA::DETECTORINDEX_t(det)) ;
-
- TDirectory * refDir = NULL ;
- TObjArray * refOCDBDir = NULL ;
- GetRefSubDir(detNameQA.Data(), taskName.Data(), refDir, refOCDBDir) ;
- if ( refDir || refOCDBDir) {
- qac->SetRefandData(refDir, refOCDBDir, taskDir) ;
- qac->Run(index) ;
- }
- }
+ 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)) ;
+ qac->Run(index, recoParam) ;
+ }
}
- AliInfo("QA performed for following detectors:") ;
- for ( Int_t det = 0; det < AliQA::kNDET; det++) {
- if (fFoundDetectors.Contains(AliQA::GetDetName(det))) {
- printf("%s, ",AliQA::GetDetName(det)) ;
- fFoundDetectors.ReplaceAll(AliQA::GetDetName(det), "") ;
+ 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))) ;
}
}
- printf("\n") ;
- rv = kTRUE ;
-
- return rv ;
-
+ AliInfo(Form("QA performed for following detectors: %s", detList.Data())) ;
+ return kTRUE ;
}
//_____________________________________________________________________________
-Bool_t AliQAChecker::Run(AliQA::DETECTORINDEX_t det, AliQA::TASKINDEX_t task, TObject * obj)
+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
+ if (det >= AliQAv1::kNDET) {
+ AliError(Form("det = %i is larger than AliQAv1::kNDET ... should never happen", det));
+ return kFALSE ;
+ }
AliQACheckerBase * qac = GetDetQAChecker(det) ;
if (qac)
- AliDebug(1, Form("QA checker found for %s", AliQA::GetDetName(det).Data())) ;
+ AliDebug(AliQAv1::GetQADebugLevel(), Form("QA checker found for %s", AliQAv1::GetDetName(det).Data())) ;
if (!qac)
- AliError(Form("QA checker not found for %s", AliQA::GetDetName(det).Data())) ;
+ 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 ;
+
+ qac->Init(det) ;
+ qac->Run(index, list, recoParam) ;
+
+ // make the image
+ qac->MakeImage(list, task, AliQAv1::Mode(task)) ;
- AliQA::ALITASK_t index = AliQA::kNULLTASK ;
- if ( task == AliQA::kRAWS )
- index = AliQA::kRAW ;
- else if ( task == AliQA::kHITS )
- index = AliQA::kSIM ;
- else if ( task == AliQA::kSDIGITS )
- index = AliQA::kSIM ;
- else if ( task == AliQA::kDIGITS )
- index = AliQA::kSIM ;
- else if ( task == AliQA::kRECPOINTS )
- index = AliQA::kREC ;
- else if ( task == AliQA::kTRACKSEGMENTS )
- index = AliQA::kREC ;
- else if ( task == AliQA::kRECPARTICLES )
- index = AliQA::kREC ;
- else if ( task == AliQA::kESDS )
- index = AliQA::kESD ;
-
- TDirectory * refDir = NULL ;
- TObjArray * refOCDBDir = NULL ;
- qac->Init(det) ;
- GetRefSubDir(AliQA::GetDetName(det), AliQA::GetTaskName(task), refDir, refOCDBDir) ;
- if ( refDir || refOCDBDir) // references found
- qac->SetRefandData(refDir, refOCDBDir) ;
-
- TString className(obj->ClassName()) ;
- if (className.Contains(TObjArray::Class()->GetName())) {
- qac->Run(index, static_cast<TObjArray *>(obj)) ;
- } else if (className.Contains(TNtupleD::Class()->GetName())) {
- qac->Run(index, static_cast<TNtupleD *>(obj)) ;
- } else {
- AliError(Form("%s class not implemented", className.Data())) ;
- return kFALSE ;
- }
return kTRUE ;
-}
\ No newline at end of file
+}
+
+//_____________________________________________________________________________
+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 ;
+
+ qac->Init(det) ;
+ qac->Run(index, list, recoParam) ;
+
+ return kTRUE ;
+}