Inheritance diagram for lar_pandora::PFParticleTrackAna:

Public Types
using	ModuleType = EDAnalyzer

Public Member Functions
	PFParticleTrackAna (fhicl::ParameterSet const &pset)
	Constructor. More...

virtual	~PFParticleTrackAna ()
	Destructor. More...

void	beginJob ()

void	endJob ()

void	analyze (const art::Event &evt)

void	reconfigure (fhicl::ParameterSet const &pset)

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
TTree *	m_pCaloTree

int	m_run

int	m_event

int	m_index

int	m_ntracks

int	m_trkid

int	m_plane

double	m_length

double	m_dEdx

double	m_dNdx

double	m_dQdx

double	m_residualRange

double	m_x

double	m_y

double	m_z

double	m_px

double	m_py

double	m_pz

bool	m_useModBox

bool	m_isCheated

std::string	m_trackModuleLabel

Detailed Description

PFParticleTrackAna class.

Definition at line 21 of file PFParticleTrackAna_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 22 of file EDAnalyzer.h.

Constructor & Destructor Documentation

lar_pandora::PFParticleTrackAna::PFParticleTrackAna ( fhicl::ParameterSet const & pset )

Constructor.

Parameters

pset

Definition at line 97 of file PFParticleTrackAna_module.cc.

References reconfigure().

                                                                       : art::EDAnalyzer(pset)
   {
     this->reconfigure(pset);
   }

lar_pandora::PFParticleTrackAna::~PFParticleTrackAna ( )

virtual

Destructor.

Definition at line 104 of file PFParticleTrackAna_module.cc.

104 {}

Member Function Documentation

void lar_pandora::PFParticleTrackAna::analyze ( const art::Event & evt )

Definition at line 148 of file PFParticleTrackAna_module.cc.

References lar_pandora::LArPandoraHelper::CollectTracks(), dir, recob::Track::DirectionAtPoint(), art::EventID::event(), art::Event::id(), recob::Track::ID(), recob::Track::Length(), recob::Track::LocationAtPoint(), m_dEdx, m_dNdx, m_dQdx, m_event, m_index, m_length, m_ntracks, m_pCaloTree, m_plane, m_px, m_py, m_pz, m_residualRange, m_run, m_trackModuleLabel, m_trkid, m_x, m_y, m_z, recob::Track::NumberTrajectoryPoints(), art::Event::run(), and track.

   {
     std::cout << " *** PFParticleTrackAna::analyze(...) *** " << std::endl;
 
     m_run = evt.run();
     m_event = evt.id().event();
     m_index = 0;
 
     m_ntracks = 0;
     m_trkid = 0;
     m_plane = 0;
     m_length = 0.0;
     m_dEdx = 0.0;
     m_dNdx = 0.0;
     m_dQdx = 0.0;
     m_residualRange = 0.0;
 
     m_x = 0.0;
     m_y = 0.0;
     m_z = 0.0;
     m_px = 0.0;
     m_py = 0.0;
     m_pz = 0.0;
 
     std::cout << "  Run: " << m_run << std::endl;
     std::cout << "  Event: " << m_event << std::endl;
 
     TrackVector trackVector;
     TracksToHits tracksToHits;
     LArPandoraHelper::CollectTracks(evt, m_trackModuleLabel, trackVector, tracksToHits);
 
     std::cout << "  Tracks: " << trackVector.size() << std::endl;
 
     // art::ServiceHandle<geo::Geometry const> theGeometry;
     // auto const* theDetector = lar::providerFrom<detinfo::DetectorPropertiesService>();
 
 
     // const double adc2eU(5.1e-3);
     // const double adc2eV(5.2e-3);
     // const double adc2eW(5.4e-3);
     // const double adc2eCheat(theDetector->ElectronsToADC());
 
     // const double tau(theDetector->ElectronLifetime());
 
     m_ntracks = trackVector.size();
 
     for (TrackVector::const_iterator iter = trackVector.begin(), iterEnd = trackVector.end();
          iter != iterEnd;
          ++iter) {
       const art::Ptr<recob::Track> track = *iter;
 
       m_trkid = track->ID();
       m_length = track->Length();
 
       m_plane = 0;
       m_dEdx = 0.0;
       m_dNdx = 0.0;
       m_dQdx = 0.0;
       m_residualRange = 0.0;
 
       m_x = 0.0;
       m_y = 0.0;
       m_z = 0.0;
       m_px = 0.0;
       m_py = 0.0;
       m_pz = 0.0;
 
       for (unsigned int p = 0; p < track->NumberTrajectoryPoints(); ++p) {
         auto pos = track->LocationAtPoint(p);
         auto dir = track->DirectionAtPoint(p);
 
         m_residualRange = track->Length(p);
 
         m_x = pos.x();
         m_y = pos.y();
         m_z = pos.z();
         m_px = dir.x();
         m_py = dir.y();
         m_pz = dir.z();
 
         /*************************************************************/
         /*                          WARNING                          */
         /*************************************************************/
         /* The dQdx information in recob::Track has been deprecated  */
         /* since 2016 and in 11/2018 the recob::Track interface was  */
         /* changed and DQdxAtPoint and NumberdQdx were removed.      */
         /* Therefore the code below is now commented out             */
         /* (note that it was most likely not functional anyways).    */
         /* For any issue please contact: larsoft-team@fnal.gov       */
         /*************************************************************/
         /*
             const double dQdxU(track->DQdxAtPoint(p, geo::kU)); // plane 0
             const double dQdxV(track->DQdxAtPoint(p, geo::kV)); // plane 1
             const double dQdxW(track->DQdxAtPoint(p, geo::kW)); // plane 2
 
             m_plane = ((dQdxU > 0.0) ? geo::kU : (dQdxV > 0.0) ? geo::kV : geo::kW);
 
             const double adc2e(m_isCheated ? adc2eCheat : (geo::kU == m_plane) ? adc2eU : (geo::kV == m_plane) ? adc2eV : adc2eW);
 
             m_dQdx = ((geo::kU == m_plane) ? dQdxU : (geo::kV == m_plane) ? dQdxV : dQdxW);
 
             // TODO: Need to include T0 information (currently assume T0 = 0)
 
             m_dNdx = ((m_dQdx / adc2e) * exp((m_x / theDetector->GetXTicksCoefficient()) * theDetector->SamplingRate() * 1.e-3 / tau));
 
             m_dEdx = (m_useModBox ? theDetector->ModBoxCorrection(m_dNdx) : theDetector->BirksCorrection(m_dNdx));
             */
         /*************************************************************/
 
         m_pCaloTree->Fill();
         ++m_index;
       }
     }
   }

void lar_pandora::PFParticleTrackAna::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 117 of file PFParticleTrackAna_module.cc.

References m_dEdx, m_dNdx, m_dQdx, m_event, m_index, m_length, m_ntracks, m_pCaloTree, m_plane, m_px, m_py, m_pz, m_residualRange, m_run, m_trkid, m_x, m_y, and m_z.

   {
     //
     art::ServiceHandle<art::TFileService const> tfs;
 
     m_pCaloTree = tfs->make<TTree>("calorimetry", "LAr Track Calo Tree");
     m_pCaloTree->Branch("run", &m_run, "run/I");
     m_pCaloTree->Branch("event", &m_event, "event/I");
     m_pCaloTree->Branch("index", &m_index, "index/I");
     m_pCaloTree->Branch("ntracks", &m_ntracks, "ntracks/I");
     m_pCaloTree->Branch("trkid", &m_trkid, "trkid/I");
     m_pCaloTree->Branch("plane", &m_plane, "plane/I");
     m_pCaloTree->Branch("length", &m_length, "length/D");
     m_pCaloTree->Branch("dEdx", &m_dEdx, "dEdx/D");
     m_pCaloTree->Branch("dNdx", &m_dNdx, "dNdx/D");
     m_pCaloTree->Branch("dQdx", &m_dQdx, "dQdx/D");
     m_pCaloTree->Branch("residualRange", &m_residualRange, "residualRange/D");
     m_pCaloTree->Branch("x", &m_x, "x/D");
     m_pCaloTree->Branch("y", &m_y, "y/D");
     m_pCaloTree->Branch("z", &m_z, "z/D");
     m_pCaloTree->Branch("px", &m_px, "px/D");
     m_pCaloTree->Branch("py", &m_py, "py/D");
     m_pCaloTree->Branch("pz", &m_pz, "pz/D");
   }

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

void lar_pandora::PFParticleTrackAna::endJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 144 of file PFParticleTrackAna_module.cc.

144 {}

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 lar_pandora::PFParticleTrackAna::reconfigure ( fhicl::ParameterSet const & pset )

Definition at line 108 of file PFParticleTrackAna_module.cc.

References fhicl::ParameterSet::get(), m_isCheated, m_trackModuleLabel, and m_useModBox.

Referenced by PFParticleTrackAna().

   {
     m_useModBox = pset.get<bool>("UeModBox", true);
     m_isCheated = pset.get<bool>("IsCheated", false);
     m_trackModuleLabel = pset.get<std::string>("TrackModule", "pandora");
   }

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