]>
Commit | Line | Data |
---|---|---|
83dbc640 | 1 | #ifndef ALIMUONTRACKK_H |
2 | #define ALIMUONTRACKK_H | |
3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
4 | * See cxx source for full Copyright notice */ | |
5 | ||
30178c30 | 6 | /*$Id$*/ |
7 | // Revision of includes 07/05/2004 | |
8 | ||
692de412 | 9 | /// \ingroup rec |
090443dd | 10 | /// \class AliMUONTrackK |
692de412 | 11 | /// \brief Kalman track in MUON arm of ALICE |
585d648a | 12 | /// |
7945aae7 | 13 | /// \author Alexander Zinchenko, JINR Dubna |
692de412 | 14 | |
cc87ebcd | 15 | class AliMUONEventRecoCombi; |
243b68c3 | 16 | class AliMUONHitForRec; |
208f139e | 17 | class AliMUONObjectPair; |
2457f726 | 18 | class AliMUONTrackReconstructorK; |
243b68c3 | 19 | #include "AliMUONTrack.h" |
83dbc640 | 20 | |
4b9e83e9 | 21 | class TArrayD; |
22 | class TClonesArray; | |
23 | class TObjArray; | |
24 | ||
25 | #include <TMatrixD.h> | |
4b9e83e9 | 26 | |
b035a148 | 27 | class AliMUONTrackK : public AliMUONTrack { |
83dbc640 | 28 | |
29 | public: | |
30 | ||
31 | AliMUONTrackK(); // Default constructor | |
32 | virtual ~AliMUONTrackK(); // Destructor | |
83dbc640 | 33 | |
2457f726 | 34 | AliMUONTrackK(AliMUONTrackReconstructorK *TrackReconstructor, TClonesArray *hitForRec); // Constructor |
208f139e | 35 | AliMUONTrackK(AliMUONObjectPair *segment); // Constructor from a segment |
83dbc640 | 36 | |
37 | // Pointer to hits on track | |
71a2d3aa | 38 | TObjArray* GetTrackHits(void) const {return fTrackHits;} ///< ptr. to hits on track |
39 | Int_t GetNTrackHits(void) const {return fNmbTrackHits;} ///< hits on track | |
40 | Double_t GetTrackQuality(void) const {return fChi2;} ///< track quality | |
41 | TMatrixD* GetTrackParameters(void) const {return fTrackPar;} ///< track parameters | |
42 | Double_t GetZ(void) const {return fPosition;} ///< Z-coordinate of track | |
43 | TMatrixD* GetCovariance(void) const {return fCovariance;} ///< covariance matrix | |
44 | Int_t GetTrackDir(void) const {return fTrackDir;} ///< get track propagation direction | |
45 | void SetTrackDir(Int_t iDir) {fTrackDir = iDir;} ///< set track propagation direction | |
46 | Bool_t GetBPFlag(void) const {return fBPFlag;} ///< get backpropagation flag | |
47 | void SetBPFlag(Bool_t BPFlag) {fBPFlag = BPFlag;} ///< set backpropagation flag | |
48 | Int_t GetRecover(void) const {return fRecover;} ///< return recover flag | |
49 | void SetRecover(Int_t iRecover) {fRecover = iRecover;} ///< set recover flag | |
50 | AliMUONObjectPair* GetStartSegment(void) const {return fStartSegment;} ///< return seed segment | |
51 | Bool_t KalmanFilter(Int_t ichamBeg, Int_t ichamEnd, Bool_t Back, Double_t zDipole1, Double_t zDipole2); ///< Kalman filter | |
83dbc640 | 52 | void StartBack(void); // start backpropagator |
53 | void SetTrackQuality(Int_t iChi2); // compute track quality or Chi2 | |
30178c30 | 54 | Bool_t KeepTrack(AliMUONTrackK* track0) const; // keep or discard track |
83dbc640 | 55 | void Kill(void); // kill track candidate |
b035a148 | 56 | Bool_t Smooth(void); // apply smoother |
57 | Double_t GetChi2PerPoint(Int_t iPoint) const; // return Chi2 at point | |
58 | void Print(FILE *lun) const; // print track information | |
71a2d3aa | 59 | void Print(const char* /*opt*/) const {return;} ///< print track information |
b035a148 | 60 | AliMUONHitForRec* GetHitLastOk(void); // get hit before the skipped one |
61 | Int_t GetStation0(void); // return seed station number | |
71a2d3aa | 62 | Int_t DebugLevel(void) const {return fgDebug;} ///< return debug level |
63 | void SetDebugLevel(Int_t iDebug) {fgDebug = iDebug;} ///< set debug level | |
b035a148 | 64 | void FillMUONTrack(void); // set track parameters as for AliMUONTrack |
65 | void SetTrackParam(AliMUONTrackParam *trackParam, TMatrixD *par, Double_t z); // fill AliMUONTrackParam object | |
83dbc640 | 66 | |
71a2d3aa | 67 | /// What is necessary for sorting TClonesArray's |
83dbc640 | 68 | Bool_t IsSortable() const { return kTRUE; } |
69 | Int_t Compare(const TObject* trackK) const; // "Compare" function for sorting | |
70 | ||
71 | ||
72 | ||
73 | protected: | |
74 | ||
30178c30 | 75 | AliMUONTrackK& operator=(const AliMUONTrackK& source); // assignment operator |
76 | ||
83dbc640 | 77 | private: |
78 | ||
829425a5 | 79 | static Int_t fgDebug; ///< debug level |
80 | static Int_t fgNOfPoints; ///< number of points in event | |
81 | //static AliMUON *fgMUON; ///< pointer to MUON module | |
2457f726 | 82 | static AliMUONTrackReconstructorK *fgTrackReconstructor; ///< pointer to event reconstructor |
829425a5 | 83 | static TClonesArray *fgHitForRec; ///< pointer to hits |
84 | static AliMUONEventRecoCombi *fgCombi; ///< pointer to combined cluster/track finder | |
85 | ||
208f139e | 86 | AliMUONObjectPair *fStartSegment; ///< seed segment |
829425a5 | 87 | Double_t fPosition; ///< Z-coordinate of track |
88 | Double_t fPositionNew; //!< Z-coordinate of track | |
89 | Double_t fChi2; ///< Chi2 of track | |
90 | TObjArray *fTrackHits; ///< pointer to hits on track | |
91 | Int_t fNmbTrackHits; ///< number of points on track | |
92 | Int_t fTrackDir; ///< track direction (+(-) towards high (low) z) | |
93 | Bool_t fBPFlag; ///< backpropagation flag (TRUE if backpropagation) | |
94 | Int_t fRecover; ///< recover flag (!=0 if recovery procedure was applied) | |
95 | AliMUONHitForRec *fSkipHit; ///< hit to skip during recovery procedure | |
96 | ||
97 | TMatrixD *fTrackPar; ///< track parameters | |
98 | TMatrixD *fTrackParNew; //!< track parameters | |
99 | TMatrixD *fCovariance; ///< covariance matrix | |
100 | TMatrixD *fWeight; //!< weight matrix (inverse of covariance) | |
83dbc640 | 101 | |
b035a148 | 102 | // For smoother |
829425a5 | 103 | TObjArray *fParExtrap; //!< extrapolated track parameters |
104 | TObjArray *fParFilter; //!< filtered track parameters | |
105 | TObjArray *fParSmooth; //!< smoothed track parameters | |
b035a148 | 106 | |
829425a5 | 107 | TObjArray *fCovExtrap; //!< extrapolated covariance matrices |
108 | TObjArray *fCovFilter; //!< filtered covariance matrices | |
b035a148 | 109 | |
829425a5 | 110 | TObjArray *fJacob; //!< Jacobian matrices |
5017ea54 | 111 | Int_t fNSteps; //!< number of the track propagation points |
112 | TArrayD *fSteps; //!< Z-coordinates of the track propagation points | |
113 | TArrayD *fChi2Array; //!< measurements' contributions to the track Chi2 | |
114 | TArrayD *fChi2Smooth;//!< measurements' contributions to the smoothed track Chi2 | |
b035a148 | 115 | |
83dbc640 | 116 | // Functions |
117 | ||
71a2d3aa | 118 | AliMUONTrackK (const AliMUONTrackK& source); ///< Not implemented copy constructor |
83dbc640 | 119 | void EvalCovariance(Double_t dZ); |
120 | void ParPropagation(Double_t zEnd); | |
b035a148 | 121 | void WeightPropagation(Double_t zEnd, Bool_t smooth); |
83dbc640 | 122 | void MSThin(Int_t sign); |
123 | void MSLine(Double_t dZ, Double_t X0); | |
cc87ebcd | 124 | Bool_t FindPoint(Int_t ichamb, Double_t zEnd, Int_t currIndx, Int_t iFB, AliMUONHitForRec *&hitAdd, Int_t iz); |
83dbc640 | 125 | void TryPoint(TMatrixD &point, const TMatrixD &pointWeight, TMatrixD &trackParTmp, Double_t &dChi2); |
126 | void SetGeantParam(Double_t *VGeant3, Int_t iFB); | |
127 | void GetFromGeantParam(Double_t *VGeant3, Int_t iFB); | |
b035a148 | 128 | Bool_t Recover(void); |
129 | void AddMatrices(AliMUONTrackK *trackK, Double_t dChi2, AliMUONHitForRec *hitAdd); | |
1ffe5558 | 130 | void CreateMatrix(TObjArray *objArray) const; |
b035a148 | 131 | void RemoveMatrices(Double_t zEnd); |
132 | void RemoveMatrices(AliMUONTrackK* trackK); | |
133 | void Outlier(); | |
134 | void SortHits(Int_t iflag, TObjArray *array); | |
135 | void DropBranches(Int_t imax, TObjArray *hits); | |
208f139e | 136 | void DropBranches(AliMUONObjectPair *segment); |
b035a148 | 137 | Bool_t ExistDouble(AliMUONHitForRec *hit); |
138 | Bool_t ExistDouble(void); | |
4fcb5a03 | 139 | void CheckBranches(TArrayD &branchChi2, Int_t nBranch); |
83dbc640 | 140 | |
343146bf | 141 | private: |
142 | // Some constants | |
829425a5 | 143 | static const Int_t fgkSize; ///< number of track parameters |
144 | static const Int_t fgkNSigma; ///< 4; acceptance window width in sigmas | |
145 | static const Double_t fgkChi2max; ///< 25; chi2 cut in smoother for outlier detection | |
146 | static const Int_t fgkTriesMax; ///< max number of attempts to find exact position during tracking | |
147 | static const Double_t fgkEpsilon; ///< tracking precision (cm) | |
343146bf | 148 | |
83dbc640 | 149 | ClassDef(AliMUONTrackK,0) // Kalman track in MUON arm of ALICE |
150 | }; | |
151 | #endif |