LArPfoHelper class. More...

#include "LArPfoHelper.h"

Public Member Functions
template<typename T >
void	SlidingFitTrajectoryImpl (const T const pT, const CartesianVector &vertexPosition, const unsigned int layerWindow, const float layerPitch, LArTrackStateVector &trackStateVector, IntVector const pIndexVector)

Static Public Member Functions
static void	GetCoordinateVector (const pandora::ParticleFlowObject *const pPfo, const pandora::HitType &hitType, pandora::CartesianPointVector &coordinateVector)
	Get a list of coordinates of a particular hit type from an input pfos. More...

static void	GetCaloHits (const pandora::PfoList &pfoList, const pandora::HitType &hitType, pandora::CaloHitList &caloHitList)
	Get a list of calo hits of a particular hit type from a list of pfos. More...

static void	GetCaloHits (const pandora::ParticleFlowObject *const pPfo, const pandora::HitType &hitType, pandora::CaloHitList &caloHitList)
	Get a list of calo hits of a particular hit type from a given pfo. More...

static void	GetIsolatedCaloHits (const pandora::PfoList &pfoList, const pandora::HitType &hitType, pandora::CaloHitList &caloHitList)
	Get a list of isolated calo hits of a particular hit type from a list of pfos. More...

static void	GetIsolatedCaloHits (const pandora::ParticleFlowObject *const pPfo, const pandora::HitType &hitType, pandora::CaloHitList &caloHitList)
	Get a list of isolated calo hits of a particular hit type from a given pfo. More...

static void	GetClusters (const pandora::PfoList &pfoList, const pandora::HitType &hitType, pandora::ClusterList &clusterList)
	Get a list of clusters of a particular hit type from a list of pfos. More...

static void	GetClusters (const pandora::ParticleFlowObject *const pPfo, const pandora::HitType &hitType, pandora::ClusterList &clusterList)
	Get a list of clusters of a particular hit type from a given pfo. More...

static void	GetTwoDClusterList (const pandora::ParticleFlowObject *const pPfo, pandora::ClusterList &clusterList)
	Get the list of 2D clusters from an input pfo. More...

static void	GetThreeDClusterList (const pandora::ParticleFlowObject *const pPfo, pandora::ClusterList &clusterList)
	Get the list of 3D clusters from an input pfo. More...

static void	GetAllConnectedPfos (const pandora::PfoList &inputPfoList, pandora::PfoList &outputPfoList)
	Get a flat list of all pfos, recursively including all daughters and parents associated with those pfos in an input list. More...

static void	GetAllConnectedPfos (const pandora::ParticleFlowObject *const pPfo, pandora::PfoList &outputPfoList)
	Get a flat list of all pfos, recursively including all daughters and parents associated with an input pfo. More...

static void	GetAllDownstreamPfos (const pandora::PfoList &inputPfoList, pandora::PfoList &outputPfoList)
	Get a flat list of all pfos, recursively, of all daughters associated with those pfos in an input list. More...

static void	GetAllDownstreamPfos (const pandora::ParticleFlowObject *const pPfo, pandora::PfoList &outputPfoList)
	Get a flat list of all pfos, recursively, of all daughters and parents associated with an input pfo. More...

static float	GetTwoDLengthSquared (const pandora::ParticleFlowObject *const pPfo)
	Calculate length of Pfo using 2D clusters. More...

static float	GetThreeDLengthSquared (const pandora::ParticleFlowObject *const pPfo)
	Calculate length of Pfo using 3D clusters. More...

static float	GetClosestDistance (const pandora::ParticleFlowObject const pPfo, const pandora::Cluster const pCluster)
	Get closest distance between Pfo and cluster. More...

static float	GetTwoDSeparation (const pandora::ParticleFlowObject const pPfo1, const pandora::ParticleFlowObject const pPfo2)
	Get distance between two Pfos using 2D clusters. More...

static float	GetThreeDSeparation (const pandora::ParticleFlowObject const pPfo1, const pandora::ParticleFlowObject const pPfo2)
	Get distance between two Pfos using 3D clusters. More...

static bool	IsTwoD (const pandora::ParticleFlowObject *const pPfo)
	Does Pfo contain 2D clusters. More...

static bool	IsThreeD (const pandora::ParticleFlowObject *const pPfo)
	Does Pfo contain 3D clusters. More...

static bool	IsTrack (const pandora::ParticleFlowObject *const pPfo)
	Return track flag based on Pfo Particle ID. More...

static bool	IsShower (const pandora::ParticleFlowObject *const pPfo)
	Return shower flag based on Pfo Particle ID. More...

static int	GetPrimaryNeutrino (const pandora::ParticleFlowObject *const pPfo)
	Get primary neutrino or antineutrino. More...

static bool	IsFinalState (const pandora::ParticleFlowObject *const pPfo)
	Whether a pfo is a primary parent particle. More...

static bool	IsNeutrinoFinalState (const pandora::ParticleFlowObject *const pPfo)
	Whether a pfo is a final-state particle from a neutrino (or antineutrino) interaction. More...

static bool	IsNeutrino (const pandora::ParticleFlowObject *const pPfo)
	Whether a pfo is a neutrino or (antineutrino) More...

static bool	IsTestBeam (const pandora::ParticleFlowObject *const pPfo)
	Whether a pfo is a test beam particle. More...

static void	GetRecoNeutrinos (const pandora::PfoList *const pPfoList, pandora::PfoList &recoNeutrinos)
	Get neutrino pfos from an input pfo list. More...

static const pandora::ParticleFlowObject *	GetParentPfo (const pandora::ParticleFlowObject *const pPfo)
	Get the primary parent pfo. More...

static const pandora::ParticleFlowObject *	GetParentNeutrino (const pandora::ParticleFlowObject *const pPfo)
	Get primary neutrino or antineutrino. More...

static const pandora::Vertex *	GetVertex (const pandora::ParticleFlowObject *const pPfo)
	Get the pfo vertex. More...

static void	GetSlidingFitTrajectory (const pandora::CartesianPointVector &pointVector, const pandora::CartesianVector &vertexPosition, const unsigned int layerWindow, const float layerPitch, LArTrackStateVector &trackStateVector, pandora::IntVector *const pIndexVector=nullptr)
	Apply 3D sliding fit to a set of 3D points and return track trajectory. More...

static void	GetSlidingFitTrajectory (const pandora::ParticleFlowObject const pPfo, const pandora::Vertex const pVertex, const unsigned int slidingFitHalfWindow, const float layerPitch, LArTrackStateVector &trackStateVector)
	Apply 3D sliding fit to Pfo and return track trajectory. More...

static LArShowerPCA	GetPrincipalComponents (const pandora::CartesianPointVector &pointVector, const pandora::CartesianVector &vertexPosition)
	Perform PCA analysis on a set of 3D points and return results. More...

static LArShowerPCA	GetPrincipalComponents (const pandora::ParticleFlowObject const pPfo, const pandora::Vertex const pVertex)
	Perform PCA analysis on Pfo and return results. More...

static bool	SortByHitProjection (const LArTrackTrajectoryPoint &lhs, const LArTrackTrajectoryPoint &rhs)
	Sort pfos by number of constituent hits. More...

static bool	SortByNHits (const pandora::ParticleFlowObject const pLhs, const pandora::ParticleFlowObject const pRhs)
	Sort pfos by number of constituent hits. More...

Static Private Member Functions
template<typename T >
static void	SlidingFitTrajectoryImpl (const T const pT, const pandora::CartesianVector &vertexPosition, const unsigned int layerWindow, const float layerPitch, LArTrackStateVector &trackStateVector, pandora::IntVector const pIndexVector=nullptr)
	Implementation of sliding fit trajectory extraction. More...

Detailed Description

LArPfoHelper class.

Definition at line 23 of file LArPfoHelper.h.

Member Function Documentation

static void lar_content::LArPfoHelper::GetAllConnectedPfos	(	const pandora::PfoList &	inputPfoList,
		pandora::PfoList &	outputPfoList
	)

static

Get a flat list of all pfos, recursively including all daughters and parents associated with those pfos in an input list.

Parameters

inputPfoList	the input pfo list
outputPfoList	to receive the output pfo list

Referenced by lar_content::EventValidationAlgorithm::FillValidationInfo(), lar_content::CheatingCosmicRayTaggingTool::FindAmbiguousPfos(), lar_content::BdtBeamParticleIdTool::SliceFeatureParameters::Initialize(), lar_pandora::LArPandoraOutput::ProduceArtOutput(), lar_content::MasterAlgorithm::RunCosmicRayHitRemoval(), and lar_content::BeamParticleIdTool::SelectOutputPfos().

static void lar_content::LArPfoHelper::GetAllConnectedPfos	(	const pandora::ParticleFlowObject *const	pPfo,
		pandora::PfoList &	outputPfoList
	)

static

Get a flat list of all pfos, recursively including all daughters and parents associated with an input pfo.

Parameters

pPfo	the input Pfo
outputPfoList	to receive the output pfo list

static void lar_content::LArPfoHelper::GetAllDownstreamPfos	(	const pandora::PfoList &	inputPfoList,
		pandora::PfoList &	outputPfoList
	)

static

Get a flat list of all pfos, recursively, of all daughters associated with those pfos in an input list.

Parameters

inputPfoList	the input pfo list
outputPfoList	to receive the output pfo list

Referenced by lar_content::LArMCParticleHelper::CollectReconstructable2DHits(), lar_content::CheatingCosmicRayTaggingTool::FindAmbiguousPfos(), lar_content::CosmicRayVertexBuildingAlgorithm::GetCosmicPfos(), lar_content::NeutrinoDaughterVerticesAlgorithm::GetDaughterPfos(), lar_content::CheatingNeutrinoDaughterVerticesAlgorithm::ProcessRecoNeutrinos(), lar_content::CheatingCosmicRayIdentificationAlg::Run(), lar_content::CheatingBeamParticleIdTool::SelectOutputPfos(), and lar_content::MasterAlgorithm::ShiftPfoHierarchy().

static void lar_content::LArPfoHelper::GetAllDownstreamPfos	(	const pandora::ParticleFlowObject *const	pPfo,
		pandora::PfoList &	outputPfoList
	)

static

Get a flat list of all pfos, recursively, of all daughters and parents associated with an input pfo.

Parameters

pPfo	the input Pfo
outputPfoList	to receive the output pfo list

static void lar_content::LArPfoHelper::GetCaloHits	(	const pandora::PfoList &	pfoList,
		const pandora::HitType &	hitType,
		pandora::CaloHitList &	caloHitList
	)

static

Get a list of calo hits of a particular hit type from a list of pfos.

Parameters

pfoList	the input list of Pfos
hitType	the cluster hit type
caloHitList	the output list of calo hits

Referenced by lar_content::CosmicRayTaggingTool::CheckIfInTime(), lar_content::NeutrinoIdTool::Collect2DHits(), lar_content::BdtBeamParticleIdTool::Collect2DHits(), lar_content::LArMCParticleHelper::CollectReconstructable2DHits(), lar_content::BdtBeamParticleIdTool::SliceFeatureParameters::Initialize(), lar_content::CosmicRayTaggingMonitoringTool::PrintPfoTable(), lar_pandora::LArPandoraOutput::ProduceArtOutput(), lar_content::TestBeamParticleCreationAlgorithm::Run(), lar_content::ThreeDOpeningAngleFeatureTool::Run(), lar_content::ThreeDPCAFeatureTool::Run(), lar_content::MasterAlgorithm::RunSlicing(), and lar_content::BeamParticleIdTool::SelectOutputPfos().

static void lar_content::LArPfoHelper::GetCaloHits	(	const pandora::ParticleFlowObject *const	pPfo,
		const pandora::HitType &	hitType,
		pandora::CaloHitList &	caloHitList
	)

static

Get a list of calo hits of a particular hit type from a given pfo.

Parameters

pPfo	the input Pfo
hitType	the cluster hit type
caloHitList	the output list of calo hits

float lar_content::LArPfoHelper::GetClosestDistance	(	const pandora::ParticleFlowObject *const	pPfo,
		const pandora::Cluster *const	pCluster
	)

static

Get closest distance between Pfo and cluster.

Parameters

pPfo	the address of the input Pfo
pCluster	the address of the input cluster

Definition at line 202 of file LArPfoHelper.cc.

References max.

 {
     const HitType hitType(LArClusterHelper::GetClusterHitType(pCluster));
 
     ClusterList clusterList;
     LArPfoHelper::GetClusters(pPfo, hitType, clusterList);
 
     if (clusterList.empty())
         throw StatusCodeException(STATUS_CODE_NOT_FOUND);
 
     float bestDistance(std::numeric_limits<float>::max());
 
     for (const Cluster *const pPfoCluster : clusterList)
     {
         const float thisDistance(LArClusterHelper::GetClosestDistance(pCluster, pPfoCluster));
 
         if (thisDistance < bestDistance)
             bestDistance = thisDistance;
     }
 
     return bestDistance;
 }

static void lar_content::LArPfoHelper::GetClusters	(	const pandora::PfoList &	pfoList,
		const pandora::HitType &	hitType,
		pandora::ClusterList &	clusterList
	)

static

Get a list of clusters of a particular hit type from a list of pfos.

Parameters

pfoList	the input list of Pfos
hitType	the cluster hit type
clusterList	the output list of clusters

Referenced by lar_content::DeltaRayMatchingAlgorithm::AddToDaughterPfo(), lar_content::CosmicRayVertexBuildingAlgorithm::BuildCosmicRayDaughter(), lar_content::CosmicRayVertexBuildingAlgorithm::BuildCosmicRayParent(), lar_content::NeutrinoDaughterVerticesAlgorithm::BuildDaughterShower(), lar_content::NeutrinoDaughterVerticesAlgorithm::BuildDaughterTrack(), lar_content::CosmicRayVertexBuildingAlgorithm::BuildPointingClusterMap(), lar_content::NeutrinoDaughterVerticesAlgorithm::BuildPointingClusterMap(), lar_content::CosmicRayVertexBuildingAlgorithm::GetCosmicPfos(), lar_content::NeutrinoDaughterVerticesAlgorithm::GetDaughterPfos(), lar_content::DeltaRayMatchingAlgorithm::GetDistanceSquaredToPfo(), lar_content::DeltaRayGrowingAlgorithm::GetListOfSeedClusters(), lar_content::DeltaRayIdentificationAlgorithm::GetTwoDSeparation(), lar_content::DeltaRayIdentificationAlgorithm::GetTwoDVertices(), lar_content::SvmPfoCharacterisationAlgorithm::IsClearTrack(), and lar_content::ThreeDChargeFeatureTool::Run().

static void lar_content::LArPfoHelper::GetClusters	(	const pandora::ParticleFlowObject *const	pPfo,
		const pandora::HitType &	hitType,
		pandora::ClusterList &	clusterList
	)

static

Get a list of clusters of a particular hit type from a given pfo.

Parameters

pPfo	the input Pfo
hitType	the cluster hit type
clusterList	the output list of clusters

void lar_content::LArPfoHelper::GetCoordinateVector	(	const pandora::ParticleFlowObject *const	pPfo,
		const pandora::HitType &	hitType,
		pandora::CartesianPointVector &	coordinateVector
	)

static

Get a list of coordinates of a particular hit type from an input pfos.

Parameters

pPfo	the address of the input Pfo
hitType	the cluster hit type
coordinateVector	the output list of coordinates

Definition at line 27 of file LArPfoHelper.cc.

 {
     ClusterList clusterList;
     LArPfoHelper::GetClusters(pPfo, hitType, clusterList);
 
     for (const Cluster *const pCluster : clusterList)
         LArClusterHelper::GetCoordinateVector(pCluster, coordinateVector);
 }

static void lar_content::LArPfoHelper::GetIsolatedCaloHits	(	const pandora::PfoList &	pfoList,
		const pandora::HitType &	hitType,
		pandora::CaloHitList &	caloHitList
	)

static

Get a list of isolated calo hits of a particular hit type from a list of pfos.

Parameters

pfoList	the input list of Pfos
hitType	the cluster hit type
caloHitList	the output list of calo hits

Referenced by lar_content::LArMCParticleHelper::CollectReconstructable2DHits().

static void lar_content::LArPfoHelper::GetIsolatedCaloHits	(	const pandora::ParticleFlowObject *const	pPfo,
		const pandora::HitType &	hitType,
		pandora::CaloHitList &	caloHitList
	)

static

Get a list of isolated calo hits of a particular hit type from a given pfo.

Parameters

pPfo	the input Pfo
hitType	the cluster hit type
caloHitList	the output list of isolated calo hits

const ParticleFlowObject * lar_content::LArPfoHelper::GetParentNeutrino ( const pandora::ParticleFlowObject *const pPfo )

static

Get primary neutrino or antineutrino.

Parameters

pPfo	the address of the Pfo

Returns: address of primary neutrino pfo

Definition at line 418 of file LArPfoHelper.cc.

Referenced by lar_content::EventValidationAlgorithm::ProcessOutput().

 {
     const ParticleFlowObject *const pParentPfo = LArPfoHelper::GetParentPfo(pPfo);
 
     if(!LArPfoHelper::IsNeutrino(pParentPfo))
         throw StatusCodeException(STATUS_CODE_NOT_FOUND);
 
     return pParentPfo;
 }

const ParticleFlowObject * lar_content::LArPfoHelper::GetParentPfo ( const pandora::ParticleFlowObject *const pPfo )

static

Get the primary parent pfo.

Parameters

pPfo	the address of the Pfo

Returns: address of the primary parent pfo

Definition at line 404 of file LArPfoHelper.cc.

 {
     const ParticleFlowObject *pParentPfo = pPfo;
 
     while (pParentPfo->GetParentPfoList().empty() == false)
     {
         pParentPfo = *(pParentPfo->GetParentPfoList().begin());
     }
 
     return pParentPfo;
 }

int lar_content::LArPfoHelper::GetPrimaryNeutrino ( const pandora::ParticleFlowObject *const pPfo )

static

Get primary neutrino or antineutrino.

Parameters

pPfo	the address of the Pfo

Returns: pdg code of neutrino (or zero, otherwise)

Definition at line 330 of file LArPfoHelper.cc.

 {
     try
     {
         const ParticleFlowObject *const pParentPfo = LArPfoHelper::GetParentNeutrino(pPfo);
         return pParentPfo->GetParticleId();
     }
     catch (const StatusCodeException &)
     {
         return 0;
     }
 }

static LArShowerPCA lar_content::LArPfoHelper::GetPrincipalComponents	(	const pandora::CartesianPointVector &	pointVector,
		const pandora::CartesianVector &	vertexPosition
	)

static

Perform PCA analysis on a set of 3D points and return results.

Parameters

pointVector	the input list of 3D positions
vertexPosition	the input vertex position

Referenced by lar_content::PcaShowerParticleBuildingAlgorithm::CreatePfo(), and lar_pandora::LArPandoraShowerCreation::produce().

static LArShowerPCA lar_content::LArPfoHelper::GetPrincipalComponents	(	const pandora::ParticleFlowObject *const	pPfo,
		const pandora::Vertex *const	pVertex
	)

static

Perform PCA analysis on Pfo and return results.

Parameters

pPfo	the address of the input Pfo
pVertex	the address of the input vertex

void lar_content::LArPfoHelper::GetRecoNeutrinos	(	const pandora::PfoList *const	pPfoList,
		pandora::PfoList &	recoNeutrinos
	)

static

Get neutrino pfos from an input pfo list.

Parameters

pPfoList	the input pfo list
recoNeutrinos	to receive the list of neutrino pfos

Definition at line 390 of file LArPfoHelper.cc.

Referenced by lar_content::CheatingNeutrinoDaughterVerticesAlgorithm::Run().

 {
     if (!pPfoList)
         return;
 
     for (const ParticleFlowObject *const pPfo : *pPfoList)
     {
         if (LArPfoHelper::IsNeutrino(pPfo))
             recoNeutrinos.push_back(pPfo);
     }
 }

static void lar_content::LArPfoHelper::GetSlidingFitTrajectory	(	const pandora::CartesianPointVector &	pointVector,
		const pandora::CartesianVector &	vertexPosition,
		const unsigned int	layerWindow,
		const float	layerPitch,
		LArTrackStateVector &	trackStateVector,
		pandora::IntVector *const	pIndexVector = `nullptr`
	)

static

Apply 3D sliding fit to a set of 3D points and return track trajectory.

Parameters

pointVector	the input list of 3D positions
vertexPosition	the input vertex position
slidingFitHalfWindow	size of half window for sliding linear fit
layerPitch	size of pitch for sliding linear fit
trackStateVector	the output track trajectory
pIndexVector	lookup vector of spacepoint indices to store trajectory point sorting

Referenced by lar_content::TrackParticleBuildingAlgorithm::CreatePfo(), and lar_pandora::LArPandoraTrackCreation::produce().

static void lar_content::LArPfoHelper::GetSlidingFitTrajectory	(	const pandora::ParticleFlowObject *const	pPfo,
		const pandora::Vertex *const	pVertex,
		const unsigned int	slidingFitHalfWindow,
		const float	layerPitch,
		LArTrackStateVector &	trackStateVector
	)

static

Apply 3D sliding fit to Pfo and return track trajectory.

Parameters

pPfo	the address of the input Pfo
pVertex	the address of the input vertex
slidingFitHalfWindow	size of half window for sliding linear fit
layerPitch	size of pitch for sliding linear fit
trackStateVector	the output track trajectory

void lar_content::LArPfoHelper::GetThreeDClusterList	(	const pandora::ParticleFlowObject *const	pPfo,
		pandora::ClusterList &	clusterList
	)

static

Get the list of 3D clusters from an input pfo.

Parameters

pPfo	the input Pfo
clusterList	the output list of clusters

Definition at line 110 of file LArPfoHelper.cc.

Referenced by lar_content::ThreeDHitCreationAlgorithm::AddThreeDHitsToPfo(), lar_content::StitchingCosmicRayMergingTool::BuildPointingClusterMaps(), lar_content::CosmicRayTaggingTool::CRCandidate::CRCandidate(), lar_content::NeutrinoIdTool::SliceFeatures::GetSpacePoints(), lar_content::EventSlicingTool::GetThreeDClusters(), lar_content::SlidingConePfoMopUpAlgorithm::GetThreeDClusters(), lar_content::SlidingConeClusterMopUpAlgorithm::GetThreeDClusters(), lar_content::CosmicRayTaggingTool::GetValid3DCluster(), lar_content::NeutrinoHierarchyAlgorithm::PfoInfo::PfoInfo(), lar_content::MasterAlgorithm::Reset(), lar_content::ThreeDOpeningAngleFeatureTool::Run(), and lar_content::ThreeDPCAFeatureTool::Run().

 {
     for (const Cluster *const pCluster : pPfo->GetClusterList())
     {
         if (TPC_3D != LArClusterHelper::GetClusterHitType(pCluster))
             continue;
 
         clusterList.push_back(pCluster);
     }
 }

float lar_content::LArPfoHelper::GetThreeDLengthSquared ( const pandora::ParticleFlowObject *const pPfo )

static

Calculate length of Pfo using 3D clusters.

Parameters

pPfo	the input Pfo

Returns: length variable

Definition at line 182 of file LArPfoHelper.cc.

References f.

Referenced by lar_content::ThreeDLinearFitFeatureTool::Run().

 {
     if (!LArPfoHelper::IsThreeD(pPfo))
         throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);
 
     float lengthSquared(0.f);
 
     for (const Cluster *const pCluster : pPfo->GetClusterList())
     {
         if (TPC_3D != LArClusterHelper::GetClusterHitType(pCluster))
             continue;
 
         lengthSquared += LArClusterHelper::GetLengthSquared(pCluster);
     }
 
     return lengthSquared;
 }

float lar_content::LArPfoHelper::GetThreeDSeparation	(	const pandora::ParticleFlowObject *const	pPfo1,
		const pandora::ParticleFlowObject *const	pPfo2
	)

static

Get distance between two Pfos using 3D clusters.

Parameters

pPfo	the address of the first Pfo
pPfo	the address of the second Pfo

Definition at line 269 of file LArPfoHelper.cc.

 {
     ClusterList clusterList1, clusterList2;
 
     LArPfoHelper::GetClusters(pPfo1, TPC_3D, clusterList1);
     LArPfoHelper::GetClusters(pPfo2, TPC_3D, clusterList2);
 
     if (clusterList1.empty() || clusterList2.empty())
         throw StatusCodeException(STATUS_CODE_NOT_FOUND);
 
     return LArClusterHelper::GetClosestDistance(clusterList1, clusterList2);
 }

void lar_content::LArPfoHelper::GetTwoDClusterList	(	const pandora::ParticleFlowObject *const	pPfo,
		pandora::ClusterList &	clusterList
	)

static

Get the list of 2D clusters from an input pfo.

Parameters

pPfo	the input Pfo
clusterList	the output list of clusters

Definition at line 97 of file LArPfoHelper.cc.

Referenced by lar_content::EventSlicingTool::CopyPfoHitsToSlices(), lar_content::LArMCParticleHelper::GetMainMCParticle(), lar_content::PfoCharacterisationBaseAlgorithm::IsClearTrack3x2D(), lar_content::LArParticleIdPlugins::LArMuonId::IsMatch(), lar_content::MasterAlgorithm::Reset(), lar_content::PfoCharacterisationBaseAlgorithm::Run(), lar_content::ThreeDLinearFitFeatureTool::Run(), and lar_content::ThreeDHitCreationAlgorithm::SeparateTwoDHits().

 {
     for (const Cluster *const pCluster : pPfo->GetClusterList())
     {
         if (TPC_3D == LArClusterHelper::GetClusterHitType(pCluster))
             continue;
 
         clusterList.push_back(pCluster);
     }
 }

float lar_content::LArPfoHelper::GetTwoDLengthSquared ( const pandora::ParticleFlowObject *const pPfo )

static

Calculate length of Pfo using 2D clusters.

Parameters

pPfo	the input Pfo

Returns: length variable

Definition at line 162 of file LArPfoHelper.cc.

References f.

Referenced by lar_content::DeltaRayMatchingAlgorithm::GetLengthFromCache(), and lar_content::DeltaRayIdentificationAlgorithm::IsAssociated().

 {
     if (!LArPfoHelper::IsTwoD(pPfo))
         throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);
 
     float lengthSquared(0.f);
 
     for (const Cluster *const pCluster : pPfo->GetClusterList())
     {
         if (TPC_3D == LArClusterHelper::GetClusterHitType(pCluster))
             continue;
 
         lengthSquared += LArClusterHelper::GetLengthSquared(pCluster);
     }
 
     return lengthSquared;
 }

float lar_content::LArPfoHelper::GetTwoDSeparation	(	const pandora::ParticleFlowObject *const	pPfo1,
		const pandora::ParticleFlowObject *const	pPfo2
	)

static

Get distance between two Pfos using 2D clusters.

Parameters

pPfo	the address of the first Pfo
pPfo	the address of the second Pfo

Definition at line 227 of file LArPfoHelper.cc.

References f.

 {
     ClusterList clusterListU1, clusterListV1, clusterListW1;
     ClusterList clusterListU2, clusterListV2, clusterListW2;
 
     LArPfoHelper::GetClusters(pPfo1, TPC_VIEW_U, clusterListU1);
     LArPfoHelper::GetClusters(pPfo1, TPC_VIEW_V, clusterListV1);
     LArPfoHelper::GetClusters(pPfo1, TPC_VIEW_W, clusterListW1);
 
     LArPfoHelper::GetClusters(pPfo2, TPC_VIEW_U, clusterListU2);
     LArPfoHelper::GetClusters(pPfo2, TPC_VIEW_V, clusterListV2);
     LArPfoHelper::GetClusters(pPfo2, TPC_VIEW_W, clusterListW2);
 
     float numViews(0.f);
     float distanceSquared(0.f);
 
     if (!clusterListU1.empty() && !clusterListU2.empty())
     {
         distanceSquared += LArClusterHelper::GetClosestDistance(clusterListU1, clusterListU2);
         numViews += 1.f;
     }
 
     if (!clusterListV1.empty() && !clusterListV2.empty())
     {
         distanceSquared += LArClusterHelper::GetClosestDistance(clusterListV1, clusterListV2);
         numViews += 1.f;
     }
 
     if (!clusterListW1.empty() && !clusterListW2.empty())
     {
         distanceSquared += LArClusterHelper::GetClosestDistance(clusterListW1, clusterListW2);
         numViews += 1.f;
     }
 
     if (numViews < std::numeric_limits<float>::epsilon())
         throw StatusCodeException(STATUS_CODE_NOT_FOUND);
 
     return std::sqrt(distanceSquared / numViews);
 }

const Vertex * lar_content::LArPfoHelper::GetVertex ( const pandora::ParticleFlowObject *const pPfo )

static

Get the pfo vertex.

Parameters

pPfo	the address of the Pfo

Returns: address of pfo vertex

Definition at line 430 of file LArPfoHelper.cc.

References f.

Referenced by lar_content::CosmicRayTaggingTool::CheckIfInTime(), lar_content::TrackParticleBuildingAlgorithm::CreatePfo(), lar_content::PcaShowerParticleBuildingAlgorithm::CreatePfo(), lar_content::EventValidationAlgorithm::ProcessOutput(), lar_content::TestBeamParticleCreationAlgorithm::Run(), lar_content::ThreeDVertexDistanceFeatureTool::Run(), lar_content::NeutrinoHierarchyAlgorithm::Run(), and lar_content::NeutrinoIdTool::SliceFeatures::SliceFeatures().

 {
     if (pPfo->GetVertexList().empty())
         throw StatusCodeException(STATUS_CODE_NOT_FOUND);
 
     if (pPfo->GetVertexList().size() != 1)
         throw StatusCodeException(STATUS_CODE_FAILURE);
 
     const Vertex *const pVertex = *(pPfo->GetVertexList().begin());
 
     if (VERTEX_3D != pVertex->GetVertexType())
         throw StatusCodeException(STATUS_CODE_FAILURE);
 
     return pVertex;
 }

bool lar_content::LArPfoHelper::IsFinalState ( const pandora::ParticleFlowObject *const pPfo )

static

Whether a pfo is a primary parent particle.

Parameters

pPfo	the address of the Pfo

Returns: boolean

Definition at line 345 of file LArPfoHelper.cc.

Referenced by lar_content::CosmicRayVertexBuildingAlgorithm::BuildCosmicRayParticles(), lar_content::TrackParticleBuildingAlgorithm::CreatePfo(), lar_content::PcaShowerParticleBuildingAlgorithm::CreatePfo(), lar_content::EventValidationAlgorithm::FillValidationInfo(), lar_content::CosmicRayVertexBuildingAlgorithm::GetCosmicPfos(), lar_content::CosmicRayTaggingTool::GetCRCandidates(), lar_content::PfoValidationAlgorithm::Run(), and lar_content::StitchingCosmicRayMergingTool::SelectPrimaryPfos().

 {
     if (pPfo->GetParentPfoList().empty() && !LArPfoHelper::IsNeutrino(pPfo))
         return true;
 
     if (LArPfoHelper::IsNeutrinoFinalState(pPfo))
         return true;
 
     return false;
 }

bool lar_content::LArPfoHelper::IsNeutrino ( const pandora::ParticleFlowObject *const pPfo )

static

Whether a pfo is a neutrino or (antineutrino)

Parameters

pPfo	the address of the Pfo

Returns: boolean

Definition at line 365 of file LArPfoHelper.cc.

Referenced by lar_content::NeutrinoDaughterVerticesAlgorithm::BuildDaughterShower(), lar_content::NeutrinoDaughterVerticesAlgorithm::BuildDaughterTrack(), lar_content::TrackParticleBuildingAlgorithm::CreatePfo(), lar_content::PcaShowerParticleBuildingAlgorithm::CreatePfo(), lar_content::NeutrinoDaughterVerticesAlgorithm::GetDaughterPfos(), lar_content::NeutrinoHierarchyAlgorithm::GetNeutrinoPfo(), lar_content::NeutrinoPropertiesAlgorithm::Run(), lar_content::TestBeamParticleCreationAlgorithm::Run(), and lar_content::CheatingBeamParticleIdTool::SelectOutputPfos().

 {
     const int absoluteParticleId(std::abs(pPfo->GetParticleId()));
 
     if ((NU_E == absoluteParticleId) || (NU_MU == absoluteParticleId) || (NU_TAU == absoluteParticleId))
         return true;
 
     return false;
 }

bool lar_content::LArPfoHelper::IsNeutrinoFinalState ( const pandora::ParticleFlowObject *const pPfo )

static

Whether a pfo is a final-state particle from a neutrino (or antineutrino) interaction.

Parameters

pPfo	the address of the Pfo

Returns: boolean

Definition at line 358 of file LArPfoHelper.cc.

Referenced by lar_content::TrackParticleBuildingAlgorithm::CreatePfo(), lar_content::PcaShowerParticleBuildingAlgorithm::CreatePfo(), and lar_content::EventValidationAlgorithm::ProcessOutput().

 {
     return ((pPfo->GetParentPfoList().size() == 1) && (LArPfoHelper::IsNeutrino(*(pPfo->GetParentPfoList().begin()))));
 }

bool lar_content::LArPfoHelper::IsShower ( const pandora::ParticleFlowObject *const pPfo )

static

Return shower flag based on Pfo Particle ID.

Parameters

pPfo	the address of the Pfo

Definition at line 320 of file LArPfoHelper.cc.

Referenced by lar_content::PcaShowerParticleBuildingAlgorithm::CreatePfo(), lar_content::SlidingConePfoMopUpAlgorithm::GetClusterMergeMap(), lar_content::SlidingConeClusterMopUpAlgorithm::GetThreeDClusters(), lar_content::ShowerHitsBaseTool::Run(), and lar_content::ThreeDHitCreationAlgorithm::Run().

 {
     const int pdg(pPfo->GetParticleId());
 
     // electron, photon
     return ((E_MINUS == std::abs(pdg)) || (PHOTON == std::abs(pdg)));
 }

bool lar_content::LArPfoHelper::IsTestBeam ( const pandora::ParticleFlowObject *const pPfo )

static

Whether a pfo is a test beam particle.

Parameters

pPfo	the address of the Pfo

Returns: boolean

Definition at line 377 of file LArPfoHelper.cc.

Referenced by lar_content::EventValidationAlgorithm::ProcessOutput().

 {
     const PropertiesMap &properties(pPfo->GetPropertiesMap());
     const PropertiesMap::const_iterator iter(properties.find("IsTestBeam"));
 
     if (iter != properties.end())
         return ((iter->second > 0.f) ? true : false);
 
     return false;
 }

bool lar_content::LArPfoHelper::IsThreeD ( const pandora::ParticleFlowObject *const pPfo )

static

Does Pfo contain 3D clusters.

Parameters

pPfo	the address of the Pfo

Definition at line 297 of file LArPfoHelper.cc.

Referenced by lar_content::SvmPfoCharacterisationAlgorithm::IsClearTrack().

 {
     for (const Cluster *const pCluster : pPfo->GetClusterList())
     {
         if (TPC_3D == LArClusterHelper::GetClusterHitType(pCluster))
             return true;
     }
 
     return false;
 }

bool lar_content::LArPfoHelper::IsTrack ( const pandora::ParticleFlowObject *const pPfo )

static

Return track flag based on Pfo Particle ID.

Parameters

pPfo	the address of the Pfo

Definition at line 310 of file LArPfoHelper.cc.

Referenced by lar_content::NeutrinoDaughterVerticesAlgorithm::BuildDaughterParticles(), lar_content::DeltaRayIdentificationAlgorithm::BuildParentDaughterLinks(), lar_content::CosmicRayVertexBuildingAlgorithm::BuildPointingClusterMap(), lar_content::NeutrinoDaughterVerticesAlgorithm::BuildPointingClusterMap(), lar_content::TrackParticleBuildingAlgorithm::CreatePfo(), lar_content::ClusterMopUpBaseAlgorithm::GetPfoClusterLists(), lar_content::SlidingConeClusterMopUpAlgorithm::GetThreeDClusters(), lar_content::SlidingConePfoMopUpAlgorithm::GetThreeDClusters(), lar_content::DeltaRayMatchingAlgorithm::GetTrackPfos(), lar_content::BranchAssociatedPfosTool::Run(), lar_content::TrackHitsBaseTool::Run(), lar_content::ThreeDHitCreationAlgorithm::Run(), and lar_content::StitchingCosmicRayMergingTool::SelectPrimaryPfos().

 {
     const int pdg(pPfo->GetParticleId());
 
     // muon, pion, proton, kaon
     return ((MU_MINUS == std::abs(pdg)) || (PI_PLUS == std::abs(pdg)) || (PROTON == std::abs(pdg)) || (K_PLUS == std::abs(pdg)));
 }

bool lar_content::LArPfoHelper::IsTwoD ( const pandora::ParticleFlowObject *const pPfo )

static

Does Pfo contain 2D clusters.

Parameters

pPfo	the address of the Pfo

Definition at line 284 of file LArPfoHelper.cc.

 {
     for (const Cluster *const pCluster : pPfo->GetClusterList())
     {
         if (TPC_3D != LArClusterHelper::GetClusterHitType(pCluster))
             return true;
     }
 
     return false;
 }

template<typename T >

static void lar_content::LArPfoHelper::SlidingFitTrajectoryImpl	(	const T *const	pT,
		const pandora::CartesianVector &	vertexPosition,
		const unsigned int	layerWindow,
		const float	layerPitch,
		LArTrackStateVector &	trackStateVector,
		pandora::IntVector *const	pIndexVector = `nullptr`
	)

staticprivate

Implementation of sliding fit trajectory extraction.

Parameters

t	the input information
pVertex	the address of the input vertex
slidingFitHalfWindow	size of half window for sliding linear fit
layerPitch	size of pitch for sliding linear fit
trackStateVector	the output track trajectory
pIndexVector	lookup vector of spacepoint indices to store trajectory point sorting

template<typename T >

void lar_content::LArPfoHelper::SlidingFitTrajectoryImpl	(	const T *const	pT,
		const CartesianVector &	vertexPosition,
		const unsigned int	layerWindow,
		const float	layerPitch,
		LArTrackStateVector &	trackStateVector,
		IntVector *const	pIndexVector
	)

Definition at line 558 of file LArPfoHelper.cc.

References f, lar_content::ThreeDSlidingFitResult::GetGlobalFitDirection(), lar_content::ThreeDSlidingFitResult::GetGlobalFitPosition(), lar_content::ThreeDSlidingFitResult::GetGlobalMaxLayerPosition(), lar_content::ThreeDSlidingFitResult::GetGlobalMinLayerPosition(), and lar_content::ThreeDSlidingFitResult::GetLongitudinalDisplacement().

 {
     CartesianPointVector pointVector;
 
     for (const auto &nextPoint : *pT)
         pointVector.push_back(LArObjectHelper::TypeAdaptor::GetPosition(nextPoint));
 
     if (pointVector.empty())
         throw StatusCodeException(STATUS_CODE_NOT_FOUND);
 
     std::sort(pointVector.begin(), pointVector.end(), LArClusterHelper::SortCoordinatesByPosition);
 
     LArTrackTrajectory trackTrajectory;
     IntVector indicesWithoutSpacePoints;
     if (pIndexVector) pIndexVector->clear();
 
     try
     {
         const ThreeDSlidingFitResult slidingFitResult(&pointVector, layerWindow, layerPitch);
         const CartesianVector minPosition(slidingFitResult.GetGlobalMinLayerPosition());
         const CartesianVector maxPosition(slidingFitResult.GetGlobalMaxLayerPosition());
 
         if ((maxPosition - minPosition).GetMagnitudeSquared() < std::numeric_limits<float>::epsilon())
             throw StatusCodeException(STATUS_CODE_NOT_FOUND);
 
         const CartesianVector seedPosition((maxPosition + minPosition) * 0.5f);
         const CartesianVector seedDirection((maxPosition - minPosition).GetUnitVector());
 
         const float scaleFactor((seedDirection.GetDotProduct(seedPosition - vertexPosition) > 0.f) ? +1.f : -1.f);
 
         int index(-1);
         for (const auto &nextPoint : *pT)
         {
             ++index;
 
             try
             {
                 const float rL(slidingFitResult.GetLongitudinalDisplacement(LArObjectHelper::TypeAdaptor::GetPosition(nextPoint)));
 
                 CartesianVector position(0.f, 0.f, 0.f);
                 const StatusCode positionStatusCode(slidingFitResult.GetGlobalFitPosition(rL, position));
 
                 if (positionStatusCode != STATUS_CODE_SUCCESS)
                     throw StatusCodeException(positionStatusCode);
 
                 CartesianVector direction(0.f, 0.f, 0.f);
                 const StatusCode directionStatusCode(slidingFitResult.GetGlobalFitDirection(rL, direction));
 
                 if (directionStatusCode != STATUS_CODE_SUCCESS)
                     throw StatusCodeException(directionStatusCode);
 
                 const float projection(seedDirection.GetDotProduct(position - seedPosition));
 
                 trackTrajectory.push_back(LArTrackTrajectoryPoint(projection * scaleFactor,
                     LArTrackState(position, direction * scaleFactor, LArObjectHelper::TypeAdaptor::GetCaloHit(nextPoint)), index));
             }
             catch (const StatusCodeException &statusCodeException1)
             {
                 indicesWithoutSpacePoints.push_back(index);
 
                 if (STATUS_CODE_FAILURE == statusCodeException1.GetStatusCode())
                     throw statusCodeException1;
             }
         }
     }
     catch (const StatusCodeException &statusCodeException2)
     {
         if (STATUS_CODE_FAILURE == statusCodeException2.GetStatusCode())
             throw statusCodeException2;
     }
 
     // Sort trajectory points by distance along track
     std::sort(trackTrajectory.begin(), trackTrajectory.end(), LArPfoHelper::SortByHitProjection);
 
     for (const LArTrackTrajectoryPoint &larTrackTrajectoryPoint : trackTrajectory)
     {
         trackStateVector.push_back(larTrackTrajectoryPoint.second);
         if (pIndexVector) pIndexVector->push_back(larTrackTrajectoryPoint.GetIndex());
     }
 
     // Store indices of spacepoints with no associated trajectory point at the end of pIndexVector
     if (pIndexVector)
     {
         for (const int index : indicesWithoutSpacePoints)
             pIndexVector->push_back(index);
     }
 }

bool lar_content::LArPfoHelper::SortByHitProjection	(	const LArTrackTrajectoryPoint &	lhs,
		const LArTrackTrajectoryPoint &	rhs
	)

static

Sort pfos by number of constituent hits.

Parameters

pLhs	address of first pfo
pRhs	address of second pfo

Definition at line 500 of file LArPfoHelper.cc.

References f.

 {
     if (lhs.first != rhs.first)
         return (lhs.first < rhs.first);
 
     // ATTN Removed to support use with CartesianVector only (no CaloHit) input
     // if (lhs.second.GetCaloHit() && rhs.second.GetCaloHit())
     //     return (lhs.second.GetCaloHit()->GetInputEnergy() > rhs.second.GetCaloHit()->GetInputEnergy());
 
     const float dx(lhs.second.GetPosition().GetX() - rhs.second.GetPosition().GetX());
     const float dy(lhs.second.GetPosition().GetY() - rhs.second.GetPosition().GetY());
     const float dz(lhs.second.GetPosition().GetZ() - rhs.second.GetPosition().GetZ());
     return (dx + dy + dz > 0.f);
 }

bool lar_content::LArPfoHelper::SortByNHits	(	const pandora::ParticleFlowObject *const	pLhs,
		const pandora::ParticleFlowObject *const	pRhs
	)

static

Sort pfos by number of constituent hits.

Parameters

pLhs	address of first pfo
pRhs	address of second pfo

Definition at line 517 of file LArPfoHelper.cc.

References f.

Referenced by lar_content::DeltaRayIdentificationAlgorithm::BuildParentDaughterLinks(), lar_content::CosmicRayTaggingMonitoringTool::CalculatePfoMetrics(), lar_content::NeutrinoHierarchyAlgorithm::DisplayPfoInfoMap(), lar_content::DeltaRayMatchingAlgorithm::GetAllPfos(), lar_content::LArMCParticleHelper::GetPfoMCParticleHitSharingMaps(), lar_content::DeltaRayGrowingAlgorithm::GetPfos(), lar_content::DeltaRayIdentificationAlgorithm::GetPfos(), lar_content::EventValidationAlgorithm::InterpretMatching(), lar_content::VertexBasedPfoMopUpAlgorithm::PfoAssociation::operator<(), lar_content::StitchingCosmicRayMergingTool::OrderPfoMerges(), lar_content::NeutrinoHierarchyAlgorithm::ProcessPfoInfoMap(), lar_pandora::LArPandoraOutput::ProduceArtOutput(), lar_content::VertexAssociatedPfosTool::Run(), lar_content::ThreeDHitCreationAlgorithm::Run(), lar_content::StitchingCosmicRayMergingTool::SelectPfoMatches(), lar_content::StitchingCosmicRayMergingTool::SelectPfoMerges(), lar_content::StitchingCosmicRayMergingTool::SelectPrimaryPfos(), lar_content::NeutrinoHierarchyAlgorithm::SeparatePfos(), lar_content::NeutrinoPropertiesAlgorithm::SetNeutrinoId(), and lar_content::StitchingCosmicRayMergingTool::StitchPfos().

 {
     unsigned int nTwoDHitsLhs(0), nThreeDHitsLhs(0); float energyLhs(0.f);
     for (ClusterList::const_iterator iter = pLhs->GetClusterList().begin(), iterEnd = pLhs->GetClusterList().end(); iter != iterEnd; ++iter)
     {
         const Cluster *const pClusterLhs = *iter;
 
         if (TPC_3D != LArClusterHelper::GetClusterHitType(pClusterLhs))
             nTwoDHitsLhs += pClusterLhs->GetNCaloHits();
         else
             nThreeDHitsLhs += pClusterLhs->GetNCaloHits();
 
         energyLhs += pClusterLhs->GetHadronicEnergy();
     }
 
     unsigned int nTwoDHitsRhs(0), nThreeDHitsRhs(0); float energyRhs(0.f);
     for (ClusterList::const_iterator iter = pRhs->GetClusterList().begin(), iterEnd = pRhs->GetClusterList().end(); iter != iterEnd; ++iter)
     {
         const Cluster *const pClusterRhs = *iter;
 
         if (TPC_3D != LArClusterHelper::GetClusterHitType(pClusterRhs))
             nTwoDHitsRhs += pClusterRhs->GetNCaloHits();
         else
             nThreeDHitsRhs += pClusterRhs->GetNCaloHits();
 
         energyRhs += pClusterRhs->GetHadronicEnergy();
     }
 
     if (nTwoDHitsLhs != nTwoDHitsRhs)
         return (nTwoDHitsLhs > nTwoDHitsRhs);
 
     if (nThreeDHitsLhs != nThreeDHitsRhs)
         return (nThreeDHitsLhs > nThreeDHitsRhs);
 
     // ATTN Need an efficient (balance with well-motivated) tie-breaker here. Pfo length, for instance, is extremely slow.
     return (energyLhs > energyRhs);
 }

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

larpandoracontent/v03_13_02/source/larpandoracontent/LArHelpers/LArPfoHelper.h
larpandoracontent/v03_13_02/source/larpandoracontent/LArHelpers/LArPfoHelper.cc

Public Member Functions

Static Public Member Functions

Static Private Member Functions

Detailed Description

Member Function Documentation