calculateVertexSharedTracks.h File Reference

#include "SimDataFormats/Associations/interface/TrackAssociation.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "SimDataFormats/TrackingAnalysis/interface/TrackingVertex.h"

Go to the source code of this file.

Functions
unsigned int	calculateVertexSharedTracks (const reco::Vertex &recoV, const TrackingVertex &simV, const reco::RecoToSimCollection &trackRecoToSimAssociation)
unsigned int	calculateVertexSharedTracks (const TrackingVertex &simV, const reco::Vertex &recoV, const reco::SimToRecoCollection &trackSimToRecoAssociation)

Function Documentation

calculateVertexSharedTracks() [1/2]

unsigned int calculateVertexSharedTracks	(	const reco::Vertex &	recoV,
		const TrackingVertex &	simV,
		const reco::RecoToSimCollection &	trackRecoToSimAssociation
	)

Definition at line 3 of file calculateVertexSharedTracks.cc.

References TrackingVertex::daughterTracks_begin(), TrackingVertex::daughterTracks_end(), edm::AssociationMap< Tag >::end(), edm::AssociationMap< Tag >::find(), newFWLiteAna::found, cmsswSequenceInfo::tp, reco::Vertex::tracks_begin(), and reco::Vertex::tracks_end().

Referenced by VertexAssociatorByPositionAndTracks::associateRecoToSim(), VertexAssociatorByPositionAndTracks::associateSimToReco(), and PrimaryVertexAnalyzer4PUSlimmed::matchReco2SimVertices().

                                                                                                   {
  unsigned int sharedTracks = 0;
  for (auto iTrack = recoV.tracks_begin(); iTrack != recoV.tracks_end(); ++iTrack) {
    auto found = trackRecoToSimAssociation.find(*iTrack);

    if (found == trackRecoToSimAssociation.end())
      continue;

    // matched TP equal to any TP of sim vertex => increase counter
    for (const auto &tp : found->val) {
      if (std::find_if(simV.daughterTracks_begin(), simV.daughterTracks_end(), [&](const TrackingParticleRef &vtp) {
            return tp.first == vtp;
          }) != simV.daughterTracks_end()) {
        sharedTracks += 1;
        break;
      }
    }
  }

  return sharedTracks;
}

calculateVertexSharedTracks() [2/2]

unsigned int calculateVertexSharedTracks	(	const TrackingVertex &	simV,
		const reco::Vertex &	recoV,
		const reco::SimToRecoCollection &	trackSimToRecoAssociation
	)

Definition at line 27 of file calculateVertexSharedTracks.cc.

References TrackingVertex::daughterTracks_begin(), TrackingVertex::daughterTracks_end(), edm::AssociationMap< Tag >::end(), edm::AssociationMap< Tag >::find(), newFWLiteAna::found, reco::Vertex::tracks_begin(), and reco::Vertex::tracks_end().

                                                                                                   {
  unsigned int sharedTracks = 0;
  for (auto iTP = simV.daughterTracks_begin(); iTP != simV.daughterTracks_end(); ++iTP) {
    auto found = trackSimToRecoAssociation.find(*iTP);

    if (found == trackSimToRecoAssociation.end())
      continue;

    // matched track equal to any track of reco vertex => increase counter
    for (const auto &tk : found->val) {
      if (std::find_if(recoV.tracks_begin(), recoV.tracks_end(), [&](const reco::TrackBaseRef &vtk) {
            return tk.first == vtk;
          }) != recoV.tracks_end()) {
        sharedTracks += 1;
        break;
      }
    }
  }

  return sharedTracks;
}