From 0b21f68696c8bdae9ab61a0e47cfb4049ed9c349 Mon Sep 17 00:00:00 2001 From: mcosenti Date: Mon, 6 Aug 2012 13:54:14 +0000 Subject: [PATCH] including another event counter --- PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx | 10 +++++++++- PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h | 1 + 2 files changed, 10 insertions(+), 1 deletion(-) diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx index cba26c6a6dc..b65d331a3a4 100644 --- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx +++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx @@ -54,6 +54,7 @@ AliAnalysisTaskEMCALIsoPhoton::AliAnalysisTaskEMCALIsoPhoton() : fESD(0), fOutputList(0), fEvtSel(0), + fNClusEt10(0), fPVtxZ(0), fCellAbsIdVsAmpl(0), fNClusHighClusE(0), @@ -82,6 +83,7 @@ AliAnalysisTaskEMCALIsoPhoton::AliAnalysisTaskEMCALIsoPhoton(const char *name) : fESD(0), fOutputList(0), fEvtSel(0), + fNClusEt10(0), fPVtxZ(0), fCellAbsIdVsAmpl(0), fNClusHighClusE(0), @@ -113,7 +115,9 @@ void AliAnalysisTaskEMCALIsoPhoton::UserCreateOutputObjects() fEvtSel = new TH1F("hEvtSel","Event selection counter (0=all trg, 1=pvz cut) ;evt cut ;dN/dcut}",2,0,2); fOutputList->Add(fEvtSel); - + + fNClusEt10 = new TH1F("hNClusEt10","# of cluster with E_{T}>10 per event;E;",101,-0.5,100.5); + fOutputList->Add(fNClusEt10); fPVtxZ = new TH1F("hPVtxZ","primary vertex Z before cut;prim-vz(cm) ;",200,-100,100); fOutputList->Add(fPVtxZ); @@ -222,6 +226,7 @@ void AliAnalysisTaskEMCALIsoPhoton::FillClusHists() if(nclus==0) return; Double_t maxE = 0; + Int_t nclus10 = 0; for(Int_t ic=0;ic(fCaloClusters->At(ic)); @@ -240,6 +245,8 @@ void AliAnalysisTaskEMCALIsoPhoton::FillClusHists() c->GetPosition(clsPos); TVector3 clsVec(clsPos); Double_t Et = c->E()*TMath::Sin(clsVec.Theta()); + if(Et>10) + nclus10++; Float_t ceiso, cephiband, cecore; Float_t triso, trphiband, trcore; Float_t alliso, allphiband, allcore; @@ -271,6 +278,7 @@ void AliAnalysisTaskEMCALIsoPhoton::FillClusHists() maxE = c->E(); } fNClusHighClusE->Fill(maxE,nclus); + fNClusEt10->Fill(nclus10); } //________________________________________________________________________ diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h index 4ba112f9843..1dd20b59a7f 100644 --- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h +++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h @@ -61,6 +61,7 @@ class AliAnalysisTaskEMCALIsoPhoton : public AliAnalysisTaskSE { TList *fOutputList; //! Output list //histograms for events with 1+ track pt>1 TH1F *fEvtSel; //!evt selection counter: 0=all trg, 1=pv cut + TH1F *fNClusEt10; //!number of clusters w/ Et>10 in the event TH1F *fPVtxZ; //!primary vertex Z before cut TH2F *fCellAbsIdVsAmpl; //!cell abs id vs cell amplitude (energy) TH2F *fNClusHighClusE; //!total number of clusters vs. highest clus energy in the event -- 2.43.0