GsfPropagatorAdapter.cc

Go to the documentation of this file.

#include "TrackingTools/GsfTools/interface/GsfPropagatorAdapter.h"

#include "TrackingTools/GsfTools/src/MultiStatePropagation.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include <atomic>

GsfPropagatorAdapter::GsfPropagatorAdapter (const Propagator& aPropagator) :
  Propagator(aPropagator.propagationDirection()),
  thePropagator(aPropagator.clone()) {}

std::pair<TrajectoryStateOnSurface,double>
GsfPropagatorAdapter::propagateWithPath (const TrajectoryStateOnSurface& tsos,
                     const Plane& plane) const {
  MultiStatePropagation<Plane> multiPropagator(*thePropagator);
  return multiPropagator.propagateWithPath(tsos,plane);
}

std::pair<TrajectoryStateOnSurface,double>
GsfPropagatorAdapter::propagateWithPath (const TrajectoryStateOnSurface& tsos,
                     const Cylinder& cylinder) const {
  MultiStatePropagation<Cylinder> multiPropagator(*thePropagator);
  return multiPropagator.propagateWithPath(tsos,cylinder);
}

std::pair<TrajectoryStateOnSurface,double>
GsfPropagatorAdapter::propagateWithPath (const FreeTrajectoryState& fts,
                     const Plane& plane) const {
  static std::atomic<int> nWarn{0};
  if ( nWarn++<5 )
    edm::LogInfo("GsfPropagatorAdapter") << "GsfPropagator used from FTS = single state mode!";
  return thePropagator->propagateWithPath(fts,plane);
}

std::pair<TrajectoryStateOnSurface,double>
GsfPropagatorAdapter::propagateWithPath (const FreeTrajectoryState& fts,
                     const Cylinder& cylinder) const {
  static std::atomic<int> nWarn{0};
  if ( nWarn++<5 )
    edm::LogInfo("GsfPropagatorAdapter") << "GsfPropagator used from FTS = single state mode!";
  return thePropagator->propagateWithPath(fts,cylinder);
}

void GsfPropagatorAdapter::setPropagationDirection (PropagationDirection dir) {
  thePropagator->setPropagationDirection(dir);
  Propagator::setPropagationDirection(dir);
}