1 #ifndef ALITPCKALMANFIT_H
2 #define ALITPCKALMANFIT_H
10 class TTreeSRedirector;
11 class AliTrackPointArray;
18 class AliTPCkalmanFit: public TNamed{
22 void InitTransformation();
23 void Add(const AliTPCkalmanFit * kalman);
25 void AddCalibration(AliTPCTransformation * calib);
26 AliTPCTransformation * GetTransformation(Int_t i){return (fCalibration)? (AliTPCTransformation *)fCalibration->At(i):0;}
28 void SetStatus(const char * name, Bool_t setOn, Bool_t isOr=kTRUE);
30 void FitTrackLinear(AliTrackPointArray& points, Int_t step=1, TTreeSRedirector *debug=0);
31 void DumpTrackLinear(AliTrackPointArray& points, TTreeSRedirector *debug);
32 void UpdateLinear(AliTrackPoint &point, TTreeSRedirector *debug=0);
34 void Propagate(TTreeSRedirector *debug=0);
35 void PropagateTime(Int_t time);
36 void Update(const AliTPCkalmanFit * kalman);
38 static AliTrackPointArray * MakePointArrayLinear(Double_t alpha, Double_t y0, Double_t z0, Double_t ky, Double_t kz, Double_t err=0.02);
39 void ApplyCalibration(AliTrackPointArray *array, Double_t csign);
40 Bool_t CheckCovariance(TMatrixD &covar, Float_t maxEl);
42 Bool_t DumpCorelation(Double_t threshold, const char *mask=0);
43 Bool_t DumpCalib(const char *mask=0);
45 Double_t GetTPCDeltaXYZ(Int_t coord, Int_t volID, Double_t x, Double_t y, Double_t z);
46 static Double_t SGetTPCDeltaXYZ(Int_t coord, Int_t volID, Double_t x, Double_t y, Double_t z);
47 AliTPCkalmanFit *Test(Int_t ntracks);
50 // Calibration parameters
52 TObjArray *fCalibration; // array of calibrations
53 TMatrixD *fCalibParam; // calibration parameters
54 TMatrixD *fCalibCovar; // calibration parameters
58 TMatrixD *fLinearParam; // linear parameters
59 TMatrixD *fLinearCovar; // linear covariance
60 THnSparse *fLinearTrackDelta[8]; // linear tracks matching residuals - delta
61 THnSparse *fLinearTrackPull[8]; // linear tracks matching residuals - pull
65 Int_t fLastTimeStamp; // last time stamp - used for propagation of parameters
66 //static AliTPCkalmanFit* Instance();
67 void SetInstance(AliTPCkalmanFit*param){fgInstance = param;}
68 static AliTPCkalmanFit* fgInstance; //! Instance of this class (singleton implementation)
70 Double_t fCurrentAlpha; //! current rotation frame
71 Double_t fCA; //! cosine of current angle
72 Double_t fSA; //! sinus of current angle
73 // AliTPCkalmanFit& operator=(const AliTPCkalmanFit&){;}// not implemented
74 // AliTPCkalmanFit(const AliTPCkalmanFit&){;} //not implemented
75 ClassDef(AliTPCkalmanFit,1);