CMS 3D CMS Logo

List of all members | Public Types | Public Member Functions | Private Member Functions | Private Attributes | Friends
PFDisplacedVertexFinder Class Reference

#include <PFDisplacedVertexFinder.h>

Public Types

enum  FitterType { F_NOTDEFINED, F_DONOTREFIT, F_KALMAN, F_ADAPTIVE }
 Fitter Type. More...
 
typedef reco::PFDisplacedVertexCollection::iterator IDV
 
typedef reco::PFDisplacedVertexSeedCollection::iterator IDVS
 -----— Useful Types -----— /// More...
 
typedef std::pair< PFTrackHitInfo, PFTrackHitInfoPFTrackHitFullInfo
 
typedef std::pair< unsigned int, unsigned int > PFTrackHitInfo
 

Public Member Functions

const std::unique_ptr< reco::PFDisplacedVertexCollection > & displacedVertices () const
 
void findDisplacedVertices ()
 -----— Main function which find vertices -----— /// More...
 
 PFDisplacedVertexFinder ()
 
void setAVFParameters (const edm::ParameterSet &ps)
 
void setDebug (bool debug)
 Sets debug printout flag. More...
 
void setEdmParameters (const MagneticField *magField, edm::ESHandle< GlobalTrackingGeometry > globTkGeomHandle, const TrackerTopology *tkerTopo, const TrackerGeometry *tkerGeom)
 Sets parameters for track extrapolation and hits study. More...
 
void setInput (const edm::Handle< reco::PFDisplacedVertexCandidateCollection > &)
 Set input collections of tracks. More...
 
void setParameters (double transvSize, double longSize, double primaryVertexCut, double tobCut, double tecCut, double minAdaptWeight, bool switchOff2TrackVertex)
 -----— Set different algo parameters -----— /// More...
 
void setPrimaryVertex (edm::Handle< reco::VertexCollection > mainVertexHandle, edm::Handle< reco::BeamSpot > beamSpotHandle)
 
void setTracksSelector (const edm::ParameterSet &ps)
 
void setVertexIdentifier (const edm::ParameterSet &ps)
 
std::unique_ptr< reco::PFDisplacedVertexCollectiontransferDisplacedVertices ()
 
 ~PFDisplacedVertexFinder ()
 

Private Member Functions

unsigned commonTracks (const reco::PFDisplacedVertex &, const reco::PFDisplacedVertex &) const
 
void findSeedsFromCandidate (const reco::PFDisplacedVertexCandidate &, reco::PFDisplacedVertexSeedCollection &)
 -----— Different steps of the finder algorithm -----— /// More...
 
bool fitVertexFromSeed (const reco::PFDisplacedVertexSeed &, reco::PFDisplacedVertex &)
 Fit one by one the vertex points with associated tracks to get displaced vertices. More...
 
std::pair< float, float > getTransvLongDiff (const GlobalPoint &, const GlobalPoint &) const
 
reco::PFDisplacedVertex::VertexTrackType getVertexTrackType (PFTrackHitFullInfo &) const
 
bool isCloseTo (const reco::PFDisplacedVertexSeed &, const reco::PFDisplacedVertexSeed &) const
 -----— Tools -----— /// More...
 
void mergeSeeds (reco::PFDisplacedVertexSeedCollection &, std::vector< bool > &bLocked)
 Sometimes two vertex candidates can be quite close and coming from the same vertex. More...
 
bool rejectAndLabelVertex (reco::PFDisplacedVertex &dv)
 
void selectAndLabelVertices (reco::PFDisplacedVertexCollection &, std::vector< bool > &)
 Remove potentially fakes displaced vertices. More...
 

Private Attributes

bool debug_
 If true, debug printouts activated. More...
 
reco::PFDisplacedVertexCandidateCollection const * displacedVertexCandidates_
 -----— Members -----— /// More...
 
std::unique_ptr< reco::PFDisplacedVertexCollectiondisplacedVertices_
 
edm::ESHandle< GlobalTrackingGeometryglobTkGeomHandle_
 Tracker geometry for discerning hit positions. More...
 
PFDisplacedVertexHelper helper_
 
PFCheckHitPattern hitPattern_
 
float longSize_
 
const MagneticFieldmagField_
 to be able to extrapolate tracks f More...
 
double minAdaptWeight_
 
double primaryVertexCut_
 
double ratio_
 
double sigmacut_
 Adaptive Vertex Fitter parameters. More...
 
bool switchOff2TrackVertex_
 
double t_ini_
 
double tecCut_
 
const TrackerGeometrytkerGeom_
 
const TrackerTopologytkerTopo_
 doc? More...
 
double tobCut_
 
float transvSize_
 -----— Parameters -----— /// More...
 

Friends

std::ostream & operator<< (std::ostream &, const PFDisplacedVertexFinder &)
 

Detailed Description

Definition at line 31 of file PFDisplacedVertexFinder.h.

Member Typedef Documentation

◆ IDV

typedef reco::PFDisplacedVertexCollection::iterator PFDisplacedVertexFinder::IDV

Definition at line 40 of file PFDisplacedVertexFinder.h.

◆ IDVS

typedef reco::PFDisplacedVertexSeedCollection::iterator PFDisplacedVertexFinder::IDVS

-----— Useful Types -----— ///

Definition at line 39 of file PFDisplacedVertexFinder.h.

◆ PFTrackHitFullInfo

Definition at line 43 of file PFDisplacedVertexFinder.h.

◆ PFTrackHitInfo

typedef std::pair<unsigned int, unsigned int> PFDisplacedVertexFinder::PFTrackHitInfo

Definition at line 42 of file PFDisplacedVertexFinder.h.

Member Enumeration Documentation

◆ FitterType

Fitter Type.

Enumerator
F_NOTDEFINED 
F_DONOTREFIT 
F_KALMAN 
F_ADAPTIVE 

Definition at line 46 of file PFDisplacedVertexFinder.h.

Constructor & Destructor Documentation

◆ PFDisplacedVertexFinder()

PFDisplacedVertexFinder::PFDisplacedVertexFinder ( )

Definition at line 30 of file PFDisplacedVertexFinder.cc.

31  : displacedVertexCandidates_(nullptr),
33  transvSize_(0.0),
34  longSize_(0.0),
35  primaryVertexCut_(0.0),
36  tobCut_(0.0),
37  tecCut_(0.0),
38  minAdaptWeight_(2.0),
39  debug_(false) {}

◆ ~PFDisplacedVertexFinder()

PFDisplacedVertexFinder::~PFDisplacedVertexFinder ( )

Definition at line 41 of file PFDisplacedVertexFinder.cc.

41 {}

Member Function Documentation

◆ commonTracks()

unsigned PFDisplacedVertexFinder::commonTracks ( const reco::PFDisplacedVertex v1,
const reco::PFDisplacedVertex v2 
) const
private

Definition at line 661 of file PFDisplacedVertexFinder.cc.

661  {
662  vector<Track> vt1 = v1.refittedTracks();
663  vector<Track> vt2 = v2.refittedTracks();
664 
665  unsigned commonTracks = 0;
666 
667  for (unsigned il1 = 0; il1 < vt1.size(); il1++) {
668  unsigned il1_idx = v1.originalTrack(vt1[il1]).key();
669 
670  for (unsigned il2 = 0; il2 < vt2.size(); il2++)
671  if (il1_idx == v2.originalTrack(vt2[il2]).key()) {
672  commonTracks++;
673  break;
674  }
675  }
676 
677  return commonTracks;
678 }

References edm::RefToBase< T >::key(), reco::Vertex::originalTrack(), and reco::Vertex::refittedTracks().

Referenced by selectAndLabelVertices().

◆ displacedVertices()

const std::unique_ptr<reco::PFDisplacedVertexCollection>& PFDisplacedVertexFinder::displacedVertices ( ) const
inline

Definition at line 104 of file PFDisplacedVertexFinder.h.

104  {
106  }

References displacedVertices_, and eostools::move().

◆ findDisplacedVertices()

void PFDisplacedVertexFinder::findDisplacedVertices ( )

-----— Main function which find vertices -----— ///

Definition at line 54 of file PFDisplacedVertexFinder.cc.

54  {
55  if (debug_)
56  cout << "========= Start Find Displaced Vertices =========" << endl;
57 
58  // The vertexCandidates have not been passed to the event
59  // So they need to be cleared is they are not empty
60  if (displacedVertices_.get())
61  displacedVertices_->clear();
62  else
63  displacedVertices_ = std::make_unique<PFDisplacedVertexCollection>();
64 
65  if (displacedVertexCandidates_ == nullptr) {
66  edm::LogInfo("EmptyVertexInput")
67  << "displacedVertexCandidates are not set or the setInput was called with invalid vertex";
68  return;
69  }
70 
71  // Prepare the collections
72  PFDisplacedVertexSeedCollection tempDisplacedVertexSeeds;
73  tempDisplacedVertexSeeds.reserve(4 * displacedVertexCandidates_->size());
74  PFDisplacedVertexCollection tempDisplacedVertices;
75  tempDisplacedVertices.reserve(4 * displacedVertexCandidates_->size());
76 
77  if (debug_)
78  cout << "1) Parsing displacedVertexCandidates into displacedVertexSeeds" << endl;
79 
80  // 1) Parsing displacedVertexCandidates into displacedVertexSeeds which would
81  // be later used for vertex fitting
82 
83  int i = -1;
84 
85  for (auto const& idvc : *displacedVertexCandidates_) {
86  i++;
87  if (debug_) {
88  cout << "Analyse Vertex Candidate " << i << endl;
89  }
90 
91  findSeedsFromCandidate(idvc, tempDisplacedVertexSeeds);
92  }
93 
94  if (debug_)
95  cout << "2) Merging Vertex Seeds" << endl;
96 
97  // 2) Some displacedVertexSeeds coming from different displacedVertexCandidates
98  // may be closed enough to be merged together. bLocked is an array which keeps the
99  // information about the seeds which are desabled.
100  vector<bool> bLockedSeeds;
101  bLockedSeeds.resize(tempDisplacedVertexSeeds.size());
102  mergeSeeds(tempDisplacedVertexSeeds, bLockedSeeds);
103 
104  if (debug_)
105  cout << "3) Fitting Vertices From Seeds" << endl;
106 
107  // 3) Fit displacedVertices from displacedVertexSeeds
108  for (unsigned idv = 0; idv < tempDisplacedVertexSeeds.size(); idv++) {
109  if (!tempDisplacedVertexSeeds[idv].isEmpty() && !bLockedSeeds[idv]) {
110  PFDisplacedVertex displacedVertex;
111  bLockedSeeds[idv] = fitVertexFromSeed(tempDisplacedVertexSeeds[idv], displacedVertex);
112  if (!bLockedSeeds[idv])
113  tempDisplacedVertices.emplace_back(displacedVertex);
114  }
115  }
116 
117  if (debug_)
118  cout << "4) Rejecting Bad Vertices and label them" << endl;
119 
120  // 4) Reject displaced vertices which may be considered as fakes
121  vector<bool> bLocked;
122  bLocked.resize(tempDisplacedVertices.size());
123  selectAndLabelVertices(tempDisplacedVertices, bLocked);
124 
125  if (debug_)
126  cout << "5) Fill the Displaced Vertices" << endl;
127 
128  // 5) Fill the displacedVertex_ which would be transfered to the producer
129  displacedVertices_->reserve(tempDisplacedVertices.size());
130 
131  for (unsigned idv = 0; idv < tempDisplacedVertices.size(); idv++)
132  if (!bLocked[idv])
133  displacedVertices_->push_back(tempDisplacedVertices[idv]);
134 
135  if (debug_)
136  cout << "========= End Find Displaced Vertices =========" << endl;
137 }

References gather_cfg::cout, debug_, displacedVertexCandidates_, displacedVertices_, findSeedsFromCandidate(), fitVertexFromSeed(), mps_fire::i, mergeSeeds(), and selectAndLabelVertices().

◆ findSeedsFromCandidate()

void PFDisplacedVertexFinder::findSeedsFromCandidate ( const reco::PFDisplacedVertexCandidate vertexCandidate,
reco::PFDisplacedVertexSeedCollection tempDisplacedVertexSeeds 
)
private

-----— Different steps of the finder algorithm -----— ///

Analyse a vertex candidate and select potential vertex point(s)

Definition at line 141 of file PFDisplacedVertexFinder.cc.

142  {
143  const PFDisplacedVertexCandidate::DistMap r2Map = vertexCandidate.r2Map();
144  bool bNeedNewCandidate = false;
145 
146  tempDisplacedVertexSeeds.push_back(PFDisplacedVertexSeed());
147 
148  IDVS idvc_current;
149 
150  for (PFDisplacedVertexCandidate::DistMap::const_iterator imap = r2Map.begin(); imap != r2Map.end(); imap++) {
151  unsigned ie1 = (*imap).second.first;
152  unsigned ie2 = (*imap).second.second;
153 
154  if (debug_)
155  cout << "ie1 = " << ie1 << " ie2 = " << ie2 << " radius = " << sqrt((*imap).first) << endl;
156 
157  GlobalPoint dcaPoint = vertexCandidate.dcaPoint(ie1, ie2);
158  if (fabs(dcaPoint.x()) > 1e9)
159  continue;
160 
161  bNeedNewCandidate = true;
162  for (idvc_current = tempDisplacedVertexSeeds.begin(); idvc_current != tempDisplacedVertexSeeds.end();
163  idvc_current++) {
164  if ((*idvc_current).isEmpty()) {
165  bNeedNewCandidate = false;
166  break;
167  }
168  const GlobalPoint vertexPoint = (*idvc_current).seedPoint();
169  std::pair<float, float> diffs = getTransvLongDiff(vertexPoint, dcaPoint);
170  if (diffs.second > longSize_)
171  continue;
172  if (diffs.first > transvSize_)
173  continue;
174  bNeedNewCandidate = false;
175  break;
176  }
177  if (bNeedNewCandidate) {
178  if (debug_)
179  cout << "create new displaced vertex" << endl;
180  tempDisplacedVertexSeeds.push_back(PFDisplacedVertexSeed());
181  idvc_current = tempDisplacedVertexSeeds.end();
182  idvc_current--;
183  }
184 
185  (*idvc_current).updateSeedPoint(dcaPoint, vertexCandidate.tref(ie1), vertexCandidate.tref(ie2));
186  }
187 }

References gather_cfg::cout, reco::PFDisplacedVertexCandidate::dcaPoint(), debug_, getTransvLongDiff(), longSize_, reco::PFDisplacedVertexCandidate::r2Map(), mathSSE::sqrt(), transvSize_, reco::PFDisplacedVertexCandidate::tref(), and PV3DBase< T, PVType, FrameType >::x().

Referenced by findDisplacedVertices().

◆ fitVertexFromSeed()

bool PFDisplacedVertexFinder::fitVertexFromSeed ( const reco::PFDisplacedVertexSeed displacedVertexSeed,
reco::PFDisplacedVertex displacedVertex 
)
private

Fit one by one the vertex points with associated tracks to get displaced vertices.

This prefit procedure allow to reduce the Warning rate from Adaptive Vertex fitter It reject also many fake tracks

Definition at line 209 of file PFDisplacedVertexFinder.cc.

210  {
211  if (debug_)
212  cout << "== Start vertexing procedure ==" << endl;
213 
214  // ---- Prepare transient track list ----
215 
216  auto const& tracksToFit = displacedVertexSeed.elements();
217  const GlobalPoint& seedPoint = displacedVertexSeed.seedPoint();
218 
219  vector<TransientTrack> transTracks;
220  vector<TransientTrack> transTracksRaw;
221  vector<TrackBaseRef> transTracksRef;
222 
223  transTracks.reserve(tracksToFit.size());
224  transTracksRaw.reserve(tracksToFit.size());
225  transTracksRef.reserve(tracksToFit.size());
226 
227  TransientVertex theVertexAdaptiveRaw;
228  TransientVertex theRecoVertex;
229 
230  // ---- 1) Clean for potentially poor seeds ------- //
231  // --------------------------------------------- //
232 
233  if (tracksToFit.size() < 2) {
234  if (debug_)
235  cout << "Only one to Fit Track" << endl;
236  return true;
237  }
238 
239  double rho = sqrt(seedPoint.x() * seedPoint.x() + seedPoint.y() * seedPoint.y());
240  double z = seedPoint.z();
241 
242  if (rho > tobCut_ || fabs(z) > tecCut_) {
243  if (debug_)
244  cout << "Seed Point out of the tracker rho = " << rho << " z = " << z << " nTracks = " << tracksToFit.size()
245  << endl;
246  return true;
247  }
248 
249  if (debug_)
250  displacedVertexSeed.Dump();
251 
252  int nStep45 = 0;
253  int nNotIterative = 0;
254 
255  // Fill vectors of TransientTracks and TrackRefs after applying preselection cuts.
256  for (auto const& ie : tracksToFit) {
257  TransientTrack tmpTk(*(ie.get()), magField_, globTkGeomHandle_);
258  transTracksRaw.emplace_back(tmpTk);
259  bool nonIt = PFTrackAlgoTools::nonIterative((ie)->algo());
260  bool step45 = PFTrackAlgoTools::step45((ie)->algo());
261  bool highQ = PFTrackAlgoTools::highQuality((ie)->algo());
262  if (step45)
263  nStep45++;
264  else if (nonIt)
265  nNotIterative++;
266  else if (!highQ) {
267  nNotIterative++;
268  nStep45++;
269  }
270  }
271 
272  if (rho > 25 && nStep45 + nNotIterative < 1) {
273  if (debug_)
274  cout << "Seed point at rho > 25 cm but no step 4-5 tracks" << endl;
275  return true;
276  }
277 
278  // ----------------------------------------------- //
279  // ---- PRESELECT GOOD TRACKS FOR FINAL VERTEX --- //
280  // ----------------------------------------------- //
281 
282  // 1)If only two track are found do not prefit
283 
284  if (transTracksRaw.size() == 2) {
285  if (debug_)
286  cout << "No raw fit done" << endl;
288  if (debug_)
289  cout << "Due to probably high pile-up conditions 2 track vertices switched off" << endl;
290  return true;
291  }
292  GlobalError globalError;
293 
294  theVertexAdaptiveRaw = TransientVertex(seedPoint, globalError, transTracksRaw, 1.);
295 
296  } else { //branch with transTracksRaw.size of at least 3
297 
298  if (debug_)
299  cout << "Raw fit done." << endl;
300 
306 
307  if (transTracksRaw.size() == 3) {
308  theVertexAdaptiveRaw = theAdaptiveFitterRaw.vertex(transTracksRaw, seedPoint);
309 
310  } else if (transTracksRaw.size() < 1000) {
313 
314  if (debug_)
315  cout << "First test with KFT" << endl;
316 
317  KalmanVertexFitter theKalmanFitter(true);
318  theVertexAdaptiveRaw = theKalmanFitter.vertex(transTracksRaw, seedPoint);
319 
320  if (!theVertexAdaptiveRaw.isValid() || theVertexAdaptiveRaw.totalChiSquared() < 0.) {
321  if (debug_)
322  cout << "Prefit failed : valid? " << theVertexAdaptiveRaw.isValid()
323  << " totalChi2 = " << theVertexAdaptiveRaw.totalChiSquared() << endl;
324  return true;
325  }
326 
327  if (debug_)
328  cout << "We use KFT instead of seed point to set up a point for AVF "
329  << " x = " << theVertexAdaptiveRaw.position().x() << " y = " << theVertexAdaptiveRaw.position().y()
330  << " z = " << theVertexAdaptiveRaw.position().z() << endl;
331 
332  // To save time: reject the Displaced vertex if it is too close to the beam pipe.
333  // Frequently it is very big vertices, with some dosens of tracks
334 
335  Vertex vtx = theVertexAdaptiveRaw;
336  rho = vtx.position().rho();
337 
338  // cout << "primary vertex cut = " << primaryVertexCut_ << endl;
339 
340  if (rho < primaryVertexCut_ || rho > 100) {
341  if (debug_)
342  cout << "KFT Vertex geometrically rejected with tracks #rho = " << rho << endl;
343  return true;
344  }
345 
346  // cout << "primary vertex cut = " << primaryVertexCut_ << " rho = " << rho << endl;
347 
348  theVertexAdaptiveRaw = theAdaptiveFitterRaw.vertex(transTracksRaw, theVertexAdaptiveRaw.position());
349 
350  } else {
351  edm::LogWarning("TooManyPFDVCandidates") << "gave up vertex reco for " << transTracksRaw.size() << " tracks";
352  }
353 
354  if (!theVertexAdaptiveRaw.isValid() || theVertexAdaptiveRaw.totalChiSquared() < 0.) {
355  if (debug_)
356  cout << "Fit failed : valid? " << theVertexAdaptiveRaw.isValid()
357  << " totalChi2 = " << theVertexAdaptiveRaw.totalChiSquared() << endl;
358  return true;
359  }
360 
361  // To save time: reject the Displaced vertex if it is too close to the beam pipe.
362  // Frequently it is very big vertices, with some dosens of tracks
363 
364  Vertex vtx = theVertexAdaptiveRaw;
365  rho = vtx.position().rho();
366 
367  if (rho < primaryVertexCut_) {
368  if (debug_)
369  cout << "Vertex "
370  << " geometrically rejected with " << transTracksRaw.size() << " tracks #rho = " << rho << endl;
371  return true;
372  }
373  }
374 
375  // ---- Remove tracks with small weight or
376  // big first (last) hit_to_vertex distance
377  // and then refit ---- //
378 
379  for (unsigned i = 0; i < transTracksRaw.size(); i++) {
380  if (debug_)
381  cout << "Raw Weight track " << i << " = " << theVertexAdaptiveRaw.trackWeight(transTracksRaw[i]) << endl;
382 
383  if (theVertexAdaptiveRaw.trackWeight(transTracksRaw[i]) > minAdaptWeight_) {
384  PFTrackHitFullInfo pattern = hitPattern_.analyze(tkerTopo_, tkerGeom_, tracksToFit[i], theVertexAdaptiveRaw);
385 
387 
388  if (vertexTrackType != PFDisplacedVertex::T_NOT_FROM_VERTEX) {
389  bool bGoodTrack = helper_.isTrackSelected(transTracksRaw[i].track(), vertexTrackType);
390 
391  if (bGoodTrack) {
392  transTracks.push_back(transTracksRaw[i]);
393  transTracksRef.push_back(tracksToFit[i]);
394  } else {
395  if (debug_)
396  cout << "Track rejected nChi2 = " << transTracksRaw[i].track().normalizedChi2()
397  << " pt = " << transTracksRaw[i].track().pt()
398  << " dxy (wrt (0,0,0)) = " << transTracksRaw[i].track().dxy()
399  << " nHits = " << transTracksRaw[i].track().numberOfValidHits() << " nOuterHits = "
400  << transTracksRaw[i].track().hitPattern().numberOfLostHits(HitPattern::MISSING_OUTER_HITS) << endl;
401  }
402  } else {
403  if (debug_) {
404  cout << "Remove track because too far away from the vertex:" << endl;
405  }
406  }
407  }
408  }
409 
410  if (debug_)
411  cout << "All Tracks " << transTracksRaw.size() << " with good weight " << transTracks.size() << endl;
412 
413  // ---- Refit ---- //
414  FitterType vtxFitter = F_NOTDEFINED;
415 
416  if (transTracks.size() < 2)
417  return true;
418  else if (transTracks.size() == 2) {
420  if (debug_)
421  cout << "Due to probably high pile-up conditions 2 track vertices switched off" << endl;
422  return true;
423  }
424  vtxFitter = F_KALMAN;
425  } else if (transTracks.size() > 2 && transTracksRaw.size() > transTracks.size())
426  vtxFitter = F_ADAPTIVE;
427  else if (transTracks.size() > 2 && transTracksRaw.size() == transTracks.size())
428  vtxFitter = F_DONOTREFIT;
429  else
430  return true;
431 
432  if (debug_)
433  cout << "Vertex Fitter " << vtxFitter << endl;
434 
435  if (vtxFitter == F_KALMAN) {
436  KalmanVertexFitter theKalmanFitter(true);
437  theRecoVertex = theKalmanFitter.vertex(transTracks, seedPoint);
438 
439  } else if (vtxFitter == F_ADAPTIVE) {
445 
446  theRecoVertex = theAdaptiveFitter.vertex(transTracks, seedPoint);
447 
448  } else if (vtxFitter == F_DONOTREFIT) {
449  theRecoVertex = theVertexAdaptiveRaw;
450  } else {
451  return true;
452  }
453 
454  // ---- Check if the fitted vertex is valid ---- //
455 
456  if (!theRecoVertex.isValid() || theRecoVertex.totalChiSquared() < 0.) {
457  if (debug_)
458  cout << "Refit failed : valid? " << theRecoVertex.isValid() << " totalChi2 = " << theRecoVertex.totalChiSquared()
459  << endl;
460  return true;
461  }
462 
463  // ---- Create vertex ---- //
464 
465  Vertex theRecoVtx = theRecoVertex;
466 
467  double chi2 = theRecoVtx.chi2();
468  double ndf = theRecoVtx.ndof();
469 
470  if (chi2 > TMath::ChisquareQuantile(0.95, ndf)) {
471  if (debug_)
472  cout << "Rejected because of chi2 = " << chi2 << " ndf = " << ndf
473  << " confid. level: " << TMath::ChisquareQuantile(0.95, ndf) << endl;
474  return true;
475  }
476 
477  // ---- Create and fill vector of refitted TransientTracks ---- //
478 
479  // -----------------------------------------------//
480  // ---- Prepare and Fill the Displaced Vertex ----//
481  // -----------------------------------------------//
482 
483  displacedVertex = theRecoVtx;
484  displacedVertex.removeTracks();
485 
486  for (unsigned i = 0; i < transTracks.size(); i++) {
487  PFTrackHitFullInfo pattern = hitPattern_.analyze(tkerTopo_, tkerGeom_, transTracksRef[i], theRecoVertex);
488 
490 
491  Track refittedTrack;
492  float weight = theRecoVertex.trackWeight(transTracks[i]);
493 
494  if (weight < minAdaptWeight_)
495  continue;
496 
497  try {
498  refittedTrack = theRecoVertex.refittedTrack(transTracks[i]).track();
499  } catch (cms::Exception& exception) {
500  continue;
501  }
502 
503  if (debug_) {
504  cout << "Vertex Track Type = " << vertexTrackType << endl;
505 
506  cout << "nHitBeforeVertex = " << pattern.first.first << " nHitAfterVertex = " << pattern.second.first
507  << " nMissHitBeforeVertex = " << pattern.first.second << " nMissHitAfterVertex = " << pattern.second.second
508  << " Weight = " << weight << endl;
509  }
510 
511  displacedVertex.addElement(transTracksRef[i], refittedTrack, pattern, vertexTrackType, weight);
512  }
513 
514  displacedVertex.setPrimaryDirection(helper_.primaryVertex());
515  displacedVertex.calcKinematics();
516 
517  if (debug_)
518  cout << "== End vertexing procedure ==" << endl;
519 
520  return false;
521 }

References reco::PFDisplacedVertex::addElement(), PFCheckHitPattern::analyze(), reco::PFDisplacedVertex::calcKinematics(), hltPixelTracks_cff::chi2, reco::Vertex::chi2(), gather_cfg::cout, debug_, reco::PFDisplacedVertexSeed::Dump(), reco::PFDisplacedVertexSeed::elements(), cppFunctionSkipper::exception, F_ADAPTIVE, F_DONOTREFIT, F_KALMAN, F_NOTDEFINED, getVertexTrackType(), globTkGeomHandle_, helper_, PFTrackAlgoTools::highQuality(), hitPattern_, mps_fire::i, PFDisplacedVertexHelper::isTrackSelected(), TransientVertex::isValid(), magField_, minAdaptWeight_, reco::Vertex::ndof(), PFTrackAlgoTools::nonIterative(), topSingleLeptonDQM_PU_cfi::pattern, TransientVertex::position(), PFDisplacedVertexHelper::primaryVertex(), primaryVertexCut_, ratio_, TransientVertex::refittedTrack(), reco::Vertex::removeTracks(), rho, reco::PFDisplacedVertexSeed::seedPoint(), reco::PFDisplacedVertex::setPrimaryDirection(), sigmacut_, mathSSE::sqrt(), PFTrackAlgoTools::step45(), switchOff2TrackVertex_, t_ini_, tecCut_, tkerGeom_, tkerTopo_, tobCut_, reco::TransientTrack::track(), HLT_FULL_cff::track, TransientVertex::trackWeight(), KalmanVertexFitter::vertex(), AdaptiveVertexFitter::vertex(), extraflags_cff::vtx, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by findDisplacedVertices().

◆ getTransvLongDiff()

std::pair< float, float > PFDisplacedVertexFinder::getTransvLongDiff ( const GlobalPoint Ref,
const GlobalPoint ToProject 
) const
private

Definition at line 630 of file PFDisplacedVertexFinder.cc.

631  {
632  const auto& vRef = Ref.basicVector();
633  const auto& vToProject = ToProject.basicVector();
634  float vRefMag2 = vRef.mag2();
635  float oneOverMag = 1.0f / sqrt(vRefMag2);
636 
637  return std::make_pair(fabs(vRef.cross(vToProject).mag() * oneOverMag),
638  fabs((vRef.dot(vToProject) - vRefMag2) * oneOverMag));
639 }

References PV3DBase< T, PVType, FrameType >::basicVector(), Basic3DVector< T >::mag2(), and mathSSE::sqrt().

Referenced by findSeedsFromCandidate(), and isCloseTo().

◆ getVertexTrackType()

reco::PFDisplacedVertex::VertexTrackType PFDisplacedVertexFinder::getVertexTrackType ( PFTrackHitFullInfo pairTrackHitInfo) const
private

Definition at line 641 of file PFDisplacedVertexFinder.cc.

642  {
643  unsigned int nHitBeforeVertex = pairTrackHitInfo.first.first;
644  unsigned int nHitAfterVertex = pairTrackHitInfo.second.first;
645 
646  unsigned int nMissHitBeforeVertex = pairTrackHitInfo.first.second;
647  unsigned int nMissHitAfterVertex = pairTrackHitInfo.second.second;
648 
649  // For the moment those definitions are given apriori a more general study would be useful to refine those criteria
650 
651  if (nHitBeforeVertex <= 1 && nHitAfterVertex >= 3 && nMissHitAfterVertex <= 1)
652  return PFDisplacedVertex::T_FROM_VERTEX;
653  else if (nHitBeforeVertex >= 3 && nHitAfterVertex <= 1 && nMissHitBeforeVertex <= 1)
654  return PFDisplacedVertex::T_TO_VERTEX;
655  else if ((nHitBeforeVertex >= 2 && nHitAfterVertex >= 3) || (nHitBeforeVertex >= 3 && nHitAfterVertex >= 2))
656  return PFDisplacedVertex::T_MERGED;
657  else
658  return PFDisplacedVertex::T_NOT_FROM_VERTEX;
659 }

Referenced by fitVertexFromSeed().

◆ isCloseTo()

bool PFDisplacedVertexFinder::isCloseTo ( const reco::PFDisplacedVertexSeed dv1,
const reco::PFDisplacedVertexSeed dv2 
) const
private

-----— Tools -----— ///

Definition at line 616 of file PFDisplacedVertexFinder.cc.

616  {
617  const GlobalPoint& vP1 = dv1.seedPoint();
618  const GlobalPoint& vP2 = dv2.seedPoint();
619 
620  std::pair<float, float> diffs = getTransvLongDiff(vP1, vP2);
621  if (diffs.second > longSize_)
622  return false;
623  if (diffs.first > transvSize_)
624  return false;
625  // if (Delta_Long < longSize_ && Delta_Transv < transvSize_) isCloseTo = true;
626 
627  return true;
628 }

References getTransvLongDiff(), longSize_, reco::PFDisplacedVertexSeed::seedPoint(), and transvSize_.

Referenced by mergeSeeds().

◆ mergeSeeds()

void PFDisplacedVertexFinder::mergeSeeds ( reco::PFDisplacedVertexSeedCollection tempDisplacedVertexSeeds,
std::vector< bool > &  bLocked 
)
private

Sometimes two vertex candidates can be quite close and coming from the same vertex.

Definition at line 189 of file PFDisplacedVertexFinder.cc.

190  {
191  // loop over displaced vertex candidates
192  // and merge them if they are close to each other
193  for (unsigned idv_mother = 0; idv_mother < tempDisplacedVertexSeeds.size(); idv_mother++) {
194  if (!bLocked[idv_mother]) {
195  for (unsigned idv_daughter = idv_mother + 1; idv_daughter < tempDisplacedVertexSeeds.size(); idv_daughter++) {
196  if (!bLocked[idv_daughter]) {
197  if (isCloseTo(tempDisplacedVertexSeeds[idv_mother], tempDisplacedVertexSeeds[idv_daughter])) {
198  tempDisplacedVertexSeeds[idv_mother].mergeWith(tempDisplacedVertexSeeds[idv_daughter]);
199  bLocked[idv_daughter] = true;
200  if (debug_)
201  cout << "Seeds " << idv_mother << " and " << idv_daughter << " merged" << endl;
202  }
203  }
204  }
205  }
206  }
207 }

References gather_cfg::cout, debug_, and isCloseTo().

Referenced by findDisplacedVertices().

◆ rejectAndLabelVertex()

bool PFDisplacedVertexFinder::rejectAndLabelVertex ( reco::PFDisplacedVertex dv)
private

◆ selectAndLabelVertices()

void PFDisplacedVertexFinder::selectAndLabelVertices ( reco::PFDisplacedVertexCollection tempDisplacedVertices,
std::vector< bool > &  bLocked 
)
private

Remove potentially fakes displaced vertices.

Definition at line 523 of file PFDisplacedVertexFinder.cc.

524  {
525  if (debug_)
526  cout << " 4.1) Reject vertices " << endl;
527 
528  for (unsigned idv = 0; idv < tempDisplacedVertices.size(); idv++) {
529  // ---- We accept a vertex only if it is not in TOB in the barrel
530  // and not in TEC in the end caps
531  // and not too much before the first pixel layer ---- //
532 
533  const float rho = tempDisplacedVertices[idv].position().rho();
534  const float z = tempDisplacedVertices[idv].position().z();
535 
536  if (rho > tobCut_ || rho < primaryVertexCut_ || fabs(z) > tecCut_) {
537  if (debug_)
538  cout << "Vertex " << idv << " geometrically rejected #rho = " << rho << " z = " << z << endl;
539 
540  bLocked[idv] = true;
541 
542  continue;
543  }
544 
545  unsigned nPrimary = tempDisplacedVertices[idv].nPrimaryTracks();
546  unsigned nMerged = tempDisplacedVertices[idv].nMergedTracks();
547  unsigned nSecondary = tempDisplacedVertices[idv].nSecondaryTracks();
548 
549  if (nPrimary + nMerged > 1) {
550  bLocked[idv] = true;
551  if (debug_)
552  cout << "Vertex " << idv << " rejected because two primary or merged tracks" << endl;
553  }
554 
555  if (nPrimary + nMerged + nSecondary < 2) {
556  bLocked[idv] = true;
557  if (debug_)
558  cout << "Vertex " << idv << " rejected because only one track related to the vertex" << endl;
559  }
560  }
561 
562  if (debug_)
563  cout << " 4.2) Check for common vertices" << endl;
564 
565  // ---- Among the remaining vertex we shall remove one
566  // of those which have two common tracks ---- //
567 
568  for (unsigned idv_mother = 0; idv_mother < tempDisplacedVertices.size(); idv_mother++) {
569  for (unsigned idv_daughter = idv_mother + 1; idv_daughter < tempDisplacedVertices.size(); idv_daughter++) {
570  if (!bLocked[idv_daughter] && !bLocked[idv_mother]) {
571  const unsigned commonTrks =
572  commonTracks(tempDisplacedVertices[idv_daughter], tempDisplacedVertices[idv_mother]);
573 
574  if (commonTrks > 1) {
575  if (debug_)
576  cout << "Vertices " << idv_daughter << " and " << idv_mother << " has many common tracks" << endl;
577 
578  // We keep the vertex vertex which contains the most of the tracks
579 
580  const int mother_size = tempDisplacedVertices[idv_mother].nTracks();
581  const int daughter_size = tempDisplacedVertices[idv_daughter].nTracks();
582 
583  if (mother_size > daughter_size)
584  bLocked[idv_daughter] = true;
585  else if (mother_size < daughter_size)
586  bLocked[idv_mother] = true;
587  else {
588  // If they have the same size we keep the vertex with the best normalised chi2
589 
590  const float mother_normChi2 = tempDisplacedVertices[idv_mother].normalizedChi2();
591  const float daughter_normChi2 = tempDisplacedVertices[idv_daughter].normalizedChi2();
592  if (mother_normChi2 < daughter_normChi2)
593  bLocked[idv_daughter] = true;
594  else
595  bLocked[idv_mother] = true;
596  }
597  }
598  }
599  }
600  }
601 
602  for (unsigned idv = 0; idv < tempDisplacedVertices.size(); idv++)
603  if (!bLocked[idv])
604  bLocked[idv] = rejectAndLabelVertex(tempDisplacedVertices[idv]);
605 }

References commonTracks(), gather_cfg::cout, debug_, primaryVertexCut_, rejectAndLabelVertex(), rho, tecCut_, tobCut_, and z.

Referenced by findDisplacedVertices().

◆ setAVFParameters()

void PFDisplacedVertexFinder::setAVFParameters ( const edm::ParameterSet ps)
inline

Definition at line 90 of file PFDisplacedVertexFinder.h.

90  {
91  sigmacut_ = ps.getParameter<double>("sigmacut");
92  t_ini_ = ps.getParameter<double>("Tini");
93  ratio_ = ps.getParameter<double>("ratio");
94  }

References edm::ParameterSet::getParameter(), ratio_, sigmacut_, and t_ini_.

◆ setDebug()

void PFDisplacedVertexFinder::setDebug ( bool  debug)
inline

Sets debug printout flag.

Definition at line 68 of file PFDisplacedVertexFinder.h.

68 { debug_ = debug; }

References debug, and debug_.

◆ setEdmParameters()

void PFDisplacedVertexFinder::setEdmParameters ( const MagneticField magField,
edm::ESHandle< GlobalTrackingGeometry globTkGeomHandle,
const TrackerTopology tkerTopo,
const TrackerGeometry tkerGeom 
)
inline

Sets parameters for track extrapolation and hits study.

Definition at line 71 of file PFDisplacedVertexFinder.h.

74  {
75  magField_ = magField;
76  globTkGeomHandle_ = globTkGeomHandle;
77  tkerTopo_ = tkerTopo;
78  tkerGeom_ = tkerGeom;
79  }

References globTkGeomHandle_, magField_, tkerGeom_, and tkerTopo_.

◆ setInput()

void PFDisplacedVertexFinder::setInput ( const edm::Handle< reco::PFDisplacedVertexCandidateCollection > &  displacedVertexCandidates)

Set input collections of tracks.

Definition at line 43 of file PFDisplacedVertexFinder.cc.

44  {
45  if (displacedVertexCandidates.isValid()) {
46  displacedVertexCandidates_ = displacedVertexCandidates.product();
47  } else {
49  }
50 }

References displacedVertexCandidates_, edm::HandleBase::isValid(), and edm::Handle< T >::product().

◆ setParameters()

void PFDisplacedVertexFinder::setParameters ( double  transvSize,
double  longSize,
double  primaryVertexCut,
double  tobCut,
double  tecCut,
double  minAdaptWeight,
bool  switchOff2TrackVertex 
)
inline

◆ setPrimaryVertex()

void PFDisplacedVertexFinder::setPrimaryVertex ( edm::Handle< reco::VertexCollection mainVertexHandle,
edm::Handle< reco::BeamSpot beamSpotHandle 
)
inline

Definition at line 85 of file PFDisplacedVertexFinder.h.

86  {
87  helper_.setPrimaryVertex(mainVertexHandle, beamSpotHandle);
88  }

References helper_, and PFDisplacedVertexHelper::setPrimaryVertex().

◆ setTracksSelector()

void PFDisplacedVertexFinder::setTracksSelector ( const edm::ParameterSet ps)
inline

◆ setVertexIdentifier()

void PFDisplacedVertexFinder::setVertexIdentifier ( const edm::ParameterSet ps)
inline

◆ transferDisplacedVertices()

std::unique_ptr<reco::PFDisplacedVertexCollection> PFDisplacedVertexFinder::transferDisplacedVertices ( )
inline
Returns
unique_ptr to collection of DisplacedVertices

Definition at line 100 of file PFDisplacedVertexFinder.h.

100  {
102  }

References displacedVertices_, and eostools::move().

Friends And Related Function Documentation

◆ operator<<

std::ostream& operator<< ( std::ostream &  out,
const PFDisplacedVertexFinder a 
)
friend

Definition at line 680 of file PFDisplacedVertexFinder.cc.

680  {
681  if (!out)
682  return out;
683  out << setprecision(3) << setw(5) << endl;
684  out << "" << endl;
685  out << " ====================================== " << endl;
686  out << " ====== Displaced Vertex Finder ======= " << endl;
687  out << " ====================================== " << endl;
688  out << " " << endl;
689 
690  a.helper_.Dump();
691  out << "" << endl
692  << " Adaptive Vertex Fitter parameters are :" << endl
693  << " sigmacut = " << a.sigmacut_ << " T_ini = " << a.t_ini_ << " ratio = " << a.ratio_ << endl
694  << endl;
695 
696  const std::unique_ptr<reco::PFDisplacedVertexCollection>& displacedVertices_ = std::move(a.displacedVertices());
697 
698  if (!displacedVertices_.get()) {
699  out << "displacedVertex already transfered" << endl;
700  } else {
701  out << "Number of displacedVertices found : " << displacedVertices_->size() << endl << endl;
702 
703  int i = -1;
704 
705  for (PFDisplacedVertexFinder::IDV idv = displacedVertices_->begin(); idv != displacedVertices_->end(); idv++) {
706  i++;
707  out << i << " ";
708  idv->Dump();
709  out << "" << endl;
710  }
711  }
712 
713  return out;
714 }

Member Data Documentation

◆ debug_

bool PFDisplacedVertexFinder::debug_
private

If true, debug printouts activated.

Definition at line 166 of file PFDisplacedVertexFinder.h.

Referenced by findDisplacedVertices(), findSeedsFromCandidate(), fitVertexFromSeed(), mergeSeeds(), selectAndLabelVertices(), and setDebug().

◆ displacedVertexCandidates_

reco::PFDisplacedVertexCandidateCollection const* PFDisplacedVertexFinder::displacedVertexCandidates_
private

-----— Members -----— ///

Definition at line 143 of file PFDisplacedVertexFinder.h.

Referenced by findDisplacedVertices(), and setInput().

◆ displacedVertices_

std::unique_ptr<reco::PFDisplacedVertexCollection> PFDisplacedVertexFinder::displacedVertices_
private

◆ globTkGeomHandle_

edm::ESHandle<GlobalTrackingGeometry> PFDisplacedVertexFinder::globTkGeomHandle_
private

Tracker geometry for discerning hit positions.

Definition at line 169 of file PFDisplacedVertexFinder.h.

Referenced by fitVertexFromSeed(), and setEdmParameters().

◆ helper_

PFDisplacedVertexHelper PFDisplacedVertexFinder::helper_
private

◆ hitPattern_

PFCheckHitPattern PFDisplacedVertexFinder::hitPattern_
private

Definition at line 178 of file PFDisplacedVertexFinder.h.

Referenced by fitVertexFromSeed().

◆ longSize_

float PFDisplacedVertexFinder::longSize_
private

Definition at line 151 of file PFDisplacedVertexFinder.h.

Referenced by findSeedsFromCandidate(), isCloseTo(), and setParameters().

◆ magField_

const MagneticField* PFDisplacedVertexFinder::magField_
private

to be able to extrapolate tracks f

Definition at line 176 of file PFDisplacedVertexFinder.h.

Referenced by fitVertexFromSeed(), and setEdmParameters().

◆ minAdaptWeight_

double PFDisplacedVertexFinder::minAdaptWeight_
private

Definition at line 155 of file PFDisplacedVertexFinder.h.

Referenced by fitVertexFromSeed(), and setParameters().

◆ primaryVertexCut_

double PFDisplacedVertexFinder::primaryVertexCut_
private

◆ ratio_

double PFDisplacedVertexFinder::ratio_
private

Definition at line 163 of file PFDisplacedVertexFinder.h.

Referenced by fitVertexFromSeed(), and setAVFParameters().

◆ sigmacut_

double PFDisplacedVertexFinder::sigmacut_
private

Adaptive Vertex Fitter parameters.

Definition at line 161 of file PFDisplacedVertexFinder.h.

Referenced by fitVertexFromSeed(), and setAVFParameters().

◆ switchOff2TrackVertex_

bool PFDisplacedVertexFinder::switchOff2TrackVertex_
private

Definition at line 157 of file PFDisplacedVertexFinder.h.

Referenced by fitVertexFromSeed(), and setParameters().

◆ t_ini_

double PFDisplacedVertexFinder::t_ini_
private

Definition at line 162 of file PFDisplacedVertexFinder.h.

Referenced by fitVertexFromSeed(), and setAVFParameters().

◆ tecCut_

double PFDisplacedVertexFinder::tecCut_
private

◆ tkerGeom_

const TrackerGeometry* PFDisplacedVertexFinder::tkerGeom_
private

Definition at line 173 of file PFDisplacedVertexFinder.h.

Referenced by fitVertexFromSeed(), and setEdmParameters().

◆ tkerTopo_

const TrackerTopology* PFDisplacedVertexFinder::tkerTopo_
private

doc?

Definition at line 172 of file PFDisplacedVertexFinder.h.

Referenced by fitVertexFromSeed(), and setEdmParameters().

◆ tobCut_

double PFDisplacedVertexFinder::tobCut_
private

◆ transvSize_

float PFDisplacedVertexFinder::transvSize_
private

-----— Parameters -----— ///

Algo parameters for the vertex finder

Definition at line 150 of file PFDisplacedVertexFinder.h.

Referenced by findSeedsFromCandidate(), isCloseTo(), and setParameters().

AdaptiveVertexFitter
Definition: AdaptiveVertexFitter.h:29
PFDisplacedVertexFinder::globTkGeomHandle_
edm::ESHandle< GlobalTrackingGeometry > globTkGeomHandle_
Tracker geometry for discerning hit positions.
Definition: PFDisplacedVertexFinder.h:169
PFTrackAlgoTools::nonIterative
bool nonIterative(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:210
PFDisplacedVertexFinder::FitterType
FitterType
Fitter Type.
Definition: PFDisplacedVertexFinder.h:46
PFDisplacedVertexFinder::getTransvLongDiff
std::pair< float, float > getTransvLongDiff(const GlobalPoint &, const GlobalPoint &) const
Definition: PFDisplacedVertexFinder.cc:630
PFDisplacedVertexFinder::longSize_
float longSize_
Definition: PFDisplacedVertexFinder.h:151
reco::Vertex::originalTrack
TrackBaseRef originalTrack(const Track &refTrack) const
Definition: Vertex.cc:96
KalmanVertexUpdator< 5 >
mps_fire.i
i
Definition: mps_fire.py:428
HLT_FULL_cff.track
track
Definition: HLT_FULL_cff.py:11779
edm::Handle::product
T const * product() const
Definition: Handle.h:70
reco::PFDisplacedVertexSeed::seedPoint
const GlobalPoint & seedPoint() const
Definition: PFDisplacedVertexSeed.h:53
particleFlowDisplacedVertex_cfi.minAdaptWeight
minAdaptWeight
Definition: particleFlowDisplacedVertex_cfi.py:34
edm::RefToBase::key
size_t key() const
Definition: RefToBase.h:219
TransientVertex::isValid
bool isValid() const
Definition: TransientVertex.h:195
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
reco::PFDisplacedVertex::isFake
bool isFake() const
Definition: PFDisplacedVertex.h:204
TransientVertex::refittedTrack
reco::TransientTrack refittedTrack(const reco::TransientTrack &track) const
Definition: TransientVertex.cc:291
PFTrackAlgoTools::step45
bool step45(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:221
PFDisplacedVertexFinder::getVertexTrackType
reco::PFDisplacedVertex::VertexTrackType getVertexTrackType(PFTrackHitFullInfo &) const
Definition: PFDisplacedVertexFinder.cc:641
PFDisplacedVertexFinder::helper_
PFDisplacedVertexHelper helper_
Definition: PFDisplacedVertexFinder.h:180
gather_cfg.cout
cout
Definition: gather_cfg.py:144
reco::PFDisplacedVertexSeed
Block of elements.
Definition: PFDisplacedVertexSeed.h:25
PFDisplacedVertexFinder::IDVS
reco::PFDisplacedVertexSeedCollection::iterator IDVS
-----— Useful Types -----— ///
Definition: PFDisplacedVertexFinder.h:39
reco::PFDisplacedVertexCandidate::tref
const TrackBaseRef & tref(unsigned ie) const
Definition: PFDisplacedVertexCandidate.h:92
PFDisplacedVertexHelper::setVertexIdentifier
void setVertexIdentifier(const edm::ParameterSet &ps)
Set Vertex identifier parameters.
Definition: PFDisplacedVertexHelper.h:26
particleFlowDisplacedVertex_cfi.transvSize
transvSize
Definition: particleFlowDisplacedVertex_cfi.py:14
particleFlowDisplacedVertex_cfi.tecCut
tecCut
Definition: particleFlowDisplacedVertex_cfi.py:29
TransientVertex::position
GlobalPoint position() const
Definition: TransientVertex.h:169
PFDisplacedVertexHelper::setTracksSelector
void setTracksSelector(const edm::ParameterSet &ps)
Set Tracks selector parameters.
Definition: PFDisplacedVertexHelper.h:23
PFDisplacedVertexHelper::setPrimaryVertex
void setPrimaryVertex(edm::Handle< reco::VertexCollection > mainVertexHandle, edm::Handle< reco::BeamSpot > beamSpotHandle)
Update the primary vertex information.
Definition: PFDisplacedVertexHelper.cc:23
reco::PFDisplacedVertexCandidate::dcaPoint
const GlobalPoint dcaPoint(unsigned ie1, unsigned ie2) const
Definition: PFDisplacedVertexCandidate.cc:126
PFDisplacedVertexFinder::tkerTopo_
const TrackerTopology * tkerTopo_
doc?
Definition: PFDisplacedVertexFinder.h:172
PFDisplacedVertexFinder::debug_
bool debug_
If true, debug printouts activated.
Definition: PFDisplacedVertexFinder.h:166
PFDisplacedVertexFinder::isCloseTo
bool isCloseTo(const reco::PFDisplacedVertexSeed &, const reco::PFDisplacedVertexSeed &) const
-----— Tools -----— ///
Definition: PFDisplacedVertexFinder.cc:616
PFTrackAlgoTools::highQuality
bool highQuality(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:193
PFDisplacedVertexFinder::fitVertexFromSeed
bool fitVertexFromSeed(const reco::PFDisplacedVertexSeed &, reco::PFDisplacedVertex &)
Fit one by one the vertex points with associated tracks to get displaced vertices.
Definition: PFDisplacedVertexFinder.cc:209
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
hltPixelTracks_cff.chi2
chi2
Definition: hltPixelTracks_cff.py:25
PFDisplacedVertexFinder::rejectAndLabelVertex
bool rejectAndLabelVertex(reco::PFDisplacedVertex &dv)
Definition: PFDisplacedVertexFinder.cc:607
PFDisplacedVertexFinder::F_KALMAN
Definition: PFDisplacedVertexFinder.h:46
reco::PFDisplacedVertexSeed::Dump
void Dump(std::ostream &out=std::cout) const
cout function
Definition: PFDisplacedVertexSeed.cc:63
Basic3DVector::mag2
T mag2() const
The vector magnitude squared. Equivalent to vec.dot(vec)
Definition: extBasic3DVector.h:113
PFDisplacedVertexFinder::displacedVertices_
std::unique_ptr< reco::PFDisplacedVertexCollection > displacedVertices_
Definition: PFDisplacedVertexFinder.h:144
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
reco::PFDisplacedVertexCollection
std::vector< PFDisplacedVertex > PFDisplacedVertexCollection
collection of PFDisplacedVertex objects
Definition: PFDisplacedVertexFwd.h:10
PFDisplacedVertexFinder::primaryVertexCut_
double primaryVertexCut_
Definition: PFDisplacedVertexFinder.h:152
PFDisplacedVertexFinder::displacedVertexCandidates_
reco::PFDisplacedVertexCandidateCollection const * displacedVertexCandidates_
-----— Members -----— ///
Definition: PFDisplacedVertexFinder.h:143
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
reco::PFDisplacedVertexCandidate::DistMap
std::map< float, std::pair< int, int > > DistMap
Definition: PFDisplacedVertexCandidate.h:49
particleFlowDisplacedVertex_cfi.tobCut
tobCut
Definition: particleFlowDisplacedVertex_cfi.py:25
PFDisplacedVertexFinder::F_ADAPTIVE
Definition: PFDisplacedVertexFinder.h:46
cmsdt::algo
algo
Definition: constants.h:164
PFDisplacedVertexFinder::switchOff2TrackVertex_
bool switchOff2TrackVertex_
Definition: PFDisplacedVertexFinder.h:157
PFDisplacedVertexFinder::F_NOTDEFINED
Definition: PFDisplacedVertexFinder.h:46
debug
#define debug
Definition: HDRShower.cc:19
PFDisplacedVertexFinder::tobCut_
double tobCut_
Definition: PFDisplacedVertexFinder.h:153
PFDisplacedVertexFinder::PFTrackHitFullInfo
std::pair< PFTrackHitInfo, PFTrackHitInfo > PFTrackHitFullInfo
Definition: PFDisplacedVertexFinder.h:43
reco::PFDisplacedVertexSeedCollection
std::vector< PFDisplacedVertexSeed > PFDisplacedVertexSeedCollection
collection of PFDisplacedVertexSeed objects
Definition: PFDisplacedVertexSeedFwd.h:10
PFCheckHitPattern::analyze
PFTrackHitFullInfo analyze(const TrackerTopology *tkerTopo, const TrackerGeometry *tkerGeom, const reco::TrackBaseRef track, const TransientVertex &vert)
Definition: PFCheckHitPattern.cc:77
KalmanVertexTrackCompatibilityEstimator
Definition: KalmanVertexTrackCompatibilityEstimator.h:24
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
reco::PFDisplacedVertexSeed::elements
const std::vector< TrackBaseRef > & elements() const
Definition: PFDisplacedVertexSeed.h:48
PFDisplacedVertexFinder::findSeedsFromCandidate
void findSeedsFromCandidate(const reco::PFDisplacedVertexCandidate &, reco::PFDisplacedVertexSeedCollection &)
-----— Different steps of the finder algorithm -----— ///
Definition: PFDisplacedVertexFinder.cc:141
DDAxes::z
reco::Track
Definition: Track.h:27
PFDisplacedVertexHelper::primaryVertex
math::XYZPoint primaryVertex() const
Set Vertex direction using the primary vertex.
Definition: PFDisplacedVertexHelper.h:39
PFDisplacedVertexFinder::magField_
const MagneticField * magField_
to be able to extrapolate tracks f
Definition: PFDisplacedVertexFinder.h:176
Point3DBase< float, GlobalTag >
reco::PFDisplacedVertexCandidate::r2Map
DistMap r2Map() const
-----— Provide useful information -----— ///
Definition: PFDisplacedVertexCandidate.cc:77
KalmanVertexSmoother
Definition: KalmanVertexSmoother.h:6
reco::Vertex::refittedTracks
const std::vector< Track > & refittedTracks() const
Returns the container of refitted tracks.
Definition: Vertex.h:183
reco::PFDisplacedVertex::setVertexType
void setVertexType(VertexType vertexType)
Set the type of this vertex.
Definition: PFDisplacedVertex.h:77
PFDisplacedVertexFinder::t_ini_
double t_ini_
Definition: PFDisplacedVertexFinder.h:162
DDAxes::rho
PFDisplacedVertexFinder::mergeSeeds
void mergeSeeds(reco::PFDisplacedVertexSeedCollection &, std::vector< bool > &bLocked)
Sometimes two vertex candidates can be quite close and coming from the same vertex.
Definition: PFDisplacedVertexFinder.cc:189
cppFunctionSkipper.exception
exception
Definition: cppFunctionSkipper.py:10
reco::PFDisplacedVertex
Block of elements.
Definition: PFDisplacedVertex.h:23
a
double a
Definition: hdecay.h:119
PFDisplacedVertexFinder::IDV
reco::PFDisplacedVertexCollection::iterator IDV
Definition: PFDisplacedVertexFinder.h:40
reco::PFDisplacedVertex::addElement
void addElement(const TrackBaseRef &r, const Track &refTrack, const PFTrackHitFullInfo &hitInfo, VertexTrackType trackType=T_NOT_FROM_VERTEX, float w=1.0)
Add a new track to the vertex.
Definition: PFDisplacedVertex.cc:12
PFDisplacedVertexFinder::ratio_
double ratio_
Definition: PFDisplacedVertexFinder.h:163
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
GeometricAnnealing
Definition: GeometricAnnealing.h:7
particleFlowDisplacedVertex_cfi.longSize
longSize
Definition: particleFlowDisplacedVertex_cfi.py:17
topSingleLeptonDQM_PU_cfi.pattern
pattern
Definition: topSingleLeptonDQM_PU_cfi.py:39
PFDisplacedVertexFinder::commonTracks
unsigned commonTracks(const reco::PFDisplacedVertex &, const reco::PFDisplacedVertex &) const
Definition: PFDisplacedVertexFinder.cc:661
GlobalErrorBase< double, ErrorMatrixTag >
PFDisplacedVertexFinder::selectAndLabelVertices
void selectAndLabelVertices(reco::PFDisplacedVertexCollection &, std::vector< bool > &)
Remove potentially fakes displaced vertices.
Definition: PFDisplacedVertexFinder.cc:523
TransientVertex
Definition: TransientVertex.h:18
PFDisplacedVertexFinder::F_DONOTREFIT
Definition: PFDisplacedVertexFinder.h:46
DefaultLinearizationPointFinder
Definition: DefaultLinearizationPointFinder.h:13
reco::PFDisplacedVertex::VertexTrackType
VertexTrackType
Definition: PFDisplacedVertex.h:36
reco::Vertex::chi2
double chi2() const
chi-squares
Definition: Vertex.h:103
PV3DBase::basicVector
const BasicVectorType & basicVector() const
Definition: PV3DBase.h:53
PFDisplacedVertexFinder::tkerGeom_
const TrackerGeometry * tkerGeom_
Definition: PFDisplacedVertexFinder.h:173
particleFlowDisplacedVertex_cfi.primaryVertexCut
primaryVertexCut
Definition: particleFlowDisplacedVertex_cfi.py:21
reco::PFDisplacedVertex::calcKinematics
void calcKinematics()
Definition: PFDisplacedVertex.h:184
eostools.move
def move(src, dest)
Definition: eostools.py:511
reco::TransientTrack
Definition: TransientTrack.h:19
reco::TransientTrack::track
const Track & track() const
Definition: TransientTrack.h:117
extraflags_cff.vtx
vtx
Definition: extraflags_cff.py:18
reco::PFDisplacedVertex::VertexType
VertexType
Definition: PFDisplacedVertex.h:40
TransientVertex::trackWeight
float trackWeight(const reco::TransientTrack &track) const
Definition: TransientVertex.cc:241
PFDisplacedVertexFinder::transvSize_
float transvSize_
-----— Parameters -----— ///
Definition: PFDisplacedVertexFinder.h:150
reco::JetExtendedAssociation::Ref
edm::Ref< Container > Ref
Definition: JetExtendedAssociation.h:32
PFDisplacedVertexHelper::identifyVertex
reco::PFDisplacedVertex::VertexType identifyVertex(const reco::PFDisplacedVertex &v) const
Vertex identification tool.
Definition: PFDisplacedVertexHelper.cc:72
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
PFDisplacedVertexFinder::tecCut_
double tecCut_
Definition: PFDisplacedVertexFinder.h:154
PFDisplacedVertexFinder::sigmacut_
double sigmacut_
Adaptive Vertex Fitter parameters.
Definition: PFDisplacedVertexFinder.h:161
MillePedeFileConverter_cfg.out
out
Definition: MillePedeFileConverter_cfg.py:31
cms::Exception
Definition: Exception.h:70
reco::PFDisplacedVertex::setPrimaryDirection
void setPrimaryDirection(const math::XYZPoint &pvtx)
Definition: PFDisplacedVertex.cc:54
reco::Vertex::removeTracks
void removeTracks()
Definition: Vertex.cc:90
PFDisplacedVertexFinder::minAdaptWeight_
double minAdaptWeight_
Definition: PFDisplacedVertexFinder.h:155
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
reco::Vertex::ndof
double ndof() const
Definition: Vertex.h:110
PFDisplacedVertexFinder::hitPattern_
PFCheckHitPattern hitPattern_
Definition: PFDisplacedVertexFinder.h:178
particleFlowDisplacedVertex_cfi.switchOff2TrackVertex
switchOff2TrackVertex
Definition: particleFlowDisplacedVertex_cfi.py:38
weight
Definition: weight.py:1
reco::Vertex
Definition: Vertex.h:35
PFDisplacedVertexHelper::isTrackSelected
bool isTrackSelected(const reco::Track &trk, const reco::PFDisplacedVertex::VertexTrackType vertexTrackType) const
Select tracks tool.
Definition: PFDisplacedVertexHelper.cc:38
KalmanVertexFitter
Definition: KalmanVertexFitter.h:22