TwoTrackMinimumDistance.h

Go to the documentation of this file.

#ifndef _TwoTrackMinimumDistance_H_
#define _TwoTrackMinimumDistance_H_

#include "TrackingTools/PatternTools/interface/ClosestApproachOnHelices.h"
#include "TrackingTools/PatternTools/interface/ClosestApproachInRPhi.h"
#include "TrackingTools/PatternTools/interface/TwoTrackMinimumDistanceHelixHelix.h"
#include "TrackingTools/PatternTools/interface/TwoTrackMinimumDistanceLineLine.h"
#include "TrackingTools/PatternTools/interface/TwoTrackMinimumDistanceHelixLine.h"
#include "FWCore/Utilities/interface/Visibility.h"

class TwoTrackMinimumDistance final : public ClosestApproachOnHelices {
public:
  enum Mode { FastMode = 0, SlowMode = 1 };

  TwoTrackMinimumDistance(const Mode m = FastMode) {
    theModus = m;
    status_ = false;
  };
  ~TwoTrackMinimumDistance() override {}

  bool calculate(const TrajectoryStateOnSurface& sta, const TrajectoryStateOnSurface& stb) override;

  bool calculate(const FreeTrajectoryState& sta, const FreeTrajectoryState& stb) override;

  virtual bool calculate(const GlobalTrajectoryParameters& sta, const GlobalTrajectoryParameters& stb);

  bool status() const override { return status_; }

  std::pair<GlobalPoint, GlobalPoint> points() const override;

  GlobalPoint crossingPoint() const override;

  float distance() const override;

  TwoTrackMinimumDistance* clone() const override { return new TwoTrackMinimumDistance(*this); }

  double firstAngle() const;
  double secondAngle() const;
  std::pair<double, double> pathLength() const;

private:
  enum Charge { hh, hl, ll };
  Mode theModus;
  Charge theCharge;
  ClosestApproachInRPhi theIniAlgo;
  TwoTrackMinimumDistanceHelixHelix theTTMDhh;
  TwoTrackMinimumDistanceLineLine theTTMDll;
  TwoTrackMinimumDistanceHelixLine theTTMDhl;
  bool status_;
  std::pair<GlobalPoint, GlobalPoint> points_;

  bool pointsLineLine(const GlobalTrajectoryParameters& sta, const GlobalTrajectoryParameters& stb) dso_internal;
  bool pointsHelixLine(const GlobalTrajectoryParameters& sta, const GlobalTrajectoryParameters& stb) dso_internal;
  bool pointsHelixHelix(const GlobalTrajectoryParameters& sta, const GlobalTrajectoryParameters& stb) dso_internal;
};

#endif