Inheritance diagram for simfilter::FilterSimPhotonLiteTime:

Public Member Functions
	FilterSimPhotonLiteTime (fhicl::ParameterSet const &p, art::ProcessingFrame const &)

	FilterSimPhotonLiteTime (FilterSimPhotonLiteTime const &)=delete

	FilterSimPhotonLiteTime (FilterSimPhotonLiteTime &&)=delete

FilterSimPhotonLiteTime &	operator= (FilterSimPhotonLiteTime const &)=delete

FilterSimPhotonLiteTime &	operator= (FilterSimPhotonLiteTime &&)=delete

Private Member Functions
bool	filter (art::Event &e, art::ProcessingFrame const &) override

void	CheckTimeWindows () const

Private Attributes
std::string const	fSimPhotonsLiteCollectionLabel
	Label for the sim::SimPhotonsLite data product. More...

std::vector< std::pair< int, int > > const	fTimeWindows
	Time windows used for filtering. Units are the same as in the sim::SimPhotonsLite. More...

int const	fMinTotalPhotons
	Minimum number of photons inside a window to pass the filter. More...

bool const	fDebug
	Set to true to print (a lot of) debug information. More...

std::size_t const	fN
	Number of time winows. More...

bool const	fUseReflectedPhotons
	Whether to include reflected photons in the filter. More...

std::string const	fReflectedLabel
	Label for the reflected photons – "Reflected" by default. More...

Detailed Description

Definition at line 32 of file FilterSimPhotonLiteTime_module.cc.

Constructor & Destructor Documentation

simfilter::FilterSimPhotonLiteTime::FilterSimPhotonLiteTime	(	fhicl::ParameterSet const &	p,
		art::ProcessingFrame const &
	)

explicit

Definition at line 58 of file FilterSimPhotonLiteTime_module.cc.

   : SharedFilter{p}
   , fSimPhotonsLiteCollectionLabel(p.get<std::string>("SimPhotonsLiteCollectionLabel"))
   , fTimeWindows(p.get<std::vector<std::pair<int, int>>>("TimeWindows"))
   , fMinTotalPhotons(p.get<int>("MinTotalPhotons"))
   , fDebug(p.get<bool>("Debug", false))
   , fN(fTimeWindows.size())
   , fUseReflectedPhotons(p.get<bool>("UseReflectedPhotons", false))
   , fReflectedLabel(p.get<std::string>("ReflectedLabel", "Reflected"))
 {
   CheckTimeWindows();
 
   // For printing out debug messages, we want to serialize the
   // event-level calls so that the messages are not garbled.
   // Otherwise, this module works well for asynchronous event-level
   // calls.
   if (fDebug) { serialize(); }
   else {
     async<art::InEvent>();
   }
 }

simfilter::FilterSimPhotonLiteTime::FilterSimPhotonLiteTime ( FilterSimPhotonLiteTime const & )

delete

simfilter::FilterSimPhotonLiteTime::FilterSimPhotonLiteTime ( FilterSimPhotonLiteTime && )

delete

Member Function Documentation

void simfilter::FilterSimPhotonLiteTime::CheckTimeWindows ( ) const

private

Definition at line 82 of file FilterSimPhotonLiteTime_module.cc.

 {
 
   if (fDebug)
     std::cout << "\tFilterSimPhotonLiteTime: TimeWindows size is " << fTimeWindows.size()
               << std::endl;
 
   for (auto const& tw : fTimeWindows) {
     if (fDebug)
       std::cout << "\t\tTimeWindow "
                 << "[" << tw.first << "," << tw.second << "]" << std::endl;
 
     if (tw.first > tw.second)
       throw cet::exception("FilterSimPhotonLiteTime::CheckTimeWindows")
         << "Bad time window initialization: tw.first>tw.second. Reverse the order!" << std::endl;
   }
 }

bool simfilter::FilterSimPhotonLiteTime::filter	(	art::Event &	e,
		art::ProcessingFrame const &
	)

overrideprivate

Definition at line 101 of file FilterSimPhotonLiteTime_module.cc.

 {
   auto const& simPhotonsLiteCollection =
     *e.getValidHandle<std::vector<sim::SimPhotonsLite>>(fSimPhotonsLiteCollectionLabel);
 
   std::vector<int> sumNPhotonArray(fN, 0);
 
   const std::vector<sim::SimPhotonsLite>& simPhotonsLiteCollectionReflected =
     fUseReflectedPhotons ? *e.getValidHandle<std::vector<sim::SimPhotonsLite>>(
                              {fSimPhotonsLiteCollectionLabel, fReflectedLabel}) :
                            std::vector<sim::SimPhotonsLite>();
 
   size_t n_sim_photons = simPhotonsLiteCollection.size() + simPhotonsLiteCollectionReflected.size();
 
   if (fDebug) {
     std::cout << "New event to filter with total # sim photons: " << n_sim_photons << std::endl;
   }
 
   for (size_t i_pc = 0; i_pc < n_sim_photons; i_pc++) {
     const sim::SimPhotonsLite& simphotonslite =
       (i_pc < simPhotonsLiteCollection.size()) ?
         simPhotonsLiteCollection[i_pc] :
         simPhotonsLiteCollectionReflected[i_pc - simPhotonsLiteCollection.size()];
 
     if (fDebug)
       std::cout << "\tFilterSimPhotonLiteTime: Processing simphotonslite channel "
                 << simphotonslite.OpChannel << std::endl;
 
     for (auto const& photon_pair : simphotonslite.DetectedPhotons) {
       for (size_t i_tw = 0; i_tw < fN; i_tw++) {
         auto const& tw(fTimeWindows[i_tw]);
         if (photon_pair.first >= tw.first && photon_pair.first <= tw.second) {
 
           if (fDebug) {
             std::string photon_string =
               (i_pc < simPhotonsLiteCollection.size()) ? "Photon" : "Reflected Photon";
             std::cout << "\t\t" << photon_string << " with number " << photon_pair.second
                       << " at time " << photon_pair.first << " detected." << std::endl;
           }
 
           sumNPhotonArray[i_tw] += photon_pair.second;
 
           if (fDebug)
             std::cout << "\t\tTotal number of photons in this window (" << i_tw << ") is now "
                       << sumNPhotonArray[i_tw] << std::endl;
 
           if (sumNPhotonArray[i_tw] >= fMinTotalPhotons) return true;
         }
       }
     }
   }
 
   if (fDebug) {
     std::cout << "\tFilterSimPhotonLiteTime: Final total numbers are below min of "
               << fMinTotalPhotons << ":" << std::endl;
     for (size_t i_tw = 0; i_tw < fN; ++i_tw) {
       std::cout << "\t\tTimeWindow "
                 << "[" << fTimeWindows[i_tw].first << "," << fTimeWindows[i_tw].second
                 << "]: " << sumNPhotonArray[i_tw] << std::endl;
     }
   }
 
   return false;
 }