Message: Re: Particle Problems at G4TessellatedSolid Boundary Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Particle Problems at G4TessellatedSolid Boundary 

Forum: Geometry
Re: Question Particle Problems at G4TessellatedSolid Boundary (Alex Leatherland)
Re: None Re: Particle Problems at G4TessellatedSolid Boundary (Ivana Hrivnacova)
Re: None Re: Particle Problems at G4TessellatedSolid Boundary (Alex Leatherland)
Date: 20 Feb, 2013
From: Alex Leatherland <Alex Leatherland>

Just noticed that the end of my message was missing:

> The reason I used this is because the implementation I have available
> for writing GDML

files based on CAD designs, outputs objects defined as tesselated shapes with triangular facets. So I wanted to confirm that the G4TesselatedSolid behaved correctly first.


I have also tried the same test of G4TessallatedSolid by hardcoding the volume description into the DetectorConstruction class, and the same incorrect result emerges.





G4TesselatedSolid* solidTarget = new G4TessellatedSolid("Target");
G4double side = 5*mm;

G4ThreeVector v1(side,side,side);
G4ThreeVector v2(side,side,-side);
G4ThreeVector v3(side,-side,side);
G4ThreeVector v4(side,-side,-side);
G4ThreeVector v5(-side,side,side);
G4ThreeVector v6(-side,side,-side);
G4ThreeVector v7(-side,-side,side);
G4ThreeVector v8(-side,-side,-side);

G4TriangularFacet *f1 = new G4TriangularFacet(v1,v3,v7,ABSOLUTE)
G4TriangularFacet *f2 = new G4TriangularFacet(v1,v5,v7,ABSOLUTE)
G4TriangularFacet *f3 = new G4TriangularFacet(v6,v8,v4,ABSOLUTE)
G4TriangularFacet *f4 = new G4TriangularFacet(v6,v4,v2,ABSOLUTE)
G4TriangularFacet *f5 = new G4TriangularFacet(v1,v2,v4,ABSOLUTE)
G4TriangularFacet *f6 = new G4TriangularFacet(v1,v4,v3,ABSOLUTE)
G4TriangularFacet *f7 = new G4TriangularFacet(v7,v8,v6,ABSOLUTE)
G4TriangularFacet *f8 = new G4TriangularFacet(v7,v6,v5,ABSOLUTE)
G4TriangularFacet *f9 = new G4TriangularFacet(v1,v5,v6,ABSOLUTE)
G4TriangularFacet *f10 = new G4TriangularFacet(v1,v6,v2,ABSOLUTE)
G4TriangularFacet *f11 = new G4TriangularFacet(v4,v8,v7,ABSOLUTE)
G4TriangularFacet *f12 = new G4TriangularFacet(v4,v7,v3,ABSOLUTE)

solidTarget->AddFacet((G4VFacet*) f1);
solidTarget->AddFacet((G4VFacet*) f2);
solidTarget->AddFacet((G4VFacet*) f3);
solidTarget->AddFacet((G4VFacet*) f4);
solidTarget->AddFacet((G4VFacet*) f5);
solidTarget->AddFacet((G4VFacet*) f6);
solidTarget->AddFacet((G4VFacet*) f7);
solidTarget->AddFacet((G4VFacet*) f8);
solidTarget->AddFacet((G4VFacet*) f9);
solidTarget->AddFacet((G4VFacet*) f10);
solidTarget->AddFacet((G4VFacet*) f11);
solidTarget->AddFacet((G4VFacet*) f12);

solidTarget->SetSolidClosed(true)

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

1 None: Re: Particle Problems at G4TessellatedSolid Boundary   (Ivana Hrivnacova - 20 Feb, 2013)
(_ None: Re: Particle Problems at G4TessellatedSolid Boundary   (Alex Leatherland - 21 Feb, 2013)
(_ News: Re: Particle Problems at G4TessellatedSolid Boundary   (Alex Leatherland - 22 Feb, 2013)
(_ None: Re: Particle Problems at G4TessellatedSolid Boundary   (Ivana Hrivnacova - 27 Feb, 2013)
(_ None: Re: Particle Problems at G4TessellatedSolid Boundary   (Michael H. Kelsey - 27 Feb, 2013)
(_ Ok: Re: Particle Problems at G4TessellatedSolid Boundary   (Alex Leatherland - 07 Mar, 2013)
 Add Message Add Message
to: "Re: Particle Problems at G4TessellatedSolid Boundary"

 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 ]