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 22 of file PositiveSideGhostTrackFitter.h.

{}

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

inline

Definition at line 24 of file PositiveSideGhostTrackFitter.h.

References vertices_cff::chi2, fit(), ndof, and fftjetcommon_cfi::prior.

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

Member Function Documentation

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

inlineoverrideprivatevirtual

Implements reco::GhostTrackFitter::FitterImpl.

Definition at line 36 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_, vertices_cff::chi2, reco::GhostTrackFitter::PredictionUpdater::contribution(), JetChargeProducer_cfi::exp, mps_fire::i, reco::GhostTrackState::isValid(), reco::GhostTrackState::lambda(), reco::GhostTrackPrediction::lambda(), reco::GhostTrackState::linearize(), min(), ndof, origin_, reco::GhostTrackPrediction::rho(), reco::GhostTrackState::setWeight(), reco::GhostTrackState::weight(), and mps_merge::weight.

Referenced by trackingPlots.Iteration::modules(), and PositiveSideGhostTrackFitter().

{
    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 40 of file PositiveSideGhostTrackFitter.h.

Referenced by fit().

GlobalPoint reco::PositiveSideGhostTrackFitter::origin_

private

Definition at line 39 of file PositiveSideGhostTrackFitter.h.

Referenced by fit().