add check on PATLeptonTimeLifeInfoProducer closestState valid state

cms-sw · Apr 29, 2024 · 36effa5 · 36effa5
1 parent 43944b8
commit 36effa5
Showing 1 changed file with 9 additions and 0 deletions.
diff --git a/PhysicsTools/PatAlgos/plugins/PATLeptonTimeLifeInfoProducer.cc b/PhysicsTools/PatAlgos/plugins/PATLeptonTimeLifeInfoProducer.cc
@@ -12,6 +12,7 @@
 #include "FWCore/ParameterSet/interface/ParameterSet.h"
 #include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
 #include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
+#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
 #include "DataFormats/Common/interface/ValueMap.h"
 #include "DataFormats/PatCandidates/interface/Electron.h"
@@ -169,6 +170,14 @@ void PATLeptonTimeLifeInfoProducer<T>::produceAndFillIPInfo(const T& lepton,
     AnalyticalImpactPointExtrapolator extrapolator(transTrack.field());
     TrajectoryStateOnSurface closestState =
         extrapolator.extrapolate(transTrack.impactPointState(), RecoVertex::convertPos(pv.position()));
+    if (!closestState.isValid()) {
+      edm::LogError("PATLeptonTimeLifeInfoProducer")
+          << "closestState not valid! From:\n"
+          << "transTrack.impactPointState():\n"
+          << transTrack.impactPointState() << "RecoVertex::convertPos(pv.position()):\n"
+          << RecoVertex::convertPos(pv.position());
+      return;
+    }
     GlobalPoint pca = closestState.globalPosition();
     GlobalError pca_cov = closestState.cartesianError().position();
     GlobalVector ip_vec = GlobalVector(pca.x() - pv.x(), pca.y() - pv.y(), pca.z() - pv.z());