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.

References cutoff(), theChi2cut, theIndex, and theIsAnnealed.

Referenced by clone().

                                                              :
  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 16 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 21 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 51 of file DeterministicAnnealing.h.

References DeterministicAnnealing().

  {
    return new DeterministicAnnealing ( * this );
  };

double DeterministicAnnealing::currentTemp ( ) const

overridevirtual

Implements AnnealingSchedule.

Definition at line 64 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 59 of file DeterministicAnnealing.cc.

References mathSSE::sqrt(), and theChi2cut.

Referenced by DeterministicAnnealing().

{
  return sqrt(theChi2cut);
}

void DeterministicAnnealing::debug ( ) const

overridevirtual

Implements AnnealingSchedule.

Definition at line 79 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 69 of file DeterministicAnnealing.cc.

References theTemperatures.

{
  return theTemperatures[0];
}

bool DeterministicAnnealing::isAnnealed ( ) const

overridevirtual

is it annealed yet?

Implements AnnealingSchedule.

Definition at line 74 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 54 of file DeterministicAnnealing.cc.

References JetChargeProducer_cfi::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 48 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 31 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 59 of file DeterministicAnnealing.h.

Referenced by cutoff(), DeterministicAnnealing(), and weight().

unsigned int DeterministicAnnealing::theIndex

private

Definition at line 58 of file DeterministicAnnealing.h.

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

bool DeterministicAnnealing::theIsAnnealed

private

Definition at line 60 of file DeterministicAnnealing.h.

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

std::vector<float> DeterministicAnnealing::theTemperatures

private

Definition at line 54 of file DeterministicAnnealing.h.

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