#include "CBAlgoStartInCone.h"

Inheritance diagram for cmtool::CBAlgoStartInCone:

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

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

virtual bool	Bool (const ::cluster::ClusterParamsAlg &cluster1, const ::cluster::ClusterParamsAlg &cluster2)
	Merging Algorithm is Here. More...

void	reconfigure ()
	Method to re-configure the instance. More...

void	SetMinHits (size_t n)
	Set Minimum number of hits for cone-cluster. More...

void	SetMinLen (double l)
	Set Minimum number of hits for cone-cluster. More...

void	SetVerbose (bool verbosity)
	Set Verbosity of messages. More...

void	SetDebug (bool debug)
	Set Debug for messages. More...

void	SetAngleCompat (double deg)
	Set Angle Compatibility betweeen the clusters. More...

void	SetLengthReach (double frac)
	Set Length Reach: How for out the cone extends as percent of cluster length. More...

virtual void	Reset ()
	Function to reset the algorithm instance called within CMergeManager/CMatchManager's Reset() ... maybe implemented via child class. 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 ()

virtual void	Report ()

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

Protected Attributes
double	_wire_2_cm

double	_time_2_cm

size_t	_NhitsMin
	Conversion factors ogtten from GeometryUtilities. More...

double	_lenMin
	Larger cluster which determines cone must have this many hits. More...

bool	_verbose
	Larger cluster which determines cone must be at least this long. More...

bool	_debug

double	_angleCompat

double	_lengthReach
	Two clusters must have direction within this value of each other. More...

TFile *	_fout
	TFile pointer to an output file. More...

Detailed Description

If start point of Cluster B in Cone region of Cluster A then merge

Definition at line 28 of file CBAlgoStartInCone.h.

Constructor & Destructor Documentation

cmtool::CBAlgoStartInCone::CBAlgoStartInCone ( )

Default constructor.

Definition at line 8 of file CBAlgoStartInCone.cxx.

References _time_2_cm, _wire_2_cm, reconfigure(), SetAngleCompat(), SetDebug(), SetLengthReach(), SetMinHits(), SetMinLen(), util::GeometryUtilities::TimeToCm(), and util::GeometryUtilities::WireToCm().

   {
     // Nothing to be done in the base class
     this->reconfigure();
 
     ::util::GeometryUtilities geou;
     _wire_2_cm = geou.WireToCm();
     _time_2_cm = geou.TimeToCm();
 
     SetMinHits(40);
     SetMinLen(10);
     SetDebug(false);
     SetAngleCompat(90.);
     SetLengthReach(1.0);
 
   }

virtual cmtool::CBAlgoStartInCone::~CBAlgoStartInCone ( )

inlinevirtual

Default destructor.

Definition at line 36 of file CBAlgoStartInCone.h.

References Bool(), and reconfigure().

36 {};

Member Function Documentation

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

virtual

Merging Algorithm is Here.

Reimplemented from cmtool::CBoolAlgoBase.

Definition at line 33 of file CBAlgoStartInCone.cxx.

References _angleCompat, _debug, _lengthReach, _lenMin, _NhitsMin, and _verbose.

Referenced by ~CBAlgoStartInCone().

   {
 
     //apply filter on hit number and length immediately
     if ( ( (cluster1.GetParams().length < _lenMin) and (cluster1.GetHitVector().size() < _NhitsMin) )
          and ( (cluster2.GetParams().length < _lenMin) and (cluster2.GetHitVector().size() < _NhitsMin) ) )
       return false;
 
     //A cluster has an associated cone defined as the cone originating
     //at the start point of that cluster opening out in the direction
     //of the end point with an angle equal to the opening angle
     //If cluster A (B) has a start point within the volume defined
     //by the cone of cluster B (A), and if cluster B (A) is good enough
     // (enough hits, charge, length...) then the two are merged.
 
     double w_start1 = cluster1.GetParams().start_point.w;// * _wire_2_cm;
     double t_start1 = cluster1.GetParams().start_point.t;// * _time_2_cm;
     double angle1   = cluster1.GetParams().angle_2d;
     double opening1 = cluster1.GetParams().opening_angle * (180./3.14);
     double len1     = cluster1.GetParams().length;
 
     double w_start2 = cluster2.GetParams().start_point.w;// * _wire_2_cm;
     double t_start2 = cluster2.GetParams().start_point.t;// * _time_2_cm;
     double angle2   = cluster2.GetParams().angle_2d;
     double opening2 = cluster2.GetParams().opening_angle * (180./3.14);
     double len2     = cluster2.GetParams().length;
 
     //check for angle compatibility: the direction of the two clusters must be within
     //this range of each other
     if ( abs(angle2-angle1) > _angleCompat ){
       if (_debug) { std::cout << "Directions too different....do not merge" << std::endl; }
       return false;
     }
     
     if (_debug){
       std::cout << "Cluster 1:" << std::endl;
       std::cout << "\tStart: ( " << w_start1 << ", " << t_start1 << " )" << std::endl;
       std::cout << "\tAngle: " << angle1 << std::endl;
       std::cout << "\tOpening Angle: " << opening1 << std::endl;
       std::cout << "\tLength: " << len1 << std::endl;
       std::cout << "Cluster 2:" << std::endl;
       std::cout << "\tStart: ( " << w_start2 << ", " << t_start2 << " )" << std::endl;
       std::cout << "\tAngle: " << angle2 << std::endl;
       std::cout << "\tOpening Angle: " << opening2 << std::endl;
       std::cout << "\tLength: " << len2 << std::endl;
       std::cout << std::endl;
     }
     
     //find if start point of A (B) in cone of B (A)
     //do this by translating point A (B) such that
     //start point of B (A) is at origin and +  axis
     //towards end point of B (A)
     double w_start1_transl = w_start1 - w_start2;
     double t_start1_transl = t_start1 - t_start2;
     double w_start2_transl = w_start2 - w_start1;
     double t_start2_transl = t_start2 - t_start1;
 
     double w_start1_rot =   w_start1_transl*cos(angle2*3.14/180.) + t_start1_transl*sin(angle2*3.14/180.);
     double t_start1_rot = - w_start1_transl*sin(angle2*3.14/180.) + t_start1_transl*cos(angle2*3.14/180.);
     double w_start2_rot =   w_start2_transl*cos(angle1*3.14/180.) + t_start2_transl*sin(angle1*3.14/180.);
     double t_start2_rot = - w_start2_transl*sin(angle1*3.14/180.) + t_start2_transl*cos(angle1*3.14/180.);
 
     if (_debug){
       std::cout << "\tStart Rot 1: ( " << w_start1_rot << ", " << t_start1_rot << " )" << std::endl;
       std::cout << "\tStart Rot 2: ( " << w_start2_rot << ", " << t_start2_rot << " )" << std::endl;
     }
 
     if ( (w_start1_rot < len2*_lengthReach ) and (w_start1_rot > 0) and
          ( abs(t_start1_rot) < (abs(w_start1_rot*sin(opening2*3.14/180.))) ) and
          (len2 > _lenMin) and 
          (cluster2.GetHitVector().size() > _NhitsMin) ){
       if (_verbose) { std::cout << "Start point of Cluster 1 in cone of Cluster 2!" << std::endl; }
       return true;
     }
     if ( (w_start2_rot < len1*_lengthReach ) and (w_start2_rot > 0) and
          ( abs(t_start2_rot) < (abs(w_start2_rot*sin(opening1*3.14/180.))) ) and
          (len1 > _lenMin) and 
          (cluster1.GetHitVector().size() > _NhitsMin) ){
       if (_verbose) { std::cout << "Start point of Cluster 2 in cone of Cluster 1!" << std::endl; }
       return true;
     }
 
     return false;
   }

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

inlinevirtualinherited

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

Reimplemented in cmtool::CFAlgoArray, cmtool::CPAlgoArray, cmtool::CBAlgoArray, and cmtool::CBAlgoPolyShortestDist.

Definition at line 45 of file CMAlgoBase.h.

Referenced by cmtool::CMergeManager::EventBegin().

46 { if(clusters.size()) return; }

virtual void cmtool::CMAlgoBase::EventEnd ( )

inlinevirtualinherited

Optional function: called at the end of event ... after the last merging iteration is over.

Reimplemented in cmtool::CFAlgoArray, cmtool::CPAlgoArray, and cmtool::CBAlgoArray.

Definition at line 51 of file CMAlgoBase.h.

Referenced by cmtool::CMatchManager::EventEnd(), and cmtool::CMergeManager::EventEnd().

52 {return;}

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

inlinevirtualinherited

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.

Reimplemented in cmtool::CFAlgoArray, cmtool::CPAlgoArray, and cmtool::CBAlgoArray.

Definition at line 59 of file CMAlgoBase.h.

Referenced by cmtool::CMatchManager::EventBegin(), cmtool::CMatchManager::IterationBegin(), and cmtool::CMergeManager::IterationBegin().

60 { if(clusters.size()) return;}

virtual void cmtool::CMAlgoBase::IterationEnd ( )

inlinevirtualinherited

Optional function: called at the end of each iteration over all pairs of clusters.

Reimplemented in cmtool::CFAlgoArray, cmtool::CPAlgoArray, and cmtool::CBAlgoArray.

Definition at line 65 of file CMAlgoBase.h.

Referenced by cmtool::CMatchManager::IterationEnd(), and cmtool::CMergeManager::IterationEnd().

66 {return; }

void cmtool::CBAlgoStartInCone::reconfigure ( )

Method to re-configure the instance.

Definition at line 26 of file CBAlgoStartInCone.cxx.

Referenced by CBAlgoStartInCone(), and ~CBAlgoStartInCone().

                                      {
 
     //not sure what needs to be reset/reconfigured for this algo
     
   }//end reconfigure function

virtual void cmtool::CMAlgoBase::Report ( )

inlinevirtualinherited

Optional function: called after Bool() function is called for all possible cluster pairs by CMergeManager/CMatchManager IFF run with verbosity level kPerIteration. Maybe useful for debugging.

Reimplemented in cmtool::CFAlgoArray, cmtool::CPAlgoArray, cmtool::CBAlgoArray, cmtool::CBAlgoStartTrack, cmtool::CBAlgoPolyShortestDist, cmtool::CBAlgoProhibitAllTracks, cmtool::CBAlgoOutOfConeSeparate, cmtool::CBAlgoTrackSeparate, cmtool::CBAlgoAngleSeparate, cmtool::CFAlgo3DAngle, cmtool::CBAlgoMergeTinyWithBig, cmtool::CFAlgoShowerCompat, cmtool::CFAlgoTimeOverlap, cmtool::CBAlgoStartNearEnd, cmtool::CFAlgoVolumeOverlap, cmtool::CFAlgoWireOverlap, cmtool::CBAlgoProhibitBigClusters, cmtool::CFAlgoStartPointCompat, cmtool::CFAlgoStartPointMatch, cmtool::CFAlgoStartTimeCompat, cmtool::CFAlgoChargeDistrib, cmtool::CFAlgoQRatio, cmtool::CFAlgoZOverlap, cmtool::CPAlgoIgnoreTracks, cmtool::CBAlgoCenterOfMass, cmtool::CBAlgoCenterOfMassSmall, and cmtool::CFAlgoTimeProf.

Definition at line 73 of file CMAlgoBase.h.

Referenced by cmtool::CMergeManager::IterationEnd(), cmtool::CMatchManager::IterationProcess(), and cmtool::CMergeManager::IterationProcess().

74 {return;}

virtual void cmtool::CMAlgoBase::Reset ( )

inlinevirtualinherited

Function to reset the algorithm instance called within CMergeManager/CMatchManager's Reset() ... maybe implemented via child class.

Reimplemented in cmtool::CBAlgoArray, cmtool::CFAlgoArray, cmtool::CPAlgoArray, cmtool::CBAlgoStartTrack, cmtool::CBAlgoPolyShortestDist, cmtool::CBAlgoProhibitAllTracks, cmtool::CBAlgoOutOfConeSeparate, cmtool::CBAlgoTrackSeparate, cmtool::CBAlgoMergeTinyWithBig, cmtool::CFAlgoShowerCompat, cmtool::CFAlgo3DAngle, cmtool::CFAlgoTimeOverlap, cmtool::CBAlgoAngleSeparate, cmtool::CFAlgoVolumeOverlap, cmtool::CFAlgoWireOverlap, cmtool::CBAlgoProhibitBigClusters, cmtool::CFAlgoStartPointCompat, cmtool::CFAlgoStartPointMatch, cmtool::CFAlgoStartTimeCompat, cmtool::CFAlgoChargeDistrib, cmtool::CFAlgoQRatio, cmtool::CFAlgoZOverlap, cmtool::CPAlgoIgnoreTracks, cmtool::CBAlgoStartNearEnd, cmtool::CBAlgoCenterOfMass, cmtool::CBAlgoCenterOfMassSmall, cmtool::CFAlgoTimeProf, cmtool::CBAlgoFake, and cmtool::CBAlgoMergeAll.

Definition at line 40 of file CMAlgoBase.h.

Referenced by cmtool::CMergeManager::Reset(), cmtool::CMatchManager::Reset(), and cmtool::CMManagerBase::Reset().

40 {}

void cmtool::CMAlgoBase::SetAnaFile ( TFile * fout )

inlineinherited

Setter function for an output plot TFile pointer.

Definition at line 77 of file CMAlgoBase.h.

References cmtool::CMAlgoBase::_fout.

Referenced by cmtool::CMergeManager::EventBegin().

77 { _fout = fout; }

cmtool::CMAlgoBase::_fout

TFile * _fout

TFile pointer to an output file.

Definition: CMAlgoBase.h:85

void cmtool::CBAlgoStartInCone::SetAngleCompat ( double deg )

inline

Set Angle Compatibility betweeen the clusters.

Definition at line 58 of file CBAlgoStartInCone.h.

References _angleCompat.

Referenced by CBAlgoStartInCone().

58 { _angleCompat = deg; }

cmtool::CBAlgoStartInCone::_angleCompat

double _angleCompat

Definition: CBAlgoStartInCone.h:71

void cmtool::CBAlgoStartInCone::SetDebug ( bool debug )

inline

Set Debug for messages.

Definition at line 55 of file CBAlgoStartInCone.h.

References _debug, and tca::debug.

Referenced by CBAlgoStartInCone().

55 { _debug = debug; }

tca::debug

DebugStuff debug

Definition: DebugStruct.cxx:4

cmtool::CBAlgoStartInCone::_debug

bool _debug

Definition: CBAlgoStartInCone.h:70

void cmtool::CBAlgoStartInCone::SetLengthReach ( double frac )

inline

Set Length Reach: How for out the cone extends as percent of cluster length.

Definition at line 61 of file CBAlgoStartInCone.h.

References _lengthReach.

Referenced by CBAlgoStartInCone().

61 { _lengthReach = frac; }

cmtool::CBAlgoStartInCone::_lengthReach

double _lengthReach

Two clusters must have direction within this value of each other.

Definition: CBAlgoStartInCone.h:72

void cmtool::CBAlgoStartInCone::SetMinHits ( size_t n )

inline

Set Minimum number of hits for cone-cluster.

Definition at line 46 of file CBAlgoStartInCone.h.

References _NhitsMin, and n.

Referenced by CBAlgoStartInCone().

46 { _NhitsMin = n; }

cmtool::CBAlgoStartInCone::_NhitsMin

size_t _NhitsMin

Conversion factors ogtten from GeometryUtilities.

Definition: CBAlgoStartInCone.h:67

n

Char_t n[5]

Definition: comparison_ascii.C:55

void cmtool::CBAlgoStartInCone::SetMinLen ( double l )

inline

Set Minimum number of hits for cone-cluster.

Definition at line 49 of file CBAlgoStartInCone.h.

References _lenMin.

Referenced by CBAlgoStartInCone().

49 { _lenMin = l; }

cmtool::CBAlgoStartInCone::_lenMin

double _lenMin

Larger cluster which determines cone must have this many hits.

Definition: CBAlgoStartInCone.h:68

void cmtool::CBAlgoStartInCone::SetVerbose ( bool verbosity )

inlinevirtual

Set Verbosity of messages.

Reimplemented from cmtool::CMAlgoBase.

Definition at line 52 of file CBAlgoStartInCone.h.

References _verbose.

52 { _verbose = verbosity; }

cmtool::CBAlgoStartInCone::_verbose

bool _verbose

Larger cluster which determines cone must be at least this long.

Definition: CBAlgoStartInCone.h:69

Member Data Documentation

double cmtool::CBAlgoStartInCone::_angleCompat

protected

Definition at line 71 of file CBAlgoStartInCone.h.

Referenced by Bool(), and SetAngleCompat().

bool cmtool::CBAlgoStartInCone::_debug

protected

Definition at line 70 of file CBAlgoStartInCone.h.

Referenced by Bool(), and SetDebug().

TFile* cmtool::CMAlgoBase::_fout

protectedinherited

TFile pointer to an output file.

Definition at line 85 of file CMAlgoBase.h.

Referenced by cmtool::CMAlgoBase::CMAlgoBase(), and cmtool::CMAlgoBase::SetAnaFile().

double cmtool::CBAlgoStartInCone::_lengthReach

protected

Two clusters must have direction within this value of each other.

Definition at line 72 of file CBAlgoStartInCone.h.

Referenced by Bool(), and SetLengthReach().

double cmtool::CBAlgoStartInCone::_lenMin

protected

Larger cluster which determines cone must have this many hits.

Definition at line 68 of file CBAlgoStartInCone.h.

Referenced by Bool(), and SetMinLen().

size_t cmtool::CBAlgoStartInCone::_NhitsMin

protected

Conversion factors ogtten from GeometryUtilities.

Definition at line 67 of file CBAlgoStartInCone.h.

Referenced by Bool(), and SetMinHits().

double cmtool::CBAlgoStartInCone::_time_2_cm

protected

Definition at line 65 of file CBAlgoStartInCone.h.

Referenced by CBAlgoStartInCone().

bool cmtool::CBAlgoStartInCone::_verbose

protected

Larger cluster which determines cone must be at least this long.

Definition at line 69 of file CBAlgoStartInCone.h.

Referenced by Bool(), and SetVerbose().

double cmtool::CBAlgoStartInCone::_wire_2_cm

protected

Definition at line 65 of file CBAlgoStartInCone.h.

Referenced by CBAlgoStartInCone().

The documentation for this class was generated from the following files:

larreco/v06_64_02/source/larreco/RecoAlg/CMTool/CMTAlgMerge/CBAlgoStartInCone.h
larreco/v06_64_02/source/larreco/RecoAlg/CMTool/CMTAlgMerge/CBAlgoStartInCone.cxx

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation