fit::RootMinuitCommands< Function > Class Template Reference

#include <RootMinuitCommands.h>

Public Types
typedef RootMinuitCommand	command
typedef RootMinuit< Function >	minuit

Public Member Functions
void	add (RootMinuit< Function > &minuit, funct::Parameter &p) const
double	err (const std::string &name)
bool	fixed (const std::string &name)
void	init (const char *fileName)
double	max (const std::string &name)
double	min (const std::string &name)
double	par (const std::string &name)
	RootMinuitCommands (bool verbose=true)
	RootMinuitCommands (const char *fileName, bool verbose=true)
void	run (RootMinuit< Function > &minuit) const

Private Types
typedef boost::tokenizer< boost::char_separator< char > >	tokenizer

Private Member Functions
std::string	errorHeader () const
std::string	nextToken (typename tokenizer::iterator &i, const typename tokenizer::iterator &end) const
const parameter_t &	parameter (const std::string &name) const
double	string2double (const std::string &str) const

Private Attributes
std::vector< command >	commands_
unsigned int	lineNumber_
std::map< std::string, size_t >	parIndices_
parameterVector_t	pars_
bool	verbose_

Detailed Description

template<class Function>
class fit::RootMinuitCommands< Function >

Definition at line 50 of file RootMinuitCommands.h.

Member Typedef Documentation

template<class Function >

typedef RootMinuitCommand fit::RootMinuitCommands< Function >::command

Definition at line 53 of file RootMinuitCommands.h.

template<class Function >

typedef RootMinuit<Function> fit::RootMinuitCommands< Function >::minuit

Definition at line 52 of file RootMinuitCommands.h.

template<class Function >

typedef boost::tokenizer<boost::char_separator<char> > fit::RootMinuitCommands< Function >::tokenizer

private

Definition at line 72 of file RootMinuitCommands.h.

Constructor & Destructor Documentation

template<class Function >

fit::RootMinuitCommands< Function >::RootMinuitCommands ( bool  verbose = true )

inline

Definition at line 54 of file RootMinuitCommands.h.

: verbose_(verbose) {}

template<class Function >

fit::RootMinuitCommands< Function >::RootMinuitCommands ( const char *  fileName, bool  verbose = true  )

inline

Definition at line 55 of file RootMinuitCommands.h.

References MillePedeFileConverter_cfg::fileName, and init.

: verbose_(verbose) { init(fileName); }

Member Function Documentation

template<class Function >

void fit::RootMinuitCommands< Function >::add ( RootMinuit< Function > &  minuit, funct::Parameter &  p  ) const

inline

Definition at line 62 of file RootMinuitCommands.h.

References fit::RootMinuit< Function >::addParameter(), fit::parameter_t::err, fit::parameter_t::fixed, fit::RootMinuit< Function >::fixParameter(), fit::parameter_t::max, fit::parameter_t::min, funct::Parameter::name(), fit::RootMinuitCommand::name, pfMETsysShiftCorrections_cfi::parameter, writedatasetfile::run, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by counter.Counter::register().

                                                                    {
      const std::string& name = p.name();
      const parameter_t& par = parameter(name);
      minuit.addParameter(p, par.err, par.min, par.max);
      if (par.fixed)
        minuit.fixParameter(name);
    }

template<class Function >

double fit::RootMinuitCommands< Function >::err ( const std::string &  name )

inline

Definition at line 58 of file RootMinuitCommands.h.

References pfMETsysShiftCorrections_cfi::parameter.

{ return parameter(name).err; }

template<class Function >

std::string fit::RootMinuitCommands< Function >::errorHeader ( ) const

inlineprivate

Definition at line 94 of file RootMinuitCommands.h.

References MillePedeFileConverter_cfg::out.

                                  {
      std::ostringstream out;
      out << "RootMinuitCommands config. error, line " << lineNumber_ << ": ";
      return out.str();
    }

template<class Function >

bool fit::RootMinuitCommands< Function >::fixed ( const std::string &  name )

inline

Definition at line 61 of file RootMinuitCommands.h.

References pfMETsysShiftCorrections_cfi::parameter.

{ return parameter(name).fixed; }

template<typename Function >

void fit::RootMinuitCommands< Function >::init

(

const char *

fileName

)

Definition at line 108 of file RootMinuitCommands.h.

References runTheMatrix::cmssw_base, tkAlTrackRefitSequence_cff::commands, edm::errors::Configuration, gather_cfg::cout, ztail::d, fit::RootMinuitCommand::doubleArgs, MillePedeFileConverter_cfg::e, fit::parameter_t::err, Exception, f, FrontierConditions_GlobalTag_cff::file, fit::parameter_t::fixed, alignBH_cfg::fixed, mps_fire::i, dqmiolumiharvest::j, join(), kFix, kMigrad, kMinimize, kParameter, kPrintAll, kRelease, kSet, dqmdumpme::last, mps_splice::line, fit::parameter_t::max, fit::parameter_t::min, eostools::move(), fit::RootMinuitCommand::name, fit::RootMinuitCommand::print(), alignCSCRings::s, fit::RootMinuitCommand::stringArgs, and fit::parameter_t::val.

                                                              {
    using namespace std;
    string cmssw_release_base = std::getenv("CMSSW_RELEASE_BASE");
    string cmssw_base = std::getenv("CMSSW_BASE");
    vector<string> directories;
    directories.reserve(3);
    directories.emplace_back(".");
    if (!cmssw_release_base.empty()) {
      directories.emplace_back(cmssw_release_base + "/src");
    }
    if (!cmssw_base.empty()) {
      directories.emplace_back(cmssw_base + "/src");
    }
    ifstream file;
    for (auto const& d : directories) {
      std::ifstream f{d + "/" + fileName};
      if (f.good()) {
        file = std::move(f);
        break;
      }
    }
    if (!file.is_open()) {
      throw edm::Exception(edm::errors::Configuration)
          << "RootMinuitCommands: can't open file: " << fileName
          << " in path: " << boost::algorithm::join(directories, ":") << "\n";
    }
    if (verbose_)
      cout << ">>> configuration file: " << fileName << endl;
    string line;
    lineNumber_ = 0;
    bool commands = false;
    while (getline(file, line)) {
      ++lineNumber_;
      if (line.size() == 0)
        continue;
      char last = *line.rbegin();
      if (!(last >= '0' && last <= 'z'))
        line.erase(line.end() - 1);
      boost::char_separator<char> sep(" ");
      tokenizer tokens(line, sep);
      tokenizer::iterator i = tokens.begin(), e = tokens.end();
      if (tokens.begin() == tokens.end())
        continue;
      if (*(i->begin()) != '#') {
        if (*i == kParameter) {
          if (commands)
            throw edm::Exception(edm::errors::Configuration)
                << errorHeader() << "please, declare all parameter before all other minuit commands.\n";
          string name = nextToken(i, e);
          parameter_t par;
          par.val = string2double(nextToken(i, e));
          par.err = string2double(nextToken(i, e));
          par.min = string2double(nextToken(i, e));
          par.max = string2double(nextToken(i, e));
          tokenizer::iterator j = i;
          ++j;
          if (j != e) {
            string fixed = nextToken(i, e);
            if (fixed == "fixed")
              par.fixed = true;
            else if (fixed == "free")
              par.fixed = false;
            else
              throw edm::Exception(edm::errors::Configuration)
                  << errorHeader() << "fix parameter option unknown: " << *i << "\n"
                  << "valid options are: fixed, free.\n";
          } else {
            par.fixed = false;
          }
          pars_.push_back(std::make_pair(name, par));
          size_t s = parIndices_.size();
          parIndices_[name] = s;
          if (verbose_)
            cout << ">>> " << kParameter << " " << name << " " << par.val << " [" << par.min << ", " << par.max << "],"
                 << " err: " << par.err << endl;
        } else if (*i == kFix || *i == kRelease) {
          commands = true;
          command com;
          com.name = *i;
          string arg = nextToken(i, e);
          com.stringArgs.push_back(arg);
          commands_.push_back(com);
          if (verbose_) {
            cout << ">>> ";
            com.print(cout);
            cout << endl;
          }
        } else if (*i == kSet) {
          commands = true;
          command com;
          com.name = *i;
          string arg = nextToken(i, e);
          com.stringArgs.push_back(arg);
          com.doubleArgs.push_back(string2double(nextToken(i, e)));
          commands_.push_back(com);
          if (verbose_) {
            cout << ">>> ";
            com.print(cout);
            cout << endl;
          }
        } else if (*i == kMinimize || *i == kMigrad || *i == kPrintAll) {
          commands = true;
          command com;
          com.name = *i;
          commands_.push_back(com);
          if (verbose_) {
            cout << ">>> ";
            com.print(cout);
            cout << endl;
          }
        } else {
          throw edm::Exception(edm::errors::Configuration) << errorHeader() << "unkonwn command:: " << *i << "\n";
        }
      }
    }
    if (verbose_)
      cout << ">>> end configuration" << endl;
  }

template<class Function >

double fit::RootMinuitCommands< Function >::max ( const std::string &  name )

inline

Definition at line 60 of file RootMinuitCommands.h.

References pfMETsysShiftCorrections_cfi::parameter.

{ return parameter(name).max; }

template<class Function >

double fit::RootMinuitCommands< Function >::min ( const std::string &  name )

inline

Definition at line 59 of file RootMinuitCommands.h.

References pfMETsysShiftCorrections_cfi::parameter.

{ return parameter(name).min; }

template<class Function >

std::string fit::RootMinuitCommands< Function >::nextToken ( typename tokenizer::iterator &  i, const typename tokenizer::iterator &  end  ) const

inlineprivate

Definition at line 99 of file RootMinuitCommands.h.

References edm::errors::Configuration, Exception, and mps_fire::i.

                                                                                                    {
      ++i;
      if (i == end)
        throw edm::Exception(edm::errors::Configuration) << errorHeader() << "missing parameter\n";
      return *i;
    }

template<class Function >

double fit::RootMinuitCommands< Function >::par ( const std::string &  name )

inline

Definition at line 57 of file RootMinuitCommands.h.

References pfMETsysShiftCorrections_cfi::parameter.

{ return parameter(name).val; }

template<class Function >

const parameter_t& fit::RootMinuitCommands< Function >::parameter ( const std::string &  name ) const

inlineprivate

Definition at line 88 of file RootMinuitCommands.h.

References edm::errors::Configuration, and AlCaHLTBitMon_ParallelJobs::p.

                                                              {
      typename std::map<std::string, size_t>::const_iterator p = parIndices_.find(name);
      if (p == parIndices_.end())
        throw edm::Exception(edm::errors::Configuration) << "RootMinuit: can't find parameter " << name << "\n";
      return pars_[p->second].second;
    }

template<typename Function >

void fit::RootMinuitCommands< Function >::run

(

RootMinuit< Function > &

minuit

)

const

Definition at line 228 of file RootMinuitCommands.h.

References HltBtagPostValidation_cff::c, gather_cfg::cout, end, fit::RootMinuit< Function >::fixParameter(), kFix, kMigrad, kMinimize, kPrintAll, kRelease, kSet, fit::RootMinuit< Function >::migrad(), fit::RootMinuit< Function >::minimize(), fit::RootMinuit< Function >::printFitResults(), fit::RootMinuit< Function >::releaseParameter(), and fit::RootMinuit< Function >::setParameter().

                                                                           {
    using namespace std;
    typename vector<command>::const_iterator c = commands_.begin(), end = commands_.end();
    for (; c != end; ++c) {
      if (verbose_) {
        cout << ">>> minuit command: ";
        c->print(cout);
        cout << endl;
      }
      if (c->name == kMinimize)
        minuit.minimize();
      else if (c->name == kMigrad)
        minuit.migrad();
      else if (c->name == kPrintAll)
        minuit.printFitResults();
      else if (c->name == kFix)
        minuit.fixParameter(c->stringArgs[0]);
      else if (c->name == kRelease)
        minuit.releaseParameter(c->stringArgs[0]);
      else if (c->name == kSet)
        minuit.setParameter(c->stringArgs[0], c->doubleArgs[0]);
    }
  }

template<class Function >

double fit::RootMinuitCommands< Function >::string2double ( const std::string &  str ) const

inlineprivate

Definition at line 78 of file RootMinuitCommands.h.

References begin, edm::errors::Configuration, end, Exception, alignCSCRings::s, and heppy_batch::val.

                                                     {
      const char* begin = str.c_str();
      char* end;
      double val = strtod(begin, &end);
      size_t s = end - begin;
      if (s < str.size()) {
        throw edm::Exception(edm::errors::Configuration) << "RootMinuitCommands: invalid double value: " << str << "\n";
      }
      return val;
    }

Member Data Documentation

template<class Function >

std::vector<command> fit::RootMinuitCommands< Function >::commands_

private

Definition at line 77 of file RootMinuitCommands.h.

template<class Function >

unsigned int fit::RootMinuitCommands< Function >::lineNumber_

private

Definition at line 74 of file RootMinuitCommands.h.

template<class Function >

std::map<std::string, size_t> fit::RootMinuitCommands< Function >::parIndices_

private

Definition at line 76 of file RootMinuitCommands.h.

template<class Function >

parameterVector_t fit::RootMinuitCommands< Function >::pars_

private

Definition at line 75 of file RootMinuitCommands.h.

template<class Function >

bool fit::RootMinuitCommands< Function >::verbose_

private

Definition at line 73 of file RootMinuitCommands.h.