SimplexParameters.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_SimplexParameters
#define ROOT_Minuit2_SimplexParameters
 
#include <cassert>
 
#include "Minuit2/MnMatrix.h"
 
#include <ROOT/RSpan.hxx>
 
#include <vector>
#include <utility>
 
namespace ROOT {
 
namespace Minuit2 {
 
/**
   class describing the simplex set of points (f(x), x )  which evolve during the minimization
   iteration process.
 */
 
class SimplexParameters {
 
public:
   SimplexParameters(std::span<const std::pair<double, MnAlgebraicVector>> simpl, unsigned int jh, unsigned int jl)
      : fSimplexParameters(simpl.begin(), simpl.end()), fJHigh(jh), fJLow(jl)
   {
   }
 
   ~SimplexParameters() {}
 
   void Update(double, const MnAlgebraicVector &);
 
   const std::vector<std::pair<double, MnAlgebraicVector>> &Simplex() const { return fSimplexParameters; }
 
   const std::pair<double, MnAlgebraicVector> &operator()(unsigned int i) const
   {
      assert(i < fSimplexParameters.size());
      return fSimplexParameters[i];
   }
 
   unsigned int Jh() const { return fJHigh; }
   unsigned int Jl() const { return fJLow; }
   double Edm() const { return fSimplexParameters[Jh()].first - fSimplexParameters[Jl()].first; }
   MnAlgebraicVector Dirin() const;
 
private:
   std::vector<std::pair<double, MnAlgebraicVector>> fSimplexParameters;
   unsigned int fJHigh;
   unsigned int fJLow;
};
 
} // namespace Minuit2
 
} // namespace ROOT
 
#endif // ROOT_Minuit2_SimplexParameters