Changeset 3156 in Sophya for trunk/SophyaLib/HiStats

Timestamp:

Jan 22, 2007, 6:23:55 PM (19 years ago)

Author:

cmv

Message:

add write and read in/from ASCII file cmv 22/01/2007

Location:

trunk/SophyaLib/HiStats

Files:

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

trunk/SophyaLib/HiStats/histerr.cc

@@ -226 +226 @@
 void HistoErr::Show(ostream & os) const
 {
-  os <<"HistoErr(nmean="<<mMean<<")"<<endl
-     <<"          nx="<<nx_<<" ["<<xmin_<<","<<xmax_<<"] dx="<<dx_<<endl;
-}
+  os <<"HistoErr(nmean="<<mMean<<")"
+     <<" nx="<<nx_<<" ["<<xmin_<<","<<xmax_<<"] dx="<<dx_<<endl;
+}
+
+/********* Methode *********/
+/*!
+  Write to an ASCII file
+*/
+int HistoErr::WriteASCII(string fname)
+{
+  FILE *file = fopen(fname.c_str(),"w");
+  if(file==NULL) {
+    cout<<"HistoErr::WriteASCII_Error: error opening "<<fname<<endl;
+    return -1;
+  }
+
+  if(NBins()<=0) {
+    cout<<"HistoErr::WriteASCII_Error: nbin= "<<NBins()<<endl;
+    return -2;
+  }
+
+  fprintf(file,"%ld %.17e %.17e %d\n",(long)NBins(),XMin(),XMax(),NMean());
+  for(long i=0;i<NBins();i++) {
+      fprintf(file,"%d %.17e %.17e %.17e %.0f\n"
+             ,i,BinCenter(i),(*this)(i),Error2(i),NEntBin(i));
+  }
+
+  fclose(file);
+  return 0;
+}
+
+/*!
+  Read from an ASCII file
+*/
+#define __LENLINE_HistoErr_ReadASCII__ 2048
+int HistoErr::ReadASCII(string fname)
+{
+  FILE *file = fopen(fname.c_str(),"r");
+  if(file==NULL) {
+    cout<<"HistoErr::ReadASCII_Error: error opening "<<fname<<endl;
+    return -1;
+  }
+
+  char line[__LENLINE_HistoErr_ReadASCII__];
+  long n=0, nbin=0;
+
+  while ( fgets(line,__LENLINE_HistoErr_ReadASCII__,file) != NULL ) {
+
+    if(n==0) {
+
+      r_8 xmin,xmax; long mnmean=1;
+      sscanf(line,"%d %lf %lf %d",&nbin,&xmin,&xmax,&mnmean);
+      if(nbin<=0) {
+        cout<<"HistoErr::ReadASCII_Error: nbin= "<<nbin<<endl;
+        return -2;
+      }
+      CreateOrResize(xmin,xmax,nbin);
+      SetMean(mnmean);
+
+    } else {
+
+      long i; r_8 x,v,e2,nb;
+      sscanf(line,"%d %lf %lf %lf %lf",&i,&x,&v,&e2,&nb);
+      SetBin(i,v);
+      SetErr2(i,e2);
+      SetNentB(i,nb);
+
+    }
+
+    n++;
+  }
+
+  fclose(file);
+  return 0;
+}
+
 
 ///////////////////////////////////////////////////////////

trunk/SophyaLib/HiStats/histerr.h

@@ -144 +144 @@
   void Show() const { Show(cout); }
 
+  // Write ASCII
+  int WriteASCII(string fname);
+  int ReadASCII(string fname);
+
 protected:
   void CreateOrResize(r_8 xmin,r_8 xmax,int_4 nx);