G4ParameterisationPolyhedraZ Class Reference

G4ParameterisationPolyhedraZ represents the parameterised positioning equivalent to dividing a G4Polyhedra along Z axis. More...

#include <G4ParameterisationPolyhedra.hh>

Inheritance diagram for G4ParameterisationPolyhedraZ:

Public Member Functions
	G4ParameterisationPolyhedraZ (EAxis axis, G4int nCopies, G4double offset, G4double step, G4VSolid *pSolid, DivisionType divType)
	~G4ParameterisationPolyhedraZ () override
void	CheckParametersValidity () override
G4double	GetMaxParameter () const override
void	ComputeTransformation (const G4int copyNo, G4VPhysicalVolume *physVol) const override
void	ComputeDimensions (G4Polyhedra &phedra, const G4int copyNo, const G4VPhysicalVolume *physVol) const override
Public Member Functions inherited from G4VParameterisationPolyhedra
	G4VParameterisationPolyhedra (EAxis axis, G4int nCopies, G4double offset, G4double step, G4VSolid *pSolid, DivisionType divType)
	~G4VParameterisationPolyhedra () override
Public Member Functions inherited from G4VDivisionParameterisation
	G4VDivisionParameterisation (EAxis axis, G4int nDiv, G4double width, G4double offset, DivisionType divType, G4VSolid *motherSolid=nullptr)
	~G4VDivisionParameterisation () override
G4VSolid *	ComputeSolid (const G4int, G4VPhysicalVolume *pv) override
const G4String &	GetType () const
EAxis	GetAxis () const
G4int	GetNoDiv () const
G4double	GetWidth () const
G4double	GetOffset () const
G4VSolid *	GetMotherSolid () const
void	SetType (const G4String &type)
void	SetHalfGap (G4double hg)
G4double	GetHalfGap () const
Public Member Functions inherited from G4VPVParameterisation
	G4VPVParameterisation ()=default
virtual	~G4VPVParameterisation ()=default
virtual G4Material *	ComputeMaterial (const G4int repNo, G4VPhysicalVolume *currentVol, const G4VTouchable *parentTouch=nullptr)
virtual G4bool	IsNested () const
virtual G4VVolumeMaterialScanner *	GetMaterialScanner ()

Additional Inherited Members
Protected Member Functions inherited from G4VDivisionParameterisation
void	ChangeRotMatrix (G4VPhysicalVolume *physVol, G4double rotZ=0.0) const
G4int	CalculateNDiv (G4double motherDim, G4double width, G4double offset) const
G4double	CalculateWidth (G4double motherDim, G4int nDiv, G4double offset) const
void	CheckOffset (G4double maxPar)
void	CheckNDivAndWidth (G4double maxPar)
G4double	OffsetZ () const
Protected Attributes inherited from G4VDivisionParameterisation
G4String	ftype
EAxis	faxis
G4int	fnDiv = 0
G4double	fwidth = 0.0
G4double	foffset = 0.0
DivisionType	fDivisionType
G4VSolid *	fmotherSolid = nullptr
G4bool	fReflectedSolid = false
G4bool	fDeleteSolid = false
G4double	kCarTolerance
G4double	fhgap = 0.0
Static Protected Attributes inherited from G4VDivisionParameterisation
static G4ThreadLocal G4RotationMatrix *	fRot = nullptr
static const G4int	verbose = 5

Detailed Description

G4ParameterisationPolyhedraZ represents the parameterised positioning equivalent to dividing a G4Polyhedra along Z axis.

Definition at line 253 of file G4ParameterisationPolyhedra.hh.

Constructor & Destructor Documentation

G4ParameterisationPolyhedraZ()

G4ParameterisationPolyhedraZ::G4ParameterisationPolyhedraZ	(	EAxis	axis,
		G4int	nCopies,
		G4double	offset,
		G4double	step,
		G4VSolid *	pSolid,
		DivisionType	divType )

Initialises a parameterised polyhedra, along the Z axis.

Parameters

[in]	axis	The axis along which apply the parameterisation.
[in]	nCopies	The total number of divided slices.
[in]	offset	Potential initial offset along the axis.
[in]	step	The width of the divided slice.
[in]	pSolid	Pointer to the original shape to parameterise.
[in]	divType	String identifier for the kind of division.

Definition at line 426 of file G4ParameterisationPolyhedra.cc.

  :  G4VParameterisationPolyhedra( axis, nDiv, width, offset, msolid, divType ),
     fOrigParamMother(((G4Polyhedra*)fmotherSolid)->GetOriginalParameters())
{ 
  CheckParametersValidity();
  SetType( "DivisionPolyhedraZ" );
 
  if( divType == DivWIDTH )
  {
    fnDiv =
      CalculateNDiv(fOrigParamMother->Z_values[fOrigParamMother->Num_z_planes-1]
                    - fOrigParamMother->Z_values[0] , width, offset);
  }
  else if( divType == DivNDIV )
  {
    fwidth =
      CalculateNDiv(fOrigParamMother->Z_values[fOrigParamMother->Num_z_planes-1]
                    - fOrigParamMother->Z_values[0] , nDiv, offset);
  }
  
#ifdef G4DIVDEBUG
  if( verbose >= 1 )
  {
    G4cout << " G4ParameterisationPolyhedraZ - # divisions " << fnDiv << " = "
           << nDiv << G4endl
           << " Offset " << foffset << " = " << offset << G4endl
           << " Width " << fwidth << " = " << width << G4endl;
  }
#endif
}

~G4ParameterisationPolyhedraZ()

G4ParameterisationPolyhedraZ::~G4ParameterisationPolyhedraZ ( )

overridedefault

Default Destructor.

Member Function Documentation

CheckParametersValidity()

void G4ParameterisationPolyhedraZ::CheckParametersValidity ( )

overridevirtual

Checks the validity of parameters given in input, issuing an exception.

Reimplemented from G4VDivisionParameterisation.

Definition at line 508 of file G4ParameterisationPolyhedra.cc.

{
  G4VDivisionParameterisation::CheckParametersValidity();
 
  // Division will be following the mother polyhedra segments
  //
  if( fDivisionType == DivNDIV )
  {
    if( fOrigParamMother->Num_z_planes-1 != fnDiv )
    { 
      std::ostringstream message;
      message << "Configuration not supported." << G4endl
              << "Division along Z will be done splitting in the defined"
              << G4endl
              << "Z planes, i.e, the number of division would be :"
              << fOrigParamMother->Num_z_planes-1 << " instead of "
              << fnDiv << " !"; 
      G4Exception("G4ParameterisationPolyhedraZ::CheckParametersValidity()",
                  "GeomDiv0001", FatalException, message);
    }
  }  
 
  // Division will be done within one polyhedra segment
  // with applying given width and offset
  //
  if( fDivisionType == DivNDIVandWIDTH || fDivisionType == DivWIDTH )
  {
    // Check if divided region does not span over more
    // than one z segment
 
    G4int isegstart = -1;  // number of the segment containing start position
    G4int isegend = -1;    // number of the segment containing end position
 
    if ( !fReflectedSolid )
    {
      // The start/end position of the divided region
      //
      G4double zstart = fOrigParamMother->Z_values[0] + foffset;
      G4double zend = fOrigParamMother->Z_values[0]
                    + foffset + fnDiv*fwidth;
   
      G4int counter = 0;
      while ( isegend < 0 && counter < fOrigParamMother->Num_z_planes-1 )
      {
        // first segment
        if ( zstart >= fOrigParamMother->Z_values[counter]  &&
             zstart  < fOrigParamMother->Z_values[counter+1] )
        {
           isegstart = counter;
        }     
        // last segment
        if ( zend  > fOrigParamMother->Z_values[counter] &&
             zend <= fOrigParamMother->Z_values[counter+1] )
        {
           isegend = counter;
        }   
        ++counter;   
      }  // Loop checking, 06.08.2015, G.Cosmo
    }
    else
    {
      // The start/end position of the divided region
      //
      G4double zstart = fOrigParamMother->Z_values[0] - foffset;
      G4double zend = fOrigParamMother->Z_values[0]
                    - (foffset + fnDiv* fwidth);
   
      G4int counter = 0;
      while ( isegend < 0 && counter < fOrigParamMother->Num_z_planes-1 )
      {
        // first segment
        if ( zstart <= fOrigParamMother->Z_values[counter]  &&
             zstart  > fOrigParamMother->Z_values[counter+1] )
        {
           isegstart = counter;
        }     
        // last segment
        if ( zend  < fOrigParamMother->Z_values[counter] &&
             zend >= fOrigParamMother->Z_values[counter+1] )
        {
           isegend = counter;
        }   
        ++counter;   
      }  // Loop checking, 06.08.2015, G.Cosmo
    }
  
    if ( isegstart != isegend )
    {
      std::ostringstream message;
      message << "Configuration not supported." << G4endl
              << "Division with user defined width." << G4endl
              << "Solid " << fmotherSolid->GetName() << G4endl
              << "Divided region is not between two Z planes."; 
      G4Exception("G4ParameterisationPolyhedraZ::CheckParametersValidity()",
                  "GeomDiv0001", FatalException, message);
    }
  
    fNSegment = isegstart;
  }  
}

Referenced by G4ParameterisationPolyhedraZ().

ComputeDimensions()

void G4ParameterisationPolyhedraZ::ComputeDimensions ( G4Polyhedra & phedra, const G4int copyNo, const G4VPhysicalVolume * physVol ) const

overridevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 668 of file G4ParameterisationPolyhedra.cc.

{
  // Define division solid
  //
  G4PolyhedraHistorical origparam;
  G4int nz = 2; 
  origparam.Num_z_planes = nz;
  origparam.numSide = fOrigParamMother->numSide;
  origparam.Start_angle = fOrigParamMother->Start_angle;
  origparam.Opening_angle = fOrigParamMother->Opening_angle;
 
  // Define division solid z sections
  //
  origparam.Z_values = new G4double[nz];
  origparam.Rmin = new G4double[nz];
  origparam.Rmax = new G4double[nz];
  origparam.Z_values[0] = - fwidth/2.;
  origparam.Z_values[1] = fwidth/2.;
 
  if ( fDivisionType == DivNDIV )
  {
    // The position of the centre of copyNo-th mother polycone segment
    //
    G4double posi = ( fOrigParamMother->Z_values[copyNo]
                    + fOrigParamMother->Z_values[copyNo+1])/2;
 
    origparam.Z_values[0] = fOrigParamMother->Z_values[copyNo] - posi;
    origparam.Z_values[1] = fOrigParamMother->Z_values[copyNo+1] - posi;
    origparam.Rmin[0] = fOrigParamMother->Rmin[copyNo];
    origparam.Rmin[1] = fOrigParamMother->Rmin[copyNo+1];
    origparam.Rmax[0] = fOrigParamMother->Rmax[copyNo];
    origparam.Rmax[1] = fOrigParamMother->Rmax[copyNo+1];
  }  
 
  if ( fDivisionType == DivNDIVandWIDTH || fDivisionType == DivWIDTH )
  {
    if ( !fReflectedSolid )
    {
      origparam.Z_values[0] = -fwidth/2.;
      origparam.Z_values[1] = fwidth/2.;
 
      // The position of the centre of copyNo-th division
      //
      G4double posi = fOrigParamMother->Z_values[0]
                    + foffset + (2*copyNo + 1) * fwidth/2.;
    
      // The first and last z sides z values
      G4double zstart = posi - fwidth/2.;
      G4double zend = posi + fwidth/2.;
      origparam.Rmin[0] = GetRmin(zstart, fNSegment); 
      origparam.Rmax[0] = GetRmax(zstart, fNSegment);  
      origparam.Rmin[1] = GetRmin(zend, fNSegment); 
      origparam.Rmax[1] = GetRmax(zend, fNSegment); 
    }
    else
    {   
      origparam.Z_values[0] = fwidth/2.;
      origparam.Z_values[1] = -fwidth/2.;
 
      // The position of the centre of copyNo-th division
      //
      G4double posi = fOrigParamMother->Z_values[0]
                    - ( foffset + (2*copyNo + 1) * fwidth/2.);
    
      // The first and last z sides z values
      //
      G4double zstart = posi + fwidth/2.;
      G4double zend = posi - fwidth/2.;
      origparam.Rmin[0] = GetRmin(zstart, fNSegment); 
      origparam.Rmax[0] = GetRmax(zstart, fNSegment);  
      origparam.Rmin[1] = GetRmin(zend, fNSegment); 
      origparam.Rmax[1] = GetRmax(zend, fNSegment); 
    } 
 
    // It can happen due to rounding errors
    //
    if ( origparam.Rmin[0]    < 0.0 ) { origparam.Rmin[0] = 0.0; }
    if ( origparam.Rmin[nz-1] < 0.0 ) { origparam.Rmin[1] = 0.0; }
  }  
 
  phedra.SetOriginalParameters(&origparam);  // copy values & transfer pointers
  phedra.Reset();                            // reset to new solid parameters
 
#ifdef G4DIVDEBUG
  if( verbose >= 2 )
  {
    G4cout << "G4ParameterisationPolyhedraZ::ComputeDimensions()" << G4endl
           << "-- Parametrised phedra copy-number: " << copyNo << G4endl;
    phedra.DumpInfo();
  }
#endif
}

ComputeTransformation()

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

overridevirtual

Concrete methods implementing the parameterisation.

Implements G4VDivisionParameterisation.

Definition at line 611 of file G4ParameterisationPolyhedra.cc.

{
  G4double posi;
  if ( fDivisionType == DivNDIV )
  {
    // The position of the centre of copyNo-th mother polycone segment
    
    posi = ( fOrigParamMother->Z_values[copyNo]
                    + fOrigParamMother->Z_values[copyNo+1])/2;
    physVol->SetTranslation( G4ThreeVector(0, 0, posi) );
  }
  
  if ( fDivisionType == DivNDIVandWIDTH || fDivisionType == DivWIDTH )
  {
    // The position of the centre of copyNo-th division
 
    posi = fOrigParamMother->Z_values[0];
    
    if ( !fReflectedSolid )
    {
      posi += foffset + (2*copyNo + 1) * fwidth/2.;
    }
    else
    {
      posi -= foffset + (2*copyNo + 1) * fwidth/2.;
    }
    
    physVol->SetTranslation( G4ThreeVector(0, 0, posi) );
  }   
 
  //----- calculate rotation matrix: unit
 
#ifdef G4DIVDEBUG
  if( verbose >= 2 )
  {
    G4cout << " G4ParameterisationPolyhedraZ - position: " << posi << G4endl
           << " copyNo: " << copyNo << " - foffset: " << foffset/CLHEP::deg
           << " - fwidth: " << fwidth/CLHEP::deg << G4endl;
  }
#endif
 
  ChangeRotMatrix( physVol );
 
#ifdef G4DIVDEBUG
  if( verbose >= 2 )
  {
    G4cout << std::setprecision(8) << " G4ParameterisationPolyhedraZ "
           << copyNo << G4endl
           << " Position: (0,0,0) - Width: " << fwidth
           << " - Axis: " << faxis  << G4endl;
  }
#endif
}

GetMaxParameter()

G4double G4ParameterisationPolyhedraZ::GetMaxParameter ( ) const

overridevirtual

Returns the max width along Z.

Returns

The maximum width of the solid to divide along the Z axis.

Implements G4VDivisionParameterisation.

Definition at line 501 of file G4ParameterisationPolyhedra.cc.

{
  return std::abs(fOrigParamMother->Z_values[fOrigParamMother->Num_z_planes-1]
                 -fOrigParamMother->Z_values[0]);
}

---

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

• G4ParameterisationPolyhedra.hh
• G4ParameterisationPolyhedra.cc