/************************************************************************** * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * * * * Author: The ALICE Off-line Project. * * Contributors are mentioned in the code where appropriate. * * * * Permission to use, copy, modify and distribute this software and its * * documentation strictly for non-commercial purposes is hereby granted * * without fee, provided that the above copyright notice appears in all * * copies and that both the copyright notice and this permission notice * * appear in the supporting documentation. The authors make no claims * * about the suitability of this software for any purpose. It is * * provided "as is" without express or implied warranty. * **************************************************************************/ /* $Id$ */ // Selector base class for analysis based on ESD // Please derive your selector-based analysis from this class, if you just want to use // information from the ESD. // // The ESD is available as member fESD // // The following methods can be overrriden. Please do not forgot to call the base class function. // // Begin(): called everytime a loop on the tree starts, // a convenient place to create your histograms. // SlaveBegin(): called after Begin(), when on PROOF called only on the // slave servers. // Init(): called for each new tree. Enable/Disable branches here. // Process(): called for each event, in this function you decide what // to read and fill your histograms. // SlaveTerminate: called at the end of the loop on the tree, when on PROOF // called only on the slave servers. // Terminate(): called at the end of the loop on the tree, // a convenient place to draw/fit your histograms. // // Author: Jan.Fiete.Grosse-Oetringhaus@cern.ch #include "AliSelector.h" #include #include #include #include #include #include #include #include #include #include #include "AliLog.h" #include "AliESD.h" ClassImp(AliSelector) AliSelector::AliSelector() : TSelector(), fTree(0), fESD(0), fCountFiles(0) { // // Constructor. Initialization of pointers // } AliSelector::~AliSelector() { // // Destructor // if (fTree) fTree->ResetBranchAddresses(); if (fESD) { delete fESD; fESD = 0; } } void AliSelector::CheckOptions() { // checks the option string for the debug flag AliLog::SetClassDebugLevel(ClassName(), AliLog::kInfo); TString option = GetOption(); if (option.Contains("moredebug")) { printf("Enabling verbose debug mode for %s\n", ClassName()); AliLog::SetClassDebugLevel(ClassName(), AliLog::kDebug+1); AliInfo(Form("Called with option %s.", option.Data())); } else if (option.Contains("debug")) { printf("Enabling debug mode for %s\n", ClassName()); AliLog::SetClassDebugLevel(ClassName(), AliLog::kDebug); AliInfo(Form("Called with option %s.", option.Data())); } } void AliSelector::Begin(TTree*) { // The Begin() function is called at the start of the query. // When running with PROOF Begin() is only called on the client. // The tree argument is deprecated (on PROOF 0 is passed). CheckOptions(); AliDebug(AliLog::kDebug, "============BEGIN==========="); } void AliSelector::SlaveBegin(TTree* tree) { // The SlaveBegin() function is called after the Begin() function. // When running with PROOF SlaveBegin() is called on each slave server. // The tree argument is deprecated (on PROOF 0 is passed). CheckOptions(); AliDebug(AliLog::kDebug, "=======SLAVEBEGIN========"); AliDebug(AliLog::kDebug, Form("Hostname: %s", gSystem->HostName())); AliDebug(AliLog::kDebug, Form("Time: %s", gSystem->Now().AsString())); if (tree != 0) Init(tree); } void AliSelector::Init(TTree *tree) { // The Init() function is called when the selector needs to initialize // a new tree or chain. Typically here the branch addresses of the tree // will be set. It is normaly not necessary to make changes to the // generated code, but the routine can be extended by the user if needed. // Init() will be called many times when running with PROOF. AliDebug(AliLog::kDebug, "=========Init=========="); fTree = tree; if (fTree == 0) { AliDebug(AliLog::kError, "ERROR: tree argument is 0."); return; } // Set branch address fTree->SetBranchAddress("ESD", &fESD); if (fESD != 0) AliDebug(AliLog::kInfo, "INFO: Found ESD branch in chain."); } Bool_t AliSelector::Notify() { // The Notify() function is called when a new file is opened. This // can be either for a new TTree in a TChain or when when a new TTree // is started when using PROOF. Typically here the branch pointers // will be retrieved. It is normaly not necessary to make changes // to the generated code, but the routine can be extended by the // user if needed. AliDebug(AliLog::kDebug, "=========NOTIFY=========="); AliDebug(AliLog::kDebug, Form("Hostname: %s", gSystem->HostName())); AliDebug(AliLog::kDebug, Form("Time: %s", TTimeStamp(time(0)).AsString())); ++fCountFiles; if (fTree) { TFile *f = fTree->GetCurrentFile(); if (f) { AliDebug(AliLog::kInfo, Form("Processing %d. file %s", fCountFiles, f->GetName())); } else AliDebug(AliLog::kError, "fTree->GetCurrentFile() is 0"); } else { AliDebug(AliLog::kError, "fTree not available"); } return kTRUE; } Bool_t AliSelector::Process(Long64_t entry) { // The Process() function is called for each entry in the tree (or possibly // keyed object in the case of PROOF) to be processed. The entry argument // specifies which entry in the currently loaded tree is to be processed. // It can be passed to either TTree::GetEntry() or TBranch::GetEntry() // to read either all or the required parts of the data. When processing // keyed objects with PROOF, the object is already loaded and is available // via the fObject pointer. // // This function should contain the "body" of the analysis. It can contain // simple or elaborate selection criteria, run algorithms on the data // of the event and typically fill histograms. // WARNING when a selector is used with a TChain, you must use // the pointer to the current TTree to call GetEntry(entry). // The entry is always the local entry number in the current tree. // Assuming that fTree is the pointer to the TChain being processed, // use fTree->GetTree()->GetEntry(entry). AliDebug(AliLog::kDebug, Form("=========PROCESS========== Entry %lld", entry)); if (!fTree) { AliDebug(AliLog::kError, "ERROR: fTree is 0."); return kFALSE; } fTree->GetTree()->GetEntry(entry); if (fESD) AliDebug(AliLog::kDebug, Form("ESD: We have %d tracks.", fESD->GetNumberOfTracks())); return kTRUE; } void AliSelector::SlaveTerminate() { // The SlaveTerminate() function is called after all entries or objects // have been processed. When running with PROOF SlaveTerminate() is called // on each slave server. AliDebug(AliLog::kDebug, "=======SLAVETERMINATE======="); } void AliSelector::Terminate() { // The Terminate() function is the last function to be called during // a query. It always runs on the client, it can be used to present // the results graphically or save the results to file. AliDebug(AliLog::kDebug, "=========TERMINATE=========="); }