Inheritance diagram for tutorial::ClusterCounter4:

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

	ClusterCounter4 (ClusterCounter4 const &)=delete

	ClusterCounter4 (ClusterCounter4 &&)=delete

ClusterCounter4 &	operator= (ClusterCounter4 const &)=delete

ClusterCounter4 &	operator= (ClusterCounter4 &&)=delete

void	analyze (art::Event const &e) override

void	beginJob () override

void	endJob () 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
float	sumAdc (const std::vector< art::Ptr< recob::Hit > > &hits) const

const simb::MCParticle *	getTruthParticle (detinfo::DetectorClocksData const &clockData, const std::vector< art::Ptr< recob::Hit > > &hits, float &fraction, bool &foundEmParent) const

Private Attributes
size_t	fEvNumber

TTree *	fEventTree

size_t	fNClusters

TTree *	fClusterTree

size_t	fNHits

float	fAdcSum

float	fClean

art::InputTag	fClusterModuleLabel

size_t	fMinSize

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 35 of file ClusterCounter4_module.cc.

Constructor & Destructor Documentation

tutorial::ClusterCounter4::ClusterCounter4 ( fhicl::ParameterSet const & p )

explicit

Definition at line 75 of file ClusterCounter4_module.cc.

                                                             : EDAnalyzer(p)
 {
     fClusterModuleLabel = p.get< std::string >("ClusterModuleLabel");
     fMinSize = p.get< size_t >("MinSize");
 }

tutorial::ClusterCounter4::ClusterCounter4 ( ClusterCounter4 const & )

delete

tutorial::ClusterCounter4::ClusterCounter4 ( ClusterCounter4 && )

delete

Member Function Documentation

void tutorial::ClusterCounter4::analyze ( art::Event const & e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 96 of file ClusterCounter4_module.cc.

 {
     fEvNumber = evt.id().event();
     mf::LogVerbatim("ClusterCounter4") << "ClusterCounter4 module on event #" << fEvNumber;
 
     // use auto to make the line shorter when you remember all art types,
     // the full type of the handle is: art::ValidHandle< std::vector<recob::Cluster> >
     auto cluHandle = evt.getValidHandle< std::vector<recob::Cluster> >(fClusterModuleLabel);
 
     // this will let us look for all hits associated to selected cluster
     art::FindManyP< recob::Hit > hitsFromClusters(cluHandle, evt, fClusterModuleLabel);
 
     fNClusters = 0;
 
     auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService>()->DataFor(evt);
     // and here we go by index:
     for (size_t i = 0; i < cluHandle->size(); ++i)
     {
         fNHits = cluHandle->at(i).NHits();
         fAdcSum = sumAdc(hitsFromClusters.at(i));
 
         bool isEM = false;
         const simb::MCParticle* p = getTruthParticle(clockData, hitsFromClusters.at(i), fClean, isEM);
         if (p)
         {
             if (isEM) { mf::LogVerbatim("ClusterCounter4") << "matched mother particle PDG: " << p->PdgCode(); }
             else { mf::LogVerbatim("ClusterCounter4") << "matched particle PDG: " << p->PdgCode(); }
         }
          else { mf::LogWarning("ClusterCounter4") << "No matcched particle??"; }
 
         fClusterTree->Fill();
 
         mf::LogVerbatim("ClusterCounter4")
             << "NHits() = " << fNHits << ", assn size = " << hitsFromClusters.at(i).size()
             << " SummedADC() = " << cluHandle->at(i).SummedADC() << ", sum hits adc = " << fAdcSum;
 
         if (fNHits >= fMinSize) { ++fNClusters; }
     }
     fEventTree->Fill();
 }

void tutorial::ClusterCounter4::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 81 of file ClusterCounter4_module.cc.

 {
     art::ServiceHandle<art::TFileService> tfs; // TTree's are created in the memory managed by ROOT (you don't delete them)
 
     fEventTree = tfs->make<TTree>("EventTree", "event by event info");
     fEventTree->Branch("event", &fEvNumber, "fEvNumber/I");
     fEventTree->Branch("nclusters", &fNClusters, "fNClusters/I");
 
     fClusterTree = tfs->make<TTree>("ClusterTree", "cluster by cluster info");
     fClusterTree->Branch("event", &fEvNumber, "fEvNumber/I");
     fClusterTree->Branch("nhits", &fNHits, "fNHits/I");
     fClusterTree->Branch("adcsum", &fAdcSum, "fAdcSum/F");
     fClusterTree->Branch("clean", &fClean, "fClean/F");
 }

void tutorial::ClusterCounter4::endJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 193 of file ClusterCounter4_module.cc.

 {
     mf::LogVerbatim("ClusterCounter4") << "ClusterCounter4 finished job";
 }

const simb::MCParticle * tutorial::ClusterCounter4::getTruthParticle	(	detinfo::DetectorClocksData const &	clockData,
		const std::vector< art::Ptr< recob::Hit > > &	hits,
		float &	fraction,
		bool &	foundEmParent
	)		const

private

Definition at line 137 of file ClusterCounter4_module.cc.

 {
     const simb::MCParticle* mcParticle = 0;
     fraction = 0;
     foundEmParent = false;
 
     art::ServiceHandle<cheat::BackTrackerService> bt_serv;
     art::ServiceHandle<cheat::ParticleInventoryService> pi_serv;
     std::unordered_map<int, double> trkIDE;
     for (auto const & h : hits)
     {
         for (auto const & ide : bt_serv->HitToTrackIDEs(clockData, h)) // loop over std::vector<sim::TrackIDE>
         {
             trkIDE[ide.trackID] += ide.energy; // sum energy contribution by each track ID
         }
     }
 
     int best_id = 0;
     double tot_e = 0, max_e = 0;
     for (auto const & contrib : trkIDE)
     {
         tot_e += contrib.second;     // sum total energy in these hits
         if (contrib.second > max_e)  // find track ID corresponding to max energy
         {
             max_e = contrib.second;
             best_id = contrib.first;
         }
     }
 
     if ((max_e > 0) && (tot_e > 0)) // ok, found something reasonable
     {
         if (best_id < 0)            // NOTE: negative ID means this is EM activity
         {                           // caused by track with the same but positive ID
             best_id = -best_id;     // --> we'll find mother MCParticle of these hits
             foundEmParent = true;
         }
         mcParticle = pi_serv->TrackIdToParticle_P(best_id); // MCParticle corresponding to track ID
         fraction = max_e / tot_e;
     }
     else { mf::LogWarning("ClusterCounter4") << "No energy deposits??"; }
 
     return mcParticle;
 }

ClusterCounter4& tutorial::ClusterCounter4::operator= ( ClusterCounter4 const & )

delete

ClusterCounter4& tutorial::ClusterCounter4::operator= ( ClusterCounter4 && )

delete

float tutorial::ClusterCounter4::sumAdc ( const std::vector< art::Ptr< recob::Hit > > & hits ) const

private

Definition at line 183 of file ClusterCounter4_module.cc.

 {
     float sum = 0;
     for (auto const & h : hits)
     {
         sum += h->SummedADC();
     }
     return sum;
 }

Member Data Documentation

float tutorial::ClusterCounter4::fAdcSum

private

Definition at line 67 of file ClusterCounter4_module.cc.

float tutorial::ClusterCounter4::fClean

private

Definition at line 68 of file ClusterCounter4_module.cc.

art::InputTag tutorial::ClusterCounter4::fClusterModuleLabel

private

Definition at line 71 of file ClusterCounter4_module.cc.

TTree* tutorial::ClusterCounter4::fClusterTree

private

Definition at line 65 of file ClusterCounter4_module.cc.

TTree* tutorial::ClusterCounter4::fEventTree

private

Definition at line 62 of file ClusterCounter4_module.cc.

size_t tutorial::ClusterCounter4::fEvNumber

private

Definition at line 60 of file ClusterCounter4_module.cc.

size_t tutorial::ClusterCounter4::fMinSize

private

Definition at line 72 of file ClusterCounter4_module.cc.

size_t tutorial::ClusterCounter4::fNClusters

private

Definition at line 63 of file ClusterCounter4_module.cc.

size_t tutorial::ClusterCounter4::fNHits

private

Definition at line 66 of file ClusterCounter4_module.cc.

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

protoduneana/protoduneana/TutorialExamples/ClusterCounter4_module.cc

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation