1 #ifndef ALITPCPREPROCESSOROFFLINE_H
2 #define ALITPCPREPROCESSOROFFLINE_H
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
9 // Class to create OCDB entries - processing the results of the OFFLINE calibration
15 class AliTPCcalibTime;
16 class AliTPCcalibTimeGain;
17 class AliTPCcalibGainMult;
23 class AliTPCPreprocessorOffline:public TNamed {
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);
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 void SetMaxVDriftCorr(Double_t maxVDriftCorr=0.03) {fMaxVdriftCorr=maxVDriftCorr;};
42 Bool_t ValidateTimeDrift();
46 void CalibTimeGain(const Char_t* fileName, Int_t startRunNumber, Int_t endRunNumber, TString ocdbStorage);
47 void ReadGainGlobal(const Char_t* fileName="CalibObjectsTrain1.root");
48 void MakeQAPlot(Float_t FPtoMIPratio);
49 Bool_t AnalyzeGain(Int_t startRunNumber, Int_t endRunNumber, Int_t minEntriesGaussFit = 500, Float_t FPtoMIPratio = 1.43);
50 Bool_t AnalyzeAttachment(Int_t startRunNumber, Int_t endRunNumber, Int_t minEntriesFit = 2000);
51 Bool_t AnalyzePadRegionGain();
52 Bool_t AnalyzeGainMultiplicity();
53 Bool_t AnalyzeGainChamberByChamber();
54 void SetTimeGainRange(Double_t minGain=2.0, Double_t maxGain = 3.0)
55 {fMinGain = minGain; fMaxGain = maxGain;};
56 Bool_t ValidateTimeGain();
58 // Alignment time part
61 void MakePrimitivesTime();
62 // static void RegisterPrimitiveTimes();
63 void CreateAlignTime(TString fstring, TVectorD paramC);
65 static Double_t EvalAt(Double_t phi, Double_t refX, Double_t theta, Int_t corr, Int_t ptype);
66 static Double_t EvalAtPar(Double_t phi, Double_t snp, Double_t refX, Double_t theta, Int_t corr, Int_t ptype, Int_t nstep);
71 static void SetPadStyle(TPad *pad, Float_t mx0, Float_t mx1, Float_t my0, Float_t my1);
72 static void PrintArray(TObjArray *array);
73 TChain *GetAlignTree(){return fAlignTree;}
75 // graph filtering part
77 static TGraphErrors* FilterGraphMedianAbs(TGraphErrors * graph, Float_t cut,Double_t &medianY);
78 static TGraphErrors* FilterGraphDrift(TGraphErrors * graph, Float_t errSigmaCut, Float_t medianCutAbs);
79 static TGraphErrors * MakeGraphFilter0(THnSparse *hisN, Int_t itime, Int_t ival, Int_t minEntries, Double_t offset=0);
82 void SwitchOnValidation(Bool_t val = kTRUE) {fSwitchOnValidation = val;}
83 Bool_t IsSwitchOnValidation() { return fSwitchOnValidation; }
87 Int_t fMinEntries; // minimal number of entries for fit
88 Int_t startRun; // start Run - used to make fast selection in THnSparse
89 Int_t endRun; // end Run - used to make fast selection in THnSparse
90 Int_t startTime; // startTime - used to make fast selection in THnSparse
91 Int_t endTime; // endTime - used to make fast selection in THnSparse
92 TString ocdbStorage; // path to the OCDB storage
93 TObjArray * fVdriftArray; // array with output calibration graphs
94 AliTPCcalibTime * fTimeDrift; // input data to construct calibration graphs
95 TGraphErrors * fGraphMIP; // graph time dependence of MIP
96 TGraphErrors * fGraphCosmic; // graph time dependence at Plateu
97 TGraphErrors * fGraphAttachmentMIP; // graph time dependence of attachment (signal vs. mean driftlength)
98 AliSplineFit * fFitMIP; // fit of dependence - MIP
99 AliSplineFit * fFitCosmic; // fit of dependence - Plateu
100 TObjArray * fGainArray; // array to be stored in the OCDB
101 AliTPCcalibTimeGain * fGainMIP; // calibration component for MIP
102 AliTPCcalibTimeGain * fGainCosmic; // calibration component for cosmic
103 AliTPCcalibGainMult * fGainMult; // calibration component for pad region gain equalization and multiplicity dependence
105 TChain *fAlignTree; //alignment tree
107 Bool_t fSwitchOnValidation; // flag to switch on validation of OCDB parameters
108 Float_t fMinGain; // min gain
109 Float_t fMaxGain; // max gain
110 Float_t fMaxVdriftCorr; // max v-drift correction
113 AliTPCPreprocessorOffline& operator=(const AliTPCPreprocessorOffline&); // not implemented
114 AliTPCPreprocessorOffline(const AliTPCPreprocessorOffline&); // not implemented
115 ClassDef(AliTPCPreprocessorOffline,1)