X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PHOS%2FAliPHOSGetter.cxx;h=faf49034194611546583293df7fb75147d729c91;hb=77f353c4a4be6534d09627cc35855605eb156713;hp=f634efc1487c1f0ff92e20c6e8e41578dffb5b48;hpb=f4f43939b91b1a20cbd2ed81eceea9fe0c368dd1;p=u%2Fmrichter%2FAliRoot.git diff --git a/PHOS/AliPHOSGetter.cxx b/PHOS/AliPHOSGetter.cxx index f634efc1487..faf49034194 100644 --- a/PHOS/AliPHOSGetter.cxx +++ b/PHOS/AliPHOSGetter.cxx @@ -20,6 +20,10 @@ Everywhere reading the treese TTree->GetEvent(i) is replaced by reading the branches TBranch->GetEntry(0) */ +/* $Log: + 08.2002 Dmitri Peressounko: + +*/ //_________________________________________________________________________ // A singleton. This class should be used in the analysis stage to get @@ -35,7 +39,7 @@ // for(Int_t irecp = 0; irecp < gime->NRecParticles() ; irecp++) // AliPHOSRecParticle * part = gime->RecParticle(1) ; // ................ -// please->GetEvent(event) ; // reads new event from galice.root +// gime->Event(event) ; // reads new event from galice.root // //*-- Author: Yves Schutz (SUBATECH) & Dmitri Peressounko (RRC KI & SUBATECH) //*-- Completely redesigned by Dmitri Peressounko March 2001 @@ -44,9 +48,7 @@ //*-- systematic usage of TFolders without changing the interface ////////////////////////////////////////////////////////////////////////////// - // --- ROOT system --- - #include "TFile.h" #include "TTree.h" #include "TROOT.h" @@ -55,41 +57,43 @@ #include "TParticle.h" // --- Standard library --- -#include // --- AliRoot header files --- - #include "AliRun.h" #include "AliConfig.h" #include "AliPHOSGetter.h" #include "AliPHOS.h" #include "AliPHOSDigitizer.h" #include "AliPHOSSDigitizer.h" -#include "AliPHOSClusterizer.h" #include "AliPHOSClusterizerv1.h" -#include "AliPHOSTrackSegmentMaker.h" #include "AliPHOSTrackSegmentMakerv1.h" #include "AliPHOSTrackSegment.h" -#include "AliPHOSPID.h" #include "AliPHOSPIDv1.h" #include "AliPHOSGeometry.h" - +#include "AliPHOSRaw2Digits.h" +#include "AliPHOSCalibrationDB.h" ClassImp(AliPHOSGetter) AliPHOSGetter * AliPHOSGetter::fgObjGetter = 0 ; + TFile * AliPHOSGetter::fFile = 0 ; //____________________________________________________________________________ -AliPHOSGetter::AliPHOSGetter(const char* headerFile, const char* branchTitle ) +AliPHOSGetter::AliPHOSGetter(const char* headerFile, const char* branchTitle, const Bool_t toSplit ) { - //Initialize all lists + // This is the ctor called by GetInstance and the only one that can be used - fHeaderFile = headerFile ; - fBranchTitle = branchTitle ; - fSDigitsTitle = branchTitle ; - fDigitsTitle = branchTitle ; - fRecPointsTitle = branchTitle ; - fRecParticlesTitle = branchTitle ; - fTrackSegmentsTitle = branchTitle ; + if( fHeaderFile.Contains("_") ) { + Fatal("AliPHOSGetter", "Invalid file name (_ not allowed) %s", fHeaderFile.Data() ) ; + } + + //Initialize all data + + fFailed = kFALSE ; + fDebug = 0 ; + fAlice = 0 ; + + fToSplit = toSplit ; + fHeaderFile = headerFile ; fPrimaries = new TObjArray(1) ; @@ -102,49 +106,93 @@ AliPHOSGetter::AliPHOSGetter(const char* headerFile, const char* branchTitle ) fQAFolder = dynamic_cast(gROOT->FindObjectAny("Folders/Run/Conditions/QA")); fTasksFolder = dynamic_cast(gROOT->FindObjectAny("Folders/Tasks")) ; + //Set titles to branches and create PHOS specific folders + SetTitle(branchTitle) ; + if ( fHeaderFile != "aliroot" ) { // to call the getter without a file - //open headers file - TFile * file = static_cast(gROOT->GetFile(fHeaderFile.Data() ) ) ; + fFile = static_cast(gROOT->GetFile(fHeaderFile.Data() ) ) ; - if(file == 0){ //if file was not opened yet, read gAlice - if(fHeaderFile.Contains("rfio")) // if we read file using HPSS - file = TFile::Open(fHeaderFile.Data(),"update") ; - else - file = new TFile(fHeaderFile.Data(),"update") ; - - if (!file->IsOpen()) { - cerr << "ERROR : AliPHOSGetter::AliPHOSGetter -> Cannot open " << fHeaderFile.Data() << endl ; - abort() ; + if(!fFile) { //if file was not opened yet, read gAlice + fFile = TFile::Open(fHeaderFile.Data(),"update") ; + if (!fFile->IsOpen()) { + Error("AliPHOSGetter", "Cannot open %s", fHeaderFile.Data() ) ; + fFailed = kTRUE ; + return ; } - - gAlice = static_cast(file->Get("gAlice")) ; - } + gAlice = static_cast(fFile->Get("gAlice")) ; + } } - + if (!gAlice) { - cerr << "ERROR : AliPHOSGetter::AliPHOSGetter -> Cannot find gAlice in " << fHeaderFile.Data() << endl ; - abort() ; + Error("AliPHOSGetter", "Cannot find gAlice in %s", fHeaderFile.Data() ) ; + fFailed = kTRUE ; + return ; } if (!PHOS()) { if (fDebug) - cout << "INFO: AliPHOSGetter -> Posting PHOS to Folders" << endl ; - AliConfig * conf = AliConfig::Instance() ; - conf->Add(static_cast(gAlice->GetDetector("PHOS"))) ; - conf->Add(static_cast(gAlice->GetDetector("PHOS"))) ; + Info("AliPHOSGetter", "-> Posting PHOS to Folders") ; + if (gAlice->GetDetector("PHOS")) { + AliConfig * conf = AliConfig::Instance() ; + conf->Add(static_cast(gAlice->GetDetector("PHOS"))) ; + conf->Add(static_cast(gAlice->GetDetector("PHOS"))) ; + } + else + Error("AliPHOSGetter", "detector PHOS not found") ; } - + fDebug=0; } + //____________________________________________________________________________ AliPHOSGetter::~AliPHOSGetter(){ + if (fPrimaries) { + fPrimaries->Delete() ; + delete fPrimaries ; + } + + TFolder * phosF = dynamic_cast(fSDigitsFolder->FindObject("PHOS")) ; + TCollection * folderslist = phosF->GetListOfFolders() ; + TIter next(folderslist) ; + TFolder * folder = 0 ; + while ( (folder = static_cast(next())) ) + phosF->Remove(folder) ; + + if (fFile) { + fFile->Close() ; + delete fFile ; + fFile = 0 ; + } + fgObjGetter = 0 ; + } //____________________________________________________________________________ -void AliPHOSGetter::CreateWhiteBoard() const +void AliPHOSGetter::CloseFile() { + delete gAlice ; + gAlice = 0 ; + delete fAlice ; + fAlice = 0 ; +} +//____________________________________________________________________________ +const TFolder * AliPHOSGetter::Folder(const TString what) const { + + // returns the PHOS folder required by what + // what = hits, sdigits, digits + + if ( what == "hits" ) + return dynamic_cast(fHitsFolder->FindObject("PHOS")) ; + else if ( what == "sdigits" ) + return dynamic_cast(fSDigitsFolder->FindObject("PHOS")) ; + else if ( what == "digits" ) + return dynamic_cast(fDigitsFolder->FindObject("PHOS")) ; + else { + Error("GetFolder", " %s illegal option (hits, sdigits, digits)", what.Data() ) ; + return 0 ; + } } //____________________________________________________________________________ @@ -152,36 +200,244 @@ AliPHOSGetter * AliPHOSGetter::GetInstance() { // Returns the pointer of the unique instance already defined - AliPHOSGetter * rv = 0 ; - if ( fgObjGetter ) - rv = fgObjGetter ; - else - cout << "AliPHOSGetter::GetInstance ERROR: not yet initialized" << endl ; - - return rv ; + if ( fgObjGetter ) { + return fgObjGetter ; + } + else { + //Warning("GetInstance", "not yet initialized") ; + return 0 ; + } } //____________________________________________________________________________ AliPHOSGetter * AliPHOSGetter::GetInstance(const char* headerFile, - const char* branchTitle) + const char* branchTitle, + const Bool_t toSplit) { // Creates and returns the pointer of the unique instance - // Must be called only when the environment has changed + // Must be called only when the environment has changed + + if(!fgObjGetter){ + fgObjGetter = new AliPHOSGetter(headerFile,branchTitle,toSplit) ; + if(fgObjGetter->fFailed) + return 0; + else + return fgObjGetter ; + } - if ( fgObjGetter ) - if((fgObjGetter->fBranchTitle.CompareTo(branchTitle) == 0) && - (fgObjGetter->fHeaderFile.CompareTo(headerFile)==0)) + //First checks, if header file already opened + if(!fgObjGetter->fFile){ + fgObjGetter = new AliPHOSGetter(headerFile,branchTitle,toSplit) ; + if(fgObjGetter->fFailed) + return 0; + else return fgObjGetter ; + } + + if(fgObjGetter->fHeaderFile.CompareTo(headerFile)==0){ //Opened the same header file + if((fgObjGetter->fBranchTitle.CompareTo(branchTitle) == 0)&& //Open the same branch title + (toSplit==fgObjGetter->fToSplit)){ //Nothing should be cleaned + } + else{ //Clean all data and AliPHOS...zers + if(fgObjGetter->fToSplit) + fgObjGetter->CloseSplitFiles() ; + fgObjGetter->CleanWhiteBoard() ; + fgObjGetter->fToSplit = toSplit ; + fgObjGetter->SetTitle(branchTitle) ; + } + } + else{ //Close already opened files, clean memory and open new header file + if(gAlice) + delete gAlice ; + gAlice= 0; + if(fgObjGetter->fFile){ + fgObjGetter->fFile->Close() ; + fgObjGetter->fFile=0; + } + if(fgObjGetter->fToSplit) + fgObjGetter->CloseSplitFiles() ; + fgObjGetter->CleanWhiteBoard() ; + fgObjGetter = new AliPHOSGetter(headerFile,branchTitle,toSplit) ; + } + return fgObjGetter ; + +} + +//____________________________________________________________________________ +const Bool_t AliPHOSGetter::BranchExists(const TString recName) const +{ + //Looks in the tree Tree"name" if branch with current name olready exists + + TString filename("") ; + TString name, dataname, zername ; + if(recName == "SDigits"){ + filename=fSDigitsFileName ; + name = "TreeS0" ; + dataname = "PHOS" ; + zername = "AliPHOSSDigitizer" ; + } + else + if(recName == "Digits"){ + filename=fDigitsFileName ; + name = "TreeD0" ; + dataname = "PHOS" ; + zername = "AliPHOSDigitizer" ; + } else - fgObjGetter->~AliPHOSGetter() ; // delete it if already exists another version + if(recName == "RecPoints"){ + filename=fRecPointsFileName ; + name = "TreeR0" ; + dataname = "PHOSEmcRP" ; + zername = "AliPHOSClusterizer" ; + } + else + if(recName == "TrackSegments"){ + filename=fTrackSegmentsFileName ; + name = "TreeR0" ; + dataname = "PHOSTS" ; + zername = "AliPHOSTrackSegmentMaker" ; + } + else + if(recName == "RecParticles"){ + filename= fRecParticlesFileName ; + name = "TreeR0" ; + dataname = "PHOSRP" ; + zername = "AliPHOSPID" ; + } + else + return kFALSE ; + + TFile * file ; + TTree * tree ; + if(fToSplit){ + file = static_cast(gROOT->GetFile(filename.Data() ) ) ; + if(!file) + file = TFile::Open(fSDigitsFileName.Data(),"update"); + } + else + file = fFile ; + + tree = (TTree *)file->Get(name.Data()) ; + + if(!tree ) + return kFALSE ; + + TObjArray * lob = static_cast(tree->GetListOfBranches()) ; + TIter next(lob) ; + TBranch * branch = 0 ; + TString titleName(fBranchTitle); + titleName+=":"; + while ((branch = (static_cast(next())))) { + TString branchName(branch->GetName() ) ; + TString branchTitle(branch->GetTitle() ) ; + if ( branchName.BeginsWith(dataname) && branchTitle.BeginsWith(fBranchTitle) ){ + Warning("BranchExists", "branch %s with title %s already exists in %s", dataname.Data(), fBranchTitle.Data(), name.Data() ) ; + return kTRUE ; + } + if ( branchName.BeginsWith(zername) && branchTitle.BeginsWith(titleName) ){ + Warning("BranchExists", "branch AliPHOS... with title %s already exists in %s", branch->GetTitle(), name.Data() ) ; + return kTRUE ; + } + } + //We can't delete three if gAlice points to it... To be redisigned somehow???!!! + if(!fToSplit){ + if(name.Contains("TreeS")) + if(tree!=gAlice->TreeS()) + tree->Delete(); + if(name.Contains("TreeD")) + if(tree!=gAlice->TreeD()) + tree->Delete(); + if(name.Contains("TreeR")) + if(tree!=gAlice->TreeR()) + tree->Delete(); + } + return kFALSE ; - fgObjGetter = new AliPHOSGetter(headerFile,branchTitle) ; +} + +//____________________________________________________________________________ +void AliPHOSGetter::ListBranches(Int_t event) const +{ - // Posts a few item to the white board (folders) - // fgObjGetter->CreateWhiteBoard() ; - - return fgObjGetter ; + TBranch * branch = 0 ; + if (gAlice->GetEvent(event) == -1) + return ; + TTree * t = gAlice->TreeH() ; + if(t){ + Info("ListBranches", "-> ****** Hits : ") ; + TObjArray * lob = t->GetListOfBranches() ; + TIter next(lob) ; + while ( (branch = static_cast(next())) ) + Info("ListBranches", " %s", branch->GetName()) ; + } else + Warning("ListBranches", "TreeH not found for event %d", event ) ; + + t = gAlice->TreeS() ; + if(t){ + Info("ListBranches", "-> ****** SDigits : ") ; + TObjArray * lob = t->GetListOfBranches() ; + TIter next(lob) ; + while ( (branch = static_cast(next())) ) + Info("ListBranches", " %s %s", branch->GetName(), branch->GetTitle()) ; + } else + Warning("ListBranches", "TreeS not found for event %d", event) ; + + + t = gAlice->TreeD() ; + if(t){ + Info("ListBranches", "-> ****** Digits : ") ; + TObjArray * lob = t->GetListOfBranches() ; + TIter next(lob) ; + while ( (branch = static_cast(next())) ) + Info("ListBranches", " %s %s", branch->GetName(), branch->GetTitle()) ; + } else + Warning("ListBranches", "TreeD not found for event %d", event) ; + + + t = gAlice->TreeR() ; + if(t){ + Info("ListBranches", "-> ****** Recon : ") ; + TObjArray * lob = t->GetListOfBranches() ; + TIter next(lob) ; + while ( (branch = static_cast(next())) ) + Info("ListBranches", " %s %s", branch->GetName(), branch->GetTitle()) ; + } else + Warning("ListBranches", "TreeR not found for event %d", event) ; + +} + +//____________________________________________________________________________ +void AliPHOSGetter::NewBranch(TString name, Int_t event) +{ + fBranchTitle = fSDigitsTitle = fDigitsTitle = fRecPointsTitle = fTrackSegmentsTitle = fRecParticlesTitle = name ; + Event(event) ; +} + +//____________________________________________________________________________ +Bool_t AliPHOSGetter::NewFile(TString name) +{ + fHeaderFile = name ; + fFile->Close() ; + fFailed = kFALSE; + + fFile = static_cast(gROOT->GetFile(fHeaderFile.Data() ) ) ; + if(!fFile) { //if file was not opened yet, read gAlice + fFile = TFile::Open(fHeaderFile.Data(),"update") ; + if (!fFile->IsOpen()) { + Error("NewFile", "Cannot open %s", fHeaderFile.Data() ) ; + fFailed = kTRUE ; + return fFailed ; + } + gAlice = static_cast(fFile->Get("gAlice")) ; + } + + if (!gAlice) { + Error("AliPHOSGetter", "Cannot find gAlice in %s", fHeaderFile.Data() ) ; + fFailed = kTRUE ; + return fFailed ; + } + return fFailed ; } //____________________________________________________________________________ @@ -191,7 +447,7 @@ const AliPHOS * AliPHOSGetter::PHOS() AliPHOS * phos = dynamic_cast(fModuleFolder->FindObject("PHOS")) ; if (!phos) if (fDebug) - cout << "WARNING: AliPHOSGetter::PHOS -> PHOS module not found in Folders" << endl ; + Warning("PHOS", "-> PHOS module not found in Folders") ; return phos ; } @@ -205,7 +461,7 @@ const AliPHOSGeometry * AliPHOSGetter::PHOSGeometry() } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostPrimaries(void) const +const Bool_t AliPHOSGetter::PostPrimaries(void) const { //------- Primaries ---------------------- // the hierarchy is //Folders/RunMC/Event/Data/Primaries @@ -213,8 +469,8 @@ Bool_t AliPHOSGetter::PostPrimaries(void) const TFolder * primariesFolder = dynamic_cast(fPrimariesFolder->FindObject("Primaries")) ; if ( !primariesFolder ) { if (fDebug) { - cout << "WARNING: AliPHOSGetter::Post Primaries -> Folder //" << fPrimariesFolder << "/Primaries/ not found!" << endl; - cout << "INFO: AliPHOSGetter::Post Primaries -> Adding Folder //" << fPrimariesFolder << "/Primaries/" << endl; + Warning("PostPrimaries", "-> Folder //%s/Primaries/ not found!", fPrimariesFolder->GetName()) ; + Info("PostPrimaries", "-> Adding Folder //%s/Primaries", fPrimariesFolder->GetName()) ; } primariesFolder = fPrimariesFolder->AddFolder("Primaries", "Primaries particles from TreeK") ; } @@ -232,27 +488,24 @@ TObject** AliPHOSGetter::PrimariesRef(void) const // the hierarchy is //Folders/RunMC/Event/Data/Primaries if ( !fPrimariesFolder ) { - cerr << "ERROR: AliPHOSGetter::Post PrimariesRef -> Folder //" << fPrimariesFolder << " not found!" << endl; - return 0; + Fatal("PrimariesRef", "Folder //%s not found", fPrimariesFolder) ; } TFolder * primariesFolder = dynamic_cast(fPrimariesFolder->FindObject("Primaries")) ; if ( !primariesFolder ) { - cerr << "ERROR: AliPHOSGetter::Post PrimariesRef -> Folder //" << fPrimariesFolder << "/Primaries/ not found!" << endl; - return 0; + Fatal("PrimariesRef", "Folder //%s/Primaries/ not found", fPrimariesFolder) ; } TObject * p = primariesFolder->FindObject("Primaries") ; if(!p) { - cerr << "ERROR: AliPHOSGetter::PrimariesRef -> " << primariesFolder->GetName() << "/Primaries not found !" << endl ; - return 0 ; + Fatal("PrimariesRef","%s /Primaries not found !", primariesFolder->GetName() ) ; } - else - return primariesFolder->GetListOfFolders()->GetObjectRef(p) ; + + return primariesFolder->GetListOfFolders()->GetObjectRef(p) ; } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostHits(void) const +const Bool_t AliPHOSGetter::PostHits(void) const { //------- Hits ---------------------- // the hierarchy is //Folders/RunMC/Event/Data/PHOS/Hits @@ -260,8 +513,8 @@ Bool_t AliPHOSGetter::PostHits(void) const TFolder * phosFolder = dynamic_cast(fHitsFolder->FindObject("PHOS")) ; if ( !phosFolder ) { if (fDebug) { - cout << "WARNING: AliPHOSGetter::Post H -> Folder //" << fHitsFolder << "/PHOS/ not found!" << endl; - cout << "INFO: AliPHOSGetter::Post H -> Adding Folder //" << fHitsFolder << "/PHOS/" << endl; + Warning("PostHits", "-> Folder //%s/PHOS/ not found!", fHitsFolder) ; + Info("PostHits", "-> Adding Folder //%s/PHOS/", fHitsFolder) ; } phosFolder = fHitsFolder->AddFolder("PHOS", "Hits from PHOS") ; } @@ -279,27 +532,24 @@ TObject** AliPHOSGetter::HitsRef(void) const // the hierarchy is //Folders/RunMC/Event/Data/PHOS/Hits if ( !fHitsFolder ) { - cerr << "ERROR: AliPHOSGetter::Post H -> Folder //" << fHitsFolder << " not found!" << endl; - return 0; + Fatal("HitsRef", "Folder //%s not found !", fHitsFolder) ; } TFolder * phosFolder = dynamic_cast(fHitsFolder->FindObject("PHOS")) ; if ( !phosFolder ) { - cerr << "ERROR: AliPHOSGetter::Post HRef -> Folder //" << fHitsFolder << "/PHOS/ not found!" << endl; - return 0; + Fatal("HitsRef", "Folder //%s/PHOS/ not found !", fHitsFolder) ; } TObject * h = phosFolder->FindObject("Hits") ; if(!h) { - cerr << "ERROR: AliPHOSGetter::HRef -> " << phosFolder->GetName() << "/Hits not found !" << endl ; - return 0 ; + Fatal("HitsRef", "%s/Hits not fount !", phosFolder->GetName() ) ; } - else - return phosFolder->GetListOfFolders()->GetObjectRef(h) ; + + return phosFolder->GetListOfFolders()->GetObjectRef(h) ; } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostSDigits(const char * name, const char * headerFile) const +const Bool_t AliPHOSGetter::PostSDigits(const char * name, const char * headerFile) const { //---------- SDigits ------------------------- @@ -309,23 +559,21 @@ Bool_t AliPHOSGetter::PostSDigits(const char * name, const char * headerFile) co TFolder * phosFolder = dynamic_cast(fSDigitsFolder->FindObject("PHOS")) ; if ( !phosFolder ) { if (fDebug) { - cout << "WARNING: AliPHOSGetter::Post S -> Folder //" << fSDigitsFolder << "/PHOS/ not found!" << endl; - cout << "INFO: AliPHOSGetter::Post S -> Adding Folder //" << fHitsFolder << "/PHOS/" << endl; + Warning("PostSDigits", "-> Folder //%s/PHOS/ not found!", fSDigitsFolder) ; + Info("PostSDigits", "-> Adding Folder //%s/PHOS/", fHitsFolder) ; } phosFolder = fSDigitsFolder->AddFolder("PHOS", "SDigits from PHOS") ; } + TString subdir(headerFile) ; subdir.ReplaceAll("/","_") ; TFolder * phosSubFolder = dynamic_cast(phosFolder->FindObject(subdir)) ; if ( !phosSubFolder ) phosSubFolder = phosFolder->AddFolder(subdir, ""); + TObject * sd = phosSubFolder->FindObject(name); - if ( sd ) { - if (fDebug) - cerr <<"INFO: AliPHOSGetter::Post S -> Folder " << subdir - << " already exists!" << endl ; - }else{ + if ( !sd ) { TClonesArray * sdigits = new TClonesArray("AliPHOSDigit",1) ; sdigits->SetName(name) ; phosSubFolder->Add(sdigits) ; @@ -340,14 +588,12 @@ TObject** AliPHOSGetter::SDigitsRef(const char * name, const char * file) const // the hierarchy is //Folders/RunMC/Event/Data/PHOS/SDigits/filename/SDigits if ( !fSDigitsFolder ) { - cerr << "ERROR: AliPHOSGetter::Post SRef -> Folder //" << fSDigitsFolder << " not found!" << endl; - return 0; + Fatal("SDigitsRef", "Folder //%s not found !", fSDigitsFolder) ; } TFolder * phosFolder = dynamic_cast(fSDigitsFolder->FindObject("PHOS")) ; if ( !phosFolder ) { - cerr << "ERROR: AliPHOSGetter::Post SRef -> Folder //" << fSDigitsFolder << "/PHOS/ not found!" << endl; - return 0; + Fatal("SDigitsRef", "Folder //%s/PHOS not found !", fSDigitsFolder) ; } TFolder * phosSubFolder = 0 ; @@ -357,20 +603,20 @@ TObject** AliPHOSGetter::SDigitsRef(const char * name, const char * file) const phosSubFolder = dynamic_cast(phosFolder->FindObject(fHeaderFile)) ; if(!phosSubFolder) { - cerr << "ERROR: AliPHOSGetter::Post SRef -> Folder //Folders/RunMC/Event/Data/PHOS/" << file << "not found!" << endl; - return 0; + Fatal("SDigitsRef", "Folder //Folders/RunMC/Event/Data/PHOS/%s not found !", file) ; } TObject * dis = phosSubFolder->FindObject(name) ; - if(!dis) - return 0 ; - else - return phosSubFolder->GetListOfFolders()->GetObjectRef(dis) ; + if(!dis){ + Fatal("SDigitsRef", "object %s not found !", name) ; + } + + return phosSubFolder->GetListOfFolders()->GetObjectRef(dis) ; } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostSDigitizer(AliPHOSSDigitizer * sdigitizer) const +const Bool_t AliPHOSGetter::PostSDigitizer(AliPHOSSDigitizer * sdigitizer) const { //---------- SDigitizer ------------------------- // the hierarchy is //Folders/Tasks/SDigitizer/PHOS/sdigitsname @@ -379,14 +625,14 @@ Bool_t AliPHOSGetter::PostSDigitizer(AliPHOSSDigitizer * sdigitizer) const TTask * sd = dynamic_cast(fTasksFolder->FindObject("SDigitizer")) ; if ( !sd ) { - cerr << "ERROR: AliPHOSGetter::Post Ser -> Task //" << fTasksFolder << "/SDigitizer not found!" << endl; + Error("PostDigitizer", "Task //%s/SDigitizer not found !", fTasksFolder) ; return kFALSE ; } TTask * phos = dynamic_cast(sd->GetListOfTasks()->FindObject("PHOS")) ; if ( !phos ) { if (fDebug) { - cout <<"WARNING: AliPHOSGetter::Post Ser ->//" << fTasksFolder << "/SDigitizer/PHOS/ not found!" << endl; - cout <<"INFO: AliPHOSGetter::Post Ser -> Adding //" << fTasksFolder << "/SDigitizer/PHOS/" << endl; + Warning("PostSDigitizer", "->//%s/SDigitizer/PHOS/ not found!", fTasksFolder) ; + Info("PostSDigitizer", "-> Adding //%s/SDigitizer/PHOS/", fTasksFolder) ; } phos = new TTask("PHOS", "") ; sd->Add(phos) ; @@ -394,7 +640,7 @@ Bool_t AliPHOSGetter::PostSDigitizer(AliPHOSSDigitizer * sdigitizer) const AliPHOSSDigitizer * phossd = dynamic_cast(phos->GetListOfTasks()->FindObject( sdigitizer->GetName() )); if (phossd) { if (fDebug) - cout << "INFO: AliPHOSGetter::Post Ser -> Task " << sdigitizer->GetName() << " already exists" << endl ; + Info("PostSDigitizer", "-> Task %s already exists", sdigitizer->GetName()) ; phos->GetListOfTasks()->Remove(phossd) ; } phos->Add(sdigitizer) ; @@ -408,14 +654,12 @@ TObject** AliPHOSGetter::SDigitizerRef(const char * name) const TTask * sd = dynamic_cast(fTasksFolder->FindObject("SDigitizer")) ; if ( !sd ) { - cerr << "ERROR: AliPHOSGetter::Post SerRef -> Task //" << fTasksFolder << "/SDigitizer not found!" << endl; - abort(); + Fatal("SDigitizerRef", "Task //%s/SDigitizer not found !", fTasksFolder) ; } TTask * phos = dynamic_cast(sd->GetListOfTasks()->FindObject("PHOS")) ; if ( !phos ) { - cerr <<"ERROR: AliPHOSGetter::Post SerRef -> //" << fTasksFolder << "/SDigitizer/PHOS not found!" << endl; - abort(); + Fatal("SDigitizerRef", "//%s/SDigitizer/PHOS not found !", fTasksFolder) ; } TTask * task = dynamic_cast(phos->GetListOfTasks()->FindObject(name)) ; @@ -425,23 +669,22 @@ TObject** AliPHOSGetter::SDigitizerRef(const char * name) const } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostSDigitizer(const char * name, const char * file) const +const Bool_t AliPHOSGetter::PostSDigitizer(const char * name, const char * file) const { //---------- SDigitizer ------------------------- // the hierarchy is //Folders/Tasks/SDigitizer/PHOS/sdigitsname - TTask * sd = dynamic_cast(fTasksFolder->FindObject("SDigitizer")) ; if ( !sd ) { - cerr << "ERROR: AliPHOSGetter::Post Ser -> Task //" << fTasksFolder << "/SDigitizer not found!" << endl; + Error("PostSDigitizer", "Task //%s/SDigitizer not found !", fTasksFolder) ; return kFALSE ; } TTask * phos = dynamic_cast(sd->GetListOfTasks()->FindObject("PHOS")) ; if ( !phos ) { if (fDebug) { - cout <<"WARNING: AliPHOSGetter::Post Ser -> //" << fTasksFolder << "/SDigitizer/PHOS/ not found!" << endl; - cout <<"INFO: AliPHOSGetter::Post Ser -> Adding //" << fTasksFolder << "/SDigitizer/PHOS" << endl; + Error("PostSDigitizer", "-> //%s/SDigitizer/PHOS/ not found!", fTasksFolder) ; + Info("PostSDigitizer", "-> Adding //%s/SDigitizer/PHOS", fTasksFolder) ; } phos = new TTask("PHOS", "") ; sd->Add(phos) ; @@ -450,10 +693,11 @@ Bool_t AliPHOSGetter::PostSDigitizer(const char * name, const char * file) const TString sdname(name) ; sdname.Append(":") ; sdname.Append(file); + sdname.ReplaceAll("/","_") ; AliPHOSSDigitizer * phossd = dynamic_cast(phos->GetListOfTasks()->FindObject( sdname )); if (!phossd) { phossd = new AliPHOSSDigitizer() ; - //Note, we can not call constructor with parameters: it will call Getter and scrud up everething + //Note, we can not call constructor with parameters: it will call Getter and scew up everething phossd->SetName(sdname) ; phossd->SetTitle(file) ; phos->Add(phossd) ; @@ -463,7 +707,7 @@ Bool_t AliPHOSGetter::PostSDigitizer(const char * name, const char * file) const } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostDigits(const char * name) const +const Bool_t AliPHOSGetter::PostDigits(const char * name) const { //---------- Digits ------------------------- // the hierarchy is //Folders/Run/Event/Data/PHOS/SDigits/name @@ -472,8 +716,8 @@ Bool_t AliPHOSGetter::PostDigits(const char * name) const if ( !phosFolder ) { if (fDebug) { - cout << "WARNING: AliPHOSGetter::Post D -> Folder //" << fDigitsFolder << "/PHOS/ not found!" << endl; - cout << "INFO: AliPHOSGetter::Post D -> Adding Folder //" << fDigitsFolder << "/PHOS/" << endl; + Warning("PostDigitizer", "-> Folder //%s/PHOS/ not found!", fDigitsFolder) ; + Info("PostDigitizer", "-> Adding Folder //%s/PHOS/", fDigitsFolder) ; } phosFolder = fDigitsFolder->AddFolder("PHOS", "Digits from PHOS") ; } @@ -494,39 +738,38 @@ TObject** AliPHOSGetter::DigitsRef(const char * name) const // the hierarchy is //Folders/Run/Event/Data/PHOS/Digits/name if ( !fDigitsFolder ) { - cerr << "ERROR: AliPHOSGetter::Post DRef -> Folder //" << fDigitsFolder << " not found!" << endl; - return 0; + Fatal("DigitsRef", "Folder //%s not found !", fDigitsFolder) ; } TFolder * phosFolder = dynamic_cast(fDigitsFolder->FindObject("PHOS")) ; if ( !phosFolder ) { - cerr << "ERROR: AliPHOSGetter::DRef -> Folder //" << fDigitsFolder << "/PHOS/ not found!" << endl; - return 0; + Fatal("DigitsRef", "Folder //%s/PHOS/ not found !", fDigitsFolder) ; } TObject * d = phosFolder->FindObject(name) ; - if(!d) - return 0 ; - else - return phosFolder->GetListOfFolders()->GetObjectRef(d) ; + if(!d) { + Fatal("DigitsRef", "object %s not found !", name) ; + } + + return phosFolder->GetListOfFolders()->GetObjectRef(d) ; } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostDigitizer(AliPHOSDigitizer * digitizer) const +const Bool_t AliPHOSGetter::PostDigitizer(AliPHOSDigitizer * digitizer) const { //---------- Digitizer ------------------------- TTask * sd = dynamic_cast(fTasksFolder->FindObject("Digitizer")) ; if ( !sd ) { - cerr << "ERROR: AliPHOSGetter::Post Der -> Task //" << fTasksFolder << "/Digitizer not found!" << endl; + Error("PostDigitizer", "Task //%s/Digitizer not found !", fTasksFolder) ; return kFALSE ; } TTask * phos = dynamic_cast(sd->GetListOfTasks()->FindObject("PHOS")) ; if ( !phos ) { if (fDebug) { - cout <<"WARNING: AliPHOSGetter::Post Der -> //" << fTasksFolder << "/Digitizer/PHOS not found!" << endl; - cout <<"INFO: AliPHOSGetter::Post Der -> Adding //" << fTasksFolder << "/Digitizer/PHOS" << endl; + Error("PostDigitizer", "//%s/Digitizer/PHOS not found!", fTasksFolder) ; + Info("PostDigitizer", "Adding //%s/Digitizer/PHOS", fTasksFolder) ; } phos = new TTask("PHOS", "") ; sd->Add(phos) ; @@ -542,33 +785,41 @@ Bool_t AliPHOSGetter::PostDigitizer(AliPHOSDigitizer * digitizer) const } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostDigitizer(const char * name) const +const Bool_t AliPHOSGetter::PostDigitizer(const char * name) const { //---------- Digitizer ------------------------- // the hierarchy is //Folders/Tasks/SDigitizer/PHOS/sdigitsname TTask * d = dynamic_cast(fTasksFolder->FindObject("Digitizer")) ; if ( !d ) { - cerr << "ERROR: AliPHOSGetter::Post Der -> Task //" << fTasksFolder << "/Digitizer not found!" << endl; + Error("PostDigitizer", "Task //%s/Digitizer not found !", fTasksFolder) ; return kFALSE ; } TTask * phos = dynamic_cast(d->GetListOfTasks()->FindObject("PHOS")) ; if ( !phos ) { if (fDebug) { - cout <<"WARNING: AliPHOSGetter::Post Der -> //" << fTasksFolder << "/Digitizer/PHOS not found!" << endl; - cout <<"INFO: AliPHOSGetter::Post Der -> Adding //" << fTasksFolder << "/Digitizer/PHOS" << endl; + Warning("PostDigitizer", "//%s/Digitizer/PHOS not found!", fTasksFolder) ; + Info("PostDigitizer", "Adding //%s/Digitizer/PHOS", fTasksFolder) ; } phos = new TTask("PHOS", "") ; d->Add(phos) ; } - AliPHOSDigitizer * phosd = dynamic_cast(phos->GetListOfTasks()->FindObject(name)) ; + TTask * phosd = dynamic_cast(phos->GetListOfTasks()->FindObject(fDigitsTitle)) ; if (!phosd) { - phosd = new AliPHOSDigitizer() ; - phosd->SetName(fDigitsTitle) ; - phosd->SetTitle(fHeaderFile) ; - phos->Add(phosd) ; + if(strcmp(name, "Digitizer")==0){ + phosd = new AliPHOSDigitizer() ; + phosd->SetName(fDigitsTitle) ; + phosd->SetTitle(fHeaderFile) ; + phos->Add(phosd) ; + } + else{ + phosd = new AliPHOSRaw2Digits() ; + phosd->SetName(fDigitsTitle) ; + phosd->SetTitle(fHeaderFile) ; + phos->Add(phosd) ; + } } return kTRUE; } @@ -578,14 +829,12 @@ TObject** AliPHOSGetter::DigitizerRef(const char * name) const { TTask * sd = dynamic_cast(fTasksFolder->FindObject("Digitizer")) ; if ( !sd ) { - cerr << "ERROR: AliPHOSGetter::Post DerRef -> Task //" << fTasksFolder << "/Digitizer not found!" << endl; - abort(); + Fatal("DigitizerRef", "Task //%s/Digitizer not found !", fTasksFolder) ; } TTask * phos = dynamic_cast(sd->GetListOfTasks()->FindObject("PHOS")) ; if ( !phos ) { - cerr <<"ERROR: AliPHOSGetter::Post DerRef -> //" << fTasksFolder << "/Digitizer/PHOS" << endl; - abort(); + Fatal("DigitizerRef", "//%s/Digitizer/PHOS", fTasksFolder) ; } TTask * task = dynamic_cast(phos->GetListOfTasks()->FindObject(name)) ; @@ -595,7 +844,7 @@ TObject** AliPHOSGetter::DigitizerRef(const char * name) const } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostRecPoints(const char * name) const +const Bool_t AliPHOSGetter::PostRecPoints(const char * name) const { // -------------- RecPoints ------------------------------------------- // the hierarchy is //Folders/Run/Event/RecData/PHOS/EMCARecPoints/name @@ -605,8 +854,8 @@ Bool_t AliPHOSGetter::PostRecPoints(const char * name) const if ( !phosFolder ) { if (fDebug) { - cout << "WARNING: AliPHOSGetter::Post RPo -> Folder //" << fRecoFolder << "/PHOS/ not found!" << endl; - cout << "INFO: AliPHOSGetter::Post Rpo -> Adding Folder //" << fRecoFolder << "/PHOS/" << endl; + Warning("PostRecPoints", "-> Folder //%s/PHOS/ not found!", fRecoFolder->GetName()) ; + Info("PostRecPoints", "-> Adding Folder //%s/PHOS/", fRecoFolder->GetName()) ; } phosFolder = fRecoFolder->AddFolder("PHOS", "Reconstructed data from PHOS") ; } @@ -615,8 +864,8 @@ Bool_t AliPHOSGetter::PostRecPoints(const char * name) const TFolder * phosRPoEMCAFolder = dynamic_cast(phosFolder->FindObject("EMCARecPoints")) ; if ( !phosRPoEMCAFolder ) { if (fDebug) { - cout << "WARNING: AliPHOSGetter::Post RPo -> Folder //" << fRecoFolder << "/PHOS/EMCARecPoints/ not found!" << endl; - cout << "INFO: AliPHOSGetter::Post Rpo -> Adding Folder //" << fRecoFolder << "/PHOS/EMCARecPoints not found!" << endl; + Warning("PostRecPoints", "-> Folder //%s/PHOS/EMCARecPoints/ not found!", fRecoFolder->GetName()) ; + Info("PostRecPoints", "-> Adding Folder //%s/PHOS/EMCARecPoints", fRecoFolder->GetName()) ; } phosRPoEMCAFolder = phosFolder->AddFolder("EMCARecPoints", "EMCA RecPoints from PHOS") ; } @@ -632,8 +881,8 @@ Bool_t AliPHOSGetter::PostRecPoints(const char * name) const TFolder * phosRPoCPVFolder = dynamic_cast(phosFolder->FindObject("CPVRecPoints")) ; if ( !phosRPoCPVFolder ) { if (fDebug) { - cout << "WARNING: AliPHOSGetter::Post RPo -> Folder //" << fRecoFolder << "/PHOS/CPVRecPoints/ not found!" << endl; - cout << "INFO: AliPHOSGetter::Post Rpo -> Adding Folder //" << fRecoFolder << "/PHOS/CPVRecPoints/" << endl; + Warning("PostRecPoints", "-> Folder //%s/PHOS/CPVRecPoints/ not found!", fRecoFolder->GetName()) ; + Info("PostRecPoints", "-> Adding Folder //%s/PHOS/CPVRecPoints/", fRecoFolder->GetName()) ; } phosRPoCPVFolder = phosFolder->AddFolder("CPVRecPoints", "CPV RecPoints from PHOS") ; } @@ -654,23 +903,21 @@ TObject** AliPHOSGetter::EmcRecPointsRef(const char * name) const // the hierarchy is //Folders/Run/Event/RecData/PHOS/EMCARecPoints/name if ( !fRecoFolder ) { - cerr << "ERROR: AliPHOSGetter::EmcRecPointsRef -> Folder //" << fRecoFolder << " not found!" << endl; - return 0 ; + Fatal("EmcRecPointsRef", "Folder //%s not found !", fRecoFolder->GetName() ) ; } TFolder * phosFolder = dynamic_cast(fRecoFolder->FindObject("PHOS/EMCARecPoints")) ; if ( !phosFolder ) { - cerr << "ERROR: AliPHOSGetter::EmcRecPointsRef -> Folder //" << fRecoFolder << "/PHOS/EMCARecPoints/ not found!" << endl; - return 0; + Fatal("EmcRecPointsRef", "Folder //%s/PHOS/EMCARecPoints/ not found !", fRecoFolder->GetName() ) ; } TObject * erp = phosFolder->FindObject(name ) ; if ( !erp ) { - return 0 ; + Fatal("EmcRecPointsRef", "object %s not found !", name) ; } return phosFolder->GetListOfFolders()->GetObjectRef(erp) ; - + } //____________________________________________________________________________ @@ -680,26 +927,24 @@ TObject** AliPHOSGetter::CpvRecPointsRef(const char * name) const // the hierarchy is //Folders/Run/Event/RecData/PHOS/CPVRecPoints/name if ( !fRecoFolder ) { - cerr << "ERROR: AliPHOSGetter::EmcRecPointsRef -> Folder //" << fRecoFolder << " not found!" << endl; - return 0 ; + Fatal("CpvRecPointsRef", "Folder //%s not found !", fRecoFolder->GetName() ) ; } TFolder * phosFolder = dynamic_cast(fRecoFolder->FindObject("PHOS/CPVRecPoints")) ; if ( !phosFolder ) { - cerr << "ERROR: AliPHOSGetter::CpvRecPointsRef -> Folder //" << fRecoFolder << "/PHOS/CPVRecPoints/" << endl; - return 0; + Fatal("CpvRecPointsRef", "Folder //%s/PHOS/CPVRecPoints/ not found !", fRecoFolder->GetName() ) ; } TObject * crp = phosFolder->FindObject(name ) ; if ( !crp ) { - return 0 ; + Fatal("CpvRecPointsRef", "object %s nott found", name ) ; } return phosFolder->GetListOfFolders()->GetObjectRef(crp) ; - + } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostClusterizer(AliPHOSClusterizer * clu) const +const Bool_t AliPHOSGetter::PostClusterizer(AliPHOSClusterizer * clu) const { // ------------------ AliPHOSClusterizer ------------------------ // the hierarchy is //Folders/Tasks/Reconstructioner/PHOS/sdigitsname @@ -707,15 +952,15 @@ Bool_t AliPHOSGetter::PostClusterizer(AliPHOSClusterizer * clu) const TTask * tasks = dynamic_cast(fTasksFolder->FindObject("Reconstructioner")) ; if ( !tasks ) { - cerr << "ERROR: AliPHOSGetter::Post Rer -> Task //" << fTasksFolder << "/Reconstructioner not found!" << endl; + Error("PostClusterizer", "Task //%s/Reconstructioner not found !", fTasksFolder) ; return kFALSE ; } TTask * phos = dynamic_cast(tasks->GetListOfTasks()->FindObject("PHOS")) ; if ( !phos ) { if (fDebug) { - cout <<"WARNING: AliPHOSGetter::Post Rer -> //" << fTasksFolder << "/ReconstructionerPHOS not found!" << endl; - cout <<"INFO: AliPHOSGetter::Post Rer -> Adding //" << fTasksFolder << "/Reconstructioner/PHOS" << endl; + Warning("PostClusterizer", "//%s/Reconstructioner/PHOS not found!", fTasksFolder) ; + Info("PostClusterizer", "Adding //%s/Reconstructioner/PHOS", fTasksFolder) ; } phos = new TTask("PHOS", "") ; tasks->Add(phos) ; @@ -724,7 +969,7 @@ Bool_t AliPHOSGetter::PostClusterizer(AliPHOSClusterizer * clu) const AliPHOSClusterizer * phoscl = dynamic_cast(phos->GetListOfTasks()->FindObject(clu->GetName())) ; if (phoscl) { if (fDebug) - cout << "INFO: AliPHOSGetter::Post Rer -> Task " << clu->GetName() << " already exists" << endl ; + Info("PostClusterizer", "Task %s already exists", clu->GetName()) ; phoscl->Delete() ; phos->GetListOfTasks()->Remove(phoscl) ; } @@ -739,14 +984,12 @@ TObject** AliPHOSGetter::ClusterizerRef(const char * name) const TTask * tasks = dynamic_cast(fTasksFolder->FindObject("Reconstructioner")) ; if ( !tasks ) { - cerr << "ERROR: AliPHOSGetter::Post RerRef -> Task //" << fTasksFolder << "/Reconstructioner not found!" << endl; - return kFALSE ; + Fatal("ClusterizerRef", "Task //%s/Reconstructioner not found !", fTasksFolder->GetName() ) ; } TTask * phos = dynamic_cast(tasks->GetListOfTasks()->FindObject("PHOS")) ; if ( !phos ) { - cerr <<"WARNING: AliPHOSGetter::Post RerRef -> //" << fTasksFolder << "/Reconstructioner/PHOS" << endl; - return 0 ; + Fatal("ClusterizerRef", " //%s/Reconstructioner/PHOS not founf !", fTasksFolder->GetName() ) ; } TList * l = phos->GetListOfTasks() ; @@ -754,7 +997,7 @@ TObject** AliPHOSGetter::ClusterizerRef(const char * name) const TTask * task ; TTask * clu = 0 ; TString cluname(name) ; - cluname+=":clu-" ; + cluname+=":clu" ; while((task = static_cast(it.Next()) )){ TString taskname(task->GetName()) ; if(taskname.BeginsWith(cluname)){ @@ -763,14 +1006,16 @@ TObject** AliPHOSGetter::ClusterizerRef(const char * name) const } } - if(clu) - return l->GetObjectRef(clu) ; - else - return 0 ; + if(!clu) { + Fatal("ClusterizerRef", "Task //%s/Reconstructioner/clusterizer/%s not found", fTasksFolder->GetName(), name) ; + } + + return l->GetObjectRef(clu) ; + } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostClusterizer(const char * name) const +const Bool_t AliPHOSGetter::PostClusterizer(const char * name) const { // ------------------ AliPHOSClusterizer ------------------------ // the hierarchy is //Folders/Tasks/Reconstructioner/PHOS/sdigitsname @@ -778,31 +1023,42 @@ Bool_t AliPHOSGetter::PostClusterizer(const char * name) const TTask * tasks = dynamic_cast(fTasksFolder->FindObject("Reconstructioner")) ; if ( !tasks ) { - cerr << "ERROR: AliPHOSGetter::Post Rer -> Task//" << fTasksFolder << "/Reconstructioner not found!" << endl; + Error("PostClusterizer", "Task//%s/Reconstructioner not found !", fTasksFolder) ; return kFALSE ; } TTask * phos = dynamic_cast(tasks->GetListOfTasks()->FindObject("PHOS")) ; if ( !phos ) { if (fDebug) { - cout <<"WARNING: AliPHOSGetter::Post Rer -> //" << fTasksFolder << "/Reconstructioner/PHOS not found!" << endl; - cout <<"INFO: AliPHOSGetter::Post Rer -> Adding //" << fTasksFolder << "/Reconstructioner/PHOS" << endl; + Warning("PostClusterizer", "//%s/Reconstructioner/PHOS not found!", fTasksFolder) ; + Info("PostClusterizer", "Adding //%s/Reconstructioner/PHOS", fTasksFolder) ; } phos = new TTask("PHOS", "") ; tasks->Add(phos) ; } - AliPHOSClusterizer * phoscl = new AliPHOSClusterizerv1() ; + TList * l = phos->GetListOfTasks() ; + TIter it(l) ; TString clun(name) ; - clun+=":clu-v1" ; + clun+=":clu" ; + TTask * task ; + while((task = static_cast(it.Next()) )){ + TString taskname(task->GetName()) ; + if(taskname.BeginsWith(clun)) + return kTRUE ; + } + + AliPHOSClusterizerv1 * phoscl = new AliPHOSClusterizerv1() ; + clun+="-v1" ; phoscl->SetName(clun) ; + phoscl->SetTitle(fHeaderFile) ; phos->Add(phoscl) ; return kTRUE; } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostTrackSegments(const char * name) const +const Bool_t AliPHOSGetter::PostTrackSegments(const char * name) const { // ---------------TrackSegments ----------------------------------- // the hierarchy is //Folders/Run/Event/RecData/PHOS/TrackSegments/name @@ -811,8 +1067,8 @@ Bool_t AliPHOSGetter::PostTrackSegments(const char * name) const if ( !phosFolder ) { if (fDebug) { - cout << "WARNING: AliPHOSGetter::Post TS -> Folder //" << fRecoFolder << "/PHOS/ not found!" << endl; - cout << "INFO: AliPHOSGetter::Post TS -> Adding Folder //" << fRecoFolder << "/PHOS" << endl; + Warning("PostTrackSegments", "-> Folder //%s/PHOS/ not found", fRecoFolder->GetName()) ; + Info("PostTrackSegments", "-> Adding Folder //%s/PHOS", fRecoFolder->GetName()) ; } phosFolder = fRecoFolder->AddFolder("PHOS", "Reconstructed data from PHOS") ; } @@ -820,8 +1076,8 @@ Bool_t AliPHOSGetter::PostTrackSegments(const char * name) const TFolder * phosTSFolder = dynamic_cast(phosFolder->FindObject("TrackSegments")) ; if ( !phosTSFolder ) { if (fDebug) { - cout << "WARNING: AliPHOSGetter::Post TS -> Folder//" << fRecoFolder << "/PHOS/TrackSegments/ not found!" << endl; - cout << "INFO: AliPHOSGetter::Post TS -> Adding Folder //" << fRecoFolder << "/PHOS/TrackSegments/" << endl; + Warning("PostTrackSegments", "-> Folder //%s/PHOS/TrackSegments/ not found!", fRecoFolder->GetName() ) ; + Info("PostTrackSegments", "-> Adding Folder //%s/PHOS/TrackSegments/", fRecoFolder->GetName()) ; } phosTSFolder = phosFolder->AddFolder("TrackSegments", "TrackSegments from PHOS") ; } @@ -842,25 +1098,23 @@ TObject** AliPHOSGetter::TrackSegmentsRef(const char * name) const // the hierarchy is //Folders/Run/Event/RecData/PHOS/TrackSegments/name if ( !fRecoFolder ) { - cerr << "ERROR: AliPHOSGetter::TrackSegmentsRef -> Folder //" << fRecoFolder << "not found!" << endl; - return 0 ; + Fatal("TrackSegmentsRef", "Folder //%s not found !", fRecoFolder->GetName() ) ; } TFolder * phosFolder = dynamic_cast(fRecoFolder->FindObject("PHOS/TrackSegments")) ; if ( !phosFolder ) { - cerr << "ERROR: AliPHOSGetter::TrackSegmentsRef -> Folder //" << fRecoFolder << "/PHOS/TrackSegments/ not found!" << endl; - return 0; + Fatal("TrackSegmentsRef", "Folder //%s/PHOS/TrackSegments/ not found !", fRecoFolder->GetName() ) ; } TObject * tss = phosFolder->FindObject(name) ; if (!tss) { - return 0 ; + Fatal("TrackSegmentsRef", "object %s not found !", name) ; } return phosFolder->GetListOfFolders()->GetObjectRef(tss) ; } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostTrackSegmentMaker(AliPHOSTrackSegmentMaker * tsmaker) const +const Bool_t AliPHOSGetter::PostTrackSegmentMaker(AliPHOSTrackSegmentMaker * tsmaker) const { //------------Track Segment Maker ------------------------------ // the hierarchy is //Folders/Tasks/Reconstructioner/PHOS/sdigitsname @@ -868,15 +1122,15 @@ Bool_t AliPHOSGetter::PostTrackSegmentMaker(AliPHOSTrackSegmentMaker * tsmaker) TTask * tasks = dynamic_cast(fTasksFolder->FindObject("Reconstructioner")) ; if ( !tasks ) { - cerr << "ERROR: AliPHOSGetter::Post Ter -> Task //" << fTasksFolder << "/Reconstructioner not found!" << endl; + Error("PostTrackSegmentMaker", "Task //%s/Reconstructioner not found !", fTasksFolder) ; return kFALSE ; } TTask * phos = dynamic_cast(tasks->GetListOfTasks()->FindObject("PHOS")) ; if ( !phos ) { if (fDebug) { - cout <<"WARNING: AliPHOSGetter::Post Rer -> //" << fTasksFolder << "/Reconstructioner/PHOS not found!" << endl; - cout <<"INFO: AliPHOSGetter::Post Rer -> Adding //" << fTasksFolder << "/Reconstructioner/PHOS" << endl; + Warning("PostTrackSegmentMaker", "//%s/Reconstructioner/PHOS not found!", fTasksFolder) ; + Info("PostTrackSegmentMaker", "Adding //%s/Reconstructioner/PHOS", fTasksFolder) ; } phos = new TTask("PHOS", "") ; tasks->Add(phos) ; @@ -893,38 +1147,44 @@ Bool_t AliPHOSGetter::PostTrackSegmentMaker(AliPHOSTrackSegmentMaker * tsmaker) } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostTrackSegmentMaker(const char * name) const +const Bool_t AliPHOSGetter::PostTrackSegmentMaker(const char * name) const { //------------Track Segment Maker ------------------------------ // the hierarchy is //Folders/Tasks/Reconstructioner/PHOS/sdigitsname - TTask * tasks = dynamic_cast(fTasksFolder->FindObject("Reconstructioner")) ; if ( !tasks ) { - cerr << "ERROR: AliPHOSGetter::Post Ter -> Task //" << fTasksFolder << "/Reconstructioner not found!" << endl; + Error("PostTrackSegmentMaker", "Task //%s/Reconstructioner not found !", fTasksFolder->GetName() ) ; return kFALSE ; } TTask * phos = dynamic_cast(tasks->GetListOfTasks()->FindObject("PHOS")) ; if ( !phos ) { if (fDebug) { - cout <<"WARNING: AliPHOSGetter::Post Rer -> //" << fTasksFolder << "/Reconstructioner/PHOS not found!" << endl; - cout <<"INFO: AliPHOSGetter::Post Rer -> Adding //" << fTasksFolder << "/Reconstructioner/PHOS" << endl; + Warning("PostTrackSegmentMaker", "//%s/Reconstructioner/PHOS not found!", fTasksFolder->GetName() ) ; + Info("PostTrackSegmentMaker", "Adding //%s/Reconstructioner/PHOS", fTasksFolder->GetName()) ; } phos = new TTask("PHOS", "") ; tasks->Add(phos) ; } - AliPHOSTrackSegmentMaker * phosts = - dynamic_cast(phos->GetListOfTasks()->FindObject(name)) ; - if (!phosts) { - phosts = new AliPHOSTrackSegmentMakerv1() ; - TString tsn(name); - tsn+=":tsm-v1" ; - phosts->SetName(tsn) ; - phos->Add(phosts) ; + TList * l = phos->GetListOfTasks() ; + TIter it(l) ; + TString tsn(name); + tsn+=":tsm" ; + TTask * task ; + while((task = static_cast(it.Next()) )){ + TString taskname(task->GetName()) ; + if(taskname.BeginsWith(tsn)) + return kTRUE ; } + + AliPHOSTrackSegmentMakerv1 * phosts = new AliPHOSTrackSegmentMakerv1() ; + tsn+="-v1" ; + phosts->SetName(tsn) ; + phosts->SetTitle(fHeaderFile) ; + phos->Add(phosts) ; return kTRUE; } @@ -936,14 +1196,12 @@ TObject** AliPHOSGetter::TSMakerRef(const char * name) const TTask * tasks = dynamic_cast(fTasksFolder->FindObject("Reconstructioner")) ; if ( !tasks ) { - cerr << "ERROR: AliPHOSGetter::Post TerRef -> Task //" << fTasksFolder << "/Reconstructioner not found!" << endl; - return kFALSE ; + Fatal("TSMakerRef", "Task //%s/Reconstructioner not found !", fTasksFolder->GetName() ) ; } TTask * phos = dynamic_cast(tasks->GetListOfTasks()->FindObject("PHOS")) ; if ( !phos ) { - cerr <<"WARNING: AliPHOSGetter::Post TerRef -> //" << fTasksFolder << "/Reconstructioner/PHOS not found!" << endl; - return 0 ; + Fatal("TSMakerRef", "//%s/Reconstructioner/PHOS not found !", fTasksFolder->GetName() ) ; } TList * l = phos->GetListOfTasks() ; @@ -951,7 +1209,7 @@ TObject** AliPHOSGetter::TSMakerRef(const char * name) const TTask * task ; TTask * tsm = 0 ; TString tsmname(name) ; - tsmname+=":tsm-" ; + tsmname+=":tsm" ; while((task = static_cast(it.Next()) )){ TString taskname(task->GetName()) ; if(taskname.BeginsWith(tsmname)){ @@ -960,25 +1218,26 @@ TObject** AliPHOSGetter::TSMakerRef(const char * name) const } } - if(tsm) - return l->GetObjectRef(tsm) ; - else - return 0 ; - + if(!tsm) { + Fatal("TSMakerRef", "Task //%s/Reconstructioner/PHOS/TrackSegmentMarker/%s not found !", fTasksFolder->GetName(), name) ; + } + + return l->GetObjectRef(tsm) ; + } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostRecParticles(const char * name) const +const Bool_t AliPHOSGetter::PostRecParticles(const char * name) const { // -------------------- RecParticles ------------------------ - // the hierarchy is //Folders/Run/Event/RecData/PHOS/TrackSegments/name + // the hierarchy is //Folders/Run/Event/RecData/PHOS/RecParticles/name TFolder * phosFolder = dynamic_cast(fRecoFolder->FindObject("PHOS")) ; if ( !phosFolder ) { if (fDebug) { - cout << "WARNING: AliPHOSGetter::Post RPa -> Folder //" << fRecoFolder << "/PHOS/ not found!" << endl; - cout << "INFO: AliPHOSGetter::Post Rpa -> Adding Folder //" << fRecoFolder << "/PHOS/" << endl; + Warning("PostRecParticles", "-> Folder //%s/PHOS/ not found!", fRecoFolder->GetName()) ; + Info("PostRecParticles", "-> Adding Folder //%s/PHOS/", fRecoFolder->GetName()) ; } phosFolder = fRecoFolder->AddFolder("PHOS", "Reconstructed data from PHOS") ; } @@ -986,8 +1245,8 @@ Bool_t AliPHOSGetter::PostRecParticles(const char * name) const TFolder * phosRPaFolder = dynamic_cast(phosFolder->FindObject("RecParticles")) ; if ( !phosRPaFolder ) { if (fDebug) { - cout << "WARNING: AliPHOSGetter::Post RPa -> Folder //" << fRecoFolder << "/PHOS/RecParticles/ not found!" << endl; - cout << "INFO: AliPHOSGetter::Post RPa -> Adding Folder //" << fRecoFolder << "/PHOS/RecParticles/" << endl; + Warning("PostRecParticles", "-> Folder //%s/PHOS/RecParticles/ not found!", fRecoFolder->GetName()) ; + Info("PostRecParticles", "-> Adding Folder //%s/PHOS/RecParticles/", fRecoFolder->GetName()) ; } phosRPaFolder = phosFolder->AddFolder("RecParticles", "RecParticles from PHOS") ; } @@ -1003,44 +1262,42 @@ Bool_t AliPHOSGetter::PostRecParticles(const char * name) const //____________________________________________________________________________ TObject** AliPHOSGetter::RecParticlesRef(const char * name) const -{ // ---------------TrackSegments ----------------------------------- +{ // ---------------RecParticles ----------------------------------- // the hierarchy is //Folders/Run/Event/RecData/PHOS/TrackSegments/name if ( !fRecoFolder ) { - cerr << "ERROR: AliPHOSGetter::RecParticlesRef -> Folder//" << fRecoFolder << " not found!" << endl; - return 0 ; + Fatal("RecParticlesRef", "Folder//%s not found !", fRecoFolder->GetName() ) ; } TFolder * phosFolder = dynamic_cast(fRecoFolder->FindObject("PHOS/RecParticles")) ; if ( !phosFolder ) { - cerr << "ERROR: AliPHOSGetter::RecParticlesRef -> Folder //" << fRecoFolder << "/PHOS/RecParticles/ not found!" << endl; - return 0; + Fatal("RecParticlesRef", "Folder //%s/PHOS/RecParticles/ not found !", fRecoFolder->GetName() ) ; } TObject * tss = phosFolder->FindObject(name ) ; if (!tss) { - return 0 ; + Fatal("RecParticlesRef", "object %s not found !", name) ; } return phosFolder->GetListOfFolders()->GetObjectRef(tss) ; } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostPID(AliPHOSPID * pid) const +const Bool_t AliPHOSGetter::PostPID(AliPHOSPID * pid) const { // ------------AliPHOS PID ----------------------------- TTask * tasks = dynamic_cast(fTasksFolder->FindObject("Reconstructioner")) ; if ( !tasks ) { - cerr << "ERROR: AliPHOSGetter::Post Per -> Task //" << fTasksFolder << "/Reconstructioner not found!" << endl; + Error("PostPID", "Task //%s/Reconstructioner not found !", fTasksFolder) ; return kFALSE ; } TTask * phos = dynamic_cast(tasks->GetListOfTasks()->FindObject("PHOS")) ; if ( !phos ) { if (fDebug) { - cout <<"WARNING: AliPHOSGetter::Post Per -> //" << fTasksFolder << "/Reconstructioner/PHOS not found!" << endl; - cout <<"INFO: AliPHOSGetter::Post Per -> Adding //" << fTasksFolder << "/Reconstructioner/PHOS" << endl; + Warning("PostPID", "//%s/Reconstructioner/PHOS not found!", fTasksFolder) ; + Info("PostPID", "Adding //%s/Reconstructioner/PHOS", fTasksFolder) ; } phos = new TTask("PHOS", "") ; tasks->Add(phos) ; @@ -1049,8 +1306,7 @@ Bool_t AliPHOSGetter::PostPID(AliPHOSPID * pid) const AliPHOSPID * phospid = dynamic_cast(phos->GetListOfTasks()->FindObject(pid->GetName())) ; if (phospid) { if (fDebug) - cout << "INFO: AliPHOSGetter::Post Per -> Task " << pid->GetName() - << " already exists" << endl ; + Info("PostPID", "-> Task %s qlready exists", pid->GetName()) ; phos->GetListOfTasks()->Remove(phospid) ; } @@ -1059,22 +1315,22 @@ Bool_t AliPHOSGetter::PostPID(AliPHOSPID * pid) const } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostPID(const char * name) const +const Bool_t AliPHOSGetter::PostPID(const char * name) const { // the hierarchy is //Folders/Tasks/Reconstructioner/PHOS/sdigitsname TTask * tasks = dynamic_cast(fTasksFolder->FindObject("Reconstructioner")) ; if ( !tasks ) { - cerr << "ERROR: AliPHOSGetter::Post Per -> Task //" << fTasksFolder << "/Reconstructioner not found!" << endl; + Error("PostPID", "Task //%s/Reconstructioner not found !", fTasksFolder->GetName() ) ; return kFALSE ; } TTask * phos = dynamic_cast(tasks->GetListOfTasks()->FindObject("PHOS")) ; if ( !phos ) { if (fDebug) { - cout <<"WARNING: AliPHOSGetter::Post Per -> //" << fTasksFolder << "/Reconstructioner/PHOS not found!" << endl; - cout <<"INFO: AliPHOSGetter::Post Per -> Adding //" << fTasksFolder << "/Reconstructioner/PHOS" << endl; + Warning("PostPID", "//%s/Reconstructioner/PHOS not found!", fTasksFolder->GetName()) ; + Info("PostPID", "Adding //%s/Reconstructioner/PHOS", fTasksFolder->GetName()) ; } phos = new TTask("PHOS", "") ; tasks->Add(phos) ; @@ -1083,7 +1339,7 @@ Bool_t AliPHOSGetter::PostPID(const char * name) const TList * l = phos->GetListOfTasks() ; TIter it(l) ; TString pidname(name) ; - pidname+=":pid" ; + pidname+=":pid" ; TTask * task ; while((task = static_cast(it.Next()) )){ TString taskname(task->GetName()) ; @@ -1093,7 +1349,8 @@ Bool_t AliPHOSGetter::PostPID(const char * name) const AliPHOSPIDv1 * phospid = new AliPHOSPIDv1() ; pidname+="-v1" ; - phospid->SetName(pidname) ; + phospid->SetName(pidname) ; + phospid->SetTitle(fHeaderFile) ; phos->Add(phospid) ; return kTRUE; @@ -1106,14 +1363,12 @@ TObject** AliPHOSGetter::PIDRef(const char * name) const TTask * tasks = dynamic_cast(fTasksFolder->FindObject("Reconstructioner")) ; if ( !tasks ) { - cerr << "ERROR: AliPHOSGetter::Post PerRef -> Task //" << fTasksFolder << "/Reconstructioner not found!" << endl; - return kFALSE ; + Fatal("PIDRef", "Task //%s/Reconstructioner not found !", fTasksFolder->GetName() ) ; } TTask * phos = dynamic_cast(tasks->GetListOfTasks()->FindObject("PHOS")) ; if ( !phos ) { - cerr <<"WARNING: AliPHOSGetter::Post PerRef -> //" << fTasksFolder << "/ReconstructionerPHOS not found!" << endl; - return 0 ; + Fatal("PIDRef", "//%s/Reconstructioner/PHOS not found !", fTasksFolder->GetName() ) ; } TList * l = phos->GetListOfTasks() ; @@ -1121,7 +1376,7 @@ TObject** AliPHOSGetter::PIDRef(const char * name) const TTask * task ; TTask * pid = 0 ; TString pidname(name) ; - pidname+=":pid-" ; + pidname+=":pid" ; while((task = static_cast(it.Next()) )){ TString taskname(task->GetName()) ; if(taskname.BeginsWith(pidname)){ @@ -1130,15 +1385,15 @@ TObject** AliPHOSGetter::PIDRef(const char * name) const } } - if(pid) - return l->GetObjectRef(pid) ; - else - return 0 ; + if(!pid) { + Fatal("PIDRef", "Task //%s/Reconstructioner/PHOS/PID/%s not found !", fTasksFolder->GetName(), name) ; + } + return l->GetObjectRef(pid) ; } //____________________________________________________________________________ -Bool_t AliPHOSGetter::PostQA(void) const +const Bool_t AliPHOSGetter::PostQA(void) const { // ------------------ QA --------------------------------- // the hierarchy is //Folders/Run/Conditions/QA/PHOS/alarmsName @@ -1146,8 +1401,8 @@ Bool_t AliPHOSGetter::PostQA(void) const TFolder * phosFolder = dynamic_cast(fQAFolder->FindObject("PHOS")) ; if ( !phosFolder ) { if (fDebug) { - cout << "WARNING: AliPHOSGetter::Post Q -> Folder //" << fQAFolder << "/PHOS/ not found!" << endl; - cout << "INFO: AliPHOSGetter::Post Q -> Adding Folder //" << fQAFolder << "/PHOS/" << endl; + Warning("PostQA", "-> Folder //%s/PHOS/ not found!", fQAFolder) ; + Info("PostQA", "-> Adding Folder //%s/PHOS", fQAFolder) ; } phosFolder = fQAFolder->AddFolder("PHOS", "QA from PHOS") ; } @@ -1162,46 +1417,132 @@ TObject** AliPHOSGetter::AlarmsRef(void) const // the hierarchy is //Folders/Run/Conditions/QA/PHOS if ( !fQAFolder ) { - cerr << "ERROR: AliPHOSGetter::Post QRef -> Folder //" << fQAFolder << " not found!" << endl; - return 0; + Fatal("AlarmsRef", "Folder //%s not found !", fQAFolder) ; } TFolder * phosFolder = dynamic_cast(fQAFolder->FindObject("PHOS")) ; if ( !phosFolder ) { - cerr << "ERROR: AliPHOSGetter::Post QRef -> Folder //" << fQAFolder << "/PHOS/ not found!" << endl; - return 0; + Fatal("AlarmsRef", "Folder //%s/PHOS/ not found !", fQAFolder) ; } return fQAFolder->GetListOfFolders()->GetObjectRef(phosFolder) ; } + +//____________________________________________________________________________ +TTree * AliPHOSGetter::TreeK(TString filename) +{ + + // returns TreeK from file filename + // usefull in case of split file + + if ( filename.IsNull() ) + filename = fHeaderFile ; + + TFile * file = 0 ; + file = static_cast(gROOT->GetFile(filename.Data() ) ) ; + if (file && (filename != fHeaderFile) ) { // file already open + file->Close() ; + delete fAlice ; + } + file = TFile::Open(filename.Data(), "read") ; + fAlice = static_cast(file->Get("gAlice")) ; + TString treeName("TreeK") ; + treeName += EventNumber() ; + TTree * tree = static_cast(file->Get(treeName.Data())) ; + if (!tree && fDebug) + Warning("TreeK", "-> %s not found in %s", treeName.Data(), filename.Data()) ; + + return tree ; +} + //____________________________________________________________________________ -const TParticle * AliPHOSGetter::Primary(Int_t index) const +TTree * AliPHOSGetter::TreeH(TString filename) { - // Return primary particle numbered by - if(index < 0) - return 0 ; - TParticle * p = gAlice->Particle(index) ; - if (p->GetFirstMother() != -1 ) { - cout << "AliPHOSGetter::Primary : Not a primary " << endl ; + // returns TreeH from file filename + // usefull in case of split file + + if ( filename.IsNull() ) + filename = fHeaderFile ; + + TFile * file = 0 ; + file = static_cast(gROOT->GetFile(filename.Data() ) ) ; + if (!file) { // file not open yet + file = TFile::Open(filename.Data(), "read") ; } - return p ; + TString treeName("TreeH") ; + treeName += EventNumber() ; + TTree * tree = static_cast(file->Get(treeName.Data())) ; + if (!tree && fDebug) + Warning("TreeH", "-> %s not found in %s", treeName.Data(), filename.Data()) ; + return tree ; +} + +//____________________________________________________________________________ +TTree * AliPHOSGetter::TreeS(TString filename) +{ + + // returns TreeS from file filename + // usefull in case of split file + + if ( filename.IsNull() ) + filename = fHeaderFile ; + + TFile * file = 0 ; + file = static_cast(gROOT->GetFile(filename.Data() ) ) ; + if (!file) { // file not open yet + file = TFile::Open(filename.Data(), "read") ; + } + TString treeName("TreeS") ; + treeName += EventNumber() ; + TTree * tree = static_cast(file->Get(treeName.Data())) ; + if (!tree && fDebug) + Warning("TreeS", "-> %s not found in %s", treeName.Data(), filename.Data() ); -// Int_t primaryIndex = index % 10000000 ; -// Int_t primaryList = (Int_t ) ((index-primaryIndex)/10000000.) ; - -// if ( primaryList > 0 ) { -// if (fDebug) { -// cout << " Getter does not support currently Mixing of primary " << endl ; -// cout << " can not return primary: " << index<< " (list "<< primaryList<< " primary # " << primaryIndex << " )"<(gROOT->GetFile(filename.Data() ) ) ; + if (!file) { // file not open yet + file = TFile::Open(filename.Data(), "read") ; + } + TString treeName("TreeD") ; + treeName += EventNumber() ; + TTree * tree = static_cast(file->Get(treeName.Data())) ; + if (!tree && fDebug) + Warning("TreeD", "-> %s not found in %s", treeName.Data(), filename.Data()) ; -// return gAlice->Particle(primaryIndex) ; + return tree ; +} + +//____________________________________________________________________________ +const TParticle * AliPHOSGetter::Primary(Int_t index) const +{ + // Return primary particle numbered by + + if(index < 0) + return 0 ; + TParticle * p = 0 ; + if (fAlice) + p = fAlice->Particle(index) ; + else + p = gAlice->Particle(index) ; + return p ; + } //____________________________________________________________________________ @@ -1223,15 +1564,29 @@ const TParticle * AliPHOSGetter::Secondary(TParticle* p, Int_t index) const } //____________________________________________________________________________ -void AliPHOSGetter::ReadTreeD() +Int_t AliPHOSGetter::ReadTreeD(const Int_t event) { // Read the digit tree gAlice->TreeD() - if(gAlice->TreeD()== 0){ - cerr << "ERROR: AliPHOSGetter::ReadTreeD: can not read TreeD " << endl ; - return ; + + TTree * treeD ; + if(fToSplit){ + TFile * file = static_cast(gROOT->GetFile(fDigitsFileName)); + if(!file) + file = TFile::Open(fDigitsFileName) ; + // Get Digits Tree header from file + TString treeName("TreeD") ; + treeName += event ; + treeD = dynamic_cast(file->Get(treeName.Data())); + if(!treeD){ // TreeD not found in header file + if (fDebug) + Warning("ReadTreeD", "-> Cannot find TreeD in %s", fDigitsFileName.Data()) ; + return 1; + } } + else + treeD = gAlice->TreeD() ; - TObjArray * lob = static_cast(gAlice->TreeD()->GetListOfBranches()) ; + TObjArray * lob = static_cast(treeD->GetListOfBranches()) ; TIter next(lob) ; TBranch * branch = 0 ; TBranch * digitsbranch = 0 ; @@ -1243,19 +1598,20 @@ void AliPHOSGetter::ReadTreeD() digitsbranch = branch ; phosfound = kTRUE ; } - else if ( (strcmp(branch->GetName(), "AliPHOSDigitizer")==0) && (strcmp(branch->GetTitle(), fDigitsTitle)==0) ) { + else if ( ((strcmp(branch->GetName(), "AliPHOSDigitizer")==0)|| + (strcmp(branch->GetName(), "AliPHOSRaw2Digits")==0)) && + (strcmp(branch->GetTitle(), fDigitsTitle)==0) ) { digitizerbranch = branch ; digitizerfound = kTRUE ; } } - + if ( !phosfound || !digitizerfound ) { if (fDebug) - cout << "WARNING: AliPHOSGetter::ReadTreeD -> Cannot find Digits and/or Digitizer with name " - << fDigitsTitle << endl ; - return ; + Warning("ReadTreeD", "-> Cannot find Digits and/or Digitizer with name %s", fDigitsTitle.Data()) ; + return 2; } - + //read digits if(!Digits(fDigitsTitle) ) PostDigits(fDigitsTitle); @@ -1264,34 +1620,95 @@ void AliPHOSGetter::ReadTreeD() // read the Digitizer - if(!Digitizer(fDigitsTitle)) - PostDigitizer(fDigitsTitle) ; + if(Digitizer()){ + if(strcmp(Digitizer()->IsA()->GetName(),digitizerbranch->GetName())!=0){ + RemoveTask("D", fDigitsTitle) ; + if(strcmp(digitizerbranch->GetName(), "AliPHOSDigitizer")==0) + PostDigitizer("Digitizer") ; + else + PostDigitizer("Raw2Digits") ; + } + } + else{ + if(strcmp(digitizerbranch->GetName(), "AliPHOSDigitizer")==0) + PostDigitizer("Digitizer") ; + else + PostDigitizer("Raw2Digits") ; + } + + digitizerbranch->SetAddress(DigitizerRef(fDigitsTitle)) ; digitizerbranch->GetEntry(0) ; - -} + if((!fcdb)&&(strcmp(digitizerbranch->GetName(), "AliPHOSRaw2Digits")==0)) + ReadCalibrationDB("Primordial","beamtest.root") ; + + if(gAlice->TreeD()!=treeD) + treeD->Delete(); + + return 0 ; +} //____________________________________________________________________________ -void AliPHOSGetter::ReadTreeH() -{ - // Read the first entry of PHOS branch in hit tree gAlice->TreeH() +void AliPHOSGetter::ReadCalibrationDB(const char * database,const char * filename){ - if(gAlice->TreeH()== 0){ - cerr << "ERROR: AliPHOSGetter::ReadTreeH: -> Cannot read TreeH " << endl ; + if(fcdb && (strcmp(database,fcdb->GetTitle())==0)) + return ; + + TFile * file = gROOT->GetFile(filename) ; + if(!file) + file = TFile::Open(filename); + if(!file){ + Error ("ReadCalibrationDB", "Cannot open file %s", filename) ; return ; } + if(fcdb) + fcdb->Delete() ; + fcdb = dynamic_cast(file->Get("AliPHOSCalibrationDB")) ; + if(!fcdb) + Error ("ReadCalibrationDB", "No database %s in file %s", database, filename) ; +} + +//____________________________________________________________________________ +Int_t AliPHOSGetter::ReadTreeH() +{ + // Read the first entry of PHOS branch in hit tree gAlice->TreeH() + + TTree * treeH = gAlice->TreeH() ; + + if(!treeH) {// TreeH not found in header file + + if (fDebug) + Warning("ReadTreeH", "-> Cannot find TreeH in %s", fHeaderFile.Data() ) ; + + TString searchFileName("PHOS.Hits") ; + if((strcmp(fBranchTitle.Data(),"Default")!=0)&&(strcmp(fBranchTitle.Data(),"")!=0)){ + searchFileName+="." ; + searchFileName += fBranchTitle ; + } + searchFileName+=".root" ; + + if ( (treeH = TreeH(searchFileName)) ) { //found TreeH in the file which contains the hits + if (fDebug) + Info("ReadTreeH", "-> TreeH found in %s", searchFileName.Data()) ; + + } else { + Error("ReadTreeH", "TreeH not found") ; + return 1; + } + } - TBranch * hitsbranch = static_cast(gAlice->TreeH()->GetBranch("PHOS")) ; + TBranch * hitsbranch = static_cast(treeH->GetBranch("PHOS")) ; if ( !hitsbranch ) { if (fDebug) - cout << "WARNING: AliPHOSGetter::ReadTreeH -> Cannot find branch PHOS" << endl ; - return ; + Warning("ReadTreeH", "-> Cannot find branch PHOS") ; + return 2; } if(!Hits()) PostHits() ; if (hitsbranch->GetEntries() > 1 ) { + (dynamic_cast (*HitsRef()))->Clear() ; TClonesArray * tempo = new TClonesArray("AliPHOSHit",1000) ; TClonesArray * hits = dynamic_cast(*HitsRef()) ; hitsbranch->SetAddress(&tempo) ; @@ -1309,35 +1726,38 @@ void AliPHOSGetter::ReadTreeH() delete tempo ; } else { + (dynamic_cast (*HitsRef()))->Clear() ; hitsbranch->SetAddress(HitsRef()) ; hitsbranch->GetEntry(0) ; } + return 0 ; } //____________________________________________________________________________ -void AliPHOSGetter::Track(Int_t itrack) +void AliPHOSGetter::Track(const Int_t itrack) { // Read the first entry of PHOS branch in hit tree gAlice->TreeH() if(gAlice->TreeH()== 0){ - cerr << "ERROR: AliPHOSGetter::ReadTreeH: -> Cannot read TreeH " << endl ; + Error("Track", "Cannot read TreeH") ; return ; } TBranch * hitsbranch = dynamic_cast(gAlice->TreeH()->GetListOfBranches()->FindObject("PHOS")) ; if ( !hitsbranch ) { if (fDebug) - cout << "WARNING: AliPHOSGetter::ReadTreeH -> Cannot find branch PHOS" << endl ; + Warning("Track", "Cannot find branch PHOS") ; return ; } if(!Hits()) PostHits() ; + (dynamic_cast (*HitsRef()))->Clear() ; hitsbranch->SetAddress(HitsRef()) ; hitsbranch->GetEntry(itrack) ; - } + //____________________________________________________________________________ void AliPHOSGetter::ReadTreeQA() { @@ -1345,14 +1765,14 @@ void AliPHOSGetter::ReadTreeQA() // so far only PHOS knows about this Tree if(PHOS()->TreeQA()== 0){ - cerr << "ERROR: AliPHOSGetter::ReadTreeQA: can not read TreeQA " << endl ; + Error("ReadTreeQA", "Cannot read TreeQA") ; return ; } TBranch * qabranch = PHOS()->TreeQA()->GetBranch("PHOS") ; if (!qabranch) { if (fDebug) - cout << "WARNING: AliPHOSGetter::ReadTreeQA -> Cannot find QA Alarms for PHOS" << endl ; + Warning("ReadTreeQA", "Cannot find QA Alarms for PHOS"); return ; } @@ -1372,26 +1792,49 @@ void AliPHOSGetter::ReadTreeQA() } //____________________________________________________________________________ -void AliPHOSGetter::ReadTreeR() +Int_t AliPHOSGetter::ReadTreeR(const Int_t event) { // Read the reconstrunction tree gAlice->TreeR() - - if(gAlice->TreeR()== 0){ - cerr << "ERROR: AliPHOSGetter::ReadTreeR: can not read TreeR " << endl ; - return ; + // A particularity has been introduced here : + // if gime->Event(ievent,"R") is called branches with the current title are read, the current title + // being for example give in AliPHOSPID(fileName, title) + // if gime(Event(ievent, "RA") is called the title of the branches is not checked anymore, "A" stands for any + // This is a feature needed by PID to be able to reconstruct several times particles (each time a ther title is given) + // from a given set of TrackSegments (with a given name) + // This is why any is NOT used to read the branch of RecParticles + // any migh have become obsolete : to be checked + // See AliPHOSPIDv1 + + TTree * treeR ; + if(fToSplit){ + TFile * file = static_cast(gROOT->GetFile(fRecPointsFileName)); + if(!file) + file = TFile::Open(fRecPointsFileName) ; + // Get Digits Tree header from file + TString treeName("TreeR") ; + treeName += event ; + treeR = dynamic_cast(file->Get(treeName.Data())); + if(!treeR){ // TreeR not found in header file + if (fDebug) + Warning("ReadTreeD", "-> Cannot find TreeR in %s", fRecPointsFileName.Data()) ; + return 1; + } } + else + treeR = gAlice->TreeR() ; // RecPoints - TObjArray * lob = static_cast(gAlice->TreeR()->GetListOfBranches()) ; + TObjArray * lob = static_cast(treeR->GetListOfBranches()) ; TIter next(lob) ; TBranch * branch = 0 ; TBranch * emcbranch = 0 ; TBranch * cpvbranch = 0 ; TBranch * clusterizerbranch = 0 ; Bool_t phosemcrpfound = kFALSE, phoscpvrpfound = kFALSE, clusterizerfound = kFALSE ; + - while ( (branch = static_cast(next())) && (!phosemcrpfound || !phoscpvrpfound || !clusterizerfound) ) - if(strcmp(branch->GetTitle(), fRecPointsTitle)==0) { + while ( (branch = static_cast(next())) && (!phosemcrpfound || !phoscpvrpfound || !clusterizerfound) ) { + if(strcmp(branch->GetTitle(), fRecPointsTitle)==0 ) { if ( strcmp(branch->GetName(), "PHOSEmcRP")==0) { emcbranch = branch ; phosemcrpfound = kTRUE ; @@ -1405,51 +1848,36 @@ void AliPHOSGetter::ReadTreeR() clusterizerfound = kTRUE ; } } + } - if ( !phosemcrpfound ) { + if ( !phosemcrpfound || !phoscpvrpfound || !clusterizerfound) { if (fDebug) - cout << "WARNING: AliPHOSGetter::ReadTreeR -> Cannot find EmcRecPoints with title " - << fRecPointsTitle << endl ; - return ; - } - if ( !phoscpvrpfound ) { - if (fDebug) - cout << "WARNING: AliPHOSGetter::ReadTreeR -> Cannot find CpvRecPoints with title " - << fRecPointsTitle << endl ; - return ; - } - if ( !clusterizerfound ) { - if (fDebug) - cout << "WARNING: AliPHOSGetter::ReadTreeR -> Can not find Clusterizer with title " - << fRecPointsTitle << endl ; - return ; - } - - // Read and Post the RecPoints - if(!EmcRecPoints(fRecPointsTitle) ) - PostRecPoints(fRecPointsTitle) ; - - emcbranch->SetAddress(EmcRecPointsRef(fRecPointsTitle)) ; - emcbranch->GetEntry(0) ; - - cpvbranch->SetAddress(CpvRecPointsRef(fRecPointsTitle)) ; - cpvbranch->GetEntry(0) ; - - if(!Clusterizer(fRecPointsTitle) ) - PostClusterizer(fRecPointsTitle) ; - - clusterizerbranch->SetAddress(ClusterizerRef(fRecPointsTitle)) ; - clusterizerbranch->GetEntry(0) ; + Warning("ReadTreeR", "-> Cannot find RecPoints and/or Clusterizer with name %s", fRecPointsTitle.Data() ) ; + } else { + if(!EmcRecPoints(fRecPointsTitle) ) + PostRecPoints(fRecPointsTitle) ; + + emcbranch->SetAddress(EmcRecPointsRef(fRecPointsTitle)) ; + emcbranch->GetEntry(0) ; + + cpvbranch->SetAddress(CpvRecPointsRef(fRecPointsTitle)) ; + cpvbranch->GetEntry(0) ; + + if(!Clusterizer(fRecPointsTitle) ) + PostClusterizer(fRecPointsTitle) ; + + clusterizerbranch->SetAddress(ClusterizerRef(fRecPointsTitle)) ; + clusterizerbranch->GetEntry(0) ; + } //------------------- TrackSegments --------------------- next.Reset() ; TBranch * tsbranch = 0 ; TBranch * tsmakerbranch = 0 ; Bool_t phostsfound = kFALSE, tsmakerfound = kFALSE ; - - while ( (branch = static_cast(next())) && (!phostsfound || !tsmakerfound) ) - if(strcmp(branch->GetTitle(), fTrackSegmentsTitle)==0) { + while ( (branch = static_cast(next())) && (!phostsfound || !tsmakerfound) ) { + if(strcmp(branch->GetTitle(), fTrackSegmentsTitle)==0 ) { if ( strcmp(branch->GetName(), "PHOSTS")==0){ tsbranch = branch ; phostsfound = kTRUE ; @@ -1459,27 +1887,24 @@ void AliPHOSGetter::ReadTreeR() tsmakerfound = kTRUE ; } } - + } + if ( !phostsfound || !tsmakerfound ) { if (fDebug) - cout << "WARNING: AliPHOSGetter::ReadTreeR -> Cannot find TrackSegments and/or TrackSegmentMaker with name " - << fTrackSegmentsTitle << endl ; - return ; - } - - // Read and Post the TrackSegments - if(!TrackSegments(fTrackSegmentsTitle)) - PostTrackSegments(fTrackSegmentsTitle) ; - - tsbranch->SetAddress(TrackSegmentsRef(fTrackSegmentsTitle)) ; - tsbranch->GetEntry(0) ; - - // Read and Post the TrackSegment Maker - if(!TrackSegmentMaker(fTrackSegmentsTitle)) - PostTrackSegmentMaker(fTrackSegmentsTitle) ; - - tsmakerbranch->SetAddress(TSMakerRef(fTrackSegmentsTitle)) ; - tsmakerbranch->GetEntry(0) ; + Warning("ReadTreeR", "-> Cannot find TrackSegments and/or TrackSegmentMaker with name %s", fTrackSegmentsTitle.Data() ) ; + } else { + // Read and Post the TrackSegments + if(!TrackSegments(fTrackSegmentsTitle)) + PostTrackSegments(fTrackSegmentsTitle) ; + tsbranch->SetAddress(TrackSegmentsRef(fTrackSegmentsTitle)) ; + tsbranch->GetEntry(0) ; + + // Read and Post the TrackSegment Maker + if(!TrackSegmentMaker(fTrackSegmentsTitle)) + PostTrackSegmentMaker(fTrackSegmentsTitle) ; + tsmakerbranch->SetAddress(TSMakerRef(fTrackSegmentsTitle)) ; + tsmakerbranch->GetEntry(0) ; + } //------------ RecParticles ---------------------------- @@ -1502,63 +1927,60 @@ void AliPHOSGetter::ReadTreeR() if ( !phosrpafound || !pidfound ) { if (fDebug) - cout << "WARNING: AliPHOSGetter::ReadTreeR -> Cannot find RecParticles and/or PID with name " - << fRecParticlesTitle << endl ; - return ; - } - - // Read and Post the RecParticles - if(!RecParticles(fRecParticlesTitle)) - PostRecParticles(fRecParticlesTitle) ; - - rpabranch->SetAddress(RecParticlesRef(fRecParticlesTitle)) ; - rpabranch->GetEntry(0) ; - - // Read and Post the PID - if(!PID(fRecParticlesTitle)) - PostPID(fRecParticlesTitle) ; - - pidbranch->SetAddress(PIDRef(fRecParticlesTitle)) ; - pidbranch->GetEntry(0) ; - - + Warning("ReadTreeR", "-> Cannot find RecParticles and/or PID with name %s", fRecParticlesTitle.Data() ) ; + } else { + // Read and Post the RecParticles + if(!RecParticles(fRecParticlesTitle)) + PostRecParticles(fRecParticlesTitle) ; + rpabranch->SetAddress(RecParticlesRef(fRecParticlesTitle)) ; + rpabranch->GetEntry(0) ; + // Read and Post the PID + if(!PID(fRecParticlesTitle)) + PostPID(fRecParticlesTitle) ; + pidbranch->SetAddress(PIDRef(fRecParticlesTitle)) ; + pidbranch->GetEntry(0) ; + } + + if(gAlice->TreeR()!=treeR) + treeR->Delete(); + return 0 ; } //____________________________________________________________________________ -void AliPHOSGetter::ReadTreeS(Int_t event) +Int_t AliPHOSGetter::ReadTreeS(const Int_t event) { - // Read the summable digits tree gAlice->TreeS() - - // loop over all opened files and read their SDigits to the White Board + // Reads the SDigits treeS from all files + // Files, which should be opened are listed in phosF + // So, first get list of files TFolder * phosF = dynamic_cast(fSDigitsFolder->FindObject("PHOS")) ; if (!phosF) phosF = fSDigitsFolder->AddFolder("PHOS", "SDigits from PHOS") ; TCollection * folderslist = phosF->GetListOfFolders() ; - //Add current file to list if it is not there yet - if ( (fHeaderFile != "aliroot") && ( !folderslist->Contains(fHeaderFile) ) ){ - phosF->AddFolder(fHeaderFile, ""); - } - + // Now iterate over the list of files and read TreeS into Whiteboard TIter next(folderslist) ; TFolder * folder = 0 ; TFile * file; TTree * treeS = 0; while ( (folder = static_cast(next())) ) { - if(fHeaderFile.CompareTo(folder->GetName()) == 0 ) - treeS=gAlice->TreeS() ; - else{ - file = static_cast(gROOT->GetFile(folder->GetName())); - file->cd() ; - - // Get SDigits Tree header from file - TString treeName("TreeS") ; - treeName += event ; - treeS = dynamic_cast(gDirectory->Get(treeName.Data())); - } - if(treeS==0){ - cerr << "ERROR: AliPHOSGetter::ReadTreeS There is no SDigit Tree" << endl; - return ; + TString fileName("") ; + if(fToSplit) + fileName = folder->GetTitle() ; + else + fileName = folder->GetName() ; + fileName.ReplaceAll("_","/") ; + file = static_cast(gROOT->GetFile(fileName)); + if(!file) + file = TFile::Open(fileName) ; + // Get SDigits Tree header from file + TString treeName("TreeS") ; + treeName += event ; + treeS = dynamic_cast(file->Get(treeName.Data())); + + if(!treeS){ // TreeS not found in header file + if (fDebug) + Warning("ReadTreeS", "-> Cannot find TreeS in %s", fileName.Data()) ; + return 1; } //set address of the SDigits and SDigitizer @@ -1568,28 +1990,29 @@ void AliPHOSGetter::ReadTreeS(Int_t event) TObjArray * lob = static_cast(treeS->GetListOfBranches()) ; TIter next(lob) ; Bool_t phosfound = kFALSE, sdigitizerfound = kFALSE ; - + while ( (branch = static_cast(next())) && (!phosfound || !sdigitizerfound) ) { if ( (strcmp(branch->GetName(), "PHOS")==0) && (strcmp(branch->GetTitle(), fSDigitsTitle)==0) ) { phosfound = kTRUE ; sdigitsBranch = branch ; } - else if ( (strcmp(branch->GetName(), "AliPHOSSDigitizer")==0) && (strcmp(branch->GetTitle(), fSDigitsTitle)==0) ) { + else if ( (strcmp(branch->GetName(), "AliPHOSSDigitizer")==0) && + (strcmp(branch->GetTitle(), fSDigitsTitle)==0) ) { sdigitizerfound = kTRUE ; sdigitizerBranch = branch ; } } if ( !phosfound || !sdigitizerfound ) { if (fDebug) - cout << "WARNING: AliPHOSDigitizer::ReadSDigits -> Digits and/or Digitizer branch with name " << GetName() - << " not found" << endl ; - return ; + Warning("ReadSDigits", "-> Digits and/or Digitizer branch with name %s not found", GetName()) ; + return 2; } if ( !folder->FindObject(fSDigitsTitle) ) PostSDigits(fSDigitsTitle,folder->GetName()) ; + ((TClonesArray*) (*SDigitsRef(fSDigitsTitle,folder->GetName())))->Clear() ; sdigitsBranch->SetAddress(SDigitsRef(fSDigitsTitle,folder->GetName())) ; sdigitsBranch->GetEntry(0) ; @@ -1598,24 +2021,17 @@ void AliPHOSGetter::ReadTreeS(Int_t event) sdname+=folder->GetName() ; if(!SDigitizer(sdname) ) PostSDigitizer(fSDigitsTitle,folder->GetName()) ; - sdigitizerBranch->SetAddress(SDigitizerRef(sdname)) ; - sdigitizerBranch->GetEntry(0) ; - + sdigitizerBranch->GetEntry(0) ; + if(gAlice->TreeS()!=treeS) + treeS->Delete(); } - - // After SDigits have been read from all files, return to the first one - - next.Reset(); - folder = static_cast(next()); - if(folder){ - file = static_cast(gROOT->GetFile(folder->GetName())); - file ->cd() ; - } - + return 0 ; } + //____________________________________________________________________________ void AliPHOSGetter::ReadTreeS(TTree * treeS, Int_t input) + { // Read the summable digits fron treeS() @@ -1648,7 +2064,7 @@ void AliPHOSGetter::ReadTreeS(TTree * treeS, Int_t input) } if ( !phosfound || !sdigitizerfound ) { if (fDebug) - cout << "WARNING: AliPHOSGetter::ReadTreeS -> Digits and/or Digitizer branch not found" << endl ; + Warning("ReadTreeS", "-> Digits and/or Digitizer branch not found") ; return ; } @@ -1661,151 +2077,80 @@ void AliPHOSGetter::ReadTreeS(TTree * treeS, Int_t input) TString sdname(sdigitsBranch->GetTitle()) ; sdname+=":" ; sdname+=filename ; + if(!SDigitizer(sdigitsBranch->GetTitle()) ) PostSDigitizer(sdigitsBranch->GetTitle(),filename) ; - sdigitizerBranch->SetAddress(SDigitizerRef(sdname)) ; sdigitizerBranch->GetEntry(0) ; - + if(gAlice->TreeS()!=treeS) + treeS->Delete(); } //____________________________________________________________________________ void AliPHOSGetter::ReadPrimaries() { + // a lot simplified.... if 2 files are opened then we have a problem - - // Read the first entry of PHOS branch in hit tree gAlice->TreeH() - TTree * particleTree = gAlice->TreeK() ; - if( ! particleTree ){ - cerr << "ERROR: AliPHOSGetter::ReadTreePrimaries: -> Cannot read TreeK " << endl ; - return ; - } - - TBranch * particlesbranch = static_cast(particleTree->GetBranch("Particles")) ; - if ( !particlesbranch ) { - if (fDebug) - cout << "WARNING: AliPHOSGetter::ReadTreePrimaries -> Cannot find branch Particles" << endl ; - return ; - } - - TParticle * particle = 0; - TClonesArray * ar ; + TClonesArray * ar = 0 ; if(! (ar = Primaries()) ) { PostPrimaries() ; ar = Primaries() ; } ar->Delete() ; - particlesbranch->SetAddress(&particle) ; - Int_t index = 0 ; - for (index = 0 ; index < particleTree->GetEntries(); index++) { - particlesbranch->GetEntry(index) ; - new ((*ar)[index]) TParticle(*particle); - } - fNPrimaries= ar->GetEntries() ; + if (TreeK(fHeaderFile)) { // treeK found in header file + if (fDebug) + Info("ReadPrimaries", "-> TreeK found in %s", fHeaderFile.Data() ); + fNPrimaries = gAlice->GetNtrack() ; + fAlice = 0 ; - // Reads specific branches of primaries - // fNPrimaries = gAlice->GetNtrack(); - - //Check, is it necessary to open new files - // TArrayI* events = fDigitizer->GetCurrentEvents() ; -// TClonesArray * filenames = fDigitizer->GetHeadersFiles() ; -// Int_t input ; -// for(input = 0; input < filenames->GetEntriesFast(); input++){ - -// TObjString * filename = (TObjString *) filenames->At(input) ; - -// //Test, if this file already open -// TFile *file = (TFile*) gROOT->GetFile( filename->GetString() ) ; -// if(file == 0) -// file = new TFile( filename->GetString()) ; -// file->cd() ; - -// // Get Kine Tree from file -// char treeName[20]; -// sprintf(treeName,"TreeK%d",events->At(input)); -// TTree * treeK = (TTree*)gDirectory->Get(treeName); -// if (treeK) -// treeK->SetBranchAddress("Particles", &fParticleBuffer); -// else -// cout << "AliPHOSGetter: cannot find Kine Tree for event:" << events->At(input) << endl; - -// // Create the particle stack -// if(!fParticles) fParticles = new TClonesArray("TParticle",1000); -// // Build the pointer list -// if(fParticleMap) { <---- -// fParticleMap->Clear(); -// fParticleMap->Expand(treeK->GetEntries()); -// } else -// fParticleMap = new TObjArray(treeK->GetEntries()); + } else { // treeK not found in header file -// // From gAlice->Particle(i) - - -// if(!(*fParticleMap)[i]) { -// Int_t nentries = fParticles->GetEntries(); - -// // algorithmic way of getting entry index -// // (primary particles are filled after secondaries) -// Int_t entry; -// if (iGetEntry(fParticleFileMap[i]); -// new ((*fParticles)[nentries]) TParticle(*fParticleBuffer); -// fParticleMap->AddAt((*fParticles)[nentries],i); -// } -// return (TParticle *) (*fParticleMap)[i]; - - + Error("ReadPrimaries", "TreeK not found") ; + return ; -// // } - - -// // //scan over opened files and read corresponding TreeK## - -// return ; + } + Int_t index = 0 ; + for (index = 0 ; index < fNPrimaries; index++) { + new ((*ar)[index]) TParticle(*(Primary(index))); + } } + //____________________________________________________________________________ -void AliPHOSGetter::Event(const Int_t event, const char* opt) +void AliPHOSGetter::Event(const Int_t event, const char* opt) { // Reads the content of all Tree's S, D and R - + if (event >= gAlice->TreeE()->GetEntries() ) { - cerr << "ERROR: AliPHOSGetter::Event -> " << event << " not found in TreeE!" << endl ; + Error("Event", "%d not found in TreeE !", event) ; return ; } + + Bool_t any = kFALSE ; + if (strstr(opt,"A") ) // do not check the title of the branches + any = kTRUE; + gAlice->GetEvent(event) ; - if(strstr(opt,"H") ) - ReadTreeH() ; - - if(strstr(opt,"S") ) - ReadTreeS(event) ; + if( strstr(opt,"R") ) + ReadTreeR(event) ; if( strstr(opt,"D") ) - ReadTreeD() ; + ReadTreeD(event) ; - if( strstr(opt,"R") ) - ReadTreeR() ; + if(strstr(opt,"S") ) + ReadTreeS(event) ; + if(strstr(opt,"H") ) + ReadTreeH() ; + if( strstr(opt,"Q") ) ReadTreeQA() ; - - // if( strstr(opt,"P") || (strcmp(opt,"")==0) ) - // ReadPrimaries() ; - + + if( strstr(opt,"P") || (strcmp(opt,"")==0) ) + ReadPrimaries() ; + } //____________________________________________________________________________ @@ -1820,19 +2165,21 @@ TObject * AliPHOSGetter::ReturnO(TString what, TString name, TString file) const TFolder * folder = 0 ; TObject * phosO = 0 ; - // if ( name.IsNull() ) { if ( what.CompareTo("Primaries") == 0 ) { folder = dynamic_cast(fPrimariesFolder->FindObject("Primaries")) ; if (folder) phosO = dynamic_cast(folder->FindObject("Primaries")) ; + else + return 0 ; } else if ( what.CompareTo("Hits") == 0 ) { folder = dynamic_cast(fHitsFolder->FindObject("PHOS")) ; if (folder) phosO = dynamic_cast(folder->FindObject("Hits")) ; } - else if ( what.CompareTo("SDigits") == 0 ) { - TString path = "PHOS/" + file ; + else if ( what.CompareTo("SDigits") == 0 ) { + file.ReplaceAll("/","_") ; + TString path = "PHOS/" + file ; folder = dynamic_cast(fSDigitsFolder->FindObject(path.Data())) ; if (folder) { if (name.IsNull()) @@ -1853,7 +2200,7 @@ TObject * AliPHOSGetter::ReturnO(TString what, TString name, TString file) const if (folder) { if (name.IsNull()) name = fRecPointsTitle ; - phosO = dynamic_cast(folder->FindObject(name)) ; + phosO = dynamic_cast(folder->FindObject(name)) ; } } else if ( what.CompareTo("CpvRecPoints") == 0 ) { @@ -1877,7 +2224,7 @@ TObject * AliPHOSGetter::ReturnO(TString what, TString name, TString file) const if (folder) { if (name.IsNull()) name = fRecParticlesTitle ; - phosO = dynamic_cast(folder->FindObject(name)) ; + phosO = dynamic_cast(folder->FindObject(name)) ; } } else if ( what.CompareTo("Alarms") == 0 ){ @@ -1893,7 +2240,7 @@ TObject * AliPHOSGetter::ReturnO(TString what, TString name, TString file) const } if (!phosO) { if(fDebug) - cerr << "ERROR : AliPHOSGetter::ReturnO -> Object " << what << " not found in " << folder->GetName() << endl ; + Warning("ReturnO", "Object %s not found in PHOS", what.Data() ) ; return 0 ; } @@ -1919,13 +2266,13 @@ const TTask * AliPHOSGetter::ReturnT(TString what, TString name) const TTask * tasks = dynamic_cast(fTasksFolder->FindObject(search)) ; if (!tasks) { - cerr << "ERROR: AliPHOSGetter::ReturnT -> Task " << what << " not found!" << endl ; + Error("ReturnT", "Task %s not found !", what.Data() ) ; return 0 ; } TTask * phosT = dynamic_cast(tasks->GetListOfTasks()->FindObject("PHOS")) ; if (!phosT) { - cerr << "ERROR: AliPHOSGetter::ReturnT -> Task " << what << "/PHOS not found!" << endl ; + Error("ReturnT", "Task %s/PHOS not found !", what.Data() ) ; return 0 ; } @@ -1966,6 +2313,380 @@ const TTask * AliPHOSGetter::ReturnT(TString what, TString name) const } if(fDebug) - cout << "WARNING: AliPHOSGetter::ReturnT -> Task " << search << "/" << name << " not found!" << endl ; + Warning("ReturnT", "-> Task %s/PHOS/%s not found", search.Data(), name.Data() ) ; return 0 ; } + +//____________________________________________________________________________ +void AliPHOSGetter::RemoveTask(TString opt, TString name) const +{ + // remove a task from the folder + // path is fTasksFolder/SDigitizer/PHOS/name + + TTask * task = 0 ; + TTask * phos = 0 ; + TList * lofTasks = 0 ; + + if (opt == "S") { // SDigitizer + task = dynamic_cast(fTasksFolder->FindObject("SDigitizer")) ; + if (!task) + return ; + } + else if (opt == "D") { // Digitizer + task = dynamic_cast(fTasksFolder->FindObject("Digitizer")) ; + if (!task) + return ; + } + else if (opt == "C" || opt == "T" || opt == "P" ) { // Clusterizer, TrackSegmentMaker, PID + task = dynamic_cast(fTasksFolder->FindObject("Reconstructioner")) ; + if (!task) + return ; + } + else { + Warning("RemoveTask", "Unknown option %s", opt.Data() ); + return ; + } + phos = dynamic_cast(task->GetListOfTasks()->FindObject("PHOS")) ; + if (!phos) + return ; + lofTasks = phos->GetListOfTasks() ; + if (!lofTasks) + return ; + TObject * obj = lofTasks->FindObject(name) ; + if (obj) + lofTasks->Remove(obj) ; + +} + +//____________________________________________________________________________ +void AliPHOSGetter::RemoveObjects(TString opt, TString name) const +{ + // remove SDigits from the folder + // path is fSDigitsFolder/fHeaderFileName/name + + TFolder * phos = 0 ; + TFolder * phosmain = 0 ; + + if (opt == "H") { // Hits + phos = dynamic_cast(fHitsFolder->FindObject("PHOS")) ; + if (!phos) + return ; + name = "Hits" ; + } + + else if ( opt == "S") { // SDigits + phosmain = dynamic_cast(fSDigitsFolder->FindObject("PHOS")) ; + if (!phosmain) + return ; + phos = dynamic_cast(phosmain->FindObject(fHeaderFile)) ; + if (!phos) + return ; + } + + else if (opt == "D") { // Digits + phos = dynamic_cast(fDigitsFolder->FindObject("PHOS")) ; + if (!phos) + return ; + } + + else if (opt == "RE") { // EMCARecPoints + phos = dynamic_cast(fRecoFolder->FindObject("PHOS/EMCARecPoints")) ; + if (!phos) + return ; + } + + else if (opt == "RC") { // CPVRecPoints + phos = dynamic_cast(fRecoFolder->FindObject("PHOS/CPVRecPoints")) ; + if (!phos) + return ; + } + + else if (opt == "T") { // TrackSegments + phos = dynamic_cast(fRecoFolder->FindObject("PHOS/TrackSegments")) ; + if (!phos) + return ; + } + + else if (opt == "P") { // RecParticles + phos = dynamic_cast(fRecoFolder->FindObject("PHOS/RecParticles")) ; + if (!phos) + return ; + } + + else { + Warning("RemoveObjects", "Unknown option %s", opt.Data() ) ; + return ; + } + + TObjArray * ar = dynamic_cast(phos->FindObject(name)) ; + if (ar) { + phos->Remove(ar) ; + ar->Delete() ; + delete ar ; + } + + if (opt == "S") + phosmain->Remove(phos) ; + +} + +//____________________________________________________________________________ +void AliPHOSGetter::RemoveSDigits() const +{ + TFolder * phos= dynamic_cast(fSDigitsFolder->FindObject("PHOS")) ; + if (!phos) + return ; + + phos->SetOwner() ; + phos->Clear() ; +} + +//____________________________________________________________________________ +void AliPHOSGetter::CleanWhiteBoard(void){ + + TFolder * phosmain = 0 ; + TFolder * phos ; + TObjArray * ar ; + TList * lofTasks = 0 ; + TTask * task = 0 ; + TTask * phost = 0 ; + + // Hits + phos = dynamic_cast(fHitsFolder->FindObject("PHOS")) ; + if (phos){ + TObjArray * ar = dynamic_cast(phos->FindObject("Hits")) ; + if (ar) { + phos->Remove(ar) ; + ar->Delete() ; + delete ar ; + } + } + + // SDigits + phosmain = dynamic_cast(fSDigitsFolder->FindObject("PHOS")) ; + if (phosmain){ + phos = dynamic_cast(phosmain->FindObject(fHeaderFile)) ; + if (phos) { + ar = dynamic_cast(phos->FindObject(fSDigitsTitle)) ; + if (ar) { + phos->Remove(ar) ; + ar->Delete() ; + delete ar ; + } + } + phosmain->Remove(phos) ; + } + + + // Digits + phos = dynamic_cast(fDigitsFolder->FindObject("PHOS")) ; + if (phos){ + ar = dynamic_cast(phos->FindObject(fDigitsTitle)) ; + if (ar) { + phos->Remove(ar) ; + ar->Delete() ; + delete ar ; + } + } + + + // EMCARecPoints + phos = dynamic_cast(fRecoFolder->FindObject("PHOS/EMCARecPoints")) ; + if (phos){ + ar = dynamic_cast(phos->FindObject(fRecPointsTitle)) ; + if (ar) { + phos->Remove(ar) ; + ar->Delete() ; + delete ar ; + } + } + + + // CPVRecPoints + phos = dynamic_cast(fRecoFolder->FindObject("PHOS/CPVRecPoints")) ; + if (phos){ + ar = dynamic_cast(phos->FindObject(fRecPointsTitle)) ; + if (ar) { + phos->Remove(ar) ; + ar->Delete() ; + delete ar ; + } + } + + + // TrackSegments + phos = dynamic_cast(fRecoFolder->FindObject("PHOS/TrackSegments")) ; + if (phos) { + ar = dynamic_cast(phos->FindObject(fTrackSegmentsTitle)) ; + if (ar) { + phos->Remove(ar) ; + ar->Delete() ; + delete ar ; + } + } + + + + // RecParticles + phos = dynamic_cast(fRecoFolder->FindObject("PHOS/RecParticles")) ; + if (phos){ + ar = dynamic_cast(phos->FindObject(fRecParticlesTitle)) ; + if (ar) { + phos->Remove(ar) ; + ar->Delete() ; + delete ar ; + } + } + + + //---- Now Tasks ----------- + + TObject * obj ; + TString sdname(fSDigitsTitle); + + // Digitizer + task = dynamic_cast(fTasksFolder->FindObject("Digitizer")) ; + if (task){ + phost = dynamic_cast(task->GetListOfTasks()->FindObject("PHOS")) ; + if (phost){ + lofTasks = phost->GetListOfTasks() ; + if (lofTasks){ + obj = lofTasks->FindObject(sdname.Data()) ; + if (obj) + lofTasks->Remove(obj) ; + } + } + } + + + sdname.Append(":") ; + // Clusterizer, TrackSegmentMaker, PID + task = dynamic_cast(fTasksFolder->FindObject("Reconstructioner")) ; + if (task){ + phost = dynamic_cast(task->GetListOfTasks()->FindObject("PHOS")) ; + if (phost){ + lofTasks = phost->GetListOfTasks() ; + TIter next(lofTasks); + while((obj=next())){ + TString oname(obj->GetName()) ; + if (oname.BeginsWith(sdname)){ + lofTasks->Remove(obj) ; + } + } + } + } + + + // SDigitizer + sdname.Append(fHeaderFile) ; + task = dynamic_cast(fTasksFolder->FindObject("SDigitizer")) ; + if (task) { + phost = dynamic_cast(task->GetListOfTasks()->FindObject("PHOS")) ; + if (phost){ + lofTasks = phost->GetListOfTasks() ; + if (lofTasks){ + obj = lofTasks->FindObject(sdname.Data()) ; + if (obj) + lofTasks->Remove(obj) ; + } + } + } + +} +//____________________________________________________________________________ +void AliPHOSGetter::SetTitle(const char * branchTitle ) +{ + fBranchTitle = branchTitle ; + fSDigitsTitle = branchTitle ; + fDigitsTitle = branchTitle ; + fRecPointsTitle = branchTitle ; + fRecParticlesTitle = branchTitle ; + fTrackSegmentsTitle = branchTitle ; + if(fToSplit){ + //First - extract full path if necessary + TString sFileName(fHeaderFile) ; + Ssiz_t islash = sFileName.Last('/') ; + if(islash(fHitsFolder->FindObject("PHOS")) ; + if ( !phosFolder ) + phosFolder = fHitsFolder->AddFolder("PHOS", "Hits from PHOS") ; + + phosFolder = dynamic_cast(fSDigitsFolder->FindObject("PHOS")) ; + if ( !phosFolder ) + phosFolder = fSDigitsFolder->AddFolder("PHOS", "SDigits from PHOS") ; + + //Make folder for SDigits + TString subdir(fHeaderFile) ; + subdir.ReplaceAll("/","_") ; + phosFolder->AddFolder(subdir, fSDigitsFileName.Data()); + + + phosFolder = dynamic_cast(fDigitsFolder->FindObject("PHOS")) ; + if ( !phosFolder ) + phosFolder = fDigitsFolder->AddFolder("PHOS", "Digits from PHOS") ; + + phosFolder = dynamic_cast(fRecoFolder->FindObject("PHOS")) ; + if ( !phosFolder ) + phosFolder = fRecoFolder->AddFolder("PHOS", "Reconstructed data from PHOS") ; + + +} +//____________________________________________________________________________ +void AliPHOSGetter::CloseSplitFiles(void){ + TFile * file ; + file = static_cast(gROOT->GetFile(fSDigitsFileName.Data() ) ) ; + if(file) + file->Close() ; + file = static_cast(gROOT->GetFile(fDigitsFileName.Data() ) ) ; + if(file) + file->Close() ; + file = static_cast(gROOT->GetFile(fRecPointsFileName.Data() ) ) ; + if(file) + file->Close() ; + file = static_cast(gROOT->GetFile(fTrackSegmentsFileName.Data() ) ) ; + if(file) + file->Close() ; + file = static_cast(gROOT->GetFile(fRecParticlesFileName.Data() ) ) ; + if(file) + file->Close() ; + +}