lar_content::MissingTrackTool Class Reference

MissingTrackTool class. More...

#include <MissingTrackTool.h>

Inheritance diagram for lar_content::MissingTrackTool:

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

Additional Inherited Members
Public Types inherited from lar_content::TransverseTensorTool
typedef ThreeViewTransverseTracksAlgorithm::MatchingType::TensorType	TensorType
typedef std::vector < TensorType::ElementList::const_iterator >	IteratorList

Detailed Description

MissingTrackTool class.

Definition at line 19 of file MissingTrackTool.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.

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

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

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

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.

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.

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.

---

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

• MissingTrackTool.h
• MissingTrackTool.cc