GsfTransientTrack.h

Go to the documentation of this file.

#ifndef TrackReco_GsfTransientTrack_h
#define TrackReco_GsfTransientTrack_h
 
#include "TrackingTools/TransientTrack/interface/BasicTransientTrack.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
#include "DataFormats/GsfTrackReco/interface/GsfTrack.h"
#include "DataFormats/GsfTrackReco/interface/GsfTrackFwd.h"
#include "TrackingTools/PatternTools/interface/TSCPBuilderNoMaterial.h"
#include "TrackingTools/PatternTools/interface/TransverseImpactPointExtrapolator.h"
 
namespace reco {
 
  class GsfTransientTrack : public GsfTrack, public BasicTransientTrack {
  public:
    // constructor from persistent track
    GsfTransientTrack();
    GsfTransientTrack(const GsfTrack& tk, const MagneticField* field);
    GsfTransientTrack(const GsfTrack& tk, const double time, const double dtime, const MagneticField* field);
 
    GsfTransientTrack(const GsfTrackRef& tk, const MagneticField* field);
    GsfTransientTrack(const GsfTrackRef& tk, const double time, const double dtime, const MagneticField* field);
 
    GsfTransientTrack(const GsfTrackRef& tk,
                      const MagneticField* field,
                      const edm::ESHandle<GlobalTrackingGeometry>& trackingGeometry);
    GsfTransientTrack(const GsfTrackRef& tk,
                      const double time,
                      const double dtime,
                      const MagneticField* field,
                      const edm::ESHandle<GlobalTrackingGeometry>& trackingGeometry);
 
    GsfTransientTrack(const GsfTrack& tk,
                      const MagneticField* field,
                      const edm::ESHandle<GlobalTrackingGeometry>& trackingGeometry);
    GsfTransientTrack(const GsfTrack& tk,
                      const double time,
                      const double dtime,
                      const MagneticField* field,
                      const edm::ESHandle<GlobalTrackingGeometry>& trackingGeometry);
 
    GsfTransientTrack(const GsfTransientTrack& tt);
 
    GsfTransientTrack& operator=(const GsfTransientTrack& tt);
 
    void setES(const edm::EventSetup&) override;
 
    void setTrackingGeometry(const edm::ESHandle<GlobalTrackingGeometry>&) override;
 
    void setBeamSpot(const reco::BeamSpot& beamSpot) override;
 
    FreeTrajectoryState initialFreeState() const override { return initialFTS; }
 
    TrajectoryStateOnSurface outermostMeasurementState() const override;
 
    TrajectoryStateOnSurface innermostMeasurementState() const override;
 
    TrajectoryStateClosestToPoint trajectoryStateClosestToPoint(const GlobalPoint& point) const override;
 
    TrajectoryStateOnSurface stateOnSurface(const GlobalPoint& point) const override;
 
    TrajectoryStateClosestToPoint impactPointTSCP() const override;
 
    TrajectoryStateOnSurface impactPointState() const override;
 
    bool impactPointStateAvailable() const override { return initialTSOSAvailable; }
 
    GsfTrackRef persistentTrackRef() const { return tkr_; }
 
    TrackBaseRef trackBaseRef() const override { return TrackBaseRef(tkr_); }
 
    TrackCharge charge() const override { return GsfTrack::charge(); }
 
    const MagneticField* field() const override { return theField; }
 
    const Track& track() const override { return *this; }
 
    TrajectoryStateClosestToBeamLine stateAtBeamLine() const override;
 
    double timeExt() const override { return (hasTime ? timeExt_ : std::numeric_limits<double>::quiet_NaN()); }
 
    double dtErrorExt() const override { return (hasTime ? dtErrorExt_ : std::numeric_limits<double>::quiet_NaN()); }
 
  private:
    void init();
    void calculateTSOSAtVertex() const;
 
    GsfTrackRef tkr_;
    bool hasTime;
    double timeExt_, dtErrorExt_;
    const MagneticField* theField;
 
    FreeTrajectoryState initialFTS;
    mutable bool initialTSOSAvailable, initialTSCPAvailable, blStateAvailable;
    mutable TrajectoryStateOnSurface initialTSOS;
    mutable TrajectoryStateClosestToPoint initialTSCP;
    TSCPBuilderNoMaterial builder;
    edm::ESHandle<GlobalTrackingGeometry> theTrackingGeometry;
    reco::BeamSpot theBeamSpot;
    mutable TrajectoryStateClosestToBeamLine trajectoryStateClosestToBeamLine;
 
    TransverseImpactPointExtrapolator theTIPExtrapolator;
  };
 
}  // namespace reco
 
#endif