SlidingConePfoMopUpAlgorithm class. More...

#include "SlidingConePfoMopUpAlgorithm.h"

Inheritance diagram for lar_content::SlidingConePfoMopUpAlgorithm:

Classes
class	ClusterMerge
	ClusterMerge class. More...

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

virtual void	MergeAndDeletePfos (const pandora::ParticleFlowObject const pPfoToEnlarge, const pandora::ParticleFlowObject const pPfoToDelete) const
	Merge and delete a pair of pfos, with a specific set of conventions for cluster merging, vertex use, etc. More...

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

Static Public Member Functions
static const pandora::Cluster *	GetParentCluster (const pandora::ClusterList &clusterList, const pandora::HitType hitType)
	Select the parent cluster (same hit type and most hits) using a provided cluster list and hit type. More...

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

Private Types
typedef std::vector< ClusterMerge >	ClusterMergeList

typedef std::unordered_map< const pandora::Cluster , const pandora::ParticleFlowObject >	ClusterToPfoMap

typedef std::unordered_map< const pandora::Cluster *, ClusterMergeList >	ClusterMergeMap

typedef std::unordered_map< const pandora::Cluster *, bool >	VertexAssociationMap

typedef std::unordered_map< const pandora::Cluster , const pandora::Cluster >	ClusterReplacementMap

Private Member Functions
pandora::StatusCode	Run ()

void	GetInteractionVertex (const pandora::Vertex *&pVertex) const
	Get the neutrino interaction vertex if it is available and if the algorithm is configured to do so. More...

void	GetThreeDClusters (pandora::ClusterVector &clusters3D, ClusterToPfoMap &clusterToPfoMap) const
	Get all 3d clusters contained in the input pfo lists and a mapping from clusters to pfos. More...

void	GetClusterMergeMap (const pandora::Vertex *const pVertex, const pandora::ClusterVector &clusters3D, const ClusterToPfoMap &clusterToPfoMap, ClusterMergeMap &clusterMergeMap) const
	Get the cluster merge map describing all potential 3d cluster merges. More...

bool	IsVertexAssociated (const pandora::Cluster const pCluster, const pandora::Vertex const pVertex, VertexAssociationMap &vertexAssociationMap, const ThreeDSlidingFitResult *const pSlidingFitResult=nullptr) const
	Whether a 3D cluster is nodally associated with a provided vertex. More...

bool	IsVertexAssociated (const pandora::Cluster const pCluster, const pandora::CartesianVector &vertexPosition, const ThreeDSlidingFitResult const pSlidingFitResult=nullptr) const
	Whether a 3D cluster is nodally associated with a provided vertex. More...

bool	MakePfoMerges (const ClusterToPfoMap &clusterToPfoMap, const ClusterMergeMap &clusterMergeMap) const
	Make pfo merges based on the provided cluster merge map. More...

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

Private Attributes
pandora::StringVector	m_inputPfoListNames
	The input pfo list names. More...

bool	m_useVertex
	Whether to use the interaction vertex to select useful cone directions. More...

bool	m_legacyMode
	Whether to use the legacy cone generation. More...

unsigned int	m_maxIterations
	The maximum allowed number of algorithm iterations. More...

unsigned int	m_maxHitsToConsider3DTrack
	The maximum number of hits in a 3d track cluster to warrant inclusion in algorithm. More...

unsigned int	m_minHitsToConsider3DShower
	The minimum number of hits in a 3d shower cluster to attempt cone fits. More...

unsigned int	m_halfWindowLayers
	The number of layers to use for half-window of sliding fit. More...

unsigned int	m_nConeFitLayers
	The number of layers over which to sum fitted direction to obtain cone fit. More...

unsigned int	m_nConeFits
	The number of cone fits to perform, spread roughly uniformly along the shower length. More...

float	m_coneLengthMultiplier
	The cone length multiplier to use when calculating bounded cluster fractions. More...

float	m_maxConeLength
	The maximum allowed cone length to use when calculating bounded cluster fractions. More...

float	m_coneTanHalfAngle1
	The cone tan half angle to use when calculating bounded cluster fractions 1. More...

float	m_coneBoundedFraction1
	The minimum cluster bounded fraction for association 1. More...

float	m_coneTanHalfAngle2
	The cone tan half angle to use when calculating bounded cluster fractions 2. More...

float	m_coneBoundedFraction2
	The minimum cluster bounded fraction for association 2. More...

float	m_minVertexLongitudinalDistance
	Vertex association check: min longitudinal distance cut. More...

float	m_maxVertexTransverseDistance
	Vertex association check: max transverse distance cut. More...

Detailed Description

SlidingConePfoMopUpAlgorithm class.

Definition at line 21 of file SlidingConePfoMopUpAlgorithm.h.

Member Typedef Documentation

typedef std::vector<ClusterMerge> lar_content::SlidingConePfoMopUpAlgorithm::ClusterMergeList

private

Definition at line 81 of file SlidingConePfoMopUpAlgorithm.h.

typedef std::unordered_map<const pandora::Cluster *, ClusterMergeList> lar_content::SlidingConePfoMopUpAlgorithm::ClusterMergeMap

private

Definition at line 102 of file SlidingConePfoMopUpAlgorithm.h.

typedef std::unordered_map<const pandora::Cluster *, const pandora::Cluster *> lar_content::SlidingConePfoMopUpAlgorithm::ClusterReplacementMap

private

Definition at line 154 of file SlidingConePfoMopUpAlgorithm.h.

typedef std::unordered_map<const pandora::Cluster *, const pandora::ParticleFlowObject *> lar_content::SlidingConePfoMopUpAlgorithm::ClusterToPfoMap

private

Definition at line 92 of file SlidingConePfoMopUpAlgorithm.h.

typedef std::unordered_map<const pandora::Cluster *, bool> lar_content::SlidingConePfoMopUpAlgorithm::VertexAssociationMap

private

Definition at line 115 of file SlidingConePfoMopUpAlgorithm.h.

Constructor & Destructor Documentation

lar_content::SlidingConePfoMopUpAlgorithm::SlidingConePfoMopUpAlgorithm ( )

Default constructor.

Definition at line 25 of file SlidingConePfoMopUpAlgorithm.cc.

                                                            :
     m_useVertex(true),
     m_legacyMode(true),
     m_maxIterations(1000),
     m_maxHitsToConsider3DTrack(100),
     m_minHitsToConsider3DShower(20),
     m_halfWindowLayers(20),
     m_nConeFitLayers(20),
     m_nConeFits(5),
     m_coneLengthMultiplier(7.f),
     m_maxConeLength(126.f),
     m_coneTanHalfAngle1(0.5f),
     m_coneBoundedFraction1(0.5f),
     m_coneTanHalfAngle2(0.75f),
     m_coneBoundedFraction2(0.75f),
     m_minVertexLongitudinalDistance(-2.5f),
     m_maxVertexTransverseDistance(3.5f)
 {
 }

Member Function Documentation

void lar_content::SlidingConePfoMopUpAlgorithm::GetClusterMergeMap	(	const pandora::Vertex *const	pVertex,
		const pandora::ClusterVector &	clusters3D,
		const ClusterToPfoMap &	clusterToPfoMap,
		ClusterMergeMap &	clusterMergeMap
	)		const

private

Get the cluster merge map describing all potential 3d cluster merges.

Parameters

pVertex	the neutrino interaction vertex, if available
clusters3D	the sorted list of 3d clusters
clusterToPfoMap	the mapping from 3d cluster to pfo
clusterMergeMap	to receive the populated cluster merge map

Definition at line 125 of file SlidingConePfoMopUpAlgorithm.cc.

References lar_content::CONE_BACKWARD_ONLY, lar_content::CONE_BOTH_DIRECTIONS, lar_content::CONE_FORWARD_ONLY, f, lar_content::SlidingConePfoMopUpAlgorithm::ClusterMerge::GetBoundedFraction1(), lar_content::SlidingConePfoMopUpAlgorithm::ClusterMerge::GetBoundedFraction2(), lar_content::LArClusterHelper::GetClusterHitType(), lar_content::ThreeDSlidingFitResult::GetGlobalMaxLayerPosition(), lar_content::ThreeDSlidingFitResult::GetGlobalMinLayerPosition(), lar_content::LArPointingCluster::GetInnerVertex(), lar_content::LArPointingCluster::GetOuterVertex(), lar_content::SlidingConePfoMopUpAlgorithm::ClusterMerge::GetParentCluster(), lar_content::LArPointingCluster::Vertex::GetPosition(), lar_content::ThreeDSlidingConeFitResult::GetSimpleConeList(), lar_content::ThreeDSlidingConeFitResult::GetSlidingFitResult(), lar_content::LArGeometryHelper::GetWirePitch(), lar_content::LArPointingClusterHelper::IsNode(), lar_content::LArPfoHelper::IsShower(), IsVertexAssociated(), m_coneBoundedFraction1, m_coneBoundedFraction2, m_coneLengthMultiplier, m_coneTanHalfAngle1, m_coneTanHalfAngle2, m_halfWindowLayers, m_legacyMode, m_maxConeLength, m_maxVertexTransverseDistance, m_minHitsToConsider3DShower, m_minVertexLongitudinalDistance, m_nConeFitLayers, and m_nConeFits.

Referenced by Run().

 {
     VertexAssociationMap vertexAssociationMap;
     const float pitchU{LArGeometryHelper::GetWirePitch(this->GetPandora(), TPC_VIEW_U)};
     const float pitchV{LArGeometryHelper::GetWirePitch(this->GetPandora(), TPC_VIEW_V)};
     const float pitchW{LArGeometryHelper::GetWirePitch(this->GetPandora(), TPC_VIEW_W)};
     const float pitchMax{std::max({pitchU, pitchV, pitchW})};
 
     for (const Cluster *const pShowerCluster : clusters3D)
     {
         if ((pShowerCluster->GetNCaloHits() < m_minHitsToConsider3DShower) || !LArPfoHelper::IsShower(clusterToPfoMap.at(pShowerCluster)))
             continue;
 
         float coneLength(0.f);
         SimpleConeList simpleConeList;
         bool isShowerVertexAssociated(false);
 
         try
         {
             float layerPitch{0.f};
             const HitType view{LArClusterHelper::GetClusterHitType(pShowerCluster)};
             if (view == TPC_VIEW_U || view == TPC_VIEW_V || view == TPC_VIEW_W)
                 layerPitch = LArGeometryHelper::GetWirePitch(this->GetPandora(), view);
             else
                 layerPitch = pitchMax;
 
             const ThreeDSlidingConeFitResult slidingConeFitResult3D(pShowerCluster, m_halfWindowLayers, layerPitch);
 
             const CartesianVector &minLayerPosition(slidingConeFitResult3D.GetSlidingFitResult().GetGlobalMinLayerPosition());
             const CartesianVector &maxLayerPosition(slidingConeFitResult3D.GetSlidingFitResult().GetGlobalMaxLayerPosition());
             coneLength = std::min(m_coneLengthMultiplier * (maxLayerPosition - minLayerPosition).GetMagnitude(), m_maxConeLength);
 
             const float vertexToMinLayer(!pVertex ? 0.f : (pVertex->GetPosition() - minLayerPosition).GetMagnitude());
             const float vertexToMaxLayer(!pVertex ? 0.f : (pVertex->GetPosition() - maxLayerPosition).GetMagnitude());
             const ConeSelection coneSelection(!pVertex      ? CONE_BOTH_DIRECTIONS
                     : (vertexToMaxLayer > vertexToMinLayer) ? CONE_FORWARD_ONLY
                                                             : CONE_BACKWARD_ONLY);
 
             slidingConeFitResult3D.GetSimpleConeList(m_nConeFitLayers, m_nConeFits, coneSelection, simpleConeList, m_coneTanHalfAngle1, m_legacyMode);
             isShowerVertexAssociated =
                 this->IsVertexAssociated(pShowerCluster, pVertex, vertexAssociationMap, &(slidingConeFitResult3D.GetSlidingFitResult()));
         }
         catch (const StatusCodeException &)
         {
             continue;
         }
 
         for (const Cluster *const pNearbyCluster : clusters3D)
         {
             if (pNearbyCluster == pShowerCluster)
                 continue;
 
             ClusterMerge bestClusterMerge(nullptr, 0.f, 0.f);
 
             for (const SimpleCone &simpleCone : simpleConeList)
             {
                 const float boundedFraction1(simpleCone.GetBoundedHitFraction(pNearbyCluster, coneLength, m_coneTanHalfAngle1));
                 const float boundedFraction2(simpleCone.GetBoundedHitFraction(pNearbyCluster, coneLength, m_coneTanHalfAngle2));
                 const ClusterMerge clusterMerge(pShowerCluster, boundedFraction1, boundedFraction2);
 
                 if (clusterMerge < bestClusterMerge)
                     bestClusterMerge = clusterMerge;
             }
 
             if (isShowerVertexAssociated && this->IsVertexAssociated(pNearbyCluster, pVertex, vertexAssociationMap))
                 continue;
 
             if (bestClusterMerge.GetParentCluster() && (bestClusterMerge.GetBoundedFraction1() > m_coneBoundedFraction1) &&
                 (bestClusterMerge.GetBoundedFraction2() > m_coneBoundedFraction2))
                 clusterMergeMap[pNearbyCluster].push_back(bestClusterMerge);
         }
     }
 
     for (ClusterMergeMap::value_type &mapEntry : clusterMergeMap)
         std::sort(mapEntry.second.begin(), mapEntry.second.end());
 }

void lar_content::SlidingConePfoMopUpAlgorithm::GetInteractionVertex ( const pandora::Vertex *& pVertex ) const

private

Get the neutrino interaction vertex if it is available and if the algorithm is configured to do so.

Parameters

pVertex to receive the neutrino interaction vertex

Definition at line 79 of file SlidingConePfoMopUpAlgorithm.cc.

References m_useVertex.

Referenced by Run().

 {
     if (!m_useVertex)
         return;
 
     const VertexList *pVertexList = nullptr;
     PANDORA_THROW_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_INITIALIZED, !=, PandoraContentApi::GetCurrentList(*this, pVertexList));
 
     pVertex =
         ((pVertexList && (pVertexList->size() == 1) && (VERTEX_3D == (*(pVertexList->begin()))->GetVertexType())) ? *(pVertexList->begin()) : nullptr);
 }

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);
 }

const Cluster * lar_content::PfoMopUpBaseAlgorithm::GetParentCluster	(	const pandora::ClusterList &	clusterList,
		const pandora::HitType	hitType
	)

staticinherited

Select the parent cluster (same hit type and most hits) using a provided cluster list and hit type.

Parameters

clusterList	the cluster list
hitType	the hit type

Returns: the address of the parent cluster

Definition at line 61 of file PfoMopUpBaseAlgorithm.cc.

Referenced by lar_content::SlidingConeClusterMopUpAlgorithm::MakeClusterMerges(), lar_content::SlidingConeClusterMopUpAlgorithm::ClusterMerge::operator<(), lar_content::SlidingConePfoMopUpAlgorithm::ClusterMerge::operator<(), and lar_content::MasterAlgorithm::StitchPfos().

 {
     unsigned int mostHits(0);
     const Cluster *pBestParentCluster(nullptr);
 
     for (const Cluster *const pParentCluster : clusterList)
     {
         if (hitType != LArClusterHelper::GetClusterHitType(pParentCluster))
             continue;
 
         const unsigned int nParentHits(pParentCluster->GetNCaloHits());
 
         if (nParentHits > mostHits)
         {
             mostHits = nParentHits;
             pBestParentCluster = pParentCluster;
         }
     }
 
     return pBestParentCluster;
 }

void lar_content::SlidingConePfoMopUpAlgorithm::GetThreeDClusters	(	pandora::ClusterVector &	clusters3D,
		ClusterToPfoMap &	clusterToPfoMap
	)		const

private

Get all 3d clusters contained in the input pfo lists and a mapping from clusters to pfos.

Parameters

clusters3D	to receive the sorted list of 3d clusters
clusterToPfoMap	to receive the mapping from 3d cluster to pfo

Definition at line 93 of file SlidingConePfoMopUpAlgorithm.cc.

References lar_content::LArPfoHelper::GetThreeDClusterList(), lar_content::LArPfoHelper::IsTrack(), m_inputPfoListNames, m_maxHitsToConsider3DTrack, and lar_content::LArClusterHelper::SortByNHits().

Referenced by Run().

 {
     for (const std::string &pfoListName : m_inputPfoListNames)
     {
         const PfoList *pPfoList(nullptr);
 
         if (STATUS_CODE_SUCCESS != PandoraContentApi::GetList(*this, pfoListName, pPfoList))
             continue;
 
         for (const Pfo *const pPfo : *pPfoList)
         {
             ClusterList pfoClusters3D;
             LArPfoHelper::GetThreeDClusterList(pPfo, pfoClusters3D);
 
             for (const Cluster *const pCluster3D : pfoClusters3D)
             {
                 if (LArPfoHelper::IsTrack(pPfo) && (pCluster3D->GetNCaloHits() > m_maxHitsToConsider3DTrack))
                     continue;
 
                 if (!clusterToPfoMap.insert(ClusterToPfoMap::value_type(pCluster3D, pPfo)).second)
                     throw StatusCodeException(STATUS_CODE_ALREADY_PRESENT);
 
                 clusters3D.push_back(pCluster3D);
             }
         }
     }
 
     std::sort(clusters3D.begin(), clusters3D.end(), LArClusterHelper::SortByNHits);
 }

bool lar_content::SlidingConePfoMopUpAlgorithm::IsVertexAssociated	(	const pandora::Cluster *const	pCluster,
		const pandora::Vertex *const	pVertex,
		VertexAssociationMap &	vertexAssociationMap,
		const ThreeDSlidingFitResult *const	pSlidingFitResult = `nullptr`
	)		const

private

Whether a 3D cluster is nodally associated with a provided vertex.

Parameters

pCluster	the address of the cluster
pVertex	the address of the vertex
vertexAssociationMap	the vertex association map, which caches the results for efficiency
pSlidingFitResult	the address of a 3D sliding linear fit, if available, which aids pointing cluster construction efficiency

Returns: boolean

Referenced by GetClusterMergeMap().

bool lar_content::SlidingConePfoMopUpAlgorithm::IsVertexAssociated	(	const pandora::Cluster *const	pCluster,
		const pandora::CartesianVector &	vertexPosition,
		const ThreeDSlidingFitResult *const	pSlidingFitResult = `nullptr`
	)		const

private

Whether a 3D cluster is nodally associated with a provided vertex.

Parameters

pCluster	the address of the cluster
vertexPosition	the vertex position
pSlidingFitResult	the address of a 3D sliding linear fit, if available, which aids pointing cluster construction efficiency

Returns: boolean

bool lar_content::SlidingConePfoMopUpAlgorithm::MakePfoMerges	(	const ClusterToPfoMap &	clusterToPfoMap,
		const ClusterMergeMap &	clusterMergeMap
	)		const

private

Make pfo merges based on the provided cluster merge map.

Parameters

clusterToPfoMap	the mapping from 3d cluster to pfo
clusterMergeMap	the populated cluster merge map

Returns: whether a pfo merge has been made

Definition at line 259 of file SlidingConePfoMopUpAlgorithm.cc.

References lar_content::PfoMopUpBaseAlgorithm::MergeAndDeletePfos(), and lar_content::LArClusterHelper::SortByNHits().

Referenced by Run().

 {
     ClusterVector daughterClusters;
     for (const ClusterMergeMap::value_type &mapEntry : clusterMergeMap)
         daughterClusters.push_back(mapEntry.first);
     std::sort(daughterClusters.begin(), daughterClusters.end(), LArClusterHelper::SortByNHits);
 
     bool pfosMerged(false);
     ClusterReplacementMap clusterReplacementMap;
 
     for (ClusterVector::const_reverse_iterator rIter = daughterClusters.rbegin(), rIterEnd = daughterClusters.rend(); rIter != rIterEnd; ++rIter)
     {
         const Cluster *const pDaughterCluster(*rIter);
 
         if (clusterReplacementMap.count(pDaughterCluster))
             throw StatusCodeException(STATUS_CODE_FAILURE);
 
         const Cluster *pParentCluster(clusterMergeMap.at(pDaughterCluster).at(0).GetParentCluster());
 
         if (clusterReplacementMap.count(pParentCluster))
             pParentCluster = clusterReplacementMap.at(pParentCluster);
 
         // ATTN Sign there was a reciprocal relationship in the cluster merge map (already actioned once)
         if (pDaughterCluster == pParentCluster)
             continue;
 
         // Key book-keeping on clusters and use cluster->pfo lookup
         const Pfo *const pDaughterPfo(clusterToPfoMap.at(pDaughterCluster));
         const Pfo *const pParentPfo(clusterToPfoMap.at(pParentCluster));
         this->MergeAndDeletePfos(pParentPfo, pDaughterPfo);
         pfosMerged = true;
 
         // Simple/placeholder book-keeping for reciprocal relationships and progressive merges
         clusterReplacementMap[pDaughterCluster] = pParentCluster;
 
         for (ClusterReplacementMap::value_type &mapEntry : clusterReplacementMap)
         {
             if (pDaughterCluster == mapEntry.second)
                 mapEntry.second = pParentCluster;
         }
     }
 
     return pfosMerged;
 }

void lar_content::PfoMopUpBaseAlgorithm::MergeAndDeletePfos	(	const pandora::ParticleFlowObject *const	pPfoToEnlarge,
		const pandora::ParticleFlowObject *const	pPfoToDelete
	)		const

virtualinherited

Merge and delete a pair of pfos, with a specific set of conventions for cluster merging, vertex use, etc.

Parameters

pPfoToEnlarge	the address of the pfo to enlarge
pPfoToDelete	the address of the pfo to delete (will become a dangling pointer)

Definition at line 20 of file PfoMopUpBaseAlgorithm.cc.

Referenced by MakePfoMerges(), and lar_content::VertexBasedPfoMopUpAlgorithm::MergePfos().

 {
     if (pPfoToEnlarge == pPfoToDelete)
         throw StatusCodeException(STATUS_CODE_NOT_ALLOWED);
 
     const PfoList daughterPfos(pPfoToDelete->GetDaughterPfoList());
     const ClusterVector daughterClusters(pPfoToDelete->GetClusterList().begin(), pPfoToDelete->GetClusterList().end());
     const VertexVector daughterVertices(pPfoToDelete->GetVertexList().begin(), pPfoToDelete->GetVertexList().end());
 
     PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::Delete(*this, pPfoToDelete, this->GetListName(pPfoToDelete)));
 
     for (const ParticleFlowObject *const pDaughterPfo : daughterPfos)
     {
         PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::SetPfoParentDaughterRelationship(*this, pPfoToEnlarge, pDaughterPfo));
     }
 
     for (const Vertex *const pDaughterVertex : daughterVertices)
     {
         PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::Delete(*this, pDaughterVertex, this->GetListName(pDaughterVertex)));
     }
 
     for (const Cluster *const pDaughterCluster : daughterClusters)
     {
         const HitType daughterHitType(LArClusterHelper::GetClusterHitType(pDaughterCluster));
         const Cluster *pParentCluster(PfoMopUpBaseAlgorithm::GetParentCluster(pPfoToEnlarge->GetClusterList(), daughterHitType));
 
         if (pParentCluster)
         {
             PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=,
                 PandoraContentApi::MergeAndDeleteClusters(
                     *this, pParentCluster, pDaughterCluster, this->GetListName(pParentCluster), this->GetListName(pDaughterCluster)));
         }
         else
         {
             PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::AddToPfo(*this, pPfoToEnlarge, pDaughterCluster));
         }
     }
 }

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

privatevirtual

Reimplemented from lar_content::PfoMopUpBaseAlgorithm.

Definition at line 330 of file SlidingConePfoMopUpAlgorithm.cc.

References m_coneBoundedFraction1, m_coneBoundedFraction2, m_coneLengthMultiplier, m_coneTanHalfAngle1, m_coneTanHalfAngle2, lar_content::MopUpBaseAlgorithm::m_daughterListNames, m_halfWindowLayers, m_inputPfoListNames, m_legacyMode, m_maxConeLength, m_maxHitsToConsider3DTrack, m_maxIterations, m_maxVertexTransverseDistance, m_minHitsToConsider3DShower, m_minVertexLongitudinalDistance, m_nConeFitLayers, m_nConeFits, m_useVertex, and lar_content::PfoMopUpBaseAlgorithm::ReadSettings().

 {
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadVectorOfValues(xmlHandle, "InputPfoListNames", m_inputPfoListNames));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "UseVertex", m_useVertex));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "LegacyMode", m_legacyMode));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MaxIterations", m_maxIterations));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
         XmlHelper::ReadValue(xmlHandle, "MaxHitsToConsider3DTrack", m_maxHitsToConsider3DTrack));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
         XmlHelper::ReadValue(xmlHandle, "MinHitsToConsider3DShower", m_minHitsToConsider3DShower));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "SlidingFitHalfWindow", m_halfWindowLayers));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "NConeFitLayers", m_nConeFitLayers));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "NConeFits", m_nConeFits));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "ConeLengthMultiplier", m_coneLengthMultiplier));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MaxConeLength", m_maxConeLength));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "ConeTanHalfAngle1", m_coneTanHalfAngle1));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "ConeBoundedFraction1", m_coneBoundedFraction1));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "ConeTanHalfAngle2", m_coneTanHalfAngle2));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "ConeBoundedFraction2", m_coneBoundedFraction2));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
         XmlHelper::ReadValue(xmlHandle, "MinVertexLongitudinalDistance", m_minVertexLongitudinalDistance));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
         XmlHelper::ReadValue(xmlHandle, "MaxVertexTransverseDistance", m_maxVertexTransverseDistance));
 
     m_daughterListNames.insert(m_daughterListNames.end(), m_inputPfoListNames.begin(), m_inputPfoListNames.end());
 
     return PfoMopUpBaseAlgorithm::ReadSettings(xmlHandle);
 }

StatusCode lar_content::SlidingConePfoMopUpAlgorithm::Run ( )

private

Definition at line 47 of file SlidingConePfoMopUpAlgorithm.cc.

References GetClusterMergeMap(), GetInteractionVertex(), GetThreeDClusters(), m_maxIterations, m_useVertex, and MakePfoMerges().

 {
     const Vertex *pVertex(nullptr);
     this->GetInteractionVertex(pVertex);
 
     if (m_useVertex && !pVertex)
     {
         if (PandoraContentApi::GetSettings(*this)->ShouldDisplayAlgorithmInfo())
             std::cout << "SlidingConePfoMopUpAlgorithm - interaction vertex not available for use." << std::endl;
         return STATUS_CODE_SUCCESS;
     }
 
     unsigned int nIterations(0);
 
     while (nIterations++ < m_maxIterations)
     {
         ClusterVector clusters3D;
         ClusterToPfoMap clusterToPfoMap;
         this->GetThreeDClusters(clusters3D, clusterToPfoMap);
 
         ClusterMergeMap clusterMergeMap;
         this->GetClusterMergeMap(pVertex, clusters3D, clusterToPfoMap, clusterMergeMap);
 
         if (!this->MakePfoMerges(clusterToPfoMap, clusterMergeMap))
             break;
     }
 
     return STATUS_CODE_SUCCESS;
 }

Member Data Documentation

float lar_content::SlidingConePfoMopUpAlgorithm::m_coneBoundedFraction1

private

The minimum cluster bounded fraction for association 1.

Definition at line 168 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by GetClusterMergeMap(), and ReadSettings().

float lar_content::SlidingConePfoMopUpAlgorithm::m_coneBoundedFraction2

private

The minimum cluster bounded fraction for association 2.

Definition at line 170 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by GetClusterMergeMap(), and ReadSettings().

float lar_content::SlidingConePfoMopUpAlgorithm::m_coneLengthMultiplier

private

The cone length multiplier to use when calculating bounded cluster fractions.

Definition at line 165 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by GetClusterMergeMap(), and ReadSettings().

float lar_content::SlidingConePfoMopUpAlgorithm::m_coneTanHalfAngle1

private

The cone tan half angle to use when calculating bounded cluster fractions 1.

Definition at line 167 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by GetClusterMergeMap(), and ReadSettings().

float lar_content::SlidingConePfoMopUpAlgorithm::m_coneTanHalfAngle2

private

The cone tan half angle to use when calculating bounded cluster fractions 2.

Definition at line 169 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by GetClusterMergeMap(), and ReadSettings().

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(), ReadSettings(), and lar_content::VertexBasedPfoMopUpAlgorithm::ReadSettings().

unsigned int lar_content::SlidingConePfoMopUpAlgorithm::m_halfWindowLayers

private

The number of layers to use for half-window of sliding fit.

Definition at line 162 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by GetClusterMergeMap(), and ReadSettings().

pandora::StringVector lar_content::SlidingConePfoMopUpAlgorithm::m_inputPfoListNames

private

The input pfo list names.

Definition at line 156 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by GetThreeDClusters(), and ReadSettings().

bool lar_content::SlidingConePfoMopUpAlgorithm::m_legacyMode

private

Whether to use the legacy cone generation.

Definition at line 158 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by GetClusterMergeMap(), and ReadSettings().

float lar_content::SlidingConePfoMopUpAlgorithm::m_maxConeLength

private

The maximum allowed cone length to use when calculating bounded cluster fractions.

Definition at line 166 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by GetClusterMergeMap(), and ReadSettings().

unsigned int lar_content::SlidingConePfoMopUpAlgorithm::m_maxHitsToConsider3DTrack

private

The maximum number of hits in a 3d track cluster to warrant inclusion in algorithm.

Definition at line 160 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by GetThreeDClusters(), and ReadSettings().

unsigned int lar_content::SlidingConePfoMopUpAlgorithm::m_maxIterations

private

The maximum allowed number of algorithm iterations.

Definition at line 159 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by ReadSettings(), and Run().

float lar_content::SlidingConePfoMopUpAlgorithm::m_maxVertexTransverseDistance

private

Vertex association check: max transverse distance cut.

Definition at line 172 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by GetClusterMergeMap(), and ReadSettings().

unsigned int lar_content::SlidingConePfoMopUpAlgorithm::m_minHitsToConsider3DShower

private

The minimum number of hits in a 3d shower cluster to attempt cone fits.

Definition at line 161 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by GetClusterMergeMap(), and ReadSettings().

float lar_content::SlidingConePfoMopUpAlgorithm::m_minVertexLongitudinalDistance

private

Vertex association check: min longitudinal distance cut.

Definition at line 171 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by GetClusterMergeMap(), and ReadSettings().

unsigned int lar_content::SlidingConePfoMopUpAlgorithm::m_nConeFitLayers

private

The number of layers over which to sum fitted direction to obtain cone fit.

Definition at line 163 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by GetClusterMergeMap(), and ReadSettings().

unsigned int lar_content::SlidingConePfoMopUpAlgorithm::m_nConeFits

private

The number of cone fits to perform, spread roughly uniformly along the shower length.

Definition at line 164 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by GetClusterMergeMap(), and ReadSettings().

bool lar_content::SlidingConePfoMopUpAlgorithm::m_useVertex

private

Whether to use the interaction vertex to select useful cone directions.

Definition at line 157 of file SlidingConePfoMopUpAlgorithm.h.

Referenced by GetInteractionVertex(), ReadSettings(), and Run().

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

larpandoracontent/v04_10_00/source/larpandoracontent/LArThreeDReco/LArPfoMopUp/SlidingConePfoMopUpAlgorithm.h
larpandoracontent/v04_10_00/source/larpandoracontent/LArThreeDReco/LArPfoMopUp/SlidingConePfoMopUpAlgorithm.cc

Classes

Public Member Functions

Static Public Member Functions

Protected Attributes

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation