MuonKinkFinder Class Reference

#include <MuonKinkFinder.h>

Public Member Functions
bool	fillTrkKink (reco::MuonQuality &quality, const reco::Track &track) const
bool	fillTrkKink (reco::MuonQuality &quality, const Trajectory &trajectory) const
void	init (const edm::EventSetup &iSetup)
	MuonKinkFinder (const edm::ParameterSet &iConfig, edm::ConsumesCollector &iC)
	~MuonKinkFinder ()

Private Member Functions
void	cropAndInvert (AlgebraicSymMatrix55 &cov) const
double	getChi2 (const TrajectoryStateOnSurface &start, const TrajectoryStateOnSurface &other) const

Private Attributes
bool	diagonalOnly_
	use only on-diagonal terms of the covariance matrices More...
TrackTransformer	refitter_
	Track Transformer. More...
bool	usePosition_
	if true, use full 5x5 track state; if false, use only the track direction More...

Detailed Description

Definition at line 8 of file MuonKinkFinder.h.

Constructor & Destructor Documentation

MuonKinkFinder()

MuonKinkFinder::MuonKinkFinder	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &	iC
	)

Definition at line 5 of file MuonKinkFinder.cc.

    : diagonalOnly_(iConfig.getParameter<bool>("diagonalOnly")),
      usePosition_(iConfig.getParameter<bool>("usePosition")),
      refitter_(iConfig, iC) {}

~MuonKinkFinder()

MuonKinkFinder::~MuonKinkFinder

(

)

Definition at line 10 of file MuonKinkFinder.cc.

{}

Member Function Documentation

cropAndInvert()

void MuonKinkFinder::cropAndInvert ( AlgebraicSymMatrix55 &  cov ) const

private

Definition at line 60 of file MuonKinkFinder.cc.

                                                                  {
  if (usePosition_) {
    if (diagonalOnly_) {
      for (size_t i = 0; i < 5; ++i) {
        for (size_t j = i + 1; j < 5; ++j) {
          cov(i, j) = 0;
        }
      }
    }
    cov.Invert();
  } else {
    // get 3x3 covariance
    AlgebraicSymMatrix33 momCov = cov.Sub<AlgebraicSymMatrix33>(0, 0);  // get 3x3 matrix
    if (diagonalOnly_) {
      momCov(0, 1) = 0;
      momCov(0, 2) = 0;
      momCov(1, 2) = 0;
    }
    // invert
    momCov.Invert();
    // place it
    cov.Place_at(momCov, 0, 0);
    // zero the rest
    for (size_t i = 3; i < 5; ++i) {
      for (size_t j = i; j < 5; ++j) {
        cov(i, j) = 0;
      }
    }
  }
}

References diagonalOnly_, mps_fire::i, dqmiolumiharvest::j, and usePosition_.

Referenced by getChi2().

fillTrkKink() [1/2]

bool MuonKinkFinder::fillTrkKink	(	reco::MuonQuality &	quality,
		const reco::Track &	track
	)		const

Definition at line 14 of file MuonKinkFinder.cc.

                                                                                       {
  std::vector<Trajectory> traj = refitter_.transform(track);
  if (traj.size() != 1) {
    quality.trkKink = 999;
    quality.tkKink_position = math::XYZPoint(0, 0, 0);
    return false;
  }
  return fillTrkKink(quality, traj.front());
}

References fillTrkKink(), quality, refitter_, HLT_FULL_cff::track, and TrackTransformer::transform().

fillTrkKink() [2/2]

bool MuonKinkFinder::fillTrkKink	(	reco::MuonQuality &	quality,
		const Trajectory &	trajectory
	)		const

Definition at line 24 of file MuonKinkFinder.cc.

                                                                                             {
  const std::vector<TrajectoryMeasurement> &tms = trajectory.measurements();
  quality.trkKink = -1.0;
  quality.tkKink_position = math::XYZPoint(0, 0, 0);
  bool found = false;
  for (int itm = 3, nm = tms.size() - 3; itm < nm; ++itm) {
    TrajectoryStateOnSurface pre = tms[itm].forwardPredictedState();
    TrajectoryStateOnSurface post = tms[itm].backwardPredictedState();
    if (!pre.isValid() || !post.isValid())
      continue;
    found = true;
    double c2f = getChi2(pre, post);
    if (c2f > quality.trkKink) {
      quality.trkKink = c2f;
      GlobalPoint pos = (tms[itm].updatedState().isValid() ? tms[itm].updatedState() : pre).globalPosition();
      quality.tkKink_position = math::XYZPoint(pos.x(), pos.y(), pos.z());
    }
  }
  if (!found)
    quality.trkKink = 999;
  return found;
}

References newFWLiteAna::found, getChi2(), TrajectoryStateOnSurface::isValid(), Trajectory::measurements(), and quality.

Referenced by fillTrkKink().

getChi2()

double MuonKinkFinder::getChi2 ( const TrajectoryStateOnSurface &  start, const TrajectoryStateOnSurface &  other  ) const

private

Definition at line 47 of file MuonKinkFinder.cc.

                                                                                                                 {
  if (!start.hasError() && !other.hasError())
    throw cms::Exception("LogicError") << "At least one of the two states must have errors to make chi2s.\n";
  AlgebraicSymMatrix55 cov;
  if (start.hasError())
    cov += start.localError().matrix();
  if (other.hasError())
    cov += other.localError().matrix();
  cropAndInvert(cov);
  AlgebraicVector5 diff(start.localParameters().mixedFormatVector() - other.localParameters().mixedFormatVector());
  return ROOT::Math::Similarity(diff, cov);
}

References cropAndInvert(), change_name::diff, and trackingPlots::other.

Referenced by fillTrkKink().

init()

void MuonKinkFinder::init

(

const edm::EventSetup &

iSetup

)

Definition at line 12 of file MuonKinkFinder.cc.

{ refitter_.setServices(iSetup); }

References refitter_, and TrackTransformer::setServices().

Member Data Documentation

diagonalOnly_

bool MuonKinkFinder::diagonalOnly_

private

use only on-diagonal terms of the covariance matrices

Definition at line 24 of file MuonKinkFinder.h.

Referenced by cropAndInvert().

refitter_

TrackTransformer MuonKinkFinder::refitter_

private

Track Transformer.

Definition at line 29 of file MuonKinkFinder.h.

Referenced by fillTrkKink(), and init().

usePosition_

bool MuonKinkFinder::usePosition_

private

if true, use full 5x5 track state; if false, use only the track direction

Definition at line 26 of file MuonKinkFinder.h.

Referenced by cropAndInvert().