106void Dnn3piCylinder::RemoveAndAddPoints(
const vector<int> & indices_to_remove,
107 const vector<EtaPhi> & points_to_add,
108 vector<int> & indices_added,
109 vector<int> & indices_of_updated_neighbours) {
114 vector<int> plane_indices_to_remove;
115 for (
unsigned int i=0; i < indices_to_remove.size(); i++) {
116 MirrorVertexInfo * mvi;
117 mvi = & _mirror_info[indices_to_remove[i]];
118 plane_indices_to_remove.push_back(mvi->main_index);
119 if (mvi->mirror_index != INEXISTENT_VERTEX) {
120 plane_indices_to_remove.push_back(mvi->mirror_index);
126 vector<EtaPhi> plane_points_to_add;
127 indices_added.clear();
128 for (
unsigned int i=0; i < points_to_add.size(); i++) {
129 indices_added.push_back(_mirror_info.size());
130 _RegisterCylinderPoint(points_to_add[i], plane_points_to_add);
137 vector<int> updated_plane_neighbours, plane_indices_added;
138 _DNN->RemoveAndAddPoints(plane_indices_to_remove, plane_points_to_add,
139 plane_indices_added, updated_plane_neighbours);
145 for (i=0; i < updated_plane_neighbours.size(); i++) {
147 _cylinder_index_of_plane_vertex[updated_plane_neighbours[i]]);}
150 indices_of_updated_neighbours.clear();
151 for (set<int>::iterator iter = index_set.begin();
152 iter != index_set.end(); iter++) {
153 indices_of_updated_neighbours.push_back(*iter);