TPC/AliH2F.h

   1 #ifndef TH2FSMOOTH_H
   2 #define TH2FSMOOTH_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 /* $Id$ */
   7
   8 // include files and class forward declarations
   9
  10 #include "TH2.h"
  11 #include "TMatrix.h"
  12
  13 class TMatrix;
  14 class TH1F;
  15 class TClonesArray;
  16
  17 class AliH2F : public TH2F {
  18 public:
  19   AliH2F();
  20   AliH2F(const Text_t *name,const Text_t *title,Int_t nbinsx,
  21              Axis_t xlow,Axis_t xup,Int_t nbinsy,Axis_t ylow,Axis_t yup);
  22     ~AliH2F();
  23
  24 public:
  25   AliH2F(const AliH2F &);
  26   AliH2F & operator = (const AliH2F &);
  27   TClonesArray * FindPeaks(Float_t threshold, Float_t noise);
  28     //find peaks and write it in form of AliTPCcluster to array
  29   void ClearSpectrum();
  30   void AddGauss(Float_t x,Float_t y,Float_t sx, Float_t sy, Float_t max);
  31   void AddNoise(Float_t sn);
  32   void ClearUnderTh(Int_t threshold);
  33   void Round();
  34   //round float values to integer values
  35   //  void Smooth();
  36   AliH2F * GetSubrange2d(Float_t xmin, Float_t xmax,
  37                              Float_t ymin, Float_t ymax);
  38   //create new  2D histogram
  39   Float_t  GetOccupancy(Float_t th=1. , Float_t xmin=0, Float_t xmax=0,
  40                              Float_t ymin=0, Float_t ymax=0);
  41   //calculate ration of channel over threshold to all channels
  42   TH1F * GetAmplitudes(Float_t zmin, Float_t zmax, Float_t th=1. , Float_t xmin=0, Float_t xmax=0,
  43                              Float_t ymin=0, Float_t ymax=0);
  44                        //generate one dim histogram of amplitudes
  45   void SetSmoothSigma(Float_t sigmaX, Float_t sigmaY);
  46   void SetSmoothRange(Float_t dx,Float_t dy);
  47   TMatrix *  GetFitMatrix(Int_t irow, Int_t icolumn);
  48
  49 public:
  50
  51 protected:
  52
  53 private:
  54   void SmoothCell(Int_t irow, Int_t icolumn) ;
  55   Float_t fSigmaX2;   //sigma x  used for weight calculation in smooth proces
  56   Float_t fSigmaY2;   //sigma x  used for weight calculation in smooth proces
  57   Int_t fdx;        //dx range on which we are smoothing around 1 point
  58   Int_t fdy;        //dy  range on which we are smoothing around 1 point
  59   TMatrix fFitMatrix; //matrix with fitted parameters
  60   ClassDef(AliH2F,1)
  61 };
  62
  63 #endif /*TH2FSMOOTH_H */