Changeset 2627 in Sophya for trunk/SophyaLib/HiStats

Timestamp:

Oct 6, 2004, 6:04:27 PM (21 years ago)

Author:

cmv

Message:

petit bug dans methode ReCenterBin cmv 6/10/04

Location:

trunk/SophyaLib/HiStats

Files:

• histerr.cc (modified) (4 diffs)
• histerr.h (modified) (1 diff)

trunk/SophyaLib/HiStats/histerr.cc

@@ -84 +84 @@
  else if(numBin>=mBins) mOver += w;
  else {
-   mData[numBin] += w; mNData[numBin]++; mErr2[numBin] += e*e;
+   mData[numBin] += w; mNData[numBin] += 1.; mErr2[numBin] += e*e;
    nHist += w; nEntries++;
  }
@@ -98 +98 @@
  else if(numBin>=mBins) mOver += w;
  else {
-   mData[numBin] += w; mNData[numBin]++; mErr2[numBin] += e*e;
+   mData[numBin] += w; mNData[numBin] += 1.; mErr2[numBin] += e*e;
    nHist += w; nEntries++;
  }
@@ -149 +149 @@
 /*!
  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;
+ the CENTER of the first bin is exactly XMin and
+ the CENTER of the last bin is exactly XMax.
+ Remember that otherwise
+ XMin is the beginning of the first bin
+ and XMax is the end of the last bin
+*/
+void HistoErr::ReCenterBin(void)
+{
  if(mBins<=1) return;
  double dx = (mMax-mMin)/(mBins-1);
-  mMin -= dx/2.;
-  mMax += dx/2.;
-  binWidth = (mMax-mMin)/mBins;
+ mMin -= dx/2.;
+ mMax += dx/2.;
+ binWidth = (mMax-mMin)/mBins;
 }
 
@@ -174 +176 @@
  if(mBins<1) return;
  for(int_4 i=0;i<mBins;i++) {
-   if(mNData[i]<1) continue;
+   if(mNData[i]<1.) continue;
    mData[i] /= mNData[i];
    mErr2[i] /= mNData[i];

trunk/SophyaLib/HiStats/histerr.h

@@ -40 +40 @@
 
   //! Re-center bin abscissa
-  void ReCenterBin(r_8 xcmin,r_8 xcmax);
+  void ReCenterBin(void);
   //! Compute the correlation histogram
   void ToCorrel(void);