30 #include "fastjet/PseudoJet.hh"
31 #include "fastjet/ClusterSequence.hh"
37 FASTJET_BEGIN_NAMESPACE
48 void ClusterSequence::_really_dumb_cluster () {
52 vector<PseudoJet *> jetsp(_jets.size());
53 vector<int> indices(_jets.size());
55 for (
size_t i = 0; i<_jets.size(); i++) {
56 jetsp[i] = & _jets[i];
60 for (
int n = jetsp.size(); n > 0; n--) {
64 double ymin = jet_scale_for_algorithm(*(jetsp[0]));
66 for (
int i = 0; i < n; i++) {
67 double yiB = jet_scale_for_algorithm(*(jetsp[i]));
69 ymin = yiB; ii = i; jj = -2;}
73 for (
int i = 0; i < n-1; i++) {
74 for (
int j = i+1; j < n; j++) {
76 double y = min(jet_scale_for_algorithm(*(jetsp[i])),
77 jet_scale_for_algorithm(*(jetsp[j])))
78 * jetsp[i]->plain_distance(*jetsp[j])*_invR2;
79 if (y < ymin) {ymin = y; ii = i; jj = j;}
95 int newn = 2*jetsp.size() - n;
99 _do_ij_recombination_step(jetsp[ii]-&_jets[0],
100 jetsp[jj]-&_jets[0], ymin, nn);
103 jetsp[ii] = &_jets[nn];
106 jetsp[jj] = jetsp[n-1];
108 indices[jj] = indices[n-1];
122 _do_iB_recombination_step(jetsp[ii]-&_jets[0], ymin);
124 jetsp[ii] = jetsp[n-1];
125 indices[ii] = indices[n-1];
132 FASTJET_END_NAMESPACE