DQM flag reassigned due to new PMT's calibration
[u/mrichter/AliRoot.git] / ACORDE / AliACORDEQAChecker.cxx
CommitLineData
51504028 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15//...
16// Checks the quality assurance for ACORDE.
17// Default implementation
18// Authors:
19// Mario Rodriguez Cahuantzi <mrodrigu@mail.cern.ch> (FCFM-BUAP)
20// Luciano Diaz Gonzalez <luciano.diaz@nucleares.unam.mx> (ICN-UNAM)
21// Arturo Fernandez <afernan@mail.cern.ch> (FCFM-BUAP)
1b368017 22// Last update: Nov. 14t 2009 --> MRC <mrodrigu@mail.cern.ch> (FCFM-BUAP)
51504028 23//...
24
25// --- ROOT system ---
26#include <TClass.h>
27#include <TH1F.h>
28#include <TH1I.h>
29#include <TIterator.h>
30#include <TKey.h>
31#include <TFile.h>
dd1da44f 32#include <TPaveText.h>
33#include <TLine.h>
51504028 34// --- Standard library ---
35
36// --- AliRoot header files ---
37#include "AliLog.h"
4e25ac79 38#include "AliQAv1.h"
51504028 39#include "AliQAChecker.h"
40#include "AliACORDEQAChecker.h"
fcdfb1ff 41#include "AliCDBEntry.h"
42#include "AliQAManager.h"
51504028 43
dd1da44f 44/*************************************************************************
45 Last update: Oct. 10th 2011 from Mario RC, <mrodrigu@mail.cern.ch>
46 |-> Adding the checker class for raw and esd index
47 |-> Setting the threshold lines and box for DQM shifter
48
49*************************************************************************/
50
51504028 51ClassImp(AliACORDEQAChecker)
52
1b368017 53//____________________________________________________________________________
dd1da44f 54
55AliACORDEQAChecker::AliACORDEQAChecker() :
56AliQACheckerBase("ACORDE","ACORDE Quality Assurance Data Maker"),
57fTextDQMShifterInfo(new TPaveText(35,0.5,55,0.7,"T")),
58fMax(new TLine(0,0.9,60,0.9)),
59fMin(new TLine(0,0.3,60,0.3))
60{
61 // default constructor
62 fMax->SetLineColor(kGreen);
63 fMax->SetLineWidth(3);
64 fMin->SetLineColor(kGreen);
65 fMax->SetLineWidth(3);
66}
67//____________________________________________________________________________
68AliACORDEQAChecker::~AliACORDEQAChecker()
69{
70 // destructor
71 delete fTextDQMShifterInfo;
72 delete fMax;
73 delete fMin;
74}
75//____________________________________________________________________________
76AliACORDEQAChecker::AliACORDEQAChecker(const AliACORDEQAChecker& qac) :
77AliQACheckerBase(qac.GetName(), qac.GetTitle()),
78fTextDQMShifterInfo(new TPaveText(35,0.5,55,0.7,"T")),
79fMax(static_cast<TLine*>(qac.fMax->Clone())),
80fMin(static_cast<TLine*>(qac.fMin->Clone()))
81{
82 //
83}
84//____________________________________________________________________________
85AliACORDEQAChecker& AliACORDEQAChecker::operator = (const AliACORDEQAChecker &qac)
86{
87
88 fTextDQMShifterInfo = new TPaveText(35,0.5,55,0.7,"T");
89 fMax = static_cast<TLine*>(qac.fMax->Clone());
90 fMin = static_cast<TLine*>(qac.fMin->Clone());
91 return *this;
92}
93//____________________________________________________________________________
86017bd8 94void AliACORDEQAChecker::Check(Double_t * test, AliQAv1::ALITASK_t /*index*/)
1b368017 95{
1b368017 96 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
86017bd8 97 test[specie] = 0.0 ;
1b368017 98}
99//____________________________________________________________________________
86017bd8 100void AliACORDEQAChecker::Check(Double_t * test, AliQAv1::ALITASK_t /*index*/, TObjArray ** list, const AliDetectorRecoParam * /*recoParam*/)
f307c383 101{
1b368017 102// Close version to the final one for the ACORDE QA Checker
57acd2d2 103
1b368017 104// Loop over the run species (for specie!= cosmic by now we set QA to INFO)
105
106 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
107 {
108 if ( !AliQAv1::Instance()->IsEventSpecieSet(specie) ) continue ;
86017bd8 109 if (list[specie]->GetEntries() == 0) test[specie] = 1.; // Nothing to check
1b368017 110 else
bd97a569 111 {
1b368017 112 TIter next(list[specie]) ;
113 TH1 * hdata ; // Data created by the AliACORDEQADataMakerXXX (Sim/Rec)
114 while ( (hdata = dynamic_cast<TH1 *>(next())) )
fcdfb1ff 115 {
1b368017 116 if (hdata)
bd97a569 117 {
1b368017 118 Double_t rv = 0.0 ;
119 if(hdata->GetEntries()>0) rv=1;
120 AliDebug(AliQAv1::GetQADebugLevel(), Form("%s -> %f", hdata->GetName(), rv)) ;
121 TString hdataName = hdata->GetName();
dd1da44f 122 if (hdata->GetListOfFunctions()->GetEntries() == 0 && (hdataName.Contains("HitsSL0_DQM")||(hdataName.Contains("HitsAMU_DQM")))){
123 hdata->GetListOfFunctions()->Add(fTextDQMShifterInfo);
124 hdata->GetListOfFunctions()->Add(fMax);
125 hdata->GetListOfFunctions()->Add(fMin);
126
127 }
1b368017 128 // Here we use the QAref ACORDE data from fRefOCDBSubDir
129
130 if ( (fRefOCDBSubDir[specie]) && (hdataName.Contains("ACORDEBitPattern")) )
bd97a569 131 {
1b368017 132 TH1 * href = NULL;
133 if (fRefSubDir) href = static_cast<TH1*>(fRefSubDir->Get(hdata->GetName()));
134 else if (fRefOCDBSubDir[specie]) href = static_cast<TH1*>(fRefOCDBSubDir[specie]->FindObject(hdata->GetName()));
86017bd8 135 test[specie] = CheckAcordeRefHits(href,hdata);
dd1da44f 136 }else if (hdataName.Contains("ACORDEBitPattern")||hdataName.Contains("HitsSL0_DQM")||hdataName.Contains("HitsAMU_DQM"))
1b368017 137 // Here we use an inner QA Checher without the QAref data
bd97a569 138 {
1b368017 139 Float_t acoDataMax = hdata->GetMaximum();
dd1da44f 140 Int_t flagAcoQAChecker = 0;
1b368017 141 Int_t flagAcoQA = 0;
dd1da44f 142 if (acoDataMax == 0) continue;
143 for(Int_t i=0;i<60;i++)
bd97a569 144 {
dd1da44f 145 if ((hdata->GetBinContent(i)/acoDataMax) < 0.75) flagAcoQAChecker++;
146 if (hdataName.Contains("HitsSL0_DQM")||hdataName.Contains("HitsAMU_DQM")){
147 if (hdata->GetBinContent(i)<0.29) flagAcoQA++;
f31df1cb 148 }
bd97a569 149 }
dd1da44f 150 if (hdataName.Contains("ACORDEBitPattern")){
151 Double_t simpleFlag = 1.-flagAcoQAChecker/60.;
152 if ( (simpleFlag >= 0.90) && (simpleFlag <= 1.0) ) test[specie] = 0.75; // INFO
153 if ( (simpleFlag >= 0.70) && (simpleFlag < 0.90) ) test[specie] = 0.50; // WARNING
154 if ( (simpleFlag >= 0.25) && (simpleFlag < 0.70) ) test[specie] = 0.25; // ERROR
155 if ( (simpleFlag >= 0.0) && (simpleFlag < 0.25) ) test[specie] = -1.0; // FATAL
156 } else if (hdataName.Contains("HitsSL0_DQM")||hdataName.Contains("HitsAMU_DQM")){
2fe32e4c 157 if (flagAcoQA < 100) test[specie] = 0.75;
dd1da44f 158 else test[specie] = 0.3;
159 }
160
161
1b368017 162 }
dd1da44f 163
1b368017 164 // Setting Warning message for possible Empty Events with the ACORDE-Trigger
165
166 if (hdataName.Contains("ACORDEMultiplicity") && (hdata->GetBinContent(0)!=0)) AliWarning("Empty event registered with ACORDE Trigger !!!");
167
168
169 }else AliError("Data type cannot be processed") ;
dd1da44f 170 if (fTextDQMShifterInfo)
171 {
172 fTextDQMShifterInfo->Clear();
173 if (test[specie] > 0.3){
174 fTextDQMShifterInfo->SetFillColor(kGreen);
175 fTextDQMShifterInfo->AddText("ACORDE: O.K.");
176 }else{
177 fTextDQMShifterInfo->SetFillColor(kRed);
178 fTextDQMShifterInfo->AddText("ACORDE: Not, O.K.");
179 fTextDQMShifterInfo->AddText("CALL THE EXPERTS");
180 }
181 }
1b368017 182
183 }
184 }
86017bd8 185 if ( (specie == AliRecoParam::kHighMult) || (specie == AliRecoParam::kLowMult) || (specie == AliRecoParam::kCalib) ) test[specie] = 0.75;
1b368017 186 }
f307c383 187}
1b368017 188//____________________________________________________________________________
189Double_t AliACORDEQAChecker::CheckAcordeRefHits(const TH1 * href, const TH1 * hdata) const
fcdfb1ff 190{
86017bd8 191 Double_t test = 0.;
1b368017 192 Int_t flag=0;
193 for (Int_t i=0;i<60;i++)
fcdfb1ff 194 {
0200d7cf 195 if (TMath::Abs(href->GetBinContent(i)-hdata->GetBinContent(i))>10) flag++;
bd97a569 196 }
86017bd8 197 if ((flag>50)&&(flag<=60)) test = -1.;
198 if ((flag>30)&&(flag<=50)) test = 0.25;
199 if ((flag>10)&&(flag<=30)) test = 0.5;
200 if ((flag>0)&&(flag<=10)) test = 0.75;
201 return test;
fcdfb1ff 202}