#include <G4UnionSolid.hh>

Inheritance diagram for G4UnionSolid:

Public Member Functions
	G4UnionSolid (const G4String &pName, G4VSolid pSolidA, G4VSolid pSolidB)

	G4UnionSolid (const G4String &pName, G4VSolid pSolidA, G4VSolid pSolidB, G4RotationMatrix *rotMatrix, const G4ThreeVector &transVector)

	G4UnionSolid (const G4String &pName, G4VSolid pSolidA, G4VSolid pSolidB, const G4Transform3D &transform)

virtual	~G4UnionSolid ()

G4GeometryType	GetEntityType () const

G4VSolid *	Clone () const

	G4UnionSolid (__void__ &)

	G4UnionSolid (const G4UnionSolid &rhs)

G4UnionSolid &	operator= (const G4UnionSolid &rhs)

G4bool	CalculateExtent (const EAxis pAxis, const G4VoxelLimits &pVoxelLimit, const G4AffineTransform &pTransform, G4double &pMin, G4double &pMax) const

EInside	Inside (const G4ThreeVector &p) const

G4ThreeVector	SurfaceNormal (const G4ThreeVector &p) const

G4double	DistanceToIn (const G4ThreeVector &p, const G4ThreeVector &v) const

G4double	DistanceToIn (const G4ThreeVector &p) const

G4double	DistanceToOut (const G4ThreeVector &p, const G4ThreeVector &v, const G4bool calcNorm=false, G4bool validNorm=0, G4ThreeVector n=0) const

G4double	DistanceToOut (const G4ThreeVector &p) const

void	ComputeDimensions (G4VPVParameterisation p, const G4int n, const G4VPhysicalVolume pRep)

void	DescribeYourselfTo (G4VGraphicsScene &scene) const

G4Polyhedron *	CreatePolyhedron () const

Public Member Functions inherited from G4BooleanSolid
	G4BooleanSolid (const G4String &pName, G4VSolid pSolidA, G4VSolid pSolidB)

	G4BooleanSolid (const G4String &pName, G4VSolid pSolidA, G4VSolid pSolidB, G4RotationMatrix *rotMatrix, const G4ThreeVector &transVector)

	G4BooleanSolid (const G4String &pName, G4VSolid pSolidA, G4VSolid pSolidB, const G4Transform3D &transform)

virtual	~G4BooleanSolid ()

virtual const G4VSolid *	GetConstituentSolid (G4int no) const

virtual G4VSolid *	GetConstituentSolid (G4int no)

G4double	GetCubicVolume ()

G4double	GetSurfaceArea ()

virtual G4Polyhedron *	GetPolyhedron () const

std::ostream &	StreamInfo (std::ostream &os) const

G4int	GetCubVolStatistics () const

G4double	GetCubVolEpsilon () const

void	SetCubVolStatistics (G4int st)

void	SetCubVolEpsilon (G4double ep)

G4int	GetAreaStatistics () const

G4double	GetAreaAccuracy () const

void	SetAreaStatistics (G4int st)

void	SetAreaAccuracy (G4double ep)

G4ThreeVector	GetPointOnSurface () const

	G4BooleanSolid (__void__ &)

	G4BooleanSolid (const G4BooleanSolid &rhs)

G4BooleanSolid &	operator= (const G4BooleanSolid &rhs)

Public Member Functions inherited from G4VSolid
	G4VSolid (const G4String &name)

virtual	~G4VSolid ()

G4bool	operator== (const G4VSolid &s) const

G4String	GetName () const

void	SetName (const G4String &name)

G4double	GetTolerance () const

void	DumpInfo () const

virtual G4VisExtent	GetExtent () const

virtual const G4DisplacedSolid *	GetDisplacedSolidPtr () const

virtual G4DisplacedSolid *	GetDisplacedSolidPtr ()

	G4VSolid (__void__ &)

	G4VSolid (const G4VSolid &rhs)

G4VSolid &	operator= (const G4VSolid &rhs)

Additional Inherited Members
Protected Member Functions inherited from G4BooleanSolid
G4Polyhedron *	StackPolyhedron (HepPolyhedronProcessor &, const G4VSolid *) const

G4double	GetAreaRatio () const

Protected Member Functions inherited from G4VSolid
void	CalculateClippedPolygonExtent (G4ThreeVectorList &pPolygon, const G4VoxelLimits &pVoxelLimit, const EAxis pAxis, G4double &pMin, G4double &pMax) const

void	ClipCrossSection (G4ThreeVectorList *pVertices, const G4int pSectionIndex, const G4VoxelLimits &pVoxelLimit, const EAxis pAxis, G4double &pMin, G4double &pMax) const

void	ClipBetweenSections (G4ThreeVectorList *pVertices, const G4int pSectionIndex, const G4VoxelLimits &pVoxelLimit, const EAxis pAxis, G4double &pMin, G4double &pMax) const

void	ClipPolygon (G4ThreeVectorList &pPolygon, const G4VoxelLimits &pVoxelLimit, const EAxis pAxis) const

G4double	EstimateCubicVolume (G4int nStat, G4double epsilon) const

G4double	EstimateSurfaceArea (G4int nStat, G4double ell) const

Protected Attributes inherited from G4BooleanSolid
G4VSolid *	fPtrSolidA

G4VSolid *	fPtrSolidB

G4double	fAreaRatio

Protected Attributes inherited from G4VSolid
G4double	kCarTolerance

Detailed Description

Definition at line 52 of file G4UnionSolid.hh.

Constructor & Destructor Documentation

G4UnionSolid::G4UnionSolid	(	const G4String &	pName,
		G4VSolid *	pSolidA,
		G4VSolid *	pSolidB
	)

Definition at line 58 of file G4UnionSolid.cc.

Referenced by Clone().

   : G4BooleanSolid(pName,pSolidA,pSolidB)
 {
 }

G4UnionSolid::G4UnionSolid	(	const G4String &	pName,
		G4VSolid *	pSolidA,
		G4VSolid *	pSolidB,
		G4RotationMatrix *	rotMatrix,
		const G4ThreeVector &	transVector
	)

Definition at line 69 of file G4UnionSolid.cc.

   : G4BooleanSolid(pName,pSolidA,pSolidB,rotMatrix,transVector)
 
 {
 }

G4UnionSolid::G4UnionSolid	(	const G4String &	pName,
		G4VSolid *	pSolidA,
		G4VSolid *	pSolidB,
		const G4Transform3D &	transform
	)

Definition at line 83 of file G4UnionSolid.cc.

   : G4BooleanSolid(pName,pSolidA,pSolidB,transform)
 {
 } 

G4UnionSolid::~G4UnionSolid ( )

virtual

Definition at line 105 of file G4UnionSolid.cc.

106 {

107 }

G4UnionSolid::G4UnionSolid ( __void__ & a )

Definition at line 96 of file G4UnionSolid.cc.

   : G4BooleanSolid(a)
 {
 }

G4UnionSolid::G4UnionSolid ( const G4UnionSolid & rhs )

Definition at line 113 of file G4UnionSolid.cc.

   : G4BooleanSolid (rhs)
 {
 }

Member Function Documentation

G4bool G4UnionSolid::CalculateExtent	(	const EAxis	pAxis,
		const G4VoxelLimits &	pVoxelLimit,
		const G4AffineTransform &	pTransform,
		G4double &	pMin,
		G4double &	pMax
	)		const

virtual

Implements G4VSolid.

Definition at line 140 of file G4UnionSolid.cc.

References G4VSolid::CalculateExtent(), G4BooleanSolid::fPtrSolidA, G4BooleanSolid::fPtrSolidB, G4INCL::Math::max(), and G4INCL::Math::min().

 {
   G4bool   touchesA, touchesB, out ;
   G4double minA =  kInfinity, minB =  kInfinity, 
            maxA = -kInfinity, maxB = -kInfinity; 
 
   touchesA = fPtrSolidA->CalculateExtent( pAxis, pVoxelLimit, 
                                           pTransform, minA, maxA);
   touchesB= fPtrSolidB->CalculateExtent( pAxis, pVoxelLimit, 
                                          pTransform, minB, maxB);
   if( touchesA || touchesB )
   {
     pMin = std::min( minA, minB ); 
     pMax = std::max( maxA, maxB );
     out  = true ; 
   }
   else out = false ;
 
   return out ;  // It exists in this slice if either one does.
 }

G4VSolid * G4UnionSolid::Clone ( ) const

virtual

Reimplemented from G4VSolid.

Definition at line 457 of file G4UnionSolid.cc.

References G4UnionSolid().

 {
   return new G4UnionSolid(*this);
 }

void G4UnionSolid::ComputeDimensions	(	G4VPVParameterisation *	p,
		const G4int	n,
		const G4VPhysicalVolume *	pRep
	)

virtual

Reimplemented from G4VSolid.

Definition at line 467 of file G4UnionSolid.cc.

470 {

471 }

G4Polyhedron * G4UnionSolid::CreatePolyhedron ( ) const

virtual

Reimplemented from G4VSolid.

Definition at line 488 of file G4UnionSolid.cc.

References HepPolyhedronProcessor::execute(), processor, and G4BooleanSolid::StackPolyhedron().

 {
   HepPolyhedronProcessor processor;
   // Stack components and components of components recursively
   // See G4BooleanSolid::StackPolyhedron
   G4Polyhedron* top = StackPolyhedron(processor, this);
   G4Polyhedron* result = new G4Polyhedron(*top);
   if (processor.execute(*result)) { return result; }
   else { return 0; }
 }

void G4UnionSolid::DescribeYourselfTo ( G4VGraphicsScene & scene ) const

virtual

Implements G4VSolid.

Definition at line 478 of file G4UnionSolid.cc.

References G4VGraphicsScene::AddSolid().

 {
   scene.AddSolid (*this);
 }

G4double G4UnionSolid::DistanceToIn	(	const G4ThreeVector &	p,
		const G4ThreeVector &	v
	)		const

virtual

Implements G4VSolid.

Definition at line 252 of file G4UnionSolid.cc.

References G4VSolid::DistanceToIn(), G4BooleanSolid::fPtrSolidA, G4BooleanSolid::fPtrSolidB, G4cerr, G4cout, G4endl, Inside(), kInside, and G4INCL::Math::min().

 {
 #ifdef G4BOOLDEBUG
   if( Inside(p) == kInside )
   {
     G4cout << "WARNING - Invalid call in "
            << "G4UnionSolid::DistanceToIn(p,v)" << G4endl
            << "  Point p is inside !" << G4endl;
     G4cout << "          p = " << p << G4endl;
     G4cout << "          v = " << v << G4endl;
     G4cerr << "WARNING - Invalid call in "
            << "G4UnionSolid::DistanceToIn(p,v)" << G4endl
            << "  Point p is inside !" << G4endl;
     G4cerr << "          p = " << p << G4endl;
     G4cerr << "          v = " << v << G4endl;
   }
 #endif
 
   return std::min(fPtrSolidA->DistanceToIn(p,v),
                     fPtrSolidB->DistanceToIn(p,v) ) ;
 }

G4double G4UnionSolid::DistanceToIn ( const G4ThreeVector & p ) const

virtual

Implements G4VSolid.

Definition at line 281 of file G4UnionSolid.cc.

References G4VSolid::DistanceToIn(), G4BooleanSolid::fPtrSolidA, G4BooleanSolid::fPtrSolidB, G4cerr, G4cout, G4endl, Inside(), kInside, and G4INCL::Math::min().

 {
 #ifdef G4BOOLDEBUG
   if( Inside(p) == kInside )
   {
     G4cout << "WARNING - Invalid call in "
            << "G4UnionSolid::DistanceToIn(p)" << G4endl
            << "  Point p is inside !" << G4endl;
     G4cout << "          p = " << p << G4endl;
     G4cerr << "WARNING - Invalid call in "
            << "G4UnionSolid::DistanceToIn(p)" << G4endl
            << "  Point p is inside !" << G4endl;
     G4cerr << "          p = " << p << G4endl;
   }
 #endif
   G4double distA = fPtrSolidA->DistanceToIn(p) ;
   G4double distB = fPtrSolidB->DistanceToIn(p) ;
   G4double safety = std::min(distA,distB) ;
   if(safety < 0.0) safety = 0.0 ;
   return safety ;
 }

G4double G4UnionSolid::DistanceToOut	(	const G4ThreeVector &	p,
		const G4ThreeVector &	v,
		const G4bool	calcNorm = `false`,
		G4bool *	validNorm = `0`,
		G4ThreeVector *	n = `0`
	)		const

virtual

Implements G4VSolid.

Definition at line 308 of file G4UnionSolid.cc.

References G4VSolid::DistanceToOut(), G4BooleanSolid::fPtrSolidA, G4BooleanSolid::fPtrSolidB, G4cerr, G4cout, G4endl, Inside(), G4VSolid::Inside(), G4VSolid::kCarTolerance, kOutside, python.hepunit::mm, CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

 {
   G4double  dist = 0.0, disTmp = 0.0 ;
   G4ThreeVector normTmp;
   G4ThreeVector* nTmp= &normTmp;
 
   if( Inside(p) == kOutside )
   {
 #ifdef G4BOOLDEBUG
       G4cout << "Position:"  << G4endl << G4endl;
       G4cout << "p.x() = "   << p.x()/mm << " mm" << G4endl;
       G4cout << "p.y() = "   << p.y()/mm << " mm" << G4endl;
       G4cout << "p.z() = "   << p.z()/mm << " mm" << G4endl << G4endl;
       G4cout << "Direction:" << G4endl << G4endl;
       G4cout << "v.x() = "   << v.x() << G4endl;
       G4cout << "v.y() = "   << v.y() << G4endl;
       G4cout << "v.z() = "   << v.z() << G4endl << G4endl;
       G4cout << "WARNING - Invalid call in "
              << "G4UnionSolid::DistanceToOut(p,v)" << G4endl
              << "  Point p is outside !" << G4endl;
       G4cout << "          p = " << p << G4endl;
       G4cout << "          v = " << v << G4endl;
       G4cerr << "WARNING - Invalid call in "
              << "G4UnionSolid::DistanceToOut(p,v)" << G4endl
              << "  Point p is outside !" << G4endl;
       G4cerr << "          p = " << p << G4endl;
       G4cerr << "          v = " << v << G4endl;
 #endif
   }
   else
   {
     EInside positionA = fPtrSolidA->Inside(p) ;
     // EInside positionB = fPtrSolidB->Inside(p) ;
 
     if( positionA != kOutside )
     { 
       do
       {
         disTmp = fPtrSolidA->DistanceToOut(p+dist*v,v,calcNorm,
                                            validNorm,nTmp)        ;
         dist += disTmp ;
 
         if(fPtrSolidB->Inside(p+dist*v) != kOutside)
         { 
           disTmp = fPtrSolidB->DistanceToOut(p+dist*v,v,calcNorm,
                                             validNorm,nTmp)         ;
           dist += disTmp ;
         }
       }
       //     while( Inside(p+dist*v) == kInside ) ;
            while( fPtrSolidA->Inside(p+dist*v) != kOutside && 
                   disTmp > 0.5*kCarTolerance ) ;
     }
     else // if( positionB != kOutside )
     {
       do
       {
         disTmp = fPtrSolidB->DistanceToOut(p+dist*v,v,calcNorm,
                                            validNorm,nTmp)        ; 
         dist += disTmp ;
 
         if(fPtrSolidA->Inside(p+dist*v) != kOutside)
         { 
           disTmp = fPtrSolidA->DistanceToOut(p+dist*v,v,calcNorm,
                                             validNorm,nTmp)         ;
           dist += disTmp ;
         }
       }
       //  while( Inside(p+dist*v) == kInside ) ;
         while( (fPtrSolidB->Inside(p+dist*v) != kOutside)
             && (disTmp > 0.5*kCarTolerance) ) ;
     }
   }
   if( calcNorm )
   { 
      *validNorm = false ;
      *n         = *nTmp ;   
   }
   return dist ;
 }

G4double G4UnionSolid::DistanceToOut ( const G4ThreeVector & p ) const

virtual

Implements G4VSolid.

Definition at line 398 of file G4UnionSolid.cc.

References G4VSolid::DistanceToOut(), G4BooleanSolid::fPtrSolidA, G4BooleanSolid::fPtrSolidB, G4cerr, G4cout, G4endl, Inside(), G4VSolid::Inside(), kInside, kOutside, kSurface, and G4INCL::Math::max().

 {
   G4double distout = 0.0;
   if( Inside(p) == kOutside )
   {
 #ifdef G4BOOLDEBUG
     G4cout << "WARNING - Invalid call in "
            << "G4UnionSolid::DistanceToOut(p)" << G4endl
            << "  Point p is outside !" << G4endl;
     G4cout << "          p = " << p << G4endl;
     G4cerr << "WARNING - Invalid call in "
            << "G4UnionSolid::DistanceToOut(p)" << G4endl
            << "  Point p is outside !" << G4endl;
     G4cerr << "          p = " << p << G4endl;
 #endif
   }
   else
   {
     EInside positionA = fPtrSolidA->Inside(p) ;
     EInside positionB = fPtrSolidB->Inside(p) ;
   
     //  Is this equivalent ??
     //    if( ! (  (positionA == kOutside)) && 
     //             (positionB == kOutside))  ) 
     if((positionA == kInside  && positionB == kInside  ) ||
        (positionA == kInside  && positionB == kSurface ) ||
        (positionA == kSurface && positionB == kInside  )     )
     {     
       distout= std::max(fPtrSolidA->DistanceToOut(p),
                           fPtrSolidB->DistanceToOut(p) ) ;
     }
     else
     {
       if(positionA == kOutside)
       {
         distout= fPtrSolidB->DistanceToOut(p) ;
       }
       else
       {
         distout= fPtrSolidA->DistanceToOut(p) ;
       }
     }
   }
   return distout;
 }

G4GeometryType G4UnionSolid::GetEntityType ( ) const

virtual

Reimplemented from G4BooleanSolid.

Definition at line 448 of file G4UnionSolid.cc.

 {
   return G4String("G4UnionSolid");
 }

EInside G4UnionSolid::Inside ( const G4ThreeVector & p ) const

virtual

Implements G4VSolid.

Definition at line 171 of file G4UnionSolid.cc.

References G4BooleanSolid::fPtrSolidA, G4BooleanSolid::fPtrSolidB, G4GeometryTolerance::GetInstance(), G4GeometryTolerance::GetRadialTolerance(), G4VSolid::Inside(), kInside, kOutside, kSurface, and G4VSolid::SurfaceNormal().

Referenced by DistanceToIn(), DistanceToOut(), and SurfaceNormal().

 {
   EInside positionA = fPtrSolidA->Inside(p);
   if (positionA == kInside)  { return kInside; }
 
   EInside positionB = fPtrSolidB->Inside(p);
 
   if( positionB == kInside  ||
     ( positionA == kSurface && positionB == kSurface &&
         ( fPtrSolidA->SurfaceNormal(p) + 
           fPtrSolidB->SurfaceNormal(p) ).mag2() < 
           1000*G4GeometryTolerance::GetInstance()->GetRadialTolerance() ) )
   {
     return kInside;
   }
   else
   {
     if( ( positionB == kSurface ) || ( positionA == kSurface ) )
       { return kSurface; }
     else
       { return kOutside; } 
   }
 }

G4UnionSolid & G4UnionSolid::operator= ( const G4UnionSolid & rhs )

Definition at line 122 of file G4UnionSolid.cc.

References G4BooleanSolid::operator=().

 {
   // Check assignment to self
   //
   if (this == &rhs)  { return *this; }
 
   // Copy base class data
   //
   G4BooleanSolid::operator=(rhs);
 
   return *this;
 }  

G4ThreeVector G4UnionSolid::SurfaceNormal ( const G4ThreeVector & p ) const

virtual

Implements G4VSolid.

Definition at line 200 of file G4UnionSolid.cc.

References G4BooleanSolid::fPtrSolidA, G4BooleanSolid::fPtrSolidB, G4cerr, G4cout, G4endl, Inside(), G4VSolid::Inside(), kInside, kOutside, kSurface, and G4VSolid::SurfaceNormal().

 {
     G4ThreeVector normal;
 
 #ifdef G4BOOLDEBUG
     if( Inside(p) == kOutside )
     {
       G4cout << "WARNING - Invalid call in "
              << "G4UnionSolid::SurfaceNormal(p)" << G4endl
              << "  Point p is outside !" << G4endl;
       G4cout << "          p = " << p << G4endl;
       G4cerr << "WARNING - Invalid call in "
              << "G4UnionSolid::SurfaceNormal(p)" << G4endl
              << "  Point p is outside !" << G4endl;
       G4cerr << "          p = " << p << G4endl;
     }
 #endif
 
     if(fPtrSolidA->Inside(p) == kSurface && fPtrSolidB->Inside(p) != kInside) 
     {
        normal= fPtrSolidA->SurfaceNormal(p) ;
     }
     else if(fPtrSolidB->Inside(p) == kSurface && 
             fPtrSolidA->Inside(p) != kInside)
     {
        normal= fPtrSolidB->SurfaceNormal(p) ;
     }
     else 
     {
       normal= fPtrSolidA->SurfaceNormal(p) ;
 #ifdef G4BOOLDEBUG
       if(Inside(p)==kInside)
       {
         G4cout << "WARNING - Invalid call in "
              << "G4UnionSolid::SurfaceNormal(p)" << G4endl
              << "  Point p is inside !" << G4endl;
         G4cout << "          p = " << p << G4endl;
         G4cerr << "WARNING - Invalid call in "
              << "G4UnionSolid::SurfaceNormal(p)" << G4endl
              << "  Point p is inside !" << G4endl;
         G4cerr << "          p = " << p << G4endl;
       }
 #endif
     }
     return normal;
 }

The documentation for this class was generated from the following files:

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation