G4tgbPlaceParamCircle Class Reference

#include <G4tgbPlaceParamCircle.hh>

Inheritance diagram for G4tgbPlaceParamCircle:

Public Member Functions
	G4tgbPlaceParamCircle (G4tgrPlaceParameterisation *)
	~G4tgbPlaceParamCircle ()
void	ComputeTransformation (const G4int copyNo, G4VPhysicalVolume *physVol) const

---

Detailed Description

Definition at line 52 of file G4tgbPlaceParamCircle.hh.

---

Constructor & Destructor Documentation

G4tgbPlaceParamCircle::G4tgbPlaceParamCircle

(

G4tgrPlaceParameterisation *

)

Definition at line 52 of file G4tgbPlaceParamCircle.cc.

References G4tgbPlaceParameterisation::CheckNExtraData(), FatalException, G4cout, G4endl, G4Exception(), G4tgrPlaceParameterisation::GetExtraData(), G4tgrPlaceParameterisation::GetParamType(), G4tgrMessenger::GetVerboseLevel(), kXAxis, kYAxis, kZAxis, G4tgbPlaceParameterisation::theAxis, G4tgbPlaceParameterisation::theNCopies, and WLSIZE_EQ.

  : G4tgbPlaceParameterisation(tgrParam)
{
  //---- Get translation and rotation 
  if( tgrParam->GetParamType() == "CIRCLE" )
  {
    CheckNExtraData( tgrParam, 7, WLSIZE_EQ, "G4tgbPlaceParamCircle:");
    theCircleAxis = G4ThreeVector( tgrParam->GetExtraData()[4],
                                   tgrParam->GetExtraData()[5],
                                   tgrParam->GetExtraData()[6] );

    G4ThreeVector zaxis(0.,0.,-1.);
    if( zaxis.cross(theCircleAxis).mag() > 1.E-6 )
    {
      theDirInPlane = zaxis.cross(theCircleAxis);
    }
    else
    { 
      theDirInPlane = theCircleAxis.cross(G4ThreeVector(0.,-1.,0.));
    }
    theAxis = kZAxis;
  }
  else
  {
    CheckNExtraData( tgrParam, 4, WLSIZE_EQ, "G4tgbPlaceParamCircle:");
    if( tgrParam->GetParamType() == "CIRCLE_XY" ) {
      theCircleAxis = G4ThreeVector(0.,0.,1.);
      theDirInPlane = G4ThreeVector(1.,0.,0.);
      theAxis = kZAxis;
    } else if( tgrParam->GetParamType() == "CIRCLE_XZ" ) {
      theCircleAxis = G4ThreeVector(0.,1.,0.);
      theDirInPlane = G4ThreeVector(1.,0.,0.);
      theAxis = kYAxis;
    } else if( tgrParam->GetParamType() == "CIRCLE_YZ" ) {
      theCircleAxis = G4ThreeVector(1.,0.,0.);
      theDirInPlane = G4ThreeVector(0.,1.,0.);
      theAxis = kXAxis;
    }
  }

  if( theCircleAxis.mag() == 0. )
  {
    G4Exception("G4tgbPlaceParamCircle::G4tgbPlaceParamCircle()",
                "InvalidSetup", FatalException, "Circle axis is zero !");
  }
  theCircleAxis /= theCircleAxis.mag();

  theAxis = kZAxis;

  theNCopies = G4int(tgrParam->GetExtraData()[0]);
  theStep = tgrParam->GetExtraData()[1];
  theOffset = tgrParam->GetExtraData()[2];
  theRadius = tgrParam->GetExtraData()[3];

#ifdef G4VERBOSE
  if( G4tgrMessenger::GetVerboseLevel() >= 2 )
  {
    G4cout << " G4tgbPlaceParamCircle::G4tgbPlaceParamCircle():" << G4endl
           << " param type " << tgrParam->GetParamType() << G4endl
           << "   no copies - " << theNCopies << G4endl
           << "   step - " << theStep << G4endl
           << "   offset - " << theOffset << G4endl
           << "   radius - " << theRadius << G4endl
           << "   circle axis - " << theCircleAxis << G4endl
           << "   dir in plane - " << theDirInPlane << G4endl;
  }
#endif
}

G4tgbPlaceParamCircle::~G4tgbPlaceParamCircle

(

)

Definition at line 45 of file G4tgbPlaceParamCircle.cc.

{
}

---

Member Function Documentation

void G4tgbPlaceParamCircle::ComputeTransformation	(	const G4int	copyNo,
		G4VPhysicalVolume *	physVol
	)			const [virtual]

Reimplemented from G4tgbPlaceParameterisation.

Definition at line 124 of file G4tgbPlaceParamCircle.cc.

References G4cout, G4endl, G4VPhysicalVolume::GetName(), G4VPhysicalVolume::GetRotation(), G4tgrMessenger::GetVerboseLevel(), G4VPhysicalVolume::SetCopyNo(), G4VPhysicalVolume::SetRotation(), G4VPhysicalVolume::SetTranslation(), G4tgbPlaceParameterisation::theNCopies, and G4tgbPlaceParameterisation::theRotationMatrix.

{ 
  G4double posi = theOffset + copyNo*theStep;
  G4ThreeVector origin = theDirInPlane * theRadius;
  origin.rotate( posi, theCircleAxis );

  //----- Calculate rotation matrix (so that all volumes point to the centre)
  G4RotationMatrix rm;
  rm.rotate( -posi, theCircleAxis );

  //----- Set translation and rotation
  physVol->SetTranslation(origin);
  G4RotationMatrix* pvRm = physVol->GetRotation();
  if( pvRm == 0 )
  {
    pvRm = new G4RotationMatrix;
  }
  *pvRm  = *theRotationMatrix * rm;
  physVol->SetRotation(pvRm);
  physVol->SetCopyNo( copyNo );

#ifdef G4VERBOSE
  if( G4tgrMessenger::GetVerboseLevel() >= 3 )
  {
    G4cout << " G4tgbPlaceParamCircle::ComputeTransformation():" 
           << physVol->GetName() << G4endl
           << "   no copies - " << theNCopies  << G4endl
           << "   centre - " << origin << G4endl
           << "   rotation-matrix - " << *pvRm << G4endl;
  }
#endif
}

---

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

• G4tgbPlaceParamCircle.hh
• G4tgbPlaceParamCircle.cc

---