lar_content::NViewTrackMatchingAlgorithm< T > Class Template Reference

NViewTrackMatchingAlgorithm class. More...

#include "NViewTrackMatchingAlgorithm.h"

Inheritance diagram for lar_content::NViewTrackMatchingAlgorithm< T >:

Public Types
typedef T	MatchingType

Public Member Functions
	NViewTrackMatchingAlgorithm ()
	Default constructor. More...
virtual	~NViewTrackMatchingAlgorithm ()
	Destructor. More...
const TwoDSlidingFitResult &	GetCachedSlidingFitResult (const pandora::Cluster *const pCluster) const
	Get a sliding fit result from the algorithm cache. More...
unsigned int	GetSlidingFitWindow () const
	Get the layer window for the sliding linear fits. More...
virtual bool	MakeClusterSplits (const SplitPositionMap &splitPositionMap)
	Make cluster splits. More...
virtual bool	MakeClusterSplit (const pandora::CartesianVector &splitPosition, const pandora::Cluster *&pCurrentCluster, const pandora::Cluster *&pLowXCluster, const pandora::Cluster *&pHighXCluster) const
	Make a cluster split. More...
virtual void	UpdateForNewCluster (const pandora::Cluster *const pNewCluster)
	Update to reflect addition of a new cluster to the problem space. More...
virtual void	UpdateUponDeletion (const pandora::Cluster *const pDeletedCluster)
	Update to reflect cluster deletion. More...
virtual void	SelectInputClusters (const pandora::ClusterList *const pInputClusterList, pandora::ClusterList &selectedClusterList) const
	Select a subset of input clusters for processing in this algorithm. More...
virtual void	PrepareInputClusters (pandora::ClusterList &preparedClusterList)
	Perform any preparatory steps required on the input clusters, e.g. caching expensive fit results. More...
virtual void	SetPfoParticleId (PandoraContentApi::ParticleFlowObject::Parameters &pfoParameters) const
	Set pfo particle id. More...
const std::string &	GetClusterListName (const pandora::HitType hitType) const
	Get the cluster list name corresponding to a specified hit type. More...
const pandora::ClusterList &	GetInputClusterList (const pandora::HitType hitType) const
	Get the input cluster list corresponding to a specified hit type. More...
const pandora::ClusterList &	GetSelectedClusterList (const pandora::HitType hitType) const
	Get the selected cluster list corresponding to a specified hit type. More...
virtual void	CalculateOverlapResult (const pandora::Cluster *const pCluster1, const pandora::Cluster *const pCluster2, const pandora::Cluster *const pCluster3=nullptr)=0
	Calculate cluster overlap result and store in container. More...
virtual bool	MakeClusterMerges (const ClusterMergeMap &clusterMergeMap)
	Merge clusters together. More...
virtual bool	CreateThreeDParticles (const ProtoParticleVector &protoParticleVector)
	Create particles using findings from recent algorithm processing. More...
virtual void	SetPfoParameters (const ProtoParticle &protoParticle, PandoraContentApi::ParticleFlowObject::Parameters &pfoParameters) const
	Set Pfo properties. More...

Static Public Member Functions
static bool	SortSplitPositions (const pandora::CartesianVector &lhs, const pandora::CartesianVector &rhs)
	Sort split position cartesian vectors by increasing x coordinate. More...

Protected Member Functions
void	AddToSlidingFitCache (const pandora::Cluster *const pCluster)
	Add a new sliding fit result, for the specified cluster, to the algorithm cache. More...
void	RemoveFromSlidingFitCache (const pandora::Cluster *const pCluster)
	Remova an existing sliding fit result, for the specified cluster, from the algorithm cache. More...
virtual void	TidyUp ()
	Tidy member variables in derived class. More...
virtual pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)
MatchingType &	GetMatchingControl ()
	Get the matching control. More...
virtual void	SelectAllInputClusters ()
	Select a subset of input clusters for processing in this algorithm. More...
virtual void	PrepareAllInputClusters ()
	Perform any preparatory steps required, e.g. caching expensive fit results for clusters. More...
virtual void	PerformMainLoop ()
	Main loop over cluster combinations in order to populate the overlap container. Responsible for calling CalculateOverlapResult. More...
virtual void	ExamineOverlapContainer ()=0
	Examine contents of overlap container, collect together best-matching 2D particles and modify clusters as required. More...

Protected Attributes
MatchingType	m_matchingControl
	The matching control. More...

Private Attributes
unsigned int	m_slidingFitWindow
	The layer window for the sliding linear fits. More...
TwoDSlidingFitResultMap	m_slidingFitResultMap
	The sliding fit result map. More...
unsigned int	m_minClusterCaloHits
	The min number of hits in base cluster selection method. More...
float	m_minClusterLengthSquared
	The min length (squared) in base cluster selection method. More...

Detailed Description

template<typename T>
class lar_content::NViewTrackMatchingAlgorithm< T >

NViewTrackMatchingAlgorithm class.

Definition at line 26 of file NViewTrackMatchingAlgorithm.h.

Member Typedef Documentation

template<typename T>

typedef T lar_content::NViewMatchingAlgorithm< T >::MatchingType

inherited

Definition at line 23 of file NViewMatchingAlgorithm.h.

Constructor & Destructor Documentation

template<typename T >

lar_content::NViewTrackMatchingAlgorithm< T >::NViewTrackMatchingAlgorithm

(

)

Default constructor.

Definition at line 28 of file NViewTrackMatchingAlgorithm.cc.

                                                            :
    m_slidingFitWindow(20),
    m_minClusterCaloHits(5),
    m_minClusterLengthSquared(3.f * 3.f)
{
}

template<typename T >

lar_content::NViewTrackMatchingAlgorithm< T >::~NViewTrackMatchingAlgorithm ( )

virtual

Destructor.

Definition at line 38 of file NViewTrackMatchingAlgorithm.cc.

{
}

Member Function Documentation

template<typename T>

void lar_content::NViewTrackMatchingAlgorithm< T >::AddToSlidingFitCache ( const pandora::Cluster *const  pCluster )

protected

Add a new sliding fit result, for the specified cluster, to the algorithm cache.

Parameters

pCluster

address of the relevant cluster

Definition at line 255 of file NViewTrackMatchingAlgorithm.cc.

References lar_content::LArClusterHelper::GetClusterHitType(), lar_content::LArGeometryHelper::GetWirePitch(), lar_content::NViewTrackMatchingAlgorithm< T >::m_slidingFitResultMap, and lar_content::NViewTrackMatchingAlgorithm< T >::m_slidingFitWindow.

Referenced by lar_content::NViewTrackMatchingAlgorithm< T >::PrepareInputClusters(), and lar_content::NViewTrackMatchingAlgorithm< T >::UpdateForNewCluster().

{
    const float slidingFitPitch(LArGeometryHelper::GetWirePitch(this->GetPandora(), LArClusterHelper::GetClusterHitType(pCluster)));
    const TwoDSlidingFitResult slidingFitResult(pCluster, m_slidingFitWindow, slidingFitPitch);

    if (!m_slidingFitResultMap.insert(TwoDSlidingFitResultMap::value_type(pCluster, slidingFitResult)).second)
        throw StatusCodeException(STATUS_CODE_FAILURE);
}

virtual void lar_content::MatchingBaseAlgorithm::CalculateOverlapResult ( const pandora::Cluster *const  pCluster1, const pandora::Cluster *const  pCluster2, const pandora::Cluster *const  pCluster3 = nullptr  )

pure virtualinherited

Calculate cluster overlap result and store in container.

Parameters

pCluster1	address of cluster1
pCluster2	address of cluster2
pCluster3	address of cluster3

Implemented in lar_content::ThreeViewShowersAlgorithm, lar_content::TwoViewDeltaRayMatchingAlgorithm, lar_content::ThreeViewTrackFragmentsAlgorithm, lar_content::ThreeViewDeltaRayMatchingAlgorithm, lar_content::ThreeViewTransverseTracksAlgorithm, lar_content::TwoViewTransverseTracksAlgorithm, lar_content::ThreeViewLongitudinalTracksAlgorithm, and lar_content::ThreeViewRemnantsAlgorithm.

Referenced by lar_content::ThreeViewMatchingControl< T >::PerformMainLoop(), lar_content::TwoViewMatchingControl< T >::PerformMainLoop(), lar_content::ThreeViewMatchingControl< T >::UpdateForNewCluster(), and lar_content::TwoViewMatchingControl< T >::UpdateForNewCluster().

bool lar_content::MatchingBaseAlgorithm::CreateThreeDParticles ( const ProtoParticleVector &  protoParticleVector )

virtualinherited

Create particles using findings from recent algorithm processing.

Parameters

protoParticleVector	the proto particle vector
whether	particles were created

Definition at line 88 of file MatchingBaseAlgorithm.cc.

Referenced by lar_content::NViewDeltaRayMatchingAlgorithm< T >::CreatePfos(), lar_content::ClearRemnantsTool::CreateThreeDParticles(), lar_content::ClearLongitudinalTracksTool::CreateThreeDParticles(), lar_content::ClearTracksTool::CreateThreeDParticles(), lar_content::TwoViewClearTracksTool::CreateThreeDParticles(), lar_content::ClearTrackFragmentsTool::FindTrackFragments(), lar_content::MatchedEndPointsTool::Run(), lar_content::SimpleShowersTool::Run(), lar_content::ConnectedRemnantsTool::Run(), lar_content::MopUpRemnantsTool::Run(), lar_content::TwoViewSimpleTracksTool::Run(), lar_content::MissingTrackTool::Run(), lar_content::TracksCrossingGapsTool::Run(), lar_content::MissingTrackSegmentTool::Run(), lar_content::LongTracksTool::Run(), lar_content::TwoViewLongTracksTool::Run(), and lar_content::ClearShowersTool::Run().

{
    bool particlesMade(false);
    const PfoList *pPfoList(nullptr);
    std::string pfoListName;
    PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::CreateTemporaryListAndSetCurrent(*this, pPfoList, pfoListName));

    for (const ProtoParticle &protoParticle : protoParticleVector)
    {
        PandoraContentApi::ParticleFlowObject::Parameters pfoParameters;
        this->SetPfoParameters(protoParticle, pfoParameters);

        const ParticleFlowObject *pPfo(nullptr);
        PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::ParticleFlowObject::Create(*this, pfoParameters, pPfo));
        particlesMade = true;
    }

    if (!pPfoList->empty())
    {
        PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::SaveList<Pfo>(*this, m_outputPfoListName));
        PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::ReplaceCurrentList<Pfo>(*this, m_outputPfoListName));
    }

    return particlesMade;
}

virtual void lar_content::MatchingBaseAlgorithm::ExamineOverlapContainer ( )

protectedpure virtualinherited

Examine contents of overlap container, collect together best-matching 2D particles and modify clusters as required.

Implemented in lar_content::TwoViewDeltaRayMatchingAlgorithm, lar_content::ThreeViewShowersAlgorithm, lar_content::ThreeViewTrackFragmentsAlgorithm, lar_content::ThreeViewTransverseTracksAlgorithm, lar_content::TwoViewTransverseTracksAlgorithm, lar_content::ThreeViewLongitudinalTracksAlgorithm, lar_content::ThreeViewDeltaRayMatchingAlgorithm, and lar_content::ThreeViewRemnantsAlgorithm.

template<typename T>

const TwoDSlidingFitResult & lar_content::NViewTrackMatchingAlgorithm< T >::GetCachedSlidingFitResult

(

const pandora::Cluster *const

pCluster

)

const

Get a sliding fit result from the algorithm cache.

Parameters

pCluster

address of the relevant cluster

Definition at line 45 of file NViewTrackMatchingAlgorithm.cc.

References lar_content::NViewTrackMatchingAlgorithm< T >::m_slidingFitResultMap.

Referenced by lar_content::UndershootTracksTool::GetIteratorListModifications(), lar_content::OvershootTracksTool::GetIteratorListModifications(), lar_content::TwoViewThreeDKinkTool::GetIteratorListModifications(), lar_content::MissingTrackSegmentTool::GetSegmentOverlapMap(), lar_content::MissingTrackSegmentTool::GetSlidingFitResultMap(), lar_content::UndershootTracksTool::IsThreeDKink(), lar_content::OvershootTracksTool::IsThreeDKink(), lar_content::TwoViewThreeDKinkTool::IsThreeDKink(), lar_content::TrackSplittingTool::PassesChecks(), and lar_content::TracksCrossingGapsTool::PassesGapChecks().

{
    TwoDSlidingFitResultMap::const_iterator iter = m_slidingFitResultMap.find(pCluster);

    if (m_slidingFitResultMap.end() == iter)
        throw StatusCodeException(STATUS_CODE_NOT_FOUND);

    return iter->second;
}

template<typename T>

const std::string & lar_content::NViewMatchingAlgorithm< T >::GetClusterListName ( const pandora::HitType  hitType ) const

virtualinherited

Get the cluster list name corresponding to a specified hit type.

Parameters

hitType

the hit type

Returns

the cluster list name

Implements lar_content::MatchingBaseAlgorithm.

Definition at line 58 of file NViewMatchingAlgorithm.cc.

References lar_content::NViewMatchingAlgorithm< T >::m_matchingControl.

Referenced by lar_content::NViewDeltaRayMatchingAlgorithm< T >::AddInStrayClusters(), lar_content::NViewTrackMatchingAlgorithm< T >::MakeClusterSplits(), lar_content::DeltaRayMergeTool::MakeOneCommonViewMerges(), lar_content::DeltaRayMergeTool::MakeTwoCommonViewMerges(), lar_content::ClearTrackFragmentsTool::ProcessTensorElement(), lar_content::CosmicRayRemovalTool::ReclusterRemnant(), lar_content::TwoViewCosmicRayRemovalTool::ReclusterRemnant(), lar_content::CosmicRayRemovalTool::SplitDeltaRayCluster(), lar_content::TwoViewCosmicRayRemovalTool::SplitDeltaRayCluster(), and lar_content::DeltaRayRemovalTool::SplitMuonCluster().

{
    return m_matchingControl.GetClusterListName(hitType);
}

template<typename T>

const pandora::ClusterList & lar_content::NViewMatchingAlgorithm< T >::GetInputClusterList ( const pandora::HitType  hitType ) const

virtualinherited

Get the input cluster list corresponding to a specified hit type.

Parameters

hitType

the hit type

Returns

the input cluster list

Implements lar_content::MatchingBaseAlgorithm.

Definition at line 66 of file NViewMatchingAlgorithm.cc.

References lar_content::NViewMatchingAlgorithm< T >::m_matchingControl.

Referenced by lar_content::NViewDeltaRayMatchingAlgorithm< T >::FillStrayClusterList(), lar_content::MissingTrackSegmentTool::GetCandidateClusters(), lar_content::NViewDeltaRayMatchingAlgorithm< T >::PrepareInputClusters(), lar_content::ShowerTensorVisualizationTool::Run(), and lar_content::TransverseTensorVisualizationTool::Run().

{
    return m_matchingControl.GetInputClusterList(hitType);
}

template<typename T >

T & lar_content::NViewMatchingAlgorithm< T >::GetMatchingControl ( )

inlineprotectedinherited

Get the matching control.

Definition at line 59 of file NViewMatchingAlgorithm.h.

References lar_content::NViewMatchingAlgorithm< T >::m_matchingControl.

{
    return m_matchingControl;
}

template<typename T>

const pandora::ClusterList & lar_content::NViewMatchingAlgorithm< T >::GetSelectedClusterList ( const pandora::HitType  hitType ) const

virtualinherited

Get the selected cluster list corresponding to a specified hit type.

Parameters

hitType

the hit type

Returns

the selected cluster list

Implements lar_content::MatchingBaseAlgorithm.

Definition at line 74 of file NViewMatchingAlgorithm.cc.

References lar_content::NViewMatchingAlgorithm< T >::m_matchingControl.

{
    return m_matchingControl.GetSelectedClusterList(hitType);
}

template<typename T >

unsigned int lar_content::NViewTrackMatchingAlgorithm< T >::GetSlidingFitWindow ( ) const

inline

Get the layer window for the sliding linear fits.

Returns

the layer window for the sliding linear fits

Definition at line 118 of file NViewTrackMatchingAlgorithm.h.

References lar_content::NViewTrackMatchingAlgorithm< T >::m_slidingFitWindow.

Referenced by lar_content::MissingTrackSegmentTool::GetSlidingFitResultMap().

{
    return m_slidingFitWindow;
}

bool lar_content::MatchingBaseAlgorithm::MakeClusterMerges ( const ClusterMergeMap &  clusterMergeMap )

virtualinherited

Merge clusters together.

Parameters

clusterMergeMap

the cluster merge map

Returns

whether changes to the overlap container have been made

Definition at line 48 of file MatchingBaseAlgorithm.cc.

Referenced by lar_content::SplitShowersTool::ApplyChanges(), lar_content::ThreeDKinkBaseTool::ApplyChanges(), lar_content::TwoViewThreeDKinkTool::ApplyChanges(), lar_content::ConnectedRemnantsTool::Run(), and lar_content::MissingTrackSegmentTool::Run().

{
    ClusterSet deletedClusters;

    ClusterList parentClusters;
    for (const auto &mapEntry : clusterMergeMap)
        parentClusters.push_back(mapEntry.first);
    parentClusters.sort(LArClusterHelper::SortByNHits);

    for (const Cluster *const pParentCluster : parentClusters)
    {
        const HitType hitType(LArClusterHelper::GetClusterHitType(pParentCluster));
        const std::string &clusterListName(this->GetClusterListName(hitType));

        if (!((TPC_VIEW_U == hitType) || (TPC_VIEW_V == hitType) || (TPC_VIEW_W == hitType)))
            throw StatusCodeException(STATUS_CODE_FAILURE);

        ClusterList daughterClusters(clusterMergeMap.at(pParentCluster));
        daughterClusters.sort(LArClusterHelper::SortByNHits);

        for (const Cluster *const pDaughterCluster : daughterClusters)
        {
            if (deletedClusters.count(pParentCluster) || deletedClusters.count(pDaughterCluster))
                throw StatusCodeException(STATUS_CODE_FAILURE);

            this->UpdateUponDeletion(pDaughterCluster);
            this->UpdateUponDeletion(pParentCluster);
            PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=,
                PandoraContentApi::MergeAndDeleteClusters(*this, pParentCluster, pDaughterCluster, clusterListName, clusterListName));

            this->UpdateForNewCluster(pParentCluster);
            deletedClusters.insert(pDaughterCluster);
        }
    }

    return !(deletedClusters.empty());
}

template<typename T>

bool lar_content::NViewTrackMatchingAlgorithm< T >::MakeClusterSplit ( const pandora::CartesianVector &  splitPosition, const pandora::Cluster *&  pCurrentCluster, const pandora::Cluster *&  pLowXCluster, const pandora::Cluster *&  pHighXCluster  ) const

virtual

Make a cluster split.

Parameters

splitPosition	the split position
pCurrentCluster	the cluster to split
pLowXCluster	to receive the low x cluster
pHighXCluster	to receive the high x cluster

Returns

whether a cluster split occurred

Definition at line 105 of file NViewTrackMatchingAlgorithm.cc.

References f, lar_content::LArPointingCluster::GetInnerVertex(), lar_content::LArPointingCluster::GetOuterVertex(), and lar_content::LArPointingCluster::Vertex::GetPosition().

Referenced by lar_content::NViewTrackMatchingAlgorithm< T >::MakeClusterSplits().

{
    CartesianVector lowXEnd(0.f, 0.f, 0.f), highXEnd(0.f, 0.f, 0.f);

    try
    {
        LArPointingCluster pointingCluster(pCurrentCluster);
        const bool innerIsLowX(pointingCluster.GetInnerVertex().GetPosition().GetX() < pointingCluster.GetOuterVertex().GetPosition().GetX());
        lowXEnd = (innerIsLowX ? pointingCluster.GetInnerVertex().GetPosition() : pointingCluster.GetOuterVertex().GetPosition());
        highXEnd = (innerIsLowX ? pointingCluster.GetOuterVertex().GetPosition() : pointingCluster.GetInnerVertex().GetPosition());
    }
    catch (const StatusCodeException &)
    {
        return false;
    }

    const CartesianVector lowXUnitVector((lowXEnd - splitPosition).GetUnitVector());
    const CartesianVector highXUnitVector((highXEnd - splitPosition).GetUnitVector());

    CaloHitList caloHitList;
    pCurrentCluster->GetOrderedCaloHitList().FillCaloHitList(caloHitList);

    std::string originalListName, fragmentListName;
    const ClusterList clusterList(1, pCurrentCluster);
    PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::InitializeFragmentation(*this, clusterList, originalListName, fragmentListName));

    pLowXCluster = nullptr;
    pHighXCluster = nullptr;

    for (const CaloHit *const pCaloHit : caloHitList)
    {
        const CartesianVector unitVector((pCaloHit->GetPositionVector() - splitPosition).GetUnitVector());
        const float dotProductLowX(unitVector.GetDotProduct(lowXUnitVector));
        const float dotProductHighX(unitVector.GetDotProduct(highXUnitVector));

        const Cluster *&pClusterToModify((dotProductLowX > dotProductHighX) ? pLowXCluster : pHighXCluster);

        if (!pClusterToModify)
        {
            PandoraContentApi::Cluster::Parameters parameters;
            parameters.m_caloHitList.push_back(pCaloHit);
            PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::Cluster::Create(*this, parameters, pClusterToModify));
        }
        else
        {
            PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::AddToCluster(*this, pClusterToModify, pCaloHit));
        }
    }

    if (!pLowXCluster || !pHighXCluster)
    {
        PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::EndFragmentation(*this, originalListName, fragmentListName));
        return false;
    }

    PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::EndFragmentation(*this, fragmentListName, originalListName));
    return true;
}

template<typename T >

bool lar_content::NViewTrackMatchingAlgorithm< T >::MakeClusterSplits ( const SplitPositionMap &  splitPositionMap )

virtual

Make cluster splits.

Parameters

splitPositionMap

the split position map

Returns

whether changes to the overlap container have been made

Definition at line 58 of file NViewTrackMatchingAlgorithm.cc.

References lar_content::LArClusterHelper::GetClusterHitType(), lar_content::NViewMatchingAlgorithm< T >::GetClusterListName(), lar_content::NViewTrackMatchingAlgorithm< T >::MakeClusterSplit(), lar_content::LArClusterHelper::SortByNHits(), lar_content::NViewTrackMatchingAlgorithm< T >::UpdateForNewCluster(), and lar_content::NViewTrackMatchingAlgorithm< T >::UpdateUponDeletion().

Referenced by lar_content::ThreeDKinkBaseTool::ApplyChanges(), lar_content::TwoViewThreeDKinkTool::ApplyChanges(), and lar_content::TrackSplittingTool::Run().

{
    bool changesMade(false);

    ClusterList splitClusters;
    for (const auto &mapEntry : splitPositionMap)
        splitClusters.push_back(mapEntry.first);
    splitClusters.sort(LArClusterHelper::SortByNHits);

    for (const Cluster *pCurrentCluster : splitClusters)
    {
        CartesianPointVector splitPositions(splitPositionMap.at(pCurrentCluster));
        std::sort(splitPositions.begin(), splitPositions.end(), NViewTrackMatchingAlgorithm<T>::SortSplitPositions);

        const HitType hitType(LArClusterHelper::GetClusterHitType(pCurrentCluster));
        const std::string &clusterListName(this->GetClusterListName(hitType));

        if (!((TPC_VIEW_U == hitType) || (TPC_VIEW_V == hitType) || (TPC_VIEW_W == hitType)))
            throw StatusCodeException(STATUS_CODE_FAILURE);

        PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::ReplaceCurrentList<Cluster>(*this, clusterListName));

        for (const CartesianVector &splitPosition : splitPositions)
        {
            const Cluster *pLowXCluster(nullptr), *pHighXCluster(nullptr);
            this->UpdateUponDeletion(pCurrentCluster);

            if (this->MakeClusterSplit(splitPosition, pCurrentCluster, pLowXCluster, pHighXCluster))
            {
                changesMade = true;
                this->UpdateForNewCluster(pLowXCluster);
                this->UpdateForNewCluster(pHighXCluster);
                pCurrentCluster = pHighXCluster;
            }
            else
            {
                this->UpdateForNewCluster(pCurrentCluster);
            }
        }
    }

    return changesMade;
}

template<typename T >

void lar_content::NViewMatchingAlgorithm< T >::PerformMainLoop ( )

protectedvirtualinherited

Main loop over cluster combinations in order to populate the overlap container. Responsible for calling CalculateOverlapResult.

Implements lar_content::MatchingBaseAlgorithm.

Reimplemented in lar_content::ThreeViewTrackFragmentsAlgorithm.

Definition at line 106 of file NViewMatchingAlgorithm.cc.

References lar_content::NViewMatchingAlgorithm< T >::m_matchingControl.

{
    m_matchingControl.PerformMainLoop();
}

template<typename T >

void lar_content::NViewMatchingAlgorithm< T >::PrepareAllInputClusters ( )

protectedvirtualinherited

Perform any preparatory steps required, e.g. caching expensive fit results for clusters.

Implements lar_content::MatchingBaseAlgorithm.

Definition at line 90 of file NViewMatchingAlgorithm.cc.

References lar_content::NViewMatchingAlgorithm< T >::m_matchingControl.

{
    m_matchingControl.PrepareAllInputClusters();
}

template<typename T>

void lar_content::NViewTrackMatchingAlgorithm< T >::PrepareInputClusters ( pandora::ClusterList &  preparedClusterList )

virtual

Perform any preparatory steps required on the input clusters, e.g. caching expensive fit results.

Parameters

preparedClusterList

to receive the prepared cluster list, non const so as to be able to modify input selected list

Reimplemented from lar_content::MatchingBaseAlgorithm.

Definition at line 225 of file NViewTrackMatchingAlgorithm.cc.

References lar_content::NViewTrackMatchingAlgorithm< T >::AddToSlidingFitCache().

{
    for (ClusterList::iterator iter = preparedClusterList.begin(), iterEnd = preparedClusterList.end(); iter != iterEnd;)
    {
        try
        {
            this->AddToSlidingFitCache(*iter);
            ++iter;
        }
        catch (const StatusCodeException &statusCodeException)
        {
            preparedClusterList.erase(iter++);

            if (STATUS_CODE_FAILURE == statusCodeException.GetStatusCode())
                throw statusCodeException;
        }
    }
}

template<typename T>

StatusCode lar_content::NViewTrackMatchingAlgorithm< T >::ReadSettings ( const pandora::TiXmlHandle  xmlHandle )

protectedvirtual

Reimplemented from lar_content::NViewMatchingAlgorithm< T >.

Reimplemented in lar_content::ThreeViewTrackFragmentsAlgorithm, lar_content::ThreeViewTransverseTracksAlgorithm, lar_content::TwoViewTransverseTracksAlgorithm, and lar_content::ThreeViewLongitudinalTracksAlgorithm.

Definition at line 287 of file NViewTrackMatchingAlgorithm.cc.

References lar_content::NViewTrackMatchingAlgorithm< T >::m_minClusterCaloHits, lar_content::NViewTrackMatchingAlgorithm< T >::m_minClusterLengthSquared, lar_content::NViewTrackMatchingAlgorithm< T >::m_slidingFitWindow, and lar_content::NViewMatchingAlgorithm< T >::ReadSettings().

{
    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "SlidingFitWindow", m_slidingFitWindow));

    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MinClusterCaloHits", m_minClusterCaloHits));

    float minClusterLength = std::sqrt(m_minClusterLengthSquared);
    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MinClusterLength", minClusterLength));
    m_minClusterLengthSquared = minClusterLength * minClusterLength;

    return NViewMatchingAlgorithm<T>::ReadSettings(xmlHandle);
}

template<typename T>

void lar_content::NViewTrackMatchingAlgorithm< T >::RemoveFromSlidingFitCache ( const pandora::Cluster *const  pCluster )

protected

Remova an existing sliding fit result, for the specified cluster, from the algorithm cache.

Parameters

pCluster

address of the relevant cluster

Definition at line 267 of file NViewTrackMatchingAlgorithm.cc.

References lar_content::NViewTrackMatchingAlgorithm< T >::m_slidingFitResultMap.

Referenced by lar_content::NViewTrackMatchingAlgorithm< T >::UpdateUponDeletion().

{
    TwoDSlidingFitResultMap::iterator iter = m_slidingFitResultMap.find(pCluster);

    if (m_slidingFitResultMap.end() != iter)
        m_slidingFitResultMap.erase(iter);
}

template<typename T >

void lar_content::NViewMatchingAlgorithm< T >::SelectAllInputClusters ( )

protectedvirtualinherited

Select a subset of input clusters for processing in this algorithm.

Implements lar_content::MatchingBaseAlgorithm.

Definition at line 82 of file NViewMatchingAlgorithm.cc.

References lar_content::NViewMatchingAlgorithm< T >::m_matchingControl.

{
    m_matchingControl.SelectAllInputClusters();
}

template<typename T>

void lar_content::NViewTrackMatchingAlgorithm< T >::SelectInputClusters ( const pandora::ClusterList *const  pInputClusterList, pandora::ClusterList &  selectedClusterList  ) const

virtual

Select a subset of input clusters for processing in this algorithm.

Parameters

pInputClusterList	address of an input cluster list
selectedClusterList	to receive the selected cluster list

Reimplemented from lar_content::MatchingBaseAlgorithm.

Definition at line 205 of file NViewTrackMatchingAlgorithm.cc.

References lar_content::LArClusterHelper::GetLengthSquared(), lar_content::NViewTrackMatchingAlgorithm< T >::m_minClusterCaloHits, and lar_content::NViewTrackMatchingAlgorithm< T >::m_minClusterLengthSquared.

Referenced by lar_content::ClearTrackFragmentsTool::ProcessTensorElement().

{
    for (const Cluster *const pCluster : *pInputClusterList)
    {
        if (!pCluster->IsAvailable())
            continue;

        if (pCluster->GetNCaloHits() < m_minClusterCaloHits)
            continue;

        if (LArClusterHelper::GetLengthSquared(pCluster) < m_minClusterLengthSquared)
            continue;

        selectedClusterList.push_back(pCluster);
    }
}

void lar_content::MatchingBaseAlgorithm::SetPfoParameters ( const ProtoParticle &  protoParticle, PandoraContentApi::ParticleFlowObject::Parameters &  pfoParameters  ) const

virtualinherited

Set Pfo properties.

Parameters

protoParticle	the input proto particle
pfoParameters	the output pfo parameters

Definition at line 116 of file MatchingBaseAlgorithm.cc.

References f, and lar_content::ProtoParticle::m_clusterList.

{
    this->SetPfoParticleId(pfoParameters);
    pfoParameters.m_charge = PdgTable::GetParticleCharge(pfoParameters.m_particleId.Get());
    pfoParameters.m_mass = PdgTable::GetParticleMass(pfoParameters.m_particleId.Get());
    pfoParameters.m_energy = 0.f;
    pfoParameters.m_momentum = CartesianVector(0.f, 0.f, 0.f);
    pfoParameters.m_clusterList = protoParticle.m_clusterList;
}

template<typename T >

void lar_content::NViewTrackMatchingAlgorithm< T >::SetPfoParticleId ( PandoraContentApi::ParticleFlowObject::Parameters &  pfoParameters ) const

virtual

Set pfo particle id.

Parameters

pfoParameters

the output pfo parameters

Reimplemented from lar_content::MatchingBaseAlgorithm.

Definition at line 247 of file NViewTrackMatchingAlgorithm.cc.

{
    pfoParameters.m_particleId = MU_MINUS; // Track
}

template<typename T >

bool lar_content::NViewTrackMatchingAlgorithm< T >::SortSplitPositions ( const pandora::CartesianVector &  lhs, const pandora::CartesianVector &  rhs  )

static

Sort split position cartesian vectors by increasing x coordinate.

Parameters

lhs	the first cartesian vector
rhs	the second cartesian vector

Definition at line 168 of file NViewTrackMatchingAlgorithm.cc.

{
    return (lhs.GetX() < rhs.GetX());
}

template<typename T >

void lar_content::NViewTrackMatchingAlgorithm< T >::TidyUp ( )

protectedvirtual

Tidy member variables in derived class.

Reimplemented from lar_content::NViewMatchingAlgorithm< T >.

Definition at line 278 of file NViewTrackMatchingAlgorithm.cc.

References lar_content::NViewTrackMatchingAlgorithm< T >::m_slidingFitResultMap, and lar_content::NViewMatchingAlgorithm< T >::TidyUp().

{
    m_slidingFitResultMap.clear();
    return NViewMatchingAlgorithm<T>::TidyUp();
}

template<typename T>

void lar_content::NViewTrackMatchingAlgorithm< T >::UpdateForNewCluster ( const pandora::Cluster *const  pNewCluster )

virtual

Update to reflect addition of a new cluster to the problem space.

Parameters

pNewCluster

address of the new cluster

Reimplemented from lar_content::NViewMatchingAlgorithm< T >.

Reimplemented in lar_content::ThreeViewTrackFragmentsAlgorithm.

Definition at line 176 of file NViewTrackMatchingAlgorithm.cc.

References lar_content::NViewTrackMatchingAlgorithm< T >::AddToSlidingFitCache(), and lar_content::NViewMatchingAlgorithm< T >::UpdateForNewCluster().

Referenced by lar_content::NViewTrackMatchingAlgorithm< T >::MakeClusterSplits().

{
    try
    {
        this->AddToSlidingFitCache(pNewCluster);
    }
    catch (const StatusCodeException &statusCodeException)
    {
        if (STATUS_CODE_FAILURE == statusCodeException.GetStatusCode())
            throw statusCodeException;

        return;
    }

    NViewMatchingAlgorithm<T>::UpdateForNewCluster(pNewCluster);
}

template<typename T>

void lar_content::NViewTrackMatchingAlgorithm< T >::UpdateUponDeletion ( const pandora::Cluster *const  pDeletedCluster )

virtual

Update to reflect cluster deletion.

Parameters

pDeletedCluster

address of the deleted cluster

Reimplemented from lar_content::NViewMatchingAlgorithm< T >.

Definition at line 196 of file NViewTrackMatchingAlgorithm.cc.

References lar_content::NViewTrackMatchingAlgorithm< T >::RemoveFromSlidingFitCache(), and lar_content::NViewMatchingAlgorithm< T >::UpdateUponDeletion().

Referenced by lar_content::ClearTrackFragmentsTool::FindTrackFragments(), lar_content::NViewTrackMatchingAlgorithm< T >::MakeClusterSplits(), and lar_content::ClearTrackFragmentsTool::ProcessTensorElement().

{
    this->RemoveFromSlidingFitCache(pDeletedCluster);
    NViewMatchingAlgorithm<T>::UpdateUponDeletion(pDeletedCluster);
}

Member Data Documentation

template<typename T>

MatchingType lar_content::NViewMatchingAlgorithm< T >::m_matchingControl

protectedinherited

The matching control.

Definition at line 53 of file NViewMatchingAlgorithm.h.

Referenced by lar_content::NViewMatchingAlgorithm< T >::GetClusterListName(), lar_content::NViewMatchingAlgorithm< T >::GetInputClusterList(), lar_content::NViewMatchingAlgorithm< T >::GetMatchingControl(), lar_content::NViewMatchingAlgorithm< T >::GetSelectedClusterList(), lar_content::NViewMatchingAlgorithm< T >::PerformMainLoop(), lar_content::NViewMatchingAlgorithm< T >::PrepareAllInputClusters(), lar_content::NViewMatchingAlgorithm< T >::ReadSettings(), lar_content::NViewMatchingAlgorithm< T >::SelectAllInputClusters(), lar_content::NViewMatchingAlgorithm< T >::TidyUp(), lar_content::NViewMatchingAlgorithm< T >::UpdateForNewCluster(), and lar_content::NViewMatchingAlgorithm< T >::UpdateUponDeletion().

template<typename T>

unsigned int lar_content::NViewTrackMatchingAlgorithm< T >::m_minClusterCaloHits

private

The min number of hits in base cluster selection method.

Definition at line 111 of file NViewTrackMatchingAlgorithm.h.

Referenced by lar_content::NViewTrackMatchingAlgorithm< T >::ReadSettings(), and lar_content::NViewTrackMatchingAlgorithm< T >::SelectInputClusters().

template<typename T>

float lar_content::NViewTrackMatchingAlgorithm< T >::m_minClusterLengthSquared

private

The min length (squared) in base cluster selection method.

Definition at line 112 of file NViewTrackMatchingAlgorithm.h.

Referenced by lar_content::NViewTrackMatchingAlgorithm< T >::ReadSettings(), and lar_content::NViewTrackMatchingAlgorithm< T >::SelectInputClusters().

template<typename T>

TwoDSlidingFitResultMap lar_content::NViewTrackMatchingAlgorithm< T >::m_slidingFitResultMap

private

The sliding fit result map.

Definition at line 109 of file NViewTrackMatchingAlgorithm.h.

Referenced by lar_content::NViewTrackMatchingAlgorithm< T >::AddToSlidingFitCache(), lar_content::NViewTrackMatchingAlgorithm< T >::GetCachedSlidingFitResult(), lar_content::NViewTrackMatchingAlgorithm< T >::RemoveFromSlidingFitCache(), and lar_content::NViewTrackMatchingAlgorithm< T >::TidyUp().

template<typename T>

unsigned int lar_content::NViewTrackMatchingAlgorithm< T >::m_slidingFitWindow

private

The layer window for the sliding linear fits.

Definition at line 108 of file NViewTrackMatchingAlgorithm.h.

Referenced by lar_content::NViewTrackMatchingAlgorithm< T >::AddToSlidingFitCache(), lar_content::NViewTrackMatchingAlgorithm< T >::GetSlidingFitWindow(), and lar_content::NViewTrackMatchingAlgorithm< T >::ReadSettings().

---

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

• larpandoracontent/v04_08_01/source/larpandoracontent/LArThreeDReco/LArThreeDBase/NViewTrackMatchingAlgorithm.h
• larpandoracontent/v04_08_01/source/larpandoracontent/LArThreeDReco/LArThreeDBase/NViewTrackMatchingAlgorithm.cc