1 | #include "sopnamsp.h"
|
---|
2 | #include "machdefs.h"
|
---|
3 | #include <stdlib.h>
|
---|
4 | #include <typeinfo>
|
---|
5 | #include <iostream>
|
---|
6 | #include <string>
|
---|
7 |
|
---|
8 | #include "nomhistadapter.h"
|
---|
9 | #include "pihisto.h"
|
---|
10 | #include "sohiswrap.h"
|
---|
11 |
|
---|
12 | #include "pihisto2d.h"
|
---|
13 | #include "pipodrw.h"
|
---|
14 |
|
---|
15 | #include "servnobjm.h"
|
---|
16 | #include "strutilxx.h"
|
---|
17 |
|
---|
18 | #ifndef SANS_EVOLPLANCK
|
---|
19 | #include "objfitter.h"
|
---|
20 | // Pour les DataTable ( Depuis Avril 2005 )
|
---|
21 | #include "datatable.h"
|
---|
22 | #include "swppfdtable.h"
|
---|
23 | #endif
|
---|
24 |
|
---|
25 | //-----------------------------------------------------------------------------
|
---|
26 | // Class Adaptateur d'objet (Pour NamedObjMgr) d'objet Histo / HProf
|
---|
27 | //-----------------------------------------------------------------------------
|
---|
28 |
|
---|
29 | /* --Methode-- */
|
---|
30 | NOMAdapter_Histo::NOMAdapter_Histo(Histo* o)
|
---|
31 | : NObjMgrAdapter(o)
|
---|
32 | {
|
---|
33 | mHis = o;
|
---|
34 | }
|
---|
35 |
|
---|
36 | /* --Methode-- */
|
---|
37 | NOMAdapter_Histo::~NOMAdapter_Histo()
|
---|
38 | {
|
---|
39 | }
|
---|
40 |
|
---|
41 | /* --Methode-- */
|
---|
42 | NObjMgrAdapter* NOMAdapter_Histo::Clone(AnyDataObj* o)
|
---|
43 | {
|
---|
44 | Histo* h = dynamic_cast<Histo *>(o);
|
---|
45 | if (h) return ( new NOMAdapter_Histo(h) );
|
---|
46 | return ( new NObjMgrAdapter(o) );
|
---|
47 | }
|
---|
48 |
|
---|
49 | /* --Methode-- */
|
---|
50 | string NOMAdapter_Histo::GetDataObjType()
|
---|
51 | {
|
---|
52 | HProf * hp = dynamic_cast<HProf *>(mHis);
|
---|
53 | if(hp) return("HProf ");
|
---|
54 | else return("Histo ");
|
---|
55 | }
|
---|
56 |
|
---|
57 | /* --Methode-- */
|
---|
58 | string NOMAdapter_Histo::GetInfoString(int lev)
|
---|
59 | {
|
---|
60 | char buff[128];
|
---|
61 | string rs;
|
---|
62 | if (lev > 2) {
|
---|
63 | sprintf(buff, "Histo: NBin=%d XMin=%lg XMax=%lg\n", mHis->NBins(),
|
---|
64 | mHis->XMin(), mHis->XMax());
|
---|
65 | rs += buff;
|
---|
66 | }
|
---|
67 | sprintf(buff, "Entries= %llu \n",mHis->NEntries());
|
---|
68 | rs += buff;
|
---|
69 | if (lev > 0) {
|
---|
70 | sprintf(buff, "Overflow= %lg \n", mHis->NOver());
|
---|
71 | rs += buff;
|
---|
72 | sprintf(buff, "Underflow= %lg \n", mHis->NUnder());
|
---|
73 | rs += buff;
|
---|
74 | }
|
---|
75 | sprintf(buff, "Mean= %lg \n", mHis->Mean());
|
---|
76 | rs += buff;
|
---|
77 | sprintf(buff, "Sigma= %lg \n", mHis->Sigma());
|
---|
78 | rs += buff;
|
---|
79 | return rs;
|
---|
80 | }
|
---|
81 |
|
---|
82 | /* --Methode-- */
|
---|
83 | string NOMAdapter_Histo::GetInfoString(vector<string>& opts)
|
---|
84 | {
|
---|
85 | string blabla = "Histo1D: nbin binw mean sigma over under nentries ndata";
|
---|
86 | blabla += " xmin xmax vmin vmax imin imax sum sum2";
|
---|
87 |
|
---|
88 | if(opts.size() == 0) return GetInfoString(3);
|
---|
89 |
|
---|
90 | char buff[64];
|
---|
91 | if(opts[0] == "nbin") {
|
---|
92 | sprintf(buff, "%d",mHis->NBins());
|
---|
93 | } else if(opts[0] == "binw") {
|
---|
94 | sprintf(buff, "%lg",mHis->BinWidth());
|
---|
95 | } else if(opts[0] == "mean") {
|
---|
96 | sprintf(buff, "%lg",mHis->Mean());
|
---|
97 | } else if(opts[0] == "sigma") {
|
---|
98 | sprintf(buff, "%lg",mHis->Sigma());
|
---|
99 | } else if(opts[0] == "over") {
|
---|
100 | sprintf(buff, "%lg",mHis->NOver());
|
---|
101 | } else if(opts[0] == "under") {
|
---|
102 | sprintf(buff, "%lg",mHis->NUnder());
|
---|
103 | } else if(opts[0] == "nentries") {
|
---|
104 | sprintf(buff, "%llu",mHis->NEntries());
|
---|
105 | } else if(opts[0] == "ndata") {
|
---|
106 | sprintf(buff, "%lg",mHis->NData());
|
---|
107 | } else if(opts[0] == "xmin") {
|
---|
108 | sprintf(buff, "%lg",mHis->XMin());
|
---|
109 | } else if(opts[0] == "xmax") {
|
---|
110 | sprintf(buff, "%lg",mHis->XMax());
|
---|
111 | } else if(opts[0] == "vmin") {
|
---|
112 | sprintf(buff, "%lg",mHis->VMin());
|
---|
113 | } else if(opts[0] == "vmax") {
|
---|
114 | sprintf(buff, "%lg",mHis->VMax());
|
---|
115 | } else if(opts[0] == "imin") {
|
---|
116 | sprintf(buff, "%d",mHis->IMin());
|
---|
117 | } else if(opts[0] == "imax") {
|
---|
118 | sprintf(buff, "%d",mHis->IMax());
|
---|
119 | } else if(opts[0] == "sum") {
|
---|
120 | sprintf(buff, "%lg",mHis->Sum());
|
---|
121 | } else if(opts[0] == "sum2") {
|
---|
122 | sprintf(buff, "%lg",mHis->Sum2());
|
---|
123 | } else {
|
---|
124 | return blabla;
|
---|
125 | }
|
---|
126 | return string(buff);
|
---|
127 | }
|
---|
128 |
|
---|
129 | /* --Methode-- */
|
---|
130 | AnyDataObj* NOMAdapter_Histo::CloneDataObj(bool /*share*/)
|
---|
131 | {
|
---|
132 | mHis->UpdateHisto(); // pour le cas ou c'est un HProf
|
---|
133 | HProf * hp = dynamic_cast<HProf *>(mHis);
|
---|
134 | if(hp) return( new HProf(*hp) );
|
---|
135 | else return( new Histo(*mHis) );
|
---|
136 | }
|
---|
137 |
|
---|
138 | /* --Methode-- */
|
---|
139 | void NOMAdapter_Histo::SavePPF(POutPersist& pos, string const & nom)
|
---|
140 | {
|
---|
141 | #ifdef SANS_EVOLPLANCK
|
---|
142 | // PEIDA-EROS L'histo est lui-meme PPersist
|
---|
143 | string tag = nom; // A cause de const
|
---|
144 | mHis->Write(pos,0,tag);
|
---|
145 | #else
|
---|
146 | ObjFileIO<Histo> fio(mHis);
|
---|
147 | fio.Write(pos, nom);
|
---|
148 | #endif
|
---|
149 | }
|
---|
150 |
|
---|
151 | /* --Methode-- */
|
---|
152 | void NOMAdapter_Histo::Print(ostream& os, int lev)
|
---|
153 | {
|
---|
154 | mHis->Show(os);
|
---|
155 | if (lev > 0) mHis->Print(60);
|
---|
156 | }
|
---|
157 |
|
---|
158 | /* --Methode-- */
|
---|
159 | PIDrawer* NOMAdapter_Histo::GetDrawer(string & dopt)
|
---|
160 | {
|
---|
161 | if (typeid(*mHis) == typeid(HProf)) dopt = "fcirclemarker5 " + dopt;
|
---|
162 | else dopt = "thinline " + dopt;
|
---|
163 | HistoWrapper* hw = new HistoWrapper(mHis, false); // false: le Wrapper ne delete pas l'objet Histo mHis
|
---|
164 | PIHisto * pih = new PIHisto(hw, true); // true: PIHisto delete l'objet HistoWrapper hw
|
---|
165 | return( pih );
|
---|
166 | }
|
---|
167 |
|
---|
168 | /* --Methode-- */
|
---|
169 | NTupleInterface* NOMAdapter_Histo::GetNTupleInterface(bool& adel)
|
---|
170 | {
|
---|
171 | adel = true;
|
---|
172 | return( new NTupInt_Histo(mHis) );
|
---|
173 | }
|
---|
174 |
|
---|
175 | /* --Methode-- */
|
---|
176 | GeneralFitData* NOMAdapter_Histo::GetGeneralFitData(bool& adel
|
---|
177 | ,GeneralFitData::FitErrType errtype,double errscale,double errmin
|
---|
178 | ,int i1,int i2,int j1,int j2)
|
---|
179 | {
|
---|
180 | adel = false;
|
---|
181 | if(!mHis) return(NULL);
|
---|
182 |
|
---|
183 | int nx = mHis->NBins();
|
---|
184 | if(nx<=0) return(NULL);
|
---|
185 |
|
---|
186 | i1 = (i1<0||i1>=nx)? 0: i1;
|
---|
187 | i2 = (i2<0||i2>=nx||i2<i1)? nx-1: i2;
|
---|
188 |
|
---|
189 | GeneralFitData* mGData = new GeneralFitData(1,i2-i1+1,0);
|
---|
190 | adel = true;
|
---|
191 |
|
---|
192 | for(int i=i1;i<=i2;i++) {
|
---|
193 | double x = mHis->BinCenter(i);
|
---|
194 | double f = (*mHis)(i);
|
---|
195 | double e = (mHis->HasErrors())? mHis->Error(i) : 1.;
|
---|
196 | e = GeneralFitData::ComputeError(f,e,errtype,errscale,errmin);
|
---|
197 | mGData->AddData1(x,f,e);
|
---|
198 | }
|
---|
199 |
|
---|
200 | return mGData;
|
---|
201 | }
|
---|
202 |
|
---|
203 | AnyDataObj* NOMAdapter_Histo::FitResidusObj(GeneralFit& mfit)
|
---|
204 | {
|
---|
205 | Histo* h = NULL;
|
---|
206 | #ifdef SANS_EVOLPLANCK
|
---|
207 | h = mHis->FitResidus(mfit);
|
---|
208 | #else
|
---|
209 | h = new Histo(ObjectFitter::FitResidus(*mHis,mfit));
|
---|
210 | #endif
|
---|
211 | return h;
|
---|
212 | }
|
---|
213 |
|
---|
214 | AnyDataObj* NOMAdapter_Histo::FitFunctionObj(GeneralFit& mfit)
|
---|
215 | {
|
---|
216 | Histo* h = NULL;
|
---|
217 | #ifdef SANS_EVOLPLANCK
|
---|
218 | h = mHis->FitFunction(mfit);
|
---|
219 | #else
|
---|
220 | h = new Histo(ObjectFitter::FitFunction(*mHis,mfit));
|
---|
221 | #endif
|
---|
222 | return h;
|
---|
223 | }
|
---|
224 |
|
---|
225 | // -------------------------------------------------------------
|
---|
226 |
|
---|
227 | /* --Methode-- */
|
---|
228 | NTupInt_Histo::NTupInt_Histo(Histo* h)
|
---|
229 | {
|
---|
230 | mHis = h;
|
---|
231 | mHpr = dynamic_cast<HProf *>(h);
|
---|
232 | }
|
---|
233 |
|
---|
234 | /* --Methode-- */
|
---|
235 | NTupInt_Histo::~NTupInt_Histo()
|
---|
236 | {
|
---|
237 | }
|
---|
238 |
|
---|
239 | /* --Methode-- */
|
---|
240 | sa_size_t NTupInt_Histo::NbLines() const
|
---|
241 | {
|
---|
242 | return(mHis->NBins());
|
---|
243 | }
|
---|
244 |
|
---|
245 | /* --Methode-- */
|
---|
246 | sa_size_t NTupInt_Histo::NbColumns() const
|
---|
247 | {
|
---|
248 | return(5);
|
---|
249 | }
|
---|
250 |
|
---|
251 | /* --Methode-- */
|
---|
252 | r_8* NTupInt_Histo::GetLineD(sa_size_t k) const
|
---|
253 | {
|
---|
254 | int i;
|
---|
255 | if ((k < 0) || (k >= mHis->NBins()))
|
---|
256 | for(i=0; i<5; i++) mRet[i] = 0.;
|
---|
257 | else {
|
---|
258 | mRet[0] = k; mRet[1] = mHis->BinCenter(k);
|
---|
259 | mRet[2] = (*mHis)(k); mRet[3] = mHis->Error(k);
|
---|
260 | if(mHpr) mRet[4] = mHpr->SumW(k); else mRet[4] = 0.;
|
---|
261 | }
|
---|
262 | return(mRet);
|
---|
263 | }
|
---|
264 |
|
---|
265 | /* --Methode-- */
|
---|
266 | string NTupInt_Histo::VarList_C(const char* nx) const
|
---|
267 | {
|
---|
268 | string nomx;
|
---|
269 | if (nx) nomx = nx;
|
---|
270 | else nomx = "_xh_";
|
---|
271 | string vardec = "double i,x,val,err,nb; \n";
|
---|
272 | vardec += "i = " + nomx + "[0]; x = " + nomx + "[1]; \n";
|
---|
273 | vardec += "val = " + nomx + "[2]; err = " + nomx + "[3]; \n";
|
---|
274 | vardec += "nb = " + nomx + "[4]; \n";
|
---|
275 | return(vardec);
|
---|
276 | }
|
---|
277 |
|
---|
278 | //-------------------------------------------------------------------------
|
---|
279 | // Class Adaptateur d'objet (Pour NamedObjMgr) d'objet Histo2D
|
---|
280 | //-------------------------------------------------------------------------
|
---|
281 |
|
---|
282 |
|
---|
283 | /* --Methode-- */
|
---|
284 | NOMAdapter_Histo2D::NOMAdapter_Histo2D(Histo2D* o)
|
---|
285 | : NObjMgrAdapter(o)
|
---|
286 | {
|
---|
287 | mHis = o;
|
---|
288 | }
|
---|
289 |
|
---|
290 | /* --Methode-- */
|
---|
291 | NOMAdapter_Histo2D::~NOMAdapter_Histo2D()
|
---|
292 | {
|
---|
293 | }
|
---|
294 |
|
---|
295 | /* --Methode-- */
|
---|
296 | NObjMgrAdapter* NOMAdapter_Histo2D::Clone(AnyDataObj* o)
|
---|
297 | {
|
---|
298 | Histo2D* h = dynamic_cast<Histo2D *>(o);
|
---|
299 | if (h) return ( new NOMAdapter_Histo2D(h) );
|
---|
300 | return ( new NObjMgrAdapter(o) );
|
---|
301 | }
|
---|
302 |
|
---|
303 | /* --Methode-- */
|
---|
304 | string NOMAdapter_Histo2D::GetDataObjType()
|
---|
305 | {
|
---|
306 | return ("Histo2D ");
|
---|
307 | }
|
---|
308 |
|
---|
309 | /* --Methode-- */
|
---|
310 | AnyDataObj* NOMAdapter_Histo2D::CloneDataObj(bool /*share*/)
|
---|
311 | {
|
---|
312 | return ( new Histo2D(*mHis) );
|
---|
313 | }
|
---|
314 |
|
---|
315 | /* --Methode-- */
|
---|
316 | void NOMAdapter_Histo2D::SavePPF(POutPersist& pos, string const & nom)
|
---|
317 | {
|
---|
318 | #ifdef SANS_EVOLPLANCK
|
---|
319 | // PEIDA-EROS L'histo est lui-meme PPersist
|
---|
320 | string tag = nom; // A cause de const
|
---|
321 | mHis->Write(pos,0,tag);
|
---|
322 | #else
|
---|
323 | ObjFileIO<Histo2D> fio(mHis);
|
---|
324 | fio.Write(pos, nom);
|
---|
325 | #endif
|
---|
326 | }
|
---|
327 |
|
---|
328 | /* --Methode-- */
|
---|
329 | void NOMAdapter_Histo2D::Print(ostream& os, int lev)
|
---|
330 | {
|
---|
331 | mHis->Show(os);
|
---|
332 | if (lev > 0) mHis->Print();
|
---|
333 | }
|
---|
334 |
|
---|
335 | /* --Methode-- */
|
---|
336 | PIDrawer* NOMAdapter_Histo2D::GetDrawer(string & dopt)
|
---|
337 | {
|
---|
338 | dopt = "thinline " + dopt;
|
---|
339 | return( new PIHisto2D(new Histo2DWrapper(mHis, false), true) );
|
---|
340 | }
|
---|
341 |
|
---|
342 | /* --Methode-- */
|
---|
343 | P2DArrayAdapter* NOMAdapter_Histo2D::Get2DArray(string & dopt)
|
---|
344 | {
|
---|
345 | Histo2DWrapper* hwp = new Histo2DWrapper(mHis, false);
|
---|
346 | return hwp;
|
---|
347 | }
|
---|
348 |
|
---|
349 | /* --Methode-- */
|
---|
350 | NTupleInterface* NOMAdapter_Histo2D::GetNTupleInterface(bool& adel)
|
---|
351 | {
|
---|
352 | adel = true;
|
---|
353 | return( new NTupInt_Histo2D(mHis) );
|
---|
354 | }
|
---|
355 |
|
---|
356 |
|
---|
357 | /* --Methode-- */
|
---|
358 | GeneralFitData* NOMAdapter_Histo2D::GetGeneralFitData(bool& adel
|
---|
359 | ,GeneralFitData::FitErrType errtype,double errscale,double errmin
|
---|
360 | ,int i1,int i2,int j1,int j2)
|
---|
361 | {
|
---|
362 | adel = false;
|
---|
363 | if(!mHis) return(NULL);
|
---|
364 |
|
---|
365 | int nx = mHis->NBinX();
|
---|
366 | int ny = mHis->NBinY();
|
---|
367 | if(nx<=0 || ny<=0) return(NULL);
|
---|
368 |
|
---|
369 | i1 = (i1<0||i1>=nx)? 0: i1;
|
---|
370 | i2 = (i2<0||i2>=nx||i2<i1)? nx-1: i2;
|
---|
371 | j1 = (j1<0||j1>=ny)? 0: j1;
|
---|
372 | j2 = (j2<0||j2>=ny||j2<j1)? ny-1: j2;
|
---|
373 |
|
---|
374 | GeneralFitData* mGData = new GeneralFitData(2,(i2-i1+1)*(j2-j1+1),0);
|
---|
375 | adel = true;
|
---|
376 |
|
---|
377 | for(int i=i1;i<=i2;i++) for(int j=j1;j<=j2;j++) {
|
---|
378 | double x,y; mHis->BinCenter(i,j,x,y);
|
---|
379 | double f = (*mHis)(i,j);
|
---|
380 | double e = (mHis->HasErrors())? mHis->Error(i,j) : 1.;
|
---|
381 | e = GeneralFitData::ComputeError(f,e,errtype,errscale,errmin);
|
---|
382 | mGData->AddData2(x,y,f,e);
|
---|
383 | }
|
---|
384 |
|
---|
385 | return mGData;
|
---|
386 | }
|
---|
387 |
|
---|
388 | AnyDataObj* NOMAdapter_Histo2D::FitResidusObj(GeneralFit& mfit)
|
---|
389 | {
|
---|
390 | Histo2D* h2 = NULL;
|
---|
391 | #ifdef SANS_EVOLPLANCK
|
---|
392 | h2 = mHis->FitFunction(mfit);
|
---|
393 | #else
|
---|
394 | h2 = new Histo2D(ObjectFitter::FitResidus(*mHis,mfit));
|
---|
395 | #endif
|
---|
396 | return h2;
|
---|
397 | }
|
---|
398 |
|
---|
399 | AnyDataObj* NOMAdapter_Histo2D::FitFunctionObj(GeneralFit& mfit)
|
---|
400 | {
|
---|
401 | Histo2D* h2 = NULL;
|
---|
402 | #ifdef SANS_EVOLPLANCK
|
---|
403 | h2 = mHis->FitFunction(mfit);
|
---|
404 | #else
|
---|
405 | h2 = new Histo2D(ObjectFitter::FitFunction(*mHis,mfit));
|
---|
406 | #endif
|
---|
407 | return h2;
|
---|
408 | }
|
---|
409 |
|
---|
410 |
|
---|
411 | /* --Methode-- */
|
---|
412 | string NOMAdapter_Histo2D::GetInfoString(vector<string>& opts)
|
---|
413 | {
|
---|
414 | string blabla = "Histo2D: nbin binw nband nslice nentries ndata xmin xmax";
|
---|
415 | blabla += " ymin ymax vmin vmax ijmin ijmax";
|
---|
416 |
|
---|
417 | if(opts.size() == 0) return blabla;
|
---|
418 |
|
---|
419 | char buff[128];
|
---|
420 | if(opts[0] == "nbin") {
|
---|
421 | sprintf(buff, "%d %d",mHis->NBinX(),mHis->NBinY());
|
---|
422 | } else if(opts[0] == "binw") {
|
---|
423 | sprintf(buff, "%lg %lg",mHis->WBinX(),mHis->WBinY());
|
---|
424 | } else if(opts[0] == "nband") {
|
---|
425 | sprintf(buff, "%d %d",mHis->NBandX(),mHis->NBandY());
|
---|
426 | } else if(opts[0] == "nslice") {
|
---|
427 | sprintf(buff, "%d %d",mHis->NSliX(),mHis->NSliY());
|
---|
428 | } else if(opts[0] == "nentries") {
|
---|
429 | sprintf(buff, "%d",mHis->NEntries());
|
---|
430 | } else if(opts[0] == "ndata") {
|
---|
431 | sprintf(buff, "%lg",mHis->NData());
|
---|
432 | } else if(opts[0] == "xmin") {
|
---|
433 | sprintf(buff, "%lg",mHis->XMin());
|
---|
434 | } else if(opts[0] == "xmax") {
|
---|
435 | sprintf(buff, "%lg",mHis->XMax());
|
---|
436 | } else if(opts[0] == "ymin") {
|
---|
437 | sprintf(buff, "%lg",mHis->YMin());
|
---|
438 | } else if(opts[0] == "ymax") {
|
---|
439 | sprintf(buff, "%lg",mHis->YMax());
|
---|
440 | } else if(opts[0] == "vmin") {
|
---|
441 | sprintf(buff, "%lg",mHis->VMin());
|
---|
442 | } else if(opts[0] == "vmax") {
|
---|
443 | sprintf(buff, "%lg",mHis->VMax());
|
---|
444 | } else if(opts[0] == "ijmin") {
|
---|
445 | int_4 i,j;
|
---|
446 | mHis->IJMin(i,j);
|
---|
447 | sprintf(buff, "%d %d",i,j);
|
---|
448 | } else if(opts[0] == "ijmax") {
|
---|
449 | int_4 i,j;
|
---|
450 | mHis->IJMax(i,j);
|
---|
451 | sprintf(buff, "%d %d",i,j);
|
---|
452 | } else {
|
---|
453 | return blabla;
|
---|
454 | }
|
---|
455 | return string(buff);
|
---|
456 | }
|
---|
457 |
|
---|
458 | // -------------------------------------------------------------
|
---|
459 |
|
---|
460 | /* --Methode-- */
|
---|
461 | NTupInt_Histo2D::NTupInt_Histo2D(Histo2D* h)
|
---|
462 | {
|
---|
463 | mHis = h;
|
---|
464 | }
|
---|
465 |
|
---|
466 | /* --Methode-- */
|
---|
467 | NTupInt_Histo2D::~NTupInt_Histo2D()
|
---|
468 | {
|
---|
469 | }
|
---|
470 |
|
---|
471 | /* --Methode-- */
|
---|
472 | sa_size_t NTupInt_Histo2D::NbLines() const
|
---|
473 | {
|
---|
474 | return(mHis->NBinX()*mHis->NBinY());
|
---|
475 | }
|
---|
476 |
|
---|
477 | /* --Methode-- */
|
---|
478 | sa_size_t NTupInt_Histo2D::NbColumns() const
|
---|
479 | {
|
---|
480 | return(6);
|
---|
481 | }
|
---|
482 |
|
---|
483 | /* --Methode-- */
|
---|
484 | r_8* NTupInt_Histo2D::GetLineD(sa_size_t n) const
|
---|
485 | {
|
---|
486 | int i,j;
|
---|
487 | r_8 f2,f3;
|
---|
488 | if ((n < 0) || (n >= mHis->NBinX()*mHis->NBinY()))
|
---|
489 | for(i=0; i<6; i++) mRet[i] = 0.;
|
---|
490 | else {
|
---|
491 | i = n%mHis->NBinX(); j = n/mHis->NBinX();
|
---|
492 | mRet[0] = i; mRet[1] = j;
|
---|
493 | mHis->BinCenter(i,j,f2,f3);
|
---|
494 | mRet[2] = f2; mRet[3] = f3;
|
---|
495 | mRet[4] = (*mHis)(i,j); mRet[5] = mHis->Error(i, j);
|
---|
496 | }
|
---|
497 | return(mRet);
|
---|
498 | }
|
---|
499 |
|
---|
500 | /* --Methode-- */
|
---|
501 | string NTupInt_Histo2D::VarList_C(const char* nx) const
|
---|
502 | {
|
---|
503 | string nomx;
|
---|
504 | if (nx) nomx = nx;
|
---|
505 | else nomx = "_xh_";
|
---|
506 | string vardec = "double i,j,x,y,val,err; \n";
|
---|
507 | vardec += "i = " + nomx + "[0]; j = " + nomx + "[1]; \n";
|
---|
508 | vardec += "x = " + nomx + "[2]; y = " + nomx + "[3]; \n";
|
---|
509 | vardec += "val = " + nomx + "[4]; err = " + nomx + "[5]; \n";
|
---|
510 | return(vardec);
|
---|
511 | }
|
---|
512 |
|
---|
513 |
|
---|
514 |
|
---|
515 | //-------------------------------------------------------------------------
|
---|
516 | // Class Adaptateur d'objet (Pour NamedObjMgr) d'objet NTuple
|
---|
517 | //-------------------------------------------------------------------------
|
---|
518 |
|
---|
519 | /* --Methode-- */
|
---|
520 | NOMAdapter_NTuple::NOMAdapter_NTuple(NTuple* o)
|
---|
521 | : NObjMgrAdapter(o)
|
---|
522 | {
|
---|
523 | mNt = o;
|
---|
524 | }
|
---|
525 |
|
---|
526 | /* --Methode-- */
|
---|
527 | NOMAdapter_NTuple::~NOMAdapter_NTuple()
|
---|
528 | {
|
---|
529 | }
|
---|
530 |
|
---|
531 | /* --Methode-- */
|
---|
532 | NObjMgrAdapter* NOMAdapter_NTuple::Clone(AnyDataObj* o)
|
---|
533 | {
|
---|
534 | NTuple* nt = dynamic_cast<NTuple *>(o);
|
---|
535 | if (nt) return ( new NOMAdapter_NTuple(nt) );
|
---|
536 | return ( new NObjMgrAdapter(o) );
|
---|
537 | }
|
---|
538 |
|
---|
539 |
|
---|
540 | /* --Methode-- */
|
---|
541 | string NOMAdapter_NTuple::GetDataObjType()
|
---|
542 | {
|
---|
543 | return ("NTuple ");
|
---|
544 | }
|
---|
545 |
|
---|
546 | /* --Methode-- */
|
---|
547 | AnyDataObj* NOMAdapter_NTuple::CloneDataObj(bool /*share*/)
|
---|
548 | {
|
---|
549 | return ( new NTuple(*mNt) );
|
---|
550 | }
|
---|
551 |
|
---|
552 |
|
---|
553 |
|
---|
554 | /* --Methode-- */
|
---|
555 | void NOMAdapter_NTuple::SavePPF(POutPersist& pos, string const & nom)
|
---|
556 | {
|
---|
557 | #ifdef SANS_EVOLPLANCK
|
---|
558 | // PEIDA-EROS L'histo est lui-meme PPersist
|
---|
559 | string tag = nom; // A cause de const
|
---|
560 | mNt->Write(pos,0,tag);
|
---|
561 | #else
|
---|
562 | ObjFileIO<NTuple> fio(mNt);
|
---|
563 | fio.Write(pos, nom);
|
---|
564 | #endif
|
---|
565 | }
|
---|
566 |
|
---|
567 | /* --Methode-- */
|
---|
568 | string NOMAdapter_NTuple::GetInfoString(vector<string>& opts)
|
---|
569 | {
|
---|
570 | if (opts.size() == 0) return NObjMgrAdapter::GetInfoString(opts);
|
---|
571 | char buff[128];
|
---|
572 | if (opts[0] == "sizes") {
|
---|
573 | sprintf(buff, "%ld %ld", mNt->NEntry(), mNt->NVar());
|
---|
574 | return string(buff);
|
---|
575 | }
|
---|
576 | else if ((opts[0] == "nlines") || (opts[0] == "nentry") || (opts[0] == "nrows")) {
|
---|
577 | sprintf(buff, "%ld", mNt->NEntry());
|
---|
578 | return string(buff);
|
---|
579 | }
|
---|
580 | else if ((opts[0] == "nvar") || (opts[0] == "ncols")) {
|
---|
581 | sprintf(buff, "%ld", mNt->NVar());
|
---|
582 | return string(buff);
|
---|
583 | }
|
---|
584 | else return "NTuple.Att: nlines/nentry/nrows nvar/ncols";
|
---|
585 | }
|
---|
586 |
|
---|
587 | /* --Methode-- */
|
---|
588 | void NOMAdapter_NTuple::Print(ostream& os, int lev)
|
---|
589 | {
|
---|
590 | mNt->Show(os);
|
---|
591 | if (lev > 2) os << mNt->Info();
|
---|
592 | if (lev > 0) mNt->Print(0, 10*lev);
|
---|
593 | }
|
---|
594 |
|
---|
595 |
|
---|
596 | /* --Methode-- */
|
---|
597 | NTupleInterface* NOMAdapter_NTuple::GetNTupleInterface(bool& adel)
|
---|
598 | {
|
---|
599 | adel = false;
|
---|
600 | return(mNt);
|
---|
601 | // return( new NTupInt_NTuple(mNt) );
|
---|
602 | }
|
---|
603 |
|
---|
604 | //-------------------------------------------------------------------------
|
---|
605 | // Class Adaptateur d'objet (Pour NamedObjMgr) d'objet XNTuple
|
---|
606 | //-------------------------------------------------------------------------
|
---|
607 |
|
---|
608 | /* --Methode-- */
|
---|
609 | NOMAdapter_XNTuple::NOMAdapter_XNTuple(XNTuple* o)
|
---|
610 | : NObjMgrAdapter(o)
|
---|
611 | {
|
---|
612 | mNt = o;
|
---|
613 | }
|
---|
614 |
|
---|
615 | /* --Methode-- */
|
---|
616 | NOMAdapter_XNTuple::~NOMAdapter_XNTuple()
|
---|
617 | {
|
---|
618 | }
|
---|
619 |
|
---|
620 | /* --Methode-- */
|
---|
621 | NObjMgrAdapter* NOMAdapter_XNTuple::Clone(AnyDataObj* o)
|
---|
622 | {
|
---|
623 | XNTuple* nt = dynamic_cast<XNTuple *>(o);
|
---|
624 | if (nt) return ( new NOMAdapter_XNTuple(nt) );
|
---|
625 | return ( new NObjMgrAdapter(o) );
|
---|
626 | }
|
---|
627 |
|
---|
628 | /* --Methode-- */
|
---|
629 | string NOMAdapter_XNTuple::GetDataObjType()
|
---|
630 | {
|
---|
631 | return ("XNTuple ");
|
---|
632 | }
|
---|
633 |
|
---|
634 | /* --Methode-- */
|
---|
635 | AnyDataObj* NOMAdapter_XNTuple::CloneDataObj(bool /*share*/)
|
---|
636 | {
|
---|
637 | return ( new XNTuple(*mNt) );
|
---|
638 | }
|
---|
639 |
|
---|
640 | /* --Methode-- */
|
---|
641 | void NOMAdapter_XNTuple::SavePPF(POutPersist& pos, string const & nom)
|
---|
642 | {
|
---|
643 | #ifdef SANS_EVOLPLANCK
|
---|
644 | // PEIDA-EROS L'histo est lui-meme PPersist
|
---|
645 | string tag = nom; // A cause de const
|
---|
646 | mNt->Write(pos,0,tag);
|
---|
647 | #else
|
---|
648 | ObjFileIO<XNTuple> fio(mNt);
|
---|
649 | fio.Write(pos, nom);
|
---|
650 | #endif
|
---|
651 | }
|
---|
652 |
|
---|
653 | /* --Methode-- */
|
---|
654 | void NOMAdapter_XNTuple::Print(ostream& os, int lev)
|
---|
655 | {
|
---|
656 | // os << mNt->Info();
|
---|
657 | mNt->Show(os);
|
---|
658 | }
|
---|
659 |
|
---|
660 |
|
---|
661 | /* --Methode-- */
|
---|
662 | NTupleInterface* NOMAdapter_XNTuple::GetNTupleInterface(bool& adel)
|
---|
663 | {
|
---|
664 | adel = false;
|
---|
665 | return(mNt);
|
---|
666 | }
|
---|
667 |
|
---|
668 |
|
---|
669 | //-------------------------------------------------------------------------
|
---|
670 | // Class Adaptateur d'objet (Pour NamedObjMgr) d'objet BaseDataTable
|
---|
671 | //-------------------------------------------------------------------------
|
---|
672 |
|
---|
673 | /* --Methode-- */
|
---|
674 | NOMAdapter_DataTable::NOMAdapter_DataTable(BaseDataTable* o)
|
---|
675 | : NObjMgrAdapter(o)
|
---|
676 | {
|
---|
677 | mDT = o;
|
---|
678 | }
|
---|
679 |
|
---|
680 | /* --Methode-- */
|
---|
681 | NOMAdapter_DataTable::~NOMAdapter_DataTable()
|
---|
682 | {
|
---|
683 | }
|
---|
684 |
|
---|
685 | /* --Methode-- */
|
---|
686 | NObjMgrAdapter* NOMAdapter_DataTable::Clone(AnyDataObj* o)
|
---|
687 | {
|
---|
688 | BaseDataTable* dt = dynamic_cast<BaseDataTable *>(o);
|
---|
689 | if (dt) return ( new NOMAdapter_DataTable(dt) );
|
---|
690 | else return ( new NObjMgrAdapter(o) );
|
---|
691 | }
|
---|
692 |
|
---|
693 | /* --Methode-- */
|
---|
694 | string NOMAdapter_DataTable::GetDataObjType()
|
---|
695 | {
|
---|
696 | DataTable* dt = dynamic_cast<DataTable *>(mDT);
|
---|
697 | if (dt) return ("DataTable ");
|
---|
698 | else {
|
---|
699 | SwPPFDataTable* swdt = dynamic_cast<SwPPFDataTable *>(mDT);
|
---|
700 | if (swdt) return ("SwPPFDataTable ");
|
---|
701 | return ("BaseDataTable ");
|
---|
702 | }
|
---|
703 | }
|
---|
704 |
|
---|
705 | /* --Methode-- */
|
---|
706 | AnyDataObj* NOMAdapter_DataTable::CloneDataObj(bool share)
|
---|
707 | {
|
---|
708 | DataTable* dt = dynamic_cast<DataTable *>(mDT);
|
---|
709 | if (dt) return new DataTable(*dt, share);
|
---|
710 | else {
|
---|
711 | SwPPFDataTable* swdt = dynamic_cast<SwPPFDataTable *>(mDT);
|
---|
712 | if (swdt) cout << "NOMAdapter_DataTable::CloneDataObj() Object type SwPPFDataTable can not be cloned !" << endl;
|
---|
713 | return NULL;
|
---|
714 | }
|
---|
715 | }
|
---|
716 |
|
---|
717 |
|
---|
718 | /* --Methode-- */
|
---|
719 | void NOMAdapter_DataTable::SavePPF(POutPersist& pos, string const & nom)
|
---|
720 | {
|
---|
721 | DataTable* dt = dynamic_cast<DataTable *>(mDT);
|
---|
722 | SwPPFDataTable* swdt = dynamic_cast<SwPPFDataTable *>(mDT);
|
---|
723 | if (dt) {
|
---|
724 | ObjFileIO<BaseDataTable> fio(dt);
|
---|
725 | fio.Write(pos, nom);
|
---|
726 | }
|
---|
727 | else if (swdt) {
|
---|
728 | ObjFileIO<BaseDataTable> fio(swdt);
|
---|
729 | fio.Write(pos, nom);
|
---|
730 | }
|
---|
731 | else {
|
---|
732 | cerr << " NOMAdapter_DataTable::SavePPF() Objet pas de type DataTable/SwPPFDataTable (nom="
|
---|
733 | << nom << ")" << endl;
|
---|
734 | }
|
---|
735 | }
|
---|
736 |
|
---|
737 | /* --Methode-- */
|
---|
738 | string NOMAdapter_DataTable::GetInfoString(vector<string>& opts)
|
---|
739 | {
|
---|
740 | if (opts.size() == 0) return NObjMgrAdapter::GetInfoString(opts);
|
---|
741 | char buff[128];
|
---|
742 | if (opts[0] == "sizes") {
|
---|
743 | sprintf(buff, "%ld %ld", mDT->NEntry(), mDT->NVar());
|
---|
744 | return string(buff);
|
---|
745 | }
|
---|
746 | else if ((opts[0] == "nlines") || (opts[0] == "nentry") || (opts[0] == "nrows")) {
|
---|
747 | sprintf(buff, "%ld", mDT->NEntry());
|
---|
748 | return string(buff);
|
---|
749 | }
|
---|
750 | else if ((opts[0] == "nvar") || (opts[0] == "ncols")) {
|
---|
751 | sprintf(buff, "%ld", mDT->NVar());
|
---|
752 | return string(buff);
|
---|
753 | }
|
---|
754 | else if (opts[0] == "info") { // Acces aux valeurs stockes ds le DVList Info()
|
---|
755 | if (opts.size() < 2) return string("");
|
---|
756 | else if (mDT->Info().HasKey(opts[1]) == false) return string("");
|
---|
757 | else return mDT->Info().GetS(opts[1]);
|
---|
758 | }
|
---|
759 | else if (opts[0] == "row") { // Pour recuperer une ligne de la table
|
---|
760 | if (opts.size() < 2) return string("");
|
---|
761 | long num = atol(opts[1].c_str());
|
---|
762 | return mDT->TableRowToString(num, true);
|
---|
763 | }
|
---|
764 |
|
---|
765 | else return "BaseDataTable.Att: nlines/nentry/nrows nvar/ncols info.varname row.num ";
|
---|
766 | }
|
---|
767 |
|
---|
768 | /* --Methode-- */
|
---|
769 | void NOMAdapter_DataTable::Print(ostream& os, int lev)
|
---|
770 | {
|
---|
771 | mDT->Show(os);
|
---|
772 | if (lev < 1) return;
|
---|
773 | if (lev < 5) mDT->Print(os, 0, lev*10);
|
---|
774 | else mDT->Print(os);
|
---|
775 | }
|
---|
776 |
|
---|
777 |
|
---|
778 | /* --Methode-- */
|
---|
779 | NTupleInterface* NOMAdapter_DataTable::GetNTupleInterface(bool& adel)
|
---|
780 | {
|
---|
781 | adel = false;
|
---|
782 | return(mDT);
|
---|
783 | }
|
---|