lar_content::LArPointingClusterHelper Class Reference

LArPointingClusterHelper class. More...

#include <LArPointingClusterHelper.h>

Static Public Member Functions
static float	GetLengthSquared (const LArPointingCluster &pointingCluster)
	Calculate distance squared between inner and outer vertices of pointing cluster. More...
static float	GetLength (const LArPointingCluster &pointingCluster)
	Calculate distance squared between inner and outer vertices of pointing cluster. More...
static bool	IsNode (const pandora::CartesianVector &parentVertex, const LArPointingCluster::Vertex &daughterVertex, const float minLongitudinalDistance, const float maxTransverseDistance)
	Whether pointing vertex is adjacent to a given position. More...
static bool	IsEmission (const pandora::CartesianVector &parentVertex, const LArPointingCluster::Vertex &daughterVertex, const float minLongitudinalDistance, const float maxLongitudinalDistance, const float maxTransverseDistance, const float angularAllowance)
	Whether pointing vertex is emitted from a given position. More...
static pandora::CartesianVector	GetProjectedPosition (const pandora::CartesianVector &initialPosition, const pandora::CartesianVector &initialDirection, const pandora::Cluster *const pCluster, const float projectionAngularAllowance)
	Get projected position on a cluster from a specified position and direction. More...
static void	GetClosestVertices (const bool useX, const bool useY, const bool useZ, const LArPointingCluster &pointingClusterI, const LArPointingCluster &pointingClusterJ, LArPointingCluster::Vertex &closestVertexI, LArPointingCluster::Vertex &closestVertexJ)
	Given a pair of pointing clusters, receive the closest or farthest pair of vertices. More...
static void	GetClosestVertices (const LArPointingCluster &pointingClusterI, const LArPointingCluster &pointingClusterJ, LArPointingCluster::Vertex &closestVertexI, LArPointingCluster::Vertex &closestVertexJ)
	Given a pair of pointing clusters, find the closest pair of vertices. More...
static void	GetClosestVerticesInX (const LArPointingCluster &pointingClusterI, const LArPointingCluster &pointingClusterJ, LArPointingCluster::Vertex &closestVertexI, LArPointingCluster::Vertex &closestVertexJ)
	Given a pair of pointing clusters, find the pair of vertices with smallest x-separation. More...
static void	GetClosestVerticesInYZ (const LArPointingCluster &pointingClusterI, const LArPointingCluster &pointingClusterJ, LArPointingCluster::Vertex &closestVertexI, LArPointingCluster::Vertex &closestVertexJ)
	Given a pair of pointing clusters, find the pair of vertices with smallest yz-separation. More...
static void	GetImpactParametersInYZ (const LArPointingCluster::Vertex &pointingVertex, const LArPointingCluster::Vertex &targetVertex, float &longitudinal, float &transverse)
	Calculate impact parameters between a pair of pointing vertices using yz-coordinates. More...
static void	GetImpactParameters (const LArPointingCluster::Vertex &pointingVertex, const LArPointingCluster::Vertex &targetVertex, float &longitudinal, float &transverse)
	Calculate impact parameters between a pair of pointing vertices. More...
static void	GetImpactParameters (const LArPointingCluster::Vertex &pointingVertex, const pandora::CartesianVector &targetPosition, float &longitudinal, float &transverse)
	Calculate impact parameters between a pointing cluster vertex and a target position. More...
static void	GetImpactParameters (const pandora::CartesianVector &initialPosition, const pandora::CartesianVector &initialDirection, const pandora::CartesianVector &targetPosition, float &longitudinal, float &transverse)
	Calculate impact parameters of a specified position and direction to a target position. More...
static void	GetIntersection (const LArPointingCluster::Vertex &firstVertex, const LArPointingCluster::Vertex &secondVertex, pandora::CartesianVector &intersectPosition, float &firstDisplacement, float &secondDisplacement)
	Get intersection of two vertices. More...
static void	GetAverageDirection (const LArPointingCluster::Vertex &firstVertex, const LArPointingCluster::Vertex &secondVertex, pandora::CartesianVector &averageDirection)
	Get average direction of two vertices. More...
static void	GetIntersection (const pandora::CartesianVector &firstPosition, const pandora::CartesianVector &firstDirection, const pandora::CartesianVector &secondPosition, const pandora::CartesianVector &secondDirection, pandora::CartesianVector &intersectPosition, float &firstDisplacement, float &secondDisplacement)
	Get intersection of two vertices. More...
static void	GetIntersection (const LArPointingCluster::Vertex &vertexCluster, const pandora::Cluster *const pTargetCluster, pandora::CartesianVector &intersectPosition, float &displacementL, float &displacementT)
	Get intersection of vertex with target cluster. More...
static LArPointingCluster::Vertex	GetBestVertexEstimate (const LArPointingClusterVertexList &vertexList, const LArPointingClusterList &pointingClusterList, const float minLongitudinalDistance, const float maxLongitudinalDistance, const float maxTransverseDistance, const float angularAllowance)
	Simple and fast vertex selection, choosing best vertex from a specified list to represent a set of pointing clusters. More...

Static Private Member Functions
static void	CollectAssociatedClusters (const LArPointingCluster::Vertex &vertex, const LArPointingClusterList &inputList, const float minLongitudinalDistance, const float maxLongitudinalDistance, const float maxTransverseDistance, const float angularAllowance, LArPointingClusterVertexList &outputList)
	Collect cluster vertices, from a provided input list, associated with a specified vertex. More...
static float	GetAssociatedEnergy (const LArPointingCluster::Vertex &vertex, const LArPointingClusterVertexList &clusterVertices)
	Get an estimate of the energy associated with a specified vertex. More...

Detailed Description

LArPointingClusterHelper class.

Definition at line 21 of file LArPointingClusterHelper.h.

Member Function Documentation

void lar_content::LArPointingClusterHelper::CollectAssociatedClusters ( const LArPointingCluster::Vertex &  vertex, const LArPointingClusterList &  inputList, const float  minLongitudinalDistance, const float  maxLongitudinalDistance, const float  maxTransverseDistance, const float  angularAllowance, LArPointingClusterVertexList &  outputList  )

staticprivate

Collect cluster vertices, from a provided input list, associated with a specified vertex.

Parameters

vertex	the vertex
inputList	the input list of pointing clusters
minLongitudinalDistance	the min longitudinal distance cut
maxLongitudinalDistance	the max longitudinal distance cut
maxTransverseDistance	the max transverse distance cut
angularAllowance	the pointing angular allowance in degrees
outputList	to receive the output list of cluster vertices associated with the specified vertex

Definition at line 357 of file LArPointingClusterHelper.cc.

{
    for (LArPointingClusterList::const_iterator iter = inputList.begin(), iterEnd = inputList.end(); iter != iterEnd; ++iter)
    {
        const LArPointingCluster &pointingCluster = *iter;
        const LArPointingCluster::Vertex &innerVertex = pointingCluster.GetInnerVertex();
        const LArPointingCluster::Vertex &outerVertex = pointingCluster.GetOuterVertex();

        const float innerDistanceSquared = (innerVertex.GetPosition() - vertex.GetPosition()).GetMagnitudeSquared();
        const float outerDistanceSquared = (outerVertex.GetPosition() - vertex.GetPosition()).GetMagnitudeSquared();

        const LArPointingCluster::Vertex &chosenVertex((innerDistanceSquared < outerDistanceSquared) ? innerVertex : outerVertex);

        if (LArPointingClusterHelper::IsNode(vertex.GetPosition(), chosenVertex, minLongitudinalDistance, maxTransverseDistance) ||
            LArPointingClusterHelper::IsEmission(vertex.GetPosition(), chosenVertex, minLongitudinalDistance, maxLongitudinalDistance,
                maxTransverseDistance, angularAllowance))
        {
            outputList.push_back(chosenVertex);
        }
    }
}

float lar_content::LArPointingClusterHelper::GetAssociatedEnergy ( const LArPointingCluster::Vertex &  vertex, const LArPointingClusterVertexList &  clusterVertices  )

staticprivate

Get an estimate of the energy associated with a specified vertex.

Parameters

vertex	the vertex
clusterVertices	the list of cluster vertices associated with the specified vertex

Returns

the energy associated with a specified vertex

Definition at line 383 of file LArPointingClusterHelper.cc.

{
    float associatedEnergy(0.f);

    for (LArPointingClusterVertexList::const_iterator iter = associatedVertices.begin(), iterEnd = associatedVertices.end(); iter != iterEnd; ++iter)
    {
        const LArPointingCluster::Vertex &clusterVertex(*iter);
        const Cluster *const pCluster(clusterVertex.GetCluster());

        const float clusterEnergy(LArClusterHelper::GetEnergyFromLength(pCluster));
        const float clusterLength(LArClusterHelper::GetLength(pCluster));
        const float deltaLength(clusterVertex.GetDirection().GetDotProduct(vertex.GetPosition() - clusterVertex.GetPosition()));

        if (deltaLength < std::numeric_limits<float>::epsilon())
        {
            associatedEnergy += clusterEnergy;
        }
        else if (deltaLength < clusterLength)
        {
            associatedEnergy += clusterEnergy * (1.f - (deltaLength / clusterLength));
        }
    }

    return associatedEnergy;
}

void lar_content::LArPointingClusterHelper::GetAverageDirection ( const LArPointingCluster::Vertex &  firstVertex, const LArPointingCluster::Vertex &  secondVertex, pandora::CartesianVector &  averageDirection  )

static

Get average direction of two vertices.

Parameters

firstVertex	the first vertex
secondVertex	the second vertex
averageDirection	the average direction

Definition at line 228 of file LArPointingClusterHelper.cc.

{
    const Cluster *const pFirstCluster(firstVertex.GetCluster());
    const Cluster *const pSecondCluster(secondVertex.GetCluster());

    if (pFirstCluster == pSecondCluster)
        throw pandora::StatusCodeException(STATUS_CODE_INVALID_PARAMETER);

    const float energy1(pFirstCluster->GetHadronicEnergy());
    const float energy2(pSecondCluster->GetHadronicEnergy());

    if ((energy1 < std::numeric_limits<float>::epsilon()) || (energy2 < std::numeric_limits<float>::epsilon()))
        throw pandora::StatusCodeException(STATUS_CODE_NOT_ALLOWED);

    averageDirection = (firstVertex.GetDirection() * energy1 + secondVertex.GetDirection() * energy2).GetUnitVector();
}

LArPointingCluster::Vertex lar_content::LArPointingClusterHelper::GetBestVertexEstimate ( const LArPointingClusterVertexList &  vertexList, const LArPointingClusterList &  pointingClusterList, const float  minLongitudinalDistance, const float  maxLongitudinalDistance, const float  maxTransverseDistance, const float  angularAllowance  )

static

Simple and fast vertex selection, choosing best vertex from a specified list to represent a set of pointing clusters.

Parameters

vertexList	the candidate vertex list
pointingClusterList	the pointing cluster list
minLongitudinalDistance	the min longitudinal distance cut
maxLongitudinalDistance	the max longitudinal distance cut
maxTransverseDistance	the max transverse distance cut
angularAllowance	the pointing angular allowance in degrees

Returns

the best vertex estimate

Definition at line 325 of file LArPointingClusterHelper.cc.

{
    float bestAssociatedEnergy(0.f);
    LArPointingClusterVertexList::const_iterator bestVertexIter(vertexList.end());

    for (LArPointingClusterVertexList::const_iterator iter = vertexList.begin(), iterEnd = vertexList.end(); iter != iterEnd; ++iter)
    {
        const LArPointingCluster::Vertex &vertex(*iter);

        LArPointingClusterVertexList associatedVertices;
        LArPointingClusterHelper::CollectAssociatedClusters(vertex, pointingClusterList, minLongitudinalDistance, maxLongitudinalDistance,
            maxTransverseDistance, angularAllowance, associatedVertices);

        const float associatedEnergy(LArPointingClusterHelper::GetAssociatedEnergy(vertex, associatedVertices));

        if (associatedEnergy > bestAssociatedEnergy)
        {
            bestVertexIter = iter;
            bestAssociatedEnergy = associatedEnergy;
        }
    }

    if (vertexList.end() == bestVertexIter)
        throw StatusCodeException(STATUS_CODE_NOT_FOUND);

    return (*bestVertexIter);
}

void lar_content::LArPointingClusterHelper::GetClosestVertices ( const bool  useX, const bool  useY, const bool  useZ, const LArPointingCluster &  pointingClusterI, const LArPointingCluster &  pointingClusterJ, LArPointingCluster::Vertex &  closestVertexI, LArPointingCluster::Vertex &  closestVertexJ  )

static

Given a pair of pointing clusters, receive the closest or farthest pair of vertices.

Parameters

useX	use relative X coordinates in the calculation
useY	use relative Y coordinates in the calculation
useZ	use relative Z coordinates in the calculation
pointingClusterI	the first pointing cluster
pointingClusterJ	the second pointing cluster
closestVertexI	to receive the relevant vertex from the first pointing cluster
closestVertexJ	to receive the relevant vertex from the second pointing cluster

Definition at line 107 of file LArPointingClusterHelper.cc.

{
    if (pointingClusterI.GetCluster() == pointingClusterJ.GetCluster())
        throw StatusCodeException(STATUS_CODE_NOT_FOUND);

    if (!useX && !useY && !useZ)
        throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);

    for (unsigned int useInnerI = 0; useInnerI < 2; ++useInnerI)
    {
        const LArPointingCluster::Vertex &vtxI(useInnerI == 1 ? pointingClusterI.GetInnerVertex() : pointingClusterI.GetOuterVertex());
        const LArPointingCluster::Vertex &endI(useInnerI == 0 ? pointingClusterI.GetInnerVertex() : pointingClusterI.GetOuterVertex());

        for (unsigned int useInnerJ = 0; useInnerJ < 2; ++useInnerJ)
        {
            const LArPointingCluster::Vertex &vtxJ(useInnerJ == 1 ? pointingClusterJ.GetInnerVertex() : pointingClusterJ.GetOuterVertex());
            const LArPointingCluster::Vertex &endJ(useInnerJ == 0 ? pointingClusterJ.GetInnerVertex() : pointingClusterJ.GetOuterVertex());

            const float vtxI_vtxJ_dx(useX ? (vtxI.GetPosition().GetX() - vtxJ.GetPosition().GetX()) : 0.f);
            const float vtxI_vtxJ_dy(useY ? (vtxI.GetPosition().GetY() - vtxJ.GetPosition().GetY()) : 0.f);
            const float vtxI_vtxJ_dz(useZ ? (vtxI.GetPosition().GetZ() - vtxJ.GetPosition().GetZ()) : 0.f);
            const float vtxI_vtxJ(vtxI_vtxJ_dx * vtxI_vtxJ_dx + vtxI_vtxJ_dy * vtxI_vtxJ_dy + vtxI_vtxJ_dz * vtxI_vtxJ_dz);

            const float vtxI_endJ_dx(useX ? (vtxI.GetPosition().GetX() - endJ.GetPosition().GetX()) : 0.f);
            const float vtxI_endJ_dy(useY ? (vtxI.GetPosition().GetY() - endJ.GetPosition().GetY()) : 0.f);
            const float vtxI_endJ_dz(useZ ? (vtxI.GetPosition().GetZ() - endJ.GetPosition().GetZ()) : 0.f);
            const float vtxI_endJ(vtxI_endJ_dx * vtxI_endJ_dx + vtxI_endJ_dy * vtxI_endJ_dy + vtxI_endJ_dz * vtxI_endJ_dz);

            const float endI_vtxJ_dx(useX ? (endI.GetPosition().GetX() - vtxJ.GetPosition().GetX()) : 0.f);
            const float endI_vtxJ_dy(useY ? (endI.GetPosition().GetY() - vtxJ.GetPosition().GetY()) : 0.f);
            const float endI_vtxJ_dz(useZ ? (endI.GetPosition().GetZ() - vtxJ.GetPosition().GetZ()) : 0.f);
            const float endI_vtxJ(endI_vtxJ_dx * endI_vtxJ_dx + endI_vtxJ_dy * endI_vtxJ_dy + endI_vtxJ_dz * endI_vtxJ_dz);

            const float endI_endJ_dx(useX ? (endI.GetPosition().GetX() - endJ.GetPosition().GetX()) : 0.f);
            const float endI_endJ_dy(useY ? (endI.GetPosition().GetY() - endJ.GetPosition().GetY()) : 0.f);
            const float endI_endJ_dz(useZ ? (endI.GetPosition().GetZ() - endJ.GetPosition().GetZ()) : 0.f);
            const float endI_endJ(endI_endJ_dx * endI_endJ_dx + endI_endJ_dy * endI_endJ_dy + endI_endJ_dz * endI_endJ_dz);

            if ((vtxI_vtxJ < std::min(vtxI_endJ, std::min(endI_vtxJ, endI_endJ))) && (endI_endJ > std::max(vtxI_endJ, std::max(endI_vtxJ, vtxI_vtxJ))))
            {
                closestVertexI = vtxI;
                closestVertexJ = vtxJ;
                return;
            }
        }
    }

    throw StatusCodeException(STATUS_CODE_NOT_FOUND);
}

void lar_content::LArPointingClusterHelper::GetClosestVertices ( const LArPointingCluster &  pointingClusterI, const LArPointingCluster &  pointingClusterJ, LArPointingCluster::Vertex &  closestVertexI, LArPointingCluster::Vertex &  closestVertexJ  )

static

Given a pair of pointing clusters, find the closest pair of vertices.

Parameters

pointingClusterI	the first pointing cluster
pointingClusterJ	the second pointing cluster
closestVertexI	to receive the relevant vertex from the first pointing cluster
closestVertexJ	to receive the relevant vertex from the second pointing cluster

Definition at line 160 of file LArPointingClusterHelper.cc.

{
    return LArPointingClusterHelper::GetClosestVertices(true, true, true, pointingClusterI, pointingClusterJ, closestVertexI, closestVertexJ);
}

void lar_content::LArPointingClusterHelper::GetClosestVerticesInX ( const LArPointingCluster &  pointingClusterI, const LArPointingCluster &  pointingClusterJ, LArPointingCluster::Vertex &  closestVertexI, LArPointingCluster::Vertex &  closestVertexJ  )

static

Given a pair of pointing clusters, find the pair of vertices with smallest x-separation.

Parameters

pointingClusterI	the first pointing cluster
pointingClusterJ	the second pointing cluster
closestVertexI	to receive the relevant vertex from the first pointing cluster
closestVertexJ	to receive the relevant vertex from the second pointing cluster

Definition at line 168 of file LArPointingClusterHelper.cc.

{
    return LArPointingClusterHelper::GetClosestVertices(true, false, false, pointingClusterI, pointingClusterJ, closestVertexI, closestVertexJ);
}

void lar_content::LArPointingClusterHelper::GetClosestVerticesInYZ ( const LArPointingCluster &  pointingClusterI, const LArPointingCluster &  pointingClusterJ, LArPointingCluster::Vertex &  closestVertexI, LArPointingCluster::Vertex &  closestVertexJ  )

static

Given a pair of pointing clusters, find the pair of vertices with smallest yz-separation.

Parameters

pointingClusterI	the first pointing cluster
pointingClusterJ	the second pointing cluster
closestVertexI	to receive the relevant vertex from the first pointing cluster
closestVertexJ	to receive the relevant vertex from the second pointing cluster

Definition at line 176 of file LArPointingClusterHelper.cc.

{
    return LArPointingClusterHelper::GetClosestVertices(false, true, true, pointingClusterI, pointingClusterJ, closestVertexI, closestVertexJ);
}

void lar_content::LArPointingClusterHelper::GetImpactParameters ( const LArPointingCluster::Vertex &  pointingVertex, const LArPointingCluster::Vertex &  targetVertex, float &  longitudinal, float &  transverse  )

static

Calculate impact parameters between a pair of pointing vertices.

Parameters

pointingVertex	the pointing vertex
targetPosition	the target position
longitidinal	the longitudinal displacement
transverse	the transverse displacement

Definition at line 199 of file LArPointingClusterHelper.cc.

{
    return LArPointingClusterHelper::GetImpactParameters(
        pointingVertex.GetPosition(), pointingVertex.GetDirection(), targetVertex.GetPosition(), longitudinal, transverse);
}

static void lar_content::LArPointingClusterHelper::GetImpactParameters ( const LArPointingCluster::Vertex &  pointingVertex, const pandora::CartesianVector &  targetPosition, float &  longitudinal, float &  transverse  )

static

Calculate impact parameters between a pointing cluster vertex and a target position.

Parameters

pointingVertex	the pointing vertex
targetPosition	the target position
longitidinal	the longitudinal displacement
transverse	the transverse displacement

static void lar_content::LArPointingClusterHelper::GetImpactParameters ( const pandora::CartesianVector &  initialPosition, const pandora::CartesianVector &  initialDirection, const pandora::CartesianVector &  targetPosition, float &  longitudinal, float &  transverse  )

static

Calculate impact parameters of a specified position and direction to a target position.

Parameters

initialPosition	the initial position of the cluster
initialDirection	the initial direction of the cluster
targetPosition	the target position
longitidinal	the longitudinal displacement
transverse	the transverse displacement

void lar_content::LArPointingClusterHelper::GetImpactParametersInYZ ( const LArPointingCluster::Vertex &  pointingVertex, const LArPointingCluster::Vertex &  targetVertex, float &  longitudinal, float &  transverse  )

static

Calculate impact parameters between a pair of pointing vertices using yz-coordinates.

Parameters

pointingVertex	the pointing vertex
targetPosition	the target position
longitidinal	the longitudinal displacement
transverse	the transverse displacement

Definition at line 184 of file LArPointingClusterHelper.cc.

{
    if (std::fabs(initialVertex.GetDirection().GetX()) > 1.f - std::numeric_limits<float>::epsilon())
        throw pandora::StatusCodeException(pandora::STATUS_CODE_NOT_FOUND);

    const pandora::CartesianVector initialPosition(0.f, initialVertex.GetPosition().GetY(), initialVertex.GetPosition().GetZ());
    const pandora::CartesianVector initialDirection(0.f, initialVertex.GetDirection().GetY(), initialVertex.GetDirection().GetZ());
    const pandora::CartesianVector targetPosition(0.f, targetVertex.GetPosition().GetY(), targetVertex.GetPosition().GetZ());

    return LArPointingClusterHelper::GetImpactParameters(initialPosition, initialDirection.GetUnitVector(), targetPosition, longitudinal, transverse);
}

static void lar_content::LArPointingClusterHelper::GetIntersection ( const LArPointingCluster::Vertex &  firstVertex, const LArPointingCluster::Vertex &  secondVertex, pandora::CartesianVector &  intersectPosition, float &  firstDisplacement, float &  secondDisplacement  )

static

Get intersection of two vertices.

Parameters

firstVertex	the first vertex
secondVertex	the second vertex
intersectPosition	the position of the intersection
firstDisplacement	the displacement from first vertex to intersection
secondDisplacement	the displacement from second vertex to intersection

static void lar_content::LArPointingClusterHelper::GetIntersection ( const pandora::CartesianVector &  firstPosition, const pandora::CartesianVector &  firstDirection, const pandora::CartesianVector &  secondPosition, const pandora::CartesianVector &  secondDirection, pandora::CartesianVector &  intersectPosition, float &  firstDisplacement, float &  secondDisplacement  )

static

Get intersection of two vertices.

Parameters

firstPosition	the position of the first vertex
firstDirection	the direction of the first vertex
secondPosition	the position of the second vertex
secondDirection	the direction of the second vertex
intersectPosition	the position of the intersection
firstDisplacement	the displacement from first vertex to intersection
secondDisplacement	the displacement from second vertex to intersection

static void lar_content::LArPointingClusterHelper::GetIntersection ( const LArPointingCluster::Vertex &  vertexCluster, const pandora::Cluster *const  pTargetCluster, pandora::CartesianVector &  intersectPosition, float &  displacementL, float &  displacementT  )

static

Get intersection of vertex with target cluster.

Parameters

vertexCluster	the vertex
pTargetCluster	the target cluster
intersectPosition	the position of the intersection
displacementL	the longitudinal displacement
displacementT	the transverse displacement

float lar_content::LArPointingClusterHelper::GetLength ( const LArPointingCluster &  pointingCluster )

static

Calculate distance squared between inner and outer vertices of pointing cluster.

Parameters

pointingCluster

the input pointing cluster

Returns

float the distance

Definition at line 26 of file LArPointingClusterHelper.cc.

{
    return std::sqrt(LArPointingClusterHelper::GetLengthSquared(pointingCluster));
}

float lar_content::LArPointingClusterHelper::GetLengthSquared ( const LArPointingCluster &  pointingCluster )

static

Calculate distance squared between inner and outer vertices of pointing cluster.

Parameters

pointingCluster

the input pointing cluster

Returns

float the distance squared

Definition at line 17 of file LArPointingClusterHelper.cc.

{
    const LArPointingCluster::Vertex &innerVertex(pointingCluster.GetInnerVertex());
    const LArPointingCluster::Vertex &outerVertex(pointingCluster.GetOuterVertex());
    return (innerVertex.GetPosition() - outerVertex.GetPosition()).GetMagnitudeSquared();
}

CartesianVector lar_content::LArPointingClusterHelper::GetProjectedPosition ( const pandora::CartesianVector &  initialPosition, const pandora::CartesianVector &  initialDirection, const pandora::Cluster *const  pCluster, const float  projectionAngularAllowance  )

static

Get projected position on a cluster from a specified position and direction.

Parameters

initialPosition	the initial position of the cluster
initialDirection	the initial direction of the cluster
pCluster	address of the cluster
projectionAngularAllowance	the projection angular allowance

Returns

the projected position

Definition at line 66 of file LArPointingClusterHelper.cc.

{
    const CaloHit *pClosestCaloHit(nullptr);
    float closestDistanceSquared(std::numeric_limits<float>::max());
    const float minCosTheta(std::cos(M_PI * projectionAngularAllowance / 180.f));

    for (const OrderedCaloHitList::value_type &layerEntry : pCluster->GetOrderedCaloHitList())
    {
        for (const CaloHit *const pCaloHit : *layerEntry.second)
        {
            const CartesianVector hitProjection(pCaloHit->GetPositionVector() - vertexPosition);
            const float distanceSquared(hitProjection.GetMagnitudeSquared());

            if (distanceSquared > std::numeric_limits<float>::epsilon())
            {
                // TODO Try to give more weight to on-axis projections
                if (distanceSquared < closestDistanceSquared)
                {
                    if (-hitProjection.GetUnitVector().GetDotProduct(vertexDirection) > minCosTheta)
                    {
                        pClosestCaloHit = pCaloHit;
                        closestDistanceSquared = distanceSquared;
                    }
                }
            }
            else
            {
                return pCaloHit->GetPositionVector();
            }
        }
    }

    if (pClosestCaloHit)
        return pClosestCaloHit->GetPositionVector();

    throw StatusCodeException(STATUS_CODE_NOT_FOUND);
}

bool lar_content::LArPointingClusterHelper::IsEmission ( const pandora::CartesianVector &  parentVertex, const LArPointingCluster::Vertex &  daughterVertex, const float  minLongitudinalDistance, const float  maxLongitudinalDistance, const float  maxTransverseDistance, const float  angularAllowance  )

static

Whether pointing vertex is emitted from a given position.

Parameters

parentVertex	the parent vertex position
daughtervertex	the daughter pointing vertex
minLongitudinalDistance	the min longitudinal distance cut
maxLongitudinalDistance	the max longitudinal distance cut
maxTransverseDistance	the max transverse distance cut
angularAllowance	the pointing angular allowance in degrees

Returns

boolean

Definition at line 47 of file LArPointingClusterHelper.cc.

{
    float rL(0.f), rT(0.f);
    LArPointingClusterHelper::GetImpactParameters(daughterVertex.GetPosition(), daughterVertex.GetDirection(), parentVertex, rL, rT);

    if (std::fabs(rL) > std::fabs(minLongitudinalDistance) && (rL < 0 || rL > maxLongitudinalDistance))
        return false;

    const float tanSqTheta(std::pow(std::tan(M_PI * angularAllowance / 180.f), 2.0));

    if (rT * rT > maxTransverseDistance * maxTransverseDistance + rL * rL * tanSqTheta)
        return false;

    return true;
}

bool lar_content::LArPointingClusterHelper::IsNode ( const pandora::CartesianVector &  parentVertex, const LArPointingCluster::Vertex &  daughterVertex, const float  minLongitudinalDistance, const float  maxTransverseDistance  )

static

Whether pointing vertex is adjacent to a given position.

Parameters

parentVertex	the parent vertex position
daughtervertex	the daughter pointing vertex
minLongitudinalDistance	the min longitudinal distance cut
maxTransverseDistance	the max transverse distance cut

Returns

boolean

Definition at line 33 of file LArPointingClusterHelper.cc.

{
    float rL(0.f), rT(0.f);
    LArPointingClusterHelper::GetImpactParameters(daughterVertex.GetPosition(), daughterVertex.GetDirection(), parentVertex, rL, rT);

    if (std::fabs(rL) > std::fabs(minLongitudinalDistance) || rT > maxTransverseDistance)
        return false;

    return true;
}

---

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

• larpandoracontent/larpandoracontent/LArHelpers/LArPointingClusterHelper.h
• larpandoracontent/larpandoracontent/LArHelpers/LArPointingClusterHelper.cc