Message: Re: Polyhedron::SetReferences: List N is not empty Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

More Re: Polyhedron::SetReferences: List N is not empty 

Forum: Geometry
Re: Question Polyhedron::SetReferences: List N is not empty (Nicola Mori)
Re: None Re: Polyhedron::SetReferences: List N is not empty (John Allison)
Date: 05 Apr, 2015
From: Nicola Mori <Nicola Mori>

I did what you suggested, but either I don't understand the output of the dump or there is something really screwed in my geometry. To better understand the dump I produced a simple geometry which does not give the Polyhedron::SetReferences message but still I cannot understand the dump. The geometry consist of a cubic box made of two triangular facets on the upper side (at Z=1000):

Y

|   3--2
|   | /|
|   |/ |
|   1--4
|
-------- X

and 5 quadrangular facets on the other 5 sides. To obey the rule of "anticlockwise when seen from the outside" I build the two G4TriangularFacets using these ordering of points: 1-2-3 and 1-4-2. The dump shows me:

Nverteces=8, Nfacets=7

xyz(1)=-500 -500 1000

xyz(2)=500 500 1000

xyz(3)=-500 500 1000

xyz(4)=500 -500 1000

xyz(5)=-500 -500 0

xyz(6)=500 -500 0

xyz(7)=500 500 0

xyz(8)=-500 500 0

face(1)= 1/2 2/4 3/5 0/0

face(2)= 1/3 4/6 2/1 0/0

face(3)= 5/7 6/6 4/2 1/5

face(4)= 7/7 8/5 3/1 2/6

face(5)= 8/7 5/3 1/1 3/4

face(6)= 6/7 7/4 2/2 4/3

face(7)= 5/5 8/4 7/6 6/3

The first 4 points are the 1,2,3 and 4 from the picture above bounding the upper face of the cube, while the others are in the Z=0 layer. My understanding of the dump of the facets is that facet 1 has an edge going from point 1 to point 2 (correct), another edge going from point 2 to point 4 (correct start but wrong end) and another edge going from point 3 to point 5 (???). This does not make sense to me so I tend to believe that I'm misinterpreting the dump. Still I can't figure out which is the correct interpretation. By the way the geometry is visualized perfectly; I also tried to define a facet clockwise to see if I misinterpreted the anticlockwise rule, but defining it clockwise it is not drawn correctly, so I think that I can handle the anticlockwise rule correctly...

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: Polyhedron::SetReferences: List N is not empty   (John Allison - 09 Apr, 2015)
(_ Ok: Re: Polyhedron::SetReferences: List N is not empty   (Nicola Mori - 09 Apr, 2015)
 Add Message Add Message
to: "Re: Polyhedron::SetReferences: List N is not empty"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]