Use correct half-length calculation.

gputnam · gputnam · commit 8d355d05c1fe · 2025-02-03T10:44:08.000-06:00
diff --git a/sbncode/Cluster3D/SnippetHit3DBuilderSBN_tool.cc b/sbncode/Cluster3D/SnippetHit3DBuilderSBN_tool.cc
@@ -1714,7 +1714,7 @@ geo::WireID SnippetHit3DBuilderSBN::NearestWireID(const Eigen::Vector3f& positio
         mf::LogWarning("Cluster3D") << "Exception caught finding nearest wire, position - " << exc.what() << std::endl;
 
         // Assume extremum for wire number depending on z coordinate
-        if (position[2] < m_geometry->TPC().HalfLength()) wire = 0;
+        if (position[2] < m_geometry->TPC().ActiveLength()*0.5) wire = 0;
         else                                              wire = m_wireReadout->Nwires(wireIDIn.asPlaneID()) - 1;
     }
 

Original file line number	Diff line number	Diff line change
`@@ -1714,7 +1714,7 @@ geo::WireID SnippetHit3DBuilderSBN::NearestWireID(const Eigen::Vector3f& positio`
`1714`	`1714`	`mf::LogWarning("Cluster3D") << "Exception caught finding nearest wire, position - " << exc.what() << std::endl;`
`1715`	`1715`
`1716`	`1716`	`// Assume extremum for wire number depending on z coordinate`
`1717`		`- if (position[2] < m_geometry->TPC().HalfLength()) wire = 0;`
	`1717`	`+ if (position[2] < m_geometry->TPC().ActiveLength()*0.5) wire = 0;`
`1718`	`1718`	`else wire = m_wireReadout->Nwires(wireIDIn.asPlaneID()) - 1;`
`1719`	`1719`	`}`
`1720`	`1720`