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");
109 if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->StartOfDetectorCycle();
110 if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->StartOfDetectorCycle();
111 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->StartOfDetectorCycle();
114 //____________________________________________________________________________
115 void AliITSQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray** list)
118 AliInfo(Form("End of Dedetctor Cycle called for %s\n",AliQAv1::GetTaskName(task).Data() ));
120 // launch the QA checking
121 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
122 if (! IsValidEventSpecie(specie, list) ){
126 Int_t idnumber=list[specie]->GetUniqueID();
127 //printf("specie %s \t id number == %d\n",AliRecoParam::GetEventSpecieName(specie),idnumber);
130 //AliInfo(Form("No check for %s\n",AliQAv1::GetTaskName(task).Data() ))
132 } //skip kDigitsR and not filled TobjArray specie
134 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSDM instantiates checker with Run(AliQAv1::kITS, task, list)\n");
135 if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->EndOfDetectorCycle(task, list[specie]);
136 if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->EndOfDetectorCycle(task, list[specie]);
137 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
158 }//end else event specie
162 //____________________________________________________________________________
163 void AliITSQADataMakerSim::InitDigits()
166 fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(60);
167 // Initialization for Digits data
168 if(fSubDetector == 0 || fSubDetector == 1) {
169 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitDigits\n");
171 fSPDDataMaker->InitDigits();
173 if(fSubDetector == 0 || fSubDetector == 2) {
174 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitDigits\n");
176 fSDDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
177 fSDDDataMaker->InitDigits();
179 if(fSubDetector == 0 || fSubDetector == 3) {
180 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitDigits\n");
182 fSSDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
183 fSSDDataMaker->InitDigits();
187 //____________________________________________________________________________
188 void AliITSQADataMakerSim::MakeDigits()
190 // Fill QA for digits
191 if(fSubDetector == 0 || fSubDetector == 1)
192 fSPDDataMaker->MakeDigits() ;
195 if(fSubDetector == 0 || fSubDetector == 2)
196 fSDDDataMaker->MakeDigits() ;
198 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits();
201 //____________________________________________________________________________
202 void AliITSQADataMakerSim::MakeDigits(TTree * digits)
204 // Fill QA for digits
205 if(fSubDetector == 0 || fSubDetector == 1)
206 fSPDDataMaker->MakeDigits(digits) ;
208 if(fSubDetector == 0 || fSubDetector == 2)
209 fSDDDataMaker->MakeDigits(digits) ;
211 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits(digits);
214 //____________________________________________________________________________
215 void AliITSQADataMakerSim::InitSDigits()
217 fSDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(70);
218 // Initialization for RECPOINTS
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();
238 //____________________________________________________________________________
239 void AliITSQADataMakerSim::MakeSDigits()
241 // Fill QA for recpoints
242 if(fSubDetector == 0 || fSubDetector == 1)
243 fSPDDataMaker->MakeSDigits() ;
246 if(fSubDetector == 0 || fSubDetector == 2)
247 fSDDDataMaker->MakeSDigits() ;
250 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits();
253 //____________________________________________________________________________
254 void AliITSQADataMakerSim::MakeSDigits(TTree * sdigits)
256 // Fill QA for recpoints
257 if(fSubDetector == 0 || fSubDetector == 1){
258 fSPDDataMaker->MakeSDigits(sdigits) ;
261 if(fSubDetector == 0 || fSubDetector == 2){
262 fSDDDataMaker->MakeSDigits(sdigits) ;
265 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits(sdigits);
268 //____________________________________________________________________________
269 void AliITSQADataMakerSim::InitHits()
271 fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(50);
272 // Initialization for hits
273 if(fSubDetector == 0 || fSubDetector == 1) {
274 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitHits\n");
275 fSPDDataMaker->InitHits();
277 if(fSubDetector == 0 || fSubDetector == 2) {
278 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitHits\n");
280 fSDDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
281 fSDDDataMaker->InitHits();
283 if(fSubDetector == 0 || fSubDetector == 3) {
284 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitHits\n");
286 fSSDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
287 fSSDDataMaker->InitHits();
291 //____________________________________________________________________________
292 void AliITSQADataMakerSim::MakeHits()
295 if(fSubDetector == 0 || fSubDetector == 1) {
296 fSPDDataMaker->MakeHits() ;
299 if(fSubDetector == 0 || fSubDetector == 2) {
300 fSDDDataMaker->MakeHits() ;
303 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits();
306 //____________________________________________________________________________
307 void AliITSQADataMakerSim::MakeHits(TTree * hits)
310 if(fSubDetector == 0 || fSubDetector == 1) {
311 fSPDDataMaker->MakeHits(hits) ;
313 if(fSubDetector == 0 || fSubDetector == 2) {
314 fSDDDataMaker->MakeHits(hits) ;
317 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits(hits);
320 //_________________________________________________________________
321 Int_t AliITSQADataMakerSim::GetDetTaskOffset(Int_t subdet,AliQAv1::TASKINDEX_t task)
328 offset=fSPDDataMaker->GetOffset(task);
332 offset=fSDDDataMaker->GetOffset(task);
336 offset=fSSDDataMaker->GetOffset(task);
340 AliWarning("No specific subdetector (SPD, SDD, SSD) selected!! Offset set to zero \n");
348 //_________________________________________________________________
349 Int_t AliITSQADataMakerSim::GetDetTaskHisto(Int_t subdet,AliQAv1::TASKINDEX_t task)
356 histo=fSPDDataMaker->GetOffset(task);
360 histo=fSDDDataMaker->GetOffset(task);
364 histo=fSSDDataMaker->GetOffset(task);
368 AliWarning("No specific subdetector (SPD, SDD, SSD) selected!! Offset set to zero \n");