SmartPropagator.h

Go to the documentation of this file.

#ifndef GeomPropagators_SmartPropagator_H
#define GeomPropagators_SmartPropagator_H

/* Collaborating Class Declarations */
#include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "TrackingTools/GeomPropagators/interface/Propagator.h"
#include "MagneticField/Engine/interface/MagneticField.h"
#include "DataFormats/GeometrySurface/interface/ReferenceCounted.h"
#include "DataFormats/TrajectorySeed/interface/PropagationDirection.h"

class Cylinder;
class Plane;

/* Class SmartPropagator Interface */

class SmartPropagator final : public Propagator {
public:
  /* Constructor */
  SmartPropagator(const Propagator* aTkProp,
                  const Propagator* aGenProp,
                  const MagneticField* field,
                  PropagationDirection dir = alongMomentum,
                  float epsilon = 5);

  SmartPropagator(const Propagator& aTkProp,
                  const Propagator& aGenProp,
                  const MagneticField* field,
                  PropagationDirection dir = alongMomentum,
                  float epsilon = 5);

  SmartPropagator(const SmartPropagator&);

  ~SmartPropagator() override;

  SmartPropagator* clone() const override {
    return new SmartPropagator(getTkPropagator(), getGenPropagator(), magneticField());
  }

  void setPropagationDirection(PropagationDirection dir) override {
    Propagator::setPropagationDirection(dir);
    theTkProp->setPropagationDirection(dir);
    theGenProp->setPropagationDirection(dir);
  }

  using Propagator::propagate;
  using Propagator::propagateWithPath;

private:
  std::pair<TrajectoryStateOnSurface, double> propagateWithPath(const FreeTrajectoryState& fts,
                                                                const Plane& plane) const override;

  std::pair<TrajectoryStateOnSurface, double> propagateWithPath(const FreeTrajectoryState& fts,
                                                                const Cylinder& cylinder) const override;

  std::pair<TrajectoryStateOnSurface, double> propagateWithPath(const TrajectoryStateOnSurface& tsos,
                                                                const Plane& sur) const override;

  std::pair<TrajectoryStateOnSurface, double> propagateWithPath(const TrajectoryStateOnSurface& tsos,
                                                                const Cylinder& sur) const override;

public:
  bool insideTkVol(const FreeTrajectoryState& fts) const;
  bool insideTkVol(const Surface& surface) const;
  bool insideTkVol(const Cylinder& cylin) const;
  bool insideTkVol(const Plane& plane) const;

  const Propagator* getTkPropagator() const;
  const Propagator* getGenPropagator() const;
  const MagneticField* magneticField() const override { return theField; }

private:
  void initTkVolume(float epsilon);

  Propagator* theTkProp;
  Propagator* theGenProp;
  const MagneticField* theField;
  ReferenceCountingPointer<Cylinder> theTkVolume;

protected:
};

#endif  // SMARTPROPAGATOR_H