BeamHaloPropagator.h

Go to the documentation of this file.

#ifndef GeomPropagators_BeamHaloPropagator_H
#define GeomPropagators_BeamHaloPropagator_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 BeamHaloPropagator Interface */

class BeamHaloPropagator final : public Propagator {
public:
  /* Constructor */
  BeamHaloPropagator(const Propagator* aEndCapTkProp,
                     const Propagator* aCrossTkProp,
                     const MagneticField* field,
                     PropagationDirection dir = alongMomentum);

  BeamHaloPropagator(const Propagator& aEndCapTkProp,
                     const Propagator& aCrossTkProp,
                     const MagneticField* field,
                     PropagationDirection dir = alongMomentum);

  BeamHaloPropagator(const BeamHaloPropagator&);

  ~BeamHaloPropagator() override;

  BeamHaloPropagator* clone() const override {
    return new BeamHaloPropagator(
        getEndCapTkPropagator(), getCrossTkPropagator(), magneticField(), propagationDirection());
  }

  void setPropagationDirection(PropagationDirection dir) override {
    Propagator::setPropagationDirection(dir);
    theEndCapTkProp->setPropagationDirection(dir);
    theCrossTkProp->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;

  bool crossingTk(const FreeTrajectoryState& fts, const Plane& plane) const;

  const Propagator* getEndCapTkPropagator() const;
  const Propagator* getCrossTkPropagator() const;
  const MagneticField* magneticField() const override { return theField; }

private:
  void directionCheck(PropagationDirection dir);

  Propagator* theEndCapTkProp;
  Propagator* theCrossTkProp;
  const MagneticField* theField;

protected:
};

#endif