Changeset 943 in Sophya for trunk/SophyaLib


Ignore:
Timestamp:
Apr 15, 2000, 1:06:28 PM (25 years ago)
Author:
ansari
Message:

Vector,Matrix -> TVector,TMatrix<r_8> cmv 15/4/00

Location:
trunk/SophyaLib/HiStats
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • trunk/SophyaLib/HiStats/hisprof.h

    r926 r943  
    3030  // Acces a l information
    3131  //! Retourne l'histogramme de profil.
    32   inline Histo GetHisto() {if(!Ok) UpdateHisto(); return (Histo) *this;}
     32  inline Histo GetHisto()
     33               {if(!Ok) UpdateHisto(); return (Histo) *this;}
    3334  //! Retourne le contenu de la moyenne dans le vecteur v
    34   inline void GetMean(Vector& v) {if(!Ok) UpdateHisto(); Histo::GetValue(v);}
     35  inline void GetMean(TVector<r_8>& v)
     36              {if(!Ok) UpdateHisto(); Histo::GetValue(v);}
    3537  //! Retourne le contenu au carre de la dispersion/erreur dans le vecteur v
    36   inline void GetError2(Vector& v) {if(!Ok) UpdateHisto(); Histo::GetError2(v);}
     38  inline void GetError2(TVector<r_8>& v)
     39              {if(!Ok) UpdateHisto(); Histo::GetError2(v);}
    3740  //! Retourne le contenu du bin i
    38   inline float operator()(int i) const {if(!Ok) UpdateHisto(); return data[i];}
     41  inline float operator()(int i) const
     42               {if(!Ok) UpdateHisto(); return data[i];}
    3943  //! Retourne le carre de la dispersion/erreur du bin i
    40   inline float Error2(int i) const {if(!Ok) UpdateHisto(); return (float) err2[i];}
     44  inline float Error2(int i) const
     45               {if(!Ok) UpdateHisto(); return (float) err2[i];}
    4146  //! Retourne la dispersion/erreur du bin i
    4247  inline float Error(int i) const
    43          {if(!Ok) UpdateHisto(); return err2[i]>0. ? (float) sqrt(err2[i]) : 0.f;}
     48         {if(!Ok) UpdateHisto();
     49          return err2[i]>0. ? (float) sqrt(err2[i]) : 0.f;}
    4450
    4551  // Operators
     
    6268  //! Print, voir detail dans Histo::Print
    6369  inline void Print(int dyn=100,float hmin=1.,float hmax=-1.,int pflag=0,int il=1,int ih=-1)
    64                            {if(!Ok) UpdateHisto(); Histo::Print(dyn,hmin,hmax,pflag,il,ih);}
     70              {if(!Ok) UpdateHisto(); Histo::Print(dyn,hmin,hmax,pflag,il,ih);}
    6571
    6672protected:
  • trunk/SophyaLib/HiStats/histos.cc

    r926 r943  
    11//
    2 // $Id: histos.cc,v 1.3 2000-04-13 18:38:32 ansari Exp $
     2// $Id: histos.cc,v 1.4 2000-04-15 11:06:27 ansari Exp $
    33//
    44
     
    387387  Remplissage d'un tableau avec la valeur des abscisses
    388388*/
    389 void Histo::GetAbsc(Vector &v)
     389void Histo::GetAbsc(TVector<r_8> &v)
    390390{
    391391v.Realloc(bins);
     
    397397  Remplissage d'un tableau avec la valeur du contenu
    398398*/
    399 void Histo::GetValue(Vector &v)
     399void Histo::GetValue(TVector<r_8> &v)
    400400{
    401401v.Realloc(bins);
     
    407407  Remplissage d'un tableau avec la valeur des erreurs au carre
    408408*/
    409 void Histo::GetError2(Vector &v)
     409void Histo::GetError2(TVector<r_8> &v)
    410410{
    411411v.Realloc(bins);
     
    418418  Remplissage d'un tableau avec la valeur des erreurs
    419419*/
    420 void Histo::GetError(Vector &v)
     420void Histo::GetError(TVector<r_8> &v)
    421421{
    422422v.Realloc(bins);
     
    430430  Remplissage du contenu de l'histo avec les valeurs d'un vecteur
    431431*/
    432 void Histo::PutValue(Vector &v, int ierr)
     432void Histo::PutValue(TVector<r_8> &v, int ierr)
    433433{
    434434if(v.NElts()<bins) THROW(sizeMismatchErr);
     
    443443  Addition du contenu de l'histo avec les valeurs d'un vecteur
    444444*/
    445 void Histo::PutValueAdd(Vector &v, int ierr)
     445void Histo::PutValueAdd(TVector<r_8> &v, int ierr)
    446446{
    447447if(v.NElts()<bins) THROW(sizeMismatchErr);
     
    456456  Remplissage des erreurs au carre de l'histo avec les valeurs d'un vecteur
    457457*/
    458 void Histo::PutError2(Vector &v)
     458void Histo::PutError2(TVector<r_8> &v)
    459459{
    460460if(v.NElts()<bins) THROW(sizeMismatchErr);
     
    467467  Addition des erreurs au carre de l'histo avec les valeurs d'un vecteur
    468468*/
    469 void Histo::PutError2Add(Vector &v)
     469void Histo::PutError2Add(TVector<r_8> &v)
    470470{
    471471if(v.NElts()<bins) THROW(sizeMismatchErr);
     
    478478  Remplissage des erreurs de l'histo avec les valeurs d'un vecteur
    479479*/
    480 void Histo::PutError(Vector &v)
     480void Histo::PutError(TVector<r_8> &v)
    481481{
    482482if(v.NElts()<bins) THROW(sizeMismatchErr);
     
    10651065  }
    10661066
    1067   Vector xFit(nLowHigh);
    1068   Vector yFit(nLowHigh);
    1069   Vector e2Fit(nLowHigh);
    1070   Vector errcoef(1);
     1067  TVector<r_8> xFit(nLowHigh);
     1068  TVector<r_8> yFit(nLowHigh);
     1069  TVector<r_8> e2Fit(nLowHigh);
     1070  TVector<r_8> errcoef(1);
    10711071  int ii = 0;
    10721072  for (int j=iLow; j<=iHigh; j++) {
     
    13851385if(f==NULL)
    13861386  throw(NullPtrError("Histo::FitResidus: NULL pointer\n"));
    1387 Vector par = gfit.GetParm();
     1387TVector<r_8> par = gfit.GetParm();
    13881388Histo h(*this);
    13891389for(int i=0;i<NBins();i++) {
     
    14041404if(f==NULL)
    14051405  throw(NullPtrError("Histo::FitFunction: NULL pointer\n"));
    1406 Vector par = gfit.GetParm();
     1406TVector<r_8> par = gfit.GetParm();
    14071407Histo h(*this);
    14081408for(int i=0;i<NBins();i++) {
  • trunk/SophyaLib/HiStats/histos.h

    r926 r943  
    11// This may look like C code, but it is really -*- C++ -*-
    22//
    3 // $Id: histos.h,v 1.4 2000-04-13 18:38:33 ansari Exp $
     3// $Id: histos.h,v 1.5 2000-04-15 11:06:28 ansari Exp $
    44//
    55
     
    6767
    6868  // get/put dans/depuis un vector
    69   void GetAbsc(Vector& v);
    70   void GetValue(Vector& v);
    71   void GetError2(Vector& v);
    72   void GetError(Vector& v);
    73   void PutValue(Vector& v, int ierr=0);
    74   void PutValueAdd(Vector &v, int ierr=0);
    75   void PutError2(Vector& v);
    76   void PutError2Add(Vector& v);
    77   void PutError(Vector& v);
     69  void GetAbsc(TVector<r_8>& v);
     70  void GetValue(TVector<r_8>& v);
     71  void GetError2(TVector<r_8>& v);
     72  void GetError(TVector<r_8>& v);
     73  void PutValue(TVector<r_8>& v, int ierr=0);
     74  void PutValueAdd(TVector<r_8> &v, int ierr=0);
     75  void PutError2(TVector<r_8>& v);
     76  void PutError2Add(TVector<r_8>& v);
     77  void PutError(TVector<r_8>& v);
    7878
    7979  // INLINES
  • trunk/SophyaLib/HiStats/histos2.cc

    r926 r943  
    603603  Remplissage d'un tableau avec les valeurs des abscisses.
    604604*/
    605 void Histo2D::GetXCoor(Vector &v)
     605void Histo2D::GetXCoor(TVector<r_8> &v)
    606606{
    607607float x,y;
     
    614614  Remplissage d'un tableau avec les valeurs des ordonnees.
    615615*/
    616 void Histo2D::GetYCoor(Vector &v)
     616void Histo2D::GetYCoor(TVector<r_8> &v)
    617617{
    618618float x,y;
     
    625625  Remplissage d'un tableau avec les valeurs du contenu.
    626626*/
    627 void Histo2D::GetValue(Matrix &v)
     627void Histo2D::GetValue(TMatrix<r_8> &v)
    628628{
    629629v.Realloc(nx,ny);
     
    636636  Remplissage d'un tableau avec les valeurs du carre des erreurs.
    637637*/
    638 void Histo2D::GetError2(Matrix &v)
     638void Histo2D::GetError2(TMatrix<r_8> &v)
    639639{
    640640int i,j;
     
    649649  Remplissage d'un tableau avec les valeurs des erreurs.
    650650*/
    651 void Histo2D::GetError(Matrix &v)
     651void Histo2D::GetError(TMatrix<r_8> &v)
    652652{
    653653int i,j;
     
    663663  Remplissage du contenu de l'histo avec les valeurs d'un tableau.
    664664*/
    665 void Histo2D::PutValue(Matrix &v, int ierr)
     665void Histo2D::PutValue(TMatrix<r_8> &v, int ierr)
    666666{
    667667int i,j;
     
    677677  Addition du contenu de l'histo avec les valeurs d'un tableau.
    678678*/
    679 void Histo2D::PutValueAdd(Matrix &v, int ierr)
     679void Histo2D::PutValueAdd(TMatrix<r_8> &v, int ierr)
    680680{
    681681int i,j;
     
    692692  avec les valeurs d'un tableau.
    693693*/
    694 void Histo2D::PutError2(Matrix &v)
     694void Histo2D::PutError2(TMatrix<r_8> &v)
    695695{
    696696int i,j;
     
    705705  avec les valeurs d'un tableau.
    706706*/
    707 void Histo2D::PutError2Add(Matrix &v)
     707void Histo2D::PutError2Add(TMatrix<r_8> &v)
    708708{
    709709int i,j;
     
    718718  Remplissage des erreurs de l'histo avec les valeurs d'un tableau.
    719719*/
    720 void Histo2D::PutError(Matrix &v)
     720void Histo2D::PutError(TMatrix<r_8> &v)
    721721{
    722722int i,j;
     
    10711071if(f==NULL)
    10721072  throw(NullPtrError("Histo2D::FitResidus: NULL pointer\n"));
    1073 Vector par = gfit.GetParm();
     1073TVector<r_8> par = gfit.GetParm();
    10741074Histo2D h2(*this);
    10751075for(int i=0;i<NBinX();i++) for(int j=0;j<NBinY();j++) {
     
    10921092if(f==NULL)
    10931093  throw(NullPtrError("Histo2D::FitFunction: NULL pointer\n"));
    1094 Vector par = gfit.GetParm();
     1094TVector<r_8> par = gfit.GetParm();
    10951095Histo2D h2(*this);
    10961096for(int i=0;i<NBinX();i++) for(int j=0;j<NBinY();j++) {
  • trunk/SophyaLib/HiStats/histos2.h

    r926 r943  
    6464
    6565  // get/put dans/depuis une matrice / vector
    66   void GetXCoor(Vector& v);
    67   void GetValue(Matrix &v);
    68   void GetYCoor(Vector& v);
    69   void GetError2(Matrix& v);
    70   void GetError(Matrix& v);
    71   void PutValue(Matrix& v, int ierr=0);
    72   void PutValueAdd(Matrix& v, int ierr=0);
    73   void PutError2(Matrix& v);
    74   void PutError2Add(Matrix& v);
    75   void PutError(Matrix& v);
     66  void GetXCoor(TVector<r_8>& v);
     67  void GetValue(TMatrix<r_8> &v);
     68  void GetYCoor(TVector<r_8>& v);
     69  void GetError2(TMatrix<r_8>& v);
     70  void GetError(TMatrix<r_8>& v);
     71  void PutValue(TMatrix<r_8>& v, int ierr=0);
     72  void PutValueAdd(TMatrix<r_8>& v, int ierr=0);
     73  void PutError2(TMatrix<r_8>& v);
     74  void PutError2Add(TMatrix<r_8>& v);
     75  void PutError(TMatrix<r_8>& v);
    7676
    7777  // INLINES
Note: See TracChangeset for help on using the changeset viewer.