/************************************************************************** * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * * * * Author: The ALICE Off-line Project. * * Contributors are mentioned in the code where appropriate. * * * * Permission to use, copy, modify and distribute this software and its * * documentation strictly for non-commercial purposes is hereby granted * * without fee, provided that the above copyright notice appears in all * * copies and that both the copyright notice and this permission notice * * appear in the supporting documentation. The authors make no claims * * about the suitability of this software for any purpose. It is * * provided "as is" without express or implied warranty. * **************************************************************************/ /* $Id$ */ // Author: Andrei Gheata, 31/05/2006 //============================================================================== // AliAnalysysDataContainer - Container of data of arbitrary type deriving // from TObject used for analysis. A container must be connected to the // output data slot of a single analysis task (producer) , but also as // input slot for possibly several other tasks (consumers). The connected // slots must enforce the same data type as the container (or a derived type). // A container becomes the owner of the contained data once this was produced. // // Containers should be defined by the analysis module using: // // AliAnalysisModule::AddContainer(const char *name, TClass *type); // // A container should be connected to a producer: // AliAnalysisModule::ConnectOutput(AliAnalysisTask *task, // AliAnalysisDataContainer *cont) // and to its consumers: // // AliAnalysisModule::ConnectInput(AliAnalysisTask *task, Int_t islot, // AliAnalysisDataContainer *cont) // // The container will create an implicit connection between the producer task // and all consumers, which will become sub-tasks of the producer. // //============================================================================== #include #include #include #include #include #include #include #include "AliAnalysisManager.h" #include "AliAnalysisDataContainer.h" #include "AliAnalysisDataSlot.h" #include "AliAnalysisTask.h" ClassImp(AliAnalysisDataContainer) //______________________________________________________________________________ AliAnalysisDataContainer::AliAnalysisDataContainer() : TNamed(), fDataReady(kFALSE), fOwnedData(kFALSE), fFileName(), fFolderName(), fFile(NULL), fData(NULL), fType(NULL), fProducer(NULL), fConsumers(NULL) { // Dummy ctor. } //______________________________________________________________________________ AliAnalysisDataContainer::AliAnalysisDataContainer(const char *name, TClass *type) :TNamed(name,""), fDataReady(kFALSE), fOwnedData(kFALSE), fFileName(), fFolderName(), fFile(NULL), fData(NULL), fType(type), fProducer(NULL), fConsumers(NULL) { // Default constructor. SetTitle(fType->GetName()); } //______________________________________________________________________________ AliAnalysisDataContainer::AliAnalysisDataContainer(const AliAnalysisDataContainer &cont) :TNamed(cont), fDataReady(cont.fDataReady), fOwnedData(kFALSE), fFileName(cont.fFileName), fFolderName(cont.fFolderName), fFile(NULL), fData(cont.fData), fType(NULL), fProducer(cont.fProducer), fConsumers(NULL) { // Copy ctor. GetType(); if (cont.fConsumers) { fConsumers = new TObjArray(2); Int_t ncons = cont.fConsumers->GetEntriesFast(); for (Int_t i=0; iAdd(cont.fConsumers->At(i)); } } //______________________________________________________________________________ AliAnalysisDataContainer::~AliAnalysisDataContainer() { // Destructor. Deletes data ! (What happens if data is a container ???) if (fData && fOwnedData) delete fData; if (fConsumers) delete fConsumers; } //______________________________________________________________________________ AliAnalysisDataContainer &AliAnalysisDataContainer::operator=(const AliAnalysisDataContainer &cont) { // Assignment. if (&cont != this) { TNamed::operator=(cont); fDataReady = cont.fDataReady; fOwnedData = kFALSE; fFileName = cont.fFileName; fFolderName = cont.fFolderName; fFile = NULL; fData = cont.fData; GetType(); fProducer = cont.fProducer; if (cont.fConsumers) { fConsumers = new TObjArray(2); Int_t ncons = cont.fConsumers->GetEntriesFast(); for (Int_t i=0; iAdd(cont.fConsumers->At(i)); } } return *this; } //______________________________________________________________________________ void AliAnalysisDataContainer::AddConsumer(AliAnalysisTask *consumer, Int_t islot) { // Add a consumer for contained data; AliAnalysisDataSlot *slot = consumer->GetInputSlot(islot); if (!slot || !slot->GetType()) { cout<<"Consumer task "<< consumer->GetName()<<" does not have an input/type #"<GetName(),islot)); return; } if (!slot->GetType()->InheritsFrom(GetType())) { cout<<"Data type "<GetTitle()<<" for input slot "<GetName()<<" does not match container type "<GetType()->GetName(),islot,consumer->GetName(),fType->GetName())); return; } if (!fConsumers) fConsumers = new TObjArray(2); fConsumers->Add(consumer); // Add the consumer task to the list of task of the producer if (fProducer && !fProducer->GetListOfTasks()->FindObject(consumer)) fProducer->Add(consumer); } //______________________________________________________________________________ Bool_t AliAnalysisDataContainer::ClientsExecuted() const { // Check if all client tasks have executed. TIter next(fConsumers); AliAnalysisTask *task; while ((task=(AliAnalysisTask*)next())) { if (!task->HasExecuted()) return kFALSE; } return kTRUE; } //______________________________________________________________________________ void AliAnalysisDataContainer::DeleteData() { // Delete data if not needed anymore. if (!fDataReady || !ClientsExecuted()) { cout<<"Data not ready or not all clients of container "<SetType(gROOT->GetClass(fTitle.Data())); if (!fType) printf("AliAnalysisDataContainer: Unknown class: %s\n", GetTitle()); return fType; } //______________________________________________________________________________ void AliAnalysisDataContainer::GetEntry(Long64_t ientry) { // If data is ready and derives from TTree or from TBranch, this will get the // requested entry in memory if not already loaded. if (!fDataReady || !GetType()) return; Bool_t istree = fType->InheritsFrom(TTree::Class()); if (istree) { TTree *tree = (TTree*)fData; if (tree->GetReadEntry() != ientry) tree->GetEntry(ientry); return; } Bool_t isbranch = fType->InheritsFrom(TBranch::Class()); if (isbranch) { TBranch *branch = (TBranch*)fData; if (branch->GetReadEntry() != ientry) branch->GetEntry(ientry); return; } } //______________________________________________________________________________ Long64_t AliAnalysisDataContainer::Merge(TCollection *list) { // Merge a list of containers with this one. Containers in the list must have // data of the same type. if (!list || !fData) return 0; printf("Merging %d containers %s\n", list->GetSize()+1, GetName()); TMethodCall callEnv; if (fData->IsA()) callEnv.InitWithPrototype(fData->IsA(), "Merge", "TCollection*"); if (!callEnv.IsValid() && !list->IsEmpty()) { cout << "No merge interface for data stored by " << GetName() << ". Merging not possible !" << endl; return 1; } if (list->IsEmpty()) return 1; TIter next(list); AliAnalysisDataContainer *cont; // Make a list where to temporary store the data to be merged. TList *collectionData = new TList(); Int_t count = 0; // object counter while ((cont=(AliAnalysisDataContainer*)next())) { TObject *data = cont->GetData(); if (!data) continue; if (strcmp(cont->GetName(), GetName())) { cout << "Not merging containers with different names !" << endl; continue; } printf(" ... merging object %s\n", data->GetName()); collectionData->Add(data); count++; } callEnv.SetParam((Long_t) collectionData); callEnv.Execute(fData); delete collectionData; return count+1; } //______________________________________________________________________________ void AliAnalysisDataContainer::PrintContainer(Option_t *option, Int_t indent) const { // Print info about this container. TString ind; for (Int_t i=0; iGetName()); else printf("%s= No data producer",ind.Data()); printf("%s = Consumer tasks: ", ind.Data()); if (!fConsumers || !fConsumers->GetEntriesFast()) printf("-none-\n"); else printf("\n"); } if (fFolderName.Length()) printf("Filename: %s folder: %s\n", fFileName.Data(), fFolderName.Data()); else printf("Filename: %s\n", fFileName.Data()); TIter next(fConsumers); AliAnalysisTask *task; while ((task=(AliAnalysisTask*)next())) task->PrintTask(option, indent+3); } //______________________________________________________________________________ Bool_t AliAnalysisDataContainer::SetData(TObject *data, Option_t *) { // Set the data as READY only if it was published by the producer. // If there is no producer declared, this is a top level container. AliAnalysisTask *task; Bool_t init = kFALSE; Int_t i, nc; if (!fProducer) { if (data != fData) init = kTRUE; fData = data; fDataReady = kTRUE; if (fConsumers) { nc = fConsumers->GetEntriesFast(); for (i=0; iAt(i); task->CheckNotify(init); } } return kTRUE; } // Check if it is the producer who published the data if (fProducer->GetPublishedData()==data) { fData = data; fDataReady = kTRUE; if (fConsumers) { nc = fConsumers->GetEntriesFast(); for (i=0; iAt(i); task->CheckNotify(); } } return kTRUE; } else { // Ignore data posting from other than the producer // cout<<"Data for container "<GetName()<GetName())); return kFALSE; } } //______________________________________________________________________________ void AliAnalysisDataContainer::SetFileName(const char *filename) { // The filename field can be actually composed by the actual file name followed // by :dirname (optional): // filename = file_name[:dirname] // No slashes (/) allowed fFileName = filename; fFolderName = ""; Int_t index = fFileName.Index(":"); // Fill the folder name if (index >= 0) { fFolderName = fFileName(index+1, fFileName.Length()-index); fFileName.Remove(index); } if (!fFileName.Length()) Fatal("SetFileName", "Empty file name"); if (fFileName.Index("/")>=0) Fatal("SetFileName", "No slashes (/) allowed in the file name"); } //______________________________________________________________________________ void AliAnalysisDataContainer::SetProducer(AliAnalysisTask *prod, Int_t islot) { // Set the producer of data. The slot number is required for data type checking. if (fProducer) { cout<<"Data container "<GetName()<GetName())); } if (fDataReady) { cout<GetOutputSlot(islot); if (!slot) { cout<<"Producer task "<GetName()<<" does not have an output #"<GetName(),islot)); return; } if (!slot->GetType()->InheritsFrom(GetType())) { cout<<"Data type "<GetTitle()<<"for output slot "<GetName()<<" does not match container type "<GetType()->GetName(),islot,prod->GetName(),fType->GetName())); return; } fProducer = prod; // Add all consumers as daughter tasks TIter next(fConsumers); AliAnalysisTask *cons; while ((cons=(AliAnalysisTask*)next())) { if (!prod->GetListOfTasks()->FindObject(cons)) prod->Add(cons); } } //______________________________________________________________________________ AliAnalysisDataWrapper *AliAnalysisDataContainer::ExportData() const { // Wraps data for sending it through the net. AliAnalysisDataWrapper *pack = 0; if (!fData) { Error("ExportData", "Container %s - No data to be wrapped !", GetName()); return pack; } AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager(); if (mgr->GetDebugLevel() > 1) printf(" ExportData: Wrapping data %s for container %s\n", fData->GetName(),GetName()); pack = new AliAnalysisDataWrapper(fData); pack->SetName(fName.Data()); return pack; } //______________________________________________________________________________ void AliAnalysisDataContainer::ImportData(AliAnalysisDataWrapper *pack) { // Unwraps data from a data wrapper. if (pack) { fData = pack->Data(); if (!fData) { Error("ImportData", "No data was wrapped for container %s", GetName()); return; } AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager(); if (mgr->GetDebugLevel() > 1) printf(" ImportData: Unwrapping data %s for container %s\n", fData->GetName(),GetName()); fDataReady = kTRUE; // Imported wrappers do not own data anymore (AG 13-11-07) pack->SetDeleteData(kFALSE); } } ClassImp (AliAnalysisDataWrapper) //______________________________________________________________________________ AliAnalysisDataWrapper::AliAnalysisDataWrapper(TObject *data) :TNamed(), fData(data) { // Ctor. if (data) SetName(data->GetName()); } //______________________________________________________________________________ AliAnalysisDataWrapper::~AliAnalysisDataWrapper() { // Dtor. if (fData && TObject::TestBit(kDeleteData)) delete fData; } //______________________________________________________________________________ AliAnalysisDataWrapper &AliAnalysisDataWrapper::operator=(const AliAnalysisDataWrapper &other) { // Assignment. if (&other != this) { TNamed::operator=(other); fData = other.fData; } return *this; } //______________________________________________________________________________ Long64_t AliAnalysisDataWrapper::Merge(TCollection *list) { // Merge a list of containers with this one. Containers in the list must have // data of the same type. if (TH1::AddDirectoryStatus()) TH1::AddDirectory(kFALSE); if (!fData) return 0; if (!list || list->IsEmpty()) return 1; SetDeleteData(); TMethodCall callEnv; if (fData->IsA()) callEnv.InitWithPrototype(fData->IsA(), "Merge", "TCollection*"); if (!callEnv.IsValid()) { cout << "No merge interface for data stored by " << GetName() << ". Merging not possible !" << endl; return 1; } TIter next1(list); AliAnalysisDataWrapper *cont; // Make a list where to temporary store the data to be merged. TList *collectionData = new TList(); Int_t count = 0; // object counter // printf("Wrapper %s 0x%lx (data=%s) merged with:\n", GetName(), (ULong_t)this, fData->ClassName()); while ((cont=(AliAnalysisDataWrapper*)next1())) { cont->SetDeleteData(); TObject *data = cont->Data(); if (!data) continue; // printf(" - %s 0x%lx (data=%s)\n", cont->GetName(), (ULong_t)cont, data->ClassName()); collectionData->Add(data); count++; } callEnv.SetParam((Long_t) collectionData); callEnv.Execute(fData); delete collectionData; return count+1; }