Inheritance diagram for shwf::ShowerCheater:

Public Types
using	ModuleType = EDProducer

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

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

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 Member Functions
void	produce (art::Event &evt) override

Private Attributes
std::string const	fCheatedClusterLabel

std::string const	fG4ModuleLabel
	label for module running G4 and making particles, etc More...

Detailed Description

Definition at line 33 of file ShowerCheater_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

shwf::ShowerCheater::ShowerCheater ( fhicl::ParameterSet const & pset )

explicit

Definition at line 49 of file ShowerCheater_module.cc.

References fCheatedClusterLabel, and fG4ModuleLabel.

     : EDProducer{pset}
     , fCheatedClusterLabel{pset.get<std::string>("CheatedClusterLabel", "cluster")}
     , fG4ModuleLabel{pset.get<std::string>("G4ModuleLabel", "largeant")}
   {
     produces<std::vector<recob::Shower>>();
     produces<std::vector<recob::SpacePoint>>();
     produces<art::Assns<recob::Shower, recob::Cluster>>();
     produces<art::Assns<recob::Shower, recob::SpacePoint>>();
     produces<art::Assns<recob::Shower, recob::Hit>>();
     produces<art::Assns<recob::Hit, recob::SpacePoint>>();
   }

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

void shwf::ShowerCheater::produce ( art::Event & evt )

overrideprivatevirtual

Todo:: figure out the max transverse width of the shower in the x and y directions

Implements art::EDProducer.

Definition at line 63 of file ShowerCheater_module.cc.

References util::abs(), util::CreateAssn(), DEFINE_ART_MODULE, e, fCheatedClusterLabel, art::fill_ptr_vector(), art::ProductRetriever::getByLabel(), hits(), cheat::BackTrackerService::HitToXYZ(), cheat::ParticleInventoryService::ParticleList(), art::Event::put(), recob::Shower::set_direction(), recob::Shower::set_direction_err(), and recob::Shower::set_id().

   {
     art::ServiceHandle<cheat::BackTrackerService const> bt_serv;
     art::ServiceHandle<cheat::ParticleInventoryService const> pi_serv;
     art::ServiceHandle<geo::Geometry const> geo;
 
     // grab the clusters that have been reconstructed
     art::Handle<std::vector<recob::Cluster>> clustercol;
     evt.getByLabel(fCheatedClusterLabel, clustercol);
 
     art::FindManyP<recob::Hit> fmh(clustercol, evt, fCheatedClusterLabel);
 
     // make a vector of them - we aren't writing anything out to a file
     // so no need for a art::PtrVector here
     std::vector<art::Ptr<recob::Cluster>> clusters;
     art::fill_ptr_vector(clusters, clustercol);
 
     // make a map of vectors of art::Ptrs keyed by eveID values
     std::map<int, std::vector<std::pair<size_t, art::Ptr<recob::Cluster>>>> eveClusterMap;
 
     // loop over all clusters and fill in the map
     for (size_t c = 0; c < clusters.size(); ++c) {
 
       // in the ClusterCheater module we set the cluster ID to be
       // the eve particle track ID*1000 + plane*100 + tpc*10 + cryostat number.  The
       // floor function on the cluster ID / 1000 will give us
       // the eve track ID
       int eveID = floor(clusters[c]->ID() / 1000.);
 
       std::pair<size_t, art::Ptr<recob::Cluster>> indexPtr(c, clusters[c]);
 
       eveClusterMap[eveID].push_back(indexPtr);
 
     } // end loop over clusters
 
     // loop over the map and make prongs
     std::unique_ptr<std::vector<recob::Shower>> showercol(new std::vector<recob::Shower>);
     std::unique_ptr<std::vector<recob::SpacePoint>> spcol(new std::vector<recob::SpacePoint>);
     std::unique_ptr<art::Assns<recob::Shower, recob::Cluster>> scassn(
       new art::Assns<recob::Shower, recob::Cluster>);
     std::unique_ptr<art::Assns<recob::Shower, recob::Hit>> shassn(
       new art::Assns<recob::Shower, recob::Hit>);
     std::unique_ptr<art::Assns<recob::Shower, recob::SpacePoint>> sspassn(
       new art::Assns<recob::Shower, recob::SpacePoint>);
     std::unique_ptr<art::Assns<recob::Hit, recob::SpacePoint>> sphassn(
       new art::Assns<recob::Hit, recob::SpacePoint>);
 
     auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
 
     for (auto const& clusterMapItr : eveClusterMap) {
 
       // separate out the hits for each particle into the different views
       std::vector<std::pair<size_t, art::Ptr<recob::Cluster>>> const& eveClusters =
         clusterMapItr.second;
 
       size_t startSPIndx = spcol->size();
 
       std::vector<art::Ptr<recob::Cluster>> ptrvs;
       std::vector<size_t> idxs;
 
       for (auto const& idxPtr : eveClusters) {
         idxs.push_back(idxPtr.first);
         ptrvs.push_back(idxPtr.second);
 
         // need to make the space points for this prong
         // loop over the hits for this cluster and make
         // a space point for each one
         // set the SpacePoint ID to be the cluster ID*10000
         // + the hit index in the cluster PtrVector of hits
         std::vector<art::Ptr<recob::Hit>> const& hits = fmh.at(idxPtr.first);
 
         for (size_t h = 0; h < hits.size(); ++h) {
           art::Ptr<recob::Hit> hit = hits[h];
           // add up the charge from the hits on the collection plane
           std::vector<double> xyz = bt_serv->HitToXYZ(clockData, hit);
           double sperr[6] = {0.01, 0.01, 0.1, 0.001, 0.001, 0.001};
 
           // make the space point and set its ID and XYZ
           // the errors and chi^2 are set to "good" values as we know the information perfectly
           recob::SpacePoint sp(&xyz[0], sperr, 0.9, idxPtr.second->ID() * 10000 + h);
           spcol->push_back(sp);
 
           // associate the space point to the hit
           util::CreateAssn(evt, *spcol, hit, *sphassn);
 
         } // end loop over hits
       }   // end loop over pairs of index values and cluster Ptrs
 
       size_t endSPIndx = spcol->size();
 
       // is this prong electro-magnetic in nature or
       // hadronic/muonic?  EM --> shower, everything else is a track
       if (std::abs(pi_serv->ParticleList()[clusterMapItr.first]->PdgCode()) == 11 ||
           std::abs(pi_serv->ParticleList()[clusterMapItr.first]->PdgCode()) == 22 ||
           std::abs(pi_serv->ParticleList()[clusterMapItr.first]->PdgCode()) == 111) {
 
         mf::LogInfo("ShowerCheater")
           << "prong of " << clusterMapItr.first << " is a shower with pdg code "
           << pi_serv->ParticleList()[clusterMapItr.first]->PdgCode();
 
         // get the direction cosine for the eve ID particle
         // just use the same for both the start and end of the prong
         const TLorentzVector initmom = pi_serv->ParticleList()[clusterMapItr.first]->Momentum();
         TVector3 dcos(
           initmom.Px() / initmom.Mag(), initmom.Py() / initmom.Mag(), initmom.Pz() / initmom.Mag());
         TVector3 dcosErr(1.e-3, 1.e-3, 1.e-3);
 
         //double maxTransWidth[2] = { util::kBogusD, util::kBogusD };
         //double distanceMaxWidth = util::kBogusD;
 
         // add a prong to the collection.  Make the prong
         // ID be the same as the track ID for the eve particle
         recob::Shower s;
         s.set_id(showercol->size());
         s.set_direction(dcos);
         s.set_direction_err(dcosErr);
         showercol->push_back(s);
         // associate the shower with its clusters
         util::CreateAssn(evt, *showercol, ptrvs, *scassn);
 
         // get the hits associated with each cluster and associate those with the shower
         for (size_t i = 0; i < idxs.size(); ++i) {
           std::vector<art::Ptr<recob::Hit>> hits = fmh.at(i);
           util::CreateAssn(evt, *showercol, hits, *shassn);
         }
 
         // associate the shower with the space points
         util::CreateAssn(evt, *showercol, *spcol, *sspassn, startSPIndx, endSPIndx);
 
         mf::LogInfo("ShowerCheater") << "adding shower: \n"
                                      << showercol->back() << "\nto collection.";
 
       } // end if this is a shower
     }   // end loop over the map
 
     evt.put(std::move(showercol));
     evt.put(std::move(spcol));
     evt.put(std::move(scassn));
     evt.put(std::move(shassn));
     evt.put(std::move(sspassn));
     evt.put(std::move(sphassn));
   } // end produce

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

Member Data Documentation

std::string const shwf::ShowerCheater::fCheatedClusterLabel

private

label for module creating recob::Cluster objects

Definition at line 40 of file ShowerCheater_module.cc.

Referenced by produce(), and ShowerCheater().

std::string const shwf::ShowerCheater::fG4ModuleLabel

private

label for module running G4 and making particles, etc

Definition at line 42 of file ShowerCheater_module.cc.

Referenced by ShowerCheater().

The documentation for this class was generated from the following file:

larreco/v09_25_00/source/larreco/ShowerFinder/ShowerCheater_module.cc

Public Types

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation