#include "MissingTrackTool.h"

Inheritance diagram for lar_content::MissingTrackTool:

Public Types
typedef ThreeViewTransverseTracksAlgorithm::MatchingType::TensorType	TensorType

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

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

bool	Run (ThreeViewTransverseTracksAlgorithm *const pAlgorithm, TensorType &overlapTensor)
	Run the algorithm tool. More...

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

void	FindMissingTracks (const TensorType &overlapTensor, ProtoParticleVector &protoParticleVector) const
	Find missing tracks, due to merging of multiple particle deposits into single hits during hit creation. More...

Private Attributes
unsigned int	m_minMatchedSamplingPoints
	The min number of matched sampling points for the unavailable tensor element. More...

float	m_minMatchedFraction
	The min matched sampling point fraction for the unavailable tensor element. More...

float	m_maxReducedChiSquared
	The max reduced chi squared value for the unavailable tensor element. More...

float	m_minXOverlapFraction
	The min x overlap fraction for the two available clusters in the tensor element. More...

Detailed Description

MissingTrackTool class.

Definition at line 19 of file MissingTrackTool.h.

Member Typedef Documentation

typedef std::vector<TensorType::ElementList::const_iterator> lar_content::TransverseTensorTool::IteratorList

inherited

Definition at line 131 of file ThreeViewTransverseTracksAlgorithm.h.

typedef ThreeViewTransverseTracksAlgorithm::MatchingType::TensorType lar_content::TransverseTensorTool::TensorType

inherited

Definition at line 130 of file ThreeViewTransverseTracksAlgorithm.h.

Constructor & Destructor Documentation

lar_content::MissingTrackTool::MissingTrackTool ( )

Default constructor.

Definition at line 17 of file MissingTrackTool.cc.

                                    :
     m_minMatchedSamplingPoints(15),
     m_minMatchedFraction(0.95f),
     m_maxReducedChiSquared(0.707f),
     m_minXOverlapFraction(0.75f)
 {
 }

Member Function Documentation

void lar_content::MissingTrackTool::FindMissingTracks	(	const TensorType &	overlapTensor,
		ProtoParticleVector &	protoParticleVector
	)		const

private

Find missing tracks, due to merging of multiple particle deposits into single hits during hit creation.

Parameters

overlapTensor	the overlap tensor
protoParticleVector	to receive the list of proto particles

Definition at line 41 of file MissingTrackTool.cc.

References lar_content::OverlapTensor< T >::GetConnectedElements(), lar_content::OverlapTensor< T >::GetSortedKeyClusters(), lar_content::ProtoParticle::m_clusterList, m_maxReducedChiSquared, m_minMatchedFraction, m_minMatchedSamplingPoints, and m_minXOverlapFraction.

Referenced by Run().

 {
     ClusterSet usedClusters;
     ClusterVector sortedKeyClusters;
     overlapTensor.GetSortedKeyClusters(sortedKeyClusters);
 
     for (const Cluster *const pKeyCluster : sortedKeyClusters)
     {
         unsigned int nU(0), nV(0), nW(0);
         TensorType::ElementList elementList;
         overlapTensor.GetConnectedElements(pKeyCluster, false, elementList, nU, nV, nW);
 
         for (TensorType::ElementList::const_iterator eIter = elementList.begin(); eIter != elementList.end(); ++eIter)
         {
             const bool includeU(eIter->GetClusterU()->IsAvailable() && !usedClusters.count(eIter->GetClusterU()));
             const bool includeV(eIter->GetClusterV()->IsAvailable() && !usedClusters.count(eIter->GetClusterV()));
             const bool includeW(eIter->GetClusterW()->IsAvailable() && !usedClusters.count(eIter->GetClusterW()));
 
             unsigned int nAvailable(0);
             if (includeU)
                 ++nAvailable;
             if (includeV)
                 ++nAvailable;
             if (includeW)
                 ++nAvailable;
 
             if (2 != nAvailable)
                 continue;
 
             const TransverseOverlapResult &overlapResult(eIter->GetOverlapResult());
 
             if (overlapResult.GetNMatchedSamplingPoints() < m_minMatchedSamplingPoints)
                 continue;
 
             if (overlapResult.GetMatchedFraction() < m_minMatchedFraction)
                 continue;
 
             if (overlapResult.GetReducedChi2() > m_maxReducedChiSquared)
                 continue;
 
             if ((overlapResult.GetXOverlap().GetXSpanU() < std::numeric_limits<float>::epsilon()) ||
                 (overlapResult.GetXOverlap().GetXSpanV() < std::numeric_limits<float>::epsilon()) ||
                 (overlapResult.GetXOverlap().GetXSpanW() < std::numeric_limits<float>::epsilon()))
             {
                 continue;
             }
 
             const float xOverlapSpan(overlapResult.GetXOverlap().GetXOverlapSpan());
 
             if (includeU && (xOverlapSpan / overlapResult.GetXOverlap().GetXSpanU() < m_minXOverlapFraction))
                 continue;
 
             if (includeV && (xOverlapSpan / overlapResult.GetXOverlap().GetXSpanV() < m_minXOverlapFraction))
                 continue;
 
             if (includeW && (xOverlapSpan / overlapResult.GetXOverlap().GetXSpanW() < m_minXOverlapFraction))
                 continue;
 
             ProtoParticle protoParticle;
             if (includeU)
                 protoParticle.m_clusterList.push_back(eIter->GetClusterU());
             if (includeV)
                 protoParticle.m_clusterList.push_back(eIter->GetClusterV());
             if (includeW)
                 protoParticle.m_clusterList.push_back(eIter->GetClusterW());
 
             protoParticleVector.push_back(protoParticle);
             usedClusters.insert(eIter->GetClusterU());
             usedClusters.insert(eIter->GetClusterV());
             usedClusters.insert(eIter->GetClusterW());
         }
     }
 }

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

private

Definition at line 117 of file MissingTrackTool.cc.

References m_maxReducedChiSquared, m_minMatchedFraction, m_minMatchedSamplingPoints, and m_minXOverlapFraction.

 {
     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, "MinMatchedFraction", m_minMatchedFraction));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MaxReducedChiSquared", m_maxReducedChiSquared));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MinXOverlapFraction", m_minXOverlapFraction));
 
     return STATUS_CODE_SUCCESS;
 }

bool lar_content::MissingTrackTool::Run	(	ThreeViewTransverseTracksAlgorithm *const	pAlgorithm,
		TensorType &	overlapTensor
	)

virtual

Run the algorithm tool.

Parameters

pAlgorithm	address of the calling algorithm
overlapTensor	the overlap tensor

Returns: whether changes have been made by the tool

Implements lar_content::TransverseTensorTool.

Definition at line 27 of file MissingTrackTool.cc.

References lar_content::MatchingBaseAlgorithm::CreateThreeDParticles(), and FindMissingTracks().

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

Member Data Documentation

float lar_content::MissingTrackTool::m_maxReducedChiSquared

private

The max reduced chi squared value for the unavailable tensor element.

Definition at line 42 of file MissingTrackTool.h.

Referenced by FindMissingTracks(), and ReadSettings().

float lar_content::MissingTrackTool::m_minMatchedFraction

private

The min matched sampling point fraction for the unavailable tensor element.

Definition at line 41 of file MissingTrackTool.h.

Referenced by FindMissingTracks(), and ReadSettings().

unsigned int lar_content::MissingTrackTool::m_minMatchedSamplingPoints

private

The min number of matched sampling points for the unavailable tensor element.

Definition at line 40 of file MissingTrackTool.h.

Referenced by FindMissingTracks(), and ReadSettings().

float lar_content::MissingTrackTool::m_minXOverlapFraction

private

The min x overlap fraction for the two available clusters in the tensor element.

Definition at line 43 of file MissingTrackTool.h.

Referenced by FindMissingTracks(), and ReadSettings().

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

larpandoracontent/v04_08_01/source/larpandoracontent/LArThreeDReco/LArTransverseTrackMatching/MissingTrackTool.h
larpandoracontent/v04_08_01/source/larpandoracontent/LArThreeDReco/LArTransverseTrackMatching/MissingTrackTool.cc

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation