trajectoryStateTransform Namespace Reference

Functions
FreeTrajectoryState	initialFreeState (const reco::Track &tk, const MagneticField *field, bool withErr=true)
FreeTrajectoryState	innerFreeState (const reco::Track &tk, const MagneticField *field, bool withErr=true)
TrajectoryStateOnSurface	innerStateOnSurface (const reco::Track &tk, const TrackingGeometry &geom, const MagneticField *field, bool withErr=true)
FreeTrajectoryState	outerFreeState (const reco::Track &tk, const MagneticField *field, bool withErr=true)
TrajectoryStateOnSurface	outerStateOnSurface (const reco::Track &tk, const TrackingGeometry &geom, const MagneticField *field, bool withErr=true)
PTrajectoryStateOnDet	persistentState (const TrajectoryStateOnSurface &ts, unsigned int detid)
TrajectoryStateOnSurface	transientState (const PTrajectoryStateOnDet &ts, const Surface *surface, const MagneticField *field)

Function Documentation

FreeTrajectoryState trajectoryStateTransform::initialFreeState	(	const reco::Track &	tk,
		const MagneticField *	field,
		bool	withErr = true
	)

Construct a FreeTrajectoryState from the reco::Track innermost or outermost state, does not require access to tracking geometry

Definition at line 69 of file TrajectoryStateTransform.cc.

References reco::TrackBase::charge(), reco::TrackBase::covariance(), reco::TrackBase::momentum(), and reco::TrackBase::vertex().

Referenced by reco::CandidatePtrTransientTrack::CandidatePtrTransientTrack(), reco::GsfTransientTrack::GsfTransientTrack(), and reco::TrackTransientTrack::TrackTransientTrack().

  {
    Basic3DVector<float> pos( tk.vertex());
    GlobalPoint gpos( pos);
    Basic3DVector<float> mom( tk.momentum());
    GlobalVector gmom( mom);
    GlobalTrajectoryParameters par( gpos, gmom, tk.charge(), field);
    if (!withErr) return FreeTrajectoryState(par);
    CurvilinearTrajectoryError err( tk.covariance());
    return FreeTrajectoryState( par, err);
  }

FreeTrajectoryState trajectoryStateTransform::innerFreeState	(	const reco::Track &	tk,
		const MagneticField *	field,
		bool	withErr = true
	)

Definition at line 82 of file TrajectoryStateTransform.cc.

References reco::TrackBase::charge(), reco::Track::extra(), reco::Track::innerMomentum(), and reco::Track::innerPosition().

Referenced by innerStateOnSurface().

  {
    Basic3DVector<float> pos( tk.innerPosition());
    GlobalPoint gpos( pos);
    Basic3DVector<float> mom( tk.innerMomentum());
    GlobalVector gmom( mom);
    GlobalTrajectoryParameters par( gpos, gmom, tk.charge(), field);
    if (!withErr) return FreeTrajectoryState(par);
    CurvilinearTrajectoryError err( tk.extra()->innerStateCovariance());
    return FreeTrajectoryState( par, err);
  }

TrajectoryStateOnSurface trajectoryStateTransform::innerStateOnSurface	(	const reco::Track &	tk,
		const TrackingGeometry &	geom,
		const MagneticField *	field,
		bool	withErr = true
	)

Construct a TrajectoryStateOnSurface from the reco::Track innermost or outermost state, requires access to tracking geometry

Definition at line 110 of file TrajectoryStateTransform.cc.

References reco::Track::extra(), TrackingGeometry::idToDet(), innerFreeState(), and GeomDet::surface().

Referenced by reco::CandidatePtrTransientTrack::innermostMeasurementState(), and reco::TrackTransientTrack::innermostMeasurementState().

  {
    const Surface& surface = geom.idToDet( DetId( tk.extra()->innerDetId()))->surface();
    return TrajectoryStateOnSurface( innerFreeState( tk, field, withErr), surface);
  }

FreeTrajectoryState trajectoryStateTransform::outerFreeState	(	const reco::Track &	tk,
		const MagneticField *	field,
		bool	withErr = true
	)

Definition at line 96 of file TrajectoryStateTransform.cc.

References reco::TrackBase::charge(), reco::Track::extra(), reco::Track::outerMomentum(), and reco::Track::outerPosition().

Referenced by outerStateOnSurface().

  {
    Basic3DVector<float> pos( tk.outerPosition());
    GlobalPoint gpos( pos);
    Basic3DVector<float> mom( tk.outerMomentum());
    GlobalVector gmom( mom);
    GlobalTrajectoryParameters par( gpos, gmom, tk.charge(), field);
    if (!withErr) return FreeTrajectoryState(par);
    CurvilinearTrajectoryError err( tk.extra()->outerStateCovariance());
    return FreeTrajectoryState( par, err);
  }

TrajectoryStateOnSurface trajectoryStateTransform::outerStateOnSurface	(	const reco::Track &	tk,
		const TrackingGeometry &	geom,
		const MagneticField *	field,
		bool	withErr = true
	)

Definition at line 118 of file TrajectoryStateTransform.cc.

References reco::Track::extra(), TrackingGeometry::idToDet(), outerFreeState(), and GeomDet::surface().

Referenced by reco::CandidatePtrTransientTrack::outermostMeasurementState(), and reco::TrackTransientTrack::outermostMeasurementState().

  {
    const Surface& surface = geom.idToDet( DetId( tk.extra()->outerDetId()))->surface();
    return TrajectoryStateOnSurface( outerFreeState( tk, field, withErr), surface);
  }

PTrajectoryStateOnDet trajectoryStateTransform::persistentState	(	const TrajectoryStateOnSurface &	ts,
		unsigned int	detid
	)

should check if corresponds to m

should check if corresponds to m

Definition at line 16 of file TrajectoryStateTransform.cc.

References TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::hasError(), i, j, relval_2017::k, TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), visualization-live-secondInstance_cfg::m, LocalTrajectoryError::matrix(), PV3DBase< T, PVType, FrameType >::perp(), EnergyCorrector::pt, and TrajectoryStateOnSurface::surfaceSide().

Referenced by MuonSeedFromRecHits::createSeed(), SeedFromProtoTrack::init(), CosmicSeedCreator::makeSeed(), and TrajectorySeedFromMuonProducer::produce().

  {
    int surfaceSide = static_cast<int>(ts.surfaceSide());
    auto pt = ts.globalMomentum().perp();
    
    if (ts.hasError()) {
      AlgebraicSymMatrix55 const & m = ts.localError().matrix();
      
      int dim = 5; 
      float localErrors[15];
      
      int k = 0;
      for (int i=0; i<dim; i++) {
    for (int j=0; j<=i; j++) {
      localErrors[k++] = m(i,j);
    }
      }
    return PTrajectoryStateOnDet(ts.localParameters(),pt,
                 localErrors, detid,
                 surfaceSide);
    }
    return PTrajectoryStateOnDet(ts.localParameters(),pt,
                 detid,
                 surfaceSide);
  }

TrajectoryStateOnSurface trajectoryStateTransform::transientState	(	const PTrajectoryStateOnDet &	ts,
		const Surface *	surface,
		const MagneticField *	field
	)

Definition at line 44 of file TrajectoryStateTransform.cc.

References PTrajectoryStateOnDet::error(), PTrajectoryStateOnDet::hasError(), i, j, relval_2017::k, visualization-live-secondInstance_cfg::m, PTrajectoryStateOnDet::parameters(), and PTrajectoryStateOnDet::surfaceSide().

Referenced by SeedMatcher::matchRecHitCombinations(), SeedTransformer::seedTransientState(), and StartingLayerFinder::startingLayers().

  {
    AlgebraicSymMatrix55 m;
    bool errInv=true;
    if (ts.hasError()) {
      errInv = false;
      int dim = 5;
      int k = 0;
      for (int i=0; i<dim; i++) {
    for (int j=0; j<=i; j++) {
      m(i,j) = ts.error(k++);       // NOTE: here we do a cast float => double.     
    }
      }
    }
    
    
    return TrajectoryStateOnSurface( ts.parameters(),
                     errInv ? LocalTrajectoryError(InvalidError()) : LocalTrajectoryError(m),
                     *surface, field,
                     static_cast<SurfaceSide>(ts.surfaceSide()));
    
}