| b3fcfd96 |
1 | #ifndef ALITRDPREPROCESSOROFFLINE_H |
| 2 | #define ALITRDPREPROCESSOROFFLINE_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 to create OCDB entries - processing the results of the OFFLINE calibration |
| 10 | // |
| 11 | |
| 12 | |
| 13 | #include "TNamed.h" |
| 81a5aeca |
14 | #include "AliTRDCalChamberStatus.h" |
| 83d0cc79 |
15 | #include <THnSparse.h> |
| b3fcfd96 |
16 | class TObjArray; |
| 01239968 |
17 | class AliTRDCalDet; |
| b3fcfd96 |
18 | class TH2I; |
| 19 | class TProfile2D; |
| 20 | class AliTRDCalibraVdriftLinearFit; |
| a0bb5615 |
21 | class AliTRDCalibraExbAltFit; |
| b3fcfd96 |
22 | class TH1I; |
| 23 | class TH2F; |
| 01239968 |
24 | class TString; |
| b3fcfd96 |
25 | |
| b3fcfd96 |
26 | class AliTRDPreprocessorOffline:public TNamed { |
| 27 | public: |
| a0bb5615 |
28 | enum{ kGain = 0, |
| 29 | kVdriftPHDet = 1, |
| 30 | kVdriftPHPad = 2, |
| 31 | kT0PHDet = 3, |
| 32 | kT0PHPad = 4, |
| 33 | kVdriftLinear = 5, |
| 34 | kLorentzLinear = 6, |
| 35 | kChamberStatus = 7, |
| 36 | kPRF = 8, |
| 37 | kExbAlt = 9, |
| e2a1c98b |
38 | kPHQ = 10, |
| 39 | kNumCalibObjs = 11 |
| a0bb5615 |
40 | }; |
| 840ec79d |
41 | enum { kGainNotEnoughStatsButFill = 2, |
| 42 | kVdriftNotEnoughStatsButFill = 4, |
| 43 | kGainNotEnoughStatsNotFill = 8, |
| 44 | kVdriftNotEnoughStatsNotFill = 16, |
| 45 | kTimeOffsetNotEnoughStatsNotFill = 32, |
| 46 | kExBErrorRange = 64}; |
| 47 | enum { kGainErrorOld = 2, |
| 48 | kVdriftErrorOld = 4, |
| 49 | kExBErrorOld = 8, |
| 50 | kGainErrorRange = 16, |
| 51 | kVdriftErrorRange = 32, |
| 52 | kTimeOffsetErrorRange = 64, |
| 53 | kChamberStatusErrorRange = 128}; |
| 82b413fd |
54 | |
| 55 | |
| b3fcfd96 |
56 | AliTRDPreprocessorOffline(); |
| 57 | virtual ~AliTRDPreprocessorOffline(); |
| 58 | |
| 83d0cc79 |
59 | Bool_t Init(const Char_t* fileName); |
| 60 | void Process(const Char_t* file, Int_t startRunNumber, Int_t endRunNumber, TString ocdbStorage); |
| 61 | |
| 62 | // settings |
| 63 | void SetLinearFitForVdrift(Bool_t methodsecond) { fMethodSecond = methodsecond;}; |
| 64 | void SetNameList(TString nameList) { fNameList = nameList;}; |
| 65 | void SetMinStatsVdriftT0PH(Int_t minStatsVdriftT0PH) { fMinStatsVdriftT0PH = minStatsVdriftT0PH; } |
| 66 | void SetMinStatsVdriftLinear(Int_t minStatsVdriftLinear) { fMinStatsVdriftLinear = minStatsVdriftLinear; } |
| 67 | void SetMinStatsGain(Int_t minStatsGain) { fMinStatsGain = minStatsGain; } |
| 68 | void SetMinStatsPRF(Int_t minStatsPRF) { fMinStatsPRF = minStatsPRF; } |
| 69 | void SetMinStatsChamberStatus(Int_t minStatsChamberStatus) { fMinStatsChamberStatus = minStatsChamberStatus; } |
| 532be2d4 |
70 | void SetSingleMinStatsChamberStatus(Int_t minSingleStatsChamberStatus) { fMinSingleStatsChamberStatus = minSingleStatsChamberStatus; } |
| 83d0cc79 |
71 | void SetLimitValidateNoData(Int_t nodatavalidate) { fNoDataValidate = nodatavalidate; }; |
| 72 | void SetLimitValidateBadCalib(Int_t badcalibvalidate) { fBadCalibValidate = badcalibvalidate; }; |
| 73 | void SetBackCorrectGain(Bool_t backCorrectGain) { fBackCorrectGain = backCorrectGain; } |
| 74 | void SetBackCorrectVdrift(Bool_t backCorrectVdrift) { fBackCorrectVdrift = backCorrectVdrift; } |
| 75 | void SetNoExBUsedInReco(Bool_t noExBUsedInReco) { fNoExBUsedInReco = noExBUsedInReco; }; |
| 76 | void SetSwitchOnValidation(Bool_t switchOnValidation) { fSwitchOnValidation = switchOnValidation;}; |
| 77 | void SetRMSBadCalibratedGain(Double_t rms) { fRMSBadCalibratedGain = rms;}; |
| 78 | void SetRMSBadCalibratedVdrift(Double_t rms) { fRMSBadCalibratedVdrift = rms;}; |
| 79 | void SetRMSBadCalibratedExB(Double_t rms) { fRMSBadCalibratedExB = rms;}; |
| 2a1a7b36 |
80 | void SetRobustFitDriftVelocity(Bool_t robustFitDriftVelocity) { fRobustFitDriftVelocity = robustFitDriftVelocity;}; |
| 81 | void SetRobustFitExbAlt(Bool_t robustFitExbAlt) { fRobustFitExbAlt = robustFitExbAlt;}; |
| 661c7be6 |
82 | void SetAlternativeDriftVelocityFit(Bool_t alt) { fAlternativeVdrfitFit = alt;}; |
| 532be2d4 |
83 | void SetMinNbOfPointVdriftFit(Int_t minNbOfPointVdriftFit) { fMinNbOfPointVdriftFit = minNbOfPointVdriftFit;}; |
| 661c7be6 |
84 | void SetAlternativeExbAltFit(Bool_t alt) { fAlternativeExbAltFit = alt;}; |
| 595cfc12 |
85 | void SetMethodeFitGain(Int_t methodeFitGain) { fMethodeGain = methodeFitGain;}; |
| 86 | void SetOutliersFitChargeLow(Float_t outliersFitChargeLow) { fOutliersFitChargeLow = outliersFitChargeLow; } |
| 87 | void SetOutliersFitChargeHigh(Float_t outliersFitChargeHigh) { fOutliersFitChargeHigh = outliersFitChargeHigh; } |
| 88 | void SetBeginFitCharge(Float_t beginFitCharge) { fBeginFitCharge = beginFitCharge;}; |
| 89 | |
| 2a1a7b36 |
90 | |
| 83d0cc79 |
91 | |
| 92 | Bool_t GetLinearFitForVdrift() const { return fMethodSecond;}; |
| 93 | TString GetNameList() const { return fNameList;}; |
| 94 | |
| 95 | // status |
| 96 | Bool_t CheckStatus(Int_t status, Int_t bitMask) const; |
| 97 | void PrintStatus() const; |
| 98 | |
| 82b413fd |
99 | Bool_t IsGainNotEnoughStatsButFill() const |
| 100 | { return CheckStatus(fStatusNeg, kGainNotEnoughStatsButFill); }; |
| 101 | Bool_t IsGainNotEnoughStatsNotFill() const |
| 102 | { return CheckStatus(fStatusNeg, kGainNotEnoughStatsNotFill); }; |
| 103 | Bool_t IsVdriftNotEnoughStatsButFill() const |
| 104 | { return CheckStatus(fStatusNeg, kVdriftNotEnoughStatsButFill); }; |
| 105 | Bool_t IsVdriftNotEnoughStatsNotFill() const |
| 106 | { return CheckStatus(fStatusNeg, kVdriftNotEnoughStatsNotFill); }; |
| 107 | Bool_t IsTimeOffsetNotEnoughStatsNotFill() const |
| 840ec79d |
108 | { return CheckStatus(fStatusNeg, kTimeOffsetNotEnoughStatsNotFill); }; |
| 109 | Bool_t IsExBErrorRange() const |
| 110 | { return CheckStatus(fStatusNeg, kExBErrorRange); }; |
| 82b413fd |
111 | |
| 112 | Bool_t IsGainErrorOld() const |
| 113 | { return CheckStatus(fStatusPos, kGainErrorOld); }; |
| 114 | Bool_t IsVdriftErrorOld() const |
| 115 | { return CheckStatus(fStatusPos, kVdriftErrorOld); }; |
| 840ec79d |
116 | Bool_t IsExBErrorOld() const |
| 117 | { return CheckStatus(fStatusPos, kExBErrorOld); }; |
| 82b413fd |
118 | Bool_t IsGainErrorRange() const |
| 119 | { return CheckStatus(fStatusPos, kGainErrorRange); }; |
| 120 | Bool_t IsVdriftErrorRange() const |
| 121 | { return CheckStatus(fStatusPos, kVdriftErrorRange); }; |
| 122 | Bool_t IsTimeOffsetErrorRange() const |
| 123 | { return CheckStatus(fStatusPos, kTimeOffsetErrorRange); }; |
| 124 | Bool_t IsChamberStatusErrorRange() const |
| 125 | { return CheckStatus(fStatusPos, kChamberStatusErrorRange); }; |
| 126 | |
| 82b413fd |
127 | |
| 83d0cc79 |
128 | // Back corrections |
| 01239968 |
129 | void SetCalDetGain(AliTRDCalDet *calDetGainUsed) {fCalDetGainUsed = calDetGainUsed;}; |
| 840ec79d |
130 | void SetCalDetVdrift(AliTRDCalDet *calDetVdriftUsed); |
| 131 | void SetCalDetVdriftExB(AliTRDCalDet *calDetVdriftUsed,AliTRDCalDet *calDetExBUsed) {fCalDetVdriftUsed = calDetVdriftUsed; fCalDetExBUsed = calDetExBUsed;}; |
| 83d0cc79 |
132 | Bool_t SetCalDetGain(Int_t runNumber, Int_t version, Int_t subversion); |
| 133 | Bool_t SetCalDetVdriftExB(Int_t runNumber, Int_t versionv, Int_t subversionv, Int_t versionexb, Int_t subversionexb); |
| 134 | |
| 01239968 |
135 | AliTRDCalDet *GetCalDetGain() const { return fCalDetGainUsed;}; |
| 4c865c34 |
136 | AliTRDCalDet *GetCalDetVdrift() const { return fCalDetVdriftUsed;}; |
| 83d0cc79 |
137 | Int_t GetFirstRunGainUsed() const { return fFirstRunGainUsed; } |
| 138 | Int_t GetVersionGainUsed() const { return fVersionGainUsed; } |
| 139 | Int_t GetSubVersionGainUsed() const { return fSubVersionGainUsed; } |
| 140 | Int_t GetFirstRunVdriftUsed() const { return fFirstRunVdriftUsed; } |
| 141 | Int_t GetVersionVdriftUsed() const { return fVersionVdriftUsed; } |
| 142 | Int_t GetSubVersionVdriftUsed() const { return fSubVersionVdriftUsed; } |
| 143 | Int_t GetFirstRunExBUsed() const { return fFirstRunExBUsed; } |
| 144 | Int_t GetVersionExBUsed() const { return fVersionExBUsed; } |
| 145 | Int_t GetSubVersionExBUsed() const { return fSubVersionExBUsed; } |
| 146 | |
| 147 | |
| 148 | // Internal functions |
| b3fcfd96 |
149 | |
| 150 | void CalibVdriftT0(const Char_t* file, Int_t startRunNumber, Int_t endRunNumber, TString ocdbStorage=""); |
| a0bb5615 |
151 | void CalibExbAlt(const Char_t* file, Int_t startRunNumber, Int_t endRunNumber, TString ocdbStorage=""); |
| b3fcfd96 |
152 | void CalibGain(const Char_t* file, Int_t startRunNumber, Int_t endRunNumber, TString ocdbStorage=""); |
| 153 | void CalibPRF(const Char_t* file, Int_t startRunNumber, Int_t endRunNumber, TString ocdbStorage=""); |
| 83d0cc79 |
154 | void CalibChamberStatus(const Char_t* file, Int_t startRunNumber, Int_t endRunNumber, TString ocdbStorage=""); |
| e2a1c98b |
155 | void CalibPHQ(const Char_t* file, Int_t startRunNumber, Int_t endRunNumber, TString ocdbStorage); |
| 81a5aeca |
156 | |
| 83d0cc79 |
157 | Bool_t ReadStatusGlobal(const Char_t* fileName="CalibObjects.root"); |
| b3fcfd96 |
158 | Bool_t ReadGainGlobal(const Char_t* fileName="CalibObjects.root"); |
| 159 | Bool_t ReadVdriftT0Global(const Char_t* fileName="CalibObjects.root"); |
| 160 | Bool_t ReadVdriftLinearFitGlobal(const Char_t* fileName="CalibObjects.root"); |
| a0bb5615 |
161 | Bool_t ReadExbAltFitGlobal(const Char_t* fileName="CalibObjects.root"); |
| b3fcfd96 |
162 | Bool_t ReadPRFGlobal(const Char_t* fileName="CalibObjects.root"); |
| e2a1c98b |
163 | Bool_t ReadPHQGlobal(const Char_t* fileName); |
| b3fcfd96 |
164 | |
| 165 | Bool_t AnalyzeGain(); |
| 166 | Bool_t AnalyzeVdriftT0(); |
| 167 | Bool_t AnalyzeVdriftLinearFit(); |
| a0bb5615 |
168 | Bool_t AnalyzeExbAltFit(); |
| a5dcf618 |
169 | Bool_t AnalyzePRF(); |
| 170 | Bool_t AnalyzeChamberStatus(); |
| e2a1c98b |
171 | Bool_t AnalyzePHQ(Int_t startRunNumber); |
| 172 | |
| 01239968 |
173 | void CorrectFromDetGainUsed(); |
| 4c865c34 |
174 | void CorrectFromDetVdriftUsed(); |
| a5dcf618 |
175 | |
| b3fcfd96 |
176 | void UpdateOCDBT0(Int_t startRunNumber, Int_t endRunNumber, const char* storagePath); |
| 177 | void UpdateOCDBVdrift(Int_t startRunNumber, Int_t endRunNumber, const char* storagePath); |
| eec29d56 |
178 | void UpdateOCDBExB(Int_t startRunNumber, Int_t endRunNumber, const Char_t *storagePath); |
| a0bb5615 |
179 | void UpdateOCDBExBAlt(Int_t startRunNumber, Int_t endRunNumber, const Char_t *storagePath); |
| b3fcfd96 |
180 | void UpdateOCDBGain(Int_t startRunNumber, Int_t endRunNumber, const char* storagePath); |
| 181 | void UpdateOCDBPRF(Int_t startRunNumber, Int_t endRunNumber, const char* storagePath); |
| a5dcf618 |
182 | void UpdateOCDBChamberStatus(Int_t startRunNumber, Int_t endRunNumber, const Char_t *storagePath); |
| e2a1c98b |
183 | void UpdateOCDBPHQ(Int_t startRunNumber, Int_t endRunNumber, const Char_t *storagePath); |
| b3fcfd96 |
184 | |
| 54f2ff1c |
185 | Bool_t ValidateGain(); |
| 00d203b6 |
186 | Bool_t ValidateVdrift(); |
| 840ec79d |
187 | Bool_t ValidateExB(); |
| 00d203b6 |
188 | Bool_t ValidateT0(); |
| 189 | Bool_t ValidatePRF() const; |
| 82b413fd |
190 | Bool_t ValidateChamberStatus(); |
| 00d203b6 |
191 | |
| 82b413fd |
192 | Int_t GetStatus() const; |
| 193 | Int_t GetStatusPos() const { return fStatusPos; } |
| 194 | Int_t GetStatusNeg() const { return fStatusNeg; } |
| a5dcf618 |
195 | |
| 9a8b0e85 |
196 | Bool_t IsPHQon() const { return fPHQon ;}; |
| 197 | void SetPHQon(const Bool_t kphq){ fPHQon = kphq; } |
| 83d0cc79 |
198 | |
| 687aa844 |
199 | Bool_t IsDebugPHQon() const { return fDebugPHQon ;}; |
| 200 | void SetDebugPHQon(const Bool_t kphq){ fDebugPHQon = kphq; } |
| 201 | |
| 4c865c34 |
202 | private: |
| 81a5aeca |
203 | Bool_t fMethodSecond; // Second Method for drift velocity |
| 204 | TString fNameList; // Name of the list |
| 205 | AliTRDCalDet *fCalDetGainUsed; // CalDet used and to be corrected for |
| 206 | AliTRDCalDet *fCalDetVdriftUsed; // CalDet used and to be corrected for |
| 840ec79d |
207 | AliTRDCalDet *fCalDetExBUsed; // CalDet used and to be corrected for |
| 81a5aeca |
208 | TH2I *fCH2d; // Gain |
| 209 | TProfile2D *fPH2d; // Drift velocity first method |
| 210 | TProfile2D *fPRF2d; // PRF |
| 83d0cc79 |
211 | THnSparseI *fSparse; // chamberstatus |
| b3fcfd96 |
212 | AliTRDCalibraVdriftLinearFit *fAliTRDCalibraVdriftLinearFit; // Drift velocity second method |
| a0bb5615 |
213 | AliTRDCalibraExbAltFit* fAliTRDCalibraExbAltFit; //ExB alternative method |
| b3fcfd96 |
214 | TH1I *fNEvents; // Number of events |
| 215 | TH2F *fAbsoluteGain; // Absolute Gain calibration |
| 216 | TObjArray * fPlots; // array with some plots to check |
| 217 | TObjArray * fCalibObjects; // array with calibration objects |
| 83d0cc79 |
218 | Int_t fFirstRunGainUsed; // first run gain used |
| 4c865c34 |
219 | Int_t fVersionGainUsed; // VersionGainUsed |
| 220 | Int_t fSubVersionGainUsed; // SubVersionGainUsed |
| ca7e6e64 |
221 | Int_t fFirstRunVdriftUsed; // FirstRunVdrift |
| 4c865c34 |
222 | Int_t fVersionVdriftUsed; // VersionVdriftUsed |
| 223 | Int_t fSubVersionVdriftUsed; // SubVersionVdriftUsed |
| 83d0cc79 |
224 | Int_t fFirstRunExBUsed; // FirstRunExB |
| 225 | Int_t fVersionExBUsed; // VersionExBUsed |
| 226 | Int_t fSubVersionExBUsed; // SubVersionExBUsed |
| 227 | Bool_t fNoExBUsedInReco; // ExB not used yet in the reco |
| 00d203b6 |
228 | Bool_t fSwitchOnValidation; // Validation |
| 229 | Bool_t fVdriftValidated; // Vdrift validation |
| 840ec79d |
230 | Bool_t fExBValidated; // ExB validation |
| 00d203b6 |
231 | Bool_t fT0Validated; // T0 validation |
| a2a4ec8e |
232 | Int_t fMinStatsVdriftT0PH; // MinStats VdriftT0 |
| 233 | Int_t fMinStatsVdriftLinear; // MinStats Vdrift Linear |
| 234 | Int_t fMinStatsGain; // MinStats Gain |
| 235 | Int_t fMinStatsPRF; // MinStats PRF |
| 83d0cc79 |
236 | Int_t fMinStatsChamberStatus; // MinStats ChamberStatus |
| 17aa2935 |
237 | Double_t fMinSingleStatsChamberStatus; // MinStats per chamber in % of mean (ChamberStatus) |
| 54f2ff1c |
238 | Bool_t fBackCorrectGain; // Back correction afterwards gain |
| 239 | Bool_t fBackCorrectVdrift; // Back correction afterwards vdrift |
| 240 | Bool_t fNotEnoughStatisticsForTheGain;// Take the chamber per chamber distribution from the default distribution |
| 241 | Bool_t fNotEnoughStatisticsForTheVdriftLinear;// Take the chamber per chamber distribution from the default distribution |
| 82b413fd |
242 | Int_t fStatusNeg; // Info but ok |
| 243 | Int_t fStatusPos; // Problems |
| 83d0cc79 |
244 | Int_t fBadCalib[18]; // number of bad calibrated chambers per sm |
| 245 | Int_t fNoData[18]; // number of chambers w/o data per sm |
| 246 | Int_t fBadCalibValidate; // validation limit for bad calibrated chambers |
| 247 | Int_t fNoDataValidate; // validation limit for chamber w/o data (sm w/o data excluded) |
| 248 | Double_t fRMSBadCalibratedGain; // value to decide when it is bad calibrated |
| 249 | Double_t fRMSBadCalibratedVdrift; // value to decide when it is bad calibrated |
| 250 | Double_t fRMSBadCalibratedExB; // value to decide when it is bad calibrated |
| 2a1a7b36 |
251 | Bool_t fRobustFitDriftVelocity; // Robust fit for the drift velocity |
| 252 | Bool_t fRobustFitExbAlt; // Robust fit for the exb alt |
| 661c7be6 |
253 | Bool_t fAlternativeVdrfitFit; // Alternative fitting method for vdrift calibration |
| 254 | Bool_t fAlternativeExbAltFit; // Alternative fitting method for the alternative exb calibarion method |
| 532be2d4 |
255 | Int_t fMinNbOfPointVdriftFit; // Min number of points for the drift velocity calibration |
| 595cfc12 |
256 | Int_t fMethodeGain; // Methode Gain Fit |
| 257 | Float_t fOutliersFitChargeLow; // The fit starts at fOutliersFitChargeLow procent number of entries |
| 258 | Float_t fOutliersFitChargeHigh; // The fit starts at fOutliersFitChargeHigh procent number of entries |
| 259 | Float_t fBeginFitCharge; // Fit Begin Charge starts at mean/fBeginFitCharge |
| b3fcfd96 |
260 | |
| 4c865c34 |
261 | Int_t GetSubVersion(TString name) const; |
| 262 | Int_t GetVersion(TString name) const; |
| ca7e6e64 |
263 | Int_t GetFirstRun(TString name) const; |
| b3fcfd96 |
264 | |
| 265 | |
| 266 | |
| 267 | private: |
| 268 | AliTRDPreprocessorOffline& operator=(const AliTRDPreprocessorOffline&); // not implemented |
| 269 | AliTRDPreprocessorOffline(const AliTRDPreprocessorOffline&); // not implemented |
| 9a8b0e85 |
270 | |
| 271 | Bool_t fPHQon; //switch of PHQ |
| 687aa844 |
272 | Bool_t fDebugPHQon; //switch of DebugPHQ |
| 9a8b0e85 |
273 | |
| 532be2d4 |
274 | ClassDef(AliTRDPreprocessorOffline,4) |
| b3fcfd96 |
275 | }; |
| 276 | |
| 277 | #endif |
| 01239968 |
278 | |
| a5dcf618 |
279 | |
git://git.uio.no / u / mrichter / AliRoot.git / blame