| [756] | 1 | // This may look like C code, but it is really -*- C++ -*-
|
|---|
| 2 | //
|
|---|
| 3 | // Methode d'ajustement non-lineaire sur des objets
|
|---|
| 4 | // R. Ansari 02/2000
|
|---|
| 5 |
|
|---|
| 6 | #ifndef OBJFITTER_SEEN
|
|---|
| 7 | #define OBJFITTER_SEEN
|
|---|
| 8 |
|
|---|
| 9 | #include "machdefs.h"
|
|---|
| 10 |
|
|---|
| 11 | #include "tmatrix.h"
|
|---|
| 12 | #include "tvector.h"
|
|---|
| [1110] | 13 | #include "histos.h"
|
|---|
| 14 | #include "histos2.h"
|
|---|
| 15 | #include "hisprof.h"
|
|---|
| [756] | 16 | #include "generalfit.h"
|
|---|
| 17 |
|
|---|
| 18 |
|
|---|
| 19 | namespace SOPHYA {
|
|---|
| 20 |
|
|---|
| 21 | class ObjectFitter {
|
|---|
| 22 | public:
|
|---|
| 23 | // Residus et fonction fittees sur matrix
|
|---|
| 24 | static TMatrix<r_4> FitResidus(TMatrix<r_4> const & mtx, GeneralFit& gfit,
|
|---|
| 25 | double xorg=0.,double yorg=0.,double dx=1.,double dy=1.);
|
|---|
| 26 | static TMatrix<r_8> FitResidus(TMatrix<r_8> const & mtx, GeneralFit& gfit,
|
|---|
| 27 | double xorg=0.,double yorg=0.,double dx=1.,double dy=1.);
|
|---|
| 28 | static TMatrix<r_4> FitFunction(TMatrix<r_4> const & mtx, GeneralFit& gfit,
|
|---|
| 29 | double xorg=0.,double yorg=0.,double dx=1.,double dy=1.);
|
|---|
| 30 | static TMatrix<r_8> FitFunction(TMatrix<r_8> const & mtx, GeneralFit& gfit,
|
|---|
| 31 | double xorg=0.,double yorg=0.,double dx=1.,double dy=1.);
|
|---|
| 32 | // Residus et fonction fittees sur vecteurs
|
|---|
| 33 | static TVector<r_4> FitResidus(TVector<r_4> const & vec, GeneralFit& gfit,
|
|---|
| 34 | double xorg=0.,double dx=1.);
|
|---|
| 35 | static TVector<r_8> FitResidus(TVector<r_8> const & vec, GeneralFit& gfit,
|
|---|
| 36 | double xorg=0.,double dx=1.);
|
|---|
| 37 | static TVector<r_4> FitFunction(TVector<r_4> const & vec, GeneralFit& gfit,
|
|---|
| 38 | double xorg=0.,double dx=1.);
|
|---|
| 39 | static TVector<r_8> FitFunction(TVector<r_8> const & vec, GeneralFit& gfit,
|
|---|
| 40 | double xorg=0.,double dx=1.);
|
|---|
| [1110] | 41 | // Residus et fonction fittees sur Histo
|
|---|
| 42 | static Histo FitResidus(Histo const& h, GeneralFit& gfit);
|
|---|
| 43 | static Histo FitFunction(Histo const& h, GeneralFit& gfit);
|
|---|
| 44 | static int_4 Fit(Histo const & h, GeneralFit& gfit,unsigned short typ_err=0);
|
|---|
| 45 | // Residus et fonction fittees sur Histo2D
|
|---|
| 46 | static Histo2D FitResidus(Histo2D const& h, GeneralFit& gfit);
|
|---|
| 47 | static Histo2D FitFunction(Histo2D const& h, GeneralFit& gfit);
|
|---|
| 48 | static int_4 Fit(Histo2D const & h, GeneralFit& gfit,unsigned short typ_err=0);
|
|---|
| 49 | // Residus et fonction fittees sur HProf
|
|---|
| 50 | static inline Histo FitResidus(HProf const& h, GeneralFit& gfit)
|
|---|
| 51 | {h.UpdateHisto(); return FitResidus((Histo const&) h,gfit);}
|
|---|
| 52 | static inline Histo FitFunction(HProf const& h, GeneralFit& gfit)
|
|---|
| 53 | {h.UpdateHisto(); return FitFunction((Histo const&) h,gfit);}
|
|---|
| 54 | static int_4 Fit(HProf const & h, GeneralFit& gfit)
|
|---|
| 55 | {h.UpdateHisto(); return Fit((Histo const&) h,gfit,0);}
|
|---|
| 56 | // Residus et fonction fittees sur GeneralFitData
|
|---|
| 57 | static inline GeneralFitData FitResidus(GeneralFitData const& g, GeneralFit& gfit)
|
|---|
| 58 | {return g.FitResidus(gfit);}
|
|---|
| 59 | static inline GeneralFitData FitFunction(GeneralFitData const& g, GeneralFit& gfit)
|
|---|
| 60 | {return g.FitFunction(gfit);}
|
|---|
| [756] | 61 | };
|
|---|
| 62 |
|
|---|
| 63 | // Classe template pour ajustement sur tableaux (Vecteurs, Matrices)
|
|---|
| 64 | template <class T>
|
|---|
| 65 | class ArrayFitter {
|
|---|
| 66 | public:
|
|---|
| 67 | // Residus et fonction fittees sur matrix
|
|---|
| 68 | static TMatrix<T> FitResidus(TMatrix<T> const & mtx, GeneralFit& gfit,
|
|---|
| 69 | double xorg=0.,double yorg=0.,double dx=1.,double dy=1.);
|
|---|
| 70 | static TMatrix<T> FitFunction(TMatrix<T> const & mtx, GeneralFit& gfit,
|
|---|
| 71 | double xorg=0.,double yorg=0.,double dx=1.,double dy=1.);
|
|---|
| 72 | // Residus et fonction fittees sur vecteurs
|
|---|
| 73 | static TVector<T> FitResidus(TVector<T> const & vec, GeneralFit& gfit,
|
|---|
| 74 | double xorg=0.,double dx=1.);
|
|---|
| 75 | static TVector<T> FitFunction(TVector<T> const & vec, GeneralFit& gfit,
|
|---|
| 76 | double xorg=0.,double dx=1.);
|
|---|
| 77 | };
|
|---|
| 78 |
|
|---|
| 79 | } // Fin du namespace
|
|---|
| 80 |
|
|---|
| 81 |
|
|---|
| [805] | 82 | #endif // OBJFITTER_SEEN
|
|---|