55 vector<fastjet::PseudoJet> input_particles;
57 double px, py , pz, E;
58 while (cin >> px >> py >> pz >> E) {
97 unsigned int n_repeat = 3;
98 double ghost_area = 0.01;
114 vector<fastjet::PseudoJet> inclusive_jets = sorted_by_pt(clust_seq.
inclusive_jets(ptmin));
125 cout <<
"Area: " << area_def.
description() << endl << endl;
128 printf(
"%5s %15s %15s %15s %15s %15s\n",
"jet #",
"rapidity",
"phi",
"pt",
"area",
"area error");
131 for (
unsigned int i = 0; i < inclusive_jets.size(); i++) {
132 printf(
"%5u %15.8f %15.8f %15.8f %15.8f %15.8f\n", i,
133 inclusive_jets[i].rap(), inclusive_jets[i].phi(), inclusive_jets[i].perp(),
134 inclusive_jets[i].area(), inclusive_jets[i].area_error());