Little changes to make g++ version 3.2 compile the hough library.
[u/mrichter/AliRoot.git] / HLT / hough / AliL3Histogram1D.cxx
1 //$Id$
2
3 // Author: Anders Vestbo <mailto:vestbo@fi.uib.no>
4 //*-- Copyright &copy ASV
5
6 #include "AliL3StandardIncludes.h"
7
8 #include "AliL3Logging.h"
9 #include "AliL3Histogram1D.h"
10
11 #if GCCVERSION == 3
12 using namespace std;
13 #endif
14
15 //_____________________________________________________________
16 // AliL3Histogram1D
17 //
18 // 1D histogram class.
19
20 ClassImp(AliL3Histogram1D)
21
22 AliL3Histogram1D::AliL3Histogram1D()
23 {
24   fNbins = 0;
25   fNcells = 0;
26   fEntries = 0;
27   fXmin = 0;
28   fXmax = 0;
29 #ifdef use_root
30   fRootHisto = 0;
31 #endif
32   fThreshold = 0;
33   fContent = 0;
34   
35 }
36   
37 AliL3Histogram1D::AliL3Histogram1D(Char_t *name,Char_t *id,Int_t nxbin,Double_t xmin,Double_t xmax)
38
39 {
40   
41   strcpy(fName,name);
42   fNbins = nxbin;
43   fNcells = fNbins + 2;
44   fEntries = 0;
45   fXmin = xmin;
46   fXmax = xmax;
47 #ifdef use_root
48   fRootHisto = 0;
49 #endif
50   fThreshold = 0;
51   
52   fContent = new Double_t[fNcells];
53   Reset();
54 }
55
56 AliL3Histogram1D::~AliL3Histogram1D()
57 {
58   //Destructor
59   if(fContent)
60     delete [] fContent;
61 #ifdef use_root
62   if(fRootHisto)
63     delete fRootHisto;
64 #endif
65 }
66
67
68 void AliL3Histogram1D::Reset()
69 {
70   bzero(fContent,fNcells*sizeof(Double_t));
71   fEntries=0;
72 }
73
74 void AliL3Histogram1D::Fill(Double_t x,Int_t weight)
75 {
76   Int_t bin = FindBin(x);
77   AddBinContent(bin,weight);
78 }
79
80
81 Int_t AliL3Histogram1D::FindBin(Double_t x)
82 {
83   if(x < fXmin || x > fXmax)
84     return 0;
85   
86   return 1 + (Int_t)(fNbins*(x-fXmin)/(fXmax-fXmin));
87
88 }
89
90 Int_t AliL3Histogram1D::GetMaximumBin()
91 {
92   Double_t max_value=0;
93   Int_t max_bin=0;
94   for(Int_t i=0; i<fNcells; i++)
95     {
96       if(fContent[i] > max_value)
97         {
98           max_value=fContent[i];
99           max_bin = i;
100         }
101     }
102   return max_bin;
103 }
104
105 Double_t AliL3Histogram1D::GetBinContent(Int_t bin)
106 {
107   if(bin >= fNcells)
108     {
109       LOG(AliL3Log::kError,"AliL3Histogram::GetBinContent","array")<<AliL3Log::kDec<<
110         "bin out of range "<<bin<<ENDLOG;
111       return 0;
112     }
113   
114   if(fContent[bin] < fThreshold)
115     return 0;
116   return fContent[bin];
117 }
118
119
120 void AliL3Histogram1D::SetBinContent(Int_t bin,Int_t value)
121 {
122
123   if(bin >= fNcells)
124     {
125       LOG(AliL3Log::kError,"AliL3Histogram::SetBinContent","array")<<AliL3Log::kDec<<
126         "bin out of range "<<bin<<ENDLOG;
127       return;
128     }
129   if(bin == 0)
130     return;
131   fContent[bin]=value;
132   
133 }
134
135 void AliL3Histogram1D::AddBinContent(Int_t bin,Int_t weight)
136 {
137   if(bin < 0 || bin > fNcells)
138     {
139       LOG(AliL3Log::kError,"AliL3Histogram::AddBinContent","array")<<AliL3Log::kDec<<
140         "bin-value out of range "<<bin<<ENDLOG;
141       return;
142     }
143   if(bin == 0)
144     return;
145   fEntries++;
146   fContent[bin] += weight;
147 }
148
149 Double_t AliL3Histogram1D::GetBinCenter(Int_t bin)
150 {
151   
152   Double_t binwidth = (fXmax - fXmin) / fNbins;
153   return fXmin + (bin-1) * binwidth + 0.5*binwidth;
154   
155 }
156
157 #ifdef use_root
158 void AliL3Histogram1D::Draw(Char_t *option)
159 {
160   fRootHisto = new TH1F(fName,"",fNbins,fXmin,fXmax);
161   for(Int_t bin=0; bin<fNcells; bin++)
162     fRootHisto->AddBinContent(bin,GetBinContent(bin));
163   
164   fRootHisto->Draw(option);
165   
166 }
167 #endif