cmtool::CBAlgoStartInCone Class Reference

#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...
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	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...
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

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

Definition at line 26 of file CBAlgoStartInCone.h.

Constructor & Destructor Documentation

cmtool::CBAlgoStartInCone::CBAlgoStartInCone

(

)

Default constructor.

Definition at line 7 of file CBAlgoStartInCone.cxx.

  {
    // 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 34 of file CBAlgoStartInCone.h.

{};

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 32 of file CBAlgoStartInCone.cxx.

  {

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

void cmtool::CBAlgoStartInCone::reconfigure

(

)

Method to re-configure the instance.

Definition at line 25 of file CBAlgoStartInCone.cxx.

                                     {

    //not sure what needs to be reset/reconfigured for this algo

  }//end reconfigure function

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

inline

Set Angle Compatibility betweeen the clusters.

Definition at line 56 of file CBAlgoStartInCone.h.

{ _angleCompat = deg; }

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

inline

Set Debug for messages.

Definition at line 53 of file CBAlgoStartInCone.h.

{ _debug = debug; }

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

inline

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

Definition at line 59 of file CBAlgoStartInCone.h.

{ _lengthReach = frac; }

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

inline

Set Minimum number of hits for cone-cluster.

Definition at line 44 of file CBAlgoStartInCone.h.

{ _NhitsMin = n; }

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

inline

Set Minimum number of hits for cone-cluster.

Definition at line 47 of file CBAlgoStartInCone.h.

{ _lenMin = l; }

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

inlinevirtual

Set Verbosity of messages.

Reimplemented from cmtool::CMAlgoBase.

Definition at line 50 of file CBAlgoStartInCone.h.

{ _verbose = verbosity; }

Member Data Documentation

double cmtool::CBAlgoStartInCone::_angleCompat

protected

Definition at line 69 of file CBAlgoStartInCone.h.

bool cmtool::CBAlgoStartInCone::_debug

protected

Definition at line 68 of file CBAlgoStartInCone.h.

double cmtool::CBAlgoStartInCone::_lengthReach

protected

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

Definition at line 70 of file CBAlgoStartInCone.h.

double cmtool::CBAlgoStartInCone::_lenMin

protected

Larger cluster which determines cone must have this many hits.

Definition at line 66 of file CBAlgoStartInCone.h.

size_t cmtool::CBAlgoStartInCone::_NhitsMin

protected

Conversion factors ogtten from GeometryUtilities.

Definition at line 65 of file CBAlgoStartInCone.h.

double cmtool::CBAlgoStartInCone::_time_2_cm

protected

Definition at line 63 of file CBAlgoStartInCone.h.

bool cmtool::CBAlgoStartInCone::_verbose

protected

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

Definition at line 67 of file CBAlgoStartInCone.h.

double cmtool::CBAlgoStartInCone::_wire_2_cm

protected

Definition at line 63 of file CBAlgoStartInCone.h.

---

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

• larreco/larreco/RecoAlg/CMTool/CMTAlgMerge/CBAlgoStartInCone.h
• larreco/larreco/RecoAlg/CMTool/CMTAlgMerge/CBAlgoStartInCone.cxx