]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - STEER/AliESDEvent.h
Herwig event header added
[u/mrichter/AliRoot.git] / STEER / AliESDEvent.h
index fe53fcf23f46b7d7386a91a7bbca56e394447470..b1272d48d3b2a0b5e9cf69f1d2e3d1d68ae769a5 100644 (file)
@@ -9,7 +9,7 @@
 
 //-------------------------------------------------------------------------
 //                          Class AliESDEvent
-//   This is the class to deal with during the physical analysis of data.
+//   This is the class to deal with during the physics analysis of data.
 //   It also ensures the backward compatibility with the old ESD format.
 //      
 // Origin: Christian Klein-Boesing, CERN, Christian.Klein-Boesing@cern.ch 
 
 class TList;
 
-
-// some includes for delegeated methds
+#include "AliVEvent.h"
+// some includes for delegated methods
 #include "AliESDCaloTrigger.h"
 #include "AliESDRun.h"
 #include "AliESDHeader.h"
 #include "AliESDTZERO.h"
 #include "AliESDZDC.h"
+#include "AliESDACORDE.h"
+
+// AliESDtrack has to be included so that the compiler 
+// knows its inheritance tree (= that it is a AliVParticle).
+#include "AliESDtrack.h"
+// same for AliESDVertex (which is a AliVVertex)
+#include "AliESDVertex.h"
 
 class AliESDfriend;
 class AliESDVZERO;
@@ -37,10 +44,9 @@ class AliESDVertex;
 class AliESDPmdTrack;
 class AliESDFMD;
 class AliESDkink;
-class AliESDtrack;
 class AliESDCaloCluster;
+class AliESDCaloCells;
 class AliESDv0;
-class AliESDFMD;
 class AliMultiplicity;
 class AliRawDataErrorLog;
 class AliESDRun;
@@ -48,22 +54,28 @@ class AliESDTrdTrack;
 class AliESDMuonTrack;
 class AliESD;
 class AliESDcascade;
+class TRefArray;
+class AliESDACORDE;
+class AliESDHLTDecision;
 
-class AliESDEvent : public TObject {
+class AliESDEvent : public AliVEvent {
 public:
 
 
-  enum ESDListIndex_t   {kESDRun,
+  enum ESDListIndex   {kESDRun,
                       kHeader,
                       kESDZDC,
                       kESDFMD,
                       kESDVZERO,
                       kESDTZERO,
+                      kTPCVertex,
                       kSPDVertex,
                       kPrimaryVertex,
                       kSPDMult,
                       kPHOSTrigger,
                       kEMCALTrigger,
+                      kSPDPileupVertices,
+                      kTrkPileupVertices,
                       kTracks,
                       kMuonTracks,
                       kPmdTracks,
@@ -72,108 +84,172 @@ public:
                       kCascades,
                       kKinks,
                       kCaloClusters,
+                      kEMCALCells,
+                      kPHOSCells,
                       kErrorLogs,
+                       kESDACORDE,
                       kESDListN
   };
 
   AliESDEvent();
-  virtual ~AliESDEvent(); 
-
+  virtual ~AliESDEvent();
+  AliESDEvent &operator=(const AliESDEvent& source); // or make private and use only copy? 
+  virtual void Copy(TObject& obj) const;
 
   // RUN
   // move this to the UserData!!!
   const AliESDRun*    GetESDRun() const {return fESDRun;}
 
   // Delegated methods for fESDRun
-  void    SetRunNumber(Int_t n) {fESDRun->SetRunNumber(n);}
-  Int_t   GetRunNumber() const {return fESDRun->GetRunNumber();}
-  void    SetPeriodNumber(Int_t n){fESDRun->SetPeriodNumber(n);}
-  Int_t   GetPeriodNumber() const {return fESDRun->GetPeriodNumber();}
-  void    SetMagneticField(Float_t mf){fESDRun->SetMagneticField(mf);}
-  Float_t GetMagneticField() const {return fESDRun->GetMagneticField();}
-  void SetDiamond(const AliESDVertex *vertex) { fESDRun->SetDiamond(vertex);}
-  Float_t GetDiamondX() const {return fESDRun->GetDiamondX();}
-  Float_t GetDiamondY() const {return fESDRun->GetDiamondY();}
-  Float_t GetSigma2DiamondX() const {return  fESDRun->GetSigma2DiamondX();}
-  Float_t GetSigma2DiamondY() const {return  fESDRun->GetSigma2DiamondY();}
-  void GetDiamondCovXY(Float_t cov[3]) const {fESDRun->GetDiamondCovXY(cov);}   
-  
-
+  void     SetRunNumber(Int_t n) {fESDRun->SetRunNumber(n);}
+  Int_t    GetRunNumber() const {return fESDRun->GetRunNumber();}
+  void     SetPeriodNumber(UInt_t n){fESDRun->SetPeriodNumber(n);}
+  UInt_t   GetPeriodNumber() const {return fESDRun->GetPeriodNumber();}
+  void     SetMagneticField(Double_t mf){fESDRun->SetMagneticField(mf);}
+  Double_t GetMagneticField() const {return fESDRun->GetMagneticField();}
+  void     SetDiamond(const AliESDVertex *vertex) { fESDRun->SetDiamond(vertex);}
+  Double_t  GetDiamondX() const {return fESDRun->GetDiamondX();}
+  Double_t  GetDiamondY() const {return fESDRun->GetDiamondY();}
+  Double_t  GetSigma2DiamondX() const {return  fESDRun->GetSigma2DiamondX();}
+  Double_t  GetSigma2DiamondY() const {return  fESDRun->GetSigma2DiamondY();}
+  void      GetDiamondCovXY(Float_t cov[3]) const {fESDRun->GetDiamondCovXY(cov);}   
+  void     SetTriggerClass(const char*name, Int_t index) {fESDRun->SetTriggerClass(name,index);}
+  void     SetPHOSMatrix(TGeoHMatrix*matrix, Int_t i) {fESDRun->SetPHOSMatrix(matrix,i);}
+  const TGeoHMatrix* GetPHOSMatrix(Int_t i) const {return fESDRun->GetPHOSMatrix(i);}
+  void     SetEMCALMatrix(TGeoHMatrix*matrix, Int_t i) {fESDRun->SetEMCALMatrix(matrix,i);}
+  const TGeoHMatrix* GetEMCALMatrix(Int_t i) const {return fESDRun->GetEMCALMatrix(i);}
+  //
+  void        SetCurrentL3(Float_t cur)           const  {fESDRun->SetCurrentL3(cur);}
+  void        SetCurrentDip(Float_t cur)          const  {fESDRun->SetCurrentDip(cur);}
+  void        SetBeamEnergy(Float_t be)           const  {fESDRun->SetBeamEnergy(be);}
+  void        SetBeamType(const char* bt)         const  {fESDRun->SetBeamType(bt);}
+  void        SetUniformBMap(Bool_t val=kTRUE)    const  {fESDRun->SetBit(AliESDRun::kUniformBMap,val);}
+  void        SetBInfoStored(Bool_t val=kTRUE)    const  {fESDRun->SetBit(AliESDRun::kBInfoStored,val);}
+  //
+  Float_t     GetCurrentL3()                      const  {return fESDRun->GetCurrentL3();}
+  Float_t     GetCurrentDip()                     const  {return fESDRun->GetCurrentDip();}
+  Float_t     GetBeamEnergy()                     const  {return fESDRun->GetBeamEnergy();}
+  const char* GetBeamType()                       const  {return fESDRun->GetBeamType();}
+  Bool_t      IsUniformBMap()                     const  {return fESDRun->TestBit(AliESDRun::kUniformBMap);}
+  //
+  Bool_t      InitMagneticField()                 const  {return fESDRun->InitMagneticField();} 
   // HEADER
-  const AliESDHeader* GetHeader() const {return fHeader;}
+  AliESDHeader* GetHeader() const {return fHeader;}
 
   // Delegated methods for fHeader
   void      SetTriggerMask(ULong64_t n) {fHeader->SetTriggerMask(n);}
   void      SetOrbitNumber(UInt_t n) {fHeader->SetOrbitNumber(n);}
   void      SetTimeStamp(UInt_t timeStamp){fHeader->SetTimeStamp(timeStamp);}
   void      SetEventType(UInt_t eventType){fHeader->SetEventType(eventType);}
+  void      SetEventSpecie(UInt_t eventSpecie){fHeader->SetEventSpecie(eventSpecie);}
   void      SetEventNumberInFile(Int_t n) {fHeader->SetEventNumberInFile(n);}
-  //  void      SetRunNumber(Int_t n) {fHeader->SetRunNumber(n);}
+  //  void     SetRunNumber(Int_t n) {fHeader->SetRunNumber(n);}
   void      SetBunchCrossNumber(UShort_t n) {fHeader->SetBunchCrossNumber(n);}
   void      SetTriggerCluster(UChar_t n) {fHeader->SetTriggerCluster(n);}
+  
   ULong64_t GetTriggerMask() const {return fHeader->GetTriggerMask();}
+  TString   GetFiredTriggerClasses() const {return fESDRun->GetFiredTriggerClasses(fHeader->GetTriggerMask());}
+  Bool_t    IsTriggerClassFired(const char *name) const {return fESDRun->IsTriggerClassFired(fHeader->GetTriggerMask(),name);}
+  Bool_t    IsEventSelected(const char *trigExpr) const;
+  TObject*  GetHLTTriggerDecision() const;
+  TString   GetHLTTriggerDescription() const;
+  Bool_t    IsHLTTriggerFired(const char* name=NULL) const;
   UInt_t    GetOrbitNumber() const {return fHeader->GetOrbitNumber();}
   UInt_t    GetTimeStamp()  const { return fHeader->GetTimeStamp();}
   UInt_t    GetEventType()  const { return fHeader->GetEventType();}
+  UInt_t    GetEventSpecie()  const { return fHeader->GetEventSpecie();}
   Int_t     GetEventNumberInFile() const {return fHeader->GetEventNumberInFile();}
   UShort_t  GetBunchCrossNumber() const {return fHeader->GetBunchCrossNumber();}
   UChar_t   GetTriggerCluster() const {return fHeader->GetTriggerCluster();}
 
   // ZDC CKB: put this in the header?
-  const AliESDZDC*    GetESDZDC() const {return fESDZDC;}
+  AliESDZDC*    GetESDZDC() const {return fESDZDC;}
 
   // Delegated methods for fESDZDC
-  Float_t GetZDCN1Energy() const {return fESDZDC->GetZDCN1Energy();}
-  Float_t GetZDCP1Energy() const {return fESDZDC->GetZDCP1Energy();}
-  Float_t GetZDCN2Energy() const {return fESDZDC->GetZDCN2Energy();}
-  Float_t GetZDCP2Energy() const {return fESDZDC->GetZDCP2Energy();}
-  Float_t GetZDCEMEnergy() const {return fESDZDC->GetZDCEMEnergy();}
-  Int_t   GetZDCParticipants() const {return fESDZDC->GetZDCParticipants();}
-  void    SetZDC(Float_t n1Energy, Float_t p1Energy, Float_t emEnergy,
-                 Float_t n2Energy, Float_t p2Energy, Int_t participants)
-  {fESDZDC->SetZDC(n1Energy, p1Energy, emEnergy, n2Energy, p2Energy, participants);}
+  Double_t GetZDCN1Energy() const {return fESDZDC->GetZDCN1Energy();}
+  Double_t GetZDCP1Energy() const {return fESDZDC->GetZDCP1Energy();}
+  Double_t GetZDCN2Energy() const {return fESDZDC->GetZDCN2Energy();}
+  Double_t GetZDCP2Energy() const {return fESDZDC->GetZDCP2Energy();}
+  Double_t GetZDCEMEnergy(Int_t i=0) const {return fESDZDC->GetZDCEMEnergy(i);}
+  Int_t    GetZDCParticipants() const {return fESDZDC->GetZDCParticipants();}
+  void     SetZDC(Float_t n1Energy, Float_t p1Energy, Float_t em1Energy, Float_t em2Energy,
+                  Float_t n2Energy, Float_t p2Energy, Int_t participants, Int_t nPartA,
+                 Int_t nPartC, Double_t b, Double_t bA, Double_t bC, UInt_t recoflag)
+           {fESDZDC->SetZDC(n1Energy, p1Energy, em1Energy, em2Energy, n2Energy, p2Energy, 
+            participants, nPartA, nPartC, b, bA, bC,  recoflag);}
+  void     SetZDCScaler(UInt_t *counts) {fESDZDC->SetZDCScaler(counts);}
 
 
   // FMD
   void SetFMDData(AliESDFMD * obj);
-  AliESDFMD *GetFMDData(){ return fESDFMD; }
+  AliESDFMD *GetFMDData() const { return fESDFMD; }
 
 
   // TZERO CKB: put this in the header?
   const AliESDTZERO*    GetESDTZERO() const {return fESDTZERO;}
   // delegetated methods for fESDTZERO
 
-  Float_t GetT0zVertex() const {return fESDTZERO->GetT0zVertex();}
+  Double_t GetT0zVertex() const {return fESDTZERO->GetT0zVertex();}
   void SetT0zVertex(Float_t z) {fESDTZERO->SetT0zVertex(z);}
-  Float_t GetT0() const {return fESDTZERO->GetT0();}
+  Double_t GetT0() const {return fESDTZERO->GetT0();}
   void SetT0(Float_t timeStart) {fESDTZERO->SetT0(timeStart);}
-  const Float_t * GetT0time() const {return fESDTZERO->GetT0time();}
+  Float_t GetT0clock() const {return fESDTZERO->GetT0clock();}
+  void SetT0clock(Float_t timeStart) {fESDTZERO->SetT0clock(timeStart);}
+  const Double_t * GetT0time() const {return fESDTZERO->GetT0time();}
   void SetT0time(Float_t time[24]) {fESDTZERO->SetT0time(time);}
-  const Float_t * GetT0amplitude() const {return fESDTZERO->GetT0amplitude();}
+  const Double_t * GetT0amplitude() const {return fESDTZERO->GetT0amplitude();}
   void SetT0amplitude(Float_t amp[24]){fESDTZERO->SetT0amplitude(amp);}
 
   // VZERO 
   AliESDVZERO *GetVZEROData() const { return fESDVZERO; }
   void SetVZEROData(AliESDVZERO * obj);
 
+ // ACORDE
+  AliESDACORDE *GetACORDEData() const { return fESDACORDE;}
+  void SetACORDEData(AliESDACORDE * obj);
 
-  void SetESDfriend(const AliESDfriend *f);
+  void SetESDfriend(const AliESDfriend *f) const;
   void GetESDfriend(AliESDfriend *f) const;
 
 
 
-  void SetVertex(const AliESDVertex *vertex); 
-  const AliESDVertex *GetVertex() const {return fSPDVertex;}
+  void SetPrimaryVertexTPC(const AliESDVertex *vertex); 
+  const AliESDVertex *GetPrimaryVertexTPC() const {return fTPCVertex;}
 
-  void SetPrimaryVertex(const AliESDVertex *vertex);
-  const AliESDVertex *GetPrimaryVertex() const {return fPrimaryVertex;}
+  void SetPrimaryVertexSPD(const AliESDVertex *vertex); 
+  const AliESDVertex *GetPrimaryVertexSPD() const {return fSPDVertex;}
+  const AliESDVertex *GetVertex() const {
+    //For the backward compatibily only
+     return GetPrimaryVertexSPD();
+  }
+
+  void SetPrimaryVertexTracks(const AliESDVertex *vertex);
+  const AliESDVertex *GetPrimaryVertexTracks() const {return fPrimaryVertex;}
+  AliESDVertex *PrimaryVertexTracksUnconstrained() const;
+
+  const AliESDVertex *GetPrimaryVertex() const;
 
   void SetMultiplicity(const AliMultiplicity *mul);
 
   const AliMultiplicity *GetMultiplicity() const {return fSPDMult;}
 
 
-  
+  Bool_t Clean(Float_t *cleanPars);
+  Bool_t RemoveKink(Int_t i)   const;
+  Bool_t RemoveV0(Int_t i)     const;
+  Bool_t RemoveTrack(Int_t i)  const;
+
+  const AliESDVertex *GetPileupVertexSPD(Int_t i) const {
+    return (const AliESDVertex *)fSPDPileupVertices->UncheckedAt(i);
+  }
+  Char_t  AddPileupVertexSPD(const AliESDVertex *vtx);
+  Bool_t  IsPileupFromSPD(Int_t ncont=2, Double_t nSigmaDeltaZ=3., Double_t nSigmaXY=2., Int_t option=0) const;
+
+  const AliESDVertex *GetPileupVertexTracks(Int_t i) const {
+    return (const AliESDVertex *)fTrkPileupVertices->UncheckedAt(i);
+  }
+  Char_t  AddPileupVertexTracks(const AliESDVertex *vtx);
+
   AliESDtrack *GetTrack(Int_t i) const {
     return (AliESDtrack *)fTracks->UncheckedAt(i);
   }
@@ -240,12 +316,16 @@ public:
   AliESDCaloCluster *GetCaloCluster(Int_t i) const {
     return (AliESDCaloCluster *)fCaloClusters->UncheckedAt(i);
   }
+
   Int_t AddCaloCluster(const AliESDCaloCluster *c);
 
+  AliESDCaloCells *GetEMCALCells() const {return fEMCALCells; }  
+  AliESDCaloCells *GetPHOSCells() const {return fPHOSCells; }  
+
   AliRawDataErrorLog *GetErrorLog(Int_t i) const {
     return (AliRawDataErrorLog *)fErrorLogs->UncheckedAt(i);
   }
-  void  AddRawDataErrorLog(const AliRawDataErrorLog *log);
+  void  AddRawDataErrorLog(const AliRawDataErrorLog *log) const;
 
   Int_t GetNumberOfErrorLogs()   const {return fErrorLogs->GetEntriesFast();}
 
@@ -255,6 +335,12 @@ public:
   void AddEMCALTriggerPosition(TArrayF array)  { fEMCALTrigger->AddTriggerPosition(array); }
   void AddEMCALTriggerAmplitudes(TArrayF array){ fEMCALTrigger->AddTriggerAmplitudes(array); }
 
+  Int_t GetNumberOfPileupVerticesSPD() const {
+    return (fSPDPileupVertices?fSPDPileupVertices->GetEntriesFast():0);
+  }
+  Int_t GetNumberOfPileupVerticesTracks() const {
+    return (fTrkPileupVertices?fTrkPileupVertices->GetEntriesFast():0);
+  }
   Int_t GetNumberOfTracks()     const {return fTracks->GetEntriesFast();}
   Int_t GetNumberOfHLTConfMapTracks()     const {return 0;} 
   // fHLTConfMapTracks->GetEntriesFast();}
@@ -267,19 +353,29 @@ public:
   Int_t GetNumberOfV0s()      const {return fV0s->GetEntriesFast();}
   Int_t GetNumberOfCascades() const {return fCascades->GetEntriesFast();}
   Int_t GetNumberOfKinks() const {return fKinks->GetEntriesFast();}
+  
+  Int_t GetEMCALClusters(TRefArray *clusters) const;
+  Int_t GetPHOSClusters(TRefArray *clusters) const;
   Int_t GetNumberOfCaloClusters() const {return fCaloClusters->GetEntriesFast();}
 
+  void SetUseOwnList(Bool_t b){fUseOwnList = b;}
+  Bool_t GetUseOwnList() const {return fUseOwnList;}
+  
+  // Remove this stuff CKB?
+  //---------------------------------------------------
   Int_t GetNumberOfEMCALClusters() const {return fEMCALClusters;}
   void  SetNumberOfEMCALClusters(Int_t clus) {fEMCALClusters = clus;}
   Int_t GetFirstEMCALCluster() const {return fFirstEMCALCluster;}
   void  SetFirstEMCALCluster(Int_t index) {fFirstEMCALCluster = index;}
-  TArrayF *GetEMCALTriggerPosition() const {return  fEMCALTrigger->GetTriggerPosition();}
-  TArrayF *GetEMCALTriggerAmplitudes() const {return  fEMCALTrigger->GetTriggerAmplitudes();}
-
   Int_t GetNumberOfPHOSClusters() const {return fPHOSClusters;}
   void  SetNumberOfPHOSClusters(Int_t part) { fPHOSClusters = part ; }
   void  SetFirstPHOSCluster(Int_t index) { fFirstPHOSCluster = index ; } 
   Int_t GetFirstPHOSCluster() const  { return fFirstPHOSCluster ; }
+  //-------------------------------------------------------
+
+  TArrayF *GetEMCALTriggerPosition() const {return  fEMCALTrigger->GetTriggerPosition();}
+  TArrayF *GetEMCALTriggerAmplitudes() const {return  fEMCALTrigger->GetTriggerAmplitudes();}
   TArrayF *GetPHOSTriggerPosition() const {return  fPHOSTrigger->GetTriggerPosition();}
   TArrayF *GetPHOSTriggerAmplitudes() const {return  fPHOSTrigger->GetTriggerAmplitudes();}
 
@@ -290,21 +386,21 @@ public:
   void  Print(Option_t *option="") const;
 
   void AddObject(TObject* obj);
-  void ReadFromTree(TTree *tree);
+  void ReadFromTree(TTree *tree, Option_t* opt = "");
   TObject* FindListObject(const char *name);
   AliESD *GetAliESDOld(){return fESDOld;}
-  const void WriteToTree(TTree* tree) const {tree->Branch(fESDObjects);}
+  void WriteToTree(TTree* tree) const;
   void GetStdContent();
   void ResetStdContent();
   void CreateStdContent();
+  void CreateStdContent(Bool_t bUseThisList);
   void SetStdNames();
   void CopyFromOldESD();
-  TList* GetList(){return fESDObjects;}
+  TList* GetList() const {return fESDObjects;}
 
 protected:
   AliESDEvent(const AliESDEvent&);
-  AliESDEvent &operator=(const AliESDEvent& source);
-
+  static Bool_t ResetWithPlacementNew(TObject *pObject);
 
   TList *fESDObjects;             // List of esd Objects
 
@@ -314,12 +410,16 @@ protected:
   AliESDFMD       *fESDFMD;           //! FMD object containing rough multiplicity
   AliESDVZERO     *fESDVZERO;         //! VZERO object containing rough multiplicity
   AliESDTZERO     *fESDTZERO;         //! TZEROObject
+  AliESDVertex    *fTPCVertex;        //! Primary vertex estimated by the TPC
   AliESDVertex    *fSPDVertex;        //! Primary vertex estimated by the SPD
   AliESDVertex    *fPrimaryVertex;    //! Primary vertex estimated using ESD tracks
   AliMultiplicity *fSPDMult;          //! SPD tracklet multiplicity
   AliESDCaloTrigger* fPHOSTrigger;     //! PHOS Trigger information
   AliESDCaloTrigger* fEMCALTrigger;    //! PHOS Trigger information
+  AliESDACORDE    *fESDACORDE;        //! ACORDE ESD object caontaining bit pattern
 
+  TClonesArray *fSPDPileupVertices;//! Pileup primary vertices reconstructed by SPD 
+  TClonesArray *fTrkPileupVertices;//! Pileup primary vertices reconstructed using the tracks 
   TClonesArray *fTracks;           //! ESD tracks 
   TClonesArray *fMuonTracks;       //! MUON ESD tracks
   TClonesArray *fPmdTracks;        //! PMD ESD tracks
@@ -328,14 +428,18 @@ protected:
   TClonesArray *fCascades;         //! Cascade vertices
   TClonesArray *fKinks;            //! Kinks
   TClonesArray *fCaloClusters;     //! Calorimeter clusters for PHOS/EMCAL
+  AliESDCaloCells *fEMCALCells;     //! EMCAL cell info
+  AliESDCaloCells *fPHOSCells;     //! PHOS cell info
   TClonesArray *fErrorLogs;        //! Raw-data reading error messages
  
 
 
-  AliESD    *fESDOld;              //! Old esd Structure
+  AliESD       *fESDOld;           //! Old esd Structure
+  AliESDfriend *fESDFriendOld;     //! Old friend esd Structure
   Bool_t    fConnected;            //! flag if leaves are alreday connected
+  Bool_t    fUseOwnList;           //! Do not use the list from the esdTree but use the one created by this class 
 
-  static const char* fESDListName[kESDListN];
+  static const char* fgkESDListName[kESDListN]; //!
 
   // Remove this stuff CKB
   Int_t        fEMCALClusters;   // Number of EMCAL clusters (subset of caloclusters)
@@ -344,7 +448,7 @@ protected:
   Int_t        fPHOSClusters;     // Number of PHOS clusters (subset of caloclusters)
   Int_t        fFirstPHOSCluster; // First PHOS cluster in the fCaloClusters list 
 
-  ClassDef(AliESDEvent,2)  //ESDEvent class 
+  ClassDef(AliESDEvent,10)  //ESDEvent class 
 };
 #endif