#include <G4ParameterisationPolyhedra.hh>

Inheritance diagram for G4VParameterisationPolyhedra:

Public Member Functions
virtual void	ComputeDimensions (G4Box &, const G4int, const G4VPhysicalVolume *) const

virtual void	ComputeDimensions (G4Cons &, const G4int, const G4VPhysicalVolume *) const

virtual void	ComputeDimensions (G4Ellipsoid &, const G4int, const G4VPhysicalVolume *) const

virtual void	ComputeDimensions (G4Hype &, const G4int, const G4VPhysicalVolume *) const

virtual void	ComputeDimensions (G4Orb &, const G4int, const G4VPhysicalVolume *) const

virtual void	ComputeDimensions (G4Para &, const G4int, const G4VPhysicalVolume *) const

virtual void	ComputeDimensions (G4Polycone &, const G4int, const G4VPhysicalVolume *) const

virtual void	ComputeDimensions (G4Polyhedra &, const G4int, const G4VPhysicalVolume *) const

virtual void	ComputeDimensions (G4Sphere &, const G4int, const G4VPhysicalVolume *) const

virtual void	ComputeDimensions (G4Torus &, const G4int, const G4VPhysicalVolume *) const

virtual void	ComputeDimensions (G4Trap &, const G4int, const G4VPhysicalVolume *) const

virtual void	ComputeDimensions (G4Trd &, const G4int, const G4VPhysicalVolume *) const

virtual void	ComputeDimensions (G4Tubs &, const G4int, const G4VPhysicalVolume *) const

virtual G4Material *	ComputeMaterial (const G4int repNo, G4VPhysicalVolume currentVol, const G4VTouchable parentTouch=nullptr)

virtual G4VSolid *	ComputeSolid (const G4int, G4VPhysicalVolume *)

virtual void	ComputeTransformation (const G4int copyNo, G4VPhysicalVolume *physVol) const =0

	G4VParameterisationPolyhedra (EAxis axis, G4int nCopies, G4double offset, G4double step, G4VSolid *msolid, DivisionType divType)

EAxis	GetAxis () const

G4double	GetHalfGap () const

virtual G4VVolumeMaterialScanner *	GetMaterialScanner ()

G4VSolid *	GetMotherSolid () const

G4int	GetNoDiv () const

G4double	GetOffset () const

const G4String &	GetType () const

G4double	GetWidth () const

virtual G4bool	IsNested () const

void	SetHalfGap (G4double hg)

void	SetType (const G4String &type)

G4int	VolumeFirstCopyNo () const

virtual	~G4VParameterisationPolyhedra ()

Protected Member Functions
G4int	CalculateNDiv (G4double motherDim, G4double width, G4double offset) const

G4double	CalculateWidth (G4double motherDim, G4int nDiv, G4double offset) const

void	ChangeRotMatrix (G4VPhysicalVolume *physVol, G4double rotZ=0.0) const

void	CheckNDivAndWidth (G4double maxPar)

void	CheckOffset (G4double maxPar)

virtual void	CheckParametersValidity ()

virtual G4double	GetMaxParameter () const =0

G4double	OffsetZ () const

Protected Attributes
EAxis	faxis

G4bool	fDeleteSolid = false

DivisionType	fDivisionType

G4double	fhgap = 0.0

G4VSolid *	fmotherSolid = nullptr

G4int	fnDiv = 0

G4double	foffset = 0.0

G4bool	fReflectedSolid = false

G4String	ftype

G4double	fwidth = 0.0

G4double	kCarTolerance

G4int	theVoluFirstCopyNo = 1

Static Protected Attributes
static G4ThreadLocal G4RotationMatrix *	fRot = nullptr

static const G4int	verbose = 5

Private Member Functions
G4double	ConvertRadiusFactor (const G4Polyhedra &phedra) const

Detailed Description

Definition at line 58 of file G4ParameterisationPolyhedra.hh.

Constructor & Destructor Documentation

◆ G4VParameterisationPolyhedra()

G4VParameterisationPolyhedra::G4VParameterisationPolyhedra	(	EAxis	axis,
		G4int	nCopies,
		G4double	offset,
		G4double	step,
		G4VSolid *	msolid,
		DivisionType	divType
	)

Definition at line 45 of file G4ParameterisationPolyhedra.cc.

  :  G4VDivisionParameterisation( axis, nDiv, width, offset, divType, msolid )
{
  std::ostringstream message;
  /* #ifdef G4MULTITHREADED
  message << "Divisions for G4Polyhedra currently NOT supported in MT-mode."
          << G4endl
          << "Sorry! Solid: " << msolid->GetName();
  G4Exception("G4VParameterisationPolyhedra::G4VParameterisationPolyhedra()",
              "GeomDiv0001", FatalException, message);
          #endif */
 
  G4Polyhedra* msol = (G4Polyhedra*)(msolid);
  if ((msolid->GetEntityType() != "G4ReflectedSolid") && (msol->IsGeneric()))
  {
    message << "Generic construct for G4Polyhedra NOT supported." << G4endl
            << "Sorry! Solid: " << msol->GetName();
    G4Exception("G4VParameterisationPolyhedra::G4VParameterisationPolyhedra()",
                "GeomDiv0001", FatalException, message);
  }
  if (msolid->GetEntityType() == "G4ReflectedSolid")
  {
    // Get constituent solid
    //  
    G4VSolid* mConstituentSolid 
       = ((G4ReflectedSolid*)msolid)->GetConstituentMovedSolid();
    msol = (G4Polyhedra*)(mConstituentSolid);
  
    // Get parameters
    //
    G4int nofSides = msol->GetOriginalParameters()->numSide;
    G4int nofZplanes = msol->GetOriginalParameters()->Num_z_planes;
    G4double* zValues = msol->GetOriginalParameters()->Z_values;
    G4double* rminValues = msol->GetOriginalParameters()->Rmin;
    G4double* rmaxValues = msol->GetOriginalParameters()->Rmax;
 
    // Invert z values, convert radius parameters
    //
    G4double* rminValues2 = new G4double[nofZplanes];
    G4double* rmaxValues2 = new G4double[nofZplanes];
    G4double* zValuesRefl = new G4double[nofZplanes];
    for (G4int i=0; i<nofZplanes; ++i)
    {
      rminValues2[i] = rminValues[i] * ConvertRadiusFactor(*msol);
      rmaxValues2[i] = rmaxValues[i] * ConvertRadiusFactor(*msol);
      zValuesRefl[i] = - zValues[i];
    }  
    
    G4Polyhedra* newSolid
      = new G4Polyhedra(msol->GetName(),
                        msol->GetStartPhi(), 
                        msol->GetEndPhi() - msol->GetStartPhi(),
                        nofSides,
                        nofZplanes, zValuesRefl, rminValues2, rmaxValues2);
 
    delete [] rminValues2;       
    delete [] rmaxValues2;       
    delete [] zValuesRefl;       
 
    msol = newSolid;
    fmotherSolid = newSolid;
    fReflectedSolid = true;
    fDeleteSolid = true;
  }    
}

References ConvertRadiusFactor(), FatalException, G4VDivisionParameterisation::fDeleteSolid, G4VDivisionParameterisation::fmotherSolid, G4VDivisionParameterisation::fReflectedSolid, G4endl, G4Exception(), G4Polyhedra::GetEndPhi(), G4VSolid::GetEntityType(), G4VSolid::GetName(), G4Polyhedra::GetOriginalParameters(), G4Polyhedra::GetStartPhi(), G4Polyhedra::IsGeneric(), G4PolyhedraHistorical::Num_z_planes, G4PolyhedraHistorical::numSide, G4PolyhedraHistorical::Rmax, G4PolyhedraHistorical::Rmin, and G4PolyhedraHistorical::Z_values.

◆ ~G4VParameterisationPolyhedra()

G4VParameterisationPolyhedra::~G4VParameterisationPolyhedra ( )

virtual

Definition at line 115 of file G4ParameterisationPolyhedra.cc.

116{

117}

Member Function Documentation

◆ CalculateNDiv()

G4int G4VDivisionParameterisation::CalculateNDiv	(	G4double	motherDim,
		G4double	width,
		G4double	offset
	)		const

protectedinherited

Definition at line 99 of file G4VDivisionParameterisation.cc.

{
#ifdef G4DIVDEBUG
  G4cout << " G4VDivisionParameterisation::CalculateNDiv: "
         << ( motherDim - offset ) / width 
         << " Motherdim: " <<  motherDim << ", Offset: " << offset
         << ", Width: " << width << G4endl;
#endif
 
  return G4int( ( motherDim - offset ) / width );
}

References G4cout, and G4endl.

◆ CalculateWidth()

G4double G4VDivisionParameterisation::CalculateWidth	(	G4double	motherDim,
		G4int	nDiv,
		G4double	offset
	)		const

protectedinherited

Definition at line 114 of file G4VDivisionParameterisation.cc.

{ 
#ifdef G4DIVDEBUG
  G4cout << " G4VDivisionParameterisation::CalculateWidth: "
         << ( motherDim - offset ) / nDiv
         << ", Motherdim: " << motherDim << ", Offset: " << offset
         << ", Number of divisions: " << nDiv << G4endl;
#endif
 
  return ( motherDim - offset ) / nDiv;
}

References G4cout, and G4endl.

◆ ChangeRotMatrix()

void G4VDivisionParameterisation::ChangeRotMatrix	(	G4VPhysicalVolume *	physVol,
		G4double	rotZ = `0.0`
	)		const

protectedinherited

Definition at line 85 of file G4VDivisionParameterisation.cc.

{
  if (fRot == nullptr)
  {
    fRot = new G4RotationMatrix();
    G4AutoDelete::Register(fRot);
  }
  fRot->rotateZ( rotZ );
  physVol->SetRotation(fRot);
}

References G4VDivisionParameterisation::fRot, G4AutoDelete::Register(), CLHEP::HepRotation::rotateZ(), and G4VPhysicalVolume::SetRotation().

◆ CheckNDivAndWidth()

void G4VDivisionParameterisation::CheckNDivAndWidth ( G4double maxPar )

protectedinherited

Definition at line 151 of file G4VDivisionParameterisation.cc.

{
  if( (fDivisionType == DivNDIVandWIDTH)
      && (foffset + fwidth*fnDiv - maxPar > kCarTolerance ) )
  {
    std::ostringstream message;
    message << "Configuration not supported." << G4endl
            << "Division of solid " << fmotherSolid->GetName()
           << " has too big offset + width*nDiv = " << G4endl
           << "        " << foffset + fwidth*fnDiv << " > "
           << foffset << ". Width = "
           << G4endl
           << "        " << fwidth << ". nDiv = " << fnDiv << " !";
    G4Exception("G4VDivisionParameterisation::CheckNDivAndWidth()",
                "GeomDiv0001", FatalException, message);
  }
}

References DivNDIVandWIDTH, FatalException, G4VDivisionParameterisation::fDivisionType, G4VDivisionParameterisation::fmotherSolid, G4VDivisionParameterisation::fnDiv, G4VDivisionParameterisation::foffset, G4VDivisionParameterisation::fwidth, G4endl, G4Exception(), G4VSolid::GetName(), and G4VDivisionParameterisation::kCarTolerance.

Referenced by G4VDivisionParameterisation::CheckParametersValidity().

◆ CheckOffset()

void G4VDivisionParameterisation::CheckOffset ( G4double maxPar )

protectedinherited

Definition at line 136 of file G4VDivisionParameterisation.cc.

{
  if( foffset >= maxPar )
  {
    std::ostringstream message;
    message << "Configuration not supported." << G4endl
            << "Division of solid " << fmotherSolid->GetName()
            << " has too big offset = " << G4endl
            << "        " << foffset << " > " << maxPar << " !";
    G4Exception("G4VDivisionParameterisation::CheckOffset()",
                "GeomDiv0001", FatalException, message);
  }
}

References FatalException, G4VDivisionParameterisation::fmotherSolid, G4VDivisionParameterisation::foffset, G4endl, G4Exception(), and G4VSolid::GetName().

Referenced by G4VDivisionParameterisation::CheckParametersValidity().

◆ CheckParametersValidity()

void G4VDivisionParameterisation::CheckParametersValidity ( )

protectedvirtualinherited

Reimplemented in G4ParameterisationPolyconeRho, G4ParameterisationPolyconeZ, G4ParameterisationPolyhedraRho, G4ParameterisationPolyhedraPhi, and G4ParameterisationPolyhedraZ.

Definition at line 128 of file G4VDivisionParameterisation.cc.

{
  G4double maxPar = GetMaxParameter();
  CheckOffset( maxPar );
  CheckNDivAndWidth( maxPar );
}

References G4VDivisionParameterisation::CheckNDivAndWidth(), G4VDivisionParameterisation::CheckOffset(), and G4VDivisionParameterisation::GetMaxParameter().

◆ ComputeDimensions() [1/13]

virtual void G4VPVParameterisation::ComputeDimensions	(	G4Box &	,
		const	G4int,
		const G4VPhysicalVolume *
	)		const

inlinevirtualinherited

Reimplemented in G4EzVoxelParameterization, G4ParameterisationConsRho, G4ParameterisationConsPhi, G4ParameterisationConsZ, G4ParameterisationParaX, G4ParameterisationParaY, G4ParameterisationParaZ, G4ParameterisationPolyconeRho, G4ParameterisationPolyconePhi, G4ParameterisationPolyconeZ, G4ParameterisationPolyhedraRho, G4ParameterisationPolyhedraPhi, G4ParameterisationPolyhedraZ, G4ParameterisationTrdX, G4ParameterisationTrdY, G4ParameterisationTrdZ, G4ParameterisationTubsRho, G4ParameterisationTubsPhi, G4ParameterisationTubsZ, G4VNestedParameterisation, G4PhantomParameterisation, G4GDMLParameterisation, G4EzVoxelParameterization, G4ParameterisationBoxX, G4ParameterisationBoxY, and G4ParameterisationBoxZ.

Definition at line 89 of file G4VPVParameterisation.hh.

91 {}

Referenced by G4Box::ComputeDimensions(), G4Cons::ComputeDimensions(), G4Orb::ComputeDimensions(), G4Para::ComputeDimensions(), G4Sphere::ComputeDimensions(), G4Torus::ComputeDimensions(), G4Trap::ComputeDimensions(), G4Trd::ComputeDimensions(), G4Tubs::ComputeDimensions(), G4Ellipsoid::ComputeDimensions(), G4Hype::ComputeDimensions(), G4Polycone::ComputeDimensions(), G4Polyhedra::ComputeDimensions(), G4tgbGeometryDumper::DumpPVParameterised(), and G4GDMLWriteParamvol::ParametersWrite().

◆ ComputeDimensions() [2/13]

virtual void G4VPVParameterisation::ComputeDimensions	(	G4Cons &	,
		const	G4int,
		const G4VPhysicalVolume *
	)		const

inlinevirtualinherited

◆ ComputeDimensions() [3/13]

virtual void G4VPVParameterisation::ComputeDimensions	(	G4Ellipsoid &	,
		const	G4int,
		const G4VPhysicalVolume *
	)		const

inlinevirtualinherited

◆ ComputeDimensions() [4/13]

virtual void G4VPVParameterisation::ComputeDimensions	(	G4Hype &	,
		const	G4int,
		const G4VPhysicalVolume *
	)		const

inlinevirtualinherited

◆ ComputeDimensions() [5/13]

virtual void G4VPVParameterisation::ComputeDimensions	(	G4Orb &	,
		const	G4int,
		const G4VPhysicalVolume *
	)		const

inlinevirtualinherited

◆ ComputeDimensions() [6/13]

virtual void G4VPVParameterisation::ComputeDimensions	(	G4Para &	,
		const	G4int,
		const G4VPhysicalVolume *
	)		const

inlinevirtualinherited

◆ ComputeDimensions() [7/13]

virtual void G4VPVParameterisation::ComputeDimensions	(	G4Polycone &	,
		const	G4int,
		const G4VPhysicalVolume *
	)		const

inlinevirtualinherited

◆ ComputeDimensions() [8/13]

virtual void G4VPVParameterisation::ComputeDimensions	(	G4Polyhedra &	,
		const	G4int,
		const G4VPhysicalVolume *
	)		const

inlinevirtualinherited

◆ ComputeDimensions() [9/13]

virtual void G4VPVParameterisation::ComputeDimensions	(	G4Sphere &	,
		const	G4int,
		const G4VPhysicalVolume *
	)		const

inlinevirtualinherited

◆ ComputeDimensions() [10/13]

virtual void G4VPVParameterisation::ComputeDimensions	(	G4Torus &	,
		const	G4int,
		const G4VPhysicalVolume *
	)		const

inlinevirtualinherited

◆ ComputeDimensions() [11/13]

virtual void G4VPVParameterisation::ComputeDimensions	(	G4Trap &	,
		const	G4int,
		const G4VPhysicalVolume *
	)		const

inlinevirtualinherited

◆ ComputeDimensions() [12/13]

virtual void G4VPVParameterisation::ComputeDimensions	(	G4Trd &	,
		const	G4int,
		const G4VPhysicalVolume *
	)		const

inlinevirtualinherited

◆ ComputeDimensions() [13/13]

virtual void G4VPVParameterisation::ComputeDimensions	(	G4Tubs &	,
		const	G4int,
		const G4VPhysicalVolume *
	)		const

inlinevirtualinherited

◆ ComputeMaterial()

G4Material * G4VPVParameterisation::ComputeMaterial	(	const G4int	repNo,
		G4VPhysicalVolume *	currentVol,
		const G4VTouchable *	parentTouch = `nullptr`
	)

virtualinherited

Reimplemented in G4VNestedParameterisation, G4PartialPhantomParameterisation, and G4PhantomParameterisation.

Definition at line 47 of file G4VPVParameterisation.cc.

{
  return pPhysicalVol->GetLogicalVolume()->GetMaterial();
}

References G4VPhysicalVolume::GetLogicalVolume(), and G4LogicalVolume::GetMaterial().

Referenced by G4tgbGeometryDumper::DumpPVParameterised(), G4LogicalVolume::GetMass(), G4ParameterisedNavigation::LevelLocate(), and G4Region::ScanVolumeTree().

◆ ComputeSolid()

G4VSolid * G4VDivisionParameterisation::ComputeSolid	(	const G4int	i,
		G4VPhysicalVolume *	pv
	)

virtualinherited

Reimplemented from G4VPVParameterisation.

Definition at line 72 of file G4VDivisionParameterisation.cc.

{
  G4VSolid* solid = G4VPVParameterisation::ComputeSolid(i, pv);
  if (solid->GetEntityType() == "G4ReflectedSolid")
  {
    solid = ((G4ReflectedSolid*)solid)->GetConstituentMovedSolid();
  }
  return solid;
}      

References G4VPVParameterisation::ComputeSolid(), and G4VSolid::GetEntityType().

◆ ComputeTransformation()

virtual void G4VDivisionParameterisation::ComputeTransformation	(	const G4int	copyNo,
		G4VPhysicalVolume *	physVol
	)		const

pure virtualinherited

◆ ConvertRadiusFactor()

G4double G4VParameterisationPolyhedra::ConvertRadiusFactor ( const G4Polyhedra & phedra ) const

private

Definition at line 121 of file G4ParameterisationPolyhedra.cc.

{
  G4double phiTotal = phedra.GetEndPhi() - phedra.GetStartPhi();
  G4int nofSides = phedra.GetOriginalParameters()->numSide;
  
  if ( (phiTotal <=0) || (phiTotal >
        2*pi+G4GeometryTolerance::GetInstance()->GetAngularTolerance()) )
    { phiTotal = 2*pi; }
  
  return std::cos(0.5*phiTotal/nofSides);
}  

References G4Polyhedra::GetEndPhi(), G4GeometryTolerance::GetInstance(), G4Polyhedra::GetOriginalParameters(), G4Polyhedra::GetStartPhi(), G4PolyhedraHistorical::numSide, and pi.

Referenced by G4VParameterisationPolyhedra().

◆ GetAxis()

EAxis G4VDivisionParameterisation::GetAxis ( ) const

inlineinherited

Referenced by G4PVDivisionFactory::CreatePVDivision().

◆ GetHalfGap()

G4double G4VDivisionParameterisation::GetHalfGap ( ) const

inlineinherited

◆ GetMaterialScanner()

G4VVolumeMaterialScanner * G4VPVParameterisation::GetMaterialScanner ( )

virtualinherited

Reimplemented in G4VNestedParameterisation.

Definition at line 62 of file G4VPVParameterisation.cc.

{
  return nullptr;
}

Referenced by G4Region::ScanVolumeTree().

◆ GetMaxParameter()

virtual G4double G4VDivisionParameterisation::GetMaxParameter ( ) const

protectedpure virtualinherited

◆ GetMotherSolid()

G4VSolid * G4VDivisionParameterisation::GetMotherSolid ( ) const

inlineinherited

Referenced by G4ParameterisationConsZ::ComputeTransformation().

◆ GetNoDiv()

G4int G4VDivisionParameterisation::GetNoDiv ( ) const

inlineinherited

Referenced by G4ReplicatedSlice::CheckAndSetParameters(), G4PVDivision::CheckAndSetParameters(), and G4PVDivisionFactory::CreatePVDivision().

◆ GetOffset()

G4double G4VDivisionParameterisation::GetOffset ( ) const

inlineinherited

Referenced by G4PVDivisionFactory::CreatePVDivision().

◆ GetType()

const G4String & G4VDivisionParameterisation::GetType ( ) const

inlineinherited

◆ GetWidth()

G4double G4VDivisionParameterisation::GetWidth ( ) const

inlineinherited

Referenced by G4ReplicatedSlice::CheckAndSetParameters(), G4PVDivision::CheckAndSetParameters(), and G4PVDivisionFactory::CreatePVDivision().

◆ IsNested()

G4bool G4VPVParameterisation::IsNested ( ) const

virtualinherited

Reimplemented in G4VNestedParameterisation.

Definition at line 55 of file G4VPVParameterisation.cc.

{
  return false;
}

Referenced by G4GMocrenFileSceneHandler::AddSolid(), G4Navigator::SetupHierarchy(), and G4ITNavigator2::SetupHierarchy().

◆ OffsetZ()

G4double G4VDivisionParameterisation::OffsetZ ( ) const

protectedinherited

Definition at line 170 of file G4VDivisionParameterisation.cc.

{
  // take into account reflection in the offset
 
  G4double offset = foffset;
  if (fReflectedSolid)
  {
    offset = GetMaxParameter() - fwidth*fnDiv - foffset;
  }
 
  return offset;
}  