ConfigurableAnnealing.cc

Go to the documentation of this file.

#include "RecoVertex/ConfigurableVertexReco/interface/ConfigurableAnnealing.h"
#include "RecoVertex/VertexTools/interface/GeometricAnnealing.h"
#include "RecoVertex/VertexTools/interface/DeterministicAnnealing.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <string>
 
using namespace std;
 
ConfigurableAnnealing::ConfigurableAnnealing(const edm::ParameterSet& m) : theImpl(nullptr) {
  string type = m.getParameter<string>("annealing");
  // edm::LogWarning("ConfigurableAnnealing") << "below one code ist still here.";
  if (type == "below") {
    edm::LogError("ConfigurableAnnealing") << "below one annealing employed!";
    vector<float> sched;
    double final = m.getParameter<double>("Tfinal");
    sched.push_back(256.);
    sched.push_back(64.);
    sched.push_back(16.);
    sched.push_back(4.);
    sched.push_back(1.);
    sched.push_back(final);
    theImpl = new DeterministicAnnealing(sched, m.getParameter<double>("sigmacut"));
  } else if (type == "geom") {
    theImpl = new GeometricAnnealing(
        m.getParameter<double>("sigmacut"), m.getParameter<double>("Tini"), m.getParameter<double>("ratio"));
  } else {
    edm::LogError("ConfigurableAnnealing") << "annealing type " << type << " is not known.";
    exit(-1);
  }
}
 
ConfigurableAnnealing::ConfigurableAnnealing(const ConfigurableAnnealing& o) : theImpl(o.theImpl->clone()) {}
 
ConfigurableAnnealing* ConfigurableAnnealing::clone() const { return new ConfigurableAnnealing(*this); }
 
ConfigurableAnnealing::~ConfigurableAnnealing() { delete theImpl; }
 
void ConfigurableAnnealing::debug() const { theImpl->debug(); }
 
void ConfigurableAnnealing::anneal() { theImpl->anneal(); }
 
double ConfigurableAnnealing::weight(double chi2) const { return theImpl->weight(chi2); }
 
void ConfigurableAnnealing::resetAnnealing() { theImpl->resetAnnealing(); }
 
inline double ConfigurableAnnealing::phi(double chi2) const { return theImpl->phi(chi2); }
 
double ConfigurableAnnealing::cutoff() const { return theImpl->cutoff(); }
 
double ConfigurableAnnealing::currentTemp() const { return theImpl->currentTemp(); }
 
double ConfigurableAnnealing::initialTemp() const { return theImpl->initialTemp(); }
 
bool ConfigurableAnnealing::isAnnealed() const { return theImpl->isAnnealed(); }