]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ITS/AliITSQADataMakerSim.cxx
fix requested by Laurent Aphecetche to clean up the memory properly during simulation...
[u/mrichter/AliRoot.git] / ITS / AliITSQADataMakerSim.cxx
CommitLineData
04236e67 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// contained in a DB
22// -------------------------------------------------------------
8dc34e37 23// W. Ferrarese + P. Cerello Feb 2008
04236e67 24// INFN Torino
25
26// --- ROOT system ---
04236e67 27#include <TTree.h>
57acd2d2 28#include <TMath.h>
04236e67 29// --- Standard library ---
30
31// --- AliRoot header files ---
32#include "AliITSQADataMakerSim.h"
8dc34e37 33#include "AliITSQASPDDataMakerSim.h"
34#include "AliITSQASDDDataMakerSim.h"
35#include "AliITSQASSDDataMakerSim.h"
04236e67 36#include "AliLog.h"
4e25ac79 37#include "AliQAv1.h"
04236e67 38#include "AliQAChecker.h"
5dfa9b71 39#include "AliITSQAChecker.h"
8dc34e37 40#include "AliRawReader.h"
04236e67 41
42ClassImp(AliITSQADataMakerSim)
43
44//____________________________________________________________________________
8dc34e37 45AliITSQADataMakerSim::AliITSQADataMakerSim(Short_t subDet) :
4e25ac79 46AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kITS), "ITS Quality Assurance Data Maker"),
8dc34e37 47fSubDetector(subDet),
48fSPDDataMaker(NULL),
49fSDDDataMaker(NULL),
50fSSDDataMaker(NULL)
51{
52 //ctor used to discriminate OnLine-Offline analysis
53 if(fSubDetector < 0 || fSubDetector > 3) {
54 AliError("Error: fSubDetector number out of range; return\n");
55 }
56
57 // Initialization for RAW data
58 if(fSubDetector == 0 || fSubDetector == 1) {
5379c4a3 59 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SPD DataMakerSim\n");
8dc34e37 60 fSPDDataMaker = new AliITSQASPDDataMakerSim(this);
61 }
62 if(fSubDetector == 0 || fSubDetector == 2) {
5379c4a3 63 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SDD DataMakerSim\n");
8dc34e37 64 //printf("AliITSQADM::Create SDD DataMakerSim\n");
65 fSDDDataMaker = new AliITSQASDDDataMakerSim(this);
66 }
67 if(fSubDetector == 0 || fSubDetector == 3) {
5379c4a3 68 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SSD DataMakerSim\n");
8dc34e37 69 fSSDDataMaker = new AliITSQASSDDataMakerSim(this);
70 }
04236e67 71}
72
73//____________________________________________________________________________
8dc34e37 74AliITSQADataMakerSim::~AliITSQADataMakerSim(){
75 // destructor
76 if(fSPDDataMaker)delete fSPDDataMaker;
77 if(fSDDDataMaker)delete fSDDDataMaker;
78 if(fSSDDataMaker)delete fSSDDataMaker;
04236e67 79}
80
81//____________________________________________________________________________
82AliITSQADataMakerSim::AliITSQADataMakerSim(const AliITSQADataMakerSim& qadm) :
8dc34e37 83AliQADataMakerSim(),
84fSubDetector(qadm.fSubDetector),
85fSPDDataMaker(NULL),
86fSDDDataMaker(NULL),
87fSSDDataMaker(NULL)
04236e67 88{
89 //copy ctor
90 SetName((const char*)qadm.GetName()) ;
8dc34e37 91 SetTitle((const char*)qadm.GetTitle());
04236e67 92}
93
94//__________________________________________________________________
95AliITSQADataMakerSim& AliITSQADataMakerSim::operator = (const AliITSQADataMakerSim& qac )
96{
97 // Equal operator.
98 this->~AliITSQADataMakerSim();
99 new(this) AliITSQADataMakerSim(qac);
100 return *this;
101}
102
103//____________________________________________________________________________
85f5e9c2 104void AliITSQADataMakerSim::StartOfDetectorCycle()
04236e67 105{
106 //Detector specific actions at start of cycle
5379c4a3 107 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Start of ITS Cycle\n");
8dc34e37 108
109 if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->StartOfDetectorCycle();
8dc34e37 110 if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->StartOfDetectorCycle();
111 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->StartOfDetectorCycle();
04236e67 112}
113
114//____________________________________________________________________________
4e25ac79 115void AliITSQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray** list)
04236e67 116{
117 // launch the QA checking
57acd2d2 118 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
6252ceeb 119 if (! IsValidEventSpecie(specie, list) )
120 continue ;
5379c4a3 121 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSDM instantiates checker with Run(AliQAv1::kITS, task, list)\n");
57acd2d2 122 if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->EndOfDetectorCycle(task, list[specie]);
123 if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->EndOfDetectorCycle(task, list[specie]);
124 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->EndOfDetectorCycle(task, list[specie]);
04236e67 125
8bb2fdfd 126
127 AliQAChecker *qac = AliQAChecker::Instance();
128 AliITSQAChecker *qacb = (AliITSQAChecker *) qac->GetDetQAChecker(0);
129 Int_t subdet=GetSubDet();
130 qacb->SetSubDet(subdet);
614c7e17 131
8bb2fdfd 132 if(subdet== 0 ){
133 qacb->SetTaskOffset(fSPDDataMaker->GetOffset(task),fSDDDataMaker->GetOffset(task),fSSDDataMaker->GetOffset(task)); //Setting the offset for the QAChecker list
614c7e17 134 }
8bb2fdfd 135 else
136 if(subdet!=0){
137 Int_t offset=GetDetTaskOffset(subdet, task);
138 qacb->SetDetTaskOffset(subdet,offset);
139 }
4e25ac79 140 qac->Run( AliQAv1::kITS , task, list);
8bb2fdfd 141 }
04236e67 142}
143
8dc34e37 144//____________________________________________________________________________
145void AliITSQADataMakerSim::InitDigits()
146{
147 // Initialization for RAW data
148 if(fSubDetector == 0 || fSubDetector == 1) {
5379c4a3 149 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitDigits\n");
4a903927 150
8dc34e37 151 fSPDDataMaker->InitDigits();
152 }
153 if(fSubDetector == 0 || fSubDetector == 2) {
5379c4a3 154 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitDigits\n");
4a903927 155
156 fSDDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
8dc34e37 157 fSDDDataMaker->InitDigits();
158 }
159 if(fSubDetector == 0 || fSubDetector == 3) {
5379c4a3 160 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitDigits\n");
4a903927 161
162 fSSDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
8dc34e37 163 fSSDDataMaker->InitDigits();
164 }
165}
166
167//____________________________________________________________________________
6252ceeb 168void AliITSQADataMakerSim::MakeDigits()
8dc34e37 169{
eca4fa66 170 // Fill QA for digits
4a903927 171 if(fSubDetector == 0 || fSubDetector == 1)
172 fSPDDataMaker->MakeDigits() ;
173
eca4fa66 174
4a903927 175 if(fSubDetector == 0 || fSubDetector == 2)
176 fSDDDataMaker->MakeDigits() ;
eca4fa66 177
6252ceeb 178 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits();
8dc34e37 179}
180
181//____________________________________________________________________________
182void AliITSQADataMakerSim::MakeDigits(TTree * digits)
183{
eca4fa66 184 // Fill QA for digits
4a903927 185 if(fSubDetector == 0 || fSubDetector == 1)
186 fSPDDataMaker->MakeDigits(digits) ;
187
188 if(fSubDetector == 0 || fSubDetector == 2)
189 fSDDDataMaker->MakeDigits(digits) ;
eca4fa66 190
8dc34e37 191 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits(digits);
192}
193
194//____________________________________________________________________________
195void AliITSQADataMakerSim::InitSDigits()
196{
197 // Initialization for RECPOINTS
198 if(fSubDetector == 0 || fSubDetector == 1) {
5379c4a3 199 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitSDigits\n");
4a903927 200
8dc34e37 201 fSPDDataMaker->InitSDigits();
202 }
4a903927 203 if(fSubDetector == 0 || fSubDetector == 2){
5379c4a3 204 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitSDigits\n");
4a903927 205
206 fSDDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
207 fSDDDataMaker->InitSDigits();
8dc34e37 208 }
209 if(fSubDetector == 0 || fSubDetector == 3) {
5379c4a3 210 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitSDigits\n");
4a903927 211
212 fSSDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
213 fSSDDataMaker->InitSDigits();
8dc34e37 214 }
215}
216
217//____________________________________________________________________________
6252ceeb 218void AliITSQADataMakerSim::MakeSDigits()
8dc34e37 219{
220 // Fill QA for recpoints
4a903927 221 if(fSubDetector == 0 || fSubDetector == 1)
222 fSPDDataMaker->MakeSDigits() ;
223
eca4fa66 224
4a903927 225 if(fSubDetector == 0 || fSubDetector == 2)
226 fSDDDataMaker->MakeSDigits() ;
227
eca4fa66 228
6252ceeb 229 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits();
8dc34e37 230}
231
232//____________________________________________________________________________
233void AliITSQADataMakerSim::MakeSDigits(TTree * sdigits)
234{
235 // Fill QA for recpoints
eca4fa66 236 if(fSubDetector == 0 || fSubDetector == 1){
4a903927 237 fSPDDataMaker->MakeSDigits(sdigits) ;
238 }
eca4fa66 239
240 if(fSubDetector == 0 || fSubDetector == 2){
4a903927 241 fSDDDataMaker->MakeSDigits(sdigits) ;
eca4fa66 242 }
243
8dc34e37 244 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits(sdigits);
245}
246
247//____________________________________________________________________________
248void AliITSQADataMakerSim::InitHits()
249{
eca4fa66 250 // Initialization for hits
8dc34e37 251 if(fSubDetector == 0 || fSubDetector == 1) {
5379c4a3 252 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitHits\n");
8dc34e37 253 fSPDDataMaker->InitHits();
254 }
255 if(fSubDetector == 0 || fSubDetector == 2) {
5379c4a3 256 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitHits\n");
4a903927 257
258 fSDDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
259 fSDDDataMaker->InitHits();
8dc34e37 260 }
261 if(fSubDetector == 0 || fSubDetector == 3) {
5379c4a3 262 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitHits\n");
4a903927 263
264 fSSDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
265 fSSDDataMaker->InitHits();
8dc34e37 266 }
267}
268
269//____________________________________________________________________________
6252ceeb 270void AliITSQADataMakerSim::MakeHits()
8dc34e37 271{
eca4fa66 272 // Fill QA for hits
273 if(fSubDetector == 0 || fSubDetector == 1) {
4a903927 274 fSPDDataMaker->MakeHits() ;
275 }
eca4fa66 276
277 if(fSubDetector == 0 || fSubDetector == 2) {
4a903927 278 fSDDDataMaker->MakeHits() ;
eca4fa66 279 }
280
6252ceeb 281 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits();
8dc34e37 282}
283
284//____________________________________________________________________________
285void AliITSQADataMakerSim::MakeHits(TTree * hits)
286{
eca4fa66 287 // Fill QA for hits
288 if(fSubDetector == 0 || fSubDetector == 1) {
4a903927 289 fSPDDataMaker->MakeHits(hits) ;
290 }
eca4fa66 291 if(fSubDetector == 0 || fSubDetector == 2) {
4a903927 292 fSDDDataMaker->MakeHits(hits) ;
eca4fa66 293 }
294
8dc34e37 295 if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits(hits);
296}
614c7e17 297
298//_________________________________________________________________
4e25ac79 299Int_t AliITSQADataMakerSim::GetDetTaskOffset(Int_t subdet,AliQAv1::TASKINDEX_t task)
614c7e17 300{
301 switch(subdet)
302 {
303
304 Int_t offset;
305 case 1:
306 offset=fSPDDataMaker->GetOffset(task);
307 return offset;
308 break;
309 case 2:
310 offset=fSDDDataMaker->GetOffset(task);
311 return offset;
312 break;
313 case 3:
314 offset=fSSDDataMaker->GetOffset(task);
315 return offset;
316 break;
317 default:
318 AliWarning("No specific subdetector (SPD, SDD, SSD) selected!! Offset set to zero \n");
319 offset=0;
320 return offset;
321 break;
322 }
323 //return offset;
324}