#include "fastjet/config.h"
#include "run_jet_finder.hh"
#include <vector>
#include <iostream>
#include "fastjet/PseudoJet.hh"
Include dependency graph for many_algs_example.cc:
Go to the source code of this file.
Functions | |
int | main (int argc, char **argv) |
int main | ( | int | argc, | |
char ** | argv | |||
) |
Definition at line 62 of file many_algs_example.cc.
References fastjet::antikt_algorithm, fastjet::cambridge_algorithm, fastjet::kt_algorithm, read_input_particles(), and run_jet_finder().
00062 { 00063 00064 // read input particles 00065 vector<fastjet::PseudoJet> input_particles; 00066 read_input_particles(cin, input_particles); 00067 00068 // we will have four jet definitions, and the first two will be 00069 // plugins 00070 vector<fastjet::JetDefinition> jet_defs(6); 00071 vector<fastjet::JetDefinition::Plugin *> plugins(3); 00072 00073 // common parameters 00074 double jet_radius = 0.7; 00075 //double jet_radius = 1.0; 00076 double overlap_threshold = 0.5; 00077 00078 vector<fastjet::JetDefinition>::const_iterator 00079 jet_def_begin = jet_defs.begin(); 00080 00081 // set up a pxcone jet definition (if wanted -- requires f77, and you 00082 // should compile the pxcone plugin (not there by default)) 00083 #ifdef ENABLE_PLUGIN_PXCONE 00084 double min_jet_energy = 5.0; 00085 bool E_scheme_jets = false; 00086 plugins[0] = new fastjet::PxConePlugin (jet_radius, min_jet_energy, 00087 overlap_threshold, E_scheme_jets); 00088 jet_defs[0] = fastjet::JetDefinition(plugins[0]); 00089 #else 00090 plugins[0] = NULL; 00091 //jet_def_begin++; // skip first jet def 00092 #endif // ENABLE_PLUGIN_PXCONE 00093 00094 00095 00096 // set up a CDF midpoint jet definition 00097 #ifdef ENABLE_PLUGIN_CDFCONES 00098 double seed_threshold = 1.0; 00099 double cone_area_fraction = 1.0; 00100 int max_pair_size = 2; 00101 int max_iterations = 100; 00102 plugins[1] = new fastjet::CDFMidPointPlugin (seed_threshold, jet_radius, 00103 cone_area_fraction, max_pair_size, 00104 max_iterations, overlap_threshold); 00105 jet_defs[1] = fastjet::JetDefinition(plugins[1]); 00106 #else 00107 plugins[1] = NULL; 00108 #endif 00109 00110 // set up a siscone jet definition 00111 #ifdef ENABLE_PLUGIN_SISCONE 00112 int npass = 0; // do infinite number of passes 00113 double protojet_ptmin = 0.0; // use all protojets 00114 plugins[2] = new fastjet::SISConePlugin (jet_radius, overlap_threshold, 00115 npass, protojet_ptmin); 00116 jet_defs[2] = fastjet::JetDefinition(plugins[2]); 00117 #else 00118 plugins[2] = NULL; 00119 #endif 00120 00121 // set up kt and cam/aachen definitions 00122 jet_defs[3] = fastjet::JetDefinition(fastjet::kt_algorithm, jet_radius); 00123 jet_defs[4] = fastjet::JetDefinition(fastjet::cambridge_algorithm, 00124 jet_radius); 00125 jet_defs[5] = fastjet::JetDefinition(fastjet::antikt_algorithm, 00126 jet_radius); 00127 00128 // call the example jet-finding routine with each of jet definitions 00129 int index=0; 00130 for (vector<fastjet::JetDefinition>::const_iterator jd_it = jet_def_begin; 00131 jd_it != jet_defs.end(); jd_it++) { 00132 // only run if the plugin is available 00133 if ((index>=3) || (plugins[index] != NULL)) 00134 run_jet_finder(input_particles, *jd_it); 00135 index++; 00136 } 00137 00138 // clean up plugin memory. 00139 if (plugins[0] != NULL) delete plugins[0]; // pxcone may not be defined 00140 if (plugins[1] != NULL) delete plugins[1]; // cdfcones may not be defined 00141 if (plugins[2] != NULL) delete plugins[2]; // siscone may not be defined 00142 }