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, Int_t eventn=1);
  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