]>
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 | //------------------------------------------------------- | |
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: | |
afc42102 | 23 | static void Digits2Clusters(const AliTPCParam *par, TFile *of, Int_t eventn=1); |
73042f01 | 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 |