From 099de61e5d5f172f25f336ac23a844953084b0c1 Mon Sep 17 00:00:00 2001 From: gconesab Date: Mon, 8 Oct 2012 13:31:19 +0000 Subject: [PATCH] add pile-up event counter based on SPD --- .../AliAnaCaloTrackCorrMaker.cxx | 24 ++++++++++++++----- .../AliAnaCaloTrackCorrMaker.h | 1 + PWG/CaloTrackCorrBase/AliCaloTrackReader.cxx | 4 ++++ 3 files changed, 23 insertions(+), 6 deletions(-) diff --git a/PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.cxx b/PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.cxx index c2a6217a5f9..1d61a8db8ee 100755 --- a/PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.cxx +++ b/PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.cxx @@ -47,7 +47,8 @@ fOutputContainer(new TList ), fAnalysisContainer(new TList ), fMakeHisto(kFALSE), fMakeAOD(kFALSE), fAnaDebug(0), fCuts(new TList), fScaleFactor(-1), -fhNEvents(0), fhZVertex(0), +fhNEvents(0), fhNPileUpEvents(0), +fhZVertex(0), fhTrackMult(0), fhCentrality(0), fhEventPlaneAngle(0), fhNMergedFiles(0), fhScaleFactor(0) @@ -68,8 +69,11 @@ fOutputContainer(new TList()), fAnalysisContainer(new TList()), fMakeHisto(maker.fMakeHisto), fMakeAOD(maker.fMakeAOD), fAnaDebug(maker.fAnaDebug), fCuts(new TList()), fScaleFactor(maker.fScaleFactor), -fhNEvents(maker.fhNEvents), fhZVertex(maker.fhZVertex), -fhTrackMult(maker.fhTrackMult),fhCentrality(maker.fhCentrality), +fhNEvents(maker.fhNEvents), +fhNPileUpEvents(maker.fhNPileUpEvents), +fhZVertex(maker.fhZVertex), +fhTrackMult(maker.fhTrackMult), +fhCentrality(maker.fhCentrality), fhEventPlaneAngle(maker.fhEventPlaneAngle), fhNMergedFiles(maker.fhNMergedFiles), fhScaleFactor(maker.fhScaleFactor) @@ -175,6 +179,12 @@ TList *AliAnaCaloTrackCorrMaker::GetOutputContainer() fhNEvents->SetYTitle("# events"); fOutputContainer->Add(fhNEvents); + fhNPileUpEvents = new TH1I("hNPileUpEvents", "Number of events considered as pile-up", 2 , 0 , 2 ) ; + fhNPileUpEvents->SetYTitle("# events"); + fhNPileUpEvents->GetXaxis()->SetBinLabel(1 ,"SPD"); + fhNPileUpEvents->GetXaxis()->SetBinLabel(2 ,"Multi SPD"); + fOutputContainer->Add(fhNPileUpEvents); + fhZVertex = new TH1F("hZVertex", " Z vertex distribution" , 200 , -50 , 50 ) ; fhZVertex->SetXTitle("v_{z} (cm)"); fOutputContainer->Add(fhZVertex); @@ -414,16 +424,18 @@ void AliAnaCaloTrackCorrMaker::ProcessEvent(const Int_t iEntry, // Event control histograms fhNEvents ->Fill(0); // Number of events analyzed + if(fReader->IsPileUpFromSPD()) + fhNPileUpEvents->Fill(0.5); + if(fReader->GetInputEvent()->IsPileupFromSPDInMultBins()) + fhNPileUpEvents->Fill(1.5); fhTrackMult ->Fill(fReader->GetTrackMultiplicity()); fhCentrality ->Fill(fReader->GetEventCentrality ()); fhEventPlaneAngle->Fill(fReader->GetEventPlaneAngle ()); - - + Double_t v[3]; fReader->GetInputEvent()->GetPrimaryVertex()->GetXYZ(v) ; fhZVertex->Fill(v[2]); - //printf(">>>>>>>>>> AFTER >>>>>>>>>>>\n"); //gObjectTable->Print(); diff --git a/PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.h b/PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.h index 8670dd5c2d4..e5746fc9128 100755 --- a/PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.h +++ b/PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.h @@ -92,6 +92,7 @@ class AliAnaCaloTrackCorrMaker : public TObject { // Control histograms TH1I * fhNEvents; //! Number of events counter histogram + TH1I * fhNPileUpEvents; //! N events pasing pile up cut TH1F * fhZVertex; //! Vertex of accepted event TH1I * fhTrackMult; //! Number of tracks per event histogram TH1F * fhCentrality; //! Histogram with centrality bins diff --git a/PWG/CaloTrackCorrBase/AliCaloTrackReader.cxx b/PWG/CaloTrackCorrBase/AliCaloTrackReader.cxx index 02b0b9257b1..f2207ce2d96 100755 --- a/PWG/CaloTrackCorrBase/AliCaloTrackReader.cxx +++ b/PWG/CaloTrackCorrBase/AliCaloTrackReader.cxx @@ -446,6 +446,8 @@ Bool_t AliCaloTrackReader::IsPileUpFromSPD() const // Default values: (3, 0.8, 3., 2., 5.) return fInputEvent->IsPileupFromSPD((Int_t) fPileUpParam[0] , fPileUpParam[1] , fPileUpParam[2] , fPileUpParam[3] , fPileUpParam[4] ); + //printf("Param : %d, %2.2f, %2.2f, %2.2f, %2.2f\n",(Int_t) fPileUpParam[0], fPileUpParam[1], fPileUpParam[2], fPileUpParam[3], fPileUpParam[4]); + } //________________________________________________________ @@ -610,6 +612,8 @@ Bool_t AliCaloTrackReader::FillInputEvent(const Int_t iEntry, TMath::Abs(fVertex[0][2] ) < 1.e-6 ) return kFALSE; } + //printf("Reader : IsPileUp %d, Multi %d\n",IsPileUpFromSPD(),fInputEvent->IsPileupFromSPDInMultBins()); + if(fDoEventSelection) { if(!fCaloFilterPatch) -- 2.43.0