Inheritance diagram for lar_pandora::PFParticleMonitoring:

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)

Public Member Functions inherited from art::EDAnalyzer
	EDAnalyzer (fhicl::ParameterSet const &pset)

template<typename Config >
	EDAnalyzer (Table< Config > const &config)

std::string	workerType () const

Public Member Functions inherited from art::detail::Analyzer
virtual	~Analyzer () noexcept

	Analyzer (fhicl::ParameterSet const &pset)

template<typename Config >
	Analyzer (Table< Config > const &config)

void	doBeginJob (SharedResources const &resources)

void	doEndJob ()

void	doRespondToOpenInputFile (FileBlock const &fb)

void	doRespondToCloseInputFile (FileBlock const &fb)

void	doRespondToOpenOutputFiles (FileBlock const &fb)

void	doRespondToCloseOutputFiles (FileBlock const &fb)

bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)

Public Member Functions inherited from art::Observer
	~Observer () noexcept

	Observer (Observer const &)=delete

	Observer (Observer &&)=delete

Observer &	operator= (Observer const &)=delete

Observer &	operator= (Observer &&)=delete

void	registerProducts (ProductDescriptions &, ModuleDescription const &)

void	fillDescriptions (ModuleDescription const &)

fhicl::ParameterSetID	selectorConfig () const

Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept

	ModuleBase ()

ModuleDescription const &	moduleDescription () const

void	setModuleDescription (ModuleDescription const &)

std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const

void	sortConsumables (std::string const &current_process_name)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

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

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >

using	ModuleType = EDAnalyzer

Protected Member Functions inherited from art::Observer
std::string const &	processName () const

bool	wantAllEvents () const noexcept

bool	wantEvent (ScheduleID id, Event const &e) const

Handle< TriggerResults >	getTriggerResults (Event const &e) const

	Observer (fhicl::ParameterSet const &config)

	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)

Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

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.

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.

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

lar_pandora::PFParticleMonitoring::~PFParticleMonitoring ( )

virtual

Destructor.

Definition at line 308 of file PFParticleMonitoring_module.cc.

308 {}

Member Function Documentation

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

Definition at line 419 of file PFParticleMonitoring_module.cc.

   {
     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 334 of file PFParticleMonitoring_module.cc.

   {
     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 1179 of file PFParticleMonitoring_module.cc.

   {
     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 1144 of file PFParticleMonitoring_module.cc.

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

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 1408 of file PFParticleMonitoring_module.cc.

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

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 413 of file PFParticleMonitoring_module.cc.

414 {}

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 1461 of file PFParticleMonitoring_module.cc.

   {
     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 1212 of file PFParticleMonitoring_module.cc.

   {
     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 1231 of file PFParticleMonitoring_module.cc.

   {
     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 1425 of file PFParticleMonitoring_module.cc.

   {
     art::ServiceHandle<geo::Geometry const> 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");
   }

void lar_pandora::PFParticleMonitoring::reconfigure ( fhicl::ParameterSet const & pset )

Definition at line 313 of file PFParticleMonitoring_module.cc.

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

Member Data Documentation

bool lar_pandora::PFParticleMonitoring::m_addDaughterMCParticles

private

Definition at line 256 of file PFParticleMonitoring_module.cc.

bool lar_pandora::PFParticleMonitoring::m_addDaughterPFParticles

private

Definition at line 255 of file PFParticleMonitoring_module.cc.

std::string lar_pandora::PFParticleMonitoring::m_backtrackerLabel

private

Definition at line 250 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_completeness

private

Definition at line 220 of file PFParticleMonitoring_module.cc.

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.

int lar_pandora::PFParticleMonitoring::m_event

private

Definition at line 167 of file PFParticleMonitoring_module.cc.

std::string lar_pandora::PFParticleMonitoring::m_geantModuleLabel

private

Definition at line 251 of file PFParticleMonitoring_module.cc.

std::string lar_pandora::PFParticleMonitoring::m_hitfinderLabel

private

Definition at line 247 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_index

private

Definition at line 168 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_mcDirX

private

Definition at line 213 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_mcDirY

private

Definition at line 214 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_mcDirZ

private

Definition at line 215 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_mcEndX

private

Definition at line 210 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_mcEndY

private

Definition at line 211 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_mcEndZ

private

Definition at line 212 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_mcEnergy

private

Definition at line 216 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_mcIsCC

private

Definition at line 182 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_mcIsDecay

private

Definition at line 181 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_mcIsNeutrino

private

Definition at line 179 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_mcIsPrimary

private

Definition at line 180 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_mcLength

private

Definition at line 217 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_mcNuPdg

private

Definition at line 176 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_mcParentPdg

private

Definition at line 177 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_mcPdg

private

Definition at line 175 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_mcPrimaryPdg

private

Definition at line 178 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_mcStraightLength

private

Definition at line 218 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_mcVertex

private

Definition at line 206 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_mcVtxX

private

Definition at line 207 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_mcVtxY

private

Definition at line 208 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_mcVtxZ

private

Definition at line 209 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_nDaughterPfos

private

Definition at line 173 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_nMatchedHits

private

Definition at line 225 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_nMatchedHitsU

private

Definition at line 235 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_nMatchedHitsV

private

Definition at line 236 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_nMatchedHitsW

private

Definition at line 237 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_nMCHits

private

Definition at line 223 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_nMCHitsU

private

Definition at line 227 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_nMCHitsV

private

Definition at line 228 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_nMCHitsW

private

Definition at line 229 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_nMCParticles

private

Definition at line 170 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_nNeutrinoPfos

private

Definition at line 171 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_nPfoHits

private

Definition at line 224 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_nPfoHitsU

private

Definition at line 231 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_nPfoHitsV

private

Definition at line 232 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_nPfoHitsW

private

Definition at line 233 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_nPrimaryPfos

private

Definition at line 172 of file PFParticleMonitoring_module.cc.

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.

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.

std::string lar_pandora::PFParticleMonitoring::m_particleLabel

private

Definition at line 249 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_pfoDirX

private

Definition at line 200 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_pfoDirY

private

Definition at line 201 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_pfoDirZ

private

Definition at line 202 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_pfoEndX

private

Definition at line 197 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_pfoEndY

private

Definition at line 198 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_pfoEndZ

private

Definition at line 199 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_pfoIsNeutrino

private

Definition at line 188 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_pfoIsPrimary

private

Definition at line 189 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_pfoIsStitched

private

Definition at line 190 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_pfoLength

private

Definition at line 203 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_pfoNuPdg

private

Definition at line 185 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_pfoParentPdg

private

Definition at line 186 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_pfoPdg

private

Definition at line 184 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_pfoPrimaryPdg

private

Definition at line 187 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_pfoStraightLength

private

Definition at line 204 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_pfoTrack

private

Definition at line 192 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_pfoVertex

private

Definition at line 193 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_pfoVtxX

private

Definition at line 194 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_pfoVtxY

private

Definition at line 195 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_pfoVtxZ

private

Definition at line 196 of file PFParticleMonitoring_module.cc.

TTree* lar_pandora::PFParticleMonitoring::m_pRecoTree

private

Definition at line 164 of file PFParticleMonitoring_module.cc.

bool lar_pandora::PFParticleMonitoring::m_printDebug

private

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

Definition at line 259 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_purity

private

Definition at line 221 of file PFParticleMonitoring_module.cc.

bool lar_pandora::PFParticleMonitoring::m_recursiveMatching

private

Definition at line 258 of file PFParticleMonitoring_module.cc.

int lar_pandora::PFParticleMonitoring::m_run

private

Definition at line 166 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_spacepointsMaxX

private

Definition at line 245 of file PFParticleMonitoring_module.cc.

double lar_pandora::PFParticleMonitoring::m_spacepointsMinX

private

Definition at line 244 of file PFParticleMonitoring_module.cc.

std::string lar_pandora::PFParticleMonitoring::m_trackLabel

private

Definition at line 248 of file PFParticleMonitoring_module.cc.

bool lar_pandora::PFParticleMonitoring::m_useDaughterMCParticles

private

Definition at line 254 of file PFParticleMonitoring_module.cc.

bool lar_pandora::PFParticleMonitoring::m_useDaughterPFParticles

private

Definition at line 253 of file PFParticleMonitoring_module.cc.

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

larpandora/larpandora/LArPandoraAnalysis/PFParticleMonitoring_module.cc

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation