Changeset 1056 in Sophya for trunk/SophyaLib/HiStats/hisprof.h

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/hisprof.h (modified) (3 diffs)

trunk/SophyaLib/HiStats/hisprof.h

@@ -22 +22 @@
 
   // UPDATING or SETTING
-  void   Zero();
-  void   UpdateHisto() const;
-  void   SetErrOpt(bool spread = true);
-  void   Add(float x, float y, float w = 1.);
-  void   AddBin(int numBin, float y, float w = 1.);
+  void UpdateHisto() const;
+  void SetErrOpt(bool spread = true);
+  void Zero();
+  void Add(float x, float y, float w = 1.);
+  void AddBin(int numBin, float y, float w = 1.);
   
   // Acces a l information
+  //! Retourne true si l'histogramme est a jours, false sinon.
+  inline bool IsOk() const {return Ok;}
   //! Retourne l'histogramme de profil.
   inline Histo GetHisto()
                {if(!Ok) UpdateHisto(); return (Histo) *this;}
   //! Retourne le contenu de la moyenne dans le vecteur v
-  inline void GetMean(TVector<r_8>& v)
+  inline void GetValue(TVector<r_8>& v)
               {if(!Ok) UpdateHisto(); Histo::GetValue(v);}
   //! Retourne le contenu au carre de la dispersion/erreur dans le vecteur v
   inline void GetError2(TVector<r_8>& v)
               {if(!Ok) UpdateHisto(); Histo::GetError2(v);}
+  //! Retourne le contenu au carre de la dispersion/erreur dans le vecteur v
+  inline void GetError(TVector<r_8>& v)
+              {if(!Ok) UpdateHisto(); Histo::GetError(v);}
   //! Retourne le contenu du bin i
   inline float operator()(int i) const
                {if(!Ok) UpdateHisto(); return data[i];}
   //! Retourne le carre de la dispersion/erreur du bin i
-  inline float Error2(int i) const
-               {if(!Ok) UpdateHisto(); return (float) err2[i];}
+  inline double Error2(int i) const
+         {if(!Ok) UpdateHisto(); return (float) err2[i];}
   //! Retourne la dispersion/erreur du bin i
   inline float Error(int i) const
          {if(!Ok) UpdateHisto();
-          return err2[i]>0. ? (float) sqrt(err2[i]) : 0.f;}
+          return (err2[i]>0.) ? (float) sqrt(err2[i]) : 0.f;}
 
   // Operators
@@ -55 +60 @@
   // Fit
   //! Fit du profile par ``gfit''.
-  inline int    Fit(GeneralFit& gfit)
+  inline int Fit(GeneralFit& gfit)
          {if(!Ok) UpdateHisto(); return Histo::Fit(gfit,0);}
   //! Retourne l'Histogramme des residus par ``gfit''.
-  inline Histo  FitResidus(GeneralFit& gfit)
+  inline Histo FitResidus(GeneralFit& gfit)
          {if(!Ok) UpdateHisto(); return Histo::FitResidus(gfit);}
   //! Retourne l'Histogramme de la fonction fittee par ``gfit''.
-  inline Histo  FitFunction(GeneralFit& gfit)
+  inline Histo FitFunction(GeneralFit& gfit)
          {if(!Ok) UpdateHisto(); return Histo::FitFunction(gfit);}
 
   // Print
   //! Print, voir detail dans Histo::Print
-  inline void Print(int dyn=100,float hmin=1.,float hmax=-1.,int pflag=0,int il=1,int ih=-1)
-              {if(!Ok) UpdateHisto(); Histo::Print(dyn,hmin,hmax,pflag,il,ih);}
+  inline void Print(int dyn=100,float hmin=1.,float hmax=-1.
+                   ,int pflag=0,int il=1,int ih=-1)
+         {if(!Ok) UpdateHisto(); Histo::Print(dyn,hmin,hmax,pflag,il,ih);}
+  //! PrintF, voir detail dans Histo::PrintF
+  inline void PrintF(FILE * fp,int dyn=100,float hmin=1.,float hmax=-1.
+                    ,int pflag=0,int il=1,int ih=-1)
+         {if(!Ok) UpdateHisto(); Histo::PrintF(fp,dyn,hmin,hmax,pflag,il,ih);}
 
 protected:
@@ -75 +85 @@
   double*        SumY2; //!< somme des carres
   double*        SumW;  //!< somme des poids
-  bool           Ok;    //!< true isiupdate fait
+  bool           Ok;    //!< true if update fait
   float          YMin;  //!< limite minimum Y pour somme
   float          YMax;  //!< limite maximum Y pour somme