]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGLF/STRANGENESS/LambdaK0PbPb/MainStreamAnalysis/FitSpectrum.C
Split: removed dirs now in AliPhysics
[u/mrichter/AliRoot.git] / PWGLF / STRANGENESS / LambdaK0PbPb / MainStreamAnalysis / FitSpectrum.C
diff --git a/PWGLF/STRANGENESS/LambdaK0PbPb/MainStreamAnalysis/FitSpectrum.C b/PWGLF/STRANGENESS/LambdaK0PbPb/MainStreamAnalysis/FitSpectrum.C
deleted file mode 100644 (file)
index 8e4eb93..0000000
+++ /dev/null
@@ -1,119 +0,0 @@
-void FitSpectrum(const char* filename, const char * listName = "clambdak0Histo_00", const char* suffix = "test",Int_t ihist = 0, Int_t iparticle) {
-  
-  // load basic libs, needed to 
-  //gROOT->LoadMacro("run.C");
-  gROOT->LoadMacro("runQA.C");
-  InitAndLoadLibs();
-
-  // Load Lee's Macro
-  gROOT->LoadMacro("AliMassFitControl.h+g");
-  gROOT->LoadMacro("PtMassAna2.C");
-  gROOT->LoadMacro("MultYields2QA.C");
-
-  char* histName = 0;
-  switch (iparticle) {
-  case 0:
-    switch (ihist){
-          case 0:
-               histName = "h2PtVsMassK0";
-                break;
-          case 1:
-               histName = "h2DcaPosToPrimVertexK0vsMassK0";
-               break;
-          case 2:
-               histName = "h2DcaNegToPrimVertexK0vsMassK0";
-               break;
-          case 3:
-               histName = "h2RadiusV0K0vsMassK0pt1";
-               break;
-          case 4:
-               histName = "h2DecayLengthV0K0vsMassK0pt1";   
-               break;
-          case 5:
-               histName = "h2DcaV0DaughtersK0vsMassK0";
-               break;
-          case 6:
-               histName = "h2CosPointAngleK0vsMassK0";
-               break;
-          default;
-          cout << "histogram " << histName << endl;
-          }
-    cout << "Particle "<< iparticle << " to be implemented" << endl;  
-    break;
-  case 1:
-    switch (ihist){
-          case 0:
-               histName = "h2PtVsMassLambda";   
-               break;
-          case 1:
-               histName = "h2DcaPosToPrimVertexLvsMassL";   
-               break;
-          case 2:
-               histName = "h2DcaNegToPrimVertexLvsMassL";   
-               break;
-          case 3:
-               histName = "h2RadiusV0LvsMassL";   
-               break;
-          case 4:
-               histName = "h2DecayLengthV0LvsMassL";   
-               break;
-          case 5:
-               histName = "h2DcaV0DaughtersLvsMassL";
-               break;
-          case 6:
-               histName = "h2CosPointAngleLvsMassL";
-               break;
-          default;
-          cout << "histogram " << histName << endl;
-          }
-    cout << "Particle "<< iparticle << " to be implemented" << endl; 
-    break; 
-  case 2:
-    switch (ihist){
-          case 0:
-               histName = "h2PtVsMassAntiLambda";   
-               break;
-          case 1:
-               histName = "h2DcaPosToPrimVertexLvsMassL";   
-               break;
-          case 2:
-               histName = "h2DcaNegToPrimVertexLvsMassL";   
-               break;
-          case 3:
-               histName = "h2RadiusV0LvsMassL";   
-               break;
-          case 4:
-               histName = "h2DecayLengthV0LvsMassL";   
-               break;
-          case 5:
-               histName = "h2DcaV0DaughtersLvsMassL";
-               break;
-          case 6:
-               histName = "h2CosPointAngleLvsMassL";
-               break;
-          default;
-          cout << "histogram " << histName << endl;
-          }
-    cout << "Particle "<< iparticle << " to be implemented" << endl; 
-    break; 
-
-  }
-
-  TFile *file = new TFile(filename);
-  TList *list = file->Get(listName); 
-
-  TH2F * h2 = (TH2F*) list->FindObject(histName);
-  h2->Draw();
-  
-  ///// iNorm is used by MultYields2QA to normalize the distributions
-  TH1F * h1 = (TH1F*) list->FindObject("h1PrimaryVertexZ");
-  Int_t iNorm = h1->GetEntries();
-  cout << "number of entries with Zvertex < |10cm|: " << iNorm;
-
-  TString suffixFull = histName;
-  if(strlen(suffix)) suffixFull = suffixFull + "_" + suffix;
-  //MultYields3((TH3F*)h2,iparticle,0,suffixFull); // FIXME: modify MultYields2 to handle 1D histos
-  MultYields2QA(h2,iparticle,ihist,iNorm,0,suffixFull); // FIXME: modify MultYields2 to handle 1D histos
-
-  
-}