db/db8/TwoTrackMinimumDistanceHelixLine_8h_source.html

#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