DeterministicAnnealing.h

Go to the documentation of this file.

#ifndef DeterministicAnnealing_H
#define DeterministicAnnealing_H

#include "RecoVertex/VertexTools/interface/AnnealingSchedule.h"
#include <vector>

class DeterministicAnnealing : public AnnealingSchedule {

public:

  DeterministicAnnealing( float cutoff = 3.0 );
  DeterministicAnnealing( const std::vector < float > & sched, float cutoff = 3.0 );

  void anneal(); //< One annealing step. theT *= theRatio.
  void resetAnnealing(); //< theT = theT0.

  double phi ( double chi2 ) const;

  double weight ( double chi2 ) const;
  
  bool isAnnealed() const;

  void debug() const;

  // double weight ( double chi2, const vector < double > & chi2s ) const;


  double cutoff() const;
  double currentTemp() const;
  double initialTemp() const;

  DeterministicAnnealing * clone() const
  {
    return new DeterministicAnnealing ( * this );
  };

private:
  std::vector<float> theTemperatures;
  unsigned int theIndex;
  double theChi2cut;
  bool theIsAnnealed;
};

#endif