42 #include "fastjet/ClusterSequenceArea.hh"
53 vector<fastjet::PseudoJet> input_particles;
55 double px, py , pz, E;
56 while (cin >> px >> py >> pz >> E) {
95 unsigned int n_repeat = 3;
96 double ghost_area = 0.01;
112 vector<fastjet::PseudoJet> inclusive_jets = sorted_by_pt(clust_seq.
inclusive_jets(ptmin));
123 cout <<
"Area: " << area_def.
description() << endl << endl;
126 printf(
"%5s %15s %15s %15s %15s %15s\n",
"jet #",
"rapidity",
"phi",
"pt",
"area",
"area error");
129 for (
unsigned int i = 0; i < inclusive_jets.size(); i++) {
130 printf(
"%5u %15.8f %15.8f %15.8f %15.8f %15.8f\n", i,
131 inclusive_jets[i].rap(), inclusive_jets[i].phi(), inclusive_jets[i].perp(),
132 inclusive_jets[i].area(), inclusive_jets[i].area_error());