pileup rejection added M. Knichel

diff --git a/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysispPb.cxx b/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysispPb.cxx
index 0d81f59bb10419ff60d4fdd4338931e64b331b4c..2e3dbaee0f2dd40e398669ed3830a6e294453427 100644 (file)
--- a/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysispPb.cxx
+++ b/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysispPb.cxx
@@ -24,7 +24,7 @@
 // - generic correction matrices\r
 // - control histograms\r
 //\r
-// last change: 2013-06-13 by M.Knichel, J.Gronefeld\r
+// last change: 2013-06-19 by M.Knichel\r
 //\r
 // meaning of different multiplicities:\r
 // multRec      : number of reconstructed tracks, after AcceptanceCuts and TrackCuts\r
@@ -57,6 +57,7 @@
 #include "AliPhysicsSelection.h"\r
 #include "AliTriggerAnalysis.h"\r
 #include "AliCentrality.h"\r
+#include "AliAnalysisUtils.h"\r
 \r
 #include "AliPWG0Helper.h"\r
 #include "AlidNdPtHelper.h"\r
@@ -188,7 +189,9 @@ ClassImp(AlidNdPtAnalysispPb)
   fBinsZv(0),\r
   \r
   fRapidityShift(-4.65409416218532379e-01),\r
-\r
+  fUtils(0),\r
+  fIs2013pA(kTRUE),\r
+  \r
   fIsInit(kFALSE)  \r
   \r
 {\r
@@ -326,6 +329,8 @@ AlidNdPtAnalysispPb::AlidNdPtAnalysispPb(Char_t* name, Char_t* title): AlidNdPt(
   fBinsZv(0),\r
   \r
   fRapidityShift(-4.65409416218532379e-01),\r
+  fUtils(0),\r
+  fIs2013pA(kTRUE),\r
 \r
   fIsInit(kFALSE)  \r
   \r
@@ -445,6 +450,9 @@ AlidNdPtAnalysispPb::~AlidNdPtAnalysispPb() {
   if (fCentralityTrackHist) delete fCentralityTrackHist; fCentralityTrackHist=0;\r
   if (fVCentralityEvent) delete[] fVCentralityEvent; fVCentralityEvent=0;\r
   if (fVCentralityTrack) delete[] fVCentralityTrack; fVCentralityTrack=0;\r
+  \r
+  if (fUtils) delete fUtils; fUtils=0;\r
+  \r
 }\r
 \r
 \r
@@ -1246,6 +1254,8 @@ void AlidNdPtAnalysispPb::Init()
   // init folder\r
   fAnalysisFolder = CreateFolder("folderdNdPt","Analysis dNdPt Folder");\r
   \r
+  fUtils = new AliAnalysisUtils();\r
+  \r
   // set init flag\r
   fIsInit = kTRUE;\r
 }\r
@@ -1417,7 +1427,17 @@ void AlidNdPtAnalysispPb::Process(AliESDEvent *const esdEvent, AliMCEvent *const
 \r
   // the previous lines have been added for same vertex as dNdEta analysis\r
 \r
-  Bool_t isEventOK = evtCuts->AcceptEvent(esdEvent,mcEvent,vtxESD) && isRecVertex; \r
+  Bool_t isEventOK = evtCuts->AcceptEvent(esdEvent,mcEvent,vtxESD) && isRecVertex;\r
+  \r
+  Bool_t isEventSelected2013 = kTRUE;\r
+  \r
+  // selection and pileup rejection for 2013 p-A  \r
+  if (fIs2013pA) {\r
+       if (fUtils->IsFirstEventInChunk(esdEvent)) { isEventSelected2013 = kFALSE;  }\r
+       if (!fUtils->IsVertexSelected2013pA(esdEvent)) { isEventSelected2013 = kFALSE;  }\r
+       if (fUtils->IsPileUpEvent(esdEvent)) { isEventSelected2013 = kFALSE;  }\r
+  }\r
+  isEventOK = isEventOK && isEventSelected2013;\r
   //printf("isEventOK %d, isEventTriggered %d \n",isEventOK, isEventTriggered);\r
   //printf("GetAnalysisMode() %d \n",GetAnalysisMode());\r
   \r

diff --git a/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysispPb.h b/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysispPb.h
index 2b84be2726a1b80ea4bb4fe885a85d43a670c551..46f236e3284403614a2fea80b6ffb7969ddfab55 100644 (file)
--- a/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysispPb.h
+++ b/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysispPb.h
@@ -6,7 +6,7 @@
 // based on AlidNdPtAnalysis class
 // 
 // Author: J.Otwinowski 04/11/2008 
-// last change: 2013-06-13 by M.Knichel, J.Gronefeld
+// last change: 2013-06-19 by M.Knichel
 //------------------------------------------------------------------------------
 
 class iostream;
@@ -25,6 +25,7 @@ class AliESD;
 class AliESDfriend;
 class AliESDfriendTrack;
 class AlidNdPtHelper;
+class AliAnalysisUtils;
 
 #include "AlidNdPt.h"
 #include "TObjString.h"
@@ -183,6 +184,10 @@ public :
   //rapidity shift getter+setter
   void SetRapidityShift(Double_t yShift) { fRapidityShift = yShift;}
   Double_t GetRapidityShift() { return fRapidityShift; }
+  
+  void Set2013pA(Bool_t is2013 = kTRUE) { fIs2013pA = is2013; }
+  Double_t Get2013pA() { return fIs2013pA; }  
+
 
 private:
 
@@ -345,16 +350,18 @@ private:
   Double_t *fBinsZv;     //[fZvNedges]
   
   Double_t fRapidityShift; //y shift CMS vs. LAB
+  AliAnalysisUtils* fUtils;
+  Bool_t fIs2013pA; 
   
   Bool_t fIsInit;
-  
+    
   // generic function to change binning
   Bool_t CanChangeBins();
 
   AlidNdPtAnalysispPb(const AlidNdPtAnalysispPb&); // not implemented
   AlidNdPtAnalysispPb& operator=(const AlidNdPtAnalysispPb&); // not implemented
 
-  ClassDef(AlidNdPtAnalysispPb,5);
+  ClassDef(AlidNdPtAnalysispPb,6);
 };
 
 #endif