| Geant4 Cross Reference |
1 // Copyright (C) 2010, Guy Barrand. All rights 1 // Copyright (C) 2010, Guy Barrand. All rights reserved.
2 // See the file tools.license for terms. 2 // See the file tools.license for terms.
3 3
4 #ifndef tools_histo_c2d 4 #ifndef tools_histo_c2d
5 #define tools_histo_c2d 5 #define tools_histo_c2d
6 6
7 #include "base_cloud" 7 #include "base_cloud"
8 8
9 #include "../mnmx" 9 #include "../mnmx"
10 10
11 #include "h2d" 11 #include "h2d"
12 12
13 namespace tools { 13 namespace tools {
14 namespace histo { 14 namespace histo {
15 15
16 class c2d : public base_cloud { 16 class c2d : public base_cloud {
17 public: 17 public:
18 static const std::string& s_class() { 18 static const std::string& s_class() {
19 static const std::string s_v("tools::histo 19 static const std::string s_v("tools::histo::c2d");
20 return s_v; 20 return s_v;
21 } 21 }
22 public: 22 public:
23 bool set_title(const std::string&); 23 bool set_title(const std::string&);
24 unsigned int dimension() const {return 2;} 24 unsigned int dimension() const {return 2;}
25 bool reset(); 25 bool reset();
26 unsigned int entries() const { return m_hist 26 unsigned int entries() const { return m_histo ? m_histo->all_entries() : (unsigned int)m_ws.size();}
27 public: 27 public:
28 double sum_of_weights() const { return (m_hi 28 double sum_of_weights() const { return (m_histo ? m_histo->sum_bin_heights() : m_Sw);}
29 bool convert_to_histogram(); 29 bool convert_to_histogram();
30 bool is_converted() const {return m_histo ? 30 bool is_converted() const {return m_histo ? true : false;}
31 bool scale(double); 31 bool scale(double);
32 public: 32 public:
33 bool fill(double,double,double = 1); 33 bool fill(double,double,double = 1);
34 double lower_edge_x() const {return m_histo 34 double lower_edge_x() const {return m_histo ? m_histo->axis_x().lower_edge() : m_lower_x;}
35 double lower_edge_y() const { return m_histo 35 double lower_edge_y() const { return m_histo ? m_histo->axis_y().lower_edge() : m_lower_y;}
36 double upper_edge_x() const { return m_histo 36 double upper_edge_x() const { return m_histo ? m_histo->axis_x().upper_edge() : m_upper_x;}
37 double upper_edge_y() const { return m_histo 37 double upper_edge_y() const { return m_histo ? m_histo->axis_y().upper_edge() : m_upper_y;}
38 double value_x(unsigned int a_index) const { 38 double value_x(unsigned int a_index) const { return m_histo ? 0 : m_xs[a_index];}
39 double value_y(unsigned int a_index) const { 39 double value_y(unsigned int a_index) const { return m_histo ? 0 : m_ys[a_index];}
40 double weight(unsigned int a_index) const { 40 double weight(unsigned int a_index) const { return m_histo ? 0 : m_ws[a_index];}
41 double mean_x() const {return m_histo ? m_hi 41 double mean_x() const {return m_histo ? m_histo->mean_x() : (m_Sw?m_Sxw/m_Sw:0);}
42 double mean_y() const {return m_histo ? m_hi 42 double mean_y() const {return m_histo ? m_histo->mean_y() : (m_Sw?m_Syw/m_Sw:0);}
43 double rms_x() const; 43 double rms_x() const;
44 double rms_y() const; 44 double rms_y() const;
45 bool convert(unsigned int,double,double,unsi 45 bool convert(unsigned int,double,double,unsigned int,double,double);
46 bool convert(const std::vector<double>&,cons 46 bool convert(const std::vector<double>&,const std::vector<double>&);
47 const histo::h2d& histogram() const; 47 const histo::h2d& histogram() const;
48 template <class HISTO> 48 template <class HISTO>
49 bool fill_histogram(HISTO& a_histo) const { 49 bool fill_histogram(HISTO& a_histo) const {
50 size_t number = m_xs.size(); 50 size_t number = m_xs.size();
51 for(size_t index=0;index<number;index++) { 51 for(size_t index=0;index<number;index++) {
52 if(!a_histo.fill(m_xs[index],m_ys[index] 52 if(!a_histo.fill(m_xs[index],m_ys[index],m_ws[index])) return false;
53 } 53 }
54 return true; 54 return true;
55 } 55 }
56 bool set_conversion_parameters(unsigned int, 56 bool set_conversion_parameters(unsigned int,double,double,unsigned int,double,double);
57 57
58 bool set_histogram(h2d* a_histo){ //we take 58 bool set_histogram(h2d* a_histo){ //we take ownership of a_histo.
59 reset(); 59 reset();
60 m_histo = a_histo; 60 m_histo = a_histo;
61 return true; 61 return true;
62 } 62 }
63 public: 63 public:
64 c2d(); 64 c2d();
65 c2d(const std::string&,int=base_cloud::UNLIM 65 c2d(const std::string&,int=base_cloud::UNLIMITED());
66 virtual ~c2d(){delete m_histo;} 66 virtual ~c2d(){delete m_histo;}
67 public: 67 public:
68 c2d(const c2d& a_from) 68 c2d(const c2d& a_from)
69 :base_cloud(a_from) 69 :base_cloud(a_from)
70 ,m_xs(a_from.m_xs) 70 ,m_xs(a_from.m_xs)
71 ,m_ys(a_from.m_ys) 71 ,m_ys(a_from.m_ys)
72 ,m_lower_x(a_from.m_lower_x) 72 ,m_lower_x(a_from.m_lower_x)
73 ,m_upper_x(a_from.m_upper_x) 73 ,m_upper_x(a_from.m_upper_x)
74 ,m_lower_y(a_from.m_lower_y) 74 ,m_lower_y(a_from.m_lower_y)
75 ,m_upper_y(a_from.m_upper_y) 75 ,m_upper_y(a_from.m_upper_y)
76 ,m_Sxw(a_from.m_Sxw) 76 ,m_Sxw(a_from.m_Sxw)
77 ,m_Sx2w(a_from.m_Sx2w) 77 ,m_Sx2w(a_from.m_Sx2w)
78 ,m_Syw(a_from.m_Syw) 78 ,m_Syw(a_from.m_Syw)
79 ,m_Sy2w(a_from.m_Sy2w) 79 ,m_Sy2w(a_from.m_Sy2w)
80 ,m_cnv_x_num(a_from.m_cnv_x_num) 80 ,m_cnv_x_num(a_from.m_cnv_x_num)
81 ,m_cnv_x_min(a_from.m_cnv_x_min) 81 ,m_cnv_x_min(a_from.m_cnv_x_min)
82 ,m_cnv_x_max(a_from.m_cnv_x_max) 82 ,m_cnv_x_max(a_from.m_cnv_x_max)
83 ,m_cnv_y_num(a_from.m_cnv_y_num) 83 ,m_cnv_y_num(a_from.m_cnv_y_num)
84 ,m_cnv_y_min(a_from.m_cnv_y_min) 84 ,m_cnv_y_min(a_from.m_cnv_y_min)
85 ,m_cnv_y_max(a_from.m_cnv_y_max) 85 ,m_cnv_y_max(a_from.m_cnv_y_max)
86 ,m_histo(0) 86 ,m_histo(0)
87 { 87 {
88 if(a_from.m_histo) { 88 if(a_from.m_histo) {
89 m_histo = new histo::h2d(*a_from.m_histo 89 m_histo = new histo::h2d(*a_from.m_histo);
90 } 90 }
91 } 91 }
92 92
93 c2d& operator=(const c2d& a_from) { 93 c2d& operator=(const c2d& a_from) {
94 base_cloud::operator=(a_from); 94 base_cloud::operator=(a_from);
95 if(&a_from==this) return *this; 95 if(&a_from==this) return *this;
96 m_xs = a_from.m_xs; 96 m_xs = a_from.m_xs;
97 m_ys = a_from.m_ys; 97 m_ys = a_from.m_ys;
98 m_lower_x = a_from.m_lower_x; 98 m_lower_x = a_from.m_lower_x;
99 m_upper_x = a_from.m_upper_x; 99 m_upper_x = a_from.m_upper_x;
100 m_lower_y = a_from.m_lower_y; 100 m_lower_y = a_from.m_lower_y;
101 m_upper_y = a_from.m_upper_y; 101 m_upper_y = a_from.m_upper_y;
102 m_Sxw = a_from.m_Sxw; 102 m_Sxw = a_from.m_Sxw;
103 m_Sx2w = a_from.m_Sx2w; 103 m_Sx2w = a_from.m_Sx2w;
104 m_Syw = a_from.m_Syw; 104 m_Syw = a_from.m_Syw;
105 m_Sy2w = a_from.m_Sy2w; 105 m_Sy2w = a_from.m_Sy2w;
106 m_cnv_x_num = a_from.m_cnv_x_num; 106 m_cnv_x_num = a_from.m_cnv_x_num;
107 m_cnv_x_min = a_from.m_cnv_x_min; 107 m_cnv_x_min = a_from.m_cnv_x_min;
108 m_cnv_x_max = a_from.m_cnv_x_max; 108 m_cnv_x_max = a_from.m_cnv_x_max;
109 m_cnv_y_num = a_from.m_cnv_y_num; 109 m_cnv_y_num = a_from.m_cnv_y_num;
110 m_cnv_y_min = a_from.m_cnv_y_min; 110 m_cnv_y_min = a_from.m_cnv_y_min;
111 m_cnv_y_max = a_from.m_cnv_y_max; 111 m_cnv_y_max = a_from.m_cnv_y_max;
112 delete m_histo; 112 delete m_histo;
113 m_histo = 0; 113 m_histo = 0;
114 if(a_from.m_histo) { 114 if(a_from.m_histo) {
115 m_histo = new histo::h2d(*a_from.m_histo 115 m_histo = new histo::h2d(*a_from.m_histo);
116 } 116 }
117 return *this; 117 return *this;
118 } 118 }
119 public: //AIDA API 119 public: //AIDA API
120 double lowerEdgeX() const {return lower_edge 120 double lowerEdgeX() const {return lower_edge_x();}
121 double lowerEdgeY() const {return lower_edge 121 double lowerEdgeY() const {return lower_edge_y();}
122 double upperEdgeX() const {return upper_edge 122 double upperEdgeX() const {return upper_edge_x();}
123 double upperEdgeY() const {return upper_edge 123 double upperEdgeY() const {return upper_edge_y();}
124 template <class HISTO> 124 template <class HISTO>
125 bool fillHistogram(HISTO& a_histo) const {re 125 bool fillHistogram(HISTO& a_histo) const {return fill_histogram<HISTO>(a_histo);}
126 protected: 126 protected:
127 void clear(); 127 void clear();
128 protected: 128 protected:
129 std::vector<double> m_xs; 129 std::vector<double> m_xs;
130 std::vector<double> m_ys; 130 std::vector<double> m_ys;
131 double m_lower_x; 131 double m_lower_x;
132 double m_upper_x; 132 double m_upper_x;
133 double m_lower_y; 133 double m_lower_y;
134 double m_upper_y; 134 double m_upper_y;
135 double m_Sxw; 135 double m_Sxw;
136 double m_Sx2w; 136 double m_Sx2w;
137 double m_Syw; 137 double m_Syw;
138 double m_Sy2w; 138 double m_Sy2w;
139 // 139 //
140 unsigned int m_cnv_x_num; 140 unsigned int m_cnv_x_num;
141 double m_cnv_x_min; 141 double m_cnv_x_min;
142 double m_cnv_x_max; 142 double m_cnv_x_max;
143 unsigned int m_cnv_y_num; 143 unsigned int m_cnv_y_num;
144 double m_cnv_y_min; 144 double m_cnv_y_min;
145 double m_cnv_y_max; 145 double m_cnv_y_max;
146 histo::h2d* m_histo; 146 histo::h2d* m_histo;
147 }; 147 };
148 148
149 }} 149 }}
150 150
151 namespace tools { 151 namespace tools {
152 namespace histo { 152 namespace histo {
153 153
154 inline 154 inline
155 c2d::c2d() 155 c2d::c2d()
156 :base_cloud(UNLIMITED()) 156 :base_cloud(UNLIMITED())
157 ,m_lower_x(0) 157 ,m_lower_x(0)
158 ,m_upper_x(0) 158 ,m_upper_x(0)
159 ,m_lower_y(0) 159 ,m_lower_y(0)
160 ,m_upper_y(0) 160 ,m_upper_y(0)
161 ,m_Sxw(0) 161 ,m_Sxw(0)
162 ,m_Sx2w(0) 162 ,m_Sx2w(0)
163 ,m_Syw(0) 163 ,m_Syw(0)
164 ,m_Sy2w(0) 164 ,m_Sy2w(0)
165 ,m_cnv_x_num(0) 165 ,m_cnv_x_num(0)
166 ,m_cnv_x_min(0) 166 ,m_cnv_x_min(0)
167 ,m_cnv_x_max(0) 167 ,m_cnv_x_max(0)
168 ,m_cnv_y_num(0) 168 ,m_cnv_y_num(0)
169 ,m_cnv_y_min(0) 169 ,m_cnv_y_min(0)
170 ,m_cnv_y_max(0) 170 ,m_cnv_y_max(0)
171 ,m_histo(0) 171 ,m_histo(0)
172 {} 172 {}
173 173
174 inline 174 inline
175 c2d::c2d(const std::string& a_title,int aLimit 175 c2d::c2d(const std::string& a_title,int aLimit)
176 :base_cloud(aLimit) 176 :base_cloud(aLimit)
177 ,m_lower_x(0) 177 ,m_lower_x(0)
178 ,m_upper_x(0) 178 ,m_upper_x(0)
179 ,m_lower_y(0) 179 ,m_lower_y(0)
180 ,m_upper_y(0) 180 ,m_upper_y(0)
181 ,m_Sxw(0) 181 ,m_Sxw(0)
182 ,m_Sx2w(0) 182 ,m_Sx2w(0)
183 ,m_Syw(0) 183 ,m_Syw(0)
184 ,m_Sy2w(0) 184 ,m_Sy2w(0)
185 ,m_cnv_x_num(0) 185 ,m_cnv_x_num(0)
186 ,m_cnv_x_min(0) 186 ,m_cnv_x_min(0)
187 ,m_cnv_x_max(0) 187 ,m_cnv_x_max(0)
188 ,m_cnv_y_num(0) 188 ,m_cnv_y_num(0)
189 ,m_cnv_y_min(0) 189 ,m_cnv_y_min(0)
190 ,m_cnv_y_max(0) 190 ,m_cnv_y_max(0)
191 ,m_histo(0) 191 ,m_histo(0)
192 { 192 {
193 set_title(a_title); 193 set_title(a_title);
194 } 194 }
195 195
196 inline 196 inline
197 void c2d::clear(){ 197 void c2d::clear(){
198 m_lower_x = 0; 198 m_lower_x = 0;
199 m_upper_x = 0; 199 m_upper_x = 0;
200 m_lower_y = 0; 200 m_lower_y = 0;
201 m_upper_y = 0; 201 m_upper_y = 0;
202 m_Sw = 0; 202 m_Sw = 0;
203 m_Sxw = 0; 203 m_Sxw = 0;
204 m_Sx2w = 0; 204 m_Sx2w = 0;
205 m_Syw = 0; 205 m_Syw = 0;
206 m_Sy2w = 0; 206 m_Sy2w = 0;
207 m_xs.clear(); 207 m_xs.clear();
208 m_ys.clear(); 208 m_ys.clear();
209 m_ws.clear(); 209 m_ws.clear();
210 } 210 }
211 211
212 inline 212 inline
213 bool c2d::convert( 213 bool c2d::convert(
214 unsigned int a_bins_x,double a_lower_edge_x,d 214 unsigned int a_bins_x,double a_lower_edge_x,double a_upper_edge_x
215 ,unsigned int a_bins_y,double a_lower_edge_y,d 215 ,unsigned int a_bins_y,double a_lower_edge_y,double a_upper_edge_y
216 ) { 216 ) {
217 if(m_histo) return true; // Done. 217 if(m_histo) return true; // Done.
218 m_histo = new histo::h2d(base_cloud::title() 218 m_histo = new histo::h2d(base_cloud::title(),
219 a_bins_x,a_lower_ed 219 a_bins_x,a_lower_edge_x,a_upper_edge_x,
220 a_bins_y,a_lower_ed 220 a_bins_y,a_lower_edge_y,a_upper_edge_y);
221 if(!m_histo) return false; 221 if(!m_histo) return false;
222 bool status = fill_histogram(*m_histo); 222 bool status = fill_histogram(*m_histo);
223 clear(); 223 clear();
224 return status; 224 return status;
225 } 225 }
226 226
227 inline 227 inline
228 bool c2d::convert_to_histogram(){ 228 bool c2d::convert_to_histogram(){
229 if( (m_cnv_x_num<=0) || (m_cnv_x_max<=m_cnv_ 229 if( (m_cnv_x_num<=0) || (m_cnv_x_max<=m_cnv_x_min) ||
230 (m_cnv_y_num<=0) || (m_cnv_y_max<=m_cnv_ 230 (m_cnv_y_num<=0) || (m_cnv_y_max<=m_cnv_y_min) ) {
231 double dx = 0.01 * (upper_edge_x() - lower 231 double dx = 0.01 * (upper_edge_x() - lower_edge_x())/BINS();
232 double dy = 0.01 * (upper_edge_y() - lower 232 double dy = 0.01 * (upper_edge_y() - lower_edge_y())/BINS();
233 return convert(BINS(),lower_edge_x(),upper 233 return convert(BINS(),lower_edge_x(),upper_edge_x()+dx,
234 BINS(),lower_edge_y(),upper 234 BINS(),lower_edge_y(),upper_edge_y()+dy);
235 } else { 235 } else {
236 return convert(m_cnv_x_num,m_cnv_x_min,m_c 236 return convert(m_cnv_x_num,m_cnv_x_min,m_cnv_x_max,
237 m_cnv_y_num,m_cnv_y_min,m_c 237 m_cnv_y_num,m_cnv_y_min,m_cnv_y_max);
238 } 238 }
239 } 239 }
240 240
241 241
242 inline 242 inline
243 bool c2d::set_title(const std::string& a_title 243 bool c2d::set_title(const std::string& a_title){
244 m_title = a_title; 244 m_title = a_title;
245 if(m_histo) m_histo->set_title(a_title); 245 if(m_histo) m_histo->set_title(a_title);
246 return true; 246 return true;
247 } 247 }
248 248
249 inline 249 inline
250 bool c2d::scale(double a_scale) { 250 bool c2d::scale(double a_scale) {
251 if(m_histo) { 251 if(m_histo) {
252 return m_histo->scale(a_scale); 252 return m_histo->scale(a_scale);
253 } else { 253 } else {
254 size_t number = m_ws.size(); 254 size_t number = m_ws.size();
255 for(size_t index=0;index<number;index++) m 255 for(size_t index=0;index<number;index++) m_ws[index] *= a_scale;
256 m_Sw *= a_scale; 256 m_Sw *= a_scale;
257 m_Sxw *= a_scale; 257 m_Sxw *= a_scale;
258 m_Sx2w *= a_scale; 258 m_Sx2w *= a_scale;
259 m_Syw *= a_scale; 259 m_Syw *= a_scale;
260 m_Sy2w *= a_scale; 260 m_Sy2w *= a_scale;
261 return true; 261 return true;
262 } 262 }
263 } 263 }
264 264
265 inline 265 inline
266 bool c2d::reset() { 266 bool c2d::reset() {
267 clear(); 267 clear();
268 delete m_histo; 268 delete m_histo;
269 m_histo = 0; 269 m_histo = 0;
270 return true; 270 return true;
271 } 271 }
272 272
273 inline 273 inline
274 bool c2d::fill(double aX,double aY,double aW){ 274 bool c2d::fill(double aX,double aY,double aW){
275 if(!m_histo && (m_limit!=UNLIMITED()) && ((i 275 if(!m_histo && (m_limit!=UNLIMITED()) && ((int)m_xs.size()>=m_limit)){
276 convert_to_histogram(); 276 convert_to_histogram();
277 } 277 }
278 278
279 if(m_histo) { 279 if(m_histo) {
280 return m_histo->fill(aX,aY,aW); 280 return m_histo->fill(aX,aY,aW);
281 } else { 281 } else {
282 if(m_xs.size()) { 282 if(m_xs.size()) {
283 m_lower_x = mn<double>(aX,m_lower_x); 283 m_lower_x = mn<double>(aX,m_lower_x);
284 m_upper_x = mx<double>(aX,m_upper_x); 284 m_upper_x = mx<double>(aX,m_upper_x);
285 } else { 285 } else {
286 m_lower_x = aX; 286 m_lower_x = aX;
287 m_upper_x = aX; 287 m_upper_x = aX;
288 } 288 }
289 if(m_ys.size()) { 289 if(m_ys.size()) {
290 m_lower_y = mn<double>(aY,m_lower_y); 290 m_lower_y = mn<double>(aY,m_lower_y);
291 m_upper_y = mx<double>(aY,m_upper_y); 291 m_upper_y = mx<double>(aY,m_upper_y);
292 } else { 292 } else {
293 m_lower_y = aY; 293 m_lower_y = aY;
294 m_upper_y = aY; 294 m_upper_y = aY;
295 } 295 }
296 m_xs.push_back(aX); 296 m_xs.push_back(aX);
297 m_ys.push_back(aY); 297 m_ys.push_back(aY);
298 m_ws.push_back(aW); 298 m_ws.push_back(aW);
299 m_Sw += aW; 299 m_Sw += aW;
300 double xw = aX * aW; 300 double xw = aX * aW;
301 m_Sxw += xw; 301 m_Sxw += xw;
302 m_Sx2w += aX * xw; 302 m_Sx2w += aX * xw;
303 double yw = aY * aW; 303 double yw = aY * aW;
304 m_Syw += yw; 304 m_Syw += yw;
305 m_Sy2w += aY * yw; 305 m_Sy2w += aY * yw;
306 return true; 306 return true;
307 } 307 }
308 } 308 }
309 309
310 inline 310 inline
311 bool c2d::convert(const std::vector<double>& a 311 bool c2d::convert(const std::vector<double>& a_edges_x,const std::vector<double>& a_edges_y) {
312 if(m_histo) return true; 312 if(m_histo) return true;
313 m_histo = new histo::h2d(base_cloud::title() 313 m_histo = new histo::h2d(base_cloud::title(),
314 a_edges_x,a_edges_y 314 a_edges_x,a_edges_y);
315 if(!m_histo) return false; 315 if(!m_histo) return false;
316 bool status = fill_histogram(*m_histo); 316 bool status = fill_histogram(*m_histo);
317 clear(); 317 clear();
318 return status; 318 return status;
319 } 319 }
320 320
321 inline 321 inline
322 bool c2d::set_conversion_parameters( 322 bool c2d::set_conversion_parameters(
323 unsigned int aCnvXnumber,double aCnvXmin,doub 323 unsigned int aCnvXnumber,double aCnvXmin,double aCnvXmax
324 ,unsigned int aCnvYnumber,double aCnvYmin,doub 324 ,unsigned int aCnvYnumber,double aCnvYmin,double aCnvYmax
325 ){ 325 ){
326 m_cnv_x_num = aCnvXnumber; 326 m_cnv_x_num = aCnvXnumber;
327 m_cnv_x_min = aCnvXmin; 327 m_cnv_x_min = aCnvXmin;
328 m_cnv_x_max = aCnvXmax; 328 m_cnv_x_max = aCnvXmax;
329 m_cnv_y_num = aCnvYnumber; 329 m_cnv_y_num = aCnvYnumber;
330 m_cnv_y_min = aCnvYmin; 330 m_cnv_y_min = aCnvYmin;
331 m_cnv_y_max = aCnvYmax; 331 m_cnv_y_max = aCnvYmax;
332 return true; 332 return true;
333 } 333 }
334 334
335 inline 335 inline
336 const h2d& c2d::histogram() const { 336 const h2d& c2d::histogram() const {
337 if(!m_histo) const_cast<c2d&>(*this).convert 337 if(!m_histo) const_cast<c2d&>(*this).convert_to_histogram();
338 return *m_histo; 338 return *m_histo;
339 } 339 }
340 340
341 inline 341 inline
342 double c2d::rms_x() const { 342 double c2d::rms_x() const {
343 double rms = 0; //FIXME nan. 343 double rms = 0; //FIXME nan.
344 if(m_histo) { 344 if(m_histo) {
345 rms = m_histo->rms_x(); 345 rms = m_histo->rms_x();
346 } else { 346 } else {
347 if(m_Sw==0) { 347 if(m_Sw==0) {
348 } else { 348 } else {
349 double mean = m_Sxw / m_Sw; 349 double mean = m_Sxw / m_Sw;
350 rms = ::sqrt(::fabs( (m_Sx2w / m_Sw) - m 350 rms = ::sqrt(::fabs( (m_Sx2w / m_Sw) - mean * mean));
351 } 351 }
352 } 352 }
353 return rms; 353 return rms;
354 } 354 }
355 inline 355 inline
356 double c2d::rms_y() const { 356 double c2d::rms_y() const {
357 double rms = 0; //FIXME nan. 357 double rms = 0; //FIXME nan.
358 if(m_histo) { 358 if(m_histo) {
359 rms = m_histo->rms_y(); 359 rms = m_histo->rms_y();
360 } else { 360 } else {
361 if(m_Sw==0) { 361 if(m_Sw==0) {
362 } else { 362 } else {
363 double mean = m_Syw / m_Sw; 363 double mean = m_Syw / m_Sw;
364 rms = ::sqrt(::fabs( (m_Sy2w / m_Sw) - m 364 rms = ::sqrt(::fabs( (m_Sy2w / m_Sw) - mean * mean));
365 } 365 }
366 } 366 }
367 return rms; 367 return rms;
368 } 368 }
369 369
370 }} 370 }}
371 371
372 #endif 372 #endif