#include <CBAlgoTrackSeparate.h>

Inheritance diagram for cmtool::CBAlgoTrackSeparate:

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

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

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

void	SetVerbose (bool on)
	Setter function for verbosity. More...

void	SetDebug (bool on)

void	SetMinNumHits (size_t n)

void	SetMinAngleDiff (float anglesep)

void	SetMaxOpeningAngle (float maxangle)

void	SetMinLength (float minlength)

void	SetMinPolyHitDensity (float mindensity)

void	SetMaxWidth (float maxwidth)

void	SetUseEP (bool flag)

void	SetEPCutoff (float epcut)

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

virtual void	Report ()
	Function to report what's going on per merging. More...

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

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

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

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

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

virtual void	EventEnd ()

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

virtual void	IterationEnd ()

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

Protected Attributes
bool	_verbose

bool	_debug

bool	_use_EP

float	_ep_cut

size_t	_MinNumHits

float	_MinAngleDiff

float	_MaxOpeningAngle

float	_MinLength

float	_MinDensity

float	_MaxWidth

double	_wire_2_cm

double	_time_2_cm

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

Definition at line 24 of file CBAlgoTrackSeparate.h.

Constructor & Destructor Documentation

cmtool::CBAlgoTrackSeparate::CBAlgoTrackSeparate ( )

Default constructor.

Definition at line 8 of file CBAlgoTrackSeparate.cxx.

                                            : CBoolAlgoBase()
   //----------------------------------------
   {
 
     //this just sets default values
     SetVerbose(true);
     SetDebug(true);
 
     //1e9 is huge; everything will be merged
     SetMinNumHits(30);
     SetMinAngleDiff(15.); //in degrees
     SetMaxOpeningAngle(12.0); //in deg (parameter in rad!!)
     SetMinLength(10.);
     SetMinPolyHitDensity(10.);
     SetMaxWidth(10.);
 
 
     //NOTE! Using this flag means all of the other crap
     //(minNumHits, anglediff, openingangle, blah blah)
     //is totally irrelevant. if we stick with this flag as the algo,
     //we probably want to delete all of the old method
     SetUseEP(true);
     SetEPCutoff(0.99000);
     ::util::GeometryUtilities geou;
     _wire_2_cm = geou.WireToCm();
     _time_2_cm = geou.TimeToCm();
 
   }

virtual cmtool::CBAlgoTrackSeparate::~CBAlgoTrackSeparate ( )

inlinevirtual

Default destructor.

Definition at line 32 of file CBAlgoTrackSeparate.h.

32 {};

Member Function Documentation

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

virtual

Core function: given the ClusterParamsAlg input, return whether a cluster should be merged or not.

Reimplemented from cmtool::CBoolAlgoBase.

Definition at line 38 of file CBAlgoTrackSeparate.cxx.

   {
     //if you are using EP method for this algo:
     if(_use_EP){
       if(cluster1.GetParams().eigenvalue_principal > _ep_cut &&
          cluster2.GetParams().eigenvalue_principal > _ep_cut)
         return true;
     }
     //if you are using the original method for this algo:
     else{
 
       //two clusters are considered un-mergable if:
       //1) both have more than _MinNumHits
       //2) opening angle for both < _MAxOpeningAngle
       //3) diff. in direction of both < _MinAngleDiff
 
       size_t N_Hits1 = cluster1.GetHitVector().size();
       size_t N_Hits2 = cluster2.GetHitVector().size();
       auto start_point1 = cluster1.GetParams().start_point;
       auto start_point2 = cluster2.GetParams().start_point;
       double angle_2d1 = cluster1.GetParams().angle_2d;
       double angle_2d2 = cluster2.GetParams().angle_2d;
       double opening_angle1 = cluster1.GetParams().opening_angle;
       double opening_angle2 = cluster2.GetParams().opening_angle;
       Polygon2D PolyObject1 = cluster1.GetParams().PolyObject;
       Polygon2D PolyObject2 = cluster2.GetParams().PolyObject;
       double length1 = cluster1.GetParams().length;
       double length2 = cluster2.GetParams().length;
       double width1 = cluster1.GetParams().width;
       double width2 = cluster2.GetParams().width;
 
       //first filter out low hits clusters
       if ( (N_Hits1 > _MinNumHits) and
            (N_Hits2 > _MinNumHits) ) {
         if (_debug) {
           std::cout << "Cluster1 Num Hits: " << N_Hits1 << std::endl;
           std::cout << "\t Start: (" << start_point1.w << " " << start_point1.t << " )" << std::endl;
           std::cout << "\t Opening ANgle " << opening_angle1*(360/(2*3.14)) << std::endl;
           std::cout << "\t Angle2D: " << angle_2d1 << std::endl;
           std::cout << "\t Length: " << length1 << std::endl;
           std::cout << "\t Width: " << width1 << std::endl;
           std::cout << "Cluster2 Num Hits: " << N_Hits2 << std::endl;
           std::cout << "\t Start: (" << start_point2.w  << " " << start_point2.t << " )" << std::endl;
           std::cout << "\t Opening ANgle " << opening_angle2*(360/(2*3.14)) << std::endl;
           std::cout << "\t Angle2D: " << angle_2d2 << std::endl;
           std::cout << "\t Length: " << length2 << std::endl;
           std::cout << "\t Width: " << width2 << std::endl;
         }
         if ( (N_Hits1 > _MinNumHits) and
              (N_Hits2 > _MinNumHits) and
              ( abs(angle_2d1 - angle_2d2) > _MinAngleDiff ) and
              //( PolyObject1.Area()/N_Hits1 > _MinDensity ) and
              //( PolyObject2.Area()/N_Hits2 > _MinDensity ) and
              (opening_angle1 < _MaxOpeningAngle/(360/(2*3.14))) and
              (opening_angle2 < _MaxOpeningAngle/(360/(2*3.14))) and
              (width1 < _MaxWidth) and
              (width2 < _MaxWidth) and
              (length1 > _MinLength) and
              (length2 > _MinLength) ){
           if (_verbose) { std::cout << "*****************************************Separate with TrackSeparate!" << std::endl; }
           return true;
         }
       }
     }
 
     return false;
 
   }