StandAloneMuonSmoother.cc

Go to the documentation of this file.

 
#include "RecoMuon/StandAloneTrackFinder/interface/StandAloneMuonSmoother.h"
 
#include "RecoMuon/TrackingTools/interface/MuonServiceProxy.h"
 
#include "TrackingTools/KalmanUpdators/interface/Chi2MeasurementEstimator.h"
#include "TrackingTools/KalmanUpdators/interface/KFUpdator.h"
#include "TrackingTools/TrackFitters/interface/KFTrajectorySmoother.h"
 
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
using namespace edm;
using namespace std;
 
StandAloneMuonSmoother::StandAloneMuonSmoother(const ParameterSet& par, const MuonServiceProxy* service)
    : theService(service) {
  // The max allowed chi2 to accept a rechit in the fit
  theMaxChi2 = par.getParameter<double>("MaxChi2");
 
  // The errors of the trajectory state are multiplied by nSigma
  // to define acceptance of BoundPlane and maximalLocalDisplacement
  theNSigma = par.getParameter<double>("NumberOfSigma");  // default = 3.
 
  // The estimator: makes the decision wheter a measure is good or not
  // it isn't used by the updator which does the real fit. In fact, in principle,
  // a looser request onto the measure set can be requested
  // (w.r.t. the request on the accept/reject measure in the fit)
  theEstimator = new Chi2MeasurementEstimator(theMaxChi2, theNSigma);
 
  theErrorRescaling = par.getParameter<double>("ErrorRescalingFactor");
 
  thePropagatorName = par.getParameter<string>("Propagator");
 
  theUpdator = new KFUpdator();
 
  // The Kalman smoother
  theSmoother = nullptr;
}
 
StandAloneMuonSmoother::~StandAloneMuonSmoother() {
  if (theEstimator)
    delete theEstimator;
  if (theUpdator)
    delete theUpdator;
  if (theSmoother)
    delete theSmoother;
}
 
const Propagator* StandAloneMuonSmoother::propagator() const { return &*theService->propagator(thePropagatorName); }
 
void StandAloneMuonSmoother::renewTheSmoother() {
  if (theService->isTrackingComponentsRecordChanged()) {
    if (theSmoother)
      delete theSmoother;
    theSmoother = new KFTrajectorySmoother(propagator(), updator(), estimator());
  }
  if (!theSmoother)
    theSmoother = new KFTrajectorySmoother(propagator(), updator(), estimator());
}
 
StandAloneMuonSmoother::SmoothingResult StandAloneMuonSmoother::smooth(const Trajectory& inputTrajectory) {
  const string metname = "Muon|RecoMuon|StandAloneMuonSmoother";
 
  renewTheSmoother();
 
  vector<Trajectory> trajectoriesSM = smoother()->trajectories(inputTrajectory);
 
  if (trajectoriesSM.empty()) {
    LogTrace(metname) << "No Track smoothed!";
    return SmoothingResult(false, inputTrajectory);
  }
 
  Trajectory smoothed = trajectoriesSM.front();
 
  return SmoothingResult(true, smoothed);
}