MUON/AliMUONClusterFinderVS.h

   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$ */
   7 // Revision of includes 07/05/2004
   8
   9 ////////////////////////////////////////////////
  10 //  MUON Cluster Finder Class                 //
  11 ////////////////////////////////////////////////
  12
  13 #include <TObject.h>
  14
  15 class TClonesArray;
  16
  17 class AliMUONClusterInput;
  18 class AliMUONHitMapA1;
  19 class AliSegmentation;
  20 class AliMUONGeometrySegmentation;
  21 class AliMUONRawCluster;
  22 class AliMUONDigit;
  23
  24
  25 class AliMUONClusterFinderVS : public TObject
  26 {
  27  public:
  28     AliMUONClusterFinderVS();
  29     virtual ~AliMUONClusterFinderVS();
  30 // Decluster ?
  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();
  38 // Find cluster
  39     virtual void  FindCluster(Int_t i, Int_t j, Int_t cath, AliMUONRawCluster &c);
  40 // Decluster
  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(const AliMUONRawCluster& cluster);
  57 //  Set tracks for debugging
  58     virtual void SetTracks(Int_t t1, Int_t t2) {fTrack[0]=t1; fTrack[1]=t2;}
  59 //  debug level
  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;}
  64 // reset raw clusters
  65     void ResetRawClusters();
  66
  67  protected:
  68     AliMUONClusterFinderVS(const AliMUONClusterFinderVS& clusterFinder);
  69 //  Assignment operator
  70     AliMUONClusterFinderVS & operator = (const AliMUONClusterFinderVS& rhs);
  71
  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
  76
  77
  78 // Configuration
  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.
  90
  91 // Local data store
  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
 103
 104 // Current Fit
 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
 113
 114 // Selected track for debugging
 115     Int_t                    fTrack[2];        // Only digits with main contributions from these tracks are
 116                                                // considered
 117     Int_t                    fDebugLevel;      // prinout control
 118     Int_t                    fSegmentationType;// new seg
 119
 120 //  Return pointer to raw clusters
 121     ClassDef(AliMUONClusterFinderVS,2) //Class for clustering and reconstruction of space points
 122 };
 123 #endif
 124
 125
 126
 127
 128
 129
 130
 131
 132
 133
 134
 135
 136
 137
 138