Inheritance diagram for TemplatedVertexMerger< VTX >:

Public Types
typedef std::vector< VTX >	Product

Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Member Functions
void	produce (edm::Event &event, const edm::EventSetup &es) override

	TemplatedVertexMerger (const edm::ParameterSet &params)

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Member Functions
bool	trackFilter (const reco::TrackRef &track) const

Private Attributes
double	maxFraction

double	minSignificance

edm::EDGetTokenT< Product >	token_secondaryVertex

Detailed Description

template<class VTX>
class TemplatedVertexMerger< VTX >

Definition at line 22 of file VertexMerger.cc.

Member Typedef Documentation

template<class VTX >

typedef std::vector<VTX> TemplatedVertexMerger< VTX >::Product

Definition at line 24 of file VertexMerger.cc.

Constructor & Destructor Documentation

template<class VTX >

TemplatedVertexMerger< VTX >::TemplatedVertexMerger ( const edm::ParameterSet & params )

Definition at line 38 of file VertexMerger.cc.

References edm::ParameterSet::getParameter(), and TemplatedVertexMerger< VTX >::token_secondaryVertex.

     : maxFraction(params.getParameter<double>("maxFraction")),
       minSignificance(params.getParameter<double>("minSignificance")) {
   token_secondaryVertex = consumes<Product>(params.getParameter<edm::InputTag>("secondaryVertices"));
   produces<Product>();
 }

Member Function Documentation

template<class VTX >

void TemplatedVertexMerger< VTX >::produce	(	edm::Event &	event,
		const edm::EventSetup &	es
	)

override

Definition at line 46 of file VertexMerger.cc.

References vertexTools::computeSharedTracks(), RecoVertex::convertError(), RecoVertex::convertPos(), VertexDistance3D::distance(), TemplatedVertexMerger< VTX >::maxFraction, TemplatedVertexMerger< VTX >::minSignificance, eostools::move(), indexGen::s2, HLT_2018_cff::secondaryVertices, Measurement1D::significance(), pfDeepBoostedJetPreprocessParams_cfi::sv, and TemplatedVertexMerger< VTX >::token_secondaryVertex.

                                                                                  {
   using namespace reco;
 
   edm::Handle<Product> secondaryVertices;
   event.getByToken(token_secondaryVertex, secondaryVertices);
 
   VertexDistance3D dist;
   auto recoVertices = std::make_unique<Product>();
   for (typename Product::const_iterator sv = secondaryVertices->begin(); sv != secondaryVertices->end(); ++sv) {
     recoVertices->push_back(*sv);
   }
   for (typename Product::iterator sv = recoVertices->begin(); sv != recoVertices->end(); ++sv) {
     bool shared = false;
     VertexState s1(RecoVertex::convertPos(sv->position()), RecoVertex::convertError(sv->error()));
     for (typename Product::iterator sv2 = recoVertices->begin(); sv2 != recoVertices->end(); ++sv2) {
       VertexState s2(RecoVertex::convertPos(sv2->position()), RecoVertex::convertError(sv2->error()));
       double fr = vertexTools::computeSharedTracks(*sv2, *sv);
       //        std::cout << sv2-recoVertices->begin() << " vs " << sv-recoVertices->begin() << " : " << fr << " "  <<  computeSharedTracks(*sv, *sv2) << " sig " << dist.distance(*sv,*sv2).significance() << std::endl;
       //      std::cout << (fr > maxFraction) << " && " << (dist.distance(*sv,*sv2).significance() < 2)  <<  " && " <<  (sv-sv2!=0)  << " && " <<  (fr >= computeSharedTracks(*sv2, *sv))  << std::endl;
       if (fr > maxFraction && dist.distance(s1, s2).significance() < minSignificance && sv - sv2 != 0 &&
           fr >= vertexTools::computeSharedTracks(*sv, *sv2)) {
         shared = true;
         // std::cout << "shared " << sv-recoVertices->begin() << " and "  << sv2-recoVertices->begin() << " fractions: " << fr << " , "  << computeSharedTracks(*sv2, *sv) << " sig: " <<  dist.distance(*sv,*sv2).significance() <<  std::endl;
       }
     }
     if (shared) {
       sv = recoVertices->erase(sv) - 1;
     }
     //    std::cout << "it = " <<  sv-recoVertices->begin() << " new size is: " << recoVertices->size() <<   std::endl;
   }
 
   event.put(std::move(recoVertices));
 }