helixtest.C

Go to the documentation of this file.

// helixtest
// Test program to test the JBLHelix class:
// B. List 4.1.05


#include <iostream>              // - cout
using std::cout;
using std::endl;

#include "jbltools/kinfit/JBLHelix.h"
#include "jbltools/kinfit/ThreeVector.h"



int main(int argc, char** argv) {
  JBLHelix h1 ( 0.01, 0.23, 1.2, 1.5, 5.0);
  JBLHelix h2 (-0.02, 0.20, 1.4, 1.2, 6.0);
  
  for (double s = -314.; s < 314.; s+=50) {
    ThreeVector p = h1.getTrajectoryPoint (s);
    double sc = h1.getClosestS (p);
    cout << "h1: s=" << s << ", p=" << p << ", s_c=" << sc << ", diff=" << s-sc << endl;
  }
  for (double s = -157.; s < 157.; s+=10) {
    ThreeVector p = h2.getTrajectoryPoint (s);
    double sc = h2.getClosestS (p);
    cout << "h2: s=" << s << ", p=" << p << ", s_c=" << sc << ", diff=" << s-sc << endl;
  }
  
  double s01, s11, s02, s12;
  h1.getClosestApproach (h2, s01, s11, s02, s12);

  return 0;
}

---