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