TPC/AliDigits.h

   1 #ifndef ALIDIGITS_H
   2 #define ALIDIGITS_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 //  Manager class generaol Alice segment digits
  10 //  segment is for example one pad row in TPC //
  11 ////////////////////////////////////////////////
  12
  13 #include   <TArrayI.h>
  14 #include   <TArrayS.h>
  15 #include   "AliSegmentID.h"
  16 class AliH2F;
  17
  18 class AliDigits: public AliSegmentID{
  19 public:
  20   AliDigits();
  21   AliDigits(const AliDigits &digits); //copy constructor
  22   AliDigits &operator = (const AliDigits & digits); //assignment operator
  23   virtual ~AliDigits();
  24   Short_t * GetDigits(){return fElements->GetArray();}   //return row  pointer to the array digits
  25   Short_t GetDigitFast(Int_t row, Int_t column);  //return value at given row and collumn
  26   void  SetDigitFast(Short_t value,Int_t row, Int_t column);  //set value at given row and collumn
  27   Bool_t BoundsOK(const char *where, Int_t row, Int_t col) ;  //Check If Bound Ok
  28   Bool_t OutOfBoundsError(const char *where, Int_t row, Int_t column);
  29   virtual void Allocate(Int_t rows, Int_t columns);  //construct empty buffer fDigits with size rows x columns
  30   virtual Short_t GetDigit(Int_t row, Int_t column);
  31   virtual void ExpandBuffer();  //expand buffer to twodimensional array
  32   virtual void CompresBuffer(Int_t bufferType,Int_t threshold); //compres buffer according buffertype algorithm
  33   virtual Bool_t First(); //adjust  first valid current digit
  34   virtual Bool_t Next();  //addjust next valid current digit
  35   void SetThreshold(Int_t th) {fThreshold = th;} //set threshold
  36   Int_t  GetThreshold() {return fThreshold;}  //return threshold
  37   Int_t GetNRows(){return fNrows;}
  38   Int_t GetNCols(){return fNcols;}
  39   Int_t CurrentRow(){ return fCurrentRow;}  //return current row
  40   Int_t CurrentColumn(){ return fCurrentCol;} //return current column
  41   Int_t CurrentDigit() {return fElements->At(fCurrentIndex);} //return degit for current row and column
  42   void AcceptHisto(AliH2F * his);  //update buffer for - it will content histogram values
  43   AliH2F * GenerHisto();           //generate 2 dimensional histogram with digits
  44   AliH2F *DrawDigits( const char *option=0,Float_t x1=-1, Float_t x2=-1, Float_t y1=-1, Float_t y2=-1); //draw digits
  45
  46   Int_t GetSize();//return total size of object in bytes
  47   Int_t GetDigitSize(); //return total size of pure digits
  48   Int_t GetOverTh(Float_t threshold,Float_t x1=-1, Float_t x2=-1, Float_t y1=-1, Float_t y2=-1); //return number of digits over threshold
  49 protected:
  50   virtual  void Invalidate();
  51   void ExpandBuffer1(); //expand buffer of type to twodimensional array
  52   void CompresBuffer1(); //compres buffer according  algorithm 1
  53   Bool_t First0();  //first for the buffer type 0
  54   Bool_t Next0();  //next for the buffer type 0
  55   Bool_t First1(); //first for the buffer type 1
  56   Bool_t Next1();//next for the buffer type 1
  57   Short_t  GetDigit1(Int_t row, Int_t column); //return digit for given row and column
  58
  59   Int_t     fNrows;   //number of rows in Segment
  60   Int_t     fNcols; //number of collumns in Segment
  61 private:
  62   TArrayS *fElements;  //buffer of 2 bytes integers for digits
  63   TArrayI *fIndex;  //index position of column
  64   Int_t     fBufType; //type of the buffer - define compression algorithm
  65   Int_t     fThreshold; //treshold for zero suppresion
  66   Int_t     fNelems;  //total number of elements
  67   Int_t fCurrentRow;   //!current row  iteration
  68   Int_t fCurrentCol;   //!current column iteration
  69   Int_t fCurrentIndex; //!current index in field
  70
  71   ClassDef(AliDigits,2)
  72 };
  73
  74
  75
  76 inline Bool_t AliDigits::BoundsOK(const char *where, Int_t row, Int_t col)
  77 {
  78   //Check If Bound Ok
  79   if ( (col>=fNcols) || (col<0) ) return OutOfBoundsError(where,row,col);
  80   Int_t index =(*fIndex).At(col)+row;
  81   if ( (index<0) || (index>fNelems)) return OutOfBoundsError(where,row,col);
  82   return kTRUE;
  83 }
  84
  85 inline Short_t AliDigits::GetDigitFast(Int_t row, Int_t column)
  86 {
  87   //
  88   //return digit from  fDigits array
  89   //if out of range return dummy value  ( value at row = 0, collumn = 0)
  90   //
  91   return fElements->At(fIndex->At(column)+row);
  92 }
  93
  94 inline void  AliDigits::SetDigitFast(Short_t value, Int_t row, Int_t column)
  95 {
  96   //
  97   //set  digit
  98   //
  99   if ( (row<0) || (row>=fNrows)  || (column<0) || (column>=fNcols) )
 100        Error("AliDigits::SetDigitFast", "row %d  col %d out of bounds (size: %d x %d, this: 0x%08x)",
 101            row, column, fNrows, fNcols, this);
 102   (*fElements)[fIndex->At(column)+row]=value;
 103 }
 104
 105 #endif
 106