00001 //STARTHEADER 00002 // $Id: pxcone_example.cc 380 2006-12-07 09:25:46Z salam $ 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 the pxcone plugin 00033 // 00034 // Compile it with: make pxcone_example 00035 // run it with : ./pxcone_example < data/single-event.dat 00036 // 00037 // When compiling it make sure that you've modified ../../Makefile 00038 // appropriately to get the relevant fortran libraries 00039 // 00040 //---------------------------------------------------------------------- 00041 00042 #include "run_jet_finder.hh" 00043 #include "SISConePlugin.hh" 00044 #include<vector> 00045 #include<iostream> 00046 using namespace std; 00047 00048 int main(int argc, char** argv) { 00049 00050 00051 // define a plugin point 00052 fastjet::JetDefinition::Plugin * plugin; 00053 00054 // allocate a new plugin 00055 double cone_radius = 0.7; 00056 double overlap_threshold = 0.5; 00057 plugin = new fastjet::SISConePlugin (cone_radius, overlap_threshold); 00058 00059 // create a jet-definition based on the plugin 00060 fastjet::JetDefinition jet_def(plugin); 00061 00062 // call the example routine that runs with arbitrary jet definitions 00063 vector<fastjet::PseudoJet> particles; 00064 read_input_particles(cin, particles); 00065 run_jet_finder(particles, jet_def); 00066 00067 // clean up plugin memory. 00068 delete plugin; 00069 }