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