X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=TPC%2FAliTPCclustererMI.h;h=d854dc0e311f6aa4bbec772ed6c2e42d6723fe3b;hb=61098ddafdd4a3351a49259599094a55986950d0;hp=43c750dd84edb37b42a6caabb0d9a3e3446ac909;hpb=1c53abe24c5f70da7045524e204946b440fe5b21;p=u%2Fmrichter%2FAliRoot.git

diff --git a/TPC/AliTPCclustererMI.h b/TPC/AliTPCclustererMI.h
index 43c750dd84e..d854dc0e311 100644
--- a/TPC/AliTPCclustererMI.h
+++ b/TPC/AliTPCclustererMI.h
@@ -1,6 +1,13 @@
 #ifndef ALITPCCLUSTERERMI_H
 #define ALITPCCLUSTERERMI_H
 
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/* $Id$ */
+
+
+
 //-------------------------------------------------------
 //                       TPC clusterer
 //
@@ -12,55 +19,75 @@
 
 class TFile;
 class AliTPCParam;
+class AliTPCRecoParam;
 class AliTPCclusterMI;
 class AliTPCClustersRow;
+class AliRawReader;
 class AliSimDigits;
 class TTree;
+class TTreeSRedirector;
+class  AliRawEventHeaderBase;
+class AliTPCCalROC;
 
 class AliTPCclustererMI : public TObject{
 public:
-  AliTPCclustererMI();
-  virtual void Digits2Clusters(const AliTPCParam *par, Int_t eventn=1);
+  AliTPCclustererMI(const AliTPCParam* par, const AliTPCRecoParam * recoParam = 0);
+  AliTPCclustererMI(const AliTPCclustererMI &param); // copy constructor
+  AliTPCclustererMI &operator = (const AliTPCclustererMI & param); //assignment
+  virtual ~AliTPCclustererMI();
+  virtual void Digits2Clusters();
+  virtual void Digits2Clusters(AliRawReader* rawReader);
   virtual void SetInput(TTree * tree);  // set input tree with digits    
-  virtual void SetOutput(TTree * tree); //set output tree with 
+  virtual void SetOutput(TTree * tree); // set output tree with 
+  virtual void FillRow();               // fill the output container - Tree or TObjArray
+  TObjArray * GetOutputArray(){return fOutputArray;}
 private:
-  Bool_t IsMaximum(Int_t k, Int_t max, const Int_t *bins); 
-  void MakeCluster2(Int_t k,Int_t max,Int_t *bins,UInt_t m,
+  Bool_t IsMaximum(Float_t k, Int_t max, const Float_t *bins) const; 
+  void MakeCluster2(Int_t k,Int_t max,Float_t *bins,UInt_t m,
    AliTPCclusterMI &c);  
-  void MakeCluster(Int_t k,Int_t max,Int_t *bins,UInt_t m,
+  void MakeCluster(Int_t k,Int_t max,Float_t *bins,UInt_t m,
    AliTPCclusterMI &c); 
   Float_t  GetSigmaY2(Int_t iz);
   Float_t  GetSigmaZ2(Int_t iz);
   Float_t  FitMax(Float_t vmatrix[5][5], Float_t y, Float_t z, Float_t sigmay, Float_t sigmaz);
-  void AddCluster(AliTPCclusterMI &c);  // add the cluster to the array
-  void UnfoldCluster(Int_t * matrix[7], Float_t recmatrix[5][5], 
+  void AddCluster(AliTPCclusterMI &c, Float_t *matrix, Int_t pos);  // add the cluster to the array
+  void UnfoldCluster(Float_t * matrix[7], Float_t recmatrix[5][5], 
 		     Float_t & meani, Float_t & meanj, Float_t & sum, Float_t &overlap );
+  void FindClusters(AliTPCCalROC * noiseROC);
+  Double_t  ProcesSignal(Float_t * signal, Int_t nchannels, Int_t id[3], Double_t &rms, Double_t &pedestalCalib);
 
-
-
-  Int_t * fBins;       //!digits array
-  Int_t * fResBins;    //!digits array with res. after 1 finder
+  Float_t * fBins;       //!digits array
+  Int_t   * fSigBins; //!digits array containg only timebins above threshold
+  Int_t     fNSigBins;//!size of fSigBins
   Int_t fLoop;         //loop - cf in 2 loops
-  Int_t fMaxBin;
-  Int_t fMaxTime;
-  Int_t fMaxPad;
+  Int_t fMaxBin;       //current ( for current sector)  maximal bin
+  Int_t fMaxTime;      //current ( for current sector)  maximal time
+  Int_t fMaxPad;       //current ( for current sector)  maximal pad
   Int_t fSector;      //!current sector
+  Int_t fRow;         //!current row
   Float_t fSign;      //!current sign 
   Float_t fRx;        // current radius
   Float_t fPadWidth;  // the width of the pad
   Float_t fPadLength;  // the width of the pad
   Float_t fZWidth;     //the z bin width
-
+  Bool_t  fPedSubtraction; // perform pedestal subtraction or not
+  AliRawEventHeaderBase *fEventHeader; //! event header information
+  UInt_t  fTimeStamp;   // Time Stamp
+  UInt_t  fEventType;   // Event Type
   TTree * fInput;   //!input  tree with digits - object not owner
   TTree * fOutput;   //!output tree with digits - object not owner
+  TObjArray *fOutputArray;     //! output TObjArray with pointers arrays of cluster
   AliTPCClustersRow * fRowCl;  //! current cluster row
   AliSimDigits * fRowDig;      //! current digits row
   const AliTPCParam * fParam;        //! tpc parameters
   Int_t fNcluster;             // number of clusters - for given row
-  ClassDef(AliTPCclustererMI,1)  // Time Projection Chamber digits
+  TTreeSRedirector *fDebugStreamer;     //!debug streamer
+  const AliTPCRecoParam  * fRecoParam;        //! reconstruction parameters
+  Bool_t  fBDumpSignal; // dump signal flag
+  ClassDef(AliTPCclustererMI,2)  // Time Projection Chamber digits
 };
 
-inline Bool_t AliTPCclustererMI::IsMaximum(Int_t q,Int_t max,const Int_t *bins){
+inline Bool_t AliTPCclustererMI::IsMaximum(Float_t q,Int_t max,const Float_t *bins) const {
   //is this a local maximum ?
   if (bins[-max] >= q) return kFALSE;
   if (bins[-1  ] >= q) return kFALSE;