Correct sign for calculated b_y (outside measured region).
[u/mrichter/AliRoot.git] / 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