]>
Commit | Line | Data |
---|---|---|
73042f01 | 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 | //------------------------------------------------------- | |
c630aafd | 9 | // The TPC cluster finder |
73042f01 | 10 | // |
11 | // Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch | |
12 | //------------------------------------------------------- | |
c630aafd | 13 | #include <TObject.h> |
73042f01 | 14 | |
c630aafd | 15 | class TTree; |
73042f01 | 16 | class AliTPCParam; |
17 | class AliTPCcluster; | |
18 | ||
c630aafd | 19 | class AliTPCclusterer : public TObject { |
73042f01 | 20 | public: |
c630aafd | 21 | AliTPCclusterer(const AliTPCParam *par); |
22 | Int_t Digits2Clusters(TTree *dig, TTree *clu); | |
73042f01 | 23 | |
24 | private: | |
c630aafd | 25 | class AliBin { |
26 | public: | |
27 | UShort_t GetQ() const {return fQ;} | |
28 | UInt_t GetMask() const {return fMask;} | |
29 | void SetQ(UShort_t q) {fQ=q;} | |
30 | void SetMask(UInt_t m) {fMask=m;} | |
31 | private: | |
32 | UShort_t fQ; //signal | |
33 | UInt_t fMask; //peak mask | |
34 | }; | |
73042f01 | 35 | |
36 | private: | |
c630aafd | 37 | static Bool_t IsMaximum(Int_t k, Int_t max, const AliBin *bins); |
73042f01 | 38 | static void FindPeaks(Int_t k,Int_t m,AliBin*b,Int_t*idx,UInt_t*msk,Int_t&n); |
c630aafd | 39 | static void MarkPeak(Int_t k, Int_t max, AliBin *bins, UInt_t m); |
40 | static void MakeCluster(Int_t k,Int_t max,AliBin *bins,UInt_t m, | |
73042f01 | 41 | AliTPCcluster &c); |
c630aafd | 42 | |
d7cbafb5 | 43 | const AliTPCParam *fPar; //! pointer to the TPC parameters |
c630aafd | 44 | |
45 | ClassDef(AliTPCclusterer,1) // the TPC cluster finder | |
73042f01 | 46 | }; |
47 | ||
48 | ||
49 | inline Bool_t AliTPCclusterer::IsMaximum(Int_t k,Int_t max,const AliBin *bins){ | |
50 | //is this a local maximum ? | |
51 | UShort_t q=bins[k].GetQ(); | |
52 | if (q==1023) return kFALSE; | |
53 | if (bins[k-max].GetQ() > q) return kFALSE; | |
54 | if (bins[k-1 ].GetQ() > q) return kFALSE; | |
55 | if (bins[k+max].GetQ() > q) return kFALSE; | |
56 | if (bins[k+1 ].GetQ() > q) return kFALSE; | |
57 | if (bins[k-max-1].GetQ() > q) return kFALSE; | |
58 | if (bins[k+max-1].GetQ() > q) return kFALSE; | |
59 | if (bins[k+max+1].GetQ() > q) return kFALSE; | |
60 | if (bins[k-max+1].GetQ() > q) return kFALSE; | |
61 | return kTRUE; | |
62 | } | |
63 | ||
64 | //----------------------------------------------------------------- | |
65 | ||
66 | #endif | |
67 | ||
68 |