reco::TrackGhostTrackState Class Reference

#include <TrackGhostTrackState.h>

Inheritance diagram for reco::TrackGhostTrackState:

Public Member Functions
CovarianceMatrix	cartesianCovariance () const override
GlobalError	cartesianError () const override
GlobalPoint	globalPosition () const override
bool	isValid () const override
bool	linearize (const GhostTrackPrediction &pred, bool initial, double lambda) override
bool	linearize (const GhostTrackPrediction &pred, double lambda) override
void	reset () override
const TransientTrack &	track () const
	TrackGhostTrackState (const TransientTrack &track)
const TrajectoryStateOnSurface &	tsos () const
Vertex	vertexStateOnGhostTrack (const GhostTrackPrediction &pred, bool withMeasurementError) const override
Vertex	vertexStateOnMeasurement (const GhostTrackPrediction &pred, bool withGhostTrackError) const override
Public Member Functions inherited from reco::BasicGhostTrackState
double	lambda () const
void	setWeight (double weight)
double	weight () const
virtual	~BasicGhostTrackState ()

Private Member Functions
pointer	clone () const override

Private Attributes
TransientTrack	track_
TrajectoryStateOnSurface	tsos_

Additional Inherited Members
Public Types inherited from reco::BasicGhostTrackState
using	BGTS = BasicGhostTrackState
typedef math::Error< 3 >::type	CovarianceMatrix
using	pointer = Proxy::pointer
using	Proxy = ProxyBase11< BGTS >
typedef std::pair< GlobalPoint, GlobalError >	Vertex
Static Public Member Functions inherited from reco::BasicGhostTrackState
template<typename T , typename... Args>
static std::shared_ptr< BGTS >	build (Args &&...args)
Protected Attributes inherited from reco::BasicGhostTrackState
double	lambda_ = 0
double	weight_ = 1.

Detailed Description

Definition at line 19 of file TrackGhostTrackState.h.

Constructor & Destructor Documentation

reco::TrackGhostTrackState::TrackGhostTrackState ( const TransientTrack &  track )

inline

Definition at line 21 of file TrackGhostTrackState.h.

: track_(track) {}

Member Function Documentation

CovarianceMatrix reco::TrackGhostTrackState::cartesianCovariance ( ) const

inlineoverridevirtual

Implements reco::BasicGhostTrackState.

Definition at line 30 of file TrackGhostTrackState.h.

References TrajectoryStateOnSurface::cartesianError(), CartesianTrajectoryError::matrix(), and tsos_.

                                                          {
      return tsos_.cartesianError().matrix().Sub<CovarianceMatrix>(0, 0);
    }

GlobalError reco::TrackGhostTrackState::cartesianError ( ) const

inlineoverridevirtual

Implements reco::BasicGhostTrackState.

Definition at line 29 of file TrackGhostTrackState.h.

References TrajectoryStateOnSurface::cartesianError(), CartesianTrajectoryError::position(), and tsos_.

{ return tsos_.cartesianError().position(); }

pointer reco::TrackGhostTrackState::clone ( void  ) const

inlineoverrideprivatevirtual

Implements reco::BasicGhostTrackState.

Definition at line 42 of file TrackGhostTrackState.h.

{ return build<TrackGhostTrackState>(*this); }

GlobalPoint reco::TrackGhostTrackState::globalPosition ( ) const

inlineoverridevirtual

Implements reco::BasicGhostTrackState.

Definition at line 28 of file TrackGhostTrackState.h.

References TrajectoryStateOnSurface::globalPosition(), and tsos_.

{ return tsos_.globalPosition(); }

bool reco::TrackGhostTrackState::isValid ( void  ) const

inlineoverridevirtual

Reimplemented from reco::BasicGhostTrackState.

Definition at line 26 of file TrackGhostTrackState.h.

References TrajectoryStateOnSurface::isValid(), and tsos_.

Referenced by ntupleDataFormat._Object::_checkIsValid(), linearize(), core.AutoHandle.AutoHandle::ReallyLoad(), vertexStateOnGhostTrack(), and vertexStateOnMeasurement().

{ return tsos_.isValid(); }

bool TrackGhostTrackState::linearize ( const GhostTrackPrediction &  pred, bool  initial, double  lambda  )

overridevirtual

Reimplemented from reco::BasicGhostTrackState.

Definition at line 43 of file TrackGhostTrackState.cc.

References reco::GhostTrackPrediction::direction(), reco::TransientTrack::field(), TrajectoryStateOnSurface::globalPosition(), reco::TransientTrack::impactPointState(), isValid(), reco::BasicGhostTrackState::lambda_, geometryCSVtoXML::line, reco::GhostTrackPrediction::origin(), reco::GhostTrackPrediction::rho2(), track_, and tsos_.

                                                                                                  {
  AnalyticalTrajectoryExtrapolatorToLine extrap(track_.field());

  GlobalPoint origin = pred.origin();
  GlobalVector direction = pred.direction();

  if (isValid() && !initial) {
    GlobalPoint pca = origin + lambda_ * direction;
    Line line(pca, direction);
    tsos_ = extrap.extrapolate(tsos_, line);
  } else {
    GlobalPoint pca = origin + lambda * direction;
    Line line(pca, direction);
    tsos_ = extrap.extrapolate(track_.impactPointState(), line);
  }

  if (!isValid())
    return false;

  lambda_ = (tsos_.globalPosition() - origin) * direction / pred.rho2();

  return true;
}

bool TrackGhostTrackState::linearize ( const GhostTrackPrediction &  pred, double  lambda  )

overridevirtual

Reimplemented from reco::BasicGhostTrackState.

Definition at line 67 of file TrackGhostTrackState.cc.

References reco::TransientTrack::field(), reco::TransientTrack::impactPointState(), isValid(), reco::BasicGhostTrackState::lambda(), reco::BasicGhostTrackState::lambda_, point, reco::GhostTrackPrediction::position(), track_, and tsos_.

                                                                                    {
  AnalyticalImpactPointExtrapolator extrap(track_.field());

  GlobalPoint point = pred.position(lambda);

  tsos_ = extrap.extrapolate(track_.impactPointState(), point);
  if (!isValid())
    return false;

  lambda_ = lambda;

  return true;
}

void reco::TrackGhostTrackState::reset ( void  )

inlineoverridevirtual

Reimplemented from reco::BasicGhostTrackState.

Definition at line 34 of file TrackGhostTrackState.h.

References tsos_.

{ tsos_ = TrajectoryStateOnSurface(); }

const TransientTrack& reco::TrackGhostTrackState::track ( ) const

inline

Definition at line 23 of file TrackGhostTrackState.h.

References track_.

Referenced by reco::GhostTrackState::track().

{ return track_; }

const TrajectoryStateOnSurface& reco::TrackGhostTrackState::tsos ( ) const

inline

Definition at line 24 of file TrackGhostTrackState.h.

References tsos_.

Referenced by reco::GhostTrackState::tsos().

{ return tsos_; }

BasicGhostTrackState::Vertex TrackGhostTrackState::vertexStateOnGhostTrack ( const GhostTrackPrediction &  pred, bool  withMeasurementError  ) const

overridevirtual

Implements reco::BasicGhostTrackState.

Definition at line 81 of file TrackGhostTrackState.cc.

References b, reco::GhostTrackPrediction::cartesianError(), TrajectoryStateOnSurface::cartesianError(), conv, ztail::d, reco::GhostTrackPrediction::direction(), relativeConstraints::error, g, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), isValid(), cmsLHEtoEOSManager::l, reco::BasicGhostTrackState::lambda(), PV3DBase< T, PVType, FrameType >::mag(), CartesianTrajectoryError::matrix(), reco::GhostTrackPrediction::origin(), reco::GhostTrackPrediction::rho2(), sqr(), mathSSE::sqrt(), and tsos_.

                                                                                                            {
  using namespace ROOT::Math;

  if (!isValid())
    return Vertex();

  GlobalPoint origin = pred.origin();
  GlobalVector direction = pred.direction();

  double rho2 = pred.rho2();
  double rho = std::sqrt(rho2);
  double lambda = (tsos_.globalPosition() - origin) * direction / rho2;
  GlobalPoint pos = origin + lambda * direction;

  GlobalVector momentum = tsos_.globalMomentum();
  double mom = momentum.mag();

  Vector3 b = conv(direction) / rho;
  Vector3 d = conv(momentum) / mom;
  double l = Dot(b, d);
  double g = 1. / (1. - sqr(l));

  Vector3 ca = conv(pos - tsos_.globalPosition());
  Vector3 bd = b - l * d;
  b *= g;

  Matrix33 pA = TensorProd(b, bd);
  Matrix33 pB = TensorProd(b, ca);

  Matrix36 jacobian;
  jacobian.Place_at(-pA + Matrix33(SMatrixIdentity()), 0, 0);
  jacobian.Place_at(pB / rho, 0, 3);
  Matrix3S error = Similarity(jacobian, pred.cartesianError(lambda));

  if (withMeasurementError) {
    jacobian.Place_at(pA, 0, 0);
    jacobian.Place_at(-pB / mom, 0, 3);
    error += Similarity(jacobian, tsos_.cartesianError().matrix());
  }

  return Vertex(pos, error);
}

BasicGhostTrackState::Vertex TrackGhostTrackState::vertexStateOnMeasurement ( const GhostTrackPrediction &  pred, bool  withGhostTrackError  ) const

overridevirtual

Implements reco::BasicGhostTrackState.

Definition at line 125 of file TrackGhostTrackState.cc.

References b, reco::GhostTrackPrediction::cartesianError(), TrajectoryStateOnSurface::cartesianError(), conv, ztail::d, reco::GhostTrackPrediction::direction(), relativeConstraints::error, g, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), isValid(), cmsLHEtoEOSManager::l, reco::BasicGhostTrackState::lambda(), PV3DBase< T, PVType, FrameType >::mag(), CartesianTrajectoryError::matrix(), reco::GhostTrackPrediction::origin(), reco::GhostTrackPrediction::rho2(), sqr(), mathSSE::sqrt(), and tsos_.

                                                                                                            {
  using namespace ROOT::Math;

  if (!isValid())
    return Vertex();

  GlobalPoint origin = pred.origin();
  GlobalVector direction = pred.direction();

  double rho2 = pred.rho2();
  double rho = std::sqrt(rho2);
  double lambda = (tsos_.globalPosition() - origin) * direction / rho2;
  GlobalPoint pos = origin + lambda * direction;

  GlobalVector momentum = tsos_.globalMomentum();
  double mom = momentum.mag();

  Vector3 b = conv(direction) / rho;
  Vector3 d = conv(momentum) / mom;
  double l = Dot(b, d);
  double g = 1. / (1. - sqr(l));

  Vector3 ca = conv(tsos_.globalPosition() - pos);
  Vector3 bd = l * b - d;
  d *= g;

  Matrix33 pC = TensorProd(d, bd);
  Matrix33 pD = TensorProd(d, ca);

  Matrix36 jacobian;
  jacobian.Place_at(pC + Matrix33(SMatrixIdentity()), 0, 0);
  jacobian.Place_at(pD / mom, 0, 3);
  Matrix3S error = Similarity(jacobian, tsos_.cartesianError().matrix());

  if (withGhostTrackError) {
    jacobian.Place_at(-pC, 0, 0);
    jacobian.Place_at(-pD / rho, 0, 3);
    error += Similarity(jacobian, pred.cartesianError(lambda));
  }

  return Vertex(tsos_.globalPosition(), error);
}

Member Data Documentation

TransientTrack reco::TrackGhostTrackState::track_

private

Definition at line 45 of file TrackGhostTrackState.h.

Referenced by linearize(), and track().

TrajectoryStateOnSurface reco::TrackGhostTrackState::tsos_

private

Definition at line 44 of file TrackGhostTrackState.h.

Referenced by cartesianCovariance(), cartesianError(), globalPosition(), isValid(), linearize(), reset(), tsos(), vertexStateOnGhostTrack(), and vertexStateOnMeasurement().