DeterministicAnnealing Class Reference

#include <DeterministicAnnealing.h>

Inheritance diagram for DeterministicAnnealing:

Classes
class	A

Public Member Functions
void	anneal () override
DeterministicAnnealing *	clone () const override
double	currentTemp () const override
double	cutoff () const override
void	debug () const override
	DeterministicAnnealing (float cutoff=3.0)
	DeterministicAnnealing (const std::vector< float > &sched, float cutoff=3.0)
double	initialTemp () const override
bool	isAnnealed () const override
double	phi (double chi2) const override
void	resetAnnealing () override
double	weight (double chi2) const override
Public Member Functions inherited from AnnealingSchedule
virtual	~AnnealingSchedule ()

Private Attributes
double	theChi2cut
unsigned int	theIndex
bool	theIsAnnealed
std::vector< float >	theTemperatures

Detailed Description

Definition at line 7 of file DeterministicAnnealing.h.

Constructor & Destructor Documentation

DeterministicAnnealing::DeterministicAnnealing

(

float

cutoff = 3.0

)

Definition at line 10 of file DeterministicAnnealing.cc.

    : theTemperatures({256, 64, 16, 4, 2, 1}), theIndex(0), theChi2cut(cutoff * cutoff), theIsAnnealed(false) {}

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

Definition at line 13 of file DeterministicAnnealing.cc.

    : theTemperatures(sched), theIndex(0), theChi2cut(cutoff * cutoff), theIsAnnealed(false) {}

Member Function Documentation

void DeterministicAnnealing::anneal ( )

overridevirtual

Implements AnnealingSchedule.

Definition at line 16 of file DeterministicAnnealing.cc.

References theIndex, theIsAnnealed, and theTemperatures.

                                    {
  if (theIndex < (theTemperatures.size() - 1)) {
    theIndex++;
  } else {
    theIsAnnealed = true;
  };
}

DeterministicAnnealing* DeterministicAnnealing::clone ( void  ) const

inlineoverridevirtual

Implements AnnealingSchedule.

Definition at line 48 of file DeterministicAnnealing.h.

{ return new DeterministicAnnealing(*this); };

double DeterministicAnnealing::currentTemp ( ) const

overridevirtual

Implements AnnealingSchedule.

Definition at line 50 of file DeterministicAnnealing.cc.

References theIndex, and theTemperatures.

{ return theTemperatures[theIndex]; }

double DeterministicAnnealing::cutoff ( ) const

overridevirtual

Returns phi(chi2) / ( phi(cutoff^2) + sum_i { phi(chi2s[i]) } )

Implements AnnealingSchedule.

Definition at line 48 of file DeterministicAnnealing.cc.

References mathSSE::sqrt(), and theChi2cut.

{ return sqrt(theChi2cut); }

void DeterministicAnnealing::debug ( ) const

overridevirtual

Implements AnnealingSchedule.

Definition at line 56 of file DeterministicAnnealing.cc.

References gather_cfg::cout, mps_fire::i, and theTemperatures.

                                         {
  cout << "[DeterministicAnnealing] schedule=";
  for (vector<float>::const_iterator i = theTemperatures.begin(); i != theTemperatures.end(); ++i) {
    cout << *i << " ";
  };
  cout << endl;
}

double DeterministicAnnealing::initialTemp ( ) const

overridevirtual

Implements AnnealingSchedule.

Definition at line 52 of file DeterministicAnnealing.cc.

References theTemperatures.

{ return theTemperatures[0]; }

bool DeterministicAnnealing::isAnnealed ( ) const

overridevirtual

is it annealed yet?

Implements AnnealingSchedule.

Definition at line 54 of file DeterministicAnnealing.cc.

References theIsAnnealed.

{ return theIsAnnealed; }

double DeterministicAnnealing::phi ( double  chi2 ) const

inlineoverridevirtual

phi ( chi2 ) = e^( -.5*chi2 / T )

Implements AnnealingSchedule.

Definition at line 46 of file DeterministicAnnealing.cc.

References funct::exp(), theIndex, and theTemperatures.

Referenced by Particle.Particle::__str__(), ntupleDataFormat.Track::phiPull(), and weight().

{ return exp(-.5 * chi2 / theTemperatures[theIndex]); }

void DeterministicAnnealing::resetAnnealing ( )

overridevirtual

Implements AnnealingSchedule.

Definition at line 41 of file DeterministicAnnealing.cc.

References theIndex, and theIsAnnealed.

                                            {
  theIndex = 0;
  theIsAnnealed = false;
}

double DeterministicAnnealing::weight ( double  chi2 ) const

overridevirtual

Returns phi(chi2) / ( phi(cutoff^2) + phi(chi2) ),

Implements AnnealingSchedule.

Definition at line 24 of file DeterministicAnnealing.cc.

References edm::isNotFinite(), phi(), and theChi2cut.

                                                       {
  long double mphi = phi(chi2);
  /*
  if ( mphi < std::numeric_limits<double>::epsilon() ) return 0.;
  return 1. / ( 1. + phi ( theChi2cut * theChi2cut ) / mphi );
  */
  // return mphi / ( mphi + phi ( theChi2cut ) );
  long double newtmp = mphi / (mphi + phi(theChi2cut));
  if (edm::isNotFinite(newtmp)) {
    if (chi2 < theChi2cut)
      newtmp = 1.;
    else
      newtmp = 0.;
  }
  return newtmp;
}

Member Data Documentation

double DeterministicAnnealing::theChi2cut

private

Definition at line 53 of file DeterministicAnnealing.h.

Referenced by cutoff(), and weight().

unsigned int DeterministicAnnealing::theIndex

private

Definition at line 52 of file DeterministicAnnealing.h.

Referenced by anneal(), currentTemp(), phi(), and resetAnnealing().

bool DeterministicAnnealing::theIsAnnealed

private

Definition at line 54 of file DeterministicAnnealing.h.

Referenced by anneal(), isAnnealed(), and resetAnnealing().

std::vector<float> DeterministicAnnealing::theTemperatures

private

Definition at line 48 of file DeterministicAnnealing.h.

Referenced by anneal(), currentTemp(), debug(), initialTemp(), and phi().