9 #include "Pandora/AlgorithmHeaders.h" 18 TransverseTensorVisualizationTool::TransverseTensorVisualizationTool() :
19 m_minClusterConnections(1),
20 m_ignoreUnavailableClusters(true),
21 m_showEachIndividualElement(false),
30 if (PandoraContentApi::GetSettings(*pAlgorithm)->ShouldDisplayAlgorithmInfo())
31 std::cout <<
"----> Running Algorithm Tool: " << this->GetInstanceName() <<
", " << this->GetType() << std::endl;
33 ClusterSet usedKeyClusters;
37 for (
const Cluster *
const pKeyCluster : sortedKeyClusters)
42 if (usedKeyClusters.count(pKeyCluster))
45 unsigned int nU(0), nV(0), nW(0);
52 if (nU * nV * nW == 0)
56 ClusterList allClusterListU, allClusterListV, allClusterListW;
57 std::cout <<
" Connections: nU " << nU <<
", nV " << nV <<
", nW " << nW <<
", nElements " << elementList.size() << std::endl;
61 if (allClusterListU.end() == std::find(allClusterListU.begin(), allClusterListU.end(), eIter->GetClusterU()))
62 allClusterListU.push_back(eIter->GetClusterU());
63 if (allClusterListV.end() == std::find(allClusterListV.begin(), allClusterListV.end(), eIter->GetClusterV()))
64 allClusterListV.push_back(eIter->GetClusterV());
65 if (allClusterListW.end() == std::find(allClusterListW.begin(), allClusterListW.end(), eIter->GetClusterW()))
66 allClusterListW.push_back(eIter->GetClusterW());
67 usedKeyClusters.insert(eIter->GetClusterU());
69 std::cout <<
" Element " << counter++ <<
": MatchedFraction " << eIter->GetOverlapResult().GetMatchedFraction()
70 <<
", MatchedSamplingPoints " << eIter->GetOverlapResult().GetNMatchedSamplingPoints() <<
", xSpanU " 71 << eIter->GetOverlapResult().GetXOverlap().GetXSpanU() <<
", xSpanV " << eIter->GetOverlapResult().GetXOverlap().GetXSpanV()
72 <<
", xSpanW " << eIter->GetOverlapResult().GetXOverlap().GetXSpanW() <<
", xOverlapSpan " 73 << eIter->GetOverlapResult().GetXOverlap().GetXOverlapSpan() << std::endl;
77 const ClusterList clusterListU(1, eIter->GetClusterU()), clusterListV(1, eIter->GetClusterV()),
78 clusterListW(1, eIter->GetClusterW());
79 PANDORA_MONITORING_API(SetEveDisplayParameters(this->GetPandora(),
false, DETECTOR_VIEW_XZ, -1.
f, -1.
f, 1.
f));
80 PANDORA_MONITORING_API(VisualizeClusters(this->GetPandora(), &clusterListU,
"UCluster", RED));
81 PANDORA_MONITORING_API(VisualizeClusters(this->GetPandora(), &clusterListV,
"VCluster", GREEN));
82 PANDORA_MONITORING_API(VisualizeClusters(this->GetPandora(), &clusterListW,
"WCluster", BLUE));
83 PANDORA_MONITORING_API(ViewEvent(this->GetPandora()));
87 std::cout <<
" All Connected Clusters " << std::endl;
88 PANDORA_MONITORING_API(SetEveDisplayParameters(this->GetPandora(),
false, DETECTOR_VIEW_XZ, -1.
f, -1.
f, 1.
f));
89 PANDORA_MONITORING_API(VisualizeClusters(this->GetPandora(), &allClusterListU,
"AllUClusters", RED));
90 PANDORA_MONITORING_API(VisualizeClusters(this->GetPandora(), &allClusterListV,
"AllVClusters", GREEN));
91 PANDORA_MONITORING_API(VisualizeClusters(this->GetPandora(), &allClusterListW,
"AllWClusters", BLUE));
95 PANDORA_MONITORING_API(VisualizeClusters(this->GetPandora(), &(pAlgorithm->
GetInputClusterList(TPC_VIEW_U)),
"InputClusterListU", GRAY));
96 PANDORA_MONITORING_API(VisualizeClusters(this->GetPandora(), &(pAlgorithm->
GetInputClusterList(TPC_VIEW_V)),
"InputClusterListV", GRAY));
97 PANDORA_MONITORING_API(VisualizeClusters(this->GetPandora(), &(pAlgorithm->
GetInputClusterList(TPC_VIEW_W)),
"InputClusterListW", GRAY));
100 PANDORA_MONITORING_API(ViewEvent(this->GetPandora()));
110 PANDORA_RETURN_RESULT_IF_AND_IF(
111 STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle,
"MinClusterConnections",
m_minClusterConnections));
113 PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
116 PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
119 PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle,
"ShowContext",
m_showContext));
121 return STATUS_CODE_SUCCESS;
bool m_showContext
Whether to show input cluster lists to add context to tensor elements.
bool m_ignoreUnavailableClusters
Whether to ignore (skip-over) unavailable clusters in the tensor.
void GetConnectedElements(const pandora::Cluster *const pCluster, const bool ignoreUnavailable, ElementList &elementList) const
Get a list of elements connected to a specified cluster.
bool Run(ThreeViewTransverseTracksAlgorithm *const pAlgorithm, TensorType &overlapTensor)
Run the algorithm tool.
pandora::StatusCode ReadSettings(const pandora::TiXmlHandle xmlHandle)
std::vector< Element > ElementList
auto counter(T begin, T end)
Returns an object to iterate values from begin to end in a range-for loop.
const pandora::ClusterList & GetInputClusterList(const pandora::HitType hitType) const
Get the input cluster list corresponding to a specified hit type.
void GetSortedKeyClusters(pandora::ClusterVector &sortedKeyClusters) const
Get a sorted vector of key clusters (U clusters with current implementation)
bool m_showEachIndividualElement
Whether to draw each individual tensor element.
unsigned int m_minClusterConnections
The minimum number of cluster connections for display.
std::vector< art::Ptr< recob::Cluster > > ClusterVector
ThreeViewTransverseTracksAlgorithm class.
TheTensor::const_iterator const_iterator