lar_content::BranchAssociatedPfosTool Class Reference

BranchAssociatedPfosTool class. More...

#include "BranchAssociatedPfosTool.h"

Inheritance diagram for lar_content::BranchAssociatedPfosTool:

Public Member Functions
	BranchAssociatedPfosTool ()
	Default constructor. More...
void	Run (const NeutrinoHierarchyAlgorithm *const pAlgorithm, const pandora::Vertex *const pNeutrinoVertex, NeutrinoHierarchyAlgorithm::PfoInfoMap &pfoInfoMap)
	Run the algorithm tool. More...

Private Member Functions
pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)

Private Attributes
float	m_minNeutrinoVertexDistance
	Branch association: min distance from branch vertex to neutrino vertex. More...
float	m_trackBranchAdditionFraction
	Branch association: min fraction of length along parent track before association allowed. More...
float	m_maxParentClusterDistance
	Branch association: max distance from branch vertex to a hit in parent 3D cluster. More...

Detailed Description

BranchAssociatedPfosTool class.

Definition at line 19 of file BranchAssociatedPfosTool.h.

Constructor & Destructor Documentation

lar_content::BranchAssociatedPfosTool::BranchAssociatedPfosTool

(

)

Default constructor.

Definition at line 27 of file BranchAssociatedPfosTool.cc.

                                                   :
    m_minNeutrinoVertexDistance(5.f),
    m_trackBranchAdditionFraction(0.4f),
    m_maxParentClusterDistance(3.5f)
{
}

Member Function Documentation

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

private

Definition at line 105 of file BranchAssociatedPfosTool.cc.

References m_maxParentClusterDistance, m_minNeutrinoVertexDistance, and m_trackBranchAdditionFraction.

{
    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "MinNeutrinoVertexDistance", m_minNeutrinoVertexDistance));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "TrackBranchAdditionFraction", m_trackBranchAdditionFraction));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "MaxParentClusterDistance", m_maxParentClusterDistance));

    return STATUS_CODE_SUCCESS;
}

void lar_content::BranchAssociatedPfosTool::Run ( const NeutrinoHierarchyAlgorithm *const  pAlgorithm, const pandora::Vertex *const  pNeutrinoVertex, NeutrinoHierarchyAlgorithm::PfoInfoMap &  pfoInfoMap  )

virtual

Run the algorithm tool.

Parameters

pAlgorithm	address of the calling algorithm
pNeutrinoVertex	the address of the three dimensional neutrino interaction vertex
pfoInfoMap	mapping from pfos to three dimensional clusters, sliding fits, vertices, etc.

Implements lar_content::PfoRelationTool.

Definition at line 36 of file BranchAssociatedPfosTool.cc.

References lar_content::LArClusterHelper::GetClosestDistance(), lar_content::ThreeDSlidingFitResult::GetGlobalMaxLayerPosition(), lar_content::ThreeDSlidingFitResult::GetGlobalMinLayerPosition(), lar_content::LArPfoHelper::IsTrack(), m_maxParentClusterDistance, m_minNeutrinoVertexDistance, m_trackBranchAdditionFraction, and lar_content::NeutrinoHierarchyAlgorithm::SeparatePfos().

{
    if (PandoraContentApi::GetSettings(*pAlgorithm)->ShouldDisplayAlgorithmInfo())
        std::cout << "----> Running Algorithm Tool: " << this->GetInstanceName() << ", " << this->GetType() << std::endl;

    bool associationsMade(true);

    while (associationsMade)
    {
        associationsMade = false;
        PfoVector assignedPfos, unassignedPfos;
        pAlgorithm->SeparatePfos(pfoInfoMap, assignedPfos, unassignedPfos);

        if (unassignedPfos.empty())
            break;

        // ATTN May want to reconsider precise association mechanics for complex situations
        PfoSet recentlyAssigned;

        for (const ParticleFlowObject *const pParentPfo : assignedPfos)
        {
            PfoInfo *const pParentPfoInfo(pfoInfoMap.at(pParentPfo));
            const Cluster *const pParentCluster3D(pParentPfoInfo->GetCluster3D());

            const bool parentIsTrack(LArPfoHelper::IsTrack(pParentPfo));
            const ThreeDSlidingFitResult &parentFitResult(*(pParentPfoInfo->GetSlidingFitResult3D()));

            const float parentLength3D((parentFitResult.GetGlobalMinLayerPosition() - parentFitResult.GetGlobalMaxLayerPosition()).GetMagnitude());
            const CartesianVector &parentVertexPosition(pParentPfoInfo->IsInnerLayerAssociated() ? parentFitResult.GetGlobalMinLayerPosition()
                                                                                                 : parentFitResult.GetGlobalMaxLayerPosition());

            for (const ParticleFlowObject *const pPfo : unassignedPfos)
            {
                if (recentlyAssigned.count(pPfo))
                    continue;

                PfoInfo *const pPfoInfo(pfoInfoMap.at(pPfo));
                const LArPointingCluster pointingCluster(*(pPfoInfo->GetSlidingFitResult3D()));

                const float dNeutrinoVertex(std::min((pointingCluster.GetInnerVertex().GetPosition() - pNeutrinoVertex->GetPosition()).GetMagnitude(),
                    (pointingCluster.GetOuterVertex().GetPosition() - pNeutrinoVertex->GetPosition()).GetMagnitude()));

                if (dNeutrinoVertex < m_minNeutrinoVertexDistance)
                    continue;

                const float dParentVertex(std::min((pointingCluster.GetInnerVertex().GetPosition() - parentVertexPosition).GetMagnitude(),
                    (pointingCluster.GetOuterVertex().GetPosition() - parentVertexPosition).GetMagnitude()));

                if (parentIsTrack && (dParentVertex < m_trackBranchAdditionFraction * parentLength3D))
                    continue;

                const float dInnerVertex(LArClusterHelper::GetClosestDistance(pointingCluster.GetInnerVertex().GetPosition(), pParentCluster3D));
                const float dOuterVertex(LArClusterHelper::GetClosestDistance(pointingCluster.GetOuterVertex().GetPosition(), pParentCluster3D));

                if ((dInnerVertex < m_maxParentClusterDistance) || (dOuterVertex < m_maxParentClusterDistance))
                {
                    associationsMade = true;
                    pParentPfoInfo->AddDaughterPfo(pPfoInfo->GetThisPfo());
                    pPfoInfo->SetParentPfo(pParentPfoInfo->GetThisPfo());
                    pPfoInfo->SetInnerLayerAssociation(dInnerVertex < dOuterVertex);
                    recentlyAssigned.insert(pPfoInfo->GetThisPfo());
                }
            }
        }
    }
}

Member Data Documentation

float lar_content::BranchAssociatedPfosTool::m_maxParentClusterDistance

private

Branch association: max distance from branch vertex to a hit in parent 3D cluster.

Definition at line 35 of file BranchAssociatedPfosTool.h.

Referenced by ReadSettings(), and Run().

float lar_content::BranchAssociatedPfosTool::m_minNeutrinoVertexDistance

private

Branch association: min distance from branch vertex to neutrino vertex.

Definition at line 33 of file BranchAssociatedPfosTool.h.

Referenced by ReadSettings(), and Run().

float lar_content::BranchAssociatedPfosTool::m_trackBranchAdditionFraction

private

Branch association: min fraction of length along parent track before association allowed.

Definition at line 34 of file BranchAssociatedPfosTool.h.

Referenced by ReadSettings(), and Run().

---

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

• larpandoracontent/v04_08_01/source/larpandoracontent/LArThreeDReco/LArEventBuilding/BranchAssociatedPfosTool.h
• larpandoracontent/v04_08_01/source/larpandoracontent/LArThreeDReco/LArEventBuilding/BranchAssociatedPfosTool.cc