new standard cuts

diff --git a/PWG4/macros/AddTaskESDFilterPWG4Train.C b/PWG4/macros/AddTaskESDFilterPWG4Train.C
index abcbee6012d10ff97f3f6e01ce934ad641ab5003..37dec7545107a7835bfdea219dbea8e0f6079fa6 100644 (file)
--- a/PWG4/macros/AddTaskESDFilterPWG4Train.C
+++ b/PWG4/macros/AddTaskESDFilterPWG4Train.C
@@ -97,7 +97,7 @@ AliAnalysisTaskESDfilter *AddTaskESDFilterPWG4Train(Bool_t useKineFilter=kTRUE,
 
    // ITS cuts for new jet analysis 
    gROOT->LoadMacro("$ALICE_ROOT/PWG4/macros/CreateTrackCutsPWG4.C");
-   AliESDtrackCuts* esdTrackCutsHG0 = CreateTrackCutsPWG4(10001005);
+   AliESDtrackCuts* esdTrackCutsHG0 = CreateTrackCutsPWG4(10001006);
 
    // throw out tracks with too low number of clusters in
    // the first pass (be consistent with TPC only tracks)
@@ -109,11 +109,11 @@ AliAnalysisTaskESDfilter *AddTaskESDFilterPWG4Train(Bool_t useKineFilter=kTRUE,
 
 
    // the complement to the one with SPD requirement
-   AliESDtrackCuts* esdTrackCutsHG1 = CreateTrackCutsPWG4(10011005);
+   AliESDtrackCuts* esdTrackCutsHG1 = CreateTrackCutsPWG4(10011006);
 
    // the tracks that must not be taken pass this cut and
    // non HGC1 and HG
-   AliESDtrackCuts* esdTrackCutsHG2 = CreateTrackCutsPWG4(10021005);
+   AliESDtrackCuts* esdTrackCutsHG2 = CreateTrackCutsPWG4(10021006);
 
    
 
@@ -124,7 +124,7 @@ AliAnalysisTaskESDfilter *AddTaskESDFilterPWG4Train(Bool_t useKineFilter=kTRUE,
    esdTrackCutsH2->SetMaxChi2PerClusterITS(36.);
    esdTrackCutsH2->SetPtRange(0.15,1E10);
 
-   AliESDtrackCuts* esdTrackCutsGCOnly = CreateTrackCutsPWG4(10041005);
+   AliESDtrackCuts* esdTrackCutsGCOnly = CreateTrackCutsPWG4(10041006);
 
    // TPC only tracks
    AliESDtrackCuts* esdTrackCutsTPCCOnly = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();

diff --git a/PWG4/macros/AddTaskPWG4HighPtQAMC.C b/PWG4/macros/AddTaskPWG4HighPtQAMC.C
index 7bf377eada084fe1024a1e2d855596e6b4fcb316..f82b038bf65d239d3dd4c7c22d100fa17b9fdd6b 100644 (file)
--- a/PWG4/macros/AddTaskPWG4HighPtQAMC.C
+++ b/PWG4/macros/AddTaskPWG4HighPtQAMC.C
@@ -56,27 +56,27 @@ AliPWG4HighPtQAMC* AddTaskPWG4HighPtQAMC(char *prodType = "LHC10e14", Int_t trac
   }
   if(trackType==0 && cuts==1) {
     //Cuts global tracks with ITSrefit requirement and SPDrequirement for jet analysis
-    trackCuts = CreateTrackCutsPWG4(10001005);
+    trackCuts = CreateTrackCutsPWG4(10001006);
    }
   if(trackType==0 && cuts==2) {
     //Cuts global tracks with ITSrefit requirement but without SPD
-    trackCuts = CreateTrackCutsPWG4(10011005);
+    trackCuts = CreateTrackCutsPWG4(10011006);
   }
   if(trackType==7 && cuts==0) {
     // tight global tracks
-    trackCuts = CreateTrackCutsPWG4(10041005);
-    trackCutsReject = CreateTrackCutsPWG4(1005);
+    trackCuts = CreateTrackCutsPWG4(10041006);
+    trackCutsReject = CreateTrackCutsPWG4(1006);
     trackCutsReject->SetEtaRange(-0.9,0.9);
     trackCutsReject->SetPtRange(0.15, 1e10);
   }
   if(trackType==7 && cuts==1) {
     // tight global tracks
-    trackCuts = CreateTrackCutsPWG4(10011005);
+    trackCuts = CreateTrackCutsPWG4(10011006);
   }
   if(trackType==7 && cuts==2) {
     // no requirements on SPD and ITSrefit failed
-    trackCuts = CreateTrackCutsPWG4(10041005);       //no ITSrefit requirement filter 256
-    trackCutsReject = CreateTrackCutsPWG4(10001005); //ITSrefit requirement filter 16
+    trackCuts = CreateTrackCutsPWG4(10041006);       //no ITSrefit requirement filter 256
+    trackCutsReject = CreateTrackCutsPWG4(10001006); //ITSrefit requirement filter 16
     trackCutsReject->SetEtaRange(-0.9,0.9);
     trackCutsReject->SetPtRange(0.15, 1e10);
   }

diff --git a/PWG4/macros/AddTaskPWG4HighPtSpectra.C b/PWG4/macros/AddTaskPWG4HighPtSpectra.C
index 5e231372b54b5c39b66b1ba9c0fe5c7cb69fa5b3..57552abf0b6b39c724da48acbe42b8475d41f082 100644 (file)
--- a/PWG4/macros/AddTaskPWG4HighPtSpectra.C
+++ b/PWG4/macros/AddTaskPWG4HighPtSpectra.C
@@ -177,27 +177,27 @@ AliPWG4HighPtSpectra* AddTaskPWG4HighPtSpectra(char *prodType = "LHC10e14", Bool
   }
   if(trackType==0 && cuts==1) {
     //Cuts global tracks with ITSrefit requirement and SPDrequirement for jet analysis
-    trackCuts = CreateTrackCutsPWG4(10001005);
+    trackCuts = CreateTrackCutsPWG4(10001006);
   }
   if(trackType==0 && cuts==2) {
     //Cuts global tracks with ITSrefit requirement but without SPD
-    trackCuts = CreateTrackCutsPWG4(10011005);
+    trackCuts = CreateTrackCutsPWG4(10011006);
   }
   if(trackType==7 && cuts==0) {
     // no requirements on SPD and ITSrefit failed
-    trackCuts = CreateTrackCutsPWG4(10041005);   //no ITSrefit requirement
-    trackCutsReject = CreateTrackCutsPWG4(1005); //ITSrefit requirement
+    trackCuts = CreateTrackCutsPWG4(10041006);   //no ITSrefit requirement
+    trackCutsReject = CreateTrackCutsPWG4(1006); //ITSrefit requirement
     trackCutsReject->SetEtaRange(etamin,etamax);
     trackCutsReject->SetPtRange(0.15, 1e10);
   }
   if(trackType==7 && cuts==1) {
     //Cuts global tracks with ITSrefit requirement but without SPD
-    trackCuts = CreateTrackCutsPWG4(10011005);
+    trackCuts = CreateTrackCutsPWG4(10011006);
   }
   if(trackType==7 && cuts==2) {
     // no requirements on SPD and ITSrefit failed
-    trackCuts = CreateTrackCutsPWG4(10041005);       //no ITSrefit requirement filter 256
-    trackCutsReject = CreateTrackCutsPWG4(10001005); //ITSrefit requirement filter 16
+    trackCuts = CreateTrackCutsPWG4(10041006);       //no ITSrefit requirement filter 256
+    trackCutsReject = CreateTrackCutsPWG4(10001006); //ITSrefit requirement filter 16
     trackCutsReject->SetEtaRange(etamin,etamax);
     trackCutsReject->SetPtRange(0.15, 1e10);
   }

diff --git a/PWG4/macros/AddTaskPWG4HighPtTrackQA.C b/PWG4/macros/AddTaskPWG4HighPtTrackQA.C
index a90331c4ef0030fe1a050b006edd232ecb20a2e6..84256090ba6f7c9a0fc3993a6c6d502b41c9e344 100644 (file)
--- a/PWG4/macros/AddTaskPWG4HighPtTrackQA.C
+++ b/PWG4/macros/AddTaskPWG4HighPtTrackQA.C
@@ -119,27 +119,27 @@ AliPWG4HighPtTrackQA* AddTaskPWG4HighPtTrackQA(char *prodType = "LHC10e14",Bool_
   }
   if(trackType==0 && cuts==1) {
     //Cuts global tracks with ITSrefit requirement and SPDrequirement for jet analysis
-    trackCuts = CreateTrackCutsPWG4(10001005);
+    trackCuts = CreateTrackCutsPWG4(10001006);
   }
   if(trackType==0 && cuts==2) {
     //Cuts global tracks with ITSrefit requirement but without SPD
-    trackCuts = CreateTrackCutsPWG4(10011005);
+    trackCuts = CreateTrackCutsPWG4(10011006);
   }
   if(trackType==7 && cuts==0) {
     // tight global tracks
-    trackCuts = CreateTrackCutsPWG4(10041005);
-    trackCutsReject = CreateTrackCutsPWG4(1005);
+    trackCuts = CreateTrackCutsPWG4(10041006);
+    trackCutsReject = CreateTrackCutsPWG4(1006);
     trackCutsReject->SetEtaRange(-0.9,0.9);
     trackCutsReject->SetPtRange(0.15, 1e10);
   }
   if(trackType==7 && cuts==1) {
     // tight global tracks
-    trackCuts = CreateTrackCutsPWG4(10011005);
+    trackCuts = CreateTrackCutsPWG4(10011006);
   }
   if(trackType==7 && cuts==2) {
     // no requirements on SPD and ITSrefit failed
-    trackCuts = CreateTrackCutsPWG4(10041005);       //no ITSrefit requirement filter 256
-    trackCutsReject = CreateTrackCutsPWG4(10001005); //ITSrefit requirement filter 16
+    trackCuts = CreateTrackCutsPWG4(10041006);       //no ITSrefit requirement filter 256
+    trackCutsReject = CreateTrackCutsPWG4(10001006); //ITSrefit requirement filter 16
     trackCutsReject->SetEtaRange(-0.9,0.9);
     trackCutsReject->SetPtRange(0.15, 1e10);
   }

diff --git a/PWG4/macros/CreateTrackCutsPWG4.C b/PWG4/macros/CreateTrackCutsPWG4.C
index e3dc7fb04577ca1e7328e84bca7c76eae1d31d7b..332d9e305ed56351f0fdb2786f6bbf3176391f2f 100644 (file)
--- a/PWG4/macros/CreateTrackCutsPWG4.C
+++ b/PWG4/macros/CreateTrackCutsPWG4.C
@@ -31,6 +31,7 @@ AliESDtrackCuts *CreateTrackCutsPWG4(Int_t cutMode) {
     trackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.8);
     trackCuts->SetMaxChi2PerClusterITS(36);
     trackCuts->SetMaxFractionSharedTPCClusters(0.4);
+    trackCuts->SetMaxChi2TPCConstrainedGlobal(36);
 
     trackCuts->SetEtaRange(-0.9,0.9);
     trackCuts->SetPtRange(0.15, 1e10);
@@ -173,6 +174,37 @@ AliESDtrackCuts *CreateTrackCutsPWG4(Int_t cutMode) {
     tag = "Global tracks jet analysis with ITSrefit and NclsIter1=70, noSPD requirement, no upper pt cut";
 
   }
+  if(stdCutMode == 1006) {
+
+    bStdCutsDefined = kTRUE;
+
+    // TPC  
+    TFormula *f1NClustersTPCLinearPtDep = new TFormula("f1NClustersTPCLinearPtDep","70.+30./20.*x");
+    trackCuts->SetMinNClustersTPCPtDep(f1NClustersTPCLinearPtDep,20.);
+    trackCuts->SetMinNClustersTPC(70);
+    trackCuts->SetMaxChi2PerClusterTPC(4);
+    trackCuts->SetRequireTPCStandAlone(kTRUE); //cut on NClustersTPC and chi2TPC Iter1
+    trackCuts->SetAcceptKinkDaughters(kFALSE);
+    trackCuts->SetRequireTPCRefit(kTRUE);
+    trackCuts->SetMaxFractionSharedTPCClusters(0.4);
+    // ITS
+    trackCuts->SetRequireITSRefit(kTRUE);
+    //accept secondaries
+    trackCuts->SetMaxDCAToVertexXY(2.4);
+    trackCuts->SetMaxDCAToVertexZ(3.2);
+    trackCuts->SetDCAToVertex2D(kTRUE);
+    //reject fakes
+    trackCuts->SetMaxChi2PerClusterITS(36);
+    trackCuts->SetMaxChi2TPCConstrainedGlobal(36);
+
+    trackCuts->SetRequireSigmaToVertex(kFALSE);
+
+    trackCuts->SetEtaRange(-0.9,0.9);
+    trackCuts->SetPtRange(0.15, 1E+15.);
+ 
+    tag = "Global tracks jet analysis with ITSrefit and NclsIter1=PtDep, noSPD requirement, no upper pt cut, golden chi2";
+
+  }