TransverseMatrixVisualizationTool class. More...

#include <TransverseMatrixVisualizationTool.h>

Inheritance diagram for lar_content::TransverseMatrixVisualizationTool:

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...

Additional Inherited Members
Public Types inherited from lar_content::TransverseMatrixTool
typedef TwoViewTransverseTracksAlgorithm::MatchingType::MatrixType	MatrixType

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

Detailed Description

TransverseMatrixVisualizationTool class.

Definition at line 19 of file TransverseMatrixVisualizationTool.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.

 {
     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.

 {
     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;
 }