TPC/AliTPCCalibTCF.h

   1 #ifndef ALI_TPC_CALIB_TCF_H
   2 #define ALI_TPC_CALIB_TCF_H
   3
   4
   5 /* Copyright(c) 2007-08, ALICE Experiment at CERN, All rights reserved. *
   6  * See cxx source for full Copyright notice                             */
   7
   8 ///////////////////////////////////////////////////////////////////////////////
   9 //                        Class AliTPCCalibTCF                               //
  10 // Class for Extraction and test of TCF parameters needed by the ALTRO chip  //
  11 ///////////////////////////////////////////////////////////////////////////////
  12
  13 #include "TSystem.h"
  14 #include <exception>
  15
  16
  17 class TObject;
  18 class TGraph;
  19 class AliTPCRawStream;
  20 class TNtuple;
  21 class TH1F;
  22 class TH2F;
  23
  24 class AliTPCCalibTCF : public TNamed {
  25
  26 public:
  27
  28   AliTPCCalibTCF();
  29   AliTPCCalibTCF(Int_t gateWidth, Int_t Sample, Int_t pulseLength, Int_t lowPulseLim, Int_t upPulseLim, Double_t rmsLim, Double_t ratioIntLim);
  30   AliTPCCalibTCF(const AliTPCCalibTCF &sig);
  31   virtual ~AliTPCCalibTCF();
  32
  33   AliTPCCalibTCF& operator = (const  AliTPCCalibTCF &source);
  34
  35
  36   void ProcessRawFile(const char *nameRawFile, const char *nameFileOut);//Int_t *tresholds);
  37   void ProcessRawEvent(AliTPCRawStream *rawStream, const char *nameFileOut);
  38
  39   void MergeHistoPerSector(const char *nameFileIn);
  40
  41   void AnalyzeRootFile(const char *nameFileIn, Int_t minNumPulse=1, Int_t histStart=1, Int_t histEnd=1000000);
  42   Int_t AnalyzePulse(TH1F *hisIn, Double_t *coefZ, Double_t *coefP);
  43
  44   void TestTCFonRootFile(const char *nameFileIn, const char *nameFileTCF, Int_t nPulseMin=0, Int_t plotFlag=0, Int_t lowKey=1, Int_t upKey=1000000);
  45   void TestTCFonRawFile(const char *nameRawFile, const char *nameFileOut, const char *nameFileTCF, Int_t nPulseMin=0, Int_t plotFlag=0);
  46
  47   Int_t DumpTCFparamToFilePerSector(const char *nameFileTCFPerSec, const char *nameMappingFile="$ALICE_ROOT/TPC/Calib/tpcMapping.root");
  48   Int_t DumpTCFparamToFilePerPad(const char *nameFileTCFPerPad,const char *nameFileTCFPerSec, const char *nameMappingFile="$ALICE_ROOT/TPC/Calib/tpcMapping.root");
  49
  50   TH2F *PlotOccupSummary2Dhist(const char *nameFileIn, Int_t side=0);
  51   void PlotOccupSummary(const char *nameFile, Int_t side=0, Int_t nPulseMin=0);
  52
  53   void PlotQualitySummary(const char *nameFileQuality, const char *plotSpec="widthRed:maxUndershot", const char *cut="maxUndershot<0.1&&maxUndershot>-40&&widthRed>0&&widthRed<100", const char *pOpt="LEGO2Z");
  54
  55   void PrintPulseThresholds();
  56
  57   void MergeHistoPerFile(const char *fileNameIn, const char *fileSum, Int_t mode=0);
  58   void MergeToOneFile(const char *nameFileSum);
  59
  60 private:
  61
  62   // tresholds for proper pulse finder (Analyze functions)
  63   Int_t fGateWidth;     // expected Gate fluctuation length
  64   Int_t fSample;        // expected usefull signal length
  65   Int_t fPulseLength;   // needed pulselength for TC characterisation
  66   Int_t fLowPulseLim;   // lower pulse height limit
  67   Int_t fUpPulseLim;    // upper pulse height limit
  68   Double_t fRMSLim;     // signal RMS limit
  69   Double_t fRatioIntLim;// ratio of signal-integral/pulse-integral limit
  70
  71   Int_t FitPulse(TNtuple *dataTuple, Double_t *coefZ, Double_t *coefP);
  72   static void FitFcn(Int_t &nPar, Double_t *grad, Double_t &f, Double_t *par, Int_t iflag);
  73
  74   Double_t* ExtractPZValues(Double_t *param);
  75   Int_t Equalization(TNtuple *dataTuple, Double_t *coefZ, Double_t *coefP);
  76
  77   Int_t FindCorTCFparam(TH1F *hisIn, const char *nameFileTCF, Double_t *coefZ, Double_t *coefP);
  78   Double_t *GetQualityOfTCF(TH1F *hisIn, Double_t *coefZ, Double_t *coefP,Int_t plotFlag=0);
  79
  80   TNtuple *ApplyTCFilter(TH1F *hisIn, Double_t *coefZ, Double_t *coefP, Int_t plotFlag=0);
  81
  82   ClassDef(AliTPCCalibTCF,1);
  83
  84 };
  85 #endif