PWGCF/EBYE/macros/configBalanceFunctionAnalysis.C

   1 //__________________________________________________//
   2 AliBalance *GetBalanceFunctionObject(const char* analysisLevel = "AOD",
   3                                      const char* centralityName = 0x0,
   4                                      Double_t centrMin = 0.,
   5                                      Double_t centrMax = 100.,
   6                                      Bool_t bShuffle = kFALSE,
   7                                      Bool_t bHBTcut = kFALSE,
   8                                      Bool_t bConversionCut = kFALSE,
   9                                      Bool_t kUsePID = kTRUE) {
  10   //Function to setup the AliBalance object and return it
  11   AliBalance *gBalance = new AliBalance();
  12   gBalance->SetAnalysisLevel(analysisLevel);
  13   gBalance->SetShuffle(bShuffle);
  14   gBalance->SetHBTcut(bHBTcut);
  15   gBalance->SetConversionCut(bConversionCut);
  16   if(centralityName) gBalance->SetCentralityIdentifier(centralityName);
  17   gBalance->SetCentralityInterval(centrMin,centrMax);
  18
  19   //Set all analyses separately
  20   //Rapidity
  21   gBalance->SetInterval(AliBalance::kRapidity,-0.8,0.8,64,0.0,1.6);
  22   //Eta
  23   gBalance->SetInterval(AliBalance::kEta,-0.8,0.8,64,0.0,1.6);
  24   //Qlong
  25   gBalance->SetInterval(AliBalance::kQlong,-1,1,100,0.0,2.0);
  26   //Qout
  27   gBalance->SetInterval(AliBalance::kQout,-1,1,100,0.0,2.0);
  28   //Qside
  29   gBalance->SetInterval(AliBalance::kQside,-1,1,100,0.0,2.0);
  30   //Qinv
  31   gBalance->SetInterval(AliBalance::kQinv,-1,1,100,0.0,5.0);
  32   //Phi
  33   gBalance->SetInterval(AliBalance::kPhi,0.,360.,90,0.,180.0);
  34
  35   //Init the histograms
  36   gBalance->InitHistograms();
  37
  38   return gBalance;
  39 }
  40
  41 //__________________________________________________//
  42 AliESDtrackCuts *GetTrackCutsObject(Double_t ptMin, Double_t ptMax, Double_t etaMin, Double_t etaMax, Double_t maxTPCchi2, Double_t maxDCAz, Double_t maxDCAxy, Int_t minNClustersTPC) {
  43
  44   // only used for ESDs
  45   // Function to setup the AliESDtrackCuts object and return it
  46   AliESDtrackCuts *cuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
  47   cuts->SetRequireTPCStandAlone(kTRUE); // TPC only cuts!
  48
  49   // extra TPC cuts (Syst studies)
  50   if(minNClustersTPC != -1)  cuts->SetMinNClustersTPC(minNClustersTPC);
  51   else cuts->SetMinNClustersTPC(70); // standard for filter bit 128
  52
  53   if(maxTPCchi2 != -1) cuts->SetMaxChi2PerClusterTPC(maxTPCchi2);
  54
  55   // extra DCA cuts (Syst studies)
  56   if(maxDCAz!=-1 && maxDCAxy != -1){
  57     cuts->SetMaxDCAToVertexZ(maxDCAz);
  58     cuts->SetMaxDCAToVertexXY(maxDCAxy);
  59   }
  60
  61   cuts->SetPtRange(ptMin,ptMax);
  62   cuts->SetEtaRange(etaMin,etaMax);
  63   cuts->DefineHistograms(1);
  64   //cuts->SaveHistograms("trackCuts");
  65
  66   return cuts;
  67 }
  68