| 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 | } | 
|---|