#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 GlobalPoint &priorPosition, const GlobalError &priorError, const GlobalVector &direction, const GlobalError &directionError)

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

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

	GhostTrackPrediction (const Track &track)

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

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

◆ CartesianError

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

Definition at line 27 of file GhostTrackPrediction.h.

◆ Error

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

Definition at line 26 of file GhostTrackPrediction.h.

◆ Vector

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

Definition at line 25 of file GhostTrackPrediction.h.

Constructor & Destructor Documentation

◆ GhostTrackPrediction() [1/6]

reco::GhostTrackPrediction::GhostTrackPrediction ( )

inline

Definition at line 29 of file GhostTrackPrediction.h.

29 {}

◆ GhostTrackPrediction() [2/6]

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

inline

Definition at line 30 of file GhostTrackPrediction.h.

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

◆ GhostTrackPrediction() [3/6]

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

Definition at line 94 of file GhostTrackPrediction.cc.

                                                               {
   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);
 }

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

◆ GhostTrackPrediction() [4/6]

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

inline

Definition at line 37 of file GhostTrackPrediction.h.

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

References direction(), and init().

◆ GhostTrackPrediction() [5/6]

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

Definition at line 120 of file GhostTrackPrediction.cc.

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

◆ GhostTrackPrediction() [6/6]

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

◆ cartesianError()

Matrix6S GhostTrackPrediction::cartesianError ( double lambda = 0. ) const

Definition at line 146 of file GhostTrackPrediction.cc.

                                                                  {
   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());
 }

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

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

◆ cotTheta()

double reco::GhostTrackPrediction::cotTheta ( ) const

inline

Definition at line 50 of file GhostTrackPrediction.h.

50 { return prediction_[2]; }

References prediction_.

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

◆ covariance()

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

inline

Definition at line 60 of file GhostTrackPrediction.h.

60 { return covariance_; }

References covariance_.

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

◆ curvilinearError()

CurvilinearTrajectoryError GhostTrackPrediction::curvilinearError ( ) const

Definition at line 168 of file GhostTrackPrediction.cc.

                                                                         {
   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);
 }

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

Referenced by fts(), and track().

◆ direction()

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

inline

Definition at line 63 of file GhostTrackPrediction.h.

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

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

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

◆ eta()

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

inline

Definition at line 57 of file GhostTrackPrediction.h.

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

References cotTheta(), dqm-mbProfile::log, and rho().

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

◆ fts()

FreeTrajectoryState GhostTrackPrediction::fts ( const MagneticField * fieldProvider ) const

Definition at line 185 of file GhostTrackPrediction.cc.

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

References curvilinearError(), and globalTrajectory().

Referenced by transientGhostTrack().

◆ globalTrajectory()

GlobalTrajectoryParameters GhostTrackPrediction::globalTrajectory ( const MagneticField * fieldProvider ) const

Definition at line 164 of file GhostTrackPrediction.cc.

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

References direction(), and origin().

Referenced by fts().

◆ init()

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

private

Definition at line 58 of file GhostTrackPrediction.cc.

                                                                    {
   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);
 }

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

Referenced by GhostTrackPrediction().

◆ ip()

double reco::GhostTrackPrediction::ip ( ) const

inline

Definition at line 49 of file GhostTrackPrediction.h.

49 { return prediction_[1]; }

References prediction_.

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

◆ lambda()

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

inline

Definition at line 65 of file GhostTrackPrediction.h.

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

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

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

◆ origin()

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

inline

Definition at line 62 of file GhostTrackPrediction.h.

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

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().

◆ phi()

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

inline

Definition at line 51 of file GhostTrackPrediction.h.

51 { return prediction_[3]; }

References prediction_.

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

◆ position()

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

inline

Definition at line 67 of file GhostTrackPrediction.h.

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

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

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

◆ positionError()

GlobalError GhostTrackPrediction::positionError ( double lambda = 0. ) const

Definition at line 131 of file GhostTrackPrediction.cc.

                                                                    {
   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());
 }

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

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

◆ prediction()

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

inline

Definition at line 59 of file GhostTrackPrediction.h.

59 { return prediction_; }

References prediction_.

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

◆ rho()

double reco::GhostTrackPrediction::rho ( ) const

inline

Definition at line 54 of file GhostTrackPrediction.h.

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

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

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

◆ rho2()

double reco::GhostTrackPrediction::rho2 ( ) const

inline

Definition at line 53 of file GhostTrackPrediction.h.

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

References cotTheta().

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

◆ sz()

double reco::GhostTrackPrediction::sz ( ) const

inline

Definition at line 55 of file GhostTrackPrediction.h.

55 { return z() / rho(); }

References rho(), and z().

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

◆ theta()

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

inline

Definition at line 56 of file GhostTrackPrediction.h.

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

References cotTheta(), and M_PI_2.

Referenced by Tau.Tau::zImpact().

◆ track()

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

Definition at line 189 of file GhostTrackPrediction.cc.

                                                                 {
   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());
 }

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

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

◆ z()

double reco::GhostTrackPrediction::z ( ) const

inline

Definition at line 48 of file GhostTrackPrediction.h.

48 { return prediction_[0]; }

References prediction_.

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

Member Data Documentation

◆ covariance_

Error reco::GhostTrackPrediction::covariance_

private

Definition at line 86 of file GhostTrackPrediction.h.

Referenced by covariance(), and init().

◆ prediction_

Vector reco::GhostTrackPrediction::prediction_

private

Definition at line 85 of file GhostTrackPrediction.h.

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

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ CartesianError

◆ Error

◆ Vector

Constructor & Destructor Documentation

◆ GhostTrackPrediction() [1/6]

◆ GhostTrackPrediction() [2/6]

◆ GhostTrackPrediction() [3/6]

◆ GhostTrackPrediction() [4/6]

◆ GhostTrackPrediction() [5/6]

◆ GhostTrackPrediction() [6/6]

Member Function Documentation

◆ cartesianError()

◆ cotTheta()

◆ covariance()

◆ curvilinearError()

◆ direction()

◆ eta()

◆ fts()

◆ globalTrajectory()

◆ init()

◆ ip()

◆ lambda()

◆ origin()

◆ phi()

◆ position()

◆ positionError()

◆ prediction()

◆ rho()

◆ rho2()

◆ sz()

◆ theta()

◆ track()

◆ z()

Member Data Documentation

◆ covariance_

◆ prediction_