#include <ProtoDUNEEmptyEventFinder.h>

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

bool	IsEmptyEvent (const art::Event &evt) const

Private Attributes
std::string	fSpacePointLabel

unsigned int	fNHitsThreshold

std::vector< float >	fMinCoordsData

std::vector< float >	fMaxCoordsData

std::vector< float >	fMinCoordsSim

std::vector< float >	fMaxCoordsSim

Detailed Description

Definition at line 12 of file ProtoDUNEEmptyEventFinder.h.

Constructor & Destructor Documentation

protoana::ProtoDUNEEmptyEventFinder::ProtoDUNEEmptyEventFinder ( const fhicl::ParameterSet & pset )

Definition at line 12 of file ProtoDUNEEmptyEventFinder.cxx.

                                                                                           :
   fSpacePointLabel(pset.get<std::string>("SpacePointLabel","reco3d")),
   fNHitsThreshold(pset.get<unsigned int>("NHitsThreshold",10)),
   fMinCoordsData(pset.get<std::vector<float>>("MinCoordValuesData")),
   fMaxCoordsData(pset.get<std::vector<float>>("MaxCoordValuesData")),
   fMinCoordsSim(pset.get<std::vector<float>>("MinCoordValuesSim")),
   fMaxCoordsSim(pset.get<std::vector<float>>("MaxCoordValuesSim"))
 { 
 }

Member Function Documentation

bool protoana::ProtoDUNEEmptyEventFinder::IsEmptyEvent ( const art::Event & evt ) const

Definition at line 22 of file ProtoDUNEEmptyEventFinder.cxx.

 {
   // Get the space points
   auto spacePointHandle = evt.getValidHandle<std::vector<recob::SpacePoint>>(fSpacePointLabel);
   std::vector<art::Ptr<recob::SpacePoint>> spacePoints;
   if(spacePointHandle.isValid())
     art::fill_ptr_vector(spacePoints,spacePointHandle);
 
   const std::vector<float> minValues{evt.isRealData() ? fMinCoordsData : fMinCoordsSim};
   const std::vector<float> maxValues{evt.isRealData() ? fMaxCoordsData : fMaxCoordsSim};
 
   // Now count the number of space points
   unsigned int nBeamRegionSpacePoints{0};
   for (const art::Ptr<recob::SpacePoint> &sp : spacePoints)
   {
     const double *xyz{sp->XYZ()};
 
     if(xyz[0] < minValues.at(0) || xyz[0] > maxValues.at(0)) continue;
     if(xyz[1] < minValues.at(1) || xyz[1] > maxValues.at(1)) continue;
     if(xyz[2] < minValues.at(2) || xyz[2] > maxValues.at(2)) continue;
 
     ++nBeamRegionSpacePoints;
 
     // If we reach the threshold then there is no point continuing
     if (nBeamRegionSpacePoints >= fNHitsThreshold) return false;
   }
 
   // If we get here then we didn't meet the threshold
   return true;
 }