#include <GsfVertexMerger.h>

Public Member Functions
GsfVertexMerger *	clone () const

	GsfVertexMerger (const edm::ParameterSet &pSet)

CachingVertex< 5 >	merge (const CachingVertex< 5 > &vertex) const

VertexState	merge (const VertexState &vertex) const

	~GsfVertexMerger ()

Private Attributes
unsigned int	maxComponents

DeepCopyPointerByClone< MultiGaussianStateMerger< 3 > >	merger

Detailed Description

Definition at line 15 of file GsfVertexMerger.h.

Constructor & Destructor Documentation

◆ GsfVertexMerger()

GsfVertexMerger::GsfVertexMerger ( const edm::ParameterSet & pSet )

Definition at line 13 of file GsfVertexMerger.cc.

References HLT_2023v12_cff::distance, edm::ParameterSet::getParameter(), maxComponents, merger, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by clone().

                                                             {
   maxComponents = pSet.getParameter<int>("maxNbrComponents");
   std::string mergerName = pSet.getParameter<std::string>("merger");
   std::string distanceName = pSet.getParameter<std::string>("distance");
 
   if (mergerName == "CloseComponentsMerger") {
     DistanceBetweenComponents<3>* distance;
 
     if (distanceName == "KullbackLeiblerDistance")
       distance = new KullbackLeiblerDistance<3>();
     //     else if ( distanceName=="MahalanobisDistance" )
     //       distance = new MahalanobisDistance();
     else
       throw VertexException("GsfVertexMerger: Distance type " + distanceName +
                             " unknown. Check distance parameter in GsfMergerParameters PSet");
 
     merger = new CloseComponentsMerger<3>(maxComponents, distance);
     delete distance;
   }
   //   else if ( mergerName=="LargestWeightsStateMerger" )
   //     merger = new LargestWeightsStateMerger(maxComponents);
   else
     throw VertexException("GsfVertexMerger: Merger type " + mergerName +
                           " unknown. Check merger parameter in GsfMergerParameters PSet");
 
   //   else if ( mergerName=="KeepingNonZeroWeightsMerger" )
   //     merger = new KeepingNonZeroWeightsMerger();
 
   //   std::string mergerName = "test";
   //   edm::ESHandle<MultiGaussianStateMerger> mergerProducer;
   //   iRecord.get(mergerName,mergerProducer);
   //   merger = (MultiGaussianStateMerger *) mergerProducer.product();
 }

◆ ~GsfVertexMerger()

GsfVertexMerger::~GsfVertexMerger ( )

inline

Definition at line 18 of file GsfVertexMerger.h.

18 {}

Member Function Documentation

◆ clone()

GsfVertexMerger* GsfVertexMerger::clone ( void ) const

inline

Definition at line 24 of file GsfVertexMerger.h.

References GsfVertexMerger().

Referenced by GsfVertexSmoother::GsfVertexSmoother(), and GsfVertexUpdator::GsfVertexUpdator().

24 { return new GsfVertexMerger(*this); }

GsfVertexMerger::GsfVertexMerger

GsfVertexMerger(const edm::ParameterSet &pSet)

Definition: GsfVertexMerger.cc:13

◆ merge() [1/2]

CachingVertex< 5 > GsfVertexMerger::merge ( const CachingVertex< 5 > & vertex ) const

Definition at line 47 of file GsfVertexMerger.cc.

References VertexState::components(), CachingVertex< N >::hasPrior(), maxComponents, CachingVertex< N >::priorError(), CachingVertex< N >::priorPosition(), CachingVertex< N >::totalChiSquared(), CachingVertex< N >::tracks(), CachingVertex< N >::vertexState(), VertexState::weight(), and VertexState::weightTimesPosition().

Referenced by GsfVertexUpdator::add(), and GsfVertexSmoother::smooth().

                                                                                {
   if (oldVertex.vertexState().components().size() <= maxComponents)
     return oldVertex;
 
   VertexState newVertexState = merge(oldVertex.vertexState());
 
   if (oldVertex.hasPrior()) {
     return CachingVertex<5>(oldVertex.priorPosition(),
                             oldVertex.priorError(),
                             newVertexState.weightTimesPosition(),
                             newVertexState.weight(),
                             oldVertex.tracks(),
                             oldVertex.totalChiSquared());
   } else {
     return CachingVertex<5>(newVertexState, oldVertex.tracks(), oldVertex.totalChiSquared());
   }
 }

◆ merge() [2/2]

VertexState GsfVertexMerger::merge ( const VertexState & vertex ) const

Definition at line 65 of file GsfVertexMerger.cc.

References VertexState::components(), maxComponents, MultiGaussianStateMerger< N >::merge(), merger, GaussianStateConversions::multiGaussianStateFromVertex(), and GaussianStateConversions::vertexFromMultiGaussianState().

                                                                      {
   using namespace GaussianStateConversions;
 
   if (oldVertex.components().size() <= maxComponents)
     return oldVertex;
 
   MultiGaussianState<3> multiGaussianState(multiGaussianStateFromVertex(oldVertex));
   MultiGaussianState<3> finalState(merger->merge(multiGaussianState));
   return vertexFromMultiGaussianState(finalState);
 }