Put back the two bellows in front of the absorber.
[u/mrichter/AliRoot.git] / TPC / AliTPCclusterer.h
CommitLineData
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
17class TFile;
18class AliTPCParam;
19class AliTPCcluster;
20
21class AliTPCclusterer {
22public:
23 static void Digits2Clusters(const AliTPCParam *par, TFile *of);
24
25private:
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
37private:
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
46inline 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