MCMatches class. More...

#include "LArHierarchyHelper.h"

Public Member Functions
	MCMatches (const MCHierarchy::Node *pMCParticle)
	Constructor. More...

void	AddRecoMatch (const RecoHierarchy::Node *pReco, const int nSharedHits, const pandora::CaloHitList &selectedRecoHits)
	Add a reconstructed node as a match for this MC node. More...

const MCHierarchy::Node *	GetMC () const
	Retrieve the MC node. More...

const RecoHierarchy::NodeVector &	GetRecoMatches () const
	Retrieve the vector of matched reco nodes. More...

const pandora::CaloHitList	GetSelectedRecoHits (const RecoHierarchy::Node *pReco) const
	Retrieve the selected hits for the given reco node. More...

unsigned int	GetSharedHits (const RecoHierarchy::Node *pReco) const
	Retrieve the number of shared hits in the match. More...

float	GetPurity (const RecoHierarchy::Node *pReco, const bool adcWeighted=false) const
	Retrieve the purity of the match. More...

float	GetPurity (const RecoHierarchy::Node *pReco, const pandora::HitType view, const bool adcWeighted=false) const
	Retrieve the purity of the match. More...

float	GetCompleteness (const RecoHierarchy::Node *pReco, const bool adcWeighted=false) const
	Retrieve the completeness of the match. More...

float	GetCompleteness (const RecoHierarchy::Node *pReco, const pandora::HitType view, const bool adcWeighted=false) const
	Retrieve the completeness of the match. More...

size_t	GetNRecoMatches () const
	Get the number of reco nodes matched (both above and below quality cut thresholds) to the MC node. More...

bool	IsQuality (const QualityCuts &qualityCuts) const
	Get whether this match passes quality cuts. More...

Private Member Functions
float	GetPurity (const pandora::CaloHitVector &intersection, const pandora::CaloHitList &recoHits, const bool adcWeighted) const
	Core purity calculation given intersecting hits and reco hits. More...

float	GetCompleteness (const pandora::CaloHitVector &intersection, const pandora::CaloHitList &mcHits, const bool adcWeighted) const
	Core completeness calculation given intersecting hits and MC hits. More...

Private Attributes
const MCHierarchy::Node *	m_pMCParticle
	MC node associated with any matches. More...

RecoHierarchy::NodeVector	m_recoNodes
	Matched reco nodes. More...

pandora::IntVector	m_sharedHits
	Number of shared hits for each match. More...

SelectedRecoHitsMap	m_selectedRecoHitsMap
	Map storing the selected CaloHits for a given reco node. More...

Detailed Description

MCMatches class.

Definition at line 618 of file LArHierarchyHelper.h.

Constructor & Destructor Documentation

lar_content::LArHierarchyHelper::MCMatches::MCMatches ( const MCHierarchy::Node * pMCParticle )

Constructor.

Parameters

pMCParticle The MCParticle being matched

Definition at line 1093 of file LArHierarchyHelper.cc.

                                                                          :
     m_pMCParticle{pMCParticle}
 {
 }

Member Function Documentation

void lar_content::LArHierarchyHelper::MCMatches::AddRecoMatch	(	const RecoHierarchy::Node *	pReco,
		const int	nSharedHits,
		const pandora::CaloHitList &	selectedRecoHits
	)

Add a reconstructed node as a match for this MC node.

Parameters

pReco	The reconstructed node that matches this MC node
nSharedHits	The number of hits shared between reco and MC nodes
selectedRecoHits	The reco node hits that are selected from the MC particle hits

Definition at line 1100 of file LArHierarchyHelper.cc.

References m_recoNodes, m_selectedRecoHitsMap, and m_sharedHits.

Referenced by lar_content::LArHierarchyHelper::MatchInfo::Match().

 {
     m_recoNodes.emplace_back(pReco);
     m_sharedHits.emplace_back(nSharedHits);
     m_selectedRecoHitsMap.insert(std::make_pair(pReco, selectedRecoHits));
 }

float lar_content::LArHierarchyHelper::MCMatches::GetCompleteness	(	const RecoHierarchy::Node *	pReco,
		const bool	adcWeighted = `false`
	)		const

Retrieve the completeness of the match.

Parameters

pReco	The reco node to consider
adcWeighted	Whether or not to weight completeness according to the charge contribution

Returns: The completeness of the match

Definition at line 1171 of file LArHierarchyHelper.cc.

References lar_content::LArHierarchyHelper::MCHierarchy::Node::GetCaloHits(), GetPurity(), GetSelectedRecoHits(), m_pMCParticle, and m_recoNodes.

Referenced by IsQuality().

 {
     auto iter{std::find(m_recoNodes.begin(), m_recoNodes.end(), pReco)};
     if (iter == m_recoNodes.end())
         throw StatusCodeException(STATUS_CODE_NOT_FOUND);
 
     const CaloHitList recoHits = LArHierarchyHelper::MCMatches::GetSelectedRecoHits(pReco);
     const CaloHitList &mcHits{m_pMCParticle->GetCaloHits()};
     CaloHitVector intersection;
     std::set_intersection(mcHits.begin(), mcHits.end(), recoHits.begin(), recoHits.end(), std::back_inserter(intersection));
 
     return this->GetCompleteness(intersection, mcHits, adcWeighted);
 }

float lar_content::LArHierarchyHelper::MCMatches::GetCompleteness	(	const RecoHierarchy::Node *	pReco,
		const pandora::HitType	view,
		const bool	adcWeighted = `false`
	)		const

Retrieve the completeness of the match.

Parameters

pReco	The reco node to consider
view	The view for which purity should be calculated
adcWeighted	Whether or not to weight completeness according to the charge contribution

Returns: The completeness of the match

float lar_content::LArHierarchyHelper::MCMatches::GetCompleteness	(	const pandora::CaloHitVector &	intersection,
		const pandora::CaloHitList &	mcHits,
		const bool	adcWeighted
	)		const

private

Core completeness calculation given intersecting hits and MC hits.

Parameters

intersection	The intersecting reco and MC hits
mcHits	The MC hits
adcWeighted	Whether or not to weight completeness according to the charge contribution

Returns: The completeness of the match

const LArHierarchyHelper::MCHierarchy::Node * lar_content::LArHierarchyHelper::MCMatches::GetMC ( ) const

inline

Retrieve the MC node.

Returns: The MC node

Definition at line 1030 of file LArHierarchyHelper.h.

Referenced by lar_content::LArHierarchyHelper::MatchInfo::Match().

 {
     return m_pMCParticle;
 }

size_t lar_content::LArHierarchyHelper::MCMatches::GetNRecoMatches ( ) const

inline

Get the number of reco nodes matched (both above and below quality cut thresholds) to the MC node.

Returns: The number of reco nodes matched to the MC node

Definition at line 1044 of file LArHierarchyHelper.h.

 {
     return m_recoNodes.size();
 }

float lar_content::LArHierarchyHelper::MCMatches::GetPurity	(	const RecoHierarchy::Node *	pReco,
		const bool	adcWeighted = `false`
	)		const

Retrieve the purity of the match.

Parameters

pReco	The reco node to consider
adcWeighted	Whether or not to weight purity according to the charge contribution

Returns: The purity of the match

Definition at line 1132 of file LArHierarchyHelper.cc.

References lar_content::LArHierarchyHelper::MCHierarchy::Node::GetCaloHits(), GetSelectedRecoHits(), m_pMCParticle, and m_recoNodes.

Referenced by GetCompleteness(), and IsQuality().

 {
     auto iter{std::find(m_recoNodes.begin(), m_recoNodes.end(), pReco)};
     if (iter == m_recoNodes.end())
         throw StatusCodeException(STATUS_CODE_NOT_FOUND);
 
     const CaloHitList recoHits = LArHierarchyHelper::MCMatches::GetSelectedRecoHits(pReco);
     const CaloHitList &mcHits{m_pMCParticle->GetCaloHits()};
     CaloHitVector intersection;
     std::set_intersection(mcHits.begin(), mcHits.end(), recoHits.begin(), recoHits.end(), std::back_inserter(intersection));
 
     return this->GetPurity(intersection, recoHits, adcWeighted);
 }

float lar_content::LArHierarchyHelper::MCMatches::GetPurity	(	const RecoHierarchy::Node *	pReco,
		const pandora::HitType	view,
		const bool	adcWeighted = `false`
	)		const

Retrieve the purity of the match.

Parameters

pReco	The reco node to consider
view	The view for which purity should be calculated
adcWeighted	Whether or not to weight purity according to the charge contribution

Returns: The purity of the match

float lar_content::LArHierarchyHelper::MCMatches::GetPurity	(	const pandora::CaloHitVector &	intersection,
		const pandora::CaloHitList &	recoHits,
		const bool	adcWeighted
	)		const

private

Core purity calculation given intersecting hits and reco hits.

Parameters

intersection	The intersecting reco and MC hits
recoHits	The reco hits
adcWeighted	Whether or not to weight purity according to the charge contribution

Returns: The purity of the match

const LArHierarchyHelper::RecoHierarchy::NodeVector & lar_content::LArHierarchyHelper::MCMatches::GetRecoMatches ( ) const

inline

Retrieve the vector of matched reco nodes.

Returns: The vector of matched reco nodes

Definition at line 1037 of file LArHierarchyHelper.h.

 {
     return m_recoNodes;
 }

const CaloHitList lar_content::LArHierarchyHelper::MCMatches::GetSelectedRecoHits ( const RecoHierarchy::Node * pReco ) const

Retrieve the selected hits for the given reco node.

Parameters

pReco The reco node pointer

Returns: The list of selected calo reco hits

Definition at line 1109 of file LArHierarchyHelper.cc.

References lar_content::LArHierarchyHelper::RecoHierarchy::Node::GetCaloHits(), and m_selectedRecoHitsMap.

Referenced by GetCompleteness(), and GetPurity().

 {
     auto iter{m_selectedRecoHitsMap.find(pReco)};
     if (iter != m_selectedRecoHitsMap.end())
         return iter->second;
     else
         return pReco->GetCaloHits();
 }

unsigned int lar_content::LArHierarchyHelper::MCMatches::GetSharedHits ( const RecoHierarchy::Node * pReco ) const

Retrieve the number of shared hits in the match.

Parameters

pReco The reco node to consider

Returns: The number of shared hits

Definition at line 1120 of file LArHierarchyHelper.cc.

References m_recoNodes, and m_sharedHits.

 {
     auto iter{std::find(m_recoNodes.begin(), m_recoNodes.end(), pReco)};
     if (iter == m_recoNodes.end())
         throw StatusCodeException(STATUS_CODE_NOT_FOUND);
     int index = iter - m_recoNodes.begin();
 
     return static_cast<int>(m_sharedHits[index]);
 }

bool lar_content::LArHierarchyHelper::MCMatches::IsQuality ( const QualityCuts & qualityCuts ) const

Get whether this match passes quality cuts.

Parameters

qualityCuts The quality cuts to pass

Returns: Whether or not this match passes quality cuts

Definition at line 1266 of file LArHierarchyHelper.cc.

References GetCompleteness(), GetPurity(), lar_content::LArHierarchyHelper::QualityCuts::m_minCompleteness, lar_content::LArHierarchyHelper::QualityCuts::m_minPurity, and m_recoNodes.

 {
     if (m_recoNodes.empty())
         return false;
 
     int nAboveThreshold{0};
     if (m_recoNodes.size() != 1)
     {
         for (const LArHierarchyHelper::RecoHierarchy::Node *const pNode : m_recoNodes)
             if (this->GetCompleteness(pNode) > 0.1f)
                 ++nAboveThreshold;
         if (nAboveThreshold != 1)
             return false;
     }
 
     if (this->GetPurity(m_recoNodes.front()) < qualityCuts.m_minPurity)
         return false;
 
     if (this->GetCompleteness(m_recoNodes.front()) < qualityCuts.m_minCompleteness)
         return false;
 
     return true;
 }

Member Data Documentation

const MCHierarchy::Node* lar_content::LArHierarchyHelper::MCMatches::m_pMCParticle

private

MC node associated with any matches.

Definition at line 750 of file LArHierarchyHelper.h.

Referenced by GetCompleteness(), and GetPurity().

RecoHierarchy::NodeVector lar_content::LArHierarchyHelper::MCMatches::m_recoNodes

private

Matched reco nodes.

Definition at line 751 of file LArHierarchyHelper.h.

Referenced by AddRecoMatch(), GetCompleteness(), GetPurity(), GetSharedHits(), and IsQuality().

SelectedRecoHitsMap lar_content::LArHierarchyHelper::MCMatches::m_selectedRecoHitsMap

private

Map storing the selected CaloHits for a given reco node.

Definition at line 753 of file LArHierarchyHelper.h.

Referenced by AddRecoMatch(), and GetSelectedRecoHits().

pandora::IntVector lar_content::LArHierarchyHelper::MCMatches::m_sharedHits

private

Number of shared hits for each match.

Definition at line 752 of file LArHierarchyHelper.h.

Referenced by AddRecoMatch(), and GetSharedHits().

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

larpandoracontent/v04_15_01/source/larpandoracontent/LArHelpers/LArHierarchyHelper.h
larpandoracontent/v04_15_01/source/larpandoracontent/LArHelpers/LArHierarchyHelper.cc

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation