Corrected UInt_t <-> Int_t conversion
[u/mrichter/AliRoot.git] / 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 AliRawReader;
20 class AliTPCRawStream;
21 class AliTPCRawStreamV3;
22 class TNtuple;
23 class TH1F;
24 class TH2F;
25
26 class AliTPCCalibTCF : public TNamed {
27
28 public:
29
30   AliTPCCalibTCF();
31   AliTPCCalibTCF(Int_t gateWidth, Int_t Sample, Int_t pulseLength, Int_t lowPulseLim, Int_t upPulseLim, Double_t rmsLim, Double_t ratioIntLim);
32   AliTPCCalibTCF(const AliTPCCalibTCF &sig);
33   virtual ~AliTPCCalibTCF();
34   
35   AliTPCCalibTCF& operator = (const  AliTPCCalibTCF &source);
36
37  
38   void ProcessRawFileV3(const char *nameRawFile, const char *nameFileOut);
39   void ProcessRawEventV3(AliRawReader *rawReader,AliTPCRawStreamV3 *rawStream, const char *nameFileOut);
40
41   void ProcessRawFile(const char *nameRawFile, const char *nameFileOut, bool bUseHLTOUT=false);
42   void ProcessRawEvent(AliTPCRawStream *rawStream, const char *nameFileOut);
43
44   void MergeHistoPerSector(const char *nameFileIn);
45
46   void AnalyzeRootFile(const char *nameFileIn, Int_t minNumPulse=1, Int_t histStart=1, Int_t histEnd=1000000);
47   Int_t AnalyzePulse(TH1F * const hisIn, Double_t *coefZ, Double_t *coefP); 
48
49   void TestTCFonRootFile(const char *nameFileIn, const char *nameFileTCF, Int_t nPulseMin=0, Int_t plotFlag=0, Int_t lowKey=1, Int_t upKey=1000000);
50   void TestTCFonRawFile(const char *nameRawFile, const char *nameFileOut, const char *nameFileTCF, Int_t nPulseMin=0, Int_t plotFlag=0, bool bUseHLTOUT=false);
51
52   Int_t DumpTCFparamToFilePerSector(const char *nameFileTCFPerSec, const char *nameMappingFile="$ALICE_ROOT/TPC/Calib/tpcMapping.root");  
53   Int_t DumpTCFparamToFilePerPad(const char *nameFileTCFPerPad,const char *nameFileTCFPerSec, const char *nameMappingFile="$ALICE_ROOT/TPC/Calib/tpcMapping.root");  
54  
55   TH2F *PlotOccupSummary2Dhist(const char *nameFileIn, Int_t side=0);
56   void PlotOccupSummary(const char *nameFile, Int_t side=0, Int_t nPulseMin=0); 
57
58   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");
59
60   void PrintPulseThresholds();
61
62   void MergeHistoPerFile(const char *fileNameIn, const char *fileSum, Int_t mode=0);
63   void MergeToOneFile(const char *nameFileSum);
64
65 private:
66   
67   // tresholds for proper pulse finder (Analyze functions)
68   Int_t fGateWidth;     // expected Gate fluctuation length
69   Int_t fSample;        // expected usefull signal length
70   Int_t fPulseLength;   // needed pulselength for TC characterisation
71   Int_t fLowPulseLim;   // lower pulse height limit
72   Int_t fUpPulseLim;    // upper pulse height limit
73   Double_t fRMSLim;     // signal RMS limit
74   Double_t fRatioIntLim;// ratio of signal-integral/pulse-integral limit
75
76   Int_t FitPulse(TNtuple *dataTuple, Double_t *coefZ, Double_t *coefP);
77   static void FitFcn(Int_t &nPar, Double_t *grad, Double_t &f, Double_t * const par, Int_t iflag);
78
79   Double_t* ExtractPZValues(Double_t *param);
80   Int_t Equalization(TNtuple *dataTuple, Double_t *coefZ, Double_t *coefP);
81
82   Int_t FindCorTCFparam(TH1F * const hisIn, const char *nameFileTCF, Double_t *coefZ, Double_t *coefP);
83   Double_t *GetQualityOfTCF(TH1F *hisIn, Double_t *coefZ, Double_t *coefP,Int_t plotFlag=0); 
84
85   TNtuple *ApplyTCFilter(TH1F * const hisIn, Double_t * const coefZ, Double_t * const coefP, Int_t plotFlag=0);
86
87   ClassDef(AliTPCCalibTCF,1);
88
89 };
90 #endif