source: trunk/MagicSoft/Mars/mcalib/MCalibrationPix.h@ 3236

Last change on this file since 3236 was 3236, checked in by gaug, 21 years ago
*** empty log message ***
File size: 10.7 KB
Line 
1#ifndef MARS_MCalibrationPix
2#define MARS_MCalibrationPix
3
4#ifndef MARS_MParContainer
5#include "MParContainer.h"
6#endif
7
8#include "MHCalibrationPixel.h"
9
10class MCalibrationPix : public MParContainer
11{
12private:
13
14 static const Float_t gkAverageQE; // The average quantum efficieny agreed on for the first analysis
15 static const Float_t gkAverageQEErr; // The error of average quantum efficieny
16
17 static const Float_t gkConversionHiLo; // The default conversion factor HI gain - Lo Gain
18 static const Float_t gkConversionHiLoError; // The error of the default conversion factor
19
20 static const Float_t gkElectronicPedRms; // The pure electronic component of the RMS
21 static const Float_t gkErrElectronicPedRms; // The error of the pure electronic component of the RMS
22 static const Float_t gkFFactor; // The laboratory F-factor of the PMTs
23 static const Float_t gkFFactorError; // The laboratory F-factor Error of the PMTs
24 static const Float_t gkChargeLimit; // The limit (in units of PedRMS) for acceptance of the fitted mean charge
25 static const Float_t gkChargeErrLimit; // The limit (in units of PedRMS) for acceptance of the fitted charge sigma
26 static const Float_t gkChargeRelErrLimit; // The limit (in units of Error of fitted charge) for acceptance of the fitted mean
27 static const Float_t gkTimeLimit; // The limit (in units of FADC slices) for acceptance of the fitted time
28 static const Float_t gkTimeErrLimit; // The limit (in units of FADC slices) for acceptance of the fitted time sigma
29 static const Float_t gkConvFFactorRelErrorLimit; // The limit (in units of [1]) for acceptance of the rel. error of the conversion factor with the FFactor method
30
31 Int_t fPixId; // the pixel Id
32
33 UInt_t fFlags; // Flag for the set Bits
34
35 Float_t fAverageQE;
36 Float_t fAverageQEErr;
37
38 Float_t fCharge; // The mean reduced charge after the fit
39 Float_t fErrCharge; // The error of reduced mean charge after the fit
40 Float_t fSigmaCharge; // The sigma of the mean charge after the fit
41 Float_t fErrSigmaCharge; // The error of the sigma of the mean charge after the fit
42 Float_t fRSigmaCharge; // The reduced squares of sigmas after the fit
43 Float_t fErrRSigmaCharge; // The reduced squares of sigmas after the fit
44 Float_t fChargeProb; // The probability of the fit function
45
46 Float_t fPed; // The mean pedestal (from MPedestalPix)
47 Float_t fPedRms; // The pedestal RMS (from MPedestalPix)
48 Float_t fErrPedRms; // The error of the pedestal RMS (from MPedestalPix)
49
50 Float_t fAbsTimeMean;
51 Float_t fAbsTimeMeanErr;
52 Float_t fAbsTimeRms;
53
54 Byte_t fTimeFirstHiGain; // The first used FADC slice
55 Byte_t fTimeLastHiGain; // The last used FADC slice
56
57 Byte_t fTimeFirstLoGain; // The first used FADC slice
58 Byte_t fTimeLastLoGain; // The last used FADC slice
59
60 Float_t fPheFFactorMethod; // The number of Phe's calculated (F-factor method)
61 Float_t fPheFFactorMethodError; // The error on the number of Phe's calculated (F-factor method)
62
63 Float_t fMeanConversionFFactorMethod; // The conversion factor to Phe's (F-factor method)
64 Float_t fMeanConversionBlindPixelMethod; // The conversion factor to Ph's (Blind Pixel method)
65 Float_t fMeanConversionPINDiodeMethod; // The conversion factor to Ph's (PIN Diode method)
66 Float_t fMeanConversionCombinedMethod; // The conversion factor to Ph's (all methods combined)
67
68 Float_t fErrorConversionFFactorMethod; // The error of the conversion factor to Phe's (F-factor method)
69 Float_t fErrorConversionBlindPixelMethod; // The error of the conversion factor to Ph's (Blind Pixel method)
70 Float_t fErrorConversionPINDiodeMethod; // The error of the conversion factor to Ph's (PIN Diode method)
71 Float_t fErrorConversionCombinedMethod; // The error of the conversion factor to Ph's (all methods combined)
72
73 Float_t fSigmaConversionFFactorMethod; // The sigma of conversion factor to Phe's (F-factor method)
74 Float_t fSigmaConversionBlindPixelMethod; // The conversion factor to Ph's (Blind Pixel method)
75 Float_t fSigmaConversionPINDiodeMethod; // The conversion factor to Ph's (PIN Diode method)
76 Float_t fSigmaConversionCombinedMethod; // The conversion factor to Ph's (all methods combined)
77
78 Float_t fTotalFFactor; // The F-Factor of the total readout system (Sigma(out)/mean(out)*Mean(in)/sigma(in)
79 Float_t fTotalFFactorError; // The error on the F-Factor of the total readout system
80
81 Float_t fConversionHiLo; // The conversion factor between Hi Gain and Lo Gain
82 Float_t fConversionHiLoError; // The error of the conversion factor between Hi Gain and Lo Gain
83
84 Float_t fNumHiGainSamples;
85 Float_t fNumLoGainSamples;
86
87 Bool_t fFactorCalculated;
88
89 enum { kHiGainSaturation,
90 kExcluded, kExcludeQualityCheck,
91 kChargeValid, kTimeFitValid,
92 kFitted, kOscillating,
93 kBlindPixelMethodValid, kFFactorMethodValid, kPINDiodeMethodValid, kCombinedMethodValid };
94
95 MHCalibrationPixel *fHist; // Pointer to the histograms performing the fits, etc.
96
97 Bool_t CheckChargeValidity();
98 Bool_t CheckTimeFitValidity();
99
100public:
101
102 MCalibrationPix(const char *name=NULL, const char *title=NULL);
103 ~MCalibrationPix();
104
105 void Clear(Option_t *o="");
106
107 // Getter
108 MHCalibrationPixel *GetHist() const { return fHist; }
109
110 // Setter
111 void SetPedestal(const Float_t ped, const Float_t pedrms,
112 const Float_t higainsamp, const Float_t logainsamp);
113 void SetConversionHiLo( const Float_t c=gkConversionHiLo) { fConversionHiLo = c; }
114 void SetConversionHiLoError(const Float_t e=gkConversionHiLoError) { fConversionHiLoError = e; }
115 void SetAverageQE(const Float_t qe=gkAverageQE, const Float_t err=gkAverageQEErr)
116 { fAverageQE = qe; fAverageQEErr = err; }
117
118 // Setters for MC
119 void SetConversionFFactorMethod(Float_t c, Float_t err, Float_t sig);
120 void SetConversionBlindPixelMethod(Float_t c, Float_t err, Float_t sig);
121 void SetConversionPINDiodeMethod(Float_t c, Float_t err, Float_t sig);
122 void SetConversionCombinedMethod(Float_t c, Float_t err, Float_t sig);
123
124 // Bit Setters
125 void SetHiGainSaturation(Bool_t b = kTRUE);
126 void SetExcluded(Bool_t b = kTRUE);
127 void SetExcludeQualityCheck(Bool_t b = kTRUE);
128 void SetChargeValid(Bool_t b = kTRUE);
129 void SetFitted(Bool_t b = kTRUE);
130 void SetOscillating(Bool_t b = kTRUE);
131 void SetBlindPixelMethodValid(Bool_t b = kTRUE);
132 void SetFFactorMethodValid(Bool_t b = kTRUE);
133 void SetPINDiodeMethodValid(Bool_t b = kTRUE);
134 void SetAbsTimeBordersHiGain(Byte_t f, Byte_t l);
135 void SetAbsTimeBordersLoGain(Byte_t f, Byte_t l);
136
137 // Charges
138 Float_t GetCharge() const { return fCharge; }
139 Float_t GetErrCharge() const { return fErrCharge; }
140 Float_t GetChargeProb() const { return fChargeProb; }
141 Float_t GetSigmaCharge() const { return fSigmaCharge; }
142 Float_t GetErrSigmaCharge() const { return fErrSigmaCharge; }
143 Float_t GetRSigmaCharge() const { return fRSigmaCharge; }
144 Float_t GetErrRSigmaCharge() const { return fErrRSigmaCharge; }
145
146
147 Float_t GetAbsTimeMean() const { return fAbsTimeMean; }
148 Float_t GetAbsTimeMeanErr() const { return fAbsTimeMeanErr; }
149 Float_t GetAbsTimeRms() const { return fAbsTimeRms; }
150
151 // Conversion Factors
152 Float_t GetConversionHiLo() const { return fConversionHiLo; }
153 Float_t GetConversionHiLoError() const { return fConversionHiLoError; }
154
155 Float_t GetMeanConversionBlindPixelMethod() const { return fMeanConversionBlindPixelMethod ; }
156 Float_t GetErrorConversionBlindPixelMethod() const { return fErrorConversionBlindPixelMethod ; }
157 Float_t GetSigmaConversionBlindPixelMethod() const { return fSigmaConversionBlindPixelMethod ; }
158
159 Float_t GetMeanConversionFFactorMethod();
160 Float_t GetErrorConversionFFactorMethod();
161 Float_t GetSigmaConversionFFactorMethod();
162
163 Float_t GetMeanConversionPINDiodeMethod() const { return fMeanConversionPINDiodeMethod ; }
164 Float_t GetErrorConversionPINDiodeMethod() const { return fErrorConversionPINDiodeMethod ; }
165 Float_t GetSigmaConversionPINDiodeMethod() const { return fSigmaConversionPINDiodeMethod ; }
166
167 Float_t GetMeanConversionCombinedMethod() const { return fMeanConversionCombinedMethod ; }
168 Float_t GetErrorConversionCombinedMethod() const { return fErrorConversionCombinedMethod ; }
169 Float_t GetSigmaConversionCombinedMethod() const { return fSigmaConversionCombinedMethod ; }
170
171 Float_t GetPheFFactorMethod();
172 Float_t GetPheFFactorMethodError();
173
174 Int_t GetPixId() const { return fPixId; }
175
176 Float_t GetPed() const { return fPed; }
177 Float_t GetPedRms() const { return fPedRms; }
178
179 Float_t GetTotalFFactorFFactorMethod();
180 Float_t GetTotalFFactorErrorFFactorMethod();
181
182 Float_t GetTotalFFactorBlindPixelMethod();
183 Float_t GetTotalFFactorErrorBlindPixelMethod();
184
185 Float_t GetTotalFFactorPINDiodeMethod();
186 Float_t GetTotalFFactorErrorPINDiodeMethod();
187
188 Float_t GetTotalFFactorCombinedMethod();
189 Float_t GetTotalFFactorErrorCombinedMethod();
190
191 Bool_t IsExcluded() const;
192 Bool_t IsExcludeQualityCheck() const;
193 Bool_t IsHiGainSaturation() const;
194 Bool_t IsChargeValid() const;
195 Bool_t IsFitted() const;
196 Bool_t IsOscillating();
197 Bool_t IsBlindPixelMethodValid() const;
198 Bool_t IsPINDiodeMethodValid() const;
199 Bool_t IsFFactorMethodValid();
200 Bool_t IsCombinedMethodValid();
201
202 // Fill histos
203 Bool_t FillChargeHiGain(Float_t q) const { return fHist->FillChargeHiGain(q); }
204 Bool_t FillAbsTimeHiGain(Float_t t) const { return fHist->FillAbsTimeHiGain(t); }
205
206 Bool_t FillChargeLoGain(Float_t q) const { return fHist->FillChargeLoGain(q); }
207 Bool_t FillAbsTimeLoGain(Float_t t) const { return fHist->FillAbsTimeLoGain(t); }
208
209 Bool_t FillGraphs(Float_t qhi,Float_t qlo) const { return fHist->FillGraphs(qhi,qlo); }
210
211 void DefinePixId(Int_t i);
212
213 // Fits
214 Bool_t FitCharge();
215
216 // Draws
217 void Draw(Option_t *opt="") { fHist->Draw(opt); }
218 TObject *DrawClone(Option_t *opt="") const { return fHist->DrawClone(opt); }
219
220 // Miscellaneous
221 void ApplyLoGainConversion();
222 void CheckOscillations();
223 Bool_t CalcFFactorMethod();
224
225 ClassDef(MCalibrationPix, 1) // Container for Calibration of one pixel
226};
227
228#endif
229
Note: See TracBrowser for help on using the repository browser.