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

Timestamp:

Jul 11, 2000, 8:57:25 PM (25 years ago)

Author:

ansari

Message:

pour PutValues,Error etc ... on prend maintenant le mini

du nbre d'elements du vecteur/matrice et de l'histo.

cmv 11/07/00

File:

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

trunk/SophyaLib/HiStats/histos.cc

@@ -1 +1 @@
 //
-// $Id: histos.cc,v 1.8 2000-07-07 09:44:15 ansari Exp $
+// $Id: histos.cc,v 1.9 2000-07-11 18:57:25 ansari Exp $
 //
 
@@ -326 +326 @@
 void Histo::PutValue(TVector<r_8> &v, int ierr)
 {
-if(v.NElts()<(uint_4) bins) THROW(sizeMismatchErr);
-for(int i=0;i<bins;i++) {
+//if(v.NElts()<(uint_4) bins) THROW(sizeMismatchErr);
+uint_4 n = (v.NElts()<(uint_4) bins) ? v.NElts(): (uint_4) bins;
+if(n>0) for(int i=0;i<n;i++) {
   data[i] = v(i);
   if(err2&&ierr) err2[i] = fabs(v(i));
@@ -339 +340 @@
 void Histo::PutValueAdd(TVector<r_8> &v, int ierr)
 {
-if(v.NElts()<(uint_4) bins) THROW(sizeMismatchErr);
-for(int i=0;i<bins;i++) {
+//if(v.NElts()<(uint_4) bins) THROW(sizeMismatchErr);
+uint_4 n = (v.NElts()<(uint_4) bins) ? v.NElts(): (uint_4) bins;
+if(n>0) for(int i=0;i<n;i++) {
   data[i] += v(i);
   if(err2 && ierr) err2[i] += fabs(v(i));
@@ -352 +354 @@
 void Histo::PutError2(TVector<r_8> &v)
 {
-if(v.NElts()<(uint_4) bins) THROW(sizeMismatchErr);
-if(!err2) Errors();
-for(int i=0;i<bins;i++) err2[i] = v(i);
+//if(v.NElts()<(uint_4) bins) THROW(sizeMismatchErr);
+uint_4 n = (v.NElts()<(uint_4) bins) ? v.NElts(): (uint_4) bins;
+if(n>0) {
+  if(!err2) Errors();
+  for(int i=0;i<n;i++) err2[i] = v(i);
+}
 return;
 }
@@ -363 +368 @@
 void Histo::PutError2Add(TVector<r_8> &v)
 {
-if(v.NElts()<(uint_4) bins) THROW(sizeMismatchErr);
-if(!err2) Errors();
-for(int i=0;i<bins;i++) if(v(i)>0.) err2[i] += v(i);
+//if(v.NElts()<(uint_4) bins) THROW(sizeMismatchErr);
+uint_4 n = (v.NElts()<(uint_4) bins) ? v.NElts(): (uint_4) bins;
+if(n>0) {
+  if(!err2) Errors();
+  for(int i=0;i<n;i++) if(v(i)>0.) err2[i] += v(i);
+}
 return;
 }
@@ -374 +382 @@
 void Histo::PutError(TVector<r_8> &v)
 {
-if(v.NElts()<(uint_4) bins) THROW(sizeMismatchErr);
-if(!err2) Errors();
-for(int i=0;i<bins;i++)
-  if(v(i)>0.) err2[i]=v(i)*v(i); else err2[i]=-v(i)*v(i);
+//if(v.NElts()<(uint_4) bins) THROW(sizeMismatchErr);
+uint_4 n = (v.NElts()<(uint_4) bins) ? v.NElts(): (uint_4) bins;
+if(n>0) {
+  if(!err2) Errors();
+  for(int i=0;i<n;i++)
+    if(v(i)>0.) err2[i]=v(i)*v(i); else err2[i]=-v(i)*v(i);
+}
 return;
 }