#include <G4ParameterisationTrd.hh>

Inheritance diagram for G4ParameterisationTrdX:

Public Member Functions
void	ComputeDimensions (G4Trap &trd, const G4int copyNo, const G4VPhysicalVolume *pv) const

void	ComputeDimensions (G4Trd &trd, const G4int copyNo, const G4VPhysicalVolume *pv) const

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

virtual G4VSolid *	ComputeSolid (const G4int, G4VPhysicalVolume *)

void	ComputeTransformation (const G4int copyNo, G4VPhysicalVolume *physVol) const

	G4ParameterisationTrdX (EAxis axis, G4int nCopies, G4double width, G4double offset, G4VSolid *motherSolid, DivisionType divType)

EAxis	GetAxis () const

G4double	GetHalfGap () const

virtual G4VVolumeMaterialScanner *	GetMaterialScanner ()

G4double	GetMaxParameter () const

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

	~G4ParameterisationTrdX ()

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 ()

G4double	OffsetZ () const

Protected Attributes
G4bool	bDivInTrap = false

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
void	ComputeDimensions (G4Box &, const G4int, const G4VPhysicalVolume *) const

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

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

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

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

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

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

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

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

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

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

Detailed Description

Definition at line 75 of file G4ParameterisationTrd.hh.

Constructor & Destructor Documentation

◆ G4ParameterisationTrdX()

G4ParameterisationTrdX::G4ParameterisationTrdX	(	EAxis	axis,
		G4int	nCopies,
		G4double	width,
		G4double	offset,
		G4VSolid *	motherSolid,
		DivisionType	divType
	)

Definition at line 78 of file G4ParameterisationTrd.cc.

  :  G4VParameterisationTrd( axis, nDiv, width, offset, msolid, divType )
{
  CheckParametersValidity();
  SetType( "DivisionTrdX" );
 
  G4Trd* msol = (G4Trd*)(fmotherSolid);
  if( divType == DivWIDTH )
  {
    fnDiv = CalculateNDiv( msol->GetXHalfLength1()+msol->GetXHalfLength2(),
                           width, offset );
  }
  else if( divType == DivNDIV )
  {
    fwidth = CalculateWidth( msol->GetXHalfLength1()+msol->GetXHalfLength2(),
                             nDiv, offset );
  }
 
#ifdef G4DIVDEBUG
  if( verbose >= 1 )
  {
    G4cout << " G4ParameterisationTrdX - ## divisions " << fnDiv << " = "
           << nDiv << G4endl
           << " Offset " << foffset << " = " << offset << G4endl
           << " Width " << fwidth << " = " << width << G4endl;
  }
#endif
 
  G4double mpDx1 = msol->GetXHalfLength1();
  G4double mpDx2 = msol->GetXHalfLength2();
  if( std::fabs(mpDx1 - mpDx2) > kCarTolerance )
  {
    bDivInTrap = true;
  }
}

◆ ~G4ParameterisationTrdX()

G4ParameterisationTrdX::~G4ParameterisationTrdX ( )

Definition at line 118 of file G4ParameterisationTrd.cc.

119{

120}

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]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 99 of file G4ParameterisationTrd.hh.

100 {}

◆ ComputeDimensions() [2/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 97 of file G4ParameterisationTrd.hh.

98 {}

◆ ComputeDimensions() [3/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 105 of file G4ParameterisationTrd.hh.

106 {}

◆ ComputeDimensions() [4/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 111 of file G4ParameterisationTrd.hh.

112 {}

◆ ComputeDimensions() [5/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 103 of file G4ParameterisationTrd.hh.

104 {}

◆ ComputeDimensions() [6/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 109 of file G4ParameterisationTrd.hh.

110 {}

◆ ComputeDimensions() [7/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 115 of file G4ParameterisationTrd.hh.

116 {}

◆ ComputeDimensions() [8/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 117 of file G4ParameterisationTrd.hh.

118 {}

◆ ComputeDimensions() [9/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 101 of file G4ParameterisationTrd.hh.

102 {}

◆ ComputeDimensions() [10/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 107 of file G4ParameterisationTrd.hh.

108 {}

◆ ComputeDimensions() [11/13]

void G4ParameterisationTrdX::ComputeDimensions	(	G4Trap &	trd,
		const G4int	copyNo,
		const G4VPhysicalVolume *	pv
	)		const

virtual

Reimplemented from G4VPVParameterisation.

Definition at line 192 of file G4ParameterisationTrd.cc.

{
  G4Trd* msol = (G4Trd*)(fmotherSolid); 
  G4double pDy1 = msol->GetYHalfLength1();
  G4double pDy2 = msol->GetYHalfLength2();
  G4double pDz = msol->GetZHalfLength();
  //fwidth is at Z=0
  G4double pDx1 = msol->GetXHalfLength1(); 
  G4double pDx2 = msol->GetXHalfLength2();
  //  G4double pDxAVE = (pDx1+pDx2)/2.;
  G4double xChangeRatio = (pDx2-pDx1) / (pDx2+pDx1);
  G4double fWidChange = xChangeRatio*fwidth;
  G4double fWid1 = fwidth - fWidChange;
  G4double fWid2 = fwidth + fWidChange;
  G4double fOffsetChange = xChangeRatio*foffset/2.;
  G4double fOffset1 = foffset - fOffsetChange;
  G4double fOffset2 = foffset + fOffsetChange;
  G4double cxy1 = -pDx1+fOffset1 + (copyNo+0.5)*fWid1;// centre of the side at y=-pDy1
  G4double cxy2 = -pDx2+fOffset2 + (copyNo+0.5)*fWid2;// centre of the side at y=+pDy1
  G4double alp = std::atan( (cxy2-cxy1)/(pDz*2.) );
 
  pDx1 = fwidth/2. - fWidChange/2.;
  pDx2 = fwidth/2. + fWidChange/2.;
 
 
  trap.SetAllParameters ( pDz,
              alp,
              0.,
              pDy1,
              pDx1,
              pDx1,
              0.,
              pDy2,
              pDx2 - fhgap,
              pDx2 - fhgap * pDx2/pDx1,
              0.);
 
#ifdef G4DIVDEBUG
  if( verbose >= 2 )
  {
    G4cout << " G4ParameterisationTrdX::ComputeDimensions():"
           << copyNo << G4endl;
    trap.DumpInfo();
  }
#endif
}

References G4VSolid::DumpInfo(), G4VDivisionParameterisation::fhgap, G4VDivisionParameterisation::fmotherSolid, G4VDivisionParameterisation::foffset, G4VDivisionParameterisation::fwidth, G4cout, G4endl, G4Trd::GetXHalfLength1(), G4Trd::GetXHalfLength2(), G4Trd::GetYHalfLength1(), G4Trd::GetYHalfLength2(), G4Trd::GetZHalfLength(), G4Trap::SetAllParameters(), and G4VDivisionParameterisation::verbose.

◆ ComputeDimensions() [12/13]

void G4ParameterisationTrdX::ComputeDimensions	(	G4Trd &	trd,
		const G4int	copyNo,
		const G4VPhysicalVolume *	pv
	)		const

virtual

Reimplemented from G4VPVParameterisation.

Definition at line 169 of file G4ParameterisationTrd.cc.

{  
  G4Trd* msol = (G4Trd*)(fmotherSolid);
  G4double pDy1 = msol->GetYHalfLength1();
  G4double pDy2 = msol->GetYHalfLength2();
  G4double pDz = msol->GetZHalfLength();
  G4double pDx = fwidth/2. - fhgap;
 
  trd.SetAllParameters ( pDx, pDx, pDy1, pDy2, pDz );
 
#ifdef G4DIVDEBUG
  if( verbose >= 2 )
  {
    G4cout << " G4ParameterisationTrdX::ComputeDimensions():"
           << copyNo << G4endl;
    trd.DumpInfo();
  }
#endif
}

References G4VSolid::DumpInfo(), G4VDivisionParameterisation::fhgap, G4VDivisionParameterisation::fmotherSolid, G4VDivisionParameterisation::fwidth, G4cout, G4endl, G4Trd::GetYHalfLength1(), G4Trd::GetYHalfLength2(), G4Trd::GetZHalfLength(), G4Trd::SetAllParameters(), and G4VDivisionParameterisation::verbose.

◆ ComputeDimensions() [13/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 113 of file G4ParameterisationTrd.hh.

114 {}

◆ 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()

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

virtual

Implements G4VDivisionParameterisation.

Definition at line 131 of file G4ParameterisationTrd.cc.

{
  G4Trd* msol = (G4Trd*)(fmotherSolid );
  G4double mdx = ( msol->GetXHalfLength1() + msol->GetXHalfLength2() ) / 2.;
  //----- translation 
  G4ThreeVector origin(0.,0.,0.); 
  G4double posi;
  posi = -mdx + foffset + (copyNo+0.5)*fwidth;
  if( faxis == kXAxis )
  {
    origin.setX( posi ); 
  }
  else
  { 
    std::ostringstream message;
    message << "Only axes along X are allowed !  Axis: " << faxis;
    G4Exception("G4ParameterisationTrdX::ComputeTransformation()",
                "GeomDiv0002", FatalException, message);
  }
 
#ifdef G4DIVDEBUG
  if( verbose >= 2 )
  {
    G4cout << std::setprecision(8)
           << " G4ParameterisationTrdX::ComputeTransformation() "
           << copyNo << G4endl
           << " Position: " << origin << " - Axis: " << faxis << G4endl;
  }
#endif
 
   //----- set translation 
  physVol->SetTranslation( origin );
}

References FatalException, G4VDivisionParameterisation::faxis, G4VDivisionParameterisation::fmotherSolid, G4VDivisionParameterisation::foffset, G4VDivisionParameterisation::fwidth, G4cout, G4endl, G4Exception(), G4Trd::GetXHalfLength1(), G4Trd::GetXHalfLength2(), kXAxis, G4VPhysicalVolume::SetTranslation(), CLHEP::Hep3Vector::setX(), and G4VDivisionParameterisation::verbose.

◆ 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()

G4double G4ParameterisationTrdX::GetMaxParameter ( ) const

virtual

Implements G4VDivisionParameterisation.

Definition at line 123 of file G4ParameterisationTrd.cc.

{
  G4Trd* msol = (G4Trd*)(fmotherSolid);
  return (msol->GetXHalfLength1()+msol->GetXHalfLength2());
}

References G4VDivisionParameterisation::fmotherSolid, G4Trd::GetXHalfLength1(), and G4Trd::GetXHalfLength2().

◆ 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;
}  