Merge branch 'feature-movesplit'

[u/mrichter/AliRoot.git] / PWGLF / FORWARD / analysis2 / AliPoissonCalculator.cxx

diff --git a/PWGLF/FORWARD/analysis2/AliPoissonCalculator.cxx b/PWGLF/FORWARD/analysis2/AliPoissonCalculator.cxx
index db24e2a82d2c534e0bef6a13466841dd34614651..c37cd66595bb94be9e7814c7d17d740670548b01 100644 (file)
--- a/PWGLF/FORWARD/analysis2/AliPoissonCalculator.cxx
+++ b/PWGLF/FORWARD/analysis2/AliPoissonCalculator.cxx
@@ -1,5 +1,5 @@
 #include "AliPoissonCalculator.h"
-#include "AliForwardCorrectionManager.h"
+// #include "AliForwardCorrectionManager.h"
 #include <TH2D.h>
 #include <TBrowser.h>
 #include <TROOT.h>
@@ -13,7 +13,7 @@
 // using Poisson statistics. 
 //
 // The input is assumed to be binned in @f$(x,y)@f$ as described by
-// the 2D histogram passwd to the Reset member function.
+// the 2D histogram passed to the Reset member function.
 //
 // The data is grouped in to regions as defined by the parameters
 // fXLumping and fYLumping.  The total number of cells and number of
@@ -179,8 +179,8 @@ AliPoissonCalculator::Define(const TAxis& xaxis, const TAxis& yaxis)
   Reset(fBasic);
 }
 //____________________________________________________________________
-void AliPoissonCalculator::MakeOutput() {
-
+void AliPoissonCalculator::MakeOutput() 
+{
   Int_t n = fXLumping * fYLumping + 1;
   fEmptyVsTotal = new TH2D("emptyVsTotal", 
                           "# of empty # bins vs total # bins", 
@@ -216,9 +216,7 @@ void AliPoissonCalculator::MakeOutput() {
   fCorr->SetYTitle("Correction 1/(1-e^{#bar{N_{c}}})");
   fCorr->SetZTitle("Events");
   fCorr->SetOption("colz");
-  fCorr->SetDirectory(0);
-  
- 
+  fCorr->SetDirectory(0); 
 }
 //____________________________________________________________________
 void
@@ -408,6 +406,13 @@ AliPoissonCalculator::Result(Bool_t correct)
       fBasic->SetBinError(ix,iy,poissonE);
     }
   }
+  return fBasic;
+}
+  
+//____________________________________________________________________
+void
+AliPoissonCalculator::FillDiagnostics()
+{
   for (Int_t ix = 1; ix <= fEmpty->GetNbinsX(); ix++) { 
     for (Int_t iy = 1; iy <= fEmpty->GetNbinsY(); iy++) { 
       Double_t empty    = fEmpty->GetBinContent(ix, iy);
@@ -421,9 +426,7 @@ AliPoissonCalculator::Result(Bool_t correct)
       fCorr->Fill(mean, corr);
     }
   }
-  return fBasic;
 }
-  
 //____________________________________________________________________
 void
 AliPoissonCalculator::Print(const Option_t*) const