/data/refman/pasoursint/CMSSW_6_1_2_SLHC2_patch1/src/TrackingTools/PatternTools/interface/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(0), theL(0), 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:
  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() const;
  bool oneIteration(double & thePhiH, double & fct, double & derivative ) const;
  mutable GlobalPoint helixPoint, linePoint;
  mutable double tL, linePath, helixPath;
  mutable bool pointsUpdated;

};
#endif