/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/TrackingTools/PatternTools/interface/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 : public ClosestApproachOnHelices {

public:

  enum Mode { FastMode=0, SlowMode=1 };

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

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

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

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

  virtual bool status() const {return status_;}

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

  virtual GlobalPoint crossingPoint() const;

  virtual float distance() const;


  virtual TwoTrackMinimumDistance * clone() const {
    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);
  bool pointsHelixLine(const GlobalTrajectoryParameters & sta,
        const GlobalTrajectoryParameters & stb);
  bool pointsHelixHelix(const GlobalTrajectoryParameters & sta,
        const GlobalTrajectoryParameters & stb);
};

#endif