TPC/AliTPCPreprocessorOffline.h

   1 #ifndef ALITPCPREPROCESSOROFFLINE_H
   2 #define ALITPCPREPROCESSOROFFLINE_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 class TObjArray;
  15 class AliTPCcalibTime;
  16 class AliTPCcalibTimeGain;
  17 class AliTPCcalibGainMult;
  18 class AliTPCROC;
  19 class AliTPCParam;
  20 class TPad;
  21 class AliCDBRunRange;
  22
  23 class AliTPCPreprocessorOffline:public TNamed {
  24 public:
  25   AliTPCPreprocessorOffline();
  26   virtual ~AliTPCPreprocessorOffline();
  27   void UpdateOCDBDrift(Int_t ustartRun, Int_t uendRun, const char* storagePath);
  28   void UpdateOCDBGain(Int_t  startRunNumber, Int_t endRunNumber, const char* storagePath);
  29   void UpdateDriftParam(AliTPCParam *param, TObjArray *const arr, Int_t startRun);
  30
  31   //
  32   // v drift part
  33   //
  34   void GetRunRange(AliTPCcalibTime* const timeDrift);
  35   void CalibTimeVdrift(const Char_t* file, Int_t ustartRun, Int_t uendRun,TString ocdbStorage="");
  36   void AddHistoGraphs(  TObjArray * vdriftArray, AliTPCcalibTime * const timeDrift, Int_t minEntries);
  37   void AddAlignmentGraphs(  TObjArray * vdriftArray, AliTPCcalibTime * const timeDrift);
  38   void AddLaserGraphs(  TObjArray * vdriftArray, AliTPCcalibTime *timeDrift);
  39   void SetDefaultGraphDrift(TGraph *graph, Int_t color, Int_t style);
  40   void MakeDefaultPlots(TObjArray * const arr, TObjArray *picArray);
  41   Bool_t ValidateTimeDrift(Double_t maxVDriftCorr=0.03);
  42   //
  43   // Gain part
  44   //
  45   void CalibTimeGain(const Char_t* fileName, Int_t startRunNumber, Int_t endRunNumber,  TString  ocdbStorage);
  46   void ReadGainGlobal(const Char_t* fileName="CalibObjectsTrain1.root");
  47   void MakeQAPlot(Float_t  FPtoMIPratio);
  48   Bool_t AnalyzeGain(Int_t startRunNumber, Int_t endRunNumber, Int_t minEntriesGaussFit = 500, Float_t FPtoMIPratio = 1.43);
  49   Bool_t AnalyzeAttachment(Int_t startRunNumber, Int_t endRunNumber, Int_t minEntriesFit = 2000);
  50   Bool_t AnalyzePadRegionGain();
  51   Bool_t AnalyzeGainMultiplicity();
  52   Bool_t ValidateTimeGain(Double_t minGain=2.0, Double_t maxGain = 3.0);
  53   //
  54   // Alignment time part
  55   //
  56   void  MakeChainTime();
  57   void  MakePrimitivesTime();
  58   void  CreateAlignTime(TString fstring, TVectorD paramC);
  59   void  MakeFitTime();
  60   static Double_t EvalAt(Double_t phi, Double_t refX, Double_t theta, Int_t corr, Int_t ptype);
  61
  62   //
  63   // QA drawing part
  64   //
  65   static void SetPadStyle(TPad *pad, Float_t mx0, Float_t mx1, Float_t my0, Float_t my1);
  66   static void PrintArray(TObjArray *array);
  67   TChain *GetAlignTree(){return fAlignTree;}
  68   //
  69   // graph filtering part
  70   //
  71   static TGraphErrors* FilterGraphMedianAbs(TGraphErrors * graph, Float_t cut,Double_t &medianY);
  72   static TGraphErrors* FilterGraphDrift(TGraphErrors * graph, Float_t errSigmaCut, Float_t medianCutAbs);
  73   static TGraphErrors * MakeGraphFilter0(THnSparse *hisN, Int_t itime, Int_t ival, Int_t minEntries, Double_t offset=0);
  74
  75   //
  76   void SwitchOnValidation(Bool_t val = kTRUE) {fSwitchOnValidation = val;}
  77   Bool_t IsSwitchOnValidation() { return fSwitchOnValidation; }
  78
  79
  80 private:
  81   Int_t fMinEntries;                      // minimal number of entries for fit
  82   Int_t startRun;                         // start Run - used to make fast selection in THnSparse
  83   Int_t endRun;                           // end   Run - used to make fast selection in THnSparse
  84   Int_t startTime;                        // startTime - used to make fast selection in THnSparse
  85   Int_t endTime;                          // endTime   - used to make fast selection in THnSparse
  86   TString  ocdbStorage;                   // path to the OCDB storage
  87   TObjArray * fVdriftArray;               // array with output calibration graphs
  88   AliTPCcalibTime * fTimeDrift;           // input data to construct calibration graphs
  89   TGraphErrors * fGraphMIP;                // graph time dependence of MIP
  90   TGraphErrors * fGraphCosmic;             // graph time dependence at Plateu
  91   TGraphErrors * fGraphAttachmentMIP;      // graph time dependence of attachment (signal vs. mean driftlength)
  92   AliSplineFit * fFitMIP;                  // fit of dependence - MIP
  93   AliSplineFit * fFitCosmic;               // fit of dependence - Plateu
  94   TObjArray    * fGainArray;               // array to be stored in the OCDB
  95   AliTPCcalibTimeGain * fGainMIP;          // calibration component for MIP
  96   AliTPCcalibTimeGain * fGainCosmic;       // calibration component for cosmic
  97   AliTPCcalibGainMult * fGainMult;         // calibration component for pad region gain equalization and multiplicity dependence
  98
  99   TChain   *fAlignTree;        //alignment tree
 100   //
 101   Bool_t fSwitchOnValidation;  // flag to switch on validation of OCDB parameters
 102
 103 private:
 104   AliTPCPreprocessorOffline& operator=(const AliTPCPreprocessorOffline&); // not implemented
 105   AliTPCPreprocessorOffline(const AliTPCPreprocessorOffline&); // not implemented
 106   ClassDef(AliTPCPreprocessorOffline,1)
 107 };
 108
 109 #endif