]> git.uio.no Git - u/mrichter/AliRoot.git/blob - STEER/AliESDEvent.h
A getter for the number of tracks composing the vertex
[u/mrichter/AliRoot.git] / STEER / AliESDEvent.h
1 // -*- mode: C++ -*- 
2 #ifndef ALIESDEVENT_H
3 #define ALIESDEVENT_H
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 AliESDEvent
12 //   This is the class to deal with during the physics analysis of data.
13 //   It also ensures the backward compatibility with the old ESD format.
14 //      
15 // Origin: Christian Klein-Boesing, CERN, Christian.Klein-Boesing@cern.ch 
16 //-------------------------------------------------------------------------
17
18 #include <TClonesArray.h>
19 #include <TObject.h>
20 #include <TTree.h>
21 #include <TArrayF.h>
22
23 class TList;
24
25 #include "AliVEvent.h"
26 // some includes for delegated methods
27 #include "AliESDCaloTrigger.h"
28 #include "AliESDRun.h"
29 #include "AliESDHeader.h"
30 #include "AliESDTZERO.h"
31 #include "AliESDZDC.h"
32
33 // AliESDtrack has to be included so that the compiler 
34 // knows its inheritance tree (= that it is a AliVParticle).
35 #include "AliESDtrack.h"
36
37 class AliESDfriend;
38 class AliESDVZERO;
39 class AliESDHLTtrack;
40 class AliESDVertex;
41 class AliESDPmdTrack;
42 class AliESDFMD;
43 class AliESDkink;
44 class AliESDCaloCluster;
45 class AliESDCaloCells;
46 class AliESDv0;
47 class AliMultiplicity;
48 class AliRawDataErrorLog;
49 class AliESDRun;
50 class AliESDTrdTrack;
51 class AliESDMuonTrack;
52 class AliESD;
53 class AliESDcascade;
54 class TRefArray;
55
56 class AliESDEvent : public AliVEvent {
57 public:
58
59
60   enum ESDListIndex   {kESDRun,
61                        kHeader,
62                        kESDZDC,
63                        kESDFMD,
64                        kESDVZERO,
65                        kESDTZERO,
66                        kSPDVertex,
67                        kPrimaryVertex,
68                        kSPDMult,
69                        kPHOSTrigger,
70                        kEMCALTrigger,
71                        kTracks,
72                        kMuonTracks,
73                        kPmdTracks,
74                        kTrdTracks,
75                        kV0s,
76                        kCascades,
77                        kKinks,
78                        kCaloClusters,
79                        kEMCALCells,
80                        kPHOSCells,
81                        kErrorLogs,
82                        kESDListN
83   };
84
85   AliESDEvent();
86   virtual ~AliESDEvent(); 
87
88
89   // RUN
90   // move this to the UserData!!!
91   const AliESDRun*    GetESDRun() const {return fESDRun;}
92
93   // Delegated methods for fESDRun
94   void     SetRunNumber(Int_t n) {fESDRun->SetRunNumber(n);}
95   Int_t    GetRunNumber() const {return fESDRun->GetRunNumber();}
96   void     SetPeriodNumber(UInt_t n){fESDRun->SetPeriodNumber(n);}
97   UInt_t   GetPeriodNumber() const {return fESDRun->GetPeriodNumber();}
98   void     SetMagneticField(Double_t mf){fESDRun->SetMagneticField(mf);}
99   Double_t GetMagneticField() const {return fESDRun->GetMagneticField();}
100   void     SetDiamond(const AliESDVertex *vertex) { fESDRun->SetDiamond(vertex);}
101   Double_t  GetDiamondX() const {return fESDRun->GetDiamondX();}
102   Double_t  GetDiamondY() const {return fESDRun->GetDiamondY();}
103   Double_t  GetSigma2DiamondX() const {return  fESDRun->GetSigma2DiamondX();}
104   Double_t  GetSigma2DiamondY() const {return  fESDRun->GetSigma2DiamondY();}
105   void      GetDiamondCovXY(Float_t cov[3]) const {fESDRun->GetDiamondCovXY(cov);}   
106   
107
108   // HEADER
109   AliESDHeader* GetHeader() const {return fHeader;}
110
111   // Delegated methods for fHeader
112   void      SetTriggerMask(ULong64_t n) {fHeader->SetTriggerMask(n);}
113   void      SetOrbitNumber(UInt_t n) {fHeader->SetOrbitNumber(n);}
114   void      SetTimeStamp(UInt_t timeStamp){fHeader->SetTimeStamp(timeStamp);}
115   void      SetEventType(UInt_t eventType){fHeader->SetEventType(eventType);}
116   void      SetEventNumberInFile(Int_t n) {fHeader->SetEventNumberInFile(n);}
117   //  void      SetRunNumber(Int_t n) {fHeader->SetRunNumber(n);}
118   void      SetBunchCrossNumber(UShort_t n) {fHeader->SetBunchCrossNumber(n);}
119   void      SetTriggerCluster(UChar_t n) {fHeader->SetTriggerCluster(n);}
120   
121   ULong64_t GetTriggerMask() const {return fHeader->GetTriggerMask();}
122   UInt_t    GetOrbitNumber() const {return fHeader->GetOrbitNumber();}
123   UInt_t    GetTimeStamp()  const { return fHeader->GetTimeStamp();}
124   UInt_t    GetEventType()  const { return fHeader->GetEventType();}
125   Int_t     GetEventNumberInFile() const {return fHeader->GetEventNumberInFile();}
126   UShort_t  GetBunchCrossNumber() const {return fHeader->GetBunchCrossNumber();}
127   UChar_t   GetTriggerCluster() const {return fHeader->GetTriggerCluster();}
128
129   // ZDC CKB: put this in the header?
130   AliESDZDC*    GetESDZDC() const {return fESDZDC;}
131
132   // Delegated methods for fESDZDC
133   Double_t GetZDCN1Energy() const {return fESDZDC->GetZDCN1Energy();}
134   Double_t GetZDCP1Energy() const {return fESDZDC->GetZDCP1Energy();}
135   Double_t GetZDCN2Energy() const {return fESDZDC->GetZDCN2Energy();}
136   Double_t GetZDCP2Energy() const {return fESDZDC->GetZDCP2Energy();}
137   Double_t GetZDCEMEnergy(Int_t i=0) const {return fESDZDC->GetZDCEMEnergy(i);}
138   Int_t   GetZDCParticipants() const {return fESDZDC->GetZDCParticipants();}
139   void    SetZDC(Float_t n1Energy, Float_t p1Energy, Float_t em1Energy, Float_t em2Energy,
140                  Float_t n2Energy, Float_t p2Energy, Int_t participants)
141   {fESDZDC->SetZDC(n1Energy, p1Energy, em1Energy, em2Energy, n2Energy, p2Energy, participants);}
142
143
144   // FMD
145   void SetFMDData(AliESDFMD * obj);
146   AliESDFMD *GetFMDData(){ return fESDFMD; }
147
148
149   // TZERO CKB: put this in the header?
150   const AliESDTZERO*    GetESDTZERO() const {return fESDTZERO;}
151   // delegetated methods for fESDTZERO
152
153   Double_t GetT0zVertex() const {return fESDTZERO->GetT0zVertex();}
154   void SetT0zVertex(Float_t z) {fESDTZERO->SetT0zVertex(z);}
155   Double_t GetT0() const {return fESDTZERO->GetT0();}
156   void SetT0(Float_t timeStart) {fESDTZERO->SetT0(timeStart);}
157   const Double_t * GetT0time() const {return fESDTZERO->GetT0time();}
158   void SetT0time(Float_t time[24]) {fESDTZERO->SetT0time(time);}
159   const Double_t * GetT0amplitude() const {return fESDTZERO->GetT0amplitude();}
160   void SetT0amplitude(Float_t amp[24]){fESDTZERO->SetT0amplitude(amp);}
161
162   // VZERO 
163   AliESDVZERO *GetVZEROData() const { return fESDVZERO; }
164   void SetVZEROData(AliESDVZERO * obj);
165
166
167   void SetESDfriend(const AliESDfriend *f);
168   void GetESDfriend(AliESDfriend *f) const;
169
170
171
172   void SetVertex(const AliESDVertex *vertex); 
173   const AliESDVertex *GetVertex() const {return fSPDVertex;}
174
175   void SetPrimaryVertex(const AliESDVertex *vertex);
176   const AliESDVertex *GetPrimaryVertex() const {return fPrimaryVertex;}
177
178   void SetMultiplicity(const AliMultiplicity *mul);
179
180   const AliMultiplicity *GetMultiplicity() const {return fSPDMult;}
181
182
183   Bool_t Clean(Float_t *cleanPars);
184   Bool_t RemoveKink(Int_t i);
185   Bool_t RemoveV0(Int_t i);
186   Bool_t RemoveTrack(Int_t i);
187
188   AliESDtrack *GetTrack(Int_t i) const {
189     return (AliESDtrack *)fTracks->UncheckedAt(i);
190   }
191   Int_t  AddTrack(const AliESDtrack *t);
192
193   
194   AliESDHLTtrack *GetHLTConfMapTrack(Int_t /*i*/) const {
195     //    return (AliESDHLTtrack *)fHLTConfMapTracks->UncheckedAt(i);
196     return 0;
197   }
198   void AddHLTConfMapTrack(const AliESDHLTtrack */*t*/) {
199     printf("ESD:: AddHLTConfMapTrack do nothing \n");
200     //    TClonesArray &fhlt = *fHLTConfMapTracks;
201     //  new(fhlt[fHLTConfMapTracks->GetEntriesFast()]) AliESDHLTtrack(*t);
202   }
203   
204
205   AliESDHLTtrack *GetHLTHoughTrack(Int_t /*i*/) const {
206     //    return (AliESDHLTtrack *)fHLTHoughTracks->UncheckedAt(i);
207     return 0;
208   }
209   void AddHLTHoughTrack(const AliESDHLTtrack */*t*/) {
210     printf("ESD:: AddHLTHoughTrack do nothing \n");
211     //    TClonesArray &fhlt = *fHLTHoughTracks;
212     //     new(fhlt[fHLTHoughTracks->GetEntriesFast()]) AliESDHLTtrack(*t);
213   }
214   
215   AliESDMuonTrack *GetMuonTrack(Int_t i) const {
216     return (AliESDMuonTrack *)fMuonTracks->UncheckedAt(i);
217   }
218
219   void AddMuonTrack(const AliESDMuonTrack *t);
220
221   AliESDPmdTrack *GetPmdTrack(Int_t i) const {
222     return (AliESDPmdTrack *)fPmdTracks->UncheckedAt(i);
223   }
224
225   void AddPmdTrack(const AliESDPmdTrack *t);
226
227
228   AliESDTrdTrack *GetTrdTrack(Int_t i) const {
229     return (AliESDTrdTrack *)fTrdTracks->UncheckedAt(i);
230   }
231
232   
233   void AddTrdTrack(const AliESDTrdTrack *t);
234
235   AliESDv0 *GetV0(Int_t i) const {
236     return (AliESDv0*)fV0s->UncheckedAt(i);
237   }
238   Int_t AddV0(const AliESDv0 *v);
239
240   AliESDcascade *GetCascade(Int_t i) const {
241     return (AliESDcascade *)fCascades->UncheckedAt(i);
242   }
243
244   void AddCascade(const AliESDcascade *c);
245
246   AliESDkink *GetKink(Int_t i) const {
247     return (AliESDkink *)fKinks->UncheckedAt(i);
248   }
249   Int_t AddKink(const AliESDkink *c);
250
251   AliESDCaloCluster *GetCaloCluster(Int_t i) const {
252     return (AliESDCaloCluster *)fCaloClusters->UncheckedAt(i);
253   }
254
255   Int_t AddCaloCluster(const AliESDCaloCluster *c);
256
257   AliESDCaloCells *GetEMCALCells() const {return fEMCALCells; }  
258   AliESDCaloCells *GetPHOSCells() const {return fPHOSCells; }  
259
260   AliRawDataErrorLog *GetErrorLog(Int_t i) const {
261     return (AliRawDataErrorLog *)fErrorLogs->UncheckedAt(i);
262   }
263   void  AddRawDataErrorLog(const AliRawDataErrorLog *log);
264
265   Int_t GetNumberOfErrorLogs()   const {return fErrorLogs->GetEntriesFast();}
266
267     
268   void AddPHOSTriggerPosition(TArrayF array)   { fPHOSTrigger->AddTriggerPosition(array); }
269   void AddPHOSTriggerAmplitudes(TArrayF array) { fPHOSTrigger->AddTriggerAmplitudes(array);}
270   void AddEMCALTriggerPosition(TArrayF array)  { fEMCALTrigger->AddTriggerPosition(array); }
271   void AddEMCALTriggerAmplitudes(TArrayF array){ fEMCALTrigger->AddTriggerAmplitudes(array); }
272
273   Int_t GetNumberOfTracks()     const {return fTracks->GetEntriesFast();}
274   Int_t GetNumberOfHLTConfMapTracks()     const {return 0;} 
275   // fHLTConfMapTracks->GetEntriesFast();}
276   Int_t GetNumberOfHLTHoughTracks()     const {return  0;  }
277   //  fHLTHoughTracks->GetEntriesFast();  }
278
279   Int_t GetNumberOfMuonTracks() const {return fMuonTracks->GetEntriesFast();}
280   Int_t GetNumberOfPmdTracks() const {return fPmdTracks->GetEntriesFast();}
281   Int_t GetNumberOfTrdTracks() const {return fTrdTracks->GetEntriesFast();}
282   Int_t GetNumberOfV0s()      const {return fV0s->GetEntriesFast();}
283   Int_t GetNumberOfCascades() const {return fCascades->GetEntriesFast();}
284   Int_t GetNumberOfKinks() const {return fKinks->GetEntriesFast();}
285   
286   Int_t GetEMCALClusters(TRefArray *clusters) const;
287   Int_t GetPHOSClusters(TRefArray *clusters) const;
288   Int_t GetNumberOfCaloClusters() const {return fCaloClusters->GetEntriesFast();}
289
290   
291   // Remove this stuff CKB?
292   //---------------------------------------------------
293   Int_t GetNumberOfEMCALClusters() const {return fEMCALClusters;}
294   void  SetNumberOfEMCALClusters(Int_t clus) {fEMCALClusters = clus;}
295   Int_t GetFirstEMCALCluster() const {return fFirstEMCALCluster;}
296   void  SetFirstEMCALCluster(Int_t index) {fFirstEMCALCluster = index;}
297  
298   Int_t GetNumberOfPHOSClusters() const {return fPHOSClusters;}
299   void  SetNumberOfPHOSClusters(Int_t part) { fPHOSClusters = part ; }
300   void  SetFirstPHOSCluster(Int_t index) { fFirstPHOSCluster = index ; } 
301   Int_t GetFirstPHOSCluster() const  { return fFirstPHOSCluster ; }
302   //-------------------------------------------------------
303
304   TArrayF *GetEMCALTriggerPosition() const {return  fEMCALTrigger->GetTriggerPosition();}
305   TArrayF *GetEMCALTriggerAmplitudes() const {return  fEMCALTrigger->GetTriggerAmplitudes();}
306   TArrayF *GetPHOSTriggerPosition() const {return  fPHOSTrigger->GetTriggerPosition();}
307   TArrayF *GetPHOSTriggerAmplitudes() const {return  fPHOSTrigger->GetTriggerAmplitudes();}
308
309   void ResetV0s() { fV0s->Clear(); }
310   void ResetCascades() { fCascades->Clear(); }
311   void Reset();
312
313   void  Print(Option_t *option="") const;
314
315   void AddObject(TObject* obj);
316   void ReadFromTree(TTree *tree);
317   TObject* FindListObject(const char *name);
318   AliESD *GetAliESDOld(){return fESDOld;}
319   const void WriteToTree(TTree* tree) const;
320   void GetStdContent();
321   void ResetStdContent();
322   void CreateStdContent();
323   void SetStdNames();
324   void CopyFromOldESD();
325   TList* GetList() const {return fESDObjects;}
326
327 protected:
328   AliESDEvent(const AliESDEvent&);
329   AliESDEvent &operator=(const AliESDEvent& source);
330
331
332   TList *fESDObjects;             // List of esd Objects
333
334   AliESDRun       *fESDRun;           //! Run information tmp put in the Userdata
335   AliESDHeader    *fHeader;           //! ESD Event Header
336   AliESDZDC       *fESDZDC;           //! ZDC information
337   AliESDFMD       *fESDFMD;           //! FMD object containing rough multiplicity
338   AliESDVZERO     *fESDVZERO;         //! VZERO object containing rough multiplicity
339   AliESDTZERO     *fESDTZERO;         //! TZEROObject
340   AliESDVertex    *fSPDVertex;        //! Primary vertex estimated by the SPD
341   AliESDVertex    *fPrimaryVertex;    //! Primary vertex estimated using ESD tracks
342   AliMultiplicity *fSPDMult;          //! SPD tracklet multiplicity
343   AliESDCaloTrigger* fPHOSTrigger;     //! PHOS Trigger information
344   AliESDCaloTrigger* fEMCALTrigger;    //! PHOS Trigger information
345
346   TClonesArray *fTracks;           //! ESD tracks 
347   TClonesArray *fMuonTracks;       //! MUON ESD tracks
348   TClonesArray *fPmdTracks;        //! PMD ESD tracks
349   TClonesArray *fTrdTracks;        //! TRD ESD tracks (triggered)
350   TClonesArray *fV0s;              //! V0 vertices
351   TClonesArray *fCascades;         //! Cascade vertices
352   TClonesArray *fKinks;            //! Kinks
353   TClonesArray *fCaloClusters;     //! Calorimeter clusters for PHOS/EMCAL
354   AliESDCaloCells *fEMCALCells;     //! EMCAL cell info
355   AliESDCaloCells *fPHOSCells;     //! PHOS cell info
356   TClonesArray *fErrorLogs;        //! Raw-data reading error messages
357  
358
359
360   AliESD       *fESDOld;           //! Old esd Structure
361   AliESDfriend *fESDFriendOld;     //! Old friend esd Structure
362   Bool_t    fConnected;            //! flag if leaves are alreday connected
363
364   static const char* fgkESDListName[kESDListN]; //!
365
366   // Remove this stuff CKB
367   Int_t        fEMCALClusters;   // Number of EMCAL clusters (subset of caloclusters)
368   Int_t        fFirstEMCALCluster; // First EMCAL cluster in the fCaloClusters list 
369
370   Int_t        fPHOSClusters;     // Number of PHOS clusters (subset of caloclusters)
371   Int_t        fFirstPHOSCluster; // First PHOS cluster in the fCaloClusters list 
372
373   ClassDef(AliESDEvent,6)  //ESDEvent class 
374 };
375 #endif 
376