lar_content::VertexSplittingAlgorithm Class Reference

VertexSplittingAlgorithm class. More...

#include <VertexSplittingAlgorithm.h>

Inheritance diagram for lar_content::VertexSplittingAlgorithm:

Public Member Functions
	VertexSplittingAlgorithm ()
	Default constructor. More...
Public Member Functions inherited from lar_content::TwoDSlidingFitSplittingAlgorithm
	TwoDSlidingFitSplittingAlgorithm ()
	Default constructor. More...

Private Member Functions
pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)
pandora::StatusCode	FindBestSplitPosition (const TwoDSlidingFitResult &slidingFitResult, pandora::CartesianVector &splitPosition) const
	Use sliding linear fit to identify the best split position. More...

Private Attributes
float	m_splitDisplacementSquared
	Maximum displacement squared. More...
float	m_vertexDisplacementSquared
	Maximum displacement squared. More...

Additional Inherited Members
Protected Member Functions inherited from lar_content::ClusterSplittingAlgorithm
virtual pandora::StatusCode	Run ()
pandora::StatusCode	RunUsingCurrentList () const
	Run the algorithm using the current cluster list as input. More...
Protected Attributes inherited from lar_content::TwoDSlidingFitSplittingAlgorithm
unsigned int	m_slidingFitHalfWindow
float	m_minClusterLength

Detailed Description

VertexSplittingAlgorithm class.

Definition at line 21 of file VertexSplittingAlgorithm.h.

Constructor & Destructor Documentation

lar_content::VertexSplittingAlgorithm::VertexSplittingAlgorithm

(

)

Default constructor.

Definition at line 21 of file VertexSplittingAlgorithm.cc.

                                                   : m_splitDisplacementSquared(4.f * 4.f), m_vertexDisplacementSquared(1.f * 1.f)
{
    // ATTN Some default values differ from base class
    m_minClusterLength = 1.f;
}

Member Function Documentation

StatusCode lar_content::VertexSplittingAlgorithm::FindBestSplitPosition ( const TwoDSlidingFitResult &  slidingFitResult, pandora::CartesianVector &  splitPosition  ) const

privatevirtual

Use sliding linear fit to identify the best split position.

Parameters

slidingFitResult	the input sliding fit result
splitPosition	the best split position

Returns

pandora::StatusCode

Implements lar_content::TwoDSlidingFitSplittingAlgorithm.

Definition at line 29 of file VertexSplittingAlgorithm.cc.

{
    // Identify event vertex
    const VertexList *pVertexList(NULL);
    PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetCurrentList(*this, pVertexList));

    if (pVertexList->empty())
        return STATUS_CODE_NOT_INITIALIZED;

    if (pVertexList->size() != 1)
        return STATUS_CODE_OUT_OF_RANGE;

    const Cluster *const pCluster(slidingFitResult.GetCluster());
    const HitType hitType(LArClusterHelper::GetClusterHitType(pCluster));

    const Vertex *const pSelectedVertex(*(pVertexList->begin()));

    if (VERTEX_3D != pSelectedVertex->GetVertexType())
        return STATUS_CODE_INVALID_PARAMETER;

    const CartesianVector theVertex2D(LArGeometryHelper::ProjectPosition(this->GetPandora(), pSelectedVertex->GetPosition(), hitType));

    const CartesianVector innerVertex2D(slidingFitResult.GetGlobalMinLayerPosition());
    const CartesianVector outerVertex2D(slidingFitResult.GetGlobalMaxLayerPosition());

    if ((outerVertex2D - innerVertex2D).GetMagnitudeSquared() < 4.f * m_vertexDisplacementSquared)
        return STATUS_CODE_NOT_FOUND;

    bool foundSplit(false);
    const StatusCode statusCode(slidingFitResult.GetGlobalFitProjection(theVertex2D, splitPosition));

    if (STATUS_CODE_SUCCESS != statusCode)
        return statusCode;

    const float splitDisplacementSquared((splitPosition - theVertex2D).GetMagnitudeSquared());
    const float vertexDisplacementSquared(
        std::min((splitPosition - innerVertex2D).GetMagnitudeSquared(), (splitPosition - outerVertex2D).GetMagnitudeSquared()));

    if ((splitDisplacementSquared < m_splitDisplacementSquared) && (vertexDisplacementSquared > m_vertexDisplacementSquared) &&
        (splitDisplacementSquared < vertexDisplacementSquared))
    {
        foundSplit = true;
    }

    if (!foundSplit)
        return STATUS_CODE_NOT_FOUND;

    return STATUS_CODE_SUCCESS;
}

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

privatevirtual

Reimplemented from lar_content::TwoDSlidingFitSplittingAlgorithm.

Definition at line 81 of file VertexSplittingAlgorithm.cc.

{
    float splitDisplacement = std::sqrt(m_splitDisplacementSquared);
    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "SplitDisplacement", splitDisplacement));
    m_splitDisplacementSquared = splitDisplacement * splitDisplacement;

    float vertexDisplacement = std::sqrt(m_vertexDisplacementSquared);
    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "VertexDisplacement", vertexDisplacement));
    m_vertexDisplacementSquared = vertexDisplacement * vertexDisplacement;

    return TwoDSlidingFitSplittingAlgorithm::ReadSettings(xmlHandle);
}

Member Data Documentation

float lar_content::VertexSplittingAlgorithm::m_splitDisplacementSquared

private

Maximum displacement squared.

Definition at line 33 of file VertexSplittingAlgorithm.h.

float lar_content::VertexSplittingAlgorithm::m_vertexDisplacementSquared

private

Maximum displacement squared.

Definition at line 34 of file VertexSplittingAlgorithm.h.

---

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

• VertexSplittingAlgorithm.h
• VertexSplittingAlgorithm.cc