TransientTrackFromFTS.cc

Go to the documentation of this file.

#include "TrackingTools/TransientTrack/interface/TransientTrackFromFTS.h"
#include "DataFormats/Math/interface/Error.h"
#include "DataFormats/Math/interface/Vector3D.h"
#include "TrackingTools/PatternTools/interface/TransverseImpactPointExtrapolator.h"
#include "Geometry/Records/interface/GlobalTrackingGeometryRecord.h"
#include "TrackingTools/PatternTools/interface/TSCBLBuilderNoMaterial.h"
// #include <iostream>
#include "FWCore/Utilities/interface/Likely.h"

using namespace reco;

TransientTrackFromFTS::TransientTrackFromFTS()
    : hasTime(false),
      timeExt_(0.),
      dtErrorExt_(0.),
      theField(nullptr),
      initialTSOSAvailable(false),
      initialTSCPAvailable(false),
      trackAvailable(false),
      blStateAvailable(false) {}

TransientTrackFromFTS::TransientTrackFromFTS(const FreeTrajectoryState& fts)
    : initialFTS(fts),
      hasTime(false),
      timeExt_(0.),
      dtErrorExt_(0.),
      theField(&(initialFTS.parameters().magneticField())),
      initialTSOSAvailable(false),
      initialTSCPAvailable(false),
      trackAvailable(false),
      blStateAvailable(false) {}

TransientTrackFromFTS::TransientTrackFromFTS(const FreeTrajectoryState& fts, const double time, const double dtime)
    : initialFTS(fts),
      hasTime(true),
      timeExt_(time),
      dtErrorExt_(dtime),
      theField(&(initialFTS.parameters().magneticField())),
      initialTSOSAvailable(false),
      initialTSCPAvailable(false),
      trackAvailable(false),
      blStateAvailable(false) {}

TransientTrackFromFTS::TransientTrackFromFTS(const FreeTrajectoryState& fts,
                                             const edm::ESHandle<GlobalTrackingGeometry>& tg)
    : initialFTS(fts),
      hasTime(false),
      timeExt_(0.),
      dtErrorExt_(0.),
      theField(&(initialFTS.parameters().magneticField())),
      initialTSOSAvailable(false),
      initialTSCPAvailable(false),
      trackAvailable(false),
      blStateAvailable(false),
      theTrackingGeometry(tg) {}

TransientTrackFromFTS::TransientTrackFromFTS(const FreeTrajectoryState& fts,
                                             const double time,
                                             const double dtime,
                                             const edm::ESHandle<GlobalTrackingGeometry>& tg)
    : initialFTS(fts),
      hasTime(true),
      timeExt_(time),
      dtErrorExt_(dtime),
      theField(&(initialFTS.parameters().magneticField())),
      initialTSOSAvailable(false),
      initialTSCPAvailable(false),
      trackAvailable(false),
      blStateAvailable(false),
      theTrackingGeometry(tg) {}

TransientTrackFromFTS::TransientTrackFromFTS(const TransientTrackFromFTS& tt)
    : initialFTS(tt.initialFreeState()),
      hasTime(tt.hasTime),
      timeExt_(tt.timeExt_),
      dtErrorExt_(tt.dtErrorExt_),
      theField(tt.field()),
      initialTSOSAvailable(false),
      initialTSCPAvailable(false),
      trackAvailable(false) {
  if (tt.initialTSOSAvailable) {
    initialTSOS = tt.impactPointState();
    initialTSOSAvailable = true;
  }
  if (tt.initialTSCPAvailable) {
    initialTSCP = tt.impactPointTSCP();
    initialTSCPAvailable = true;
  }
}

void TransientTrackFromFTS::setES(const edm::EventSetup& setup) {
  setup.get<GlobalTrackingGeometryRecord>().get(theTrackingGeometry);
}

void TransientTrackFromFTS::setTrackingGeometry(const edm::ESHandle<GlobalTrackingGeometry>& tg) {
  theTrackingGeometry = tg;
}

void TransientTrackFromFTS::setBeamSpot(const BeamSpot& beamSpot) {
  theBeamSpot = beamSpot;
  blStateAvailable = false;
}

TrajectoryStateOnSurface TransientTrackFromFTS::impactPointState() const {
  if
    UNLIKELY(!initialTSOSAvailable) calculateTSOSAtVertex();
  return initialTSOS;
}

TrajectoryStateClosestToPoint TransientTrackFromFTS::impactPointTSCP() const {
  if
    UNLIKELY(!initialTSCPAvailable) {
      initialTSCP = builder(initialFTS, initialFTS.position());
      initialTSCPAvailable = true;
    }
  return initialTSCP;
}

TrajectoryStateOnSurface TransientTrackFromFTS::outermostMeasurementState() const {
  throw cms::Exception("LogicError") << "TransientTrack built from a FreeTrajectoryState (TransientTrackFromFTS) can "
                                        "not have an outermostMeasurementState";
}

TrajectoryStateOnSurface TransientTrackFromFTS::innermostMeasurementState() const {
  throw cms::Exception("LogicError") << "TransientTrack built from a FreeTrajectoryState (TransientTrackFromFTS) can "
                                        "not have an innermostMeasurementState";
}

void TransientTrackFromFTS::calculateTSOSAtVertex() const {
  TransverseImpactPointExtrapolator tipe(theField);
  initialTSOS = tipe.extrapolate(initialFTS, initialFTS.position());
  initialTSOSAvailable = true;
}

TrajectoryStateOnSurface TransientTrackFromFTS::stateOnSurface(const GlobalPoint& point) const {
  TransverseImpactPointExtrapolator tipe(theField);
  return tipe.extrapolate(initialFTS, point);
}

const Track& TransientTrackFromFTS::track() const {
  if
    UNLIKELY(!trackAvailable) {
      GlobalPoint v = initialFTS.position();
      math::XYZPoint pos(v.x(), v.y(), v.z());
      GlobalVector p = initialFTS.momentum();
      math::XYZVector mom(p.x(), p.y(), p.z());

      theTrack = Track(0., 0., pos, mom, initialFTS.charge(), initialFTS.curvilinearError());
      trackAvailable = true;
    }
  return theTrack;
}

TrajectoryStateClosestToBeamLine TransientTrackFromFTS::stateAtBeamLine() const {
  if
    UNLIKELY(!blStateAvailable) {
      TSCBLBuilderNoMaterial blsBuilder;
      trajectoryStateClosestToBeamLine = blsBuilder(initialFTS, theBeamSpot);
      blStateAvailable = true;
    }
  return trajectoryStateClosestToBeamLine;
}