VariableMetricBuilder.h

Go to the documentation of this file.

// @(#)root/minuit2:$Id$
// Authors: M. Winkler, F. James, L. Moneta, A. Zsenei   2003-2005
 
/**********************************************************************
 *                                                                    *
 * Copyright (c) 2005 LCG ROOT Math team,  CERN/PH-SFT                *
 *                                                                    *
 **********************************************************************/
 
#ifndef ROOT_Minuit2_VariableMetricBuilder
#define ROOT_Minuit2_VariableMetricBuilder
 
#include "Minuit2/MnConfig.h"
#include "Minuit2/MinimumBuilder.h"
#include "Minuit2/VariableMetricEDMEstimator.h"
#include "Minuit2/DavidonErrorUpdator.h"
#include "Minuit2/BFGSErrorUpdator.h"
 
#include <vector>
#include <memory>
 
namespace ROOT {
 
namespace Minuit2 {
 
/**
   Build (find) function minimum using the Variable Metric method (MIGRAD)
   Two possible error updators can be chosen
    - Davidon : this is the standard formula used in Migrad
    - BFGS this is the new formula based on BFGS algorithm
      (see Broyden–Fletcher–Goldfarb–Shanno algorithm
      https://en.wikipedia.org/wiki/Broyden–Fletcher–Goldfarb–Shanno_algorithm )
 */
class VariableMetricBuilder : public MinimumBuilder {
 
public:
   enum ErrorUpdatorType { kDavidon, kBFGS };
 
   VariableMetricBuilder(ErrorUpdatorType type = kDavidon) : fEstimator(VariableMetricEDMEstimator())
   {
      if (type == kBFGS)
         fErrorUpdator = std::unique_ptr<MinimumErrorUpdator>(new BFGSErrorUpdator());
      else
         fErrorUpdator = std::unique_ptr<MinimumErrorUpdator>(new DavidonErrorUpdator());
   }
 
   FunctionMinimum Minimum(const MnFcn &, const GradientCalculator &, const MinimumSeed &, const MnStrategy &,
                                   unsigned int, double) const override;
 
   FunctionMinimum Minimum(const MnFcn &, const GradientCalculator &, const MinimumSeed &, std::vector<MinimumState> &,
                           unsigned int, double) const;
 
   const VariableMetricEDMEstimator &Estimator() const { return fEstimator; }
   const MinimumErrorUpdator &ErrorUpdator() const { return *fErrorUpdator; }
 
   void AddResult(std::vector<MinimumState> &result, const MinimumState &state) const;
 
private:
   VariableMetricEDMEstimator fEstimator;
   std::shared_ptr<MinimumErrorUpdator> fErrorUpdator;
};
 
} // namespace Minuit2
 
} // namespace ROOT
 
#endif // ROOT_Minuit2_VariableMetricBuilder