cmtool::CBAlgoTrackSeparate Class Reference

#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...
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
TFile *	_fout
	TFile pointer to an output file. More...

Detailed Description

Definition at line 25 of file CBAlgoTrackSeparate.h.

Constructor & Destructor Documentation

cmtool::CBAlgoTrackSeparate::CBAlgoTrackSeparate

(

)

Default constructor.

Definition at line 9 of file CBAlgoTrackSeparate.cxx.

References _time_2_cm, _wire_2_cm, SetDebug(), SetEPCutoff(), SetMaxOpeningAngle(), SetMaxWidth(), SetMinAngleDiff(), SetMinLength(), SetMinNumHits(), SetMinPolyHitDensity(), SetUseEP(), SetVerbose(), util::GeometryUtilities::TimeToCm(), and util::GeometryUtilities::WireToCm().

                                           : 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 33 of file CBAlgoTrackSeparate.h.

References Bool().

{};

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 39 of file CBAlgoTrackSeparate.cxx.

References _debug, _ep_cut, _MaxOpeningAngle, _MaxWidth, _MinAngleDiff, _MinLength, _MinNumHits, _use_EP, and _verbose.

Referenced by ~CBAlgoTrackSeparate().

  {
    //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;

  }

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

    { 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().

    {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().

    { 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().

    {return; }

void cmtool::CBAlgoTrackSeparate::Report ( )

virtual

Function to report what's going on per merging.

Reimplemented from cmtool::CMAlgoBase.

Definition at line 111 of file CBAlgoTrackSeparate.cxx.

Referenced by Reset().

  {
  }

virtual void cmtool::CBAlgoTrackSeparate::Reset ( )

inlinevirtual

Function to reset the algorithm instance ... maybe implemented via child class.

Reimplemented from cmtool::CMAlgoBase.

Definition at line 63 of file CBAlgoTrackSeparate.h.

References Report().

{}

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

{ _fout = fout; }

void cmtool::CBAlgoTrackSeparate::SetDebug ( bool  on )

inline

Definition at line 44 of file CBAlgoTrackSeparate.h.

References _debug.

Referenced by CBAlgoTrackSeparate(), and cluster::SimpleClusterMerger::SimpleClusterMerger().

{ _debug = on; }

void cmtool::CBAlgoTrackSeparate::SetEPCutoff ( float  epcut )

inline

Definition at line 60 of file CBAlgoTrackSeparate.h.

References _ep_cut.

Referenced by CBAlgoTrackSeparate().

{ _ep_cut = epcut; }

void cmtool::CBAlgoTrackSeparate::SetMaxOpeningAngle ( float  maxangle )

inline

Definition at line 50 of file CBAlgoTrackSeparate.h.

References _MaxOpeningAngle.

Referenced by CBAlgoTrackSeparate(), and cluster::SimpleClusterMerger::SimpleClusterMerger().

{ _MaxOpeningAngle = maxangle; }

void cmtool::CBAlgoTrackSeparate::SetMaxWidth ( float  maxwidth )

inline

Definition at line 56 of file CBAlgoTrackSeparate.h.

References _MaxWidth.

Referenced by CBAlgoTrackSeparate().

{ _MaxWidth = maxwidth; }

void cmtool::CBAlgoTrackSeparate::SetMinAngleDiff ( float  anglesep )

inline

Definition at line 48 of file CBAlgoTrackSeparate.h.

References _MinAngleDiff.

Referenced by CBAlgoTrackSeparate(), and cluster::SimpleClusterMerger::SimpleClusterMerger().

{ _MinAngleDiff = anglesep; }

void cmtool::CBAlgoTrackSeparate::SetMinLength ( float  minlength )

inline

Definition at line 52 of file CBAlgoTrackSeparate.h.

References _MinLength.

Referenced by CBAlgoTrackSeparate(), and cluster::SimpleClusterMerger::SimpleClusterMerger().

{ _MinLength = minlength; }

void cmtool::CBAlgoTrackSeparate::SetMinNumHits ( size_t  n )

inline

Definition at line 46 of file CBAlgoTrackSeparate.h.

References _MinNumHits, and n.

Referenced by CBAlgoTrackSeparate(), and cluster::SimpleClusterMerger::SimpleClusterMerger().

{ _MinNumHits = n; }

void cmtool::CBAlgoTrackSeparate::SetMinPolyHitDensity ( float  mindensity )

inline

Definition at line 54 of file CBAlgoTrackSeparate.h.

References _MinDensity.

Referenced by CBAlgoTrackSeparate().

{ _MinDensity = mindensity; }

void cmtool::CBAlgoTrackSeparate::SetUseEP ( bool  flag )

inline

Definition at line 58 of file CBAlgoTrackSeparate.h.

References _use_EP.

Referenced by CBAlgoTrackSeparate().

{ _use_EP = flag; }

void cmtool::CBAlgoTrackSeparate::SetVerbose ( bool  doit )

inlinevirtual

Setter function for verbosity.

Reimplemented from cmtool::CMAlgoBase.

Definition at line 42 of file CBAlgoTrackSeparate.h.

References _verbose.

Referenced by CBAlgoTrackSeparate(), and cluster::SimpleClusterMerger::SimpleClusterMerger().

{ _verbose = on; }

Member Data Documentation

bool cmtool::CBAlgoTrackSeparate::_debug

protected

Definition at line 71 of file CBAlgoTrackSeparate.h.

Referenced by Bool(), and SetDebug().

float cmtool::CBAlgoTrackSeparate::_ep_cut

protected

Definition at line 73 of file CBAlgoTrackSeparate.h.

Referenced by Bool(), and SetEPCutoff().

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

float cmtool::CBAlgoTrackSeparate::_MaxOpeningAngle

protected

Definition at line 76 of file CBAlgoTrackSeparate.h.

Referenced by Bool(), and SetMaxOpeningAngle().

float cmtool::CBAlgoTrackSeparate::_MaxWidth

protected

Definition at line 79 of file CBAlgoTrackSeparate.h.

Referenced by Bool(), and SetMaxWidth().

float cmtool::CBAlgoTrackSeparate::_MinAngleDiff

protected

Definition at line 75 of file CBAlgoTrackSeparate.h.

Referenced by Bool(), and SetMinAngleDiff().

float cmtool::CBAlgoTrackSeparate::_MinDensity

protected

Definition at line 78 of file CBAlgoTrackSeparate.h.

Referenced by SetMinPolyHitDensity().

float cmtool::CBAlgoTrackSeparate::_MinLength

protected

Definition at line 77 of file CBAlgoTrackSeparate.h.

Referenced by Bool(), and SetMinLength().

size_t cmtool::CBAlgoTrackSeparate::_MinNumHits

protected

Definition at line 74 of file CBAlgoTrackSeparate.h.

Referenced by Bool(), and SetMinNumHits().

double cmtool::CBAlgoTrackSeparate::_time_2_cm

protected

Definition at line 81 of file CBAlgoTrackSeparate.h.

Referenced by CBAlgoTrackSeparate().

bool cmtool::CBAlgoTrackSeparate::_use_EP

protected

Definition at line 72 of file CBAlgoTrackSeparate.h.

Referenced by Bool(), and SetUseEP().

bool cmtool::CBAlgoTrackSeparate::_verbose

protected

Definition at line 70 of file CBAlgoTrackSeparate.h.

Referenced by Bool(), and SetVerbose().

double cmtool::CBAlgoTrackSeparate::_wire_2_cm

protected

Definition at line 80 of file CBAlgoTrackSeparate.h.

Referenced by CBAlgoTrackSeparate().

---

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

• larreco/v06_64_02/source/larreco/RecoAlg/CMTool/CMTAlgMerge/CBAlgoTrackSeparate.h
• larreco/v06_64_02/source/larreco/RecoAlg/CMTool/CMTAlgMerge/CBAlgoTrackSeparate.cxx