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 **************************************************************************/
17 ///////////////////////////////////////////////////////////////////////////////
19 // class for running the QA makers //
21 // AliQADataMakerSteer qas; //
22 // qas.Run(AliQA::kRAWS, rawROOTFileName); //
23 // qas.Run(AliQA::kHITS); //
24 // qas.Run(AliQA::kSDIGITS); //
25 // qas.Run(AliQA::kDIGITS); //
26 // qas.Run(AliQA::kRECPOINTS); //
27 // qas.Run(AliQA::kESDS); //
29 ///////////////////////////////////////////////////////////////////////////////
33 #include <TFileMerger.h>
34 #include <TPluginManager.h>
39 #include "AliCDBManager.h"
40 #include "AliCDBEntry.h"
42 #include "AliCDBMetaData.h"
43 #include "AliESDEvent.h"
44 #include "AliHeader.h"
46 #include "AliModule.h"
48 #include "AliQADataMakerRec.h"
49 #include "AliQADataMakerSim.h"
50 #include "AliQADataMakerSteer.h"
51 #include "AliRawReaderDate.h"
52 #include "AliRawReaderFile.h"
53 #include "AliRawReaderRoot.h"
55 #include "AliRunLoader.h"
57 ClassImp(AliQADataMakerSteer)
59 //_____________________________________________________________________________
60 AliQADataMakerSteer::AliQADataMakerSteer(const char* gAliceFilename, const char * name, const char * title) :
69 fGAliceFileName(gAliceFilename),
71 fNumberOfEvents(999999),
74 fRawReaderDelete(kTRUE),
78 fMaxEvents = fNumberOfEvents ;
79 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
80 if (IsSelected(AliQA::GetDetName(iDet))) {
81 fLoader[iDet] = NULL ;
82 fQADataMaker[iDet] = NULL ;
83 fQACycles[iDet] = 999999 ;
88 //_____________________________________________________________________________
89 AliQADataMakerSteer::AliQADataMakerSteer(const AliQADataMakerSteer & qas) :
91 fCurrentEvent(qas.fCurrentEvent),
93 fDetectors(qas.fDetectors),
94 fDetectorsW(qas.fDetectorsW),
98 fGAliceFileName(qas.fGAliceFileName),
99 fMaxEvents(qas.fMaxEvents),
100 fNumberOfEvents(qas.fNumberOfEvents),
101 fRunNumber(qas.fRunNumber),
103 fRawReaderDelete(kTRUE),
107 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
108 fLoader[iDet] = qas.fLoader[iDet] ;
109 fQADataMaker[iDet] = qas.fQADataMaker[iDet] ;
110 fQACycles[iDet] = qas.fQACycles[iDet] ;
114 //_____________________________________________________________________________
115 AliQADataMakerSteer & AliQADataMakerSteer::operator = (const AliQADataMakerSteer & qas)
117 // assignment operator
118 this->~AliQADataMakerSteer() ;
119 new(this) AliQADataMakerSteer(qas) ;
123 //_____________________________________________________________________________
124 AliQADataMakerSteer::~AliQADataMakerSteer()
127 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
128 if (IsSelected(AliQA::GetDetName(iDet))) {
129 fLoader[iDet] = NULL;
130 if (fQADataMaker[iDet]) {
131 (fQADataMaker[iDet])->Finish() ;
132 delete fQADataMaker[iDet] ;
133 fQADataMaker[iDet] = NULL ;
138 if (fRawReaderDelete) {
145 //_____________________________________________________________________________
146 Bool_t AliQADataMakerSteer::DoIt(const AliQA::TASKINDEX taskIndex, const char * mode)
148 // Runs all the QA data Maker for every detector
151 // Fill QA data in event loop
152 for (UInt_t iEvent = 0 ; iEvent < (UInt_t)fMaxEvents ; iEvent++) {
155 if ( iEvent%10 == 0 )
156 AliInfo(Form("processing event %d", iEvent));
157 if ( taskIndex == AliQA::kRAWS ) {
158 if ( !fRawReader->NextEvent() )
160 } else if ( taskIndex == AliQA::kESDS ) {
161 if ( fESDTree->GetEntry(iEvent) == 0 )
164 if ( fRunLoader->GetEvent(iEvent) != 0 )
167 // loop over detectors
168 TObjArray* detArray = NULL ;
169 if (fRunLoader) // check if RunLoader exists
170 if ( fRunLoader->GetAliRun() ) // check if AliRun exists in gAlice.root
171 detArray = fRunLoader->GetAliRun()->Detectors() ;
172 for (UInt_t iDet = 0 ; iDet < fgkNDetectors ; iDet++) {
174 AliModule* det = static_cast<AliModule*>(detArray->FindObject(AliQA::GetDetName(iDet))) ;
175 if (!det || !det->IsActive())
178 if (!IsSelected(AliQA::GetDetName(iDet)))
180 AliQADataMaker * qadm = GetQADataMaker(iDet, mode) ;
184 if ( qadm->IsCycleDone() ) {
185 qadm->EndOfCycle(AliQA::kRAWS) ;
186 qadm->StartOfCycle(AliQA::kRAWS) ;
191 qadm->Exec(taskIndex, fRawReader) ;
194 GetLoader(iDet)->LoadHits() ;
195 data = GetLoader(iDet)->TreeH() ;
197 AliWarning(Form(" Hit Tree not found for %s", AliQA::GetDetName(iDet))) ;
199 qadm->Exec(taskIndex, data) ;
202 case AliQA::kSDIGITS :
203 GetLoader(iDet)->LoadSDigits() ;
204 data = GetLoader(iDet)->TreeS() ;
206 AliWarning(Form(" SDigit Tree not found for %s", AliQA::GetDetName(iDet))) ;
208 qadm->Exec(taskIndex, data) ;
211 case AliQA::kDIGITS :
212 GetLoader(iDet)->LoadDigits() ;
213 data = GetLoader(iDet)->TreeD() ;
215 AliWarning(Form(" Digit Tree not found for %s", AliQA::GetDetName(iDet))) ;
217 qadm->Exec(taskIndex, data) ;
220 case AliQA::kRECPOINTS :
221 GetLoader(iDet)->LoadRecPoints() ;
222 data = GetLoader(iDet)->TreeR() ;
224 AliWarning(Form("RecPoints not found for %s", AliQA::GetDetName(iDet))) ;
226 qadm->Exec(taskIndex, data) ;
229 case AliQA::kTRACKSEGMENTS :
231 case AliQA::kRECPARTICLES :
234 qadm->Exec(taskIndex, fESD) ;
236 case AliQA::kNTASKINDEX :
239 //qadm->Increment() ;
243 // // Save QA data for all detectors
244 rv = Finish(taskIndex, mode) ;
249 //_____________________________________________________________________________
250 Bool_t AliQADataMakerSteer::Finish(const AliQA::TASKINDEX taskIndex, const char * mode)
252 // write output to file for all detectors
253 for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
254 if (IsSelected(AliQA::GetDetName(iDet))) {
255 AliQADataMaker * qadm = GetQADataMaker(iDet, mode) ;
257 qadm->EndOfCycle(taskIndex) ;
264 //_____________________________________________________________________________
265 TObjArray * AliQADataMakerSteer::GetFromOCDB(AliQA::DETECTORINDEX det, AliQA::TASKINDEX task) const
267 // Retrieve the list of QA data for a given detector and a given task
268 TObjArray * rv = NULL ;
269 TString tmp(AliQA::GetQARefStorage()) ;
270 if ( tmp.IsNull() ) {
271 AliError("No storage defined, use AliQA::SetQARefStorage") ;
274 AliCDBManager* man = AliCDBManager::Instance() ;
275 if ( ! man->IsDefaultStorageSet() ) {
276 man->SetDefaultStorage(AliQA::GetQARefDefaultStorage()) ;
277 man->SetSpecificStorage(Form("%s/*", AliQA::GetQAOCDBDirName()), AliQA::GetQARefStorage()) ;
279 char detOCDBDir[10] ;
280 sprintf(detOCDBDir, "%s/%s/%s", AliQA::GetQAOCDBDirName(), AliQA::GetDetName((Int_t)det), AliQA::GetRefOCDBDirName()) ;
281 AliInfo(Form("Retrieving reference data from %s/%s for %s", AliQA::GetQARefStorage(), detOCDBDir, AliQA::GetTaskName(task).Data())) ;
282 AliCDBEntry* entry = man->Get(detOCDBDir, 0) ; //FIXME 0 --> Run Number
283 TList * listDetQAD = dynamic_cast<TList *>(entry->GetObject()) ;
285 rv = dynamic_cast<TObjArray *>(listDetQAD->FindObject(AliQA::GetTaskName(task))) ;
289 //_____________________________________________________________________________
290 AliLoader * AliQADataMakerSteer::GetLoader(Int_t iDet)
292 // get the loader for a detector
294 TString detName = AliQA::GetDetName(iDet) ;
295 fLoader[iDet] = fRunLoader->GetLoader(detName + "Loader");
297 return fLoader[iDet] ;
299 // load the QA data maker object
300 TPluginManager* pluginManager = gROOT->GetPluginManager() ;
301 TString loaderName = "Ali" + detName + "Loader" ;
303 AliLoader * loader = NULL ;
304 // first check if a plugin is defined for the quality assurance data maker
305 TPluginHandler* pluginHandler = pluginManager->FindHandler("AliLoader", detName) ;
306 // if not, add a plugin for it
307 if (!pluginHandler) {
308 AliDebug(1, Form("defining plugin for %s", loaderName.Data())) ;
309 TString libs = gSystem->GetLibraries() ;
310 if (libs.Contains("lib" + detName + "base.so") || (gSystem->Load("lib" + detName + "base.so") >= 0)) {
311 pluginManager->AddHandler("AliQADataMaker", detName, loaderName, detName + "loader", loaderName + "()") ;
313 pluginManager->AddHandler("AliLoader", detName, loaderName, detName, loaderName + "()") ;
315 pluginHandler = pluginManager->FindHandler("AliLoader", detName) ;
317 if (pluginHandler && (pluginHandler->LoadPlugin() == 0)) {
318 loader = (AliLoader *) pluginHandler->ExecPlugin(0) ;
321 fLoader[iDet] = loader ;
325 //_____________________________________________________________________________
326 AliQADataMaker * AliQADataMakerSteer::GetQADataMaker(const Int_t iDet, const char * mode )
328 // get the quality assurance data maker for a detector
330 if (fQADataMaker[iDet])
331 return fQADataMaker[iDet] ;
333 AliQADataMaker * qadm = NULL ;
336 // load the QA data maker object
337 TPluginManager* pluginManager = gROOT->GetPluginManager() ;
338 TString detName = AliQA::GetDetName(iDet) ;
340 if (tmp.Contains("sim"))
341 tmp.ReplaceAll("s", "S") ;
342 else if (tmp.Contains("rec"))
343 tmp.ReplaceAll("r", "R") ;
344 TString qadmName = "Ali" + detName + "QADataMaker" + tmp ;
346 // first check if a plugin is defined for the quality assurance data maker
347 TPluginHandler* pluginHandler = pluginManager->FindHandler("AliQADataMaker", detName) ;
348 // if not, add a plugin for it
349 if (!pluginHandler) {
350 AliDebug(1, Form("defining plugin for %s", qadmName.Data())) ;
351 TString libs = gSystem->GetLibraries() ;
352 if (libs.Contains("lib" + detName + mode + ".so") || (gSystem->Load("lib" + detName + mode + ".so") >= 0)) {
353 pluginManager->AddHandler("AliQADataMaker", detName, qadmName, detName + "qadm", qadmName + "()") ;
355 pluginManager->AddHandler("AliQADataMaker", detName, qadmName, detName, qadmName + "()") ;
357 pluginHandler = pluginManager->FindHandler("AliQADataMaker", detName) ;
359 if (pluginHandler && (pluginHandler->LoadPlugin() == 0)) {
360 qadm = (AliQADataMaker *) pluginHandler->ExecPlugin(0) ;
363 fQADataMaker[iDet] = qadm ;
368 //_____________________________________________________________________________
369 Bool_t AliQADataMakerSteer::Init(const AliQA::TASKINDEX taskIndex, const char * mode, const char * input )
371 // Initialize the event source and QA data makers
373 if (taskIndex == AliQA::kRAWS) {
375 TString fileName(input);
376 if (fileName.EndsWith("/")) {
377 fRawReader = new AliRawReaderFile(fileName);
378 } else if (fileName.EndsWith(".root")) {
379 fRawReader = new AliRawReaderRoot(fileName);
380 } else if (!fileName.IsNull()) {
381 fRawReader = new AliRawReaderDate(fileName);
382 fRawReader->SelectEvents(7);
387 fRawReaderDelete = kTRUE ;
388 fRawReader->NextEvent() ;
389 fRunNumber = fRawReader->GetRunNumber() ;
390 AliCDBManager::Instance()->SetRun(fRunNumber) ;
391 fRawReader->RewindEvents();
392 fNumberOfEvents = 999999 ;
393 if ( fMaxEvents < 0 )
394 fMaxEvents = fNumberOfEvents ;
395 } else if (taskIndex == AliQA::kESDS) {
396 if (!gSystem->AccessPathName("AliESDs.root")) { // AliESDs.root exists
397 TFile * esdFile = TFile::Open("AliESDs.root") ;
398 fESDTree = dynamic_cast<TTree *> (esdFile->Get("esdTree")) ;
400 AliError("esdTree not found") ;
403 fESD = new AliESDEvent() ;
404 fESD->ReadFromTree(fESDTree) ;
405 fESDTree->GetEntry(0) ;
406 fRunNumber = fESD->GetRunNumber() ;
407 fNumberOfEvents = fESDTree->GetEntries() ;
408 if ( fMaxEvents < 0 )
409 fMaxEvents = fNumberOfEvents ;
412 AliError("AliESDs.root not found") ;
416 if ( !InitRunLoader() ) {
417 AliWarning("No Run Loader not found") ;
419 fNumberOfEvents = fRunLoader->GetNumberOfEvents() ;
420 if ( fMaxEvents < 0 )
421 fMaxEvents = fNumberOfEvents ;
425 // Initialize all QA data makers for all detectors
426 for (UInt_t iDet = 0; iDet < fgkNDetectors ; iDet++) {
428 if (IsSelected(AliQA::GetDetName(iDet))) {
429 AliQADataMaker * qadm = GetQADataMaker(iDet, mode) ;
431 AliError(Form("AliQADataMaker not found for %s", AliQA::GetDetName(iDet))) ;
432 fDetectorsW.ReplaceAll(AliQA::GetDetName(iDet), "") ;
434 AliDebug(1, Form("Data Maker found for %s", qadm->GetName())) ;
435 qadm->Init(taskIndex, fRunNumber, GetQACycles(iDet)) ;
436 qadm->StartOfCycle(taskIndex, fCycleSame) ;
444 //_____________________________________________________________________________
445 Bool_t AliQADataMakerSteer::InitRunLoader()
447 // get or create the run loader
453 if (!gSystem->AccessPathName(fGAliceFileName.Data())) { // galice.root exists
454 // load all base libraries to get the loader classes
455 TString libs = gSystem->GetLibraries() ;
456 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
457 if (!IsSelected(AliQA::GetDetName(iDet)))
459 TString detName = AliQA::GetDetName(iDet) ;
460 if (detName == "HLT")
462 if (libs.Contains("lib" + detName + "base.so"))
464 gSystem->Load("lib" + detName + "base.so");
466 fRunLoader = AliRunLoader::Open(fGAliceFileName.Data());
468 AliError(Form("no run loader found in file %s", fGAliceFileName.Data()));
471 fRunLoader->CdGAFile();
472 if (fRunLoader->LoadgAlice() == 0) {
473 gAlice = fRunLoader->GetAliRun();
477 AliError(Form("no gAlice object found in file %s", fGAliceFileName.Data()));
481 } else { // galice.root does not exist
482 AliError(Form("the file %s does not exist", fGAliceFileName.Data()));
489 //_____________________________________________________________________________
490 Bool_t AliQADataMakerSteer::IsSelected(const char * det)
492 // check whether detName is contained in detectors
493 // if yes, it is removed from detectors
495 const TString detName(det) ;
496 // check if all detectors are selected
497 if ((fDetectors.CompareTo("ALL") == 0) ||
498 fDetectors.BeginsWith("ALL ") ||
499 fDetectors.EndsWith(" ALL") ||
500 fDetectors.Contains(" ALL ")) {
505 // search for the given detector
507 //AliInfo(Form("SSSSSSSSSSSSS fd = %s det = %s ", fDetectors.Data(), det)) ;
508 if ((fDetectors.CompareTo(detName) == 0) ||
509 fDetectors.BeginsWith(detName+" ") ||
510 fDetectors.EndsWith(" "+detName) ||
511 fDetectors.Contains(" "+detName+" ")) {
512 // fDetectors.ReplaceAll(detName, "");
516 // clean up the detectors string
517 // while (fDetectors.Contains(" "))
518 // fDetectors.ReplaceAll(" ", " ");
519 // while (fDetectors.BeginsWith(" "))
520 // fDetectors.Remove(0, 1);
521 // while (fDetectors.EndsWith(" "))
522 // fDetectors.Remove(fDetectors.Length()-1, 1);
527 //_____________________________________________________________________________
528 Bool_t AliQADataMakerSteer::Merge(const Int_t runNumber) const
530 // Merge all the cycles from all detectors in one single file per run
532 if ( runNumber == -1 )
533 sprintf(cmd, ".! ls *%s*.*.*.root > tempo.txt", AliQA::GetQADataFileName()) ;
535 sprintf(cmd, ".! ls *%s*.%d.*.root > tempo.txt", AliQA::GetQADataFileName(), runNumber) ;
536 gROOT->ProcessLine(cmd) ;
537 ifstream in("tempo.txt") ;
538 const Int_t runMax = 10 ;
539 TString file[AliQA::kNDET*runMax] ;
540 Int_t run[AliQA::kNDET*runMax] ;
547 AliInfo(Form("index = %d file = %s", index, (file[index]).Data())) ;
552 AliError(Form("run number %d not found", runNumber)) ;
557 Int_t runIndexMax = 0 ;
559 sprintf(stmp, ".%s.", AliQA::GetQADataFileName()) ;
560 for (Int_t ifile = 0 ; ifile < index ; ifile++) {
561 TString tmp(file[ifile]) ;
562 tmp.ReplaceAll(".root", "") ;
563 TString det = tmp(0, tmp.Index(".")) ;
564 tmp.Remove(0, tmp.Index(stmp)+4) ;
565 TString ttmp = tmp(0, tmp.Index(".")) ;
566 Int_t newRun = ttmp.Atoi() ;
567 for (Int_t irun = 0; irun <= runIndexMax; irun++) {
568 if (newRun == run[irun])
570 run[runIndex] = newRun ;
573 runIndexMax = runIndex ;
574 ttmp = tmp(tmp.Index(".")+1, tmp.Length()) ;
575 Int_t cycle = ttmp.Atoi() ;
576 AliDebug(1, Form("%s : det = %s run = %d cycle = %d \n", file[ifile].Data(), det.Data(), newRun, cycle)) ;
578 for (Int_t irun = 0 ; irun < runIndexMax ; irun++) {
580 char outFileName[20] ;
581 sprintf(outFileName, "Merged.%s.%d.root", AliQA::GetQADataFileName(), run[irun]) ;
582 merger.OutputFile(outFileName) ;
583 for (Int_t ifile = 0 ; ifile < index-1 ; ifile++) {
585 sprintf(pattern, "%s.%d.", AliQA::GetQADataFileName(), run[irun]) ;
586 TString tmp(file[ifile]) ;
587 if (tmp.Contains(pattern)) {
588 merger.AddFile(tmp) ;
597 //_____________________________________________________________________________
598 void AliQADataMakerSteer::Reset(const Bool_t sameCycle)
600 // Reset the default data members
601 for (UInt_t iDet = 0; iDet < fgkNDetectors; iDet++) {
602 if (IsSelected(AliQA::GetDetName(iDet))) {
603 fLoader[iDet] = NULL;
604 if (fQADataMaker[iDet]) {
605 (fQADataMaker[iDet])->Reset(sameCycle) ;
606 //delete fQADataMaker[iDet] ;
607 //fQADataMaker[iDet] = NULL ;
612 if (fRawReaderDelete) {
617 fCycleSame = sameCycle ;
621 fNumberOfEvents = 999999 ;
624 //_____________________________________________________________________________
625 TString AliQADataMakerSteer::Run(const char * detectors, AliRawReader * rawReader)
627 //Runs all the QA data Maker for Raws only
629 fRawReader = rawReader ;
630 fDetectors = detectors ;
631 fDetectorsW = detectors ;
633 if ( !Init(AliQA::kRAWS, "rec") )
635 fRawReaderDelete = kFALSE ;
637 DoIt(AliQA::kRAWS, "rec") ;
641 //_____________________________________________________________________________
642 TString AliQADataMakerSteer::Run(const char * detectors, const char * fileName)
644 //Runs all the QA data Maker for Raws only
646 fDetectors = detectors ;
647 fDetectorsW = detectors ;
649 if ( !Init(AliQA::kRAWS, "rec", fileName) )
652 DoIt(AliQA::kRAWS, "rec") ;
656 //_____________________________________________________________________________
657 TString AliQADataMakerSteer::Run(const char * detectors, const AliQA::TASKINDEX taskIndex, const char * fileName )
659 // Runs all the QA data Maker for every detector
661 fDetectors = detectors ;
662 fDetectorsW = detectors ;
665 if ( (taskIndex == AliQA::kHITS) || (taskIndex == AliQA::kSDIGITS) || (taskIndex == AliQA::kDIGITS) )
667 else if ( (taskIndex == AliQA::kRAWS) || (taskIndex == AliQA::kRECPOINTS) || (taskIndex == AliQA::kESDS) )
670 AliError(Form("%s not implemented", AliQA::GetTaskName(taskIndex).Data())) ;
674 if ( !Init(taskIndex, mode, fileName) )
677 DoIt(taskIndex, mode) ;
683 //_____________________________________________________________________________
684 Bool_t AliQADataMakerSteer::Save2OCDB(const Int_t runNumber, const Int_t cycleNumber, const char * detectors) const
686 // take the locasl QA data merge into a single file and save in OCDB
688 TString tmp(AliQA::GetQARefStorage()) ;
689 if ( tmp.IsNull() ) {
690 AliError("No storage defined, use AliQA::SetQARefStorage") ;
693 if ( !(tmp.Contains(AliQA::GetLabLocalOCDB()) || tmp.Contains(AliQA::GetLabAliEnOCDB())) ) {
694 AliError(Form("%s is a wrong storage, use %s or %s", AliQA::GetQARefStorage(), AliQA::GetLabLocalOCDB().Data(), AliQA::GetLabAliEnOCDB().Data())) ;
697 TString sdet(detectors) ;
700 if ( sdet.Contains("ALL") ) {
701 rv = Merge(runNumber) ;
704 char inputFileName[20] ;
705 sprintf(inputFileName, "Merged.%s.%d.root", AliQA::GetQADataFileName(), runNumber) ;
706 inputFile = TFile::Open(inputFileName) ;
707 rv = SaveIt2OCDB(runNumber, inputFile) ;
709 for (Int_t index = 0; index < AliQA::kNDET; index++) {
710 if (sdet.Contains(AliQA::GetDetName(index))) {
711 char inputFileName[20] ;
712 sprintf(inputFileName, "%s.%s.%d.%d.root", AliQA::GetDetName(index), AliQA::GetQADataFileName(), runNumber, cycleNumber) ;
713 inputFile = TFile::Open(inputFileName) ;
714 rv *= SaveIt2OCDB(runNumber, inputFile) ;
721 //_____________________________________________________________________________
722 Bool_t AliQADataMakerSteer::SaveIt2OCDB(const Int_t runNumber, TFile * inputFile) const
724 // reads the TH1 from file and adds it to appropriate list before saving to OCDB
726 AliInfo(Form("Saving TH1s in %s to %s", inputFile->GetName(), AliQA::GetQARefStorage())) ;
727 AliCDBManager* man = AliCDBManager::Instance() ;
728 if ( ! man->IsDefaultStorageSet() ) {
729 man->SetDefaultStorage(AliQA::GetQARefDefaultStorage()) ;
730 man->SetSpecificStorage(Form("%s/*", AliQA::GetQAOCDBDirName()), AliQA::GetQARefStorage()) ;
732 if(man->GetRun() < 0)
733 man->SetRun(runNumber);
735 for ( Int_t detIndex = 0 ; detIndex < AliQA::kNDET ; detIndex++) {
736 TDirectory * detDir = inputFile->GetDirectory(AliQA::GetDetName(detIndex)) ;
738 AliInfo(Form("Entering %s", detDir->GetName())) ;
739 char detOCDBDir[20] ;
740 sprintf(detOCDBDir, "%s/%s/%s", AliQA::GetQAOCDBDirName(), AliQA::GetDetName(detIndex), AliQA::GetRefOCDBDirName()) ;
741 AliCDBId idr(detOCDBDir, runNumber, 999999999) ;
742 TList * listDetQAD = new TList() ;
744 sprintf(listName, "%s QA data Reference", AliQA::GetDetName(detIndex)) ;
745 listDetQAD->SetName(listName) ;
746 TList * taskList = detDir->GetListOfKeys() ;
747 TIter nextTask(taskList) ;
749 while ( (taskKey = dynamic_cast<TKey*>(nextTask())) ) {
750 TDirectory * taskDir = detDir->GetDirectory(taskKey->GetName()) ;
751 AliInfo(Form("Saving %s", taskDir->GetName())) ;
752 TObjArray * listTaskQAD = new TObjArray(100) ;
753 listTaskQAD->SetName(taskKey->GetName()) ;
754 listDetQAD->Add(listTaskQAD) ;
755 TList * histList = taskDir->GetListOfKeys() ;
756 TIter nextHist(histList) ;
758 while ( (histKey = dynamic_cast<TKey*>(nextHist())) ) {
759 TObject * odata = taskDir->Get(histKey->GetName()) ;
760 if ( odata->IsA()->InheritsFrom("TH1") ) {
761 AliInfo(Form("Adding %s", histKey->GetName())) ;
762 TH1 * hdata = static_cast<TH1*>(odata) ;
763 listTaskQAD->Add(hdata) ;
768 man->Put(listDetQAD, idr, &mdr) ;