Fitter.cxx

Go to the documentation of this file.

#include "Fitter.h"

#include "StatedFCN.h"

#include "FitterException.h"

using std::string;
using std::vector;

Fitter::
Fitter ( const char * name )
  : m_name ( name ),
    m_fcn ( 0 ),
    m_max_iterations ( 100 )
{
}

Fitter::
Fitter ( const Fitter & fitter )
  : m_name ( fitter.m_name ),
    m_max_iterations ( fitter.m_max_iterations )
{
  if ( fitter.m_fcn != 0 ) m_fcn = fitter.m_fcn -> clone ();
}

Fitter::
~Fitter ()
{
  if ( m_fcn != 0 ) delete m_fcn;
}

void
Fitter::
setFCN ( StatedFCN * fcn )
{
  if ( m_fcn != 0 ) delete m_fcn;

  m_fcn = fcn;
}

StatedFCN *
Fitter::
getFCN ()
{
  return m_fcn;
}

const std::string &
Fitter::
name () const
{
  return m_name;
}

void
Fitter::
fillFreeParameters ( std::vector < double > & free_parms) const
{
  return m_fcn -> fillFreeParameters ( free_parms );
}

void
Fitter::
setFixedFlags ( const std::vector < int > & flags )
{
  m_fcn -> setFixedFlags ( flags );
}

const vector < int > &
Fitter::
getFixedFlags ( ) const
{
  return m_fcn -> getFixedFlags ();
}

void
Fitter::
setLimits ( unsigned int i, double lower, double upper )
{
  string what ( "This " );
  what += m_name;
  what += " minimizer does not support limits on parameters";
  throw FitterException ( what );
}

unsigned int
Fitter::
getParameterIndex ( const std::string & name )
{
  unsigned int index = UINT_MAX;
  const vector < string > & names = m_fcn -> getParmNames ();
  unsigned int size = names.size();
  for ( unsigned int i = 0; i < size; i++ ) {
    const string & pname = names [i];
    if ( pname == name ) {
      index = i;
      break;
    }
  }
  if ( index == UINT_MAX ) {
    string what ( "No parameter named `" );
    what += name;
    what += "' for this function.";
    throw FitterException ( what );
  }

  return index;
}

void
Fitter::
setLimits ( const std::string & name, double lower, double upper )
{
//   unsigned int index = UINT_MAX;
//   const vector < string > & names = m_fcn -> getParmNames ();
//   unsigned int size = names.size();
//   for ( unsigned int i = 0; i < size; i++ ) {
//     const string & pname = names [i];
//     if ( pname == name ) {
//       index = i;
//       break;
//     }
//   }
  unsigned int index = getParameterIndex ( name );
  setLimits ( index, lower, upper );
}

void
Fitter::
setStepSize ( unsigned int i, double size )
{
  string what ( "This " );
  what += m_name;
  what += " minimizer does not support setting step size.";
  throw FitterException ( what );
}

void
Fitter::
setStepSize ( const std::string & name, double size )
{
  unsigned int index = getParameterIndex ( name );
  setStepSize ( index, size );
}

double
Fitter::objectiveValue () const
{
  return m_fcn -> objectiveValue ();
}

int
Fitter::
calcDegreesOfFreedom () const
{
  return m_fcn -> degreesOfFreedom();
}

int
Fitter::
calcCovariance ( std::vector < std::vector < double > > & cov )
{
  return EXIT_FAILURE;
}

void
Fitter::
setFitCut ( TupleCut * cut )
{
  m_fcn -> setFitCut ( cut );
}

void
Fitter::
setFitRange ( bool yes )
{
  m_fcn -> setFitRange ( yes );
}

---