source: huonglan/ANA2011/simpletest.C

#include <iostream>
#include <math.h>
#include <utility>
#include "TLorentzVector.h"
#include "TVector2.h"
#include "TVector3.h"
#include "Thrust_func_collection.h"

using namespace std;

simpletest()
{
  TLorentzVector v1(10,0,0,0);
  TLorentzVector v2(-10,0,0,0);
  TLorentzVector v3(0,1,0,0);

  std::vector<TLorentzVector> vec;
  vec.push_back(v1);
  vec.push_back(v2);
  vec.push_back(v3);

  Thrust_func_collection obj;

  std::vector<double> result = obj.Calc_Thrust(vec);

  std::cout << "result1 = " << result.at(0) << std::endl;
  std::cout << "result2 = " << result.at(1) << std::endl;
  
}