TransverseTrackHitsBaseTool class. More...

#include <TransverseTrackHitsBaseTool.h>

Inheritance diagram for lar_content::TransverseTrackHitsBaseTool:

Protected Member Functions
virtual void	GetTransverseTrackHit3D (const MatchedSlidingFitMap &matchedSlidingFitMap, ProtoHit &protoHit) const =0
	Get the three dimensional position using a provided two dimensional calo hit and sliding linear fits in the other two views. More...

virtual void	GetTrackHits3D (const pandora::CaloHitVector &inputTwoDHits, const MatchedSlidingFitMap &matchedSlidingFitMap, ProtoHitVector &protoHitVector) const
	Calculate 3D hits from an input list of 2D hits. More...

virtual void	AddTransverseChi2 (const MatchedSlidingFitMap &matchedSlidingFitMap, ProtoHit &protoHit) const
	Calculate an additional contribution to the chi-squared based on the steepness of the track. More...

Protected Member Functions inherited from lar_content::TrackHitsBaseTool
virtual void	BuildSlidingFitMap (const pandora::ParticleFlowObject *const pPfo, MatchedSlidingFitMap &matchedSlidingFitMap) const
	Calculate sliding fit results for clusters from each view. More...

virtual pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)

Protected Member Functions inherited from lar_content::HitCreationBaseTool
virtual void	GetBestPosition3D (const pandora::HitType hitType1, const pandora::HitType hitType2, const pandora::CartesianPointVector &fitPositionList1, const pandora::CartesianPointVector &fitPositionList2, ProtoHit &protoHit) const
	Get the three dimensional position using a provided two dimensional calo hit and candidate fit positions from the other two views. More...

virtual void	GetBestPosition3D (const pandora::HitType hitType1, const pandora::HitType hitType2, const pandora::CartesianVector &fitPosition1, const pandora::CartesianVector &fitPosition2, ProtoHit &protoHit) const
	Get the three dimensional position using a provided two dimensional calo hit and candidate fit positions from the other two views. More...

virtual void	GetBestPosition3D (const pandora::HitType hitType, const pandora::CartesianVector &fitPosition, ProtoHit &protoHit) const
	Get the three dimensional position using a provided two dimensional calo hit and a candidate fit position from another view. More...

Private Member Functions
double	GetTransverseChi2 (const pandora::CartesianVector &position2D, const TwoDSlidingFitResult &fitResult) const
	Calculate an additional contribution to the chi-squared based on the steepness of the track. More...

Additional Inherited Members
Public Types inherited from lar_content::HitCreationBaseTool
typedef ThreeDHitCreationAlgorithm::ProtoHit	ProtoHit

typedef ThreeDHitCreationAlgorithm::ProtoHitVector	ProtoHitVector

typedef ThreeDHitCreationAlgorithm::TrajectorySample	TrajectorySample

Public Member Functions inherited from lar_content::TrackHitsBaseTool
	TrackHitsBaseTool ()
	Default constructor. More...

virtual void	Run (ThreeDHitCreationAlgorithm const pAlgorithm, const pandora::ParticleFlowObject const pPfo, const pandora::CaloHitVector &inputTwoDHits, ProtoHitVector &protoHitVector)
	Run the algorithm tool. More...

Public Member Functions inherited from lar_content::HitCreationBaseTool
	HitCreationBaseTool ()
	Default constructor. More...

virtual	~HitCreationBaseTool ()
	Destructor. More...

Protected Types inherited from lar_content::TrackHitsBaseTool
typedef std::map < pandora::HitType, TwoDSlidingFitResult >	MatchedSlidingFitMap

Protected Attributes inherited from lar_content::TrackHitsBaseTool
unsigned int	m_minViews
	The minimum number of views required for building hits. More...

unsigned int	m_slidingFitWindow
	The layer window for the sliding linear fits. More...

Protected Attributes inherited from lar_content::HitCreationBaseTool
double	m_sigmaX2
	The sigmaX squared value, for calculation of chi2 deltaX term. More...

double	m_chiSquaredCut
	The chi squared cut (accept only values below the cut value) More...

Detailed Description

TransverseTrackHitsBaseTool class.

Definition at line 19 of file TransverseTrackHitsBaseTool.h.

Member Function Documentation

void lar_content::TransverseTrackHitsBaseTool::AddTransverseChi2	(	const MatchedSlidingFitMap &	matchedSlidingFitMap,
		ProtoHit &	protoHit
	)		const

protectedvirtual

Calculate an additional contribution to the chi-squared based on the steepness of the track.

Parameters

matchedSlidingFitMap	map of sliding fit results from each view
protoHit	to receive the modified proto hit

Definition at line 43 of file TransverseTrackHitsBaseTool.cc.

 {
     // TODO Develop a proper treatment of the |dz/dx| * sigmaX uncertainty
     double chiSquared(protoHit.GetChi2());
     const HitType inputHitType(protoHit.GetParentCaloHit2D()->GetHitType());
     const CartesianVector &inputPosition3D(protoHit.GetPosition3D());
 
     for (const MatchedSlidingFitMap::value_type &mapEntry : matchedSlidingFitMap)
     {
         if (mapEntry.first == inputHitType)
             continue;
 
         const CartesianVector inputPosition2D(LArGeometryHelper::ProjectPosition(this->GetPandora(), inputPosition3D, mapEntry.first));
         chiSquared += this->GetTransverseChi2(inputPosition2D, mapEntry.second);
     }
 
     protoHit.SetPosition3D(inputPosition3D, chiSquared);
 }

void lar_content::TransverseTrackHitsBaseTool::GetTrackHits3D	(	const pandora::CaloHitVector &	inputTwoDHits,
		const MatchedSlidingFitMap &	matchedSlidingFitMap,
		ProtoHitVector &	protoHitVector
	)		const

protectedvirtual

Calculate 3D hits from an input list of 2D hits.

Parameters

pAlgorithm	the hit creation algorithm
inputTwoDHits	the input vector of 2D hits
matchedSlidingFitMap	the group of sliding fit results
protoHitVector	to receive the new three dimensional proto hits

Implements lar_content::TrackHitsBaseTool.

Definition at line 21 of file TransverseTrackHitsBaseTool.cc.

 {
     for (const CaloHit *const pCaloHit2D : inputTwoDHits)
     {
         try
         {
             ProtoHit protoHit(pCaloHit2D);
             this->GetTransverseTrackHit3D(matchedSlidingFitMap, protoHit);
             this->AddTransverseChi2(matchedSlidingFitMap, protoHit);
 
             if (protoHit.IsPositionSet() && (protoHit.GetChi2() < m_chiSquaredCut))
                 protoHitVector.push_back(protoHit);
         }
         catch (StatusCodeException &)
         {
         }
     }
 }

double lar_content::TransverseTrackHitsBaseTool::GetTransverseChi2	(	const pandora::CartesianVector &	position2D,
		const TwoDSlidingFitResult &	fitResult
	)		const

private

Calculate an additional contribution to the chi-squared based on the steepness of the track.

Parameters

position2D	the calculated two dimensional position
fitResult	the sliding fit to the track

Definition at line 64 of file TransverseTrackHitsBaseTool.cc.

 {
     const float minLayerX(fitResult.GetGlobalMinLayerPosition().GetX());
     const float maxLayerX(fitResult.GetGlobalMaxLayerPosition().GetX());
 
     const float minX(std::min(minLayerX, maxLayerX));
     const float maxX(std::max(minLayerX, maxLayerX));
 
     if (((position2D.GetX() - minX) > -std::numeric_limits<float>::epsilon()) && ((position2D.GetX() - maxX) < +std::numeric_limits<float>::epsilon()))
         return 0.;
 
     const float minLayerDistanceSquared((position2D - fitResult.GetGlobalMinLayerPosition()).GetMagnitudeSquared());
     const float maxLayerDistanceSquared((position2D - fitResult.GetGlobalMaxLayerPosition()).GetMagnitudeSquared());
 
     double px(0.), pz(0.);
 
     if (minLayerDistanceSquared < maxLayerDistanceSquared)
     {
         px = fitResult.GetGlobalMinLayerDirection().GetX();
         pz = fitResult.GetGlobalMinLayerDirection().GetZ();
     }
     else if (maxLayerDistanceSquared < minLayerDistanceSquared)
     {
         px = fitResult.GetGlobalMaxLayerDirection().GetX();
         pz = fitResult.GetGlobalMaxLayerDirection().GetZ();
     }
 
     if (std::fabs(px) > std::numeric_limits<double>::epsilon())
     {
         return (0.5 * (pz * pz) / (px * px)); // TODO: Figure out the appropriate scale factor
     }
 
     throw StatusCodeException(STATUS_CODE_NOT_FOUND);
 }

virtual void lar_content::TransverseTrackHitsBaseTool::GetTransverseTrackHit3D	(	const MatchedSlidingFitMap &	matchedSlidingFitMap,
		ProtoHit &	protoHit
	)		const

protectedpure virtual

Get the three dimensional position using a provided two dimensional calo hit and sliding linear fits in the other two views.

Parameters

matchedSlidingFitMap	map of sliding fit results from each view
protoHit	to receive the populated proto hit

Implemented in lar_content::ClearTransverseTrackHitsTool, and lar_content::MultiValuedTransverseTrackHitsTool.

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

Protected Member Functions

Private Member Functions

Additional Inherited Members

Detailed Description

Member Function Documentation