Inheritance diagram for tutorial::ClusterCounter:

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
	ClusterCounter (Parameters const &config)

	ClusterCounter (ClusterCounter const &)=delete

	ClusterCounter (ClusterCounter &&)=delete

ClusterCounter &	operator= (ClusterCounter const &)=delete

ClusterCounter &	operator= (ClusterCounter &&)=delete

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

void	beginJob () override

void	endJob () override

	ClusterCounter (fhicl::ParameterSet const &p)

	ClusterCounter (ClusterCounter const &)=delete

	ClusterCounter (ClusterCounter &&)=delete

ClusterCounter &	operator= (ClusterCounter const &)=delete

ClusterCounter &	operator= (ClusterCounter &&)=delete

void	analyze (art::Event const &e) 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
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 34 of file ClusterCounter5_module.cc.

Member Typedef Documentation

using tutorial::ClusterCounter::Parameters = art::EDAnalyzer::Table<Config>

Definition at line 51 of file ClusterCounter5_module.cc.

Constructor & Destructor Documentation

tutorial::ClusterCounter::ClusterCounter ( Parameters const & config )

explicit

Definition at line 91 of file ClusterCounter5_module.cc.

                                                        : art::EDAnalyzer(config),
     fClusterModuleLabel(config().ClusterModuleLabel()),
         fMinSize(config().MinSize())
 {
 }

tutorial::ClusterCounter::ClusterCounter ( ClusterCounter const & )

delete

tutorial::ClusterCounter::ClusterCounter ( ClusterCounter && )

delete

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

explicit

Definition at line 37 of file ClusterCounter_module.cc.

37 : EDAnalyzer(p) {}

art::EDAnalyzer::EDAnalyzer

EDAnalyzer(fhicl::ParameterSet const &pset)

Definition: EDAnalyzer.h:25

test.p

p

Definition: test.py:223

tutorial::ClusterCounter::ClusterCounter ( ClusterCounter const & )

delete

tutorial::ClusterCounter::ClusterCounter ( ClusterCounter && )

delete

Member Function Documentation

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

overridevirtual

Implements art::EDAnalyzer.

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

overridevirtual

Implements art::EDAnalyzer.

Definition at line 112 of file ClusterCounter5_module.cc.

 {
     fEvNumber = evt.id().event();
     mf::LogVerbatim("ClusterCounter") << "ClusterCounter 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("ClusterCounter") << "matched mother particle PDG: " << p->PdgCode(); }
             else { mf::LogVerbatim("ClusterCounter") << "matched particle PDG: " << p->PdgCode(); }
         }
         else { mf::LogWarning("ClusterCounter") << "No matcched particle??"; }
 
         fClusterTree->Fill();
 
         mf::LogVerbatim("ClusterCounter")
             << "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::ClusterCounter::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 97 of file ClusterCounter5_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::ClusterCounter::endJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 209 of file ClusterCounter5_module.cc.

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

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

private

Definition at line 153 of file ClusterCounter5_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("ClusterCounter") << "No energy deposits??"; }
 
     return mcParticle;
 }

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

delete

ClusterCounter& tutorial::ClusterCounter::operator= ( ClusterCounter && )

delete

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

delete

ClusterCounter& tutorial::ClusterCounter::operator= ( ClusterCounter && )

delete

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

private

Definition at line 199 of file ClusterCounter5_module.cc.

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

Member Data Documentation

float tutorial::ClusterCounter::fAdcSum

private

Definition at line 83 of file ClusterCounter5_module.cc.

float tutorial::ClusterCounter::fClean

private

Definition at line 84 of file ClusterCounter5_module.cc.

art::InputTag tutorial::ClusterCounter::fClusterModuleLabel

private

Definition at line 87 of file ClusterCounter5_module.cc.

TTree* tutorial::ClusterCounter::fClusterTree

private

Definition at line 81 of file ClusterCounter5_module.cc.

TTree* tutorial::ClusterCounter::fEventTree

private

Definition at line 78 of file ClusterCounter5_module.cc.

size_t tutorial::ClusterCounter::fEvNumber

private

Definition at line 76 of file ClusterCounter5_module.cc.

size_t tutorial::ClusterCounter::fMinSize

private

Definition at line 88 of file ClusterCounter5_module.cc.

size_t tutorial::ClusterCounter::fNClusters

private

Definition at line 79 of file ClusterCounter5_module.cc.

size_t tutorial::ClusterCounter::fNHits

private

Definition at line 82 of file ClusterCounter5_module.cc.

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

protoduneana/protoduneana/TutorialExamples/ClusterCounter5_module.cc
protoduneana/protoduneana/TutorialExamples/ClusterCounter_module.cc

Classes

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation