]> git.uio.no Git - u/mrichter/AliRoot.git/blame - TPC/TPCbase/AliTPCcalibDButil.h
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
doxy: TPC/TPCbase converted
[u/mrichter/AliRoot.git] / TPC / TPCbase / AliTPCcalibDButil.h
CommitLineData
8166d5d7 1#ifndef ALITPCCALIBDBUTIL_H
2#define ALITPCCALIBDBUTIL_H
892226be 3
4/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
6
7d855b04 7/// \class AliTPCcalibDButil
8/// \brief Class providing the calculation of derived quantities (mean,rms,fits,...) of calibration entries
892226be 9
10#include <TObject.h>
ea09e1dd 11#include <TArrayI.h>
1be09c9e 12#include <TVectorDfwd.h>
13#include <TVectorFfwd.h>
14#include <TMatrixDfwd.h>
892226be 15
949d8707 16class TGraph;
17class TMap;
892226be 18class AliDCSSensorArray;
19class AliTPCcalibDB;
20class AliTPCCalPad;
37f2177e 21class AliTPCCalROC;
892226be 22class AliTPCmapper;
6e7d7dc4 23class AliTPCCalibRaw;
949d8707 24class AliCDBEntry;
37f2177e 25class AliDCSSensor;
26class AliDCSSensorArray;
27class AliTPCSensorTempArray;
949d8707 28class AliTPCdataQA;
a23ba1c3 29class TGraphErrors;
30class TTreeSRedirector;
37f2177e 31class AliTPCCalROC;
a3b590cf 32class AliSplineFit;
8166d5d7 33
892226be 34class AliTPCcalibDButil : public TObject
35{
36public:
37 AliTPCcalibDButil();
38 virtual ~AliTPCcalibDButil();
39
40 void UpdateFromCalibDB();
41 //data processing functions
7390f655 42 void ProcessCEdata(const char* fitFormula, TVectorD &fitResultsA, TVectorD &fitResultsC,
8166d5d7 43 Int_t &noutliersCE, Double_t & chi2A, Double_t &chi2C, AliTPCCalPad * const outCE=0);
892226be 44 void ProcessCEgraphs(TVectorD &vecTEntries, TVectorD &vecTMean, TVectorD &vecTRMS, TVectorD &vecTMedian,
45 TVectorD &vecQEntries, TVectorD &vecQMean, TVectorD &vecQRMS, TVectorD &vecQMedian,
46 Float_t &driftTimeA, Float_t &driftTimeC );
47 void ProcessNoiseData(TVectorD &vNoiseMean, TVectorD &vNoiseMeanSenRegions,
48 TVectorD &vNoiseRMS, TVectorD &vNoiseRMSSenRegions,
93425263 49 Int_t &nonMaskedZero, Int_t &nNaN);
abb20887 50 void ProcessQAData(TVectorD &vQaOcc, TVectorD &vQaQtot, TVectorD &vQaQmax);
892226be 51 void ProcessPulser(TVectorD &vMeanTime);
52 void ProcessALTROConfig(Int_t &nMasked);
53 void ProcessGoofie(TVectorD & vecEntries, TVectorD & vecMedian, TVectorD &vecMean, TVectorD &vecRMS);
949d8707 54
7390f655 55 //processing functions using reference data
56 void ProcessPedestalVariations(TVectorF &pedestalDeviations);
57 void ProcessNoiseVariations(TVectorF &noiseDeviations);
58 void ProcessPulserVariations(TVectorF &pulserQdeviations, Float_t &varQMean, Int_t &npadsOutOneTB, Int_t &npadsOffAdd);
949d8707 59
892226be 60 //getter preprocess information
61 Int_t GetNPulserOutliers() const {return fNpulserOutliers;}
8166d5d7 62 Float_t GetMeanAltro(const AliTPCCalROC *roc, const Int_t row, const Int_t pad, AliTPCCalROC * const rocOut=0x0);
892226be 63 AliTPCCalPad *GetPulserOutlierMap() const {return fPulserOutlier;}
949d8707 64
65 //getters ref data
66 TMap *GetReferenceMap() const {return fRefMap;}
c95e0219 67 Int_t GetReferenceRun(const char* type) const;
949d8707 68 const char* GetRefValidity() const {return fRefValidity.Data();}
69
70 AliTPCCalPad* GetRefPadNoise() const {return fRefPadNoise;}
71 AliTPCCalPad* GetRefPedestals() const {return fRefPedestals;}
72 AliTPCCalPad* GetRefPedestalMasked() const {return fRefPedestalMasked;}
73 AliTPCCalPad* GetRefPulserTmean() const {return fRefPulserTmean;}
74 AliTPCCalPad* GetRefPulserTrms() const {return fRefPulserTrms;}
75 AliTPCCalPad* GetRefPulserQmean() const {return fRefPulserQmean;}
76 AliTPCCalPad* GetRefPulserOutlier() const {return fRefPulserOutlier;}
77 AliTPCCalPad* GetRefPulserMasked() const {return fRefPulserMasked;}
78 AliTPCCalPad* GetRefCETmean() const {return fRefCETmean;}
79 AliTPCCalPad* GetRefCETrms() const {return fRefCETrms;}
80 AliTPCCalPad* GetRefCEQmean() const {return fRefCEQmean;}
81 AliTPCCalPad* GetRefCEMasked() const {return fRefCEMasked;}
82 AliTPCCalPad* GetRefALTROFPED() const {return fRefALTROFPED;}
83 AliTPCCalPad* GetRefALTROZsThr() const {return fRefALTROZsThr;}
84 AliTPCCalPad* GetRefALTROAcqStart() const {return fRefALTROAcqStart;}
85 AliTPCCalPad* GetRefALTROAcqStop() const {return fRefALTROAcqStop;}
86 AliTPCCalPad* GetRefALTROMasked() const {return fRefALTROMasked;}
87
88
892226be 89 //setters for pad by pad information
8166d5d7 90 void SetPulserData(AliTPCCalPad * const tmean, AliTPCCalPad * const trms=0x0, AliTPCCalPad * const qmean=0x0)
892226be 91 {fPulserTmean=tmean; fPulserTrms=trms; fPulserQmean=qmean;}
98a4cc77 92 void SetCEData(AliTPCCalPad *const tmean, AliTPCCalPad *const trms=0x0, AliTPCCalPad *const qmean=0x0)
892226be 93 {fCETmean=tmean; fCETrms=trms; fCEQmean=qmean;}
98a4cc77 94 void SetNoisePedestal(AliTPCCalPad *const noise, AliTPCCalPad *const pedestal=0x0)
892226be 95 {fPadNoise=noise; fPedestals=pedestal;}
98a4cc77 96 void SetALTROData(AliTPCCalPad *const masked)
892226be 97 {fALTROMasked=masked;}
98a4cc77 98 void SetGoofieArray(AliDCSSensorArray *const arr) {fGoofieArray=arr;}
949d8707 99
7390f655 100 //setters for pad by pad information
101 void SetRefFile(const char* filename);
949d8707 102 void SetReferenceRun(Int_t run=-1);
103 void UpdateRefDataFromOCDB();
98a4cc77 104 void SetRefPulserData(AliTPCCalPad *const tmean, AliTPCCalPad *const trms=0x0, AliTPCCalPad *const qmean=0x0)
7390f655 105 {fRefPulserTmean=tmean; fRefPulserTrms=trms; fRefPulserQmean=qmean;}
98a4cc77 106 void SetRefCEData(AliTPCCalPad *const tmean, AliTPCCalPad *const trms=0x0, AliTPCCalPad *const qmean=0x0)
7390f655 107 {fRefCETmean=tmean; fRefCETrms=trms; fRefCEQmean=qmean;}
98a4cc77 108 void SetRefNoisePedestal(AliTPCCalPad *const noise, AliTPCCalPad *const pedestal=0x0)
7390f655 109 {fRefPadNoise=noise; fRefPedestals=pedestal;}
98a4cc77 110 void SetRefALTROData(AliTPCCalPad *const masked)
7390f655 111 {fRefALTROMasked=masked;}
112
892226be 113 //creation of derived pad by pad calibration data
2cb269df 114 AliTPCCalPad *CreatePadTime0(Int_t model, Double_t &gyA, Double_t &gyC, Double_t &chi2A, Double_t &chi2C);
892226be 115 //
2cb269df 116 // create outlyer maps
117 //
8166d5d7 118 AliTPCCalPad *CreateCEOutlyerMap(Int_t &noutliersCE, AliTPCCalPad * const ceOut=0, Float_t minSignal=10, Float_t cutTrmsMin=0.9, Float_t cutTrmsMax=1.2, Float_t cutMaxDistT=0.7);
119 AliTPCCalPad *CreatePulserOutlyerMap(Int_t &noutliersPulser, AliTPCCalPad * const pulserOut=0, Float_t cutTime=3, Float_t cutnRMSQ=5, Float_t cutnRMSrms=5);
2cb269df 120 //
121 AliTPCCalPad *CreatePadTime0CE(TVectorD &fitResultsA, TVectorD&fitResultsC, Int_t &nOut, Double_t &chi2A, Double_t &chi2C, const char *dumpfile=0);
122 //
123
892226be 124 void UpdatePulserOutlierMap();
7390f655 125 void UpdateRefPulserOutlierMap();
126 void PulserOutlierMap(AliTPCCalPad *pulOut, const AliTPCCalPad *pulT, const AliTPCCalPad *pulQ);
817766d5 127
949d8707 128 const char* GetGUIRefTreeDefaultName();
129
130 Bool_t CreateGUIRefTree(const char* filename="");
817766d5 131 //
132 // graph tools
133 //
a23ba1c3 134 static Double_t GetLaserTime0(Int_t run, Int_t timeStamp, Int_t deltaT, Int_t side);
8166d5d7 135 static TGraph* FilterGraphMedian(TGraph * const graph, Float_t sigmaCut, Double_t &medianY);
a23ba1c3 136 static TGraph* FilterGraphMedianAbs(TGraph * graph, Float_t cut, Double_t &medianY);
137 static TGraphErrors* FilterGraphMedianErr(TGraphErrors * graph, Float_t sigmaCut,Double_t &medianY);
138 //
139 static void Sort(TGraph *graph);
140 static void SmoothGraph(TGraph *graph, Double_t delta);
817766d5 141 static Int_t GetNearest(TGraph *graph, Double_t xref, Double_t &dx, Double_t &y);
8166d5d7 142 static Double_t EvalGraphConst(TGraph * const graph, Double_t xref);
108a67db 143 static Double_t EvalGraphConst(AliSplineFit *graph, Double_t xref);
144
1e722a63 145 //
146 // Filter sensors
147 //
148 static Float_t FilterSensor(AliDCSSensor * sensor, Double_t ymin, Double_t ymax, Double_t maxdy, Double_t sigmaCut);
a980538f 149 //
150 // Filter AliRelAlignmentKalman - Alignment/Drift velocity
151 //
8166d5d7 152 static TMatrixD* MakeStatRelKalman(TObjArray * const array, Float_t minFraction, Int_t minStat, Float_t maxvd);
153 static TObjArray *SmoothRelKalman(TObjArray * const array,const TMatrixD & stat, Bool_t direction, Float_t sigmaCut);
154 static TObjArray *SmoothRelKalman(TObjArray * const arrayP, TObjArray * const arrayM);
155 static void FilterCE(Double_t deltaT=100, Double_t cutAbs=10, Double_t cutSigma=4., TTreeSRedirector * const pcstream=0);
156 static void FilterTracks(Int_t run, Double_t cutSigma=20., TTreeSRedirector * const pcstream=0);
1e722a63 157 static Float_t FilterTemperature(AliTPCSensorTempArray *tempArray, Double_t ymin=15, Double_t ymax=22, Double_t sigmaCut=5);
a23ba1c3 158
8166d5d7 159 static void FilterGoofie(AliDCSSensorArray * goofieArray, Double_t deltaT=2, Double_t cutSigma=4., Double_t minVdn=8.5, Double_t maxVdn=9.05, TTreeSRedirector * const pcstream=0);
a23ba1c3 160 static Double_t GetTriggerOffsetTPC(Int_t run, Int_t timeStamp, Double_t deltaT=86400, Double_t deltaTLaser=3600, Int_t valType=0);
161 static Double_t GetVDriftTPC(Double_t &dist, Int_t run, Int_t timeStamp, Double_t deltaT=86400, Double_t deltaTLaser=3600, Int_t valType=0);
162 static Double_t GetVDriftTPCLaserTracks(Double_t &dist,Int_t run, Int_t timeStamp, Double_t deltaT=43200, Int_t side=2);
78f17711 163 static Double_t GetVDriftTPCLaserTracksOnline(Double_t &dist, Int_t /*run*/, Int_t timeStamp, Double_t deltaT=43200, Int_t side=2);
a23ba1c3 164 static Double_t GetVDriftTPCCE(Double_t &dist, Int_t run, Int_t timeStamp, Double_t deltaT=43200, Int_t side=2);
cc65e4f5 165 static Double_t GetVDriftTPCITS(Double_t &dist, Int_t run, Int_t timeStamp);
166 static Double_t GetTime0TPCITS(Double_t &dist, Int_t run, Int_t timeStamp);
a23ba1c3 167 Int_t MakeRunList(Int_t startRun, Int_t stopRun); // find the list of usable runs
168 Int_t FindRunTPC(Int_t itime, Bool_t debug=kFALSE);
b5db1171 169
170 AliTPCCalPad* GetRefCalPad(AliCDBEntry *entry, const char* objName);
171 AliTPCCalPad* GetRefCalPad(AliCDBEntry *entry);
172 AliTPCCalPad* GetAltroMasked(const char* cdbPath, const char* name);
173 Bool_t HasRefChanged(const char *cdbPath);
174 Int_t GetCurrentReferenceRun(const char* type) const;
175 AliCDBEntry* GetRefEntry(const char* cdbPath);
c4e4498c 176
177 static TTree* ConnectGainTrees(TString baseDir);
06a262f6 178 static TTree* ConnectPulserTrees(TString baseDir, TTree * tree=0);
179 //static TTree* ConnectNoiseTrees(TString baseDir, TTree * tree=0);
525f22d2 180 static TTree* ConnectDistortionTrees(TString baseDir, TString selection, TTree *tMain);
085597e9 181 static TTree* ConnectCalPadTrees(TString baseDir, TString pattern, TTree *tMain, Bool_t checkAliases);
c4f17a0d 182
892226be 183private:
7d855b04 184 AliTPCcalibDB *fCalibDB; ///< pointer to calibDB object
185 AliTPCCalPad *fPadNoise; ///< noise information
186 AliTPCCalPad *fPedestals; ///< pedestal information
187 AliTPCCalPad *fPulserTmean; ///< pulser mean time information
188 AliTPCCalPad *fPulserTrms; ///< pulser rms time information
189 AliTPCCalPad *fPulserQmean; ///< pulser mean q information
190 AliTPCCalPad *fPulserOutlier; ///< pulser outlier map
191 AliTPCCalPad *fCETmean; ///< central electrode mean time information
192 AliTPCCalPad *fCETrms; ///< central electrode rms time information
193 AliTPCCalPad *fCEQmean; ///< central electrode mean q information
194 AliTPCCalPad *fALTROMasked; ///< ALTRO masked channels information
892226be 195 //
7d855b04 196 AliTPCCalibRaw *fCalibRaw; ///< raw calibration object
949d8707 197 //
7d855b04 198 AliTPCdataQA *fDataQA; ///< data qa
7390f655 199 //reference data
7d855b04 200 TMap *fRefMap; ///< latest map to reference information
201 TMap *fCurrentRefMap; ///< reference data map of entries currently loaded
202 TString fRefValidity; ///< validity range of reference data
949d8707 203 //
7d855b04 204 AliTPCCalPad *fRefPadNoise; ///< Reference noise information
205 AliTPCCalPad *fRefPedestals; ///< Reference pedestal information
206 AliTPCCalPad *fRefPedestalMasked; ///< Reference masked channels in pedestal run
207 AliTPCCalPad *fRefPulserTmean; ///< Reference pulser mean time information
208 AliTPCCalPad *fRefPulserTrms; ///< Reference pulser rms time information
209 AliTPCCalPad *fRefPulserQmean; ///< Reference pulser mean q information
210 AliTPCCalPad *fRefPulserOutlier; ///< Reference pulser outlier map
211 AliTPCCalPad *fRefPulserMasked; ///< Reference masked channels in pulser run
212 AliTPCCalPad *fRefCETmean; ///< Reference central electrode mean time information
213 AliTPCCalPad *fRefCETrms; ///< Reference central electrode rms time information
214 AliTPCCalPad *fRefCEQmean; ///< Reference central electrode mean q information
215 AliTPCCalPad *fRefCEMasked; ///< Reference masked channels in laser run
216 AliTPCCalPad *fRefALTROFPED; ///< Reference fixed pedestal value
217 AliTPCCalPad *fRefALTROZsThr; ///< Reference zero suppression threshol
218 AliTPCCalPad *fRefALTROAcqStart; ///< Reference accquistion start time bin
219 AliTPCCalPad *fRefALTROAcqStop; ///< Reference accquistion stop time bin
220 AliTPCCalPad *fRefALTROMasked; ///< Reference ALTRO masked channels information
7390f655 221 //
7d855b04 222 AliTPCCalibRaw *fRefCalibRaw; ///< Reference raw calibration object
949d8707 223 //
7d855b04 224 AliTPCdataQA *fRefDataQA; ///< Reference data QA
6e7d7dc4 225 //
7d855b04 226 AliDCSSensorArray* fGoofieArray; ///< Goofie Data
892226be 227 //
7d855b04 228 AliTPCmapper *fMapper; ///< TPC mapping handler
229 Int_t fNpulserOutliers; ///< number of outliers from Pulser calibration
892226be 230
7d855b04 231 Float_t fIrocTimeOffset; ///< timing offset between IROC and OROC in timebins
232 Float_t fCETmaxLimitAbs; ///< maximum variation in CE data before pads will be treated as outliers
233 Float_t fPulTmaxLimitAbs; ///< maximum variation of Pulser Signals (time) before pads will be treated as outliers
234 Float_t fPulQmaxLimitAbs; ///< maximum variation of Pulser Signals (charge) before pads will be treated as outliers
235 Float_t fPulQminLimit; ///< minimum charge value for Pulser Signals before pads will be treated as outliers
a23ba1c3 236
237 //
238 // helpers to get the run number for given time stamps
239 //
240 // filters
892226be 241
7d855b04 242 TArrayI fRuns; ///< run list with OCDB info
243 TArrayI fRunsStart; ///< start time for given run
244 TArrayI fRunsStop; ///< stop time for given run
8166d5d7 245
892226be 246 AliTPCcalibDButil (const AliTPCcalibDButil& );
247 AliTPCcalibDButil& operator= (const AliTPCcalibDButil& );
248
949d8707 249
78f17711 250 static Double_t GetVDriftTPCLaserTracksCommon(Double_t &dist, Int_t timeStamp, Double_t deltaT, Int_t side, TObjArray * const array);
251
252
7d855b04 253 /// \cond CLASSIMP
892226be 254 ClassDef(AliTPCcalibDButil,0)
7d855b04 255 /// \endcond
892226be 256};
257
258
259#endif