HepPolyhedronEllipsoid Class Reference

#include <HepPolyhedron.h>

Inheritance diagram for HepPolyhedronEllipsoid:

Public Member Functions
	HepPolyhedronEllipsoid (G4double dx, G4double dy, G4double dz, G4double zcut1, G4double zcut2)
virtual	~HepPolyhedronEllipsoid ()
Public Member Functions inherited from HepPolyhedron
	HepPolyhedron ()
	HepPolyhedron (const HepPolyhedron &from)
virtual	~HepPolyhedron ()
HepPolyhedron &	operator= (const HepPolyhedron &from)
G4int	GetNoVertices () const
G4int	GetNoFacets () const
HepPolyhedron &	Transform (const G4Transform3D &t)
G4bool	GetNextVertexIndex (G4int &index, G4int &edgeFlag) const
G4Point3D	GetVertex (G4int index) const
G4bool	GetNextVertex (G4Point3D &vertex, G4int &edgeFlag) const
G4bool	GetNextVertex (G4Point3D &vertex, G4int &edgeFlag, G4Normal3D &normal) const
G4bool	GetNextEdgeIndeces (G4int &i1, G4int &i2, G4int &edgeFlag, G4int &iface1, G4int &iface2) const
G4bool	GetNextEdgeIndeces (G4int &i1, G4int &i2, G4int &edgeFlag) const
G4bool	GetNextEdge (G4Point3D &p1, G4Point3D &p2, G4int &edgeFlag) const
G4bool	GetNextEdge (G4Point3D &p1, G4Point3D &p2, G4int &edgeFlag, G4int &iface1, G4int &iface2) const
void	GetFacet (G4int iFace, G4int &n, G4int *iNodes, G4int *edgeFlags=0, G4int *iFaces=0) const
void	GetFacet (G4int iFace, G4int &n, G4Point3D *nodes, G4int *edgeFlags=0, G4Normal3D *normals=0) const
G4bool	GetNextFacet (G4int &n, G4Point3D *nodes, G4int *edgeFlags=0, G4Normal3D *normals=0) const
G4Normal3D	GetNormal (G4int iFace) const
G4Normal3D	GetUnitNormal (G4int iFace) const
G4bool	GetNextNormal (G4Normal3D &normal) const
G4bool	GetNextUnitNormal (G4Normal3D &normal) const
HepPolyhedron	add (const HepPolyhedron &p) const
HepPolyhedron	subtract (const HepPolyhedron &p) const
HepPolyhedron	intersect (const HepPolyhedron &p) const
G4double	GetSurfaceArea () const
G4double	GetVolume () const
G4int	createTwistedTrap (G4double Dz, const G4double xy1[][2], const G4double xy2[][2])
G4int	createPolyhedron (G4int Nnodes, G4int Nfaces, const G4double xyz[][3], const G4int faces[][4])

Additional Inherited Members
Static Public Member Functions inherited from HepPolyhedron
static G4int	GetNumberOfRotationSteps ()
static void	SetNumberOfRotationSteps (G4int n)
static void	ResetNumberOfRotationSteps ()
Protected Member Functions inherited from HepPolyhedron
void	AllocateMemory (G4int Nvert, G4int Nface)
G4int	FindNeighbour (G4int iFace, G4int iNode, G4int iOrder) const
G4Normal3D	FindNodeNormal (G4int iFace, G4int iNode) const
void	CreatePrism ()
void	RotateEdge (G4int k1, G4int k2, G4double r1, G4double r2, G4int v1, G4int v2, G4int vEdge, G4bool ifWholeCircle, G4int ns, G4int &kface)
void	SetSideFacets (G4int ii[4], G4int vv[4], G4int *kk, G4double *r, G4double dphi, G4int ns, G4int &kface)
void	RotateAroundZ (G4int nstep, G4double phi, G4double dphi, G4int np1, G4int np2, const G4double *z, G4double *r, G4int nodeVis, G4int edgeVis)
void	SetReferences ()
void	InvertFacets ()
Protected Attributes inherited from HepPolyhedron
G4int	nvert
G4int	nface
G4Point3D *	pV
G4Facet *	pF
Static Protected Attributes inherited from HepPolyhedron
static G4ThreadLocal G4int	fNumberOfRotationSteps = DEFAULT_NUMBER_OF_STEPS

Detailed Description

Definition at line 497 of file HepPolyhedron.h.

Constructor & Destructor Documentation

HepPolyhedronEllipsoid::HepPolyhedronEllipsoid	(	G4double	dx,
		G4double	dy,
		G4double	dz,
		G4double	zcut1,
		G4double	zcut2
	)

Definition at line 2058 of file HepPolyhedron.cc.

References test::a, FatalException, G4Exception(), HepPolyhedron::GetNumberOfRotationSteps(), HepPolyhedron::nvert, python.hepunit::pi, HepPolyhedron::pV, HepPolyhedron::RotateAroundZ(), HepPolyhedron::SetReferences(), HepGeom::BasicVector3D< T >::setX(), HepGeom::BasicVector3D< T >::setY(), python.hepunit::twopi, HepGeom::BasicVector3D< T >::x(), and HepGeom::BasicVector3D< T >::y().

{
  //   C H E C K   I N P U T   P A R A M E T E R S

  if (zCut1 >= cz || zCut2 <= -cz || zCut1 > zCut2) {
    std::cerr << "HepPolyhedronEllipsoid: wrong zCut1 = " << zCut1
           << " zCut2 = " << zCut2
           << " for given cz = " << cz << std::endl;
    return;
  }
  if (cz <= 0.0) {
    std::cerr << "HepPolyhedronEllipsoid: bad z semi-axis: cz = " << cz
      << std::endl;
    return;
  }

  G4double dthe;
  G4double sthe;
  G4int cutflag;
  cutflag= 0;
  if (zCut2 >= cz)
    {
      sthe= 0.0;
    }
  else
    {
      sthe= std::acos(zCut2/cz);
      cutflag++;
    }
  if (zCut1 <= -cz)
    {
      dthe= pi - sthe;
    }
  else
    {
      dthe= std::acos(zCut1/cz)-sthe;
      cutflag++;
    }

  //   P R E P A R E   T W O   P O L Y L I N E S
  //   generate sphere of radius cz first, then rescale x and y later

  G4int nds = (GetNumberOfRotationSteps() + 1) / 2;
  G4int np1 = G4int(dthe*nds/pi) + 2 + cutflag;

  G4double *zz, *rr;
  zz = new G4double[np1+1];
  rr = new G4double[np1+1];
  if (!zz || !rr)
    {
      G4Exception("HepPolyhedronEllipsoid::HepPolyhedronEllipsoid",
          "greps1002", FatalException, "Out of memory");
    }

  G4double a = dthe/(np1-cutflag-1);
  G4double cosa, sina;
  G4int j=0;
  if (sthe > 0.0)
    {
      zz[j]= zCut2;
      rr[j]= 0.;
      j++;
    }
  for (G4int i=0; i<np1-cutflag; i++) {
    cosa  = std::cos(sthe+i*a);
    sina  = std::sin(sthe+i*a);
    zz[j] = cz*cosa;
    rr[j] = cz*sina;
    j++;
  }
  if (j < np1)
    {
      zz[j]= zCut1;
      rr[j]= 0.;
      j++;
    }
  if (j > np1)
    {
      std::cerr << "Logic error in HepPolyhedronEllipsoid, memory corrupted!"
                << std::endl;
    }
  if (j < np1)
    {
      std::cerr << "Warning: logic error in HepPolyhedronEllipsoid."
                << std::endl;
      np1= j;
    }
  zz[j] = 0.;
  rr[j] = 0.;

  
  //   R O T A T E    P O L Y L I N E S

  RotateAroundZ(0, 0.0, twopi, np1, 1, zz, rr, -1, 1); 
  SetReferences();

  delete [] zz;
  delete [] rr;

  // rescale x and y vertex coordinates
  {
    G4Point3D * p= pV;
    for (G4int i=0; i<nvert; i++, p++) {
      p->setX( p->x() * ax/cz );
      p->setY( p->y() * by/cz );
    }
  }
}

HepPolyhedronEllipsoid::~HepPolyhedronEllipsoid ( )

virtual

Definition at line 2183 of file HepPolyhedron.cc.

{}

---

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

• HepPolyhedron.h
• HepPolyhedron.cc