ZDC centrality selection (Jochen and Michael)

diff --git a/PWG0/dNdPt/AlidNdPtAnalysisPbPb.cxx b/PWG0/dNdPt/AlidNdPtAnalysisPbPb.cxx
index 682453c1b2f6f45e9a49f8fa2a8dd3d1f47a49f1..e25ec1f52140f48c978b92f8e54c6dda60767722 100644 (file)
--- a/PWG0/dNdPt/AlidNdPtAnalysisPbPb.cxx
+++ b/PWG0/dNdPt/AlidNdPtAnalysisPbPb.cxx
@@ -104,7 +104,9 @@ ClassImp(AlidNdPtAnalysisPbPb)
   fMCMultRecTrackHist1(0), 
 
   // rec. track control histograms
-  fRecTrackHist3(0)
+  fRecTrackHist3(0),
+  
+  fTriggerAnalysis(0)
 {
   // default constructor
   for(Int_t i=0; i<AlidNdPtHelper::kCutSteps; i++) { 
@@ -170,7 +172,9 @@ AlidNdPtAnalysisPbPb::AlidNdPtAnalysisPbPb(Char_t* name, Char_t* title): AlidNdP
   fMCMultRecTrackHist1(0), 
 
   // rec. track control histograms
-  fRecTrackHist3(0)
+  fRecTrackHist3(0),
+  
+  fTriggerAnalysis(0)
 {
   //
   // constructor
@@ -237,6 +241,8 @@ AlidNdPtAnalysisPbPb::~AlidNdPtAnalysisPbPb() {
   if(fRecTrackHist3) delete fRecTrackHist3; fRecTrackHist3=0; 
   //
   if(fAnalysisFolder) delete fAnalysisFolder; fAnalysisFolder=0;
+  
+  if (fTriggerAnalysis) delete fTriggerAnalysis;  fTriggerAnalysis = 0;
 }
 
 //_____________________________________________________________________________
@@ -616,6 +622,9 @@ void AlidNdPtAnalysisPbPb::Init(){
 
   // init folder
   fAnalysisFolder = CreateFolder("folderdNdPt","Analysis dNdPt Folder");
+  
+  // init trigger analysis (for zdc cut)
+  fTriggerAnalysis = new AliTriggerAnalysis;
 }
 
 //_____________________________________________________________________________
@@ -628,7 +637,23 @@ void AlidNdPtAnalysisPbPb::Process(AliESDEvent *const esdEvent, AliMCEvent *cons
     AliDebug(AliLog::kError, "esdEvent not available");
     return;
   }
+  
+   // zdc cut
+   if (!fTriggerAnalysis->ZDCTimeTrigger(esdEvent)) {
+     return; 
+   }
 
+  // track cuts from Jochen
+  const AliESDVertex* vtxESDTPC = esdEvent->GetPrimaryVertexTPC();
+  if( vtxESDTPC->GetNContributors() < 1 ) {
+    return;
+  }
+
+   // francesco prino cut
+  const AliMultiplicity* multESD = esdEvent->GetMultiplicity();
+  if( vtxESDTPC->GetNContributors() < (-10.+0.25*multESD->GetNumberOfITSClusters(0)) ) {
+    return;
+  } 
 
   // get selection cuts
   AlidNdPtEventCuts *evtCuts = GetEventCuts(); 

diff --git a/PWG0/dNdPt/AlidNdPtAnalysisPbPb.h b/PWG0/dNdPt/AlidNdPtAnalysisPbPb.h
index 15297043f3b3264b5d4f1786f6448ad94e7ab670..583540f6d6a5ae645e261976c016407a5735d28f 100644 (file)
--- a/PWG0/dNdPt/AlidNdPtAnalysisPbPb.h
+++ b/PWG0/dNdPt/AlidNdPtAnalysisPbPb.h
@@ -25,6 +25,7 @@ class AliESD;
 class AliESDfriend;
 class AliESDfriendTrack;
 class AlidNdPtHelper;
+class AliTriggerAnalysis;
 
 #include "AlidNdPt.h"
 
@@ -181,6 +182,8 @@ private:
   // track control histograms
   THnSparseF *fRecTrackHist3;  //-> nclust:chi2:Pt:Eta:Phi
 
+  AliTriggerAnalysis *fTriggerAnalysis; //! trigger analysis object;
+
   AlidNdPtAnalysisPbPb(const AlidNdPtAnalysisPbPb&); // not implemented
   AlidNdPtAnalysisPbPb& operator=(const AlidNdPtAnalysisPbPb&); // not implemented