# -*- coding: iso-8859-1 -*-

###
### This file is generated automatically by SALOME v6.3.0 with dump python functionality
###

import sys
import salome

salome.salome_init()
theStudy = salome.myStudy

import GEOM
import geompy
import math
import SALOMEDS


geompy.init_geom(theStudy)

OZ = geompy.MakeVectorDXDYDZ(0,0,1)
Cylinder_1 = geompy.MakeCylinderRH(10, 300)
Cylinder_2 = geompy.MakeCylinderRH(15, 30)
Vector_1 = geompy.MakeVectorDXDYDZ(0, 1, 0)
geompy.TranslateDXDYDZ(Cylinder_2, 0, 0, -15)
geompy.Rotate(Cylinder_2, Vector_1, 90*math.pi/180.0)
Cut_1 = geompy.MakeCut(Cylinder_1, Cylinder_2)
tool = geompy.MakeBox( 0,-50,-50, 50,50,350 )
NewCyl = geompy.MakeCut( Cut_1, tool )
[ToolFace,Bottom,Top,Face_1] = geompy.ExtractShapes(NewCyl, geompy.ShapeType["FACE"], True)



geompy.addToStudy( Cut_1, 'Cut_1' )
geompy.addToStudy( NewCyl, 'NewCyl' )
geompy.addToStudyInFather( NewCyl, Bottom, 'Bottom' )
geompy.addToStudyInFather( NewCyl, Top, 'Top' )


###
### SMESH component
###

import smesh, SMESH, SALOMEDS, math

smesh.SetCurrentStudy(theStudy)
Mesh_1        = smesh.Mesh(NewCyl)
Regular_1D    = Mesh_1.Segment()
Nb_Segments_1 = Regular_1D.NumberOfSegments(15)
Quadrangle_2D = Mesh_1.Quadrangle()
Prism_3D      = Mesh_1.Prism()
MEFISTO_2D    = Mesh_1.Triangle(geom=Top)
Projection_2D = Mesh_1.Projection2D(geom=Bottom)
Source_Face_1 = Projection_2D.SourceFace(Top)
isDone        = Mesh_1.Compute()
Mesh_1.Rotate([],OZ,math.pi,Copy=True)

