lar_content::VertexAssociatedPfosTool Class Reference

VertexAssociatedPfosTool class. More...

#include <VertexAssociatedPfosTool.h>

Inheritance diagram for lar_content::VertexAssociatedPfosTool:

Public Member Functions
	VertexAssociatedPfosTool ()
	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_minVertexLongitudinalDistance
	Vertex association check: min longitudinal distance cut. More...
float	m_maxVertexLongitudinalDistance
	Vertex association check: max longitudinal distance cut. More...
float	m_maxVertexTransverseDistance
	Vertex association check: max transverse distance cut. More...
float	m_vertexAngularAllowance
	Vertex association check: pointing angular allowance in degrees. More...

Detailed Description

VertexAssociatedPfosTool class.

Definition at line 19 of file VertexAssociatedPfosTool.h.

Constructor & Destructor Documentation

lar_content::VertexAssociatedPfosTool::VertexAssociatedPfosTool

(

)

Default constructor.

Definition at line 27 of file VertexAssociatedPfosTool.cc.

                                                   :
    m_minVertexLongitudinalDistance(-2.5f),
    m_maxVertexLongitudinalDistance(20.f),
    m_maxVertexTransverseDistance(3.5f),
    m_vertexAngularAllowance(3.f)
{
}

Member Function Documentation

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

private

Definition at line 74 of file VertexAssociatedPfosTool.cc.

{
    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "MinVertexLongitudinalDistance", m_minVertexLongitudinalDistance));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "MaxVertexLongitudinalDistance", m_maxVertexLongitudinalDistance));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "MaxVertexTransverseDistance", m_maxVertexTransverseDistance));

    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "VertexAngularAllowance", m_vertexAngularAllowance));

    return STATUS_CODE_SUCCESS;
}

void lar_content::VertexAssociatedPfosTool::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 37 of file VertexAssociatedPfosTool.cc.

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

    const CartesianVector &neutrinoVertex(pNeutrinoVertex->GetPosition());

    PfoVector sortedPfos;
    for (const auto &mapEntry : pfoInfoMap)
        sortedPfos.push_back(mapEntry.first);
    std::sort(sortedPfos.begin(), sortedPfos.end(), LArPfoHelper::SortByNHits);

    for (const Pfo *const pPfo : sortedPfos)
    {
        PfoInfo *const pPfoInfo(pfoInfoMap.at(pPfo));

        if (pPfoInfo->IsNeutrinoVertexAssociated() || pPfoInfo->GetParentPfo())
            continue;

        const LArPointingCluster pointingCluster(*(pPfoInfo->GetSlidingFitResult3D()));
        const bool useInner((pointingCluster.GetInnerVertex().GetPosition() - neutrinoVertex).GetMagnitudeSquared() <
                            (pointingCluster.GetOuterVertex().GetPosition() - neutrinoVertex).GetMagnitudeSquared());

        const LArPointingCluster::Vertex &daughterVertex(useInner ? pointingCluster.GetInnerVertex() : pointingCluster.GetOuterVertex());

        if (LArPointingClusterHelper::IsNode(neutrinoVertex, daughterVertex, m_minVertexLongitudinalDistance, m_maxVertexTransverseDistance) ||
            LArPointingClusterHelper::IsEmission(neutrinoVertex, daughterVertex, m_minVertexLongitudinalDistance,
                m_maxVertexLongitudinalDistance, m_maxVertexTransverseDistance, m_vertexAngularAllowance))
        {
            pPfoInfo->SetNeutrinoVertexAssociation(true);
            pPfoInfo->SetInnerLayerAssociation(useInner);
        }
    }
}

Member Data Documentation

float lar_content::VertexAssociatedPfosTool::m_maxVertexLongitudinalDistance

private

Vertex association check: max longitudinal distance cut.

Definition at line 34 of file VertexAssociatedPfosTool.h.

float lar_content::VertexAssociatedPfosTool::m_maxVertexTransverseDistance

private

Vertex association check: max transverse distance cut.

Definition at line 35 of file VertexAssociatedPfosTool.h.

float lar_content::VertexAssociatedPfosTool::m_minVertexLongitudinalDistance

private

Vertex association check: min longitudinal distance cut.

Definition at line 33 of file VertexAssociatedPfosTool.h.

float lar_content::VertexAssociatedPfosTool::m_vertexAngularAllowance

private

Vertex association check: pointing angular allowance in degrees.

Definition at line 36 of file VertexAssociatedPfosTool.h.

---

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

• VertexAssociatedPfosTool.h
• VertexAssociatedPfosTool.cc