Changeset 943 in Sophya
- Timestamp:
- Apr 15, 2000, 1:06:28 PM (25 years ago)
- Location:
- trunk/SophyaLib/HiStats
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaLib/HiStats/hisprof.h
r926 r943 30 30 // Acces a l information 31 31 //! 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;} 33 34 //! 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);} 35 37 //! 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);} 37 40 //! 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];} 39 43 //! 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];} 41 46 //! Retourne la dispersion/erreur du bin i 42 47 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;} 44 50 45 51 // Operators … … 62 68 //! Print, voir detail dans Histo::Print 63 69 inline void Print(int dyn=100,float hmin=1.,float hmax=-1.,int pflag=0,int il=1,int ih=-1) 64 70 {if(!Ok) UpdateHisto(); Histo::Print(dyn,hmin,hmax,pflag,il,ih);} 65 71 66 72 protected: -
trunk/SophyaLib/HiStats/histos.cc
r926 r943 1 1 // 2 // $Id: histos.cc,v 1. 3 2000-04-13 18:38:32ansari Exp $2 // $Id: histos.cc,v 1.4 2000-04-15 11:06:27 ansari Exp $ 3 3 // 4 4 … … 387 387 Remplissage d'un tableau avec la valeur des abscisses 388 388 */ 389 void Histo::GetAbsc( Vector&v)389 void Histo::GetAbsc(TVector<r_8> &v) 390 390 { 391 391 v.Realloc(bins); … … 397 397 Remplissage d'un tableau avec la valeur du contenu 398 398 */ 399 void Histo::GetValue( Vector&v)399 void Histo::GetValue(TVector<r_8> &v) 400 400 { 401 401 v.Realloc(bins); … … 407 407 Remplissage d'un tableau avec la valeur des erreurs au carre 408 408 */ 409 void Histo::GetError2( Vector&v)409 void Histo::GetError2(TVector<r_8> &v) 410 410 { 411 411 v.Realloc(bins); … … 418 418 Remplissage d'un tableau avec la valeur des erreurs 419 419 */ 420 void Histo::GetError( Vector&v)420 void Histo::GetError(TVector<r_8> &v) 421 421 { 422 422 v.Realloc(bins); … … 430 430 Remplissage du contenu de l'histo avec les valeurs d'un vecteur 431 431 */ 432 void Histo::PutValue( Vector&v, int ierr)432 void Histo::PutValue(TVector<r_8> &v, int ierr) 433 433 { 434 434 if(v.NElts()<bins) THROW(sizeMismatchErr); … … 443 443 Addition du contenu de l'histo avec les valeurs d'un vecteur 444 444 */ 445 void Histo::PutValueAdd( Vector&v, int ierr)445 void Histo::PutValueAdd(TVector<r_8> &v, int ierr) 446 446 { 447 447 if(v.NElts()<bins) THROW(sizeMismatchErr); … … 456 456 Remplissage des erreurs au carre de l'histo avec les valeurs d'un vecteur 457 457 */ 458 void Histo::PutError2( Vector&v)458 void Histo::PutError2(TVector<r_8> &v) 459 459 { 460 460 if(v.NElts()<bins) THROW(sizeMismatchErr); … … 467 467 Addition des erreurs au carre de l'histo avec les valeurs d'un vecteur 468 468 */ 469 void Histo::PutError2Add( Vector&v)469 void Histo::PutError2Add(TVector<r_8> &v) 470 470 { 471 471 if(v.NElts()<bins) THROW(sizeMismatchErr); … … 478 478 Remplissage des erreurs de l'histo avec les valeurs d'un vecteur 479 479 */ 480 void Histo::PutError( Vector&v)480 void Histo::PutError(TVector<r_8> &v) 481 481 { 482 482 if(v.NElts()<bins) THROW(sizeMismatchErr); … … 1065 1065 } 1066 1066 1067 VectorxFit(nLowHigh);1068 VectoryFit(nLowHigh);1069 Vectore2Fit(nLowHigh);1070 Vectorerrcoef(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); 1071 1071 int ii = 0; 1072 1072 for (int j=iLow; j<=iHigh; j++) { … … 1385 1385 if(f==NULL) 1386 1386 throw(NullPtrError("Histo::FitResidus: NULL pointer\n")); 1387 Vectorpar = gfit.GetParm();1387 TVector<r_8> par = gfit.GetParm(); 1388 1388 Histo h(*this); 1389 1389 for(int i=0;i<NBins();i++) { … … 1404 1404 if(f==NULL) 1405 1405 throw(NullPtrError("Histo::FitFunction: NULL pointer\n")); 1406 Vectorpar = gfit.GetParm();1406 TVector<r_8> par = gfit.GetParm(); 1407 1407 Histo h(*this); 1408 1408 for(int i=0;i<NBins();i++) { -
trunk/SophyaLib/HiStats/histos.h
r926 r943 1 1 // This may look like C code, but it is really -*- C++ -*- 2 2 // 3 // $Id: histos.h,v 1. 4 2000-04-13 18:38:33ansari Exp $3 // $Id: histos.h,v 1.5 2000-04-15 11:06:28 ansari Exp $ 4 4 // 5 5 … … 67 67 68 68 // 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); 78 78 79 79 // INLINES -
trunk/SophyaLib/HiStats/histos2.cc
r926 r943 603 603 Remplissage d'un tableau avec les valeurs des abscisses. 604 604 */ 605 void Histo2D::GetXCoor( Vector&v)605 void Histo2D::GetXCoor(TVector<r_8> &v) 606 606 { 607 607 float x,y; … … 614 614 Remplissage d'un tableau avec les valeurs des ordonnees. 615 615 */ 616 void Histo2D::GetYCoor( Vector&v)616 void Histo2D::GetYCoor(TVector<r_8> &v) 617 617 { 618 618 float x,y; … … 625 625 Remplissage d'un tableau avec les valeurs du contenu. 626 626 */ 627 void Histo2D::GetValue( Matrix&v)627 void Histo2D::GetValue(TMatrix<r_8> &v) 628 628 { 629 629 v.Realloc(nx,ny); … … 636 636 Remplissage d'un tableau avec les valeurs du carre des erreurs. 637 637 */ 638 void Histo2D::GetError2( Matrix&v)638 void Histo2D::GetError2(TMatrix<r_8> &v) 639 639 { 640 640 int i,j; … … 649 649 Remplissage d'un tableau avec les valeurs des erreurs. 650 650 */ 651 void Histo2D::GetError( Matrix&v)651 void Histo2D::GetError(TMatrix<r_8> &v) 652 652 { 653 653 int i,j; … … 663 663 Remplissage du contenu de l'histo avec les valeurs d'un tableau. 664 664 */ 665 void Histo2D::PutValue( Matrix&v, int ierr)665 void Histo2D::PutValue(TMatrix<r_8> &v, int ierr) 666 666 { 667 667 int i,j; … … 677 677 Addition du contenu de l'histo avec les valeurs d'un tableau. 678 678 */ 679 void Histo2D::PutValueAdd( Matrix&v, int ierr)679 void Histo2D::PutValueAdd(TMatrix<r_8> &v, int ierr) 680 680 { 681 681 int i,j; … … 692 692 avec les valeurs d'un tableau. 693 693 */ 694 void Histo2D::PutError2( Matrix&v)694 void Histo2D::PutError2(TMatrix<r_8> &v) 695 695 { 696 696 int i,j; … … 705 705 avec les valeurs d'un tableau. 706 706 */ 707 void Histo2D::PutError2Add( Matrix&v)707 void Histo2D::PutError2Add(TMatrix<r_8> &v) 708 708 { 709 709 int i,j; … … 718 718 Remplissage des erreurs de l'histo avec les valeurs d'un tableau. 719 719 */ 720 void Histo2D::PutError( Matrix&v)720 void Histo2D::PutError(TMatrix<r_8> &v) 721 721 { 722 722 int i,j; … … 1071 1071 if(f==NULL) 1072 1072 throw(NullPtrError("Histo2D::FitResidus: NULL pointer\n")); 1073 Vectorpar = gfit.GetParm();1073 TVector<r_8> par = gfit.GetParm(); 1074 1074 Histo2D h2(*this); 1075 1075 for(int i=0;i<NBinX();i++) for(int j=0;j<NBinY();j++) { … … 1092 1092 if(f==NULL) 1093 1093 throw(NullPtrError("Histo2D::FitFunction: NULL pointer\n")); 1094 Vectorpar = gfit.GetParm();1094 TVector<r_8> par = gfit.GetParm(); 1095 1095 Histo2D h2(*this); 1096 1096 for(int i=0;i<NBinX();i++) for(int j=0;j<NBinY();j++) { -
trunk/SophyaLib/HiStats/histos2.h
r926 r943 64 64 65 65 // 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); 76 76 77 77 // INLINES
Note:
See TracChangeset
for help on using the changeset viewer.