361b291b8ad339aaf85d9f7390af18e502f5b2ad
[u/mrichter/AliRoot.git] / TPC / Base / AliTPCCalROC.h
1 #ifndef ALITPCCALROC_H
2 #define ALITPCCALROC_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4  * See cxx source for full Copyright notice                               */
5
6 /* $Id: AliTPCCalROC.h -1   */
7
8 //////////////////////////////////////////////////
9 //                                              //
10 //  TPC calibration base class for one ROC      //
11 //                                              //
12 //////////////////////////////////////////////////
13 #include <TObject.h>
14 #include <TMath.h>
15 #include <AliTPCROC.h>
16 #include "TLinearFitter.h"
17
18 class TH1F;
19 class TH2F;
20 class TArrayI;
21 class TF1;
22 //_____________________________________________________________________________
23 class AliTPCCalROC : public TNamed {
24
25  public:
26   
27   AliTPCCalROC();
28   AliTPCCalROC(UInt_t sector);
29   AliTPCCalROC(const AliTPCCalROC &c);
30  AliTPCCalROC &operator = (const AliTPCCalROC & param);
31   virtual           ~AliTPCCalROC();  
32   UInt_t        GetSector() const { return fSector;}
33   UInt_t        GetNrows() const               { return fNRows;};
34   UInt_t        GetNchannels()       const     { return fNChannels;};
35   Bool_t        IsInRange(UInt_t row, UInt_t pad) { return (row<fNRows)? (fkIndexes[row]+pad)<fNChannels:kFALSE;}
36   UInt_t        GetNPads(UInt_t row)  const     { return (row<fNRows)? AliTPCROC::Instance()->GetNPads(fSector,row):0;};
37   Float_t      GetValue(UInt_t row, UInt_t pad) const { return ( (row<fNRows) && (fkIndexes[row]+pad)<fNChannels)? fData[fkIndexes[row]+pad]: 0; };
38   Float_t      GetValue(UInt_t channel) const { return  fData[channel]; };
39   void         SetValue(UInt_t row, UInt_t pad, Float_t vd) { if ( row<fNRows && (fkIndexes[row]+pad)<fNChannels)fData[fkIndexes[row]+pad]= vd; }
40   void         SetValue(UInt_t channel, Float_t vd) {fData[channel]= vd; }
41   void         Reset();
42   virtual void Draw(Option_t* option = "");
43   //
44   Bool_t MedianFilter(Int_t deltaRow, Int_t deltaPad, AliTPCCalROC*outlierROC=0, Bool_t doEdge=kTRUE);
45   Bool_t LTMFilter(Int_t deltaRow, Int_t deltaPad, Float_t fraction, Int_t type,  AliTPCCalROC*outlierROC=0, Bool_t doEdge=kTRUE);
46   Bool_t Convolute(Double_t sigmaPad, Double_t sigmaRow,  AliTPCCalROC*outlierPad=0, TF1 *fpad=0, TF1 *frow=0 );
47   //
48   // algebra
49   void Add(Float_t c1); // add c1 to each channel of the ROC
50   void Multiply(Float_t c1); // multiply each channel of the ROC with c1
51   void Add(const AliTPCCalROC * roc, Double_t c1 = 1);  // multiply AliTPCCalROC roc by c1 and add each channel to the coresponing channel in the ROC
52   void Multiply(const AliTPCCalROC * roc);   // multiply each channel of the ROC with the coresponding channel of 'roc'
53   void Divide(const AliTPCCalROC * roc);   // divide each channel of the ROC by the coresponding value of 'roc'
54   // statistic
55   //
56   Double_t GetMean(AliTPCCalROC *const outlierROC = 0) const;
57   Double_t GetRMS(AliTPCCalROC *const outlierROC = 0) const;
58   Double_t GetMedian(AliTPCCalROC *const outlierROC = 0) const;
59   Double_t GetLTM(Double_t *const sigma=0, Double_t fraction=0.9, AliTPCCalROC *const outlierROC = 0);
60   TH1F * MakeHisto1D(Float_t min=4, Float_t max=-4, Int_t type=0);     
61   TH2F * MakeHisto2D(Float_t min=4, Float_t max=-4, Int_t type=0);   
62   TH2F * MakeHistoOutliers(Float_t delta=4, Float_t fraction=0.7, Int_t mode=0);
63
64   AliTPCCalROC * LocalFit(Int_t rowRadius, Int_t padRadius, AliTPCCalROC* ROCoutliers = 0, Bool_t robust = kFALSE, Double_t chi2Threshold = 5, Double_t robustFraction = 0.7);
65   void GlobalFit(const AliTPCCalROC* ROCoutliers, Bool_t robust, TVectorD &fitParam, TMatrixD &covMatrix, Float_t & chi2, Int_t fitType = 1, Double_t chi2Threshold = 5, Double_t robustFraction = 0.7, Double_t err=1);
66   
67   static AliTPCCalROC* CreateGlobalFitCalROC(TVectorD &fitParam, Int_t sector);
68   
69   static void Test();
70  protected:
71   
72   Double_t GetNeighbourhoodValue(TLinearFitter* fitterQ, Int_t row, Int_t pad, Int_t rRadius, Int_t pRadius, AliTPCCalROC *const ROCoutliers, Bool_t robust, Double_t chi2Threshold, Double_t robustFraction);
73   void GetNeighbourhood(TArrayI* &rowArray, TArrayI* &padArray, Int_t row, Int_t pad, Int_t rRadius, Int_t pRadius);
74   
75   UInt_t     fSector;          // sector number
76   UInt_t     fNChannels;       // number of channels
77   UInt_t     fNRows;           // number of rows
78   const UInt_t* fkIndexes;      //!indexes
79   Float_t  *fData;            //[fNChannels] Data
80   ClassDef(AliTPCCalROC,2)    //  TPC ROC calibration class
81
82 };
83
84 #endif