43 #include "fastjet/PseudoJet.hh"
44 #include "fastjet/ClusterSequenceArea.hh"
45 #include "fastjet/ClusterSequencePassiveArea.hh"
48 #include "fastjet/config.h"
50 #ifdef ENABLE_PLUGIN_SISCONE
51 #include "fastjet/SISConePlugin.hh"
61 void print_jets (
const vector<fastjet::PseudoJet> &);
66 vector<fastjet::PseudoJet> input_particles;
69 double px, py , pz, E;
70 while (cin >> px >> py >> pz >> E) {
87 bool use_voronoi =
false;
89 double ghost_etamax = 6.0;
90 double ghost_area = 0.01;
91 int active_area_repeats = 1;
100 double effective_Rfact = 1.0;
118 cout <<
"Jet definition was: " << jet_def.
description() << endl;
119 cout <<
"Area definition was: " << area_def.
description() << endl;
120 cout <<
"Strategy adopted by FastJet was "<<
125 vector<fastjet::PseudoJet> inclusive_jets = clust_seq.
inclusive_jets(ptmin);
128 cout <<
"Printing inclusive jets with pt > "<< ptmin<<
" GeV\n";
129 cout <<
"---------------------------------------\n";
134 cout <<
"Number of unclustered particles: "
143 void print_jets (
const vector<fastjet::PseudoJet> & unsorted_jets) {
146 vector<fastjet::PseudoJet> jets =
sorted_by_pt(unsorted_jets);
148 printf(
" ijet rap phi Pt area +- err\n");
149 for (
unsigned int j = 0; j < jets.size(); j++) {
151 double area = jets[j].area();
152 double area_error = jets[j].area_error();
154 printf(
"%5u %9.5f %8.5f %10.3f %8.3f +- %6.3f\n",j,jets[j].rap(),
155 jets[j].phi(),jets[j].perp(), area, area_error);