#include <G4ParameterisationTrd.hh>

Inheritance diagram for G4ParameterisationTrdY:

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

	~G4ParameterisationTrdY ()

void	CheckParametersValidity ()

G4double	GetMaxParameter () const

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

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

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=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 G4VParameterisationTrd
G4bool	bDivInTrap

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

Constructor & Destructor Documentation

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

Definition at line 291 of file G4ParameterisationTrd.cc.

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

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

G4ParameterisationTrdY::~G4ParameterisationTrdY ( )

Definition at line 323 of file G4ParameterisationTrd.cc.

324 {

325 }

Member Function Documentation

void G4ParameterisationTrdY::CheckParametersValidity ( )

virtual

Reimplemented from G4VDivisionParameterisation.

Definition at line 397 of file G4ParameterisationTrd.cc.

References G4VDivisionParameterisation::CheckParametersValidity(), FatalException, G4VDivisionParameterisation::fmotherSolid, G4endl, G4Exception(), G4Trd::GetYHalfLength1(), G4Trd::GetYHalfLength2(), and G4VDivisionParameterisation::kCarTolerance.

Referenced by G4ParameterisationTrdY().

 {
   G4VDivisionParameterisation::CheckParametersValidity();
 
   G4Trd* msol = (G4Trd*)(fmotherSolid);
 
   G4double mpDy1 = msol->GetYHalfLength1();
   G4double mpDy2 = msol->GetYHalfLength2();
 
   if( std::fabs(mpDy1 - mpDy2) > kCarTolerance )
   {
     std::ostringstream message;
     message << "Invalid solid specification. NOT supported." << G4endl
             << "Making a division of a TRD along axis Y while" << G4endl
             << "the Y half lengths are not equal is not (yet)" << G4endl
             << "supported. It will result in non-equal" << G4endl
             << "division solids.";
     G4Exception("G4ParameterisationTrdY::CheckParametersValidity()",
                 "GeomDiv0001", FatalException, message);
   }
 }

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

virtual

Reimplemented from G4VPVParameterisation.

Definition at line 374 of file G4ParameterisationTrd.cc.

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

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

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

virtual

Implements G4VDivisionParameterisation.

Definition at line 337 of file G4ParameterisationTrd.cc.

References FatalException, G4VDivisionParameterisation::faxis, G4VDivisionParameterisation::fmotherSolid, G4VDivisionParameterisation::foffset, G4VDivisionParameterisation::fwidth, G4cout, G4endl, G4Exception(), G4Trd::GetYHalfLength1(), kYAxis, G4VPhysicalVolume::SetTranslation(), CLHEP::Hep3Vector::setY(), and G4VDivisionParameterisation::verbose.

 {
   G4Trd* msol = (G4Trd*)(fmotherSolid );
   G4double mdy = msol->GetYHalfLength1();
 
   //----- 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 );
 }

G4double G4ParameterisationTrdY::GetMaxParameter ( ) const

virtual

Implements G4VDivisionParameterisation.

Definition at line 328 of file G4ParameterisationTrd.cc.

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

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

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