ShowerRegionFeatureTool to calculate variables related to the shower region. More...

#include "ConnectionPathwayFeatureTool.h"

Inheritance diagram for lar_content::ShowerRegionFeatureTool:

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

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

Public Member Functions
	ShowerRegionFeatureTool ()
	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)

void	GetViewShowerRegionVariables (const pandora::Algorithm const pAlgorithm, const pandora::ParticleFlowObject const pShowerPfo, const pandora::CartesianVector &nuVertex3D, const ProtoShowerMatch &protoShowerMatch, const pandora::HitType hitType, const pandora::CartesianVector &showerStart3D, float &nHits, float &foundHitRatio, float &scatterAngle, float &openingAngle, float &nuVertexEnergyAsymmetry, float &nuVertexEnergyWeightedMeanRadialDistance, float &showerStartEnergyAsymmetry, float &showerStartMoliereRadius)
	Calculate the shower region variables for the input view. More...

void	BuildViewShower (const pandora::ParticleFlowObject *const pShowerPfo, const TwoDSlidingFitResult &spineFit, const pandora::HitType hitType, const pandora::CartesianVector &showerStart2D, const pandora::CartesianVector &nuVertex2D, pandora::CaloHitList &postShowerHitList, pandora::CartesianPointVector &postShowerPositions)
	Collect the shower region hits in a given view. More...

void	GetShowerHitVariables (const pandora::CaloHitList &spineHitList, const pandora::CaloHitList &postShowerHitList, const pandora::ParticleFlowObject *const pShowerPfo, const pandora::HitType hitType, float &nHits, float &foundHitRatio)
	Evaluate the variables associated with the shower region hit multiplicity. More...

void	CalculateViewScatterAngle (const pandora::CartesianVector &nuVertex2D, const TwoDSlidingFitResult &spineFitResult, const pandora::CartesianVector &showerStart2D, const TwoDSlidingFitResult &showerFitResult, float &scatterAngle)
	Calculate the connection pathway-shower region scatter angle. More...

void	CalculateViewOpeningAngle (const TwoDSlidingFitResult &showerFitResult, const pandora::CaloHitList &postShowerHitList, const pandora::CartesianVector &showerStart2D, float &openingAngle)
	Calculate the opening angle of the shower region. More...

void	CalculateViewNuVertexConsistencyVariables (const TwoDSlidingFitResult &spineFitResult, const pandora::CaloHitList &postShowerHitList, const bool isDownstream, const pandora::CartesianVector &nuVertex2D, float &nuVertexEnergyAsymmetry, float &nuVertexEnergyWeightedMeanRadialDistance)
	Evaluate the neutrino vertex consistency variables. More...

void	CalculateViewShowerStartConsistencyVariables (const TwoDSlidingFitResult &showerFitResult, const pandora::CaloHitList &postShowerHitList, const bool isShowerDownstream, float &showerStartEnergyAsymmetry, float &showerStartMoliereRadius)
	Evaluate the shower start consistency variables. More...

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

float	m_defaultRatio
	Default float value for ratios. More...

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

float	m_showerRadius
	The max. separation distance between a shower region hit and the shower core. More...

unsigned int	m_showerFitWindow
	The shower fit window. More...

float	m_edgeStep
	The binning of the shower boundaries. More...

float	m_moliereFraction
	The energy fraction which corresponds to minShowerStartMoliereRadius. More...

float	m_maxNHitsLimit
	maxNHits max. limit More...

float	m_maxFoundHitRatioLimit
	maxFoundHitRatio max. limit More...

float	m_maxScatterAngleLimit
	maxScatterAngle max. limit More...

float	m_maxOpeningAngleLimit
	maxOpeningAngle max. limit More...

float	m_maxNuVertexEnergyWeightedMeanRadialDistanceLimit
	maxNuVertexEnergyWeightedMeanRadialDistance max. limit More...

float	m_minShowerStartMoliereRadiusLimit
	minShowerStartMoliereRadius max. limit More...

Detailed Description

ShowerRegionFeatureTool to calculate variables related to the shower region.

Definition at line 128 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::ShowerRegionFeatureTool::ShowerRegionFeatureTool ( )

Default constructor.

Definition at line 335 of file ConnectionPathwayFeatureTool.cc.

References lar_content::LArConnectionPathwayHelper::GetMinMiddleMax(), GetViewShowerRegionVariables(), m_defaultFloat, m_defaultRatio, m_maxFoundHitRatioLimit, m_maxNHitsLimit, m_maxNuVertexEnergyWeightedMeanRadialDistanceLimit, m_maxOpeningAngleLimit, m_maxScatterAngleLimit, m_minShowerStartMoliereRadiusLimit, and Run().

                                                  :
     m_defaultFloat(-10.f),
     m_defaultRatio(-0.5f),
     m_spineFitWindow(20),
     m_showerRadius(14.f),
     m_showerFitWindow(1000),
     m_edgeStep(2.f),
     m_moliereFraction(0.9f),
     m_maxNHitsLimit(2000.f),
     m_maxFoundHitRatioLimit(1.5f),
     m_maxScatterAngleLimit(40.f),
     m_maxOpeningAngleLimit(20.f),
     m_maxNuVertexEnergyWeightedMeanRadialDistanceLimit(20.f),
     m_minShowerStartMoliereRadiusLimit(10.f)
 {
 }

Member Function Documentation

void lar_content::ShowerRegionFeatureTool::BuildViewShower	(	const pandora::ParticleFlowObject *const	pShowerPfo,
		const TwoDSlidingFitResult &	spineFit,
		const pandora::HitType	hitType,
		const pandora::CartesianVector &	showerStart2D,
		const pandora::CartesianVector &	nuVertex2D,
		pandora::CaloHitList &	postShowerHitList,
		pandora::CartesianPointVector &	postShowerPositions
	)

private

Collect the shower region hits in a given view.

Parameters

pShowerPfo	the shower pfo
spineFit	the shower spine fit
hitType	the 2D view
showerStart2D	the 2D shower start position
nuVertex2D	the 2D neutrino vertex position
postShowerHitList	the collected shower region hit list
postShowerPositions	the collected shower region hit position vector

Definition at line 579 of file ConnectionPathwayFeatureTool.cc.

References f, lar_content::LArPfoHelper::GetCaloHits(), lar_content::TwoDSlidingFitResult::GetGlobalDirection(), lar_content::TwoDSlidingFitResult::GetLayer(), lar_content::TwoDSlidingFitResult::GetLayerFitResultMap(), lar_content::TwoDSlidingFitResult::GetLocalPosition(), and m_showerRadius.

Referenced by GetViewShowerRegionVariables().

 {
     // Find initial shower direction from spine fit
     float lShowerStart(0.f), tShowerStart(0.f);
     spineFit.GetLocalPosition(showerStart2D, lShowerStart, tShowerStart);
 
     const LayerFitResultMap &layerFitResultMap(spineFit.GetLayerFitResultMap());
     const int showerStartLayer(spineFit.GetLayer(lShowerStart));
     int closestLayer(std::numeric_limits<int>::max());
 
     for (const auto &entry : layerFitResultMap)
     {
         if (std::fabs(entry.first - showerStartLayer) < std::fabs(entry.first - closestLayer))
             closestLayer = entry.first;
     }
 
     const float gradient(layerFitResultMap.at(closestLayer).GetGradient());
     CartesianVector showerDirection(0.f, 0.f, 0.f);
 
     spineFit.GetGlobalDirection(gradient, showerDirection);
 
     // Now find the shower
     const bool isDownstream(showerStart2D.GetZ() > nuVertex2D.GetZ());
     CaloHitList caloHitList;
 
     LArPfoHelper::GetCaloHits(pShowerPfo, hitType, caloHitList);
 
     for (const CaloHit *const pCaloHit : caloHitList)
     {
         const CartesianVector &hitPosition(pCaloHit->GetPositionVector());
         const CartesianVector displacement(hitPosition - showerStart2D);
         const float l(showerDirection.GetDotProduct(displacement));
         const float t(showerDirection.GetCrossProduct(displacement).GetMagnitude());
 
         if (((isDownstream && (l > 0.f)) || (!isDownstream && (l < 0.f))) && (t < m_showerRadius))
         {
             postShowerHitList.push_back(pCaloHit);
             postShowerPositions.push_back(pCaloHit->GetPositionVector());
         }
     }
 }

void lar_content::ShowerRegionFeatureTool::CalculateViewNuVertexConsistencyVariables	(	const TwoDSlidingFitResult &	spineFitResult,
		const pandora::CaloHitList &	postShowerHitList,
		const bool	isDownstream,
		const pandora::CartesianVector &	nuVertex2D,
		float &	nuVertexEnergyAsymmetry,
		float &	nuVertexEnergyWeightedMeanRadialDistance
	)

private

Evaluate the neutrino vertex consistency variables.

Parameters

spineHitList	the shower spine hit list
postShowerHitList	the collected shower region hit list
isDownstream	whether the shower direction is downstream (in Z) of the neutrino vertex
nuVertex2D	the 2D neutrino vertex
nuVertexEnergyAsymmetry	the output neutrino vertex energy asymmetry
nuVertexEnergyWeightedMeanRadialDistance	the output neutrino vertex energy weighted mean radial distance

Definition at line 725 of file ConnectionPathwayFeatureTool.cc.

References f, lar_content::TwoDSlidingFitResult::GetGlobalMaxLayerDirection(), lar_content::TwoDSlidingFitResult::GetGlobalMinLayerDirection(), m_defaultFloat, and m_defaultRatio.

Referenced by GetViewShowerRegionVariables().

 {
     const CartesianVector &directionAxis(isDownstream ? spineFitResult.GetGlobalMinLayerDirection() : spineFitResult.GetGlobalMaxLayerDirection());
     const CartesianVector orthoAxis = directionAxis.GetCrossProduct(CartesianVector(0.f, 1.f, 0.f));
 
     float totalEnergy(0.f);
 
     // nuVertexEnergyAsymmetry
     nuVertexEnergyAsymmetry = 0.f;
 
     for (const CaloHit *const pCaloHit : postShowerHitList)
     {
         const float hitEnergy(std::fabs(pCaloHit->GetElectromagneticEnergy()));
 
         totalEnergy += hitEnergy;
 
         const CartesianVector position(pCaloHit->GetPositionVector() - nuVertex2D);
         const float thisL(orthoAxis.GetDotProduct(position));
 
         nuVertexEnergyAsymmetry += (thisL < 0.f) ? (-1.f * hitEnergy) : hitEnergy;
     }
 
     nuVertexEnergyAsymmetry = (totalEnergy < std::numeric_limits<float>::epsilon()) ? m_defaultRatio : (nuVertexEnergyAsymmetry / totalEnergy);
     nuVertexEnergyAsymmetry = std::fabs(nuVertexEnergyAsymmetry);
 
     // nuVertexEnergyWeightedMeanRadialDistance
     nuVertexEnergyWeightedMeanRadialDistance = 0.f;
 
     for (const CaloHit *const pCaloHit : postShowerHitList)
     {
         const CartesianVector position(pCaloHit->GetPositionVector() - nuVertex2D);
         const float hitEnergy(std::fabs(pCaloHit->GetElectromagneticEnergy()));
 
         nuVertexEnergyWeightedMeanRadialDistance += (directionAxis.GetCrossProduct(position).GetMagnitude() * hitEnergy);
     }
 
     nuVertexEnergyWeightedMeanRadialDistance =
         (totalEnergy < std::numeric_limits<float>::epsilon()) ? m_defaultFloat : nuVertexEnergyWeightedMeanRadialDistance / totalEnergy;
 }

void lar_content::ShowerRegionFeatureTool::CalculateViewOpeningAngle	(	const TwoDSlidingFitResult &	showerFitResult,
		const pandora::CaloHitList &	postShowerHitList,
		const pandora::CartesianVector &	showerStart2D,
		float &	openingAngle
	)

private

Calculate the opening angle of the shower region.

Parameters

showerFitResult	the shower region fit
postShowerHitList	the collected shower region hit list
showerStart2D	the 2D shower start position
openingAngle	the output opening angle

Definition at line 661 of file ConnectionPathwayFeatureTool.cc.

References f, lar_content::TwoDSlidingFitResult::GetGlobalMinLayerDirection(), lar_content::TwoDSlidingFitResult::GetLayerPitch(), m_edgeStep, and m_showerFitWindow.

Referenced by GetViewShowerRegionVariables().

 {
     try
     {
         const CartesianVector directionAxis = showerFitResult.GetGlobalMinLayerDirection();
         const CartesianVector orthoAxis = directionAxis.GetCrossProduct(CartesianVector(0.f, 1.f, 0.f));
 
         std::map<int, float> positiveEdges, negativeEdges;
 
         for (const CaloHit *const pCaloHit : postShowerHitList)
         {
             const CartesianVector position(pCaloHit->GetPositionVector() - showerStart2D);
             const float thisT(directionAxis.GetCrossProduct(position).GetMagnitude());
             const float thisL(directionAxis.GetDotProduct(position));
             const float orthoL(orthoAxis.GetDotProduct(position));
 
             std::map<int, float> &edgeMap(orthoL > 0.f ? positiveEdges : negativeEdges);
 
             const int lIndex(std::floor(thisL / m_edgeStep));
 
             edgeMap[lIndex] = (edgeMap.find(lIndex) == edgeMap.end() ? thisT : std::max(edgeMap[lIndex], thisT));
         }
 
         CartesianPointVector positiveEdgePositions, negativeEdgePositions;
 
         for (auto &entry : positiveEdges)
             positiveEdgePositions.push_back(CartesianVector(entry.second, 0.f, entry.first));
 
         for (auto &entry : negativeEdges)
             negativeEdgePositions.push_back(CartesianVector(entry.second, 0.f, entry.first));
 
         const TwoDSlidingFitResult positiveEdgeFit(&positiveEdgePositions, m_showerFitWindow, showerFitResult.GetLayerPitch());
         const TwoDSlidingFitResult negativeEdgeFit(&negativeEdgePositions, m_showerFitWindow, showerFitResult.GetLayerPitch());
 
         const CartesianVector positiveMinLayer(positiveEdgeFit.GetGlobalMinLayerPosition());
         const CartesianVector positiveMaxLayer(positiveEdgeFit.GetGlobalMaxLayerPosition());
         const CartesianVector negativeMinLayer(negativeEdgeFit.GetGlobalMinLayerPosition());
         const CartesianVector negativeMaxLayer(negativeEdgeFit.GetGlobalMaxLayerPosition());
 
         const CartesianVector globalPositiveMinLayer(
             showerStart2D + (directionAxis * positiveMinLayer.GetZ()) + (orthoAxis * positiveMinLayer.GetX()));
         const CartesianVector globalPositiveMaxLayer(
             showerStart2D + (directionAxis * positiveMaxLayer.GetZ()) + (orthoAxis * positiveMaxLayer.GetX()));
         const CartesianVector globalNegativeMinLayer(
             showerStart2D + (directionAxis * negativeMinLayer.GetZ()) - (orthoAxis * negativeMinLayer.GetX()));
         const CartesianVector globalNegativeMaxLayer(
             showerStart2D + (directionAxis * negativeMaxLayer.GetZ()) - (orthoAxis * negativeMaxLayer.GetX()));
 
         const CartesianVector positiveEdgeVector((globalPositiveMaxLayer - globalPositiveMinLayer).GetUnitVector());
         const CartesianVector negativeEdgeVector((globalNegativeMaxLayer - globalNegativeMinLayer).GetUnitVector());
 
         const float positiveOpeningAngle = directionAxis.GetOpeningAngle(positiveEdgeVector) * 180.f / M_PI;
         const float negativeOpeningAngle = directionAxis.GetOpeningAngle(negativeEdgeVector) * 180.f / M_PI;
 
         openingAngle = std::max(positiveOpeningAngle, negativeOpeningAngle);
     }
     catch (...)
     {
     }
 }

void lar_content::ShowerRegionFeatureTool::CalculateViewScatterAngle	(	const pandora::CartesianVector &	nuVertex2D,
		const TwoDSlidingFitResult &	spineFitResult,
		const pandora::CartesianVector &	showerStart2D,
		const TwoDSlidingFitResult &	showerFitResult,
		float &	scatterAngle
	)

private

Calculate the connection pathway-shower region scatter angle.

Parameters

nuVertex2D	the 2D neutrino vertex
spineFitResult	the shower spine fit
showerStart2D	the 2D shower start position
showerFitResult	the shower region fit
scatterAngle	the output scatter angle

Definition at line 644 of file ConnectionPathwayFeatureTool.cc.

References lar_content::TwoDSlidingFitResult::GetGlobalMaxLayerDirection(), lar_content::TwoDSlidingFitResult::GetGlobalMaxLayerPosition(), lar_content::TwoDSlidingFitResult::GetGlobalMinLayerDirection(), and lar_content::TwoDSlidingFitResult::GetGlobalMinLayerPosition().

Referenced by GetViewShowerRegionVariables().

 {
     const bool isDownstream(showerStart2D.GetZ() > nuVertex2D.GetZ());
     const CartesianVector streamCorrectedConnectionPathwayDirection(
         isDownstream ? spineFitResult.GetGlobalMinLayerDirection() : spineFitResult.GetGlobalMaxLayerDirection() * (-1.f));
 
     const bool isShowerDownstream((showerStart2D - showerFitResult.GetGlobalMinLayerPosition()).GetMagnitude() <
         (showerStart2D - showerFitResult.GetGlobalMaxLayerPosition()).GetMagnitude());
     const CartesianVector streamCorrectedShowerDirection(
         isShowerDownstream ? showerFitResult.GetGlobalMinLayerDirection() : showerFitResult.GetGlobalMaxLayerDirection() * (-1.f));
 
     scatterAngle = streamCorrectedConnectionPathwayDirection.GetOpeningAngle(streamCorrectedShowerDirection) * 180.f / M_PI;
 }

void lar_content::ShowerRegionFeatureTool::CalculateViewShowerStartConsistencyVariables	(	const TwoDSlidingFitResult &	showerFitResult,
		const pandora::CaloHitList &	postShowerHitList,
		const bool	isShowerDownstream,
		float &	showerStartEnergyAsymmetry,
		float &	showerStartMoliereRadius
	)

private

Evaluate the shower start consistency variables.

Parameters

showerFitResult	the shower fit result
postShowerHitList	the collected shower region hit list
isShowerDownstream	whether the shower direction is downstream (in Z) of the neutrino vertex
showerStartEnergyAsymmetry	the output shower start energy asymmetry
showerStartMoliereRadius	the output shower start moliere radius

Definition at line 768 of file ConnectionPathwayFeatureTool.cc.

References f, lar_content::TwoDSlidingFitResult::GetGlobalMaxLayerDirection(), lar_content::TwoDSlidingFitResult::GetGlobalMaxLayerPosition(), lar_content::TwoDSlidingFitResult::GetGlobalMinLayerDirection(), lar_content::TwoDSlidingFitResult::GetGlobalMinLayerPosition(), m_defaultFloat, m_defaultRatio, and m_moliereFraction.

Referenced by GetViewShowerRegionVariables().

 {
     const CartesianVector fitShowerStart(isShowerDownstream ? showerFitResult.GetGlobalMinLayerPosition() : showerFitResult.GetGlobalMaxLayerPosition());
     const CartesianVector directionAxis(
         isShowerDownstream ? showerFitResult.GetGlobalMinLayerDirection() : showerFitResult.GetGlobalMaxLayerDirection());
     const CartesianVector orthoAxis = directionAxis.GetCrossProduct(CartesianVector(0.f, 1.f, 0.f));
 
     float totalEnergy(0.f);
 
     // showerStartEnergyAsymmetry
     showerStartEnergyAsymmetry = 0.f;
 
     for (const CaloHit *const pCaloHit : postShowerHitList)
     {
         const float hitEnergy(std::fabs(pCaloHit->GetElectromagneticEnergy()));
 
         totalEnergy += hitEnergy;
 
         const CartesianVector position(pCaloHit->GetPositionVector() - fitShowerStart);
         const float thisL(orthoAxis.GetDotProduct(position));
 
         showerStartEnergyAsymmetry += (thisL < 0.f) ? (-1.f * hitEnergy) : hitEnergy;
     }
 
     showerStartEnergyAsymmetry = (totalEnergy < std::numeric_limits<float>::epsilon()) ? m_defaultRatio : (showerStartEnergyAsymmetry / totalEnergy);
     showerStartEnergyAsymmetry = std::fabs(showerStartEnergyAsymmetry);
 
     // showerStartMoliereRadius
     showerStartMoliereRadius = m_defaultFloat;
 
     CaloHitVector showerStartPostShowerHitVector(postShowerHitList.begin(), postShowerHitList.end());
 
     std::sort(showerStartPostShowerHitVector.begin(), showerStartPostShowerHitVector.end(),
         [&fitShowerStart, &directionAxis](const CaloHit *const pCaloHitA, const CaloHit *const pCaloHitB) -> bool
         {
             const CartesianVector positionA(pCaloHitA->GetPositionVector() - fitShowerStart);
             const CartesianVector positionB(pCaloHitB->GetPositionVector() - fitShowerStart);
 
             const float tA(directionAxis.GetCrossProduct(positionA).GetMagnitude());
             const float tB(directionAxis.GetCrossProduct(positionB).GetMagnitude());
 
             return tA < tB;
         });
 
     float showerStartRunningEnergySum(0.f);
 
     for (const CaloHit *const pCaloHit : showerStartPostShowerHitVector)
     {
         const float hitEnergy(std::fabs(pCaloHit->GetElectromagneticEnergy()));
         showerStartRunningEnergySum += hitEnergy;
 
         if ((showerStartRunningEnergySum / totalEnergy) > m_moliereFraction)
         {
             const CartesianVector position(pCaloHit->GetPositionVector() - fitShowerStart);
             showerStartMoliereRadius = directionAxis.GetCrossProduct(position).GetMagnitude();
             break;
         }
     }
 }

void lar_content::ShowerRegionFeatureTool::GetShowerHitVariables	(	const pandora::CaloHitList &	spineHitList,
		const pandora::CaloHitList &	postShowerHitList,
		const pandora::ParticleFlowObject *const	pShowerPfo,
		const pandora::HitType	hitType,
		float &	nHits,
		float &	foundHitRatio
	)

private

Evaluate the variables associated with the shower region hit multiplicity.

Parameters

spineHitList	the shower spine hit list
postShowerHitList	the collected shower region hit list
pShowerPfo	the shower pfo
hitType	the 2D view
nHits	the output number of shower region hits
foundHitRatio	the output found hit ratio

Definition at line 624 of file ConnectionPathwayFeatureTool.cc.

References lar_content::LArPfoHelper::GetCaloHits().

Referenced by GetViewShowerRegionVariables().

 {
     int foundHits(spineHitList.size());
 
     for (const CaloHit *const pCaloHit : postShowerHitList)
     {
         if (std::find(spineHitList.begin(), spineHitList.end(), pCaloHit) == spineHitList.end())
             ++foundHits;
     }
 
     CaloHitList viewHitList;
     LArPfoHelper::GetCaloHits(pShowerPfo, hitType, viewHitList);
 
     foundHitRatio = static_cast<float>(foundHits) / static_cast<float>(viewHitList.size());
     nHits = postShowerHitList.size();
 }

void lar_content::ShowerRegionFeatureTool::GetViewShowerRegionVariables	(	const pandora::Algorithm *const	pAlgorithm,
		const pandora::ParticleFlowObject *const	pShowerPfo,
		const pandora::CartesianVector &	nuVertex3D,
		const ProtoShowerMatch &	protoShowerMatch,
		const pandora::HitType	hitType,
		const pandora::CartesianVector &	showerStart3D,
		float &	nHits,
		float &	foundHitRatio,
		float &	scatterAngle,
		float &	openingAngle,
		float &	nuVertexEnergyAsymmetry,
		float &	nuVertexEnergyWeightedMeanRadialDistance,
		float &	showerStartEnergyAsymmetry,
		float &	showerStartMoliereRadius
	)

private

Calculate the shower region variables for the input view.

Parameters

pAlgorithm	the algorithm
pShowerPfo	the shower pfo
nuVertex3D	the 3D neutrino vertex
protoShowerMatch	the ProtoShower match
hitType	the 2D view
showerStart3D	the 3D shower start position
nHits	the output number of shower region hits
foundHitRatio	the output found hit ratio
scatterAngle	the output scatter angle
openingAngle	the output opening angle
nuVertexEnergyAsymmetry	the output neutrino vertex energy asymmetry
nuVertexEnergyWeightedMeanRadialDistance	the output neutrino vertex energy weighted mean radial distance
showerStartEnergyAsymmetry	the output shower start energy asymmetry
showerStartMoliereRadius	the output shower start moliere radius

Definition at line 502 of file ConnectionPathwayFeatureTool.cc.

References BuildViewShower(), CalculateViewNuVertexConsistencyVariables(), CalculateViewOpeningAngle(), CalculateViewScatterAngle(), CalculateViewShowerStartConsistencyVariables(), f, lar_content::LArPfoHelper::GetCaloHits(), lar_content::TwoDSlidingFitResult::GetGlobalMaxLayerPosition(), lar_content::TwoDSlidingFitResult::GetGlobalMinLayerDirection(), lar_content::TwoDSlidingFitResult::GetGlobalMinLayerPosition(), lar_content::ProtoShowerMatch::GetProtoShowerU(), lar_content::ProtoShowerMatch::GetProtoShowerV(), lar_content::ProtoShowerMatch::GetProtoShowerW(), GetShowerHitVariables(), lar_content::ProtoShower::GetSpineHitList(), lar_content::LArGeometryHelper::GetWirePitch(), m_showerFitWindow, m_showerRadius, m_spineFitWindow, and lar_content::LArGeometryHelper::ProjectPosition().

Referenced by ShowerRegionFeatureTool().

 {
     CaloHitList viewHitList;
     LArPfoHelper::GetCaloHits(pShowerPfo, hitType, viewHitList);
 
     const CartesianVector nuVertex2D(LArGeometryHelper::ProjectPosition(pAlgorithm->GetPandora(), nuVertex3D, hitType));
     const CartesianVector showerStart2D(LArGeometryHelper::ProjectPosition(pAlgorithm->GetPandora(), showerStart3D, hitType));
     const bool isDownstream(showerStart2D.GetZ() > nuVertex2D.GetZ());
 
     try
     {
         // Fit the spine to get the initial shower direction
         const CaloHitList &spineHitList(hitType == TPC_VIEW_U ? protoShowerMatch.GetProtoShowerU().GetSpineHitList()
                                                               : (hitType == TPC_VIEW_V ? protoShowerMatch.GetProtoShowerV().GetSpineHitList()
                                                                                        : protoShowerMatch.GetProtoShowerW().GetSpineHitList()));
 
         CartesianPointVector spinePositions;
         for (const CaloHit *const pCaloHit : spineHitList)
             spinePositions.push_back(pCaloHit->GetPositionVector());
 
         const TwoDSlidingFitResult spineFitResult(&spinePositions, m_spineFitWindow, LArGeometryHelper::GetWirePitch(pAlgorithm->GetPandora(), hitType));
 
         // Now can build the shower
         CaloHitList postShowerHitList;
         CartesianPointVector postShowerPositions;
 
         this->BuildViewShower(pShowerPfo, spineFitResult, hitType, showerStart2D, nuVertex2D, postShowerHitList, postShowerPositions);
 
         this->GetShowerHitVariables(spineHitList, postShowerHitList, pShowerPfo, hitType, nHits, foundHitRatio);
 
         // Fit the shower
         TwoDSlidingFitResult showerFitResult(
             &postShowerPositions, m_showerFitWindow, LArGeometryHelper::GetWirePitch(pAlgorithm->GetPandora(), hitType));
 
         const bool isShowerDownstream((showerStart2D - showerFitResult.GetGlobalMinLayerPosition()).GetMagnitude() <
             (showerStart2D - showerFitResult.GetGlobalMaxLayerPosition()).GetMagnitude());
 
         // Collect more hits
         for (const CaloHit *const pCaloHit : viewHitList)
         {
             if (std::find(postShowerHitList.begin(), postShowerHitList.end(), pCaloHit) != postShowerHitList.end())
                 continue;
 
             const CartesianVector &hitPosition(pCaloHit->GetPositionVector());
             const CartesianVector displacement(hitPosition - showerStart2D);
             const float l(showerFitResult.GetGlobalMinLayerDirection().GetDotProduct(displacement));
             const float t(showerFitResult.GetGlobalMinLayerDirection().GetCrossProduct(displacement).GetMagnitude());
 
             if (((isDownstream && (l > 0.f)) || (!isDownstream && (l < 0.f))) && (t < m_showerRadius))
             {
                 postShowerHitList.push_back(pCaloHit);
                 postShowerPositions.push_back(pCaloHit->GetPositionVector());
             }
         }
 
         this->GetShowerHitVariables(spineHitList, postShowerHitList, pShowerPfo, hitType, nHits, foundHitRatio);
 
         this->CalculateViewScatterAngle(nuVertex2D, spineFitResult, showerStart2D, showerFitResult, scatterAngle);
 
         this->CalculateViewOpeningAngle(showerFitResult, postShowerHitList, showerStart2D, openingAngle);
 
         this->CalculateViewNuVertexConsistencyVariables(
             spineFitResult, postShowerHitList, isDownstream, nuVertex2D, nuVertexEnergyAsymmetry, nuVertexEnergyWeightedMeanRadialDistance);
 
         this->CalculateViewShowerStartConsistencyVariables(
             showerFitResult, postShowerHitList, isShowerDownstream, showerStartEnergyAsymmetry, showerStartMoliereRadius);
     }
     catch (...)
     {
     }
 }

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

private

Definition at line 831 of file ConnectionPathwayFeatureTool.cc.

References m_defaultFloat, m_defaultRatio, m_edgeStep, m_maxFoundHitRatioLimit, m_maxNHitsLimit, m_maxNuVertexEnergyWeightedMeanRadialDistanceLimit, m_maxOpeningAngleLimit, m_maxScatterAngleLimit, m_minShowerStartMoliereRadiusLimit, m_moliereFraction, m_showerFitWindow, m_showerRadius, 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, "DefaultRatio", m_defaultRatio));
 
     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, "ShowerRadius", m_showerRadius));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "ShowerFitWindow", m_showerFitWindow));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "EdgeStep", m_edgeStep));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MoliereFraction", m_moliereFraction));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MaxNHitsLimit", m_maxNHitsLimit));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MaxFoundHitRatioLimit", m_maxFoundHitRatioLimit));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MaxScatterAngleLimit", m_maxScatterAngleLimit));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MaxOpeningAngleLimit", m_maxOpeningAngleLimit));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
         XmlHelper::ReadValue(xmlHandle, "MaxNuVertexEnergyWeightedMeanRadialDistanceLimit", m_maxNuVertexEnergyWeightedMeanRadialDistanceLimit));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
         XmlHelper::ReadValue(xmlHandle, "MinShowerStartMoliereRadiusLimit", m_minShowerStartMoliereRadiusLimit));
 
     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::ShowerRegionFeatureTool::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 ShowerRegionFeatureTool().

void lar_content::ShowerRegionFeatureTool::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::ShowerRegionFeatureTool::m_defaultFloat

private

Default float value.

Definition at line 248 of file ConnectionPathwayFeatureTool.h.

Referenced by CalculateViewNuVertexConsistencyVariables(), CalculateViewShowerStartConsistencyVariables(), ReadSettings(), and ShowerRegionFeatureTool().

float lar_content::ShowerRegionFeatureTool::m_defaultRatio

private

Default float value for ratios.

Definition at line 249 of file ConnectionPathwayFeatureTool.h.

Referenced by CalculateViewNuVertexConsistencyVariables(), CalculateViewShowerStartConsistencyVariables(), ReadSettings(), and ShowerRegionFeatureTool().

float lar_content::ShowerRegionFeatureTool::m_edgeStep

private

The binning of the shower boundaries.

Definition at line 253 of file ConnectionPathwayFeatureTool.h.

Referenced by CalculateViewOpeningAngle(), and ReadSettings().

float lar_content::ShowerRegionFeatureTool::m_maxFoundHitRatioLimit

private

maxFoundHitRatio max. limit

Definition at line 256 of file ConnectionPathwayFeatureTool.h.

Referenced by ReadSettings(), and ShowerRegionFeatureTool().

float lar_content::ShowerRegionFeatureTool::m_maxNHitsLimit

private

maxNHits max. limit

Definition at line 255 of file ConnectionPathwayFeatureTool.h.

Referenced by ReadSettings(), and ShowerRegionFeatureTool().

float lar_content::ShowerRegionFeatureTool::m_maxNuVertexEnergyWeightedMeanRadialDistanceLimit

private

maxNuVertexEnergyWeightedMeanRadialDistance max. limit

Definition at line 259 of file ConnectionPathwayFeatureTool.h.

Referenced by ReadSettings(), and ShowerRegionFeatureTool().

float lar_content::ShowerRegionFeatureTool::m_maxOpeningAngleLimit

private

maxOpeningAngle max. limit

Definition at line 258 of file ConnectionPathwayFeatureTool.h.

Referenced by ReadSettings(), and ShowerRegionFeatureTool().

float lar_content::ShowerRegionFeatureTool::m_maxScatterAngleLimit

private

maxScatterAngle max. limit

Definition at line 257 of file ConnectionPathwayFeatureTool.h.

Referenced by ReadSettings(), and ShowerRegionFeatureTool().

float lar_content::ShowerRegionFeatureTool::m_minShowerStartMoliereRadiusLimit

private

minShowerStartMoliereRadius max. limit

Definition at line 260 of file ConnectionPathwayFeatureTool.h.

Referenced by ReadSettings(), and ShowerRegionFeatureTool().

float lar_content::ShowerRegionFeatureTool::m_moliereFraction

private

The energy fraction which corresponds to minShowerStartMoliereRadius.

Definition at line 254 of file ConnectionPathwayFeatureTool.h.

Referenced by CalculateViewShowerStartConsistencyVariables(), and ReadSettings().

unsigned int lar_content::ShowerRegionFeatureTool::m_showerFitWindow

private

The shower fit window.

Definition at line 252 of file ConnectionPathwayFeatureTool.h.

Referenced by CalculateViewOpeningAngle(), GetViewShowerRegionVariables(), and ReadSettings().

float lar_content::ShowerRegionFeatureTool::m_showerRadius

private

The max. separation distance between a shower region hit and the shower core.

Definition at line 251 of file ConnectionPathwayFeatureTool.h.

Referenced by BuildViewShower(), GetViewShowerRegionVariables(), and ReadSettings().

unsigned int lar_content::ShowerRegionFeatureTool::m_spineFitWindow

private

The spine fit window.

Definition at line 250 of file ConnectionPathwayFeatureTool.h.

Referenced by GetViewShowerRegionVariables(), 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