Inheritance diagram for opdet::OpFlashMCTruthAna:

Public Types
using	ModuleType = EDAnalyzer

Public Member Functions
	OpFlashMCTruthAna (const fhicl::ParameterSet &)

void	analyze (const art::Event &)

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 Attributes
std::string	fFlashInputModule

std::string	fTruthInputModule

TTree *	fAnalysisTree

TTree *	fPerEventTree

Int_t	fEventID

Int_t	fNFlashes

Int_t	fNTruths

Float_t	fFlashY

Float_t	fFlashZ

Float_t	fFlashU

Float_t	fFlashV

Float_t	fFlashT

Float_t	fFlashPE

Float_t	fFlashFastToTotal

Float_t	fFlashWidthY

Float_t	fFlashWidthZ

Float_t	fFlashWidthU

Float_t	fFlashWidthV

Float_t	fVertexX

Float_t	fVertexY

Float_t	fVertexZ

Float_t	fVertexT

Float_t	fTrueE

Int_t	fTruePDG

Float_t	fCenterX

Float_t	fCenterY

Float_t	fCenterZ

Float_t	fDistFlashCenter

Float_t	fDistFlashVertex

Float_t	fDistFlashCenterNorm

Float_t	fDistFlashVertexNorm

Detailed Description

Definition at line 31 of file OpFlashMCTruthAna_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 22 of file EDAnalyzer.h.

Constructor & Destructor Documentation

opdet::OpFlashMCTruthAna::OpFlashMCTruthAna ( const fhicl::ParameterSet & )

Definition at line 81 of file OpFlashMCTruthAna_module.cc.

References fAnalysisTree, fCenterX, fCenterY, fCenterZ, fDistFlashCenter, fDistFlashVertex, fEventID, fFlashFastToTotal, fFlashInputModule, fFlashPE, fFlashT, fFlashU, fFlashV, fFlashWidthU, fFlashWidthV, fFlashWidthY, fFlashWidthZ, fFlashY, fFlashZ, fNFlashes, fPerEventTree, fTrueE, fTruePDG, fTruthInputModule, fVertexX, fVertexY, fVertexZ, and fhicl::ParameterSet::get().

                                                                     : EDAnalyzer(pset)
   {
 
     // Indicate that the Input Module comes from .fcl
     fFlashInputModule = pset.get<std::string>("FlashInputModule");
     fTruthInputModule = pset.get<std::string>("TruthInputModule");
 
     art::ServiceHandle<art::TFileService const> tfs;
 
     fPerEventTree = tfs->make<TTree>("PerEventTree", "PerEventTree");
 
     fPerEventTree->Branch("EventID", &fEventID, "EventID/I");
     fPerEventTree->Branch("NFlashes", &fNFlashes, "NFlashes/I");
 
     fPerEventTree->Branch("VertexX", &fVertexX, "VertexX/F");
     fPerEventTree->Branch("VertexY", &fVertexY, "VertexY/F");
     fPerEventTree->Branch("VertexZ", &fVertexZ, "VertexZ/F");
 
     fPerEventTree->Branch("TrueE", &fTrueE, "TrueE/F");
     fPerEventTree->Branch("TruePDG", &fTruePDG, "TruePDG/I");
 
     fPerEventTree->Branch("CenterX", &fCenterX, "CenterX/F");
     fPerEventTree->Branch("CenterY", &fCenterY, "CenterY/F");
     fPerEventTree->Branch("CenterZ", &fCenterZ, "CenterZ/F");
 
     fAnalysisTree = tfs->make<TTree>("AnalysisTree", "AnalysisTree");
 
     fAnalysisTree->Branch("EventID", &fEventID, "EventID/I");
     fAnalysisTree->Branch("NFlashes", &fNFlashes, "NFlashes/I");
 
     fAnalysisTree->Branch("FlashY", &fFlashY, "FlashY/F");
     fAnalysisTree->Branch("FlashZ", &fFlashZ, "FlashZ/F");
     fAnalysisTree->Branch("FlashU", &fFlashU, "FlashU/F");
     fAnalysisTree->Branch("FlashV", &fFlashV, "FlashV/F");
 
     fAnalysisTree->Branch("FlashWidthY", &fFlashWidthY, "FlashWidthY/F");
     fAnalysisTree->Branch("FlashWidthZ", &fFlashWidthZ, "FlashWidthZ/F");
     fAnalysisTree->Branch("FlashWidthU", &fFlashWidthU, "FlashWidthU/F");
     fAnalysisTree->Branch("FlashWidthV", &fFlashWidthV, "FlashWidthV/F");
 
     fAnalysisTree->Branch("FlashT", &fFlashT, "FlashT/F");
     fAnalysisTree->Branch("FlashPE", &fFlashPE, "FlashPE/F");
     fAnalysisTree->Branch("FlashFastToTotal", &fFlashFastToTotal, "FlashFastToTotal/F");
 
     fAnalysisTree->Branch("VertexX", &fVertexX, "VertexX/F");
     fAnalysisTree->Branch("VertexY", &fVertexY, "VertexY/F");
     fAnalysisTree->Branch("VertexZ", &fVertexZ, "VertexZ/F");
 
     fAnalysisTree->Branch("TrueE", &fTrueE, "TrueE/F");
     fAnalysisTree->Branch("TruePDG", &fTruePDG, "TruePDG/I");
 
     fAnalysisTree->Branch("CenterX", &fCenterX, "CenterX/F");
     fAnalysisTree->Branch("CenterY", &fCenterY, "CenterY/F");
     fAnalysisTree->Branch("CenterZ", &fCenterZ, "CenterZ/F");
 
     fAnalysisTree->Branch("DistFlashCenter", &fDistFlashCenter, "DistFlashCenter/F");
     fAnalysisTree->Branch("DistFlashVertex", &fDistFlashVertex, "DistFlashVertex/F");
   }

Member Function Documentation

void opdet::OpFlashMCTruthAna::analyze ( const art::Event & evt )

Definition at line 141 of file OpFlashMCTruthAna_module.cc.

References simb::MCParticle::E(), simb::MCParticle::EndX(), simb::MCParticle::EndY(), simb::MCParticle::EndZ(), art::EventID::event(), fAnalysisTree, recob::OpFlash::FastToTotal(), fCenterX, fCenterY, fCenterZ, fDistFlashCenter, fDistFlashCenterNorm, fDistFlashVertex, fDistFlashVertexNorm, fEventID, fFlashFastToTotal, fFlashInputModule, fFlashPE, fFlashT, fFlashU, fFlashV, fFlashWidthU, fFlashWidthV, fFlashWidthY, fFlashWidthZ, fFlashY, fFlashZ, fNFlashes, fNTruths, fPerEventTree, fTrueE, fTruePDG, fTruthInputModule, fVertexT, fVertexX, fVertexY, fVertexZ, art::ProductRetriever::getByLabel(), art::Event::id(), simb::MCParticle::PdgCode(), simb::MCParticle::T(), recob::OpFlash::Time(), recob::OpFlash::TotalPE(), simb::MCParticle::Vx(), simb::MCParticle::Vy(), simb::MCParticle::Vz(), recob::OpFlash::WireCenters(), recob::OpFlash::WireWidths(), recob::OpFlash::YCenter(), recob::OpFlash::YWidth(), recob::OpFlash::ZCenter(), and recob::OpFlash::ZWidth().

   {
 
     // Create a handles
     art::Handle<std::vector<recob::OpFlash>> FlashHandle;
     art::Handle<std::vector<simb::MCTruth>> TruthHandle;
 
     // Read in data
     evt.getByLabel(fFlashInputModule, FlashHandle);
     evt.getByLabel(fTruthInputModule, TruthHandle);
 
     fEventID = evt.id().event();
 
     fNTruths = TruthHandle->at(0).NParticles();
     fNFlashes = FlashHandle->size();
 
     std::cout << "Size of truth collection : " << TruthHandle->size() << std::endl;
     std::cout << "We found " << fNTruths << " truth particles and " << fNFlashes << " flashes"
               << std::endl;
 
     for (int iPart = 0; iPart != fNTruths; ++iPart) {
       const simb::MCParticle ThisPart = TruthHandle->at(0).GetParticle(iPart);
       fTruePDG = ThisPart.PdgCode();
       fTrueE = ThisPart.E(0);
 
       fVertexX = ThisPart.Vx(0);
       fVertexY = ThisPart.Vy(0);
       fVertexZ = ThisPart.Vz(0);
       fVertexT = ThisPart.T(0);
 
       fCenterX = (ThisPart.Vx(0) + ThisPart.EndX()) / 2.;
       fCenterY = (ThisPart.Vy(0) + ThisPart.EndY()) / 2.;
       fCenterZ = (ThisPart.Vz(0) + ThisPart.EndZ()) / 2.;
 
       for (unsigned int i = 0; i < FlashHandle->size(); ++i) {
 
         // Get OpFlash
         art::Ptr<recob::OpFlash> TheFlashPtr(FlashHandle, i);
 
         fFlashT = TheFlashPtr->Time();
         fFlashPE = TheFlashPtr->TotalPE();
         fFlashFastToTotal = TheFlashPtr->FastToTotal();
 
         fFlashY = TheFlashPtr->YCenter();
         fFlashZ = TheFlashPtr->ZCenter();
         fFlashU = TheFlashPtr->WireCenters().at(0);
         fFlashV = TheFlashPtr->WireCenters().at(1);
 
         fFlashWidthY = TheFlashPtr->YWidth();
         fFlashWidthZ = TheFlashPtr->ZWidth();
         fFlashWidthU = TheFlashPtr->WireWidths().at(0);
         fFlashWidthV = TheFlashPtr->WireWidths().at(1);
 
         fDistFlashCenter = pow(pow(fCenterY - fFlashY, 2) + pow(fCenterZ - fFlashZ, 2), 0.5);
 
         fDistFlashVertex = pow(pow(fVertexY - fFlashY, 2) + pow(fVertexZ - fFlashZ, 2), 0.5);
 
         fDistFlashCenterNorm = pow(pow((fCenterY - fFlashY) / fFlashWidthY, 2) +
                                      pow((fCenterZ - fFlashZ) / fFlashWidthZ, 2),
                                    0.5);
 
         fDistFlashVertexNorm = pow(pow((fVertexY - fFlashY) / fFlashWidthY, 2) +
                                      pow((fVertexZ - fFlashZ) / fFlashWidthZ, 2),
                                    0.5);
 
         fAnalysisTree->Fill();
       }
     }
 
     fPerEventTree->Fill();
   }

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