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);}
118 void SetEMCALMatrix(TGeoHMatrix*matrix, Int_t i) {fESDRun->SetEMCALMatrix(matrix,i);}
119 const TGeoHMatrix* GetEMCALMatrix(Int_t i) const {return fESDRun->GetEMCALMatrix(i);}
121 void SetCurrentL3(Float_t cur) const {fESDRun->SetCurrentL3(cur);}
122 void SetCurrentDip(Float_t cur) const {fESDRun->SetCurrentDip(cur);}
123 void SetBeamEnergy(Float_t be) const {fESDRun->SetBeamEnergy(be);}
124 void SetBeamType(const char* bt) const {fESDRun->SetBeamType(bt);}
125 void SetUniformBMap(Bool_t val=kTRUE) const {fESDRun->SetBit(AliESDRun::kUniformBMap,val);}
126 void SetBInfoStored(Bool_t val=kTRUE) const {fESDRun->SetBit(AliESDRun::kBInfoStored,val);}
128 Float_t GetCurrentL3() const {return fESDRun->GetCurrentL3();}
129 Float_t GetCurrentDip() const {return fESDRun->GetCurrentDip();}
130 Float_t GetBeamEnergy() const {return fESDRun->GetBeamEnergy();}
131 const char* GetBeamType() const {return fESDRun->GetBeamType();}
132 Bool_t IsUniformBMap() const {return fESDRun->TestBit(AliESDRun::kUniformBMap);}
134 Bool_t InitMagneticField() const {return fESDRun->InitMagneticField();}
136 AliESDHeader* GetHeader() const {return fHeader;}
138 // Delegated methods for fHeader
139 void SetTriggerMask(ULong64_t n) {fHeader->SetTriggerMask(n);}
140 void SetOrbitNumber(UInt_t n) {fHeader->SetOrbitNumber(n);}
141 void SetTimeStamp(UInt_t timeStamp){fHeader->SetTimeStamp(timeStamp);}
142 void SetEventType(UInt_t eventType){fHeader->SetEventType(eventType);}
143 void SetEventNumberInFile(Int_t n) {fHeader->SetEventNumberInFile(n);}
144 // void SetRunNumber(Int_t n) {fHeader->SetRunNumber(n);}
145 void SetBunchCrossNumber(UShort_t n) {fHeader->SetBunchCrossNumber(n);}
146 void SetTriggerCluster(UChar_t n) {fHeader->SetTriggerCluster(n);}
148 ULong64_t GetTriggerMask() const {return fHeader->GetTriggerMask();}
149 TString GetFiredTriggerClasses() const {return fESDRun->GetFiredTriggerClasses(fHeader->GetTriggerMask());}
150 Bool_t IsTriggerClassFired(const char *name) const {return fESDRun->IsTriggerClassFired(fHeader->GetTriggerMask(),name);}
151 Bool_t IsEventSelected(const char *trigExpr) const;
152 TObject* GetHLTTriggerDecision() const;
153 TString GetHLTTriggerDescription() const;
154 Bool_t IsHLTTriggerFired(const char* name=NULL) const;
155 UInt_t GetOrbitNumber() const {return fHeader->GetOrbitNumber();}
156 UInt_t GetTimeStamp() const { return fHeader->GetTimeStamp();}
157 UInt_t GetEventType() const { return fHeader->GetEventType();}
158 Int_t GetEventNumberInFile() const {return fHeader->GetEventNumberInFile();}
159 UShort_t GetBunchCrossNumber() const {return fHeader->GetBunchCrossNumber();}
160 UChar_t GetTriggerCluster() const {return fHeader->GetTriggerCluster();}
162 // ZDC CKB: put this in the header?
163 AliESDZDC* GetESDZDC() const {return fESDZDC;}
165 // Delegated methods for fESDZDC
166 Double_t GetZDCN1Energy() const {return fESDZDC->GetZDCN1Energy();}
167 Double_t GetZDCP1Energy() const {return fESDZDC->GetZDCP1Energy();}
168 Double_t GetZDCN2Energy() const {return fESDZDC->GetZDCN2Energy();}
169 Double_t GetZDCP2Energy() const {return fESDZDC->GetZDCP2Energy();}
170 Double_t GetZDCEMEnergy(Int_t i=0) const {return fESDZDC->GetZDCEMEnergy(i);}
171 Int_t GetZDCParticipants() const {return fESDZDC->GetZDCParticipants();}
172 void SetZDC(Float_t n1Energy, Float_t p1Energy, Float_t em1Energy, Float_t em2Energy,
173 Float_t n2Energy, Float_t p2Energy, Int_t participants, Int_t nPartA,
174 Int_t nPartC, Double_t b, Double_t bA, Double_t bC, UInt_t recoflag)
175 {fESDZDC->SetZDC(n1Energy, p1Energy, em1Energy, em2Energy, n2Energy, p2Energy,
176 participants, nPartA, nPartC, b, bA, bC, recoflag);}
180 void SetFMDData(AliESDFMD * obj);
181 AliESDFMD *GetFMDData(){ return fESDFMD; }
184 // TZERO CKB: put this in the header?
185 const AliESDTZERO* GetESDTZERO() const {return fESDTZERO;}
186 // delegetated methods for fESDTZERO
188 Double_t GetT0zVertex() const {return fESDTZERO->GetT0zVertex();}
189 void SetT0zVertex(Float_t z) {fESDTZERO->SetT0zVertex(z);}
190 Double_t GetT0() const {return fESDTZERO->GetT0();}
191 void SetT0(Float_t timeStart) {fESDTZERO->SetT0(timeStart);}
192 Float_t GetT0clock() const {return fESDTZERO->GetT0clock();}
193 void SetT0clock(Float_t timeStart) {fESDTZERO->SetT0clock(timeStart);}
194 const Double_t * GetT0time() const {return fESDTZERO->GetT0time();}
195 void SetT0time(Float_t time[24]) {fESDTZERO->SetT0time(time);}
196 const Double_t * GetT0amplitude() const {return fESDTZERO->GetT0amplitude();}
197 void SetT0amplitude(Float_t amp[24]){fESDTZERO->SetT0amplitude(amp);}
200 AliESDVZERO *GetVZEROData() const { return fESDVZERO; }
201 void SetVZEROData(AliESDVZERO * obj);
204 AliESDACORDE *GetACORDEData() const { return fESDACORDE;}
205 void SetACORDEData(AliESDACORDE * obj);
207 void SetESDfriend(const AliESDfriend *f) const;
208 void GetESDfriend(AliESDfriend *f) const;
212 void SetPrimaryVertexTPC(const AliESDVertex *vertex);
213 const AliESDVertex *GetPrimaryVertexTPC() const {return fTPCVertex;}
215 void SetPrimaryVertexSPD(const AliESDVertex *vertex);
216 const AliESDVertex *GetPrimaryVertexSPD() const {return fSPDVertex;}
217 const AliESDVertex *GetVertex() const {
218 //For the backward compatibily only
219 return GetPrimaryVertexSPD();
222 void SetPrimaryVertexTracks(const AliESDVertex *vertex);
223 const AliESDVertex *GetPrimaryVertexTracks() const {return fPrimaryVertex;}
224 AliESDVertex *PrimaryVertexTracksUnconstrained() const;
226 const AliESDVertex *GetPrimaryVertex() const;
228 void SetMultiplicity(const AliMultiplicity *mul);
230 const AliMultiplicity *GetMultiplicity() const {return fSPDMult;}
233 Bool_t Clean(Float_t *cleanPars);
234 Bool_t RemoveKink(Int_t i) const;
235 Bool_t RemoveV0(Int_t i) const;
236 Bool_t RemoveTrack(Int_t i) const;
238 const AliESDVertex *GetPileupVertexSPD(Int_t i) const {
239 return (const AliESDVertex *)fSPDPileupVertices->UncheckedAt(i);
241 Char_t AddPileupVertexSPD(const AliESDVertex *vtx);
242 Bool_t IsPileupFromSPD(Int_t ncont=2, Double_t nSigmaDeltaZ=3., Double_t nSigmaXY=2., Int_t option=0) const;
244 const AliESDVertex *GetPileupVertexTracks(Int_t i) const {
245 return (const AliESDVertex *)fTrkPileupVertices->UncheckedAt(i);
247 Char_t AddPileupVertexTracks(const AliESDVertex *vtx);
249 AliESDtrack *GetTrack(Int_t i) const {
250 return (AliESDtrack *)fTracks->UncheckedAt(i);
252 Int_t AddTrack(const AliESDtrack *t);
255 AliESDHLTtrack *GetHLTConfMapTrack(Int_t /*i*/) const {
256 // return (AliESDHLTtrack *)fHLTConfMapTracks->UncheckedAt(i);
259 void AddHLTConfMapTrack(const AliESDHLTtrack */*t*/) {
260 printf("ESD:: AddHLTConfMapTrack do nothing \n");
261 // TClonesArray &fhlt = *fHLTConfMapTracks;
262 // new(fhlt[fHLTConfMapTracks->GetEntriesFast()]) AliESDHLTtrack(*t);
266 AliESDHLTtrack *GetHLTHoughTrack(Int_t /*i*/) const {
267 // return (AliESDHLTtrack *)fHLTHoughTracks->UncheckedAt(i);
270 void AddHLTHoughTrack(const AliESDHLTtrack */*t*/) {
271 printf("ESD:: AddHLTHoughTrack do nothing \n");
272 // TClonesArray &fhlt = *fHLTHoughTracks;
273 // new(fhlt[fHLTHoughTracks->GetEntriesFast()]) AliESDHLTtrack(*t);
276 AliESDMuonTrack *GetMuonTrack(Int_t i) const {
277 return (AliESDMuonTrack *)fMuonTracks->UncheckedAt(i);
280 void AddMuonTrack(const AliESDMuonTrack *t);
282 AliESDPmdTrack *GetPmdTrack(Int_t i) const {
283 return (AliESDPmdTrack *)fPmdTracks->UncheckedAt(i);
286 void AddPmdTrack(const AliESDPmdTrack *t);
289 AliESDTrdTrack *GetTrdTrack(Int_t i) const {
290 return (AliESDTrdTrack *)fTrdTracks->UncheckedAt(i);
294 void AddTrdTrack(const AliESDTrdTrack *t);
296 AliESDv0 *GetV0(Int_t i) const {
297 return (AliESDv0*)fV0s->UncheckedAt(i);
299 Int_t AddV0(const AliESDv0 *v);
301 AliESDcascade *GetCascade(Int_t i) const {
302 return (AliESDcascade *)fCascades->UncheckedAt(i);
305 void AddCascade(const AliESDcascade *c);
307 AliESDkink *GetKink(Int_t i) const {
308 return (AliESDkink *)fKinks->UncheckedAt(i);
310 Int_t AddKink(const AliESDkink *c);
312 AliESDCaloCluster *GetCaloCluster(Int_t i) const {
313 return (AliESDCaloCluster *)fCaloClusters->UncheckedAt(i);
316 Int_t AddCaloCluster(const AliESDCaloCluster *c);
318 AliESDCaloCells *GetEMCALCells() const {return fEMCALCells; }
319 AliESDCaloCells *GetPHOSCells() const {return fPHOSCells; }
321 AliRawDataErrorLog *GetErrorLog(Int_t i) const {
322 return (AliRawDataErrorLog *)fErrorLogs->UncheckedAt(i);
324 void AddRawDataErrorLog(const AliRawDataErrorLog *log) const;
326 Int_t GetNumberOfErrorLogs() const {return fErrorLogs->GetEntriesFast();}
329 void AddPHOSTriggerPosition(TArrayF array) { fPHOSTrigger->AddTriggerPosition(array); }
330 void AddPHOSTriggerAmplitudes(TArrayF array) { fPHOSTrigger->AddTriggerAmplitudes(array);}
331 void AddEMCALTriggerPosition(TArrayF array) { fEMCALTrigger->AddTriggerPosition(array); }
332 void AddEMCALTriggerAmplitudes(TArrayF array){ fEMCALTrigger->AddTriggerAmplitudes(array); }
334 Int_t GetNumberOfPileupVerticesSPD() const {
335 return (fSPDPileupVertices?fSPDPileupVertices->GetEntriesFast():0);
337 Int_t GetNumberOfPileupVerticesTracks() const {
338 return (fTrkPileupVertices?fTrkPileupVertices->GetEntriesFast():0);
340 Int_t GetNumberOfTracks() const {return fTracks->GetEntriesFast();}
341 Int_t GetNumberOfHLTConfMapTracks() const {return 0;}
342 // fHLTConfMapTracks->GetEntriesFast();}
343 Int_t GetNumberOfHLTHoughTracks() const {return 0; }
344 // fHLTHoughTracks->GetEntriesFast(); }
346 Int_t GetNumberOfMuonTracks() const {return fMuonTracks->GetEntriesFast();}
347 Int_t GetNumberOfPmdTracks() const {return fPmdTracks->GetEntriesFast();}
348 Int_t GetNumberOfTrdTracks() const {return fTrdTracks->GetEntriesFast();}
349 Int_t GetNumberOfV0s() const {return fV0s->GetEntriesFast();}
350 Int_t GetNumberOfCascades() const {return fCascades->GetEntriesFast();}
351 Int_t GetNumberOfKinks() const {return fKinks->GetEntriesFast();}
353 Int_t GetEMCALClusters(TRefArray *clusters) const;
354 Int_t GetPHOSClusters(TRefArray *clusters) const;
355 Int_t GetNumberOfCaloClusters() const {return fCaloClusters->GetEntriesFast();}
357 void SetUseOwnList(Bool_t b){fUseOwnList = b;}
358 Bool_t GetUseOwnList(){return fUseOwnList;}
360 // Remove this stuff CKB?
361 //---------------------------------------------------
362 Int_t GetNumberOfEMCALClusters() const {return fEMCALClusters;}
363 void SetNumberOfEMCALClusters(Int_t clus) {fEMCALClusters = clus;}
364 Int_t GetFirstEMCALCluster() const {return fFirstEMCALCluster;}
365 void SetFirstEMCALCluster(Int_t index) {fFirstEMCALCluster = index;}
367 Int_t GetNumberOfPHOSClusters() const {return fPHOSClusters;}
368 void SetNumberOfPHOSClusters(Int_t part) { fPHOSClusters = part ; }
369 void SetFirstPHOSCluster(Int_t index) { fFirstPHOSCluster = index ; }
370 Int_t GetFirstPHOSCluster() const { return fFirstPHOSCluster ; }
371 //-------------------------------------------------------
373 TArrayF *GetEMCALTriggerPosition() const {return fEMCALTrigger->GetTriggerPosition();}
374 TArrayF *GetEMCALTriggerAmplitudes() const {return fEMCALTrigger->GetTriggerAmplitudes();}
375 TArrayF *GetPHOSTriggerPosition() const {return fPHOSTrigger->GetTriggerPosition();}
376 TArrayF *GetPHOSTriggerAmplitudes() const {return fPHOSTrigger->GetTriggerAmplitudes();}
378 void ResetV0s() { fV0s->Clear(); }
379 void ResetCascades() { fCascades->Clear(); }
382 void Print(Option_t *option="") const;
384 void AddObject(TObject* obj);
385 void ReadFromTree(TTree *tree, Option_t* opt = "");
386 TObject* FindListObject(const char *name);
387 AliESD *GetAliESDOld(){return fESDOld;}
388 void WriteToTree(TTree* tree) const;
389 void GetStdContent();
390 void ResetStdContent();
391 void CreateStdContent();
392 void CreateStdContent(Bool_t bUseThisList);
394 void CopyFromOldESD();
395 TList* GetList() const {return fESDObjects;}
398 AliESDEvent(const AliESDEvent&);
399 static Bool_t ResetWithPlacementNew(TObject *pObject);
401 TList *fESDObjects; // List of esd Objects
403 AliESDRun *fESDRun; //! Run information tmp put in the Userdata
404 AliESDHeader *fHeader; //! ESD Event Header
405 AliESDZDC *fESDZDC; //! ZDC information
406 AliESDFMD *fESDFMD; //! FMD object containing rough multiplicity
407 AliESDVZERO *fESDVZERO; //! VZERO object containing rough multiplicity
408 AliESDTZERO *fESDTZERO; //! TZEROObject
409 AliESDVertex *fTPCVertex; //! Primary vertex estimated by the TPC
410 AliESDVertex *fSPDVertex; //! Primary vertex estimated by the SPD
411 AliESDVertex *fPrimaryVertex; //! Primary vertex estimated using ESD tracks
412 AliMultiplicity *fSPDMult; //! SPD tracklet multiplicity
413 AliESDCaloTrigger* fPHOSTrigger; //! PHOS Trigger information
414 AliESDCaloTrigger* fEMCALTrigger; //! PHOS Trigger information
415 AliESDACORDE *fESDACORDE; //! ACORDE ESD object caontaining bit pattern
417 TClonesArray *fSPDPileupVertices;//! Pileup primary vertices reconstructed by SPD
418 TClonesArray *fTrkPileupVertices;//! Pileup primary vertices reconstructed using the tracks
419 TClonesArray *fTracks; //! ESD tracks
420 TClonesArray *fMuonTracks; //! MUON ESD tracks
421 TClonesArray *fPmdTracks; //! PMD ESD tracks
422 TClonesArray *fTrdTracks; //! TRD ESD tracks (triggered)
423 TClonesArray *fV0s; //! V0 vertices
424 TClonesArray *fCascades; //! Cascade vertices
425 TClonesArray *fKinks; //! Kinks
426 TClonesArray *fCaloClusters; //! Calorimeter clusters for PHOS/EMCAL
427 AliESDCaloCells *fEMCALCells; //! EMCAL cell info
428 AliESDCaloCells *fPHOSCells; //! PHOS cell info
429 TClonesArray *fErrorLogs; //! Raw-data reading error messages
433 AliESD *fESDOld; //! Old esd Structure
434 AliESDfriend *fESDFriendOld; //! Old friend esd Structure
435 Bool_t fConnected; //! flag if leaves are alreday connected
436 Bool_t fUseOwnList; //! Do not use the list from the esdTree but use the one created by this class
438 static const char* fgkESDListName[kESDListN]; //!
440 // Remove this stuff CKB
441 Int_t fEMCALClusters; // Number of EMCAL clusters (subset of caloclusters)
442 Int_t fFirstEMCALCluster; // First EMCAL cluster in the fCaloClusters list
444 Int_t fPHOSClusters; // Number of PHOS clusters (subset of caloclusters)
445 Int_t fFirstPHOSCluster; // First PHOS cluster in the fCaloClusters list
447 ClassDef(AliESDEvent,10) //ESDEvent class