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"

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;
  mutable Charge theCharge;
  ClosestApproachInRPhi theIniAlgo;
  mutable TwoTrackMinimumDistanceHelixHelix theTTMDhh;
  mutable TwoTrackMinimumDistanceLineLine theTTMDll;
  mutable 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