MINT::Minimiser Class Reference

#include <Minimiser.h>

Inheritance diagram for MINT::Minimiser:

Public Types
enum	FitStatus { NOTCALCULATED = 0, NOTACCURATE, NOTPOSDEF, CONVERGED }

Public Member Functions
	Minimiser (IMinimisable *fitFunction=0, const double errdef=1.)
virtual	~Minimiser ()
bool	attachFunction (IMinimisable *fcn)
IMinimisable *	theFunction ()
const IMinimisable *	theFunction () const
const MinuitParameterSet *	parSet () const
MinuitParameterSet *	parSet ()
unsigned int	nPars () const
IMinuitParameter *	getParPtr (unsigned int i)
const IMinuitParameter *	getParPtr (unsigned int i) const
bool	OK () const
bool	parsOK () const
bool	fcnOK () const
double	getFCNVal ()
void	FCNGradient (std::vector< double > &grad)
bool	initialiseVariables ()
bool	setPrintLevel (int level=-1)
bool	temporarilyQuiet ()
bool	resetPrintLevel ()
bool	SetSomeMinuitOptions ()
bool	CallMigrad ()
bool	CallMinos ()
bool	CallSeek (int maxCalls=100, int devs=5)
bool	CallSimplex (int maxCalls=300, double tolerance=1.)
bool	CallImprove (int maxCalls=1500, int searches=5)
bool	prepFit ()
bool	doFit ()
bool	doMinosFit ()
bool	doSeekFit (int maxCalls=100, int devs=5)
bool	doSimplexFit (int maxCalls=300, double tolerance=1.)
bool	scanMarked ()
bool	scanAll ()
TGraph *	scan (int i, double from=0, double to=0)
TGraph *	scan (IMinuitParameter &fp, double from=0, double to=0)
void	setMaxCalls (int maxCalls)
int	getMaxCalls () const
void	printResultVsInput (std::ostream &os=std::cout) const
TMatrixTSym< double >	covMatrix ()
TMatrixTSym< double >	covMatrixFull ()
int	status ()
bool	isStatus (Minimiser::FitStatus)
bool	isStatusAtLeast (Minimiser::FitStatus)
int	printStatus ()
bool	isConverged ()
TGraph *	contour (unsigned, unsigned, float nsigma=1., unsigned npoints=40)

Static Public Member Functions
static Minimiser *	getDefaultMinimiser ()

Protected Member Functions
bool	init ()
bool	MakeSpace (int needSpace)
void	TMinInit ()
bool	updateFitParameters (Double_t *p)
bool	setParametersToResult ()
bool	endOfFit ()
Int_t	Eval (Int_t npar, Double_t *grad, Double_t &fval, Double_t *par, Int_t flag)

Protected Attributes
bool	_useAnalyticGradient
MinuitParameterSet *	_parSet
IMinimisable *	_theFunction
int	_maxCalls
int	_printLevel
double	_errdef

Static Protected Attributes
static Minimiser *	_defaultMinimiser =0
static int	_defaultMaxCalls =100000

Private Attributes
Double_t	arglist [10]
Int_t	ierflg

Detailed Description

Definition at line 22 of file Minimiser.h.

Member Enumeration Documentation

FitStatus

enum MINT::Minimiser::FitStatus

Fit status enumerators.

Enumerator
NOTCALCULATED	FAILED: Covariance matrix not calculated at all.
NOTACCURATE	FAILED: Covariance matrix approximation only, not accurate.
NOTPOSDEF	FAILED: Full covariance matrix, but forced positive-definite.
CONVERGED	CONVERGED: Full accurate covariance matrix.

Definition at line 121 of file Minimiser.h.

                 {
    NOTCALCULATED = 0, 
    NOTACCURATE, 
    NOTPOSDEF, 
    CONVERGED 
  } ;

Constructor & Destructor Documentation

Minimiser()

Minimiser::Minimiser	(	IMinimisable *	fitFunction = 0,
		const double	errdef = 1.
	)

Definition at line 32 of file Minimiser.cpp.

  : TMinuit(0)
  , ierflg(0)
  , _useAnalyticGradient(false)
  , _theFunction(fitFunction)
  , _maxCalls(_defaultMaxCalls)
  , _printLevel(3)
  , _errdef(errdef)
{
  if(0 != theFunction()){
    init();
  }
}

~Minimiser()

Minimiser::~Minimiser ( )

virtual

Definition at line 46 of file Minimiser.cpp.

                     {
}

Member Function Documentation

attachFunction()

bool Minimiser::attachFunction

(

IMinimisable *

fcn

)

Definition at line 139 of file Minimiser.cpp.

                                               {
  if(0==fcn) return false;
  //  detachParameters();
  _theFunction = fcn;
  return init();
}

CallImprove()

bool Minimiser::CallImprove	(	int	maxCalls = 1500,
		int	searches = 5
	)

Definition at line 367 of file Minimiser.cpp.

                                                      {
    bool dbThis=true;
    bool success=true;
    bool useAnalyticGradient = _useAnalyticGradient;
    _useAnalyticGradient = false;
    arglist[0] = maxCalls; arglist[1] = searches;
    if(dbThis) cout << "calling IMPROVE" << endl;
    TMinuit::mnexcm("IMPROVE", arglist ,2,ierflg);
    if(dbThis) cout << "did that. How did I do? ierflg=" << ierflg << endl;
    success &= (! ierflg);
    _useAnalyticGradient = useAnalyticGradient;
    return success;
}

CallMigrad()

bool Minimiser::CallMigrad

(

)

Definition at line 317 of file Minimiser.cpp.

                          {
  bool dbThis=true;
  bool success=true;
  arglist[0] = _maxCalls; arglist[1] = 1.e-2;
  if(dbThis) cout << "calling MIGRAD" << endl;
  TMinuit::mnexcm("MIGRAD", arglist ,2,ierflg);
  if(dbThis) cout << "did that. How did I do? ierflg=" << ierflg << endl;
  success &= (! ierflg);
  return success;
}

CallMinos()

bool Minimiser::CallMinos

(

)

Definition at line 327 of file Minimiser.cpp.

                         {
  // thanks to Daniel Johnson for contributing this.
  bool dbThis=false;
  bool success=true;
  arglist[0] = _maxCalls; arglist[1] = 1.e-2;
  if(dbThis) cout << "calling MINOS" << endl;
  TMinuit::mnexcm("MINOS", arglist ,2,ierflg);
  if(dbThis) cout << "did that. How did I do? ierflg=" << ierflg << endl;
  success &= (! ierflg);
  return success;
}

CallSeek()

bool Minimiser::CallSeek	(	int	maxCalls = 100,
		int	devs = 5
	)

Definition at line 353 of file Minimiser.cpp.

                                              {
    bool dbThis=true;
    bool success=true;
    bool useAnalyticGradient = _useAnalyticGradient;
    _useAnalyticGradient = false;
    arglist[0] = maxCalls; arglist[1] = devs; 
    if(dbThis) cout << "calling SEEK" << endl;
    TMinuit::mnexcm("SEEK", arglist ,2,ierflg);
    if(dbThis) cout << "did that. How did I do? ierflg=" << ierflg << endl;
    success &= (! ierflg);
    _useAnalyticGradient = useAnalyticGradient;
    return success;
}

CallSimplex()

bool Minimiser::CallSimplex	(	int	maxCalls = 300,
		double	tolerance = 1.
	)

Definition at line 339 of file Minimiser.cpp.

                                                          {
    bool dbThis=true;
    bool success=true;
    bool useAnalyticGradient = _useAnalyticGradient;
    _useAnalyticGradient = false;
    arglist[0] = maxCalls; arglist[1] = tolerance;
    if(dbThis) cout << "calling SIMPLEX" << endl;
    TMinuit::mnexcm("SIMPLEX", arglist ,2,ierflg);
    if(dbThis) cout << "did that. How did I do? ierflg=" << ierflg << endl;
    success &= (! ierflg);
    _useAnalyticGradient = useAnalyticGradient;
    return success;
}

contour()

TGraph * Minimiser::contour	(	unsigned	iparx,
		unsigned	ipary,
		float	nsigma = 1.,
		unsigned	npoints = 40
	)

Make the nsigma contour plot for the given parameters - note that parameter numbers run from 0-N.

Definition at line 596 of file Minimiser.cpp.

                                                                                         {
  // Cache the initial values and errors.
  deque<deque<double> > initvals ;
  for(unsigned i = 0 ; i < parSet()->size() ; ++i){
    FitParameter* par = (FitParameter*)getParPtr(i) ;
    initvals.push_back(deque<double>()) ;
    deque<double>& parvals = initvals.back() ;
    parvals.push_back(par->mean()) ;
    parvals.push_back(par->err()) ;
    parvals.push_back(par->errPos()) ;
    parvals.push_back(par->errNeg()) ;
  }

  FitParameter* parx = (FitParameter*)getParPtr(iparx) ;
  FitParameter* pary = (FitParameter*)getParPtr(ipary) ;

  float errdef = _errdef * nsigma * nsigma ;
  SetErrorDef(errdef) ;
  TGraph* graph = (TGraph*)Contour(npoints, iparx, ipary) ;
  SetErrorDef(_errdef) ;
  ostringstream title ;
  title << pary->name() << "_vs_" << parx->name() << "_" << nsigma << "_sigma_contour" ;
  string titlestr(title.str()) ;
  graph->SetName(titlestr.c_str()) ;
  

  for(unsigned i = 0 ; i < parSet()->size() ; ++i)
    getParPtr(i)->setResult(initvals[i][0], initvals[i][1], initvals[i][2], initvals[i][3]) ;
  theFunction()->parametersChanged() ;

  return graph ;
}

covMatrix()

TMatrixTSym< double > Minimiser::covMatrix

(

)

Definition at line 506 of file Minimiser.cpp.

                                         {
  // only entries for variable parameters.
  // if you want rows and columns to correspond
  // to parameter numbers in MinuitParameterSets
  // use covMatrixFull.
  unsigned int internalPars = fNpar;

  TwoDArray<double> m1(internalPars, internalPars, 0.0);

  this->mnemat(&m1[0][0], internalPars);
  TMatrixTSym<double> matrix(internalPars);
  for(unsigned int i=0; i < internalPars; i++){
    for(unsigned int j=i; j < internalPars; j++){
      matrix(i,j) = matrix(j,i) = m1[i][j];
    }
  }
  return matrix;
}

covMatrixFull()

TMatrixTSym< double > Minimiser::covMatrixFull

(

)

Definition at line 525 of file Minimiser.cpp.

                                             {
  // rows and columns correspond to external parameter numbers
  // (if you have many fixed parameters, you'll get a lot of zeros)

  bool dbThis=false;

  unsigned int internalPars = fNpar;

  TwoDArray<double> m1(internalPars, internalPars, 0.0);

  this->mnemat(&m1[0][0], internalPars);
  TMatrixTSym<double> matrix(nPars());
  for(unsigned int i=0; i < internalPars; i++){
    for(unsigned int j=i; j < internalPars; j++){
      int ex_i = fNexofi[i] -1;
      int ex_j = fNexofi[j] -1;
      if(dbThis){
    cout << "Minimiser::covMatrixFull(): doing: "
         << "matrix(" << ex_i<< ", " << ex_j << " )"
         << " = matrix(" << ex_j << ", " << ex_i << ")"
         << " = m1[" << i << "][" << j<< "];" 
         << " = " << m1[i][j]
         << endl;
      }
      matrix(ex_i,ex_j) = matrix(ex_j,ex_i) = m1[i][j];
    }
  }
  return matrix;
}

doFit()

bool Minimiser::doFit

(

)

Definition at line 391 of file Minimiser.cpp.

                     {
  bool dbThis=false;
  bool success = true;
  if(dbThis) cout << "Minimiser::doFit() called" << endl;
  success &= prepFit();
  if(dbThis) cout << "... called prepFit " << success << endl;
  success &= CallMigrad();
  if(dbThis) cout << "... called MIGRAD" << success << endl;
  scanMarked();
  if(dbThis) cout << "... scanned" << success << endl;
  success &= this->endOfFit();
  return success;
}

doMinosFit()

bool Minimiser::doMinosFit

(

)

Definition at line 404 of file Minimiser.cpp.

                          {
  bool dbThis=false;
  bool success = true;
  success &= prepFit();
  success &= CallMinos();
  if(dbThis) cout << "called MINOS" << endl;
  scanMarked();
  success &= this->endOfFit();
  return success;
}

doSeekFit()

bool Minimiser::doSeekFit	(	int	maxCalls = 100,
		int	devs = 5
	)

Definition at line 415 of file Minimiser.cpp.

                                               {
    bool dbThis=false;
    bool success = true;
    success &= prepFit();
    success &= CallSeek(maxCalls, devs);
    if(dbThis) cout << "called SEEK" << endl;
    //scanMarked();
    success &= this->endOfFit();
    return success;
}

doSimplexFit()

bool Minimiser::doSimplexFit	(	int	maxCalls = 300,
		double	tolerance = 1.
	)

Definition at line 426 of file Minimiser.cpp.

                                                          {
    bool dbThis=false;
    bool success = true;
    success &= prepFit();
    success &= CallSimplex(maxCalls, tolerance);
    if(dbThis) cout << "called SIMPLEX" << endl;
    //scanMarked();
    success &= this->endOfFit();
    return success;
}

endOfFit()

bool Minimiser::endOfFit ( )

protected

Definition at line 250 of file Minimiser.cpp.

                        {
  if(! parsOK())return false;
  setParametersToResult();
  theFunction()->endFit();
  parSet()->setCovMatrix(covMatrixFull()) ;
  return true;
}

Eval()

Int_t Minimiser::Eval ( Int_t  npar, Double_t *  grad, Double_t &  fval, Double_t *  par, Int_t  flag  )

protected

Definition at line 56 of file Minimiser.cpp.

               {
  bool dbThis=false;
  if(dbThis) cout << "Eval got called " << endl;
  if(! this->OK()){
    cout << "ERROR in Minimiser::Eval (file: Minimiser.C)"
     << " Got called although I'm not OK" << endl;
    fval = -9999;
    return -1;
  }
  this->updateFitParameters(par);
  fval = this->getFCNVal();
  if (flag == 4) {
    if(_useAnalyticGradient){
      //calculate GRAD, the first derivatives of FVAL
      vector<double> gradient(npar);
      this->FCNGradient(gradient);
      for(unsigned int i = 0; i < static_cast<unsigned int>(npar); i++){
    grad[i]=gradient[i];
      }
    }
  }
  return -1;
  (void)npar;
}

FCNGradient()

void Minimiser::FCNGradient

(

std::vector< double > &

grad

)

Definition at line 281 of file Minimiser.cpp.

                                               {
    if(! this->OK()){
        cout << "ERROR IN Minimiser::FCNGradient()"
        << " I'm not OK!!" << endl;
    }
    return theFunction()->Gradient(grad);
}

fcnOK()

bool Minimiser::fcnOK

(

)

const

Definition at line 263 of file Minimiser.cpp.

                           {
  if(0 == theFunction()) return false;
  return true;
}

getDefaultMinimiser()

Minimiser * Minimiser::getDefaultMinimiser ( )

static

Definition at line 23 of file Minimiser.cpp.

                                         {

  if(0 == _defaultMinimiser){
    _defaultMinimiser = new Minimiser;
  }
  return _defaultMinimiser;
}

getFCNVal()

double Minimiser::getFCNVal

(

)

Definition at line 272 of file Minimiser.cpp.

                           {
  if(! this->OK()){
    cout << "ERROR IN Minimiser::getFCNVal()"
     << " I'm not OK!!" << endl;
    return -9999.0;
  }
  return theFunction()->getVal();
}

getMaxCalls()

int Minimiser::getMaxCalls

(

)

const

Definition at line 52 of file Minimiser.cpp.

                                {
  return _maxCalls;
}

getParPtr() [1/2]

IMinuitParameter * Minimiser::getParPtr

(

unsigned int

i

)

Definition at line 149 of file Minimiser.cpp.

                                                    {
  if(0 == _parSet) return 0;
  if(i >= nPars()) return 0;
  return _parSet->getParPtr(i);
}

getParPtr() [2/2]

const IMinuitParameter * Minimiser::getParPtr

(

unsigned int

i

)

const

Definition at line 154 of file Minimiser.cpp.

                                                                {
  if(0 == _parSet) return 0;
  if(i >= nPars()) return 0;
  return _parSet->getParPtr(i);
}

init()

bool Minimiser::init ( )

protected

Definition at line 105 of file Minimiser.cpp.

                    {
  bool dbThis=true;
  if(dbThis) cout << "Minimiser::init(): you called me" << endl;
  if(0 == theFunction()){
    std::cout << "ERROR IN Minimiser::init():"
          << " the function-ptr is empty."
          << " I won't do anything for now."
          << std::endl;
    return false;
  }
  if(dbThis) cout << " ... calling theFunction()->beginFit()" << endl;
  theFunction()->beginFit();

  _parSet = theFunction()->getParSet();
  if(0 == _parSet){
    std::cout << "ERROR IN Minimiser::init():"
          << " theFunction()->getParSet() returned empty pointer"
          << " I won't do anything for now."
          << std::endl;
    return false;
  }

  MakeSpace(nPars());

  if(dbThis) cout << "... made space, now initialising variables" << endl;
  initialiseVariables();
  if(dbThis) cout << " initialised variables." << endl;

  _useAnalyticGradient = theFunction()->useAnalyticGradient();  
    
  if(dbThis) cout << "Minimiser::init(): returning true" << endl;
  return true;
}

initialiseVariables()

bool Minimiser::initialiseVariables

(

)

Definition at line 170 of file Minimiser.cpp.

                                   {
  bool dbThis=false;
  if(! parsOK()) return false;
  
  bool success=true;
  if(dbThis) {
    cout << "Minimiser::initialiseVariables() called. " << endl;
    cout << "\n\t(declaring them to MINUIT)" << endl;
  } 

  //remove hidden MinuitParameter
  for(unsigned int i=0; i < nPars(); i++){
    if(getParPtr(i)->hidden()){
        if(dbThis) cout << i << ")" << getParPtr(i)->name() << endl;
    parSet()->unregister(getParPtr(i));
    i--;
     }
  }

  //temporarilyQuiet();
  for(unsigned int i=0; i < nPars(); i++){
    int ierflag=0;
    //if(! getParPtr(i)->hidden()){        //hidden MinuitParamter already removed
      if(dbThis) cout << i << ")" << getParPtr(i)->name() << endl;
      double step = getParPtr(i)->stepInit();
      if(getParPtr(i)->iFixInit()) step=0;
      this->mnparm( i
            , getParPtr(i)->name().c_str()
            , getParPtr(i)->meanInit()
            , step
            , getParPtr(i)->minInit()
            , getParPtr(i)->maxInit()
            , ierflag);
      //if(getParPtr(i)->iFixInit()) FixParameter(i);
    //}
    //getParPtr(i)->associate(this, i);
    success &= ! ierflag;
  }
  if(dbThis){
    cout << "Minimiser::initialiseVariables():"
     << "done declaring variables, now:"
     << "theFunction()->parametersChanged();"
     << endl;
  }
  theFunction()->parametersChanged();
  if(dbThis) cout << "Minimiser::initialiseVariables():"
          << " all done - returning " << success << endl;
  //resetPrintLevel();
  return success;
}

isConverged()

bool Minimiser::isConverged

(

)

Check if the fit status is CONVERGED

Definition at line 592 of file Minimiser.cpp.

                            {
  return isStatus(Minimiser::CONVERGED) ;
}

isStatus()

bool Minimiser::isStatus

(

Minimiser::FitStatus

stat

)

Check if the fit status is of the given type.

Definition at line 562 of file Minimiser.cpp.

                                                {
  return status() == stat ;
}

isStatusAtLeast()

bool Minimiser::isStatusAtLeast

(

Minimiser::FitStatus

stat

)

Check if the fit status is >= the given status type.

Definition at line 566 of file Minimiser.cpp.

                                                       {
  return status() >= stat ;
}

MakeSpace()

bool Minimiser::MakeSpace ( int  needSpace )

protected

Definition at line 97 of file Minimiser.cpp.

                                      {
  //return true;

  this->DeleteArrays();
  this->BuildArrays(needSpace + 2); // +2 for good luck
  TMinInit();
  return true;
}

nPars()

unsigned int Minimiser::nPars

(

)

const

Definition at line 145 of file Minimiser.cpp.

                                   {
  if(0 == _parSet) return 0;
  return _parSet->size();
}

OK()

bool Minimiser::OK

(

)

const

Definition at line 268 of file Minimiser.cpp.

                        {
  return parsOK() && fcnOK();
}

parSet() [1/2]

const MinuitParameterSet* MINT::Minimiser::parSet ( ) const

inline

Definition at line 63 of file Minimiser.h.

{ return _parSet;}

parSet() [2/2]

MinuitParameterSet* MINT::Minimiser::parSet ( )

inline

Definition at line 64 of file Minimiser.h.

{ return _parSet;}

parsOK()

bool Minimiser::parsOK

(

)

const

Definition at line 258 of file Minimiser.cpp.

                            {
  if(0 == _parSet) return false;
  if(0 == nPars()) return false;
  return true;
}

prepFit()

bool Minimiser::prepFit

(

)

Definition at line 381 of file Minimiser.cpp.

                       {
  bool dbThis=false;
  bool success=true;
  success &= SetSomeMinuitOptions();
  theFunction()->beginFit();
  success &= this->initialiseVariables();
  if(dbThis) cout << "re-initialised variables" << endl;
  return success;
}

printResultVsInput()

void Minimiser::printResultVsInput

(

std::ostream &

os = std::cout

)

const

Definition at line 501 of file Minimiser.cpp.

                                                      {
  if(0 == _parSet) return;
  _parSet->printResultVsInput(os);
}

printStatus()

int Minimiser::printStatus

(

)

Print the fit status.

Definition at line 570 of file Minimiser.cpp.

                           {
  int stat = status() ;
  cout << "Fit status: " << stat << endl ;
  switch(stat) {
    case NOTCALCULATED :
      cout << "FAILED: Covariance matrix not calculated at all." << endl ;
      break ;
    case NOTACCURATE :
      cout << "FAILED: Covariance matrix approximation only, not accurate." << endl ;
      break ;
    case NOTPOSDEF :
      cout << "FAILED: Full covariance matrix, but forced positive-definite." << endl ;
      break ;
    case CONVERGED :
      cout << "CONVERGED: Full accurate covariance matrix." << endl ;
      break ;
    default :
      cout << "UNKNOWN." << endl ;
    }
  return stat ;
}

resetPrintLevel()

bool Minimiser::resetPrintLevel

(

)

Definition at line 300 of file Minimiser.cpp.

                               {
  return setPrintLevel(-9999);
}

scan() [1/2]

TGraph * Minimiser::scan	(	int	i,
		double	from = 0,
		double	to = 0
	)

note that the index i is the one in the parameter list which goes from 0 to n-1 (i.e. C-style). This corresponds to Minuit's parameter number i+1. (so if you want Minuit's fit parameter 1, pass it 1-1=0)

Definition at line 473 of file Minimiser.cpp.

                                                    {
  // note that the index i is the one in the parameter list
  // which goes from 0 to n-1 (i.e. C-style).
  // This corresponds to Minuit's parameter number i+1.
  // (so if you want Minuit's fit parameter 1, pass it 1-1=0)
  IMinuitParameter* p = getParPtr(i);
  if(0 == p) return 0;

  double points = 100;
  string fname = "scan_" + p->name() + ".root";
  
  Double_t arglist[10] = {0};
  Int_t ierflg=0;

  arglist[0]=i+1; arglist[1]=points; arglist[2]=from; arglist[3]=to;
  this->mnexcm("SCAN", arglist, 4, ierflg);
  TGraph *gr = (TGraph*) this->GetPlot();
  if(0 == gr){
    cout << " didn't get plot " << endl;
  }else{
    TFile fscan(fname.c_str(), "RECREATE");
    fscan.cd();
    gr->Write();
    fscan.Close();
  }
  return gr;
}

scan() [2/2]

TGraph * Minimiser::scan	(	IMinuitParameter &	fp,
		double	from = 0,
		double	to = 0
	)

Definition at line 459 of file Minimiser.cpp.

                                                                   {
  if(fp.parSet() != this->parSet()){
    cout << "ERROR Minimiser::scanParameter: "
     << " You are trying to scan a parameter that is"
     << " associated to a different fit (different fit parameter set)"
     << " will not do it."
     << " The parameter you tried was: " << endl;
    fp.print();
    cout <<  endl;
    return 0;
  }
  return scan(fp.parSetIndex(), from, to);
}

scanAll()

bool Minimiser::scanAll

(

)

Definition at line 448 of file Minimiser.cpp.

                       {
  bool sc=true;
  if(0 == nPars()) return false;
  for(unsigned int i=0; i < nPars(); i++){
    IMinuitParameter* p=getParPtr(i);
    if(0 == p) continue;
    sc &= (bool) scan(i, p->scanMin(), p->scanMax());
  }
  return sc;
}

scanMarked()

bool Minimiser::scanMarked

(

)

Definition at line 438 of file Minimiser.cpp.

                          {
  bool sc=true;
  if(0 == nPars()) return false;
  for(unsigned int i=0; i < nPars(); i++){
    IMinuitParameter* p=getParPtr(i);
    if(0 == p || (! p->scan()) ) continue;
    sc &= (bool) scan(i, p->scanMin(), p->scanMax());
  }
  return sc;
}

setMaxCalls()

void Minimiser::setMaxCalls

(

int

maxCalls

)

Definition at line 49 of file Minimiser.cpp.

                                       {
  _maxCalls = maxCalls;
}

setParametersToResult()

bool Minimiser::setParametersToResult ( )

protected

Definition at line 221 of file Minimiser.cpp.

                                     {
  if(! parsOK())return false;

  Double_t mean, err, errN, errP;
  Double_t gcc;
  for(unsigned int i=0; i < nPars(); i++){
    //if(! getParPtr(i)->hidden()){
      this->mnerrs(i, errP, errN, err, gcc);
      this->GetParameter(i, mean, err);
      mean += this->getParPtr(i)->blinding();
      this->getParPtr(i)->setResult(mean, err, errP, errN);
    //}
  }
  theFunction()->parametersChanged();
  return true;
}

setPrintLevel()

bool Minimiser::setPrintLevel

(

int

level = -1

)

Definition at line 289 of file Minimiser.cpp.

                                      {
  if(level >=0 )_printLevel=level;
  arglist[0] = _printLevel;
  TMinuit::mnexcm("SET PRINTOUT", arglist , 1, ierflg);
  return (! ierflg);
}

SetSomeMinuitOptions()

bool Minimiser::SetSomeMinuitOptions

(

)

Definition at line 303 of file Minimiser.cpp.

                                    {
  bool success = true;
  success &= setPrintLevel();
  arglist[0] = _errdef ;
  TMinuit::mnexcm("SET ERR", arglist , 1, ierflg);
  success &= (! ierflg);
  arglist[0] = 1;
  TMinuit::mnexcm("SET STRATEGY", arglist , 1, ierflg);
  success &= (! ierflg);
  if(_useAnalyticGradient)TMinuit::mnexcm("SET GRADIENT", arglist , 1, ierflg);
  else TMinuit::mnexcm("SET NOGRADIENT", arglist , 1, ierflg);
  return success;
}

status()

int Minimiser::status

(

)

Get the fit status

Definition at line 555 of file Minimiser.cpp.

                      {
  double fmin, fedm, errdef ;
  int npari, nparx, status ;
  mnstat( fmin, fedm, errdef, npari, nparx, status ) ;
  return status ;
}

temporarilyQuiet()

bool Minimiser::temporarilyQuiet

(

)

Definition at line 295 of file Minimiser.cpp.

                                {
  arglist[0] = -1;
  TMinuit::mnexcm("SET PRINTOUT", arglist , 1, ierflg);
  return (! ierflg);
}

theFunction() [1/2]

IMinimisable* MINT::Minimiser::theFunction ( )

inline

Definition at line 60 of file Minimiser.h.

{ return _theFunction;}

theFunction() [2/2]

const IMinimisable* MINT::Minimiser::theFunction ( ) const

inline

Definition at line 61 of file Minimiser.h.

{ return _theFunction;}

TMinInit()

void Minimiser::TMinInit ( )

protected

Definition at line 86 of file Minimiser.cpp.

                        {
   fStatus       = 0;
   fEmpty        = 0;
   fObjectFit    = 0;
   fMethodCall   = 0;
   fPlot         = 0;
   fGraphicsMode = kTRUE;
   SetMaxIterations();
   mninit(5,6,7);
}

updateFitParameters()

bool Minimiser::updateFitParameters ( Double_t *  p )

protected

Definition at line 238 of file Minimiser.cpp.

                                                {
  if(! parsOK())return false;

  for(unsigned int i=0; i < nPars(); i++){
    //if(! getParPtr(i)->hidden()){
      //getParPtr(i)->setCurrentFitVal(par[i]);
       getParPtr(i)->setCurrentFitVal(par[i] + this->getParPtr(i)->blinding());
    //}
  }
  theFunction()->parametersChanged();
  return true;
}

Member Data Documentation

_defaultMaxCalls

int Minimiser::_defaultMaxCalls =100000

staticprotected

Definition at line 28 of file Minimiser.h.

_defaultMinimiser

Minimiser * Minimiser::_defaultMinimiser =0

staticprotected

Definition at line 27 of file Minimiser.h.

_errdef

double MINT::Minimiser::_errdef

protected

Definition at line 37 of file Minimiser.h.

_maxCalls

int MINT::Minimiser::_maxCalls

protected

Definition at line 34 of file Minimiser.h.

_parSet

MinuitParameterSet* MINT::Minimiser::_parSet

protected

Definition at line 30 of file Minimiser.h.

_printLevel

int MINT::Minimiser::_printLevel

protected

Definition at line 36 of file Minimiser.h.

_theFunction

IMinimisable* MINT::Minimiser::_theFunction

protected

Definition at line 33 of file Minimiser.h.

_useAnalyticGradient

bool MINT::Minimiser::_useAnalyticGradient

protected

Definition at line 26 of file Minimiser.h.

arglist

Double_t MINT::Minimiser::arglist[10]

mutableprivate

Definition at line 23 of file Minimiser.h.

ierflg

Int_t MINT::Minimiser::ierflg

mutableprivate

Definition at line 24 of file Minimiser.h.

---

The documentation for this class was generated from the following files:

• Mint/Minimiser.h
• src/Mint/Minimiser/Minimiser.cpp