G4TwistedTrd Class Reference

G4TwistedTrd is a twisted trapezoid with the X and Y dimensions varying along Z. More...

#include <G4TwistedTrd.hh>

Inheritance diagram for G4TwistedTrd:

Public Member Functions
	G4TwistedTrd (const G4String &pName, G4double pDx1, G4double pDx2, G4double pDy1, G4double pDy2, G4double pDz, G4double pPhiTwist)
	~G4TwistedTrd () override=default
G4double	GetX1HalfLength () const
G4double	GetX2HalfLength () const
G4double	GetY1HalfLength () const
G4double	GetY2HalfLength () const
G4double	GetZHalfLength () const
G4double	GetPhiTwist () const
G4GeometryType	GetEntityType () const override
G4double	GetCubicVolume () override
G4double	GetSurfaceArea () override
G4VSolid *	Clone () const override
std::ostream &	StreamInfo (std::ostream &os) const override
	G4TwistedTrd (__void__ &)
	G4TwistedTrd (const G4TwistedTrd &rhs)
G4TwistedTrd &	operator= (const G4TwistedTrd &rhs)
Public Member Functions inherited from G4VTwistedFaceted
	G4VTwistedFaceted (const G4String &pname, G4double PhiTwist, G4double pDz, G4double pTheta, G4double pPhi, G4double pDy1, G4double pDx1, G4double pDx2, G4double pDy2, G4double pDx3, G4double pDx4, G4double pAlph)
	~G4VTwistedFaceted () override
void	ComputeDimensions (G4VPVParameterisation *, const G4int, const G4VPhysicalVolume *) override
void	BoundingLimits (G4ThreeVector &pMin, G4ThreeVector &pMax) const override
G4bool	CalculateExtent (const EAxis pAxis, const G4VoxelLimits &pVoxelLimit, const G4AffineTransform &pTransform, G4double &pMin, G4double &pMax) const override
G4double	DistanceToIn (const G4ThreeVector &p, const G4ThreeVector &v) const override
G4double	DistanceToIn (const G4ThreeVector &p) const override
G4double	DistanceToOut (const G4ThreeVector &p, const G4ThreeVector &v, const G4bool calcnorm=false, G4bool *validnorm=nullptr, G4ThreeVector *n=nullptr) const override
G4double	DistanceToOut (const G4ThreeVector &p) const override
EInside	Inside (const G4ThreeVector &p) const override
G4ThreeVector	SurfaceNormal (const G4ThreeVector &p) const override
G4ThreeVector	GetPointOnSurface () const override
G4double	GetCubicVolume () override
G4double	GetSurfaceArea () override
void	DescribeYourselfTo (G4VGraphicsScene &scene) const override
G4Polyhedron *	CreatePolyhedron () const override
G4Polyhedron *	GetPolyhedron () const override
G4VisExtent	GetExtent () const override
std::ostream &	StreamInfo (std::ostream &os) const override
G4double	GetTwistAngle () const
G4double	GetDx1 () const
G4double	GetDx2 () const
G4double	GetDx3 () const
G4double	GetDx4 () const
G4double	GetDy1 () const
G4double	GetDy2 () const
G4double	GetDz () const
G4double	GetPhi () const
G4double	GetTheta () const
G4double	GetAlpha () const
G4double	Xcoef (G4double u, G4double phi, G4double ftg) const
G4double	GetValueA (G4double phi) const
G4double	GetValueB (G4double phi) const
G4double	GetValueD (G4double phi) const
G4GeometryType	GetEntityType () const override
	G4VTwistedFaceted (__void__ &)
	G4VTwistedFaceted (const G4VTwistedFaceted &rhs)
G4VTwistedFaceted &	operator= (const G4VTwistedFaceted &rhs)
Public Member Functions inherited from G4VSolid
	G4VSolid (const G4String &name)
virtual	~G4VSolid ()
	G4VSolid (const G4VSolid &rhs)
G4VSolid &	operator= (const G4VSolid &rhs)
G4bool	operator== (const G4VSolid &s) const
G4String	GetName () const
void	SetName (const G4String &name)
G4double	GetTolerance () const
virtual G4int	GetNumOfConstituents () const
virtual G4bool	IsFaceted () const
void	DumpInfo () const
virtual const G4VSolid *	GetConstituentSolid (G4int no) const
virtual G4VSolid *	GetConstituentSolid (G4int no)
virtual const G4DisplacedSolid *	GetDisplacedSolidPtr () const
virtual G4DisplacedSolid *	GetDisplacedSolidPtr ()
	G4VSolid (__void__ &)
G4double	EstimateCubicVolume (G4int nStat, G4double epsilon) const
G4double	EstimateSurfaceArea (G4int nStat, G4double epsilon) const

Additional Inherited Members
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
Protected Attributes inherited from G4VTwistedFaceted
G4bool	fRebuildPolyhedron = false
G4Polyhedron *	fpPolyhedron = nullptr
G4double	fCubicVolume = 0.0
G4double	fSurfaceArea = 0.0
Protected Attributes inherited from G4VSolid
G4double	kCarTolerance

Detailed Description

G4TwistedTrd is a twisted trapezoid with the X and Y dimensions varying along Z.

Definition at line 55 of file G4TwistedTrd.hh.

Constructor & Destructor Documentation

G4TwistedTrd() [1/3]

G4TwistedTrd::G4TwistedTrd	(	const G4String &	pName,
		G4double	pDx1,
		G4double	pDx2,
		G4double	pDy1,
		G4double	pDy2,
		G4double	pDz,
		G4double	pPhiTwist )

Constructs a twisted Trd, given its parameters.

Parameters

[in]	pName	The solid name.
[in]	pDx1	Half-length along X at the surface positioned at -dz.
[in]	pDx2	Half-length along X at the surface positioned at +dz.
[in]	pDy1	Half-length along Y at the surface positioned at -dz.
[in]	pDy2	Half-length along Y at the surface positioned at +dz.
[in]	pDz	Half-length along Z axis.
[in]	pPhiTwist	The twist angle.

Definition at line 45 of file G4TwistedTrd.cc.

  : G4VTwistedFaceted( pName, pPhiTwist,pDz,0.,0.,
                       pDy1, pDx1, pDx1, pDy2, pDx2, pDx2,0.)
{
}

Referenced by Clone(), G4TwistedTrd(), and operator=().

~G4TwistedTrd()

G4TwistedTrd::~G4TwistedTrd ( )

overridedefault

Default destructor.

G4TwistedTrd() [2/3]

G4TwistedTrd::G4TwistedTrd

(

__void__ &

a

)

Fake default constructor for usage restricted to direct object persistency for clients requiring preallocation of memory for persistifiable objects.

Definition at line 61 of file G4TwistedTrd.cc.

  : G4VTwistedFaceted(a)
{
}

G4TwistedTrd() [3/3]

G4TwistedTrd::G4TwistedTrd

(

const G4TwistedTrd &

rhs

)

Copy constructor and assignment operator.

Definition at line 69 of file G4TwistedTrd.cc.

  : G4VTwistedFaceted(rhs)
{
  fpPolyhedron = GetPolyhedron();
}

Member Function Documentation

Clone()

G4VSolid * G4TwistedTrd::Clone ( ) const

overridevirtual

Makes a clone of the object for use in multi-treading.

Returns

A pointer to the new cloned allocated solid.

Reimplemented from G4VSolid.

Definition at line 127 of file G4TwistedTrd.cc.

{
  return new G4TwistedTrd(*this);
}

GetCubicVolume()

double G4TwistedTrd::GetCubicVolume ( )

overridevirtual

Returning an estimation of the solid volume (capacity) and surface area, in internal units.

Reimplemented from G4VSolid.

Definition at line 135 of file G4TwistedTrd.cc.

{
  if (fCubicVolume == 0)
  {
    G4AutoLock l(&twtrdMutex);
    G4double x1 = GetX1HalfLength();
    G4double x2 = GetX2HalfLength();
    G4double y1 = GetY1HalfLength();
    G4double y2 = GetY2HalfLength();
    G4double h = 2.*GetZHalfLength();
    fCubicVolume = h*((x1 + x2)*(y1 + y2) + (x2 - x1)*(y2 - y1)/3.);
    l.unlock();
  }
  return fCubicVolume;
}

GetEntityType()

G4GeometryType G4TwistedTrd::GetEntityType ( ) const

overridevirtual

Returns the type ID, "G4TwistedTrd" of the solid.

Implements G4VSolid.

Definition at line 119 of file G4TwistedTrd.cc.

{
  return {"G4TwistedTrd"};
}

GetPhiTwist()

G4double G4TwistedTrd::GetPhiTwist ( ) const

inline

Definition at line 90 of file G4TwistedTrd.hh.

{ return GetTwistAngle() ; }

Referenced by G4tgbGeometryDumper::GetSolidParams(), GetSurfaceArea(), StreamInfo(), and G4GDMLWriteSolids::TwistedtrdWrite().

GetSurfaceArea()

double G4TwistedTrd::GetSurfaceArea ( )

overridevirtual

Returns an estimation of the solid surface area in internal units. This method may be overloaded by derived classes to compute the exact geometrical quantity for solids where this is possible, or anyway to cache the computed value. Note: the computed value is NOT cached.

Reimplemented from G4VSolid.

Definition at line 154 of file G4TwistedTrd.cc.

{
  if (fSurfaceArea == 0)
  {
    G4AutoLock l(&twtrdMutex);
    G4double ang = GetPhiTwist();
    G4double x1 = GetX1HalfLength();
    G4double x2 = GetX2HalfLength();
    G4double y1 = GetY1HalfLength();
    G4double y2 = GetY2HalfLength();
    G4double h = 2.*GetZHalfLength();
    G4double hh = h*h;
    G4double delX = x2 - x1;
    G4double delY = y2 - y1;
    if (ang == 0.)
    {
      G4double hx = std::sqrt(delY*delY + hh);
      G4double hy = std::sqrt(delX*delX + hh);
      return fSurfaceArea =
        2.*(x1 + x2)*hx + 2.*(y1 + y2)*hy + 4.*(x1*y1 + x2*y2);
    }
 
    // compute area of x-faces
    G4double U1, U2, V1, V2;
    G4double areaX = 0.;
    U1 = delY + x1*ang;
    U2 = delY + x2*ang;
    V1 = delY - x1*ang;
    V2 = delY - x2*ang;
    if (std::abs(delX) < kCarTolerance) // case x1 == x2
    {
      areaX = (U1*std::sqrt(hh + U1*U1) + hh*std::asinh(U1/h) -
               V1*std::sqrt(hh + V1*V1) - hh*std::asinh(V1/h))/ang;
    }
    else
    {
      // U contribution
      areaX += ((hh + U2*U2)*std::sqrt(hh + U2*U2) -
                (hh + U1*U1)*std::sqrt(hh + U1*U1))/3.
        + hh*(U2*std::asinh(U2/h) - U1*std::asinh(U1/h))
        - hh*(std::sqrt(hh + U2*U2) - std::sqrt(hh + U1*U1));
      // V contribution
      areaX += ((hh + V2*V2)*std::sqrt(hh + V2*V2) -
                (hh + V1*V1)*std::sqrt(hh + V1*V1))/3.
        + hh*(V2*std::asinh(V2/h) - V1*std::asinh(V1/h))
        - hh*(std::sqrt(hh + V2*V2) - std::sqrt(hh + V1*V1));
      areaX /= delX*ang*ang;
    }
 
    // compute area of y-faces
    G4double areaY = 0.;
    U1 = delX + y1*ang;
    U2 = delX + y2*ang;
    V1 = delX - y1*ang;
    V2 = delX - y2*ang;
    if (std::abs(delY) < kCarTolerance) // case y1 == y2
    {
      areaY = (U1*std::sqrt(hh + U1*U1) + hh*std::asinh(U1/h) -
               V1*std::sqrt(hh + V1*V1) - hh*std::asinh(V1/h))/ang;
    }
    else
    {
      // U contribution
      areaY += ((hh + U2*U2)*std::sqrt(hh + U2*U2) -
                (hh + U1*U1)*std::sqrt(hh + U1*U1))/3.
        + hh*(U2*std::asinh(U2/h) - U1*std::asinh(U1/h))
        - hh*(std::sqrt(hh + U2*U2) - std::sqrt(hh + U1*U1));
      // V contribution
      areaY += ((hh + V2*V2)*std::sqrt(hh + V2*V2) -
                (hh + V1*V1)*std::sqrt(hh + V1*V1))/3.
        + hh*(V2*std::asinh(V2/h) - V1*std::asinh(V1/h))
        - hh*(std::sqrt(hh + V2*V2) - std::sqrt(hh + V1*V1));
      areaY /= delY*ang*ang;
    }
    fSurfaceArea = areaX + areaY + 4.*(x1*y1 + x2*y2);
    l.unlock();
  }
  return fSurfaceArea;
}

GetX1HalfLength()

G4double G4TwistedTrd::GetX1HalfLength ( ) const

inline

Accessors.

Definition at line 85 of file G4TwistedTrd.hh.

{ return GetDx1() ; }

Referenced by GetCubicVolume(), G4tgbGeometryDumper::GetSolidParams(), GetSurfaceArea(), StreamInfo(), and G4GDMLWriteSolids::TwistedtrdWrite().

GetX2HalfLength()

G4double G4TwistedTrd::GetX2HalfLength ( ) const

inline

Definition at line 86 of file G4TwistedTrd.hh.

{ return GetDx3() ; }

Referenced by GetCubicVolume(), G4tgbGeometryDumper::GetSolidParams(), GetSurfaceArea(), StreamInfo(), and G4GDMLWriteSolids::TwistedtrdWrite().

GetY1HalfLength()

G4double G4TwistedTrd::GetY1HalfLength ( ) const

inline

Definition at line 87 of file G4TwistedTrd.hh.

{ return GetDy1() ; }

Referenced by GetCubicVolume(), G4tgbGeometryDumper::GetSolidParams(), GetSurfaceArea(), StreamInfo(), and G4GDMLWriteSolids::TwistedtrdWrite().

GetY2HalfLength()

G4double G4TwistedTrd::GetY2HalfLength ( ) const

inline

Definition at line 88 of file G4TwistedTrd.hh.

{ return GetDy2() ; }

Referenced by GetCubicVolume(), G4tgbGeometryDumper::GetSolidParams(), GetSurfaceArea(), StreamInfo(), and G4GDMLWriteSolids::TwistedtrdWrite().

GetZHalfLength()

G4double G4TwistedTrd::GetZHalfLength ( ) const

inline

Definition at line 89 of file G4TwistedTrd.hh.

{ return GetDz()  ; }

Referenced by GetCubicVolume(), G4tgbGeometryDumper::GetSolidParams(), GetSurfaceArea(), StreamInfo(), and G4GDMLWriteSolids::TwistedtrdWrite().

operator=()

G4TwistedTrd & G4TwistedTrd::operator=

(

const G4TwistedTrd &

rhs

)

Definition at line 78 of file G4TwistedTrd.cc.

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

StreamInfo()

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

overridevirtual

Streams the object contents to an output stream.

Implements G4VSolid.

Definition at line 95 of file G4TwistedTrd.cc.

{
  //
  // Stream object contents to an output stream
  //
  os << "-----------------------------------------------------------\n"
     << "    *** Dump for solid - " << GetName() << " ***\n"
     << "    ===================================================\n"
     << " Solid type: G4TwistedTrd\n"
     << " Parameters: \n"
     << "    pDx1 = " << GetX1HalfLength()/cm << " cm" << G4endl
     << "    pDx2 = " << GetX2HalfLength()/cm << " cm" << G4endl
     << "    pDy1 = " << GetY1HalfLength()/cm << " cm" << G4endl
     << "    pDy2 = " << GetY2HalfLength()/cm << " cm" << G4endl
     << "    pDz = "  << GetZHalfLength()/cm << " cm" << G4endl
     << "    pPhiTwist = " << GetPhiTwist()/degree << " deg" << G4endl
     << "-----------------------------------------------------------\n";
 
  return os;
}

---

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

• G4TwistedTrd.hh
• G4TwistedTrd.cc