MultiValuedLongitudinalTrackHitsTool class. More...

#include "MultiValuedLongitudinalTrackHitsTool.h"

Inheritance diagram for lar_content::MultiValuedLongitudinalTrackHitsTool:

Public Types
typedef ThreeDHitCreationAlgorithm::ProtoHit	ProtoHit

typedef ThreeDHitCreationAlgorithm::ProtoHitVector	ProtoHitVector

typedef ThreeDHitCreationAlgorithm::TrajectorySample	TrajectorySample

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

Protected Types
typedef std::map< pandora::HitType, TwoDSlidingFitResult >	MatchedSlidingFitMap

Protected Member Functions
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 pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)

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

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...

Protected Attributes
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...

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...

Private Member Functions
void	GetLongitudinalTrackHit3D (const MatchedSlidingFitMap &matchedSlidingFitMap, const pandora::CartesianVector &vtx3D, const pandora::CartesianVector &end3D, ProtoHit &protoHit) const
	Get the three dimensional position using a provided two dimensional calo hit and sliding linear fits in the other two views. More...

Detailed Description

MultiValuedLongitudinalTrackHitsTool class.

Definition at line 19 of file MultiValuedLongitudinalTrackHitsTool.h.

Member Typedef Documentation

typedef std::map<pandora::HitType, TwoDSlidingFitResult> lar_content::TrackHitsBaseTool::MatchedSlidingFitMap

protectedinherited

Definition at line 35 of file TrackHitsBaseTool.h.

typedef ThreeDHitCreationAlgorithm::ProtoHit lar_content::HitCreationBaseTool::ProtoHit

inherited

Definition at line 24 of file HitCreationBaseTool.h.

typedef ThreeDHitCreationAlgorithm::ProtoHitVector lar_content::HitCreationBaseTool::ProtoHitVector

inherited

Definition at line 25 of file HitCreationBaseTool.h.

typedef ThreeDHitCreationAlgorithm::TrajectorySample lar_content::HitCreationBaseTool::TrajectorySample

inherited

Definition at line 26 of file HitCreationBaseTool.h.

Member Function Documentation

void lar_content::TrackHitsBaseTool::BuildSlidingFitMap	(	const pandora::ParticleFlowObject *const	pPfo,
		MatchedSlidingFitMap &	matchedSlidingFitMap
	)		const

protectedvirtualinherited

Calculate sliding fit results for clusters from each view.

Parameters

pPfo	the input particle flow object
matchedSlidingFitMap	the group of sliding fit results

Definition at line 57 of file TrackHitsBaseTool.cc.

References lar_content::LArClusterHelper::GetClusterHitType(), lar_content::LArGeometryHelper::GetWirePitch(), lar_content::TrackHitsBaseTool::m_slidingFitWindow, and lar_content::LArClusterHelper::SortByNHits().

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

 {
     const ClusterList &pfoClusterList(pPfo->GetClusterList());
 
     ClusterVector pfoClusterVector;
     pfoClusterVector.insert(pfoClusterVector.end(), pfoClusterList.begin(), pfoClusterList.end());
     std::sort(pfoClusterVector.begin(), pfoClusterVector.end(), LArClusterHelper::SortByNHits);
 
     for (const Cluster *const pCluster : pfoClusterVector)
     {
         const HitType hitType(LArClusterHelper::GetClusterHitType(pCluster));
 
         if (TPC_3D == hitType)
             continue;
 
         if (matchedSlidingFitMap.end() != matchedSlidingFitMap.find(hitType))
             continue;
 
         try
         {
             const float slidingFitPitch(LArGeometryHelper::GetWirePitch(this->GetPandora(), LArClusterHelper::GetClusterHitType(pCluster)));
             const TwoDSlidingFitResult slidingFitResult(pCluster, m_slidingFitWindow, slidingFitPitch);
 
             if (!matchedSlidingFitMap.insert(MatchedSlidingFitMap::value_type(hitType, slidingFitResult)).second)
                 throw StatusCodeException(STATUS_CODE_FAILURE);
         }
         catch (StatusCodeException &statusCodeException)
         {
             if (STATUS_CODE_FAILURE == statusCodeException.GetStatusCode())
                 throw statusCodeException;
         }
     }
 }

virtual void lar_content::HitCreationBaseTool::GetBestPosition3D	(	const pandora::HitType	hitType1,
		const pandora::HitType	hitType2,
		const pandora::CartesianPointVector &	fitPositionList1,
		const pandora::CartesianPointVector &	fitPositionList2,
		ProtoHit &	protoHit
	)		const

protectedvirtualinherited

Get the three dimensional position using a provided two dimensional calo hit and candidate fit positions from the other two views.

Parameters

hitType1	the hit type identifying the first view
hitType2	the hit type identifying the second view
fitPositionList1	the candidate sliding fit position in the first view
fitPositionList2	the candidate sliding fit position in the second view
protoHit	to receive the populated proto hit

Referenced by lar_content::ThreeViewShowerHitsTool::GetShowerHit3D(), and lar_content::HitCreationBaseTool::~HitCreationBaseTool().

virtual void lar_content::HitCreationBaseTool::GetBestPosition3D	(	const pandora::HitType	hitType1,
		const pandora::HitType	hitType2,
		const pandora::CartesianVector &	fitPosition1,
		const pandora::CartesianVector &	fitPosition2,
		ProtoHit &	protoHit
	)		const

protectedvirtualinherited

Get the three dimensional position using a provided two dimensional calo hit and candidate fit positions from the other two views.

Parameters

hitType1	the hit type identifying the first view
hitType2	the hit type identifying the second view
fitPosition1	the candidate sliding fit position in the first view
fitPosition2	the candidate sliding fit position in the second view
protoHit	to receive the populated proto hit

virtual void lar_content::HitCreationBaseTool::GetBestPosition3D	(	const pandora::HitType	hitType,
		const pandora::CartesianVector &	fitPosition,
		ProtoHit &	protoHit
	)		const

protectedvirtualinherited

Get the three dimensional position using a provided two dimensional calo hit and a candidate fit position from another view.

Parameters

hitType	the hit type identifying the other view
fitPosition	the candidate sliding fit position in the other view
protoHit	to receive the populated proto hit

void lar_content::MultiValuedLongitudinalTrackHitsTool::GetLongitudinalTrackHit3D	(	const MatchedSlidingFitMap &	matchedSlidingFitMap,
		const pandora::CartesianVector &	vtx3D,
		const pandora::CartesianVector &	end3D,
		ProtoHit &	protoHit
	)		const

privatevirtual

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
vtx3D	the 3D vertex position
end3D	the 3D end position
protoHit	to receive the populated proto hit

Implements lar_content::LongitudinalTrackHitsBaseTool.

Definition at line 20 of file MultiValuedLongitudinalTrackHitsTool.cc.

References f, lar_content::TwoDSlidingFitResult::GetGlobalFitProjection(), and lar_content::ThreeDHitCreationAlgorithm::ProtoHit::GetParentCaloHit2D().

 {
     const CaloHit *const pCaloHit2D(protoHit.GetParentCaloHit2D());
     const HitType hitType(pCaloHit2D->GetHitType());
     const HitType hitType1((TPC_VIEW_U == hitType) ? TPC_VIEW_V : (TPC_VIEW_V == hitType) ? TPC_VIEW_W : TPC_VIEW_U);
     const HitType hitType2((TPC_VIEW_U == hitType) ? TPC_VIEW_W : (TPC_VIEW_V == hitType) ? TPC_VIEW_U : TPC_VIEW_V);
 
     const CartesianVector vtx2D(LArGeometryHelper::ProjectPosition(this->GetPandora(), vtx3D, hitType));
     const CartesianVector end2D(LArGeometryHelper::ProjectPosition(this->GetPandora(), end3D, hitType));
 
     if ((end2D - vtx2D).GetMagnitudeSquared() < std::numeric_limits<float>::epsilon())
         throw StatusCodeException(STATUS_CODE_NOT_FOUND);
 
     const float frac((end2D - vtx2D).GetDotProduct(pCaloHit2D->GetPositionVector() - vtx2D) / (end2D - vtx2D).GetMagnitudeSquared());
     const CartesianVector projection3D(vtx3D + (end3D - vtx3D) * frac);
 
     CartesianPointVector fitPositionList1, fitPositionList2;
 
     MatchedSlidingFitMap::const_iterator fIter1 = matchedSlidingFitMap.find(hitType1);
     if (matchedSlidingFitMap.end() != fIter1)
     {
         const TwoDSlidingFitResult &fitResult1 = fIter1->second;
         const CartesianVector position2D(LArGeometryHelper::ProjectPosition(this->GetPandora(), projection3D, hitType1));
 
         CartesianVector position1(0.f, 0.f, 0.f);
         const StatusCode statusCode(fitResult1.GetGlobalFitProjection(position2D, position1));
 
         if ((STATUS_CODE_SUCCESS != statusCode) && (STATUS_CODE_NOT_FOUND != statusCode))
             throw StatusCodeException(statusCode);
 
         if (STATUS_CODE_SUCCESS == statusCode)
             fitPositionList1.push_back(position1);
     }
 
     MatchedSlidingFitMap::const_iterator fIter2 = matchedSlidingFitMap.find(hitType2);
     if (matchedSlidingFitMap.end() != fIter2)
     {
         const TwoDSlidingFitResult &fitResult2 = fIter2->second;
         const CartesianVector position2D(LArGeometryHelper::ProjectPosition(this->GetPandora(), projection3D, hitType2));
 
         CartesianVector position2(0.f, 0.f, 0.f);
         const StatusCode statusCode(fitResult2.GetGlobalFitProjection(position2D, position2));
 
         if ((STATUS_CODE_SUCCESS != statusCode) && (STATUS_CODE_NOT_FOUND != statusCode))
             throw StatusCodeException(statusCode);
 
         if (STATUS_CODE_SUCCESS == statusCode)
             fitPositionList2.push_back(position2);
     }
 
     unsigned int nViews(1);
     if (fitPositionList1.size() > 0)
         ++nViews;
     if (fitPositionList2.size() > 0)
         ++nViews;
 
     if (nViews < m_minViews)
         throw StatusCodeException(STATUS_CODE_NOT_FOUND);
 
     this->GetBestPosition3D(hitType1, hitType2, fitPositionList1, fitPositionList2, protoHit);
 }

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

protectedvirtualinherited

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 29 of file LongitudinalTrackHitsBaseTool.cc.

References f, lar_content::ThreeDHitCreationAlgorithm::ProtoHit::GetChi2(), lar_content::LongitudinalTrackHitsBaseTool::GetLongitudinalTrackHit3D(), lar_content::LongitudinalTrackHitsBaseTool::GetVertexAndEndPositions(), lar_content::ThreeDHitCreationAlgorithm::ProtoHit::IsPositionSet(), and lar_content::HitCreationBaseTool::m_chiSquaredCut.

 {
     MatchedSlidingFitMap matchedSlidingFitMap;
     CartesianVector vtx3D(0.f, 0.f, 0.f), end3D(0.f, 0.f, 0.f);
     this->GetVertexAndEndPositions(inputSlidingFitMap, matchedSlidingFitMap, vtx3D, end3D);
 
     for (const CaloHit *const pCaloHit2D : inputTwoDHits)
     {
         try
         {
             ProtoHit protoHit(pCaloHit2D);
             this->GetLongitudinalTrackHit3D(matchedSlidingFitMap, vtx3D, end3D, protoHit);
 
             if (protoHit.IsPositionSet() && (protoHit.GetChi2() < m_chiSquaredCut))
                 protoHitVector.push_back(protoHit);
         }
         catch (StatusCodeException &)
         {
         }
     }
 }

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

protectedvirtualinherited

Reimplemented from lar_content::TrackHitsBaseTool.

Definition at line 214 of file LongitudinalTrackHitsBaseTool.cc.

References lar_content::LongitudinalTrackHitsBaseTool::m_minTrackLengthSquared, lar_content::LongitudinalTrackHitsBaseTool::m_vtxDisplacementCutSquared, and lar_content::TrackHitsBaseTool::ReadSettings().

 {
     float vtxDisplacementCut = std::sqrt(m_vtxDisplacementCutSquared);
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "VertexDisplacementCut", vtxDisplacementCut));
     m_vtxDisplacementCutSquared = vtxDisplacementCut * vtxDisplacementCut;
 
     float minTrackLength = std::sqrt(m_minTrackLengthSquared);
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MinTrackLength", minTrackLength));
     m_minTrackLengthSquared = minTrackLength * minTrackLength;
 
     return TrackHitsBaseTool::ReadSettings(xmlHandle);
 }

void lar_content::TrackHitsBaseTool::Run	(	ThreeDHitCreationAlgorithm *const	pAlgorithm,
		const pandora::ParticleFlowObject *const	pPfo,
		const pandora::CaloHitVector &	inputTwoDHits,
		ProtoHitVector &	protoHitVector
	)

virtualinherited

Run the algorithm tool.

Parameters

pAlgorithm	address of the calling algorithm
pPfo	the address of the pfo
inputTwoDHits	the vector of input two dimensional hits
protoHitVector	to receive the new three dimensional proto hits

Implements lar_content::HitCreationBaseTool.

Definition at line 31 of file TrackHitsBaseTool.cc.

References lar_content::TrackHitsBaseTool::BuildSlidingFitMap(), lar_content::TrackHitsBaseTool::GetTrackHits3D(), and lar_content::LArPfoHelper::IsTrack().

 {
     if (PandoraContentApi::GetSettings(*pAlgorithm)->ShouldDisplayAlgorithmInfo())
         std::cout << "----> Running Algorithm Tool: " << this->GetInstanceName() << ", " << this->GetType() << std::endl;
 
     try
     {
         if (!LArPfoHelper::IsTrack(pPfo))
             return;
 
         MatchedSlidingFitMap matchedSlidingFitMap;
         this->BuildSlidingFitMap(pPfo, matchedSlidingFitMap);
 
         if (matchedSlidingFitMap.size() < 2)
             return;
 
         this->GetTrackHits3D(inputTwoDHits, matchedSlidingFitMap, protoHitVector);
     }
     catch (StatusCodeException &)
     {
     }
 }

Member Data Documentation

double lar_content::HitCreationBaseTool::m_chiSquaredCut

protectedinherited

The chi squared cut (accept only values below the cut value)

Definition at line 86 of file HitCreationBaseTool.h.

Referenced by lar_content::ShowerHitsBaseTool::GetShowerHits3D(), lar_content::LongitudinalTrackHitsBaseTool::GetTrackHits3D(), and lar_content::HitCreationBaseTool::ReadSettings().

unsigned int lar_content::TrackHitsBaseTool::m_minViews

protectedinherited

The minimum number of views required for building hits.

Definition at line 58 of file TrackHitsBaseTool.h.

Referenced by lar_content::TrackHitsBaseTool::ReadSettings().

double lar_content::HitCreationBaseTool::m_sigmaX2

protectedinherited

The sigmaX squared value, for calculation of chi2 deltaX term.

Definition at line 85 of file HitCreationBaseTool.h.

Referenced by lar_content::HitCreationBaseTool::ReadSettings(), and lar_content::HitCreationBaseTool::~HitCreationBaseTool().

unsigned int lar_content::TrackHitsBaseTool::m_slidingFitWindow

protectedinherited

The layer window for the sliding linear fits.

Definition at line 59 of file TrackHitsBaseTool.h.

Referenced by lar_content::TrackHitsBaseTool::BuildSlidingFitMap(), and lar_content::TrackHitsBaseTool::ReadSettings().

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

larpandoracontent/v04_08_01/source/larpandoracontent/LArThreeDReco/LArHitCreation/MultiValuedLongitudinalTrackHitsTool.h
larpandoracontent/v04_08_01/source/larpandoracontent/LArThreeDReco/LArHitCreation/MultiValuedLongitudinalTrackHitsTool.cc

Public Types

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Member Typedef Documentation

Member Function Documentation

Member Data Documentation