40 ntup =
new TNtuple(
"hitntup",
"G4Hits",
"detid:lyr:slat:x0:y0:z0:x1:y1:z1:edep");
42 TH1 *
h1 =
new TH1F(
"edep1GeV",
"edep 0-1GeV", 1000, 0, 1);
44 h1 =
new TH1F(
"edep100GeV",
"edep 0-100GeV", 1000, 0, 100);
51 ostringstream nodename;
52 set<string>::const_iterator iter;
53 vector<TH1 *>::const_iterator eiter;
58 nodename <<
"G4HIT_" << *iter;
69 esum += hit_iter->second->get_edep();
71 hit_iter->second->get_layer(),
72 hit_iter->second->get_scint_id(),
73 hit_iter->second->get_x(0),
74 hit_iter->second->get_y(0),
75 hit_iter->second->get_z(0),
76 hit_iter->second->get_x(1),
77 hit_iter->second->get_y(1),
78 hit_iter->second->get_z(1),
79 hit_iter->second->get_edep());
81 for (eiter =
eloss.begin(); eiter !=
eloss.end(); ++eiter)