Changeset 1315 for trunk/source/geometry/solids/Boolean/src/G4BooleanSolid.cc

Timestamp:

Jun 18, 2010, 11:42:07 AM (14 years ago)

Author:

garnier

Message:

update geant4-09-04-beta-cand-01 interfaces-V09-03-09 vis-V09-03-08

File:

• trunk/source/geometry/solids/Boolean/src/G4BooleanSolid.cc (modified) (3 diffs)

trunk/source/geometry/solids/Boolean/src/G4BooleanSolid.cc

@@ -25 +25 @@
 //
 //
-// $Id: G4BooleanSolid.cc,v 1.21 2006/10/19 15:34:49 gcosmo Exp $
-// GEANT4 tag $Name: geant4-09-03 $
+// $Id: G4BooleanSolid.cc,v 1.23 2010/05/11 15:03:45 gcosmo Exp $
+// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
 //
 // Implementation for the abstract base class for solids created by boolean 
@@ -40 +40 @@
 #include "G4VSolid.hh"
 #include "G4Polyhedron.hh"
+#include "HepPolyhedronProcessor.h"
 #include "Randomize.hh"
 
@@ -231 +232 @@
   return fpPolyhedron;
 }
+
+//////////////////////////////////////////////////////////////////////////
+//
+// Stacks polyhedra for processing. Returns top polyhedron.
+
+G4Polyhedron*
+G4BooleanSolid::StackPolyhedron(HepPolyhedronProcessor& processor,
+                                const G4VSolid* solid) const
+{
+  HepPolyhedronProcessor::Operation operation;
+  const G4String& type = solid->GetEntityType();
+  if (type == "G4UnionSolid")
+    { operation = HepPolyhedronProcessor::UNION; }
+  else if (type == "G4IntersectionSolid")
+    { operation = HepPolyhedronProcessor::INTERSECTION; }
+  else if (type == "G4SubtractionSolid")
+    { operation = HepPolyhedronProcessor::SUBTRACTION; }
+  else
+  {
+    std::ostringstream oss;
+    oss << "Solid - " << solid->GetName()
+        << " - Unrecognised composite solid"
+        << "\n  Returning NULL !";
+    G4Exception("StackPolyhedron()", "InvalidSetup",
+                JustWarning, oss.str().c_str());
+    return 0;
+  }
+
+  G4Polyhedron* top = 0;
+  const G4VSolid* solidA = solid->GetConstituentSolid(0);
+  const G4VSolid* solidB = solid->GetConstituentSolid(1);
+
+  if (solidA->GetConstituentSolid(0))
+  {
+    top = StackPolyhedron(processor, solidA);
+  }
+  else
+  {
+    top = solidA->GetPolyhedron();
+  }
+  G4Polyhedron* operand = solidB->GetPolyhedron();
+  processor.push_back (operation, *operand);
+
+  return top;
+}