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

Public Member Functions inherited from art::SharedFilter
	SharedFilter (fhicl::ParameterSet const &pset)

template<typename Config >
	SharedFilter (Table< Config > const &config)

std::string	workerType () const

Public Member Functions inherited from art::detail::Filter
virtual	~Filter () noexcept

	Filter (fhicl::ParameterSet const &)

	Filter (Filter const &)=delete

	Filter (Filter &&)=delete

Filter &	operator= (Filter const &)=delete

Filter &	operator= (Filter &&)=delete

void	doBeginJob (SharedResources const &resources)

void	doEndJob ()

void	doRespondToOpenInputFile (FileBlock const &fb)

void	doRespondToCloseInputFile (FileBlock const &fb)

void	doRespondToOpenOutputFiles (FileBlock const &fb)

void	doRespondToCloseOutputFiles (FileBlock const &fb)

bool	doBeginRun (RunPrincipal &rp, ModuleContext const &)

bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)

Public Member Functions inherited from art::Modifier
	~Modifier () noexcept

	Modifier ()

	Modifier (Modifier const &)=delete

	Modifier (Modifier &&)=delete

Modifier &	operator= (Modifier const &)=delete

Modifier &	operator= (Modifier &&)=delete

Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept

	ModuleBase ()

ModuleDescription const &	moduleDescription () const

void	setModuleDescription (ModuleDescription const &)

std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const

void	sortConsumables (std::string const &current_process_name)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Public Member Functions inherited from art::detail::SharedModule
	SharedModule ()

	SharedModule (std::string const &moduleLabel)

hep::concurrency::SerialTaskQueueChain *	serialTaskQueueChain () const

std::set< std::string > const &	sharedResources () const

void	createQueues (SharedResources const &resources)

template<BranchType , typename... T>
void	serialize (T const &...resources)

template<BranchType , typename... T>
void	serializeExternal (T const &...resources)

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

Additional Inherited Members
Public Types inherited from art::SharedFilter
using	ModuleType = SharedFilter

using	WorkerType = WorkerT< SharedFilter >

Public Types inherited from art::detail::Filter
template<typename UserConfig >
using	Table = Modifier::Table< UserConfig >

Public Types inherited from art::Modifier
template<typename UserConfig , typename UserKeysToIgnore = void>
using	Table = ProducerTable< UserConfig, detail::ModuleConfig, UserKeysToIgnore >

Static Public Attributes inherited from art::detail::Filter
static constexpr bool	Pass {true}

static constexpr bool	Fail {false}

Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Protected Member Functions inherited from art::detail::SharedModule
template<BranchType BT = InEvent, typename... T>
void	serialize (T const &...)

template<BranchType BT = InEvent, typename... T>
void	serializeExternal (T const &...)

template<BranchType BT = InEvent>
void	async ()

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 &
	)

overrideprivatevirtual

Implements art::SharedFilter.

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;
 }