4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
10 //-------------------------------------------------------------------------
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.
15 // Origin: Christian Klein-Boesing, CERN, Christian.Klein-Boesing@cern.ch
16 //-------------------------------------------------------------------------
18 #include <TClonesArray.h>
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 #include "AliESDACORDE.h"
34 // AliESDtrack has to be included so that the compiler
35 // knows its inheritance tree (= that it is a AliVParticle).
36 #include "AliESDtrack.h"
37 // same for AliESDVertex (which is a AliVVertex)
38 #include "AliESDVertex.h"
47 class AliESDCaloCluster;
48 class AliESDCaloCells;
50 class AliMultiplicity;
51 class AliRawDataErrorLog;
54 class AliESDMuonTrack;
60 class AliESDEvent : public AliVEvent {
64 enum ESDListIndex {kESDRun,
94 virtual ~AliESDEvent();
95 AliESDEvent &operator=(const AliESDEvent& source); // or make private and use only copy?
96 virtual void Copy(TObject& obj) const;
99 // move this to the UserData!!!
100 const AliESDRun* GetESDRun() const {return fESDRun;}
102 // Delegated methods for fESDRun
103 void SetRunNumber(Int_t n) {fESDRun->SetRunNumber(n);}
104 Int_t GetRunNumber() const {return fESDRun->GetRunNumber();}
105 void SetPeriodNumber(UInt_t n){fESDRun->SetPeriodNumber(n);}
106 UInt_t GetPeriodNumber() const {return fESDRun->GetPeriodNumber();}
107 void SetMagneticField(Double_t mf){fESDRun->SetMagneticField(mf);}
108 Double_t GetMagneticField() const {return fESDRun->GetMagneticField();}
109 void SetDiamond(const AliESDVertex *vertex) { fESDRun->SetDiamond(vertex);}
110 Double_t GetDiamondX() const {return fESDRun->GetDiamondX();}
111 Double_t GetDiamondY() const {return fESDRun->GetDiamondY();}
112 Double_t GetSigma2DiamondX() const {return fESDRun->GetSigma2DiamondX();}
113 Double_t GetSigma2DiamondY() const {return fESDRun->GetSigma2DiamondY();}
114 void GetDiamondCovXY(Float_t cov[3]) const {fESDRun->GetDiamondCovXY(cov);}
115 void SetTriggerClass(const char*name, Int_t index) {fESDRun->SetTriggerClass(name,index);}
116 void SetPHOSMatrix(TGeoHMatrix*matrix, Int_t i) {fESDRun->SetPHOSMatrix(matrix,i);}
117 const TGeoHMatrix* GetPHOSMatrix(Int_t i) const {return fESDRun->GetPHOSMatrix(i);}
121 AliESDHeader* GetHeader() const {return fHeader;}
123 // Delegated methods for fHeader
124 void SetTriggerMask(ULong64_t n) {fHeader->SetTriggerMask(n);}
125 void SetOrbitNumber(UInt_t n) {fHeader->SetOrbitNumber(n);}
126 void SetTimeStamp(UInt_t timeStamp){fHeader->SetTimeStamp(timeStamp);}
127 void SetEventType(UInt_t eventType){fHeader->SetEventType(eventType);}
128 void SetEventNumberInFile(Int_t n) {fHeader->SetEventNumberInFile(n);}
129 // void SetRunNumber(Int_t n) {fHeader->SetRunNumber(n);}
130 void SetBunchCrossNumber(UShort_t n) {fHeader->SetBunchCrossNumber(n);}
131 void SetTriggerCluster(UChar_t n) {fHeader->SetTriggerCluster(n);}
133 ULong64_t GetTriggerMask() const {return fHeader->GetTriggerMask();}
134 TString GetFiredTriggerClasses() const {return fESDRun->GetFiredTriggerClasses(fHeader->GetTriggerMask());}
135 Bool_t IsTriggerClassFired(const char *name) const {return fESDRun->IsTriggerClassFired(fHeader->GetTriggerMask(),name);}
136 TObject* GetHLTTriggerDecision() const;
137 TString GetHLTTriggerDescription() const;
138 Bool_t IsHLTTriggerFired(const char* name=NULL) const;
139 UInt_t GetOrbitNumber() const {return fHeader->GetOrbitNumber();}
140 UInt_t GetTimeStamp() const { return fHeader->GetTimeStamp();}
141 UInt_t GetEventType() const { return fHeader->GetEventType();}
142 Int_t GetEventNumberInFile() const {return fHeader->GetEventNumberInFile();}
143 UShort_t GetBunchCrossNumber() const {return fHeader->GetBunchCrossNumber();}
144 UChar_t GetTriggerCluster() const {return fHeader->GetTriggerCluster();}
146 // ZDC CKB: put this in the header?
147 AliESDZDC* GetESDZDC() const {return fESDZDC;}
149 // Delegated methods for fESDZDC
150 Double_t GetZDCN1Energy() const {return fESDZDC->GetZDCN1Energy();}
151 Double_t GetZDCP1Energy() const {return fESDZDC->GetZDCP1Energy();}
152 Double_t GetZDCN2Energy() const {return fESDZDC->GetZDCN2Energy();}
153 Double_t GetZDCP2Energy() const {return fESDZDC->GetZDCP2Energy();}
154 Double_t GetZDCEMEnergy(Int_t i=0) const {return fESDZDC->GetZDCEMEnergy(i);}
155 Int_t GetZDCParticipants() const {return fESDZDC->GetZDCParticipants();}
156 void SetZDC(Float_t n1Energy, Float_t p1Energy, Float_t em1Energy, Float_t em2Energy,
157 Float_t n2Energy, Float_t p2Energy, Int_t participants, Int_t nPartA,
158 Int_t nPartC, Double_t b, Double_t bA, Double_t bC, UInt_t recoflag)
159 {fESDZDC->SetZDC(n1Energy, p1Energy, em1Energy, em2Energy, n2Energy, p2Energy,
160 participants, nPartA, nPartC, b, bA, bC, recoflag);}
164 void SetFMDData(AliESDFMD * obj);
165 AliESDFMD *GetFMDData(){ return fESDFMD; }
168 // TZERO CKB: put this in the header?
169 const AliESDTZERO* GetESDTZERO() const {return fESDTZERO;}
170 // delegetated methods for fESDTZERO
172 Double_t GetT0zVertex() const {return fESDTZERO->GetT0zVertex();}
173 void SetT0zVertex(Float_t z) {fESDTZERO->SetT0zVertex(z);}
174 Double_t GetT0() const {return fESDTZERO->GetT0();}
175 void SetT0(Float_t timeStart) {fESDTZERO->SetT0(timeStart);}
176 Float_t GetT0clock() const {return fESDTZERO->GetT0clock();}
177 void SetT0clock(Float_t timeStart) {fESDTZERO->SetT0clock(timeStart);}
178 const Double_t * GetT0time() const {return fESDTZERO->GetT0time();}
179 void SetT0time(Float_t time[24]) {fESDTZERO->SetT0time(time);}
180 const Double_t * GetT0amplitude() const {return fESDTZERO->GetT0amplitude();}
181 void SetT0amplitude(Float_t amp[24]){fESDTZERO->SetT0amplitude(amp);}
184 AliESDVZERO *GetVZEROData() const { return fESDVZERO; }
185 void SetVZEROData(AliESDVZERO * obj);
188 AliESDACORDE *GetACORDEData() const { return fESDACORDE;}
189 void SetACORDEData(AliESDACORDE * obj);
191 void SetESDfriend(const AliESDfriend *f) const;
192 void GetESDfriend(AliESDfriend *f) const;
196 void SetPrimaryVertexTPC(const AliESDVertex *vertex);
197 const AliESDVertex *GetPrimaryVertexTPC() const {return fTPCVertex;}
199 void SetPrimaryVertexSPD(const AliESDVertex *vertex);
200 const AliESDVertex *GetPrimaryVertexSPD() const {return fSPDVertex;}
201 const AliESDVertex *GetVertex() const {
202 //For the backward compatibily only
203 return GetPrimaryVertexSPD();
206 void SetPrimaryVertexTracks(const AliESDVertex *vertex);
207 const AliESDVertex *GetPrimaryVertexTracks() const {return fPrimaryVertex;}
208 AliESDVertex *PrimaryVertexTracksUnconstrained() const;
210 const AliESDVertex *GetPrimaryVertex() const;
212 void SetMultiplicity(const AliMultiplicity *mul);
214 const AliMultiplicity *GetMultiplicity() const {return fSPDMult;}
217 Bool_t Clean(Float_t *cleanPars);
218 Bool_t RemoveKink(Int_t i) const;
219 Bool_t RemoveV0(Int_t i) const;
220 Bool_t RemoveTrack(Int_t i) const;
222 const AliESDVertex *GetPileupVertexSPD(Int_t i) const {
223 return (const AliESDVertex *)fSPDPileupVertices->UncheckedAt(i);
225 Char_t AddPileupVertexSPD(const AliESDVertex *vtx);
227 const AliESDVertex *GetPileupVertexTracks(Int_t i) const {
228 return (const AliESDVertex *)fTrkPileupVertices->UncheckedAt(i);
230 Char_t AddPileupVertexTracks(const AliESDVertex *vtx);
232 AliESDtrack *GetTrack(Int_t i) const {
233 return (AliESDtrack *)fTracks->UncheckedAt(i);
235 Int_t AddTrack(const AliESDtrack *t);
238 AliESDHLTtrack *GetHLTConfMapTrack(Int_t /*i*/) const {
239 // return (AliESDHLTtrack *)fHLTConfMapTracks->UncheckedAt(i);
242 void AddHLTConfMapTrack(const AliESDHLTtrack */*t*/) {
243 printf("ESD:: AddHLTConfMapTrack do nothing \n");
244 // TClonesArray &fhlt = *fHLTConfMapTracks;
245 // new(fhlt[fHLTConfMapTracks->GetEntriesFast()]) AliESDHLTtrack(*t);
249 AliESDHLTtrack *GetHLTHoughTrack(Int_t /*i*/) const {
250 // return (AliESDHLTtrack *)fHLTHoughTracks->UncheckedAt(i);
253 void AddHLTHoughTrack(const AliESDHLTtrack */*t*/) {
254 printf("ESD:: AddHLTHoughTrack do nothing \n");
255 // TClonesArray &fhlt = *fHLTHoughTracks;
256 // new(fhlt[fHLTHoughTracks->GetEntriesFast()]) AliESDHLTtrack(*t);
259 AliESDMuonTrack *GetMuonTrack(Int_t i) const {
260 return (AliESDMuonTrack *)fMuonTracks->UncheckedAt(i);
263 void AddMuonTrack(const AliESDMuonTrack *t);
265 AliESDPmdTrack *GetPmdTrack(Int_t i) const {
266 return (AliESDPmdTrack *)fPmdTracks->UncheckedAt(i);
269 void AddPmdTrack(const AliESDPmdTrack *t);
272 AliESDTrdTrack *GetTrdTrack(Int_t i) const {
273 return (AliESDTrdTrack *)fTrdTracks->UncheckedAt(i);
277 void AddTrdTrack(const AliESDTrdTrack *t);
279 AliESDv0 *GetV0(Int_t i) const {
280 return (AliESDv0*)fV0s->UncheckedAt(i);
282 Int_t AddV0(const AliESDv0 *v);
284 AliESDcascade *GetCascade(Int_t i) const {
285 return (AliESDcascade *)fCascades->UncheckedAt(i);
288 void AddCascade(const AliESDcascade *c);
290 AliESDkink *GetKink(Int_t i) const {
291 return (AliESDkink *)fKinks->UncheckedAt(i);
293 Int_t AddKink(const AliESDkink *c);
295 AliESDCaloCluster *GetCaloCluster(Int_t i) const {
296 return (AliESDCaloCluster *)fCaloClusters->UncheckedAt(i);
299 Int_t AddCaloCluster(const AliESDCaloCluster *c);
301 AliESDCaloCells *GetEMCALCells() const {return fEMCALCells; }
302 AliESDCaloCells *GetPHOSCells() const {return fPHOSCells; }
304 AliRawDataErrorLog *GetErrorLog(Int_t i) const {
305 return (AliRawDataErrorLog *)fErrorLogs->UncheckedAt(i);
307 void AddRawDataErrorLog(const AliRawDataErrorLog *log) const;
309 Int_t GetNumberOfErrorLogs() const {return fErrorLogs->GetEntriesFast();}
312 void AddPHOSTriggerPosition(TArrayF array) { fPHOSTrigger->AddTriggerPosition(array); }
313 void AddPHOSTriggerAmplitudes(TArrayF array) { fPHOSTrigger->AddTriggerAmplitudes(array);}
314 void AddEMCALTriggerPosition(TArrayF array) { fEMCALTrigger->AddTriggerPosition(array); }
315 void AddEMCALTriggerAmplitudes(TArrayF array){ fEMCALTrigger->AddTriggerAmplitudes(array); }
317 Int_t GetNumberOfPileupVerticesSPD() const {
318 return (fSPDPileupVertices?fSPDPileupVertices->GetEntriesFast():0);
320 Int_t GetNumberOfPileupVerticesTracks() const {
321 return (fTrkPileupVertices?fTrkPileupVertices->GetEntriesFast():0);
323 Int_t GetNumberOfTracks() const {return fTracks->GetEntriesFast();}
324 Int_t GetNumberOfHLTConfMapTracks() const {return 0;}
325 // fHLTConfMapTracks->GetEntriesFast();}
326 Int_t GetNumberOfHLTHoughTracks() const {return 0; }
327 // fHLTHoughTracks->GetEntriesFast(); }
329 Int_t GetNumberOfMuonTracks() const {return fMuonTracks->GetEntriesFast();}
330 Int_t GetNumberOfPmdTracks() const {return fPmdTracks->GetEntriesFast();}
331 Int_t GetNumberOfTrdTracks() const {return fTrdTracks->GetEntriesFast();}
332 Int_t GetNumberOfV0s() const {return fV0s->GetEntriesFast();}
333 Int_t GetNumberOfCascades() const {return fCascades->GetEntriesFast();}
334 Int_t GetNumberOfKinks() const {return fKinks->GetEntriesFast();}
336 Int_t GetEMCALClusters(TRefArray *clusters) const;
337 Int_t GetPHOSClusters(TRefArray *clusters) const;
338 Int_t GetNumberOfCaloClusters() const {return fCaloClusters->GetEntriesFast();}
340 void SetUseOwnList(Bool_t b){fUseOwnList = b;}
341 Bool_t GetUseOwnList(){return fUseOwnList;}
343 // Remove this stuff CKB?
344 //---------------------------------------------------
345 Int_t GetNumberOfEMCALClusters() const {return fEMCALClusters;}
346 void SetNumberOfEMCALClusters(Int_t clus) {fEMCALClusters = clus;}
347 Int_t GetFirstEMCALCluster() const {return fFirstEMCALCluster;}
348 void SetFirstEMCALCluster(Int_t index) {fFirstEMCALCluster = index;}
350 Int_t GetNumberOfPHOSClusters() const {return fPHOSClusters;}
351 void SetNumberOfPHOSClusters(Int_t part) { fPHOSClusters = part ; }
352 void SetFirstPHOSCluster(Int_t index) { fFirstPHOSCluster = index ; }
353 Int_t GetFirstPHOSCluster() const { return fFirstPHOSCluster ; }
354 //-------------------------------------------------------
356 TArrayF *GetEMCALTriggerPosition() const {return fEMCALTrigger->GetTriggerPosition();}
357 TArrayF *GetEMCALTriggerAmplitudes() const {return fEMCALTrigger->GetTriggerAmplitudes();}
358 TArrayF *GetPHOSTriggerPosition() const {return fPHOSTrigger->GetTriggerPosition();}
359 TArrayF *GetPHOSTriggerAmplitudes() const {return fPHOSTrigger->GetTriggerAmplitudes();}
361 void ResetV0s() { fV0s->Clear(); }
362 void ResetCascades() { fCascades->Clear(); }
365 void Print(Option_t *option="") const;
367 void AddObject(TObject* obj);
368 void ReadFromTree(TTree *tree, Option_t* opt = "");
369 TObject* FindListObject(const char *name);
370 AliESD *GetAliESDOld(){return fESDOld;}
371 void WriteToTree(TTree* tree) const;
372 void GetStdContent();
373 void ResetStdContent();
374 void CreateStdContent();
375 void CreateStdContent(Bool_t bUseThisList);
377 void CopyFromOldESD();
378 TList* GetList() const {return fESDObjects;}
381 AliESDEvent(const AliESDEvent&);
382 static Bool_t ResetWithPlacementNew(TObject *pObject);
384 TList *fESDObjects; // List of esd Objects
386 AliESDRun *fESDRun; //! Run information tmp put in the Userdata
387 AliESDHeader *fHeader; //! ESD Event Header
388 AliESDZDC *fESDZDC; //! ZDC information
389 AliESDFMD *fESDFMD; //! FMD object containing rough multiplicity
390 AliESDVZERO *fESDVZERO; //! VZERO object containing rough multiplicity
391 AliESDTZERO *fESDTZERO; //! TZEROObject
392 AliESDVertex *fTPCVertex; //! Primary vertex estimated by the TPC
393 AliESDVertex *fSPDVertex; //! Primary vertex estimated by the SPD
394 AliESDVertex *fPrimaryVertex; //! Primary vertex estimated using ESD tracks
395 AliMultiplicity *fSPDMult; //! SPD tracklet multiplicity
396 AliESDCaloTrigger* fPHOSTrigger; //! PHOS Trigger information
397 AliESDCaloTrigger* fEMCALTrigger; //! PHOS Trigger information
398 AliESDACORDE *fESDACORDE; //! ACORDE ESD object caontaining bit pattern
400 TClonesArray *fSPDPileupVertices;//! Pileup primary vertices reconstructed by SPD
401 TClonesArray *fTrkPileupVertices;//! Pileup primary vertices reconstructed using the tracks
402 TClonesArray *fTracks; //! ESD tracks
403 TClonesArray *fMuonTracks; //! MUON ESD tracks
404 TClonesArray *fPmdTracks; //! PMD ESD tracks
405 TClonesArray *fTrdTracks; //! TRD ESD tracks (triggered)
406 TClonesArray *fV0s; //! V0 vertices
407 TClonesArray *fCascades; //! Cascade vertices
408 TClonesArray *fKinks; //! Kinks
409 TClonesArray *fCaloClusters; //! Calorimeter clusters for PHOS/EMCAL
410 AliESDCaloCells *fEMCALCells; //! EMCAL cell info
411 AliESDCaloCells *fPHOSCells; //! PHOS cell info
412 TClonesArray *fErrorLogs; //! Raw-data reading error messages
416 AliESD *fESDOld; //! Old esd Structure
417 AliESDfriend *fESDFriendOld; //! Old friend esd Structure
418 Bool_t fConnected; //! flag if leaves are alreday connected
419 Bool_t fUseOwnList; //! Do not use the list from the esdTree but use the one created by this class
421 static const char* fgkESDListName[kESDListN]; //!
423 // Remove this stuff CKB
424 Int_t fEMCALClusters; // Number of EMCAL clusters (subset of caloclusters)
425 Int_t fFirstEMCALCluster; // First EMCAL cluster in the fCaloClusters list
427 Int_t fPHOSClusters; // Number of PHOS clusters (subset of caloclusters)
428 Int_t fFirstPHOSCluster; // First PHOS cluster in the fCaloClusters list
430 ClassDef(AliESDEvent,10) //ESDEvent class