]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGLF/SPECTRA/AntiprotonToProton/configProtonAnalysisBaseObject.C
Split: removed dirs now in AliPhysics
[u/mrichter/AliRoot.git] / PWGLF / SPECTRA / AntiprotonToProton / configProtonAnalysisBaseObject.C
diff --git a/PWGLF/SPECTRA/AntiprotonToProton/configProtonAnalysisBaseObject.C b/PWGLF/SPECTRA/AntiprotonToProton/configProtonAnalysisBaseObject.C
deleted file mode 100644 (file)
index 32b9135..0000000
+++ /dev/null
@@ -1,138 +0,0 @@
-AliProtonAnalysisBase *GetProtonAnalysisBaseObject(const char* analysisLevel = "ESD",
-                                                  Bool_t kAnalyzeMC = kTRUE,
-                                                  const char* esdAnalysisType = "Hybrid",
-                                                  const char* pidMode = "Bayesian",
-                                                  Bool_t kUseOnlineTrigger = kFALSE,
-                                                  Bool_t kUseOfflineTrigger = kFALSE,
-                                                  Bool_t kRunQA = kFALSE) {
-  //Function to setup the AliProtonAnalysisBase object and return it
-  AliProtonAnalysisBase *baseAnalysis = new AliProtonAnalysisBase();
-  //baseAnalysis->SetDebugMode();
-  if(kRunQA) baseAnalysis->SetRunQA();
-  baseAnalysis->SetAnalysisLevel(analysisLevel);
-  if(analysisLevel == "ESD") {
-    if(kAnalyzeMC)
-      baseAnalysis->SetTriggerMode(AliProtonAnalysisBase::kMB2);
-    //use the offline trigger
-    if(kUseOnlineTrigger) baseAnalysis->UseOnlineTrigger();
-
-    //use the offline trigger
-    if(kUseOfflineTrigger) baseAnalysis->OfflineTriggerInit();
-
-    baseAnalysis->SetMinTPCClusters(80);
-    baseAnalysis->SetMaxChi2PerTPCCluster(3.5);
-    /*baseAnalysis->SetMaxCov11(2.0);
-    baseAnalysis->SetMaxCov22(2.0);
-    baseAnalysis->SetMaxCov33(0.5);
-    baseAnalysis->SetMaxCov44(0.5);
-    baseAnalysis->SetMaxCov55(2.0);*/
-    baseAnalysis->SetMinTPCdEdxPoints(80);
-    switch(esdAnalysisType) {
-    case "TPC":
-      baseAnalysis->SetAnalysisMode(AliProtonAnalysisBase::kTPC);
-      baseAnalysis->SetPhaseSpace(10, -0.5, 0.5, 16, 0.5, 0.9);
-      baseAnalysis->SetTPCpid();
-      baseAnalysis->SetMaxSigmaToVertexTPC(2.0);
-      //baseAnalysis->SetMaxDCAXYTPC(1.5);
-      //baseAnalysis->SetMaxDCAZTPC(1.5);
-      break;
-    case "Hybrid":
-      baseAnalysis->SetAnalysisMode(AliProtonAnalysisBase::kHybrid);
-      baseAnalysis->SetPhaseSpace(9, -0.9, 0.9, 6, 0.45, 1.05);
-      //baseAnalysis->SetPhaseSpace(18, -0.9, 0.9, 32, 0.5, 1.3);
-      baseAnalysis->SetTPCpid();
-      //baseAnalysis->SetMaxSigmaToVertex(3.0);
-      //baseAnalysis->SetMaxDCAXY(0.5);
-      //baseAnalysis->SetMaxDCAZ(0.7);
-      baseAnalysis->SetMaxDCA3D(2.0);
-      //baseAnalysis->SetPointOnITSLayer6();
-      //baseAnalysis->SetPointOnITSLayer5();
-      //baseAnalysis->SetPointOnITSLayer4();
-      //baseAnalysis->SetPointOnITSLayer3();
-      //baseAnalysis->SetPointOnITSLayer2();
-      //baseAnalysis->SetPointOnITSLayer1();
-      baseAnalysis->SetPointOnSPDLayers();
-      baseAnalysis->SetMinITSClusters(2);
-      break;
-    case "FullHybrid":
-      baseAnalysis->SetAnalysisMode(AliProtonAnalysisBase::kFullHybrid);
-      baseAnalysis->SetPhaseSpace(9, -0.9, 0.9, 6, 0.45, 1.05);
-      //baseAnalysis->SetPhaseSpace(18, -0.9, 0.9, 32, 0.5, 1.3);
-      baseAnalysis->SetTPCpid();
-      //baseAnalysis->SetMaxSigmaToVertex(3.0);
-      //baseAnalysis->SetMaxDCAXY(0.2);
-      //baseAnalysis->SetMaxDCAZ(0.7);
-      //baseAnalysis->SetMaxDCA3D(0.2);
-      baseAnalysis->SetPtDependentDCAxy(5,2.89575e+02,6.62161e+01,1.99085e+00);
-      //baseAnalysis->SetPointOnITSLayer6();
-      //baseAnalysis->SetPointOnITSLayer5();
-      //baseAnalysis->SetPointOnITSLayer4();
-      //baseAnalysis->SetPointOnITSLayer3();
-      //baseAnalysis->SetPointOnITSLayer2();
-      //baseAnalysis->SetPointOnITSLayer1();
-      baseAnalysis->SetPointOnSPDLayers();
-      baseAnalysis->SetMinITSClusters(2);
-      break;
-    case "Global":
-      baseAnalysis->SetAnalysisMode(AliProtonAnalysisBase::kGlobal);
-      baseAnalysis->SetPhaseSpace(20, -1.0, 1.0, 48, 0.3, 1.5);
-      baseAnalysis->SetMaxSigmaToVertex(2.0);
-      //baseAnalysis->SetMaxDCAXY(2.0);
-      //baseAnalysis->SetMaxDCAZ(2.0);
-      baseAnalysis->SetTPCRefit();
-      baseAnalysis->SetPointOnITSLayer1();
-      baseAnalysis->SetPointOnITSLayer2();
-      //baseAnalysis->SetPointOnITSLayer3();
-      //baseAnalysis->SetPointOnITSLayer4();
-      baseAnalysis->SetPointOnITSLayer5();
-      baseAnalysis->SetPointOnITSLayer6();
-      baseAnalysis->SetMinITSClusters(5);
-      baseAnalysis->SetITSRefit();
-      baseAnalysis->SetESDpid();
-      baseAnalysis->SetTOFpid();
-      break;
-    default:
-      break;
-    }
-    baseAnalysis->SetAcceptedVertexDiamond(1.,1.,10.);
-    baseAnalysis->SetMinNumOfContributors(0);
-    //baseAnalysis->SetEtaMode();
-    switch(pidMode) {
-    case "Bayesian":
-      baseAnalysis->SetPIDMode(AliProtonAnalysisBase::kBayesian);
-      //Momentum dependent priors
-      /*TFile *f = TFile::Open("$ALICE_ROOT/PWG2/data/PriorProbabilities.root ");
-       TF1 *fitElectrons = (TF1 *)f->Get("fitElectrons");
-       TF1 *fitMuons = (TF1 *)f->Get("fitMuons");
-       TF1 *fitPions = (TF1 *)f->Get("fitPions");
-       TF1 *fitKaons = (TF1 *)f->Get("fitKaons");
-       TF1 *fitProtons = (TF1 *)f->Get("fitProtons");
-       baseAnalysis->SetPriorProbabilityFunctions(fitElectrons,
-       fitMuons,
-       fitPions,
-       fitKaons,
-       fitProtons);*/
-      //Fixed prior probabilities
-      Double_t partFrac[5] = {0.01, 0.01, 0.85, 0.10, 0.05};
-      if(!baseAnalysis->IsPriorProbabilityFunctionUsed())
-       baseAnalysis->SetPriorProbabilities(partFrac);
-      break;
-    case "Ratio":
-      baseAnalysis->SetPIDMode(AliProtonAnalysisBase::kRatio);
-      baseAnalysis->SetRatio(-0.2);
-      break;
-    case "Sigma":
-      baseAnalysis->SetPIDMode(AliProtonAnalysisBase::kSigma1);
-      baseAnalysis->SetNSigma(4);
-      break;
-    default:
-      break;
-    }//PID mode
-  }//ESD
-  if(analysisLevel == "MC") 
-    baseAnalysis->SetPhaseSpace(10, -0.5, 0.5, 16, 0.5, 0.9);
-  if(analysisLevel == "AOD")
-    baseAnalysis->SetPhaseSpace(10, -0.5, 0.5, 16, 0.5, 0.9);
-
-  return baseAnalysis;
-}