From 08e4c457014aadde2f210e47d594b4e7738adffb Mon Sep 17 00:00:00 2001 From: hristov Date: Tue, 7 Aug 2007 20:54:04 +0000 Subject: [PATCH] Add inheritance from AliVEvent. Fix standard naming scheme. Fix to be able to ReadFromTree in case of TChains. Bug fixes concerning muon tracks (Markus) --- STEER/AliAODEvent.cxx | 83 +++++++++++++++++++++++++++++++++---------- STEER/AliAODEvent.h | 82 +++++++++++++++++++++++++++++++----------- 2 files changed, 126 insertions(+), 39 deletions(-) diff --git a/STEER/AliAODEvent.cxx b/STEER/AliAODEvent.cxx index 01d0f6ae731..ddaede88b03 100644 --- a/STEER/AliAODEvent.cxx +++ b/STEER/AliAODEvent.cxx @@ -28,8 +28,16 @@ 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; +} diff --git a/STEER/AliAODEvent.h b/STEER/AliAODEvent.h index b52ab93d425..15d4eec1726 100644 --- a/STEER/AliAODEvent.h +++ b/STEER/AliAODEvent.h @@ -16,6 +16,8 @@ #include #include +#include "AliVEvent.h" +#include "AliVParticle.h" #include "AliAODHeader.h" #include "AliAODTrack.h" #include "AliAODVertex.h" @@ -25,9 +27,18 @@ 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 -- 2.39.3