many_algs_example.cc

Go to the documentation of this file.
00001 //STARTHEADER
00002 // $Id: many_algs_example.cc 947 2007-11-09 21:06:32Z soyez $
00003 //
00004 // Copyright (c) 2005-2006, Matteo Cacciari and Gavin Salam
00005 //
00006 //----------------------------------------------------------------------
00007 // This file is part of FastJet.
00008 //
00009 //  FastJet is free software; you can redistribute it and/or modify
00010 //  it under the terms of the GNU General Public License as published by
00011 //  the Free Software Foundation; either version 2 of the License, or
00012 //  (at your option) any later version.
00013 //
00014 //  The algorithms that underlie FastJet have required considerable
00015 //  development and are described in hep-ph/0512210. If you use
00016 //  FastJet as part of work towards a scientific publication, please
00017 //  include a citation to the FastJet paper.
00018 //
00019 //  FastJet is distributed in the hope that it will be useful,
00020 //  but WITHOUT ANY WARRANTY; without even the implied warranty of
00021 //  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
00022 //  GNU General Public License for more details.
00023 //
00024 //  You should have received a copy of the GNU General Public License
00025 //  along with FastJet; if not, write to the Free Software
00026 //  Foundation, Inc.:
00027 //      59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
00028 //----------------------------------------------------------------------
00029 //ENDHEADER
00030 
00031 //----------------------------------------------------------------------
00032 // example program to illustrate the use of multiple jet-definitions,
00033 // some of which are plugin based.
00034 //
00035 // Compile it with: make many_algs_example
00036 // run it with    : ./many_algs_example < data/single-event.dat
00037 //----------------------------------------------------------------------
00038 
00039 #include "fastjet/config.h"
00040 #include "run_jet_finder.hh"
00041 
00042 // get info on how fastjet was configured
00043 #include "fastjet/config.h"
00044 
00045 #ifdef ENABLE_PLUGIN_SISCONE
00046 #  include "fastjet/SISConePlugin.hh"
00047 #endif
00048 #ifdef ENABLE_PLUGIN_CDFCONES
00049 #  include "fastjet/CDFMidPointPlugin.hh"
00050 #  include "fastjet/CDFJetCluPlugin.hh"
00051 #endif
00052 #ifdef ENABLE_PLUGIN_PXCONE
00053 #  include "fastjet/PxConePlugin.hh"
00054 #endif
00055 
00056 #include<vector>
00057 #include<iostream>
00058 #include "fastjet/PseudoJet.hh"
00059 
00060 using namespace std;
00061 
00062 int main(int argc, char** argv) {
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:02 2007 for fastjet by  doxygen 1.5.2