#include "fastjet/config.h"#include "run_jet_finder.hh"#include "fastjet/SISConePlugin.hh"#include "fastjet/CDFMidPointPlugin.hh"#include "fastjet/CDFJetCluPlugin.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) |
|
||||||||||||
|
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 }
|
1.4.2