Re: TGeoArb8 and check overlap

I forgot to mention that I'm using root v4.03/05 last updated from cvs on 4/25/2005.
-Sue

Sue Kasahara wrote:

> Hi root team,
> I'm attempting to insert a volume with a TGeoArb8 shape as
> a node inside a mother volume with a TGeoArb8 shape, and I find
> that CheckOverlaps run on the mother volume detects an extrusion
> where there should be none. I'm wondering if this is a bug or if I'm
> misunderstanding
> some part of the code usage. I'm attaching a simple example to
> illustrate the problem.
> Thanks in advance for your help!
> -Sue
>
> {
> // usage:
> // root[0] .x thisscript.C
> //
> gSystem->Load("libGeom.so");
>
> new TGeoManager("testGeo","test overlaps");
>
> TGeoArb8* outShape = new TGeoArb8("F0001_0_tVa",0.0055);
> for ( int iz = 0; iz < 2; iz++ ) {
> outShape -> SetVertex(0+iz*4,-1.7715,-0.57575);
> outShape -> SetVertex(1+iz*4,+1.7715,-0.57575);
> outShape -> SetVertex(2+iz*4,+2.906398,+0.57565);
> outShape -> SetVertex(3+iz*4,-2.906398,+0.57565);
> }
> TGeoVolume* outVolume = new TGeoVolume(outShape->GetName(),outShape);
>
>
> TGeoArb8* inShape = new TGeoArb8("F0001_0_tVa_Air",0.005);
> for ( int iz = 0; iz < 2; iz++ ) {
> inShape -> SetVertex(0+iz*4,-1.7710,-0.57525);
> inShape -> SetVertex(1+iz*4,+1.7710,-0.57525);
> inShape -> SetVertex(2+iz*4,+2.905898,+0.57515);
> inShape -> SetVertex(3+iz*4,-2.905898,+0.57515);
> }
> TGeoVolume* inVolume = new TGeoVolume(inShape->GetName(),inShape);
>
> outVolume->AddNode(inVolume,1,gGeoIdentity);
>
> gGeoManager->SetTopVolume(outVolume);
> gGeoManager->CloseGeometry();
>
> gGeoManager->CheckOverlaps();
> gGeoManager->PrintOverlaps();
>
> }
>
Received on Tue May 31 2005 - 03:45:57 MEST