ConnectionRegionFeatureTool to calculate variables related to the connection pathway region. More...

#include "ConnectionPathwayFeatureTool.h"

Inheritance diagram for lar_content::ConnectionRegionFeatureTool:

Public Types
typedef std::vector< MvaFeatureTool< Ts... > * >	FeatureToolVector

typedef std::map< std::string, MvaFeatureTool< Ts... > * >	FeatureToolMap

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

void	Run (LArMvaHelper::MvaFeatureVector &featureVector, const pandora::Algorithm const pAlgorithm, const pandora::ParticleFlowObject const pShowerPfo, const pandora::CartesianVector &nuVertex3D, const ProtoShowerMatch &protoShowerMatch, const pandora::CartesianPointVector &showerStarts3D)

void	Run (LArMvaHelper::MvaFeatureMap &featureMap, pandora::StringVector &featureOrder, const std::string &featureToolName, const pandora::Algorithm const pAlgorithm, const pandora::ParticleFlowObject const pShowerPfo, const pandora::CartesianVector &nuVertex3D, const ProtoShowerMatch &protoShowerMatch, const pandora::CartesianPointVector &showerStarts3D)

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

virtual void	Run (MvaTypes::MvaFeatureMap &featureMap, pandora::StringVector &featureOrder, const std::string &featureToolName, Ts...)

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

float	Get2DKink (const pandora::Algorithm *const pAlgorithm, const ProtoShowerMatch &protoShowerMatch, const pandora::CartesianVector &showerStart3D) const
	Obtain a cautious estimate of the largest 2D deflection of the connection pathway. More...

float	GetLargest2DKinkFromView (const pandora::Algorithm *const pAlgorithm, const TwoDSlidingFitResult &spineFit, const pandora::HitType hitType, const pandora::CartesianVector &showerStart3D) const
	Obtain a cautious estimate of the largest 2D deflection of a connection pathway in a given view. More...

Private Attributes
float	m_defaultFloat
	Default float value. More...

unsigned int	m_spineFitWindow
	The spine fit window. More...

int	m_nLayersHalfWindow
	The half window of each segment. More...

float	m_pathwayLengthLimit
	pathwayLengthLimit max. limit More...

float	m_pathwayScatteringAngle2DLimit
	pathwayScatteringAngle2DLimit max. limit More...

Detailed Description

ConnectionRegionFeatureTool to calculate variables related to the connection pathway region.

Definition at line 72 of file ConnectionPathwayFeatureTool.h.

Member Typedef Documentation

template<typename... Ts>

typedef std::map<std::string, MvaFeatureTool<Ts...> *> lar_content::MvaFeatureTool< Ts >::FeatureToolMap

inherited

Definition at line 37 of file LArMvaHelper.h.

template<typename... Ts>

typedef std::vector<MvaFeatureTool<Ts...> *> lar_content::MvaFeatureTool< Ts >::FeatureToolVector

inherited

Definition at line 36 of file LArMvaHelper.h.

Constructor & Destructor Documentation

lar_content::ConnectionRegionFeatureTool::ConnectionRegionFeatureTool ( )

Default constructor.

Definition at line 135 of file ConnectionPathwayFeatureTool.cc.

References Get2DKink(), m_pathwayLengthLimit, m_pathwayScatteringAngle2DLimit, and Run().

                                                          :
     m_defaultFloat(-10.f),
     m_spineFitWindow(10),
     m_nLayersHalfWindow(5),
     m_pathwayLengthLimit(30.f),
     m_pathwayScatteringAngle2DLimit(10.f)
 {
 }

Member Function Documentation

float lar_content::ConnectionRegionFeatureTool::Get2DKink	(	const pandora::Algorithm *const	pAlgorithm,
		const ProtoShowerMatch &	protoShowerMatch,
		const pandora::CartesianVector &	showerStart3D
	)		const

private

Obtain a cautious estimate of the largest 2D deflection of the connection pathway.

Parameters

pAlgorithm	the algorithm
protoShowerMatch	the ProtoShower match
showerStart3D	the 3D shower start position

Returns: a cautious estimate of the largest 2D deflection of the connection pathway

Definition at line 187 of file ConnectionPathwayFeatureTool.cc.

References GetLargest2DKinkFromView(), lar_content::LArConnectionPathwayHelper::GetMinMiddleMax(), lar_content::ProtoShowerMatch::GetProtoShowerU(), lar_content::ProtoShowerMatch::GetProtoShowerV(), lar_content::ProtoShowerMatch::GetProtoShowerW(), lar_content::ProtoShower::GetSpineHitList(), lar_content::LArGeometryHelper::GetWirePitch(), m_defaultFloat, and m_spineFitWindow.

Referenced by ConnectionRegionFeatureTool().

 {
     try
     {
         CartesianPointVector spinePositionsU, spinePositionsV, spinePositionsW;
 
         for (HitType hitType : {TPC_VIEW_U, TPC_VIEW_V, TPC_VIEW_W})
         {
             const ProtoShower &protoShower(hitType == TPC_VIEW_U
                     ? protoShowerMatch.GetProtoShowerU()
                     : (hitType == TPC_VIEW_V ? protoShowerMatch.GetProtoShowerV() : protoShowerMatch.GetProtoShowerW()));
             CartesianPointVector &spinePositions(hitType == TPC_VIEW_U ? spinePositionsU : (hitType == TPC_VIEW_V ? spinePositionsV : spinePositionsW));
 
             for (const CaloHit *const pCaloHit : protoShower.GetSpineHitList())
                 spinePositions.push_back(pCaloHit->GetPositionVector());
         }
 
         const TwoDSlidingFitResult spineFitU(&spinePositionsU, m_spineFitWindow, LArGeometryHelper::GetWirePitch(pAlgorithm->GetPandora(), TPC_VIEW_U));
         const TwoDSlidingFitResult spineFitV(&spinePositionsV, m_spineFitWindow, LArGeometryHelper::GetWirePitch(pAlgorithm->GetPandora(), TPC_VIEW_V));
         const TwoDSlidingFitResult spineFitW(&spinePositionsW, m_spineFitWindow, LArGeometryHelper::GetWirePitch(pAlgorithm->GetPandora(), TPC_VIEW_W));
 
         const float scatterAngleU(this->GetLargest2DKinkFromView(pAlgorithm, spineFitU, TPC_VIEW_U, showerStart3D));
         const float scatterAngleV(this->GetLargest2DKinkFromView(pAlgorithm, spineFitV, TPC_VIEW_V, showerStart3D));
         const float scatterAngleW(this->GetLargest2DKinkFromView(pAlgorithm, spineFitW, TPC_VIEW_W, showerStart3D));
 
         float minScatterAngle(m_defaultFloat), middleScatterAngle(m_defaultFloat), maxScatterAngle(m_defaultFloat);
         LArConnectionPathwayHelper::GetMinMiddleMax(scatterAngleU, scatterAngleV, scatterAngleW, minScatterAngle, middleScatterAngle, maxScatterAngle);
 
         return middleScatterAngle;
     }
     catch (...)
     {
     }
 
     return m_defaultFloat;
 }

float lar_content::ConnectionRegionFeatureTool::GetLargest2DKinkFromView	(	const pandora::Algorithm *const	pAlgorithm,
		const TwoDSlidingFitResult &	spineFit,
		const pandora::HitType	hitType,
		const pandora::CartesianVector &	showerStart3D
	)		const

private

Obtain a cautious estimate of the largest 2D deflection of a connection pathway in a given view.

Parameters

pAlgorithm	the algorithm
spineFit	the shower spine fit
hitType	the 2D view
showerStart3D	the 3D shower start position

Returns: a cautious estimate of the largest 2D deflection of a connection pathway in a given view

Definition at line 227 of file ConnectionPathwayFeatureTool.cc.

References f, lar_content::TwoDSlidingFitResult::GetGlobalFitPosition(), lar_content::TwoDSlidingFitResult::GetL(), lar_content::TwoDSlidingFitResult::GetLayer(), lar_content::TwoDSlidingFitResult::GetLayerFitResultMap(), lar_content::TwoDSlidingFitResult::GetLocalPosition(), m_defaultFloat, m_nLayersHalfWindow, and lar_content::LArGeometryHelper::ProjectPosition().

Referenced by Get2DKink().

 {
     const LayerFitResultMap &layerFitResultMap(spineFit.GetLayerFitResultMap());
     const int minLayer(layerFitResultMap.begin()->first), maxLayer(layerFitResultMap.rbegin()->first);
     const int nLayersSpanned(1 + maxLayer - minLayer);
 
     if (nLayersSpanned <= 2 * m_nLayersHalfWindow)
         return m_defaultFloat;
 
     const CartesianVector projectedShowerStart(LArGeometryHelper::ProjectPosition(pAlgorithm->GetPandora(), showerStart3D, hitType));
 
     float showerStartL(0.f), showerStartT(0.f);
     spineFit.GetLocalPosition(projectedShowerStart, showerStartL, showerStartT);
 
     float maxCentralLayer(spineFit.GetLayer(showerStartL) - m_nLayersHalfWindow);
     float minCentralLayer(layerFitResultMap.begin()->first + m_nLayersHalfWindow + 1);
 
     float highestOpeningAngle(m_defaultFloat);
     CartesianVector kinkPosition(0.f, 0.f, 0.f);
     CartesianVector highestKinkPosition(0.f, 0.f, 0.f);
 
     for (LayerFitResultMap::const_iterator iter = layerFitResultMap.begin(), iterEnd = layerFitResultMap.end(); iter != iterEnd; ++iter)
     {
         const int layer(iter->first);
 
         if (layer < minCentralLayer)
             continue;
 
         if (layer > maxCentralLayer)
             continue;
 
         bool found(false);
         float openingAngle2D(std::numeric_limits<float>::max());
 
         float thisHighestOpeningAngle(m_defaultFloat);
         CartesianVector thisHighestKinkPosition(0.f, 0.f, 0.f);
 
         for (int i = 0; i < m_nLayersHalfWindow; ++i)
         {
             const int testLayer(layer + i);
             const float rL(spineFit.GetL(testLayer));
             const float rL1(spineFit.GetL(testLayer - m_nLayersHalfWindow));
             const float rL2(spineFit.GetL(testLayer + m_nLayersHalfWindow));
 
             CartesianVector firstPosition(0.f, 0.f, 0.f), centralPosition(0.f, 0.f, 0.f), secondPosition(0.f, 0.f, 0.f);
 
             if ((STATUS_CODE_SUCCESS != spineFit.GetGlobalFitPosition(rL1, firstPosition)) ||
                 (STATUS_CODE_SUCCESS != spineFit.GetGlobalFitPosition(rL, centralPosition)) ||
                 (STATUS_CODE_SUCCESS != spineFit.GetGlobalFitPosition(rL2, secondPosition)))
             {
                 continue;
             }
 
             const CartesianVector firstDirection((centralPosition - firstPosition).GetUnitVector());
             const CartesianVector secondDirection((secondPosition - centralPosition).GetUnitVector());
 
             const float testOpeningAngle2D(secondDirection.GetOpeningAngle(firstDirection) * 180.f / M_PI);
 
             if (testOpeningAngle2D < openingAngle2D)
             {
                 openingAngle2D = testOpeningAngle2D;
                 found = true;
             }
 
             if (testOpeningAngle2D > thisHighestOpeningAngle)
             {
                 thisHighestOpeningAngle = openingAngle2D;
                 thisHighestKinkPosition = centralPosition;
             }
         }
 
         if (found)
         {
             if (openingAngle2D > highestOpeningAngle)
             {
                 highestOpeningAngle = openingAngle2D;
                 highestKinkPosition = thisHighestKinkPosition;
             }
         }
     }
 
     return highestOpeningAngle;
 }

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

private

Definition at line 314 of file ConnectionPathwayFeatureTool.cc.

References m_defaultFloat, m_nLayersHalfWindow, m_pathwayLengthLimit, m_pathwayScatteringAngle2DLimit, and m_spineFitWindow.

 {
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "DefaultFloat", m_defaultFloat));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "SpineFitWindow", m_spineFitWindow));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "NLayersHalfWindow", m_nLayersHalfWindow));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "PathwayLengthLimit", m_pathwayLengthLimit));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
         XmlHelper::ReadValue(xmlHandle, "PathwayScatteringAngle2DLimit", m_pathwayScatteringAngle2DLimit));
 
     return STATUS_CODE_SUCCESS;
 }

template<typename... Ts>

virtual void lar_content::MvaFeatureTool< Ts >::Run	(	MvaTypes::MvaFeatureVector &	featureVector,
		Ts...	args
	)

pure virtualinherited

Run the algorithm tool.

Parameters

featureVector	the vector of features to append
args	arguments to pass to the tool

template<typename... Ts>

virtual void lar_content::MvaFeatureTool< Ts >::Run	(	MvaTypes::MvaFeatureMap &	featureMap,
		pandora::StringVector &	featureOrder,
		const std::string &	featureToolName,
		Ts...
	)

inlinevirtualinherited

Definition at line 51 of file LArMvaHelper.h.

     {
         (void)featureMap;
         (void)featureOrder;
         (void)featureToolName;
         return;
     };

void lar_content::ConnectionRegionFeatureTool::Run	(	LArMvaHelper::MvaFeatureVector &	featureVector,
		const pandora::Algorithm *const	pAlgorithm,
		const pandora::ParticleFlowObject *const	pShowerPfo,
		const pandora::CartesianVector &	nuVertex3D,
		const ProtoShowerMatch &	protoShowerMatch,
		const pandora::CartesianPointVector &	showerStarts3D
	)

Referenced by ConnectionRegionFeatureTool().

void lar_content::ConnectionRegionFeatureTool::Run	(	LArMvaHelper::MvaFeatureMap &	featureMap,
		pandora::StringVector &	featureOrder,
		const std::string &	featureToolName,
		const pandora::Algorithm *const	pAlgorithm,
		const pandora::ParticleFlowObject *const	pShowerPfo,
		const pandora::CartesianVector &	nuVertex3D,
		const ProtoShowerMatch &	protoShowerMatch,
		const pandora::CartesianPointVector &	showerStarts3D
	)

Member Data Documentation

float lar_content::ConnectionRegionFeatureTool::m_defaultFloat

private

Default float value.

Definition at line 116 of file ConnectionPathwayFeatureTool.h.

Referenced by Get2DKink(), GetLargest2DKinkFromView(), and ReadSettings().

int lar_content::ConnectionRegionFeatureTool::m_nLayersHalfWindow

private

The half window of each segment.

Definition at line 118 of file ConnectionPathwayFeatureTool.h.

Referenced by GetLargest2DKinkFromView(), and ReadSettings().

float lar_content::ConnectionRegionFeatureTool::m_pathwayLengthLimit

private

pathwayLengthLimit max. limit

Definition at line 119 of file ConnectionPathwayFeatureTool.h.

Referenced by ConnectionRegionFeatureTool(), and ReadSettings().

float lar_content::ConnectionRegionFeatureTool::m_pathwayScatteringAngle2DLimit

private

pathwayScatteringAngle2DLimit max. limit

Definition at line 120 of file ConnectionPathwayFeatureTool.h.

Referenced by ConnectionRegionFeatureTool(), and ReadSettings().

unsigned int lar_content::ConnectionRegionFeatureTool::m_spineFitWindow

private

The spine fit window.

Definition at line 117 of file ConnectionPathwayFeatureTool.h.

Referenced by Get2DKink(), and ReadSettings().

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

larpandoracontent/v04_08_01/source/larpandoracontent/LArShowerRefinement/ConnectionPathwayFeatureTool.h
larpandoracontent/v04_08_01/source/larpandoracontent/LArShowerRefinement/ConnectionPathwayFeatureTool.cc

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation