TwoTrackMinimumDistanceHelixHelix.h

Go to the documentation of this file.

#ifndef _Tracker_TwoTrackMinimumDistanceHelixHelix_H_
#define _Tracker_TwoTrackMinimumDistanceHelixHelix_H_

#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
// #include <string>
// #include <sstream>
#include <utility>

class GlobalTrajectoryParameters;

class TwoTrackMinimumDistanceHelixHelix {

public:
  TwoTrackMinimumDistanceHelixHelix();
  ~TwoTrackMinimumDistanceHelixHelix();

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

  std::pair <GlobalPoint, GlobalPoint> points() const {
    if (!pointsUpdated) finalPoints();
    return std::pair<GlobalPoint, GlobalPoint> (pointG, pointH);
  }

  std::pair <double, double> pathLength() const {
    if (!pointsUpdated) finalPoints();
    return std::pair <double, double> ( pathG, pathH);
  }



  double firstAngle() const {return thepG;}
  double secondAngle() const {return thepH;}

private:
  bool updateCoeffs( const GlobalPoint & , const GlobalPoint & );
  bool oneIteration ( double &, double & ) const;

// bool parallelTracks () const;
  void finalPoints() const;

private:
  GlobalTrajectoryParameters const *theH, *theG;
  // the 'GH-track data' (constants)
  double thea, theb, thec1, thec2, thed1, thed2, thee1, thee2, theg, theh;
  // double thelambdaG, thelambdaH;
  double thetanlambdaG, thetanlambdaH;
  double thesinpG0, thecospG0;
  double thesinpH0, thecospH0;
  double thepG0, thepH0;

  // the variable stuff
  // = the point we are currently looking at.
  mutable double thepG, thepH;
  mutable double thesinpG, thesinpH;
  mutable double thecospG, thecospH;
  mutable GlobalPoint pointG, pointH;
  mutable double pathG, pathH;
  mutable bool pointsUpdated;

  double themaxjump, thesingjacI;
  int themaxiter;

};
#endif