]>
Commit | Line | Data |
---|---|---|
75d8233f | 1 | #ifndef ALITPCCALIBCE_H |
2 | #define ALITPCCALIBCE_H | |
3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
4 | * See cxx source for full Copyright notice */ | |
5 | ||
6 | #include <TVectorT.h> | |
7 | class TObjArray; | |
8 | class TH1S; | |
9 | class TH2S; | |
10 | class TH1F; | |
11 | class TTreeSRedirector; | |
12 | class AliTPCCalPad; | |
13 | class AliTPCROC; | |
14 | class AliTPCCalROC; | |
15 | class AliTPCParam; | |
16 | class AliRawReader; | |
17 | class AliTPCRawStream; | |
18 | class TGraph; | |
19 | struct eventHeaderStruct; | |
20 | ||
21 | class AliTPCCalibCE : public TObject { | |
22 | ||
23 | public: | |
24 | AliTPCCalibCE(); | |
25 | AliTPCCalibCE(const AliTPCCalibCE &sig); | |
26 | virtual ~AliTPCCalibCE(); | |
27 | ||
28 | AliTPCCalibCE& operator = (const AliTPCCalibCE &source); | |
29 | ||
30 | ||
31 | Bool_t ProcessEvent(AliTPCRawStream *rawStream); | |
32 | Bool_t ProcessEvent(AliRawReader *rawReader); | |
33 | Bool_t ProcessEvent(eventHeaderStruct *event); | |
34 | ||
35 | Int_t Update(const Int_t isector, const Int_t iRow, const Int_t iPad, | |
36 | const Int_t iTimeBin, const Float_t signal); | |
37 | void Analyse(); | |
38 | // | |
39 | AliTPCCalROC* GetCalRocT0 (Int_t sector, Bool_t force=kFALSE); // get calibration object - sector | |
40 | AliTPCCalROC* GetCalRocQ (Int_t sector, Bool_t force=kFALSE); // get calibration object - sector | |
41 | AliTPCCalROC* GetCalRocRMS(Int_t sector, Bool_t force=kFALSE); // get calibration object - sector | |
42 | AliTPCCalROC* GetCalRocOutliers(Int_t sector, Bool_t force=kFALSE); // get calibration object - sector | |
43 | ||
44 | const TObjArray* GetCalPadT0() { return &fCalRocArrayT0; } // get calibration object | |
45 | const TObjArray* GetCalPadQ() { return &fCalRocArrayQ; } // get calibration object | |
46 | const TObjArray* GetCalPadRMS(){ return &fCalRocArrayRMS;} // get calibration object | |
47 | const TObjArray* GetCalPadOutliers(){ return &fCalRocArrayOutliers;} // get calibration object | |
48 | ||
49 | TH2S* GetHistoQ (Int_t sector, Bool_t force=kFALSE); // get refernce histogram | |
50 | TH2S* GetHistoT0 (Int_t sector, Bool_t force=kFALSE); // get refernce histogram | |
51 | TH2S* GetHistoRMS(Int_t sector, Bool_t force=kFALSE); // get refernce histogram | |
52 | ||
53 | TH1S* GetHistoTmean(Int_t sector, Bool_t force=kFALSE); // get refernce histogram | |
54 | ||
7fb602b1 | 55 | //needed here to merge ClibCE objects |
56 | TObjArray* GetParamArrayPol1(Int_t sector, Bool_t force=kFALSE); | |
57 | TObjArray* GetParamArrayPol2(Int_t sector, Bool_t force=kFALSE); | |
58 | ||
59 | // TObjArray* GetTMeanArrayEvent(){ return &fTMeanArrayEvent; } | |
60 | // TObjArray* GetQMeanArrayEvent(){ return &fQMeanArrayEvent; } | |
61 | TVectorF* GetTMeanEvents(Int_t sector, Bool_t force=kFALSE); | |
62 | TVectorF* GetQMeanEvents(Int_t sector, Bool_t force=kFALSE); | |
63 | ||
64 | TVectorD* GetEventTimes() { return &fVEventTime; } | |
65 | TVectorD* GetEventIds() { return &fVEventNumber; } | |
66 | ||
75d8233f | 67 | Short_t GetDebugLevel() const { return fDebugLevel; } |
68 | // | |
69 | void SetRangeTime (Int_t firstTimeBin, Int_t lastTimeBin) { fFirstTimeBin=firstTimeBin; fLastTimeBin=lastTimeBin; } //Set range in which the pulser signal is expected | |
70 | // | |
71 | void SetRangeRefQ (Int_t nBins, Float_t xMin, Float_t xMax){ fNbinsQ = nBins; fXminQ = xMin; fXmaxQ = xMax; } //Set range for Q reference histograms | |
72 | void SetRangeRefT0 (Int_t nBins, Float_t xMin, Float_t xMax){ fNbinsT0 = nBins; fXminT0 = xMin; fXmaxT0 = xMax; } //Set range for T0 reference histograms | |
73 | void SetRangeRefRMS(Int_t nBins, Float_t xMin, Float_t xMax){ fNbinsRMS = nBins; fXminRMS = xMin; fXmaxRMS = xMax; } //Set range for T0 reference histograms | |
74 | // | |
7fb602b1 | 75 | void SetTimeStampEvent(Double_t timestamp){ fTimeStamp = timestamp; } |
76 | void SetRunNumber(Double_t eventnumber){ fRunNumber = eventnumber; } | |
77 | ||
78 | void SetEventInfo(Double_t runNumber, Double_t timestamp, Double_t eventId){ fRunNumber=runNumber; fTimeStamp=timestamp; fEventId=eventId;} | |
75d8233f | 79 | |
80 | void SetOldRCUformat(Bool_t format=kTRUE){ fOldRCUformat = format; } | |
81 | ||
82 | void SetDebugLevel(Short_t debug=1){ fDebugLevel = debug;} | |
83 | ||
bf57d87d | 84 | void SetPedestalDatabase(AliTPCCalPad *pedestalTPC, AliTPCCalPad *padNoiseTPC) {fPedestalTPC = pedestalTPC; fPadNoiseTPC = padNoiseTPC;} |
85 | ||
75d8233f | 86 | Int_t GetFirstTimeBin() const { return fFirstTimeBin; } |
87 | Int_t GetLastTimeBin() const { return fLastTimeBin; } | |
88 | ||
89 | Int_t GetNeventsProcessed() { return fNevents; } | |
90 | ||
7fb602b1 | 91 | void Merge(AliTPCCalibCE *ce); |
92 | ||
75d8233f | 93 | TGraph *MakeGraphTimeCE(Int_t sector, Int_t xVariable=0, Int_t fitType=0, Int_t fitParameter=0); |
94 | ||
95 | void DumpToFile(const Char_t *filename, const Char_t *dir="", Bool_t append=kFALSE); | |
96 | ||
bf57d87d | 97 | private: |
75d8233f | 98 | Int_t fFirstTimeBin; // First Time bin needed for analysis |
99 | Int_t fLastTimeBin; // Last Time bin needed for analysis | |
100 | ||
101 | // reference histogram ranges | |
102 | Int_t fNbinsT0; // Number of bins for T0 reference histogram | |
103 | Float_t fXminT0; // xmin of T0 reference histogram | |
104 | Float_t fXmaxT0; // xmax of T0 reference histogram | |
105 | Int_t fNbinsQ; // Number of bins for T0 reference histogram | |
106 | Float_t fXminQ; // xmin of T0 reference histogram | |
107 | Float_t fXmaxQ; // xmax of T0 reference histogram | |
108 | Int_t fNbinsRMS; // Number of bins for T0 reference histogram | |
109 | Float_t fXminRMS; // xmin of T0 reference histogram | |
110 | Float_t fXmaxRMS; // xmax of T0 reference histogram | |
111 | ||
112 | Int_t fLastSector; //! Last sector processed | |
113 | ||
114 | Bool_t fOldRCUformat; //! Should we use the old RCU format for data reading | |
115 | ||
116 | AliTPCROC *fROC; //! ROC information | |
117 | AliTPCParam *fParam; //! TPC information | |
118 | ||
119 | AliTPCCalPad *fPedestalTPC; //! Pedestal Information whole TPC | |
120 | AliTPCCalPad *fPadNoiseTPC; //! Pad noise Information whole TPC | |
121 | AliTPCCalROC *fPedestalROC; //! Pedestal Information for current ROC | |
122 | AliTPCCalROC *fPadNoiseROC; //! Pad noise Information for current ROC | |
75d8233f | 123 | |
124 | TObjArray fCalRocArrayT0; // Array of AliTPCCalROC class for Time0 calibration | |
125 | TObjArray fCalRocArrayQ; // Array of AliTPCCalROC class for Charge calibration | |
126 | TObjArray fCalRocArrayRMS; // Array of AliTPCCalROC class for signal width calibration | |
127 | TObjArray fCalRocArrayOutliers; // Array of AliTPCCalROC class for signal outliers | |
128 | ||
129 | TObjArray fHistoQArray; // Calibration histograms for Charge distribution | |
130 | TObjArray fHistoT0Array; // Calibration histograms for Time0 distribution | |
131 | TObjArray fHistoRMSArray; // Calibration histograms for signal width distribution | |
132 | ||
133 | TObjArray fHistoTmean; //! Calibration histograms of the mean CE position for all sectors | |
134 | ||
75d8233f | 135 | TObjArray fParamArrayEventPol1; // Store mean arrival time parameters for each sector event by event from global plane fit |
136 | TObjArray fParamArrayEventPol2; // Store mean arrival time parameters for each sector event by event from global parabola fit | |
bf57d87d | 137 | TObjArray fTMeanArrayEvent; // Store mean arrival time for each sector event by event |
138 | TObjArray fQMeanArrayEvent; // Store mean arrival Charge for each sector event by event | |
75d8233f | 139 | TVectorD fVEventTime; // Timestamps of the events |
140 | TVectorD fVEventNumber; // Eventnumbers of the events | |
141 | Int_t fNevents; // Event counter | |
7fb602b1 | 142 | Double_t fTimeStamp; //! Timestamp of the current event |
143 | Double_t fEventId; //! Event Id of the current event | |
144 | Double_t fRunNumber; //! Run Number of the current event | |
145 | Double_t fOldRunNumber; //! Old Run Number | |
75d8233f | 146 | |
147 | TObjArray fPadTimesArrayEvent; //! Pad Times for the event, before mean Time0 corrections | |
148 | TObjArray fPadQArrayEvent; //! Charge for the event, only needed for debugging streamer | |
149 | TObjArray fPadRMSArrayEvent; //! Signal width for the event, only needed for debugging streamer | |
150 | TObjArray fPadPedestalArrayEvent; //! Signal width for the event, only needed for debugging streamer | |
151 | ||
152 | Int_t fCurrentChannel; //! current channel processed | |
153 | Int_t fCurrentSector; //! current sector processed | |
154 | Int_t fCurrentRow; //! current row processed | |
155 | Float_t fMaxPadSignal; //! maximum bin of current pad | |
156 | Int_t fMaxTimeBin; //! time bin with maximum value | |
157 | TVectorF fPadSignal; //! signal of current Pad | |
158 | Float_t fPadPedestal; //! Pedestal Value of current pad | |
159 | Float_t fPadNoise; //! Noise Value of current pad | |
160 | ||
161 | TVectorD fVTime0Offset; //! Time0 Offset for each sector; | |
162 | TVectorD fVTime0OffsetCounter; //! Time0 Offset counter for each sector; | |
bf57d87d | 163 | TVectorD fVMeanQ; //! Mean Q for each sector; |
164 | TVectorD fVMeanQCounter; //! Mean Q counter for each sector; | |
75d8233f | 165 | //debugging |
166 | Int_t fEvent; | |
167 | TTreeSRedirector *fDebugStreamer; //! debug streamer | |
168 | ||
169 | Short_t fDebugLevel; | |
170 | //! debugging | |
171 | ||
172 | void FindPedestal(Float_t part=.6); | |
173 | void FindCESignal(TVectorD ¶m, Float_t &qSum, const TVectorF maxima); | |
174 | void FindLocalMaxima(TVectorF &maxima); | |
175 | Bool_t IsPeak(Int_t pos, Int_t tminus, Int_t tplus); | |
176 | ||
177 | TH2S* GetHisto(Int_t sector, TObjArray *arr, | |
178 | Int_t nbinsY, Float_t ymin, Float_t ymax, | |
179 | Char_t *type, Bool_t force); | |
180 | TH1S* GetHisto(Int_t sector, TObjArray *arr, | |
181 | Char_t *type, Bool_t force); | |
182 | ||
183 | AliTPCCalROC* GetCalRoc(Int_t sector, TObjArray* arr, Bool_t force); | |
184 | ||
7fb602b1 | 185 | TVectorF* GetVectSector(Int_t sector, TObjArray *arr, UInt_t size, Bool_t force=kFALSE); |
75d8233f | 186 | TVectorF* GetPadTimesEvent(Int_t sector, Bool_t force=kFALSE); |
187 | ||
188 | TObjArray* GetParamArray(Int_t sector, TObjArray *arr, Bool_t force=kFALSE); | |
75d8233f | 189 | |
190 | void ResetEvent(); | |
191 | void ResetPad(); | |
192 | void ProcessPad(); | |
193 | void EndEvent(); | |
194 | ||
195 | ||
196 | //debug | |
75d8233f | 197 | TVectorF* GetPadQEvent(Int_t sector, Bool_t force=kFALSE); |
198 | TVectorF* GetPadRMSEvent(Int_t sector, Bool_t force=kFALSE); | |
199 | TVectorF* GetPadPedestalEvent(Int_t sector, Bool_t force=kFALSE); | |
200 | ||
201 | ||
202 | public: | |
203 | ||
204 | ||
7fb602b1 | 205 | ClassDef(AliTPCCalibCE,2) //Implementation of the TPC Central Electrode calibration |
206 | ||
75d8233f | 207 | }; |
208 | ||
209 | ||
210 | ||
211 | #endif | |
212 |