src/PartSet/PartSet_Tools.cpp

   1 // Copyright (C) 2014-20xx CEA/DEN, EDF R&D
   2
   3 // File:        PartSet_Tools.cpp
   4 // Created:     28 Apr 2014
   5 // Author:      Natalia ERMOLAEVA
   6
   7 #include <PartSet_Tools.h>
   8 #include <PartSet_Module.h>
   9 #include <PartSet_SketcherMgr.h>
  10
  11 #include <ModelAPI_Data.h>
  12 #include <ModelAPI_AttributeDouble.h>
  13 #include <ModelAPI_AttributeRefList.h>
  14 #include <ModelAPI_Document.h>
  15 #include <ModelAPI_Session.h>
  16 #include <ModelAPI_ResultConstruction.h>
  17 #include <ModelAPI_Events.h>
  18 #include <ModelAPI_Validator.h>
  19
  20 #include <Events_Loop.h>
  21
  22 #include <SketcherPrs_Tools.h>
  23
  24 #include <XGUI_ModuleConnector.h>
  25 #include <XGUI_Displayer.h>
  26 #include <XGUI_Workshop.h>
  27 #include <XGUI_SelectionMgr.h>
  28 #include <XGUI_Selection.h>
  29
  30 #include <GeomDataAPI_Point.h>
  31 #include <GeomDataAPI_Dir.h>
  32 #include <GeomDataAPI_Point2D.h>
  33 #include <GeomAPI_Pln.h>
  34 #include <GeomAPI_Pnt2d.h>
  35 #include <GeomAPI_Pnt.h>
  36 #include <GeomAPI_Edge.h>
  37 #include <GeomAPI_Vertex.h>
  38
  39 #include <GeomAPI_Dir.h>
  40 #include <GeomAPI_XYZ.h>
  41
  42 #include <SketchPlugin_Feature.h>
  43 #include <SketchPlugin_Sketch.h>
  44 #include <SketchPlugin_ConstraintCoincidence.h>
  45 #include <SketchPlugin_ConstraintDistance.h>
  46 #include <SketchPlugin_ConstraintLength.h>
  47 #include <SketchPlugin_ConstraintRadius.h>
  48 #include <SketchPlugin_ConstraintRigid.h>
  49 #include <SketchPlugin_Constraint.h>
  50 #include <SketchPlugin_Circle.h>
  51 #include <SketchPlugin_Arc.h>
  52 #include <SketchPlugin_Line.h>
  53 #include <SketchPlugin_Point.h>
  54
  55 #include <ModuleBase_IWorkshop.h>
  56 #include <ModuleBase_ViewerPrs.h>
  57 #include <ModuleBase_Tools.h>
  58
  59 #include <V3d_View.hxx>
  60 #include <gp_Pln.hxx>
  61 #include <gp_Circ.hxx>
  62 #include <ProjLib.hxx>
  63 #include <ElSLib.hxx>
  64 #include <Geom_Line.hxx>
  65 #include <GeomAPI_ProjectPointOnCurve.hxx>
  66 #include <BRep_Tool.hxx>
  67 #include <TopoDS.hxx>
  68 #include <TopoDS_Edge.hxx>
  69 #include <TopoDS_Vertex.hxx>
  70 #include <AIS_InteractiveObject.hxx>
  71 #include <StdSelect_BRepOwner.hxx>
  72 #include <SelectMgr_IndexedMapOfOwner.hxx>
  73
  74 #ifdef _DEBUG
  75 #include <QDebug>
  76 #endif
  77
  78 const double PRECISION_TOLERANCE = 0.000001;
  79 const int AIS_DEFAULT_WIDTH = 2;
  80
  81 int PartSet_Tools::getAISDefaultWidth()
  82 {
  83   return AIS_DEFAULT_WIDTH;
  84 }
  85
  86 gp_Pnt PartSet_Tools::convertClickToPoint(QPoint thePoint, Handle(V3d_View) theView)
  87 {
  88   if (theView.IsNull())
  89     return gp_Pnt();
  90
  91   V3d_Coordinate XEye, YEye, ZEye, XAt, YAt, ZAt;
  92   theView->Eye(XEye, YEye, ZEye);
  93
  94   theView->At(XAt, YAt, ZAt);
  95   gp_Pnt EyePoint(XEye, YEye, ZEye);
  96   gp_Pnt AtPoint(XAt, YAt, ZAt);
  97
  98   gp_Vec EyeVector(EyePoint, AtPoint);
  99   gp_Dir EyeDir(EyeVector);
 100
 101   gp_Pln PlaneOfTheView = gp_Pln(AtPoint, EyeDir);
 102   Standard_Real X, Y, Z;
 103   theView->Convert(thePoint.x(), thePoint.y(), X, Y, Z);
 104   gp_Pnt ConvertedPoint(X, Y, Z);
 105
 106   gp_Pnt2d ConvertedPointOnPlane = ProjLib::Project(PlaneOfTheView, ConvertedPoint);
 107   gp_Pnt ResultPoint = ElSLib::Value(ConvertedPointOnPlane.X(), ConvertedPointOnPlane.Y(),
 108                                      PlaneOfTheView);
 109   return ResultPoint;
 110 }
 111
 112 void PartSet_Tools::convertTo2D(const gp_Pnt& thePoint, FeaturePtr theSketch,
 113 Handle(V3d_View) theView,
 114                                 double& theX, double& theY)
 115 {
 116   if (!theSketch)
 117     return;
 118
 119   AttributeDoublePtr anAttr;
 120   std::shared_ptr<ModelAPI_Data> aData = theSketch->data();
 121
 122   std::shared_ptr<GeomDataAPI_Point> anOrigin = std::dynamic_pointer_cast<GeomDataAPI_Point>(
 123       aData->attribute(SketchPlugin_Sketch::ORIGIN_ID()));
 124
 125   std::shared_ptr<GeomDataAPI_Dir> aX = std::dynamic_pointer_cast<GeomDataAPI_Dir>(
 126       aData->attribute(SketchPlugin_Sketch::DIRX_ID()));
 127   std::shared_ptr<GeomDataAPI_Dir> aNorm = std::dynamic_pointer_cast<GeomDataAPI_Dir>(
 128       aData->attribute(SketchPlugin_Sketch::NORM_ID()));
 129   std::shared_ptr<GeomAPI_XYZ> anY = aNorm->xyz()->cross(aX->xyz());
 130
 131   gp_Pnt anOriginPnt(anOrigin->x(), anOrigin->y(), anOrigin->z());
 132   gp_Vec aVec(anOriginPnt, thePoint);
 133
 134   if (!theView.IsNull()) {
 135     V3d_Coordinate XEye, YEye, ZEye, XAt, YAt, ZAt;
 136     theView->Eye(XEye, YEye, ZEye);
 137
 138     theView->At(XAt, YAt, ZAt);
 139     gp_Pnt EyePoint(XEye, YEye, ZEye);
 140     gp_Pnt AtPoint(XAt, YAt, ZAt);
 141
 142     gp_Vec anEyeVec(EyePoint, AtPoint);
 143     anEyeVec.Normalize();
 144
 145     std::shared_ptr<GeomDataAPI_Dir> aNormal = std::dynamic_pointer_cast<GeomDataAPI_Dir>(
 146         aData->attribute(SketchPlugin_Sketch::NORM_ID()));
 147     gp_Vec aNormalVec(aNormal->x(), aNormal->y(), aNormal->z());
 148
 149     double aDen = anEyeVec * aNormalVec;
 150     double aLVec = aDen != 0 ? aVec * aNormalVec / aDen : DBL_MAX;
 151
 152     gp_Vec aDeltaVec = anEyeVec * aLVec;
 153     aVec = aVec - aDeltaVec;
 154   }
 155   theX = aVec.X() * aX->x() + aVec.Y() * aX->y() + aVec.Z() * aX->z();
 156   theY = aVec.X() * anY->x() + aVec.Y() * anY->y() + aVec.Z() * anY->z();
 157 }
 158
 159 std::shared_ptr<GeomAPI_Pnt2d> PartSet_Tools::convertTo2D(FeaturePtr theSketch,
 160                                                     const std::shared_ptr<GeomAPI_Pnt>& thePnt)
 161 {
 162   std::shared_ptr<GeomAPI_Pnt2d> aRes;
 163   if (theSketch->getKind() != SketchPlugin_Sketch::ID())
 164     return aRes;
 165   std::shared_ptr<GeomDataAPI_Point> aC = std::dynamic_pointer_cast<GeomDataAPI_Point>(
 166       theSketch->data()->attribute(SketchPlugin_Sketch::ORIGIN_ID()));
 167   std::shared_ptr<GeomDataAPI_Dir> aNorm = std::dynamic_pointer_cast<GeomDataAPI_Dir>(
 168       theSketch->data()->attribute(SketchPlugin_Sketch::NORM_ID()));
 169   std::shared_ptr<GeomDataAPI_Dir> aX = std::dynamic_pointer_cast<GeomDataAPI_Dir>(
 170       theSketch->data()->attribute(SketchPlugin_Sketch::DIRX_ID()));
 171   std::shared_ptr<GeomAPI_Dir> aY(new GeomAPI_Dir(aNorm->dir()->cross(aX->dir())));
 172   return thePnt->to2D(aC->pnt(), aX->dir(), aY);
 173 }
 174
 175
 176 std::shared_ptr<GeomAPI_Pnt> PartSet_Tools::convertTo3D(const double theX, const double theY, FeaturePtr theSketch)
 177 {
 178   std::shared_ptr<ModelAPI_Data> aData = theSketch->data();
 179
 180   std::shared_ptr<GeomDataAPI_Point> aC = std::dynamic_pointer_cast<GeomDataAPI_Point>(
 181       aData->attribute(SketchPlugin_Sketch::ORIGIN_ID()));
 182   std::shared_ptr<GeomDataAPI_Dir> aX = std::dynamic_pointer_cast<GeomDataAPI_Dir>(
 183       aData->attribute(SketchPlugin_Sketch::DIRX_ID()));
 184   std::shared_ptr<GeomDataAPI_Dir> aNorm = std::dynamic_pointer_cast<GeomDataAPI_Dir>(
 185       aData->attribute(SketchPlugin_Sketch::NORM_ID()));
 186   std::shared_ptr<GeomAPI_Dir> aY(new GeomAPI_Dir(aNorm->dir()->cross(aX->dir())));
 187
 188   std::shared_ptr<GeomAPI_Pnt2d> aPnt2d =
 189     std::shared_ptr<GeomAPI_Pnt2d>(new GeomAPI_Pnt2d(theX, theY));
 190
 191   return aPnt2d->to3D(aC->pnt(), aX->dir(), aY);
 192 }
 193
 194 std::shared_ptr<ModelAPI_Document> PartSet_Tools::document()
 195 {
 196   return ModelAPI_Session::get()->moduleDocument();
 197 }
 198
 199 void PartSet_Tools::setFeatureValue(FeaturePtr theFeature, double theValue,
 200                                     const std::string& theAttribute)
 201 {
 202   if (!theFeature)
 203     return;
 204   std::shared_ptr<ModelAPI_Data> aData = theFeature->data();
 205   AttributeDoublePtr anAttribute = std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(
 206       aData->attribute(theAttribute));
 207   if (anAttribute)
 208     anAttribute->setValue(theValue);
 209 }
 210
 211 double PartSet_Tools::featureValue(FeaturePtr theFeature, const std::string& theAttribute,
 212                                    bool& isValid)
 213 {
 214   isValid = false;
 215   double aValue = 0;
 216   if (theFeature) {
 217     std::shared_ptr<ModelAPI_Data> aData = theFeature->data();
 218     AttributeDoublePtr anAttribute = std::dynamic_pointer_cast<ModelAPI_AttributeDouble>(
 219         aData->attribute(theAttribute));
 220     if (anAttribute) {
 221       aValue = anAttribute->value();
 222       isValid = true;
 223     }
 224   }
 225   return aValue;
 226 }
 227
 228 FeaturePtr PartSet_Tools::feature(FeaturePtr theFeature, const std::string& theAttribute,
 229                                   const std::string& theKind)
 230 {
 231   FeaturePtr aFeature;
 232   if (!theFeature)
 233     return aFeature;
 234
 235   std::shared_ptr<ModelAPI_Data> aData = theFeature->data();
 236   std::shared_ptr<ModelAPI_AttributeRefAttr> anAttr = std::dynamic_pointer_cast<
 237       ModelAPI_AttributeRefAttr>(aData->attribute(theAttribute));
 238   if (anAttr) {
 239     aFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(anAttr->object());
 240     if (!theKind.empty() && aFeature && aFeature->getKind() != theKind) {
 241       aFeature = FeaturePtr();
 242     }
 243   }
 244   return aFeature;
 245 }
 246
 247 void PartSet_Tools::createConstraint(CompositeFeaturePtr theSketch,
 248                                      std::shared_ptr<GeomDataAPI_Point2D> thePoint1,
 249                                      std::shared_ptr<GeomDataAPI_Point2D> thePoint2)
 250 {
 251   FeaturePtr aFeature;
 252   if (theSketch) {
 253     aFeature = theSketch->addFeature(SketchPlugin_ConstraintCoincidence::ID());
 254   } else {
 255     std::shared_ptr<ModelAPI_Document> aDoc = document();
 256     aFeature = aDoc->addFeature(SketchPlugin_ConstraintCoincidence::ID());
 257   }
 258
 259   std::shared_ptr<ModelAPI_Data> aData = aFeature->data();
 260
 261   std::shared_ptr<ModelAPI_AttributeRefAttr> aRef1 = std::dynamic_pointer_cast<
 262       ModelAPI_AttributeRefAttr>(aData->attribute(SketchPlugin_Constraint::ENTITY_A()));
 263   aRef1->setAttr(thePoint1);
 264
 265   std::shared_ptr<ModelAPI_AttributeRefAttr> aRef2 = std::dynamic_pointer_cast<
 266       ModelAPI_AttributeRefAttr>(aData->attribute(SketchPlugin_Constraint::ENTITY_B()));
 267   aRef2->setAttr(thePoint2);
 268
 269   // we need to flush created signal in order to coincidence is processed by solver
 270   Events_Loop::loop()->flush(Events_Loop::eventByName(EVENT_OBJECT_CREATED));
 271 }
 272
 273 /*std::shared_ptr<GeomDataAPI_Point2D> PartSet_Tools::
 274   findAttributePoint(CompositeFeaturePtr theSketch, double theX, double theY,
 275   double theTolerance, const QList<FeaturePtr>& theIgnore)
 276 {
 277   std::shared_ptr<GeomAPI_Pnt2d> aClickedPoint = std::shared_ptr<GeomAPI_Pnt2d>(
 278       new GeomAPI_Pnt2d(theX, theY));
 279
 280   std::list<std::shared_ptr<ModelAPI_Attribute> > anAttiributes;
 281   for (int i = 0; i < theSketch->numberOfSubs(); i++) {
 282     FeaturePtr aFeature = theSketch->subFeature(i);
 283     if (!theIgnore.contains(aFeature)) {
 284       anAttiributes = aFeature->data()->attributes(GeomDataAPI_Point2D::typeId());
 285
 286       std::list<std::shared_ptr<ModelAPI_Attribute> >::const_iterator anIt;
 287       for (anIt = anAttiributes.cbegin(); anIt != anAttiributes.cend(); ++anIt) {
 288         std::shared_ptr<GeomDataAPI_Point2D> aCurPoint =
 289           std::dynamic_pointer_cast<GeomDataAPI_Point2D>(*anIt);
 290         double x = aCurPoint->x();
 291         double y = aCurPoint->y();
 292         if (aCurPoint && (aCurPoint->pnt()->distance(aClickedPoint) < theTolerance)) {
 293           return aCurPoint;
 294         }
 295       }
 296     }
 297   }
 298   return std::shared_ptr<GeomDataAPI_Point2D>();
 299 }*/
 300
 301
 302 std::shared_ptr<GeomDataAPI_Point2D> PartSet_Tools::findFirstEqualPointInArgumentFeatures(
 303                   const FeaturePtr& theFeature, const std::shared_ptr<GeomAPI_Pnt2d>& thePoint)
 304 {
 305   std::shared_ptr<GeomDataAPI_Point2D> aFeaturePoint;
 306
 307   // may be feature is not updated yet, execute is not performed and references features
 308   // are not created. Case: rectangle macro feature
 309   ModuleBase_Tools::flushUpdated(theFeature);
 310
 311   std::list<AttributePtr> anAttributes = theFeature->data()->attributes(
 312                                           ModelAPI_AttributeRefList::typeId());
 313   std::list<AttributePtr>::const_iterator anIt = anAttributes.begin(), aLast = anAttributes.end();
 314   for (; anIt != aLast && !aFeaturePoint.get(); anIt++) {
 315     std::shared_ptr<ModelAPI_AttributeRefList> aCurSelList =
 316                                       std::dynamic_pointer_cast<ModelAPI_AttributeRefList>(*anIt);
 317     for (int i = 0, aNb = aCurSelList->size(); i < aNb && !aFeaturePoint.get(); i++) {
 318       ObjectPtr anObject = aCurSelList->object(i);
 319       FeaturePtr aFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(anObject);
 320       if (aFeature.get())
 321         aFeaturePoint = PartSet_Tools::findFirstEqualPoint(aFeature, thePoint);
 322     }
 323   }
 324   return aFeaturePoint;
 325 }
 326
 327 std::shared_ptr<GeomDataAPI_Point2D> PartSet_Tools::findFirstEqualPoint(const FeaturePtr& theFeature,
 328                                                       const std::shared_ptr<GeomAPI_Pnt2d>& thePoint)
 329 {
 330   std::shared_ptr<GeomDataAPI_Point2D> aFPoint;
 331
 332   // find the given point in the feature attributes
 333   std::list<std::shared_ptr<ModelAPI_Attribute> > anAttiributes =
 334                                     theFeature->data()->attributes(GeomDataAPI_Point2D::typeId());
 335   std::list<std::shared_ptr<ModelAPI_Attribute> >::const_iterator anIt = anAttiributes.begin(),
 336       aLast = anAttiributes.end();
 337   ModelAPI_ValidatorsFactory* aValidators = ModelAPI_Session::get()->validators();
 338
 339   for (; anIt != aLast && !aFPoint; anIt++) {
 340     std::shared_ptr<GeomDataAPI_Point2D> aCurPoint =
 341                                              std::dynamic_pointer_cast<GeomDataAPI_Point2D>(*anIt);
 342     if (aCurPoint && aCurPoint->isInitialized() &&
 343         aValidators->isCase(theFeature, aCurPoint->id()) &&
 344         (aCurPoint->pnt()->distance(thePoint) < Precision::Confusion())) {
 345       aFPoint = aCurPoint;
 346       break;
 347     }
 348   }
 349   return aFPoint;
 350 }
 351
 352 void PartSet_Tools::setConstraints(CompositeFeaturePtr theSketch, FeaturePtr theFeature,
 353                                    const std::string& theAttribute, double theClickedX,
 354                                    double theClickedY)
 355 {
 356   if (!theFeature.get())
 357     return;
 358
 359   std::shared_ptr<GeomAPI_Pnt2d> aClickedPoint = std::shared_ptr<GeomAPI_Pnt2d>(
 360       new GeomAPI_Pnt2d(theClickedX, theClickedY));
 361
 362   // find a feature point by the selection mode
 363   std::shared_ptr<GeomDataAPI_Point2D> aFeaturePoint;
 364   if (theFeature->isMacro()) {
 365     // the macro feature will be removed after the operation is stopped, so we need to build
 366     // coicidence to possible sub-features
 367     aFeaturePoint = PartSet_Tools::findFirstEqualPointInArgumentFeatures(theFeature, aClickedPoint);
 368   }
 369   else {
 370     aFeaturePoint = std::dynamic_pointer_cast<
 371         GeomDataAPI_Point2D>(theFeature->data()->attribute(theAttribute));
 372   }
 373   if (!aFeaturePoint)
 374     return;
 375
 376   // get all sketch features. If the point with the given coordinates belong to any sketch feature,
 377   // the constraint is created between the feature point and the found sketch point
 378   std::shared_ptr<ModelAPI_Data> aData = theSketch->data();
 379   std::shared_ptr<ModelAPI_AttributeRefList> aRefList = std::dynamic_pointer_cast<
 380       ModelAPI_AttributeRefList>(aData->attribute(SketchPlugin_Sketch::FEATURES_ID()));
 381
 382   std::list<ObjectPtr> aFeatures = aRefList->list();
 383   std::list<ObjectPtr>::const_iterator anIt = aFeatures.begin(), aLast = aFeatures.end();
 384   std::list<std::shared_ptr<ModelAPI_Attribute> > anAttiributes;
 385   for (; anIt != aLast; anIt++) {
 386     FeaturePtr aFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(*anIt);
 387     if (!aFeature.get() || theFeature == aFeature)
 388       continue;
 389     std::shared_ptr<GeomDataAPI_Point2D> aFPoint = PartSet_Tools::findFirstEqualPoint(aFeature,
 390                                                                                 aClickedPoint);
 391     if (aFPoint)
 392       PartSet_Tools::createConstraint(theSketch, aFPoint, aFeaturePoint);
 393   }
 394 }
 395
 396 std::shared_ptr<GeomAPI_Pln> PartSet_Tools::sketchPlane(CompositeFeaturePtr theSketch)
 397 {
 398   std::shared_ptr<GeomDataAPI_Point> anOrigin = std::dynamic_pointer_cast<GeomDataAPI_Point>(
 399       theSketch->data()->attribute(SketchPlugin_Sketch::ORIGIN_ID()));
 400   std::shared_ptr<GeomDataAPI_Dir> aNorm = std::dynamic_pointer_cast<GeomDataAPI_Dir>(
 401       theSketch->data()->attribute(SketchPlugin_Sketch::NORM_ID()));
 402
 403   if (!anOrigin || !aNorm)
 404     return std::shared_ptr<GeomAPI_Pln>();
 405
 406   return std::shared_ptr<GeomAPI_Pln>(new GeomAPI_Pln(anOrigin->pnt(), aNorm->dir()));
 407 }
 408
 409 std::shared_ptr<GeomAPI_Pnt> PartSet_Tools::point3D(std::shared_ptr<GeomAPI_Pnt2d> thePoint2D,
 410                                                       CompositeFeaturePtr theSketch)
 411 {
 412   std::shared_ptr<GeomAPI_Pnt> aPoint;
 413   if (!theSketch || !thePoint2D)
 414     return aPoint;
 415
 416   DataPtr aData = theSketch->data();
 417   std::shared_ptr<GeomDataAPI_Point> aC = std::dynamic_pointer_cast<GeomDataAPI_Point>(
 418       aData->attribute(SketchPlugin_Sketch::ORIGIN_ID()));
 419   std::shared_ptr<GeomDataAPI_Dir> aX = std::dynamic_pointer_cast<GeomDataAPI_Dir>(
 420       aData->attribute(SketchPlugin_Sketch::DIRX_ID()));
 421   std::shared_ptr<GeomDataAPI_Dir> aNorm = std::dynamic_pointer_cast<GeomDataAPI_Dir>(
 422       aData->attribute(SketchPlugin_Sketch::NORM_ID()));
 423   std::shared_ptr<GeomAPI_Dir> aY(new GeomAPI_Dir(aNorm->dir()->cross(aX->dir())));
 424
 425   return thePoint2D->to3D(aC->pnt(), aX->dir(), aY);
 426 }
 427
 428 ResultPtr PartSet_Tools::findFixedObjectByExternal(const TopoDS_Shape& theShape,
 429                                                    const ObjectPtr& theObject,
 430                                                    CompositeFeaturePtr theSketch)
 431 {
 432   ResultPtr aResult;
 433   if (theShape.ShapeType() == TopAbs_EDGE) {
 434     // Check that we already have such external edge
 435     std::shared_ptr<GeomAPI_Edge> aInEdge = std::shared_ptr<GeomAPI_Edge>(new GeomAPI_Edge());
 436     aInEdge->setImpl(new TopoDS_Shape(theShape));
 437     aResult = findExternalEdge(theSketch, aInEdge);
 438   }
 439   if (theShape.ShapeType() == TopAbs_VERTEX) {
 440     std::shared_ptr<GeomAPI_Vertex> aInVert = std::shared_ptr<GeomAPI_Vertex>(new GeomAPI_Vertex());
 441     aInVert->setImpl(new TopoDS_Shape(theShape));
 442     aResult = findExternalVertex(theSketch, aInVert);
 443   }
 444   return aResult;
 445 }
 446
 447 ResultPtr PartSet_Tools::createFixedObjectByExternal(const TopoDS_Shape& theShape,
 448                                                      const ObjectPtr& theObject,
 449                                                      CompositeFeaturePtr theSketch,
 450                                                      const bool theTemporary)
 451 {
 452   if (theShape.ShapeType() == TopAbs_EDGE) {
 453     Standard_Real aStart, aEnd;
 454     Handle(V3d_View) aNullView;
 455     FeaturePtr aMyFeature;
 456
 457     Handle(Geom_Curve) aCurve = BRep_Tool::Curve(TopoDS::Edge(theShape), aStart, aEnd);
 458     GeomAdaptor_Curve aAdaptor(aCurve);
 459     std::shared_ptr<GeomAPI_Edge> anEdge = std::shared_ptr<GeomAPI_Edge>(new GeomAPI_Edge);
 460     anEdge->setImpl(new TopoDS_Shape(theShape));
 461     if (aAdaptor.GetType() == GeomAbs_Line) {
 462       // Create line
 463       aMyFeature = theSketch->addFeature(SketchPlugin_Line::ID());
 464       if (!theObject.get()) {
 465         // There is no selected result
 466         std::shared_ptr<GeomAPI_Pnt> aPnt1 = anEdge->firstPoint();
 467         std::shared_ptr<GeomAPI_Pnt> aPnt2 = anEdge->lastPoint();
 468         std::shared_ptr<GeomAPI_Pnt2d> aPnt2d1 = convertTo2D(theSketch, aPnt1);
 469         std::shared_ptr<GeomAPI_Pnt2d> aPnt2d2 = convertTo2D(theSketch, aPnt2);
 470
 471         std::shared_ptr<ModelAPI_Data> aData = aMyFeature->data();
 472         std::shared_ptr<GeomDataAPI_Point2D> aPoint1 =
 473           std::dynamic_pointer_cast<GeomDataAPI_Point2D>(aData->attribute(SketchPlugin_Line::START_ID()));
 474         std::shared_ptr<GeomDataAPI_Point2D> aPoint2 =
 475           std::dynamic_pointer_cast<GeomDataAPI_Point2D>(aData->attribute(SketchPlugin_Line::END_ID()));
 476
 477         aPoint1->setValue(aPnt2d1);
 478         aPoint2->setValue(aPnt2d2);
 479
 480         // If this is an axis then its name has to be changed correspondently
 481         std::string aSuffix = "";
 482         bool aXdir = fabs(aPnt1->x() - aPnt2->x()) > Precision::Confusion();
 483         bool aYdir = fabs(aPnt1->y() - aPnt2->y()) > Precision::Confusion();
 484         bool aZdir = fabs(aPnt1->z() - aPnt2->z()) > Precision::Confusion();
 485         if (aXdir && (!aYdir) && (!aZdir))
 486           aSuffix = "X";
 487         else if ((!aXdir) && aYdir && (!aZdir))
 488           aSuffix = "Y";
 489         else if ((!aXdir) && (!aYdir) && aZdir)
 490           aSuffix = "Z";
 491         if (aSuffix.length() > 0)
 492           aData->setName("Axis_" + aSuffix);
 493         aMyFeature->execute();
 494
 495       }
 496     } else if (aAdaptor.GetType() == GeomAbs_Circle) {
 497       if (anEdge->isArc()) {
 498         // Create arc
 499         aMyFeature = theSketch->addFeature(SketchPlugin_Arc::ID());
 500       }
 501       else {
 502         // Create circle
 503         aMyFeature = theSketch->addFeature(SketchPlugin_Circle::ID());
 504       }
 505     }
 506     if (aMyFeature) {
 507       DataPtr aData = aMyFeature->data();
 508       AttributeSelectionPtr anAttr =
 509         std::dynamic_pointer_cast<ModelAPI_AttributeSelection>
 510         (aData->attribute(SketchPlugin_SketchEntity::EXTERNAL_ID()));
 511
 512       ResultPtr aRes = std::dynamic_pointer_cast<ModelAPI_Result>(theObject);
 513       // selection shape has no result owner => the trihedron axis
 514       // TODO: make reference to the real axes when they are implemented in the initialization plugin
 515       if (!aRes.get()) {
 516         ObjectPtr aPointObj = ModelAPI_Session::get()->moduleDocument()->objectByName(
 517           ModelAPI_ResultConstruction::group(), "Origin");
 518         if (aPointObj.get()) { // if initialization plugin performed well
 519           aRes = std::dynamic_pointer_cast<ModelAPI_Result>(aPointObj);
 520         }
 521       }
 522       if (!aRes.get()) {
 523         aRes = aMyFeature->firstResult();
 524       }
 525       if (anAttr.get() && aRes.get()) {
 526         std::shared_ptr<GeomAPI_Shape> anEdge(new GeomAPI_Shape);
 527         anEdge->setImpl(new TopoDS_Shape(theShape));
 528
 529         anAttr->setValue(aRes, anEdge);
 530         //if (!theTemporary) {
 531           aMyFeature->execute();
 532
 533           // fix this edge
 534           FeaturePtr aFix = theSketch->addFeature(SketchPlugin_ConstraintRigid::ID());
 535           aFix->data()->refattr(SketchPlugin_Constraint::ENTITY_A())->
 536             setObject(aMyFeature->lastResult());
 537         //}
 538         return aMyFeature->lastResult();
 539       }
 540     }
 541   }
 542   if (theShape.ShapeType() == TopAbs_VERTEX) {
 543     FeaturePtr aMyFeature = theSketch->addFeature(SketchPlugin_Point::ID());
 544
 545     if (aMyFeature) {
 546       DataPtr aData = aMyFeature->data();
 547       AttributeSelectionPtr anAttr =
 548         std::dynamic_pointer_cast<ModelAPI_AttributeSelection>
 549         (aData->attribute(SketchPlugin_SketchEntity::EXTERNAL_ID()));
 550
 551       ResultPtr aRes = std::dynamic_pointer_cast<ModelAPI_Result>(theObject);
 552       // if there is no object, it means that this is the origin point: search it in the module document
 553       if (!aRes.get()) {
 554         ObjectPtr aPointObj = ModelAPI_Session::get()->moduleDocument()->objectByName(
 555           ModelAPI_ResultConstruction::group(), "Origin");
 556         if (aPointObj.get()) { // if initialization plugin performed well
 557           aRes = std::dynamic_pointer_cast<ModelAPI_Result>(aPointObj);
 558         }
 559       }
 560       // reference to itself with name "Origin" (but this may cause the infinitive cycling)
 561       if (!aRes.get()) {
 562         // If the point is selected not from Result object
 563         std::shared_ptr<GeomAPI_Shape> aShape =
 564           std::shared_ptr<GeomAPI_Shape>(new GeomAPI_Shape());
 565         aShape->setImpl(new TopoDS_Shape(theShape));
 566
 567         std::shared_ptr<GeomAPI_Vertex> aVertex =
 568           std::shared_ptr<GeomAPI_Vertex>(new GeomAPI_Vertex(aShape));
 569         std::shared_ptr<GeomAPI_Pnt> aPnt = aVertex->point();
 570
 571         std::shared_ptr<GeomAPI_Pnt2d> aPnt2d = convertTo2D(theSketch, aPnt);
 572         std::shared_ptr<GeomDataAPI_Point2D> aPoint =
 573           std::dynamic_pointer_cast<GeomDataAPI_Point2D>(aData->attribute(SketchPlugin_Point::COORD_ID()));
 574         aPoint->setValue(aPnt2d);
 575         if ((aPnt->x() < Precision::Confusion()) &&
 576             (aPnt->y() < Precision::Confusion()) &&
 577             (aPnt->z() < Precision::Confusion()))
 578           aData->setName("Origin");
 579
 580         aMyFeature->execute();
 581         aRes = aMyFeature->firstResult();
 582       }
 583       if (anAttr.get() && aRes.get()) {
 584         std::shared_ptr<GeomAPI_Shape> aVert(new GeomAPI_Shape);
 585         aVert->setImpl(new TopoDS_Shape(theShape));
 586
 587         anAttr->setValue(aRes, aVert);
 588         //if (theTemporary) {
 589           aMyFeature->execute();
 590
 591           // fix this edge
 592           FeaturePtr aFix = theSketch->addFeature(SketchPlugin_ConstraintRigid::ID());
 593           aFix->data()->refattr(SketchPlugin_Constraint::ENTITY_A())->
 594             setObject(aMyFeature->lastResult());
 595         //}
 596         return aMyFeature->lastResult();
 597       }
 598     }
 599   }
 600   return ResultPtr();
 601 }
 602
 603 bool PartSet_Tools::isContainPresentation(const QList<ModuleBase_ViewerPrsPtr>& theSelected,
 604                                           const ModuleBase_ViewerPrsPtr& thePrs)
 605 {
 606   foreach (ModuleBase_ViewerPrsPtr aPrs, theSelected) {
 607     if (aPrs->object() == thePrs->object())
 608       return true;
 609   }
 610   return false;
 611 }
 612
 613 ResultPtr PartSet_Tools::findExternalEdge(CompositeFeaturePtr theSketch, std::shared_ptr<GeomAPI_Edge> theEdge)
 614 {
 615   for (int i = 0; i < theSketch->numberOfSubs(); i++) {
 616     FeaturePtr aFeature = theSketch->subFeature(i);
 617     std::shared_ptr<SketchPlugin_Feature> aSketchFea =
 618       std::dynamic_pointer_cast<SketchPlugin_Feature>(aFeature);
 619     if (aSketchFea) {
 620       if (aSketchFea->isExternal()) {
 621         std::list<ResultPtr> aResults = aSketchFea->results();
 622         std::list<ResultPtr>::const_iterator aIt;
 623         for (aIt = aResults.cbegin(); aIt != aResults.cend(); ++aIt) {
 624           ResultConstructionPtr aRes =
 625             std::dynamic_pointer_cast<ModelAPI_ResultConstruction>(*aIt);
 626           if (aRes) {
 627             std::shared_ptr<GeomAPI_Shape> aShape = aRes->shape();
 628             if (aShape) {
 629               if (theEdge->isEqual(aShape))
 630                 return aRes;
 631             }
 632           }
 633         }
 634       }
 635     }
 636   }
 637   return ResultPtr();
 638 }
 639
 640
 641 ResultPtr PartSet_Tools::findExternalVertex(CompositeFeaturePtr theSketch, std::shared_ptr<GeomAPI_Vertex> theVert)
 642 {
 643   for (int i = 0; i < theSketch->numberOfSubs(); i++) {
 644     FeaturePtr aFeature = theSketch->subFeature(i);
 645     std::shared_ptr<SketchPlugin_Feature> aSketchFea =
 646       std::dynamic_pointer_cast<SketchPlugin_Feature>(aFeature);
 647     if (aSketchFea) {
 648       if (aSketchFea->isExternal()) {
 649         std::list<ResultPtr> aResults = aSketchFea->results();
 650         std::list<ResultPtr>::const_iterator aIt;
 651         for (aIt = aResults.cbegin(); aIt != aResults.cend(); ++aIt) {
 652           ResultConstructionPtr aRes =
 653             std::dynamic_pointer_cast<ModelAPI_ResultConstruction>(*aIt);
 654           if (aRes) {
 655             std::shared_ptr<GeomAPI_Shape> aShape = aRes->shape();
 656             if (aShape) {
 657               if (theVert->isEqual(aShape))
 658                 return aRes;
 659             }
 660           }
 661         }
 662       }
 663     }
 664   }
 665   return ResultPtr();
 666 }
 667
 668
 669 bool PartSet_Tools::hasVertexShape(const ModuleBase_ViewerPrsPtr& thePrs, FeaturePtr theSketch,
 670                                    Handle_V3d_View theView, double& theX, double& theY)
 671 {
 672   bool aHasVertex = false;
 673
 674   const GeomShapePtr& aShape = thePrs->shape();
 675   if (aShape.get() && !aShape->isNull() && aShape->shapeType() == GeomAPI_Shape::VERTEX)
 676   {
 677     const TopoDS_Shape& aTDShape = aShape->impl<TopoDS_Shape>();
 678     const TopoDS_Vertex& aVertex = TopoDS::Vertex(aTDShape);
 679     if (!aVertex.IsNull())
 680     {
 681       gp_Pnt aPoint = BRep_Tool::Pnt(aVertex);
 682       PartSet_Tools::convertTo2D(aPoint, theSketch, theView, theX, theY);
 683       aHasVertex = true;
 684     }
 685   }
 686
 687   return aHasVertex;
 688 }
 689
 690 GeomShapePtr PartSet_Tools::findShapeBy2DPoint(const AttributePtr& theAttribute,
 691                                                        ModuleBase_IWorkshop* theWorkshop)
 692 {
 693   GeomShapePtr aShape;
 694   XGUI_ModuleConnector* aConnector = dynamic_cast<XGUI_ModuleConnector*>(theWorkshop);
 695   XGUI_Displayer* aDisplayer = aConnector->workshop()->displayer();
 696
 697   // 2. find visualized vertices of the attribute and if the attribute of the vertex is
 698   // the same, return it
 699   FeaturePtr anAttributeFeature = ModelAPI_Feature::feature(theAttribute->owner());
 700   // 2.1 get visualized results of the feature
 701   const std::list<ResultPtr>& aResList = anAttributeFeature->results();
 702   std::list<ResultPtr>::const_iterator anIt = aResList.begin(), aLast = aResList.end();
 703   for (; anIt != aLast; anIt++) {
 704     AISObjectPtr aAISObj = aDisplayer->getAISObject(*anIt);
 705     if (aAISObj.get() != NULL) {
 706       Handle(AIS_InteractiveObject) anAISIO = aAISObj->impl<Handle(AIS_InteractiveObject)>();
 707       // 2.2 find selected owners of a visualizedd object
 708       SelectMgr_IndexedMapOfOwner aSelectedOwners;
 709       aConnector->workshop()->selector()->selection()->entityOwners(anAISIO, aSelectedOwners);
 710       for (Standard_Integer i = 1, n = aSelectedOwners.Extent(); i <= n; i++) {
 711         Handle(SelectMgr_EntityOwner) anOwner = aSelectedOwners(i);
 712         if (!anOwner.IsNull()) {
 713           Handle(StdSelect_BRepOwner) aBRepOwner = Handle(StdSelect_BRepOwner)::DownCast(anOwner);
 714           if (!aBRepOwner.IsNull() && aBRepOwner->HasShape()) {
 715             const TopoDS_Shape& aBRepShape = aBRepOwner->Shape();
 716             if (aBRepShape.ShapeType() == TopAbs_VERTEX) {
 717               // 2.3 if the owner is vertex and an attribute of the vertex is equal to the initial
 718               // attribute, returns the shape
 719               PartSet_Module* aModule = dynamic_cast<PartSet_Module*>(theWorkshop->module());
 720               PartSet_SketcherMgr* aSketchMgr = aModule->sketchMgr();
 721               AttributePtr aPntAttr = PartSet_Tools::findAttributeBy2dPoint(anAttributeFeature,
 722                                                         aBRepShape, aSketchMgr->activeSketch());
 723               if (aPntAttr.get() != NULL && aPntAttr == theAttribute) {
 724                 aShape = std::shared_ptr<GeomAPI_Shape>(new GeomAPI_Shape);
 725                 aShape->setImpl(new TopoDS_Shape(aBRepShape));
 726                 break;
 727               }
 728             }
 729           }
 730         }
 731       }
 732     }
 733   }
 734   return aShape;
 735 }
 736
 737 std::shared_ptr<GeomAPI_Pnt2d> PartSet_Tools::getPoint(std::shared_ptr<ModelAPI_Feature>& theFeature,
 738                                                        const std::string& theAttribute)
 739 {
 740   std::shared_ptr<GeomDataAPI_Point2D> aPointAttr;
 741
 742   if (!theFeature->data())
 743     return std::shared_ptr<GeomAPI_Pnt2d>();
 744
 745   FeaturePtr aFeature;
 746   std::shared_ptr<ModelAPI_AttributeRefAttr> anAttr = std::dynamic_pointer_cast<
 747       ModelAPI_AttributeRefAttr>(theFeature->data()->attribute(theAttribute));
 748   if (!anAttr)
 749     return std::shared_ptr<GeomAPI_Pnt2d>();
 750
 751   aFeature = ModelAPI_Feature::feature(anAttr->object());
 752
 753   if (aFeature && aFeature->getKind() == SketchPlugin_Point::ID())
 754     aPointAttr = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(
 755         aFeature->data()->attribute(SketchPlugin_Point::COORD_ID()));
 756
 757   else if (anAttr->attr()) {
 758     aPointAttr = std::dynamic_pointer_cast<GeomDataAPI_Point2D>(anAttr->attr());
 759   }
 760   if (aPointAttr.get() != NULL)
 761     return aPointAttr->pnt();
 762   return std::shared_ptr<GeomAPI_Pnt2d>();
 763 }
 764
 765 FeaturePtr findFirstCoincidenceByData(const DataPtr& theData, std::shared_ptr<GeomAPI_Pnt2d> thePoint)
 766 {
 767   FeaturePtr aCoincident;
 768
 769   const std::set<AttributePtr>& aRefsList = theData->refsToMe();
 770   std::set<AttributePtr>::const_iterator aIt;
 771   for (aIt = aRefsList.cbegin(); aIt != aRefsList.cend(); ++aIt) {
 772     std::shared_ptr<ModelAPI_Attribute> aAttr = (*aIt);
 773     FeaturePtr aConstrFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(aAttr->owner());
 774     if (aConstrFeature->getKind() == SketchPlugin_ConstraintCoincidence::ID()) {
 775       std::shared_ptr<GeomAPI_Pnt2d> a2dPnt =
 776         PartSet_Tools::getPoint(aConstrFeature, SketchPlugin_ConstraintCoincidence::ENTITY_A());
 777       if (a2dPnt.get() && thePoint->isEqual(a2dPnt)) {
 778         aCoincident = aConstrFeature;
 779         break;
 780       } else {
 781         a2dPnt = PartSet_Tools::getPoint(aConstrFeature,
 782                                           SketchPlugin_ConstraintCoincidence::ENTITY_B());
 783         if (a2dPnt.get() && thePoint->isEqual(a2dPnt)) {
 784           aCoincident = aConstrFeature;
 785           break;
 786         }
 787       }
 788     }
 789   }
 790   return aCoincident;
 791 }
 792
 793 FeaturePtr PartSet_Tools::findFirstCoincidence(const FeaturePtr& theFeature,
 794                                                std::shared_ptr<GeomAPI_Pnt2d> thePoint)
 795 {
 796   FeaturePtr aCoincident;
 797   if (theFeature.get() == NULL)
 798     return aCoincident;
 799
 800   const std::set<AttributePtr>& aRefsList = theFeature->data()->refsToMe();
 801   std::set<AttributePtr>::const_iterator aIt;
 802   for (aIt = aRefsList.cbegin(); aIt != aRefsList.cend(); ++aIt) {
 803     std::shared_ptr<ModelAPI_Attribute> aAttr = (*aIt);
 804     FeaturePtr aConstrFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(aAttr->owner());
 805     if (aConstrFeature->getKind() == SketchPlugin_ConstraintCoincidence::ID()) {
 806       std::shared_ptr<GeomAPI_Pnt2d> a2dPnt =
 807         PartSet_Tools::getPoint(aConstrFeature, SketchPlugin_ConstraintCoincidence::ENTITY_A());
 808       if (a2dPnt.get() && thePoint->isEqual(a2dPnt)) {
 809         aCoincident = aConstrFeature;
 810         break;
 811       } else {
 812         a2dPnt = PartSet_Tools::getPoint(aConstrFeature,
 813                                           SketchPlugin_ConstraintCoincidence::ENTITY_B());
 814         if (a2dPnt.get() && thePoint->isEqual(a2dPnt)) {
 815           aCoincident = aConstrFeature;
 816           break;
 817         }
 818       }
 819     }
 820   }
 821   /// Find by result
 822   if (!aCoincident.get()) {
 823     std::list<ResultPtr> aResults = theFeature->results();
 824     std::list<ResultPtr>::const_iterator aIt;
 825     for (aIt = aResults.cbegin(); aIt != aResults.cend(); ++aIt) {
 826       ResultPtr aResult = *aIt;
 827       aCoincident = findFirstCoincidenceByData(aResult->data(), thePoint);
 828       if (aCoincident.get())
 829         break;
 830     }
 831   }
 832   return aCoincident;
 833 }
 834
 835 void PartSet_Tools::findCoincidences(FeaturePtr theStartCoin, QList<FeaturePtr>& theList,
 836                                      QList<FeaturePtr>& theCoincidencies,
 837                                      std::string theAttr)
 838 {
 839   std::shared_ptr<GeomAPI_Pnt2d> aOrig = getCoincedencePoint(theStartCoin);
 840   if (aOrig.get() == NULL)
 841     return;
 842
 843   AttributeRefAttrPtr aPnt = theStartCoin->refattr(theAttr);
 844   if (!aPnt)
 845     return;
 846   ObjectPtr aObj = aPnt->object();
 847   FeaturePtr aFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(aObj);
 848   if (aFeature.get()) {
 849     if (!theList.contains(aFeature)) {
 850       theList.append(aFeature);
 851       theCoincidencies.append(theStartCoin);
 852       const std::set<AttributePtr>& aRefsList = aFeature->data()->refsToMe();
 853       std::set<AttributePtr>::const_iterator aIt;
 854       for (aIt = aRefsList.cbegin(); aIt != aRefsList.cend(); ++aIt) {
 855         std::shared_ptr<ModelAPI_Attribute> aAttr = (*aIt);
 856         FeaturePtr aConstrFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(aAttr->owner());
 857         if (aConstrFeature->getKind() == SketchPlugin_ConstraintCoincidence::ID()) {
 858           if (!theCoincidencies.contains(aConstrFeature)) {
 859             std::shared_ptr<GeomAPI_Pnt2d> aPnt = getCoincedencePoint(aConstrFeature);
 860             if (aPnt.get() && aOrig->isEqual(aPnt)) {
 861               findCoincidences(aConstrFeature, theList, theCoincidencies,
 862                 SketchPlugin_ConstraintCoincidence::ENTITY_A());
 863               findCoincidences(aConstrFeature, theList, theCoincidencies,
 864                 SketchPlugin_ConstraintCoincidence::ENTITY_B());
 865             }
 866           }
 867         }
 868       }
 869     }
 870   } else {
 871     // Find by Results
 872     ResultConstructionPtr aResult = std::dynamic_pointer_cast<ModelAPI_ResultConstruction>(aObj);
 873     if (aResult.get()) {
 874       FeaturePtr aFeature = ModelAPI_Feature::feature(aPnt->object());
 875       if (!theList.contains(aFeature))
 876         theList.append(aFeature);
 877       theCoincidencies.append(theStartCoin);
 878
 879       const std::set<AttributePtr>& aRefsList = aResult->data()->refsToMe();
 880       std::set<AttributePtr>::const_iterator aIt;
 881       for (aIt = aRefsList.cbegin(); aIt != aRefsList.cend(); ++aIt) {
 882         std::shared_ptr<ModelAPI_Attribute> aAttr = (*aIt);
 883         FeaturePtr aConstrFeature = std::dynamic_pointer_cast<ModelAPI_Feature>(aAttr->owner());
 884         if (aConstrFeature->getKind() == SketchPlugin_ConstraintCoincidence::ID()) {
 885           if (!theCoincidencies.contains(aConstrFeature)) {
 886             std::shared_ptr<GeomAPI_Pnt2d> aPnt = getCoincedencePoint(aConstrFeature);
 887             if (aPnt.get() && aOrig->isEqual(aPnt)) {
 888               findCoincidences(aConstrFeature, theList, theCoincidencies,
 889                 SketchPlugin_ConstraintCoincidence::ENTITY_A());
 890               findCoincidences(aConstrFeature, theList, theCoincidencies,
 891                 SketchPlugin_ConstraintCoincidence::ENTITY_B());
 892             }
 893           }
 894         }
 895       }
 896     }
 897   }
 898 }
 899
 900 std::shared_ptr<GeomAPI_Pnt2d> PartSet_Tools::getCoincedencePoint(FeaturePtr theStartCoin)
 901 {
 902   std::shared_ptr<GeomAPI_Pnt2d> aPnt = SketcherPrs_Tools::getPoint(theStartCoin.get(),
 903                                                                     SketchPlugin_Constraint::ENTITY_A());
 904   if (aPnt.get() == NULL)
 905     aPnt = SketcherPrs_Tools::getPoint(theStartCoin.get(), SketchPlugin_Constraint::ENTITY_B());
 906   return aPnt;
 907 }
 908
 909 AttributePtr PartSet_Tools::findAttributeBy2dPoint(ObjectPtr theObj,
 910                                                    const TopoDS_Shape theShape,
 911                                                    FeaturePtr theSketch)
 912 {
 913
 914   AttributePtr anAttribute;
 915   FeaturePtr aFeature = ModelAPI_Feature::feature(theObj);
 916   if (aFeature) {
 917     if (theShape.ShapeType() == TopAbs_VERTEX) {
 918       const TopoDS_Vertex& aVertex = TopoDS::Vertex(theShape);
 919       if (!aVertex.IsNull())  {
 920         gp_Pnt aPoint = BRep_Tool::Pnt(aVertex);
 921         std::shared_ptr<GeomAPI_Pnt> aValue = std::shared_ptr<GeomAPI_Pnt>(
 922             new GeomAPI_Pnt(aPoint.X(), aPoint.Y(), aPoint.Z()));
 923
 924         // find the given point in the feature attributes
 925         std::list<AttributePtr> anAttiributes =
 926           aFeature->data()->attributes(GeomDataAPI_Point2D::typeId());
 927         std::list<AttributePtr>::const_iterator anIt = anAttiributes.begin(),
 928                                                 aLast = anAttiributes.end();
 929         for (; anIt != aLast && !anAttribute; anIt++) {
 930           std::shared_ptr<GeomDataAPI_Point2D> aCurPoint =
 931             std::dynamic_pointer_cast<GeomDataAPI_Point2D>(*anIt);
 932
 933           std::shared_ptr<GeomAPI_Pnt> aPnt = convertTo3D(aCurPoint->x(), aCurPoint->y(), theSketch);
 934           if (aPnt && (aPnt->distance(aValue) < Precision::Confusion())) {
 935             anAttribute = aCurPoint;
 936             break;
 937           }
 938         }
 939       }
 940     }
 941   }
 942   return anAttribute;
 943 }
 944
 945 void PartSet_Tools::sendSubFeaturesEvent(const CompositeFeaturePtr& theComposite,
 946                                          const Events_ID theEventId)
 947 {
 948   if (!theComposite.get())
 949     return;
 950
 951   static Events_Loop* aLoop = Events_Loop::loop();
 952   for (int i = 0; i < theComposite->numberOfSubs(); i++) {
 953     FeaturePtr aSubFeature = theComposite->subFeature(i);
 954     static const ModelAPI_EventCreator* aECreator = ModelAPI_EventCreator::get();
 955     aECreator->sendUpdated(aSubFeature, theEventId);
 956   }
 957   Events_Loop::loop()->flush(theEventId);
 958 }
 959
 960 bool PartSet_Tools::isAuxiliarySketchEntity(const ObjectPtr& theObject)
 961 {
 962   bool isAuxiliaryFeature = false;
 963
 964   FeaturePtr anObjectFeature = ModelAPI_Feature::feature(theObject);
 965   std::string anAuxiliaryAttribute = SketchPlugin_SketchEntity::AUXILIARY_ID();
 966   AttributeBooleanPtr anAuxiliaryAttr = std::dynamic_pointer_cast<ModelAPI_AttributeBoolean>(
 967                                     anObjectFeature->data()->attribute(anAuxiliaryAttribute));
 968   if (anAuxiliaryAttr.get())
 969     isAuxiliaryFeature = anAuxiliaryAttr->value();
 970
 971
 972   return isAuxiliaryFeature;
 973 }