Bug fix. Removed delete statement
[u/mrichter/AliRoot.git] / STEER / AliQADataMakerSim.cxx
CommitLineData
04236e67 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
16
17/* $Id$ */
18
202374b1 19//
20// Base Class
21// Produces the data needed to calculate the quality assurance.
22// All data must be mergeable objects.
23// Y. Schutz CERN July 2007
24//
04236e67 25
26// --- ROOT system ---
04236e67 27#include <TFile.h>
04236e67 28#include <TTree.h>
29#include <TClonesArray.h>
30
31// --- Standard library ---
32
33// --- AliRoot header files ---
34#include "AliLog.h"
35#include "AliQADataMakerSim.h"
04236e67 36
37ClassImp(AliQADataMakerSim)
38
39//____________________________________________________________________________
40AliQADataMakerSim::AliQADataMakerSim(const char * name, const char * title) :
41 AliQADataMaker(name, title),
eca4fa66 42 fDigitsQAList(NULL),
43 fHitsQAList(NULL),
6252ceeb 44 fSDigitsQAList(NULL),
45 fHitsArray(NULL),
46 fSDigitsArray(NULL)
04236e67 47{
48 // ctor
49 fDetectorDirName = GetName() ;
50}
51
52//____________________________________________________________________________
53AliQADataMakerSim::AliQADataMakerSim(const AliQADataMakerSim& qadm) :
54 AliQADataMaker(qadm.GetName(), qadm.GetTitle()),
55 fDigitsQAList(qadm.fDigitsQAList),
56 fHitsQAList(qadm.fHitsQAList),
6252ceeb 57 fSDigitsQAList(qadm.fSDigitsQAList),
58 fHitsArray(NULL),
59 fSDigitsArray(NULL)
04236e67 60{
61 //copy ctor
62 fDetectorDirName = GetName() ;
63}
64
63c6f8ae 65//____________________________________________________________________________
66AliQADataMakerSim::~AliQADataMakerSim()
67{
68 //dtor: delete the TObjArray and thei content
7ff8385d 69 if ( fDigitsQAList ) {
57acd2d2 70 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
71 if ( fDigitsQAList[specie]->IsOwner() )
72 fDigitsQAList[specie]->Delete() ;
73 }
eca4fa66 74 delete[] fDigitsQAList ;
75 }
7ff8385d 76 if ( fHitsQAList ) {
57acd2d2 77 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
78 if ( fHitsQAList[specie]->IsOwner() )
79 fHitsQAList[specie]->Delete() ;
80 }
eca4fa66 81 delete[] fHitsQAList ;
82 }
7ff8385d 83 if ( fSDigitsQAList ) {
57acd2d2 84 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
85 if ( fSDigitsQAList[specie]->IsOwner() )
86 fSDigitsQAList[specie]->Delete() ;
87 }
eca4fa66 88 delete[] fSDigitsQAList ;
89 }
6252ceeb 90 if (fHitsArray) {
91 fHitsArray->Clear() ;
92 delete fHitsArray ;
93 }
94 if (fSDigitsArray) {
95 fSDigitsArray->Clear() ;
96 delete fSDigitsArray ;
97 }
63c6f8ae 98}
99
04236e67 100//__________________________________________________________________
101AliQADataMakerSim& AliQADataMakerSim::operator = (const AliQADataMakerSim& qadm )
102{
103 // Assignment operator.
104 this->~AliQADataMakerSim();
105 new(this) AliQADataMakerSim(qadm);
106 return *this;
107}
108
109//____________________________________________________________________________
930e6e3e 110void AliQADataMakerSim::EndOfCycle()
111{
112 // Finishes a cycle of QA for all tasks
4e25ac79 113 EndOfCycle(AliQAv1::kHITS) ;
114 EndOfCycle(AliQAv1::kSDIGITS) ;
115 EndOfCycle(AliQAv1::kDIGITS) ;
930e6e3e 116 ResetCycle() ;
117}
118
119//____________________________________________________________________________
4e25ac79 120void AliQADataMakerSim::EndOfCycle(AliQAv1::TASKINDEX_t task)
04236e67 121{
122 // Finishes a cycle of QA data acquistion
57acd2d2 123 TObjArray ** list = NULL ;
04236e67 124
4e25ac79 125 if ( task == AliQAv1::kHITS )
04236e67 126 list = fHitsQAList ;
4e25ac79 127 else if ( task == AliQAv1::kSDIGITS )
04236e67 128 list = fSDigitsQAList ;
4e25ac79 129 else if ( task == AliQAv1::kDIGITS )
04236e67 130 list = fDigitsQAList ;
930e6e3e 131
132 if ( ! list )
133 return ;
04236e67 134 EndOfDetectorCycle(task, list) ;
eca4fa66 135 fDetectorDir = fOutput->GetDirectory(GetDetectorDirName()) ;
136 if (!fDetectorDir)
137 fDetectorDir = fOutput->mkdir(GetDetectorDirName()) ;
138 TDirectory * subDir = fDetectorDir->GetDirectory(AliQAv1::GetTaskName(task)) ;
139 if (!subDir)
140 subDir = fDetectorDir->mkdir(AliQAv1::GetTaskName(task)) ;
141 subDir->cd() ;
142 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
143 if (! AliQAv1::Instance(AliQAv1::GetDetIndex(GetName()))->IsEventSpecieSet(AliRecoParam::ConvertIndex(specie)) )
144 continue ;
145 TDirectory * eventSpecieDir = subDir->GetDirectory(AliRecoParam::GetEventSpecieName(specie)) ;
146 if (!eventSpecieDir)
147 eventSpecieDir = subDir->mkdir(AliRecoParam::GetEventSpecieName(specie)) ;
148 eventSpecieDir->cd() ;
149 TIter next(list[specie]) ;
150 TObject * obj ;
151 while ( (obj = next()) ) {
152 if (!obj->TestBit(AliQAv1::GetExpertBit()))
153 obj->Write() ;
154 }
155 if (WriteExpert()) {
156 TDirectory * expertDir = eventSpecieDir->GetDirectory(AliQAv1::GetExpert()) ;
157 if (!expertDir)
158 expertDir = eventSpecieDir->mkdir(AliQAv1::GetExpert()) ;
159 expertDir->cd() ;
160 next.Reset() ;
161 while ( (obj = next()) ) {
162 if (!obj->TestBit(AliQAv1::GetExpertBit()))
163 continue ;
164 obj->Write() ;
165 }
6e65319c 166 }
57acd2d2 167 fOutput->Save() ;
6e65319c 168 }
04236e67 169}
eca4fa66 170
04236e67 171//____________________________________________________________________________
4e25ac79 172void AliQADataMakerSim::Exec(AliQAv1::TASKINDEX_t task, TObject * data)
04236e67 173{
174 // creates the quality assurance data for the various tasks (Hits, SDigits, Digits, ESDs)
eca4fa66 175
4e25ac79 176 if ( task == AliQAv1::kHITS ) {
5379c4a3 177 AliDebug(AliQAv1::GetQADebugLevel(), "Processing Hits QA") ;
eca4fa66 178 if (strcmp(data->ClassName(), "TClonesArray") == 0) {
6252ceeb 179 fHitsArray = static_cast<TClonesArray *>(data) ;
180 MakeHits() ;
eca4fa66 181 } else if (strcmp(data->ClassName(), "TTree") == 0) {
182 TTree * tree = static_cast<TTree *>(data) ;
183 MakeHits(tree) ;
184 } else {
185 AliWarning("data are neither a TClonesArray nor a TTree") ;
186 }
4e25ac79 187 } else if ( task == AliQAv1::kSDIGITS ) {
5379c4a3 188 AliDebug(AliQAv1::GetQADebugLevel(), "Processing SDigits QA") ;
eca4fa66 189 if (strcmp(data->ClassName(), "TClonesArray") == 0) {
6252ceeb 190 fSDigitsArray = static_cast<TClonesArray *>(data) ;
191 MakeSDigits() ;
eca4fa66 192 } else if (strcmp(data->ClassName(), "TTree") == 0) {
193 TTree * tree = static_cast<TTree *>(data) ;
194 MakeSDigits(tree) ;
195 } else {
196 AliWarning("data are neither a TClonesArray nor a TTree") ;
197 }
198 } else if ( task == AliQAv1::kDIGITS ) {
5379c4a3 199 AliDebug(AliQAv1::GetQADebugLevel(), "Processing Digits QA") ;
eca4fa66 200 if (strcmp(data->ClassName(), "TClonesArray") == 0) {
6252ceeb 201 fDigitsArray = static_cast<TClonesArray *>(data) ;
202 MakeDigits() ;
eca4fa66 203 } else if (strcmp(data->ClassName(), "TTree") == 0) {
204 TTree * tree = static_cast<TTree *>(data) ;
205 MakeDigits(tree) ;
206 } else {
207 AliWarning("data are neither a TClonesArray nor a TTree") ;
208 }
209 }
04236e67 210}
211
212//____________________________________________________________________________
4e25ac79 213TObjArray ** AliQADataMakerSim::Init(AliQAv1::TASKINDEX_t task, Int_t cycles)
04236e67 214{
215 // general intialisation
216
04236e67 217 if (cycles > 0)
218 SetCycle(cycles) ;
57acd2d2 219 TObjArray ** rv = NULL ;
4e25ac79 220 if ( task == AliQAv1::kHITS ) {
63c6f8ae 221 if ( ! fHitsQAList ) {
57acd2d2 222 fHitsQAList = new TObjArray *[AliRecoParam::kNSpecies] ;
223 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
49466ea2 224 fHitsQAList[specie] = new TObjArray(AliQAv1::GetMaxQAObj()) ;
4e25ac79 225 fHitsQAList[specie]->SetName(Form("%s_%s_%s", GetName(), AliQAv1::GetTaskName(task).Data(), AliRecoParam::GetEventSpecieName(specie))) ;
57acd2d2 226 }
63c6f8ae 227 }
04236e67 228 rv = fHitsQAList ;
4e25ac79 229 } else if ( task == AliQAv1::kSDIGITS ) {
63c6f8ae 230 if ( ! fSDigitsQAList ) {
57acd2d2 231 fSDigitsQAList = new TObjArray *[AliRecoParam::kNSpecies] ;
232 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
49466ea2 233 fSDigitsQAList[specie] = new TObjArray(AliQAv1::GetMaxQAObj()) ;
4e25ac79 234 fSDigitsQAList[specie]->SetName(Form("%s_%s_%s", GetName(), AliQAv1::GetTaskName(task).Data(), AliRecoParam::GetEventSpecieName(specie))) ;
57acd2d2 235 }
63c6f8ae 236 }
04236e67 237 rv = fSDigitsQAList ;
4e25ac79 238 } else if ( task == AliQAv1::kDIGITS ) {
63c6f8ae 239 if ( ! fDigitsQAList ) {
57acd2d2 240 fDigitsQAList = new TObjArray *[AliRecoParam::kNSpecies] ;
241 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
49466ea2 242 fDigitsQAList[specie] = new TObjArray(AliQAv1::GetMaxQAObj()) ;
4e25ac79 243 fDigitsQAList[specie]->SetName(Form("%s_%s_%s", GetName(), AliQAv1::GetTaskName(task).Data(), AliRecoParam::GetEventSpecieName(specie))) ;
57acd2d2 244 }
63c6f8ae 245 }
04236e67 246 rv = fDigitsQAList ;
247 }
248
249 return rv ;
eca4fa66 250}
04236e67 251
252//____________________________________________________________________________
4e25ac79 253void AliQADataMakerSim::Init(AliQAv1::TASKINDEX_t task, TObjArray ** list, Int_t run, Int_t cycles)
04236e67 254{
255 // Intialisation by passing the list of QA data booked elsewhere
256
257 fRun = run ;
258 if (cycles > 0)
259 SetCycle(cycles) ;
260
4e25ac79 261 if ( task == AliQAv1::kHITS ) {
04236e67 262 fHitsQAList = list ;
4e25ac79 263 } else if ( task == AliQAv1::kSDIGITS) {
04236e67 264 fSDigitsQAList = list ;
4e25ac79 265 } else if ( task == AliQAv1::kDIGITS ) {
04236e67 266 fDigitsQAList = list ;
267 }
268}
269
270//____________________________________________________________________________
930e6e3e 271void AliQADataMakerSim::StartOfCycle(Int_t run)
272{
273 // Finishes a cycle of QA for all tasks
274 Bool_t samecycle = kFALSE ;
4e25ac79 275 StartOfCycle(AliQAv1::kHITS, run, samecycle) ;
930e6e3e 276 samecycle = kTRUE ;
4e25ac79 277 StartOfCycle(AliQAv1::kSDIGITS, run, samecycle) ;
278 StartOfCycle(AliQAv1::kDIGITS, run, samecycle) ;
930e6e3e 279}
280
281//____________________________________________________________________________
4e25ac79 282void AliQADataMakerSim::StartOfCycle(AliQAv1::TASKINDEX_t task, Int_t run, const Bool_t sameCycle)
04236e67 283{
284 // Finishes a cycle of QA data acquistion
930e6e3e 285 if ( run > 0 )
286 fRun = run ;
04236e67 287 if ( !sameCycle || fCurrentCycle == -1) {
288 ResetCycle() ;
289 if (fOutput)
290 fOutput->Close() ;
4e25ac79 291 fOutput = AliQAv1::GetQADataFile(GetName(), fRun) ;
04236e67 292 }
293
5379c4a3 294 AliDebug(AliQAv1::GetQADebugLevel(), Form(" Run %d Cycle %d task %s file %s",
4e25ac79 295 fRun, fCurrentCycle, AliQAv1::GetTaskName(task).Data(), fOutput->GetName() )) ;
04236e67 296
eca4fa66 297 //fDetectorDir = fOutput->GetDirectory(GetDetectorDirName()) ;
298// if (!fDetectorDir)
299// fDetectorDir = fOutput->mkdir(GetDetectorDirName()) ;
300//
301// TDirectory * subDir = fDetectorDir->GetDirectory(AliQAv1::GetTaskName(task)) ;
302// if (!subDir)
303// subDir = fDetectorDir->mkdir(AliQAv1::GetTaskName(task)) ;
304//
305// for ( Int_t index = AliRecoParam::kDefault ; index < AliRecoParam::kNSpecies ; index++ ) {
306// TDirectory * eventSpecieDir = subDir->GetDirectory(AliRecoParam::GetEventSpecieName(index)) ;
307// if (!eventSpecieDir)
308// eventSpecieDir = subDir->mkdir(AliRecoParam::GetEventSpecieName(index)) ;
309// TDirectory * expertDir = eventSpecieDir->GetDirectory(AliQAv1::GetExpert()) ;
310// if (!expertDir)
311// expertDir = eventSpecieDir->mkdir(AliQAv1::GetExpert()) ;
312// }
04236e67 313 StartOfDetectorCycle() ;
314}