Prints all the hits associated to the specified shower. More...

Inheritance diagram for lar::test::AssnsChainTest:

Classes
struct	Config

Public Types
using	Parameters = art::EDAnalyzer::Table< Config >

Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >

using	ModuleType = EDAnalyzer

Public Member Functions
	AssnsChainTest (Parameters const &config)

virtual void	analyze (art::Event const &event) override

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

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

std::string	workerType () const

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

	Analyzer (fhicl::ParameterSet const &pset)

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

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)

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

	Observer (Observer const &)=delete

	Observer (Observer &&)=delete

Observer &	operator= (Observer const &)=delete

Observer &	operator= (Observer &&)=delete

void	registerProducts (ProductDescriptions &, ModuleDescription const &)

void	fillDescriptions (ModuleDescription const &)

fhicl::ParameterSetID	selectorConfig () const

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)

Private Member Functions
template<typename ShowerHandle >
void	printAssociatedPFOs (art::Event const &event, ShowerHandle const &showers) const

template<typename ShowerHandle >
void	printAssociatedClusters (art::Event const &event, ShowerHandle const &showers) const

template<typename ShowerHandle >
void	printAssociatedHits (art::Event const &event, ShowerHandle const &showers) const

Private Attributes
art::InputTag	showerTag
	Label of the input collection of showers. More...

unsigned int	nObjectsPerLine
	Number of objects to print on one line. More...

unsigned int	nShowers
	Total number of expected showers. More...

unsigned int	nPFOs
	Total number of expected particles. More...

unsigned int	nClusters
	Total number of expected clusters. More...

unsigned int	nHits
	Total number of expected hits. More...

Additional Inherited Members
Protected Member Functions inherited from art::Observer
std::string const &	processName () const

bool	wantAllEvents () const noexcept

bool	wantEvent (ScheduleID id, Event const &e) const

Handle< TriggerResults >	getTriggerResults (Event const &e) const

	Observer (fhicl::ParameterSet const &config)

	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)

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

Detailed Description

Prints all the hits associated to the specified shower.

The hits are searched transversing associations from showers to particle flow objects to clusters to hits.

Configuration parameters

showers (input tag, required): the shower collection
hitsPerLine (positive integer, default: 10): how many associated hits to print per line
nShowers _(positive integer, mandatory)_: total number of expected showers
nParticles _(positive integer, mandatory)_: total number of expected particle flow objects
nClusters _(positive integer, mandatory)_: total number of expected clusters
nHits _(positive integer, mandatory)_: total number of expected hits

Definition at line 65 of file AssnsChainTest_module.cc.

Member Typedef Documentation

using lar::test::AssnsChainTest::Parameters = art::EDAnalyzer::Table<Config>

Definition at line 105 of file AssnsChainTest_module.cc.

Constructor & Destructor Documentation

lar::test::AssnsChainTest::AssnsChainTest ( Parameters const & config )

explicit

Definition at line 141 of file AssnsChainTest_module.cc.

   : art::EDAnalyzer(config)
   , showerTag(config().showers())
   , nObjectsPerLine(config().hitsPerLine())
   , nShowers(config().nShowers())
   , nPFOs(config().nParticles())
   , nClusters(config().nClusters())
   , nHits(config().nHits())
   {}

Member Function Documentation

void lar::test::AssnsChainTest::analyze ( art::Event const & event )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 153 of file AssnsChainTest_module.cc.

                                                          {
 
   //
   // read the input collection
   //
   auto showers = event.getValidHandle<std::vector<recob::Shower>>(showerTag);
   mf::LogVerbatim("AssnsChainTest") << event.id() << " contains "
     << showers->size() << " showers from '" << showerTag.encode() << "'";
 
   if (showers->size() != nShowers) {
     throw cet::exception("AssnsChainTest")
       << "Data product '" << showerTag.encode() << "' contains "
       << showers->size() << " showers, " << nShowers << " were expected.\n";
   }
 
   mf::LogVerbatim("AssnsChainTest") << "\nPrinting: shower particle";
   printAssociatedPFOs(event, showers);
 
   mf::LogVerbatim("AssnsChainTest") << "\nPrinting: shower clusters";
   printAssociatedClusters(event, showers);
 
   mf::LogVerbatim("AssnsChainTest") << "\nPrinting: shower hits";
   printAssociatedHits(event, showers);
 
 } // lar::test::AssnsChainTest::analyze()

template<typename ShowerHandle >

void lar::test::AssnsChainTest::printAssociatedClusters	(	art::Event const &	event,
		ShowerHandle const &	showers
	)		const

private

Definition at line 261 of file AssnsChainTest_module.cc.

 {
   std::string const objectsDesc = "clusters";
 
   //
   // get the associated objects
   //
   lar::FindManyInChainP<recob::Cluster, recob::PFParticle> showerObjects
     (showers, event, showerTag);
   assert(showerObjects.size() == showers->size());
 
   using ObjectPtr_t = decltype(showerObjects)::TargetPtr_t;
 
   //
   // print the associated objects (just the art pointer so far)
   //
   unsigned int nDuplicates = 0;
   std::set<art::Ptr<recob::Cluster>> foundObjects; // all objects found
   std::set<art::ProductID> foundObjectProducts;
 
   unsigned int const pageSize = nObjectsPerLine; // objects per line
   mf::LogVerbatim log("AssnsChainTest");
   for (std::size_t iShower = 0; iShower < showers->size(); ++iShower) {
     auto const& objects = showerObjects.at(iShower);
     log << "\n #" << iShower << ": " << objects.size() << " " << objectsDesc;
     if (objects.empty()) continue;
 
     unsigned int objectsLeft = pageSize;
     for (ObjectPtr_t const& object: objects) {
       if (foundObjects.count(object) == 0) {
         foundObjects.insert(object);
         foundObjectProducts.insert(object.id());
       }
       else {
         mf::LogProblem("AssnsChainTest")
           << "ERROR: Cluster " << object << " appears in more than one shower!";
         ++nDuplicates;
       }
       if (objectsLeft-- == 0) {
         objectsLeft = pageSize;
         log << "\n  ";
       }
       log << " " << object; // just print the art pointer
     } // for objects
 
   } // for iShower
 
   mf::LogVerbatim("AssnsChainTest")
     << foundObjects.size() << " " << objectsDesc << " collected for "
     << showers->size() << " showers ('"
     << showers.provenance()->inputTag().encode() << "') from "
     << foundObjectProducts.size() << " data products:";
   for (art::ProductID const& PID: foundObjectProducts) {
     art::Handle<std::vector<recob::Cluster>> objects;
     if (event.get(PID, objects)) {
       mf::LogVerbatim("AssnsChainTest") << " - '"
         << objects.provenance()->inputTag().encode() << "' (contains "
         << objects->size() << " " << objectsDesc << ")";
     }
     else {
       mf::LogVerbatim("AssnsChainTest") << " - <" << PID << "> (not found!)";
     }
   } // for PIDs
 
   if (nDuplicates > 0) {
     throw cet::exception("AssnsChainTest")
       << "Test failed: " << nDuplicates
       << " clusters appear in more than one shower.\n";
   }
   if (foundObjects.size() != nClusters) {
     throw cet::exception("AssnsChainTest")
       << "Test failed: counted " << foundObjects.size()
       << " clusters, expected " << nClusters << ".\n";
   }
 
 } // lar::test::AssnsChainTest::printAssociatedClusters()

template<typename ShowerHandle >

void lar::test::AssnsChainTest::printAssociatedHits	(	art::Event const &	event,
		ShowerHandle const &	showers
	)		const

private

Definition at line 183 of file AssnsChainTest_module.cc.

 {
 
   //
   // get the associated hits
   //
   lar::FindManyInChainP<recob::Hit, recob::Cluster, recob::PFParticle>
     showerHits(showers, event, showerTag);
   assert(showerHits.size() == showers->size());
 
   //
   // print the associated hits (just the art pointer so far)
   //
   unsigned int nDuplicates = 0;
   std::set<art::Ptr<recob::Hit>> foundHits; // all hits found
   std::set<art::ProductID> foundHitProducts;
 
   unsigned int const pageSize = nObjectsPerLine; // hits per line
   mf::LogVerbatim log("AssnsChainTest");
   for (std::size_t iShower = 0; iShower < showers->size(); ++iShower) {
     auto const& hits = showerHits.at(iShower);
     log << "\n #" << iShower << ": " << hits.size() << " hits";
     if (hits.empty()) continue;
 
     unsigned int hitsLeft = pageSize;
     for (art::Ptr<recob::Hit> const& hit: hits) {
       if (foundHits.count(hit) == 0) {
         foundHits.insert(hit);
         foundHitProducts.insert(hit.id());
       }
       else {
         mf::LogProblem("AssnsChainTest")
           << "ERROR: Hit " << hit << " appears in more than one shower!";
         ++nDuplicates;
       }
       if (hitsLeft-- == 0) {
         hitsLeft = pageSize;
         log << "\n  ";
       }
       log << " " << hit; // just print the art pointer
     } // for hits
 
   } // for iShower
 
   mf::LogVerbatim("AssnsChainTest")
     << foundHits.size() << " hits collected for "
     << showers->size() << " showers ('"
     << showers.provenance()->inputTag().encode() << "') from "
     << foundHitProducts.size() << " data products:";
   for (art::ProductID const& PID: foundHitProducts) {
     art::Handle<std::vector<recob::Hit>> hits;
     if (event.get(PID, hits)) {
       mf::LogVerbatim("AssnsChainTest") << " - '"
         << hits.provenance()->inputTag().encode() << "' (contains "
         << hits->size() << " hits)";
     }
     else {
       mf::LogVerbatim("AssnsChainTest") << " - <" << PID << "> (not found!)";
     }
   } // for PIDs
 
   if (nDuplicates > 0) {
     throw cet::exception("AssnsChainTest")
       << "Test failed: " << nDuplicates
       << " hits appear in more than one shower.\n";
   }
   if (foundHits.size() != nHits) {
     throw cet::exception("AssnsChainTest")
       << "Test failed: counted " << foundHits.size()
       << " hits, expected " << nHits << ".\n";
   }
 
 } // lar::test::AssnsChainTest::printAssociatedHits()

template<typename ShowerHandle >

void lar::test::AssnsChainTest::printAssociatedPFOs	(	art::Event const &	event,
		ShowerHandle const &	showers
	)		const

private

Definition at line 342 of file AssnsChainTest_module.cc.

 {
   std::string const objectsDesc = "particle flow objects";
 
   //
   // get the associated objects
   //
   lar::FindManyInChainP<recob::PFParticle> showerObjects
     (showers, event, showerTag);
   assert(showerObjects.size() == showers->size());
 
   using ObjectPtr_t = decltype(showerObjects)::TargetPtr_t;
 
   //
   // print the associated objects (just the art pointer so far)
   //
   unsigned int nDuplicates = 0;
   std::set<art::Ptr<recob::PFParticle>> foundObjects; // all objects found
   std::set<art::ProductID> foundObjectProducts;
 
   unsigned int const pageSize = nObjectsPerLine; // objects per line
   mf::LogVerbatim log("AssnsChainTest");
   for (std::size_t iShower = 0; iShower < showers->size(); ++iShower) {
     auto const& objects = showerObjects.at(iShower);
     log << "\n #" << iShower << ": " << objects.size() << " " << objectsDesc;
 
     unsigned int objectsLeft = pageSize;
     for (ObjectPtr_t const& object: objects) {
       if (foundObjects.count(object) == 0) {
         foundObjects.insert(object);
         foundObjectProducts.insert(object.id());
       }
       else {
         mf::LogProblem("AssnsChainTest")
           << "ERROR: Particle " << object
           << " appears in more than one shower!";
         ++nDuplicates;
       }
       if (objectsLeft-- == 0) {
         objectsLeft = pageSize;
         log << "\n  ";
       }
       log << " " << object; // just print the art pointer
     } // for objects
 
     if (objects.size() != 1) {
       throw cet::exception("AssnsChainTest")
         << "all showers are expected to have one PFParticle associated, while #"
         << iShower << " has " << objects.size() << "\n";
     }
 
   } // for iShower
 
   mf::LogVerbatim("AssnsChainTest")
     << foundObjects.size() << " " << objectsDesc << " collected for "
     << showers->size() << " showers ('"
     << showers.provenance()->inputTag().encode() << "') from "
     << foundObjectProducts.size() << " data products:";
   for (art::ProductID const& PID: foundObjectProducts) {
     art::Handle<std::vector<recob::PFParticle>> objects;
     if (event.get(PID, objects)) {
       mf::LogVerbatim("AssnsChainTest") << " - '"
         << objects.provenance()->inputTag().encode() << "' (contains "
         << objects->size() << " " << objectsDesc << ")";
     }
     else {
       mf::LogVerbatim("AssnsChainTest") << " - <" << PID << "> (not found!)";
     }
   } // for PIDs
 
   if (nDuplicates > 0) {
     throw cet::exception("AssnsChainTest")
       << "Test failed: " << nDuplicates
       << " particle flow objects appear in more than one shower.\n";
   }
   if (foundObjects.size() != nPFOs) {
     throw cet::exception("AssnsChainTest")
       << "Test failed: counted " << foundObjects.size()
       << " particle flow objects, expected " << nPFOs << ".\n";
   }
 
 } // lar::test::AssnsChainTest::printAssociatedPFOs()

Member Data Documentation

unsigned int lar::test::AssnsChainTest::nClusters

private

Total number of expected clusters.

Definition at line 116 of file AssnsChainTest_module.cc.

unsigned int lar::test::AssnsChainTest::nHits

private

Total number of expected hits.

Definition at line 117 of file AssnsChainTest_module.cc.

unsigned int lar::test::AssnsChainTest::nObjectsPerLine

private

Number of objects to print on one line.

Definition at line 113 of file AssnsChainTest_module.cc.

unsigned int lar::test::AssnsChainTest::nPFOs

private

Total number of expected particles.

Definition at line 115 of file AssnsChainTest_module.cc.

unsigned int lar::test::AssnsChainTest::nShowers

private

Total number of expected showers.

Definition at line 114 of file AssnsChainTest_module.cc.

art::InputTag lar::test::AssnsChainTest::showerTag

private

Label of the input collection of showers.

Definition at line 112 of file AssnsChainTest_module.cc.

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

lardata/test/Utilities/testAssnsChainUtils/AssnsChainTest_module.cc

Classes

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Configuration parameters

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation