TwoTrackMinimumDistanceHelixLine.h

Go to the documentation of this file.

#ifndef _Tracker_TwoTrackMinimumDistanceHelixLine_H_
#define _Tracker_TwoTrackMinimumDistanceHelixLine_H_

#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/GeometryVector/interface/GlobalVector.h"
#include <utility>
class GlobalTrajectoryParameters;

class TwoTrackMinimumDistanceHelixLine {
public:
  TwoTrackMinimumDistanceHelixLine() : theH(nullptr), theL(nullptr), themaxiter(12), pointsUpdated(false) {}
  ~TwoTrackMinimumDistanceHelixLine() {}

  bool calculate(const GlobalTrajectoryParameters &,
                 const GlobalTrajectoryParameters &,
                 const float qual = .0001);  // retval=true? error occured.

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

  double firstAngle() const;
  double secondAngle() const;

private:
  const GlobalTrajectoryParameters *theH, *theL, *firstGTP, *secondGTP;
  GlobalVector posDiff;
  GlobalVector theLp;
  double X, Y, Z, px, py, pz, px2, py2, pz2, baseFct, baseDer;
  double theh, thePhiH0, thesinPhiH0, thecosPhiH0, thetanlambdaH;
  double thePhiH;
  double Hn, Ln;
  double aa, bb, cc, dd, ee, ff;

  int themaxiter;
  bool updateCoeffs();
  void finalPoints();
  bool oneIteration(double &thePhiH, double &fct, double &derivative) const;
  GlobalPoint helixPoint, linePoint;
  double tL, linePath, helixPath;
  bool pointsUpdated;
};
#endif