ATO-97 - addopt to the new naming covnevention for branches in calibration trees
[u/mrichter/AliRoot.git] / TPC / Base / AliTPCCalPad.h
CommitLineData
07627591 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
8///////////////////////////////////////////////////////////////////////////////
9// //
10// TPC calibration class for parameters which are saved per pad //
11// //
12///////////////////////////////////////////////////////////////////////////////
13
14#include "TNamed.h"
a6d2bd0c 15//#include <TMath.h>
16//#include <AliTPCROC.h>
17#include "TLinearFitter.h"
5312f439 18#include "TVectorD.h"
a6d2bd0c 19//#include <iostream>
20
21
07627591 22
23class AliTPCCalROC;
24class AliTPCCalDet;
184bcc16 25class TObjArray;
26class TGraph;
200be8a6 27class TH2F;
90127643 28class TH1F;
4486a91f 29class TCanvas;
30class TTree;
af6a50bb 31class TH2;
07627591 32
33class AliTPCCalPad : public TNamed {
07627591 34 public:
07627591 35 enum { kNsec = 72 };
07627591 36 AliTPCCalPad();
37 AliTPCCalPad(const Text_t* name, const Text_t* title);
38 AliTPCCalPad(const AliTPCCalPad &c);
184bcc16 39 AliTPCCalPad(TObjArray *arrayROC);
07627591 40 virtual ~AliTPCCalPad();
41 AliTPCCalPad &operator=(const AliTPCCalPad &c);
42 virtual void Copy(TObject &c) const;
a6d2bd0c 43 AliTPCCalROC *GetCalROC(Int_t sector) const {return fROC[sector]; };
4486a91f 44 void SetCalROC(AliTPCCalROC* roc, Int_t sector = -1);
45 virtual void Draw(Option_t* option = "");
7b83b139 46 // TTree functions
47 static AliTPCCalPad *MakePadFromTree(TTree * tree, const char *query, const char*name=0);
48 void AddFriend(TTree * tree, const char *friendName, const char *fname=0);
49 //
50 // convolution
5e24c3ce 51 Bool_t MedianFilter(Int_t deltaRow, Int_t deltaPad, AliTPCCalPad*outlierPad=0, Bool_t doEdge=kTRUE);
52 Bool_t LTMFilter(Int_t deltaRow, Int_t deltaPad, Float_t fraction, Int_t type, AliTPCCalPad*outlierPad=0, Bool_t doEdge=kTRUE);
184bcc16 53 //
90127643 54 // algebra
a6d2bd0c 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'
184bcc16 60 //
a6d2bd0c 61 Double_t GetMeanRMS(Double_t &rms); // Calculates mean and RMS of all ROCs
62 Double_t GetMean(AliTPCCalPad* outlierPad = 0); // return mean of the mean of all ROCs
63 Double_t GetRMS(AliTPCCalPad* outlierPad = 0) ; // return mean of the RMS of all ROCs
64 Double_t GetMedian(AliTPCCalPad* outlierPad = 0) ; // return mean of the median of all ROCs
65 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
90127643 66 TGraph *MakeGraph(Int_t type=0, Float_t ratio=0.7);
67 TH2F *MakeHisto2D(Int_t side=0);
4486a91f 68 TH1F *MakeHisto1D(Float_t min=4, Float_t max=-4, Int_t type=0, Int_t side=0);
a6d2bd0c 69
72d0ab7e 70 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;
6d8681ef 71 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);
5312f439 72
732e90a8 73 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);
74
75 static AliTPCCalPad* CreateCalPadFit(const char* fitFormula, const TVectorD &fitParamSideA, const TVectorD &fitParamSideC);
76
77 static TObjArray *CreateFormulaArray(const char *fitFormula);
78 static void EvalFormulaArray(const TObjArray &arrFitFormulas, TVectorD &results,
79 const Int_t sec, const Int_t row, const Int_t pad);
4486a91f 80 //
81 // default report
82 //
83 static TCanvas * MakeReportPadSector(TTree *chain, const char* varName, const char*varTitle, const char *axisTitle, Float_t min, Float_t max, const char * cutUser="");
84 static TCanvas * MakeReportPadSector2D(TTree *chain, const char* varName, const char*varTitle, const char *axisTitle, Float_t min, Float_t max, const char *cutUser="");
af6a50bb 85 static AliTPCCalPad *MakeCalPadFromHistoRPHI(TH2 * hisA, TH2* hisC);
07627591 86 protected:
87 AliTPCCalROC *fROC[kNsec]; // Array of ROC objects which contain the values per pad
88 ClassDef(AliTPCCalPad,1) // TPC calibration class for parameters which are saved per pad
07627591 89};
90
91#endif