From 1be4eb33b917c5d0b528ed9add0b4571ab41d288 Mon Sep 17 00:00:00 2001
From: jgrosseo <jgrosseo@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Thu, 23 Jun 2011 13:14:27 +0000
Subject: [PATCH] updated track cut systematic study

---
 PWG4/JetTasks/AliAnalyseLeadingTrackUE.cxx | 29 +++++++++++++++++++---
 1 file changed, 25 insertions(+), 4 deletions(-)

diff --git a/PWG4/JetTasks/AliAnalyseLeadingTrackUE.cxx b/PWG4/JetTasks/AliAnalyseLeadingTrackUE.cxx
index f1ab13d7899..18d32a58e1b 100644
--- a/PWG4/JetTasks/AliAnalyseLeadingTrackUE.cxx
+++ b/PWG4/JetTasks/AliAnalyseLeadingTrackUE.cxx
@@ -112,16 +112,37 @@ Bool_t AliAnalyseLeadingTrackUE::ApplyCuts(TObject* track)
 
 
 //____________________________________________________________________
-void AliAnalyseLeadingTrackUE::DefineESDCuts(Int_t /*filterbit*/){
+void AliAnalyseLeadingTrackUE::DefineESDCuts(Int_t filterbit) {
   
   // Reproduces the cuts of the corresponding bit in the ESD->AOD filtering
   // (see $ALICE_ROOT/ANALYSIS/macros/AddTaskESDFilter.C)
+  
+  if (filterbit == -1)
+    filterbit = fFilterBit;
 
-  if (fFilterBit == 128)
+  if (filterbit == 128)
   {
     fEsdTrackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
     fEsdTrackCuts->SetMinNClustersTPC(70);
   }
+  else if (filterbit == 256)
+  {
+    // syst study
+    fEsdTrackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
+    fEsdTrackCuts->SetMinNClustersTPC(80);
+    fEsdTrackCuts->SetMaxChi2PerClusterTPC(3);
+    fEsdTrackCuts->SetMaxDCAToVertexZ(2.7);
+    fEsdTrackCuts->SetMaxDCAToVertexXY(1.9);
+  }
+  else if (filterbit == 512)
+  {
+    // syst study
+    fEsdTrackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
+    fEsdTrackCuts->SetMinNClustersTPC(60);
+    fEsdTrackCuts->SetMaxChi2PerClusterTPC(5);
+    fEsdTrackCuts->SetMaxDCAToVertexZ(3.7);
+    fEsdTrackCuts->SetMaxDCAToVertexXY(2.9);
+  }
   else
   {
     fEsdTrackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010();
@@ -181,7 +202,7 @@ TObjArray* AliAnalyseLeadingTrackUE::GetAcceptedParticles(TObject* obj, TObject*
   
   // for TPC only tracks
   Bool_t hasOwnership = kFALSE;
-  if (fFilterBit == 128 && obj->InheritsFrom("AliESDEvent"))
+  if ((fFilterBit == 128 || fFilterBit == 256 || fFilterBit == 512) && obj->InheritsFrom("AliESDEvent"))
     hasOwnership = kTRUE;
   
   if (!arrayMC)
@@ -451,7 +472,7 @@ AliVParticle*  AliAnalyseLeadingTrackUE::ParticleWithCuts(TObject* obj, Int_t ip
 	if (!( ApplyCuts(part)) )
 	 return 0; 
 	
-	if (fFilterBit == 128)
+	if (fFilterBit == 128 || fFilterBit == 256 || fFilterBit == 512)
 	{
 	  // create TPC only tracks constrained to the SPD vertex
 
-- 
2.43.5