TPC/AliTPCclustererMI.h

   1 #ifndef ALITPCCLUSTERERMI_H
   2 #define ALITPCCLUSTERERMI_H
   3
   4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   5  * See cxx source for full Copyright notice                               */
   6
   7 /* $Id$ */
   8
   9
  10
  11 //-------------------------------------------------------
  12 //                       TPC clusterer
  13 //
  14 //   Origin: Marian Ivanov
  15 //-------------------------------------------------------
  16 #include <Rtypes.h>
  17 #include <TObject.h>
  18 #define kMAXCLUSTER 2500
  19
  20 class TFile;
  21 class AliTPCParam;
  22 class AliTPCRecoParam;
  23 class AliTPCclusterMI;
  24 class AliTPCClustersRow;
  25 class AliRawReader;
  26 class AliSimDigits;
  27 class TTree;
  28 class TTreeSRedirector;
  29 class  AliRawEventHeaderBase;
  30 class AliTPCCalROC;
  31 class TVirtualFFT;
  32
  33 class AliTPCclustererMI : public TObject{
  34 public:
  35   AliTPCclustererMI(const AliTPCParam* par, const AliTPCRecoParam * recoParam = 0);
  36   AliTPCclustererMI(const AliTPCclustererMI &param); // copy constructor
  37   AliTPCclustererMI &operator = (const AliTPCclustererMI & param); //assignment
  38   virtual ~AliTPCclustererMI();
  39   virtual void Digits2Clusters();
  40   virtual void Digits2Clusters(AliRawReader* rawReader);
  41   virtual void SetOldRCUFormat(Bool_t rcuFormat = kFALSE)
  42     { fIsOldRCUFormat = rcuFormat; };
  43   virtual void SetInput(TTree * tree);  // set input tree with digits
  44   virtual void SetOutput(TTree * tree); //set output tree with
  45 private:
  46   Bool_t IsMaximum(Float_t k, Int_t max, const Float_t *bins) const;
  47   void MakeCluster2(Int_t k,Int_t max,Float_t *bins,UInt_t m,
  48    AliTPCclusterMI &c);
  49   void MakeCluster(Int_t k,Int_t max,Float_t *bins,UInt_t m,
  50    AliTPCclusterMI &c);
  51   Float_t  GetSigmaY2(Int_t iz);
  52   Float_t  GetSigmaZ2(Int_t iz);
  53   Float_t  FitMax(Float_t vmatrix[5][5], Float_t y, Float_t z, Float_t sigmay, Float_t sigmaz);
  54   void AddCluster(AliTPCclusterMI &c, Float_t *matrix, Int_t pos);  // add the cluster to the array
  55   void UnfoldCluster(Float_t * matrix[7], Float_t recmatrix[5][5],
  56                      Float_t & meani, Float_t & meanj, Float_t & sum, Float_t &overlap );
  57   void FindClusters(AliTPCCalROC * noiseROC);
  58   Double_t  ProcesSignal(Float_t * signal, Int_t nchannels, Int_t id[3], Double_t &rms, Double_t &pedestalCalib);
  59   void DumpHistos();
  60   Int_t  TransformFFT(Float_t *input, Float_t threshold, Bool_t locMax, Float_t *freq, Float_t *re, Float_t *im, Float_t *mag, Float_t *phi);
  61
  62   Float_t * fBins;       //!digits array
  63   Float_t * fResBins;    //!digits array with res. after 1 finder
  64   Int_t fLoop;         //loop - cf in 2 loops
  65   Int_t fMaxBin;       //current ( for current sector)  maximal bin
  66   Int_t fMaxTime;      //current ( for current sector)  maximal time
  67   Int_t fMaxPad;       //current ( for current sector)  maximal pad
  68   Int_t fSector;      //!current sector
  69   Int_t fRow;         //!current row
  70   Float_t fSign;      //!current sign
  71   Float_t fRx;        // current radius
  72   Float_t fPadWidth;  // the width of the pad
  73   Float_t fPadLength;  // the width of the pad
  74   Float_t fZWidth;     //the z bin width
  75
  76   Bool_t  fPedSubtraction; // perform pedestal subtraction or not
  77   Bool_t  fIsOldRCUFormat; // assume old RCU raw data format
  78   AliRawEventHeaderBase *fEventHeader; //! event header information
  79   UInt_t  fTimeStamp;   // Time Stamp
  80   UInt_t  fEventType;   // Event Type
  81   TTree * fInput;   //!input  tree with digits - object not owner
  82   TTree * fOutput;   //!output tree with digits - object not owner
  83   AliTPCClustersRow * fRowCl;  //! current cluster row
  84   AliSimDigits * fRowDig;      //! current digits row
  85   const AliTPCParam * fParam;        //! tpc parameters
  86   Int_t fNcluster;             // number of clusters - for given row
  87   TObjArray * fAmplitudeHisto;          //! array of histograms of amplitudes
  88   TTreeSRedirector *fDebugStreamer;     //!debug streamer
  89   const AliTPCRecoParam  * fRecoParam;        //! reconstruction parameters
  90   TVirtualFFT *fFFTr2c;                 //! Fast Furier transform object
  91   ClassDef(AliTPCclustererMI,1)  // Time Projection Chamber digits
  92 };
  93
  94 inline Bool_t AliTPCclustererMI::IsMaximum(Float_t q,Int_t max,const Float_t *bins) const {
  95   //is this a local maximum ?
  96   if (bins[-max] >= q) return kFALSE;
  97   if (bins[-1  ] >= q) return kFALSE;
  98   if (bins[+max] > q) return kFALSE;
  99   if (bins[+1  ] > q) return kFALSE;
 100   if (bins[-max-1] >= q) return kFALSE;
 101   if (bins[+max-1] >= q) return kFALSE;
 102   if (bins[+max+1] > q) return kFALSE;
 103   if (bins[-max+1] >= q) return kFALSE;
 104   return kTRUE;
 105 }
 106
 107
 108
 109 //-----------------------------------------------------------------
 110
 111 #endif
 112
 113