#include <GeometricAnnealing.h>

Inheritance diagram for GeometricAnnealing:

Classes
class	A

Public Member Functions
void	anneal ()

GeometricAnnealing *	clone () const

double	currentTemp () const

double	cutoff () const

void	debug () const

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

double	initialTemp () const

bool	isAnnealed () const

double	phi (double chi2) const

void	resetAnnealing ()

double	weight (double chi2) const

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	(	const double	cutoff = `3.0`,
		const double	T = `256.0`,
		const double	annealing_ratio = `0.25`
	)

Definition at line 7 of file GeometricAnnealing.cc.

Referenced by clone().

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

Member Function Documentation

void GeometricAnnealing::anneal ( )

virtual

Implements AnnealingSchedule.

Definition at line 12 of file GeometricAnnealing.cc.

References theRatio, and theT.

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

GeometricAnnealing* GeometricAnnealing::clone ( void ) const

inlinevirtual

Implements AnnealingSchedule.

Definition at line 45 of file GeometricAnnealing.h.

References GeometricAnnealing().

   {
     return new GeometricAnnealing ( * this );
   };

double GeometricAnnealing::currentTemp ( ) const

virtual

Implements AnnealingSchedule.

Definition at line 45 of file GeometricAnnealing.cc.

References theT.

 {
   return theT;
 }

double GeometricAnnealing::cutoff ( ) const

virtual

Implements AnnealingSchedule.

Definition at line 39 of file GeometricAnnealing.cc.

References mathSSE::sqrt(), and theChi2cut.

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

void GeometricAnnealing::debug ( ) const

virtual

Implements AnnealingSchedule.

Definition at line 60 of file GeometricAnnealing.cc.

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

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

double GeometricAnnealing::initialTemp ( ) const

virtual

Implements AnnealingSchedule.

Definition at line 50 of file GeometricAnnealing.cc.

References theT0.

 {
   return theT0;
 }

bool GeometricAnnealing::isAnnealed ( ) const

virtual

is it annealed yet?

Implements AnnealingSchedule.

Definition at line 55 of file GeometricAnnealing.cc.

References theT.

 {
   return ( theT < 1.02 );
 }

double GeometricAnnealing::phi ( double chi2 ) const

virtual

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

Implements AnnealingSchedule.

Definition at line 34 of file GeometricAnnealing.cc.

References JetChargeProducer_cfi::exp, and theT.

Referenced by Particle.Particle::__str__(), and weight().

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

void GeometricAnnealing::resetAnnealing ( )

virtual

Implements AnnealingSchedule.

Definition at line 29 of file GeometricAnnealing.cc.

References theT, and theT0.

 {
   theT=theT0;
 }

double GeometricAnnealing::weight ( double chi2 ) const

virtual

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

Implements AnnealingSchedule.

Definition at line 17 of file GeometricAnnealing.cc.

References 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;
 }