42 for (titer = trange.first; titer != trange.second; ++titer)
44 if (titer->first < imax)
49 cout <<
"min index: " << imax << endl;
50 std::pair< std::map<int,int>::const_iterator, std::map<int,int>::const_iterator > embtrk_b_e = truthcont->
GetEmbeddedTrkIds();
51 std::map<int,int>::const_iterator embiter;
52 for (embiter = embtrk_b_e.first; embiter != embtrk_b_e.second; ++embiter)
54 cout <<
"embedded trkid: " << embiter->first << endl;
56 PHG4HitContainer *ghit = findNode::getClass<PHG4HitContainer>(topNode,
"G4HIT_CEMC_E");
62 for (hit = hit_begin_end.first; hit != hit_begin_end.second; ++hit)
65 int trkid = hit->second->get_trkid();
70 cout <<
"could not locate geant particle " << trkid <<
" in G4HIT_CEMC_E" << endl;
78 if (printpart.find(primary_id) == printpart.end())
80 cout <<
"primary id " << primary_id <<
" is embedded" << endl;
81 printpart.insert(primary_id);
89 ghit = findNode::getClass<PHG4HitContainer>(topNode,
"G4HIT_SVTX");
94 for (hit = hit_begin_end.first; hit != hit_begin_end.second; ++hit)
96 int trkid = hit->second->get_trkid();
100 cout <<
"could not locate geant particle " << trkid <<
" in G4HIT_SVTX" << endl;