]>
Commit | Line | Data |
---|---|---|
a9e2aefa | 1 | #ifndef ALIMUONCLUSTERFINDERVS_H |
2 | #define ALIMUONCLUSTERFINDERVS_H | |
3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
4 | * See cxx source for full Copyright notice */ | |
5 | ||
6 | /* $Id$ */ | |
30178c30 | 7 | // Revision of includes 07/05/2004 |
a9e2aefa | 8 | |
692de412 | 9 | /// \ingroup rec |
10 | /// \class AliMUONClusterFinderVS | |
11 | /// \brief Class for clustering and reconstruction of space points | |
c1a185bf | 12 | |
30aaba74 | 13 | #include <TObject.h> |
9825400f | 14 | |
30178c30 | 15 | class TClonesArray; |
16 | ||
c1a185bf | 17 | class AliMUONClusterInput; |
5289cf2f | 18 | class AliMUONDigitMapA1; |
fed772f3 | 19 | class AliMUONGeometrySegmentation; |
30aaba74 | 20 | class AliMUONRawCluster; |
21 | class AliMUONDigit; | |
22 | ||
23 | ||
24 | class AliMUONClusterFinderVS : public TObject | |
a9e2aefa | 25 | { |
26 | public: | |
a9e2aefa | 27 | AliMUONClusterFinderVS(); |
4da78c65 | 28 | virtual ~AliMUONClusterFinderVS(); |
30aaba74 | 29 | // Decluster ? |
30 | virtual void SetDeclusterFlag(Int_t flag=1) {fDeclusterFlag =flag;} | |
31 | // Set max. cluster size ; bigger clusters will deconvoluted | |
32 | virtual void SetClusterSize(Int_t clsize=5) {fClusterSize = clsize;} | |
33 | // Set max. number of pads per local cluster | |
34 | virtual void SetNperMax(Int_t npermax=5) {fNperMax = npermax;} | |
a9e2aefa | 35 | // Search for raw clusters |
30aaba74 | 36 | virtual void FindRawClusters(); |
a9e2aefa | 37 | // Find cluster |
38 | virtual void FindCluster(Int_t i, Int_t j, Int_t cath, AliMUONRawCluster &c); | |
39 | // Decluster | |
40 | virtual void Decluster(AliMUONRawCluster *cluster); | |
41 | // Perform split by local maxima | |
42 | virtual void SplitByLocalMaxima(AliMUONRawCluster *cluster); | |
43 | virtual void FindLocalMaxima(AliMUONRawCluster *cluster); | |
44 | virtual void Split(AliMUONRawCluster * cluster); | |
a9e2aefa | 45 | // Perform Double Mathieson Fit |
30aaba74 | 46 | Bool_t DoubleMathiesonFit(AliMUONRawCluster *c, Int_t cath); |
a9e2aefa | 47 | Float_t CombiDoubleMathiesonFit(AliMUONRawCluster *c); |
48 | Float_t SingleMathiesonFit(AliMUONRawCluster *c, Int_t cath); | |
49 | Float_t CombiSingleMathiesonFit(AliMUONRawCluster *c); | |
50 | // Build up full cluster information | |
51 | virtual void FillCluster(AliMUONRawCluster *cluster, Int_t flag, Int_t cath); | |
52 | virtual void FillCluster(AliMUONRawCluster *cluster, Int_t cath); | |
30aaba74 | 53 | virtual void FillCluster(AliMUONRawCluster *cluster) {FillCluster(cluster,1,0);} |
54 | // Add a new raw cluster | |
b137f8b9 | 55 | virtual void AddRawCluster(AliMUONRawCluster& cluster); |
30aaba74 | 56 | // Set tracks for debugging |
57 | virtual void SetTracks(Int_t t1, Int_t t2) {fTrack[0]=t1; fTrack[1]=t2;} | |
07cfabcf | 58 | void SetGhostChi2Cut(Float_t cut) {fGhostChi2Cut = cut;} |
4da78c65 | 59 | // get raw cluster pointer |
60 | TClonesArray* GetRawClusters() {return fRawClusters;} | |
61 | // reset raw clusters | |
62 | void ResetRawClusters(); | |
30aaba74 | 63 | |
64 | protected: | |
30178c30 | 65 | AliMUONClusterFinderVS(const AliMUONClusterFinderVS& clusterFinder); |
66 | // Assignment operator | |
67 | AliMUONClusterFinderVS & operator = (const AliMUONClusterFinderVS& rhs); | |
68 | ||
829425a5 | 69 | AliMUONClusterInput* fInput; //!< AliMUONClusterInput instance |
70 | AliMUONDigitMapA1* fDigitMap[2]; ///< Hit Maps for cathode 1 and 2 | |
71 | AliMUONGeometrySegmentation* fSeg2[2]; ///< New Segmentations for cathode 1 and 2 | |
fed772f3 | 72 | |
f0d86bc4 | 73 | |
30aaba74 | 74 | // Configuration |
829425a5 | 75 | Int_t fDeclusterFlag; ///< flag for declusterin |
76 | Int_t fClusterSize; ///< cluster size | |
77 | Int_t fNperMax; ///< Maximum number of pads per peak | |
78 | Float_t fGhostChi2Cut; ///< \brief Cut in charge matching chi2 | |
79 | /// (2 degrees of freedom) | |
80 | /// Used by ghost removal | |
b137f8b9 | 81 | // Current decluster result |
829425a5 | 82 | Int_t fMul[2]; ///< current multiplicity |
83 | Int_t fNPeaks; ///< number of local maxima | |
84 | Int_t fNRawClusters; ///< Number of Raw Clusters | |
85 | TClonesArray* fRawClusters; ///< array of cluster per ch. | |
4da78c65 | 86 | |
b137f8b9 | 87 | // Local data store |
829425a5 | 88 | AliMUONDigit* fDig[100][2]; ///< current list of digits |
89 | Int_t fIx[100][2]; ///< current list of x-pad-coord. | |
90 | Int_t fIy[100][2]; ///< current list of y-pad-coord. | |
91 | Float_t fX[100][2]; ///< current list of x-coord. | |
92 | Float_t fY[100][2]; ///< current list of y-coord. | |
93 | Float_t fZ[100][2]; ///< current list of z-coord. | |
94 | Int_t fIndLocal[100][2]; ///< indices of local maxima | |
95 | Int_t fNLocal[2]; ///< Number of local maxima | |
96 | Int_t fQ[100][2]; ///< current list of charges | |
97 | Float_t fZPlane; ///< currenz z-plane position | |
98 | Int_t fSector; ///< current sector | |
f0d86bc4 | 99 | |
b137f8b9 | 100 | // Current Fit |
829425a5 | 101 | Double_t fXFit[2]; ///< x-coordinate |
102 | Double_t fYFit[2]; ///< y-coordinate | |
103 | Double_t fQrFit[2]; ///< charge ratio | |
104 | Float_t fChi2[2]; ///< chi2 of fit | |
105 | Float_t fXInit[2]; ///< start values | |
106 | Float_t fYInit[2]; ///< start values | |
107 | Float_t fQrInit[2]; ///< start values | |
108 | Int_t fFitStat; ///< status of fit | |
a9e2aefa | 109 | |
b137f8b9 | 110 | // Selected track for debugging |
829425a5 | 111 | Int_t fTrack[2]; ///< Only digits with main contributions from these tracks are |
b137f8b9 | 112 | // considered |
b137f8b9 | 113 | |
114 | // Return pointer to raw clusters | |
4da78c65 | 115 | ClassDef(AliMUONClusterFinderVS,2) //Class for clustering and reconstruction of space points |
b137f8b9 | 116 | }; |
a9e2aefa | 117 | #endif |
118 | ||
119 | ||
120 | ||
121 | ||
122 | ||
123 | ||
124 | ||
125 | ||
126 | ||
127 | ||
30aaba74 | 128 | |
129 | ||
130 | ||
131 | ||
132 |