4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
10 //-------------------------------------------------------------------------
13 // Origin: Andreas.Morsch, CERN, andreas.morsch@cern.ch
14 //-------------------------------------------------------------------------
18 #include <TRefArray.h>
19 #include <TClonesArray.h>
21 #include <AliVEvent.h>
22 #include "AliVHeader.h"
23 #include "AliVParticle.h"
24 #include "AliVVertex.h"
25 #include "AliMCParticle.h"
31 class AliGenEventHeader;
36 class AliMCEvent : public AliVEvent {
41 virtual ~AliMCEvent() {;}
42 AliMCEvent(const AliMCEvent& mcEvnt);
43 AliMCEvent& operator=(const AliMCEvent& mcEvnt);
45 // Methods implementing the interface
48 virtual void AddObject(TObject* /*obj*/) {;}
49 virtual TObject* FindListObject(const char */*name*/) const {return 0;}
50 virtual TList* GetList() const {return 0;}
51 virtual void CreateStdContent() {;}
52 virtual void GetStdContent() {;}
53 virtual void ReadFromTree(TTree * /*tree*/, Option_t* /*opt*/) {;}
54 virtual void WriteToTree(TTree* /*tree*/) const {;}
56 virtual void SetStdNames() {;}
57 virtual void Print(Option_t */*option=""*/) const {;}
58 virtual void PreReadAll();
59 virtual void Reset() {;}
62 virtual AliVHeader* GetHeader() const {return 0;}
64 // Delegated methods for fESDRun or AODHeader
66 virtual void SetRunNumber(Int_t /*n*/) {;}
67 virtual void SetPeriodNumber(UInt_t /*n*/) {;}
68 virtual void SetMagneticField(Double_t /*mf*/) {;}
71 virtual Int_t GetRunNumber() const {return 0;}
72 virtual UInt_t GetPeriodNumber() const {return 0;}
73 virtual Double_t GetMagneticField() const {return 0.;}
76 virtual void SetOrbitNumber(UInt_t /*n*/) {;}
77 virtual void SetBunchCrossNumber(UShort_t /*n*/) {;}
78 virtual void SetEventType(UInt_t /*eventType*/) {;}
79 virtual void SetTriggerMask(ULong64_t /*n*/) {;}
80 virtual void SetTriggerCluster(UChar_t /*n*/) {;}
82 virtual UInt_t GetOrbitNumber() const {return 0;}
83 virtual UShort_t GetBunchCrossNumber() const {return 0;}
85 virtual UInt_t GetEventType() const {return 0;}
87 virtual ULong64_t GetTriggerMask() const {return 0;}
88 virtual UChar_t GetTriggerCluster() const {return 0;}
89 virtual TString GetFiredTriggerClasses()const {return ("");}
90 virtual Double_t GetZDCN1Energy() const {return 0.;}
91 virtual Double_t GetZDCP1Energy() const {return 0.;}
92 virtual Double_t GetZDCN2Energy() const {return 0.;}
93 virtual Double_t GetZDCP2Energy() const {return 0.;}
94 virtual Double_t GetZDCEMEnergy(Int_t /*i*/)
97 virtual AliVParticle *GetTrack(Int_t i) const;
98 virtual Int_t GetNumberOfTracks() const {return fNparticles;}
99 virtual Int_t GetNumberOfV0s() const {return -1;}
100 virtual Int_t GetNumberOfCascades() const {return -1;}
102 virtual const AliVVertex *GetPrimaryVertex() const;
105 // MC Specific methods
108 AliStack* Stack() {return fStack;}
109 AliHeader* Header() {return fHeader;}
110 AliGenEventHeader* GenEventHeader() const;
112 virtual void ConnectTreeE (TTree* tree);
113 virtual void ConnectTreeK (TTree* tree);
114 virtual void ConnectTreeTR(TTree* tree);
115 virtual void Clean();
116 virtual void InitEvent();
117 virtual void FinishEvent();
118 virtual Int_t GetParticleAndTR(Int_t i, TParticle*& particle, TClonesArray*& trefs);
119 virtual void DrawCheck(Int_t i, Int_t search);
120 virtual void AddSubsidiaryEvent(AliMCEvent* event);
121 virtual Int_t GetNumberOfPrimaries() {return fNprimaries;}
122 virtual Int_t GetPrimaryOffset() const {return fPrimaryOffset;}
123 virtual Int_t GetSecondaryOffset() const {return fSecondaryOffset;}
124 virtual void SetPrimaryOffset(Int_t ioff) {fPrimaryOffset = ioff;}
125 virtual void SetSecondaryOffset(Int_t ioff) {fSecondaryOffset = ioff;}
126 virtual Bool_t IsPhysicalPrimary(Int_t i) const;
128 virtual Int_t BgLabelToIndex(Int_t label);
129 static Int_t BgLabelOffset() {return fgkBgLabelOffset;}
130 virtual Bool_t IsFromBGEvent(Int_t index);
131 virtual Bool_t IsSecondaryFromWeakDecay(Int_t index);
132 virtual Bool_t IsSecondaryFromMaterial(Int_t index);
133 // External particle array
134 virtual void SetParticleArray(TClonesArray* mcParticles)
135 {fMCParticles = mcParticles; fNparticles = fMCParticles->GetEntries(); fExternal = kTRUE;}
137 virtual AliGenEventHeader *FindHeader(Int_t ipart);
138 //Following needed only for mixed event
139 virtual Int_t EventIndex(Int_t) const {return 0;}
140 virtual Int_t EventIndexForCaloCluster(Int_t) const {return 0;}
141 virtual Int_t EventIndexForPHOSCell(Int_t) const {return 0;}
142 virtual Int_t EventIndexForEMCALCell(Int_t) const {return 0;}
143 AliCentrality* GetCentrality() {return 0;}
144 AliEventplane* GetEventplane() {return 0;}
146 virtual AliVVZERO *GetVZEROData() const {return 0;}
147 virtual AliVZDC *GetZDCData() const {return 0;}
149 // #ifdef MFT_UPGRADE
150 // virtual AliVMFT *GetMFTData() const {return 0;}
152 virtual AliVMFT *GetMFTData() const {return 0;} // AU
155 virtual void ReorderAndExpandTreeTR();
156 virtual Int_t FindIndexAndEvent(Int_t oldidx, AliMCEvent*& event) const;
159 // Stanndard implementation for ESD production
160 AliStack *fStack; // Current pointer to stack
161 TClonesArray *fMCParticles; // Pointer to list of particles
162 TObjArray *fMCParticleMap; // Map of MC Particles
163 AliHeader *fHeader; // Current pointer to header
164 TClonesArray *fTRBuffer; // Track reference buffer
165 TClonesArray *fTrackReferences; // Array of track references
166 TTree *fTreeTR; // Pointer to Track Reference Tree
167 TTree *fTmpTreeTR; // Temporary tree TR to read old format
168 TFile *fTmpFileTR; // Temporary file with TreeTR to read old format
169 Int_t fNprimaries; // Number of primaries
170 Int_t fNparticles; // Number of particles
171 TList *fSubsidiaryEvents; // List of possible subsidiary events (for example merged underlying event)
172 Int_t fPrimaryOffset; // Offset for primaries
173 Int_t fSecondaryOffset; // Offset for secondaries
174 Bool_t fExternal; // True if external particle array
175 static Int_t fgkBgLabelOffset; // Standard branch name
176 mutable AliVVertex* fVertex; // MC Vertex
177 Int_t fNBG; //! Background particles in current event
179 ClassDef(AliMCEvent, 2) // AliVEvent realisation for MC data