#include "DeltaRayShowerHitsTool.h"

Inheritance diagram for lar_content::DeltaRayShowerHitsTool:

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

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

Protected Attributes
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	CreateDeltaRayShowerHits3D (const pandora::CaloHitVector &inputTwoDHits, const pandora::CaloHitVector &parentHits3D, ProtoHitVector &protoHitVector) const
	Create three dimensional hits, using a list of input two dimensional hits and the 3D hits from the parent particle. More...

Detailed Description

DeltaRayShowerHitsTool class.

Definition at line 19 of file DeltaRayShowerHitsTool.h.

Member Typedef Documentation

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::DeltaRayShowerHitsTool::CreateDeltaRayShowerHits3D	(	const pandora::CaloHitVector &	inputTwoDHits,
		const pandora::CaloHitVector &	parentHits3D,
		ProtoHitVector &	protoHitVector
	)		const

private

Create three dimensional hits, using a list of input two dimensional hits and the 3D hits from the parent particle.

Parameters

inputTwoDHits	the vector of input two dimensional hits
parentHits3D	the vector of 3D hits from the parent particle
protoHitVector	to receive the new three dimensional proto hits

Definition at line 54 of file DeltaRayShowerHitsTool.cc.

References f, and lar_content::ThreeDHitCreationAlgorithm::ProtoHit::IsPositionSet().

 {
     for (const CaloHit *const pCaloHit2D : inputTwoDHits)
     {
         try
         {
             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);
 
             bool foundClosestPosition(false);
             float closestDistanceSquared(std::numeric_limits<float>::max());
             CartesianVector closestPosition3D(0.f, 0.f, 0.f);
 
             for (const CaloHit *const pCaloHit3D : parentHits3D)
             {
                 const CartesianVector thisPosition3D(pCaloHit3D->GetPositionVector());
                 const CartesianVector thisPosition2D(LArGeometryHelper::ProjectPosition(this->GetPandora(), thisPosition3D, hitType));
                 const float thisDistanceSquared((pCaloHit2D->GetPositionVector() - thisPosition2D).GetMagnitudeSquared());
 
                 if (thisDistanceSquared < closestDistanceSquared)
                 {
                     foundClosestPosition = true;
                     closestDistanceSquared = thisDistanceSquared;
                     closestPosition3D = thisPosition3D;
                 }
             }
 
             if (!foundClosestPosition)
                 continue;
 
             const CartesianVector position1(LArGeometryHelper::ProjectPosition(this->GetPandora(), closestPosition3D, hitType1));
             const CartesianVector position2(LArGeometryHelper::ProjectPosition(this->GetPandora(), closestPosition3D, hitType2));
 
             ProtoHit protoHit(pCaloHit2D);
             this->GetBestPosition3D(hitType1, hitType2, position1, position2, protoHit);
 
             if (protoHit.IsPositionSet())
                 protoHitVector.push_back(protoHit);
         }
         catch (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

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

protectedvirtualinherited

Reimplemented in lar_content::TrackHitsBaseTool, lar_content::ShowerHitsBaseTool, lar_content::LongitudinalTrackHitsBaseTool, and lar_content::ThreeViewShowerHitsTool.

Definition at line 141 of file HitCreationBaseTool.cc.

References lar_content::HitCreationBaseTool::m_chiSquaredCut, and lar_content::HitCreationBaseTool::m_sigmaX2.

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

 {
     double sigmaX(std::sqrt(m_sigmaX2));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "SigmaX", sigmaX));
 
     m_sigmaX2 = sigmaX * sigmaX;
 
     if (m_sigmaX2 < std::numeric_limits<double>::epsilon())
     {
         std::cout << "HitCreationBaseTool - Invalid parameter, SigmaX: " << sigmaX << std::endl;
         return STATUS_CODE_INVALID_PARAMETER;
     }
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "ChiSquaredCut", m_chiSquaredCut));
 
     return STATUS_CODE_SUCCESS;
 }

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

virtual

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 23 of file DeltaRayShowerHitsTool.cc.

 {
     if (PandoraContentApi::GetSettings(*pAlgorithm)->ShouldDisplayAlgorithmInfo())
         std::cout << "----> Running Algorithm Tool: " << this->GetInstanceName() << ", " << this->GetType() << std::endl;
 
     try
     {
         if (!LArPfoHelper::IsShower(pPfo) || (1 != pPfo->GetParentPfoList().size()))
             return;
 
         const ParticleFlowObject *const pParentPfo(pPfo->GetParentPfoList().front());
 
         CaloHitList parentHitList3D;
         LArPfoHelper::GetCaloHits(pParentPfo, TPC_3D, parentHitList3D);
 
         if (parentHitList3D.empty())
             return;
 
         CaloHitVector parentHitVector3D(parentHitList3D.begin(), parentHitList3D.end());
         std::sort(parentHitVector3D.begin(), parentHitVector3D.end(), LArClusterHelper::SortHitsByPosition);
 
         this->CreateDeltaRayShowerHits3D(inputTwoDHits, parentHitVector3D, 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().

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

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

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

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Member Typedef Documentation

Member Function Documentation

Member Data Documentation