pile up cuts (Leonardo)

diff --git a/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx b/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx
index 85c7ab266516591791baf013e7ea632f4a06dfcf..b186797ae4f33f9361e758ec44f1c7f5122910fd 100644 (file)
--- a/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx
+++ b/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx
@@ -55,6 +55,7 @@
 #include "AliESDtrackCuts.h"
 
 #include "AliHelperPID.h"
+#include "AliAnalysisUtils.h"
 
 ////////////////////////////////////////////////////////////////////////
 //
@@ -97,6 +98,7 @@ fSkipTrigger(kFALSE),
 fInjectedSignals(kFALSE),
 // pointers to UE classes
 fHelperPID(0x0),
+fAnalysisUtils(0x0),
 fAnalyseUE(0x0),
 fHistos(0x0),
 fHistosMixed(0),
@@ -867,7 +869,7 @@ void  AliAnalysisTaskPhiCorrelations::AnalyseDataMode()
   // skip fast cluster events here if requested
   if (fSkipFastCluster && (fInputHandler->IsEventSelected() & AliVEvent::kFastOnly))
     return;
-
+ 
   // Support for ESD and AOD based analysis
   AliVEvent* inputEvent = fAOD;
   if (!inputEvent)
@@ -958,10 +960,19 @@ void  AliAnalysisTaskPhiCorrelations::AnalyseDataMode()
   
   // Trigger selection ************************************************
   if (!fSkipTrigger && !fAnalyseUE->TriggerSelection(fInputHandler)) return;
-  
+
   // Fill the "event-counting-container", it is needed to get the number of events remaining after each event-selection cut
   fHistos->FillEvent(centrality, AliUEHist::kCFStepTriggered);
   
+  // Pileup selection ************************************************
+  if (fAnalysisUtils && fAnalysisUtils->IsPileUpEvent(inputEvent)) 
+  {
+    // count the removed events
+    fHistos->FillEvent(centrality, AliUEHist::kCFStepAnaTopology);
+
+    return;
+  }
+  
   // Vertex selection *************************************************
   if(!fAnalyseUE->VertexSelection(inputEvent, fnTracksVertex, fZVertex)) return;
   

diff --git a/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.h b/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.h
index 9deac802d56ff1ee8487a7c389817368c00787e1..457407d7ca99024478af78008e24902501abbf7a 100644 (file)
--- a/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.h
+++ b/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.h
@@ -45,6 +45,7 @@ class TObjArray;
 class AliEventPoolManager;
 class AliESDEvent;
 class AliHelperPID;
+class AliAnalysisUtils;
 
 
 class  AliAnalysisTaskPhiCorrelations : public AliAnalysisTask
@@ -111,6 +112,9 @@ class  AliAnalysisTaskPhiCorrelations : public AliAnalysisTask
     
     AliHelperPID* GetHelperPID() { return fHelperPID; }
     void   SetHelperPID(AliHelperPID* pid){ fHelperPID = pid; }
+   
+    AliAnalysisUtils* GetAnalysisUtils() { return fAnalysisUtils; }
+    void   SetAnalysisUtils(AliAnalysisUtils* utils){ fAnalysisUtils = utils; }
     
   private:
     AliAnalysisTaskPhiCorrelations(const  AliAnalysisTaskPhiCorrelations &det);
@@ -141,6 +145,7 @@ class  AliAnalysisTaskPhiCorrelations : public AliAnalysisTask
     
     // Pointers to external UE classes
     AliHelperPID*     fHelperPID;      // points to class containing common analysis algorithms
+    AliAnalysisUtils*     fAnalysisUtils;      // points to class containing common analysis algorithms
     AliAnalyseLeadingTrackUE*     fAnalyseUE;      //! points to class containing common analysis algorithms
     AliUEHistograms*  fHistos;       //! points to class to handle histograms/containers  
     AliUEHistograms*  fHistosMixed;       //! points to class to handle mixed histograms/containers  
@@ -196,7 +201,7 @@ class  AliAnalysisTaskPhiCorrelations : public AliAnalysisTask
     
     Bool_t fFillpT;                // fill sum pT instead of number density
     
-    ClassDef( AliAnalysisTaskPhiCorrelations, 31); // Analysis task for delta phi correlations
+    ClassDef( AliAnalysisTaskPhiCorrelations, 32); // Analysis task for delta phi correlations
   };
 
 class AliDPhiBasicParticle : public AliVParticle