3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
6 /* $Id: AliTPCCalROC.h -1 */
8 //////////////////////////////////////////////////
10 // TPC calibration base class for one ROC //
12 //////////////////////////////////////////////////
15 #include <AliTPCROC.h>
16 #include "TLinearFitter.h"
22 //_____________________________________________________________________________
23 class AliTPCCalROC : public TNamed {
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; }
42 virtual void Draw(Option_t* option = "");
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 );
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'
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);
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);
67 static AliTPCCalROC* CreateGlobalFitCalROC(TVectorD &fitParam, Int_t sector);
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);
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