Inheritance diagram for pdsp::MDMAna:

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

	MDMAna (MDMAna const &)=delete

	MDMAna (MDMAna &&)=delete

MDMAna &	operator= (MDMAna const &)=delete

MDMAna &	operator= (MDMAna &&)=delete

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

void	beginJob () 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 Attributes
TTree *	mdmtree

int	run

int	subrun

int	event

std::vector< int >	trackid

std::vector< double >	tracklen

std::vector< int >	trackorg

std::vector< int >	trackpdg

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 MDMAna_module.cc.

Constructor & Destructor Documentation

pdsp::MDMAna::MDMAna ( fhicl::ParameterSet const & p )

explicit

Definition at line 68 of file MDMAna_module.cc.

   : EDAnalyzer{p}  // ,
   // More initializers here.
 {
   // Call appropriate consumes<>() for any products to be retrieved by this module.
 }

pdsp::MDMAna::MDMAna ( MDMAna const & )

delete

pdsp::MDMAna::MDMAna ( MDMAna && )

delete

Member Function Documentation

void pdsp::MDMAna::analyze ( art::Event const & e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 75 of file MDMAna_module.cc.

 {
   run = e.run();
   subrun = e.subRun();
   event = e.id().event();
   trackid.clear();
   tracklen.clear();
   trackorg.clear();
   trackpdg.clear();
 
   //Services
   art::ServiceHandle<cheat::BackTrackerService> bt_serv;
   art::ServiceHandle<cheat::ParticleInventoryService> pi_serv;
 
    // Reconstruciton information
   std::vector < art::Ptr < recob::Track > > trackList;
   auto trackListHandle = e.getHandle < std::vector < recob::Track > >("pandoraTrack");
   if (trackListHandle) {
     art::fill_ptr_vector(trackList, trackListHandle);
   }
   else return;
 
   //Get hits associated with track
   art::FindManyP < recob::Hit > hitsFromTrack(trackListHandle, e, "pandoraTrack");
 
   auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService>()->DataFor(e);
   for (auto const& track : trackList){
     int this_trackid = track.key();
     double this_tracklen = track->Length();
     int this_trackorg = -1;
     int this_trackpdg = 0;
 
     auto const & allHits = hitsFromTrack.at(track.key());
 
     if (!e.isRealData()){
       // Find true particle for reconstructed track
       int TrackID = 0;
       std::map<int,double> trkide;
       for(size_t h = 0; h < allHits.size(); ++h){
         art::Ptr<recob::Hit> hit = allHits[h];
         std::vector<sim::TrackIDE> TrackIDs = bt_serv->HitToTrackIDEs(clockData, hit);
         for(size_t e = 0; e < TrackIDs.size(); ++e){
           trkide[TrackIDs[e].trackID] += TrackIDs[e].energy;
         }           
       }
       // Work out which IDE despoited the most charge in the hit if there was more than one.
       double maxe = -1;
       double tote = 0;
       for (std::map<int,double>::iterator ii = trkide.begin(); ii!=trkide.end(); ++ii){
         tote += ii->second;
         if ((ii->second)>maxe){
           maxe = ii->second;
           TrackID = ii->first;
         }
       }
       // Now have trackID, so get PdG code and T0 etc.
       if (TrackID){
         const simb::MCParticle *particle = pi_serv->TrackIdToParticle_P(TrackID);
         if (particle){
           this_trackpdg = particle->PdgCode();
           auto & truth = pi_serv->ParticleToMCTruth_P(particle);
           this_trackorg = truth->Origin();
         }
       }
     }
     trackid.push_back(this_trackid);
     tracklen.push_back(this_tracklen);
     trackorg.push_back(this_trackorg);
     trackpdg.push_back(this_trackpdg);
   }
   if (!trackid.empty()) mdmtree->Fill();
 }

void pdsp::MDMAna::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 148 of file MDMAna_module.cc.

 {
   art::ServiceHandle<art::TFileService> fileServiceHandle;
   mdmtree = fileServiceHandle->make<TTree>("mdmtree", "MDM info");
   mdmtree->Branch("run", &run, "run/I");
   mdmtree->Branch("subrun", &subrun, "subrun/I");
   mdmtree->Branch("event", &event, "event/I");
   mdmtree->Branch("trackid", &trackid);
   mdmtree->Branch("tracklen", &tracklen);
   mdmtree->Branch("trackorg", &trackorg);
   mdmtree->Branch("trackpdg", &trackpdg);
 
 }

MDMAna& pdsp::MDMAna::operator= ( MDMAna const & )

delete

MDMAna& pdsp::MDMAna::operator= ( MDMAna && )

delete

Member Data Documentation

int pdsp::MDMAna::event

private

Definition at line 59 of file MDMAna_module.cc.

TTree* pdsp::MDMAna::mdmtree

private

Definition at line 56 of file MDMAna_module.cc.

int pdsp::MDMAna::run

private

Definition at line 57 of file MDMAna_module.cc.

int pdsp::MDMAna::subrun

private

Definition at line 58 of file MDMAna_module.cc.

std::vector<int> pdsp::MDMAna::trackid

private

Definition at line 60 of file MDMAna_module.cc.

std::vector<double> pdsp::MDMAna::tracklen

private

Definition at line 61 of file MDMAna_module.cc.

std::vector<int> pdsp::MDMAna::trackorg

private

Definition at line 62 of file MDMAna_module.cc.

std::vector<int> pdsp::MDMAna::trackpdg

private

Definition at line 63 of file MDMAna_module.cc.

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

protoduneana/protoduneana/singlephase/MDMAna/MDMAna_module.cc

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation