lar_pandora::PFParticleMonitoring Class Reference

PFParticleMonitoring class. More...

Inheritance diagram for lar_pandora::PFParticleMonitoring:

Public Types
using	WorkerType = WorkerT< EDAnalyzer >
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)
std::string	workerType () const
bool	modifiesEvent () const
void	registerProducts (MasterProductRegistry &, ProductDescriptions &, ModuleDescription const &)
std::string const &	processName () const
bool	wantAllEvents () const
bool	wantEvent (Event const &e)
fhicl::ParameterSetID	selectorConfig () const
art::Handle< art::TriggerResults >	getTriggerResults (Event const &e) const
template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)
template<typename T , art::BranchType BT>
art::ProductToken< T >	consumes (InputTag const &it)
template<typename T , BranchType = InEvent>
void	consumesMany ()
template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)
template<typename T , art::BranchType BT>
art::ViewToken< T >	consumesView (InputTag const &it)
template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)
template<typename T , art::BranchType BT>
art::ProductToken< T >	mayConsume (InputTag const &it)
template<typename T , BranchType = InEvent>
void	mayConsumeMany ()
template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)
template<typename T , art::BranchType BT>
art::ViewToken< T >	mayConsumeView (InputTag const &it)
base_engine_t &	createEngine (seed_t seed)
base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make)
base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make, label_t const &engine_label)
seed_t	get_seed_value (fhicl::ParameterSet const &pset, char const key[]="seed", seed_t const implicit_seed=-1)

Static Public Member Functions
static cet::exempt_ptr< Consumer >	non_module_context ()

Protected Member Functions
CurrentProcessingContext const *	currentContext () const
detail::CachedProducts &	cachedProducts ()
void	validateConsumedProduct (BranchType const bt, ProductInfo const &pi)
void	prepareForJob (fhicl::ParameterSet const &pset)
void	showMissingConsumes () const

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

Detailed Description

PFParticleMonitoring class.

Definition at line 25 of file PFParticleMonitoring_module.cc.

Member Typedef Documentation

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

private

Definition at line 48 of file PFParticleMonitoring_module.cc.

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

private

Definition at line 49 of file PFParticleMonitoring_module.cc.

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 39 of file EDAnalyzer.h.

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

private

Definition at line 47 of file PFParticleMonitoring_module.cc.

using art::EDAnalyzer::WorkerType = WorkerT<EDAnalyzer>

inherited

Definition at line 38 of file EDAnalyzer.h.

Constructor & Destructor Documentation

lar_pandora::PFParticleMonitoring::PFParticleMonitoring

(

fhicl::ParameterSet const &

pset

)

Constructor.

Parameters

pset

Definition at line 281 of file PFParticleMonitoring_module.cc.

References reconfigure().

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

lar_pandora::PFParticleMonitoring::~PFParticleMonitoring ( )

virtual

Destructor.

Definition at line 288 of file PFParticleMonitoring_module.cc.

{
}

Member Function Documentation

void lar_pandora::PFParticleMonitoring::analyze

(

const art::Event &

evt

)

Definition at line 397 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_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, max, min, simb::MCTruth::NeutrinoSet(), simb::MCParticle::P(), recob::PFParticle::PdgCode(), simb::MCParticle::PdgCode(), simb::MCParticle::Process(), simb::MCParticle::Px(), simb::MCParticle::Py(), simb::MCParticle::Pz(), R, 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 (!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 313 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, m_spacepointsMinX, and art::TFileDirectory::make().

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

    //
    art::ServiceHandle<art::TFileService> 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 1119 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 1090 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);
            }
        }
    }
}

detail::CachedProducts& art::EventObserverBase::cachedProducts ( )

inlineprotectedinherited

Definition at line 79 of file EventObserverBase.h.

References art::EventObserverBase::selectors_.

Referenced by art::EDAnalyzer::doEvent(), and art::OutputModule::doWriteEvent().

  {
    return selectors_;
  }

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::consumes ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::consumes ( InputTag const &  it )

inherited

Definition at line 147 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

{
  if (!moduleContext_)
    return ProductToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::Product,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ProductToken<T>{it};
}

template<typename T , art::BranchType BT>

void art::Consumer::consumesMany ( )

inherited

Definition at line 162 of file Consumer.h.

{
  if (!moduleContext_)
    return;

  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
}

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::consumesView ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::consumesView ( InputTag const &  it )

inherited

Definition at line 172 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

{
  if (!moduleContext_)
    return ViewToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ViewToken<T>{it};
}

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

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed )

inherited

Definition at line 26 of file EngineCreator.cc.

References art::EngineCreator::rng().

Referenced by evgen::CosmicsGen::CosmicsGen(), rndm::NuRandomService::createEngine(), cluster::fuzzyCluster::fuzzyCluster(), cluster::HoughLineFinder::HoughLineFinder(), art::MixFilter< T >::initEngine_(), larg4::LArG4::LArG4(), evgen::LightSource::LightSource(), evgen::NeutronOsc::NeutronOsc(), evgen::NucleonDecay::NucleonDecay(), opdet::OpMCDigi::OpMCDigi(), opdet::OptDetDigitizer::OptDetDigitizer(), phot::PhotonLibraryPropagation::PhotonLibraryPropagation(), detsim::SimDriftElectrons::SimDriftElectrons(), evgen::SingleGen::SingleGen(), evgen::SNNueAr40CCGen::SNNueAr40CCGen(), ToyOneShowerGen::ToyOneShowerGen(), and trkf::Track3DKalman::Track3DKalman().

{
  return rng()->createEngine(placeholder_schedule_id(), seed);
}

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed, std::string const &  kind_of_engine_to_make  )

inherited

Definition at line 32 of file EngineCreator.cc.

References art::EngineCreator::rng().

{
  return rng()->createEngine(
    placeholder_schedule_id(), seed, kind_of_engine_to_make);
}

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed, std::string const &  kind_of_engine_to_make, label_t const &  engine_label  )

inherited

Definition at line 40 of file EngineCreator.cc.

References art::EngineCreator::rng().

{
  return rng()->createEngine(
    placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
}

CurrentProcessingContext const * art::EDAnalyzer::currentContext ( ) const

protectedinherited

Definition at line 114 of file EDAnalyzer.cc.

References art::EDAnalyzer::current_context_.

Referenced by art::EDAnalyzer::workerType().

  {
    return current_context_.get();
  }

void lar_pandora::PFParticleMonitoring::endJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 391 of file PFParticleMonitoring_module.cc.

{
}

EngineCreator::seed_t EngineCreator::get_seed_value ( fhicl::ParameterSet const &  pset, char const  key[] = "seed", seed_t const  implicit_seed = -1  )

inherited

Definition at line 49 of file EngineCreator.cc.

References fhicl::ParameterSet::get().

Referenced by art::MixFilter< T >::initEngine_().

{
  auto const& explicit_seeds = pset.get<std::vector<int>>(key, {});
  return explicit_seeds.empty() ? implicit_seed : explicit_seeds.front();
}

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 1366 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 1147 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 1157 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 1329 of file PFParticleMonitoring_module.cc.

References e, max, simb::MCParticle::NumberTrajectoryPoints(), geo::GeometryCore::PositionToTPC(), simb::MCParticle::Vx(), simb::MCParticle::Vy(), and simb::MCParticle::Vz().

Referenced by analyze().

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

    bool foundStartPosition(false);

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

    for (int nt = 0; nt < numTrajectoryPoints; ++nt)
    {
        try
        {
            double pos[3] = {particle->Vx(nt), particle->Vy(nt), particle->Vz(nt)};
            unsigned int which_tpc(std::numeric_limits<unsigned int>::max());
            unsigned int which_cstat(std::numeric_limits<unsigned int>::max());
            theGeometry->PositionToTPC(pos, which_tpc, which_cstat);

            // TODO: Apply fiducial cut due to readout window

            endT = nt;
            if (!foundStartPosition)
            {
                startT = endT;
                foundStartPosition = true;
            }
        }
        catch (cet::exception &e){
            continue;
        }
    }

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

art::Handle<art::TriggerResults> art::EventObserverBase::getTriggerResults ( Event const &  e ) const

inlineinherited

Definition at line 61 of file EventObserverBase.h.

References art::detail::CachedProducts::getOneTriggerResults(), and art::EventObserverBase::selectors_.

Referenced by art::OutputModule::doWriteEvent().

  {
    return selectors_.getOneTriggerResults(e);
  }

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::mayConsume ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::mayConsume ( InputTag const &  it )

inherited

Definition at line 190 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

{
  if (!moduleContext_)
    return ProductToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::Product,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ProductToken<T>{it};
}

template<typename T , art::BranchType BT>

void art::Consumer::mayConsumeMany ( )

inherited

Definition at line 205 of file Consumer.h.

{
  if (!moduleContext_)
    return;

  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
}

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::mayConsumeView ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::mayConsumeView ( InputTag const &  it )

inherited

Definition at line 215 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

{
  if (!moduleContext_)
    return ViewToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ViewToken<T>{it};
}

bool art::EventObserverBase::modifiesEvent ( ) const

inlineinherited

Definition at line 25 of file EventObserverBase.h.

  {
    return false;
  }

cet::exempt_ptr< art::Consumer > art::Consumer::non_module_context ( )

staticinherited

Definition at line 76 of file Consumer.cc.

Referenced by art::RootOutput::beginSubRun(), art::OutputModule::doBeginRun(), art::OutputModule::doBeginSubRun(), art::OutputModule::doEndRun(), art::OutputModule::doEndSubRun(), art::ProducingService::doPostReadEvent(), art::ProducingService::doPostReadRun(), art::ProducingService::doPostReadSubRun(), art::OutputModule::doWriteEvent(), art::ProcessPackage< L >::postScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::Run >::End::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::End::postScheduleSignal(), art::ProcessPackage< L >::preScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::preScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::preScheduleSignal(), art::EventProcessor::readEvent(), art::EventProcessor::readRun(), art::EmptyEvent::readRun_(), art::EventProcessor::readSubRun(), and art::EmptyEvent::readSubRun_().

{
  static Consumer invalid{InvalidTag{}};
  return &invalid;
}

void art::Consumer::prepareForJob ( fhicl::ParameterSet const &  pset )

protectedinherited

Definition at line 89 of file Consumer.cc.

References fhicl::ParameterSet::get_if_present().

Referenced by art::EDProducer::doBeginJob(), art::EDFilter::doBeginJob(), and art::EDAnalyzer::doBeginJob().

{
  if (!moduleContext_)
    return;

  pset.get_if_present("errorOnMissingConsumes", requireConsumes_);
  for (auto& consumablesPerBranch : consumables_) {
    cet::sort_all(consumablesPerBranch);
  }
}

std::string const& art::EventObserverBase::processName ( ) const

inlineinherited

Definition at line 41 of file EventObserverBase.h.

References art::EventObserverBase::process_name_.

Referenced by art::FileDumperOutput::printPrincipal(), and art::RootOutput::RootOutput().

  {
    return process_name_;
  }

void lar_pandora::PFParticleMonitoring::reconfigure

(

fhicl::ParameterSet const &

pset

)

Definition at line 294 of file PFParticleMonitoring_module.cc.

References fhicl::ParameterSet::get(), m_addDaughterMCParticles, m_addDaughterPFParticles, m_backtrackerLabel, 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);
}

void art::EventObserverBase::registerProducts ( MasterProductRegistry &  , ProductDescriptions &  , ModuleDescription const &    )

inlineinherited

Definition at line 33 of file EventObserverBase.h.

  {}

fhicl::ParameterSetID art::EventObserverBase::selectorConfig ( ) const

inlineinherited

Definition at line 56 of file EventObserverBase.h.

References art::EventObserverBase::selector_config_id_.

Referenced by art::RootOutputFile::writeOne().

  {
    return selector_config_id_;
  }

void art::Consumer::showMissingConsumes ( ) const

protectedinherited

Definition at line 125 of file Consumer.cc.

Referenced by art::EDProducer::doEndJob(), art::EDFilter::doEndJob(), art::EDAnalyzer::doEndJob(), and art::RootOutput::endJob().

{
  if (!moduleContext_)
    return;

  // If none of the branches have missing consumes statements, exit early.
  if (std::all_of(cbegin(missingConsumes_),
                  cend(missingConsumes_),
                  [](auto const& perBranch) { return perBranch.empty(); }))
    return;

  constexpr cet::HorizontalRule rule{60};
  mf::LogPrint log{"MTdiagnostics"};
  log << '\n'
      << rule('=') << '\n'
      << "The following consumes (or mayConsume) statements are missing from\n"
      << module_context(moduleDescription_) << '\n'
      << rule('-') << '\n';

  cet::for_all_with_index(
    missingConsumes_, [&log](std::size_t const i, auto const& perBranch) {
      for (auto const& pi : perBranch) {
        log << "  "
            << assemble_consumes_statement(static_cast<BranchType>(i), pi)
            << '\n';
      }
    });
  log << rule('=');
}

void art::Consumer::validateConsumedProduct ( BranchType const  bt, ProductInfo const &  pi  )

protectedinherited

Definition at line 101 of file Consumer.cc.

References art::errors::ProductRegistrationFailure.

{
  // Early exits if consumes tracking has been disabled or if the
  // consumed product is an allowed consumable.
  if (!moduleContext_)
    return;

  if (cet::binary_search_all(consumables_[bt], pi))
    return;

  if (requireConsumes_) {
    throw Exception(errors::ProductRegistrationFailure,
                    "Consumer: an error occurred during validation of a "
                    "retrieved product\n\n")
      << "The following consumes (or mayConsume) statement is missing from\n"
      << module_context(moduleDescription_) << ":\n\n"
      << "  " << assemble_consumes_statement(bt, pi) << "\n\n";
  }

  missingConsumes_[bt].insert(pi);
}

bool art::EventObserverBase::wantAllEvents ( ) const

inlineinherited

Definition at line 46 of file EventObserverBase.h.

References art::EventObserverBase::wantAllEvents_.

Referenced by art::EDAnalyzer::doEvent(), art::OutputModule::doEvent(), art::OutputModule::doWriteEvent(), art::RootOutput::RootOutput(), and art::OutputWorker::wantAllEvents().

  {
    return wantAllEvents_;
  }

bool art::EventObserverBase::wantEvent ( Event const &  e )

inlineinherited

Definition at line 51 of file EventObserverBase.h.

References art::EventObserverBase::selectors_, and art::detail::CachedProducts::wantEvent().

Referenced by art::EDAnalyzer::doEvent(), art::OutputModule::doEvent(), and art::OutputModule::doWriteEvent().

  {
    return selectors_.wantEvent(e);
  }

std::string art::EDAnalyzer::workerType ( ) const

inlineinherited

Definition at line 109 of file EDAnalyzer.h.

References art::EDAnalyzer::currentContext().

    {
      return "WorkerT<EDAnalyzer>";
    }

Member Data Documentation

bool lar_pandora::PFParticleMonitoring::m_addDaughterMCParticles

private

Definition at line 238 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

bool lar_pandora::PFParticleMonitoring::m_addDaughterPFParticles

private

Definition at line 237 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

std::string lar_pandora::PFParticleMonitoring::m_backtrackerLabel

private

Definition at line 232 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

double lar_pandora::PFParticleMonitoring::m_completeness

private

Definition at line 204 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_event

private

Definition at line 151 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

std::string lar_pandora::PFParticleMonitoring::m_geantModuleLabel

private

Definition at line 233 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

std::string lar_pandora::PFParticleMonitoring::m_hitfinderLabel

private

Definition at line 229 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

int lar_pandora::PFParticleMonitoring::m_index

private

Definition at line 152 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcDirX

private

Definition at line 197 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcDirY

private

Definition at line 198 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcDirZ

private

Definition at line 199 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcEndX

private

Definition at line 194 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcEndY

private

Definition at line 195 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcEndZ

private

Definition at line 196 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcEnergy

private

Definition at line 200 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcIsCC

private

Definition at line 166 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcIsDecay

private

Definition at line 165 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcIsNeutrino

private

Definition at line 163 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcIsPrimary

private

Definition at line 164 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcLength

private

Definition at line 201 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcNuPdg

private

Definition at line 160 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcParentPdg

private

Definition at line 161 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcPdg

private

Definition at line 159 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcPrimaryPdg

private

Definition at line 162 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcStraightLength

private

Definition at line 202 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_mcVertex

private

Definition at line 190 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcVtxX

private

Definition at line 191 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcVtxY

private

Definition at line 192 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_mcVtxZ

private

Definition at line 193 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nDaughterPfos

private

Definition at line 157 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMatchedHits

private

Definition at line 209 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMatchedHitsU

private

Definition at line 219 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMatchedHitsV

private

Definition at line 220 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMatchedHitsW

private

Definition at line 221 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMCHits

private

Definition at line 207 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMCHitsU

private

Definition at line 211 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMCHitsV

private

Definition at line 212 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMCHitsW

private

Definition at line 213 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nMCParticles

private

Definition at line 154 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nNeutrinoPfos

private

Definition at line 155 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nPfoHits

private

Definition at line 208 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nPfoHitsU

private

Definition at line 215 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nPfoHitsV

private

Definition at line 216 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nPfoHitsW

private

Definition at line 217 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_nPrimaryPfos

private

Definition at line 156 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 224 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 223 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

std::string lar_pandora::PFParticleMonitoring::m_particleLabel

private

Definition at line 231 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

double lar_pandora::PFParticleMonitoring::m_pfoDirX

private

Definition at line 184 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoDirY

private

Definition at line 185 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoDirZ

private

Definition at line 186 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoEndX

private

Definition at line 181 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoEndY

private

Definition at line 182 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoEndZ

private

Definition at line 183 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoIsNeutrino

private

Definition at line 172 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoIsPrimary

private

Definition at line 173 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoIsStitched

private

Definition at line 174 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoLength

private

Definition at line 187 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoNuPdg

private

Definition at line 169 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoParentPdg

private

Definition at line 170 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoPdg

private

Definition at line 168 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoPrimaryPdg

private

Definition at line 171 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoStraightLength

private

Definition at line 188 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoTrack

private

Definition at line 176 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

int lar_pandora::PFParticleMonitoring::m_pfoVertex

private

Definition at line 177 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoVtxX

private

Definition at line 178 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoVtxY

private

Definition at line 179 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_pfoVtxZ

private

Definition at line 180 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

TTree* lar_pandora::PFParticleMonitoring::m_pRecoTree

private

Definition at line 148 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 241 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

double lar_pandora::PFParticleMonitoring::m_purity

private

Definition at line 205 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

bool lar_pandora::PFParticleMonitoring::m_recursiveMatching

private

Definition at line 240 of file PFParticleMonitoring_module.cc.

Referenced by GetRecoToTrueMatches(), and reconfigure().

int lar_pandora::PFParticleMonitoring::m_run

private

Definition at line 150 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_spacepointsMaxX

private

Definition at line 227 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

double lar_pandora::PFParticleMonitoring::m_spacepointsMinX

private

Definition at line 226 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and beginJob().

std::string lar_pandora::PFParticleMonitoring::m_trackLabel

private

Definition at line 230 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

bool lar_pandora::PFParticleMonitoring::m_useDaughterMCParticles

private

Definition at line 236 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

bool lar_pandora::PFParticleMonitoring::m_useDaughterPFParticles

private

Definition at line 235 of file PFParticleMonitoring_module.cc.

Referenced by analyze(), and reconfigure().

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/larpandora/v07_06_01/source/larpandora/LArPandoraAnalysis/PFParticleMonitoring_module.cc