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 #include "AliESDACORDE.h"
  33
  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"
  39
  40 class AliESDfriend;
  41 class AliESDVZERO;
  42 class AliESDHLTtrack;
  43 class AliESDVertex;
  44 class AliESDPmdTrack;
  45 class AliESDFMD;
  46 class AliESDkink;
  47 class AliESDCaloCluster;
  48 class AliESDCaloCells;
  49 class AliESDv0;
  50 class AliMultiplicity;
  51 class AliRawDataErrorLog;
  52 class AliESDRun;
  53 class AliESDTrdTrack;
  54 class AliESDMuonTrack;
  55 class AliESD;
  56 class AliESDcascade;
  57 class TRefArray;
  58 class AliESDACORDE;
  59
  60 class AliESDEvent : public AliVEvent {
  61 public:
  62
  63
  64   enum ESDListIndex   {kESDRun,
  65                        kHeader,
  66                        kESDZDC,
  67                        kESDFMD,
  68                        kESDVZERO,
  69                        kESDTZERO,
  70                        kTPCVertex,
  71                        kSPDVertex,
  72                        kPrimaryVertex,
  73                        kSPDMult,
  74                        kPHOSTrigger,
  75                        kEMCALTrigger,
  76                        kSPDPileupVertices,
  77                        kTrkPileupVertices,
  78                        kTracks,
  79                        kMuonTracks,
  80                        kPmdTracks,
  81                        kTrdTracks,
  82                        kV0s,
  83                        kCascades,
  84                        kKinks,
  85                        kCaloClusters,
  86                        kEMCALCells,
  87                        kPHOSCells,
  88                        kErrorLogs,
  89                        kESDACORDE,
  90                        kESDListN
  91   };
  92
  93   AliESDEvent();
  94   virtual ~AliESDEvent();
  95   AliESDEvent &operator=(const AliESDEvent& source); // or make private and use only copy?
  96   virtual void Copy(TObject& obj) const;
  97
  98   // RUN
  99   // move this to the UserData!!!
 100   const AliESDRun*    GetESDRun() const {return fESDRun;}
 101
 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);}
 120
 121   // HEADER
 122   AliESDHeader* GetHeader() const {return fHeader;}
 123
 124   // Delegated methods for fHeader
 125   void      SetTriggerMask(ULong64_t n) {fHeader->SetTriggerMask(n);}
 126   void      SetOrbitNumber(UInt_t n) {fHeader->SetOrbitNumber(n);}
 127   void      SetTimeStamp(UInt_t timeStamp){fHeader->SetTimeStamp(timeStamp);}
 128   void      SetEventType(UInt_t eventType){fHeader->SetEventType(eventType);}
 129   void      SetEventNumberInFile(Int_t n) {fHeader->SetEventNumberInFile(n);}
 130   //  void     SetRunNumber(Int_t n) {fHeader->SetRunNumber(n);}
 131   void      SetBunchCrossNumber(UShort_t n) {fHeader->SetBunchCrossNumber(n);}
 132   void      SetTriggerCluster(UChar_t n) {fHeader->SetTriggerCluster(n);}
 133
 134   ULong64_t GetTriggerMask() const {return fHeader->GetTriggerMask();}
 135   TString   GetFiredTriggerClasses() const {return fESDRun->GetFiredTriggerClasses(fHeader->GetTriggerMask());}
 136   Bool_t    IsTriggerClassFired(const char *name) const {return fESDRun->IsTriggerClassFired(fHeader->GetTriggerMask(),name);}
 137   Bool_t    IsEventSelected(const char *trigExpr) const;
 138   TObject*  GetHLTTriggerDecision() const;
 139   TString   GetHLTTriggerDescription() const;
 140   Bool_t    IsHLTTriggerFired(const char* name=NULL) const;
 141   UInt_t    GetOrbitNumber() const {return fHeader->GetOrbitNumber();}
 142   UInt_t    GetTimeStamp()  const { return fHeader->GetTimeStamp();}
 143   UInt_t    GetEventType()  const { return fHeader->GetEventType();}
 144   Int_t     GetEventNumberInFile() const {return fHeader->GetEventNumberInFile();}
 145   UShort_t  GetBunchCrossNumber() const {return fHeader->GetBunchCrossNumber();}
 146   UChar_t   GetTriggerCluster() const {return fHeader->GetTriggerCluster();}
 147
 148   // ZDC CKB: put this in the header?
 149   AliESDZDC*    GetESDZDC() const {return fESDZDC;}
 150
 151   // Delegated methods for fESDZDC
 152   Double_t GetZDCN1Energy() const {return fESDZDC->GetZDCN1Energy();}
 153   Double_t GetZDCP1Energy() const {return fESDZDC->GetZDCP1Energy();}
 154   Double_t GetZDCN2Energy() const {return fESDZDC->GetZDCN2Energy();}
 155   Double_t GetZDCP2Energy() const {return fESDZDC->GetZDCP2Energy();}
 156   Double_t GetZDCEMEnergy(Int_t i=0) const {return fESDZDC->GetZDCEMEnergy(i);}
 157   Int_t   GetZDCParticipants() const {return fESDZDC->GetZDCParticipants();}
 158   void    SetZDC(Float_t n1Energy, Float_t p1Energy, Float_t em1Energy, Float_t em2Energy,
 159                  Float_t n2Energy, Float_t p2Energy, Int_t participants, Int_t nPartA,
 160                  Int_t nPartC, Double_t b, Double_t bA, Double_t bC, UInt_t recoflag)
 161   {fESDZDC->SetZDC(n1Energy, p1Energy, em1Energy, em2Energy, n2Energy, p2Energy,
 162    participants, nPartA, nPartC, b, bA, bC,  recoflag);}
 163
 164
 165   // FMD
 166   void SetFMDData(AliESDFMD * obj);
 167   AliESDFMD *GetFMDData(){ return fESDFMD; }
 168
 169
 170   // TZERO CKB: put this in the header?
 171   const AliESDTZERO*    GetESDTZERO() const {return fESDTZERO;}
 172   // delegetated methods for fESDTZERO
 173
 174   Double_t GetT0zVertex() const {return fESDTZERO->GetT0zVertex();}
 175   void SetT0zVertex(Float_t z) {fESDTZERO->SetT0zVertex(z);}
 176   Double_t GetT0() const {return fESDTZERO->GetT0();}
 177   void SetT0(Float_t timeStart) {fESDTZERO->SetT0(timeStart);}
 178   Float_t GetT0clock() const {return fESDTZERO->GetT0clock();}
 179   void SetT0clock(Float_t timeStart) {fESDTZERO->SetT0clock(timeStart);}
 180   const Double_t * GetT0time() const {return fESDTZERO->GetT0time();}
 181   void SetT0time(Float_t time[24]) {fESDTZERO->SetT0time(time);}
 182   const Double_t * GetT0amplitude() const {return fESDTZERO->GetT0amplitude();}
 183   void SetT0amplitude(Float_t amp[24]){fESDTZERO->SetT0amplitude(amp);}
 184
 185   // VZERO
 186   AliESDVZERO *GetVZEROData() const { return fESDVZERO; }
 187   void SetVZEROData(AliESDVZERO * obj);
 188
 189  // ACORDE
 190   AliESDACORDE *GetACORDEData() const { return fESDACORDE;}
 191   void SetACORDEData(AliESDACORDE * obj);
 192
 193   void SetESDfriend(const AliESDfriend *f) const;
 194   void GetESDfriend(AliESDfriend *f) const;
 195
 196
 197
 198   void SetPrimaryVertexTPC(const AliESDVertex *vertex);
 199   const AliESDVertex *GetPrimaryVertexTPC() const {return fTPCVertex;}
 200
 201   void SetPrimaryVertexSPD(const AliESDVertex *vertex);
 202   const AliESDVertex *GetPrimaryVertexSPD() const {return fSPDVertex;}
 203   const AliESDVertex *GetVertex() const {
 204     //For the backward compatibily only
 205      return GetPrimaryVertexSPD();
 206   }
 207
 208   void SetPrimaryVertexTracks(const AliESDVertex *vertex);
 209   const AliESDVertex *GetPrimaryVertexTracks() const {return fPrimaryVertex;}
 210   AliESDVertex *PrimaryVertexTracksUnconstrained() const;
 211
 212   const AliESDVertex *GetPrimaryVertex() const;
 213
 214   void SetMultiplicity(const AliMultiplicity *mul);
 215
 216   const AliMultiplicity *GetMultiplicity() const {return fSPDMult;}
 217
 218
 219   Bool_t Clean(Float_t *cleanPars);
 220   Bool_t RemoveKink(Int_t i)   const;
 221   Bool_t RemoveV0(Int_t i)     const;
 222   Bool_t RemoveTrack(Int_t i)  const;
 223
 224   const AliESDVertex *GetPileupVertexSPD(Int_t i) const {
 225     return (const AliESDVertex *)fSPDPileupVertices->UncheckedAt(i);
 226   }
 227   Char_t  AddPileupVertexSPD(const AliESDVertex *vtx);
 228
 229   const AliESDVertex *GetPileupVertexTracks(Int_t i) const {
 230     return (const AliESDVertex *)fTrkPileupVertices->UncheckedAt(i);
 231   }
 232   Char_t  AddPileupVertexTracks(const AliESDVertex *vtx);
 233
 234   AliESDtrack *GetTrack(Int_t i) const {
 235     return (AliESDtrack *)fTracks->UncheckedAt(i);
 236   }
 237   Int_t  AddTrack(const AliESDtrack *t);
 238
 239
 240   AliESDHLTtrack *GetHLTConfMapTrack(Int_t /*i*/) const {
 241     //    return (AliESDHLTtrack *)fHLTConfMapTracks->UncheckedAt(i);
 242     return 0;
 243   }
 244   void AddHLTConfMapTrack(const AliESDHLTtrack */*t*/) {
 245     printf("ESD:: AddHLTConfMapTrack do nothing \n");
 246     //    TClonesArray &fhlt = *fHLTConfMapTracks;
 247     //  new(fhlt[fHLTConfMapTracks->GetEntriesFast()]) AliESDHLTtrack(*t);
 248   }
 249
 250
 251   AliESDHLTtrack *GetHLTHoughTrack(Int_t /*i*/) const {
 252     //    return (AliESDHLTtrack *)fHLTHoughTracks->UncheckedAt(i);
 253     return 0;
 254   }
 255   void AddHLTHoughTrack(const AliESDHLTtrack */*t*/) {
 256     printf("ESD:: AddHLTHoughTrack do nothing \n");
 257     //    TClonesArray &fhlt = *fHLTHoughTracks;
 258     //     new(fhlt[fHLTHoughTracks->GetEntriesFast()]) AliESDHLTtrack(*t);
 259   }
 260
 261   AliESDMuonTrack *GetMuonTrack(Int_t i) const {
 262     return (AliESDMuonTrack *)fMuonTracks->UncheckedAt(i);
 263   }
 264
 265   void AddMuonTrack(const AliESDMuonTrack *t);
 266
 267   AliESDPmdTrack *GetPmdTrack(Int_t i) const {
 268     return (AliESDPmdTrack *)fPmdTracks->UncheckedAt(i);
 269   }
 270
 271   void AddPmdTrack(const AliESDPmdTrack *t);
 272
 273
 274   AliESDTrdTrack *GetTrdTrack(Int_t i) const {
 275     return (AliESDTrdTrack *)fTrdTracks->UncheckedAt(i);
 276   }
 277
 278
 279   void AddTrdTrack(const AliESDTrdTrack *t);
 280
 281   AliESDv0 *GetV0(Int_t i) const {
 282     return (AliESDv0*)fV0s->UncheckedAt(i);
 283   }
 284   Int_t AddV0(const AliESDv0 *v);
 285
 286   AliESDcascade *GetCascade(Int_t i) const {
 287     return (AliESDcascade *)fCascades->UncheckedAt(i);
 288   }
 289
 290   void AddCascade(const AliESDcascade *c);
 291
 292   AliESDkink *GetKink(Int_t i) const {
 293     return (AliESDkink *)fKinks->UncheckedAt(i);
 294   }
 295   Int_t AddKink(const AliESDkink *c);
 296
 297   AliESDCaloCluster *GetCaloCluster(Int_t i) const {
 298     return (AliESDCaloCluster *)fCaloClusters->UncheckedAt(i);
 299   }
 300
 301   Int_t AddCaloCluster(const AliESDCaloCluster *c);
 302
 303   AliESDCaloCells *GetEMCALCells() const {return fEMCALCells; }
 304   AliESDCaloCells *GetPHOSCells() const {return fPHOSCells; }
 305
 306   AliRawDataErrorLog *GetErrorLog(Int_t i) const {
 307     return (AliRawDataErrorLog *)fErrorLogs->UncheckedAt(i);
 308   }
 309   void  AddRawDataErrorLog(const AliRawDataErrorLog *log) const;
 310
 311   Int_t GetNumberOfErrorLogs()   const {return fErrorLogs->GetEntriesFast();}
 312
 313
 314   void AddPHOSTriggerPosition(TArrayF array)   { fPHOSTrigger->AddTriggerPosition(array); }
 315   void AddPHOSTriggerAmplitudes(TArrayF array) { fPHOSTrigger->AddTriggerAmplitudes(array);}
 316   void AddEMCALTriggerPosition(TArrayF array)  { fEMCALTrigger->AddTriggerPosition(array); }
 317   void AddEMCALTriggerAmplitudes(TArrayF array){ fEMCALTrigger->AddTriggerAmplitudes(array); }
 318
 319   Int_t GetNumberOfPileupVerticesSPD() const {
 320     return (fSPDPileupVertices?fSPDPileupVertices->GetEntriesFast():0);
 321   }
 322   Int_t GetNumberOfPileupVerticesTracks() const {
 323     return (fTrkPileupVertices?fTrkPileupVertices->GetEntriesFast():0);
 324   }
 325   Int_t GetNumberOfTracks()     const {return fTracks->GetEntriesFast();}
 326   Int_t GetNumberOfHLTConfMapTracks()     const {return 0;}
 327   // fHLTConfMapTracks->GetEntriesFast();}
 328   Int_t GetNumberOfHLTHoughTracks()     const {return  0;  }
 329   //  fHLTHoughTracks->GetEntriesFast();  }
 330
 331   Int_t GetNumberOfMuonTracks() const {return fMuonTracks->GetEntriesFast();}
 332   Int_t GetNumberOfPmdTracks() const {return fPmdTracks->GetEntriesFast();}
 333   Int_t GetNumberOfTrdTracks() const {return fTrdTracks->GetEntriesFast();}
 334   Int_t GetNumberOfV0s()      const {return fV0s->GetEntriesFast();}
 335   Int_t GetNumberOfCascades() const {return fCascades->GetEntriesFast();}
 336   Int_t GetNumberOfKinks() const {return fKinks->GetEntriesFast();}
 337
 338   Int_t GetEMCALClusters(TRefArray *clusters) const;
 339   Int_t GetPHOSClusters(TRefArray *clusters) const;
 340   Int_t GetNumberOfCaloClusters() const {return fCaloClusters->GetEntriesFast();}
 341
 342   void SetUseOwnList(Bool_t b){fUseOwnList = b;}
 343   Bool_t GetUseOwnList(){return fUseOwnList;}
 344
 345   // Remove this stuff CKB?
 346   //---------------------------------------------------
 347   Int_t GetNumberOfEMCALClusters() const {return fEMCALClusters;}
 348   void  SetNumberOfEMCALClusters(Int_t clus) {fEMCALClusters = clus;}
 349   Int_t GetFirstEMCALCluster() const {return fFirstEMCALCluster;}
 350   void  SetFirstEMCALCluster(Int_t index) {fFirstEMCALCluster = index;}
 351
 352   Int_t GetNumberOfPHOSClusters() const {return fPHOSClusters;}
 353   void  SetNumberOfPHOSClusters(Int_t part) { fPHOSClusters = part ; }
 354   void  SetFirstPHOSCluster(Int_t index) { fFirstPHOSCluster = index ; }
 355   Int_t GetFirstPHOSCluster() const  { return fFirstPHOSCluster ; }
 356   //-------------------------------------------------------
 357
 358   TArrayF *GetEMCALTriggerPosition() const {return  fEMCALTrigger->GetTriggerPosition();}
 359   TArrayF *GetEMCALTriggerAmplitudes() const {return  fEMCALTrigger->GetTriggerAmplitudes();}
 360   TArrayF *GetPHOSTriggerPosition() const {return  fPHOSTrigger->GetTriggerPosition();}
 361   TArrayF *GetPHOSTriggerAmplitudes() const {return  fPHOSTrigger->GetTriggerAmplitudes();}
 362
 363   void ResetV0s() { fV0s->Clear(); }
 364   void ResetCascades() { fCascades->Clear(); }
 365   void Reset();
 366
 367   void  Print(Option_t *option="") const;
 368
 369   void AddObject(TObject* obj);
 370   void ReadFromTree(TTree *tree, Option_t* opt = "");
 371   TObject* FindListObject(const char *name);
 372   AliESD *GetAliESDOld(){return fESDOld;}
 373   void WriteToTree(TTree* tree) const;
 374   void GetStdContent();
 375   void ResetStdContent();
 376   void CreateStdContent();
 377   void CreateStdContent(Bool_t bUseThisList);
 378   void SetStdNames();
 379   void CopyFromOldESD();
 380   TList* GetList() const {return fESDObjects;}
 381
 382 protected:
 383   AliESDEvent(const AliESDEvent&);
 384   static Bool_t ResetWithPlacementNew(TObject *pObject);
 385
 386   TList *fESDObjects;             // List of esd Objects
 387
 388   AliESDRun       *fESDRun;           //! Run information tmp put in the Userdata
 389   AliESDHeader    *fHeader;           //! ESD Event Header
 390   AliESDZDC       *fESDZDC;           //! ZDC information
 391   AliESDFMD       *fESDFMD;           //! FMD object containing rough multiplicity
 392   AliESDVZERO     *fESDVZERO;         //! VZERO object containing rough multiplicity
 393   AliESDTZERO     *fESDTZERO;         //! TZEROObject
 394   AliESDVertex    *fTPCVertex;        //! Primary vertex estimated by the TPC
 395   AliESDVertex    *fSPDVertex;        //! Primary vertex estimated by the SPD
 396   AliESDVertex    *fPrimaryVertex;    //! Primary vertex estimated using ESD tracks
 397   AliMultiplicity *fSPDMult;          //! SPD tracklet multiplicity
 398   AliESDCaloTrigger* fPHOSTrigger;     //! PHOS Trigger information
 399   AliESDCaloTrigger* fEMCALTrigger;    //! PHOS Trigger information
 400   AliESDACORDE    *fESDACORDE;        //! ACORDE ESD object caontaining bit pattern
 401
 402   TClonesArray *fSPDPileupVertices;//! Pileup primary vertices reconstructed by SPD
 403   TClonesArray *fTrkPileupVertices;//! Pileup primary vertices reconstructed using the tracks
 404   TClonesArray *fTracks;           //! ESD tracks
 405   TClonesArray *fMuonTracks;       //! MUON ESD tracks
 406   TClonesArray *fPmdTracks;        //! PMD ESD tracks
 407   TClonesArray *fTrdTracks;        //! TRD ESD tracks (triggered)
 408   TClonesArray *fV0s;              //! V0 vertices
 409   TClonesArray *fCascades;         //! Cascade vertices
 410   TClonesArray *fKinks;            //! Kinks
 411   TClonesArray *fCaloClusters;     //! Calorimeter clusters for PHOS/EMCAL
 412   AliESDCaloCells *fEMCALCells;     //! EMCAL cell info
 413   AliESDCaloCells *fPHOSCells;     //! PHOS cell info
 414   TClonesArray *fErrorLogs;        //! Raw-data reading error messages
 415
 416
 417
 418   AliESD       *fESDOld;           //! Old esd Structure
 419   AliESDfriend *fESDFriendOld;     //! Old friend esd Structure
 420   Bool_t    fConnected;            //! flag if leaves are alreday connected
 421   Bool_t    fUseOwnList;           //! Do not use the list from the esdTree but use the one created by this class
 422
 423   static const char* fgkESDListName[kESDListN]; //!
 424
 425   // Remove this stuff CKB
 426   Int_t        fEMCALClusters;   // Number of EMCAL clusters (subset of caloclusters)
 427   Int_t        fFirstEMCALCluster; // First EMCAL cluster in the fCaloClusters list
 428
 429   Int_t        fPHOSClusters;     // Number of PHOS clusters (subset of caloclusters)
 430   Int_t        fFirstPHOSCluster; // First PHOS cluster in the fCaloClusters list
 431
 432   ClassDef(AliESDEvent,10)  //ESDEvent class
 433 };
 434 #endif
 435