AdaptiveGsfVertexFitter.cc

Go to the documentation of this file.

#include "RecoVertex/GaussianSumVertexFit/interface/AdaptiveGsfVertexFitter.h"
#include "RecoVertex/VertexTools/interface/DummyVertexSmoother.h"
#include "RecoVertex/GaussianSumVertexFit/interface/GsfVertexSmoother.h"
#include "RecoVertex/GaussianSumVertexFit/interface/MultiPerigeeLTSFactory.h"
#include "RecoVertex/GaussianSumVertexFit/interface/GsfVertexTrackCompatibilityEstimator.h"
#include "RecoVertex/VertexTools/interface/GeometricAnnealing.h"



AdaptiveGsfVertexFitter::AdaptiveGsfVertexFitter(const edm::ParameterSet& pSet,
    const LinearizationPointFinder & linP )
{

  bool limitComponents_ = pSet.getParameter<bool>("limitComponents");

  DeepCopyPointerByClone<GsfVertexMerger> theMerger;

  if (limitComponents_) {
    theMerger = new GsfVertexMerger(pSet.getParameter<edm::ParameterSet>("GsfMergerParameters"));
  }

  theFitter = new AdaptiveVertexFitter(
      GeometricAnnealing(),
      linP,
      GsfVertexUpdator(limitComponents_, &*theMerger),
      GsfVertexTrackCompatibilityEstimator(),
      GsfVertexSmoother(limitComponents_, &*theMerger),
      MultiPerigeeLTSFactory() );
  theFitter->gsfIntermediarySmoothing(true);

  theFitter->setParameters ( pSet.getParameter<double>("maxshift"),
                    pSet.getParameter<double>("maxlpshift"),
                    pSet.getParameter<int>("maxstep"),
                    pSet.getParameter<double>("weightthreshold") );

}

AdaptiveGsfVertexFitter::AdaptiveGsfVertexFitter(const AdaptiveGsfVertexFitter & original)
{
  theFitter = original.theFitter->clone();
}

AdaptiveGsfVertexFitter::~AdaptiveGsfVertexFitter()
{
  delete theFitter;
}