Inheritance diagram for cluster::HoughLineFinderAna:

Public Types
using	ModuleType = EDAnalyzer

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

	~HoughLineFinderAna ()

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)

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
std::string const &	processName () const

bool	wantAllEvents () const noexcept

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

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

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	analyze (const art::Event &)

void	beginJob ()

Private Attributes
std::string	fHoughModuleLabel

std::string	fDigitModuleLabel

std::string	fHitsModuleLabel

std::string	fDBScanModuleLabel

TTree *	ftree

int	fm_run

unsigned long int	fm_run_timestamp

int	fm_event

int	fm_plane

int	fm_dbsize

int	fm_clusterid

int	fm_wirespan

int	fm_sizeClusterZ

int	fm_sizeHitZ

float	fm_clusterslope

float	fm_clusterintercept

int *	fm_wireZ

int *	fm_hitidZ

float *	fm_mipZ

float *	fm_drifttimeZ

float *	fm_widthZ

float *	fm_upadcZ

Detailed Description

Definition at line 39 of file HoughLineFinderAna_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 22 of file EDAnalyzer.h.

Constructor & Destructor Documentation

cluster::HoughLineFinderAna::HoughLineFinderAna ( fhicl::ParameterSet const & pset )

explicit

Definition at line 79 of file HoughLineFinderAna_module.cc.

     : EDAnalyzer(pset)
     , fHoughModuleLabel(pset.get<std::string>("HoughModuleLabel"))
     , fDigitModuleLabel(pset.get<std::string>("DigitModuleLabel"))
     , fHitsModuleLabel(pset.get<std::string>("HitsModuleLabel"))
     , fDBScanModuleLabel(pset.get<std::string>("DBScanModuleLabel"))
     , fm_run(0)
     , fm_event(0)
     , fm_plane(0)
     , fm_dbsize(0)
     , fm_clusterid(0)
     , fm_wirespan(0)
     , fm_sizeClusterZ(10000)
     , fm_sizeHitZ(10000)
     , fm_clusterslope(0)
     , fm_clusterintercept(0)
   {}

cluster::HoughLineFinderAna::~HoughLineFinderAna ( )

Definition at line 98 of file HoughLineFinderAna_module.cc.

References fm_drifttimeZ, fm_hitidZ, fm_mipZ, fm_upadcZ, fm_widthZ, and fm_wireZ.

   {
     delete fm_hitidZ;
     delete fm_mipZ;
     delete fm_drifttimeZ;
     delete fm_widthZ;
     delete fm_upadcZ;
     delete fm_wireZ;
   }

Member Function Documentation

void cluster::HoughLineFinderAna::analyze ( const art::Event & evt )

private

Definition at line 139 of file HoughLineFinderAna_module.cc.

References DEFINE_ART_MODULE, art::EventID::event(), fDBScanModuleLabel, fHitsModuleLabel, fHoughModuleLabel, fm_clusterid, fm_clusterintercept, fm_clusterslope, fm_dbsize, fm_drifttimeZ, fm_event, fm_hitidZ, fm_mipZ, fm_plane, fm_run, fm_run_timestamp, fm_sizeClusterZ, fm_sizeHitZ, fm_upadcZ, fm_widthZ, fm_wirespan, fm_wireZ, ftree, art::ProductRetriever::getByLabel(), art::Event::id(), MF_LOG_VERBATIM, art::PtrVector< T >::push_back(), art::EventID::run(), art::PtrVector< T >::size(), art::Event::time(), art::Timestamp::value(), and geo::GeometryCore::Views().

   {
 
     art::Handle<std::vector<recob::Cluster>> hlfListHandle;
     evt.getByLabel(fHoughModuleLabel, hlfListHandle);
     art::Handle<std::vector<recob::Hit>> hitListHandle;
     evt.getByLabel(fHitsModuleLabel, hitListHandle);
     art::Handle<std::vector<recob::Cluster>> dbscanListHandle;
     evt.getByLabel(fDBScanModuleLabel, dbscanListHandle);
 
     art::FindManyP<recob::Hit> fmh(dbscanListHandle, evt, fDBScanModuleLabel);
     art::FindManyP<recob::Hit> fmhhl(hlfListHandle, evt, fHoughModuleLabel);
 
     art::PtrVector<recob::Cluster> clusters;
     art::PtrVector<recob::Cluster> dbclusters;
     //   art::PtrVector<recob::Hit> hits;// unused, as yet. EC, 5-Oct-2010.
 
     for (size_t ii = 0; ii < hlfListHandle->size(); ++ii) {
       art::Ptr<recob::Cluster> cluster(hlfListHandle, ii);
       clusters.push_back(cluster);
     }
 
     for (size_t ii = 0; ii < dbscanListHandle->size(); ++ii) {
       art::Ptr<recob::Cluster> dbcluster(dbscanListHandle, ii);
       dbclusters.push_back(dbcluster);
     }
 
     MF_LOG_VERBATIM("HoughLineFinderAna") << "run    : " << evt.id().run();
     //std::cout << "subrun : " << evt.subRun() << std::endl;
     MF_LOG_VERBATIM("HoughLineFinderAna") << "event  : " << evt.id().event();
     fm_run = evt.id().run();
     fm_event = evt.id().event();
     fm_run_timestamp = evt.time().value(); // won't cast, EC, 7-Oct-2010.
     unsigned int firstwire = 0;
     unsigned int lastwire = 0;
     fm_sizeClusterZ = 0;
     fm_sizeHitZ = 0;
     fm_dbsize = 0;
     art::ServiceHandle<geo::Geometry const> geo;
 
     for (auto view : geo->Views()) {
 
       fm_dbsize = 0;
       fm_sizeClusterZ = clusters.size();
 
       for (size_t j = 0; j < dbclusters.size(); ++j) {
         if (dbclusters[j]->View() == view) {
           std::vector<art::Ptr<recob::Hit>> _dbhits = fmh.at(j);
           fm_dbsize += _dbhits.size();
           if (_dbhits.size() > 0) fm_plane = _dbhits.at(0)->WireID().Plane;
         }
       }
 
       for (size_t j = 0; j < clusters.size(); ++j) {
         if (clusters[j]->View() == view) {
           fm_clusterid = clusters[j]->ID();
           std::vector<art::Ptr<recob::Hit>> _hits = fmhhl.at(j);
           fm_clusterslope = (double)std::tan(clusters[j]->StartAngle());
           fm_clusterintercept = (double)clusters[j]->StartTick();
           if (_hits.size() != 0) {
             fm_plane = _hits.at(0)->WireID().Plane;
             firstwire = _hits[0]->WireID().Wire;
             lastwire = _hits[_hits.size() - 1]->WireID().Wire;
             fm_wirespan = lastwire - firstwire;
             fm_sizeHitZ = _hits.size();
 
             for (unsigned int i = 0; i < _hits.size(); ++i) {
 
               fm_hitidZ[i] = i;
               fm_wireZ[i] = _hits[i]->WireID().Wire;
               fm_mipZ[i] = (double)_hits[i]->Integral();
               fm_drifttimeZ[i] = (double)_hits[i]->PeakTime();
               fm_widthZ[i] = (double)(2. * _hits[i]->RMS());
               fm_upadcZ[i] = (double)_hits[i]->Integral();
             }
 
             ftree->Fill();
           }
         } //end if in the correct view
       }   // end loop over clusters
     }     // end loop over views
   }

void cluster::HoughLineFinderAna::beginJob ( )

privatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 109 of file HoughLineFinderAna_module.cc.

References fm_clusterid, fm_clusterintercept, fm_clusterslope, fm_dbsize, fm_drifttimeZ, fm_event, fm_hitidZ, fm_mipZ, fm_plane, fm_run, fm_run_timestamp, fm_sizeClusterZ, fm_sizeHitZ, fm_upadcZ, fm_widthZ, fm_wirespan, fm_wireZ, and ftree.

   {
 
     // get access to the TFile service
     art::ServiceHandle<art::TFileService const> tfs;
     ftree = tfs->make<TTree>("HoughTree", "HoughTree");
     fm_hitidZ = new int[fm_sizeHitZ];
     fm_mipZ = new float[fm_sizeHitZ];
     fm_drifttimeZ = new float[fm_sizeHitZ];
     fm_widthZ = new float[fm_sizeHitZ];
     fm_upadcZ = new float[fm_sizeHitZ];
     fm_wireZ = new int[fm_sizeHitZ];
     ftree->Branch("run", &fm_run, "run/I");
     ftree->Branch("run_timestamp", &fm_run_timestamp, "run_timestamp/l"); //l is for ULong64_t
     ftree->Branch("event", &fm_event, "event/I");
     ftree->Branch("plane", &fm_plane, "plane/I");
     ftree->Branch("dbsize", &fm_dbsize, "dbsize/I");
     ftree->Branch("clusterid", &fm_clusterid, "clusterid/I");
     ftree->Branch("clusterslope", &fm_clusterslope, "clusterslope/F");
     ftree->Branch("clusterintercept", &fm_clusterintercept, "clusterintecept/F");
     ftree->Branch("wirespan", &fm_wirespan, "wirespan/I");
     ftree->Branch("numberHits", &fm_sizeHitZ, "numberHits/I");
     ftree->Branch("numberClusters", &fm_sizeClusterZ, "numberClusters/I");
     ftree->Branch("hitidZ", fm_hitidZ, "hitidZ[numberHits]/I");
     ftree->Branch("wireZ", fm_wireZ, "wireZ[numberHits]/I");
     ftree->Branch("mipZ", fm_mipZ, "mipZ[numberHits]/F");
     ftree->Branch("drifttimeZ", fm_drifttimeZ, "drifttitmeZ[numberHits]/F");
     ftree->Branch("widthZ", fm_widthZ, "widthZ[numberHits]/F");
   }

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::Analyzer::doBeginJob ( SharedResources const & resources )

inherited

Definition at line 25 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

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

inherited

Definition at line 68 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     beginRunWithFrame(std::as_const(rp).makeRun(mc), frame);
     return true;
   }

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

inherited

Definition at line 84 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     beginSubRunWithFrame(std::as_const(srp).makeSubRun(mc), frame);
     return true;
   }

void art::detail::Analyzer::doEndJob ( )

inherited

Definition at line 33 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

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

inherited

Definition at line 76 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     endRunWithFrame(std::as_const(rp).makeRun(mc), frame);
     return true;
   }

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

inherited

Definition at line 92 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     endSubRunWithFrame(std::as_const(srp).makeSubRun(mc), frame);
     return true;
   }

bool art::detail::Analyzer::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 100 of file Analyzer.cc.

References e, and art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     auto const e = std::as_const(ep).makeEvent(mc);
     if (wantEvent(mc.scheduleID(), e)) {
       ++counts_run;
       ProcessingFrame const frame{mc.scheduleID()};
       analyzeWithFrame(e, frame);
       ++counts_passed;
     }
     return true;
   }

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

inherited

Definition at line 47 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

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

inherited

Definition at line 61 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

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

inherited

Definition at line 40 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

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

inherited

Definition at line 54 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

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

Handle< TriggerResults > art::Observer::getTriggerResults ( Event const & e ) const

protectedinherited

Definition at line 75 of file Observer.cc.

References art::ProductRetriever::get(), and art::Observer::selectors_.

Referenced by art::OutputModule::doWriteEvent(), and art::Observer::wantAllEvents().

   {
     if (selectors_) {
       return selectors_->getOneTriggerResults(e);
     }
 
     // The following applies for cases where no SelectEvents entries
     // exist.
     Handle<TriggerResults> h;
     if (e.get(empty_process_name, h)) {
       return h;
     }
     return Handle<TriggerResults>{};
   }

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

string const & art::Observer::processName ( ) const

protectedinherited

Definition at line 57 of file Observer.cc.

References art::Observer::process_name_.

Referenced by art::FileDumperOutput::printPrincipal().

   {
     return process_name_;
   }

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

bool art::Observer::wantAllEvents ( ) const

inlineprotectednoexceptinherited

Definition at line 31 of file Observer.h.

References e, art::Observer::getTriggerResults(), art::Observer::wantAllEvents_, and art::Observer::wantEvent().

     {
       return wantAllEvents_;
     }

bool art::Observer::wantEvent	(	ScheduleID	id,
		Event const &	e
	)		const

protectedinherited

Definition at line 63 of file Observer.cc.

References art::Observer::rejectors_, art::Observer::selectors_, and art::Observer::wantAllEvents_.

Referenced by art::OutputModule::doEvent(), art::OutputModule::doWriteEvent(), and art::Observer::wantAllEvents().

   {
     if (wantAllEvents_) {
       return true;
     }
     bool const select_event = selectors_ ? selectors_->matchEvent(id, e) : true;
     bool const reject_event =
       rejectors_ ? rejectors_->matchEvent(id, e) : false;
     return select_event and not reject_event;
   }