#include <AlignmentTracksFromVertexSelector.h>

Public Types
typedef std::vector< const reco::Track * >	Tracks

Public Member Functions
	AlignmentTrackFromVertexSelector (const edm::ParameterSet &cfg, edm::ConsumesCollector &iC)
	constructor More...

const reco::Vertex *	findClosestVertex (const reco::TrackCollection &leptonTracks, const reco::VertexCollection *vertices, const edm::EventSetup &setup) const

Tracks	select (const edm::Handle< reco::TrackCollection > &tc, const edm::Event &evt, const edm::EventSetup &setup) const
	select tracks More...

	~AlignmentTrackFromVertexSelector ()
	destructor More...

Private Attributes
edm::EDGetTokenT< reco::TrackCollection >	diLeptonToken_

const edm::ESGetToken< TransientTrackBuilder, TransientTrackRecord >	ttbESToken_

bool	useClosestVertex_

unsigned int	vertexIndex_

edm::EDGetTokenT< reco::VertexCollection >	vertexToken_

Detailed Description

Definition at line 21 of file AlignmentTracksFromVertexSelector.h.

Member Typedef Documentation

◆ Tracks

typedef std::vector<const reco::Track*> AlignmentTrackFromVertexSelector::Tracks

Definition at line 23 of file AlignmentTracksFromVertexSelector.h.

Constructor & Destructor Documentation

◆ AlignmentTrackFromVertexSelector()

AlignmentTrackFromVertexSelector::AlignmentTrackFromVertexSelector	(	const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &	iC
	)

constructor

Definition at line 17 of file AlignmentTracksFromVertexSelector.cc.

     : ttbESToken_(
           iC.esConsumes<TransientTrackBuilder, TransientTrackRecord>(edm::ESInputTag("", "TransientTrackBuilder"))),
       vertexToken_(iC.consumes<reco::VertexCollection>(cfg.getParameter<edm::InputTag>("vertices"))),
       diLeptonToken_(iC.consumes<reco::TrackCollection>(cfg.getParameter<edm::InputTag>("leptonTracks"))),
       useClosestVertex_(cfg.getParameter<bool>("useClosestVertexToDilepton")),
       vertexIndex_(cfg.getParameter<unsigned int>("vertexIndex")) {}

◆ ~AlignmentTrackFromVertexSelector()

AlignmentTrackFromVertexSelector::~AlignmentTrackFromVertexSelector ( )

destructor

Definition at line 27 of file AlignmentTracksFromVertexSelector.cc.

27 {}

Member Function Documentation

◆ findClosestVertex()

const reco::Vertex * AlignmentTrackFromVertexSelector::findClosestVertex	(	const reco::TrackCollection &	leptonTracks,
		const reco::VertexCollection *	vertices,
		const edm::EventSetup &	setup
	)		const

Definition at line 30 of file AlignmentTracksFromVertexSelector.cc.

References TransientTrackBuilder::build(), counter, VertexDistance3D::distance(), TransientVertex::isValid(), singleTopDQM_cfi::setup, HLT_2022v12_cff::track, ttbESToken_, Measurement1D::value(), KalmanVertexFitter::vertex(), AlignmentTracksFromVertexSelector_cfi::vertices, and extraflags_cff::vtx.

Referenced by select().

                                                                                                           {
   reco::Vertex* defaultVtx = nullptr;
 
   // fill the transient track collection with the lepton tracks
   const TransientTrackBuilder* theB = &setup.getData(ttbESToken_);
   std::vector<reco::TransientTrack> tks;
   for (const auto& track : leptonTracks) {
     reco::TransientTrack trajectory = theB->build(track);
     tks.push_back(trajectory);
   }
 
   // compute the secondary vertex
   TransientVertex aTransVtx;
   KalmanVertexFitter kalman(true);
   aTransVtx = kalman.vertex(tks);
 
   if (!aTransVtx.isValid())
     return defaultVtx;
 
   // find the closest vertex to the secondary vertex in 3D
   VertexDistance3D vertTool3D;
   float minD = 9999.;
   int closestVtxIndex = 0;
   int counter = 0;
   for (const auto& vtx : *vertices) {
     double dist3D = vertTool3D.distance(aTransVtx, vtx).value();
     if (dist3D < minD) {
       minD = dist3D;
       closestVtxIndex = counter;
     }
     counter++;
   }
   if ((*vertices).at(closestVtxIndex).isValid()) {
     return &(vertices->at(closestVtxIndex));
   } else {
     return defaultVtx;
   }
 }

◆ select()

AlignmentTrackFromVertexSelector::Tracks AlignmentTrackFromVertexSelector::select	(	const edm::Handle< reco::TrackCollection > &	tc,
		const edm::Event &	evt,
		const edm::EventSetup &	setup
	)		const

select tracks

Definition at line 72 of file AlignmentTracksFromVertexSelector.cc.

References diLeptonToken_, spr::find(), findClosestVertex(), edm::Event::get(), edm::Event::getHandle(), edm::HandleBase::isValid(), crabWrapper::key, edm::Ref< C, T, F >::key(), LogDebug, edm::Handle< T >::product(), mps_fire::result, singleTopDQM_cfi::setup, jetUpdater_cfi::sort, reco::Vertex::tracks_begin(), reco::Vertex::tracks_end(), useClosestVertex_, vertexIndex_, vertexToken_, and AlignmentTracksFromVertexSelector_cfi::vertices.

Referenced by TrackFromVertexConfigSelector::select().

                                                                                                      {
   Tracks result;
 
   std::vector<unsigned int> thePVkeys;
 
   // get collection of reconstructed vertices from event
   edm::Handle<reco::VertexCollection> vertexHandle = evt.getHandle(vertexToken_);
 
   // get collection of the di-lepton traxks
   const auto& leptonTracks = evt.get(diLeptonToken_);
 
   // fill the vector of keys
   if (vertexHandle.isValid()) {
     const reco::VertexCollection* vertices = vertexHandle.product();
     const reco::Vertex* theVtx = nullptr;
 
     if (useClosestVertex_) {
       theVtx = findClosestVertex(leptonTracks, vertices, setup);
     } else {
       if ((*vertices).at(vertexIndex_).isValid()) {
         theVtx = &(vertices->at(vertexIndex_));
       }
     }
 
     if (theVtx) {
       for (auto tv = theVtx->tracks_begin(); tv != theVtx->tracks_end(); tv++) {
         if (tv->isNonnull()) {
           const reco::TrackRef trackRef = tv->castTo<reco::TrackRef>();
           thePVkeys.push_back(trackRef.key());
         }
       }
     }
   }
 
   std::sort(thePVkeys.begin(), thePVkeys.end());
 
   LogDebug("AlignmentTrackFromVertexSelector")
       << "after collecting the PV keys: thePVkeys.size()" << thePVkeys.size() << std::endl;
   for (const auto& key : thePVkeys) {
     LogDebug("AlignmentTrackFromVertexSelector") << key << ", ";
   }
   LogDebug("AlignmentTrackFromVertexSelector") << std::endl;
 
   if (tc.isValid()) {
     int indx(0);
     // put the track in the collection is it was used for the vertex
     for (reco::TrackCollection::const_iterator tk = tc->begin(); tk != tc->end(); ++tk, ++indx) {
       reco::TrackRef trackRef = reco::TrackRef(tc, indx);
       if (std::find(thePVkeys.begin(), thePVkeys.end(), trackRef.key()) != thePVkeys.end()) {
         LogDebug("AlignmentTrackFromVertexSelector") << "track index: " << indx << "filling result vector" << std::endl;
         result.push_back(&(*tk));
       }  // if a valid key is found
     }    // end loop over tracks
   }      // if the handle is valid
   return result;
 }