From 6300a3af6a104f7675722dab8d9b322b6119e306 Mon Sep 17 00:00:00 2001 From: mcosenti Date: Fri, 12 Sep 2014 10:37:34 -0300 Subject: [PATCH] including a setter to turn on pile-up rejection with one of the SPD methods (not the multiplicity dependent one) --- PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx | 7 +++++++ PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h | 2 ++ 2 files changed, 9 insertions(+) diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx index db5e57734a4..6fc8afe7170 100644 --- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx +++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx @@ -90,6 +90,7 @@ AliAnalysisTaskEMCALIsoPhoton::AliAnalysisTaskEMCALIsoPhoton() : fNCuts(5), fTrCoreRem(kFALSE), fClusTDiff(30e-9), + fPileUpRejSPD(kFALSE), fESD(0), fAOD(0), fVEvent(0), @@ -197,6 +198,7 @@ AliAnalysisTaskEMCALIsoPhoton::AliAnalysisTaskEMCALIsoPhoton(const char *name) : fNCuts(5), fTrCoreRem(kFALSE), fClusTDiff(30e-9), + fPileUpRejSPD(kFALSE), fESD(0), fAOD(0), fVEvent(0), @@ -541,6 +543,11 @@ void AliAnalysisTaskEMCALIsoPhoton::UserExec(Option_t *) printf("passed vertex cut\n"); fEvtSel->Fill(1); + if(fVEvent->IsPileupFromSPD(3, 0.8, 3., 2., 5.) && fPileUpRejSPD){ + if(fDebug) + printf("Event is SPD pile-up!***\n"); + return; + } if(fESD) fTracks = dynamic_cast(InputEvent()->FindListObject("Tracks")); if(fAOD) diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h index 9605a7492e1..4df15e4e1c4 100644 --- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h +++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h @@ -75,6 +75,7 @@ class AliAnalysisTaskEMCALIsoPhoton : public AliAnalysisTaskSE { void SetMinIsoClusE(Double_t emin) { fMinIsoClusE = emin; } void SetTrCoreRemoval(Bool_t b) { fTrCoreRem = b; } void SetClusTDiff(Double_t diff) { fClusTDiff = diff; } + void SetPileUpRejSPD() { fPileUpRejSPD = kTRUE; } protected: TObjArray *fESDClusters; //!pointer to EMCal clusters TObjArray *fAODClusters; //!pointer to EMCal clusters @@ -121,6 +122,7 @@ class AliAnalysisTaskEMCALIsoPhoton : public AliAnalysisTaskSE { Int_t fNCuts; // number of cuts (QA purposes) Bool_t fTrCoreRem; // flag to set the removal of the core in track isolation (true removes it, default) Double_t fClusTDiff; // variable to hold the time diff between the candidate cluster and the isolation clusters + Bool_t fPileUpRejSPD; // flag to set pile-up rejection via SPD (multiple vertices) private: AliESDEvent *fESD; //! ESD object -- 2.39.3