]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TPC/AliTPCclustererMI.h
undo accidental commit of CMake and LinkDef file without PhiCorrelation Task
[u/mrichter/AliRoot.git] / TPC / AliTPCclustererMI.h
index 588085bd65af56169eda982f4bcc99cfd5567ccf..01b1c3804401097c6ad58bae0bad9efc2dea0ae6 100644 (file)
@@ -15,6 +15,7 @@
 //-------------------------------------------------------
 #include <Rtypes.h>
 #include <TObject.h>
+#include <AliTPCRecoParam.h>
 #define kMAXCLUSTER 2500
 
 class TFile;
@@ -37,8 +38,9 @@ public:
   AliTPCclustererMI &operator = (const AliTPCclustererMI & param); //assignment
   virtual ~AliTPCclustererMI();
   virtual void Digits2Clusters();
+  virtual void Digits2ClustersOld(AliRawReader* rawReader);
   virtual void Digits2Clusters(AliRawReader* rawReader);
-  virtual void SetInput(TTree * tree);  // set input tree with digits    
+  virtual void SetInput(TTree * tree);  // set input tree with digits
   virtual void SetOutput(TTree * tree); // set output tree with 
   virtual void FillRow();               // fill the output container - Tree or TObjArray
   TObjArray * GetOutputArray(){return fOutputArray;}
@@ -61,7 +63,8 @@ private:
   void FindClusters(AliTPCCalROC * noiseROC);
   Bool_t AcceptCluster(AliTPCclusterMI*c);
   Double_t  ProcesSignal(Float_t * signal, Int_t nchannels, Int_t id[3], Double_t &rms, Double_t &pedestalCalib);
-
+  void ProcessSectorData();
+  
   Float_t * fBins;       //!digits array
   Int_t   * fSigBins; //!digits array containg only timebins above threshold
   Int_t     fNSigBins;//!size of fSigBins
@@ -94,6 +97,13 @@ private:
   Bool_t  fBDumpSignal; // dump signal flag
   Bool_t  fBClonesArray; // output clusters stored in TClonesArray 
 
+  // Non-persistent arrays
+
+  Float_t** fAllBins; //! All sector bins
+  Int_t** fAllSigBins;//! All signal bins in a sector
+  Int_t*  fAllNSigBins;//! Number of signal bins in a sector
+
+
   ClassDef(AliTPCclustererMI,2)  // Time Projection Chamber digits
 };
 
@@ -107,6 +117,17 @@ inline Bool_t AliTPCclustererMI::IsMaximum(Float_t q,Int_t max,const Float_t *bi
   if (bins[+max-1] >= q) return kFALSE; 
   if (bins[+max+1] > q) return kFALSE; 
   if (bins[-max+1] >= q) return kFALSE;
+  //
+  //
+  if (fRecoParam->GetClusterMaxRange(1)>0){  //local maxima in z on more than 1 time bin 
+    if (bins[-2] > q) return kFALSE;
+    if (bins[ 2] > q) return kFALSE;
+  }
+  if (fRecoParam->GetClusterMaxRange(0)>0){  //local maxima in y on more than pad 
+    if (bins[-2*max] > q) return kFALSE;
+    if (bins[ 2*max] > q) return kFALSE;
+  }
+
   return kTRUE; 
 }