#include <AsymmetryFeatureBaseTool.h>

Inheritance diagram for lar_content::AsymmetryFeatureBaseTool:

Public Member Functions
	AsymmetryFeatureBaseTool ()
	Default constructor. More...

void	Run (LArMvaHelper::MvaFeatureVector &featureVector, const VertexSelectionBaseAlgorithm const pAlgorithm, const pandora::Vertex const pVertex, const VertexSelectionBaseAlgorithm::SlidingFitDataListMap &slidingFitDataListMap, const VertexSelectionBaseAlgorithm::ClusterListMap &, const VertexSelectionBaseAlgorithm::KDTreeMap &, const VertexSelectionBaseAlgorithm::ShowerClusterListMap &showerClusterListMap, const float, float &)
	Run the tool. More...

Public Member Functions inherited from lar_content::MvaFeatureTool< Ts >
	MvaFeatureTool ()=default
	Default constructor. More...

virtual void	Run (MvaTypes::MvaFeatureVector &featureVector, Ts...args)=0
	Run the algorithm tool. More...

Protected Member Functions
pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)

virtual float	GetAsymmetryForView (const pandora::CartesianVector &vertexPosition2D, const VertexSelectionBaseAlgorithm::SlidingFitDataList &slidingFitDataList, const VertexSelectionBaseAlgorithm::ShowerClusterList &showerClusterList) const =0
	Get the asymmetry feature for a given view. More...

void	IncrementAsymmetryParameters (const float weight, const pandora::CartesianVector &clusterDirection, pandora::CartesianVector &localWeightedDirectionSum) const
	Increment the asymmetry parameters. More...

virtual float	CalculateAsymmetry (const bool useEnergyMetrics, const pandora::CartesianVector &vertexPosition2D, const pandora::ClusterVector &asymmetryClusters, const pandora::CartesianVector &localWeightedDirectionSum) const
	Calculate the asymmetry feature. More...

Protected Attributes
float	m_maxAsymmetryDistance
	The max distance between cluster (any hit) and vertex to calculate asymmetry score. More...

Additional Inherited Members
Public Types inherited from lar_content::MvaFeatureTool< Ts >
typedef std::vector< MvaFeatureTool< Ts... > * >	FeatureToolVector

Detailed Description

AsymmetryFeatureBaseTool class.

Definition at line 19 of file AsymmetryFeatureBaseTool.h.

Constructor & Destructor Documentation

lar_content::AsymmetryFeatureBaseTool::AsymmetryFeatureBaseTool ( )

Default constructor.

Definition at line 19 of file AsymmetryFeatureBaseTool.cc.

                                                    : m_maxAsymmetryDistance(5.f)
 {
 }

Member Function Documentation

float lar_content::AsymmetryFeatureBaseTool::CalculateAsymmetry	(	const bool	useEnergyMetrics,
		const pandora::CartesianVector &	vertexPosition2D,
		const pandora::ClusterVector &	asymmetryClusters,
		const pandora::CartesianVector &	localWeightedDirectionSum
	)		const

protectedvirtual

Calculate the asymmetry feature.

Parameters

useEnergyMetrics	whether to use energy-based metrics instead of hit-counting-based metrics
vertexPosition2D	the vertex position in this view
asymmetryCluster	the vector of cluster objects to be used in the asymmetry calculation
localWeightedDirectionSum	the local event axis

Returns: the asymmetry feature

Reimplemented in lar_content::EnergyDepositionAsymmetryFeatureTool.

Definition at line 69 of file AsymmetryFeatureBaseTool.cc.

 {
     // Project every hit onto local event axis direction and record side of the projected vtx position on which it falls
     float beforeVtxHitEnergy(0.f), afterVtxHitEnergy(0.f);
     unsigned int beforeVtxHitCount(0), afterVtxHitCount(0);
 
     const CartesianVector localWeightedDirection(localWeightedDirectionSum.GetUnitVector());
     const float evtProjectedVtxPos(vertexPosition2D.GetDotProduct(localWeightedDirection));
 
     for (const Cluster *const pCluster : asymmetryClusters)
     {
         CaloHitList caloHitList;
         pCluster->GetOrderedCaloHitList().FillCaloHitList(caloHitList);
 
         CaloHitVector caloHitVector(caloHitList.begin(), caloHitList.end());
         std::sort(caloHitVector.begin(), caloHitVector.end(), LArClusterHelper::SortHitsByPosition);
 
         for (const CaloHit *const pCaloHit : caloHitVector)
         {
             if (pCaloHit->GetPositionVector().GetDotProduct(localWeightedDirection) < evtProjectedVtxPos)
             {
                 beforeVtxHitEnergy += pCaloHit->GetElectromagneticEnergy();
                 ++beforeVtxHitCount;
             }
             else
             {
                 afterVtxHitEnergy += pCaloHit->GetElectromagneticEnergy();
                 ++afterVtxHitCount;
             }
         }
     }
 
     // Use energy metrics if possible, otherwise fall back on hit counting.
     const float totHitEnergy(beforeVtxHitEnergy + afterVtxHitEnergy);
     const unsigned int totHitCount(beforeVtxHitCount + afterVtxHitCount);
 
     if (useEnergyMetrics && (totHitEnergy > std::numeric_limits<float>::epsilon()))
         return std::fabs((afterVtxHitEnergy - beforeVtxHitEnergy)) / totHitEnergy;
 
     if (0 == totHitCount)
         throw StatusCodeException(STATUS_CODE_FAILURE);
 
     return std::fabs((static_cast<float>(afterVtxHitCount) - static_cast<float>(beforeVtxHitCount))) / static_cast<float>(totHitCount);
 }

virtual float lar_content::AsymmetryFeatureBaseTool::GetAsymmetryForView	(	const pandora::CartesianVector &	vertexPosition2D,
		const VertexSelectionBaseAlgorithm::SlidingFitDataList &	slidingFitDataList,
		const VertexSelectionBaseAlgorithm::ShowerClusterList &	showerClusterList
	)		const

protectedpure virtual

Get the asymmetry feature for a given view.

Parameters

vertexPosition2D	the vertex position projected into this view
slidingFitDataList	the list of sliding fit data objects for this view
showerClusterList	the list of shower cluster objects for this view

Returns: the asymmetry feature

Implemented in lar_content::ShowerAsymmetryFeatureTool, lar_content::GlobalAsymmetryFeatureTool, and lar_content::LocalAsymmetryFeatureTool.

void lar_content::AsymmetryFeatureBaseTool::IncrementAsymmetryParameters	(	const float	weight,
		const pandora::CartesianVector &	clusterDirection,
		pandora::CartesianVector &	localWeightedDirectionSum
	)		const

protected

Increment the asymmetry parameters.

Parameters

weight	the weight to assign to this vector
clusterDirection	the direction of the cluster
localWeightedDirectionSum	the current energy-weighted local cluster direction vector

Definition at line 52 of file AsymmetryFeatureBaseTool.cc.

 {
     // If the new axis direction is at an angle of greater than 90 deg to the current axis direction, flip it 180 degs.
     CartesianVector newDirection(clusterDirection);
 
     if (localWeightedDirectionSum.GetMagnitudeSquared() > std::numeric_limits<float>::epsilon())
     {
         if (localWeightedDirectionSum.GetCosOpeningAngle(clusterDirection) < 0.f)
             newDirection *= -1.f;
     }
 
     localWeightedDirectionSum += newDirection * weight;
 }

StatusCode lar_content::AsymmetryFeatureBaseTool::ReadSettings ( const pandora::TiXmlHandle xmlHandle )

protected

Definition at line 117 of file AsymmetryFeatureBaseTool.cc.

 {
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MaxAsymmetryDistance", m_maxAsymmetryDistance));
 
     return STATUS_CODE_SUCCESS;
 }

void lar_content::AsymmetryFeatureBaseTool::Run	(	LArMvaHelper::MvaFeatureVector &	featureVector,
		const VertexSelectionBaseAlgorithm *const	pAlgorithm,
		const pandora::Vertex *const	pVertex,
		const VertexSelectionBaseAlgorithm::SlidingFitDataListMap &	slidingFitDataListMap,
		const VertexSelectionBaseAlgorithm::ClusterListMap &	,
		const VertexSelectionBaseAlgorithm::KDTreeMap &	,
		const VertexSelectionBaseAlgorithm::ShowerClusterListMap &	showerClusterListMap,
		const float	,
		float &
	)

Run the tool.

Parameters

pAlgorithm	address of the calling algorithm
pVertex	address of the vertex
slidingFitDataListMap	map of the sliding fit data lists
showerClusterListMap	map of the shower cluster lists

Returns: the asymmetry feature

Definition at line 25 of file AsymmetryFeatureBaseTool.cc.

 {
     if (PandoraContentApi::GetSettings(*pAlgorithm)->ShouldDisplayAlgorithmInfo())
         std::cout << "----> Running Algorithm Tool: " << this->GetInstanceName() << ", " << this->GetType() << std::endl;
 
     float asymmetry(0.f);
 
     asymmetry += this->GetAsymmetryForView(LArGeometryHelper::ProjectPosition(this->GetPandora(), pVertex->GetPosition(), TPC_VIEW_U),
         slidingFitDataListMap.at(TPC_VIEW_U),
         showerClusterListMap.empty() ? VertexSelectionBaseAlgorithm::ShowerClusterList() : showerClusterListMap.at(TPC_VIEW_U));
 
     asymmetry += this->GetAsymmetryForView(LArGeometryHelper::ProjectPosition(this->GetPandora(), pVertex->GetPosition(), TPC_VIEW_V),
         slidingFitDataListMap.at(TPC_VIEW_V),
         showerClusterListMap.empty() ? VertexSelectionBaseAlgorithm::ShowerClusterList() : showerClusterListMap.at(TPC_VIEW_V));
 
     asymmetry += this->GetAsymmetryForView(LArGeometryHelper::ProjectPosition(this->GetPandora(), pVertex->GetPosition(), TPC_VIEW_W),
         slidingFitDataListMap.at(TPC_VIEW_W),
         showerClusterListMap.empty() ? VertexSelectionBaseAlgorithm::ShowerClusterList() : showerClusterListMap.at(TPC_VIEW_W));
 
     featureVector.push_back(asymmetry);
 }

Member Data Documentation

float lar_content::AsymmetryFeatureBaseTool::m_maxAsymmetryDistance

protected

The max distance between cluster (any hit) and vertex to calculate asymmetry score.

Definition at line 81 of file AsymmetryFeatureBaseTool.h.

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

larpandoracontent/larpandoracontent/LArVertex/AsymmetryFeatureBaseTool.h
larpandoracontent/larpandoracontent/LArVertex/AsymmetryFeatureBaseTool.cc

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation