Merge pull request #39873 from swagata87/gsfCrashFixT0

[TsosGaussianStateConversions] Require GSF states to be positive definite [12_6_X]
cms-sw · Oct 31, 2022 · d624d70 · d624d70
2 parents 5932f9b + 57a8f0a
commit d624d70
Showing 1 changed file with 9 additions and 6 deletions.
diff --git a/TrackingTools/GsfTracking/src/TsosGaussianStateConversions.cc b/TrackingTools/GsfTracking/src/TsosGaussianStateConversions.cc
@@ -40,12 +40,15 @@ namespace GaussianStateConversions {
     std::vector<TrajectoryStateOnSurface> components;
     components.reserve(singleStates.size());
     for (auto const& ic : singleStates) {
-      components.emplace_back((*ic).weight(),
-                              LocalTrajectoryParameters((*ic).mean(), pzSign, charged),
-                              LocalTrajectoryError((*ic).covariance()),
-                              surface,
-                              field,
-                              side);
+      //require states to be positive-definite
+      if (double det = 0; (*ic).covariance().Det2(det) && det > 0) {
+        components.emplace_back((*ic).weight(),
+                                LocalTrajectoryParameters((*ic).mean(), pzSign, charged),
+                                LocalTrajectoryError((*ic).covariance()),
+                                surface,
+                                field,
+                                side);
+      }
     }
     return TrajectoryStateOnSurface((BasicTrajectoryState*)new BasicMultiTrajectoryState(components));
   }