AliACORDEQAChecker::AliACORDEQAChecker() :
AliQACheckerBase("ACORDE","ACORDE Quality Assurance Data Maker"),
-fTextDQMShifterInfo(new TPaveText(2.2,45,2.9,50,"T")),
-fMax(new TLine(1,30,3,30))
+fTextDQMShifterInfo(new TPaveText(2.2,0.90,2.9,0.99,"T")),
+fMax(new TLine(1,0.5,3,0.5))
{
// default constructor
fMax->SetLineColor(kGreen);
//____________________________________________________________________________
AliACORDEQAChecker::AliACORDEQAChecker(const AliACORDEQAChecker& qac) :
AliQACheckerBase(qac.GetName(), qac.GetTitle()),
-fTextDQMShifterInfo(new TPaveText(2.2,45,2.9,50,"T")),
+fTextDQMShifterInfo(new TPaveText(2.2,0.90,2.9,0.99,"T")),
fMax(static_cast<TLine*>(qac.fMax->Clone()))
{
//
for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
{
+ test[specie] = 1.0;
if ( !AliQAv1::Instance()->IsEventSpecieSet(specie) ) continue ;
if (list[specie]->GetEntries() == 0) test[specie] = 1.; // Nothing to check
else
TH1 * hdata ; // Data created by the AliACORDEQADataMakerXXX (Sim/Rec)
while ( (hdata = dynamic_cast<TH1 *>(next())) )
{
- if (hdata)
+ //if (hdata)
+ if (hdata->GetEntries()!=0)
{
Double_t rv = 0.0 ;
if(hdata->GetEntries()>0) rv=1;
AliDebug(AliQAv1::GetQADebugLevel(), Form("%s -> %f", hdata->GetName(), rv)) ;
TString hdataName = hdata->GetName();
//if (hdata->GetListOfFunctions()->GetEntries() == 0) continue;
- if (hdataName.Contains("fhACORDEStatusAMU_DQM")){
+
+ if (hdataName.Contains("fhACORDEStatusAMU_DQM"))
+ {
if (!hdata->GetListOfFunctions()->Contains(fTextDQMShifterInfo))
- hdata->GetListOfFunctions()->Add(fTextDQMShifterInfo);
+ hdata->GetListOfFunctions()->Add(fTextDQMShifterInfo);
if (!hdata->GetListOfFunctions()->Contains(fMax))
- hdata->GetListOfFunctions()->Add(fMax);
-
+ hdata->GetListOfFunctions()->Add(fMax);
}
- // Here we use the QAref ACORDE data from fRefOCDBSubDir
-
+ // check with OCDB ref. data
if ( (fRefOCDBSubDir[specie]) && (hdataName.Contains("ACORDEBitPattern")) )
{
TH1 * href = NULL;
if (fRefSubDir) href = static_cast<TH1*>(fRefSubDir->Get(hdata->GetName()));
else if (fRefOCDBSubDir[specie]) href = static_cast<TH1*>(fRefOCDBSubDir[specie]->FindObject(hdata->GetName()));
test[specie] = CheckAcordeRefHits(href,hdata);
- }else if (hdataName.Contains("ACORDEBitPattern")||hdataName.Contains("fhACORDEStatusAMU_DQM"))
- // Here we use an inner QA Checher without the QAref data
+ }else
{
- Float_t acoDataMax = hdata->GetMaximum();
- Int_t flagAcoQAChecker = 0;
- if (acoDataMax == 0) continue;
- for(Int_t i=0;i<60;i++)
+ if (hdataName.Contains("ACORDEBitPattern"))
+ // Here we use an inner QA Checher without the QAref data
{
- if ((hdata->GetBinContent(i)/acoDataMax) < 0.75) flagAcoQAChecker++;
- // if (hdataName.Contains("fhACORDEStatusSL0_DQM")||hdataName.Contains("fhACORDEStatusAMU_DQM")){
- // if (hdata->GetBinContent(i)<0.29) flagAcoQA++;
- // }
- }
- if (hdataName.Contains("ACORDEBitPattern")){
+ Float_t acoDataMax = hdata->GetMaximum();
+ Int_t flagAcoQAChecker = 0;
+ if (acoDataMax == 0) continue;
+ for(Int_t i=0;i<60;i++)
+ {
+ if ((hdata->GetBinContent(i)/acoDataMax) < 0.75) flagAcoQAChecker++;
+ }
Double_t simpleFlag = 1.-flagAcoQAChecker/60.;
if ( (simpleFlag >= 0.90) && (simpleFlag <= 1.0) ) test[specie] = 0.75; // INFO
if ( (simpleFlag >= 0.70) && (simpleFlag < 0.90) ) test[specie] = 0.50; // WARNING
if ( (simpleFlag >= 0.25) && (simpleFlag < 0.70) ) test[specie] = 0.25; // ERROR
if ( (simpleFlag >= 0.0) && (simpleFlag < 0.25) ) test[specie] = -1.0; // FATAL
+ if (test[specie]>0) rv = test[specie];
+ else rv = 1.0;
}
- /*if (hdataName.Contains("fhACORDEStatusSL0_DQM"))
- {
- if (hdata->GetBinContent(1) != 0)
- {
- if (hdata->GetBinContent(2) < 31)
- test[specie] = 0.75;
- else test[specie] = 0.3;
- } else continue;
- }*/
if (hdataName.Contains("fhACORDEStatusAMU_DQM"))
{
- if (hdata->GetBinContent(1) != 0)
- {
- if (hdata->GetBinContent(2) < 31)
+ Float_t goodModules = hdata->GetBinContent(1);
+ Float_t badModules = hdata->GetBinContent(2);
+ if (goodModules >= badModules)
test[specie] = 0.75;
- else test[specie] = 0.3;
- } else continue;
- }
-
- }
-
- // Setting Warning message for possible Empty Events with the ACORDE-Trigger
-
- if (hdataName.Contains("ACORDEMultiplicity") && (hdata->GetBinContent(0)!=0)) AliWarning("Empty event registered with ACORDE Trigger !!!");
-
-
- }else AliError("Data type cannot be processed") ;
- if (fTextDQMShifterInfo)
- {
- fTextDQMShifterInfo->Clear();
- if (test[specie] > 0.3){
- fTextDQMShifterInfo->SetFillColor(kGreen);
- fTextDQMShifterInfo->AddText("ACORDE: O.K.");
- }else{
- fTextDQMShifterInfo->SetFillColor(kRed);
- fTextDQMShifterInfo->AddText("ACORDE: Not, O.K.");
- }
- }
-
- }
- }
+ else test[specie] = 0.30;
+
+ } // check for DQM condition
+ rv = test[specie];
+ if (fTextDQMShifterInfo) // DQM-flag asignment
+ {
+ fTextDQMShifterInfo->Clear();
+ if (rv > 0.4){
+ fTextDQMShifterInfo->SetFillColor(kGreen);
+ fTextDQMShifterInfo->AddText("ACORDE: O.K.");
+ }else{
+ fTextDQMShifterInfo->SetFillColor(kRed);
+ fTextDQMShifterInfo->AddText("ACORDE: Not, O.K.");
+ }
+ } // end DQM flag asignment
+ }// NO OCDB checking
+ } // hdata with entries
+ else AliError("Data type cannot be processed") ;
+ } // loop over histograms from QADataMaker
+ } // end checker
if ( (specie == AliRecoParam::kHighMult) || (specie == AliRecoParam::kLowMult) || (specie == AliRecoParam::kCalib) ) test[specie] = 0.75;
- }
+ } // end of species
}
//____________________________________________________________________________
Double_t AliACORDEQAChecker::CheckAcordeRefHits(const TH1 * href, const TH1 * hdata) const
\r
Float_t meanHitsSL0 = 0.;\r
Float_t meanHitsAMU = 0.;\r
- if ((indexActiveModuleSL0==0) || (indexActiveModuleAMU == 0)) continue;\r
+ //if ((indexActiveModuleSL0==0) || (indexActiveModuleAMU == 0)) continue;\r
+ if ((indexActiveModuleAMU == 0)) continue;\r
\r
meanHitsSL0 = h0->Integral()/indexActiveModuleSL0;\r
meanHitsAMU = h2->Integral()/indexActiveModuleAMU;\r
\r
TH1* h6 = (TH1*)harr[6];\r
- TH1* h7 = (TH1*)harr[7];\r
+ //TH1* h7 = (TH1*)harr[7];\r
\r
- Int_t goodModulesSL0 = 0;\r
- Int_t badModulesSL0 = 0;\r
+ //Int_t goodModulesSL0 = 0;\r
+ //Int_t badModulesSL0 = 0;\r
Int_t goodModulesAMU = 0;\r
Int_t badModulesAMU = 0;\r
\r
for (Int_t imod = 0; imod < 60; imod++)\r
{\r
- if (TMath::Abs(h0->GetBinContent(imod)/meanHitsSL0-1) < 0.65) goodModulesSL0++;\r
- else badModulesSL0++;\r
- if (TMath::Abs(h2->GetBinContent(imod)/meanHitsAMU-1) < 0.65) goodModulesAMU++;\r
- else badModulesAMU++;\r
+ //if (TMath::Abs(h0->GetBinContent(imod)/meanHitsSL0-1) < 0.65) goodModulesSL0++;\r
+ //else badModulesSL0++;\r
+ if (meanHitsAMU!=0)\r
+ {\r
+ if (TMath::Abs(h2->GetBinContent(imod)/meanHitsAMU-1) <= 0.90) goodModulesAMU++;\r
+ else badModulesAMU++;\r
+ }\r
\r
}\r
- h6->SetMaximum(60);\r
- h6->SetMinimum(0);\r
- h6->Fill(1,goodModulesSL0);\r
- h6->Fill(2,badModulesSL0);\r
+ h6->Fill(1,goodModulesAMU);\r
+ h6->Fill(2,badModulesAMU);\r
\r
- h7->SetMaximum(60);\r
- h7->SetMinimum(0);\r
- h7->Fill(1,goodModulesAMU);\r
- h7->Fill(2,badModulesAMU);\r
+ Float_t maxMod = 1*h6->GetMaximum();\r
+ Float_t scaleMaxMod = 0;\r
+ if (maxMod!=0) scaleMaxMod = 1./maxMod;\r
+ else scaleMaxMod = 0;\r
+ h6->Scale(scaleMaxMod);\r
\r
TH1* h4 = (TH1*)harr[4];\r
TH1* h5 = (TH1*)harr[5];\r
\r
\r
\r
- TH1F * fhACORDEStatusSL0_DQM = new TH1F("fhACORDEStatusSL0_DQM","Status of rate for ACORDE's modules (SL0 mode)",2,1,3);\r
TH1F * fhACORDEStatusAMU_DQM = new TH1F("fhACORDEStatusAMU_DQM","Status of rate for ACORDE's modules (AMU mode)",2,1,3);\r
for (Int_t i=0;i<2;i++) \r
{\r
- fhACORDEStatusSL0_DQM->GetXaxis()->SetBinLabel(i+1,acoStatus[i]); \r
fhACORDEStatusAMU_DQM->GetXaxis()->SetBinLabel(i+1,acoStatus[i]); \r
}\r
- fhACORDEStatusSL0_DQM->SetYTitle("No. of modules");\r
- fhACORDEStatusAMU_DQM->SetYTitle("No. of modules");\r
+ fhACORDEStatusAMU_DQM->SetYTitle("No. of modules (Norm.)");\r
\r
- Add2RawsList(fhACORDEStatusSL0_DQM,6,expert,image,!saveCorr);\r
- Add2RawsList(fhACORDEStatusAMU_DQM,7,!expert,image,!saveCorr);\r
+ Add2RawsList(fhACORDEStatusAMU_DQM,6,!expert,image,!saveCorr);\r
\r
- fhACORDEStatusSL0_DQM->SetFillColor(kAzure-5);\r
fhACORDEStatusAMU_DQM->SetFillColor(kAzure-7);\r
\r
\r