1 #ifndef ALITPCCALIBDBUTIL_H
2 #define ALITPCCALIBDBUTIL_H
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
7 ///////////////////////////////////////////////////////////////////////////////
9 // Class providing the calculation of derived quantities (mean,rms,fits,...) //
10 // of calibration entries //
12 ///////////////////////////////////////////////////////////////////////////////
16 #include <TVectorDfwd.h>
17 #include <TVectorFfwd.h>
18 #include <TMatrixDfwd.h>
22 class AliDCSSensorArray;
30 class AliDCSSensorArray;
31 class AliTPCSensorTempArray;
34 class TTreeSRedirector;
38 class AliTPCcalibDButil : public TObject
42 virtual ~AliTPCcalibDButil();
44 void UpdateFromCalibDB();
45 //data processing functions
46 void ProcessCEdata(const char* fitFormula, TVectorD &fitResultsA, TVectorD &fitResultsC,
47 Int_t &noutliersCE, Double_t & chi2A, Double_t &chi2C, AliTPCCalPad * const outCE=0);
48 void ProcessCEgraphs(TVectorD &vecTEntries, TVectorD &vecTMean, TVectorD &vecTRMS, TVectorD &vecTMedian,
49 TVectorD &vecQEntries, TVectorD &vecQMean, TVectorD &vecQRMS, TVectorD &vecQMedian,
50 Float_t &driftTimeA, Float_t &driftTimeC );
51 void ProcessNoiseData(TVectorD &vNoiseMean, TVectorD &vNoiseMeanSenRegions,
52 TVectorD &vNoiseRMS, TVectorD &vNoiseRMSSenRegions,
53 Int_t &nonMaskedZero, Int_t &nNaN);
54 void ProcessPulser(TVectorD &vMeanTime);
55 void ProcessALTROConfig(Int_t &nMasked);
56 void ProcessGoofie(TVectorD & vecEntries, TVectorD & vecMedian, TVectorD &vecMean, TVectorD &vecRMS);
58 //processing functions using reference data
59 void ProcessPedestalVariations(TVectorF &pedestalDeviations);
60 void ProcessNoiseVariations(TVectorF &noiseDeviations);
61 void ProcessPulserVariations(TVectorF &pulserQdeviations, Float_t &varQMean, Int_t &npadsOutOneTB, Int_t &npadsOffAdd);
63 //getter preprocess information
64 Int_t GetNPulserOutliers() const {return fNpulserOutliers;}
65 Float_t GetMeanAltro(const AliTPCCalROC *roc, const Int_t row, const Int_t pad, AliTPCCalROC * const rocOut=0x0);
66 AliTPCCalPad *GetPulserOutlierMap() const {return fPulserOutlier;}
69 TMap *GetReferenceMap() const {return fRefMap;}
70 Int_t GetReferenceRun(const char* type) const;
71 const char* GetRefValidity() const {return fRefValidity.Data();}
73 AliTPCCalPad* GetRefPadNoise() const {return fRefPadNoise;}
74 AliTPCCalPad* GetRefPedestals() const {return fRefPedestals;}
75 AliTPCCalPad* GetRefPedestalMasked() const {return fRefPedestalMasked;}
76 AliTPCCalPad* GetRefPulserTmean() const {return fRefPulserTmean;}
77 AliTPCCalPad* GetRefPulserTrms() const {return fRefPulserTrms;}
78 AliTPCCalPad* GetRefPulserQmean() const {return fRefPulserQmean;}
79 AliTPCCalPad* GetRefPulserOutlier() const {return fRefPulserOutlier;}
80 AliTPCCalPad* GetRefPulserMasked() const {return fRefPulserMasked;}
81 AliTPCCalPad* GetRefCETmean() const {return fRefCETmean;}
82 AliTPCCalPad* GetRefCETrms() const {return fRefCETrms;}
83 AliTPCCalPad* GetRefCEQmean() const {return fRefCEQmean;}
84 AliTPCCalPad* GetRefCEMasked() const {return fRefCEMasked;}
85 AliTPCCalPad* GetRefALTROFPED() const {return fRefALTROFPED;}
86 AliTPCCalPad* GetRefALTROZsThr() const {return fRefALTROZsThr;}
87 AliTPCCalPad* GetRefALTROAcqStart() const {return fRefALTROAcqStart;}
88 AliTPCCalPad* GetRefALTROAcqStop() const {return fRefALTROAcqStop;}
89 AliTPCCalPad* GetRefALTROMasked() const {return fRefALTROMasked;}
92 //setters for pad by pad information
93 void SetPulserData(AliTPCCalPad * const tmean, AliTPCCalPad * const trms=0x0, AliTPCCalPad * const qmean=0x0)
94 {fPulserTmean=tmean; fPulserTrms=trms; fPulserQmean=qmean;}
95 void SetCEData(AliTPCCalPad *const tmean, AliTPCCalPad *const trms=0x0, AliTPCCalPad *const qmean=0x0)
96 {fCETmean=tmean; fCETrms=trms; fCEQmean=qmean;}
97 void SetNoisePedestal(AliTPCCalPad *const noise, AliTPCCalPad *const pedestal=0x0)
98 {fPadNoise=noise; fPedestals=pedestal;}
99 void SetALTROData(AliTPCCalPad *const masked)
100 {fALTROMasked=masked;}
101 void SetGoofieArray(AliDCSSensorArray *const arr) {fGoofieArray=arr;}
103 //setters for pad by pad information
104 void SetRefFile(const char* filename);
105 void SetReferenceRun(Int_t run=-1);
106 void UpdateRefDataFromOCDB();
107 void SetRefPulserData(AliTPCCalPad *const tmean, AliTPCCalPad *const trms=0x0, AliTPCCalPad *const qmean=0x0)
108 {fRefPulserTmean=tmean; fRefPulserTrms=trms; fRefPulserQmean=qmean;}
109 void SetRefCEData(AliTPCCalPad *const tmean, AliTPCCalPad *const trms=0x0, AliTPCCalPad *const qmean=0x0)
110 {fRefCETmean=tmean; fRefCETrms=trms; fRefCEQmean=qmean;}
111 void SetRefNoisePedestal(AliTPCCalPad *const noise, AliTPCCalPad *const pedestal=0x0)
112 {fRefPadNoise=noise; fRefPedestals=pedestal;}
113 void SetRefALTROData(AliTPCCalPad *const masked)
114 {fRefALTROMasked=masked;}
116 //creation of derived pad by pad calibration data
117 AliTPCCalPad *CreatePadTime0(Int_t model, Double_t &gyA, Double_t &gyC, Double_t &chi2A, Double_t &chi2C);
119 // create outlyer maps
121 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);
122 AliTPCCalPad *CreatePulserOutlyerMap(Int_t &noutliersPulser, AliTPCCalPad * const pulserOut=0, Float_t cutTime=3, Float_t cutnRMSQ=5, Float_t cutnRMSrms=5);
124 AliTPCCalPad *CreatePadTime0CE(TVectorD &fitResultsA, TVectorD&fitResultsC, Int_t &nOut, Double_t &chi2A, Double_t &chi2C, const char *dumpfile=0);
127 void UpdatePulserOutlierMap();
128 void UpdateRefPulserOutlierMap();
129 void PulserOutlierMap(AliTPCCalPad *pulOut, const AliTPCCalPad *pulT, const AliTPCCalPad *pulQ);
131 const char* GetGUIRefTreeDefaultName();
133 Bool_t CreateGUIRefTree(const char* filename="");
137 static Double_t GetLaserTime0(Int_t run, Int_t timeStamp, Int_t deltaT, Int_t side);
138 static TGraph* FilterGraphMedian(TGraph * const graph, Float_t sigmaCut, Double_t &medianY);
139 static TGraph* FilterGraphMedianAbs(TGraph * graph, Float_t cut, Double_t &medianY);
140 static TGraphErrors* FilterGraphMedianErr(TGraphErrors * graph, Float_t sigmaCut,Double_t &medianY);
142 static void Sort(TGraph *graph);
143 static void SmoothGraph(TGraph *graph, Double_t delta);
144 static Int_t GetNearest(TGraph *graph, Double_t xref, Double_t &dx, Double_t &y);
145 static Double_t EvalGraphConst(TGraph * const graph, Double_t xref);
149 static Float_t FilterSensor(AliDCSSensor * sensor, Double_t ymin, Double_t ymax, Double_t maxdy, Double_t sigmaCut);
151 // Filter AliRelAlignmentKalman - Alignment/Drift velocity
153 static TMatrixD* MakeStatRelKalman(TObjArray * const array, Float_t minFraction, Int_t minStat, Float_t maxvd);
154 static TObjArray *SmoothRelKalman(TObjArray * const array,const TMatrixD & stat, Bool_t direction, Float_t sigmaCut);
155 static TObjArray *SmoothRelKalman(TObjArray * const arrayP, TObjArray * const arrayM);
156 static void FilterCE(Double_t deltaT=100, Double_t cutAbs=10, Double_t cutSigma=4., TTreeSRedirector * const pcstream=0);
157 static void FilterTracks(Int_t run, Double_t cutSigma=20., TTreeSRedirector * const pcstream=0);
158 static Float_t FilterTemperature(AliTPCSensorTempArray *tempArray, Double_t ymin=15, Double_t ymax=22, Double_t sigmaCut=5);
160 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);
161 static Double_t GetTriggerOffsetTPC(Int_t run, Int_t timeStamp, Double_t deltaT=86400, Double_t deltaTLaser=3600, Int_t valType=0);
162 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);
163 static Double_t GetVDriftTPCLaserTracks(Double_t &dist,Int_t run, Int_t timeStamp, Double_t deltaT=43200, Int_t side=2);
164 static Double_t GetVDriftTPCCE(Double_t &dist, Int_t run, Int_t timeStamp, Double_t deltaT=43200, Int_t side=2);
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);
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);
170 AliTPCcalibDB *fCalibDB; //pointer to calibDB object
171 AliTPCCalPad *fPadNoise; //noise information
172 AliTPCCalPad *fPedestals; //pedestal information
173 AliTPCCalPad *fPulserTmean; //pulser mean time information
174 AliTPCCalPad *fPulserTrms; //pulser rms time information
175 AliTPCCalPad *fPulserQmean; //pulser mean q information
176 AliTPCCalPad *fPulserOutlier; //pulser outlier map
177 AliTPCCalPad *fCETmean; //central electrode mean time information
178 AliTPCCalPad *fCETrms; //central electrode rms time information
179 AliTPCCalPad *fCEQmean; //central electrode mean q information
180 AliTPCCalPad *fALTROMasked; //ALTRO masked channels information
182 AliTPCCalibRaw *fCalibRaw; //raw calibration object
184 AliTPCdataQA *fDataQA; //data qa
186 TMap *fRefMap; // latest map to reference information
187 TMap *fCurrentRefMap; // reference data map of entries currently loaded
188 TString fRefValidity; // validity range of reference data
190 AliTPCCalPad *fRefPadNoise; //Reference noise information
191 AliTPCCalPad *fRefPedestals; //Reference pedestal information
192 AliTPCCalPad *fRefPedestalMasked; //Reference masked channels in pedestal run
193 AliTPCCalPad *fRefPulserTmean; //Reference pulser mean time information
194 AliTPCCalPad *fRefPulserTrms; //Reference pulser rms time information
195 AliTPCCalPad *fRefPulserQmean; //Reference pulser mean q information
196 AliTPCCalPad *fRefPulserOutlier; //Reference pulser outlier map
197 AliTPCCalPad *fRefPulserMasked; //Reference masked channels in pulser run
198 AliTPCCalPad *fRefCETmean; //Reference central electrode mean time information
199 AliTPCCalPad *fRefCETrms; //Reference central electrode rms time information
200 AliTPCCalPad *fRefCEQmean; //Reference central electrode mean q information
201 AliTPCCalPad *fRefCEMasked; //Reference masked channels in laser run
202 AliTPCCalPad *fRefALTROFPED; //Reference fixed pedestal value
203 AliTPCCalPad *fRefALTROZsThr; //Reference zero suppression threshol
204 AliTPCCalPad *fRefALTROAcqStart; //Reference accquistion start time bin
205 AliTPCCalPad *fRefALTROAcqStop; //Reference accquistion stop time bin
206 AliTPCCalPad *fRefALTROMasked; //Reference ALTRO masked channels information
208 AliTPCCalibRaw *fRefCalibRaw; //Reference raw calibration object
210 AliTPCdataQA *fRefDataQA; //Reference data QA
212 AliDCSSensorArray* fGoofieArray; //Goofie Data
214 AliTPCmapper *fMapper; //TPC mapping handler
215 Int_t fNpulserOutliers; //number of outliers from Pulser calibration
217 Float_t fIrocTimeOffset; //timing offset between IROC and OROC in timebins
218 Float_t fCETmaxLimitAbs; //maximum variation in CE data before pads will be treated as outliers
219 Float_t fPulTmaxLimitAbs; //maximum variation of Pulser Signals (time) before pads will be treated as outliers
220 Float_t fPulQmaxLimitAbs; //maximum variation of Pulser Signals (charge) before pads will be treated as outliers
221 Float_t fPulQminLimit; //minimum charge value for Pulser Signals before pads will be treated as outliers
224 // helpers to get the run number for given time stamps
228 TArrayI fRuns; // run list with OCDB info
229 TArrayI fRunsStart; // start time for given run
230 TArrayI fRunsStop; // stop time for given run
232 AliTPCcalibDButil (const AliTPCcalibDButil& );
233 AliTPCcalibDButil& operator= (const AliTPCcalibDButil& );
235 AliTPCCalPad* GetRefCalPad(AliCDBEntry *entry, const char* objName);
236 AliTPCCalPad* GetRefCalPad(AliCDBEntry *entry);
237 AliTPCCalPad* GetAltroMasked(const char* cdbPath, const char* name);
238 Bool_t HasRefChanged(const char *cdbPath);
239 Int_t GetCurrentReferenceRun(const char* type) const;
240 AliCDBEntry* GetRefEntry(const char* cdbPath);
242 ClassDef(AliTPCcalibDButil,0)