Skip to content

Commit 325c13e

Browse files
committed
Protection against acos argument > 1
1 parent 71e1e50 commit 325c13e

File tree

1 file changed

+5
-2
lines changed
  • Detectors/ITSMFT/MFT/tracking/include/MFTTracking

1 file changed

+5
-2
lines changed

Detectors/ITSMFT/MFT/tracking/include/MFTTracking/Tracker.h

Lines changed: 5 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -212,8 +212,11 @@ inline const Float_t Tracker<T>::getCellDeviation(const Cell& cell1, const Cell&
212212
Float_t cell2mod = std::sqrt(cell2dx * cell2dx + cell2dy * cell2dy + cell2dz * cell2dz);
213213

214214
Float_t cosAngle = (cell1dx * cell2dx + cell1dy * cell2dy + cell1dz * cell2dz) / (cell1mod * cell2mod);
215-
216-
return std::acos(cosAngle);
215+
if (std::abs(cosAngle) < o2::constants::math::Almost1) {
216+
return std::acos(cosAngle);
217+
} else {
218+
return cosAngle > 0 ? 0. : -o2::constants::math::PI;
219+
}
217220
}
218221

219222
//_________________________________________________________________________________________________

0 commit comments

Comments
 (0)