52 idsmap = std::vector<std::vector<size_t>>(
planes.size(), std::vector<size_t>());
53 for (
auto h : hitlist) {
54 idsmap[h->View()].push_back(h.key());
57 vector<int32_t> hit_table_hit_id_data;
58 vector<int32_t> hit_table_local_plane_data;
59 vector<float> hit_table_local_time_data;
60 vector<int32_t> hit_table_local_wire_data;
61 vector<float> hit_table_integral_data;
62 vector<float> hit_table_rms_data;
63 vector<int32_t> spacepoint_table_spacepoint_id_data;
64 vector<int32_t> spacepoint_table_hit_id_u_data;
65 vector<int32_t> spacepoint_table_hit_id_v_data;
66 vector<int32_t> spacepoint_table_hit_id_y_data;
69 for (
auto h : hitlist) {
70 hit_table_hit_id_data.push_back(h.key());
71 hit_table_local_plane_data.push_back(h->View());
72 hit_table_local_time_data.push_back(h->PeakTime());
73 hit_table_local_wire_data.push_back(h->WireID().Wire);
74 hit_table_integral_data.push_back(h->Integral());
75 hit_table_rms_data.push_back(h->RMS());
80 vector<art::Ptr<recob::SpacePoint>> splist;
83 vector<vector<art::Ptr<recob::Hit>>> sp2Hit(splist.size());
84 if (splist.size() > 0) {
86 for (
size_t spIdx = 0; spIdx < sp2Hit.size(); ++spIdx) {
87 sp2Hit[spIdx] = fmp.at(spIdx);
92 for (
size_t i = 0; i < splist.size(); ++i) {
93 spacepoint_table_spacepoint_id_data.push_back(i);
94 spacepoint_table_hit_id_u_data.push_back(-1);
95 spacepoint_table_hit_id_v_data.push_back(-1);
96 spacepoint_table_hit_id_y_data.push_back(-1);
97 for (
size_t j = 0; j < sp2Hit[i].size(); ++j) {
98 if (sp2Hit[i][j]->View() == 0) spacepoint_table_hit_id_u_data.back() = sp2Hit[i][j].key();
99 if (sp2Hit[i][j]->View() == 1) spacepoint_table_hit_id_v_data.back() = sp2Hit[i][j].key();
100 if (sp2Hit[i][j]->View() == 2) spacepoint_table_hit_id_y_data.back() = sp2Hit[i][j].key();
104 inputs.emplace_back(
"hit_table_hit_id", hit_table_hit_id_data);
105 inputs.emplace_back(
"hit_table_local_plane", hit_table_local_plane_data);
106 inputs.emplace_back(
"hit_table_local_time", hit_table_local_time_data);
107 inputs.emplace_back(
"hit_table_local_wire", hit_table_local_wire_data);
108 inputs.emplace_back(
"hit_table_integral", hit_table_integral_data);
109 inputs.emplace_back(
"hit_table_rms", hit_table_rms_data);
111 inputs.emplace_back(
"spacepoint_table_spacepoint_id", spacepoint_table_spacepoint_id_data);
112 inputs.emplace_back(
"spacepoint_table_hit_id_u", spacepoint_table_hit_id_u_data);
113 inputs.emplace_back(
"spacepoint_table_hit_id_v", spacepoint_table_hit_id_v_data);
114 inputs.emplace_back(
"spacepoint_table_hit_id_y", spacepoint_table_hit_id_y_data);
bool getByLabel(std::string const &label, std::string const &instance, Handle< PROD > &result) const
void fill_ptr_vector(std::vector< Ptr< T >> &ptrs, H const &h)