GeometricAnnealing.h

Go to the documentation of this file.

#ifndef GeometricAnnealing_H
#define GeometricAnnealing_H

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

class GeometricAnnealing : public AnnealingSchedule {

public:

  GeometricAnnealing( const double cutoff=3.0, const double T=256.0,
     const double annealing_ratio=0.25 );

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

  double phi ( double chi2 ) const;

  double weight ( double chi2 ) const;

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

  bool isAnnealed() const;

  void debug() const;

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

private:
  double theT0;
  double theT;
  double theChi2cut;
  double theRatio;

};

#endif