X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=STEER%2FAliQAChecker.cxx;h=0a109643c3a6b8eb0b967bbc9a4cee31f26b03ff;hb=cd507f9c8fa4e4cd6ce5d6efe7d567d1cb0c5da8;hp=33b8e338bb05a6c83d55497f30945fb6cc71281f;hpb=31316597b5b63998c7282f4ef74da8ee811a3373;p=u%2Fmrichter%2FAliRoot.git diff --git a/STEER/AliQAChecker.cxx b/STEER/AliQAChecker.cxx index 33b8e338bb0..0a109643c3a 100644 --- a/STEER/AliQAChecker.cxx +++ b/STEER/AliQAChecker.cxx @@ -22,16 +22,17 @@ /////////////////////////////////////////////////////////////////////////////// #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 #include @@ -54,10 +55,12 @@ AliQAChecker::AliQAChecker(const char* name, const char* title) : 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 ; } @@ -68,11 +71,13 @@ AliQAChecker::AliQAChecker(const AliQAChecker& qac) : 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 ; } @@ -92,8 +97,9 @@ AliQAChecker::~AliQAChecker() // clean up if (fRunInfo) delete fRunInfo ; - delete [] fCheckers ; - AliQA::Close() ; + for (Int_t det=0; detGetPluginManager() ; TString qacName = "Ali" + detName + "QAChecker" ; @@ -121,7 +127,7 @@ AliQAChecker::~AliQAChecker() 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)) @@ -141,56 +147,6 @@ AliQAChecker::~AliQAChecker() 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(entry->GetObject()) ; - if ( listDetQAD ) - dirOCDB = dynamic_cast(listDetQAD->FindObject(task)) ; - } - } -} //_____________________________________________________________________________ AliQAChecker * AliQAChecker::Instance() @@ -205,173 +161,250 @@ 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(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(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(entry->GetObject()); // new GRP entry + } + + entry->SetOwner(0); + AliCDBManager::Instance()->UnloadFromCache("GRP/GRP/Data"); } - TObjString *lhcState= - dynamic_cast(data->GetValue("fLHCState")); - if (!lhcState) { - AliWarning(Form("%s entry: missing value for the LHC state ! Using UNKNOWN", AliQA::GetGRPPath().Data())); + + if (!grpObject) { + AliFatal("No GRP entry found in OCDB!"); } - TObjString *beamType= - dynamic_cast(data->GetValue("fAliceBeamType")); - if (!beamType) { - AliWarning(Form("%s entry: missing value for the LHC state ! Using UNKNOWN", AliQA::GetGRPPath().Data())); + + 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"; } - TObjString *beamEnergyStr= - dynamic_cast(data->GetValue("fAliceBeamEnergy")); - if (!beamEnergyStr) { - AliWarning(Form("%s entry: missing value for the beam energy ! Using 0", 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 *runType= - dynamic_cast(data->GetValue("fRunType")); - if (!runType) { - AliWarning(Form("%s entry: missing value for the run type ! Using UNKNOWN", AliQA::GetGRPPath().Data())); } - TObjString *activeDetectors= - dynamic_cast(data->GetValue("fDetectorMask")); - if (!activeDetectors) { - AliWarning(Form("%s entry: missing value for the detector mask ! Using 1074790399", 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"; + } + + 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(nextd()) ) ) { - AliDebug(1, Form("Found %s", detKey->GetName())) ; + while ( (detKey = static_cast(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(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(obj)) ; - } else if (className.Contains(TNtupleD::Class()->GetName())) { - qac->Run(index, static_cast(obj)) ; - } else { - AliError(Form("%s class not implemented", className.Data())) ; - return kFALSE ; - } 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 ; + + qac->Init(det) ; + qac->Run(index, list, recoParam) ; + + return kTRUE ; +}