StandAloneMuonRefitter Class Reference

#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

RefitResult

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

Definition at line 23 of file StandAloneMuonRefitter.h.

Constructor & Destructor Documentation

StandAloneMuonRefitter()

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

Constructor.

Definition at line 22 of file StandAloneMuonRefitter.cc.

    : 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");
}

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

~StandAloneMuonRefitter()

StandAloneMuonRefitter::~StandAloneMuonRefitter ( )

virtual

Destructor.

Definition at line 34 of file StandAloneMuonRefitter.cc.

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

References LogDebug.

Member Function Documentation

refit()

StandAloneMuonRefitter::RefitResult StandAloneMuonRefitter::refit

(

const Trajectory &

trajectory

)

Definition at line 63 of file StandAloneMuonRefitter.cc.

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

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

Referenced by StandAloneMuonTrajectoryBuilder::trajectories().

singleRefit()

StandAloneMuonRefitter::RefitResult StandAloneMuonRefitter::singleRefit

(

const Trajectory &

trajectory

)

Refit.

Definition at line 41 of file StandAloneMuonRefitter.cc.

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

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

Referenced by refit().

Member Data Documentation

errorRescale

double StandAloneMuonRefitter::errorRescale

private

Definition at line 46 of file StandAloneMuonRefitter.h.

Referenced by singleRefit(), and StandAloneMuonRefitter().

isForceAllIterations

bool StandAloneMuonRefitter::isForceAllIterations

private

Definition at line 44 of file StandAloneMuonRefitter.h.

Referenced by refit(), and StandAloneMuonRefitter().

theFitter

edm::ESHandle<TrajectoryFitter> StandAloneMuonRefitter::theFitter

private

Definition at line 41 of file StandAloneMuonRefitter.h.

Referenced by singleRefit().

theFitterName

std::string StandAloneMuonRefitter::theFitterName

private

Definition at line 42 of file StandAloneMuonRefitter.h.

Referenced by singleRefit(), and StandAloneMuonRefitter().

theMaxFractionOfLostHits

double StandAloneMuonRefitter::theMaxFractionOfLostHits

private

Definition at line 45 of file StandAloneMuonRefitter.h.

Referenced by refit(), and StandAloneMuonRefitter().

theNumberOfIterations

unsigned int StandAloneMuonRefitter::theNumberOfIterations

private

Definition at line 43 of file StandAloneMuonRefitter.h.

Referenced by refit(), and StandAloneMuonRefitter().

theService

const MuonServiceProxy* StandAloneMuonRefitter::theService

private

Definition at line 40 of file StandAloneMuonRefitter.h.

Referenced by singleRefit().