#include <G4ParameterisationTrd.hh>

Inheritance diagram for G4ParameterisationTrdY:

Public Member Functions
void	ComputeDimensions (G4Trap &, const G4int, const G4VPhysicalVolume *) 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

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

	~G4ParameterisationTrdY ()

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 122 of file G4ParameterisationTrd.hh.

Constructor & Destructor Documentation

◆ G4ParameterisationTrdY()

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

Definition at line 242 of file G4ParameterisationTrd.cc.

  : G4VParameterisationTrd( axis, nDiv, width, offset, msolid, divType )
{
  CheckParametersValidity();
  SetType( "DivisionTrdY" );
 
  G4Trd* msol = (G4Trd*)(fmotherSolid);
  if( divType == DivWIDTH )
  {
    fnDiv = CalculateNDiv( 2*msol->GetYHalfLength1(),
                           width, offset );
  }
  else if( divType == DivNDIV )
  {
    fwidth = CalculateWidth( 2*msol->GetYHalfLength1(),
                             nDiv, offset );
  }
 
#ifdef G4DIVDEBUG
  if( verbose >= 1 )
  {
     G4cout << " G4ParameterisationTrdY no divisions " << fnDiv
            << " = " << nDiv << G4endl
            << " Offset " << foffset << " = " << offset << G4endl
            << " width " << fwidth << " = " << width << G4endl;
  }
#endif
}

References G4VDivisionParameterisation::CalculateNDiv(), G4VDivisionParameterisation::CalculateWidth(), G4VDivisionParameterisation::CheckParametersValidity(), DivNDIV, DivWIDTH, G4VDivisionParameterisation::fmotherSolid, G4VDivisionParameterisation::fnDiv, G4VDivisionParameterisation::foffset, G4VDivisionParameterisation::fwidth, G4cout, G4endl, G4Trd::GetYHalfLength1(), G4VDivisionParameterisation::SetType(), and G4VDivisionParameterisation::verbose.

◆ ~G4ParameterisationTrdY()

G4ParameterisationTrdY::~G4ParameterisationTrdY ( )

Definition at line 275 of file G4ParameterisationTrd.cc.

276{

277}

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 146 of file G4ParameterisationTrd.hh.

147 {}

◆ ComputeDimensions() [2/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 144 of file G4ParameterisationTrd.hh.

145 {}

◆ ComputeDimensions() [3/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 152 of file G4ParameterisationTrd.hh.

153 {}

◆ ComputeDimensions() [4/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 158 of file G4ParameterisationTrd.hh.

159 {}

◆ ComputeDimensions() [5/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 150 of file G4ParameterisationTrd.hh.

151 {}

◆ ComputeDimensions() [6/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 156 of file G4ParameterisationTrd.hh.

157 {}

◆ ComputeDimensions() [7/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 162 of file G4ParameterisationTrd.hh.

163 {}

◆ ComputeDimensions() [8/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 164 of file G4ParameterisationTrd.hh.

165 {}

◆ ComputeDimensions() [9/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 148 of file G4ParameterisationTrd.hh.

149 {}

◆ ComputeDimensions() [10/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 154 of file G4ParameterisationTrd.hh.

155 {}

◆ ComputeDimensions() [11/13]

void G4ParameterisationTrdY::ComputeDimensions	(	G4Trap &	trap,
		const G4int	copyNo,
		const G4VPhysicalVolume *
	)		const

virtual

Reimplemented from G4VPVParameterisation.

Definition at line 350 of file G4ParameterisationTrd.cc.

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

References CLHEP::deg, 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 G4ParameterisationTrdY::ComputeDimensions	(	G4Trd &	trd,
		const G4int	copyNo,
		const G4VPhysicalVolume *	pv
	)		const

virtual

Reimplemented from G4VPVParameterisation.

Definition at line 325 of file G4ParameterisationTrd.cc.

{
  //---- The division along Y of a Trd will result a Trd, only 
  //--- if Y at -Z and +Z are equal, else use the G4Trap version
  G4Trd* msol = (G4Trd*)(fmotherSolid);
  
  G4double pDx1 = msol->GetXHalfLength1();
  G4double pDx2 = msol->GetXHalfLength2();
  G4double pDz = msol->GetZHalfLength();
  G4double pDy = fwidth/2. - fhgap;
 
  trd.SetAllParameters ( pDx1, pDx2, pDy, pDy, pDz );
 
#ifdef G4DIVDEBUG
  if( verbose >= 2 )
  {
    G4cout << " G4ParameterisationTrdY::ComputeDimensions():" << G4endl;
    trd.DumpInfo();
  }
#endif
}

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

◆ ComputeDimensions() [13/13]

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

inlineprivatevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 160 of file G4ParameterisationTrd.hh.

161 {}

◆ 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 G4ParameterisationTrdY::ComputeTransformation	(	const G4int	copyNo,
		G4VPhysicalVolume *	physVol
	)		const

virtual

Implements G4VDivisionParameterisation.

Definition at line 288 of file G4ParameterisationTrd.cc.

{
  G4Trd* msol = (G4Trd*)(fmotherSolid );
  G4double mdy =  ( msol->GetYHalfLength1() + msol->GetYHalfLength2() ) / 2.;
 
  //----- translation 
  G4ThreeVector origin(0.,0.,0.); 
  G4double posi = -mdy + foffset + (copyNo+0.5)*fwidth;
 
  if( faxis == kYAxis )
  {
    origin.setY( posi ); 
  }
  else
  { 
    std::ostringstream message;
    message << "Only axes along Y are allowed !  Axis: " << faxis;
    G4Exception("G4ParameterisationTrdY::ComputeTransformation()",
                "GeomDiv0002", FatalException, message);
  }
 
#ifdef G4DIVDEBUG
  if( verbose >= 2 )
  {
    G4cout << std::setprecision(8)
           << " G4ParameterisationTrdY::ComputeTransformation " << copyNo
           << " pos " << origin << " rot mat " << " axis " << faxis << G4endl;
  }
#endif
 
   //----- set translation 
  physVol->SetTranslation( origin );
}

References FatalException, G4VDivisionParameterisation::faxis, G4VDivisionParameterisation::fmotherSolid, G4VDivisionParameterisation::foffset, G4VDivisionParameterisation::fwidth, G4cout, G4endl, G4Exception(), G4Trd::GetYHalfLength1(), G4Trd::GetYHalfLength2(), kYAxis, G4VPhysicalVolume::SetTranslation(), CLHEP::Hep3Vector::setY(), 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 G4ParameterisationTrdY::GetMaxParameter ( ) const

virtual

Implements G4VDivisionParameterisation.

Definition at line 280 of file G4ParameterisationTrd.cc.

{
  G4Trd* msol = (G4Trd*)(fmotherSolid);
  return (msol->GetYHalfLength1()+msol->GetYHalfLength2());
}

References G4VDivisionParameterisation::fmotherSolid, G4Trd::GetYHalfLength1(), and G4Trd::GetYHalfLength2().

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