adding momentum difference cut (suppression of HBT effects)

[u/mrichter/AliRoot.git] / PWGCF / EBYE / macros / configBalanceFunctionPsiAnalysis.C

//__________________________________________________//\r
AliBalancePsi *GetBalanceFunctionObject(const char* analysisLevel = "ESD", \r
					const char* centralityName = 0x0,\r
					Double_t centrMin = 0.,\r
					Double_t centrMax = 100.,\r
					Bool_t bShuffle = kFALSE,\r
					Bool_t bResonancesCut = kFALSE,\r
					Bool_t bHBTCut = kFALSE,\r
					Bool_t bConversionCut = kFALSE,\r
					Bool_t bMomentumDifferenceCut = kFALSE,\r
					TString fArgEventClass = "EventPlane",\r
					Double_t deltaEtaMax = 2.0,\r
					Bool_t bVertexBinning = kFALSE) {\r
  //Function to setup the AliBalance object and return it\r
  AliBalancePsi *gBalance = new AliBalancePsi();\r
  gBalance->SetAnalysisLevel(analysisLevel);\r
  gBalance->SetShuffle(bShuffle);\r
  if(bResonancesCut) gBalance->UseResonancesCut();\r
  if(bHBTCut) gBalance->UseHBTCut();\r
  if(bConversionCut) gBalance->UseConversionCut();\r
  if(bMomentumDifferenceCut) gBalance->UseMomentumDifferenceCut();\r
  if(centralityName) gBalance->SetCentralityIdentifier(centralityName);\r
  if(bVertexBinning) gBalance->SetVertexZBinning();\r
  gBalance->SetCentralityInterval(centrMin,centrMax);\r
  gBalance->SetEventClass(fArgEventClass);\r
  gBalance->SetDeltaEtaMax(deltaEtaMax);\r
\r
  //Set all analyses separately\r
  //Rapidity\r
  //gBalance->SetInterval(AliBalance::kRapidity,-0.8,0.8,32,-1.6,1.6,15.);  \r
  //Eta\r
  //gBalance->SetInterval(AliBalance::kEta,-0.8,0.8,32,-1.6,1.6,15);\r
  //Qlong\r
  //gBalance->SetInterval(AliBalance::kQlong,-1,1,200,0.0,4.0,15);\r
  //Qout\r
  //gBalance->SetInterval(AliBalance::kQout,-1,1,200,0.0,4.0,15);\r
  //Qside\r
  //gBalance->SetInterval(AliBalance::kQside,-1,1,200,0.0,4.0,15);\r
  //Qinv\r
  //gBalance->SetInterval(AliBalance::kQinv,-1,1,200,0.0,4.0,15);\r
  //Phi\r
  //gBalance->SetInterval(AliBalance::kPhi,0.,360.,90,-180.,180.0,15);\r
\r
  //Init the histograms\r
  gBalance->InitHistograms();\r
  \r
  return gBalance;\r
}\r
\r
//__________________________________________________//\r
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) {\r
\r
  // only used for ESDs\r
  // Function to setup the AliESDtrackCuts object and return it\r
  AliESDtrackCuts *cuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();\r
  cuts->SetRequireTPCStandAlone(kTRUE); // TPC only cuts!  \r
\r
  // extra TPC cuts (Syst studies)\r
  if(minNClustersTPC != -1)  cuts->SetMinNClustersTPC(minNClustersTPC);\r
  else cuts->SetMinNClustersTPC(70); // standard for filter bit 128\r
  \r
  if(maxTPCchi2 != -1) cuts->SetMaxChi2PerClusterTPC(maxTPCchi2);\r
\r
  // extra DCA cuts (Syst studies)  \r
  if(maxDCAz!=-1 && maxDCAxy != -1){\r
    cuts->SetMaxDCAToVertexZ(maxDCAz);\r
    cuts->SetMaxDCAToVertexXY(maxDCAxy);\r
  }\r
\r
  cuts->SetPtRange(ptMin,ptMax);\r
  cuts->SetEtaRange(etaMin,etaMax);\r
  cuts->DefineHistograms(1);\r
  //cuts->SaveHistograms("trackCuts");\r
\r
  return cuts;\r
}\r
\r