PWG0/CreateStandardCuts.C

   1 /* $Id: CreateCuts.C,v 1.5 2008/01/11 08:28:52 jgrosseo Exp $ */
   2
   3 // this macro creates the track cuts used in this analysis
   4
   5 AliESDtrackCuts* CreateTrackCuts(AliPWG0Helper::AnalysisMode analysisMode,  Bool_t hists = kTRUE,  Float_t ptMin = 0,  Float_t etacut =1e10)
   6 {
   7   AliESDtrackCuts* esdTrackCuts = 0;
   8
   9   // see https://twiki.cern.ch/twiki/bin/view/ALICE/SelectionOfPrimaryTracksForPp2009DataAnalysis
  10
  11   if (analysisMode & AliPWG0Helper::kTPC)
  12   {
  13     esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts");
  14
  15     TString tag("TPC-only tracking");
  16
  17     esdTrackCuts->SetMaxDCAToVertexZ(3.2);
  18     esdTrackCuts->SetMaxDCAToVertexXY(2.4);
  19     esdTrackCuts->SetDCAToVertex2D(kTRUE);
  20
  21     esdTrackCuts->SetRequireTPCRefit(kTRUE);
  22     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
  23     esdTrackCuts->SetMinNClustersTPC(70);
  24     esdTrackCuts->SetMaxChi2PerClusterTPC(4);
  25   }
  26
  27   if (analysisMode & AliPWG0Helper::kTPCITS);
  28   {
  29     esdTrackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2009(kTRUE);
  30     esdTrackCuts->SetPtRange(ptMin);  // adding pt cut
  31     esdTrackCuts->SetEtaRange(-etacut,etacut);
  32     TString tag("Global tracking");
  33   }
  34   if ( analysisMode & AliPWG0Helper::kTPCSPD) {
  35
  36     esdTrackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2009(kFALSE);
  37     TString tag("Global tracking+tracklets");
  38
  39   }
  40
  41   if (hists)
  42     esdTrackCuts->DefineHistograms(1);
  43
  44   // cuts for data without field
  45   if (!(analysisMode & AliPWG0Helper::kFieldOn))
  46   {
  47     tag += " without field";
  48   }
  49
  50   Printf("Created track cuts for: %s", tag.Data());
  51
  52   return esdTrackCuts;
  53 }
  54