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 **************************************************************************/
16 /* $Id: AliQAManager.cxx 30894 2009-02-05 13:46:48Z schutz $ */
17 ///////////////////////////////////////////////////////////////////////////////
19 // class for running the QA makers //
21 // AliQAManager qas; //
22 // qas.Run(AliQAv1::kRAWS, rawROOTFileName); //
23 // qas.Run(AliQAv1::kHITS); //
24 // qas.Run(AliQAv1::kSDIGITS); //
25 // qas.Run(AliQAv1::kDIGITS); //
26 // qas.Run(AliQAv1::kRECPOINTS); //
27 // qas.Run(AliQAv1::kESDS); //
29 ///////////////////////////////////////////////////////////////////////////////
34 #include <TFileMerger.h>
36 #include <TGridCollection.h>
37 #include <TGridResult.h>
38 #include <TPluginManager.h>
42 #include <TStopwatch.h>
44 #include "AliCDBManager.h"
45 #include "AliCDBEntry.h"
47 #include "AliCDBMetaData.h"
48 #include "AliCodeTimer.h"
49 #include "AliCorrQADataMakerRec.h"
50 #include "AliDetectorRecoParam.h"
51 #include "AliESDEvent.h"
52 #include "AliGeomManager.h"
53 #include "AliGlobalQADataMaker.h"
54 #include "AliHeader.h"
56 #include "AliModule.h"
58 #include "AliQAChecker.h"
59 #include "AliQACheckerBase.h"
60 #include "AliQADataMakerRec.h"
61 #include "AliQADataMakerSim.h"
62 #include "AliQAManager.h"
63 #include "AliRawReaderDate.h"
64 #include "AliRawReaderFile.h"
65 #include "AliRawReaderRoot.h"
67 #include "AliRunLoader.h"
68 #include "AliRunTag.h"
71 ClassImp(AliQAManager)
72 AliQAManager* AliQAManager::fgQAInstance = 0x0;
74 //_____________________________________________________________________________
75 AliQAManager::AliQAManager() :
88 fNumberOfEvents(999999),
92 fRawReaderDelete(kTRUE),
95 fEventSpecie(AliRecoParam::kDefault),
100 fMaxEvents = fNumberOfEvents ;
101 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
102 if (IsSelected(AliQAv1::GetDetName(iDet))) {
103 fLoader[iDet] = NULL ;
104 fQADataMaker[iDet] = NULL ;
105 fQACycles[iDet] = 999999 ;
111 //_____________________________________________________________________________
112 AliQAManager::AliQAManager(AliQAv1::MODE_t mode, const Char_t* gAliceFilename) :
121 fGAliceFileName(gAliceFilename),
124 fMode(AliQAv1::GetModeName(mode)),
125 fNumberOfEvents(999999),
129 fRawReaderDelete(kTRUE),
132 fEventSpecie(AliRecoParam::kDefault),
137 fMaxEvents = fNumberOfEvents ;
138 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
139 if (IsSelected(AliQAv1::GetDetName(iDet))) {
140 fLoader[iDet] = NULL ;
141 fQADataMaker[iDet] = NULL ;
142 fQACycles[iDet] = 999999 ;
148 //_____________________________________________________________________________
149 AliQAManager::AliQAManager(const AliQAManager & qas) :
151 fCurrentEvent(qas.fCurrentEvent),
153 fDetectors(qas.fDetectors),
154 fDetectorsW(qas.fDetectorsW),
158 fGAliceFileName(qas.fGAliceFileName),
159 fFirstEvent(qas.fFirstEvent),
160 fMaxEvents(qas.fMaxEvents),
162 fNumberOfEvents(qas.fNumberOfEvents),
164 fRunNumber(qas.fRunNumber),
166 fRawReaderDelete(kTRUE),
169 fEventSpecie(qas.fEventSpecie),
170 fPrintImage(qas.fPrintImage),
171 fSaveData(qas.fSaveData)
175 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
176 fLoader[iDet] = qas.fLoader[iDet] ;
177 fQADataMaker[iDet] = qas.fQADataMaker[iDet] ;
178 fQACycles[iDet] = qas.fQACycles[iDet] ;
179 fQAWriteExpert[iDet] = qas.fQAWriteExpert[iDet] ;
183 //_____________________________________________________________________________
184 AliQAManager & AliQAManager::operator = (const AliQAManager & qas)
186 // assignment operator
187 this->~AliQAManager() ;
188 new(this) AliQAManager(qas) ;
192 //_____________________________________________________________________________
193 AliQAManager::~AliQAManager()
196 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
197 if (IsSelected(AliQAv1::GetDetName(iDet))) {
198 fLoader[iDet] = NULL;
199 if (fQADataMaker[iDet]) {
200 (fQADataMaker[iDet])->Finish() ;
201 delete fQADataMaker[iDet] ;
205 if (fRawReaderDelete) {
211 //_____________________________________________________________________________
212 Bool_t AliQAManager::DoIt(const AliQAv1::TASKINDEX_t taskIndex)
214 // Runs all the QA data Maker for every detector
217 // Fill QA data in event loop
218 for (UInt_t iEvent = fFirstEvent ; iEvent < (UInt_t)fMaxEvents ; iEvent++) {
221 if ( iEvent%10 == 0 )
222 AliDebug(AliQAv1::GetQADebugLevel(), Form("processing event %d", iEvent));
223 if ( taskIndex == AliQAv1::kRAWS ) {
224 if ( !fRawReader->NextEvent() )
226 } else if ( taskIndex == AliQAv1::kESDS ) {
227 if ( fESDTree->GetEntry(iEvent) == 0 )
230 if ( fRunLoader->GetEvent(iEvent) != 0 )
233 // loop over active loaders
236 detList = GetEventInfo()->GetTriggerCluster() ;
237 for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
238 if (!detList.IsNull() && !detList.Contains(AliQAv1::GetDetName(iDet)))
240 if (IsSelected(AliQAv1::GetDetName(iDet)) ){
241 AliQADataMaker * qadm = GetQADataMaker(iDet) ;
242 if (!qadm) continue; // This detector doesn't have any QA (for example, HLT)
243 if ( qadm->IsCycleDone() ) {
244 qadm->EndOfCycle(taskIndex) ;
246 TTree * data = NULL ;
247 AliLoader* loader = GetLoader(qadm->GetUniqueID());
249 case AliQAv1::kNULLTASKINDEX :
251 case AliQAv1::kRAWS :
252 qadm->Exec(taskIndex, fRawReader) ;
254 case AliQAv1::kHITS :
257 data = loader->TreeH() ;
259 AliWarning(Form(" Hit Tree not found for %s", AliQAv1::GetDetName(iDet))) ;
262 qadm->Exec(taskIndex, data) ;
265 case AliQAv1::kSDIGITS :
268 TString fileName(Form("%s.SDigits.root", AliQAv1::GetDetName(iDet))) ;
269 if (gSystem->FindFile("./", fileName)) {
271 loader->LoadSDigits() ;
272 data = loader->TreeS() ;
274 AliWarning(Form(" SDigit Tree not found for %s", AliQAv1::GetDetName(iDet))) ;
277 qadm->Exec(taskIndex, data) ;
282 case AliQAv1::kDIGITS :
284 loader->LoadDigits() ;
285 data = loader->TreeD() ;
287 AliWarning(Form(" Digit Tree not found for %s", AliQAv1::GetDetName(iDet))) ;
290 qadm->Exec(taskIndex, data) ;
293 case AliQAv1::kDIGITSR :
295 loader->LoadDigits() ;
296 data = loader->TreeD() ;
298 AliWarning(Form(" Digit Tree not found for %s", AliQAv1::GetDetName(iDet))) ;
301 qadm->Exec(taskIndex, data) ;
304 case AliQAv1::kRECPOINTS :
306 loader->LoadRecPoints() ;
307 data = loader->TreeR() ;
309 AliWarning(Form("RecPoints not found for %s", AliQAv1::GetDetName(iDet))) ;
312 qadm->Exec(taskIndex, data) ;
315 case AliQAv1::kTRACKSEGMENTS :
317 case AliQAv1::kRECPARTICLES :
319 case AliQAv1::kESDS :
320 qadm->Exec(taskIndex, fESD) ;
322 case AliQAv1::kNTASKINDEX :
327 Increment(taskIndex) ;
329 // Save QA data for all detectors
333 if ( taskIndex == AliQAv1::kRAWS )
334 fRawReader->RewindEvents() ;
339 //_____________________________________________________________________________
340 Bool_t AliQAManager::Finish(const AliQAv1::TASKINDEX_t taskIndex)
342 // write output to file for all detectors
344 AliQAChecker::Instance()->SetRunNumber(fRunNumber) ;
346 for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
347 if (IsSelected(AliQAv1::GetDetName(iDet))) {
348 AliQADataMaker * qadm = GetQADataMaker(iDet) ;
350 qadm->EndOfCycle(taskIndex) ;
356 //_____________________________________________________________________________
357 TObjArray * AliQAManager::GetFromOCDB(AliQAv1::DETECTORINDEX_t det, AliQAv1::TASKINDEX_t task, const Char_t * year) const
359 // Retrieve the list of QA data for a given detector and a given task
360 TObjArray * rv = NULL ;
361 if ( !strlen(AliQAv1::GetQARefStorage()) ) {
362 AliError("No storage defined, use AliQAv1::SetQARefStorage") ;
365 if ( ! IsDefaultStorageSet() ) {
366 TString tmp(AliQAv1::GetQARefDefaultStorage()) ;
369 Instance()->SetDefaultStorage(tmp.Data()) ;
370 Instance()->SetSpecificStorage(Form("%s/*", AliQAv1::GetQAName()), AliQAv1::GetQARefStorage()) ;
372 TString detOCDBDir(Form("%s/%s/%s", AliQAv1::GetQAName(), AliQAv1::GetDetName((Int_t)det), AliQAv1::GetRefOCDBDirName())) ;
373 AliDebug(AliQAv1::GetQADebugLevel(), Form("Retrieving reference data from %s/%s for %s", AliQAv1::GetQARefStorage(), detOCDBDir.Data(), AliQAv1::GetTaskName(task).Data())) ;
374 AliCDBEntry* entry = QAManager()->Get(detOCDBDir.Data(), 0) ; //FIXME 0 --> Run Number
375 TList * listDetQAD = static_cast<TList *>(entry->GetObject()) ;
377 rv = static_cast<TObjArray *>(listDetQAD->FindObject(AliQAv1::GetTaskName(task))) ;
381 //_____________________________________________________________________________
382 TCanvas ** AliQAManager::GetImage(Char_t * detName)
384 // retrieves QA Image for the given detector
385 TCanvas ** rv = NULL ;
386 Int_t detIndex = AliQAv1::GetDetIndex(detName) ;
387 if ( detIndex != AliQAv1::kNULLDET) {
388 AliQACheckerBase * qac = AliQAChecker::Instance()->GetDetQAChecker(detIndex) ;
389 rv = qac->GetImage() ;
394 //_____________________________________________________________________________
395 AliLoader * AliQAManager::GetLoader(Int_t iDet)
397 // get the loader for a detector
399 if ( !fRunLoader || iDet == AliQAv1::kCORR || iDet == AliQAv1::kGLOBAL )
402 TString detName = AliQAv1::GetDetName(iDet) ;
403 fLoader[iDet] = fRunLoader->GetLoader(detName + "Loader");
405 return fLoader[iDet] ;
407 // load the QA data maker object
408 TPluginManager* pluginManager = gROOT->GetPluginManager() ;
409 TString loaderName = "Ali" + detName + "Loader" ;
411 AliLoader * loader = NULL ;
412 // first check if a plugin is defined for the quality assurance data maker
413 TPluginHandler* pluginHandler = pluginManager->FindHandler("AliLoader", detName) ;
414 // if not, add a plugin for it
415 if (!pluginHandler) {
416 AliDebug(AliQAv1::GetQADebugLevel(), Form("defining plugin for %s", loaderName.Data())) ;
417 TString libs = gSystem->GetLibraries() ;
418 if (libs.Contains("lib" + detName + "base.so") || (gSystem->Load("lib" + detName + "base.so") >= 0)) {
419 pluginManager->AddHandler("AliQADataMaker", detName, loaderName, detName + "loader", loaderName + "()") ;
421 pluginManager->AddHandler("AliLoader", detName, loaderName, detName, loaderName + "()") ;
423 pluginHandler = pluginManager->FindHandler("AliLoader", detName) ;
425 if (pluginHandler && (pluginHandler->LoadPlugin() == 0)) {
426 loader = (AliLoader *) pluginHandler->ExecPlugin(0) ;
429 fLoader[iDet] = loader ;
433 //_____________________________________________________________________________
434 AliQAv1 * AliQAManager::GetQA(UInt_t run, UInt_t evt)
436 // retrieves the QA object stored in a file named "Run{run}.Event{evt}_1.ESD.tag.root"
437 Char_t * fileName = Form("Run%d.Event%d_1.ESD.tag.root", run, evt) ;
438 TFile * tagFile = TFile::Open(fileName) ;
440 AliError(Form("File %s not found", fileName)) ;
443 TTree * tagTree = static_cast<TTree *>(tagFile->Get("T")) ;
445 AliError(Form("Tree T not found in %s", fileName)) ;
449 AliRunTag * tag = new AliRunTag ;
450 tagTree->SetBranchAddress("AliTAG", &tag) ;
451 tagTree->GetEntry(evt) ;
452 AliQAv1 * qa = AliQAv1::Instance(tag->GetQALength(), tag->GetQAArray(), tag->GetESLength(), tag->GetEventSpecies()) ;
457 //_____________________________________________________________________________
458 AliQADataMaker * AliQAManager::GetQADataMaker(const Int_t iDet)
460 // get the quality assurance data maker for a detector
462 AliQADataMaker * qadm = fQADataMaker[iDet] ;
466 qadm->SetEventSpecie(fEventSpecie) ;
467 if ( qadm->GetRecoParam() )
468 if ( AliRecoParam::Convert(qadm->GetRecoParam()->GetEventSpecie()) != AliRecoParam::kDefault)
469 qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ;
471 } else if (iDet == AliQAv1::kGLOBAL && strcmp(GetMode(), AliQAv1::GetModeName(AliQAv1::kRECMODE)) == 0) { //Global QA
473 qadm = new AliGlobalQADataMaker();
474 qadm->SetName(AliQAv1::GetDetName(iDet));
475 qadm->SetUniqueID(iDet);
476 fQADataMaker[iDet] = qadm;
477 qadm->SetEventSpecie(fEventSpecie) ;
478 if ( qadm->GetRecoParam() )
479 if ( AliRecoParam::Convert(qadm->GetRecoParam()->GetEventSpecie()) != AliRecoParam::kDefault)
480 qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ;
482 } else if (iDet == AliQAv1::kCORR && strcmp(GetMode(), AliQAv1::GetModeName(AliQAv1::kRECMODE)) == 0 ) { //the data maker for correlations among detectors
483 qadm = new AliCorrQADataMakerRec(fQADataMaker) ;
484 qadm->SetName(AliQAv1::GetDetName(iDet));
485 qadm->SetUniqueID(iDet);
486 fQADataMaker[iDet] = qadm;
487 qadm->SetEventSpecie(fEventSpecie) ;
488 if ( qadm->GetRecoParam() )
489 if ( AliRecoParam::Convert(qadm->GetRecoParam()->GetEventSpecie()) != AliRecoParam::kDefault)
490 qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ;
492 } else if ( iDet < AliQAv1::kGLOBAL ) {
493 TString smode(GetMode()) ;
494 if (smode.Contains(AliQAv1::GetModeName(AliQAv1::kQAMODE)))
495 smode = AliQAv1::GetModeName(AliQAv1::kRECMODE) ;
496 // load the QA data maker object
497 TPluginManager* pluginManager = gROOT->GetPluginManager() ;
498 TString detName = AliQAv1::GetDetName(iDet) ;
499 TString qadmName = "Ali" + detName + "QADataMaker" + smode ;
501 // first check if a plugin is defined for the quality assurance data maker
502 TPluginHandler* pluginHandler = pluginManager->FindHandler("AliQADataMaker", detName) ;
503 // if not, add a plugin for it
504 if (!pluginHandler) {
505 AliDebug(AliQAv1::GetQADebugLevel(), Form("defining plugin for %s", qadmName.Data())) ;
506 TString libs = gSystem->GetLibraries() ;
507 TString temp(smode) ;
509 if (libs.Contains("lib" + detName + smode + ".so") || (gSystem->Load("lib" + detName + temp.Data() + ".so") >= 0)) {
510 if ( iDet == AliQAv1::kMUON ) {
511 pluginManager->AddHandler("AliQADataMaker", detName, qadmName, detName + "qadm", qadmName + "(Bool_t,Bool_t)");
513 pluginManager->AddHandler("AliQADataMaker", detName, qadmName, detName + "qadm", qadmName + "()");
516 pluginManager->AddHandler("AliQADataMaker", detName, qadmName, detName, qadmName + "()");
518 pluginHandler = pluginManager->FindHandler("AliQADataMaker", detName) ;
520 if (pluginHandler && (pluginHandler->LoadPlugin() == 0)) {
521 if ( iDet == AliQAv1::kMUON ) {
522 Bool_t mch = fActiveOnlineDetectors.Contains("MUONTRK");
523 Bool_t mtr = fActiveOnlineDetectors.Contains("MUONTRG");
524 qadm = (AliQADataMaker *) pluginHandler->ExecPlugin(2,mch,mtr);
527 qadm = (AliQADataMaker *) pluginHandler->ExecPlugin(0);
531 qadm->SetName(AliQAv1::GetDetName(iDet));
532 qadm->SetUniqueID(iDet);
533 fQADataMaker[iDet] = qadm ;
534 qadm->SetEventSpecie(fEventSpecie) ;
535 if ( qadm->GetRecoParam() )
536 if ( AliRecoParam::Convert(qadm->GetRecoParam()->GetEventSpecie()) != AliRecoParam::kDefault)
537 qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ;
543 //_____________________________________________________________________________
544 void AliQAManager::EndOfCycle(TObjArray * detArray)
546 // End of cycle QADataMakers
548 AliQAChecker::Instance()->SetRunNumber(fRunNumber) ;
551 fakeCanvas.Print(Form("%s%s%d.%s[", AliQAv1::GetImageFileName(), GetMode(), fRunNumber, AliQAv1::GetImageFileFormat()), "ps") ;
552 for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
553 if (IsSelected(AliQAv1::GetDetName(iDet))) {
554 AliQADataMaker * qadm = GetQADataMaker(iDet) ;
557 // skip non active detectors
559 AliModule* det = static_cast<AliModule*>(detArray->FindObject(AliQAv1::GetDetName(iDet))) ;
560 if (!det || !det->IsActive())
563 AliQACheckerBase * qac = AliQAChecker::Instance()->GetDetQAChecker(iDet) ;
565 qac->SetPrintImage(fPrintImage) ;
566 for (UInt_t taskIndex = 0; taskIndex < AliQAv1::kNTASKINDEX; taskIndex++) {
567 if ( fTasks.Contains(Form("%d", taskIndex)) )
568 qadm->EndOfCycle(AliQAv1::GetTaskIndex(AliQAv1::GetTaskName(taskIndex))) ;
574 fakeCanvas.Print(Form("%s%s%d.%s]", AliQAv1::GetImageFileName(), GetMode(), fRunNumber, AliQAv1::GetImageFileFormat()), "ps");
577 //_____________________________________________________________________________
578 void AliQAManager::EndOfCycle(TString detectors)
580 // End of cycle QADataMakers
582 AliQAChecker::Instance()->SetRunNumber(fRunNumber) ;
585 fakeCanvas.Print(Form("%s%s%d.%s[", AliQAv1::GetImageFileName(), GetMode(), fRunNumber, AliQAv1::GetImageFileFormat()), "ps") ;
586 for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
587 if (IsSelected(AliQAv1::GetDetName(iDet))) {
588 AliQADataMaker * qadm = GetQADataMaker(iDet) ;
591 // skip non active detectors
592 if (!detectors.Contains(AliQAv1::GetDetName(iDet)))
594 AliQACheckerBase * qac = AliQAChecker::Instance()->GetDetQAChecker(iDet) ;
596 qac->SetPrintImage(fPrintImage) ;
597 for (UInt_t taskIndex = 0; taskIndex < AliQAv1::kNTASKINDEX; taskIndex++) {
598 if ( fTasks.Contains(Form("%d", taskIndex)) )
599 qadm->EndOfCycle(AliQAv1::GetTaskIndex(AliQAv1::GetTaskName(taskIndex))) ;
605 fakeCanvas.Print(Form("%s%s%d.%s]", AliQAv1::GetImageFileName(), GetMode(), fRunNumber, AliQAv1::GetImageFileFormat()), "ps");
608 //_____________________________________________________________________________
609 AliRecoParam::EventSpecie_t AliQAManager::GetEventSpecieFromESD()
611 AliRecoParam::EventSpecie_t runtype = AliRecoParam::kDefault ;
612 if (!gSystem->AccessPathName("AliESDs.root")) { // AliESDs.root exists
613 TFile * esdFile = TFile::Open("AliESDs.root") ;
614 TTree * esdTree = static_cast<TTree *> (esdFile->Get("esdTree")) ;
616 AliError("esdTree not found") ;
618 AliESDEvent * esd = new AliESDEvent() ;
619 esd->ReadFromTree(esdTree) ;
620 esdTree->GetEntry(0) ;
621 runtype = AliRecoParam::Convert(esd->GetEventType()) ;
624 AliError("AliESDs.root not found") ;
629 //_____________________________________________________________________________
630 void AliQAManager::Increment(const AliQAv1::TASKINDEX_t taskIndex)
632 // Increments the cycle counter for all QA Data Makers
633 static AliQAv1::TASKINDEX_t currentTask = AliQAv1::kNTASKINDEX ;
634 if ( (currentTask == taskIndex) && taskIndex != AliQAv1::kNULLTASKINDEX )
637 currentTask = taskIndex ;
638 for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
639 if (IsSelected(AliQAv1::GetDetName(iDet))) {
640 AliQADataMaker * qadm = GetQADataMaker(iDet) ;
647 //_____________________________________________________________________________
648 Bool_t AliQAManager::InitQA(const AliQAv1::TASKINDEX_t taskIndex, const Char_t * input )
650 // Initialize the event source and QA data makers
652 fTasks += Form("%d", taskIndex) ;
654 if (taskIndex == AliQAv1::kRAWS) {
656 fRawReader = AliRawReader::Create(input);
660 fRawReaderDelete = kTRUE ;
661 fRawReader->NextEvent() ;
662 fRunNumber = fRawReader->GetRunNumber() ;
664 fRawReader->RewindEvents();
665 fNumberOfEvents = 999999 ;
666 if ( fMaxEvents < 0 )
667 fMaxEvents = fNumberOfEvents ;
668 } else if (taskIndex == AliQAv1::kESDS) {
669 fTasks = AliQAv1::GetTaskName(AliQAv1::kESDS) ;
670 if (!gSystem->AccessPathName("AliESDs.root")) { // AliESDs.root exists
671 TFile * esdFile = TFile::Open("AliESDs.root") ;
672 fESDTree = static_cast<TTree *> (esdFile->Get("esdTree")) ;
674 AliError("esdTree not found") ;
677 fESD = new AliESDEvent() ;
678 fESD->ReadFromTree(fESDTree) ;
679 fESDTree->GetEntry(0) ;
680 fRunNumber = fESD->GetRunNumber() ;
681 fNumberOfEvents = fESDTree->GetEntries() ;
682 if ( fMaxEvents < 0 )
683 fMaxEvents = fNumberOfEvents ;
686 AliError("AliESDs.root not found") ;
690 if ( !InitRunLoader() ) {
691 AliWarning("No Run Loader not found") ;
693 fNumberOfEvents = fRunLoader->GetNumberOfEvents() ;
694 if ( fMaxEvents < 0 )
695 fMaxEvents = fNumberOfEvents ;
700 TObjArray* detArray = NULL ;
701 if (fRunLoader) // check if RunLoader exists
702 if ( fRunLoader->GetAliRun() ) { // check if AliRun exists in gAlice.root
703 detArray = fRunLoader->GetAliRun()->Detectors() ;
704 fRunNumber = fRunLoader->GetHeader()->GetRun() ;
707 // Initialize all QA data makers for all detectors
708 fRunNumber = AliCDBManager::Instance()->GetRun() ;
709 if ( ! AliGeomManager::GetGeometry() )
710 AliGeomManager::LoadGeometry() ;
712 InitQADataMaker(fRunNumber, detArray) ; //, fCycleSame, kTRUE, detArray) ;
717 while (timer.CpuTime()<5) {
719 gSystem->ProcessEvents();
721 fakeCanvas.Print(Form("%s%s%d.%s[", AliQAv1::GetImageFileName(), GetMode(), fRunNumber, AliQAv1::GetImageFileFormat()), "ps") ;
726 //_____________________________________________________________________________
727 void AliQAManager::InitQADataMaker(UInt_t run, TObjArray * detArray)
729 // Initializes The QADataMaker for all active detectors and for all active tasks
731 for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
732 if (IsSelected(AliQAv1::GetDetName(iDet))) {
733 AliQADataMaker * qadm = GetQADataMaker(iDet) ;
735 AliError(Form("AliQADataMaker not found for %s", AliQAv1::GetDetName(iDet))) ;
736 fDetectorsW.ReplaceAll(AliQAv1::GetDetName(iDet), "") ;
738 if (fQAWriteExpert[iDet])
739 qadm->SetWriteExpert() ;
740 AliDebug(AliQAv1::GetQADebugLevel(), Form("Data Maker found for %s %d", qadm->GetName(), qadm->WriteExpert())) ;
741 // skip non active detectors
743 AliModule* det = static_cast<AliModule*>(detArray->FindObject(AliQAv1::GetDetName(iDet))) ;
744 if (!det || !det->IsActive())
747 // Set default reco params
748 Bool_t sameCycle = kFALSE ;
749 for (UInt_t taskIndex = 0; taskIndex < AliQAv1::kNTASKINDEX; taskIndex++) {
750 if ( fTasks.Contains(Form("%d", taskIndex)) ) {
751 qadm->Init(AliQAv1::GetTaskIndex(AliQAv1::GetTaskName(taskIndex)), GetQACycles(qadm->GetUniqueID())) ;
752 qadm->StartOfCycle(AliQAv1::GetTaskIndex(AliQAv1::GetTaskName(taskIndex)), run, sameCycle) ;
762 //_____________________________________________________________________________
763 Bool_t AliQAManager::InitRunLoader()
765 // get or create the run loader
769 if (!gSystem->AccessPathName(fGAliceFileName.Data())) { // galice.root exists
770 // load all base libraries to get the loader classes
771 TString libs = gSystem->GetLibraries() ;
772 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
773 if (!IsSelected(AliQAv1::GetDetName(iDet)))
775 TString detName = AliQAv1::GetDetName(iDet) ;
776 if (detName == "HLT")
778 if (libs.Contains("lib" + detName + "base.so"))
780 gSystem->Load("lib" + detName + "base.so");
782 fRunLoader = AliRunLoader::Open(fGAliceFileName.Data());
784 AliError(Form("no run loader found in file %s", fGAliceFileName.Data()));
787 fRunLoader->CdGAFile();
788 if (fRunLoader->LoadgAlice() == 0) {
789 gAlice = fRunLoader->GetAliRun();
793 AliError(Form("no gAlice object found in file %s", fGAliceFileName.Data()));
797 } else { // galice.root does not exist
798 AliError(Form("the file %s does not exist", fGAliceFileName.Data()));
804 fRunLoader->LoadHeader();
805 fRunNumber = fRunLoader->GetHeader()->GetRun() ;
810 //_____________________________________________________________________________
811 Bool_t AliQAManager::IsSelected(const Char_t * det)
813 // check whether detName is contained in detectors
814 // if yes, it is removed from detectors
817 const TString detName(det) ;
818 // always activates Correlation
819 // if ( detName.Contains(AliQAv1::GetDetName(AliQAv1::kCORR)) || detName.Contains(AliQAv1::GetDetName(AliQAv1::kGLOBAL))) {
822 // check if all detectors are selected
823 if (fDetectors.Contains("ALL")) {
826 } else if ((fDetectors.CompareTo(detName) == 0) ||
827 fDetectors.BeginsWith(detName+" ") ||
828 fDetectors.EndsWith(" "+detName) ||
829 fDetectors.Contains(" "+detName+" ")) {
836 //_____________________________________________________________________________
837 Bool_t AliQAManager::Merge(Int_t runNumber, const char *fileName) const
839 // Merge data from all detectors from a given run in one single file
840 // Merge the QA results from all the data chunks in one run
841 // The 'fileName' is name of the output file with merged QA data
842 if ( runNumber == -1)
843 runNumber = fRunNumber ;
844 Bool_t rv = MergeData(runNumber,fileName) ;
845 //rv *= MergeResults(runNumber) ; // not needed for the time being
849 //______________________________________________________________________
850 Bool_t AliQAManager::MergeXML(const Char_t * collectionFile, const Char_t * subFile, const Char_t * outFile)
852 // merges files listed in a xml collection
853 // usage Merge(collection, outputFile))
854 // collection: is a xml collection
858 if ( strstr(collectionFile, ".xml") == 0 ) {
859 AliError("Input collection file must be an \".xml\" file\n") ;
864 TGrid::Connect("alien://");
868 // Open the file collection
869 AliInfoClass(Form("*** Create Collection ***\n*** Wk-Dir = |%s| \n*** Coll = |%s| \n",gSystem->WorkingDirectory(), collectionFile));
871 TGridCollection * collection = (TGridCollection*)gROOT->ProcessLine(Form("TAlienCollection::Open(\"%s\")",collectionFile));
872 TGridResult* result = collection->GetGridResult("", 0, 0);
875 const Char_t * turl ;
876 TFileMerger merger(kFALSE) ;
878 TString tempo(collectionFile) ;
880 tempo.ReplaceAll(".xml", subFile) ;
882 tempo.ReplaceAll(".xml", "_Merged.root") ;
883 outFile = tempo.Data() ;
885 merger.OutputFile(outFile) ;
887 while ( (turl = result->GetKey(index, "turl")) ) {
890 file = Form("%s#%s", turl, subFile) ;
892 file = Form("%s", turl) ;
894 AliDebug(AliQAv1::GetQADebugLevel(), Form("%s\n", file)) ;
895 merger.AddFile(file) ;
902 AliDebug(AliQAv1::GetQADebugLevel(), Form("Files merged into %s\n", outFile)) ;
908 //_____________________________________________________________________________
909 void AliQAManager::MergeCustom() const
911 // Custom Merge of QA data from all detectors for all runs in one single file
912 // search all the run numbers
913 // search all the run numbers
914 gROOT->ProcessLine(".! ls *QA*.root > QAtempo.txt") ;
916 FILE * theQAfiles = fopen("QAtempo.txt", "r") ;
919 TIter nextRun(&srunList) ;
920 TObjString * srun = NULL ;
921 Int_t loRun = 999999999 ;
923 while ( theQAfile.Gets(theQAfiles) ) {
924 Bool_t runExist = kFALSE ;
925 TString srunNew(theQAfile(theQAfile.Index("QA.")+3, theQAfile.Index(".root")-(theQAfile.Index("QA.")+3))) ;
926 Int_t cuRun = srunNew.Atoi() ;
931 while ( (srun = static_cast<TObjString *> (nextRun())) ) {
932 if ( cuRun == (srun->String()).Atoi() ) {
939 srunList.Add(new TObjString(srunNew.Data()));
942 Int_t runNumber = 0 ;
943 TFile mergedFile(Form("Merged.%s.Data.root", AliQAv1::GetQADataFileName()), "RECREATE") ;
944 TH1I * hisRun = new TH1I("hLMR", "List of merged runs", hiRun-loRun+10, loRun, hiRun+10) ;
945 // create the structure into the merged file
946 for (Int_t iDet = 0; iDet < AliQAv1::kNDET ; iDet++) {
947 TDirectory * detDir = mergedFile.mkdir(AliQAv1::GetDetName(iDet)) ;
948 for (Int_t taskIndex = 0; taskIndex < AliQAv1::kNTASKINDEX; taskIndex++) {
950 TDirectory * taskDir = gDirectory->mkdir(AliQAv1::GetTaskName(taskIndex)) ;
951 for (Int_t es = 0 ; es < AliRecoParam::kNSpecies ; es++) {
953 TDirectory * esDir = gDirectory->mkdir(AliRecoParam::GetEventSpecieName(es)) ;
955 gDirectory->mkdir(AliQAv1::GetExpert()) ;
959 while ( (srun = static_cast<TObjString *> (nextRun())) ) {
960 runNumber = (srun->String()).Atoi() ;
961 hisRun->Fill(runNumber) ;
962 AliDebug(AliQAv1::GetQADebugLevel(), Form("Merging run number %d", runNumber)) ;
963 // search all QA files for runNumber in the current directory
964 Char_t * fileList[AliQAv1::kNDET] ;
966 for (Int_t iDet = 0; iDet < AliQAv1::kNDET ; iDet++) {
967 Char_t * file = gSystem->Which(gSystem->WorkingDirectory(), Form("%s.%s.%d.root", AliQAv1::GetDetName(iDet), AliQAv1::GetQADataFileName(), runNumber));
969 fileList[index++] = file ;
972 AliError("No QA data file found\n") ;
975 for ( Int_t i = 0 ; i < index ; i++) {
976 TFile * inFile = TFile::Open(fileList[i]) ;
977 TList * listOfKeys =inFile->GetListOfKeys() ;
978 TIter nextkey(listOfKeys) ;
980 TString dirName("") ;
981 while ( (obj1 = nextkey()) ) {
982 TDirectory * directoryDet = inFile->GetDirectory(obj1->GetName()) ;
983 if ( directoryDet ) {
984 AliDebug(AliQAv1::GetQADebugLevel(), Form("%s dir = %s", inFile->GetName(), directoryDet->GetName())) ;
985 dirName += Form("%s/", directoryDet->GetName() ) ;
987 TList * listOfTasks = directoryDet->GetListOfKeys() ;
988 TIter nextTask(listOfTasks) ;
990 while ( (obj2 = nextTask()) ) {
991 TDirectory * directoryTask = directoryDet->GetDirectory(obj2->GetName()) ;
992 if ( directoryTask ) {
993 dirName += Form("%s", obj2->GetName()) ;
994 AliDebug(AliQAv1::GetQADebugLevel(), Form("%s", dirName.Data())) ;
995 directoryTask->cd() ;
996 TList * listOfEventSpecie = directoryTask->GetListOfKeys() ;
997 TIter nextEventSpecie(listOfEventSpecie) ;
999 while ( (obj3 = nextEventSpecie()) ) {
1000 TDirectory * directoryEventSpecie = directoryTask->GetDirectory(obj3->GetName()) ;
1001 if ( directoryEventSpecie ) {
1002 dirName += Form("/%s/", obj3->GetName()) ;
1003 AliDebug(AliQAv1::GetQADebugLevel(), Form("%s\n", dirName.Data())) ;
1004 directoryEventSpecie->cd() ;
1005 // histograms are here
1006 TDirectory * mergedDirectory = mergedFile.GetDirectory(dirName.Data()) ;
1007 TList * listOfData = directoryEventSpecie->GetListOfKeys() ;
1008 TIter nextData(listOfData) ;
1010 while ( (key = static_cast<TKey *>(nextData())) ) {
1011 TString className(key->GetClassName()) ;
1012 if ( className.Contains("TH") || className.Contains("TProfile") ) {
1013 TH1 * histIn = static_cast<TH1*> (key->ReadObj()) ;
1014 TH1 * histOu = static_cast<TH1*> (mergedDirectory->FindObjectAny(histIn->GetName())) ;
1015 AliDebug(AliQAv1::GetQADebugLevel(), Form("%s %p %p\n", key->GetName(), histIn, histOu)) ;
1016 mergedDirectory->cd() ;
1020 histOu->Add(histIn) ;
1021 histOu->Write(histOu->GetName(), kOverwrite) ;
1024 else if ( className.Contains("TDirectoryFile") ) {
1025 TDirectory * dirExpert = directoryEventSpecie->GetDirectory(key->GetName()) ;
1027 TDirectory * mergedDirectoryExpert = mergedDirectory->GetDirectory(dirExpert->GetName()) ;
1028 TList * listOfExpertData = dirExpert->GetListOfKeys() ;
1029 TIter nextExpertData(listOfExpertData) ;
1031 while ( (keykey = static_cast<TKey *>(nextExpertData())) ) {
1032 TString classNameExpert(keykey->GetClassName()) ;
1033 if (classNameExpert.Contains("TH")) {
1034 TH1 * histInExpert = static_cast<TH1*> (keykey->ReadObj()) ;
1035 TH1 * histOuExpert = static_cast<TH1*> (mergedDirectory->FindObjectAny(histInExpert->GetName())) ;
1036 mergedDirectoryExpert->cd() ;
1037 if ( ! histOuExpert ) {
1038 histInExpert->Write() ;
1040 histOuExpert->Add(histInExpert) ;
1041 histOuExpert->Write(histOuExpert->GetName(), kOverwrite) ;
1046 AliError(Form("No merge done for this object %s in %s", key->GetName(), dirName.Data())) ;
1049 dirName.ReplaceAll(Form("/%s/",obj3->GetName()), "") ;
1052 dirName.ReplaceAll(obj2->GetName(), "") ;
1062 mergedFile.Close() ;
1066 //_____________________________________________________________________________
1067 Bool_t AliQAManager::MergeData(const Int_t runNumber, const char *fileName) const
1069 // Merge QA data from all detectors for a given run in one single file
1071 TFileMerger merger(kFALSE) ;
1072 TString outFileName = fileName;
1073 if (outFileName.IsNull()) outFileName.Form("Merged.%s.Data.root",AliQAv1::GetQADataFileName());
1074 merger.OutputFile(outFileName.Data()) ;
1075 for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
1076 Char_t * file = gSystem->Which(gSystem->WorkingDirectory(), Form("%s.%s.%d.root", AliQAv1::GetDetName(iDet), AliQAv1::GetQADataFileName(), runNumber));
1078 merger.AddFile(file);
1085 //_____________________________________________________________________________
1086 Bool_t AliQAManager::MergeResults(const Int_t runNumber) const
1088 // Merge the QA result from all the data chunks in a run
1089 // to be revised whwn it will be used (see MergeData)
1091 cmd = Form(".! ls %s*.root > tempo.txt", AliQAv1::GetQADataFileName()) ;
1092 gROOT->ProcessLine(cmd.Data()) ;
1093 ifstream in("tempo.txt") ;
1094 const Int_t chunkMax = 100 ;
1095 TString fileList[chunkMax] ;
1100 in >> fileList[index] ;
1103 AliDebug(AliQAv1::GetQADebugLevel(), Form("index = %d file = %s", index, (fileList[index].Data()))) ;
1108 AliError("No QA Result File found") ;
1112 TFileMerger merger ;
1113 TString outFileName ;
1114 if (runNumber != -1)
1115 outFileName = Form("Merged.%s.Result.%d.root",AliQAv1::GetQADataFileName(),runNumber);
1117 outFileName = Form("Merged.%s.Result.root",AliQAv1::GetQADataFileName());
1118 merger.OutputFile(outFileName.Data()) ;
1119 for (Int_t ifile = 0 ; ifile < index ; ifile++) {
1120 TString file = fileList[ifile] ;
1121 merger.AddFile(file) ;
1128 //_____________________________________________________________________________
1129 void AliQAManager::Reset(const Bool_t sameCycle)
1131 // Reset the default data members
1133 for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
1134 if (IsSelected(AliQAv1::GetDetName(iDet))) {
1135 AliQADataMaker * qadm = GetQADataMaker(iDet);
1140 if (fRawReaderDelete) {
1145 fCycleSame = sameCycle ;
1149 fNumberOfEvents = 999999 ;
1152 //_____________________________________________________________________________
1153 void AliQAManager::ResetDetectors(AliQAv1::TASKINDEX_t task, AliQAv1::DETECTORINDEX_t det)
1155 //calls ResetDetector of specified or all detectors
1157 UInt_t iDetMax = fgkNDetectors ;
1158 if ( det != AliQAv1::kNULLDET ) {
1163 for (iDet = 0; iDet < iDetMax ; iDet++) {
1164 if (IsSelected(AliQAv1::GetDetName(iDet))) {
1165 AliQADataMaker * qadm = GetQADataMaker(iDet);
1166 qadm->ResetDetector(task);
1171 //_____________________________________________________________________________
1172 AliQAManager * AliQAManager::QAManager(AliQAv1::MODE_t mode, TMap *entryCache, Int_t run)
1174 // returns AliQAManager instance (singleton)
1176 if (!fgQAInstance) {
1177 if ( (mode != AliQAv1::kSIMMODE) && (mode != AliQAv1::kRECMODE) && (mode != AliQAv1::kQAMODE) ) {
1178 AliWarningClass("You must specify kSIMMODE or kRECMODE or kQAMODE") ;
1181 fgQAInstance = new AliQAManager(mode) ;
1183 fgQAInstance->Init();
1185 fgQAInstance->InitFromCache(entryCache,run);
1187 return fgQAInstance;
1190 //_____________________________________________________________________________
1191 AliQAManager * AliQAManager::QAManager(AliQAv1::TASKINDEX_t task)
1193 // returns AliQAManager instance (singleton)
1194 return QAManager(AliQAv1::Mode(task)) ;
1197 //_____________________________________________________________________________
1198 TString AliQAManager::Run(const Char_t * detectors, AliRawReader * rawReader, const Bool_t sameCycle)
1200 //Runs all the QA data Maker for Raws only
1202 fCycleSame = sameCycle ;
1203 fRawReader = rawReader ;
1204 fDetectors = detectors ;
1205 fDetectorsW = detectors ;
1207 AliCDBManager* man = AliCDBManager::Instance() ;
1209 if ( man->GetRun() == -1 ) {// check if run number not set previously and set it from raw data
1210 rawReader->NextEvent() ;
1211 man->SetRun(fRawReader->GetRunNumber()) ;
1212 rawReader->RewindEvents() ;
1216 if ( !InitQA(AliQAv1::kRAWS) )
1218 fRawReaderDelete = kFALSE ;
1220 DoIt(AliQAv1::kRAWS) ;
1221 return fDetectorsW ;
1224 //_____________________________________________________________________________
1225 TString AliQAManager::Run(const Char_t * detectors, const Char_t * fileName, const Bool_t sameCycle)
1227 //Runs all the QA data Maker for Raws only
1229 fCycleSame = sameCycle ;
1230 fDetectors = detectors ;
1231 fDetectorsW = detectors ;
1233 AliCDBManager* man = AliCDBManager::Instance() ;
1234 if ( man->GetRun() == -1 ) { // check if run number not set previously and set it from AliRun
1235 AliRunLoader * rl = AliRunLoader::Open("galice.root") ;
1237 AliFatal("galice.root file not found in current directory") ;
1241 if ( ! rl->GetAliRun() ) {
1242 AliFatal("AliRun not found in galice.root") ;
1245 man->SetRun(rl->GetHeader()->GetRun());
1251 if ( !InitQA(AliQAv1::kRAWS, fileName) )
1254 DoIt(AliQAv1::kRAWS) ;
1255 return fDetectorsW ;
1258 //_____________________________________________________________________________
1259 TString AliQAManager::Run(const Char_t * detectors, const AliQAv1::TASKINDEX_t taskIndex, Bool_t const sameCycle, const Char_t * fileName )
1261 // Runs all the QA data Maker for every detector
1263 fCycleSame = sameCycle ;
1264 fDetectors = detectors ;
1265 fDetectorsW = detectors ;
1267 AliCDBManager* man = AliCDBManager::Instance() ;
1268 if ( man->GetRun() == -1 ) { // check if run number not set previously and set it from AliRun
1269 AliRunLoader * rl = AliRunLoader::Open("galice.root") ;
1271 AliFatal("galice.root file not found in current directory") ;
1275 if ( ! rl->GetAliRun() ) {
1276 AliDebug(AliQAv1::GetQADebugLevel(), "AliRun not found in galice.root") ;
1279 man->SetRun(rl->GetHeader()->GetRun()) ;
1283 if ( taskIndex == AliQAv1::kNULLTASKINDEX) {
1284 for (UInt_t task = 0; task < AliQAv1::kNTASKINDEX; task++) {
1285 if ( fTasks.Contains(Form("%d", task)) ) {
1287 if ( !InitQA(AliQAv1::GetTaskIndex(AliQAv1::GetTaskName(task)), fileName) )
1289 DoIt(AliQAv1::GetTaskIndex(AliQAv1::GetTaskName(task))) ;
1294 if ( !InitQA(taskIndex, fileName) )
1298 return fDetectorsW ;
1301 //_____________________________________________________________________________
1302 void AliQAManager::RunOneEvent(AliRawReader * rawReader)
1304 //Runs all the QA data Maker for Raws only and on one event only (event loop done by calling method)
1308 if (fTasks.Contains(Form("%d", AliQAv1::kRAWS))){
1310 if ( GetEventInfo())
1311 detList = GetEventInfo()->GetTriggerCluster() ;
1312 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
1313 if (!IsSelected(AliQAv1::GetDetName(iDet)) || (!detList.IsNull() && !detList.Contains(AliQAv1::GetDetName(iDet))))
1315 AliQADataMaker *qadm = GetQADataMaker(iDet);
1318 if ( qadm->IsCycleDone() ) {
1319 qadm->EndOfCycle() ;
1321 qadm->SetEventSpecie(fEventSpecie) ;
1322 if ( qadm->GetRecoParam() )
1323 if ( AliRecoParam::Convert(qadm->GetRecoParam()->GetEventSpecie()) != AliRecoParam::kDefault)
1324 qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ;
1325 qadm->Exec(AliQAv1::kRAWS, rawReader) ;
1330 //_____________________________________________________________________________
1331 void AliQAManager::RunOneEvent(AliESDEvent *& esd, AliESDEvent *& hltesd)
1333 //Runs all the QA data Maker for ESDs only and on one event only (event loop done by calling method)
1335 if (fTasks.Contains(Form("%d", AliQAv1::kESDS))) {
1337 if ( GetEventInfo())
1338 detList = GetEventInfo()->GetTriggerCluster() ;
1339 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
1340 if (!IsSelected(AliQAv1::GetDetName(iDet)) || (!detList.IsNull() && !detList.Contains(AliQAv1::GetDetName(iDet))))
1342 AliQADataMaker *qadm = GetQADataMaker(iDet);
1345 qadm->SetEventSpecie(fEventSpecie) ;
1346 if ( qadm->GetRecoParam() )
1347 if ( AliRecoParam::Convert(qadm->GetRecoParam()->GetEventSpecie()) != AliRecoParam::kDefault)
1348 qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ;
1349 if ( qadm->IsCycleDone() ) {
1350 qadm->EndOfCycle() ;
1352 if (iDet == AliQAv1::kHLT) {
1355 esdarray.Add(hltesd);
1356 qadm->Exec(AliQAv1::kESDS, &esdarray);
1358 qadm->Exec(AliQAv1::kESDS, esd) ;
1364 //_____________________________________________________________________________
1365 void AliQAManager::RunOneEventInOneDetector(Int_t det, TTree * tree)
1367 // Runs all the QA data Maker for ESDs only and on one event only (event loop done by calling method)
1370 if ( GetEventInfo())
1371 detList = GetEventInfo()->GetTriggerCluster() ;
1372 if (!detList.IsNull() && !detList.Contains(AliQAv1::GetDetName(det)))
1375 TString test(tree->GetName()) ;
1376 if (fTasks.Contains(Form("%d", AliQAv1::kRECPOINTS))) {
1377 if (IsSelected(AliQAv1::GetDetName(det))) {
1378 AliQADataMaker *qadm = GetQADataMaker(det);
1380 qadm->SetEventSpecie(fEventSpecie) ;
1381 if ( qadm->GetRecoParam() ) {
1382 if ( AliRecoParam::Convert(qadm->GetRecoParam()->GetEventSpecie()) != AliRecoParam::kDefault)
1383 qadm->SetEventSpecie(qadm->GetRecoParam()->GetEventSpecie()) ;
1385 AliError(Form("%d defined by %s is not an event specie", qadm->GetRecoParam()->GetEventSpecie(), qadm->GetName())) ;
1387 if ( qadm->IsCycleDone() ) {
1388 qadm->EndOfCycle() ;
1390 if (test.Contains("TreeD")) {
1391 qadm->Exec(AliQAv1::kDIGITSR, tree) ;
1392 } else if (test.Contains("TreeR")) {
1393 qadm->Exec(AliQAv1::kRECPOINTS, tree) ;
1400 //_____________________________________________________________________________
1401 Bool_t AliQAManager::Save2OCDB(const Int_t runNumber, AliRecoParam::EventSpecie_t es, const Char_t * year, const Char_t * detectors) const
1403 // take the locasl QA data merge into a single file and save in OCDB
1405 TString tmp(AliQAv1::GetQARefStorage()) ;
1406 if ( tmp.IsNull() ) {
1407 AliError("No storage defined, use AliQAv1::SetQARefStorage") ;
1410 if ( !(tmp.Contains(AliQAv1::GetLabLocalOCDB()) || tmp.Contains(AliQAv1::GetLabAliEnOCDB())) ) {
1411 AliError(Form("%s is a wrong storage, use %s or %s", AliQAv1::GetQARefStorage(), AliQAv1::GetLabLocalOCDB().Data(), AliQAv1::GetLabAliEnOCDB().Data())) ;
1414 TString sdet(detectors) ;
1417 if ( sdet.Contains("ALL") ) {
1418 rv = Merge(runNumber) ;
1421 TString inputFileName(Form("Merged.%s.Data.%d.root", AliQAv1::GetQADataFileName(), runNumber)) ;
1422 inputFile = TFile::Open(inputFileName.Data()) ;
1423 rv = SaveIt2OCDB(runNumber, inputFile, year, es) ;
1425 for (Int_t index = 0; index < AliQAv1::kNDET; index++) {
1426 if (sdet.Contains(AliQAv1::GetDetName(index))) {
1427 TString inputFileName(Form("%s.%s.%d.root", AliQAv1::GetDetName(index), AliQAv1::GetQADataFileName(), runNumber)) ;
1428 inputFile = TFile::Open(inputFileName.Data()) ;
1429 rv *= SaveIt2OCDB(runNumber, inputFile, year, es) ;
1436 //_____________________________________________________________________________
1437 Bool_t AliQAManager::SaveIt2OCDB(const Int_t runNumber, TFile * inputFile, const Char_t * year, AliRecoParam::EventSpecie_t es) const
1439 // reads the TH1 from file and adds it to appropriate list before saving to OCDB
1441 AliDebug(AliQAv1::GetQADebugLevel(), Form("Saving TH1s in %s to %s", inputFile->GetName(), AliQAv1::GetQARefStorage())) ;
1442 if ( ! IsDefaultStorageSet() ) {
1443 TString tmp( AliQAv1::GetQARefStorage() ) ;
1444 if ( tmp.Contains(AliQAv1::GetLabLocalOCDB()) )
1445 Instance()->SetDefaultStorage(AliQAv1::GetQARefStorage()) ;
1447 TString tmp1(AliQAv1::GetQARefDefaultStorage()) ;
1449 tmp1.Append("?user=alidaq") ;
1450 Instance()->SetDefaultStorage(tmp1.Data()) ;
1453 Instance()->SetSpecificStorage("*", AliQAv1::GetQARefStorage()) ;
1455 Instance()->SetRun(runNumber);
1457 AliCDBMetaData mdr ;
1458 mdr.SetResponsible("yves schutz");
1460 for ( Int_t detIndex = 0 ; detIndex < AliQAv1::kNDET ; detIndex++) {
1461 TDirectory * detDir = inputFile->GetDirectory(AliQAv1::GetDetName(detIndex)) ;
1463 AliDebug(AliQAv1::GetQADebugLevel(), Form("Entering %s", detDir->GetName())) ;
1464 AliQAv1::SetQARefDataDirName(es) ;
1465 TString detOCDBDir(Form("%s/%s/%s", AliQAv1::GetDetName(detIndex), AliQAv1::GetRefOCDBDirName(), AliQAv1::GetRefDataDirName())) ;
1466 AliCDBId idr(detOCDBDir.Data(), runNumber, AliCDBRunRange::Infinity()) ;
1467 TList * listDetQAD = new TList() ;
1468 TString listName(Form("%s QA data Reference", AliQAv1::GetDetName(detIndex))) ;
1469 mdr.SetComment(Form("%s QA stuff", AliQAv1::GetDetName(detIndex)));
1470 listDetQAD->SetName(listName) ;
1471 TList * taskList = detDir->GetListOfKeys() ;
1472 TIter nextTask(taskList) ;
1474 while ( (taskKey = static_cast<TKey*>(nextTask())) ) {
1475 TDirectory * taskDir = detDir->GetDirectory(taskKey->GetName()) ;
1476 TDirectory * esDir = taskDir->GetDirectory(AliRecoParam::GetEventSpecieName(es)) ;
1477 AliDebug(AliQAv1::GetQADebugLevel(), Form("Saving %s", esDir->GetName())) ;
1478 TObjArray * listTaskQAD = new TObjArray(100) ;
1479 listTaskQAD->SetName(Form("%s/%s", taskKey->GetName(), AliRecoParam::GetEventSpecieName(es))) ;
1480 listDetQAD->Add(listTaskQAD) ;
1481 TList * histList = esDir->GetListOfKeys() ;
1482 TIter nextHist(histList) ;
1484 while ( (histKey = static_cast<TKey*>(nextHist())) ) {
1485 TObject * odata = esDir->Get(histKey->GetName()) ;
1487 AliError(Form("%s in %s/%s returns a NULL pointer !!", histKey->GetName(), detDir->GetName(), taskDir->GetName())) ;
1489 if ( AliQAv1::GetExpert() == histKey->GetName() ) {
1490 TDirectory * expertDir = esDir->GetDirectory(histKey->GetName()) ;
1491 TList * expertHistList = expertDir->GetListOfKeys() ;
1492 TIter nextExpertHist(expertHistList) ;
1493 TKey * expertHistKey ;
1494 while ( (expertHistKey = static_cast<TKey*>(nextExpertHist())) ) {
1495 TObject * expertOdata = expertDir->Get(expertHistKey->GetName()) ;
1496 if ( !expertOdata ) {
1497 AliError(Form("%s in %s/%s/Expert returns a NULL pointer !!", expertHistKey->GetName(), detDir->GetName(), taskDir->GetName())) ;
1499 AliDebug(AliQAv1::GetQADebugLevel(), Form("Adding %s", expertHistKey->GetName())) ;
1500 if ( expertOdata->IsA()->InheritsFrom("TH1") ) {
1501 AliDebug(AliQAv1::GetQADebugLevel(), Form("Adding %s", expertHistKey->GetName())) ;
1502 TH1 * hExpertdata = static_cast<TH1*>(expertOdata) ;
1503 listTaskQAD->Add(hExpertdata) ;
1508 AliDebug(AliQAv1::GetQADebugLevel(), Form("Adding %s", histKey->GetName())) ;
1509 if ( odata->IsA()->InheritsFrom("TH1") ) {
1510 AliDebug(AliQAv1::GetQADebugLevel(), Form("Adding %s", histKey->GetName())) ;
1511 TH1 * hdata = static_cast<TH1*>(odata) ;
1512 listTaskQAD->Add(hdata) ;
1517 Instance()->Put(listDetQAD, idr, &mdr) ;
1523 //_____________________________________________________________________________
1525 void AliQAManager::SetCheckerExternParam(AliQAv1::DETECTORINDEX_t detIndex, TList * parameterList)
1527 // set the external parameters list for the detector checkers
1528 AliQACheckerBase * qac = AliQAChecker::Instance()->GetDetQAChecker(detIndex) ;
1529 qac->SetExternParamlist(parameterList) ;
1530 qac->PrintExternParam() ;
1533 //_____________________________________________________________________________
1534 void AliQAManager::SetEventSpecie(AliRecoParam::EventSpecie_t es)
1536 // set the current event specie and inform AliQAv1 that this event specie has been encountered
1538 AliQAv1::Instance()->SetEventSpecie(es) ;
1541 //_____________________________________________________________________________
1542 void AliQAManager::SetRecoParam(const Int_t det, const AliDetectorRecoParam *par)
1544 // Set custom reconstruction parameters for a given detector
1545 // Single set of parameters for all the events
1546 GetQADataMaker(det)->SetRecoParam(par) ;
1549 //_____________________________________________________________________________
1550 void AliQAManager::SetWriteExpert()
1552 // enable the writing of QA expert data
1553 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
1554 if (IsSelected(AliQAv1::GetDetName(iDet)))
1555 fQAWriteExpert[iDet] = kTRUE ;
1559 //_____________________________________________________________________________
1560 void AliQAManager::Destroy() {
1561 // delete AliQAManager instance and
1562 // all associated objects
1565 delete fgQAInstance ;
1566 fgQAInstance = NULL ;
1570 //_____________________________________________________________________________
1571 void AliQAManager::ShowQA() {
1572 // Show the result of the QA checking
1573 // for all detectors
1574 for ( Int_t detIndex = 0 ; detIndex < AliQAv1::kNDET ; detIndex++)
1575 if ( IsSelected(AliQAv1::GetDetName(detIndex)) )
1576 AliQAv1::Instance(AliQAv1::GetDetIndex(AliQAv1::GetDetName(detIndex)))->Show() ;