]>
Commit | Line | Data |
---|---|---|
95a00d93 | 1 | //$Id$ |
2 | ||
b1886074 | 3 | // Author: Anders Vestbo <mailto:vestbo@fi.uib.no> |
4 | //*-- Copyright © ASV | |
f7473b79 | 5 | |
e06900d5 | 6 | #include "AliL3StandardIncludes.h" |
7 | ||
f7473b79 | 8 | #include "AliL3Logging.h" |
9 | #include "AliL3Histogram1D.h" | |
10 | ||
e06900d5 | 11 | #if GCCVERSION == 3 |
12 | using namespace std; | |
13 | #endif | |
14 | ||
b1886074 | 15 | //_____________________________________________________________ |
16 | // AliL3Histogram1D | |
17 | // | |
18 | // 1D histogram class. | |
f7473b79 | 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 | } | |
f7473b79 | 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; | |
95a00d93 | 47 | #ifdef use_root |
f7473b79 | 48 | fRootHisto = 0; |
95a00d93 | 49 | #endif |
f7473b79 | 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 | ||
b1886074 | 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 | ||
f7473b79 | 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 |