CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Types | Static Private Member Functions | Private Attributes
SecondaryVertexProducer Class Reference
Inheritance diagram for SecondaryVertexProducer:
edm::EDProducer edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

virtual void produce (edm::Event &event, const edm::EventSetup &es) override
 
 SecondaryVertexProducer (const edm::ParameterSet &params)
 
 ~SecondaryVertexProducer ()
 
- Public Member Functions inherited from edm::EDProducer
 EDProducer ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDProducer ()
 
- Public Member Functions inherited from edm::ProducerBase
 ProducerBase ()
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription
const &)> 		registrationCallback () const
 used by the fwk to register list of products More...
 
virtual ~ProducerBase ()
 
- Public Member Functions inherited from edm::EDConsumerBase
 EDConsumerBase ()
 
ProductHolderIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
std::vector
< ProductHolderIndexAndSkipBit >
const & 		itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Private Types

enum  ConstraintType {
  CONSTRAINT_NONE = 0, CONSTRAINT_BEAMSPOT, CONSTRAINT_PV_BEAMSPOT_SIZE, CONSTRAINT_PV_BS_Z_ERRORS_SCALED,
  CONSTRAINT_PV_ERROR_SCALED, CONSTRAINT_PV_PRIMARIES_IN_FIT
}
 

Static Private Member Functions

static ConstraintType getConstraintType (const std::string &name)
 

Private Attributes

ConstraintType constraint
 
double constraintScaling
 
double extSVDeltaRToJet
 
double minTrackWeight
 
TrackIPTagInfo::SortCriteria sortCriterium
 
edm::EDGetTokenT< reco::BeamSpot > token_BeamSpot
 
edm::EDGetTokenT
< reco::VertexCollection > 		token_extSVCollection
 
edm::EDGetTokenT
< reco::TrackIPTagInfoCollection > 		token_trackIPTagInfo
 
TrackSelector trackSelector
 
bool useExternalSV
 
bool useGhostTrack
 
VertexFilter vertexFilter
 
VertexSorting vertexSorting
 
edm::ParameterSet vtxRecoPSet
 
bool withPVError
 

Additional Inherited Members

- Public Types inherited from edm::EDProducer
typedef EDProducer ModuleType
 
- Public Types inherited from edm::ProducerBase
typedef
ProductRegistryHelper::TypeLabelList 		TypeLabelList
 
- Static Public Member Functions inherited from edm::EDProducer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 65 of file SecondaryVertexProducer.cc.

Member Enumeration Documentation

enum SecondaryVertexProducer::ConstraintType
private
Enumerator
CONSTRAINT_NONE 
CONSTRAINT_BEAMSPOT 
CONSTRAINT_PV_BEAMSPOT_SIZE 
CONSTRAINT_PV_BS_Z_ERRORS_SCALED 
CONSTRAINT_PV_ERROR_SCALED 
CONSTRAINT_PV_PRIMARIES_IN_FIT 

Definition at line 73 of file SecondaryVertexProducer.cc.

73  {
74  CONSTRAINT_NONE = 0,
75  CONSTRAINT_BEAMSPOT,
76  CONSTRAINT_PV_BEAMSPOT_SIZE,
77  CONSTRAINT_PV_BS_Z_ERRORS_SCALED,
78  CONSTRAINT_PV_ERROR_SCALED,
79  CONSTRAINT_PV_PRIMARIES_IN_FIT
80  };

Constructor & Destructor Documentation

SecondaryVertexProducer::SecondaryVertexProducer ( const edm::ParameterSet &  params)
explicit

Definition at line 140 of file SecondaryVertexProducer.cc.

References constraint, CONSTRAINT_BEAMSPOT, CONSTRAINT_PV_BEAMSPOT_SIZE, CONSTRAINT_PV_BS_Z_ERRORS_SCALED, CONSTRAINT_PV_ERROR_SCALED, CONSTRAINT_PV_PRIMARIES_IN_FIT, constraintScaling, edm::ParameterSet::existsAs(), extSVDeltaRToJet, edm::ParameterSet::getParameter(), token_BeamSpot, token_extSVCollection, token_trackIPTagInfo, and useExternalSV.

141  :
142  sortCriterium(TrackSorting::getCriterium(params.getParameter<std::string>("trackSort"))),
143  trackSelector(params.getParameter<edm::ParameterSet>("trackSelection")),
144  constraint(getConstraintType(params.getParameter<std::string>("constraint"))),
145  constraintScaling(1.0),
146  vtxRecoPSet(params.getParameter<edm::ParameterSet>("vertexReco")),
147  useGhostTrack(vtxRecoPSet.getParameter<std::string>("finder") == "gtvr"),
148  withPVError(params.getParameter<bool>("usePVError")),
149  minTrackWeight(params.getParameter<double>("minimumTrackWeight")),
150  vertexFilter(params.getParameter<edm::ParameterSet>("vertexCuts")),
151  vertexSorting(params.getParameter<edm::ParameterSet>("vertexSelection"))
152 {
153  token_trackIPTagInfo = consumes<reco::TrackIPTagInfoCollection>(params.getParameter<edm::InputTag>("trackIPTagInfos"));
154  if (constraint == CONSTRAINT_PV_ERROR_SCALED ||
155  constraint == CONSTRAINT_PV_BS_Z_ERRORS_SCALED)
156  constraintScaling = params.getParameter<double>("pvErrorScaling");
157 
158  if (constraint == CONSTRAINT_PV_BEAMSPOT_SIZE ||
159  constraint == CONSTRAINT_PV_BS_Z_ERRORS_SCALED ||
160  constraint == CONSTRAINT_BEAMSPOT ||
161  constraint == CONSTRAINT_PV_PRIMARIES_IN_FIT )
162  token_BeamSpot = consumes<reco::BeamSpot>(params.getParameter<edm::InputTag>("beamSpotTag"));
163  useExternalSV = false;
164  if(params.existsAs<bool>("useExternalSV")) useExternalSV = params.getParameter<bool> ("useExternalSV");
165  if(useExternalSV) {
166  token_extSVCollection = consumes<reco::VertexCollection>(params.getParameter<edm::InputTag>("extSVCollection"));
167  extSVDeltaRToJet = params.getParameter<double>("extSVDeltaRToJet");
168  }
169  produces<SecondaryVertexTagInfoCollection>();
170 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::existsAs
bool existsAs(std::string const &parameterName, bool trackiness=true) const
checks if a parameter exists as a given type
Definition: ParameterSet.h:184
TrackSorting::getCriterium
reco::TrackIPTagInfo::SortCriteria getCriterium(const std::string &name)
Definition: TrackSorting.cc:11
SecondaryVertexProducer::getConstraintType
static ConstraintType getConstraintType(const std::string &name)
Definition: SecondaryVertexProducer.cc:102
SecondaryVertexProducer::token_extSVCollection
edm::EDGetTokenT< reco::VertexCollection > token_extSVCollection
Definition: SecondaryVertexProducer.cc:98
SecondaryVertexProducer::sortCriterium
TrackIPTagInfo::SortCriteria sortCriterium
Definition: SecondaryVertexProducer.cc:86
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
SecondaryVertexProducer::token_BeamSpot
edm::EDGetTokenT< reco::BeamSpot > token_BeamSpot
Definition: SecondaryVertexProducer.cc:84
SecondaryVertexProducer::token_trackIPTagInfo
edm::EDGetTokenT< reco::TrackIPTagInfoCollection > token_trackIPTagInfo
Definition: SecondaryVertexProducer.cc:85
SecondaryVertexProducer::~SecondaryVertexProducer ( )

Definition at line 172 of file SecondaryVertexProducer.cc.

173 {
174 }

Member Function Documentation

SecondaryVertexProducer::ConstraintType SecondaryVertexProducer::getConstraintType ( const std::string &  name)
staticprivate

Definition at line 102 of file SecondaryVertexProducer.cc.

References edm::hlt::Exception.

103 {
104  if (name == "None")
105  return CONSTRAINT_NONE;
106  else if (name == "BeamSpot")
107  return CONSTRAINT_BEAMSPOT;
108  else if (name == "BeamSpot+PVPosition")
109  return CONSTRAINT_PV_BEAMSPOT_SIZE;
110  else if (name == "BeamSpotZ+PVErrorScaledXY")
111  return CONSTRAINT_PV_BS_Z_ERRORS_SCALED;
112  else if (name == "PVErrorScaled")
113  return CONSTRAINT_PV_ERROR_SCALED;
114  else if (name == "BeamSpot+PVTracksInFit")
115  return CONSTRAINT_PV_PRIMARIES_IN_FIT;
116  else
117  throw cms::Exception("InvalidArgument")
118  << "SecondaryVertexProducer: ``constraint'' parameter "
119  "value \"" << name << "\" not understood."
120  << std::endl;
121 }
void SecondaryVertexProducer::produce ( edm::Event &  event,
const edm::EventSetup &  es 
)
overridevirtual

Implements edm::EDProducer.

Definition at line 211 of file SecondaryVertexProducer.cc.

References SiPixelRawToDigiRegional_cfi::beamSpot, edm::RefVector< C, T, F >::begin(), constraint, CONSTRAINT_BEAMSPOT, CONSTRAINT_NONE, CONSTRAINT_PV_BEAMSPOT_SIZE, CONSTRAINT_PV_BS_Z_ERRORS_SCALED, CONSTRAINT_PV_ERROR_SCALED, CONSTRAINT_PV_PRIMARIES_IN_FIT, constraintScaling, RecoVertex::convertError(), RecoVertex::convertPos(), reco::Vertex::covariance(), Geom::deltaR(), reco::SecondaryVertex::dist2d(), reco::SecondaryVertex::dist3d(), edm::RefVector< C, T, F >::end(), reco::Vertex::error(), edm::hlt::Exception, extSVDeltaRToJet, spr::find(), edm::EventSetup::get(), getGhostTrackFitType(), edm::ParameterSet::getParameter(), ghostTrackES_cfi::ghostTrack, i, customizeTrackingMonitorSeedNumber::idx, getHLTprescales::index, edm::HandleBase::isValid(), j, reco::GhostTrackState::linearize(), python.multivaluedict::map(), minTrackWeight, reco::LeafCandidate::momentum(), reco::Vertex::p4(), reco::Vertex::position(), edm::RefVector< C, T, F >::push_back(), reco::GhostTrackState::setWeight(), sortCriterium, step1_ZMM_7Tev::Status, AlCaHLTBitMon_QueryRunRegistry::string, token_BeamSpot, token_extSVCollection, token_trackIPTagInfo, reco::Vertex::tracks_begin(), reco::Vertex::tracks_end(), trackSelector, reco::Vertex::trackWeight(), Unknown, useExternalSV, useGhostTrack, findQualityFiles::v, vertexFilter, vertexSorting, vtxRecoPSet, withPVError, reco::Vertex::x(), reco::Vertex::y(), and reco::Vertex::z().

213 {
214  typedef std::map<TrackBaseRef, TransientTrack,
215  RefToBaseLess<Track> > TransientTrackMap;
216 
217  edm::ESHandle<TransientTrackBuilder> trackBuilder;
218  es.get<TransientTrackRecord>().get("TransientTrackBuilder",
219  trackBuilder);
220 
221  edm::Handle<TrackIPTagInfoCollection> trackIPTagInfos;
222  event.getByToken(token_trackIPTagInfo, trackIPTagInfos);
223 
224  // External Sec Vertex collection (e.g. for IVF usage)
225  edm::Handle<reco::VertexCollection> extSecVertex;
226  if(useExternalSV) event.getByToken(token_extSVCollection,extSecVertex);
227 
228 
229  edm::Handle<BeamSpot> beamSpot;
230  unsigned int bsCovSrc[7] = { 0, };
231  double sigmaZ = 0.0, beamWidth = 0.0;
232  switch(constraint) {
233  case CONSTRAINT_PV_BEAMSPOT_SIZE:
234  event.getByToken(token_BeamSpot,beamSpot);
235  bsCovSrc[3] = bsCovSrc[4] = bsCovSrc[5] = bsCovSrc[6] = 1;
236  sigmaZ = beamSpot->sigmaZ();
237  beamWidth = beamSpot->BeamWidthX();
238  break;
239 
240  case CONSTRAINT_PV_BS_Z_ERRORS_SCALED:
241  event.getByToken(token_BeamSpot,beamSpot);
242  bsCovSrc[0] = bsCovSrc[1] = 2;
243  bsCovSrc[3] = bsCovSrc[4] = bsCovSrc[5] = 1;
244  sigmaZ = beamSpot->sigmaZ();
245  break;
246 
247  case CONSTRAINT_PV_ERROR_SCALED:
248  bsCovSrc[0] = bsCovSrc[1] = bsCovSrc[2] = 2;
249  break;
250 
251  case CONSTRAINT_BEAMSPOT:
252  case CONSTRAINT_PV_PRIMARIES_IN_FIT:
253  event.getByToken(token_BeamSpot,beamSpot);
254  break;
255 
256  default:
257  /* nothing */;
258  }
259 
260  std::auto_ptr<ConfigurableVertexReconstructor> vertexReco;
261  std::auto_ptr<GhostTrackVertexFinder> vertexRecoGT;
262  if (useGhostTrack)
263  vertexRecoGT.reset(new GhostTrackVertexFinder(
264  vtxRecoPSet.getParameter<double>("maxFitChi2"),
265  vtxRecoPSet.getParameter<double>("mergeThreshold"),
266  vtxRecoPSet.getParameter<double>("primcut"),
267  vtxRecoPSet.getParameter<double>("seccut"),
268  getGhostTrackFitType(vtxRecoPSet.getParameter<std::string>("fitType"))));
269  else
270  vertexReco.reset(
271  new ConfigurableVertexReconstructor(vtxRecoPSet));
272 
273  TransientTrackMap primariesMap;
274 
275  // result secondary vertices
276 
277  std::auto_ptr<SecondaryVertexTagInfoCollection>
278  tagInfos(new SecondaryVertexTagInfoCollection);
279 
280  for(TrackIPTagInfoCollection::const_iterator iterJets =
281  trackIPTagInfos->begin(); iterJets != trackIPTagInfos->end();
282  ++iterJets) {
283  std::vector<SecondaryVertexTagInfo::IndexedTrackData> trackData;
284 // std::cout << "Jet " << iterJets-trackIPTagInfos->begin() << std::endl;
285 
286  const Vertex &pv = *iterJets->primaryVertex();
287 
288  std::set<TransientTrack> primaries;
289  if (constraint == CONSTRAINT_PV_PRIMARIES_IN_FIT) {
290  for(Vertex::trackRef_iterator iter = pv.tracks_begin();
291  iter != pv.tracks_end(); ++iter) {
292  TransientTrackMap::iterator pos =
293  primariesMap.lower_bound(*iter);
294 
295  if (pos != primariesMap.end() &&
296  pos->first == *iter)
297  primaries.insert(pos->second);
298  else {
299  TransientTrack track =
300  trackBuilder->build(
301  iter->castTo<TrackRef>());
302  primariesMap.insert(pos,
303  std::make_pair(*iter, track));
304  primaries.insert(track);
305  }
306  }
307  }
308 
309  edm::RefToBase<Jet> jetRef = iterJets->jet();
310 
311  GlobalVector jetDir(jetRef->momentum().x(),
312  jetRef->momentum().y(),
313  jetRef->momentum().z());
314 
315  std::vector<std::size_t> indices =
316  iterJets->sortedIndexes(sortCriterium);
317 
318  TrackRefVector trackRefs = iterJets->sortedTracks(indices);
319 
320  const std::vector<TrackIPTagInfo::TrackIPData> &ipData =
321  iterJets->impactParameterData();
322 
323  // build transient tracks used for vertex reconstruction
324 
325  std::vector<TransientTrack> fitTracks;
326  std::vector<GhostTrackState> gtStates;
327  std::auto_ptr<GhostTrackPrediction> gtPred;
328  if (useGhostTrack)
329  gtPred.reset(new GhostTrackPrediction(
330  *iterJets->ghostTrack()));
331 
332  for(unsigned int i = 0; i < indices.size(); i++) {
333  typedef SecondaryVertexTagInfo::IndexedTrackData IndexedTrackData;
334 
335  const TrackRef &trackRef = trackRefs[i];
336 
337  trackData.push_back(IndexedTrackData());
338  trackData.back().first = indices[i];
339 
340  // select tracks for SV finder
341 
342  if (!trackSelector(*trackRef, ipData[indices[i]], *jetRef,
343  RecoVertex::convertPos(
344  pv.position()))) {
345  trackData.back().second.svStatus =
346  SecondaryVertexTagInfo::TrackData::trackSelected;
347  continue;
348  }
349 
350  TransientTrackMap::const_iterator pos =
351  primariesMap.find(
352  TrackBaseRef(trackRef));
353  TransientTrack fitTrack;
354  if (pos != primariesMap.end()) {
355  primaries.erase(pos->second);
356  fitTrack = pos->second;
357  } else
358  fitTrack = trackBuilder->build(trackRef);
359  fitTracks.push_back(fitTrack);
360 
361  trackData.back().second.svStatus =
362  SecondaryVertexTagInfo::TrackData::trackUsedForVertexFit;
363 
364  if (useGhostTrack) {
365  GhostTrackState gtState(fitTrack);
366  GlobalPoint pos =
367  ipData[indices[i]].closestToGhostTrack;
368  gtState.linearize(*gtPred, true,
369  gtPred->lambda(pos));
370  gtState.setWeight(ipData[indices[i]].ghostTrackWeight);
371  gtStates.push_back(gtState);
372  }
373  }
374 
375  std::auto_ptr<GhostTrack> ghostTrack;
376  if (useGhostTrack)
377  ghostTrack.reset(new GhostTrack(
378  GhostTrackPrediction(
379  RecoVertex::convertPos(pv.position()),
380  RecoVertex::convertError(pv.error()),
381  GlobalVector(
382  iterJets->ghostTrack()->px(),
383  iterJets->ghostTrack()->py(),
384  iterJets->ghostTrack()->pz()),
385  0.05),
386  *gtPred, gtStates,
387  iterJets->ghostTrack()->chi2(),
388  iterJets->ghostTrack()->ndof()));
389 
390  // perform actual vertex finding
391 
392 
393  std::vector<reco::Vertex> extAssoCollection;
394  std::vector<TransientVertex> fittedSVs;
395  std::vector<SecondaryVertex> SVs;
396  if(!useExternalSV){
397  switch(constraint) {
398  case CONSTRAINT_NONE:
399  if (useGhostTrack)
400  fittedSVs = vertexRecoGT->vertices(
401  pv, *ghostTrack);
402  else
403  fittedSVs = vertexReco->vertices(fitTracks);
404  break;
405 
406  case CONSTRAINT_BEAMSPOT:
407  if (useGhostTrack)
408  fittedSVs = vertexRecoGT->vertices(
409  pv, *beamSpot, *ghostTrack);
410  else
411  fittedSVs = vertexReco->vertices(fitTracks,
412  *beamSpot);
413  break;
414 
415  case CONSTRAINT_PV_BEAMSPOT_SIZE:
416  case CONSTRAINT_PV_BS_Z_ERRORS_SCALED:
417  case CONSTRAINT_PV_ERROR_SCALED: {
418  BeamSpot::CovarianceMatrix cov;
419  for(unsigned int i = 0; i < 7; i++) {
420  unsigned int covSrc = bsCovSrc[i];
421  for(unsigned int j = 0; j < 7; j++) {
422  double v=0.0;
423  if (!covSrc || bsCovSrc[j] != covSrc)
424  v = 0.0;
425  else if (covSrc == 1)
426  v = beamSpot->covariance(i, j);
427  else if (j<3 && i<3)
428  v = pv.covariance(i, j) *
429  constraintScaling;
430  cov(i, j) = v;
431  }
432  }
433 
434  BeamSpot bs(pv.position(), sigmaZ,
435  beamSpot.isValid() ? beamSpot->dxdz() : 0.,
436  beamSpot.isValid() ? beamSpot->dydz() : 0.,
437  beamWidth, cov, BeamSpot::Unknown);
438 
439  if (useGhostTrack)
440  fittedSVs = vertexRecoGT->vertices(
441  pv, bs, *ghostTrack);
442  else
443  fittedSVs = vertexReco->vertices(fitTracks, bs);
444  } break;
445 
446  case CONSTRAINT_PV_PRIMARIES_IN_FIT: {
447  std::vector<TransientTrack> primaries_(
448  primaries.begin(), primaries.end());
449  if (useGhostTrack)
450  fittedSVs = vertexRecoGT->vertices(
451  pv, *beamSpot, primaries_,
452  *ghostTrack);
453  else
454  fittedSVs = vertexReco->vertices(
455  primaries_, fitTracks,
456  *beamSpot);
457  } break;
458  }
459 
460  // build combined SV information and filter
461 
462  SVBuilder svBuilder(pv, jetDir, withPVError);
463  std::remove_copy_if(boost::make_transform_iterator(
464  fittedSVs.begin(), svBuilder),
465  boost::make_transform_iterator(
466  fittedSVs.end(), svBuilder),
467  std::back_inserter(SVs),
468  SVFilter(vertexFilter, pv, jetDir));
469 
470  // clean up now unneeded collections
471  }else{
472  for(size_t iExtSv = 0; iExtSv < extSecVertex->size(); iExtSv++){
473  const reco::Vertex & extVertex = (*extSecVertex)[iExtSv];
474 // GlobalVector vtxDir = GlobalVector(extVertex.p4().X(),extVertex.p4().Y(),extVertex.p4().Z());
475 // if(Geom::deltaR(extVertex.position() - pv.position(), vtxDir)>0.2) continue; //pointing angle
476 // std::cout << " dR " << iExtSv << " " << Geom::deltaR( ( extVertex.position() - pv.position() ), jetDir ) << "eta: " << ( extVertex.position() - pv.position()).eta() << " vs " << jetDir.eta() << " phi: " << ( extVertex.position() - pv.position()).phi() << " vs " << jetDir.phi() << std::endl;
477  if( Geom::deltaR( ( extVertex.position() - pv.position() ), jetDir ) > extSVDeltaRToJet || extVertex.p4().M() < 0.3)
478  continue;
479 // std::cout << " SV added " << iExtSv << std::endl;
480  extAssoCollection.push_back( extVertex );
481  }
482  SVBuilder svBuilder(pv, jetDir, withPVError);
483  std::remove_copy_if(boost::make_transform_iterator( extAssoCollection.begin(), svBuilder),
484  boost::make_transform_iterator(extAssoCollection.end(), svBuilder),
485  std::back_inserter(SVs),
486  SVFilter(vertexFilter, pv, jetDir));
487 
488 
489  }
490 // std::cout << "size: " << SVs.size() << std::endl;
491  gtPred.reset();
492  ghostTrack.reset();
493  gtStates.clear();
494  fitTracks.clear();
495  fittedSVs.clear();
496  extAssoCollection.clear();
497 
498  // sort SVs by importance
499 
500  std::vector<unsigned int> vtxIndices = vertexSorting(SVs);
501 
502  std::vector<SecondaryVertexTagInfo::VertexData> svData;
503 
504  svData.resize(vtxIndices.size());
505  for(unsigned int idx = 0; idx < vtxIndices.size(); idx++) {
506  const SecondaryVertex &sv = SVs[vtxIndices[idx]];
507 
508  svData[idx].vertex = sv;
509  svData[idx].dist2d = sv.dist2d();
510  svData[idx].dist3d = sv.dist3d();
511  svData[idx].direction =
512  GlobalVector(sv.x() - pv.x(),
513  sv.y() - pv.y(),
514  sv.z() - pv.z());
515 
516  // mark tracks successfully used in vertex fit
517 
518  for(Vertex::trackRef_iterator iter = sv.tracks_begin();
519  iter != sv.tracks_end(); iter++) {
520  if (sv.trackWeight(*iter) < minTrackWeight)
521  continue;
522 
523  TrackRefVector::const_iterator pos =
524  std::find(trackRefs.begin(), trackRefs.end(),
525  iter->castTo<TrackRef>());
526 
527  if (pos == trackRefs.end() ) {
528  if(!useExternalSV)
529  throw cms::Exception("TrackNotFound")
530  << "Could not find track from secondary "
531  "vertex in original tracks."
532  << std::endl;
533  } else {
534  unsigned int index = pos - trackRefs.begin();
535  trackData[index].second.svStatus =
536  (SecondaryVertexTagInfo::TrackData::Status)
537  ((unsigned int)SecondaryVertexTagInfo::TrackData::trackAssociatedToVertex + idx);
538  }
539  }
540  }
541 
542  // fill result into tag infos
543 
544  tagInfos->push_back(
545  SecondaryVertexTagInfo(
546  trackData, svData, SVs.size(),
547  TrackIPTagInfoRef(trackIPTagInfos,
548  iterJets - trackIPTagInfos->begin())));
549  }
550 
551  event.put(tagInfos);
552 }
reco::BeamSpot::CovarianceMatrix
math::Error< dimension >::type CovarianceMatrix
Definition: BeamSpot.h:31
RecoVertex::convertPos
reco::Vertex::Point convertPos(const GlobalPoint &p)
Definition: ConvertToFromReco.h:7
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
i
int i
Definition: DBlmapReader.cc:9
SiPixelRawToDigiRegional_cfi.beamSpot
tuple beamSpot
Definition: SiPixelRawToDigiRegional_cfi.py:11
reco::Vertex::tracks_end
trackRef_iterator tracks_end() const
last iterator over tracks
Definition: Vertex.cc:44
reco::SecondaryVertex::dist3d
Measurement1D dist3d() const
Definition: SecondaryVertex.h:23
reco::Vertex::y
double y() const
y coordinate
Definition: Vertex.h:96
RecoVertex::convertError
reco::Vertex::Error convertError(const GlobalError &ge)
Definition: ConvertError.h:8
SecondaryVertexProducer::token_extSVCollection
edm::EDGetTokenT< reco::VertexCollection > token_extSVCollection
Definition: SecondaryVertexProducer.cc:98
reco::Vertex::covariance
double covariance(int i, int j) const
(i, j)-th element of error matrix, i, j = 0, ... 2
Definition: Vertex.h:109
edm::RefVector::end
const_iterator end() const
Termination of iteration.
Definition: RefVector.h:249
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:7
reco::Vertex::position
const Point & position() const
position
Definition: Vertex.h:92
SecondaryVertexProducer::sortCriterium
TrackIPTagInfo::SortCriteria sortCriterium
Definition: SecondaryVertexProducer.cc:86
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
edm::RefVector::begin
const_iterator begin() const
Initialize an iterator over the RefVector.
Definition: RefVector.h:244
SecondaryVertexProducer::token_BeamSpot
edm::EDGetTokenT< reco::BeamSpot > token_BeamSpot
Definition: SecondaryVertexProducer.cc:84
reco::TrackBaseRef
edm::RefToBase< reco::Track > TrackBaseRef
persistent reference to a Track, using views
Definition: TrackFwd.h:22
reco::Vertex::trackWeight
float trackWeight(const TrackBaseRef &r) const
returns the weight with which a Track has contributed to the vertex-fit.
j
int j
Definition: DBlmapReader.cc:9
reco::Vertex::z
double z() const
y coordinate
Definition: Vertex.h:98
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
reco::SecondaryVertex::dist2d
Measurement1D dist2d() const
Definition: SecondaryVertex.h:22
step1_ZMM_7Tev.Status
tuple Status
Definition: step1_ZMM_7Tev.py:112
SecondaryVertexProducer::token_trackIPTagInfo
edm::EDGetTokenT< reco::TrackIPTagInfoCollection > token_trackIPTagInfo
Definition: SecondaryVertexProducer.cc:85
reco::Vertex::x
double x() const
x coordinate
Definition: Vertex.h:94
reco::SecondaryVertexTagInfo::IndexedTrackData
std::pair< unsigned int, TrackData > IndexedTrackData
Definition: SecondaryVertexTagInfo.h:44
customizeTrackingMonitorSeedNumber.idx
tuple idx
DEBUGGING if hasattr(process,&quot;trackMonIterativeTracking2012&quot;): print &quot;trackMonIterativeTracking2012 D...
Definition: customizeTrackingMonitorSeedNumber.py:15
reco::Vertex::p4
math::XYZTLorentzVectorD p4(float mass=0.13957018, float minWeight=0.5) const
Returns the four momentum of the sum of the tracks, assuming the given mass for the decay products...
Definition: Vertex.cc:112
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
reco::LeafCandidate::momentum
virtual Vector momentum() const GCC11_FINAL
spatial momentum vector
Definition: LeafCandidate.h:127
reco::Vertex::error
Error error() const
return SMatrix
Definition: Vertex.h:115
Geom::deltaR
double deltaR(const Vector1 &v1, const Vector2 &v2)
Definition: VectorUtil.h:84
reco::Vertex::trackRef_iterator
std::vector< TrackBaseRef >::const_iterator trackRef_iterator
The iteratator for the vector&lt;TrackRef&gt;
Definition: Vertex.h:37
edm::RefVector::push_back
void push_back(value_type const &ref)
Add a Ref&lt;C, T&gt; to the RefVector.
Definition: RefVector.h:64
reco::Vertex::tracks_begin
trackRef_iterator tracks_begin() const
first iterator over tracks
Definition: Vertex.cc:39
ghostTrackES_cfi.ghostTrack
tuple ghostTrack
Definition: ghostTrackES_cfi.py:5
getGhostTrackFitType
static GhostTrackVertexFinder::FitType getGhostTrackFitType(const std::string &name)
Definition: SecondaryVertexProducer.cc:124
GlobalVector
Global3DVector GlobalVector
Definition: GlobalVector.h:10

Member Data Documentation

ConstraintType SecondaryVertexProducer::constraint
private

Definition at line 88 of file SecondaryVertexProducer.cc.

Referenced by produce(), and SecondaryVertexProducer().

double SecondaryVertexProducer::constraintScaling
private

Definition at line 89 of file SecondaryVertexProducer.cc.

Referenced by produce(), and SecondaryVertexProducer().

double SecondaryVertexProducer::extSVDeltaRToJet
private

Definition at line 97 of file SecondaryVertexProducer.cc.

Referenced by produce(), and SecondaryVertexProducer().

double SecondaryVertexProducer::minTrackWeight
private

Definition at line 93 of file SecondaryVertexProducer.cc.

Referenced by produce().

TrackIPTagInfo::SortCriteria SecondaryVertexProducer::sortCriterium
private

Definition at line 86 of file SecondaryVertexProducer.cc.

Referenced by produce().

edm::EDGetTokenT<reco::BeamSpot> SecondaryVertexProducer::token_BeamSpot
private

Definition at line 84 of file SecondaryVertexProducer.cc.

Referenced by produce(), and SecondaryVertexProducer().

edm::EDGetTokenT<reco::VertexCollection> SecondaryVertexProducer::token_extSVCollection
private

Definition at line 98 of file SecondaryVertexProducer.cc.

Referenced by produce(), and SecondaryVertexProducer().

edm::EDGetTokenT<reco::TrackIPTagInfoCollection> SecondaryVertexProducer::token_trackIPTagInfo
private

Definition at line 85 of file SecondaryVertexProducer.cc.

Referenced by produce(), and SecondaryVertexProducer().

TrackSelector SecondaryVertexProducer::trackSelector
private

Definition at line 87 of file SecondaryVertexProducer.cc.

Referenced by produce().

bool SecondaryVertexProducer::useExternalSV
private

Definition at line 96 of file SecondaryVertexProducer.cc.

Referenced by produce(), and SecondaryVertexProducer().

bool SecondaryVertexProducer::useGhostTrack
private

Definition at line 91 of file SecondaryVertexProducer.cc.

Referenced by produce().

VertexFilter SecondaryVertexProducer::vertexFilter
private

Definition at line 94 of file SecondaryVertexProducer.cc.

Referenced by produce().

VertexSorting SecondaryVertexProducer::vertexSorting
private

Definition at line 95 of file SecondaryVertexProducer.cc.

Referenced by produce().

edm::ParameterSet SecondaryVertexProducer::vtxRecoPSet
private

Definition at line 90 of file SecondaryVertexProducer.cc.

Referenced by produce().

bool SecondaryVertexProducer::withPVError
private

Definition at line 92 of file SecondaryVertexProducer.cc.

Referenced by produce().

Generated for CMSSW Reference Manual by doxygen 1.8.5