ClsqSolver.hh~

#ifndef CLSQSOLVER_HH
#define CLSQSOLVER_hh

// Solver for constrained linear least squares

#include "Constraints.hh"
#include "ConstraintFunction.hh"

#include "TVectorD.h"
#include "TMatrixD.h"

#include <string>
#include <map>
#include <vector>

class ClsqSolver {

public:

  // Ctor: pass input data, constraint function object and optional
  // parameter name maps mapping parameter number (0-n) to names:
  ClsqSolver( const TVectorD& d, 
	      const TMatrixDSym& c, 
	      const TVectorD& u, 
	      const ConstraintFunction& cfun, 
	      const std::map<int,std::string>& upn=std::map<int,std::string>(),
	      const std::map<int,std::string>& mpn=std::map<int,std::string>(), 
	      Int_t nd=0, Int_t mi=100, 
	      Double_t eps=0.0001, Double_t dc=0.0001 );
  
  // Accessors:
  std::vector<std::string> getMParNames() const;
  std::vector<std::string> getUParNames() const;

private:
  
  // Helper to setup parameter names
  std::vector<std::string> 
  setParameterNames( const std::map<int,std::string>& parnames,
		     Int_t npar );
			  
  // Instance variables:
  TVectorD data;
  TMatrixDSym covm;
  TVectorD upar;
  TVectorD mpar;
  Constraints constraints;
  std::vector<std::string> uparnames;
  std::vector<std::string> mparnames;
  TMatrixDSym invm;
  Int_t ndof;
  Int_t maxiterations;
  Double_t epsilon;
  Double_t deltachi2;
  Int_t niterations;

};


#endif