Changeset 307 in Sophya for trunk/SophyaLib/NTools/hisprof.cc

Timestamp:

May 19, 1999, 5:58:06 PM (26 years ago)

Author:

ansari

Message:

FIO_... + grosses modifs cmv 19/5/99

File:

• trunk/SophyaLib/NTools/hisprof.cc (modified) (3 diffs)

trunk/SophyaLib/NTools/hisprof.cc

@@ -7 +7 @@
 #include "cvector.h"
 
+using namespace PlanckDPC;
+
 //++
 // Class        HProf
@@ -47 +49 @@
   Histo::Errors();
   Zero();
-  END_CONSTRUCTOR
-}
-
-/********* Methode *********/
-//++
-HProf::HProf(char *flnm)
-//
-//      Constructeur a partir d'un fichier ppersist.
-//--
-: Histo()
-, SumY(NULL), SumY2(NULL), SumW(NULL), Ok(false), YMin(1.), YMax(-1.), Opt(0)
-{
-  PInPersist s(flnm);
-  Read(s);
   END_CONSTRUCTOR
 }
@@ -261 +249 @@
 //--
 {
-if (b.bins!=a.bins) THROW(sizeMismatchErr);
+if (b.NBins()!=a.NBins()) THROW(sizeMismatchErr);
 HProf c(a);
 return (c += b);
 }
 
-/********* Methode *********/
-//++
-void  HProf::WriteSelf(POutPersist& s)  const
-//
-//      Ecriture ppersist
-//--
+// Rappel des inlines functions pour commentaires
+//++
+//   inline Histo GetHisto()
+//      Retourne l'histogramme de profil.
+//--
+//++
+// inline void GetMean(Vector& v)
+//      Retourne le contenu de la moyenne dans le vecteur v
+//--
+//++
+// inline void GetError2(Vector& v)
+//      Retourne le contenu au carre de la dispersion/erreur dans le vecteur v
+//--
+//++
+// inline float operator()(int i) const
+//      Retourne le contenu du bin i
+//--
+//++
+// inline float Error2(int i) const
+//      Retourne le carre de la dispersion/erreur du bin i
+//--
+//++
+// inline float Error(int i) const
+//      Retourne la dispersion/erreur du bin i
+//--
+//++
+// inline int Fit(GeneralFit& gfit)
+//      Fit du profile par ``gfit''.
+//--
+//++
+// inline Histo* FitResidus(GeneralFit& gfit)
+//      Retourne l'Histogramme des residus par ``gfit''.
+//--
+//++
+// inline Histo* FitFunction(GeneralFit& gfit)
+//      Retourne l'Histogramme de la fonction fittee par ``gfit''.
+//--
+//++
+// inline void Print(int dyn,float hmin,float hmax,int pflag,int il,int ih)
+//      Print, voir detail dans Histo::Print
+//--
+
+
+///////////////////////////////////////////////////////////
+// --------------------------------------------------------
+//   Les objets delegues pour la gestion de persistance 
+// --------------------------------------------------------
+///////////////////////////////////////////////////////////
+
+FIO_HProf::FIO_HProf()
+{
+dobj=new HProf;
+ownobj=true;
+}
+
+FIO_HProf::FIO_HProf(string const & filename) 
+{
+dobj=new HProf;
+ownobj=true; 
+Read(filename);
+}
+
+FIO_HProf::FIO_HProf(const HProf & obj) 
+{ 
+dobj = new HProf(obj);
+ownobj=true; 
+}
+
+FIO_HProf::FIO_HProf(HProf * obj) 
+{ 
+dobj = obj;
+ownobj=false; 
+}
+
+FIO_HProf::~FIO_HProf()
+{
+if (ownobj && dobj) delete dobj;
+}
+
+AnyDataObj* FIO_HProf::DataObj()
+{
+return(dobj);
+}
+
+void FIO_HProf::ReadSelf(PInPersist& is)
 {
 char strg[256];
 
-UpdateHisto();
+if(dobj==NULL) dobj=new HProf;
+  else         dobj->Delete();
+
+// Lecture entete
+is.GetLine(strg,255);
+
+// Ecriture des valeurs
+is.Get(dobj->bins);
+is.Get(dobj->YMin);
+is.Get(dobj->YMax);
+is.Get(dobj->Opt);
+dobj->Ok = true;
+
+// Ecriture des donnees propres a l'histogramme de profil.
+is.GetLine(strg,255);
+dobj->SumY  = new double[dobj->bins];
+dobj->SumY2 = new double[dobj->bins];
+dobj->SumW  = new double[dobj->bins];
+is.Get(dobj->SumY,  dobj->bins);
+is.Get(dobj->SumY2, dobj->bins);
+is.Get(dobj->SumW,  dobj->bins);
+
+// Ecriture de l'histogramme
+FIO_Histo fio_h((Histo&)dobj);
+fio_h.Read(is);
+
+return;
+}
+
+void FIO_HProf::WriteSelf(POutPersist& os) const
+{
+if (dobj == NULL)   return;
+char strg[256];
+
+dobj->UpdateHisto();
 
 // Ecriture entete pour identifier facilement
-sprintf(strg,"HProf: YMin=%f  YMax=%f  Opt=%1d",YMin,YMax,Opt);
-s.PutLine(strg);
+sprintf(strg,"HProf: YMin=%f  YMax=%f  Opt=%1d",dobj->YMin,dobj->YMax,dobj->Opt);
+os.PutLine(strg);
 
 // Ecriture des valeurs
-s.PutI4(bins);
-s.PutR4(YMin);
-s.PutR4(YMax);
-s.PutU2(Opt);
+os.Put(dobj->bins);
+os.Put(dobj->YMin);
+os.Put(dobj->YMax);
+os.Put(dobj->Opt);
 
 // Ecriture des donnees propres a l'histogramme de profil.
 sprintf(strg,"HProf: SumY SumY2 SumW");
-s.PutLine(strg);
-s.PutR8s(SumY,  bins);
-s.PutR8s(SumY2, bins);
-s.PutR8s(SumW,  bins);
+os.PutLine(strg);
+os.Put(dobj->SumY,  dobj->bins);
+os.Put(dobj->SumY2, dobj->bins);
+os.Put(dobj->SumW,  dobj->bins);
 
 // Ecriture de l'histogramme
-Histo::WriteSelf(s);
+FIO_Histo fio_h((Histo&)dobj);
+fio_h.Write(os);
 
 return;
 }
-
-/********* Methode *********/
-//++
-void  HProf::ReadSelf(PInPersist& s)
-//
-//      Lecture ppersist
-//--
-{
-char strg[256];
-
-Delete();
-
-// Lecture entete
-s.GetLine(strg,255);
-
-// Ecriture des valeurs
-s.GetI4(bins);
-s.GetR4(YMin);
-s.GetR4(YMax);
-s.GetU2(Opt);
-Ok = true;
-
-// Ecriture des donnees propres a l'histogramme de profil.
-s.GetLine(strg,255);
-SumY  = new double[bins];
-SumY2 = new double[bins];
-SumW  = new double[bins];
-s.GetR8s(SumY,  bins);
-s.GetR8s(SumY2, bins);
-s.GetR8s(SumW,  bins);
-
-// Ecriture de l'histogramme
-Histo::ReadSelf(s);
-
-return;
-}
-
-
-// Rappel des inlines functions pour commentaires
-//++
-//   inline Histo GetHisto()
-//      Retourne l'histogramme de profil.
-//--
-//++
-// inline void GetMean(Vector& v)
-//      Retourne le contenu de la moyenne dans le vecteur v
-//--
-//++
-// inline void GetError2(Vector& v)
-//      Retourne le contenu au carre de la dispersion/erreur dans le vecteur v
-//--
-//++
-// inline float operator()(int i) const
-//      Retourne le contenu du bin i
-//--
-//++
-// inline float Error2(int i) const
-//      Retourne le carre de la dispersion/erreur du bin i
-//--
-//++
-// inline float Error(int i) const
-//      Retourne la dispersion/erreur du bin i
-//--
-//++
-// inline int Fit(GeneralFit& gfit)
-//      Fit du profile par ``gfit''.
-//--
-//++
-// inline Histo* FitResidus(GeneralFit& gfit)
-//      Retourne l'Histogramme des residus par ``gfit''.
-//--
-//++
-// inline Histo* FitFunction(GeneralFit& gfit)
-//      Retourne l'Histogramme de la fonction fittee par ``gfit''.
-//--
-//++
-// inline void Print(int dyn,float hmin,float hmax,int pflag,int il,int ih)
-//      Print, voir detail dans Histo::Print
-//--