#include <TwoViewLongTracksTool.h>

Inheritance diagram for lar_content::TwoViewLongTracksTool:

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

bool	Run (TwoViewTransverseTracksAlgorithm *const pAlgorithm, MatrixType &overlapMatrix)
	Run the algorithm tool. More...

Static Public Member Functions
static bool	HasLongDirectConnections (IteratorList::const_iterator iIter, const IteratorList &iteratorList)
	Whether a long element shares clusters with any other long elements. More...

static bool	IsLongerThanDirectConnections (IteratorList::const_iterator iIter, const MatrixType::ElementList &elementList, const unsigned int minMatchedSamplingPointRatio, const pandora::ClusterSet &usedClusters)
	Whether a long element is significantly longer that other elements with which it shares a cluster. More...

Private Member Functions
pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)

void	FindLongTracks (const MatrixType &overlapMatrix, ProtoParticleVector &protoParticleVector) const
	Find long tracks, hidden by simple ambiguities in the matrix. More...

void	SelectLongElements (const MatrixType::ElementList &elementList, const pandora::ClusterSet &usedClusters, IteratorList &iteratorList) const
	Select a list of long track-like elements from a set of connected matrix elements. More...

Private Attributes
float	m_minMatchedFraction
	The min matched sampling point fraction for particle creation. More...

float	m_minMatchingScore
	The min global matching score for particle creation. More...

unsigned int	m_minMatchedSamplingPoints
	The min number of matched sampling points for particle creation. More...

float	m_minXOverlapFraction
	The min x overlap fraction (in each view) for particle creation. More...

unsigned int	m_minMatchedSamplingPointRatio
	The min ratio between 1st and 2nd highest msps for simple ambiguity resolution. More...

Additional Inherited Members
Public Types inherited from lar_content::TransverseMatrixTool
typedef TwoViewTransverseTracksAlgorithm::MatchingType::MatrixType	MatrixType

typedef std::vector < MatrixType::ElementList::const_iterator >	IteratorList

Detailed Description

TwoViewLongTracksTool class.

Definition at line 19 of file TwoViewLongTracksTool.h.

Constructor & Destructor Documentation

lar_content::TwoViewLongTracksTool::TwoViewLongTracksTool ( )

Default constructor.

Definition at line 18 of file TwoViewLongTracksTool.cc.

                                              :
     m_minMatchedFraction(0.3f),
     m_minMatchingScore(0.95f),
     m_minMatchedSamplingPoints(15),
     m_minXOverlapFraction(0.9f),
     m_minMatchedSamplingPointRatio(2)
 {
 }

Member Function Documentation

void lar_content::TwoViewLongTracksTool::FindLongTracks	(	const MatrixType &	overlapMatrix,
		ProtoParticleVector &	protoParticleVector
	)		const

private

Find long tracks, hidden by simple ambiguities in the matrix.

Parameters

overlapMatrix	the overlap matrix
protoParticleVector	to receive the list of proto particles

Definition at line 84 of file TwoViewLongTracksTool.cc.

 {
     ClusterSet usedClusters;
     ClusterVector sortedKeyClusters;
     overlapMatrix.GetSortedKeyClusters(sortedKeyClusters);
 
     for (const Cluster *const pKeyCluster : sortedKeyClusters)
     {
         if (!pKeyCluster->IsAvailable())
             continue;
 
         unsigned int n0(0), n1(0);
         MatrixType::ElementList elementList;
         overlapMatrix.GetConnectedElements(pKeyCluster, true, elementList, n0, n1);
 
         IteratorList iteratorList;
         this->SelectLongElements(elementList, usedClusters, iteratorList);
 
         // Check that elements are not directly connected and are significantly longer than any other directly connected elements
         for (IteratorList::const_iterator iIter = iteratorList.begin(), iIterEnd = iteratorList.end(); iIter != iIterEnd; ++iIter)
         {
             if (TwoViewLongTracksTool::HasLongDirectConnections(iIter, iteratorList))
                 continue;
 
             if (!TwoViewLongTracksTool::IsLongerThanDirectConnections(iIter, elementList, m_minMatchedSamplingPointRatio, usedClusters))
                 continue;
 
             ProtoParticle protoParticle;
             protoParticle.m_clusterList.push_back((*iIter)->GetCluster1());
             protoParticle.m_clusterList.push_back((*iIter)->GetCluster2());
             protoParticleVector.push_back(protoParticle);
 
             usedClusters.insert((*iIter)->GetCluster1());
             usedClusters.insert((*iIter)->GetCluster2());
         }
     }
 }

bool lar_content::TwoViewLongTracksTool::HasLongDirectConnections	(	IteratorList::const_iterator	iIter,
		const IteratorList &	iteratorList
	)

static

Whether a long element shares clusters with any other long elements.

Parameters

iIter	specifies the long element under consideration
iteratorList	list of iterators to other long elements

Returns: boolean

Definition at line 29 of file TwoViewLongTracksTool.cc.

 {
     for (IteratorList::const_iterator iIter2 = iteratorList.begin(), iIter2End = iteratorList.end(); iIter2 != iIter2End; ++iIter2)
     {
         if (iIter == iIter2)
             continue;
 
         if (((*iIter)->GetCluster1() == (*iIter2)->GetCluster1()) || ((*iIter)->GetCluster2() == (*iIter2)->GetCluster2()))
             return true;
     }
 
     return false;
 }

bool lar_content::TwoViewLongTracksTool::IsLongerThanDirectConnections	(	IteratorList::const_iterator	iIter,
		const MatrixType::ElementList &	elementList,
		const unsigned int	minMatchedSamplingPointRatio,
		const pandora::ClusterSet &	usedClusters
	)

static

Whether a long element is significantly longer that other elements with which it shares a cluster.

Parameters

iIter	specifies the long element under consideration
elementList	the full list of connected matrix elements
minMatchedSamplingPointRatio	the min ratio between 1st and 2nd highest msps for simple ambiguity resolution
usedClusters	the list of clusters already marked as to be added to a pfo

Definition at line 45 of file TwoViewLongTracksTool.cc.

 {
     const unsigned int nMatchedReUpsampledSamplingPoints((*iIter)->GetOverlapResult().GetNMatchedReUpsampledSamplingPoints());
 
     for (MatrixType::ElementList::const_iterator eIter = elementList.begin(); eIter != elementList.end(); ++eIter)
     {
         if ((*iIter) == eIter)
             continue;
 
         if (usedClusters.count(eIter->GetCluster1()) || usedClusters.count(eIter->GetCluster2()))
             continue;
 
         if (((*iIter)->GetCluster1() != eIter->GetCluster1()) && ((*iIter)->GetCluster2() != eIter->GetCluster2()))
             continue;
 
         if (nMatchedReUpsampledSamplingPoints < minMatchedSamplingPointRatio * eIter->GetOverlapResult().GetNMatchedReUpsampledSamplingPoints())
             return false;
     }
 
     return true;
 }

StatusCode lar_content::TwoViewLongTracksTool::ReadSettings ( const pandora::TiXmlHandle xmlHandle )

private

Definition at line 152 of file TwoViewLongTracksTool.cc.

 {
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MinMatchedFraction", m_minMatchedFraction));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MinMatchingScore", m_minMatchingScore));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
         XmlHelper::ReadValue(xmlHandle, "MinMatchedSamplingPoints", m_minMatchedSamplingPoints));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MinXOverlapFraction", m_minXOverlapFraction));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
         XmlHelper::ReadValue(xmlHandle, "MinMatchedSamplingPointRatio", m_minMatchedSamplingPointRatio));
 
     return STATUS_CODE_SUCCESS;
 }

bool lar_content::TwoViewLongTracksTool::Run	(	TwoViewTransverseTracksAlgorithm *const	pAlgorithm,
		MatrixType &	overlapMatrix
	)

virtual

Run the algorithm tool.

Parameters

pAlgorithm	address of the calling algorithm
overlapMatrix	the overlap matrix

Returns: whether changes have been made by the tool

Implements lar_content::TransverseMatrixTool.

Definition at line 70 of file TwoViewLongTracksTool.cc.

 {
     if (PandoraContentApi::GetSettings(*pAlgorithm)->ShouldDisplayAlgorithmInfo())
         std::cout << "----> Running Algorithm Tool: " << this->GetInstanceName() << ", " << this->GetType() << std::endl;
 
     ProtoParticleVector protoParticleVector;
     this->FindLongTracks(overlapMatrix, protoParticleVector);
 
     const bool particlesMade(pAlgorithm->CreateThreeDParticles(protoParticleVector));
     return particlesMade;
 }

void lar_content::TwoViewLongTracksTool::SelectLongElements	(	const MatrixType::ElementList &	elementList,
		const pandora::ClusterSet &	usedClusters,
		IteratorList &	iteratorList
	)		const

private

Select a list of long track-like elements from a set of connected matrix elements.

Parameters

elementList	the full list of connected matrix elements
usedClusters	the list of clusters already marked as to be added to a pfo
iteratorList	to receive a list of iterators to long track-like elements

Definition at line 124 of file TwoViewLongTracksTool.cc.

 {
     for (MatrixType::ElementList::const_iterator eIter = elementList.begin(); eIter != elementList.end(); ++eIter)
     {
         if (usedClusters.count(eIter->GetCluster1()) || usedClusters.count(eIter->GetCluster2()))
             continue;
 
         if (eIter->GetOverlapResult().GetLocallyMatchedFraction() < m_minMatchedFraction)
             continue;
 
         if (eIter->GetOverlapResult().GetMatchingScore() < m_minMatchingScore)
             continue;
 
         if (eIter->GetOverlapResult().GetNMatchedReUpsampledSamplingPoints() < m_minMatchedSamplingPoints)
             continue;
 
         const TwoViewXOverlap &xOverlap(eIter->GetOverlapResult().GetTwoViewXOverlap());
 
         if ((xOverlap.GetXOverlapFraction0() > m_minXOverlapFraction) && (xOverlap.GetXOverlapFraction1() > m_minXOverlapFraction))
         {
             iteratorList.push_back(eIter);
         }
     }
 }

Member Data Documentation

float lar_content::TwoViewLongTracksTool::m_minMatchedFraction

private

The min matched sampling point fraction for particle creation.

Definition at line 70 of file TwoViewLongTracksTool.h.

unsigned int lar_content::TwoViewLongTracksTool::m_minMatchedSamplingPointRatio

private

The min ratio between 1st and 2nd highest msps for simple ambiguity resolution.

Definition at line 74 of file TwoViewLongTracksTool.h.

unsigned int lar_content::TwoViewLongTracksTool::m_minMatchedSamplingPoints

private

The min number of matched sampling points for particle creation.

Definition at line 72 of file TwoViewLongTracksTool.h.

float lar_content::TwoViewLongTracksTool::m_minMatchingScore

private

The min global matching score for particle creation.

Definition at line 71 of file TwoViewLongTracksTool.h.

float lar_content::TwoViewLongTracksTool::m_minXOverlapFraction

private

The min x overlap fraction (in each view) for particle creation.

Definition at line 73 of file TwoViewLongTracksTool.h.

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

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation