[u/mrichter/AliRoot.git] / ANALYSIS / AliAnalysisManager.h

#ifndef ALIANALYSISMANAGER_H
#define ALIANALYSISMANAGER_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

/* $Id$ */
// Author: Andrei Gheata, 31/05/2006

//==============================================================================
//   AliAnalysysManager - Manager analysis class. Allows creation of several
// analysis tasks and data containers storing their input/output. Allows 
// connecting/chaining tasks via shared data containers. Serializes the current
// event for all tasks depending only on initial input data.
//==============================================================================

#ifndef ROOT_TNamed
#include <TNamed.h>
#endif
#ifndef ROOT_THashTable
#include <THashTable.h>
#endif

class TClass;
class TTree;
class TFile;
class TStopwatch;
class AliAnalysisSelector;
class AliAnalysisDataContainer;
class AliAnalysisTask;
class AliVEventHandler;
class AliVEventPool;
class AliAnalysisGrid;


class AliAnalysisManager : public TNamed {

public:

enum EAliAnalysisContType {
   kExchangeContainer  = 0,   // use to exchange data between tasks
   kInputContainer     = 1,   // use for the task private input data
   kOutputContainer    = 2,   // use for the task private output data
   kParamContainer     = 3    // use for storing task private parameters/cuts per analysis session
};   

enum EAliAnalysisExecMode {
   kLocalAnalysis    = 0,
   kProofAnalysis    = 1,
   kGridAnalysis     = 2,
   kMixingAnalysis   = 3
};

enum EAliAnalysisFlags {
   kEventLoop        = BIT(14),
   kDisableBranches  = BIT(15),
   kUseDataSet       = BIT(16),
   kSaveCanvases     = BIT(17),
   kExternalLoop     = BIT(18),
   kSkipTerminate    = BIT(19),
   kUseProgressBar   = BIT(20)
};   

   AliAnalysisManager(const char *name = "mgr", const char *title="");
   virtual            ~AliAnalysisManager();

   AliAnalysisManager(const AliAnalysisManager& other);
   AliAnalysisManager& operator=(const AliAnalysisManager& other);
   
   // Event loop control
   virtual Int_t       GetEntry(Long64_t entry, Int_t getall = 0);
   virtual Bool_t      Init(TTree *tree);   
   virtual Bool_t      Notify();
   virtual Bool_t      ProcessCut(Long64_t entry) {return Process(entry);}
   virtual Bool_t      Process(Long64_t entry);
   TFile              *OpenProofFile(AliAnalysisDataContainer *cont, const char *option);
   static TFile       *OpenFile(AliAnalysisDataContainer *cont, const char *option, Bool_t ignoreProof=kFALSE);
   void                PackOutput(TList *target);
   void                RegisterExtraFile(const char *fname);
   Long64_t            StartAnalysis(const char *type="local", TTree * const tree=0, Long64_t nentries=1234567890, Long64_t firstentry=0);
   Long64_t            StartAnalysis(const char *type, const char *dataset, Long64_t nentries=1234567890, Long64_t firstentry=0);
   virtual void        SlaveBegin(TTree *tree);
   virtual void        Terminate();
   void                UnpackOutput(TList *source);

   // Getters/Setters
   static AliAnalysisManager *GetAnalysisManager() {return fgAnalysisManager;}
   EAliAnalysisExecMode 
                       GetAnalysisType() const    {return fMode;}
   void                GetAnalysisTypeString(TString &type) const;                    
   Bool_t              GetAutoBranchLoading() const {return fAutoBranchHandling;} 
   static const char  *GetCommonFileName()        {return fgCommonFileName.Data();}
   AliAnalysisDataContainer *
                       GetCommonInputContainer() const  {return fCommonInput;}
   AliAnalysisDataContainer *
                       GetCommonOutputContainer() const {return fCommonOutput;}
   TObjArray          *GetContainers() const      {return fContainers;}
   Long64_t            GetCurrentEntry() const    {return fCurrentEntry;}
   UInt_t              GetDebugLevel() const      {return fDebug;}
   TString             GetExtraFiles() const      {return fExtraFiles;}
   AliVEventPool*      GetEventPool()  const      {return fEventPool;}
   Bool_t              GetFileFromWrapper(const char *filename, const TList *source);
   AliAnalysisGrid*    GetGridHandler()           {return fGridHandler;}
   TObjArray          *GetInputs() const          {return fInputs;}
   AliVEventHandler*   GetInputEventHandler() const   {return fInputEventHandler;}
   AliVEventHandler*   GetMCtruthEventHandler() const {return fMCtruthEventHandler;}
   AliVEventHandler*   GetOutputEventHandler() const  {return fOutputEventHandler;}
   TObjArray          *GetOutputs() const         {return fOutputs;}
   TObjArray          *GetParamOutputs() const    {return fParamCont;}
   TObjArray          *GetTasks() const           {return fTasks;}
   TObjArray          *GetTopTasks() const        {return fTopTasks;}
   TTree              *GetTree() const            {return fTree;}
   TObjArray          *GetZombieTasks() const     {return fZombies;}
   Bool_t              IsProofMode() const        {return (fMode==kProofAnalysis)?kTRUE:kFALSE;}
   Bool_t              IsRemote() const           {return fIsRemote;}
   Bool_t              IsUsingDataSet() const     {return TObject::TestBit(kUseDataSet);}
   void                LoadBranch(const char *n)  { if(fAutoBranchHandling) return; DoLoadBranch(n); }
   void                SetAnalysisType(EAliAnalysisExecMode mode) {fMode = mode;}
   void                SetAutoBranchLoading(Bool_t b) { fAutoBranchHandling = b; }
   void                SetCurrentEntry(Long64_t entry)            {fCurrentEntry = entry;}
   void                SetCollectSysInfoEach(Int_t nevents=0)     {fNSysInfo = nevents;}
   static void         SetCommonFileName(const char *name)        {fgCommonFileName = name;}
   void                SetDebugLevel(UInt_t level)                {fDebug = level;}
   void                SetDisableBranches(Bool_t disable=kTRUE)   {TObject::SetBit(kDisableBranches,disable);}
   void                SetExternalLoop(Bool_t flag)               {TObject::SetBit(kExternalLoop,flag);}
   void                SetEventPool(AliVEventPool* const epool)   {fEventPool = epool;}
   void                SetGridHandler(AliAnalysisGrid * const handler) {fGridHandler = handler;}
   void                SetInputEventHandler(AliVEventHandler* const handler);
   void                SetMCtruthEventHandler(AliVEventHandler* const handler) {fMCtruthEventHandler = handler;}
   void                SetNSysInfo(Long64_t nevents)              {fNSysInfo = nevents;}
   void                SetOutputEventHandler(AliVEventHandler* const handler);
   void                SetSelector(AliAnalysisSelector * const sel)      {fSelector = sel;}
   void                SetSaveCanvases(Bool_t flag=kTRUE)         {TObject::SetBit(kSaveCanvases,flag);}
   void                SetSkipTerminate(Bool_t flag)              {TObject::SetBit(kSkipTerminate,flag);}
   void                SetUseProgressBar(Bool_t flag)             {TObject::SetBit(kUseProgressBar,flag);}
   void                SetSpecialOutputLocation(const char *loc)  {fSpecialOutputLocation = loc;}

   // Container handling
   AliAnalysisDataContainer *CreateContainer(const char *name, TClass *datatype, 
                       EAliAnalysisContType type     = kExchangeContainer, 
                       const char          *filename = NULL);
   
   // Including tasks and getting them
   void                 AddTask(AliAnalysisTask *task);
   AliAnalysisTask     *GetTask(const char *name) const;
   
   // Connecting data containers to task inputs/outputs
   Bool_t               ConnectInput(AliAnalysisTask *task, Int_t islot,
                                     AliAnalysisDataContainer *cont);
   Bool_t               ConnectOutput(AliAnalysisTask *task, Int_t islot,
                                     AliAnalysisDataContainer *cont);
   // Garbage collection
   void                 CleanContainers();
   
   // Analysis initialization and execution, status
   Bool_t               InitAnalysis();
   Bool_t               IsInitialized() const {return fInitOK;}
   Bool_t               IsExternalLoop() const {return TObject::TestBit(kExternalLoop);}
   Bool_t               IsEventLoop() const {return TObject::TestBit(kEventLoop);}
   Bool_t               IsSkipTerminate() const {return TObject::TestBit(kSkipTerminate);}
   void                 ResetAnalysis();
   void                 ExecAnalysis(Option_t *option="");
   void                 PrintStatus(Option_t *option="all") const;
   void                 ProfileTask(const char *name, const char *option="VM") const;
   void                 ProfileTask(Int_t itop, const char *option="VM") const;
   static void          ProgressBar(const char *opname, Long64_t current, Long64_t size, TStopwatch * const watch=0, Bool_t last=kFALSE, Bool_t refresh=kFALSE);
   Bool_t               ValidateOutputFiles() const;

protected:
   void                 ImportWrappers(TList *source);
   void                 SetEventLoop(Bool_t flag=kTRUE) {TObject::SetBit(kEventLoop,flag);}
   void                 DoLoadBranch(const char *name);

private:
   TTree                  *fTree;                //! Input tree in case of TSelector model
   AliVEventHandler       *fInputEventHandler;   //  Optional common input  event handler
   AliVEventHandler       *fOutputEventHandler;  //  Optional common output event handler
   AliVEventHandler       *fMCtruthEventHandler; //  Optional common MC Truth event handler
   AliVEventPool          *fEventPool;           //  Event pool for mixing analysis
   Long64_t                fCurrentEntry;        //! Current processed entry in the tree
   Long64_t                fNSysInfo;            // Event frequency for collecting system information
   EAliAnalysisExecMode    fMode;                // Execution mode
   Bool_t                  fInitOK;              // Initialisation done
   Bool_t                  fIsRemote;            //! Flag is set for remote analysis
   UInt_t                  fDebug;               // Debug level
   TString                 fSpecialOutputLocation; // URL/path where the special outputs will be copied
   TObjArray              *fTasks;               // List of analysis tasks
   TObjArray              *fTopTasks;            // List of top tasks
   TObjArray              *fZombies;             // List of zombie tasks
   TObjArray              *fContainers;          // List of all containers
   TObjArray              *fInputs;              // List of containers with input data
   TObjArray              *fOutputs;             // List of containers with results
   TObjArray              *fParamCont;           // List of containers with results
   AliAnalysisDataContainer *fCommonInput;       // Common input container
   AliAnalysisDataContainer *fCommonOutput;      // Common output container
   AliAnalysisSelector    *fSelector;            //! Current selector
   AliAnalysisGrid        *fGridHandler;         //! Grid handler plugin
   TString                 fExtraFiles;          // List of extra files to be merged
   Bool_t                  fAutoBranchHandling;  // def=kTRUE, turn off if you use LoadBranch
   THashTable              fTable;               // keep branch ptrs in case of manual branch loading

   static TString          fgCommonFileName;     //! Common output file name (not streamed)
   static AliAnalysisManager *fgAnalysisManager; //! static pointer to object instance
   ClassDef(AliAnalysisManager,8)  // Analysis manager class
};   
#endif
Commit	Line	Data
d3106602	1	#ifndef ALIANALYSISMANAGER_H
	2	#define ALIANALYSISMANAGER_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	/* $Id$ */
	7	// Author: Andrei Gheata, 31/05/2006
	8
	9	//==============================================================================
	10	// AliAnalysysManager - Manager analysis class. Allows creation of several
	11	// analysis tasks and data containers storing their input/output. Allows
	12	// connecting/chaining tasks via shared data containers. Serializes the current
	13	// event for all tasks depending only on initial input data.
	14	//==============================================================================
	15
c52c2132	16	#ifndef ROOT_TNamed
c52c2132	17	#include <TNamed.h>
d3106602	18	#endif
012e169c	19	#ifndef ROOT_THashTable
	20	#include <THashTable.h>
	21	#endif
d3106602	22
d3106602	23	class TClass;
c52c2132	24	class TTree;
8d7d3b59	25	class TFile;
cd11251e	26	class TStopwatch;
8d7d3b59	27	class AliAnalysisSelector;
d3106602	28	class AliAnalysisDataContainer;
d3106602	29	class AliAnalysisTask;
d2f1d9ef	30	class AliVEventHandler;
c2922515	31	class AliVEventPool;
c57f56b7	32	class AliAnalysisGrid;
8c9485b2	33
d3106602	34
c52c2132	35	class AliAnalysisManager : public TNamed {
d3106602	36
	37	public:
	38
	39	enum EAliAnalysisContType {
6cd44ee0	40	kExchangeContainer = 0, // use to exchange data between tasks
	41	kInputContainer = 1, // use for the task private input data
	42	kOutputContainer = 2, // use for the task private output data
	43	kParamContainer = 3 // use for storing task private parameters/cuts per analysis session
d3106602	44	};
c52c2132	45
	46	enum EAliAnalysisExecMode {
	47	kLocalAnalysis = 0,
	48	kProofAnalysis = 1,
aee5ee44	49	kGridAnalysis = 2,
aee5ee44	50	kMixingAnalysis = 3
c52c2132	51	};
c52c2132	52
efd53803	53	enum EAliAnalysisFlags {
b1310ef5	54	kEventLoop = BIT(14),
d86ed856	55	kDisableBranches = BIT(15),
a0e2e8b1	56	kUseDataSet = BIT(16),
06a59280	57	kSaveCanvases = BIT(17),
90a4b3ee	58	kExternalLoop = BIT(18),
ff0d33e8	59	kSkipTerminate = BIT(19),
ff0d33e8	60	kUseProgressBar = BIT(20)
efd53803	61	};
efd53803	62
6e14aaa2	63	AliAnalysisManager(const char name = "mgr", const char title="");
d3106602	64	virtual ~AliAnalysisManager();
c52c2132	65
	66	AliAnalysisManager(const AliAnalysisManager& other);
	67	AliAnalysisManager& operator=(const AliAnalysisManager& other);
d3106602	68
06a59280	69	// Event loop control
06a59280	70	virtual Int_t GetEntry(Long64_t entry, Int_t getall = 0);
2d626244	71	virtual Bool_t Init(TTree *tree);
327eaf46	72	virtual Bool_t Notify();
327eaf46	73	virtual Bool_t ProcessCut(Long64_t entry) {return Process(entry);}
d3106602	74	virtual Bool_t Process(Long64_t entry);
84fcd93f	75	TFile OpenProofFile(AliAnalysisDataContainer cont, const char *option);
84fcd93f	76	static TFile OpenFile(AliAnalysisDataContainer cont, const char *option, Bool_t ignoreProof=kFALSE);
c52c2132	77	void PackOutput(TList *target);
06a59280	78	void RegisterExtraFile(const char *fname);
61505f8b	79	Long64_t StartAnalysis(const char type="local", TTree const tree=0, Long64_t nentries=1234567890, Long64_t firstentry=0);
1f87e9fb	80	Long64_t StartAnalysis(const char type, const char dataset, Long64_t nentries=1234567890, Long64_t firstentry=0);
06a59280	81	virtual void SlaveBegin(TTree *tree);
d3106602	82	virtual void Terminate();
06a59280	83	void UnpackOutput(TList *source);
d3106602	84
c52c2132	85	// Getters/Setters
c52c2132	86	static AliAnalysisManager *GetAnalysisManager() {return fgAnalysisManager;}
c52c2132	87	EAliAnalysisExecMode
06a59280	88	GetAnalysisType() const {return fMode;}
d29168d6	89	void GetAnalysisTypeString(TString &type) const;
012e169c	90	Bool_t GetAutoBranchLoading() const {return fAutoBranchHandling;}
06a59280	91	static const char *GetCommonFileName() {return fgCommonFileName.Data();}
06a59280	92	AliAnalysisDataContainer *
61505f8b	93	GetCommonInputContainer() const {return fCommonInput;}
06a59280	94	AliAnalysisDataContainer *
61505f8b	95	GetCommonOutputContainer() const {return fCommonOutput;}
06a59280	96	TObjArray *GetContainers() const {return fContainers;}
	97	Long64_t GetCurrentEntry() const {return fCurrentEntry;}
	98	UInt_t GetDebugLevel() const {return fDebug;}
	99	TString GetExtraFiles() const {return fExtraFiles;}
61505f8b	100	AliVEventPool* GetEventPool() const {return fEventPool;}
61505f8b	101	Bool_t GetFileFromWrapper(const char filename, const TList source);
06a59280	102	AliAnalysisGrid* GetGridHandler() {return fGridHandler;}
06a59280	103	TObjArray *GetInputs() const {return fInputs;}
61505f8b	104	AliVEventHandler* GetInputEventHandler() const {return fInputEventHandler;}
	105	AliVEventHandler* GetMCtruthEventHandler() const {return fMCtruthEventHandler;}
	106	AliVEventHandler* GetOutputEventHandler() const {return fOutputEventHandler;}
06a59280	107	TObjArray *GetOutputs() const {return fOutputs;}
ff07ec61	108	TObjArray *GetParamOutputs() const {return fParamCont;}
06a59280	109	TObjArray *GetTasks() const {return fTasks;}
	110	TObjArray *GetTopTasks() const {return fTopTasks;}
	111	TTree *GetTree() const {return fTree;}
	112	TObjArray *GetZombieTasks() const {return fZombies;}
d791fad6	113	Bool_t IsProofMode() const {return (fMode==kProofAnalysis)?kTRUE:kFALSE;}
8e1f0465	114	Bool_t IsRemote() const {return fIsRemote;}
06a59280	115	Bool_t IsUsingDataSet() const {return TObject::TestBit(kUseDataSet);}
012e169c	116	void LoadBranch(const char *n) { if(fAutoBranchHandling) return; DoLoadBranch(n); }
c52c2132	117	void SetAnalysisType(EAliAnalysisExecMode mode) {fMode = mode;}
012e169c	118	void SetAutoBranchLoading(Bool_t b) { fAutoBranchHandling = b; }
06a59280	119	void SetCurrentEntry(Long64_t entry) {fCurrentEntry = entry;}
	120	void SetCollectSysInfoEach(Int_t nevents=0) {fNSysInfo = nevents;}
	121	static void SetCommonFileName(const char *name) {fgCommonFileName = name;}
	122	void SetDebugLevel(UInt_t level) {fDebug = level;}
	123	void SetDisableBranches(Bool_t disable=kTRUE) {TObject::SetBit(kDisableBranches,disable);}
	124	void SetExternalLoop(Bool_t flag) {TObject::SetBit(kExternalLoop,flag);}
61505f8b	125	void SetEventPool(AliVEventPool* const epool) {fEventPool = epool;}
	126	void SetGridHandler(AliAnalysisGrid * const handler) {fGridHandler = handler;}
	127	void SetInputEventHandler(AliVEventHandler* const handler);
	128	void SetMCtruthEventHandler(AliVEventHandler* const handler) {fMCtruthEventHandler = handler;}
06a59280	129	void SetNSysInfo(Long64_t nevents) {fNSysInfo = nevents;}
61505f8b	130	void SetOutputEventHandler(AliVEventHandler* const handler);
61505f8b	131	void SetSelector(AliAnalysisSelector * const sel) {fSelector = sel;}
06a59280	132	void SetSaveCanvases(Bool_t flag=kTRUE) {TObject::SetBit(kSaveCanvases,flag);}
90a4b3ee	133	void SetSkipTerminate(Bool_t flag) {TObject::SetBit(kSkipTerminate,flag);}
ff0d33e8	134	void SetUseProgressBar(Bool_t flag) {TObject::SetBit(kUseProgressBar,flag);}
06a59280	135	void SetSpecialOutputLocation(const char *loc) {fSpecialOutputLocation = loc;}
d3106602	136
	137	// Container handling
	138	AliAnalysisDataContainer CreateContainer(const char name, TClass *datatype,
c52c2132	139	EAliAnalysisContType type = kExchangeContainer,
c52c2132	140	const char *filename = NULL);
d3106602	141
	142	// Including tasks and getting them
	143	void AddTask(AliAnalysisTask *task);
	144	AliAnalysisTask GetTask(const char name) const;
	145
	146	// Connecting data containers to task inputs/outputs
	147	Bool_t ConnectInput(AliAnalysisTask *task, Int_t islot,
	148	AliAnalysisDataContainer *cont);
	149	Bool_t ConnectOutput(AliAnalysisTask *task, Int_t islot,
	150	AliAnalysisDataContainer *cont);
	151	// Garbage collection
	152	void CleanContainers();
	153
	154	// Analysis initialization and execution, status
	155	Bool_t InitAnalysis();
	156	Bool_t IsInitialized() const {return fInitOK;}
90a4b3ee	157	Bool_t IsExternalLoop() const {return TObject::TestBit(kExternalLoop);}
efd53803	158	Bool_t IsEventLoop() const {return TObject::TestBit(kEventLoop);}
90a4b3ee	159	Bool_t IsSkipTerminate() const {return TObject::TestBit(kSkipTerminate);}
d3106602	160	void ResetAnalysis();
d3106602	161	void ExecAnalysis(Option_t *option="");
d3106602	162	void PrintStatus(Option_t *option="all") const;
096b5a2e	163	void ProfileTask(const char name, const char option="VM") const;
096b5a2e	164	void ProfileTask(Int_t itop, const char *option="VM") const;
61505f8b	165	static void ProgressBar(const char opname, Long64_t current, Long64_t size, TStopwatch const watch=0, Bool_t last=kFALSE, Bool_t refresh=kFALSE);
923e2ca5	166	Bool_t ValidateOutputFiles() const;
d3106602	167
d3106602	168	protected:
981f2614	169	void ImportWrappers(TList *source);
efd53803	170	void SetEventLoop(Bool_t flag=kTRUE) {TObject::SetBit(kEventLoop,flag);}
012e169c	171	void DoLoadBranch(const char *name);
d3106602	172
5daf9fd2	173	private:
6bb2b24f	174	TTree *fTree; //! Input tree in case of TSelector model
54cff064	175	AliVEventHandler *fInputEventHandler; // Optional common input event handler
d2f1d9ef	176	AliVEventHandler *fOutputEventHandler; // Optional common output event handler
d2f1d9ef	177	AliVEventHandler *fMCtruthEventHandler; // Optional common MC Truth event handler
c2922515	178	AliVEventPool *fEventPool; // Event pool for mixing analysis
6bb2b24f	179	Long64_t fCurrentEntry; //! Current processed entry in the tree
8c0ab8e8	180	Long64_t fNSysInfo; // Event frequency for collecting system information
6bb2b24f	181	EAliAnalysisExecMode fMode; // Execution mode
6bb2b24f	182	Bool_t fInitOK; // Initialisation done
8e1f0465	183	Bool_t fIsRemote; //! Flag is set for remote analysis
6bb2b24f	184	UInt_t fDebug; // Debug level
13ef3bb0	185	TString fSpecialOutputLocation; // URL/path where the special outputs will be copied
6bb2b24f	186	TObjArray *fTasks; // List of analysis tasks
	187	TObjArray *fTopTasks; // List of top tasks
	188	TObjArray *fZombies; // List of zombie tasks
	189	TObjArray *fContainers; // List of all containers
	190	TObjArray *fInputs; // List of containers with input data
	191	TObjArray *fOutputs; // List of containers with results
ff07ec61	192	TObjArray *fParamCont; // List of containers with results
60a04972	193	AliAnalysisDataContainer *fCommonInput; // Common input container
60a04972	194	AliAnalysisDataContainer *fCommonOutput; // Common output container
8d7d3b59	195	AliAnalysisSelector *fSelector; //! Current selector
c57f56b7	196	AliAnalysisGrid *fGridHandler; //! Grid handler plugin
c9e39043	197	TString fExtraFiles; // List of extra files to be merged
012e169c	198	Bool_t fAutoBranchHandling; // def=kTRUE, turn off if you use LoadBranch
012e169c	199	THashTable fTable; // keep branch ptrs in case of manual branch loading
5daf9fd2	200
84fcd93f	201	static TString fgCommonFileName; //! Common output file name (not streamed)
c52c2132	202	static AliAnalysisManager *fgAnalysisManager; //! static pointer to object instance
8e1f0465	203	ClassDef(AliAnalysisManager,8) // Analysis manager class
d3106602	204	};
d3106602	205	#endif