1 /**************************************************************************
2 * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
18 // *************************************************************
19 // Checks the quality assurance
20 // by comparing with reference data
22 // -------------------------------------------------------------
23 // W. Ferrarese + P. Cerello Feb 2008
26 // --- ROOT system ---
29 // --- Standard library ---
31 // --- AliRoot header files ---
32 #include "AliITSQADataMakerSim.h"
33 #include "AliITSQASPDDataMakerSim.h"
34 #include "AliITSQASDDDataMakerSim.h"
35 #include "AliITSQASSDDataMakerSim.h"
38 #include "AliQAChecker.h"
39 #include "AliITSQAChecker.h"
40 #include "AliRawReader.h"
42 ClassImp(AliITSQADataMakerSim)
44 //____________________________________________________________________________
45 AliITSQADataMakerSim::AliITSQADataMakerSim(Short_t subDet) :
46 AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kITS), "ITS Quality Assurance Data Maker"),
52 //ctor used to discriminate OnLine-Offline analysis
53 if(fSubDetector < 0 || fSubDetector > 3) {
54 AliError("Error: fSubDetector number out of range; return\n");
57 // Initialization for each subdetector
58 if(fSubDetector == 0 || fSubDetector == 1) {
59 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SPD DataMakerSim\n");
60 fSPDDataMaker = new AliITSQASPDDataMakerSim(this);
62 if(fSubDetector == 0 || fSubDetector == 2) {
63 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SDD DataMakerSim\n");
64 //printf("AliITSQADM::Create SDD DataMakerSim\n");
65 fSDDDataMaker = new AliITSQASDDDataMakerSim(this);
67 if(fSubDetector == 0 || fSubDetector == 3) {
68 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SSD DataMakerSim\n");
69 fSSDDataMaker = new AliITSQASSDDataMakerSim(this);
73 //____________________________________________________________________________
74 AliITSQADataMakerSim::~AliITSQADataMakerSim(){
76 if(fSPDDataMaker)delete fSPDDataMaker;
77 if(fSDDDataMaker)delete fSDDDataMaker;
78 if(fSSDDataMaker)delete fSSDDataMaker;
81 //____________________________________________________________________________
82 AliITSQADataMakerSim::AliITSQADataMakerSim(const AliITSQADataMakerSim& qadm) :
84 fSubDetector(qadm.fSubDetector),
90 SetName((const char*)qadm.GetName()) ;
91 SetTitle((const char*)qadm.GetTitle());
94 //__________________________________________________________________
95 AliITSQADataMakerSim& AliITSQADataMakerSim::operator = (const AliITSQADataMakerSim& qac )
98 this->~AliITSQADataMakerSim();
99 new(this) AliITSQADataMakerSim(qac);
103 //____________________________________________________________________________
104 void AliITSQADataMakerSim::StartOfDetectorCycle()
106 //Detector specific actions at start of cycle
107 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Start of ITS Cycle\n");
108 ResetEventTrigClasses(); // reset triggers list to select all histos
111 if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->StartOfDetectorCycle();
112 if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->StartOfDetectorCycle();
113 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->StartOfDetectorCycle();
116 //____________________________________________________________________________
117 void AliITSQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray** list)
120 AliInfo(Form("End of Dedetctor Cycle called for %s\n",AliQAv1::GetTaskName(task).Data() ));
121 ResetEventTrigClasses(); // reset triggers list to select all histos
123 // launch the QA checking
124 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
125 if (! IsValidEventSpecie(specie, list) ) continue;
126 SetEventSpecie(AliRecoParam::ConvertIndex(specie));
128 Int_t idnumber=list[specie]->GetUniqueID();
129 //printf("specie %s \t id number == %d\n",AliRecoParam::GetEventSpecieName(specie),idnumber);
131 //AliInfo(Form("No check for %s\n",AliQAv1::GetTaskName(task).Data() ))
133 } //skip kDigitsR and not filled TobjArray specie
135 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSDM instantiates checker with Run(AliQAv1::kITS, task, list)\n");
136 if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->EndOfDetectorCycle(task, list); //[specie]);
137 if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->EndOfDetectorCycle(task, list); //[specie]);
138 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->EndOfDetectorCycle(task, list); //[specie]);
140 AliQAChecker *qac = AliQAChecker::Instance();
141 AliITSQAChecker *qacb = (AliITSQAChecker *) qac->GetDetQAChecker(0);
142 Int_t subdet=GetSubDet();
143 qacb->SetSubDet(subdet);
146 qacb->SetTaskOffset(fSPDDataMaker->GetOffset(task,specie),fSDDDataMaker->GetOffset(task,specie),fSSDDataMaker->GetOffset(task,specie)); //Setting the offset for the QAChecker list
147 qacb->SetHisto(fSPDDataMaker->GetTaskHisto(task), fSDDDataMaker->GetTaskHisto(task), fSSDDataMaker->GetTaskHisto(task));
151 Int_t offset=GetDetTaskOffset(subdet, task);
152 qacb->SetDetTaskOffset(subdet,offset);
153 Int_t histo=GetDetTaskHisto(subdet, task);
154 qacb->SetDetHisto(subdet,histo);
156 qac->Run( AliQAv1::kITS , task, list);
157 }//end else unique id
161 //____________________________________________________________________________
162 void AliITSQADataMakerSim::InitDigits()
165 // Initialization for Digits data
166 fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(60);
167 if(fSubDetector == 0 || fSubDetector == 1) {
168 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitDigits\n");
170 fSPDDataMaker->InitDigits();
172 if(fSubDetector == 0 || fSubDetector == 2) {
173 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitDigits\n");
175 fSDDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
176 fSDDDataMaker->InitDigits();
178 if(fSubDetector == 0 || fSubDetector == 3) {
179 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitDigits\n");
181 fSSDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
182 fSSDDataMaker->InitDigits();
185 ClonePerTrigClass(AliQAv1::kDIGITS); // this should be the last line
188 //____________________________________________________________________________
189 void AliITSQADataMakerSim::MakeDigits()
191 // Fill QA for digits
192 if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeDigits();
193 if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeDigits();
194 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits();
196 IncEvCountCycleDigits();
197 IncEvCountTotalDigits();
201 //____________________________________________________________________________
202 void AliITSQADataMakerSim::MakeDigits(TTree * digits)
204 // Fill QA for digits
205 if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeDigits(digits);
206 if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeDigits(digits);
207 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits(digits);
209 IncEvCountCycleDigits();
210 IncEvCountTotalDigits();
214 //____________________________________________________________________________
215 void AliITSQADataMakerSim::InitSDigits()
217 // Initialization for SDigits
218 fSDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(70);
219 if(fSubDetector == 0 || fSubDetector == 1) {
220 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitSDigits\n");
222 fSPDDataMaker->InitSDigits();
224 if(fSubDetector == 0 || fSubDetector == 2){
225 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitSDigits\n");
227 fSDDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
228 fSDDDataMaker->InitSDigits();
230 if(fSubDetector == 0 || fSubDetector == 3) {
231 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitSDigits\n");
233 fSSDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
234 fSSDDataMaker->InitSDigits();
237 ClonePerTrigClass(AliQAv1::kSDIGITS); // this should be the last line
240 //____________________________________________________________________________
241 void AliITSQADataMakerSim::MakeSDigits()
243 // Fill QA for sdigits
244 if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeSDigits();
245 if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeSDigits();
246 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits();
248 IncEvCountCycleSDigits();
249 IncEvCountTotalSDigits();
253 //____________________________________________________________________________
254 void AliITSQADataMakerSim::MakeSDigits(TTree * sdigits)
256 // Fill QA for recpoints
257 if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeSDigits(sdigits);
258 if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeSDigits(sdigits);
259 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits(sdigits);
261 IncEvCountCycleSDigits();
262 IncEvCountTotalSDigits();
266 //____________________________________________________________________________
267 void AliITSQADataMakerSim::InitHits()
269 // Initialization for hits
270 fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(50);
271 if(fSubDetector == 0 || fSubDetector == 1) {
272 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitHits\n");
273 fSPDDataMaker->InitHits();
275 if(fSubDetector == 0 || fSubDetector == 2) {
276 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitHits\n");
278 fSDDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
279 fSDDDataMaker->InitHits();
281 if(fSubDetector == 0 || fSubDetector == 3) {
282 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitHits\n");
284 fSSDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
285 fSSDDataMaker->InitHits();
288 ClonePerTrigClass(AliQAv1::kHITS); // this should be the last line
291 //____________________________________________________________________________
292 void AliITSQADataMakerSim::MakeHits()
295 if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeHits();
296 if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeHits();
297 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits();
299 IncEvCountCycleHits();
300 IncEvCountTotalHits();
304 //____________________________________________________________________________
305 void AliITSQADataMakerSim::MakeHits(TTree * hits)
308 if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeHits(hits);
309 if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeHits(hits);
310 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits(hits);
312 IncEvCountCycleHits();
313 IncEvCountTotalHits();
317 //_________________________________________________________________
318 Int_t AliITSQADataMakerSim::GetDetTaskOffset(Int_t subdet,AliQAv1::TASKINDEX_t task)
321 //return the offset for each subdetector
326 offset=fSPDDataMaker->GetOffset(task);
330 offset=fSDDDataMaker->GetOffset(task);
334 offset=fSSDDataMaker->GetOffset(task);
338 AliWarning("No specific subdetector (SPD, SDD, SSD) selected!! Offset set to zero \n");
346 //_________________________________________________________________
347 Int_t AliITSQADataMakerSim::GetDetTaskHisto(Int_t subdet,AliQAv1::TASKINDEX_t task)
349 //return of the number of histograms for each task and for each sub detector
355 histo=fSPDDataMaker->GetOffset(task);
359 histo=fSDDDataMaker->GetOffset(task);
363 histo=fSSDDataMaker->GetOffset(task);
367 AliWarning("No specific subdetector (SPD, SDD, SSD) selected!! Offset set to zero \n");