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 */
7 // Revision of includes 07/05/2004
10 /// \class AliMUONClusterFinderVS
11 /// \brief Class for clustering and reconstruction of space points
17 class AliMUONClusterInput;
18 class AliMUONDigitMapA1;
19 class AliMUONGeometrySegmentation;
20 class AliMUONRawCluster;
24 class AliMUONClusterFinderVS : public TObject
27 AliMUONClusterFinderVS();
28 virtual ~AliMUONClusterFinderVS();
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;}
35 // Search for raw clusters
36 virtual void FindRawClusters();
38 virtual void FindCluster(Int_t i, Int_t j, Int_t cath, AliMUONRawCluster &c);
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);
45 // Perform Double Mathieson Fit
46 Bool_t DoubleMathiesonFit(AliMUONRawCluster *c, Int_t cath);
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);
53 virtual void FillCluster(AliMUONRawCluster *cluster) {FillCluster(cluster,1,0);}
54 // Add a new raw cluster
55 virtual void AddRawCluster(AliMUONRawCluster& cluster);
56 // Set tracks for debugging
57 virtual void SetTracks(Int_t t1, Int_t t2) {fTrack[0]=t1; fTrack[1]=t2;}
58 void SetGhostChi2Cut(Float_t cut) {fGhostChi2Cut = cut;}
59 // get raw cluster pointer
60 TClonesArray* GetRawClusters() {return fRawClusters;}
62 void ResetRawClusters();
66 AliMUONClusterInput* fInput; //!< AliMUONClusterInput instance
67 AliMUONDigitMapA1* fDigitMap[2]; ///< Hit Maps for cathode 1 and 2
68 AliMUONGeometrySegmentation* fSeg2[2]; ///< New Segmentations for cathode 1 and 2
72 Int_t fDeclusterFlag; ///< flag for declusterin
73 Int_t fClusterSize; ///< cluster size
74 Int_t fNperMax; ///< Maximum number of pads per peak
75 Float_t fGhostChi2Cut; ///< \brief Cut in charge matching chi2
76 /// (2 degrees of freedom)
77 /// Used by ghost removal
78 // Current decluster result
79 Int_t fMul[2]; ///< current multiplicity
80 Int_t fNPeaks; ///< number of local maxima
81 Int_t fNRawClusters; ///< Number of Raw Clusters
82 TClonesArray* fRawClusters; ///< array of cluster per ch.
85 AliMUONDigit* fDig[100][2]; ///< current list of digits
86 Int_t fIx[100][2]; ///< current list of x-pad-coord.
87 Int_t fIy[100][2]; ///< current list of y-pad-coord.
88 Float_t fX[100][2]; ///< current list of x-coord.
89 Float_t fY[100][2]; ///< current list of y-coord.
90 Float_t fZ[100][2]; ///< current list of z-coord.
91 Int_t fIndLocal[100][2]; ///< indices of local maxima
92 Int_t fNLocal[2]; ///< Number of local maxima
93 Float_t fQ[100][2]; ///< current list of charges
94 Float_t fZPlane; ///< currenz z-plane position
95 Int_t fSector; ///< current sector
98 Double_t fXFit[2]; ///< x-coordinate
99 Double_t fYFit[2]; ///< y-coordinate
100 Double_t fQrFit[2]; ///< charge ratio
101 Float_t fChi2[2]; ///< chi2 of fit
102 Float_t fXInit[2]; ///< start values
103 Float_t fYInit[2]; ///< start values
104 Float_t fQrInit[2]; ///< start values
105 Int_t fFitStat; ///< status of fit
107 // Selected track for debugging
108 Int_t fTrack[2]; ///< Only digits with main contributions from these tracks are
112 AliMUONClusterFinderVS(const AliMUONClusterFinderVS& clusterFinder);
113 // Assignment operator
114 AliMUONClusterFinderVS & operator = (const AliMUONClusterFinderVS& rhs);
116 ClassDef(AliMUONClusterFinderVS,3) //Class for clustering and reconstruction of space points