Updating the parameters

[u/mrichter/AliRoot.git] / PWG2 / SPECTRA / macros / configProtonAnalysisBaseObject.C

diff --git a/PWG2/SPECTRA/macros/configProtonAnalysisBaseObject.C b/PWG2/SPECTRA/macros/configProtonAnalysisBaseObject.C
index 3d3273a974a3ed289dd12e303e53a58e96e13d32..0cf5e673eac81b87256222c18b3b0a80fc9f6a36 100644 (file)
--- a/PWG2/SPECTRA/macros/configProtonAnalysisBaseObject.C
+++ b/PWG2/SPECTRA/macros/configProtonAnalysisBaseObject.C
@@ -1,20 +1,29 @@
 AliProtonAnalysisBase *GetProtonAnalysisBaseObject(const char* analysisLevel = "ESD",\r
+                                                  Bool_t kAnalyzeMC = kTRUE,\r
                                                   const char* esdAnalysisType = "Hybrid",\r
-                                                  const char* pidMode = "Bayesian") {\r
-  \r
+                                                  const char* pidMode = "Bayesian",\r
+                                                  Bool_t kUseOnlineTrigger = kFALSE,\r
+                                                  Bool_t kUseOfflineTrigger = kFALSE) {\r
   //Function to setup the AliProtonAnalysisBase object and return it\r
   AliProtonAnalysisBase *baseAnalysis = new AliProtonAnalysisBase();\r
   //baseAnalysis->SetDebugMode();\r
   baseAnalysis->SetAnalysisLevel(analysisLevel);\r
-  if(analysisLevel == "ESD") {  \r
-    baseAnalysis->SetTriggerMode(AliProtonAnalysisBase::kMB2);\r
-    baseAnalysis->SetMinTPCClusters(110);\r
-    baseAnalysis->SetMaxChi2PerTPCCluster(2.2);\r
-    baseAnalysis->SetMaxCov11(0.5);\r
-    baseAnalysis->SetMaxCov22(0.5);\r
+  if(analysisLevel == "ESD") {\r
+    if(kAnalyzeMC)\r
+      baseAnalysis->SetTriggerMode(AliProtonAnalysisBase::kMB2);\r
+    //use the offline trigger\r
+    if(kUseOnlineTrigger) baseAnalysis->UseOnlineTrigger();\r
+\r
+    //use the offline trigger\r
+    if(kUseOfflineTrigger) baseAnalysis->OfflineTriggerInit();\r
+\r
+    baseAnalysis->SetMinTPCClusters(80);\r
+    baseAnalysis->SetMaxChi2PerTPCCluster(3.5);\r
+    /*baseAnalysis->SetMaxCov11(2.0);\r
+    baseAnalysis->SetMaxCov22(2.0);\r
     baseAnalysis->SetMaxCov33(0.5);\r
     baseAnalysis->SetMaxCov44(0.5);\r
-    baseAnalysis->SetMaxCov55(0.5);\r
+    baseAnalysis->SetMaxCov55(2.0);*/\r
     baseAnalysis->SetMinTPCdEdxPoints(80);\r
     switch(esdAnalysisType) {\r
     case "TPC":\r
@@ -27,19 +36,39 @@ AliProtonAnalysisBase *GetProtonAnalysisBaseObject(const char* analysisLevel = "
       break;\r
     case "Hybrid":\r
       baseAnalysis->SetAnalysisMode(AliProtonAnalysisBase::kHybrid);\r
-      baseAnalysis->SetPhaseSpace(10, -0.5, 0.5, 16, 0.5, 0.9);\r
+      baseAnalysis->SetPhaseSpace(9, -0.9, 0.9, 6, 0.45, 1.05);\r
       //baseAnalysis->SetPhaseSpace(18, -0.9, 0.9, 32, 0.5, 1.3);\r
       baseAnalysis->SetTPCpid();\r
-      baseAnalysis->SetMaxSigmaToVertex(2.0);\r
-      /*baseAnalysis->SetMaxDCAXY(1.5);\r
-       baseAnalysis->SetMaxDCAZ(1.5);*/\r
-      baseAnalysis->SetPointOnITSLayer6();\r
-      baseAnalysis->SetPointOnITSLayer5();\r
+      //baseAnalysis->SetMaxSigmaToVertex(3.0);\r
+      //baseAnalysis->SetMaxDCAXY(0.5);\r
+      //baseAnalysis->SetMaxDCAZ(0.7);\r
+      baseAnalysis->SetMaxDCA3D(2.0);\r
+      //baseAnalysis->SetPointOnITSLayer6();\r
+      //baseAnalysis->SetPointOnITSLayer5();\r
       //baseAnalysis->SetPointOnITSLayer4();\r
       //baseAnalysis->SetPointOnITSLayer3();\r
-      baseAnalysis->SetPointOnITSLayer2();\r
-      baseAnalysis->SetPointOnITSLayer1();\r
-      baseAnalysis->SetMinITSClusters(4);\r
+      //baseAnalysis->SetPointOnITSLayer2();\r
+      //baseAnalysis->SetPointOnITSLayer1();\r
+      baseAnalysis->SetPointOnSPDLayers();\r
+      baseAnalysis->SetMinITSClusters(2);\r
+      break;\r
+    case "FullHybrid":\r
+      baseAnalysis->SetAnalysisMode(AliProtonAnalysisBase::kFullHybrid);\r
+      baseAnalysis->SetPhaseSpace(9, -0.9, 0.9, 6, 0.45, 1.05);\r
+      //baseAnalysis->SetPhaseSpace(18, -0.9, 0.9, 32, 0.5, 1.3);\r
+      baseAnalysis->SetTPCpid();\r
+      //baseAnalysis->SetMaxSigmaToVertex(3.0);\r
+      //baseAnalysis->SetMaxDCAXY(0.2);\r
+      //baseAnalysis->SetMaxDCAZ(0.7);\r
+      baseAnalysis->SetMaxDCA3D(0.2);\r
+      //baseAnalysis->SetPointOnITSLayer6();\r
+      //baseAnalysis->SetPointOnITSLayer5();\r
+      //baseAnalysis->SetPointOnITSLayer4();\r
+      //baseAnalysis->SetPointOnITSLayer3();\r
+      //baseAnalysis->SetPointOnITSLayer2();\r
+      //baseAnalysis->SetPointOnITSLayer1();\r
+      baseAnalysis->SetPointOnSPDLayers();\r
+      baseAnalysis->SetMinITSClusters(2);\r
       break;\r
     case "Global":\r
       baseAnalysis->SetAnalysisMode(AliProtonAnalysisBase::kGlobal);\r
@@ -62,8 +91,9 @@ AliProtonAnalysisBase *GetProtonAnalysisBaseObject(const char* analysisLevel = "
     default:\r
       break;\r
     }\r
-    baseAnalysis->SetAcceptedVertexDiamond(5.,5.,15.);\r
-    baseAnalysis->SetEtaMode();\r
+    baseAnalysis->SetAcceptedVertexDiamond(1.,1.,10.);\r
+    baseAnalysis->SetMinNumOfContributors(0);\r
+    //baseAnalysis->SetEtaMode();\r
     switch(pidMode) {\r
     case "Bayesian":\r
       baseAnalysis->SetPIDMode(AliProtonAnalysisBase::kBayesian);\r
@@ -86,16 +116,11 @@ AliProtonAnalysisBase *GetProtonAnalysisBaseObject(const char* analysisLevel = "
       break;\r
     case "Ratio":\r
       baseAnalysis->SetPIDMode(AliProtonAnalysisBase::kRatio);\r
+      baseAnalysis->SetRatio(-0.2);\r
       break;\r
-    case "Sigma1":\r
+    case "Sigma":\r
       baseAnalysis->SetPIDMode(AliProtonAnalysisBase::kSigma1);\r
-      baseAnalysis->SetNSigma(3);\r
-      baseAnalysis->SetdEdxBandInfo("$ALICE_ROOT/PWG2/data/protonsdEdxInfo.dat");\r
-      break;\r
-    case "Sigma2":\r
-      baseAnalysis->SetPIDMode(AliProtonAnalysisBase::kSigma2);\r
-      baseAnalysis->SetNSigma(3);\r
-      baseAnalysis->SetdEdxBandInfo("$ALICE_ROOT/PWG2/data/protonsdEdxInfo.dat");\r
+      baseAnalysis->SetNSigma(4);\r
       break;\r
     default:\r
       break;\r