testwwfit.C

Go to the documentation of this file.

// testwwfit
// Test program to test fitter classes 
// J. Boehme, B. List 3.4.04


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

#include "jbltools/kinfit/JetFitObject.h"
#include "jbltools/kinfit/PxConstraint.h"
#include "jbltools/kinfit/PyConstraint.h"
#include "jbltools/kinfit/WWFitter.h"


int main(int argc, char** argv) {
 
   
  JetFitObject j1 (44., 1.571, 0.087, 5.0, 0.087, 0.1, 0.);
  JetFitObject j2 (46., 1.571, 3.220, 5.0, 0.087, 0.1, 0.);
  cout << "start four-vector of first  jet: " << j1 << endl;
  cout << "start four-vector of second jet: " << j2 << endl;
  
  PxConstraint pxc;
  pxc.addToFOList (j1);
  pxc.addToFOList (j2);
  
  PyConstraint pyc;
  pyc.addToFOList (j1);
  pyc.addToFOList (j2);
  
  WWFitter fitter;
  fitter.addFitObject (j1);
  fitter.addFitObject (j2);
  fitter.addConstraint (pxc);
  fitter.addConstraint (pyc);
  
  double prob = fitter.fit();
  cout << "fit probability = " << prob << endl;
  cout << "error code: " << fitter.getError() << endl;
  cout << "final four-vector of first  jet: " << j1 << endl;
  cout << "final four-vector of second jet: " << j2 << endl;
  
  return 0;
  
  
   
}

---