30#include "siscone/momentum.h"
31#include "siscone/siscone.h"
38using namespace siscone;
41 vector<Cmomentum> particles;
50 flux = fopen(
"events/single-event.dat",
"r");
52 cerr <<
"cannot read event" << endl;
57 while (fgets(fline, 512, flux)!=NULL){
59 if (sscanf(fline,
"%le%le%le%le", &px, &py, &pz, &E)==4){
60 particles.push_back(
Cmomentum(px, py, pz, E));
63 cout <<
"error in reading event file Giving up." << endl;
74 cout <<
" " << i <<
" jets found in multi-pass run" << endl;
78 cout <<
" " << i <<
" jets found with alternative f" << endl;
82 cout <<
" " << i <<
" jets found in single-pass run" << endl;
85 vector<Cjet>::iterator it_j;
87 fprintf(stdout,
"# pT eta phi px py pz E \n");
88 for (it_j = siscone.
jets.begin(), i1=0 ;
89 it_j != siscone.
jets.end() ; it_j++, i1++){
90 fprintf(stdout,
"Jet %3d: %10.3f %8.3f %8.3f %10.3f %10.3f %10.3f %10.3f\n",
91 i1, it_j->v.perp(), it_j->v.eta, it_j->v.phi, it_j->v.px, it_j->v.py, it_j->v.pz, it_j->v.E);
base class for dynamic coordinates management
final class: gather everything to compute the jet contents.
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 recompute_jets(double _f, double _ptmin=0.0, Esplit_merge_scale _split_merge_scale=SM_pttilde)
recompute the jets with a different overlap parameter.
std::vector< Cjet > jets
list of jets