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);
vector< PseudoJet > sorted_by_pt(const vector< PseudoJet > &jets)
return a vector of jets sorted into decreasing kt2
General class for user to obtain ClusterSequence with additional area information.
std::vector< PseudoJet > inclusive_jets(const double &ptmin=0.0) const
return a vector of all jets (in the sense of the inclusive algorithm) with pt >= ptmin.
the longitudinally invariant kt algorithm
class that holds a generic area definition
std::vector< PseudoJet > unclustered_particles() const
return the set of particles that have not been clustered.
int main()
an example program showing how to use fastjet
Specification for the computation of the Voronoi jet area.
std::string description() const
return a textual description of the current jet definition
std::string description() const
return a description of the current area definition
void print_jets(const vector< fastjet::PseudoJet > &)
a function that pretty prints a list of jets
automatic selection of the best (based on N)
Parameters to configure the computation of jet areas using ghosts.
Strategy
the various options for the algorithmic strategy to adopt in clustering events with kt and cambridge ...
Class to contain pseudojets, including minimal information of use to jet-clustering routines...
class that is intended to hold a full definition of the jet clusterer
std::string strategy_string() const
return the name of the strategy used to cluster the event