[u/mrichter/AliRoot.git] / STEER / AliMCEventHandler.h

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


/* $Id$ */

//-------------------------------------------------------------------------
//                          Class AliMCEvent
// This class gives access to MC truth during the analysis.
// Monte Carlo truth is contained in the kinematics tree (produced particles) and 
// the tree of reference hits.
//      
// Origin: Andreas Morsch, CERN, andreas.morsch@cern.ch 
//-------------------------------------------------------------------------
#include "AliVEventHandler.h"
#include "AliHeader.h"
#include <TExMap.h>

class TFile;
class TTree;
class TList;

class TParticle;
class TString;
class TClonesArray;
class TDirectoryFile;

class AliMCEvent;


class AliMCEventHandler : public AliVEventHandler
{
public:
    AliMCEventHandler();
    AliMCEventHandler(const char* name, const char* title);
    virtual ~AliMCEventHandler();
    virtual void         SetOutputFileName(const char* /* fname */) {;}
    virtual const char*  GetOutputFileName() {return 0;}
    virtual void         SetInputPath(const char* fname); 
    virtual void         SetInputTree(TTree* /*tree*/) {;}
    virtual TString*     GetInputPath() const {return fPathName;}
    virtual Bool_t       Init(Option_t* opt);
    virtual Bool_t       InitIO(Option_t* opt) {return Init(opt);};
    virtual Bool_t       Init(TTree* /*tree*/, Option_t* /*opt*/) {return kTRUE;}
    virtual Bool_t       BeginEvent(Long64_t entry);
    virtual Bool_t       Notify() { return AliVEventHandler::Notify(); };
    virtual Bool_t       Notify(const char* path);
    virtual Bool_t       FinishEvent();
    virtual Bool_t       Terminate();
    virtual Bool_t       TerminateIO();
    virtual void         ResetIO();
    virtual Bool_t       GetEvent(Int_t iev);
    virtual void         SetReadTR(Bool_t flag) { fReadTR = flag; }
    virtual void         AddSubsidiaryHandler(AliMCEventHandler* handler);
    virtual void         SetNumberOfEventsInContainer(Int_t nev) {fEventsInContainer = nev;}
    //
    AliMCEvent* MCEvent() const {return fMCEvent;} 
    TTree*      TreeTR()  const {return fTreeTR;}
    TTree*      TreeK()   const {return fTreeK;}
    Int_t       GetParticleAndTR(Int_t i, TParticle*& particle, TClonesArray*& trefs);
    void        DrawCheck(Int_t i, Int_t search=0);
    Bool_t      InitOk() {return fInitOk;}
    // Label manipulation
    void   SelectParticle(Int_t i);
    Bool_t IsParticleSelected(Int_t i);
    void   CreateLabelMap();
    Int_t  GetNewLabel(Int_t i);

private:
    Bool_t      OpenFile(Int_t i);
    void  VerifySelectedParticles();
    AliMCEventHandler(const AliMCEventHandler& handler);             
    AliMCEventHandler& operator=(const AliMCEventHandler& handler);  
private:
    AliMCEvent            *fMCEvent;            //! MC Event
    TFile                 *fFileE;              //! File with TreeE
    TFile                 *fFileK;              //! File with TreeK
    TFile                 *fFileTR;             //! File with TreeTR
    TTree                 *fTreeE;              //! TreeE  (Event Headers)
    TTree                 *fTreeK;              //! TreeK  (kinematics tree)
    TTree                 *fTreeTR;             //! TreeTR (track references tree)
    TDirectoryFile        *fDirK;               //! Directory for Kine Tree
    TDirectoryFile        *fDirTR;              //! Directory for TR Tree
    TExMap                 fParticleSelected;   //! List of selected MC particles for t
    TExMap                 fLabelMap;           //! Stores the Map of MC (ESDLabel,AODlabel)  
    Int_t                  fNEvent;             //! Number of events
    Int_t                  fEvent;              //! Current event
    TString               *fPathName;           //! Input file path 
    const Char_t          *fExtension;          //! File name extension 
    Int_t                  fFileNumber;         //! Input file number
    Int_t                  fEventsPerFile;      //! Number of events per file
    Bool_t                 fReadTR;             // determines if TR shall be read
    Bool_t                 fInitOk;             // Initialization ok
    TList                 *fSubsidiaryHandlers; //! List of subsidiary MC handlers (for example for Background)
    Int_t                  fEventsInContainer;  //! Number of events in container class 
    ClassDef(AliMCEventHandler,1)  //MC Truth EventHandler class
};
#endif
Commit	Line	Data
969c7896	1	// -- mode: C++ --
5fe09262	2	#ifndef ALIMCEVENTHANDLER_H
5fe09262	3	#define ALIMCEVENTHANDLER_H
033693d3	4	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	5	* See cxx source for full Copyright notice */
	6
	7
	8	/* $Id$ */
	9
	10	//-------------------------------------------------------------------------
	11	// Class AliMCEvent
	12	// This class gives access to MC truth during the analysis.
415d9f5c	13	// Monte Carlo truth is contained in the kinematics tree (produced particles) and
033693d3	14	// the tree of reference hits.
	15	//
	16	// Origin: Andreas Morsch, CERN, andreas.morsch@cern.ch
	17	//-------------------------------------------------------------------------
d2f1d9ef	18	#include "AliVEventHandler.h"
47b95647	19	#include "AliHeader.h"
da97a08a	20	#include <TExMap.h>
da97a08a	21
033693d3	22	class TFile;
033693d3	23	class TTree;
93836e1b	24	class TList;
93836e1b	25
033693d3	26	class TParticle;
0a05cd41	27	class TString;
033693d3	28	class TClonesArray;
5efedd31	29	class TDirectoryFile;
5efedd31	30
415d9f5c	31	class AliMCEvent;
415d9f5c	32
033693d3	33
033693d3	34
276be602	35	class AliMCEventHandler : public AliVEventHandler
033693d3	36	{
033693d3	37	public:
5fe09262	38	AliMCEventHandler();
	39	AliMCEventHandler(const char* name, const char* title);
	40	virtual ~AliMCEventHandler();
0931e76a	41	virtual void SetOutputFileName(const char* /* fname */) {;}
	42	virtual const char* GetOutputFileName() {return 0;}
	43	virtual void SetInputPath(const char* fname);
933fd60f	44	virtual void SetInputTree(TTree* /tree/) {;}
d0916b2f	45	virtual TString* GetInputPath() const {return fPathName;}
300d5701	46	virtual Bool_t Init(Option_t* opt);
	47	virtual Bool_t InitIO(Option_t* opt) {return Init(opt);};
	48	virtual Bool_t Init(TTree* /tree/, Option_t* /opt/) {return kTRUE;}
ed97dc98	49	virtual Bool_t BeginEvent(Long64_t entry);
24e577e5	50	virtual Bool_t Notify() { return AliVEventHandler::Notify(); };
890126ab	51	virtual Bool_t Notify(const char* path);
033693d3	52	virtual Bool_t FinishEvent();
	53	virtual Bool_t Terminate();
	54	virtual Bool_t TerminateIO();
5fe09262	55	virtual void ResetIO();
9aea8469	56	virtual Bool_t GetEvent(Int_t iev);
969c7896	57	virtual void SetReadTR(Bool_t flag) { fReadTR = flag; }
93836e1b	58	virtual void AddSubsidiaryHandler(AliMCEventHandler* handler);
93836e1b	59	virtual void SetNumberOfEventsInContainer(Int_t nev) {fEventsInContainer = nev;}
033693d3	60	//
2bb794ba	61	AliMCEvent* MCEvent() const {return fMCEvent;}
d0916b2f	62	TTree* TreeTR() const {return fTreeTR;}
d0916b2f	63	TTree* TreeK() const {return fTreeK;}
415d9f5c	64	Int_t GetParticleAndTR(Int_t i, TParticle& particle, TClonesArray& trefs);
415d9f5c	65	void DrawCheck(Int_t i, Int_t search=0);
b72029a3	66	Bool_t InitOk() {return fInitOk;}
	67	// Label manipulation
	68	void SelectParticle(Int_t i);
	69	Bool_t IsParticleSelected(Int_t i);
	70	void CreateLabelMap();
	71	Int_t GetNewLabel(Int_t i);
da97a08a	72
5fe09262	73	private:
415d9f5c	74	Bool_t OpenFile(Int_t i);
da97a08a	75	void VerifySelectedParticles();
87d395bd	76	AliMCEventHandler(const AliMCEventHandler& handler);
87d395bd	77	AliMCEventHandler& operator=(const AliMCEventHandler& handler);
033693d3	78	private:
93836e1b	79	AliMCEvent *fMCEvent; //! MC Event
	80	TFile *fFileE; //! File with TreeE
	81	TFile *fFileK; //! File with TreeK
	82	TFile *fFileTR; //! File with TreeTR
	83	TTree *fTreeE; //! TreeE (Event Headers)
	84	TTree *fTreeK; //! TreeK (kinematics tree)
	85	TTree *fTreeTR; //! TreeTR (track references tree)
	86	TDirectoryFile *fDirK; //! Directory for Kine Tree
	87	TDirectoryFile *fDirTR; //! Directory for TR Tree
	88	TExMap fParticleSelected; //! List of selected MC particles for t
	89	TExMap fLabelMap; //! Stores the Map of MC (ESDLabel,AODlabel)
	90	Int_t fNEvent; //! Number of events
	91	Int_t fEvent; //! Current event
	92	TString *fPathName; //! Input file path
	93	const Char_t *fExtension; //! File name extension
	94	Int_t fFileNumber; //! Input file number
	95	Int_t fEventsPerFile; //! Number of events per file
	96	Bool_t fReadTR; // determines if TR shall be read
	97	Bool_t fInitOk; // Initialization ok
	98	TList *fSubsidiaryHandlers; //! List of subsidiary MC handlers (for example for Background)
	99	Int_t fEventsInContainer; //! Number of events in container class
969c7896	100	ClassDef(AliMCEventHandler,1) //MC Truth EventHandler class
033693d3	101	};
	102	#endif
	103