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
9 ////////////////////////////////////////////////
10 // MUON Cluster Finder Class //
11 ////////////////////////////////////////////////
17 class AliMUONClusterInput;
18 class AliMUONHitMapA1;
19 class AliSegmentation;
20 class AliMUONGeometrySegmentation;
21 class AliMUONRawCluster;
25 class AliMUONClusterFinderVS : public TObject
28 AliMUONClusterFinderVS();
29 virtual ~AliMUONClusterFinderVS();
31 virtual void SetDeclusterFlag(Int_t flag=1) {fDeclusterFlag =flag;}
32 // Set max. cluster size ; bigger clusters will deconvoluted
33 virtual void SetClusterSize(Int_t clsize=5) {fClusterSize = clsize;}
34 // Set max. number of pads per local cluster
35 virtual void SetNperMax(Int_t npermax=5) {fNperMax = npermax;}
36 // Search for raw clusters
37 virtual void FindRawClusters();
39 virtual void FindCluster(Int_t i, Int_t j, Int_t cath, AliMUONRawCluster &c);
41 virtual void Decluster(AliMUONRawCluster *cluster);
42 // Perform split by local maxima
43 virtual void SplitByLocalMaxima(AliMUONRawCluster *cluster);
44 virtual void FindLocalMaxima(AliMUONRawCluster *cluster);
45 virtual void Split(AliMUONRawCluster * cluster);
46 // Perform Double Mathieson Fit
47 Bool_t DoubleMathiesonFit(AliMUONRawCluster *c, Int_t cath);
48 Float_t CombiDoubleMathiesonFit(AliMUONRawCluster *c);
49 Float_t SingleMathiesonFit(AliMUONRawCluster *c, Int_t cath);
50 Float_t CombiSingleMathiesonFit(AliMUONRawCluster *c);
51 // Build up full cluster information
52 virtual void FillCluster(AliMUONRawCluster *cluster, Int_t flag, Int_t cath);
53 virtual void FillCluster(AliMUONRawCluster *cluster, Int_t cath);
54 virtual void FillCluster(AliMUONRawCluster *cluster) {FillCluster(cluster,1,0);}
55 // Add a new raw cluster
56 virtual void AddRawCluster(AliMUONRawCluster& cluster);
57 // Set tracks for debugging
58 virtual void SetTracks(Int_t t1, Int_t t2) {fTrack[0]=t1; fTrack[1]=t2;}
60 void SetDebugLevel(Int_t level) {fDebugLevel = level;}
61 void SetGhostChi2Cut(Float_t cut) {fGhostChi2Cut = cut;}
62 // get raw cluster pointer
63 TClonesArray* GetRawClusters() {return fRawClusters;}
65 void ResetRawClusters();
68 AliMUONClusterFinderVS(const AliMUONClusterFinderVS& clusterFinder);
69 // Assignment operator
70 AliMUONClusterFinderVS & operator = (const AliMUONClusterFinderVS& rhs);
72 AliMUONClusterInput* fInput; // ! AliMUONClusterInput instance
73 AliMUONHitMapA1* fHitMap[2]; // Hit Maps for cathode 1 and 2
74 AliSegmentation* fSeg[2]; // Segmentations for cathode 1 and 2
75 AliMUONGeometrySegmentation* fSeg2[2]; // New Segmentations for cathode 1 and 2
79 Int_t fDeclusterFlag; // flag for declusterin
80 Int_t fClusterSize; // cluster size
81 Int_t fNperMax; // Maximum number of pads per peak
82 Float_t fGhostChi2Cut; // Cut in charge matching chi2
83 // (2 degrees of freedom)
84 // Used by ghost removal
85 // Current decluster result
86 Int_t fMul[2]; // current multiplicity
87 Int_t fNPeaks; // number of local maxima
88 Int_t fNRawClusters; // Number of Raw Clusters
89 TClonesArray* fRawClusters; // array of cluster per ch.
92 AliMUONDigit* fDig[100][2]; // current list of digits
93 Int_t fIx[100][2]; // current list of x-pad-coord.
94 Int_t fIy[100][2]; // current list of y-pad-coord.
95 Float_t fX[100][2]; // current list of x-coord.
96 Float_t fY[100][2]; // current list of y-coord.
97 Float_t fZ[100][2]; // current list of z-coord.
98 Int_t fIndLocal[100][2]; // indices of local maxima
99 Int_t fNLocal[2]; // Number of local maxima
100 Int_t fQ[100][2]; // current list of charges
101 Float_t fZPlane; // currenz z-plane position
102 Int_t fSector; // current sector
105 Double_t fXFit[2]; // x-coordinate
106 Double_t fYFit[2]; // y-coordinate
107 Double_t fQrFit[2]; // charge ratio
108 Float_t fChi2[2]; // chi2 of fit
109 Float_t fXInit[2]; // start values
110 Float_t fYInit[2]; // start values
111 Float_t fQrInit[2]; // start values
112 Int_t fFitStat; // status of fit
114 // Selected track for debugging
115 Int_t fTrack[2]; // Only digits with main contributions from these tracks are
117 Int_t fDebugLevel; // prinout control
118 Int_t fSegmentationType;// new seg
120 // Return pointer to raw clusters
121 ClassDef(AliMUONClusterFinderVS,2) //Class for clustering and reconstruction of space points