trkParamErr upper cut

cms-sw · Jul 27, 2024 · 6f3a1c9 · 6f3a1c9
1 parent b770785
commit 6f3a1c9
Show file tree

Hide file tree

Showing 3 changed files with 6 additions and 3 deletions.
diff --git a/RecoParticleFlow/PFTracking/interface/PFDisplacedVertexCandidateFinder.h b/RecoParticleFlow/PFTracking/interface/PFDisplacedVertexCandidateFinder.h
@@ -47,6 +47,7 @@ class PFDisplacedVertexCandidateFinder {
     pt_min_ = ps_trk.getParameter<double>("pt_min");
     pt_min_prim_ = ps_trk.getParameter<double>("pt_min_prim");
     dxy_ = ps_trk.getParameter<double>("dxy");
+    qoverpError_max_ = ps_trk.getParameter<double>("qoverpError_max");
   }
 
   /// sets debug printout flag
@@ -137,7 +138,8 @@ class PFDisplacedVertexCandidateFinder {
 
   double pt_min_prim_;
   double dxy_;
-
+  double qoverpError_max_;
+
   /// Max number of expected vertexCandidates in the event
   /// Used to allocate the memory and avoid multiple copy
   unsigned vertexCandidatesSize_;

diff --git a/RecoParticleFlow/PFTracking/python/particleFlowDisplacedVertexCandidate_cfi.py b/RecoParticleFlow/PFTracking/python/particleFlowDisplacedVertexCandidate_cfi.py
@@ -39,7 +39,7 @@
     # PFDisplacedVertex timing
         pt_min_prim = cms.double(.8),
         dxy = cms.double(.2),
-
+        qoverpError_max =cms.double(1.0e+7)
     )                                                     
 
 )

diff --git a/RecoParticleFlow/PFTracking/src/PFDisplacedVertexCandidateFinder.cc b/RecoParticleFlow/PFTracking/src/PFDisplacedVertexCandidateFinder.cc
@@ -314,10 +314,11 @@ bool PFDisplacedVertexCandidateFinder::goodPtResolution(const TrackBaseRef& trac
   double dxy = trackref->dxy(pvtx_);
 
   double pt_error = dpt / pt * 100;
+  double qoverpError = trackref->qoverpError();
 
   LogDebug("PFDisplacedVertexCandidateFinder")
       << " PFDisplacedVertexFinder: PFrecTrack->Track Pt= " << pt << " dPt/Pt = " << pt_error << "% nChi2 = " << nChi2;
-  if (nChi2 > nChi2_max_ || pt < pt_min_) {
+  if (nChi2 > nChi2_max_ || pt < pt_min_ || qoverpError > qoverpError_max_) {
     LogDebug("PFDisplacedVertexCandidateFinder") << " PFBlockAlgo: skip badly measured or low pt track"
                                                  << " nChi2_cut = " << 5 << " pt_cut = " << 0.2;
     return false;