#include <StandAloneMuonRefitter.h>

Public Types
typedef std::pair< bool, Trajectory >	RefitResult

Public Member Functions
RefitResult	refit (const Trajectory &)

RefitResult	singleRefit (const Trajectory &)
	Refit. More...

	StandAloneMuonRefitter (const edm::ParameterSet &par, const MuonServiceProxy *service)
	Constructor. More...

virtual	~StandAloneMuonRefitter ()
	Destructor. More...

Private Attributes
double	errorRescale

bool	isForceAllIterations

edm::ESHandle< TrajectoryFitter >	theFitter

std::string	theFitterName

double	theMaxFractionOfLostHits

unsigned int	theNumberOfIterations

const MuonServiceProxy *	theService

Detailed Description

Class ti interface the muon system rechits with the standard KF tools.

Authors: R. Bellan - INFN Torino ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch D. Trocino - INFN Torino danie.nosp@m.le.t.nosp@m.rocin.nosp@m.o@to.nosp@m..infn.nosp@m..it

Class ti interface the muon system rechits with the standard KF tools.

Authors: R. Bellan - INFN Torino ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch, D. Trocino - INFN Torino danie.nosp@m.le.t.nosp@m.rocin.nosp@m.o@to.nosp@m..infn.nosp@m..it

Definition at line 21 of file StandAloneMuonRefitter.h.

Member Typedef Documentation

typedef std::pair<bool, Trajectory> StandAloneMuonRefitter::RefitResult

Definition at line 23 of file StandAloneMuonRefitter.h.

Constructor & Destructor Documentation

StandAloneMuonRefitter::StandAloneMuonRefitter	(	const edm::ParameterSet &	par,
		const MuonServiceProxy *	service
	)

Constructor.

Definition at line 22 of file StandAloneMuonRefitter.cc.

References errorRescale, edm::ParameterSet::getParameter(), isForceAllIterations, LogDebug, theFitterName, theMaxFractionOfLostHits, and theNumberOfIterations.

     : theService(service) {
   LogDebug("Muon|RecoMuon|StandAloneMuonRefitter") << "Constructor called." << endl;
 
   theFitterName = par.getParameter<string>("FitterName");
   theNumberOfIterations = par.getParameter<unsigned int>("NumberOfIterations");
   isForceAllIterations = par.getParameter<bool>("ForceAllIterations");
   theMaxFractionOfLostHits = par.getParameter<double>("MaxFractionOfLostHits");
   errorRescale = par.getParameter<double>("RescaleError");
 }

StandAloneMuonRefitter::~StandAloneMuonRefitter ( )

virtual

Destructor.

Definition at line 34 of file StandAloneMuonRefitter.cc.

References LogDebug.

                                                 {
   LogDebug("Muon|RecoMuon|StandAloneMuonRefitter") << "Destructor called." << endl;
 }

Member Function Documentation

StandAloneMuonRefitter::RefitResult StandAloneMuonRefitter::refit ( const Trajectory & trajectory )

Definition at line 63 of file StandAloneMuonRefitter.cc.

References isForceAllIterations, dqmiolumiharvest::j, LogDebug, Trajectory::recHits(), singleRefit(), theMaxFractionOfLostHits, and theNumberOfIterations.

Referenced by StandAloneMuonTrajectoryBuilder::trajectories().

                                                                                             {
   LogDebug("Muon|RecoMuon|StandAloneMuonRefitter") << "---------------------------------" << endl;
   LogDebug("Muon|RecoMuon|StandAloneMuonRefitter") << "Starting refitting loop:" << endl;
 
   unsigned int nSuccess = 0;
   unsigned int nOrigHits = trajectory.recHits().size();
   Trajectory lastFitted = trajectory;
   bool allIter = true;
   bool enoughRH = true;
 
   for (unsigned int j = 0; j < theNumberOfIterations; ++j) {
     StandAloneMuonRefitter::RefitResult singleRefitResult = singleRefit(lastFitted);
     lastFitted = singleRefitResult.second;
     unsigned int nLastHits = lastFitted.recHits().size();
 
     if (!singleRefitResult.first) {
       allIter = false;
       LogDebug("Muon|RecoMuon|StandAloneMuonRefitter") << "  refit n. " << nSuccess + 1 << ": failed" << endl;
       break;
     }
 
     double lostFract = 1 - double(nLastHits) / nOrigHits;
     if (lostFract > theMaxFractionOfLostHits) {
       enoughRH = false;
       LogDebug("Muon|RecoMuon|StandAloneMuonRefitter") << "  refit n. " << nSuccess + 1 << ": too many RH lost" << endl;
       LogDebug("Muon|RecoMuon|StandAloneMuonRefitter")
           << "     Survived RecHits: " << nLastHits << "/" << nOrigHits << endl;
       break;
     }
 
     nSuccess++;
     LogDebug("Muon|RecoMuon|StandAloneMuonRefitter") << "  refit n. " << nSuccess << ": OK" << endl;
     LogDebug("Muon|RecoMuon|StandAloneMuonRefitter")
         << "     Survived RecHits: " << nLastHits << "/" << nOrigHits << endl;
 
   }  // end for
 
   LogDebug("Muon|RecoMuon|StandAloneMuonRefitter") << nSuccess << " successful refits!" << endl;
 
   // if isForceAllIterations==true  =>   3 successful refits: (true, refitted trajectory)
   //                                    <3 successful refits: (false, original trajectory)
   // if isForceAllIterations==false =>  >0 successful refits: (true, last refitted trajectory)
   //                                     0 successful refits: (false, original trajectory)
   if (!enoughRH)
     return RefitResult(false, trajectory);
   else if (isForceAllIterations)
     return allIter ? RefitResult(allIter, lastFitted) : RefitResult(allIter, trajectory);
   else
     return nSuccess == 0 ? RefitResult(false, trajectory) : RefitResult(true, lastFitted);
 }

StandAloneMuonRefitter::RefitResult StandAloneMuonRefitter::singleRefit ( const Trajectory & trajectory )

Refit.

Definition at line 41 of file StandAloneMuonRefitter.cc.

References errorRescale, TrajectoryFitter::fit(), Trajectory::lastMeasurement(), Trajectory::recHits(), TrajectoryStateOnSurface::rescaleError(), groupFilesInBlocks::reverse, Trajectory::seed(), theFitter, theFitterName, theService, and TrajectoryMeasurement::updatedState().

Referenced by refit().

                                                                                                   {
   theService->eventSetup().get<TrajectoryFitter::Record>().get(theFitterName, theFitter);
 
   vector<Trajectory> refitted;
 
   const TrajectoryMeasurement& lastTM = trajectory.lastMeasurement();
 
   TrajectoryStateOnSurface firstTsos(lastTM.updatedState());
 
   // Rescale errors before refit, not to bias the result
   firstTsos.rescaleError(errorRescale);
 
   TransientTrackingRecHit::ConstRecHitContainer trajRH = trajectory.recHits();
   reverse(trajRH.begin(), trajRH.end());
   refitted = theFitter->fit(trajectory.seed(), trajRH, firstTsos);
 
   if (!refitted.empty())
     return RefitResult(true, refitted.front());
   else
     return RefitResult(false, trajectory);
 }