NearbyClusterMopUpAlgorithm class. More...

#include "NearbyClusterMopUpAlgorithm.h"

Inheritance diagram for lar_content::NearbyClusterMopUpAlgorithm:

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

template<typename T >
const std::string	GetListName (const T *const pT) const
	Find the name of the list hosting a specific object. More...

Protected Types
typedef std::unordered_map< const pandora::Cluster *, float >	AssociationDetails

typedef std::unordered_map< const pandora::Cluster *, AssociationDetails >	ClusterAssociationMap

Protected Member Functions
virtual pandora::StatusCode	Run ()

virtual void	GetPfoClusterLists (pandora::ClusterList &clusterListU, pandora::ClusterList &clusterListV, pandora::ClusterList &clusterListW) const
	Get the two dimensional clusters contained in the input pfo list, divided into three different views. More...

virtual void	GetDaughterClusterLists (pandora::ClusterList &clusterListU, pandora::ClusterList &clusterListV, pandora::ClusterList &clusterListW) const
	Get the two dimensional clusters contained in the input remant/daughter cluster lists, divided into three different views. More...

virtual void	GetClusterLists (const pandora::ClusterList &inputClusterList, const bool availabilityFlag, pandora::ClusterList &clusterListU, pandora::ClusterList &clusterListV, pandora::ClusterList &clusterListW) const
	Get the two dimensional clusters contained in an input cluster list, divided into three different views. More...

virtual void	MakeClusterMerges (const ClusterAssociationMap &clusterAssociationMap) const
	Make the cluster merges specified in the cluster association map, using list name information in the cluster list name map. More...

Protected Attributes
pandora::StringVector	m_pfoListNames
	The list of pfo list names. More...

bool	m_excludePfosContainingTracks
	Whether to exclude any pfos containing clusters flagged as fixed tracks. More...

pandora::StringVector	m_daughterListNames
	The list of potential daughter object list names. More...

Private Member Functions
void	ClusterMopUp (const pandora::ClusterList &pfoClusters, const pandora::ClusterList &remnantClusters) const
	Cluster mop up for a single view. This function is responsible for instructing pandora to make cluster alterations. More...

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

Private Attributes
unsigned int	m_minHitsInCluster
	Minimum number of hits in order to consider a cluster. More...

float	m_vertexProximity
	Distance between cluster inner/outer centroid and vtx to declare cluster vtx associated. More...

float	m_minClusterSeparation
	Minimum distance between parent and daughter clusters to declare clusters associated. More...

float	m_touchingDistance
	Threshold (small) distance below which parent and daughter clusters are declated touching. More...

Detailed Description

NearbyClusterMopUpAlgorithm class.

Definition at line 22 of file NearbyClusterMopUpAlgorithm.h.

Member Typedef Documentation

typedef std::unordered_map<const pandora::Cluster*, float> lar_content::ClusterMopUpBaseAlgorithm::AssociationDetails

protectedinherited

Definition at line 71 of file ClusterMopUpBaseAlgorithm.h.

typedef std::unordered_map<const pandora::Cluster*, AssociationDetails> lar_content::ClusterMopUpBaseAlgorithm::ClusterAssociationMap

protectedinherited

Definition at line 72 of file ClusterMopUpBaseAlgorithm.h.

Constructor & Destructor Documentation

lar_content::NearbyClusterMopUpAlgorithm::NearbyClusterMopUpAlgorithm ( )

Default constructor.

Definition at line 21 of file NearbyClusterMopUpAlgorithm.cc.

                                                          :
     m_minHitsInCluster(5),
     m_vertexProximity(5.f),
     m_minClusterSeparation(2.5f),
     m_touchingDistance(0.001f)
 {
 }

Member Function Documentation

void lar_content::NearbyClusterMopUpAlgorithm::ClusterMopUp	(	const pandora::ClusterList &	pfoClusters,
		const pandora::ClusterList &	remnantClusters
	)		const

privatevirtual

Cluster mop up for a single view. This function is responsible for instructing pandora to make cluster alterations.

Parameters

pfoClusters	the list of pfo clusters
remnantClusters	the list of remnant clusters

Implements lar_content::ClusterMopUpBaseAlgorithm.

Definition at line 31 of file NearbyClusterMopUpAlgorithm.cc.

References f, lar_content::LArClusterHelper::GetClosestDistance(), lar_content::LArClusterHelper::GetClusterHitType(), m_minClusterSeparation, m_minHitsInCluster, m_touchingDistance, m_vertexProximity, lar_content::ClusterMopUpBaseAlgorithm::MakeClusterMerges(), min, lar_content::LArGeometryHelper::ProjectPosition(), and lar_content::LArClusterHelper::SortByNHits().

 {
     ClusterAssociationMap clusterAssociationMap;
 
     const VertexList *pVertexList(NULL);
     PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetCurrentList(*this, pVertexList));
     const Vertex *const pVertex(((pVertexList->size() == 1) && (VERTEX_3D == (*(pVertexList->begin()))->GetVertexType())) ? *(pVertexList->begin()) : NULL);
 
     ClusterVector sortedPfoClusters(pfoClusters.begin(), pfoClusters.end());
     std::sort(sortedPfoClusters.begin(), sortedPfoClusters.end(), LArClusterHelper::SortByNHits);
 
     ClusterVector sortedRemnantClusters(remnantClusters.begin(), remnantClusters.end());
     std::sort(sortedRemnantClusters.begin(), sortedRemnantClusters.end(), LArClusterHelper::SortByNHits);
 
     for (const Cluster *const pClusterP : sortedPfoClusters)
     {
         const HitType hitType(LArClusterHelper::GetClusterHitType(pClusterP));
         const CartesianVector vertexPosition2D(!pVertex ? CartesianVector(0.f, 0.f, 0.f) :
             LArGeometryHelper::ProjectPosition(this->GetPandora(), pVertex->GetPosition(), hitType));
 
         const float innerPV((vertexPosition2D - pClusterP->GetCentroid(pClusterP->GetInnerPseudoLayer())).GetMagnitude());
         const float outerPV((vertexPosition2D - pClusterP->GetCentroid(pClusterP->GetOuterPseudoLayer())).GetMagnitude());
 
         for (const Cluster *const pClusterR : sortedRemnantClusters)
         {
             if (pClusterR->GetNCaloHits() < m_minHitsInCluster)
                 continue;
 
             const float innerRV((vertexPosition2D - pClusterR->GetCentroid(pClusterR->GetInnerPseudoLayer())).GetMagnitude());
             const float outerRV((vertexPosition2D - pClusterR->GetCentroid(pClusterR->GetOuterPseudoLayer())).GetMagnitude());
 
             // ATTN Could use pointing clusters here, for consistency with other vertex association mechanics
             if (pVertex && (((innerPV < m_vertexProximity) || (outerPV < m_vertexProximity)) && ((innerRV < m_vertexProximity) || (outerRV < m_vertexProximity))))
                 continue;
 
             const float innerRP(LArClusterHelper::GetClosestDistance(pClusterR->GetCentroid(pClusterR->GetInnerPseudoLayer()), pClusterP));
             const float outerRP(LArClusterHelper::GetClosestDistance(pClusterR->GetCentroid(pClusterR->GetOuterPseudoLayer()), pClusterP));
 
             const float minSeparation(std::min(innerRP, outerRP));
 
             if (minSeparation > m_minClusterSeparation)
                 continue;
 
             // ATTN Use of ClusterMopUp base algorithm assumes bigger figure of merit means better association
             if (m_touchingDistance < std::numeric_limits<float>::epsilon())
                 throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);
 
             AssociationDetails &associationDetails(clusterAssociationMap[pClusterR]);
             const float figureOfMerit((minSeparation < m_touchingDistance) ? 1.f / m_touchingDistance : 1.f / minSeparation);
 
             if (!associationDetails.insert(AssociationDetails::value_type(pClusterP, figureOfMerit)).second)
                 throw StatusCodeException(STATUS_CODE_ALREADY_PRESENT);
         }
     }
 
     this->MakeClusterMerges(clusterAssociationMap);
 }

void lar_content::ClusterMopUpBaseAlgorithm::GetClusterLists	(	const pandora::ClusterList &	inputClusterList,
		const bool	availabilityFlag,
		pandora::ClusterList &	clusterListU,
		pandora::ClusterList &	clusterListV,
		pandora::ClusterList &	clusterListW
	)		const

protectedvirtualinherited

Get the two dimensional clusters contained in an input cluster list, divided into three different views.

Parameters

inputClusterList	the input cluster list
availabilityFlag	only clusters with matching availability will be considered
clusterListU	to receive the list of clusters in the u view
clusterListV	to receive the list of clusters in the v view
clusterListW	to receive the list of clusters in the w view

Definition at line 82 of file ClusterMopUpBaseAlgorithm.cc.

References lar_content::LArClusterHelper::GetClusterHitType(), and target.

Referenced by lar_content::ClusterMopUpBaseAlgorithm::GetDaughterClusterLists(), and lar_content::ClusterMopUpBaseAlgorithm::GetPfoClusterLists().

 {
     for (ClusterList::const_iterator cIter = inputClusterList.begin(), cIterEnd = inputClusterList.end(); cIter != cIterEnd; ++cIter)
     {
         const Cluster *const pCluster(*cIter);
 
         if (availabilityFlag != pCluster->IsAvailable())
             continue;
 
         const HitType hitType(LArClusterHelper::GetClusterHitType(pCluster));
 
         if ((TPC_VIEW_U != hitType) && (TPC_VIEW_V != hitType) && (TPC_VIEW_W != hitType))
             continue;
 
         ClusterList &target((TPC_VIEW_U == hitType) ? clusterListU : (TPC_VIEW_V == hitType) ? clusterListV : clusterListW);
         target.push_back(pCluster);
     }
 }

void lar_content::ClusterMopUpBaseAlgorithm::GetDaughterClusterLists	(	pandora::ClusterList &	clusterListU,
		pandora::ClusterList &	clusterListV,
		pandora::ClusterList &	clusterListW
	)		const

protectedvirtualinherited

Get the two dimensional clusters contained in the input remant/daughter cluster lists, divided into three different views.

Parameters

clusterListU	to receive the list of clusters in the u view
clusterListV	to receive the list of clusters in the v view
clusterListW	to receive the list of clusters in the w view

Definition at line 67 of file ClusterMopUpBaseAlgorithm.cc.

References lar_content::ClusterMopUpBaseAlgorithm::GetClusterLists(), and lar_content::MopUpBaseAlgorithm::m_daughterListNames.

Referenced by lar_content::ClusterMopUpBaseAlgorithm::Run().

 {
     for (const std::string &daughterListName : m_daughterListNames)
     {
         const ClusterList *pClusterList(nullptr);
 
         if (STATUS_CODE_SUCCESS != PandoraContentApi::GetList(*this, daughterListName, pClusterList))
             continue;
 
         this->GetClusterLists(*pClusterList, true, clusterListU, clusterListV, clusterListW);
     }
 }

template<typename T >

template const std::string lar_content::MopUpBaseAlgorithm::GetListName ( const T *const pT ) const

inherited

Find the name of the list hosting a specific object.

Parameters

pT	the address of the object

Returns: the name of the list

Definition at line 19 of file MopUpBaseAlgorithm.cc.

Referenced by lar_content::IsolatedClusterMopUpAlgorithm::DissolveClustersToHits(), lar_content::ClusterMopUpBaseAlgorithm::MakeClusterMerges(), and lar_content::SlidingConeClusterMopUpAlgorithm::MakeClusterMerges().

 {
     std::string currentListName;
     const MANAGED_CONTAINER<const T*> *pCurrentList(nullptr);
     (void) PandoraContentApi::GetCurrentList(*this, pCurrentList, currentListName);
 
     if (pCurrentList && (pCurrentList->end() != std::find(pCurrentList->begin(), pCurrentList->end(), pT)))
         return currentListName;
 
     for (const std::string &listName : m_daughterListNames)
     {
         const MANAGED_CONTAINER<const T*> *pList(nullptr);
         (void) PandoraContentApi::GetList(*this, listName, pList);
 
         if (pList && (pList->end() != std::find(pList->begin(), pList->end(), pT)))
             return listName;
     }
 
     throw StatusCodeException(STATUS_CODE_NOT_FOUND);
 }

void lar_content::ClusterMopUpBaseAlgorithm::GetPfoClusterLists	(	pandora::ClusterList &	clusterListU,
		pandora::ClusterList &	clusterListV,
		pandora::ClusterList &	clusterListW
	)		const

protectedvirtualinherited

Get the two dimensional clusters contained in the input pfo list, divided into three different views.

Parameters

clusterListU	to receive the list of clusters in the u view
clusterListV	to receive the list of clusters in the v view
clusterListW	to receive the list of clusters in the w view

Definition at line 45 of file ClusterMopUpBaseAlgorithm.cc.

References lar_content::ClusterMopUpBaseAlgorithm::GetClusterLists(), lar_content::LArPfoHelper::IsTrack(), lar_content::ClusterMopUpBaseAlgorithm::m_excludePfosContainingTracks, and lar_content::ClusterMopUpBaseAlgorithm::m_pfoListNames.

Referenced by lar_content::ClusterMopUpBaseAlgorithm::Run().

 {
     for (StringVector::const_iterator sIter = m_pfoListNames.begin(), sIterEnd = m_pfoListNames.end(); sIter != sIterEnd; ++sIter)
     {
         const PfoList *pPfoList = NULL;
         if (STATUS_CODE_SUCCESS != PandoraContentApi::GetList(*this, *sIter, pPfoList))
             continue;
 
         for (PfoList::const_iterator pIter = pPfoList->begin(), pIterEnd = pPfoList->end(); pIter != pIterEnd; ++pIter)
         {
             const ParticleFlowObject *const pPfo = *pIter;
 
             if (m_excludePfosContainingTracks && LArPfoHelper::IsTrack(pPfo))
                 continue;
 
             this->GetClusterLists(pPfo->GetClusterList(), false, clusterListU, clusterListV, clusterListW);
         }
     }
 }

void lar_content::ClusterMopUpBaseAlgorithm::MakeClusterMerges ( const ClusterAssociationMap & clusterAssociationMap ) const

protectedvirtualinherited

Make the cluster merges specified in the cluster association map, using list name information in the cluster list name map.

Parameters

clusterAssociationMap the cluster association map

Definition at line 104 of file ClusterMopUpBaseAlgorithm.cc.

References lar_content::MopUpBaseAlgorithm::GetListName(), max, and lar_content::LArClusterHelper::SortByNHits().

Referenced by lar_content::ConeClusterMopUpAlgorithm::ClusterMopUp(), ClusterMopUp(), and lar_content::BoundedClusterMopUpAlgorithm::ClusterMopUp().

 {
     ClusterVector sortedRemnantClusters;
     for (const auto &remnantMapEntry : clusterAssociationMap) sortedRemnantClusters.push_back(remnantMapEntry.first);
     std::sort(sortedRemnantClusters.begin(), sortedRemnantClusters.end(), LArClusterHelper::SortByNHits);
 
     for (const Cluster *const pRemnantCluster : sortedRemnantClusters)
     {
         const AssociationDetails &associationDetails(clusterAssociationMap.at(pRemnantCluster));
         const Cluster *pBestPfoCluster(nullptr);
         float bestFigureOfMerit(-std::numeric_limits<float>::max());
 
         ClusterVector sortedPfoClusters;
         for (const auto &pfoMapEntry : associationDetails) sortedPfoClusters.push_back(pfoMapEntry.first);
         std::sort(sortedPfoClusters.begin(), sortedPfoClusters.end(), LArClusterHelper::SortByNHits);
 
         for (const Cluster *const pPfoCluster : sortedPfoClusters)
         {
             const float figureOfMerit(associationDetails.at(pPfoCluster));
 
             if (figureOfMerit > bestFigureOfMerit)
             {
                 pBestPfoCluster = pPfoCluster;
                 bestFigureOfMerit = figureOfMerit;
             }
         }
 
         if (!pBestPfoCluster)
             continue;
 
         PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::MergeAndDeleteClusters(*this, pBestPfoCluster, pRemnantCluster,
             this->GetListName(pBestPfoCluster), this->GetListName(pRemnantCluster)));
     }
 }

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

privatevirtual

Reimplemented from lar_content::ClusterMopUpBaseAlgorithm.

Definition at line 91 of file NearbyClusterMopUpAlgorithm.cc.

References m_minClusterSeparation, m_minHitsInCluster, m_touchingDistance, m_vertexProximity, and lar_content::ClusterMopUpBaseAlgorithm::ReadSettings().

 {
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle,
         "MinHitsInCluster", m_minHitsInCluster));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle,
         "VertexProximity", m_vertexProximity));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle,
         "MinClusterSeparation", m_minClusterSeparation));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle,
         "TouchingDistance", m_touchingDistance));
 
     return ClusterMopUpBaseAlgorithm::ReadSettings(xmlHandle);
 }

StatusCode lar_content::ClusterMopUpBaseAlgorithm::Run ( )

protectedvirtualinherited

Definition at line 28 of file ClusterMopUpBaseAlgorithm.cc.

References lar_content::ClusterMopUpBaseAlgorithm::ClusterMopUp(), lar_content::ClusterMopUpBaseAlgorithm::GetDaughterClusterLists(), and lar_content::ClusterMopUpBaseAlgorithm::GetPfoClusterLists().

 {
     ClusterList pfoClusterListU, pfoClusterListV, pfoClusterListW;
     this->GetPfoClusterLists(pfoClusterListU, pfoClusterListV, pfoClusterListW);
 
     ClusterList daughterClusterListU, daughterClusterListV, daughterClusterListW;
     this->GetDaughterClusterLists(daughterClusterListU, daughterClusterListV, daughterClusterListW);
 
     this->ClusterMopUp(pfoClusterListU, daughterClusterListU);
     this->ClusterMopUp(pfoClusterListV, daughterClusterListV);
     this->ClusterMopUp(pfoClusterListW, daughterClusterListW);
 
     return STATUS_CODE_SUCCESS;
 }

Member Data Documentation

pandora::StringVector lar_content::MopUpBaseAlgorithm::m_daughterListNames

protectedinherited

The list of potential daughter object list names.

Definition at line 35 of file MopUpBaseAlgorithm.h.

Referenced by lar_content::SlidingConeClusterMopUpAlgorithm::GetAvailableTwoDClusters(), lar_content::ClusterMopUpBaseAlgorithm::GetDaughterClusterLists(), lar_content::SlidingConePfoMopUpAlgorithm::ReadSettings(), and lar_content::VertexBasedPfoMopUpAlgorithm::ReadSettings().

bool lar_content::ClusterMopUpBaseAlgorithm::m_excludePfosContainingTracks

protectedinherited

Whether to exclude any pfos containing clusters flagged as fixed tracks.

Definition at line 84 of file ClusterMopUpBaseAlgorithm.h.

Referenced by lar_content::ClusterMopUpBaseAlgorithm::GetPfoClusterLists(), lar_content::IsolatedClusterMopUpAlgorithm::IsolatedClusterMopUpAlgorithm(), and lar_content::ClusterMopUpBaseAlgorithm::ReadSettings().

float lar_content::NearbyClusterMopUpAlgorithm::m_minClusterSeparation

private

Minimum distance between parent and daughter clusters to declare clusters associated.

Definition at line 37 of file NearbyClusterMopUpAlgorithm.h.

Referenced by ClusterMopUp(), and ReadSettings().

unsigned int lar_content::NearbyClusterMopUpAlgorithm::m_minHitsInCluster

private

Minimum number of hits in order to consider a cluster.

Definition at line 35 of file NearbyClusterMopUpAlgorithm.h.

Referenced by ClusterMopUp(), and ReadSettings().

pandora::StringVector lar_content::ClusterMopUpBaseAlgorithm::m_pfoListNames

protectedinherited

The list of pfo list names.

Definition at line 83 of file ClusterMopUpBaseAlgorithm.h.

Referenced by lar_content::ClusterMopUpBaseAlgorithm::GetPfoClusterLists(), and lar_content::ClusterMopUpBaseAlgorithm::ReadSettings().

float lar_content::NearbyClusterMopUpAlgorithm::m_touchingDistance

private

Threshold (small) distance below which parent and daughter clusters are declated touching.

Definition at line 38 of file NearbyClusterMopUpAlgorithm.h.

Referenced by ClusterMopUp(), and ReadSettings().

float lar_content::NearbyClusterMopUpAlgorithm::m_vertexProximity

private

Distance between cluster inner/outer centroid and vtx to declare cluster vtx associated.

Definition at line 36 of file NearbyClusterMopUpAlgorithm.h.

Referenced by ClusterMopUp(), and ReadSettings().

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

/cvmfs/larsoft.opensciencegrid.org/products/larpandoracontent/v03_14_03/source/larpandoracontent/LArTwoDReco/LArClusterMopUp/NearbyClusterMopUpAlgorithm.h
/cvmfs/larsoft.opensciencegrid.org/products/larpandoracontent/v03_14_03/source/larpandoracontent/LArTwoDReco/LArClusterMopUp/NearbyClusterMopUpAlgorithm.cc

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation