---

StandAloneMuonRefitter.cc

Go to the documentation of this file.

#include "RecoMuon/StandAloneTrackFinder/interface/StandAloneMuonRefitter.h"

#include "RecoMuon/TrackingTools/interface/MuonServiceProxy.h"

#include "TrackingTools/PatternTools/interface/TrajectoryFitter.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "TrackingTools/Records/interface/TrackingComponentsRecord.h"
#include "TrackingTools/TrackFitters/interface/TrajectoryStateWithArbitraryError.h"

using namespace edm;
using namespace std;

StandAloneMuonRefitter::StandAloneMuonRefitter(const ParameterSet& par, const MuonServiceProxy* service):theService(service){
  theFitterName = par.getParameter<string>("FitterName");
  theTEMPORARYoption = par.getParameter<int>("Option");
}

StandAloneMuonRefitter::~StandAloneMuonRefitter(){

}

  // Operations

StandAloneMuonRefitter::RefitResult StandAloneMuonRefitter::refit(const Trajectory& trajectory){
  
  theService->eventSetup().get<TrackingComponentsRecord>().get(theFitterName, theFitter);

  vector<Trajectory> refitted;

  if(theTEMPORARYoption == 1)
    refitted = theFitter->fit(trajectory);
  
  else if(theTEMPORARYoption == 2 || theTEMPORARYoption == 3){
    TrajectoryMeasurement lastTM = trajectory.lastMeasurement();
    TrajectoryStateOnSurface firstTsos = lastTM.updatedState();
    if (theTEMPORARYoption == 3) 
      firstTsos = TrajectoryStateWithArbitraryError()(lastTM.updatedState());
    TransientTrackingRecHit::ConstRecHitContainer trajRH = trajectory.recHits();
    reverse(trajRH.begin(),trajRH.end());

    vector<Trajectory> refitted = theFitter->fit(trajectory.seed(), trajRH, firstTsos);
  }
  
  if(!refitted.empty()) return RefitResult(true,refitted.front());
  else return RefitResult(false,trajectory);
}


// {
//   TransientTrackingRecHit::ConstRecHitContainer trajRH = trajectory.recHits();
//   for()

// }

---