G4ParameterisationTrdX Class Reference

#include <G4ParameterisationTrd.hh>

Inheritance diagram for G4ParameterisationTrdX:

Public Member Functions
	G4ParameterisationTrdX (EAxis axis, G4int nCopies, G4double width, G4double offset, G4VSolid *motherSolid, DivisionType divType)
	~G4ParameterisationTrdX ()
void	CheckParametersValidity ()
G4double	GetMaxParameter () const
void	ComputeTransformation (const G4int copyNo, G4VPhysicalVolume *physVol) const
void	ComputeDimensions (G4Trd &trd, const G4int copyNo, const G4VPhysicalVolume *pv) const
void	ComputeDimensions (G4Trap &trd, const G4int copyNo, const G4VPhysicalVolume *pv) const
G4VSolid *	ComputeSolid (const G4int, G4VPhysicalVolume *)
Public Member Functions inherited from G4VParameterisationTrd
	G4VParameterisationTrd (EAxis axis, G4int nCopies, G4double offset, G4double step, G4VSolid *msolid, DivisionType divType)
virtual	~G4VParameterisationTrd ()
Public Member Functions inherited from G4VDivisionParameterisation
	G4VDivisionParameterisation (EAxis axis, G4int nDiv, G4double width, G4double offset, DivisionType divType, G4VSolid *motherSolid=nullptr)
virtual	~G4VDivisionParameterisation ()
virtual G4VSolid *	ComputeSolid (const G4int, G4VPhysicalVolume *)
virtual void	ComputeTransformation (const G4int copyNo, G4VPhysicalVolume *physVol) const =0
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 void	ComputeTransformation (const G4int, G4VPhysicalVolume *) const =0
virtual G4VSolid *	ComputeSolid (const G4int, G4VPhysicalVolume *)
virtual G4Material *	ComputeMaterial (const G4int repNo, G4VPhysicalVolume *currentVol, const G4VTouchable *parentTouch=nullptr)
virtual G4bool	IsNested () const
virtual G4VVolumeMaterialScanner *	GetMaterialScanner ()
virtual void	ComputeDimensions (G4Box &, const G4int, const G4VPhysicalVolume *) const
virtual void	ComputeDimensions (G4Tubs &, const G4int, const G4VPhysicalVolume *) const
virtual void	ComputeDimensions (G4Trd &, const G4int, const G4VPhysicalVolume *) const
virtual void	ComputeDimensions (G4Trap &, const G4int, const G4VPhysicalVolume *) const
virtual void	ComputeDimensions (G4Cons &, const G4int, const G4VPhysicalVolume *) const
virtual void	ComputeDimensions (G4Sphere &, const G4int, const G4VPhysicalVolume *) const
virtual void	ComputeDimensions (G4Orb &, const G4int, const G4VPhysicalVolume *) const
virtual void	ComputeDimensions (G4Ellipsoid &, const G4int, const G4VPhysicalVolume *) const
virtual void	ComputeDimensions (G4Torus &, const G4int, const G4VPhysicalVolume *) const
virtual void	ComputeDimensions (G4Para &, const G4int, const G4VPhysicalVolume *) const
virtual void	ComputeDimensions (G4Polycone &, const G4int, const G4VPhysicalVolume *) const
virtual void	ComputeDimensions (G4Polyhedra &, const G4int, const G4VPhysicalVolume *) const
virtual void	ComputeDimensions (G4Hype &, const G4int, const G4VPhysicalVolume *) const

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)
virtual G4double	GetMaxParameter () const =0
G4double	OffsetZ () const
Protected Attributes inherited from G4VParameterisationTrd
G4bool	bDivInTrap = false
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 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.

{
}

Member Function Documentation

CheckParametersValidity()

void G4ParameterisationTrdX::CheckParametersValidity ( )

virtual

Reimplemented from G4VDivisionParameterisation.

Definition at line 255 of file G4ParameterisationTrd.cc.

{
  G4VDivisionParameterisation::CheckParametersValidity();
/*
  G4Trd* msol = (G4Trd*)(fmotherSolid);
 
  G4double mpDx1 = msol->GetXHalfLength1();
  G4double mpDx2 = msol->GetXHalfLength2();
  bDivInTrap = false;
 
  if( std::fabs(mpDx1 - mpDx2) > kCarTolerance )
  {
    std::ostringstream message;
    message << "Invalid solid specification. NOT supported." << G4endl
            << "Making a division of a TRD along axis X," << G4endl
            << "while the X half lengths are not equal," << G4endl
            << "is not (yet) supported. It will result" << G4endl
            << "in non-equal division solids.";
    G4Exception("G4ParameterisationTrdX::CheckParametersValidity()",
                "GeomDiv0001", FatalException, message);
  }
*/
}

Referenced by G4ParameterisationTrdX().

ComputeDimensions() [1/2]

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

virtual

Reimplemented from G4VPVParameterisation.

Definition at line 216 of file G4ParameterisationTrd.cc.

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

ComputeDimensions() [2/2]

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

virtual

Reimplemented from G4VPVParameterisation.

Definition at line 178 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
}

ComputeSolid()

G4VSolid * G4ParameterisationTrdX::ComputeSolid ( const G4int  i, G4VPhysicalVolume *  pv  )

virtual

Reimplemented from G4VDivisionParameterisation.

Definition at line 200 of file G4ParameterisationTrd.cc.

{
  if( bDivInTrap ) 
  {
    return G4VDivisionParameterisation::ComputeSolid(i, pv);
  } 
  else 
  {
    return fmotherSolid;
  }
}

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;
  if( !bDivInTrap )
  {
    posi = -mdx + foffset + (copyNo+0.5)*fwidth;
  }
  else
  {
    G4double aveHL = (msol->GetXHalfLength1()+msol->GetXHalfLength2())/2.;
    posi = - aveHL + foffset + (copyNo+0.5)*aveHL/fnDiv*2;
  }
  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 );
}

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

---

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

• G4ParameterisationTrd.hh
• G4ParameterisationTrd.cc