30 #ifndef DROP_CGAL // in case we do not have the code for CGAL
32 #include "fastjet/internal/Dnn4piCylinder.hh"
35 FASTJET_BEGIN_NAMESPACE
39 Dnn4piCylinder::Dnn4piCylinder(
40 const vector<EtaPhi> & input_points,
const bool & verbose) {
43 vector<EtaPhi> copied_points(input_points.size());
44 for (
unsigned int i=0; i < input_points.size(); i++) {
45 double phi = input_points[i].second;
46 assert(phi >= 0.0 && phi < 2*pi);
47 copied_points[i] = _remap_phi(input_points[i]);
50 if (_verbose) cout <<
"============== Preparing _DNN1" << endl;
51 _DNN1 =
new DnnPlane(input_points, verbose);
52 if (_verbose) cout <<
"============== Preparing _DNN2" << endl;
53 _DNN2 =
new DnnPlane(copied_points, verbose);
59 void Dnn4piCylinder::RemoveAndAddPoints(
const vector<int> & indices_to_remove,
60 const vector<EtaPhi> & points_to_add,
61 vector<int> & indices_added,
62 vector<int> & indices_of_updated_neighbours) {
64 vector<int> indices1, indices2;
66 _DNN1->RemoveAndAddPoints(indices_to_remove,points_to_add,
67 indices_added,indices1);
70 vector<EtaPhi> remapped_points(points_to_add.size());
71 for (
size_t i = 0; i < points_to_add.size(); i++) {
72 remapped_points[i] = _remap_phi(points_to_add[i]);
74 _DNN2->RemoveAndAddPoints(indices_to_remove, remapped_points,
75 indices_added,indices2);
81 for (i=0; i < indices1.size(); i++) {index_set.insert(indices1[i]);}
82 for (i=0; i < indices2.size(); i++) {index_set.insert(indices2[i]);}
84 indices_of_updated_neighbours.clear();
85 for (set<int>::iterator iter = index_set.begin();
86 iter != index_set.end(); iter++) {
87 indices_of_updated_neighbours.push_back(*iter);