TwoViewCosmicRayRemovalTool class. More...

#include <TwoViewCosmicRayRemovalTool.h>

Inheritance diagram for lar_content::TwoViewCosmicRayRemovalTool:

Public Types
typedef std::vector < pandora::HitType >	HitTypeVector

Public Types inherited from lar_content::DeltaRayMatrixTool
typedef TwoViewDeltaRayMatchingAlgorithm::MatchingType::MatrixType	MatrixType

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

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

Private Member Functions
bool	Run (TwoViewDeltaRayMatchingAlgorithm *const pAlgorithm, MatrixType &overlapMatrix)
	Run the algorithm tool. More...

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

bool	RemoveCosmicRayHits (const MatrixType::ElementList &elementList) const
	Remove hits from delta ray clusters that belong to the parent muon. More...

bool	PassElementChecks (const MatrixType::Element &element, const pandora::HitType hitType) const
	Determine whether element satifies simple checks. More...

bool	IsMuonEndpoint (const MatrixType::Element &element, const pandora::HitType hitType) const
	Determine whether the matched clusters suggest that the delta ray is at the endpoint of the cosmic ray (and is likely to be a michel) More...

bool	IsContaminated (const MatrixType::Element &element, const pandora::HitType hitType) const
	Determine whether the cluster under investigation has muon contamination. More...

bool	IsCloseToLine (const pandora::CartesianVector &hitPosition, const pandora::CartesianVector &lineStart, const pandora::CartesianVector &lineEnd, const float distanceToLine) const
	Whether a given position is close to a defined line. More...

bool	IsInLineSegment (const pandora::CartesianVector &lowerBoundary, const pandora::CartesianVector &upperBoundary, const pandora::CartesianVector &point) const
	Whether the projection of a given position lies on a defined line. More...

bool	IsBestElement (const MatrixType::Element &element, const pandora::HitType hitType, const MatrixType::ElementList &elementList) const
	Determine whether the input element is the best to use to modify the contaminated cluster (best is defined by the total hit count) More...

void	CreateSeed (const MatrixType::Element &element, const pandora::HitType hitType, pandora::CaloHitList &collectedHits) const
	Collect hits in the delta ray cluster that lie close to calculated projected positions forming a seed to later grow. More...

pandora::StatusCode	GrowSeed (const MatrixType::Element &element, const pandora::HitType hitType, pandora::CaloHitList &collectedHits, pandora::CaloHitList &deltaRayRemantHits) const
	Examine remaining hits in the delta ray cluster adding them to the delta ray seed or parent muon if appropriate and a separate list if not. More...

void	CollectHitsFromDeltaRay (const pandora::CartesianVector &positionOnMuon, const pandora::CartesianVector &muonDirection, const pandora::Cluster *const pDeltaRayCluster, const float &minDistanceFromMuon, const bool demandCloserToCollected, const pandora::CaloHitList &protectedHits, pandora::CaloHitList &collectedHits) const
	Collect hits from the delta ray cluster to either keep (demandCloserToCollected == true) or separate into a new shower (demandCloserToCollected == false) More...

void	SplitDeltaRayCluster (const MatrixType::Element &element, const pandora::HitType hitType, pandora::CaloHitList &collectedHits, pandora::CaloHitList &deltaRayRemnantHits) const
	Fragment the delta ray cluster, refining the matched delta ray cluster perhaps creating significant remnants in the process. More...

void	ReclusterRemnant (const pandora::HitType hitType, const pandora::Cluster const pMuonCluster, const pandora::Cluster const pDeltaRayRemnant, pandora::ClusterVector &clusterVector, pandora::PfoVector &pfoVector) const
	Reculster the remnant cluster, merging insignificant created clusters into the parent muon (if proximity checks pass) More...

pandora::StatusCode	ProjectDeltaRayPositions (const MatrixType::Element &element, const pandora::HitType hitType, pandora::CartesianPointVector &projectedPositions) const
	Use two views of a delta ray pfo to calculate projected positions in a given third view. More...

Private Attributes
float	m_minSeparation
	The minimum delta ray - parent muon cluster separation required to investigate delta ray cluster. More...

unsigned int	m_slidingFitWindow
	The sliding fit window used in cosmic ray parameterisations. More...

float	m_distanceToLine
	The maximum perpendicular distance of a position to a line for it to be considered close. More...

float	m_minContaminationLength
	The minimum projected length of a delta ray cluster onto the muon track for it to be considered contaminated. More...

float	m_maxDistanceToHit
	The maximum distance of a hit from the cosmic ray track for it to be added into the CR. More...

unsigned int	m_minRemnantClusterSize
	The minimum hit number of a remnant cluster for it to be considered significant. More...

float	m_maxDistanceToTrack
	The minimum distance of an insignificant remnant cluster from the cosmic ray track for it to be merged into the CR. More...

Additional Inherited Members
Public Attributes inherited from lar_content::DeltaRayMatrixTool
TwoViewDeltaRayMatchingAlgorithm *	m_pParentAlgorithm
	Address of the parent matching algorithm. More...

Detailed Description

TwoViewCosmicRayRemovalTool class.

Definition at line 18 of file TwoViewCosmicRayRemovalTool.h.

Member Typedef Documentation

typedef std::vector<pandora::HitType> lar_content::TwoViewCosmicRayRemovalTool::HitTypeVector

Definition at line 21 of file TwoViewCosmicRayRemovalTool.h.

Constructor & Destructor Documentation

lar_content::TwoViewCosmicRayRemovalTool::TwoViewCosmicRayRemovalTool ( )

Default constructor.

Definition at line 22 of file TwoViewCosmicRayRemovalTool.cc.

                                                          :
     m_minSeparation(2.f),
     m_slidingFitWindow(10000),
     m_distanceToLine(0.5f),
     m_minContaminationLength(3.f),
     m_maxDistanceToHit(1.f),
     m_minRemnantClusterSize(3),
     m_maxDistanceToTrack(2.f)
 {
 }

Member Function Documentation

void lar_content::TwoViewCosmicRayRemovalTool::CollectHitsFromDeltaRay	(	const pandora::CartesianVector &	positionOnMuon,
		const pandora::CartesianVector &	muonDirection,
		const pandora::Cluster *const	pDeltaRayCluster,
		const float &	minDistanceFromMuon,
		const bool	demandCloserToCollected,
		const pandora::CaloHitList &	protectedHits,
		pandora::CaloHitList &	collectedHits
	)		const

private

Collect hits from the delta ray cluster to either keep (demandCloserToCollected == true) or separate into a new shower (demandCloserToCollected == false)

Parameters

positionOnMuon	the parameterised muon position
muonDirection	the parameterised muon direction
pDeltaRayCluster	the delta ray cluster under investigation
minDistanceFromMuon	the minimum distance of a hit from the muon track for it to not belong to the muon
demandCloserToCollected	whether to demand a hit be closer to the collected hits than to the muon hits for it to be collected
protectedHits	the hits that are protected from being collected
collectedHits	the output list of collected hits

Definition at line 374 of file TwoViewCosmicRayRemovalTool.cc.

 {
     CaloHitList deltaRayHitList;
     pDeltaRayCluster->GetOrderedCaloHitList().FillCaloHitList(deltaRayHitList);
 
     bool hitsAdded(false);
 
     do
     {
         hitsAdded = false;
 
         for (const CaloHit *const pCaloHit : deltaRayHitList)
         {
             if (std::find(protectedHits.begin(), protectedHits.end(), pCaloHit) != protectedHits.end())
                 continue;
 
             if (std::find(collectedHits.begin(), collectedHits.end(), pCaloHit) != collectedHits.end())
                 continue;
 
             const float distanceToCollectedHits(collectedHits.empty()
                                                     ? std::numeric_limits<float>::max()
                                                     : LArClusterHelper::GetClosestDistance(pCaloHit->GetPositionVector(), collectedHits));
             const float distanceToMuonHits(muonDirection.GetCrossProduct(pCaloHit->GetPositionVector() - positionOnMuon).GetMagnitude());
 
             if ((distanceToMuonHits < minDistanceFromMuon) || (demandCloserToCollected && (distanceToCollectedHits > distanceToMuonHits)))
                 continue;
 
             collectedHits.push_back(pCaloHit);
             hitsAdded = true;
         }
     } while (hitsAdded);
 }

void lar_content::TwoViewCosmicRayRemovalTool::CreateSeed	(	const MatrixType::Element &	element,
		const pandora::HitType	hitType,
		pandora::CaloHitList &	collectedHits
	)		const

private

Collect hits in the delta ray cluster that lie close to calculated projected positions forming a seed to later grow.

Parameters

element	the matrix element
hitType	the hit type of the cluster under investigation
collectedHits	the output list of identified delta ray hits

Definition at line 302 of file TwoViewCosmicRayRemovalTool.cc.

 {
     // To avoid fluctuations, parameterise the muon track
     CartesianVector positionOnMuon(0.f, 0.f, 0.f), muonDirection(0.f, 0.f, 0.f);
     if (m_pParentAlgorithm->ParameteriseMuon(element.GetOverlapResult().GetCommonMuonPfoList().front(),
             m_pParentAlgorithm->GetCluster(element, hitType), positionOnMuon, muonDirection) != STATUS_CODE_SUCCESS)
         return;
 
     CartesianPointVector deltaRayProjectedPositions;
     if (this->ProjectDeltaRayPositions(element, hitType, deltaRayProjectedPositions) != STATUS_CODE_SUCCESS)
         return;
 
     CaloHitList deltaRayHitList;
     m_pParentAlgorithm->GetCluster(element, hitType)->GetOrderedCaloHitList().FillCaloHitList(deltaRayHitList);
 
     CaloHitVector deltaRayHitVector(deltaRayHitList.begin(), deltaRayHitList.end());
     std::sort(deltaRayHitVector.begin(), deltaRayHitVector.end(), LArClusterHelper::SortHitsByPosition);
 
     for (const CaloHit *const pCaloHit : deltaRayHitVector)
     {
         const CartesianVector &position(pCaloHit->GetPositionVector());
 
         for (const CartesianVector &projectedPosition : deltaRayProjectedPositions)
         {
             const float distanceToProjectionSquared((position - projectedPosition).GetMagnitudeSquared());
 
             if (distanceToProjectionSquared < m_maxDistanceToHit * m_maxDistanceToHit)
             {
                 // To prevent gappy seeds
                 if ((!collectedHits.empty()) && (LArClusterHelper::GetClosestDistance(position, collectedHits) > m_maxDistanceToHit))
                     continue;
 
                 const float distanceToMuonHitsSquared(muonDirection.GetCrossProduct(position - positionOnMuon).GetMagnitudeSquared());
 
                 if (distanceToMuonHitsSquared < m_maxDistanceToHit * m_maxDistanceToHit)
                     continue;
 
                 collectedHits.push_back(pCaloHit);
 
                 break;
             }
         }
     }
 }

StatusCode lar_content::TwoViewCosmicRayRemovalTool::GrowSeed	(	const MatrixType::Element &	element,
		const pandora::HitType	hitType,
		pandora::CaloHitList &	collectedHits,
		pandora::CaloHitList &	deltaRayRemantHits
	)		const

private

Examine remaining hits in the delta ray cluster adding them to the delta ray seed or parent muon if appropriate and a separate list if not.

Parameters

element	the matrix element
hitType	the hit type of the cluster under investigation
collectedHits	the list of identified delta ray hits
deltaRayRemantHits	the list of remainder hits

Returns: whether the muon projection mechanics were successful - abort process if not

Definition at line 349 of file TwoViewCosmicRayRemovalTool.cc.

 {
     const Cluster *const pDeltaRayCluster(m_pParentAlgorithm->GetCluster(element, hitType)), *pMuonCluster(nullptr);
 
     if (m_pParentAlgorithm->GetMuonCluster(element.GetOverlapResult().GetCommonMuonPfoList(), hitType, pMuonCluster) != STATUS_CODE_SUCCESS)
         return STATUS_CODE_NOT_FOUND;
 
     // To avoid fluctuatuions, parameterise the muon track
     CartesianVector positionOnMuon(0.f, 0.f, 0.f), muonDirection(0.f, 0.f, 0.f);
     if (m_pParentAlgorithm->ParameteriseMuon(element.GetOverlapResult().GetCommonMuonPfoList().front(), pDeltaRayCluster, positionOnMuon,
             muonDirection) != STATUS_CODE_SUCCESS)
         return STATUS_CODE_NOT_FOUND;
 
     // Identify delta ray hits
     this->CollectHitsFromDeltaRay(positionOnMuon, muonDirection, pDeltaRayCluster, m_maxDistanceToHit, true, deltaRayHits, deltaRayHits);
 
     // Identify remnant hits
     this->CollectHitsFromDeltaRay(positionOnMuon, muonDirection, pDeltaRayCluster, m_maxDistanceToHit, false, deltaRayHits, remnantHits);
 
     return STATUS_CODE_SUCCESS;
 }

bool lar_content::TwoViewCosmicRayRemovalTool::IsBestElement	(	const MatrixType::Element &	element,
		const pandora::HitType	hitType,
		const MatrixType::ElementList &	elementList
	)		const

private

Determine whether the input element is the best to use to modify the contaminated cluster (best is defined by the total hit count)

Parameters

element	the matrix element
hitType	the hit type of the cluster under investigation
elementList	the matrix element list

Returns: whether the input element is the best element

Definition at line 274 of file TwoViewCosmicRayRemovalTool.cc.

 {
     const float chiSquared(element.GetOverlapResult().GetReducedChiSquared());
     const unsigned int hitSum(element.GetCluster1()->GetNCaloHits() + element.GetCluster2()->GetNCaloHits());
 
     for (const MatrixType::Element &testElement : elementList)
     {
         if (m_pParentAlgorithm->GetCluster(testElement, hitType) != m_pParentAlgorithm->GetCluster(element, hitType))
             continue;
 
         if ((testElement.GetCluster1() == element.GetCluster1()) && (testElement.GetCluster2() == element.GetCluster2()))
             continue;
 
         const float testChiSquared(testElement.GetOverlapResult().GetReducedChiSquared());
         const unsigned int testHitSum(testElement.GetCluster1()->GetNCaloHits() + testElement.GetCluster2()->GetNCaloHits());
 
         if ((testHitSum < hitSum) || ((testHitSum == hitSum) && (testChiSquared > chiSquared)))
             continue;
 
         if (this->PassElementChecks(testElement, hitType))
             return false;
     }
 
     return true;
 }

bool lar_content::TwoViewCosmicRayRemovalTool::IsCloseToLine	(	const pandora::CartesianVector &	hitPosition,
		const pandora::CartesianVector &	lineStart,
		const pandora::CartesianVector &	lineEnd,
		const float	distanceToLine
	)		const

private

Whether a given position is close to a defined line.

Parameters

hitPosition	the input position
lineStart	the start position of the line
lineEnd	the end position of the line
distanceToLine	the definition of close

Returns: whether the position is close to the definied line

Definition at line 234 of file TwoViewCosmicRayRemovalTool.cc.

 {
     CartesianVector lineDirection(lineStart - lineEnd);
     lineDirection = lineDirection.GetUnitVector();
 
     const float transverseDistanceFromLine(lineDirection.GetCrossProduct(hitPosition - lineStart).GetMagnitude());
 
     if (transverseDistanceFromLine > distanceToLine)
         return false;
 
     return true;
 }

bool lar_content::TwoViewCosmicRayRemovalTool::IsContaminated	(	const MatrixType::Element &	element,
		const pandora::HitType	hitType
	)		const

private

Determine whether the cluster under investigation has muon contamination.

Parameters

element	the matrix element
hitType	the hit type of the cluster under investigation

Returns: whether the cluster contains muon hits to remove

Definition at line 177 of file TwoViewCosmicRayRemovalTool.cc.

 {
     const Cluster *pMuonCluster(nullptr), *const pDeltaRayCluster(m_pParentAlgorithm->GetCluster(element, hitType));
 
     if (m_pParentAlgorithm->GetMuonCluster(element.GetOverlapResult().GetCommonMuonPfoList(), hitType, pMuonCluster) != STATUS_CODE_SUCCESS)
         return false;
 
     const float slidingFitPitch(LArGeometryHelper::GetWireZPitch(this->GetPandora()));
     const TwoDSlidingFitResult slidingFitResult(pMuonCluster, m_slidingFitWindow, slidingFitPitch);
 
     CartesianVector muonDirection(0.f, 0.f, 0.f);
     slidingFitResult.GetGlobalDirection(slidingFitResult.GetLayerFitResultMap().begin()->second.GetGradient(), muonDirection);
 
     CartesianVector deltaRayVertex(0.f, 0.f, 0.f), muonVertex(0.f, 0.f, 0.f);
     LArClusterHelper::GetClosestPositions(pDeltaRayCluster, pMuonCluster, deltaRayVertex, muonVertex);
 
     // Find furthest point in DR cluster that lies on the projected muon track
     float furthestSeparation(0.f);
     CartesianVector extendedPoint(0.f, 0.f, 0.f);
 
     CaloHitList deltaRayHitList;
     pDeltaRayCluster->GetOrderedCaloHitList().FillCaloHitList(deltaRayHitList);
 
     for (const CaloHit *const pCaloHit : deltaRayHitList)
     {
         const CartesianVector &position(pCaloHit->GetPositionVector());
         const float separation((position - muonVertex).GetMagnitude());
 
         if (separation > furthestSeparation)
         {
             if (this->IsCloseToLine(position, muonVertex, muonVertex + muonDirection, m_distanceToLine))
             {
                 furthestSeparation = separation;
                 extendedPoint = position;
             }
         }
     }
 
     // Check if delta ray has significant track contamination
     if (furthestSeparation < m_minContaminationLength)
         return false;
 
     // ATTN: Avoid cases where opening angle is small - it is easy to make mistakes in these instances
     CaloHitList muonHitList;
     pMuonCluster->GetOrderedCaloHitList().FillCaloHitList(muonHitList);
 
     for (const CaloHit *const pCaloHit : muonHitList)
     {
         if (this->IsInLineSegment(deltaRayVertex, extendedPoint, pCaloHit->GetPositionVector()))
             return false;
     }
 
     return true;
 }

bool lar_content::TwoViewCosmicRayRemovalTool::IsInLineSegment	(	const pandora::CartesianVector &	lowerBoundary,
		const pandora::CartesianVector &	upperBoundary,
		const pandora::CartesianVector &	point
	)		const

private

Whether the projection of a given position lies on a defined line.

Parameters

lowerBoundary	the start position of the line
upperBoundary	the end position of the line
point	the input position

Returns: whether the position lies between the two points

Definition at line 250 of file TwoViewCosmicRayRemovalTool.cc.

 {
     const float segmentBoundaryGradient = (-1.f) * (upperBoundary.GetX() - lowerBoundary.GetX()) / (upperBoundary.GetZ() - lowerBoundary.GetZ());
     const float xPointOnUpperLine((point.GetZ() - upperBoundary.GetZ()) / segmentBoundaryGradient + upperBoundary.GetX());
     const float xPointOnLowerLine((point.GetZ() - lowerBoundary.GetZ()) / segmentBoundaryGradient + lowerBoundary.GetX());
 
     if (std::fabs(xPointOnUpperLine - point.GetX()) < std::numeric_limits<float>::epsilon())
         return true;
 
     if (std::fabs(xPointOnLowerLine - point.GetX()) < std::numeric_limits<float>::epsilon())
         return true;
 
     if ((point.GetX() > xPointOnUpperLine) && (point.GetX() > xPointOnLowerLine))
         return false;
 
     if ((point.GetX() < xPointOnUpperLine) && (point.GetX() < xPointOnLowerLine))
         return false;
 
     return true;
 }

bool lar_content::TwoViewCosmicRayRemovalTool::IsMuonEndpoint	(	const MatrixType::Element &	element,
		const pandora::HitType	hitType
	)		const

private

Determine whether the matched clusters suggest that the delta ray is at the endpoint of the cosmic ray (and is likely to be a michel)

Parameters

element	the matrix element
ignoreHitType	whether to ignore the cluster under investigation
hitTypeToIgnore	the hit type to ignore

Returns: whether the delta ray is at the endpoint of the cosmic ray

Definition at line 141 of file TwoViewCosmicRayRemovalTool.cc.

 {
     for (const HitType &otherHitType : m_pParentAlgorithm->GetHitTypeVector())
     {
         if (otherHitType == hitType)
             continue;
 
         const Cluster *pMuonCluster(nullptr), *const pDeltaRayCluster(m_pParentAlgorithm->GetCluster(element, hitType));
 
         if (m_pParentAlgorithm->GetMuonCluster(element.GetOverlapResult().GetCommonMuonPfoList(), hitType, pMuonCluster) != STATUS_CODE_SUCCESS)
             return false;
 
         float xMinDR(+std::numeric_limits<float>::max()), xMaxDR(-std::numeric_limits<float>::max());
         pDeltaRayCluster->GetClusterSpanX(xMinDR, xMaxDR);
 
         float xMinCR(-std::numeric_limits<float>::max()), xMaxCR(+std::numeric_limits<float>::max());
         pMuonCluster->GetClusterSpanX(xMinCR, xMaxCR);
 
         if ((xMinDR < xMinCR) || (xMaxDR > xMaxCR))
             return false;
 
         float zMinDR(+std::numeric_limits<float>::max()), zMaxDR(-std::numeric_limits<float>::max());
         pDeltaRayCluster->GetClusterSpanZ(xMinDR, xMaxDR, zMinDR, zMaxDR);
 
         float zMinCR(-std::numeric_limits<float>::max()), zMaxCR(+std::numeric_limits<float>::max());
         pMuonCluster->GetClusterSpanZ(xMinCR, xMaxCR, zMinCR, zMaxCR);
 
         if ((zMinDR < zMinCR) || (zMaxDR > zMaxCR))
             return false;
     }
 
     return true;
 }

bool lar_content::TwoViewCosmicRayRemovalTool::PassElementChecks	(	const MatrixType::Element &	element,
		const pandora::HitType	hitType
	)		const

private

Determine whether element satifies simple checks.

Parameters

element	the matrix element
hitType	the hit type of the cluster under investigation

Returns: whether the checks pass

Definition at line 123 of file TwoViewCosmicRayRemovalTool.cc.

 {
     // ATTN: Avoid endpoints, topological michel reconstruction is very ambiguous
     if (this->IsMuonEndpoint(element, hitType))
         return false;
 
     const Cluster *pMuonCluster(nullptr), *const pDeltaRayCluster(m_pParentAlgorithm->GetCluster(element, hitType));
 
     if (m_pParentAlgorithm->GetMuonCluster(element.GetOverlapResult().GetCommonMuonPfoList(), hitType, pMuonCluster) != STATUS_CODE_SUCCESS)
         return false;
 
     const float separation(LArClusterHelper::GetClosestDistance(pDeltaRayCluster, pMuonCluster));
 
     return separation <= m_minSeparation;
 }

StatusCode lar_content::TwoViewCosmicRayRemovalTool::ProjectDeltaRayPositions	(	const MatrixType::Element &	element,
		const pandora::HitType	hitType,
		pandora::CartesianPointVector &	projectedPositions
	)		const

private

Use two views of a delta ray pfo to calculate projected positions in a given third view.

Parameters

element	the matrix element
hitType	the view to be projected into
projectedPositions	the output list of projected positions

Returns: a status code reflecting whether the procedure ran smoothly and if the outcome is good

Definition at line 511 of file TwoViewCosmicRayRemovalTool.cc.

 {
     HitTypeVector hitTypes({TPC_VIEW_U, TPC_VIEW_V, TPC_VIEW_W});
 
     hitTypes.erase(std::find(hitTypes.begin(), hitTypes.end(), hitType));
 
     const Cluster *const pCluster1(m_pParentAlgorithm->GetCluster(element, hitTypes[0]));
     const Cluster *const pCluster2(m_pParentAlgorithm->GetCluster(element, hitTypes[1]));
 
     if (!pCluster1 || !pCluster2)
         return STATUS_CODE_NOT_FOUND;
 
     return m_pParentAlgorithm->GetProjectedPositions(pCluster1, pCluster2, projectedPositions);
 }

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

private

Definition at line 529 of file TwoViewCosmicRayRemovalTool.cc.

 {
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MinSeparation", m_minSeparation));
 
     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, "DistanceToLine", m_distanceToLine));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MinContaminationLength", m_minContaminationLength));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MaxDistanceToHit", m_maxDistanceToHit));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MinRemnantClusterSize", m_minRemnantClusterSize));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MaxDistanceToTrack", m_maxDistanceToTrack));
 
     return STATUS_CODE_SUCCESS;
 }

void lar_content::TwoViewCosmicRayRemovalTool::ReclusterRemnant	(	const pandora::HitType	hitType,
		const pandora::Cluster *const	pMuonCluster,
		const pandora::Cluster *const	pDeltaRayRemnant,
		pandora::ClusterVector &	clusterVector,
		pandora::PfoVector &	pfoVector
	)		const

private

Reculster the remnant cluster, merging insignificant created clusters into the parent muon (if proximity checks pass)

Parameters

hitType	the hit type of the cluster under investigation
pMuonCluster	the address of the parent muon cluster
pDeltaRayRemnant	the address of the delta ray remnant
clusterVector	a vector containing the address of created/modified clusters for bookeeping purposes
pfoVector	a vector containing the address of the pfo to which a modified muon cluster belongs for bookeeping purposes

Definition at line 471 of file TwoViewCosmicRayRemovalTool.cc.

 {
     std::string caloHitListName(hitType == TPC_VIEW_U ? "CaloHitListU" : hitType == TPC_VIEW_V ? "CaloHitListV" : "CaloHitListW");
 
     PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::ReplaceCurrentList<CaloHit>(*m_pParentAlgorithm, caloHitListName));
     PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=,
         PandoraContentApi::ReplaceCurrentList<Cluster>(*m_pParentAlgorithm, m_pParentAlgorithm->GetClusterListName(hitType)));
     PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::Delete<Cluster>(*m_pParentAlgorithm, pDeltaRayRemnant));
 
     const ClusterList *pClusterList(nullptr);
     std::string newClusterListName;
     PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=,
         PandoraContentApi::RunClusteringAlgorithm(*m_pParentAlgorithm, m_pParentAlgorithm->GetClusteringAlgName(), pClusterList, newClusterListName));
 
     const ClusterList remnantClusters(pClusterList->begin(), pClusterList->end());
 
     PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=,
         PandoraContentApi::SaveList<Cluster>(*m_pParentAlgorithm, newClusterListName, m_pParentAlgorithm->GetClusterListName(hitType)));
     PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=,
         PandoraContentApi::ReplaceCurrentList<Cluster>(*m_pParentAlgorithm, m_pParentAlgorithm->GetClusterListName(hitType)));
 
     for (const Cluster *const pRemnant : remnantClusters)
     {
         if (pRemnant->GetNCaloHits() < m_minRemnantClusterSize)
         {
             if (LArClusterHelper::GetClosestDistance(pRemnant, pMuonCluster) < m_maxDistanceToTrack)
             {
                 PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::MergeAndDeleteClusters(*m_pParentAlgorithm, pMuonCluster, pRemnant));
                 continue;
             }
         }
 
         clusterVector.push_back(pRemnant);
         pfoVector.push_back(nullptr);
     }
 }

bool lar_content::TwoViewCosmicRayRemovalTool::RemoveCosmicRayHits ( const MatrixType::ElementList & elementList ) const

private

Remove hits from delta ray clusters that belong to the parent muon.

Parameters

elementList the matrix element list

Returns: whether any clusters have been modified

Definition at line 69 of file TwoViewCosmicRayRemovalTool.cc.

 {
     ClusterSet modifiedClusters, checkedClusters;
 
     for (const MatrixType::Element &element : elementList)
     {
         for (const HitType hitType : m_pParentAlgorithm->GetHitTypeVector())
         {
             const Cluster *const pDeltaRayCluster(m_pParentAlgorithm->GetCluster(element, hitType));
 
             if (checkedClusters.count(pDeltaRayCluster))
                 continue;
 
             // ATTN: The underlying matrix will update during this loop, do not proceed if element has been modified
             if ((modifiedClusters.count(element.GetCluster1())) || (modifiedClusters.count(element.GetCluster2())))
                 continue;
 
             if (!this->PassElementChecks(element, hitType))
                 continue;
 
             if (!this->IsContaminated(element, hitType))
                 continue;
 
             if (!this->IsBestElement(element, hitType, elementList))
                 continue;
 
             checkedClusters.insert(pDeltaRayCluster);
 
             // Attempt to pull delta ray hits out of cluster
             CaloHitList deltaRayHits;
             this->CreateSeed(element, hitType, deltaRayHits);
 
             if (deltaRayHits.empty())
                 continue;
 
             // ATTN: If seed cannot be grown, abort
             CaloHitList deltaRayRemnantHits;
             if (this->GrowSeed(element, hitType, deltaRayHits, deltaRayRemnantHits) != STATUS_CODE_SUCCESS)
                 continue;
 
             if (deltaRayHits.size() == pDeltaRayCluster->GetNCaloHits())
                 continue;
 
             modifiedClusters.insert(pDeltaRayCluster);
 
             this->SplitDeltaRayCluster(element, hitType, deltaRayHits, deltaRayRemnantHits);
         }
     }
 
     return !modifiedClusters.empty();
 }

bool lar_content::TwoViewCosmicRayRemovalTool::Run	(	TwoViewDeltaRayMatchingAlgorithm *const	pAlgorithm,
		MatrixType &	matrixTensor
	)

privatevirtual

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

Definition at line 35 of file TwoViewCosmicRayRemovalTool.cc.

 {
     m_pParentAlgorithm = pAlgorithm;
 
     if (PandoraContentApi::GetSettings(*m_pParentAlgorithm)->ShouldDisplayAlgorithmInfo())
         std::cout << "----> Running Algorithm Tool: " << this->GetInstanceName() << ", " << this->GetType() << std::endl;
 
     bool changesMade(false);
 
     ClusterVector sortedKeyClusters;
     overlapMatrix.GetSortedKeyClusters(sortedKeyClusters);
 
     ClusterSet usedKeyClusters;
     for (const Cluster *const pKeyCluster : sortedKeyClusters)
     {
         if (usedKeyClusters.count(pKeyCluster))
             continue;
 
         MatrixType::ElementList elementList;
         overlapMatrix.GetConnectedElements(pKeyCluster, true, elementList);
 
         for (const MatrixType::Element &element : elementList)
             usedKeyClusters.insert(element.GetCluster1());
 
         const bool changesMadeInIteration = this->RemoveCosmicRayHits(elementList);
 
         changesMade = (changesMade ? changesMade : changesMadeInIteration);
     }
 
     return changesMade;
 }

void lar_content::TwoViewCosmicRayRemovalTool::SplitDeltaRayCluster	(	const MatrixType::Element &	element,
		const pandora::HitType	hitType,
		pandora::CaloHitList &	collectedHits,
		pandora::CaloHitList &	deltaRayRemnantHits
	)		const

private

Fragment the delta ray cluster, refining the matched delta ray cluster perhaps creating significant remnants in the process.

Parameters

element	the matrix element
hitType	the hit type of the cluster under investigation
collectedHits	the list of identified delta ray hits
deltaRayRemantHits	the list of remainder hits

Definition at line 411 of file TwoViewCosmicRayRemovalTool.cc.

 {
     const Cluster *const pDeltaRayCluster(m_pParentAlgorithm->GetCluster(element, hitType)), *pMuonCluster(nullptr);
 
     if (m_pParentAlgorithm->GetMuonCluster(element.GetOverlapResult().GetCommonMuonPfoList(), hitType, pMuonCluster) != STATUS_CODE_SUCCESS)
         return;
 
     m_pParentAlgorithm->UpdateUponDeletion(pMuonCluster);
     m_pParentAlgorithm->UpdateUponDeletion(pDeltaRayCluster);
 
     std::string originalListName, fragmentListName;
     ClusterList originalClusterList(1, pDeltaRayCluster);
 
     PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=,
         PandoraContentApi::ReplaceCurrentList<Cluster>(*m_pParentAlgorithm, m_pParentAlgorithm->GetClusterListName(hitType)));
     PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=,
         PandoraContentApi::InitializeFragmentation(*m_pParentAlgorithm, originalClusterList, originalListName, fragmentListName));
 
     CaloHitList deltaRayHitList;
     pDeltaRayCluster->GetOrderedCaloHitList().FillCaloHitList(deltaRayHitList);
 
     const Cluster *pDeltaRay(nullptr), *pDeltaRayRemnant(nullptr);
 
     for (const CaloHit *const pCaloHit : deltaRayHitList)
     {
         const bool isDeltaRay(std::find(collectedHits.begin(), collectedHits.end(), pCaloHit) != collectedHits.end());
         const bool isDeltaRayRemnant(
             !isDeltaRay && (std::find(deltaRayRemnantHits.begin(), deltaRayRemnantHits.end(), pCaloHit) != deltaRayRemnantHits.end()));
         const Cluster *&pCluster(isDeltaRay ? pDeltaRay : isDeltaRayRemnant ? pDeltaRayRemnant : pMuonCluster);
 
         if (!pCluster)
         {
             PandoraContentApi::Cluster::Parameters parameters;
             parameters.m_caloHitList.push_back(pCaloHit);
             PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::Cluster::Create(*m_pParentAlgorithm, parameters, pCluster));
         }
         else
         {
             PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::AddToCluster(*m_pParentAlgorithm, pCluster, pCaloHit));
         }
     }
 
     PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::EndFragmentation(*m_pParentAlgorithm, fragmentListName, originalListName));
 
     ClusterVector clusterVector;
     PfoVector pfoVector;
     if (pDeltaRayRemnant)
         this->ReclusterRemnant(hitType, pMuonCluster, pDeltaRayRemnant, clusterVector, pfoVector);
 
     clusterVector.push_back(pMuonCluster);
     pfoVector.push_back(element.GetOverlapResult().GetCommonMuonPfoList().front());
     clusterVector.push_back(pDeltaRay);
     pfoVector.push_back(nullptr);
 
     m_pParentAlgorithm->UpdateForNewClusters(clusterVector, pfoVector);
 }