Changeset 3184 in Sophya for trunk/SophyaLib/HiStats/hist2err.cc

Timestamp:

Feb 13, 2007, 3:28:00 PM (19 years ago)

Author:

cmv

Message:

write Histo2DErr into ASCII file cmv 13/02/2007

File:

• trunk/SophyaLib/HiStats/hist2err.cc (modified) (1 diff)

trunk/SophyaLib/HiStats/hist2err.cc

@@ -261 +261 @@
 }
 
+/********* Methode *********/
+/*!
+  Write to an ASCII file
+*/
+int Histo2DErr::WriteASCII(string fname)
+{
+  FILE *file = fopen(fname.c_str(),"w");
+  if(file==NULL) {
+    cout<<"Histo2DErr::WriteASCII_Error: error opening "<<fname<<endl;
+    return -1;
+  }
+
+  if(NBinX()<=0 || NBinY()<=0) {
+    cout<<"Histo2DErr::WriteASCII_Error: wrong number of bins"<<endl;
+    return -2;
+  }
+
+  fprintf(file,"%ld %.17e %.17e %.17e %ld %.17e %.17e %.17e %d\n"
+         ,(long)NBinX(),XMin(),XMax(),WBinX()
+         ,(long)NBinY(),YMin(),YMax(),WBinY()
+         ,NMean());
+  for(long i=0;i<NBinX();i++) for(long j=0;j<NBinY();j++) {
+    // ligne = i*NY+j
+    fprintf(file,"%d %d %.17e %.17e %.0f\n"
+           ,i,j,(*this)(i,j),Error2(i,j),NEntBin(i,j));
+  }
+
+  fclose(file);
+  return 0;
+}
+
+/*!
+  Read from an ASCII file
+*/
+#define __LENLINE_Histo2DErr_ReadASCII__ 2048
+int Histo2DErr::ReadASCII(string fname)
+{
+  FILE *file = fopen(fname.c_str(),"r");
+  if(file==NULL) {
+    cout<<"Histo2DErr::ReadASCII_Error: error opening "<<fname<<endl;
+    return -1;
+  }
+
+  char line[__LENLINE_Histo2DErr_ReadASCII__];
+  long n=0, nbinx=0, nbiny=0;
+
+  while ( fgets(line,__LENLINE_Histo2DErr_ReadASCII__,file) != NULL ) {
+
+    if(n==0) {
+
+      r_8 xmin,xmax,wx, ymin,ymax,wy; long mnmean=1;
+      sscanf(line,"%d %lf %lf %lf %d %lf %lf %lf %d"
+            ,&nbinx,&xmin,&xmax,&wx
+            ,&nbiny,&ymin,&ymax,&wy
+            ,&mnmean);
+      if(nbinx<=0 || nbiny<=0) {
+        cout<<"Histo2Err::ReadASCII_Error: wrong number of bins"<<endl;
+        return -2;
+      }
+      CreateOrResize(xmin,xmax,nbinx,ymin,ymax,nbiny);
+      SetMean(mnmean);
+
+    } else {
+
+      long i,j; r_8 v,e2,nb;
+      sscanf(line,"%d %d %lf %lf %lf",&i,&j,&v,&e2,&nb);
+      SetBin(i,j,v);
+      SetErr2(i,j,e2);
+      SetNentB(i,j,nb);
+
+    }
+
+    n++;
+  }
+
+  fclose(file);
+  return 0;
+}
+
 ///////////////////////////////////////////////////////////
 // --------------------------------------------------------