lar_pandora::PFParticleMonitoring Class Reference

PFParticleMonitoring class. More...

Inheritance diagram for lar_pandora::PFParticleMonitoring:

Public Types
using	ModuleType = EDAnalyzer

Public Member Functions
	PFParticleMonitoring (fhicl::ParameterSet const &pset)
	Constructor. More...
virtual	~PFParticleMonitoring ()
	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 Types
typedef std::set< art::Ptr< recob::PFParticle > >	PFParticleSet
typedef std::set< art::Ptr< simb::MCParticle > >	MCParticleSet
typedef std::set< art::Ptr< simb::MCTruth > >	MCTruthSet

Private Member Functions
void	BuildTrueNeutrinoHitMaps (const MCTruthToMCParticles &truthToParticles, const MCParticlesToHits &trueParticlesToHits, MCTruthToHits &trueNeutrinosToHits, HitsToMCTruth &trueHitsToNeutrinos) const
	Build mapping from true neutrinos to hits. More...
void	BuildRecoNeutrinoHitMaps (const PFParticleMap &recoParticleMap, const PFParticlesToHits &recoParticlesToHits, PFParticlesToHits &recoNeutrinosToHits, HitsToPFParticles &recoHitsToNeutrinos) const
	Build mapping from reconstructed neutrinos to hits. More...
void	GetRecoToTrueMatches (const PFParticlesToHits &recoNeutrinosToHits, const HitsToMCTruth &trueHitsToNeutrinos, MCTruthToPFParticles &matchedNeutrinos, MCTruthToHits &matchedNeutrinoHits) const
	Perform matching between true and reconstructed neutrino events. More...
void	GetRecoToTrueMatches (const PFParticlesToHits &recoNeutrinosToHits, const HitsToMCTruth &trueHitsToNeutrinos, MCTruthToPFParticles &matchedNeutrinos, MCTruthToHits &matchedNeutrinoHits, PFParticleSet &recoVeto, MCTruthSet &trueVeto) const
	Perform matching between true and reconstructed neutrino events. More...
void	GetRecoToTrueMatches (const PFParticlesToHits &recoParticlesToHits, const HitsToMCParticles &trueHitsToParticles, MCParticlesToPFParticles &matchedParticles, MCParticlesToHits &matchedHits) const
	Perform matching between true and reconstructed particles. More...
void	GetRecoToTrueMatches (const PFParticlesToHits &recoParticlesToHits, const HitsToMCParticles &trueHitsToParticles, MCParticlesToPFParticles &matchedParticles, MCParticlesToHits &matchedHits, PFParticleSet &recoVeto, MCParticleSet &trueVeto) const
	Perform matching between true and reconstructed particles. More...
int	CountHitsByType (const int view, const HitVector &hitVector) const
	Count the number of reconstructed hits in a given wire plane. More...
void	GetStartAndEndPoints (const art::Ptr< simb::MCParticle > trueParticle, int &startT, int &endT) const
	Find the start and end points of the true particle in the active region of detector. More...
double	GetLength (const art::Ptr< simb::MCParticle > trueParticle, const int startT, const int endT) const
	Find the length of the true particle trajectory through the active region of the detector. More...

Private Attributes
TTree *	m_pRecoTree
int	m_run
int	m_event
int	m_index
int	m_nMCParticles
int	m_nNeutrinoPfos
int	m_nPrimaryPfos
int	m_nDaughterPfos
int	m_mcPdg
int	m_mcNuPdg
int	m_mcParentPdg
int	m_mcPrimaryPdg
int	m_mcIsNeutrino
int	m_mcIsPrimary
int	m_mcIsDecay
int	m_mcIsCC
int	m_pfoPdg
int	m_pfoNuPdg
int	m_pfoParentPdg
int	m_pfoPrimaryPdg
int	m_pfoIsNeutrino
int	m_pfoIsPrimary
int	m_pfoIsStitched
int	m_pfoTrack
int	m_pfoVertex
double	m_pfoVtxX
double	m_pfoVtxY
double	m_pfoVtxZ
double	m_pfoEndX
double	m_pfoEndY
double	m_pfoEndZ
double	m_pfoDirX
double	m_pfoDirY
double	m_pfoDirZ
double	m_pfoLength
double	m_pfoStraightLength
int	m_mcVertex
double	m_mcVtxX
double	m_mcVtxY
double	m_mcVtxZ
double	m_mcEndX
double	m_mcEndY
double	m_mcEndZ
double	m_mcDirX
double	m_mcDirY
double	m_mcDirZ
double	m_mcEnergy
double	m_mcLength
double	m_mcStraightLength
double	m_completeness
double	m_purity
int	m_nMCHits
int	m_nPfoHits
int	m_nMatchedHits
int	m_nMCHitsU
int	m_nMCHitsV
int	m_nMCHitsW
int	m_nPfoHitsU
int	m_nPfoHitsV
int	m_nPfoHitsW
int	m_nMatchedHitsU
int	m_nMatchedHitsV
int	m_nMatchedHitsW
int	m_nTrueWithoutRecoHits
	True hits which don't belong to any reconstructed particle - "available". More...
int	m_nRecoWithoutTrueHits
	Reconstructed hits which don't belong to any true particle - "missing". More...
double	m_spacepointsMinX
double	m_spacepointsMaxX
std::string	m_hitfinderLabel
std::string	m_trackLabel
std::string	m_particleLabel
std::string	m_backtrackerLabel
std::string	m_geantModuleLabel
bool	m_useDaughterPFParticles
bool	m_useDaughterMCParticles
bool	m_addDaughterPFParticles
bool	m_addDaughterMCParticles
bool	m_recursiveMatching
bool	m_printDebug
	switch for print statements (TODO: use message service!) More...
bool	m_disableRealDataCheck
	Whether to check if the input file contains real data before accessing MC information. More...

Detailed Description

PFParticleMonitoring class.

Definition at line 24 of file PFParticleMonitoring_module.cc.

Member Typedef Documentation

typedef std::set<art::Ptr<simb::MCParticle> > lar_pandora::PFParticleMonitoring::MCParticleSet

private

Definition at line 45 of file PFParticleMonitoring_module.cc.

typedef std::set<art::Ptr<simb::MCTruth> > lar_pandora::PFParticleMonitoring::MCTruthSet

private

Definition at line 46 of file PFParticleMonitoring_module.cc.

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 22 of file EDAnalyzer.h.

typedef std::set<art::Ptr<recob::PFParticle> > lar_pandora::PFParticleMonitoring::PFParticleSet

private

Definition at line 44 of file PFParticleMonitoring_module.cc.

Constructor & Destructor Documentation

lar_pandora::PFParticleMonitoring::PFParticleMonitoring

(

fhicl::ParameterSet const &

pset

)

Constructor.

Parameters

pset

Definition at line 300 of file PFParticleMonitoring_module.cc.

References reconfigure().

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

lar_pandora::PFParticleMonitoring::~PFParticleMonitoring ( )

virtual

Destructor.

Definition at line 308 of file PFParticleMonitoring_module.cc.

{}

Member Function Documentation

void lar_pandora::PFParticleMonitoring::analyze

(

const art::Event &

evt

)

Definition at line 414 of file PFParticleMonitoring_module.cc.

References lar_pandora::LArPandoraHelper::BuildMCParticleHitMaps(), lar_pandora::LArPandoraHelper::BuildMCParticleMap(), lar_pandora::LArPandoraHelper::BuildPFParticleHitMaps(), lar_pandora::LArPandoraHelper::BuildPFParticleMap(), BuildRecoNeutrinoHitMaps(), BuildTrueNeutrinoHitMaps(), lar_pandora::LArPandoraHelper::CollectHits(), lar_pandora::LArPandoraHelper::CollectMCParticles(), lar_pandora::LArPandoraHelper::CollectPFParticles(), lar_pandora::LArPandoraHelper::CollectSpacePoints(), lar_pandora::LArPandoraHelper::CollectT0s(), lar_pandora::LArPandoraHelper::CollectTracks(), lar_pandora::LArPandoraHelper::CollectVertices(), CountHitsByType(), e, simb::MCParticle::E(), recob::Track::End(), art::EventID::event(), lar_pandora::LArPandoraHelper::GetFinalStateMCParticle(), lar_pandora::LArPandoraHelper::GetFinalStatePFParticle(), GetLength(), simb::MCTruth::GetNeutrino(), lar_pandora::LArPandoraHelper::GetParentMCParticle(), lar_pandora::LArPandoraHelper::GetParentNeutrino(), lar_pandora::LArPandoraHelper::GetParentPFParticle(), GetRecoToTrueMatches(), GetStartAndEndPoints(), art::Event::id(), lar_pandora::LArPandoraHelper::IsFinalState(), lar_pandora::LArPandoraHelper::IsNeutrino(), art::Event::isRealData(), lar_pandora::LArPandoraHelper::kAddDaughters, simb::kCC, lar_pandora::LArPandoraHelper::kIgnoreDaughters, geo::kU, lar_pandora::LArPandoraHelper::kUseDaughters, geo::kV, geo::kW, recob::Track::Length(), m_addDaughterMCParticles, m_addDaughterPFParticles, m_backtrackerLabel, m_completeness, m_disableRealDataCheck, m_event, m_geantModuleLabel, m_hitfinderLabel, m_index, m_mcDirX, m_mcDirY, m_mcDirZ, m_mcEndX, m_mcEndY, m_mcEndZ, m_mcEnergy, m_mcIsCC, m_mcIsDecay, m_mcIsNeutrino, m_mcIsPrimary, m_mcLength, m_mcNuPdg, m_mcParentPdg, m_mcPdg, m_mcPrimaryPdg, m_mcStraightLength, m_mcVertex, m_mcVtxX, m_mcVtxY, m_mcVtxZ, m_nDaughterPfos, m_nMatchedHits, m_nMatchedHitsU, m_nMatchedHitsV, m_nMatchedHitsW, m_nMCHits, m_nMCHitsU, m_nMCHitsV, m_nMCHitsW, m_nMCParticles, m_nNeutrinoPfos, m_nPfoHits, m_nPfoHitsU, m_nPfoHitsV, m_nPfoHitsW, m_nPrimaryPfos, m_nRecoWithoutTrueHits, m_nTrueWithoutRecoHits, m_particleLabel, m_pfoDirX, m_pfoDirY, m_pfoDirZ, m_pfoEndX, m_pfoEndY, m_pfoEndZ, m_pfoIsNeutrino, m_pfoIsPrimary, m_pfoIsStitched, m_pfoLength, m_pfoNuPdg, m_pfoParentPdg, m_pfoPdg, m_pfoPrimaryPdg, m_pfoStraightLength, m_pfoTrack, m_pfoVertex, m_pfoVtxX, m_pfoVtxY, m_pfoVtxZ, m_pRecoTree, m_printDebug, m_purity, m_run, m_spacepointsMaxX, m_spacepointsMinX, m_trackLabel, m_useDaughterMCParticles, m_useDaughterPFParticles, simb::MCTruth::NeutrinoSet(), simb::MCParticle::P(), recob::PFParticle::PdgCode(), simb::MCParticle::PdgCode(), simb::MCParticle::Process(), simb::MCParticle::Px(), simb::MCParticle::Py(), simb::MCParticle::Pz(), art::Event::run(), lar_pandora::LArPandoraHelper::SelectNeutrinoPFParticles(), recob::Track::Vertex(), recob::Track::VertexDirection(), simb::MCParticle::Vx(), simb::MCParticle::Vy(), simb::MCParticle::Vz(), X, recob::SpacePoint::XYZ(), and recob::Vertex::XYZ().

  {
    if (m_printDebug) std::cout << " *** PFParticleMonitoring::analyze(...) *** " << std::endl;

    m_run = evt.run();
    m_event = evt.id().event();
    m_index = 0;

    m_nMCParticles = 0;
    m_nNeutrinoPfos = 0;
    m_nPrimaryPfos = 0;
    m_nDaughterPfos = 0;

    m_mcPdg = 0;
    m_mcNuPdg = 0;
    m_mcParentPdg = 0;
    m_mcPrimaryPdg = 0;
    m_mcIsNeutrino = 0;
    m_mcIsPrimary = 0;
    m_mcIsDecay = 0;
    m_mcIsCC = 0;

    m_pfoPdg = 0;
    m_pfoNuPdg = 0;
    m_pfoParentPdg = 0;
    m_pfoPrimaryPdg = 0;
    m_pfoIsNeutrino = 0;
    m_pfoIsPrimary = 0;
    m_pfoIsStitched = 0;
    m_pfoTrack = 0;
    m_pfoVertex = 0;
    m_pfoVtxX = 0.0;
    m_pfoVtxY = 0.0;
    m_pfoVtxZ = 0.0;
    m_pfoEndX = 0.0;
    m_pfoEndY = 0.0;
    m_pfoEndZ = 0.0;
    m_pfoDirX = 0.0;
    m_pfoDirY = 0.0;
    m_pfoDirZ = 0.0;
    m_pfoLength = 0.0;
    m_pfoStraightLength = 0.0;

    m_mcVertex = 0;
    m_mcVtxX = 0.0;
    m_mcVtxY = 0.0;
    m_mcVtxZ = 0.0;
    m_mcEndX = 0.0;
    m_mcEndY = 0.0;
    m_mcEndZ = 0.0;
    m_mcDirX = 0.0;
    m_mcDirY = 0.0;
    m_mcDirZ = 0.0;
    m_mcEnergy = 0.0;
    m_mcLength = 0.0;
    m_mcStraightLength = 0.0;

    m_completeness = 0.0;
    m_purity = 0.0;

    m_nMCHits = 0;
    m_nPfoHits = 0;
    m_nMatchedHits = 0;
    m_nMCHitsU = 0;
    m_nMCHitsV = 0;
    m_nMCHitsW = 0;
    m_nPfoHitsU = 0;
    m_nPfoHitsV = 0;
    m_nPfoHitsW = 0;
    m_nMatchedHitsU = 0;
    m_nMatchedHitsV = 0;
    m_nMatchedHitsW = 0;

    m_nTrueWithoutRecoHits = 0;
    m_nRecoWithoutTrueHits = 0;

    m_spacepointsMinX = 0.0;
    m_spacepointsMaxX = 0.0;

    if (m_printDebug) {
      std::cout << "  Run: " << m_run << std::endl;
      std::cout << "  Event: " << m_event << std::endl;
    }

    // Collect Hits
    // ============
    HitVector hitVector;
    LArPandoraHelper::CollectHits(evt, m_hitfinderLabel, hitVector);

    if (m_printDebug) std::cout << "  Hits: " << hitVector.size() << std::endl;

    // Collect SpacePoints and SpacePoint <-> Hit Associations
    // =======================================================
    SpacePointVector spacePointVector;
    SpacePointsToHits spacePointsToHits;
    HitsToSpacePoints hitsToSpacePoints;
    LArPandoraHelper::CollectSpacePoints(
      evt, m_particleLabel, spacePointVector, spacePointsToHits, hitsToSpacePoints);

    if (m_printDebug) std::cout << "  SpacePoints: " << spacePointVector.size() << std::endl;

    // Collect Tracks and PFParticle <-> Track Associations
    // ====================================================
    TrackVector recoTrackVector;
    PFParticlesToTracks recoParticlesToTracks;
    LArPandoraHelper::CollectTracks(evt, m_trackLabel, recoTrackVector, recoParticlesToTracks);

    if (m_printDebug) std::cout << "  Tracks: " << recoTrackVector.size() << std::endl;

    // Collect TOs and PFParticle <-> T0 Associations
    // ==============================================
    T0Vector t0Vector;
    PFParticlesToT0s particlesToT0s;
    LArPandoraHelper::CollectT0s(evt, m_particleLabel, t0Vector, particlesToT0s);

    // Collect Vertices and PFParticle <-> Vertex Associations
    // =======================================================
    VertexVector recoVertexVector;
    PFParticlesToVertices recoParticlesToVertices;
    LArPandoraHelper::CollectVertices(
      evt, m_particleLabel, recoVertexVector, recoParticlesToVertices);

    if (m_printDebug) std::cout << "  Vertices: " << recoVertexVector.size() << std::endl;

    // Collect PFParticles and match Reco Particles to Hits
    // ====================================================
    PFParticleVector recoParticleVector;
    PFParticleVector recoNeutrinoVector;
    PFParticlesToHits recoParticlesToHits;
    HitsToPFParticles recoHitsToParticles;

    LArPandoraHelper::CollectPFParticles(evt, m_particleLabel, recoParticleVector);
    LArPandoraHelper::SelectNeutrinoPFParticles(recoParticleVector, recoNeutrinoVector);
    LArPandoraHelper::BuildPFParticleHitMaps(
      evt,
      m_particleLabel,
      recoParticlesToHits,
      recoHitsToParticles,
      (m_useDaughterPFParticles ? (m_addDaughterPFParticles ? LArPandoraHelper::kAddDaughters :
                                                              LArPandoraHelper::kUseDaughters) :
                                  LArPandoraHelper::kIgnoreDaughters));

    if (m_printDebug) {
      std::cout << "  RecoNeutrinos: " << recoNeutrinoVector.size() << std::endl;
      std::cout << "  RecoParticles: " << recoParticleVector.size() << std::endl;
    }

    // Collect MCParticles and match True Particles to Hits
    // ====================================================
    MCParticleVector trueParticleVector;
    MCTruthToMCParticles truthToParticles;
    MCParticlesToMCTruth particlesToTruth;
    MCParticlesToHits trueParticlesToHits;
    HitsToMCParticles trueHitsToParticles;

    if (m_disableRealDataCheck || !evt.isRealData()) {
      LArPandoraHelper::CollectMCParticles(evt, m_geantModuleLabel, trueParticleVector);
      LArPandoraHelper::CollectMCParticles(
        evt, m_geantModuleLabel, truthToParticles, particlesToTruth);

      LArPandoraHelper::BuildMCParticleHitMaps(
        evt,
        m_geantModuleLabel,
        hitVector,
        trueParticlesToHits,
        trueHitsToParticles,
        (m_useDaughterMCParticles ? (m_addDaughterMCParticles ? LArPandoraHelper::kAddDaughters :
                                                                LArPandoraHelper::kUseDaughters) :
                                    LArPandoraHelper::kIgnoreDaughters));

      if (trueHitsToParticles.empty()) {
        if (m_backtrackerLabel.empty())
          throw cet::exception("LArPandora") << " PFParticleMonitoring::analyze - no sim channels "
                                                "found, backtracker module must be set in FHiCL "
                                             << std::endl;

        LArPandoraHelper::BuildMCParticleHitMaps(
          evt,
          m_geantModuleLabel,
          m_hitfinderLabel,
          m_backtrackerLabel,
          trueParticlesToHits,
          trueHitsToParticles,
          (m_useDaughterMCParticles ? (m_addDaughterMCParticles ? LArPandoraHelper::kAddDaughters :
                                                                  LArPandoraHelper::kUseDaughters) :
                                      LArPandoraHelper::kIgnoreDaughters));
      }
    }

    if (m_printDebug) {
      std::cout << "  TrueParticles: " << particlesToTruth.size() << std::endl;
      std::cout << "  TrueEvents: " << truthToParticles.size() << std::endl;
      std::cout << "  MatchedParticles: " << trueParticlesToHits.size() << std::endl;
    }

    if (trueParticlesToHits.empty()) {
      m_pRecoTree->Fill();
      return;
    }

    // Build Reco and True Particle Maps (for Parent/Daughter Navigation)
    // =================================================================
    MCParticleMap trueParticleMap;
    PFParticleMap recoParticleMap;

    LArPandoraHelper::BuildMCParticleMap(trueParticleVector, trueParticleMap);
    LArPandoraHelper::BuildPFParticleMap(recoParticleVector, recoParticleMap);

    m_nMCParticles = trueParticlesToHits.size();
    m_nNeutrinoPfos = 0;
    m_nPrimaryPfos = 0;
    m_nDaughterPfos = 0;

    // Count reconstructed particles
    for (PFParticleVector::const_iterator iter = recoParticleVector.begin(),
                                          iterEnd = recoParticleVector.end();
         iter != iterEnd;
         ++iter) {
      const art::Ptr<recob::PFParticle> recoParticle = *iter;

      if (LArPandoraHelper::IsNeutrino(recoParticle)) { m_nNeutrinoPfos++; }
      else if (LArPandoraHelper::IsFinalState(recoParticleMap, recoParticle)) {
        m_nPrimaryPfos++;
      }
      else {
        m_nDaughterPfos++;
      }
    }

    // Match Reco Neutrinos to True Neutrinos
    // ======================================
    PFParticlesToHits recoNeutrinosToHits;
    HitsToPFParticles recoHitsToNeutrinos;
    HitsToMCTruth trueHitsToNeutrinos;
    MCTruthToHits trueNeutrinosToHits;
    this->BuildRecoNeutrinoHitMaps(
      recoParticleMap, recoParticlesToHits, recoNeutrinosToHits, recoHitsToNeutrinos);
    this->BuildTrueNeutrinoHitMaps(
      truthToParticles, trueParticlesToHits, trueNeutrinosToHits, trueHitsToNeutrinos);

    MCTruthToPFParticles matchedNeutrinos;
    MCTruthToHits matchedNeutrinoHits;
    this->GetRecoToTrueMatches(
      recoNeutrinosToHits, trueHitsToNeutrinos, matchedNeutrinos, matchedNeutrinoHits);

    for (MCTruthToHits::const_iterator iter = trueNeutrinosToHits.begin(),
                                       iterEnd = trueNeutrinosToHits.end();
         iter != iterEnd;
         ++iter) {
      const art::Ptr<simb::MCTruth> trueEvent = iter->first;
      const HitVector& trueHitVector = iter->second;

      if (trueHitVector.empty()) continue;

      if (!trueEvent->NeutrinoSet()) continue;

      const simb::MCNeutrino trueNeutrino(trueEvent->GetNeutrino());
      const simb::MCParticle trueParticle(trueNeutrino.Nu());

      m_mcIsCC = ((simb::kCC == trueNeutrino.CCNC()) ? 1 : 0);
      m_mcPdg = trueParticle.PdgCode();
      m_mcNuPdg = m_mcPdg;
      m_mcParentPdg = 0;
      m_mcPrimaryPdg = 0;
      m_mcIsNeutrino = 1;
      m_mcIsPrimary = 0;
      m_mcIsDecay = 0;

      m_mcVertex = 1;
      m_mcVtxX = trueParticle.Vx();
      m_mcVtxY = trueParticle.Vy();
      m_mcVtxZ = trueParticle.Vz();
      m_mcEndX = m_mcVtxX;
      m_mcEndY = m_mcVtxY;
      m_mcEndZ = m_mcVtxZ;
      m_mcDirX = trueParticle.Px() / trueParticle.P();
      m_mcDirY = trueParticle.Py() / trueParticle.P();
      m_mcDirZ = trueParticle.Pz() / trueParticle.P();
      m_mcEnergy = trueParticle.E();
      m_mcLength = 0.0;
      m_mcStraightLength = 0.0;

      m_nMCHits = trueHitVector.size();
      m_nMCHitsU = this->CountHitsByType(geo::kU, trueHitVector);
      m_nMCHitsV = this->CountHitsByType(geo::kV, trueHitVector);
      m_nMCHitsW = this->CountHitsByType(geo::kW, trueHitVector);

      m_pfoPdg = 0;
      m_pfoNuPdg = 0;
      m_pfoParentPdg = 0;
      m_pfoPrimaryPdg = 0;
      m_pfoIsNeutrino = 0;
      m_pfoIsPrimary = 0;
      m_pfoIsStitched = 0;
      m_pfoTrack = 0;
      m_pfoVertex = 0;
      m_pfoVtxX = 0.0;
      m_pfoVtxY = 0.0;
      m_pfoVtxZ = 0.0;
      m_pfoEndX = 0.0;
      m_pfoEndY = 0.0;
      m_pfoEndZ = 0.0;
      m_pfoDirX = 0.0;
      m_pfoDirY = 0.0;
      m_pfoDirZ = 0.0;
      m_pfoLength = 0.0;
      m_pfoStraightLength = 0.0;

      m_nPfoHits = 0;
      m_nPfoHitsU = 0;
      m_nPfoHitsV = 0;
      m_nPfoHitsW = 0;

      m_nMatchedHits = 0;
      m_nMatchedHitsU = 0;
      m_nMatchedHitsV = 0;
      m_nMatchedHitsW = 0;

      m_nTrueWithoutRecoHits = 0;
      m_nRecoWithoutTrueHits = 0;

      m_spacepointsMinX = 0.0;
      m_spacepointsMaxX = 0.0;

      m_completeness = 0.0;
      m_purity = 0.0;

      for (HitVector::const_iterator hIter1 = trueHitVector.begin(),
                                     hIterEnd1 = trueHitVector.end();
           hIter1 != hIterEnd1;
           ++hIter1) {
        if (recoHitsToNeutrinos.find(*hIter1) == recoHitsToNeutrinos.end())
          ++m_nTrueWithoutRecoHits;
      }

      MCTruthToPFParticles::const_iterator pIter1 = matchedNeutrinos.find(trueEvent);
      if (matchedNeutrinos.end() != pIter1) {
        const art::Ptr<recob::PFParticle> recoParticle = pIter1->second;
        m_pfoPdg = recoParticle->PdgCode();
        m_pfoNuPdg = m_pfoPdg;
        m_pfoParentPdg = m_pfoPdg;
        m_pfoPrimaryPdg = 0;
        m_pfoIsNeutrino = 1;
        m_pfoIsPrimary = 0;

        if (!LArPandoraHelper::IsNeutrino(recoParticle))
          std::cout << " Warning: Found neutrino with an invalid PDG code " << std::endl;

        PFParticlesToHits::const_iterator pIter2 = recoNeutrinosToHits.find(recoParticle);
        if (recoParticlesToHits.end() != pIter2) {
          const HitVector& recoHitVector = pIter2->second;

          for (HitVector::const_iterator hIter2 = recoHitVector.begin(),
                                         hIterEnd2 = recoHitVector.end();
               hIter2 != hIterEnd2;
               ++hIter2) {
            if (trueHitsToNeutrinos.find(*hIter2) == trueHitsToNeutrinos.end())
              ++m_nRecoWithoutTrueHits;
          }

          MCTruthToHits::const_iterator pIter3 = matchedNeutrinoHits.find(trueEvent);
          if (matchedNeutrinoHits.end() != pIter3) {
            const HitVector& matchedHitVector = pIter3->second;

            m_nPfoHits = recoHitVector.size();
            m_nPfoHitsU = this->CountHitsByType(geo::kU, recoHitVector);
            m_nPfoHitsV = this->CountHitsByType(geo::kV, recoHitVector);
            m_nPfoHitsW = this->CountHitsByType(geo::kW, recoHitVector);

            m_nMatchedHits = matchedHitVector.size();
            m_nMatchedHitsU = this->CountHitsByType(geo::kU, matchedHitVector);
            m_nMatchedHitsV = this->CountHitsByType(geo::kV, matchedHitVector);
            m_nMatchedHitsW = this->CountHitsByType(geo::kW, matchedHitVector);

            PFParticlesToVertices::const_iterator pIter4 =
              recoParticlesToVertices.find(recoParticle);
            if (recoParticlesToVertices.end() != pIter4) {
              const VertexVector& vertexVector = pIter4->second;
              if (!vertexVector.empty()) {
                if (vertexVector.size() != 1 && m_printDebug)
                  std::cout << " Warning: Found particle with more than one associated vertex "
                            << std::endl;

                const art::Ptr<recob::Vertex> recoVertex = *(vertexVector.begin());
                double xyz[3] = {0.0, 0.0, 0.0};
                recoVertex->XYZ(xyz);

                m_pfoVertex = 1;
                m_pfoVtxX = xyz[0];
                m_pfoVtxY = xyz[1];
                m_pfoVtxZ = xyz[2];
              }
            }
          }
        }
      }

      m_purity =
        ((m_nPfoHits == 0) ? 0.0 :
                             static_cast<double>(m_nMatchedHits) / static_cast<double>(m_nPfoHits));
      m_completeness =
        ((m_nPfoHits == 0) ? 0.0 :
                             static_cast<double>(m_nMatchedHits) / static_cast<double>(m_nMCHits));

      if (m_printDebug)
        std::cout << "    MCNeutrino [" << m_index << "]"
                  << "  trueNu=" << m_mcNuPdg << ", truePdg=" << m_mcPdg
                  << ", recoNu=" << m_pfoNuPdg << ", recoPdg=" << m_pfoPdg
                  << ", mcHits=" << m_nMCHits << ", pfoHits=" << m_nPfoHits
                  << ", matchedHits=" << m_nMatchedHits
                  << ", availableHits=" << m_nTrueWithoutRecoHits << std::endl;

      m_pRecoTree->Fill();
      ++m_index; // Increment index number
    }

    // Match Reco Particles to True Particles
    // ======================================
    MCParticlesToPFParticles matchedParticles;
    MCParticlesToHits matchedParticleHits;
    this->GetRecoToTrueMatches(
      recoParticlesToHits, trueHitsToParticles, matchedParticles, matchedParticleHits);

    // Compare true and reconstructed particles
    for (MCParticlesToHits::const_iterator iter = trueParticlesToHits.begin(),
                                           iterEnd = trueParticlesToHits.end();
         iter != iterEnd;
         ++iter) {
      const art::Ptr<simb::MCParticle> trueParticle = iter->first;
      const HitVector& trueHitVector = iter->second;

      if (trueHitVector.empty()) continue;

      m_mcPdg = trueParticle->PdgCode();
      m_mcNuPdg = 0;
      m_mcParentPdg = 0;
      m_mcPrimaryPdg = 0;
      m_mcIsNeutrino = 0;
      m_mcIsPrimary = 0;
      m_mcIsDecay = 0;
      m_mcIsCC = 0;

      m_pfoPdg = 0;
      m_pfoNuPdg = 0;
      m_pfoParentPdg = 0;
      m_pfoPrimaryPdg = 0;
      m_pfoIsNeutrino = 0;
      m_pfoIsPrimary = 0;
      m_pfoIsStitched = 0;
      m_pfoTrack = 0;
      m_pfoVertex = 0;
      m_pfoVtxX = 0.0;
      m_pfoVtxY = 0.0;
      m_pfoVtxZ = 0.0;
      m_pfoEndX = 0.0;
      m_pfoEndY = 0.0;
      m_pfoEndZ = 0.0;
      m_pfoDirX = 0.0;
      m_pfoDirY = 0.0;
      m_pfoDirZ = 0.0;
      m_pfoLength = 0.0;
      m_pfoStraightLength = 0.0;

      m_mcVertex = 0;
      m_mcVtxX = 0.0;
      m_mcVtxY = 0.0;
      m_mcVtxZ = 0.0;
      m_mcEndX = 0.0;
      m_mcEndY = 0.0;
      m_mcEndZ = 0.0;
      m_mcDirX = 0.0;
      m_mcDirY = 0.0;
      m_mcDirZ = 0.0;
      m_mcEnergy = 0.0;
      m_mcLength = 0.0;
      m_mcStraightLength = 0.0;

      m_completeness = 0.0;
      m_purity = 0.0;

      m_nMCHits = 0;
      m_nMCHitsU = 0;
      m_nMCHitsV = 0;
      m_nMCHitsW = 0;

      m_nPfoHits = 0;
      m_nPfoHitsU = 0;
      m_nPfoHitsV = 0;
      m_nPfoHitsW = 0;

      m_nMatchedHits = 0;
      m_nMatchedHitsU = 0;
      m_nMatchedHitsV = 0;
      m_nMatchedHitsW = 0;

      m_nTrueWithoutRecoHits = 0;
      m_nRecoWithoutTrueHits = 0;

      m_spacepointsMinX = 0.0;
      m_spacepointsMaxX = 0.0;

      // Set true properties
      try {
        int startT(-1);
        int endT(-1);
        this->GetStartAndEndPoints(trueParticle, startT, endT);

        // vertex and end positions
        m_mcVertex = 1;
        m_mcVtxX = trueParticle->Vx(startT);
        m_mcVtxY = trueParticle->Vy(startT);
        m_mcVtxZ = trueParticle->Vz(startT);
        m_mcEndX = trueParticle->Vx(endT);
        m_mcEndY = trueParticle->Vy(endT);
        m_mcEndZ = trueParticle->Vz(endT);

        const double dx(m_mcEndX - m_mcVtxX);
        const double dy(m_mcEndY - m_mcVtxY);
        const double dz(m_mcEndZ - m_mcVtxZ);

        m_mcStraightLength = std::sqrt(dx * dx + dy * dy + dz * dz);
        m_mcLength = this->GetLength(trueParticle, startT, endT);

        // energy and momentum
        const double Ptot(trueParticle->P(startT));

        if (Ptot > 0.0) {
          m_mcDirX = trueParticle->Px(startT) / Ptot;
          m_mcDirY = trueParticle->Py(startT) / Ptot;
          m_mcDirZ = trueParticle->Pz(startT) / Ptot;
          m_mcEnergy = trueParticle->E(startT);
        }
      }
      catch (cet::exception& e) {
      }

      // Get the true parent neutrino
      MCParticlesToMCTruth::const_iterator nuIter = particlesToTruth.find(trueParticle);
      if (particlesToTruth.end() == nuIter)
        throw cet::exception("LArPandora") << " PFParticleMonitoring::analyze --- Found a true "
                                              "particle without any ancestry information ";

      const art::Ptr<simb::MCTruth> trueEvent = nuIter->second;

      if (trueEvent->NeutrinoSet()) {
        const simb::MCNeutrino neutrino(trueEvent->GetNeutrino());
        m_mcNuPdg = neutrino.Nu().PdgCode();
        m_mcIsCC = ((simb::kCC == neutrino.CCNC()) ? 1 : 0);
      }

      // Get the true 'parent' and 'primary' particles
      try {
        const art::Ptr<simb::MCParticle> parentParticle(
          LArPandoraHelper::GetParentMCParticle(trueParticleMap, trueParticle));
        const art::Ptr<simb::MCParticle> primaryParticle(
          LArPandoraHelper::GetFinalStateMCParticle(trueParticleMap, trueParticle));
        m_mcParentPdg = ((parentParticle != trueParticle) ? parentParticle->PdgCode() : 0);
        m_mcPrimaryPdg = primaryParticle->PdgCode();
        m_mcIsPrimary = (primaryParticle == trueParticle);
        m_mcIsDecay = ("Decay" == trueParticle->Process());
      }
      catch (cet::exception& e) {
      }

      // Find min and max X positions of space points
      bool foundSpacePoints(false);

      for (HitVector::const_iterator hIter1 = trueHitVector.begin(),
                                     hIterEnd1 = trueHitVector.end();
           hIter1 != hIterEnd1;
           ++hIter1) {
        const art::Ptr<recob::Hit> hit = *hIter1;

        HitsToSpacePoints::const_iterator hIter2 = hitsToSpacePoints.find(hit);
        if (hitsToSpacePoints.end() == hIter2) continue;

        const art::Ptr<recob::SpacePoint> spacepoint = hIter2->second;
        const double X(spacepoint->XYZ()[0]);

        if (!foundSpacePoints) {
          m_spacepointsMinX = X;
          m_spacepointsMaxX = X;
          foundSpacePoints = true;
        }
        else {
          m_spacepointsMinX = std::min(m_spacepointsMinX, X);
          m_spacepointsMaxX = std::max(m_spacepointsMaxX, X);
        }
      }

      // Count number of available hits
      for (HitVector::const_iterator hIter1 = trueHitVector.begin(),
                                     hIterEnd1 = trueHitVector.end();
           hIter1 != hIterEnd1;
           ++hIter1) {
        if (recoHitsToParticles.find(*hIter1) == recoHitsToParticles.end())
          ++m_nTrueWithoutRecoHits;
      }

      // Match true and reconstructed hits
      m_nMCHits = trueHitVector.size();
      m_nMCHitsU = this->CountHitsByType(geo::kU, trueHitVector);
      m_nMCHitsV = this->CountHitsByType(geo::kV, trueHitVector);
      m_nMCHitsW = this->CountHitsByType(geo::kW, trueHitVector);

      MCParticlesToPFParticles::const_iterator pIter1 = matchedParticles.find(trueParticle);
      if (matchedParticles.end() != pIter1) {
        const art::Ptr<recob::PFParticle> recoParticle = pIter1->second;
        m_pfoPdg = recoParticle->PdgCode();
        m_pfoNuPdg = LArPandoraHelper::GetParentNeutrino(recoParticleMap, recoParticle);
        m_pfoIsPrimary = LArPandoraHelper::IsFinalState(recoParticleMap, recoParticle);

        const art::Ptr<recob::PFParticle> parentParticle =
          LArPandoraHelper::GetParentPFParticle(recoParticleMap, recoParticle);
        m_pfoParentPdg = parentParticle->PdgCode();

        const art::Ptr<recob::PFParticle> primaryParticle =
          LArPandoraHelper::GetFinalStatePFParticle(recoParticleMap, recoParticle);
        m_pfoPrimaryPdg = primaryParticle->PdgCode();

        PFParticlesToHits::const_iterator pIter2 = recoParticlesToHits.find(recoParticle);
        if (recoParticlesToHits.end() == pIter2)
          throw cet::exception("LArPandora")
            << " PFParticleMonitoring::analyze --- Found a reco particle without any hits ";

        const HitVector& recoHitVector = pIter2->second;

        for (HitVector::const_iterator hIter2 = recoHitVector.begin(),
                                       hIterEnd2 = recoHitVector.end();
             hIter2 != hIterEnd2;
             ++hIter2) {
          if (trueHitsToParticles.find(*hIter2) == trueHitsToParticles.end())
            ++m_nRecoWithoutTrueHits;
        }

        MCParticlesToHits::const_iterator pIter3 = matchedParticleHits.find(trueParticle);
        if (matchedParticleHits.end() == pIter3)
          throw cet::exception("LArPandora") << " PFParticleMonitoring::analyze --- Found a "
                                                "matched true particle without matched hits ";

        const HitVector& matchedHitVector = pIter3->second;

        m_nPfoHits = recoHitVector.size();
        m_nPfoHitsU = this->CountHitsByType(geo::kU, recoHitVector);
        m_nPfoHitsV = this->CountHitsByType(geo::kV, recoHitVector);
        m_nPfoHitsW = this->CountHitsByType(geo::kW, recoHitVector);

        m_nMatchedHits = matchedHitVector.size();
        m_nMatchedHitsU = this->CountHitsByType(geo::kU, matchedHitVector);
        m_nMatchedHitsV = this->CountHitsByType(geo::kV, matchedHitVector);
        m_nMatchedHitsW = this->CountHitsByType(geo::kW, matchedHitVector);

        PFParticlesToVertices::const_iterator pIter4 = recoParticlesToVertices.find(recoParticle);
        if (recoParticlesToVertices.end() != pIter4) {
          const VertexVector& vertexVector = pIter4->second;
          if (!vertexVector.empty()) {
            if (vertexVector.size() != 1 && m_printDebug)
              std::cout << " Warning: Found particle with more than one associated vertex "
                        << std::endl;

            const art::Ptr<recob::Vertex> recoVertex = *(vertexVector.begin());
            double xyz[3] = {0.0, 0.0, 0.0};
            recoVertex->XYZ(xyz);

            m_pfoVertex = 1;
            m_pfoVtxX = xyz[0];
            m_pfoVtxY = xyz[1];
            m_pfoVtxZ = xyz[2];
          }
        }

        PFParticlesToTracks::const_iterator pIter5 = recoParticlesToTracks.find(recoParticle);
        if (recoParticlesToTracks.end() != pIter5) {
          const TrackVector& trackVector = pIter5->second;
          if (!trackVector.empty()) {
            if (trackVector.size() != 1 && m_printDebug)
              std::cout << " Warning: Found particle with more than one associated track "
                        << std::endl;

            const art::Ptr<recob::Track> recoTrack = *(trackVector.begin());
            const auto& vtxPosition = recoTrack->Vertex();
            const auto& endPosition = recoTrack->End();
            const auto& vtxDirection = recoTrack->VertexDirection();

            m_pfoTrack = 1;
            m_pfoVtxX = vtxPosition.x();
            m_pfoVtxY = vtxPosition.y();
            m_pfoVtxZ = vtxPosition.z();
            m_pfoEndX = endPosition.x();
            m_pfoEndY = endPosition.y();
            m_pfoEndZ = endPosition.z();
            m_pfoDirX = vtxDirection.x();
            m_pfoDirY = vtxDirection.y();
            m_pfoDirZ = vtxDirection.z();
            m_pfoStraightLength = (endPosition - vtxPosition).R();
            m_pfoLength = recoTrack->Length();
          }
        }

        m_pfoIsStitched = (particlesToT0s.end() != particlesToT0s.find(recoParticle));
      }

      m_purity =
        ((m_nPfoHits == 0) ? 0.0 :
                             static_cast<double>(m_nMatchedHits) / static_cast<double>(m_nPfoHits));
      m_completeness =
        ((m_nPfoHits == 0) ? 0.0 :
                             static_cast<double>(m_nMatchedHits) / static_cast<double>(m_nMCHits));

      if (m_printDebug)
        std::cout << "    MCParticle [" << m_index << "]"
                  << "  trueNu=" << m_mcNuPdg << ", truePdg=" << m_mcPdg
                  << ", recoNu=" << m_pfoNuPdg << ", recoPdg=" << m_pfoPdg
                  << ", mcHits=" << m_nMCHits << ", pfoHits=" << m_nPfoHits
                  << ", matchedHits=" << m_nMatchedHits
                  << ", availableHits=" << m_nTrueWithoutRecoHits << std::endl;

      m_pRecoTree->Fill();
      ++m_index; // Increment index number
    }
  }

void lar_pandora::PFParticleMonitoring::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 332 of file PFParticleMonitoring_module.cc.

References m_completeness, m_event, m_index, m_mcDirX, m_mcDirY, m_mcDirZ, m_mcEndX, m_mcEndY, m_mcEndZ, m_mcEnergy, m_mcIsCC, m_mcIsDecay, m_mcIsNeutrino, m_mcIsPrimary, m_mcLength, m_mcNuPdg, m_mcParentPdg, m_mcPdg, m_mcPrimaryPdg, m_mcStraightLength, m_mcVertex, m_mcVtxX, m_mcVtxY, m_mcVtxZ, m_nDaughterPfos, m_nMatchedHits, m_nMatchedHitsU, m_nMatchedHitsV, m_nMatchedHitsW, m_nMCHits, m_nMCHitsU, m_nMCHitsV, m_nMCHitsW, m_nMCParticles, m_nNeutrinoPfos, m_nPfoHits, m_nPfoHitsU, m_nPfoHitsV, m_nPfoHitsW, m_nPrimaryPfos, m_nRecoWithoutTrueHits, m_nTrueWithoutRecoHits, m_pfoDirX, m_pfoDirY, m_pfoDirZ, m_pfoEndX, m_pfoEndY, m_pfoEndZ, m_pfoIsNeutrino, m_pfoIsPrimary, m_pfoIsStitched, m_pfoLength, m_pfoNuPdg, m_pfoParentPdg, m_pfoPdg, m_pfoPrimaryPdg, m_pfoStraightLength, m_pfoTrack, m_pfoVertex, m_pfoVtxX, m_pfoVtxY, m_pfoVtxZ, m_pRecoTree, m_purity, m_run, m_spacepointsMaxX, and m_spacepointsMinX.

  {
    mf::LogDebug("LArPandora") << " *** PFParticleMonitoring::beginJob() *** " << std::endl;

    //
    art::ServiceHandle<art::TFileService const> tfs;

    m_pRecoTree = tfs->make<TTree>("pandora", "LAr Reco vs True");
    m_pRecoTree->Branch("run", &m_run, "run/I");
    m_pRecoTree->Branch("event", &m_event, "event/I");
    m_pRecoTree->Branch("index", &m_index, "index/I");
    m_pRecoTree->Branch("nMCParticles", &m_nMCParticles, "nMCParticles/I");
    m_pRecoTree->Branch("nNeutrinoPfos", &m_nNeutrinoPfos, "nNeutrinoPfos/I");
    m_pRecoTree->Branch("nPrimaryPfos", &m_nPrimaryPfos, "nPrimaryPfos/I");
    m_pRecoTree->Branch("nDaughterPfos", &m_nDaughterPfos, "nDaughterPfos/I");
    m_pRecoTree->Branch("mcPdg", &m_mcPdg, "mcPdg/I");
    m_pRecoTree->Branch("mcNuPdg", &m_mcNuPdg, "mcNuPdg/I");
    m_pRecoTree->Branch("mcParentPdg", &m_mcParentPdg, "mcParentPdg/I");
    m_pRecoTree->Branch("mcPrimaryPdg", &m_mcPrimaryPdg, "mcPrimaryPdg/I");
    m_pRecoTree->Branch("mcIsNeutrino", &m_mcIsNeutrino, "mcIsNeutrino/I");
    m_pRecoTree->Branch("mcIsPrimary", &m_mcIsPrimary, "mcIsPrimary/I");
    m_pRecoTree->Branch("mcIsDecay", &m_mcIsDecay, "mcIsDecay/I");
    m_pRecoTree->Branch("mcIsCC", &m_mcIsCC, "mcIsCC/I");
    m_pRecoTree->Branch("pfoPdg", &m_pfoPdg, "pfoPdg/I");
    m_pRecoTree->Branch("pfoNuPdg", &m_pfoNuPdg, "pfoNuPdg/I");
    m_pRecoTree->Branch("pfoParentPdg", &m_pfoParentPdg, "pfoParentPdg/I");
    m_pRecoTree->Branch("pfoPrimaryPdg", &m_pfoPrimaryPdg, "pfoPrimaryPdg/I");
    m_pRecoTree->Branch("pfoIsNeutrino", &m_pfoIsNeutrino, "pfoIsNeutrino/I");
    m_pRecoTree->Branch("pfoIsPrimary", &m_pfoIsPrimary, "pfoIsPrimary/I");
    m_pRecoTree->Branch("pfoIsStitched", &m_pfoIsStitched, "pfoIsStitched/I");
    m_pRecoTree->Branch("pfoTrack", &m_pfoTrack, "pfoTrack/I");
    m_pRecoTree->Branch("pfoVertex", &m_pfoVertex, "pfoVertex/I");
    m_pRecoTree->Branch("pfoVtxX", &m_pfoVtxX, "pfoVtxX/D");
    m_pRecoTree->Branch("pfoVtxY", &m_pfoVtxY, "pfoVtxY/D");
    m_pRecoTree->Branch("pfoVtxZ", &m_pfoVtxZ, "pfoVtxZ/D");
    m_pRecoTree->Branch("pfoEndX", &m_pfoEndX, "pfoEndX/D");
    m_pRecoTree->Branch("pfoEndY", &m_pfoEndY, "pfoEndY/D");
    m_pRecoTree->Branch("pfoEndZ", &m_pfoEndZ, "pfoEndZ/D");
    m_pRecoTree->Branch("pfoDirX", &m_pfoDirX, "pfoDirX/D");
    m_pRecoTree->Branch("pfoDirY", &m_pfoDirY, "pfoDirY/D");
    m_pRecoTree->Branch("pfoDirZ", &m_pfoDirZ, "pfoDirZ/D");
    m_pRecoTree->Branch("pfoLength", &m_pfoLength, "pfoLength/D");
    m_pRecoTree->Branch("pfoStraightLength", &m_pfoStraightLength, "pfoStraightLength/D");
    m_pRecoTree->Branch("mcVertex", &m_mcVertex, "mcVertex/I");
    m_pRecoTree->Branch("mcVtxX", &m_mcVtxX, "mcVtxX/D");
    m_pRecoTree->Branch("mcVtxY", &m_mcVtxY, "mcVtxY/D");
    m_pRecoTree->Branch("mcVtxZ", &m_mcVtxZ, "mcVtxZ/D");
    m_pRecoTree->Branch("mcEndX", &m_mcEndX, "mcEndX/D");
    m_pRecoTree->Branch("mcEndY", &m_mcEndY, "mcEndY/D");
    m_pRecoTree->Branch("mcEndZ", &m_mcEndZ, "mcEndZ/D");
    m_pRecoTree->Branch("mcDirX", &m_mcDirX, "mcDirX/D");
    m_pRecoTree->Branch("mcDirY", &m_mcDirY, "mcDirY/D");
    m_pRecoTree->Branch("mcDirZ", &m_mcDirZ, "mcDirZ/D");
    m_pRecoTree->Branch("mcEnergy", &m_mcEnergy, "mcEnergy/D");
    m_pRecoTree->Branch("mcLength", &m_mcLength, "mcLength/D");
    m_pRecoTree->Branch("mcStraightLength", &m_mcStraightLength, "mcStraightLength/D");
    m_pRecoTree->Branch("completeness", &m_completeness, "completeness/D");
    m_pRecoTree->Branch("purity", &m_purity, "purity/D");
    m_pRecoTree->Branch("nMCHits", &m_nMCHits, "nMCHits/I");
    m_pRecoTree->Branch("nPfoHits", &m_nPfoHits, "nPfoHits/I");
    m_pRecoTree->Branch("nMatchedHits", &m_nMatchedHits, "nMatchedHits/I");
    m_pRecoTree->Branch("nMCHitsU", &m_nMCHitsU, "nMCHitsU/I");
    m_pRecoTree->Branch("nMCHitsV", &m_nMCHitsV, "nMCHitsV/I");
    m_pRecoTree->Branch("nMCHitsW", &m_nMCHitsW, "nMCHitsW/I");
    m_pRecoTree->Branch("nPfoHitsU", &m_nPfoHitsU, "nPfoHitsU/I");
    m_pRecoTree->Branch("nPfoHitsV", &m_nPfoHitsV, "nPfoHitsV/I");
    m_pRecoTree->Branch("nPfoHitsW", &m_nPfoHitsW, "nPfoHitsW/I");
    m_pRecoTree->Branch("nMatchedHitsU", &m_nMatchedHitsU, "nMatchedHitsU/I");
    m_pRecoTree->Branch("nMatchedHitsV", &m_nMatchedHitsV, "nMatchedHitsV/I");
    m_pRecoTree->Branch("nMatchedHitsW", &m_nMatchedHitsW, "nMatchedHitsW/I");
    m_pRecoTree->Branch("nTrueWithoutRecoHits", &m_nTrueWithoutRecoHits, "nTrueWithoutRecoHits/I");
    m_pRecoTree->Branch("nRecoWithoutTrueHits", &m_nRecoWithoutTrueHits, "nRecoWithoutTrueHits/I");
    m_pRecoTree->Branch("spacepointsMinX", &m_spacepointsMinX, "spacepointsMinX/D");
    m_pRecoTree->Branch("spacepointsMaxX", &m_spacepointsMaxX, "spacepointsMaxX/D");
  }

void lar_pandora::PFParticleMonitoring::BuildRecoNeutrinoHitMaps ( const PFParticleMap &  recoParticleMap, const PFParticlesToHits &  recoParticlesToHits, PFParticlesToHits &  recoNeutrinosToHits, HitsToPFParticles &  recoHitsToNeutrinos  ) const

private

Build mapping from reconstructed neutrinos to hits.

Parameters

recoParticleMap	the input mapping from reconstructed particle and particle ID
recoParticlesToHits	the input mapping from reconstructed particles to hits
recoNeutrinosToHits	the output mapping from reconstructed particles to hits
recoHitsToNeutrinos	the output mapping from reconstructed hits to particles

Definition at line 1172 of file PFParticleMonitoring_module.cc.

References lar_pandora::LArPandoraHelper::GetParentPFParticle(), and lar_pandora::LArPandoraHelper::IsNeutrino().

Referenced by analyze().

  {
    for (PFParticleMap::const_iterator iter1 = recoParticleMap.begin(),
                                       iterEnd1 = recoParticleMap.end();
         iter1 != iterEnd1;
         ++iter1) {
      const art::Ptr<recob::PFParticle> recoParticle = iter1->second;
      const art::Ptr<recob::PFParticle> recoNeutrino =
        LArPandoraHelper::GetParentPFParticle(recoParticleMap, recoParticle);

      if (!LArPandoraHelper::IsNeutrino(recoNeutrino)) continue;

      const PFParticlesToHits::const_iterator iter2 = recoParticlesToHits.find(recoParticle);
      if (recoParticlesToHits.end() == iter2) continue;

      const HitVector& hitVector = iter2->second;

      for (HitVector::const_iterator iter3 = hitVector.begin(), iterEnd3 = hitVector.end();
           iter3 != iterEnd3;
           ++iter3) {
        const art::Ptr<recob::Hit> hit = *iter3;
        recoHitsToNeutrinos[hit] = recoNeutrino;
        recoNeutrinosToHits[recoNeutrino].push_back(hit);
      }
    }
  }

void lar_pandora::PFParticleMonitoring::BuildTrueNeutrinoHitMaps ( const MCTruthToMCParticles &  truthToParticles, const MCParticlesToHits &  trueParticlesToHits, MCTruthToHits &  trueNeutrinosToHits, HitsToMCTruth &  trueHitsToNeutrinos  ) const

private

Build mapping from true neutrinos to hits.

Parameters

truthToParticles	the input mapping from true event to true particles
trueParticlesToHits	the input mapping from true particles to hits
trueNeutrinosToHits	the output mapping from trues event to hits
trueHitsToNeutrinos	the output mappign from hits to true events

Definition at line 1138 of file PFParticleMonitoring_module.cc.

Referenced by analyze().

  {
    for (MCTruthToMCParticles::const_iterator iter1 = truthToParticles.begin(),
                                              iterEnd1 = truthToParticles.end();
         iter1 != iterEnd1;
         ++iter1) {
      const art::Ptr<simb::MCTruth> trueNeutrino = iter1->first;
      const MCParticleVector& trueParticleVector = iter1->second;

      for (MCParticleVector::const_iterator iter2 = trueParticleVector.begin(),
                                            iterEnd2 = trueParticleVector.end();
           iter2 != iterEnd2;
           ++iter2) {
        const MCParticlesToHits::const_iterator iter3 = trueParticlesToHits.find(*iter2);
        if (trueParticlesToHits.end() == iter3) continue;

        const HitVector& hitVector = iter3->second;

        for (HitVector::const_iterator iter4 = hitVector.begin(), iterEnd4 = hitVector.end();
             iter4 != iterEnd4;
             ++iter4) {
          const art::Ptr<recob::Hit> hit = *iter4;
          trueHitsToNeutrinos[hit] = trueNeutrino;
          trueNeutrinosToHits[trueNeutrino].push_back(hit);
        }
      }
    }
  }

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

int lar_pandora::PFParticleMonitoring::CountHitsByType ( const int  view, const HitVector &  hitVector  ) const

private

Count the number of reconstructed hits in a given wire plane.

Parameters

view	the wire plane ID
hitVector	the input vector of reconstructed hits

Definition at line 1396 of file PFParticleMonitoring_module.cc.

References recob::Hit::View().

Referenced by analyze().

  {
    int nHits(0);

    for (HitVector::const_iterator iter = hitVector.begin(), iterEnd = hitVector.end();
         iter != iterEnd;
         ++iter) {
      const art::Ptr<recob::Hit> hit = *iter;
      if (hit->View() == view) ++nHits;
    }

    return nHits;
  }

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::PFParticleMonitoring::endJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 410 of file PFParticleMonitoring_module.cc.

{}

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

double lar_pandora::PFParticleMonitoring::GetLength ( const art::Ptr< simb::MCParticle >  trueParticle, const int  startT, const int  endT  ) const

private

Find the length of the true particle trajectory through the active region of the detector.

Parameters

trueParticle	the input true particle
startT	the true start point
endT	the true end point

Definition at line 1440 of file PFParticleMonitoring_module.cc.

References simb::MCParticle::Vx(), simb::MCParticle::Vy(), and simb::MCParticle::Vz().

Referenced by analyze().

  {
    if (endT <= startT) return 0.0;

    double length(0.0);

    for (int nt = startT; nt < endT; ++nt) {
      const double dx(particle->Vx(nt + 1) - particle->Vx(nt));
      const double dy(particle->Vy(nt + 1) - particle->Vy(nt));
      const double dz(particle->Vz(nt + 1) - particle->Vz(nt));
      length += sqrt(dx * dx + dy * dy + dz * dz);
    }

    return length;
  }

void lar_pandora::PFParticleMonitoring::GetRecoToTrueMatches ( const PFParticlesToHits &  recoNeutrinosToHits, const HitsToMCTruth &  trueHitsToNeutrinos, MCTruthToPFParticles &  matchedNeutrinos, MCTruthToHits &  matchedNeutrinoHits  ) const

private

Perform matching between true and reconstructed neutrino events.

Parameters

recoNeutrinosToHits	the mapping from reconstructed neutrino events to hits
trueHitsToNeutrinos	the mapping from hits to true neutrino events
matchedNeutrinos	the output matches between reconstructed and true neutrinos
matchedNeutrinoHits	the output matches between reconstructed neutrinos and hits

Definition at line 1204 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and GetRecoToTrueMatches().

  {
    PFParticleSet recoVeto;
    MCTruthSet trueVeto;

    this->GetRecoToTrueMatches(recoNeutrinosToHits,
                               trueHitsToNeutrinos,
                               matchedNeutrinos,
                               matchedNeutrinoHits,
                               recoVeto,
                               trueVeto);
  }

void lar_pandora::PFParticleMonitoring::GetRecoToTrueMatches ( const PFParticlesToHits &  recoNeutrinosToHits, const HitsToMCTruth &  trueHitsToNeutrinos, MCTruthToPFParticles &  matchedNeutrinos, MCTruthToHits &  matchedNeutrinoHits, PFParticleSet &  recoVeto, MCTruthSet &  trueVeto  ) const

private

Perform matching between true and reconstructed neutrino events.

Parameters

recoNeutrinosToHits	the mapping from reconstructed neutrino events to hits
trueHitsToNeutrinos	the mapping from hits to true neutrino events
matchedNeutrinos	the output matches between reconstructed and true neutrinos
matchedNeutrinoHits	the output matches between reconstructed neutrinos and hits
recoVeto	the veto list for reconstructed particles
trueVeto	the veto list for true particles

Definition at line 1222 of file PFParticleMonitoring_module.cc.

References GetRecoToTrueMatches(), and m_recursiveMatching.

  {
    bool foundMatches(false);

    for (PFParticlesToHits::const_iterator iter1 = recoNeutrinosToHits.begin(),
                                           iterEnd1 = recoNeutrinosToHits.end();
         iter1 != iterEnd1;
         ++iter1) {
      const art::Ptr<recob::PFParticle> recoNeutrino = iter1->first;
      if (vetoReco.count(recoNeutrino) > 0) continue;

      const HitVector& hitVector = iter1->second;

      MCTruthToHits truthContributionMap;

      for (HitVector::const_iterator iter2 = hitVector.begin(), iterEnd2 = hitVector.end();
           iter2 != iterEnd2;
           ++iter2) {
        const art::Ptr<recob::Hit> hit = *iter2;

        HitsToMCTruth::const_iterator iter3 = trueHitsToNeutrinos.find(hit);
        if (trueHitsToNeutrinos.end() == iter3) continue;

        const art::Ptr<simb::MCTruth> trueNeutrino = iter3->second;
        if (vetoTrue.count(trueNeutrino) > 0) continue;

        truthContributionMap[trueNeutrino].push_back(hit);
      }

      MCTruthToHits::const_iterator mIter = truthContributionMap.end();

      for (MCTruthToHits::const_iterator iter4 = truthContributionMap.begin(),
                                         iterEnd4 = truthContributionMap.end();
           iter4 != iterEnd4;
           ++iter4) {
        if ((truthContributionMap.end() == mIter) ||
            (iter4->second.size() > mIter->second.size())) {
          mIter = iter4;
        }
      }

      if (truthContributionMap.end() != mIter) {
        const art::Ptr<simb::MCTruth> trueNeutrino = mIter->first;

        MCTruthToHits::const_iterator iter5 = matchedNeutrinoHits.find(trueNeutrino);

        if ((matchedNeutrinoHits.end() == iter5) || (mIter->second.size() > iter5->second.size())) {
          matchedNeutrinos[trueNeutrino] = recoNeutrino;
          matchedNeutrinoHits[trueNeutrino] = mIter->second;
          foundMatches = true;
        }
      }
    }

    if (!foundMatches) return;

    for (MCTruthToPFParticles::const_iterator pIter = matchedNeutrinos.begin(),
                                              pIterEnd = matchedNeutrinos.end();
         pIter != pIterEnd;
         ++pIter) {
      vetoTrue.insert(pIter->first);
      vetoReco.insert(pIter->second);
    }

    if (m_recursiveMatching)
      this->GetRecoToTrueMatches(recoNeutrinosToHits,
                                 trueHitsToNeutrinos,
                                 matchedNeutrinos,
                                 matchedNeutrinoHits,
                                 vetoReco,
                                 vetoTrue);
  }

void lar_pandora::PFParticleMonitoring::GetRecoToTrueMatches ( const PFParticlesToHits &  recoParticlesToHits, const HitsToMCParticles &  trueHitsToParticles, MCParticlesToPFParticles &  matchedParticles, MCParticlesToHits &  matchedHits  ) const

private

Perform matching between true and reconstructed particles.

Parameters

recoParticlesToHits	the mapping from reconstructed particles to hits
trueHitsToParticles	the mapping from hits to true particles
matchedParticles	the output matches between reconstructed and true particles
matchedHits	the output matches between reconstructed particles and hits

void lar_pandora::PFParticleMonitoring::GetRecoToTrueMatches ( const PFParticlesToHits &  recoParticlesToHits, const HitsToMCParticles &  trueHitsToParticles, MCParticlesToPFParticles &  matchedParticles, MCParticlesToHits &  matchedHits, PFParticleSet &  recoVeto, MCParticleSet &  trueVeto  ) const

private

Perform matching between true and reconstructed particles.

Parameters

recoParticlesToHits	the mapping from reconstructed particles to hits
trueHitsToParticles	the mapping from hits to true particles
matchedParticles	the output matches between reconstructed and true particles
matchedHits	the output matches between reconstructed particles and hits
recoVeto	the veto list for reconstructed particles
trueVeto	the veto list for true particles

void lar_pandora::PFParticleMonitoring::GetStartAndEndPoints ( const art::Ptr< simb::MCParticle >  trueParticle, int &  startT, int &  endT  ) const

private

Find the start and end points of the true particle in the active region of detector.

Parameters

trueParticle	the input true particle
startT	the true start point
endT	the true end point

Definition at line 1412 of file PFParticleMonitoring_module.cc.

References simb::MCParticle::NumberTrajectoryPoints(), geo::GeometryCore::PositionToTPCptr(), simb::MCParticle::Vx(), simb::MCParticle::Vy(), and simb::MCParticle::Vz().

Referenced by analyze().

  {
    art::ServiceHandle<geo::Geometry const> theGeometry;

    bool foundStartPosition(false);

    const int numTrajectoryPoints(static_cast<int>(particle->NumberTrajectoryPoints()));

    for (int nt = 0; nt < numTrajectoryPoints; ++nt) {
      geo::Point_t const pos{particle->Vx(nt), particle->Vy(nt), particle->Vz(nt)};
      if (theGeometry->PositionToTPCptr(pos) == nullptr) { continue; }

      // TODO: Apply fiducial cut due to readout window

      endT = nt;
      if (!foundStartPosition) {
        startT = endT;
        foundStartPosition = true;
      }
    }

    if (!foundStartPosition) throw cet::exception("LArPandora");
  }

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::PFParticleMonitoring::reconfigure

(

fhicl::ParameterSet const &

pset

)

Definition at line 312 of file PFParticleMonitoring_module.cc.

References fhicl::ParameterSet::get(), m_addDaughterMCParticles, m_addDaughterPFParticles, m_backtrackerLabel, m_disableRealDataCheck, m_geantModuleLabel, m_hitfinderLabel, m_particleLabel, m_printDebug, m_recursiveMatching, m_trackLabel, m_useDaughterMCParticles, and m_useDaughterPFParticles.

Referenced by PFParticleMonitoring().

  {
    m_trackLabel = pset.get<std::string>("TrackModule", "pandoraTracks");
    m_particleLabel = pset.get<std::string>("PFParticleModule", "pandora");
    m_hitfinderLabel = pset.get<std::string>("HitFinderModule", "gaushit");
    m_backtrackerLabel = pset.get<std::string>("BackTrackerModule", "gaushitTruthMatch");
    m_geantModuleLabel = pset.get<std::string>("GeantModule", "largeant");

    m_useDaughterPFParticles = pset.get<bool>("UseDaughterPFParticles", false);
    m_useDaughterMCParticles = pset.get<bool>("UseDaughterMCParticles", true);
    m_addDaughterPFParticles = pset.get<bool>("AddDaughterPFParticles", true);
    m_addDaughterMCParticles = pset.get<bool>("AddDaughterMCParticles", true);

    m_recursiveMatching = pset.get<bool>("RecursiveMatching", false);
    m_printDebug = pset.get<bool>("PrintDebug", false);
    m_disableRealDataCheck = pset.get<bool>("DisableRealDataCheck", false);
  }

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

Member Data Documentation

bool lar_pandora::PFParticleMonitoring::m_addDaughterMCParticles

private

Definition at line 256 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

bool lar_pandora::PFParticleMonitoring::m_addDaughterPFParticles

private

Definition at line 255 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

std::string lar_pandora::PFParticleMonitoring::m_backtrackerLabel

private

Definition at line 250 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

double lar_pandora::PFParticleMonitoring::m_completeness

private

Definition at line 220 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

bool lar_pandora::PFParticleMonitoring::m_disableRealDataCheck

private

Whether to check if the input file contains real data before accessing MC information.

Definition at line 261 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

int lar_pandora::PFParticleMonitoring::m_event

private

Definition at line 167 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

std::string lar_pandora::PFParticleMonitoring::m_geantModuleLabel

private

Definition at line 251 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

std::string lar_pandora::PFParticleMonitoring::m_hitfinderLabel

private

Definition at line 247 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

int lar_pandora::PFParticleMonitoring::m_index

private

Definition at line 168 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcDirX

private

Definition at line 213 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcDirY

private

Definition at line 214 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcDirZ

private

Definition at line 215 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcEndX

private

Definition at line 210 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcEndY

private

Definition at line 211 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcEndZ

private

Definition at line 212 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcEnergy

private

Definition at line 216 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcIsCC

private

Definition at line 182 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcIsDecay

private

Definition at line 181 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcIsNeutrino

private

Definition at line 179 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcIsPrimary

private

Definition at line 180 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcLength

private

Definition at line 217 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcNuPdg

private

Definition at line 176 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcParentPdg

private

Definition at line 177 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcPdg

private

Definition at line 175 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcPrimaryPdg

private

Definition at line 178 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcStraightLength

private

Definition at line 218 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcVertex

private

Definition at line 206 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcVtxX

private

Definition at line 207 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcVtxY

private

Definition at line 208 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcVtxZ

private

Definition at line 209 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nDaughterPfos

private

Definition at line 173 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMatchedHits

private

Definition at line 225 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMatchedHitsU

private

Definition at line 235 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMatchedHitsV

private

Definition at line 236 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMatchedHitsW

private

Definition at line 237 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMCHits

private

Definition at line 223 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMCHitsU

private

Definition at line 227 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMCHitsV

private

Definition at line 228 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMCHitsW

private

Definition at line 229 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMCParticles

private

Definition at line 170 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nNeutrinoPfos

private

Definition at line 171 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nPfoHits

private

Definition at line 224 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nPfoHitsU

private

Definition at line 231 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nPfoHitsV

private

Definition at line 232 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nPfoHitsW

private

Definition at line 233 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nPrimaryPfos

private

Definition at line 172 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nRecoWithoutTrueHits

private

Reconstructed hits which don't belong to any true particle - "missing".

Definition at line 242 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nTrueWithoutRecoHits

private

True hits which don't belong to any reconstructed particle - "available".

Definition at line 240 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

std::string lar_pandora::PFParticleMonitoring::m_particleLabel

private

Definition at line 249 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

double lar_pandora::PFParticleMonitoring::m_pfoDirX

private

Definition at line 200 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoDirY

private

Definition at line 201 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoDirZ

private

Definition at line 202 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoEndX

private

Definition at line 197 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoEndY

private

Definition at line 198 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoEndZ

private

Definition at line 199 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoIsNeutrino

private

Definition at line 188 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoIsPrimary

private

Definition at line 189 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoIsStitched

private

Definition at line 190 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoLength

private

Definition at line 203 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoNuPdg

private

Definition at line 185 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoParentPdg

private

Definition at line 186 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoPdg

private

Definition at line 184 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoPrimaryPdg

private

Definition at line 187 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoStraightLength

private

Definition at line 204 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoTrack

private

Definition at line 192 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoVertex

private

Definition at line 193 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoVtxX

private

Definition at line 194 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoVtxY

private

Definition at line 195 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoVtxZ

private

Definition at line 196 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

TTree* lar_pandora::PFParticleMonitoring::m_pRecoTree

private

Definition at line 164 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

bool lar_pandora::PFParticleMonitoring::m_printDebug

private

switch for print statements (TODO: use message service!)

Definition at line 259 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

double lar_pandora::PFParticleMonitoring::m_purity

private

Definition at line 221 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

bool lar_pandora::PFParticleMonitoring::m_recursiveMatching

private

Definition at line 258 of file PFParticleMonitoring_module.cc.

Referenced by GetRecoToTrueMatches(), and reconfigure().

int lar_pandora::PFParticleMonitoring::m_run

private

Definition at line 166 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_spacepointsMaxX

private

Definition at line 245 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_spacepointsMinX

private

Definition at line 244 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

std::string lar_pandora::PFParticleMonitoring::m_trackLabel

private

Definition at line 248 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

bool lar_pandora::PFParticleMonitoring::m_useDaughterMCParticles

private

Definition at line 254 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

bool lar_pandora::PFParticleMonitoring::m_useDaughterPFParticles

private

Definition at line 253 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

---

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

• larpandora/v09_22_05/source/larpandora/LArPandoraAnalysis/PFParticleMonitoring_module.cc