]> git.uio.no Git - u/mrichter/AliRoot.git/blame - STEER/STEER/AliQACheckerBase.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
added number of crossed rows (Markus Koehler)
[u/mrichter/AliRoot.git] / STEER / STEER / AliQACheckerBase.cxx
CommitLineData
421ab0fb 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
17/* $Id$ */
18
202374b1 19//
20// Base class for detectors quality assurance checkers
21// Compares Data made by QADataMakers with reference data
22// Y. Schutz CERN August 2007
23//
421ab0fb 24
25// --- ROOT system ---
634696f5 26#include <TCanvas.h>
421ab0fb 27#include <TClass.h>
28#include <TH1F.h>
29#include <TH1I.h>
30#include <TIterator.h>
31#include <TKey.h>
32#include <TFile.h>
a4976ef3 33#include <TList.h>
a2b64fbd 34#include <TNtupleD.h>
11a18a03 35#include <TParameter.h>
634696f5 36#include <TPaveText.h>
421ab0fb 37
38// --- Standard library ---
39
40// --- AliRoot header files ---
78decf65 41#include "AliCDBEntry.h"
421ab0fb 42#include "AliLog.h"
4e25ac79 43#include "AliQAv1.h"
2e42b4d4 44#include "AliQAChecker.h"
45#include "AliQACheckerBase.h"
46#include "AliQADataMaker.h"
78decf65 47#include "AliQAManager.h"
c26cb20d 48#include "AliDetectorRecoParam.h"
421ab0fb 49
2e42b4d4 50ClassImp(AliQACheckerBase)
421ab0fb 51
52
53//____________________________________________________________________________
2e42b4d4 54AliQACheckerBase::AliQACheckerBase(const char * name, const char * title) :
421ab0fb 55 TNamed(name, title),
a5fa6165 56 fDataSubDir(0x0),
4edbc5bc 57 fRefSubDir(0x0),
aa2beccb 58 fRefOCDBSubDir(new TObjArray*[AliRecoParam::kNSpecies]),
a04582a1 59 fLowTestValue(new Float_t[AliQAv1::kNBIT]),
60 fUpTestValue(new Float_t[AliQAv1::kNBIT]),
634696f5 61 fImage(new TCanvas*[AliRecoParam::kNSpecies]),
11a18a03 62 fPrintImage(kTRUE),
a80d46d0 63 fExternParamList(new TList())
421ab0fb 64{
65 // ctor
4e25ac79 66 fLowTestValue[AliQAv1::kINFO] = 0.5 ;
67 fUpTestValue[AliQAv1::kINFO] = 1.0 ;
68 fLowTestValue[AliQAv1::kWARNING] = 0.002 ;
69 fUpTestValue[AliQAv1::kWARNING] = 0.5 ;
70 fLowTestValue[AliQAv1::kERROR] = 0.0 ;
71 fUpTestValue[AliQAv1::kERROR] = 0.002 ;
72 fLowTestValue[AliQAv1::kFATAL] = -1.0 ;
73 fUpTestValue[AliQAv1::kFATAL] = 0.0 ;
57acd2d2 74
5379c4a3 75 AliDebug(AliQAv1::GetQADebugLevel(), "Default setting is:") ;
76 if ( AliDebugLevel() == AliQAv1::GetQADebugLevel() ) {
eca4fa66 77 const Char_t * text= Form(" INFO -> %1.5f < value < %1.5f WARNING -> %1.5f < value <= %1.5f \n ERROR -> %1.5f < value <= %1.5f \n FATAL -> %1.5f <= value < %1.5f \n",
78 fLowTestValue[AliQAv1::kINFO], fUpTestValue[AliQAv1::kINFO],
79 fLowTestValue[AliQAv1::kWARNING], fUpTestValue[AliQAv1::kWARNING],
80 fLowTestValue[AliQAv1::kERROR], fUpTestValue[AliQAv1::kERROR],
81 fLowTestValue[AliQAv1::kFATAL], fUpTestValue[AliQAv1::kFATAL]) ;
82 AliInfo(Form("%s", text)) ;
5379c4a3 83 }
634696f5 84
aa2beccb 85 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
634696f5 86 fImage[specie] = NULL ;
aa2beccb 87 fRefOCDBSubDir[specie] = NULL ;
88 }
421ab0fb 89}
90
91//____________________________________________________________________________
4edbc5bc 92AliQACheckerBase::AliQACheckerBase(const AliQACheckerBase& qac) :
93 TNamed(qac.GetName(), qac.GetTitle()),
94 fDataSubDir(qac.fDataSubDir),
95 fRefSubDir(qac.fRefSubDir),
57acd2d2 96 fRefOCDBSubDir(qac.fRefOCDBSubDir),
a04582a1 97 fLowTestValue(new Float_t[AliQAv1::kNBIT]),
98 fUpTestValue(new Float_t[AliQAv1::kNBIT]),
99 fImage(new TCanvas*[AliRecoParam::kNSpecies]),
11a18a03 100 fPrintImage(kTRUE),
a80d46d0 101 fExternParamList(new TList())
421ab0fb 102{
103 //copy ctor
4e25ac79 104 for (Int_t index = 0 ; index < AliQAv1::kNBIT ; index++) {
57acd2d2 105 fLowTestValue[index] = qac.fLowTestValue[index] ;
106 fUpTestValue[index] = qac.fUpTestValue[index] ;
107 }
aa2beccb 108 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
109 fImage[specie] = qac.fImage[specie] ;
110 fRefOCDBSubDir[specie] = qac.fRefOCDBSubDir[specie] ;
111 }
11a18a03 112 if (qac.fExternParamList) {
11a18a03 113 TIter next(qac.fExternParamList) ;
114 TParameter<double> * p ;
115 while ( (p = (TParameter<double>*)next()) )
116 fExternParamList->Add(p) ;
117 }
421ab0fb 118}
119
120//____________________________________________________________________________
634696f5 121AliQACheckerBase& AliQACheckerBase::operator = (const AliQACheckerBase& qac )
421ab0fb 122{
123 // Equal operator.
2e42b4d4 124 this->~AliQACheckerBase();
634696f5 125 new(this) AliQACheckerBase(qac);
421ab0fb 126 return *this;
127}
128
57acd2d2 129//____________________________________________________________________________
130AliQACheckerBase::~AliQACheckerBase()
131{
132 delete [] fLowTestValue ;
133 delete [] fUpTestValue ;
a4e305ee 134 DeleteImages();
6b5405c1 135 delete[] fImage ;
aa2beccb 136 delete[] fRefOCDBSubDir ;
5d1a2a7d 137 AliQAv1::GetQAResultFile()->Close() ;
11a18a03 138 if (fExternParamList) {
139 fExternParamList->Clear() ;
140 delete fExternParamList ;
141 }
57acd2d2 142}
143
a5fa6165 144//____________________________________________________________________________
75a8b78a 145void AliQACheckerBase::PrivateCheck(Double_t * test, AliQAv1::ALITASK_t index, const AliDetectorRecoParam * recoParam)
a5fa6165 146{
147 // Performs a basic checking
148 // Compares all the histograms stored in the directory
4edbc5bc 149 // With reference histograms either in a file of in OCDB
a5fa6165 150
ed448a70 151 TObjArray ** list = new TObjArray *[AliRecoParam::kNSpecies] ;
c39ee44c 152 Int_t specie ;
153 for (specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
ed448a70 154 list[specie] = new TObjArray(AliQAv1::GetMaxQAObj()) ;
4e25ac79 155 if ( !AliQAv1::Instance()->IsEventSpecieSet(specie) )
57acd2d2 156 continue ;
c39ee44c 157 if (fDataSubDir) {
ed448a70 158 TList * keyList = fDataSubDir->GetListOfKeys() ;
159 TIter next(keyList) ;
160 TKey * key ;
161 while ( (key = static_cast<TKey *>(next())) ) {
162 TDirectory * specieDir = fDataSubDir->GetDirectory(key->GetName()) ;
163 TList * keykeyList = specieDir->GetListOfKeys() ;
164 TIter next2(keykeyList) ;
165 TKey * keykey ;
166 while ( (keykey = static_cast<TKey *>(next2())) ) {
167 TObject * odata = specieDir->Get(keykey->GetName()) ;
168 if ( odata->IsA()->InheritsFrom("TH1") ) {
169 TH1 * hdata = static_cast<TH1*>(odata) ;
170 list[specie]->Add(hdata) ;
171 } else if (!odata->IsA()->InheritsFrom("TDirectory")) // skip the expert directory
172 AliError(Form("%s Is a Classname that cannot be processed", key->GetClassName())) ;
c39ee44c 173 }
57acd2d2 174 }
57acd2d2 175 }
176 }
ed448a70 177
86017bd8 178 Check(test, index, list, recoParam) ;
ed448a70 179
180 delete[] list ;
181
a5fa6165 182}
183
a4976ef3 184//____________________________________________________________________________
a7e91ced 185void AliQACheckerBase::Check(Double_t * test, AliQAv1::ALITASK_t task, TObjArray ** list, const AliDetectorRecoParam * recoParam)
a4976ef3 186{
187 // Performs a basic checking
188 // Compares all the histograms in the list
189
92664bc8 190 Int_t count[AliRecoParam::kNSpecies] = { 0 };
57acd2d2 191
78decf65 192 GetRefSubDir(GetName(), AliQAv1::GetTaskName(task), fRefSubDir, fRefOCDBSubDir) ;
193 // SetRefandData(refDir, refOCDBDir) ;
194
57acd2d2 195 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
196 test[specie] = 1.0 ;
92664bc8 197 if ( !AliQAv1::Instance()->IsEventSpecieSet(specie)) continue ;
57acd2d2 198 if (list[specie]->GetEntries() == 0)
199 test[specie] = 0. ; // nothing to check
200 else {
201 if (!fRefSubDir && !fRefOCDBSubDir)
202 test[specie] = -1 ; // no reference data
203 else {
204 TIter next(list[specie]) ;
205 TH1 * hdata ;
206 count[specie] = 0 ;
eca4fa66 207 while ( (hdata = static_cast<TH1 *>(next())) ) {
ab1bb842 208 if ( hdata->IsA()->InheritsFrom("TH1") ) {
7775195b 209 if ( hdata->TestBit(AliQAv1::GetExpertBit()) ) // does not perform the test for expert data
210 continue ;
92664bc8 211 //
212 // First try to find the ref histo with exact name (with possible trigger clon ending)
213 TString hname = hdata->GetName();
214 TH1 * href = NULL ;
215 if (fRefSubDir) href = static_cast<TH1*>(fRefSubDir->Get(hname.Data())) ;
216 else if (fRefOCDBSubDir[specie]) href = static_cast<TH1*>(fRefOCDBSubDir[specie]->FindObject(hname.Data()));
217 //
218 if (!href && hdata->TestBit(AliQAv1::GetClonedBit())) { // try to find the histo for the base name (w/o trigger ending
219 int ind = hname.Index(AliQADataMaker::GetTriggerPrefix());
220 if (ind>0) {
221 hname.Resize(ind);
222 if (fRefSubDir) href = static_cast<TH1*>(fRefSubDir->Get(hname.Data())) ;
223 else if (fRefOCDBSubDir[specie]) href = static_cast<TH1*>(fRefOCDBSubDir[specie]->FindObject(hname.Data()));
224 }
225 }
226 //
7775195b 227 if (!href)
228 test[specie] = -1 ; // no reference data ;
229 else {
230 Double_t rv = DiffK(hdata, href) ;
231 AliDebug(AliQAv1::GetQADebugLevel(), Form("%s ->Test = %f", hdata->GetName(), rv)) ;
232 test[specie] += rv ;
233 count[specie]++ ;
57acd2d2 234 }
c39ee44c 235 } else
57acd2d2 236 AliError("Data type cannot be processed") ;
c39ee44c 237 if (count[specie] != 0)
238 test[specie] /= count[specie] ;
57acd2d2 239 }
57acd2d2 240 }
241 }
242 }
a4976ef3 243}
244
57acd2d2 245
a4e305ee 246//____________________________________________________________________________
247void AliQACheckerBase::DeleteImages()
248{
6b5405c1 249 // clean images
250 for (Int_t esIndex = 0 ; esIndex < AliRecoParam::kNSpecies ; esIndex++) {
251 if ( fImage[esIndex] ) {delete fImage[esIndex]; fImage[esIndex] = 0;}
252 if ( fRefOCDBSubDir[esIndex] ) {delete fRefOCDBSubDir[esIndex]; fRefOCDBSubDir[esIndex] = 0;}
253 }
a4e305ee 254}
255
421ab0fb 256//____________________________________________________________________________
f12d42ce 257Double_t AliQACheckerBase::DiffC(const TH1 * href, const TH1 * hin) const
421ab0fb 258{
259 // compares two histograms using the Chi2 test
260 if ( hin->Integral() == 0 ) {
5379c4a3 261 AliDebug(AliQAv1::GetQADebugLevel(), Form("Spectrum %s is empty", hin->GetName())) ;
421ab0fb 262 return 0. ;
263 }
264
265 return hin->Chi2Test(href) ;
266}
267
268//____________________________________________________________________________
f12d42ce 269Double_t AliQACheckerBase::DiffK(const TH1 * href, const TH1 * hin) const
421ab0fb 270{
271 // compares two histograms using the Kolmogorov test
6252ceeb 272 if ( hin->Integral() == 0 || href->Integral() == 0) {
273 AliDebug(AliQAv1::GetQADebugLevel(), Form("Spectrum %s or its reference is empty", hin->GetName())) ;
421ab0fb 274 return 0. ;
275 }
276
277 return hin->KolmogorovTest(href) ;
278}
279
78decf65 280 //_____________________________________________________________________________
281void AliQACheckerBase::GetRefSubDir(const char * det, const char * task, TDirectory *& dirFile, TObjArray **& dirOCDB)
282{
283 // Opens and returns the file with the reference data
284 dirFile = NULL ;
285 TString refStorage(AliQAv1::GetQARefStorage()) ;
286 if (!refStorage.Contains(AliQAv1::GetLabLocalOCDB()) && !refStorage.Contains(AliQAv1::GetLabAliEnOCDB())) {
287 AliError(Form("%s is not a valid location for reference data", refStorage.Data())) ;
288 return ;
289 } else {
290 AliQAManager* manQA = AliQAManager::QAManager(AliQAv1::GetTaskIndex(task)) ;
a04582a1 291 // dirOCDB = new TObjArray*[AliRecoParam::kNSpecies] ;
78decf65 292 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
293 dirOCDB[specie] = NULL ;
294 if ( !AliQAv1::Instance()->IsEventSpecieSet(specie) )
295 continue ;
296 AliQAv1::SetQARefDataDirName(specie) ;
297 if ( ! manQA->GetLock() ) {
298 manQA->SetDefaultStorage(AliQAv1::GetQARefStorage()) ;
299 manQA->SetSpecificStorage("*", AliQAv1::GetQARefStorage()) ;
300 manQA->SetRun(AliCDBManager::Instance()->GetRun()) ;
301 manQA->SetLock() ;
302 }
303 char * detOCDBDir = Form("%s/%s/%s", det, AliQAv1::GetRefOCDBDirName(), AliQAv1::GetRefDataDirName()) ;
304 AliCDBEntry * entry = manQA->Get(detOCDBDir, manQA->GetRun()) ;
305 if (entry) {
306 TList * listDetQAD =static_cast<TList *>(entry->GetObject()) ;
1c871ee3 307 if ( listDetQAD && strcmp(listDetQAD->ClassName(), "TList") != 0 ) {
78decf65 308 AliError(Form("Expected a Tlist and found a %s for detector %s", listDetQAD->ClassName(), det)) ;
309 listDetQAD = NULL ;
310 continue ;
311 }
312 if ( listDetQAD ) {
313 TIter next(listDetQAD) ;
d90c4827 314 while ( (TObjArray*)next() )
78decf65 315 dirOCDB[specie] = static_cast<TObjArray *>(listDetQAD->FindObject(Form("%s/%s", task, AliRecoParam::GetEventSpecieName(specie)))) ;
316 }
317 }
318 }
319 }
320}
321
c39ee44c 322//____________________________________________________________________________
11a18a03 323void AliQACheckerBase::PrintExternParam()
324{
325 // Print the list of external parameter list
326 TIter next(fExternParamList) ;
327 TParameter<double> *pp ;
328 TString printit("\n") ;
329 while( (pp = (TParameter<double>*)next()) )
330 printit += Form("%s = %f\n", pp->GetName(), pp->GetVal());
331 AliInfo(Form("%s", printit.Data())) ;
332}
333
334//____________________________________________________________________________
9025307b 335void AliQACheckerBase::Run(AliQAv1::ALITASK_t index, const AliDetectorRecoParam * recoParam)
c39ee44c 336{
9025307b 337 //Run the checker for all kind of species
92664bc8 338 AliDebug(AliQAv1::GetQADebugLevel(), Form("Processing %s", AliQAv1::GetAliTaskName(index))) ;
c39ee44c 339
92664bc8 340 Double_t * rv = new Double_t[AliRecoParam::kNSpecies] ;
255c1c92 341 for (int i=AliRecoParam::kNSpecies;i--;) rv[i] = 0.0;
75a8b78a 342 PrivateCheck(rv, index, recoParam) ;
92664bc8 343 SetQA(index, rv) ;
344
c39ee44c 345 AliDebug(AliQAv1::GetQADebugLevel(), Form("Test result of %s", AliQAv1::GetAliTaskName(index))) ;
92664bc8 346
d9e8b603 347 delete [] rv ;
c39ee44c 348 Finish() ;
349}
350
421ab0fb 351//____________________________________________________________________________
9025307b 352void AliQACheckerBase::Run(AliQAv1::ALITASK_t index, TObjArray ** list, const AliDetectorRecoParam * recoParam)
421ab0fb 353{
92664bc8 354 // RS: perform check for all trigger classes in loop
355 Double_t * rv = new Double_t[AliRecoParam::kNSpecies] ;
356 //
357 TObjArray ** listTrig = new TObjArray *[AliRecoParam::kNSpecies];
358 //
359 for (int itc=-1;itc<AliQADataMaker::GetNTrigClasses();itc++) {
360 //
361 // RS: fetch the histograms for each specie and for given trigger
362 //AliInfo(Form("Processing %s for trigger: %s", AliQAv1::GetAliTaskName(index),AliQADataMaker::GetTrigClassName(itc)));
363
364 for (int specie=0;specie<AliRecoParam::kNSpecies;specie++) {
365 listTrig[specie] = 0;
366 if ( !AliQAv1::Instance()->IsEventSpecieSet(specie) || !list[specie]) continue;
367 listTrig[specie] = new TObjArray( list[specie]->GetSize() ); // destination for clones of this trigger
368 AliQADataMaker::GetDataOfTrigClass(list[specie],itc, listTrig[specie]);
369 }
370 AliDebug(AliQAv1::GetQADebugLevel(), Form("Processing %s for trigger: %s", AliQAv1::GetAliTaskName(index),AliQADataMaker::GetTrigClassName(itc)));
371 Check(rv, index, listTrig, recoParam) ;
372 SetQA(index, rv) ;
373 AliDebug(AliQAv1::GetQADebugLevel(), Form("Test result of %s for trigger: %s", AliQAv1::GetAliTaskName(index),AliQADataMaker::GetTrigClassName(itc)));
374 //
375 for (int specie=0;specie<AliRecoParam::kNSpecies;specie++) if (listTrig[specie]) delete listTrig[specie]; // clean temporary container
376 }
86017bd8 377 delete [] rv ;
385b26a8 378 delete [] listTrig;
a2b64fbd 379 Finish() ;
421ab0fb 380}
381
382//____________________________________________________________________________
2e42b4d4 383void AliQACheckerBase::Finish() const
421ab0fb 384{
92664bc8 385 // wrap up and save QA in proper file
386 AliQAv1::GetQAResultFile() ;
387 AliQAv1 * qa = AliQAv1::Instance() ;
388 qa->Write(AliQAv1::GetQAName(), kWriteDelete) ;
421ab0fb 389}
46b8a84d 390
634696f5 391//____________________________________________________________________________
392void AliQACheckerBase::MakeImage( TObjArray ** list, AliQAv1::TASKINDEX_t task, AliQAv1::MODE_t mode)
393{
394 // makes the QA image for sim and rec
92664bc8 395 TObjArray tmpArr; // array to store flat version of original array (which may contain clones)
396 //
397 for (Int_t esIndex = 0; esIndex < AliRecoParam::kNSpecies; esIndex++) {
398 if (! AliQAv1::Instance(AliQAv1::GetDetIndex(GetName()))->IsEventSpecieSet(AliRecoParam::ConvertIndex(esIndex)) || list[esIndex]->GetEntries() == 0) continue;
399 Int_t nImages = 0;
400 TIter next(list[esIndex]);
401 TObject* hdata = NULL;
402 tmpArr.Clear();
403 while ( (hdata=(next())) ) { // count histos and transfere to flat array
404 if (hdata->InheritsFrom(TH1::Class()) && hdata->TestBit(AliQAv1::GetImageBit()) ) { // histo, not cloned
405 nImages++;
406 tmpArr.AddLast(hdata);
407 continue;
94b94b5f 408 }
92664bc8 409 if (!hdata->TestBit(AliQAv1::GetClonedBit())) continue; // not an array of clones, unknown object
410 TIter nextCl((TObjArray*)hdata); // array of histo clones
411 TObject* hcl = 0;
412 while ((hcl=nextCl())) if (hcl->InheritsFrom(TH1::Class()) && hcl->TestBit(AliQAv1::GetImageBit())) {tmpArr.AddLast(hcl); nImages++;}
413 }
414 //
415 if ( nImages == 0 ) {
416 AliDebug(AliQAv1::GetQADebugLevel(), Form("No histogram will be plotted for %s %s %s\n", GetName(), AliQAv1::GetTaskName(task).Data(), AliRecoParam::GetEventSpecieName(esIndex)));
417 continue;
634696f5 418 }
92664bc8 419 AliDebug(AliQAv1::GetQADebugLevel(), Form("%d histograms will be plotted for %s %s %s\n", nImages, GetName(), AliQAv1::GetTaskName(task).Data(),AliRecoParam::GetEventSpecieName(esIndex)));
420 //
421 const Char_t * title = Form("QA_%s_%s_%s", GetName(), AliQAv1::GetTaskName(task).Data(), AliRecoParam::GetEventSpecieName(esIndex));
422 //
423 if ( !fImage[esIndex] ) fImage[esIndex] = new TCanvas(title, title);
424 //
425 fImage[esIndex]->Clear();
426 fImage[esIndex]->SetTitle(title);
427 fImage[esIndex]->cd();
428 TPaveText someText(0.015, 0.015, 0.98, 0.98);
429 someText.AddText(title);
430 someText.Draw();
8386cd34 431 fImage[esIndex]->Print(Form("%s%s%d.%s", AliQAv1::GetImageFileName(), AliQAv1::GetModeName(mode), AliQAChecker::Instance()->GetRunNumber(), AliQAv1::GetImageFileFormat()));
92664bc8 432 fImage[esIndex]->Clear();
433 Int_t nx = TMath::Nint(TMath::Sqrt(nImages));
434 Int_t ny = nx;
435 if (nx < TMath::Sqrt(nImages)) ny++;
436 //
437 fImage[esIndex]->Divide(nx, ny);
438 TIter nexthist(&tmpArr);
439 Int_t npad = 1;
440 fImage[esIndex]->cd(npad);
441 TH1* histo = 0;
442 while ( (histo=(TH1*)nexthist()) ) { // tmpArr is guaranteed to contain only plottable histos, no checks needed
443 TString opts = histo->GetDrawOption();
444 if (opts.Contains("logy",TString::kIgnoreCase)) gPad->SetLogy();
445 if (opts.Contains("logx",TString::kIgnoreCase)) gPad->SetLogx();
446 histo->DrawCopy();
447 fImage[esIndex]->cd(++npad);
448 }
449 fImage[esIndex]->Print(Form("%s%s%d.%s", AliQAv1::GetImageFileName(), AliQAv1::GetModeName(mode), AliQAChecker::Instance()->GetRunNumber(), AliQAv1::GetImageFileFormat()), "ps");
450 }
634696f5 451}
452
46b8a84d 453//____________________________________________________________________________
57acd2d2 454void AliQACheckerBase::SetHiLo(Float_t * hiValue, Float_t * lowValue)
46b8a84d 455{
5379c4a3 456 AliDebug(AliQAv1::GetQADebugLevel(), "Previous setting was:") ;
457 if ( AliDebugLevel() == AliQAv1::GetQADebugLevel() ) {
eca4fa66 458 const Char_t * text= Form(" INFO -> %1.5f < value < %1.5f WARNING -> %1.5f < value <= %1.5f \n ERROR -> %1.5f < value <= %1.5f \n FATAL -> %1.5f <= value < %1.5f \n",
459 fLowTestValue[AliQAv1::kINFO], fUpTestValue[AliQAv1::kINFO],
460 fLowTestValue[AliQAv1::kWARNING], fUpTestValue[AliQAv1::kWARNING],
461 fLowTestValue[AliQAv1::kERROR], fUpTestValue[AliQAv1::kERROR],
462 fLowTestValue[AliQAv1::kFATAL], fUpTestValue[AliQAv1::kFATAL]) ;
463 AliInfo(Form("%s", text)) ;
5379c4a3 464 }
57acd2d2 465
4e25ac79 466 for (Int_t index = 0 ; index < AliQAv1::kNBIT ; index++) {
57acd2d2 467 fLowTestValue[index] = lowValue[index] ;
eca4fa66 468 fUpTestValue[index] = hiValue[index] ;
57acd2d2 469 }
5379c4a3 470 AliDebug(AliQAv1::GetQADebugLevel(), "Current setting is:") ;
471 if ( AliDebugLevel() == AliQAv1::GetQADebugLevel() ) {
eca4fa66 472 const Char_t * text= Form(" INFO -> %1.5f < value < %1.5f WARNING -> %1.5f < value <= %1.5f \n ERROR -> %1.5f < value <= %1.5f \n FATAL -> %1.5f <= value < %1.5f \n",
473 fLowTestValue[AliQAv1::kINFO], fUpTestValue[AliQAv1::kINFO],
474 fLowTestValue[AliQAv1::kWARNING], fUpTestValue[AliQAv1::kWARNING],
475 fLowTestValue[AliQAv1::kERROR], fUpTestValue[AliQAv1::kERROR],
476 fLowTestValue[AliQAv1::kFATAL], fUpTestValue[AliQAv1::kFATAL]) ; AliInfo(Form("%s", text)) ;
5379c4a3 477 }
57acd2d2 478}
46b8a84d 479
57acd2d2 480//____________________________________________________________________________
4e25ac79 481void AliQACheckerBase::SetQA(AliQAv1::ALITASK_t index, Double_t * value) const
57acd2d2 482{
483 // sets the QA according the return value of the Check
46b8a84d 484
4e25ac79 485 AliQAv1 * qa = AliQAv1::Instance(index) ;
d5b176f5 486
57acd2d2 487 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
d5b176f5 488 if (! qa->IsEventSpecieSet(AliRecoParam::ConvertIndex(specie)))
489 continue ;
57acd2d2 490 if ( value == NULL ) { // No checker is implemented, set all QA to Fatal
4e25ac79 491 qa->Set(AliQAv1::kFATAL, specie) ;
57acd2d2 492 } else {
4e25ac79 493 if ( value[specie] >= fLowTestValue[AliQAv1::kFATAL] && value[specie] < fUpTestValue[AliQAv1::kFATAL] )
b872fee3 494 qa->Set(AliQAv1::kFATAL, AliRecoParam::ConvertIndex(specie)) ;
4e25ac79 495 else if ( value[specie] > fLowTestValue[AliQAv1::kERROR] && value[specie] <= fUpTestValue[AliQAv1::kERROR] )
b872fee3 496 qa->Set(AliQAv1::kERROR, AliRecoParam::ConvertIndex(specie)) ;
4e25ac79 497 else if ( value[specie] > fLowTestValue[AliQAv1::kWARNING] && value[specie] <= fUpTestValue[AliQAv1::kWARNING] )
b872fee3 498 qa->Set(AliQAv1::kWARNING, AliRecoParam::ConvertIndex(specie)) ;
4e25ac79 499 else if ( value[specie] > fLowTestValue[AliQAv1::kINFO] && value[specie] <= fUpTestValue[AliQAv1::kINFO] )
b872fee3 500 qa->Set(AliQAv1::kINFO, AliRecoParam::ConvertIndex(specie)) ;
57acd2d2 501 }
502 }
46b8a84d 503}