Inheritance diagram for simfilter::FilterSimPhotonTime:

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

	FilterSimPhotonTime (FilterSimPhotonTime const &)=delete

	FilterSimPhotonTime (FilterSimPhotonTime &&)=delete

FilterSimPhotonTime &	operator= (FilterSimPhotonTime const &)=delete

FilterSimPhotonTime &	operator= (FilterSimPhotonTime &&)=delete

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

void	CheckTimeWindows () const

Private Attributes
std::string const	fSimPhotonsCollectionLabel

std::vector< std::pair< float, float > > const	fTimeWindows

float const	fMinTotalEnergy

float const	fMinPhotonEnergy

bool const	fDebug

std::size_t const	fN

bool const	fUseReflectedPhotons

std::string const	fReflectedLabel

Detailed Description

Definition at line 28 of file FilterSimPhotonTime_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 53 of file FilterSimPhotonTime_module.cc.

   : SharedFilter{p}
   , fSimPhotonsCollectionLabel(p.get<std::string>("SimPhotonsCollectionLabel"))
   , fTimeWindows(p.get<std::vector<std::pair<float, float>>>("TimeWindows"))
   , fMinTotalEnergy(p.get<float>("MinTotalEnergy", 0.0))
   , fMinPhotonEnergy(p.get<float>("MinPhotonEnergy", -1))
   , 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::FilterSimPhotonTime::FilterSimPhotonTime ( FilterSimPhotonTime const & )

delete

simfilter::FilterSimPhotonTime::FilterSimPhotonTime ( FilterSimPhotonTime && )

delete

Member Function Documentation

void simfilter::FilterSimPhotonTime::CheckTimeWindows ( ) const

private

Definition at line 78 of file FilterSimPhotonTime_module.cc.

 {
 
   if (fDebug)
     std::cout << "\tFilterSimPhotonTime: 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("FilterSimPhotonTime::CheckTimeWindows")
         << "Bad time window initialization: tw.first>tw.second. Reverse the order!" << std::endl;
   }
 }

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

overrideprivate

Definition at line 96 of file FilterSimPhotonTime_module.cc.

 {
   auto const& simPhotonsCollection =
     *e.getValidHandle<std::vector<sim::SimPhotons>>(fSimPhotonsCollectionLabel);
 
   std::vector<double> sumEnergyArray(fN, 0.0);
 
   const std::vector<sim::SimPhotons>& simPhotonsCollectionReflected =
     fUseReflectedPhotons ? *e.getValidHandle<std::vector<sim::SimPhotons>>(
                              {fSimPhotonsCollectionLabel, fReflectedLabel}) :
                            std::vector<sim::SimPhotons>();
 
   size_t n_sim_photons = simPhotonsCollection.size() + simPhotonsCollectionReflected.size();
 
   for (size_t i_pc = 0; i_pc < n_sim_photons; i_pc++) {
     const sim::SimPhotons& simphotons =
       (i_pc < simPhotonsCollection.size()) ?
         simPhotonsCollection[i_pc] :
         simPhotonsCollectionReflected[i_pc - simPhotonsCollection.size()];
 
     if (fDebug)
       std::cout << "\tFilterSimPhotonTime: Processing simphotons channel " << simphotons.OpChannel()
                 << std::endl;
 
     for (auto const& photon : simphotons)
       for (size_t i_tw = 0; i_tw < fN; i_tw++) {
         auto const& tw(fTimeWindows[i_tw]);
         if (photon.Time >= tw.first && photon.Time <= tw.second &&
             photon.Energy > fMinPhotonEnergy) {
 
           if (fDebug) {
             std::string photon_string =
               (i_pc < simPhotonsCollection.size()) ? "Photon" : "Reflected Photon";
             std::cout << "\t\t" << photon_string << " with time " << photon.Time << " detected. "
                       << "Energy is  " << photon.Energy << "." << std::endl;
           }
 
           sumEnergyArray[i_tw] += photon.Energy;
 
           if (fDebug)
             std::cout << "\t\tTotal energy in this window (" << i_tw << ") is now "
                       << sumEnergyArray[i_tw] << std::endl;
 
           if (sumEnergyArray[i_tw] > fMinTotalEnergy) return true;
         }
       }
   }
 
   if (fDebug) {
     std::cout << "\tFilterSimPhotonTime: Final total energies are below min of " << fMinTotalEnergy
               << ":" << 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
                 << "]: " << sumEnergyArray[i_tw] << std::endl;
     }
   }
 
   return false;
 }