56 vector<fastjet::PseudoJet> input_particles;
58 double px, py , pz, E;
59 while (cin >> px >> py >> pz >> E) {
81 vector<fastjet::PseudoJet> inclusive_jets = sorted_by_pt(clust_seq.
inclusive_jets(ptmin));
93 printf(
"%5s %15s %15s %15s\n",
"jet #",
"rapidity",
"phi",
"pt");
96 for (
unsigned int i = 0; i < inclusive_jets.size(); i++) {
97 printf(
"%5u %15.8f %15.8f %15.8f\n",
98 i, inclusive_jets[i].rap(), inclusive_jets[i].phi(),
99 inclusive_jets[i].perp());