TPC/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,v */
   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 //_____________________________________________________________________________
  22 class AliTPCCalROC : public TObject {
  23
  24  public:
  25
  26   AliTPCCalROC();
  27   AliTPCCalROC(UInt_t sector);
  28   AliTPCCalROC(const AliTPCCalROC &c);
  29  AliTPCCalROC &operator = (const AliTPCCalROC & param);
  30   virtual           ~AliTPCCalROC();
  31   UInt_t        GetSector() const { return fSector;}
  32   UInt_t        GetNrows() const               { return fNRows;};
  33   UInt_t        GetNchannels()       const     { return fNChannels;};
  34   UInt_t        GetNPads(UInt_t row)  const     { return (row<fNRows)? AliTPCROC::Instance()->GetNPads(fSector,row):0;};
  35   Float_t      GetValue(UInt_t row, UInt_t pad) const { return ( (row<fNRows) && (fIndexes[row]+pad)<fNChannels)? fData[fIndexes[row]+pad]: 0; };
  36   Float_t      GetValue(UInt_t channel) const { return  fData[channel]; };
  37   void         SetValue(UInt_t row, UInt_t pad, Float_t vd) { if ( row<fNRows && (fIndexes[row]+pad)<fNChannels)fData[fIndexes[row]+pad]= vd; };
  38   void         SetValue(UInt_t channel, Float_t vd) {fData[channel]= vd; };
  39   virtual void Draw(Option_t* option = "");
  40   //
  41   // algebra
  42   void Add(Float_t c1); // add c1 to each channel of the ROC
  43   void Multiply(Float_t c1); // multiply each channel of the ROC with c1
  44   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
  45   void Multiply(const AliTPCCalROC * roc);   // multiply each channel of the ROC with the coresponding channel of 'roc'
  46   void Divide(const AliTPCCalROC * roc);   // divide each channel of the ROC by the coresponding value of 'roc'
  47   // statistic
  48   //
  49   Double_t GetMean(AliTPCCalROC* outlierROC = 0);
  50   Double_t GetRMS(AliTPCCalROC* outlierROC = 0);
  51   Double_t GetMedian(AliTPCCalROC* outlierROC = 0) ;
  52   Double_t GetLTM(Double_t *sigma=0, Double_t fraction=0.9, AliTPCCalROC* outlierROC = 0);
  53   TH1F * MakeHisto1D(Float_t min=4, Float_t max=-4, Int_t type=0);
  54   TH2F * MakeHisto2D(Float_t min=4, Float_t max=-4, Int_t type=0);
  55   TH2F * MakeHistoOutliers(Float_t delta=4, Float_t fraction=0.7, Int_t mode=0);
  56
  57   AliTPCCalROC * LocalFit(Int_t rowRadius, Int_t padRadius, AliTPCCalROC* ROCoutliers = 0, Bool_t robust = kFALSE, Double_t chi2Threshold = 5, Double_t robustFraction = 0.7);
  58   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);
  59
  60   static AliTPCCalROC* CreateGlobalFitCalROC(TVectorD &fitParam, Int_t sector);
  61
  62   static void Test();
  63  protected:
  64
  65   Double_t GetNeighbourhoodValue(TLinearFitter* fitterQ, Int_t row, Int_t pad, Int_t rRadius, Int_t pRadius, AliTPCCalROC* ROCoutliers, Bool_t robust, Double_t chi2Threshold, Double_t robustFraction);
  66   void GetNeighbourhood(TArrayI* &rowArray, TArrayI* &padArray, Int_t row, Int_t pad, Int_t rRadius, Int_t pRadius);
  67
  68   UInt_t     fSector;          // sector number
  69   UInt_t     fNChannels;       // number of channels
  70   UInt_t     fNRows;           // number of rows
  71   const UInt_t* fIndexes;      //!indexes
  72   Float_t  *fData;            //[fNChannels] Data
  73   ClassDef(AliTPCCalROC,1)    //  TPC ROC calibration class
  74
  75 };
  76
  77 #endif