30 #include "siscone/momentum.h"
31 #include "siscone/siscone.h"
35 using namespace siscone;
37 int main(
int argc,
char *argv[]){
38 vector<Cmomentum> particles;
65 if (opts.
verbose_flag) cout <<
"reading particles" << endl;
66 flux = fopen(opts.
ev_name,
"r");
68 cerr <<
"cannot read event '" << opts.
ev_name <<
"'" << endl;
69 cerr <<
"specify the event to read using the -e option" << endl;
74 fpart = fopen(
"particles.dat",
"w+");
75 while ((opts.
N_stop!=0) && (fgets(fline, 512, flux)!=NULL)){
77 if (sscanf(fline,
"%le%le%le%le", &px, &py, &pz, &E)==4){
78 particles.push_back(
Cmomentum(px, py, pz, E));
79 fprintf(fpart,
"%e\t%e\n", particles[N].eta, particles[N].phi);
83 cout <<
"error in reading event file Giving up." << endl;
93 cout <<
" working with " << N <<
" particles" << endl;
96 if (opts.
verbose_flag) cout <<
"computing jet contents" << endl;
101 cout <<
" pass " << pass <<
" found " << siscone.
protocones_list[pass].size()
102 <<
" stable cones" << endl;
103 cout <<
" Final result: " << i <<
" jets found" << endl;
108 cout <<
"saving result" << endl;
109 flux = fopen(
"jets.dat",
"w+");
114 cout <<
"bye..." << endl;