//_________________________________________________________________________
//*--
//*-- Author: Gines Martinez & Yves Schutz (SUBATECH)
-//*-- Compleetely redisigned by Dmitri Peressounko (SUBATECH & RRC KI) March 2001
+//*-- Compleetely redesigned by Dmitri Peressounko (SUBATECH & RRC KI) March 2001
/////////////////////////////////////////////////////////////////////////////////////
// Wrapping class for reconstruction. Allows to produce reconstruction from
// different steps: from previously produced hits,sdigits, etc. Each new reconstruction
// Thanks to getters, one can set
// parameters to reconstruction briks. The full set of parameters is saved in the
// corresponding branch: e.g. parameters of clusterizer are stored in branch
-// TreeR::AliPHOSClusterizer with the same title as the branch containing the RecPoints.
-// TTree does not support overwriting, therefore one can not produce several
+// TreeR::AliPHOSClusterizer with the same title as the branch containing the RecPoints. // TTree does not support overwriting, therefore one can not produce several
// branches with the same names and titles - use different titles.
//
// Use case:
#include "TFile.h"
// --- Standard library ---
-#include <iostream.h>
// --- AliRoot header files ---
#include "AliRun.h"
+#include "AliRunLoader.h"
#include "AliPHOSReconstructioner.h"
#include "AliPHOSClusterizerv1.h"
#include "AliPHOSDigitizer.h"
#include "AliPHOSPIDv1.h"
#include "AliPHOSFastRecParticle.h"
#include "AliPHOSCpvRecPoint.h"
+#include "AliPHOSLoader.h"
ClassImp(AliPHOSReconstructioner)
fTSMaker = 0 ;
fPID = 0 ;
fSDigitizer = 0 ;
- fHeaderFileName = "galice.root" ;
fIsInitialized = kFALSE ;
}
//____________________________________________________________________________
-AliPHOSReconstructioner::AliPHOSReconstructioner(const char* headerFile):TTask("AliPHOSReconstructioner","")
+AliPHOSReconstructioner::AliPHOSReconstructioner(const char* evFoldName,const char * branchName):
+TTask("AliPHOSReconstructioner",evFoldName)
{
// ctor
-
- fHeaderFileName = headerFile ;
-
- fSDigitsBranch="" ;
- fSDigitizer = new AliPHOSSDigitizer(fHeaderFileName.Data(),fSDigitsBranch.Data()) ;
- Add(fSDigitizer) ;
-
- fDigitsBranch="" ;
- fDigitizer = new AliPHOSDigitizer(fHeaderFileName.Data(),fDigitsBranch.Data()) ;
+ AliRunLoader* rl = AliRunLoader::GetRunLoader(evFoldName);
+ if (rl == 0x0)
+ {
+ Fatal("AliPHOSReconstructioner","Can not get Run Loader from folder %s.",evFoldName);
+ }
+ if (rl->GetAliRun() == 0x0)
+ {
+ delete gAlice;
+ gAlice = 0x0;
+ rl->LoadgAlice();
+ gAlice = rl->GetAliRun();
+ }
+
+ AliPHOSLoader* gime = dynamic_cast<AliPHOSLoader*>(rl->GetLoader("PHOSLoader"));
+ if (gime == 0x0)
+ {
+ Error("AliPHOSReconstructioner","Can not get PHOS Loader");
+ return;
+ }
+
+ TString galicefn = rl->GetFileName();
+ TString method("AliPHOSReconstructioner::AliPHOSReconstructioner(");
+ method = (((method + evFoldName)+",")+branchName)+"): ";
+
+ fSDigitsBranch= branchName;
+
+ //P.Skowronski remark
+ // Tasks has default fixed names
+ // other tasks can be added, even runtime
+ // with arbitrary name. See AliDataLoader::
+ cout<<"\n\n\n";
+ cout<<method<<"\n\nCreating SDigitizer\n";
+ fSDigitizer = new AliPHOSSDigitizer(galicefn,GetTitle());
+ Add(fSDigitizer);
+ gime->PostSDigitizer(fSDigitizer);
+
+ fDigitsBranch=branchName ;
+ cout<<"\n\n\n";
+ cout<<method<<"\n\nCreating Digitizer\n";
+ fDigitizer = new AliPHOSDigitizer(galicefn,GetTitle()) ;
Add(fDigitizer) ;
+ gime->PostDigitizer(fDigitizer);
-
- fRecPointBranch="" ;
- fClusterizer = new AliPHOSClusterizerv1(fHeaderFileName.Data(),fRecPointBranch.Data()) ;
- Add(fClusterizer) ;
+ fRecPointBranch=branchName ;
+ cout<<"\n\n\n";
+ cout<<method<<"Creating Clusterizer\n";
+ fClusterizer = new AliPHOSClusterizerv1(galicefn,GetTitle());
+ Add(fClusterizer);
+ gime->PostReconstructioner(fClusterizer);
-
- fTSBranch="" ;
- fTSMaker = new AliPHOSTrackSegmentMakerv1(fHeaderFileName.Data(),fTSBranch.Data()) ;
+ fTSBranch=branchName ;
+ fTSMaker = new AliPHOSTrackSegmentMakerv1(galicefn,GetTitle());
Add(fTSMaker) ;
+ gime->PostTracker(fTSMaker);
+
-
- fRecPartBranch="" ;
- fPID = new AliPHOSPIDv1(fHeaderFileName.Data(),fRecPartBranch.Data()) ;
- Add(fPID) ;
+ fRecPartBranch=branchName ;
+ cout<<"\n\n\n";
+ cout<<method<<"Creating PID\n";
+ fPID = new AliPHOSPIDv1(galicefn,GetTitle());
+ Add(fPID);
+ cout<<"\nFINISHED \n\n"<<method;
fIsInitialized = kTRUE ;
-
}
//____________________________________________________________________________
void AliPHOSReconstructioner::Exec(Option_t *option)
//existing
if(!fIsInitialized)
Init() ;
-
- gAlice->GetEvent(0) ;
-
- if(fSDigitizer->IsActive()&& gAlice->TreeS()){ //Will produce SDigits
-
- TBranch * sdigitsBranch = 0;
- TBranch * sdigitizerBranch = 0;
-
- TObjArray * branches = gAlice->TreeS()->GetListOfBranches() ;
- Int_t ibranch;
- Bool_t phosNotFound = kTRUE ;
- Bool_t sdigitizerNotFound = kTRUE ;
-
- for(ibranch = 0;ibranch <branches->GetEntries();ibranch++){
- if(phosNotFound){
- sdigitsBranch=(TBranch *) branches->At(ibranch) ;
- if(( strcmp("PHOS",sdigitsBranch->GetName())==0 ) &&
- (fSDigitsBranch.CompareTo(sdigitsBranch->GetTitle())== 0 ))
- phosNotFound = kFALSE ;
- }
- if(sdigitizerNotFound){
- sdigitizerBranch = (TBranch *) branches->At(ibranch) ;
- if(( strcmp(sdigitizerBranch->GetName(),"AliPHOSSDigitizer") == 0) &&
- (fSDigitsBranch.CompareTo(sdigitizerBranch->GetTitle())== 0 ) )
- sdigitizerNotFound = kFALSE ;
- }
- }
-
- if(!(sdigitizerNotFound && phosNotFound)){
- cout << "AliPHOSReconstructioner error: "<< endl ;
- cout << " Branches ''PHOS'' or ''AliPHOSSDigitizer'' with title ``" << fSDigitsBranch.Data() << "''" << endl ;
- cout << " already exist in TreeS. ROOT does not allow updating/overwriting." << endl ;
- cout << " Specify another title for branches or use ''StartFrom()'' method" << endl ;
-
- //mark all tasks as inactive
- TIter next(fTasks);
- TTask *task;
- while((task=(TTask*)next()))
- task->SetActive(kFALSE) ;
-
- return ;
- }
- }
-
- if(fDigitizer->IsActive() && gAlice->TreeD()){ //Will produce Digits
- TBranch * digitsBranch = 0;
- TBranch * digitizerBranch = 0;
-
- TObjArray * branches = gAlice->TreeD()->GetListOfBranches() ;
- Int_t ibranch;
- Bool_t phosNotFound = kTRUE ;
- Bool_t digitizerNotFound = kTRUE ;
-
- for(ibranch = 0;ibranch <branches->GetEntries();ibranch++){
- if(phosNotFound){
- digitsBranch=(TBranch *) branches->At(ibranch) ;
- if(( strcmp("PHOS",digitsBranch->GetName())==0 ) &&
- (fDigitsBranch.CompareTo(digitsBranch->GetTitle())== 0 ))
- phosNotFound = kFALSE ;
- }
- if(digitizerNotFound){
- digitizerBranch = (TBranch *) branches->At(ibranch) ;
- if(( strcmp(digitizerBranch->GetName(),"AliPHOSDigitizer") == 0) &&
- (fDigitsBranch.CompareTo(digitizerBranch->GetTitle())== 0 ) )
- digitizerNotFound = kFALSE ;
- }
- }
-
- if(!(digitizerNotFound && phosNotFound)){
- cout << "AliPHOSReconstructioner error: "<< endl ;
- cout << " Branches ''PHOS'' or ''AliPHOSDigitizer'' with title ``" << fDigitsBranch.Data() << "''" << endl ;
- cout << " already exist in TreeD. ROOT does not allow updating/overwriting." << endl ;
- cout << " Specify another title for branches or use ''StartFrom()'' method" << endl ;
-
- //mark all tasks as inactive
- TIter next(fTasks);
- TTask *task;
- while((task=(TTask*)next()))
- task->SetActive(kFALSE) ;
-
- return ;
- }
- }
-
- if(fClusterizer->IsActive() && gAlice->TreeR()){ //Will produce RecPoints
- TBranch * emcBranch = 0;
- TBranch * cpvBranch = 0;
- TBranch * clusterizerBranch = 0;
-
- TObjArray * branches = gAlice->TreeR()->GetListOfBranches() ;
- Int_t ibranch;
- Bool_t emcNotFound = kTRUE ;
- Bool_t cpvNotFound = kTRUE ;
- Bool_t clusterizerNotFound = kTRUE ;
-
- for(ibranch = 0;ibranch <branches->GetEntries();ibranch++){
-
- if(emcNotFound){
- emcBranch=(TBranch *) branches->At(ibranch) ;
- if(fRecPointBranch.CompareTo(emcBranch->GetTitle())==0 )
- if( strcmp(emcBranch->GetName(),"PHOSEmcRP") == 0)
- emcNotFound = kFALSE ;
- }
- if(cpvNotFound){
- cpvBranch=(TBranch *) branches->At(ibranch) ;
- if(fRecPointBranch.CompareTo(cpvBranch->GetTitle())==0 )
- if( strcmp(cpvBranch->GetName(),"PHOSCpvRP") == 0)
- cpvNotFound = kFALSE ;
- }
- if(clusterizerNotFound){
- clusterizerBranch = (TBranch *) branches->At(ibranch) ;
- if( fRecPointBranch.CompareTo(clusterizerBranch->GetTitle()) == 0)
- if( strcmp(clusterizerBranch->GetName(),"AliPHOSClusterizer") == 0)
- clusterizerNotFound = kFALSE ;
- }
- }
-
- if(!(clusterizerNotFound && emcNotFound && cpvNotFound)){
- cout << "AliPHOSReconstructioner error: "<< endl ;
- cout << " Branches ''PHOSEmcRP'', ''PHOSCpvRP'' or ''AliPHOSClusterizer'' with title ``"
- << fRecPointBranch.Data() << "''" << endl ;
- cout << " already exist in TreeR. ROOT does not allow updating/overwriting." << endl ;
- cout << " Specify another title for branches or use ''StartFrom()'' method" << endl ;
-
- //mark all tasks as inactive
- TIter next(fTasks);
- TTask *task;
- while((task=(TTask*)next()))
- task->SetActive(kFALSE) ;
- return ;
- }
- }
-
- if(fTSMaker->IsActive() && gAlice->TreeR()){ //Produce TrackSegments
-
- TBranch * tsMakerBranch = 0;
- TBranch * tsBranch = 0;
-
- TObjArray * branches = gAlice->TreeR()->GetListOfBranches() ;
- Int_t ibranch;
- Bool_t tsMakerNotFound = kTRUE ;
- Bool_t tsNotFound = kTRUE ;
-
- for(ibranch = 0;(ibranch <branches->GetEntries())&&(tsMakerNotFound||tsNotFound);ibranch++){
- if(tsMakerNotFound){
- tsMakerBranch=(TBranch *) branches->At(ibranch) ;
- if( fTSBranch.CompareTo(tsMakerBranch->GetTitle())==0 )
- if( strcmp(tsMakerBranch->GetName(),"AliPHOSTrackSegmentMaker") == 0)
- tsMakerNotFound = kFALSE ;
- }
- if(tsNotFound){
- tsBranch=(TBranch *) branches->At(ibranch) ;
- if( fTSBranch.CompareTo(tsBranch->GetTitle())==0 )
- if( strcmp(tsBranch->GetName(),"PHOSTS") == 0)
- tsNotFound = kFALSE ;
- }
- }
-
- if(!(tsMakerNotFound &&tsNotFound) ){
- cout << "AliPHOSReconstructioner error: "<< endl ;
- cout << " Branches ''PHOSTS'' or ''AliPHOSTrackSegmentMaker'' with title ``"
- << fTSBranch.Data() << "''" << endl ;
- cout << " already exist in TreeR. ROOT does not allow updating/overwriting." << endl ;
- cout << " Specify another title for branches or use ''StartFrom()'' method" << endl ;
-
- //mark all tasks as inactive
- TIter next(fTasks);
- TTask *task;
- while((task=(TTask*)next()))
- task->SetActive(kFALSE) ;
- return ;
-
- }
-
- }
-
- if(fPID->IsActive() && gAlice->TreeR()){ //Produce RecParticles
- TBranch * pidBranch = 0;
- TBranch * rpBranch = 0;
-
- TObjArray * branches = gAlice->TreeR()->GetListOfBranches() ;
- Int_t ibranch;
- Bool_t pidNotFound = kTRUE ;
- Bool_t rpNotFound = kTRUE ;
-
- for(ibranch = 0;(ibranch <branches->GetEntries()) && pidNotFound && rpNotFound ;ibranch++){
- if(pidNotFound){
- pidBranch=(TBranch *) branches->At(ibranch) ;
- if( (strcmp(fRecPartBranch,pidBranch->GetTitle())==0 ) &&
- (strcmp(pidBranch->GetName(),"AliPHOSPID") == 0) )
- pidNotFound = kFALSE ;
- }
- if(rpNotFound){
- rpBranch=(TBranch *) branches->At(ibranch) ;
- if( (strcmp(fRecPartBranch,rpBranch->GetTitle())==0 ) &&
- (strcmp(rpBranch->GetName(),"PHOSRP") == 0) )
- rpNotFound = kFALSE ;
- }
- }
-
- if(!pidNotFound || !rpNotFound ){
- cout << "AliPHOSReconstructioner error: "<< endl ;
- cout << " Branches ''PHOSRP'' or ''AliPHOSPID'' with title ``"
- << fRecPartBranch.Data() << "''" << endl ;
- cout << " already exist in TreeR. ROOT does not allow updating/overwriting." << endl ;
- cout << " Specify another title for branches." << endl ;
-
- //mark all tasks as inactive
- TIter next(fTasks);
- TTask *task;
- while((task=(TTask*)next()))
- task->SetActive(kFALSE) ;
- return ;
- }
-
- }
}
//____________________________________________________________________________
void AliPHOSReconstructioner::Init()
if(!fIsInitialized){
// Initialisation
- fSDigitsBranch="" ;
- fSDigitizer = new AliPHOSSDigitizer(fHeaderFileName.Data(),fSDigitsBranch.Data()) ;
+ fSDigitsBranch="Default" ;
+ fSDigitizer = new AliPHOSSDigitizer(GetTitle(),fSDigitsBranch.Data()) ;
Add(fSDigitizer) ;
- fDigitsBranch="" ;
- fDigitizer = new AliPHOSDigitizer(fHeaderFileName.Data(),fDigitsBranch.Data()) ;
+ fDigitsBranch="Default" ;
+ fDigitizer = new AliPHOSDigitizer(GetTitle(),fDigitsBranch.Data());
Add(fDigitizer) ;
- fRecPointBranch="" ;
- fClusterizer = new AliPHOSClusterizerv1(fHeaderFileName.Data(),fRecPointBranch.Data()) ;
+ fRecPointBranch="Default" ;
+ fClusterizer = new AliPHOSClusterizerv1(GetTitle(),fRecPointBranch.Data());
Add(fClusterizer) ;
- fTSBranch="" ;
- fTSMaker = new AliPHOSTrackSegmentMakerv1(fHeaderFileName.Data(),fTSBranch.Data()) ;
+ fTSBranch="Default" ;
+ fTSMaker = new AliPHOSTrackSegmentMakerv1(GetTitle(),fTSBranch.Data());
Add(fTSMaker) ;
- fRecPartBranch="" ;
- fPID = new AliPHOSPIDv1(fHeaderFileName.Data(),fRecPartBranch.Data()) ;
+ fRecPartBranch="Default";
+ fPID = new AliPHOSPIDv1(GetTitle(),fRecPartBranch.Data()) ;
Add(fPID) ;
fIsInitialized = kTRUE ;
+
}
}
//____________________________________________________________________________
AliPHOSReconstructioner::~AliPHOSReconstructioner()
{
// Delete data members if any
-
- if(fSDigitizer)
- delete fSDigitizer ;
-
- if(fDigitizer)
- delete fDigitizer ;
-
- if(fClusterizer)
- delete fClusterizer ;
-
- if(fTSMaker)
- delete fTSMaker ;
-
- if(fPID)
- delete fPID ;
-
-// TFile * file = (TFile*) gROOT->GetFile(fHeaderFileName.Data()) ;
-
-// if(file != 0) {
-// file->Close();
-// delete file;
-// printf("File %s is closed\n",fHeaderFileName.Data());
-// }
-
}
-//____________________________________________________________________________
-void AliPHOSReconstructioner::SetBranchTitle(const char* branch, const char * title)
-{
- //Diverge correcpoinding branch to the file "title"
-
- if(strcmp(branch,"SDigits") == 0){
- fSDigitizer->SetSDigitsBranch(title) ;
- fDigitizer->SetSDigitsBranch(title) ;
- fSDigitsBranch = title ;
- return ;
- }
-
- if(strcmp(branch,"Digits") == 0){
- fDigitizer->SetDigitsBranch(title) ;
- fClusterizer->SetDigitsBranch(title) ;
- fDigitsBranch = title ;
- return ;
- }
-
- if(strcmp(branch,"RecPoints") == 0){
- fClusterizer->SetRecPointsBranch(title) ;
- fTSMaker->SetRecPointsBranch(title) ;
- fRecPointBranch = title ;
- return ;
- }
-
- if(strcmp(branch,"TrackSegments") == 0){
- fTSMaker->SetTrackSegmentsBranch(title) ;
- fPID->SetTrackSegmentsBranch(title) ;
- fTSBranch = title ;
- return ;
- }
-
- if(strcmp(branch,"RecParticles") == 0){
- fPID->SetRecParticlesBranch(title) ;
- fRecPartBranch = title ;
- return ;
- }
-
- cout << "There is no branch " << branch << "!"<< endl ;
- cout << "Available branches `SDigits', `Digits', `RecPoints', `TrackSegments' and `RecParticles' " << endl ;
-
-}
-//____________________________________________________________________________
-void AliPHOSReconstructioner::StartFrom(char * module,char* title)
-{
- // in the next pass of reconstruction (call ExecuteTask()) reconstruction will
- // start from the module "module", and in the case of non zero title all
- // pruduced branches will have title "title". The following "modules" are recognized
- // "SD" - AliPHOSSDigitizer,
- // "D" - AliPHOSDigitizer
- // "C" - AliPHOSClusterizer
- // "TS" - AliPHOSTrackSegmentMaker
- // "RP" - AliPHOSPID
-
- if(!fIsInitialized)
- Init() ;
-
- char * moduleName = new char[30];
- if(strstr(module,"SD"))
- sprintf(moduleName,"AliPHOSSDigitizer") ;
- else
- if(strstr(module,"D") )
- sprintf(moduleName,"AliPHOSDigitizer") ;
- else
- if(strstr(module,"C") || strstr(module,"RecPoint") )
- sprintf(moduleName,"AliPHOSClusterizer") ;
- else
- if(strstr(module,"TS") || strstr(module,"Track") )
- sprintf(moduleName,"AliPHOSTrackSegmentMaker") ;
- else
- if(strstr(module,"PID") || strstr(module,"Particle") || strstr(module,"RP") )
- sprintf(moduleName,"AliPHOSPID") ;
- else{
- cout << "Do not know such a module / Rec Object " << endl;
- return ;
- }
-
- TIter next(fTasks);
- TTask *task;
- Bool_t active = kFALSE ;
- while((task=(TTask*)next())){
- if (strcmp(moduleName,task->GetName())==0)
- active = kTRUE;
- task->SetActive(active) ;
- if(active && title){ // set title to branches
- switch(strlen(task->GetName()) ) {
- case 17: // "AliPHOSSDigitizer"
- fSDigitizer->SetSDigitsBranch(title) ;
- fDigitizer->SetSDigitsBranch(title) ;
- fSDigitsBranch = title ;
- break ;
- case 16: //"AliPHOSDigitizer"
- fDigitizer->SetDigitsBranch(title) ;
- fClusterizer->SetDigitsBranch(title) ;
- fDigitsBranch = title ;
- break ;
- case 18: //"AliPHOSClusterizer"
- fClusterizer->SetRecPointsBranch(title) ;
- fTSMaker->SetRecPointsBranch(title) ;
- fRecPointBranch = title ;
- break ;
- case 24: //"AliPHOSTrackSegmentMaker"
- fTSMaker->SetTrackSegmentsBranch(title) ;
- fPID->SetTrackSegmentsBranch(title) ;
- fTSBranch = title ;
- break ;
- case 10: // "AliPHOSPID"
- fPID->SetRecParticlesBranch(title) ;
- fRecPartBranch = title ;
- break ;
- }
-
- }
- }
-
- delete moduleName;
-}
-//____________________________________________________________________________
void AliPHOSReconstructioner::Print(Option_t * option)const {
// Print reconstructioner data
- cout << "-----------------AliPHOSReconstructioner---------------" << endl ;
- cout << " Reconstruction of the header file " <<fHeaderFileName.Data() << endl ;
- cout << " with the following modules: " << endl ;
+ TString message ;
+ message = "-----------------AliPHOSReconstructioner---------------\n" ;
+ message += " Reconstruction of the header file %s\n" ;
+ message += " with the following modules:\n" ;
if(fSDigitizer->IsActive()){
- cout << " (+) " << fSDigitizer->GetName() << " to branch : " << fSDigitsBranch.Data() << endl ;
- cout << endl ;
+ message += " (+) %s to branch %s\n" ;
}
if(fDigitizer->IsActive()){
- cout << " (+) " << fDigitizer->GetName() << " to branch : " << fDigitsBranch.Data() << endl ;
- cout << endl ;
+ message += " (+) %s to branch %s\n" ;
}
if(fClusterizer->IsActive()){
- cout << " (+) " <<fClusterizer->GetName() << " to branch : " <<fRecPointBranch.Data() << endl ;
- cout << endl ;
+ message += " (+) %s to branch %s\n" ;
}
if(fTSMaker->IsActive()){
- cout << " (+) " << fTSMaker->GetName() << " to branch : " << fTSBranch.Data() << endl ;
- cout << endl ;
+ message += " (+) %s to branch %s\n" ;
}
-
if(fPID->IsActive()){
- cout << " (+) " << fPID->GetName() << " to branch : " <<fRecPartBranch.Data() << endl ;
- cout << endl ;
- }
-
-
+ message += " (+) %s to branch %s\n" ;
+ }
+ Info("Print", message.Data(),
+ GetTitle(),
+ fSDigitizer->GetName(), fSDigitsBranch.Data(),
+ fDigitizer->GetName(), fDigitsBranch.Data() ,
+ fClusterizer->GetName(), fRecPointBranch.Data(),
+ fTSMaker->GetName(), fTSBranch.Data() ,
+ fPID->GetName(), fRecPartBranch.Data() ) ;
}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff