reco::GhostTrackPrediction Class Reference

#include <GhostTrackPrediction.h>

Public Types
typedef ROOT::Math::SMatrix < double, 6, 6, ROOT::Math::MatRepSym < double, 6 > >	CartesianError
typedef ROOT::Math::SMatrix < double, 4, 4, ROOT::Math::MatRepSym < double, 4 > >	Error
typedef ROOT::Math::SVector < double, 4 >	Vector

Public Member Functions
CartesianError	cartesianError (double lambda=0.) const
double	cotTheta () const
const Error &	covariance () const
CurvilinearTrajectoryError	curvilinearError () const
const GlobalVector	direction () const
double	eta () const
FreeTrajectoryState	fts (const MagneticField *fieldProvider) const
	GhostTrackPrediction ()
	GhostTrackPrediction (const Vector &prediction, const Error &error)
	GhostTrackPrediction (const GlobalPoint &priorPosition, const GlobalError &priorError, const GlobalVector &direction, double coneRadius)
	GhostTrackPrediction (const GlobalPoint &priorPosition, const GlobalError &priorError, const GlobalVector &direction, const GlobalError &directionError)
	GhostTrackPrediction (const GlobalTrajectoryParameters &trajectory, const CurvilinearTrajectoryError &error)
	GhostTrackPrediction (const Track &track)
GlobalTrajectoryParameters	globalTrajectory (const MagneticField *fieldProvider) const
double	ip () const
double	lambda (const GlobalPoint &point) const
const GlobalPoint	origin () const
double	phi () const
GlobalPoint	position (double lambda=0.) const
GlobalError	positionError (double lambda=0.) const
const Vector &	prediction () const
double	rho () const
double	rho2 () const
double	sz () const
double	theta () const
Track	track (double ndof=0., double chi2=0.) const
double	z () const

Private Member Functions
void	init (const GlobalPoint &priorPosition, const GlobalError &priorError, const GlobalVector &direction, const GlobalError &directionError)

Private Attributes
Error	covariance_
Vector	prediction_

Detailed Description

Definition at line 21 of file GhostTrackPrediction.h.

Member Typedef Documentation

typedef ROOT::Math::SMatrix<double, 6, 6, ROOT::Math::MatRepSym<double, 6> > reco::GhostTrackPrediction::CartesianError

Definition at line 27 of file GhostTrackPrediction.h.

typedef ROOT::Math::SMatrix<double, 4, 4, ROOT::Math::MatRepSym<double, 4> > reco::GhostTrackPrediction::Error

Definition at line 26 of file GhostTrackPrediction.h.

typedef ROOT::Math::SVector<double, 4> reco::GhostTrackPrediction::Vector

Definition at line 25 of file GhostTrackPrediction.h.

Constructor & Destructor Documentation

reco::GhostTrackPrediction::GhostTrackPrediction ( )

inline

Definition at line 29 of file GhostTrackPrediction.h.

{}

reco::GhostTrackPrediction::GhostTrackPrediction ( const Vector &  prediction, const Error &  error  )

inline

Definition at line 30 of file GhostTrackPrediction.h.

: prediction_(prediction), covariance_(error) {}

GhostTrackPrediction::GhostTrackPrediction	(	const GlobalPoint &	priorPosition,
		const GlobalError &	priorError,
		const GlobalVector &	direction,
		double	coneRadius
	)

Definition at line 94 of file GhostTrackPrediction.cc.

References init(), AlCaHLTBitMon_ParallelJobs::p, perp(), alignCSCRings::r, diffTwoXMLs::r2, mathSSE::sqrt(), and submitPVValidationJobs::t.

                                                              {
  double dTheta = std::cosh((double)direction.eta()) * coneRadius;

  double r2 = direction.mag2();
  double r = std::sqrt(r2);
  double perp = direction.perp();
  double P = direction.x() / perp;
  double p = direction.y() / perp;
  double T = direction.z() / r;
  double t = perp / r;
  double h2 = dTheta * dTheta;
  double d2 = coneRadius * coneRadius;

  GlobalError cov(r2 * (T * T * P * P * h2 + t * t * p * p * d2),
                  r2 * p * P * (T * T * h2 - t * t * d2),
                  r2 * (T * T * p * p * h2 + t * t * P * P * d2),
                  -r2 * t * T * P * h2,
                  -r2 * t * T * p * h2,
                  r2 * t * t * h2);

  init(priorPosition, priorError, direction, cov);
}

reco::GhostTrackPrediction::GhostTrackPrediction ( const GlobalPoint &  priorPosition, const GlobalError &  priorError, const GlobalVector &  direction, const GlobalError &  directionError  )

inline

Definition at line 37 of file GhostTrackPrediction.h.

References init().

                                                            {
      init(priorPosition, priorError, direction, directionError);
    }

GhostTrackPrediction::GhostTrackPrediction	(	const GlobalTrajectoryParameters &	trajectory,
		const CurvilinearTrajectoryError &	error
	)

Definition at line 120 of file GhostTrackPrediction.cc.

    : prediction_(convert(trajectory)), covariance_(convert(prediction_, error)) {}

GhostTrackPrediction::GhostTrackPrediction

(

const Track &

track

)

Definition at line 124 of file GhostTrackPrediction.cc.

    : prediction_(convert(GlobalTrajectoryParameters(GlobalPoint(track.vx(), track.vy(), track.vz()),
                                                     GlobalVector(track.px(), track.py(), track.pz()),
                                                     0,
                                                     nullptr))),
      covariance_(convert(prediction_, track.covariance())) {}

Member Function Documentation

Matrix6S GhostTrackPrediction::cartesianError

(

double

lambda = 0.

)

const

Definition at line 146 of file GhostTrackPrediction.cc.

References funct::cos(), covariance(), ip(), lambda(), phi(), and funct::sin().

Referenced by reco::VertexGhostTrackState::vertexStateOnGhostTrack(), reco::TrackGhostTrackState::vertexStateOnGhostTrack(), and reco::TrackGhostTrackState::vertexStateOnMeasurement().

                                                                 {
  double x = std::cos(phi());
  double y = std::sin(phi());

  Matrix64 jacobian;
  jacobian(0, 1) = -y;
  jacobian(0, 3) = -y * lambda - x * ip();
  jacobian(1, 1) = x;
  jacobian(1, 3) = x * lambda - y * ip();
  jacobian(2, 0) = 1.;
  jacobian(2, 2) = lambda;
  jacobian(3, 3) = -y;
  jacobian(4, 3) = x;
  jacobian(5, 2) = 1.;

  return ROOT::Math::Similarity(jacobian, covariance());
}

double reco::GhostTrackPrediction::cotTheta ( ) const

inline

Definition at line 50 of file GhostTrackPrediction.h.

References prediction_.

Referenced by curvilinearError(), direction(), eta(), rho2(), and theta().

{ return prediction_[2]; }

const Error& reco::GhostTrackPrediction::covariance ( void  ) const

inline

Definition at line 60 of file GhostTrackPrediction.h.

References covariance_.

Referenced by cartesianError(), curvilinearError(), positionError(), and reco::KalmanGhostTrackUpdater::update().

{ return covariance_; }

CurvilinearTrajectoryError GhostTrackPrediction::curvilinearError

(

)

const

Definition at line 168 of file GhostTrackPrediction.cc.

References cotTheta(), covariance(), submitPVValidationJobs::err, rho2(), mathSSE::sqrt(), and z().

Referenced by fts(), and track().

                                                                        {
  double rho2I = 1. / rho2();
  double rhoI = std::sqrt(rho2I);

  Matrix54 jacobian;
  jacobian(1, 2) = rho2I;
  jacobian(2, 3) = 1.;
  jacobian(3, 1) = 1.;
  jacobian(4, 0) = rhoI;
  jacobian(4, 2) = -z() * rhoI * cotTheta() * rho2I;

  AlgebraicSymMatrix55 err = Similarity(jacobian, covariance());
  err(0, 0) = 1.;

  return CurvilinearTrajectoryError(err);
}

const GlobalVector reco::GhostTrackPrediction::direction ( ) const

inline

Definition at line 63 of file GhostTrackPrediction.h.

References funct::cos(), cotTheta(), phi(), and funct::sin().

Referenced by reco::GhostTrackState::axisDistance(), globalTrajectory(), lambda(), reco::GhostTrackState::lambdaError(), reco::TrackGhostTrackState::linearize(), position(), IPProducer< Container, Base, Helper >::produce(), track(), reco::VertexGhostTrackState::vertexStateOnGhostTrack(), reco::TrackGhostTrackState::vertexStateOnGhostTrack(), and reco::TrackGhostTrackState::vertexStateOnMeasurement().

{ return GlobalVector(std::cos(phi()), std::sin(phi()), cotTheta()); }

double reco::GhostTrackPrediction::eta ( void  ) const

inline

Definition at line 57 of file GhostTrackPrediction.h.

References cotTheta(), log, and rho().

Referenced by Particle.Particle::__str__(), Jet.Jet::jetID(), Jet.Jet::puJetId(), and reco::SequentialGhostTrackFitter::stable().

{ return -std::log(rho() - cotTheta()); }

FreeTrajectoryState GhostTrackPrediction::fts

(

const MagneticField *

fieldProvider

)

const

Definition at line 185 of file GhostTrackPrediction.cc.

References curvilinearError(), and globalTrajectory().

Referenced by transientGhostTrack().

                                                                                      {
  return FreeTrajectoryState(globalTrajectory(fieldProvider), curvilinearError());
}

GlobalTrajectoryParameters GhostTrackPrediction::globalTrajectory

(

const MagneticField *

fieldProvider

)

const

Definition at line 164 of file GhostTrackPrediction.cc.

References direction(), and origin().

Referenced by fts().

                                                                                                          {
  return GlobalTrajectoryParameters(origin(), direction(), 0, fieldProvider);
}

void GhostTrackPrediction::init ( const GlobalPoint &  priorPosition, const GlobalError &  priorError, const GlobalVector &  direction, const GlobalError &  directionError  )

private

Definition at line 58 of file GhostTrackPrediction.cc.

References covariance_, DeadROC_duringRun::dir, cmsLHEtoEOSManager::l, GlobalErrorBase< T, ErrorWeightType >::matrix(), perp2(), prediction_, mathSSE::sqrt(), and HLT_FULL_cff::tip.

Referenced by GhostTrackPrediction().

                                                                   {
  double perp2 = direction.perp2();
  GlobalVector dir = direction / std::sqrt(perp2);
  double tip = priorPosition.y() * dir.x() - priorPosition.x() * dir.y();
  double l = priorPosition.x() * dir.x() + priorPosition.y() * dir.y();

  prediction_[0] = priorPosition.z() - dir.z() * l;
  prediction_[1] = tip;
  prediction_[2] = dir.z();
  prediction_[3] = dir.phi();

  Matrix46 jacobian;
  jacobian(0, 0) = -dir.x() * dir.z();
  jacobian(1, 0) = -dir.y();
  jacobian(0, 1) = -dir.y() * dir.z();
  jacobian(1, 1) = dir.x();
  jacobian(0, 2) = 1.;
  jacobian(0, 3) = -dir.z() * priorPosition.x();
  jacobian(1, 3) = priorPosition.y();
  jacobian(3, 3) = -dir.y();
  jacobian(0, 4) = -dir.z() * priorPosition.y();
  jacobian(1, 4) = -priorPosition.x();
  jacobian(3, 4) = dir.x();
  jacobian(0, 5) = -l;
  jacobian(2, 5) = 1.;

  Matrix6S origCov;
  origCov.Place_at(priorError.matrix(), 0, 0);
  origCov.Place_at(directionError.matrix() / perp2, 3, 3);

  covariance_ = ROOT::Math::Similarity(jacobian, origCov);
}

double reco::GhostTrackPrediction::ip ( ) const

inline

Definition at line 49 of file GhostTrackPrediction.h.

References prediction_.

Referenced by cartesianError(), origin(), positionError(), and reco::SequentialGhostTrackFitter::stable().

{ return prediction_[1]; }

double reco::GhostTrackPrediction::lambda ( const GlobalPoint &  point ) const

inline

Definition at line 65 of file GhostTrackPrediction.h.

References direction(), origin(), and rho2().

Referenced by cartesianError(), reco::PositiveSideGhostTrackFitter::fit(), reco::GhostTrackFitter::fit(), reco::GhostTrack::GhostTrack(), position(), positionError(), IPProducer< Container, Base, Helper >::produce(), and reco::GhostTrackVertexFinder::reassignTracks().

{ return (point - origin()) * direction() / rho2(); }

const GlobalPoint reco::GhostTrackPrediction::origin ( ) const

inline

Definition at line 62 of file GhostTrackPrediction.h.

References funct::cos(), ip(), phi(), funct::sin(), and z().

Referenced by reco::GhostTrackState::axisDistance(), globalTrajectory(), lambda(), reco::TrackGhostTrackState::linearize(), position(), track(), reco::VertexGhostTrackState::vertexStateOnGhostTrack(), reco::TrackGhostTrackState::vertexStateOnGhostTrack(), and reco::TrackGhostTrackState::vertexStateOnMeasurement().

{ return GlobalPoint(-std::sin(phi()) * ip(), std::cos(phi()) * ip(), z()); }

double reco::GhostTrackPrediction::phi ( void  ) const

inline

Definition at line 51 of file GhostTrackPrediction.h.

References prediction_.

Referenced by Particle.Particle::__str__(), cartesianError(), direction(), origin(), ntupleDataFormat.Track::phiPull(), positionError(), and reco::SequentialGhostTrackFitter::stable().

{ return prediction_[3]; }

GlobalPoint reco::GhostTrackPrediction::position ( double  lambda = 0. ) const

inline

Definition at line 67 of file GhostTrackPrediction.h.

References direction(), lambda(), and origin().

Referenced by reco::TrackGhostTrackState::linearize(), IPProducer< Container, Base, Helper >::produce(), and vertexAtState().

{ return origin() + lambda * direction(); }

GlobalError GhostTrackPrediction::positionError

(

double

lambda = 0.

)

const

Definition at line 131 of file GhostTrackPrediction.cc.

References funct::cos(), covariance(), ip(), lambda(), phi(), and funct::sin().

Referenced by IPProducer< Container, Base, Helper >::produce(), and vertexAtState().

                                                                   {
  double x = std::cos(phi());
  double y = std::sin(phi());

  Matrix34 jacobian;
  jacobian(0, 1) = -y;
  jacobian(0, 3) = -y * lambda - x * ip();
  jacobian(1, 1) = x;
  jacobian(1, 3) = x * lambda - y * ip();
  jacobian(2, 0) = 1.;
  jacobian(2, 2) = lambda;

  return ROOT::Math::Similarity(jacobian, covariance());
}

const Vector& reco::GhostTrackPrediction::prediction ( ) const

inline

Definition at line 59 of file GhostTrackPrediction.h.

References prediction_.

Referenced by IPProducer< Container, Base, Helper >::produce(), and reco::KalmanGhostTrackUpdater::update().

{ return prediction_; }

double reco::GhostTrackPrediction::rho ( ) const

inline

Definition at line 54 of file GhostTrackPrediction.h.

References rho2(), and mathSSE::sqrt().

Referenced by Lepton.Lepton::absIsoFromEA(), Muon.Muon::absIsoWithFSR(), eta(), reco::PositiveSideGhostTrackFitter::fit(), and sz().

{ return std::sqrt(rho2()); }

double reco::GhostTrackPrediction::rho2 ( ) const

inline

Definition at line 53 of file GhostTrackPrediction.h.

References cotTheta().

Referenced by curvilinearError(), lambda(), reco::GhostTrackState::lambdaError(), reco::TrackGhostTrackState::linearize(), rho(), reco::VertexGhostTrackState::vertexStateOnGhostTrack(), reco::TrackGhostTrackState::vertexStateOnGhostTrack(), and reco::TrackGhostTrackState::vertexStateOnMeasurement().

{ return cotTheta() * cotTheta() + 1.; }

double reco::GhostTrackPrediction::sz ( ) const

inline

Definition at line 55 of file GhostTrackPrediction.h.

References rho(), and z().

Referenced by reco::SequentialGhostTrackFitter::stable().

{ return z() / rho(); }

double reco::GhostTrackPrediction::theta ( void  ) const

inline

Definition at line 56 of file GhostTrackPrediction.h.

References cotTheta(), and M_PI_2.

Referenced by Tau.Tau::zImpact().

{ return M_PI_2 - std::atan(cotTheta()); }

Track GhostTrackPrediction::track	(	double	ndof = 0.,
		double	chi2 = 0.
	)		const

Definition at line 189 of file GhostTrackPrediction.cc.

References curvilinearError(), DeadROC_duringRun::dir, direction(), origin(), point, and trackerHitRTTI::vector.

Referenced by reco::GhostTrack::operator Track().

                                                                {
  GlobalPoint origin = this->origin();
  GlobalVector dir = direction().unit();

  Track::Point point(origin.x(), origin.y(), origin.z());
  Track::Vector vector(dir.x(), dir.y(), dir.z());

  return Track(chi2, ndof, point, vector, 0, curvilinearError());
}

double reco::GhostTrackPrediction::z ( ) const

inline

Definition at line 48 of file GhostTrackPrediction.h.

References prediction_.

Referenced by curvilinearError(), origin(), geometryXMLparser.Alignable::pos(), ntupleDataFormat._HitObject::r3D(), and sz().

{ return prediction_[0]; }

Member Data Documentation

Error reco::GhostTrackPrediction::covariance_

private

Definition at line 86 of file GhostTrackPrediction.h.

Referenced by covariance(), and init().

Vector reco::GhostTrackPrediction::prediction_

private

Definition at line 85 of file GhostTrackPrediction.h.

Referenced by cotTheta(), init(), ip(), phi(), prediction(), and z().