00001 //STARTHEADER 00002 // $Id: cdfmidpoint_example.cc 385 2006-12-12 17:00:33Z 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 cdf midpoint plugin 00033 // 00034 // Compile it with: make cdfmidpoint_example 00035 // run it with : ./cdfmidpoint_example < data/single-event.dat 00036 // 00037 //---------------------------------------------------------------------- 00038 00039 #include "run_jet_finder.hh" 00040 #include "CDFMidPointPlugin.hh" 00041 #include<vector> 00042 #include<iostream> 00043 using namespace std; 00044 00045 int main(int argc, char** argv) { 00046 00047 00048 // define a plugin pointer 00049 fastjet::JetDefinition::Plugin * plugin; 00050 00051 // allocate a new plugin 00052 double seed_threshold = 1.0; 00053 double cone_radius = 0.7; 00054 double cone_area_fraction = 1.0; 00055 int max_pair_size = 2; 00056 int max_iterations = 100; 00057 double overlap_threshold = 0.5; 00058 plugin = new fastjet::CDFMidPointPlugin (seed_threshold, cone_radius, 00059 cone_area_fraction, max_pair_size, 00060 max_iterations, overlap_threshold); 00061 00062 // create a jet-definition based on the plugin 00063 fastjet::JetDefinition jet_def(plugin); 00064 00065 // call the example routine that runs with arbitrary jet definitions 00066 vector<fastjet::PseudoJet> particles; 00067 read_input_particles(cin, particles); 00068 run_jet_finder(particles, jet_def); 00069 00070 // clean up plugin memory. 00071 delete plugin; 00072 }