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; }