AliRoot Core  v5-06-30 (35d6c57)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AliTPCCalPad.h
Go to the documentation of this file.
1 #ifndef ALITPCCALPAD_H
2 #define ALITPCCALPAD_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 
11 
12 #include "TNamed.h"
13 //#include <TMath.h>
14 //#include <AliTPCROC.h>
15 #include "TLinearFitter.h"
16 #include "TVectorD.h"
17 //#include <iostream>
18 
19 
20 
21 class AliTPCCalROC;
22 class AliTPCCalDet;
23 class TObjArray;
24 class TGraph;
25 class TH2F;
26 class TH1F;
27 class TCanvas;
28 class TTree;
29 class TH2;
30 class TF1;
31 
32 class AliTPCCalPad : public TNamed {
33  public:
34  enum { kNsec = 72 };
35  AliTPCCalPad();
36  AliTPCCalPad(const Text_t* name, const Text_t* title);
37  AliTPCCalPad(const AliTPCCalPad &c);
38  AliTPCCalPad(TObjArray *arrayROC);
39  virtual ~AliTPCCalPad();
41  virtual void Copy(TObject &c) const;
42  AliTPCCalROC *GetCalROC(Int_t sector) const {return fROC[sector]; };
43  void SetCalROC(AliTPCCalROC* roc, Int_t sector = -1);
44  virtual void Draw(Option_t* option = "");
45  // TTree functions
46  static AliTPCCalPad *MakePadFromTree(TTree * tree, const char *query, const char*name=0, Bool_t doFast=kFALSE);
47  void AddFriend(TTree * tree, const char *friendName, const char *fname=0);
48  //
49  // convolution
50  Bool_t MedianFilter(Int_t deltaRow, Int_t deltaPad, AliTPCCalPad*outlierPad=0, Bool_t doEdge=kTRUE);
51  Bool_t LTMFilter(Int_t deltaRow, Int_t deltaPad, Float_t fraction, Int_t type, AliTPCCalPad*outlierPad=0, Bool_t doEdge=kTRUE);
52  Bool_t Convolute(Double_t sigmaPad, Double_t sigmaRow, AliTPCCalPad*outlierPad=0, TF1 *fpad=0, TF1 *frow=0 );
53  //
54  // algebra
55  void Add(Float_t c1); // add constant c1 to all channels of all ROCs
56  void Multiply(Float_t c1); // multiply each channel of all ROCs with c1
57  void Add(const AliTPCCalPad * roc, Double_t c1 = 1); // multiply AliTPCCalPad 'pad' by c1 and add each channel to the coresponing channel in all ROCs
58  void Multiply(const AliTPCCalPad * pad); // multiply each channel of all ROCs with the coresponding channel of 'pad'
59  void Divide(const AliTPCCalPad * pad); // divide each channel of all ROCs by the coresponding channel of 'pad'
60  void Reset();
61  //
62  Double_t GetMeanRMS(Double_t &rms); // Calculates mean and RMS of all ROCs
63  Double_t GetMean(AliTPCCalPad* outlierPad = 0); // return mean of the mean of all ROCs
64  Double_t GetRMS(AliTPCCalPad* outlierPad = 0) ; // return mean of the RMS of all ROCs
65  Double_t GetMedian(AliTPCCalPad* outlierPad = 0) ; // return mean of the median of all ROCs
66  Double_t GetLTM(Double_t *sigma=0, Double_t fraction=0.9, AliTPCCalPad* outlierPad = 0); // return mean of the LTM and sigma of all ROCs
67  TGraph *MakeGraph(Int_t type=0, Float_t ratio=0.7);
68  TH2F *MakeHisto2D(Int_t side=0);
69  TH1F *MakeHisto1D(Float_t min=4, Float_t max=-4, Int_t type=0, Int_t side=0);
70 
71  AliTPCCalPad* LocalFit(const char* padName, Int_t rowRadius, Int_t padRadius, AliTPCCalPad* Padoutliers = 0, Bool_t robust = kFALSE, Double_t chi2Threshold = 5, Double_t robustFraction = 0.7, Bool_t printCurrentSector = kFALSE) const;
72  AliTPCCalPad* GlobalFit(const char* padName, AliTPCCalPad* Padoutliers = 0, Bool_t robust = kFALSE, Int_t fitType = 1, Double_t chi2Threshold = 5, Double_t robustFraction = 0.7, Double_t err=1, TObjArray *fitParArr=0x0, TObjArray *fitCovArr=0x0);
73 
74  void GlobalSidesFit(const AliTPCCalPad* PadOutliers, const char* fitFormula, TVectorD &fitParamSideA, TVectorD &fitParamSideC, TMatrixD &covMatrixSideA, TMatrixD &covMatrixSideC, Float_t &chi2SideA, Float_t &chi2SideC, AliTPCCalPad *pointError=0, Bool_t robust = kFALSE, Double_t robustFraction = 0.7);
75 
76  static AliTPCCalPad* CreateCalPadFit(const char* fitFormula, const TVectorD &fitParamSideA, const TVectorD &fitParamSideC);
77 
78  static TObjArray *CreateFormulaArray(const char *fitFormula);
79  static void EvalFormulaArray(const TObjArray &arrFitFormulas, TVectorD &results,
80  const Int_t sec, const Int_t row, const Int_t pad);
81  //
82  // default report
83  //
84  static TCanvas * MakeReportPadSector(TTree *chain, const char* varName, const char*varTitle, const char *axisTitle, Float_t min, Float_t max, const char * cutUser="");
85  static TCanvas * MakeReportPadSector2D(TTree *chain, const char* varName, const char*varTitle, const char *axisTitle, Float_t min, Float_t max, const char *cutUser="");
86  static AliTPCCalPad *MakeCalPadFromHistoRPHI(TH2 * hisA, TH2* hisC);
87  protected:
89  ClassDef(AliTPCCalPad,1) // TPC calibration class for parameters which are saved per pad
92 };
93 
94 #endif
static TCanvas * MakeReportPadSector(TTree *chain, const char *varName, const char *varTitle, const char *axisTitle, Float_t min, Float_t max, const char *cutUser="")
Double_t GetMean(AliTPCCalPad *outlierPad=0)
AliTPCCalPad * GlobalFit(const char *padName, AliTPCCalPad *Padoutliers=0, Bool_t robust=kFALSE, Int_t fitType=1, Double_t chi2Threshold=5, Double_t robustFraction=0.7, Double_t err=1, TObjArray *fitParArr=0x0, TObjArray *fitCovArr=0x0)
AliTPCCalROC * GetCalROC(Int_t sector) const
Definition: AliTPCCalPad.h:42
virtual void Draw(Option_t *option="")
Double_t GetMedian(AliTPCCalPad *outlierPad=0)
#define TObjArray
TGraph * MakeGraph(Int_t type=0, Float_t ratio=0.7)
void SetCalROC(AliTPCCalROC *roc, Int_t sector=-1)
static void EvalFormulaArray(const TObjArray &arrFitFormulas, TVectorD &results, const Int_t sec, const Int_t row, const Int_t pad)
TCut cutUser
void Divide(const AliTPCCalPad *pad)
TChain * chain
Double_t GetRMS(AliTPCCalPad *outlierPad=0)
TTree * tree
static TCanvas * MakeReportPadSector2D(TTree *chain, const char *varName, const char *varTitle, const char *axisTitle, Float_t min, Float_t max, const char *cutUser="")
static AliTPCCalPad * MakeCalPadFromHistoRPHI(TH2 *hisA, TH2 *hisC)
static AliTPCCalPad * CreateCalPadFit(const char *fitFormula, const TVectorD &fitParamSideA, const TVectorD &fitParamSideC)
static AliTPCCalPad * MakePadFromTree(TTree *tree, const char *query, const char *name=0, Bool_t doFast=kFALSE)
Double_t GetLTM(Double_t *sigma=0, Double_t fraction=0.9, AliTPCCalPad *outlierPad=0)
AliTPCCalPad * LocalFit(const char *padName, Int_t rowRadius, Int_t padRadius, AliTPCCalPad *Padoutliers=0, Bool_t robust=kFALSE, Double_t chi2Threshold=5, Double_t robustFraction=0.7, Bool_t printCurrentSector=kFALSE) const
TH1F * MakeHisto1D(Float_t min=4, Float_t max=-4, Int_t type=0, Int_t side=0)
Bool_t MedianFilter(Int_t deltaRow, Int_t deltaPad, AliTPCCalPad *outlierPad=0, Bool_t doEdge=kTRUE)
static TObjArray * CreateFormulaArray(const char *fitFormula)
void Add(Float_t c1)
Double_t GetMeanRMS(Double_t &rms)
TPC calibration base class for one ROC.
Definition: AliTPCCalROC.h:19
TH2F * MakeHisto2D(Int_t side=0)
Bool_t LTMFilter(Int_t deltaRow, Int_t deltaPad, Float_t fraction, Int_t type, AliTPCCalPad *outlierPad=0, Bool_t doEdge=kTRUE)
Bool_t Convolute(Double_t sigmaPad, Double_t sigmaRow, AliTPCCalPad *outlierPad=0, TF1 *fpad=0, TF1 *frow=0)
AliTPCCalPad & operator=(const AliTPCCalPad &c)
virtual void Copy(TObject &c) const
void GlobalSidesFit(const AliTPCCalPad *PadOutliers, const char *fitFormula, TVectorD &fitParamSideA, TVectorD &fitParamSideC, TMatrixD &covMatrixSideA, TMatrixD &covMatrixSideC, Float_t &chi2SideA, Float_t &chi2SideC, AliTPCCalPad *pointError=0, Bool_t robust=kFALSE, Double_t robustFraction=0.7)
virtual ~AliTPCCalPad()
void AddFriend(TTree *tree, const char *friendName, const char *fname=0)
char * fname
void Multiply(Float_t c1)
AliTPCCalROC * fROC[kNsec]
Definition: AliTPCCalPad.h:88