#include "AliPoissonCalculator.h"
-#include "AliForwardCorrectionManager.h"
+// #include "AliForwardCorrectionManager.h"
#include <TH2D.h>
#include <TBrowser.h>
#include <TROOT.h>
// 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
Reset(fBasic);
}
//____________________________________________________________________
-void AliPoissonCalculator::MakeOutput() {
-
+void AliPoissonCalculator::MakeOutput()
+{
Int_t n = fXLumping * fYLumping + 1;
fEmptyVsTotal = new TH2D("emptyVsTotal",
"# of empty # bins vs total # bins",
fCorr->SetYTitle("Correction 1/(1-e^{#bar{N_{c}}})");
fCorr->SetZTitle("Events");
fCorr->SetOption("colz");
- fCorr->SetDirectory(0);
-
-
+ fCorr->SetDirectory(0);
}
//____________________________________________________________________
void
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);
Double_t corr = CalculateCorrection(empty, total);
fEmptyVsTotal->Fill(total, empty);
fMean->Fill(mean);
- fOcc->Fill(100 * (1 - empty/total));
+ if (total > 1e-6) fOcc->Fill(100 * (1 - empty/total));
//Old fOcc->Fill(100 * (1 - TMath::PoissonI(0,mean)));
fCorr->Fill(mean, corr);
}
}
- return fBasic;
}
-
//____________________________________________________________________
void
AliPoissonCalculator::Print(const Option_t*) const