#include "fastjet/PseudoJet.hh"
#include "fastjet/ClusterSequence.hh"
#include <iostream>
#include <sstream>
#include <vector>
Include dependency graph for run_jet_finder.cc:
Go to the source code of this file.
Functions | |
void | print_jets (const fastjet::ClusterSequence &, const vector< fastjet::PseudoJet > &) |
a function that pretty prints a list of jets | |
void | run_jet_finder (const vector< fastjet::PseudoJet > &input_particles, const fastjet::JetDefinition &jet_def) |
subroutine that runs fastjet (on input particles from cin) using an arbitrary jet definition supplied as an argument | |
void | read_input_particles (istream &input, vector< fastjet::PseudoJet > &input_particles) |
a function that reads particles from the supplied istream |
|
a function that pretty prints a list of jets
|
|
a function that reads particles from the supplied istream
Definition at line 73 of file run_jet_finder.cc. Referenced by main(). 00074 { 00075 00076 // read in input particles 00077 double px, py , pz, E; 00078 string line; 00079 while (getline(input, line)) { 00080 if (line.substr(0,1) == "#") {continue;} 00081 istringstream linestream(line); 00082 linestream >> px >> py >> pz >> E; 00083 // create a fastjet::PseudoJet with these components and put it onto 00084 // back of the input_particles vector 00085 input_particles.push_back(fastjet::PseudoJet(px,py,pz,E)); 00086 } 00087 00088 }
|
|
subroutine that runs fastjet (on input particles from cin) using an arbitrary jet definition supplied as an argument
Definition at line 47 of file run_jet_finder.cc. References fastjet::JetDefinition::description(), fastjet::ClusterSequence::inclusive_jets(), and print_jets(). Referenced by main(). 00048 { 00049 00050 // run the jet clustering with the above jet definition 00051 fastjet::ClusterSequence clust_seq(input_particles, jet_def); 00052 00053 // tell the user what was done 00054 cout << "Ran " << jet_def.description() << endl; 00055 00056 // extract the inclusive jets with pt > 5 GeV 00057 double ptmin = 5.0; 00058 vector<fastjet::PseudoJet> inclusive_jets = clust_seq.inclusive_jets(ptmin); 00059 00060 // print them out 00061 cout << "Printing inclusive jets with pt > "<< ptmin<<" GeV\n"; 00062 cout << "---------------------------------------\n"; 00063 print_jets(clust_seq, inclusive_jets); 00064 cout << endl; 00065 00066 // print out unclustered stuff 00067 cout << clust_seq.unclustered_particles().size() << " particles unclustered" << endl << endl; 00068 }
|