Inheritance diagram for cluster::HoughLineFinderAna:

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

	~HoughLineFinderAna ()

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

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

using	ModuleType = EDAnalyzer

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

Definition at line 39 of file HoughLineFinderAna_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 84 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 104 of file HoughLineFinderAna_module.cc.

   {
     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 147 of file HoughLineFinderAna_module.cc.

   {
 
     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 115 of file HoughLineFinderAna_module.cc.

   {
 
 
     // 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");
   }