TGeoArb8 and check overlap

From: Sue Kasahara <schubert_at_hep.umn.edu>
Date: Mon, 30 May 2005 19:41:21 -0500

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 - 02:42:05 MEST