Inheritance diagram for G4InfoReducer:

Public Types
using	ModuleType = EDProducer

template<typename UserConfig , typename KeysToIgnore = void>
using	Table = Modifier::Table< UserConfig, KeysToIgnore >

Public Member Functions
	G4InfoReducer (fhicl::ParameterSet const &p)

	G4InfoReducer (G4InfoReducer const &)=delete

	G4InfoReducer (G4InfoReducer &&)=delete

G4InfoReducer &	operator= (G4InfoReducer const &)=delete

G4InfoReducer &	operator= (G4InfoReducer &&)=delete

void	produce (art::Event &e) override

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

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)

void	fillProductDescriptions ()

void	registerProducts (ProductDescriptions &productsToRegister)

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)

std::unique_ptr< Worker >	makeWorker (WorkerParams const &wp)

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

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

Protected Member Functions
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 ()

Private Attributes
art::InputTag	fSedLabel
	module making the SimEnergyDeposit More...

double	fMinX

double	fMinY

double	fMinZ
	bottom left coordinate of union of all TPC active volumes More...

double	fVoxelSizeX

double	fVoxelSizeY

double	fVoxelSizeZ
	size of a voxel (cm) More...

double	fElectronDriftVel
	electron drift velocity (cm/us) More...

bool	fUseOrigTrackID

const geo::GeometryCore &	fGeometry

Detailed Description

Definition at line 32 of file G4InfoReducer_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 17 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::detail::Producer::Table = Modifier::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 26 of file Producer.h.

Constructor & Destructor Documentation

G4InfoReducer::G4InfoReducer ( fhicl::ParameterSet const & p )

explicit

Definition at line 58 of file G4InfoReducer_module.cc.

References fElectronDriftVel, fGeometry, fMinX, fMinY, fMinZ, fSedLabel, fUseOrigTrackID, fVoxelSizeX, fVoxelSizeY, fVoxelSizeZ, and geo::Iterable< IterationPolicy, Transform >::Iterate().

   : EDProducer{p} // ,
   , fGeometry(*lar::providerFrom<geo::Geometry>())
 {
   fSedLabel = p.get<art::InputTag>("SimEnergyDepositLabel", "largeant:TPCActive");
 
   // Prepare for voxelization
   fVoxelSizeX = p.get<double>("VoxelSizeX", 0.3);
   fVoxelSizeY = p.get<double>("VoxelSizeY", 0.3);
   fVoxelSizeZ = p.get<double>("VoxelSizeZ", 0.3);
 
   //Ionization electrons travel at ~1.6mm/us at 500 V/cm electric field --> 1.6/10 cm/us
   fElectronDriftVel = p.get<double>("ElectronDriftVel", 0.16); //cm/us
 
   //Use orig track id
   fUseOrigTrackID = p.get<bool>("useOrigTrackID", true);
 
   if (fVoxelSizeX <= 0. || fVoxelSizeY <= 0. || fVoxelSizeZ <= 0.) {
     std::cerr << "Voxel size must be strictly greater than zero." << std::endl;
     throw std::exception();
   }
 
   double min_x = std::numeric_limits<double>::max();
   double min_y = std::numeric_limits<double>::max();
   double min_z = std::numeric_limits<double>::max();
   for (geo::TPCGeo const& tpc : fGeometry.Iterate<geo::TPCGeo>()) {
     auto const& tpcabox = tpc.ActiveBoundingBox();
     min_x = std::min(min_x, tpcabox.MinX());
     min_y = std::min(min_y, tpcabox.MinY());
     min_z = std::min(min_z, tpcabox.MinZ());
   }
   auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataForJob();
   // Take into account TPC readout window size
   fMinX = min_x + clockData.TriggerOffsetTPC() * fElectronDriftVel;
   fMinY = min_y;
   fMinZ = min_z;
   //std::cout << "G4InfoReducer " << fMinX << " " << fMinY << " " << fMinZ << std::endl;
   //std::cout << clockData.TriggerOffsetTPC() << std::endl;
 
   // Call appropriate produces<>() functions here.
   produces<std::vector<sim::SimEnergyDepositLite>>();
 
   // Call appropriate consumes<>() for any products to be retrieved by this module.
   consumes<std::vector<sim::SimEnergyDeposit>>(fSedLabel);
 }

G4InfoReducer::G4InfoReducer ( G4InfoReducer const & )

delete

G4InfoReducer::G4InfoReducer ( G4InfoReducer && )

delete

Member Function Documentation

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::consumes ( InputTag const & tag )

protectedinherited

Definition at line 61 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumes().

   {
     return collector_.consumes<T, BT>(tag);
   }

ConsumesCollector & art::ModuleBase::consumesCollector ( )

protectedinherited

Definition at line 57 of file ModuleBase.cc.

References art::ModuleBase::collector_.

   {
     return collector_;
   }

template<typename T , BranchType BT>

void art::ModuleBase::consumesMany ( )

protectedinherited

Definition at line 75 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesMany().

   {
     collector_.consumesMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::consumesView ( InputTag const & )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::consumesView ( InputTag const & tag )

inherited

Definition at line 68 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesView().

   {
     return collector_.consumesView<T, BT>(tag);
   }

void art::detail::Producer::doBeginJob ( SharedResources const & resources )

inherited

Definition at line 22 of file Producer.cc.

References art::detail::Producer::beginJobWithFrame(), and art::detail::Producer::setupQueues().

   {
     setupQueues(resources);
     ProcessingFrame const frame{ScheduleID{}};
     beginJobWithFrame(frame);
   }

bool art::detail::Producer::doBeginRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 65 of file Producer.cc.

References art::detail::Producer::beginRunWithFrame(), art::RangeSet::forRun(), art::RunPrincipal::makeRun(), r, art::RunPrincipal::runID(), and art::ModuleContext::scheduleID().

   {
     auto r = rp.makeRun(mc, RangeSet::forRun(rp.runID()));
     ProcessingFrame const frame{mc.scheduleID()};
     beginRunWithFrame(r, frame);
     r.commitProducts();
     return true;
   }

bool art::detail::Producer::doBeginSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 85 of file Producer.cc.

References art::detail::Producer::beginSubRunWithFrame(), art::RangeSet::forSubRun(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::SubRunPrincipal::subRunID().

   {
     auto sr = srp.makeSubRun(mc, RangeSet::forSubRun(srp.subRunID()));
     ProcessingFrame const frame{mc.scheduleID()};
     beginSubRunWithFrame(sr, frame);
     sr.commitProducts();
     return true;
   }

void art::detail::Producer::doEndJob ( )

inherited

Definition at line 30 of file Producer.cc.

References art::detail::Producer::endJobWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     endJobWithFrame(frame);
   }

bool art::detail::Producer::doEndRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 75 of file Producer.cc.

References art::detail::Producer::endRunWithFrame(), art::RunPrincipal::makeRun(), r, art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

   {
     auto r = rp.makeRun(mc, rp.seenRanges());
     ProcessingFrame const frame{mc.scheduleID()};
     endRunWithFrame(r, frame);
     r.commitProducts();
     return true;
   }

bool art::detail::Producer::doEndSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 95 of file Producer.cc.

References art::detail::Producer::endSubRunWithFrame(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

   {
     auto sr = srp.makeSubRun(mc, srp.seenRanges());
     ProcessingFrame const frame{mc.scheduleID()};
     endSubRunWithFrame(sr, frame);
     sr.commitProducts();
     return true;
   }

bool art::detail::Producer::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
	)

inherited

Definition at line 105 of file Producer.cc.

References art::detail::Producer::checkPutProducts_, e, art::EventPrincipal::makeEvent(), art::detail::Producer::produceWithFrame(), and art::ModuleContext::scheduleID().

   {
     auto e = ep.makeEvent(mc);
     ++counts_run;
     ProcessingFrame const frame{mc.scheduleID()};
     produceWithFrame(e, frame);
     e.commitProducts(checkPutProducts_, &expectedProducts<InEvent>());
     ++counts_passed;
     return true;
   }

void art::detail::Producer::doRespondToCloseInputFile ( FileBlock const & fb )

inherited

Definition at line 44 of file Producer.cc.

References art::detail::Producer::respondToCloseInputFileWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseInputFileWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToCloseOutputFiles ( FileBlock const & fb )

inherited

Definition at line 58 of file Producer.cc.

References art::detail::Producer::respondToCloseOutputFilesWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseOutputFilesWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToOpenInputFile ( FileBlock const & fb )

inherited

Definition at line 37 of file Producer.cc.

References art::detail::Producer::respondToOpenInputFileWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenInputFileWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToOpenOutputFiles ( FileBlock const & fb )

inherited

Definition at line 51 of file Producer.cc.

References art::detail::Producer::respondToOpenOutputFilesWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenOutputFilesWithFrame(fb, frame);
   }

void art::Modifier::fillProductDescriptions ( )

inherited

Definition at line 10 of file Modifier.cc.

References art::ProductRegistryHelper::fillDescriptions(), and art::ModuleBase::moduleDescription().

   {
     ProductRegistryHelper::fillDescriptions(moduleDescription());
   }

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

inherited

Definition at line 43 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::getConsumables().

   {
     return collector_.getConsumables();
   }

std::unique_ptr< Worker > art::ModuleBase::makeWorker ( WorkerParams const & wp )

inherited

Definition at line 37 of file ModuleBase.cc.

References art::ModuleBase::doMakeWorker(), and art::NumBranchTypes.

   {
     return doMakeWorker(wp);
   }

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::mayConsume ( InputTag const & tag )

protectedinherited

Definition at line 82 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsume().

   {
     return collector_.mayConsume<T, BT>(tag);
   }

template<typename T , BranchType BT>

void art::ModuleBase::mayConsumeMany ( )

protectedinherited

Definition at line 96 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeMany().

   {
     collector_.mayConsumeMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::mayConsumeView ( InputTag const & )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::mayConsumeView ( InputTag const & tag )

inherited

Definition at line 89 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeView().

   {
     return collector_.mayConsumeView<T, BT>(tag);
   }

ModuleDescription const & art::ModuleBase::moduleDescription ( ) const

inherited

Definition at line 13 of file ModuleBase.cc.

References art::errors::LogicError.

Referenced by art::OutputModule::doRespondToOpenInputFile(), art::OutputModule::doWriteEvent(), art::Modifier::fillProductDescriptions(), art::OutputModule::makePlugins_(), art::OutputWorker::OutputWorker(), reco::shower::LArPandoraModularShowerCreation::produce(), art::Modifier::registerProducts(), and art::OutputModule::registerProducts().

   {
     if (md_.has_value()) {
       return *md_;
     }
 
     throw Exception{errors::LogicError,
                     "There was an error while calling moduleDescription().\n"}
       << "The moduleDescription() base-class member function cannot be called\n"
          "during module construction.  To determine which module is "
          "responsible\n"
          "for calling it, find the '<module type>:<module "
          "label>@Construction'\n"
          "tag in the message prefix above.  Please contact artists@fnal.gov\n"
          "for guidance.\n";
   }

G4InfoReducer& G4InfoReducer::operator= ( G4InfoReducer const & )

delete

G4InfoReducer& G4InfoReducer::operator= ( G4InfoReducer && )

delete

void G4InfoReducer::produce ( art::Event & e )

overridevirtual

Implements art::EDProducer.

Definition at line 104 of file G4InfoReducer_module.cc.

References DEFINE_ART_MODULE, sim::SimEnergyDepositLite::E(), fMinX, fMinY, fMinZ, fSedLabel, fUseOrigTrackID, fVoxelSizeX, fVoxelSizeY, fVoxelSizeZ, art::ProductRetriever::getValidHandle(), art::Event::put(), sim::SimEnergyDepositLite::T(), sim::SimEnergyDepositLite::TrackID(), x, sim::SimEnergyDepositLite::X(), y, sim::SimEnergyDepositLite::Y(), sim::SimEnergyDepositLite::Z(), and z.

 {
   // Implementation of required member function here.
 
   // Collect SimEnergyDeposit
   auto handle = e.getValidHandle<std::vector<sim::SimEnergyDeposit>>(fSedLabel);
   if (!handle.isValid()) {
     // throw some error
     std::cerr << "SimEnergyDeposit not found" << std::endl;
     throw std::exception();
   }
 
   struct comp {
     bool operator()(sim::SimEnergyDepositLite lhs, sim::SimEnergyDepositLite rhs) const
     {
       if (lhs.X() < rhs.X()) return true;
       if (lhs.X() > rhs.X()) return false;
       if (lhs.Y() < rhs.Y()) return true;
       if (lhs.Y() > rhs.Y()) return false;
       if (lhs.Z() < rhs.Z()) return true;
       if (lhs.Z() > rhs.Z()) return false;
       return lhs.TrackID() < rhs.TrackID();
     }
   };
   std::set<sim::SimEnergyDepositLite, comp> sedlite_v2;
   sim::SimEnergyDepositLite sed_lite;
   sim::SimEnergyDepositLite new_sed_lite;
 
   auto const& sed_v = *handle;
 
   /*double total_edep = 0.;
   for (size_t idx = 0; idx < sed_v.size(); ++idx) {
     total_edep += sed_v[idx].E();
   }
   std::cout << "total edep = " << total_edep << " count = " << sed_v.size() << std::endl;*/
 
   for (size_t idx = 0; idx < sed_v.size(); ++idx) {
     auto const& sed = sed_v[idx];
     // Voxelize coordinates to closest coordinate using fVoxelSize
     double x = sed.X() - std::fmod(sed.X() - fMinX, fVoxelSizeX) + fVoxelSizeX / 2.;
     double y = sed.Y() - std::fmod(sed.Y() - fMinY, fVoxelSizeY) + fVoxelSizeY / 2.;
     double z = sed.Z() - std::fmod(sed.Z() - fMinZ, fVoxelSizeZ) + fVoxelSizeZ / 2.;
 
     // Copy info to SimEnergyDepositLite
     if (fUseOrigTrackID) {
       sed_lite =
         sim::SimEnergyDepositLite(sed.E(), geo::Point_t(x, y, z), sed.T(), sed.OrigTrackID());
     }
     else {
       sed_lite = sim::SimEnergyDepositLite(sed.E(), geo::Point_t(x, y, z), sed.T(), sed.TrackID());
     }
     auto it = sedlite_v2.find(sed_lite);
     // Attempt to insert, if already exist we sum up energy deposit
 
     if (it != sedlite_v2.end()) {
       double new_energy = sed_lite.E() + it->E();
       double new_time = std::min(sed_lite.T(), it->T());
       sedlite_v2.erase(it);
       if (fUseOrigTrackID) {
         new_sed_lite =
           sim::SimEnergyDepositLite(new_energy, geo::Point_t(x, y, z), new_time, sed.OrigTrackID());
       }
       else {
         new_sed_lite =
           sim::SimEnergyDepositLite(new_energy, geo::Point_t(x, y, z), new_time, sed.TrackID());
       }
       sedlite_v2.insert(new_sed_lite);
     }
     else {
       sedlite_v2.insert(std::move(sed_lite));
     }
   }
 
   /*double new_total_edep = 0.;
   int counts = 0;
 
   for (auto it : sedlite_v2) {
     new_total_edep += it.E();
     counts += 1;
   }
   std::cout << "new total edep = " << new_total_edep << " with counts " << counts << " " << sedlite_v2.size() << std::endl;
   */
 
   // Create vector for SEDLite
   std::unique_ptr<std::vector<sim::SimEnergyDepositLite>> sedlite_v(
     new std::vector<sim::SimEnergyDepositLite>(sedlite_v2.begin(), sedlite_v2.end()));
 
   // Store SEDLite in event
   e.put(std::move(sedlite_v));
 }

void art::Modifier::registerProducts ( ProductDescriptions & productsToRegister )

inherited

Definition at line 16 of file Modifier.cc.

References art::ModuleBase::moduleDescription(), and art::ProductRegistryHelper::registerProducts().

   {
     ProductRegistryHelper::registerProducts(productsToRegister,
                                             moduleDescription());
   }

void art::ModuleBase::setModuleDescription ( ModuleDescription const & md )

inherited

Definition at line 31 of file ModuleBase.cc.

References art::ModuleBase::md_.

   {
     md_ = md;
   }

void art::ModuleBase::sortConsumables ( std::string const & current_process_name )

inherited

Definition at line 49 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::sortConsumables().

   {
     // Now that we know we have seen all the consumes declarations,
     // sort the results for fast lookup later.
     collector_.sortConsumables(current_process_name);
   }