end-of-line normalization

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

diff --git a/PWGCF/EBYE/macros/configBalanceFunctionAnalysis.C b/PWGCF/EBYE/macros/configBalanceFunctionAnalysis.C
index d5ed4dbee3d702c5a29bc638ccf75fe5efc87300..8a020aff84018ed46bdb55dc37d24c4cee44fa9d 100644 (file)
--- a/PWGCF/EBYE/macros/configBalanceFunctionAnalysis.C
+++ b/PWGCF/EBYE/macros/configBalanceFunctionAnalysis.C
@@ -1,67 +1,68 @@
-//__________________________________________________//\r
-AliBalance *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 bHBTcut = kFALSE,\r
-                                    Bool_t bConversionCut = kFALSE) {\r
-  //Function to setup the AliBalance object and return it\r
-  AliBalance *gBalance = new AliBalance();\r
-  gBalance->SetAnalysisLevel(analysisLevel);\r
-  gBalance->SetShuffle(bShuffle);\r
-  gBalance->SetHBTcut(bHBTcut);\r
-  gBalance->SetConversionCut(bConversionCut);\r
-  if(centralityName) gBalance->SetCentralityIdentifier(centralityName);\r
-  gBalance->SetCentralityInterval(centrMin,centrMax);\r
-\r
-  //Set all analyses separately\r
-  //Rapidity\r
-  gBalance->SetInterval(AliBalance::kRapidity,-0.8,0.8,64,0.0,1.6);  \r
-  //Eta\r
-  gBalance->SetInterval(AliBalance::kEta,-0.8,0.8,64,0.0,1.6);\r
-  //Qlong\r
-  gBalance->SetInterval(AliBalance::kQlong,-1,1,100,0.0,2.0);\r
-  //Qout\r
-  gBalance->SetInterval(AliBalance::kQout,-1,1,100,0.0,2.0);\r
-  //Qside\r
-  gBalance->SetInterval(AliBalance::kQside,-1,1,100,0.0,2.0);\r
-  //Qinv\r
-  gBalance->SetInterval(AliBalance::kQinv,-1,1,100,0.0,2.0);\r
-  //Phi\r
-  gBalance->SetInterval(AliBalance::kPhi,0.,360.,90,0.,180.0);\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
+//__________________________________________________//
+AliBalance *GetBalanceFunctionObject(const char* analysisLevel = "AOD", 
+                                    const char* centralityName = 0x0,
+                                    Double_t centrMin = 0.,
+                                    Double_t centrMax = 100.,
+                                    Bool_t bShuffle = kFALSE,
+                                    Bool_t bHBTcut = kFALSE,
+                                    Bool_t bConversionCut = kFALSE,
+                                    Bool_t kUsePID = kTRUE) {
+  //Function to setup the AliBalance object and return it
+  AliBalance *gBalance = new AliBalance();
+  gBalance->SetAnalysisLevel(analysisLevel);
+  gBalance->SetShuffle(bShuffle);
+  gBalance->SetHBTcut(bHBTcut);
+  gBalance->SetConversionCut(bConversionCut);
+  if(centralityName) gBalance->SetCentralityIdentifier(centralityName);
+  gBalance->SetCentralityInterval(centrMin,centrMax);
+
+  //Set all analyses separately
+  //Rapidity
+  gBalance->SetInterval(AliBalance::kRapidity,-0.8,0.8,64,0.0,1.6);  
+  //Eta
+  gBalance->SetInterval(AliBalance::kEta,-0.8,0.8,64,0.0,1.6);
+  //Qlong
+  gBalance->SetInterval(AliBalance::kQlong,-1,1,100,0.0,2.0);
+  //Qout
+  gBalance->SetInterval(AliBalance::kQout,-1,1,100,0.0,2.0);
+  //Qside
+  gBalance->SetInterval(AliBalance::kQside,-1,1,100,0.0,2.0);
+  //Qinv
+  gBalance->SetInterval(AliBalance::kQinv,-1,1,100,0.0,5.0);
+  //Phi
+  gBalance->SetInterval(AliBalance::kPhi,0.,360.,90,0.,180.0);
+
+  //Init the histograms
+  gBalance->InitHistograms();
+  
+  return gBalance;
+}
+
+//__________________________________________________//
+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) {
+
+  // only used for ESDs
+  // Function to setup the AliESDtrackCuts object and return it
+  AliESDtrackCuts *cuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
+  cuts->SetRequireTPCStandAlone(kTRUE); // TPC only cuts!  
+
+  // extra TPC cuts (Syst studies)
+  if(minNClustersTPC != -1)  cuts->SetMinNClustersTPC(minNClustersTPC);
+  else cuts->SetMinNClustersTPC(70); // standard for filter bit 128
+  
+  if(maxTPCchi2 != -1) cuts->SetMaxChi2PerClusterTPC(maxTPCchi2);
+
+  // extra DCA cuts (Syst studies)  
+  if(maxDCAz!=-1 && maxDCAxy != -1){
+    cuts->SetMaxDCAToVertexZ(maxDCAz);
+    cuts->SetMaxDCAToVertexXY(maxDCAxy);
+  }
+
+  cuts->SetPtRange(ptMin,ptMax);
+  cuts->SetEtaRange(etaMin,etaMax);
+  cuts->DefineHistograms(1);
+  //cuts->SaveHistograms("trackCuts");
+
+  return cuts;
+}
+