]> git.uio.no Git - u/mrichter/AliRoot.git/blob - TPC/AliTPCcalibDButil.h
Missing protection - check the presence of array
[u/mrichter/AliRoot.git] / TPC / AliTPCcalibDButil.h
1 #ifndef AliTPCcalibDButil_H
2 #define AliTPCcalibDButil_H
3
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5  * See cxx source for full Copyright notice                               */
6
7 ///////////////////////////////////////////////////////////////////////////////
8 //                                                                           //
9 // Class providing the calculation of derived quantities (mean,rms,fits,...) //
10 //       of calibration entries                                              //
11 //                                                                           //
12 ///////////////////////////////////////////////////////////////////////////////
13
14 #include <TObject.h>
15
16 class AliDCSSensorArray;
17 class AliTPCcalibDB;
18 class AliTPCCalPad;
19 class AliTPCmapper;
20 class AliTPCCalibRaw;
21
22 class AliTPCcalibDButil : public TObject
23 {
24 public:
25   AliTPCcalibDButil();
26   virtual ~AliTPCcalibDButil();
27
28   void UpdateFromCalibDB();
29   //data processing functions
30   void ProcessCEdata(const char* fitFormula, TVectorD &fitResultsA, TVectorD &fitResultsC,
31                      Int_t &noutliersCE, Double_t & chi2A, Double_t &chi2C, AliTPCCalPad *outCE=0);
32   void ProcessCEgraphs(TVectorD &vecTEntries, TVectorD &vecTMean, TVectorD &vecTRMS, TVectorD &vecTMedian,
33                        TVectorD &vecQEntries, TVectorD &vecQMean, TVectorD &vecQRMS, TVectorD &vecQMedian,
34                        Float_t &driftTimeA, Float_t &driftTimeC );
35   void ProcessNoiseData(TVectorD &vNoiseMean, TVectorD &vNoiseMeanSenRegions,
36                         TVectorD &vNoiseRMS, TVectorD &vNoiseRMSSenRegions,
37                         Int_t &nonMaskedZero);
38   void ProcessPulser(TVectorD &vMeanTime);
39   void ProcessALTROConfig(Int_t &nMasked);
40   void ProcessGoofie(TVectorD & vecEntries, TVectorD & vecMedian, TVectorD &vecMean, TVectorD &vecRMS);
41   //processing functions using reference data
42   void ProcessPedestalVariations(TVectorF &pedestalDeviations);
43   void ProcessNoiseVariations(TVectorF &noiseDeviations);
44   void ProcessPulserVariations(TVectorF &pulserQdeviations, Float_t &varQMean, Int_t &npadsOutOneTB, Int_t &npadsOffAdd);
45   //getter preprocess information
46   Int_t GetNPulserOutliers() const {return fNpulserOutliers;}
47   Float_t GetMeanAltro(const AliTPCCalROC *roc, const Int_t row, const Int_t pad, AliTPCCalROC *rocOut=0x0);
48   AliTPCCalPad *GetPulserOutlierMap() const {return fPulserOutlier;}
49   //setters for pad by pad information
50   void SetPulserData(AliTPCCalPad *tmean, AliTPCCalPad *trms=0x0, AliTPCCalPad *qmean=0x0)
51                 {fPulserTmean=tmean; fPulserTrms=trms; fPulserQmean=qmean;}
52   void SetCEData(AliTPCCalPad *tmean, AliTPCCalPad *trms=0x0, AliTPCCalPad *qmean=0x0)
53                 {fCETmean=tmean; fCETrms=trms; fCEQmean=qmean;}
54   void SetNoisePedestal(AliTPCCalPad *noise, AliTPCCalPad *pedestal=0x0)
55                 {fPadNoise=noise; fPedestals=pedestal;}
56   void SetALTROData(AliTPCCalPad *masked)
57                 {fALTROMasked=masked;}
58   void SetGoofieArray(AliDCSSensorArray *arr) {fGoofieArray=arr;}
59   //setters for pad by pad information
60   void SetRefFile(const char* filename);
61   void SetRefPulserData(AliTPCCalPad *tmean, AliTPCCalPad *trms=0x0, AliTPCCalPad *qmean=0x0)
62                 {fRefPulserTmean=tmean; fRefPulserTrms=trms; fRefPulserQmean=qmean;}
63   void SetRefCEData(AliTPCCalPad *tmean, AliTPCCalPad *trms=0x0, AliTPCCalPad *qmean=0x0)
64                 {fRefCETmean=tmean; fRefCETrms=trms; fRefCEQmean=qmean;}
65   void SetRefNoisePedestal(AliTPCCalPad *noise, AliTPCCalPad *pedestal=0x0)
66                 {fRefPadNoise=noise; fRefPedestals=pedestal;}
67   void SetRefALTROData(AliTPCCalPad *masked)
68                 {fRefALTROMasked=masked;}
69   
70   //creation of derived pad by pad calibration data
71   AliTPCCalPad *CreatePadTime0(Int_t model, Double_t &gyA, Double_t &gyC, Double_t &chi2A, Double_t &chi2C);
72   //
73   // create outlyer maps
74   //
75   AliTPCCalPad *CreateCEOutlyerMap(Int_t &noutliersCE, AliTPCCalPad *ceOut=0, Float_t minSignal=10, Float_t cutTrmsMin=0.9, Float_t cutTrmsMax=1.2, Float_t cutMaxDistT=0.7);
76   AliTPCCalPad *CreatePulserOutlyerMap(Int_t &noutliersPulser, AliTPCCalPad *pulserOut=0, Float_t cutTime=3, Float_t cutnRMSQ=5, Float_t cutnRMSrms=5);
77   //
78   AliTPCCalPad *CreatePadTime0CE(TVectorD &fitResultsA, TVectorD&fitResultsC, Int_t &nOut, Double_t &chi2A, Double_t &chi2C, const char *dumpfile=0);
79   //
80
81   void UpdatePulserOutlierMap();
82   void UpdateRefPulserOutlierMap();
83   void PulserOutlierMap(AliTPCCalPad *pulOut, const AliTPCCalPad *pulT, const AliTPCCalPad *pulQ);
84 private:
85   AliTPCcalibDB *fCalibDB;            //pointer to calibDB object
86   AliTPCCalPad  *fPadNoise;           //noise information
87   AliTPCCalPad  *fPedestals;          //pedestal information
88   AliTPCCalPad  *fPulserTmean;        //pulser mean time information 
89   AliTPCCalPad  *fPulserTrms;         //pulser rms time information
90   AliTPCCalPad  *fPulserQmean;        //pulser mean q information
91   AliTPCCalPad  *fPulserOutlier;      //pulser outlier map
92   AliTPCCalPad  *fCETmean;            //central electrode mean time information
93   AliTPCCalPad  *fCETrms;             //central electrode rms time information
94   AliTPCCalPad  *fCEQmean;            //central electrode mean q information
95   AliTPCCalPad  *fALTROMasked;        //ALTRO masked channels information
96   //
97   AliTPCCalibRaw *fCalibRaw;          //raw calibration object
98   //reference data
99   AliTPCCalPad  *fRefPadNoise;           //Reference noise information
100   AliTPCCalPad  *fRefPedestals;          //Reference pedestal information
101   AliTPCCalPad  *fRefPulserTmean;        //Reference pulser mean time information
102   AliTPCCalPad  *fRefPulserTrms;         //Reference pulser rms time information
103   AliTPCCalPad  *fRefPulserQmean;        //Reference pulser mean q information
104   AliTPCCalPad  *fRefPulserOutlier;      //Reference pulser outlier map
105   AliTPCCalPad  *fRefCETmean;            //Reference central electrode mean time information
106   AliTPCCalPad  *fRefCETrms;             //Reference central electrode rms time information
107   AliTPCCalPad  *fRefCEQmean;            //Reference central electrode mean q information
108   AliTPCCalPad  *fRefALTROMasked;        //Reference ALTRO masked channels information
109   //
110   AliTPCCalibRaw *fRefCalibRaw;          //Reference raw calibration object
111   
112   //
113   AliDCSSensorArray* fGoofieArray;    //Goofie Data
114   //
115   AliTPCmapper  *fMapper;             //TPC mapping handler
116   Int_t fNpulserOutliers;             //number of outliers from Pulser calibration
117   
118   Float_t fIrocTimeOffset;               //timing offset between IROC and OROC in timebins
119   Float_t fCETmaxLimitAbs;               //maximum variation in CE data before pads will be treated as outliers
120   Float_t fPulTmaxLimitAbs;              //maximum variation of Pulser Signals (time) before pads will be treated as outliers
121   Float_t fPulQmaxLimitAbs;              //maximum variation of Pulser Signals (charge) before pads will be treated as outliers
122   Float_t fPulQminLimit;                 //minimum charge value for Pulser Signals before pads will be treated as outliers
123   
124   AliTPCcalibDButil (const AliTPCcalibDButil& );
125   AliTPCcalibDButil& operator= (const AliTPCcalibDButil& );
126
127     
128   ClassDef(AliTPCcalibDButil,0)
129 };
130
131
132 #endif