#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	(	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 ( )

overridevirtual

Implements AnnealingSchedule.

Definition at line 12 of file GeometricAnnealing.cc.

References theRatio, and theT.

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

GeometricAnnealing* GeometricAnnealing::clone ( void ) const

inlineoverridevirtual

Implements AnnealingSchedule.

Definition at line 45 of file GeometricAnnealing.h.

References GeometricAnnealing().

   {
     return new GeometricAnnealing ( * this );
   };

double GeometricAnnealing::currentTemp ( ) const

overridevirtual

Implements AnnealingSchedule.

Definition at line 45 of file GeometricAnnealing.cc.

References theT.

 {
   return theT;
 }

double GeometricAnnealing::cutoff ( ) const

overridevirtual

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

overridevirtual

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

overridevirtual

Implements AnnealingSchedule.

Definition at line 50 of file GeometricAnnealing.cc.

References theT0.

 {
   return theT0;
 }

bool GeometricAnnealing::isAnnealed ( ) const

overridevirtual

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

overridevirtual

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__(), ntupleDataFormat.Track::phiPull(), and weight().

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

void GeometricAnnealing::resetAnnealing ( )

overridevirtual

Implements AnnealingSchedule.

Definition at line 29 of file GeometricAnnealing.cc.

References theT, and theT0.

 {
   theT=theT0;
 }

double GeometricAnnealing::weight ( double chi2 ) const

overridevirtual

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