lar_content::LArClusterHelper Class Reference

LArClusterHelper class. More...

#include "LArClusterHelper.h"

Static Public Member Functions
static pandora::HitType	GetClusterHitType (const pandora::Cluster *const pCluster)
	Get the hit type associated with a two dimensional cluster. More...
static void	GetClustersUVW (const pandora::ClusterList &inputClusters, pandora::ClusterList &clusterListU, pandora::ClusterList &clusterListV, pandora::ClusterList &clusterListW)
	Divide an input cluster list into separate u, v and w lists (exception raised if alternative hit type encountered) More...
static void	GetClustersByHitType (const pandora::ClusterList &inputClusters, const pandora::HitType hitType, pandora::ClusterList &clusterList)
	Get the subset of clusters, from a provided list, that match the specified hit type. More...
static float	GetLengthSquared (const pandora::Cluster *const pCluster)
	Get length squared of cluster. More...
static float	GetLength (const pandora::Cluster *const pCluster)
	Get length of cluster. More...
static float	GetEnergyFromLength (const pandora::Cluster *const pCluster)
	Get energy of cluster, based on length. More...
static unsigned int	GetLayerSpan (const pandora::Cluster *const pCluster)
	Get number of layers spanned by cluster (1+Last-First) More...
static float	GetLayerOccupancy (const pandora::Cluster *const pCluster)
	Fraction of occupied layers in cluster. More...
static float	GetLayerOccupancy (const pandora::Cluster *const pCluster1, const pandora::Cluster *const pCluster2)
	Fraction of occupied layers in a pair of clusters. More...
static float	GetClosestDistance (const pandora::ClusterList &clusterList1, const pandora::ClusterList &clusterList2)
	Get closest distance between clusters in a pair of cluster lists. More...
static float	GetClosestDistance (const pandora::Cluster *const pCluster, const pandora::ClusterList &clusterList)
	Get closest distance between a specified cluster and list of clusters. More...
static float	GetClosestDistance (const pandora::Cluster *const pCluster1, const pandora::Cluster *const pCluster2)
	Get closest distance between a pair of clusters. More...
static float	GetClosestDistance (const pandora::CartesianVector &position, const pandora::ClusterList &clusterList)
	Get closest distance between a specified position and list of clusters. More...
static float	GetClosestDistance (const pandora::CartesianVector &position, const pandora::Cluster *const pCluster)
	Get closest distance between a specified position vector and the hits in a specified cluster. More...
static pandora::CartesianVector	GetClosestPosition (const pandora::CartesianVector &position, const pandora::ClusterList &clusterList)
	Get closest position in a list of clusters to a specified input position vector. More...
static pandora::CartesianVector	GetClosestPosition (const pandora::CartesianVector &position, const pandora::Cluster *const pCluster)
	Get closest position on a cluster to a specified input position vector. More...
static void	GetClosestPositions (const pandora::Cluster *const pCluster1, const pandora::Cluster *const pCluster2, pandora::CartesianVector &position1, pandora::CartesianVector &position2)
	Get pair of closest positions for a pair of clusters. More...
static void	GetExtremalCoordinates (const pandora::ClusterList &clusterList, pandora::CartesianVector &innerCoordinate, pandora::CartesianVector &outerCoordinate)
	Get positions of the two most distant calo hits in a list of cluster (ordered by Z) More...
static void	GetExtremalCoordinates (const pandora::Cluster *const pCluster, pandora::CartesianVector &innerCoordinate, pandora::CartesianVector &outerCoordinate)
	Get positions of the two most distant calo hits in a cluster (ordered by Z) More...
static void	GetExtremalCoordinates (const pandora::OrderedCaloHitList &orderedCaloHitList, pandora::CartesianVector &innerCoordinate, pandora::CartesianVector &outerCoordinate)
	Get positions of the two most distant calo hits in an ordered calo hit list (ordered by Z) More...
static void	GetExtremalCoordinates (const pandora::CartesianPointVector &coordinateVector, pandora::CartesianVector &innerCoordinate, pandora::CartesianVector &outerCoordinate)
	Get positions of the two most distant points in a provided list (ordered by Z) More...
static void	GetClusterBoundingBox (const pandora::Cluster *const pCluster, pandora::CartesianVector &minimumCoordinate, pandora::CartesianVector &maximumCoordinate)
	Get minimum and maximum X, Y and Z positions of the calo hits in a cluster. More...
static void	GetClusterSpanX (const pandora::Cluster *const pCluster, float &xmin, float &xmax)
	Get minimum and maximum X positions of the calo hits in a cluster. More...
static void	GetClusterSpanZ (const pandora::Cluster *const pCluster, const float xmin, const float xmax, float &zmin, float &zmax)
	Get upper and lower Z positions of the calo hits in a cluster in range xmin to xmax. More...
static void	GetCoordinateVector (const pandora::Cluster *const pCluster, pandora::CartesianPointVector &coordinateVector)
	Get vector of hit coordinates from an input cluster. More...
static pandora::StatusCode	GetAverageZ (const pandora::Cluster *const pCluster, const float xmin, const float xmax, float &averageZ)
	Get average Z positions of the calo hits in a cluster in range xmin to xmax. More...
static bool	SortByNOccupiedLayers (const pandora::Cluster *const pLhs, const pandora::Cluster *const pRhs)
	Sort clusters by number of occupied layers, and by inner layer, then energy in event of a tie. More...
static bool	SortByNHits (const pandora::Cluster *const pLhs, const pandora::Cluster *const pRhs)
	Sort clusters by number of hits, then layer span, then inner layer, then position, then pulse-height. More...
static bool	SortByLayerSpan (const pandora::Cluster *const pLhs, const pandora::Cluster *const pRhs)
	Sort clusters by layer span, then inner layer, then position, then pulse-height. More...
static bool	SortByInnerLayer (const pandora::Cluster *const pLhs, const pandora::Cluster *const pRhs)
	Sort clusters by inner layer, then position, then pulse-height. More...
static bool	SortByPosition (const pandora::Cluster *const pLhs, const pandora::Cluster *const pRhs)
	Sort clusters by position, then pulse-height. More...
static bool	SortByPulseHeight (const pandora::Cluster *const pLhs, const pandora::Cluster *const pRhs)
	Sort clusters by pulse-height. More...
static bool	SortHitsByPosition (const pandora::CaloHit *const pLhs, const pandora::CaloHit *const pRhs)
	Sort calo hits by their position (use Z, followed by X, followed by Y) More...
static bool	SortHitsByPulseHeight (const pandora::CaloHit *const pLhs, const pandora::CaloHit *const pRhs)
	Sort calo hits by their pulse height. More...
static bool	SortCoordinatesByPosition (const pandora::CartesianVector &lhs, const pandora::CartesianVector &rhs)
	Sort cartesian vectors by their position (use Z, followed by X, followed by Y) More...

Detailed Description

LArClusterHelper class.

Definition at line 19 of file LArClusterHelper.h.

Member Function Documentation

StatusCode lar_content::LArClusterHelper::GetAverageZ ( const pandora::Cluster *const  pCluster, const float  xmin, const float  xmax, float &  averageZ  )

static

Get average Z positions of the calo hits in a cluster in range xmin to xmax.

Parameters

pCluster	address of the cluster
xmin	for range in x
xmax	for range in x
averageZ	to receive the average Z position

Returns

status code, faster than throwing in regular use-cases

Definition at line 403 of file LArClusterHelper.cc.

References f, max, and min.

Referenced by lar_content::ThreeDRemnantsAlgorithm::CalculateOverlapResult(), lar_content::ParticleRecoveryAlgorithm::CheckConsistency(), and lar_content::CosmicRayShowerMatchingAlgorithm::CheckMatchedClusters3D().

{
    averageZ = std::numeric_limits<float>::max();

    if (xmin > xmax)
        return STATUS_CODE_INVALID_PARAMETER;

    const OrderedCaloHitList &orderedCaloHitList(pCluster->GetOrderedCaloHitList());

    float zsum(0.f);
    int count(0);

    for (OrderedCaloHitList::const_iterator ochIter = orderedCaloHitList.begin(), ochIterEnd = orderedCaloHitList.end(); ochIter != ochIterEnd; ++ochIter)
    {
        for (CaloHitList::const_iterator hIter = ochIter->second->begin(), hIterEnd = ochIter->second->end(); hIter != hIterEnd; ++hIter)
        {
            const CaloHit *const pCaloHit = *hIter;
            const CartesianVector &hit(pCaloHit->GetPositionVector());

            if (hit.GetX() < xmin || hit.GetX() > xmax)
                continue;

            zsum += hit.GetZ();
            ++count;
        }
    }

    if (count == 0)
        return STATUS_CODE_NOT_FOUND;

    averageZ = zsum / static_cast<float>(count);
    return STATUS_CODE_SUCCESS;
}

static float lar_content::LArClusterHelper::GetClosestDistance ( const pandora::ClusterList &  clusterList1, const pandora::ClusterList &  clusterList2  )

static

Get closest distance between clusters in a pair of cluster lists.

Parameters

clusterList1	the first cluster list
clusterList2	the second cluster list

Returns

the closest distance

Referenced by lar_content::ShowerGrowingAlgorithm::AreClustersAssociated(), lar_content::NeutrinoDaughterVerticesAlgorithm::BuildDaughterTrack(), lar_content::OvershootSplittingAlgorithm::BuildIntersectionMap(), lar_content::SplitShowersTool::CheckClusterProximities(), lar_content::CosmicRayTrackMatchingAlgorithm::CheckMatchedClusters3D(), lar_content::NearbyClusterMopUpAlgorithm::ClusterMopUp(), lar_content::CosmicRaySplittingAlgorithm::ConfirmSplitPosition(), lar_content::DeltaRaySplittingAlgorithm::FindBestSplitPosition(), lar_content::CrossedTrackSplittingAlgorithm::FindBestSplitPosition(), lar_content::DeltaRayIdentificationAlgorithm::GetClosestDistance(), lar_content::DeltaRayMatchingAlgorithm::GetDistanceSquaredToPfo(), lar_content::GlobalAsymmetryFeatureTool::GetGlobalAsymmetryForView(), lar_content::TransverseExtensionAlgorithm::GetListOfCleanClusters(), lar_content::DeltaRayGrowingAlgorithm::GetListOfSeedClusters(), lar_content::LocalAsymmetryFeatureTool::GetLocalAsymmetryForView(), lar_content::CutClusterCharacterisationAlgorithm::GetVertexDistance(), lar_content::SimpleClusterMergingAlgorithm::IsAssociated(), lar_content::ConnectedRemnantsTool::IsConnected(), lar_content::ClusterGrowingAlgorithm::PopulateClusterMergeMap(), lar_content::BranchAssociatedPfosTool::Run(), lar_content::CosmicRayTrackMatchingAlgorithm::SelectCleanClusters(), and lar_content::ShowerAsymmetryFeatureTool::ShouldUseShowerCluster().

static float lar_content::LArClusterHelper::GetClosestDistance ( const pandora::Cluster *const  pCluster, const pandora::ClusterList &  clusterList  )

static

Get closest distance between a specified cluster and list of clusters.

Parameters

pCluster	address of the input cluster
clusterList	list of input clusters

Returns

the closest distance

static float lar_content::LArClusterHelper::GetClosestDistance ( const pandora::Cluster *const  pCluster1, const pandora::Cluster *const  pCluster2  )

static

Get closest distance between a pair of clusters.

Parameters

pCluster1	address of the first cluster
pCluster2	address of the second cluster

Returns

the closest distance

static float lar_content::LArClusterHelper::GetClosestDistance ( const pandora::CartesianVector &  position, const pandora::ClusterList &  clusterList  )

static

Get closest distance between a specified position and list of clusters.

Parameters

position	the position vector
clusterList	list of input clusters

Returns

the closest distance

static float lar_content::LArClusterHelper::GetClosestDistance ( const pandora::CartesianVector &  position, const pandora::Cluster *const  pCluster  )

static

Get closest distance between a specified position vector and the hits in a specified cluster.

Parameters

position	the position vector
pCluster	address of the cluster

Returns

the closest distance

static pandora::CartesianVector lar_content::LArClusterHelper::GetClosestPosition ( const pandora::CartesianVector &  position, const pandora::ClusterList &  clusterList  )

static

Get closest position in a list of clusters to a specified input position vector.

Parameters

position	the position vector
clusterList	list of input clusters

Returns

the closest position

Referenced by lar_content::NeutrinoDaughterVerticesAlgorithm::BuildDaughterShower(), lar_content::CrossedTrackSplittingAlgorithm::FindCandidateSplitPositions(), lar_content::CosmicRaySplittingAlgorithm::GetCaloHitListToMove(), lar_content::DeltaRayExtensionAlgorithm::GetExtremalCoordinatesFromCache(), lar_content::CosmicRaySplittingAlgorithm::IdentifyCrossedTracks(), lar_content::TransverseAssociationAlgorithm::IsAssociated(), lar_content::TransverseAssociationAlgorithm::IsOverlapping(), and lar_content::TrackConsolidationAlgorithm::TrackConsolidationAlgorithm().

static pandora::CartesianVector lar_content::LArClusterHelper::GetClosestPosition ( const pandora::CartesianVector &  position, const pandora::Cluster *const  pCluster  )

static

Get closest position on a cluster to a specified input position vector.

Parameters

position	the position vector
pCluster	address of the cluster

Returns

the closest position

void lar_content::LArClusterHelper::GetClosestPositions ( const pandora::Cluster *const  pCluster1, const pandora::Cluster *const  pCluster2, pandora::CartesianVector &  position1, pandora::CartesianVector &  position2  )

static

Get pair of closest positions for a pair of clusters.

Parameters

pCluster1	the address of the first cluster
pCluster2	the address of the second cluster
the	closest position in the first cluster
the	closest position in the second cluster

Definition at line 267 of file LArClusterHelper.cc.

References f, and max.

Referenced by lar_content::CosmicRayVertexBuildingAlgorithm::BuildCosmicRayDaughter(), lar_content::NeutrinoDaughterVerticesAlgorithm::BuildDaughterShower(), and lar_content::EndAssociatedPfosTool::IsCloseToParentEndpoint().

{
    bool distanceFound(false);
    float minDistanceSquared(std::numeric_limits<float>::max());

    CartesianVector closestPosition1(0.f, 0.f, 0.f);
    CartesianVector closestPosition2(0.f, 0.f, 0.f);

    const OrderedCaloHitList &orderedCaloHitList1(pCluster1->GetOrderedCaloHitList());
    const OrderedCaloHitList &orderedCaloHitList2(pCluster2->GetOrderedCaloHitList());

    // Loop over hits in cluster 1
    for (OrderedCaloHitList::const_iterator iter1 = orderedCaloHitList1.begin(), iter1End = orderedCaloHitList1.end(); iter1 != iter1End; ++iter1)
    {
        for (CaloHitList::const_iterator hitIter1 = iter1->second->begin(), hitIter1End = iter1->second->end(); hitIter1 != hitIter1End; ++hitIter1)
        {
            const CartesianVector &positionVector1((*hitIter1)->GetPositionVector());

            // Loop over hits in cluster 2
            for (OrderedCaloHitList::const_iterator iter2 = orderedCaloHitList2.begin(), iter2End = orderedCaloHitList2.end(); iter2 != iter2End; ++iter2)
            {
                for (CaloHitList::const_iterator hitIter2 = iter2->second->begin(), hitIter2End = iter2->second->end(); hitIter2 != hitIter2End; ++hitIter2)
                { 
                    const CartesianVector &positionVector2((*hitIter2)->GetPositionVector());

                    const float distanceSquared((positionVector1 - positionVector2).GetMagnitudeSquared());

                    if (distanceSquared < minDistanceSquared)
                    {
                        minDistanceSquared = distanceSquared;
                        closestPosition1 = positionVector1;
                        closestPosition2 = positionVector2;
                        distanceFound = true;
                    }
                }
            }
        }
    }

    if (!distanceFound)
        throw StatusCodeException(STATUS_CODE_NOT_FOUND);

    outputPosition1 = closestPosition1;
    outputPosition2 = closestPosition2;
}

void lar_content::LArClusterHelper::GetClusterBoundingBox ( const pandora::Cluster *const  pCluster, pandora::CartesianVector &  minimumCoordinate, pandora::CartesianVector &  maximumCoordinate  )

static

Get minimum and maximum X, Y and Z positions of the calo hits in a cluster.

Parameters

pCluster	address of the cluster
the	minimum positions (x,y,z)
the	maximum positions (x,y,z)

Definition at line 316 of file LArClusterHelper.cc.

References max, and min.

Referenced by lar_content::ThreeDTrackFragmentsAlgorithm::CheckMatchedClusters(), lar_content::SvmVertexSelectionAlgorithm::GetEventShapeFeatures(), and lar_content::CosmicRayTrackRecoveryAlgorithm::SelectCleanClusters().

{
    const OrderedCaloHitList &orderedCaloHitList(pCluster->GetOrderedCaloHitList());

    float xmin(std::numeric_limits<float>::max());
    float ymin(std::numeric_limits<float>::max());
    float zmin(std::numeric_limits<float>::max());
    float xmax(-std::numeric_limits<float>::max());
    float ymax(-std::numeric_limits<float>::max());
    float zmax(-std::numeric_limits<float>::max());

    for (OrderedCaloHitList::const_iterator ochIter = orderedCaloHitList.begin(), ochIterEnd = orderedCaloHitList.end(); ochIter != ochIterEnd; ++ochIter)
    {
        for (CaloHitList::const_iterator hIter = ochIter->second->begin(), hIterEnd = ochIter->second->end(); hIter != hIterEnd; ++hIter)
        {
            const CaloHit *const pCaloHit = *hIter;
            const CartesianVector &hit(pCaloHit->GetPositionVector());
            xmin = std::min(hit.GetX(), xmin);
            xmax = std::max(hit.GetX(), xmax);
            ymin = std::min(hit.GetY(), ymin);
            ymax = std::max(hit.GetY(), ymax);
            zmin = std::min(hit.GetZ(), zmin);
            zmax = std::max(hit.GetZ(), zmax);
        }
    }

    minimumCoordinate.SetValues(xmin, ymin, zmin);
    maximumCoordinate.SetValues(xmax, ymax, zmax);
}

HitType lar_content::LArClusterHelper::GetClusterHitType ( const pandora::Cluster *const  pCluster )

static

Get the hit type associated with a two dimensional cluster.

Parameters

pCluster

the address of the cluster

Returns

the cluster hit type

Definition at line 20 of file LArClusterHelper.cc.

Referenced by lar_content::ParticleRecoveryAlgorithm::SimpleOverlapTensor::AddAssociation(), lar_content::DeltaRayMatchingAlgorithm::AddToDaughterPfo(), lar_content::DeltaRayMatchingAlgorithm::AreClustersMatched(), lar_content::EventSlicingTool::AssignRemainingHitsToSlices(), lar_content::CrossGapsExtensionAlgorithm::BuildPointingClusterList(), lar_content::TrackHitsBaseTool::BuildSlidingFitMap(), lar_content::TwoDLinearFitFeatureTool::CalculateVariablesSlidingLinearFit(), lar_content::ThreeDLinearFitFeatureTool::CalculateVariablesSlidingLinearFit(), lar_content::SplitShowersTool::CheckClusterVertexRelations(), lar_content::CosmicRayShowerMatchingAlgorithm::CheckMatchedClusters3D(), lar_content::CosmicRayTrackMatchingAlgorithm::CheckMatchedClusters3D(), lar_content::TracksCrossingGapsTool::CheckXPositionInGap(), lar_content::ConeClusterMopUpAlgorithm::ClusterMopUp(), lar_content::NearbyClusterMopUpAlgorithm::ClusterMopUp(), lar_content::EventSlicingTool::CopyPfoHitsToSlices(), lar_content::CandidateVertexCreationAlgorithm::CreateEndpointCandidates(), lar_content::CandidateVertexCreationAlgorithm::CreateEndpointVertex(), lar_content::EventSlicingTool::CreateSlices(), lar_content::CrossGapsExtensionAlgorithm::FillClusterAssociationMatrix(), lar_content::VertexSplittingAlgorithm::FindBestSplitPosition(), lar_content::ShowerGrowingAlgorithm::GetAllVertexSeedCandidates(), lar_content::VertexBasedPfoRecoveryAlgorithm::GetAvailableClusters(), lar_content::VertexBasedPfoMopUpAlgorithm::GetClusterAssociation(), lar_content::ClusterMopUpBaseAlgorithm::GetClusterLists(), lar_content::VertexSelectionBaseAlgorithm::GetClusterLists(), lar_content::SlidingConeClusterMopUpAlgorithm::GetClusterMergeMap(), lar_content::ParticleRecoveryAlgorithm::SimpleOverlapTensor::GetConnectedElements(), lar_content::DeltaRayMatchingAlgorithm::GetDistanceSquaredToPfo(), lar_content::ShowerGrowingAlgorithm::GetFigureOfMerit(), lar_content::ParticleRecoveryAlgorithm::GetInputClusters(), lar_content::VertexBasedPfoMopUpAlgorithm::GetInputPfos(), lar_content::DeltaRayGrowingAlgorithm::GetListOfSeedClusters(), lar_content::CheatingPfoCreationAlgorithm::GetNHitTypesAboveThreshold(), lar_content::NeutrinoPropertiesAlgorithm::GetNTwoDHitsInPfoChain(), lar_pandora::LArPandoraOutput::GetPandoraToArtHitMap(), lar_content::VertexBasedPfoMopUpAlgorithm::GetPfoAssociations(), lar_content::ThreeDTrackFragmentsAlgorithm::GetProjectedPositions(), lar_content::CutClusterCharacterisationAlgorithm::GetVertexDistance(), lar_content::CrossGapsAssociationAlgorithm::IsAssociated(), lar_content::PfoCharacterisationBaseAlgorithm::IsClearTrack3x2D(), lar_content::ParticleRecoveryAlgorithm::IsOverlap(), lar_content::UndershootTracksTool::IsThreeDKink(), lar_content::OvershootTracksTool::IsThreeDKink(), lar_content::ThreeDBaseAlgorithm< T >::MakeClusterMerges(), lar_content::SlidingConeClusterMopUpAlgorithm::MakeClusterMerges(), lar_content::ThreeDTracksBaseAlgorithm< T >::MakeClusterSplits(), lar_content::CosmicRayTrackRecoveryAlgorithm::MatchClusters(), lar_content::VertexBasedPfoRecoveryAlgorithm::MatchThreeViews(), lar_content::VertexBasedPfoRecoveryAlgorithm::MatchTwoViews(), lar_content::ThreeDChargeFeatureTool::OrderCaloHitsByDistanceToVertex(), lar_content::CosmicRayBaseMatchingAlgorithm::Particle::Particle(), lar_content::MissingTrackSegmentTool::Particle::Particle(), lar_content::VertexBasedPfoRecoveryAlgorithm::Particle::Particle(), lar_content::DeltaRayMatchingAlgorithm::Particle::Particle(), lar_content::TrackSplittingTool::PassesChecks(), lar_content::CandidateVertexCreationAlgorithm::SelectClusters(), lar_content::VertexBasedPfoRecoveryAlgorithm::SelectClusters(), lar_content::VertexBasedPfoRecoveryAlgorithm::SelectVertexClusters(), lar_content::ThreeDHitCreationAlgorithm::SeparateTwoDHits(), lar_content::OvershootTracksTool::SetSplitPosition(), lar_content::MasterAlgorithm::StitchPfos(), lar_content::ThreeDTrackFragmentsAlgorithm::UpdateForNewCluster(), and lar_content::ThreeDBaseAlgorithm< T >::UpdateForNewCluster().

{
    if (0 == pCluster->GetNCaloHits())
        throw StatusCodeException(STATUS_CODE_NOT_INITIALIZED);

    // TODO Find a way of confirming single hit-type clustering mode; currently only checked in ListPreparation algorithm
    // if (!pandora->GetSettings()->SingleHitTypeClusteringMode())
    //     throw StatusCodeException(STATUS_CODE_FAILURE);

    return (*(pCluster->GetOrderedCaloHitList().begin()->second->begin()))->GetHitType();
}

void lar_content::LArClusterHelper::GetClustersByHitType ( const pandora::ClusterList &  inputClusters, const pandora::HitType  hitType, pandora::ClusterList &  clusterList  )

static

Get the subset of clusters, from a provided list, that match the specified hit type.

Parameters

inputClusters	the input cluster list
hitType	the specified hit type
clusterList	to receive the clusters

Definition at line 49 of file LArClusterHelper.cc.

Referenced by lar_content::CosmicRayTrackRecoveryAlgorithm::MergeClusters().

{
    for (ClusterList::const_iterator iter = inputClusters.begin(), iterEnd = inputClusters.end(); iter != iterEnd; ++iter)
    {
        if (hitType == LArClusterHelper::GetClusterHitType(*iter))
            clusterList.push_back(*iter);
    }
}

void lar_content::LArClusterHelper::GetClusterSpanX ( const pandora::Cluster *const  pCluster, float &  xmin, float &  xmax  )

static

Get minimum and maximum X positions of the calo hits in a cluster.

Parameters

pCluster	address of the cluster
the	minimum position of x
the	maximum position of x

Definition at line 348 of file LArClusterHelper.cc.

References max, and min.

Referenced by lar_content::DeltaRayMatchingAlgorithm::AreClustersMatched(), lar_content::ThreeDRemnantsAlgorithm::CalculateOverlapResult(), lar_content::ParticleRecoveryAlgorithm::CheckConsistency(), lar_content::CosmicRayTaggingTool::CheckIfInTime(), lar_content::CosmicRayShowerMatchingAlgorithm::CheckMatchedClusters3D(), lar_content::ThreeDTrackFragmentsAlgorithm::GetProjectedPositions(), lar_content::ParticleRecoveryAlgorithm::IsOverlap(), lar_content::CosmicRayShowerMatchingAlgorithm::MatchClusters(), lar_content::CosmicRayTrackMatchingAlgorithm::MatchClusters(), lar_content::ThreeDTrackFragmentsAlgorithm::PerformMainLoop(), and lar_content::ParticleRecoveryAlgorithm::StandardClusterSelection().

{
    const OrderedCaloHitList &orderedCaloHitList(pCluster->GetOrderedCaloHitList());
    xmin = std::numeric_limits<float>::max();
    xmax = -std::numeric_limits<float>::max();

    for (OrderedCaloHitList::const_iterator ochIter = orderedCaloHitList.begin(), ochIterEnd = orderedCaloHitList.end(); ochIter != ochIterEnd; ++ochIter)
    {
        for (CaloHitList::const_iterator hIter = ochIter->second->begin(), hIterEnd = ochIter->second->end(); hIter != hIterEnd; ++hIter)
        {
            const CaloHit *const pCaloHit = *hIter;
            const CartesianVector &hit(pCaloHit->GetPositionVector());
            xmin = std::min(hit.GetX(), xmin);
            xmax = std::max(hit.GetX(), xmax);
        }
    }
}

void lar_content::LArClusterHelper::GetClusterSpanZ ( const pandora::Cluster *const  pCluster, const float  xmin, const float  xmax, float &  zmin, float &  zmax  )

static

Get upper and lower Z positions of the calo hits in a cluster in range xmin to xmax.

Parameters

pCluster	address of the cluster
xmin	for range in x
xmax	for range in x
zmin	the lower z for this range of x
zmax	the upper z for this range in x

Definition at line 368 of file LArClusterHelper.cc.

References max, and min.

Referenced by lar_content::DeltaRayMatchingAlgorithm::AreClustersMatched().

{
    if (xmin > xmax)
        throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);

    const OrderedCaloHitList &orderedCaloHitList(pCluster->GetOrderedCaloHitList());

    zmin = std::numeric_limits<float>::max();
    zmax = -std::numeric_limits<float>::max();

    bool foundHits(false);
    
    for (OrderedCaloHitList::const_iterator ochIter = orderedCaloHitList.begin(), ochIterEnd = orderedCaloHitList.end(); ochIter != ochIterEnd; ++ochIter)
    {
        for (CaloHitList::const_iterator hIter = ochIter->second->begin(), hIterEnd = ochIter->second->end(); hIter != hIterEnd; ++hIter)
        {
            const CaloHit *const pCaloHit = *hIter;
            const CartesianVector &hit(pCaloHit->GetPositionVector());

            if (hit.GetX() < xmin || hit.GetX() > xmax)
                continue;

            zmin = std::min(hit.GetZ(), zmin);
            zmax = std::max(hit.GetZ(), zmax);
            foundHits = true;
        }
    }

    if (!foundHits)
        throw StatusCodeException(STATUS_CODE_NOT_FOUND);
}

void lar_content::LArClusterHelper::GetClustersUVW ( const pandora::ClusterList &  inputClusters, pandora::ClusterList &  clusterListU, pandora::ClusterList &  clusterListV, pandora::ClusterList &  clusterListW  )

static

Divide an input cluster list into separate u, v and w lists (exception raised if alternative hit type encountered)

Parameters

inputClusters	the input cluster list
clusterListU	to receive the u clusters
clusterListV	to receive the v clusters
clusterListW	to receive the w clusters

Definition at line 34 of file LArClusterHelper.cc.

{
    for (ClusterList::const_iterator iter = inputClusters.begin(), iterEnd = inputClusters.end(); iter != iterEnd; ++iter)
    {
        const HitType hitType(LArClusterHelper::GetClusterHitType(*iter));

        if (TPC_VIEW_U == hitType) clusterListU.push_back(*iter);
        else if (TPC_VIEW_V == hitType) clusterListV.push_back(*iter);
        else if (TPC_VIEW_W == hitType) clusterListW.push_back(*iter);
        else throw StatusCodeException(STATUS_CODE_NOT_FOUND);
    }
}

void lar_content::LArClusterHelper::GetCoordinateVector ( const pandora::Cluster *const  pCluster, pandora::CartesianPointVector &  coordinateVector  )

static

Get vector of hit coordinates from an input cluster.

Parameters

pCluster	address of the cluster
coordinateVector

Definition at line 590 of file LArClusterHelper.cc.

Referenced by lar_content::VertexSelectionBaseAlgorithm::ShowerCluster::GetClusterListCoordinateVector(), lar_content::ThreeDSlidingFitResult::GetGlobalDirection(), lar_content::TwoDSlidingShowerFitResult::GetShowerEdges(), lar_content::CandidateVertexCreationAlgorithm::GetSpacepoints(), and lar_content::TwoDSlidingFitResult::TwoDSlidingFitResult().

{
    for (const OrderedCaloHitList::value_type &layerEntry : pCluster->GetOrderedCaloHitList())
    {
        for (const CaloHit *const pCaloHit : *layerEntry.second)
            coordinateVector.push_back(pCaloHit->GetPositionVector());
    }

    std::sort(coordinateVector.begin(), coordinateVector.end(), LArClusterHelper::SortCoordinatesByPosition);
}

float lar_content::LArClusterHelper::GetEnergyFromLength ( const pandora::Cluster *const  pCluster )

static

Get energy of cluster, based on length.

Parameters

pCluster

address of the cluster

Returns

the energy

Definition at line 99 of file LArClusterHelper.cc.

References f.

{
    const float dEdX(0.002f); // approximately 2 MeV/cm
    return (dEdX * LArClusterHelper::GetLength(pCluster));
}

static void lar_content::LArClusterHelper::GetExtremalCoordinates ( const pandora::ClusterList &  clusterList, pandora::CartesianVector &  innerCoordinate, pandora::CartesianVector &  outerCoordinate  )

static

Get positions of the two most distant calo hits in a list of cluster (ordered by Z)

Parameters

clusterList	the input cluster list
the	inner extremal position
the	outer extremal position

Referenced by lar_content::CosmicRayVertexBuildingAlgorithm::BuildCosmicRayParent(), lar_content::NeutrinoDaughterVerticesAlgorithm::BuildDaughterTrack(), lar_content::CosmicRayTrackMatchingAlgorithm::CheckMatchedClusters3D(), lar_content::DeltaRayExtensionAlgorithm::GetExtremalCoordinatesFromCache(), lar_content::TransverseAssociationAlgorithm::GetExtremalCoordinatesXZ(), lar_content::TransverseExtensionAlgorithm::GetListOfCleanClusters(), lar_content::DeltaRayIdentificationAlgorithm::GetTwoDVertices(), lar_content::CosmicRaySplittingAlgorithm::IdentifyCrossedTracks(), lar_content::CosmicRayTrackMatchingAlgorithm::MatchClusters(), lar_content::CosmicRayTrackMatchingAlgorithm::SelectCleanClusters(), and lar_content::ShowerGrowingAlgorithm::SortClusters().

static void lar_content::LArClusterHelper::GetExtremalCoordinates ( const pandora::Cluster *const  pCluster, pandora::CartesianVector &  innerCoordinate, pandora::CartesianVector &  outerCoordinate  )

static

Get positions of the two most distant calo hits in a cluster (ordered by Z)

Parameters

pCluster	the input cluster
the	inner extremal position
the	outer extremal position

static void lar_content::LArClusterHelper::GetExtremalCoordinates ( const pandora::OrderedCaloHitList &  orderedCaloHitList, pandora::CartesianVector &  innerCoordinate, pandora::CartesianVector &  outerCoordinate  )

static

Get positions of the two most distant calo hits in an ordered calo hit list (ordered by Z)

Parameters

orderedCaloHitList	the ordered calo hit list
the	inner extremal position
the	outer extremal position

static void lar_content::LArClusterHelper::GetExtremalCoordinates ( const pandora::CartesianPointVector &  coordinateVector, pandora::CartesianVector &  innerCoordinate, pandora::CartesianVector &  outerCoordinate  )

static

Get positions of the two most distant points in a provided list (ordered by Z)

Parameters

coordinateVector	the hit list
the	inner extremal position
the	outer extremal position

static float lar_content::LArClusterHelper::GetLayerOccupancy ( const pandora::Cluster *const  pCluster )

static

Fraction of occupied layers in cluster.

Parameters

pCluster

address of the cluster

Returns

float

Referenced by lar_content::LongitudinalExtensionAlgorithm::FillClusterAssociationMatrix(), lar_content::LongitudinalExtensionAlgorithm::GetListOfCleanClusters(), and lar_content::LArParticleIdPlugins::LArMuonId::LArMuonId().

static float lar_content::LArClusterHelper::GetLayerOccupancy ( const pandora::Cluster *const  pCluster1, const pandora::Cluster *const  pCluster2  )

static

Fraction of occupied layers in a pair of clusters.

Parameters

pCluster1	address of the first cluster
pCluster2	address of the second cluster

Returns

float

unsigned int lar_content::LArClusterHelper::GetLayerSpan ( const pandora::Cluster *const  pCluster )

static

Get number of layers spanned by cluster (1+Last-First)

Parameters

pCluster

address of the cluster

Returns

the layer span

Definition at line 107 of file LArClusterHelper.cc.

References f, max, and min.

{
    return (1 + pCluster->GetOuterPseudoLayer() - pCluster->GetInnerPseudoLayer());
}

float lar_content::LArClusterHelper::GetLength ( const pandora::Cluster *const  pCluster )

static

Get length of cluster.

Parameters

pCluster

address of the cluster

Returns

the length

Definition at line 92 of file LArClusterHelper.cc.

Referenced by lar_content::SvmVertexSelectionAlgorithm::IsClusterShowerLike().

{
    return std::sqrt(LArClusterHelper::GetLengthSquared(pCluster));
}

float lar_content::LArClusterHelper::GetLengthSquared ( const pandora::Cluster *const  pCluster )

static

Get length squared of cluster.

Parameters

pCluster

address of the cluster

Returns

the length squared

Definition at line 60 of file LArClusterHelper.cc.

References max, maxY, min, and minY.

Referenced by lar_content::CrossGapsExtensionAlgorithm::FillClusterAssociationMatrix(), lar_content::DeltaRayMatchingAlgorithm::FindBestParentPfo(), lar_content::OvershootSplittingAlgorithm::GetListOfCleanClusters(), lar_content::LongitudinalExtensionAlgorithm::GetListOfCleanClusters(), lar_content::DeltaRayExtensionAlgorithm::GetListOfCleanClusters(), lar_content::CrossGapsExtensionAlgorithm::GetListOfCleanClusters(), lar_content::CosmicRayExtensionAlgorithm::GetListOfCleanClusters(), lar_content::TransverseExtensionAlgorithm::GetListOfCleanClusters(), lar_content::CosmicRaySplittingAlgorithm::GetListOfCleanClusters(), lar_content::TwoDSlidingFitSplittingAndSwitchingAlgorithm::GetListOfCleanClusters(), lar_content::TwoDSlidingFitSplittingAndSplicingAlgorithm::GetListOfCleanClusters(), lar_content::DeltaRayGrowingAlgorithm::GetListOfSeedClusters(), lar_content::CosmicRayTrackMatchingAlgorithm::SelectCleanClusters(), lar_content::CosmicRayTrackRecoveryAlgorithm::SelectCleanClusters(), lar_content::CandidateVertexCreationAlgorithm::SelectClusters(), lar_content::ThreeDShowersAlgorithm::SelectInputClusters(), lar_content::ThreeDTracksBaseAlgorithm< T >::SelectInputClusters(), lar_content::TwoDSlidingFitConsolidationAlgorithm::SortInputClusters(), lar_content::ParticleRecoveryAlgorithm::StandardClusterSelection(), lar_content::TrackConsolidationAlgorithm::TrackConsolidationAlgorithm(), and lar_content::TwoDSlidingFitSplittingAlgorithm::TwoDSlidingFitSplittingAlgorithm().

{
    const OrderedCaloHitList &orderedCaloHitList(pCluster->GetOrderedCaloHitList());

    if (orderedCaloHitList.empty())
        throw StatusCodeException(STATUS_CODE_NOT_INITIALIZED);

    // ATTN In 2D case, we will actually calculate the quadrature sum of deltaX and deltaU/V/W
    float minX(std::numeric_limits<float>::max()), maxX(-std::numeric_limits<float>::max());
    float minY(std::numeric_limits<float>::max()), maxY(-std::numeric_limits<float>::max());
    float minZ(std::numeric_limits<float>::max()), maxZ(-std::numeric_limits<float>::max());

    for (OrderedCaloHitList::const_iterator iter = orderedCaloHitList.begin(), iterEnd = orderedCaloHitList.end(); iter != iterEnd; ++iter)
    {
        for (CaloHitList::const_iterator hitIter = iter->second->begin(), hitIterEnd = iter->second->end(); hitIter != hitIterEnd; ++hitIter)
        {
            const CartesianVector &hitPosition((*hitIter)->GetPositionVector());
            minX = std::min(hitPosition.GetX(), minX);
            maxX = std::max(hitPosition.GetX(), maxX);
            minY = std::min(hitPosition.GetY(), minY);
            maxY = std::max(hitPosition.GetY(), maxY);
            minZ = std::min(hitPosition.GetZ(), minZ);
            maxZ = std::max(hitPosition.GetZ(), maxZ);
        }
    }

    const float deltaX(maxX - minX), deltaY(maxY - minY), deltaZ(maxZ - minZ);
    return (deltaX * deltaX + deltaY * deltaY + deltaZ * deltaZ);
}

bool lar_content::LArClusterHelper::SortByInnerLayer ( const pandora::Cluster *const  pLhs, const pandora::Cluster *const  pRhs  )

static

Sort clusters by inner layer, then position, then pulse-height.

Parameters

pLhs	address of first cluster
pRhs	address of second cluster

Definition at line 642 of file LArClusterHelper.cc.

Referenced by lar_content::CrossGapsAssociationAlgorithm::GetListOfCleanClusters(), and lar_content::LongitudinalAssociationAlgorithm::GetListOfCleanClusters().

{
    const unsigned int innerLayerLhs(pLhs->GetInnerPseudoLayer());
    const unsigned int innerLayerRhs(pRhs->GetInnerPseudoLayer());

    if (innerLayerLhs != innerLayerRhs)
      return (innerLayerLhs < innerLayerRhs);

    return SortByPosition(pLhs, pRhs);
}

bool lar_content::LArClusterHelper::SortByLayerSpan ( const pandora::Cluster *const  pLhs, const pandora::Cluster *const  pRhs  )

static

Sort clusters by layer span, then inner layer, then position, then pulse-height.

Parameters

pLhs	address of first cluster
pRhs	address of second cluster

Definition at line 629 of file LArClusterHelper.cc.

{
    const unsigned int layerSpanLhs(pLhs->GetOuterPseudoLayer() - pLhs->GetInnerPseudoLayer());
    const unsigned int layerSpanRhs(pRhs->GetOuterPseudoLayer() - pRhs->GetInnerPseudoLayer());

    if (layerSpanLhs != layerSpanRhs)
        return (layerSpanLhs > layerSpanRhs);

    return SortByInnerLayer(pLhs, pRhs);
}

bool lar_content::LArClusterHelper::SortByNHits ( const pandora::Cluster *const  pLhs, const pandora::Cluster *const  pRhs  )

static

Sort clusters by number of hits, then layer span, then inner layer, then position, then pulse-height.

Parameters

pLhs	address of first cluster
pRhs	address of second cluster

Definition at line 616 of file LArClusterHelper.cc.

Referenced by lar_content::TwoDSlidingFitConsolidationAlgorithm::AddHitsToClusters(), lar_content::ClusterAssociationAlgorithm::AmbiguousPropagation(), lar_content::SplitShowersTool::ApplyChanges(), lar_content::ThreeDKinkBaseTool::ApplyChanges(), lar_content::EventSlicingTool::AssignRemainingHitsToSlices(), lar_content::OvershootSplittingAlgorithm::BuildIntersectionMap(), lar_content::TrackHitsBaseTool::BuildSlidingFitMap(), lar_content::OvershootSplittingAlgorithm::BuildSortedIntersectionMap(), lar_content::VertexSelectionBaseAlgorithm::CalculateClusterSlidingFits(), lar_content::ConeClusterMopUpAlgorithm::ClusterMopUp(), lar_content::NearbyClusterMopUpAlgorithm::ClusterMopUp(), lar_content::BoundedClusterMopUpAlgorithm::ClusterMopUp(), lar_pandora::LArPandoraOutput::CollectClusters(), lar_content::EventSlicingTool::CopyPfoHitsToSlices(), lar_content::ParticleRecoveryAlgorithm::ExamineTensor(), lar_content::CrossGapsExtensionAlgorithm::FillClusterMergeMap(), lar_content::DeltaRayExtensionAlgorithm::FillClusterMergeMap(), lar_content::LongitudinalExtensionAlgorithm::FillClusterMergeMap(), lar_content::CosmicRayExtensionAlgorithm::FillClusterMergeMap(), lar_content::TransverseExtensionAlgorithm::FillClusterMergeMap(), lar_content::TransverseAssociationAlgorithm::FillReducedAssociationMap(), lar_content::TransverseAssociationAlgorithm::FillSymmetricAssociationMap(), lar_content::ClearTrackFragmentsTool::GetAffectedKeyClusters(), lar_content::TransverseAssociationAlgorithm::GetAssociatedClusters(), lar_content::CosmicRayTrackRecoveryAlgorithm::GetAvailableClusters(), lar_content::VertexBasedPfoRecoveryAlgorithm::GetAvailableClusters(), lar_content::SlidingConeClusterMopUpAlgorithm::GetAvailableTwoDClusters(), lar_content::DeltaRayMatchingAlgorithm::GetClusters(), lar_content::EventSlicingTool::GetClusterSliceList(), lar_content::ShowerGrowingAlgorithm::GetFigureOfMerit(), lar_content::EventSlicingTool::GetKDTreeEntries3D(), lar_content::OvershootSplittingAlgorithm::GetListOfCleanClusters(), lar_content::CrossGapsExtensionAlgorithm::GetListOfCleanClusters(), lar_content::SimpleClusterGrowingAlgorithm::GetListOfCleanClusters(), lar_content::DeltaRayExtensionAlgorithm::GetListOfCleanClusters(), lar_content::LongitudinalExtensionAlgorithm::GetListOfCleanClusters(), lar_content::SimpleClusterMergingAlgorithm::GetListOfCleanClusters(), lar_content::DeltaRayGrowingAlgorithm::GetListOfCleanClusters(), lar_content::TransverseExtensionAlgorithm::GetListOfCleanClusters(), lar_content::CosmicRayExtensionAlgorithm::GetListOfCleanClusters(), lar_content::CosmicRaySplittingAlgorithm::GetListOfCleanClusters(), lar_content::TwoDSlidingFitSplittingAndSwitchingAlgorithm::GetListOfCleanClusters(), lar_content::TransverseAssociationAlgorithm::GetListOfCleanClusters(), lar_content::TwoDSlidingFitSplittingAndSplicingAlgorithm::GetListOfCleanClusters(), lar_content::ClusterGrowingAlgorithm::GetListOfNonSeedClusters(), lar_content::DeltaRayGrowingAlgorithm::GetListOfSeedClusters(), lar_content::SimpleClusterGrowingAlgorithm::GetListOfSeedClusters(), lar_content::ThreeDTrackFragmentsAlgorithm::GetMatchedClusters(), lar_content::MissingTrackSegmentTool::GetSegmentOverlapMap(), lar_content::SlidingConePfoMopUpAlgorithm::GetThreeDClusters(), lar_content::SlidingConeClusterMopUpAlgorithm::GetThreeDClusters(), lar_content::CrossGapsAssociationAlgorithm::IsExtremalCluster(), lar_content::LongitudinalAssociationAlgorithm::IsExtremalCluster(), lar_content::TransverseAssociationAlgorithm::IsExtremalCluster(), lar_content::ClusterMopUpBaseAlgorithm::MakeClusterMerges(), lar_content::ThreeDBaseAlgorithm< T >::MakeClusterMerges(), lar_content::SlidingConeClusterMopUpAlgorithm::MakeClusterMerges(), lar_content::ThreeDTracksBaseAlgorithm< T >::MakeClusterSplits(), lar_content::MissingTrackSegmentTool::MakeDecisions(), lar_content::SlidingConePfoMopUpAlgorithm::MakePfoMerges(), lar_content::ClusterGrowingAlgorithm::MergeClusters(), lar_content::SlidingConeClusterMopUpAlgorithm::ClusterMerge::operator<(), lar_content::SlidingConePfoMopUpAlgorithm::ClusterMerge::operator<(), lar_content::ThreeDTrackFragmentsAlgorithm::PerformMainLoop(), lar_content::ThreeDBaseAlgorithm< T >::PerformMainLoop(), lar_content::OvershootSplittingAlgorithm::PopulateSplitPositionMap(), lar_content::ShowerGrowingAlgorithm::ProcessSeedAssociationDetails(), lar_content::ClearTrackFragmentsTool::ProcessTensorElement(), lar_content::TwoDSlidingFitSplittingAndSplicingAlgorithm::PruneClusterExtensionList(), lar_content::TwoDSlidingFitConsolidationAlgorithm::RebuildClusters(), lar_content::TwoDSlidingFitConsolidationAlgorithm::RemoveHitsFromClusters(), lar_content::ThreeDBaseAlgorithm< T >::RemoveUnavailableTensorElements(), lar_content::CandidateVertexCreationAlgorithm::SelectClusters(), lar_content::TwoDSlidingFitConsolidationAlgorithm::SortInputClusters(), lar_content::TransverseAssociationAlgorithm::SortInputClusters(), lar_content::TwoDSlidingFitMultiSplitAlgorithm::SplitClusters(), lar_content::ThreeDTrackFragmentsAlgorithm::UpdateForNewCluster(), and lar_content::ThreeDBaseAlgorithm< T >::UpdateForNewCluster().

{
    const unsigned int nHitsLhs(pLhs->GetNCaloHits());
    const unsigned int nHitsRhs(pRhs->GetNCaloHits());

    if (nHitsLhs != nHitsRhs)
        return (nHitsLhs > nHitsRhs);

    return SortByLayerSpan(pLhs, pRhs);
}

bool lar_content::LArClusterHelper::SortByNOccupiedLayers ( const pandora::Cluster *const  pLhs, const pandora::Cluster *const  pRhs  )

static

Sort clusters by number of occupied layers, and by inner layer, then energy in event of a tie.

Parameters

pLhs	address of first cluster
pRhs	address of second cluster

Definition at line 603 of file LArClusterHelper.cc.

{
    const unsigned int nOccupiedLayersLhs(pLhs->GetOrderedCaloHitList().size());
    const unsigned int nOccupiedLayersRhs(pRhs->GetOrderedCaloHitList().size());

    if (nOccupiedLayersLhs != nOccupiedLayersRhs)
        return (nOccupiedLayersLhs > nOccupiedLayersRhs);

    return SortByNHits(pLhs, pRhs);
}

bool lar_content::LArClusterHelper::SortByPosition ( const pandora::Cluster *const  pLhs, const pandora::Cluster *const  pRhs  )

static

Sort clusters by position, then pulse-height.

Parameters

pLhs	address of first cluster
pRhs	address of second cluster

Definition at line 655 of file LArClusterHelper.cc.

{
    const CartesianVector deltaPositionIL(pRhs->GetCentroid(pRhs->GetInnerPseudoLayer()) - pLhs->GetCentroid(pLhs->GetInnerPseudoLayer()));

    if (std::fabs(deltaPositionIL.GetZ()) > std::numeric_limits<float>::epsilon())
        return (deltaPositionIL.GetZ() > std::numeric_limits<float>::epsilon());

    if (std::fabs(deltaPositionIL.GetX()) > std::numeric_limits<float>::epsilon())
        return (deltaPositionIL.GetX() > std::numeric_limits<float>::epsilon());

    if (std::fabs(deltaPositionIL.GetY()) > std::numeric_limits<float>::epsilon())
        return (deltaPositionIL.GetY() > std::numeric_limits<float>::epsilon());

    const CartesianVector deltaPositionOL(pRhs->GetCentroid(pRhs->GetOuterPseudoLayer()) - pLhs->GetCentroid(pLhs->GetOuterPseudoLayer()));

    if (std::fabs(deltaPositionOL.GetZ()) > std::numeric_limits<float>::epsilon())
        return (deltaPositionOL.GetZ() > std::numeric_limits<float>::epsilon());

    if (std::fabs(deltaPositionOL.GetX()) > std::numeric_limits<float>::epsilon())
        return (deltaPositionOL.GetX() > std::numeric_limits<float>::epsilon());

    if (std::fabs(deltaPositionOL.GetY()) > std::numeric_limits<float>::epsilon())
        return (deltaPositionOL.GetY() > std::numeric_limits<float>::epsilon());

    // Use pulse height to resolve ties
    return SortByPulseHeight(pLhs, pRhs);
}

bool lar_content::LArClusterHelper::SortByPulseHeight ( const pandora::Cluster *const  pLhs, const pandora::Cluster *const  pRhs  )

static

Sort clusters by pulse-height.

Parameters

pLhs	address of first cluster
pRhs	address of second cluster

Definition at line 685 of file LArClusterHelper.cc.

{
    return (pLhs->GetHadronicEnergy() > pRhs->GetHadronicEnergy());
}

bool lar_content::LArClusterHelper::SortCoordinatesByPosition ( const pandora::CartesianVector &  lhs, const pandora::CartesianVector &  rhs  )

static

Sort cartesian vectors by their position (use Z, followed by X, followed by Y)

Parameters

lhs	first point
rhs	second point

Definition at line 719 of file LArClusterHelper.cc.

Referenced by lar_content::CandidateVertexCreationAlgorithm::GetSpacepoints().

{
    const CartesianVector deltaPosition(rhs - lhs);

    if (std::fabs(deltaPosition.GetZ()) > std::numeric_limits<float>::epsilon())
        return (deltaPosition.GetZ() > std::numeric_limits<float>::epsilon());

    if (std::fabs(deltaPosition.GetX()) > std::numeric_limits<float>::epsilon())
        return (deltaPosition.GetX() > std::numeric_limits<float>::epsilon());

    return (deltaPosition.GetY() > std::numeric_limits<float>::epsilon());
}

bool lar_content::LArClusterHelper::SortHitsByPosition ( const pandora::CaloHit *const  pLhs, const pandora::CaloHit *const  pRhs  )

static

Sort calo hits by their position (use Z, followed by X, followed by Y)

Parameters

pLhs	address of first calo hit
pRhs	address of second calo hit

Definition at line 692 of file LArClusterHelper.cc.

Referenced by lar_content::TrackClusterCreationAlgorithm::AddFilteredCaloHits(), lar_content::ShowerAsymmetryFeatureTool::CalculateAsymmetryParameters(), lar_content::GlobalAsymmetryFeatureTool::CalculateGlobalAsymmetry(), lar_content::LocalAsymmetryFeatureTool::CalculateLocalAsymmetry(), lar_content::IsolatedClusterMopUpAlgorithm::ClusterMopUp(), lar_pandora::LArPandoraOutput::Collect3DHits(), lar_content::SimpleClusterCreationAlgorithm::CollectAssociatedHits(), lar_content::SimpleClusterCreationAlgorithm::CreateClusters(), lar_content::TrackClusterCreationAlgorithm::CreateClusters(), lar_content::TrackClusterCreationAlgorithm::FilterCaloHits(), lar_content::CrossedTrackSplittingAlgorithm::FindCandidateSplitPositions(), lar_content::ThreeDTrackFragmentsAlgorithm::GetFragmentOverlapResult(), lar_pandora::LArPandoraOutput::GetHitsInCluster(), lar_content::ThreeDTrackFragmentsAlgorithm::GetMatchedHits(), lar_content::SvmVertexSelectionAlgorithm::GetShowerLikeClusterEndPoints(), lar_content::TrackClusterCreationAlgorithm::IdentifyJoins(), lar_content::TrackClusterCreationAlgorithm::MakePrimaryAssociations(), lar_content::TrackClusterCreationAlgorithm::MakeSecondaryAssociations(), and lar_content::ThreeDHitCreationAlgorithm::SeparateTwoDHits().

{
    const CartesianVector deltaPosition(pRhs->GetPositionVector() - pLhs->GetPositionVector());
    
    if (std::fabs(deltaPosition.GetZ()) > std::numeric_limits<float>::epsilon())
        return (deltaPosition.GetZ() > std::numeric_limits<float>::epsilon());

    if (std::fabs(deltaPosition.GetX()) > std::numeric_limits<float>::epsilon())
        return (deltaPosition.GetX() > std::numeric_limits<float>::epsilon());

    if (std::fabs(deltaPosition.GetY()) > std::numeric_limits<float>::epsilon())
        return (deltaPosition.GetY() > std::numeric_limits<float>::epsilon());

    // Use pulse height to resolve ties
    return SortHitsByPulseHeight(pLhs, pRhs);
}

bool lar_content::LArClusterHelper::SortHitsByPulseHeight ( const pandora::CaloHit *const  pLhs, const pandora::CaloHit *const  pRhs  )

static

Sort calo hits by their pulse height.

Parameters

pLhs	address of first calo hit
pRhs	address of second calo hit

Definition at line 711 of file LArClusterHelper.cc.

{
    // TODO: Think about the correct energy to use here (should they ever be different)
    return (pLhs->GetHadronicEnergy() > pRhs->GetHadronicEnergy());
}

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/larpandoracontent/v03_14_03/source/larpandoracontent/LArHelpers/LArClusterHelper.h
• /cvmfs/larsoft.opensciencegrid.org/products/larpandoracontent/v03_14_03/source/larpandoracontent/LArHelpers/LArClusterHelper.cc