TRD/AliTRDPreprocessorOffline.h

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