]>
Commit | Line | Data |
---|---|---|
075a0e70 | 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 | // --- ROOT system --- | |
075a0e70 | 18 | #include <TH1F.h> |
075a0e70 | 19 | #include <TIterator.h> |
fdc38bb2 | 20 | #include <TString.h> |
891ddd89 | 21 | #include <TPaveText.h> |
22 | #include <TObjArray.h> | |
23 | #include "TList.h" | |
075a0e70 | 24 | |
25 | // --- Standard library --- | |
26 | ||
27 | // --- AliRoot header files --- | |
28 | #include "AliLog.h" | |
075a0e70 | 29 | #include "AliZDCQAChecker.h" |
30 | ||
31 | ClassImp(AliZDCQAChecker) | |
32 | ||
075a0e70 | 33 | //____________________________________________________________________________ |
a42ceb0e | 34 | void AliZDCQAChecker::Check(Double_t * test, AliQAv1::ALITASK_t index, TObjArray ** list, |
486788fc | 35 | const AliDetectorRecoParam * /*recoParam*/) |
075a0e70 | 36 | { |
57acd2d2 | 37 | // Checks the QA histograms on the input list: |
38 | // | |
4e25ac79 | 39 | const char* taskName = AliQAv1::GetAliTaskName(index); |
2d9c70ab | 40 | //printf("\n\tAliZDCQAChecker -> checking QA histos for task %s\n",taskName); |
57acd2d2 | 41 | // |
c473b43d | 42 | Int_t ihitHisto=0, idigHisto=0; |
43 | Int_t irecHisto=0, irawHisto=0, esdInd=0; | |
44 | ||
891ddd89 | 45 | TObjArray messages; |
46 | messages.SetOwner(kTRUE); | |
47 | ||
f5e4ee59 | 48 | for(Int_t specie = 0; specie<AliRecoParam::kNSpecies; specie++){ |
a42ceb0e | 49 | Int_t count = 0; |
c473b43d | 50 | if(!AliQAv1::Instance(AliQAv1::GetDetIndex(GetName()))->IsEventSpecieSet(AliRecoParam::ConvertIndex(specie))) |
51 | continue ; | |
52 | //printf("\tAliZDCQAChecker -> specie %d, AliRecoParam::ConvertIndex(specie) %d, AliRecoParam::kLowMult %d, IsEventSpecieSet(specie) %d\n", | |
53 | // specie, AliRecoParam::ConvertIndex(specie) ,AliRecoParam::kLowMult, | |
54 | // AliQAv1::Instance(AliQAv1::GetDetIndex(GetName()))->IsEventSpecieSet(AliRecoParam::ConvertIndex(specie))); | |
f5e4ee59 | 55 | |
57acd2d2 | 56 | // ==================================================================== |
57 | // Checks for p-p events | |
58 | // ==================================================================== | |
c473b43d | 59 | if(AliRecoParam::ConvertIndex(specie) == AliRecoParam::kLowMult){ |
57acd2d2 | 60 | if(list[specie]->GetEntries()==0){ |
213a7d34 | 61 | AliWarning("\t The list to be checked is empty!"); // nothing to check |
c473b43d | 62 | return; |
12b9729c | 63 | } |
5379c4a3 | 64 | //AliDebug(AliQAv1::GetQADebugLevel(), Form("\n\tAliZDCQAChecker-> checking QA histograms for task %s\n\n",taskName)); |
57acd2d2 | 65 | TIter next(list[specie]); |
a42ceb0e | 66 | count = 0; |
57acd2d2 | 67 | TH1 * hdata; |
c473b43d | 68 | |
331464fb | 69 | Float_t res=0., percentageDiff=0.30; |
c473b43d | 70 | Float_t meanZNA=0., meanZNC=0., meanZPA=0., meanZPC=0.; |
71 | Float_t pmCZNA=0., pmCZNC=0., pmCZPA=0., pmCZPC=0.; | |
72 | Float_t pmQZNA=0., pmQZNC=0., pmQZPA=0., pmQZPC=0.; | |
73 | Float_t sumADCZNA=0., sumADCZNC=0., sumADCZPA=0., sumADCZPC=0.; | |
74 | Float_t adcCZNA=0., adcCZNC=0., adcCZPA=0., adcCZPC=0.; | |
75 | Float_t adcQZNA=0., adcQZNC=0., adcQZPA=0., adcQZPC=0.; | |
4998d47e | 76 | Int_t nentries = -99; |
c473b43d | 77 | |
57acd2d2 | 78 | while((hdata = dynamic_cast<TH1 *>(next()))){ |
79 | if(hdata){ | |
80 | // ------------------------------------------------------------------- | |
4e25ac79 | 81 | if(index == AliQAv1::kSIM){ |
5379c4a3 | 82 | //AliDebug(AliQAv1::GetQADebugLevel(), Form("\tAliZDCQAChecker-> checking histo %s",hdata->GetName())); |
57acd2d2 | 83 | // Check HITS histos |
213a7d34 | 84 | // |
85 | if(!(strncmp(hdata->GetName(),"hHits",5))){ | |
57acd2d2 | 86 | if(hdata->GetEntries()>0){ |
c473b43d | 87 | if(ihitHisto==0) meanZNC = hdata->GetMean(); |
88 | else if(ihitHisto==1) meanZNA = hdata->GetMean(); | |
89 | else if(ihitHisto==2) meanZPC = hdata->GetMean(); | |
90 | else if(ihitHisto==3) meanZPA = hdata->GetMean(); | |
213a7d34 | 91 | else if(ihitHisto==4) pmQZNC = hdata->GetMean(); |
92 | else if(ihitHisto==5) pmQZNA = hdata->GetMean(); | |
93 | else if(ihitHisto==6) pmQZPC = hdata->GetMean(); | |
94 | else if(ihitHisto==7) pmQZPA = hdata->GetMean(); | |
95 | else if(ihitHisto==8) pmCZNC = hdata->GetMean(); | |
96 | else if(ihitHisto==9) pmCZNA = hdata->GetMean(); | |
97 | else if(ihitHisto==10) pmCZPC = hdata->GetMean(); | |
98 | else if(ihitHisto==11) pmCZPA = hdata->GetMean(); | |
99 | } | |
100 | // | |
101 | // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff | |
102 | if(ihitHisto==11){ | |
c473b43d | 103 | if(TMath::Abs(meanZNC)>1.e-10){ |
213a7d34 | 104 | if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) |
105 | res=1.; | |
106 | else | |
107 | res=.5; | |
108 | test[specie] += res; | |
a42ceb0e | 109 | count++; |
213a7d34 | 110 | } |
c473b43d | 111 | if(TMath::Abs(meanZNA)>1.e-10){ |
213a7d34 | 112 | if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) |
113 | res=1.; | |
114 | else percentageDiff= | |
115 | res=.5; | |
116 | test[specie] += res; | |
a42ceb0e | 117 | count++; |
213a7d34 | 118 | } |
c473b43d | 119 | if(TMath::Abs(meanZPC)>1.e-10){ |
213a7d34 | 120 | if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) |
121 | res=1.; | |
122 | else | |
123 | res=.5; | |
124 | test[specie] += res; | |
a42ceb0e | 125 | count++; |
213a7d34 | 126 | } |
c473b43d | 127 | if(TMath::Abs(meanZPA)>1.e-10){ |
213a7d34 | 128 | if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) |
129 | res=1.; | |
130 | else | |
131 | res=.5; | |
132 | test[specie] += res; | |
a42ceb0e | 133 | count++; |
213a7d34 | 134 | } |
135 | } | |
136 | ihitHisto++; | |
fdc38bb2 | 137 | } |
213a7d34 | 138 | // Check DIGIT HIGH GAIN CHAIN histos |
139 | else if(!(strncmp(hdata->GetName(),"hDig",4))){ | |
140 | if(hdata->GetEntries()>0){ | |
141 | if(idigHisto==0) sumADCZNC = hdata->GetMean(); | |
142 | else if(idigHisto==1) sumADCZNA = hdata->GetMean(); | |
143 | else if(idigHisto==2) sumADCZPC = hdata->GetMean(); | |
144 | else if(idigHisto==3) sumADCZPA = hdata->GetMean(); | |
145 | else if(idigHisto==4) pmQZNC = hdata->GetMean(); | |
146 | else if(idigHisto==5) pmQZNA = hdata->GetMean(); | |
147 | else if(idigHisto==6) pmQZPC = hdata->GetMean(); | |
148 | else if(idigHisto==7) pmQZPA = hdata->GetMean(); | |
149 | else if(idigHisto==8) pmCZNC = hdata->GetMean(); | |
150 | else if(idigHisto==9) pmCZNA = hdata->GetMean(); | |
151 | else if(idigHisto==10) pmCZPC = hdata->GetMean(); | |
152 | else if(idigHisto==11) pmCZPA = hdata->GetMean(); | |
153 | } | |
154 | // | |
155 | // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff | |
156 | if(idigHisto==11){ | |
c473b43d | 157 | if(TMath::Abs(sumADCZNC)>1.e-10){ |
213a7d34 | 158 | if((TMath::Abs(adcQZNC-adcCZNC)/adcCZNC)<percentageDiff) |
159 | res=1.; | |
160 | else | |
161 | res=.5; | |
162 | test[specie] += res; | |
a42ceb0e | 163 | count++; |
213a7d34 | 164 | } |
c473b43d | 165 | if(TMath::Abs(sumADCZNA)>1.e-10){ |
213a7d34 | 166 | if((TMath::Abs(adcQZNA-adcCZNA)/adcCZNA)<percentageDiff) |
167 | res=1.; | |
168 | else | |
169 | res=.5; | |
170 | test[specie] += res; | |
a42ceb0e | 171 | count++; |
213a7d34 | 172 | } |
c473b43d | 173 | if(TMath::Abs(sumADCZPC)>1.e-10){ |
213a7d34 | 174 | if((TMath::Abs(adcQZPC-adcCZPC)/adcCZPC)<percentageDiff) |
175 | res=1.; | |
176 | else | |
177 | res=.5; | |
178 | test[specie] += res; | |
a42ceb0e | 179 | count++; |
213a7d34 | 180 | } |
c473b43d | 181 | if(TMath::Abs(sumADCZPA)>1.e-10){ |
213a7d34 | 182 | if((TMath::Abs(adcQZPA-adcCZPA)/adcCZPA)<percentageDiff) |
183 | res=1.; | |
184 | else | |
185 | res=.5; | |
186 | test[specie] += res; | |
a42ceb0e | 187 | count++; |
213a7d34 | 188 | } |
189 | } | |
190 | idigHisto++; | |
fdc38bb2 | 191 | } |
213a7d34 | 192 | } |
193 | // ------------------------------------------------------------------- | |
194 | else if(index == AliQAv1::kRAW) { | |
213a7d34 | 195 | // |
196 | // Check RAW HIGH GAIN CHAIN histos | |
891ddd89 | 197 | |
198 | messages.Clear(); | |
199 | ||
200 | Bool_t iDetPM = kTRUE; | |
201 | // --- Checks | |
4998d47e | 202 | if(irawHisto==16) nentries = Int_t (hdata->GetEntries()); |
92664bc8 | 203 | if(irawHisto==18){ |
891ddd89 | 204 | Float_t resADC=0.; |
205 | for(int ibin=1; ibin<=hdata->GetNbinsX(); ibin++){ | |
206 | if((hdata->GetBinContent(ibin))>10.){ | |
207 | res=1.; | |
208 | } | |
17ffe44a | 209 | else if((hdata->GetBinContent(ibin))<10.){ |
891ddd89 | 210 | res=0.5; |
17ffe44a | 211 | if(ibin==1 || ibin==6 || ibin==11 || ibin==12 || ibin==13 || ibin==18){ |
212 | iDetPM = kFALSE; | |
213 | } | |
891ddd89 | 214 | } |
215 | // | |
216 | resADC += res; | |
217 | test[specie] += res; | |
218 | count++; | |
219 | } | |
4998d47e | 220 | if(nentries != -99) messages.Add(new TObjString(Form("#entries %d",nentries))); |
221 | else messages.Add(new TObjString("#entries not known")); | |
222 | // | |
891ddd89 | 223 | Float_t rv=1.; |
224 | if(hdata->GetNbinsX() != 0) rv = resADC/hdata->GetNbinsX(); | |
225 | if(rv == 1.) messages.Add(new TObjString("ADCs are OK!")); | |
226 | else if(iDetPM==kFALSE){ | |
4998d47e | 227 | messages.Add(new TObjString("Problem with some ADC!")); |
228 | messages.Add(new TObjString("IF THIS IS A PHYSICS RUN")); | |
891ddd89 | 229 | } |
230 | else if(iDetPM==kTRUE) messages.Add(new TObjString("Minor problem with ADCs")); | |
231 | SetupHisto(messages, *hdata, rv); | |
213a7d34 | 232 | } |
92664bc8 | 233 | else if(irawHisto==19){ |
ce9bf216 | 234 | // Reference values must be inserted in the order: |
235 | // ZNC, ZPC, ZNA, ZPA, ZEM1, ZEM2 | |
f33cb8ca | 236 | // Reference values from RUN 145456 |
ce9bf216 | 237 | //Double_t refTDCs[6] = {-322.5,-319.1,-320.9,-319.2,-319.7,-319.2}; |
238 | // 2012 -> Reference values from RUN 177399 | |
239 | Double_t refTDCs[6] = {-322.7,-321.4,-321.6,-321.7,-316.2,-315.4}; | |
891ddd89 | 240 | Float_t resTDC=0.; |
414b1eb6 | 241 | for(int ibin=1; ibin<=hdata->GetNbinsX(); ibin++){ |
f33cb8ca | 242 | if(TMath::Abs((hdata->GetBinContent(ibin))-refTDCs[ibin-1])<3.){ |
891ddd89 | 243 | res=1.; |
244 | } | |
f33cb8ca | 245 | else if(TMath::Abs((hdata->GetBinContent(ibin))-refTDCs[ibin-1])<5.){ |
246 | res=0.8; | |
891ddd89 | 247 | } |
f33cb8ca | 248 | else res=0.5; |
891ddd89 | 249 | // |
250 | resTDC += res; | |
251 | test[specie] += res; | |
252 | count++; | |
253 | } | |
254 | Float_t rv=1.; | |
414b1eb6 | 255 | if(hdata->GetNbinsX() != 0) rv = resTDC/hdata->GetNbinsX(); |
891ddd89 | 256 | if(rv == 1.) messages.Add(new TObjString("TDCs are OK!")); |
f33cb8ca | 257 | else if(rv<1. && rv>0.75) messages.Add(new TObjString("Minor problem with TDCs")); |
891ddd89 | 258 | else{ |
259 | messages.Add(new TObjString("Serious problem in ZDC timing")); | |
4998d47e | 260 | messages.Add(new TObjString("IF THIS IS A PHYSICS RUN")); |
891ddd89 | 261 | } |
262 | SetupHisto(messages, *hdata, rv); | |
263 | } | |
891ddd89 | 264 | irawHisto++; |
265 | ||
213a7d34 | 266 | } |
267 | // ------------------------------------------------------------------- | |
f5e4ee59 | 268 | else if(index == AliQAv1::kREC) { |
f5e4ee59 | 269 | // |
270 | // Check REC HIGH GAIN CHAIN histos | |
271 | if(hdata->GetEntries()>0){ | |
c473b43d | 272 | if(irecHisto==0) meanZNC = hdata->GetMean(); |
273 | else if(irecHisto==1) meanZNA = hdata->GetMean(); | |
274 | else if(irecHisto==2) meanZPC = hdata->GetMean(); | |
275 | else if(irecHisto==3) meanZPA = hdata->GetMean(); | |
f5e4ee59 | 276 | else if(irecHisto==4) pmQZNC = hdata->GetMean(); |
277 | else if(irecHisto==5) pmQZNA = hdata->GetMean(); | |
278 | else if(irecHisto==6) pmQZPC = hdata->GetMean(); | |
279 | else if(irecHisto==7) pmQZPA = hdata->GetMean(); | |
280 | else if(irecHisto==8) pmCZNC = hdata->GetMean(); | |
281 | else if(irecHisto==9) pmCZNA = hdata->GetMean(); | |
282 | else if(irecHisto==10) pmCZPC = hdata->GetMean(); | |
283 | else if(irecHisto==11) pmCZPA = hdata->GetMean(); | |
284 | } | |
285 | // | |
286 | // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff | |
287 | if(irecHisto==11){ | |
c473b43d | 288 | if(TMath::Abs(meanZNC)>1.e-10){ |
f5e4ee59 | 289 | if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) |
290 | res=1.; | |
291 | else | |
292 | res=.5; | |
293 | test[specie] += res; | |
a42ceb0e | 294 | count++; |
f5e4ee59 | 295 | } |
c473b43d | 296 | if(TMath::Abs(meanZNA)>1.e-10){ |
f5e4ee59 | 297 | if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) |
298 | res=1.; | |
299 | else | |
300 | res=.5; | |
301 | test[specie] += res; | |
a42ceb0e | 302 | count++; |
f5e4ee59 | 303 | } |
c473b43d | 304 | if(TMath::Abs(meanZPC)>1.e-10){ |
f5e4ee59 | 305 | if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) |
306 | res=1.; | |
307 | else | |
308 | res=.5; | |
309 | test[specie] += res; | |
a42ceb0e | 310 | count++; |
f5e4ee59 | 311 | } |
c473b43d | 312 | if(TMath::Abs(meanZPA)>1.e-10){ |
f5e4ee59 | 313 | if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) |
314 | res=1.; | |
315 | else | |
316 | res=.5; | |
317 | test[specie] += res; | |
a42ceb0e | 318 | count++; |
f5e4ee59 | 319 | } |
320 | } | |
321 | irecHisto++; | |
322 | } | |
323 | // ------------------------------------------------------------------- | |
213a7d34 | 324 | else if(index == AliQAv1::kESD) { |
213a7d34 | 325 | // |
326 | // Check ESD HIGH GAIN CHAIN histos | |
c473b43d | 327 | if(hdata->GetEntries()>0){ |
92664bc8 | 328 | if(esdInd==0) sumADCZNC = hdata->GetMean(); |
329 | else if(esdInd==1) sumADCZNA = hdata->GetMean(); | |
330 | else if(esdInd==2) sumADCZPC = hdata->GetMean(); | |
331 | else if(esdInd==3) sumADCZPA = hdata->GetMean(); | |
332 | else if(esdInd==6) pmQZNC = hdata->GetMean(); | |
333 | else if(esdInd==7) pmQZNA = hdata->GetMean(); | |
334 | else if(esdInd==8) pmQZPC = hdata->GetMean(); | |
335 | else if(esdInd==9) pmQZPA = hdata->GetMean(); | |
336 | else if(esdInd==10) pmCZNC = hdata->GetMean(); | |
337 | else if(esdInd==11) pmCZNA = hdata->GetMean(); | |
338 | else if(esdInd==12) pmCZPC = hdata->GetMean(); | |
339 | else if(esdInd==13) pmCZPA = hdata->GetMean(); | |
213a7d34 | 340 | } |
341 | // | |
342 | // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff | |
c473b43d | 343 | if(esdInd==15){ |
344 | if(TMath::Abs(sumADCZNC)>1.e-10){ | |
213a7d34 | 345 | if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) |
346 | res=1.; | |
347 | else | |
348 | res=.5; | |
349 | test[specie] += res; | |
a42ceb0e | 350 | count++; |
213a7d34 | 351 | } |
c473b43d | 352 | if(TMath::Abs(sumADCZNA)>1.e-10){ |
213a7d34 | 353 | if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) |
354 | res=1.; | |
355 | else | |
356 | res=.5; | |
357 | test[specie] += res; | |
a42ceb0e | 358 | count++; |
213a7d34 | 359 | } |
c473b43d | 360 | if(TMath::Abs(sumADCZPC)>1.e-10){ |
213a7d34 | 361 | if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) |
362 | res=1.; | |
363 | else | |
364 | res=.5; | |
365 | test[specie] += res; | |
a42ceb0e | 366 | count++; |
213a7d34 | 367 | } |
c473b43d | 368 | if(TMath::Abs(sumADCZPA)>1.e-10){ |
213a7d34 | 369 | if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) |
370 | res=1.; | |
371 | else | |
372 | res=.5; | |
373 | test[specie] += res; | |
a42ceb0e | 374 | count++; |
213a7d34 | 375 | } |
fdc38bb2 | 376 | } |
57acd2d2 | 377 | esdInd++; |
213a7d34 | 378 | } |
379 | else { | |
57acd2d2 | 380 | AliWarning(Form("\n\t No ZDC QA for %s task\n",taskName)); |
a42ceb0e | 381 | return ; |
57acd2d2 | 382 | } |
213a7d34 | 383 | }//if(hdata) |
384 | else AliError("AliZDCQAChecker-> No histos!!!\n"); | |
075a0e70 | 385 | } |
213a7d34 | 386 | } // LowMult (p-p) |
57acd2d2 | 387 | // ==================================================================== |
388 | // Checks for A-A events | |
389 | // ==================================================================== | |
4998d47e | 390 | if (AliRecoParam::ConvertIndex(specie) == AliRecoParam::kHighMult) { |
57acd2d2 | 391 | if(list[specie]->GetEntries()==0){ |
213a7d34 | 392 | AliWarning("\t The list to be checked is empty!"); |
a42ceb0e | 393 | return ; |
57acd2d2 | 394 | } |
5379c4a3 | 395 | //AliDebug(AliQAv1::GetQADebugLevel(), Form("\n\tAliZDCQAChecker-> checking QA histograms for task %s\n\n",taskName)); |
213a7d34 | 396 | // |
57acd2d2 | 397 | TIter next(list[specie]); |
a42ceb0e | 398 | count = 0; |
57acd2d2 | 399 | TH1 * hdata; |
c473b43d | 400 | |
213a7d34 | 401 | Float_t res=0., percentageDiff=0.10; |
c473b43d | 402 | Float_t meanZNA=0., meanZNC=0., meanZPA=0., meanZPC=0.; |
403 | Float_t pmCZNA=0., pmCZNC=0., pmCZPA=0., pmCZPC=0.; | |
404 | Float_t pmQZNA=0., pmQZNC=0., pmQZPA=0., pmQZPC=0.; | |
405 | Float_t sumADCZNA=0., sumADCZNC=0., sumADCZPA=0., sumADCZPC=0.; | |
406 | Float_t adcCZNA=0., adcCZNC=0., adcCZPA=0., adcCZPC=0.; | |
407 | Float_t adcQZNA=0., adcQZNC=0., adcQZPA=0., adcQZPC=0.; | |
4998d47e | 408 | Int_t nentries=-99; |
c473b43d | 409 | |
57acd2d2 | 410 | while((hdata = dynamic_cast<TH1 *>(next()))){ |
411 | if(hdata){ | |
5379c4a3 | 412 | //AliDebug(AliQAv1::GetQADebugLevel(), Form("\tAliZDCQAChecker-> checking histo %s",hdata->GetName())); |
57acd2d2 | 413 | // ------------------------------------------------------------------- |
4e25ac79 | 414 | if(index == AliQAv1::kSIM){ |
213a7d34 | 415 | // Check HITS histos |
416 | if (!(strncmp(hdata->GetName(),"hHits",5))){ | |
417 | if(hdata->GetEntries()>0){ | |
c473b43d | 418 | if(ihitHisto==0) meanZNC = hdata->GetMean(); |
419 | else if(ihitHisto==1) meanZNA = hdata->GetMean(); | |
420 | else if(ihitHisto==2) meanZPC = hdata->GetMean(); | |
421 | else if(ihitHisto==3) meanZPA = hdata->GetMean(); | |
422 | else if(ihitHisto==4) pmQZNC = hdata->GetMean(); | |
423 | else if(ihitHisto==5) pmQZNA = hdata->GetMean(); | |
424 | else if(ihitHisto==6) pmQZPC = hdata->GetMean(); | |
425 | else if(ihitHisto==7) pmQZPA = hdata->GetMean(); | |
426 | else if(ihitHisto==8) pmCZNC = hdata->GetMean(); | |
427 | else if(ihitHisto==9) pmCZNA = hdata->GetMean(); | |
428 | else if(ihitHisto==10) pmCZPC = hdata->GetMean(); | |
429 | else if(ihitHisto==11) pmCZPA = hdata->GetMean(); | |
213a7d34 | 430 | } |
431 | // | |
432 | // --- Check whether 2*|Mean ZNA - Mean ZNC|/(Mean ZNA + Mean ZNC) < percentageDiff | |
433 | // --- and 2*|Mean ZPA - Mean ZPC|/(Mean ZPA + Mean ZPC) < 2*percentageDiff | |
c473b43d | 434 | if(ihitHisto==3){ |
435 | if(TMath::Abs(meanZNC)>1.e-10 && TMath::Abs(meanZNA)>1.e-10){ | |
213a7d34 | 436 | if((2*TMath::Abs(meanZNC-meanZNA)/(meanZNA+meanZNC))<percentageDiff) |
437 | res=1.; | |
438 | else | |
439 | res=.5; | |
440 | test[specie] += res; | |
a42ceb0e | 441 | count++; |
213a7d34 | 442 | } |
c473b43d | 443 | if(TMath::Abs(meanZPC)>1.e-10 && TMath::Abs(meanZPA)>1.e-10){ |
213a7d34 | 444 | if((TMath::Abs(meanZPC-meanZPA)/(meanZPA+meanZPC))<percentageDiff) |
445 | res=1.; | |
446 | else | |
447 | res=.5; | |
448 | test[specie] += res; | |
a42ceb0e | 449 | count++; |
213a7d34 | 450 | } |
57acd2d2 | 451 | } |
213a7d34 | 452 | // --- Check whether (mean PMQi - PMC)/PMC < percentageDiff |
c473b43d | 453 | if(ihitHisto==11){ |
454 | if(TMath::Abs(meanZNC)>1.e-10){ | |
213a7d34 | 455 | if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) |
57acd2d2 | 456 | res=1.; |
457 | else | |
213a7d34 | 458 | res=.5; |
57acd2d2 | 459 | test[specie] += res; |
a42ceb0e | 460 | count++; |
213a7d34 | 461 | } |
c473b43d | 462 | if(TMath::Abs(meanZNA)>1.e-10){ |
213a7d34 | 463 | if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) |
57acd2d2 | 464 | res=1.; |
465 | else | |
213a7d34 | 466 | res=.5; |
57acd2d2 | 467 | test[specie] += res; |
a42ceb0e | 468 | count++; |
213a7d34 | 469 | } |
c473b43d | 470 | if(TMath::Abs(meanZPC)>1.e-10){ |
213a7d34 | 471 | if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) |
57acd2d2 | 472 | res=1.; |
473 | else | |
213a7d34 | 474 | res=.5; |
57acd2d2 | 475 | test[specie] += res; |
a42ceb0e | 476 | count++; |
213a7d34 | 477 | } |
c473b43d | 478 | if(TMath::Abs(meanZPA)>1.e-10){ |
213a7d34 | 479 | if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) |
57acd2d2 | 480 | res=1.; |
481 | else | |
213a7d34 | 482 | res=.5; |
57acd2d2 | 483 | test[specie] += res; |
a42ceb0e | 484 | count++; |
213a7d34 | 485 | } |
486 | } | |
c473b43d | 487 | ihitHisto++; |
fdc38bb2 | 488 | } |
213a7d34 | 489 | // Check DIGITS histos |
490 | else if (!(strncmp(hdata->GetName(),"hDig",4))){ | |
491 | if(hdata->GetEntries()>0){ | |
c473b43d | 492 | if(idigHisto==0) sumADCZNC = hdata->GetMean(); |
493 | else if(idigHisto==1) sumADCZNA = hdata->GetMean(); | |
494 | else if(idigHisto==2) sumADCZPC = hdata->GetMean(); | |
495 | else if(idigHisto==3) sumADCZPA = hdata->GetMean(); | |
496 | else if(idigHisto==4) adcQZNC = hdata->GetMean(); | |
497 | else if(idigHisto==5) adcQZNA = hdata->GetMean(); | |
498 | else if(idigHisto==6) adcQZPC = hdata->GetMean(); | |
499 | else if(idigHisto==7) adcQZPA = hdata->GetMean(); | |
500 | else if(idigHisto==8) adcCZNC = hdata->GetMean(); | |
501 | else if(idigHisto==9) adcCZNA = hdata->GetMean(); | |
502 | else if(idigHisto==10) adcCZPC = hdata->GetMean(); | |
503 | else if(idigHisto==11) adcCZPA = hdata->GetMean(); | |
213a7d34 | 504 | } |
505 | // | |
506 | // --- Check whether 2*|Mean ZNA - Mean ZNC|/(Mean ZNA + Mean ZNC) < percentageDiff | |
507 | // --- and 2*|Mean ZPA - Mean ZPC|/(Mean ZPA + Mean ZPC) < 2*percentageDiff | |
c473b43d | 508 | if(idigHisto==3){ |
509 | if(TMath::Abs(sumADCZNC)>1.e-10 && TMath::Abs(sumADCZNA)>1.e-10){ | |
213a7d34 | 510 | if((2*TMath::Abs(sumADCZNC-sumADCZNA)/(sumADCZNA+sumADCZNC))<percentageDiff) |
57acd2d2 | 511 | res=1.; |
512 | else | |
513 | res=.5; | |
514 | test[specie] += res; | |
a42ceb0e | 515 | count++; |
213a7d34 | 516 | } |
c473b43d | 517 | if(TMath::Abs(sumADCZPC)>1.e-10 && TMath::Abs(sumADCZPA)>1.e-10){ |
213a7d34 | 518 | if((TMath::Abs(sumADCZPC-sumADCZPA)/(sumADCZPA+sumADCZPC))<percentageDiff) |
57acd2d2 | 519 | res=1.; |
520 | else | |
521 | res=.5; | |
522 | test[specie] += res; | |
a42ceb0e | 523 | count++; |
213a7d34 | 524 | } |
57acd2d2 | 525 | } |
213a7d34 | 526 | // --- Check whether (sumADC PMQi - PMC)/PMC < percentageDiff |
c473b43d | 527 | if(idigHisto==11){ |
528 | if(TMath::Abs(sumADCZNC)>1.e-10){ | |
213a7d34 | 529 | if((TMath::Abs(adcQZNC-adcCZNC)/adcCZNC)<percentageDiff) |
57acd2d2 | 530 | res=1.; |
531 | else | |
213a7d34 | 532 | res=.5; |
57acd2d2 | 533 | test[specie] += res; |
a42ceb0e | 534 | count++; |
213a7d34 | 535 | } |
c473b43d | 536 | if(TMath::Abs(sumADCZNA)>1.e-10){ |
213a7d34 | 537 | if((TMath::Abs(adcQZNA-adcCZNA)/adcCZNA)<percentageDiff) |
57acd2d2 | 538 | res=1.; |
539 | else | |
213a7d34 | 540 | res=.5; |
57acd2d2 | 541 | test[specie] += res; |
a42ceb0e | 542 | count++; |
213a7d34 | 543 | } |
c473b43d | 544 | if(TMath::Abs(sumADCZPC)>1.e-10){ |
213a7d34 | 545 | if((TMath::Abs(adcQZPC-adcCZPC)/adcCZPC)<percentageDiff) |
57acd2d2 | 546 | res=1.; |
547 | else | |
213a7d34 | 548 | res=.5; |
57acd2d2 | 549 | test[specie] += res; |
a42ceb0e | 550 | count++; |
213a7d34 | 551 | } |
c473b43d | 552 | if(TMath::Abs(sumADCZPA)>1.e-10){ |
213a7d34 | 553 | if((TMath::Abs(adcQZPA-adcCZPA)/adcCZPA)<percentageDiff) |
57acd2d2 | 554 | res=1.; |
555 | else | |
213a7d34 | 556 | res=.5; |
57acd2d2 | 557 | test[specie] += res; |
a42ceb0e | 558 | count++; |
213a7d34 | 559 | } |
560 | } | |
c473b43d | 561 | idigHisto++; |
fdc38bb2 | 562 | } |
563 | } | |
213a7d34 | 564 | // ------------------------------------------------------------------- |
565 | else if(index == AliQAv1::kRAW){ | |
213a7d34 | 566 | // |
567 | // Check RAW HIGH GAIN CHAIN histos | |
891ddd89 | 568 | |
569 | messages.Clear(); | |
570 | ||
571 | Bool_t iDetPM = kTRUE; | |
572 | // --- Checks | |
4998d47e | 573 | if(irawHisto==16) nentries = Int_t (hdata->GetEntries()); |
92664bc8 | 574 | if(irawHisto==18){ |
891ddd89 | 575 | Float_t resADC=0.; |
576 | for(int ibin=1; ibin<=hdata->GetNbinsX(); ibin++){ | |
577 | if((hdata->GetBinContent(ibin))>10.){ | |
578 | res=1.; | |
579 | } | |
580 | else if((hdata->GetBinContent(ibin))<10. && | |
581 | (ibin==1 || ibin==6 || ibin==11 || ibin==12 || ibin==13 || ibin==18)){ | |
582 | res=0.5; | |
583 | iDetPM = kFALSE; | |
584 | } | |
585 | // | |
586 | resADC += res; | |
587 | test[specie] += res; | |
588 | count++; | |
589 | } | |
4998d47e | 590 | if(nentries != -99) messages.Add(new TObjString(Form("#entries %d",nentries))); |
591 | else messages.Add(new TObjString("#entries not known")); | |
592 | // | |
891ddd89 | 593 | Float_t rv=1.; |
594 | if(hdata->GetNbinsX() != 0) rv = resADC/hdata->GetNbinsX(); | |
4998d47e | 595 | if(rv > 0.98) messages.Add(new TObjString("ADCs are OK!")); |
891ddd89 | 596 | else if(iDetPM==kFALSE){ |
597 | messages.Add(new TObjString("Problem with ADCs!")); | |
4998d47e | 598 | messages.Add(new TObjString("IF THIS IS A PHYSICS RUN")); |
891ddd89 | 599 | } |
600 | else if(iDetPM==kTRUE) messages.Add(new TObjString("Minor problem with ADCs")); | |
601 | SetupHisto(messages, *hdata, rv); | |
213a7d34 | 602 | } |
92664bc8 | 603 | else if(irawHisto==19){ |
4998d47e | 604 | // Reference values from RUN 137161 |
605 | Double_t refTDCs[6] = {-320.7,-319.0,-318.6,-319.9,-321.3,-320.8}; | |
891ddd89 | 606 | Float_t resTDC=0.; |
607 | for(int ibin=1; ibin<=hdata->GetNbinsX(); ibin++){ | |
4998d47e | 608 | if(TMath::Abs((hdata->GetBinContent(ibin))-refTDCs[ibin-1])<3.){ |
891ddd89 | 609 | res=1.; |
610 | } | |
4998d47e | 611 | else if(TMath::Abs((hdata->GetBinContent(ibin))-refTDCs[ibin-1])<4.){ |
612 | res=0.8; | |
891ddd89 | 613 | } |
4998d47e | 614 | else res=0.5; |
891ddd89 | 615 | // |
616 | resTDC += res; | |
617 | test[specie] += res; | |
618 | count++; | |
619 | } | |
620 | Float_t rv=1.; | |
621 | if(hdata->GetNbinsX() != 0) rv = resTDC/hdata->GetNbinsX(); | |
622 | if(rv == 1.) messages.Add(new TObjString("TDCs are OK!")); | |
4998d47e | 623 | else if(rv<1 && rv>0.75) messages.Add(new TObjString("Minor problem with TDCs")); |
891ddd89 | 624 | else{ |
625 | messages.Add(new TObjString("Serious problem in ZDC timing")); | |
4998d47e | 626 | messages.Add(new TObjString("IF THIS IS A PHYSICS RUN")); |
891ddd89 | 627 | } |
628 | SetupHisto(messages, *hdata, rv); | |
629 | } | |
891ddd89 | 630 | irawHisto++; |
631 | ||
213a7d34 | 632 | } |
633 | // ------------------------------------------------------------------- | |
f5e4ee59 | 634 | else if(index == AliQAv1::kREC){ |
f5e4ee59 | 635 | // |
636 | // Check RAW HIGH GAIN CHAIN histos | |
637 | if(hdata->GetEntries()>0){ | |
c473b43d | 638 | if(irecHisto==0) meanZNC = hdata->GetMean(); |
639 | else if(irecHisto==1) meanZNA = hdata->GetMean(); | |
640 | else if(irecHisto==2) meanZPC = hdata->GetMean(); | |
641 | else if(irecHisto==3) meanZPA = hdata->GetMean(); | |
642 | else if(irecHisto==4) pmQZNC = hdata->GetMean(); | |
643 | else if(irecHisto==5) pmQZNA = hdata->GetMean(); | |
644 | else if(irecHisto==6) pmQZPC = hdata->GetMean(); | |
645 | else if(irecHisto==7) pmQZPA = hdata->GetMean(); | |
646 | else if(irecHisto==8) pmCZNC = hdata->GetMean(); | |
647 | else if(irecHisto==9) pmCZNA = hdata->GetMean(); | |
648 | else if(irecHisto==10) pmCZPC = hdata->GetMean(); | |
649 | else if(irecHisto==11) pmCZPA = hdata->GetMean(); | |
f5e4ee59 | 650 | } |
651 | // | |
652 | // --- Check whether 2*|Mean ZNA - Mean ZNC|/(Mean ZNA + Mean ZNC) < percentageDiff | |
653 | // --- and 2*|Mean ZPA - Mean ZPC|/(Mean ZPA + Mean ZPC) < 2*percentageDiff | |
c473b43d | 654 | if(irecHisto==3){ |
655 | if(TMath::Abs(meanZNC)>1.e-10 && TMath::Abs(meanZNA)>1.e-10){ | |
656 | if((2*TMath::Abs(meanZNC-meanZNA)/(meanZNA+meanZNC))<percentageDiff) | |
f5e4ee59 | 657 | res=1.; |
658 | else | |
659 | res=.5; | |
660 | test[specie] += res; | |
a42ceb0e | 661 | count++; |
f5e4ee59 | 662 | } |
c473b43d | 663 | if(TMath::Abs(meanZPC)>1.e-10 && TMath::Abs(meanZPA)>1.e-10){ |
664 | if((TMath::Abs(meanZPC-meanZPA)/(meanZPA+meanZPC))<percentageDiff) | |
f5e4ee59 | 665 | res=1.; |
666 | else | |
667 | res=.5; | |
668 | test[specie] += res; | |
a42ceb0e | 669 | count++; |
f5e4ee59 | 670 | } |
671 | } | |
672 | // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff | |
c473b43d | 673 | if(irecHisto==11){ |
9e05925b | 674 | if((TMath::Abs(meanZNC)>1.e-10) && (pmCZNC>1.e-10)){ |
f5e4ee59 | 675 | if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) |
676 | res=1.; | |
677 | else | |
678 | res=.5; | |
679 | test[specie] += res; | |
a42ceb0e | 680 | count++; |
f5e4ee59 | 681 | } |
9e05925b | 682 | if((TMath::Abs(meanZNA)>1.e-10) && (pmCZNA>1.e-10)){ |
f5e4ee59 | 683 | if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) |
684 | res=1.; | |
685 | else | |
686 | res=.5; | |
687 | test[specie] += res; | |
a42ceb0e | 688 | count++; |
f5e4ee59 | 689 | } |
9e05925b | 690 | if((TMath::Abs(meanZPC)>1.e-10) && (pmCZPC>1.e-10)){ |
f5e4ee59 | 691 | if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) |
692 | res=1.; | |
693 | else | |
694 | res=.5; | |
695 | test[specie] += res; | |
a42ceb0e | 696 | count++; |
f5e4ee59 | 697 | } |
9e05925b | 698 | if((TMath::Abs(meanZPA)>1.e-10) && (pmCZPA>1.e-10)){ |
f5e4ee59 | 699 | if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) |
700 | res=1.; | |
701 | else | |
702 | res=.5; | |
703 | test[specie] += res; | |
a42ceb0e | 704 | count++; |
f5e4ee59 | 705 | } |
706 | } | |
c473b43d | 707 | irecHisto++; |
f5e4ee59 | 708 | } |
709 | // ------------------------------------------------------------------- | |
213a7d34 | 710 | else if(index == AliQAv1::kESD){ |
213a7d34 | 711 | // |
712 | // Check ESD HIGH GAIN CHAIN histos | |
c473b43d | 713 | if(hdata->GetEntries()>0){ |
92664bc8 | 714 | if(esdInd==0) sumADCZNC = hdata->GetMean(); |
715 | else if(esdInd==1) sumADCZNA = hdata->GetMean(); | |
716 | else if(esdInd==2) sumADCZPC = hdata->GetMean(); | |
717 | else if(esdInd==3) sumADCZPA = hdata->GetMean(); | |
718 | else if(esdInd==6) pmQZNC = hdata->GetMean(); | |
719 | else if(esdInd==7) pmQZNA = hdata->GetMean(); | |
720 | else if(esdInd==8) pmQZPC = hdata->GetMean(); | |
721 | else if(esdInd==9) pmQZPA = hdata->GetMean(); | |
722 | else if(esdInd==10) pmCZNC = hdata->GetMean(); | |
723 | else if(esdInd==11) pmCZNA = hdata->GetMean(); | |
724 | else if(esdInd==12) pmCZPC = hdata->GetMean(); | |
725 | else if(esdInd==13) pmCZPA = hdata->GetMean(); | |
213a7d34 | 726 | } |
727 | // | |
728 | // --- Check whether 2*|Mean ZNA - Mean ZNC|/(Mean ZNA + Mean ZNC) < percentageDiff | |
729 | // --- and 2*|Mean ZPA - Mean ZPC|/(Mean ZPA + Mean ZPC) < 2*percentageDiff | |
730 | if(esdInd==5){ | |
c473b43d | 731 | if(TMath::Abs(sumADCZNC)>1.e-10 && TMath::Abs(sumADCZNA)>1.e-10){ |
213a7d34 | 732 | if((2*TMath::Abs(sumADCZNC-sumADCZNA)/(sumADCZNA+sumADCZNC))<percentageDiff) |
733 | res=1.; | |
734 | else | |
735 | res=.5; | |
736 | test[specie] += res; | |
a42ceb0e | 737 | count++; |
213a7d34 | 738 | } |
c473b43d | 739 | if(TMath::Abs(sumADCZPC)>1.e-10 && TMath::Abs(sumADCZPA)>1.e-10){ |
213a7d34 | 740 | if((TMath::Abs(sumADCZPC-sumADCZPA)/(sumADCZPA+sumADCZPC))<percentageDiff) |
741 | res=1.; | |
742 | else | |
743 | res=.5; | |
744 | test[specie] += res; | |
a42ceb0e | 745 | count++; |
213a7d34 | 746 | } |
747 | } | |
748 | // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff | |
c473b43d | 749 | if(esdInd==15){ |
750 | if(TMath::Abs(sumADCZNC)>1.e-10){ | |
213a7d34 | 751 | if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) |
752 | res=1.; | |
753 | else | |
754 | res=.5; | |
755 | test[specie] += res; | |
a42ceb0e | 756 | count++; |
213a7d34 | 757 | } |
c473b43d | 758 | if(TMath::Abs(sumADCZNA)>1.e-10){ |
213a7d34 | 759 | if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) |
760 | res=1.; | |
761 | else | |
762 | res=.5; | |
763 | test[specie] += res; | |
a42ceb0e | 764 | count++; |
213a7d34 | 765 | } |
32e2fda5 | 766 | /* if(TMath::Abs(sumADCZPC)>1.e-10){ |
213a7d34 | 767 | if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) |
768 | res=1.; | |
769 | else | |
770 | res=.5; | |
771 | test[specie] += res; | |
a42ceb0e | 772 | count++; |
213a7d34 | 773 | } |
c473b43d | 774 | if(TMath::Abs(sumADCZPA)>1.e-10){ |
213a7d34 | 775 | if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) |
776 | res=1.; | |
777 | else | |
778 | res=.5; | |
779 | test[specie] += res; | |
a42ceb0e | 780 | count++; |
213a7d34 | 781 | } |
32e2fda5 | 782 | */ |
213a7d34 | 783 | } |
784 | esdInd++; | |
785 | } | |
786 | else { | |
787 | AliWarning(Form("\n\t No ZDC QA for %s task\n",taskName)); | |
a42ceb0e | 788 | return ; |
213a7d34 | 789 | } |
790 | }//if(hdata) | |
791 | else AliError("\t No histos found for ZDC!!!\n"); | |
792 | } | |
793 | } // HighMult (Pb-Pb) | |
f5e4ee59 | 794 | // ==================================================================== |
795 | // Checks for Calibration events | |
796 | // ==================================================================== | |
4998d47e | 797 | if (AliRecoParam::ConvertIndex(specie) == AliRecoParam::kCalib) { |
798 | AliWarning(Form("\n\t No check implemented in ZDC QA for %s task in CALIBRATION events\n",taskName)); | |
a42ceb0e | 799 | return ; |
f5e4ee59 | 800 | } // Calibration |
801 | // ==================================================================== | |
802 | // Checks for cosmic events | |
803 | // ==================================================================== | |
c473b43d | 804 | else if (AliRecoParam::ConvertIndex(specie) == AliRecoParam::kCosmic) { |
4998d47e | 805 | AliWarning(Form("\n\t No check needed in ZDC QA for %s task in COSMIC events\n",taskName)); |
a42ceb0e | 806 | return ; |
f5e4ee59 | 807 | } // Cosmic |
c473b43d | 808 | if(TMath::Abs(count)>1.e-10) test[specie] = test[specie]/count; |
213a7d34 | 809 | AliDebug(AliQAv1::GetQADebugLevel(), Form("\n\t ZDC QA check result = %1.2f\n",test[specie])); |
a42ceb0e | 810 | } // Loop on species |
075a0e70 | 811 | } |
891ddd89 | 812 | |
813 | //___________________________________________________________________ | |
814 | void AliZDCQAChecker::SetupHisto(const TObjArray& messages, TH1& histo, Float_t& code) | |
815 | { | |
816 | // | |
817 | /// Add text to histos | |
818 | // | |
819 | ||
9d397f52 | 820 | TPaveText* text = new TPaveText(0.70,0.70,0.99,0.99,"NDC"); |
891ddd89 | 821 | |
822 | TIter next(&messages); | |
823 | TObjString* str; | |
824 | ||
825 | while ( ( str = static_cast<TObjString*>(next()) ) ){ | |
826 | text->AddText(str->String()); | |
827 | } | |
828 | ||
829 | TString defaultText = ""; | |
830 | ||
831 | Int_t color = 0; | |
832 | if(code==1.){ | |
833 | color = kGreen; | |
834 | defaultText = "Everything is fine!"; | |
835 | } | |
4998d47e | 836 | else if(code<1. && code>=0.85){ |
891ddd89 | 837 | color = kYellow; |
838 | defaultText = "To be monitored in next runs"; | |
839 | } | |
4998d47e | 840 | else if(code<0.85 && code>=0.6){ |
891ddd89 | 841 | color = kOrange; |
842 | defaultText = "notify the expert DURING THE DAY!"; | |
843 | } | |
9d397f52 | 844 | else if(code<0.6){ |
891ddd89 | 845 | color = kRed; |
4998d47e | 846 | defaultText = "PLEASE CALL THE EXPERT!!!!"; |
891ddd89 | 847 | } |
848 | ||
849 | ||
850 | text->AddText(defaultText.Data()); | |
851 | text->SetFillColor(color); | |
852 | ||
853 | //histo.SetFillStyle(1001); | |
854 | //histo.SetFillColor(color); | |
855 | ||
856 | histo.SetStats(kFALSE); | |
321e378b | 857 | |
858 | TList* lst = histo.GetListOfFunctions(); | |
444b012a | 859 | if(!lst){ |
860 | printf(" No list found\n"); | |
861 | } | |
862 | else { | |
321e378b | 863 | TObject *stats = lst->FindObject("stats"); |
864 | lst->Remove(stats); | |
865 | TObject *obj; | |
866 | while ((obj = lst->First())) { | |
867 | while(lst->Remove(obj)) { } | |
868 | delete obj; | |
869 | } | |
870 | if(stats) lst->Add(stats); | |
444b012a | 871 | lst->Add(text); |
07e38cef | 872 | } |
891ddd89 | 873 | } |