From 2ce7791be78edc544a2f48fc5be4321f5badcf03 Mon Sep 17 00:00:00 2001
From: rsadek <rita.sadek@cern.ch>
Date: Wed, 27 Oct 2021 19:03:11 +0200
Subject: [PATCH] Fwd Secondary Vertexing: Minor change on search range for
 dichotomy method + cleaned unnecessary cout

---
 .../Vertexing/include/DetectorsVertexing/FwdDCAFitterN.h | 9 ++++-----
 1 file changed, 4 insertions(+), 5 deletions(-)

diff --git a/Detectors/Vertexing/include/DetectorsVertexing/FwdDCAFitterN.h b/Detectors/Vertexing/include/DetectorsVertexing/FwdDCAFitterN.h
index dbb01f04e77a2..c6fad6f71e808 100644
--- a/Detectors/Vertexing/include/DetectorsVertexing/FwdDCAFitterN.h
+++ b/Detectors/Vertexing/include/DetectorsVertexing/FwdDCAFitterN.h
@@ -767,13 +767,13 @@ void FwdDCAFitterN<N, Args...>::findZatXY_mid(int mCurHyp)
 {
   // look into dXY of T0 - T1 between 2 points(0,40cm); the one with the highest dXY is moved to mid
 
-  double startPoint = 0.;
-  double endPoint = 20.; // first disk
+  double startPoint = -40.;
+  double endPoint = 50.;
   double midPoint = 0.5 * (startPoint + endPoint);
 
   double z[2][2] = {{startPoint, endPoint}, {startPoint, endPoint}}; // z for tracks 0/1 on starting poing and endpoint
 
-  double DeltaZ = endPoint - startPoint;
+  double DeltaZ = std::abs(endPoint - startPoint);
 
   double newX[2][2];
   double newY[2][2];
@@ -810,7 +810,7 @@ void FwdDCAFitterN<N, Args...>::findZatXY_mid(int mCurHyp)
     dstXY[1] = (newX[0][1] - newX[1][1]) * (newX[0][1] - newX[1][1]) +
                (newY[0][1] - newY[1][1]) * (newY[0][1] - newY[1][1]);
 
-    DeltaZ = endPoint - startPoint;
+    DeltaZ = std::abs(endPoint - startPoint);
 
     if (DeltaZ < epsilon) {
       finalZ = 0.5 * (startPoint + endPoint);
@@ -1197,7 +1197,6 @@ bool FwdDCAFitterN<N, Args...>::minimizeChi2NoErr()
       return false;
     }
     VecND dz = mD2Chi2Dz2 * mDChi2Dz;
-    std::cout << "the dz: " << dz << std::endl;
 
     if (!FwdcorrectTracks(dz)) {
       return false;