]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ITS/AliITSQASDDChecker.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Add Reaction Plane
[u/mrichter/AliRoot.git] / ITS / AliITSQASDDChecker.cxx
CommitLineData
5b716935 1/**************************************************************************
2 * Copyright(c) 2007-2009, 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/* $Id$ */
17
18// *****************************************
19// Checks the quality assurance
20// by comparing with reference data
21// P. Cerello Apr 2008
22// INFN Torino
23
24// --- ROOT system ---
3647765c 25#include "TH1.h"
80b9610c 26#include <TH1D.h>
27#include <TH2.h>
8db5c9f9 28#include "TCanvas.h"
29#include "TPaveText.h"
30#include "TPad.h"
31//#include "TPaletteAxis.h"
5b716935 32// --- AliRoot header files ---
f09c0b42 33#include "AliITSQADataMakerRec.h"
5b716935 34#include "AliITSQASDDChecker.h"
c71529b0 35#include "AliLog.h"
3647765c 36#include "AliCDBEntry.h"
f09c0b42 37#include "AliCDBManager.h"
80b9610c 38#include "AliITSCalibrationSDD.h"
39#include "AliITSgeomTGeo.h"
8db5c9f9 40#include "AliQAManager.h"
41#include "AliQAv1.h"
42#include "AliQAChecker.h"
43#include "AliQACheckerBase.h"
f09c0b42 44
5b716935 45
3647765c 46ClassImp(AliITSQASDDChecker)
5b716935 47//__________________________________________________________________
48AliITSQASDDChecker& AliITSQASDDChecker::operator = (const AliITSQASDDChecker& qac )
49{
50 // Equal operator.
51 this->~AliITSQASDDChecker();
52 new(this) AliITSQASDDChecker(qac);
53 return *this;
54}
55
80b9610c 56AliITSQASDDChecker::~AliITSQASDDChecker()
57{
f09c0b42 58
59 //destructor
80b9610c 60 if(fStepBitSDD)
61 {
62 delete[] fStepBitSDD ;
63 fStepBitSDD = NULL;
64 }
65 if(fLowSDDValue)
66 {
67 delete[]fLowSDDValue;
68 fLowSDDValue=NULL;
69 }
70 if(fHighSDDValue)
71 {
72 delete[]fHighSDDValue;
73 fHighSDDValue=NULL;
74 }
75 if(fCalibration)
76 {
77 delete fCalibration;
78 fCalibration=NULL;
79 }
8db5c9f9 80 if(fImage)
81 {
82 delete []fImage;
83 fImage=NULL;
84 }
80b9610c 85} // dtor
86
5b716935 87//__________________________________________________________________
f09c0b42 88Double_t AliITSQASDDChecker::Check(AliQAv1::ALITASK_t index, const TObjArray * list, const AliDetectorRecoParam * /*recoparam*/)
89{
90 //check histograms of the different lists
80b9610c 91 AliInfo(Form("AliITSQASDDChecker called with offset: %d\n", fSubDetOffset) );
92
3647765c 93 AliDebug(1,Form("AliITSQASDDChecker called with offset: %d\n", fSubDetOffset));
3647765c 94
62b7a70f 95 Double_t sddQACheckerValue = 0.;
80b9610c 96 TH1 *hdata=NULL;
97 Double_t entries=0.;
98 Double_t entries2[2];
99 for(Int_t i=0;i<2;i++)entries2[i]=0.;
80118b05 100
101 if(!fCalibration){
102 AliCDBEntry *calibSDD = AliCDBManager::Instance()->Get("ITS/Calib/CalibSDD");
103 Bool_t cacheStatus = AliCDBManager::Instance()->GetCacheFlag();
104 if(!calibSDD)
105 {
106 AliError("Calibration object retrieval failed! SDD will not be processed");
107 fCalibration = NULL;
62b7a70f 108 sddQACheckerValue= fHighSDDValue[AliQAv1::kWARNING];
80118b05 109 }
d2f99642 110 else{
111 fCalibration = (TObjArray *)calibSDD->GetObject();
112
113 if(!cacheStatus)calibSDD->SetObject(NULL);
114 calibSDD->SetOwner(kTRUE);
115 if(!cacheStatus)
116 {
117 delete calibSDD;
118 }
119 }//end calibsdd
120 }//end f calibration
80118b05 121
78fb0b3d 122 AliInfo("Calib SDD Created\n ");
80118b05 123
78fb0b3d 124 TIter next(list);
8a55554e 125 TString results1;
126 TString results2;
127 Int_t color=1;
78fb0b3d 128
6bf3d857 129 switch(index) {
130 case AliQAv1::kRAW:{
80b9610c 131 AliInfo(Form("Check on %s\n",AliQAv1::GetAliTaskName(index)));
6bf3d857 132// if(fRawModulePattern) { delete fRawModulePattern; fRawModulePattern = 0; }
133// if(fRawL3Pattern) { delete fRawL3Pattern; fRawL3Pattern = 0; }
134// if(fRawL4Pattern) { delete fRawL4Pattern; fRawL4Pattern = 0; }
62b7a70f 135 if (list->GetEntries() == 0){AliError("Raw List for SDD is empty \n");sddQACheckerValue += fHighSDDValue[AliQAv1::kFATAL]; break;}
6bf3d857 136 TH1 *hmodule=NULL;
137 TH2 *hlayer[2];
d2f99642 138 hdata=NULL;
8a55554e 139 for(Int_t i=0;i<2;i++)hlayer[i]=NULL;
140
141 //check counters
142 Int_t emptymodules[2], filledmodules[2],emptyladders[2],filledladders[2],emptydriftregion[2], filleddriftregion[2], excludedmoduleperlayer[2], excludeddrperlayer[2], activemoduleperlayer[2],activedrperlayer[2],exactivemoduleperlayer[2],exactivedrperlayer[2];
143 Int_t excluded=0; //excluded modules
144 Int_t excludeddriftregion=0; //excluded single drift region
145 Int_t active=0; //active modules
146 Int_t activedriftregion=0; //active single drift region
147 Int_t exactive=0; //excluded modules but taking data
148 Int_t exactivedriftregion=0; //excluded single drift region but taking data
149 Int_t empty=0;
150 Int_t filled=0;
151 Int_t emptydr=0;
152 Int_t filleddr=0;
153 //Int_t emptyactivemodule=0;
154 Int_t emptyactivemoduleperlayer[2];
155 //Int_t emptydractivemodule=0;
156 Int_t emptyactivedrperlayer[2];
157 Int_t emptysum=0;
158 Int_t emptydiff=0;
159 Int_t emptydrsum=0;
160 Int_t emptydrdiff=0;
161
162 for(Int_t i=0;i<2;i++)
163 {
164 emptymodules[i]=0;
165 filledmodules[i]=0;
166 emptyladders[i]=0;
167 filledladders[i]=0;
168 emptydriftregion[i]=0;
169 filleddriftregion[i]=0;
170 excludedmoduleperlayer[i]=0;
171 excludeddrperlayer[i]=0;
172 activemoduleperlayer[i]=0;
173 activedrperlayer[i]=0;
174 exactivemoduleperlayer[i]=0;
175 exactivedrperlayer[i]=0;
176 emptyactivemoduleperlayer[i]=0;
177 emptyactivedrperlayer[i]=0;
178 }
179
180 Int_t neventsraw=0;
181
182 //take the number of events
183
184 while( (hdata = dynamic_cast<TH1* >(next())) ){
185 if (hdata){
186 TString hname=hdata->GetName();
187
188 if(hname.Contains("SDDRawDataCheck"))
189 {
190 neventsraw=(Int_t)hdata->GetBinContent(1);
191 //break;
192 }
193 else continue;
194 }//end if hdata
195 }//end while
196
197 next.Begin();
78fb0b3d 198 while( (hdata = dynamic_cast<TH1* >(next())) ){
6bf3d857 199 if (hdata){TString hname=hdata->GetName();
78fb0b3d 200 if(hname.Contains("SDDchargeMap"))continue;
8a55554e 201 if(hname.Contains("SDDDDLPattern"))continue;
202 if(hname.Contains("SDDEventSize"))continue;
203 if(hname.Contains("_RelativeOccupancy"))continue;
78fb0b3d 204 if(hname.Contains("SDDModPattern")){
6bf3d857 205 if(hname.Contains("NORM")) continue;
206 hmodule=(TH1*)hdata->Clone();
207 entries= hdata->GetEntries();
62b7a70f 208 if(AliITSQADataMakerRec::AreEqual(entries,0.)){AliWarning(Form("===================>>>>>> No entries in %s \n",hname.Data()));sddQACheckerValue += fStepBitSDD[AliQAv1::kFATAL];}//endif entries
8a55554e 209 else{
210 int modmax=hdata->GetNbinsX();
6bf3d857 211 Double_t content=0;
8a55554e 212 Int_t llay=0;
213 for(Int_t i=1;i<=modmax;i++)
214 {
215 if(i<85)llay=0;
216 else llay=1;
217 content=hdata->GetBinContent(i);
218 if(AliITSQADataMakerRec::AreEqual(content,0.))
219 {
220 empty++;
221 emptymodules[llay]++;
222 }
223 else
224 {
225 filled++;
226 filledmodules[llay]++;
227 }
228 }//end for
229 //output of the check at the level of the modules. Drift region in the following checks
230
231 AliInfo(Form(" %s : empty modules %i \t filled modules %i",hname.Data(), empty, filled));
232 AliInfo(Form(" %s : Layer 3 empty modules %i \t filled modules %i",hname.Data(), emptymodules[0], filledmodules[0]));
233 AliInfo(Form(" %s : Layer 4 empty modules %i \t filled modules %i",hname.Data(), emptymodules[1], filledmodules[1]));
234 }//end else pattern entries !=0
235 } //modpattern (1d histogram)
78fb0b3d 236 if(hname.Contains("_RelativeOccupancy")) {
6bf3d857 237 //fRawModulePattern = (TH1F *) hdata;
238 Float_t threshold = hdata->GetMean() + 4*hdata->GetRMS();
239 if(hname.Contains("L3")) AliInfo(Form("SDD check number 1: L3 mean: %f, rms: ,%f",hdata->GetMean(),hdata->GetRMS()));
240 if(hname.Contains("L4")) AliInfo(Form("SDD check number 2: L4 mean: %f, rms: ,%f",hdata->GetMean(),hdata->GetRMS()));
78fb0b3d 241 Int_t aboveThreshold = 0;
6bf3d857 242 for(Int_t k=0; k<= hdata->GetNbinsX(); k++) {if(hdata->GetBinLowEdge(k) > threshold) aboveThreshold += (int)(hdata->GetBinContent(k));}
6753619d 243 Float_t fractionAboveThreshold=0.;
244 if(hdata->GetEntries()>0.)fractionAboveThreshold=((Float_t) aboveThreshold)/hdata->GetEntries();
78fb0b3d 245 if(hname.Contains("L3")) AliInfo(Form("SDD check number 1, L3: Raw fractionAboveThreshold: %f",fractionAboveThreshold));
246 if(hname.Contains("L4")) AliInfo(Form("SDD check number 2, L4: Raw fractionAboveThreshold: %f",fractionAboveThreshold));
62b7a70f 247 if(fractionAboveThreshold > fThresholdForRelativeOccupancy) {sddQACheckerValue=fHighSDDValue[AliQAv1::kWARNING];
78fb0b3d 248 if(hname.Contains("L3")) AliInfo(Form("SDD check number 1: Set Warning (L3 Raw)"));
8a55554e 249 if(hname.Contains("L4")) AliInfo(Form("SDD check number 2: Set Warning (L4 Raw)")); } }//relativeoccupancy
250
6bf3d857 251 if(hname.Contains("SDDphizL3") || hname.Contains("SDDphizL4")){if(hname.Contains("NORM"))continue;
6bf3d857 252 Int_t layer=0;
253 if(hname.Contains("3"))layer=0;
254 else if(hname.Contains("4"))layer=1;
255 entries2[layer]=hdata->GetEntries();
256 if(entries2[layer]==0){AliWarning(Form("===================>>>>>> No entries in %s \n",hname.Data()));
62b7a70f 257 sddQACheckerValue += fStepBitSDD[AliQAv1::kFATAL];}//end if getentries
6bf3d857 258 else{
259 Int_t layer1=0;
260 if(hname.Contains("3"))layer1=0;
261 else if(hname.Contains("4"))layer1=1;
262 TH2* htemp=dynamic_cast<TH2*>(hdata);
d2f99642 263 if(htemp){
264 hlayer[layer1]=(TH2*)htemp->Clone();
265 hlayer[layer1]->SetName(Form("%s_copy",hname.Data()));
8a55554e 266 //hlayer[layer1]->RebinX(2);
d2f99642 267 int modmay=hlayer[layer1]->GetNbinsY();
268 TH1D* hproj= hlayer[layer1]->ProjectionY();
269 Double_t ladcontent=0;
270 for(Int_t i=1;i<=modmay;i++) {//loop on the ladders
271 ladcontent=hproj->GetBinContent(i);
272 if(AliITSQADataMakerRec::AreEqual(ladcontent,0.)) emptyladders[layer1]++;
273 else filledladders[layer1]++;}//end for
8a55554e 274 AliInfo(Form(" %s : empty ladders %i \t filled ladders %i\n",hname.Data(), emptyladders[layer1], filledladders[layer1]));//end else layer 3
d2f99642 275 delete hproj;
276 hproj=NULL;
277 }//end if htemp
278 }//end else entries !=0
6bf3d857 279 }//end check on phiz
78fb0b3d 280 }//end if hdata
281 }//end while
8a55554e 282 for(Int_t ii=0;ii<2;ii++)
283 {
284 filledmodules[ii]=0;
285 emptymodules[ii]=0;
286 }
287 filled=0;
288 empty=0;
289 if(AliITSQADataMakerRec::AreEqual(entries,0.)&& AliITSQADataMakerRec::AreEqual(entries2[0],0.)&& AliITSQADataMakerRec::AreEqual(entries2[1],0.)) break;
78fb0b3d 290 //else{
291 if(hmodule || (hlayer[0] && hlayer[1])){
78fb0b3d 292 for(Int_t imod=0;imod<fgknSDDmodules;imod++){
293 Int_t lay=0;
294 Int_t lad=0;
295 Int_t det=0;
296 Int_t module=0;
297 module=imod+fgkmodoffset;
298 AliITSCalibrationSDD * cal=(AliITSCalibrationSDD*)fCalibration->At(imod);
d2f99642 299 if(cal==0) { continue;}
78fb0b3d 300 AliITSgeomTGeo::GetModuleId(module,lay,lad,det);
301 if (cal->IsBad()){
302 excluded++;
8a55554e 303 excludedmoduleperlayer[lay-3]++;
78fb0b3d 304 Double_t content=0.;
305 Double_t contentlayer[2];
306 for(Int_t i=0;i<2;i++)contentlayer[i]=0.;
8a55554e 307 if(hmodule)content=hmodule->GetBinContent(imod+1);//if expert bit is active the histogram will be stored in the QA file otherwise the histogram will not be written on the logbook
308 // if(hlayer[lay-3]) contentlayer[lay-3]=hlayer[lay-3]->GetBinContent(det,lad);
309 if(AliITSQADataMakerRec::AreEqual(content,0.)== kFALSE) {
78fb0b3d 310 filledmodules[lay-3]++;
8a55554e 311 filled++;
312 AliError(Form("The module %d (layer %i, ladder %i det %i ) excluded from the acquisition, took data \n ",module,lay,lad,det));
78fb0b3d 313 exactive++;
8a55554e 314 exactivemoduleperlayer[lay-3]++;
315 } else if(AliITSQADataMakerRec::AreEqual(content,0.))
316 {
317 emptymodules[lay-3]++; //it has to be empty
318 empty++;
319 }
78fb0b3d 320 } else {
8a55554e 321 Int_t totside=0;
322 Int_t totactiveside=0;
323 //Int_t totbadside=0;
78fb0b3d 324 Double_t contentgood=0.;
8a55554e 325
326 for(Int_t i=0;i<2;i++){
327 if(hlayer[lay-3]) contentgood=hlayer[lay-3]->GetBinContent(2*det+i-1,lad);
328 if(cal->IsWingBad(i))
329 {
330 excludeddriftregion++;
331 excludeddrperlayer[lay-3]++;
332 if(AliITSQADataMakerRec::AreEqual(contentgood,0.)==kFALSE){
333 AliError(Form("The side %d of the module %d (layer %i, ladder %i det %i ) excluded from the acquisition, took data \n ",i,module,lay,lad,det));
334 exactivedriftregion++;
335 exactivedrperlayer[lay-3]++;
336 filleddr++;
337 filleddriftregion[lay-3]++;
338 }
339 }//end wingbad
340 else{
341 if(AliITSQADataMakerRec::AreEqual(contentgood,0.)==kTRUE)
342 {
343 AliWarning(Form("The side %d of the module %d (layer %i, ladder %i det %i ) is in acquisition, but it didn't take data\n ",i,module, lay, lad, det));
344 }
345 else
346 {
347 totside++;
348 }
349 totactiveside++;
350 }
351 }//end for
352 if(totside==0){
353 AliWarning(Form("The module %d (layer %i, ladder %i det %i ) is in acquisition, but it didn't take data\n ",module, lay, lad, det));
78fb0b3d 354 emptymodules[lay-3]++;
8a55554e 355 empty++;
356
357 }
358 else
359 if(totside==2){
360 filledmodules[lay-3]++;
361 filled++;
362 }
363 else
364 if(totside==1)
365 {
366 // emptydr++;
367 //emptydriftregion[lay-3]++; //it has to be empty
368 emptydriftregion[lay-3]++;
369 emptydr++;
370 filleddr++;
371 filleddriftregion[lay-3]++;
372 }
373 if(totactiveside==1)
374 {
375 activedriftregion++;
376 activedrperlayer[lay-3]++;
377 }else if(totactiveside==2)
378 {
379 active++;
380 activemoduleperlayer[lay-3]++;
381 }
382
78fb0b3d 383 }
384 }//end for
8a55554e 385 AliInfo(Form("In total %d modules and %d single drift regions took data.\n ",filled, filleddr));
386 AliInfo(Form("In total %d modules and %d single drift regions were empty\n",empty, emptydr));
387 for(Int_t i=0;i<2;i++)
388 {
389 AliInfo(Form("Layer %i \tempty modules %i \t filled modules %i\n", i+3,emptymodules[i], filledmodules[i]));
390 AliInfo(Form("Layer %i \tempty single drift regions %i \t filled single drift regions %i\n",i+3,emptydriftregion[i], filleddriftregion[i]));
391 }//end else layers
392 emptysum=emptymodules[0]+emptymodules[1];
393 emptydiff=emptysum-excluded;
394 emptyactivemoduleperlayer[0]=emptymodules[0]- excludedmoduleperlayer[0];
395 emptyactivemoduleperlayer[1]=emptymodules[1]- excludedmoduleperlayer[1];
396
397 emptydrsum=emptydriftregion[0]+emptydriftregion[1];
398 emptydrdiff=emptydrsum-excludeddriftregion;
399 emptyactivedrperlayer[0]=emptydriftregion[0]- excludeddrperlayer[0];
400 emptyactivedrperlayer[1]=emptydriftregion[1]- excludeddrperlayer[1];
401
62b7a70f 402 Int_t numlimit=1000;
8a55554e 403
62b7a70f 404 if(emptysum>excluded||emptydrsum>excludeddriftregion){
405 AliWarning(Form(" %i good module(s) and %i good single drift regions didn't take data! \n",emptydiff,emptydrdiff));
406 AliWarning(Form(" Layer 3: %i good module(s) and %i good single drift regions didn't take data! \n",emptyactivemoduleperlayer[0] ,emptyactivedrperlayer[0] ));
407 AliWarning(Form(" Layer 4: %i good module(s) and %i good single drift regions didn't take data! \n",emptyactivemoduleperlayer[1] ,emptyactivedrperlayer[1] ));
408 //printf("========================= %d",AliQAv1::Instance(AliQAv1::GetDetIndex(GetName()))->IsEventSpecieSet(AliRecoParam::kCosmic));
409 if((AliQAv1::Instance(AliQAv1::GetDetIndex("ITS")))->IsEventSpecieSet(AliRecoParam::kCosmic)==kFALSE){
410
411 results1.Form("%i good module(s) and %i good drift regions didn't take data!",emptydiff,emptydrdiff);
412 if(neventsraw<numlimit)
413 {
414 results2.Form(" Events %d .Too few events.DO NOT CALL the Expert ",neventsraw);
415 color=5;
416 sddQACheckerValue=fHighSDDValue[AliQAv1::kWARNING];
417 }
418 else
419 {
420 results2.Form(" Events %d .Follow the TWiki instruction and call the Expert ",neventsraw);
421 color=2;
422 sddQACheckerValue=fHighSDDValue[AliQAv1::kERROR];
423 }
424 }
8a55554e 425 else
62b7a70f 426 if((AliQAv1::Instance(AliQAv1::GetDetIndex("ITS")))->IsEventSpecieSet(AliRecoParam::kCosmic)==kTRUE)
427 {
428 numlimit=10000;
429 if(neventsraw<numlimit)
430 {
431 AliWarning(Form("This is a cosmic run. Some module and drift region are empty but all is OK. "));
432 results1.Form("OK. Thi is a cosmic run. you need a lot of events");
433 results2.Form("%i good module(s) and %i good drift are empty! DO NOT CALL THE EXPERT",emptydiff,emptydrdiff);
434 color=5;
435 sddQACheckerValue=fHighSDDValue[AliQAv1::kWARNING];
436 }
437 else
438 {
439 results1.Form("%i good module(s) and %i good drift region(s) have not recpoints!",emptydiff,emptydrdiff);
440 results2.Form(" Events %d .Follow the TWiki instruction and call the Expert ",neventsraw);
441 color=2;
442 sddQACheckerValue=fHighSDDValue[AliQAv1::kERROR];
443 }
444 }
8a55554e 445 }
62b7a70f 446
8a55554e 447 if(exactive==0 && emptydiff==0 && exactivedriftregion==0 && emptydrdiff==0){
448 AliInfo(Form("All the active modules (%i) and single drift regions (%i) are in acquisition. The number of excluded modules are %i and the excluded single drift regions are %i\n",active,activedriftregion,excluded,excludeddriftregion));
449 results1.Form("OK.");
450 results2.Form(" All active modules and drift regions in acquisition");
451 color=3;
62b7a70f 452 sddQACheckerValue=fHighSDDValue[AliQAv1::kINFO];
3647765c 453 }
8a55554e 454 if(exactive!=0||exactivedriftregion!=0){
455 AliError(Form("%i modules and %i single drift regions excluded from the acquisition took data. Active modules%i single drift region %i \n ",exactive,exactivedriftregion,active,activedriftregion));
456 AliError(Form("Layer 3: %i modules and %i single drift regions excluded from the acquisition took data. Active modules%i single drift region %i \n ",exactivemoduleperlayer[0],exactivedrperlayer[0],activemoduleperlayer[0],activedrperlayer[0]));
457 AliError(Form("Layer 3: %i modules and %i single drift regions excluded from the acquisition took data. Active modules%i single drift region %i \n ",exactivemoduleperlayer[1],exactivedrperlayer[1],activemoduleperlayer[1],activedrperlayer[1]));
458 results1.Form("%i modules and %i drift region excluded from the acquisition took data",exactive,exactivedriftregion);
459 results2.Form("Follow the TWiki instructions and Call the SDD expert ");
460 color=2;
62b7a70f 461 sddQACheckerValue=fHighSDDValue[AliQAv1::kERROR];
78fb0b3d 462 }
8a55554e 463 if(excluded==exactive||excludeddriftregion==exactivedriftregion){
464 AliError(Form("All the modules (%d) or single drift regions (%d) excluded from the acquisition took data!\n Active modules %i \t Active drfift regions %i\n",excluded,excludeddriftregion,active,activedriftregion));
465 results1.Form("All the modules (%d) or drift regions (%d) excluded from the acquisition took data!",excluded,excludeddriftregion );
466 results2.Form("Follow the TWiki instructions and Call the SDD expert ");
467 color=6;
62b7a70f 468 sddQACheckerValue=fHighSDDValue[AliQAv1::kFATAL];
78fb0b3d 469 }
8a55554e 470 if(active==0||activedriftregion==0){
471 AliError(Form("No modules or single drift regions took data: excluded %i \t excluded active %i \n\t\t excluded single drift regions %i \t excluded active drift regions %i \n", excluded, exactive, excludeddriftregion, exactivedriftregion));
472 results1.Form("No modules or drift region took data: excluded modules %i excluded drift regions %i ", excluded, excludeddriftregion );
473 results2.Form("Follow the TWiki instructions and Call the SDD expert ");
474 color=6;
62b7a70f 475 sddQACheckerValue=fHighSDDValue[AliQAv1::kFATAL];
78fb0b3d 476 }
8a55554e 477
478 TPaveText *pave[2];
479 next.Begin();
480
481 while( (hdata=dynamic_cast<TH1* >(next())) )
482 {
483 if (hdata){
484 TString hname=hdata->GetName();
485 if(hname.Contains("SDDphizL3") || hname.Contains("SDDphizL4")){
486 if(hname.Contains("NORM"))continue;
487 //AliInfo("========================================Found histo 11\n");
488 Int_t lay=0;
489 if(hname.Contains("3"))lay=0;
490 else if(hname.Contains("4"))lay=1;
491 pave[lay]=new TPaveText(0.3,0.88,0.9,0.99,"NDC");
492 pave[lay]->AddText(results1.Data());
493 pave[lay]->AddText(results2.Data());
494 pave[lay]->SetFillColor(color);
495 pave[lay]->SetBorderSize(1);
496 pave[lay]->SetLineWidth(1);
497 hdata->GetListOfFunctions()->Add(pave[lay]);
498 }
499 else
500 if(hname.Contains("SDDRawDataCheck"))
501 {
502
503 //AliInfo("========================================Found histo\n");
504 ((TH1F*)hdata)->SetBinContent(5,active);
505 ((TH1F*)hdata)->SetBinContent(6,filled);
506 ((TH1F*)hdata)->SetBinContent(7,activedriftregion);
507 ((TH1F*)hdata)->SetBinContent(8,filleddr);
508 ((TH1F*)hdata)->SetBinContent(9,excluded);
509 ((TH1F*)hdata)->SetBinContent(10,empty);
510 ((TH1F*)hdata)->SetBinContent(11,excludeddriftregion);
511 ((TH1F*)hdata)->SetBinContent(12,emptydr);
512 ((TH1F*)hdata)->SetBinContent(13,exactive);
513 ((TH1F*)hdata)->SetBinContent(14,emptydiff);
514 ((TH1F*)hdata)->SetBinContent(15,exactivedriftregion);
515 ((TH1F*)hdata)->SetBinContent(16,emptydr);
516
517 //layer 3
518 ((TH1F*)hdata)->SetBinContent(19,activemoduleperlayer[0]);
519 ((TH1F*)hdata)->SetBinContent(20,filledmodules[0]);
520 ((TH1F*)hdata)->SetBinContent(21,activedrperlayer[0]);
521 ((TH1F*)hdata)->SetBinContent(22,filleddriftregion[0]);
522 ((TH1F*)hdata)->SetBinContent(23,excludedmoduleperlayer[0]);
523 ((TH1F*)hdata)->SetBinContent(24,emptymodules[0]);
524 ((TH1F*)hdata)->SetBinContent(25,excludeddrperlayer[0]);
525 ((TH1F*)hdata)->SetBinContent(26,emptydriftregion[0]);
526 ((TH1F*)hdata)->SetBinContent(27,exactivemoduleperlayer[0]);
527 ((TH1F*)hdata)->SetBinContent(28,emptyactivemoduleperlayer[0]);
528 ((TH1F*)hdata)->SetBinContent(29,activedrperlayer[0]);
529 ((TH1F*)hdata)->SetBinContent(30,emptyactivedrperlayer[0]);
530
531 //layer 4
532 ((TH1F*)hdata)->SetBinContent(33,activemoduleperlayer[1]);
533 ((TH1F*)hdata)->SetBinContent(34,filledmodules[1]);
534 ((TH1F*)hdata)->SetBinContent(35,activedrperlayer[1]);
535 ((TH1F*)hdata)->SetBinContent(36,filleddriftregion[1]);
536 ((TH1F*)hdata)->SetBinContent(37,excludedmoduleperlayer[1]);
537 ((TH1F*)hdata)->SetBinContent(38,emptymodules[1]);
538 ((TH1F*)hdata)->SetBinContent(39,excludeddrperlayer[1]);
539 ((TH1F*)hdata)->SetBinContent(40,emptydriftregion[1]);
540 ((TH1F*)hdata)->SetBinContent(41,exactivemoduleperlayer[1]);
541 ((TH1F*)hdata)->SetBinContent(42,emptyactivemoduleperlayer[1]);
542 ((TH1F*)hdata)->SetBinContent(43,activedrperlayer[1]);
543 ((TH1F*)hdata)->SetBinContent(44,emptyactivedrperlayer[1]);
544 //break;
545 }
546 }//if hdata
547
548 }//end while
6bf3d857 549
78fb0b3d 550 }//end else
6bf3d857 551 delete hmodule;
552 hmodule=NULL;
553 for(Int_t i=0;i<2;i++) {
554 delete hlayer[i];
555 hlayer[i]=NULL;
8a55554e 556 }//end for
78fb0b3d 557
8a55554e 558 }//end raw
6bf3d857 559
80b9610c 560 break;
6bf3d857 561
562 case AliQAv1::kNULLTASK:{
563 AliInfo(Form("No Check on %s\n",AliQAv1::GetAliTaskName(index)));
62b7a70f 564 sddQACheckerValue=1.;
6bf3d857 565 }
566 break;
567
568 case AliQAv1::kREC:
569 {
8db5c9f9 570 Int_t uidrec=list->GetUniqueID();
80b9610c 571 AliInfo(Form("Check on %s\n",AliQAv1::GetAliTaskName(index)));
8db5c9f9 572 if(uidrec==20){
573 //recpoints
574 if (list->GetEntries() == 0){ //check if the list is empty
62b7a70f 575 //printf("sddQACheckerValue = %f \t value %f\n",sddQACheckerValue,fHighSDDValue[AliQAv1::kFATAL]);
576 sddQACheckerValue=fHighSDDValue[AliQAv1::kFATAL];
8a55554e 577 //break;
8db5c9f9 578 }//end if getentries
6bf3d857 579
8a55554e 580
581 TH1 *hmodule=NULL;
582 TH2 *hlayer[2];
583 hdata=NULL;
584 for(Int_t i=0;i<2;i++)hlayer[i]=NULL;
585
586 //check counters
587 Int_t emptymodules[2], filledmodules[2],emptyladders[2],filledladders[2],emptydriftregion[2], filleddriftregion[2], excludedmoduleperlayer[2], excludeddrperlayer[2], activemoduleperlayer[2],activedrperlayer[2],exactivemoduleperlayer[2],exactivedrperlayer[2];
588 Int_t excluded=0; //excluded modules
589 Int_t excludeddriftregion=0; //excluded single drift region
590 Int_t active=0; //active modules
591 Int_t activedriftregion=0; //active single drift region
592 Int_t exactive=0; //excluded modules but taking data
593 Int_t exactivedriftregion=0; //excluded single drift region but taking data
594 Int_t empty=0;
595 Int_t filled=0;
596 Int_t emptydr=0;
597 Int_t filleddr=0;
598 //Int_t emptyactivemodule=0;
599 Int_t emptyactivemoduleperlayer[2];
600 //Int_t emptydractivemodule=0;
601 Int_t emptyactivedrperlayer[2];
602 Int_t emptysum=0;
603 Int_t emptydiff=0;
604 Int_t emptydrsum=0;
605 Int_t emptydrdiff=0;
606
607 for(Int_t i=0;i<2;i++)
608 {
609 emptymodules[i]=0;
610 filledmodules[i]=0;
611 emptyladders[i]=0;
612 filledladders[i]=0;
613 emptydriftregion[i]=0;
614 filleddriftregion[i]=0;
615 excludedmoduleperlayer[i]=0;
616 excludeddrperlayer[i]=0;
617 activemoduleperlayer[i]=0;
618 activedrperlayer[i]=0;
619 exactivemoduleperlayer[i]=0;
620 exactivedrperlayer[i]=0;
621 emptyactivemoduleperlayer[i]=0;
622 emptyactivedrperlayer[i]=0;
623 }
624
625 Int_t neventsrecpoints=0;
626
627 while( (hdata = dynamic_cast<TH1* >(next())) ){
628 if (hdata){
629 TString hname=hdata->GetName();
630
631 if(hname.Contains("SDDRecPointCheck"))
632 {
633 neventsrecpoints=(Int_t)hdata->GetBinContent(1);
634 //break;
635 }
636 else{continue;}
637 }//end if hdata
638 }//end while
639
640 next.Begin();
641
8db5c9f9 642 while((hdata=dynamic_cast<TH1* >(next()))){
643 if (hdata){
644 TString hname=hdata->GetName();
645 if(hname.Contains("_RelativeOccupancy")) {
646 Float_t threshold = hdata->GetMean() + 4*hdata->GetRMS();
647 if(hname.Contains("L3")) AliInfo(Form("SDD check number 3: L3 mean: %f, rms: ,%f",hdata->GetMean(),hdata->GetRMS()));
648 if(hname.Contains("L4")) AliInfo(Form("SDD check number 4: L4 mean: %f, rms: ,%f",hdata->GetMean(),hdata->GetRMS()));
649 Int_t aboveThreshold = 0;
650 for(Int_t k=0; k<= ((Int_t)hdata->GetNbinsX()); k++) {
651 if(hdata->GetBinLowEdge(k) > threshold) aboveThreshold += (Int_t)(hdata->GetBinContent(k));
652 }
6753619d 653 Float_t fractionAboveThreshold=0.;
654 if(hdata->GetEntries()>0.) fractionAboveThreshold = ((Float_t) aboveThreshold)/hdata->GetEntries();
8db5c9f9 655 if(hname.Contains("L3")) AliInfo(Form("SDD check number 3, L3: RecPoints fractionAboveThreshold: %f",fractionAboveThreshold));
656 if(hname.Contains("L4")) AliInfo(Form("SDD check number 4, L4: RecPoints fractionAboveThreshold: %f",fractionAboveThreshold));
657 if(fractionAboveThreshold > fThresholdForRelativeOccupancy) {
62b7a70f 658 sddQACheckerValue=fHighSDDValue[AliQAv1::kWARNING];
8db5c9f9 659 if(hname.Contains("L3")) AliInfo(Form("SDD check number 3: Set Warning (L3 RecPoints)"));
660 if(hname.Contains("L4")) AliInfo(Form("SDD check number 4: Set Warning (L4 RecPoints)"));
661 }
78fb0b3d 662 }
8db5c9f9 663 if(hname.Contains("Rec2Raw") && !hname.Contains("2D")) {
664 //Float_t threshold = 0.;
665 if(hname.Contains("L3")) AliInfo(Form("SDD check number 5: L3 R2R mean: %f, rms: ,%f",((TH1F *) hdata)->GetMean(),((TH1F *) hdata)->GetRMS()));
666 if(hname.Contains("L4")) AliInfo(Form("SDD check number 6: L4 R2R mean: %f, rms: ,%f",((TH1F *) hdata)->GetMean(),((TH1F *) hdata)->GetRMS()));
667 Int_t belowThreshold = 0;
668 for(Int_t k=0; k<=((TH1F *)hdata)->GetNbinsX(); k++) {
669 if(((TH1F *) hdata)->GetBinLowEdge(k) < fThresholdForRecToRawRatio) belowThreshold += ((Int_t)((TH1F *) hdata)->GetBinContent(k));
670 }
671 Double_t fractionBelowThreshold =0.;
672 Double_t entries3=((TH1F *)hdata)->GetEntries();
673 if(entries3>0.001)fractionBelowThreshold = ((Double_t)(belowThreshold))/entries3;
674 else{ AliWarning(Form("No entries on %s. The check will retuns zero.\n",hdata->GetName() )); }
675 if(hname.Contains("L3")) AliInfo(Form("SDD check number 5, L3: RecPoints2Raws fractionBelowThreshold: %f",fractionBelowThreshold));
676 if(hname.Contains("L4")) AliInfo(Form("SDD check number 6, L4: RecPoints2Raws fractionBelowThreshold: %f",fractionBelowThreshold));
677 if(fractionBelowThreshold > fThresholdForRelativeOccupancy) {
62b7a70f 678 sddQACheckerValue=fHighSDDValue[AliQAv1::kWARNING];
8db5c9f9 679 if(hname.Contains("L3")) AliInfo(Form("SDD check number 5: Set Warning (L3 RecPoints2Raws)"));
680 if(hname.Contains("L4")) AliInfo(Form("SDD check number 6: Set Warning (L4 RecPoints2Raws)"));
681 }
78fb0b3d 682 }
8db5c9f9 683 if(hname.Contains("dedx")) {
684 if(hname.Contains("L3")) AliInfo(Form("SDD check number 7: L3 average charge: %f, rms: ,%f",hdata->GetMean(),hdata->GetRMS()));
685 if(hname.Contains("L4")) AliInfo(Form("SDD check number 8: L4 average charge: %f, rms: ,%f",hdata->GetMean(),hdata->GetRMS()));
8a55554e 686 }//end if dedx
687 if(hname.Contains("SDDModPatternRP")){
688 if(hname.Contains("NORM")) continue;
689 hmodule=(TH1*)hdata->Clone();
690 entries= hdata->GetEntries();
62b7a70f 691 if(AliITSQADataMakerRec::AreEqual(entries,0.)){AliWarning(Form("===================>>>>>> No entries in %s \n",hname.Data()));sddQACheckerValue += fStepBitSDD[AliQAv1::kFATAL];}//endif entries
8a55554e 692 else{
693 int modmax=hdata->GetNbinsX();
694 Double_t content=0;
695 Int_t llay=0;
696 for(Int_t i=1;i<=modmax;i++)
697 {
698 if(i<85)llay=0;
699 else llay=1;
700 content=hdata->GetBinContent(i);
701 if(AliITSQADataMakerRec::AreEqual(content,0.))
702 {
703 empty++;
704 emptymodules[llay]++;
705 }
706 else
707 {
708 filled++;
709 filledmodules[llay]++;
710 }
711 }//end for
712 //output of the check at the level of the modules. Drift region in the following checks
713
714 AliInfo(Form(" %s : empty modules %i \t filled modules %i",hname.Data(), empty, filled));
715 AliInfo(Form(" %s : Layer 3 empty modules %i \t filled modules %i",hname.Data(), emptymodules[0], filledmodules[0]));
716 AliInfo(Form(" %s : Layer 4 empty modules %i \t filled modules %i",hname.Data(), emptymodules[1], filledmodules[1]));
717 }//end else pattern entries !=0
718 } //modpattern (1d histogram)
719
720 if(hname.Contains("SDDModPatternL3RP") || hname.Contains("SDDModPatternL4RP")){if(hname.Contains("NORM"))continue;
721 Int_t layer=0;
722 if(hname.Contains("3"))layer=0;
723 else if(hname.Contains("4"))layer=1;
724 entries2[layer]=hdata->GetEntries();
725 if(entries2[layer]==0){AliWarning(Form("===================>>>>>> No entries in %s \n",hname.Data()));
62b7a70f 726 sddQACheckerValue += fStepBitSDD[AliQAv1::kFATAL];}//end if getentries
8a55554e 727 else{
728 Int_t layer1=0;
729 if(hname.Contains("3"))layer1=0;
730 else if(hname.Contains("4"))layer1=1;
731 TH2* htemp=dynamic_cast<TH2*>(hdata);
732 if(htemp){
733 hlayer[layer1]=(TH2*)htemp->Clone();
734 hlayer[layer1]->SetName(Form("%s_copy",hname.Data()));
735 //hlayer[layer1]->RebinX(2);
736 int modmay=hlayer[layer1]->GetNbinsY();
737 TH1D* hproj= hlayer[layer1]->ProjectionY();
738 Double_t ladcontent=0;
739 for(Int_t i=1;i<=modmay;i++) {//loop on the ladders
740 ladcontent=hproj->GetBinContent(i);
741 if(AliITSQADataMakerRec::AreEqual(ladcontent,0.)) emptyladders[layer1]++;
742 else filledladders[layer1]++;}//end for
743 AliInfo(Form(" %s : empty ladders %i \t filled ladders %i\n",hname.Data(), emptyladders[layer1], filledladders[layer1]));//end else layer 3
744 delete hproj;
745 hproj=NULL;
746 }//end if htemp
747 }//end else entries !=0
748 }//end check on phiz
749 }//end if hdata
750 }//end while
751 for(Int_t ii=0;ii<2;ii++)
752 {
753 filledmodules[ii]=0;
754 emptymodules[ii]=0;
755 }
756 filled=0;
757 empty=0;
758 if(AliITSQADataMakerRec::AreEqual(entries,0.)&& AliITSQADataMakerRec::AreEqual(entries2[0],0.)&& AliITSQADataMakerRec::AreEqual(entries2[1],0.)) break;
759 //else{
760 if(hmodule || (hlayer[0] && hlayer[1])){
761 for(Int_t imod=0;imod<fgknSDDmodules;imod++){
762 Int_t lay=0;
763 Int_t lad=0;
764 Int_t det=0;
765 Int_t module=0;
766 module=imod+fgkmodoffset;
767 AliITSCalibrationSDD * cal=(AliITSCalibrationSDD*)fCalibration->At(imod);
768 if(cal==0) { continue;}
769 AliITSgeomTGeo::GetModuleId(module,lay,lad,det);
770 if (cal->IsBad()){
771 excluded++;
772 excludedmoduleperlayer[lay-3]++;
773 Double_t content=0.;
774 Double_t contentlayer[2];
775 for(Int_t i=0;i<2;i++)contentlayer[i]=0.;
776 if(hmodule)content=hmodule->GetBinContent(imod+1);//if expert bit is active the histogram will be stored in the QA file otherwise the histogram will not be written on the logbook
777 // if(hlayer[lay-3]) contentlayer[lay-3]=hlayer[lay-3]->GetBinContent(det,lad);
778 if(AliITSQADataMakerRec::AreEqual(content,0.)== kFALSE) {
779 filledmodules[lay-3]++;
780 filled++;
781 AliError(Form("The module %d (layer %i, ladder %i det %i ) excluded from the acquisition,has recpoints \n ",module,lay,lad,det));
782 exactive++;
783 exactivemoduleperlayer[lay-3]++;
784 } else if(AliITSQADataMakerRec::AreEqual(content,0.))
785 {
786 emptymodules[lay-3]++; //it has to be empty
787 empty++;
788 }
789 } else {
790 Int_t totside=0;
791 Int_t totactiveside=0;
792 //Int_t totbadside=0;
793 Double_t contentgood=0.;
794
795 for(Int_t i=0;i<2;i++){
796 if(hlayer[lay-3]) contentgood=hlayer[lay-3]->GetBinContent(2*det+i-1,lad);
797 if(cal->IsWingBad(i))
798 {
799 excludeddriftregion++;
800 excludeddrperlayer[lay-3]++;
801 if(AliITSQADataMakerRec::AreEqual(contentgood,0.)==kFALSE){
802 AliError(Form("The side %d of the module %d (layer %i, ladder %i det %i ) excluded from the acquisition, has recpoints \n ",i,module,lay,lad,det));
803 exactivedriftregion++;
804 exactivedrperlayer[lay-3]++;
805 filleddr++;
806 filleddriftregion[lay-3]++;
807 }
808 }//end wingbad
809 else{
810 if(AliITSQADataMakerRec::AreEqual(contentgood,0.)==kTRUE)
811 {
812 AliWarning(Form("The side %d of the module %d (layer %i, ladder %i det %i ) is in acquisition, but no recpoints are present\n ",i,module, lay, lad, det));
813 }
814 else
815 {
816 totside++;
817 }
818 totactiveside++;
819 }
820 }//end for
821 if(totside==0){
822 AliWarning(Form("The module %d (layer %i, ladder %i det %i ) is in acquisition, but no recpoints are present \n ",module, lay, lad, det));
823 emptymodules[lay-3]++;
824 empty++;
825
78fb0b3d 826 }
8a55554e 827 else
828 if(totside==2){
829 filledmodules[lay-3]++;
830 filled++;
831 }
832 else
833 if(totside==1)
834 {
835 // emptydr++;
836 //emptydriftregion[lay-3]++; //it has to be empty
837 emptydriftregion[lay-3]++;
838 emptydr++;
839 filleddr++;
840 filleddriftregion[lay-3]++;
841 }
842 if(totactiveside==1)
843 {
844 activedriftregion++;
845 activedrperlayer[lay-3]++;
846 }else if(totactiveside==2)
847 {
848 active++;
849 activemoduleperlayer[lay-3]++;
850 }
851
78fb0b3d 852 }
8a55554e 853 }//end for
854 AliInfo(Form("In total %d modules and %d single drift regions have recpoints.\n ",filled, filleddr));
855 AliInfo(Form("In total %d modules and %d single drift regions are empty\n",empty, emptydr));
856 for(Int_t i=0;i<2;i++)
857 {
858 AliInfo(Form("Layer %i \tempty modules %i \t filled modules %i\n", i+3,emptymodules[i], filledmodules[i]));
859 AliInfo(Form("Layer %i \tempty single drift regions %i \t filled single drift regions %i\n",i+3,emptydriftregion[i], filleddriftregion[i]));
860 }//end else layers
861 emptysum=emptymodules[0]+emptymodules[1];
862 emptydiff=emptysum-excluded;
863 emptyactivemoduleperlayer[0]=emptymodules[0]- excludedmoduleperlayer[0];
864 emptyactivemoduleperlayer[1]=emptymodules[1]- excludedmoduleperlayer[1];
865
866 emptydrsum=emptydriftregion[0]+emptydriftregion[1];
867 emptydrdiff=emptydrsum-excludeddriftregion;
868 emptyactivedrperlayer[0]=emptydriftregion[0]- excludeddrperlayer[0];
869 emptyactivedrperlayer[1]=emptydriftregion[1]- excludeddrperlayer[1];
870
62b7a70f 871 Int_t numlimit=1000;
8a55554e 872 if(emptysum>excluded||emptydrsum>excludeddriftregion){
62b7a70f 873 AliWarning(Form(" %i good module(s) and %i good single drift regions have not recpoints! \n",emptydiff,emptydrdiff));
874 AliWarning(Form(" Layer 3: %i good module(s) and %i good single drift regions have not recpoints! \n",emptyactivemoduleperlayer[0] ,emptyactivedrperlayer[0] ));
875 AliWarning(Form(" Layer 4: %i good module(s) and %i good single drift regions have not recpoints! \n",emptyactivemoduleperlayer[1] ,emptyactivedrperlayer[1] ));
8a55554e 876
62b7a70f 877 //sddQACheckerValue=fHighSDDValue[AliQAv1::kWARNING];
878 Int_t numlimits=1000;
879
8a55554e 880 results1.Form("%i good module(s) and %i good drift region(s) have not recpoints!",emptydiff,emptydrdiff);
62b7a70f 881 if((AliQAv1::Instance(AliQAv1::GetDetIndex("ITS")))->IsEventSpecieSet(AliRecoParam::kCosmic)==kFALSE)
8a55554e 882 {
62b7a70f 883 if(neventsrecpoints<numlimits)
884 {
885 results2.Form(" Events %d .Too few events.DO NOT CALL the Expert ",neventsrecpoints);
886 color=5;
887 sddQACheckerValue=fHighSDDValue[AliQAv1::kWARNING];
888 }
889 else
890 {
891 results2.Form(" Events %d .Follow the TWiki instruction and call the Expert ",neventsrecpoints);
892 color=2;
893 sddQACheckerValue=fHighSDDValue[AliQAv1::kERROR];
894 }
8a55554e 895 }
896 else
62b7a70f 897 if((AliQAv1::Instance(AliQAv1::GetDetIndex("ITS")))->IsEventSpecieSet(AliRecoParam::kCosmic)==kTRUE)
898 {
899 numlimit=10000;
900 if(neventsrecpoints<numlimit)
901 {
902 AliWarning(Form("This is a cosmic run. Some module and drift region are empty but all is OK. "));
903 results1.Form("OK. Thi is a cosmic run. You need a lot of events");
904 results2.Form("%i good module(s) and %i good drift are empty! DO NOT CALL THE EXPERT",emptydiff,emptydrdiff);
905 color=5;
906 sddQACheckerValue=fHighSDDValue[AliQAv1::kWARNING];
907 }
908 else
909 {
910 results1.Form("%i good module(s) and %i good drift region(s) have not recpoints!",emptydiff,emptydrdiff);
911 results2.Form(" Events %d .Follow the TWiki instruction and call the Expert ",neventsrecpoints);
912 color=2;
913 sddQACheckerValue=fHighSDDValue[AliQAv1::kERROR];
914 }
915 }
8a55554e 916 }
917
6bf3d857 918
8a55554e 919
920 if(exactive==0 && emptydiff==0 && exactivedriftregion==0 && emptydrdiff==0){
921 AliInfo(Form("All the active modules (%i) and single drift regions (%i) are in acquisition. The number of excluded modules are %i and the excluded single drift regions are %i\n",active,activedriftregion,excluded,excludeddriftregion));
922 results1.Form("OK.");
923 results2.Form(" All active modules have recpoints");
924 color=3;
62b7a70f 925 sddQACheckerValue=fHighSDDValue[AliQAv1::kINFO];
8a55554e 926 }
927 if(exactive!=0||exactivedriftregion!=0){
928 AliError(Form("%i modules and %i single drift regions excluded from the acquisition have recpoints. Active modules %i single drift region %i \n ",exactive,exactivedriftregion,active,activedriftregion));
929 AliError(Form("Layer 3: %i modules and %i single drift regions excluded from the acquisition have recpoints. Active modules %i single drift region %i \n ",exactivemoduleperlayer[0],exactivedrperlayer[0],activemoduleperlayer[0],activedrperlayer[0]));
930 AliError(Form("Layer 3: %i modules and %i single drift regions excluded from the acquisition have recpoints. Active modules %i single drift region %i \n ",exactivemoduleperlayer[1],exactivedrperlayer[1],activemoduleperlayer[1],activedrperlayer[1]));
931 results1.Form("%i modules and %i drift region excluded from the acquisition have recpoints",exactive,exactivedriftregion);
932 results2.Form("Follow the TWiki instructions and Call the SDD expert ");
933 color=2;
62b7a70f 934 sddQACheckerValue=fHighSDDValue[AliQAv1::kERROR];
8a55554e 935 }
936 if(excluded==exactive||excludeddriftregion==exactivedriftregion){
937 AliError(Form("All the modules (%d) or single drift regions (%d) excluded from the acquisition have recpoints!\n Active modules %i \t Active drfift regions %i\n",excluded,excludeddriftregion,active,activedriftregion));
938 results1.Form("All the modules (%d) or drift regions (%d) excluded from the acquisition have recpoints!",excluded,excludeddriftregion );
939 results2.Form("Follow the TWiki instructions and Call the SDD expert ");
940 color=6;
62b7a70f 941 sddQACheckerValue=fHighSDDValue[AliQAv1::kFATAL];
8a55554e 942 }
943 if(active==0||activedriftregion==0){
944 AliError(Form("No modules or single drift regions have recpoints: excluded %i \t excluded active %i \n\t\t excluded single drift regions %i \t excluded active drift regions %i \n", excluded, exactive, excludeddriftregion, exactivedriftregion));
945 results1.Form("No modules or drift region have recpoints: excluded modules %i excluded drift regions %i ", excluded, excludeddriftregion );
946 results2.Form("Follow the TWiki instructions and Call the SDD expert ");
947 color=6;
62b7a70f 948 sddQACheckerValue=fHighSDDValue[AliQAv1::kFATAL];
8a55554e 949 }
950
951 TPaveText *pave[2];
952 next.Begin();
953
954 while( (hdata=dynamic_cast<TH1* >(next())) )
955 {
956 if (hdata){
957 TString hname=hdata->GetName();
958 if(hname.Contains("SDDModPatternL3RP") || hname.Contains("SDDModPatternL4RP")){
959 if(hname.Contains("NORM"))continue;
960 //AliInfo("========================================Found histo 11\n");
961 Int_t lay=0;
962 if(hname.Contains("3"))lay=0;
963 else if(hname.Contains("4"))lay=1;
964 pave[lay]=new TPaveText(0.3,0.88,0.9,0.99,"NDC");
965 pave[lay]->AddText(results1.Data());
966 pave[lay]->AddText(results2.Data());
967 pave[lay]->SetFillColor(color);
968 pave[lay]->SetBorderSize(1);
969 pave[lay]->SetLineWidth(1);
970 hdata->GetListOfFunctions()->Add(pave[lay]);
971 }
972 else
973 if(hname.Contains("SDDRecPointCheck"))
974 {
975
976 //AliInfo("========================================Found histo\n");
977 ((TH1F*)hdata)->SetBinContent(5,active);
978 ((TH1F*)hdata)->SetBinContent(6,filled);
979 ((TH1F*)hdata)->SetBinContent(7,activedriftregion);
980 ((TH1F*)hdata)->SetBinContent(8,filleddr);
981 ((TH1F*)hdata)->SetBinContent(9,excluded);
982 ((TH1F*)hdata)->SetBinContent(10,empty);
983 ((TH1F*)hdata)->SetBinContent(11,excludeddriftregion);
984 ((TH1F*)hdata)->SetBinContent(12,emptydr);
985 ((TH1F*)hdata)->SetBinContent(13,exactive);
986 ((TH1F*)hdata)->SetBinContent(14,emptydiff);
987 ((TH1F*)hdata)->SetBinContent(15,exactivedriftregion);
988 ((TH1F*)hdata)->SetBinContent(16,emptydr);
989
990 //layer 3
991 ((TH1F*)hdata)->SetBinContent(19,activemoduleperlayer[0]);
992 ((TH1F*)hdata)->SetBinContent(20,filledmodules[0]);
993 ((TH1F*)hdata)->SetBinContent(21,activedrperlayer[0]);
994 ((TH1F*)hdata)->SetBinContent(22,filleddriftregion[0]);
995 ((TH1F*)hdata)->SetBinContent(23,excludedmoduleperlayer[0]);
996 ((TH1F*)hdata)->SetBinContent(24,emptymodules[0]);
997 ((TH1F*)hdata)->SetBinContent(25,excludeddrperlayer[0]);
998 ((TH1F*)hdata)->SetBinContent(26,emptydriftregion[0]);
999 ((TH1F*)hdata)->SetBinContent(27,exactivemoduleperlayer[0]);
1000 ((TH1F*)hdata)->SetBinContent(28,emptyactivemoduleperlayer[0]);
1001 ((TH1F*)hdata)->SetBinContent(29,activedrperlayer[0]);
1002 ((TH1F*)hdata)->SetBinContent(30,emptyactivedrperlayer[0]);
1003
1004 //layer 4
1005 ((TH1F*)hdata)->SetBinContent(35,activemoduleperlayer[1]);
1006 ((TH1F*)hdata)->SetBinContent(36,filledmodules[1]);
1007 ((TH1F*)hdata)->SetBinContent(37,activedrperlayer[1]);
1008 ((TH1F*)hdata)->SetBinContent(38,filleddriftregion[1]);
1009 ((TH1F*)hdata)->SetBinContent(39,excludedmoduleperlayer[1]);
1010 ((TH1F*)hdata)->SetBinContent(40,emptymodules[1]);
1011 ((TH1F*)hdata)->SetBinContent(41,excludeddrperlayer[1]);
1012 ((TH1F*)hdata)->SetBinContent(42,emptydriftregion[1]);
1013 ((TH1F*)hdata)->SetBinContent(43,exactivemoduleperlayer[1]);
1014 ((TH1F*)hdata)->SetBinContent(44,emptyactivemoduleperlayer[1]);
1015 ((TH1F*)hdata)->SetBinContent(45,activedrperlayer[1]);
1016 ((TH1F*)hdata)->SetBinContent(46,emptyactivedrperlayer[1]);
1017
1018 }
1019 }//if hadata
1020
1021 }//end while
1022
1023 }//end else
1024 delete hmodule;
1025 hmodule=NULL;
1026 for(Int_t i=0;i<2;i++) {
1027 delete hlayer[i];
1028 hlayer[i]=NULL;
1029 }//end for
1030
1031
62b7a70f 1032 //sddQACheckerValue=1.;
8a55554e 1033 }//end recpoint list uid = 20
8db5c9f9 1034 else if(uidrec==40)
1035 {
1036 //digitsr
1037 if (list->GetEntries() == 0){
62b7a70f 1038 sddQACheckerValue=fHighSDDValue[AliQAv1::kFATAL];
8db5c9f9 1039 break;
1040 } else{
1041
1042 while( (hdata = dynamic_cast<TH1* >(next())) ){
1043 if (hdata){
62b7a70f 1044 if(hdata->GetEntries()==0)sddQACheckerValue += fStepBitSDD[AliQAv1::kFATAL];
8db5c9f9 1045 else {
1046 TString hname=hdata->GetName();
1047 if(hname.Contains("SDD DIGITS Module Pattern")) {
1048 //see raws
1049
62b7a70f 1050 sddQACheckerValue += fStepBitSDD[AliQAv1::kINFO];
8db5c9f9 1051 } else if(hname.Contains("SDD Anode Distribution")) {
62b7a70f 1052 sddQACheckerValue += fStepBitSDD[AliQAv1::kINFO];
8db5c9f9 1053 } else if(hname.Contains("SDD Tbin Distribution")) {
1054 //to do as rp
62b7a70f 1055 sddQACheckerValue += fStepBitSDD[AliQAv1::kINFO];
8db5c9f9 1056 } else if(hname.Contains("SDD ADC Counts Distribution")) {
62b7a70f 1057 sddQACheckerValue += fStepBitSDD[AliQAv1::kINFO];
8db5c9f9 1058 }//end adc counts
1059
1060 }//end entries !=0
1061 }//end hdata
1062 }//end while
1063 }//end else
62b7a70f 1064 sddQACheckerValue=1.;
8db5c9f9 1065 }
1066
6bf3d857 1067 }
1068 break;
1069 case AliQAv1::kANA:
1070 {
80118b05 1071 AliInfo(Form("===================> No Check on %s\n",AliQAv1::GetAliTaskName(index)));
62b7a70f 1072 sddQACheckerValue=1.;
6bf3d857 1073 }
1074 break;
1075 case AliQAv1::kESD:
1076 {
80118b05 1077 AliInfo(Form("==================> No Check on %s\n",AliQAv1::GetAliTaskName(index)));
62b7a70f 1078 sddQACheckerValue=1.;
6bf3d857 1079 }
1080 break;
1081 case AliQAv1::kNTASK:{
1082 AliInfo(Form("==================> No Check on %s\n",AliQAv1::GetAliTaskName(index)));
62b7a70f 1083 sddQACheckerValue=1.;
6bf3d857 1084 }
1085 break;
1086 case AliQAv1::kSIM:{
1087 AliInfo(Form("Check on %s\n",AliQAv1::GetAliTaskName(index)));
1088 Int_t uid=list->GetUniqueID();
1089 if(uid==60) {
1090 //digits
1091 if (list->GetEntries() == 0){
62b7a70f 1092 sddQACheckerValue=fHighSDDValue[AliQAv1::kFATAL];
6bf3d857 1093 break;
1094 } else{
1095
1096 while( (hdata = dynamic_cast<TH1* >(next())) ){
1097 if (hdata){
62b7a70f 1098 if(hdata->GetEntries()==0)sddQACheckerValue += fStepBitSDD[AliQAv1::kFATAL];
6bf3d857 1099 else {
1100 TString hname=hdata->GetName();
1101 if(hname.Contains("SDDDIGITSModulePattern")) {
1102 //see raws
1103
62b7a70f 1104 sddQACheckerValue += fStepBitSDD[AliQAv1::kINFO];
6bf3d857 1105 } else if(hname.Contains("SDDAnodeDistribution")) {
62b7a70f 1106 sddQACheckerValue += fStepBitSDD[AliQAv1::kINFO];
6bf3d857 1107 } else if(hname.Contains("SDDTbinDistribution")) {
1108 //to do as rp
62b7a70f 1109 sddQACheckerValue += fStepBitSDD[AliQAv1::kINFO];
6bf3d857 1110 } else if(hname.Contains("SDDADCCountsDistribution")) {
62b7a70f 1111 sddQACheckerValue += fStepBitSDD[AliQAv1::kINFO];
6bf3d857 1112 }//end adc counts
1113
1114 }//end entries !=0
1115 }//end hdata
1116 }//end while
1117 }//end else
1118 } else if(uid==50)
1119 {
1120 //hits
1121 if (list->GetEntries() == 0){
62b7a70f 1122 sddQACheckerValue=fHighSDDValue[AliQAv1::kFATAL];
6bf3d857 1123 break;
1124 }
1125 else{
1126
1127 while( (hdata = dynamic_cast<TH1* >(next())) ){
1128 if (hdata){
62b7a70f 1129 if(hdata->GetEntries()==0)sddQACheckerValue += fStepBitSDD[AliQAv1::kFATAL];
6bf3d857 1130 else {
1131 TString hname=hdata->GetName();
1132 if(hname.Contains("SDDHITSModulePattern")) {
1133 //to do as raws
62b7a70f 1134 sddQACheckerValue += fStepBitSDD[AliQAv1::kINFO];
6bf3d857 1135 } else if(hname.Contains("SDDHITlenghtalonglocalYCoord")) {
62b7a70f 1136 sddQACheckerValue += fStepBitSDD[AliQAv1::kINFO];
6bf3d857 1137 } else if(hname.Contains("SDDHITlenghtalonglocalYCoordZoom")) {
62b7a70f 1138 sddQACheckerValue += fStepBitSDD[AliQAv1::kINFO];
6bf3d857 1139 } else if(hname.Contains("SDDDepositedEnergyDistribution")) {
62b7a70f 1140 sddQACheckerValue += fStepBitSDD[AliQAv1::kINFO];
6bf3d857 1141 }//end deposited energy
1142
1143 }//end entries !=0
1144 }//end hdata
1145 }//end while
1146 }//end else
1147 } else if(uid==70)
1148 {
1149 //sdigits
78fb0b3d 1150 if (list->GetEntries() == 0){
62b7a70f 1151 sddQACheckerValue=fHighSDDValue[AliQAv1::kFATAL];
6bf3d857 1152 break;
78fb0b3d 1153 } else{
6bf3d857 1154
78fb0b3d 1155 while( (hdata = dynamic_cast<TH1* >(next())) ){
1156 if (hdata){
62b7a70f 1157 if(hdata->GetEntries()==0)sddQACheckerValue += fStepBitSDD[AliQAv1::kFATAL];
78fb0b3d 1158 else {
1159 TString hname=hdata->GetName();
6bf3d857 1160 if(hname.Contains("SDDSDIGITSModulePattern")) {
1161 //to do as raws
62b7a70f 1162 sddQACheckerValue += fStepBitSDD[AliQAv1::kINFO];
78fb0b3d 1163 } else if(hname.Contains("SDDAnodeDistribution")) {
62b7a70f 1164 sddQACheckerValue += fStepBitSDD[AliQAv1::kINFO];
78fb0b3d 1165 } else if(hname.Contains("SDDTbinDistribution")) {
1166 //to do as rp
62b7a70f 1167 sddQACheckerValue += fStepBitSDD[AliQAv1::kINFO];
78fb0b3d 1168 } else if(hname.Contains("SDDADCCountsDistribution")) {
62b7a70f 1169 sddQACheckerValue += fStepBitSDD[AliQAv1::kINFO];
6bf3d857 1170 }//end adc counts bindistribution
78fb0b3d 1171 }//end entries !=0
1172 }//end hdata
1173 }//end while
1174 }//end else
6bf3d857 1175 }//end sdigits
62b7a70f 1176 sddQACheckerValue=1.;
6bf3d857 1177 }
1178 break;
1179
1180 }//end switch
1181
80118b05 1182 fCalibration=NULL;
d2f99642 1183 if(hdata) delete hdata;
6bf3d857 1184
1185
62b7a70f 1186 return sddQACheckerValue;
5b716935 1187}
6bf3d857 1188
3647765c 1189//__________________________________________________________________
80b9610c 1190void AliITSQASDDChecker::SetTaskOffset(Int_t taskoffset)
1191{
f09c0b42 1192 //set the number of the histograms already present in the list before the SDD histograms
80b9610c 1193 fSubDetOffset = taskoffset;
1194}
1195
1196
1197//__________________________________________________________________
f09c0b42 1198void AliITSQASDDChecker::SetStepBit(const Double_t *steprange)
80b9610c 1199{
f09c0b42 1200 //set the values of the step bit for each QA bit range calculated in the AliITSQAChecker class
80b9610c 1201 fStepBitSDD = new Double_t[AliQAv1::kNBIT];
1202 for(Int_t bit=0;bit<AliQAv1::kNBIT;bit++)
1203 {
1204 fStepBitSDD[bit]=steprange[bit];
1205 }
1206}
1207
1208//__________________________________________________________________
f09c0b42 1209void AliITSQASDDChecker::SetSDDLimits(const Float_t *lowvalue, const Float_t * highvalue)
3647765c 1210{
f09c0b42 1211 //set the low and high values in for each QA bit range calculated in the AliITSQAChecker class
80b9610c 1212 fLowSDDValue = new Float_t[AliQAv1::kNBIT];
1213 fHighSDDValue= new Float_t[AliQAv1::kNBIT];
1214
1215 for(Int_t bit=0;bit<AliQAv1::kNBIT;bit++)
1216 {
1217 fLowSDDValue[bit]=lowvalue[bit];
1218 fHighSDDValue[bit]= highvalue[bit];
1219 }
1220
3647765c 1221}
8db5c9f9 1222//__________________________________________________________________
1223Bool_t AliITSQASDDChecker::MakeSDDImage( TObjArray ** list, AliQAv1::TASKINDEX_t task, AliQAv1::MODE_t mode)
1224{
62b7a70f 1225 //create the image for raws and recpoints. In the other case, the default methodof CheckerBase class will be used
1226 //
8db5c9f9 1227 Bool_t rval=kFALSE;
1228 fImage=(TCanvas**)AliQAChecker::Instance()->GetDetQAChecker(0)->GetImage();
62b7a70f 1229
8db5c9f9 1230 switch(task)
1231 {
1232 case AliQAv1::kRAWS:{
1233 rval=MakeSDDRawsImage(list, task,mode);
1234 }
1235 break;
1236 case AliQAv1::kRECPOINTS:{ rval=MakeSDDRecPointsImage(list, task,mode); }
1237 break;
1238 case AliQAv1::kHITS:; case AliQAv1::kESDS:; case AliQAv1::kDIGITS:;case AliQAv1::kDIGITSR:;case AliQAv1::kSDIGITS:;case AliQAv1::kTRACKSEGMENTS:;case AliQAv1::kRECPARTICLES:; default:
1239 {
1240 rval=kFALSE;
1241 //AliQAChecker::Instance()->GetDetQAChecker(0)->MakeImage(list,task,mode);
1242 }
1243 break;
1244 case AliQAv1::kNULLTASKINDEX:; case AliQAv1::kNTASKINDEX:
1245 {AliWarning(Form("No histograms for this task ( %s ) \n", AliQAv1::GetTaskName(task).Data())); rval=kFALSE;}
1246 break;
1247 }
1248return rval;
1249}
1250
1251
1252//_______________________________________________________________________
1253Bool_t AliITSQASDDChecker::MakeSDDRawsImage(TObjArray ** list, AliQAv1::TASKINDEX_t task, AliQAv1::MODE_t mode )
1254{
62b7a70f 1255 // MakeSDDRawsImage: raw data QA plots
8db5c9f9 1256
1257 for (Int_t esIndex = 0 ; esIndex < AliRecoParam::kNSpecies ; esIndex++) {
1258 //printf("-------------------------> %i \n", esIndex);
1259 if (! AliQAv1::Instance(AliQAv1::GetDetIndex(GetName()))->IsEventSpecieSet(AliRecoParam::ConvertIndex(esIndex)) || list[esIndex]->GetEntries() == 0)
1260 {//printf ("Nothing for %s \n", AliRecoParam::GetEventSpecieName(esIndex));
1261 continue;
1262 }
1263 else{
1264 const Char_t * title = Form("QA_%s_%s_%s", GetName(), AliQAv1::GetTaskName(task).Data(), AliRecoParam::GetEventSpecieName(esIndex)) ;
1265 if ( !fImage[esIndex] ) {
1266 fImage[esIndex] = new TCanvas(title, title,1280,980) ;
1267 }
1268
1269 fImage[esIndex]->Clear() ;
1270 fImage[esIndex]->SetTitle(title) ;
1271 fImage[esIndex]->cd();
1272
1273 TPaveText someText(0.015, 0.015, 0.98, 0.98);
1274 someText.AddText(title);
1275 someText.Draw();
1276 fImage[esIndex]->Print(Form("%s%s%d.%s", AliQAv1::GetImageFileName(), AliQAv1::GetModeName(mode), AliQAChecker::Instance()->GetRunNumber(), AliQAv1::GetImageFileFormat()), "ps") ;
1277 fImage[esIndex]->Clear() ;
1278 Int_t nx =2; //TMath::Nint(TMath::Sqrt(nImages));
8a55554e 1279 Int_t ny =2; // nx ;
8db5c9f9 1280 //if (nx < TMath::Sqrt(nImages))
1281 //ny++ ;
1282 fImage[esIndex]->Divide(nx, ny) ;
1283 TIter nexthist(list[esIndex]) ;
1284 TH1* hist = NULL ;
1285 Int_t npad = 1 ;
1286 fImage[esIndex]->cd(npad);
1287 fImage[esIndex]->cd(npad)->SetBorderMode(0) ;
1288 while ( (hist=static_cast<TH1*>(nexthist())) ) {
1289 //gPad=fImage[esIndex]->cd(npad)->GetPad(npad);
1290 TString cln(hist->ClassName()) ;
1291 if ( ! cln.Contains("TH") )
1292 continue ;
1293
1294 if(hist->TestBit(AliQAv1::GetImageBit())) {
1295 hist->GetXaxis()->SetTitleSize(0.02);
1296 hist->GetYaxis()->SetTitleSize(0.02);
1297 hist->GetXaxis()->SetLabelSize(0.02);
1298 hist->GetYaxis()->SetLabelSize(0.02);
1299 if(cln.Contains("TH2"))
1300 {
1301 gPad->SetRightMargin(0.15);
1302 gPad->SetLeftMargin(0.05);
1303 hist->SetStats(0);
1304 hist->SetOption("colz") ;
1305 //hist->GetListOfFunctions()->FindObject("palette")->SetLabelSize(0.025);
1306 //gPad->Update();
1307 }
1308 hist->DrawCopy() ;
1309 fImage[esIndex]->cd(++npad) ;
1310 fImage[esIndex]->cd(npad)->SetBorderMode(0) ;
1311 }
1312 }
1313 fImage[esIndex]->Print(Form("%s%s%d.%s", AliQAv1::GetImageFileName(), AliQAv1::GetModeName(mode), AliQAChecker::Instance()->GetRunNumber(), AliQAv1::GetImageFileFormat()), "ps") ;
1314 }
1315 }
1316 return kTRUE;
1317}
1318
1319
1320
1321
1322//_______________________________________________________________________
1323Bool_t AliITSQASDDChecker::MakeSDDRecPointsImage(TObjArray ** list, AliQAv1::TASKINDEX_t task, AliQAv1::MODE_t mode )
1324{
62b7a70f 1325 // MakeSDDRecPointsImage: rec point QA plots
80b9610c 1326
8db5c9f9 1327 for (Int_t esIndex = 0 ; esIndex < AliRecoParam::kNSpecies ; esIndex++) {
1328 if (! AliQAv1::Instance(AliQAv1::GetDetIndex(GetName()))->IsEventSpecieSet(AliRecoParam::ConvertIndex(esIndex)) || list[esIndex]->GetEntries() == 0)
1329 {
1330 //printf ("Nothing for %s \n", AliQAv1::GetTaskName(task).Data());
1331 continue;
1332 }
1333 const Char_t * title = Form("QA_%s_%s_%s", GetName(), AliQAv1::GetTaskName(task).Data(), AliRecoParam::GetEventSpecieName(esIndex)) ;
1334 if ( !fImage[esIndex] ) {
1335 fImage[esIndex] = new TCanvas(title, title,1280,980) ;
1336 }
1337 fImage[esIndex]->Clear() ;
1338 fImage[esIndex]->SetTitle(title) ;
1339 fImage[esIndex]->cd();
1340 fImage[esIndex]->SetBorderMode(0) ;
1341 TPaveText someText(0.015, 0.015, 0.98, 0.98);
1342 someText.AddText(title);
1343 someText.Draw();
1344 fImage[esIndex]->Print(Form("%s%s%d.%s", AliQAv1::GetImageFileName(), AliQAv1::GetModeName(mode), AliQAChecker::Instance()->GetRunNumber(), AliQAv1::GetImageFileFormat()), "ps") ;
1345 fImage[esIndex]->Clear() ;
1346 Int_t nx =2; //TMath::Nint(TMath::Sqrt(nImages));
1347 Int_t ny =4; // nx ;
1348 //if (nx < TMath::Sqrt(nImages))
1349 //ny++ ;
1350 fImage[esIndex]->Divide(nx, ny) ;
1351 TIter nexthist(list[esIndex]) ;
1352 TH1* hist = NULL ;
1353 Int_t npad = 1 ;
1354 fImage[esIndex]->cd(npad) ;
1355 fImage[esIndex]->cd(npad)->SetBorderMode(0) ;
1356 while ( (hist=static_cast<TH1*>(nexthist())) ) {
1357 //gPad=fImage[esIndex]->cd(npad)->GetPad(npad);
1358 TString cln(hist->ClassName()) ;
1359 //printf("=====================> Class name %s \n",cln.Data());
1360 if ( ! cln.Contains("TH") )
1361 continue ;
1362 if(hist->TestBit(AliQAv1::GetImageBit())) {
1363 hist->GetXaxis()->SetTitleSize(0.02);
1364 hist->GetYaxis()->SetTitleSize(0.02);
1365 hist->GetXaxis()->SetLabelSize(0.02);
1366 hist->GetYaxis()->SetLabelSize(0.02);
1367 if(cln.Contains("TH1"))
1368 {
1369 hist->SetFillColor(kOrange+7);
1370 //SetFrameFillColor(kAzure-9);
1371 //hist->DrawCopy() ;
1372 }
1373 if(cln.Contains("TH2"))
1374 {
1375 gPad->SetRightMargin(0.15);
1376 gPad->SetLeftMargin(0.05);
1377 hist->SetStats(0);
1378 hist->SetOption("colz") ;
1379 // TPaletteAxis *paletta =(TPaletteAxis*)hist->GetListOfFunctions()->FindObject("palette");
1380 //paletta->SetLabelSize(0.025);
1381 //gPad->Update();
1382 }
1383 hist->DrawCopy();
1384 fImage[esIndex]->cd(++npad) ;
1385 fImage[esIndex]->cd(npad)->SetBorderMode(0) ;
1386 }
1387 }
1388 fImage[esIndex]->Print(Form("%s%s%d.%s", AliQAv1::GetImageFileName(), AliQAv1::GetModeName(mode), AliQAChecker::Instance()->GetRunNumber(), AliQAv1::GetImageFileFormat()), "ps") ;
1389 }
1390 // }
1391 return kTRUE;
1392}