38 #include <TMatrixFfwd.h>
40 #include <TMatrixTUtils.h>
57 auto dstNode =
dynamic_cast<PHCompositeNode *
>(iter.findFirst(
"PHCompositeNode",
"DST"));
60 auto DetNode =
dynamic_cast<PHCompositeNode *
>(dstiter.findFirst(
"PHCompositeNode",
"TRKR"));
64 dstNode->addNode(DetNode);
73 DetNode->addNode(newNode);
76 m_clusters = findNode::getClass<TrkrClusterContainer>(
m_topNode,
"TRKR_TRUTHCLUSTERCONTAINER");
81 DetNode->addNode(newNode);
87 std::cout <<
PHWHERE <<
" PHG4TruthInfoContainer node not found on node tree" << std::endl;
97 int new_trkid = (hit==
nullptr) ? -1 : hit->
get_trkid();
99 if (
m_verbosity>5) std::cout <<
PHWHERE << std::endl <<
" -> Checking status of PHG4Hit. Track id("<<new_trkid<<
")" << std::endl;
114 for (
auto cluster = range.first; cluster != range.second; ++cluster) {
116 m_clusters->addClusterSpecifyKey(ckey, cluster->second);
138 hit = hitsetit->second->getHit(hitkey);
143 hitsetit->second->addHitSpecificKey(hitkey, hit);
150 std::cout <<
PHWHERE <<
": content of clusters " << std::endl;
151 auto& tmap = m_truthtracks->getMap();
152 std::cout <<
" Number of tracks: " << tmap.size() << std::endl;
153 for (
auto& _pair : tmap) {
154 auto& track = _pair.second;
156 printf(
"id(%2i) phi:eta:pt(", (
int)track->getTrackid());
157 std::cout <<
"phi:eta:pt(";
158 printf(
"%5.2f:%5.2f:%5.2f", track->getPhi(), track->getPseudoRapidity(), track->getPt());
161 std::cout <<
") nclusters(" << track->getClusters().size() <<
") ";
165 for (
auto cluskey : track->getClusters()) {
170 if (nclusprint > 0 && nclus >= nclusprint) {
171 std::cout <<
" ... ";
177 std::cout <<
PHWHERE <<
" ----- end of clusters " << std::endl;