|
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
The
SISCone project has been developed by
Gavin Salam and
Gregory Soyez
Documentation generated on Tue Jun 20 2023 18:08:37 for SISCone by
Doxygen 1.9.4
|