#include "fastjet/ClusterSequence.hh"
#include "UWEvent.hh"
#include <iostream>

using namespace std;
using namespace fastjet;

int main() {
  double R = 0.6;
  JetDefinition jet_def(cambridge_algorithm, R);

  vector<PseudoJet> event;
  int iev = 0;

  while (readUWEvent(cin, event) && ++iev <= 5) {
    cout << "Event " << iev << ", size = " << event.size() << endl;

    // run the jet clustering
    ClusterSequence cs(event, jet_def);
    double ptmin = 500.0;
    vector<PseudoJet> jets = sorted_by_pt(cs.inclusive_jets(ptmin));

    // loop over the jets
    for (unsigned i = 0; i < jets.size(); i++) {
      const PseudoJet & jet = jets[i];
      cout << "  jet " << i << ":" 
	   << "pt = " << jet.perp() << ", "
	   << "phi = " << jet.phi() << ", "
	   << "m = "  << jet.m() << endl;
    }
  }
}