#include <cstdlib>
#include <TClass.h>
#include <TFile.h>
#include <TSystem.h>
#include <TROOT.h>
#include "AliLog.h"
#include "AliQA.h"
ClassImp(AliQA)
AliQA * AliQA::fgQA = 0x0 ;
TFile * AliQA::fgQADataFile = 0x0 ;
TString AliQA::fgQADataFileName = "QA" ;
TFile * AliQA::fgQARefFile = 0x0 ;
TString AliQA::fgQARefDirName = "" ;
TString AliQA::fgQARefFileName = "QA.root" ;
TFile * AliQA::fgQAResultFile = 0x0 ;
TString AliQA::fgQAResultDirName = "" ;
TString AliQA::fgQAResultFileName = "QA.root" ;
TString AliQA::fgDetNames[] = {"ITS", "TPC", "TRD", "TOF", "PHOS", "HMPID", "EMCAL", "MUON", "FMD",
"ZDC", "PMD", "T0", "VZERO", "ACORDE", "HLT", "Global", "CORR","MFT"};
TString AliQA::fgGRPPath = "GRP/GRP/Data" ;
TString AliQA::fgRTNames[] = {"UNKNOWN", "AUTO_TEST", "CALIBRATION", "CALIBRATION_PULSER", "CHANNEL_DELAY_TUNING", "COSMIC",
"COSMICS", "DAQ_FO_UNIF_SCAN", "DAQ_GEN_DAC_SCAN", "DAQ_MEAN_TH_SCAN", "DAQ_MIN_TH_SCAN",
"DAQ_NOISY_PIX_SCAN", "DAQ_PIX_DELAY_SCAN", "DAQ_UNIFORMITY_SCAN", "DCS_FO_UNIF_SCAN",
"DCS_MEAN_TH_SCAN", "DCS_MIN_TH_SCAN", "DCS_PIX_DELAY_SCAN", "DCS_UNIFORMITY_SCAN",
"DDL_TEST", "GAIN", "PEDESTAL", "INJECTOR", "LASER", "MONTECARLO", "NOISE", "NOISY_PIX_SCAN",
"PHYSICS", "PULSER", "STANDALONE", "STANDALONE_BC", "STANDALONE_CENTRAL", "STANDALONE_COSMIC",
"STANDALONE_EMD", "STANDALONE_LASER", "STANDALONE_MB", "STANDALONE_PEDESTAL",
"STANDALONE_SEMICENTRAL", "STANDALONE_PULSER" } ;
TString AliQA::fgTaskNames[] = {"Raws", "Hits", "SDigits", "Digits", "RecPoints", "TrackSegments", "RecParticles", "ESDs"} ;
const TString AliQA::fgkLabLocalFile = "file://" ;
const TString AliQA::fgkLabLocalOCDB = "local://" ;
const TString AliQA::fgkLabAliEnOCDB = "alien://" ;
const TString AliQA::fgkRefFileName = "QA.root" ;
const TString AliQA::fgkQAName = "QA" ;
const TString AliQA::fgkQACorrNtName = "CorrQA" ;
const TString AliQA::fgkRefOCDBDirName = "Ref" ;
TString AliQA::fgkRefDataDirName = "" ;
const TString AliQA::fgkQARefOCDBDefault = "alien://folder=/alice/QA/20" ;
const TString AliQA::fgkExpert = "Expert" ;
const UInt_t AliQA::fgkExpertBit = 16 ;
const UInt_t AliQA::fgkQABit = 17 ;
AliQA::AliQA() :
TNamed("", ""),
fNdet(kNDET),
fQA(new ULong_t[fNdet]),
fDet(kNULLDET),
fTask(kNULLTASK)
{
for (Int_t index = 0 ; index < fNdet ; index++)
fQA[index] = 0 ;
}
AliQA::AliQA(const AliQA& qa) :
TNamed(qa),
fNdet(qa.fNdet),
fQA(new ULong_t[qa.fNdet]),
fDet(qa.fDet),
fTask(qa.fTask)
{
memcpy(fQA, qa.fQA, sizeof(ULong_t)*fNdet);
}
AliQA& AliQA::operator = (const AliQA& qa)
{
if(this != &qa) {
TNamed::operator=(qa);
fNdet = qa.fNdet;
delete [] fQA;
fQA = new ULong_t[fNdet];
memcpy(fQA, qa.fQA, sizeof(ULong_t)*fNdet);
fDet = qa.fDet;
fTask = qa.fTask;
}
return *this;
}
AliQA::AliQA(const DETECTORINDEX_t det) :
TNamed("QA", "Quality Assurance status"),
fNdet(kNDET),
fQA(new ULong_t[fNdet]),
fDet(det),
fTask(kNULLTASK)
{
if (! CheckRange(det) ) {
fDet = kNULLDET ;
return ;
}
Int_t index ;
for (index = 0; index < fNdet; index++)
fQA[index] = 0 ;
}
AliQA::AliQA(const ALITASK_t tsk) :
TNamed("QA", "Quality Assurance status"),
fNdet(kNDET),
fQA(new ULong_t[fNdet]),
fDet(kNULLDET),
fTask(tsk)
{
if (! CheckRange(tsk) ) {
fTask = kNULLTASK ;
return ;
}
Int_t index ;
for (index = 0; index < fNdet; index++)
fQA[index] = 0 ;
}
AliQA::~AliQA()
{
delete[] fQA ;
}
void AliQA::Close()
{
if (fgQADataFile)
if (fgQADataFile->IsOpen())
fgQADataFile->Close() ;
if (fgQAResultFile)
if (fgQAResultFile->IsOpen())
fgQAResultFile->Close() ;
if (fgQARefFile)
if (fgQARefFile->IsOpen())
fgQARefFile->Close() ;
}
Bool_t AliQA::CheckFatal() const
{
Bool_t rv = kFALSE ;
Int_t index ;
for (index = 0; index < kNDET ; index++)
rv = rv || IsSet(DETECTORINDEX_t(index), fTask, kFATAL) ;
return rv ;
}
Bool_t AliQA::CheckRange(DETECTORINDEX_t det) const
{
Bool_t rv = ( det < 0 || det > kNDET ) ? kFALSE : kTRUE ;
if (!rv)
AliFatal(Form("Detector index %d is out of range: 0 <= index <= %d", det, kNDET)) ;
return rv ;
}
Bool_t AliQA::CheckRange(ALITASK_t task) const
{
Bool_t rv = ( task < kRAW || task > kNTASK ) ? kFALSE : kTRUE ;
if (!rv)
AliFatal(Form("Module index %d is out of range: 0 <= index <= %d", task, kNTASK)) ;
return rv ;
}
Bool_t AliQA::CheckRange(QABIT_t bit) const
{
Bool_t rv = ( bit < 0 || bit > kNBIT ) ? kFALSE : kTRUE ;
if (!rv)
AliFatal(Form("Status bit %d is out of range: 0 <= bit <= %d", bit, kNBIT)) ;
return rv ;
}
char * AliQA::GetAliTaskName(ALITASK_t tsk)
{
TString tskName ;
switch (tsk) {
case kNULLTASK:
break ;
case kRAW:
tskName = "RAW" ;
break ;
case kSIM:
tskName = "SIM" ;
break ;
case kREC:
tskName = "REC" ;
break ;
case kESD:
tskName = "ESD" ;
break ;
case kANA:
tskName = "ANA" ;
break ;
default:
tsk = kNULLTASK ;
break ;
}
return Form("%s", tskName.Data()) ;
}
char * AliQA::GetBitName(QABIT_t bit) const
{
TString bitName ;
switch (bit) {
case kNULLBit:
break ;
case kINFO:
bitName = "INFO" ;
break ;
case kWARNING:
bitName = "WARNING" ;
break ;
case kERROR:
bitName = "ERROR" ;
break ;
case kFATAL:
bitName = "FATAL" ;
break ;
default:
bit = kNULLBit ;
break ;
}
return Form("%s", bitName.Data()) ;
}
AliQA::DETECTORINDEX_t AliQA::GetDetIndex(const char * name)
{
TString sname(name) ;
DETECTORINDEX_t rv = kNULLDET ;
for (Int_t det = 0; det < kNDET ; det++) {
if ( GetDetName(det) == sname ) {
rv = DETECTORINDEX_t(det) ;
break ;
}
}
return rv ;
}
const char * AliQA::GetDetName(Int_t det)
{
if ( det >= 0 && det < kNDET)
return (fgDetNames[det]).Data() ;
else
return NULL ;
}
TFile * AliQA::GetQADataFile(const char * name, const Int_t run)
{
const char * temp = Form("%s.%s.%d.root", name, fgQADataFileName.Data(), run) ;
TString opt ;
if (! fgQADataFile ) {
if (gSystem->AccessPathName(temp))
opt = "NEW" ;
else
opt = "UPDATE" ;
fgQADataFile = TFile::Open(temp, opt.Data()) ;
} else {
if ( strcmp(temp, fgQADataFile->GetName()) != 0 ) {
fgQADataFile = dynamic_cast<TFile *>(gROOT->FindObject(temp)) ;
if ( !fgQADataFile ) {
if (gSystem->AccessPathName(temp))
opt = "NEW" ;
else
opt = "UPDATE" ;
fgQADataFile = TFile::Open(temp, opt.Data()) ;
}
}
}
return fgQADataFile ;
}
TFile * AliQA::GetQADataFile(const char * fileName)
{
if (!fgQADataFile)
if (!fileName)
fileName = AliQA::GetQADataFileName() ;
if (!gSystem->AccessPathName(fileName)) {
fgQADataFile = TFile::Open(fileName) ;
} else {
AliFatalClass(Form("File %s not found", fileName)) ;
}
return fgQADataFile ;
}
TFile * AliQA::GetQAResultFile()
{
if (fgQAResultFile)
fgQAResultFile->Close() ;
fgQAResultFile = 0x0 ;
TString dirName(fgQAResultDirName) ;
if ( dirName.Contains(fgkLabLocalFile))
dirName.ReplaceAll(fgkLabLocalFile, "") ;
TString fileName(dirName + fgQAResultFileName) ;
TString opt("") ;
if ( !gSystem->AccessPathName(fileName) )
opt = "UPDATE" ;
else {
if ( gSystem->AccessPathName(dirName) )
gSystem->mkdir(dirName) ;
opt = "NEW" ;
}
fgQAResultFile = TFile::Open(fileName, opt) ;
return fgQAResultFile ;
}
const TString AliQA::GetRunTypeName(RUNTYPE_t rt)
{
TString rv("Invalid Run Type") ;
if ( rt == kNULLTYPE ) {
rv = "Known RUN_TYPE are: \n" ;
for (Int_t index = 0 ; index < kNTYPE; index++) {
rv += Form("%2d -- %s\n", index, fgRTNames[index].Data()) ;
}
AliErrorClass(Form("%s", rv.Data())) ;
return "" ;
}
else {
if ( rt > kNULLTYPE && rt < kNTYPE )
rv = fgRTNames[rt] ;
}
return rv ;
}
AliQA::TASKINDEX_t AliQA::GetTaskIndex(const char * name)
{
TString sname(name) ;
TASKINDEX_t rv = kNULLTASKINDEX ;
for (Int_t tsk = 0; tsk < kNTASKINDEX ; tsk++) {
if ( GetTaskName(tsk) == sname ) {
rv = TASKINDEX_t(tsk) ;
break ;
}
}
return rv ;
}
Bool_t AliQA::IsSet(DETECTORINDEX_t det, ALITASK_t tsk, QABIT_t bit) const
{
CheckRange(det) ;
CheckRange(tsk) ;
CheckRange(bit) ;
ULong_t offset = Offset(tsk) ;
ULong_t status = GetStatus(det) ;
offset+= bit ;
status = (status & 1 << offset) != 0 ;
return status ;
}
Bool_t AliQA::IsSetAny(DETECTORINDEX_t det, ALITASK_t tsk) const
{
CheckRange(det) ;
CheckRange(tsk) ;
ULong_t offset = Offset(tsk) ;
ULong_t status = GetStatus(det) ;
UShort_t st = 0 ;
for ( Int_t bit = 0 ; bit < kNBIT ; bit++) {
offset+= bit ;
st += (status & 1 << offset) != 0 ;
}
if ( st == 0 )
return kFALSE ;
else
return kTRUE ;
}
Bool_t AliQA::IsSetAny(DETECTORINDEX_t det) const
{
CheckRange(det) ;
ULong_t status = GetStatus(det) ;
UShort_t st = 0 ;
for ( Int_t tsk = 0 ; tsk < kNTASK ; tsk++) {
ULong_t offset = Offset(ALITASK_t(tsk)) ;
for ( Int_t bit = 0 ; bit < kNBIT ; bit++) {
offset+= bit ;
st += (status & 1 << offset) != 0 ;
}
}
if ( st == 0 )
return kFALSE ;
else
return kTRUE ;
}
AliQA * AliQA::Instance()
{
return fgQA ;
}
AliQA * AliQA::Instance(const DETECTORINDEX_t det)
{
if ( ! fgQA) {
TFile * f = GetQAResultFile() ;
fgQA = dynamic_cast<AliQA *>(f->Get("QA")) ;
if ( ! fgQA )
fgQA = new AliQA(det) ;
}
fgQA->Set(det) ;
return fgQA ;
}
AliQA * AliQA::Instance(const ALITASK_t tsk)
{
if ( ! fgQA)
switch (tsk) {
case kNULLTASK:
break ;
case kRAW:
fgQA = new AliQA(tsk) ;
break ;
case kSIM:
fgQA = new AliQA(tsk) ;
break ;
case kREC:
AliInfoClass("fgQA = gAlice->GetQA()") ;
break ;
case kESD:
AliInfoClass("fgQA = dynamic_cast<AliQA *> (esdFile->Get(\"QA\")") ;
break ;
case kANA:
AliInfoClass("fgQA = dynamic_cast<AliQA *> (esdFile->Get(\"QA\")") ;
break ;
case kNTASK:
break ;
}
if (fgQA)
fgQA->Set(tsk) ;
return fgQA ;
}
AliQA * AliQA::Instance(const TASKINDEX_t tsk)
{
ALITASK_t index = kNULLTASK ;
if ( tsk == kRAWS )
index = kRAW ;
else if (tsk < kDIGITS)
index = kSIM ;
else if (tsk < kRECPARTICLES)
index = kREC ;
else if (tsk == kESDS)
index = kESD ;
return Instance(index) ;
}
Long64_t AliQA::Merge(const TCollection * list) {
Long64_t nmerge=0;
for (Int_t det = 0 ; det < kNDET ; det++) {
Set(DETECTORINDEX_t(det)) ;
for (Int_t task = 0 ; task < kNTASK ; task++) {
Set(ALITASK_t(task)) ;
for (Int_t bit = 0 ; bit < kNBIT ; bit++) {
TIter next(list) ;
AliQA * qa ;
while ( (qa = (AliQA*)next() ) ) {
++nmerge;
qa->IsSet(DETECTORINDEX_t(det), ALITASK_t(task), QABIT_t(bit)) ;
Set(QABIT_t(bit)) ;
}
}
}
}
return nmerge;
}
ULong_t AliQA::Offset(ALITASK_t tsk) const
{
CheckRange(tsk) ;
ULong_t offset = 0 ;
switch (tsk) {
case kNULLTASK:
break ;
case kRAW:
offset+= 0 ;
break ;
case kSIM:
offset+= 4 ;
break ;
case kREC:
offset+= 8 ;
break ;
case kESD:
offset+= 12 ;
break ;
case kANA:
offset+= 16 ;
break ;
case kNTASK:
break ;
}
return offset ;
}
void AliQA::Set(QABIT_t bit)
{
SetStatusBit(fDet, fTask, bit) ;
}
void AliQA::SetQARefStorage(const char * name)
{
fgQARefDirName = name ;
if ( fgQARefDirName.Contains(fgkLabLocalFile) )
fgQARefFileName = fgkRefFileName ;
else if ( fgQARefDirName.Contains(fgkLabLocalOCDB) )
fgQARefFileName = fgkQAName ;
else if ( fgQARefDirName.Contains(fgkLabAliEnOCDB) )
fgQARefFileName = fgkQAName ;
else {
AliErrorClass(Form("%s is an invalid storage definition\n", name)) ;
fgQARefDirName = "" ;
fgQARefFileName = "" ;
}
TString tmp(fgQARefDirName) ;
AliInfoClass(Form("AliQA::SetQARefDir: QA references are in %s\n", tmp.Data()) ) ;
}
void AliQA::SetQARefDataDirName(const char * name)
{
TString test(name) ;
RUNTYPE_t rt = kNULLTYPE ;
for (Int_t index = 0; index < kNTYPE; index++) {
if (test.CompareTo(fgRTNames[index]) == 0) {
rt = (RUNTYPE_t) index ;
break ;
}
}
if (rt == kNULLTYPE) {
AliWarningClass(Form("AliQA::SetQARefDataDirName: %s is an unknown RUN TYPE name\n", name)) ;
return ;
}
SetQARefDataDirName(rt) ;
}
void AliQA::SetQAResultDirName(const char * name)
{
fgQAResultDirName.Prepend(name) ;
AliInfoClass(Form("AliQA::SetQAResultDirName: QA results are in %s\n", fgQAResultDirName.Data())) ;
if ( fgQAResultDirName.Contains(fgkLabLocalFile))
fgQAResultDirName.ReplaceAll(fgkLabLocalFile, "") ;
fgQAResultFileName.Prepend(fgQAResultDirName) ;
}
void AliQA::SetStatusBit(DETECTORINDEX_t det, ALITASK_t tsk, QABIT_t bit)
{
CheckRange(det) ;
CheckRange(tsk) ;
CheckRange(bit) ;
ULong_t offset = Offset(tsk) ;
ULong_t status = GetStatus(det) ;
offset+= bit ;
status = status | 1 << offset ;
SetStatus(det, status) ;
}
void AliQA::ShowAll() const
{
Int_t index ;
for (index = 0 ; index < kNDET ; index++) {
for (Int_t tsk = kRAW ; tsk < kNTASK ; tsk++) {
ShowStatus(DETECTORINDEX_t(index), ALITASK_t(tsk)) ;
}
}
}
void AliQA::ShowStatus(DETECTORINDEX_t det, ALITASK_t tsk) const
{
CheckRange(det) ;
ULong_t status = GetStatus(det) ;
ULong_t tskStatus[kNTASK] ;
tskStatus[kRAW] = status & 0x0000f ;
tskStatus[kSIM] = status & 0x000f0 ;
tskStatus[kREC] = status & 0x00f00 ;
tskStatus[kESD] = status & 0x0f000 ;
tskStatus[kANA] = status & 0xf0000 ;
AliInfo(Form("====> QA Status for %8s raw =0x%lx, sim=0x%lx, rec=0x%lx, esd=0x%lx, ana=0x%lx", GetDetName(det).Data(),
tskStatus[kRAW], tskStatus[kSIM], tskStatus[kREC], tskStatus[kESD], tskStatus[kANA] )) ;
if (tsk == kNULLTASK) {
for (Int_t itsk = kRAW ; itsk < kNTASK ; itsk++) {
ShowASCIIStatus(det, ALITASK_t(itsk), tskStatus[itsk]) ;
}
} else {
ShowASCIIStatus(det, tsk, tskStatus[tsk]) ;
}
}
void AliQA::ShowASCIIStatus(DETECTORINDEX_t det, ALITASK_t tsk, const ULong_t status) const
{
TString text;
for (Int_t bit = kINFO ; bit < kNBIT ; bit++) {
if (IsSet(det, tsk, QABIT_t(bit))) {
text = GetBitName(QABIT_t(bit)) ;
text += " " ;
}
}
if (! text.IsNull())
AliWarningClass(Form(" %8s %4s 0x%4lx, Problem signalled: %8s \n", GetDetName(det).Data(), GetAliTaskName(tsk), status, text.Data())) ;
}
void AliQA::UnSet(QABIT_t bit)
{
UnSetStatusBit(fDet, fTask, bit) ;
}
void AliQA::UnSetStatusBit(DETECTORINDEX_t det, ALITASK_t tsk, QABIT_t bit)
{
CheckRange(det) ;
CheckRange(tsk) ;
CheckRange(bit) ;
ULong_t offset = Offset(tsk) ;
ULong_t status = GetStatus(det) ;
offset+= bit ;
status = status & 0 << offset ;
SetStatus(det, status) ;
}