implement integrated flow in scalar product method and subevents in flowevent

[u/mrichter/AliRoot.git] / PWG0 / CreateStandardCuts.C

diff --git a/PWG0/CreateStandardCuts.C b/PWG0/CreateStandardCuts.C
index c801069d015c9b99292f623bef47a21094c05104..fb6c760593298992566e416d342471b328429d76 100644 (file)
--- a/PWG0/CreateStandardCuts.C
+++ b/PWG0/CreateStandardCuts.C
@@ -15,17 +15,24 @@ AliESDtrackCuts* CreateTrackCuts(AliPWG0Helper::AnalysisMode analysisMode, Bool_
   Double_t cov3 = 0.5;
   Double_t cov4 = 0.5;
   Double_t cov5 = 2;
-  Double_t nSigma = 3;
+  Double_t nSigma = 4;
+
+  Bool_t tpcRefit = kTRUE;
+  Bool_t sigmaToVertex = kTRUE;
 
   TString tag("Global tracking");
 
   // TPC-only cuts
-  if (analysisMode == AliPWG0Helper::kTPC) 
+  if (analysisMode == AliPWG0Helper::kTPC)
   {
-    // beta cuts (still under investigation)
-    cov1 = 4;
-    cov2 = 4;
-    nSigma = 4;
+    cov1 = 9;
+    cov2 = 9;
+    cov3 = 1e10;
+    cov4 = 1e10;
+    cov5 = 1e10;
+    
+    tpcRefit = kFALSE;
+    sigmaToVertex = kFALSE;
     
     tag = "TPC-only tracking";
   }
@@ -36,14 +43,23 @@ AliESDtrackCuts* CreateTrackCuts(AliPWG0Helper::AnalysisMode analysisMode, Bool_
     cov5 = 1e10;
     tag += " without field";
   }
-  
+
   esdTrackCuts->SetMaxCovDiagonalElements(cov1, cov2, cov3, cov4, cov5);
 
-  esdTrackCuts->SetMinNsigmaToVertex(nSigma);
-  esdTrackCuts->SetRequireSigmaToVertex(kTRUE);
+  esdTrackCuts->SetRequireSigmaToVertex(sigmaToVertex);
+
+  if (sigmaToVertex) {
+    esdTrackCuts->SetMaxNsigmaToVertex(nSigma);
+  }
+  else{
+    
+    esdTrackCuts->SetMaxDCAToVertexZ(3.2);
+    esdTrackCuts->SetMaxDCAToVertexXY(2.4);
+    esdTrackCuts->SetDCAToVertex2D(kTRUE);
+  }
 
-  esdTrackCuts->SetRequireTPCRefit(kTRUE);
-  esdTrackCuts->SetAcceptKingDaughters(kFALSE);
+  esdTrackCuts->SetRequireTPCRefit(tpcRefit);
+  esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
 
   esdTrackCuts->SetMinNClustersTPC(50);
   esdTrackCuts->SetMaxChi2PerClusterTPC(3.5);