Changeset 3053 in Sophya for trunk/SophyaLib/HiStats/histos2.cc

Timestamp:

Aug 12, 2006, 7:29:29 PM (19 years ago)

Author:

cmv

Message:

correction bug + addaptation pour ecriture fits cmv 12/8/2006

File:

• trunk/SophyaLib/HiStats/histos2.cc (modified) (7 diffs)

trunk/SophyaLib/HiStats/histos2.cc

@@ -39 +39 @@
 
 ///////////////////////////////////////////////////////////////////
+
+/*!
+  Constructeur par defaut.
+*/
+Histo2D::Histo2D()
+: mData(NULL), mErr2(NULL)
+, nHist(0), nEntries(0)
+, mNx(0), mNy(0), mNxy(0), mXmin(0), mXmax(0), mYmin(0), mYmax(0), mWBinx(0), mWBiny(0)
+, mHprojx(NULL), mHprojy(NULL)
+{
+ for(int_4 i=0;i<3;i++) for(int_4 j=0;j<3;j++) mOver[i][j]=0.;
+ mB_s.H = NULL;
+}
+
 /*!
   Createur d'un histogramme 2D ayant nxBin,nyBin bins
@@ -44 +58 @@
 */
 Histo2D::Histo2D(r_8 xMin,r_8 xMax,int_4 nxBin,r_8 yMin,r_8 yMax,int_4 nyBin)
-      : mData(new r_8[nxBin*nyBin]), mErr2(NULL)
-      , nHist(0), nEntries(0)
-      , mNx(nxBin), mNy(nyBin), mNxy(nxBin*nyBin)
-      , mXmin(xMin), mXmax(xMax), mYmin(yMin), mYmax(yMax)
-      , mWBinx((xMax - xMin)/nxBin), mWBiny((yMax - yMin)/nyBin)
-      , mHprojx(NULL), mHprojy(NULL)
-{
-ASSERT(nxBin>0 && nyBin>0 && xMin<xMax && yMin<yMax);
-for(int_4 i=0;i<3;i++) for(int_4 j=0;j<3;j++) mOver[i][j]=0.;
-Zero();
-mB_s.H = NULL;
+: mData(NULL), mErr2(NULL), mHprojx(NULL), mHprojy(NULL)
+{
+  CreateOrResize(xMin,xMax,nxBin,yMin,yMax,nyBin);
 }
 
@@ -62 +68 @@
 */
 Histo2D::Histo2D(r_4 xMin,r_4 xMax,int_4 nxBin,r_4 yMin,r_4 yMax,int_4 nyBin)
-      : mData(new r_8[nxBin*nyBin]), mErr2(NULL)
-      , nHist(0), nEntries(0)
-      , mNx(nxBin), mNy(nyBin), mNxy(nxBin*nyBin)
-      , mXmin((r_8)xMin), mXmax((r_8)xMax), mYmin((r_8)yMin), mYmax((r_8)yMax)
-      , mWBinx((xMax - xMin)/nxBin), mWBiny((yMax - yMin)/nyBin)
-      , mHprojx(NULL), mHprojy(NULL)
-{
-ASSERT(nxBin>0 && nyBin>0 && xMin<xMax && yMin<yMax);
-for(int_4 i=0;i<3;i++) for(int_4 j=0;j<3;j++) mOver[i][j]=0.;
-Zero();
-mB_s.H = NULL;
+: mData(NULL), mErr2(NULL), mHprojx(NULL), mHprojy(NULL)
+{
+  CreateOrResize((r_8)xMin,(r_8)xMax,nxBin,(r_8)yMin,(r_8)yMax,nyBin);
 }
 
@@ -142 +140 @@
 
 /*!
-  Constructeur par defaut.
-*/
-Histo2D::Histo2D()
-      : mData(NULL), mErr2(NULL)
-      , nHist(0), nEntries(0)
-      , mNx(0), mNy(0), mNxy(0)
-      , mXmin(0), mXmax(0), mYmin(0), mYmax(0)
-      , mWBinx(0), mWBiny(0)
-      , mHprojx(NULL), mHprojy(NULL)
-{
-for(int_4 i=0;i<3;i++) for(int_4 j=0;j<3;j++) mOver[i][j]=0.;
-mB_s.H = NULL;
-}
-
-///////////////////////////////////////////////////////////////////
+  Destructeur.
+*/
+Histo2D::~Histo2D()
+{
+  Delete();
+}
+
+///////////////////////////////////////////////////////////////////
+/*!
+  allocation de la place de l'histogramme (fct privee).
+*/
+void Histo2D::CreateOrResize(r_8 xMin,r_8 xMax,int_4 nxBin,r_8 yMin,r_8 yMax,int_4 nyBin)
+{
+ int_8 n = nxBin*nyBin;
+ if(mData!=NULL) {delete[] mData; mData = NULL;}
+ if(n>0) {mData = new r_8[n]; memset(mData, 0, n*sizeof(r_8));}
+ mNx = nxBin; mNy = nyBin; mNxy = nxBin*nyBin;
+ mXmin = xMin; mXmax = xMax; mYmin = yMin; mYmax = yMax;
+ mWBinx = (nxBin>0) ? (xMax - xMin)/nxBin: 0.;
+ mWBiny = (nyBin>0) ? (yMax - yMin)/nyBin: 0.;
+ nHist = 0; nEntries = 0;
+ for(int_4 i=0;i<3;i++) for(int_4 j=0;j<3;j++) mOver[i][j]=0.;
+ mB_s.H = NULL;
+}
+
 /*!
   Desallocation de la place de l'histogramme (fct privee).
@@ -163 +171 @@
 {
   if( mData != NULL ) { delete[] mData; mData = NULL;}
-
   if( mErr2 != NULL ) { delete[] mErr2; mErr2 = NULL;}
 
@@ -180 +187 @@
   mWBinx = 0; mWBiny = 0;
   for(int_4 i=0;i<3;i++) for(int_4 j=0;j<3;j++) mOver[i][j]=0.;
-  mB_s.H = NULL;
-}
-
-/*!
-  Destructeur.
-*/
-Histo2D::~Histo2D()
-{
-Delete();
+  mB_s.H = NULL;  // Ne pas faire de delete!
 }
 
@@ -219 +218 @@
 void Histo2D::Errors()
 {
- if( mNxy > 0 ) {
-   if(mErr2==NULL) mErr2 = new r_8[mNxy];
-   memset(mErr2, 0, mNxy*sizeof(r_8));
- }
+ if(mErr2!=NULL) {delete[] mErr2; mErr2 = NULL;}
+ if(mNxy <= 0) return;
+ mErr2 = new r_8[mNxy];
+ memset(mErr2,0,mNxy*sizeof(r_8));
 }