Add inheritance from AliVEvent. Fix standard naming scheme. Fix to be able to ReadFro...
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 7 Aug 2007 20:54:04 +0000 (20:54 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 7 Aug 2007 20:54:04 +0000 (20:54 +0000)
STEER/AliAODEvent.cxx
STEER/AliAODEvent.h

index 01d0f6a..ddaede8 100644 (file)
 
 ClassImp(AliAODEvent)
 
+// definition of std AOD member names
+  const char* AliAODEvent::fAODListName[kAODListN] = {"header",
+                                                     "tracks",
+                                                     "vertices",
+                                                     "clusters",
+                                                     "jets",
+                                                     "tracklets"};
 //______________________________________________________________________________
 AliAODEvent::AliAODEvent() :
+  AliVEvent(),
   fAODObjects(new TList()),
   fHeader(0),
   fTracks(0),
@@ -59,7 +67,7 @@ void AliAODEvent::AddObject(TObject* obj)
 }
 
 //______________________________________________________________________________
-TObject *AliAODEvent::GetObject(const char *objName) const 
+TObject *AliAODEvent::FindListObject(const char *objName)
 {
   // Return the pointer to the object with the given name.
 
@@ -79,28 +87,56 @@ void AliAODEvent::CreateStdContent()
   AddObject(new TClonesArray("AliAODJet", 0));
   AddObject(new AliAODTracklets());
 
+  // set names
+  SetStdNames();
+
   // read back pointers
   GetStdContent();
 
-  // set names
-  fTracks->SetName("tracks");
-  fVertices->SetName("vertices");
-  fClusters->SetName("clusters");
-  fJets->SetName("jets");
-  fTracklets->SetName("tracklets");    
+  return;
 }
 
 //______________________________________________________________________________
+void AliAODEvent::SetStdNames()
+{
+  // introduce the standard naming
+
+  if(fAODObjects->GetEntries()==kAODListN){
+    for(int i = 0;i < fAODObjects->GetEntries();i++){
+      TObject *fObj = fAODObjects->At(i);
+      if(fObj->InheritsFrom("TNamed")){
+       ((TNamed*)fObj)->SetName(fAODListName[i]);
+      }
+      else if(fObj->InheritsFrom("TClonesArray")){
+       ((TClonesArray*)fObj)->SetName(fAODListName[i]);
+      }
+    }
+  }
+  else{
+    printf("%s:%d SetStdNames() Wrong number of Std Entries \n",(char*)__FILE__,__LINE__);
+  }
+} 
+
+//______________________________________________________________________________
 void AliAODEvent::GetStdContent()
 {
   // set pointers for standard content
 
+  fHeader    = (AliAODHeader*)fAODObjects->FindObject("header");
+  fTracks    = (TClonesArray*)fAODObjects->FindObject("tracks");
+  fVertices  = (TClonesArray*)fAODObjects->FindObject("vertices");
+  fClusters  = (TClonesArray*)fAODObjects->FindObject("clusters");
+  fJets      = (TClonesArray*)fAODObjects->FindObject("jets");
+  fTracklets = (AliAODTracklets*)fAODObjects->FindObject("tracklets");
+  
+  /* old implementation
   fHeader    = (AliAODHeader*)fAODObjects->At(0);
   fTracks    = (TClonesArray*)fAODObjects->At(1);
   fVertices  = (TClonesArray*)fAODObjects->At(2);
   fClusters  = (TClonesArray*)fAODObjects->At(3);
   fJets      = (TClonesArray*)fAODObjects->At(4);
   fTracklets = (AliAODTracklets*)fAODObjects->At(5);
+  */
 }
 
 //______________________________________________________________________________
@@ -135,7 +171,7 @@ Int_t AliAODEvent::GetMuonTracks(TRefArray *muonTracks) const
 
   AliAODTrack *track = 0;
   for (Int_t iTrack = 0; iTrack < GetNTracks(); iTrack++) {
-    if ((track = GetTrack(iTrack))->GetMostProbablePID() == AliAODTrack::kMuon) {
+    if ((track = GetTrack(iTrack))->IsMuonTrack()) {
       muonTracks->Add(track);
     }
   }
@@ -144,18 +180,27 @@ Int_t AliAODEvent::GetMuonTracks(TRefArray *muonTracks) const
 }
 
 
-
 void AliAODEvent::ReadFromTree(TTree *tree)
 {
-    // connects aod event to tree
-
-    fAODObjects = (TList*)((AliAODEvent*)tree->GetTree()->GetUserInfo()->FindObject("AliAODEvent"))->GetList(); 
-    TIter next(fAODObjects);
-    TNamed *el;
-    while((el=(TNamed*)next())){
-       TString bname(el->GetName());
-       tree->SetBranchAddress(bname.Data(),fAODObjects->GetObjectRef(el));
-    }
-    GetStdContent();
+  // connects aod event to tree
+  
+  // load the TTree
+  tree->LoadTree(0);
+
+  fAODObjects = (TList*)((AliAODEvent*)tree->GetTree()->GetUserInfo()->FindObject("AliAODEvent"))->GetList(); 
+  TIter next(fAODObjects);
+  TNamed *el;
+  while((el=(TNamed*)next())){
+    TString bname(el->GetName());
+    tree->SetBranchAddress(bname.Data(),fAODObjects->GetObjectRef(el));
+  }
+  GetStdContent();
 }
 
+//______________________________________________________________________________
+void AliAODEvent::Print(Option_t *) const
+{
+  // Something meaningful should be implemented here.
+  
+  return;
+}
index b52ab93..15d4eec 100644 (file)
@@ -16,6 +16,8 @@
 #include <TTree.h>
 #include <TNamed.h>
 
+#include "AliVEvent.h"
+#include "AliVParticle.h"
 #include "AliAODHeader.h"
 #include "AliAODTrack.h"
 #include "AliAODVertex.h"
 
 class TTree;
 
-class AliAODEvent : public TObject {
+class AliAODEvent : public AliVEvent {
 
  public :
+  enum AODListIndex_t {kAODHeader,
+                      kAODTracks,
+                      kAODVertices,
+                      kAODClusters,
+                      kAODJets,
+                      kAODTracklets,
+                      kAODListN
+  };
+
   AliAODEvent();
   virtual ~AliAODEvent();
 
@@ -35,62 +46,91 @@ class AliAODEvent : public TObject {
   //AliAODEvent& operator=(const AliAODEvent& aodevent);  // not implemented
 
   void          AddObject(TObject *obj);
-  TObject      *GetObject(const char *objName) const;
+  TObject      *FindListObject(const char *objName) ;
   TList        *GetList()                const { return fAODObjects; }
 
   // -- Header
-  AliAODHeader *GetHeader()              const { return (AliAODHeader*)fHeader; }
+  AliAODHeader *GetHeader()              const { return fHeader; }
   void          AddHeader(const AliAODHeader* hdx)
     {
        delete fHeader; fHeader = new AliAODHeader(*hdx);
        (fAODObjects->FirstLink())->SetObject(fHeader);
     }
 
+  // setters and getters for header information
+  void     SetRunNumber(Int_t n) {if (fHeader) fHeader->SetRunNumber(n);}
+  void     SetPeriodNumber(UInt_t n){if (fHeader) fHeader->SetPeriodNumber(n);}
+  void     SetOrbitNumber(UInt_t n) {if (fHeader) fHeader->SetOrbitNumber(n);}
+  void     SetBunchCrossNumber(UShort_t n) {if (fHeader) fHeader->SetBunchCrossNumber(n);}
+  void     SetMagneticField(Double_t mf){if (fHeader) fHeader->SetMagneticField(mf);}
+
+  Int_t    GetRunNumber() const {return fHeader ? fHeader->GetRunNumber() : -999;}
+  UInt_t   GetPeriodNumber() const {return fHeader ? fHeader->GetPeriodNumber() : 0;}
+  UInt_t   GetOrbitNumber() const {return fHeader ? fHeader->GetOrbitNumber() : 0;}
+  UShort_t GetBunchCrossNumber() const {return fHeader ? fHeader->GetBunchCrossNumber() : 0;}
+  Double_t GetMagneticField() const {return fHeader ? fHeader->GetMagneticField() : -999.;}
+  
+  void      SetEventType(UInt_t eventType){fHeader->SetEventType(eventType);}
+  void      SetTriggerMask(ULong64_t n) {fHeader->SetTriggerMask(n);}
+  void      SetTriggerCluster(UChar_t n) {fHeader->SetTriggerCluster(n);}
+
+  UInt_t    GetEventType()  const { return fHeader ? fHeader->GetEventType() : 0;}
+  ULong64_t GetTriggerMask() const {return fHeader ? fHeader->GetTriggerMask() : 0;}
+  UChar_t   GetTriggerCluster() const {return fHeader ? fHeader->GetTriggerCluster() : 0;}
+
+  Double_t  GetZDCN1Energy()        const { return fHeader ? fHeader->GetZDCN1Energy() : -999.; }
+  Double_t  GetZDCP1Energy()        const { return fHeader ? fHeader->GetZDCP1Energy() : -999.; }
+  Double_t  GetZDCN2Energy()        const { return fHeader ? fHeader->GetZDCN2Energy() : -999.; }
+  Double_t  GetZDCP2Energy()        const { return fHeader ? fHeader->GetZDCP2Energy() : -999.; }
+  Double_t  GetZDCEMEnergy()        const { return fHeader ? fHeader->GetZDCEMEnergy() : -999.; }
+
   // -- Tracks
   TClonesArray *GetTracks()              const { return fTracks; }
   Int_t         GetNTracks()             const { return fTracks->GetEntriesFast(); }
-  AliAODTrack  *GetTrack(Int_t nTrack)   const { return (AliAODTrack*)fTracks->At(nTrack); }
-  void          AddTrack(const AliAODTrack* trk)
-    {new((*fTracks)[fTracks->GetEntries()]) AliAODTrack(*trk);}
+  Int_t         GetNumberOfTracks()      const { return GetNTracks(); }
+  AliAODTrack  *GetTrack(Int_t nTrack)   const { return (AliAODTrack*)fTracks->UncheckedAt(nTrack); }
+  Int_t         AddTrack(const AliAODTrack* trk)
+    {new((*fTracks)[fTracks->GetEntriesFast()]) AliAODTrack(*trk); return fTracks->GetEntriesFast()-1;}
   Int_t         GetMuonTracks(TRefArray *muonTracks) const;
 
   // -- Vertex
   TClonesArray *GetVertices()            const { return fVertices; }
   Int_t         GetNVertices()           const { return fVertices->GetEntriesFast(); }
-  AliAODVertex *GetVertex(Int_t nVertex) const { return (AliAODVertex*)fVertices->At(nVertex); }
-  void          AddVertex(const AliAODVertex* vtx)
-    {new((*fVertices)[fVertices->GetEntries()]) AliAODVertex(*vtx);}
+  AliAODVertex *GetVertex(Int_t nVertex) const { return (AliAODVertex*)fVertices->UncheckedAt(nVertex); }
+  Int_t         AddVertex(const AliAODVertex* vtx)
+  {new((*fVertices)[fVertices->GetEntriesFast()]) AliAODVertex(*vtx); return fVertices->GetEntriesFast()-1;}
   virtual AliAODVertex *GetPrimaryVertex() const { return GetVertex(0); }
   
 
   // -- Cluster
   TClonesArray *GetClusters()            const { return fClusters; }
   Int_t         GetNClusters()           const { return fClusters->GetEntriesFast(); }
-  AliAODCluster *GetCluster(Int_t nCluster) const { return (AliAODCluster*)fClusters->At(nCluster); }
-  void          AddCluster(const AliAODCluster* vtx)
-    {new((*fClusters)[fClusters->GetEntries()]) AliAODCluster(*vtx);}
+  AliAODCluster *GetCluster(Int_t nCluster) const { return (AliAODCluster*)fClusters->UncheckedAt(nCluster); }
+  Int_t         AddCluster(const AliAODCluster* vtx)
+  {new((*fClusters)[fClusters->GetEntriesFast()]) AliAODCluster(*vtx); return fClusters->GetEntriesFast()-1;}
 
   // -- Jet
   TClonesArray *GetJets()            const { return fJets; }
   Int_t         GetNJets()           const { return fJets->GetEntriesFast(); }
-  AliAODJet *GetJet(Int_t nJet) const { return (AliAODJet*)fJets->At(nJet); }
-  void          AddJet(const AliAODJet* vtx)
-    {new((*fJets)[fJets->GetEntries()]) AliAODJet(*vtx);}
+  AliAODJet    *GetJet(Int_t nJet) const { return (AliAODJet*)fJets->UncheckedAt(nJet); }
+  Int_t         AddJet(const AliAODJet* vtx)
+    {new((*fJets)[fJets->GetEntriesFast()]) AliAODJet(*vtx); return fJets->GetEntriesFast()-1;}
 
   // -- Tracklets
   AliAODTracklets *GetTracklets() const { return fTracklets; }
 
   // -- Services
   void    CreateStdContent();
+  void    SetStdNames();
   void    GetStdContent();
   void    ResetStd(Int_t trkArrSize = 0, Int_t vtxArrSize = 0);
   void    ClearStd();
   void    ReadFromTree(TTree *tree);
-  void    WriteToTree(TTree* tree) {tree->Branch(fAODObjects);}
- private :
+  const void WriteToTree(TTree* tree) const {tree->Branch(fAODObjects);}
+
+  void  Print(Option_t *option="") const;
 
-  AliAODEvent(const AliAODEvent&); // Not implemented
-  AliAODEvent& operator=(const AliAODEvent&); // Not implemented
+ private :
 
   TList *fAODObjects; //  list of AODObjects
   
@@ -102,7 +142,9 @@ class AliAODEvent : public TObject {
   TClonesArray    *fJets;      //! jets
   AliAODTracklets *fTracklets; //! SPD tracklets
 
-  ClassDef(AliAODEvent,2);
+  static const char* fAODListName[kAODListN]; //!
+
+  ClassDef(AliAODEvent,3);
 };
 
 #endif