New macro "MUONreco.C" for using the event reconstruction package in C++
[u/mrichter/AliRoot.git] / TPC / AliH2F.h
CommitLineData
cc80f89e 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
13class TMatrix;
14class TH1F;
15class TClonesArray;
16
17class AliH2F : public TH2F {
18public:
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
24public:
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
49public:
50
51protected:
52
53private:
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 */