11 |
#include <TArrayF.h> |
#include <TArrayF.h> |
12 |
#include <TClonesArray.h> |
#include <TClonesArray.h> |
13 |
|
|
14 |
|
#include <math.h> // EMILIANO |
15 |
|
|
16 |
#include <ToFStruct.h> |
#include <ToFStruct.h> |
17 |
|
|
18 |
|
|
48 |
ToFPMT(const ToFPMT&); |
ToFPMT(const ToFPMT&); |
49 |
// |
// |
50 |
ToFPMT* GetToFPMT(){return this;}; |
ToFPMT* GetToFPMT(){return this;}; |
51 |
void Clear(); |
void Clear(Option_t *t=""); |
52 |
|
|
53 |
|
|
54 |
|
|
84 |
TArrayI tdcflag; ///<flag for artificial TDC, "0" if normal TDC value |
TArrayI tdcflag; ///<flag for artificial TDC, "0" if normal TDC value |
85 |
|
|
86 |
/** |
/** |
87 |
* \brief beta, 12 measurements for the 12 combinations, beta[13] is weighted mean |
* \brief beta, 12 measurements for the 12 combinations, beta[13] is modified weighted mean |
88 |
* |
* |
89 |
* The 12 measurements are S11-S31, S11-S32, S12-S31, S12-S32, and then analogue for |
* The 12 measurements are S11-S31, S11-S32, S12-S31, S12-S32, and then analogue for |
90 |
* S2-S3 and S1-S2. |
* S2-S3 and S1-S2. |
91 |
* In the moment all measurements are taken and the weighted mean is calculated. |
* The calculation of beta[13] is now modified: |
92 |
* Note that the weights are just simple overall results for S1-S3, S2-S3, and S1-S2. |
* We check the individual weights for artificial TDC values, then calculate |
93 |
* Artificial measurments are not treated correct, (since there is only one real |
* am mean beta for the first time. In a second step we loop again through |
94 |
* measurment the weight should be different then for two meassurments). |
* the single measurements, checking for the residual from the mean |
95 |
* The beta calculation will be improved in the next release. |
* The cut on the residual reject measurements > "x"-sigma. A chi2 value is |
96 |
|
* calculated, furthermore a "quality" value by adding the weights which |
97 |
|
* are finally used. If all measurements are taken, "quality" will be = 505. |
98 |
|
* A chi2 cut around 3-4 and a quality-cut > 400 is needed for clean beta |
99 |
|
* The Level2 beta[12] which is derived in the fortran routines uses: 10.,200.,20. |
100 |
|
* This is not a very high quality measurement. One can re-calculate a new beta[13] |
101 |
|
* using the L2-method "CalcBeta" |
102 |
*/ |
*/ |
103 |
Float_t beta[13]; |
Float_t beta[13]; |
104 |
// |
// |
117 |
ToFTrkVar(const ToFTrkVar&); |
ToFTrkVar(const ToFTrkVar&); |
118 |
|
|
119 |
ToFTrkVar* GetToFTrkVar(){return this;}; |
ToFTrkVar* GetToFTrkVar(){return this;}; |
120 |
void Clear(); |
void Clear(Option_t *t=""); |
121 |
|
|
122 |
ClassDef(ToFTrkVar,1); |
ClassDef(ToFTrkVar,1); |
123 |
// |
// |
141 |
Int_t default_calib; ///< one if the default calibration has been used to process the data, zero otherwise |
Int_t default_calib; ///< one if the default calibration has been used to process the data, zero otherwise |
142 |
// |
// |
143 |
Float_t GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl); // gf Apr 07 |
Float_t GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl); // gf Apr 07 |
144 |
|
|
145 |
|
Float_t CalcBeta(Int_t notrack, Float_t resmax, Float_t qualitycut, Float_t chi2cut); // wm feb 08 |
146 |
|
|
147 |
|
// |
148 |
|
// Float_t CalcBeta(Int_t notrack, Float_t resmax, Float_t chi2cut, Float_t qualitycut); // wm feb 08 |
149 |
// |
// |
150 |
// methods to make life simplier during the analysis, returns a pointer to the ToFTrkVar class containing track related variables |
// methods to make life simplier during the analysis, returns a pointer to the ToFTrkVar class containing track related variables |
151 |
// |
// |
175 |
void GetPaddlePlane(Int_t padid, Int_t &plane, Int_t &paddle); |
void GetPaddlePlane(Int_t padid, Int_t &plane, Int_t &paddle); |
176 |
Int_t GetNPaddle(Int_t plane); |
Int_t GetNPaddle(Int_t plane); |
177 |
// |
// |
178 |
|
// |
179 |
|
|
180 |
// |
// |
181 |
// constructor |
// constructor |
182 |
// |
// |
183 |
ToFLevel2(); |
ToFLevel2(); |
184 |
~ToFLevel2(){Delete();}; //ELENA |
~ToFLevel2(){Delete();}; //ELENA |
185 |
void Delete(); //ELENA |
void Delete(Option_t *t=""); //ELENA |
186 |
void Set();//ELENA |
void Set();//ELENA |
187 |
// |
// |
188 |
// |
// |
220 |
Int_t GetToFPlaneIndex(Int_t plane_id); |
Int_t GetToFPlaneIndex(Int_t plane_id); |
221 |
Bool_t HitPaddle(Int_t ,Int_t); |
Bool_t HitPaddle(Int_t ,Int_t); |
222 |
Int_t GetNHitPaddles(Int_t plane); |
Int_t GetNHitPaddles(Int_t plane); |
223 |
void Clear(); |
void Clear(Option_t *t=""); |
224 |
// |
// |
225 |
ClassDef(ToFLevel2,3); |
ClassDef(ToFLevel2,3); |
226 |
}; |
}; |