1 /**************************************************************************
2 * Copyright(c) 1998-1999, 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 **************************************************************************/
19 #include <TPluginManager.h>
24 #include "AliESDEvent.h"
27 #include "AliQADataMaker.h"
28 #include "AliQADataMakerSteer.h"
29 #include "AliRawReaderRoot.h"
31 #include "AliRunLoader.h"
33 ClassImp(AliQADataMakerSteer)
35 //_____________________________________________________________________________
36 AliQADataMakerSteer::AliQADataMakerSteer(const char* gAliceFilename, const char * name, const char * title) :
42 fGAliceFileName(gAliceFilename),
49 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
50 fLoader[iDet] = NULL ;
51 fQADataMaker[iDet] = NULL ;
52 fQACycles[iDet] = 999999 ;
56 //_____________________________________________________________________________
57 AliQADataMakerSteer::AliQADataMakerSteer(const AliQADataMakerSteer & qas) :
63 fGAliceFileName(qas.fGAliceFileName),
64 fRunNumber(qas.fRunNumber),
65 fNumberOfEvents(qas.fNumberOfEvents),
70 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
71 fLoader[iDet] = qas.fLoader[iDet] ;
72 fQADataMaker[iDet] = qas.fQADataMaker[iDet] ;
73 fQACycles[iDet] = qas.fQACycles[iDet] ;
77 //_____________________________________________________________________________
78 AliQADataMakerSteer & AliQADataMakerSteer::operator = (const AliQADataMakerSteer & qas)
80 // assignment operator
81 this->~AliQADataMakerSteer() ;
82 new(this) AliQADataMakerSteer(qas) ;
86 //_____________________________________________________________________________
87 AliQADataMakerSteer::~AliQADataMakerSteer()
90 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
92 if (fQADataMaker[iDet]) {
93 (fQADataMaker[iDet])->Finish() ;
94 delete fQADataMaker[iDet] ;
95 fQADataMaker[iDet] = NULL ;
104 //_____________________________________________________________________________
105 AliLoader * AliQADataMakerSteer::GetLoader(Int_t iDet)
107 // get the loader for a detector
109 TString detName = AliQA::GetDetName(iDet) ;
110 fLoader[iDet] = fRunLoader->GetLoader(detName + "Loader");
112 return fLoader[iDet] ;
114 // load the QA data maker object
115 TPluginManager* pluginManager = gROOT->GetPluginManager() ;
116 TString loaderName = "Ali" + detName + "Loader" ;
118 AliLoader * loader = NULL ;
119 // first check if a plugin is defined for the quality assurance data maker
120 TPluginHandler* pluginHandler = pluginManager->FindHandler("AliLoader", detName) ;
121 // if not, add a plugin for it
122 if (!pluginHandler) {
123 AliDebug(1, Form("defining plugin for %s", loaderName.Data())) ;
124 TString libs = gSystem->GetLibraries() ;
125 if (libs.Contains("lib" + detName + "base.so") || (gSystem->Load("lib" + detName + "base.so") >= 0)) {
126 pluginManager->AddHandler("AliQADataMaker", detName, loaderName, detName + "loader", loaderName + "()") ;
128 pluginManager->AddHandler("AliLoader", detName, loaderName, detName, loaderName + "()") ;
130 pluginHandler = pluginManager->FindHandler("AliLoader", detName) ;
132 if (pluginHandler && (pluginHandler->LoadPlugin() == 0)) {
133 loader = (AliLoader *) pluginHandler->ExecPlugin(0) ;
136 fLoader[iDet] = loader ;
140 //_____________________________________________________________________________
141 AliQADataMaker * AliQADataMakerSteer::GetQADataMaker(Int_t iDet)
143 // get the quality assurance data maker for a detector
145 if (fQADataMaker[iDet])
146 return fQADataMaker[iDet] ;
148 AliQADataMaker * qadm = NULL ;
151 // load the QA data maker object
152 TPluginManager* pluginManager = gROOT->GetPluginManager() ;
153 TString detName = AliQA::GetDetName(iDet) ;
154 TString qadmName = "Ali" + detName + "QADataMaker" ;
156 // first check if a plugin is defined for the quality assurance data maker
157 TPluginHandler* pluginHandler = pluginManager->FindHandler("AliQADataMaker", detName) ;
158 // if not, add a plugin for it
159 if (!pluginHandler) {
160 AliDebug(1, Form("defining plugin for %s", qadmName.Data())) ;
161 TString libs = gSystem->GetLibraries() ;
162 if (libs.Contains("lib" + detName + "base.so") || (gSystem->Load("lib" + detName + "base.so") >= 0)) {
163 pluginManager->AddHandler("AliQADataMaker", detName, qadmName, detName + "qadm", qadmName + "()") ;
165 pluginManager->AddHandler("AliQADataMaker", detName, qadmName, detName, qadmName + "()") ;
167 pluginHandler = pluginManager->FindHandler("AliQADataMaker", detName) ;
169 if (pluginHandler && (pluginHandler->LoadPlugin() == 0)) {
170 qadm = (AliQADataMaker *) pluginHandler->ExecPlugin(0) ;
173 fQADataMaker[iDet] = qadm ;
178 //_____________________________________________________________________________
179 Bool_t AliQADataMakerSteer::Init(const AliQA::TASKINDEX taskIndex, const char * fileName )
181 // Initialize the event source and QA data makers
183 if (taskIndex == AliQA::kRAWS) {
184 fRawReader = new AliRawReaderRoot(fileName) ;
187 fRawReader->NextEvent() ;
188 fRunNumber = fRawReader->GetRunNumber() ;
189 fRawReader->RewindEvents();
190 fNumberOfEvents = 999999 ;
192 } else if (taskIndex == AliQA::kESDS) {
193 if (!gSystem->AccessPathName("AliESDs.root")) { // AliESDs.root exists
194 TFile * esdFile = TFile::Open("AliESDs.root") ;
195 fESDTree = dynamic_cast<TTree *> (esdFile->Get("esdTree")) ;
196 fESD = new AliESDEvent() ;
197 fESD->ReadFromTree(fESDTree) ;
198 fESDTree->GetEntry(0) ;
199 fRunNumber = fESD->GetRunNumber() ;
200 fNumberOfEvents = fESDTree->GetEntries() ;
202 AliError("AliESDs.root not found") ;
207 if ( !InitRunLoader() ) {
208 AliError("Problems in getting the Run Loader") ;
211 fRunNumber = fRunLoader->GetAliRun()->GetRunNumber() ;
212 fNumberOfEvents = fRunLoader->GetNumberOfEvents() ;
216 // Initialize all QA data makers for all detectors
217 for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
218 AliQADataMaker * qadm = GetQADataMaker(iDet) ;
220 AliWarning(Form("AliQADataMaker not found for %s", AliQA::GetDetName(iDet))) ;
222 AliInfo(Form("Data Maker found for %s", qadm->GetName())) ;
223 qadm->Init(taskIndex, fRunNumber, GetQACycles(iDet)) ;
224 qadm->StartOfCycle(taskIndex, fCycleSame) ;
231 //_____________________________________________________________________________
232 Bool_t AliQADataMakerSteer::InitRunLoader()
234 // get or create the run loader
240 if (!gSystem->AccessPathName(fGAliceFileName.Data())) { // galice.root exists
241 // load all base libraries to get the loader classes
242 TString libs = gSystem->GetLibraries() ;
243 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
244 TString detName = AliQA::GetDetName(iDet) ;
245 if (detName == "HLT")
247 if (libs.Contains("lib" + detName + "base.so"))
249 gSystem->Load("lib" + detName + "base.so");
251 fRunLoader = AliRunLoader::Open(fGAliceFileName.Data());
253 AliError(Form("no run loader found in file %s", fGAliceFileName.Data()));
256 fRunLoader->CdGAFile();
257 if (fRunLoader->LoadgAlice() == 0) {
258 gAlice = fRunLoader->GetAliRun();
261 AliError(Form("no gAlice object found in file %s", fGAliceFileName.Data()));
265 } else { // galice.root does not exist
266 AliError(Form("the file %s does not exist", fGAliceFileName.Data()));
273 //_____________________________________________________________________________
274 Bool_t AliQADataMakerSteer::Finish(const AliQA::TASKINDEX taskIndex)
276 // write output to file for all detectors
277 for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
278 AliQADataMaker * qadm = GetQADataMaker(iDet) ;
280 qadm->EndOfCycle(taskIndex) ;
286 //_____________________________________________________________________________
287 void AliQADataMakerSteer::Reset()
289 // Reset the default data members
290 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
291 fLoader[iDet] = NULL;
292 if (fQADataMaker[iDet]) {
293 (fQADataMaker[iDet])->Reset() ;
294 // delete fQADataMaker[iDet] ;
295 // fQADataMaker[iDet] = NULL ;
303 fCycleSame = kFALSE ;
307 fNumberOfEvents = 0 ;
310 //_____________________________________________________________________________
311 Bool_t AliQADataMakerSteer::Run(const AliQA::TASKINDEX taskIndex, const char * fileName )
313 // Runs all the QA data Maker for every detector
316 if ( !Init(taskIndex, fileName) )
319 // Fill QA data in event loop
320 for (UInt_t iEvent = 0 ; iEvent < fNumberOfEvents ; iEvent++) {
322 AliInfo(Form("processing event %d", iEvent));
323 if ( taskIndex == AliQA::kRAWS ) {
324 if ( !fRawReader->NextEvent() )
326 } else if ( taskIndex == AliQA::kESDS ) {
327 fESDTree->GetEntry(iEvent) ;
329 fRunLoader->GetEvent(iEvent);
331 // loop over detectors
332 for (UInt_t iDet = 0 ; iDet < fgkNDetectors ; iDet++) {
333 AliQADataMaker * qadm = GetQADataMaker(iDet) ;
337 if ( qadm->IsCycleDone() ) {
338 qadm->EndOfCycle(AliQA::kRAWS) ;
339 qadm->StartOfCycle(AliQA::kRAWS) ;
344 qadm->Exec(taskIndex, fRawReader) ;
347 GetLoader(iDet)->LoadHits() ;
348 data = GetLoader(iDet)->TreeH() ;
350 AliWarning(Form(" Hit Tree not found for %s", AliQA::GetDetName(iDet))) ;
352 qadm->Exec(taskIndex, data) ;
355 case AliQA::kSDIGITS :
356 GetLoader(iDet)->LoadSDigits() ;
357 data = GetLoader(iDet)->TreeS() ;
359 AliWarning(Form(" SDigit Tree not found for %s", AliQA::GetDetName(iDet))) ;
361 qadm->Exec(taskIndex, data) ;
364 case AliQA::kDIGITS :
365 GetLoader(iDet)->LoadDigits() ;
366 data = GetLoader(iDet)->TreeD() ;
368 AliWarning(Form(" Digit Tree not found for %s", AliQA::GetDetName(iDet))) ;
370 qadm->Exec(taskIndex, data) ;
373 case AliQA::kRECPOINTS :
374 GetLoader(iDet)->LoadRecPoints() ;
375 data = GetLoader(iDet)->TreeR() ;
377 AliWarning(Form("RecPoints not found for %s", AliQA::GetDetName(iDet))) ;
379 qadm->Exec(taskIndex, data) ;
382 case AliQA::kTRACKSEGMENTS :
384 case AliQA::kRECPARTICLES :
387 qadm->Exec(taskIndex, fESD) ;
394 // Save QA data for all detectors
395 rv = Finish(taskIndex) ;