a85475b06aa3725cffe4a75e68a8dafbf674fce4
[u/mrichter/AliRoot.git] / TPC / AliTPCclusterer.h
1 #ifndef ALITPCCLUSTERER_H
2 #define ALITPCCLUSTERER_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
8 //-------------------------------------------------------
9 //                       TPC clusterer
10 //
11 //   Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch 
12 //-------------------------------------------------------
13 #include <Rtypes.h>
14
15 #define kMAXCLUSTER 2500
16
17 class TFile;
18 class AliTPCParam;
19 class AliTPCcluster;
20
21 class AliTPCclusterer {
22 public:
23    static void Digits2Clusters(const AliTPCParam *par, TFile *of);
24
25 private:
26    class AliBin {
27    public:
28      UShort_t GetQ()    const   {return fQ;}
29      UInt_t   GetMask() const   {return fMask;}
30      void     SetQ(UShort_t q)  {fQ=q;}
31      void     SetMask(UInt_t m) {fMask=m;}
32    private:
33      UShort_t fQ;  //signal
34      UInt_t fMask; //peak mask
35    };
36
37 private:
38    static Bool_t IsMaximum(Int_t k, Int_t max, const AliBin *bins); 
39   static void FindPeaks(Int_t k,Int_t m,AliBin*b,Int_t*idx,UInt_t*msk,Int_t&n);
40    static void MarkPeak(Int_t k, Int_t max, AliBin *bins, UInt_t m);
41    static void MakeCluster(Int_t k,Int_t max,AliBin *bins,UInt_t m,
42    AliTPCcluster &c);
43 };
44
45
46 inline Bool_t AliTPCclusterer::IsMaximum(Int_t k,Int_t max,const AliBin *bins){
47   //is this a local maximum ?
48   UShort_t q=bins[k].GetQ();
49   if (q==1023) return kFALSE;
50   if (bins[k-max].GetQ() > q) return kFALSE;
51   if (bins[k-1  ].GetQ() > q) return kFALSE; 
52   if (bins[k+max].GetQ() > q) return kFALSE; 
53   if (bins[k+1  ].GetQ() > q) return kFALSE; 
54   if (bins[k-max-1].GetQ() > q) return kFALSE;
55   if (bins[k+max-1].GetQ() > q) return kFALSE; 
56   if (bins[k+max+1].GetQ() > q) return kFALSE; 
57   if (bins[k-max+1].GetQ() > q) return kFALSE;
58   return kTRUE; 
59 }
60
61 //-----------------------------------------------------------------
62
63 #endif
64
65