Inheritance diagram for cluster::SmallClusterFilter:

Public Member Functions
	SmallClusterFilter (fhicl::ParameterSet const &pset)

Private Member Functions
bool	filter (art::Event &evt) override

Private Attributes
std::string	fHitFinderModuleLabel
	label of module making hits More...

std::vector< std::size_t >	fMaxHitsByPlane
	maximum hits on each plane More...

std::size_t	fMaxTotalHits
	maximum number of hits allowed More...

std::size_t	fNPlanes
	number of planes More...

Detailed Description

Definition at line 36 of file SmallClusterFilter_module.cc.

Constructor & Destructor Documentation

cluster::SmallClusterFilter::SmallClusterFilter ( fhicl::ParameterSet const & pset )

explicit

Definition at line 51 of file SmallClusterFilter_module.cc.

   : EDFilter{pset}
   , fHitFinderModuleLabel{pset.get<std::string>("HitFinderModuleLabel")}
   , fMaxHitsByPlane{pset.get<std::vector<std::size_t>>("MaxHitsByPlane")}
   , fMaxTotalHits{pset.get<std::size_t>("MaxTotalHits")}
   , fNPlanes{art::ServiceHandle<geo::Geometry const>()->Nplanes()}
 {
   if (size(fMaxHitsByPlane) != fNPlanes) {
     throw art::Exception{art::errors::Configuration}
       << "Mismatch in number of planes specified in 'MaxHitsByPlane' (" << size(fMaxHitsByPlane)
       << ") and the number of planes (" << fNPlanes << ")\n";
   }
 }

Member Function Documentation

bool cluster::SmallClusterFilter::filter ( art::Event & evt )

overrideprivate

Definition at line 66 of file SmallClusterFilter_module.cc.

 {
   auto const& hits = *evt.getValidHandle<std::vector<recob::Hit>>(fHitFinderModuleLabel);
 
   mf::LogVerbatim("SmallClusterFilter")
     << " ++++ Hitsreceived received " << size(hits) << " +++++ ";
 
   if (empty(hits)) {
     mf::LogWarning("SmallClusterFilter") << " no hits received! exiting ";
     return false;
   }
 
   if (size(hits) > fMaxTotalHits) {
     mf::LogWarning("SmallClusterFinder") << "Not an empty event, exiting.";
     return false;
   }
 
   bool collFound = false;
 
   std::map<unsigned int, std::size_t> hitsPerPlane;
 
   for (auto const& hit : hits) {
 
     // Skip crazy hits:
     if (hit.Integral() > 500) continue;
 
     ++hitsPerPlane[hit.WireID().Plane];
 
     if (hit.SignalType() == geo::kCollection) collFound = true;
 
   }
 
   for (unsigned int i = 0; i < fNPlanes; i++) {
     if (hitsPerPlane[i] > fMaxHitsByPlane[i]) return false;
   }
 
   //check that there is at least 1 hit on collection:
   return collFound;
 }