G4ParameterisationConsRho Class Reference

#include <G4ParameterisationCons.hh>

Inheritance diagram for G4ParameterisationConsRho:

Public Member Functions
	G4ParameterisationConsRho (EAxis axis, G4int nCopies, G4double offset, G4double step, G4VSolid *motherSolid, DivisionType divType)
	~G4ParameterisationConsRho () override
G4double	GetMaxParameter () const override
void	ComputeTransformation (const G4int copyNo, G4VPhysicalVolume *physVol) const override
void	ComputeDimensions (G4Cons &tubs, const G4int copyNo, const G4VPhysicalVolume *physVol) const override
Public Member Functions inherited from G4VParameterisationCons
	G4VParameterisationCons (EAxis axis, G4int nCopies, G4double offset, G4double step, G4VSolid *msolid, DivisionType divType)
	~G4VParameterisationCons () 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 *) 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)
G4int	VolumeFirstCopyNo () const
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
virtual void	CheckParametersValidity ()
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
G4int	theVoluFirstCopyNo = 1
G4double	kCarTolerance
G4double	fhgap = 0.0
Static Protected Attributes inherited from G4VDivisionParameterisation
static G4ThreadLocal G4RotationMatrix *	fRot = nullptr
static const G4int	verbose = 5

Detailed Description

Definition at line 70 of file G4ParameterisationCons.hh.

Constructor & Destructor Documentation

G4ParameterisationConsRho()

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

Definition at line 76 of file G4ParameterisationCons.cc.

  :  G4VParameterisationCons( axis, nDiv, width, offset, msolid, divType )
{
  CheckParametersValidity();
  SetType( "DivisionConsRho" );
 
  auto msol = (G4Cons*)(fmotherSolid);
  if( msol->GetInnerRadiusPlusZ() == 0. )
  {
    std::ostringstream message;
    message << "OuterRadiusMinusZ = 0" << G4endl 
            << "Width is calculated as that of OuterRadiusMinusZ !";
    G4Exception("G4ParameterisationConsRho::G4ParameterisationConsRho()",
                "GeomDiv1001", JustWarning, message);
  } 
 
  if( divType == DivWIDTH )
  {
    fnDiv = CalculateNDiv( msol->GetOuterRadiusMinusZ()
                         - msol->GetInnerRadiusMinusZ(), width, offset );
  }
  else if( divType == DivNDIV )
  {
    auto mconsol = (G4Cons*)(msolid);
    fwidth = CalculateWidth( mconsol->GetOuterRadiusMinusZ()
                           - mconsol->GetInnerRadiusMinusZ(), nDiv, offset );
  }
 
#ifdef G4DIVDEBUG
  if( verbose >= 1 )
  {
    G4cout << " G4ParameterisationConsRho - no divisions " << fnDiv << " = "
           << nDiv << G4endl
           << " Offset " << foffset << " = " << offset
           << " - Width " << fwidth << " = " << width << G4endl;
  }
#endif
}

~G4ParameterisationConsRho()

G4ParameterisationConsRho::~G4ParameterisationConsRho ( )

overridedefault

Member Function Documentation

ComputeDimensions()

void G4ParameterisationConsRho::ComputeDimensions ( G4Cons & tubs, const G4int copyNo, const G4VPhysicalVolume * physVol ) const

overridevirtual

Reimplemented from G4VPVParameterisation.

Definition at line 163 of file G4ParameterisationCons.cc.

{
  auto msol = (G4Cons*)(fmotherSolid);
 
  G4double pRMin1 = msol->GetInnerRadiusMinusZ() + foffset + fwidth*copyNo;
  G4double pRMax1 = msol->GetInnerRadiusMinusZ() + foffset + fwidth*(copyNo+1);
 
  //width at Z Plus
  //- G4double fwidthPlus =
  //   fwidth * ( msol->GetOuterRadiusPlusZ()/ msol->GetInnerRadiusPlusZ())
  //-         / ( msol->GetOuterRadiusMinusZ() - msol->GetInnerRadiusMinusZ());
  G4double fwidthPlus = CalculateWidth( msol->GetOuterRadiusPlusZ()
                           - msol->GetInnerRadiusPlusZ(), fnDiv, foffset );
  G4double pRMin2 = msol->GetInnerRadiusPlusZ()
                  + foffset + fwidthPlus * copyNo;
  G4double pRMax2 = msol->GetInnerRadiusPlusZ()
                  + foffset + fwidthPlus * (copyNo+1);
  G4double pDz = msol->GetZHalfLength();
 
  G4double d_half_gap = fhgap * pRMax2 / pRMax1;
  //- already rotated  double pSR = foffset + copyNo*fwidth;
  G4double pSPhi = msol->GetStartPhiAngle();
  G4double pDPhi = msol->GetDeltaPhiAngle();;
 
  cons.SetInnerRadiusMinusZ( pRMin1 + fhgap );
  cons.SetOuterRadiusMinusZ( pRMax1 - fhgap );
  cons.SetInnerRadiusPlusZ( pRMin2 + d_half_gap );
  cons.SetOuterRadiusPlusZ( pRMax2 - d_half_gap );
  cons.SetZHalfLength( pDz );
  cons.SetStartPhiAngle( pSPhi, false );
  cons.SetDeltaPhiAngle( pDPhi );
 
#ifdef G4DIVDEBUG
  if( verbose >= 2 )
  {
    G4cout << " G4ParameterisationConsRho::ComputeDimensions()" << G4endl
           << " pRMin: " << pRMin1 << " - pRMax: " << pRMax1 << G4endl;
    if( verbose >= 4 ) cons.DumpInfo();
  }
#endif
}

ComputeTransformation()

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

overridevirtual

Implements G4VDivisionParameterisation.

Definition at line 130 of file G4ParameterisationCons.cc.

{
  //----- translation 
  G4ThreeVector origin(0.,0.,0.); 
  //----- set translation 
  physVol->SetTranslation( origin );
 
  //----- calculate rotation matrix: unit
 
#ifdef G4DIVDEBUG
  if( verbose >= 2 )
  {
    G4cout <<  " G4ParameterisationConsRho " << G4endl
           << " Offset: " << foffset
           << " - Width: " << fwidth << G4endl;
  }
#endif
 
  ChangeRotMatrix( physVol );
 
#ifdef G4DIVDEBUG
  if( verbose >= 2 )
  {
    G4cout << std::setprecision(8) << " G4ParameterisationConsRho" << G4endl
           << " Position: " << origin << " - Width: " << fwidth
           << " - Axis: " << faxis  << G4endl;
  }
#endif
}

GetMaxParameter()

G4double G4ParameterisationConsRho::GetMaxParameter ( ) const

overridevirtual

Implements G4VDivisionParameterisation.

Definition at line 122 of file G4ParameterisationCons.cc.

{
  auto msol = (G4Cons*)(fmotherSolid);
  return msol->GetOuterRadiusMinusZ() - msol->GetInnerRadiusMinusZ();
}

---

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

• G4ParameterisationCons.hh
• G4ParameterisationCons.cc