lar_content::TransverseMatrixVisualizationTool Class Reference

TransverseMatrixVisualizationTool class. More...

#include "TransverseMatrixVisualizationTool.h"

Inheritance diagram for lar_content::TransverseMatrixVisualizationTool:

Public Types
typedef TwoViewTransverseTracksAlgorithm::MatchingType::MatrixType	MatrixType
typedef std::vector< MatrixType::ElementList::const_iterator >	IteratorList

Public Member Functions
	TransverseMatrixVisualizationTool ()
	Default constructor. More...
bool	Run (TwoViewTransverseTracksAlgorithm *const pAlgorithm, MatrixType &overlapMatrix)
	Run the algorithm tool. More...

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

Private Attributes
unsigned int	m_minClusterConnections
	The minimum number of cluster connections for display. More...
bool	m_ignoreUnavailableClusters
	Whether to ignore (skip-over) unavailable clusters in the matrix. More...
bool	m_showEachIndividualElement
	Whether to draw each individual matrix element. More...
bool	m_showOnlyTrueMatchIndividualElements
	Whether to draw only truly matching individual matrix elements. More...

Detailed Description

TransverseMatrixVisualizationTool class.

Definition at line 19 of file TransverseMatrixVisualizationTool.h.

Member Typedef Documentation

typedef std::vector<MatrixType::ElementList::const_iterator> lar_content::TransverseMatrixTool::IteratorList

inherited

Definition at line 103 of file TwoViewTransverseTracksAlgorithm.h.

typedef TwoViewTransverseTracksAlgorithm::MatchingType::MatrixType lar_content::TransverseMatrixTool::MatrixType

inherited

Definition at line 102 of file TwoViewTransverseTracksAlgorithm.h.

Constructor & Destructor Documentation

lar_content::TransverseMatrixVisualizationTool::TransverseMatrixVisualizationTool

(

)

Default constructor.

Definition at line 18 of file TransverseMatrixVisualizationTool.cc.

                                                                     :
    m_minClusterConnections(1),
    m_ignoreUnavailableClusters(true),
    m_showEachIndividualElement(false),
    m_showOnlyTrueMatchIndividualElements(false)
{
}

Member Function Documentation

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

private

Definition at line 137 of file TransverseMatrixVisualizationTool.cc.

References m_ignoreUnavailableClusters, m_minClusterConnections, m_showEachIndividualElement, and m_showOnlyTrueMatchIndividualElements.

{
    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MinClusterConnections", m_minClusterConnections));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "IgnoreUnavailableClusters", m_ignoreUnavailableClusters));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "ShowEachIndividualElement", m_showEachIndividualElement));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "ShowOnlyTrueMatchIndividualElements", m_showOnlyTrueMatchIndividualElements));

    return STATUS_CODE_SUCCESS;
}

bool lar_content::TransverseMatrixVisualizationTool::Run ( TwoViewTransverseTracksAlgorithm *const  pAlgorithm, MatrixType &  overlapMatrix  )

virtual

Run the algorithm tool.

Parameters

pAlgorithm	address of the calling algorithm
overlapMatrix	the overlap matrix

Returns

whether changes have been made by the tool

Implements lar_content::TransverseMatrixTool.

Definition at line 28 of file TransverseMatrixVisualizationTool.cc.

References util::counter(), f, lar_content::OverlapMatrix< T >::GetConnectedElements(), lar_content::OverlapMatrix< T >::GetSortedKeyClusters(), m_ignoreUnavailableClusters, m_minClusterConnections, m_showEachIndividualElement, and m_showOnlyTrueMatchIndividualElements.

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

    ClusterSet usedKeyClusters;
    ClusterVector sortedKeyClusters;
    overlapMatrix.GetSortedKeyClusters(sortedKeyClusters);

    for (const Cluster *const pKeyCluster : sortedKeyClusters)
    {
        if (m_ignoreUnavailableClusters && !pKeyCluster->IsAvailable())
            continue;

        if (usedKeyClusters.count(pKeyCluster))
            continue;

        unsigned int n1(0), n2(0);
        MatrixType::ElementList elementList;
        overlapMatrix.GetConnectedElements(pKeyCluster, m_ignoreUnavailableClusters, elementList, n1, n2);

        if ((n1 < m_minClusterConnections) && (n2 < m_minClusterConnections))
            continue;

        if (n1 * n2 == 0)
            continue;

        int counter(0);
        ClusterList allClusterList1, allClusterList2;
        std::cout << " Connections: n1 " << n1 << ", n2 " << n2 << ", nElements " << elementList.size() << std::endl;

        for (MatrixType::ElementList::const_iterator eIter = elementList.begin(); eIter != elementList.end(); ++eIter)
        {
            if (allClusterList1.end() == std::find(allClusterList1.begin(), allClusterList1.end(), eIter->GetCluster1()))
                allClusterList1.push_back(eIter->GetCluster1());
            if (allClusterList2.end() == std::find(allClusterList2.begin(), allClusterList2.end(), eIter->GetCluster2()))
                allClusterList2.push_back(eIter->GetCluster2());
            usedKeyClusters.insert(eIter->GetCluster1());
        }

        for (MatrixType::ElementList::const_iterator eIter = elementList.begin(); eIter != elementList.end(); ++eIter)
        {
            int pdg0(0);
            int pdg1(0);
            bool isPrimary0(false);
            bool isPrimary1(false);
            bool sameParticle(false);
            try
            {
                const MCParticle *particle0(MCParticleHelper::GetMainMCParticle(eIter->GetCluster1()));
                const MCParticle *particle1(MCParticleHelper::GetMainMCParticle(eIter->GetCluster2()));
                pdg0 = (particle0->GetParticleId());
                isPrimary0 = (particle0->IsRootParticle());
                pdg1 = (particle1->GetParticleId());
                isPrimary1 = (particle1->IsRootParticle());
                sameParticle = (particle0->GetUid() == particle1->GetUid());
            }
            catch (const StatusCodeException &)
            {
            };

            if (m_showOnlyTrueMatchIndividualElements && !sameParticle)
                continue;

            std::cout << " Element " << counter++ << std::endl;
            std::cout << " ---True PDG 0: " << pdg0 << std::endl;
            std::cout << " ---True PDG 1: " << pdg1 << std::endl;
            std::cout << " ---True is primary 0: " << isPrimary0 << std::endl;
            std::cout << " ---True is primary 1: " << isPrimary1 << std::endl;
            std::cout << " ---True is same particle: " << sameParticle << std::endl;
            std::cout << " ---Is cluster 0 available: " << eIter->GetCluster1()->IsAvailable() << std::endl;
            std::cout << " ---Is cluster 1 available: " << eIter->GetCluster2()->IsAvailable() << std::endl;
            std::cout << " ---XOverlap: " << eIter->GetOverlapResult().GetTwoViewXOverlap().GetTwoViewXOverlapSpan() << std::endl;
            std::cout << " ---XOverlap fraction view0: " << eIter->GetOverlapResult().GetTwoViewXOverlap().GetXOverlapFraction0() << std::endl;
            std::cout << " ---XOverlap fraction view1: " << eIter->GetOverlapResult().GetTwoViewXOverlap().GetXOverlapFraction1() << std::endl;
            std::cout << " ---Matching score: " << eIter->GetOverlapResult().GetMatchingScore() << std::endl;
            std::cout << " ---N. sampling points: " << eIter->GetOverlapResult().GetNSamplingPoints() << std::endl;
            std::cout << " ---N. matched sampling points: " << eIter->GetOverlapResult().GetNMatchedSamplingPoints() << std::endl;
            std::cout << " ---N. (re-)upsampled sampling points: " << eIter->GetOverlapResult().GetNReUpsampledSamplingPoints() << std::endl;
            std::cout << " ---N. (re-)upsampled matched sampling points: " << eIter->GetOverlapResult().GetNMatchedReUpsampledSamplingPoints()
                      << std::endl;
            std::cout << " ---Correlation coeff.: " << eIter->GetOverlapResult().GetCorrelationCoefficient() << std::endl;
            std::cout << " ---Locally matched fraction: " << eIter->GetOverlapResult().GetLocallyMatchedFraction() << std::endl;

            if (m_showEachIndividualElement)
            {
                const ClusterList clusterList1(1, eIter->GetCluster1()), clusterList2(1, eIter->GetCluster2());
                PANDORA_MONITORING_API(SetEveDisplayParameters(this->GetPandora(), false, DETECTOR_VIEW_XZ, -1.f, -1.f, 1.f));
                PANDORA_MONITORING_API(VisualizeClusters(this->GetPandora(), &allClusterList1, "AllClusters1", LIGHTORANGE));
                PANDORA_MONITORING_API(VisualizeClusters(this->GetPandora(), &allClusterList2, "AllClusters2", LIGHTYELLOW));
                PANDORA_MONITORING_API(VisualizeClusters(this->GetPandora(), &clusterList1, "Cluster1", RED));
                PANDORA_MONITORING_API(VisualizeClusters(this->GetPandora(), &clusterList2, "Cluster2", GREEN));
                PANDORA_MONITORING_API(ViewEvent(this->GetPandora()));
            }
        }

        std::cout << " All Connected Clusters " << std::endl;
        PANDORA_MONITORING_API(SetEveDisplayParameters(this->GetPandora(), false, DETECTOR_VIEW_XZ, -1.f, -1.f, 1.f));
        PANDORA_MONITORING_API(VisualizeClusters(this->GetPandora(), &allClusterList1, "AllClusters1", RED));
        PANDORA_MONITORING_API(VisualizeClusters(this->GetPandora(), &allClusterList2, "AllClusters2", GREEN));

        PANDORA_MONITORING_API(ViewEvent(this->GetPandora()));
    }

    return false;
}

Member Data Documentation

bool lar_content::TransverseMatrixVisualizationTool::m_ignoreUnavailableClusters

private

Whether to ignore (skip-over) unavailable clusters in the matrix.

Definition at line 33 of file TransverseMatrixVisualizationTool.h.

Referenced by ReadSettings(), and Run().

unsigned int lar_content::TransverseMatrixVisualizationTool::m_minClusterConnections

private

The minimum number of cluster connections for display.

Definition at line 32 of file TransverseMatrixVisualizationTool.h.

Referenced by ReadSettings(), and Run().

bool lar_content::TransverseMatrixVisualizationTool::m_showEachIndividualElement

private

Whether to draw each individual matrix element.

Definition at line 34 of file TransverseMatrixVisualizationTool.h.

Referenced by ReadSettings(), and Run().

bool lar_content::TransverseMatrixVisualizationTool::m_showOnlyTrueMatchIndividualElements

private

Whether to draw only truly matching individual matrix elements.

Definition at line 35 of file TransverseMatrixVisualizationTool.h.

Referenced by ReadSettings(), and Run().

---

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

• larpandoracontent/v04_08_01/source/larpandoracontent/LArMonitoring/TransverseMatrixVisualizationTool.h
• larpandoracontent/v04_08_01/source/larpandoracontent/LArMonitoring/TransverseMatrixVisualizationTool.cc