#include <G4ParameterisationPolyhedra.hh>

Inheritance diagram for G4ParameterisationPolyhedraZ:

Public Member Functions
	G4ParameterisationPolyhedraZ (EAxis axis, G4int nCopies, G4double offset, G4double step, G4VSolid *motherSolid, DivisionType divType)

	~G4ParameterisationPolyhedraZ ()

void	CheckParametersValidity ()

G4double	GetMaxParameter () const

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

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

Public Member Functions inherited from G4VParameterisationPolyhedra
	G4VParameterisationPolyhedra (EAxis axis, G4int nCopies, G4double offset, G4double step, G4VSolid *msolid, DivisionType divType)

virtual	~G4VParameterisationPolyhedra ()

Public Member Functions inherited from G4VDivisionParameterisation
	G4VDivisionParameterisation (EAxis axis, G4int nDiv, G4double width, G4double offset, DivisionType divType, G4VSolid *motherSolid=0)

virtual	~G4VDivisionParameterisation ()

virtual G4VSolid *	ComputeSolid (const G4int, G4VPhysicalVolume *)

const G4String &	GetType () const

EAxis	GetAxis () const

G4int	GetNoDiv () const

G4double	GetWidth () const

G4double	GetOffset () const

G4VSolid *	GetMotherSolid () const

void	SetType (const G4String &type)

G4int	VolumeFirstCopyNo () const

void	SetHalfGap (G4double hg)

G4double	GetHalfGap () const

Public Member Functions inherited from G4VPVParameterisation
	G4VPVParameterisation ()

virtual	~G4VPVParameterisation ()

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

virtual G4bool	IsNested () const

virtual G4VVolumeMaterialScanner *	GetMaterialScanner ()

Additional Inherited Members
Protected Member Functions inherited from G4VDivisionParameterisation
void	ChangeRotMatrix (G4VPhysicalVolume *physVol, G4double rotZ=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

G4double	fwidth

G4double	foffset

DivisionType	fDivisionType

G4VSolid *	fmotherSolid

G4bool	fReflectedSolid

G4bool	fDeleteSolid

G4int	theVoluFirstCopyNo

G4double	kCarTolerance

G4double	fhgap

Static Protected Attributes inherited from G4VDivisionParameterisation
static G4ThreadLocal G4int	verbose = 5

Detailed Description

Definition at line 188 of file G4ParameterisationPolyhedra.hh.

Constructor & Destructor Documentation

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

Definition at line 426 of file G4ParameterisationPolyhedra.cc.

References G4VDivisionParameterisation::CalculateNDiv(), CheckParametersValidity(), DivNDIV, DivWIDTH, G4VDivisionParameterisation::fnDiv, G4VDivisionParameterisation::foffset, G4VDivisionParameterisation::fwidth, G4cout, G4endl, G4PolyhedraHistorical::Num_z_planes, G4VDivisionParameterisation::SetType(), G4VDivisionParameterisation::verbose, and G4PolyhedraHistorical::Z_values.

   :  G4VParameterisationPolyhedra( axis, nDiv, width, offset, msolid, divType ),
      fNSegment(0),
      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 ( )

Definition at line 461 of file G4ParameterisationPolyhedra.cc.

462 {

463 }

Member Function Documentation

void G4ParameterisationPolyhedraZ::CheckParametersValidity ( )

virtual

Reimplemented from G4VDivisionParameterisation.

Definition at line 510 of file G4ParameterisationPolyhedra.cc.

References G4VDivisionParameterisation::CheckParametersValidity(), DivNDIV, DivNDIVandWIDTH, DivWIDTH, FatalException, G4VDivisionParameterisation::fDivisionType, G4VDivisionParameterisation::fmotherSolid, G4VDivisionParameterisation::fnDiv, G4VDivisionParameterisation::foffset, G4VDivisionParameterisation::fReflectedSolid, G4VDivisionParameterisation::fwidth, G4endl, G4Exception(), G4VSolid::GetName(), G4PolyhedraHistorical::Num_z_planes, and G4PolyhedraHistorical::Z_values.

Referenced by G4ParameterisationPolyhedraZ().

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

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

virtual

Reimplemented from G4VPVParameterisation.

Definition at line 649 of file G4ParameterisationPolyhedra.cc.

References DivNDIV, DivNDIVandWIDTH, DivWIDTH, G4VSolid::DumpInfo(), G4VDivisionParameterisation::fDivisionType, G4VDivisionParameterisation::foffset, G4VDivisionParameterisation::fReflectedSolid, G4VDivisionParameterisation::fwidth, G4cout, G4endl, G4PolyhedraHistorical::Num_z_planes, G4PolyhedraHistorical::numSide, G4PolyhedraHistorical::Opening_angle, G4Polyhedra::Reset(), G4PolyhedraHistorical::Rmax, G4PolyhedraHistorical::Rmin, G4Polyhedra::SetOriginalParameters(), G4PolyhedraHistorical::Start_angle, G4VDivisionParameterisation::verbose, and G4PolyhedraHistorical::Z_values.

 {
   // 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
 }

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

virtual

Implements G4VDivisionParameterisation.

Definition at line 600 of file G4ParameterisationPolyhedra.cc.

References G4VDivisionParameterisation::ChangeRotMatrix(), python.hepunit::deg, DivNDIV, DivNDIVandWIDTH, DivWIDTH, G4VDivisionParameterisation::faxis, G4VDivisionParameterisation::fDivisionType, G4VDivisionParameterisation::foffset, G4VDivisionParameterisation::fReflectedSolid, G4VDivisionParameterisation::fwidth, G4cout, G4endl, G4VPhysicalVolume::SetTranslation(), G4VDivisionParameterisation::verbose, and G4PolyhedraHistorical::Z_values.

 {
   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;
     physVol->SetTranslation( G4ThreeVector(0, 0, posi) );
   }
   
   if ( fDivisionType == DivNDIVandWIDTH || fDivisionType == DivWIDTH ) {
     // The position of the centre of copyNo-th division
 
     G4double 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/deg
            << " - fwidth: " << fwidth/deg << G4endl;
   }
 #endif
 
   ChangeRotMatrix( physVol );
 
 #ifdef G4DIVDEBUG
   if( verbose >= 2 )
   {
     G4cout << std::setprecision(8) << " G4ParameterisationPolyhedraZ "
            << copyNo << G4endl
            << " Position: " << origin << " - Width: " << fwidth
            << " - Axis: " << faxis  << G4endl;
   }
 #endif
 }

G4double G4ParameterisationPolyhedraZ::GetMaxParameter ( ) const

virtual

Implements G4VDivisionParameterisation.

Definition at line 503 of file G4ParameterisationPolyhedra.cc.

References G4PolyhedraHistorical::Num_z_planes, and G4PolyhedraHistorical::Z_values.

 {
   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:

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation