9 #include "Pandora/AlgorithmHeaders.h" 22 OvershootTracksTool::OvershootTracksTool() :
25 m_maxVertexXSeparation(2.
f),
26 m_cosThetaCutForKinkSearch(0.94
f)
44 const unsigned int nMatchedSamplingPoints1((*iIter1)->GetOverlapResult().GetNMatchedSamplingPoints());
45 const unsigned int nMatchedSamplingPoints2((*iIter2)->GetOverlapResult().GetNMatchedSamplingPoints());
49 Particle particle(*(*iIterA), *(*iIterB));
64 const bool isA1LowestInX(this->
IsALowestInX(pointingClusterA1, pointingClusterB1));
65 const bool isA2LowestInX(this->
IsALowestInX(pointingClusterA2, pointingClusterB2));
80 const bool vertex1AIsLowX(vertexA1.GetPosition().GetX() < vertexB1.GetPosition().GetX());
85 const bool vertex2AIsLowX(vertexA2.GetPosition().GetX() < vertexB2.GetPosition().GetX());
97 modificationList.push_back(modification);
99 catch (StatusCodeException &statusCodeException)
101 if (STATUS_CODE_FAILURE == statusCodeException.GetStatusCode())
102 throw statusCodeException;
114 float longitudinalAB(-std::numeric_limits<float>::max()), transverseAB(std::numeric_limits<float>::max());
117 float longitudinalBA(-std::numeric_limits<float>::max()), transverseBA(std::numeric_limits<float>::max());
131 bool splitAtElementA(
false), splitAtElementB(
false);
135 splitAtElementA =
true;
139 splitAtElementB =
true;
142 if (!splitAtElementA && !splitAtElementB)
143 throw StatusCodeException(STATUS_CODE_NOT_FOUND);
148 CartesianVector splitPosition(0.
f, 0.
f, 0.
f);
149 float chiSquared(std::numeric_limits<float>::max());
180 if ((STATUS_CODE_SUCCESS != lowXFitResult1.GetGlobalFitPositionAtX(minusX, minus1)) ||
181 (STATUS_CODE_SUCCESS != highXFitResult1.GetGlobalFitPositionAtX(plusX, plus1)) ||
182 (STATUS_CODE_SUCCESS != lowXFitResult2.GetGlobalFitPositionAtX(minusX, minus2)) ||
183 (STATUS_CODE_SUCCESS != highXFitResult2.GetGlobalFitPositionAtX(plusX, plus2)) ||
184 (STATUS_CODE_SUCCESS != fitResultCommon3.GetGlobalFitPositionAtX(minusX, minus3)) ||
185 (STATUS_CODE_SUCCESS != fitResultCommon3.GetGlobalFitPositionAtX(plusX, plus3)))
195 CartesianVector minus(0.
f, 0.
f, 0.
f), split(0.
f, 0.
f, 0.
f), plus(0.
f, 0.
f, 0.
f);
196 float chi2Minus(std::numeric_limits<float>::max()), chi2Split(std::numeric_limits<float>::max()),
197 chi2Plus(std::numeric_limits<float>::max());
203 const CartesianVector minusToSplit((split - minus).GetUnitVector());
204 const CartesianVector splitToPlus((plus - split).GetUnitVector());
205 const float dotProduct(minusToSplit.GetDotProduct(splitToPlus));
210 catch (StatusCodeException &)
221 m_splitPosition(0.
f, 0.
f, 0.
f),
222 m_splitPosition1(0.
f, 0.
f, 0.
f),
223 m_splitPosition2(0.
f, 0.
f, 0.
f)
225 const HitType commonView((elementA.GetClusterU() == elementB.GetClusterU()) ? TPC_VIEW_U
226 : (elementA.GetClusterV() == elementB.GetClusterV()) ? TPC_VIEW_V
227 : (elementA.GetClusterW() == elementB.GetClusterW()) ? TPC_VIEW_W
230 if (HIT_CUSTOM == commonView)
231 throw StatusCodeException(STATUS_CODE_FAILURE);
234 : (TPC_VIEW_V == commonView) ? elementA.GetClusterV()
235 : elementA.GetClusterW();
236 m_pClusterA1 = (TPC_VIEW_U == commonView) ? elementA.GetClusterV() : elementA.GetClusterU();
237 m_pClusterA2 = (TPC_VIEW_U == commonView) ? elementA.GetClusterW()
238 : (TPC_VIEW_V == commonView) ? elementA.GetClusterW()
239 : elementA.GetClusterV();
240 m_pClusterB1 = (TPC_VIEW_U == commonView) ? elementB.GetClusterV() : elementB.GetClusterU();
241 m_pClusterB2 = (TPC_VIEW_U == commonView) ? elementB.GetClusterW()
242 : (TPC_VIEW_V == commonView) ? elementB.GetClusterW()
243 : elementB.GetClusterV();
246 throw StatusCodeException(STATUS_CODE_FAILURE);
254 PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle,
"SplitMode",
m_splitMode));
256 PANDORA_RETURN_RESULT_IF_AND_IF(
257 STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle,
"MaxVertexXSeparation",
m_maxVertexXSeparation));
259 PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
pandora::ClusterList m_affectedClusters
The list of affected clusters.
std::vector< Modification > ModificationList
float GetXSamplingPoint(const pandora::CartesianVector &splitPosition1, const bool isForwardInX, const TwoDSlidingFitResult &fitResult1, const TwoDSlidingFitResult &fitResult2, const TwoDSlidingFitResult &fitResult3) const
Get a sampling point in x that is common to sliding linear fit objects in three views.
pandora::CartesianVector m_splitPosition2
The candidate split position in view 2.
bool IsThreeDKink(ThreeViewTransverseTracksAlgorithm *const pAlgorithm, const Particle &particle, const bool isA1LowestInX, const bool isA2LowestInX) const
Whether the provided particle is consistent with being a kink, when examined in three dimensions at t...
std::vector< TensorType::ElementList::const_iterator > IteratorList
void GetIteratorListModifications(ThreeViewTransverseTracksAlgorithm *const pAlgorithm, const IteratorList &iteratorList, ModificationList &modificationList) const
Get modification objects for a specific elements of the tensor, identifying required splits and merge...
bool m_splitMode
Whether to run in cluster splitting mode, as opposed to cluster merging mode.
const TwoDSlidingFitResult & GetCachedSlidingFitResult(const pandora::Cluster *const pCluster) const
Get a sliding fit result from the algorithm cache.
static void GetImpactParameters(const LArPointingCluster::Vertex &pointingVertex, const LArPointingCluster::Vertex &targetVertex, float &longitudinal, float &transverse)
Calculate impact parameters between a pair of pointing vertices.
const pandora::Cluster * m_pCommonCluster
Address of the common cluster.
void SetSplitPosition(const LArPointingCluster::Vertex &vertexA1, const LArPointingCluster::Vertex &vertexA2, const LArPointingCluster::Vertex &vertexB1, const LArPointingCluster::Vertex &vertexB2, Particle &particle) const
Set split position for a provided particle.
LArPointingCluster class.
static pandora::HitType GetClusterHitType(const pandora::Cluster *const pCluster)
Get the hit type associated with a two dimensional cluster.
ClusterMergeMap m_clusterMergeMap
The cluster merge map.
pandora::StatusCode ReadSettings(const pandora::TiXmlHandle xmlHandle)
const pandora::Cluster * m_pClusterA1
Address of cluster in element A, view 1.
Header file for the geometry helper class.
float m_maxVertexXSeparation
The max separation between accompanying clusters vertex x positions to make split.
Header file for the cluster helper class.
pandora::StatusCode ReadSettings(const pandora::TiXmlHandle xmlHandle)
static float MergeTwoPositions(const pandora::Pandora &pandora, const pandora::HitType view1, const pandora::HitType view2, const float position1, const float position2)
Merge two views (U,V) to give a third view (Z).
static void MergeThreePositions3D(const pandora::Pandora &pandora, const pandora::HitType view1, const pandora::HitType view2, const pandora::HitType view3, const pandora::CartesianVector &position1, const pandora::CartesianVector &position2, const pandora::CartesianVector &position3, pandora::CartesianVector &position3D, float &chiSquared)
Merge 2D positions from three views to give unified 3D position.
const pandora::Cluster * m_pClusterB2
Address of cluster in element B, view 2.
static bool IsALowestInX(const LArPointingCluster &pointingClusterA, const LArPointingCluster &pointingClusterB)
Whether pointing cluster labelled A extends to lowest x positions (as opposed to that labelled B) ...
float m_cosThetaCutForKinkSearch
The cos theta cut used for the kink search in three dimensions.
pandora::CartesianVector m_splitPosition
The candidate split position for the common cluster.
Particle(const TensorType::Element &elementA, const TensorType::Element &elementB)
Constructor.
bool PassesVertexCuts(const LArPointingCluster::Vertex &vertexA, const LArPointingCluster::Vertex &vertexB) const
Whether a pair of vertices pass longitudinal projection cuts.
pandora::CartesianVector m_splitPosition1
The candidate split position in view 1.
ThreeDKinkBaseTool class.
static void GetClosestVerticesInX(const LArPointingCluster &pointingClusterI, const LArPointingCluster &pointingClusterJ, LArPointingCluster::Vertex &closestVertexI, LArPointingCluster::Vertex &closestVertexJ)
Given a pair of pointing clusters, find the pair of vertices with smallest x-separation.
bool m_majorityRulesMode
Whether to run in majority rules mode (always split overshoots, always merge undershoots) ...
float m_minLongitudinalImpactParameter
The min longitudinal impact parameter for connecting accompanying clusters.
const pandora::Cluster * m_pClusterA2
Address of cluster in element A, view 2.
ThreeViewTransverseTracksAlgorithm class.
const pandora::Cluster * m_pClusterB1
Address of cluster in element B, view 1.
const pandora::CartesianVector & GetPosition() const
Get the vertex position.
TwoDSlidingFitResult class.
SplitPositionMap m_splitPositionMap
The split position map.