many_algs_example.cc File Reference

#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:

run_jet_finder.hhfastjet/PseudoJet.hhfastjet/JetDefinition.hhfastjet/internal/numconsts.hhfastjet/internal/base.hh

Go to the source code of this file.

Functions

int main (int argc, char **argv)


Function Documentation

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 }


Generated on Tue Dec 18 17:05:30 2007 for fastjet by  doxygen 1.5.2