30 #include "siscone/momentum.h"
31 #include "siscone/siscone.h"
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;
base class for dynamic coordinates management
double f
split/merge threshold
int compute_jets(std::vector< Cmomentum > &_particles, double _radius, double _f, int _n_pass_max=0, double _ptmin=0.0, Esplit_merge_scale _split_merge_scale=SM_pttilde)
compute the jets from a given particle set.
int save_contents(FILE *flux)
save final jets
options for the 'cone' sample
int verbose_flag
do we need to print the help message
int version_flag
do we need to print the version description
double ptmin
minimal pT for jet candidates
std::vector< std::vector< Cmomentum > > protocones_list
list of protocones found pass-by-pass
siscone::Esplit_merge_scale SM_var
variable for split-merge
int npass
number of passes (0 for )
int print_version()
print program version
char * ev_name
event to read
int parse_options(int argc, char **argv)
parse oprions
int help_flag
do we need to print the help message
int print_help()
print the help message
final class: gather everything to compute the jet contents.
int N_stop
maximum number of particle