LArMuonLeadingHelper class. More...

#include <LArMuonLeadingHelper.h>

Classes
class	ValidationParameters
	ValidationParameters class. More...

Public Types
typedef std::map< const pandora::MCParticle *, pandora::CaloHitList >	LeadingMCParticleToPostBremsstrahlungHitList

Static Public Member Functions
static MCProcess	GetLeadingProcess (const pandora::MCParticle *const pMCParticle)
	Return the MCProcess of the leading particle (tier 1) in the delta ray/michel hierarchy. More...

static bool	IsDeltaRay (const pandora::MCParticle *const pMCParticle)
	Return true if input MCParticle is a child of a cosmic ray and has 'delta ray' process tag. More...

static bool	IsMichel (const pandora::MCParticle *const pMCParticle)
	Return true if input MCParticle is a child of a cosmic ray and has 'decay' process tag. More...

static bool	IsMuonLeading (const pandora::MCParticle *const pMCParticle)
	Return true if input MCParticle is in tier 1 of the cosmic ray hierarchy. More...

static const pandora::MCParticle *	GetLeadingParticle (const pandora::MCParticle *const pMCParticle)
	Return leading particle in the delta ray/michel hierarchy i.e. tier below cosmic ray. More...

static void	SelectReconstructableLeadingParticles (const pandora::MCParticleList pMCParticleList, const pandora::CaloHitList pCaloHitList, const ValidationParameters &parameters, const pandora::CaloHitList &recoMuonHitList, LArMCParticleHelper::MCContributionMap &selectedMCParticlesToHitsMap)
	Select target, reconstructable mc particles in the cosmic ray hierarchy. More...

static void	GetPfoMatchContamination (const pandora::MCParticle *const pLeadingParticle, const pandora::CaloHitList &matchedPfoHitList, pandora::CaloHitList &parentTrackHits, pandora::CaloHitList &otherTrackHits, pandora::CaloHitList &otherShowerHits)
	Separate a leading pfo hit list according to the true owner of the hit e.g. other shower. More...

static void	GetMuonPfoContaminationContribution (const pandora::CaloHitList &cosmicRayPfoHitList, const pandora::CaloHitList &leadingMCHitList, pandora::CaloHitList &leadingHitsInParentCosmicRay)
	Determine the leading MCParticle hits within a cosmic ray pfo hit list. More...

static float	GetClosestDistance (const pandora::Cluster *const pCluster, const pandora::CartesianPointVector &cartesianPointVector)
	Get closest distance between a specified cluster and list of positions. More...

static float	GetClosestDistance (const pandora::CaloHit *const pCaloHit, const pandora::CartesianPointVector &cartesianPointVector)
	Get closest distance between a specified calo hit and list of positions. More...

static pandora::StatusCode	GetClosestPosition (const pandora::CartesianVector &referencePoint, const pandora::CartesianPointVector &cartesianPointVector, const pandora::Cluster *const pCluster, const float maxDistanceToCluster, const float maxDistanceToReferencePoint, pandora::CartesianVector &closestPosition)
	Get the closest position from an input list of projected positions that lies close to both a reference point and an input cluster. More...

static void	GetClosestPositions (const pandora::CartesianPointVector &cartesianPointVector1, const pandora::Cluster *const pCluster2, pandora::CartesianVector &outputPosition1, pandora::CartesianVector &outputPosition2)
	Get the closest positions between a list of positions and a cluster. More...

Static Private Member Functions
static void	GetMCToLeadingMap (const pandora::MCParticleList *const pMCParticleList, LArMCParticleHelper::MCRelationMap &mcToLeadingMap)
	Construct the hierarchy folding map (cosmic rays folded to themselves, delta ray/michel hierarchy folded to leading particle) More...

static void	SelectCaloHits (const pandora::CaloHitList *const pCaloHitList, const LArMCParticleHelper::MCRelationMap &mcToTargetMCMap, pandora::CaloHitList &selectedCaloHitList, const bool selectInputHits, const float minHitSharingFraction, const pandora::CaloHitList &recoMuonHitList, LeadingMCParticleToPostBremsstrahlungHitList &leadingMCParticleToPostBremsstrahlungHitList)
	Select a subset of calo hits representing those that represent "reconstructable" regions of the event. More...

static bool	RejectBremsstrahlungHits (const pandora::CaloHit *const pCaloHit, LeadingMCParticleToPostBremsstrahlungHitList &leadingMCParticleToPostBremsstrahlungHitList)
	Identify and record the hits that are post-bremsstralung radiation in a cosmic ray hierarchy. More...

static void	AddInPostBremsstrahlungHits (const LeadingMCParticleToPostBremsstrahlungHitList &leadingMCParticleToPostBremsstrahlungHitList, const float maxBremsstrahlungSeparation, LArMCParticleHelper::MCContributionMap &leadingMCToTrueHitListMap)
	Identify the reconstructable post-bremsstrahlung radiation hits. More...

static void	AddInPostBremsstrahlungHits (const pandora::MCParticle *const pLeadingParticle, const LeadingMCParticleToPostBremsstrahlungHitList &leadingMCParticleToPostBremsstrahlungHitList, const float maxBremsstrahlungSeparation, LArMCParticleHelper::MCContributionMap &leadingMCToTrueHitListMap, const pandora::HitType tpcView)
	Identify the reconstructable post-bremsstrahlung radiation hits. More...

static void	SelectLeadingMCParticles (const pandora::MCParticleList *pMCParticleList, pandora::MCParticleVector &selectedParticles)
	Select all tier 0 and tier 1 MCParticles in cosmic ray hierarchies from an input list. More...

Detailed Description

LArMuonLeadingHelper class.

Definition at line 20 of file LArMuonLeadingHelper.h.

Member Typedef Documentation

typedef std::map<const pandora::MCParticle *, pandora::CaloHitList> lar_content::LArMuonLeadingHelper::LeadingMCParticleToPostBremsstrahlungHitList

Definition at line 37 of file LArMuonLeadingHelper.h.

Member Function Documentation

void lar_content::LArMuonLeadingHelper::AddInPostBremsstrahlungHits	(	const LeadingMCParticleToPostBremsstrahlungHitList &	leadingMCParticleToPostBremsstrahlungHitList,
		const float	maxBremsstrahlungSeparation,
		LArMCParticleHelper::MCContributionMap &	leadingMCToTrueHitListMap
	)

staticprivate

Identify the reconstructable post-bremsstrahlung radiation hits.

Parameters

leadingMCParticleToPostBremsstrahlungHitList	the mapping of leading MCParticles to post-bremsstrahlung hits
maxBremsstrahlungSeparation	the maximum separation of a reconstructable post-bremsstrahlung hit from the pre-radiation hits
leadingMCToTrueHitListMap	the mapping of cosmic ray and leading MCParticles to their reconstructable hits

Definition at line 265 of file LArMuonLeadingHelper.cc.

 {
     MCParticleVector leadingMCParticleVector;
     for (auto &entry : leadingMCParticleToPostBremsstrahlungHitList)
         leadingMCParticleVector.push_back(entry.first);
 
     for (const MCParticle *const pLeadingMCParticle : leadingMCParticleVector)
     {
         // Do not add in hits for which there are no main particle hits
         if (leadingMCToTrueHitListMap.find(pLeadingMCParticle) == leadingMCToTrueHitListMap.end())
             continue;
 
         LArMuonLeadingHelper::AddInPostBremsstrahlungHits(pLeadingMCParticle, leadingMCParticleToPostBremsstrahlungHitList,
             maxBremsstrahlungSeparation, leadingMCToTrueHitListMap, TPC_VIEW_U);
         LArMuonLeadingHelper::AddInPostBremsstrahlungHits(pLeadingMCParticle, leadingMCParticleToPostBremsstrahlungHitList,
             maxBremsstrahlungSeparation, leadingMCToTrueHitListMap, TPC_VIEW_V);
         LArMuonLeadingHelper::AddInPostBremsstrahlungHits(pLeadingMCParticle, leadingMCParticleToPostBremsstrahlungHitList,
             maxBremsstrahlungSeparation, leadingMCToTrueHitListMap, TPC_VIEW_W);
     }
 }

static void lar_content::LArMuonLeadingHelper::AddInPostBremsstrahlungHits	(	const pandora::MCParticle *const	pLeadingParticle,
		const LeadingMCParticleToPostBremsstrahlungHitList &	leadingMCParticleToPostBremsstrahlungHitList,
		const float	maxBremsstrahlungSeparation,
		LArMCParticleHelper::MCContributionMap &	leadingMCToTrueHitListMap,
		const pandora::HitType	tpcView
	)

staticprivate

Identify the reconstructable post-bremsstrahlung radiation hits.

Parameters

pLeadingParticle	the address of the input leading MCParticle
leadingMCParticleToPostBremsstrahlungHitList	the mapping of leading MCParticles to post-bremsstrahlung hits
maxBremsstrahlungSeparation	the maximum separation of a reconstructable post-bremsstrahlung hit from the pre-radiation hits
leadingMCToTrueHitListMap	the mapping of cosmic ray and leading MCParticles to their reconstructable hits
tpcView	the TPC view of the hits to be collected

static float lar_content::LArMuonLeadingHelper::GetClosestDistance	(	const pandora::Cluster *const	pCluster,
		const pandora::CartesianPointVector &	cartesianPointVector
	)

static

Get closest distance between a specified cluster and list of positions.

Parameters

pCluster	address of the input cluster
cartesianPointVector	the list of input positions

Returns: the closest distance

static float lar_content::LArMuonLeadingHelper::GetClosestDistance	(	const pandora::CaloHit *const	pCaloHit,
		const pandora::CartesianPointVector &	cartesianPointVector
	)

static

Get closest distance between a specified calo hit and list of positions.

Parameters

pCaloHit	the address of the input calo hit
cartesianPointVector	the list of input positions

Returns: the closest distance

StatusCode lar_content::LArMuonLeadingHelper::GetClosestPosition	(	const pandora::CartesianVector &	referencePoint,
		const pandora::CartesianPointVector &	cartesianPointVector,
		const pandora::Cluster *const	pCluster,
		const float	maxDistanceToCluster,
		const float	maxDistanceToReferencePoint,
		pandora::CartesianVector &	closestPosition
	)

static

Get the closest position from an input list of projected positions that lies close to both a reference point and an input cluster.

Parameters

referencePoint	the input reference point
cartesianPointVector	the input list of projected positions
pCluster	the input cluster
maxDistanceToCluster	the maximum distance to the cluster
maxDistanceToReferencePoint	the maximum distance to the reference point
closestPosition	to receive the closest position if found

Returns: whether a closest position could be found

Definition at line 445 of file LArMuonLeadingHelper.cc.

 {
     bool found(false);
     float shortestDistanceSquared(std::numeric_limits<float>::max());
 
     for (const CartesianVector &testPosition : cartesianPointVector)
     {
         if (LArClusterHelper::GetClosestDistance(testPosition, pCluster) > maxDistanceToCluster)
             continue;
 
         const float separationSquared((testPosition - referencePoint).GetMagnitude());
 
         if (separationSquared > maxDistanceToReferencePoint)
             continue;
 
         if (separationSquared < shortestDistanceSquared)
         {
             shortestDistanceSquared = separationSquared;
             closestPosition = testPosition;
             found = true;
         }
     }
 
     return found ? STATUS_CODE_SUCCESS : STATUS_CODE_NOT_FOUND;
 }

void lar_content::LArMuonLeadingHelper::GetClosestPositions	(	const pandora::CartesianPointVector &	cartesianPointVector1,
		const pandora::Cluster *const	pCluster2,
		pandora::CartesianVector &	outputPosition1,
		pandora::CartesianVector &	outputPosition2
	)

static

Get the closest positions between a list of positions and a cluster.

Parameters

cartesianPointVector1	the input list of positions
pCluster2	the address of the input cluster
outputPosition1	the closest position in the list of positions
outputPosition2	the closest position in the cluster

Definition at line 474 of file LArMuonLeadingHelper.cc.

 {
     bool distanceFound(false);
     float minDistanceSquared(std::numeric_limits<float>::max());
 
     CartesianVector closestPosition1(0.f, 0.f, 0.f);
     CartesianVector closestPosition2(0.f, 0.f, 0.f);
 
     CaloHitList caloHitList2;
     pCluster2->GetOrderedCaloHitList().FillCaloHitList(caloHitList2);
 
     for (const CartesianVector &positionVector1 : cartesianPointVector1)
     {
         for (const CaloHit *const pCaloHit : caloHitList2)
         {
             const CartesianVector &positionVector2(pCaloHit->GetPositionVector());
 
             const float distanceSquared((positionVector1 - positionVector2).GetMagnitudeSquared());
 
             if (distanceSquared < minDistanceSquared)
             {
                 minDistanceSquared = distanceSquared;
                 closestPosition1 = positionVector1;
                 closestPosition2 = positionVector2;
                 distanceFound = true;
             }
         }
     }
 
     if (!distanceFound)
         throw StatusCodeException(STATUS_CODE_NOT_FOUND);
 
     outputPosition1 = closestPosition1;
     outputPosition2 = closestPosition2;
 }

const MCParticle * lar_content::LArMuonLeadingHelper::GetLeadingParticle ( const pandora::MCParticle *const pMCParticle )

static

Return leading particle in the delta ray/michel hierarchy i.e. tier below cosmic ray.

Parameters

pMCParticle the input MCParticle

Definition at line 73 of file LArMuonLeadingHelper.cc.

 {
     if (!pMCParticle)
         throw StatusCodeException(STATUS_CODE_FAILURE);
 
     if (LArMCParticleHelper::GetHierarchyTier(pMCParticle) == 0)
         throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);
 
     const MCParticle *pParentMCParticle = pMCParticle;
 
     while (LArMCParticleHelper::GetHierarchyTier(pParentMCParticle) != 1)
     {
         const MCParticleList &parentList(pParentMCParticle->GetParentList());
 
         if (parentList.size() != 1)
             throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);
 
         pParentMCParticle = parentList.front();
     }
 
     return pParentMCParticle;
 }

MCProcess lar_content::LArMuonLeadingHelper::GetLeadingProcess ( const pandora::MCParticle *const pMCParticle )

static

Return the MCProcess of the leading particle (tier 1) in the delta ray/michel hierarchy.

Parameters

pMCParticle the address of the input MCParticle

Returns: the MCProcess of the leading particle (tier 1) in the delta ray/michel hierarchy

Definition at line 34 of file LArMuonLeadingHelper.cc.

 {
     if (LArMCParticleHelper::GetHierarchyTier(pMCParticle) == 0)
         return MC_PROC_UNKNOWN;
 
     const MCParticle *const pLeadingParticle(LArMuonLeadingHelper::GetLeadingParticle(pMCParticle));
     const LArMCParticle *const pLArMCParticle(dynamic_cast<const LArMCParticle *>(pLeadingParticle));
 
     if (!pLArMCParticle)
         throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);
 
     return pLArMCParticle->GetProcess();
 }

void lar_content::LArMuonLeadingHelper::GetMCToLeadingMap	(	const pandora::MCParticleList *const	pMCParticleList,
		LArMCParticleHelper::MCRelationMap &	mcToLeadingMap
	)

staticprivate

Construct the hierarchy folding map (cosmic rays folded to themselves, delta ray/michel hierarchy folded to leading particle)

Parameters

pMCParticleList	the address of the list of MCParticles
mcToLeadingMap	the hierarchy folding map

Definition at line 98 of file LArMuonLeadingHelper.cc.

 {
     for (const MCParticle *const pMCParticle : *pMCParticleList)
     {
         const MCParticle *const pParentMCParticle(LArMCParticleHelper::GetParentMCParticle(pMCParticle));
 
         if (!LArMCParticleHelper::IsCosmicRay(pParentMCParticle))
             continue;
 
         // For the CRs: fold hits to themselves, for the DRs: fold hits to the leading MCParticle
         if (pMCParticle == pParentMCParticle)
         {
             mcToLeadingMap[pMCParticle] = pMCParticle;
         }
         else
         {
             const MCParticle *const pLeadingMCParticle(LArMuonLeadingHelper::GetLeadingParticle(pMCParticle));
             mcToLeadingMap[pMCParticle] = pLeadingMCParticle;
         }
     }
 }

void lar_content::LArMuonLeadingHelper::GetMuonPfoContaminationContribution	(	const pandora::CaloHitList &	cosmicRayPfoHitList,
		const pandora::CaloHitList &	leadingMCHitList,
		pandora::CaloHitList &	leadingHitsInParentCosmicRay
	)

static

Determine the leading MCParticle hits within a cosmic ray pfo hit list.

Parameters

cosmicRayPfoHitList	the cosmic ray pfo hit list
leadingMCHitList	the list of hits that belong to the leading MCParticle
leadingHitsInParentCosmicRay	the output list of 'stolen' leading MCParticle hits

Definition at line 395 of file LArMuonLeadingHelper.cc.

 {
     for (const CaloHit *const pCaloHit : cosmicRayPfoHitList)
     {
         if (std::find(leadingMCHitList.begin(), leadingMCHitList.end(), pCaloHit) != leadingMCHitList.end())
             leadingHitsInParentCosmicRay.push_back(pCaloHit);
     }
 }

void lar_content::LArMuonLeadingHelper::GetPfoMatchContamination	(	const pandora::MCParticle *const	pLeadingParticle,
		const pandora::CaloHitList &	matchedPfoHitList,
		pandora::CaloHitList &	parentTrackHits,
		pandora::CaloHitList &	otherTrackHits,
		pandora::CaloHitList &	otherShowerHits
	)

static

Separate a leading pfo hit list according to the true owner of the hit e.g. other shower.

Parameters

pLeadingParticle	the address of the input leading MCParticle
matchedPfoHitList	the input leading pfo hit list
parentTrackHits	the output list of hits that belong to the parent cosmic ray
otherTrackHits	the output list of hits that belong to a cosmic ray that is not the parent
otherShowerHits	the output list of hits that belong to a different shower hierarchy

Definition at line 370 of file LArMuonLeadingHelper.cc.

 {
     const MCParticle *const pParentCosmicRay(LArMCParticleHelper::GetParentMCParticle(pLeadingParticle));
 
     for (const CaloHit *const pCaloHit : matchedPfoHitList)
     {
         const MCParticle *const pHitParticle(MCParticleHelper::GetMainMCParticle(pCaloHit));
 
         if (LArMCParticleHelper::IsCosmicRay(pHitParticle))
         {
             (pHitParticle == pParentCosmicRay) ? parentTrackHits.push_back(pCaloHit) : otherTrackHits.push_back(pCaloHit);
         }
         else
         {
             const MCParticle *const pHitLeadingParticle(LArMuonLeadingHelper::GetLeadingParticle(pHitParticle));
 
             if (pHitLeadingParticle != pLeadingParticle)
                 otherShowerHits.push_back(pCaloHit);
         }
     }
 }

bool lar_content::LArMuonLeadingHelper::IsDeltaRay ( const pandora::MCParticle *const pMCParticle )

static

Return true if input MCParticle is a child of a cosmic ray and has 'delta ray' process tag.

Parameters

pMCParticle the input MCParticle

Definition at line 50 of file LArMuonLeadingHelper.cc.

 {
     return (LArMuonLeadingHelper::GetLeadingProcess(pMCParticle) == MC_PROC_MU_IONI);
 }

bool lar_content::LArMuonLeadingHelper::IsMichel ( const pandora::MCParticle *const pMCParticle )

static

Return true if input MCParticle is a child of a cosmic ray and has 'decay' process tag.

Parameters

pMCParticle the input MCParticle

Definition at line 57 of file LArMuonLeadingHelper.cc.

 {
     return (LArMuonLeadingHelper::GetLeadingProcess(pMCParticle) == MC_PROC_DECAY);
 }

bool lar_content::LArMuonLeadingHelper::IsMuonLeading ( const pandora::MCParticle *const pMCParticle )

static

Return true if input MCParticle is in tier 1 of the cosmic ray hierarchy.

Parameters

pMCParticle the input MCParticle

Definition at line 64 of file LArMuonLeadingHelper.cc.

 {
     const MCParticle *const pParentMCParticle(LArMCParticleHelper::GetParentMCParticle(pMCParticle));
 
     return ((LArMCParticleHelper::GetHierarchyTier(pMCParticle) == 1) && (LArMCParticleHelper::IsCosmicRay(pParentMCParticle)));
 }

bool lar_content::LArMuonLeadingHelper::RejectBremsstrahlungHits	(	const pandora::CaloHit *const	pCaloHit,
		LeadingMCParticleToPostBremsstrahlungHitList &	leadingMCParticleToPostBremsstrahlungHitList
	)

staticprivate

Identify and record the hits that are post-bremsstralung radiation in a cosmic ray hierarchy.

Parameters

pCaloHit	the address of the input calo hit
leadingMCParticleToPostBremsstrahlungHitList	the mapping of leading MCParticles to post-bremsstrahlung hits

Returns: whether the hit lies post-bremsstrahlung radiation in a cosmic ray hierarchy

Definition at line 231 of file LArMuonLeadingHelper.cc.

 {
     const MCParticle *const pHitMCParticle(MCParticleHelper::GetMainMCParticle(pCaloHit));
 
     MCParticleList ancestorMCParticleList;
     LArMCParticleHelper::GetAllAncestorMCParticles(pHitMCParticle, ancestorMCParticleList);
 
     bool isPostBremsstrahlung(false);
     const MCParticle *leadingMCParticle(nullptr);
 
     for (const MCParticle *const pAncestorMCParticle : ancestorMCParticleList)
     {
         if (LArMCParticleHelper::GetHierarchyTier(pAncestorMCParticle) == 1)
         {
             if (LArMuonLeadingHelper::IsMuonLeading(pAncestorMCParticle))
                 leadingMCParticle = pAncestorMCParticle;
         }
 
         if (pAncestorMCParticle->GetParticleId() == PHOTON)
             isPostBremsstrahlung = true;
     }
 
     if (isPostBremsstrahlung && leadingMCParticle)
     {
         leadingMCParticleToPostBremsstrahlungHitList[leadingMCParticle].push_back(pCaloHit);
         return true;
     }
 
     return false;
 }

void lar_content::LArMuonLeadingHelper::SelectCaloHits	(	const pandora::CaloHitList *const	pCaloHitList,
		const LArMCParticleHelper::MCRelationMap &	mcToTargetMCMap,
		pandora::CaloHitList &	selectedCaloHitList,
		const bool	selectInputHits,
		const float	minHitSharingFraction,
		const pandora::CaloHitList &	recoMuonHitList,
		LeadingMCParticleToPostBremsstrahlungHitList &	leadingMCParticleToPostBremsstrahlungHitList
	)

staticprivate

Select a subset of calo hits representing those that represent "reconstructable" regions of the event.

Parameters

pCaloHitList	the address of the input calo hit list
mcToTargetMCMap	the MCParticle to leading MCParticle map
selectedCaloHitList	to receive the populated selected calo hit list
selectInputHits	whether to select input hits
minHitSharingFraction	the minimum required charge share of the hit
recoMuonHitList	the list of reconstructed cosmic ray hits
leadingMCParticleToPostBremsstrahlungHitList	the mapping of leading MCParticles to post-bremsstrahlung hits

Definition at line 154 of file LArMuonLeadingHelper.cc.

 {
     if (!selectInputHits)
     {
         selectedCaloHitList.insert(selectedCaloHitList.end(), pCaloHitList->begin(), pCaloHitList->end());
         return;
     }
 
     for (const CaloHit *const pCaloHit : *pCaloHitList)
     {
         try
         {
             const MCParticle *const pHitParticle(MCParticleHelper::GetMainMCParticle(pCaloHit));
 
             if (mcToTargetMCMap.find(pHitParticle) == mcToTargetMCMap.end())
                 continue;
 
             // Remove delta ray hits that have been 'stolen' by the muon
             if (!LArMCParticleHelper::IsCosmicRay(mcToTargetMCMap.at(pHitParticle)))
             {
                 if (std::find(recoMuonHitList.begin(), recoMuonHitList.end(), pCaloHit) != recoMuonHitList.end())
                     continue;
             }
 
             MCParticleVector mcParticleContributionVector;
             for (const auto &mapEntry : pCaloHit->GetMCParticleWeightMap())
                 mcParticleContributionVector.push_back(mapEntry.first);
 
             std::sort(mcParticleContributionVector.begin(), mcParticleContributionVector.end(), PointerLessThan<MCParticle>());
 
             MCParticleWeightMap targetWeightMap;
             for (const MCParticle *const pMCParticle : mcParticleContributionVector)
             {
                 const float weight(pCaloHit->GetMCParticleWeightMap().at(pMCParticle));
                 LArMCParticleHelper::MCRelationMap::const_iterator mcIter = mcToTargetMCMap.find(pMCParticle);
 
                 if (mcToTargetMCMap.end() != mcIter)
                     targetWeightMap[mcIter->second] += weight;
             }
 
             MCParticleVector mcTargetContributionVector;
             for (const auto &mapEntry : targetWeightMap)
                 mcTargetContributionVector.push_back(mapEntry.first);
             std::sort(mcTargetContributionVector.begin(), mcTargetContributionVector.end(), PointerLessThan<MCParticle>());
 
             float bestTargetWeight(0.f), targetWeightSum(0.f);
 
             for (const MCParticle *const pTargetMCParticle : mcTargetContributionVector)
             {
                 const float targetWeight(targetWeightMap.at(pTargetMCParticle));
                 targetWeightSum += targetWeight;
 
                 if (targetWeight > bestTargetWeight)
                 {
                     bestTargetWeight = targetWeight;
                 }
             }
 
             if ((targetWeightSum < std::numeric_limits<float>::epsilon()) || ((bestTargetWeight / targetWeightSum) < minHitSharingFraction))
                 continue;
 
             // Remove and record post bremsstrahlung hits
             if (LArMuonLeadingHelper::RejectBremsstrahlungHits(pCaloHit, leadingMCParticleToPostBremsstrahlungHitList))
                 continue;
 
             selectedCaloHitList.push_back(pCaloHit);
         }
         catch (const StatusCodeException &)
         {
         }
     }
 }

void lar_content::LArMuonLeadingHelper::SelectLeadingMCParticles	(	const pandora::MCParticleList *	pMCParticleList,
		pandora::MCParticleVector &	selectedParticles
	)

staticprivate

Select all tier 0 and tier 1 MCParticles in cosmic ray hierarchies from an input list.

Parameters

pMCParticle	the address of the input MCParticle list
selectedParticles	the output vector of selected MCParticles

Definition at line 345 of file LArMuonLeadingHelper.cc.

 {
     for (const MCParticle *const pMCParticle : *pMCParticleList)
     {
         const MCParticle *const pParentMCParticle(LArMCParticleHelper::GetParentMCParticle(pMCParticle));
 
         if (!LArMCParticleHelper::IsCosmicRay(pParentMCParticle))
             continue;
 
         if (pMCParticle == pParentMCParticle)
         {
             selectedParticles.push_back(pMCParticle);
         }
         else
         {
             if (LArMuonLeadingHelper::IsMuonLeading(pMCParticle))
                 selectedParticles.push_back(pMCParticle);
         }
     }
 
     std::sort(selectedParticles.begin(), selectedParticles.end(), LArMCParticleHelper::SortByMomentum);
 }

void lar_content::LArMuonLeadingHelper::SelectReconstructableLeadingParticles	(	const pandora::MCParticleList *	pMCParticleList,
		const pandora::CaloHitList *	pCaloHitList,
		const ValidationParameters &	parameters,
		const pandora::CaloHitList &	recoMuonHitList,
		LArMCParticleHelper::MCContributionMap &	selectedMCParticlesToHitsMap
	)

static

Select target, reconstructable mc particles in the cosmic ray hierarchy.

Parameters

pMCParticleList	the address of the list of MCParticles
pCaloHitList	the address of the list of CaloHits
parameters	validation parameters to decide when an MCParticle is considered reconstructable
recoMuonHitList	the list of reconstructed cosmic ray hits
selectedMCParticlesToHitsMap	the output mapping from selected MCParticles to their hits

Definition at line 122 of file LArMuonLeadingHelper.cc.

 {
     // Obtain hierarchy folding map:
     LArMCParticleHelper::MCRelationMap mcToLeadingMCMap;
     LArMuonLeadingHelper::GetMCToLeadingMap(pMCParticleList, mcToLeadingMCMap);
 
     // Select reconstructable hits, e.g. remove delta ray hits 'stolen' by the cosmic rays
     CaloHitList selectedCaloHitList;
     LeadingMCParticleToPostBremsstrahlungHitList leadingMCParticleToPostBremsstrahlungHitList;
     LArMuonLeadingHelper::SelectCaloHits(pCaloHitList, mcToLeadingMCMap, selectedCaloHitList, parameters.m_selectInputHits,
         parameters.m_minHitSharingFraction, recoMuonHitList, leadingMCParticleToPostBremsstrahlungHitList);
 
     // Obtain maps: [hit -> leading MCParticle], [leading MCParticle -> list of hits]
     LArMCParticleHelper::CaloHitToMCMap trueHitToLeadingMCMap;
     LArMCParticleHelper::MCContributionMap leadingMCToTrueHitListMap;
     LArMCParticleHelper::GetMCParticleToCaloHitMatches(&selectedCaloHitList, mcToLeadingMCMap, trueHitToLeadingMCMap, leadingMCToTrueHitListMap);
 
     // Add in close post bremsstrahlung hits
     LArMuonLeadingHelper::AddInPostBremsstrahlungHits(
         leadingMCParticleToPostBremsstrahlungHitList, parameters.m_maxBremsstrahlungSeparation, leadingMCToTrueHitListMap);
 
     // Obtain vector: all mc particles
     MCParticleVector leadingMCVector;
     LArMuonLeadingHelper::SelectLeadingMCParticles(pMCParticleList, leadingMCVector);
 
     // Ensure the MCParticles have enough "good" hits to be reconstructed
     LArMCParticleHelper::SelectParticlesByHitCount(leadingMCVector, leadingMCToTrueHitListMap, mcToLeadingMCMap, parameters, selectedMCParticlesToHitsMap);
 }

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

larpandoracontent/larpandoracontent/LArHelpers/LArMuonLeadingHelper.h
larpandoracontent/larpandoracontent/LArHelpers/LArMuonLeadingHelper.cc

Classes

Public Types

Static Public Member Functions

Static Private Member Functions

Detailed Description

Member Typedef Documentation

Member Function Documentation