]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TRD/AliTRDclusterizer.h
AliTPCcalibDB.cxx.difF Add possibility to get RCU configuration, RCU Trigger...
[u/mrichter/AliRoot.git] / TRD / AliTRDclusterizer.h
index 0db02de592b30240e25094bfa9fb0b9049f10d52..624cbd5847cd01ce0b23faa9d295602fc82e03f5 100644 (file)
@@ -31,6 +31,7 @@ class AliTRDtransform;
 class AliTRDCalROC;
 class AliTRDReconstructor;
 class AliTRDCalSingleChamberStatus;
+class AliTRDrawStreamBase;
 
 class AliTRDclusterizer : public TNamed 
 {
@@ -40,26 +41,26 @@ class AliTRDclusterizer : public TNamed
   // steering flags
   enum{
     kTrOwner= BIT(13)  //  toggle online tracklets ownership
-   ,kClOwner= BIT(14)  //  toggle cluster ownership
-   ,kLabels = BIT(15)  //  toggle MC labels for clusters
-   ,kHLT    = BIT(16)  //  HLT mode
-   ,kLUT    = BIT(17)  //  using look up table for cluster's r-phi position
-   ,kGAUS   = BIT(18)  //  using gauss approx. for cluster's r-phi position
+    ,kClOwner= BIT(14)  //  toggle cluster ownership
+    ,kLabels = BIT(15)  //  toggle MC labels for clusters
+    ,kHLT    = BIT(16)  //  HLT mode
+    ,kLUT    = BIT(17)  //  using look up table for cluster's r-phi position
+    ,kGAUS   = BIT(18)  //  using gauss approx. for cluster's r-phi position
+    ,knewDM  = BIT(19)  //  was the digitsmanger created by raw2clusters?
   };
 
   struct MaxStruct
   {
-    Int_t       Row;
-    Int_t       Col;
-    Int_t       Time;
-    UChar_t     padStatus;
-    Short_t     Signals[3];
-    Bool_t      FivePad;
-  MaxStruct():Row(-1),Col(0),Time(0),padStatus(0),FivePad(kFALSE)
-      {memset(Signals, 0, 3*sizeof(Short_t));}
+    Int_t       row;           // row of the current cluster candidate
+    Int_t       col;           // col of the current cluster candidate
+    Int_t       time;          // time -"-
+    Short_t     signals[3];    // signals of the maximum pad and it's twon neigbours
+    UChar_t     padStatus;     // padStatus of the current cluster candidate
+    Bool_t      fivePad;       // is this cluster candidate part of a 5 pad cluster (two overlaping clusters)?
+  MaxStruct():row(-1),col(0),time(0),padStatus(0),fivePad(kFALSE)
+      {memset(signals, 0, 3*sizeof(Short_t));}
     MaxStruct &operator=(const MaxStruct &a)
-    {Row=a.Row; Col=a.Col; Time=a.Time; padStatus=a.padStatus; FivePad=a.FivePad;
-     memcpy(Signals, a.Signals, 3*sizeof(Short_t)); return *this;}
+    {memcpy(this, &a, sizeof(MaxStruct)); return *this;}
   };
   
   AliTRDclusterizer(const AliTRDReconstructor *const rec = 0x0);
@@ -73,7 +74,7 @@ class AliTRDclusterizer : public TNamed
   Bool_t   Open(const Char_t *name, Int_t nEvent = 0);
   Bool_t   OpenInput(Int_t nEvent = 0);
   Bool_t   OpenOutput();
-  Bool_t   OpenOutput(TTree *clusterTree);
+  Bool_t   OpenOutput(TTree *const clusterTree);
   Bool_t   OpenTrackletOutput();
 
   Bool_t   ReadDigits();
@@ -83,7 +84,7 @@ class AliTRDclusterizer : public TNamed
   Bool_t   WriteClusters(Int_t det);
   void     ResetRecPoints();
   virtual TClonesArray    *RecPoints();
-  TClonesArray    *TrackletsArray();
+  virtual TClonesArray    *TrackletsArray();
   Bool_t   WriteTracklets(Int_t det);
 
   Bool_t   Raw2Clusters(AliRawReader *rawReader);
@@ -97,7 +98,7 @@ class AliTRDclusterizer : public TNamed
   void     SetRawVersion(const Int_t iver) { fRawVersion = iver; } // set the expected raw data version
   void             SetReconstructor(const AliTRDReconstructor *rec) {fReconstructor = rec;}
   static UChar_t   GetStatus(Short_t &signal);
-  Int_t            GetAddedClusters() {return fNoOfClusters;}
+  Int_t            GetAddedClusters() const {return fNoOfClusters;}
 
   Bool_t   IsClustersOwner() const {return TestBit(kClOwner);}
   virtual void     SetClustersOwner(Bool_t own=kTRUE) {SetBit(kClOwner, own); if(!own) {fRecPoints = 0x0; fNoOfClusters=0;} }
@@ -109,19 +110,23 @@ protected:
                             ,const Int_t nTimeTotal, const Int_t nexp);
   void             TailCancelation();
 
-  Float_t  Unfold(Double_t eps, Int_t layer, Double_t *padSignal) const;
+  Float_t  Unfold(Double_t eps, Int_t layer, const Double_t *const padSignal) const;
   
-  void             SetPadStatus(const UChar_t status, UChar_t &encoding);
+  void             SetPadStatus(const UChar_t status, UChar_t &encoding) const;
   UChar_t          GetPadStatus(UChar_t encoding) const;
   Int_t            GetCorruption(UChar_t encoding) const;
 
   Bool_t           IsMaximum(const MaxStruct &Max, UChar_t &padStatus, Short_t *const Signals);       //for const correctness reasons not const parameters are given separately
   Bool_t           FivePadCluster(MaxStruct &ThisMax, MaxStruct &NeighbourMax);
   void             CreateCluster(const MaxStruct &Max); 
+
+  virtual void     AddClusterToArray(AliTRDcluster* cluster);
+  virtual void     AddTrackletsToArray();
+
+private:
   inline void      CalcAdditionalInfo(const MaxStruct &Max, Short_t *const signals, Int_t &nPadCount);
-  virtual void     AddClusterToArray(AliTRDcluster *cluster);
-  inline void      AddTrackletsToArray();
 
+protected:
   const AliTRDReconstructor *fReconstructor; //! reconstructor
   AliRunLoader        *fRunLoader;           //! Run Loader
   TTree               *fClusterTree;         //! Tree with the cluster
@@ -140,7 +145,6 @@ protected:
 
   AliTRDarrayADC      *fDigits;               // Array holding the digits
   AliTRDSignalIndex   *fIndexes;              // Array holding the indexes to the digits
-  Float_t              fADCthresh;            // ADC thresholds: There is no ADC threshold anymore, and simParam should not be used in clusterizer. KO
   Float_t              fMaxThresh;            // Threshold value for the maximum
   Float_t              fSigThresh;            // Threshold value for the digit signal
   Float_t              fMinMaxCutSigma;       // Threshold value for the maximum (cut noise)
@@ -158,6 +162,8 @@ protected:
   Int_t                fClusterROC;           // The index to the first cluster of a given ROC
   Int_t                firstClusterROC;       // The number of cluster in a given ROC
   Int_t                fNoOfClusters;         // Number of Clusters already processed and still owned by the clusterizer
+  Int_t                fBaseline;             // Baseline of the ADC values
+  AliTRDrawStreamBase *fRawStream;            // Currently used RawStream
 
   ClassDef(AliTRDclusterizer,6)               //  TRD clusterfinder