305 std::vector<PiZeroProcess>
result;
308 auto showerHandle = evt.
getHandle<std::vector<recob::Shower>>(showerLabel);
309 if (!showerHandle || showerHandle->size() < 2)
return result;
312 double showerDistThreshold = 50;
313 std::vector<std::pair<const recob::Shower*, const recob::Shower*>> showerPairs;
314 for(
unsigned si = 0; si < showerHandle->size()-1; ++si) {
315 for(
unsigned ssi = si+1; ssi < showerHandle->size(); ++ssi) {
319 showerPairs.push_back(std::make_pair(shower1, shower2));
323 std::cout << showerPairs.size() <<
" shower pairs.\n";
const TVector3 & ShowerStart() const
Handle< PROD > getHandle(SelectorBase const &) const