]> git.uio.no Git - u/mrichter/AliRoot.git/blame - T0/AliT0QAChecker.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
ana update salvatore
[u/mrichter/AliRoot.git] / T0 / AliT0QAChecker.cxx
CommitLineData
c03351ac 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 **************************************************************************/
c03351ac 15//...
16// Checks the quality assurance.
17// By comparing with reference data
18// Skeleton for T0
f6c2d5e2 19//---------------------------------------------
20//checkig without reference data:
21//for RAW QA all histograms should have approximatly the same
22//number of entries as RefPoint
23//for Rec Points checks
24// - amplitude measured by 2 methos
25// - online and offline T0 measurements
26// for ESD quality of reconstruction ( and measurements):
27// RMS of vertex and T0 less than 75ps
28//
29// Alla.Maevskaya@cern.ch
c03351ac 30//...
31
32// --- ROOT system ---
b09247a2 33#include <Riostream.h>
c03351ac 34#include <TClass.h>
35#include <TH1F.h>
be3ada9f 36#include <TF1.h>
37#include <TFitResultPtr.h>
0d849919 38#include <TH2.h>
c03351ac 39#include <TIterator.h>
40#include <TKey.h>
41#include <TFile.h>
446d6ec4 42#include <TMath.h>
0d849919 43#include <TString.h>
be3ada9f 44#include <TPaveText.h>
c03351ac 45
46// --- Standard library ---
47
48// --- AliRoot header files ---
49#include "AliLog.h"
4e25ac79 50#include "AliQAv1.h"
c03351ac 51#include "AliQAChecker.h"
394c1a6d 52#include "AliCDBEntry.h"
53#include "AliQAManager.h"
c03351ac 54#include "AliT0QAChecker.h"
55
56ClassImp(AliT0QAChecker)
be3ada9f 57//____________________________________________________________________________
58AliT0QAChecker::AliT0QAChecker() :
59AliQACheckerBase("T0","T0 Quality Assurance Checker")
60
61{
62 // Standard constructor
63
64}
65
49e33ff1 66//____________________________________________________________________________
67AliT0QAChecker::AliT0QAChecker(const AliT0QAChecker& qac):
68 AliQACheckerBase(qac.GetName(), qac.GetTitle())
69{
70 // copy constructor
71 AliError("Copy should not be used with this class\n");
72}
73//____________________________________________________________________________
74AliT0QAChecker& AliT0QAChecker::operator=(const AliT0QAChecker& qac){
75 // assignment operator
76 this->~AliT0QAChecker();
77 new(this)AliT0QAChecker(qac);
78 return *this;
79}
80
81
be3ada9f 82//____________________________________________________________________________
83AliT0QAChecker::~AliT0QAChecker(){
84 // destructor
85
86}
c03351ac 87
4e67e173 88//__________________________________________________________________
a42ceb0e 89void AliT0QAChecker::Check(Double_t * test, AliQAv1::ALITASK_t index, TObjArray ** list, const AliDetectorRecoParam * /*recoParam*/)
3685c21c 90{
3685c21c 91 // Super-basic check on the QA histograms on the input list:
92 // look whether they are empty!
a42ceb0e 93
be3ada9f 94 char * detOCDBDir = Form("T0/%s/%s", AliQAv1::GetRefOCDBDirName(), AliQAv1::GetRefDataDirName()) ;
8825d145 95
394c1a6d 96 AliCDBEntry *QARefRec = AliQAManager::QAManager()->Get(detOCDBDir);
97 // QARefRec->Dump();
98 if( !QARefRec){
99 AliInfo("QA reference data NOT retrieved for Reconstruction check. No T0 reference distribution");
100 }
c724dd64 101
49e33ff1 102
103 for(Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++){
104 test[specie] = 1.0; //FK// initiate qa flag for the whole set of histograms as good
105 }
106
107
108 for(Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
109 if(!(AliQAv1::Instance()->IsEventSpecieSet(specie) && list[specie]) || list[specie]->GetEntries() == 0) {
92664bc8 110 continue;
57acd2d2 111 }
1271d24f 112
49e33ff1 113 if(index == AliQAv1::kRAW){
114
115 if(AliRecoParam::ConvertIndex(specie) == AliRecoParam::kCalib){// if (index == AliQAv1::kRAW )
116 //check laser data efficiencies
117 Double_t qaFlag = CheckLaser(list[specie]);
118 if(qaFlag < test[specie]) test[specie] = qaFlag;
119 }
120
121 if(AliRecoParam::ConvertIndex(specie) == AliRecoParam::kCalib ||
7b51ee39 122 AliRecoParam::ConvertIndex(specie) == AliRecoParam::kHighMult ||
1271d24f 123 AliRecoParam::ConvertIndex(specie) == AliRecoParam::kLowMult){
124 //AliRecoParam::ConvertIndex(specie) == AliRecoParam::kDefault ||
49e33ff1 125
126 //check BCID
127 Double_t qaFlag = CheckBCID(list[specie]);
128 if(qaFlag < test[specie]) test[specie] = qaFlag;
129 }
130
7b51ee39 131 if(AliRecoParam::ConvertIndex(specie) == AliRecoParam::kHighMult ||
1271d24f 132 AliRecoParam::ConvertIndex(specie) == AliRecoParam::kLowMult){
133 //AliRecoParam::ConvertIndex(specie) == AliRecoParam::kDefault ||
49e33ff1 134 //check physics
135 Double_t qaFlag = CheckRaw(list[specie]);
136 if(qaFlag < test[specie]) test[specie] = qaFlag;
92664bc8 137 }
49e33ff1 138 }
139
140 if(index == AliQAv1::kESD && AliRecoParam::Convert(specie) != AliRecoParam::kCalib){
92664bc8 141 test[specie] = CheckESD(list[specie]);
49e33ff1 142 }
be3ada9f 143 }
144}
394c1a6d 145
be3ada9f 146//--------------------------------------------------------------------------
49e33ff1 147Double_t AliT0QAChecker::CheckLaser(TObjArray *listrec) const {
be3ada9f 148
49e33ff1 149 TH1 *hdata;
150 TH1 *fhRawEff[10];
151 Int_t nEffHistos=0;
152
153 //thresholds for warning and error on efficiencies
154 Float_t thrWarning = 0.5; //FK// warning level
155 Float_t thrError = 0.2; //FK// error level
156
157 const int kNumberOfHistos = 3;
158 Int_t consecutiveHistoNumber[kNumberOfHistos] = { 207, 208, 209}; //Checked histos fhCDFeff, hEffLED, hEffQTC
159 Int_t qualityFlag[kNumberOfHistos]; //quality flag for a given histogram
160
161 for(Int_t ir=0; ir<kNumberOfHistos; ir++){
162 qualityFlag[ir] = kT0Info; //init quality flag for a given histogram
163
164 hdata = (TH1*) listrec->UncheckedAt(consecutiveHistoNumber[ir]);
165 if(hdata){
166 fhRawEff[nEffHistos] = hdata;
167 nEffHistos++;
168 }
169 }
170
171 TLine linelowyellow(0, thrWarning, 24, thrWarning);
172 linelowyellow.SetLineColor(5);
173 linelowyellow.SetLineStyle(3);
174 linelowyellow.SetLineWidth(4);
175 TLine linelowred(0, thrError, 24, thrError);
176 linelowred.SetLineColor(2);
177 linelowred.SetLineStyle(3);
178 linelowred.SetLineWidth(4);
179
180 Bool_t bEffHistosNotEmpty = kFALSE; //check if all histograms have some counts
181
182 for(Int_t ih = 0; ih < nEffHistos; ih++){
183
184 EraseOldMessages((TH1*) fhRawEff[ih]);// clean objects added at previous checks
be3ada9f 185
49e33ff1 186 fhRawEff[ih]->SetLineWidth(2);
187 fhRawEff[ih]->SetMaximum(2.);
188 fhRawEff[ih]->SetMinimum(0.);
189 fhRawEff[ih]->GetListOfFunctions()->Add((TLine*)linelowyellow.Clone());
190 fhRawEff[ih]->GetListOfFunctions()->Add((TLine*)linelowred.Clone());
394c1a6d 191
49e33ff1 192 if(fhRawEff[ih]->Integral()>0) bEffHistosNotEmpty = kTRUE; //this histo does have some counts in it
193
194 Int_t nbins= fhRawEff[ih]->GetNbinsX();
195 for(Int_t ib=1; ib<=nbins; ib++){ //loop over bins and check if the efficiency is above level
196
197 Float_t chcont = fhRawEff[ih]->GetBinContent(ib);
198 if(chcont < thrWarning && qualityFlag[ih] > kT0Error ) qualityFlag[ih] = kT0Warning;//Warning level
199 if(chcont < thrError) qualityFlag[ih] = kT0Error;//Error level
be3ada9f 200 }
49e33ff1 201
202 if(qualityFlag[ih] == kT0Info ){
203 AliDebug(AliQAv1::GetQADebugLevel(), Form("T0 efficiency %s is good", fhRawEff[ih]->GetName() ));
204 }else if(qualityFlag[ih] == kT0Warning){
205 AliDebug(AliQAv1::GetQADebugLevel(), Form("T0 efficiency %s is not so good", fhRawEff[ih]->GetName() ));
206 }else if(qualityFlag[ih] == kT0Error){
207 AliDebug(AliQAv1::GetQADebugLevel(), Form("T0 efficiency %s is not good", fhRawEff[ih]->GetName() ));
208 }
209 }
210
211 //executive summary
212 int lowestQualityFlag = (int) kT0Info;
213 for(Int_t ih = 0; ih < nEffHistos; ih++){
214
215 if(!bEffHistosNotEmpty){ //all laser efficiency plots are empty
7b51ee39 216
49e33ff1 217 TPaveText text(0.20,0.50,0.99,0.99,"NDC");
218 text.AddText(Form("1) T0 is in BEAMTUNIG: empty plots are ok"));
219 text.AddText(Form("2) T0 is in READY: check calibriation trigger"));
220 text.AddText(Form("if also physics data are empty report"));
221 text.AddText(Form("readout problem to the T0 on-call expert"));
7b51ee39 222 text.SetBorderSize(0);
223 text.SetFillStyle(0);
49e33ff1 224 fhRawEff[ih]->GetListOfFunctions()->Add((TPaveText*)text.Clone());
225 }
226
227 if( qualityFlag[ih] <lowestQualityFlag ) lowestQualityFlag = qualityFlag[ih];
228 }
229
230 return ConvertQualityFlagToDouble(lowestQualityFlag);
231}
232//--------------------------------------------------------------------------
233Double_t AliT0QAChecker::CheckBCID(TObjArray *listrec) const {
234
235 Int_t qualityFlagBCID = kT0Info; //init quality flag for a given histogram;
236
237 TH2F *hBCID = (TH2F*) listrec->UncheckedAt(224); //BCID versus TRM BCID
238
be3ada9f 239
49e33ff1 240 // clean objects added at previous checks
241 EraseOldMessages((TH1*)hBCID);
242
243 if(hBCID->Integral()>0){
244 //BCID does have some counts in it
245
246 Int_t nbinsX = hBCID->GetNbinsX();
247 Int_t nbinsY = hBCID->GetNbinsY();
248 double entriesOnDiagonal = 0; //count diagonal and off diagonal entries
249 double entriesOffDiagonal = 0;
7b51ee39 250 Int_t offset = 37;
251 Int_t period = 3564;
252
253 for(Int_t itrm=1; itrm<=nbinsX; itrm++){ //BCID TRM
254 for(Int_t ibcid=1; ibcid<=nbinsY; ibcid++){ //BCID
255 if(TMath::Abs(itrm-ibcid)<6 || TMath::Abs( (itrm+offset)%period - ibcid)<2 ){
256 entriesOnDiagonal += hBCID->GetBinContent(itrm,ibcid); //On Diagonal
257 }else{
258 entriesOffDiagonal += hBCID->GetBinContent(itrm,ibcid); //Off Diagonal
259 }
57acd2d2 260 }
49e33ff1 261 }
7b51ee39 262 if(entriesOnDiagonal<1 || entriesOffDiagonal>20){
263 qualityFlagBCID = kT0Warning; //no entries on diagonal
49e33ff1 264 AliDebug(AliQAv1::GetQADebugLevel(), Form("T0 %s is not diagonal", hBCID->GetName() ));
265
266 TPaveText text(0.20,0.50,0.99,0.99,"NDC");
267 text.AddText(Form("Check if entries are on a diagonal."));
268 text.AddText(Form("Report readout problem to the T0 on-call expert"));
7b51ee39 269 text.SetBorderSize(0);
270 text.SetFillStyle(0);
49e33ff1 271 hBCID->GetListOfFunctions()->Add((TPaveText*)text.Clone());
272 }
273 }else{ //BCID empty
274
275 qualityFlagBCID = kT0Error;
276 AliDebug(AliQAv1::GetQADebugLevel(), Form("T0 : %s has NO entries", hBCID->GetName() ));
277
278 TPaveText text(0.20,0.50,0.99,0.99,"NDC");
279 text.AddText(Form("NO ENTRIES!!!"));
280 text.AddText(Form("If T0 is READY report"));
7b51ee39 281 text.AddText(Form("readout problem to the T0 on-call expert"));
282 text.SetBorderSize(0);
283 text.SetFillStyle(0);
49e33ff1 284 hBCID->GetListOfFunctions()->Add((TPaveText*)text.Clone());
285 }
286
287 //executive summary
288 int lowestQualityFlag = (int) qualityFlagBCID;
289
290 return ConvertQualityFlagToDouble(lowestQualityFlag);
be3ada9f 291
3685c21c 292}
5ed41460 293
49e33ff1 294//--------------------------------------------------------------------------
295Double_t AliT0QAChecker::CheckRaw(TObjArray *listrec) const {
296
297
298 Int_t qualityFlagTrigger = kT0Info; //init quality flag for a given histogram;
299
300 TH1F *hTrigger = (TH1F*) listrec->UncheckedAt(169);//hRawTrigger
301
302
303 // clean objects added at previous checks
304 EraseOldMessages((TH1*) hTrigger);
305
306 if(hTrigger->Integral()>0){
307 //trigger plot does have some counts in it
308 //are Mean, ORA and ORC not empty?
309 if( hTrigger->GetBinContent(1)<0.001 || hTrigger->GetBinContent(3)<0.001 || hTrigger->GetBinContent(4)<0.001){
310 qualityFlagTrigger = kT0Error; //no entries on diagonal
311 AliDebug(AliQAv1::GetQADebugLevel(), Form("T0: too little ORA and ORC in %s", hTrigger->GetName() ));
312
313 TPaveText text(0.20,0.50,0.99,0.99,"NDC");
314 text.AddText(Form("Check ORA and ORC"));
315 text.AddText(Form("Report problem to the T0 on-call expert"));
7b51ee39 316 text.SetBorderSize(0);
317 text.SetFillStyle(0);
49e33ff1 318 hTrigger->GetListOfFunctions()->Add((TPaveText*)text.Clone());
319 }
1271d24f 320
49e33ff1 321 }else{ //Trigger histo empty
322
323 qualityFlagTrigger = kT0Error;
324 AliDebug(AliQAv1::GetQADebugLevel(), Form("T0 histogram %s has NO entries", hTrigger->GetName() ));
325
326 TPaveText text(0.20,0.50,0.99,0.99,"NDC");
327 text.AddText(Form("NO ENTRIES!!!"));
328 text.AddText(Form("If T0 is READY report"));
329 text.AddText(Form("readout problem to the T0 on-call expert"));
7b51ee39 330 text.SetBorderSize(0);
331 text.SetFillStyle(0);
49e33ff1 332 hTrigger->GetListOfFunctions()->Add((TPaveText*)text.Clone());
1271d24f 333
49e33ff1 334 }
335
1271d24f 336 //--------- timing plots ---------
337 Int_t qualityFlagMeanBest = kT0Info; //init quality flag for a given histogram;
338 TH1F *hMeanBest = (TH1F*) listrec->UncheckedAt(223); //fhMeanBest time
339 EraseOldMessages((TH1*) hMeanBest);
340
341 if(hMeanBest->Integral()<1){
342 qualityFlagMeanBest = kT0Error;
343 AliDebug(AliQAv1::GetQADebugLevel(), Form("T0 histogram %s has NO entries", hMeanBest->GetName() ));
344
345 TPaveText text(0.20,0.50,0.99,0.99,"NDC");
346 text.AddText(Form("NO ENTRIES!!!"));
347 text.AddText(Form("If T0 is READY and beam is on report"));
7b51ee39 348 text.AddText(Form("the problem to the T0 on-call expert"));
349 text.SetBorderSize(0);
350 text.SetFillStyle(0);
1271d24f 351 hMeanBest->GetListOfFunctions()->Add((TPaveText*)text.Clone());
352 }
353
354
355
356
357 //----------- vetext plots ---------
358
359 Int_t qualityFlagVetrexFirst = kT0Info; //init quality flag for a given histogram;
360 Int_t qualityFlagVetrexBest = kT0Info; //init quality flag for a given histogram;
361
362 TH1F *hVertexFirst = (TH1F*) listrec->UncheckedAt(225);//fhVertex1st
363 TH1F *hVertexBest = (TH1F*) listrec->UncheckedAt(226);//fhVertexBest
364 TH1F *hOrCminOrATvdcOn = (TH1F*) listrec->UncheckedAt(217);//vertex with TVDC on
365
366 // clean objects added at previous checks
367 EraseOldMessages((TH1*) hVertexFirst);
368 EraseOldMessages((TH1*) hVertexBest);
7b51ee39 369 float errorLevelDifference = 400;
370 float warningLevelDifference = 200;
1271d24f 371
372 if(hVertexFirst->Integral()>0 && hOrCminOrATvdcOn->Integral()>0){
7b51ee39 373 Float_t thr = 0.05 * hVertexFirst->GetBinContent(hVertexFirst->GetMaximumBin());
374
375 float meanVertexFirst = GetMeanAboveThreshold( hVertexFirst, thr);
376 float meanVertexTVDC = GetMeanAboveThreshold( hOrCminOrATvdcOn,thr);
1271d24f 377 float diff = TMath::Abs(meanVertexFirst - meanVertexTVDC);
378 if(diff > errorLevelDifference){
7b51ee39 379 qualityFlagVetrexFirst = kT0Warning; //init quality flag for a given histogram;
1271d24f 380
381 TPaveText text(0.20,0.50,0.99,0.99,"NDC");
7b51ee39 382 text.AddText(Form("large diff. between TVDC vertex"));
383 text.AddText(Form("and the first vertex."));
384 text.SetBorderSize(0);
385 text.SetFillStyle(0);
1271d24f 386 hVertexFirst->GetListOfFunctions()->Add((TPaveText*)text.Clone());
387
388 }else if(diff > warningLevelDifference){
389 qualityFlagVetrexFirst = kT0Warning; //init quality flag for a given histogram;
390
391 TPaveText text(0.20,0.50,0.99,0.99,"NDC");
7b51ee39 392 text.AddText(Form("large diff. between TVDC vertex"));
393 text.AddText(Form("and the first vertex."));
394 text.SetBorderSize(0);
395 text.SetFillStyle(0);
1271d24f 396 hVertexFirst->GetListOfFunctions()->Add((TPaveText*)text.Clone());
397
398 }
399 }else{
400 qualityFlagVetrexFirst = kT0Error;
401 AliDebug(AliQAv1::GetQADebugLevel(), Form("T0 histogram %s has NO entries", hVertexFirst->GetName() ));
402
403 TPaveText text(0.20,0.50,0.99,0.99,"NDC");
404 text.AddText(Form("NO ENTRIES!!!"));
405 text.AddText(Form("If T0 is READY and beam is on report"));
7b51ee39 406 text.AddText(Form("the problem to the T0 on-call expert"));
407 text.SetBorderSize(0);
408 text.SetFillStyle(0);
1271d24f 409 hVertexFirst->GetListOfFunctions()->Add((TPaveText*)text.Clone());
410 }
411
412 if(hVertexBest->Integral()>0 && hOrCminOrATvdcOn->Integral()>0){
7b51ee39 413
414 Float_t thr = 0.05 * hVertexBest->GetBinContent(hVertexBest->GetMaximumBin());
415 float meanVertexBest = GetMeanAboveThreshold(hVertexBest,thr);
416 float meanVertexTVDC = GetMeanAboveThreshold(hOrCminOrATvdcOn,thr);
1271d24f 417 float diff = TMath::Abs(meanVertexBest - meanVertexTVDC);
418 if(diff > errorLevelDifference){
7b51ee39 419 qualityFlagVetrexBest = kT0Warning; //init quality flag for a given histogram;
1271d24f 420
421 TPaveText text(0.20,0.50,0.99,0.99,"NDC");
7b51ee39 422 text.AddText(Form("Large. diff. between TVDC vertex"));
423 text.AddText(Form("and the best vertex."));
424 text.SetBorderSize(0);
425 text.SetFillStyle(0);
1271d24f 426 hVertexBest->GetListOfFunctions()->Add((TPaveText*)text.Clone());
427
428 }else if(diff > warningLevelDifference){
429 qualityFlagVetrexBest = kT0Warning; //init quality flag for a given histogram;
430
431 TPaveText text(0.20,0.50,0.99,0.99,"NDC");
7b51ee39 432 text.AddText(Form("Large. diff. between TVDC vertex"));
433 text.AddText(Form("and the best vertex."));
434 text.SetBorderSize(0);
435 text.SetFillStyle(0);
1271d24f 436 hVertexBest->GetListOfFunctions()->Add((TPaveText*)text.Clone());
437 }
438 }else{
439 qualityFlagVetrexBest = kT0Error;
440 AliDebug(AliQAv1::GetQADebugLevel(), Form("T0 histogram %s has NO entries", hVertexBest->GetName() ));
441
442 TPaveText text(0.20,0.50,0.99,0.99,"NDC");
443 text.AddText(Form("NO ENTRIES!!!"));
444 text.AddText(Form("If T0 is READY and beam is on report"));
7b51ee39 445 text.AddText(Form("the problem to the T0 on-call expert"));
446 text.SetBorderSize(0);
447 text.SetFillStyle(0);
1271d24f 448 hVertexBest->GetListOfFunctions()->Add((TPaveText*)text.Clone());
449 }
450
451
49e33ff1 452 //executive summary
453 int lowestQualityFlag = (int) qualityFlagTrigger;
1271d24f 454 if(qualityFlagVetrexBest<lowestQualityFlag) lowestQualityFlag = qualityFlagVetrexBest;
455 if(qualityFlagVetrexFirst<lowestQualityFlag) lowestQualityFlag = qualityFlagVetrexFirst;
456 if(qualityFlagMeanBest <lowestQualityFlag) lowestQualityFlag = qualityFlagMeanBest;
457
49e33ff1 458
459 return ConvertQualityFlagToDouble(lowestQualityFlag);
460
461}
be3ada9f 462
394c1a6d 463//--------------------------------------------------------------------------
be3ada9f 464Double_t AliT0QAChecker::CheckESD(TObjArray *listrec ) const
394c1a6d 465{
394c1a6d 466 Float_t checkr = 0;
be3ada9f 467 TH1 *fhESD;
468
469 fhESD = (TH1*) listrec->UncheckedAt(2);
470 if(fhESD){
471 AliDebug(AliQAv1::GetQADebugLevel(), Form("count %s ", fhESD->GetName()) );
472 TF1 *f1 = new TF1("f1","gaus",-1,1);
473 fhESD->Fit("f1","R","Q", -1,1);
474 Double_t par[3];
475 f1->GetParameters(&par[0]);
476
105ff51c 477 TPaveText text(0.30,0.50,0.99,0.99,"NDC");
be3ada9f 478
105ff51c 479 text.AddText(Form("T0 RUN %d ",AliCDBManager::Instance()->GetRun()));
be3ada9f 480
db33b25a 481 AliDebug(AliQAv1::GetQADebugLevel(), Form("numentries %d mean %f #sigma %f", (int)fhESD->GetEntries(),par[1], par[2]));
be3ada9f 482
483
484 if (par[2] > 0.07 && par[2] < 1.) {
485 checkr=0.5;
105ff51c 486 text.AddText(Form("not good resolution :\n %f ns\n", par[2] ));
487 text.SetFillColor(5);
be3ada9f 488 printf("T0 detector resolution is not good enouph: %f ns\n",par[2] );
489 }
490 if(TMath::Abs(par[1])>0.05) {
491 checkr = 0.5;
105ff51c 492 text.AddText(Form(" Check clock shift on %f ns", par[1]));
493 text.SetFillColor(5);
be3ada9f 494 }
495 if (par[2] > 1. || TMath::Abs(par[1])>0.1) {
496 checkr = 0.25;
105ff51c 497 text.AddText(Form(" Bad resolution:\n mean %f ns sigma %f ns", par[1], par[2]));
498 text.SetFillColor(2);
499 { // RS Clean previous additions
500 TList* lstF = fhESD->GetListOfFunctions();
501 if (lstF) {
502 TObject *stats = lstF->FindObject("stats");
503 lstF->Remove(stats);
504 TObject *obj;
505 while ((obj = lstF->First())) {
506 while(lstF->Remove(obj)) { }
507 delete obj;
508 }
509 if (stats) lstF->Add(stats);
510 }
511 }
512 fhESD->GetListOfFunctions()->Add((TPaveText*)text.Clone());
513 AliDebug(AliQAv1::GetQADebugLevel(),
514 Form("Please, check calibration: shift= %f resolution %f test=%f\n",
515 par[1], par[2], checkr) ) ;
be3ada9f 516 }
394c1a6d 517 }
be3ada9f 518 else
519 {
520 AliDebug(AliQAv1::GetQADebugLevel(),
521 Form("No ESD QA histogram found, nothing to check"));
522 checkr=0;
523 }
524
332b9234 525
394c1a6d 526 return checkr;
527}
49e33ff1 528
529
530//--------------------------------------------------------------------------
531void AliT0QAChecker::EraseOldMessages(TH1* h) const
532{
533 //erase the old captions
534 TList* lstF = h->GetListOfFunctions();
535 if(lstF){
536 TObject *stats = lstF->FindObject("stats");
537 lstF->Remove(stats);
538 TObject *obj;
539 while ((obj = lstF->First())) {
540 while(lstF->Remove(obj)) { }
541 delete obj;
542 }
543 if (stats) lstF->Add(stats);
544 }
545}
546//--------------------------------------------------------------------------
547Double_t AliT0QAChecker::ConvertQualityFlagToDouble(int qualityFlag) const
548{
549 //covert quality flag to double
550 Double_t checkr=1.0;
551
552 switch ( qualityFlag ){
553 case kT0Info:
554 checkr = 1.0; break;
555 case kT0Warning:
556 checkr = 0.75; break;
557 case kT0Error:
558 checkr = 0.25; break;
559 case kT0Fatal:
560 checkr = -1.0; break;
561 default:
562 AliError("Invalid ecc value. FIXME !");
563 checkr = 0.25; break;
564 };
565
566 return checkr;
7b51ee39 567}
568
569//--------------------------------------------------------------------------
570Float_t AliT0QAChecker::GetMeanAboveThreshold(TH1F* hV, Float_t thr) const{
571 //caculate mean value of histo bins above threshold
572 Int_t nBins = hV->GetNbinsX();
573 Int_t nBinsAboveThr = 0;
574 Float_t sum = 0;
575
576 for(Int_t ib=1;ib<=nBins;ib++){
577 Float_t val = hV->GetBinContent(ib);
578 if(val<thr){
579 sum+=val;
580 nBinsAboveThr++;
581 }
582 }
583
584 if(nBinsAboveThr>0) return sum/nBinsAboveThr;
585 else return hV->GetMean();
586}