Changeset 1056 in Sophya for trunk/SophyaLib/HiStats/histos.cc

Timestamp:

Jul 7, 2000, 10:41:11 AM (25 years ago)

Author:

ansari

Message:

cmv 7/7/2000

hisprof.cc,h:

fonction IsOk()
GetMean -> GetValue (mauvais nom)
float Error2() -> double Error2()
nouveau: GetError(TVector<r_8>& v)
HProf::PrintF() sur-ecriture de Histo::PrintF
protection dans createur par copie dans alloc

SumY... pour le cas ou H.bins==0

protection dans UpdateHisto pour HProf cree par defaut (bins==0)
Dans WriteSelf UpdateHisto appele que si necessaire (IsOk()?)

histos2.cc : le Print() dit si les erreurs ont ete demandees
histos.cc,h:

protection dans createur par copie dans alloc

pour le cas ou H.bins==0

protection dans Zero() pour Histo cree par defaut (bins==0)
protection dans operator=() pour Histo cree par defaut (bins==0)
le Print() dit si les erreurs ont ete demandees

cmv 7/7/2000

File:

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

trunk/SophyaLib/HiStats/histos.cc

@@ -1 +1 @@
 //
-// $Id: histos.cc,v 1.6 2000-06-30 13:21:52 ansari Exp $
+// $Id: histos.cc,v 1.7 2000-07-07 08:41:11 ansari Exp $
 //
 
@@ -45 +45 @@
 /*! Constructeur par copie */
 Histo::Histo(const Histo& H)
-: data(new float[H.bins]), err2(NULL),
+: data((H.bins>0)? new float[H.bins] : NULL),
+  err2(NULL),
   under(H.under), over(H.over), nHist(H.nHist), nEntries(H.nEntries),
   bins(H.bins), min(H.min), max(H.max),
   binWidth(H.binWidth)
 {
-  memcpy(data, H.data, bins*sizeof(float));
-  if( H.err2 != NULL ) {
-    err2 = new double[bins];
-    memcpy(err2, H.err2, bins*sizeof(double));
+  if(bins>0) {
+    memcpy(data, H.data, bins*sizeof(float));
+    if( H.err2 != NULL ) {
+      err2 = new double[bins];
+      memcpy(err2, H.err2, bins*sizeof(double));
+    }
   }
   END_CONSTRUCTOR
@@ -82 +85 @@
 void Histo::Zero()
 {
+  if(bins<=0) return;
   memset(data, 0, bins*sizeof(float));
   under = over = 0;
@@ -95 +99 @@
 void Histo::Errors()
 {
- if( bins > 0 ) {
-   if(err2==NULL) err2 = new double[bins];
-   memset(err2, 0, bins*sizeof(double));
- }
-}
-
-/********* Methode *********/
-/*!
-  Operateur egal
+ if(bins<=0) return;
+ if(err2==NULL) err2 = new double[bins];
+ memset(err2, 0, bins*sizeof(double));
+}
+
+/********* Methode *********/
+/*!
+  Operateur egal Histo = Histo
 */
 Histo& Histo::operator = (const Histo& h)
 {
   if(this == &h) return *this;
+  if( h.bins <= 0 ) {Delete(); return *this;}
   if( h.bins > bins ) Delete();
+  if(!h.err2 && err2 ) { delete [] err2; err2=NULL;}
   if(!data) data = new float[h.bins];
-  if( !h.err2 && err2 ) { delete [] err2; err2=NULL;}
-  if( h.err2 && !err2 ) err2 = new double[h.bins];
+  if(h.err2 && !err2 ) err2 = new double[h.bins];
 
   under = h.under;
@@ -817 +821 @@
   if( nbinew <= 0 ) return;
 
-  // memoraisation de l'histogramme original
+  // memorisation de l'histogramme original
   Histo H(*this);
 
@@ -1324 +1328 @@
               int il, int ih)
 {
-
  
  if( il > ih ) { il = 0; ih = bins-1; }
@@ -1342 +1345 @@
       << "  min=" << min << "  max=" << max
       << "  binWidth=" << binWidth << endl;
- cout << "       mean=" << Mean() << " r.m.s=" << Sigma() << endl;
+ cout << "       mean=" << Mean() << " r.m.s=" << Sigma() 
+      << "  Errors="<<HasErrors()<< endl;
 
  if(hdyn<0 || pflag<0 ) return;