Changeset 2608 in Sophya for trunk/SophyaLib/HiStats/histerr.cc

Timestamp:

Sep 7, 2004, 12:50:35 PM (21 years ago)

Author:

cmv

Message:

correct. erreur persistance HistoErr, methode ReCenterBin cmv (sur les conseils de rz) 7/9/04

File:

• trunk/SophyaLib/HiStats/histerr.cc (modified) (8 diffs)

trunk/SophyaLib/HiStats/histerr.cc

@@ -23 +23 @@
 /*! Constructeur d'un histo de nBin bins allant de xMin a xMax */
 HistoErr::HistoErr(r_8 xMin, r_8 xMax, int_4 nBin)
-: Histo(xMin,xMax,nBin), mNData((nBin>0) ? new uint_4[nBin] : NULL)
+: Histo(xMin,xMax,nBin), mNData(NULL)
 {
  this->Errors();
- if(mNData) memset(mNData,0,mBins*sizeof(uint_4));
+ allocate_mNData(nBin);
+ Zero();
 }
 
@@ -32 +33 @@
 /*! Constructeur par copie */
 HistoErr::HistoErr(const HistoErr& H)
-: Histo(H), mNData((H.mBins>0)? new uint_4[H.mBins] : NULL)
-{
+: Histo(H), mNData(NULL)
+{
+ allocate_mNData(H.mBins);
  if(mBins>0) memcpy(mNData,H.mNData,mBins*sizeof(uint_4));
 }
@@ -42 +44 @@
 {
  if(mNData) {delete [] mNData; mNData=NULL;}
+}
+/********* Methode *********/
+/*! Allocation du tableau mNData */
+void HistoErr::allocate_mNData(int nbin)
+{
+ if(nbin<=0) return;
+ if(mNData) {delete [] mNData; mNData=NULL;}
+ mNData = new uint_4[nbin];
 }
 
@@ -126 +136 @@
 }
 
+/********* Methode *********/
+/*!
+ Recompute XMin and XMax so that
+ the center of the first bin is exactly xcmin and
+ the center of the last bin is exactly xcmax
+*/
+void HistoErr::ReCenterBin(r_8 xcmin,r_8 xcmax)
+{
+ if(xcmax<=xcmin) return;
+ if(mBins<=1) return;
+ double dx = (mMax-mMin)/(mBins-1);
+  mMin -= dx/2.;
+  mMax += dx/2.;
+  binWidth = (mMax-mMin)/mBins;
+}
 
 /********* Methode *********/
@@ -155 +180 @@
 void ObjFileIO<HistoErr>::ReadSelf(PInPersist& is)
 {
-char strg[256];
+string strg;
 
 if(dobj==NULL) dobj=new HistoErr;
@@ -161 +186 @@
 
 // Lecture entete
-is.GetLine(strg, 255);
+is.GetStr(strg);
 
 // Lecture des donnees HistoErr
 is.Get(dobj->mBins);
-is.GetLine(strg, 255);
+  dobj->allocate_mNData(dobj->mBins);
 is.Get(dobj->mNData, dobj->mBins);
 
@@ -178 +203 @@
 {
 if(dobj == NULL)   return;
-char strg[256];
-
-sprintf(strg,"HistErr");
-os.PutLine(strg);
+string strg;
+
+strg = "HistErr";
+os.PutStr(strg);
 
 // Ecriture des valeurs
@@ -187 +212 @@
 
 // Ecriture des donnees HistoErr nombre d entree par bin
-sprintf(strg,"Tableau du nombre d entree par bin");
-os.PutLine(strg);
 os.Put(dobj->mNData, dobj->mBins);