From 3f4073ba02599990b0d83b1073c03adeef88e79f Mon Sep 17 00:00:00 2001 From: mcosenti Date: Thu, 12 Jul 2012 14:09:08 +0000 Subject: [PATCH] updato on how to read the tracks --- PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx | 13 ++++++++++--- PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h | 1 + 2 files changed, 11 insertions(+), 3 deletions(-) diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx index f9ec72619df..99d4a0409fb 100644 --- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx +++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx @@ -39,6 +39,7 @@ AliAnalysisTaskEMCALIsoPhoton::AliAnalysisTaskEMCALIsoPhoton() : AliAnalysisTaskSE(), fCaloClusters(0), fSelPrimTracks(0), + fTracks(0), fEMCalCells(0), fPrTrCuts(0), fGeom(0x0), @@ -66,6 +67,7 @@ AliAnalysisTaskEMCALIsoPhoton::AliAnalysisTaskEMCALIsoPhoton(const char *name) : AliAnalysisTaskSE(name), fCaloClusters(0), fSelPrimTracks(0), + fTracks(0), fEMCalCells(0), fPrTrCuts(0), fGeom(0x0), @@ -174,9 +176,14 @@ void AliAnalysisTaskEMCALIsoPhoton::UserExec(Option_t *) fEvtSel->Fill(1); + if (!fTracks) + fTracks = dynamic_cast(InputEvent()->FindListObject("Tracks")); // Track loop to fill a pT spectrum - for (Int_t iTracks = 0; iTracks < fESD->GetNumberOfTracks(); iTracks++) { - AliESDtrack* track = (AliESDtrack*)fESD->GetTrack(iTracks); + const Int_t Ntracks = fTracks->GetEntriesFast(); + for (Int_t iTracks = 0; iTracks < Ntracks; ++iTracks) { + // for (Int_t iTracks = 0; iTracks < fESD->GetNumberOfTracks(); iTracks++) { + //AliESDtrack* track = (AliESDtrack*)fESD->GetTrack(iTracks); + AliVTrack *track = static_cast(fTracks->At(iTracks)); if (!track) continue; if (fPrTrCuts && fPrTrCuts->IsSelected(track)){ @@ -322,7 +329,7 @@ void AliAnalysisTaskEMCALIsoPhoton::GetTrIso(TVector3 vec, Float_t &iso, Float_t if(phicl<0) phicl+=TMath::TwoPi(); for(int itrack=0;itrack (fSelPrimTracks->At(itrack)); + AliVTrack *track = static_cast (fSelPrimTracks->At(itrack)); if(!track) continue; Double_t dphi = TMath::Abs(track->Phi()-phicl); diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h index f9f6620ed3d..37237f51cb1 100644 --- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h +++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h @@ -43,6 +43,7 @@ class AliAnalysisTaskEMCALIsoPhoton : public AliAnalysisTaskSE { protected: TRefArray *fCaloClusters; //!pointer to EMCal clusters TObjArray *fSelPrimTracks; //!pointer to ESD primary tracks + TClonesArray *fTracks; //!track input array AliESDCaloCells *fEMCalCells; //!pointer to EMCal cells AliESDtrackCuts *fPrTrCuts; //!pointer to hold the prim track cuts AliEMCALGeometry *fGeom; // geometry utils -- 2.43.0