reco::PositiveSideGhostTrackFitter Class Reference

#include <PositiveSideGhostTrackFitter.h>

Inheritance diagram for reco::PositiveSideGhostTrackFitter:

Public Member Functions
GhostTrackPrediction	fit (const GhostTrackFitter::PredictionUpdater &updater, const GhostTrackPrediction &prior, std::vector< GhostTrackState > &states, double &ndof, double &chi2) override
	PositiveSideGhostTrackFitter (const GlobalPoint &origin, const GhostTrackFitter::FitterImpl &actualFitter)
	PositiveSideGhostTrackFitter (const PositiveSideGhostTrackFitter &orig)
	~PositiveSideGhostTrackFitter () override
Public Member Functions inherited from reco::GhostTrackFitter::FitterImpl
virtual	~FitterImpl ()

Private Member Functions
FitterImpl *	clone () const override

Private Attributes
std::unique_ptr < GhostTrackFitter::FitterImpl >	actualFitter_
GlobalPoint	origin_

Detailed Description

Definition at line 16 of file PositiveSideGhostTrackFitter.h.

Constructor & Destructor Documentation

reco::PositiveSideGhostTrackFitter::PositiveSideGhostTrackFitter ( const GlobalPoint &  origin, const GhostTrackFitter::FitterImpl &  actualFitter  )

inline

Definition at line 18 of file PositiveSideGhostTrackFitter.h.

Referenced by clone().

        : origin_(origin), actualFitter_(actualFitter.clone()) {}

reco::PositiveSideGhostTrackFitter::~PositiveSideGhostTrackFitter ( )

inlineoverride

Definition at line 20 of file PositiveSideGhostTrackFitter.h.

{}

reco::PositiveSideGhostTrackFitter::PositiveSideGhostTrackFitter ( const PositiveSideGhostTrackFitter &  orig )

inline

Definition at line 22 of file PositiveSideGhostTrackFitter.h.

        : origin_(orig.origin_), actualFitter_(orig.actualFitter_->clone()) {}

Member Function Documentation

FitterImpl* reco::PositiveSideGhostTrackFitter::clone ( void  ) const

inlineoverrideprivatevirtual

Implements reco::GhostTrackFitter::FitterImpl.

Definition at line 32 of file PositiveSideGhostTrackFitter.h.

References PositiveSideGhostTrackFitter().

{ return new PositiveSideGhostTrackFitter(*this); }

GhostTrackPrediction PositiveSideGhostTrackFitter::fit ( const GhostTrackFitter::PredictionUpdater &  updater, const GhostTrackPrediction &  prior, std::vector< GhostTrackState > &  states, double &  ndof, double &  chi2  )

overridevirtual

Implements reco::GhostTrackFitter::FitterImpl.

Definition at line 15 of file PositiveSideGhostTrackFitter.cc.

References actualFitter_, HLT_FULL_cff::chi2, reco::GhostTrackFitter::PredictionUpdater::contribution(), fileCollector::done, funct::exp(), mps_fire::i, reco::GhostTrackState::isValid(), reco::GhostTrackState::lambda(), reco::GhostTrackPrediction::lambda(), reco::GhostTrackState::linearize(), SiStripPI::min, ndof, origin_, reco::GhostTrackPrediction::rho(), reco::GhostTrackState::setWeight(), histoStyle::weight, and reco::GhostTrackState::weight().

Referenced by trackingPlots.Iteration::modules().

                                                                     {
  double rho = prior.rho();
  for (unsigned int i = 0; i < states.size(); i++) {
    GhostTrackState &state = states[i];
    state.linearize(prior, true, .5 / rho);
  }

  GhostTrackPrediction pred = actualFitter_->fit(updater, prior, states, ndof, chi2);

  double origin = pred.lambda(origin_);
  bool done = true;
  for (unsigned int i = 0; i < states.size(); i++) {
    GhostTrackState &state = states[i];
    double lambda = state.lambda();
    if (lambda < origin && (origin - lambda) < 3.5) {
      GhostTrackState testState = state;
      testState.linearize(pred, 2. * origin - lambda);
      double ndof, chi2;

      if (testState.isValid()) {
        updater.contribution(prior, testState, ndof, chi2, true);
        if (ndof > 0. && chi2 < 10.) {
          state = testState;
          if (state.weight() != 1.)
            state.setWeight(3.);
          done = false;
        }
      } else {
        edm::LogError("InvalidGhostTrackState") << "Invalid GhostTrackState encountered!";
      }
    }
  }

  if (!done) {
    for (unsigned int i = 0; i < states.size(); i++) {
      GhostTrackState &state = states[i];
      double lambda = state.lambda();
      if (state.weight() != 1. && lambda < origin) {
        double weight = std::exp(10. * (origin - lambda) - .1);
        state.setWeight(std::min(state.weight(), weight));
      }
    }

    pred = actualFitter_->fit(updater, prior, states, ndof, chi2);
  }

  return pred;
}

Member Data Documentation

std::unique_ptr<GhostTrackFitter::FitterImpl> reco::PositiveSideGhostTrackFitter::actualFitter_

private

Definition at line 35 of file PositiveSideGhostTrackFitter.h.

Referenced by fit().

GlobalPoint reco::PositiveSideGhostTrackFitter::origin_

private

Definition at line 34 of file PositiveSideGhostTrackFitter.h.

Referenced by fit().