GeometricAnnealing Class Reference

#include <GeometricAnnealing.h>

Inheritance diagram for GeometricAnnealing:

Classes
class	A

Public Member Functions
void	anneal () override
GeometricAnnealing *	clone () const override
double	currentTemp () const override
double	cutoff () const override
void	debug () const override
	GeometricAnnealing (const double cutoff=3.0, const double T=256.0, const double annealing_ratio=0.25)
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
double	theRatio
double	theT
double	theT0

Detailed Description

Definition at line 7 of file GeometricAnnealing.h.

Constructor & Destructor Documentation

GeometricAnnealing()

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

Definition at line 7 of file GeometricAnnealing.cc.

    : theT0(T), theT(T), theChi2cut(cutoff * cutoff), theRatio(ratio) {}

Member Function Documentation

anneal()

void GeometricAnnealing::anneal ( )

overridevirtual

Implements AnnealingSchedule.

Definition at line 10 of file GeometricAnnealing.cc.

References theRatio, and theT.

{ theT = 1 + (theT - 1) * theRatio; }

clone()

GeometricAnnealing* GeometricAnnealing::clone ( void  ) const

inlineoverridevirtual

Implements AnnealingSchedule.

Definition at line 42 of file GeometricAnnealing.h.

{ return new GeometricAnnealing(*this); };

currentTemp()

double GeometricAnnealing::currentTemp ( ) const

overridevirtual

Implements AnnealingSchedule.

Definition at line 33 of file GeometricAnnealing.cc.

References theT.

{ return theT; }

cutoff()

double GeometricAnnealing::cutoff ( ) const

overridevirtual

Implements AnnealingSchedule.

Definition at line 28 of file GeometricAnnealing.cc.

References mathSSE::sqrt(), and theChi2cut.

                                        {
  // std::cout << "[GeometricAnnealing] cutoff called!" << std::endl;
  return sqrt(theChi2cut);
}

debug()

void GeometricAnnealing::debug ( ) const

overridevirtual

Implements AnnealingSchedule.

Definition at line 39 of file GeometricAnnealing.cc.

References gather_cfg::cout, theChi2cut, theRatio, and theT0.

                                     {
  std::cout << "[GeometricAnnealing] chi2_cut=" << theChi2cut << ", Tini=" << theT0 << ", ratio=" << theRatio
            << std::endl;
}

initialTemp()

double GeometricAnnealing::initialTemp ( ) const

overridevirtual

Implements AnnealingSchedule.

Definition at line 35 of file GeometricAnnealing.cc.

References theT0.

{ return theT0; }

isAnnealed()

bool GeometricAnnealing::isAnnealed ( ) const

overridevirtual

is it annealed yet?

Implements AnnealingSchedule.

Definition at line 37 of file GeometricAnnealing.cc.

References theT.

{ return (theT < 1.02); }

phi()

double GeometricAnnealing::phi ( double  chi2 ) const

overridevirtual

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

Implements AnnealingSchedule.

Definition at line 26 of file GeometricAnnealing.cc.

References nano_mu_local_reco_cff::chi2, JetChargeProducer_cfi::exp, and theT.

Referenced by Particle.Particle::__str__(), datamodel.Object::DeltaR(), datamodel.Object::p4(), ntupleDataFormat.Track::phiPull(), and weight().

{ return exp(-.5 * chi2 / theT); }

resetAnnealing()

void GeometricAnnealing::resetAnnealing ( )

overridevirtual

Implements AnnealingSchedule.

Definition at line 24 of file GeometricAnnealing.cc.

References theT, and theT0.

{ theT = theT0; }

weight()

double GeometricAnnealing::weight ( double  chi2 ) const

overridevirtual

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

Implements AnnealingSchedule.

Definition at line 12 of file GeometricAnnealing.cc.

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

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

Member Data Documentation

theChi2cut

double GeometricAnnealing::theChi2cut

private

Definition at line 47 of file GeometricAnnealing.h.

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

theRatio

double GeometricAnnealing::theRatio

private

Definition at line 48 of file GeometricAnnealing.h.

Referenced by anneal(), and debug().

theT

double GeometricAnnealing::theT

private

Definition at line 46 of file GeometricAnnealing.h.

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

theT0

double GeometricAnnealing::theT0

private

Definition at line 42 of file GeometricAnnealing.h.

Referenced by debug(), initialTemp(), and resetAnnealing().