#include <CBAlgoPolyShortestDist.h>

Inheritance diagram for cmtool::CBAlgoPolyShortestDist:

Public Member Functions
	CBAlgoPolyShortestDist ()
	Default constructor. More...

virtual	~CBAlgoPolyShortestDist ()
	Default destructor. More...

virtual void	EventBegin (const std::vector< cluster::ClusterParamsAlg > &clusters)

virtual bool	Bool (const ::cluster::ClusterParamsAlg &cluster1, const ::cluster::ClusterParamsAlg &cluster2)

virtual void	Report ()

virtual void	Reset ()
	Function to reset the algorithm instance ... maybe implemented via child class. More...

void	SetMinNumHits (size_t nhits)

void	SetMaxNumHits (int nhits)

void	SetMinDistSquared (double dist)

void	SetDebug (bool flag)

Public Member Functions inherited from cmtool::CBoolAlgoBase
	CBoolAlgoBase ()
	Default constructor. More...

virtual	~CBoolAlgoBase ()
	Default destructor. More...

Public Member Functions inherited from cmtool::CMAlgoBase
	CMAlgoBase ()

virtual	~CMAlgoBase ()=default

virtual void	EventEnd ()

virtual void	IterationBegin (const std::vector< cluster::ClusterParamsAlg > &)

virtual void	IterationEnd ()

void	SetAnaFile (TFile *fout)
	Setter function for an output plot TFile pointer. More...

virtual void	SetVerbose (bool doit=true)
	Setter function for verbosity. More...

Private Attributes
size_t	_min_hits

size_t	_max_hits

double	_dist_sqrd_cut

bool	_debug

double	tmp_min_dist

Additional Inherited Members
Protected Attributes inherited from cmtool::CMAlgoBase
TFile *	_fout
	TFile pointer to an output file. More...

bool	_verbose
	Boolean to choose verbose mode. Turned on if CMergeManager/CMatchManager's verbosity level is >= kPerMerging. More...

Detailed Description

User implementation for CBoolAlgoBase class doxygen documentation!

Definition at line 28 of file CBAlgoPolyShortestDist.h.

Constructor & Destructor Documentation

cmtool::CBAlgoPolyShortestDist::CBAlgoPolyShortestDist ( )

Default constructor.

Definition at line 8 of file CBAlgoPolyShortestDist.cxx.

                                                  : CBoolAlgoBase()
   //-------------------------------------------------------
   {
     SetMinDistSquared(1.);
     SetMinNumHits(0);
     SetMaxNumHits(99999);
     SetDebug(false);
   }

virtual cmtool::CBAlgoPolyShortestDist::~CBAlgoPolyShortestDist ( )

inlinevirtual

Default destructor.

Definition at line 36 of file CBAlgoPolyShortestDist.h.

36 {};

Member Function Documentation

bool cmtool::CBAlgoPolyShortestDist::Bool	(	const ::cluster::ClusterParamsAlg &	cluster1,
		const ::cluster::ClusterParamsAlg &	cluster2
	)

virtual

Optional function: called at the end of event ... after the last merging iteration is over. Optional function: called at the beggining of each iteration over all pairs of clusters. This provides all clusters' information in case the algorithm need them. Note this is called per iteration which may be more than once per event. Optional function: called at the end of each iteration over all pairs of clusters. Core function: given the CPAN input, return whether a cluster should be merged or not.

Reimplemented from cmtool::CBoolAlgoBase.

Definition at line 54 of file CBAlgoPolyShortestDist.cxx.

   {
     if( (cluster1.GetHitVector().size() < _min_hits) ||
         (cluster2.GetHitVector().size() < _min_hits) )
       return false;
 
     if( (cluster1.GetHitVector().size() > _max_hits) ||
         (cluster2.GetHitVector().size() > _max_hits) )
       return false;
 
     //if either has < 3 sides do not merge!
     if ( (cluster1.GetParams().PolyObject.Size() < 2) or
          (cluster2.GetParams().PolyObject.Size() < 2) ){
       return false;
     }
 
     //loop over the points on the first polygon and calculate
     //distance to each point on the second polygon
     //if any two points are close enough to each other,
     //merge the two clusters
 
     unsigned int npoints1 = cluster1.GetParams().PolyObject.Size();
     unsigned int npoints2 = cluster2.GetParams().PolyObject.Size();
     //loop over points on first polygon
     for(unsigned int i = 0; i < npoints1; ++i){
       float pt1w = cluster1.GetParams().PolyObject.Point(i).first;
       float pt1t = cluster1.GetParams().PolyObject.Point(i).second;
       //loop over points on second polygon
       for(unsigned int j = 0; j < npoints2; ++j){
         float pt2w = cluster2.GetParams().PolyObject.Point(j).first;
         float pt2t = cluster2.GetParams().PolyObject.Point(j).second;
         double distsqrd = pow(pt2w-pt1w,2)+pow(pt2t-pt1t,2);
 
         if(distsqrd < tmp_min_dist) tmp_min_dist = distsqrd;
 
         if(_debug){
           std::cout<<"two polygon points dist2 is "<<distsqrd<<std::endl;
           std::cout<<"minimum dist was "<<tmp_min_dist<<std::endl;
         }
         if(distsqrd<_dist_sqrd_cut)
           return true;
       }
 
     }
 
     return false;
   }

void cmtool::CBAlgoPolyShortestDist::EventBegin ( const std::vector< cluster::ClusterParamsAlg > & clusters )

virtual

Optional function: called at the beginning of 1st iteration. This is called per event.

Reimplemented from cmtool::CMAlgoBase.

Definition at line 25 of file CBAlgoPolyShortestDist.cxx.

   {
     if(clusters.size())
       tmp_min_dist = 99999;
   }

void cmtool::CBAlgoPolyShortestDist::Report ( )

virtual

Optional function: called after each Merge() function call by CMergeManager IFF CMergeManager is run with verbosity level kPerMerging. Maybe useful for debugging.

Reimplemented from cmtool::CMAlgoBase.

Definition at line 105 of file CBAlgoPolyShortestDist.cxx.