/************************************************************************** * 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. * **************************************************************************/ /* Checks the quality assurance. By comparing with reference data Based on PHOS code written by Y. Schutz CERN July 2007 For the moment we only implement the checking of raw data QA. The checked for ESD and RecPoints will be implemented later. */ // --- ROOT system --- #include #include #include #include #include #include #include #include #include #include #include #include // --- Standard library --- // --- AliRoot header files --- #include "AliLog.h" #include "AliQAv1.h" #include "AliQAChecker.h" #include "AliEMCALQAChecker.h" ClassImp(AliEMCALQAChecker) //__________________________________________________________________ AliEMCALQAChecker::AliEMCALQAChecker() : AliQACheckerBase("EMCAL","EMCAL Quality Assurance Data Maker"), fTextSM(new TText*[fgknSM]), fLineCol(new TLine(47.5,-0.5,47.5,119.5)), fText(new TPaveText(0.2,0.7,0.8,0.9,"NDC")) { // ctor fLineCol->SetLineColor(1); fLineCol->SetLineWidth(2); fTextSM[0]= new TText(20, 12, "SM A0"); fTextSM[1]= new TText(20, 36, "SM A1"); fTextSM[2]= new TText(20, 60, "SM A2"); fTextSM[3]= new TText(20, 84, "SM A3"); fTextSM[4]= new TText(20, 108,"SM A4"); fTextSM[5]= new TText(20, 132,"SM A5"); fTextSM[6]= new TText(64, 12, "SM C0"); fTextSM[7]= new TText(64, 36, "SM C1"); fTextSM[8]= new TText(64, 60, "SM C2"); fTextSM[9]= new TText(64, 84, "SM C3"); fTextSM[10]= new TText(64, 108,"SM C4"); fTextSM[11]= new TText(64, 132,"SM C5"); for(int i = 0; i < 5; i++) { fLineRow[i] = new TLine(-0.5,23.5+(24*i),95.5,23.5+(24*i)); fLineRow[i]->SetLineColor(1); fLineRow[i]->SetLineWidth(2); } for(int i = 0; i < 3; i++) { fTextL1[i] = new TPaveText(0.2,0.7,0.8,0.9,"NDC"); } } //__________________________________________________________________ AliEMCALQAChecker::~AliEMCALQAChecker() { /// dtor delete [] fTextSM ; delete fLineCol ; for (Int_t i=0; i<5; ++i) delete fLineRow[i] ; delete fText ; } //______________________________________________________________________________ void AliEMCALQAChecker::Check(Double_t * test, AliQAv1::ALITASK_t index, TObjArray ** list, const AliDetectorRecoParam * /*recoParam*/) { /// Check objects in list if ( index == AliQAv1::kRAW ) { CheckRaws(test, list); printf ("checkers for task %d \n", index) ; } if ( index == AliQAv1::kREC) { CheckRecPoints(test, list); } if ( index == AliQAv1::kESD ) { CheckESD(test, list); } AliWarning(Form("Checker for task %d not implement for the moment",index)); } //______________________________________________________________________________ TH1* AliEMCALQAChecker::GetHisto(TObjArray* list, const char* hname, Int_t specie) const { /// Get a given histo from the list TH1* h = static_cast(list->FindObject(Form("%s_%s",AliRecoParam::GetEventSpecieName(specie),hname))); if (!h) { AliError(Form("Did not find expected histo %s",hname)); } return h; } //______________________________________________________________________________ Double_t AliEMCALQAChecker::MarkHisto(TH1& histo, Double_t value) const { /// Mark histo as originator of some QA error/warning if ( value != 1.0 ) { histo.SetBit(AliQAv1::GetQABit()); } return value; } //______________________________________________________________________________ void AliEMCALQAChecker::CheckRaws(Double_t * test, TObjArray ** list) { // Check RAW QA histograms // -- Yaxian Mao, CCNU/CERN/LPSC //adding new checking method: 25/04/2010, Yaxian Mao //Comparing the amplitude from current run to the reference run, if the ratio in the range [0.8, .12], count as a good tower. //If more than 90% towers are good, EMCAL works fine, otherwise experts should be contacted. Int_t nTowersPerSM = 24*48; // number of towers in a SuperModule; 24x48 Double_t nTot = fgknSM * nTowersPerSM ; TList *lstF = 0; Int_t calibSpecieId = (Int_t)TMath::Log2( AliRecoParam::kCalib ); for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) { test[specie] = 0.0 ; if ( !AliQAv1::Instance()->IsEventSpecieSet(specie)) continue ; if (list[specie]->GetEntries() == 0) test[specie] = 0. ; // nothing to check else { //get calib histos TH2F * hdata = (TH2F*)list[specie]->At(k2DRatioAmp) ; TH1F * ratio = (TH1F*)list[specie]->At(kRatioDist) ; //get L1 histos TH2F *hL1GammaPatch = (TH2F*)list[specie]->At(kGL1); TH2F *hL1JetPatch = (TH2F*)list[specie]->At(kJL1); TH1I *hFrameR = (TH1I*)list[specie]->At(kSTUTRU); //calib histo checker first: if(hdata->GetEntries()!=0 && ratio->GetEntries()!=0) { //adding the lines to distinguish different SMs lstF = hdata->GetListOfFunctions(); CleanListOfFunctions(lstF); lstF->Add(fLineCol->Clone()); for(Int_t iLine = 0; iLine < 4; iLine++) { lstF->Add(fLineRow[iLine]->Clone()); } //Now adding the text to for each SM for(Int_t iSM = 0 ; iSM < fgknSM ; iSM++){ //number of SMs loop start lstF->Add(fTextSM[iSM]->Clone()); } // lstF = ratio->GetListOfFunctions(); CleanListOfFunctions(lstF); // //now check the ratio histogram Double_t binContent = 0. ; Int_t nGoodTower = 0 ; Double_t rv = 0. ; for(Int_t ix = 1; ix <= hdata->GetNbinsX(); ix++) { for(Int_t iy = 1; iy <= hdata->GetNbinsY(); iy++) { binContent = hdata->GetBinContent(ix, iy) ; if (binContent < 1.2 && binContent > 0.8) nGoodTower++ ; } } rv = nGoodTower/nTot ; printf("%2.2f %% towers out of range [0.8, 1.2]\n", (1-rv)*100); if(fText){ lstF->Add(fText->Clone()) ; fText->Clear() ; fText->AddText(Form("%2.2f %% towers out of range [0.8, 1.2]", (1-rv)*100)); if (rv < 0.9) { test[specie] = 0.9 ; // 2 lines text info for quality fText->SetFillColor(2) ; fText->AddText(Form("EMCAL = NOK, CALL EXPERTS!!!")); } else { test[specie] = 0.1 ; fText->SetFillColor(3) ; fText->AddText(Form("EMCAL = OK, ENJOY...")); } }//fText } // calib histo checking done //now L1 checks: if (specie != calibSpecieId) { if(hL1GammaPatch->GetEntries() !=0 ) { lstF = hL1GammaPatch->GetListOfFunctions(); CleanListOfFunctions(lstF); //Checker for L1GammaPatch (if a patch triggers > sigmaG * mean value (1/#patch positions total) says "hot spot !") Double_t dL1GmeanTrig = 1/2961.; Int_t sigmaG = 5; // deviation from mean value Double_t dL1GEntries = hL1GammaPatch->GetEntries(); Int_t badL1G[48][64] = {{0}} ; Int_t nBadL1G = 0; for(Int_t ix = 1; ix <= hL1GammaPatch->GetNbinsX(); ix++) { for(Int_t iy = 1; iy <= hL1GammaPatch->GetNbinsY(); iy++) { Double_t binContent = hL1GammaPatch->GetBinContent(ix, iy) ; if (binContent != 0) { if ((double)binContent/(double)dL1GEntries > sigmaG*dL1GmeanTrig) { badL1G[ix-1][iy-1] += 1; nBadL1G += 1; } } } } if(fTextL1[0]){ lstF->Add(fTextL1[0]->Clone()) ; fTextL1[0]->Clear() ; if (nBadL1G == 0) { fTextL1[0]->SetFillColor(3) ; fTextL1[0]->AddText(Form("L1 GAMMA TRIGGER = OK, ENJOY...")); } else { fTextL1[0]->SetFillColor(2) ; fTextL1[0]->AddText(Form("HOT SPOT IN L1 GAMMA TRIGGER = CALL EXPERT!!")); /* for(Int_t ix = 1; ix <= hL1GammaPatch->GetNbinsX(); ix++) { for(Int_t iy = 1; iy <= hL1GammaPatch->GetNbinsY(); iy++) { if(badL1G[ix-1][iy-1] != 0) printf("L1 Gamma patch with position x = %d, y = %d is out of range\n",ix,iy); } } */ } }//fTextL1[0] }// L1 gamma patch checking done if(hL1JetPatch->GetEntries() !=0) { lstF = hL1JetPatch->GetListOfFunctions(); CleanListOfFunctions(lstF); //Checker for L1JetPatch (if a patch triggers > sigmaJ * mean value (1/#patch positions total) says "hot spot !") Double_t dL1JmeanTrig = 1/126.; Int_t sigmaJ = 5; // deviation from mean value Double_t dL1JEntries = hL1JetPatch->GetEntries(); Int_t badL1J[12][16] = {{0}} ; Int_t nBadL1J = 0; for(Int_t ix = 1; ix <= hL1JetPatch->GetNbinsX(); ix++) { for(Int_t iy = 1; iy <= hL1JetPatch->GetNbinsY(); iy++) { Double_t binContent = hL1JetPatch->GetBinContent(ix, iy) ; if (binContent != 0) { if ((double)binContent/(double)dL1JEntries > sigmaJ*dL1JmeanTrig) { badL1J[ix-1][iy-1] += 1 ; nBadL1J += 1; } } } } if(fTextL1[1]){ lstF->Add(fTextL1[1]->Clone()) ; fTextL1[1]->Clear() ; if (nBadL1J == 0) { fTextL1[1]->SetFillColor(3) ; fTextL1[1]->AddText(Form("L1 JET TRIGGER = OK, ENJOY...")); } else { fTextL1[1]->SetFillColor(2) ; fTextL1[1]->AddText(Form("HOT SPOT IN L1 JET TRIGGER = CALL EXPERT!!")); /* for(Int_t ix = 1; ix <= hL1JetPatch->GetNbinsX(); ix++) { for(Int_t iy = 1; iy <= hL1JetPatch->GetNbinsY(); iy++) { if(badL1J[ix-1][iy-1] != 0) printf("L1 Jet patch with position x = %d, y = %d is out of range\n",(4*ix-4),(4*iy-4)); } } */ } }//fTextL1[1] } // L1 Jet patch checking done } // if (specie != calibSpecieId) .. if(hFrameR->GetEntries() !=0) { lstF = hFrameR->GetListOfFunctions(); CleanListOfFunctions(lstF); Int_t badLink[32] = {0}; Int_t nBadLink = 0; for(Int_t ix = 1; ix <= hFrameR->GetNbinsX(); ix++) { Double_t binContent = hFrameR->GetBinContent(ix) ; if (binContent == 0) { badLink[ix-1] += 1; nBadLink += 1; } } if(fTextL1[2]){ lstF->Add(fTextL1[2]->Clone()) ; fTextL1[2]->Clear() ; if (nBadLink == 0) { fTextL1[2]->SetFillColor(3) ; fTextL1[2]->AddText(Form("LINK TRU-STU = OK, ENJOY...")); } else { fTextL1[2]->SetFillColor(2) ; fTextL1[2]->AddText(Form("PROBLEM WITH TRU-STU LINK = CALL EXPERT!!")); /* for(Int_t ix = 0; ix <= hFrameR->GetNbinsX(); ix++) { if(badLink[ix] != 0) printf("STU link with TRU %d is out\n",ix); } */ } }//fTextL1[2] } // Checker for link TRU-STU done } } // specie } //______________________________________________________________________________ void AliEMCALQAChecker::Init(const AliQAv1::DETECTORINDEX_t det) { /// intialises QA and QA checker settings AliQAv1::Instance(det) ; Float_t hiValue[AliQAv1::kNBIT] ; Float_t lowValue[AliQAv1::kNBIT] ; lowValue[AliQAv1::kINFO] = 0.0 ; hiValue[AliQAv1::kINFO] = 0.1 ; lowValue[AliQAv1::kWARNING] = 0.1 ; hiValue[AliQAv1::kWARNING] = 0.5 ; lowValue[AliQAv1::kERROR] = 0.5 ; hiValue[AliQAv1::kERROR] = 0.8 ; lowValue[AliQAv1::kFATAL] = 0.8 ; hiValue[AliQAv1::kFATAL] = 1.0 ; SetHiLo(&hiValue[0], &lowValue[0]) ; } //______________________________________________________________________________ void AliEMCALQAChecker::CleanListOfFunctions(TList *list) { // clean up if (list) { TObject *stats = list->FindObject("stats"); list->Remove(stats); TObject *obj; while ((obj = list->First())) { while(list->Remove(obj)) { } delete obj; } if (stats) list->Add(stats); } else { AliWarning(Form("Checker : empty list of data functions; returning")); return; } }