testtrackfit.C

Go to the documentation of this file.

// testtrackfit
// Test program to test fitter classes for tracks 
// B. List, J. List, 17.11.04


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

#include "jbltools/kinfit/ChargedParticleTrack.h"
#include "jbltools/kinfit/TTVertexConstraint.h"
#include "jbltools/kinfit/ThreeVector.h"
#include "jbltools/kinfit/FourVector.h"
#include "jbltools/kinfit/WWFitter.h"


int main(int argc, char** argv) {

  TrackFitObject::setBfield (11.4);
   
  ChargedParticleTrack t1 ("pi1", 0.01, 1.57, 1.1, 1.0, 5.0, 0.135, 1, 0.0, 100.);
  
  cout << t1 << endl;
  
  for (double s = 0; s < 50; s += 1.0) {
    cout << "s=" << s << ", track=" << t1.getTrajectoryPoint (s) << endl;
  
  }

  ChargedParticleTrack t2 ("pi2", -0.015, 1.7, 1.3, 0.5, 4.0, 0.135, 1, 0.0, 100.);
  
  TTVertexConstraint vc_x (t1, 0, t2, 0, 0);
  TTVertexConstraint vc_y (t1, 0, t2, 0, 1);
  TTVertexConstraint vc_z (t1, 0, t2, 0, 2);
  
  cout << "t1 = " << t1 << "\nt2 = " << t2 << endl;
  cout << "p1 = " << t1.getMomentum(0) 
       << "\np2 = " << t2.getMomentum(0) << endl;
  
  WWFitter fitter;
  fitter.addFitObject (t1);
  fitter.addFitObject (t2);
  fitter.addConstraint (vc_x);
  fitter.addConstraint (vc_y);
  fitter.addConstraint (vc_z);
    
  double prob = fitter.fit();
  cout << "fit probability = " << prob << endl;
  cout << "error code: " << fitter.getError() << endl;
  cout << "t1 = " << t1 << "\nt2 = " << t2 << endl;
  cout << "t1.startvertex = " << t1.getVertex(0)
       << "\nt2.startvertex = " << t2.getVertex(0) << endl;
  cout << "p1 = " << t1.getMomentum(0) 
       << "\np2 = " << t2.getMomentum(0) << endl;
  
  return 0;
  
   
}

---