30 , _truth_container(NULL)
43 _truth_container = findNode::getClass<PHG4TruthInfoContainer>(topNode,
"G4TruthInfo");
46 cout <<
PHWHERE <<
" PHG4TruthInfoContainer node not found on node tree"
50 ostringstream nodename;
55 PHG4HitContainer *cylinder_hits = findNode::getClass<PHG4HitContainer>(topNode, nodename.str());
66 ostringstream nodename;
73 PHG4HitContainer *cylinder_hits = findNode::getClass<PHG4HitContainer>(topNode, nodename.str().c_str());
76 cerr <<
"can't find " << nodename.str() << endl;
83 itr != itr_range.second; ++itr)
86 if (!particle)
continue;
90 layerit != hits->
getLayers().second; layerit++)
96 if (hit_iter->second->get_trkid() != trkid)
continue;
97 if (hit_iter->second->get_hit_type() == 1)
99 else if (hit_iter->second->get_hit_type() == 2)
100 _hit_Z = hit_iter->second;
115 if (h1 ==
nullptr || h2 ==
nullptr)
return nullptr;
127 merged->
set_x(0, R_in * TMath::Cos(phi_measured));
128 merged->
set_y(0, R_in * TMath::Sin(phi_measured));
130 merged->
set_x(1, R_out * TMath::Cos(phi_measured));
131 merged->
set_y(1, R_out * TMath::Sin(phi_measured));
133 float t0 = h1->
get_t(0);
137 merged->
set_t(0, t0);
138 merged->
set_t(1, t1);