AliPoissonCalculator.h

Go to the documentation of this file.

#ifndef ALIPOISSONCALCULATOR_H
#define ALIPOISSONCALCULATOR_H
#include <TNamed.h>
class TH2D;
class TH1D;
class TBrowser;
class TAxis;

class AliPoissonCalculator : public TNamed
{
public:
  AliPoissonCalculator(); 
  AliPoissonCalculator(const char*/*, UShort_t d, Char_t r*/);
  AliPoissonCalculator(const AliPoissonCalculator& o);

  virtual ~AliPoissonCalculator();
  AliPoissonCalculator& operator=(const AliPoissonCalculator& o);
  void SetLumping(UShort_t nx, UShort_t ny);
  void SetXLumping(UShort_t nx) { SetLumping(nx, fYLumping); } //*MENU*
  void SetYLumping(UShort_t ny) { SetLumping(fYLumping, ny); } //*MENU*
  void Init(Int_t xLumping=-1, Int_t yLumping=-1);
  
  void Define(const TAxis& xaxis, const TAxis& yaxis);
  void MakeOutput();
  void Output(TList* d);
  void Reset(const TH2D* base);
  void Fill(UShort_t strip, UShort_t sec, Bool_t hit, Double_t weight=1);
  TH2D* Result(Bool_t correct=true);
  void FillDiagnostics();
  Bool_t IsFolder() const { return kTRUE; }
  void Print(const Option_t* option="") const;
  void Browse(TBrowser* b);

  TH2D* GetEmptyVsTotal() const { return fEmptyVsTotal; }
  TH1D* GetMean() const { return fMean; }
  TH1D* GetOccupancy() const { return fOcc; }
  TH2D* GetCorrection() const { return fCorr; }

  Int_t GetReducedXBin(Int_t ix) const;
  Int_t GetReducedXBin(Double_t x) const;
  Int_t GetReducedYBin(Int_t iy) const;
  Int_t GetReducedYBin(Double_t y) const;

protected:
  Int_t CheckLumping(char which, Int_t nBins, Int_t lumping) const;
  void CleanUp();
  Double_t CalculateMean(Double_t empty, Double_t total) const;
  Double_t CalculateCorrection(Double_t empty, Double_t total) const;
  UShort_t fXLumping;   // Grouping of eta bins 
  UShort_t fYLumping;   // Grouping of phi bins 
  TH2D*    fTotal;        // Total number of strips in a region
  TH2D*    fEmpty;        // Total number of strips in a region
  TH2D*    fBasic;        // Total number basic hits in a region
  TH2D*    fEmptyVsTotal; // Empty versus total cells 
  TH1D*    fMean;         // Mean calculated by poisson method 
  TH1D*    fOcc;          // Histogram of occupancies 
  TH2D*    fCorr;         // Correction as a function of mean 
  ClassDef(AliPoissonCalculator,3) // Calculate N_ch using Poisson
};

#endif
// Local Variables:
//   mode: C++
// End: