#include <CBAlgoStartTrack.h>

Inheritance diagram for cmtool::CBAlgoStartTrack:

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

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

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

bool	IsStartTrack (const ::cluster::ClusterParamsAlg &cluster)

bool	IsOverlappingBlob (const ::cluster::ClusterParamsAlg &cluster)

void	SetMinWidth (double value)

void	SetMinOpeningAngle (double value)

void	SetMinEP (double value)

void	SetMinHits (size_t value)

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	EventBegin (const std::vector< cluster::ClusterParamsAlg > &)

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

Protected Attributes
size_t	_min_hits

double	_min_width

double	_min_opening_angle

double	_min_EP

bool	_debug

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 CBAlgoStartTrack.h.

Constructor & Destructor Documentation

cmtool::CBAlgoStartTrack::CBAlgoStartTrack ( )

Default constructor.

Definition at line 6 of file CBAlgoStartTrack.cxx.

                                      : CBoolAlgoBase()
   //-------------------------------------------------------
   {
 
     SetMinHits(10);
     SetMinWidth(1.00000);
     SetMinOpeningAngle(0.15000);
     SetMinEP(0.99000);
     SetDebug(false);
 
   }

virtual cmtool::CBAlgoStartTrack::~CBAlgoStartTrack ( )

inlinevirtual

Default destructor.

Definition at line 36 of file CBAlgoStartTrack.h.

36 {};

Member Function Documentation

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

virtual

Optional function: called at the beginning of 1st iteration. This is called per event. 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 CBAlgoStartTrack.cxx.

   {
 
     if(cluster1.GetHitVector().size() < _min_hits ||
        cluster2.GetHitVector().size() < _min_hits )
       return false;
 
     bool FirstTrackSecondBlob = false;
     bool SecondTrackFirstBlob = false;
 
     if(IsStartTrack(cluster1) && IsOverlappingBlob(cluster2))
       FirstTrackSecondBlob = true;
     if(IsStartTrack(cluster2) && IsOverlappingBlob(cluster1))
       SecondTrackFirstBlob = true;
 
     if(_debug){
       if(FirstTrackSecondBlob){
         std::cout<<"*********************************************"<<std::endl;
         std::cout<<"First track second blob! First has params:"<<std::endl;
         std::cout<<"Angle 2D = "<<cluster1.GetParams().angle_2d<<std::endl;
         std::cout<<"Opening Angle = "<<cluster1.GetParams().opening_angle<<std::endl;
         std::cout<<"Sum charge = "<<cluster1.GetParams().sum_charge<<std::endl;
         std::cout<<"Length = "<<cluster1.GetParams().length<<std::endl;
         std::cout<<"Width = "<<cluster1.GetParams().width<<std::endl;
         std::cout<<"N_Hits = "<<cluster1.GetHitVector().size()<<std::endl;
         std::cout<<"eigenvalue_principal = "<<cluster1.GetParams().eigenvalue_principal<<std::endl;
         std::cout<<"Second has params:"<<std::endl;
         std::cout<<"Angle 2D = "<<cluster2.GetParams().angle_2d<<std::endl;
         std::cout<<"Opening Angle = "<<cluster2.GetParams().opening_angle<<std::endl;
         std::cout<<"Sum charge = "<<cluster2.GetParams().sum_charge<<std::endl;
         std::cout<<"Length = "<<cluster2.GetParams().length<<std::endl;
         std::cout<<"Width = "<<cluster2.GetParams().width<<std::endl;
         std::cout<<"N_Hits = "<<cluster2.GetHitVector().size()<<std::endl;
         std::cout<<"eigenvalue_principal = "<<cluster2.GetParams().eigenvalue_principal<<std::endl;
 
       }
 
 
     }
 
     //if cluster1 looks like a track and cluster2 looks like a blob
     if (FirstTrackSecondBlob)
       {
         std::pair<float,float> trackEndPoint =
           std::pair<float,float>( cluster1.GetParams().end_point.w, cluster1.GetParams().end_point.t );
 
         //if the tracklike cluster's endpoint is inside polygon of blob
         if(cluster2.GetParams().PolyObject.PointInside(trackEndPoint))
           return true;
         else
           return false;
       }
 
     //if cluster2 looks like a track and cluster1 looks like a blob
     else if (SecondTrackFirstBlob)
       {
         std::pair<float,float> trackEndPoint =
           std::pair<float,float>( cluster2.GetParams().end_point.w, cluster2.GetParams().end_point.t );
 
         //if the tracklike cluster's endpoint is inside polygon of blob
         if(cluster1.GetParams().PolyObject.PointInside(trackEndPoint))
           return true;
         else
           return false;
       }
 
     else
       return false;
 
     //the reason for not just using (one track one blob) and
     //polyoverlap for this is that
     //this was merging BNB events where the track is clearly in a different
     //direction than the blob
 
     //this will fail when the tracklike cluster is reconstructed backwards
     //(looks like it's a track going into a blob, rather than exiting)
 
 
   }

bool cmtool::CBAlgoStartTrack::IsOverlappingBlob ( const ::cluster::ClusterParamsAlg & cluster )

Definition at line 153 of file CBAlgoStartTrack.cxx.

   {
     if(cluster.GetParams().width > _min_width &&
        cluster.GetParams().opening_angle > _min_opening_angle &&
        cluster.GetParams().eigenvalue_principal < _min_EP)
       return true;
     else
       return false;
   }

bool cmtool::CBAlgoStartTrack::IsStartTrack ( const ::cluster::ClusterParamsAlg & cluster )

Definition at line 144 of file CBAlgoStartTrack.cxx.

   {
     if(cluster.GetParams().eigenvalue_principal > _min_EP)
       return true;
     else
       return false;
 
   }

void cmtool::CBAlgoStartTrack::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 137 of file CBAlgoStartTrack.cxx.